Dark Matter’s Hidden Boost: How Unstable Particles Amplify Annihilation

Author: Denis Avetisyan


New research reveals that interactions between the products of dark matter annihilation can significantly increase the rate at which dark matter particles disappear, impacting our understanding of its abundance in the universe.

The annihilation cross section of <span class="katex-eq" data-katex-display="false">\chi_{1}\bar{\chi}_{1}</span> particles-calculated with parameters including <span class="katex-eq" data-katex-display="false">m_{1} = 1~\text{TeV}</span>, <span class="katex-eq" data-katex-display="false">m_{2}/m_{1} = 1.01</span>, and <span class="katex-eq" data-katex-display="false">\alpha = 0.2</span>-exhibits a pronounced threshold effect at <span class="katex-eq" data-katex-display="false">E_{2} = 0</span>, with the inclusion of final-state Sommerfeld enhancements-modeled via a cutoff method and defined by equation (29)-significantly altering the cross section compared to calculations neglecting these enhancements and demonstrating a clear dependence on the strength of the Coulomb potential.
The annihilation cross section of \chi_{1}\bar{\chi}_{1} particles-calculated with parameters including m_{1} = 1~\text{TeV}, m_{2}/m_{1} = 1.01, and \alpha = 0.2-exhibits a pronounced threshold effect at E_{2} = 0, with the inclusion of final-state Sommerfeld enhancements-modeled via a cutoff method and defined by equation (29)-significantly altering the cross section compared to calculations neglecting these enhancements and demonstrating a clear dependence on the strength of the Coulomb potential.

This study demonstrates a final-state Sommerfeld enhancement effect arising from long-range interactions, altering predicted relic densities and offering new avenues for dark matter detection.

The standard calculation of dark matter annihilation often assumes instantaneous interactions, potentially overlooking crucial final-state effects. This paper, ‘Sommerfeld enhancement from unstable final-state particles in dark matter annihilation’, investigates the impact of long-range interactions between unstable annihilation products on the overall process. We demonstrate that these interactions can significantly enhance the annihilation cross-section via resonant effects and, importantly, naturally incorporate contributions from off-shell final states. How do these final-state interactions refine our understanding of dark matter’s relic abundance and its potential detection signals?


The Unseen Universe: Defining the Dark Matter Abundance

Observations across multiple scales – from galactic rotation curves to the cosmic microwave background – consistently indicate that roughly 85% of the universe’s matter content is non-baryonic, existing as what scientists term “dark matter”. Despite this overwhelming evidence for its existence, the fundamental composition of dark matter remains one of the most profound mysteries in modern cosmology. It doesn’t interact with light, making direct observation impossible, and its interactions with ordinary matter appear extraordinarily weak. Leading candidates range from weakly interacting massive particles (WIMPs) and axions to primordial black holes, but decades of dedicated searches have yet to definitively identify the particle or object comprising this invisible mass. This elusiveness drives ongoing research, pushing the boundaries of both theoretical physics and experimental techniques in the quest to unveil the true nature of dark matter and its role in the universe’s formation and evolution.

Establishing the precise relic abundance of dark matter – the amount remaining from the early universe – serves as a critical benchmark for theoretical physicists. These models, attempting to define the nature of this invisible substance, predict a specific quantity of dark matter that should have survived the intense conditions following the Big Bang. A calculated abundance that drastically diverges from observational estimates would invalidate the underlying assumptions of the model, necessitating a re-evaluation of its fundamental parameters or even its core principles. Conversely, a strong match between theoretical predictions and observed abundance provides compelling evidence supporting a particular dark matter candidate and guides future research efforts towards refining and testing that specific hypothesis. Essentially, the relic abundance acts as a powerful constraint, shaping the landscape of dark matter research and pushing the boundaries of our understanding of the cosmos.

