Testing the isothermal Jeans model for self-interacting dark matter halos in the collapse phase

arXiv:2603.01772 · astro-ph.CO, astro-ph.GA · Submitted 2026-03-02 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: We're looking at the paper: "Testing the isothermal Jeans model for self-interacting dark matter halos in the collapse phase".

Jocelyn: Here's the gist of the paper. The isothermal Jeans model provides a reliable and computationally efficient description of Self-Interacting Dark Matter (SIDM) halo evolution, successfully benchmarking its predictions against high-resolution N-body simulations across different stages of collapse.

Vera: The gist is that the isothermal Jeans model offers a fast and broadly reliable framework for modeling the gravothermal evolution of SIDM halos, including the core-collapse phase <ref:2603.01772#pg11>.

Jocelyn: How it works The isothermal Jeans model is built on the assumption that self-interactions in the dense, central region of a SIDM halo lead to an isothermal core determined by solving coupled Jeans and Poisson equations <ref:2603.01772#pg4>.

Vera: The density profile of this inner core is uniquely defined by two parameters: its central density (ρ0) and its constant velocity dispersion (ν0) <ref:2603.01772#pg5>.

Title and authors: Vera: We started by looking at the title of "Testing the isothermal Jeans model for self-interacting dark matter halos in the collapse phase" and what that actually means for us. It's about testing how well this specific simple model performs when you look at SIDM halos during their collapse.

Jocelyn: Basically, it’s a direct comparison between the predictions of this semi-analytical model against those high-resolution Nbody simulations we use to describe what happens in these halos.

Subrahmanyan: It suggests that if you want to test SIDM predictions across all stages of halo evolution, you need a robust way to model the full gravothermal evolution, and this paper offers that framework.

Vera: The authors are trying to show that their Jeans model can handle the complexity of going from initial halo state all the way into deep core collapse.

Jocelyn: They are showing how they use Nbody simulations to benchmark this model against those full particle trajectories, which gives them confidence in their results.

Subrahmanyan: The big implication here is that we need these kinds of complementary approaches because Nbody simulations are too slow for systematic exploration of large parameter spaces and long timescales.

Vera: So, they aren't just proposing a single model; they are showing that this Jeans approach is one piece of a set needed to handle the whole picture.

Jocelyn: It’s about finding methods that can cover the full evolution without needing a super expensive calculation every single time.

Subrahmanyan: This connects directly to how we can constrain SIDM models, and if we want to test them, we need these fast tools alongside the slower simulations.

Vera: So the paper is positioning this isothermal Jeans model as a reliable framework for modeling SIDM halos across their entire lifecycle, from start to finish.

Jocelyn: It’s about showing that even with its simplifications, it still gives us useful physical insight into how these halos evolve gravitationally.

The paper's summary: Vera: Now let’s break down the summary they give of "Testing the isothermal Jeans model for self-interacting dark matter halos in the collapse phase" and what that actually means on the ground.

Jocelyn: The core idea is that self-interactions in a dense central region of a SIDM halo lead to an isothermal core determined by solving coupled Jeans and Poisson equations.

Subrahmanyan: They assume this inner core structure is uniquely defined by two parameters: the central density rho0 and a constant velocity dispersion nu0.

Vera: Outside that core, where self-interactions are negligible, they assume the density follows a standard Navarro-Frenk-White NFW profile.

Jocelyn: This whole setup is smoothly joined at a stitching radius r1, which marks the boundary where the local scattering rate is high enough for a particle to undergo one interaction on average over time t.

Subrahmanyan: They also introduce an extension into the core-collapse regime using mirroring prescriptions proposed by S. Yang et al., which exploits two branches of matched isothermal solutions at a merger time tmerge.

Vera: This mirrors the high-density branch about tmerge to provide a continuous mapping onto the subsequent collapse evolution, enabling them to construct SIDM halo profiles up to the latest collapse stages accessible within this framework.

