Showing 1–50 of 90 results for author: Hix, W R

Hydrodynamic Simulations of Oxygen-Neon Classical Novae as Galactic $^7$Li Producers and Potential Accretion Induced Collapse Progenitors

Authors: Sumner Starrfield, Maitrayee Bose, Christian Iliadis, W. Raphael Hix, Charles E. Woodward, R. Mark Wagner

Abstract: We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on oxygen-neon (ONe) white dwarfs (WDs). Using NOVA, a one-dimensional hydrodynamic computer code, we accrete Solar matter until the TNR is ongoing and then switch to a mixed composition. This approach is guided by the results of multi-dimensional studies of TNRs in WDs which find… ▽ More We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on oxygen-neon (ONe) white dwarfs (WDs). Using NOVA, a one-dimensional hydrodynamic computer code, we accrete Solar matter until the TNR is ongoing and then switch to a mixed composition. This approach is guided by the results of multi-dimensional studies of TNRs in WDs which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment of the CNONeMg elements are reached that agree with observations of CN ejecta abundances. Because the amount of accreted material is inversely proportional to the oxygen abundance, by first accreting Solar matter, the amount of accreted material is larger than in those simulations with an initially enriched composition. We vary the mass of the WD (from 0.6 Msun to 1.35 Msun) and the composition of the mixed materials. Our results show large enrichments of 7Be in the ejected gases implying that ONe CNe and CO CNe (Starrfield et al. 2020) may be responsible for a significant fraction (about 100 Msun) of the galactic 7Li ( about 1000 Msun). The production of 22Na and 26Al in CN explosions and the gamma-ray emission predicted by our simulations is discussed. The WDs in all our simulations eject less material than they accrete and we predict that the WD is growing in mass as a consequence of the CN outburst. ONe CNe, therefore, may be an important channel for accretion induced collapse (AIC) events. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 50 pages, 15 figures, 9 Tables accepted to the Astrophysical Journal This article is a follow on to Starrfield et al. ApJ, 895, 70, 2020 [arXiv:1910.00575]. In that paper we presented simulations for CO novae in this paper we consider ONe novae. The outline is the same but in this case we consider evolution to accretion induced collapse instead of SN Ia explosions

Impact of Pycnonuclear Fusion Uncertainties on the Cooling of Accreting Neutron Star Crusts

Authors: R. Jain, E. F. Brown, H. Schatz, A. V. Afanasjev, M. Beard, L. R. Gasques, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Lau, P. Moller, W. J. Ong, M. Wiescher, Y. Xu

Abstract: The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncerta… ▽ More The observation of X-rays during quiescence from transiently accreting neutron stars provides unique clues about the nature of dense matter. This, however, requires extensive modeling of the crusts and matching the results to observations. The pycnonuclear fusion reaction rates implemented in these models are theoretically calculated by extending phenomenological expressions and have large uncertainties spanning many orders of magnitude. We present the first sensitivity studies of these pycnonuclear fusion reactions in realistic network calculations. We also couple the reaction network with the thermal evolution code dStar to further study their impact on the neutron star cooling curves in quiescence. Varying the pycnonuclear fusion reaction rates alters the depth at which nuclear heat is deposited although the total heating remains constant. The enhancement of the pycnonuclear fusion reaction rates leads to an overall shallower deposition of nuclear heat. The impurity factors are also altered depending on the type of ashes deposited on the crust. These total changes correspond to a variation of up to 9 eV in the modeled cooling curves. While this is not sufficient to explain the shallow heat source, it is comparable to the observational uncertainties and can still be important for modeling the neutron star crust. △ Less

Submitted 1 October, 2023; v1 submitted 20 June, 2023; originally announced June 2023.

Comments: AASTeX63, 11 pages with 9 figures

Journal ref: ApJ 955 51 (2023)

Tracer particles for core-collapse supernova nucleosynthesis: The advantages of moving backward

Authors: Andre Sieverding, Preston G. Waldrop, J. Austin Harris, W. Raphael Hix, Eric J. Lentz, Stephen W. Bruenn, O. E. Bronson Messer

Abstract: After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multi-dimensional simulations. Obtaining nucleosynthesis yields based on such multi-dimensional core-collapse supernova (CCSN) simulations, however, is not straightforward and frequently tracer particles are employed. Tracer particles may b… ▽ More After decades, the theoretical study of core-collapse supernova explosions is moving from parameterized, spherically symmetric models to increasingly realistic multi-dimensional simulations. Obtaining nucleosynthesis yields based on such multi-dimensional core-collapse supernova (CCSN) simulations, however, is not straightforward and frequently tracer particles are employed. Tracer particles may be tracked in situ during the simulation, but often they are reconstructed in a post-processing step based on the information saved during the hydrodynamics simulation. Reconstruction can be done in a number of ways and here we compare the approaches of backward and forward integration of the equations of motion to the results based on inline particle trajectories. We find that both methods agree reasonably well with the inline results for isotopes for which a large number of particles contribute. However, for rarer isotopes that are produced only by a small number of particle trajectories, deviations can be large. For our setup, we find that backward integration leads to a better agreement with the inline particles by more accurately reproducing the conditions following freeze-out from nuclear statistical equilibrium, because the establishment of nuclear statistical equilibrium erases the need for detailed trajectories at earlier times. Based on our results, if inline tracers are unavailable, we recommend backward reconstruction, to the point when nuclear statistical equilibrium last applied, with an interval between simulation snapshots of at most 1 ms for nucleosynthesis post-processing. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: 23 pages, 18 figures, sumitted to ApJ

Comparison of the Core-Collapse Evolution of Two Nearly Equal Mass Progenitors

Authors: Stephen W. Bruenn, Andre Sieverding, Eric J. Lentz, Tuguldur Sukhbold, W. Raphael Hix, Leah N. Huk, J. Austin Harris, O. E. Bronson Messer, Anthony Mezzacappa

Abstract: We compare the core-collapse evolution of a pair of 15.8 $M_\odot$ stars with significantly different internal structures, a consequence of bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 $M_\odot $ progenitors have core masses of 1.47 and 1.78 $M_\odot$ and compactness parameters $ξ_{1.75}$ of 0.302 and 0.604. The core collapse simulations… ▽ More We compare the core-collapse evolution of a pair of 15.8 $M_\odot$ stars with significantly different internal structures, a consequence of bimodal variability exhibited by massive stars during their late evolutionary stages. The 15.78 and 15.79 $M_\odot $ progenitors have core masses of 1.47 and 1.78 $M_\odot$ and compactness parameters $ξ_{1.75}$ of 0.302 and 0.604. The core collapse simulations are carried out in 2D to nearly 3 s post-bounce and show substantial differences in the times of shock revival and explosion energies. The 15.78 $M_\odot$ model explodes promptly at 120 ms post-bounce when a strong density decrement at the Si--Si/O shell interface encounters the stalled shock. The 15.79 $M_\odot$ model, which lacks the density decrement, takes 100 ms longer to explode but ultimately produces a more powerful explosion. Larger mass accretion rate of the 15.79 $M_\odot$ model during the first 0.8 s post-bounce results in larger $ν_{e}$/$\bar ν_{e}$ luminosities and rms energies. The $ν_{e}$/$\bar ν_{e}$ luminosities and rms energies arising from the inner core are also larger in the 15.79 $M_\odot$ model throughout due to the larger negative temperature gradient of this core due to greater adiabatic compression. Larger luminosities and rms energies in the 15.79 $M_\odot$ model and a flatter and higher density heating region, result in more energy deposition behind the shock and more ejected matter with higher enthalpy. We find the ejected $^{56}$Ni mass of the 15.79 $M_\odot$ model is more than double that of the 15.78 $M_\odot$ model. Most of the ejecta in both models is moderately proton-rich, though counterintuitively the highest electron fraction ($Y_e=0.61$) ejecta in either model is in the less energetic 15.78 $M_\odot$ model while the lowest electron fraction ($Y_e=0.45$) ejecta in either model is in the 15.79 $M_\odot$ model. △ Less

Submitted 10 February, 2023; v1 submitted 22 November, 2022; originally announced November 2022.

Comments: 24 pages; Accepted for publication in ApJ

Journal ref: ApJ, 2023, 497, 35

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

Core Collapse Supernova Gravitational Wave Emission for Progenitors of 9.6, 15, and 25 Solar Masses

Authors: Anthony Mezzacappa, Pedro Marronetti, Ryan E. Landfield, Eric J. Lentz, W. Raphael Hix, J. Austin Harris, Stephen W. Bruenn, John M. Blondin, O. E. Bronson Messer, Jordi Casanova, Luke L. Kronzer

Abstract: We present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25 Solar masses, are all initially non-rotating, and are of two metallicities: zero and Solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross… ▽ More We present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and 25 Solar masses, are all initially non-rotating, and are of two metallicities: zero and Solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: Submitted to Physical Review D

The Impact of Resolution on Double-Detonation Models for Type Ia Supernovae

