---
title: LUX Formalism in Dark Matter Experiments
url: https://www.emergentmind.com/topics/lux-formalism
type: topic
---

# LUX Formalism in Dark Matter Experiments

Searching arXiv for the specified LUX papers and closely related formalism papers.
The LUX formalism is the integrated analysis framework used by the Large Underground Xenon experiment to map dark-matter interaction models into observable distributions in a dual-phase liquid-xenon time projection chamber, calibrate those observables with in situ electron-recoil and nuclear-recoil data, and derive exclusion limits through profile-likelihood inference under explicit astrophysical and detector-response assumptions. In the LUX literature, the formalism appears in several closely related forms: a canonical elastic WIMP-search pipeline for spin-independent and spin-dependent scattering [1310.8214; 1512.03506; 1710.03572], a nuclear-recoil calibration and microphysics framework anchored by D–D neutron scattering [1608.05381], and an extension to sub-GeV dark matter using Bremsstrahlung and Migdal electron-recoil signatures when the nuclear recoil itself is below threshold [1811.11241].

## 1. Detector basis and observable space

LUX is a dual-phase xenon TPC in which an interaction in liquid xenon produces prompt scintillation, denoted S1, and ionization electrons that are drifted upward and extracted into the gas to generate proportional scintillation, denoted S2 [1710.03572]. The detector records both signals with 122 PMTs arranged in top and bottom arrays, and reconstructs depth from the S1–S2 drift time while the transverse position is inferred from the S2 light pattern on the top array [1710.03572]. In the first-results configuration, the active liquid xenon target was 250 kg within a TPC of 47 cm diameter and 48 cm drift height, with photodetection via two arrays of 61 PMTs [1310.8214].

A central feature of the LUX formalism is that analysis is performed directly in detector observables rather than in recoil energy alone. In the early 2013 analysis, the unbinned extended profile-likelihood ratio used four observables, radius, depth, S1, and S2\(_b\), where S2\(_b\) denotes S2 measured from the bottom PMT array alone [1310.8214]. In the 2013 reanalysis, the signal and background PDFs were constructed in \((cS1, cS2, r, z)\), where weekly \(^{83\mathrm{m}}\)Kr calibrations were used to define corrected observables that equalize detector response throughout the active volume [1512.03506]. In the combined WS2013 and WS2014–16 analysis, the observable spaces were \(\{r_{\mathrm{ver}}, z_{\mathrm{ver}}, S1, S2\}\) for WS2013 and \(\{r_{S2}, \phi_{S2}, \tau_d, S1, S2\}\) for WS2014–16, reflecting the more complicated field geometry of the later run [1710.03572].

The observables are expressed in units of detected photons (phd), not photoelectrons, to account for double photoelectron emission [1710.03572]. LUX also used photon counting for low S1 to improve resolution [1710.03572]. This detector-level representation is foundational: it is the space in which both calibration data and dark-matter signal models are compared to the search sample.

## 2. Response model, gains, and energy reconstruction

The LUX formalism connects deposited energy to measured signals through a quanta model. For electron recoils, the reconstructed energy is obtained by combining scintillation and ionization according to  
\[
E = W\left(n_{\gamma}+n_e\right)=W\left(\frac{S1}{g_1}+\frac{S2}{g_2}\right),
\]
with \(W = 13.7 \pm 0.2\) eV per quanta [1811.11241]. For the WS2013 sub-GeV analysis, the gains were \(g_1 = 0.117\) phd/photon and \(g_2 = 12.2\) phd/electron, with electron extraction efficiency \(49\% \pm 3\%\) [1811.11241]. The combined WIMP-search formalism reports closely related WS2013 values \(g_1 = 0.117 \pm 0.003\) phd/photon and \(g_2 = 12.1 \pm 0.8\) phd/electron [1710.03572], while the 2013 reanalysis quotes \(g_1 = 0.117 \pm 0.003\) phd per scintillation photon and \(g_2 = 12.1 \pm 0.8\) phd per extracted electron, with anticorrelation \(\rho = -0.6\) [1512.03506].

