Among the reactions involved in the production and destruction of deuterium during Big Bang Nucleosynthesis, the deuterium-burning D(p,gamma)3He reaction has the largest uncertainty and limits the precision of theoretical estimates of primordial deuterium abundance. Here we report the results of a careful commissioning of the experimental setup used to measure the cross-section of the D(p,gamma)3He reaction at the Laboratory for Underground Nuclear Astrophysics of the Gran Sasso Laboratory (Italy). The commissioning was aimed at minimising all sources of systematic uncertainty in the measured cross sections. The overall systematic error achieved (< 3 %) will enable improved predictions of BBN deuterium abundance.

The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.

We present a detailed description of the electromagnetic filter for the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Starting with an initial estimate for the orbital magnetic moment, the higher-order drift process of ExB is configured to balance the gradient-B drift motion of the electron in such a way as to guide the trajectory into the standing voltage potential along the mid-plane of the filter. As a function of drift distance along the length of the filter, the filter zooms in with exponentially increasing precision on the transverse velocity component of the electron kinetic energy. This yields a linear dimension for the total filter length that is exceptionally compact compared to previous techniques for electromagnetic filtering. The parallel velocity component of the electron kinetic energy oscillates in an electrostatic harmonic trap as the electron drifts along the length of the filter. An analysis of the phase-space volume conservation validates the expected behavior of the filter from the adiabatic invariance of the orbital magnetic moment and energy conservation following Liouville's theorem for Hamiltonian systems.

We propose to achieve the proof-of-principle of the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Each of the technological challenges described in [1,2] will be targeted and hopefully solved by the use of the latest experimental developments and profiting from the low background environment provided by the LNGS underground site. The first phase will focus on the graphene technology for a tritium target and the demonstration of TES microcalorimetry with an energy resolution of better than 0.05 eV for low energy electrons. These technologies will be evaluated using the PTOLEMY prototype, proposed for underground installation, using precision HV controls to step down the kinematic energy of endpoint electrons to match the calorimeter dynamic range and rate capabilities. The second phase will produce a novel implementation of the EM filter that is scalable to the full target size and which demonstrates intrinsic triggering capability for selecting endpoint electrons. Concurrent with the CNB program, we plan to exploit and develop the unique properties of graphene to implement an intermediate program for direct directional detection of MeV dark matter [3,4]. This program will evaluate the radio-purity and scalability of the graphene fabrication process with the goal of using recently identified ultra-high radio-purity CO2 sources. The direct detection of the CNB is a snapshot of early universe dynamics recorded by the thermal relic neutrino yield taken at a time that predates the epochs of Big Bang Nucleosynthesis, the Cosmic Microwave Background and the recession of galaxies (Hubble Expansion). Big Bang neutrinos are believed to have a central role in the evolution of the Universe and a direct measurement with PTOLEMY will unequivocally establish the extent to which these predictions match present-day neutrino densities.

The experimental study of nuclear reactions of astrophysical interest is greatly facilitated by a low-background, high-luminosity setup. The Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator offers ultra-low cosmic-ray induced background due to its location deep underground in the Gran Sasso National Laboratory (INFN-LNGS), Italy, and high intensity, 250-500 $\mu$A, proton and $\alpha$ ion beams. In order to fully exploit these features, a high-purity, recirculating gas target system for isotopically enriched gases is coupled to a high-efficiency, six-fold optically segmented bismuth germanate (BGO) $\gamma$-ray detector. The beam intensity is measured with a beam calorimeter with constant temperature gradient. Pressure and temperature measurements have been carried out at several positions along the beam path, and the resultant gas density profile has been determined. Calibrated $\gamma$-intensity standards and the well-known $E_p$ = 278 keV $\mathrm{^{14}N(p,\gamma)^{15}O}$ resonance were used to determine the $\gamma$-ray detection efficiency and to validate the simulation of the target and detector setup. As an example, the recently measured resonance at $E_p$ = 189.5 keV in the $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction has been investigated with high statistics, and the $\gamma$-decay branching ratios of the resonance have been determined.

The astrophysical $S$-factor for the radiative capture $d(p,\gamma)^3$He in the energy-range of interest for Big Bang Nucleosynthesis (BBN) is calculated using an {\it ab-initio} approach. The nuclear Hamiltonian retains both two- and three-nucleon interactions - the Argonne $v_{18}$ and the Urbana IX, respectively. Both one- and many-body contributions to the nuclear current operator are included. The former retain for the first time, besides the $1/m$ leading order contribution ($m$ is the nucleon mass), also the next-to-leading order term, proportional to $1/m^3$. The many-body currents are constructed in order to satisfy the current conservation relation with the adopted Hamiltonian model. The hyperspherical harmonics technique is applied to solve the $A=3$ bound and scattering states. A particular attention is used in this second case in order to obtain, in the energy range of BBN, an uncertainty on the astrophysical $S$-factor of the order or below $\sim$1 %. Then, in this energy range, the $S$-factor is found to be $\sim$10 % larger than the currently adopted values.Part of this increase (1-3 %) is due to the $1/m^3$ one-body operator, while the remaining is due to the new more accurate scattering wave functions. We have studied the implication of this new determination for the $d(p,\gamma)^3$He $S$-factor on deuterium primordial abundance. We find that the predicted theoretical value for $^2$H/H is in excellent agreement with its experimental determination, using the most recent determination of baryon density of Planck experiment, and with a standard number of relativistic degrees of freedom $N_{\rm eff}=3.046$ during primordial nucleosynthesis.

