Detecting Light Dark Matter with Magnons Tanner Trickle,1, 2, 3 Zhengkang Zhang,1, 2, 3 and Kathryn M. Zurek3, 2, 1 1Department of Physics, University of California, Berkeley, CA 94720, USA 2Theoretical Physics Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 3Walter Burke Institute for Theoretical Physics, California Institute of Technology, Pasadena, CA 91125, USA Scattering of light dark matter with sub-eV energy deposition can be detected with collective excitations in condensed matter systems. When dark matter has spin-independent couplings to atoms or ions, it has been shown to efficiently excite phonons. Here we show that, if dark matter couples to the electron spin, magnon excitations in materials with magnetic dipole order offer a promising detection path. We derive general formulae for single magnon excitation rates from dark matter scattering, and demonstrate as a proof of principle the projected reach of a yttrium iron garnet target for several dark matter models with spin-dependent interactions. This highlights the complementarity of various collective excitations in probing different dark matter interactions. Introduction | Direct detection of dark matter weaker optical response, have excellent reach [25, 30, 31]. (DM) has undergone a dramatic expansion of scope in Similarly, collective excitations can arise from different recent years. Well-motivated theories where sub-GeV degrees of freedom, such as charge or spin, and some DM arises in a hidden sector/hidden valley, with new excitations may be advantageous over others for certain weakly or strongly coupled dynamics (see e.g. Refs. [1{ types of DM couplings. Therefore, in order to identify 8] for early examples), have given impetus to new ideas the broadest DM detection strategy, it is important to to search for light DM. Conventional nuclear recoils, consider collective excitations of all types. well-matched kinematically to search for weak-scale DM, From this perspective, previous proposals via phonon are not effective for light DM | once the DM mass excitations are aimed at probing SI responses. While drops below the target nucleus mass, the fraction of the they cover many simple DM models, including those with DM's kinetic energy that can be deposited on the target a dark photon or scalar mediator, there are other scenar- falls. Beyond nuclear recoils, better DM-target kinematic ios that are equally plausible, where the leading DM-SM matching allows us to probe qualitatively new parame- interactions lead to stronger SD responses. For exam- ter space, through lighter targets (e.g. electrons) with ple, in dark photon mediated models, the DM may in ∼ eV (as in semiconductors and atoms [9{18] as well as fact be charge neutral, but couple to the dark photon molecules [19{21]) or ∼ meV (as in superconductors [22{ via a higher multipole, e.g. magnetic dipole or anapole 24] and Dirac materials [25]) energy gaps. Reading out [33, 35{45]. Also, a spin-0 mediator may dominantly cou- such small energy depositions is achieved through im- ple to the pseudoscalar (rather than scalar) current of provements to cryogenic superconducting calorimeters, SM fermions. In these scenarios, summarized in TableI, such as transition edge sensors (TES) and microwave SI responses are suppressed compared to the previously kinetic inductance devices (MKIDs). Collective excita- considered cases, and ideas of detecting SD responses are tions, such as phonons in superfluid helium [26{29] and needed. More generally, SI and SD couplings can co- crystals [30, 31], open new avenues for good kinematic exist, so it is desirable to pursue detection channels for matching. For example, the presence of O(10-100) meV both in order to have a more complete picture of DM gapped optical phonons in some systems facilitates the interactions. extraction of a large fraction of DM's kinetic energy for In this Letter, we propose a novel detection path DM as light as ∼10 keV. for spin-dependent light DM-electron interactions via arXiv:1905.13744v2 [hep-ph] 15 May 2020 Beyond kinematics, there is also a dynamics aspect magnon excitations. Magnons are quanta of collective of the problem | depending on how the DM couples spin wave excitations in condensed matter systems that to Standard Model (SM) particles, different target re- exhibit magnetic dipole order in the ground state. They sponses are relevant. A familiar example from nuclear can be thought of as the SD counterpart of phonons for recoils is the presence of several nuclear responses { spin- DM detection with similar kinematics. We demonstrate independent (SI), spin-dependent (SD), etc. { which can as a proof of principle that single magnon excitations can probe different DM-SM interactions [32{35]. Together probe interesting DM scenarios through scattering, thus they provide broad coverage of the DM theory space, broadening the coverage of the DM theory space. In fu- with various target nuclei offering complementary infor- ture work we will pursue DM (in particular axion DM) mation. Another example is dark photon mediated DM: absorption through magnon excitations. a material with a strong optical response, such as a super- conductor, has weak reach since the effective coupling of Magnons in magnetically ordered materials | Mag- the dark photon is suppressed due to in-medium effects, netic order can arise in solid state systems due to the while Dirac materials and polar crystals, which have interplay between electron-electron interactions, electron 2 α β 2 2 2 2 gχ µν µ ^α 4gχge αβ q q ^β gχge 6mχ+me Magnetic dipole DM L = χσ¯ χ Vµν + geeγ¯ e Vµ Oχ = δ − 2 Sχ σ¯e = 2 2 Λχ Λχme q π Λχ(mχ+me) 2 2 2 2 gχ µ 5 ν µ ^α 2gχge αβγ β ^γ gχge 3α µχe Anapole DM L = 2 χγ¯ γ χ @ Vµν + geeγ¯ e Vµ Oχ = 2 iq Sχ σ¯e = 4 Λχ Λχme π 2Λχ 2 2 2 5 ^α gχge α gχge µχe Pseudo-mediated DM L = gχχχφ¯ + gee¯ iγ e φ Oχ = − 2 iq 1χ σ¯e = 2 4 q me 4π α me TABLE I. Dark matter models, having Lagrangian L, with SD interactions considered in this work; these models are particularly well-motivated when DM does not carry a charge of any type, see e.g. Refs. [33, 35{45]. χ is a spin-1/2 DM particle, and V , φ are ultralight (typically eV) spin-1, spin-0 mediators, respectively. gχ, ge are dimensionless couplings, and Λχ is the ^α effective theory cutoff. In the nonrelativistic limit, these Lagrangians reduce to the operators Oχ (with Cartesian coordinates ^α α α = 1; 2; 3), as in Eq. (4). q ≡ jqj is the momentum transfer, and Sχ = σ =2 is the DM spin operator. σe is the reference cross section defined in Eq. (11) that we will use to present the reach. kinetic energy and Pauli exclusion (see e.g. Refs. [46, 47]). the nonrelativistic effective Lagrangian takes the form Such systems are usually described by a spin lattice 3 model, e.g. the Heisenberg model, X ^α ^α L = − Oχ (q)Se ; (4) N n α=1 1 X X H = J 0 0 S · S 0 0 : (1) ll jj lj l j where α denotes the Cartesian coordinates, and q is the 2 0 0 l;l =1 j;j =1 momentum transfer from the DM to the target. The op- ^ Here l; l0 label the magnetic unit cells, and j; j0 label the erators Oχ that follow from the three Lagrangians we magnetic atoms/ions inside the unit cell. Depending on consider are listed in TableI. Focusing on transitions from the ground state to single magnon states jν; ki, we the sign of the exchange coupling Jll0jj0 , the spins Slj obtain the matrix element as (see the Supplemental Ma- and Sl0j0 tend to align or anti-align. The low energy ex- citations are obtained by applying the Holstein-Primakoff terial for details) transformation to expand the spins around the ordered 3 ground state in terms of bosonic creation and annihila- sisf 1 X α α M (q) = δq;k+G p hsf jO^ (q)jsii ; (5) y ν;k χ ν;k;G tion operatorsa ^ ; a^. The quadratic part of the Hamilto- NΩ α=1 nian can then be diagonalized via a Bogoliubov transfor- mation (see the Supplemental Material for details), where Ω is the volume of the magnetic unit cell, G de- notes a reciprocal lattice vector, and jsi;f i are the initial ! ! ! a^ U V ^b and final DM spin states. ν;k;G is the analog of polar- j;k = jν;k jν;k ν;k ; (2) y ∗ ∗ ^y ization vectors for the magnon modes, a^j;−k Vjν;−k Ujν;−k bν;−k n r n S X j ∗ ∗ iG·xj X X ^y ^ ν;k;G = Vjν;−krj + Ujν;krj e ; (6) H = !ν;kb bν;k ; (3) 2 ν;k j=1 ν=1 k21BZ where rα ≡ Rα1+iRα2 parameterize the spin orientations ^y ^ j j j so that b ; b are creation and annihilation operators of the in the ground state, canonical magnon modes, which are collective excitations of the spins. For a system with N magnetic unit cells α X αβ 0β 01 02 03 Slj = Rj Slj ; fhSlj i; hSlj i; hSlj ig = f0; 0;Sjg ; and n magnetic atoms/ions in the unit cell, there are n β magnon branches, labeled by ν, with N modes on each (7) branch, labeled by momentum vectors k within the first and xj ≡ xlj − xl is the position of the jth site within a (magnetic) Brillouin zone (1BZ). The n × n matrices U, magnetic unit cell. As a simple example, a ferromagnet V can be calculated for each k.