For nuclear recoils, the response model is expressed through the light yield \(L_y(E_{nr},F)\) and charge yield \(Q_y(E_{nr},F)\), which determine the number of photons and electrons:
\[
N_{\mathrm{ph}} = L_y(E_{nr},F)\,E_{nr},\qquad N_e = Q_y(E_{nr},F)\,E_{nr},
\]
and then
\[
S1 = g_1\,N_{\mathrm{ph}},\qquad S2 = g_2\,N_e
\]
[1608.05381]. The D–D calibration paper gives \(g_1 = 0.115 \pm 0.004\) phd/photon and \(g_2 = 11.5 \pm 0.9\) phd/electron for the calibration at \(F \approx 180\) V/cm, with extraction efficiency \(0.48 \pm 0.04\) [1608.05381].

The response model in later analyses was implemented with NEST, tuned to calibration data. In WS2014–16, variable electric fields necessitated splitting the exposure into 16 subdatasets, each with a unique NEST response model constrained by tritiated methane and D–D neutron calibrations [1710.03572]. The proceedings summary states that LUX used NEST to generate NR and ER Monte Carlo signal and background models for each of the 16 time/z bins, with yield formulae updated using LUX’s own tritium and D–D calibrations [1611.05525].

A plausible implication is that “LUX formalism” denotes not one equation set but a calibrated forward model: quanta production, transport, gain, correction, and fluctuation modeling are all embedded before any likelihood evaluation.

## 3. Nuclear-recoil microphysics and in situ calibration

The D–D neutron calibration established the nuclear-recoil branch of the LUX formalism by reconstructing recoil energy from double-scatter neutron kinematics [1608.05381]. The analysis defines the recoil energy using the center-of-mass scattering angle and employs the approximation
\[
E_{nr} \simeq \frac{2\,m_n\,M}{(m_n+M)^2}\,E_n\,(1-\cos\theta_{\mathrm{lab}})
\]
with better than \(2\%\) accuracy for all angles [1608.05381]. This kinematic reconstruction enabled absolute in situ measurements of \(Q_y\) from \(0.70\) to \(24.2\) keV\(_{nr}\) plus an endpoint at \(74\) keV\(_{nr}\), and \(L_y\) from \(1.08\) to \(12.8\) keV\(_{nr}\) plus an endpoint at \(74\) keV\(_{nr}\) [1608.05381].

The calibration demonstrated measurable signals down to \(\mathcal{O}(1\ \mathrm{keV}_{nr})\), with representative values near \(1.1\) keV\(_{nr}\) of \(Q_y \approx 7.4\ e^{-}/\mathrm{keV}_{nr}\) and \(L_y \approx 4.9\ \mathrm{ph}/\mathrm{keV}_{nr}\) [1608.05381]. At \(10.9\) keV\(_{nr}\), the reported representative values are \(Q_y \approx 5.9\ e^{-}/\mathrm{keV}_{nr}\) and \(L_y \approx 8.1\ \mathrm{ph}/\mathrm{keV}_{nr}\) [1608.05381]. The paper states that these results extend measurable NR response down to \(\sim 1\) keV\(_{nr}\) and improved low-mass WIMP sensitivity by a factor of \(7\) at \(7\ \mathrm{GeV}/c^2\), reducing the lowest accessible WIMP mass from \(5.2\) to \(3.3\ \mathrm{GeV}/c^2\) [1608.05381].

The microphysical energy scale uses the \(W\)-value and an electronic energy fraction \(L(E)\):
\[
E_{nr} = \frac{W\,(N_e+N_{\mathrm{ph}})}{L(E_{nr})}
\]
[1608.05381]. LUX fit both a Lindhard-based model and a Bezrukov/Ziegler-based parameterization to the data, with a Lindhard best-fit parameter \(k = 0.1735 \pm 0.0060\) and a biexcitonic quenching parameter \(\eta = 13.2 \pm 2.3\) [1608.05381]. In the 2013 reanalysis, the mean fraction of energy going to quanta for NRs was parameterized by the Lindhard model with \(k = 0.174 \pm 0.006\), treated as a Gaussian-constrained nuisance parameter in the likelihood [1512.03506].

These calibrations replaced earlier conservative assumptions. The 2013 reanalysis explicitly states that the previous analysis modeled the signal only above \(3\) keV minimum energy, whereas the revised analysis truncated the WIMP signal only below the lowest D–D calibration point of \(1.1\) keV [1512.03506]. This change dominates the low-mass improvement [1512.03506].

## 4. Astrophysical and scattering-rate framework

