## Electron Ionization via Dark Matter-Electron Scattering and the Migdal Effect

Daniel Baxter, 1, 2, ∗ Yonatan Kahn, 1, 3, † and Gordan Krnjaic 4, ‡

1 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL USA

2 Enrico Fermi Institute, University of Chicago, Chicago, IL USA

3 University of Illinois at Urbana-Champaign, Urbana, IL USA

4 Fermi National Accelerator Laboratory, Batavia, IL USA

(Dated: November 24, 2020)

There are currently several existing and proposed experiments designed to probe sub-GeV dark matter (DM) using electron ionization in various materials. The projected signal rates for these experiments assume that this ionization yield arises only from DM scattering directly off electron targets, ignoring secondary ionization contributions from DM scattering off nuclear targets. We investigate the validity of this assumption and show that if sub-GeV DM couples with comparable strength to both protons and electrons, as would be the case for a dark photon mediator, the ionization signal from atomic scattering via the Migdal effect scales with the atomic number Z and 3-momentum transfer q as Z 2 q 2 . The result is that the Migdal effect is always subdominant to electron scattering when the mediator is light, but that Migdal-induced ionization can dominate over electron scattering for heavy mediators and DM masses in the hundreds of MeV range. We put these two ionization processes on identical theoretical footing, address some theoretical uncertainties in the choice of atomic wavefunctions used to compute rates, and discuss the implications for DM scenarios where the Migdal process dominates, including for XENON10, XENON100, and the recent XENON1T results on light DM scattering.

## I. INTRODUCTION

Although the evidence for dark matter (DM) is overwhelming, its microscopic properties remain unknown and motivate various experimental techniques to uncover its possible non-gravitational interactions [1]. In recent years, there have been several novel experimental techniques introduced to directly detect hitherto inaccessible DMcandidates below the GeV scale [2]. One particularly promising strategy involves detecting single electron ionization from DM-electron scattering [3-27].

Since momentum transfer from the DM to a target particle T is most efficient when m DM &gt; m T , a bound atomic electron can capture an order-1 fraction of the DM kinetic energy for m DM &gt; m e and be ionized. Similar reasoning would suggest that DM lighter than an atomic nucleus cannot efficiently transfer momentum to the nucleus, which is why experiments searching for nuclear recoils are typically insensitive to m DM &lt; m p . However, in a bound atomic system, momentum transfer to the entire atom will be redistributed among electrons and the nucleus through the electronic binding energy. This is known as the Migdal effect [28-34] and can also result in a final state with an ionized electron and a recoiling atom. Until now, the Migdal effect has solely been used to set limits on nucleon coupling for low mass weakly interacting massive particle models [35-38].

The main result of this paper is the following: in models where sub-GeV DM couples comparably to electrons and protons, the ratio of the differential ionization rate

∗ dbaxter@kicp.uchicago.edu

† yfkahn@illinois.edu

‡ krnjaicg@fnal.gov

dR M /d q due to the Migdal effect (which we will refer to as 'Migdal scattering' for brevity) to the corresponding direct electron scattering rate dR /d e q satisfies

$$\frac { d R _ { M } / d \mathbf q } { d R _ { e } / d \mathbf q } > Z ^ { 2 } \left ( \frac { m _ { e } } { m _ { N } } \right ) ^ { 2 } ( \mathbf q r _ { a } ) ^ { 2 }, \text{ \quad \ } ( 1 )$$

where m N is the mass of the target nucleus, Z its atomic number, q the 3-momentum transferred from the DM to the atom, and r a is an effective atomic radius which we will define more precisely in Sec. III B. 1 As we will show, the rate computation for both electron scattering and Migdal scattering involves identical atomic ingredients because the initial hard scatter factorizes from the dynamics of ionization in bound atoms. However, the scattering probability in the latter case is enhanced by Z 2 due to coherent scattering off the nucleus, and simultaneously suppressed by the small electron mass compared to the heavy nucleus, though this suppression is mitigated somewhat when the momentum transfer to the atom is large.

Due to these competing effects, the Migdal scattering rate in heavy atoms such as Xe is generically dominated by the largest kinematically-permitted momentum transfers, typically hundreds of keV, which are small on nuclear scales but large on electron scales; by contrast, direct DM-electron scattering is dominated by the smallest momentum transfers. Thus, the direct DM-electron scattering rate generically dominates over Migdal scattering for light-mediator exchange, which favors small momentum transfers, whereas Migdal scattering can dominate for heavier mediators and heavier DM which imparts

1 If DM couples equally to protons and neutrons, Z should be replaced by the mass number A .

larger momentum transfers to the target system. When Eq. (1) is integrated over the momentum transfer q , the total rate R M will then dominate over R e for sufficiently heavy DM.

This paper is organized as follows. In Sec. II we define our reference model of DM coupling to both electrons and protons. In Sec. III we develop in parallel the formalisms for electron and Migdal scattering, illustrating their similarities and differences. In Sec. IV we discuss the conversion from electron recoil spectra to observed ionization spectra in Xe and highlight the importance of the electron binding energies and wavefunctions in obtaining accurate exclusion curves. We conclude in Sec. V with a comparison of Migdal and electron exclusion curves for XENON10 and XENON100 data [39, 40], as well as the recent XENON1T results on light DM scattering [41]. We emphasize throughout that considerable theoretical uncertainty exists as to the correct choice of wavefunctions to use in computing limits on both Migdal and electron scattering. Consequently, our limits presented here should be considered provisional pending a dedicated analysis of relativistic and electron correlation effects in heavy atomic systems.

## II. REFERENCE MODEL

Our benchmark model consists of a DM candidate χ , which scatters off both electrons and protons through the exchange of a massive dark photon A ′ [42, 43]. The A ′ kinetically mixes with the visible photon, and after rotating away the kinetic mixing term glyph[epsilon1] 2 F µν F ′ µν , the Lagrangian for this scenario contains

$$\mathcal { L } = - \frac { 1 } { 4 } F _ { \mu \nu } ^ { \prime } F ^ { \prime \, \mu \nu } + \frac { m _ { A ^ { \prime } } ^ { 2 } } { 2 } A _ { \mu } ^ { \prime } A ^ { \prime \, \mu } + A _ { \mu } ^ { \prime } \left ( g _ { D } J _ { D } ^ { \mu } + \epsilon e J _ { \text{EM} } ^ { \mu } \right ), ( 2 )$$

where glyph[epsilon1] is the kinetic mixing parameter, m A ′ is the A ′ mass, J µ EM = ∑ f Q fγ f f ¯ µ is the electromagnetic current of all Standard Model fermions f with charges Q f , g D ≡ √ 4 πα D is the dark photon gauge coupling, and J µ D = ¯ χγ µ χ or iχ ∂ ∗ µ χ is the DM current for a Dirac fermion or complex scalar DM candidate, respectively. The fiducial non-relativistic cross section for χ scattering off a free-particle target T is defined at a reference 3-momentum transfer q 0 as

$$\overline { \sigma } _ { T } = \frac { 1 6 \pi \epsilon ^ { 2 } \alpha \alpha _ { D } \mu _ { \chi T } ^ { 2 } } { ( m _ { A ^ { \prime } } ^ { 2 } + | { \mathbf q } _ { 0 } | ^ { 2 } ) ^ { 2 } }, \quad \quad ( 3 ) \quad \ \ ( n \epsilon \\ b u )$$

where µ χT is the χ T -reduced mass; by coincidence this same parametric expression hold for both complex scalar and Dirac fermion DM candidates.

This popular scenario features comparable mediator couplings to electrons and protons, so it serves as a good benchmark for comparing DM-induced ionization from direct electron scattering and Migdal scattering. Indeed, both processes will always be present, so for the remainder of this paper, we will only consider this model. However, comparable quark and lepton couplings are by no means unique to dark photons. This feature applies to most anomaly-free U (1) extensions to the SM (e.g. gauged B -L ) whose gauge bosons couple to DM [44]; it is also generic for (pseudo)scalar-mediated DM scattering to feature comparable electron and proton couplings [45, 46].

We note for completeness that, in the general case with arbitrary mediator couplings to different particle species, the relationship between Migdal and electron scattering in Eq. (1) becomes

