Source author record

Christian Iliadis

Christian Iliadis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

23 published item(s)

preprint2024arXiv

Hydrodynamic Simulations of Oxygen-Neon Classical Novae as Galactic $^7$Li Producers and Potential Accretion Induced Collapse Progenitors

We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on oxygen-neon (ONe) white dwarfs (WDs). Using NOVA, a one-dimensional hydrodynamic computer code, we accrete Solar matter until the TNR is ongoing and then switch to a mixed composition. This approach is guided by the results of multi-dimensional studies of TNRs in WDs which find that sufficient mixing with WD core material occurs after the TNR is well underway, and levels of enrichment of the CNONeMg elements are reached that agree with observations of CN ejecta abundances. Because the amount of accreted material is inversely proportional to the oxygen abundance, by first accreting Solar matter, the amount of accreted material is larger than in those simulations with an initially enriched composition. We vary the mass of the WD (from 0.6 Msun to 1.35 Msun) and the composition of the mixed materials. Our results show large enrichments of 7Be in the ejected gases implying that ONe CNe and CO CNe (Starrfield et al. 2020) may be responsible for a significant fraction (about 100 Msun) of the galactic 7Li ( about 1000 Msun). The production of 22Na and 26Al in CN explosions and the gamma-ray emission predicted by our simulations is discussed. The WDs in all our simulations eject less material than they accrete and we predict that the WD is growing in mass as a consequence of the CN outburst. ONe CNe, therefore, may be an important channel for accretion induced collapse (AIC) events.

preprint2022arXiv

Hydrogen burning of $^{29}$Si and its impact on presolar stardust grains from classical novae

Presolar stardust grains found in primitive meteorites are believed to retain the isotopic composition of stellar outflows at the time of grain condensation. Therefore, laboratory measurements of their isotopic ratios represent sensitive probes for investigating open questions related to stellar evolution, stellar explosions, nucleosynthesis, mixing mechanisms, dust formation, and galactic chemical evolution. For a few selected presolar grains, classical novae have been discussed as a potential source. For SiC, silicate, and graphite presolar grains, the association is based on the observation of small $N(^{12}$C)/$N(^{13}$C) and $N(^{14}$N)/$N(^{15}$N) number abundance ratios compared to solar values, and abundance excesses in $^{30}$Si relative to $^{29}$Si, as previously predicted by models of classical novae. We report on a direct measurement of the $^{29}$Si(p,$γ$)$^{30}$P reaction, which strongly impacts simulated $δ^{29}$Si values from classical novae. Our new experimental $^{29}$Si(p,$γ$)$^{30}$P thermonuclear reaction rate differs from previous results by up to 50\% in the classical nova temperature range ($T$ $=$ $100$ $-$ $400$~MK), while the rate uncertainty is reduced by up to a factor of $3$. Using our new reaction rate in Monte Carlo reaction network and hydrodynamic simulations of classical novae, we estimate $δ^{29}$Si values with much reduced uncertainties. Our results establish $δ^{29}$Si values measured in presolar grains as a sensitive probe for assessing their classical nova paternity. We also demonstrate that $δ^{30}$Si values from nova simulations are presently not a useful diagnostic tool unless the large uncertainty of the $^{30}$P(p,$γ$)$^{31}$S reaction rate can be significantly reduced.

preprint2022arXiv

Thermonuclear reaction rate of $^{29}$Si(p,$γ$)$^{30}$P

The thermonuclear rate of the $^{29}$Si(p,$γ$)$^{30}$P reaction impacts the $^{29}$Si abundance in classical novae. A reliable reaction rate is essential for testing the nova paternity of presolar stardust grains. At present, the fact that no classical nova grains have been unambiguously identified in primitive meteorites among thousands of grains studied is puzzling, considering that classical novae are expected to be prolific producers of dust grains. We investigated the $^{29}$Si $+$ $p$ reaction at center-of-mass energies of $200$ $-$ $420$~keV, and present improved values for resonance energies, level excitation energies, resonance strengths, and branching ratios. One new resonance was found at a center-of-mass energy of $303$ keV. For an expected resonance at $215$~keV, an experimental upper limit could be determined for the strength. We evaluated the level structure near the proton threshold, and present new reaction rates based on all the available experimental information. Our new reaction rates have much reduced uncertainties compared to previous results at temperatures of $T$ $\ge$ $140$~MK, which are most important for classical nova nucleosynthesis. Future experiments to improve the reaction rates at lower temperatures are discussed.

