No Arabic abstract
A Fuzzy Dark Matter (FDM) halo consists of a soliton core close to the center and an NFW-like density profile in the outer region. Previous investigations found that the soliton core exhibits temporal oscillations and random walk excursions around the halo center. Analyzing a set of numerical simulations, we show that both phenomena can be understood as the results of wave interference -- a suitable superposition of the ground (solitonic) state and excited states in a fixed potential suffices to account for the main features of these phenomena. Such an eigenmode analysis can shed light on the evolution of a satellite halo undergoing tidal disruption. As the outer halo is stripped away, reducing the amplitudes of the excited states, the ground state evolves adiabatically. This suggests diminished soliton oscillations and random walk excursions, an effect to consider in deducing constraints from stellar heating.
Fuzzy dark matter (FDM) has been a promising alternative to standard cold dark matter. The model consists of ultralight bosons with mass $m_b sim 10^{-22}$ eV and features a quantum-pressure-supported solitonic core that oscillates. In this work, we show that the soliton density oscillations persist even after significant tidal stripping of the outer halo. We report two intrinsic yet distinct timescales associated, respectively, with the ground-state soliton wavefunction $tau_{00}$ and the soliton density oscillations $tau_text{soliton}$, obeying $tau_text{soliton} /tau_{00} simeq 2.3$. The central star cluster (SC) in Eridanus II has a characteristic timescale $tau_text{soliton} / tau_text{SC} sim 2$ to $3$ that deviates substantially from unity. As a result, we demonstrate, both analytically and numerically with three-dimensional self-consistent FDM simulations, that the gravitational heating of the SC owing to soliton density oscillations is negligible irrespective of $m_b$. We also show that the subhalo mass function to form Eridanus II does not place a strong constraint on $m_b$. These results are contrary to the previous findings by Marsh & Niemeyer (2019).
For idealized (spherical, smooth) dark matter halos described by single-parameter density profiles (such as the NFW profile) there exists a one-to-one mapping between the energy of the halo and the scale radius of its density profile. The energy therefore uniquely determines the concentration parameter of such halos. We exploit this fact to predict the concentrations of dark matter halos via a random walk in halo energy space. Given a full merger tree for a halo, the total internal energy of each halo in that tree is determined by summing the internal and orbital energies of progenitor halos. We show that, when calibrated, this model can accurately reproduce the mean of the concentration--mass relation measured in N-body simulations, and reproduces more of the scatter in that relation than previous models. We further test this model by examining both the autocorrelation of scale radii across time, and the correlations between halo concentration and spin, and comparing to results measured from cosmological N-body simulations. In both cases we find that our model closely matches the N-body results. Our model is implemented within the open source Galacticus toolkit.
Simulations of ultralight, $sim 10^{-22},rm eV$, bosonic dark matter exhibit rich wave-like structure, including a soliton core within a surrounding halo that continuously self-interferes on the de Broglie scale. We show here that as an inherent consequence, the soliton undergoes a confined random walk at the base of the halo potential. This is significant for the fate of the ancient central star cluster in Eridanus II, as the agitated soliton gravitationally shakes the star cluster in and out of the soliton on a time scale of $sim 100,rm Myr$, so complete tidal disruption of the star cluster can occur within $sim 1,rm Gyr$. This destructive effect can be mitigated by tidal stripping of the halo of Eridanus II, thereby reducing the agitation, depending on its orbit around the Milky Way. Our simulations show the Milky Way tide affects the halo much more than the soliton, so the star cluster in Eridanus II can survive for over $5,rm Gyr$ within the soliton if it formed after significant halo stripping.
We extend the random-walk model of Vitvitska et al. for predicting the spins of dark matter halos from their merger histories. Using updated merger rates, orbital parameter distributions, and N-body constraints we show that this model can accurately reproduce the distribution of spin parameters measured in N-body simulations when we include a weak correlation between the spins of halos and the angular momenta of infalling subhalos. We further show that this model is in approximate agreement with the correlation of the spin magnitude over time as determined from N-body simulations, while it slightly underpredicts the correlation in the direction of the spin vector measured from the same simulations. This model is useful for predicting spins from merger histories derived from non-N-body sources, thereby circumventing the need for very high resolution simulations to permit accurate measurements of spins. It may be particularly relevant to modeling systems which accumulate angular momentum from halos over time (such as galactic discs) - we show that this model makes small but significant changes in the distribution of galactic disc sizes computed using the Galacticus semi-analytic galaxy formation model.
We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $msimmathcal{O}(10^{-22}),$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that emph{(1)}~for FDM masses $mgtrsim 10^{-21}, {rm eV}, c^{-2}$, the timing problem of the Fornax dwarf spheroidals globular clusters is no longer solved and emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{rm rel}$ should require detailed numerical simulations for $left(M/10^9~M_{odot}right) left(m/10^{-22}~{rm eV}right)left(100~{rm km}~{rm s}^{-1}/v_{rm rel}right) sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.