skip to main content

Title: The impact of (n,γ) reaction rate uncertainties on the predicted abundances of i-process elements with 32 ≤ Z ≤ 48 in the metal-poor star HD94028

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 more » 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
 ;  ;  ;  ;  
Award ID(s):
1913554 1927130
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 5179-5187
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 ofmore »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.« less
  2. 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
  3. Abstract Despite the advances in discovering new nuclei, modeling microscopic nuclear structure, nuclear reactors, and stellar nucleosynthesis, we still lack a systemic tool, such as a network approach, to understand the structure and dynamics of over 70 thousands reactions compiled in JINA REACLIB. To this end, we develop an analysis framework, under which it is simple to know which reactions generally are possible and which are not, by counting neutrons and protons incoming to and outgoing from any target nucleus. Specifically, we assemble here a nuclear reaction network in which a node represents a nuclide, and a link represents amore »direct reaction between nuclides. Interestingly, the degree distribution of nuclear network exhibits a bimodal distribution that significantly deviates from the common power-law distribution of scale-free networks and Poisson distribution of random networks. Based on the dynamics from the cross section parameterizations in REACLIB, we surprisingly find that the distribution is universal for reactions with a rate below the threshold, λ < e − T γ , where T is the temperature and γ ≈ 1.05. Moreover, we discover three rules that govern the structure pattern of nuclear reaction network: (i) reaction-type is determined by linking choices, (ii) network distances between the reacting nuclides on 2D grid of Z vs N of nuclides are short, and (iii) each node in- and out-degrees are close to each other. By incorporating these three rules, our model universally unveils the underlying nuclear reaction patterns hidden in a large and dense nuclear reaction network regardless of nuclide chart expansions. It enables us to predict missing links that represent possible new nuclear reactions not yet discovered.« less
  4. Abstract. Photoacoustic spectroscopy (PAS) has become a popular technique for measuringabsorption of light by atmospheric aerosols in both the laboratory andfield campaigns. It has low detection limits, measures suspended aerosols,and is insensitive to scattering. But PAS requires rigorous calibration to beapplied quantitatively. Often, a PAS instrument is either filled with a gasof known concentration and absorption cross section, such that the absorptionin the cell can be calculated from the product of the two, or the absorptionis measured independently with a technique such as cavity ring-downspectroscopy. Then, the PAS signal can be regressed upon the known absorptionto determine a calibration slopemore »that reflects the sensitivity constant ofthe cell and microphone. Ozone has been used for calibrating PAS instrumentsdue to its well-known UV–visible absorption spectrum and the ease with whichit can be generated. However, it is known to photodissociate up toapproximately 1120nm via the O3 + hν(&gt;1.1eV)O2(3Σg-) + O(3P) pathway, which is likely tolead to inaccuracies in aerosol measurements. Two recent studies haveinvestigated the use of O3 for PAS calibration but have reachedseemingly contradictory conclusions with one finding that it results in asensitivity that is a factor of 2 low and the other concluding that it isaccurate. The present work is meant to add to this discussion by exploringthe extent to which O3 photodissociates in the PAS cell and the rolethat the identity of the bath gas plays in determining the PAS sensitivity.We find a 5% loss in PAS signal attributable to photodissociation at 532nmin N2 but no loss in a 5% mixture of O2 in N2.Furthermore, we discovered a dramatic increase of more than a factor of 2in the PAS sensitivity as we increased the O2 fraction in the bathgas, which reached an asymptote near 100% O2 that nearly matched thesensitivity measured with both NO2 and nigrosin particles. Weinterpret this dependence with a kinetic model that suggests the reason forthe observed results is a more efficient transfer of energy from excitedO3 to O2 than to N2 by a factor of 22–55 depending onexcitation wavelength. Notably, the two prior studies on this topic useddifferent bath gas compositions, and although the results presented here donot fully resolve the differences in their results, they may at leastpartially explain them.

    « less
  5. 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