Comprehensive neutrino light curves and spectra: from pre-supernova evolution to early supernova phase

This paper presents the first systematic study of neutrino emissions from massive stars (10–40 $M_\odot$) spanning pre-supernova evolution to the early core-collapse phase, revealing that neutrino luminosities and spectra strongly correlate with progenitor compactness and core mass, thereby offering a robust method to infer internal stellar structure through joint observational analysis.

The Gist

Imagine a massive star as a giant, cosmic pressure cooker. For millions of years, it cooks up heavier and heavier elements in its core, layer by layer, like an onion made of nuclear fire. Eventually, the fuel runs out, the pressure cooker can't hold the weight anymore, and it collapses in on itself, triggering a supernova—a stellar explosion so bright it can outshine an entire galaxy.

For decades, scientists have tried to predict exactly how these stars behave right before they explode. But it's like trying to guess the contents of a sealed, black box just by looking at the outside. We know the star's mass, but we don't know the exact structure of its deep interior.

This paper is like a super-advanced X-ray machine for dying stars. The authors, a team of astrophysicists, have created a massive simulation that tracks a star's neutrino emissions (ghostly particles that pass through everything) from hundreds of years before the explosion, all the way through the first split-second of the blast.

Here is the story of their discovery, broken down into simple concepts:

1. The "Ghostly" Messenger: Neutrinos

Neutrinos are like invisible messengers. They are produced in the star's core and escape immediately, carrying information about what's happening deep inside. Unlike light, which gets trapped and scattered by the star's outer layers, neutrinos tell us the truth about the core's temperature, density, and structure.

2. The Two-Act Play: Before and After the Crash

The authors looked at the neutrino signal in two distinct acts:

• Act I: The Pre-Supernova (The "Winding Up"):
  Imagine the star is a clock winding up. As it gets closer to the end, the core gets hotter and denser. The authors found that the total number of neutrinos emitted during the final days depends on how "compact" the star's core is.

  • The Analogy: Think of the core as a sponge. If the sponge is very dense and tight (high compactness), it squeezes out neutrinos in a specific pattern. If it's fluffier, the pattern is different. By counting the neutrinos in the final hours, we can tell if the star's core was a tight sponge or a loose one.

