A New Evaluation of the 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{11}$$\end{document}B(p,α)αα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha )\alpha \alpha$$\end{document} Reaction Rates

Reaction rates of the 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{11}$$\end{document}B(p,α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document})αα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \alpha$$\end{document} process have been evaluated on the basis of a data set spanning incident proton energies Ep\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_p$$\end{document} from 0.15 to 3.8 MeV. A previously published analysis (Spraker et al. in J Fusion Energy 31(4):357, 2012) of these data provided the number of outgoing α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-particles in a restricted range of the detected α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-energy spectrum, making it unsuitable for the evaluation of the reaction rates. The present work takes advantage of a calculation of the α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-energy spectrum based on a sequential model of the reaction and the assumption that the primary α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-particles are emitted with ℓ=3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\ell =3$$\end{document}. A full description of this ansatz, which has been shown to reproduce the essential features of the observed α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-energy spectra, can be found in Stave et al. (Phys Lett B 696:26, 2011). The accuracy of these calculated spectra has made it possible to reliably extrapolate the new data to zero-energy α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-particles. In the ensuing calculation of the cross section, the total measured α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-yield is then divided by a fixed factor of three at all incident proton energies. In addition, this technique has enabled a treatment of the α0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha _0$$\end{document} channel where the 12C nuclei decay to the ground state of 8Be via emission of an α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}-particle. This channel contributes at incident proton energies above 2 MeV. The new cross section data have then been used to evaluate the 11\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{11}$$\end{document}B(p,α\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document})αα\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha \alpha$$\end{document} reaction rates. The new evaluation is ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}10–15 % higher than the currently accepted result (Angulo et al. in Nucl Phys A 656(1):3, 1999) at temperatures between 200 and 600 keV (2–7 ×\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document} 109\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^9$$\end{document} K). The inclusion of a narrow, low-lying resonance at Ep=0.162\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_p=0.162$$\end{document} MeV in the evaluation is found to have a minimal effect on the reaction rate above 100 keV (1.2 ×\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document} 109\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^9$$\end{document} K), and a higher-lying state at Ep=3.75\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_p=3.75$$\end{document} MeV is shown to enhance the reaction rates by only ∼\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sim$$\end{document}15 % above 400 keV (4.6 ×\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\times$$\end{document} 109\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^9$$\end{document} K).


Introduction
The 11 B(p,a)aa reaction at incident proton energies E p \4 MeV has been studied since the 1930s [4][5][6][7][8]. The reaction proceeds through a two-step process involving emission of a primary a-particle followed by two secondary a-particles emitted from the decay of 8 Be. Based on cross section data collected previously at the Triangle Universities Nuclear Laboratory (TUNL) [1] and a model of the reaction mechanism developed to account for the observed a-energy spectra [2], the reaction can be described via its behavior at resonances at E p ¼ 0:675 and E p ¼ 2:64 MeV. At E p ¼ 0:675 MeV, the primary a-particles produced from the decay of the 2 À state in 12 C have ' ¼ 3, forming the 2 þ first-excited state of 8 Be. This then decays to two a-particles with ' 0 ¼ 2. The 3 À state in 12 C formed at E p ¼ 2:64 MeV has two possible decay modes. Primary aparticles can be emitted with either ' ¼ 1 to the first-excited state of 8 Be, or with ' ¼ 3 to the 0 þ ground state. The residual 8 Be nucleus then decays, emitting two secondary a-particles.
Interest in the 11 B(p,a)aa reaction is driven by studies in stellar evolution, where relative abundances of 11

B, Li and
Be provide insight into stellar processes [9]. This reaction is also being investigated for use in practical nuclear fusion reactors, exploiting the aneutronic nature of this fusion reaction [10]. Total cross sections and reaction rates are key input parameters in the design of such reactors.
The full observed a-energy spectra must be integrated over all outgoing energies and the yield divided by 3 to account for the final-state multiplicity in order to extract absolute cross sections. Practically, this cannot be done using only the collected data for two reasons. First, the aenergy spectra contain a dominant peak at the nominal beam energy from protons elastically scattered from the target. The other reason is the low-energy truncation of the energy spectra due to hardware thresholds required to reduce background rates which would overwhelm the data acquisition system. Consequently, the analysis presented in [1] calculated yields by integrating over a 2.75 MeV wide window centered on the dominant a-particle peak observed in each detector, circumventing the ambiguity in the number of a-particles located in the integration region by reporting results in units of Counts/Luminosity ( X), defined as where N t is the target density, N p is the number of incident protons, and dX is the detector solid angle. To convert X to a differential cross section, the number of a-particles in the energy window used must be divided out. This number varies with proton energy and is difficult to determine. The present work addresses this complication by utilizing the model described in Ref. [2] to extrapolate the entire a-energy spectra. Yields are then obtained from these extrapolations to extract angular distributions and total cross sections, which are then used to evaluate the reaction rate.

