I review measurements of the most important reactions involved in the creation and destruction of Li and Be during big bang nucleosynthesis (BBN) as well as their uncertainties and the relative contributions they make to the uncertainty in the primordial $^7$Li abundance ($^7$Li/H). Examining the sensitivity of calculated $^7$Li/H to these reactions as predicted by different BBN codes I find no significant differences. I compare my calculation of primordial $^7$Li/H to some recently published values and conclude that in the absence of a major undetected experimental blunder, nuclear physics uncertainties cannot account for the cosmological Li problem. With an estimated 13% uncertainty in the calculated abundance, the discrepancy with observation amounts to some $4.6\sigma$.

In a recent Letter entitled "Measurement of d + 7Be Cross Sections for Big-Bang Nucleosynthesis", Rijal et al. describe the results of an experiment that put the d + 7Be reaction rate on firmer ground. The erratum of Rijal et al. claims that including this reaction rate in big bang nucleosynthesis brings a significant reduction of up to 8.1% of lithium production relative to calculations that exclude this reaction. In contrast, using several well tested BBN codes, we confirm earlier studies showing that the effect is less than 1%, so that this measurement cannot alleviate the cosmological lithium problem.

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

The astrophysical $S$-factor for the radiative proton capture reaction on $^7$Be ($S_{17}$) at low energies is affected by the $s$-wave scattering lengths. We report the measurement of elastic and inelastic scattering cross sections for the $^7$Be+p system in the center-of-mass energy range 0.474 - 2.740 MeV and center-of-mass angular range of 70$^\circ$- 150$^\circ$. A radioactive $^7$Be beam produced at Oak Ridge National Laboratory's (ORNL) Holifield Radioactive Ion Beam Facility was accelerated and bombarded a thin polypropylene (CH$_{2}$)$_\text n$ target. Scattered ions were detected in the segmented Silicon Detector Array. Using an $\textit{R}$-matrix analysis of ORNL and Louvain-la-Neuve cross section data, the $s$-wave scattering lengths for channel spins 1 and 2 were determined to be 17.34$^{+1.11}_{-1.33}$ and -3.18$^{+0.55}_{-0.50}$ fm, respectively. The uncertainty in the $s$-wave scattering lengths reported in this work is smaller by a factor of 5-8 compared to the previous measurement, which may reduce the overall uncertainty in $S_{17}$ at zero energy. The level structure of $^8$B is discussed based upon the results from this work. Evidence for the existence of 0$^+$ and 2$^+$ levels in $^8$B at 1.9 and 2.21 MeV, respectively, is observed.

According to sensitivity studies, the $^{38}\mathrm{K}\left( p, \gamma \right){}^{39}\mathrm{Ca}$ reaction has a significant influence on $\mathrm{Ar}$, $\mathrm{K}$, and $\mathrm{Ca}$ production in classical novae. In order to constrain the rate of this reaction, we have performed a direct measurement of the strengths of three candidate $\ell = 0$ resonances within the Gamow window, at $386 \pm 10~\mathrm{keV}$, $515 \pm 10~\mathrm{keV}$, and $689 \pm 10~\mathrm{keV}$. The experiment was performed in inverse kinematics using a beam of unstable $^{38}\mathrm{K}$ impinged on a windowless $\mathrm{H}_2$ target. The $^{39}\mathrm{Ca}$ recoils and prompt $\gamma$ rays from $^{38}\mathrm{K}\left( p, \gamma \right){}^{39}\mathrm{Ca}$ reactions were detected in coincidence using a recoil mass separator and a BGO array, respectively. For the $689$ keV resonance, we observed a clear recoil-$\gamma$ coincidence signal and extracted resonance strength and energy values of $120^{+50}_{-30}~\mathrm{(stat.)}^{+20}_{-60}~\mathrm{(sys.)}~\mathrm{meV}$ and $679^{+2}_{-1}~\mathrm{(stat.)} \pm 1~\mathrm{(sys.)}~\mathrm{keV}$, respectively. We also performed a singles analysis, extracting a resonance strength of $120 \pm 20~\mathrm{(stat.)} \pm 15~\mathrm{(sys.)}~\mathrm{meV}$, consistent with the coincidence result. For the $386$ keV and $515$ keV resonances, we extract $90\%$ confidence level upper limits of $2.54$ meV and $18.4$ meV, respectively. We have established a new recommended $^{38}\mathrm{K}(p, \gamma){}^{39}\mathrm{Ca}$ rate based on experimental information, which reduces overall uncertainties near the peak temperatures of nova burning by a factor of ${\sim} 250$. Using the rate obtained in this work in model calculations of the hottest oxygen-neon novae reduces overall uncertainties on $\mathrm{Ar}$, $\mathrm{K}$, and $\mathrm{Ca}$ synthesis to factors of $15$ or less in all cases.

