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 present the first ab initio calculations for p-shell single-Lambda hypernuclei. For the solution of the many-baryon problem, we develop two variants of the no-core shell model with explicit $\Lambda$ and $\Sigma^+$, $\Sigma^0$, $\Sigma^-$ hyperons including $\Lambda$-$\Sigma$ conversion, optionally supplemented by a similarity renormalization group transformation to accelerate model-space convergence. In addition to state-of-the-art chiral two- and three-nucleon interactions, we use leading-order chiral hyperon-nucleon interactions and a recent meson-exchange hyperon-nucleon interaction. We validate the approach for s-shell hypernuclei and apply it to p-shell hypernuclei, in particular to $^7_\Lambda$Li, $^9_\Lambda$Be and $^{13}_\Lambda$C. We show that the chiral hyperon-nucleon interactions provide ground-state and excitation energies that agree with experiment within the cutoff dependence. At the same time we demonstrate that hypernuclear spectroscopy provides tight constraints on the hyperon-nucleon interactions and we discuss the impact of induced hyperon-nucleon-nucleon interactions.