Foraging for dark matter in large volume liquid scintillator neutrino detectors with multiscatter events

We show that dark matter with a per-nucleon scattering cross section $\gtrsim 10^{-28}~{\rm cm^2}$ could be discovered by liquid scintillator neutrino detectors like BOREXINO, SNO+, and JUNO. Due to the large dark matter fluxes admitted, these detectors could find dark matter with masses up to $10^{21}$ GeV, surpassing the mass sensitivity of current direct detection experiments (such as XENON1T and PICO) by over two orders of magnitude. We derive the spin-independent and spin-dependent cross section sensitivity of these detectors using existing selection triggers, and propose an improved trigger program that enhances this sensitivity by two orders of magnitude. We interpret these sensitivities in terms of three dark matter scenarios: (1) effective contact operators for scattering, (2) QCD-charged dark matter, and (3) a recently proposed model of Planck-mass baryon-charged dark matter. We calculate the flux attenuation of dark matter at these detectors due to the earth overburden, taking into account the earth's density profile and elemental composition, and nuclear spins.


I. INTRODUCTION
Despite ongoing experimental inquiry, the characteristics and couplings of dark matter remain a compelling puzzle. In this document we show how dark matter may be discovered by looking for multiple-scatter nuclear recoil events at large volume liquid scintillator experiments, which are traditionally used to detect solar neutrinos with MeV energies.
Ongoing dark matter direct detection experiments typically seek dark matter that scatters in underground detectors at most once. These are often called wimps, or weakly interacting massive particles, and their signature is a recoiling nucleus with ∼ keV kinetic energy. As some of us recently showed, traditional dark matter direct detection experiments are also capable of seeking dark matter candidates that scatter multiple times as they transit these detectors [1]. Because they have a large cross section with Standard Model particles, these mimps, or multiply interacting massive particles, leave a signature different from wimps. Mimps recoil against nuclei multiple times as they travel through dark matter direct detection experiments, leaving a characteristic collinear trail of deposited nuclear recoil energy. Mimp searches extend the reach of dark matter direct detection searches to the largest dark matter masses; the maximum dark matter mass an experiment can search for is determined by the total dark matter flux admitted by the detector, which in turn is determined by the area of the detector and the observation time. Current xenon-based and bubble chamber experiments may eventually probe dark matter up to Planck masses (M Pl ).
On the other hand, as this document will show, dark matter with mass greater than M Pl can be readily sought in a number of existing and imminent large volume neutrino detectors, such as borexino, sno+, and juno. While the energy threshold of neutrino detectors is typically larger than dedicated dark matter experiments, their larger volumes collect over three orders of magnitude more dark matter flux. Of course, neutrino detectors are unsuitable for a standard wimp search, where a single O(10 keV) recoil energy deposition falls well short of O(MeV) thresholds that trigger these experiments. However, dark matter with sufficiently many scattering interactions will impart enough aggregate kinetic energy to trigger these detectors, rendering them well suited for mimp searches. Their cross section vs mass reach is given by the region above and to the left of the blue lines in Fig. 1. We will calculate this reach in the following.
Other proposed searches for dark matter that deposits enough energy to exceed the ∼ MeV threshold of neutrino detectors include dark matter that destroys target baryons [2][3][4][5], that yields annihilation or decay products detectable in these experiments [6][7][8][9][10][11][12][13][14][15][16], that deposits its entire (mass + kinetic) energy [17], is produced at high-intensity accelerators or radioactive sources [18], or is bounced off energetic cosmic rays [19,20]. All these approaches (excepting the last one) require specific models with a number of massive dark sector states. In contrast, our proposed search is model-independent and sensitive to any single dark matter state with a large nuclear scattering cross section over a range of masses. Another important difference is that in our study,Čerenkov light-based detectors like sno and Super-Kamiokande will not be effective, since mimp scattering does not excite electrons/nuclei to relativistic speeds. Future argon-based detectors such as Protodune and dune, if operated in single-phase mode with particle energies measured via electronic readout, may also be unsuitable to our study, since the energy threshold is around 5 MeV for even a single scatter. Lowering this threshold by operating in dual-phase may make these detectors more interesting for mimp searches.
The remainder of this document proceeds as follows. In Sec. II, we detail the dark  Figure 1. Sensitivities, as derived in Sec. II, of liquid scintillator detectors to multiscatter-inducing dark matter undergoing spin-independent (top) or spin-dependent (bottom) scattering, illustrated with borexino. The mass reach is determined by the total flux of dark matter transiting the detector, while the detection threshold cross section is set by the minimum rate of energy deposition required for the photomultiplier tubes to tell significantly between signal scintillation and dark counts. Due to similar dimensions and detection technology, sno+ would have very similar sensitivities; due to larger dimensions, juno would reach greater dark matter masses. The vertical green line is the maximum mass reach of xenont after 1 year [1]. The black lines indicate upper bounds (derived in Sec. IV) on the cross section from the slowing down of dark matter by the earth's atmosphere and rock. The dot-dashed (dashed) lines indicate the halving (tenthing) of total dark matter flux by this earth overburden. Also shown are spin-independent scattering bounds from molecular gas cloud cooling [21], Skylab [22][23][24], and searches for dark matter tracks in ancient mica [25,26]. Similar bounds should apply to spin-dependent scattering; this would require a re-analysis of these constraints accounting for the fraction of spin-containing nuclei in gas clouds, Skylab detectors, and ancient mica. The horizontal brown line denotes the maximum cross section imposed by unitarity at zero angular momentum; this cross section could be higher for higher modes.
matter-nucleon scattering cross sections that could be uncovered using large volume liquid scintillator detectors. We consider both a conservative search, which requires dark matter passing the triggers defined at these experiments, and an aggressive future search, where the irreducible background at large volume liquid scintillator experiments is estimated using the dark count rate in photomultiplier tubes (pmts) as reported by borexino. In Sec. III we discuss three scenarios of super-heavy dark matter that scatters multiple times in and could be discovered by large volume liquid scintillator experiments. In Sec. IV we provide overburden calculations, determining the fraction of multiply interacting dark matter that is blocked from detection by scattering with the earth. In Sec. V we conclude.