Total Cross Sections
The aforementioned model was used to generate simulated a-energy spectra at each measurement angle and beam energy over 0.15 MeV \E p \ 3.8 MeV. The reaction was taken to proceed through excitation of the 2 À state in 12 C for E p 2.0 MeV, while the 3 À state was used for E p [ 2.0 MeV. Each simulated spectrum was normalized to the corresponding measured a-energy spectrum using the previously discussed 2.75 MeV wide summing region. The number of counts observed in the data over the integration region was compared to the simulated yield. This ratio defined a scaling factor that was then applied to the total simulated a-energy spectrum. The process is illustrated in Fig. 1 at E p ¼ 0:65 MeV at h lab a ¼ 90 . By definition, the scaled simulated a-energy spectrum has the same number of counts in the original summing region as the data. The final total yield was then obtained by integrating the scaled simulated a-energy spectrum down to 0.
As noted in Sect. 1, the two decay modes of the 3 À state in 12 C can lead to contributions from both the a 0 and a 1 channels from decays to the ground state or the first excited state of 8 Be, respectively. For incident proton energies exciting this state, the same procedure was used to obtain The total a-contribution calculated using the model was normalized to the data over the integration window used in the previous analysis, indicated above. The a-energy spectrum was then extended to E a ¼ 0 by using the normalized model output  Comparison of observed and calculated a-energy spectra at h lab a ¼ 90 and E p ¼ 2:6 MeV. The total a 1 -contribution calculated using the model was normalized to the data over the integration window used in the previous analysis, indicated above. The a-energy spectrum was then extended to E a ¼ 0 by using the normalized model output J Fusion Energ (2016) 35:538-543 539 the a 1 yields as described for the resonance at E p ¼ 0:675 MeV, shown in Fig. 2. The same technique was then applied separately to the a 0 channel, as demonstrated in Fig. 3. The primary a 0 peak is clearly separated from the a 1 contribution, providing a convenient reference for normalizing the modeled primary a 0 yield. The scaling factor obtained in this way is then applied to the calculated yield of secondary a 0 particles, which have an energy distribution lying below the proton elastic scattering peak. As in [1], the extrapolated a 0 and a 1 spectra were then integrated down to E a ¼ 0, and the resulting Counts/Luminosity were used to obtain angular distributions which were fit to an expansion in terms of Legendre polynomials The data obtained at E p \1:4 MeV were fit with a Legendre expansion up to and including i ¼ 2, while for E p ! 1:4 MeV the expansion included up to i ¼ 4. Similarly, the a 0 angular distributions went up to i ¼ 4. More than 95 % of the individual data points were within 3 % of the corresponding fit result. The expansion coefficients are listed in Tables 1, 2 and 3 in the Appendix. The A 0 coefficients extracted from the fits were then multiplied by 4p=3 (with the factor of 3 accounting for the total number of final-state a-particles) to give total cross sections. The results, which include the a 0 channel where applicable, are displayed as a function of E p in Fig. 4.
A 2 þ state in 12 C with a full width at half maximum of C ¼ 5 keV is excited at E p ¼ 0:162 MeV with a maximum cross section of 100 mb. Measuring this state experimentally was not possible due to the narrow width relative to the energy resolution of the proton beam. In order to account for the effect of this state on the reaction rates, cross section values were generated using an existing parameterization of the astrophysical factor S(E) based on a fit to experimental data [7,11]. The cross section as a function of center of mass energy E could then be obtained using

Reaction Rates
Using the total cross section values from Sect. 2, the Maxwellian-averaged reaction rates [3] were computed as a function of temperature using where l is the reduced mass, k B is the Boltzmann constant, T is temperature, r is the total cross section, m is the relative velocity, and E is the center of mass energy. The reaction rate as a function of temperature is shown in Fig. 5. The dot-dashed lines were obtained by varying the total cross section by the combined statistical ( $ 0.1 %) and systematic (3.4 %) uncertainties. Also shown is a previous evaluation of the reaction rate [3].
The previous result was based on an interpretation of the reaction described in [7] which assumed that primary aparticles emitted from the decay of the 2 À resonance at E p ¼ 0:675 MeV have ' ¼ 1 rather than ' ¼ 3. In this ansatz, the simulated kinematics of the secondary a-particles indicated that the observed a 1 yields should be divided by 2 based on the chosen integration region. This interpretation, adapted in [3] in evaluating the reaction rate displayed in Fig. 5, leads to rates that are 10-15 % lower than the new results over a temperature range of 200-600 keV. As shown in [2], this ansatz does not give an accurate description of the new data set.
To determine the importance of including the narrow resonance at E p ¼ 0:162 MeV, the reaction rates were evaluated both with and without inclusion of the cross section values calculated with Eq. 3. As the comparison shown in Fig. 6 demonstrates, inclusion of these generated values more than doubles the reaction rates at temperatures near 20 keV. At temperatures greater than 100 keV, the enhancement becomes negligible, as the reaction mechanism is dominated by the much stronger resonances at E p ¼ 0:675 and E p ¼ 2:64 MeV.
The current analysis has only considered the resonances occuring at E p ¼ 0:675 MeV and E p ¼ 2:64 MeV in modeling the energy and angular distributions of the outgoing a-particles. The range of incident proton energies extends high enough to excite a higher-lying state in 12 C at E p ¼ 3:75 MeV (C ¼ 480 keV) with a possible spin assignment of J p ¼ 2 À . The effect on the reaction rate of neglecting the contribution from this state can be assesed by truncating the cross section data set at various values of E p , as shown in Fig. 7. The effect is minimal below 400 keV and and is less than $ 15 % at higher temperatures, indicating that the two resonances included in the simulation provide a reasonable description of the data even at higher energies.  Conclusion A data set built on previously published cross section data of the 11 B(p,a)aa reaction [1] and extended with a detailed model of the reaction mechanism [2] has been used to generate new angular distributions and total cross section results. The new results cover the entire a-energy spectra and are now reported as true cross sections rather than in terms of Counts/Luminosity. These new results provide the basis for a new evaluation of the reaction rates, which are 10-15 % higher than a previous evaluation [3] at temperatures between 200 and 600 keV (2-7 Â 10 9 K). These results will provide useful inputs to astrophysical studies and will have application in the design of novel aneutronic fusion reactors.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.