{"id":"a08cb3b2-2433-49c1-b8b4-03b23767720e","arxiv_id":"2504.21129","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"From 68 fast radio bursts and 1048 supernovae, the authors find a photon rest mass of roughly 18 to 29 times 10^-51 kilograms depending on the assumed intergalactic gas fraction, with zero mass still allowed at 2 to 3 sigma.","lead":"Astronomers used 68 fast radio bursts and supernova data to place limits on whether the photon has a tiny mass, without assuming a particular model of cosmic expansion. The best fit is around 10^-50 kilograms, but the result depends strongly on assumptions about gas filling intergalactic space.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model-independent rewrite of DM_IGM and DM_γ (Eqs. 4.5, 4.6) silently assumes flat FLRW via d_L=(1+z)c∫dz/H; nonzero spatial curvature changes the integrals and can bias m_γ.","rationale":"The reader's weakest_assumption matches my own scan of the paper. The derivation of Eqs. (4.5) and (4.6) is the bridge that turns raw FRB DMs and SNe d_L data into constraints on m_γ; if that bridge is only valid in flat space, the headline 'cosmological-model independent' is not supported and the numerical limits are conditional on k=0. I checked the integration by parts and confirmed that both equations use d_L/(1+z) = c∫dz'/H, which is the flat-space distance relation. Other candidate concerns (post hoc sample pruning, fixed H_0/Ω_b/M_B, DM_halo = 50 pc/cm^3) are real but either are explicitly documented, affect nuisance parameters rather than the geometric kernel, or are secondary to the derivation. The paper deserves credit for not claiming a detection, for treating f_IGM both fixed and free, and for a transparent data list; none of these rescues the unstated flatness premise. The proposed test quantifies whether the premise matters numerically; if it does not, the verdict can be upgraded to acceptance with a clarifying assumption statement.","tokens_in":15935,"tokens_out":16121,"duration_ms":175534,"concrete_test":"Re-derive Eqs. (4.5) and (4.6) in a curved FLRW metric, then for each of the 68 FRB redshifts compute the exact DM_IGM and DM_γ values for Ω_K = ±0.01 and ±0.1 (holding H_0 = 74.03 and Ω_b h^2 fixed as in the paper) and compare them with the flat-formula outputs. If the fractional differences exceed the fractional statistical uncertainties in DM_IGM/DM_γ, re-run the MCMC with Ω_K added as a free parameter (or marginalized over a Planck/CMB prior) and report the shift in m_γ; if the shift is within the 1σ width, the flatness assumption is not numerically load-bearing, but it should still be stated explicitly.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the replacement of the DM_IGM and DM_γ path integrals by luminosity-distance terms in Eqs. (4.5) and (4.6). Both identities use the flat-FLRW relation d_L(z) = (1+z)c ∫_0^z dz'/H(z'). For DM_IGM, Eq. (4.5) follows from writing d_L/(1+z) = c∫dz'/H and integrating by parts; for DM_γ, Eq. (4.6) follows from d_L/(1+z)^3 = c/(1+z)^2 ∫dz'/H. In a curved FLRW metric, d_L/(1+z) is the transverse comoving distance D_M(z) = S_k(χ), whose derivative is not c/H, so both expressions acquire Ω_k-dependent correction factors. The manuscript never states this geometric assumption; the headline 'cosmological-model independent' claim therefore rests on an unstated flatness premise. Since DM_γ enters as m_γ^2 and is degenerate with f_IGM,0 and DM_host,0, a curvature correction shifts the inferred m_γ and also modifies the uncertainty propagation through Eq. (4.2). This is a correctness risk in the derivation itself, not a mere labeling issue: the central limit is obtained from integrals that are proven only for k=0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper constrains the photon rest mass mγ by combining 68 well-localized fast radio bursts with luminosity distances reconstructed from the 1048-object Pantheon supernova sample. Starting from the massive-photon dispersion relation, the authors define a photon-mass contribution DMγ(z) to the dispersion measure and rewrite both DM_IGM and DM_γ in Eqs. (4.5)–(4.6) as combinations of the luminosity distance and redshift integrals, avoiding any assumed H(z). An MCMC is then run for two scenarios: f_IGM fixed at 0.83 and f_IGM free. The free case yields mγ = (29.4^{+5.8}_{−15.5}) × 10^{-51} kg, f_IGM,0 = 0.902^{+0.034}_{−0.631}, and DM_host,0 = 90^{+19}_{−22} pc/cm^3 at 1σ, together with a strong anticorrelation between f_IGM and mγ. The fixed case gives mγ = (18.2^{+2.7}_{−5.9}) × 10^{-51} kg.","tokens_in":16209,"tokens_out":9985,"duration_ms":103910,"significance":"The intended contribution is a photon-mass limit that avoids specifying the dark-energy expansion history and therefore avoids the ΛCDM circularity present in earlier FRB-based limits. The paper's algebraic transformation of the DM integrals into d_L-based terms is transparent and, under a flat FLRW geometry, internally consistent. It also makes efficient use of the growing sample of localized FRBs and tabulates the individual events. The main limitation is that the 'cosmological-model independent' phrasing overstates the result: the key identities assume flat spatial geometry, and the error budget omits some external uncertainties. With these points addressed, the method is a useful and falsifiable approach to photon-mass constraints.","major_comments":[{"comment":"The derivation of DM_IGM and DM_γ is done by integration by parts using the flat-FLRW relation d_L(z) = (1+z)c∫_0^z dz'/H(z'). For nonvanishing spatial curvature, d_L/(1+z) is the transverse comoving distance D_M(z) = S_k(χ), whose derivative is not c/H, so the identities in Eqs. (4.5) and (4.6) acquire Ω_k-dependent correction factors. The paper never states that flatness is assumed. This matters because DMγ enters as mγ² and is degenerate with f_IGM,0 and DM_host,0 through Eq. (3.12), so a curvature correction shifts the inferred mγ and its uncertainty. Please state the flatness assumption explicitly, or generalize the derivation to S_k, or quantify the sensitivity to Ω_k; as written, the 'cosmological-model independent' claim rests on an unstated geometric premise.","section":"§4.3, Eqs. (4.5)–(4.6)"},{"comment":"The quoted external uncertainties H0 = 74.03 ± 1.4 km/s/Mpc and Ω_b h² = 0.02235 ± 0.00037 are not propagated into σ_IGM or into the MCMC likelihood. The prefactor A = 3cΩ_b H0²/(8πG m_p) in Eq. (4.5) contains H0² and Ω_b, so its fractional uncertainty is roughly 4%, yet Eq. (4.2) includes only σ_dL and σ_I. Because DM_IGM is proportional to A f_IGM,0 and mγ is found to be anticorrelated with f_IGM,0, the omission can underestimate the error bars on mγ. Please propagate the H0 and Ω_b uncertainties, or show explicitly that their effect is negligible for the reported 1σ intervals.","section":"§4.3, Eqs. (4.1)–(4.2)"},{"comment":"The uncertainty σ_I of the integral in Eq. (4.6) is not defined. The integral is computed from GP-reconstructed d_L values at a set of redshifts, and those values are strongly correlated in the GP posterior. If σ_I is computed from the diagonal errors only, the likelihood overstates the information content of the data and the reported uncertainties are not reliable. Please specify how σ_I is calculated, including whether the full GP covariance matrix is propagated.","section":"§4.3, Eq. (4.2)"}],"minor_comments":[{"comment":"The phrase 'a discrepancy value between them at 11.2×10^-31 kg' should read 10^-51 kg, since the constraints are quoted in units of 10^-51 kg.","section":"§5, first paragraph"},{"comment":"The likelihood function and the priors used in the emcee runs are not stated; please include the exact likelihood expression and the prior ranges or distributions for mγ, f_IGM,0, and DM_host,0 for reproducibility.","section":"§4.3"},{"comment":"The Pantheon systematic covariance matrix is not used in the GP reconstruction of d_L; the authors should state whether this is intentional and whether the supernova systematics are reflected in σ_dL.","section":"§4.2"},{"comment":"The fixed 30 pc/cm^3 host-galaxy uncertainty enters Eq. (4.1) in quadrature, but the treatment of the host-galaxy scatter and its relation to the fitted DM_host,0 is not discussed; a brief justification would help.","section":"§4.1"},{"comment":"The comparison with Ref. [45] mentions sample size and number of free parameters, but not the different treatment of f_IGM; adding a sentence noting that [45] fixes f_IGM while the present work fits it would make the comparison clearer.","section":"§6"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the flatness issue is the main substantive concern; it is correctable, but it touches the paper's central claim of model independence, so I recommend major revision rather than outright acceptance. I do not see grounds for rejection: the derivation is transparent, the data selection is documented, and the method can be fixed by either an explicit flatness caveat or a curvature-generalized treatment."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things. This is the same group's model-independent FRB+SNe machinery, now pointed at the photon mass; the new pieces are the luminosity-distance form of the massive-photon DM integral (Eq. 4.6) and a constraint from 68 FRBs plus Pantheon. The result is modest: the free-f_IGM case gives m_γ = (29.4^{+5.8}_{-15.5})×10^{-51} kg at 1σ, consistent with zero at ~2σ, and about ten times weaker than the existing cosmology-independent limit from Ran, Wang & Wei (≤3.5×10^{-51} kg). That is stated plainly.\n\nThe authors do a lot right. The data cuts are itemized with reasons, the f_IGM–m_γ degeneracy is acknowledged and displayed, and there is no detection overclaim. The method does avoid the ΛCDM circularity it criticizes: the photon mass is fitted, not assumed. The citation pattern is healthy—closely related work is cited and compared, and the self-reference to [27] is legitimate since this is their own machinery.\n\nSoft spots, in order of how much they matter. The flatness assumption: the stress-test is right that Eqs. 4.5–4.6 use d_L = (1+z)c∫dz/H, the flat-FLRW identity, and the paper never says so. The practical bias is small—with |Ω_k| ≲ 0.02 the curvature correction is of order Ω_k (∫_0^z H_0/H dz')^2/6, around 10^{-3}—so it will not change the answer, but the label 'cosmological-model independent' should be qualified to 'flat', and a sentence plus a curvature robustness run would settle it.\n\nThen three smaller things. Eq. 4.6 as printed has a units problem in the first term: it should be d_L/(c(1+z)^3) to match the (2/c)∫ term; likely a typo, but it needs fixing. Eq. 4.5 drops the χ = 7/8 factor that Eq. 3.14 defines; the fitted f_IGM,0 then absorbs that factor, which changes the comparison with external baryon-fraction estimates and can shift the fixed-case m_γ by ~10%. And the MCMC priors are not reported, while H0 and Ω_b enter with quoted errors that are not marginalized over; the first is a reproducibility issue, the second is an easy improvement.\n\nThe single global DM_host,0 for all 68 FRBs is crude, and the pruning from 97 to 68 events is justified but deserves a robustness check. The fixed-f_IGM case actually excludes m_γ = 0 at ~3σ; given the missing 7/8 factor, that tension should be examined before anyone quotes it as a hint.\n\nBottom line: a publishable-after-revision paper, not a breakthrough. It advances the program incrementally, is honestly presented, and the central method survives the flatness critique numerically. If you work on FRB dispersion or photon-mass bounds, it is worth a read as a cross-check; if not, it can wait. It deserves a serious referee—the derivation issue is real but fixable, and the transparency makes referee work tractable. My recommendation: send it to peer review, and ask for the flatness statement, the two equation fixes, and the robustness checks.","headline":"Honest, incremental extension of the authors' own model-independent FRB+SNe method to photon mass, with an unstated flatness assumption behind Eq. 4.6 that is real but numerically minor; the paper deserves a serious referee.","tokens_in":16805,"tokens_out":26095,"would_cite":false,"duration_ms":235439,"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":"The paper shows that photon rest mass can be constrained without assuming a cosmological model, by using fast radio burst dispersion measures together with supernova distances.","keywords":["photon rest mass","fast radio bursts","dispersion measure","type Ia supernovae","intergalactic medium","model-independent cosmology","Pantheon","baryon fraction"],"falsifier":"A single well-localized FRB with an independently measured luminosity distance (for example from a gravitational-wave standard siren or a Cepheid-calibrated host) would settle the flatness question: if the DM_IGM and DM_gamma residuals computed with the true d_L disagree with the fitted m_gamma at more than the quoted 1-sigma level, the massive-photon interpretation fails. More directly, a future sample of hundreds of FRBs should show residuals around the model that are statistically consistent with the fitted m_gamma; a significant excess that correlates with redshift and not with host properties would falsify the model.","tokens_in":15694,"feed_emoji":"🔭","tokens_out":7271,"duration_ms":69694,"temperature":0.7,"pith_summary":"This paper tries to bound the rest mass of the photon without committing to a cosmological model. It combines the dispersion measures of 68 well-localized fast radio bursts with luminosity distances of 1048 type Ia supernovae from the Pantheon catalog, rewriting both the intergalactic-medium and the massive-photon dispersion integrals as functions of the supernova luminosity distance. In the case where the IGM baryon fraction is treated as free, the photon mass is found to be $m_\\gamma = (29.4^{+5.8}_{-15.5}) \\times 10^{-51}$ kg at $1\\sigma$, with $f_{\\mathrm{IGM,0}} = 0.902^{+0.034}_{-0.631}$. The result is consistent with a zero photon mass at about $2\\sigma$, and it shows a strong anticorrelation between $m_\\gamma$ and $f_{\\mathrm{IGM}}$, which the authors read as a warning that baryon-fraction assumptions dominate the constraint.","feed_headline":"68 FRBs put photon mass near 3e-50 kg","feed_subtitle":"A model-independent match of 68 FRBs to 1048 supernovae hints at a nonzero photon mass while flagging a baryon-fraction degeneracy.","key_machinery":"The central object is the decomposition of the observed extragalactic dispersion measure into host, IGM, and photon-mass terms, $\\mathrm{DM}_{\\mathrm{ext}}^{\\mathrm{th}}(z) = \\mathrm{DM}_{\\mathrm{host}}(z) + \\mathrm{DM}_{\\mathrm{IGM}}(z) + \\mathrm{DM}_\\gamma(z)$. The load-bearing move is the integration-by-parts identities $$\\mathrm{DM}_{\\mathrm{IGM}}(z) = A\\,f_{\\mathrm{IGM,0}}\\left[ \\frac{d_L(z)}{c} - \\frac{1}{c}\\int_0^z \\frac{d_L(z')}{1+z'}\\,dz' \\right]$$ and the displayed expression for $\\mathrm{DM}_\\gamma(z)$, which convert the unknown Hubble-parameter integrals into luminosity distances that the Pantheon SNe supply through a Gaussian-process reconstruction. All the cosmology is thereby concentrated in $d_L(z)$, so the comparison of theory to the measured DM of 68 FRBs constrains $m_\\gamma$, $f_{\\mathrm{IGM,0}}$, and $\\mathrm{DM}_{\\mathrm{host,0}}$ without specifying a dark-energy model.","core_discovery":"The paper's central claim is that the photon rest mass can be constrained from FRB dispersion measures without assuming $\\Lambda$CDM by expressing the massive-photon contribution as $$\\mathrm{DM}_\\gamma(z) = B\\,m_\\$gamma^{2}$\\left[ \\frac{d_L(z)}{(1+z)^3} + \\frac{2}{c}\\int_0^z \\frac{d_L(z')}{(1+z')^4}\\,dz' \\right],$$ where $d_L$ comes directly from SNe data. This expression follows from integration by parts on the massive-photon dispersion integral using the luminosity-distance relation, avoiding the Hubble-parameter integrals that would require a cosmological model. With 68 well-localized FRBs and Pantheon SNe, the analysis yields $m_\\gamma = (29.4^{+5.8}_{-15.5}) \\times 10^{-51}$ kg ($1\\sigma$) when $f_{\\mathrm{IGM}}$ is free, and $m_\\gamma = (18.2^{+2.7}_{-5.9}) \\times 10^{-51}$ kg when $f_{\\mathrm{IGM}} = 0.83$ is fixed; the host-galaxy term is constrained to around 90--100 pc/cm$^3$. The paper also claims a tight anticorrelation between $m_\\gamma$ and $f_{\\mathrm{IGM,0}}$.","pith_inferences":["Going beyond the paper, the unstated flat-geometry assumption may be the weakest link: the luminosity-distance identities used in Eqs. (4.5)--(4.6) rely on $d_L(z) = (1+z)c\\int_0^z dz'/H(z')$, which holds only for zero spatial curvature; a non-flat universe would introduce extra geometric factors and shift the inferred $m_\\gamma$.","A direct testable extension would measure the frequency-dependent arrival-time residual within individual FRBs (or among sub-bursts) and compare the implied $m_\\gamma$ with the ensemble value; a mismatch would expose the ensemble modeling of $\\mathrm{DM}_{\\mathrm{host}}$ and $f_{\\mathrm{IGM}}$.","The same integration-by-parts machinery could be applied to other frequency-dependent propagation effects, such as axion-photon mixing or Lorentz-invariance violation, replacing $m_\\gamma$ with the relevant coupling."],"forward_implications":["If the free-$f_{\\mathrm{IGM}}$ result is right, the photon rest mass is not zero at the roughly $2\\sigma$ level, with a central value around $3\\times10^{-50}$ kg, a factor of several above the previous model-independent limit from 32 FRBs ($m_\\gamma \\le 3.5\\times10^{-51}$ kg).","The strong anticorrelation between $m_\\gamma$ and $f_{\\mathrm{IGM,0}}$ means that any assumed baryon fraction directly biases the inferred photon mass; fixing $f_{\\mathrm{IGM,0}}$ can move the central value by roughly $11\\times10^{-51}$ kg.","The host-galaxy dispersion measure is constrained to $\\mathrm{DM}_{\\mathrm{host,0}} \\sim 90$--$100$ pc/cm$^3$, consistent with current FRB host-galaxy modeling.","The method is agnostic to the expansion history, so it can be applied as the sample of well-localized FRBs grows without needing to revisit the cosmological model."],"supporting_citations":[{"why":"Supplies the cosmological-model-independent integration-by-parts method that rewrites DM_IGM and DM_gamma in terms of the luminosity distance.","marker":"[27]"},{"why":"Provides the 1048 Pantheon type Ia supernovae from which the luminosity distances at FRB redshifts are reconstructed.","marker":"[93]"},{"why":"Gives the massive-photon dispersion relation and group velocity from which DM_gamma is derived.","marker":"[48]"},{"why":"Provides the average IGM electron-density expression that yields DM_IGM(z).","marker":"[49]"},{"why":"Supplies the assumed DM_halo = 50 pc/cm^3 and the host-galaxy/IGM fluctuation prescriptions used in the uncertainty budget.","marker":"[30]"},{"why":"Provides the trapezoidal numerical integration used to evaluate the d_L integrals in Eqs. (4.5) and (4.6).","marker":"[96]"},{"why":"Gives the earlier model-independent photon-mass upper limit from 32 FRBs that this work compares against.","marker":"[45]"},{"why":"Presents a LambdaCDM-based joint m_gamma-f_IGM analysis without the anticorrelation, used to contrast the model-independent result.","marker":"[44]"},{"why":"Fixes the Hubble constant H0 used in the DM_IGM prefactor.","marker":"[94]"},{"why":"Supplies the BBN value of Omega_b h^2 used in the IGM electron-density prefactor.","marker":"[99]"}],"fun_headline_variants":["Photon mass bound via FRB-SNe without cosmic model","FRBs + supernovae fix photon mass near 3e-50 kg","Model-independent photon mass from 68 FRBs and SNe","Cosmology-free limits on photon mass from FRB data"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation rests on the flat-universe relation between luminosity distance and expansion rate, $d_L(z) = (1+z)c\\int_0^z dz'/H(z')$, which the paper does not state explicitly; if spatial curvature is nonzero, the dispersion-measure terms acquire extra geometric factors and the inferred photon mass shifts.","fun_headline_variants_meta":{"raw":{"variants":["Photon mass bound via FRB-SNe without cosmic model","FRBs + supernovae fix photon mass near 3e-50 kg","Model-independent photon mass from 68 FRBs and SNe","Cosmology-free limits on photon mass from FRB data"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000626,"raw_usage":{"total_tokens":2975,"prompt_tokens":1099,"completion_tokens":1876,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":715,"completion_tokens_details":{"reasoning_tokens":1813}},"tokens_in":715,"tokens_out":1876,"duration_ms":15457,"temperature":1.0,"reasoning_tokens":1813,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:13:42.039889+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A single well-localized FRB with an independently measured luminosity distance (for example from a gravitational-wave standard siren or a Cepheid-calibrated host) would settle the flatness question: if the DM_IGM and DM_gamma residuals computed with the true d_L disagree with the fitted m_gamma at more than the quoted 1-sigma level, the massive-photon interpretation fails. More directly, a future sample of hundreds of FRBs should show residuals around the model that are statistically consistent with the fitted m_gamma; a significant excess that correlates with redshift and not with host properties would falsify the model.","supporting_citations":[],"review_version":1}