II. LIQUID SCINTILLATOR SEARCHES FOR DARK MATTER
Large volume liquid scintillator experiments like borexino, sno+, and juno can be used to detect strongly interacting, heavy dark matter. The dark matter scattering multiplicity must be large enough that enough energy is deposited to exceed the trigger threshold of these experiments, which are primarily designed to detect neutrinos with MeV energies. The dark matter must be heavy enough that its initial kinetic energy greatly exceeds the energy lost scattering with the intervening atmosphere and earth as it travels to these experiments situated a few kilometers underground. The advantage of large liquid scintillator experiments is that they are a few orders of magnitude larger than traditional direct detection experiments and therefore intercept an (a few) 2 orders of magnitude larger flux of dark matter.
To explore sensitivities to dark matter at these experiments, we begin with the canonical, momentum-independent spin-independent and spin-dependent dark matter-nucleon scattering cross sections. In Sec. III we discuss a number of explicit dark matter models which permit large masses and sizable nuclear scattering cross sections. For a target with atomic number A, dark matter-nucleus reduced mass µ Tχ , and dark matter-nucleon reduced mass µ nχ , the spin-independent si and spin-dependent sd per-nucleon cross sections in terms of the per-nucleus cross section are given by where the proportionality is in the m χ m T limit. The effect of nuclear form factors is negligible in our treatment, hence they have been ignored here. We will consider isospinindependent scattering and set a n = a p = 1, although it would be possible to extend these results to the case a n = a p . The above formula only applies to scenarios where perturbativity is valid; for non-perturbative scenarios, such as when dark matter is comprised of large bound states, only per-nuclear cross sections are relevant.
Because dark matter will only have to scatter (multiple times) to be detected at sno+, borexino, and juno, a determination of which dark matter parameters can be probed by these detectors requires computing just two quantities: (i) the maximum dark matter mass to which the experiment is sensitive, and (ii) the minimum cross section for the dark matter to exceed the energy thresholds of the experiments. We will estimate these two quantities in the following, assuming uniform dark matter velocity v χ = 10 −3 c. This is a good approximation to a full estimation of these quantities assuming a Maxwell-Boltzmann speed distribution, which we will encounter in Sec. IV.

