We construct an effective field theory (EFT) description of the hard photon spectrum for heavy WIMP annihilation. This facilitates precision predictions relevant for line searches, and allows the incorporation of non-trivial energy resolution effects. Our framework combines techniques from non-relativistic EFTs and soft-collinear effective theory (SCET), as well as its multi-scale extensions that have been recently introduced for studying jet substructure. We find a number of interesting features, including the simultaneous presence of SCET$_{\text{I}}$ and SCET$_{\text{II}}$ modes, as well as collinear-soft modes at the electroweak scale. We derive a factorization formula that enables both the resummation of the leading large Sudakov double logarithms that appear in the perturbative spectrum, and the inclusion of Sommerfeld enhancement effects. Consistency of this factorization is demonstrated to leading logarithmic order through explicit calculation. Our final result contains both the exclusive and the inclusive limits, thereby providing a unifying description of these two previously-considered approximations. We estimate the impact on experimental sensitivity, focusing for concreteness on an SU(2)$_{W}$ triplet fermion dark matter - the pure wino - where the strongest constraints are due to a search for gamma-ray lines from the Galactic Center. We find numerically significant corrections compared to previous results, thereby highlighting the importance of accounting for the photon spectrum when interpreting data from current and future indirect detection experiments.

WIMP-nucleon scattering is analyzed at order $1/M$ in Heavy WIMP Effective Theory. The $1/M$ power corrections, where $M\gg m_W$ is the WIMP mass, distinguish between different underlying UV models with the same universal limit and their impact on direct detection rates can be enhanced relative to naive expectations due to generic amplitude-level cancellations at leading order. The necessary one- and two-loop matching calculations onto the low-energy effective theory for WIMP interactions with Standard Model quarks and gluons are performed for the case of an electroweak SU(2) triplet WIMP, considering both the cases of elementary fermions and composite scalars. The low-velocity WIMP-nucleon scattering cross section is evaluated and compared with current experimental limits and projected future sensitivities. Our results provide the most robust prediction for electroweak triplet Majorana fermion dark matter direct detection rates; for this case, a cancellation between two sources of power corrections yields a small total $1/M$ correction, and a total cross section close to the universal limit for $M \gtrsim {\rm few} \times 100\,{\rm GeV}$. For the SU(2) composite scalar, the $1/M$ corrections introduce dependence on underlying strong dynamics. Using a leading chiral logarithm evaluation, the total $1/M$ correction has a larger magnitude and uncertainty than in the fermionic case, with a sign that further suppresses the total cross section. These examples provide definite targets for future direct detection experiments and motivate large scale detectors capable of probing to the neutrino floor in the TeV mass regime.

Primordial non-Gaussianities enhanced at small wavevectors can induce a power spectrum of the galaxy overdensity that differs greatly from that of the matter overdensity at large length scales. In previous work, it was shown that "squeezed" three-point and "collapsed" four-point functions of the curvature perturbation $\zeta$ can generate these non-Gaussianities and give rise to so-called scale-dependent and stochastic bias in the galaxy overdensity power spectrum. We explore a third way to generate non-Gaussianities enhanced at small wavevectors: the infrared behavior of quantum loop contributions to the four-point correlations of $\zeta$. We show that these loop effects lead to stochastic bias, which can be observable in the context of quasi-single field inflation.

We apply an orthogonalization procedure on the effective field theory of large scale structure (EFT of LSS) shapes, relevant for the angle-averaged bispectrum and non-Gaussian covariance of the matter power spectrum at one loop. Assuming natural-sized EFT parameters, this identifies a linear combination of EFT shapes - referred to as the principal shape - that gives the dominant contribution for the whole kinematic plane, with subdominant combinations suppressed by a few orders of magnitude. For the covariance, our orthogonal transformation is in excellent agreement with a principal component analysis applied to available data. Additionally we find that, for both observables, the coefficients of the principal shapes are well approximated by the EFT coefficients appearing in the squeezed limit, and are thus measurable from power spectrum response functions. Employing data from N-body simulations for the growth-only response, we measure the single EFT coefficient describing the angle-averaged bispectrum with $\mathcal{O}(10\%)$ precision. These methods of shape orthogonalization and measurement of coefficients from response functions are valuable tools for developing the EFT of LSS framework, and can be applied to more general observables.