For elastic scattering, LUX adopts the standard differential-rate formalism under the Standard Halo Model. The 2013 first-results paper writes
\[
\frac{dR}{dE_R} = \frac{\rho_0}{m_\chi}\int_{v>v_{\min}(E_R)} \frac{f(\mathbf{v})}{v}\,\frac{d\sigma}{dE_R}(v,E_R)\,d^3v
\]
with
\[
\frac{d\sigma}{dE_R} = \frac{m_A\,\sigma_A^{\mathrm{SI}}}{2\,\mu_A^2\,v^2}F^2(E_R),
\]
where \(F(E_R)\) is the Helm form factor, and
\[
\sigma_A^{\mathrm{SI}} = \sigma_n\,A^2\left(\frac{\mu_A}{\mu_n}\right)^2
\]
for isospin-invariant spin-independent coupling [1310.8214]. The minimum speed is
\[
v_{\min}(E_R)=\sqrt{\frac{m_A E_R}{2\mu_A^2}}
\]
[1310.8214].

The halo assumptions are stated explicitly in the early WIMP-search literature. The first-results paper uses
\[
v_0 = 220~\mathrm{km/s},\qquad v_{\mathrm{esc}} = 544~\mathrm{km/s},\qquad \rho_0 = 0.3~\mathrm{GeV/cm^3},\qquad v_E = 245~\mathrm{km/s}
\]
[1310.8214]. The reanalysis preserves the same Standard Halo Model structure and again lists \(\rho_0 = 0.3\ \mathrm{GeV/cm^3}\), \(v_0 = 220\ \mathrm{km\,s^{-1}}\), \(v_{\mathrm{esc}} = 544\ \mathrm{km\,s^{-1}}\), and \(v_E = 245\ \mathrm{km\,s^{-1}}\) [1512.03506]. The combined WIMP-search paper states that LUX adopts the same Standard Halo Model assumptions used by XENON and CDMS, with \(v_E \approx 245\ \mathrm{km/s}\) for WS2013 and \(v_E \approx 230\ \mathrm{km/s}\) for WS2014–16, while other SHM parameters are standard in the cited works but not explicitly listed there [1710.03572].

For spin-dependent scattering, the combined analysis summarizes the common LUX SD formalism in terms of nuclear structure functions \(S_{ij}(q)\):
\[
\sigma_A^{SD}(q) = \frac{4\mu_A^2}{2J+1}\left[S_{00}(q)a_0^2 + S_{11}(q)a_1^2 + S_{01}(q)a_0 a_1\right]
\]
[1710.03572]. The same source notes that neutron-only and proton-only limits are quoted as \(\sigma_n^{SD}\) and \(\sigma_p^{SD}\) [1710.03572].

The proceedings summary presents the scattering-rate formalism in a compact form and emphasizes that LUX combined standard astrophysical assumptions, isospin-invariant couplings, and in situ detector calibrations in the same pipeline [1611.05525]. This suggests that the formalism’s defining feature is not any single interaction model, but rather the reproducible passage from halo model and cross section to observable-space PDFs.

## 5. Profile-likelihood inference, backgrounds, and datasets

The statistical core of the LUX formalism is an extended profile-likelihood ratio analysis with nuisance-parameter constraints. In the first-results paper, the extended likelihood for \(N\) observed events with signal mean \(\mu\) and background means \(b_k\) is written schematically as
\[
L(\mu,\{b_k\},\boldsymbol{\theta}) = e^{-(\mu + \sum_k b_k)}\prod_{i=1}^N\left[\,\mu\,s(x_i\mid\boldsymbol{\theta}) + \sum_k b_k\,b_k(x_i\mid\boldsymbol{\theta})\right]\times \prod_{j}\,G(b_j\mid \bar{b}_j,\sigma_j)
\]
[1310.8214]. One-sided \(90\%\) confidence limits were then set on the SI WIMP–nucleon cross section by profiling the likelihood over nuisance parameters [1310.8214].

The 2013 reanalysis likewise used a double-sided profile-likelihood ratio. For each WIMP mass, an unbinned extended likelihood over \((cS1, cS2, r, z)\) was built as a mixture of signal and background PDFs with Gaussian-constrained nuisance parameters, including the NR-response parameters \(k\) and \(g_{2,\mathrm{dd}}/g_{2,\mathrm{ws}}\), as well as background normalizations [1512.03506]. Monte Carlo pseudo-experiments via RooStats were used to construct \(90\%\) confidence intervals, and a power constraint at the median expected limit was applied to avoid over-exclusion due to downward background fluctuations [1512.03506].

