---
title: Component-wise Baryonification (BFC)
url: https://www.emergentmind.com/topics/component-wise-baryonification-bfc
type: topic
---

# Component-wise Baryonification (BFC)

to=arxiv_search 《凤凰大参考json
{"query":"Component-wise Baryonification BFC baryonification 2507.07991 2507.07892 2601.18784","max_results":10,"sort_by":"submittedDate","sort_order":"descending"}【อ่านข้อความเต็มjson
[{"arxiv_id":"2601.18784","title":"Baryonification III: An accurate analytical model for the dispersion measure probability density function of fast radio bursts","authors":["Matthieu Schaller","João M. Zorrilla Matilla","Lucas Lombriser"],"abstract":"We develop a fully analytical framework for predicting the one-point probability distribution function (PDF) of dispersion measures (DM) for fast radio bursts (FRBs) using 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 across a range of redshifts up 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 (Mc), mass-dependent slope (μ), and outer truncation (δ), 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."},{"arxiv_id":"2507.07991","title":"Baryonification II: Constraining feedback with X-ray and kinematic Sunyaev-Zel'dovich observations","authors":["Alessandro R. Murk","Martijn M. S. L. Brouwer","João M. Zorrilla Matilla","Lukas Schaller","Lucas Lombriser"],"abstract":"Baryonic feedback alters the matter distribution on small and intermediate scales, posing a challenge for precision cosmology. The new, component-wise baryonification (BFC) approach provides a self-consistent framework to model feedback effects for different observables. In this paper we use this framework to fit kinematic Sunyaev-Zel'dovich (kSZ) observations from the Atacama Cosmology Telescope (ACT) alongside halo X-ray gas fractions from eROSITA, investigating baryonic feedback in a cosmological context. We first show that the kSZ data from ACT is consistent with the gas fractions from eROSITA, both suggesting a feedback model that is stronger than what is assumed in most hydrodynamical simulations. This finding is in contrast to older, pre-eROSITA gas fraction measurements that point towards weaker feedback in tension with the kSZ results. We suspect these discrepancies to be due to selection bias in the pre-eROSITA sample, or differences in halo mass estimation between the two data sets. In a further step, we use the BFC model to predict the baryonic suppression of the matter power spectrum. Based on our combined fit to data from ACT and eROSITA, we find a power spectrum suppression that exceeds the percent-level at modes above k=0.3-0.6 h/Mpc, growing to 2-8 percent at k=1 h/Mpc, and to 20-25 percent at k=5 h/Mpc, consistent with strong-feedback hydrodynamical simulations. Finally, we compare our best-fitting model to the observed gas density and pressure profiles of massive galaxy clusters from the X-COP sample, finding excellent agreement. These results show that BFC provides a self-consistent picture of feedback across mass- and length scales as well as different cosmological observables, thus making it promising for applications to multiwavelength studies to jointly constrain cosmology and baryonic effects."},{"arxiv_id":"2507.07892","title":"Baryonification: An alternative to hydrodynamical simulations for cosmological studies","authors":["João M. Zorrilla Matilla","Lukas Schaller","Anders C. Bygott","Martijn M. S. L. Brouwer","Alessandro R. Murk","Viraj Pandya","Cora Uhlemann","Lucas Lombriser"],"abstract":"We present an improved baryonification (BFC) model that modifies dark-matter-only N-body simulations to generate particle-level outputs for gas, dark matter, and stars. Unlike previous implementations, our approach first splits each simulation particle into separate dark matter and baryonic components, which are then displaced individually using the BFC technique. By applying the hydrostatic and ideal gas equations, we assign pressure and temperature values to individual gas particles. The model is validated against hydrodynamical simulations from the FLAMINGO and TNG suites (which feature varied feedback prescriptions) showing good agreement at the level of density and pressure profiles across a wide range of halo masses. As a further step, we calibrate the BFC model parameters to gas and stellar mass ratio profiles from the hydrodynamical simulations. Based on these calibrations, we baryonify N-body simulations and compare the resulting total matter power spectrum suppressions to the ones from the same hydrodynamical simulation. Carrying out this test of the BFC method at each redshift individually, we obtain a 2 percent agreement up to k=5 h/Mpc across all tested feedback scenarios. We also define a reduced, 2+1 parameter BFC model that simultaneously accounts for feedback variations (2 parameters) and redshift evolution (1 parameter). The 2+1 parameter model agrees with the hydrodynamical simulations to better than 2.5 percent over the scales and redshifts relevant for cosmological surveys. Finally, we present a map-level comparison between a baryonified N-body simulation and a full hydrodynamical run from the TNG simulation suite. Visual inspection of dark matter, gas, and stellar density fields, along with the integrated pressure map, shows promising agreement. Further work is needed to quantify the accuracy at the level of observables."}]
Component-wise Baryonification (BFC) is a baryonification framework that modifies gravity-only or dark-matter-only \(N\)-body outputs to emulate the effects of gas, stars, and baryonic feedback by treating dark matter, gas, and stars as independent fields at the particle level [2507.07892]. In its modern form, BFC displaces particles so that the final mass distribution is a physically motivated superposition of dark matter, hot and cold gas, and stellar components, while also providing a self-consistent description of hot-gas thermodynamics that can be mapped to kinematic Sunyaev-Zel'dovich (kSZ), X-ray, and related observables [2507.07991]. The framework has been developed as an alternative to hydrodynamical simulations for cosmological studies, validated against FLAMINGO and TNG, used to jointly constrain feedback with ACT and eROSITA data, and extended to a fully analytical model for the dispersion-measure probability density function of fast radio bursts (FRBs) [2507.07892] [2507.07991] [2601.18784].

