Gravitational scalar production with a generic reheating scenario

This paper investigates the gravitational production of decoupled scalars across various reheating scenarios, deriving constraints on their properties and demonstrating that while universal gravity effects do not necessarily invalidate non-thermal dark matter models for specific inflaton potentials, the reheating dynamics significantly influence relic abundance and impose stringent limits on scalar parameters.

The Gist

The Big Picture: The Universe's "Ghost" Particles

Imagine the early Universe as a giant, expanding balloon. Inside this balloon, there is a mysterious, invisible substance called Dark Matter. Scientists have a theory about how this Dark Matter was created: it didn't come from a hot soup of particles (like the rest of us), but rather "froze in" from nothing, slowly appearing over time.

However, there is a problem. Even if these particles don't talk to normal matter, gravity itself might be creating too many of them. If gravity creates too many, the Universe would be too heavy and would collapse or look completely different than it does today.

This paper asks: How much gravity-made Dark Matter is created, and does it ruin our theories? The answer depends entirely on how the Universe "cooled down" after its initial rapid expansion (a period called inflation).

---

The Cast of Characters

1. The Inflaton (The Inflator): A field that caused the Universe to expand incredibly fast. After it finished, it started vibrating like a plucked guitar string.
2. The Scalar (The Ghost): The candidate for Dark Matter. It is "decoupled," meaning it ignores almost everything else. It only interacts via gravity.
3. Reheating (The Cooling Phase): The period after inflation where the energy of the inflaton turns into the hot soup of particles we know today. This is the "reheating" phase.

---

The Two Ways Ghosts Are Born

The paper looks at two different ways these "ghost" particles are created by gravity:

1. The "Stochastic Noise" (During Inflation)

Imagine the inflaton field is a calm lake. During inflation, the rapid expansion of the Universe creates "waves" (quantum fluctuations) on the surface. Even if the ghost particles are heavy, these waves can push them into existence.

• The Analogy: Think of a snowstorm. If the wind (inflation) blows hard enough, it blows snowflakes (ghost particles) into the air. Once the wind stops, the snowflakes start to fall.
• The Catch: If the snow piles up too high, it crushes the house (over-produces Dark Matter). The paper calculates exactly how much snow piles up based on how the house cools down afterward.

2. The "Quantum Gravity Operators" (During Reheating)

Even if the ghost particles are invisible, the laws of quantum gravity suggest they might have tiny, hidden connections to the inflaton field.

• The Analogy: Imagine the inflaton is a drum being beaten. Even if the ghost is in a soundproof room, the vibrations of the drum might be strong enough to rattle the ghost's cage and shake a few out.
• The Catch: These "shakes" happen most violently right at the start of the drumming. The paper calculates how many ghosts get shaken out based on the rhythm of the drum.

---

The Critical Factor: The "Cooling Rate" (Reheating)

The most important discovery in this paper is that how fast the Universe cools down changes everything.

The authors tested different "cooling scenarios" (Reheating), which are like different ways to let a hot oven cool down:

• Scenario A: Instant Cooling (The Open Window)
  The oven cools instantly. The paper finds that in this case, the constraints on the ghost particles are very strict. There isn't much room for error.

• Scenario B: The "Slow Cool" (The Insulated Oven)
  The oven cools down slowly over a long time.

  • If the inflaton is a "Soft" oscillator (Power $k < 4$): Think of this as a gentle, slow vibration. If the cooling is slow, the "ghosts" created early on get diluted. It's like adding a lot of water to a cup of coffee; the coffee (ghosts) is still there, but it's weak.
    • Result: This is good news! It means we can have stronger interactions between the ghost and the normal world without breaking the Universe.
  • If the inflaton is a "Hard" oscillator (Power $k > 4$): Think of this as a violent, jerky vibration. If the cooling is slow, the "ghosts" actually get concentrated. It's like squeezing a sponge; the water (ghosts) gets pushed together.
    • Result: This is bad news. It creates too many ghosts, making the theory impossible unless the interactions are incredibly weak.

• Scenario C: Multi-Stage Cooling (The Staircase)
  What if the cooling happens in steps? First, it cools like a soft oscillator, then switches to a hard one?

  • The Finding: The paper shows that the effects factorize. This means the first stage creates the ghosts, and the subsequent stages just act as a "dilution" or "concentration" filter. You can calculate the total effect by multiplying the steps together.

---

The Bottom Line

The paper concludes that gravity is a double-edged sword.

1. For "Soft" Universes ($k < 4$): A long, slow cooling period acts as a safety valve. It dilutes the extra particles created by gravity, allowing for more flexible and interesting theories about Dark Matter.
2. For "Hard" Universes ($k > 4$): A long cooling period acts like a pressure cooker. It amplifies the particle production, making the theory very fragile and requiring extremely weak interactions to avoid disaster.

In simple terms: Whether the Universe is a "good" place for these invisible particles to exist depends entirely on the "music" the Universe played while it was cooling down. If the music was the right kind of slow and gentle, the ghosts stay hidden and manageable. If the music was too jerky or the cooling too slow in the wrong way, the ghosts take over.

Technical Summary

Technical Summary: Gravitational Scalar Production with a Generic Reheating Scenario

