Christian Iliadis

Light nuclei in galactic globular clusters: constraints on the self-enrichment scenario from nucleosynthesis

Aims. Hydrogen-burning is the root cause of the star-to-star abundance variations of light nuclei in Galactic globular clusters (GC). In the present work we constrain the physical conditions that gave rise to the observed abundance patterns of Li, C, N, O, Na, Mg, Al, as well as Mg isotopes in the typical case of NGC 6752. Methods. We perform nucleosynthesis calculations at constant temperature, adopting realistic initial abundances for the proto-cluster gas. We use a detailed nuclear reaction network and state-of-the-art nuclear reaction rates. Results. Although simplistic, our analysis provides original results and new constraints on the self-enrichment scenario for GCs. Our parametric calculations allow us to determine a narrow range of temperature where the observed extreme abundances of all light elements and isotopes in NGC 6752 are nicely reproduced simultaneously. This agreement is obtained after mixing of the H-processed material with ∼30% of unprocessed gas. We show that the observed C-N, O-Na, Mg-Al, Li-Na and F-Na anticorrelations as well as the behaviour of the Mg isotopes can be recovered by assuming mixing with various dilution factors. Li production by the stars that build up the other abundance anomalies is not mandatory in the case of NGC 6752. Conclusions. Observations of O, Na, Mg and Al constrain the temperature range for H-burning; such temperatures are encountered in the two main candidate “polluters” proposed for GCs, namely massive AGBs and the most massive main-sequence stars. Furthermore, observations require dilution of H-burning processed material with pristine one. They provide no clue, however, as to the nature of the unprocessed material required for mixing. The complementary observations of the fragile Li and F clearly point to ISM origin for the mixed material.

Proton-induced Thermonuclear Reaction Rates for A = 20–40 Nuclei

Proton-induced reaction rates on 26 stable and 29 unstable target nuclei in the mass A = 20–40 region have been evaluated and compiled. Recommended reaction rates, assuming that all interacting nuclei are in the ground state, are presented in tabular form on a temperature grid in the range T = 0.01–10.0 GK. Most reaction rates involving stable targets were normalized to a set of measured standard resonance strengths in the sd shell. For the majority of reaction rates, experimental information from transfer reaction studies has been used consistently. Our results are compared with recent statistical model (Hauser-Feshbach) calculations. Reaction rate uncertainties are presented and amount to several orders of magnitude for many of the reactions. Several of these reaction rates and/or their corresponding uncertainties deviate from results of previous compilations. In most cases, the deviations are explained by the fact that new experimental information became available recently. Examples are given for calculating reaction rates and reverse reaction rates for thermally excited nuclei from the present results. The survey of literature for this review was concluded in 2000 August.

The Effects of Thermonuclear Reaction-Rate Variations on Nova Nucleosynthesis: A Sensitivity Study

We investigate the effects of thermonuclear reaction-rate uncertainties on nova nucleosynthesis. One-zone nucleosynthesis calculations have been performed by adopting temperature-density-time profiles of the hottest hydrogen-burning zone (i.e., the region in which most of the nucleosynthesis takes place). We obtain our profiles from seven different, recently published, hydrodynamic nova simulations covering peak temperatures in the range from Tpeak = 0.145 to 0.418 GK. For each of these profiles, we individually varied the rates of 175 reactions within their associated errors and analyzed the resulting abundance changes of 142 isotopes in the mass range below A = 40. In total, we performed ≈7350 nuclear reaction network calculations. We use the most recent thermonuclear reaction-rate evaluations for the mass ranges A = 1-20 and 20-40. For the theoretical astrophysicist, our results indicate the extent to which nova nucleosynthesis calculations depend on currently uncertain nuclear physics input, while for the experimental nuclear physicist, our results represent at least a qualitative guide for future measurements at stable and radioactive ion beam facilities. We find that present reaction-rate estimates are reliable for predictions of Li, Be, C, and N abundances in nova nucleosynthesis. However, rate uncertainties of several reactions have to be reduced significantly in order to predict more reliable O, F, Ne, Na, Mg, Al, Si, S, Cl, and Ar abundances. Results are presented in tabular form for each adopted nova simulation.

Charged-Particle Thermonuclear Reaction Rates: II. Tables and Graphs of Reaction Rates and Probability Density Functions