Estimating the amount of dark matter present in the universe – its relic abundance – fundamentally depends on understanding how effectively dark matter particles annihilate each other. This annihilation rate is quantified by the annihilation cross section, a measure of the probability that two dark matter particles will collide and destroy each other, producing standard model particles. Current calculations assume a specific annihilation cross section, and even slight inaccuracies in this value can dramatically alter the predicted dark matter abundance. Consequently, refining the precision of the annihilation cross section – through both direct detection experiments and indirect searches for annihilation products – is paramount for confirming or refuting various dark matter theories. The challenge lies in the fact that the predicted cross sections vary wildly depending on the specific dark matter model, creating a significant uncertainty in determining the true relic abundance and, ultimately, the nature of this mysterious substance.

For the attractive Hulthén potential with <span class="katex-eq" data-katex-display="false">m_1 = 1~\text{TeV}</span>, <span class="katex-eq" data-katex-display="false">\alpha = 0.2</span>, and <span class="katex-eq" data-katex-display="false">a = 10^{-7}~\text{GeV}^{-2}</span>, the relic abundance <span class="katex-eq" data-katex-display="false">\Omega h^2</span> varies with the mass ratio <span class="katex-eq" data-katex-display="false">m_2/m_1</span>, exhibiting values consistent with observed dark matter density (represented by the blue line) for <span class="katex-eq" data-katex-display="false">m_<i> = 0.1~\text{TeV}</span> (top) and <span class="katex-eq" data-katex-display="false">m_</i> = 0.3~\text{TeV}</span> (bottom), with a magnified view highlighting the region <span class="katex-eq" data-katex-display="false">1 \leq m_2/m_1 \leq 1.01</span>.
For the attractive Hulthén potential with m_1 = 1~\text{TeV}, \alpha = 0.2, and a = 10^{-7}~\text{GeV}^{-2}, the relic abundance \Omega h^2 varies with the mass ratio m_2/m_1, exhibiting values consistent with observed dark matter density (represented by the blue line) for m_<i> = 0.1~\text{TeV} (top) and m_</i> = 0.3~\text{TeV} (bottom), with a magnified view highlighting the region 1 \leq m_2/m_1 \leq 1.01.

Beyond Simplification: Refining the Annihilation Calculation

The Cutoff Method, a common simplification in calculating annihilation rates, introduces significant inaccuracies by artificially limiting the momentum transfer in the interaction. This limitation effectively truncates the integral representing the annihilation process, neglecting contributions from particles with higher energies and momenta. Consequently, the calculated annihilation rate is consistently lower than the true value, with the degree of underestimation dependent on the specific particles involved and the chosen cutoff value. More accurate calculations necessitate including the full momentum space integration, accounting for the complete range of possible interaction energies, and avoiding such arbitrary limitations on the phase space.

Accurate calculation of annihilation rates necessitates the inclusion of off-shell particles and Sommerfeld enhancement effects, particularly when dealing with non-relativistic particles. Traditional calculations often assume particles are ‘on-shell’ – meaning they satisfy the energy-momentum relation E^2 = p^2c^2 + m^2c^4. However, intermediate virtual particles in the annihilation process can be off-shell, altering the cross-section. Sommerfeld enhancement, a non-perturbative effect, arises from the long-range nature of interactions – such as the Coulomb interaction – between the annihilating particles. This effect significantly boosts the annihilation rate at low velocities by increasing the phase space volume available for annihilation, and requires detailed consideration of the initial state wave function and binding energy to accurately model the enhancement factor.

Sommerfeld Enhancement, a non-perturbative effect, significantly increases the annihilation rate of dark matter particles by accounting for the distortion of wavefunctions due to long-range interactions, specifically the Coulomb interaction. This effect is prominent at low velocities where particles can be deflected into bound states before annihilation. Accurate calculation necessitates detailed modeling of the initial state, including the momentum distribution and the potential describing the interaction; simplified treatments that ignore the full momentum spectrum or potential shape will underestimate the enhancement. The magnitude of Sommerfeld Enhancement is dependent on the strength of the interaction and can range from a modest increase to orders of magnitude enhancement of the expected annihilation cross-section, crucially impacting indirect detection rates of dark matter.

