⚛️ phenomenology

Solar Neutrino Signals and Capture-Annihilation Equilibrium of SU(5) Neutralino Dark Matter

Motivated by IceCube solar searches, this paper investigates neutralino dark matter in supersymmetric SU(5) with nonuniversal gaugino masses to establish a consistent framework linking local abundance, composition, and capture-annihilation equilibrium to solar neutrino and muon signals, thereby providing a tool for interpreting subdominant dark matter candidates against experimental limits.

Original authors: M. Adeel Ajaib, Abdul Rehman

Published 2026-10-20

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: M. Adeel Ajaib, Abdul Rehman

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

Deep in the universe, vast clouds of invisible matter drift between the stars. Scientists call this dark matter, and while we cannot see it directly, its gravity shapes the galaxies we observe. One leading idea suggests that this dark matter is made of tiny, heavy particles called neutralinos. These particles are their own antiparticles, meaning that when two of them collide, they can annihilate each other and release a burst of energy. If these particles drift into our Sun, the Sun's immense gravity can trap them. Over billions of years, they sink to the core, collide, and vanish, producing a stream of ghostly particles called neutrinos that fly out into space.

Detecting these neutrinos is a primary goal for the IceCube Neutrino Observatory, a massive detector buried deep in the ice of Antarctica. By looking for neutrinos coming from the direction of the Sun, scientists hope to catch a glimpse of dark matter in action. However, a crucial question remains: how do we interpret what we see? If the dark matter particles in our neighborhood are fewer than the maximum amount theory allows, does the signal from the Sun simply get weaker in direct proportion? Or does the physics of how these particles gather and collide inside the Sun change the story entirely?

A team of researchers at Pennsylvania State University and the University of Delaware set out to answer this question by building a detailed model of the Sun's dark matter population. They focused on a specific theoretical framework known as supersymmetric SU(5), which predicts a rich variety of possible dark matter particles. Using powerful computer simulations, they tracked how these particles would behave if they were captured by the Sun, how long it would take for them to reach a steady state of collisions, and what kind of signal they would eventually produce. Their work reveals that the relationship between the amount of dark matter and the signal we see is far more complex than a simple scaling rule.

The researchers began by generating a selected sample of possible versions of the universe, each with different properties for the dark matter particles. In some scenarios, the particles are mostly one type, while in others, they are a mix of different components. They calculated how many of these particles would be captured by the Sun and how quickly they would annihilate. A key finding emerged when they compared two different assumptions about the local density of dark matter. In the first scenario, they assumed the Sun is surrounded by the maximum possible density of these particles. In the second, they assumed the density was lower, reflecting a situation where dark matter makes up only a fraction of the total mass in our region of the galaxy.

When the density is high, the particles capture quickly and reach a steady state where the rate of capture equals the rate of annihilation. In this state, the signal is strong and predictable. However, when the density is lower, the capture rate drops. This reduction has a double effect: fewer particles are caught, and the time it takes for the remaining particles to reach that steady state of collisions becomes much longer. Because the Sun is only about 4.6 billion years old, many of these lower-density scenarios have not had enough time to reach that steady state yet. Consequently, the signal from the Sun is not just weaker because there are fewer particles; it is significantly weaker because the particles have not had time to build up a high rate of collisions.

The study shows that for many of the models they tested, simply assuming a linear drop in the signal based on the drop in density would lead to a massive overestimation of what IceCube should see. In some cases, the actual signal was less than half of what a simple linear calculation would predict. This is particularly true for models where the dark matter particles are dominated by a specific component known as a wino. For these particles, the delay in reaching equilibrium is so pronounced that the signal is heavily suppressed. Other types of particles, such as those dominated by a component called a bino, behave differently and stay closer to the linear prediction, but even they do not follow a simple rule across the board.

To make these findings concrete, the researchers selected five specific examples, or benchmarks, that represent different types of dark matter behavior. One example involved a particle that is almost entirely a wino. In this case, even though the particle is captured by the Sun, the rate of collisions is so slow that the signal remains very faint, far below what one might expect if the density were simply scaled down. Another example featured a particle that is a mix of different components, which showed a moderate suppression of the signal. A third example involved a particle that is almost entirely a Higgsino, which behaved almost exactly as a simple linear scaling would predict, reaching equilibrium quickly enough to maintain a stronger signal relative to its density.

