Showing 1–50 of 126 results for author: Timmes, F

An Expanded Set of Los Alamos OPLIB Tables in MESA: Type-1 Rosseland-mean Opacities and Solar Models

Authors: Ebraheem Farag, Christopher J. Fontes, F. X. Timmes, Earl P. Bellinger, Joyce A. Guzik, Evan B. Bauer, Suzannah R. Wood, Katie Mussack, Peter Hakel, James Colgan, David P. Kilcrease, Manolo E. Sherrill, Tryston C. Raecke, Morgan T. Chidester

Abstract: We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10… ▽ More We present a set of 1194 Type-1 Rosseland-mean opacity tables for four different metallicity mixtures. These new Los Alamos OPLIB atomic radiative opacity tables are an order of magnitude larger in number than any previous opacity table release, and span regimes where previous opacity tables have not existed. For example, the new set of opacity tables expands the metallicity range to $Z$\,=\,10$^{-6}$ to $Z$\,=\,0.2 which allows improved accuracy of opacities at low and high metallicity, increases the table density in the metallicity range $Z$\,=\,10$^{-4}$ to $Z$\,=\,0.1 to enhance the accuracy of opacities drawn from interpolations across neighboring metallicities, and adds entries for hydrogen mass fractions between $X$\,=\,0 and $X$\,=\,0.1 including $X$\,=\,$10^{-2}, 10^{-3}, 10^{-4}, 10^{-5}, 10^{-6}$ that can improve stellar models of hydrogen deficient stars. We implement these new OPLIB radiative opacity tables in \MESA, and find that calibrated solar models agree broadly with previously published helioseismic and solar neutrino results. We find differences between using the new 1194 OPLIB opacity tables and the 126 OPAL opacity tables range from $\approx$\,20--80\% across individual chemical mixtures, up to $\approx$\,8\% and $\approx$\,15\% at the bottom and top of the solar convection zone respectively, and $\approx$\,7\% in the solar core. We also find differences between standard solar models using different opacity table sources that are on par with altering the initial abundance mixture. We conclude that this new, open-access set of OPLIB opacity tables does not solve the solar modeling problem, and suggest the investigation of physical mechanisms other than the atomic radiative opacity. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Comments: 28 pages, 12 (13) figures. Accepted for Publication in ApJ

Photons from neutrinos: the gamma ray echo of a supernova neutrino burst

Authors: Cecilia Lunardini, Joshua Loeffler, Mainak Mukhopadhyay, Matthew J. Hurley, Ebraheem Farag, F. X. Timmes

Abstract: When a star undergoes core collapse, a vast amount of energy is released in a ~10 s long burst of neutrinos of all species. Inverse beta decay in the star's hydrogen envelope causes an electromagnetic cascade which ultimately results in a flare of gamma rays - an "echo" of the neutrino burst - at the characteristic energy of 0.511 MeV. We study the phenomenology and detectability of this flare. It… ▽ More When a star undergoes core collapse, a vast amount of energy is released in a ~10 s long burst of neutrinos of all species. Inverse beta decay in the star's hydrogen envelope causes an electromagnetic cascade which ultimately results in a flare of gamma rays - an "echo" of the neutrino burst - at the characteristic energy of 0.511 MeV. We study the phenomenology and detectability of this flare. Its luminosity curve is characterized by a fast, seconds-long, rise and an equally fast decline, with a minute- or hour-long plateau in between. For a near-Earth star (distance D<1 kpc) the echo will be observable at near future gamma ray telescopes with an effective area of 10^3 cm^2 or larger. Its observation will inform us on the envelope size and composition. In conjunction with the direct detection of the neutrino burst, it will also give information on the neutrino emission away from the line of sight and will enable tests of neutrino propagation effects between the stellar surface and Earth. △ Less

Submitted 9 July, 2024; v1 submitted 20 December, 2023; originally announced December 2023.

Comments: 10 pages, 5 figures and 1 table. Description of the microphysics has been expanded; minor edits were made for clarity. This version matches the one to be published in the Astrophysical Journal

Stellar Neutrino Emission Across The Mass-Metallicity Plane

Authors: Ebraheem Farag, F. X. Timmes, Morgan T. Chidester, Samalka Anandagoda, Dieter H. Hartmann

Abstract: We explore neutrino emission from nonrotating, single star models across six initial metallicities and seventy initial masses from the zero-age main sequence to the final fate. Overall, across the mass spectrum, we find metal-poor stellar models tend to have denser, hotter and more massive cores with lower envelope opacities, larger surface luminosities, and larger effective temperatures than thei… ▽ More We explore neutrino emission from nonrotating, single star models across six initial metallicities and seventy initial masses from the zero-age main sequence to the final fate. Overall, across the mass spectrum, we find metal-poor stellar models tend to have denser, hotter and more massive cores with lower envelope opacities, larger surface luminosities, and larger effective temperatures than their metal-rich counterparts. Across the mass-metallicity plane we identify the sequence (initial CNO $\rightarrow$ $^{14}$N $\rightarrow$ $^{22}$Ne $\rightarrow$ $^{25}$Mg $\rightarrow$ $^{26}$Al $\rightarrow$ $^{26}$Mg $\rightarrow$ $^{30}$P $\rightarrow$ $^{30}$Si) as making primary contributions to the neutrino luminosity at different phases of evolution. For the low-mass models we find neutrino emission from the nitrogen flash and thermal pulse phases of evolution depend strongly on the initial metallicity. For the high-mass models, neutrino emission at He-core ignition and He-shell burning depends strongly on the initial metallicity. Anti-neutrino emission during C, Ne, and O burning shows a strong metallicity dependence with $^{22}$Ne($α$,$n$)$^{25}$Mg providing much of the neutron excess available for inverse-$β$ decays. We integrate the stellar tracks over an initial mass function and time to investigate the neutrino emission from a simple stellar population. We find average neutrino emission from simple stellar populations to be 0.5--1.2 MeV electron neutrinos. Lower metallicity stellar populations produce slightly larger neutrino luminosities and average $β$ decay energies. This study can provide targets for neutrino detectors from individual stars and stellar populations. We provide convenient fitting formulae and open access to the photon and neutrino tracks for more sophisticated population synthesis models. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: 37 pages, 23 figures, accepted for publication in ApJs

MESA-Web: A cloud resource for stellar evolution in astronomy curricula

Authors: Carl E. Fields, Richard H. D. Townsend, A. L. Dotter, Michael Zingale, F. X. Timmes

Abstract: We present MESA-Web, a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MES… ▽ More We present MESA-Web, a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MESA-Web can be used as an educational tool for stars in the classroom or for scientific investigations. We report on new capabilities of MESA-Web introduced since its 2015 release including learner-supplied nuclear reaction rates, custom stopping conditions, and an expanded selection of input parameters. To foster collaboration we have created a Zenodo MESA-Web community hub where instructors can openly share examples of using MESA-Web in the classroom. We discuss two examples in the current community hub. The first example is a lesson module on Red Giant Branch stars that includes a suite of exercises designed to fit a range of learners and a Jupyter workbook for additional analysis. The second example is lesson materials for an upper-level Astronomy majors course in Stars and Radiation that includes an assignment verifying some of the expected trends that are presented in a popular stellar physics textbook. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: Accepted for publication in the Astronomy Education Journal. We encourage users of MESA-Web to share their educational materials with the larger community using our Zenodo hub at https://zenodo.org/communities/mesa-web/

Seismic Signatures of the $^{12}$C($α$, $γ$)$^{16}$O Reaction Rate in White Dwarf Models with Overshooting

Authors: Morgan T. Chidester, F. X. Timmes, Ebraheem Farag

Abstract: We consider the combined effects that overshooting and the $^{12}$C($α$, $γ$)$^{16}$O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution. These signatures ho… ▽ More We consider the combined effects that overshooting and the $^{12}$C($α$, $γ$)$^{16}$O reaction rate have on variable white dwarf stellar models. We find that carbon-oxygen white dwarf models continue to yield pulsation signatures of the current experimental $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution. These signatures hold because the resonating mantle region, encompassing $\simeq$\,0.2\,\Msun\ in a typical $\simeq$\,0.6\,\Msun\ white dwarf model, still undergoes radiative helium burning during the evolution to a white dwarf. Our specific models show two potential low-order adiabatic g-modes, $g_2$ and $g_6$, that signalize the $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function. Both g-mode signatures induce average relative period shifts of $ΔP/P = 0.44 \%$ and $ΔP/P = 1.33\%$ for $g_2$ and $g_6$ respectively. We find that $g_6$ is a trapped mode, and the $g_2$ period signature is inversely proportional to the $^{12}$C($α$, $γ$)$^{16}$O reaction rate. The $g_6$ period signature generally separates the slower and faster reaction rates, and has a maximum relative period shift of $ΔP/P = 3.45\%$. We conclude that low-order g-mode periods from carbon-oxygen white dwarfs may still serve as viable probes for the $^{12}$C($α$, $γ$)$^{16}$O reaction rate probability distribution function when overshooting is included in the evolution. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: 18 pages, 10 total figures, 2 online interactive figures, zenodo data set link included. Accepted to the Astrophysical Journal

A Highly Magnified Star at Redshift 6.2

Authors: Brian Welch, Dan Coe, Jose M. Diego, Adi Zitrin, Erik Zackrisson, Paola Dimauro, Yolanda Jimenez-Teja, Patrick Kelly, Guillaume Mahler, Masamune Oguri, F. X. Timmes, Rogier Windhorst, Michael Florian, S. E. DeMink, Roberto J. Avila, Jay Anderson, Larry Bradley, Keren Sharon, Anton Vikaeus, Stephan McCandliss, Marusa Bradac, Jane Rigby, Brenda Frye, Sune Toft, Victoria Strait , et al. (4 additional authors not shown)

