Improved supernova bounds on CP-even scalars: cooling and decay constraints

This paper presents significantly improved supernova constraints on CP-even scalars mixing with the Higgs boson by combining updated cooling calculations with new decay-based limits, thereby probing mixing angles down to $10^{-9}$and Yukawa couplings to$10^{-10}$, which extends sensitivity by over five orders of magnitude beyond existing collider bounds.

The Gist

Imagine the universe as a giant, cosmic detective story. For decades, scientists have been looking for "ghost particles"—new, invisible particles that might explain dark matter or why the universe is the way it is. Usually, we try to catch these ghosts by smashing atoms together in giant machines like the Large Hadron Collider (LHC). But if these ghosts are very light and interact very weakly with normal matter, our machines are like trying to catch a whisper with a megaphone; they just aren't sensitive enough.

So, this paper suggests we look elsewhere: Supernovae.

Think of a supernova (the explosion of a dying star) as the universe's most powerful, natural particle accelerator. It's a place so hot and dense that it's like a cosmic pressure cooker. If these ghostly particles exist, they should be born in the thousands inside the exploding star.

The authors of this paper are like detectives who have just upgraded their magnifying glass. They are looking at a specific type of ghost particle called a CP-even scalar (a fancy name for a particle that mixes with the Higgs boson). Here is what they did, explained simply:

1. The "Leaky Bucket" Problem (Cooling)

Imagine a hot cup of coffee (the core of the dying star) sitting in a cold room. Normally, it cools down by radiating heat (neutrinos). But if you poke a hole in the cup, it cools down too fast.

In 1987, we saw a supernova (SN1987a) and watched its "heat" (neutrinos) escape over about 10 seconds. This matched our predictions perfectly.

• The Old Theory: Previous scientists thought these ghost particles might escape, but they calculated the "leak" (how many particles are made) using a rough estimate.
• The New Discovery: These authors did a much more precise calculation. They realized that when you look closely at how these particles are made inside the star, the "leak" is actually 10 times bigger than we thought for light particles.
• The Result: If the leak is that big, the star would have cooled down way too fast, and the neutrinos would have arrived at Earth much sooner than they actually did. Since they didn't, we can now say with high confidence: "These ghost particles cannot exist if they are too weakly coupled." They have tightened the net around these particles significantly.

2. The "Positron Rain" (Decay)

Now, imagine some of these ghost particles escape the star and fly across the galaxy. Eventually, they might decay (break apart) into pairs of electrons and their antimatter twins, positrons.

• The Analogy: If too many of these particles decay, it would be like a sudden, massive rain of positrons hitting the center of our galaxy.
• The Evidence: We have satellites (like the INTEGRAL satellite) that watch the center of the Milky Way. They see a specific glow (511 keV gamma rays) caused by positrons hitting electrons.
• The Constraint: The authors calculated that if these ghost particles were too common, they would create a "positron storm" that would make the galaxy glow much brighter than it actually does. Since the galaxy isn't glowing that bright, these particles must be even rarer than we thought.

3. The "Energy Overload" (Low-Energy Supernovae)

Some supernovae are "low-energy" explosions. They are like a firecracker compared to a nuclear bomb. They don't have much extra energy to spare.

• The Scenario: If ghost particles are created and then decay inside the star's outer layers, they dump their energy back into the star.
• The Constraint: For a low-energy supernova, dumping extra energy is like adding a heavy weight to a balloon that's barely floating. It would change the explosion's shape and brightness. By looking at the few low-energy supernovae we've seen, the authors found that these ghost particles can't be dumping too much energy, which sets another limit on how common they can be.

4. The "Special Guest" (Hadrophilic Scalars)

The paper also looked at a special version of these particles that only talks to "nuclear matter" (protons and neutrons) and ignores electrons.

• The Twist: Because these particles don't decay into electrons, they live much longer and travel further before breaking apart. This makes them harder to catch via the "positron rain" method, but the "cooling" method still works very well to rule them out.