preprint2021arXiv

On the analysis of signal peaks in pulse-height spectra

The estimation of the signal location and intensity of a peak in a pulse height spectrum is important for x-ray and $γ$-ray spectroscopy, charged-particle spectrometry, liquid chromatography, and many other subfields. However, both the "centroid" and "signal intensity" of a peak in a pulse-height spectrum are ill-defined quantities and different methods of analysis will yield different numerical results. Here, we apply three methods of analysis. Method A is based on simple count summation and is likely the technique most frequently applied in practice. The analysis is straightforward and fast, and does not involve any statistical modeling. We find that it provides reliable results only for high signal-to-noise data, but has severe limitations in all other cases. Method B employs a Bayesian model to extract signal counts and centroid from the measured total and background counts. The resulting values are derived from the respective posteriors and, therefore, have a rigorous statistical meaning. The method makes no assumptions about the peak shape. It yields reliable and relatively small centroid uncertainties. However, it provides relatively large signal count uncertainties. Method C makes a strong assumption regarding the peak shape by fitting a Gaussian function to the data. The fit is based again on a Bayesian model. Although Method C requires careful consideration of the Gaussian width (usually given by the detector resolution) used in the fitting, it provides reliable values and relatively small uncertainties both for the signal counts and the centroid.

preprint2020arXiv

Hierarchical Bayesian Thermonuclear Rate for the $^7$Be(n,p)$^7$Li Big Bang Nucleosynthesis Reaction

Big bang nucleosynthesis provides the earliest probe of standard model physics, at a time when the universe was less than a thousand seconds old. It determines the abundances of the lightest nuclides, which give rise to the subsequent history of the visible matter in the Universe. This work derives new $^7$Be(n,p)$^7$Li thermonuclear reaction rates based on all available experimental information. This reaction sensitively impacts the primordial abundances of $^{7}$Be and $^7$Li during big bang nucleosynthesis. We critically evaluate all available data and disregard experimental results that are questionable. For the nuclear model, we adopt an incoherent sum of single-level, two-channel R-matrix approximation expressions, which are implemented into a hierarchical Bayesian model, to analyze the remaining six data sets we deem most reliable. In the fitting of the data, we consistently model all known sources of uncertainty, including discrepant absolute normalizations of different data sets, and also take the variation of the neutron and proton channel radii into account, hence providing less biased estimates of the $^7$Be(n,p)$^7$Li thermonuclear rates. From the resulting posteriors, we extract R-matrix parameters ($E_r$, $γ^2_n$, $γ^2_p$) and derive excitation energies, partial and total widths. Our fit is sensitive to the contributions of the first three levels above the neutron threshold. Reaction rates were computed by integrating 10,000 samples of the reduced cross section. Our $^7$Be(n,p)$^7$Li thermonuclear rates have uncertainties between 1.5% and 2.0% at temperatures of $\leq$1 GK. We compare our rates to previous results and find that the $^7$Be(n,p)$^7$Li rates most commonly used in big bang simulations have too optimistic uncertainties.

preprint2019arXiv

Carbon-Oxygen Classical Novae are Galactic $^7$Li Producers as well as Potential Supernova Ia Progenitors

