skip to main content


Title: The impact of (n,γ) reaction rate uncertainties of unstable isotopes on the i -process nucleosynthesis of the elements from Ba to W
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
Award ID(s):
1927130 1913554
NSF-PAR ID:
10296881
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
503
Issue:
3
ISSN:
0035-8711
Page Range / eLocation ID:
3913 to 3925
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Several anomalous elemental abundance ratios have been observed in the metal-poor star HD94028. We assume that its high [As/Ge] ratio is a product of a weak intermediate (i) neutron-capture process. Given that observational errors are usually smaller than predicted nuclear physics uncertainties, we have first set-up a benchmark one-zone i-process nucleosynthesis simulation results of which provide the best fit to the observed abundances. We have then performed Monte Carlo simulations in which 113 relevant (n,γ) reaction rates of unstable species were randomly varied within Hauser–Feshbach model uncertainty ranges for each reaction to estimate the impact on the predicted stellar abundances. One of the interesting results of these simulations is a double-peaked distribution of the As abundance, which is caused by the variation of the 75Ga (n,γ) cross-section. This variation strongly anticorrelates with the predicted As abundance, confirming the necessity for improved theoretical or experimental bounds on this cross-section. The 66Ni (n,γ) reaction is found to behave as a major bottleneck for the i-process nucleosynthesis. Our analysis finds the Pearson product–moment correlation coefficient rP > 0.2 for all of the i-process elements with 32 ≤ Z ≤ 42, with significant changes in their predicted abundances showing up when the rate of this reaction is reduced to its theoretically constrained lower bound. Our results are applicable to any other stellar nucleosynthesis site with the similar i-process conditions, such as Sakurai’s object (V4334 Sagittarii) or rapidly accreting white dwarfs.

     
    more » « less
  2. 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
  3. Abstract

    Accurate nuclear reaction rates for26P(p,γ)27S are pivotal for a comprehensive understanding of therp-process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of26P(p,γ)27S because of the lack of nuclear mass and energy level structure information for27S. We reevaluate this reaction rate using the experimentally constrained27S mass, together with the shell model predicted level structure. It is found that the26P(p,γ)27S reaction rate is dominated by a direct capture reaction mechanism despite the presence of three resonances atE= 1.104, 1.597, and 1.777 MeV above the proton threshold in27S. The new rate is overall smaller than the other previous rates from the Hauser–Feshbach statistical model by at least 1 order of magnitude in the temperature range of X-ray burst interest. In addition, we consistently update the photodisintegration rate using the new27S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in therp-process is explored by postprocessing nucleosynthesis calculations. The final abundance ratio of27S/26P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show that the reaction path26P(p,γ)27S(β+,ν)27P is not as important as previously thought for producing27P. The adoption of the new reaction rates for26P(p,γ)27S only reduces the final production of aluminum by 7.1% and has no discernible impact on the yield of other elements.

     
    more » « less
  4. 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
  5. Abstract

    We present results from high-resolution (R∼ 40,000) spectroscopic observations of over 200 metal-poor stars, mostly selected from the RAVE survey, using the Southern African Large Telescope. We were able to derive stellar parameters for a total of 108 stars; an additional sample of 50 stars from this same effort was previously reported on by Rasmussen et al. Among our newly reported observations, we identify 84 very metal-poor (VMP; [Fe/H] < −2.0, 53 newly identified) stars and three extremely metal-poor (EMP; [Fe/H] < −3.0, one newly identified) stars. The elemental abundances were measured for carbon, as well as several otherα-elements (Mg, Ca, Sc, and Ti), iron-peak elements (Mn, Co, Ni, and Zn), and neutron-capture elements (Sr, Ba, and Eu). Based on these measurements, the stars are classified by their carbon and neutron-capture abundances into carbon-enhanced metal-poor (CEMP; [C/Fe] > +0.70), CEMP subclasses, and by the level of theirr-process abundances. A total of 17 are classified as CEMP stars. There are 11 CEMP-rstars (eight newly identified), one CEMP-sstar (newly identified), two possible CEMP-istars (one newly identified), and three CEMP-no stars (all newly identified) in this work. We found 11 stars (eight newly identified) that are strongly enhanced inr-process elements (r-II; [Eu/Fe] > +0.70), 38 stars (31 newly identified) that are moderately enhanced inr-process elements (r-I; +0.30 < [Eu/Fe] ≤ + 0.70), and one newly identified limited-rstar.

     
    more » « less