Supernova production of axion-like particles coupling to electrons, reloaded

This paper revisits the production of axion-like particles coupled to electrons in supernovae by incorporating previously neglected processes and thermal corrections to derive analytical emissivity rates, ultimately updating astrophysical constraints from SN 1987A across various coupling regimes.

The Gist

Imagine a supernova not just as a giant explosion, but as a cosmic pressure cooker. Inside this cooker, the matter is so hot and dense that it behaves like a super-charged fluid. For decades, physicists have been trying to figure out what happens to the energy inside this cooker. They suspect that invisible, ghostly particles called Axion-Like Particles (ALPs) might be sneaking out, stealing energy and cooling the star down faster than expected.

This paper is a "reloaded" version of a previous study. The authors, Damiano Fiorillo, Tetyana Pitik, and Edoardo Vitagliano, went back to the kitchen to check the recipe. They found that the old recipe was missing a few key ingredients and had some incorrect measurements. Here is what they discovered, explained simply:

1. The "Ghost" Particles and the Electron Dance

ALPs are hypothetical particles that interact very weakly with normal matter. In this study, the authors focus on how ALPs interact with electrons.

Think of the electrons in the supernova as a crowded dance floor.

• The Old View: Scientists previously thought the main way ALPs were produced was like a "bremsstrahlung" (braking radiation) event: an electron bumps into a heavy nucleus (like a dancer bumping into a wall) and emits an ALP.
• The New Discovery: The authors realized there is a much more common dance move they ignored: Semi-Compton Scattering. Imagine an electron (dancer) colliding with a photon (a flash of light) and, in the chaos, spitting out an ALP.
  • The Analogy: It's like realizing that while dancers bumping into walls do create noise, the dancers bumping into each other while spinning in the light is actually creating way more noise. This "Semi-Compton" process turns out to be the dominant way ALPs are made in these extreme environments, surpassing the old "wall-bumping" method.

2. Fixing the "Thermal Mass" Confusion

In a super-hot plasma, particles don't act like they do in a vacuum. They gain a "thermal mass" (they get "heavier" because of the heat).

• The Mistake: Previous calculations treated these hot electrons as if they were just heavy versions of cold electrons.
• The Correction: The authors explain that the math for these hot electrons is more subtle. You can't just slap a "heavy" label on them; you have to change how they move and interact fundamentally. They fixed the math to ensure they aren't overestimating or underestimating how many ALPs are being produced.

3. The "Fireball" and the "Trapping" Effect

The paper also looks at what happens if the ALPs interact too strongly.

• The Trap: If the coupling is too strong, the ALPs get stuck inside the star. They bounce around like a pinball in a machine, unable to escape. This is called the "trapping regime."
• The Fireball: If they do escape but are still interacting heavily, they might create a "fireball" of new particles (electrons and positrons) that glow in X-rays.
• The New Rule: The authors developed a new, more precise way to calculate how much energy gets trapped versus how much escapes. This is crucial because if we get the math wrong, we might think a star is cooling too fast when it's actually just holding onto its heat.

4. The "Decay" Detective Work

Once ALPs escape the supernova, they might decay (break apart) into other particles. The authors looked at three ways this happens:

1. Into Positrons: These create a specific signal (511 keV line) that we can look for in the galaxy.
2. Into Two Photons: A direct flash of light.
3. Into an Electron, Positron, and Photon: This is the new big find. The authors realized that this specific three-part decay produces a lot of gamma rays that previous studies missed.
   • The Impact: This new channel is actually the strongest "smoking gun" for detecting these particles at low interaction strengths. It's like finding a new type of fingerprint at a crime scene that everyone else overlooked.

5. The Final Verdict: New Limits on the Invisible

By combining all these new calculations, the authors drew a new map of where these ALPs can and cannot exist.

• Small Couplings: If the ALPs interact very weakly, the new "three-part decay" (electron + positron + photon) sets the strictest limits.
• Large Couplings: If they interact strongly, the "energy deposition" (how much heat they dump back into the star) sets the limits.
• The Fireball Zone: There is a middle ground where the ALPs create a fireball that turns into X-rays. This region is now ruled out by observations from the Pioneer Venus Orbiter.