Abstract: Galaxy clusters magnify background objects through strong gravitational lensing. Typical magnifications for lensed galaxies are factors of a few but can also be as high as tens or hundreds, stretching galaxies into giant arcs. Individual stars can attain even higher magnifications given fortuitous alignment with the lensing cluster. Recently, several individual stars at redshift $z \sim 1 - 1.5$ h… ▽ More Galaxy clusters magnify background objects through strong gravitational lensing. Typical magnifications for lensed galaxies are factors of a few but can also be as high as tens or hundreds, stretching galaxies into giant arcs. Individual stars can attain even higher magnifications given fortuitous alignment with the lensing cluster. Recently, several individual stars at redshift $z \sim 1 - 1.5$ have been discovered, magnified by factors of thousands, temporarily boosted by microlensing. Here we report observations of a more distant and persistent magnified star at redshift $z_{\rm phot} = 6.2 \pm 0.1$, 900 Myr after the Big Bang. This star is magnified by a factor of thousands by the foreground galaxy cluster lens WHL0137--08 ($z = 0.566$), as estimated by four independent lens models. Unlike previous lensed stars, the magnification and observed brightness (AB mag 27.2) have remained roughly constant over 3.5 years of imaging and follow-up. The delensed absolute UV magnitude $M_{UV} = -10 \pm 2$ is consistent with a star of mass $M > 50 M_{\odot}$. Confirmation and spectral classification are forthcoming from approved observations with the James Webb Space Telescope △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: 50 pages, 11 figures (3 main text, 8 extended data). Published in Nature

Journal ref: Nature 603, 815-818 (2022)

Flash-X, a multiphysics simulation software instrument

Authors: Anshu Dubey, Klaus Weide, Jared O'Neal, Akash Dhruv, Sean Couch, J. Austin Harris, Tom Klosterman, Rajeev Jain, Johann Rudi, Bronson Messer, Michael Pajkos, Jared Carlson, Ran Chu, Mohamed Wahib, Saurabh Chawdhary, Paul M. Ricker, Dongwook Lee, Katie Antypas, Katherine M. Riley, Christopher Daley, Murali Ganapathy, Francis X. Timmes, Dean M. Townsley, Marcos Vanella, John Bachan , et al. (6 additional authors not shown)

Abstract: Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Fla… ▽ More Flash-X is a highly composable multiphysics software system that can be used to simulate physical phenomena in several scientific domains. It derives some of its solvers from FLASH, which was first released in 2000. Flash-X has a new framework that relies on abstractions and asynchronous communications for performance portability across a range of increasingly heterogeneous hardware platforms. Flash-X is meant primarily for solving Eulerian formulations of applications with compressible and/or incompressible reactive flows. It also has a built-in, versatile Lagrangian framework that can be used in many different ways, including implementing tracers, particle-in-cell simulations, and immersed boundary methods. △ Less

Submitted 24 August, 2022; originally announced August 2022.

Comments: 16 pages, 5 Figures, published open access in SoftwareX

Journal ref: SoftwareX, Volume 19, 2022, 101168,ISSN 2352-7110

Resolving The Peak Of The Black Hole Mass Spectrum

Authors: Ebraheem Farag, Mathieu Renzo, Robert Farmer, Morgan T. Chidester, F. X. Timmes

Abstract: Gravitational wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernova (PISN). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections from merging BHs. We use \MESA\ to evolve single, no… ▽ More Gravitational wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernova (PISN). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections from merging BHs. We use \MESA\ to evolve single, non-rotating, massive helium cores with a metallicity of $Z = 10^{-5}$ until they either collapse to form a BH or explode as a PISN without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the $\pm 3σ$ uncertainty in our high resolution tabulated $^{12}$C($α$,$γ$)$^{16}$O reaction rate probability distribution function. We extensively test the temporal and mass resolution to resolve the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower resolution models reach convergence. We establish a new lower edge of the upper mass gap as M\textsubscript{lower} $\simeq$\,60$^{+32}_{-14}$\,\Msun\ from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16}\text{O}$ rate. We explore the effect of a larger 3-$α$ rate on the lower edge of the upper mass gap, finding M\textsubscript{lower} $\simeq$\,69$^{+34}_{-18}$\,\Msun. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog. △ Less

Submitted 20 August, 2022; originally announced August 2022.

Comments: 22 pages, 10 figures, accepted for publication in the astrophysical journal

JWST Imaging of Earendel, the Extremely Magnified Star at Redshift $z=6.2$

Authors: Brian Welch, Dan Coe, Erik Zackrisson, S. E. de Mink, Swara Ravindranath, Jay Anderson, Gabriel Brammer, Larry Bradley, Jinmi Yoon, Patrick Kelly, Jose M. Diego, Rogier Windhorst, Adi Zitrin, Paola Dimauro, Yolanda Jimenez-Teja, Abdurro'uf, Mario Nonino, Ana Acebron, Felipe Andrade-Santos, Roberto J. Avila, Matthew B. Bayliss, Alex Benitez, Tom Broadhurst, Rachana Bhatawdekar, Marusa Bradac , et al. (38 additional authors not shown)

Abstract: The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μ$m. In these higher resolution images, Earendel remains a single unresolved point… ▽ More The gravitationally lensed star WHL0137-LS, nicknamed Earendel, was identified with a photometric redshift $z_{phot} = 6.2 \pm 0.1$ based on images taken with the Hubble Space Telescope. Here we present James Webb Space Telescope (JWST) Near Infrared Camera (NIRCam) images of Earendel in 8 filters spanning 0.8--5.0$μ$m. In these higher resolution images, Earendel remains a single unresolved point source on the lensing critical curve, increasing the lower limit on the lensing magnification to $μ> 4000$ and restricting the source plane radius further to $r < 0.02$ pc, or $\sim 4000$ AU. These new observations strengthen the conclusion that Earendel is best explained by an individual star or multiple star system, and support the previous photometric redshift estimate. Fitting grids of stellar spectra to our photometry yields a stellar temperature of $T_{\mathrm{eff}} \simeq 13000$--16000 K assuming the light is dominated by a single star. The delensed bolometric luminosity in this case ranges from $\log(L) = 5.8$--6.6 $L_{\odot}$, which is in the range where one expects luminous blue variable stars. Follow-up observations, including JWST NIRSpec scheduled for late 2022, are needed to further unravel the nature of this object, which presents a unique opportunity to study massive stars in the first billion years of the universe. △ Less

Submitted 9 November, 2022; v1 submitted 18 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL. Data products, lens models, and analysis code will be available online at https://cosmic-spring.github.io

Modules for Experiments in Stellar Astrophysics (MESA): Time-Dependent Convection, Energy Conservation, Automatic Differentiation, and Infrastructure

Authors: Adam S. Jermyn, Evan B. Bauer, Josiah Schwab, R. Farmer, Warrick H. Ball, Earl P. Bellinger, Aaron Dotter, Meridith Joyce, Pablo Marchant, Joey S. G. Mombarg, William M. Wolf, Tin Long Sunny Wong, Giulia C. Cinquegrana, Eoin Farrell, R. Smolec, Anne Thoul, Matteo Cantiello, Falk Herwig, Odette Toloza, Lars Bildsten, Richard H. D. Townsend, F. X. Timmes

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MES… ▽ More We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). The new auto_diff module implements automatic differentiation in MESA, an enabling capability that alleviates the need for hard-coded analytic expressions or finite difference approximations. We significantly enhance the treatment of the growth and decay of convection in MESA with a new model for time-dependent convection, which is particularly important during late-stage nuclear burning in massive stars and electron degenerate ignition events. We strengthen MESA's implementation of the equation of state, and we quantify continued improvements to energy accounting and solver accuracy through a discussion of different energy equation features and enhancements. To improve the modeling of stars in MESA we describe key updates to the treatment of stellar atmospheres, molecular opacities, Compton opacities, conductive opacities, element diffusion coefficients, and nuclear reaction rates. We introduce treatments of starspots, an important consideration for low-mass stars, and modifications for superadiabatic convection in radiation-dominated regions. We describe new approaches for increasing the efficiency of calculating monochromatic opacities and radiative levitation, and for increasing the efficiency of evolving the late stages of massive stars with a new operator split nuclear burning mode. We close by discussing major updates to MESA's software infrastructure that enhance source code development and community engagement. △ Less

Submitted 30 December, 2022; v1 submitted 7 August, 2022; originally announced August 2022.

Comments: 50 pages, 29 figures; Accepted to ApJS

On Trapped Modes In Variable White Dwarfs As Probes Of The $^{12}$C($α, γ$)$^{16}$O Reaction Rate

Authors: Morgan T. Chidester, Ebraheem Farag, F. X. Timmes

Abstract: We seek signatures of the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it f… ▽ More We seek signatures of the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function in the pulsation periods of carbon-oxygen white dwarf models. We find that adiabatic g-modes trapped by the interior carbon-rich layer offer potentially useful signatures of this reaction rate probability distribution function. Probing the carbon-rich region is relevant because it forms during the evolution of low-mass stars under radiative helium burning conditions, mitigating the impact of convective mixing processes. We make direct quantitative connections between the pulsation periods of the identified trapped g-modes in variable WD models and the current experimental $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function. We find an average spread in relative period shifts of $ΔP/P \simeq \pm$ 2\% for the identified trapped g-modes over the $\pm$ 3$σ$ uncertainty in the $^{12}$C$(α,γ)^{16}$O reaction rate probability distribution function -- across the effective temperature range of observed DAV and DBV white dwarfs and for different white dwarf masses, helium shell masses, and hydrogen shell masses. The g-mode pulsation periods of observed white dwarfs are typically given to 6-7 significant figures of precision. This suggests that an astrophysical constraint on the $^{12}$C$(α,γ)^{16}$O reaction rate could, in principle, be extractable from the period spectrum of observed variable white dwarfs. △ Less

Submitted 5 July, 2022; originally announced July 2022.

Comments: This manuscript has been accepted for publication to the Astrophysical Journal. 13 pages, 9 figures. Presented at the AAS 240th summer conference in Pasadena. Zenodo link for our figures, data, etc. found here: 10.5281/zenodo.6779983

Horizons: Nuclear Astrophysics in the 2020s and Beyond