The Big Picture

The authors combined these three detective methods (Cooling, Positron Rain, and Energy Overload) to create a map of where these particles cannot exist.

• The Achievement: They have pushed the search for these particles into a territory that is 100,000 times more sensitive than what our best particle colliders on Earth can currently do.
• The Metaphor: If the particle colliders are looking for a needle in a haystack with a metal detector, this paper is like using a super-magnet to scan the entire barn, finding the needle even if it's hidden under a layer of dust.

In short: By doing better math on how these particles are born in exploding stars and how they behave afterward, the authors have proven that if these specific "ghost particles" exist, they are even more elusive than we hoped. They have ruled out a huge chunk of the "hunting ground" for dark matter models, forcing physicists to look in new places or rethink their theories.

Technical Summary

Here is a detailed technical summary of the paper "Improved supernova bounds on CP-even scalars: cooling and decay constraints."

1. Problem Statement

The search for physics Beyond the Standard Model (BSM) has shifted focus toward light, weakly-coupled "dark sector" particles, particularly in the MeV mass range, as collider experiments lose sensitivity in this regime. A specific class of interest is the CP-even scalar ($S$) that mixes with the Standard Model (SM) Higgs boson (the "Higgs portal").

Existing constraints on these scalars rely heavily on:

• Collider bounds: Limited to relatively large mixing angles ($\sin \theta$).
• Supernova (SN) cooling bounds: Previously derived for scalars, but often relying on approximations that may underestimate production rates.

The paper addresses the need for improved supernova constraints by:

1. Refining the calculation of scalar production rates in the proto-neutron star (PNS) core.
2. Incorporating new constraints based on the decay of these scalars into visible SM particles (positrons and energy deposition), which were previously unexplored for Higgs-portal scalars.
3. Extending the analysis to hadrophilic scalars (coupling only to nucleons), relevant for dark matter models.

2. Methodology

A. Scalar Production (Improved Calculation)

The dominant production mechanism in the SN core is nucleon-nucleon bremsstrahlung ($NN \to NNS$).

• Approach: The authors utilize the One-Pion Exchange (OPE) approximation in a vacuum.
• Key Improvement: Previous works often truncated the non-relativistic expansion of the nucleon propagator at leading order. The authors perform a consistent non-relativistic expansion including higher-order terms in the nucleon and scalar momenta ($p_i$ and $k_s$).
• Result: In the isospin-symmetric limit, leading-order amplitudes cancel out. By retaining next-to-leading-order (NLO) terms, the authors find the production rate increases by an order of magnitude for scalars with mass $m_s < 100$ MeV.
• Correction: They also correct the scalar-pion interaction vertex, which differs from previous literature (Ref. [34]), significantly altering contributions from diagrams with internal pion propagators.

B. Absorption and Decay

The authors model the fate of produced scalars using the optical depth ($\tau$), accounting for:

• Decay: $S \to e^+e^-$, $\mu^+\mu^-$, $\pi^+\pi^-$, $\pi^0\pi^0$. The decay width depends on the mixing angle $\sin \theta$ and mass $m_s$.
• Re-absorption: Inverse bremsstrahlung ($NNS \to NN$).
• Regimes:
  • Free-streaming: Low coupling; scalars escape the core.
  • Trapping: High coupling; scalars are re-absorbed or decay before escaping, effectively thermalizing. The authors use a black-body approximation at a "decoupling radius" ($R_b$) where $\tau \approx 2/3$ to model this regime.

C. Constraints Applied

Three distinct observational bounds are applied to the parameter space ($m_s$ vs. $\sin \theta$):

1. SN1987a Cooling Bound (Raffelt Criterion):
   • New particles must not carry away energy faster than neutrinos.
   • Limit: $L_{NP} \leq 3 \times 10^{52}$ erg/s.
   • Applied to both free-streaming and trapped regimes.
