Are carbon deflagration supernovae triggered by dark matter ?

Collisions between stellar remnants and dark matter in the Galactic bulge are frequent, and the kinetic energy of a primordial black hole incident on a white dwarf, if it is all thermalized, will raise the degenerate core’s temperature, by at least a degree in the case of a lunar mass black hole. This is an underestimate in two ways: the specific heat is less than 3k/2 per particle, and the incoming object is accelerated by gravitational focusing. Detailed physical models have recently been made of this triggering event. Present observational data are equivocal as to whether the radial distribution of type Ia supernovae in galaxies follows the starlight in the galaxies, or is more concentrated towards the center, as collisional triggering would suggest. But future samples of millions of supernovae from the Rubin telescope will change that.

💡 Research Summary

The manuscript titled “Are carbon deflagration supernovae triggered by dark matter?” investigates the provocative hypothesis that primordial black holes (PBHs), a candidate component of dark matter, can collide with white dwarfs (WDs) in galactic bulges and ignite carbon‑deflagration Type Ia supernovae (SNe Ia). The authors combine elementary kinetic‑energy arguments, collision‑rate estimates, semi‑analytic galaxy modeling, and existing supernova surveys to assess whether this mechanism can account for a non‑negligible fraction of observed SNe Ia.

Physical mechanism
A PBH of lunar mass (~10⁻⁷ M⊙) falling into a WD releases a kinetic energy that, if fully thermalized, raises the degenerate core temperature by at least 1 K. The authors note two under‑estimates: (i) the specific heat of a degenerate electron gas is lower than the classical 3k/2 per particle, and (ii) gravitational focusing dramatically increases the impact velocity, with a focusing factor g ≈ 1 + (v_esc²/v²) ≈ 4.8 × 10⁵ for typical bulge velocity dispersions. Consequently, the temperature rise could be substantially larger, potentially pushing the core over the carbon‑ignition threshold.

Collision rate calculation
The rate per unit volume is expressed as ρ = n₁ n₂ σ v, where n₁ and n₂ are the number densities of WDs and PBHs, σ the geometric cross‑section, and v the relative velocity. Using a dark‑matter density in the Galactic bulge of 0.01–0.1 M⊙ pc⁻³, a PBH fraction f of the total dark matter, and a 1‑D velocity dispersion of 110 km s⁻¹, the authors obtain a PBH‑WD encounter rate of 13.8–138 f × 10² Myr⁻¹ for a lunar‑mass PBH. Scaling to the cosmic volume of L* galaxies (≈3.9 × 10⁶ Gpc⁻³) yields a global SNe Ia production rate of 0.9–9 × 10⁴ f Gpc⁻³ yr⁻¹. This matches the observed volumetric rate of (2.55 ± 0.12) × 10⁴ yr⁻¹ Gpc⁻³ h⁻³ if f ≈ 0.2, i.e., if roughly 20 % of dark matter is in PBHs of this mass. The authors acknowledge that this fraction sits near current microlensing limits, and that variations in PBH mass, WD mass distribution, and the minimum WD mass that can explode (≈0.7 M⊙) could shift the rate by factors of a few.

Observational tests
The paper examines two data sets: (1) the Dark Energy Survey (DES) and OzDES, which together provide 1,829 spectroscopically confirmed SNe Ia with measured offsets from host‑galaxy centers; and (2) a historical catalog of SNe Ia from the Central Bureau for Astronomical Telegrams (CBAT) spanning 1930–2015. For the DES sample, the authors compute the second radial moment d_DLR and compare it to theoretical expectations for an exponential disk (Freeman 1970) and a de Vaucouleurs r¹⁄⁴ law. The observed distribution lies between the two curves but shows a modest excess in the innermost radial bin, which the authors attribute to possible selection effects (nuclear SNe being missed).

For the historical sample, they convert angular offsets to dimensionless radii by dividing by the D25 isophotal diameter (corrected for (1+z)⁴ surface‑brightness dimming). The resulting histogram follows closely the Freeman disk prediction, with no central excess. A similar analysis using 2MASS K‑band isophotal radii yields consistent results. The discrepancy between the DES and historical samples suggests that current data are insufficient to definitively detect a central concentration of SNe Ia.

