{"id":"62b62bc3-f26c-4721-9994-cccf9758589d","arxiv_id":"2411.17682","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A baryonification-based analytic model reproduces the host-galaxy dispersion measure distribution in CAMELS simulations and links it to baryonic feedback.","lead":"A new analytic model calculates how much of a fast radio burst's signal delay comes from electrons in its host galaxy. It matches hydrodynamic simulations and could make FRBs more reliable probes of cosmic structure.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Host DM PDF reproduction and mass/redshift scalings are contingent on the assumed FRB spatial profile; the paper's own Fig. 11 degeneracy analysis shows the gas parameters and predictions shift under alternative FRB priors, so the central claim is not yet robust.","rationale":"The reader identified p_FRB proportional to rho_star as the weakest assumption, and I agree it is the most load-bearing condition for the central claim. The paper's own analysis (Fig. 11) shows that the host-DM PDF only constrains the product of the FRB and gas profiles, not each independently. Consequently, the fitted BCM parameters and the predicted mass/redshift scaling are prior-dependent. This is not an internal inconsistency; the paper is transparent about the degeneracy and lists the simplistic stellar profile as a limitation (Sec. 5.7). However, it means the central claim as stated in the abstract, that the model reproduces the PDFs and their scaling, is established for one specific FRB model and should not be read as a general property of the BCM. The concrete test with alternative FRB distributions would determine whether the mass/redshift predictions are stable. The single-CV0 comparison and Poisson-only error model are additional weaknesses, but they are secondary to the FRB-prior dependence because the model's flexibility could absorb a different PDF shape. The verdict CONDITIONAL remains appropriate: the paper's conclusions are sound for the assumed prior but require confirmation under other FRB placement models before use in FRB cosmology.","tokens_in":21582,"tokens_out":18517,"duration_ms":160847,"concrete_test":"Refit the BCM to the SIMBA and TNG CV0 host-DM PDFs with three alternative FRB sampling distributions: p_FRB proportional to star-formation-rate density from the simulations, p_FRB uniform inside r_vir, and p_FRB proportional to rho_star with r_cut fixed to 0.1 r_vir. For each case, record the best-fit gas parameters and check (i) whether the simulated PDFs remain within the Poisson error band, (ii) whether the parameters stay inside the Giri and Schneider (2021) emulator prior, and (iii) whether the predicted mean DM-mass and DM-redshift relations still match Theis et al. (2024). If the gas parameters or the mass/redshift scaling change materially, the central claim must be explicitly restricted to the stellar-tracing FRB prior.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the BCM model reproduces the CAMELS host-DM PDF and its mass/redshift dependence. This is shown for p_FRB proportional to rho_star (Eq. 20). However, Section 5.6 and Figure 11 demonstrate a strong degeneracy: a more compact FRB profile is compensated by a shallower gas profile, so the same PDF is consistent with very different gas parameters depending on the FRB prior. Therefore the best-fit parameters in Table 1, the inferred feedback differences (SIMBA beta around 1.3 vs TNG around 0.3), and the predicted mean DM-mass/redshift curves (Figs. 7-8) are all conditional on this assumption. The paper acknowledges 'we cannot distinguish between these two cases without strong priors on the FRB distribution' (Sec. 5.6). Since real FRB populations may trace star-forming regions, globular clusters, or other stellar components, the claimed mass/redshift scaling and the feedback interpretation are not robust properties of the BCM; they are outputs of one chosen FRB prior. The validation against a single CV0 realization with a Poisson-only error model (Sec. 4) further weakens the strength of the 'reproduction' claim, but the FRB prior is the more fundamental dependence.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the baryonification model (BCM) to analytically model the host-galaxy contribution to fast radio burst dispersion measures. Assuming that FRBs trace the stellar density profile, the authors construct a PDF for the host DM by sampling sightlines and integrating the gas density profile. They fit the model's free parameters (with alpha fixed) to the host DM PDFs of the SIMBA and IllustrisTNG CV0 runs of CAMELS taken from Theis et al. (2024), and show the resulting fit in Figure 2. They then use the fitted parameters to predict the matter power spectrum suppression with the Giri & Schneider (2021) BCM emulator, finding agreement within the CAMELS CV scatter. The paper also presents the dependence of the host DM PDF on halo mass, redshift, and all BCM parameters, and highlights a degeneracy between the compactness of the FRB profile and the shallowness of the gas profile. The authors conclude that the model provides a flexible, physically motivated description of the host DM contribution and could enable joint modeling with baryonic feedback.","tokens_in":21896,"tokens_out":7367,"duration_ms":60939,"significance":"The proposed model is a useful conceptual step: it provides an analytic, computationally inexpensive framework for a quantity that is currently a major systematic in FRB cosmology. The parameter-dependence study (Figures 9-10) and the discussion of the FRB/gas degeneracy in Section 5.6 are informative and will be valuable for the community. The model's ability to jointly describe host DM PDFs and baryonic power-spectrum suppression within the same framework is attractive. However, the validation is currently limited to a single simulation realization and the inference is degenerate with the assumed FRB spatial prior, so the significance of the specific parameter values and the mass/redshift scalings is not yet established. Strengths of the manuscript: the model is explicitly specified, the fitting procedure is described, and the limitations are acknowledged in Section 5.7, which aids reproducibility.","major_comments":[{"comment":"The degeneracy between the FRB spatial distribution and the gas profile, demonstrated in Figure 11, is a load-bearing caveat. The best-fit BCM parameters in Table 1 and the resulting feedback interpretation (SIMBA beta≈1.3 vs TNG beta≈0.3) were obtained with p_FRB ∝ rho_star (Eq. 20). Figure 11 shows that a more compact FRB profile can be compensated by a shallower gas profile while still matching the same simulation PDFs, and the paper itself notes that 'we cannot distinguish between these two cases without strong priors on the FRB distribution' (Sec. 5.6). Since the true FRB distribution (e.g., star-forming regions or globular clusters) is unknown, the claimed reproduction of the mass- and redshift-dependence in Figures 7 and 8 is conditional on this prior and may not be robust.","section":"Section 5.6, Figure 11"},{"comment":"The validation is based on a single CV0 realization, and the error model in Equation (26) includes only Poisson scatter in the halo mass function, not cosmic variance across the full CV set. The grey bands in Figure 2 therefore underestimate the true simulation uncertainty. A fit to the ensemble of CV realizations, or at least a comparison against the PDF scatter among CV runs, is needed to support the claim that the BCM reproduces the host DM PDF. Without this, the excellent agreement in Figure 2 may be partly a consequence of the narrow error model.","section":"Section 4, Figure 2"},{"comment":"The power spectrum suppression check is not an independent prediction: the BCM emulator of Giri & Schneider (2021) is calibrated to hydrodynamic simulations with similar feedback prescriptions, and the comparison in Figure 3 uses the same CAMELS runs on which the PDF fit is based. The result that the BCM prediction lies within the grey scatter is therefore a consistency check rather than a validation of the model. The abstract and Section 4 should make this distinction explicit.","section":"Section 4, Figure 3"},{"comment":"The mass-dependence claim is stronger than the evidence. Figure 7 shows that the BCM does not reproduce the non-monotonic feature around 8×10^12 h^-1 M_sun seen in the SIMBA measurements, and the simulation points are shown without Poisson error bars. In addition, the comparison in Figure 8 is qualitative, with no cosmic variance estimate. The text should either quantify the level of agreement or soften the claim that the mass and redshift dependence is reproduced.","section":"Section 5.3, Figure 7"}],"minor_comments":[{"comment":"The text states that the model has 'eight free parameters', but Table 1 lists seven varied parameters plus alpha which is fixed to 2; please correct the count.","section":"Section 6"},{"comment":"The text says 'we simply measure the distribution function of halos, p_halo(M), in SIMBA' but Figure 2 also shows a TNG comparison; please clarify whether p_halo(M) is separately measured in IllustrisTNG or taken from Theis et al. (2024).","section":"Section 4"},{"comment":"The symbol 'B' appears to be a typographical artifact for a definition operator; it should be replaced with '≡' or '='.","section":"Equations (8) and (16)"},{"comment":"The phrase 'carried out up to a fiducial value of three (m=3)' should refer to ξ=3, not m.","section":"Section 5.1"},{"comment":"The caption uses 'viral mass' instead of 'virial mass'.","section":"Figure 4 caption"},{"comment":"The statement that 'strong feedback cannot be compensated by increasingly compact stellar distributions' is not quantitatively supported by Figure 11; please either provide a metric for the breakdown of the degeneracy or label it as a qualitative inference.","section":"Section 5.6"},{"comment":"The sentence 'if uses matches spherically symmetric profiles from simulations to the BCM' contains a typo; it should likely read 'if one matches'.","section":"Section 3.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript overlaps with Theis et al. (2024) through the shared use of the CAMELS CV0 data and overlapping authorship. This is not by itself a problem, but the validation dependence on that single realization should be addressed before publication. The paper's central claim is defensible, but the FRB-profile degeneracy and the single-realization validation are load-bearing and require additional work. The power-spectrum consistency check should also be described as such, not as an independent prediction."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The first application of baryonification to the FRB host DM PDF is genuinely new, and the paper does a solid job showing that a simple analytic model can reproduce the CAMELS PDFs and the matter power spectrum suppression within cosmic variance. The degeneracy result, that compact FRB profiles can be compensated by shallower gas profiles, is a real insight, and the authors are upfront about its consequence: without strong priors on the FRB distribution, the inferred gas parameters and feedback interpretation are not unique.\n\nThe validation rests on a single CV0 realization with an assumed Poisson error model, and the comparison PDFs come from a companion paper with shared authorship. That is not a fatal issue, but it does weaken the strength of the 'reproduction' claim. The bigger soft spot is the FRB profile prior: the best-fit parameters in Table 1, the SIMBA vs TNG feedback difference, and the mass/redshift curves all shift under alternative choices of p_FRB. The authors know this and say it in Sec. 5.6, so it is a stated limitation rather than a hidden flaw. The power-spectrum check is a consistency test, not fully independent validation, because it uses the same fitted parameters, though it is still a useful cross-check.\n\nThe transfer to survey-relevant halo masses is uncertain, and the authors are honest about that too. The two-halo term is shown to be negligible, which is good to see. The geometry and sampling are straightforward and the comparison to simulations is reproducible from the public BCM emulator.\n\nI would send this to peer review. The model is useful for the FRB community, and the limitations are identifiable and addressable. A referee should ask for a more systematic exploration of the FRB profile prior and a more realistic error model, but the central framework is sound.","headline":"A genuinely new and honest application of baryonification to FRB host DMs; the central degeneracy is clearly stated and the paper deserves a serious referee.","tokens_in":22399,"tokens_out":2995,"would_cite":true,"duration_ms":28110,"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":"An analytic baryonification model reproduces the probability distribution of host-galaxy dispersion measures for fast radio bursts, matching two hydrodynamic simulation suites and tying the host correction to baryonic feedback.","keywords":["fast radio bursts","dispersion measure","baryonification","baryonic feedback","host galaxies","circumgalactic medium","CAMELS simulations","halo model"],"falsifier":"Fit the same model to the host DM PDFs of the third CAMELS suite, ASTRID, or of a larger-box hydrodynamic simulation with different subgrid feedback: if no parameter set within the BCM prior can simultaneously reproduce its host DM PDFs and its power-spectrum suppression, the claimed self-consistency holds only for SIMBA and IllustrisTNG. On the data side, a sample of a few hundred localized FRBs with host halo mass estimates would test the predicted $\\langle \\mathrm{DM_{host}} \\rangle(M)$ scaling and the long high-mass tail of the PDF, while measured FRB offsets from their host centers would test the $p_{\\mathrm{FRB}} \\propto \\rho_{\\star}$ assumption directly.","tokens_in":21409,"feed_emoji":"📡","tokens_out":9800,"duration_ms":79864,"temperature":0.7,"pith_summary":"This paper claims that the host-galaxy contribution to the dispersion measure (DM) of fast radio bursts — currently one of the largest unknowns in FRB cosmology — can be described by a simple analytic model based on baryonification. The model assigns each halo physically motivated gas and stellar profiles, draws FRB sightlines from the stellar distribution, and integrates the free-electron column along each line of sight to build a full probability distribution of host DMs. That distribution reproduces the host DM PDFs measured in the SIMBA and IllustrisTNG runs of the CAMELS hydrodynamic simulations, including their halo-mass and redshift dependence, and the same fitted parameters predict the baryonic suppression of the matter power spectrum within cosmic variance. The paper's central physical finding is that the shape of the host DM distribution is governed by how compact the FRB population is relative to the gas: more concentrated FRB profiles must be paired with shallower gas profiles, a degeneracy that is real but bounded.","feed_headline":"One analytic model reproduces FRB host dispersion measures","feed_subtitle":"It matches two hydrodynamic simulation suites and ties the host correction to baryonic feedback.","key_machinery":"The load-bearing object is the baryonification (BCM) halo model, which remaps a dark-matter-only halo into dark-matter-baryon profiles with analytic gas, stellar, and collisionless components. The gas follows a cored double-power-law profile whose shape is set by a core radius $\\theta_{\\mathrm{co}}$, an ejection radius $\\theta_{\\mathrm{ej}}$, an outer slope $\\delta$, a transition slope $\\gamma$, and a mass-dependent inner slope $\\beta(M)$ controlled by the parameters $M_c$ and $\\mu$. The stellar profile — a truncated NFW profile with exponential cutoff governed by a cutoff radius $r_{\\mathrm{cut}}$ and slope $\\alpha$ — doubles as the FRB sampling distribution via $p_{\\mathrm{FRB}}(r) \\propto \\rho_{\\star}(r)$. The host DM probability distribution is built by sampling roughly $10^4$ sightlines from that distribution and integrating the electron column along each; the two-halo term is shown in an appendix to be negligible. The machinery converts a handful of feedback-related parameters into a concrete, testable prediction for the full host DM distribution, which is what makes the comparison to simulations and to future FRB samples possible.","core_discovery":"On its own terms, the paper establishes that the baryonification model (BCM) — a framework previously built to describe baryonic effects on the matter power spectrum — can be extended into a statistical model of FRB host DMs. The extension adapts the stellar profile to low-mass halos as a truncated NFW profile with an exponential cutoff and assumes FRBs are sampled proportionally to that stellar density, $p_{\\mathrm{FRB}}(r) \\propto \\rho_{\\star}(r)$. With eight free parameters fitted to the host DM PDFs of SIMBA and IllustrisTNG, the model reproduces the distributions, their growth in mean and width with halo mass, and their mild evolution from high to low redshift. As a consistency test, the same fitted parameters, passed through a BCM emulator, yield a matter power spectrum suppression consistent with the same simulations within cosmic variance. The paper further argues that the long-tailed, quasi-log-normal DM distributions seen in simulations arise naturally from the interplay of a compact stellar (FRB) profile and an extended gas profile, and that this interplay contains a bounded degeneracy: steeper FRB profiles demand shallower gas profiles, but not at arbitrary strength, so the two components are not fully interchangeable.","pith_inferences":["If future FRB localizations measure burst offsets from host centers for a large sample, the degenerate pair of profiles could be disentangled, turning the host DM PDF into a direct probe of the gas profile shape — a measurement the paper's own framework makes available but does not perform.","The model's strong sensitivity to the stellar cutoff radius suggests host DM statistics could constrain the FRB progenitor population itself, distinguishing star-formation-tracing bursts from those in old stellar environments.","The BCM parameters preferred by low-mass FRB hosts differ from those constrained by cosmic shear at cluster scales; a joint analysis would test whether feedback parameters must depend on halo mass, a question the paper leaves open.","Adding the partially neutral cold-gas component the paper flags as missing could turn host DM PDFs into a probe of the cold circumgalactic medium and its ionization state."],"forward_implications":["The same BCM parameters predict both the host DM PDF and the matter power spectrum suppression, so FRB cosmology analyses can marginalize over the host contribution in a way that is consistent with baryonic feedback constraints from other probes.","The host DM PDF is sensitive enough to feedback parameters that it can serve as a complementary probe of baryonic physics on galactic scales, where cosmic shear has little leverage.","Because the mean host DM grows with halo mass and the PDF broadens toward log-normal shape, the host correction cannot be approximated as a constant or a fixed scatter in FRB samples spanning different host populations.","The bounded nature of the FRB-gas degeneracy means that combining DM PDFs with even weak priors on the FRB distribution can break the degeneracy, allowing a joint measurement of the gas profile and the stellar profile."],"supporting_citations":[{"why":"Original baryonification method whose halo remapping this paper extends to FRB host DMs.","marker":"Schneider & Teyssier (2015)"},{"why":"Extended BCM with analytic gas, stellar, and collisionless profiles plus adiabatic relaxation.","marker":"Schneider et al. (2019)"},{"why":"BCM emulator used to predict the matter power spectrum suppression from the fitted host-DM parameters.","marker":"Giri & Schneider (2021)"},{"why":"Source of the CAMELS host DM PDF measurements that the model is fitted to and compared against.","marker":"Theis et al. (2024)"},{"why":"The CAMELS simulation suite and its cosmology, providing the hydrodynamic runs used for validation.","marker":"Villaescusa-Navarro et al. (2021)"},{"why":"Modified NFW gas profile used as a benchmark for the electron distribution in halos.","marker":"Prochaska & Zheng (2019)"},{"why":"Concentration-mass relation used to construct the halo density profiles in the model.","marker":"Dutton & Macciò (2014)"},{"why":"Abundance-matching based double power law that sets the stellar fraction in low-mass halos.","marker":"Moster et al. (2012)"}],"fun_headline_variants":["Analytic model reproduces FRB host DM distributions from simulations","Baryonification model links FRB host DMs to gas feedback","FRB host DM shapes trace baryonic feedback in halos","One model explains FRB host DMs and baryonic effects"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes fast radio bursts live where the stars are — FRB sightlines are drawn from the stellar density profile — and every inferred gas property, as well as the central compactness degeneracy, shifts if the true FRB population instead follows star-forming gas, globular clusters, or some other spatial distribution.","fun_headline_variants_meta":{"raw":{"variants":["Analytic model reproduces FRB host DM distributions from simulations","Baryonification model links FRB host DMs to gas feedback","FRB host DM shapes trace baryonic feedback in halos","One model explains FRB host DMs and baryonic effects"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000665,"raw_usage":{"total_tokens":3125,"prompt_tokens":1124,"completion_tokens":2001,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":740,"completion_tokens_details":{"reasoning_tokens":1928}},"tokens_in":740,"tokens_out":2001,"duration_ms":13257,"temperature":1.0,"reasoning_tokens":1928,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T11:50:38.720350+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same model to the host DM PDFs of the third CAMELS suite, ASTRID, or of a larger-box hydrodynamic simulation with different subgrid feedback: if no parameter set within the BCM prior can simultaneously reproduce its host DM PDFs and its power-spectrum suppression, the claimed self-consistency holds only for SIMBA and IllustrisTNG. On the data side, a sample of a few hundred localized FRBs with host halo mass estimates would test the predicted $\\langle \\mathrm{DM_{host}} \\rangle(M)$ scaling and the long high-mass tail of the PDF, while measured FRB offsets from their host centers would test the $p_{\\mathrm{FRB}} \\propto \\rho_{\\star}$ assumption directly.","supporting_citations":[],"review_version":1}