A. Mass reach
The maximum dark matter mass that large liquid scintillator detectors are sensitive to, m max χ , depends upon the detectable dark matter flux passing through the detector. Because the dark matter density around earth is constant, the largest dark matter mass probed by the detector will correspond to a flux that delivers a few detectable dark matter events over the observation time of the detector, , v χ and m χ are respectively the local density, average velocity, and mass of incident dark matter, A det is the effective area of the fiducial volume of the detector, and t exp is the observation time. From this, we obtain Here we have normalized A det to the area of the spherical, 8.5 meter diameter, "inner volume" region of borexino [28], which yields A det = πr 2 5 × 10 5 cm 2 . We have indicated the borexino 10-year sensitivity with a vertical blue line in Fig. 1. At sno+, A det 10 6 cm 2 , hence m max χ 2.8 × 10 21 GeV after 10 years. Similarly, at juno, A det 10 7 cm 2 , hence m max χ 2.8 × 10 22 GeV after 10 years. For a sufficiently large dark matter-nucleon scattering cross section, the mass reach of large volume scintillator experiments will diminish, since a portion of the dark matter flux will become undetectably slow after scattering many times with nuclei in the earth. This diminished mass reach is also shown in Fig. 1. We discuss this further in Sec. IV.

B. Cross section reach with existing triggers
Next we determine the minimum dark matter-nucleon cross section required to trigger large volume liquid scintillatior neutrino detectors. At some of these detectors, this threshold is quoted in terms of the energy deposited per ∼100 ns window. At borexino, the detector is triggered if E trig 50 keV of energy is deposited within a 100 ns timing window [28]. In this study, we consider a dark matter flux which travels 3 − 10 meters through these detectors over a time 10 − 30 µs. We note that assuming the dark matter travels at least 3 meters alters the flux calculations given above by less than 10%. The signal of such a mimp would be a continuous emission of photons across a succession of many O(100) ns timing windows. The typical lab frame recoil energy of a carbon nucleus struck by heavy dark matter = m T v 2 χ in the limit m χ m T , where m T is the mass of the target nucleus. The density of the scintillator pseudocumene (C 9 H 12 ) is 0.88 g/cm 3 . Therefore, using the mean scattering rate formula for scattering on carbon with number density n C (Γ s = n C σ Cχ v χ ), the 50 keV per 100 ns trigger threshold amounts to a per-carbon-nucleus cross section of σ Cχ = 3.7 × 10 −23 cm 2 . From Eq. (1), this corresponds to a per-nucleon scattering cross section σ (si,min, 100% eff.) nχ = 1.8×10 −27 cm 2 , set by scattering on 12 C. The threshold for a spindependent cross section is set by scattering on 13 C. Accounting for its natural abundance = 1.1%, and from Eq. (1), this threshold is given by σ  [29]. The effect of sd scattering on (the more abundant) 1 H is relatively negligible, since the energy deposition is small, and since Eq. (1) implies poor sensitivity to small per-nucleon cross sections for A=1. Note that this is because we had assumed that a p = a n in Eq. (1); were a n = 0, the sd scattering sensitivity would be set by 1 H recoils. However, these cross sections have assumed a 100% efficiency for conversion of carbon nuclear recoil energy into photons detectable by pmts. In reality, we expect only about 10% of the carbon recoil energy to be converted to detectable photons. In a study of liquid scintillation nuclear recoils [30], the photoelectron yield for carbon and hydrogen nuclear recoils were determined for recoil energies as small as 40 keV. Extrapolating their result using their modified Birks' formula, the light yield efficiency for 10 − 20 keV nuclear recoils is approximately 10%. This results in the scaling up of the minimum cross sections quoted above by a factor of ten.
Therefore, the borexino trigger outline above implies that heavy dark matter can be detected with a minimum per-nucleon spin-independent scattering cross section given by and a minimum spin-dependent scattering cross section given by where E trig is the 100 ns trigger energy of the detector. We have plotted these threshold cross sections in Fig. 1 with solid blue horizontal lines. We have also shown existing constraints on spin-independent scattering from molecular gas cloud cooling [21], as well as from searches for dark matter tracks in plastic etch detectors on board Skylab [22][23][24], and in ancient muscovite mica [23,25]; see Appendix A for details of our recasting of mica bounds. Also shown in Fig. 1 is the maximum cross section dictated by unitarity, 4π(2 + 1)/q 2 C , for the partial wave mode = 0, with q C the typical momentum transferred to a carbon nucleus. This cross section could of course be higher if scattering proceeds in higher partial wave modes. The black and gray lines are cross sections above which the earth overburden significantly degrades dark matter kinetic energies; we will treat this more carefully in Sec. IV.
Note that, the efficiencies obtained in [30] imply that on average one photoelectron will be emitted per ∼ 10 keV carbon-atom recoil. This means, with a 10%-efficient light yield, the borexino trigger is activated if a minimum of 50 photoelectrons are detected in a 100 ns window, or 5000 in a 10 µs window. In the next subsection we argue that a photoelectron count that is O(100) smaller can be tolerated in mimp searches, with the use of new trigger software.

