🔭 astrophysics

Probing the potential high-energy messengers of the anticipated T Coronae Borealis outburst

This paper employs 3D hydrodynamical simulations coupled with diffusive shock acceleration models to predict the gamma-ray and neutrino emission from the anticipated T Coronae Borealis outburst, revealing that early high-energy signals are dominated by ejecta and the accretion disk, with potential PeV particle acceleration and feasible neutrino detection in high-energy, high-density scenarios.

Original authors: O. Petruk, T. Kuzyo, S. Orlando, L. Chomiuk, F. Bocchino, M. Miceli, S. Ustamujic

Published 2026-03-17

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: O. Petruk, T. Kuzyo, S. Orlando, L. Chomiuk, F. Bocchino, M. Miceli, S. Ustamujic

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a cosmic fireworks display that hasn't happened in nearly 80 years. This is T Coronae Borealis (T CrB), a "recurring nova" located relatively close to Earth. It's a binary star system where a dense, dead star (a white dwarf) is greedily stealing gas from its bloated, aging neighbor (a red giant). Every few decades, the white dwarf gets so full of stolen gas that it explodes in a thermonuclear blast.

This paper is a crystal ball for astronomers. The authors are using supercomputers to predict exactly what this explosion will look like not just in visible light, but in the invisible, high-energy "messengers" of the universe: Gamma-rays (extremely energetic light) and Neutrinos (ghostly particles that pass through everything).

Here is a breakdown of their findings using simple analogies:

1. The Setting: A Cosmic Obstacle Course

The explosion doesn't happen in empty space. The white dwarf is surrounded by a complex "circumbinary medium" (CBM), which is like a messy room filled with different obstacles:

The Red Giant's Wind: A steady breeze of gas blowing out from the neighbor star.
The Equatorial Density Enhancement (EDE): A thick, dense ring of gas around the star's equator, like a heavy belt.
The Accretion Disk: A flat, spinning disk of gas falling onto the white dwarf, like a whirlpool.

When the nova explodes, it sends a shockwave (a giant wall of compressed gas) flying outward at thousands of kilometers per second. The authors simulated how this shockwave crashes into these different obstacles.

2. The Accelerator: A Cosmic Particle Smasher

The main goal of the paper is to see how this shockwave acts as a particle accelerator. Think of the shockwave as a giant cosmic pinball machine.

The Bumpers: As the shockwave hits the gas, it creates turbulent magnetic fields (invisible force fields).
The Balls: Protons and electrons get trapped in these fields, bouncing back and forth across the shockwave. Every time they cross, they get a kick of energy, speeding up to near the speed of light.
The Twist: The authors discovered something new. Usually, scientists assume the gas behind the shockwave flows smoothly. But in their 3D simulations, they found the flow is bumpy and uneven.
- Analogy: Imagine a river flowing over rocks. The water doesn't flow in a straight line; it swirls and speeds up in some spots. The authors found that these "swirls" (velocity gradients) actually help particles get accelerated to much higher energies than previously thought—potentially reaching PeV (Peta-electronvolt) scales, which is the energy of the most powerful cosmic rays known!

3. The Messengers: What Will We See?

Gamma-Rays (The Flash)

When these super-fast protons smash into other gas particles, they create a flash of high-energy gamma-rays.

The Early Show: In the first few hours, the explosion is dominated by the ejecta (the debris from the explosion itself) and the accretion disk. The disk acts like a booster seat, making the gamma-rays much brighter and more energetic during this short window.
The Duration: The paper predicts we will be able to see these gamma-rays with telescopes like Fermi-LAT (which sees lower energy) for about one month, and with ground-based telescopes like H.E.S.S. or CTAO (which see higher energy) for about one week.
The Shape: Because the environment is messy (with the disk and the ring), the gamma-ray light curve won't be a smooth slide down; it will have bumps and breaks, acting like a fingerprint of the environment.

Neutrinos (The Ghosts)

Neutrinos are created in the same collisions that make gamma-rays, but they are much harder to catch because they rarely interact with anything.

The Catch: Detecting neutrinos from T CrB is a long shot, but not impossible.
The Condition: It depends heavily on the density of the gas the shockwave hits.
- If the accretion disk is very dense and the explosion is very powerful (like in their "RUN10" model), we might catch a few neutrinos with detectors like IceCube or KM3NeT.
- If the disk is thin or the explosion is weaker, the signal will be too faint to detect.
The Verdict: The authors are more optimistic than previous studies. They believe that if the explosion is strong enough, neutrino detectors might actually see something, which would be a historic first for a nova.

