Search

We explore the properties of Milky Way subhalos in self-interacting dark matter models for moderate cross sections of 1 to 5 cm$^2$g$^{-1}$ using high-resolution zoom-in N-body simulations. We include the gravitational potential of a baryonic disk and bulge matched to the Milky Way, which is critical for getting accurate predictions. The predicted number and distribution of subhalos within the host halo are similar for 1 and 5 cm$^2$g$^{-1}$ models, and they agree with observations of Milky Way satellite galaxies only if subhalos with peak circular velocity over all time > 4.5 km/s are able to form galaxies. We do not find distinctive signatures in the pericenter distribution of the subhalos that could help distinguish the models. Using an analytic model to extend the simulation results, we are able to show that subhalos in models with cross sections between 1 and 5 cm$^2$g$^{-1}$ are not dense enough to match the densest ultra-faint and classical dwarf spheroidal galaxies in the Milky Way. This motivates velocity-dependent cross sections with values larger than 5 cm$^2$g$^{-1}$ at the velocities relevant for the satellites such that core collapse would occur in some of the ultra-faint and classical dwarf spheroidals., Version accepted for publication by MNRAS. 19 pages, 16 figures, 2 tables, 5 appendices

We investigate cosmological structure formation in Fuzzy Dark Matter (FDM) with an attractive self-interaction (SI) with numerical simulations. Such a SI would arise if the FDM boson were an ultra-light axion, which has a strong CP symmetry-breaking scale (decay constant). Although weak, the attractive SI may be strong enough to counteract the quantum 'pressure' and alter structure formation. We find in our simulations that the SI can enhance small-scale structure formation, and soliton cores above a critical mass undergo a phase transition, transforming from dilute to dense solitons., 10 pages, 3 figures, submitted to mnras

Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes:gizmo and arepo. We begin with identical initial conditions for an isolated 1010 M⊙ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in gizmo and arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2\, g^{-1}$. Our tests are restricted to the core expansion phase, where the core density decreases and core radius increases with time. We find better than 30 per cent agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10 per cent outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2\, g^{-1}$, but finer distinctions would require further investigation.

In recent years, the scalar field dark matter (SFDM), also called ultralight bosonic dark matter, has received considerable attention due to the number of problems it might help to solve. Among these are the cusp-core problem and the abundance of small structures of the standard cold dark matter (CDM) model. In this paper we show that multi-state solutions of the low energy and weak gravitational field limit of field equations, interpreted as galactic halo density profiles, can provide a possible explanation to the anisotropic distribution of satellite galaxies observed in the Milky Way, M31 and Centaurus A, where satellites trajectories seem to concentrate on planes close to the poles of the galaxies instead of following homogeneously distributed trajectories. The core hypothesis is that multi-state solutions of the equations describing the dynamics of this dark matter candidate, namely, the Gross-Pitaevskii-Poisson equations, with monopolar and dipolar contributions, can possibly explain the anisotropy of satellite trajectories. In order to construct a proof of concept, we study the trajectories of a number of test particles traveling on top of the gravitational potential due to a multi-state halo with modes (1,0,0)+(2,1,0). The result is that particles accumulate asymptotically in time on planes passing close to the poles. Satellite galaxies are not test particles but interpreted as such, our results indicate that in the asymptotic time their trajectories do not distribute isotropically, instead they prefer to have orbital poles accumulating near the equatorial plane of the multistate halo. The concentration of orbital poles depends on whether the potential is monopolar or dipolar dominated., This version has been accepted for publication in Physical Review D