Authors: Fernando Rivas, J. Austin Harris, W. Raphael Hix, O. E. Bronson Messer

Abstract: Thermonuclear supernovae are the result of the violent unbinding of a white dwarf, but the precise nature of the explosion mechanism(s) is a matter of active debate. To this end, several specific scenarios have been proposed to explain the observable traits of SNe Ia. A promising pathway is the double-detonation scenario, where a white dwarf accretes a shell of helium-rich material from a companio… ▽ More Thermonuclear supernovae are the result of the violent unbinding of a white dwarf, but the precise nature of the explosion mechanism(s) is a matter of active debate. To this end, several specific scenarios have been proposed to explain the observable traits of SNe Ia. A promising pathway is the double-detonation scenario, where a white dwarf accretes a shell of helium-rich material from a companion and a detonation in the resulting helium shell is the primary cause of the explosion. Through a set of two-dimensional grid-based simulations of this scenario we clearly distinguish three phases of evolution: external helium-rich detonation, core compressive heating, and a final core carbon burn. Though final disruption of the whole system is achieved at all resolutions, only models with minimum resolutions of 4~km and better exhibit all three phases. Particularly, core compression heating is only observed for higher resolutions, producing qualitatively different nucleosynthetic outcomes. We identify the effect of finer spatial resolution on the mixing of hot silicon at the interface between the detonating helium layer and the underlying C/O WD as a primary driver of these dynamic differences. △ Less

Submitted 18 August, 2022; v1 submitted 24 May, 2022; originally announced May 2022.

Comments: 10 figures, 12 pages

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

The role of low-lying resonances for $^{10}$Be(p,$α$) reaction rate and implications for the formation of the Solar System

Authors: A. Sieverding, J. S. Randhawa, D. Zetterberg, R. J. deBoer, T. Ahn, R. Mancino, G. Martinez-Pinedo, W. R. Hix

Abstract: Evidence for the presence of short-lived radioactive isotopes when the Solar System formed is preserved in meteorites, providing insights into the conditions at the birth of our Sun. A low-mass core-collapse supernova had been postulated as a candidate for the origin of $^{10}$Be, reinforcing the idea that a supernova triggered the formation of the Solar System. We present a detailed study of the… ▽ More Evidence for the presence of short-lived radioactive isotopes when the Solar System formed is preserved in meteorites, providing insights into the conditions at the birth of our Sun. A low-mass core-collapse supernova had been postulated as a candidate for the origin of $^{10}$Be, reinforcing the idea that a supernova triggered the formation of the Solar System. We present a detailed study of the production of $^{10}$Be by the $ν$ process in supernovae, which is very sensitive to the reaction rate of the major destruction channel, $^{10}$Be(p,$α$)$^{7}$Li. With data from recent nuclear experiments that show the presence of a resonant state in $^{11}$B at $\approx$~193 keV, we derive new values for the $^{10}$Be(p,$α$)$^{7}$Li reaction rate which are significantly higher than previous estimates. We show that, with the new $^{10}$Be(p,$α$)$^{7}$Li reaction rate, a low mass CCSN is unlikely to produce enough $^{10}$Be to explain the observed $^{10}$Be/ $^{9}$Be ratio in meteorites, even for a wide range of neutrino spectra considered in our models. These findings point towards spallation reactions induced by solar energetic particles in the early Solar System as the origin of $^{10}$Be. △ Less

Submitted 13 June, 2022; v1 submitted 12 March, 2022; originally announced March 2022.

Comments: 12 pages, 9 Figures, 4 Tables. Submitted to Physical Review C

The Impact of Neutron Transfer Reactions on Heating and Cooling of Accreted Neutron Star Crusts

Authors: H. Schatz, Z. Meisel, E. F. Brown, S. S. Gupta, G. W. Hitt, W. R. Hix, R. Jain, R. Lau, P. Möller, W. -J. Ong, P. S. Shternin, Y. Xu, M. Wiescher

Abstract: Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now incl… ▽ More Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and of long-term cooling in transiently accreting systems. It was recently suggested that previously neglected neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions, and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust. △ Less

Submitted 21 December, 2021; originally announced December 2021.

Comments: 8 pages, 10 figures, accepted for Astrophysical Journal

Three Dimensional Core-Collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

Authors: Michael A. Sandoval, W. Raphael Hix, O. E. Bronson Messer, Eric J. Lentz, J. Austin Harris

Abstract: We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino-radiation hydrodynamic simulations of the neutrino-driven phase performed with the CHIMERA code. We consider a 9.6-$M_{\odot}$ zero-metallicity progenitor starting from both 2D and 3D CHIMERA models, and a 10-… ▽ More We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino-radiation hydrodynamic simulations of the neutrino-driven phase performed with the CHIMERA code. We consider a 9.6-$M_{\odot}$ zero-metallicity progenitor starting from both 2D and 3D CHIMERA models, and a 10-$M_{\odot}$ solar-metallicity progenitor starting from a 2D CHIMERA model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh-Taylor (R-T) "bullets" determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum $^{56}\rm{Ni}$ velocities of ${\sim} 1950~\rm{km~s}^{-1}$ and ${\sim} 1750~\rm{km~s}^{-1}$ at shock breakout from 2D and 3D 9.6-$M_{\odot}$ CHIMERA models, respectively, due to the bullets' ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher velocity structures is suppressed when the 2D CHIMERA model and 3D FLASH model meshes are aligned. The development of faster growing spherical-bubble structures, as opposed to the slower growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10-$M_{\odot}$ model, which achieves maximum $^{56}\rm{Ni}$ velocities of ${\sim} 2500~\rm{km~s}^{-1}$ at shock breakout. △ Less

Submitted 12 November, 2021; v1 submitted 2 June, 2021; originally announced June 2021.

Comments: 31 pages, 23 figures, published in ApJ

Journal ref: ApJ 921 113 (2021)

Evaluation of Experimental Constraints on the $^{44}$Ti($α$,p)$^{47}$V Reaction Cross Section Relevant for Supernovae

Authors: Kelly A. Chipps, Philip Adsley, Manoel Couder, W. Raphael Hix, Zachary Meisel, Konrad Schmidt

Abstract: Due to its importance as an astronomical observable in core-collapse supernovae (CCSNe), the reactions producing and destroying $^{44}$Ti must be well constrained. Generally, statistical model calculations such as Hauser-Feshbach are employed when experimental cross sections are not available, but the variation in such adopted rates can be large. Here, data from the literature is compared with sta… ▽ More Due to its importance as an astronomical observable in core-collapse supernovae (CCSNe), the reactions producing and destroying $^{44}$Ti must be well constrained. Generally, statistical model calculations such as Hauser-Feshbach are employed when experimental cross sections are not available, but the variation in such adopted rates can be large. Here, data from the literature is compared with statistical model calculations of the $^{44}$Ti($α$,p)$^{47}$V reaction cross section and used to constrain the possible reaction rate variation over the temperatures relevant to CCSNe. Suggestions for targeted future measurements are given. △ Less

Submitted 30 November, 2020; v1 submitted 31 July, 2020; originally announced August 2020.

Comments: 7 pages, 4 figures. PRC, published

Journal ref: Phys. Rev. C 102, 035806 (2020)

Gravitational-wave signal of a core-collapse supernova explosion of a 15 Solar mass star

Authors: Anthony Mezzacappa, Pedro Marronetti, Ryan E. Landfield, Eric J. Lentz, Konstantin N. Yakunin, Stephen W. Bruenn, W. Raphael Hix, O. E. Bronson Messer, Eirik Endeve, John M. Blondin, J. Austin Harris

Abstract: We report on the gravitational wave signal computed in the context of a three-dimensional simulation of a core collapse supernova explosion of a 15 Solar mass star. The simulation was performed with our neutrino hydrodynamics code Chimera. We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding sp… ▽ More We report on the gravitational wave signal computed in the context of a three-dimensional simulation of a core collapse supernova explosion of a 15 Solar mass star. The simulation was performed with our neutrino hydrodynamics code Chimera. We detail the gravitational wave strains as a function of time, for both polarizations, and discuss their physical origins. We also present the corresponding spectral signatures. Gravitational wave emission in our model has two key features: low-frequency emission (< 200 Hz) emanates from the gain layer as a result of neutrino-driven convection and the SASI and high-frequency emission (> 600 Hz) emanates from the proto-neutron star due to Ledoux convection within it. The high-frequency emission dominates the gravitational wave emission in our model and emanates largely from the convective layer itself, not from the convectively stable layer above it, due to convective overshoot. Moreover, the low-frequency emission emanates from the gain layer itself, not from the proto-neutron star, due to accretion onto it. We provide evidence of the SASI in our model and demonstrate that the peak of our low-frequency gravitational wave emission spectrum corresponds to it. Given its origin in the gain layer, we classify the SASI emission in our model as p-mode emission and assign a purely acoustic origin, not a vortical-acoustic origin, to it. Our dominant proto-neutron star gravitational wave emission is not well characterized by emission from surface g-modes, complicating the relationship between peak frequencies observed and the mass and radius of the proto-neutron star expressed by analytic estimates under the assumption of surface g-mode emission. We present our frequency normalized characteristic strain along with the sensitivity curves of current- and next-generation gravitational wave detectors. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Journal ref: Physical Review D 102, 023027 (2020)

