Search NASA⌕ Search

Engineering topics

Hix, William Raphael

Publications and source records attributed to Hix, William Raphael.

Nuclear Reactions in the Crusts of Accreting Neutron Stars

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 $\rho =2\times {10}^{12}\,{\rm{g}}\,{\mathrm{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 $\rho =2\times {10}^{12}\,{\rm{g}}\,{\mathrm{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 Q imp ≈ 20 owing 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. As a result, 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.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrodynamic Simulations of Classical Novae: Accretion onto CO White Dwarfs as SN Ia Progenitors

We report on our continuing studies of Classical Nova explosions by following the evolution of thermonuclear runaways (TNRs) on carbon-oxygen (CO) white dwarfs (WDs). We have varied both the mass of the WD and the composition of the accreted material. Rather than assuming that the material has mixed from the beginning, we now rely on the results of the multidimensional (multi-D) studies of mixing as a consequence of the TNRs in WDs that accreted only Solar matter. The multi-D studies find that mixing with the core occurs after the TNR is well underway and reach enrichment levels in agreement with observations of the ejecta abundances. We report on 3 studies in this paper. First, simulations in which we accrete only Solar matter with NOVA (our 1-D, fully implicit, hydro code). Second, we use MESA for similar studies in which we accrete only Solar material and compare the results. Third, we accrete Solar matter until the TNR is initiated and then switch the composition in the accreted layers to a mixed composition: either 25% core and 75% Solar or 50% core and 50% Solar. The amount of accreted material is inversely proportional to the initial 12 C abundance. Thus, accreting Solar material results in more material to fuel the outburst - much larger than in the earlier studies where mixed materials were used from the beginning. We tabulate the amount of ejected gases, their velocities, and abundances. We predict the amount of 7 Li and 7 Be produced and ejected by the explosion and compare our predictions to our Large Binocular Telescope (LBT) high dispersion spectra which determined the abundance of 7 Li in nova V5668 Sgr. Finally, many of these simulations eject significantly less mass than accreted and, therefore, the WD is growing in mass toward the Chandrasekhar Limit.

79 ASTRONOMY AND ASTROPHYSICS↗

White paper on nuclear astrophysics and low energy nuclear physics Part 1: Nuclear astrophysics

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 also 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). Our white paper is informed 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. Answers to long standing key questions are well within reach in the coming decade because of the developments outlined in this white paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The development of explosions in axisymmetric ab initio core-collapse supernova simulations of 12–25 M{sub ⊙} STARS

We present four ab initio axisymmetric core-collapse supernova simulations initiated from 12, 15, 20, and 25 M{sub ⊙} zero-age main sequence progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 s after core bounce and 1 s after material first becomes unbound. These simulations were computed with our Chimera code employing RbR spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 s after core bounce allows the explosions to develop more fully and the processes involved in powering the explosions to become more clearly evident. We compute explosion energy estimates, including the negative gravitational binding energy of the stellar envelope outside the expanding shock, of 0.34, 0.88, 0.38, and 0.70 Bethe (B ≡ 10{sup 51} erg) and increasing at 0.03, 0.15, 0.19, and 0.52 B s{sup −1}, respectively, for the 12, 15, 20, and 25 M{sub ⊙} models at the endpoint of this report. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We discuss how the explosion energies may be subject to stochastic variations as exemplfied by the effect of the explosion geometry of the 20 M{sub ⊙} model in reducing its explosion energy. We compute the proto-neutron star masses and kick velocities. We compare our results for the explosion energies and ejected {sup 56}Ni masses against some observational standards despite the large error bars in both models and observations.

79 ASTRONOMY AND ASTROPHYSICS↗

The Multi-Dimensional Character of Core-Collapse Supernovae

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.

Hix, William Raphael↗

Multidimensional simulations of core-collapse supernovae and

Core-collapse supernovae (CCSNe), the culmination of massive stellar evolution, are 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 which is highly hydrodynamically-unstable. Simulations of increasing sophistication show 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 three dimensional fluid flows. Regrettably, much of our understanding of the nucleosynthesis that occurs in these explosions, and their impact on galactic chemical evolution, is based on spherically symmetric simulations with parameterized explosions, ignoring much that has been learned about the central engine of these supernovae over the past two decades. Here in this paper we discuss recent results from two-dimensional CCSN simulations using our CHIMERA code, as well as ongoing three-dimensional simulations, and discuss how the multidimensional character of the explosions directly impacts the nucleosynthesis and other observables of core-collapse supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

White Paper on Nuclear Astrophysics and Low Energy Nuclear Physics Part 1: Nuclear Astrophysics

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 AM 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.

Nuclear astrophysics; White paper; Nucleosynthesis↗

Essential ingredients in core-collapse supernovae

Carrying 10{sup 44} joules of kinetic energy and a rich mix of newly synthesized atomic nuclei, core-collapse supernovae are the preeminent foundries of the nuclear species which make up our solar system and ourselves. Signaling the inevitable death of a massive star, and the birth of a neutron star or black hole, core-collapse supernovae combine physics over a wide range in spatial scales, from kilometer-sized hydrodynamic motions (eventually growing to gigameter scale) down to femtometer-scale nuclear reactions. We will discuss our emerging understanding of the convectively-unstable, neutrino-driven explosion mechanism, based on increasingly realistic neutrino radiation hydrodynamic simulations that include progressively better nuclear and particle physics. Multi-dimensional models with spectral neutrino transport from several research groups, which slowly develop successful explosions for a range of progenitors, have recently motivated changes in our understanding of the neutrino reheating mechanism. In a similar fashion, improvements in nuclear physics, most notably explorations of weak interactions on nuclei and the nuclear equation of state, continue to refine our understanding of the births of neutron stars and the supernovae that result. Recent progress on both the macroscopic and microscopic effects that affect core-collapse supernovae are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

AXISYMMETRIC AB INITIO CORE-COLLAPSE SUPERNOVA SIMULATIONS OF 12-25 M{sub Sun} STARS

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 and Heger progenitors of mass 12, 15, 20, and 25 M{sub Sun }. All four models exhibit shock revival over {approx}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{sub Sun} model and the standing accretion shock instability appearing first in the 25 M{sub Sun} model. Three of the models have developed pronounced prolate morphologies (the 20 M{sub Sun} model has remained approximately spherical). By 500 ms after bounce the mean shock radii in all four models exceed 3000 km and the diagnostic explosion energies are 0.33, 0.66, 0.65, and 0.70 Bethe (B = 10{sup 51} erg) for the 12, 15, 20, and 25 M{sub Sun} 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{sub Sun} diagnostic explosion energy has saturated at 0.38 B, and our estimate for the final kinetic energy of the ejecta is {approx}0.3 B, which is comparable to observations for lower mass progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗

Strong neutrino cooling by cycles of electron capture and decay in neutron star crusts

The temperature in the crust of an accreting neutron star, which comprises its outermost kilometre, 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 metres, 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 dwarfs13 and type Ia supernovae, but hitherto was not considered in neutron stars, because previous models1, 2 computed the crust reactions using a zero-temperature approximation and assumed that only a single nuclear species was present at any given depth. The 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.

Schatz, Hendrik↗

Explicit Integration of Extremely Stiff Reaction Networks: Asymptotic Methods

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 stabilizing algebra differs for systems well removed from equilibrium and those 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 [1] examines quasi-steady-state methods as an alternative to asymptotic methods in systems well away from equilibrium and a third paper [2] 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 explicit methods can execute a timestep faster than an implicit method, our results imply that algebraically stabilized explicit algorithms may offer a means to integration of larger networks than have been feasible previously in various disciplines.

Guidry, Mike W↗

Explicit Integration of Extremely Stiff Reaction Networks: Partial Equilibrium Methods

In two preceding papers [1,2] 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 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 algebraically stabilized 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 feasible previously in a variety of fields.

Guidry, Mike W↗

Core-collapse Supernovae

Marking the inevitable death of a massive star, and the birth of a neutron star or black hole, core-collapse supernovae bring together physics at a wide range in spatial scales, from kilometer-sized hydrodynamic motions (growing to gigameter scale) down to femtometer scale nuclear reactions. Carrying 10$^{51}$ ergs of kinetic energy and a rich-mix of newly synthesized atomic nuclei, core-collapse supernovae are the preeminent foundries of the nuclear species which make up ourselves and our solar system. We will discuss our emerging understanding of the convectively unstable, neutrino-driven explosion mechanism, based on increasingly realistic neutrino-radiation hydrodynamic simulations that include progressively better nuclear and particle physics. Recent multi-dimensional models with spectral neutrino transport from several research groups, which slowly develop successful explosions for a range of progenitors, have motivated changes in our understanding of the neutrino reheating mechanism. In a similar fashion, improvements in nuclear physics, most notably explorations of weak interactions on nuclei and the nuclear equation of state, continue to refine our understanding of how supernovae explode. Recent progress on both the macroscopic and microscopic effects that affect core-collapse supernovae are discussed.

Hix, William Raphael↗

Nucleosynthesis in Thermonuclear Supernovae

We review our understanding of the nucleosynthesis that occurs in thermonuclear supernovae and their contribution to Galactic Chemical evolution. We discuss the prospects to improve the modeling of the nucleosynthesis within simulations of these events.

Claudia, Travaglio↗