W-H Liu,Y-N Fu,Z-G Deng,J-H. Huang
Abstract:In this paper a prescription for generating an equilibrium spherical system depicted with cuspy density profiles is extended to a core density profile, the Burkert profile, which is observationally more suitable to dwarfs. By using a time-saving Monte Carlo method instead of N-body simulations, we show that the Burkert halo initialized with the distribution function that depends only on energy is in a practically stable equilibrium. The one generated using the local Maxwellian approximation is unstable, where the flat core density structure tends to steepen. This is fundamentally different from the previous study on a halo with cuspy density profile. The deviation of the unstable "Burkert" from the initial Burkert profile is found to be closely related to taking a Gaussian as the velocity distribution at any given point. The significance of not using the local Maxwellian approximation is further demonstrated by exploring the dynamical evolution of compact super star clusters in the Burkert halo. In particular, the local Maxwellian approximation will result in underestimating the dynamical friction and overestimating sinking time scale. Accordingly, it will lead to lower probability of forming massive bulges and young nuclear star clusters in bulgeless galaxies.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is related to the construction of initial conditions for dark matter halos and their stability research. Specifically, the authors focus on how to generate a stable spherically symmetric dark matter halo with a Burkert density profile (i.e., a cored density profile), and explore the effects of using different initialization methods (local Maxwellian approximation and steady - state distribution function) on the stability of dark matter halos.
### Research Background and Problem
In cosmology and galactic dynamics, constructing initial conditions for systems with specific density profiles (such as dark matter halos) is an important but difficult task. For dark matter halos with sharp central density profiles (cuspy profiles), previous studies have shown that halos initialized using the local Maxwellian approximation may lead to unstable evolution, which in turn affects the results of subsequent dynamical simulations. However, for dark matter halos with flat - core density profiles (such as Burkert profiles), there have been fewer previous studies.
### Specific Problems
1. **Generating Stable Core Dark Matter Halos**: How to generate a stable dark matter halo with a Burkert density profile?
2. **Effects of Different Initialization Methods**: Compare the stability differences of Burkert halos initialized using the local Maxwellian approximation and the steady - state distribution function (DF).
3. **Effects of Dynamical Evolution**: Explore the effects of different initialization methods on the dynamical friction and sinking time - scale of compact super - star clusters (SSC) in dark matter halos.
### Solutions
- **Mathematical Model**: By extending the research of Kazantzidis, Magorrian & Moore (2004), the authors propose a method for generating dark matter halos with Burkert density profiles and use the Monte Carlo method for numerical simulations.
- **Formula Representation**:
- The Burkert density profile formula is:
\[
\rho=\frac{\rho_s}{(1 + r/r_s)[1+(r/r_s)^2]}
\]
where \(\rho_s\) is the central density and \(r_s\) is the scale radius.
- At the outer boundary, an exponential truncation form is adopted:
\[
\rho=\frac{\rho_s}{(1 + c)(1 + c^2)}\left(\frac{r}{r_{\text{vir}}}\right)^\alpha\exp\left[-\frac{r - r_{\text{vir}}}{r_{\text{decay}}}\right]\quad\text{for}\quad(r > r_{\text{vir}})
\]
where \(c = r_{\text{vir}}/r_s\) is the concentration parameter, and \(\alpha\) is determined by the continuity condition.
- **Initialization Methods**:
- Initialize using the exact distribution function.
- Initialize using the local Maxwellian approximation.
### Results and Discussion
- **Stability Test**: Through Monte Carlo simulations, it is found that Burkert halos initialized using the exact distribution function remain stable, while halos initialized using the local Maxwellian approximation exhibit unstable behavior, especially with a steeper density in the core region and a decreasing outer - layer density.
- **Dynamical Friction**: For compact super - star clusters (SSC), in halos initialized using the exact distribution function, SSC experiences stronger dynamical friction and has a faster sinking speed; while in halos initialized using the local Maxwellian approximation, SSC has a slower sinking speed.
### Conclusions
- **Importance**: Initializing dark matter halos using the steady - state distribution function is crucial for obtaining accurate dynamical evolution results, especially for astrophysical problems such as studying the sinking history of compact super - star clusters, the formation of nuclear star clusters, and the evolution of massive black hole seeds.
- **Future Work**: It is recommended to give priority to using the steady - state in high - resolution simulations.