We report on studies of Classical Nova (CN) explosions where we follow the evolution of thermonuclear runaways (TNRs) on Carbon Oxygen (CO) white dwarfs (WDs). We vary both the mass of the WD (from 0.6 M$_\odot$ to 1.35 M$_\odot$) and the composition of the accreted material. Our simulations are guided by the results of multi-dimensional studies of TNRs in WDs that find sufficient mixing with WD core material occurs after the TNR is well underway, reaching levels of enrichment that agree with observations of CN ejecta abundances. We use NOVA (our 1-dimensional hydrodynamic code) to accrete solar matter until the TNR is ongoing and then switch to a mixed composition (either 25% WD material and 75% solar or 50% WD material and 50% solar). Because the amount of accreted material is inversely proportional to the initial $^{12}$C abundance, by first accreting solar matter the amount of material taking part in the outburst is larger than in those simulations where we assume a mixed composition from the beginning. Our results show large enrichments of $^7$Be in the ejected gases implying that CO CNe may be responsible for a significant fraction ($\sim$ 100 M$_\odot$) of the $^7$Li in the galaxy ($\sim$1000 M$_\odot$). In addition, although the ejected gases are enriched in WD material, the WDs in these simulations eject less material than they accrete. We predict that the WD is growing in mass as a consequence of the accretion-outburst-accretion cycle and CO CNe may be an important channel of Supernova Ia progenitors.

preprint2019arXiv

Thermonuclear Reaction Rate of $^{30}$Si(p,$γ$)$^{31}$P

Silicon synthesis in high-temperature hydrogen burning environments presents one possible avenue for the study of abundance anomalies in globular clusters. This was suggested in a previous study, which found that the large uncertainties associated with the $^{30}$Si(p,$γ$)$^{31}$P reaction rate preclude a firm understanding of the stellar conditions that give rise to the Mg-K anti-correlation observed in the globular cluster NGC 2419. In an effort to improve the reaction rate, we present new strength measurements of the $E_r^{lab} = 435$ keV and $E_r^{lab} = 501$ keV resonances in $^{30}$Si(p,$γ$)$^{31}$P. For the former, which was previously unobserved, we obtain a resonance strength of $ωγ= (1.28 \pm 0.25$) $\times 10^{-4}$ eV. For the latter, we obtain a value of $ωγ= (1.88 \pm 0.14)$ $\times 10^{-1}$ eV, which has a smaller uncertainty compared to previously measured strengths. Based on these results, the thermonuclear reaction rate has been re-evaluated. The impact of the new measurements is to lower the reaction rate by a factor of $\approx$10 at temperatures important to the study of NGC 2419. The rate uncertainty at these temperatures has also been reduced significantly.

preprint2016arXiv

Bayesian Estimation of Thermonuclear Reaction Rates

The problem of estimating non-resonant astrophysical S-factors and thermonuclear reaction rates, based on measured nuclear cross sections, is of major interest for nuclear energy generation, neutrino physics, and element synthesis. Many different methods have been applied in the past to this problem, almost all of them based on traditional statistics. Bayesian methods, on the other hand, are now in widespread use in the physical sciences. In astronomy, for example, Bayesian statistics is applied to the observation of extra-solar planets, gravitational waves, and type Ia supernovae. However, nuclear physics, in particular, has been slow to adopt Bayesian methods. We present astrophysical S-factors and reaction rates based on Bayesian statistics. We develop a framework that incorporates robust parameter estimation, systematic effects, and non-Gaussian uncertainties in a consistent manner. The method is applied to the d(p,$γ$)$^3$He, $^3$He($^3$He,2p)$^4$He, and $^3$He($α$,$γ$)$^7$Be reactions, important for deuterium burning, solar neutrinos, and big bang nucleosynthesis.

preprint2016arXiv

On Potassium and Other Abundance Anomalies of Red Giants in NGC 2419

Globular clusters are of paramount importance for testing theories of stellar evolution and early galaxy formation. Strong evidence for multiple populations of stars in globular clusters derives from observed abundance anomalies. A puzzling example is the recently detected Mg-K anticorrelation in NGC 2419. We perform Monte Carlo nuclear reaction network calculations to constrain the temperature-density conditions that gave rise to the elemental abundances observed in this elusive cluster. We find a correlation between stellar temperature and density values that provide a satisfactory match between simulated and observed abundances in NGC 2419 for all relevant elements (Mg, Si, K, Ca, Sc, Ti, and V). Except at the highest densities ($ρ\gtrsim 10^8$~g/cm$^3$), the acceptable conditions range from $\approx$ $100$~MK at $\approx$ $10^8$~g/cm$^3$ to $\approx$ $200$~MK at $\approx$ $10^{-4}$~g/cm$^3$. This result accounts for uncertainties in nuclear reaction rates and variations in the assumed initial composition. We review hydrogen burning sites and find that low-mass stars, AGB stars, massive stars, or supermassive stars cannot account for the observed abundance anomalies in NGC 2419. Super-AGB stars could be viable candidates for the polluter stars if stellar model parameters can be fine-tuned to produce higher temperatures. Novae, either involving CO or ONe white dwarfs, could be interesting polluter candidates, but a current lack of low-metallicity nova models precludes firmer conclusions. We also discuss if additional constraints for the first-generation polluters can be obtained by future measurements of oxygen, or by evolving models of second-generation low-mass stars with a non-canonical initial composition.