The thermally averaged annihilation cross section, plotted as a function of <span class="katex-eq" data-katex-display="false">x = m_1/T</span>, demonstrates strong agreement between the full result from equation (39), the free result with <span class="katex-eq" data-katex-display="false">S_f(E_2, \Gamma) = 1</span>, and the cutoff method using equation (41) with <span class="katex-eq" data-katex-display="false">v_{cut} = 10^{-2}</span>.
The thermally averaged annihilation cross section, plotted as a function of x = m_1/T, demonstrates strong agreement between the full result from equation (39), the free result with S_f(E_2, \Gamma) = 1, and the cutoff method using equation (41) with v_{cut} = 10^{-2}.

Resonances Unveiled: The Impact of Final State Interactions

The Final State Enhancement (FSE) addresses a phenomenon where the annihilation cross section is modified by forces acting between the particles created in the annihilation process. Unlike traditional calculations which treat annihilation products as free particles immediately after creation, FSE accounts for residual, long-range interactions – typically electromagnetic or strong force interactions – that can bind or otherwise influence the final state particles. These interactions are not necessarily strong enough to form a stable, fully bound state, but they can create a temporary, resonant state which significantly alters the probability of the annihilation occurring at a given energy. The effect is particularly pronounced when the final state particles have a low relative velocity, increasing the time they spend interacting before fully separating.

Determination of the Final State Enhancement (FSE) requires a quantum mechanical treatment of the interacting annihilation products via the Schrödinger Equation. Specifically, solving the time-independent Schrödinger Equation, H\Psi = E\Psi, yields the bound state energies (E) and corresponding wavefunctions (Ψ) of the final state particles. From these wavefunctions, the decay width (Γ) – inversely proportional to the lifetime of the bound state – can be calculated using time-dependent perturbation theory or related methods. Accurate determination of both the bound state energy and decay width is essential for quantifying the resonant enhancement of the annihilation cross section; the magnitude of the enhancement is directly related to the proximity of a pole to the real axis in the complex momentum space, which is determined by these parameters.

Accurate modeling of final state interactions requires the use of realistic potential functions beyond simple approximations. Potentials like the Hulthen potential, which incorporates screening effects, provide a more faithful representation of the forces between annihilation products. This is because the inclusion of these potentials allows for the formation of bound states, leading to resonance phenomena that significantly alter the annihilation cross section. Specifically, the annihilation cross section can be enhanced by a factor approaching ~10 when these bound state resonances occur, a result not achievable with less sophisticated potential models. The precise magnitude of this enhancement is directly dependent on the parameters of the chosen potential and the specific interacting particles.

Annihilation cross sections of <span class="katex-eq" data-katex-display="false">\chi_1 \bar{\chi}_1</span> as a function of <span class="katex-eq" data-katex-display="false">E_2</span> under the attractive Hulthén potential, calculated with varying <span class="katex-eq" data-katex-display="false">m_*</span> values (0.1 TeV left, 0.3 TeV right) and showing the impact of final-state suppression effects on the cross section.
Annihilation cross sections of \chi_1 \bar{\chi}_1 as a function of E_2 under the attractive Hulthén potential, calculated with varying m_* values (0.1 TeV left, 0.3 TeV right) and showing the impact of final-state suppression effects on the cross section.

Constraining the Invisible: Implications for Dark Matter Models

The predicted abundance of dark matter remaining in the universe today, known as the relic abundance, is a cornerstone for evaluating the plausibility of different theoretical models. Calculations that omit final-state enhancement (FSE) – the boosting of annihilation rates when dark matter particles decay into multiple final-state products – can yield dramatically inaccurate predictions. Incorporating FSE effects often leads to a substantial alteration of the relic abundance, potentially shifting a previously viable dark matter candidate into an excluded region, or conversely, rescuing a seemingly improbable model. This is particularly true when the masses of the dark matter particle and its annihilation products are closely matched, creating resonant enhancement of the annihilation cross-section and influencing the overall dark matter density. Consequently, a precise accounting of FSE is not merely a refinement, but a fundamental necessity for distinguishing between competing dark matter hypotheses and accurately mapping the landscape of possible solutions to the mystery of dark matter.