These differences matter deeply for the search for dark matter. The IceCube observatory has set strict limits on how many neutrinos it can see coming from the Sun. If scientists assume a simple linear relationship between density and signal, they might incorrectly conclude that certain types of dark matter are ruled out because the signal is too weak to be seen. However, this new analysis shows that the signal is naturally weaker for many models due to the time it takes to build up. Therefore, a lack of detection does not necessarily mean the dark matter doesn't exist; it might just mean that the particles are still gathering speed, so to speak, and have not yet reached their peak collision rate.

The researchers also compared their predictions with data from direct detection experiments, which look for dark matter particles bumping into atomic nuclei in underground detectors. They found that the rules for these two types of experiments are different. In direct detection, the signal scales linearly with density, so the limits are straightforward. But for the Sun, the signal depends on the history of the particle population. This means that a model could be consistent with direct detection limits but still produce a very faint signal in the Sun, or vice versa. The study provides a new map that connects the local abundance of dark matter to the time it takes to reach equilibrium, allowing scientists to interpret IceCube's data with much greater precision.

Ultimately, this work refines how we listen for the universe's invisible inhabitants. It demonstrates that the Sun is not just a passive trap for dark matter but a dynamic system where the history of capture and the speed of collision matter just as much as the number of particles present. By accounting for the time it takes for these particles to settle into a steady rhythm of annihilation, scientists can avoid overestimating what they should see. This clarity is essential for the next generation of searches, ensuring that when we look at the Sun and see nothing, we understand exactly what that silence means for the nature of dark matter.

Technical Summary: Solar Neutrino Signals and Capture–Annihilation Equilibrium of SU(5) Neutralino Dark Matter

Problem Statement
The paper addresses the interpretation of solar neutrino and muon signals from dark matter (DM) annihilation, specifically within the context of the IceCube neutrino telescope. While the lightest neutralino is a well-motivated DM candidate in supersymmetric (SUSY) theories with conserved R-parity, its abundance and detection signatures depend heavily on its composition (bino, wino, Higgsino) and the surrounding SUSY spectrum.

A critical issue arises when considering neutralinos that are thermally underabundant (i.e., their relic density $\Omega_\chi h^2$ is less than the observed cosmological density). Such particles may exist as a subcomponent of the local DM density. Previous studies often assumed that the solar signal from these subcomponents scales linearly with their local abundance fraction ( $\xi$ ). However, this linear scaling relies on the assumption that the captured neutralino population in the Sun has reached capture–annihilation equilibrium. For subdominant components, the reduced capture rate lengthens the equilibration timescale, potentially pushing the system out of equilibrium. Consequently, the annihilation rate may scale closer to $\xi^2$ rather than $\xi$ , leading to significant overestimates of the expected signal if linear rescaling is applied blindly.

Methodology
The authors investigate this problem using a supersymmetric $SU(5)$ Grand Unified Theory framework with nonuniversal gaugino masses.

Model Setup: The study adopts conventional $SU(5)$ matter assignments at the unification scale, with independent gaugino masses ( $M_1, M_2, M_3$ ) and specific sign conventions ( $M_1, M_2, \mu > 0$ ; $M_3 < 0$ ). The scan covers scalar masses up to 30 TeV and gaugino masses up to 5 TeV.
Numerical Tools:
- ISAJET/ISASUGRA: Used to generate SUSY spectra.
- micrOMEGAs: Calculated the thermal relic abundance ( $\Omega_\chi h^2$ ) and classified neutralino composition (bino-like, wino-like, Higgsino-like, or mixed) based on 90% purity.
- DarkSUSY (v6.4.3): Computed solar capture rates ( $C_\odot$ ), annihilation coefficients ( $C_A$ ), and the resulting neutrino and muon fluxes.
Population Dynamics: The evolution of the number of neutralinos in the Sun ( $N$ ) is governed by the standard population equation: $dN/dt = C_\odot - C_A N^2 - C_E N$ . The authors neglect evaporation ( $C_E$ ) for the retained mass range and assume no self-capture. The solution yields an annihilation rate $\Gamma_A = \frac{1}{2}C_\odot \tanh^2(t_\odot/\tau)$ , where $\tau = (C_\odot C_A)^{-1/2}$ is the equilibration timescale.
Rescaling Procedure: To model subdominant DM, the local density is scaled by $\xi = \Omega_\chi h^2 / 0.12$ . The capture rate scales linearly ( $C_\odot(\xi) = \xi C_\odot(1)$ ), while the timescale scales as $\tau(\xi) = \tau(1)/\sqrt{\xi}$ . The resulting annihilation rate ratio is derived as:
$\frac{\Gamma_A(\xi)}{\Gamma_A(1)} = \xi \frac{\tanh^2[\sqrt{\xi} t_\odot/\tau(1)]}{\tanh^2[t_\odot/\tau(1)]}$
This formulation allows the authors to map the departure from linear scaling based on the full-density equilibration timescale and the thermal fraction $\xi$ .
Constraints: The sample is filtered by experimental constraints, including the Higgs mass ($123-127$ GeV), flavor observables ( $b \to s\gamma$ , $B_s \to \mu^+\mu^-$ ), and direct detection limits (LZ graphical reference). The study explicitly notes that Sommerfeld corrections are omitted, which affects wino-like populations.

