RNAAS
Einasto gravitational potentials have difficulty to hold spherically
symmetric stellar systems with cores
Abstract
It was known that an ideal spherically symmetric stellar system with isotropic velocities and an inner core cannot reside in a Navarro, Frenk, and White (NFW) gravitational potential. The incompatibility can be pinned down to the radial gradient of the NFW potential in the very center of the system, which differs from zero. The gradient is identically zero in an Einasto potential, also an alternative representation of the dark matter (DM) halos created by the kind of cold DM (CDM) defining the current cosmological model. Here we show that, despite the inner gradient being zero, stellar cores are also inconsistent with Einasto potentials. This result may have implications to constrain the nature of DM through interpreting the stellar cores often observed in dwarf galaxies.
1 Discusion
Sánchez Almeida et al. (2023) showed that an ideal spherical stellar distribution having an inner core, i.e., having
(1) |
is not consistent with a NFW potential (named after Navarro et al., 1997). The symbol stands for the stellar mass volume density, with the distance to the center of the gravitating system. This result would be only a curiosity if it were not for the fact that stellar cores are very common in dwarf galaxies (e.g., Moskowitz & Walker, 2020; Richstein et al., 2024; Montes et al., 2024). Even if it is an academic result for ideal systems, the existence of cores seems to suggest the gravitational potential to differ from the NFW potential characteristic of the CDM particles defining the current cosmological model. Deviations from the NFW potential are well known from kinematical measurements (e.g., Oh et al., 2015) and they are ascribed to the effect of stellar feedback processes modifying the overall gravitational potential (e.g., Governato et al., 2010). This baryon driven feedback is expected to be inoperative for galaxies with stellar mass (e.g., Chan et al., 2015; Jackson et al., 2021), therefore, the existence of stellar cores in tiny galaxies would reflect the need to go beyond the CDM paradigm.
Sánchez Almeida et al. (2023) proved the inconsistency using the Eddington Inversion Method (EIM; e.g., Binney & Tremaine, 2008) which provides the distribution function in the phase space needed to reproduce a stellar distribution in an assumed gravitational potential. For a cored stellar profile to reside in a NFW potential has to be negative, which is unphysical. The proof employs an intermediate step in the EIM formalism, which reads,
(2) |
with the relative gravitational potential and the total energy per unit mass of a star. In the case of a NFW potential,
(3) |
which together with Eq. (1) implies the integral in the right-hand-side of Eq. (2) to be zero when . For the integral to be zero must be negative somewhere in the phase space, which is unphysical and evidences the incompatibility.
The CDM halos giving rise to a NFW potential are equally well represented by Einasto profiles (e.g., Wang et al., 2020), where the derivative is zero (e.g., Retana-Montenegro et al., 2012)
(4) |
and, consequently, for which the above argument of incompatibility does not hold. The question arises as to whether the incompatibility of stellar cores is specific to the NFW potential or if it more generally affects all CDM-produced halos. Recall that the NFW potential stems from a mass density profile scaling as , therefore, with an artificial singularity when . The singularity is not present in the Einasto profile which scales as with and two positive constants. This Note is meant to show that the inconsistency still remains in the case of an Einasto potential. In terms of Eq. (2), the Einasto profile has
(5) |
for values of where
(6) |
and this condition is enough to force the need of .
The argument is illustrated with the simulation in Fig. 1, which shows a cored stellar profile (the black line; a Plummer profile) with size and central density chosen to mimic a dwarf galaxy. The red line shows the NFW profile used to compute the potential giving rise to the distribution function shown in Fig. 1b (the red line). It was evaluated with the numerical tools developed by Sánchez Almeida et al. (2023). As expected, there is a range of values for which . The same exercise is carried out with an Einasto profile having the same characteristic radius and density as the NFW profile, as in CDM simulations (e.g., Wang et al., 2020). Firstly, we note that the corresponding also has negative parts (the blue dashed line in Fig. 1b), making also the Einasto potential inconsistent with the stellar core. Secondly, and in contrast with the NFW potential, the distribution function for the Einasto potential becomes positive near the largest energy. Such energies are reached only by the stars moving close to the center of the potential at all times. Thus, the incompatibility is not at but at a finite radius within the core (Eq. [6]) but where the potential keeps changing (Eq. [5]).
2 Conclusions
The original inconsistency happens for spherically symmetric systems having isotropic velocities, and representing the CDM with NFW profiles. However, none of these simplifying assumptions seem to be responsible for the incompatibility. Sánchez Almeida et al. (2023) showed that the incompatibility remains for other potentials with an inner slope much shallower than the value of characteristic of the NFW profiles. They also showed that the assumption of isotropic velocities can be relaxed since the incompatibility remains for radial orbits and for isotropic orbits that turn into radial in the outskirts, as expected from cosmological numerical simulations (e.g., El-Badry et al., 2017; Orkney et al., 2023). On the other hand, Sánchez Almeida et al. (2024) proved that the incompatibility persists for axi-symmetric systems, which implies that it goes beyond the assumption of spherical symmetry. This Note adds that the incompatibility remains for Einasto potentials, which are also a good representation of the CDM halos and, in contrast to the NFW profiles, do not diverge at the center of the system. Thus, the incompatibility is not artificially set by the mathematical singularity at present in the NFW profiles.
References
- Binney & Tremaine (2008) Binney, J., & Tremaine, S. 2008, Galactic Dynamics: Second Edition
- Chan et al. (2015) Chan, T. K., Kereš, D., Oñorbe, J., et al. 2015, MNRAS, 454, 2981
- El-Badry et al. (2017) El-Badry, K., Wetzel, A. R., Geha, M., et al. 2017, ApJ, 835, 193
- Governato et al. (2010) Governato, F., Brook, C., Mayer, L., et al. 2010, Nature, 463, 203
- Jackson et al. (2021) Jackson, R. A., Martin, G., Kaviraj, S., et al. 2021, MNRAS, 502, 4262
- Montes et al. (2024) Montes, M., Trujillo, I., Karunakaran, A., et al. 2024, A&A, 681, A15
- Moskowitz & Walker (2020) Moskowitz, A. G., & Walker, M. G. 2020, ApJ, 892, 27
- Navarro et al. (1997) Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493
- Oh et al. (2015) Oh, S.-H., Hunter, D. A., Brinks, E., et al. 2015, AJ, 149, 180
- Orkney et al. (2023) Orkney, M. D. A., Taylor, E., Read, J. I., et al. 2023, MNRAS, 525, 3516
- Retana-Montenegro et al. (2012) Retana-Montenegro, E., van Hese, E., Gentile, G., Baes, M., & Frutos-Alfaro, F. 2012, A&A, 540, A70
- Richstein et al. (2024) Richstein, H., Kallivayalil, N., Simon, J. D., et al. 2024, ApJ, 967, 72
- Sánchez Almeida et al. (2023) Sánchez Almeida, J., Plastino, A. R., & Trujillo, I. 2023, ApJ, 954, 153
- Sánchez Almeida et al. (2024) —. 2024, A&A, submitted, …
- Virtanen et al. (2020) Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261
- Wang et al. (2020) Wang, J., Bose, S., Frenk, C. S., et al. 2020, Nature, 585, 39