The combined WS2013 and WS2014–16 analysis used an unbinned PLR across 17 exposure segments, with WS2013 as a 17th segment and WS2014–16 divided into 16 time–drift bins [1710.03572]. The generic likelihood form is given as
\[
L(\mu_s,\{\mu_{b,k}\},\theta) = \prod_{i=1}^{N}\left[\mu_s S(x_i|\theta) + \sum_k \mu_{b,k} B_k(x_i|\theta)\right] e^{-(\mu_s+\sum_k \mu_{b,k})}\prod_j G(\theta_j;\theta_{j,0},\sigma_{\theta_j})
\]
[1710.03572].

The background model is a major part of the formalism. In the 2013 reanalysis, the ER background model included gamma rays, beta decays, \(^{127}\)Xe, \(^{37}\)Ar, and a new empirical model for wall-originating events that enabled an enlarged fiducial radius to 20 cm [1512.03506]. The fitted background normalizations are reported explicitly, including \(\mu_{\gamma,\mathrm{bottom}} = 165 \pm 16\), \(\mu_{\gamma,\mathrm{rest}} = 228 \pm 19\), \(\mu_\beta = 84 \pm 15\), \(\mu_{\mathrm{Xe\text{-}127}} = 78 \pm 12\), \(\mu_{\mathrm{Ar\text{-}37}} = 12 \pm 8\), and \(\mu_{\mathrm{wall}} = 22 \pm 4\) [1512.03506]. In the sub-GeV analysis, six background components were constrained with Gaussian priors: low-z gamma rays, other gamma rays, beta particles, \(^{127}\)Xe, \(^{37}\)Ar, and wall events [1811.11241].

Several concrete datasets define the canonical LUX implementations. The first-results search used \(85.3\) live days with a fiducial mass \(118.3 \pm 6.5\) kg and yielded 160 events after cuts in the ROI [1310.8214]. The 2013 reanalysis used \(95.0\) live days, fiducial mass \(145.4 \pm 1.3\) kg, and 591 events after cuts [1512.03506]. The combined WIMP-search paper adds the WS2014–16 run with 332 live days and 1,221 events after cuts, for a total WS2014–16 exposure of \(95.8\) kg·live-years [1710.03572].

## 6. Extension to sub-GeV dark matter: Bremsstrahlung and Migdal channels

The sub-GeV extension of the LUX formalism addresses a regime in which ordinary elastic nuclear recoils are often invisible in xenon because the deposited NR energy is below threshold [1811.11241]. The key observation is that the tree-level DM–nucleus scattering diagram can be accompanied by Bremsstrahlung photon emission or by the Migdal effect, producing an electron-recoil component at higher observable energy than the NR alone [1811.11241].

LUX modeled these channels by starting from the standard SI DM–nucleus scattering framework and weighting the differential NR cross section by the probability of an accompanying ER process [1811.11241]. The schematic rate is given as
\[
R = \frac{\rho_\chi}{m_\chi}\int d^3v\, f(\mathbf{v})\, v \int dE_R\, \frac{d\sigma_{\mathrm{NR}}}{dE_R}(v,E_R)\,P_{\text{sig}}(E_{\mathrm{obs}}\mid E_R,\text{process})
\]
[1811.11241].

For Bremsstrahlung, the calculation follows Kouvaris and Pradler, with
\[
\frac{d\sigma_{\mathrm{br}}}{d\omega\, dE_R} = \frac{d\sigma_{\mathrm{NR}}}{dE_R}(v,E_R)\times \mathcal{B}(Z,q,\omega)
\]
in a soft-photon approximation [1811.11241]. For the Migdal effect, LUX followed Ibe et al. and included only ionization, omitting excitation and excluding valence-electron shells \(n=5\) because liquid-phase effects may perturb valence spectra [1811.11241]. The Migdal observable energy is
\[
E_{\mathrm{ER}} = E_e + E_{\mathrm{binding}}
\]
and the ionization probability density is summed shell by shell [1811.11241].

