pith. sign in

arxiv: 2606.19428 · v1 · pith:GCLBVB3Vnew · submitted 2026-06-17 · 🌌 astro-ph.CO · astro-ph.GA· hep-ph

Cooling, conduction, compact objects: Gravothermal evolution of dissipative self-interacting dark matter halos

Pith reviewed 2026-06-26 19:35 UTC · model grok-4.3

classification 🌌 astro-ph.CO astro-ph.GAhep-ph
keywords self-interacting dark matterdissipative SIDMgravothermal evolutiondark matter halosstrong gravitational lensingN-body simulationsheat conduction
0
0 comments X

The pith

Strong central cooling in dissipative SIDM halos inverts heat conduction and keeps it directed inward throughout evolution.

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

The paper investigates how radiative dissipation from self-interactions modifies the gravothermal evolution of isolated dark matter halos. It finds that sufficient cooling suppresses the usual isothermal core, so conduction flows inward instead of outward and outer regions cool rather than heat. Density profiles then show a broader infall region and reduced indentation between core and halo. Weak dissipation also reproduces the properties of the observed strong lens perturber in JVAS B1938+666 with shorter times or smaller cross sections than the elastic case. The work extends N-body methods for frequent small-angle interactions to include dissipation and validates results against a fluid model.

Core claim

Dissipation qualitatively changes gravothermal evolution beyond accelerating collapse. Sufficiently strong central cooling inverts the usual role of heat conduction by suppressing isothermal core formation, so conduction remains directed inward throughout. Outer halo regions beyond the scale radius cool efficiently rather than being heated, producing a larger region of mass infall and a less pronounced indentation in the final density profile. These effects depend strongly on the cooling rate but are insensitive to the angular dependence of the cross section. Weakly dissipative self-interactions explain the JVAS B1938+666 perturber with significantly shorter evolution times or smaller cross

What carries the argument

The first extension of the N-body formalism for frequent small-angle self-interactions (fSIDM) to include effective dissipation, compared against a dissipative gravothermal fluid model.

If this is right

  • Outer halo regions cool efficiently, increasing the scale of mass infall.
  • Density profiles show a less pronounced indentation between the core and outer halo.
  • Conduction stays inward for the entire evolution when cooling is strong enough.
  • Observed compact objects can be matched with smaller self-interaction cross sections or shorter times.

Where Pith is reading between the lines

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

  • Dissipative models may relax upper bounds on elastic cross sections while still fitting halo data.
  • The same cooling mechanism could affect the structure of other compact dark-matter objects beyond the reported lens.
  • Varying the cooling rate in future runs could map the transition between standard and inverted conduction regimes.

Load-bearing premise

The extended N-body formalism with dissipation, together with the fluid model comparison, accurately captures radiative processes without major artifacts that would alter the reported qualitative changes in halo evolution.

What would settle it

A density profile from simulation or observation that still develops a clear isothermal core despite strong central cooling rates would falsify the claimed inversion of conduction direction.

Figures

Figures reproduced from arXiv: 2606.19428 by Ludwig D. Schmidt, Mathias Garny, Moritz S. Fischer.