C. Cross section reach with improved triggers
It should be possible to improve the cross section sensitivity of large volume liquid scintillator experiments with the addition of a dedicated trigger. As long as this trigger is sensitive to excessive photoelectrons detected within a certain time window, mimps may be sought. As discussed, the window of time for a dark matter particle to transit these detectors is ∼ 10 µs. The pmt dark count rate reported by borexino [28] is roughly 10 photoelectrons per 10 µs. Using this, one may now determine the minimum cross section -or minimum multiplicity of scatters -potentially accessible at a similar future experiment, such as sno+. Assuming 10 µs timing windows, we will demand that the aggressive trigger threshold admit no more than a few randomly generated dark count events over the course of the experiment. Said differently, we demand that the dark matter cross section sought with the aggressive trigger be just large enough that the maximum dark counts in any 10 µs window over the lifetime of the experiment, will only produce about one background event. Then, assuming Poisson statistics the required multiplicity N c is obtained by solving where N bg = 10 is the dark count in a 10 µs window as mentioned above, and t life is the run-time of the experiment. For t life = 10 years, we get N c = 42, which corresponds to a minimum detectable spin-independent cross section of and a minimum detectable spin-dependent cross section of We have plotted these improved threshold cross sections in Fig. 1 with dashed blue horizontal lines. 14 C background. The β decay of carbon-14, which produces ∼ 20 − 150 keV electrons, has the highest rate of background processes observed at underground liquid scintillator experiments (see e.g. Section XI.3 in [28]). Using liquid scintillator refined from petroleum extracted deep underground, an abundance of 14 C/ 12 C as low as 2 × 10 −18 has been obtained at borexino [31]. With this purity, a kiloton of liquid scintillator will see 14 C decays at a rate of about one kHz. Additionally, we note that pmt hits from β decays last <∼100 ns [28]. Therefore, we should expect a 14 C β decay event in around one out of every hundred 10-µs mimp search windows. However, we should not expect more than 4 14 C events in any 10 µs window over the lifetime of the detector, meaning that we expect at least 10 − (4 × .01) = 9.6 µs of every 10 µs window to be background-free (except for dark counts, which we have discussed above). Using this, and the fact that the observed energy deposition rate for 14 C β decays is a few orders of magnitude higher than the energy deposition for a mimp in our aggressive search region (σ (si min, trig+) nχ 10 −26 cm 2 ), it should be straightforward to design a suitable trigger for the aggressive mimp search region that rejects 10 µs windows with any 14 C β decays. In principle, it should also be possible to design the aggressive trigger to "scrub" 14 C β decays, by excluding 100 ns windows with pmt hits occurring in excess of a certain threshold. The particulars of the aggressive trigger will depend on the experimental capabilities of e.g. borexino , sno+, juno, and we leave a detailed study to future work.

