{"id":"7c9d59ae-ebde-4592-ac8e-8c6675853cf6","arxiv_id":"2505.21256","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Collision rate simulations between asteroid-mass primordial black holes and white dwarfs reproduce the observed rate, host, and brightness distributions of type Ia supernovae, and the author argues these dark matter impacts are favored over binary companion ignitions.","lead":"A single-author simulation study argues that asteroid-mass black holes, which could make up dark matter, crash into white dwarf stars often enough to set off type Ia supernovae, the explosions used as cosmic distance beacons. The author claims these dark matter impacts reproduce observed supernova rates and brightnesses and beat the standard binary-companion picture in a statistical comparison.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim inherits an unverified prompt-detonation ignition cross section Rm from the author's prior PRL; if the contested ignition mechanism is wrong, the fitted PBH parameters and the 21-sigma preference lose meaning.","rationale":"The paper is a serious end-to-end forward-modeling effort: it propagates WD masses, SFHs, halo profiles, and PBH mass functions through a collision-rate equation and produces six rate distributions in one framework. The agreement with the shape of the ZTF NMD, the DTD, and the host-offset distributions is genuinely interesting, and the comparison exposes a real weakness of the double-degenerate NMD. However, the headline claim is only as strong as the ignition physics in eq. (2). That physics is not independently established in this Letter; it is carried over from PRL [15], and the contested status of the prompt-detonation mechanism is acknowledged only obliquely via [14] and [21]. The 21-sigma language is also statistically nonstandard (converting AIC evidence ratios to Gaussian significances is not a valid hypothesis test), and the absolute reduced chi-square of the PBH fits is poor (about 4 to 9), but those issues could be repaired by reframing the comparison. The ignition cross section cannot be repaired by reframing: if the collision does not produce a prompt detonation, the scenario is simply not viable for normal SNe Ia. I therefore agree with the reader's weakest-assumption identification. The appropriate verdict remains CONDITIONAL—the modeling is worth publishing, but the ignition mechanism must be independently verified or clearly presented as an assumption, and the sigma-language should be replaced by a standard model-comparison statistic. No change to the reader's verdict is needed.","tokens_in":27762,"tokens_out":5279,"duration_ms":62871,"concrete_test":"Run an independent 3D hydrodynamical simulation (e.g., Castro or FLASH) of a PBH transiting a 1.0 solar-mass CO white dwarf with XC = 0.3 for log m_bh = 21, using the same nuclear network and initial conditions as [15]; measure whether a self-sustained supersonic detonation wave forms at the radius Rm predicted by eq. (2). Compare the outcome with the deflagration-quenching result of Montero-Camacho et al. [14]. If no prompt detonation occurs, the rate distributions, the fitted PBH parameters, and the claimed preference over the DD scenario are invalidated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that PBH ignition is favored over double-degenerate ignition at >21 sigma and that the six observed rate distributions are consistent—flows through eq. (2), where the per-WD collision rate Gamma_Ia is proportional to the ignition cross-section radius Rm(mw, m_bh, XC) inherited from the author's prior PRL [15]. That cross section is not independently verified in this Letter: it assumes a quasi-nuclear-statistical-equilibrium velocity threshold, an ignition efficiency alpha=5 (justified in the Supplemental Material only by matching the deflagration ignition cross section in the over-satisfied regime), and a CO-to-ONe composition transition. The mechanism is contested: Montero-Camacho et al. [14] argued that the burning is a quenched subsonic deflagration, and the paper itself cites Kushnir et al. [21] as tension with the observed t0-MNi relation. The Supplemental Material verifies only the induction-length scaling xi_perp proportional to XC^-2 (Fig. S4), not the existence of a self-sustained detonation. If Rm is wrong, every predicted distribution, the fitted log mc and sigma_c in eq. (15), and the AIC comparison lose meaning. This is the single most load-bearing assumption: all six distributions and the headline significance pass through it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that collisions between white dwarfs (WDs) and asteroid-mass primordial black holes (PBHs), which trigger prompt detonations, can explain the observed rate distributions of normal type Ia supernovae (SNe Ia). The author builds a forward model that convolves WD populations, galactic dark-matter densities, and PBH mass functions to predict the 56Ni-mass distribution, cluster delay-time distribution, volumetric rate history, host-offset distributions, and host-mass distribution. These predictions are compared with literature data for eight combinations of modeling assumptions, and the PBH scenario is compared with a double-degenerate (DD) scenario using chi-square minimization and the Akaike information criterion. The paper claims that the PBH rate distributions are consistent with observations and that PBH ignition is favored over DD ignition at more than 21 sigma confidence, with a best-fit log-normal PBH mass function peaking at log mc = 21.13 +/- 0.12 g and width sigma_c = 0.40 +/- 0.20.","tokens_in":27887,"tokens_out":5751,"duration_ms":65533,"significance":"If the central claims were fully supported, the paper would be highly significant: it would identify a dark-matter mechanism for SN Ia ignition, populate the asteroid-mass PBH window, and resolve a long-standing progenitor problem. The modeling framework is ambitious and physically transparent: it covers WD mass, host stellar mass, galactocentric offset, and cosmic time in a single simulation, and it tests eight systematic variations of star formation history, IMF, IFMR, Iax contribution, and 56Ni dispersion. The use of a wide set of literature rate distributions as benchmarks is a strength. However, the significance is conditional: the two headline claims rest on an unverified ignition cross-section from the author's prior work, and the absolute quality of the fits is poor even for the best PBH models. The paper is therefore better viewed as a proof-of-concept for a PBH-trigger scenario than as a demonstrated resolution of the SN Ia progenitor problem.","major_comments":[{"comment":"The claim in the Summary and discussion that 'rate distributions are consistent with observations' is not supported by the reported fits. For the restricted sample, the best PBH model (PBH2d) has total chi-square 204.1 for 47 bins, i.e. reduced chi-square about 4.4 after accounting for two fitted parameters. Per-distribution values are also poor: HOD contributes 32.2 for 8 bins, HND 15.1 for 5 bins, and HMD 43.6 for 13 bins. The AIC comparison therefore shows only that the PBH model is less bad than the DD model, not that the model is an adequate description of the data. The wording 'consistent' should be replaced by a quantitative statement of goodness-of-fit, and the main result should be reframed accordingly.","section":"Supplemental Material, Table SII"},{"comment":"The entire rate prediction and the fitted parameters (log mc, sigma_c in Eq. 15) flow through the ignition cross-section radius Rm(mw, m_bh, XC) inherited from the author's prior PRL [15]. The Supplemental Material verifies only the scaling xi_perp proportional to XC^-2 (Fig. S4) and justifies alpha = 5 by matching the deflagration ignition cross-section (Fig. S2); it does not demonstrate that a self-sustained prompt detonation actually occurs. The mechanism is contested by Montero-Camacho et al. [14], and the paper itself cites Kushnir et al. [21] for tension with the observed t0-MNi56 relation. If Rm is incorrect, the fitted parameters, the predicted distributions, and the 21-sigma preference all lose their quantitative meaning. The paper should either provide independent verification of the cross-section or explicitly present all results as conditional on this unverified assumption, with a sensitivity study treating detonation and deflagration as separate cases.","section":"Eq. (2) and Supplemental Material SI"},{"comment":"The abstract's claim that the brightness distribution 'comes out without finetuning' is not supported, because the 56Ni-mass distribution is one of the six distributions entering the joint chi-square in Eq. (S7), and the two PBH parameters (log mc, sigma_c) are fitted to all six distributions simultaneously. The NMD is therefore a fitted output, not an a priori prediction. A genuine test would fix the PBH parameters using, for example, the VRH, DTD, HOD, HND, and HMD data and then compare the predicted NMD. In addition, the mNi(mw) relation used in Eq. (5) is itself a fit to hydrodynamical simulation results (Fig. S5), adding a further fitted ingredient. The 'without finetuning' phrasing should be removed or explicitly qualified.","section":"Abstract and Eq. (S7)"},{"comment":"The conversion from AIC differences to 'more than 21 sigma' via p = exp(Delta AIC / 2) conflates relative model comparison with absolute goodness of fit. With best-fit reduced chi-square near 4, the very large Delta AIC is driven primarily by the DD model being catastrophically bad, not by the PBH model being statistically acceptable. The reported significance also does not incorporate systematic uncertainties in the data binning, survey selection, or the modeling assumptions, despite the large spread of chi-square values across the eight models in Table SII. The authors should report absolute goodness-of-fit statistics and use a model-comparison framework that accounts for model inadequacy, or at minimum state clearly that the significance is a relative likelihood under the assumed models.","section":"Parameter fitting and scenario comparison"}],"minor_comments":[{"comment":"In the cosmology line, 'H0 = 70 km s^-1 Mpc^-3' should read 'H0 = 70 km s^-1 Mpc^-1'.","section":"Simulation setup"},{"comment":"The phrase 'the evolution equation maybe setup as' should be 'the evolution equation may be set up as'.","section":"Simulation setup"},{"comment":"The legend label 'Friedman+18' is inconsistent with reference [6], which is 'Friedmann and Maoz'; please standardize.","section":"Fig. 1b"},{"comment":"In the HOD section, 'Palomer Transient Facility' should be 'Palomar Transient Facility'.","section":"Supplemental Material, Data analysis"},{"comment":"The caption contains the typo 'white dwar mass'; it should be 'white dwarf mass'.","section":"Supplemental Material, Fig. S5 caption"},{"comment":"The funding agency name appears as 'F APES/CAPES'; this should be 'FAPES/CAPES'.","section":"Acknowledgments"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's headline claims are substantially stronger than the internal evidence supports: the absolute fits are poor (Table SII), and the central mechanism is inherited from the author's own prior PRL [15] without independent verification. The accompanying companion paper [23] supplies the SFH framework, so the new contribution here is the rate-distribution comparison, but that comparison is currently presented as a definitive resolution rather than a conditional exploration. If the editor wishes to maintain a high evidentiary bar for claims about SN Ia progenitors, the paper needs major revision with reframed claims; in its present form the '21-sigma' and 'without finetuning' statements are disproportionate to the reported chi-square values."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the end-to-end forward model: collision rates between WDs and PBHs, integrated over WD mass, host stellar mass, galactocentric radius, and cosmic time, giving six observed rate distributions. That is a real step beyond the order-of-magnitude estimate in the author's earlier PRL. The qualitative simultaneity—NMD, DTD, VRH, HOD, HND, and HMD all coming out roughly right, while the DD model visibly fails on the nickel mass distribution—is worth taking seriously.\n\nBut the paper overclaims. The \"more than 21 sigma\" number comes from converting an AIC difference into Gaussian significance via p = exp(-ΔAIC/2). That is a relative likelihood, not a p-value, and the conversion to sigma is not a valid hypothesis test. Worse, the comparison is between two bad models: the PBH best fit has reduced chi-square of roughly 4–9, so both scenarios are formally excluded by the data. Saying one is preferred at 21 sigma when neither is a good fit is misleading. Also, the NMD advertised as \"without finetuning\" is one of the six fitted targets in eq. (S7); its agreement is not a prediction.\n\nThe load-bearing piece, the ignition cross section Rm, is inherited from the author's contested PRL, and the Supplemental Material only checks the induction-length scaling, not the existence of a self-sustained detonation. If Montero-Camacho et al. are right that the burning is a quenched deflagration, the fitted PBH parameters and the AIC comparison lose their meaning. The discussion acknowledges this, but the abstract does not.\n\nThe restricted-sample cuts (removing offsets beyond 20 kpc and NMD edge bins) are justified, but they make \"consistency\" relative and remove exactly the regions where the model fails.\n\nWhat I would want before accepting the central claim: hold out one distribution as a genuine prediction, give the DD scenario a free normalization (it gets zero free parameters while PBH gets two), replace the sigma-language with a standard model comparison or likelihood ratio test, release code and data, and engage the ignition-mechanism criticism head-on.\n\nThis is a serious paper that deserves a serious referee. The framework is coherent, the systematic coverage is careful, and the qualitative DD problem with the nickel mass distribution is a useful result. But the headline significance should be scaled back, and the fitted-to-predicted distinction needs to be clear.","headline":"Real forward-modeling step, but the 21-sigma headline is not supported by the paper's own statistics.","tokens_in":28624,"tokens_out":2940,"would_cite":false,"duration_ms":38402,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that collisions between white dwarfs and asteroid-mass primordial black holes reproduce all major rate distributions of normal type Ia supernovae, and that this dark-matter ignition channel is favored over binary double…","keywords":["type Ia supernovae","primordial black holes","white dwarf detonation","dark matter","nickel-56 mass distribution","delay time distribution","asteroid-mass black holes","SN Ia progenitors"],"falsifier":"Run a multidimensional hydrodynamical simulation of a $\\sim 10^{21}$ g black hole crossing a $\\sim 1\\,M_\\odot$ carbon-oxygen white dwarf and resolve the shocked region: a quenched deflagration would invalidate the central mechanism. Independently, a gamma-ray burst lensing parallax survey sensitive to $\\sim 10^{21}$ g objects would directly test whether the inferred dark-matter population exists.","tokens_in":27356,"feed_emoji":"💥","tokens_out":9702,"duration_ms":92132,"temperature":0.7,"pith_summary":"This paper argues that normal type Ia supernovae are ignited not by a binary companion but by the passage of an asteroid-mass primordial black hole through a white dwarf. The author builds collision-rate simulations spanning white-dwarf masses, host galaxy stellar masses, galactocentric offsets, and cosmic time, and compares the predicted rate distributions with literature observations. The central result is that a two-parameter log-normal black-hole mass function centered near $\\log m_c = 21.13$ (with $m_c$ in grams) reproduces the observed nickel-56 mass distribution, delay-time distribution, volumetric rate history, host-offset distributions, and host-mass distribution simultaneously. The same data disfavor the standard double-detonation binary scenario at more than $21\\sigma$ under all eight modeling-assumption combinations tested. If the paper is right, the long-standing progenitor problem is resolved by making SN Ia ignition a dark-matter phenomenon, at the price of a black-hole mass scale that is not predicted from first principles.","feed_headline":"Asteroid-mass black holes may ignite type Ia supernovae","feed_subtitle":"One dark-matter mass function fits nickel yields, delay times, cosmic rates, and host offsets together.","key_machinery":"The load-bearing object is the primordial-black-hole--white-dwarf prompt-detonation ignition cross section, a radius $R_m(m_w, m_\\bullet, X_C)$ inside which the shock from the infalling compact object exceeds the quasi-nuclear statistical equilibrium detonation velocity of degenerate carbon-oxygen matter. This cross section enters a rate equation for the number of white dwarfs per mass and galactocentric shell, so the collision rate scales with the local dark-matter density and with the log-normal black-hole mass function $\\phi(m_\\bullet)$. The same rate machinery, using a literature double-detonation delay-time distribution, generates the comparison scenario. The ignition efficiency is fixed at $\\alpha = 5$ by matching the deflagration ignition cross section in massive white dwarfs.","core_discovery":"The paper claims that the long-unexplained distribution of peak brightnesses of normal type Ia supernovae, together with their delay-time distribution, cosmic rate history, host-offset distributions, and host-mass distribution, can all be reproduced by one mechanism: a white dwarf in a dark-matter halo is struck by a compact object of roughly $10^{21}$ grams, and the resulting shock crosses the detonation threshold so the star burns supersonically. The mass function that does this work is log-normal, with $\\log m_c = 21.13 \\pm 0.12\\,(\\mathrm{sys}) \\pm 0.05\\,(\\mathrm{stat})$ and width $\\sigma_c = 0.40 \\pm 0.20\\,(\\mathrm{sys}) \\pm 0.10\\,(\\mathrm{stat})$. The same machinery, fed with the literature double-detonation delay-time distribution, fits the observations far worse; across eight combinations of star-formation history, initial mass function, initial-final mass relation, Iax contamination, and nickel-mass dispersion, the primordial-black-hole scenario is preferred at more than $21\\sigma$. The paper identifies its main result as the consistency of the rate distributions with observations, with the only severe contradiction being that the inferred black-hole mass scale is unpredicted from first principles.","pith_inferences":["Not in the paper: if the fitted mass function is taken at face value, SN Ia catalogs could be cross-correlated with dark-matter halos to map halo density profiles, because the host-offset distributions computed here are exactly the observable that carries that signal.","Not in the paper: the paper's interpretation of the low-mass tail of the nickel distribution as lighter white dwarfs needing heavier black holes predicts a measurable anti-correlation between SN Ia brightness and local dark-matter density that volumetric samples could test.","Not in the paper: the best-fit width $\\sigma_c = 0.40$ brackets the critical-collapse value 0.26, so a sharper measurement of the nickel-distribution shape could distinguish inflationary log-normal production from critical collapse.","Not in the paper: since the inferred peak mass shifts by 0.18 dex per unit change in the ignition efficiency $\\alpha$, an independent multidimensional hydrodynamical determination of $\\alpha$ would either tighten or dissolve the claimed statistical preference."],"forward_implications":["If the paper is right, normal type Ia supernovae need no binary companion at all, explaining the absence of companion ejecta interaction signatures in observations.","Supernova rate surveys become a probe of the asteroid-mass primordial black hole mass function and of dark-matter density inside galaxies, because the predicted rate is proportional to the local dark-matter density.","The width-luminosity relation and the skewed nickel-56 mass distribution follow from white-dwarf mass and ignition-threshold physics rather than from binary population statistics, removing the need for fine tuning.","All known binary ignition channels, represented here by the double-detonation scenario, are disfavored by the rate-distribution data at more than $21\\sigma$ across the modeling assumptions tested."],"supporting_citations":[{"why":"supplies the semi-analytical prompt-detonation ignition cross section that drives the collision rate in eq. (2).","marker":"[15]"},{"why":"provides the ZTF-BTS nickel-56 mass distribution used as the main brightness benchmark.","marker":"[10]"},{"why":"provides cluster supernova rates used for the delay-time distribution comparison.","marker":"[6]"},{"why":"compiles the volumetric rate history data over cosmic time.","marker":"[41]"},{"why":"supplies the ZTF host-offset and host-normalized offset distributions.","marker":"[39]"},{"why":"gives the host stellar-mass distribution from DES, PTF, and PS1.","marker":"[42]"},{"why":"provides the double-detonation delay-time distribution used for the DD scenario.","marker":"[8]"},{"why":"supplies the double-detonation nickel-mass distribution fit that the comparison disfavors.","marker":"[12]"},{"why":"argues the burning is a quenched subsonic deflagration, the competing ignition picture the prompt-detonation cross section must exclude.","marker":"[14]"},{"why":"models galaxy star formation histories in radial zones, needed because encounter rates scale with dark-matter density.","marker":"[23]"}],"fun_headline_variants":["Primordial black holes spark supernova brightness puzzle","Asteroid-mass black holes trigger type Ia supernovae","Dark matter black holes could solve supernova ignition","One PBH mass fits all supernova rate distributions","PBH encounters explain type Ia supernova brightness"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire rate prediction inherits the prompt-detonation ignition cross section from the author's earlier work, including the adopted ignition efficiency $\\alpha = 5$; if the burning actually stalls as a quenched subsonic deflagration rather than detonating supersonically, the fitted black-hole masses and the $21\\sigma$ preference lose their basis.","fun_headline_variants_meta":{"raw":{"variants":["Primordial black holes spark supernova brightness puzzle","Asteroid-mass black holes trigger type Ia supernovae","Dark matter black holes could solve supernova ignition","One PBH mass fits all supernova rate distributions","PBH encounters explain type Ia supernova brightness"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000211,"raw_usage":{"total_tokens":1409,"prompt_tokens":932,"completion_tokens":477,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":548,"completion_tokens_details":{"reasoning_tokens":402}},"tokens_in":548,"tokens_out":477,"duration_ms":5495,"temperature":1.0,"reasoning_tokens":402,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T13:35:22.869218+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a multidimensional hydrodynamical simulation of a $\\sim 10^{21}$ g black hole crossing a $\\sim 1\\,M_\\odot$ carbon-oxygen white dwarf and resolve the shocked region: a quenched deflagration would invalidate the central mechanism. Independently, a gamma-ray burst lensing parallax survey sensitive to $\\sim 10^{21}$ g objects would directly test whether the inferred dark-matter population exists.","supporting_citations":[{"cited_title":"The ZTF-BTS Type Ia supernovae luminosity function is consistent with a single progenitor channel for the explosions","cited_arxiv_id":"2109.06219","evidence_quote":"provides the ZTF-BTS nickel-56 mass distribution used as the main brightness benchmark."},{"cited_title":"The rate of Type-Ia supernovae in galaxy clusters and the delay-time distribution out to redshift 1.75","cited_arxiv_id":"1803.04421","evidence_quote":"provides cluster supernova rates used for the delay-time distribution comparison."},{"cited_title":"Delay Time Distributions of Type Ia Supernovae From Galaxy and Cosmic Star Formation Histories","cited_arxiv_id":"2001.05967","evidence_quote":"compiles the volumetric rate history data over cosmic time."},{"cited_title":"The Zwicky Transient Facility Census of the Local Universe I: Systematic search for Calcium rich gap transients reveal three related spectroscopic sub-classes","cited_arxiv_id":"2004.09029","evidence_quote":"supplies the ZTF host-offset and host-normalized offset distributions."},{"cited_title":"Projected distances to host galaxy reduce SNIa dispersion","cited_arxiv_id":"1612.04417","evidence_quote":"gives the host stellar-mass distribution from DES, PTF, and PS1."},{"cited_title":"Delay Times and Rates for Type Ia Supernovae and Thermonuclear Explosions from Double-detonation Sub-Chandrasekhar Mass Models","cited_arxiv_id":"1011.1407","evidence_quote":"provides the double-detonation delay-time distribution used for the DD scenario."},{"cited_title":"Confronting double-detonation sub-Chandrasekhar models with the low-luminosity suppression of Type Ia supernovae","cited_arxiv_id":"2110.13169","evidence_quote":"supplies the double-detonation nickel-mass distribution fit that the comparison disfavors."},{"cited_title":"Dark-matter-induced transients over cosmic time: The role of star formation history profiles","cited_arxiv_id":"2505.21260","evidence_quote":"models galaxy star formation histories in radial zones, needed because encounter rates scale with dark-matter density."}],"review_version":1}