## 1. Concept and defining characteristics

BFC is a post-processing tool that modifies gravity-only \(N\)-body outputs to emulate the effects of gas, stars, and baryonic feedback, but with a crucial change relative to earlier baryonification schemes: it treats dark matter, gas, and stars as independent fields at the particle level [2507.07892]. Each original DMO simulation particle is duplicated into a dark-matter particle and a baryonic particle; the two start at the same position but undergo different radial displacements around halo centers, with particle masses renormalized to match the cosmic fractions \(f_{\rm dm}=\Omega_{\rm dm}/\Omega_m\) and \(f_{\rm bar}=\Omega_b/\Omega_m\) [2507.07892].

The framework differs from earlier baryonification methods in three key ways. First, it performs component-wise modeling across all matter species: gas and stars are decomposed into physically distinct subcomponents, each with its own mass fraction \(f_i(M,z)\) and profile \(\rho_i(r|M,z)\), enabling simultaneous predictions for gas- and star-sensitive observables alongside total-matter effects such as matter power suppression [2507.07991]. Second, it implements self-consistent hot-gas thermodynamics by solving for the total pressure via hydrostatic equilibrium and then partitioning it into thermal and non-thermal components with an explicit model for non-thermal pressure support [2507.07991]. Third, it provides a unified cosmology–feedback forward model in which the same halo component fractions and profiles generate kSZ and X-ray gas-fraction predictions and the dark-matter response that suppresses the matter power spectrum, enabling joint constraints from multiwavelength data without re-calibration to hydrodynamical simulations [2507.07991].

A recurring misconception is to treat BFC as merely a fitted correction to the matter power spectrum. The published formulation is broader: it yields particle-level outputs for three components—DM, gas, and stars—and, for gas particles, additionally assigns a thermal pressure and temperature using hydrostatic equilibrium plus an empirical non-thermal pressure prescription and the ideal-gas law [2507.07892]. This suggests that BFC is intended as a physically organized field-level construction rather than only a summary transfer function for clustering.

## 2. Mass decomposition, profiles, and transport map

BFC starts from a DMO halo described by a truncated NFW profile and a two-halo term [2507.07991]. For a halo with \(M_{200{\rm c}}\) and concentration \(c_{200}\), \(r_{200}\equiv r_{200{\rm c}}\) is the radius enclosing 200 times the critical density [2507.07991]. The matter content is split into five components: dark matter (\({\rm dm}\)), hot bound gas (\({\rm hga}\)), cold/inner gas (\({\rm iga}\)), central galaxy stars (\({\rm cga}\)), and satellite galaxy stars (\({\rm sga}\)) [2507.07991].

The initial and final cumulative radial mass profiles define the radial displacements that map the initial to final density profiles while enforcing mass conservation,
\[
M_{i,f}(r)=\int_0^r ds\, s^2 \rho_{i,f}(s).
\]
Particles are displaced so that the final density equals the prescribed component sum [2507.07991]. In the formulation used for the transport map, the initial and final “one-halo + two-halo” forms are
\[
\rho_i(r)=f_{\rm dm}\,[\rho_{\rm nfw}(r)+\rho_{2{\rm h}}(r)],\qquad
\rho_f(r)=\rho_{\rm dm}(r)+f_{\rm dm}\rho_{2{\rm h}}(r),
\]
for dark matter, and
\[
\rho_i(r)=f_{\rm bar}\,[\rho_{\rm nfw}(r)+\rho_{2{\rm h}}(r)],\qquad
\rho_f(r)=\rho_{\rm gas}(r)+\rho_{\rm stars}(r)+f_{\rm bar}\rho_{2{\rm h}}(r),
\]
for baryons, where \(\rho_{\rm stars}=\rho_{\rm cga}+\rho_{\rm sga}\) and \(\rho_{\rm gas}=\rho_{\rm hga}+\rho_{\rm iga}\) [2507.07991].

The initial halo profile is a truncated NFW,
\[
\rho_{\rm nfw}(r)=\rho_{{\rm nfw},0}\,[x_s(1+x_s)^2]^{-1}[1+x_t^2]^{-2},
\]
with \(x_s\equiv r/r_s\), \(r_s=r_{200}/c_{200}\), \(x_t\equiv r/r_t\), \(r_t=\epsilon r_{200}\), and truncation parameter \(\epsilon=\epsilon_0+\epsilon_1\nu\) with \(\epsilon_0=4.0\) and \(\epsilon_1=0.5\) [2507.07991]. The normalization \(\rho_{{\rm nfw},0}\) ensures \(\int \rho_{\rm nfw}\,d^3r=M_{200{\rm c}}\) [2507.07991].

The hot bound gas profile is modeled with a smooth core, a mass-dependent interior slope, and an outer steepening or truncation,
\[
\rho_{\rm hga}(r)=\rho_{{\rm hga},0}\,[1+x_c^\alpha]^{-\beta(M)/\alpha}[1+x_t^\gamma]^{-\delta/\gamma},
\]
with \(x_c\equiv r/r_c\), \(r_c=\theta_c r_{200}\), \(x_t\equiv r/r_t\), \(r_t=\epsilon r_{200}\), and
\[
\beta(M)=\frac{3(M/M_c)^\mu}{1+(M/M_c)^\mu}.
\]
The mass dependence \(\beta(M)\) captures the relative flattening of gas profiles at lower halo masses where feedback is more efficient [2507.07991]. In the implementation paper, \(\alpha=1\) and \(\gamma=3/2\) are fixed in the hot-gas profile [2507.07892].

The cold or inner gas fraction is tied to the central galaxy fraction, \(f_{\rm iga}=c_{\rm iga}f_{\rm cga}\), and is centrally concentrated; in practice, it is a small fraction of the baryons at group and cluster scales and does not affect kSZ significantly [2507.07991]. The central galaxy stellar profile is
\[
\rho_{\rm cga}(r)=\frac{f_{\rm cga}M_{\rm tot}}{4\pi R_h}\,r^{-2}\exp(-r/R_h),
\qquad R_h=0.03\,r_{200},
\]
while satellite stars are assumed to trace a rescaled NFW cumulative mass,
\[
\rho_{\rm sga}(r)=\frac{f_{\rm sga}}{4\pi r^2}\frac{d}{dr}M_{\rm nfw}(\xi r),
\]
where \(\xi\) encodes net contraction or expansion effects [2507.07991].

At the particle-assignment stage, the displaced baryonic particle is stochastically assigned to become a gas or a star particle with radius-dependent probabilities
\[
P_i(r)=\frac{\rho_i(r)}{\rho_{\rm gas}(r)+\rho_{\rm star}(r)},\qquad i\in\{{\rm gas,star}\},
\]
stars are assumed to reside within the virial radius, and baryons not belonging to any halo after displacement are assigned to gas [2507.07892]. A practical rule preserves bound substructure: if a baryonic particle belongs to a neighboring halo, it is displaced using the dark-matter prescription so baryons and DM in that satellite are moved coherently [2507.07892].

## 3. Thermodynamics and dark-matter response

A distinctive element of BFC is that hot-gas thermodynamics is not imposed by a fixed, simulation-calibrated pressure profile. Instead, the total pressure profile \(P_{\rm tot}\) is obtained from hydrostatic equilibrium,
\[
\frac{dP_{\rm tot}}{dr}=-\rho_{\rm gas}(r)\frac{G\,M(<r)}{r^2},
\]
and is then split into thermal and non-thermal contributions through a radial and redshift-dependent non-thermal fraction based on the Shaw et al. model [2507.07991]. In the observational analysis,
\[
\frac{P_{\rm nt}}{P_{\rm tot}}(r,z)=\alpha(z)\left[\frac{r}{R_{500}}\right]^{n_{\rm nt}},
\qquad
\alpha(z)=\alpha_{0,{\rm nt}}\,f(z),
\]
with
\[
f(z)=\min\left[(1+z)^\beta,\,(f_{\max}-1)\tanh(\beta z)+1\right],
\]
\(\beta=0.5\), and \(\alpha_{0,{\rm nt}}=0.1\) fixed in that paper; the cap \(f_{\max}=4^{-n_{\rm nt}/\alpha_0}\) prevents \(P_{\rm nt}>P_{\rm tot}\) [2507.07991]. The thermal pressure is then \(P_{\rm th}=P_{\rm tot}-P_{\rm nt}\) [2507.07991]. In the implementation paper, the non-thermal correction is given as
\[
P_{\rm th}(r)=\left[1-\alpha_{\rm nth}(z)(r/r_{200})^{n_{\rm nth}}\right]P_{\rm tot},
\]
with \(\alpha_{{\rm nth},0}=0.18\), \(\beta_{\rm nth}=0.5\), and \(n_{\rm nth}=0.8\), and \(P_{\rm th}=0\) where the bracket becomes negative in the outer halo [2507.07892].

Given the gas mass density, the electron number density used for kSZ and X-ray applications is
\[
n_e(r)=\frac{X_H+1}{2m_{\rm amu}}\,\rho_{\rm gas}(r),
\]
with hydrogen mass fraction \(X_H=0.76\) [2507.07991]. The electron pressure is
\[
P_e(r)=n_e(r)k_B T_e(r)=\frac{n_e(r)}{n_g(r)}P_{\rm th}(r)=\frac{\mu}{\mu_e}\frac{\rho_{\rm gas}(r)}{m_p}k_B T_e(r),
\]
where \(m_p\) is the proton mass and \(\mu,\mu_e\) are mean molecular weights for total gas and electrons, respectively [2507.07991]. In the implementation paper, the ideal-gas temperature is
\[
T(r)=\frac{m_P\mu_m P_{\rm th}(r)}{k_B\rho_{\rm hga}(r)},
\]
with \(\mu_m=0.6125\) [2507.07892].

BFC also includes an explicit dark-matter back-reaction model. The empirical mapping is written in terms of \(\xi=r_i/r_f\), with terms in the relation producing contraction in the inner halo from stars and cold gas and expansion from ejected hot gas [2507.07892]. The calibrated values \(Q_0=0.075\) and \(Q_1=0.25\) are fixed across suites, while \(Q_2\) depends on simulation resolution, taking the values \(0.5\) for FLAMINGO m9, \(0.7\) for FLAMINGO m8, and \(0.8\) for TNG-300 [2507.07892]. This component is central to the framework’s claim that the same baryonic redistribution governs both gas-sensitive observables and the suppression of the nonlinear matter power spectrum.

## 4. Parameterization, calibration, and reduced models