$$\underset { \mathbb { Z } \text{ in} } { \overset { \text{or} \text{-} } { \text{or} \text{-} } } \quad \frac { d R _ { M } / d \mathbf q } { d R _ { e } / d \mathbf q } \gtrsim \left ( \frac { Z g _ { p } + ( A - Z ) g _ { n } } { g _ { e } } \right ) ^ { 2 } \left ( \frac { m _ { e } } { m _ { N } } \right ) ^ { 2 } ( \mathbf q r _ { a } ) ^ { 2 }, \, ( 4 )$$

where g e and g p,n are the mediator's couplings to electrons, protons, and neutrons respectively. Thus, unless g p,n glyph[lessmuch] g e , there is a kinematic regime for which Migdal scattering can predominate over electron scattering.

## III. COMPARISON OF ELECTRON SCATTERING AND MIGDAL SCATTERING

In this section, we carefully define the kinematics and dynamics relevant for sub-GeV DM interacting with atoms through electron and Migdal scattering. We will work in the framework of the dark photon model described above in Sec. II, where the dark photon mediates DM-SM interactions and couples equally to electrons and protons with strength glyph[epsilon1] | e | , but our results are applicable to any model where the momentum dependence of DMelectron and DM-proton scattering (that is, the form factor F DM ( q ) defined below Eq. (16)) is identical.

## A. Kinematics

In both Migdal and electron scattering, the incoming and outgoing states are the same: a DM particle plus a bound atom and a DM particle plus an ionized atom plus an unbound electron, respectively. The incoming DM is assumed to be a plane wave, which is both an energy eigenstate and a momentum eigenstate. The incoming atom (at rest in the lab frame) is an energy eigenstate, and is also a momentum eigenstate for the total momentum of the atom p A = p N + ∑ Z i =1 p i , where the sum runs over the Z electrons in the electron cloud of the (neutral) atom. The outgoing DM is also a plane wave, but the outgoing atom can either be treated as an ionized atom with a separate ionized electron, or an atom in an excited state where the ionized electron belongs to the continuum spectrum of the atomic Hamiltonian. Following Ref. [32], we will take the second perspective where energy-momentum conservation is more transparent, in which case the entire atom recoils with velocity v A and has momentum p A = m A v A , where m A = m N + Zm e is the nominal mass of the atom neglecting binding energy. In all cases we will consider, it is appropriate to approximate m A by m N since the nucleus is so much heavier

than the electron cloud. The energetics of the ionized electron are accounted for by treating it as an excited state of the electron cloud.

To summarize, when treating the atom as a composite system of electrons and nucleus with a spectrum of internal energy levels, both energy and momentum are conserved in DM-atom interactions. For DM with mass m χ , incoming velocity v , and outgoing momentum p ′ χ , momentum conservation requires

$$\mathbf q \equiv m _ { \chi } \mathbf v - \mathbf p _ { \chi } ^ { \prime } = m _ { N } \mathbf v _ { A } \quad \quad ( 5 ) \quad \ \text{is a}$$

and energy conservation requires

$$\Delta E _ { e } = \frac { 1 } { 2 } m _ { \chi } v ^ { 2 } - \frac { | m _ { \chi } \mathbf v - \mathbf q | ^ { 2 } } { 2 m _ { \chi } } - \frac { \mathbf q ^ { 2 } } { 2 m _ { N } } = \mathbf q \cdot \mathbf v - \frac { \mathbf q ^ { 2 } } { 2 \mu _ { \chi N } } ( 6 ) \quad \stackrel { \mathbf \iota \infty } { \text{Hamilt} }$$

where µ χN = m m / m χ N ( χ + m N ) is the DM-nucleus reduced mass and ∆ E e ≡ E e,f -E e,i is the energy transferred to the scattered electron.

We emphasize that these kinematics are identical for electron scattering and Migdal scattering, provided the ionized electron is treated as a scattering state of the electron cloud Hamiltonian. In thinking of the nucleus and electrons as a single many-particle system in this formalism, it helps to regard q as simply the momentum transferred from the DM, rather than as a momentum transferred to any particular constituent of the target system. However, since the nucleus makes up the vast majority of the mass of the atom, one may think of q as the nuclear recoil momentum, as shown in Eq. (5).

## B. Dynamics

While the kinematics of Migdal and electron scattering are identical, their dynamics differ in a crucial way depending on whether DM interacts directly with electrons or nuclei. In the language of nonrelativistic quantum mechanics, the perturbing Hamiltonian for the DM-atom interaction in the case of electron scattering is

$$H _ { \text{int}, e } = - \int \frac { d ^ { 3 } { \mathbf q } } { ( 2 \pi ) ^ { 3 } } e ^ { i { \mathbf q } \cdot ( - { \mathbf x } _ { \chi } + \sum _ { i = 1 } ^ { z } { \mathbf x } _ { i } ) } \frac { \mathcal { M } _ { e \chi } ( q ) } { 4 m _ { \chi } m _ { e } }, \ \ ( 7 )$$

where M eχ ( q ) is the Lorentz-invariant matrix element for DM scattering off a free electron through 4-momentum transfer q ≈ (0 , q ), and x χ and x i are the position operators for the DM and electrons, respectively. Because the DM interacts directly with electrons, we can ignore the nuclear part of the atomic Hamiltonian, and the rate will be proportional to [3]

