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 more » 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. « less
; ; ;
Award ID(s):
1927130 1913554
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
3913 to 3925
Sponsoring Org:
National Science Foundation
More Like this

    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 stellarmore »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.

    « 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 aremore »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.« less
  3. 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 ofmore »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.« less
  4. Abstract We demonstrate that using up to seven stellar abundance ratios can place observational constraints on the star formation histories (SFHs) of Local Group dSphs, using Sculptor dSph as a test case. We use a one-zone chemical evolution model to fit the overall abundance patterns of α elements (which probe the core-collapse supernovae that occur shortly after star formation), s -process elements (which probe AGB nucleosynthesis at intermediate delay times), and iron-peak elements (which probe delayed Type Ia supernovae). Our best-fit model indicates that Sculptor dSph has an ancient SFH, consistent with previous estimates from deep photometry. However, we derivemore »a total star formation duration of ∼0.9 Gyr, which is shorter than photometrically derived SFHs. We explore the effect of various model assumptions on our measurement and find that modifications to these assumptions still produce relatively short SFHs of duration ≲1.4 Gyr. Our model is also able to compare sets of predicted nucleosynthetic yields for supernovae and AGB stars, and can provide insight into the nucleosynthesis of individual elements in Sculptor dSph. We find that observed [Mn/Fe] and [Ni/Fe] trends are most consistent with sub- M Ch Type Ia supernova models, and that a combination of “prompt” (delay times similar to core-collapse supernovae) and “delayed” (minimum delay times ≳50 Myr) r -process events may be required to reproduce observed [Ba/Mg] and [Eu/Mg] trends.« less
  5. Context. Young open clusters (ages of less than 200 Myr) have been observed to exhibit several peculiarities in their chemical compositions. These anomalies include a slightly sub-solar iron content, super-solar abundances of some atomic species (e.g. ionised chromium), and atypical enhancements of [Ba/Fe], with values up to ~0.7 dex. Regarding the behaviour of the other s -process elements like yttrium, zirconium, lanthanum, and cerium, there is general disagreement in the literature: some authors claim that they follow the same trend as barium, while others find solar abundances at all ages. Aims. In this work we expand upon our previous analysismore »of a sample of five young open clusters (IC 2391, IC 2602, IC 4665, NGC 2516, and NGC 2547) and one star-forming region (NGC 2264), with the aim of determining abundances of different neutron-capture elements, mainly Cu  I , Sr  I , Sr  II , Y  II , Zr  II , Ba  II , La  II , and Ce  II . For NGC 2264 and NGC 2547 we present the measurements of these elements for the first time. Methods. We analysed high-resolution, high signal-to-noise spectra of 23 solar-type stars observed within the Gaia -ESO survey. After a careful selection, we derived abundances of isolated and clean lines via spectral synthesis computations and in a strictly differential way with respect to the Sun. Results. We find that our clusters have solar [Cu/Fe] within the uncertainties, while we confirm that [Ba/Fe] is super-solar, with values ranging from +0.22 to +0.64 dex. Our analysis also points to a mild enhancement of Y, with [Y/Fe] ratios covering values between 0 and +0.3 dex. For the other s -process elements we find that [X/Fe] ratios are solar at all ages. Conclusions. It is not possible to reconcile the anomalous behaviour of Ba and Y at young ages with standard stellar yields and Galactic chemical evolution model predictions. We explore different possible scenarios related to the behaviour of spectral lines, from the dependence on the different ionisation stages and the sensitivity to the presence of magnetic fields (through the Landé factor) to the first ionisation potential effect. We also investigate the possibility that they may arise from alterations of the structure of the stellar photosphere due to the increased levels of stellar activity that affect the spectral line formation, and consequently the derived abundances. These effects seem to be stronger in stars at ages of less than ~ 100 Myr. However, we are still unable to explain these enhancements, and the Ba puzzle remains unsolved. With the present study we suggest that other elements, for example Sr, Zr, La, and Ce, might be more reliable tracer of the s -process at young ages, and we strongly encourage further critical observations.« less