pith. sign in

arxiv: 2607.02151 · v1 · pith:DPH7UUFDnew · submitted 2026-07-02 · 🌌 astro-ph.CO

Spherically Symmetric Fluid Simulations of Black Hole Accretion in Self-Interacting Dark Matter Halos

Pith reviewed 2026-07-03 06:46 UTC · model grok-4.3

classification 🌌 astro-ph.CO
keywords black hole accretionself-interacting dark matterhydrodynamic simulationsspherically symmetricsingular isothermal sphereNavarro-Frenk-White profilethermal conductiondark matter halos
0
0 comments X

The pith

Black holes in SIS-like self-interacting dark matter halos grow from 100 to 10,000 solar masses in 2 million years.

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

The paper develops a one-dimensional spherically symmetric hydrodynamic code to simulate black hole accretion inside self-interacting dark matter halos treated as self-gravitating fluids with thermal conduction. It compares singular isothermal sphere and Navarro-Frenk-White density profiles and shows that accretion is set by the balance between gravity-driven inflow and heat transport from dark matter scattering. In the denser, SIS-like case this balance allows a 100 solar mass seed to reach 10,000 solar masses inside 2 million years. The same runs show that bigger seeds, steeper profiles, and larger scattering cross sections all speed growth. The work supplies a fluid-dynamical description of how early black holes can grow quickly inside certain dark matter structures.

Core claim

Black hole growth in SIDM halos is regulated by the competition between gravity-driven inflow and SIDM heat transport. An SIS-like environment facilitates rapid accretion, allowing a 100 M⊙ seed to grow to 10^4 M⊙ within 2 Myr. Larger initial black hole masses, steeper density profiles, and higher scattering cross sections significantly enhance the accretion rate.

What carries the argument

The self-gravitating fluid model with thermal conduction simulated by a one-dimensional spherically symmetric hydrodynamic code that uses operator-splitting finite-volume methods to track the balance between inflow and heat transport.

If this is right

  • Larger initial black hole masses produce higher accretion rates.
  • Steeper initial density profiles increase the growth rate.
  • Higher SIDM scattering cross sections accelerate black hole growth.
  • The SIS profile yields faster growth than the NFW profile under otherwise identical conditions.

Where Pith is reading between the lines

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

  • If SIDM halos with SIS-like cores are common at high redshift, the mechanism could help explain the existence of massive black holes observed at early cosmic times.
  • Direct comparison of the predicted accretion rates against future 3D SIDM simulations would test the spherical-symmetry assumption.
  • Observational upper limits on early black hole masses could be used to place bounds on the SIDM cross section once the model is extended to cosmological volumes.

Load-bearing premise

The self-gravitating fluid model with thermal conduction and the assumption of perfect spherical symmetry accurately capture the accretion dynamics in SIDM halos.

What would settle it

A three-dimensional simulation of the identical initial conditions that produces accretion rates differing by more than a factor of a few due to non-spherical flows or non-fluid particle behavior would falsify the central result.

Figures

Figures reproduced from arXiv: 2607.02151 by Bin Hu, Bocheng Zhu, Fan Zhou, Liang Gao, Rong-Gen Cai, Tan Chen, Zhe Meng.