Hydrodynamic Simulations of Classical Novae; CO and ONe White Dwarfs are Supernova Ia Progenitors

Authors: S. Starrfield, M. Bose, C. Iliadis, W. R. Hix, C. E. Woodward, R. M. Wagner

Abstract: Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of th… ▽ More Cataclysmic Variables (CVs) and Symbiotic Binaries are close (or not so close) binary star systems which contain both a white dwarf (WD) primary and a larger cooler secondary star that typically fills its Roche Lobe. The cooler star is losing mass through the inner Lagrangian point of the binary and a fraction of this material is accreted by the WD. Here we report on our hydrodynamic studies of the thermonuclear runaway (TNR) in the accreted material that ends in a Classical Nova explosion. We have followed the evolution of the TNRs on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs. We report on 3 studies in this paper. First, simulations in which we accrete only solar matter using NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar matter and compare the results. Third, we accrete solar matter until the TNR is ongoing and then switch the composition in the accreted layers to a mixed composition: either 25% WD and 75% solar or 50% WD and 50% Solar matter. We find that the amount of accreted material is inversely proportional to the initial 12C abundance (as expected). Thus, accreting solar matter results in a larger amount of accreted material to fuel the outburst; much larger than in earlier studies where a mixed composition was assumed from the beginning of the simulation. Our most important result is that all these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit. △ Less

Submitted 2 June, 2020; originally announced June 2020.

Comments: To appear in "The Golden Age of Cataclysmic Variables and Related Objects V (GOLDEN2019)" 2-7 September 2019, Palermo, Italy. Proceedings of Science, in press

On the character of turbulence in self-consistent models of core-collapse supernovae

Authors: Jordi Casanova, Eirik Endeve, Eric J. Lentz, O. E. Bronson Messer, W. Raphael Hix, J. Austin Harris, Stephen W. Bruenn

Abstract: Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities and the development of turbulent convection will determine the morphology of… ▽ More Neutrino-driven convection plays a crucial role in the development of core-collapse supernova (CCSN) explosions. However, the complex mechanism that triggers the shock revival and the subsequent explosion has remained inscrutable for many decades. Multidimensional simulations suggest that the growth of fluid instabilities and the development of turbulent convection will determine the morphology of the explosion. We have performed 3D simulations using spherical-polar coordinates covering a reduced angular extent (90 degree computational domain), and with angular resolutions of 2 degrees, 1 degree, 1/2 degree, and 1/4 degree, to study the development of turbulence in core-collapse supernova explosions on a time scale of order 100 ms. We have employed the multi-physics Chimera code that includes detailed nuclear physics and spectral neutrino transport. Coarse resolution models do not develop an inertial range, presumably due to the bottleneck effect, such that the energy is prevented from cascading down to small scales and tends to accumulate at large scales. High-resolution models instead, start to recover the k^{-5/3} scaling of Kolmogorov's theory. Stochasticity and few simulation samples limit our ability to predict the development of explosions. Over the simulated time period, our models show no clear trend in improving (or diminishing) conditions for explosion as the angular resolution is increased. However, we find that turbulence provides an effective pressure behind the shock (approx. 40 - 50 % of the thermal pressure), which can contribute to the shock revival and be conducive for the development of the explosion. Finally, we show that the turbulent energy power spectrum of reduced angular extent and full 4 pi models are consistent, thus indicating that a 90 degree computational domain is an adequate configuration to study the character of turbulence in CCSNe. △ Less

Submitted 4 April, 2020; originally announced April 2020.

Comments: Accepted Special Issue Article in Physica Scripta

Carbon-Oxygen Classical Novae are Galactic $^7$Li Producers as well as Potential Supernova Ia Progenitors

Authors: Sumner Starrfield, Maitrayee Bose, Christian Iliadis, W. Raphael Hix, Charles E. Woodward, R. Mark Wagner