4. The "Gotchas" (Absorption)

There is a catch. The white dwarf is so bright and hot right after the explosion that it acts like a giant lightbulb.

The Fog: This intense light can actually "eat" some of the gamma-rays trying to escape, turning them into electron-positron pairs.
The Result: This might dim the very highest energy gamma-rays (TeV range) during the first few days, making them harder to see, while the lower energy ones (GeV range) pass through fine.

Summary: Why Does This Matter?

This paper is a roadmap for the next few years.

Timing is everything: Astronomers need to be ready to look at T CrB immediately after it erupts (expected around 2026). The first 24 hours are critical to see the interaction with the accretion disk.
Multi-Messenger Astronomy: By looking at light (gamma-rays) and ghosts (neutrinos) together, we can figure out exactly how much energy the explosion had and how dense the environment was.
New Physics: The discovery that "bumpy" flows in the gas can boost particle energy to PeV scales changes how we understand how cosmic rays are born in the universe.

In short, the authors have built a detailed simulation of a cosmic explosion to tell us exactly what to look for, when to look, and what we might learn about the most energetic processes in our galaxy.

1. Problem and Motivation

T Coronae Borealis (T CrB) is a nearby recurrent nova (a symbiotic binary system consisting of a $1.37 M_\odot$ white dwarf and a $0.69 M_\odot$ red giant) expected to undergo a thermonuclear outburst around 2026. While novae are known to emit GeV gamma-rays, the mechanisms of particle acceleration and the potential for neutrino emission in such systems remain under-constrained.

The Gap: Previous studies relied on semi-analytical or 1D models. There is a lack of detailed 3D hydrodynamical simulations that account for the complex circumbinary medium (CBM)—specifically the red giant wind, equatorial density enhancement (EDE), and accretion disk—and their impact on diffusive shock acceleration (DSA).
The Goal: To predict the high-energy (gamma-ray and neutrino) emission from T CrB during the first months of its outburst, assess its detectability by current observatories (Fermi-LAT, CTA, IceCube, KM3NeT), and investigate the role of flow non-uniformity in particle acceleration.

2. Methodology

The authors employed a multi-stage numerical approach combining hydrodynamics with particle acceleration physics:

3D Hydrodynamical (HD) Simulations:
- Utilized the PLUTO code to simulate the nova outburst based on models from a companion paper (Paper I).
- Environment: The CBM includes a red giant wind with an equatorial density enhancement and an accretion disk around the white dwarf.
- Models: Four specific models were analyzed (Table 1), varying in explosion energy ( $E_{bw}$ $E_{b w}$ ), EDE density ( $n_{ede}$ $n_{e d e}$ ), and the presence of an accretion disk.
  - RUN04: Reference model ( $3 \times 10^{43}$ erg, disk present).
  - RUN03: High density EDE ( $10 \times$ denser).
  - RUN10: High explosion energy ( $10^{44}$ erg) and high density.
  - RUN14: No accretion disk.
Particle Acceleration Model:
- Diffusive Shock Acceleration (DSA): Applied to the forward shock.
- Magnetic Field (MF): Post-processed using the Bell (2004) non-resonant instability to estimate amplified MF strength ( $\delta B$ ).
- Maximum Momentum ( $p_{max}$ ): Calculated as the minimum of three limits: time-limited acceleration ( $p_{m1}$ ), Bell instability growth time ( $p_{m2}$ ), and radiative losses ( $p_{m3}$ ).
- Key Innovation (Flow Gradients): Unlike standard models assuming uniform downstream flow, this study incorporates spatially varying flow velocity gradients ( $\mu = du/dx$ ) in the post-shock region. This modifies the particle momentum distribution from a simple power law to a function involving an exponential cutoff and a gradient-dependent term.
Emission Synthesis:
- Multi-zone Approach: Emissivities were calculated for every grid cell and summed.
- Hadronic: Proton-proton ($pp$) interactions producing gamma-rays and neutrinos (using the Naima/Pythia8 framework).
- Leptonic: Inverse-Compton (IC) scattering off the Cosmic Microwave Background (CMB) and Red Giant (RG) wind photons.
- Absorption: Considered $\gamma$ - $\gamma$ absorption by the nova's own radiation field and $\gamma$ - $p$ absorption in the dense accretion disk.