preprint2016arXiv

White Paper on Nuclear Astrophysics

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

preprint2014arXiv

Statistical Methods for Thermonuclear Reaction Rates and Nucleosynthesis Simulations

Rigorous statistical methods for estimating thermonuclear reaction rates and nucleosynthesis are becoming increasingly established in nuclear astrophysics. The main challenge being faced is that experimental reaction rates are highly complex quantities derived from a multitude of different measured nuclear parameters (e.g., astrophysical S-factors, resonance energies and strengths, particle and gamma-ray partial widths). We discuss the application of the Monte Carlo method to two distinct, but related, questions. First, given a set of measured nuclear parameters, how can one best estimate the resulting thermonuclear reaction rates and associated uncertainties? Second, given a set of appropriate reaction rates, how can one best estimate the abundances from nucleosynthesis (i.e., reaction network) calculations? The techniques described here provide probability density functions that can be used to derive statistically meaningful reaction rates and final abundances for any desired coverage probability. Examples are given for applications to s-process neutron sources, core-collapse supernovae, classical novae, and big bang nucleosynthesis.

preprint2013arXiv

Nuclear Mixing Meters for Classical Novae

Classical novae are caused by mass transfer episodes from a main sequence star onto a white dwarf via Roche lobe overflow. This material forms an accretion disk around the white dwarf. Ultimately, a fraction of this material spirals in and piles up on the white dwarf surface under electron-degenerate conditions. The subsequently occurring thermonuclear runaway reaches hundreds of megakelvin and explosively ejects matter into the interstellar medium. The exact peak temperature strongly depends on the underlying white dwarf mass, the accreted mass and metallicity, and the initial white dwarf luminosity. Observations of elemental abundance enrichments in these classical nova events imply that the ejected matter consists not only of processed solar material from the main sequence partner but also of material from the outer layers of the underlying white dwarf. This indicates that white dwarf and accreted matter mix prior to the thermonuclear runaway. The processes by which this mixing occurs require further investigation to be understood. In this work, we analyze elemental abundances ejected from hydrodynamic nova models in search of elemental abundance ratios that are useful indicators of the total amount of mixing. We identify the abundance ratios $Σ$CNO/H, Ne/H, Mg/H, Al/H, and Si/H as useful mixing meters in ONe novae. The impact of thermonuclear reaction rate uncertainties on the mixing meters is investigated using Monte Carlo post-processing network calculations with temperature-density evolutions of all mass zones computed by the hydrodynamic models. We find that the current uncertainties in the $^{30}$P($p$,$γ$)$^{31}$S rate influence the Si/H abundance ratio, but overall the mixing meters found here are robust against nuclear physics uncertainties. A comparison of our results with observations of ONe novae provides strong constraints for classical nova models.

preprint2012arXiv

Nuclear Thermometers for Classical Novae

