Properties of carbon-oxygen white dwarfs from Monte Carlo stellar models

We investigate properties of carbon-oxygen white dwarfs with respect to the composite uncertainties in the reaction rates using the stellar evolution toolkit, Modules for Experiments in Stellar Astrophysics (MESA) and the probability density functions in the reaction rate library STARLIB. These are the first Monte Carlo stellar evolution studies that use complete stellar models. Focusing on 3 models evolved from the pre main-sequence to the first thermal pulse, we survey the remnant core mass, composition, and structure properties as a function of 26 STARLIB reaction rates covering hydrogen and helium burning using a Principal Component Analysis and Spearman Rank-Order Correlation. Relative to the arithmetic mean value, we find the width of the 95% confidence interval to be ΔM_1TP≈0.019 M_⊙ for the core mass at the first thermal pulse, Δ_t1TP ≈ 12.50 Myr for the age, Δlog(T_c/K)≈0.013 for the central temperature, Δlog(ρ_c/g cm^-3)≈0.060 for the central density, ΔY_e,c≈2.6 × 10^-5 for the central electron fraction, ΔX_c(²²Ne)≈5.8 × 10^-4, ΔX_c(12^C)≈0.392, and ΔX_c(¹⁶O)≈0.392. Uncertainties in the experimental ¹²C(α,γ)¹⁶O triple-α, and ¹⁴N(p,γ)¹⁵O reaction rates dominate these variations. We also consider a grid of 1-6 M_⊙ models evolved from the pre main-sequence to the final white dwarf to probe the sensitivity of the initial-final mass relation to experimental uncertainties in the hydrogen and helium reaction rates.