Abstract: We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on Carbon Oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M$_\odot$ to 1.35 M$_\odot$) and the composition of the accreted material. Our simulations are guided by the results of multi-dimensional studies of TNRs in WDs that find sufficient mixing with WD c… ▽ More We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on Carbon Oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M$_\odot$ to 1.35 M$_\odot$) and the composition of the accreted material. Our simulations are guided by the results of multi-dimensional studies of TNRs in WDs that find sufficient mixing with WD core material occurs after the TNR is well underway, reaching levels of enrichment that agree with observations of CN ejecta abundances. We use NOVA (our 1-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial $^{12}$C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of $^7$Be in the ejected gases implying that CO CNe may be responsible for a significant fraction ($\sim$ 100 M$_\odot$) of the $^7$Li in the galaxy ($\sim$1000 M$_\odot$). In addition, although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion-outburst-accretion cycle and CO CNe may be an important channel of Supernova Ia progenitors. △ Less

Submitted 1 October, 2019; originally announced October 2019.

Comments: Submitted to the Astrophysical Journal on 24 April 2019 (30 pages, 12 Figures, 6 tables)

Journal ref: ApJ 895 70 (2020)

Chimera: A massively parallel code for core-collapse supernova simulation

Authors: Stephen W. Bruenn, John M. Blondin, W. Raphael Hix, Eric J. Lentz, O. E. Bronson Messer, Anthony Mezzacappa, Eirik Endeve, J. Austin Harris, Pedro Marronetti, Reuben D. Budiardja, Merek A. Chertkow, Ching-Tsai Lee

Abstract: We provide a detailed description of the Chimera code, a code developed to model core collapse supernovae in multiple spatial dimensions. The core collapse supernova explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core flu… ▽ More We provide a detailed description of the Chimera code, a code developed to model core collapse supernovae in multiple spatial dimensions. The core collapse supernova explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the foundational physics of the neutrino-stellar core weak interactions and the equations of state of all stellar core constituents -particularly, the nuclear equation of state associated with nucleons, both free and bound in nuclei. Chimera, by incorporating detailed neutrino transport, realistic neutrino-matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with core collapse supernovae. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. Chimera is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications. △ Less

Submitted 25 February, 2020; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: Accepted for publication in ApJS

Journal ref: ApJS (2020), 248, 11

Nuclear Reactions in the Crusts of Accreting Neutron Stars

Authors: R. Lau, M. Beard, S. S. Gupta, H. Schatz, A. V. Afanasjev, E. F. Brown, A. Deibel, L. R. Gasques, G. W. Hitt, W. R. Hix, L. Keek, P. Möller, P. S. Shternin, A. Steiner, M. Wiescher, Y. Xu

Abstract: X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nucl… ▽ More X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models. △ Less

Submitted 21 April, 2018; v1 submitted 10 March, 2018; originally announced March 2018.

Comments: 25 Pages, accepted for publication in Ap. J

Gravitational wave signals from multi-dimensional core-collapse supernova explosion simulations

Authors: K. N. Yakunin, E. Endeve, A. Mezzacappa, M. Szczpanczyk, M. Zanolin, P. Marronetti, E. J. Lentz, S. W. Bruenn, W. R. Hix, O. E. B. Messer, J. M. Blondin, J. A. Harris

Abstract: In this work we report briefly on the gravitational wave (GW) signal computed in the context of a self-consistent, 3D simulation of a core-collapse supernova (CCSN) explosion of a 15M$_\odot$ progenitor star. We present a short overview of the GW signal, including signal amplitude, frequency distribution, and the energy emitted in the form of GWs for each phase of explosion, along with neutrino lu… ▽ More In this work we report briefly on the gravitational wave (GW) signal computed in the context of a self-consistent, 3D simulation of a core-collapse supernova (CCSN) explosion of a 15M$_\odot$ progenitor star. We present a short overview of the GW signal, including signal amplitude, frequency distribution, and the energy emitted in the form of GWs for each phase of explosion, along with neutrino luminosities, and discuss correlations between them. △ Less

Submitted 23 October, 2017; originally announced October 2017.

Comments: 5 pages, 2 figures, submitted in the proceedings of the 52nd Rencontres de Moriond

Implications for Post-Processing Nucleosynthesis of Core-Collapse Supernova Models with Lagrangian Particles

Authors: J. Austin Harris, W. Raphael Hix, Merek A. Chertkow, C. -T. Lee, Eric J. Lentz, O. E. Bronson Messer

Abstract: We investigate core-collapse supernova (CCSN) nucleosynthesis with self-consistent, axisymmetric (2D) simulations performed using the radiation-hydrodynamics code Chimera. Computational costs have traditionally constrained the evolution of the nuclear composition within multidimensional CCSN models to, at best, a 14-species $α$-network capable of tracking only $(α,γ)$ reactions from $^{4}$He to… ▽ More We investigate core-collapse supernova (CCSN) nucleosynthesis with self-consistent, axisymmetric (2D) simulations performed using the radiation-hydrodynamics code Chimera. Computational costs have traditionally constrained the evolution of the nuclear composition within multidimensional CCSN models to, at best, a 14-species $α$-network capable of tracking only $(α,γ)$ reactions from $^{4}$He to $^{60}$Zn. Such a simplified network limits the ability to accurately evolve detailed composition and neutronization or calculate the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. However, limitations such as poor spatial resolution of the tracer particles, inconsistent thermodynamic evolution, including misestimation of expansion timescales, and uncertain determination of the multidimensional mass-cut at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of such uncertainties for four self-consistent axisymmetric CCSN models initiated from stellar metallicity, non-rotating progenitors of 12 $M_\odot$, 15 $M_\odot$, 20 $M_\odot$, and 25 $M_\odot$ and evolved with the smaller $α$-network to more than 1 s after the launch of an explosion. △ Less

Submitted 7 June, 2017; v1 submitted 30 January, 2017; originally announced January 2017.

Comments: 25 pages, 17 figures, accepted for publication in ApJ

Journal ref: ApJ, 843, 21 (2017)

The Gravitational Wave Signal of a Core Collapse Supernova Explosion of a 15M$_\odot$ Star

Authors: Konstantin N. Yakunin, Anthony Mezzacappa, Pedro Marronetti, Eric J. Lentz, Stephen W. Bruenn, W. Raphael Hix, O. E. Bronson Messer, Eirik Endeve, John M. Blondin, J. Austin Harris

Abstract: In this Letter, we report on the gravitational wave signal computed in the context of an $ab$ $initio$, three-dimensional simulation of a core collapse supernova explosion, beginning with a 15M$_\odot$ star and using state-of-the-art weak interactions. The simulation was performed with our neutrino hydrodynamics code Chimera. We discuss the potential for detection of our predicted gravitational si… ▽ More In this Letter, we report on the gravitational wave signal computed in the context of an $ab$ $initio$, three-dimensional simulation of a core collapse supernova explosion, beginning with a 15M$_\odot$ star and using state-of-the-art weak interactions. The simulation was performed with our neutrino hydrodynamics code Chimera. We discuss the potential for detection of our predicted gravitational signal by the current generation of gravitational wave detectors. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: 6 pages, 4 figures

The Thermonuclear Runaway and the Classical Nova Outburst

Authors: S. Starrfield, C. Iliadis, W. R. Hix

Abstract: Nova explosions occur on the white dwarf component of a Cataclysmic Variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a Classical Nova explosion. We describe both the recent advances in our understanding of the progress of the outburs… ▽ More Nova explosions occur on the white dwarf component of a Cataclysmic Variable binary stellar system that is accreting matter lost by its companion. When sufficient material has been accreted by the white dwarf, a thermonuclear runaway occurs and ejects material in what is observed as a Classical Nova explosion. We describe both the recent advances in our understanding of the progress of the outburst and outline some of the puzzles that are still outstanding. We report on the effects of improving both the nuclear reaction rate library and including a modern nuclear reaction network in our one-dimensional, fully implicit, hydrodynamic computer code. In addition, there has been progress in observational studies of Supernovae Ia with implications about the progenitors and we discuss that in this review. △ Less

Submitted 13 May, 2016; originally announced May 2016.

Comments: 18 pages, 7 figures

Journal ref: Publications of the Astronomical Society of the Pacific, 128:051001, 2016

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

The Multi-Dimensional Character of Core-Collapse Supernovae

Authors: W. R. Hix, E. J. Lentz, S. W. Bruenn, A. Mezzacappa, O. E. B. Messer, E. Endeve, J. M. Blondin, J. A. Harris, P. Marronett, K. N. Yakunin

Abstract: Core-collapse supernovae, the culmination of massive stellar evolution, are spectacular astronomical events and the principle actors in the story of our elemental origins. Our understanding of these events, while still incomplete, centers around a neutrino-driven central engine that is highly hydrodynamically unstable. Increasingly sophisticated simulations reveal a shock that stalls for hundreds… ▽ More Core-collapse supernovae, the culmination of massive stellar evolution, are spectacular astronomical events and the principle actors in the story of our elemental origins. Our understanding of these events, while still incomplete, centers around a neutrino-driven central engine that is highly hydrodynamically unstable. Increasingly sophisticated simulations reveal a shock that stalls for hundreds of milliseconds before reviving. Though brought back to life by neutrino heating, the development of the supernova explosion is inextricably linked to multi-dimensional fluid flows. In this paper, the outcomes of three-dimensional simulations that include sophisticated nuclear physics and spectral neutrino transport are juxtaposed to learn about the nature of the three dimensional fluid flow that shapes the explosion. Comparison is also made between the results of simulations in spherical symmetry from several groups, to give ourselves confidence in the understanding derived from this juxtaposition. △ Less

Submitted 17 February, 2016; originally announced February 2016.

Comments: 13 pages, 4 figures, to appear in proceedings of XXXIV Mazurian Lakes Conference on Physics to be published in Acta Physica Polonica B

Journal ref: Acta Physica Polonica B Vol. 47 No. 3 Pg 645 (2016)

Multimessengers from 3D Core-Collapse Supernovae

Authors: Konstantin N. Yakunin, Pedro Marronetti, Anthony Mezzacappa, O. E. Bronson Messer, Eric J. Lentz, Stephen W. Bruenn, W. Raphael Hix, J. Austin Harris

Abstract: We present gravitational wave and neutrino signatures obtained in our first principle 3D core-collapse supernova simulation of 15M non-rotating progenitor with Chimera code. Observations of neutrinos emitted by the forming neutron star, and gravitational waves, which are produced by hydrodynamic instabilities is the only way to get direct information about the supernova engine. Both GW and neutrin… ▽ More We present gravitational wave and neutrino signatures obtained in our first principle 3D core-collapse supernova simulation of 15M non-rotating progenitor with Chimera code. Observations of neutrinos emitted by the forming neutron star, and gravitational waves, which are produced by hydrodynamic instabilities is the only way to get direct information about the supernova engine. Both GW and neutrino signals show different phases of supernova evolution. △ Less

Submitted 21 July, 2015; originally announced July 2015.

Comments: 2 pages, 3 figures, Proceedings, 50th Rencontres de Moriond on Gravitation: La Thuile, Italy, March 21-28, 2014

A Neutrino-Driven Core Collapse Supernova Explosion of a 15 M Star

Authors: Anthony Mezzacappa, Eric J. Lentz, Stephen W. Bruenn, W. Raphael Hix, O. E. Bronson Messer, Eirik Endeve, John M. Blondin, J. Austin Harris, Pedro Marronetti, Konstantin N. Yakunin, Eric J. Lingerfelt

Abstract: We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the "ray-by-ray" approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward tra… ▽ More We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the "ray-by-ray" approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward trajectory at approximately 275 ms after bounce, giving the first indication of a developing explosion in the model. The onset of this shock expansion is delayed relative to our two-dimensional counterpart model, which begins at approximately 200 ms after core bounce. At a time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional model has reached 751 km. Further quantitative analysis of the outcomes in this model must await further development of the post-bounce dynamics and a simulation that will extend well beyond 1 s after stellar core bounce, based on the results for the same progenitor in the context of our two-dimensional, counterpart model. This more complete analysis will determine whether or not the explosion is robust and whether or not observables such as the explosion energy, 56Ni mass, etc. are in agreement with observations. Nonetheless, the onset of explosion in our ab initio three-dimensional multi-physics model with multi-frequency neutrino transport and general relativity is encouraging. △ Less

Submitted 20 July, 2015; originally announced July 2015.

Comments: To appear in the Proceedings of ASTRONUM 2014. arXiv admin note: text overlap with arXiv:1405.7075

Gravitational Wave Signatures of Ab Initio Two-Dimensional Core Collapse Supernova Explosion Models for 12-25 Solar Masses Stars

Authors: Konstantin N. Yakunin, Anthony Mezzacappa, Pedro Marronetti, Shin'ichirou Yoshida, Stephen W. Bruenn, W. Raphael Hix, Eric J. Lentz, O. E. Bronson Messer, J. Austin Harris, Eirik Endeve, John M. Blondin, Eric J. Lingerfelt

Abstract: We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the Chimera code for progenitor masses between 12 and 25 solar masses. All models employ multi-frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutri… ▽ More We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the Chimera code for progenitor masses between 12 and 25 solar masses. All models employ multi-frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and SASI activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway. △ Less

Submitted 21 May, 2015; originally announced May 2015.

Comments: 12 pages, 7 figures

Three-dimensional core-collapse supernova simulated using a 15 $M_\odot$ progenitor

Authors: Eric J. Lentz, Stephen W. Bruenn, W. Raphael Hix, Anthony Mezzacappa, O. E. Bronson Messer, Eirik Endeve, John M. Blondin, J. Austin Harris, Pedro Marronetti, Konstantin N. Yakunin

Abstract: We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 $M_\odot$ progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed by $\sim$100 ms in 3D relative to the 2D counterpart and the growth of the dia… ▽ More We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 $M_\odot$ progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed by $\sim$100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the $\sim$100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions. △ Less

Submitted 16 June, 2015; v1 submitted 19 May, 2015; originally announced May 2015.

Comments: 7 pages, accepted for publication in ApJ Letters, animation available at ChimeraSN.org

Recent Progress on Ascertaining the Core Collapse Supernova Explosion Mechanism

Authors: Anthony Mezzacappa, Stephen W. Bruenn, Eric J. Lentz, W. Raphael Hix, J. Austin Harris, O. E. Bronson Messer, Eirik Endeve, Merek A. Chertkow, John M. Blondin, Pedro Marronetti, Konstantin N. Yakunin

Abstract: We have been working within the fundamental paradigm that core collapse supernovae (CCSNe) may be neutrino driven, since the first suggestion of this by Colgate and White nearly five decades ago. Computational models have become increasingly sophisticated, first in one spatial dimension assuming spherical symmetry, then in two spatial dimensions assuming axisymmetry, and now in three spatial dimen… ▽ More We have been working within the fundamental paradigm that core collapse supernovae (CCSNe) may be neutrino driven, since the first suggestion of this by Colgate and White nearly five decades ago. Computational models have become increasingly sophisticated, first in one spatial dimension assuming spherical symmetry, then in two spatial dimensions assuming axisymmetry, and now in three spatial dimensions with no imposed symmetries. The increase in the number of spatial dimensions has been accompanied by an increase in the physics included in the models, and an increase in the sophistication with which this physics has been modeled. Computation has played an essential role in the development of CCSN theory, not simply for the obvious reason that such multidimensional, multi-physics, nonlinear events cannot possibly be fully captured analytically, but for its role in discovery. In particular, the discovery of the standing accretion shock instability (SASI) through computation about a decade ago has impacted all simulations performed since then. Today, we appear to be at a threshold, where neutrinos, neutrino-driven convection, and the SASI, working together over time scales significantly longer than had been anticipated in the past, are able to generate explosions, and in some cases, robust explosions, in a number of axisymmetric models. But how will this play out in three dimensions? Early results from the first three-dimensional (3D), multi-physics simulation of the "Oak Ridge" group are promising. I will discuss the essential components of today's models and the requirements of realistic CCSN modeling, present results from our one-, two-, and three-dimensional models, place our models in context with respect to other efforts around the world, and discuss short- and long-term next steps. △ Less

Submitted 7 January, 2015; originally announced January 2015.

Advancing Nucleosynthesis in Self-consistent, Multidimensional Models of Core-Collapse Supernovae

Authors: J. Austin Harris, W. Raphael Hix, Merek A. Chertkow, Stephen W. Bruenn, Eric J. Lentz, O. E. Bronson Messer, Anthony Mezzacappa, John M. Blondin, Pedro Marronetti, Konstantin N. Yakunin

Abstract: We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $α$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutr… ▽ More We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $α$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutronization and the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. Limitations such as poor spatial resolution of the tracer particles, estimation of the expansion timescales, and determination of the "mass-cut" at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of these uncertainties on post-processing nucleosynthesis calculations and implications for future models. △ Less

Submitted 31 October, 2014; originally announced November 2014.

Comments: Proceedings of the 13th Symposium on Nuclei in the Cosmos. 7-11 July 2014. Debrecen, Hungary

The Development of Explosions in Axisymmetric Ab Initio Core-Collapse Supernova Simulations of 12-25 $M_\odot$ Stars

Authors: Stephen W. Bruenn, Eric J. Lentz, W. Raphael Hix, Anthony Mezzacappa, James Austin Harris, O. E. Bronson Messer, Eirik Endeve, John M. Blondin, Merek Austin Chertkow, Eric J. Lingerfelt, Pedro Marronetti, Konstantin N. Yakunin

Abstract: We present four ab initio axisymmetric core-collapse supernova simulations for 12, 15, 20, and 25 $M_\odot$ progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 seconds after core bounce and 1 second after material first becomes unbound. Simulations were computed with our Chimera code employing spectral neutrino transport, special and general relativistic tra… ▽ More We present four ab initio axisymmetric core-collapse supernova simulations for 12, 15, 20, and 25 $M_\odot$ progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 seconds after core bounce and 1 second after material first becomes unbound. Simulations were computed with our Chimera code employing spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 second allows explosions to develop more fully and the processes powering the explosions to become more clearly evident. We compute explosion energy estimates, including the binding energy of the stellar envelope outside the shock, of 0.34, 0.88, 0.38, and 0.70 B ($10^{51}$ ergs) and increasing at 0.03, 0.15, 0.19, and 0.52 B s$^{-1}$, respectively, for the 12, 15, 20, and 25 $M_\odot$ models. Three models developed pronounced prolate shock morphologies, while the 20 $M_\odot$ model, though exhibiting lobes and accretion streams like the other models, develops an approximately spherical, off-center shock as the explosion begins and then becomes moderately prolate $\sim$600 ms after bounce. This reduces the explosion energy relative to the other models by reducing mass accretion during the critical explosion power-up phase. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We find that the 12 and 20 $M_\odot$ models have explosion energies comparable to that of the lower range of observed explosion energies while the 15 and 25 $M_\odot$ models are within the range of observed explosion energies, particularly considering the rate at which their explosion energies are increasing. The ejected $^{56}$Ni masses given by our models are all within observational limits as are the proto-neutron star masses and kick velocities. (Truncated) △ Less

Submitted 3 November, 2015; v1 submitted 19 September, 2014; originally announced September 2014.

Comments: 39 pages, 29 figures, accepted for publicatoin in The Astrophysical Journal, additional animations available at ChimeraSN.org

Journal ref: ApJ, 818, 123 (2016)

Two- and Three-Dimensional Multi-Physics Simulations of Core Collapse Supernovae: A Brief Status Report and Summary of Results from the "Oak Ridge" Group

Authors: Anthony Mezzacappa, Stephen W. Bruenn, Eric J. Lentz, W. Raphael Hix, O. E. Bronson Messer, J. Austin Harris, Eric J. Lingerfelt, Eirik Endeve, Konstantin N. Yakunin, John M. Blondin, Pedro Marronetti

Abstract: We summarize the results of core collapse supernova theory from one-, two-, and three-dimensional models and provide a snapshot of the field at this time. We also present results from the "Oak Ridge" group in this context. Studies in both one and two spatial dimensions define the necessary} physics that must be included in core collapse supernova models: a general relativistic treatment of gravity… ▽ More We summarize the results of core collapse supernova theory from one-, two-, and three-dimensional models and provide a snapshot of the field at this time. We also present results from the "Oak Ridge" group in this context. Studies in both one and two spatial dimensions define the necessary} physics that must be included in core collapse supernova models: a general relativistic treatment of gravity (at least an approximate one), spectral neutrino transport, including relativistic effects such as gravitational redshift, and a complete set of neutrino weak interactions that includes state-of-the-art electron capture on nuclei and energy-exchanging scattering on electrons and nucleons. Whether or not the necessarily approximate treatment of this physics in current models that include it is sufficient remains to be determined in the context of future models that remove the approximations. We summarize the results of the Oak Ridge group's two-dimensional supernova models. In particular, we demonstrate that robust neutrino-driven explosions can be obtained. We also demonstrate that our predictions of the explosion energies and remnant neutron star masses are in agreement with observations, although a much larger number of models must be developed before more confident conclusions can be made. We provide preliminary results from our ongoing three dimensional model with the same physics. Finally, we speculate on future outcomes and directions. △ Less

