Millicharged Particle Production During Late-Stage Stellar Evolution

This paper computes and provides semi-analytical fits for millicharged particle energy-loss rates in pre-supernova stars, identifying three distinct production regimes dominated by plasmon decay and electron-positron annihilation under high-temperature, low-plasma-frequency conditions.

The Gist

Imagine a star not just as a giant ball of fire, but as a bustling, super-hot factory. Inside this factory, atoms are smashing together, creating immense heat and pressure. Usually, the only things that can escape this factory are light (photons) and ghostly particles called neutrinos, which zip right through everything without stopping.

But what if there were secret escape artists hiding inside the star?

This paper is about a hypothetical type of particle called a Millicharged Particle (MCP). Think of these as "ghosts with a tiny bit of static electricity." They are so light and interact so weakly with normal matter that they can slip out of the star's core like a thief sneaking out of a vault, carrying away energy.

If these particles exist, they would act like a leaky radiator for the star. By stealing energy, they would cool the star down faster than expected, changing how the star ages and, eventually, how it dies (exploding as a supernova).

Here is a simple breakdown of what the scientists in this paper did:

1. The Setting: The Star's "Final Exam"

The authors focused on massive stars in their late stages of life, just before they explode.

• The Environment: The core is incredibly hot (like a million degrees) and dense.
• The Problem: Previous studies looked at how these particles escape from cooler stars (like our Sun) or very dense, cold stars (like white dwarfs). But no one had done the math for these specific "pre-explosion" stars, where the conditions are a unique mix of high heat and high density.

2. The Three Escape Routes

The paper calculates exactly how these particles escape the star's core. Depending on the particle's weight (mass) and the star's temperature, they use one of three "doors":

• Door A: The "Plasmon Decay" (The Pop)

  • Analogy: Imagine the star's core is filled with a fog of charged particles. Sometimes, a ripple in this fog (called a "plasmon") spontaneously pops and splits into two millicharged particles.
  • When it happens: This is the main escape route for very light particles. It's like a bubble bursting in a soda.

• Door B: The "Compton Scatter" (The Bump)

  • Analogy: Imagine a billiard ball (an electron) getting hit by a fast-moving cue ball (a photon). Instead of just bouncing, the collision is so energetic that it creates two new ghost particles that fly off.
  • When it happens: This happens when the particles are heavier and the star is hot, but not too hot. It's like a high-speed collision creating new debris.

• Door C: The "Pair Annihilation" (The Crash)

  • Analogy: In the hottest parts of the star, matter and antimatter (electrons and positrons) are constantly being created and destroyed. When they crash into each other, they usually vanish into energy. But sometimes, instead of vanishing, they turn into two millicharged particles.
  • When it happens: This is the dominant route when the star is extremely hot (hotter than the mass of an electron).

3. The "Recipe Book" for Astronomers

The most important part of this paper isn't just the theory; it's the math.
The authors didn't just say "particles escape." They wrote down precise formulas (like a recipe) that tell astronomers exactly how much energy is lost for every type of particle and every type of star condition.

• Why does this matter? Astronomers use computer simulations to predict how stars live and die. Before this paper, they didn't have the right "ingredients" to add millicharged particles into their simulations. Now, they can plug these new formulas into their code.

4. The Big Picture: Catching the Ghosts

The authors found that for certain types of heavy millicharged particles, the "leak" they cause in massive stars could be huge—bigger than the leak caused by neutrinos.

• The Detective Work: If we observe massive stars and they seem to be cooling down or evolving differently than our current models predict, it might be a clue that these "ghosts" are real.
• The New Frontier: This paper opens a new window to test physics. We might not be able to see these particles in a lab on Earth, but the stars themselves are giant laboratories testing the limits of our universe.

Summary

Think of this paper as a manual for a new type of star thermostat. The authors figured out how "millicharged ghosts" steal heat from dying stars. By providing the exact math for this heat loss, they give astronomers a new tool to either find these particles or prove they don't exist, helping us understand the final moments of the universe's most massive stars.

Technical Summary

1. Problem Statement

Millicharged particles (MCPs), hypothetical particles with mass $m_\chi$ and a tiny electric charge $q e$ (where $q \ll 1$), are predicted by various Beyond the Standard Model (BSM) scenarios (e.g., kinetic mixing with a dark photon). Stars act as natural laboratories for detecting such particles because their cores produce them via thermal interactions, leading to energy loss that can alter stellar evolution.

While MCP production rates have been calculated for various stellar environments (e.g., the Sun, red giants, white dwarfs), there was a significant gap in the literature regarding late-stage massive stars (pre-supernova objects). In these environments, temperatures reach $T \sim 10\text{--}100$ keV and densities are high, but the plasma frequency $\omega_{pl}$ is typically much smaller than the temperature ($\omega_{pl} \ll T$). Previous calculations often relied on approximations valid for degenerate or lower-temperature regimes, or failed to correctly bridge the transition between non-relativistic and ultra-relativistic plasma conditions. The authors aim to compute precise energy-loss rates for MCPs in these specific pre-supernova conditions to constrain MCP parameters.

2. Methodology

The authors perform a rigorous quantum field theoretic calculation of MCP production rates in a hot, non-degenerate, non-relativistic to ultra-relativistic plasma. They focus on three dominant production mechanisms (labeled A, C, and D), analogous to neutrino production processes:

1. Plasmon Decay (D-process): $\gamma_{L,T} \to \chi \bar{\chi}$
   • Dominant when $m_\chi < \omega_{pl}/2$.
   • The authors derive the emissivity for both longitudinal and transverse plasmon modes, accounting for medium-induced dispersion relations and wavefunction renormalization.