$$R _ { e } \in \left | \langle \Psi _ { f } | H _ { \text{int} } | \Psi _ { i } \rangle \right | ^ { 2 } \sim | \langle \psi _ { f } | e ^ { i \mathbf q \cdot \mathbf x } | \psi _ { i } \rangle | ^ { 2 }, \quad ( 8 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where we have made the approximation that the initialand final-state electron cloud wavefunctions (Ψ i and Ψ , f respectively) factorize such that only a single electron (with position operator x ) participates in a transition between the single-electron states ψ i and ψ f . The scale of atomic wavefunctions is parametrically the size of the atom: to make this precise, we define an effective atomic radius r a as r a ≡ 1 / | q | max , where | q | max is the momentum transfer at which the matrix element in Eq. (8) is maximized, for a given choice of initial and final states. Thus we may expect the electronic matrix element in Eq. (8) to be as large as O (1) when | q | r a = O (1). Typical r a are on the scale of the Bohr radius a 0 , which is ∼ 1 / keV in natural units, and hence q · x glyph[greaterorsimilar] O (1) where the atomic wavefunctions have support. Note that for DM heavier than 1 MeV, a momentum transfer of 1 /r a is always kinematically allowed, since the DM carries momentum of at least 1 keV.

For Migdal scattering, where the fundamental DMatom interaction is with the nucleus, the interaction Hamiltonian is

$$H _ { \text{int}, N } = - \int \frac { d ^ { 3 } { \mathbf q } } { ( 2 \pi ) ^ { 3 } } e ^ { i { \mathbf q } \cdot ( { \mathbf x } _ { N } - { \mathbf x } _ { \mathbf x } ) } \frac { { \mathcal { M } } _ { N _ { \mathbf x } } ( q ) } { 4 m _ { \mathbf x } m _ { N } }, \quad ( 9 )$$

where x N is the position operator for the nucleus. H int ,N does not involve the electron position operators x i , so by itself, it cannot induce electronic transitions. However, the light-crossing time of the electron cloud is ∼ nm /c ∼ 5 keV -1 , so the timescale of momentum transfer to the nucleus is 'fast' as long as the mediator mass m A ′ satisfies m A ′ glyph[greatermuch] keV. 2 In this regime, the entire atom suddenly acquires velocity v A but leaves behind its stationary electrostatic potential; the electrons of the moving atom are no longer in energy eigenstates of the old electron cloud Hamiltonian. As a result, electronic transitions can arise, but not through a matrix element with the perturbing Hamiltonian H int ,N .

Rather, following Ref. [32], we construct approximate energy eigenstates of the moving atom by applying a Galilean transformation with velocity parameter v A , which results in a final-state atomic wavefunction containing a phase exp( i ∑ Z i =1 q e · x i ) multiplying the full wavefunction of the atom at rest, where q e ≡ m e v A . 3 Consequently, as shown in Ref. [32], the matrix element of the nuclear wavefunction with H int ,N results in a factor of M eN ( q ) times the overlap of the electronic wavefunctions, and thus

$$R _ { M } \in | \langle \Psi _ { \nu _ { A } } | \Psi _ { i } \rangle | ^ { 2 } \sim | \langle \psi _ { f } | e ^ { i { \mathbf q } _ { e } \cdot \mathbf x } | \psi _ { i } \rangle | ^ { 2 }, \quad ( 1 0 )$$

2 By construction, this always holds for contact interactions, F DM = 1, but it is not clear to us if the formalism developed in Ref. [32] remains valid for ultralight mediators with F DM ∝ q -2 , where for sufficiently small momentum transfers, the timescale for momentum transfer can be slow enough that the atomic state changes adiabatically.

3 It is important to note that q e is not to be interpreted as the momentum of the outgoing electron, but rather the effective momentum which appears in the matrix element due to the Galilean transformation to account for the atomic recoil velocity v A . The advantage of defining this quantity is to explicitly enable the dipole approximation of the matrix element for the masses considered here. Care should be taken not to place too much physical interpretation in this quantity, and so we translate our final result in Eq. (12) below back in terms of the physical momentum of the atomic system q .

where Ψ v A is the Galilean transformation of the initial state Ψ i of the electron cloud, with velocity parameter v A . We have made the same approximation as in Eq. (8) that only a single electronic transition contributes; note that q e instead of q appears in the exponent.

By momentum conservation from Eq. (5),

$$\mathbf q _ { e } = \frac { m _ { e } } { m _ { N } } \mathbf q. \quad \quad \quad ( 1 1 ) \quad \ \ s y \quad \\ \ \ \ \ \ \ \ k ^ { \prime } _ { \dots }$$

Unlike the case for electron scattering, where q · x glyph[greaterorsimilar] 1 for all DM masses greater than an MeV, q e · x glyph[lessmuch] 1 for all sub-GeV DM because m /m e N glyph[lessmuch] 1. Indeed, r a ∼ 4 a 0 for xenon, and q e &lt; 1 /r a as long as m χ &lt; 100 GeV. As q e → 0, the matrix element in Eq. (10) must vanish because ψ f and ψ i are energy eigenstates of the same Hamiltonian with different energy eigenvalues, by assumption. Hence the leading order term in the Taylor expansion of the exponential is linear in q e , and R M scales as q 2 e = q 2 ( m /m e N ) 2 . There are also additional selection rules now that the matrix element has a dipole form, 〈 ψ f | x | ψ i 〉 , but in general, for a given q and a choice of initial- and final-state wavefunctions, the ratio of Migdal and electron scattering spectra for each i → f transition scales as

$$\frac { d R _ { M } / d \mathbf q } { d R _ { e } / d \mathbf q } > \mathbf q ^ { 2 } r _ { a } ^ { 2 } \left ( \frac { m _ { e } } { m _ { N } } \right ) ^ { 2 } \frac { | \mathcal { M } _ { \chi N } ( q ) / m _ { N } | ^ { 2 } } { | \mathcal { M } _ { \chi e } ( q ) / m _ { e } | ^ { 2 } }, \ \ ( 1 2 )$$

where the appearance of r a arises from the expectation value of x , which is parametrically of order r a . We have written the above relation as an inequality because, for sufficiently large momentum transfers ( | q | glyph[greatermuch] keV), the exponential in Eq. (8) will oscillate rapidly and R e will become suppressed, thereby enhancing the Migdal rate relative to the electron scattering rate. On the other hand, the Migdal spectrum may always be approximated by a dipole matrix element which scales with q 2 r 2 a , since q e &lt; /r 1 a for all targets relevant for sub-GeV DM.

## C. Spectra and rates

To compute the ionization rate for both processes, we must integrate over the momentum transfer q and the DM velocity v , and sum over the final electronic states ψ f , weighted by a delta function enforcing energy conservation (3-momentum conservation has already been enforced in the definitions of q and q e above). 4 We perform the integral over v by approximating the DM velocity distribution as spherically symmetric, f ( v ) = f ( v ), such that the total rate between initial state i and final state f is

$$R _ { i \rightarrow f } = \frac { \rho _ { \chi } } { m _ { \chi } } \int d ^ { 3 } v \, f ( v ) \, \sigma v _ { i \rightarrow f }, \quad \ \ ( 1 3 )$$

4 Note that integrating over q is equivalent to integrating over the nuclear recoil energy E R ≈ q 2 / (2 m N ), since in this paper we are concerned only with the electronic energy spectrum.

where ρ χ is the local DM density. For the sum over final states, we choose the normalization [3, 16]

$$\sum _ { f } = \frac { 1 } { 2 } \sum _ { l ^ { \prime } m ^ { \prime } } \int \frac { k ^ { \prime 3 } d \ln E _ { e } } { ( 2 \pi ) ^ { 3 } },$$

which is appropriate for scattering states in a sphericallysymmetric potential which have asymptotic momentum k ′ = √ 2 m E e e and angular momentum quantum numbers l ′ and m ′ . Here, E e is the recoil energy of the ionized electron asymptotically far away from the ionized atom; from now on our final state f will always be a scattered electron of energy E e , and the initial state i will be a bound state of (negative) energy E nl indexed by principal quantum number n and angular momentum quantum number l , appropriate for a spherically-symmetric atom ignoring spin-orbit coupling and relativistic effects. The only difference between Migdal and electron scattering in the above procedure is the expression for σv i → f .

To perform the integral over q and compute σv i → f we must specify the free-particle matrix elements. In the dark photon model, we can define a spin-averaged fiducial cross section for DM χ scattering off an isolated target T with charge | e | as in Eq. (3). For T = p, e , these fiducial cross sections satisfy

$$\frac { \overline { \sigma } _ { e } } { \mu _ { \chi e } ^ { 2 } } = \frac { \overline { \sigma } _ { p } } { \mu _ { \chi p } ^ { 2 } },$$

so σ e and σ p are proportional by a factor which only depends on the DM mass χ . To emphasize the point that σ e and σ p are related in this model, we shall refer to σ e as simply σ .

The fiducial cross section defined in Eq. (3) is related to the free-particle scattering matrix element as

$$| \mathcal { M } ( q ) | ^ { 2 } = \frac { 1 6 \pi m _ { \chi } ^ { 2 } m _ { T } ^ { 2 } \overline { \sigma } _ { T } } { \mu _ { \chi T } ^ { 2 } } | F _ { \text{DM} } ( q ) | ^ { 2 } \, \quad ( 1 6 )$$

where we have assumed that the appropriate electron and DM spins have been summed and/or averaged. Here, F DM ( q ) is the DM form factor which parametrizes all momentum dependence in the free-particle matrix element: if m A ′ glyph[lessmuch] m v χ , F DM ( q ) ∝ 1 /q 2 , while if m A ′ glyph[greatermuch] m v χ , F DM ( q ) = 1. Note that in the dark photon model with equal proton and electron couplings,

$$\frac { | \mathcal { M } _ { \chi N } ( q ) | ^ { 2 } } { m _ { N } ^ { 2 } } = Z ^ { 2 } | F _ { N } ( q ) | ^ { 2 } \frac { | \mathcal { M } _ { \chi e } ( q ) | ^ { 2 } } { m _ { e } ^ { 2 } } \quad ( 1 7 )$$

in Eq. (12), where F N is the form factor of the nucleus; this relation between the matrix elements gives Eq. (1).

Putting all the pieces together, the electron recoil spectrum per unit detector mass for both electron and Migdal scattering is

$$\frac { d R _ { e, M } } { d \ln E _ { e } } = N _ { T } \frac { \rho _ { \chi } } { m _ { \chi } } \frac { \overline { \sigma } } { 8 \mu _ { \chi e } ^ { 2 } } I _ { e, M } ( E _ { e } ), \quad \ \ ( 1 8 )$$

where N T is the number of atomic targets and

$$\L \, \text{we are} \quad I _ { e, M } ( E _ { e } ) = \, \int d | \mathbf q | | \mathbf q | | F _ { \text{DM} } ( q ) | ^ { 2 } \eta ( v _ { \min } ) | f _ { e, M } ( E _ { e }, \mathbf q ) | ^ { 2 }. \, ( 1 9 )$$

Here, we have solved the delta function for energy conservation, δ E ( e -E nl + q 2 2 µ χN -q · v ), to perform the integral over the DM velocity distribution, resulting in a factor of η v ( min ) ≡ 〈 v -1 θ v ( -v min ) 〉 , the mean inverse DM speed in the lab frame, as a function of

$$v _ { \min } = \frac { \Delta E _ { e } } { | \mathbf q | } + \frac { | \mathbf q | } { 2 \mu _ { \chi N } } = \frac { | E _ { n l } | + E _ { e } } { | \mathbf q | } + \frac { | \mathbf q | } { 2 \mu _ { \chi N } }, \ ( 2 0 ) \quad \text{$\text{$\text{$kiner}$} } \\ \text{and $\epsilon$} \\ \text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text{$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$\text$} }$$

