We study the process of dark matter particles scattering off $^{3,4}$He with nuclear wave functions computed using an ab initio many-body framework. We employ realistic nuclear interactions from chiral effective field theory at next-to-next-to-leading order (NNLO) and develop an ab initio scheme to compute a general set of different nuclear response functions. In particular, we then perform an accompanying uncertainty quantification on these quantities and study error propagation to physical observables. We find a rich structure of allowed nuclear responses with significant uncertainties for certain spin-dependent interactions. The approach and results that are presented in this Paper establish a new framework for nuclear structure calculations and uncertainty quantification in the context of direct and (certain) indirect searches for dark matter.

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%.

Spectra of the 9Be(e,e') reaction have been measured at the S-DALINAC at an electron energy E_0 = 73 MeV and scattering angles of 93{\deg} and 141{\deg} with high energy resolution up to excitation energies E_x = 8 MeV. The astrophysically relevant resonance parameters of the first excited 1/2+ state of 9Be have been extracted in a one-level approximation of R-matrix theory resulting in a resonance energy E_R = 1.748(6) MeV and width Gamma_R = 274(8) keV in good agreement with the latest 9Be(gamma,n) experiment but with considerably improved uncertainties. However, the reduced B(E1) transition strength deduced from an extrapolation of the (e,e') data to the photon point is a factor of two smaller. Implications of the new results for a possible production of 12C in neutron-rich astrophysical scenarios are discussed.