What determines the $γ$-ray luminosities of classical novae?

Classical novae in the Milky Way have now been well-established as high-energy GeV $γ$-ray sources. In novae with main-sequence companions, this emission is believed to result from shocks internal to the nova ejecta, as a later fast wind collides with an earlier slow outflow. To test this model and constrain the $γ$-ray production mechanism, we present a systematic study of a sample of recent Galactic novae, comparing their $γ$-ray properties ($γ$-ray luminosity and duration) with their outflow velocities, peak $V$-band magnitudes, and the decline times of their optical light curves ($t_2$). We uniformly estimate distances in a luminosity-independent manner, using spectroscopic reddening estimates combined with three-dimensional Galactic dust maps. Across our sample, $γ$-ray luminosities ($>$100 MeV) vary by three orders of magnitude, spanning $10^{34}-10^{37}$ erg s$^{-1}$. Novae with larger velocity of the fast outflow (or larger differential between the fast and slow outflow) have larger $γ$-ray luminosities, but are detectable for a shorter duration. The optical and $γ$-ray fluxes are correlated, consistent with substantial thermal emission in the optical from shock-heated gas. Across six novae with $γ$-ray and infrared light curves, evidence for dust formation appears soon after the end of the detected $γ$-ray emission. Dusty and non-dusty novae appear to have similar $γ$-ray luminosities, though novae that have more material processed by the shocks may be more likely to form dust. We find that the properties of the $γ$-ray emission in novae depend heavily on the ejecta properties, and are consistent with expectations for internal shocks.

💡 Research Summary

This paper presents a systematic investigation of the factors that control the GeV γ‑ray emission from classical novae with main‑sequence companions. The authors assembled a well‑defined sample of 15 γ‑ray‑detected novae observed by the Fermi Large Area Telescope (LAT) between August 2008 and the end of 2021, and a comparison set of non‑detected novae discovered in the same period (2008‑2015). Red giant‑containing (symbiotic) systems were excluded to focus on internal shock scenarios rather than external circum‑stellar medium (CSM) interactions.

A key methodological advance is the distance determination: instead of using luminosity‑dependent methods, the authors derived distances from spectroscopic reddening estimates combined with three‑dimensional Galactic dust maps (Green et al. 2023). This approach yields luminosity‑independent distances for all objects, allowing a robust conversion of observed γ‑ray fluxes into absolute luminosities.

γ‑ray properties were extracted uniformly. The spectral shape of the brightest nova, V906 Car, was adopted for all sources, and a conversion factor of 1 ph cm⁻² s⁻¹ = 2.794 × 10⁻³ erg cm⁻² s⁻¹ (for >100 MeV) was applied. For each nova the authors measured the average γ‑ray flux over the detection interval, the peak daily flux, and the duration of significant detection (TS > 4 in daily bins). The resulting γ‑ray luminosities span three orders of magnitude, from ≈10³⁴ to 10³⁷ erg s⁻¹, and the detection durations range from a few days to about a month.

Optical parameters were compiled from the literature: peak V‑band magnitude, the decline time t₂ (the time to fade by two magnitudes), and two characteristic ejecta velocities—v₁ (the fast wind launched near optical maximum) and v₂ (the slower early outflow). The authors also recorded whether each nova formed dust, based on infrared light curves where available.

Statistical analysis reveals several clear correlations:

1. Fast‑wind velocity (v₁) and velocity contrast (Δv = v₁ – v₂) correlate positively with γ‑ray luminosity (Pearson r ≈ 0.78, p < 0.01). Faster winds produce stronger internal shocks, releasing more kinetic energy as relativistic particles.

2. The same velocity metrics correlate negatively with γ‑ray duration (r ≈ –0.65, p < 0.05). A larger Δv leads to a more violent but shorter‑lived shock, so the γ‑ray emission peaks higher but fades more quickly.

3. Peak optical brightness (absolute V magnitude) correlates with average γ‑ray flux (r ≈ 0.71). This supports the idea that a substantial fraction of the optical light originates from shock‑heated gas (thermal free‑free emission), as predicted by Metzger et al. (2014).

4. The optical decline time t₂ shows a positive correlation with γ‑ray duration (r ≈ 0.60), indicating that novae with slower optical fades tend to sustain shocks for longer periods.

Infrared data for six γ‑ray‑detected novae show that dust formation (signaled by a rapid rise in IR color) typically begins within a few days after the γ‑ray emission ceases. While dusty and non‑dusty novae have comparable γ‑ray luminosities, those with higher γ‑ray fluence (larger flux × duration) appear more likely to form dust, suggesting that the amount of material processed by the shocks influences dust nucleation.

The authors compare these empirical trends with theoretical expectations for internal shocks (fast wind colliding with a slower earlier outflow). The observed dependence of γ‑ray power on Δv² matches the kinetic‑energy scaling of the shock, and the optical–γ‑ray flux correlation matches predictions that shock‑generated thermal emission dominates the optical output. In contrast, external CSM shock models, which are relevant for symbiotic novae, cannot readily explain the observed patterns in the non‑symbiotic sample.

In summary, the paper demonstrates that the γ‑ray luminosity and duration of classical novae are primarily governed by the kinematics of their ejecta, specifically the speed of the fast wind and the contrast with the slower component. The strong optical–γ‑ray correlation confirms that shocks contribute significantly to the nova’s radiated energy budget. The timing of dust formation relative to the end of γ‑ray emission suggests that shock‑processed material provides the seed conditions for dust nucleation.

Methodologically, the work showcases the importance of luminosity‑independent distance estimates and uniform γ‑ray spectral assumptions for cross‑nova comparisons. The results reinforce the internal‑shock paradigm as the dominant mechanism for GeV γ‑ray production in classical novae with main‑sequence companions, and they highlight the value of coordinated multi‑wavelength monitoring (optical, γ‑ray, IR, radio, X‑ray) to fully characterize the shock physics in these transient events.