Authors: H. Schatz, A. D. Becerril Reyes, A. Best, E. F. Brown, K. Chatziioannou, K. A. Chipps, C. M. Deibel, R. Ezzeddine, D. K. Galloway, C. J. Hansen, F. Herwig, A. P. Ji, M. Lugaro, Z. Meisel, D. Norman, J. S. Read, L. F. Roberts, A. Spyrou, I. Tews, F. X. Timmes, C. Travaglio, N. Vassh, C. Abia, P. Adsley, S. Agarwal , et al. (140 additional authors not shown)

Abstract: Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilit… ▽ More Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 96 pages. Submitted to Journal of Physics G

Report number: LA-UR-22-23997

Localised thermonuclear bursts from accreting magnetic white dwarfs

Authors: S. Scaringi, P. J. Groot, C. Knigge, A. J. Bird, E. Breedt, D. A. H. Buckley, Y. Cavecchi, N. D. Degenaar, D. de Martino, C. Done, M. Fratta, K. Ilkiewicz, E. Koerding, J. -P. Lasota, C. Littlefield, C. F. Manara, M. O'Brien, P. Szkody, F. X. Timmes

Abstract: Nova explosions are caused by global thermonuclear runaways triggered in the surface layers of accreting white dwarfs. It has been predicted that localised thermonuclear bursts on white dwarfs can also take place, similar to Type I X-ray bursts observed in accreting neutron stars. Unexplained rapid bursts from the binary system TV Columbae, in which mass is accreted onto a moderately-strong magnet… ▽ More Nova explosions are caused by global thermonuclear runaways triggered in the surface layers of accreting white dwarfs. It has been predicted that localised thermonuclear bursts on white dwarfs can also take place, similar to Type I X-ray bursts observed in accreting neutron stars. Unexplained rapid bursts from the binary system TV Columbae, in which mass is accreted onto a moderately-strong magnetised white dwarf from a low-mass companion, have been observed on several occasions in the past $\approx40$ years. During these bursts the optical/UV luminosity increases by a factor of $>3$ in less than an hour and fades over $\approx10$ hours. Fast outflows have been observed in UV spectral lines, with velocities $>3500$ km s$^{-1}$, comparable to the escape velocity from the white dwarf surface. Here we report on optical bursts observed in TV Columbae as well as in two additional accreting systems, EI Ursae Majoris and ASASSN-19bh. The bursts have a total energy $\approx~10^{-6}$ those of classical nova explosions ("micronovae"), and bear a strong resemblance to Type I X-ray bursts. We exclude accretion or stellar magnetic reconnection events as their origin and suggest thermonuclear runaway events in magnetically-confined accretion columns as a viable explanation. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 11 pages, 7 figures. Submitted on 4 October 2021. Accepted for publication in Nature on 1 February 2022

Stable nickel production in Type Ia supernovae: A smoking gun for the progenitor mass?

Authors: Stéphane Blondin, Eduardo Bravo, Frank Timmes, Luc Dessart, D. John Hillier

Abstract: At present, there are strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit (MCh ~ 1.4 Msun) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between MCh and sub-MCh events. In particular, it is thought that the higher-density… ▽ More At present, there are strong indications that white dwarf (WD) stars with masses well below the Chandrasekhar limit (MCh ~ 1.4 Msun) contribute a significant fraction of SN Ia progenitors. The relative fraction of stable iron-group elements synthesized in the explosion has been suggested as a possible discriminant between MCh and sub-MCh events. In particular, it is thought that the higher-density ejecta of MCh WDs, which favours the synthesis of stable isotopes of nickel, results in prominent [Ni II] lines in late-time spectra. We study the explosive nucleosynthesis of stable nickel in SNe Ia resulting from MCh and sub-MCh progenitors. We explore the potential for lines of [Ni II] at 7378 Å and 1.94 microns in late-time spectra to serve as a diagnostic of the exploding WD mass, using nonlocal thermodynamic equilibrium radiative-transfer simulations with the CMFGEN code. We find that the radiative proton-capture reaction 57Co(p,gamma)58Ni is the dominant production mode for 58Ni in both MCh and sub-MCh models, while the alpha-capture reaction on 54Fe has a negligible impact on the final 58Ni yield. More importantly, we demonstrate that the lack of [Ni II] lines in late-time spectra of sub-MCh events is not always due to an under-abundance of stable Ni; rather, it results from the higher ionization of Ni in the inner ejecta. Conversely, the strong [Ni II] lines predicted in our 1D MCh models are completely suppressed when 56Ni is sufficiently mixed with the innermost layers, which are rich in stable iron-group elements. [Ni II] lines in late-time SN Ia spectra have a complex dependency on the abundance of stable Ni, which limits their use in distinguishing among MCh and sub-MCh progenitors. However, we argue that a low-luminosity SN Ia displaying strong [Ni II] lines would most likely result from a Chandrasekhar-mass progenitor. [Abridged] △ Less

Submitted 23 March, 2022; v1 submitted 28 September, 2021; originally announced September 2021.

Comments: Accepted for publication in A&A, replaced with accepted version (+ corrected a typo in the conclusions: "overabundance" replaced with "over abundance"). 20 pages, 10 figures. Model spectra available at https://zenodo.org/record/5528088

Journal ref: A&A 660, A96 (2022)

Observing intermediate-mass black holes and the upper--stellar-mass gap with LIGO and Virgo

Authors: Ajit Kumar Mehta, Alessandra Buonanno, Jonathan Gair, M. Coleman Miller, Ebraheem Farag, R. J. deBoer, M. Wiescher, F. X. Timmes

Abstract: Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses… ▽ More Using ground-based gravitational-wave detectors, we probe the mass function of intermediate-mass black holes (IMBHs) wherein we also include BHs in the upper mass gap $\sim 60-130~M_\odot$. Employing the projected sensitivity of the upcoming LIGO and Virgo fourth observing (O4) run, we perform Bayesian analysis on quasi-circular non-precessing, spinning IMBH binaries (IMBHBs) with total masses $50\mbox{--} 500\, M_\odot$, mass ratios 1.25, 4, and 10, and dimensionless spins up to 0.95, and estimate the precision with which the source-frame parameters can be measured. We find that, at $2σ$, the mass of the heavier component of IMBHBs can be constrained with an uncertainty of $\sim 10-40\%$ at a signal-to-noise ratio of $20$. Focusing on the stellar-mass gap with new tabulations of the $^{12}\text{C}(α, γ)^{16} \text{O}$ reaction rate and its uncertanties, we evolve massive helium core stars using \MESA\, to establish the lower and upper edge of the mass gap as $\simeq$\,59$^{+34}_{-13}$\,$M_{\odot}$ and $\simeq$\,139$^{+30}_{-14}$\,$M_{\odot}$ respectively, where the error bars give the mass range that follows from the $\pm 3σ$ uncertainty in the $^{12}\text{C}(α, γ) ^{16} \text{O}$ nuclear reaction rate. We find that high resolution of the tabulated reaction rate and fine temporal resolution are necessary to resolve the peak of the BH mass spectrum. We then study IMBHBs with components lying in the mass gap and show that the O4 run will be able to robustly identify most such systems. Finally, we re-analyse GW190521 with a state-of-the-art aligned-spin waveform model, finding that the primary mass lies in the mass gap with 90\% credibility. △ Less

Submitted 19 October, 2021; v1 submitted 13 May, 2021; originally announced May 2021.

Comments: 29 pages, 18 figures

Journal ref: The Astrophysical Journal, Volume 924, Number 1 (2022)

Skye: A Differentiable Equation of State

Authors: Adam S. Jermyn, Josiah Schwab, Evan Bauer, F. X. Timmes, Alexander Y. Potekhin

Abstract: Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Cou… ▽ More Stellar evolution and numerical hydrodynamics simulations depend critically on access to fast, accurate, thermodynamically consistent equations of state. We present Skye, a new equation of state for fully-ionized matter. Skye includes the effects of positrons, relativity, electron degeneracy, Coulomb interactions, non-linear mixing effects, and quantum corrections. Skye determines the point of Coulomb crystallization in a self-consistent manner, accounting for mixing and composition effects automatically. A defining feature of this equation of state is that it uses analytic free energy terms and provides thermodynamic quantities using automatic differentiation machinery. Because of this, Skye is easily extended to include new effects by simply writing new terms in the free energy. We also introduce a novel thermodynamic extrapolation scheme for extending analytic fits to the free energy beyond the range of the fitting data while preserving desirable properties like positive entropy and sound speed. We demonstrate Skye in action in the MESA stellar evolution software instrument by computing white dwarf cooling curves. △ Less

Submitted 1 April, 2021; originally announced April 2021.

Comments: Accepted in ApJ. 27 pages, 19 figures

Journal ref: Astrophys. J., 913, 72 (2021)

On The Impact Of 22Ne On The Pulsation Periods Of Carbon-Oxygen White Dwarfs With Helium Dominated Atmospheres

Authors: Morgan T. Chidester, F. X. Timmes, Josiah Schwab, Richard H. D. Townsend, Ebraheem Farag, Anne Thoul, C. E. Fields, Evan B. Bauer, Michael H. Montgomery

Abstract: We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models wi… ▽ More We explore changes in the adiabatic low-order g-mode pulsation periods of 0.526, 0.560, and 0.729 M$_\odot$ carbon-oxygen white dwarf models with helium-dominated envelopes due to the presence, absence, and enhancement of $^{22}$Ne in the interior. The observed g-mode pulsation periods of such white dwarfs are typically given to 6$-$7 significant figures of precision. Usually white dwarf models without $^{22}$Ne are fit to the observed periods and other properties. The root-mean-square residuals to the $\simeq$ 150$-$400 s low-order g-mode periods are typically in the range of $σ_{rms}$ $\lesssim$ 0.3 s, for a fit precision of $σ_{rms}/ P$ $\lesssim$ 0.3 %. We find average relative period shifts of $ΔP/P$ $\simeq$ $\pm$ 0.5 % for the low-order dipole and quadrupole g-mode pulsations within the observed effective temperature window, with the range of $ΔP/P$ depending on the specific g-mode, abundance of $^{22}$Ne, effective temperature, and mass of the white dwarf model. This finding suggests a systematic offset may be present in the fitting process of specific white dwarfs when $^{22}$Ne is absent. As part of the fitting processes involves adjusting the composition profiles of a white dwarf model, our study on the impact of $^{22}$Ne can provide new inferences on the derived interior mass fraction profiles. We encourage routinely including $^{22}$Ne mass fraction profiles, informed by stellar evolution models, to future generations of white dwarf model fitting processes. △ Less