3. Key Contributions

3D Hydrodynamic Modeling of High-Energy Messengers: First application of full 3D HD simulations to synthesize gamma-ray and neutrino light curves and spectra for T CrB, capturing the complex geometry of the EDE and accretion disk.
Impact of Flow Non-Uniformity: Demonstrated that accounting for downstream velocity gradients ( $\mu \neq 0$ ) significantly alters the particle spectrum. Negative gradients (found in ejecta protrusions/bubbles) lead to harder spectra and allow particles to reach energies an order of magnitude higher than predicted by uniform flow models (reaching the PeV scale in high-energy scenarios).
Component Decomposition: Detailed the temporal evolution of emission sources, distinguishing between the shocked ejecta, the accretion disk, and the ambient CBM.
Neutrino Predictions: Provided specific neutrino detection estimates for IceCube and KM3NeT based on 3D models, challenging previous pessimistic estimates.

4. Key Results

A. Particle Acceleration and Spectra

Maximum Energies: In the reference model (RUN04), protons can reach $\sim 2$ TeV at 153 days. In high-energy models (RUN10), protons can reach the PeV scale ( $10^{15}$ eV) in the first few hours due to high densities and amplified magnetic fields.
Gradient Effect: Including flow gradients ( $\mu \neq 0$ ) increases the effective maximum energy ( $E_{max,eff}$ ) by a factor of $\sim 10$ compared to the uniform assumption ( $\mu=0$ ).
Spectral Features: The multi-zone nature of the remnant creates spectral breaks (e.g., a break from $E^{-2}$ to $E^{-2.8}$ at a few TeV) and bumps in the integrated spectra.

B. Gamma-Ray Emission

Dominance: Early emission (first days/weeks) is dominated by hadronic processes (proton-proton interactions) within the ejecta.
Disk Influence: The accretion disk significantly boosts the gamma-ray flux and particle normalization during the first hours. Once the shock exits the disk (approx. 8 hours in RUN04), the flux drops sharply.
Detectability:
- GeV Range: Detectable by Fermi-LAT for approximately one month.
- TeV Range: Detectable by Cherenkov telescopes (H.E.S.S., VERITAS, CTAO) for approximately one week.
- PeV Range: Potentially detectable by LHAASO in high-energy models (RUN10) for a very short duration.
Absorption: During the first day, the intense radiation field from the nuclear-burning white dwarf causes significant $\gamma$ - $\gamma$ absorption of TeV photons and enhances IC emission, potentially making leptonic emission dominant for the first week.

C. Neutrino Emission

Flux: Neutrino flux is directly linked to the explosion energy and ambient density.
Detection Probability:
- RUN10 (High Energy): High probability of detection by IceCube and KM3NeT (tens of events in the first 4 hours).
- RUN04/RUN03: Lower flux, but still potentially detectable by IceCube Main Array.
- RUN14 (No Disk): Negligible flux; detection unlikely.
Conclusion: Neutrino detection is feasible primarily in models with high explosion energy ( $>10^{44}$ erg) and high ambient density (presence of a dense accretion disk).

D. Spatial Distribution

Gamma-ray images show that emission is not uniform. The highest energies are achieved in "ejecta protrusions" (bubbles) where flow gradients are negative and in regions where the shock interacts with the dense accretion disk.

5. Significance and Implications

Multi-Messenger Strategy: The paper argues that multi-messenger observations must prioritize the first 24 hours to capture the critical interaction between the shock and the accretion disk. This phase is essential for distinguishing between different CBM configurations.
Re-evaluating Nova Potential: The study suggests that novae, despite having lower explosion energies than supernovae, can accelerate particles to near-PeV energies if the environment is dense and magnetic fields are highly amplified. This positions recurrent novae as potential contributors to the Galactic cosmic ray population.
Observational Guidance: The models provide specific light curves and spectra to guide observations by Fermi-LAT, CTA, IceCube, and KM3NeT for the anticipated 2026 T CrB outburst.
Methodological Advance: The inclusion of flow velocity gradients in DSA models offers a more realistic framework for high-energy emission synthesis in transient astrophysical events, applicable to other novae and interacting supernovae.

In summary, the paper predicts that the 2026 T CrB outburst will be a bright, multi-messenger event, with a high probability of neutrino detection if the explosion is energetic and the accretion disk is dense, offering a unique laboratory to study particle acceleration in complex, non-uniform environments.