The nucleosynthesis of elements beyond iron is dominated by neutron captures in the s and r processes. However, 32 stable, proton-rich isotopes cannot be formed during those processes, because they are shielded from the s-process flow and r-process beta-decay chains. These nuclei are attributed to the p and rp process. For all those processes, current research in nuclear astrophysics addresses the need for more precise reaction data involving radioactive isotopes. Depending on the particular reaction, direct or inverse kinematics, forward or time-reversed direction are investigated to determine or at least to constrain the desired reaction cross sections. The Facility for Antiproton and Ion Research (FAIR) will offer unique, unprecedented opportunities to investigate many of the important reactions. The high yield of radioactive isotopes, even far away from the valley of stability, allows the investigation of isotopes involved in processes as exotic as the r or rp processes.

The 33S(p,\gamma)34Cl reaction is important for constraining predictions of certain isotopic abundances in oxygen-neon novae. Models currently predict as much as 150 times the solar abundance of 33S in oxygen-neon nova ejecta. This overproduction factor may, however, vary by orders of magnitude due to uncertainties in the 33S(p,\gamma)34Cl reaction rate at nova peak temperatures. Depending on this rate, 33S could potentially be used as a diagnostic tool for classifying certain types of presolar grains. Better knowledge of the 33S(p,\gamma)34Cl rate would also aid in interpreting nova observations over the S-Ca mass region and contribute to the firm establishment of the maximum endpoint of nova nucleosynthesis. Additionally, the total S elemental abundance which is affected by this reaction has been proposed as a thermometer to study the peak temperatures of novae. Previously, the 33S(p,\gamma)34Cl reaction rate had only been studied directly down to resonance energies of 432 keV. However, for nova peak temperatures of 0.2-0.4 GK there are 7 known states in 34Cl both below the 432 keV resonance and within the Gamow window that could play a dominant role. Direct measurements of the resonance strengths of these states were performed using the DRAGON recoil separator at TRIUMF. Additionally two new states within this energy region are reported. Several hydrodynamic simulations have been performed, using all available experimental information for the 33S(p,\gamma)34Cl rate, to explore the impact of the remaining uncertainty in this rate on nucleosynthesis in nova explosions. These calculations give a range of ~ 20-150 for the expected 33S overproduction factor, and a range of ~ 100-450 for the 32S/33S ratio expected in ONe novae.

Data from a recent 9Be(3He,t)9B measurement are used to rule out a possible solution to the cosmological lithium problem based on conventional nuclear physics.

We present a Monte Carlo calculation of the astrophysical rate of the 15O(alpha,gamma)19Ne reaction based on an evaluation of published experimental data. By considering the likelihood distributions of individual resonance parameters derived from measurements, estimates of upper and lower limits on the reaction rate at the 99.73% confidence level are derived in addition to the recommended, median value. These three reaction rates are used as input for three separate calculations of Type I x-ray bursts using spherically symmetric, hydrodynamic simulations of an accreting neutron star. In this way the influence of the 15O(alpha,gamma)19Ne reaction rate on the peak luminosity, recurrence time, and associated nucleosynthesis in models of Type I x-ray bursts is studied. Contrary to previous findings, no substantial effect on any of these quantities is observed in a sequence of four bursts when varying the reaction rate between its lower and upper limits. Rather, the differences in these quantities are comparable to the burst-to-burst variations with a fixed reaction rate, indicating that uncertainties in the 15O(alpha,gamma)19Ne reaction rate do not strongly affect the predictions of this Type I x-ray burst model.

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

