Tachyonic gravitational dark matter production after inflation

This paper proposes a novel non-thermal dark matter production mechanism driven by curvature-induced tachyonic instabilities of a spectator scalar field coupled to the Gauss-Bonnet invariant, which triggers explosive particle creation after inflation and successfully reproduces the observed dark matter relic density across a wide range of masses and inflationary scales.

The Gist

The Big Picture: A Gravity-Only Factory

Imagine the universe right after the Big Bang. For a long time, scientists have been trying to figure out what Dark Matter is. We know it's there because it holds galaxies together with its gravity, but we can't see it, and we haven't been able to catch it in a lab.

Most theories suggest Dark Matter is a particle that interacts with normal matter (like us) in some way. But this paper proposes a different idea: What if Dark Matter was created entirely by gravity, with no help from other forces?

The authors suggest a mechanism where the changing shape of space and time itself acts like a factory, manufacturing Dark Matter particles out of nothing but the energy of the universe's expansion.

The Setup: The "Spectator" and the "Inflaton"

To understand the story, we need two characters:

1. The Inflaton: This is the field that drove the rapid expansion of the universe (inflation) in the very beginning. Think of it as the engine that pushed the universe apart.
2. The Spectator (Dark Matter): This is a field that just sat there watching. It didn't drive the expansion, and it didn't talk to normal matter. It was a "spectator."

In standard physics, a spectator field usually just sits quietly. But the authors found a way to make it wake up and explode into existence.

The Trigger: The "Tachyonic" Switch

The key to the story is a moment right after inflation stopped. The universe was expanding rapidly, then suddenly slowed down to enter the "Radiation Domination" era (a hot soup of particles).

The authors propose that this sudden change in the universe's speed caused a curvature-induced instability.

The Analogy: The Tightrope Walker
Imagine the Dark Matter field is a tightrope walker balancing on a wire.

• During Inflation: The wire is perfectly straight and stable. The walker stands still at the center.
• The Switch: When inflation stops and the universe changes its expansion rate, the wire suddenly flips upside down. The center of the wire becomes a peak, and the sides become valleys.
• The Tachyonic Instability: In physics, when a field's "mass" becomes negative (like the walker falling off the peak), it's called "tachyonic." The walker doesn't just fall; they fall explosively fast.

Because the wire flipped, the spectator field couldn't stay at zero. It had to roll down into the valleys. As it rolled, it didn't just move; it amplified. Tiny, invisible quantum fluctuations (ripples in the field) were stretched and blown up into massive waves.

The Process: From Ripples to Particles

This "explosive rolling" created a huge amount of energy in the Dark Matter field.

• The Amplification: Just like shaking a rug creates a big wave, the changing curvature of space shook the Dark Matter field, turning tiny quantum whispers into a loud roar.
• The Result: These amplified waves eventually settled down and behaved like a swarm of particles. These particles are our Dark Matter.

The paper uses a specific mathematical tool called the Gauss-Bonnet invariant to describe this. Think of this as a specific type of "curvature sensor" built into the fabric of space. It's special because it only reacts to the shape of space, not to other forces, ensuring that this Dark Matter is truly "gravitational" and doesn't interact with light or normal matter.

The Simulation: Watching the Movie

The authors didn't just guess this would work; they ran complex computer simulations (using a "3+1 classical lattice," which is like a 3D grid of space evolving over time).

They watched the simulation play out:

1. The Flip: The universe transitions from inflation to radiation.
2. The Explosion: The Dark Matter field goes unstable. The energy density spikes.
3. The Cooling: The field stops growing and starts oscillating (wiggling back and forth).
4. The Transformation: At first, these wiggles act like light (radiation), moving fast. But as the universe expands and cools, they slow down and start acting like heavy, slow-moving matter (dust).

This transition from "fast light" to "slow matter" is crucial. It explains why we have the right amount of Dark Matter today.

The Conclusion: A Robust Recipe

The paper concludes that this mechanism is very robust. It works across a wide range of masses and energy scales.

• No Fine-Tuning: You don't need to set the dials perfectly. The mechanism works naturally because of how the universe expands.
• Pure Gravity: It requires no new forces or interactions with the Standard Model. It is a "pure gravity" solution.
• The Fit: They found a simple mathematical formula (a "fitting function") that predicts exactly how much Dark Matter would be created based on the mass of the particle and the energy of the early universe.

In summary: The universe's expansion history acted like a switch that flipped the stability of a hidden field. This caused the field to roll down a hill, amplifying tiny ripples into a sea of Dark Matter particles, all driven solely by the curvature of spacetime.

Technical Summary

Technical Summary: Tachyonic Gravitational Dark Matter Production after Inflation

Problem Statement
The particle nature of dark matter (DM) remains unknown despite extensive experimental searches for Weakly Interacting Massive Particles (WIMPs) and axions. This work addresses the problem of generating a non-thermal DM relic abundance through purely gravitational mechanisms, avoiding direct couplings to the Standard Model (SM). Specifically, the authors investigate whether the rapid reorganization of spacetime curvature at the end of inflation can trigger a tachyonic instability in a spectator scalar field, leading to explosive particle production without requiring fine-tuned initial conditions or non-minimal couplings to the Ricci scalar alone.

Methodology
The authors employ an Effective Field Theory (EFT) framework for gravity coupled to a real scalar spectator field $\chi$. The methodology combines analytical estimates with 3+1 classical lattice simulations.

