Search for: All records

A few per cent of red giants are enriched in lithium with $A(\mathrm{Li}) \gt 1.5$. Their evolutionary status has remained uncertain because these Li-rich giants can 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. 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. We propose that the RGB stars are getting enriched in Li by the RGB extra mixing that is getting enhanced and begins to produce Li, instead of destroying it, when the stars are approaching the RGB tip. After a discussion of several mechanisms of the RGB extra mixing, including the joint operation of rotation-driven meridional circulation and turbulent diffusion, the azimuthal magnetorotational instability (AMRI), thermohaline convection, buoyancy of magnetic flux tubes, and internal gravity waves, and based on results of (magneto-) hydrodynamics simulations and asteroseismology observations, we are inclined to conclude that it is the mechanism of the AMRI or magnetically enhanced thermohaline convection, that is most likely to support our hypothesis.

 

We present 3D hydrodynamical simulations of core convection with a stably stratified envelope of a 25M_⊙star in the early phase of the main sequence. We use the explicit gas-dynamics codePPMstar, 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. We perform very long simulations on 896³grids accelerated by luminosity boost factors of 1000, 3162 and 10,000. 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 1D method to predict the extent and character of penetrative convection beyond the Schwarzschild boundary. The 1D model is based on a spherically averaged reduced entropy equation that takes the turbulent dissipation as input from the 3D hydrodynamic simulation and takes buoyancy and all other energy sources and sinks into account. This 1D 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.

 

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\,\,{\rm{M}_\odot }$ star performed with the PPMstar code. The model star has a convective core mass of $\approx 12\,\,{\rm{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 of 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 the 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.

 

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$\, {\rm M}_{\odot }$ of a 3 M⊙ CHeB star. Simulations are performed using different Cartesian grid resolutions (7683, 11523, and 17283) and heating rates. We use two different initial states, one based on mesas'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.

 

We performed 3D hydrodynamic simulations of the inner $\approx 50{{\ \rm per\ cent}}$ radial extent of a $25\,\,\mathrm{\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{\mathrm{M}_\odot \, {\mathrm{yr}}^{-1}}$), 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 N2-peak region, mixing is due to internal gravity waves and follows the scaling relation DIGW-hydro ∝ L4/3 over a $\gtrapprox 2\,\,\mathrm{\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^{-1}}$ at nominal heating in the N2-peak region above the convective core.

 

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\, {\rm M}_{\odot }$ upper RGB star. We perform simulations for different grid resolutions (7683, 15363, and 28803), 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 set-up 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.

 

Supermassive stars are Population III stars with masses exceeding $10^4\, {\rm 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\, {\rm M}_{\odot }$ supermassive main-sequence star. These are the first three-dimensional hydrodynamics simulations of core convection in supermassive stars. We perform a series of 10 simulations at different heating rates and on Cartesian grids with resolutions of 7683, 11523, and 17283. 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 one-dimensional models of supermassive stars.

 

Asymptotic Giant Branch (AGB) stars play a key role in the chemical evolution of galaxies. These stars are the fundamental stellar site for the production of light elements such as C, N and F, and half of the elements heavier than Fe via the slow neutron capture process (s-process). Hence, detailed computational models of AGB stars’ evolution and nucleosynthesis are essential for galactic chemical evolution. In this work, we discuss the progress in updating the NuGrid data set of AGB stellar models and abundance yields. All stellar models have been computed using the MESA stellar evolution code, coupled with the post-processing mppnp code to calculate the full nucleosynthesis. The final data set will include the initial masses Mini/M⊙ = 1, 1.65, 2, 3, 4, 5, 6 and 7 for initial metallicities Z = 0.0001, 0.001, 0.006, 0.01, 0.02 and 0.03. Observed s-process abundances on the surfaces of evolved stars as well as the typical light elements in the composition of H-deficient post-AGB stars are reproduced. A key short-term goal is to complete and expand the AGB stars data set for the full metallicity range. Chemical yield tables are provided for the available models. 

ABSTRACT The abundances of neutron (n)-capture elements in the carbon-enhanced metal-poor (CEMP)-r/s stars agree with predictions of intermediate n-density nucleosynthesis, at Nn ∼ 1013–1015 cm−3, in rapidly accreting white dwarfs (RAWDs). We have performed Monte Carlo simulations of this intermediate-process (i-process) nucleosynthesis to determine the impact of (n,γ) reaction rate uncertainties of 164 unstable isotopes, from 131I to 189Hf, on the predicted abundances of 18 elements from Ba to W. The impact study is based on two representative one-zone models with constant values of Nn = 3.16 × 1014 and 3.16 × 1013 cm−3 and on a multizone model based on a realistic stellar evolution simulation of He-shell convection entraining H in a RAWD model with [Fe/H] = −2.6. For each of the selected elements, we have identified up to two (n,γ) reactions having the strongest correlations between their rate variations constrained by Hauser–Feshbach computations and the predicted abundances, with the Pearson product–moment correlation coefficients |rP| > 0.15. We find that the discrepancies between the predicted and observed abundances of Ba and Pr in the CEMP-i star CS 31062−050 are significantly diminished if the rate of 137Cs(n,γ)138Cs is reduced and the rates of 141Ba(n,γ)142Ba or 141La(n,γ)142La increased. The uncertainties of temperature-dependent β-decay rates of the same unstable isotopes have a negligible effect on the predicted abundances. One-zone Monte Carlo simulations can be used instead of computationally time-consuming multizone Monte Carlo simulations in reaction rate uncertainty studies if they use comparable values of Nn. We discuss the key challenges that RAWD simulations of i process for CEMP-i stars meet by contrasting them with recently published low-Z asymptotic giant branch (AGB) i process.