Jocelyn: It’s essentially building a mathematical bridge across different evolutionary phases of the halo structure.

The paper's improvements: Vera: The paper points out some specific ways they improved their approach to make it more robust, and we need to talk about those so we can see how they made this model better.

Jocelyn: They used the mirroring prescription from S. Yang et al., which is a key improvement because it lets them exploit two branches of matched isothermal solutions at a merger time tmerge.

Subrahmanyan: This specific technique allows for that continuous mapping onto the subsequent collapse evolution, which enables construction of SIDM halo profiles up to the latest collapse stages accessible within this framework.

Vera: It’s about giving them that continuous mapping, which means they can map out more of the evolution than before.

Jocelyn: That gives them a better way to construct the density profile for those later stages, and it's less arbitrary than just picking an arbitrary point on a curve.

Conclusion: Vera: So to wrap up this paper, we’ve seen that the main conclusion is that this isothermal Jeans model provides a reliable description of SIDM halo evolution through core formation and into the early stages of collapse.

Subrahmanyan: The model successfully captures the core-formation phase but it does not fully capture that rapid rise in central velocity dispersion during the later collapse phase.

Shubo Li, Moritz S. Fischer, Zixiang Jia, Fangzhou Jiang, Ran Li

School of Physics and Astronomy, Beijing Normal University · National Astronomical Observatories, Chinese Academy of Sciences · School of Astronomy and Space Science, University of Chinese Academy of Sciences · Donostia International Physics Center (DIPC) · Department of Astronomy, Peking University · Kavli Institute for Astronomy and Astrophysics

astro-ph.CO, astro-ph.GA

Submitted: 2026-03-02

Updated: 2026-10-18

Comments: 12 pages, 6 figures. Published in The Astrophysical Journal

DOI: 10.3847/1538-4357/aea78a

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 78/100

The gist: The isothermal Jeans model provides a reliable and computationally efficient description of Self-Interacting Dark Matter (SIDM) halo evolution, successfully benchmarking its predictions against

Key concepts

Isothermal Jeans Model: This is a mathematical framework used to model the gravothermal evolution of dark matter halos. It assumes that self-interactions in the dense center create an isothermal core characterized by a constant central density and velocity dispersion, while outer regions follow a standard NFW profile.
Core-Collapse Phase: This refers to the stage in halo evolution where self-interactions drive the halo toward a dense, centrally concentrated state. The model attempts to describe this transition using coupled Jeans and Poisson equations, specifically looking at how the core density ($ ho_0$) and velocity dispersion ($ u_0$) change over time.
NFW Profile: The Navarro-Frenk-White profile is a standard mathematical description used for the outer regions of dark matter halos. It describes how the density of dark matter gradually decreases as you move farther away from the halo's center, typically following a specific power law.
Gravothermal Evolution: This describes the process where self-interactions cause dark matter halos to evolve dynamically over time, often leading to core formation. The Jeans model is used to track this evolution, showing how the halo structure changes from its initial state through core formation and into collapse.

Terminology

Benchmark — 1× in this paper · explained in 2 episodes
Core — 1× in this paper · explained in 2 episodes
Gravitational Lensing — 1× in this paper · explained in 2 episodes
Velocity Dispersion — 1× in this paper · explained in 2 episodes

Summary

The isothermal Jeans model provides a reliable and computationally efficient description of Self-Interacting Dark Matter (SIDM) halo evolution, successfully benchmarking its predictions against high-resolution N-body simulations across different stages of collapse. The gist is that the isothermal Jeans model offers a fast and broadly reliable framework for modeling the gravothermal evolution of SIDM halos, including the core-collapse phase <ref:2603.01772#pg11>.

How it works