• Act II: The Early Supernova (The "Bounce"):
  When the core finally collapses, it hits a point of no return and "bounces" back, sending a shockwave outward. This happens in a fraction of a second.

  • The Analogy: Imagine a rubber ball hitting the floor. The initial bounce is universal—it happens the same way for almost any ball. The authors found that the very first burst of neutrinos (the "neutronization burst") is surprisingly similar for all stars. However, the accretion phase (when the star's outer layers rain down onto the new core) acts like a different rainstorm for every star. The heavier the rain (accretion), the hotter the core gets, and the more neutrinos are released.

3. The Secret Code: Compactness and Core Mass

The big breakthrough in this paper is finding a "Rosetta Stone" to decode the star's interior. They found two main keys:

• Compactness ($\xi_{2.5}$): How tightly packed the core is.
• Carbon-Oxygen Core Mass ($M_{CO}$): How big the star's "kitchen" was before it started cooking the final ingredients.

They discovered that if you count the neutrinos over a long period (years), the number tells you about the size of the kitchen ($M_{CO}$). If you count them only in the final day, the number tells you how tight the core is ($\xi_{2.5}$).

4. The Real-World Test: Can We Actually See This?

The authors didn't just run simulations; they asked, "If a star exploded nearby, could our detectors actually see this?"

• They simulated what detectors like Super-Kamiokande (a giant tank of water in Japan) and JUNO (a massive liquid scintillator in China) would see.
• The Result: Yes! Even with the noise of background radiation (like cosmic rays), if a star like Betelgeuse (which is relatively close to us) were to explode, our detectors could spot the "pre-supernova" neutrinos days in advance.
• The "False Alarm" Check: They used a statistical method (like checking if a noise in a room is just a cough or a gunshot) to ensure we wouldn't mistake background noise for a real star dying. They found that for nearby stars, we could get a "heads up" alert hours or even days before the explosion.

5. Why This Matters: The "Synergy"

The most exciting part is the combination of the two phases.

• Pre-Supernova Neutrinos tell us about the star's structure before the explosion gets messy.
• Early Supernova Neutrinos tell us about the explosion dynamics after the core bounces.

If we observe both, we can cross-check our theories. If the pre-explosion data says the core is "tight," but the explosion data suggests it was "loose," we know our physics models are wrong. It's like checking a suspect's alibi (pre-explosion) against the crime scene evidence (explosion). If they don't match, we have to rewrite the story of how stars die.

The Bottom Line

This paper is a roadmap for the future of astronomy. It tells us that by listening to the "ghostly whispers" (neutrinos) of a dying star, we can:

1. Predict when a star is about to explode (giving us time to point our telescopes at it).
2. Diagnose the internal structure of the star, which we can never see with normal light.
3. Test our understanding of physics under the most extreme conditions in the universe.

It turns the terrifying, chaotic death of a star into a readable, understandable story, written in the language of neutrinos.

Technical Summary

Here is a detailed technical summary of the paper "Comprehensive Neutrino Light Curves and Spectra: From Pre-supernova Evolution to Early Supernova Phase" by Kato et al.

1. Problem Statement

Massive stars ($M \gtrsim 8 M_\odot$) undergo complex evolution leading to core-collapse supernovae (CCSNe). However, the final internal structure of these progenitors (specifically the iron core and surrounding shells) is subject to large uncertainties due to poorly constrained input physics, such as mass-loss rates, binary interactions, and mixing processes.

• The Gap: While neutrino emission is the primary messenger from the stellar interior, previous studies have largely treated the pre-supernova (preSN) phase and the post-bounce supernova (SN) phase separately. There is a lack of systematic, continuous modeling linking the late evolutionary stages (hundreds of years prior to collapse) to the early explosion phase (first 200 ms).
• The Challenge: Conventional electromagnetic observations cannot probe the deep interior. Neutrino detection offers a direct window, but to utilize this, theoretical models must consistently predict neutrino luminosities and spectra across the entire timeline to establish robust correlations between observable neutrino signals and progenitor structural parameters (e.g., compactness, core mass).

2. Methodology

The authors performed a systematic study using 30 progenitor models with Zero-Age Main Sequence (ZAMS) masses ranging from 10 to $40 M_\odot$ at solar metallicity. The methodology involved a two-step approach:

A. Stellar Evolution and Collapse (PreSN Phase)

• Codes: Two distinct stellar evolution codes were used to test robustness:
  • HOSHI: A 1D hydrodynamics code using the Ledoux criterion for convection and mixing-length theory. It defines core collapse onset when the central temperature $T_c \sim 10^{9.9}$ K.
  • MESA: A public stellar evolution code. It defines collapse onset when the iron core material reaches velocities $> 100$ km/s.
• Physics: Both codes utilized a 300-species nuclear reaction network. The HOSHI models included wind mass loss, while MESA models (for simplicity in this suite) did not.
• Neutrino Calculation (PreSN): Neutrino luminosities and spectra were calculated in a post-processing step assuming a stationary state (neutrino timescale $\ll$ dynamical timescale). Six emission processes were included: pair annihilation, $\beta^\pm$ decay, electron/positron capture on nuclei, and capture on free protons.

B. Dynamical Evolution (Early SN Phase)

• Code: A general-relativistic Boltzmann radiation-hydrodynamics code was used for the post-collapse phase.
• Physics: Spherically symmetric (1D) simulations were run up to 200 ms post-bounce. The code solves the Boltzmann equations for three neutrino species ($\nu_e, \bar{\nu}_e, \nu_x$) with full energy and angular dependence.
• Limitations: The study is restricted to 1D to manage computational costs; multi-dimensional effects (SASI, convection) and collective neutrino oscillations are acknowledged as future work but are noted to have limited impact on the early (first 200 ms) universal features.

C. Observational Simulation

• Detectors: Event rates were estimated for SK, HK, KamLAND, JUNO, and DUNE.
• Oscillations: Vacuum and MSW effects were included (assuming Normal and Inverted mass orderings).
• Detection Logic: A False Alarm Rate (FAR) approach was used to determine the "first detection time" (alert time) based on sliding time windows (12h/24h) and background rates.

3. Key Contributions

1. First Continuous Modeling: The paper provides the first systematic, continuous calculation of neutrino light curves and spectra spanning from $\sim 10^9$ seconds (hundreds of years) before collapse to 200 ms after bounce.
2. Robust Correlations: It identifies and validates strong correlations between neutrino observables and key progenitor structural parameters ($\xi_{2.5}$ and $M_{CO}$) that persist across different stellar evolution codes.
3. Observational Feasibility: It evaluates the detectability of these correlations under realistic conditions (including oscillations, energy thresholds, and background noise) using a FAR-based alert strategy.
4. Public Data Release: All calculated time-series data for 30 models are made publicly available to facilitate future multi-messenger studies.

4. Key Results

A. Neutrino Emission Characteristics

• PreSN Phase:
  • $\nu_x$ and $\bar{\nu}_e$: Emission is dominated by pair annihilation. Luminosity correlates positively with the size of the core (extended hot cores like $40 M_\odot$ emit more).
  • $\nu_e$: Emission shifts from pair annihilation to electron capture (EC) as the core degenerates. Just prior to collapse, models with compact, high-density cores (high electron chemical potential $\mu_e$) exhibit the highest $\nu_e$ luminosity due to efficient EC.
  • Energy: Average energies increase monotonically. $\nu_e$ average energy is highly sensitive to $\mu_e$ (and thus core compactness), while $\bar{\nu}_e$ and $\nu_x$ are driven by temperature.
• Early SN Phase (< 200 ms):
  • Universality: The $\nu_e$ neutronization burst is largely universal across models, governed by the self-similar collapse of the inner iron core.
  • Divergence: $\bar{\nu}_e$ and $\nu_x$ luminosities diverge significantly based on the progenitor's outer core structure. The accretion rate onto the Proto-Neutron Star (PNS) drives compressive heating, making these emissions highly sensitive to the density gradients of the progenitor.

B. Correlations with Structural Parameters

The study defines two key parameters:

1. Compactness ($\xi_{2.5}$): Mass enclosed within a specific radius at $2.5 M_\odot$.
2. CO Core Mass ($M_{CO}$): Mass coordinate at the base of the oxygen-rich region.

• PreSN Correlations:
  • Integrating neutrino numbers over the final day ($t_i \sim -10^5$ s) yields a strong correlation with $\xi_{2.5}$.
  • Integrating over longer durations ($t_i \sim -10^9$ s) yields a strong correlation with $M_{CO}$.
  • Physical Reason: Short-term emission reflects the immediate core structure (compactness), while long-term emission integrates the evolutionary history (core mass growth).
• Early SN Correlations:
  • The peak luminosities of $\bar{\nu}_e$ and $\nu_x$ and the width of the $\nu_e$ burst show a robust correlation with $\xi_{2.5}$.
  • This correlation is insensitive to the specific stellar evolution code (HOSHI vs. MESA) during the first 200 ms, as the inner core structure converges (Mazurek's law) and the early accretion phase is dominated by universal physics.

C. Observational Prospects

• Alert Times: For a nearby progenitor (200 pc), JUNO (liquid scintillator, low threshold) provides the earliest alert (days to hours before bounce), followed by KamLAND. Water Cherenkov detectors (SK, HK) have higher thresholds and detect later.
• Event Counts:
  • At 200 pc: JUNO detects $\sim 10^3-10^4$ preSN events; SK/HK detect $\sim 10^6-10^7$ SN events.
  • At 1 kpc: JUNO can still detect $\sim 10-100$ preSN events and issue alerts, though statistical reconstruction becomes harder.
• Survival of Correlations: Even after applying realistic constraints (oscillations, energy thresholds, alert time windows), the theoretical correlation between the time-integrated preSN neutrino count and $\xi_{2.5}$ remains largely intact.

5. Significance

• Independent Probes: This work demonstrates that preSN neutrinos and early SN neutrinos act as independent probes of the progenitor's internal structure. PreSN neutrinos constrain the core compactness and mass before explosion uncertainties set in, while early SN neutrinos confirm these constraints immediately post-bounce.
• Model Validation: A discrepancy between the $\xi_{2.5}$ derived from preSN neutrinos and that derived from early SN neutrinos would immediately expose flaws in either the stellar evolution models (e.g., mixing, mass loss) or the explosion physics.
• Multi-Messenger Strategy: The paper establishes a framework for using future neutrino alerts (from JUNO, HK, etc.) to trigger electromagnetic and gravitational wave observations, transforming a single supernova detection into a comprehensive diagnostic of massive star evolution.
• Data Utility: The public release of time-series data for 30 models provides a critical resource for the community to refine detection algorithms and interpret future real-world observations.

In conclusion, the paper argues that continuous neutrino monitoring from the late preSN phase through the early SN phase offers a unique, high-precision method to constrain the internal structure of massive stars, bypassing the uncertainties inherent in modeling the explosion mechanism itself.