🔬 physics

Superheavy dark matter and stepped dark radiation for Hubble Tension

This paper proposes a superheavy nonthermal dark matter framework where gravity-induced collapse into quantum-gravitational bound states converts an overabundance of dark matter into dark radiation, thereby increasing $N_{\text{eff}}$ to resolve the Hubble tension while remaining consistent with cosmological observations.

Original authors: Zhijie Xu

Published 2026-07-21

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Zhijie Xu

Original paper licensed under CC BY 4.0 (https://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

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to measure exactly how fast this balloon is inflating. They have two main ways to do this: looking at the "baby pictures" of the universe (the Cosmic Microwave Background) and measuring the speed of nearby exploding stars (supernovae). The problem is, these two methods give different answers. It's like asking two people to measure the same room with a tape measure, but one says it's 10 feet wide and the other says 12. This disagreement is called the "Hubble Tension," and it's one of the biggest headaches in modern physics because it suggests our understanding of the universe's history might be missing a crucial piece of the puzzle.

To solve this, scientists often look for "dark" ingredients—things we can't see but know must be there because of their gravity. Two popular suspects are "Dark Matter" (the invisible glue holding galaxies together) and "Dark Radiation" (invisible, fast-moving particles that act like extra heat). The paper you're about to read explores a wild new idea: what if the invisible glue is actually made of super-heavy, non-thermal particles that behave very differently from the standard suspects we usually imagine?

The Case of the Super-Heavy Ghosts

In this paper, author Zhijie Xu proposes a solution to the Hubble Tension that involves a very strange type of dark matter. Let's call these particles "Super-Heavy Ghosts." Unlike the standard dark matter particles (known as WIMPs) that are thought to be light and thermal (like hot soup cooling down), these ghosts are incredibly massive—about $10^{12}$ GeV. To put that in perspective, if a standard atom were the size of a house, these particles would be the size of a mountain.

Here is the story of how these ghosts save the day, told in three acts.

Act 1: The Freezing Out

Usually, scientists think dark matter was created in a hot, dense soup and then "froze out" as the universe cooled, leaving a steady amount behind. But Xu suggests these Super-Heavy Ghosts were created differently. They were born non-thermal (not from a hot soup) and were initially produced in massive excess—way more than we see today.

As the universe expanded, these heavy ghosts cooled down incredibly fast. They became "ultra-cold," moving so slowly that their thermal motion (wiggling around due to heat) became negligible. In fact, they became so cold that their movement was dominated not by heat, but by the gentle tug of gravity from large structures in the universe. This is the first key step: the ghosts go from being a hot, chaotic crowd to a frozen, silent congregation.

Act 2: The Great Collapse and Conversion

This is where things get weird. Because there were so many of these ghosts initially, and they were so cold, they started to clump together under their own gravity. Imagine a crowd of people so dense and still that they collapse into a single, tight huddle.

In this paper, the author suggests that these clumps formed tiny, quantum-gravitational bound states. Think of it like two magnets snapping together, but on a scale so small ( $10^{-13}$ meters) that quantum mechanics takes over. Inside these tiny, super-dense "droplets," the particles interacted in a special way. A tiny, hidden connection (called a "portal") allowed these bound pairs to instantly transform into something else: Dark Radiation.

This process is called a "cold freeze-out." It's like a magical factory where the heavy ghosts walk in, get squeezed into a tiny box, and pop out as invisible, fast-moving radiation. The paper calculates that this factory was incredibly efficient. For every billion ghosts that started the process, only one survived as a ghost. The other 999,999,999 were converted into Dark Radiation.

Act 3: Solving the Puzzle

So, why does this matter? The universe is filled with this newly created Dark Radiation. This extra radiation changes the expansion history of the early universe. Specifically, it shrinks the "sound horizon"—a cosmic ruler used to measure distances in the early universe.

By shrinking this ruler, the model allows the "baby picture" of the universe (the CMB) to predict a faster expansion rate today. This brings the two conflicting measurements of the Hubble constant into agreement. The paper suggests that this conversion process produces just the right amount of extra radiation, increasing the effective number of neutrino species by about $\Delta N_{eff} \simeq 0.4$ . This number is in the "sweet spot" that could resolve the tension without breaking other rules of physics.

The Blueprint: How It All Fits Together

The paper doesn't just guess; it builds a mathematical framework to support this story.

The Mass: The author identifies a unique mass scale of $10^{12}$ GeV. This isn't a random guess; it's the specific mass where the "free streaming" (how far a particle can travel) equals the mass of the particle itself. It's the boundary where Newtonian quantum gravity becomes the main player.
The Timing: This conversion happens very early, around $t_X \simeq 10^{-6}$ seconds after the Big Bang. This is before the universe gets hot enough for nuclear fusion (Big Bang Nucleosynthesis), so it doesn't mess up the formation of elements like hydrogen and helium.
The Physics: The model uses a "stepped" dark radiation framework. It suggests that the dark radiation doesn't just appear all at once; it has a "step" in its energy density. The Super-Heavy Ghosts provide the initial overabundance, and their conversion creates the step.
The Drag: The model also includes a "drag" effect. The dark matter and dark radiation interact with each other (like two fluids rubbing together) before the step happens. This interaction helps smooth out the clumping of matter on small scales, which might also help solve a different problem called the " $S_8$ tension" (where the universe seems less clumpy than we expect).

The Fine Print: What We Don't Know Yet

While the story is compelling, the paper is careful to note that this is a theoretical proposal, not a confirmed fact.

It's a Suggestion: The authors suggest this mechanism could work. They have solved the equations (the Boltzmann equation) to show that the numbers add up, but they haven't observed these particles yet.
The Assumptions: The model relies on a few big assumptions, such as the existence of a specific type of supersymmetry (a theory about extra dimensions and partner particles) and a "no-scale" supergravity setup. These are complex theories that are still being tested by physicists.
The Limits: The paper admits that the exact details of how the particles collapse and convert are simplified. Real simulations of this "Jeans collapse" (the gravitational clumping) would be needed to confirm the efficiency of the conversion.

The Takeaway

In simple terms, this paper offers a creative new chapter for the story of the universe. It proposes that the universe was once filled with a massive, overabundant population of super-heavy dark matter ghosts. These ghosts froze, clumped together, and then magically converted into invisible radiation. This extra radiation tweaked the universe's expansion speed just enough to make our measurements of the Hubble constant agree.

It's a bold idea that links the microscopic world of quantum gravity with the macroscopic mystery of the universe's expansion. While it remains a hypothesis waiting for more data and simulations, it provides a fresh, mathematically consistent path to solving one of cosmology's most persistent riddles.

Technical Summary: Superheavy Dark Matter and Stepped Dark Radiation for Hubble Tension

Problem Statement
The paper addresses the persistent $\sim 5\sigma$ Hubble tension, the discrepancy between the Hubble constant ( $H_0$ ) inferred from Cosmic Microwave Background (CMB) data within the $\Lambda$ CDM model and local measurements using calibrated supernovae. Resolutions that modify the pre-recombination expansion history to reduce the sound horizon ( $r_s$ ) are particularly compelling. While Early Dark Energy (EDE) and interacting Dark Radiation (DR) are proposed solutions, they face tight constraints from CMB phase shifts, damping tails, and large-scale structure (specifically the $S_8$ tension). The paper seeks a mechanism that generates a specific amount of interacting dark radiation ( $\Delta N_{\text{eff}} \approx 0.4$ ) with a "step-like" injection history, consistent with both CMB and structure data, while linking this macroscopic cosmological effect to microscopic dark matter (DM) physics.

Methodology and Framework
The authors propose a framework based on superheavy nonthermal dark matter ( $X$ ) with a mass scale $m_X \sim 10^{12}$ GeV. Unlike thermal WIMPs, the abundance of $X$ is not set by chemical freeze-out but by a subsequent "cold" depletion stage. The methodology proceeds through three distinct physical stages:

Free Streaming and Critical Mass (Step 1): The authors perform a free-streaming analysis to identify a critical particle mass where the free-streaming mass ( $M_{fs}$ ) equals the individual particle mass ( $m_X$ ). For nonthermal relics, $M_{fs} \propto m_X^{-3}$ . The crossing condition $M_{fs}(m_X) = m_X$ yields a unique critical mass $m_{Xc} \approx 10^{12}$ GeV. At this scale, Newtonian quantum gravity becomes relevant for single-particle binding, delineating the regime where the smallest possible dark matter structures form.
Gravity-Induced Velocity (Step 2): The paper analyzes the evolution of DM peculiar velocity, decomposing it into a transverse thermal component ( $V_\perp$ ) and a longitudinal gravity-induced component ( $V_\parallel$ ). In the ultra-cold regime of superheavy $X$ , $V_\perp$ redshifts rapidly to negligible values ( $\sim 10^{-3}$ m/s), while $V_\parallel$ remains "frozen" at a constant scale of $\sim$ km/s during radiation domination, sourced by large-scale density perturbations. This velocity hierarchy enables the particles to undergo gravitational collapse despite their extreme coldness.
Cold Freeze-Out and Bound-State Conversion (Step 3): When the thermal speed drops below $V_\parallel$ , the overproduced $X$ population (with initial density $\rho_{ini} \gg \rho_{\infty}$ ) becomes comparable to the radiation density. This triggers gravitational collapse into compact, quantum-gravitational bound states where quantum pressure balances gravity.
- Microphysics: The bound state size is determined by virial equilibrium and angular momentum quantization ( $m_X v_X r_X = \hbar$ ), yielding a characteristic radius $r_X \sim 10^{-13}$ m and an internal orbital velocity $v_X \sim 10^{-7}$ m/s.
- Depletion Mechanism: A Planck-suppressed "portal" coupling ( $y_K \sim 10^{-14}$ ) induces the conversion of these bound states into dark radiation. This process is modeled as an effective annihilation channel $X + X \to \text{Dark Radiation}$ with an effective cross-section $\langle \sigma v \rangle \approx 10^{-21} \text{ m}^3 \text{s}^{-1}$ .
- Boltzmann Solution: Solving the non-equilibrium Boltzmann equation reveals that the depletion is highly efficient, reducing the initial number density $n_i$ to a relic density $n_\infty = \gamma n_i$ with a survival fraction $\gamma \sim 10^{-9}$ . The converted energy thermalizes into a dark radiation bath.

Model Realization (WZDR+)
The cosmological scenario is embedded in a minimal dark sector based on the Stepped Wess–Zumino Dark Radiation (WZDR) framework:

Dark Sector Content: A chiral superfield $S = (\phi, \psi)$ (scalar $\phi$ , fermion $\psi$ ) and a dark matter superfield $X = (\tilde{\chi}, \chi)$ .
Interactions:
- A renormalizable Yukawa coupling $\lambda_X \phi X X$ ( $\lambda_X \sim 10^{-4}$ ) mediates DM-DR drag, active before the step and shutting off when $\phi$ becomes non-relativistic. This addresses the $S_8$ tension by suppressing power on 10–20 Mpc scales.
- A Planck-suppressed spurion portal $(Z^2/M_{Pl}^2) S X^2$ induces the tiny coupling $y_K$ governing the bound-state conversion.
Supersymmetry Breaking: The model utilizes no-scale supergravity with a high SUSY breaking scale $\sqrt{F} \sim 10^{12}$ GeV. This scale fixes the DM mass $m_X$ and, due to the cancellation of tree-level gravity-mediated soft masses, naturally generates an eV-scale mass for the scalar $\phi$ , enabling the "step" in DR energy density at redshift $z_t \approx 10^4$ .

Key Results

Mass Scale: The framework predicts a unique superheavy DM mass $m_X \approx 10^{12}$ GeV, derived independently from both free-streaming analysis and the cold freeze-out dynamics.
Relic Abundance: The "cold" freeze-out mechanism converts $\sim 99.9999999\%$ of the initial overabundance into dark radiation, leaving a relic fraction $\gamma \sim 10^{-9}$ .
Dark Radiation: The process generates an effective number of relativistic species $\Delta N_{\text{eff}} \approx 0.4$ . This value is sufficient to reduce the sound horizon and raise the CMB-inferred $H_0$ , potentially resolving the Hubble tension.
Timing: The conversion occurs at $t_X \sim 10^{-6}$ s (temperature $T \sim 1$ GeV), well before Big Bang Nucleosynthesis (BBN). The pre-step radiation density $N_{\text{UV}} \approx 0.4$ satisfies BBN constraints, while the post-step density $N_{\text{IR}} \approx 0.5$ (after entropy injection from $\phi$ decay) aligns with CMB constraints.

Significance and Claims
The paper claims to provide a candidate mechanism that links macroscopic cold freeze-out to microscopic couplings, offering a unified solution to the Hubble and $S_8$ tensions. By tying the origin and timing of dark radiation to the microphysics of superheavy nonthermal DM, the model avoids the fine-tuning issues of ad-hoc EDE models. The authors emphasize that their framework:

Alleviates Hubble Tension: Through the generation of $\Delta N_{\text{eff}} \approx 0.4$ .
Mitigates $S_8$ Tension: Via DM-DR drag interactions that suppress clustering on small scales.
Respects Constraints: The "step" nature of the DR injection differentially affects low- and high- $\ell$ CMB multipoles, softening phase-shift constraints, while the timing ensures compliance with BBN.

Limitations and Future Work
The authors acknowledge several limitations and areas for future development:

Numerical Uncertainties: The effective cross-section relies on idealized assumptions of complete virialization and prompt conversion; finite- $N$ and incomplete-virialization corrections require dedicated simulations.
Cosmological Likelihood: While the predicted $\Delta N_{\text{eff}}$ falls within favored ranges, a dedicated CMB+BAO+Large Scale Structure (LSS) likelihood analysis has not yet been performed to confirm full consistency with precision data.
UV Completion: The no-scale supergravity embedding provides a parametric argument for the mass scales and couplings, but a complete treatment involving explicit Kähler potentials, moduli stabilization, and radiative stability against loop corrections remains to be developed.