Code dependencies of pre-supernova evolution and nucleosynthesis in massive stars: Evolution to the end of core helium burning

Massive stars are key sources of radiative, kinetic and chemical feedback in the Universe. Grids of massive star models computed by different groups each using their own codes, input physics choices and numerical approximations, however, lead to inconsistent results for the same stars. We use three of these 1D codes – genec, kepler and mesa – to compute non-rotating stellar models of 15, 20 and 25 M⊙ and compare their nucleosynthesis. We follow the evolution from the main sequence until the end of core helium burning. The genec and kepler models hold physics assumptions used in large grids of published models. The mesa code was set up to use convective core overshooting such that the CO core masses are consistent with those obtained by genec. For all models, full nucleosynthesis is computed using the NuGrid post-processing tool mppnp. We find that the surface abundances predicted by the models are in reasonable agreement. In the helium core, the standard deviation of the elemental overproduction factors for Fe to Mo is less than 30 per cent – smaller than the impact of the present nuclear physics uncertainties. For our three initial masses, the three stellar evolution codes yield consistent results. Differences in key properties of the models, e.g. helium and CO core masses and the time spent as a red supergiant, are traced back to the treatment of convection and, to a lesser extent, mass loss. The mixing processes in stars remain the key uncertainty in stellar modelling. Better constrained prescriptions are thus necessary to improve the predictive power of stellar evolution models

MESA and NuGrid simulations of classical novae: CO and ONe nova nucleosynthesis

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates (\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ \dot{M} = 10^{-10}\,M_\odot/\mbox{yr}, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)Comment: 19 pages, 23 figures, 2 tables, accepted to publication by MNRA

i-process Nucleosynthesis and Mass Retention Efficiency in He-shell Flash Evolution of Rapidly Accreting White Dwarfs

© 2017. The American Astronomical Society. All rights reserved. Based on stellar evolution simulations, we demonstrate that rapidly accreting white dwarfs (WDs) in close binary systems are an astrophysical site for the intermediate neutron-capture process. During recurrent and very strong He-shell flashes in the stable H-burning accretion regime H-rich material enters the He-shell flash convection zone. 12 C(p, γ) 13 N reactions release enough energy to potentially impact convection, and i process is activated through the 13 C(α, n) 16 O reaction. The H-ingestion flash may not cause a split of the convection zone as it was seen in simulations of He-shell flashes in post-AGB and low-Z asymptotic giant branch (AGB) stars. We estimate that for the production of first-peak heavy elements this site can be of similar importance for galactic chemical evolution as the s-process production by low-mass AGB stars. The He-shell flashes result in the expansion and, ultimately, ejection of the accreted and then i-process enriched material, via super-Eddington-luminosity winds or Roche-lobe overflow. The WD models do not retain any significant amount of the accreted mass, with a He retention efficiency of ≲ 10% depending on mass and convective boundary mixing assumptions. This makes the evolutionary path of such systems to supernova Ia explosion highly unlikely

Stellar neutron capture cross sections of ⁴¹K and ⁴⁵Sc