The isothermal Jeans model is built on the assumption that self-interactions in the dense, central region of a SIDM halo lead to an isothermal core determined by solving coupled Jeans and Poisson equations <ref:2603.01772#pg4>. The density profile of this inner core is uniquely defined by two parameters: its central density (ρ0) and its constant velocity dispersion (ν0) <ref:2603.01772#pg5>. In the outer regions, where self-interactions are negligible, the density follows a standard Navarro-Frenk-White (NFW) profile <ref:2603.01772#pg6>. The inner isothermal core and the outer NFW halo are smoothly joined at a stitching radius, r1, which marks the boundary where the local scattering rate is high enough for a particle to undergo one interaction on average over time <ref:2603.01772#pg6>.

To extend this framework into the core-collapse regime, the paper adopts the mirroring prescription proposed by S. Yang et al. (2024), which exploits two branches of matched isothermal solutions: a low-density branch and a high-density branch that approach each other at a merger time, tmerge <ref:2603.01772#pg4>. The high-density branch is then mirrored about tmerge to provide a continuous mapping onto the subsequent collapse evolution, enabling the construction of SIDM halo profiles up to the latest collapse stages accessible within this framework <ref:2603.01772#pg5>.

The simulation and comparison methodology

The study uses a state-of-the-art, high-resolution isolated N-body simulation from M. S. Fischer et al. (2025) as the reference for validation <ref:2603.01772#pg6>. This simulation is based on an initial halo density profile described by an NFW profile with a characteristic density ρs = 4.42 × 107 M⊙ kpc−3 and a scale radius rs = 1.28 kpc <ref:2603.01772#pg6>. The simulation assumes isotropic, velocity-independent scattering with a total cross section per physical dark matter particle mass of σm = 80 cm2g−1 and a gravitational softening length of ϵ = 1.0 pc <ref:2603.01772#pg6>.

The comparison involves timeresolved comparisons of three-dimensional density and velocity-dispersion profiles, as well as projected lensing observables <ref:2603.01772#pg6>. The Jeans model is tested by tracking the halo’s evolution from 0.001 Gyr up to the formal endpoint of the Jeans model, tcoll ≡ 2 × tmerge <ref:2603.01772#pg5>. The central density ρ0 and velocity dispersion ν0 are sampled at 1500 logarithmically spaced time steps to construct interpolation functions, ρ0(t) and ν0(t). The radial velocity-dispersion squared, σ2r(r), is derived by solving the spherical Jeans equation assuming an isotropic velocity distribution by setting the anisotropy parameter β = 0.

Key findings in lensing observables

The comparison of projected lensing observables, specifically the projected logarithmic density slope, γ2D versus the enclosed mass, M2D, reveals a clear two-phase behavior. Starting from the initial NFW profile at t = 0, both the simulation and the Jeans model yield nearly identical results. The evolution exhibits a clear two-phase behavior where the projected slope rapidly flattens while the enclosed mass decreases, reflecting core-formation. Following this stage, the projected slope steepens and the enclosed mass increases, marking the transition into gravothermal core collapse.

The evolutionary tracks from both methods agree well, with two notable discrepancies. The first appears at the point of maximum core size where the Jeans model tends to predict a more extended core than the simulation. The second discrepancy arises during the core-collapse phase when the Jeans model fails to achieve the same level of central concentration as in the simulation, leading to a lower enclosed mass and a shallower projected density slope in inner regions. These differences diminish with increasing projected radius since self-interactions play a significant role only in the halo’s central region.

Internal structure and velocity dispersion discrepancies

When comparing the underlying three-dimensional density profiles, the core central density, ρ0, is used as a reference quantity for comparison. Over most of the evolutionary history, the density profiles predicted by both methods agree remarkably well. However, noticeable differences emerge only in the late stages of core collapse when ρ0 ≳ 200 ρs. In this regime, the halo predicted by the Jeans model develops a less concentrated central region than the simulated halo, manifesting as a shallower projected logarithmic slope and a lower enclosed mass.

