Source author record

Richard Longland

Richard Longland 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

13works
6topics
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

13 published item(s)

preprint2026arXiv

Propagating Uncertainties from Nuclear Physics to Gamma-rays in Core Collapse Supernovae

Nuclear yields are powerful probes of supernova explosions, their engines and their progenitors. In addition, as we improve our understanding of these explosions, we can use nuclear yields to probe dense matter and neutrino physics, both of which play a critical role in the central supernova engine. Especially with upcoming gamma-ray detectors that can directly detect radioactive isotopes out to increasing distances from gamma-rays emitted during their decay, nuclear yields have the potential to provide some of the most direct probes of supernova engines and stellar burning. To utilize these probes, we must understand and limit the uncertainties in their production. Uncertainties in the nuclear physics can be minimized by combining both laboratory experiments and nuclear theory. Similarly, astrophysical uncertainties caused by simplified explosion trajectories can be minimized by higher-fidelity stellar-evolution and supernova-engine models. This paper reviews the physics and astrophysics uncertainties in modeling nucleosynthetic yields, identifying the key areas of study needed to maximize the potential of supernova yields as probes of astrophysical transients and dense-matter physics.

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.

preprint2015arXiv

New reaction rates for improved primordial D/H calculation and the cosmic evolution of deuterium

Primordial or big bang nucleosynthesis (BBN) is one of the three historical strong evidences for the big bang model. Standard BBN is now a parameter free theory, since the baryonic density of the Universe has been deduced with an unprecedented precision from observations of the anisotropies of the cosmic microwave background (CMB) radiation. There is a good agreement between the primordial abundances of 4He, D, 3He and 7Li deduced from observations and from primordial nucleosynthesis calculations. However, the 7Li calculated abundance is significantly higher than the one deduced from spectroscopic observations and remains an open problem. In addition, recent deuterium observations have drastically reduced the uncertainty on D/H, to reach a value of 1.6%. It needs to be matched by BBN predictions whose precision is now limited by thermonuclear reaction rate uncertainties. This is especially important as many attempts to reconcile Li observations with models lead to an increased D prediction. Here, we re-evaluates the D(p,g)3He, D(d,n)3He and D(d,p)3H reaction rates that govern deuterium destruction, incorporating new experimental data and carefully accounting for systematic uncertainties. Contrary to previous evaluations, we use theoretical ab initio models for the energy dependence of the S-factors. As a result, these rates increase at BBN temperatures, leading to a reduced value of D/H = (2.45$\pm0.10)\times10^{-5}$ (2$σ$), in agreement with observations.

preprint2014arXiv

Performance Improvements for Nuclear Reaction Network Integration

Aims: The aim of this work is to compare the performance of three reaction network integration methods used in stellar nucleosynthesis calculations. These are the Gear's backward differentiation method, Wagoner's method (a 2nd-order Runge-Kutta method), and the Bader-Deuflehard semi-implicit multi-step method. Methods: To investigate the efficiency of each of the integration methods considered here, a test suite of temperature and density versus time profiles is used. This suite provides a range of situations ranging from constant temperature and density to the dramatically varying conditions present in white dwarf mergers, novae, and x-ray bursts. Some of these profiles are obtained separately from full hydrodynamic calculations. The integration efficiencies are investigated with respect to input parameters that constrain the desired accuracy and precision. Results: Gear's backward differentiation method is found to improve accuracy, performance, and stability in integrating nuclear reaction networks. For temperature-density profiles that vary strongly with time, it is found to outperform the Bader-Deuflehard method (although that method is very powerful for more smoothly varying profiles). Wagoner's method, while relatively fast for many scenarios, exhibits hard-to-predict inaccuracies for some choices of integration parameters owing to its lack of error estimations.

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.

preprint2012arXiv

Constraints on the mass and radius of the accreting neutron star in the Rapid Burster

The Rapid Burster (MXB 1730-335) is a unique object, showing both type I and type II X-ray bursts. A type I burst of the Rapid Burster was observed with Swift/XRT on 2009 March 5, showing photospheric radius expansion for the first time in this source. We report here on the mass and radius determination from this photospheric radius expansion burst using a Bayesian approach. After marginalization over the likely distance of the system (5.8-10 kpc) we obtain M=1.1+/-0.3 M_sun and R=9.6+/-1.5 km (1-sigma uncertainties) for the compact object, ruling out the stiffest equations of state for the neutron star. We study the sensitivity of the results to the distance, the color correction factor, and the hydrogen mass fraction in the envelope. We find that only the distance plays a crucial role.

preprint2012arXiv

Lithium production in the merging of white dwarf stars

The origin of R Coronae Borealis stars has been elusive for over 200 years. Currently, two theories for their formation have been presented. These are the Final Flash scenario, in which a dying asymptotic giant branch (AGB) star throws off its atmosphere to reveal the hydrogen poor, heavily processed material underneath, and the double degenerate scenario, in which two white dwarfs merge to produce a new star with renewed vigour. Some theories predict that the temperatures reached during the latter scenario would destroy any lithium originally present in the white dwarfs. The observed lithium content of some R Coronae Borealis stars, therefore, is often interpreted as an indication that the Final Flash scenario best describes their formation. In this paper, it is shown that lithium production can, indeed, occur in the merging of a helium white dwarf with a carbon-oxygen white dwarf if their chemical composition, particularly that of 3He, is fully considered. The production mechanism is described in detail, and the sensitivity of lithium production to the merger environment is investigated. Nucleosynthesis post-processing calculations of smoothed-particle hydrodynamics (SPH) tracer particles are performed to show that any lithium produced in these environments will be concentrated towards the cloud of material surrounding the R CrB star. Measurements of the lithium content of these clouds would, therefore, provide a valuable insight into the formation mechanism of R CrB stars.

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.

preprint2012arXiv

Recommendations for Monte Carlo nucleosynthesis sampling (Research Note)

Context: Recent reaction rate evaluations include reaction rate uncertainties that have been determined in a statistically meaningful manner. Furthermore, reaction rate probability density distributions have been determined and published in the form of lognormal parameters with the specific goal of pursuing Monte Carlo nucleosynthesis studies. Aims: To test and assess different methods of randomly sampling over reaction rate probability densities and to determine the most accurate method for estimating elemental abundance uncertainties. Methods: Experimental Monte Carlo reaction rates are first computed for the 22Ne+alpha, 20Ne(p,g)21Na, 25Mg(p,g)26Al, and 18F(p,alpha)15O reactions, which are used to calculate reference nucleosynthesis yields for 16 nuclei affected by nucleosynthesis in massive stars and classical novae. Five different methods of randomly sampling over these reaction rate probability distributions are then developed, tested, and compared with the reference nucleosynthesis yields. Results: Given that the reaction rate probability density distributions can be described accurately with a lognormal distribution, Monte Carlo nucleosynthesis variations arising from the parametrised estimates for the reaction rate variations agree remarkably well with those obtained from the true rate samples. Most significantly, the most simple parametrisation agrees within just a few percent, meaning that Monte Carlo nucleosynthesis studies can be performed reliably using lognormal parametrisations of reaction rate probability density functions.

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).