We compare a suite of four simulated dwarf galaxies formed in 10$^{10} M_{\odot}$ haloes of collisionless Cold Dark Matter (CDM) with galaxies simulated in the same haloes with an identical galaxy formation model but a non-zero cross-section for dark matter self-interactions. These cosmological zoom-in simulations are part of the Feedback In Realistic Environments (FIRE) project and utilize the FIRE-2 model for hydrodynamics and galaxy formation physics. We find the stellar masses of the galaxies formed in Self-Interacting Dark Matter (SIDM) with $\sigma/m= 1\, cm^2/g$ are very similar to those in CDM (spanning $M_{\star} \approx 10^{5.7 - 7.0} M_{\odot}$) and all runs lie on a similar stellar mass -- size relation. The logarithmic dark matter density slope ($\alpha=d\log \rho / d\log r$) in the central $250-500$ pc remains steeper than $\alpha= -0.8$ for the CDM-Hydro simulations with stellar mass $M_{\star} \sim 10^{6.6} M_{\odot}$ and core-like in the most massive galaxy. In contrast, every SIDM hydrodynamic simulation yields a flatter profile, with $\alpha >-0.4$. Moreover, the central density profiles predicted in SIDM runs without baryons are similar to the SIDM runs that include FIRE-2 baryonic physics. Thus, SIDM appears to be much more robust to the inclusion of (potentially uncertain) baryonic physics than CDM on this mass scale, suggesting SIDM will be easier to falsify than CDM using low-mass galaxies. Our FIRE simulations predict that galaxies less massive than $M_{\star} < 3 \times 10^6 M_{\odot}$ provide potentially ideal targets for discriminating models, with SIDM producing substantial cores in such tiny galaxies and CDM producing cusps., Comment: 10 Pages, 7 figures, submitted to MNRAS

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec (kpc) scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter we investigate how the first galaxies are assembled in a `fuzzy' dark matter (FDM) cosmology where dark matter is an ultralight $\sim 10^{-22}$~eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical soliton-like cores which are unstable under gravity and collapse into kpc-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter -- a smoking gun signature of FDM., Comment: 6 pages, 2 figures, accepted to Phys. Rev. Lett

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like ("WDM") model often used as a proxy for BECDM. Our simulations confirm that "WDM" is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/"WDM" compared to the CDM case: in BECDM/"WDM" first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in "WDM", is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM., Comment: 19 pages, 12 figures, submitted to MNRAS

Building upon results of cosmological simulations of ultra-light scalar field dark matter (SFDM), we present a comprehensive model for the density profiles of SFDM haloes as a function of halo virial mass $M_{\rm h}$ and scalar field mass $m$. The central regions of SFDM haloes are dominated by solitons with characteristic densities that increase with increasing halo mass and asymptote to CDM-like profiles at large radii. For scalar field masses $m \sim 10^{-22}$ eV, consistent with large-scale structure observations, $M_{\rm h} \sim 10^{10} \,M_\odot$ haloes have lower core densities than their Cold Dark Matter (CDM) counterparts and this alleviates the Too Big to Fail problem (TBTF) in a regime where feedback is less effective. However, higher-mass SFDM haloes with $M_{\rm h} \sim 10^{11} \,M_\odot$ are denser than their CDM counterparts at small, observationally relevant radii. We use rotation curves of $V \sim 100$ km s$^{-1}$ galaxies from the SPARC database to show that SFDM exacerbates the cusp/core and central density problems seen in CDM at this scale. We conclude that if the conventional cosmological SFDM scaling relations are correct, then baryonic feedback is required to lower densities in SFDM haloes even more so than in CDM. This motivates cosmological and self-consistent hydrodynamic simulations of SFDM to determine whether central soliton structure can be altered by realistic feedback implementations., Comment: 13 pages, 8 figures, submitted to MNRAS

The quenching `maintenance' and `cooling flow' problems are important from the Milky Way through massive cluster elliptical galaxies. Previous work has shown that some source of energy beyond that from stars and pure magnetohydrodynamic processes is required, perhaps from AGN, but even the qualitative form of this energetic input remains uncertain. Different scenarios include thermal `heating,' direct wind or momentum injection, cosmic ray heating or pressure support, or turbulent `stirring' of the intra-cluster medium (ICM). We investigate these in $10^{12}-10^{14}\,{\rm M}_{\odot}$ halos using high-resolution non-cosmological simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, including simplified toy energy-injection models, where we arbitrarily vary the strength, injection scale, and physical form of the energy. We explore which scenarios can quench without violating observational constraints on energetics or ICM gas. We show that turbulent stirring in the central $\sim100\,$kpc, or cosmic-ray injection, can both maintain a stable low-SFR halo for $>$Gyr timescales with modest energy input, by providing a non-thermal pressure which stably lowers the core density and cooling rates. In both cases, associated thermal-heating processes are negligible. Turbulent stirring preserves cool-core features while mixing condensed core gas into the hotter halo and is by far the most energy efficient model. Pure thermal heating or nuclear isotropic momentum injection require vastly larger energy, are less efficient in lower-mass halos, easily over-heat cores, and require fine-tuning to avoid driving unphysical temperature gradients or gas expulsion from the halo center., 22 pages, 16 figures

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schr\"odinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic $r^{-3}$ NFW-like profile. We find a fundamental relation of the core=to-halo mass with the dimensionless invariant $\Xi \equiv \lvert E \rvert/M^3/(Gm/\hbar)^2$ or $M_{\rm c}/M \simeq 2.6 \Xi^{1/3}$, linking the soliton to global halo properties. For $r \geq 3.5 \,r_{\rm c}$ core radii, we find equipartition between potential, classical kinetic, and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behavior driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a $k^{-1.1}$ power-law. This suggests the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a $k^{-5/3}$ Kolmogorov power-law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying the soliton-sized granules carry most of the turbulent energy in BECDM haloes., Comment: 14 pages, 7 figures, MNRAS accepted; a movie rendering of Figure 1 can be found here: https://youtu.be/UBf9ghy2lPI