We apply effective field theory methods to compute bino-nucleon scattering, in the case where tree-level interactions are suppressed and the leading contribution is at loop order via heavy flavor squarks or sleptons. We find that leading log corrections to fixed-order calculations can increase the bino mass reach of direct detection experiments by a factor of two in some models. These effects are particularly large for the bino-sbottom coannihilation region, where bino dark matter as heavy as 5-10 TeV may be detected by near future experiments. For the case of stop- and selectron-loop mediated scattering, an experiment reaching the neutrino background will probe thermal binos as heavy as 500 and 300 GeV, respectively. We present three key examples that illustrate in detail the framework for determining weak scale coefficients, and for mapping onto a low energy theory at hadronic scales, through a sequence of effective theories and renormalization group evolution. For the case of a squark degenerate with the bino, we extend the framework to include a squark degree of freedom at low energies using heavy particle effective theory, thus accounting for large logarithms through a "heavy-light current." Benchmark predictions for scattering cross sections are evaluated, including complete leading order matching onto quark and gluon operators, and a systematic treatment of perturbative and hadronic uncertainties.

We compute the non-Gaussian contribution to the covariance of the matter power spectrum at one-loop order in Standard Perturbation Theory (SPT), and using the framework of the effective field theory (EFT) of large scale structure (LSS). The complete one-loop contributions are evaluated for the first time, including the leading EFT corrections that involve seven independent operators, of which four appear in the power spectrum and bispectrum. We compare the non-Gaussian part of the one-loop covariance computed with both SPT and EFT of LSS to two separate simulations. In one simulation, we find that the one-loop prediction from SPT reproduces the simulation well to $k_i + k_j \sim$ 0.25 h/Mpc, while in the other simulation we find a substantial improvement of EFT of LSS (with one free parameter) over SPT, more than doubling the range of $k$ where the theory accurately reproduces the simulation. The disagreement between these two simulations points to unaccounted for systematics, highlighting the need for improved numerical and analytic understanding of the covariance.

We compute the connected four point correlation function (the trispectrum in Fourier space) of cosmological density perturbations at one-loop order in Standard Perturbation Theory (SPT) and the Effective Field Theory of Large Scale Structure (EFT of LSS). This paper is a companion to our earlier work on the non-Gaussian covariance of the matter power spectrum, which corresponds to a particular wavenumber configuration of the trispectrum. In the present calculation, we highlight and clarify some of the subtle aspects of the EFT framework that arise at third order in perturbation theory for general wavenumber configurations of the trispectrum. We consistently incorporate vorticity and non-locality in time into the EFT counterterms and lay out a complete basis of building blocks for the stress tensor. We show predictions for the one-loop SPT trispectrum and the EFT contributions, focusing on configurations which have particular relevance for using LSS to constrain primordial non-Gaussianity.

The discovery of a Standard Model-like Higgs boson and the hitherto absence of evidence for other new states may indicate that if WIMPs comprise cosmological dark matter, they are heavy compared to electroweak scale particles, $M \gg m_{W^\pm}, m_{Z^0}$. In this limit, the absolute cross section for a WIMP of given electroweak quantum numbers to scatter from a nucleon becomes computable in terms of Standard Model parameters. We develop effective field theory techniques to analyze the heavy WIMP limit of WIMP-nucleon scattering, and present the first complete calculation of the leading spin-independent cross section in Standard Model extensions consisting of one or two electroweak $SU(2)_W \times U(1)_Y$ multiplets.

Particles that are heavy compared to the electroweak scale ($M \gg m_W$), and that are charged under electroweak SU(2) gauge interactions display universal properties such as a characteristic fine structure in the mass spectrum induced by electroweak symmetry breaking, and an approximately universal cross section for scattering on nuclear targets. The heavy particle effective theory framework is developed to compute these properties. As illustration, the spin independent cross section for low-velocity scattering on a nucleon is evaluated in the limit $M \gg m_W$, including complete leading-order matching onto quark and gluon operators, renormalization analysis, and systematic treatment of perturbative and hadronic-input uncertainties.