The p(p,e^+ \nu_e)^2H reaction rate is an essential ingredient for theoretical computations of stellar models. In the past several values of the corresponding S-factor have been made available by different authors. Prompted by a recent evaluation of S(E), we analysed the effect of the adoption of different proton-proton reaction rates on stellar models, focusing, in particular, on the age of mid and old stellar clusters (1-12 Gyr) and on standard solar model predictions. By comparing different widely adopted p(p,e^+ \nu_e)^2H reaction rates, we found a maximum difference in the temperature regimes typical of main sequence hydrogen-burning stars (5x10^6 - 3x10^7 K) of about 3%. Such a variation translates into a change of cluster age determination lower than 1%. A slightly larger effect is observed in the predicted solar neutrino fluxes with a maximum difference, in the worst case, of about 8%. Finally we also notice that the uncertainty evaluation of the present proton-proton rate is at the level of few \permil, thus the p(p,e^+ \nu_e)^2H reaction rate does not constitute anymore a significant uncertainty source in stellar models.

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.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the threshold $S$-factors for the solar $pp$ (proton-fusion) and $hep$ processes in an effective field theory that {\it combines} the merits of the standard nuclear physics method and systematic chiral expansion. The power of the EFT adopted here is that one can correlate in a unified formalism the weak-current matrix elements of two-, three- and four-nucleon systems. Using the tritium $\beta$-decay rate as an input to fix the only unknown parameter in the theory, we can evaluate the threshold $S$ factors with drastically improved precision; the results are $S_{pp}(0) = 3.94\times(1 \pm 0.004) \times 10^{-25} {MeV-b}$ and $S_{hep}(0) = (8.6\pm 1.3)\times 10^{-20} {keV-b}$. The dependence of the calculated $S$-factors on the momentum cutoff parameter $\Lambda$ has been examined for a physically reasonable range of $\Lambda$. This dependence is found to be extremely small for the $pp$ process, and to be within acceptable levels for the $hep$ process, substantiating the consistency of our calculational scheme.

Using effective field theory, we calculate the S-factor for the hep process in a totally parameter-free formulation. The transition operators are organized according to chiral counting, and their matrix elements are evaluated using the realistic nuclear wave functions obtained in the correlated-hyperspherical-harmonics method. Terms of up to next-to-next-to-next-to-leading order in heavy-baryon chiral perturbation theory are considered. Fixing the only parameter in the theory by fitting the tritium \beta-decay rate, we predict the hep S-factor with accuracy better than \sim 20 %.

Spurred by the recent complete determination of the weak currents in two-nucleon systems up to ${\cal O}(Q^3)$ in heavy-baryon chiral perturbation theory, we carry out a parameter-free calculation of the solar proton fusion rate in an effective field theory that combines the merits of the standard nuclear physics method and systematic chiral expansion. Using the tritium beta-decay rate as an input to fix the only unknown parameter in the effective Lagrangian, we can evaluate with drastically improved precision the ratio of the two-body contribution to the well established one-body contribution; the ratio is determined to be (0.86\pm 0.05) %. This result is essentially independent of the cutoff parameter for a wide range of its variation (500 MeV \le \Lambda \le 800 MeV), a feature that substantiates the consistency of the calculation.

The astrophysical S-factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to realistic Hamiltonians consisting of the Argonne v14 or Argonne v18 two-nucleon and Urbana-VIII or Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p 3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p 3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b with the AV18/UIX Hamiltonian, a factor of about 4.5 larger than the value adopted in the standard solar model. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor. The energy dependence is rather weak: the AV18/UIX zero-energy S-factor is 9.64 10^{-20} keV b, only 5 % smaller than the 10 keV result quoted above. The model dependence is also found to be weak: the zero-energy S-factor is calculated to be 10.2 10^{-20} keV b with the older AV14/UVIII model, only 6 % larger than the AV18/UIX result. Our best estimate for the S-factor at 10 keV is therefore (10.1 \pm 0.6) 10^{-20} keV b, when the theoretical uncertainty due to the model dependence is included. This value for the calculated S-factor is not as large as determined in fits to the Super-Kamiokande data in which the hep flux normalization is free. However, the precise calculation of the S-factor and the consequent absolute prediction for the hep neutrino flux will allow much greater discrimination among proposed solar neutrino oscillation solutions.

The astrophysical factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p-3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p-3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b, a factor of five larger than the standard-solar-model value. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor.