arxiv: 2601.18784 · v2 · submitted 2026-01-26 · 🌌 astro-ph.CO · astro-ph.HE

Recognition: no theorem link

Baryonification III: An accurate analytical model for the dispersion measure probability density function of fast radio bursts

MohammadReza Torkamani , Robert Reischke , Michael Kova\v{c} , Andrina Nicola , Jozef Bucko , Alexandre Refregier , Sambit K. Giri , Aurel Schneider

show 1 more author

Steffen Hagstotz

Authors on Pith no claims yet

Pith reviewed 2026-05-16 10:25 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.HE

keywords fast radio burstsdispersion measureprobability distribution functionbaryonificationhalo mass functionbaryonic feedbackcosmological simulations

0 comments

The pith

An analytical model using baryonification predicts the dispersion measure PDF of fast radio bursts and matches hydrodynamical simulations from z=0 to 5.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an analytical framework to compute the one-point PDF of dispersion measures for fast radio bursts by convolving contributions from halos of varying masses and redshifts, using the halo mass function and bias inside the baryonification model. This supplies a fast alternative to full hydrodynamical runs while reproducing the DM statistics seen in IllustrisTNG across the full redshift range examined. A sympathetic reader would care because the approach directly ties adjustable gas-profile parameters to observable DM distributions, allowing future FRB data to constrain baryonic feedback without repeated expensive simulations. The work also shows that the PDF shape is controlled mainly by the parameters Mc, mu and delta and that a log-normal form approximates the distribution well for samples of a few hundred bursts.

Core claim

We develop an analytical framework to predict the one-point probability distribution function of dispersion measures for fast radio bursts within the baryonification model. By applying the halo mass function and halo bias we convolve contributions from individual halos across masses and redshifts to obtain the large-scale structure part of the DM PDF. Validation against consistency-check simulations and direct comparison with IllustrisTNG shows excellent agreement from z=0 to z=5. The model remains self-consistent when gas profiles are fitted or when the PDF is predicted, with the parameters Mc, mu and delta identified as the primary drivers of PDF shape.

What carries the argument

Convolution of halo contributions via the halo mass function and bias, which folds the baryonification gas density profiles into the cosmological DM PDF.

If this is right

The gas-profile parameters Mc, mu and delta primarily control the shape of the predicted DM PDF.
The analytical model yields consistent results whether it is used to fit gas profiles or to predict the PDF.
A log-normal distribution provides an adequate description of the DM PDF for samples containing only a few hundred FRBs.
FRB observations can constrain baryonic feedback processes through the self-consistent link between gas profiles and integrated DM statistics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Large future FRB catalogs could be used to measure the baryonification parameters directly from the observed DM distribution.
The demonstrated agreement with IllustrisTNG suggests the model can forecast DM statistics under varied cosmological parameters or feedback strengths.
The same convolution approach could be applied to predict probability distributions for other line-of-sight integrals such as the thermal Sunyaev-Zeldovich effect.

Load-bearing premise

The baryonification model's chosen parametrization of gas profiles inside halos correctly captures the true baryon distribution across the relevant range of halo masses and redshifts.

What would settle it

A direct measurement of the DM PDF from a large sample of FRBs with known redshifts showing a clear mismatch in the high-DM tail or overall shape compared with the analytical prediction would falsify the accuracy of the halo convolution or the underlying gas profiles.

read the original abstract

We develop an analytical framework to predict the one-point probability distribution function (PDF) of dispersion measures (DMs) for fast radio bursts (FRBs) within the baryonification (BFC) model. BFC provides a computationally efficient alternative to expensive hydrodynamical simulations for modelling baryonic effects on cosmological scales. By applying the halo mass function and halo bias, we convolve contributions from individual halos across a range of masses and redshifts to derive the large-scale structure contribution to the DM PDF. We validate our analytical predictions against consistency-check simulations and compare them with the IllustrisTNG hydrodynamical simulation over the redshift range $ z = 0$ to $z = 5$, demonstrating excellent agreement. We demonstrate that our model produces consistent results when fitting gas profiles and predicting the PDF, and vice versa. We show that the BFC parameters controlling the gas profile, particularly the halo mass scale ($M_\mathrm{c}$), mass-dependent slope ($\mu$), and outer truncation ($\delta$), are the primary drivers of the PDF shape. Additionally, we investigate the validity of the log-normal approximation commonly used for DM distributions, finding that it provides a sufficient description for a few hundred FRBs. Our work provides a self-consistent model that links gas density profiles to integrated DM statistics, enabling future constraints on baryonic feedback processes from FRB observations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper builds a workable analytical convolution for the DM PDF of FRBs inside the baryonification model and reports solid agreement with IllustrisTNG from z=0 to 5.

read the letter

The main point to take away is that the authors have constructed an analytical model for the probability density function of dispersion measures from fast radio bursts by convolving individual halo contributions using the halo mass function and bias, all within their baryonification framework, and this matches the IllustrisTNG simulation results from redshift zero out to five. This extends their previous baryonification work, which was more about two-point functions or average profiles, to the one-point statistics of the DM. The approach is efficient because it avoids running new hydro simulations for each parameter set. They demonstrate that the PDF can be predicted from the gas profile parameters and that fitting the profiles to simulations gives back a consistent PDF prediction. The parameters Mc, mu, and delta control the shape, with Mc being particularly important. They also check that the log-normal distribution often assumed for DMs is a decent approximation when you have a few hundred bursts. The strengths are the direct comparison to a major hydro run and the internal consistency test between profile fitting and PDF calculation. This makes the model more credible for use in cosmology. The convolution itself seems to follow standard large-scale structure techniques without obvious gaps. The softer parts are that the accuracy still hinges on how faithfully the three-parameter gas profile represents real baryonic physics across all halo masses and redshifts. If there are residuals in the profile fit at high or low masses, they could affect the tails of the DM PDF, which matter for rare events. The paper reports excellent agreement, but details on the quantitative level of that agreement, such as chi-squared values or Kolmogorov-Smirnov statistics, would help judge how close it really is. Since the parameters are tuned to simulations, there's some built-in dependence, though the bidirectional check mitigates it. This work is relevant for researchers studying baryonic feedback with FRB observations or those building mock catalogs for next-generation surveys. It provides a faster alternative to full simulations for exploring parameter space. I would send this to peer review. The central claim holds up based on the validation described, and the method is transparent enough for referees to assess.

Referee Report

2 major / 3 minor

Summary. The manuscript develops an analytical framework to predict the one-point PDF of FRB dispersion measures within the baryonification (BFC) model. It convolves individual halo contributions to the DM using the halo mass function and halo bias, with gas profiles controlled by the three parameters M_c, μ, and δ. The model is validated against consistency-check simulations and compared to IllustrisTNG hydrodynamical results over z=0 to z=5, reporting excellent agreement. It further demonstrates internal consistency when gas-profile parameters are fitted to simulations and then used to predict the PDF (and vice versa), identifies these parameters as the primary drivers of PDF shape, and tests the adequacy of the log-normal approximation for samples of a few hundred FRBs.

Significance. If the central results hold, the work supplies a computationally efficient, self-consistent link between baryonic gas profiles and FRB DM statistics that could enable observational constraints on feedback processes without requiring full hydrodynamical runs. The reported agreement with IllustrisTNG across a wide redshift range and the bidirectional consistency checks are concrete strengths that would make the framework useful for cosmological analyses of upcoming FRB catalogs.

major comments (2)

[Results section on IllustrisTNG comparison] Validation against IllustrisTNG: the claim of 'excellent agreement' over z=0–5 requires quantitative metrics (e.g., Kolmogorov-Smirnov statistics, binned residuals, or χ² values) rather than qualitative statements; without these, it is difficult to assess how well the analytical PDF reproduces the simulation tails and peak location.
[Consistency-checks subsection] Parameter consistency: the bidirectional fitting (gas profiles ↔ PDF) uses the same three BFC parameters (M_c, μ, δ) for both steps; an explicit test that fixes parameters from one simulation suite and predicts the PDF in an independent suite would reduce the risk of circularity noted in the abstract.

minor comments (3)

[Analytical framework section] The convolution integral that combines the halo mass function, bias, and gas profile should be written explicitly (including redshift integration limits) to allow readers to reproduce the PDF derivation without ambiguity.
[Figures showing PDF comparisons] Figure captions for the PDF comparisons should state the number of FRBs per realization and the binning scheme used for the histograms so that the visual agreement can be evaluated quantitatively.
[Log-normal approximation subsection] The statement that the log-normal approximation is 'sufficient' for a few hundred FRBs would benefit from a quantitative threshold (e.g., the sample size at which the Kolmogorov-Smirnov distance exceeds a chosen significance level).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and positive recommendation for minor revision. We address the major comments point by point below.

read point-by-point responses

Referee: [Results section on IllustrisTNG comparison] Validation against IllustrisTNG: the claim of 'excellent agreement' over z=0–5 requires quantitative metrics (e.g., Kolmogorov-Smirnov statistics, binned residuals, or χ² values) rather than qualitative statements; without these, it is difficult to assess how well the analytical PDF reproduces the simulation tails and peak location.

Authors: We agree that quantitative metrics would provide a more rigorous validation. In the revised manuscript we will add Kolmogorov-Smirnov statistics, binned residuals, and χ² values comparing the analytical PDF to the IllustrisTNG results over z=0–5, with particular attention to the peak location and high-DM tails. revision: yes
Referee: [Consistency-checks subsection] Parameter consistency: the bidirectional fitting (gas profiles ↔ PDF) uses the same three BFC parameters (M_c, μ, δ) for both steps; an explicit test that fixes parameters from one simulation suite and predicts the PDF in an independent suite would reduce the risk of circularity noted in the abstract.

Authors: The bidirectional test is performed on independent consistency-check simulations that are distinct from the IllustrisTNG runs. To address the concern, we will add explicit text clarifying the independence of the simulation suites and the fitting procedure. We will also include a cross-validation exercise using parameters fitted on one consistency-check realization to predict the PDF on a second, held-out realization from the same suite. revision: partial

Circularity Check

1 steps flagged

Moderate circularity from bidirectional fitting of BFC gas parameters to simulations before PDF prediction

specific steps

fitted input called prediction [Abstract]
"We demonstrate that our model produces consistent results when fitting gas profiles and predicting the PDF, and vice versa."

The BFC parameters (Mc, mu, delta) are adjusted to match simulated gas profiles; the analytical PDF is then obtained by convolving those same parameters with the halo mass function and bias. Bidirectional consistency therefore reduces the 'prediction' to a re-expression of the fit rather than an independent test.

full rationale

The derivation convolves HMF and bias using the three-parameter BFC gas profile (Mc, mu, delta) to obtain the DM PDF, then validates by showing the same parameters produce consistent results when fitted directly to gas profiles versus when used for PDF prediction. This creates a fitted-input-called-prediction loop where agreement with IllustrisTNG is partly enforced by the shared parameters rather than emerging from an independent first-principles calculation. No self-citation chain or self-definitional reduction is load-bearing for the core convolution step itself, and external simulation checks provide some independent grounding, keeping the overall circularity moderate rather than forcing the result by construction.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard cosmological halo statistics and the BFC parametrization of gas profiles. Free parameters are the primary adjustable elements; no new physical entities are introduced.

free parameters (3)

M_c
Halo mass scale controlling the gas profile, identified as a primary driver of PDF shape.
mu
Mass-dependent slope of the gas profile, primary driver of PDF shape.
delta
Outer truncation parameter for gas profiles, primary driver of PDF shape.

axioms (2)

domain assumption Halo mass function and halo bias accurately describe the distribution and clustering of halos across masses and redshifts
Invoked to convolve individual halo contributions into the large-scale DM PDF.
domain assumption Baryonification model provides a computationally efficient and accurate approximation to baryonic effects on cosmological scales
Basis for deriving gas profiles that are then integrated into the DM PDF.

pith-pipeline@v0.9.0 · 5592 in / 1637 out tokens · 51817 ms · 2026-05-16T10:25:19.434823+00:00 · methodology