We calculate the $p$-wave contribution to the proton-proton fusion $S$ factor and its energy derivative in pionless effective field theory (EFT) up to next-to-leading order. The leading contributions are given by a recoil piece from the Gamow-Teller and Fermi operators, and from relativistic $1/m$ suppressed weak interaction operators. We obtain the value of $(2.5\pm0.3 )\times 10^{-28}~\mathrm{MeV\ fm^2}$ for the $S$ factor and $(2.2\pm0.2) \times 10^{-26}~\mathrm{fm^2}$ for its energy derivative at threshold. These are smaller than the results of a prior study that employed chiral EFT by several orders of magnitude. We conclude that, contrary to what has been previously reported, the $p$-wave contribution does not need to be considered in a high-precision determination of the $S$ factor at astrophysical energies. Combined with the chiral EFT calculation of Acharya {\it et al.} [Phys. Lett. B \bf{760}, 584 (2016)] for the $s$-wave channel, this gives a total threshold $S$ factor of $S(0) = (4.047^{+0.024}_{-0.032}) \times 10^{-23}~{\rm MeV~fm}^2$.

We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($\chi$EFT) up to next-to-next-to-leading order (NNLO) and perform a rigorous uncertainty analysis of the results. We quantify the uncertainties due to (i) the computational method used to compute the $pp$ cross section in momentum space, (ii) the statistical uncertainties in the low-energy coupling constants of $\chi$EFT, (iii) the systematic uncertainty due to the $\chi$EFT cutoff, and (iv) systematic variations in the database used to calibrate the nucleon-nucleon interaction. We also examine the robustness of the polynomial extrapolation procedure, which is commonly used to extract the threshold $S$-factor and its energy-derivatives. By performing a statistical analysis of the polynomial fit of the energy-dependent $S$-factor at several different energy intervals, we eliminate a systematic uncertainty that can arise from the choice of the fit interval in our calculations. In addition, we explore the statistical correlations between the $S$-factor and few-nucleon observables such as the binding energies and point-proton radii of $^{2,3}$H and $^3$He as well as the $D$-state probability and quadrupole moment of $^2$H, and the $\beta$-decay of $^{3}$H. We find that, with the state-of-the-art optimization of the nuclear Hamiltonian, the statistical uncertainty in the threshold $S$-factor cannot be reduced beyond 0.7%.

We study the impact on the primordial abundances of light elements created by a variation of the quark masses at the time of Big Bang nucleosynthesis (BBN). In order to navigate through the particle and nuclear physics required to connect quark masses to binding energies and reaction rates in a model-independent way, we use lattice QCD data and a hierarchy of effective field theories. We find that the measured Helium-4 abundances put a bound of -1 % <~ d m_q/m_q <~ 0.7 % on a possible variation of quark masses. The effect of quark mass variations on the deuterium abundances can be largely compensated by changes of the baryon-to-photon ratio eta. Including bounds on the variation of eta coming from WMAP results and adding some additional assumptions further narrows the range of allowed values of d m_q/m_q.