⚛️ phenomenology

Coherence from interference: a solvable model of sub-GeV dark matter-nucleus scattering

This paper employs an exactly solvable 1D crystal lattice model to demonstrate that the transition from coherent to incoherent dark matter-nucleus scattering is governed by the diminishing importance of crystal momentum conservation as multiphonon production increases, thereby validating the use of incoherent approximations for sub-GeV dark matter detection in realistic 3D crystals.

Original authors: Lynn Lin, Tongyan Lin, Momei Fang

Published 2026-07-29

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Lynn Lin, Tongyan Lin, Momei Fang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Invisible Ghost and the Crystal Ball

Imagine the universe is filled with a mysterious, invisible substance called dark matter. We know it's there because its gravity holds galaxies together, but we've never actually seen a single particle of it. For decades, scientists have been trying to catch these particles by building giant, ultra-sensitive detectors deep underground. They wait for a dark matter particle to bump into an atom in the detector, hoping to see a tiny flash of energy. This works great if the dark matter is heavy, like a bowling ball hitting a pin. But what if the dark matter is incredibly light, like a tiny speck of dust? A speck of dust hitting a pin doesn't knock it over; it just makes the pin vibrate.

This is the challenge for "light" dark matter. When these tiny particles hit a crystal detector, they don't knock the whole atom flying. Instead, they make the entire crystal lattice—a perfectly ordered grid of atoms—wiggle. In physics, these wiggles are called "phonons," which are basically sound waves traveling through the solid material. The big question scientists are trying to answer is: How do we calculate exactly how much the crystal wiggles? Do we treat the crystal as a single, unified drum that vibrates in perfect harmony (coherent scattering), or do we treat it as a bag of individual marbles that bump into each other randomly (incoherent scattering)? Getting this math right is crucial because if we get it wrong, we might miss the dark matter signal entirely or mistake background noise for a discovery.

The Crystal Orchestra vs. The Soloist

In this paper, the authors tackle this tricky math problem using a clever trick: they shrink the universe down to a one-dimensional line of atoms. Imagine a long row of identical beads connected by springs, like a necklace. They use this simple "toy model" to simulate how a dark matter particle would scatter off the crystal. Because the model is so simple, they can solve the equations exactly, without needing to make messy guesses. This allows them to see the "interference" effects—the way the waves from different atoms overlap and either amplify or cancel each other out.

The authors discovered a fascinating rule about how these waves behave. When a dark matter particle hits the crystal and creates just one phonon (a single ripple), the atoms must act in perfect unison. It's like a choir singing a single note; every singer must hit the exact right pitch and timing, or the sound doesn't work. In physics terms, this is called "coherent scattering," and it requires strict conservation of "crystal momentum." The atoms are locked together by the rules of the crystal lattice, and they can only vibrate in specific, allowed patterns.

However, the story changes when the dark matter particle hits hard enough to create two or more phonons. Suddenly, the strict rules loosen up. It's as if the choir breaks into a chaotic jam session. The particles can create ripples in many different combinations, and the strict requirement for everyone to be perfectly synchronized starts to fade. The authors found that when you have multiple phonons, the complex interference patterns that make the math so hard start to cancel each other out. The result? You can stop worrying about the perfect harmony of the whole crystal and just treat the atoms as individual, independent bumpers. This is called the "incoherent approximation," and it is much, much easier to calculate.

The "Hybrid" Solution

The paper's main finding is a validation of a "hybrid" strategy for scientists. They showed that for the hardest part of the calculation—where the dark matter is light and only creates a single ripple—you must use the complex, exact math that accounts for the crystal's perfect harmony. But as soon as the dark matter creates two or more ripples, or if the energy is high enough, you can safely switch to the simpler "incoherent" math.

To prove this, they ran simulations on their 1D crystal model. They compared the "exact" calculation (the hard way) against the "incoherent" approximation (the easy way) and a new "hybrid" method that switches between them. The results were clear:

For single phonons: The simple approximation failed miserably, missing the specific energy peaks where the crystal actually vibrates.
For two or more phonons: The simple approximation was spot-on. The complex interference effects had washed out, and the easy math worked perfectly.
The Hybrid Method: By using the hard math for single phonons and the easy math for everything else, they created a method that is both fast and accurate.

They tested this against different types of dark matter interactions (some with heavy "messengers" and some with massless ones) and found that this hybrid approach reproduces the exact results with very little error—usually less than 30%, and often just a few percent. This is a huge win because calculating the exact math for real, 3D crystals is computationally impossible for many scenarios. By proving that the "easy" way works for most of the action, the authors have given experimentalists a reliable, fast tool to predict what their detectors should see.

Why It Matters

This isn't just about solving a math puzzle; it's about keeping the hunt for dark matter on track. As detectors become more sensitive, they are starting to look for these tiny, single-phonon signals. If scientists used the wrong math, they might look in the wrong place or misinterpret the data. This paper provides a clear roadmap: use the detailed, complex model when the crystal is singing a single note, but switch to the simpler, faster model when the crystal starts jamming with multiple notes. It turns a computationally prohibitive nightmare into a manageable calculation, bringing us one step closer to finally catching that invisible ghost.