Submitted 27 May, 2014; originally announced May 2014.

Comments: To appear in the proceedings of ASTRONUM 2013

Strong neutrino cooling by cycles of electron capture and $β^-$ decay in neutron star crusts

Authors: Hendrik Schatz, Sanjib Gupta, Peter Möller, Mary Beard, Edward F. Brown, Alex T. Deibel, Leandro R. Gasques, William Raphael Hix, Laurens Keek, Rita Lau, Andrew W. Steiner, Michael Wiescher

Abstract: The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron cap… ▽ More The temperature in the crust of an accreting neutron star, which comprises its outermost kilometer, is set by heating from nuclear reactions at large densities, neutrino cooling, and heat transport from the interior. The heated crust has been thought to affect observable phenomena at shallower depths, such as thermonuclear bursts in the accreted envelope. Here we report that cycles of electron capture and its inverse, $β^-$ decay, involving neutron-rich nuclei at a typical depth of about 150 m, cool the outer neutron star crust by emitting neutrinos while also thermally decoupling the surface layers from the deeper crust. This Urca mechanism has been studied in the context of white dwarfs and Type Ia supernovae, but hitherto was not considered in neutron stars, because previous models computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. This thermal decoupling means that X-ray bursts and other surface phenomena are largely independent of the strength of deep crustal heating. The unexpectedly short recurrence times, of the order of years, observed for very energetic thermonuclear superbursts are therefore not an indicator of a hot crust, but may point instead to an unknown local heating mechanism near the neutron star surface. △ Less