which is the minimum DM velocity required to ionize the target electron through a momentum transfer | q | . The lab frame velocity distribution is cut off at v E + v esc . , where v E ∼ 240 km s is the average speed of the Earth / relative to the DM halo, and v esc . = 544 km s is the / galactic escape velocity (these parameters are chosen to facilitate comparisons with Ref. [16]).

The differences between the Migdal and electron scattering processes are entirely contained in the ionization form factors | f e,M ( E , e q ) | 2 , which are independent of all DM properties and depend only on the electronic and nuclear structure of the target. We will discuss in some detail in Sec. IV B the issues with accurately computing the atomic wavefunctions required for these ionization form factors. For electron scattering,

$$| f _ { e } ( E _ { e }, { \mathbf q } ) | ^ { 2 } = \frac { k ^ { \prime 3 } } { 4 \pi ^ { 3 } } \times 2 \sum _ { n, l, l ^ { \prime }, m ^ { \prime } } | \langle \psi _ { E _ { e } } ^ { f } | e ^ { i { \mathbf q } \cdot { \mathbf x } } | \psi _ { E _ { n l } } ^ { i } \rangle | ^ { 2 } ( 2 1 ) \quad \stackrel { E _ { e } } { | f _ { e, M } } \\ \text{plot } \varsigma \\ \text{bv } n ( \cdot }$$

where ψ i E nl is a bound orbital of energy E nl with unit norm (the factor of 2 accounts for the approximate spin degeneracy of occupied states) and ψ f E e is an unbound electronic state of energy E e (indexed by the continuously-valued energy E e and the angular momentum quantum numbers l ′ and m ′ ), normalized to

$$\langle \psi _ { E _ { 1 }, l, m } | \psi _ { E _ { 2 }, l ^ { \prime }, m ^ { \prime } } \rangle = \frac { ( 2 \pi ) ^ { 3 } } { k _ { 1 } ^ { 2 } } \delta ( k _ { 1 } - k _ { 2 } ) \delta _ { l l ^ { \prime } } \delta _ { m m ^ { \prime } } \, ( 2 2 ) \quad \text{A co} \\ \text{$\text{$\text{$atom}$} }$$

where k 1 2 , = √ 2 m E e 1 2 , . For Migdal scattering, the analogous ionization form factor is

$$| f _ { M } ( E _ { e }, \mathbf q ) | ^ { 2 } & = \frac { k ^ { \prime 3 } } { 4 \pi ^ { 3 } } Z ^ { 2 } | F _ { N } ( q ) | ^ { 2 } & \begin{array} { c } \begin{array} { c } \begin{array} { c } \begin{array} { c } \end{array} \\ \end{array} \\ \end{array} \\ & \times 2 \sum _ { n, l, l ^ { \prime }, m ^ { \prime } } | \langle \psi _ { E _ { e } } ^ { f } | e ^ { i \mathbf q _ { e } \cdot \mathbf x } | \psi _ { E _ { n l } } ^ { i } \end{array} \rangle | ^ { 2 }. & \begin{array} { c } \begin{array} { c } \begin{array} { c } \end{array} \\ \end{array} \\ \end{array} \begin{array} { c } \begin{array} { c } \begin{array} { c } \begin{array} { c } \begin{array} { c } \end{array} \\ \end{array} \\ \end{array}$$

The differences with respect to electron scattering are the appearance of Z 2 from coherent scattering off all the protons in the nucleus, a nuclear form factor F N parametrizing loss of coherence at large momentum transfers (which is largely irrelevant for the sub-MeV momentum transfers typical of sub-GeV DM), and the appearance of q e instead of q in the matrix element between initial and final states.

The key quantity controlling the relative size of Migdal and electron scattering rates is q . From Eq. (20), the smallest allowed | q | is

$$| \mathbf q | _ { \min } = \frac { E _ { b } } { v _ { \max } } \, \quad \ \ ( 2 4 ) \quad \ \ \.$$

where E b is the first ionization energy (positive by convention) of the atom in question, and v max is the largest possible DM speed, which is the Galactic escape velocity in the lab frame. Note that | q | min is independent of the DM mass: for xenon with E b ∼ 12 eV, and v max ∼ 770 km s, / | q | min ∼ 5 keV. Thus for all kinematically-allowed momentum transfers, | q | r a &gt; 1, and electron scattering is dominated by the smallest possible q before f e ( E , e q ) is suppressed by the quicklyoscillating exponential in the matrix element. On the other hand, the largest allowed | q | is

$$| { \mathbf q } | _ { \max } = 2 \mu _ { \chi N } v _ { \max } \sim 5 \ k e V \left ( \frac { m _ { \chi } } { \text{MeV} } \right ), \quad ( 2 5 ) \\ \dots \quad \dots \quad \dots \quad \dots \quad \dots \quad \dots$$

which grows with DM mass and can be as large as hundreds of keV for m χ = O (100 MeV). For these momentum transfers, | f M ( E , e q ) | 2 still does not feel any suppression from the nuclear form factor F N , which is still ∼ 1 for | q | glyph[lessorsimilar] MeV, and likewise is still in the regime of small q e and so grows with q 2 . Thus the Migdal ionization form factor is largest when q is the largest, and Migdal scattering is dominated by the largest kinematically-allowed momentum transfers. 5

Weillustrate this behavior in Fig. 1 for m χ = 300 MeV, E e = 10 eV, and F DM = 1. The left plot shows | f e,M ( E , e q ) | 2 evaluated at E e = 10 eV, and the right plot shows the integrand of Eq. (19) which is weighted by η v ( min ). The q values plotted span the kinematicallyallowed range between | q | min and | q | max , as can be seen from the right plot where the velocity distribution cuts off the integrand at small and large q . To compute f e,M for both electron and Migdal scattering from the same set of wavefunctions, the orthogonality of ψ f and ψ i is crucial, and to ensure this, both bound and free wavefunctions must be constructed from the same atomic Hamiltonian. A complete treatment would require a full numerical solution to the many-body Schr¨ odinger equation for the atom in question, but here we capture the essential features by using hydrogenic wavefunctions for the 5p shell of Xe and matching the effective nuclear charge to the binding energy of the 5p state, with scattering states constructed from the same hydrogenic potential. This unphysical choice of wavefunctions is for illustrative purposes only; the wavefunctions used in the remainder of this paper are discussed in detail in Sec. IV B. We note that for a DM form factor proportional to 1 /q 2 , as would be the case for an ultralight dark photon mediator, the spectrum integrand is weighted by 1 /q 4 ∼ 1 / q 4 which heavily suppresses the Migdal spectrum compared to the electron spectrum for all DM masses. 6 For these form

5 In principle, there should be interference between electron and Migdal scattering, which have identical final states, but the distinct kinematics of these two processes should minimize these effects.

6 As noted in Footnote 2, long-range interactions may result in adiabatic rather than sudden changes in atomic states during Migdal scattering, which would further suppress electronic transitions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-791f9c87-24b8-4519-887b-6ecaebfe96db.jpg)