The influence of final-state enhancement (FSE) becomes critically important when considering dark matter models built upon composite particles, such as Accidental Composite Dark Matter. These models posit that dark matter isn’t elementary, but rather bound states formed from new, strongly interacting particles. The presence of these internal constituents dramatically alters the decay and annihilation pathways of the dark matter candidate, and FSE – the boosting of annihilation rates due to the production of resonant final states – can significantly enhance these processes. Without accurately accounting for FSE, predictions for the relic abundance of these composite candidates can be severely skewed, potentially leading to the dismissal of viable models. The effect is particularly pronounced when the mass ratio of the constituent particles is close to unity, leading to resonant behavior and a substantial alteration in the predicted dark matter density.

Determining the true nature of dark matter necessitates increasingly precise calculations of its relic abundance-the amount remaining from the early universe. Recent work highlights the crucial role of incorporating both Sommerfeld enhancements, arising from long-range forces between dark matter particles, and final-state enhancement (FSE), which accounts for the impact of particle decays or annihilations into nearly identical final states. The inclusion of FSE can dramatically alter predicted relic abundances, particularly for composite dark matter models, revealing substantial deviations from calculations that neglect these final-state effects. Specifically, calculations indicate a significant resonance-a bound state-occurs when the mass ratio of the interacting particles is approximately 1.001, exhibiting a binding energy of -2.6 \text{ GeV}. This resonance significantly boosts the annihilation rate, influencing the predicted dark matter density and providing a critical benchmark for distinguishing between viable and excluded dark matter candidates.

The dark matter relic abundance <span class="katex-eq" data-katex-display="false">\Omega h^2</span> is shown as a function of the mass ratio <span class="katex-eq" data-katex-display="false">m_2/m_1</span>, demonstrating that the calculated values closely align with the observed dark matter energy density (blue) when using a cutoff of <span class="katex-eq" data-katex-display="false">v_{cut} = 10^{-2}</span> or a damping factor of <span class="katex-eq" data-katex-display="false">\Gamma/m_2 = 10^{-4}</span>, particularly within the mass ratio range of <span class="katex-eq" data-katex-display="false">1 \leq m_2/m_1 \leq 1.01</span>.
The dark matter relic abundance \Omega h^2 is shown as a function of the mass ratio m_2/m_1, demonstrating that the calculated values closely align with the observed dark matter energy density (blue) when using a cutoff of v_{cut} = 10^{-2} or a damping factor of \Gamma/m_2 = 10^{-4}, particularly within the mass ratio range of 1 \leq m_2/m_1 \leq 1.01.

The study meticulously strips away extraneous assumptions regarding dark matter annihilation, revealing a more nuanced understanding of relic abundance. It posits that final-state interactions, previously considered negligible, exert a demonstrable influence on the annihilation cross-section. This echoes a sentiment articulated by Georg Wilhelm Friedrich Hegel: “The truth is the whole.” The research doesn’t seek to add complexity, but to account for all contributing factors-even those initially obscured by simplification. By considering the Sommerfeld enhancement arising from unstable final-state particles, the work moves closer to a complete picture of dark matter’s properties, a truth revealed not through accretion, but through careful subtraction of falsehoods.

Where Do We Go From Here?

This work clarifies a simple point: final-state interactions matter. The Sommerfeld enhancement, once a calculation of convenience, reveals itself as a consequence of physics, not merely a mathematical trick. Abstractions age, principles don’t. The relic abundance, that cornerstone of dark matter searches, demands scrutiny beyond simplistic assumptions.

Limitations remain. The treatment of bound states, while improved, still relies on approximations. Every complexity needs an alibi. A full, non-perturbative understanding of these final-state effects is necessary. Future work must address the interplay between these interactions and the broader cosmological landscape.

The path forward isn’t more parameters, but sharper principles. A deeper investigation into non-relativistic scattering and the formation of long-lived resonances is essential. This isn’t about finding a dark matter candidate. It’s about understanding the fundamental physics that governs its annihilation.


Original article: https://arxiv.org/pdf/2603.02647.pdf

Contact the author: https://www.linkedin.com/in/avetisyan/

See also:

2026-03-05 05:49