Dynamical evolutions in globular clusters and dwarf galaxies: Conduction fluid simulations
Phys. Rev. D 112, 023040 – Published 28 July, 2025
DOI: https://doi.org/10.1103/z1fy-2p8s
Abstract
We present a new two-fluid conduction scheme to simulate the evolution of an isolated, self-gravitating, equilibrium cluster of stars and collisionless dark matter on secular (gravothermal) timescales. We integrate the equations in Lagrangian coordinates via a second-order, semi-implicit algorithm, which is unconditionally stable when the mass of the lighter species is much less than that of the heavier species. The method can be straightforwardly generalized to handle a multispecies system with a population of stars or components beyond collisionless dark matter and stars. We apply the method to simulate the dynamical evolution of stellar-dark matter systems, exploring the consequences of mass segregation and gravothermal core collapse, and assessing those effects for observed globular clusters and dwarf galaxies in the Local Volume.