FIG. 1. Qualitative comparison of electron (green) and Migdal (purple) ionization form factors and spectrum integrands for m χ = 300 MeV, E e = 10 eV, and F DM = 1. The form factors from Eqs. (21) and (23) are plotted in the left panel, and the integrand from Eq. (19) (which is weighted by the inverse mean DM speed) is plotted in the right panel. The electron scattering rate is dominated by small q while the Migdal scattering rate is dominated by large q . These plots are computed using initial- and final-state electronic wavefunctions constructed from a hydrogenic potential for Xe, but the qualitative features are independent of the choice of wavefunctions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-44328742-035e-49ff-8963-d4a43b661ed0.jpg)

factors, electron scattering always dominates over Migdal scattering by several orders of magnitude, and as such, for the remainder of this paper we will focus on the case F DM = 1.

light produced from energy deposited in liquid xenon.

We can confirm the relative strength of Migdal and electron scattering by considering the full electron recoil spectrum, as shown in Fig. 2. Here, to facilitate comparison with the literature, the electron spectra are calculated using the wavefunctions and binding energies of Ref. [16], and the Migdal spectra are calculated using f M as tabulated in Ref. [32]. Despite some differences in these wavefunctions and binding energies (which we discuss further in Sec. IV B), the intuition developed above holds well: for sufficiently large DM masses and equal couplings to protons and electrons, the Migdal spectrum dominates over the electron spectrum for all electron recoil energies.

## IV. NUMERICAL MODELING AND SYSTEMATIC UNCERTAINTIES

To apply the formalism of the previous section to experimental data, we must choose a model for generating ionization spectra from recoil spectra, as well as a set of atomic wavefunctions. We choose the previously published low-threshold analysis [16] of the XENON10 [47, 48] and XENON100 [49] detectors, as well as the newly-released S2-only analysis [41] from the XENON1T detector [50]. These data are chosen for containing a relatively large exposure (for this mass range) at extremely low thresholds. All of these experiments use xenon time projection chambers to measure charge and

An interaction in the xenon will create some number of xenon ions N i and initially excited xenon atoms N glyph[star] . These atoms will form dimer states in the xenon which will release energy in the form of charge and UV scintillation photons. Scintillation photons produced at this step are immediately detected in what is referred to as the S1 signal. Before they can recombine, emitted electrons are drifted in a ∼ 300 V/cm electric field to a liquidgas interface where a stronger ( ∼ 10 kV/cm) extraction field is used to accelerate the electrons into the gas, producing a second (amplified) burst of light, referred to as the S2 signal. It has been well measured that (relatively speaking) interactions with xenon nuclei preferentially deposit energy via S1, whereas interactions with electrons in the detector preferentially deposit energy via S2 [51-54]. In the case of sub-GeV DM interacting with a xenon atom through either electron or Migdal scattering, the momentum transfer to the recoiling atom is sufficiently small that the S1 signal is expected to be effectively zero. Thus, we consider an S2-only analysis using a 1(4) electron threshold for XENON10(100) [16] and a 5 electron threshold for XENON1T [41].

## A. Ionization model and quantization

To compare with data, we must quantize the calculated recoil spectra in terms of the number of electrons extracted. For this, we adopt the ionization model from Refs. [3, 16] to determine the number of electrons produced ( n e ) from an initial energy transfer ∆ E e which

FIG. 2. Comparison of electron (dashed) and Migdal (solid) spectra for reference values σ = 5 × 10 -39 cm 2 and m χ = 100 MeV (blue) and 300 MeV (red). For both spectra, we show the inclusive rates summed over all E nl → E e transitions in xenon, where contributions from nl = 5p, 5s, and 4d dominate. Migdal spectra are computed using the wavefunctions and binding energies from Ref. [32], while electron spectra are computed using wavefunctions and binding energies from Ref. [16]. The differences between these choices are irrelevant for our qualitative argument here and are discussed further in Sec. IV B.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-fb118ff6-180d-4199-9f39-725d6dc241f3.jpg)

ionizes an electron from a specified electron shell with binding energy E nl to the continuum with energy E e . We begin by considering the ejected electron, which has a probability f R to recombine (avoiding detection). According to the Thomas-Imel recombination model, f R is determined to be very small at low energies [51, 55] in good agreement with measurement [56] and is assumed to be zero for this analysis. We can thus write the probability of observing an initially produced electron as

$$f _ { 0 } = \frac { 1 - f _ { R } } { 1 + ( N _ { * } / N _ { i } ) } \approx 0. 8 3, \quad \ \ ( 2 6 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

where the ratio of initially excited atoms to initially ionized atoms satisfies N /N glyph[star] i ≈ 0 2 at high energies [57, 58]. .

At high energies, the average energy W required to produce one charge quantum in xenon is measured to be W = 13 7 eV [51]. . To convert E e and E nl into an expected quantized signal, we consider n t trials of a binomial process with probability of success, f 0 , which satisfy

$$n _ { t } = \text{floor} \left ( \frac { E _ { e } } { W } \right ) + \text{floor} \left ( \frac { | E _ { n l } | - E _ { b } } { W } \right ), \quad ( 2 7 ) \quad \text{spin}.$$

where | E nl | -E b is the available de-excitation energy. Thus, we can write

$$n _ { e } = ( 1 - f _ { R } ) + \binom { n _ { t } } { f _ { 0 } }, \quad \ \ ( 2 8 ) \quad \ \ _ { b } ^ { 7 }$$

TABLE I. Comparison of the ionization energy E b and electron binding energies | E nl | (eV) of the 5p, 5s, and 4d electron shells of xenon from calculations using the formalisms of Refs. [16, 32]. The measured values [62, 63] tend to fall somewhere in the middle.

| Method                |   E b |   | E 5p | |   | E 5s | |   | E 4d | |
|-----------------------|-------|------------|------------|------------|
| Ibe et al. [32]       |   9.8 |        9.8 |       21   |       61   |
| Essig et al. [16, 61] |  12.4 |       12.4 |       25.7 |       75.6 |
| Measured [62, 63]     |  12.1 |       12.8 |       23.3 |       68.5 |

where we assume that the number n e of quanta produced is equivalent to the number of electrons extracted (observed), as the extraction efficiency should be ∼ 100% [59, 60]. Examples of quantized spectra for the F DM = 1 case and different DM masses are shown in Fig. 3. This simple ionization model is sufficient for comparing electron and Migdal scattering here, but a more robust model would be needed to correctly interpret a signal through either channel.

Upper limits on the number of events for each value of n e in XENON10(100) have been determined in Ref. [16] to be r 1 &lt; 15 18, . r 2 &lt; 3 37, . r 3 &lt; 0 95, . and r 4 &lt; 0 35(0 17) counts kg . . -1 day -1 . Upper limits from XENON1T [41] are not given directly. Instead, an upper bound of 22.5 events is reported in the range 165-275 photoelectrons in 15 ton-days of exposure 7 . We take the measured ratio of photoelectrons detected per n e (single electron gain) to be ∼ 33 [41] and conservatively assume that all events in the reported range are for the lowest encompassed bin n e = 5, to obtain r 5 &lt; 0 0015 counts . kg -1 day -1 . These rates already account for analysis and detector efficiencies and thus can be directly compared against our quantized spectra.

## B. Electron Binding Energies and Wavefunctions

The dominant quantity controlling the sensitivity for small DM masses ( glyph[lessorsimilar] 100 MeV) is the outer-shell binding energy of xenon, or equivalently its first ionization energy. This quantity is well-measured experimentally with a value of 12.1 eV [62]. Following Ref. [32], we define the electron shell binding energies E nl as an average over angular momentum states κ ,

$$E _ { n l } = \frac { 1 } { 2 } \sum _ { \kappa } \delta _ { l, | \kappa + 1 / 2 | - 1 / 2 } E _ { n \kappa }, \text{ \quad \ \ } ( 2 9 )$$

but define the ionization energy E b as the minimum binding energy for all spin states in the atom. As a result, if spin-orbit coupling is ignored, E b = | E 5p | , but the measured values show a ∼ 5% difference between the two, as can be seen in Table I.

7 The average exposure for the range 165-198 photoelectrons has been determined from Figs. 1 and 4 of Ref. [41].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-b8301af3-9430-4cda-9f1d-e5d1ca35f596.jpg)