Figure 1
Figure 1. Figure 1: FIG. 1: Comparison of our numerical results with those obtained from the gravothermal fluid code of Nishikawa et [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Simulation results for an SIS halo ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Comparison of the SIS halo evolution with and without heat flow. Parameters are identical to those in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Simulation results for an NFW halo ( [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Comparison of the NFW profile evolution with and without heat flow. Parameters are identical to those in [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Evolution of physical quantities under various parameter configurations. The baseline setup is an NFW [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Time evolution of the black hole mass (left panel) and accretion rate (right panel) under various parameter [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
read the original abstract

We investigate black hole accretion in self-interacting dark matter (SIDM) halos using a self-gravitating fluid model with thermal conduction. We develop a robust one-dimensional spherically symmetric hydrodynamic code based on an operator-splitting finite-volume method. Simulating both Singular Isothermal Sphere (SIS) and Navarro-Frenk-White (NFW) profiles, we find that black hole growth is regulated by the competition between gravity-driven inflow and SIDM heat transport. Our results demonstrate that an SIS-like environment facilitates rapid accretion, allowing a $100\,\mathrm{M_{\odot}}$ seed to grow to $10^4\,\mathrm{M_{\odot}}$ within $2\,\mathrm{Myr}$. Furthermore, we show that larger initial black hole masses, steeper density profiles, and higher scattering cross sections significantly enhance the accretion rate. This study provides a comprehensive fluid-dynamical picture of black hole growth in SIDM halos.

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

3 major / 1 minor

Summary. The paper develops a one-dimensional spherically symmetric hydrodynamic code based on an operator-splitting finite-volume method to simulate black hole accretion in self-interacting dark matter (SIDM) halos, modeled as self-gravitating fluids with thermal conduction. It compares Singular Isothermal Sphere (SIS) and Navarro-Frenk-White (NFW) profiles and reports that an SIS-like environment enables rapid accretion, allowing a 100 M⊙ seed to reach 10^4 M⊙ within 2 Myr, with accretion rates enhanced by larger initial black hole masses, steeper density profiles, and higher scattering cross sections.

Significance. If the numerical results are robust, the work supplies a fluid-dynamical picture of how SIDM heat transport competes with gravity-driven inflow to regulate black hole growth, potentially informing seed black hole evolution in alternative dark matter models. The explicit demonstration of profile-dependent growth rates is a concrete contribution, though its weight depends on validation of the 1D setup.

major comments (3)
  1. [Abstract] Abstract: the headline result (100 M⊙ seed growing to 10^4 M⊙ in 2 Myr under SIS-like conditions) is stated without any reported validation tests, convergence checks, error bars, or comparison to analytic limits, leaving the numerical support for the central claim unassessable from the provided information.
  2. [model development paragraph] Model development paragraph: the self-gravitating fluid model with thermal conduction is implemented under the assumption of perfect spherical symmetry, yet no justification or test is supplied for why net angular momentum or non-radial motions can be neglected in SIDM halos; a run with even a modest specific angular momentum floor would be required to show that centrifugal support does not throttle the reported inflow.
  3. [Abstract] Abstract and results description: the claim that larger initial masses, steeper profiles, and higher cross sections enhance accretion is presented as a simulation outcome, but without quantitative sensitivity analysis or resolution studies it is unclear whether the reported factor-of-100 growth is converged or sensitive to the operator-splitting treatment of conduction and self-gravity.
minor comments (1)
  1. Notation for the scattering cross section and conduction coefficient should be defined explicitly at first use to avoid ambiguity between particle-physics and fluid parameters.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive feedback, which highlights areas where the presentation of our numerical results can be strengthened. We address each major comment below and indicate the revisions we will make.

read point-by-point responses
  1. Referee: [Abstract] Abstract: the headline result (100 M⊙ seed growing to 10^4 M⊙ in 2 Myr under SIS-like conditions) is stated without any reported validation tests, convergence checks, error bars, or comparison to analytic limits, leaving the numerical support for the central claim unassessable from the provided information.

    Authors: We agree that the abstract and main text would benefit from explicit references to validation. In the revised manuscript we will add a dedicated methods subsection summarizing grid-convergence tests (doubling and halving the radial resolution), direct comparisons of the conduction-free limit to the analytic Bondi accretion rate, and quantitative error estimates on the integrated mass growth. These additions will make the numerical support for the reported growth rates assessable. revision: yes

  2. Referee: [model development paragraph] Model development paragraph: the self-gravitating fluid model with thermal conduction is implemented under the assumption of perfect spherical symmetry, yet no justification or test is supplied for why net angular momentum or non-radial motions can be neglected in SIDM halos; a run with even a modest specific angular momentum floor would be required to show that centrifugal support does not throttle the reported inflow.

    Authors: The one-dimensional formulation is adopted to isolate the radial competition between gravity-driven inflow and SIDM heat transport that constitutes the central physical question of the study. We will expand the model section with a brief justification referencing the rapid velocity isotropization expected in SIDM halos on timescales shorter than the accretion window considered. We nevertheless recognize that a quantitative test with a specific-angular-momentum floor lies outside the present scope; we will explicitly list this as a limitation and defer such explorations to future multidimensional work. revision: partial

  3. Referee: [Abstract] Abstract and results description: the claim that larger initial masses, steeper profiles, and higher cross sections enhance accretion is presented as a simulation outcome, but without quantitative sensitivity analysis or resolution studies it is unclear whether the reported factor-of-100 growth is converged or sensitive to the operator-splitting treatment of conduction and self-gravity.

    Authors: The manuscript already contains a parameter survey (varying seed mass, initial density slope, and scattering cross-section) whose outcomes are shown in the figures and text. To address the concern directly we will insert a new table in the results section that tabulates accretion rates at multiple resolutions together with the corresponding operator-split time-step errors, thereby demonstrating convergence of the reported growth factor within the 1D framework. revision: yes

Circularity Check

0 steps flagged

No circularity: results obtained via direct numerical integration of fluid equations

full rationale

The paper reports outcomes from a custom 1D spherically symmetric finite-volume hydrodynamic code applied to self-gravitating fluid equations with thermal conduction. No analytic derivation chain exists whose steps could reduce to fitted parameters, self-definitions, or self-citations. The reported accretion rates (e.g., 100 M⊙ seed to 10^4 M⊙ in 2 Myr under SIS conditions) are generated outputs of the numerical integration under stated initial profiles and boundary conditions, not predictions obtained by renaming or fitting subsets of the same data. The spherical-symmetry assumption is an explicit modeling choice whose validity is external to any internal reduction; it does not create a self-referential loop within the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; the fluid model and spherical symmetry are implicit modeling choices whose validity is not quantified here.

pith-pipeline@v0.9.1-grok · 5707 in / 1158 out tokens · 27742 ms · 2026-07-03T06:46:58.971887+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

63 extracted references · 51 canonical work pages · 22 internal anchors

  1. [1]

    The MC slope limiter In computational fluid dynamics (CFD), slope lim- iters are used to reconstruct interfacial values from cell- averaged quantities [52]. Given the cell-averaged state vectorU i, the left and right reconstructed states at the interfacei+ 1/2, denoted byU L i+1/2 andU R i+1/2, can be obtained through linear reconstruction as follows: UL ...

  2. [2]

    The Roe Riemann solver In computational fluid dynamics, a Riemann problem is an initial-value problem for a hyperbolic system of conservation laws with piecewise-constant initial data. When the computational domain is discretized into cells, a local Riemann problem is defined at each interface i+ 1/2 by the reconstructed left and right states,U L i+1/2 12...

  3. [3]

    The Waning of the WIMP? A Review of Models, Searches, and Constraints

    G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mam- brini, M. Pierre, S. Profumo, and F. S. Queiroz, The waning of the WIMP? A review of models, searches, and constraints, Eur. Phys. J. C78, 203 (2018), arXiv:1703.07364 [hep-ph]

  4. [4]

    Dark Matter Search Results from 4.2 Tonne-Years of Exposure of the LUX-ZEPLIN (LZ) Experiment

    J. Aalberset al.(LZ), Dark Matter Search Results from 4.2 Tonne-Years of Exposure of the LUX-ZEPLIN (LZ) Experiment, Phys. Rev. Lett.135, 011802 (2025), arXiv:2410.17036 [hep-ex]

  5. [5]

    Aprileet al.(XENON), First Search for Light Dark Matter in the Neutrino Fog with XENONnT, Phys

    E. Aprileet al.(XENON), First Search for Light Dark Matter in the Neutrino Fog with XENONnT, Phys. Rev. Lett.134, 111802 (2025), arXiv:2409.17868 [hep-ex]

  6. [6]

    D. N. Spergel and P. J. Steinhardt, Observational ev- idence for selfinteracting cold dark matter, Phys. Rev. Lett.84, 3760 (2000), arXiv:astro-ph/9909386

  7. [7]

    Subhaloes in Self-Interacting Galactic Dark Matter Haloes

    M. Vogelsberger, J. Zavala, and A. Loeb, Subhaloes in Self-Interacting Galactic Dark Matter Haloes, Mon. Not. Roy. Astron. Soc.423, 3740 (2012), arXiv:1201.5892 [astro-ph.CO]

  8. [8]

    Dark Matter Self-interactions and Small Scale Structure

    S. Tulin and H.-B. Yu, Dark Matter Self-interactions and Small Scale Structure, Phys. Rept.730, 1 (2018), arXiv:1705.02358 [hep-ph]

  9. [9]

    Jianget al., A semi-analytic study of self-interacting dark-matter haloes with baryons, Mon

    F. Jianget al., A semi-analytic study of self-interacting dark-matter haloes with baryons, Mon. Not. Roy. Astron. Soc.521, 4630 (2023), arXiv:2206.12425 [astro-ph.CO]

  10. [10]

    D. Yang, E. O. Nadler, H.-B. Yu, and Y.-M. Zhong, A parametric model for self-interacting dark matter halos, JCAP02(2024), 032, arXiv:2305.16176 [astro-ph.CO]

  11. [11]

    Zhong, D

    Y.-M. Zhong, D. Yang, and H.-B. Yu, The impact of baryonic potentials on the gravothermal evolution of self- interacting dark matter haloes, Mon. Not. Roy. Astron. Soc.526, 758 (2023), arXiv:2306.08028 [astro-ph.CO]

  12. [12]

    Hu, R.-G

    L. Hu, R.-G. Cai, and S.-J. Wang, Distinctive GWBs from eccentric inspiraling SMBH binaries with a DM spike, JCAP02(2025), 067, arXiv:2312.14041 [gr-qc]

  13. [13]

    Z. Jia, F. Jiang, S. Li, R. Li, J. Wang, and L. Zhu, An Enhanced Isothermal Jeans Approach to Constraining Dark Matter Self-Interactions from Galactic Kinematics (2026), arXiv:2601.17118 [astro-ph.GA]

  14. [14]

    S. Li, M. S. Fischer, Z. Jia, F. Jiang, R. Li, and H.- B. Yu, Testing the isothermal Jeans model for self- interacting dark matter halos in the collapse phase (2026), arXiv:2603.01772 [astro-ph.CO]

  15. [15]

    Dark Matter Halos as Particle Colliders: A Unified Solution to Small-Scale Structure Puzzles from Dwarfs to Clusters

    M. Kaplinghat, S. Tulin, and H.-B. Yu, Dark Matter Halos as Particle Colliders: Unified Solution to Small- Scale Structure Puzzles from Dwarfs to Clusters, Phys. Rev. Lett.116, 041302 (2016), arXiv:1508.03339 [astro- ph.CO]

  16. [16]

    J. S. Bullock and M. Boylan-Kolchin, Small-Scale Chal- lenges to the ΛCDM Paradigm, Ann. Rev. Astron. As- trophys.55, 343 (2017), arXiv:1707.04256 [astro-ph.CO]

  17. [17]

    W. J. G. de Blok, The Core-Cusp Problem, Advances in 14 Astronomy2010, 789293 (2010), arXiv:0910.3538 [astro- ph.CO]

  18. [18]

    Self-Interacting Dark Matter with Mass Segregation: A Unified Explanation of Dwarf Cores and Small-Scale Lenses

    D. Yang, Y.-Z. Fan, S. Hou, and Y.-L. S. Tsai, Self- interacting dark matter with mass segregation: a unified explanation of dwarf cores and small-scale lenses, Sci. Bull.71, 1349 (2026), arXiv:2506.14898 [astro-ph.CO]

  19. [19]

    Too big to fail? The puzzling darkness of massive Milky Way subhaloes

    M. Boylan-Kolchin, J. S. Bullock, and M. Kaplinghat, Too big to fail? The puzzling darkness of massive Milky Way subhaloes, Mon. Not. Roy. Astron. Soc.415, L40 (2011), arXiv:1103.0007 [astro-ph.CO]

  20. [20]

    Constraining Self-Interacting Dark Matter with the Milky Way's dwarf spheroidals

    J. Zavala, M. Vogelsberger, and M. G. Walker, Constrain- ing Self-Interacting Dark Matter with the Milky Way’s dwarf spheroidals, Mon. Not. Roy. Astron. Soc.431, L20 (2013), arXiv:1211.6426 [astro-ph.CO]

  21. [21]

    K. A. Omanet al., The unexpected diversity of dwarf galaxy rotation curves, Mon. Not. Roy. Astron. Soc.452, 3650 (2015), arXiv:1504.01437 [astro-ph.GA]

  22. [22]

    How the Self-Interacting Dark Matter Model Explains the Diverse Galactic Rotation Curves

    A. Kamada, M. Kaplinghat, A. B. Pace, and H.-B. Yu, How the Self-Interacting Dark Matter Model Explains the Diverse Galactic Rotation Curves, Phys. Rev. Lett. 119, 111102 (2017), arXiv:1611.02716 [astro-ph.GA]

  23. [23]

    Kaplinghat, T

    M. Kaplinghat, T. Ren, and H.-B. Yu, Dark Matter Cores and Cusps in Spiral Galaxies and their Explanations, JCAP06(2020), 027, arXiv:1911.00544 [astro-ph.GA]

  24. [24]

    M. G. Roberts, M. Kaplinghat, M. Valli, and H.-B. Yu, Gravothermal collapse and the diversity of galac- tic rotation curves, Phys. Rev. D111, 103041 (2025), arXiv:2407.15005 [astro-ph.GA]

  25. [25]

    Yang and H.-B

    D. Yang and H.-B. Yu, Self-interacting dark matter and small-scale gravitational lenses in galaxy clusters, Phys. Rev. D104, 103031 (2021), arXiv:2102.02375 [astro- ph.GA]

  26. [26]

    Yang and H.-B

    D. Yang and H.-B. Yu, Gravothermal evolution of dark matter halos with differential elastic scattering, JCAP09 (2022), 077, arXiv:2205.03392 [astro-ph.CO]

  27. [27]

    Gilman, Y.-M

    D. Gilman, Y.-M. Zhong, and J. Bovy, Constrain- ing resonant dark matter self-interactions with strong gravitational lenses, Phys. Rev. D107, 103008 (2023), arXiv:2207.13111 [astro-ph.CO]

  28. [28]

    S. Hou, D. Yang, N. Li, and G. Li, A universal analytic model for gravitational lensing by self-interacting dark matter halos, JCAP08(2025), 048, arXiv:2502.14964 [astro-ph.CO]

  29. [29]

    Little Dark Dot

    S. Liet al., The “Little Dark Dot”: Evidence for Self-interacting Dark Matter in the Strong Lens SDSS J0946+1006?, Astrophys. J.994, 201 (2025), arXiv:2504.11800 [astro-ph.GA]

  30. [30]

    Lei, Y.-Y

    L. Lei, Y.-Y. Wang, Q. Li, J. Dong, Z.-F. Wang, W.-L. Lin, Y.-P. Shu, X.-Y. Cao, D.-N. Yang, and Y.-Z. Fan, A Dense Dark Matter Core of the Subhalo in the Strong Lensing System JVAS B1938+666, Astrophys. J. Lett. 991, L27 (2025), arXiv:2509.07808 [astro-ph.CO]

  31. [31]

    D. Kong, E. O. Nadler, and H.-B. Yu, Strong Lens- ing Perturbers from the SIDM Concerto Suite (2025), arXiv:2510.01491 [astro-ph.CO]

  32. [32]

    Yu, Three Birds with One Stone: Core-Collapsed SIDM Halos as the Common Origin of Dense Per- turbers in Lenses, Streams, and Satellites (2025), arXiv:2510.11006 [astro-ph.GA]

    H.-B. Yu, Three Birds with One Stone: Core-Collapsed SIDM Halos as the Common Origin of Dense Per- turbers in Lenses, Streams, and Satellites (2025), arXiv:2510.11006 [astro-ph.GA]

  33. [33]

    Matthee, R

    J. Mattheeet al., Little Red Dots: An Abundant Popu- lation of Faint Active Galactic Nuclei at z∼5 Revealed by the EIGER and FRESCO JWST Surveys, Astrophys. J.963, 129 (2024), arXiv:2306.05448 [astro-ph.GA]

  34. [34]

    Z. Li, K. Inayoshi, K. Chen, K. Ichikawa, and L. C. Ho, Little Red Dots: Rapidly Growing Black Holes Reddened by Extended Dusty Flows, Astrophys. J.980, 36 (2025)

  35. [35]

    Feng, H.-B

    W.-X. Feng, H.-B. Yu, and Y.-M. Zhong, Dark Bondi Accretion Aided by Baryons and the Origin of JWST Little Red Dots (2025), arXiv:2506.17641 [astro-ph.GA]

  36. [36]

    Inayoshi, Little Red Dots as the Very First Activity of Black Hole Growth, Astrophys

    K. Inayoshi, Little Red Dots as the Very First Activity of Black Hole Growth, Astrophys. J. Lett.988, L22 (2025)

  37. [37]

    Jiang, Z

    F. Jiang, Z. Jia, H. Zheng, L. C. Ho, K. Inayoshi, X. Shen, M. Vogelsberger, and W.-X. Feng, Formation of the Lit- tle Red Dots from the Core Collapse of Self-interacting Dark Matter Halos, Astrophys. J. Lett.996, L19 (2026), arXiv:2503.23710 [astro-ph.GA]

  38. [38]

    J. I. Read and G. Gilmore, Can supermassive black holes alter cold dark matter cusps through accretion?, Mon. Not. Roy. Astron. Soc.339, 949 (2003), arXiv:astro- ph/0210658

  39. [39]

    S. L. Shapiro, Spikes and accretion of unbound, collision- less matter around black holes, Phys. Rev. D108, 083037 (2023), arXiv:2310.13739 [astro-ph.GA]

  40. [40]

    Self-Interacting Dark Matter Halos and the Gravothermal Catastrophe

    S. Balberg, S. L. Shapiro, and S. Inagaki, Self-interacting dark matter halos and the gravothermal catastrophe, As- trophys. J.568, 475 (2002), arXiv:astro-ph/0110561

  41. [41]

    Feng, H.-B

    W.-X. Feng, H.-B. Yu, and Y.-M. Zhong, Dynamical in- stability of collapsed dark matter halos, JCAP05(2022), 036, arXiv:2108.11967 [astro-ph.CO]

  42. [42]

    Feng, H.-B

    W.-X. Feng, H.-B. Yu, and Y.-M. Zhong, Seeding Su- permassive Black Holes with Self-interacting Dark Mat- ter: A Unified Scenario with Baryons, Astrophys. J. Lett. 914, L26 (2021), arXiv:2010.15132 [astro-ph.CO]

  43. [43]

    H.-P. Gu, F. Jiang, X. Chen, and R. Li, Non-Equilibrium Relativistic Core Collapse of Self-Interacting Dark Mat- ter Halos – Limits On Seed Black Hole Mass (2026), arXiv:2601.17117 [astro-ph.CO]

  44. [44]

    J. P. Ostriker, Collisional dark matter and the origin of massive black holes, Phys. Rev. Lett.84, 5258 (2000), arXiv:astro-ph/9912548

  45. [45]

    J. F. Hennawi and J. P. Ostriker, Observational con- straints on the self interacting dark matter scenario and the growth of supermassive black holes, Astrophys. J. 572, 41 (2002), arXiv:astro-ph/0108203

  46. [46]

    J. Hu, Y. Shen, Y.-Q. Lou, and S. Zhang, Forming su- permassive black holes by accreting dark and baryon matter, Mon. Not. Roy. Astron. Soc.365, 345 (2006), arXiv:astro-ph/0510222

  47. [47]

    V. M. Sabarish, M. Br¨ uggen, K. Schmidt-Hoberg, and M. S. Fischer, Accretion of self-interacting dark matter onto supermassive black holes, Astron. Astrophys.703, A142 (2025), arXiv:2505.14779 [astro-ph.CO]

  48. [48]

    S. L. Shapiro and V. Paschalidis, Self-interacting dark matter cusps around massive black holes, Phys. Rev. D 89, 023506 (2014), arXiv:1402.0005 [astro-ph.CO]

  49. [49]

    J. F. Navarro, C. S. Frenk, and S. D. M. White, The Structure of cold dark matter halos, Astrophys. J.462, 563 (1996), arXiv:astro-ph/9508025

  50. [50]

    Formation and Evolution of the Self-Interacting Dark Matter Halos

    K.-J. Ahn and P. R. Shapiro, Formation and evolution of the self-interacting dark matter halos, Mon. Not. Roy. Astron. Soc.363, 1092 (2005), arXiv:astro-ph/0412169

  51. [51]

    Gravothermal collapse of isolated self-interacting dark matter haloes: N-body simulation versus the fluid model

    J. Koda and P. R. Shapiro, Gravothermal collapse of iso- lated self-interacting dark matter haloes: N-body simula- tion versus the fluid model, Mon. Not. Roy. Astron. Soc. 415, 1125 (2011), arXiv:1101.3097 [astro-ph.CO]

  52. [52]

    Nishikawa, K

    H. Nishikawa, K. K. Boddy, and M. Kaplinghat, Accel- erated core collapse in tidally stripped self-interacting dark matter halos, Phys. Rev. D101, 063009 (2020), 15 arXiv:1901.00499 [astro-ph.GA]

  53. [53]

    Blazek,Computational fluid dynamics: principles and applications(Butterworth-Heinemann, 2015)

    J. Blazek,Computational fluid dynamics: principles and applications(Butterworth-Heinemann, 2015)

  54. [54]

    R. J. LeVeque,Finite volume methods for hyperbolic problems, Vol. 31 (Cambridge university press, 2002)

  55. [55]

    P. L. Roe, Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes, Journal of Computa- tional Physics43, 357 (1981)

  56. [56]

    E. F. Toro,Riemann solvers and numerical methods for fluid dynamics: a practical introduction(Springer Sci- ence & Business Media, 2013)

  57. [57]

    van Leer, Towards the Ultimate Conservative Differ- ence Scheme

    B. van Leer, Towards the Ultimate Conservative Differ- ence Scheme. IV. A New Approach to Numerical Convec- tion, Journal of Computational Physics23, 276 (1977)

  58. [58]

    N. J. Outmezguine, K. K. Boddy, S. Gad-Nasr, M. Kaplinghat, and L. Sagunski, Universal gravother- mal evolution of isolated self-interacting dark matter ha- los for velocity-dependent cross-sections, Mon. Not. Roy. Astron. Soc.523, 4786 (2023), arXiv:2204.06568 [astro- ph.GA]

  59. [59]

    Gad-Nasr, K

    S. Gad-Nasr, K. K. Boddy, M. Kaplinghat, N. J. Out- mezguine, and L. Sagunski, On the late-time evolution of velocity-dependent self-interacting dark matter halos, JCAP05(2024), 131, arXiv:2312.09296 [astro-ph.GA]

  60. [60]

    D. Yang, E. O. Nadler, and H.-B. Yu, Strong Dark Mat- ter Self-interactions Diversify Halo Populations within and surrounding the Milky Way, Astrophys. J.949, 67 (2023), arXiv:2211.13768 [astro-ph.GA]

  61. [61]

    Godunov, Finite difference methods for the computa- tion of discontinuous solutions of the equations of fluid dynamics, Mat

    S. Godunov, Finite difference methods for the computa- tion of discontinuous solutions of the equations of fluid dynamics, Mat. Sb.47, 271 (1959)

  62. [62]

    Harten, P

    A. Harten, P. D. Lax, and B. van Leer, On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws, SIAM Rev.25, 35 (1983)

  63. [63]

    E. F. Toro, M. Spruce, and W. Speares, Restoration of the contact surface in the HLL-Riemann solver, Shock Waves4, 25 (1994)