This extension also introduced explicit mediator classes. Four mediator classes were tested: heavy or light, each with scalar or vector couplings [1811.11241]. Scalar mediator couplings scale coherently as \(A^2\), while vector mediator couplings scale as \(Z^2\) [1811.11241]. The mediator form factor is parameterized as
\[
F_{\mathrm{med}} \simeq 1 \quad (m_{\mathrm{med}} \gg q),
\]
and
\[
F_{\mathrm{med}} = \frac{q_{\mathrm{ref}}^4}{q^4}\quad (m_{\mathrm{med}} \ll q),\qquad q_{\mathrm{ref}} = 1~\mathrm{MeV}
\]
[1811.11241].

The experimental motivation is directly tied to thresholds. LUX had \(50\%\) detection efficiency at \(1.24\) keV for ERs versus \(3.3\) keV for NRs, so ER-side channels opened sensitivity to \(m_\chi\) below about \(5\ \mathrm{GeV}/c^2\) where NR-only searches are inefficient [1811.11241]. The sub-GeV search used the 2013 WS2013 dataset with \(95\) live days and search exposure \(1.4\times 10^4\) kg·day, and set constraints on SI DM–nucleon scattering for masses \(0.4\)–\(5\ \mathrm{GeV}/c^2\) [1811.11241]. The observed event counts were consistent with the background-only hypothesis for all tested masses and mediator classes [1811.11241].

A plausible implication is that this sub-GeV analysis did not replace the canonical LUX formalism; it reused the same detector-response and PLR machinery while substituting ER-band signal models for otherwise invisible NR interactions.

## 7. Results, scope, and interpretive boundaries

Across its main WIMP analyses, the LUX formalism yielded successively stronger \(90\%\) C.L. limits under standard halo assumptions. The first-results paper reported a minimum SI upper limit of \(7.6 \times 10^{-46}\ \mathrm{cm}^2\) at \(m_\chi = 33\ \mathrm{GeV}/c^2\) [1310.8214]. The 2013 reanalysis improved this to \(0.6\) zb, explicitly \(6\times 10^{-46}\ \mathrm{cm}^2\), again at \(33\ \mathrm{GeV}/c^2\) [1512.03506]. The combined WS2013+WS2014–16 analysis reported a minimum SI limit of \(0.11\times 10^{-45}\ \mathrm{cm}^2\), equivalently \(1.1\times 10^{-46}\ \mathrm{cm}^2\), at \(50\ \mathrm{GeV}/c^2\), and spin-dependent minima of \(1.6\times 10^{-41}\ \mathrm{cm}^2\) for neutron-only coupling and \(5.0\times 10^{-40}\ \mathrm{cm}^2\) for proton-only coupling at \(35\ \mathrm{GeV}/c^2\) [1710.03572]. The sub-GeV search extended sensitivity to \(0.4\)–\(5\ \mathrm{GeV}/c^2\) using Bremsstrahlung and Migdal ER signatures [1811.11241].

The formalism’s scope is defined by explicit assumptions and limitations. LUX’s principal WIMP analyses adopted the Standard Halo Model and did not profile astrophysical uncertainties in the first-results paper [1310.8214]. The reanalysis states that uncertainties in the NR-response nuisance parameters changed the limit by less than \(20\%\) relative to fixing the best-fit model [1512.03506]. The sub-GeV paper notes that variations in \(\rho_\chi\), \(v_0\), and \(v_{\mathrm{esc}}\) can shift limits modestly, that light-mediator definitions can differ slightly between experiments, and that omitting valence-shell Migdal contributions makes the limits conservative [1811.11241].

The literature also contains reinterpretive analyses that stress the dependence of low-mass exclusions on low-energy liquid-xenon response assumptions. “Dark Matter in Light of LUX” studies exothermic scattering, isospin-dependent couplings, halo-independent comparisons, and variations in LXe scintillation assumptions, concluding that only a highly tuned isospin-dependent scenario remains viable unless LXe scintillation properties are dramatically different from those assumed by LUX [1401.0216]. This does not revise the collaboration’s formalism, but it identifies where sensitivity to modeling assumptions is concentrated.

Taken together, the LUX formalism is best understood as a detector-calibrated, likelihood-based translation layer between dark-matter theory and xenon-TPC observables. Its stable components are the S1/S2 signal model, in situ ER and NR calibration, nuisance-constrained PLR inference, and Standard Halo Model rate calculation. Its adaptable components are the interaction model and the choice of signal channel: elastic NRs for canonical WIMPs, or ER-band Bremsstrahlung and Migdal signatures for sub-GeV dark matter [1710.03572; 1608.05381; 1811.11241].

Source: https://www.emergentmind.com/topics/lux-formalism