FIG. 3. Quantized rate spectra for DM-electron scattering (left) and DM-Migdal scattering (right) in xenon per number of electrons observed in the case of a heavy mediator ( F DM = 1) for various DM masses between 10 and 1000 MeV. We normalize these spectra to an exposure of 1000 kg-yr and a fiducial cross section of σ = 5 × 10 -39 cm . 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-d50d59ec-1d38-4295-9109-31ff61a88309.jpg)

FIG. 4. Effect of wavefunction and binding energy choices on 90% CL limits on σ for direct DM-electron scattering in XENON10 with n e = 1 (blue), 2 (green), and 3 (purple) and in XENON100 with n e = 4 (red). Solid curves are published limits from Ref. [16] and dashed curves use binding energies from Ref. [32] to construct hydrogenic final-state wavefunctions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-93902da2-429a-4ce9-b982-95a835ff037f.jpg)

in Ref. [16] takes the calculated binding energy from the Roothaan-Hartree-Fock atomic wavefunctions tabulated in Ref. [61] and constructs outgoing wavefunctions from a hydrogenic potential with a shell-dependent effective nuclear charge which reproduces the appropriate binding energy for each electron shell. As the binding energies from Ref. [61] are closer to the observed values, this procedure retains the physical binding energies at the cost of losing orthogonality between initial and final electronic states, as well as any electron correlation effects. This orthogonality is crucial in order to obtain the behavior of f M as a function of q , as discussed above, so it is not possible to compute Migdal scattering rates using these wavefunctions. However, as noted in Sec. III C, electron scattering is dominated by the region where q · x glyph[greaterorsimilar] 1, so the form factor never probes the region where the matrix element must vanish as q → 0; thus, the overall kinematic features of electron scattering are sufficiently captured by this formalism.

Ideally, one wishes to compute the ionization spectrum using atomic wavefunctions with energy eigenvalues matching the observed binding energies. The binding energies used in the two analyses in Refs. [16, 32] are compared in Table I. As can be seen from this table, obtaining accurate binding energies is not entirely trivial, as the FAC code [64] used in Ref. [32] gives an outer-shell binding energy of 9.8 eV, a significant difference of 20% from the observed value. On the other hand, the procedure used

While our interest in this paper is primarily a qualitative comparison of electron and Migdal scattering, we estimate one source of systematic error which can affect both processes by computing electron scattering rates using hydrogenic final-state wavefunctions (as in Ref. [16]), but constructed from the systematically-lower binding energies used in Ref. [32]. The main difference from the illustrative example computed in Sec. III C is that the initial-state wavefunctions are now taken from Ref. [61] instead of using a crude hydrogenic potential, in order to isolate the effects of binding energies and final-state wavefunctions. Fig. 4 compares the cross section limits obtained from these binding energies to the published electron scattering limits from Ref. [16]. The systematic error on the electron scattering case associated with this procedure is less than an order of magnitude over the

FIG. 5. Limits on σ (defined below Eq. (15); to compare with standard conventions, the reference cross section for electron scattering is σ e = σ , and the reference cross section for proton scattering is σ p = σ × µ 2 χp /µ 2 χe ) for heavy A ′ ( F DM = 1) mediated DM-electron (dashed purple) and DM-Migdal (dashed red) scattering are shown for n e = 4 for XENON100 (left) and for n e = 5 for XENON1T (right). In the mass and coupling range shown on the plot, the XENON10 limits are sub-dominant. For comparison, we show the published electron scattering limits [16, 41] computed with hydrogenic final-state wavefunctions and binding energies from Ref. [16] (solid purple); our electron scattering results use the smaller (unphysical) binding energies from [32] to facilitate a comparison with Migdal scattering using the same binding energies (see Sec. IV B). The thick blue curve is the complex scalar DM freeze-out target (particle-antiparticle symmetric DM population). Points along this curve account for the full DM abundance as long as m A ′ glyph[greatermuch] m χ ; near resonance at m A ′ ≈ 2 m χ this target moves down in the parameter space, but is otherwise robust [65, 66]. The thin blue curve is the looser asymmetric Dirac fermion DM target. Any points above this line can account for the full DM abundance, but with different particle-antiparticle asymmetries [2, 65]; points below this curve are excluded by Planck limits on cosmic microwave background energy injection from the annihilation of the symmetric component [67]. The dotted blue curve taken from Ref. [2] represents sensitivity targets for ELDER DM [68]; points above this curve correspond to SIMP DM models with the same A ′ mediator considered here [69]. Shaded regions represent an envelope of exclusions from beam dump searches (LSND [70], E137 [71, 72], and MiniBooNE [73, 74]), nuclear recoil direct detection limits from CRESST II [75], and the BaBar monophoton search for invisibly decaying dark photons [76-78].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-e6506669-bb18-41fd-974a-54eaa7e49123.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00012v3-figure-b5a92710-e14d-477a-b9c7-5ee0330e2156.jpg)

full DM mass range, but as expected, smaller binding energies lead to stronger cross section limits. The same procedure cannot be directly applied to Migdal scattering as the wavefunctions from Ref. [16] are not orthogonal, but the magnitude of the difference should be comparable (about a factor of 2 for masses above 100 MeV where Migdal scattering dominates).

We leave to future work a precise determination of experimental limits on Migdal scattering using more sophisticated quantum chemistry codes which correctly reproduce the observed binding energies with orthogonal wavefunctions. Indeed, recent progress for electron scattering has already been made by incorporating relativistic effects and electron-electron interactions into a manybody calculation, using the observed ionization energies from photoabsorption data as a figure of merit for the quality of the wavefunctions [79]; the result is that at large n e , the spectrum differs significantly from that obtained with hydrogenic final-state wavefunctions, potentially affecting the limits at large DM masses by an order of magnitude. In particular, the size of relativistic effects will grow as the atomic number of the atom increases, so this may be a significant source of systematic uncertainty for Migdal scattering in xenon.

## V. RESULTS AND CONCLUSION

In this paper we have placed sub-GeV DM detection via electron and Migdal scattering on equivalent theoretical footing. Intriguingly, we have found that if DM couples comparably to electrons and protons through a contact interaction ( F DM = 1), the Migdal rate can dominate for masses above ∼ 100 MeV. Thus, all existing limits for electron scattering in such models (such as dark photonmediated scenarios), including those from XENON10, XENON100, and XENON1T [16, 41], have omitted the dominant signal component at higher DM masses. In Fig. 5, we recalculate the full signal for DM-xenon scattering in XENON100 and XENON1T and extract upper bounds on σ which include both electron and Migdal scattering. It is clear that by exploiting the combination of electron and Migdal scattering, experiments with exposures and background rates comparable to XENON1T can start to probe the target parameter space for complex scalar DM freezing out through a heavy dark photon,

m A ′ glyph[greatermuch] m χ , but as we have emphasized, a definitive conclusion requires a more careful treatment of the atomic wavefunctions than has been used previously in the literature. In the event of a signal in the DM mass range where electron and Migdal rates are within a few orders of magnitude of each other, the unique spectral shapes can be used as a discriminant, though interference effects should be carefully considered.

Although our treatment here has focused on scattering from isolated atoms, we note that additional ionization from Migdal scattering should also contribute in semiconductor targets (mainly Si and Ge), for which low electronic band gaps represent the next frontier in electronionization direct detection. This additional signal channel can be probed by numerous future and ongoing experiments, including DAMIC at SNOLAB [80], SENSEI [26], SuperCDMS [24], and DAMIC-M [81]. However, a proper comparison of electron scattering and Migdal scattering in such materials is beyond the scope of the present work and deserves a dedicated study. At a minimum, the formalism for Migdal scattering must incorporate the nontrivial harmonic potential between neighboring ions, which may result in some portion of the DM energy loss appearing as phonons rather than electronic excitations.

Despite its robust theoretical underpinnings, Migdal scattering has not yet been experimentally observed. We emphasize that the same theoretical uncertainties which apply to Migdal scattering are present for the case of DMelectron scattering. We take advantage of this fact to show that the systematic uncertainty in the rate calculations due to different computations of the xenon binding energies can be estimated to be a factor of a few at high masses, where the Migdal scattering rate is dominant, but other sources of theoretical uncertainty due to relativistic effects and electron correlations may be equally