The velocity-dispersion profiles show good agreement up to approximately 4 Gyr. At later times, the Jeans-model predictions increasingly deviate from the simulation, leading to a systematic underestimation of the central velocity dispersion. The largest discrepancy, reaching up to ∼ 10 km s−1, occurs at the final simulation snapshot. This stems from a limitation where the region inside the stitching radius r1 is enforced to follow a single isothermal solution smoothly matched to an outer NFW profile, causing it to overestimate the size of the isothermal region and underestimate the true central velocity dispersion ν0.

Timescale interpretation and model limitations

The timescale for core collapse appears somewhat longer in the Jeans model than in the isolated simulation. This slower evolution is consistent with cosmological Pippin SIDM simulations and reflects the NFW outer boundary adopted in the Jeans framework. The discrepancy between isolated N-body and isothermal Jeans models should not be interpreted as a failure of the semi-analytic approach but rather as a reflection of the more rapid gravothermal evolution expected in isolated systems where environmental effects like mergers are absent.

The parametric gravothermal framework, developed by D. Yang et al. (2024, 2025), provides a complementary fast approach that describes the full gravothermal evolution with an analytic density profile whose parameters are given by fitting functions calibrated to high-resolution SIDM simulations. While the parametric model shows good agreement at small projected radii, deviations become more significant at larger projected radii during the deep core-collapse phase where the analytic parameterization lacks flexibility. The isothermal Jeans model provides a fast and broadly reliable framework for modeling SIDM halos, including the collapse phase.

Conclusion

In summary, the isothermal Jeans model provides a reliable description of SIDM halo evolution through core formation and into the early stages of collapse. The model accurately reproduces projected lensing observables and density profiles over most of the evolution, with deviations remaining within observational uncertainties at comparable mass scales. The model successfully captures the core-formation phase but does not fully capture the rapid rise in central velocity dispersion during the later collapse phase. Thus, it is an efficient and practical tool for modeling SIDM halos with broad potential applications.

How it works

The Jeans model provides a fast and broadly reliable framework for modeling the gravothermal evolution of SIDM halos, including the collapse phase. The model accurately reproduces the projected lensing observables of the simulation throughout the gravothermal evolution, with deviations remaining well within current observational uncertainties at comparable mass scales. The Jeans model reproduces three-dimensional density profiles of isolated simulation with excellent accuracy over most of the halo evolution while noticeable deviations emerge only in the deep core-collapse phase.

The Jeans model predicts a more extended core and a delayed collapse time relative to the isolated simulation, which is consistent with cosmological Pippin SIDM simulations and reflects the NFW outer boundary adopted in the Jeans framework. The model provides a reliable description of the velocity-dispersion profile through core formation and early stages of collapse although it does not fully capture the rapid rise in central velocity dispersion during later collapse phase.

Improvements for AI systems

Bold header: Enhanced SIDM Halo Evolution Modeling

The improved AI system can accurately track the full gravothermal evolution of SIDM halos, including the core-collapse phase by integrating state-of-the-art Nbody simulation results with the isothermal Jeans model, as demonstrated in Section 3.

Bold header: Robust Lensing Observable Prediction

The system will be capable of predicting projected logarithmic density slope, γ2D versus the enclosed mass, M2D with high fidelity across evolutionary stages, leveraging the comparison shown in Figure 1 to assess model deviations.

Bold header: Accurate Central Density and Velocity Dispersion Tracking

The improved AI can provide reliable predictions for the central density evolution as a function of time t, specifically identifying the minimum central density predicted by the Jeans model is approximately half that of the simulation and tracking velocity dispersion profiles to capture structural limitations like underestimation of the true central velocity dispersion ν0.

Bold header: Automated Cross-Section Inference Calibration

The system can be equipped to infer physical parameters, such as the self-interaction cross section, by analyzing the timescale discrepancy, potentially using the mapping described in Section 4.1 to relate simulation results to Jeans model predictions.