About 30% to 40% of classical novae produce dust 20-100 days after the outburst, but no presolar stardust grains from classical novae have been unambiguously identified yet. Although several studies claimed a nova paternity for certain grains, the measured and simulated isotopic ratios could only be reconciled assuming that the grains condensed after the nova ejecta mixed with a much larger amount of close-to-solar matter. However, the source and mechanism of this potential post-explosion dilution of the ejecta remains a mystery. A major problem with previous studies is the small number of simulations performed and the implied poor exploration of the large nova parameter space. We report the results of a different strategy, based on a Monte Carlo technique, that involves the random sampling over the most important nova model parameters: the white dwarf composition; the mixing of the outer white dwarf layers with the accreted material before the explosion; the peak temperature and density; the explosion time scales; and the possible dilution of the ejecta after the outburst. We discuss and take into account the systematic uncertainties for both the presolar grain measurements and the simulation results. Only those simulations that are consistent with all measured isotopic ratios of a given grain are accepted for further analysis. We also present the numerical results of the model parameters. We identify 18 presolar grains with measured isotopic signatures consistent with a CO nova origin, without assuming any dilution of the ejecta. Among these, the grains G270 2, M11-334-2, G278, M11-347-4, M11-151-4, and Ag2 6 have the highest probability of a CO nova paternity.

Classical novae are caused by mass transfer episodes from a main sequence star onto a white dwarf via Roche lobe overflow. This material forms an accretion disk around the white dwarf. Ultimately, a fraction of this material spirals in and piles up on the white dwarf surface under electron-degenerate conditions. The subsequently occurring thermonuclear runaway reaches hundreds of megakelvin and explosively ejects matter into the interstellar medium. The exact peak temperature strongly depends on the underlying white dwarf mass, the accreted mass and metallicity, and the initial white dwarf luminosity. Observations of elemental abundance enrichments in these classical nova events imply that the ejected matter consists not only of processed solar material from the main sequence partner but also of material from the outer layers of the underlying white dwarf. This indicates that white dwarf and accreted matter mix prior to the thermonuclear runaway. The processes by which this mixing occurs require further investigation to be understood. In this work, we analyze elemental abundances ejected from hydrodynamic nova models in search of elemental abundance ratios that are useful indicators of the total amount of mixing. We identify the abundance ratios $\Sigma$CNO/H, Ne/H, Mg/H, Al/H, and Si/H as useful mixing meters in ONe novae. The impact of thermonuclear reaction rate uncertainties on the mixing meters is investigated using Monte Carlo post-processing network calculations with temperature-density evolutions of all mass zones computed by the hydrodynamic models. We find that the current uncertainties in the $^{30}$P($p$,$\gamma$)$^{31}$S rate influence the Si/H abundance ratio, but overall the mixing meters found here are robust against nuclear physics uncertainties. A comparison of our results with observations of ONe novae provides strong constraints for classical nova models.

Classical novae are stellar explosions occurring in binary systems, consisting of a white dwarf and a main sequence companion. Thermonuclear runaways on the surface of massive white dwarfs, consisting of oxygen and neon, are believed to reach peak temperatures of several hundred million kelvin. These temperatures are strongly correlated with the underlying white dwarf mass. The observational counterparts of such models are likely associated with outbursts that show strong spectral lines of neon in their shells (neon novae). The goals of this work are to investigate how useful elemental abundances are for constraining the peak temperatures achieved during these outbursts and determine how robust "nova thermometers" are with respect to uncertain nuclear physics input. We present updated observed abundances in neon novae and perform a series of hydrodynamic simulations for several white dwarf masses. We find that the most useful thermometers, N/O, N/Al, O/S, S/Al, O/Na, Na/Al, O/P, and P/Al, are those with the steepest monotonic dependence on peak temperature. The sensitivity of these thermometers to thermonuclear reaction rate variations is explored using post-processing nucleosynthesis simulations. The ratios N/O, N/Al, O/Na, and Na/Al are robust, meaning they are minimally affected by uncertain rates. However, their dependence on peak temperature is relatively weak. The ratios O/S, S/Al, O/P, and P/Al reveal strong dependences on temperature and the poorly known 30P(p,g)31S rate. We compare our model predictions to neon nova observations and obtain the following estimates for the underlying white dwarf masses: 1.34-1.35 solar masses (V838 Her), 1.18-1.21 solar masses (V382 Vel), <1.3 solar masses (V693 CrA), <1.2 solar masses (LMC 1990#1), and <1.2 solar masses (QU Vul).