In both the Sun and the early universe, the 3He(alpha,gamma)7Be reaction plays a key role. The rate of this reaction is the least certain nuclear input needed to calculate both the primordial 7Li abundance in big bang nucleosynthesis (BBN) and the solar neutrino flux. Taking advantage of several recent highly precise experiments, we analyse modern 3He(alpha,gamma)7Be data using a robust and minimally model dependent approach capable of handling discrepant data sets dominated by systematic rather than statistical errors. We find S34(0)=0.580 pm 0.043(0.054) keV b at the 68.3(95.4)% confidence level.

This paper suffers from conceptual difficulties and unjustified approximations that render its conclusions dubious.

One of the two breakout reactions from the hot CNOcycle is 15O(a,g)19Ne, which at low temperatures depends strongly on the resonance strength of the 4.033 MeV state in 19Ne. An experimental upper limit has been placed on its strength, but the lower limit on the resonance strength and thereby the astrophysical reaction rate is unconstrained experimentally. However, this breakout reaction is crucial to the thermonuclear runaway which causes type I X-ray bursts on accreting neutron stars. In this paper we exploit astronomical observations in an attempt to constrain the relevant nuclear physics and deduce a lower limit on the reaction rate. Our sensitivity study implies that if the rate were sufficiently small, accreting material would burn stably without bursts. The existence of type I X-ray bursts and superbursts consequently suggests a lower limit on the 15O(a,g)19Ne reaction rate at low temperatures.

The experimental landscape for the 7Be+p radiative capture reaction is rapidly changing as new high precision data become available. We present an evaluation of existing data, detailing the treatment of systematic errors and discrepancies, and show how they constrain the astrophysical S factor (S17), independent of any nuclear structure model. With theoretical models robustly determining the behavior of the sub-threshold pole, the extrapolation error can be reduced and a constraint placed on the slope of S17. Using only radiative capture data, we find S17(0) = 20.7 +/- 0.6 (stat) +/- 1.0 (syst) eV b if data sets are completely independent, while if data sets are completely correlated we find S17(0) = 21.4 +/- 0.5 (stat) +/- 1.4 (syst) eV b. The truth likely lies somewhere in between these two limits. Although we employ a formalism capable of treating discrepant data, we note that the central value of the S factor is dominated by the recent high precision data of Junghans et al., which imply a substantially higher value than other radiative capture and indirect measurements. Therefore we conclude that further progress will require new high precision data with a detailed error budget.

We present S factor data obtained from the Coulomb dissociation of 83 MeV/nucleon 8B, and analyze 7Be longitudinal momentum distributions measured at 44 and 81 MeV/nucleon using a potential model, first-order perturbation theory, and dynamical solution of the time-dependent Schroedinger equation. Comparing our results with independent continuum-discretized coupled channels calculations, we study the reaction model and beam energy dependence of the E2 contribution to the dissociation cross section. By fitting radiative capture and Coulomb breakup data taken below relative energies of 400 keV with potential models constrained by 7Li + n and 7Be + p elastic scattering data, we examine the mutual consistency of recent S17 measurements and obtain a recommended value for S17(0) of 18.6 +/- 0.4 (experimental) +/- 1.1 (extrapolation) eV b (1 sigma). This result is in good agreement with recent experimental determinations of the asymptotic normalization coefficient of the valence proton wave function in 8B.

Populating states in 22Mg via the (p,t) reaction in inverse kinematics with a 55 MeV/nucleon 24Mg beam, we have measured the proton-decay branching ratios of the levels at 5.96 MeV and 6.05 MeV and obtained an experimental upper limit on the branching ratio of the 5.71 MeV state. On the basis of the present and previous measurements, we assign spins and parities to the 5.96 MeV and 6.05 MeV states. We combine our branching ratios with independent measurements of the lifetimes of these states or their 22Ne analogs to compute the resonance strengths and thereby the astrophysical rate of the 21Na(p,g)22Mg reaction. We perform hydrodynamic calculations of nova outbursts with this new rate and analyze its impact on 22Na yields.