- [1] G. Bertone and D. Hooper, Rev. Mod. Phys. 90 , 045002 (2018), arXiv:1605.04909 [astro-ph.CO].
- [2] M. Battaglieri et al. , in U.S. Cosmic Visions: New Ideas in Dark Matter College Park, MD, USA, March 23-25, 2017 (2017) arXiv:1707.04591 [hep-ph].
- [3] R. Essig, J. Mardon, and T. Volansky, Phys. Rev. D85 , 076007 (2012), arXiv:1108.5383 [hep-ph].
- [4] R. Essig, A. Manalaysay, J. Mardon, P. Sorensen, and T. Volansky, Phys. Rev. Lett. 109 , 021301 (2012), arXiv:1206.2644 [astro-ph.CO].
- [5] P. W. Graham, D. E. Kaplan, S. Rajendran, and M. T. Walters, Phys. Dark Univ. 1 , 32 (2012), arXiv:1203.2531 [hep-ph].
- [6] H. An, M. Pospelov, J. Pradler, and A. Ritz, Phys. Lett. B747 , 331 (2015), arXiv:1412.8378 [hep-ph].
- [7] S. K. Lee, M. Lisanti, S. Mishra-Sharma, and B. R. Safdi, Phys. Rev. D92 , 083517 (2015), arXiv:1508.07361 [hep-ph].
- [8] R. Essig, M. Fernandez-Serra, J. Mardon, A. Soto, T. Volansky, and T.-T. Yu, JHEP 05 , 046 (2016),

important. The uncertainty due to the ionization model does not significantly affect the relative comparison of our calculated electron and Migdal scattering limits, but does explain the bulk of the difference between our limits and those in Ref. [38] when interpreted in the dark photon model.

This analysis further highlights the importance of developing low-energy calibration techniques. We have shown that Migdal and electron scattering processes probe the same atomic wavefunctions, but in different kinematic regimes. As noted in [79], atomic many-body effects are crucial for understanding DM-atom interactions. Calibrations of both ionization processes, Migdal scattering and electron scattering, would help to resolve the theoretical uncertainty in the wavefunctions, which we believe has been substantially underappreciated by the sub-GeV DM community.

Acknowledgments. We thank Rouven Essig, Noah Kurinsky, Masahiro Ibe, Josef Pradler, Lucas Wagner, and Tien-Tien Yu for helpful conversations. We additionally thank Evan Shockley for helpful clarification regarding the recent XENON1T results. Fermilab is operated by Fermi Research Alliance, LLC, under Contract No. DE-AC02-07CH11359 with the US Department of Energy. This work was supported in part by the Kavli Institute for Cosmological Physics at the University of Chicago through an endowment from the Kavli Foundation and its founder Fred Kavli. We thank the KICP pheno journal club for providing a stimulating environment for discussion where this work originated. This research was supported by the Munich Institute for Astroand Particle Physics (MIAPP) of the DFG cluster of excellence 'Origin and Structure of the Universe.' GK thanks MIAPP for hospitality while part of this work was completed.

arXiv:1509.01598 [hep-ph].

- [9] Y. Hochberg, Y. Zhao, and K. M. Zurek, Phys. Rev. Lett. 116 , 011301 (2016), arXiv:1504.07237 [hep-ph].
- [10] Y. Hochberg, M. Pyle, Y. Zhao, and K. M. Zurek, JHEP 08 , 057 (2016), arXiv:1512.04533 [hep-ph].
- [11] G. Cavoto, E. N. M. Cirillo, F. Cocina, J. Ferretti, and A. D. Polosa, Eur. Phys. J. C76 , 349 (2016), arXiv:1602.03216 [physics.ins-det].
- [12] S. Derenzo, R. Essig, A. Massari, A. Soto, and T.-T. Yu, Phys. Rev. D96 , 016026 (2017), arXiv:1607.01009 [hep-ph].
- [13] R. Essig, J. Mardon, O. Slone, and T. Volansky, Phys. Rev. D95 , 056011 (2017), arXiv:1608.02940 [hep-ph].
- [14] Y. Hochberg, Y. Kahn, M. Lisanti, C. G. Tully, and K. M. Zurek, Phys. Lett. B772 , 239 (2017), arXiv:1606.08849 [hep-ph].
- [15] C. Kouvaris and J. Pradler, Phys. Rev. Lett. 118 , 031803 (2017), arXiv:1607.01789 [hep-ph].
- [16] R. Essig, T. Volansky, and T.-T. Yu, Phys. Rev. D96 , 043017 (2017), arXiv:1703.00910 [hep-ph].