D. Directionality of recoiling nuclei in MIMP searches
Multiply interacting dark matter, with a per-nucleon cross section above the threshold cross sections identified in Eqs. (3), (4), (6) and (7), will scatter at least 40 times, yielding a minimum of 40 detectable photoelectrons (see discussion in Sec. II B), as it transits a large volume liquid scintillator-filled neutrino detector. The sites of these nuclear recoil interactions will lie along an extremely straight line pointing through the detector, since the deflection of a heavy particle by a carbon nucleus is, in the small-angle scattering limit, This collinearity of mimp events can be used to validate candidate super-heavy dark matter events, by using the location and timing of pmt hits, and thereby reconstructing the track of the mimp as it transits the detector. While it is not feasible to resolve the directionality of the photons themselves, it has been demonstrated that by collecting the arrival time and pmt location of all photon hits, directional information can be extracted; see e.g. [32], which showed how particle tracks in liquid scintillator can be reconstructed for the purposes of identifying neutrinos from dark matter that annihilates in the sun. The precise details of such a directional analysis will depend on the data acquisition and timing capabilities of a given experiment.

III. SURVEY OF MULTIPLY INTERACTING DARK MATTER MODELS
There are a number of theories of dark matter that may be discovered at large volume liquid scintillator detectors. We will outline three scenarios here: an effective field theory description of dark matter, quantum chromodynamic (qcd)-charged dark matter, and heavy dark matter mediated by a light vector carrying baryon number. A number of works demonstrate the production of super-heavy dark matter in the early universe [33][34][35][36].