Key Contributions

Abundance–Timescale Map: The authors construct a general map relating the thermal abundance fraction ( $\xi$ ) and the full-density equilibration timescale ( $t_\odot/\tau$ ) to the suppression factor of the solar signal. This map demonstrates that the signal suppression is a function of both the relic fraction and the specific composition-dependent capture/annihilation rates.
Composition-Dependent Analysis: The study quantifies how different neutralino compositions respond to abundance rescaling. It highlights that wino-like models, which often dominate the selected sample, are particularly susceptible to moving away from equilibrium when $\xi$ is reduced, leading to strong signal suppression.
Benchmarking: The paper provides five representative benchmarks (B1–B5) illustrating distinct behaviors:
- Wino-like (B1, B2): Show significant deviation from linear scaling due to long equilibration times at reduced density.
- Bino-like (B3): Remains far from equilibrium regardless of density.
- Higgsino-like (B5): Exhibits nearly linear scaling because it is already in equilibrium at full density and remains so at thermal fractions.
IceCube Comparison: The integrated neutrino ( $\nu_\mu + \bar{\nu}_\mu$ ) and through-going muon fluxes are calculated and compared against published IceCube reference limits (both 3-year and 10-year analyses). The authors clarify that these comparisons serve as context rather than definitive exclusion statements due to differences in annihilation spectra and detector response assumptions.

Results

Non-Linear Scaling: For the selected $SU(5)$ sample, linear rescaling substantially overestimates the solar signal. The median ratio of the actual thermal signal to the linearly scaled signal is approximately 0.0522.
Equilibrium Status: While 19.9% of the baseline sample is near equilibrium ( $2\Gamma_A/C_\odot \ge 0.95$ ) at full density, this drops to 1.9% when scaled to thermal fractions.
Composition Effects:
- Wino-like: 92% of these models fall below a factor of 0.5 of the linear estimate; only 1.4% remain near equilibrium at thermal fractions.
- Bino-like: 61.7% fall below the 0.5 threshold; none are near equilibrium.
- Higgsino-like: 82.1% are near equilibrium at full density, and 57.1% remain near equilibrium at thermal fractions, exhibiting behavior closer to linear scaling.
Direct Detection vs. Solar Signals: The study confirms that the departure from linear solar scaling persists even in the subset of models that satisfy the graphical LZ direct detection limits. Unlike direct detection rates, which scale linearly with local density, solar signals do not.

Significance and Claims
The paper claims to provide a consistent framework for interpreting solar searches for subdominant neutralino dark matter. Its primary significance lies in correcting the common assumption that solar signals from underabundant DM scale linearly with the local density fraction. The authors demonstrate that:

Equilibrium at full density does not guarantee linear signal rescaling when the abundance is reduced.
The response is jointly determined by the neutralino composition, the relic fraction, and the equilibration timescale.
For many models, particularly wino-like candidates, the expected signal is significantly lower than linear extrapolations suggest, potentially rendering them undetectable by current IceCube limits even if they are viable DM candidates.

The authors maintain a modest tone, noting that their numerical predictions are conditional on the specific particle calculations used (e.g., omission of Sommerfeld corrections) and that model-specific exclusion requires mixed-channel spectra and detailed detector response modeling not fully implemented here. The work serves as a necessary step to refine the interpretation of IceCube solar neutrino data in the context of specific SUSY scenarios.

Technical Summary: Solar Neutrino Signals and Capture–Annihilation Equilibrium of SU(5) Neutralino Dark Matter

More like this