2. Low-Energy Supernova (LE-SN) Bound:
   • Targets Type II-P supernovae with low explosion energies ($\sim 10^{50}$ erg).
   • Constraint: Energy deposited by decaying scalars into the progenitor mantle must not exceed the observed explosion energy.
   • Accounts for the fraction of decay energy ($f_{dep}$) that thermalizes in the mantle (vs. escaping as neutrinos).
3. Galactic Positron Bound:
   • Targets the 511 keV gamma-ray flux from the galactic center (observed by INTEGRAL).
   • Constraint: The total rate of positrons produced by scalar decays from galactic supernovae must not exceed $\sim 4 \times 10^{43} \text{ s}^{-1}$.
   • This constrains the region where scalars escape the SN but decay before reaching the interstellar medium.

D. Hadrophilic Scalar Extension

The analysis is extended to a model where the scalar couples to nucleons via a dimension-5 operator ($S N G^2$) but has negligible coupling to leptons. This changes the decay channels (no $e^+e^-$ or $\mu^+\mu^-$) and decay widths, leading to different trapping dynamics.

3. Key Contributions

1. Order-of-Magnitude Improvement in Production Rate: By correctly expanding the nucleon propagator to include momentum-dependent terms, the authors demonstrate that scalar production is significantly higher than previously estimated for $m_s < 100$ MeV. This tightens the cooling bound by more than an order of magnitude.
2. Novel Decay Constraints:
   • Positron Bound: First application of the galactic 511 keV flux limit to Higgs-portal scalars. This excludes a region of parameter space where cooling is allowed, but decays produce too many positrons.
   • LE-SN Bound: Utilizes low-energy supernova observations to constrain the upper limit of the coupling (trapping regime) where energy deposition in the mantle would be excessive.
3. Hadrophilic Scalar Bounds: Provides the first comprehensive SN cooling bounds for a simple hadrophilic scalar model, showing that the absence of leptonic decays leads to stronger constraints in the trapping regime due to longer lifetimes.

4. Results

• Mixing Angle Sensitivity: The combined constraints probe mixing angles as small as $\sin \theta \sim 10^{-9}$ for MeV-scale scalars.
• Comparison to Colliders: This sensitivity is more than five orders of magnitude below existing collider bounds.
• Parameter Space Coverage: The combination of astrophysical (SN) and collider probes covers over nine orders of magnitude in coupling strength, effectively testing the majority of the parameter space motivated by dark matter models.
• Specific Features:
  • A sharp cutoff in the positron bound occurs at $m_s = 2m_\mu$ (approx. 211 MeV) because the decay channel shifts to muons, drastically reducing the scalar lifetime and causing decays to occur inside the SN envelope rather than in the galaxy.
  • For hadrophilic scalars, the bound is stronger in the trapping regime because the lack of leptonic decays increases the mean free path, allowing more energy to be trapped and re-radiated.

5. Significance

• Robustness of Astrophysical Probes: The paper demonstrates that supernovae remain the most powerful tool for probing weakly coupled light scalars, outperforming current and near-future terrestrial experiments.
• Theoretical Precision: The work highlights the importance of higher-order terms in non-relativistic expansions for accurate astrophysical constraints. Previous approximations (soft-radiation or leading-order expansion) were insufficient for the dense SN environment.
• Dark Matter Implications: Since many dark matter models (e.g., thermal freeze-out via Higgs portal) rely on these specific couplings, the results significantly narrow the viable parameter space for such models.
• Complementarity: The study establishes a clear synergy between different astrophysical signals (cooling, energy deposition, and decay products) and collider searches, providing a multi-faceted approach to constraining BSM physics.

Uncertainties: The authors note that theoretical uncertainties (e.g., OPE approximation validity, SN profile variations between "hot" and "cold" models) introduce an order-of-magnitude uncertainty in the bounds. However, the overlap of different constraints and the use of conservative profiles suggest the results are robust to within a factor of 2–3 in the coupling.