-
Enhanced extra mixing in low-mass stars approaching the RGB tip and the problem of Li-rich red-clump stars
Authors:
Pavel A. Denissenkov,
Simon Blouin,
Falk Herwig,
Jacob Stott,
Paul R. Woodward
Abstract:
A few percent of red giants are enriched in Lithium with $A(\mathrm{Li}) > 1.5$. The evolutionary phase of the Li-rich red giants has remained uncertain because they could be placed both on the red-giant branch (RGB) near the bump luminosity and in the red clump (RC) region. However, thanks to asteroseismology, it has been found that most of them are actually RC stars. Starting at the bump luminos…
▽ More
A few percent of red giants are enriched in Lithium with $A(\mathrm{Li}) > 1.5$. The evolutionary phase of the Li-rich red giants has remained uncertain because they could be placed both on the red-giant branch (RGB) near the bump luminosity and in the red clump (RC) region. However, thanks to asteroseismology, it has been found that most of them are actually RC stars. Starting at the bump luminosity, RGB progenitors of the RC stars experience extra mixing in the radiative zone separating the H-burning shell from the convective envelope followed by a series of convective He-shell flashes at the RGB tip, known as the He-core flash. Therefore, the He-core flash was proposed to cause fast extra mixing in the stars at the RGB tip that is needed for the Cameron-Fowler mechanism to produce Li. Alternatively, we propose that the RGB stars are getting enriched in Li by the same extra mixing that starts at the bump luminosity and initially leads to a decrease of the surface Li abundance but that is getting enhanced and begins to produce Li when the stars are approaching the RGB tip. We discuss five mechanisms of the RGB extra mixing, namely, the joint operation of rotation-driven meridional circulation and turbulent diffusion, the Azimuthal Magneto-Rotational Instability (AMRI), thermohaline convection, buoyancy of magnetic flux tubes, and internal gravity waves, and, based on results of (magneto-) hydrodynamics simulations, come to the conclusion that it is the mechanism of the AMRI that is most likely to support our hypothesis.
△ Less
Submitted 8 September, 2023;
originally announced September 2023.
-
3D hydrodynamics simulations of a 3 $M_{\odot}$ core-helium burning star
Authors:
Simon Blouin,
Falk Herwig,
Huaqing Mao,
Pavel Denissenkov,
Paul R. Woodward
Abstract:
The inner structure of core-helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective laye…
▽ More
The inner structure of core-helium burning (CHeB) stars remains uncertain due to the yet unknown nature of mixing at the boundary of their cores. Large convective cores beyond a bare Schwarzschild model are favoured both from theoretical arguments and from asteroseismological constraints. However, the exact nature of this extra mixing, and in particular the possible presence of semiconvective layers, is still debated. In this work, we approach this problem through a new avenue by performing the first full-sphere 3D hydrodynamics simulations of the interiors of CHeB stars. We use the PPMstar explicit gas dynamics code to simulate the inner 0.45 $M_{\odot}$ of a 3 $M_{\odot}$ CHeB star. Simulations are performed using different Cartesian grid resolutions (768$^3$, 1152$^3$ and 1728$^3$) and heating rates. We use two different initial states, one based on MESA's predictive mixing scheme (which significantly extends the core beyond the Schwarzschild boundary) and one based on the convective premixing approach (which exhibits a semiconvective interface). The general behaviour of the flow in the convective core and in the stable envelope (where internal gravity waves are observed) is consistent with our recent simulations of core convection in massive main-sequence stars, and so are the various luminosity scaling relations. The semiconvective layers are dominated by strong internal gravity waves that do not produce measurable species mixing, but overshooting motions from the convective core gradually homogenize the semiconvective interface. This process can possibly completely erase the semiconvective layers, which would imply that CHeB stars do not harbour a semiconvection zone.
△ Less
Submitted 17 November, 2023; v1 submitted 21 August, 2023;
originally announced August 2023.
-
3D hydrodynamic simulations of massive main-sequence stars. III. The effect of radiation pressure and diffusion leading to a 1D equilibrium model
Authors:
Huaqing Mao,
Paul Woodward,
Falk Herwig,
Pavel A. Denissenkov,
Simon Blouin,
William Thompson
Abstract:
We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 $\mathrm{M}_\odot$ star in the early phase of the main-sequence. We use the explicit gas-dynamics code $\texttt{PPMstar}$ which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities…
▽ More
We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25 $\mathrm{M}_\odot$ star in the early phase of the main-sequence. We use the explicit gas-dynamics code $\texttt{PPMstar}$ which tracks two fluids and includes radiation pressure and radiative diffusion. Multiple series of simulations with different luminosities and radiative thermal conductivities are presented. The entrainment rate at the convective boundary, internal gravity waves in and above the boundary region, and the approach to dynamical equilibrium shortly after a few convective turnovers are investigated. From the results of these simulations we extrapolate to find the entrainment rate at the nominal heating rate and thermal diffusion given by the $\texttt{MESA}$ stellar evolution model on which the 3D stratification is based. Further, to study the effect of radiative diffusion on the thermal timescale, we perform very long simulations accelerated by 10000 times their nominal luminosities. In these simulations the growing penetrative convection reduces the initially unrealistically large entrainment. This reduction is enabled by a spatial separation that develops between the entropy gradient and the composition gradient. The convective boundary moves outward much more slowly at the end of these simulations. Finally, we present a method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. This method is intended to be ultimately deployed in 1D stellar evolution calculations and is based on the properties of penetrative convection in our simulations carried forward through the local thermal timescale.
△ Less
Submitted 20 April, 2023;
originally announced April 2023.
-
3D hydrodynamics simulations of core convection in supermassive main-sequence stars
Authors:
Simon Blouin,
Huaqing Mao,
Tyrone E. Woods,
Pavel Denissenkov,
Paul R. Woodward,
Falk Herwig
Abstract:
Supermassive stars are Population III stars with masses exceeding $10^4\,M_{\odot}$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code PPMstar to simulate the hydrogen-burning core of a $10^4\,M_{\odot}$ supermassive main-sequence star. The…
▽ More
Supermassive stars are Population III stars with masses exceeding $10^4\,M_{\odot}$ that could be the progenitors of the first supermassive black holes. Their interiors are in a regime where radiation pressure dominates the equation of state. In this work, we use the explicit gas dynamics code PPMstar to simulate the hydrogen-burning core of a $10^4\,M_{\odot}$ supermassive main-sequence star. These are the first 3D hydrodynamics simulations of core convection in supermassive stars. We perform a series of ten simulations at different heating rates and on Cartesian grids with resolutions of $768^3$, $1152^3$ and $1728^3$. We examine different properties of the convective flow, including its large-scale morphology, its velocity spectrum and its mixing properties. We conclude that the radiation pressure-dominated nature of the interior does not noticeably affect the behaviour of convection compared to the case of core convection in a massive main-sequence star where gas pressure dominates. Our simulations also offer support for the use of mixing-length theory in 1D models of supermassive stars.
△ Less
Submitted 17 March, 2023;
originally announced March 2023.
-
3D hydrodynamics simulations of internal gravity waves in red giant branch stars
Authors:
Simon Blouin,
Huaqing Mao,
Falk Herwig,
Pavel Denissenkov,
Paul R. Woodward,
William R. Thompson
Abstract:
We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform sim…
▽ More
We present the first 3D hydrodynamics simulations of the excitation and propagation of internal gravity waves (IGWs) in the radiative interiors of low-mass stars on the red giant branch (RGB). We use the PPMstar explicit gas dynamics code to simulate a portion of the convective envelope and all the radiative zone down to the hydrogen-burning shell of a 1.2$M_{\odot}$ upper RGB star. We perform simulations for different grid resolutions (768$^3$, 1536$^3$ and 2880$^3$), a range of driving luminosities, and two different stratifications (corresponding to the bump luminosity and the tip of the RGB). Our RGB tip simulations can be directly performed at the nominal luminosity, circumventing the need for extrapolations to lower luminosities. A rich, continuous spectrum of IGWs is observed, with a significant amount of total power contained at high wavenumbers. By following the time evolution of a passive dye in the stable layers, we find that IGW mixing in our simulations is weaker than predicted by a simple analytical prescription based on shear mixing and not efficient enough to explain the missing RGB extra mixing. However, we may be underestimating the efficiency of IGW mixing given that our simulations include a limited portion of the convective envelope. Quadrupling its radial extent compared to our fiducial setup increases convective velocities by up to a factor 2 and IGW velocities by up to a factor 4. We also report the formation of a $\sim 0.2\,H_P$ penetration zone and evidence that IGWs are excited by plumes that overshoot into the stable layers.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
3D hydrodynamic simulations of massive main-sequence stars II. Convective excitation and spectra of internal gravity waves
Authors:
William Thompson,
Falk Herwig,
Paul R. Woodward,
Huaqing Mao,
Pavel Denissenkov,
Dominic M. Bowman,
Simon Blouin
Abstract:
Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a 25 $\mathrm{M}_\odot$ star performed with the PPMStar code. The model star has a convective core mass of…
▽ More
Recent photometric observations of massive stars have identified a low-frequency power excess which appears as stochastic low-frequency variability in light curve observations. We present the oscillation properties of high resolution hydrodynamic simulations of a 25 $\mathrm{M}_\odot$ star performed with the PPMStar code. The model star has a convective core mass of $\approx\, 12\, \mathrm{M}_\odot$ and approximately half of the envelope simulated. From this simulation, we extract light curves from several directions, average them over each hemisphere, and process them as if they were real photometric observations. We show how core convection excites waves with a similar frequency as the convective time scale in addition to significant power across a forest of low and high angular degree $l$ modes. We find that the coherence of these modes is relatively low as a result of their stochastic excitation by core convection, with lifetimes on the order of 10s of days. Thanks to the still significant power at higher $l$ and this relatively low coherence, we find that integrating over a hemisphere produces a power spectrum that still contains measurable power up to the Brunt--Väisälä frequency. These power spectra extracted from the stable envelope are qualitatively similar to observations, with same order of magnitude yet lower characteristic frequency. This work further shows the potential of long-duration, high-resolution hydrodynamic simulations for connecting asteroseismic observations to the structure and dynamics of core convection and the convective boundary.
△ Less
Submitted 29 April, 2024; v1 submitted 10 March, 2023;
originally announced March 2023.
-
3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves
Authors:
Falk Herwig,
Paul R. Woodward,
Huaqing Mao,
William R. Thompson,
Pavel Denissenkov,
Josh Lau,
Simon Blouin,
Robert Andrassy,
Adam Paul
Abstract:
We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications.…
▽ More
We performed 3D hydrodynamic simulations of the inner $\approx 50\%$ radial extent of a $25\ \mathrm{M_\odot}$ star in the early phase of the main sequence and investigate core convection and internal gravity waves in the core-envelope boundary region. Simulations for different grid resolutions and driving luminosities establish scaling relations to constrain models of mixing for 1D applications. As in previous works, the turbulent mass entrainment rate extrapolated to nominal heating is unrealistically high ($1.58\times 10^{-4}\ \mathrm{M_\odot/yr}$), which is discussed in terms of the non-equilibrium response of the simulations to the initial stratification. We measure quantitatively the effect of mixing due to internal gravity waves excited by core convection interacting with the boundary in our simulations. The wave power spectral density as a function of frequency and wavelength agrees well with the GYRE eigenmode predictions based on the 1D spherically averaged radial profile. A diffusion coefficient profile that reproduces the spherically averaged abundance distribution evolution is determined for each simulation. Through a combination of eigenmode analysis and scaling relations it is shown that in the $N^2$-peak region, mixing is due to internal gravity waves and follows the scaling relation $D_\mathrm{IGW-hydro} \propto L^{4/3}$ over a $\gtrapprox 2\ \mathrm{dex}$ range of heating factors. Different extrapolations of the mixing efficiency down to nominal heating are discussed. If internal gravity wave mixing is due to thermally-enhanced shear mixing, an upper limit is $D_\mathrm{IGW} \lessapprox 2$ to $3\times 10^4\ \mathrm{cm^2/s}$ at nominal heating in the $N^2$-peak region above the convective core.
△ Less
Submitted 19 July, 2023; v1 submitted 9 March, 2023;
originally announced March 2023.
-
Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
Authors:
R. Andrassy,
J. Higl,
H. Mao,
M. Mocák,
D. G. Vlaykov,
W. D. Arnett,
I. Baraffe,
S. W. Campbell,
T. Constantino,
P. V. F. Edelmann,
T. Goffrey,
T. Guillet,
F. Herwig,
R. Hirschi,
L. Horst,
G. Leidi,
C. Meakin,
J. Pratt,
F. Rizzuti,
F. K. Roepke,
P. Woodward
Abstract:
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o…
▽ More
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.
△ Less
Submitted 26 January, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
3D1D hydro-nucleosynthesis simulations. I. Advective-reactive post-processing method and its application to H ingestion into He-shell flash convection in rapidly accreting white dwarfs
Authors:
David Stephens,
Falk Herwig,
Paul Woodward,
Pavel Denissenkov,
Robert Andrassy,
Huaqing Mao
Abstract:
We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simula…
▽ More
We present two mixing models for post-processing of 3D hydrodynamic simulations applied to convective-reactive i-process nucleosynthesis in a rapidly accreting white dwarf (RAWD) with $\mathrm{[Fe/H]} = -2.6$, in which H is ingested into a convective He shell. A 1D advective two-stream model adopts physically motivated radial and horizontal mixing coefficients constrained by 3D hydrodynamic simulations. A simpler approach uses diffusion coefficients calculated from the same simulations. All 3D simulations include the energy feedback of the $^{12}$C(p,$γ$)$^{13}$N reaction from the H entrainment. Global oscillations of shell H ingestion in two of the RAWD simulations cause bursts of entrainment of H and non-radial hydrodynamic feedback. With the same nuclear network as in the 3D simulations, the 1D advective two-stream model reproduces the rate and location of the H burning within the He shell closely matching the 3D simulation predictions, as well as qualitatively displaying the asymmetry of the $X_{\mathrm{H}}$ profiles between the up- and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the up- and downstream. For example, $^{89}$Kr and $^{90}$Kr with half-lives of 3.18 min and 32.3 s differ by a factor 2-10 in the two streams. In this particular application the diffusion approach provides globally the same abundance distribution as the advective two-stream mixing model. The resulting i-process yields are in excellent agreement with observations of the exemplary CEMP-r/s star CS31062-050.
△ Less
Submitted 17 February, 2021; v1 submitted 29 January, 2020;
originally announced January 2020.
-
Simulating 3-D Stellar Hydrodynamics using PPM and PPB Multifluid Gas Dynamics on CPU and CPU+GPU Nodes
Authors:
Paul R. Woodward,
Pei-Hung Lin,
Huaqing Mao,
Robert Andrassy,
Falk Herwig
Abstract:
The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated…
▽ More
The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary.
△ Less
Submitted 31 October, 2018;
originally announced October 2018.
-
H-He Shell Interactions and Nucleosynthesis in Massive Population III Stars
Authors:
Ondrea Clarkson,
Falk Herwig,
Robert Andrassy,
Marco Pignatari,
Paul Woodward,
Huaqing Mao
Abstract:
We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M…
▽ More
We report on our ongoing investigation into the nucleosynthetic and hydrodynamic nature of mixing at the interface between the H- and He-convection zones in massive Pop III stars. Studying a grid of 26 1D stellar evolution simulations with different mixing assumptions, we find that H-He interactions occur in 23/26 cases. We demonstrate the nucleosynthesis expected in a H-He interaction in an 80M$_\odot$. Finally, we describe our progress in simulating a Pop III double convection zone in the PPMStar hydrodynamics code.
△ Less
Submitted 29 October, 2018;
originally announced October 2018.
-
The i-process yields of rapidly-accreting white dwarfs from multicycle He-shell flash stellar evolution models with mixing parameterizations from 3D hydrodynamics simulations
Authors:
Pavel Denissenkov,
Falk Herwig,
Paul Woodward,
Robert Andrassy,
Marco Pignatari,
Samuel Jones
Abstract:
We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash in…
▽ More
We have modelled the multicycle evolution of rapidly-accreting CO white dwarfs (RAWDs) with stable H burning intermittent with strong He-shell flashes on their surfaces for $0.7\leq M_\mathrm{RAWD}/M_\odot\leq 0.75$ and [Fe/H] ranging from $0$ to $-2.6$. We have also computed the i-process nucleosynthesis yields for these models. The i process occurs when convection driven by the He-shell flash ingests protons from the accreted H-rich surface layer, which results in maximum neutron densities $N_\mathrm{n,max}\approx 10^{13}$-$10^{15}\ \mathrm{cm}^{-3}$. The H-ingestion rate and the convective boundary mixing (CBM) parameter $f_\mathrm{top}$ adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through 3D hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter $f_\mathrm{top}$ have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($η< 10\%$). A new result is that RAWDs with [Fe/H]$< -2$ have $η> 20\%$, therefore their masses may reach the Chandrasekhar limit and they may eventually explode as SNeIa. This result and the good fits of the i-process yields from the metal-poor RAWDs to the observed chemical composition of the CEMP-r/s stars suggest that some of the present-day CEMP-r/s stars could be former distant members of triple systems, orbiting close binary systems with RAWDs that may have later exploded as SNeIa.
△ Less
Submitted 9 July, 2019; v1 submitted 10 September, 2018;
originally announced September 2018.
-
3D hydrodynamic simulations of C ingestion into a convective O shell
Authors:
Robert Andrassy,
Falk Herwig,
Paul Woodward,
Christian Ritter
Abstract:
Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, includi…
▽ More
Interactions between convective shells in evolved massive stars have been linked to supernova impostors, to the production of the odd-Z elements Cl, K, and Sc, and they might also help generate the large-scale asphericities that are known to facilitate shock revival in supernova explosion models. We investigate the process of ingestion of C-shell material into a convective O-burning shell, including the hydrodynamic feedback from the nuclear burning of the ingested material. Our 3D hydrodynamic simulations span almost 3 dex in the total luminosity $L_\mathrm{tot}$. All but one of the simulations reach a quasi-stationary state with the entrainment rate and convective velocity proportional to $L_\mathrm{tot}$ and $L_\mathrm{tot}^{1/3}$, respectively. Carbon burning provides $14\,$--$\,33\%$ of the total luminosity, depending on the set of reactions considered. Equivalent simulations done on $768^3$ and $1152^3$ grids are in excellent quantitative agreement. The flow is dominated by a few large-scale convective cells. An instability leading to large-scale oscillations with Mach numbers in excess of $0.2$ develops in an experimental run with the energy yield from C burning increased by a factor of 10. This run represents most closely the conditions expected in a violent O-C shell merger, which is a potential production site for odd-Z elements such as K and Sc and which may seed asymmetries in the supernova progenitor. 1D simulations may underestimate the energy generation from the burning of ingested material by as much as a factor two owing to their missing the effect of clumpiness of entrained material on the nuclear reaction rate.
△ Less
Submitted 21 October, 2019; v1 submitted 12 August, 2018;
originally announced August 2018.
-
Cyberhubs: Virtual Research Environments for Astronomy
Authors:
Falk Herwig,
Robert Andrassy,
Nic Annau,
Ondrea Clarkson,
Benoit Cote,
Aaron D'Sa,
Sam Jones,
Belaid Moa,
Jericho O'Connell,
David Porter,
Christian Ritter,
Paul Woodward
Abstract:
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to med…
▽ More
Collaborations in astronomy and astrophysics are faced with numerous cyber infrastructure challenges, such as large data sets, the need to combine heterogeneous data sets, and the challenge to effectively collaborate on those large, heterogeneous data sets with significant processing requirements and complex science software tools. The cyberhubs system is an easy-to-deploy package for small to medium-sized collaborations based on the Jupyter and Docker technology, that allows web-browser enabled, remote, interactive analytic access to shared data. It offers an initial step to address these challenges. The features and deployment steps of the system are described, as well as the requirements collection through an account of the different approaches to data structuring, handling and available analytic tools for the NuGrid and PPMstar collaborations. NuGrid is an international collaboration that creates stellar evolution and explosion physics and nucleosynthesis simulation data. The PPMstar collaboration performs large-scale 3D stellar hydrodynamics simulation of interior convection in the late phases of stellar evolution. Examples of science that is presently performed on cyberhubs, in the areas 3D stellar hydrodynamic simulations, stellar evolution and nucleosynthesis and Galactic chemical evolution, are presented.
△ Less
Submitted 6 February, 2018;
originally announced February 2018.
-
Convective-reactive nucleosynthesis of K, Sc, Cl and p-process isotopes in O-C shell mergers
Authors:
C. Ritter,
R. Andrassy,
B. Côté,
F. Herwig,
P. R. Woodward,
M. Pignatari,
S. Jones
Abstract:
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1…
▽ More
We address the deficiency of odd-Z elements P, Cl, K and Sc in galactic chemical evolution models through an investigation of the nucleosynthesis of interacting convective O- and C shells in massive stars. 3D hydrodynamic simulations of O-shell convection with moderate C-ingestion rates show no dramatic deviation from spherical symmetry. We derive a spherically averaged diffusion coefficient for 1D nucleosynthesis simulations which show that such convective-reactive ingestion events can be a production site for P, Cl, K and Sc. An entrainment rate of $10^{-3}{{\rm M}_\odot}$/s features overproduction factors $OP_\mathrm{s} \approx 7$. Full O-C shell mergers in our 1D stellar evolution massive star models have overproduction factors $OP_\mathrm{m}>1 \mathrm{dex}$ but for such cases 3D hydrodynamic simulations suggest deviations from spherical symmetry. p-process species can be produced with overproduction factors of $OP_\mathrm{m}>1 \mathrm{dex}$, e.g. for $^{130,132}Ba$. Using the uncertain prediction of the $15{{\rm M}_\odot}$, $Z=0.02$ massive star model ($OP_\mathrm{m} \approx 15$) as representative for merger or entrainment convective-reactive events involving O- and C-burning shells, and assume that such events occur in more than 50% of all stars, our chemical evolution models reproduce the observed Galactic trends of the odd-Z elements.
△ Less
Submitted 25 September, 2017; v1 submitted 19 April, 2017;
originally announced April 2017.
-
Idealised hydrodynamic simulations of turbulent oxygen-burning shell convection in 4π geometry
Authors:
Samuel Jones,
Robert Andrassy,
Stou Sandalski,
Austin Davis,
Paul Woodward,
Falk Herwig
Abstract:
This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4π$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The…
▽ More
This work investigates the properties of convection in stars with particular emphasis on entrainment across the upper convective boundary (CB). Idealised simulations of turbulent convection in the O-burning shell of a massive star are performed in $4π$ geometry on $768^3$ and $1536^3$ grids, driven by a representative heating rate. A heating series is also performed on the $768^3$ grid. The $1536^3$ simulation exhibits an entrainment rate at the upper CB of $1.33\times10^{-6}~M_\odot~\mathrm{s}^{-1}$. The $768^3$ simulation with the same heating rate agrees within 17 per cent. The entrainment rate at the upper convective boundary is found to scale linearly with the driving luminosity and with the cube of the shear velocity at the upper boundary, while the radial RMS fluid velocity scales with the cube root of the driving luminosity, as expected. The mixing is analysed in a 1D diffusion framework, resulting in a simple model for CB mixing. The analysis confirms previous findings that limiting the MLT mixing length to the distance to the CB in 1D simulations better represents the spherically-averaged radial velocity profiles from the 3D simulations and provides an improved determination of the reference diffusion coefficient $D_0$ for the exponential diffusion CB mixing model in 1D. From the 3D simulation data we adopt as the convective boundary the location of the maximum gradient in the horizontal velocity component which has $2σ$ spatial fluctuations of $\approx0.17 H_P$ . The exponentially decaying diffusion CB mixing model with $f = 0.03$ reproduces the spherically-averaged 3D abundance profiles.
△ Less
Submitted 2 November, 2016; v1 submitted 12 May, 2016;
originally announced May 2016.
-
The i-process and CEMP-r/s stars
Authors:
L. Dardelet,
C. Ritter,
P. Prado,
E. Heringer,
C. Higgs,
S. Sandalski,
S. Jones,
P. Denissenkov,
K. Venn,
M. Bertolli,
M. Pignatari,
P. Woodward,
F. Herwig
Abstract:
We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-s+r (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple s…
▽ More
We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-s+r (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars, super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found that single-exposure one-zone simulations of the i-process reproduce the abundance patterns in some of the CEMP-r/s stars much better than the model that assumes a superposition of yields from s- and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-process revealed that they could have a significant impact on the i-process yields obtained in our idealized one-zone calculations, leading, for example, to ~0.7dex uncertainty in our predicted [Ba/La] ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGB Sakurai's object have discovered a new mode of non-radial instabilities: the Global Oscillation of Shell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulations cannot be used to accurately model physical conditions for the i-process.
△ Less
Submitted 26 July, 2015; v1 submitted 20 May, 2015;
originally announced May 2015.
-
Global non-spherical oscillations in 3-D $4π$ simulations of the H-ingestion flash
Authors:
Falk Herwig,
Paul R. Woodward,
Pei-Hung Lin,
Mike Knox,
Chris Fryer
Abstract:
We performed 3-D simulations of proton-rich material entrainment into \czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurai's object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We inc…
▽ More
We performed 3-D simulations of proton-rich material entrainment into \czw-rich He-shell flash convection and the subsequent H-ingestion flash that took place in the post-AGB star Sakurai's object. Observations of the transient nature and anomalous abundance features are available to validate our method and assumptions, with the aim to apply them to very low metallicity stars in the future. We include nuclear energy feedback from H burning and cover the full $4π$ geometry of the shell. Runs on $768^3$ and $1536^3$ grids agree well with each other and have been followed for $1500\mathrm{min}$ and $1200\mathrm{min}$. After a $850\mathrm{min}$ long quiescent entrainment phase the simulations enter into a global non-spherical oscillation that is launched and sustained by individual ignition events of H-rich fluid pockets. Fast circumferential flows collide at the antipode and cause the formation and localized ignition of the next H-overabundant pocket. The cycle repeats for more than a dozen times while its amplitude decreases. During the global oscillation the entrainment rate increases temporarily by a factor $\approx 100$. Entrained entropy quenches convective motions in the upper layer until the burning of entrained H establishes a separate convection zone. The lower-resolution run hints at the possibility that another global oscillation, perhaps even more violent will follow. The location of the H-burning convection zone agrees with a 1-D model in which the mixing efficiency is calibrated to reproduce the light curve. The simulations have been performed at the NSF Blue Waters supercomputer at NCSA.
△ Less
Submitted 25 September, 2014; v1 submitted 17 October, 2013;
originally announced October 2013.
-
Hydrodynamic simulations of H entrainment at the top of He-shell flash convection
Authors:
Paul R. Woodward,
Falk Herwig,
Pei-Hung Lin
Abstract:
We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4π$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment pro…
▽ More
We present the first 3-dimensional, fully compressible gas-dynamics simulations in $4π$ geometry of He-shell flash convection with proton-rich fuel entrainment at the upper boundary. This work is motivated by the insufficiently understood observed consequences of the H-ingestion flash in post-AGB stars (Sakurai's object) and metal-poor AGB stars. Our investigation is focused on the entrainment process at the top convection boundary and on the subsequent advection of H-rich material into deeper layers, and we therefore ignore the burning of the proton-rich fuel in this study. We find that, for our deep convection zone, coherent convective motions of near global scale appear to dominate the flow. At the top boundary convective shear flows are stable against Kelvin-Helmholtz instabilities. However, such shear instabilities are induced by the boundary-layer separation in large-scale, opposing flows. This links the global nature of thick shell convection with the entrainment process. We establish the quantitative dependence of the entrainment rate on grid resolution. With our numerical technique simulations with $1024^3$ cells or more are required to reach a numerical fidelity appropriate for this problem. However, only the result from the $1536^3$ simulation provides a clear indication that we approach convergence with regard to the entrainment rate. Our results demonstrate that our method, which is described in detail, can provide quantitative results related to entrainment and convective boundary mixing in deep stellar interior environments with veryvstiff convective boundaries. For the representative case we study in detail, we find an entrainment rate of $4.38 \pm 1.48 \times 10^{-13}M_\odot \mathrm{/s}$.
△ Less
Submitted 25 September, 2014; v1 submitted 15 July, 2013;
originally announced July 2013.
-
Convective-reactive proton-C12 combustion in Sakurai's object (V4334 Sagittarii) and implications for the evolution and yields from the first generations of stars
Authors:
Falk Herwig,
Marco Pignatari,
Paul R. Woodward,
David H. Porter,
Gabriel Rockefeller,
Chris L. Fryer,
Michael Bennett,
Raphael Hirschi
Abstract:
Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-…
▽ More
Depending on mass and metallicity as well as evolutionary phase, stars occasionally experience convective-reactive nucleosynthesis episodes. We specifically investigate the situation when nucleosynthetically unprocessed, H-rich material is convectively mixed with a He-burning zone, for example in convectively unstable shell on top of electron-degenerate cores in AGB stars, young white dwarfs or X-ray bursting neutron stars. Such episodes are frequently encountered in stellar evolution models of stars of extremely low or zero metal content [...] We focus on the convective-reactive episode in the very-late thermal pulse star Sakurai's object (V4334 Sagittarii). Asplund etal. (1999) determined the abundances of 28 elements, many of which are highly non-solar, ranging from H, He and Li all the way to Ba and La, plus the C isotopic ratio. Our simulations show that the mixing evolution according to standard, one-dimensional stellar evolution models implies neutron densities in the He that are too low to obtain a significant neutron capture nucleosynthesis on the heavy elements. We have carried out 3D hydrodynamic He-shell flash convection [...] we assume that the ingestion process of H into the He-shell convection zone leads only after some delay time to a sufficient entropy barrier that splits the convection zone [...] we obtain significantly higher neutron densities (~few 10^15 1/cm^3) and reproduce the key observed abundance trends found in Sakurai's object. These include an overproduction of Rb, Sr and Y by about 2 orders of magnitude higher than the overproduction of Ba and La. Such a peculiar nucleosynthesis signature is impossible to obtain with the mixing predictions in our one-dimensional stellar evolution models. [...] We determine how our results depend on uncertainties of nuclear reaction rates, for example for the C13(α, n)O16 reaction.
△ Less
Submitted 4 January, 2011; v1 submitted 10 February, 2010;
originally announced February 2010.
-
The Hydrodynamic Environment for the s Process in the He-Shell Flash of AGB Stars
Authors:
Paul R. Woodward,
David H. Porter,
Falk Herwig,
Marco Pignatari,
Jagan Jayaraj,
Pei-Hung Lin
Abstract:
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodyn…
▽ More
The He-shell flash convection in AGB stars is the site for the high-temperature component of the s-process in low- and intermediate mass giants, driven by the Ne22 neutron source. [...] The upper convection boundary plays a critical role during the H-ingestion episode that may lead to neutron-bursts in the most metal-poor AGB stars. We address these problems through global 3-dimensional hydrodynamic simulations including the entire spherical He-shell flash convection zone (as oposed to the 3D box-in-a-star simulations). An important aspect of our current effort is to establish the feasibility of our appoach. We explain why we favour the explicit treatment over the anelastic approximation for this problem. The simulations presented in this paper use a Cartesian grid of 512^3 cells and have been run on four 8-core workstations for four days to simulate ~5000s, which corresponds to almost ten convective turn-over times. The convection layer extends radially at the simulated point in the flash evolution over 7 H_p pressure scale-heights and exceeds the size of the underlying core. Convection is dominated by large convective cells that fill more than an entire octant. [...]
△ Less
Submitted 11 January, 2009;
originally announced January 2009.
-
Nuclear burning and mixing in the first stars: entrainment at a convective boundary using the PPB advection scheme
Authors:
Paul Woodward,
Falk Herwig,
David Porter,
Tyler Fuchs,
Anthony Nowatzki,
Marco Pignatari
Abstract:
The evolution of the first generations of stars at zero or extremly low metallicity, and especially some crucial properties like the primary N14 production, is charactarized by convective-reactive mixing events that are mostly absent from similar evolution phases at solar-like metallicity. These episodes occur when unprocessed H-rich material is mixed accross a convective boundary into C12 rich…
▽ More
The evolution of the first generations of stars at zero or extremly low metallicity, and especially some crucial properties like the primary N14 production, is charactarized by convective-reactive mixing events that are mostly absent from similar evolution phases at solar-like metallicity. These episodes occur when unprocessed H-rich material is mixed accross a convective boundary into C12 rich He-burning material, as for example in He-shell flashes of extremely-low metallicity AGB stars. In this paper we describe the astrophysical context of such convective-reactive events, including the difficulty of current one-dimensional stellar evolution models to correctly simulate these evolutionary phases. We then describe the requirements and current state of modeling convective-reactive processes in the first stars environment. We demonstrate some of the new concepts that we are applying to this problem, i.e. the highly accurate PPB advection scheme in the framework of PPM hydrodynamic simulations of mixing accross a very stiff convective boundary. We show initial results of such simulations that address the first non-reactive step of this problem, which is the entrainment of H at the top boundary of the He-shell flash convection zone.
△ Less
Submitted 13 November, 2007;
originally announced November 2007.
-
Convective and non-convective mixing in AGB stars
Authors:
Falk Herwig,
Bernd Freytag,
Tyler Fuchs,
James P. Hansen,
Robert M. Hueckstaedt,
David H. Porter,
Francis X. Timmes,
Paul R. Woodward
Abstract:
We review the current state of modeling convective mixing in AGB stars. The focus is on results obtained through multi-dimensional hydrodynamic simulations of AGB convection, both in the envelope and the unstable He-shell. Using two different codes and a wide range of resolutions and modeling assumptions we find that mixing across convective boundaries is significant for He-shell flash convectio…
▽ More
We review the current state of modeling convective mixing in AGB stars. The focus is on results obtained through multi-dimensional hydrodynamic simulations of AGB convection, both in the envelope and the unstable He-shell. Using two different codes and a wide range of resolutions and modeling assumptions we find that mixing across convective boundaries is significant for He-shell flash convection. We present a preliminary quantitative analysis of this convectively induced extra mixing, based on a sub-set of our simulations. Other non-standard mixing will be discussed briefly.
△ Less
Submitted 3 September, 2007;
originally announced September 2007.