A recoil coincidence technique has been applied to measure the alpha-decay branching ratios of near-threshold states in 19Ne. Populating these states using the (p,t) reaction in inverse kinematics, we detected the recoils and their decay products with 100% geometric efficiency using a magnetic spectrometer. Combining our branching ratio measurements with independent determinations of the radiative widths of these states, we calculate the astrophysical rate of 15O(alpha,gamma)19Ne. Using this reaction rate, we perform hydrodynamic calculations of nova outbursts and conclude that no significant breakout from the hot CNO cycles into the rp process occurs in novae via 15O(alpha,gamma)19Ne.

The 15O(alpha,gamma)19Ne reaction is one of two known routes for breakout from the hot CNO cycles into the rp process. Its astrophysical rate depends on the decay properties of excited states in 19Ne lying just above the 15O + alpha threshold. We have measured the alpha-decay branching ratios for these states using the p(21Ne,t)19Ne reaction at 43 MeV/u. Combining our branching ratio measurements with previous determinations of the radiative widths of these states, we calculate the astrophysical rate of 15O(alpha,gamma)19Ne. Using this reaction rate, we perform hydrodynamic calculations of nova outbursts and conclude that no significant breakout from the hot CNO cycles into the rp process occurs in classical novae via 15O(alpha,gamma)19Ne.

The 15O(alpha,gamma)19Ne reaction is one of two routes for breakout from the hot CNO cycles into the rp process in accreting neutron stars. Its astrophysical rate depends critically on the decay properties of excited states in 19Ne lying just above the 15O + alpha threshold. We have measured the alpha-decay branching ratios for these states using the p(21Ne,t)19Ne reaction at 43 MeV/u. Combining our measurements with previous determinations of the radiative widths of these states, we conclude that no significant breakout from the hot CNO cycle into the rp process in novae is possible via 15O(alpha,gamma)19Ne, assuming current models accurately represent their temperature and density conditions.

In an effort to better determine the 7Be(p,gamma)8B reaction rate, we have performed inclusive and exclusive measurements of the Coulomb dissociation of 8B. The former was a study of longitudinal momentum distributions of 7Be fragments emitted in the Coulomb breakup of intermediate energy 8B beams on Pb and Ag targets. Analysis of these data yielded the E2 contribution to the breakup cross section. In the exclusive measurement, we determined the cross section for the Coulomb breakup of 8B on Pb at low relative energies in order to infer the astrophysical S factor for the 7Be(p,gamma)8B reaction. Interpreting the measurements with 1st-order perturbation theory, we obtained SE2/SE1 = 4.7 (+ 2.0,- 1.3) times 10^-4 at Erel = 0.6 MeV, and S17(0) = 17.8 (+ 1.4,- 1.2) eV b. Semiclassical 1st-order perturbation theory and fully quantum mechanical continuum-discretized coupled channels analyses yield nearly identical results for the E1 strength relevant to solar neutrino flux calculations, suggesting that theoretical reaction mechanism uncertainties need not limit the precision of Coulomb breakup determinations of the 7Be(p,gamma)8B S factor. A recommended value of S17(0) based on a weighted average of this and other measurements is presented.

A kinematically complete measurement was made of the Coulomb dissociation of 8B nuclei on a Pb target at 83 MeV/nucleon. The cross section was measured at low relative energies in order to infer the astrophysical S factor for the 7Be(p,gamma)8B reaction. A first-order perturbation theory analysis of the reaction dynamics including E1, E2, and M1 transitions was employed to extract the E1 strength relevant to neutrino-producing reactions in the solar interior. By fitting the measured cross section from Erel = 130 keV to 400 keV, we find S17(0) = 17.8 (+1.4, -1.2) eV b.