2. Compton-like Scattering (C-process): $e^- + \gamma \to e^- + \chi \bar{\chi}$
   • Dominant when $m_\chi > \omega_{pl}/2$ and $T \lesssim 0.5$ MeV.
   • The authors calculate the cross-section for off-shell photon decay into MCP pairs. They explicitly address the treatment of longitudinal plasmons, correcting previous literature (Refs. [13, 14]) which used an incorrect self-energy expression for longitudinal modes in the off-shell regime.
   • They derive separate limits for the Non-Relativistic (NR) regime ($T, m_\chi \ll m_e$) and the Ultra-Relativistic (UR) regime ($T \gg m_e$).
3. Pair Annihilation (A-process): $e^+ + e^- \to \chi \bar{\chi}$
   • Dominant when $m_\chi > \omega_{pl}/2$ and $T \gtrsim m_e$ (where positron abundance is significant).
   • Calculated using the standard QED cross-section for $e^+e^- \to \mu^+\mu^-$, rescaled by the MCP charge factor $q^2$.

Key Technical Steps:

• Plasma Conditions: They model the stellar core as a non-degenerate plasma where $\omega_{pl} \ll T$. They verify that for the relevant temperatures, longitudinal plasmons decouple from incoherent processes in the Compton regime, allowing them to focus on transverse contributions.
• Derivation of Emissivity: They compute the energy-momentum loss rate ($dE/dVdt$) by integrating matrix elements over phase space, utilizing relativistic invariance to simplify integrals.
• Semi-Analytical Fits: To facilitate implementation in stellar evolution codes, the authors derive semi-analytical fitting functions for the emissivity in the NR and UR limits and construct an interpolating function valid across the entire temperature range ($T \ll m_e$ to $T \gg m_e$).

3. Key Contributions

• First Calculation for Pre-Supernova Stars: This is the first dedicated computation of MCP energy-loss rates specifically for the late evolutionary stages of massive stars (post-main sequence, pre-core collapse), where $T \sim 10\text{--}100$ keV and $\omega_{pl} \ll T$.
• Correction of Previous Literature: The authors identify and correct an error in previous works regarding the self-energy of longitudinal plasmons in the Compton scattering regime. They demonstrate that for off-shell photons, the longitudinal contribution behaves differently than previously assumed, though the transverse contribution remains dominant.
• Unified Interpolation Formula: They provide a robust, semi-analytical formula (Eq. 69) that smoothly interpolates between the non-relativistic and ultra-relativistic limits for Compton production, covering the transition region around $T \sim m_e$.
• Regime Mapping: They map out the dominant production mechanism as a function of temperature ($T$), density ($\rho$), and MCP mass ($m_\chi$), identifying three distinct regimes:
  1. Low Mass ($m_\chi < \omega_{pl}/2$): Plasmon decay dominates.
  2. Intermediate Mass/Temp ($m_\chi > \omega_{pl}/2, T \lesssim 0.5$ MeV): Compton scattering dominates.
  3. High Temp ($T \gtrsim m_e$): Pair annihilation becomes competitive or dominant.

4. Results

• Emissivity Fits: The paper provides explicit fitting functions for the emissivity $\epsilon$ (energy loss per unit volume per unit time) for all three processes. These fits are accurate to within $\sim 10\%$ across the relevant parameter space.
  • Plasmon Decay: Dominated by transverse modes; emissivity scales roughly as $T^3$ in the $\omega_{pl} \ll T$ limit.
  • Compton Scattering: The derived emissivity corrects previous estimates, particularly in the longitudinal sector, and provides a continuous function for $T$ ranging from keV to MeV.
  • Pair Annihilation: Emissivity scales as $T^5$ and is significant only when $T \gtrsim m_e$.
• Dominance Regions:
  • For light MCPs ($m_\chi = 1$ keV), plasmon decay is the dominant channel in high-density regions ($\rho \gtrsim 5 \times 10^3$ g cm$^{-3}$). Compton scattering takes over only at lower densities where plasmon decay is kinematically forbidden.
  • For heavier MCPs ($m_\chi = 100$ keV), Compton scattering dominates at lower temperatures, while pair annihilation dominates at high temperatures ($T \gtrsim m_e$).
• Comparison with Neutrinos: The authors compare the MCP energy loss ($Q_\chi$) to standard neutrino energy loss ($Q_\nu$). They find that for specific parameter spaces (e.g., $m_\chi \sim 100$ keV and $q \sim 10^{-10}$), MCP emission can exceed neutrino emission during the Helium-burning phase of massive stars.

5. Significance

• New Astrophysical Constraints: The results open a new window for constraining MCPs in a previously unexplored mass range ($10\text{--}100$ keV) and charge range ($q \sim 10^{-10}$). These parameters were not previously constrained by red giant or white dwarf observations.
• Impact on Stellar Evolution: If MCPs exist with these parameters, their enhanced energy loss could significantly alter the evolution of massive stars, potentially affecting the black hole mass gap or the timing of core collapse.
• Tool for Future Research: The provided semi-analytical fits are designed for direct implementation in stellar evolution codes. The authors indicate that forthcoming papers will utilize these rates to simulate the specific effects of MCPs on the late stages of massive star evolution, potentially offering new observational signatures for astrophysicists.

In summary, this paper provides the necessary theoretical foundation and computational tools to rigorously test the existence of millicharged particles using the extreme conditions found in the cores of massive stars before they explode as supernovae.