Why Does This Matter?

Think of the universe as a giant puzzle. Supernovae are the pieces that fit together to tell us about the laws of physics. If we use the wrong recipe (old math) to calculate how these pieces fit, we might think a piece is missing when it's actually just hidden in a different spot.

This paper corrects the recipe. It tells us:

1. We missed a major production line (Semi-Compton scattering).
2. We missed a major decay channel (the three-part decay).
3. We have a better way to measure the "trap."

By fixing these errors, the authors have tightened the noose around these hypothetical particles. They haven't found the ALPs yet, but they have told us exactly where not to look, narrowing the search for new physics significantly. They even made their data public (like a shared spreadsheet) so other scientists can use their new, more accurate numbers to test their own theories.

Technical Summary

1. Problem Statement

Axion-like particles (ALPs) coupling to electrons ($g_{ae}$) are a candidate for new physics in the MeV–GeV mass range. Core-collapse supernovae (SNe) provide extreme environments (temperatures $T \sim 30\text{--}60$ MeV, electron chemical potentials $\mu_e \sim 100\text{--}500$ MeV) ideal for producing such particles. However, previous calculations of ALP emissivity in these environments suffered from several critical limitations:

• Neglected Channels: The semi-Compton scattering process ($\gamma e^- \to a e^-$) was largely ignored or considered subdominant, despite potentially dominating the production rate.
• Inconsistent Medium Effects: Previous works often treated the electron propagator in the medium inconsistently, sometimes using vacuum masses in loop calculations or incorrectly applying effective masses to spinor numerators.
• Loop-Induced Dominance: Some studies suggested that loop-induced ALP-photon couplings (leading to Primakoff production or photon coalescence) could dominate over direct electron couplings.
• Incomplete Constraints: Constraints derived from SN 1987A often relied on simplified cooling arguments or neglected specific decay channels (like $a \to e^+ e^- \gamma$) and the "fireball" formation regime.

The paper aims to provide a rigorous, first-principles recalculation of ALP production rates in a relativistic plasma and derive updated, comprehensive constraints on $g_{ae}$.

2. Methodology

The authors employ finite-temperature field theory to model the dense plasma of a protoneutron star (PNS).

• Medium Properties: They rigorously treat electrons and photons as having thermal masses ($m_{th}$ and $m_\gamma$) derived from the plasma frequency and self-energies. They demonstrate that for ultra-relativistic particles, the wavefunction renormalization $Z \to 1$, and the plasmino excitation is negligible. Crucially, they clarify that while the denominator of the electron propagator contains the thermal mass, the numerator (spinor structure) should not be treated as $\not{p} + m_{th}$; rather, medium corrections in the numerator are negligible for the processes considered.
• Coupling Equivalence: They prove the equivalence between pseudoscalar ($g_{ae} a \bar{e} \gamma_5 e$) and derivative couplings in the medium, noting that the equivalence relies on the vacuum mass $m_e$ in the transformation, not the thermal mass.
• Production Channels: The paper calculates emissivities for four main channels:
  1. Semi-Compton Scattering ($\gamma e^- \to a e^-$): Derived analytically and numerically, identifying it as the dominant channel in many regimes.
  2. Pair Annihilation ($e^+ e^- \to a \gamma$): Evaluated for high-temperature regions.
  3. Bremsstrahlung ($e^- N \to e^- N a$): Reduced from a 4D to a 2D numerical integral by analytically solving angular integrals, accounting for weak/strong screening and nucleon degeneracy.
  4. Coalescence ($e^+ e^- \to a$): Dominant for heavy ALPs ($m_a \gtrsim 2m_{th}$).
