We show how the relativistic effective field theory for the superfluid phase of helium-4 can replace the standard methods used to compute the production rates of low momentum excitations due to the interaction with an external probe. This is done by studying the scattering problem of a light dark matter particle in the superfluid, and comparing to some existing results. We show that the rate of emission of two phonons, the Goldstone modes of the effective theory, gets strongly suppressed for sub-MeV dark matter particles due to a fine cancellation between two different tree-level diagrams in the limit of small exchanged momenta. This phenomenon is found to be a consequence of the particular choice of the potential felt by the dark matter particle in helium. The predicted rates can vary by orders of magnitude if this potential is changed. We prove that the dominant contribution to the total emission rate is provided by the phonons. Finally, we analyze the angular distributions for the emissions of one and two phonons, and discuss how they can be used to measure the mass of the hypothetical dark matter particle hitting the helium target.

We study the response of a He-4 detector to the interaction of sub-GeV dark matter using an effective field theory for the superfluid. We compute the lifetime of the phonon, which agrees with what known from standard techniques, hence providing an important check of the effective field theory. We then study the process of emission of two phonons, and show how its rate is much more suppressed than the phase space expectations; this is a consequence of the conservation of the current associated to the superfluid symmetries.

In this paper, we verify the large scale structure consistency relations using N-body simulations, including modes in the highly non-linear regime. These relations (pointed out by Kehagias & Riotto and Peloso & Pietroni) follow from the symmetry of the dynamics under a shift of the Newtonian potential by a constant and a linear gradient, and predict the absence of certain poles in the ratio between the (equal time) squeezed bispectrum and power spectrum. The consistency relations, as symmetry statements, are exact, but have not been previously checked beyond the perturbative regime. Our test using N-body simulations not only offers a non-perturbative check, but also serves as a warm-up exercise for applications to observational data. A number of subtleties arise when taking the squeezed limit of the bispectrum--we show how to circumvent or address them. An interesting by-product of our investigation is an explicit demonstration that the linear-gradient symmetry is unaffected by the periodic boundary condition of the simulations. Lastly, we verify using simulations that the consistency relations are violated when the initial conditions are non-gaussian (of the local fNL type). The methodology developed here paves the way for constraining primordial non-gaussianity using large scale structure data, including (numerous) highly non-linear modes that are otherwise hard to interpret and utilize.