Recent high-quality observations of dwarf and low surface brightness (LSB) galaxies have shown that their dark matter (DM) halos prefer flat central density profiles. On the other hand the standard cold dark matter model simulations predict a more cuspy behavior. Feedback from star formation has been widely used to reconcile simulations with observations, this might be successful in field dwarf galaxies but its success in low mass galaxies remains uncertain. One model that have received much attention is the scalar field dark matter model. Here the dark matter is a self-interacting ultra light scalar field that forms a cosmological Bose-Einstein condensate, a mass of $10^{-22}$eV/c$^2$ is consistent with flat density profiles in the centers of dwarf spheroidal galaxies, reduces the abundance of small halos, might account for the rotation curves even to large radii in spiral galaxies and has an early galaxy formation. The next generation of telescopes will provide better constraints to the model that will help to distinguish this particular alternative to the standard model of cosmology shedding light into the nature of the mysterious dark matter., Comment: 6 pages, to appear in: Proceedings of the Fourteenth Marcel Grossman Meeting on General Relativity

Using the scalar field dark matter (SFDM) model, it is proposed that galaxies form by con- densation of a scalar field (SF) very early in the Universe, forming Bose-Einstein condensate (BEC) drops (i.e. in this model, the haloes of galaxies are gigantic drops of SF). Here, as in thecold dark matter (LCDM) model, large structures form by hierarchy, and thus all the predictions of the LCDM model at large scales are reproduced by the SFDM model. This model predicts that all galaxies must be very similar and must exist for larger redshifts than in the LCDM model. In this paper, we show that BEC dark matter haloes fit the high-resolution rotation curves of a sample of 13 low-surface-brightness galaxies. We compare our fits to those obtained using Navarro-Frenk-White and pseudo-isothermal (PI) profiles. We have found bet- ter agreement with the SFDM and PI profiles. The mean value of the logarithmic inner density slopes is α =− 0.27 ± 0.18. As a second result, we find a natural way to define the core radius with the advantage of being model-independent. Using this new definition in the BEC density profile, we find that the recent observation of the constant dark matter central surface density can be reproduced. We conclude that, in light of the difficulties that the standard model is currently facing, the SFDM model could be a worthy alternative to enable us to continue exploring further.

We review the work done so far aimed at modeling in an alternative way the dark matter in the Universe: the scalar field/ Bose-Einstein condensate dark matter (SFDM/BEC) model. We discuss a number of important achievements and characteristics of the model. We also describe some of our most recent results and predictions of the model compared to those of the standard model of $\Lambda$CDM.

We investigate the gravitational constraints imposed to dark matter halos in the context of finite temperature scalar field dark matter. We find constraints to produce multiple images by dark matter only, we show that there are differences with respect to the full Bose Einstein condensate halo when the temperature of the scalar field in dark matter halos is taken into account. The non zero temperature allows the scalar field to be in excited states and recently, their existence has proved to be necessary to fit rotation curves of dark matter dominated galaxies of all sizes, it also explained the non universality of the halo density profiles. Therefore, we expect that combining observations of rotation curves and strong lensing systems can give us a clue to the nature of dark matter. Finally, we propose a method to identify the excited state of a strong lens halo, knowing various halo excited states can provide information of the scalar field dark matter halo evolution which can be tested using numerical simulations., Comment: 4 pages, 2 figures