The full BFC parameter set comprises gas parameters \(\theta_c\), \(M_c\), \(\mu\), and \(\delta\), and stellar or cold-gas parameters \(c_{\rm IGA}\), \(\eta\), \(d\eta\), and \(N_{\rm star}\); back-reaction parameters \(Q_0\), \(Q_1\), and \(Q_2\) are fixed as above [2507.07892]. The stellar fractions follow a double power-law inspired by abundance matching,
\[
f_{\rm star}(M)=N_{\rm star}\left[\left(\frac{M}{M_{\rm star}}\right)^{-\zeta_{\rm star}}+\left(\frac{M}{M_{\rm star}}\right)^{-\eta}\right]^{-1},
\]
\[
f_{\rm cga}(M)=N_{\rm star}\left[\left(\frac{M}{M_{\rm star}}\right)^{-\zeta_{\rm star}}+\left(\frac{M}{M_{\rm star}}\right)^{-(\eta+d\eta)}\right]^{-1},
\qquad
f_{\rm sga}=f_{\rm star}-f_{\rm cga},
\]
with the cold-gas fraction given by \(f_{\rm iga}=c_{\rm iga}f_{\rm cga}\), and the hot-gas fraction from cosmic closure, \(f_{\rm gas}=f_{\rm bar}-f_{\rm star}-f_{\rm iga}\) [2507.07991]. In the implementation paper, \(M_{\rm star}=2.5\times10^{11}\,M_\odot/h\) and \(\zeta=1.376\) are specified in the stellar-fraction model [2507.07892].

Calibration targets in the implementation paper are gas and stellar mass-ratio profiles,
\[
R_{\rm gas}(r)=\frac{M_{\rm gas}(r)}{M_{\rm tot}(r)},\qquad
R_{\rm star}(r)=\frac{M_{\rm star}(r)}{M_{\rm tot}(r)},
\]
measured from the FLAMINGO runs m8, m9, m9 Jet, m9 \(f_{\rm gas}\)-8\(\sigma\), and from TNG-300 [2507.07892]. Carrying out the test of the BFC method at each redshift individually, the authors obtain a 2 percent agreement up to \(k=5\,h/{\rm Mpc}\) across all tested feedback scenarios; at \(z=0\), the full 8-parameter and reduced 2-parameter BFC agree with hydro to \(\lesssim 2\%\) for \(0.1<k<10\,h/{\rm Mpc}\) [2507.07892]. The reduced “2+1 parameter” BFC model takes \(M_c\) and \(\delta\) as free feedback parameters and introduces one redshift-evolution parameter \(\alpha_c\) through
\[
\log\left[\frac{M_c}{M_{c,0}}\right]=-\alpha_c z,
\]
with \(\delta\) approximated constant with redshift [2507.07892]. The 2+1 parameter model agrees with the hydrodynamical simulations to better than 2.5 percent over the scales and redshifts relevant for cosmological surveys [2507.07892].

For survey applications, the implementation paper recommends using the reduced 2-parameter BFC at each redshift with \(\theta_c(z)=0.3(1+z)^{1/2}\) and \(c_{\rm IGA}(z)=0.1(1+z)^{3/2}\), while fixing stellar parameters via abundance matching or prior fits to survey stellar mass–halo mass relations [2507.07892]. For a redshift-spanning fit, the same paper recommends the 2+1 parameter model with \(M_c(z)\) following the redshift law above and \(\delta\) constant [2507.07892].

## 5. Mapping to observables and empirical constraints