The neutron capture cross sections of light nuclei (

High-resolution abundance analysis of red giants in the globular cluster NGC 6522

The [Sr/Ba] and [Y/Ba] scatter observed in some galactic halo stars that are very metal-poor stars and in a few individual stars of the oldest known Milky Way globular cluster NGC 6522,have been interpreted as evidence of early enrichment by massive fast-rotating stars (spinstars). Because NGC 6522 is a bulge globular cluster, the suggestion was that not only the very-metal poor halo stars, but also bulge stars at [Fe/H]~-1 could be used as probes of the stellar nucleosynthesis signatures from the earlier generations of massive stars, but at much higher metallicity. For the bulge the suggestions were based on early spectra available for stars in NGC 6522, with a medium resolution of R~22,000 and a moderate signal-to-noise ratio. The main purpose of this study is to re-analyse the NGC 6522 stars previously reported using new high-resolution (R~45,000) and high signal-to-noise spectra (S/N>100). We aim at re-deriving their stellar parameters and elemental ratios, in particular the abundances of the neutron-capture s-process-dominated elements such as Sr, Y, Zr, La, and Ba, and of the r-element Eu. High-resolution spectra of four giants belonging to the bulge globular cluster NGC 6522 were obtained at the 8m VLT UT2-Kueyen telescope with the UVES spectrograph in FLAMES-UVESconfiguration. The spectroscopic parameters were derived based on the excitation and ionization equilibrium of \ion{Fe}{I} and \ion{Fe}{II}. Our analysis confirms a metallicity [Fe/H] = -0.95+-0.15 for NGC 6522, and the overabundance of the studied stars in Eu (with +~0.2 < [Eu/Fe] < +~0.4) and alpha-elements O and Mg. The neutron-capture s-element-dominated Sr, Y, Zr, Ba, La now show less pronounced variations from star to star. Enhancements are in the range 0.0 < [Sr/Fe] < +0.4, +0.23 < [Y/Fe] < +0.43, 0.0 < [Zr/Fe] < +0.4, 0.0 < [La/Fe] < +0.35,and 0.05 < [Ba/Fe] < +0.55.Comment: date of acceptation: 31/07/2014, in press, 24 pages, 19 figures,Astronomy & Astrophysics, 201

Convective overshooting and production of s-nuclei in massive stars during their core He-burning phase

With the "post-processing" technique we explore the role of the convective overshooting on the production of s-nuclei in stellar models of different initial mass and metallicity ($15 \leq M_{ZAMS}/M_{\odot} \leq 25$; $10^{-4} \leq Z \leq 0.02$), considering a range of values for the parameter $f$, which determines the overall efficiency of convective overshooting.We find enhancements in the production of s-nuclei until a factor $\sim 6$ (measured as the average overproduction factor of the 6 s-only nuclear species with $60\lesssim A\lesssim90$) in all our models of different initial mass and metallicity with $f$ in the range $0.01{-}0.035$ (i.e. models with overshooting) compared to the production obtained with "no-overshooting" models (i.e. models with the same initial mass and metallicity, but $f=10^{-5}$). Moreover the results indicate that the link between the overshooting parameter $f$ and the s-process efficiency is essentially monotonic in all our models of different initial mass and metallicity. Also evident is the higher s-process efficiency when we progressively increase for a given f value both the mass of the models from 15 M$_\odot$ to 25 M$_\odot$ and the Z value from 10$^{-4}$ to 0.02. We also briefly discuss the possible consequences of these results for some open questions linked to the s-process weak component efficiency, as well as a "rule of thumb" to evaluate the impact of the convective overshooting on the yields of a generation of stars.Comment: 12 pages, 6 figures, A&A accepted (corrected typos plus minor changes in order to fulfill the guidelines for A&A manuscripts

Stellar (n,γ) cross sections of ²³Na

The cross section of the ²³Na(n,γ)²⁴Na reaction has been measured via the activation method at the Karlsruhe 3.7 MV Van de Graaff accelerator. NaCl samples were exposed to quasistellar neutron spectra at kT = 5.1 and 25 keV produced via the ¹⁸O(p,n)¹⁸F and ⁷Li(p,n)⁷Be reactions, respectively. The derived capture cross sections (σ)kT=5keV = 9.1 ± 0.3mb and (σ)kT=25keV = 2.03 ± 0.05 mb are significantly lower than reported in literature. These results were used to substantially revise the radiative width of the first ²³Na resonance and to establish an improved set of Maxwellian average cross sections. The implications of the lower capture cross section for current models of s-process nucleosynthesis are discussed

Stellar neutron capture cross sections of ²⁰ ²¹ ²²Ne

The stellar (n,γ) cross sections of the Ne isotopes are important for a number of astrophysical quests, i.e., for the interpretation of abundance patterns in presolar material or with respect to the s-process neutron balance in red giant stars. This paper presents resonance studies of experimental data in the keV range, which had not been fully analyzed before. The analyses were carried out with the R-matrix code sammy. With these results for the resonant part and by adding the components due to direct radiative capture, improved Maxwellian-averaged cross sections (MACS) could be determined. At kT=30keV thermal energy we obtain MACS values of 240±29,1263±160, and 53.2±2.7 μbarn for ²⁰Ne,²¹Ne, and ²²Ne, respectively. In earlier work the stellar rates of ²⁰Ne and ²¹Ne had been grossly overestimated. ²²Ne and ²⁰Ne are significant neutron poisons for the s process in stars because their very small MACS values are compensated by their large abundances