Abstract Numerical values of charged-particle thermonuclear reaction rates for nuclei in the A = 14 to 40 region are tabulated. The results are obtained using a method, based on Monte Carlo techniques, that has been described in the preceding paper of this issue (Paper I). We present a low rate, median rate and high rate which correspond to the 0.16, 0.50 and 0.84 quantiles, respectively, of the cumulative reaction rate distribution. The meaning of these quantities is in general different from the commonly reported, but statistically meaningless expressions, “lower limit”, “nominal value” and “upper limit” of the total reaction rate. In addition, we approximate the Monte Carlo probability density function of the total reaction rate by a lognormal distribution and tabulate the lognormal parameters μ and σ at each temperature. We also provide a quantitative measure (Anderson–Darling test statistic) for the reliability of the lognormal approximation. The user can implement the approximate lognormal reaction rate probability density functions directly in a stellar model code for studies of stellar energy generation and nucleosynthesis. For each reaction, the Monte Carlo reaction rate probability density functions, together with their lognormal approximations, are displayed graphically for selected temperatures in order to provide a visual impression. Our new reaction rates are appropriate for bare nuclei in the laboratory. The nuclear physics input used to derive our reaction rates is presented in the subsequent paper of this issue (Paper III). In the fourth paper of this issue (Paper IV) we compare our new reaction rates to previous results.

THE EFFECTS OF THERMONUCLEAR REACTION RATE VARIATIONS ON 26Al PRODUCTION IN MASSIVE STARS: A SENSITIVITY STUDY

We investigate the effects of thermonuclear reaction rate variations on 26Al production in massive stars. The dominant production sites in such events were recently investigated by using stellar model calculations: explosive neon-carbon burning, convective shell carbon burning, and convective core hydrogen burning. Post-processing nucleosynthesis calculations are performed for each of these sites by adopting temperature-density-time profiles from recent stellar evolution models. For each profile, we individually multiplied the rates of all relevant reactions by factors of 10, 2, 0.5, and 0.1, and analyzed the resulting abundance changes of 26Al. In total, we performed ≈900 nuclear reaction network calculations. Our simulations are based on a next-generation nuclear physics library, called STARLIB, which contains a recent evaluation of Monte Carlo reaction rates. Particular attention is paid to quantifying the rate uncertainties of those reactions that most sensitively influence 26Al production. For stellar modelers our results indicate to what degree predictions of 26Al nucleosynthesis depend on currently uncertain nuclear physics input, while for nuclear experimentalists our results represent a guide for future measurements. We also investigate equilibration effects of 26Al. In all previous massive star investigations, either a single species or two species of 26Al were taken into account, depending on whether thermal equilibrium was achieved or not. These are two extreme assumptions, and in a hot stellar plasma the ground and isomeric states may communicate via γ-ray transitions involving higher-lying 26Al levels. We tabulate the results of our reaction rate sensitivity study for each of the three distinct massive star sites referred to above. It is found that several current reaction rate uncertainties influence the production of 26Al. Particularly important reactions are 26Al(n,p)26Mg, 25Mg(α,n)28Si, 24Mg(n,γ)25Mg, and 23Na(α,p)26Mg. These reactions should be prime targets for future measurements. Overall, we estimate that the nuclear physics uncertainty of the 26Al yield predicted by the massive star models explored here amounts to about a factor of three. We also find that taking the equilibration of 26Al levels explicitly into account in any of the massive star sites investigated here has only minor effects on the predicted 26Al yields. Furthermore, we provide for the interested reader detailed comments regarding the current status of certain reactions, including 12C(12C,n)23Mg, 23Na(α,p)26Mg, 25Mg(α,n)28Si, 26Al m (p,γ)27Si, 26Al(n,p)26Mg, and 26Al(n,α)23Na.

The effects of variations in nuclear processes on type I X-ray burst nucleosynthesis

Type I X-ray bursts are violent stellar events that take place in the H/He-rich envelopes of accreting neutron stars. We have investigated the role played by uncertainties in nuclear processes on the nucleosynthesis accompanying these explosive phenomena. Two different approaches have been adopted, in the framework of postprocessing calculations. In the first one, nuclear rates are varied individually within uncertainties. Ten different models, covering the characteristic parameter space for these stellar events, have been considered. The second, somewhat complementary approach involves a Monte Carlo code in which all nuclear rates are randomly varied within uncertainty limits simultaneously. All in all, about 50,000 postprocessing calculations, with a network containing 606 nuclides (H to 113Xe) and more than 3500 nuclear processes, have been performed in this work. A brief comparison between both procedures is outlined together with an overall account of the key nuclear reactions whose uncertainties have the largest impact in our X-ray burst nucleosynthesis studies.

Starlib: A next-generation reaction-rate library for nuclear astrophysics