Submitted 20 January, 2021; originally announced January 2021.

Comments: 23 pages, 19 figures total. Accepted to the ApJ. Presented results at 2021 AAS January meeting

Presupernova neutrinos: directional sensitivity and prospects for progenitor identification

Authors: Mainak Mukhopadhyay, Cecilia Lunardini, F. X. Timmes, Kai Zuber

Abstract: We explore the potential of current and future liquid scintillator neutrino detectors of O (10) kt mass to localize a pre-supernova neutrino signal in the sky. In the hours preceding the core collapse of a nearby star (at distance D < 1 kpc), tens to hundreds of inverse beta decay events will be recorded, and their reconstructed topology in the detector can be used to estimate the direction to the… ▽ More We explore the potential of current and future liquid scintillator neutrino detectors of O (10) kt mass to localize a pre-supernova neutrino signal in the sky. In the hours preceding the core collapse of a nearby star (at distance D < 1 kpc), tens to hundreds of inverse beta decay events will be recorded, and their reconstructed topology in the detector can be used to estimate the direction to the star. Although the directionality of inverse beta decay is weak (~8% forward-backward asymmetry for currently available liquid scintillators), we find that for a fiducial signal of 200 events (which is realistic for Betelgeuse), a positional error of ~60 degrees can be achieved, resulting in the possibility to narrow the list of potential stellar candidates to less than ten, typically. For a configuration with improved forward-backward asymmetry (~40%, as expected for a lithium-loaded liquid scintillator), the angular sensitivity improves to ~15 degrees, and - when a distance upper limit is obtained from the overall event rate - it is in principle possible to uniquely identify the progenitor star. Any localization information accompanying an early supernova alert will be useful to multi-messenger observations and to particle physics tests using collapsing stars. △ Less

Submitted 6 October, 2020; v1 submitted 4 April, 2020; originally announced April 2020.

Comments: LaTeX, 16 pages, 9 figures, 7 tables. Updated to reflect published version. Minor corrections and clarifications, results unchanged

On Stellar Evolution In A Neutrino Hertzsprung-Russell Diagram

Authors: Ebraheem Farag, F. X. Timmes, Morgan Taylor, Kelly M. Patton, R. Farmer