Technical Summary: Coherence from Interference: A Solvable Model of Sub-GeV Dark Matter-Nucleus Scattering

Problem Statement
Direct detection of sub-GeV dark matter (DM) requires a theoretical framework capable of describing scattering off crystal targets where momentum transfer ( $q$ ) and energy deposition ( $\omega$ ) are small. In this regime, the approximation of elastic free nuclear recoils fails, and DM interacts with the crystal lattice to produce collective excitations, specifically phonons. While single-phonon scattering is well-understood, the transition to multiphonon excitations and the eventual limit of incoherent nuclear recoils remains computationally challenging. Existing approaches often rely on the "incoherent approximation," which neglects interference terms between different lattice sites to simplify calculations of the dynamic structure factor $S(q, \omega)$ . However, the accuracy of this approximation, particularly regarding the transition from coherent to incoherent regimes and its validity at low $q$ , has not been rigorously quantified.

Methodology
The authors employ a 1D $N$ -site crystal lattice model to analytically and numerically compute the full dynamic structure factor exactly. This model allows for the explicit separation of interference terms (coherent scattering) from auto-correlation terms (incoherent scattering).

Toy Model (2 Sites): The authors first analyze a two-site system to demonstrate how interference terms manifest as oscillatory $\cos(qa)$ factors in the structure factor. They show that for single-phonon processes, these terms enforce momentum conservation, suppressing scattering into optical modes at low $q$ . For higher-order processes ( $n \ge 2$ ), alternating phase factors in the interference terms lead to cancellations when summing over states.
Exact 1D $N$ -Site Solution: The model is extended to an $N$ -site periodic lattice. The authors derive the exact structure factor by summing over all phonon modes and lattice sites. They utilize a recursive relation (Appendix A) to efficiently compute $n$ -phonon contributions for large $N$ .
Comparison of Regimes: The full structure factor is compared against the incoherent approximation (where cross-terms $\ell \neq \ell'$ are dropped) and the impulse approximation (valid at very high $q$ ).
Scattering Rate Calculation: The authors compute differential scattering rates ($dR/dq$ and $dR/d\omega$ ) for both massive and massless scalar mediators, applying the derived structure factors to a 1D model serving as a proxy for isotropic 3D crystals.

Key Contributions and Results

Mechanism of Coherence: The paper establishes that the sole difference between the full (coherent) structure factor and the incoherent approximation is the enforcement of crystal momentum conservation via a delta function constraint ( $q = \sum k_i + G$ ). In the incoherent approximation, this constraint is removed, and the response is smeared over the phonon density of states.
Transition to Incoherence: The study demonstrates that as the number of produced phonons ( $n$ ) increases, the momentum conservation constraint becomes a weak restriction on the possible combinations of phonon energies and momenta. Consequently, the interference terms (which enforce this constraint) average out, and the full structure factor converges to the incoherent approximation.
Validity of Approximations:
- Single Phonon ( $n=1$ ): The incoherent approximation fails significantly at low $q$ ( $q < q_{BZ}$ ), missing the sharp dispersion relation and kinematic thresholds. It incorrectly predicts a non-zero response below the single-phonon energy threshold.
- Multiphonon ( $n \ge 2$ ): For two or more phonons, the incoherent approximation becomes increasingly accurate. Numerical results show that the oscillatory features of the full structure factor dampen rapidly with increasing $n$ , and the integrated rates converge.
Hybrid Scheme (Inc + LW): The authors propose and validate a "Combined Incoherent + Long-Wavelength" (Inc + LW) approximation. This scheme uses the exact (or long-wavelength) coherent single-phonon result for $q < q_{BZ}$ and switches to the computationally efficient incoherent approximation for $n \ge 2$ or $q > q_{BZ}$ .
Quantitative Validation: Using the Inc + LW scheme, the authors calculate the DM-nucleon cross-section required to achieve 3 events/kg-year. They find that for massive mediators, the maximum error relative to the full calculation does not exceed 30% (and is often percent-level) across most of the sub-GeV mass range. For massless mediators, the agreement is similarly robust, except near kinematic thresholds.

Significance and Claims
The paper claims to provide the first rigorous, solvable model that explicitly identifies the role of interference in the transition from coherent to incoherent dark matter scattering. By isolating momentum conservation as the defining feature of coherence, the authors generalize this intuition beyond 1D to realistic 3D crystals.

The primary significance lies in validating the use of the incoherent approximation for multiphonon processes, which are computationally prohibitive to calculate exactly in 3D. The authors assert that their hybrid Inc + LW scheme offers a robust, accurate, and computationally efficient method for predicting scattering rates in direct detection experiments targeting sub-GeV dark matter. This allows experimentalists to model detector sensitivities for MeV-scale dark matter candidates without needing to perform exact, expensive multiphonon calculations for the entire parameter space, provided they correctly treat the single-phonon coherent regime at low momentum transfer. The paper concludes that while the analytical derivation is specific to a 1D model, the fundamental kinematic behaviors regarding the suppression of interference in multiphonon regimes are expected to persist in 3D targets.

The Invisible Ghost and the Crystal Ball

The Crystal Orchestra vs. The Soloist

The "Hybrid" Solution

Why It Matters

More like this