STARLIB is a next-generation, all-purpose nuclear reaction-rate library. For the first time, this library provides the rate probability density at all temperature grid points for convenient implementation in models of stellar phenomena. The recommended rate and its associated uncertainties are also included. Currently, uncertainties are absent from all other rate libraries, and, although estimates have been attempted in previous evaluations and compilations, these are generally not based on rigorous statistical definitions. A common standard for deriving uncertainties is clearly warranted. STARLIB represents a first step in addressing this deficiency by providing a tabular, up-to-date database that supplies not only the rate and its uncertainty but also its distribution. Because a majority of rates are lognormally distributed, this allows the construction of rate probability densities from the columns of STARLIB. This structure is based on a recently suggested Monte Carlo method to calculate reaction rates, where uncertainties are rigorously defined. In STARLIB, experimental rates are supplemented with: (1) theoretical TALYS rates for reactions for which no experimental input is available, and (2) laboratory and theoretical weak rates. STARLIB includes all types of reactions of astrophysical interest to Z = 83, such as (p, γ), (p, α), (α, n), and corresponding reverse rates. Strong rates account for thermal target excitations. Here, we summarize our Monte Carlo formalism, introduce the library, compare methods of correcting rates for stellar environments, and discuss how to implement our library in Monte Carlo nucleosynthesis studies. We also present a method for accessing STARLIB on the Internet and outline updated Monte Carlo-based rates.

THE EFFECTS OF THE pep NUCLEAR REACTION AND OTHER IMPROVEMENTS IN THE NUCLEAR REACTION RATE LIBRARY ON SIMULATIONS OF THE CLASSICAL NOVA OUTBURST

Nova explosions occur on the white dwarf (WD) component of a cataclysmic variable binary stellar system which is accreting matter lost by its companion. When sufficient material has been accreted by the WD, a thermonuclear runaway (TNR) occurs and ejects material in what is observed as a classical nova (CN) explosion. We have continued our studies of TNRs on 1.25 M ☉ and 1.35 M ☉ WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our one-dimensional, fully implicit, hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. Because the evolutionary sequences on the 1.35 M ☉ WD reach higher temperatures, the effects of library changes are more important for this mass. In addition, as a result of our improvements, we discovered that the pep reaction (p + e – + p → d + ν) was not included in our previous studies of CN explosions (or to the best of our knowledge those of other investigators). Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed 104 g cm–3 and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. Since it is only the p-p chain that is operating during most of the accretion phase prior to the final rise to the TNR, the effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper.

Charged-particle thermonuclear reaction rates: I. Monte Carlo method and statistical distributions

Abstract A method based on Monte Carlo techniques is presented for evaluating thermonuclear reaction rates. We begin by reviewing commonly applied procedures and point out that reaction rates that have been reported up to now in the literature have no rigorous statistical meaning. Subsequently, we associate each nuclear physics quantity entering in the calculation of reaction rates with a specific probability density function, including Gaussian, lognormal and chi-squared distributions. Based on these probability density functions the total reaction rate is randomly sampled many times until the required statistical precision is achieved. This procedure results in a median (Monte Carlo) rate which agrees under certain conditions with the commonly reported recommended “classical” rate. In addition, we present at each temperature a low rate and a high rate, corresponding to the 0.16 and 0.84 quantiles of the cumulative reaction rate distribution. These quantities are in general different from the statistically meaningless “minimum” (or “lower limit”) and “maximum” (or “upper limit”) reaction rates which are commonly reported. Furthermore, we approximate the output reaction rate probability density function by a lognormal distribution and present, at each temperature, the lognormal parameters μ and σ . The values of these quantities will be crucial for future Monte Carlo nucleosynthesis studies. Our new reaction rates, appropriate for bare nuclei in the laboratory , are tabulated in the second paper of this issue (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this issue (Paper III). In the fourth paper of this issue (Paper IV) we compare our new reaction rates to previous results.

Explosive hydrogen burning of 27Si, 31S, 35Ar, and 39Ca in novae and X-ray bursts

Stellar reaction rates for proton captures on the nuclei 27Si, 31S, 35Ar, and 39Ca are estimated from the most recent nuclear structure information available. Reliable mirror-state correspondences are found by using the isobaric multiplet mass equation. An improved method for calculating proton partial widths is applied. Systematic comparisons of excitation energies, spectroscopic factors, proton partial widths, and γ-ray partial widths for states of the same isospin multiplet are presented. Stellar reaction-rate uncertainties are deduced, and our reaction rates are compared to previous estimates. Reaction network calculations are performed to investigate implications of the new reaction rates for nucleosynthesis in novae and X-ray bursts. Our recommended reaction rates are varied within their assigned uncertainties, and we find only minor effects on the nuclear energy generation and the final abundances after the outbursts. Thus, contrary to previous claims, we find no compelling reason for measuring the proton capture reactions on 27Si, 31S, 35Ar, and 39Ca by using radioactive ion beams.