2017/12/17 by C. E. Fields, F. X. Timmes, R. Farmer +3
Physics and Astronomy · #Astro and Planetary Science #Gamma-ray bursts and supernovae #Isotope #Monte Carlo method #Nuclear physics research studies #Nuclear reaction #Nucleosynthesis #Reaction rate #Stellar evolution #Supernova #Thermonuclear fusion #astro-ph.HE #astro-ph.SR #nucl-ex #nucl-th
paper · pdf · doi:10.3847/1538-4365/aaa29b
Accepted for publication in ApJS, 33 pages, 21 figures
arxiv created 2017/12/17 · openalex created_date 2018/01/05 · openalex publication_date 2018/01/25 · arxiv updated 2018/02/07 · openalex updated_date 2026/08/05
Abstract We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the probability density functions of the reaction rates provided by the STARLIB reaction rate library with MESA stellar models. We evolve 1000 models of 15 from the pre-main sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a MESA in situ reaction network that follows 127 isotopes from 1 H to 64 Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman rank-order correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in the reaction rates of , triple- α , , 12 C( 12 C, p) 23 Na, 12 C( 16 O, p) 27 Al, 16 O( 16 O, n) 31 S, 16 O( 16 O, p) 31 P, and 16 O( 16 O, α ) 28 Si dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion they are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation, causing uncertainty in various properties of the stellar model in the evolution toward iron core-collapse.