BFC was designed to map component profiles into observables without changing the underlying physical parameterization. For kSZ, the temperature shift is
\[
\frac{\Delta T_{\rm kSZ}}{T_{\rm CMB}}=-\frac{\sigma_T}{c}\int_{\rm los} dl\, n_e(r)\,[v(l)\cdot \hat n],
\]
and, in the ACT stacking analysis used in the observational paper, the gas is assumed to move with the halo bulk velocity and the signal is modeled using the RMS line-of-sight velocity of the sample, so that
\[
\frac{\Delta T_{\rm kSZ}}{T_{\rm CMB}}\approx \frac{\sigma_T}{c}\,v_r\int_{\rm los} dl\, n_e(l),
\]
for \(\tau\ll 1\) in the relevant redshift range [2507.07991]. The optical depth profile is
\[
\tau(\theta)=\sigma_T\int_{\rm los} dl\, n_e\!\left(\sqrt{l^2+d_A(z)^2\theta^2}\right),
\]
and the ACT modeling pipeline computes \(n_e\) from \(\rho_{\rm gas}\), projects to 2D, convolves with Gaussian beams with FWHM \(2.1'\) at 98 GHz and \(1.3'\) at 150 GHz, applies the compensated aperture-photometry filter \(W_{\theta_d}\), and then computes \(\mathcal{T}(\theta_d)=\int d^2\theta\,\Delta T(\theta)W_{\theta_d}(\theta)\) [2507.07991].

For X-ray gas fractions, the quantity within overdensity radius \(R_\Delta\) is
\[
f_{{\rm gas},\Delta}(M,z)=\frac{\int_0^{R_\Delta}dr\,r^2\rho_{\rm gas}(r)}{\int_0^{R_\Delta}dr\,r^2\rho_{\rm tot}(r)},
\]
with
\[
\rho_{\rm tot}(r)=\rho_{\rm gas}+\rho_{\rm cga}+\rho_{\rm sga}+\rho_{\rm dm}+\rho_{2{\rm h}}.
\]
Feedback modifies \(\rho_{\rm gas}\) and the dark-matter response, changing \(f_{{\rm gas},500}(M,z)\), and the model predicts the full \(f_{\rm gas}\)–mass relation [2507.07991].

The same baryonified fields are used to define matter power-spectrum suppression,
\[
S(k,z)\equiv \frac{P_{\rm bary}(k,z)}{P_{\rm DMO}(k,z)},
\]
where \(P_{\rm bary}(k,z)=\langle |\delta_{m,{\rm bary}}(k,z)|^2\rangle\) and \(P_{\rm DMO}(k,z)=\langle |\delta_{m,{\rm DMO}}(k,z)|^2\rangle\) [2507.07991]. By construction, \(S(k\rightarrow 0)\rightarrow 1\); the suppression grows with \(k\) as gas is pushed to larger radii and stars concentrate at the center [2507.07991].

The main observational constraint to date combines ACT DR5 kSZ stacked profiles following Schaan et al. (2021) with halo X-ray gas fractions from eROSITA, specifically eFEDS \(\times\) GAMA [2507.07991]. The ACT observable is the velocity-reconstruction-weighted stacked kSZ temperature profile at 98 and 150 GHz for BOSS CMASS galaxies, over redshift range \(z\approx 0.4\)–0.6 with mean \(\bar z\approx 0.55\), halo mass range \(M_{200}\sim10^{13}\)–\(10^{14}\,M_\odot\), and mean \(\bar M_{200}\approx 3\times10^{13}\,M_\odot\) [2507.07991]. The eROSITA observable is stacked X-ray-inferred hot-gas fractions in groups and clusters over \(M_{500}\sim10^{12}\)–\(10^{15}\,M_\odot\), with the analysis excluding \(M_{500}<10^{13}\,M_\odot\) due to uncertain mass calibration [2507.07991].

A joint fit to ACT kSZ and eROSITA \(f_{\rm gas}\) is described as excellent, with \(\chi^2=24.57\) and \({\rm PTE}=0.41\); the kSZ-only and \(f_{\rm gas}\)-only fits at the joint best-fit point have PTEs of \(0.22\) and \(0.82\), respectively [2507.07991]. The kSZ data strongly constrain \(\delta\), while eROSITA \(f_{\rm gas}\) tightens \(M_c\) and \(\mu\); the combined data prefer stronger-than-standard feedback, with flattened inner gas profiles and steeper outskirts [2507.07991]. Even when left free, the best-fit \(\bar M_{200,{\rm ksz}}\) matches independent estimates for CMASS halos, specifically the McCarthy et al. (2024) estimate \(\log_{10}M_{500{\rm c}}=13.34\pm0.04\), implying \(\log_{10}M_{200{\rm c}}\approx 13.52\) [2507.07991].

