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.

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.

We calculate the energy per particle of symmetric nuclear matter and pure neutron matter using the microscopic many-body Brueckner-Hartree-Fock (BHF) approach and employing the Argonne V18 (AV18) nucleon-nucleon (NN) potential supplemented with two different three-nucleon force models recently constructed to reproduce the binding energy of $^3$H, $^3$He and $^4$He nuclei as well as the neutron-deuteron doublet scattering length. We find that none of these new three-nucleon force models is able to reproduce simultaneously the empirical saturation point of symmetric nuclear matter and the properties of three- and four-nucleon systems.

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.