© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. We investigate whether the anomalous elemental abundance patterns in some of the C-enhanced metal-poor-r/s (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from the i-process, a neutron-capture regime at neutron densities intermediate between those typical for the slow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met at multiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars, super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found that single-exposure one-zone simulations of the i-process reproduce the abundance patterns in some of the CEMP-r/s stars much better than the model that assumes a superposition of yields from s and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-process revealed that they could have a significant impact on the i-process yields obtained in our idealized one-zone calculations, leading, for example, to ∼ 0:7dex uncertainty in our predicted [Ba/La] ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGB Sakurai's object have discovered a new mode of non-radial instabilities: the Global Oscillation of Shell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulations cannot be used to accurately model physical conditions for the i-process

Bertolli, M.

Bertolli, Michael

Dardelet, L.

Dardelet, Laurent

Denisenkov, P A

Denissenkov, P.

Heringer, E.

Herwig, F.

Higgs, C.

Higgs, Clare

Jones, S.

Jones, Samuel

Masikiv Heringer, Epson Thiago

Pignatari, M.

Prado, P.

Prado, Pablo

Ritter, C.

Ritter, Christian

Sandalski, S.

Sandalski, Stou

Venn, K.

Venn, Kim

Woodward, P.

English

Repository@Hull - Worktribe

The i-process and CEMP-r/s starsL. Dardeleta, C. Ritterb,c,d , P. Pradob, E. Heringerb, C. Higgsb, S. Sandalskih,c,S. Jonesb,d , P. Denissenkovb,c,d , K.Vennb, M. Bertollie, f ,d , M. Pignatarig,d , P.Woodwardh,c, F. Herwigb,c,daDepartment of Physics, École Normale Supérieure, 45 rue d’Ulm, 75005 Paris, FrancebDepartment of Physics and Astronomy, University of Victoria, Victoria, BC, V8P5C2, CanadacJoint Institute for Nuclear Astrophysics, Center for the Evolution of the Elements, MichiganState University, 640 South Shaw Lane, East Lansing, MI 48824, USAdNuGrid collaboration, http://www.nugridstars.orgeOak Ridge National Laboratory, P.O. Box 2008,Oak Ridge, TN 37831, USAf Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996, USAgDepartment of Physics, University of Basel, Klingelbergstrasse 82, CH-4056 Basel, SwitzerlandhLCSE and Department of Astronomy, University of Minnesota, Minneapolis, MN 55455, USAE-mail: critter@uvic.caWe investigate whether the anomalous elemental abundance patterns in some of the C-enhancedmetal-poor-r/s (CEMP-r/s) stars are consistent with predictions of nucleosynthesis yields from thei-process, a neutron-capture regime at neutron densities intermediate between those typical for theslow (s) and rapid (r) processes. Conditions necessary for the i-process are expected to be met atmultiple stellar sites, such as the He-core and He-shell flashes in low-metallicity low-mass stars,super-AGB and post-AGB stars, as well as low-metallicity massive stars. We have found thatsingle-exposure one-zone simulations of the i-process reproduce the abundance patterns in someof the CEMP-r/s stars much better than the model that assumes a superposition of yields from s-and r-process sources. Our previous study of nuclear data uncertainties relevant to the i-processrevealed that they could have a significant impact on the i-process yields obtained in our idealizedone-zone calculations, leading, for example, to ∼ 0.7dex uncertainty in our predicted [Ba/La]ratio. Recent 3D hydrodynamic simulations of convection driven by a He-shell flash in post-AGBSakurai’s object have discovered a new mode of non-radial instabilities: the Global Oscillation ofShell H-ingestion. This has demonstrated that spherically symmetric stellar evolution simulationscannot be used to accurately model physical conditions for the i-process.XIII Nuclei in the Cosmos,7-11 July, 2014Debrecen, Hungaryc© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. http://pos.sissa.it/arXiv:1505.05500v2  [astro-ph.SR]  26 Jul 2015The i-process and CEMP-r/s stars1. IntroductionThe fraction of C-Enhanced Metal Poor (CEMP) stars (objects with [C/Fe]>1) increases with adecrease in metallicity, reaching 20% at [Fe/H]=-2 [1, 3]. CEMP-r/s stars form a subclass char-acterized by simultaneous enhancements of both s- (e.g., Ba) and r-process (e.g., Eu) elements.[17, 21, 22]. The prevailing model of a CEMP-r/s star assumes that it is a component of a binarysystem initially pre-enriched with r-process material that has additionally accreted s-process mate-rial from its heavier AGB-star companion, which leads to a superposition of both processes [17].However, this model fails to explain CEMP-r/s stars with large [hs/ls]1 and low [La/Eu] ratios. Anexample of such a star is CS 31062-050 shown in Fig. 26 in [5].The intermediate neutron-capture process (hereafter, the i-process, see [7]) operates at neutrondensities between those characteristic of the s- and r-process, when H is convectively entrainedand advected into a He-burning zone. When the advected fluid parcel has reached a depth withT ≈ 1.5×108K, 13N is produced in the reaction 12C(p,γ)13N. The beta-decay of 13N to 13C followedby the reaction 13C(α ,n)16O result in a production of neutrons with a density Nn ≈ 1015cm−3.Possible astrophysical sites for the i-process have been identified in stellar evolution computations.These are the He-shell thermal pulses in the low- and zero-metallicity AGB stars [9], the very latethermal pulse (VLTP) in post-AGB stars (e.g., Sakurai’s object, [13]) and the He-core flash in low-metallicity low-mass stars ([6]). The i-process in VLTP models with mixing assumptions motivatedby 3D hydrodynamic simulations was found to reproduce the observed abundance distribution inSakurai’s object [13]. It may also explain the excess of 32S and anomalous Ti isotopic ratios inpresolar A+B grains [10, 15].In this paper, we investigate if the i-process nucleosynthesis can explain the observed abundancepatterns in peculiar CEMP-r/s stars, such as CS 31062-050. We also briefly discuss the impactof nuclear physics uncertainties on the calculated i-process yields, as well as the results of 3Dhydrodynamic simulations of nuclear burning and mixing under conditions relevant to the i-processin Sakurai’s object recently reported in [14, 24].2. A simple i-proces model for the CEMP-r/s starsWe compare the abundance pattern in CS 31062-050 with results of our nucleosynthesis simu-lations of the i-process in one zone with constant temperature and density corresponding to thephysical conditions in a He pulse-driven convective zone (PDCZ) from 1D stellar evolution modelsof AGB stars. Our assumed values of T = 2× 108 K and ρ = 104 gcm−3 prevent destruction of13N via the (p,γ) channel, but they allow the neutron release via 13C(α,n). In the He PDCZ ofa real star, with H entrained by convection from the surrounding H-rich envelope, the reactions12C(p,γ)13N and 13C(α ,n)16O are spatially separated, the first occurring close to the upper and thesecond close to the lower convective boundary. 13N produced in the first reaction decays to 13Cwhile being transported down by convection. The one-zone model reproduces the characteristici-process neutron density (Nn ∼ 1015cm−3) and neutron exposures of 10 . . . 50 mbarn−1. The initialabundances are taken from the solar abundance distribution [11], with isotopic ratios from [20],1hs and ls are the total abundances of selected elements representing the second (high) and first (light) s-processpeaks in the solar abundance distribution.2The i-process and CEMP-r/s starsFigure 1: The evolution of the neutron density Nn with time for the initial H mass fractions 0.2, 0.1 and 0.05(upper panel). The best fit of the abundance pattern in the CEMP-r/s star CS 31062-050 with the i-processmodel (lower panel, with observational data from [16]). Shown in the bottom-right corner are the fittedvalues of the neutron density Nn and exposure (tau).scaled down to the metallicity Z = 10−3. The abundances of C and O are additionally modified,so that they are close to those found in the He PDCZ, i.e. X(12C) = 0.5 and X(16O) = 0.05. Aspecified fraction of H, assumed to be ingested from the H-rich envelope, is added to the mixture,with the combined mass fraction of C+H being always set to 0.7. The i-process nucleosynthesis isfollowed with the NuGrid single-zone PPN code [12].During the first second, 12C is almost all transformed into 13N which then decays into 13C onthe timescale of 9.97min. After the decay of 13N, the reaction 13C(α ,n)16O releases neutronswith a high neutron density. The upper panel of Fig. 1 shows the evolution of the neutron densityNn with time for different initial mass fractions of H. As the simulation proceeds, the neutronexposure τ =∫ t0 NnvT dt′ (vT is the neutron thermal velocity) increases with time, and the abundance3The i-process and CEMP-r/s starsdistribution is shifted to heavier elements. The model results shown in the lower panel of Fig. 1give the best representation of the abundance pattern in CS 31062-050. Its corresponding neutrondensity and exposure are indicated in the panel’s bottom-right corner.Hydrodynamic simulations suggest that the i-process nucleosynthesis in the stellar environmentis terminated when the energy released through the combustion of ingested H leads to a signifi-cant perturbation of the spherically symmetric convective shell structure. This self-quenching willdepend on the exact configuration of the convective shells at the time of H-ingestion and can beexpected to vary. Therefore, a range of the termination time is expected for the i-process, whichmotivates us to consider it as a free parameter in this simple model. However, in spite of its simplic-ity, our one-zone i-process model reproduces surprisingly well the entire heavy-element abundancedistribution from Y to Ir within the observational errors in CS 31062-050 (Fig. 1). Specifically, thisincludes the large ratio of [hs/ls]≈ 1, the [La/Eu] ratios between 0.0 and 0.5, which are between theratios predicted for the s- and r-processes, as well as high abundances in the Er-W region comparedto those in the Os-Ir region, which are characteristic signatures of the i-process nucleosynthesis athigh neutron density, and which are all present in the abundance pattern of CS 31062-050. On theother hand, the one-zone model is not suited to correctly describe the evolution of the Fe groupelements and Pb, because of its unrealistic depletion of the n-capture seed elements. The observedabundance distributions in the CEMP-r/s stars HE 0338-3945 [17] and CS 22898-027 [2] are alsowell reproduced by the one-zone i-process model (Fig. 2).Figure 2: Fit of the results of the one-zone i-process model to the abundance distributions in the CEMP-r/sstars HE 0338-3945 (left panel, with data from [17]) and CS 22898-027 (right panel, with data from [2]) forthe indicated neutron densities (Nn) and exposures (tau).3. Nuclear physics uncertainties relevant to the i-processWhen comparing predicted nucleosynthesis yields from the i-process with observations, it is im-portant to understand what and how strongly nuclear physics uncertainties can affect the calculatedabundances of isotopes along the i-process nucleosynthesis path, that goes four to five species offthe valley of stability, as shown in Fig. 3. The first i-process uncertainty analysis was done in4The i-process and CEMP-r/s stars[4], where propagating systematic uncertainties of nuclear reaction cross sections from differenttheoretical models were investigated for elements of the second s-process peak. The analysis con-sidered Ba, La and Eu. The same i-process one-zone model described in Section 2 was used. Fig.6 in [4] shows changes of the [La/Eu] versus [Ba/La] ratio for cases when nuclear reaction ratesare estimated using three different theoretical Hauser-Feshbach codes: NON-SMOKER [23] (withrates from JINA REACLIB [8]), CoH 3.3 [18] and TALYS 1.4 [19]. The changes were found to bestrongly model-dependent. The resulting [Ba/La] and [La/Eu] ratios varied, at least, by factors of∼ 0.7 dex and ∼ 0.3 dex, respectively, between the models. The conclusions are that the nuclearphysics uncertainties strongly limit, at present, the predictive power of i-process simulations [4]and that it would be very important to further study the impact of nuclear physics uncertainties onthe i-process element production at various sites.Figure 3: The i-process nucleosynthesis path in the chart of isotopes with the color-coded mass fraction (X)of each species from our one-zone simulations.4. The i-process conditions in Sakurai’s object3D hydrodynamic simulations of mixing and H entrainment in the He PDCZ of the post-AGBstar Sakurai’s object are now possible at fine enough grids, so that the entrainment process isnumerically converged [24]. Combustion of ingested protons in the reaction 12C(p,γ)13N and itsenergy feedback on mixing was taken into account [14]. Two runs with 7683 and 15363 gridresolutions were performed and showed quantitatively quite similar behaviour. The continued Hingestion led to its accumulation in the upper part of the convective zone. At some point, a largeamount of H fuel ignited there and a resulting burning pocket rose and collided with the stiffconvective boundary. The mass-conservation law then made the pocket to spread and move alongthe boundary to the antipode, where the two fronts collided causing a strong downdraft of thematerial, and the process repeated several times. That new mode of non-radial instabilities, coinedby the authors “The Global Oscillation of Shell H-ingestion (GOSH)”, first appeared after nearly850 min in the 15363 run. It is shown in Fig. 4 at 1149 min [14, 24]. Such a GOSH event becomes5The i-process and CEMP-r/s starsFigure 4: A hemisphere with mixtures of entrained H-rich gas and He+C-rich gas of the He PDCZ. Theenergy release rate from the burning of ingested H is shown in very dark blue, yellow, and white.more and more violent as it is repeated for about a dozen times in a row, leading to entrainmentof entropy from the stable layer. In the run with 7683 grid resolution, a new upper convectiveboundary forms, while the run with 15363 cells has not been followed to this point yet. It is notclear yet how the GOSH can be taken into account in 1D stellar evolution calculations.5. ConclusionSimple one-zone simulations of i-process nucleosynthesis give yields that can fit almost all of theobserved heavy-element abundances from Y to Ir, within their observational errors, in some ofthe CEMP-r/s stars. The observed high enrichments with s- and r-process material in these starscannot be explained by the currently prevailing model that assumes a superposition of yields fromdifferent s- and r-process sources. These CEMP-r/s stars could therefore represent an i-processsite. Nuclear physics uncertainties are a main obstacle to obtaining reliable yields from the i-process nucleosynthesis calculations, as found by [4]. Mixing and H entrainment relevant to thei-process in the post-AGB star Sakurai’s object have recently been investigated, which has revealeda new mode of non-radial instabilities, the GOSH [14].Acknowledgments. This material is based upon work supported by the National Science Founda-tion under Grant No. PHY-1430152 (JINA Center for the Evolution of the Elements).References[1] D. M. Allen, S. G. Ryan, S. Rossi, T. C. Beers, S. A. Tsangarides, 2012, Elemental abundances andclassification of carbon-enhanced metal-poor stars, A&A,548, AA34[2] W. Aoki, J. E. Norris, S. G. Ryan, T. C. Beers, H. Ando, 2002, The Chemical Composition ofCarbon-rich, Very Metal Poor Stars: A New Class of Mildly Carbon Rich Objects without Excess ofNeutron-Capture Elements, APJ, 567, 11666The i-process and CEMP-r/s stars[3] T. C. Beers, N. Christlieb, 2005, The Discovery and Analysis of Very Metal-Poor Stars in the Galaxy,ARA&A , 43, 531[4] M. G. Bertolli, F. Herwig, M. Pignatari, T. Kawano, 2013, Systematic and correlated nuclearuncertainties in the i-process at the neutron shell closure N = 82, arXiv:1310.4578[5] S. Bisterzo, R. Gallino, O. Straniero, S. Cristallo, F. Kaeppeler, 2012, The s-process in low-metallicitystars - III. Individual analysis of CEMP-s and CEMP-s/r with asymptotic giant branch models,MNRAS, 422, 849[6] S. W. Campbell, M. Lugaro, A. I. Karakas, 2010, Evolution and nucleosynthesis of extremelymetal-poor and metal-free low- and intermediate-mass stars. II. s-process nucleosynthesis during thecore He flash, A&A, 522, LL6[7] J. J. Cowan, W. K. Rose, 1977, Production of C-14 and neutrons in red giants, APJ, 212, 149[8] R. H. Cyburt, A. M. Amthor, R. Ferguson, et al., 2010, The JINA REACLIB Database: Its RecentUpdates and Impact on Type-I X-ray Bursts, APJ, 189, 240[9] M. Y. Fujimoto, Y. Ikeda, Iben, I., Jr., 2000, The Origin of Extremely Metal-poor Carbon Stars andthe Search for Population III, ApJ, 529, L25[10] W. Fujiya, P. Hoppe, E. Zinner, M. Pignatari, F. Herwig, 2013, Evidence for Radiogenic Sulfur-32 inType AB Presolar Silicon Carbide Grains?, APJ, 776, LL29[11] N. Grevesse, A. Noels, 1993, Origin and Evolution of the Elements Cambridge University Press,Cambridge, 1993[12] F. Herwig, S. Diehl, C. L. Fryer, et al., 2008, Nucleosynthesis simulations for a wide range of nuclearproduction sites from NuGrid, in proceedings of Nuclei in the Cosmos (NIC X), PoS(NIC X)023[13] F. Herwig, M. Pignatari, P. R. Woodward, et al., 2011, Convective-reactive Proton-12C Combustion inSakurai’s Object (V4334 Sagittarii) and Implications for the Evolution and Yields from the FirstGenerations of Stars, APJ, 727, 89[14] F. Herwig, P. R. Woodward, P.-H. Lin, M. Knox, C. Fryer, 2014, Global Non-spherical Oscillations inThree-dimensional 4pi Simulations of the H-ingestion Flash, APJ, 792, LL3[15] M. Jadhav, M. Pignatari, F. Herwig, et al., 2013, Relics of Ancient Post-AGB Stars in a PrimitiveMeteorite, APJ, 777, LL27[16] J. A. Johnson, M. Bolte, 2004, The s-Process in Metal-Poor Stars: Abundances for 22Neutron-Capture Elements in CS 31062-050, APJ, 605, 462[17] K. Jonsell, P. S. Barklem, B. Gustafsson, et al., 2006, The Hamburg/ESO R-process enhanced starsurvey (HERES). III. HE 0338-3945 and the formation of the r + s stars, A&A, 451, 651[18] T. Kawano, S. Chiba, T. Maruyama, et al., 2004, Evaluation of Neutron Reaction Cross sections forAstrophysics, Proceedings of the 2003 Symposium on Nuclear Data, pp. 196-201, edited by T.Ohsawa and T. Fukahori,[19] A. Koning, S. Hilaire, S. Goriely, 2011, TALYS-1.4 (NRG, 2011),http://www.talys.eu/fileadmin/talys/user/docs/talys1.4.pdf[20] K. Lodders, 2003, Solar system abundances and condensation temperatures of the elements APJ, 591,1220[21] M. Lugaro, A. I. Karakas, R. J. Stancliffe, C. Rijs, 2012, The s-process in Asymptotic Giant BranchStars of Low Metallicity and the Composition of Carbon-enhanced Metal-poor Stars, APJ, 747, 27The i-process and CEMP-r/s stars[22] T. Masseron, J. A. Johnson, B. Plez, et al., 2010, A holistic approach to carbon-enhanced metal-poorstars, A&A,509, AA93[23] T. Rauscher, F.-K. Thielemann, 2000, Astrophysical Reaction Rates From Statistical ModelCalculations, Atomic Data and Nuclear Data Tables 75, 1[24] P. R. Woodward, F. Herwig, P.-H. Lin, 2013, Hydrodynamic Simulations of H Entrainment at the Topof He-shell Flash Convection, APJ, 798, 498

The i-process and CEMP-r/s stars

https://hull-repository.worktribe.com/preview/383738/i-process.pdf

The i-process and CEMP-r/s stars

Abstract

Similar works

Full text

Available Versions

Repository@Hull - Worktribe