Quantitatively, the joint fit predicts a matter power suppression that exceeds the percent level for \(k\gtrsim 0.3\)–\(0.6\,h\,{\rm Mpc}^{-1}\), reaches 2–8 percent at \(k=1\,h\,{\rm Mpc}^{-1}\), and 20–25 percent at \(k=5\,h\,{\rm Mpc}^{-1}\), closely tracking the strong-feedback FLAMINGO \(f_{\rm gas}\)-8\(\sigma\) model and lying well below the fiducial FLAMINGO m8 run [2507.07991]. Using the kSZ+eROSITA best-fit parameters, the model reproduces the observed X-COP cluster electron density and pressure profiles, with mean \(M_{200}\approx 9.7\times10^{14}\,M_\odot\) at \(z\approx 0.065\), within \(1\sigma\) across radii \(\sim 0.02\)–\(1.5\,R_{200}\), despite not fitting X-COP [2507.07991].

## 6. Extensions, limitations, and current scope

BFC has been extended beyond clustering, X-ray, and SZ applications to a fully analytical framework for predicting the one-point probability distribution function of dispersion measures for FRBs [2601.18784]. In that application, the dominant ionized component is the hot circumgalactic or collapsed gas, and the hot-gas profile parameters \(M_c\), \(\mu\), and \(\delta\) are identified as the primary drivers of the PDF shape [2601.18784]. The analytical FRB model uses the halo mass function and halo bias to convolve contributions from individual halos across a range of masses and redshifts, includes halo clustering via linear bias and Gaussian averaging, and compares favorably with IllustrisTNG across redshifts up to \(z=5\) [2601.18784]. The paper further finds that the log-normal approximation commonly used for DM distributions provides a sufficient description for a few hundred FRBs, while departures can become important as samples grow larger [2601.18784].

The framework’s limitations are stated explicitly. Small-scale gas physics such as metal-dependent cooling, detailed multiphase structure, cold clumps, and feedback mode switching are not explicitly modeled; the inner gas profile \(\rho_{\rm iga}\) is phenomenological and mainly improves small-\(r\) upturns at higher redshift [2507.07892]. Non-thermal pressure and hydrostatic bias are treated with an empirical correction, and deviations from equilibrium and halo-to-halo scatter are not modeled beyond this [2507.07892]. Dark-matter back-reaction is calibrated to several suites, but \(Q_2\) shows resolution dependence, and at high redshift the smallest-scale upturn in \(S(k)\) is slightly underpredicted [2507.07892]. Low-mass halos outside the calibration bins contribute more to \(S(k,z)\) at higher redshift, so extrapolation can degrade small-scale accuracy [2507.07892].

The observational program has also exposed a substantive dataset-level tension. Joint fits to ACT kSZ and pre-eROSITA gas fractions from Lovisari et al. (2015), Gonzalez et al. (2013), and Sun et al. (2009) are acceptable but exhibit clear parameter tension: the X-ray-only fit prefers higher gas fractions at group masses, and combining with kSZ forces extreme parameter values and a very steep downturn below the lowest X-ray mass points [2507.07991]. The likely sources of tension identified in the paper are selection biases in pre-eROSITA samples, which favor X-ray-bright, gas-rich systems, and hydrostatic mass bias and calibration systematics in total mass estimates [2507.07991]. By contrast, eROSITA uses stacked profiles of optically selected groups with eROSITA X-rays and produces lower and more representative \(f_{\rm gas}\) at group masses that align with ACT kSZ and with a strong-feedback scenario [2507.07991].

Within those stated assumptions, BFC occupies a specific methodological niche. It is presented as an alternative to hydrodynamical simulations for cosmological studies, with multi-probe readiness for weak lensing, galaxy clustering, X-ray, SZ, and FRB statistics, and with physically interpretable parameters linked to halo-scale gas and stellar content and feedback strength [2507.07892] [2601.18784]. A plausible implication is that its main value lies in carrying baryonic structure, thermodynamics, and dark-matter response through a single forward model across observables and across the mass scales relevant for Stage-IV survey analyses.

Source: https://www.emergentmind.com/topics/component-wise-baryonification-bfc