• Loop Suppression: They argue that loop-induced ALP-photon processes are suppressed in the medium. The energy cost to excite a virtual electron is the thermal mass $m_{th} \gg m_e$, not the vacuum mass. This suppresses the loop-induced emissivity by a factor of $(m_e/m_{th})^2 \sim 0.01$, rendering these channels subdominant.
• Absorption and Trapping: Using the principle of detailed balance, they derive absorption rates to determine the "trapping regime" where ALPs are thermalized within the star. They apply a new prescription for energy deposition in the trapping regime (based on radiative transfer theory) rather than simple luminosity subtraction.
• Observables: Constraints are derived from:
  • Cooling: Duration of the SN 1987A neutrino burst (Raffelt criterion).
  • Energy Deposition: Heating of the progenitor star by decaying ALPs (affecting low-energy SNe).
  • Decay Products: Positron injection (511 keV line), two-photon decay ($a \to \gamma \gamma$), and the previously neglected three-body decay ($a \to e^+ e^- \gamma$).
  • Fireball Formation: Thermalization of decay products into an X-ray emitting fireball (constrained by Pioneer Venus Orbiter data).

3. Key Contributions

1. Discovery of Dominant Channel: The paper establishes that semi-Compton scattering is the dominant production mechanism for ALPs coupled to electrons across a wide range of masses and plasma parameters, superseding bremsstrahlung in many relevant SN regions. This channel was previously neglected.
2. Theoretical Consistency: They correct the treatment of the electron propagator in the medium, showing that effective mass corrections belong in the denominator only, and proving that loop-induced processes are negligible due to thermal mass suppression.
3. Analytical and Numerical Efficiency: They derived new analytical expressions for bremsstrahlung emissivity in weak and strong screening limits and reduced the numerical integration dimensionality, enabling a complete scan of the $(T, \mu_e)$ parameter space.
4. Comprehensive Constraints: They provide the first comprehensive set of constraints on $g_{ae}$ for $m_a \gtrsim 1$ MeV, incorporating the $a \to e^+ e^- \gamma$ decay channel and the fireball formation regime.

4. Results

• Emissivity:
  • For light ALPs ($m_a \ll T$), semi-Compton scattering dominates, scaling roughly as $T^7/\mu_e^2$ in degenerate regimes.
  • For heavy ALPs ($m_a \gtrsim 2m_{th}$), electron-positron coalescence becomes the primary source.
  • Bremsstrahlung is generally subdominant due to nucleon degeneracy and screening effects.
  • The total emissivity is made available in a public repository as a function of $m_a$, $T$, and $\mu_e$.
• Constraints (SN 1987A):
  • Small Couplings ($g_{ae} \lesssim 10^{-10}$): The strongest bounds come from the non-observation of gamma-rays from the $a \to e^+ e^- \gamma$ decay. This channel provides constraints significantly stronger than the 511 keV line or $a \to \gamma \gamma$ bounds.
  • Large Couplings: Bounds are dominated by energy deposition arguments in low-energy supernovae. The new trapping regime prescription significantly alters the exclusion limits compared to previous "cooling-only" arguments.
  • Fireball Region: For very large couplings, decay products thermalize, forming a fireball that emits X-rays. This region is excluded by Pioneer Venus Orbiter (PVO) data, creating a distinct "fireball" exclusion zone in the parameter space.
  • Mass Dependence: The $a \to e^+ e^- \gamma$ bounds are nearly independent of the ALP mass, whereas $a \to \gamma \gamma$ bounds scale as $g_{ae} \propto m_a^{-1/2}$.

5. Significance

This work fundamentally updates the theoretical framework for ALP production in supernovae. By identifying semi-Compton scattering as the dominant channel and rigorously suppressing loop-induced processes via thermal mass effects, the authors provide a more accurate baseline for astrophysical constraints.

The derivation of new, tighter constraints—particularly from the $a \to e^+ e^- \gamma$ decay and the energy deposition in the trapping regime—narrows the viable parameter space for ALPs coupling to electrons in the MeV–GeV range. This has direct implications for:

• Particle Physics: Ruling out specific models of ALPs and dark sectors.
• Astrophysics: Refining our understanding of SN 1987A energetics and the explosion mechanism.
• Future Observations: Providing precise predictions for upcoming supernova neutrino and gamma-ray observations (e.g., from Hyper-Kamiokande or CTA).

The authors emphasize that previous bounds may have been over- or under-estimated due to the omission of the semi-Compton channel and the incorrect treatment of medium effects in loop diagrams. Their results serve as the new standard for SN-based ALP searches.