skip to main content


Title: The i-process yields of rapidly accreting white dwarfs from multicycle He-shell flash stellar evolution models with mixing parametrizations from 3D hydrodynamics simulations
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 ≤ MRAWD/M⊙ ≤ 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 Nn, max ≈ 1013–1015 cm−3. The H-ingestion rate and the convective boundary mixing (CBM) parameter ftop adopted in the one-dimensional nucleosynthesis and stellar evolution models are constrained through three-dimensional (3D) hydrodynamic simulations. The mass ingestion rate and, for the first time, the scaling laws for the CBM parameter ftop have been determined from 3D hydrodynamic simulations. We confirm our previous result that the high-metallicity RAWDs have a low mass retention efficiency ($\eta \lesssim 10{{\ \rm per\ cent}}$). A new result is that RAWDs with [Fe/H] $\lesssim -2$ have $\eta \gtrsim 20{{\ \rm per\ cent}}$; 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.  more » « less
Award ID(s):
1814181
NSF-PAR ID:
10166285
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
488
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
4258 to 4270
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    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 [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 12C(p, γ)13N 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 XH profiles between the upstream and downstream. With a full i-process network the advective mixing model captures the difference in the n-capture nucleosynthesis in the upstream and downstream. For example, 89Kr and 90Kr with half-lives of $3.18\,\,\mathrm{\mathrm{min}}$ and $32.3\,\,\mathrm{\mathrm{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. 
    more » « less
  2. null (Ed.)
    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. 
    more » « less
  3. null (Ed.)
    ABSTRACT Carbon enhanced metal poor (CEMP)-no stars, a subset of CEMP stars ($\rm [C/Fe]\ge 0.7$ and $\rm [Fe/H]\lesssim -1$) have been discovered in ultra-faint dwarf (UFD) galaxies, with $M_{\rm vir}\approx 10^8{\, \mathrm{ M}_\odot }$ and $M_{\ast }\approx 10^3-10^4{\, \mathrm{ M}_\odot }$ at z = 0, as well as in the halo of the Milky Way (MW). These CEMP-no stars are local fossils that may reflect the properties of the first (Pop III) and second (Pop II) generation of stars. However, cosmological simulations have struggled to reproduce the observed level of carbon enhancement of the known CEMP-no stars. Here, we present new cosmological hydrodynamic zoom-in simulations of isolated UFDs that achieve a gas mass resolution of $m_{\rm gas}\approx 60{\, \mathrm{ M}_\odot }$. We include enrichment from Pop III faint supernovae (SNe), with ESN = 0.6 × 1051 erg, to understand the origin of CEMP-no stars. We confirm that Pop III and Pop II stars are mainly responsible for the formation of CEMP and C-normal stars, respectively. New to this study, we find that a majority of CEMP-no stars in the observed UFDs and the MW halo can be explained by Pop III SNe with normal explosion energy (ESN = 1.2 × 1051 erg) and Pop II enrichment, but faint SNe might also be needed to produce CEMP-no stars with $\rm [C/Fe]\gtrsim 2$, corresponding to the absolute carbon abundance of $\rm A(C)\gtrsim 6.0$. Furthermore, we find that while we create CEMP-no stars with high carbon ratio $\rm [C/Fe]\approx 3-4$, by adopting faint SNe, it is still challenging to reproduce CEMP-no stars with extreme level of carbon abundance of $\rm A(C)\approx 7.0-7.5$, observed both in the MW halo and UFDs. 
    more » « less
  4. ABSTRACT

    We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $\lesssim 1~{{\ \rm per\ cent}}$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers.

     
    more » « less
  5. ABSTRACT

    We present detailed chemical compositions of four stars on the first-ascent red giant branch that are classified as chemically peculiar, but lack comprehensive analyses at high spectral resolution. For BD+03°2688, HE 0457−1805, HE 1255−2324, and HE 2207−1746, we derived metallicities [Fe/H] = −1.21, −0.19, −0.31, and −0.55, respectively, indicating a range in Galactic population membership. In addition to atmospheric parameters, we extracted elemental abundances for 28 elements, including the evolutionary-sensitive CNO group and 12C/13C ratios. Novel results are also presented for the heavy elements tungsten and thallium. All four stars have very large enhancements of neutron-capture elements, with high [La/Eu] ratios indicating enrichments from the slow neutron capture (s-process). To interpret these abundances, all indicative of [s/Fe] >1.0, we compared our results with data from literature, as well as with predictions from the Monash and fruity s-process nucleosynthesis models. BD+03°2688, HE 1255−2324, and HE 2207−1746 show C/O >1, while HE 0457−1805 has C/O <1. Since HE 0457−1805 and HE 1255−2324 are binary stars, their peculiarities are attributable to mass transfer. We identified HE 0457−1805 as a new barium giant star, and HE 1255−2324 as a new CH star, in fact a higher metallicity analogue CEMP-r/s star; the single object reported in literature so far with similar characteristics is the barium star HD 100503 ([Fe/H] = −0.72). A systematic monitoring is needed to confirm the binary nature of BD+03°2688 and HE 2207−1746, which are probably CH stars.

     
    more » « less