Problem Statement
The paper addresses a fundamental challenge in non-thermal dark matter (DM) scenarios, specifically the freeze-in mechanism. While freeze-in relies on the assumption of a negligible initial DM abundance and feeble interactions with the Standard Model (SM), gravitational particle production during the early Universe can generate significant relic abundances even for decoupled scalars. This gravitational production occurs via two primary channels: (1) stochastic fluctuations during inflation (Starobinsky-Yokoyama mechanism) and (2) production driven by Planck-suppressed operators induced by quantum gravity during the reheating phase. If these gravitational contributions are too large, they lead to the over-production of DM, potentially invalidating freeze-in models or requiring unnaturally small coupling constants. The authors investigate how the dynamics of the reheating phase—specifically the form of the inflaton potential and the duration of the reheating epoch—modulate these gravitational production rates.

Methodology
The authors perform a comprehensive analytical study of a decoupled scalar field $s$ with mass $m_s$ and self-interaction coupling $\lambda_s$. The analysis is divided into two main sections:

1. Production During Inflation:

   • The scalar field acquires a variance $\langle s^2 \rangle_{\text{end}}$ due to stochastic noise from horizon-crossing modes.
   • The evolution of the resulting condensate is tracked from the end of inflation through various reheating scenarios. The condensate initially oscillates around a quartic potential (radiation-like, $\rho \propto a^{-4}$) and later transitions to a quadratic potential (matter-like, $\rho \propto a^{-3}$) when the bare mass dominates.
   • The yield $Y = n_s/s_{SM}$ is calculated at the onset of the matter-dominated phase, accounting for entropy conservation and the scaling of the Hubble parameter $H$ and scale factor $a$.

2. Production During Reheating:

   • The authors consider Planck-suppressed operators arising from quantum gravity effects: $\mathcal{O}_1 = C_1 M_P^{-2} (\partial s)^2 \phi^2$ and $\mathcal{O}_2 = C_2 M_P^{-2} \phi^4 s^2$, where $\phi$ is the inflaton.
   • The inflaton is modeled as oscillating around a generic power-law potential $V_{\text{inf}} \propto \phi^k$.
   • Using the Boltzmann equation, the production rate is integrated over the reheating epoch. The analysis utilizes the Fourier decomposition of the oscillating inflaton field to compute interaction amplitudes.
   • The study extends from single-stage reheating to multi-stage scenarios where the inflaton potential changes form (e.g., from $k_1$ to $k_2$) during the reheating phase.

Key Scenarios Analyzed
The paper systematically evaluates five reheating scenarios:

1. Instantaneous Reheating: Immediate transition to radiation domination.
2. Matter-Dominated Reheating: Oscillation around a quadratic potential ($k=2$).
3. Generic Power-Law Reheating: Oscillation around $V \propto \phi^k$.
4. Two-Stage Reheating: A transition from potential $k_1$ to $k_2$.
5. Multi-Stage ($m$-stage) Reheating: A sequence of $m$ potential transitions.

Key Results

• Dependence on Reheating Dynamics ($k$):

  • $k < 4$: The reheating phase acts as a dilution mechanism. A prolonged reheating epoch (lower reheating temperature $T_{\text{reh}}$) significantly relaxes constraints on the scalar mass $m_s$ and self-interaction $\lambda_s$. This allows for viable freeze-in scenarios with larger couplings than in high-temperature scenarios.
  • $k = 4$: The constraints become independent of the reheating temperature, recovering the results of the instantaneous reheating scenario.
  • $k > 4$: The reheating phase enhances the relic abundance compared to the radiation-dominated case. Consequently, constraints on the scalar become more stringent as the reheating temperature decreases (or the reheating duration increases).

• Multi-Stage Reheating Factorization:

  • In multi-stage scenarios, the final relic abundance factorizes. The production is dominated by the earliest stage (immediately after inflation). Subsequent stages act purely as dilution or enhancement factors depending on their respective power-law indices $k_i$.
  • For inflationary production, the constraint coefficients depend only on the initial power $k_1$.
  • For reheating production via operators $\mathcal{O}_1$ and $\mathcal{O}_2$, the production is UV-dominated (peaking at the earliest times) for $k < 7$. Thus, the first stage determines the total yield, while subsequent stages modulate it.

• Constraints on Parameters:

  • Self-Interaction ($\lambda_s$): The paper derives upper bounds on $\lambda_s$ as a function of $m_s$ and $T_{\text{reh}}$. For $k < 4$, low $T_{\text{reh}}$ opens up parameter space for GeV–TeV mass scalars.
  • Wilson Coefficients ($C_1, C_2$): Constraints on the quantum gravity-induced operators are derived. For natural values $|C_1| \sim \mathcal{O}(0.01-1)$, the GeV–TeV mass range remains accessible. However, constraints on $|C_2|$ are generally more stringent, often requiring very low reheating temperatures to avoid over-production for keV–TeV scalars.

Significance and Claims
The paper claims to demonstrate that the predictivity of non-thermal dark matter scenarios is highly sensitive to the details of the reheating history.

• Dilution vs. Enhancement: The authors show that universal gravity effects do not necessarily spoil non-thermal DM scenarios with $k < 4$ and low reheating temperatures, as the extended reheating phase efficiently dilutes gravitationally produced relics. Conversely, for $k > 4$, the enhancement of relics during reheating imposes stringent constraints.
• Factorization: A key theoretical finding is that in multi-stage reheating, the effects of subsequent stages factorize into simple dilution/enhancement terms, allowing for a compact description of complex reheating histories.
• Viability of Low-T Reheating: The analysis supports the viability of low-reheating temperature regimes, showing they allow for larger DM-SM couplings (up to $\mathcal{O}(1)$) while avoiding thermalization, making these scenarios experimentally testable.

The authors conclude that incorporating realistic post-inflationary dynamics is essential when confronting feebly-interacting DM scenarios with cosmological constraints, as the reheating dynamics can either rescue or rule out specific model parameters.