# The Pirate Solution

## WIMP Dark Matter Direct Detection

All direct WIMP dark matter detectors aim to observe nuclear recoils produced by WIMP-nucleon scattering. WIMPs interact with all target materials in exactly the same way as neutrons and generally scatter off nuclei rather than electrons due to their relatively large mass. There are a few different methods that can be used to detect nuclear recoils, including collecting ionization, scintillation or thermal energy deposition data. The energy spectra of nuclear recoils for different target materials and WIMP masses can be simulated and calculated, then used to compare with real data from detectors to aid in the identification, or in placing an upper limit on characteristics (usually mass and WIMP-nucleon scattering cross-section) of WIMPs.

### Rates for Direct Detection of Dark Matter

#### Direct Detection Formalism

In this exercise we want to compute the cross section relevant for “direct” detection of Dark Matter WIMPs which is based on the elastic scattering of WIMPs off nuclei in a detector. In the case of spin-independent interactions, the procedure is fairly straightforward. The nucleus can be approximated as a collection of nucleons with overlapping wavefunctions, so an effective WIMP-nucleus interaction may be written:

are the coupling strengths of the dark matter particle to neutrons and protons respectively. are calculated from a coherent sum over the couplings to the quark model constituents of the nucleon. is the mass number of the nucleus. Integrating this interaction-Lagrangian density over space to get the matrix element for a transition from a nucleus at rest to a nucleus with momentum gives

is the probability of a given nucleon being at averaged over all nucleons. is the Fourier transform of which is called the “nuclear form factor” because it expresses the dependence of the interaction on the shape of the nucleus, given by is the proton number of the nucleus. Under the assumption that electric charge is distributed in the same way as the nucleons (i.e., there is no segregation of neutrons and protons), low-energy lepton-nucleus elastic scattering can be used to determine and It is observed that low-energy scattering is isotropic, so and Calculation of the differential cross section from the above matrix element for elastic → scattering and Fermi's Golden Rule gives

The factor arises from final-state density of states and the standard Golden Rule factor. This may be written in the form

where is the reduced mass of the nucleus and WIMP and is the relative WIMP-nucleus velocity. gives the "pointlike" total cross section if there were no form-factor suppression which corrects for the above formula being strictly correct only as the momentum transfer →

Direct DM search collaborations quote constraints on WIMP-nuclei cross sections normalized to the WIMP-nucleon cross section. The total WIMP-nucleon cross section is given by

The tiny mass difference between a proton and a neutron has been neglected. The above definitions yield : The WIMP-nucleus cross section is usually written differentially in the recoil energy , i.e., is the kinetic energy of the struck nucleus in the lab frame:

if (which is true in most WIMP models, up to small corrections, e.g., for the lightest supersymmetric neutralino, and for all WIMPs which interact primarily through Higgs exchange) and are assumed. The differential WIMP-nucleus scattering cross section depends on astrophysical, particle physics, and nuclear physics inputs.
We calculated the differential cross section for the scattering of a WIMP with given energy on a nucleus. Here we want to calculate the differential scattering *rate*, for given cross section and given WIMP distribution function. Recall that the scattering rate is the product of cross section and incident flux.
In general, the rate of interactions per unit mass of target material of a particle with interaction cross section (particle physics input) is given by

where is the number of target nuclei per unit mass (detector physics input), is the dark matter particle velocity in the lab frame (astrophysics input), and is the differential particle density. Including the normalization factor gives

with

and

- Here is the mean dark matter particle number density ( for dark matter particle mass [particle physics input], local dark matter halo density [astrophysics input]), is velocity onto the (Earth-borne) target, is Earth (target) velocity relative to the dark matter distribution, and is the local Galactic escape velocity (astrophysics inputs); is the particle density of dark matter particles with relative velocities within velocity space volume element about Then, the total event rate is

where is the average dark matter particle velocity in the lab frame. The differential rate per unit detector mass in nuclear recoil energy is given by

The quantity is the normalization with the dependence on shown explicitly. The WIMP velocity distribution is given by the WIMP phase-space distribution. The full phase-space distribution function for an isothermal halo composed of particles of mass is

which is a Boltzmann distribution with the gravitational potential included. At a given point, the position-dependent factor is fixed and can be included in the normalization, leaving a simple Maxwellian velocity distribution: : From the observational side, the most relevant piece of information coming from astrophysics is related to the rotational velocity of objects bounded to the Galaxy: The quantity is characteristic of the WIMP kinetic energy It corresponds to the most probable velocity and is known as the velocity dispersion. For The velocity distribution entering the differential rate is evaluated in the frame of the detector, therefore we need to transform the galactic velocity distribution to the Earth frame via a Galilean transformation. "Galactic" coordinates are usually used. and are simply related through: : where (). and denote the velocities of the Earth and the Sun in the Galactic rest frame and is the Earth’s orbital velocity around the Sun. is the sum of the Sun's circular velocity relative to the galactic center, and peculiar velocity relative to Sun's circular velocity. Projecting in the galactic plane, one gets the Earth’s speed relative to the galactic halo: : where is the inclination of the plane of rotation with respect to the galactic one and corresponds to the day when the Earth’s velocity is at its maximum. and . We set the most probable dark matter speed in the galactic frame to . The Maxwellian distribution is cut off at by the halo escape velocity. Note that the cutoff is isotropic in the galactocentric WIMP velocity not in the Earth-centric Heuristically, we may write : The total number of nuclear recoil events in a recoil energy range between and is

where the sum is over each nuclear species in the detector. is the effective exposure of species expressed in kg-days. We can parametrize where is the time of exposure, is the target mass of species and is a detection efficiency.

*Sheldon was calculating the rate for the wrong target material, xenon not sodium. Some experiments competing to find dark matter by observing directly the extremely small amount of energy a dark matter particle might deposit in a detector as they pass through Earth use Sodium (which has an atomic mass 23) and others use Xenon (with atomic mass 131). This is why Raj erases 131 and changes it to 23.*

In The Pirate Solution, British astronomer Professor Laughlin's first name is never stated. However, the plaque on his door reads: PROF. ___ LAUGHLIN, PhD. (first line) ASTRONOMY (second line). From this Blu Ray HD image, it appears that the obscured name on the plaque is IAN. If this is the case, then he has the same name as a Marvel Comics colorist.