1. Theoretical Framework:

   • The spectator field $\chi$ is coupled quadratically to spacetime curvature invariants via dimension-six operators suppressed by a scale $\Lambda$. The action includes the Gauss-Bonnet topological invariant ($G = R^2 - 4R_{\mu\nu}R^{\mu\nu} + R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma}$) as a benchmark, chosen because it preserves second-order equations of motion and avoids introducing ghost degrees of freedom.
   • The effective mass squared of the spectator field is derived as $M^2_{\text{eff}} = M^2 + \xi R + \frac{2}{\Lambda^2}(\alpha R^2 + \beta R_{\mu\nu}R^{\mu\nu} + \gamma R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma})$.
   • The authors analyze the transition from inflation ($w \approx -1$) to radiation domination (RD, $w = 1/3$). In the Gauss-Bonnet case, the curvature contribution to the mass is positive during inflation but becomes negative during RD, triggering a spontaneous symmetry breaking and a transient tachyonic instability.

2. Analytical Analysis:

   • The system is first treated in a homogeneous approximation to identify the timescale of symmetry restoration.
   • Quantum fluctuations are then quantized. The authors utilize the "uniform approximation" (mapping the mode equation to an Airy differential equation) to solve for the evolution of Fourier modes across the turning point where the effective mass changes sign.
   • This analysis identifies an instability band in momentum space, showing that infrared (IR) modes are exponentially amplified while ultraviolet (UV) modes remain unaffected. The authors demonstrate that these amplified fluctuations undergo a quantum-to-classical transition, justifying the use of classical lattice simulations.

3. Numerical Simulations:

   • To capture the full non-linear, stochastic evolution, the authors perform 3+1 dimensional classical lattice simulations using the CosmoLattice library on a fixed FLRW background.
   • They scan the parameter space defined by the DM bare mass ($M$), the EFT cutoff ($\Lambda$), and the Hubble scale at the end of inflation ($H_*$).
   • Simulations track the redistribution of energy among kinetic, gradient, interaction, and potential terms, and monitor the evolution of the equation of state from radiation-like ($w \approx 1/3$) to matter-like ($w \approx 0$).

Key Contributions and Results

• Mechanism Identification: The paper establishes that the transition to radiation domination naturally induces a sign flip in the effective mass of a spectator field coupled to the Gauss-Bonnet invariant. This triggers a short-lived phase of tachyonic instability and spontaneous symmetry breaking, leading to efficient, non-thermal particle production.
• IR Dominance and Classicalization: Analytical and numerical results confirm that the instability is dominated by IR modes. These modes grow exponentially and rapidly classicalize, allowing the system to be described as a classical stochastic field.
• Relic Density Calculation:
  • The authors derive a parametric estimate for the present-day DM abundance $\Omega_{\text{DM}}^0$, showing a strong dependence on the curvature-induced instability parameters.
  • Lattice simulations (approx. 200 runs) validate the analytical scaling and provide a precise fitting function:
    $$\Omega_{\text{DM}}^0(H_*, M, \Lambda) \sim 10^{24.31} \left(\frac{H_*}{M_P}\right)^{2.50} \left(\frac{M}{H_*}\right)^{0.97} \nu^{2.86} \exp[0.38 \nu]$$
    where $\nu = \chi_*/H_*$ and $\chi_* = \sqrt{48}H_*^2/\Lambda$.
  • The mechanism successfully reproduces the observed DM relic density ($\Omega_{\text{DM}}^0 h^2 \approx 0.12$) over a wide range of masses ($M \sim 10 \text{ GeV} - 10^5 \text{ GeV}$) and inflationary scales ($H_* \sim 10^{10} - 10^{13} \text{ GeV}$).
• Non-Thermal Nature: The produced DM is inherently non-thermal. The momentum distribution is set by the tachyonic amplification and preserved via redshifting, rather than thermalization. This results in a residual velocity dispersion that could mimic warm DM behavior depending on the parameter space, though the model generally approaches cold DM behavior if the bare mass becomes dominant early enough.

Significance and Claims
The paper claims to provide a "robust and predictive mechanism" for gravitational DM production that does not rely on direct couplings to the visible sector or fine-tuned cancellations. By focusing on the Gauss-Bonnet invariant, the authors isolate the curvature-driven instability as a property of the spectator sector on a standard Einsteinian background, avoiding the pathologies of modified gravity theories.

The significance of the work lies in:

1. Theoretical Control: Demonstrating that curvature-induced symmetry breaking is a generic feature of quadratic curvature couplings in EFTs, specifically realized here via the Gauss-Bonnet term which maintains theoretical consistency (no ghosts).
2. Viability: Showing that this mechanism can account for the observed DM abundance across a broad, well-motivated region of parameter space, offering an alternative to WIMP paradigms.
3. Practical Application: Providing a lattice-independent fitting function (Eq. 5.1) that allows for the application of these results to broader phenomenological studies without requiring new simulations.
4. Future Prospects: The authors note that the violent, out-of-equilibrium dynamics may source stochastic gravitational wave backgrounds, offering a potential observational signature for future high-frequency detectors, though a quantitative prediction of this spectrum is left for future work.

The authors remain modest regarding the immediate experimental verification, noting that the mechanism generates non-thermal DM with specific transfer functions that require further study against small-scale structure constraints (e.g., Lyman-$\alpha$ forest), which is identified as a direction for future research.