Abstract: We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{ν,\odot}$ = 0.02398 $\cdot$ $L_{γ,\odot}$ = 9.1795 $\times$ 10… ▽ More We explore the evolution of a select grid of solar metallicity stellar models from their pre-main sequence phase to near their final fates in a neutrino Hertzsprung-Russell diagram, where the neutrino luminosity replaces the traditional photon luminosity. Using a calibrated \MESA\ solar model for the solar neutrino luminosity ($L_{ν,\odot}$ = 0.02398 $\cdot$ $L_{γ,\odot}$ = 9.1795 $\times$ 10$^{31}$ erg s$^{-1}$) as a normalization, we identify $\simeq$ 0.3 MeV electron neutrino emission from helium burning during the helium flash (peak $L_ν / L_{ν,\odot} \simeq$ 10$^4$, flux $Φ_{ν, {\rm He \ flash}} \simeq$ 170 (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$ for a star located at a distance of $d$ parsec, timescale $\simeq$ 3 days) and the thermal pulse (peak $L_ν / L_{ν,\odot} \simeq$ 10$^9$, flux $Φ_{ν, {\rm TP}} \simeq$ 1.7$\times$10$^7$ (10 pc/$d$)$^{2}$ cm$^{-2}$ s$^{-1}$, timescale $\simeq$ 0.1 yr) phases of evolution in low mass stars as potential probes for stellar neutrino astronomy. We also delineate the contribution of neutrinos from nuclear reactions and thermal processes to the total neutrino loss along the stellar tracks in a neutrino Hertzsprung-Russell diagram. We find, broadly but with exceptions, that neutrinos from nuclear reactions dominate whenever hydrogen and helium burn, and that neutrinos from thermal processes dominate otherwise. △ Less

Submitted 12 March, 2020; originally announced March 2020.

Comments: 11 pages, 5 figures, accepted for publication in the astrophysical journal

Laminar Flame Speeds in Degenerate Oxygen-Neon Mixtures

Authors: Josiah Schwab, R. Farmer, F. X. Timmes

Abstract: The collapse of degenerate oxygen-neon cores (i.e., electron-capture supernovae or accretion-induced collapse) proceeds through a phase in which a deflagration wave ("flame") forms at or near the center and propagates through the star. In models, the assumed speed of this flame influences whether this process leads to an explosion or to the formation of a neutron star. We calculate the laminar fla… ▽ More The collapse of degenerate oxygen-neon cores (i.e., electron-capture supernovae or accretion-induced collapse) proceeds through a phase in which a deflagration wave ("flame") forms at or near the center and propagates through the star. In models, the assumed speed of this flame influences whether this process leads to an explosion or to the formation of a neutron star. We calculate the laminar flame speeds in degenerate oxygen-neon mixtures with compositions motivated by detailed stellar evolution models. These mixtures include trace amounts of carbon and have a lower electron fraction than those considered in previous work. We find that trace carbon has little effect on the flame speeds, but that material with electron fraction $Y_e \approx 0.48-0.49$ has laminar flame speeds that are $\approx 2$ times faster than those at $Y_e = 0.5$. We provide tabulated flame speeds and a corresponding fitting function so that the impact of this difference can be assessed via full star hydrodynamical simulations of the collapse process. △ Less

Submitted 21 January, 2020; originally announced January 2020.

Comments: 12 pages, 12 figures; Accepted to ApJ

Digital Infrastructure in Astrophysics

Authors: Frank Timmes, Rich Townsend, Lars Bildsten

Abstract: Astronomy, as a field, has long encouraged the development of free, open digital infrastructure (e.g., National Research Council 2010, 2011). Examples range from simple scripts that enable individual scientific research, through software instruments for entire communities, to data reduction pipelines for telescope operations at national facilities. As with the digital infrastructure of our larger… ▽ More Astronomy, as a field, has long encouraged the development of free, open digital infrastructure (e.g., National Research Council 2010, 2011). Examples range from simple scripts that enable individual scientific research, through software instruments for entire communities, to data reduction pipelines for telescope operations at national facilities. As with the digital infrastructure of our larger society today (e.g., Eghbal 2016), nearly all astronomical research relies on free, open source software (FOSS) written and maintained by a small number of developers. And like the physical infrastructure of roads or bridges, digital infrastructure needs regular upkeep and maintenance (e.g., Eghbal 2016). In astronomy, financial support for maintaining existing digital infrastructure is generally much harder to secure than funding for developing new digital infrastructures that promise new science. Sustaining astronomy's digital infrastructure is a new topic for many, the sustainability challenges are not always widely known... △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 21 pages; report for the Ford and Sloan Foundation's Digital Infrastructure Research Program

On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits

Authors: Rogier A. Windhorst, M. Alpaslan, S. Andrews, T. Ashcraft, T. Broadhurst, D. Coe, S. Cohen, C. Conselice, J. Diego, M. Dijkstra, S. Driver, K. Duncan, S. Finkelstein, B. Frye, A. Griffiths, N. Grogin, N. Hathi, A. Hopkins, R. Jansen, B. Joshi, A. Kashlinsky, W. Keel, P. Kelly, D. Kim, A. Koekemoer , et al. (22 additional authors not shown)

Abstract: Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stella… ▽ More Recent near-IR power-spectra and panchromatic Extragalactic Background Light measurements provide upper limits on the near-IR surface brightness (SB>31 mag/arcsec^2) that may come from Pop III stars and accretion disks around resulting stellar-mass black holes (BHs) in the epoch of First Light (z=7-17). Physical parameters for zero metallicity Pop III stars at z>7 can be estimated from MESA stellar evolution models through helium-depletion, and for BH accretion disks from quasar microlensing results and multicolor accretion models. Second-generation stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions in their AGB stage. The near-IR SB constraints can be used to calculate the number of caustic transits behind lensing clusters that JWST and the 25~39 m ground-based telescopes may detect for both Pop III stars and stellar mass BH accretion disks. Because Pop III stars and stellar mass BH accretion disks have sizes of a few x 10^-11 arcsec at z>7, typical caustic magnifications can be mu=10^4~10^5, with rise times of hours and decline times of < 1 year for cluster transverse velocities of v<1000 km/s. Microlensing by intracluster medium objects can modify transit magnifications, and lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Pop III stars. To observe Pop III caustic transits directly may require monitoring 3~30 lensing clusters to AB< 29 mag over a decade or more. Such a program must be started with JWST in Cycle 1, and -- depending on the role of microlensing in the Intra Cluster Light -- should be continued for decades with the GMT and TMT, where JWST and the ground-based telescopes each will play a unique and strongly complementary role. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 6 pages + references, 3 figures, science white paper submitted to the Astro2020 decadal survey. arXiv admin note: substantial text overlap with arXiv:1801.03584

Cyberinfrastructure Requirements to Enhance Multi-messenger Astrophysics

Authors: Philip Chang, Gabrielle Allen, Warren Anderson, Federica B. Bianco, Joshua S. Bloom, Patrick R. Brady, Adam Brazier, S. Bradley Cenko, Sean M. Couch, Tyce DeYoung, Ewa Deelman, Zachariah B Etienne, Ryan J. Foley, Derek B Fox, V. Zach Golkhou, Darren R Grant, Chad Hanna, Kelly Holley-Bockelmann, D. Andrew Howell, E. A. Huerta, Margaret W. G. Johnson, Mario Juric, David L. Kaplan, Daniel S. Katz, Azadeh Keivani , et al. (17 additional authors not shown)

Abstract: The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up,… ▽ More The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 7 pages, astro2020 white paper

Near-Earth Supernova Explosions: Evidence, Implications, and Opportunities

Authors: Brian D. Fields, John R. Ellis, Walter R. Binns, Dieter Breitschwerdt, Georgia A. de Nolfo, Roland Diehl, Vikram V. Dwarkadas, Adrienne Ertel, Thomas Faestermann, Jenny Feige, Caroline Fitoussi, Priscilla Frisch, David Graham, Brian Haley, Alexander Heger, Wolfgang Hillebrandt, Martin H. Israel, Thomas Janka, Michael Kachelriess, Gunther Korschinek, Marco Limongi, Maria Lugaro, Franciole Marinho, Adrian Melott, Richard A. Mewaldt , et al. (14 additional authors not shown)

Abstract: There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact o… ▽ More There is now solid experimental evidence of at least one supernova explosion within 100 pc of Earth within the last few million years, from measurements of the short-lived isotope 60Fe in widespread deep-ocean samples, as well as in the lunar regolith and cosmic rays. This is the first established example of a specific dated astrophysical event outside the Solar System having a measurable impact on the Earth, offering new probes of stellar evolution, nuclear astrophysics, the astrophysics of the solar neighborhood, cosmic-ray sources and acceleration, multi-messenger astronomy, and astrobiology. Interdisciplinary connections reach broadly to include heliophysics, geology, and evolutionary biology. Objectives for the future include pinning down the nature and location of the established near-Earth supernova explosions, seeking evidence for others, and searching for other short-lived isotopes such as 26Al and 244Pu. The unique information provided by geological and lunar detections of radioactive 60Fe to assess nearby supernova explosions make now a compelling time for the astronomy community to advocate for supporting multi-disciplinary, cross-cutting research programs. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 11 pages, 2 figures. Astro2020 Science White Paper submitted to the 2020 Decadal Survey on Astronomy and Astrophysics

Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation

Authors: Bill Paxton, R. Smolec, Josiah Schwab, A. Gautschy, Lars Bildsten, Matteo Cantiello, Aaron Dotter, R. Farmer, Jared A. Goldberg, Adam S. Jermyn, S. M. Kanbur, Pablo Marchant, Anne Thoul, Richard H. D. Townsend, William M. Wolf, Michael Zhang, F. X. Timmes

Abstract: We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For exam… ▽ More We update the capabilities of the open-knowledge software instrument Modules for Experiments in Stellar Astrophysics (MESA). RSP is a new functionality in MESAstar that models the non-linear radial stellar pulsations that characterize RR Lyrae, Cepheids, and other classes of variable stars. We significantly enhance numerical energy conservation capabilities, including during mass changes. For example, this enables calculations through the He flash that conserve energy to better than 0.001 %. To improve the modeling of rotating stars in MESA, we introduce a new approach to modifying the pressure and temperature equations of stellar structure, and a formulation of the projection effects of gravity darkening. A new scheme for tracking convective boundaries yields reliable values of the convective-core mass, and allows the natural emergence of adiabatic semiconvection regions during both core hydrogen- and helium-burning phases. We quantify the parallel performance of MESA on current generation multicore architectures and demonstrate improvements in the computational efficiency of radiative levitation. We report updates to the equation of state and nuclear reaction physics modules. We briefly discuss the current treatment of fallback in core-collapse supernova models and the thermodynamic evolution of supernova explosions. We close by discussing the new MESA Testhub software infrastructure to enhance source-code development. △ Less

Submitted 16 May, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: 57 pages, 57 figures; Accepted to ApJS

Catching Element Formation In The Act

Authors: Chris L. Fryer, Frank Timmes, Aimee L. Hungerford, Aaron Couture, Fred Adams, Wako Aoki, Almudena Arcones, David Arnett, Katie Auchettl, Melina Avila, Carles Badenes, Eddie Baron, Andreas Bauswein, John Beacom, Jeff Blackmon, Stephane Blondin, Peter Bloser, Steve Boggs, Alan Boss, Terri Brandt, Eduardo Bravo, Ed Brown, Peter Brown, Steve Bruenn. Carl Budtz-Jorgensen, Eric Burns , et al. (194 additional authors not shown)

Abstract: Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-ray… ▽ More Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions. △ Less

Submitted 7 February, 2019; originally announced February 2019.

Comments: 14 pages including 3 figures

Report number: LA-UR-18-29748

Perspectives on Astrophysics Based on Atomic, Molecular, and Optical (AMO) Techniques

Authors: Daniel Wolf Savin, James F. Babb, Paul M. Bellan, Crystal Brogan, Jan Cami, Paola Caselli, Lia Corrales, Gerardo Dominguez, Steven R. Federman, Chris J. Fontes, Richard Freedman, Brad Gibson, Leon Golub, Thomas W. Gorczyca, Michael Hahn, Sarah M. Hörst, Reggie L. Hudson, Jeffrey Kuhn, James E. Lawler, Maurice A. Leutenegger, Joan P. Marler, Michael C. McCarthy, Brett A. McGuire, Stefanie N. Milam, Nicholas A. Murphy , et al. (13 additional authors not shown)

Abstract: About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides i… ▽ More About two generations ago, a large part of AMO science was dominated by experimental high energy collision studies and perturbative theoretical methods. Since then, AMO science has undergone a transition and is now dominated by quantum, ultracold, and ultrafast studies. But in the process, the field has passed over the complexity that lies between these two extremes. Most of the Universe resides in this intermediate region. We put forward that the next frontier for AMO science is to explore the AMO complexity that describes most of the Cosmos. △ Less

Submitted 14 November, 2018; originally announced November 2018.

Comments: White paper submission to the Decadal Assessment and Outlook Report on Atomic, Molecular, and Optical (AMO) Science (AMO 2020)

Common Envelope Evolution of Massive Stars

Authors: Paul M. Ricker, Frank X. Timmes, Ronald E. Taam, Ronald F. Webbink

Abstract: The discovery via gravitational waves of binary black hole systems with total masses greater than $60M_\odot$ has raised interesting questions for stellar evolution theory. Among the most promising formation channels for these systems is one involving a common envelope binary containing a low metallicity, core helium burning star with mass $\sim 80-90M_\odot$ and a black hole with mass… ▽ More The discovery via gravitational waves of binary black hole systems with total masses greater than $60M_\odot$ has raised interesting questions for stellar evolution theory. Among the most promising formation channels for these systems is one involving a common envelope binary containing a low metallicity, core helium burning star with mass $\sim 80-90M_\odot$ and a black hole with mass $\sim 30-40M_\odot$. For this channel to be viable, the common envelope binary must eject more than half the giant star's mass and reduce its orbital separation by as much as a factor of 80. We discuss issues faced in numerically simulating the common envelope evolution of such systems and present a 3D AMR simulation of the dynamical inspiral of a low-metallicity red supergiant with a massive black hole companion. △ Less

Submitted 8 November, 2018; originally announced November 2018.

Comments: 6 pages, 3 figures, oral contribution: IAU Symposium 346 "High Mass X-ray Binaries: illuminating the passage from massive binaries to merging compact objects", Vienna, Austria, 27-31 August 2018

Journal ref: Proc. IAU 14 (2018) 449-454

The Impact of White Dwarf Luminosity Profiles on Oscillation Frequencies

Authors: F. X. Timmes, Richard H. D. Townsend, Evan B. Bauer, Anne Thoul, C. E. Fields, William M. Wolf

Abstract: KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportion… ▽ More KIC 08626021 is a pulsating DB white dwarf of considerable recent interest, and first of its class to be extensively monitored by Kepler for its pulsation properties. Fitting the observed oscillation frequencies of KIC 08626021 to a model can yield insights into its otherwise-hidden internal structure. Template-based white dwarf models choose a luminosity profile where the luminosity is proportional to the enclosed mass, $L_r \propto M_r$, independent of the effective temperature $T_{\rm eff}$. Evolutionary models of young white dwarfs with $T_{\rm eff} \gtrsim$ 25,000 K suggest neutrino emission gives rise to luminosity profiles with $L_r$ $\not\propto$ $M_r$. We explore this contrast by comparing the oscillation frequencies between two nearly identical white dwarf models: one with an enforced $L_r \propto M_r$ luminosity profile and the other with a luminosity profile determined by the star's previous evolution history. We find the low order g-mode frequencies differ by up to $\simeq$ 70 $μ$Hz over the range of Kepler observations for KIC 08626021. This suggests that by neglecting the proper thermal structure of the star (e.g., accounting for the effect of plasmon neutrino losses), the model frequencies calculated by using an $L_r \propto M_r$ profile may have uncorrected, effectively-random errors at the level of tens of $μ$Hz. A mean frequency difference of 30 $μ$Hz, based on linearly extrapolating published results, suggests a template model uncertainty in the fit precision of $\simeq$ 12% in white dwarf mass, $\simeq$ 9% in the radius, and $\simeq$ 3% in the central oxygen mass fraction. △ Less

Submitted 8 October, 2018; originally announced October 2018.

Comments: 7 pages, 5 figures; Accepted for publication in ApJ Letters

Quantifying How Density Gradients and Front Curvature Affect Carbon Detonation Strength During Type Ia Supernovae

Authors: Broxton J. Miles, Dean M. Townsley, Ken J. Shen, F. X. Timmes, Kevin Moore

Abstract: Accurately reproducing the physics behind the detonations of Type Ia supernovae and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full star simulations in the present and near future. In the lower density r… ▽ More Accurately reproducing the physics behind the detonations of Type Ia supernovae and the resultant nucleosynthetic yields is important for interpreting observations of spectra and remnants. The scales of the processes involved span orders of magnitudes, making the problem computationally impossible to ever fully resolve in full star simulations in the present and near future. In the lower density regions of the star, the curvature of the detonation front will slow the detonation, affecting the production of intermediate mass elements. We find that shock strengthening due to the density gradient present in the outer layers of the progenitor is essential for understanding the nucleosynthesis there, with burning extending well below the density at which a steady-state detonation is extinct. We show that a complete reaction network is not sufficient to obtain physical detonations at high densities and modest resolution due to numerical mixing at the unresolved reaction front. At low densities, below 6$\times$10$^{5}$ g cm$^{-3}$, it is possible to achieve high enough resolution to separate the shock and the reaction region,and the abundance structure predicted by fully resolved quasi-steady-state calculations is obtained. For our best current benchmark yields, we utilize a method in which the unresolved portion of Lagrangian histories are reconstructed based on fully resolved quasi-steady-state detonation calculations. These computations demonstrate that under-resolved simulations agree approximately, $\sim$10\% in post-shock values of temperature, pressure, density, and abundances, with expected detonation structures sufficiently far from the under-resolved region, but that there is still room for some improvement in the treatment of subgrid reactions in the hydrodynamics to before better than 1$\%$ can be achieved at all densities. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: Submitted to ApJ

On the observability of individual Population III stars and their stellar-mass black hole accretion disks through cluster caustic transits

Authors: Rogier A. Windhorst, F. X. Timmes, J. Stuart B. Wyithe, Mehmet Alpaslan, Stephen K. Andrews, Daniel Coe, Jose M. Diego, Mark Dijkstra, Simon P. Driver, Patrick L. Kelly, Duho Kim

Abstract: We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter const… ▽ More We summarize panchromatic Extragalactic Background Light data to place upper limits on the integrated near-infrared surface brightness (SB) that may come from Population III stars and possible accretion disks around their stellar-mass black holes (BHs) in the epoch of First Light, broadly taken from z$\simeq$7-17. Theoretical predictions and recent near-infrared power-spectra provide tighter constraints on their sky-signal. We outline the physical properties of zero metallicity Population III stars from MESA stellar evolution models through helium-depletion and of BH accretion disks at z$\gtrsim$7. We assume that second-generation non-zero metallicity stars can form at higher multiplicity, so that BH accretion disks may be fed by Roche-lobe overflow from lower-mass companions. We use these near-infrared SB constraints to calculate the number of caustic transits behind lensing clusters that the James Webb Space Telescope and the next generation ground-based telescopes may observe for both Population III stars and their BH accretion disks. Typical caustic magnifications can be $μ$$\simeq$10$^4$-10$^5$, with rise times of hours and decline times of $\lesssim$1 year for cluster transverse velocities of $v_{T}$$\lesssim$1000 km s$^{-1}$. Microlensing by intracluster medium objects can modify transit magnifications, but lengthen visibility times. Depending on BH masses, accretion-disk radii and feeding efficiencies, stellar-mass BH accretion-disk caustic transits could outnumber those from Population III stars. To observe Population III caustic transits directly may require to monitor 3-30 lensing clusters to AB$\lesssim$29 mag over a decade. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: 53 pages, 5 figures, Accepted for publication in ApJS

The Impact of Nuclear Reaction Rate Uncertainties On The Evolution of Core-Collapse Supernova Progenitors

Authors: C. E. Fields, F. X. Timmes, R. Farmer, I. Petermann, William M. Wolf, S. M. Couch

Abstract: We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicit… ▽ More We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models 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 $\texttt{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 $^{14}$N$(p,γ)^{15}$O, triple-$α$, $^{12}$C$(α,γ)^{16}$O, $^{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 reaction rates 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 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 towards iron core-collapse. △ Less

Submitted 17 December, 2017; originally announced December 2017.

Comments: Accepted for publication in ApJS, 33 pages, 21 figures

Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions

Authors: Bill Paxton, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Sergei Blinnikov, Paul Duffell, R. Farmer, Jared A. Goldberg, Pablo Marchant, Elena Sorokina, Anne Thoul, Richard H. D. Townsend, F. X. Timmes

Abstract: We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white… ▽ More We update the capabilities of the software instrument Modules for Experiments in Stellar Astrophysics (MESA) and enhance its ease of use and availability. Our new approach to locating convective boundaries is consistent with the physics of convection, and yields reliable values of the convective core mass during both hydrogen and helium burning phases. Stars with $M<8\,{\rm M_\odot}$ become white dwarfs and cool to the point where the electrons are degenerate and the ions are strongly coupled, a realm now available to study with MESA due to improved treatments of element diffusion, latent heat release, and blending of equations of state. Studies of the final fates of massive stars are extended in MESA by our addition of an approximate Riemann solver that captures shocks and conserves energy to high accuracy during dynamic epochs. We also introduce a 1D capability for modeling the effects of Rayleigh-Taylor instabilities that, in combination with the coupling to a public version of the STELLA radiation transfer instrument, creates new avenues for exploring Type II supernovae properties. These capabilities are exhibited with exploratory models of pair-instability supernova, pulsational pair-instability supernova, and the formation of stellar mass black holes. The applicability of MESA is now widened by the capability of importing multi-dimensional hydrodynamic models into MESA. We close by introducing software modules for handling floating point exceptions and stellar model optimization, and four new software tools -- MESAWeb, MESA-Docker, pyMESA, and mesastar.org -- to enhance MESA's education and research impact. △ Less

Submitted 3 January, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

Comments: 64 pages, 61 figures; Accepted to AAS Journals

Neutrinos from beta processes in a presupernova: probing the isotopic evolution of a massive star

Authors: Kelly M. Patton, Cecilia Lunardini, Robert J. Farmer, F. X. Timmes

Abstract: We present a new calculation of the neutrino flux received at Earth from a massive star in the $\sim 24$ hours of evolution prior to its explosion as a supernova (presupernova). Using the stellar evolution code MESA, the neutrino emissivity in each flavor is calculated at many radial zones and time steps. In addition to thermal processes, neutrino production via beta processes is modeled in detail… ▽ More We present a new calculation of the neutrino flux received at Earth from a massive star in the $\sim 24$ hours of evolution prior to its explosion as a supernova (presupernova). Using the stellar evolution code MESA, the neutrino emissivity in each flavor is calculated at many radial zones and time steps. In addition to thermal processes, neutrino production via beta processes is modeled in detail, using a network of 204 isotopes. We find that the total produced $ν_{e}$ flux has a high energy spectrum tail, at $E \gtrsim 3 - 4$ MeV, which is mostly due to decay and electron capture on isotopes with $A = 50 - 60$. In a tentative window of observability of $E \gtrsim 0.5$ MeV and $t < 2$ hours pre-collapse, the contribution of beta processes to the $ν_{e}$ flux is at the level of $\sim90\%$ . For a star at $D=1$ kpc distance, a 17 kt liquid scintillator detector would typically observe several tens of events from a presupernova, of which up to $\sim 30\%$ due to beta processes. These processes dominate the signal at a liquid argon detector, thus greatly enhancing its sensitivity to a presupernova. △ Less

Submitted 6 September, 2017; originally announced September 2017.

Comments: 14 pages, 5 figures

Report number: INT-PUB-17-037

Constraining The Single-Degenerate Channel of Type Ia Supernovae With Stable Iron-Group Elements in SNR 3C 397

Authors: Pranav Dave, Rahul Kashyap, Robert Fisher, Frank Timmes, Dean Townsley, Chris Byrohl

Abstract: Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenito… ▽ More Recent Suzaku X-ray spectra of SNR 3C 397 indicate enhanced stable iron-group element abundances of Ni, Mn, Cr, and Fe. Seeking to address key questions about the progenitor and explosion mechanism of 3C 397, we compute nucleosynthetic yields from a suite of multidimensional hydrodynamics models in the near-Chandrasekhar mass, single-degenerate paradigm for supernova Type Ia. Varying the progenitor white dwarf internal structure, composition, ignition, and explosion mechanism, we find the best match to the observed iron-peak elements of 3C 397 are dense (central density $\ge$ 6$\times$10$^{9}$ g cm$^{-3}$), low-carbon white dwarfs that undergo a weak, centrally-ignited deflagration, followed by a subsequent detonation. The amount of $^{56}$Ni produced is consistent with a normal or bright normal supernova Type Ia. A pure deflagration of a centrally-ignited, low central density ($\simeq$ 2$\times$10$^{9}$ g cm$^{-3}$) progenitor white dwarf, frequently considered in the literature, is also found to produce good agreement with 3C 397 nucleosynthetic yields, but leads to a subluminous SN Ia event, in conflict with X-ray linewidth data. Additionally, in contrast to prior work which suggested a large super-solar metallicity for the white dwarf progenitor for SNR 3C 397, we find satisfactory agreement for solar and sub-solar metallicity progenitors. We discuss a range of implications our results have for the single-degenerate channel. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: 22 pages and 8 figures. Accepted for publication in The Astrophysical Journal

Observational evidence for high neutronization in supernova remnants: implications for Type Ia supernova progenitors

Authors: Héctor Martínez-Rodríguez, Carles Badenes, Hiroya Yamaguchi, Eduardo Bravo, F. X. Timmes, Broxton J. Miles, Dean M. Townsley, Anthony L. Piro, Hideyuki Mori, Brett Andrews, Sangwook Park

Abstract: The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using C… ▽ More The physical process whereby a carbon--oxygen white dwarf explodes as a Type Ia supernova (SN Ia) remains highly uncertain. The degree of neutronization in SN Ia ejecta holds clues to this process because it depends on the mass and the metallicity of the stellar progenitor, and on the thermodynamic history prior to the explosion. We report on a new method to determine ejecta neutronization using Ca and S lines in the X-ray spectra of Type Ia supernova remnants (SNRs). Applying this method to \textit{Suzaku} data of Tycho, Kepler, 3C 397 and G337.2$-$0.7 in the Milky Way, and N103B in the Large Magellanic Cloud, we find that the neutronization of the ejecta in N103B is comparable to that of Tycho and Kepler, which suggests that progenitor metallicity is not the only source of neutronization in SNe Ia. We then use a grid of SN Ia explosion models to infer the metallicities of the stellar progenitors of our SNRs. The implied metallicities of 3C 397, G337.2$-$0.7, and N103B are major outliers compared to the local stellar metallicity distribution functions, indicating that progenitor metallicity can be ruled out as the origin of neutronization for these SNRs. Although the relationship between ejecta neutronization and equivalent progenitor metallicity is subject to uncertainties stemming from the $^{12}$C$\,$+$^{16}$O reaction rate, which affects the Ca/S mass ratio, our main results are not sensitive to these details. △ Less

Submitted 28 June, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 12 pages, 7 figures, accepted for publication in the ApJ

Journal ref: The Astrophysical Journal, 843, 35 (2017)

On Variations Of Pre-Supernova Model Properties

Authors: R. Farmer, C. E. Fields, I. Petermann, Luc Dessart, M. Cantiello, B. Paxton, F. X. Timmes

Abstract: We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear b… ▽ More We explore the variation in single star 15-30 $\rm{M}_{\odot}$, non-rotating, solar metallicity, pre-supernova MESA models due to changes in the number of isotopes in a fully-coupled nuclear reaction network and adjustments in the mass resolution. Within this two-dimensional plane we quantitatively detail the range of core masses at various stages of evolution, mass locations of the main nuclear burning shells, electron fraction profiles, mass fraction profiles, burning lifetimes, stellar lifetimes, and compactness parameter at core-collapse for models with and without mass loss. Up to carbon burning we generally find mass resolution has a larger impact on the variations than the number of isotopes, while the number of isotopes plays a more significant role in determining the span of the variations for neon, oxygen and silicon burning. Choice of mass resolution dominates the variations in the structure of the intermediate convection zone and secondary convection zone during core and shell hydrogen burning respectively, where we find a minimum mass resolution of $\approx$ 0.01 $\rm{M}_{\odot}$ is necessary to achieve convergence in the helium core mass at the $\approx$5% level. On the other hand, at the onset of core-collapse we find $\approx$30% variations in the central electron fraction and mass locations of the main nuclear burning shells, a minimum of $\approx$127 isotopes is needed to attain convergence of these values at the $\approx$10% level. △ Less

Submitted 3 November, 2016; originally announced November 2016.

Comments: 32 pages, 1 table, 25 figures, Accepted ApJS

Bayesian Estimation of Thermonuclear Reaction Rates

Authors: Christian Iliadis, Kevin Anderson, Alain Coc, Frank Timmes, Sumner Starrfield

Abstract: The problem of estimating non-resonant astrophysical S-factors and thermonuclear reaction rates, based on measured nuclear cross sections, is of major interest for nuclear energy generation, neutrino physics, and element synthesis. Many different methods have been applied in the past to this problem, almost all of them based on traditional statistics. Bayesian methods, on the other hand, are now i… ▽ More The problem of estimating non-resonant astrophysical S-factors and thermonuclear reaction rates, based on measured nuclear cross sections, is of major interest for nuclear energy generation, neutrino physics, and element synthesis. Many different methods have been applied in the past to this problem, almost all of them based on traditional statistics. Bayesian methods, on the other hand, are now in widespread use in the physical sciences. In astronomy, for example, Bayesian statistics is applied to the observation of extra-solar planets, gravitational waves, and type Ia supernovae. However, nuclear physics, in particular, has been slow to adopt Bayesian methods. We present astrophysical S-factors and reaction rates based on Bayesian statistics. We develop a framework that incorporates robust parameter estimation, systematic effects, and non-Gaussian uncertainties in a consistent manner. The method is applied to the d(p,$γ$)$^3$He, $^3$He($^3$He,2p)$^4$He, and $^3$He($α$,$γ$)$^7$Be reactions, important for deuterium burning, solar neutrinos, and big bang nucleosynthesis. △ Less

Submitted 28 August, 2016; v1 submitted 20 August, 2016; originally announced August 2016.

Comments: 8 Figures

The Importance of Computation in Astronomy Education

Authors: M. Zingale, F. X. Timmes, R. Fisher, B. W. O'Shea

Abstract: Computational skills are required across all astronomy disciplines. Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research. We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergra… ▽ More Computational skills are required across all astronomy disciplines. Many students enter degree programs without sufficient skills to solve computational problems in their core classes or contribute immediately to research. We recommend advocacy for computational literacy, familiarity with fundamental software carpentry skills, and mastery of basic numerical methods by the completion of an undergraduate degree in Astronomy. We recommend the AAS Education Task Force advocate for a significant increase in computational literacy. We encourage the AAS to modestly fund efforts aimed at providing Open Education Resources (OER) that will significantly impact computational literacy in astronomy education. △ Less

Submitted 5 June, 2016; originally announced June 2016.

Comments: white paper submitted to the AAS Education Taskforce call (https://aas.org/posts/opportunity/2016/04/aas-task-force-education-begins-its-work)

A Tracer Method for Computing Type Ia Supernova Yields: Burning Model Calibration, Reconstruction of Thickened Flames, and Verification for Planar Detonations

Authors: Dean M. Townsley, Broxton J. Miles, F. X. Timmes, Alan C. Calder, Edward F. Brown

Abstract: We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing… ▽ More We refine our previously introduced parameterized model for explosive carbon-oxygen fusion during thermonuclear supernovae (SN Ia) by adding corrections to post-processing of recorded Lagrangian fluid element histories to obtain more accurate isotopic yields. Deflagration and detonation products are verified for propagation in a uniform density medium. A new method is introduced for reconstructing the temperature-density history within the artificially thick model deflagration front. We obtain better than 5\% consistency between the electron capture computed by the burning model and yields from post-processing. For detonations, we compare to a benchmark calculation of the structure of driven steady-state planar detonations performed with a large nuclear reaction network and error-controlled integration. We verify that, for steady-state planar detonations down to a density of 5x10^6 g/cc, our post processing matches the major abundances in the benchmark solution typically to better than 10% for times greater than 0.01 s after the shock front passage. As a test case to demonstrate the method, presented here with post-processing for the first time, we perform a two dimensional simulation of a SN Ia in the Chandrasekhar-mass deflagration-detonation transition (DDT) scenario. We find that reconstruction of deflagration tracks leads to slightly more complete silicon burning than without reconstruction. The resulting abundance structure of the ejecta is consistent with inferences from spectroscopic studies of observed SNe Ia. We confirm the absence of a central region of stable Fe-group material for the multi-dimensional DDT scenario. Detailed isotopic yields are tabulated and only change modestly when using deflagration reconstruction. △ Less

Submitted 16 May, 2016; originally announced May 2016.

Comments: 28 pages, 16 figures, Accepted to the Astrophysical Journal Supplemental

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

Authors: Daniel Lecoanet, Josiah Schwab, Eliot Quataert, Lars Bildsten, F. X. Timmes, Keaton J. Burns, Geoffrey M. Vasil, Jeffrey S. Oishi, Benjamin P. Brown

Abstract: It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much l… ▽ More It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution. △ Less

Submitted 25 September, 2016; v1 submitted 29 March, 2016; originally announced March 2016.

Comments: Accepted to ApJ

Properties of Carbon-Oxygen White Dwarfs From Monte Carlo Stellar Models

Authors: C. E. Fields, R. Farmer, I. Petermann, C. Iliadis, F. X. Timmes

Abstract: 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 M… ▽ More 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 M$_{\odot}$ 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_{\rm 1TP}$ $\approx$ 0.019 M$_{\odot}$ for the core mass at the first thermal pulse, $Δ$$t_{\rm{1TP}}$ $\approx$ 12.50 Myr for the age, $Δ\log(T_{\rm c}/{\rm K}) \approx$ 0.013 for the central temperature, $Δ\log(ρ_{\rm c}/{\rm g \ cm}^{-3}) \approx$ 0.060 for the central density, $ΔY_{\rm{e,c}} \approx$ 2.6$\times$10$^{-5}$ for the central electron fraction, $ΔX_{\rm c}(^{22}\rm{Ne}) \approx$ 5.8$\times$10$^{-4}$, $ΔX_{\rm c}(^{12}\rm{C}) \approx$ 0.392, and $ΔX_{\rm c}(^{16}\rm{O}) \approx$ 0.392. Uncertainties in the experimental $^{12}$C($α,γ)^{16}\rm{O}$, triple-$α$, and $^{14}$N($p,γ)^{15}\rm{O}$ reaction rates dominate these variations. We also consider a grid of 1 to 6 M$_{\odot}$ 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. △ Less

Submitted 21 March, 2016; originally announced March 2016.

Comments: Accepted for publication in The Astrophysical Journal; 19 Pages, 23 Figures, 5 Tables

White Paper on Nuclear Astrophysics

Authors: Almudena Arcones, Dan W. Bardayan, Timothy C. Beers, Lee A. Berstein, Jeffrey C. Blackmon, Bronson Messer, B. Alex Brown, Edward F. Brown, Carl R. Brune, Art E. Champagne, Alessandro Chieffi, Aaron J. Couture, Pawel Danielewicz, Roland Diehl, Mounib El-Eid, Jutta Escher, Brian D. Fields, Carla Fröhlich, Falk Herwig, William Raphael Hix, Christian Iliadis, William G. Lynch, Gail C. McLaughlin, Bradley S. Meyer, Anthony Mezzacappa , et al. (18 additional authors not shown)

Abstract: This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation… ▽ More This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade. △ Less

Submitted 24 March, 2016; v1 submitted 4 March, 2016; originally announced March 2016.

Comments: 131 pages, community white paper based on 2012 JINA Town Meeting in Detroit, MI, and 2014 APS Town Meeting in College Station, TX, corrected author names

Convective Properties of Rotating Two-Dimensional Core-Collapse Supernova Progenitors

Authors: Emmanouil Chatzopoulos, Sean M. Couch, W. David Arnett, F. X. Timmes

Abstract: We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow… ▽ More We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions. △ Less

Submitted 25 January, 2016; v1 submitted 21 January, 2016; originally announced January 2016.

Comments: 10 pages, 9 figures

Modeling Astrophysical Explosions with Sustained Exascale Computing

Authors: M. Zingale, A. C. Calder, C. M. Malone, F. X. Timmes

Abstract: Our understanding of stars and their fates is based on coupling observations to theoretical models. Unlike laboratory physicists, we cannot perform experiments on stars, but rather must patiently take what nature allows us to observe. Simulation offers a means of virtual experimentation, enabling a detailed understanding of the most violent ongoing explosions in the Universe---the deaths of stars. Our understanding of stars and their fates is based on coupling observations to theoretical models. Unlike laboratory physicists, we cannot perform experiments on stars, but rather must patiently take what nature allows us to observe. Simulation offers a means of virtual experimentation, enabling a detailed understanding of the most violent ongoing explosions in the Universe---the deaths of stars. △ Less

Submitted 16 October, 2015; originally announced October 2015.

Comments: Response to RFI NOT-GM-15-122: Science Drivers Requiring Capable Exascale High Performance Computing

On Measuring the Metallicity of a Type Ia Supernova's Progenitor

Authors: Broxton J. Miles, Daniel R. van Rossum, Dean M. Townsley, F. X. Timmes, Aaron P. Jackson, Alan C. Calder, Edward F. Brown

Abstract: In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single de… ▽ More In Type Ia Supernovae (\sneia), the relative abundances of chemical elements are affected by the neutron excess in the composition of the progenitor white dwarf. Since these products leave signatures in the spectra near maximum light, spectral features may be used to constrain the composition of the progenitor. We calculate the nucleosynthetic yields for three \snia simulations, assuming single degenerate, Chandrasekhar mass progenitors, for a wide range of progenitor metallicities, and calculate synthetic light curves and spectra to explore correlations between progenitor metallicity and the strength of spectral features. We use two 2D simulations of the deflagration-detonation-transition scenario with different $^{56}$Ni yields and the W7 simulation to control for differences between explosion models and total yields. While the overall yields of intermediate mass elements (16 $<$ A $\leq$ 40) differ between the three cases, trends in the yields are similar. With increasing metallicity, $^{28}$Si yields remain nearly constant, $^{40}$Ca yields decline, and Ti and $^{54}$Fe yields increase. In the synthetic spectra, we identify two features at 30 days post explosion that appear to deepen with progenitor metallicity: a Ti feature around 4200\,Å and a Fe feature around 5200\,Å\@. In all three simulations, their pseudo equivalent widths show a systematic trend with progenitor metallicity. This suggests that these two features may allow differentiation among progenitor metallicities of observed \sneia and potentially help reduce the intrinsic Hubble scatter. △ Less

Submitted 3 May, 2016; v1 submitted 24 August, 2015; originally announced August 2015.

Comments: 20 pages, 18 figures, accepted by the Astrophysical Journal

Carbon-rich presolar grains from massive stars. Subsolar 12C/13C and 14N/15N ratios and the mystery of 15N

Authors: M. Pignatari, E. Zinner, P. Hoppe, C. J. Jordan, B. K. Gibson, R. Trappitsch, F. Herwig, C. Fryer, R. Hirschi, F. X. Timmes

Abstract: Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 1… ▽ More Carbon-rich grains with isotopic anomalies compared to the Sun are found in primitive meteorites. They were made by stars, and carry the original stellar nucleosynthesis signature. Silicon carbide grains of Type X and C, and low-density graphites condensed in the ejecta of core-collapse supernovae. We present a new set of models for the explosive He shell and compare them with the grains showing 12C/13C and 14N/15N ratios lower than solar. In the stellar progenitor H was ingested into the He shell and not fully destroyed before the explosion. Different explosion energies and H concentrations are considered. If the SN shock hits the He-shell region with some H still present, the models can reproduce the C and N isotopic signatures in C-rich grains. Hot-CNO cycle isotopic signatures are obtained, including a large production of 13C and 15N. The short-lived radionuclides 22Na and 26Al are increased by orders of magnitude. The production of radiogenic 22Ne from the decay of 22Na in the He shell might solve the puzzle of the Ne-E(L) component in low-density graphite grains. This scenario is attractive for the SiC grains of type AB with 14N/15N ratios lower than solar, and provides an alternative solution for SiC grains originally classified as nova grains. Finally, this process may contribute to the production of 14N and 15N in the Galaxy, helping to produce the 14N/15N ratio in the solar system. △ Less

Submitted 30 June, 2015; originally announced June 2015.

Comments: 13 pages, 4 figures, 1 table, The Astrophysical Journal Letters accepted

Modules for Experiments in Stellar Astrophysics (MESA): Binaries, Pulsations, and Explosions

Authors: Bill Paxton, Pablo Marchant, Josiah Schwab, Evan B. Bauer, Lars Bildsten, Matteo Cantiello, Luc Dessart, R. Farmer, H. Hu, N. Langer, R. H. D. Townsend, Dean M. Townsley, F. X. Timmes

Abstract: We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of… ▽ More We substantially update the capabilities of the open-source software instrument Modules for Experiments in Stellar Astrophysics (MESA). MESA can now simultaneously evolve an interacting pair of differentially rotating stars undergoing transfer and loss of mass and angular momentum, greatly enhancing the prior ability to model binary evolution. New MESA capabilities in fully coupled calculation of nuclear networks with hundreds of isotopes now allow MESA to accurately simulate advanced burning stages needed to construct supernova progenitor models. Implicit hydrodynamics with shocks can now be treated with MESA, enabling modeling of the entire massive star lifecycle, from pre-main sequence evolution to the onset of core collapse and nucleosynthesis from the resulting explosion. Coupling of the GYRE non-adiabatic pulsation instrument with MESA allows for new explorations of the instability strips for massive stars while also accelerating the astrophysical use of asteroseismology data. We improve treatment of mass accretion, giving more accurate and robust near-surface profiles. A new MESA capability to calculate weak reaction rates "on-the-fly" from input nuclear data allows better simulation of accretion induced collapse of massive white dwarfs and the fate of some massive stars. We discuss the ongoing challenge of chemical diffusion in the strongly coupled plasma regime, and exhibit improvements in MESA that now allow for the simulation of radiative levitation of heavy elements in hot stars. We close by noting that the MESA software infrastructure provides bit-for-bit consistency for all results across all the supported platforms, a profound enabling capability for accelerating MESA's development. △ Less

Submitted 8 January, 2017; v1 submitted 9 June, 2015; originally announced June 2015.

Comments: Published in the ApJ Supplement Series; updated to account for published erratum

Journal ref: ApJS, 220, 15 (2015)

On Carbon Burning in Super Asymptotic Giant Branch Stars

Authors: R. Farmer, C. E. Fields, F. X. Timmes

Abstract: We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters su… ▽ More We explore the detailed and broad properties of carbon burning in Super Asymptotic Giant Branch (SAGB) stars with 2755 MESA stellar evolution models. The location of first carbon ignition, quenching location of the carbon burning flames and flashes, angular frequency of the carbon core, and carbon core mass are studied as a function of the ZAMS mass, initial rotation rate, and mixing parameters such as convective overshoot, semiconvection, thermohaline and angular momentum transport. In general terms, we find these properties of carbon burning in SAGB models are not a strong function of the initial rotation profile, but are a sensitive function of the overshoot parameter. We quasi-analytically derive an approximate ignition density, $ρ_{ign} \approx 2.1 \times 10^6$ g cm$^{-3}$, to predict the location of first carbon ignition in models that ignite carbon off-center. We also find that overshoot moves the ZAMS mass boundaries where off-center carbon ignition occurs at a nearly uniform rate of $ΔM_{\rm ZAMS}$/$Δf_{\rm{ov}}\approx$ 1.6 $M_{\odot}$. For zero overshoot, $f_{\rm{ov}}$=0.0, our models in the ZAMS mass range $\approx$ 8.9 to 11 $M_{\odot}$ show off-center carbon ignition. For canonical amounts of overshooting, $f_{\rm{ov}}$=0.016, the off-center carbon ignition range shifts to $\approx$ 7.2 to 8.8 $M_{\odot}$. Only systems with $f_{\rm{ov}}$ $\geq 0.01$ and ZAMS mass $\approx$ 7.2-8.0 $M_{\odot}$ show carbon burning is quenched a significant distance from the center. These results suggest a careful assessment of overshoot modeling approximations on claims that carbon burning quenches an appreciable distance from the center of the carbon core. △ Less

Submitted 2 June, 2015; originally announced June 2015.

Comments: Accepted ApJ; 23 pages, 21 figures, 5 tables

Constraints on explosive silicon burning in core-collapse supernovae from measured Ni/Fe ratios

Authors: A. Jerkstrand, F. X. Timmes, G. Magkotsios, S. A. Sim, C. Fransson, J. Spyromilio, A. Heger, B. Müller, J. Sollerman, S. J. Smartt

Abstract: Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54F… ▽ More Measurements of explosive nucleosynthesis yields in core-collapse supernovae provide tests for explosion models. We investigate constraints on explosive conditions derivable from measured amounts of nickel and iron after radioactive decays using nucleosynthesis networks with parameterized thermodynamic trajectories. The Ni/Fe ratio is for most regimes dominated by the production ratio of 58Ni/(54Fe + 56Ni), which tends to grow with higher neutron excess and with higher entropy. For SN 2012ec, a supernova that produced a Ni/Fe ratio of $3.4\pm1.2$ times solar, we find that burning of a fuel with neutron excess $η\approx 6\times 10^{-3}$ is required. Unless the progenitor metallicity is over 5 times solar, the only layer in the progenitor with such a neutron excess is the silicon shell. Supernovae producing large amounts of stable nickel thus suggest that this deep-lying layer can be, at least partially, ejected in the explosion. We find that common spherically symmetric models of $M_{\rm ZAMS} \lesssim 13$ Msun stars exploding with a delay time of less than one second ($M_{\rm cut} < 1.5$ Msun) are able to achieve such silicon-shell ejection. Supernovae that produce solar or sub-solar Ni/Fe ratios, such as SN 1987A, must instead have burnt and ejected only oxygen-shell material, which allows a lower limit to the mass cut to be set. Finally, we find that the extreme Ni/Fe value of 60-75 times solar derived for the Crab cannot be reproduced by any realistic-entropy burning outside the iron core, and neutrino-neutronization obtained in electron-capture models remains the only viable explanation. △ Less

Submitted 20 May, 2015; originally announced May 2015.

Comments: 13 pages, 9 figures, accepted for publication in ApJ