Classical novae are stellar explosions occurring in binary systems, consisting of a white dwarf and a main sequence companion. Thermonuclear runaways on the surface of massive white dwarfs, consisting of oxygen and neon, are believed to reach peak temperatures of several hundred million kelvin. These temperatures are strongly correlated with the underlying white dwarf mass. The observational counterparts of such models are likely associated with outbursts that show strong spectral lines of neon in their shells (neon novae). The goals of this work are to investigate how useful elemental abundances are for constraining the peak temperatures achieved during these outbursts and determine how robust "nova thermometers" are with respect to uncertain nuclear physics input. We present updated observed abundances in neon novae and perform a series of hydrodynamic simulations for several white dwarf masses. We find that the most useful thermometers, N/O, N/Al, O/S, S/Al, O/Na, Na/Al, O/P, and P/Al, are those with the steepest monotonic dependence on peak temperature. The sensitivity of these thermometers to thermonuclear reaction rate variations is explored using post-processing nucleosynthesis simulations. The ratios N/O, N/Al, O/Na, and Na/Al are robust, meaning they are minimally affected by uncertain rates. However, their dependence on peak temperature is relatively weak. The ratios O/S, S/Al, O/P, and P/Al reveal strong dependences on temperature and the poorly known 30P(p,g)31S rate. We compare our model predictions to neon nova observations and obtain the following estimates for the underlying white dwarf masses: 1.34-1.35 solar masses (V838 Her), 1.18-1.21 solar masses (V382 Vel), <1.3 solar masses (V693 CrA), <1.2 solar masses (LMC 1990#1), and <1.2 solar masses (QU Vul).

preprint2012arXiv

Reaction rates for the s-process neutron source 22Ne+α

The 22Ne(α,n)25Mg reaction is an important source of neutrons for the s-process. In massive stars responsible for the weak component of the s-process, 22Ne(α,n)25Mg is the dominant source of neutrons, both during core helium burning and in shell carbon burning. For the main s-process component produced in Asymptotic Giant Branch (AGB) stars, the 13C(α,n)16O reaction is the dominant source of neutrons operating during the interpulse period, with the 22Ne+α source affecting mainly the s-process branchings during a thermal pulse. Rate uncertainties in the competing 22Ne(α,n)25Mg and 22Ne(α,γ)26Mg reactions result in large variations of s-process nucleosynthesis. Here, we present up-to-date and statistically rigorous 22Ne+α reaction rates using recent experimental results and Monte Carlo sampling. Our new rates are used in post-processing nucleosynthesis calculations both for massive stars and AGB stars. We demonstrate that the nucleosynthesis uncertainties arising from the new rates are dramatically reduced in comparison to previously published results, but several ambiguities in the present data must still be addressed. Recommendations for further study to resolve these issues are provided.

preprint2011arXiv

Nuclear astrophysics: the unfinished quest for the origin of the elements

Half a century has passed since the foundation of nuclear astrophysics. Since then, this discipline has reached its maturity. Today, nuclear astrophysics constitutes a multidisciplinary crucible of knowledge that combines the achievements in theoretical astrophysics, observational astronomy, cosmochemistry and nuclear physics. New tools and developments have revolutionized our understanding of the origin of the elements: supercomputers have provided astrophysicists with the required computational capabilities to study the evolution of stars in a multidimensional framework; the emergence of high-energy astrophysics with space-borne observatories has opened new windows to observe the Universe, from a novel panchromatic perspective; cosmochemists have isolated tiny pieces of stardust embedded in primitive meteorites, giving clues on the processes operating in stars as well as on the way matter condenses to form solids; and nuclear physicists have measured reactions near stellar energies, through the combined efforts using stable and radioactive ion beam facilities. This review provides comprehensive insight into the nuclear history of the Universe and related topics: starting from the Big Bang, when the ashes from the primordial explosion were transformed to hydrogen, helium, and few trace elements, to the rich variety of nucleosynthesis mechanisms and sites in the Universe. Particular attention is paid to the hydrostatic processes governing the evolution of low-mass stars, red giants and asymptotic giant-branch stars, as well as to the explosive nucleosynthesis occurring in core-collapse and thermonuclear supernovae, gamma-ray bursts, classical novae, X-ray bursts, superbursts, and stellar mergers.

preprint2011arXiv

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. 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 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(alpha,n)28Si, 24Mg(n,gamma)25Mg and 23Na(alpha,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 3.

preprint2011arXiv

Theoretical Evaluation of the Reaction Rates for 26Al(n,p)26Mg and 26Al(n,a)23Na

The reactions that destroy 26Al in massive stars have significance in a number of astrophysical contexts. We evaluate the reaction rates of 26Al(n,p)26Mg and 26Al(n,a)23Na using cross sections obtained from the codes EMPIRE and TALYS. These have been compared to the published rates obtained from the non-smoker code and to some experimental data. We show that the results obtained from EMPIRE and TALYS are comparable to those from non-smoker. We also show how the theoretical results vary with respect to changes in the input parameters. Finally, we present recommended rates for these reactions using the available experimental data and our new theoretical results.

preprint2010arXiv

Charged-Particle Thermonuclear Reaction Rates: I. Monte Carlo Method and Statistical Distributions

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 miu and sigma. 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 series (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

preprint2010arXiv

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

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 series (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 series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

preprint2010arXiv

Charged-Particle Thermonuclear Reaction Rates: III. Nuclear Physics Input

The nuclear physics input used to compute the Monte Carlo reaction rates and probability density functions that are tabulated in the second paper of this series (Paper II) is presented. Specifically, we publish the input files to the Monte Carlo reaction rate code RatesMC, which is based on the formalism presented in the first paper of this series (Paper I). This data base contains overwhelmingly experimental nuclear physics information. The survey of literature for this review was concluded in November 2009.

preprint2010arXiv

Charged-Particle Thermonuclear Reaction Rates: IV. Comparison to Previous Work

We compare our Monte Carlo reaction rates (see Paper II of this series) to previous results that were obtained by using the classical method of computing thermonuclear reaction rates. For each reaction, the comparison is presented using two types of graphs: the first shows the change in reaction rate uncertainties, while the second displays our new results normalized to the previously recommended reaction rate. We find that the rates have changed significantly for almost all reactions considered here. The changes are caused by (i) our new Monte Carlo method of computing reaction rates (see Paper I of this series), and (ii) newly available nuclear physics information (see Paper III of this series).

preprint2010arXiv

Hydrodynamic Models of Type I X-Ray Bursts: Metallicity Effects

Type I X-ray bursts are thermonuclear stellar explosions driven by charged-particle reactions. In the regime for combined H/He-ignition, the main nuclear flow is dominated by the rp-process (rapid proton-captures and beta+ decays), the 3 alpha-reaction, and the alpha-p-process (a suite of (alpha,p) and (p,gamma) reactions). The main flow is expected to proceed away from the valley of stability, eventually reaching the proton drip-line beyond A = 38. Detailed analysis of the relevant reactions along the main path has only been scarcely addressed, mainly in the context of parameterized one-zone models. In this paper, we present a detailed study of the nucleosynthesis and nuclear processes powering type I X-ray bursts. The reported 11 bursts have been computed by means of a spherically symmetric (1D), Lagrangian, hydrodynamic code, linked to a nuclear reaction network that contains 325 isotopes (from 1H to 107Te), and 1392 nuclear processes. These evolutionary sequences, followed from the onset of accretion up to the explosion and expansion stages, have been performed for 2 different metallicities to explore the dependence between the extension of the main nuclear flow and the initial metal content. We carefully analyze the dominant reactions and the products of nucleosynthesis, together with the the physical parameters that determine the light curve (including recurrence times, ratios between persistent and burst luminosities, or the extent of the envelope expansion). Results are in qualitative agreement with the observed properties of some well-studied bursting sources. Leakage from the predicted SbSnTe-cycle cannot be discarded in some of our models. Production of 12C (and implications for the mechanism that powers superbursts), light p-nuclei, and the amount of H left over after the bursting episodes will also be discussed.

preprint2009arXiv

Introduction to nuclear astrophysics

In the first lecture of this volume, we will present the basic fundamental ideas regarding nuclear processes occurring in stars. We start from stellar observations, will then elaborate on some important quantum-mechanical phenomena governing nuclear reactions, continue with how nuclear reactions proceed in a hot stellar plasma and, finally, we will provide an overview of stellar burning stages. At the end, the current knowledge regarding the origin of the elements is briefly summarized. This lecture is directed towards the student of nuclear astrophysics. Our intention is to present seemingly unrelated phenomena of nuclear physics and astrophysics in a coherent framework.