A. Effective operators
In Fig. 1 and Eqs. (3), (4), (6) and (7), the minimum cross sections required to trigger these detectors are a few orders of magnitude below those corresponding to qcd hadronic scattering. We should therefore expect that a mediator transmitting these dark matter interactions has a minimum mass around the qcd scale 100 MeV.
This intuition is borne out by the effective field theory cutoff for dark matter scattering interactions that just exceed the thresholds of large volume liquid scintillator experiments. Taking dark matter as a Majorana fermion χ, we may characterize with effective operators its parton-level interactions with nucleons. For illustration we pick two operators: the scalarscalar (ss) and axial-axial (aa) operators, describing respectively si and sd scattering in the non-relativistic limit: Following the usual convention (see [37,38]), in order to evade stringent constraints from flavor observables, we assume the principle of minimal flavor violation [39], whereby the U (3) 5 flavor symmetry of the Standard Model is broken only by Higgs Yukawa spurions. Consequently, the scalar bilinears for quarks are proportional to the corresponding Higgs Yukawa coupling (or quark mass), and the axial bilinears are flavor universal.
We can now recast the threshold cross sections derived in Sec. II, specifically assuming the new triggers proposed in Eqs. (6) and (7). The per-nucleon cross section for the above operators is given by [37,38] where the effective nucleonic couplings are with v = 246 GeV the Higgs expectation value, and the nucleon form factors at zero momentum are taken as [40,41] {f Using this, we find that large volume liquid scintillator detectors like borexino and sno+ are potentially sensitive to the cutoffs Λ ss ≤ 59.8 MeV , Λ aa ≤ 55.9 MeV .
Note that for the effective operator treatment to be valid, Λ 11 MeV, the typical momentum transfer for dark matter scattering scattering off carbon.

B. QCD-charged dark matter
The low-energy nucleon-on-nucleon elastic scattering cross section in the Standard Model (∼ 10 −26 cm 2 ) is similar to the dark matter cross section sensitivity thresholds found for large volume liquid scintillator experiments in Sec. II. Indeed, upon further investigation, dark matter charged under the Standard Model qcd gauge group SU (3) C provides a generic and viable class of models that can be sought out by these experiments.
Dark matter charged under SU (3) C has been considered by a number of authors, e.g. [42][43][44][45]. We will find that the super-massive ( 10 17 GeV) dark matter considered here will likely need to transform as an octet of SU (3) C , once constraints on electromagnetically (em) charged particles from the macro experiment [46] are taken into account. The logic is as follows. Consider two color-charged dark matter states, a stable electroweak-singlet quark-like X 3 , and a stable gluon-like X 8 . Formally, these states transform as a triplet (X 3 ) and an octet (X 8 ) of SU (3) C . After the qcd phase transition, X 3 will form color-neutral states such as X 3q and X 3 qq, where q are Standard Model quarks. Notice that while these X 3 composite states are color-neutral, they carry O(1) em charge. The macro experiment, by not observing minimum-ionizing particles in liquid scintillator tanks in the Gran Sasso mine, set a flux limit on heavy em charged relic particles (see Fig. 1 of [46]): Assuming that this search was sensitive to O(1) em-charged particles (see [47]), the macro search excludes such particles with masses 10 8 GeV m X 3 10 20 GeV. This mass bound is an order of magnitude greater than quoted in [47], possibly due to a transcription error in [47] that quotes the flux limit an order of magnitude higher than Eq. (12) when deriving a bound from macro. While we expect the lightest X 3 -containing states (X 3q and X 3 qq) to be em-charged, we expect the lightest X 8 -containing states to be em-neutral. These states are X 8 qq and X 8 g, where q is an up or down quark and g is the Standard Model gluon. It is reasonable to assume that the lightest X 8 -containing states are em-neutral, as is the case for the lightest pion in qcd, and the X 8 g state, which could be even lighter, is manifestly em-neutral. Detailed lattice computations will be necessary to establish the exact mass spectrum for X 8 qq and X 8 g states. The elastic scattering cross section for X 8 qq and X 8 g, which will be surrounded by a cloud of quarks and gluons after confinement, should fall in the range where this is the typical range for nucleon-nucleon cross sections at low energy. We defer to future work a complete cosmological treatment of superheavy qcd-charged dark matter, as well as a lattice calculation to determine the structure and mass spectrum of X 8 qq and X 8 g bound states.
C. Baryo-charged dark matter with a sub-GeV mass mediator The effective operator treatment in Sec. III A, which found operator cutoffs of ∼ 100 MeV for the cross section thresholds quoted in Sec. II, indicates that dark matter scattering with GeV mediators may be detected at large volume liquid scintillator experiments. Recently, Ref. [48] discussed models of Planck-mass mimps mediated by a light vector boson that carries gauged baryon number U (1) B , under which dark matter is charged. While such a light baryon-charged mediator is subject to a number of phenomenological constraints [48][49][50][51][52][53], they nevertheless allow for mimps that could be detected by large volume liquid scintillator experiments. Assuming dark matter carries baryon charge Q χ , for a U (1) B gauge boson with coupling strength α d and mass m A d , and normalizing with parameters unconstrained by meson decays, the si per-nucleon scattering cross section is showing that borexino and sno+ are capable of probing this scenario.

IV. EARTH OVERBURDEN FLUX ATTENUATION
If dark matter's cross section for scattering on Standard Model particles is sufficiently large, then dark matter particles will be slowed by repeated scattering while passing through Left. Geometric schematic of a dark matter particle traveling along a straight path to a dark matter detector through the Earth's atmosphere and rock. This diagram is not to scale. Right. Contours in the cos θ-m χ plane, of the spin-independent per-nucleon cross sections (in cm 2 ) required for the Earth overburden to slow down dark matter below detector thresholds. This illustrative calculation assumes an isotropic dark matter distribution with a uniform speed = 220 km/s. the earth's atmosphere and crust. For non-relativistic scattering, the fraction of kinetic energy depleted per scatter is zβ ≡ 4z(m T m χ )/(m T + m χ ) 2 , where z ∈ [0, 1] is a kinematic factor that depends on the scattering angle (and m χ , m T are respectively the dark matter and target masses). Now consider a dark matter particle traversing a straight line path of length L through the earth atmosphere and rock on its way to a detector, where r det is the distance between the detector and the center of the earth, R e is the radius of the earth (including its atmosphere), and θ is the angle between the incoming particle velocity and the line from the center of the earth to the detector; see the left-hand side of Figure 2. At some distance s along the path of length L that the dark matter is traversing, the distance to the center of the earth r is given by After traversing a path L, the final kinetic energy of the dark matter is given by where n A (r) is the number density of nuclei of type A at distance r from the center of the Earth and we take an average value of z = 1 2 in computations. The Earth and its atmosphere are assumed to be spherically symmetric, and our model of the density of Earth's interior follows the Preliminary Reference Earth Model [54]. See also [1] for the dimensions, densities and composition of the core, mantle and crust. The S p , S n values of spincarrying elements are taken from [29].
A good picture of the effects of the overburden may be gained by assuming an isotropic dark matter distribution with a uniform speed of 220 km/s. The length of dark matter's path through the Earth to the detector, related to cos θ by Eq. (16), determines the cross section σ nχ above which E f is too low to trigger the detector. For concreteness, we set the detector energy threshold to 1 keV. In the right hand side of Fig. 2 we show, in the cos θ-m χ plane, contours of σ nχ /cm 2 required for such a slow-down via spin-independent scatters. The visible change of slope near cos θ = 0.1 and −0.8 is due to the transit of dark matter through the mantle and core, which are (progressively) denser than the crust. The overburden cross sections at cos θ = 1 and cos θ = −1 are seen to agree with those derived in [1]. For overburden cross sections between cos θ = 1 and cos θ = −1, reconstructing directionality of dark matter in the detector (as discussed in Sec. II D) would provide a handle on the dark matter mass and cross section; see [1] for more details. Above cos θ = 0.2 the slope of the contours is steep, indicating that the detectable dark matter flux becomes exponentially sensitive to the crustal overburden cross section.
A more detailed calculation of the overburden bounds must take into account the dark matter velocity distribution and statistical considerations. To this end, we generated a Monte Carlo sample of dark matter particles assuming an isotropic 1 population that follows the Maxwell-Boltzmann speed distribution (as seen by an observer moving through the halo) [55]: where v = 220 km/s is the speed of the sun through the halo, the velocity dispersion v disp = 220 km/s, and the galactic escape velocity v esc = 533 km/s [56], and c min = max[−1, (v 2 − (v 2 esc −v 2 ))/2vv ]; Θ denotes the Heaviside theta function. Bounds were obtained by locating (m χ , σ nχ ) pairs that yield 2.3 expected mimps over 10 years, with enough kinetic energy to deposit an average of 11 keV (the typical energy transfer to a carbon nucleus) in each carbon nuclear recoil as they transit borexino. These are presented in Fig. 1 as solid black lines. Due to the exponential sensitivity of the detectable flux to the overburden cross section (Eq. (17)), the relation between this cross section and dark matter mass is practically linear, and can be approximated by This approximate linearity in Fig. 1 One sees that σ SI,Φ/2 nχ corresponds to the cross section contours near cos θ = 0 in the right hand side of Fig. 2. In presenting these overburden bounds, we have neglected the effect of "saturated overburden scattering" [1], i.e. the insufficient slowing down of dark matter due to the finiteness of the number of terrestrial nuclei on its path to the detector. This effect applies only to short-range interactions, whereas our treatment here is more generic; when this effect does apply, cross sections above the overburden lines in Fig. 1 may also be probed by liquid scintillator neutrino detectors.

V. CONCLUSIONS AND DISCUSSION
It is sometimes the case that well-designed experiments are well-suited for novel physics searches beyond their intended scope. This document has shown that large volume liquid scintillator neutrino detectors stand at the vanguard of super-heavy dark matter detection. Because they are the largest terrestrial detectors capable of identifying a handful of nuclear recoil interactions over ambient backgrounds, experiments like borexino, sno+, and juno provide unmatched high mass dark matter sensitivity. After estimating the sensitivity of these experiments to high mass dark matter in Sec. II, we proceeded in Sec. III to explore how canonical and newly proposed models of dark matter would first be discovered by large volume liquid scintillator experiments.
We have also identified a major trigger improvement for future large volume liquid scintillator experiments, enabling the detection of dark matter-nucleon scattering cross sections over two orders of magnitude smaller than would be accessible with current triggers. We have shown that by recording experimental data whenever there are more than 40 pmt hits distributed over a 10 microsecond window, where this number of pmt hits should never arise from randomly distributed pmt misfires ("dark counts"), future liquid scintillator experiments can greatly increase their sensitivity to heavy dark matter.
It may be possible to further extend the sensitivity of large volume liquid scintillator experiments to even smaller nucleon-scattering cross sections, by incorporating directional information in these detectors. In Sec. II D, we demonstrated that multiply interacting dark matter traversing large volume liquid scintillator detectors will leave a characteristic, highly collinear trail of nuclear recoil sites, which will produce a pattern of scintillated photoelectrons depending on the timing and pmt positions. Prior work on reconstructing the directionality of neutrinos from dark matter annihilating in the sun, using pmt-instrumented liquid scintillator, indicates that this may be feasible at experiments like borexino, sno+, and juno. If directional capabilities are established for these detectors, daily and annual modulation of multiply scattering dark matter events can also be used to validate putative signals.