- [17] G. Cavoto, F. Luchetta, and A. D. Polosa, Phys. Lett. B776 , 338 (2018), arXiv:1706.02487 [hep-ph].
- [18] S. Fichet, Phys. Rev. Lett. 120 , 131801 (2018), arXiv:1705.10331 [hep-ph].
- [19] J. Tiffenberg, M. Sofo-Haro, A. Drlica-Wagner, R. Essig, Y. Guardincerri, S. Holland, T. Volansky, and T.T. Yu (SENSEI), Phys. Rev. Lett. 119 , 131802 (2017), arXiv:1706.00028 [physics.ins-det].
- [20] Y. Hochberg, Y. Kahn, M. Lisanti, K. M. Zurek, A. G. Grushin, R. Ilan, S. M. Griffin, Z.-F. Liu, S. F. Weber, and J. B. Neaton, Phys. Rev. D97 , 015004 (2018), arXiv:1708.08929 [hep-ph].
- [21] R. K. Romani et al. , Appl. Phys. Lett. 112 , 043501 (2018), arXiv:1710.09335 [physics.ins-det].
- [22] P. Agnes et al. (DarkSide), Phys. Rev. Lett. 121 , 111303 (2018), arXiv:1802.06998 [astro-ph.CO].
- [23] M. Crisler, R. Essig, J. Estrada, G. Fernandez, J. Tiffenberg, M. Sofo haro, T. Volansky, and T.-T. Yu (SENSEI), Phys. Rev. Lett. 121 , 061803 (2018), arXiv:1804.00088 [hep-ex].
- [24] R. Agnese et al. (SuperCDMS), Phys. Rev. Lett. 121 , 051301 (2018), [erratum: Phys. Rev. Lett.122,no.6,069901(2019)], arXiv:1804.10697 [hep-ex].
- [25] R. Essig, M. Sholapurkar, and T.-T. Yu, Phys. Rev. D97 , 095029 (2018), arXiv:1801.10159 [hep-ph].
- [26] O. Abramoff et al. (SENSEI), Phys. Rev. Lett. 122 , 161801 (2019), arXiv:1901.10478 [hep-ex].
- [27] T. Emken, R. Essig, C. Kouvaris, and M. Sholapurkar, (2019), arXiv:1905.06348 [hep-ph].
- [28] A. Migdal, J. Phys. (USSR) (1941).
- [29] J. D. Vergados and H. Ejiri, Phys. Lett. B606 , 313 (2005), arXiv:hep-ph/0401151 [hep-ph].
- [30] C. C. Moustakidis, J. D. Vergados, and H. Ejiri, Nucl. Phys. B727 , 406 (2005), arXiv:hep-ph/0507123 [hep-ph].
- [31] R. Bernabei et al. , Int. J. Mod. Phys. A22 , 3155 (2007), arXiv:0706.1421 [astro-ph].
- [32] M. Ibe, W. Nakano, Y. Shoji, and K. Suzuki, JHEP 03 , 194 (2018), arXiv:1707.07258 [hep-ph].
- [33] M. J. Dolan, F. Kahlhoefer, and C. McCabe, Phys. Rev. Lett. 121 , 101801 (2018), arXiv:1711.09906 [hep-ph].
- [34] N. F. Bell, J. B. Dent, J. L. Newstead, S. Sabharwale, and T. J. Weiler, (2019), arXiv:1905.00046 [hep-ph].
- [35] D. S. Akerib et al. (LUX), Phys. Rev. Lett. 122 , 131301 (2019), arXiv:1811.11241 [astro-ph.CO].
- [36] E. Armengaud et al. (EDELWEISS Collaboration), Phys. Rev. D 99 , 082003 (2019).
- [37] Z. Z. Liu et al. (CDEX), (2019), arXiv:1905.00354 [hepex].
- [38] E. Aprile et al. , (2019), arXiv:1907.12771 [hep-ex].
- [39] J. Angle et al. (XENON10), Phys. Rev. Lett. 107 , 051301 (2011), [Erratum: Phys. Rev. Lett.110,249901(2013)], arXiv:1104.3088 [astro-ph.CO].
- [40] E. Aprile et al. (XENON100), J. Phys. G41 , 035201 (2014), arXiv:1311.1088 [physics.ins-det].
- [41] E. Aprile et al. , (2019), arXiv:1907.11485 [hep-ex].
- [42] B. Holdom, Phys. Lett. B166 , 196 (1986).
- [43] L. Okun, Sov. Phys. JETP 56 , 502 (1982).
- [44] M. Bauer, P. Foldenauer, and J. Jaeckel, (2018), arXiv:1803.05466 [hep-ph].
- [45] G. Krnjaic, Phys. Rev. D94 , 073009 (2016), arXiv:1512.04119 [hep-ph].
- [46] S. Knapen, T. Lin, and K. M. Zurek, Phys. Rev. D96 , 115021 (2017), arXiv:1709.07882 [hep-ph].
- [47] J. Angle et al. (XENON10 Collaboration), Phys. Rev.
- Lett. 107 , 051301 (2011).
- [48] R. Essig, A. Manalaysay, J. Mardon, P. Sorensen, and T. Volansky, Phys. Rev. Lett. 109 , 021301 (2012).
- [49] E. Aprile et al. (XENON Collaboration), Phys. Rev. D 94 , 092001 (2016).
- [50] E. Aprile et al. , (2019), arXiv:1906.04717 [physics.insdet].
- [51] C. E. Dahl, Ph.D. thesis, Princeton University (2009).
- [52] D. S. Akerib et al. (LUX), Phys. Rev. D 93 , 072009 (2016), arXiv:1512.03133 [physics.ins-det].
- [53] D. S. Akerib et al. (LUX), (2016), arXiv:1608.05381 [physics.ins-det].
- [54] E. Aprile et al. (XENON), Phys. Rev. D 97 , 092007 (2018), arXiv:1709.10149 [astro-ph.IM].
- [55] P. Sorensen and C. E. Dahl, Phys. Rev. D 83 , 063501 (2011).
- [56] D. S. Akerib et al. (LUX), Phys. Rev. D96 , 112011 (2017), arXiv:1709.00800 [physics.ins-det].
- [57] T. Doke, A. Hitachi, J. Kikuchi, K. Masuda, H. Okada, and E. Shibamura, Japanese Journal of Applied Physics 41 , 1538 (2002).
- [58] E. Aprile, K. L. Giboni, P. Majewski, K. Ni, and M. Yamashita, Phys. Rev. B 76 , 014115 (2007).
- [59] A. K. E. Gushchin and I. Obodovskii, Sov. Phys. JETP 55 , 860 (1982).
- [60] J. Xu, S. Pereverzev, B. Lenardo, J. Kingston, D. Naim, A. Bernstein, K. Kazkaz, and M. Tripathi, Phys. Rev. D 99 , 103024 (2019), arXiv:1904.02885 [physics.ins-det].
- [61] C. F. Bunge, J. A. Barrientos, and A. V. Bunge, Atomic Data and Nuclear Data Tables 53 , 113 (1993).
- [62] F. Brandi, I. Velchev, W. Hogervorst, and W. Ubachs, Physical Review A 64 , 032505 (2001).
- [63] M. Cardona and L. Ley, Photoemission in Solids I: General Principles (Springer-Verlag, Berlin, 1978).
- [64] M. F. Gu, Canadian Journal of Physics 86 , 675 (2008).
- [65] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, Phys. Rev. Lett. 115 , 251301 (2015), arXiv:1505.00011 [hepph].
- [66] J. L. Feng and J. Smolinsky, Phys. Rev. D96 , 095022 (2017), arXiv:1707.03835 [hep-ph].
- [67] N. Aghanim et al. (Planck), (2018), arXiv:1807.06209 [astro-ph.CO].
- [68] E. Kuflik, M. Perelstein, N. R.-L. Lorier, and Y.-D. Tsai, Phys. Rev. Lett. 116 , 221302 (2016), arXiv:1512.04545 [hep-ph].
- [69] Y. Hochberg, E. Kuflik, T. Volansky, and J. G. Wacker, Phys. Rev. Lett. 113 , 171301 (2014), arXiv:1402.5143 [hep-ph].
- [70] P. deNiverville, M. Pospelov, and A. Ritz, Phys. Rev. D84 , 075020 (2011), arXiv:1107.4580 [hep-ph].
- [71] J. D. Bjorken, S. Ecklund, W. R. Nelson, A. Abashian, C. Church, B. Lu, L. W. Mo, T. A. Nunamaker, and P. Rassmann, Phys. Rev. D38 , 3375 (1988).
- [72] B. Batell, R. Essig, and Z. Surujon, Phys. Rev. Lett. 113 , 171802 (2014), arXiv:1406.2698 [hep-ph].
- [73] R. Dharmapalan et al. (MiniBooNE), (2012), arXiv:1211.2258 [hep-ex].
- [74] A. A. Aguilar-Arevalo et al. (MiniBooNE DM), Phys. Rev. D98 , 112004 (2018), arXiv:1807.06137 [hep-ex].
- [75] G. Angloher et al. (CRESST), Eur. Phys. J. C76 , 25 (2016), arXiv:1509.01515 [astro-ph.CO].
- [76] J. P. Lees et al. (BaBar), Phys. Rev. Lett. 119 , 131804 (2017), arXiv:1702.03327 [hep-ex].
- [77] E. Izaguirre, G. Krnjaic, P. Schuster, and N. Toro, Phys.

- Rev. D88 , 114015 (2013), arXiv:1307.6554 [hep-ph]. [78] R. Essig, J. Mardon, M. Papucci, T. Volansky, and Y.-M. Zhong, JHEP 11 , 167 (2013), arXiv:1309.5084 [hep-ph]. [79] M. K. Pandey, L. Singh, C.-P. Wu, J.-W. Chen, H.-C. (2018),
- Chi, C.-C. Hsieh, C. P. Liu, and H. T. Wong, arXiv:1812.11759 [hep-ph].
- [80] A. Aguilar-Arevalo et al. (DAMIC), (2019), arXiv:1907.12628 [astro-ph.CO].
- [81] M. Settimo (DAMIC), in 53rd Rencontres de Moriond on QCD and High Energy Interactions (Moriond QCD 2018) La Thuile, Italy, March 17-24, 2018 (2018) arXiv:1805.10001 [astro-ph.IM].
Daniel Baxter Yonatan Kahn Gordan Krnjaic
Electron Ionization via Dark Matter-Electron Scattering and the Migdal Effect
There are currently several existing and proposed experiments designed to
probe sub-GeV dark matter (DM) using electron ionization in various materials.
The projected signal rates for these experiments assume that this ionization
yield arises only from DM scattering directly off electron targets, ignoring
secondary ionization contributions from DM scattering off nuclear targets. We
investigate the validity of this assumption and show that if sub-GeV DM couples
with comparable strength to both protons and electrons, as would be the case
for a dark photon mediator, the ionization signal from atomic scattering via
the Migdal effect scales with the atomic number $Z$ and 3-momentum transfer
$\mathbf{q}$ as $Z^2 \mathbf{q}^2$. The result is that the Migdal effect is
always subdominant to electron scattering when the mediator is light, but that
Migdal-induced ionization can dominate over electron scattering for heavy
mediators and DM masses in the hundreds of MeV range. We put these two
ionization processes on identical theoretical footing, address some theoretical
uncertainties in the choice of atomic wavefunctions used to compute rates, and
discuss the implications for DM scenarios where the Migdal process dominates,
including for XENON10, XENON100, and the recent XENON1T results on light DM
scattering.