(Not as) Big as a Barn: Upper Bounds on Dark Matter-Nucleus Cross Sections

Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $\sigma_{\chi A} \propto A^2 \mu_A^2 \sigma_{\chi N}\propto A^4 \sigma_{\chi N}$, where the $A^2$ comes from coherence and the $\mu_A^2\simeq A^2 m_N^2$ from kinematics for $m_\chi\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $\sigma_{\chi A} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $\sigma_{\chi N}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $\sigma_{\chi N}>10^{-32}\;\text{cm}^2$ with $m_\chi\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.


I. INTRODUCTION
The nature of dark matter is one of the most pressing problems in both fundamental physics and cosmology. Decades of observations indicate that dark matter makes up the vast majority of matter in our Universe, yet increasingly advanced experiments have yet to determine its physical nature. Once discovered, the particle properties of dark matter will be a guidepost to physics beyond the Standard Model as well as to an improved understanding of galaxies and cosmic structure [1][2][3][4][5][6][7].
For these searches, the two most common benchmarks for the performance of dark matter detection experiments are the dark matter self-annihilation cross section and the spin-independent dark matter-nucleon elastic scattering * digman.12@osu.edu Claimed constraints on the spin-independent dark matter-nucleon cross section [56,[59][60][61][62]. Those from cosmology directly probe scattering with protons, but all others here are based on scattering with nuclei, and thus require the use of "model-independent" scaling relations. Below, we show that assumptions used to derive these results are invalid over most of the plane.
cross section, the simplest case (for more general treatments, see, e.g., Refs. [63][64][65][66][67][68][69]). These benchmarks allow constraints set by different experiments to be scaled to each other. Here, we focus on spin-independent elas-tic dark matter-nucleus scattering for dark matter with m χ 1 GeV. For generality, we do not require that dark matter be a thermal relic. Most direct-detection searches focus on pushing sensitivity to small cross sections, but it is also important to consider constraints on large cross sections [59][60][61][62][70][71][72][73][74][75][76][77][78][79]. Direct-detection experiments are typically located beneath the atmosphere, rock, and detector shielding, such that dark matter with too large of a cross section loses too much energy above the detector. Energy loss in the detector overburden may open a window where strongly interacting dark matter is allowed [70]. Figure 1 summarizes prior claimed constraints. The "IMP+IMAX+SKYLAB" region is based on atmospheric and space-based detectors and is dashed because the results are commonly cited but are not based on detailed analyses in peer-reviewed papers [59]. The X-ray Quantum Calorimeter (XQC) experiment is rocket based [80]. There are several similar proposed XQC regions [59,60,62]; we adopt that of Ref. [62]. The "Underground Detectors" region is taken directly from the summary plot in Ref. [61]. For the "Cosmology" region, we plot the strongest constraint that depends only on dark matterproton scattering [56] (including helium would make the constraints somewhat stronger [48]). The details of which constraints are plotted do not affect our conclusions.
Direct-detection searches for spin-independent interactions benefit from an essentially model-independent A 2 coherent enhancement, as well as a kinematic factor of µ 2 A , such that σ χA is related to the dark matternucleon elastic scattering cross section σ χN by σ χA ∝ A 2 µ 2 A σ χN . For m χ m A , the dark matter-nucleus reduced mass µ A Am N , such that the scaling becomes σ χA ∝ A 4 σ χN . This straightforward scaling allows constraints on dark matter-nucleus scattering to be related to each other and to the cross section on nucleons. This scaling is model independent in the sense that it is independent of the detailed shape of the potential. In Fig. 1, all constraints except the one labeled "Cosmology" deal with nuclear targets with A > 1, and hence assume this scaling relation.
How large of cross sections are allowed before the defining assumptions are violated? Here we systematically calculate the theoretical upper limits on dark matternucleon cross sections. We show that most of the parameter space of Fig. 1 is beyond the point where the simple scaling relations above are valid, or where pointlike dark matter is even allowed. Our results are based first on generic considerations of theoretically allowed cross sections for short-range interactions with nuclei, and second on classes of models where we consider light mediators as a mechanism to obtain large cross sections. As far as we are aware, this is the first systematic exploration of these issues for dark matter-nucleus scattering (for related considerations in strongly self-interacting dark matter sectors, see, e.g., Ref. [81]). Our results will require the reinterpretation of a large and varied body of work, e.g. Refs. [43, 46, 48, 49, 59-62, 70-79, 82-92].
In Sec. II, we review the nonrelativistic scattering theory used to obtain the model-independent scaling relations. In Sec. III, we examine the various ways that scaling relations can break down for contact interactions. In Sec. IV, we examine the possibility of achieving a larger cross section with a light mediator in light of present constraints on light mediators. In Sec. V, we briefly discuss the possibility that dark matter itself could have a nonzero physical extent. In Sec. VI, we discuss the implications for existing constraints and future experiments. Finally, we summarize our results and the outlook for future work in Sec. VII.

II. DARK MATTER SCATTERING THEORY
We briefly review the basic nonrelativistic scattering theory required to derive the model-independent scaling relation for the spin-independent elastic scattering cross section. We also discuss how some of the key assumptions may break down. Throughout, we set = c = 1.
A. Overview of basic assumptions 1. Single particle: Dark matter is primarily a single unknown particle. The number density of dark matter is then determined only by its mass and the local dark matter density.

2.
Pointlike: Dark matter is a pointlike particle with no excitation spectrum.
3. Electrically neutral: It is typically assumed that dark matter is electrically neutral. Millicharged dark matter has different dynamics and is too strongly constrained to produce large cross sections. 4. Equal coupling to all nucleons: For simplicity, we assume that dark matter has equal coupling to both protons and neutrons, although this assumption is not essential to any of our conclusions.

Local:
The interaction is assumed to be local, 6. Energy-independent potential: The potential for the interaction is assumed to be energy independent, such that the cross section for the interaction is also energy independent up to a form factor. For a spin-independent interaction, the potential must also be independent of the incident angular momentum l.

7.
Elastic: For laboratory experiments, dark matternucleus scattering is assumed to occur at typical Milky Way virial velocities, v ∼ 10 −3 c. Typical recoil energies of O(1 keV) are not sufficient to produce Standard-Model particles, or to excite internal degrees of freedom of nuclei. Therefore, elastic scattering is the dominant interaction channel. In any case, all physical scattering processes have at least some elastic component [93]. 8. Coherence: Closely related to the assumption of purely elastic scattering is the assumption of coherence. For coherence to hold, it must be a good approximation to treat the dark matter as interacting with the nucleus as a whole, rather than with individual nucleons. Coherence is typically a good approximation provided the momentum transfer q is insufficient to excite internal degrees of freedom in a nucleus, which is true provided 1/q is large compared to the characteristic nuclear radius r A , qr A 1 [94]. The breakdown of coherence can be parametrized by including a momentum-dependent form factor in the differential cross section. 9. No bulk effects: The scattering should be well approximated as being with a single nucleus, such that initial and final state effects in the bulk medium can be ignored. This approximation is good as long as the characteristic momentum transfer q is large compared to the characteristic interatomic spacing, which is typical.
The rest of this paper deals with the failure of the following additional assumptions: 10. S-wave scattering For s-wave (l = 0) scattering, the scattering is isotropic in the center-ofmomentum frame. As shown in Sec. II C 2, assuming l = 0 is required to derive the modelindependent A 2 µ 2 A scaling relation. However, real interactions may deviate significantly from isotropic scattering, and we do not require l = 0 in this analysis. 11. Weak Interaction: For A 2 µ 2 A scaling to hold, the interaction must be weak enough for the Born approximation to hold. We discuss this assumption in Sec. II C 1.
To be detectable, dark matter must have some kind of interaction with ordinary matter in a detector, written here as a potential V (r). We specialize to spinindependent interactions, and restrict our analysis to spherically symmetric potentials V (r) = V (r) that fall off faster than r −1 as r → ∞. In the center-of-momentum frame, the time-independent Schrödinger equation giving the evolution of a nonrelativistic two-particle system with wave function ψ(r) and reduced mass µ is given as As shown in Appendix A, far from the potential the solution of Eq. (1) may be written as where f (k i , k f ) = f (k, θ) is the scattering amplitude, ψ 0 (r) ≡ (2π) −3/2 e iki·r , and k i ≡ kẑ and k f are the initial and final dark matter momenta, respectively. From the scattering amplitude, we obtain the differential cross section and the total elastic scattering cross section: If the scattering is isotropic, f (k, θ) = f (k), Eq. (4) is proportional to the rate of detectable scattering events in a detector. However, to be detectable, a collision must deposit sufficient energy into the detector. If the scattering angle is peaked close to θ = 0, very little momentum is transferred, and hence insufficient energy is deposited in the detector. Therefore, it is sometimes more useful to weight the integral in Eq. (4) by the momentum transfer to obtain the momentum-transfer cross section, For isotropic scattering, σ mt χA = σ χA . For a potential with characteristic radius r A , isotropic scattering is generically a good approximation at low energies, kr A 1, as discussed further in Sec. II C 2. Forward scattering is a major concern for light mediators (Sec. IV); in the Coulomb scattering limit where the mediator mass m φ → 0, then σ χA → ∞, while σ mt χA remains finite.

C. Derivation of model-independent scaling
Now we discuss approximation methods for f (k, θ). The two approaches we consider here are the Born approximation and the partial wave expansion. Both approaches allow us to derive the σ χA = A 2 µ 2 A /µ 2 N σ χN scaling with nuclear mass number A. The reduced masses are defined as where the mass of a nucleus with mass number A is related to the mass of a single nucleon m N using m A = Am N . We begin with the Born approximation because it is the simple and familiar derivation. Because the partial wave expansion is valid even when the Born approximation fails, it allows us to more concretely show the behavior at large scattering cross sections.

Born approximation
Inspecting Eq. (2), a natural first approach to obtaining f (k, θ) is to solve for ψ(r) by iteration, which is the Born approximation, as demonstrated in Appendix B. The first Born approximation to f (k, θ) is simply the Fourier transform of the potential: where q = |q| = 2k sin θ/2 is the momentum transfer. Now, assume that the potential has some maximum radius r A , and we have low energy scattering, kr A 1. Then we can approximate sin(qr ) ≈ qr and integrate only up to the maximum radius r A : Equation (7) is a remarkable result. If the required approximations are valid, f (1) (k, θ) depends only on the volume integral of the potential; it contains no information at all about the shape. Provided the volume integral of the potential is proportional to the nuclear mass number A, we have the scaling Plugging into Eq. (4): which can be recast more precisely in terms of the dark matter-nucleon reduced mass µ N and scattering cross section σ (1) χN : Eq. (10) is the famous model-independent scaling relation for the spin-independent elastic scattering cross section. Provided the potential falls off faster than 1/r, this scaling relation is generally a good approximation at sufficiently low energies, so long as the first Born approximation reasonably approximates f (k, θ). However, we must examine when the first Born approximation fails.
We discuss the validity of the Born approximation in Appendix B. A useful condition for the validity of the first Born approximation is [99] We can simplify Eq. (11) using our assumption of a maximum range r A and kr A 1: Equation (12) is equivalent to the statement that the potential is much too weak to form a bound state even if V (r) was purely attractive [101]. While Eq. (7) is a volume integral, Eq. (12) is an area integral of the potential. Therefore, the first Born approximation is valid when some potential-weighted effective area is small. The effective area in question is, in fact, the elastic scattering cross section, as shown for a contact interaction in Sec. III.

Partial wave expansion
To investigate what happens when the Born approximation fails, the first step is to expand the scattered wave function in terms of Legendre polynomials and calculate the phase shift of each contribution. The phase shifts may be found by numerically integrating the Schrödinger equation, as described in Appendix C. The elastic scattering cross section may be written in terms of the phase shifts δ l (k): The momentum-transfer cross section in Eq. (5) may also be written in terms of partial wave phase shifts: The mathematical decomposition in terms of partial waves is valid even beyond interactions that can be described in nonrelativistic potential scattering theory. However, when the number of partial waves becomes too large, it may be impractical to compute the phase shifts individually, and semiclassical approximations become useful [93,98]. Physically, the sum over partial waves in Eq. (13) is equivalent to the classical operation of averaging over all possible impact factors b = l(l + 1)/k [93,100]. Classically, for a potential with maximum range r A , there would be no collisions for b > r A . Therefore, a useful approximate upper limit on the highest partial wave that can meaningfully contribute to the sum in nonrelativistic quantum scattering is l max ≈ kr A , and contributions from higher l > l max fall off quickly [93]. Our derivation of the model-independent scaling of Eq. (10) in the Born approximation assumes kr A 1, which is equivalent to saying only the l = 0 (s-wave) term contributes a nonvanishing phase shift.
Using the same iterative procedure as for the Born approximation, we can obtain the model-independent form of the s-wave phase shift [100]: where the required approximation is δ 0 (k) 1. Plugging Eq. (15) into Eq. (13) and again using δ 0 (k) 1, we obtain precisely the same expression for the scattering amplitude we obtained for the Born approximation in Eq. (7), as expected. Again, the requirement δ 0 (k) 1 places an upper bound on the maximum σ χA where the relation can apply.
If the cross section were instead the maximum allowed by unitarity, δ 0 (k) = π/2, we would obtain which decreases as 1/A 2 with increasing A, assuming k ∝ A, rather than increasing as A 4 . Higher partial waves necessarily scale differently with k Higher δ l (k) also contain information about the shape of the potential. Therefore, we do not expect any special model-independent scaling when higher partial waves contribute.

III. CONTACT INTERACTIONS
In this section, we consider the limits on cross sections that can be obtained through a contact interaction, and how the scaling relations break down. A contact interaction is useful as an illustrative case because we do not need to consider the specific mechanism that produces the interaction.

A. Contact interaction with Born approximation
As a simple case, we consider a contact interaction with a nucleus, as could be produced by a heavy mediator. We roughly approximate the nuclear charge density as having a top-hat shape with radius r A : We assume the maximum charge density is roughly independent of atomic mass number A, such that We use this toy model with a sharp cutoff because both the Born approximation and the partial wave phase shift δ l (k) can be found analytically. The effect of using a more realistic charge distribution is discussed in Sec. III C 1.
Fourier transforming Eq. (17) using Eq. (6) gives The total elastic scattering cross section in the first Born approximation is then: In the limit kr A 1, Eq. (19) becomes For scattering with a nucleon, Eq. (20) would become σ we recover the required scaling relation of Eq. (10): In the kr A 1 limit, the condition for validity of the first Born approximation in Eq. (12) is simply Comparing to Eq. (20), we can rewrite the condition Eq. (22) as: Equation (23) has a simple physical interpretation. The first Born approximation is only applicable for elastic scattering cross sections much smaller than the geometric cross section of the nucleus. Using r N ≈ 1.2 fm in Eq. (23), the Born approximation result only applies for σ (1) χN 10 −25 cm 2 . Going to higher orders in the Born approximation does not unlock cross sections significantly exceeding the geometric limit of the potential. For σ (1) χA > 16 9 πr 2 A , the Born series is not even guaranteed to converge for all energies [93].
However, it may still be possible to obtain a meaningful cross section in regimes where the Born approximation fails using partial wave analysis. We explore this technique below.

B. Contact interaction with partial waves
For r < r A in Eq. (17), the radial wave function decomposed in partial waves has an analytic solution in terms of partial waves, u l (r) = C l rj l (k r), where k ≡ k 2 − 2µ A V 0 1/2 could be either pure real or pure imaginary. First, we consider the s-wave cross section, with l = 0. Expanding in the limit where V 0 and k are small, The corresponding s-wave cross section is which is identical to Eq. (20) to lowest order in |V 0 |, as anticipated in Sec. II C 2. Now we can see how the model-independent A scaling fails as the coupling strength gets stronger. The secondorder term in Eq. (25) scales ∝ µ 3 A /µ 3 N A 8/3 , and either reduces the cross section for a repulsive potential V 0 > 0 or increases it for an attractive one. For V 0 large enough for the second-order and higher corrections to matter, there is, therefore, no model-independent scaling with A because σ l=0 χA depends on the details of the potential. The breakdown of the scaling as a function of A is shown in Fig. 2. Once details of the potential begin to matter, corrections from a more realistic charge distribution would also become important, as discussed in Sec. III C 1.
To illustrate further, if we instead considered the strong coupling limit, and so the repulsive cross section completely saturates at 4 times the geometric cross section. The saturation is plotted as a function of |V 0 | in Fig. 3. Physically, the well has simply become an impenetrable hard sphere with a fixed radius. Therefore, we have obtained a physical maximum cross section for the repulsive contact interaction at small k which is only 9/4 times larger than the maximum we obtained using the first Born approximation.

Higher partial waves
While the s-wave cross section is limited to the geometric cross section, it is natural to wonder if the contributions from higher partial waves could allow a larger total cross section. For kr A 1, we can approximate semiclassically l max ≈ kr A , as described in Sec. II C 2. Of course, in the quantum case, it is possible for higher partial waves to contribute, but their contributions fall off rapidly for l > l max . Assuming a sharp cutoff is adequate for deriving an approximate maximum physically achievable cross section. The maximum possible cross section is achieved by saturating partial wave unitarity, i.e., taking δ l≤lmax (k) = π/2: Now, we see that the saturation at approximately the geometric cross section, found for kr A 1 in Eq. (28), also holds for kr A 1. In fact, for a very strong repulsive contact interaction, the phase shifts for l ≤ l max approach δ l≤lmax (k) ≈ lπ/2 − kr A [93], such that σ χA ≈ 2πr 2 A . Therefore, a repulsive hard sphere almost saturates the unitarity limit of Eq. (31).
Including higher partial waves is therefore not a useful way of increasing the cross section, because the potential remains limited by the characteristic radius r A . Figure 4 shows the breakdown of the A 4 scaling for several example nuclei, fully taking into account contributions from higher partial waves. Direct-detection constraints for underground detectors are affected by the breakdown of scaling at the O(1) level for σ χN 10 −32 cm 2 .

C. Attractive resonances
A final possible approach would be to saturate unitarity at δ l = π/2 while kr A 1, such that can become large. The limit δ l = π/2 at kr A 1 can be achieved through resonances, which occur when an attractive potential becomes strong enough to support a bound state.
In reality, the resonant scattering cross section would achieve large values only for a narrow range of k relative to the incident dark matter velocity distribution [93,102]. Because k = µ A v, the resonances are generically at different incident dark matter velocities for different elements, which guarantees that there are not any useful modelindependent scaling relations relating the observed cross sections for strongly attractive potentials between different target materials.
In Fig. 2, we show the behavior as a function of A for two different values of V 0 . When resonances are possible, the scaling with A need not be monotonic. The behavior is fairly complex for even the simple rectangular well of Eq. (17). Realistic scaling is likely to be even more complicated, because the nuclear charge distribution changes as a function of A. Even two different nuclei with the same A but different atomic numbers could have different charge distributions, and hence different resonant cross sections. Scaling with A for strong attractive couplings is therefore highly model-dependent. In Fig. 3, we show the saturation of the s-wave cross section as a function of the coupling strength, as well as the resonant behavior which occurs once the potential becomes strong enough to support a quasibound state. The resonances for A = 4 are fairly narrow, as for a nucleon the scattering is still well approximated in the low-k limit. However, for A = 131, the low-k limit is no longer a good approximation, and resonances are broadened to the point that they do not We include partial waves up to lmax = 8, which is sufficient to converge σχA to ∼ 10 −16 precision. Attractive and repulsive interactions scale similarly, although the scaling deviates from A 4 at high A due to form-factor suppression, accounted for here by including the contributions from higher partial waves. Bottom: Same as above, but with |V0| = 1.18 × 10 −3 GeV, which corresponds to the 'scaling relations unreliable for A > 12' line in Fig. 6. Repulsive and attractive interactions no longer scale the same way, and both saturate close to 4πR 2 A . The attractive potential shows resonances with A, which are sensitive to the specific choice of potential. For cross sections approaching the geometric cross section, any scaling with A is highly model dependent.
significantly increase the scattering cross section. Additionally, there are many more resonances, from multiple partial waves. Note we have not implemented any velocity dispersion for this plot; the spreading of resonances is entirely due to broadening of peaks and overlapping contributions from multiple partial waves at finite k. Applying a realistic dark matter velocity distribution would smooth the peaks. For heavy nuclei with multiple naturally occurring isotopes (e.g. xenon), averaging over a distribution of isotopes would smooth the peaks even further.
Overall, even if a carefully tuned resonance could achieve a large cross section for a single light nucleus, other nuclei would not necessarily have correspondingly large cross sections. Scaling relations between specific The cross sections are computed using the analytic partial wave results. For attractive potentials, once |V0| becomes large enough to support quasibound states, resonances can increase the cross section by several orders of magnitude, but only in a narrow range. Bottom: Same as above, but with A = 131 (xenon). A larger number of partial waves contribute due to the larger k ∝ A. There are many resonances, but they are not large enough to meaningfully increase the cross section above the geometric limit. Additionally, the resonances are not at the same values of |V0|, which prevents resonances from achieving a large cross section which scales predictably with A, as shown in Fig. 2. nuclear cross sections would also be highly model dependent, such that constraints from different types of nuclei would be difficult to compare because the full resonance structure would not be known. For example, using more realistic charge distributions, such as an exponential potential for A = 4 and a Woods-Saxon potential for A = 131 [103] would shift the positions of the resonances somewhat.

More realistic charge distributions
The rectangular barrier potential in Eq. (17) is a toy model. Realistic nuclear charge distributions have a smooth cutoff, and an exponential tail to larger radii [103], as in a Woods-Saxon potential. Because there is not a sharp cutoff, allowing the interaction strength to be arbitrarily large would cause the potential to grow logarithmically with |V 0 |. However, there are limits to how strong |V 0 | can be. For |V 0 | 10 −1 GeV, QCD corrections break the simple nonrelativistic contact interaction picture. For |V 0 | 2 GeV, the interaction may be strong enough to pull proton-antiproton pairs out of the vacuum.
We have verified by numerically computing partial wave amplitudes that, for |V 0 | 10 −1 GeV, using a Woods-Saxon potential increases the maximum σ χA by a factor of 10. By definition, an increase in the computed cross section due to a different potential can only appear for |V 0 | strong enough that the model-independent form of the cross section has already significantly broken down. Therefore, using a more realistic charge distribution cannot significantly change our conclusion that the cross section for a contact interaction cannot be much larger than the geometric cross section of the nucleus.
Realistic potentials are also not perfectly spherically symmetric; however, the same basic picture of geometrical limitations still applies, and resonances are still possible.

D. Beyond contact interactions
We have established that a contact interaction with a nucleus cannot achieve cross sections much larger than the geometrical cross section of the nucleus. The case where large cross sections might be achieved, a strongly attractive potential, produces resonances that are sensitive to the detailed structure of the potential and is far too model dependent to possess any simple scaling relation relating the cross sections at different A. To circumvent these problems, we need an interaction with a larger characteristic range. One possible way to achieve a larger characteristic range is to insert a light mediator for the interaction, as discussed in Sec. IV. Another possibility is composite dark matter with an intrinsic radius, discussed in Sec. V.

IV. LIGHT MEDIATOR
A simple approach to achieving a larger characteristic radius is to insert a light mediator, of mass m φ = 1/r φ , which generically results in a potential of the form where λ χ and λ A = Aλ N are the coupling strengths of the particle φ to the dark matter and nucleus respectively. To achieve cross sections much larger than a nucleus, we should have r φ 1 fm. The dark matter and target nucleus are distinguishable particles, so the Yukawa potential can be either attractive or repulsive. We assume the mediator is a scalar, although the general form of the potential would be similar for other light mediator candidates. The scattering amplitude of Eq. (33) is easily calculated using Eq. (6): which gives the total elastic scattering cross section: Because the characteristic radius is larger than the geometric radius of a nucleon, Eq. (35) can, in principle, achieve larger cross sections within the domain of validity of the Born approximation than a contact interaction could. Because k ∝ µ A , the scaling with A is now more complicated due to the k 2 r 2 φ term in the denominator. Assuming m χ m A such that µ A ≈ Am N , we have two limits: For direct detection, k N ∼ 0.005 fm −1 is set by Milky Way halo velocities v 10 −3 c and the mass of a single nucleon, k A Ak N . For 131 Xe, k A ≈ 0.7 fm −1 , such that k A r φ > 1 occurs for r φ 1.4 fm. Therefore, A 4 scaling is, at best, marginal for heavy nuclei for any r φ that could conceivably produce a cross section σ χA 10 −25 cm 2 . Additionally, as established for a general potential with a characteristic radius r φ in Sec. II C 1, increasing the coupling strengths λ χ or λ N at fixed r φ causes the Born approximation, and therefore the A 4 scaling, to fail before cross sections larger than the geometric cross section are achieved. Therefore, for a light mediator, the A 4 scaling is only possible if σ (1) χA 10 −25 cm 2 . The failure of A 4 scaling occurs without even considering the constraints on the existence of light mediators discussed in Sec. IV B. The A 4 scaling is preserved to somewhat larger cross sections than for the contact interaction shown in Fig. 2.

A. Momentum-transfer cross section
In fact, even the A 2 scaling is too optimistic for the detectable momentum transfer in a detector with high A. Inspection of Eq. (34) shows that for kr φ 1, the scattering becomes strongly peaked at θ = 0. Therefore, it is more useful to consider the momentum-transfer cross section, Eq. (5). Using the Born approximation, we can calculate Eq. (5) analytically for the Yukawa potential: For kr φ 1, Eq. (37) simplifies to σ mt,(1) χA ≈ σ (1) χA as expected for isotropic scattering. However, for kr φ 1, we have: Equation (38) grows only ∝ log(A), such that, for a fixed total detector mass, the total energy deposited in the detector would be larger for nuclei with smaller A. Direct detection experiments that focus on protons and other light nuclei, such as Refs. [32,[106][107][108], may therefore be effective ways of constraining the landscape for modeldependent direct detection.

B. Existing limits on light mediators
If there were no other constraints on r φ or λ A , Eq. (38) would allow σ For m φ < 10 −9 GeV, fifth-force constraints become many orders of magnitude stronger and dominate other constraints [105]. Bottom: Same as above, but for the momentum-transfer cross section. The largest cross sections achieved in unconstrained regions are σ mt χN 10 −32 cm 2 . such a particle. Figure 5 shows the maximum achievable σ χN and σ mt χN for a repulsive Yukawa potential, conservatively using the perturbativity limit λ χ = 4π, µ N ≈ m p , and k = 0.005 fm −1 . When Eq. (12) is > 10 −4 , Fig. 5 uses the results from a numerical partial wave expansion with adaptive l max , rather than the Born approximation. In practice, the Born approximation is adequate in the entire unconstrained region.
Including all such constraints, we have σ χN 10 −27 cm 2 and σ mt χN 10 −32 cm 2 . Constraints that rely on lower relative velocities, such as the cosmological constraints discussed in Sec. VI C, could achieve larger cross sections, but their constraints would need to be scaled correctly to compare them to direct-detection constraints. The momentum-transfer cross section is restricted to be σ mt χN 10 −25 cm 2 even for velocities as low as 10 −6 c.
It is also possible to produce the light mediator in a collision [109]. Particle production is an inelastic scattering process and beyond the scope of this paper, but it could be another avenue to transfer momentum between dark matter and a detector.
The detailed constraints in Fig. 5 could be different for different types of mediators. For example, for a vector mediator, the BBN constraints would be stronger [104]. Other constraints might be weaker. However, σ χN 10 −27 cm 2 is already smaller than the geometric cross section of the nucleus, and circumventing individual constraints is unlikely to drastically change the overall conclusion that light mediators do not appear to be a promising approach to achieving large cross sections.

V. COMPOSITE DARK MATTER
Another mechanism for achieving a larger characteristic interaction radius is dark matter that is not a point particle, but instead has a finite physical extent [90,[110][111][112][113][114][115][116][117][118][119][120][121][122][123][124]. Such dark matter could take the form of a composite particle. Because such dark matter would likely require an entire dark sector, any conclusions about the largest possible cross section with composite dark matter would be intrinsically model dependent. Because the largest physical scale in the problem is no longer related to a property of the target nucleus, the cross section need not scale with A at all.
The actual scaling with A could only be determined by examining the particular model of composite dark matter. Additionally, achieving cross sections significantly larger than a nucleus with composite dark matter will always require kr dm 1 for typical Milky Way virial velocities, so constraints on composite dark matter will need to be computed with a specific dark matter form factor in mind. See Sec. VI C for a discussion of limits at the lower velocities relevant to cosmological limits. Analyses setting constraints on specific form factors at large cross sections should consider whether their specific choice of form factor can be achieved at the cross sections they are constraining in a physically realistic model. Therefore, limits on composite dark matter need to be calculated in specific models. Calculation of constraints on specific models of composite dark matter is left to future work. Figure 6 summarizes the approximate limits for the repulsive contact-interaction cross sections discussed in Sec. III. In the colored regions, the Born approximation begins to break down when the proton cross section is scaled to heavier nuclei, ultimately failing even for light nuclei. For pointlike dark matter with a contact interaction, cross sections much larger than the geometric cross section are completely forbidden. As discussed in Sec. IV, the limits for a light mediator are similarly below the geometric cross section. For m χ 10 16 GeV, the entire (small) exclusion region for underground detectors is affected by the failure of scaling relations. Future improvements to constraints could change the region where the entire exclusion region would fail. Additionally, all detectors' computed ceilings are affected by the breakdown of scaling relations.

A. Scaling constraints
In the regime where scaling relations are unreliable, it becomes more difficult to compare constraints between experiments. When the scaling relations fail, scaling constraints from different nuclei to the dark matter-nucleon cross section using the A 4 is no longer meaningful.
For both contact interactions and light mediators, as the cross section begins to saturate, the momentumtransfer cross section scales less than linearly with A. Therefore, at fixed total detector mass, there is more detectable momentum transfer into the detector for lighter target nuclei. The failure of the scaling relations also occurs at larger cross sections for smaller A. For example, a 12 C-based detector would be able to use the Born approximation, and therefore the scaling relations, up to about 3000 times larger dark matter-nucleon cross section than a 131 Xe-based detector. Therefore, robustly covering the large cross-section regime may be best accomplished by detectors using light nuclei, e.g. [32,[106][107][108].
One option is to simply not scale constraints at large cross sections. While with resonances it could be possible for heavy nuclei to have smaller cross sections than a single nucleon, broadening by the dark matter velocity dispersion may limit the effect of narrow resonances on the overall detectable signature. Therefore, a relatively conservative approach could be to plot the actual momentum-transfer cross-section constraints obtained from different nuclei on the same scale. In fact, if composite dark matter as discussed in Sec. V is indeed the most plausible strongly interacting dark matter candidate, disregarding scaling with A may be the most correct way of plotting constraints.

B. Detection ceilings
Now we briefly consider if the detection ceilings (i.e., the largest cross sections that can be probed by a given detector based on the detector's overburden) shown in Fig. 1 are preserved. In our simple model in Sec. III, cross sections simply saturate at 4 times the geometric cross section for heavier nuclei. Even if all nuclei in the detector overburden have an elastic scattering cross section equal to their geometric cross section, dark matter cannot be stopped by the overburden above some m χ [82].   Fig. 4. The breakdown is purely on theoretical grounds. Also shown is the maximum allowed momentum-transfer cross section for a m φ = 10 −4 GeV light mediator using the constraints shown in Fig. 5, coincidentally at a comparable scale. For mχ 10 4 GeV we have applied a conservative self-interaction constraint σχχ/mχ < 10 cm 2 /g [125]. For σχN 10 −25 cm 2 , no viable pointlike dark matter candidates exist.
Because all currently computed detector ceilings exist at cross sections where the breakdown of the A 4 scaling is severe, correctly calculated detector ceilings must be specialized to a specific model. For basic energy-independent cross-section scaling, the weakened ceilings likely lead to stronger direct-detection constraints for m χ 10 16 GeV. For such models, direct detection may even have exhausted the parameter space for cross sections up to the largest cross sections achievable with pointlike dark matter. For other dark matter form factors, the behavior around the ceiling could be more complicated.
Further work is required to make detailed adjustments to existing constraint contours to determine what dark matter parameter space has been constrained at large cross sections.

C. Dark matter-proton scattering constraints
Constraints that rely only on dark matter scattering directly with protons are not directly affected by the breakdown of scaling relations with A. These are primarily constraints from cosmology and astrophysics, al- though at least one laboratory experiment uses proton targets [108]. Astrophysics constraints (e.g. disk stability, stars, cosmic ray interactions, gas clouds, etc) are typically assumed to occur at galactic virial velocities, as for direct detection. Cosmology constraints, such as CMB and structure formation constraints, typically assume that collisions occur at smaller relative velocities. As shown in Fig. 6, the cross sections of interest for cosmological or astrophysical constraints are too large to be pointlike dark matter. Therefore, they should be reinterpreted as constraints on specific models of composite dark matter with a specified form factor, as discussed in Sec. V.
For cosmology constraints set at lower relative velocities, the suppression of the cross section by the form factor of dark matter is not as severe. One consequence is that it is possible to achieve somewhat larger cross sections for pointlike dark matter with a light mediator than those shown in Fig. 5, although even for velocities as low as v 0.3 km/s, existing constraints would still require σ mt χN 10 −25 cm 2 . However, invoking such a model would require additional caution, as direct-detection constraints would not be scaled correctly relative to the cosmology constraints, such that it would no longer be appropriate to plot cosmology and direct-detection constraints on the same axes, as done in Fig. 1.
Cosmological and astrophysical constraints set at masses m χ < 1 GeV, discussed in Sec. VI D, are at lower cross sections, and may still be meaningful constraints on pointlike dark matter. However, analyses at lower masses should either directly investigate how high their limits can be extrapolated, or make it much clearer that there are caveats in extrapolating their results to much larger masses.

D. Low-mass dark matter
Because for m χ 1 GeV, µ A m χ , low-mass dark matter constraints benefit only from a single factor of A 2 from coherence. Therefore, the loss of the A 2 scaling at large cross sections will be orders of magnitude less severe than the impact from the loss of A 4 scaling at larger masses. The momentum transfer is also smaller, so the loss of coherence due to an assumed form factor for the dark matter would be less severe. Contact interactions are still limited by the geometric size of the nucleus, but constraints on light mediators will become a function of m χ [104].
We leave a detailed assessment of the impact of our considerations at low mass to future work. However, we reiterate our caution that constraints set at low masses should carefully state the limitations on extrapolating their constraints to m χ 1 GeV.

VII. CONCLUSIONS
How do dark matter particles interact with matter? One of the most commonly considered cases to probe is the spin-independent interactions of m χ > 1 GeV pointlike dark matter with nuclei. In the literature, a vast array of constraints -based on astrophysical and cosmological tests, as well as direct-detection searches with a wide range of nuclei and overburdens -are all compared to each other in simple plots of the dark matter-nucleon cross section and dark matter mass. Comparing searches in this way requires the assumption of scaling relations, e.g., σ χA ∝ A 4 σ χN for m χ m A , that are widely assumed to be model independent.
We systematically examine the validity of the assumptions used to derive these relations, calculating where model independence ends. Figure 7 summarizes our results. We find the following: 1. For small cross sections, σ χN 10 −32 cm 2 , the usual scaling relations are valid, and multiple reasonable models can produce the same scaling relation.
2. For 10 −32 cm 2 σ χN 10 −25 cm 2 , the assumed A 4 scaling for a contact interaction progressively fails for all nuclear targets as cross sections for heavier nuclei begin to saturate at their geometric cross sections. Experimental constraints on the existence of light mediators prevent simple light mediator models from achieving cross sections in this range at all, such that constraints set in this range of cross sections should be specialized to a model. 3. For σ χN > 10 −25 cm 2 , dark matter cannot be pointlike. Contact interactions cannot achieve cross sections larger than the geometric cross section σ χA 4πr 2 A , and simple light mediators are strongly ruled out. Dark matter with cross sections in this range must be composite.
The failure of the scaling relations should influence the design of future dark matter searches. For interactions with cross sections that scale less than linearly with A, such as some models of composite dark matter, dark matter detectors with lighter nuclei are more efficient per unit detector mass. As a result, future direct-detection searches for strongly interacting dark matter may benefit from constructing detectors with light nuclei.
Constraints on dark matter parameter space are most useful if they can be compared between different experiments. Where the A 4 scaling is not reliable, results need to be recast in terms of specific models. A comprehensive analysis should include clear statements about the mass ranges their results can reasonably be extrapolated to. Because constraints will not be the same for different models, plots including cross sections σ χN 10 −32 cm 2 must specify a model, whether it involves a contact interaction, light mediator, composite dark matter, or something else.

ERRATUM
This paper focused on calculating where model independence ends for the scaling relations relating the spinindependent elastic-scattering cross sections of dark matter (m χ 1 GeV) between different nuclei. We also made some broader statements about regions of parameter space that are ruled out on theoretical grounds for pointlike dark matter. In the main text, these statements are accompanied by appropriate caveats; however, in some places -mostly the the abstract, introduction, and conclusions -the caveats were not always present in isolated passages. In this erratum, we note some of those missing caveats more explicitly. None of our main results or conclusions are affected.
An example of a type of statement that needs clarification is given in the Conclusions: "For σ χN > 10 −25 cm 2 , dark matter cannot be pointlike. Contact interactions cannot achieve cross sections larger than the geometric cross section σ χA 4πr 2 A , and simple light mediators are strongly ruled out. Dark matter with cross sections in this range must be composite." For repulsive interactions, this statement is correct. However, there are exceptions for attractive interactions, as shown in, e.g., Fig. 3, which shows resonances that can exceed the geometric limit. Though resonances are a generic feature of strongly attractive potentials (see, e.g., Refs. [96,126,127]), the details of resonances are inherently model dependent, as discussed in Sec. III C.
When the incident momentum is small, as for very light nuclei or DM with 1 GeV m χ m A , and the potential is attractive, we acknowledge that the cross section can exceed the geometric limit by 3-5 orders of magnitude at typical halo velocities. We show this for helium in Fig. 3, which assumes m χ m A . But when m χ m A , no nuclei heavier than neon can have attractive cross sections even one order of magnitude larger than the geometric limit. The requirement of m χ m A is appropriate for nearly all of the mass range considered (e.g., the mass range in Fig. 7 spans eighteen orders of magnitude), but we should have been more clear than there are some subtleties in the narrow but important range where m χ m A .
For clarity, we also reiterate a point that was mentioned in the "Overview of basic assumptions" in Sec. II A but perhaps not made clear enough; that unless otherwise specified, we focus on Milky Way virial velocities, v ∼ 10 −3 c, relevant for direct detection. As we discuss in Sec. VI C, attractive s-wave cross sections could be even larger for velocities lower than v ∼ 10 −3 c, as is relevant for dwarf galaxies and especially for early universe cosmology. However, the limits on light mediator masses and couplings discussed in Sec. IV B do not depend on the interaction velocity and are so restrictive that, in the unconstrained regions, the upper limits on the cross section would essentially be the same for both attractive and repulsive interactions. For a model finely tuned to evade some of the light mediator constraints in Fig. 5, or an interaction potential generated by a different mechanism entirely, attractive interactions at lower relative velocities could potentially achieve much larger cross sections. Such cases would be model-dependent by construction.
We encourage further work on the many interesting possible model-dependent interactions affecting dark matter scattering with nuclei. We also reiterate that the scope of the paper is explicitly restricted to m χ 1 GeV, so that we strongly discourage any attempt to extrapolate any of our results to lower masses without a careful supporting analysis. Such analyses will be an important complement to our work on where model independence ends. We thank Eric Braaten and Glennys Farrar for helpful discussions on these issues.
at r = R. We can obtain the wave function for r < R by directly integrating the one-dimensional Schrödinger equation, where we have defined u l (r) ≡ rA l (r). Note that because we are matching the phase shift with the logarithmic derivative of the wave function, the overall normalization of A l (r) is irrelevant for our purposes and can be chosen arbitrarily. We can then obtain A l (R) for any arbitrary potential V (r) by analytically or numerically evaluating Eq. (C6) from r = 0 to r = R.
Once we have A l (R), we can obtain δ l using tan(δ l ) = kRA l j l (kR) − (dA l /d ln r)| r=R j l (kR) kRA l n l (kR) − (dA l /d ln r)| r=R n l (kR) .
We have avoided canceling A l (r) in order to preserve signs, to ensure we obtain the correct quadrant for δ l (k).
The boundary condition at r = 0 should properly be u l (0) = 0, but for numerical solutions taking the boundary to be at r = 0 is inconvenient because of the 1/r 2 centrifugal term. Instead we can take advantage of the arbitrary normalization of A l (r) and fix the boundary conditions u l (r min ) = 1, du l /dr(r min ) = (l + 1)/r min , where r min is some small minimum radius.
As k → 0, it can be shown [93] that generically δ l (k) ∝ k 2l+1 , except in special cases where it is possible to achieve δ l (k) ∝ k 2l−1 for a specific value of l. Inspecting Eq. (C2), we can see the contributions from l = 0 and l = 1 are the only values of l which can be nonvanishing as k → 0. The l = 0 cross section is called the s-wave cross section.