Simulation‑based predictions
To model the spatial distribution expected from PBH‑triggered SNe Ia, the authors use the Millennium simulation’s halo catalog and assume an NFW dark‑matter profile with a concentration ratio (virial radius / scale radius) of 12. They select galaxies with stellar masses 10⁹–10¹⁰·⁵ M⊙, exclude satellites, and separate disk‑dominated (bulge mass < disk mass) from bulge‑dominated systems. Integrating the product of stellar and dark‑matter densities along the line of sight yields a surface‑density map that is more centrally peaked than the Freeman or de Vaucouleurs expectations. The plotted curves in Figure 1 (not reproduced here) illustrate this enhanced central concentration.

Critical assessment
The paper’s strength lies in its interdisciplinary approach: it ties a speculative dark‑matter candidate to a concrete astrophysical observable, and it provides a clear roadmap for testing the idea with upcoming surveys. However, several uncertainties limit the robustness of the conclusions:

1. PBH mass function – The analysis assumes a monochromatic lunar‑mass PBH population. Realistic formation scenarios predict extended mass spectra; lighter PBHs would increase the encounter rate but may lack sufficient energy to ignite a WD, while heavier PBHs could swallow the WD without triggering a thermonuclear runaway.

2. White‑dwarf physics – The ignition condition depends sensitively on the WD’s central density, composition, and temperature profile. The simple “1 K temperature rise” estimate ignores the steep temperature dependence of carbon fusion rates and the possibility that a modest heating may be insufficient to overcome the high degeneracy pressure. Detailed hydrodynamic simulations (e.g., Leung et al. 2025) are required to confirm that a PBH impact can indeed trigger a deflagration.

3. Selection effects – Both the DES and historical samples suffer from incompleteness near galaxy nuclei due to bright backgrounds and fiber‑collision constraints. The authors acknowledge this but do not model the bias quantitatively. A forward‑modeling approach, injecting synthetic SNe Ia into real images and recovering them with the survey pipelines, would provide a more reliable correction.

4. Cosmic variance – The rate calculation scales linearly with the assumed PBH fraction f and the dark‑matter density in bulges. Variations in bulge properties across galaxy types, as well as uncertainties in the Milky Way’s bulge mass, could shift the predicted global rate by an order of magnitude.

5. Alternative triggers – Other mechanisms (e.g., double‑degenerate mergers, single‑degenerate accretion, or dynamical interactions in dense stellar clusters) also predict a modest central concentration of SNe Ia. Disentangling a PBH contribution from these backgrounds will require additional observables, such as the presence of a faint, high‑velocity remnant or specific nucleosynthetic signatures.

Future prospects with Rubin (LSST)
The authors correctly point out that the Rubin Observatory’s Legacy Survey of Space and Time will deliver millions of SNe Ia with precise host‑galaxy photometry. By measuring the radial offset distribution for different host morphologies, stellar masses, and redshift bins, LSST can test whether a statistically significant central excess exists beyond what is expected from stellar light alone. Moreover, LSST’s deep imaging will reduce nuclear masking, and its cadence will enable early‑time light‑curve characterization that might reveal subtle differences between PBH‑triggered and conventional SNe Ia (e.g., rise time, color evolution).

Conclusion
The manuscript presents a plausible, though speculative, pathway for dark matter—specifically primordial black holes—to trigger carbon‑deflagration supernovae. The calculated encounter rates can, under optimistic assumptions, reproduce the observed SNe Ia volumetric rate, and the predicted spatial distribution is more centrally concentrated than stellar light. Current observational data are inconclusive, but the upcoming Rubin LSST survey offers a realistic avenue to confirm or refute the hypothesis. Further work should focus on (i) detailed hydrodynamic simulations of PBH‑WD impacts, (ii) rigorous modeling of selection biases in existing SN surveys, and (iii) exploration of complementary signatures (e.g., remnants, nucleosynthesis) that could uniquely identify a PBH‑triggered event.