Attention:The NSF Public Access Repository (PAR) system and access will be unavailable from 10:00 PM ET on Thursday, July 16 until 12:00 AM ET on Friday 17 due to maintenance. We apologize for the inconvenience.


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
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 Presolar grains are stardust particles that condensed in the ejecta or in the outflows of dying stars and can today be extracted from meteorites. They recorded the nucleosynthetic fingerprint of their parent stars and thus serve as valuable probes of these astrophysical sites. The most common types of presolar silicon carbide grains (called mainstream SiC grains) condensed in the outflows of asymptotic giant branch stars. Their measured silicon isotopic abundances are not significantly influenced by nucleosynthesis within the parent star but rather represent the pristine stellar composition. Silicon isotopes can thus be used as a proxy for galactic chemical evolution (GCE). However, the measured correlation of29Si/28Si versus30Si/28Si does not agree with any current chemical evolution model. Here, we use a Monte Carlo model to vary nuclear reaction rates within their theoretical or experimental uncertainties and process them through stellar nucleosynthesis and GCE models to study the variation of silicon isotope abundances based on these nuclear reaction rate uncertainties. We find that these uncertainties can indeed be responsible for the discrepancy between measurements and models and that the slope of the silicon isotope correlation line measured in mainstream SiC grains agrees with chemical evolution models within the nuclear reaction rate uncertainties. Our result highlights the importance of future precision reaction rate measurements for resolving the apparent data–model discrepancy. 
    more » « less
  4. Context. Barium (Ba) stars are characterised by an abundance of heavy elements made by the slow neutron capture process ( s -process). This peculiar observed signature is due to the mass transfer from a stellar companion, bound in a binary stellar system, to the Ba star observed today. The signature is created when the stellar companion is an asymptotic giant branch (AGB) star. Aims. We aim to analyse the abundance pattern of 169 Ba stars using machine learning techniques and the AGB final surface abundances predicted by the F RUITY and Monash stellar models. Methods. We developed machine learning algorithms that use the abundance pattern of Ba stars as input to classify the initial mass and metallicity of each Ba star’s companion star using stellar model predictions. We used two algorithms. The first exploits neural networks to recognise patterns, and the second is a nearest-neighbour algorithm that focuses on finding the AGB model that predicts the final surface abundances closest to the observed Ba star values. In the second algorithm, we included the error bars and observational uncertainties in order to find the best-fit model. The classification process was based on the abundances of Fe, Rb, Sr, Zr, Ru, Nd, Ce, Sm, and Eu. We selected these elements by systematically removing s -process elements from our AGB model abundance distributions and identifying the elements whose removal had the biggest positive effect on the classification. We excluded Nb, Y, Mo, and La. Our final classification combined the output of both algorithms to identify an initial mass and metallicity range for each Ba star companion. Results. With our analysis tools, we identified the main properties for 166 of the 169 Ba stars in the stellar sample. The classifications based on both stellar sets of AGB final abundances show similar distributions, with an average initial mass of M = 2.23 M ⊙ and 2.34 M ⊙ and an average [Fe/H] = −0.21 and −0.11, respectively. We investigated why the removal of Nb, Y, Mo, and La improves our classification and identified 43 stars for which the exclusion had the biggest effect. We found that these stars have statistically significant and different abundances for these elements compared to the other Ba stars in our sample. We discuss the possible reasons for these differences in the abundance patterns. 
    more » « less
  5. Abstract Accurate nuclear reaction rates for 26 P( p , γ ) 27 S are pivotal for a comprehensive understanding of the rp -process nucleosynthesis path in the region of proton-rich sulfur and phosphorus isotopes. However, large uncertainties still exist in the current rate of 26 P( p , γ ) 27 S because of the lack of nuclear mass and energy level structure information for 27 S. We reevaluate this reaction rate using the experimentally constrained 27 S mass, together with the shell model predicted level structure. It is found that the 26 P( p , γ ) 27 S reaction rate is dominated by a direct capture reaction mechanism despite the presence of three resonances at E = 1.104, 1.597, and 1.777 MeV above the proton threshold in 27 S. 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 new 27 S mass. The influence of new rates of forward and reverse reaction in the abundances of isotopes produced in the rp -process is explored by postprocessing nucleosynthesis calculations. The final abundance ratio of 27 S/ 26 P obtained using the new rates is only 10% of that from the old rate. The abundance flow calculations show that the reaction path 26 P( p , γ ) 27 S( β + , ν ) 27 P is not as important as previously thought for producing 27 P. The adoption of the new reaction rates for 26 P( p , γ ) 27 S only reduces the final production of aluminum by 7.1% and has no discernible impact on the yield of other elements. 
    more » « less