Figure 1
Figure 1. Figure 1: Left: Schematic scattering of numerical particles. Both decel￾erate by ∆vdrag and experience opposite ∆vrand in a random direction perpendicular to the relative velocity ∆vi j. Right: Zoom-in to the level of physical particles for a dissipative process like Eq. (6). The physical particles undergo many such scatterings in a single numerical interac￾tion. Therefore, the radiated momenta k average out and the… view at source ↗
Figure 3
Figure 3. Figure 3: Core formation times (circles) and collapse times (squares) as a function of the dissipation strength σT/mχ(rdiss − 1) at fixed conduc￾tion strength σTrdiss ≃ 50 cm2 g −1 . The colors are consistent with [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 2
Figure 2. Figure 2: Time evolution of the central density (top), velocity dispersion (middle), and total dissipated energy (bottom) of an isolated halo with different initial dissipation rates. Solid lines (foreground) show the N￾body results, while faint lines (background) indicate the corresponding fluid-model predictions. The dashed horizontal line marks five times the initial central density, which is used to define the c… view at source ↗
Figure 4
Figure 4. Figure 4: Upper row: Density profiles at t = tcore, t = 0.9tcoll, and t = tcoll for different dissipation strengths. The initial NFW halo is shown in gray for reference, and corresponding fluid-model profiles are shown as faint lines. Lower row: Corresponding velocity dispersion profiles. and exhibits a shallower outer slope for strong dissipation. We discuss the underlying reason for this behavior in Sect. 4. The f… view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of central density (top), velocity dispersion (mid￾dle), and total dissipated energy (bottom) as in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Density (upper row) and velocity dispersion profiles (lower row) as in [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Specific power, i.e., energy change per mass bin, due to dissipation ˙qdiss, conduction ˙qcond, and mechanical work ˙w, together with their sum, the specific internal energy change ˙u, as a function of radius. The columns show different times relative to the core formation time tcore and collapse time tcoll. The upper row was obtained from an elastic fluid model run with σT/mχ = 50 cm2 g −1 , the lower row… view at source ↗
Figure 9
Figure 9. Figure 9: Evolutionary tracks of all simulations in the plane of projected enclosed mass and the logarithmic slope of the surface density profile. The upper row shows the simulations with varying dissipation (Sect. 3.2), while the lower row shows the simulations with varying heat conduction (Sect. 3.3). Colors are consistent with the legends in these sections, with faint background lines representing the fluid model… view at source ↗
Figure 10
Figure 10. Figure 10: Ratio of projected enclosed masses M(r < r2D)/M(r < 4.5r2D) as a function of the normalized inner radius r2D/rs for our elastic (left) and least dissipative (right) simulation. Darker colors correspond to later snapshots. The dashed horizontal line indicates the mass ratio of the exotic compact object in JVAS B1938+666 inferred by Vegetti et al. (2026), with the shaded region showing the 3σ confidence int… view at source ↗
read the original abstract

Many proposed self-interacting dark matter (SIDM) models give rise to radiative processes that can dissipate energy. Understanding their impact on astrophysical objects through simulations and comparing the results with observations may thus constrain SIDM models. In this work, we systematically investigate how dissipation alters the gravothermal evolution of isolated SIDM halos by independently varying dissipation and heat conduction and identify potential observational signatures. To this end, we present the first extension of the $N$-body formalism for frequent small-angle self-interactions (fSIDM) to include effective dissipation. We compare all results for isolated halos with a dissipative gravothermal fluid model to assess its validity and limitations. We find that dissipation qualitatively changes the gravothermal evolution of SIDM halos beyond simply accelerating collapse. Sufficiently strong central cooling can invert the usual role of heat conduction: the formation of an isothermal core is suppressed such that conduction remains directed inward throughout the evolution. Outer halo regions beyond the scale radius can cool efficiently rather than being heated by conduction, resulting in a larger region of mass infall and a less pronounced indentation between the core and the outer halo in the final density profile. These effects depend strongly on the cooling rate but are comparatively insensitive to the angular dependence of the self-interaction cross section. We further show that weakly dissipative self-interactions can explain the properties of the recently observed strong lens perturber in JVAS~B1938+666 with significantly shorter evolution times or, equivalently, smaller cross sections compared to the elastic case. Our results open a new route to connecting halo structure and recently reported compact objects to dark-sector microphysics.

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.

Referee Report

2 major / 3 minor

Summary. The manuscript extends the N-body formalism for frequent small-angle self-interacting dark matter (fSIDM) to include effective dissipation, compares results for isolated halos against a dissipative gravothermal fluid model, and reports that sufficiently strong central cooling inverts the usual role of heat conduction: isothermal core formation is suppressed, conduction remains inward-directed, outer regions cool efficiently, and the final density profile shows a larger mass-infall region with reduced core-outer indentation. These effects depend strongly on cooling rate but weakly on angular dependence of the cross section. The work further claims that weakly dissipative self-interactions can reproduce the observed properties of the strong-lens perturber in JVAS B1938+666 with shorter evolution times or smaller cross sections than the elastic case.

Significance. If the reported inversion and its mapping to the JVAS B1938+666 perturber hold, the results supply a concrete route from dark-sector radiative processes to observable halo structure and compact objects, extending beyond elastic SIDM. The dual-method comparison (N-body extension plus fluid model) and the identification of qualitative changes that are not merely accelerated collapse constitute clear strengths.

major comments (2)
  1. [results section (N-body vs. fluid comparison)] The central claim that conduction remains directed inward throughout the evolution (abstract and strongest claim) is load-bearing for the inversion result; the manuscript must show explicit time series of the conductive heat flux or temperature gradient (e.g., in the results section comparing N-body and fluid runs) to demonstrate that the inversion is sustained rather than transient.
  2. [discussion of observational application] The quantitative statement that weakly dissipative interactions explain the JVAS B1938+666 perturber with 'significantly shorter evolution times or equivalently smaller cross sections' requires tabulated values of the adopted cooling rate, cross section, and resulting collapse timescale for both dissipative and elastic cases so that the magnitude of the reduction can be assessed directly.
minor comments (3)
  1. [methods] Notation for the effective dissipation term in the extended fSIDM N-body scheme should be defined explicitly (e.g., the relation between the dissipation parameter and the radiative cooling function) to allow direct comparison with the fluid-model equations.
  2. [figures] Figure captions for the density-profile evolution plots should state the number of independent realizations and the resolution used, to clarify the robustness of the reported qualitative changes in core indentation.
  3. [abstract] The abstract states that effects are 'comparatively insensitive to the angular dependence'; a brief sentence in the results summarizing the tested angular regimes (e.g., isotropic vs. forward-peaked) would make this statement self-contained.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment and constructive comments. We address the two major comments below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [results section (N-body vs. fluid comparison)] The central claim that conduction remains directed inward throughout the evolution (abstract and strongest claim) is load-bearing for the inversion result; the manuscript must show explicit time series of the conductive heat flux or temperature gradient (e.g., in the results section comparing N-body and fluid runs) to demonstrate that the inversion is sustained rather than transient.

    Authors: We agree that explicit time series are needed to confirm the inversion is sustained. We will add plots of the conductive heat flux and temperature gradient versus time (and radius) in the results section for representative N-body and fluid runs, demonstrating that inward conduction persists for strong cooling cases rather than being transient. revision: yes

  2. Referee: [discussion of observational application] The quantitative statement that weakly dissipative interactions explain the JVAS B1938+666 perturber with 'significantly shorter evolution times or equivalently smaller cross sections' requires tabulated values of the adopted cooling rate, cross section, and resulting collapse timescale for both dissipative and elastic cases so that the magnitude of the reduction can be assessed directly.

    Authors: We will add a table in the observational discussion section listing the cooling rates, velocity-dependent cross sections, and resulting collapse timescales for both the dissipative and elastic cases applied to the JVAS B1938+666 perturber, allowing direct comparison of the reduction factors. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper extends N-body methods for dissipative fSIDM and compares results to an independent gravothermal fluid model. Central claims (conduction inversion under strong cooling, shorter evolution times for JVAS B1938+666) are outputs of these direct simulations and model comparisons rather than quantities defined by the paper's own fitted parameters or self-citation chains. No load-bearing step reduces by construction to its inputs; the derivation remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, specific free parameters, axioms, and invented entities cannot be identified. The work presumably relies on standard SIDM modeling assumptions such as the validity of the fSIDM small-angle approximation and the effective treatment of dissipation as an added cooling term.

pith-pipeline@v0.9.1-grok · 5837 in / 1263 out tokens · 29472 ms · 2026-06-26T19:35:43.740853+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

72 extracted references · 5 canonical work pages

  1. [1]

    doi:10.48550/arXiv.2603.19362 , archivePrefix =

    Survival of the most compact: the life and death of satellite halos in self-interacting dark matter. doi:10.48550/arXiv.2603.19362 , archivePrefix =. 2603.19362 , primaryClass =

  2. [2]

    2026 , month = may, number =

    , title =. 2026 , month = may, number =. arXiv , doi =:2504.13004 , keywords =

  3. [3]

    2026 , month = apr, number =

    , title =. 2026 , month = apr, number =. arXiv , doi =:2510.11006 , keywords =

  4. [4]

    2026 , month = mar, pages =

    Nature Astronomy , title =. 2026 , month = mar, pages =

  5. [5]

    2026 , month = mar, number =

    , title =. 2026 , month = mar, number =. arXiv , doi =:2509.12317 , keywords =

  6. [6]

    arXiv , doi =:2603.10107 , keywords =

    2026 , month = mar, adsnote =. arXiv , doi =:2603.10107 , keywords =

  7. [7]

    2025 , month = nov, pages =

    , title =. 2025 , month = nov, pages =. arXiv , doi =:2506.06269 , keywords =

  8. [8]

    2025 , month = nov, pages =

    Nature Astronomy , title =. 2025 , month = nov, pages =. arXiv , doi =:2510.07382 , keywords =

  9. [9]

    2025 , month = oct, number =

    Reviews of Modern Physics , title =. 2025 , month = oct, number =. arXiv , doi =:2207.10638 , keywords =

  10. [10]

    2025 , month = aug, pages =

    , title =. 2025 , month = aug, pages =. arXiv , doi =:2504.12393 , keywords =

  11. [11]

    2025 , month = apr, number =

    , title =. 2025 , month = apr, number =. arXiv , doi =:2408.15317 , keywords =

  12. [12]

    2025 , month = feb, pages =

    , title =. 2025 , month = feb, pages =. arXiv , doi =:2410.07175 , keywords =

  13. [13]

    2024 , month = dec, number =

    , title =. 2024 , month = dec, number =. arXiv , doi =:2402.01604 , keywords =

  14. [14]

    2024 , month = may, number =

    , title =. 2024 , month = may, number =. arXiv , doi =:2206.05327 , keywords =

  15. [15]

    2024 , month = may, number =

    , title =. 2024 , month = may, number =. arXiv , doi =:2312.09296 , keywords =

  16. [16]

    2024 , month = apr, number =

    , title =. 2024 , month = apr, number =. arXiv , doi =:2310.07750 , keywords =

  17. [17]

    2023 , month = sep, number =

    , title =. 2023 , month = sep, number =. arXiv , doi =:2304.09878 , keywords =

  18. [18]

    2023 , month = aug, number =

    , title =. 2023 , month = aug, number =. arXiv , doi =:2204.06568 , keywords =

  19. [19]

    2022 , month = nov, number =

    , title =. 2022 , month = nov, number =. arXiv , doi =:2209.00064 , keywords =

  20. [20]

    2022 , month = sep, number =

    , title =. 2022 , month = sep, number =. arXiv , doi =:2205.03392 , keywords =

  21. [21]

    2022 , month = feb, number =

    , title =. 2022 , month = feb, number =. arXiv , doi =:2012.06611 , keywords =

  22. [22]

    2021 , month = sep, number =

    , title =. 2021 , month = sep, number =. arXiv , doi =:2102.09580 , keywords =

  23. [23]

    2021 , month = jul, number =

    , title =. 2021 , month = jul, number =. arXiv , doi =:2012.10277 , keywords =

  24. [24]

    2021 , month = jul, number =

    , title =. 2021 , month = jul, number =. arXiv , doi =:2103.13407 , keywords =

  25. [25]

    2021 , month = jan, number =

    , title =. 2021 , month = jan, number =. arXiv , doi =:2006.12515 , keywords =

  26. [26]

    2020 , month = sep, pages =

    , title =. 2020 , month = sep, pages =. arXiv , doi =:1807.06209 , keywords =

  27. [27]

    2020 , month = jun, number =

    , title =. 2020 , month = jun, number =. arXiv , doi =:1912.06757 , keywords =

  28. [28]

    2020 , month = mar, number =

    , title =. 2020 , month = mar, number =. arXiv , doi =:1901.00499 , keywords =

  29. [29]

    2019 , month = sep, number =

    , title =. 2019 , month = sep, number =. arXiv , doi =:1809.01144 , keywords =

  30. [30]

    2019 , month = jul, number =

    , title =. 2019 , month = jul, number =. arXiv , doi =:1812.05088 , keywords =

  31. [31]

    2019 , month = feb, number =

    , title =. 2019 , month = feb, number =. arXiv , doi =:1809.07326 , keywords =

  32. [32]

    2018 , month = dec, number =

    , title =. 2018 , month = dec, number =. arXiv , doi =:1712.04512 , keywords =

  33. [33]

    2018 , month = nov, number =

    , title =. 2018 , month = nov, number =. arXiv , doi =:1805.04512 , keywords =

  34. [34]

    2018 , month = feb, pages =

    , title =. 2018 , month = feb, pages =. arXiv , doi =:1705.02358 , keywords =

  35. [35]

    2017 , month = aug, number =

    , title =. 2017 , month = aug, number =. arXiv , doi =:1707.04256 , keywords =

  36. [36]

    2017 , month = jun, number =

    , title =. 2017 , month = jun, number =. arXiv , doi =:1612.03906 , keywords =

  37. [37]

    2016 , month = aug, number =

    International Journal of Modern Physics A , title =. 2016 , month = aug, number =. arXiv , doi =:1604.04627 , keywords =

  38. [38]

    2016 , month = may, pages =

    Advances in Parallel Computing , title =. 2016 , month = may, pages =. arXiv , doi =:1810.09898 , keywords =

  39. [39]

    2016 , month = jan, pages =

    Saas-Fee Advanced Course , title =. 2016 , month = jan, pages =. arXiv , doi =:1412.5187 , keywords =

  40. [40]

    2015 , month = oct, number =

    , title =. 2015 , month = oct, number =. arXiv , doi =:1412.1477 , keywords =

  41. [41]

    2015 , month = sep, number =

    , title =. 2015 , month = sep, number =. arXiv , doi =:1412.1839 , keywords =

  42. [42]

    2014 , month = sep, number =

    , title =. 2014 , month = sep, number =. arXiv , doi =:1407.4121 , keywords =

  43. [43]

    2014 , month = jul, number =

    , title =. 2014 , month = jul, number =. arXiv , doi =:1402.7073 , keywords =

  44. [44]

    2014 , month = jan, number =

    , title =. 2014 , month = jan, number =. arXiv , doi =:1308.3419 , keywords =

  45. [45]

    2014 , month = jan, number =

    , title =. 2014 , month = jan, number =. arXiv , doi =:1311.2937 , keywords =

  46. [46]

    2013 , month = sep, number =

    Physics of the Dark Universe , title =. 2013 , month = sep, number =. arXiv , doi =:1306.1305 , keywords =

  47. [47]

    2013 , month = sep, number =

    Physics of the Dark Universe , title =. 2013 , month = sep, number =. arXiv , doi =:1303.1521 , keywords =

  48. [48]

    2013 , month = aug, number =

    , title =. 2013 , month = aug, number =. arXiv , doi =:1301.3137 , keywords =

  49. [49]

    2013 , month = jun, number =

    , title =. 2013 , month = jun, number =. arXiv , doi =:1302.3898 , keywords =

  50. [50]

    2012 , month = sep, number =

    , title =. 2012 , month = sep, number =. arXiv , doi =:1205.2729 , keywords =

  51. [51]

    2011 , month = aug, number =

    , title =. 2011 , month = aug, number =. arXiv , doi =:1101.3097 , keywords =

  52. [52]

    2010 , month = may, number =

    , title =. 2010 , month = may, number =. arXiv , doi =:0909.0753 , keywords =

  53. [53]

    2008 , month = jun, number =

    , title =. 2008 , month = jun, number =. arXiv , doi =:0704.0261 , keywords =

  54. [54]

    2005 , month = dec, number =

    , title =. 2005 , month = dec, number =. arXiv , doi =:astro-ph/0505010 , keywords =

  55. [55]

    2004 , month = jan, number =

    International Journal of Modern Physics D , title =. 2004 , month = jan, number =. arXiv , doi =:astro-ph/0407623 , keywords =

  56. [56]

    2003 , month = jan, number =

    , title =. 2003 , month = jan, number =. arXiv , doi =:astro-ph/0201544 , keywords =

  57. [57]

    2002 , month = dec, number =

    , title =. 2002 , month = dec, number =. arXiv , doi =:astro-ph/0205322 , keywords =

  58. [58]

    2002 , month = apr, number =

    , title =. 2002 , month = apr, number =. arXiv , doi =:astro-ph/0110561 , keywords =

  59. [59]

    2001 , month = nov, number =

    , title =. 2001 , month = nov, number =. arXiv , doi =:astro-ph/0010436 , keywords =

  60. [60]

    2001 , month = feb, number =

    , title =. 2001 , month = feb, number =. arXiv , doi =:astro-ph/0006218 , keywords =

  61. [61]

    2000 , month = nov, number =

    , title =. 2000 , month = nov, number =. arXiv , doi =:astro-ph/0003483 , keywords =

  62. [62]

    2000 , month = may, number =

    , title =. 2000 , month = may, number =. arXiv , doi =:astro-ph/0002409 , keywords =

  63. [63]

    2000 , month = apr, number =

    , title =. 2000 , month = apr, number =. arXiv , doi =:astro-ph/9909386 , keywords =

  64. [64]

    1997 , month = dec, number =

    , title =. 1997 , month = dec, number =. arXiv , doi =:astro-ph/9611107 , keywords =

  65. [65]
  66. [66]

    1968 , month = jan, pages =

    , title =. 1968 , month = jan, pages =

  67. [67]

    , keywords =

    Simulations of galaxy cluster mergers with velocity-dependent, rare, and frequent self-interactions. , keywords =. doi:10.1093/mnras/stae664 , archivePrefix =. 2310.07769 , primaryClass =

  68. [68]

    , keywords =

    On the consequences of the gravothermal catastrophe. , keywords =. doi:10.1093/mnras/191.3.483 , adsurl =

  69. [69]

    , keywords =

    Anisotropy in stellar dynamics. , keywords =

  70. [70]

    SIDM and CDM interpretations of the million-solar-mass lensing perturber JVAS B1938+666- V

  71. [71]

    R., Millman, K

    Charles R. Harris and K. Jarrod Millman and St. Array programming with. 2020 , month = sep, journal =. doi:10.1038/s41586-020-2649-2 , publisher =

  72. [72]

    Hunter, J. D. , Title =. Computing in Science & Engineering , Volume =