Submitted 9 December, 2013; originally announced December 2013.

Comments: 15 pages with 5 figures. Nature advance online publication, December 1, 2013

Two- and three-dimensional simulations of core-collapse supernovae with CHIMERA

Authors: Eric J. Lentz, Stephen W. Bruenn, J. Austin Harris, Merek Austin Chertkow, W. Raphael Hix, Anthony Mezzacappa, O. E. Bronson Messer, John M. Bondin, Pedro Marronetti, Christopher M. Mauney, Konstantin N. Yakunin

Abstract: Ascertaining the core-collapse supernova mechanism is a complex, and yet unsolved, problem dependent on the interaction of general relativity, hydrodynamics, neutrino transport, neutrino-matter interactions, and nuclear equations of state and reaction kinetics. Ab initio modeling of core-collapse supernovae and their nucleosynthetic outcomes requires care in the coupling and approximations of the… ▽ More Ascertaining the core-collapse supernova mechanism is a complex, and yet unsolved, problem dependent on the interaction of general relativity, hydrodynamics, neutrino transport, neutrino-matter interactions, and nuclear equations of state and reaction kinetics. Ab initio modeling of core-collapse supernovae and their nucleosynthetic outcomes requires care in the coupling and approximations of the physical components. We have built our multi-physics CHIMERA code for supernova modeling in 1-, 2-, and 3-D, using ray-by-ray neutrino transport, approximate general relativity, and detailed neutrino and nuclear physics. We discuss some early results from our current series of exploding 2D simulations and our work to perform computationally tractable simulations in 3D using the "Yin-Yang" grid. △ Less

Submitted 10 January, 2013; v1 submitted 7 January, 2013; originally announced January 2013.

Comments: Proceedings of the 12th Symposium on Nuclei in the Cosmos. 5-12 August 2012. Cairns, Australia. Published online at http://pos.sissa.it/archive/conferences/146/208/NIC%20XII_208.pdf Corrected typo

Axisymmetric Ab Initio Core-Collapse Supernova Simulations of 12-25 M_sol Stars

Authors: Stephen W. Bruenn, Anthony Mezzacappa, W. Raphael Hix, Eric J. Lentz, O. E. Bronson Messer, Eric J. Lingerfelt, John M. Blondin, Eirik Endeve, Pedro Marronetti, Konstantin N. Yakunin

Abstract: We present an overview of four ab initio axisymmetric core-collapse supernova simulations employing detailed spectral neutrino transport computed with our CHIMERA code and initiated from Woosley & Heger (2007) progenitors of mass 12, 15, 20, and 25 M_sol. All four models exhibit shock revival over \sim 200 ms (leading to the possibility of explosion), driven by neutrino energy deposition. Hydrodyn… ▽ More We present an overview of four ab initio axisymmetric core-collapse supernova simulations employing detailed spectral neutrino transport computed with our CHIMERA code and initiated from Woosley & Heger (2007) progenitors of mass 12, 15, 20, and 25 M_sol. All four models exhibit shock revival over \sim 200 ms (leading to the possibility of explosion), driven by neutrino energy deposition. Hydrodynamic instabilities that impart substantial asymmetries to the shock aid these revivals, with convection appearing first in the 12 M_sol model and the standing accretion shock instability (SASI) appearing first in the 25 M_sol model. Three of the models have developed pronounced prolate morphologies (the 20 M_sol model has remained approximately spherical). By 500 ms after bounce the mean shock radii in all four models exceed 3,000 km and the diagnostic explosion energies are 0.33, 0.66, 0.65, and 0.70 Bethe (B = $10^{51}$ ergs) for the 12, 15, 20, and 25 M_sol models, respectively, and are increasing. The three least massive of our models are already sufficiently energetic to completely unbind the envelopes of their progenitors (i.e., to explode), as evidenced by our best estimate of their explosion energies, which first become positive at 320, 380, and 440 ms after bounce. By 850 ms the 12 M_sol diagnostic explosion energy has saturated at 0.38 B, and our estimate for the final kinetic energy of the ejecta is \sim 0.3 B, which is comparable to observations for lower-mass progenitors. △ Less

Submitted 12 March, 2013; v1 submitted 7 December, 2012; originally announced December 2012.

Comments: Accepted for publication in The Astrophysical Journal Letters

Hydrodynamic Studies of the Evolution of Recurrent, Symbiotic, and Dwarf Novae: The White Dwarf Components are Growing in Mass