Recent high-quality observations of low surface brightness (LSB) galaxies have shown that their dark matter (DM) halos prefer flat central density profiles. On the other hand the standard cold dark matter model simulations predict a more cuspy behavior. Feedback from star formation has been widely used to reconcile simulations with observations, this might be successful in field dwarf galaxies but its success in high mass LSB galaxies remains unclear. Additionally, including too much feedback in the simulations is a double-edged sword, in order to obtain a cored DM distribution from an initially cuspy one, feedback recipes require to remove a large quantity of baryons from the center of galaxies, however, other feedback recipes produce twice more satellite galaxies of a given luminosity and with much smaller mass to light ratios from those that are observed. Therefore, one DM profile that produces cores naturally and that does not require large amounts of feedback would be preferable. We find both requirements to be satisfied in the scalar field dark matter model. Here, we consider that the dark matter is an auto-interacting real scalar field in a thermal bath of temperature T with an initial Z 2 symmetric potential, as the universe expands the temperature drops so that the Z 2 symmetry is spontaneously broken and the field rolls down to a new minimum. We give an exact analytic solution to the Newtonian limit of this system and show both, that it satisfies the two desired requirements and that the rotation curve profile is not longer universal.

The cold dark matter (CDM) model has two unsolved issues: simulations overpredict the satellite abundance around the Milky Way (MW) and it disagrees with observations of the central densities of dwarf galaxies which prefer constant density (core) profiles.One alternative explanation known as the scalar field dark matter (SFDM) model, assumes that the dark matter is a scalar field of mass($\sim 10^{-22}$ eV/$c^2$); this model can reduce the overabundance issue due to the lack of halo formation below a mass scale of $\sim 10^8$M$_{\odot}$ and successfully fits the density distribution in dwarfs. One of the attractive features of the model is predicting core profiles in halos, although the determination of the core sizes is set by fitting the observational data. We perform \textit{N}-body simulations to explore the influence of tidal forces over a stellar distribution embedded in a SFDM halo orbiting a MW-like SFDM host halo with a disk. Our simulations intend to test the viability of SFDM as an alternative model by comparing the tidal effects that result in this paradigm with those obtained in CDM for similar mass halos. We found that galaxies in subhalos with core profiles and high central densities survive for 10 Gyr. The same occurs for galaxies in low density subhalos located far from the host disk influence, whereas satellites in low density DM halos and in tight orbits can eventually be stripped of stars. We conclude that SFDM shows consistency with results from CDM for dwarf galaxies, but naturally offer a possibility to solve the missing satellite problem., Comment: 17 pages, 6 figures, matches the accepted version in ApJ

Using high-resolution simulations with explicit treatment of stellar feedback physics based on the FIRE (Feedback In Realistic Environments) project, we study how galaxy formation and the interstellar medium (ISM) are affected by magnetic fields, anisotropic Spitzer–Braginskii conduction and viscosity, and sub-grid metal diffusion from unresolved turbulence. We consider controlled simulations of isolated (non-cosmological) galaxies but also a limited set of cosmological ‘zoom-in’ simulations. Although simulations have shown significant effects from these physics with weak or absent stellar feedback, the effects are much weaker than those of stellar feedback when the latter is modelled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs). In contrast, removing stellar feedback leads to SFRs being overpredicted by factors of ∼10–100. Without feedback, neither galactic winds nor volume-filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scaleheights with unphysically dense clumps congregating at the galactic centre. With stellar feedback, a multi-phase, turbulent medium with galactic fountains and winds is established. At currently achievable resolutions and for the investigated halo mass range 10^(10)–10^(13) M⊙, the additional physics investigated here (magnetohydrodynamic, conduction, viscosity, metal diffusion) have only weak (∼10 per cent-level) effects on regulating SFR and altering the balance of phases, outflows or the energy in ISM turbulence, consistent with simple equipartition arguments. We conclude that galactic star formation and the ISM are primarily governed by a combination of turbulence, gravitational instabilities and feedback. We add the caveat that active galactic nucleus feedback is not included in the present work.