Authors: S. Starrfield, F. X. Timmes, C. Iliadis, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely lum… ▽ More Symbiotic binaries are systems containing white dwarfs (WDs) and red giants. Symbiotic novae are those systems in which thermonuclear eruptions occur on the WD components. These are to be distinguished from events driven by accretion disk instabilities analogous to dwarf novae eruptions in cataclysmic variable outbursts. Another class of symbiotic systems are those in which the WD is extremely luminous and it seems likely that quiescent nuclear burning is ongoing on the accreting WD. A fundamental question is the secular evolution of the WD. Do the repeated outbursts or quiescent burning in these accreting systems cause the WD to gain or lose mass? If it is gaining mass, can it eventually reach the Chandrasekhar Limit and become a supernova (a SN Ia if it can hide the hydrogen and helium in the system)? In order to better understand these systems, we have begun a new study of the evolution of Thermonuclear Runaways (TNRs) in the accreted envelopes of WDs using a variety of initial WD masses, luminosities and mass accretion rates. We use our 1-D hydro code, NOVA, which includes the new convective algorithm of Arnett, Meakin and Young, the Hix and Thielemann nuclear reaction solver, the Iliadis reaction rate library, the Timmes equation of state, and the OPAL opacities. We assume a solar composition (Lodders abundance distribution) and do not allow any mixing of accreted material with core material. This assumption strongly influences our results. We report here (1) that the WD grows in mass for all simulations so that canonical `steady burning' does not occur, and (2) that only a small fraction of the accreted matter is ejected in some (but not all) simulations. We also find that the accreting systems, before thermonuclear runaway, are too cool to be seen in X-ray searches for SN Ia progenitors. △ Less

Submitted 26 November, 2012; originally announced November 2012.

Comments: 11 pages 10 figures in the 2011 Asiago Symbiotic Workshop

Journal ref: Baltic Astronomy, Vol 21, pages 76 - 87, 2012

Theoretical Studies of Accretion of Matter onto White Dwarfs and the Single Degenerate Scenario for Supernovae of Type Ia

Authors: S. Starrfield, C. Iliadis, F. X. Timmes, W. R. Hix, W. D. Arnett, C. Meakin, W. M. Sparks

Abstract: We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and s… ▽ More We present a brief summary of the Single Degenerate Scenario for the progenitors of Type Ia Supernovae in which it is assumed that a low mass carbon-oxygen white dwarf is growing in mass as a result of accretion from a secondary star in a close binary system. Recent hydrodynamic simulations of accretion of solar material onto white dwarfs without mixing always produce a thermonuclear runaway and steady burning does not occur. For a broad range in WD mass (0.4 Solar masses to 1.35 Solar Masses), the maximum ejected material occurs for the 1.25 Solar Mass sequences and then decreases as the white dwarf mass decreases. Therefore, the white dwarfs are growing in mass as a consequence of the accretion of solar material and as long as there is no mixing of accreted material with core material. In contrast, a thermonuclear runaway in the accreted hydrogen-rich layers on the low luminosity WDs in close binary systems where mixing of core matter with accreted material has occurred is the outburst mechanism for Classical, Recurrent, and Symbiotic novae. The differences in characteristics of these systems is likely the WD mass and mass accretion rate. The high levels of enrichment of CN ejecta in elements ranging from carbon to sulfur confirm that there is dredge-up of matter from the core of the WD and enable them to contribute to the chemical enrichment of the interstellar medium. Therefore, studies of CNe can lead to an improved understanding of Galactic nucleosynthesis, some sources of pre-solar grains, and the Extragalactic distance scale. The characteristics of the outburst depend on the white dwarf mass, luminosity, mass accretion rate, and the chemical composition of both the accreting material and WD material. The properties of the outburst also depends on when, how, and if the accreted layers are mixed with the WD core and the mixing mechanism is still unknown. △ Less

Submitted 22 October, 2012; originally announced October 2012.

Comments: 25 Pages, Bulletin of the Astronomical Society of India (BASI) in press

Interplay of Neutrino Opacities in Core-collapse Supernova Simulations

Authors: Eric J. Lentz, Anthony Mezzacappa, O. E. Bronson Messer, W. Raphael Hix, Stephen W. Bruenn

Abstract: We have conducted a series of numerical experiments using spherically symmetric, general relativistic, neutrino radiation hydrodynamics with the code Agile-BOLTZTRAN to examine the effects of modern neutrino opacities on the development of supernova simulations. We test the effects of opacities by removing opacities or by undoing opacity improvements for individual opacities and groups of opacitie… ▽ More We have conducted a series of numerical experiments using spherically symmetric, general relativistic, neutrino radiation hydrodynamics with the code Agile-BOLTZTRAN to examine the effects of modern neutrino opacities on the development of supernova simulations. We test the effects of opacities by removing opacities or by undoing opacity improvements for individual opacities and groups of opacities. We find that improvements to electron capture (EC) on nuclei, namely EC on an ensemble of nuclei using modern nuclear structure models rather than the simpler independent-particle approximation (IPA) for EC on a mean nucleus, plays the most important role during core collapse of all tested neutrino opacities. Low-energy neutrinos emitted by modern nuclear EC preferentially escape during collapse without the energy downscattering on electrons required to enhance neutrino escape and deleptonization for the models with IPA nuclear EC. During shock breakout the primary influence on the emergent neutrinos arises from NIS on electrons. For the accretion phase, non-isoenergetic scattering on free nucleons and pair emission by $e^+e^-$ annihilation have the largest impact on the neutrino emission and shock evolution. Other opacities evaluated, including nucleon--nucleon bremsstrahlung and especially neutrino--positron scattering, have little measurable impact on neutrino emission or shock dynamics. Modern treatments of nuclear electron capture, $e^+e^-$-annihilation pair emission, and non-isoenergetic scattering on electrons and free nucleons are critical elements of core-collapse simulations of all dimensionality. △ Less

Submitted 2 October, 2012; v1 submitted 5 June, 2012; originally announced June 2012.

Comments: 12 pages, 15 figures, revised, accepted to ApJ

Journal ref: ApJ, 2012, 760, 94

Explicit Integration of Extremely-Stiff Reaction Networks: Partial Equilibrium Methods

Authors: M. W. Guidry, J. J. Billings, W. R. Hix

Abstract: In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard implicit methods in accuracy and speed for extremely stiff systems. However, we also showed that these explicit methods remain accurate but are no longer competitive… ▽ More In two preceding papers we have shown that, when reaction networks are well-removed from equilibrium, explicit asymptotic and quasi-steady-state approximations can give algebraically-stabilized integration schemes that rival standard implicit methods in accuracy and speed for extremely stiff systems. However, we also showed that these explicit methods remain accurate but are no longer competitive in speed as the network approaches equilibrium. In this paper we analyze this failure and show that it is associated with the presence of fast equilibration timescales that neither asymptotic nor quasi-steady-state approximations are able to remove efficiently from the numerical integration. Based on this understanding, we develop a partial equilibrium method to deal effectively with the approach to equilibrium and show that explicit asymptotic methods, combined with the new partial equilibrium methods, give an integration scheme that plausibly can deal with the stiffest networks, even in the approach to equilibrium, with accuracy and speed competitive with that of implicit methods. Thus we demonstrate that such explicit methods may offer alternatives to implicit integration of even extremely stiff systems, and that these methods may permit integration of much larger networks than have been possible before in a number of fields. △ Less

Submitted 20 December, 2011; v1 submitted 20 December, 2011; originally announced December 2011.

Comments: Updated references

Journal ref: Comput. Sci. Disc. 6, 015003 (2013)

Explicit Integration of Extremely-Stiff Reaction Networks: Asymptotic Methods

Authors: M. W. Guidry, R. Budiardja, E. Feger, J. J. Billings, W. R. Hix, O. E. B. Messer, K. J. Roche, E. McMahon, M. He

Abstract: We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends on whether the system is well-removed from equilibrium or near equilibrium. This paper introduces a quantitative distinction between these two regimes and addr… ▽ More We show that, even for extremely stiff systems, explicit integration may compete in both accuracy and speed with implicit methods if algebraic methods are used to stabilize the numerical integration. The required stabilizing algebra depends on whether the system is well-removed from equilibrium or near equilibrium. This paper introduces a quantitative distinction between these two regimes and addresses the former case in depth, presenting explicit asymptotic methods appropriate when the system is extremely stiff but only weakly equilibrated. A second paper examines quasi-steady-state methods as an alternative to asymptotic methods in systems well away from equilibrium and a third paper extends these methods to equilibrium conditions in extremely stiff systems using partial equilibrium methods. All three papers present systematic evidence for timesteps competitive with implicit methods. Because an explicit method can execute a timestep faster than an implicit method, algebraically-stabilized explicit algorithms might permit integration of larger networks than have been feasible before in various disciplines. △ Less

Submitted 20 December, 2011; v1 submitted 20 December, 2011; originally announced December 2011.

Comments: Updated references

Journal ref: Comput. Sci. Disc. 6, 015001 (2013)

On the Requirements for Realistic Modeling of Neutrino Transport in Simulations of Core-Collapse Supernovae

Authors: Eric J. Lentz, Anthony Mezzacappa, O. E. Bronson Messer, Matthias Liebendörfer, W. Raphael Hix, Stephen W. Bruenn

Abstract: We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-BOLTZTRAN to examine the effects of several approximations used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of these approximations quantitatively by removing, or substituting for, the pieces of… ▽ More We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-BOLTZTRAN to examine the effects of several approximations used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of these approximations quantitatively by removing, or substituting for, the pieces of supernova physics of interest. These approximations include: (1) using Newtonian versus general relativistic gravity, hydrodynamics, and transport; (2) using a reduced set of weak interactions, including the omission of non-isoenergetic neutrino scattering, versus the current state-of-the-art; and (3) omitting the velocity-dependent terms, or observer corrections, from the neutrino Boltzmann kinetic equation. We demonstrate that each of these changes has noticeable effects on the outcomes of our simulations. Of these, we find that the omission of observer corrections is particularly detrimental to the potential for neutrino-driven explosions and exhibits a failure to conserve lepton number. Finally, we discuss the impact of these results on our understanding of current, and the requirements for future, multidimensional models. △ Less

Submitted 15 December, 2011; originally announced December 2011.

Comments: 12 pages, 9 figures, Accepted for publication in ApJ

Journal ref: 2012, ApJ, 747, 73

Time-of-flight mass measurements for nuclear processes in neutron star crusts

Authors: A. Estrade, M. Matos, H. Schatz, A. M. Amthor, D. Bazin, M. Beard, A. Becerril, E. F. Brown, R. Cyburt, T. Elliot, A. Gade, D. Galaviz, S. George, S. S. Gupta, W. R. Hix, R. Lau, G. Lorusso, P. Moller, J. Pereira, M. Portillo, A. M. Rogers, D. Shapira, E. Smith, A. Stolz, M. Wallace , et al. (1 additional authors not shown)

Abstract: The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were dete… ▽ More The location of electron capture heat sources in the crust of accreting neutron stars depends on the masses of extremely neutron-rich nuclei. We present first results from a new implementation of the time-of-flight technique to measure nuclear masses of rare isotopes at the National Superconducting Cyclotron Laboratory. The masses of 16 neutron-rich nuclei in the scandium -- nickel range were determined simultaneously, improving the accuracy compared to previous data in 12 cases. The masses of $^{61}${V}, $^{63}${Cr}, $^{66}${Mn}, and $^{74}${Ni} were measured for the first time with mass excesses of $-30.510(890)$ MeV, $-35.280(650)$ MeV, $-36.900(790)$ MeV, and $-49.210(990)$ MeV, respectively. With the measurement of the $^{66}$Mn mass, the locations of the two dominant electron capture heat sources in the outer crust of accreting neutron stars that exhibit superbursts are now experimentally constrained. We find that the location of the $^{66}$Fe$\rightarrow^{66}$Mn electron capture transition occurs significantly closer to the surface than previously assumed because our new experimental Q-value is 2.1 MeV (2.6$σ$) smaller than predicted by the FRDM mass model. △ Less

Submitted 23 September, 2011; originally announced September 2011.

Comments: 4 pages, 3 figures

Evaluating nuclear physics inputs in core-collapse supernova models

Authors: Eric J. Lentz, W. Raphael Hix, Mark L. Baird, O. E. Bronson Messer, Anthony Mezzacappa

Abstract: Core-collapse supernova models depend on the details of the nuclear and weak interaction physics inputs just as they depend on the details of the macroscopic physics (transport, hydrodynamics, etc.), numerical methods, and progenitors. We present preliminary results from our ongoing comparison studies of nuclear and weak interaction physics inputs to core collapse supernova models using the spheri… ▽ More Core-collapse supernova models depend on the details of the nuclear and weak interaction physics inputs just as they depend on the details of the macroscopic physics (transport, hydrodynamics, etc.), numerical methods, and progenitors. We present preliminary results from our ongoing comparison studies of nuclear and weak interaction physics inputs to core collapse supernova models using the spherically-symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-Boltztran. We focus on comparisons of the effects of the nuclear EoS and the effects of improving the opacities, particularly neutrino--nucleon interactions. △ Less

Submitted 30 December, 2010; originally announced January 2011.

Comments: 5 pages, 2 figures, presented at Nuclei in the Cosmos XI, Heidelberg, Germany

Core Collapse Supernovae using CHIMERA: Gravitational Radiation from Non-Rotating Progenitors

Authors: Konstantin Yakunin, Pedro Marronetti, Anthony Mezzacappa, Stephen W Bruenn, Ching-Tsai Lee, Merek A Chertkow, W Raphael Hix, John M Blondin, Eric J Lentz, O E Bronson Messer, Shin Yoshida

Abstract: The CHIMERA code is a multi-dimensional multi-physics engine dedicated primarily to the simulation of core collapse supernova explosions. One of the most important aspects of these explosions is their capacity to produce gravitational radiation that is detectable by Earth-based laser-interferometric gravitational wave observatories such as LIGO and VIRGO. We present here preliminary gravitational… ▽ More The CHIMERA code is a multi-dimensional multi-physics engine dedicated primarily to the simulation of core collapse supernova explosions. One of the most important aspects of these explosions is their capacity to produce gravitational radiation that is detectable by Earth-based laser-interferometric gravitational wave observatories such as LIGO and VIRGO. We present here preliminary gravitational signatures of two-dimensional models with non-rotating progenitors. These simulations exhibit explosions, which are followed for more than half a second after stellar core bounce. △ Less

Submitted 16 November, 2010; originally announced November 2010.

Comments: 4 pages, 1 figure, SciDAC 2010 Conference Proceedings

Gravitational Waves from Core Collapse Supernovae

Authors: Konstantin N Yakunin, Pedro Marronetti, Anthony Mezzacappa, Stephen W Bruenn, Ching-Tsai Lee, Merek A Chertkow, W Raphael Hix, John M Blondin, Eric J Lentz, O E Bronson Messer, Shin'ichirou Yoshida

Abstract: We present the gravitational wave signatures for a suite of axisymmetric core collapse supernova models with progenitors masses between 12 and 25 solar masses. These models are distinguished by the fact they explode and contain essential physics (in particular, multi-frequency neutrino transport and general relativity) needed for a more realistic description. Thus, we are able to compute complete… ▽ More We present the gravitational wave signatures for a suite of axisymmetric core collapse supernova models with progenitors masses between 12 and 25 solar masses. These models are distinguished by the fact they explode and contain essential physics (in particular, multi-frequency neutrino transport and general relativity) needed for a more realistic description. Thus, we are able to compute complete waveforms (i.e., through explosion) based on non-parameterized, first-principles models. This is essential if the waveform amplitudes and time scales are to be computed more precisely. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models. △ Less

Submitted 14 July, 2010; v1 submitted 5 May, 2010; originally announced May 2010.

Comments: 10 pages, 5 figures

Journal ref: Class.Quant.Grav.27:194005,2010

2D and 3D Core-Collapse Supernovae Simulation Results Obtained with the CHIMERA Code

Authors: S. W. Bruenn, A. Mezzacappa, W. R. Hix, J. M. Blondin, P. Marronetti, O. E. B. Messer, C. J. Dirk, S. Yoshida

Abstract: Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with enough realism that the explosion mechanism, long a mystery, may soon be delineated. We briefly describe the CHIMERA code, a supernova code we have developed… ▽ More Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with enough realism that the explosion mechanism, long a mystery, may soon be delineated. We briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We then describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 20,000 km. We also briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor. △ Less

Submitted 25 February, 2010; originally announced February 2010.

Comments: 5 pages, 3 figures

Journal ref: SciDAC 2009 Scientific Discovery through Advanced Computing, 14-18 June, San Diego, California, USA. AIP Conference Proceedings, Volume 180 (2009) 012018, pp. 1-5

Mechanisms of Core-Collapse Supernovae & Simulation Results from the CHIMERA Code

Authors: S. W. Bruenn, A. Mezzacappa, W. R. Hix, J. M. Blondin, P. Marronetti, O. E. B. Messer, C. J. Dirk, S. Yoshida

Abstract: Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to s… ▽ More Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to some key insights which may clarify the picture in the not too distant future. Here we briefly review the current status of the three explosion mechanisms (acoustic, MHD, and neutrino heating) that are currently under active investigation, concentrating on the neutrino heating mechanism as the one most likely responsible for producing explosions from progenitors in the mass range ~10 to ~25 solar masses. We then briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We finally describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 10,000 km. We finally very briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor. △ Less

Submitted 25 February, 2010; originally announced February 2010.

Comments: 10 pages, 4 figures

Journal ref: PROBING STELLAR POPULATIONS OUT TO THE DISTANT UNIVERSE: CEFALU 2008, Proceedings of the International Conference. AIP Conference Proceedings, Volume 1111, pp. 593-601 (2009)

R-process nucleosynthesis calculations with complete nuclear physics input

Authors: I. Petermann, A. Arcones, A. Kelić, K. Langanke, G. Martínez-Pinedo, K. -H. Schmidt, W. R. Hix, I. Panov, T. Rauscher, F. -K. Thielemann, N. Zinner

Abstract: The r-process constitutes one of the major challenges in nuclear astrophysics. Its astrophysical site has not yet been identified but there is observational evidence suggesting that at least two possible sites should contribute to the solar system abundance of r-process elements and that the r-process responsible for the production of elements heavier than Z=56 operates quite robustly producing… ▽ More The r-process constitutes one of the major challenges in nuclear astrophysics. Its astrophysical site has not yet been identified but there is observational evidence suggesting that at least two possible sites should contribute to the solar system abundance of r-process elements and that the r-process responsible for the production of elements heavier than Z=56 operates quite robustly producing always the same relative abundances. From the nuclear-physics point of view the r-process requires the knowledge of a large number of reaction rates involving exotic nuclei. These include neutron capture rates, beta-decays and fission rates, the latter for the heavier nuclei produced in the r-process. We have developed for the first time a complete database of reaction rates that in addition to neutron-capture rates and beta-decay half-lives includes all possible reactions that can induce fission (neutron-capture, beta-decay and spontaneous fission) and the corresponding fission yields. In addition, we have implemented these reaction rates in a fully implicit reaction network. We have performed r-process calculations for the neutrino-driven wind scenario to explore whether or not fission can contribute to provide a robust r-process pattern. △ Less

Submitted 4 December, 2008; originally announced December 2008.

Journal ref: Phys.Rev.C81:014308,2010

Tidally-induced thermonuclear Supernovae

Authors: S. Rosswog, E. Ramirez-Ruiz, W. R. Hix

Abstract: We discuss the results of 3D simulations of tidal disruptions of white dwarfs by moderate-mass black holes as they may exist in the cores of globular clusters or dwarf galaxies. Our simulations follow self-consistently the hydrodynamic and nuclear evolution from the initial parabolic orbit over the disruption to the build-up of an accretion disk around the black hole. For strong enough encounter… ▽ More We discuss the results of 3D simulations of tidal disruptions of white dwarfs by moderate-mass black holes as they may exist in the cores of globular clusters or dwarf galaxies. Our simulations follow self-consistently the hydrodynamic and nuclear evolution from the initial parabolic orbit over the disruption to the build-up of an accretion disk around the black hole. For strong enough encounters (pericentre distances smaller than about 1/3 of the tidal radius) the tidal compression is reversed by a shock and finally results in a thermonuclear explosion. These explosions are not restricted to progenitor masses close to the Chandrasekhar limit, we find exploding examples throughout the whole white dwarf mass range. There is, however, a restriction on the masses of the involved black holes: black holes more massive than $2\times 10^5$ M$_\odot$ swallow a typical 0.6 M$_\odot$ dwarf before their tidal forces can overwhelm the star's self-gravity. Therefore, this mechanism is characteristic for black holes of moderate masses. The material that remains bound to the black hole settles into an accretion disk and produces an X-ray flare close to the Eddington limit of $L_{\rm Edd} \simeq 10^{41} {\rm erg/s} M_{\rm bh}/1000 M$_\odot$), typically lasting for a few months. The combination of a peculiar thermonuclear supernova together with an X-ray flare thus whistle-blows the existence of such moderate-mass black holes. The next generation of wide field space-based instruments should be able to detect such events. △ Less

Submitted 13 November, 2008; originally announced November 2008.

Comments: 8 pages, 2 figures, EuroWD08

Journal ref: J.Phys.Conf.Ser.172:012036,2009

The Effects of the pep Nuclear Reaction and Other Improvements in the Nuclear Reaction Rate Library on Simulations of the Classical Nova Outburst

Authors: S. Starrfield, C. Iliadis, W. R. Hix, F. X. Timmes, W. M. Sparks

Abstract: We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations ove… ▽ More We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. In addition, as a result of our improvements, we discovered that the pep reaction was not included in our previous studies of CN explosions. Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed $10^4$ gm cm$^{-3}$ and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. The effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper. △ Less

Submitted 2 November, 2008; originally announced November 2008.

Comments: Accepted to the Astrophysical Journal

Journal ref: Astrophys.J.692:1532-1542,2009