Search

We present the first study on the gravitational impact of supernova feedback in an isolated soliton and a spherically symmetric dwarf SFDM halo of virial mass $1\times 10^{10}\mathrm{M_\odot}$. We use a boson mass $m=10^{-22}\mathrm{eV/c^2}$ and a soliton core $r_c \approx 0.7$kpc, comparable to typical half-light radii of Local Group dwarf galaxies. We simulate the rapid gas removal from the center of the soliton by a concentric external time-dependent Hernquist potential. We explore two scenarios of feedback blowouts: i) a massive single burst, and ii) multiple consecutive blowouts injecting the same total energy to the system, including various magnitudes for the blowouts in both scenarios. In all cases, we find one single blowout has a stronger effect on reducing the soliton central density. Feedback leads to central soliton densities that oscillate quasi-periodically for an isolated soliton and stochastically for a SFDM halo. The range in the density amplitude depends on the strength of the blowout, however we observe typical variations of a factor of $\geqslant$2. One important consequence of the stochastic fluctuating densities is that, if we had no prior knowledge of the system evolution, we can only know the configuration profile at a specific time within some accuracy. By fitting soliton profiles at different times to our simulated structures, we found the (1-$\sigma$) scatter of their time-dependent density profiles. For configurations within the 1$\sigma$ range, we find the inferred boson mass is typically less than 20\% different from the real value used in our simulations. Finally, we compare the observed dynamical masses of field dwarf galaxies in our Local Group with the implied range of viable solitons from our simulations and find good agreement., Comment: 11 pages, 6 figures, submitted to MNRAS

Fuzzy Dark Matter (FDM), consisting of ultralight bosons, is an intriguing alternative to Cold Dark Matter. Numerical simulations solving the Schr\"odinger-Poisson (SP) equation, which governs FDM dynamics, show that FDM halos consist of a central solitonic core (representing the ground state of the SP equation), surrounded by a large envelope of excited states. Wave interference gives rise to order unity density fluctuations throughout the envelope and causes the soliton to undergo density oscillations and execute a confined random walk in the central region of the halo. The resulting gravitational potential perturbations are an efficient source of dynamical heating. Using high-resolution numerical simulations of a $6.6 \times 10^{9} \rm M_{\odot}$ FDM halo with boson mass, $m_{\rm b}=8 \times 10^{-23} \ \rm eV$, we investigate the impact of this dynamical heating on the structure and kinematics of spheroidal dwarf galaxies of a fixed mass but different initial sizes and ellipticities. The galaxies are set up in equilibrium in the time-and-azimuthally averaged halo potential and evolved for $10 \ \rm Gyr$ in the live FDM halo. We find that they continuously increase their sizes and central velocity dispersions. In addition, their kinematic structures become strongly radially anisotropic, especially in the outskirts. Dynamical heating also causes initially ellipsoidal galaxies to become more spherical over time from the inside out and gives rise to distorted, non-concentric isodensity contours. These tell-tale characteristics of dynamical heating of dwarf galaxies in FDM halos can potentially be used to constrain the boson mass., Comment: 18 pages, 11 figures, accepted for publication in ApJ

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., Comment: 10 pages, 3 figures, submitted to mnras

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., Comment: Version accepted for publication by MNRAS. 19 pages, 16 figures, 2 tables, 5 appendices

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 $10^{10}$ M$_\odot$ 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% 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% 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.

Fuzzy Dark Matter (FDM), consisting of ultralight bosons ($m_{\rm b} \sim 10^{-22}\ \rm eV$), is an intriguing alternative to Cold Dark Matter. Numerical simulations that solve the Schr\"odinger-Poisson (SP) equation show that FDM halos consist of a central solitonic core, which is the ground state of the SP equation, surrounded by an envelope of interfering excited states. These excited states also interfere with the soliton, causing it to oscillate and execute a confined random walk with respect to the halo center of mass. Using high-resolution numerical simulations of a $6.6 \times 10^9 M_{\odot}$ FDM halo with $m_{\rm b} = 8 \times 10^{-23}\ \rm eV$ in isolation, we demonstrate that the wobbling, oscillating soliton gravitationally perturbs nuclear objects, such as supermassive black holes or dense star clusters, causing them to diffuse outwards. In particular, we show that, on average, objects with mass $\lesssim 0.3 \%$ of the soliton mass ($M_{\rm sol}$) are expelled from the soliton in $\sim 3\ \rm Gyr$, after which they continue their outward diffusion due to gravitational interactions with the soliton and the halo granules. More massive objects ($\gtrsim 1 \% M_{\rm sol}$), while executing a random walk, remain largely confined to the soliton due to dynamical friction. We also present an effective treatment of the diffusion, based on kinetic theory, that accurately reproduces the outward motion of low mass objects and briefly discuss how the observed displacements of star clusters and active galactic nuclei from the centers of their host galaxies can be used to constrain FDM., Comment: 20 pages, 9 figures, minor revisions, accepted for publication in ApJ

We present a suite of baryonic cosmological zoom-in simulations of self-interacting dark matter (SIDM) haloes within the ``Feedback In Realistic Environment'' (FIRE) project. The three simulated haloes have virial masses of $\sim 10^{12}\, \text{M}_\odot$ at $z=0$, and we study velocity-independent self-interaction cross sections of 1 and 10 ${\rm cm^2 \, g^{-1}}$. We study star formation rates and the shape of dark matter density profiles of the parent haloes in both cold dark matter (CDM) and SIDM models. Galaxies formed in the SIDM haloes have higher star formation rates at $z\leq1$, resulting in more massive galaxies compared to the CDM simulations. While both CDM and SIDM simulations show diverse shape of the dark matter density profiles, the SIDM haloes can reach higher and more steep central densities within few kpcs compared to the CDM haloes. We identify a correlation between the build-up of the stars within the half-mass radii of the galaxies and the growth in the central dark matter densities. The thermalization process in the SIDM haloes is enhanced in the presence of a dense stellar component. Hence, SIDM haloes with highly concentrated baryonic profiles are predicted to have higher central dark matter densities than the CDM haloes. Overall, the SIDM haloes are more responsive to the presence of a massive baryonic distribution than their CDM counterparts., Comment: 10 pages, 5 figures. Accepted for publication in MNRAS

Using the phat-ELVIS suite of Milky Way-size halo simulations, we show that subhalo orbital pericenters, $r_{\rm peri}$, correlate with their dark matter halo structural properties. Specifically, at fixed maximum circular velocity, $V_{\rm max}$, subhalos with smaller $r_{\rm peri}$ are more concentrated (have smaller $r_{\rm max}$ values) and have lost more mass, with larger peak circular velocities, $V_{\rm peak}$, prior to infall. These trends provide information that can tighten constraints on the inferred $V_{\rm max}$ and $V_{\rm peak}$ values for known Milky Way satellites. We illustrate this using published pericenter estimates enabled by Gaia for the nine classical Milky Way dwarf spheroidal satellites. The two densest dSph satellites (Draco and Ursa Minor) have relatively small pericenters, and this pushes their inferred $r_{\rm max}$ and $V_{\rm max}$ values lower than they would have been without pericenter information. For Draco, we infer $V_{\rm max} = 23.5 \, \pm 3.3$ km s$^{-1}$ (compared to $27.3 \, \pm 7.1$ km s$^{-1}$ without pericenter information). Such a shift exacerbates the traditional Too Big to Fail problem. Draco's peak circular velocity range prior to infall narrows from $V_{\rm peak} = 21 - 49$ km s$^{-1}$ without pericenter information to $V_{\rm peak} = 25-37$ km s$^{-1}$ with the constraint. Over the full population of classical dwarf spheroidals, we find no correlation between $V_{\rm peak}$ and stellar mass today, indicative of a high level of stochasticity in galaxy formation at stellar masses below $\sim 10^7$ M$_\odot$. As proper motion measurements for dwarf satellites become more precise, they should enable useful priors on the expected structure and evolution of their host dark matter subhalos., Comment: 6 Pages, 4 figures, submitted to MNRAS

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., Comment: This version has been accepted for publication in Physical Review D

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

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

We perform high-resolution simulations of a MW-like galaxy in a self-interacting cold dark matter model with elastic cross section over mass of $1~\rm cm^2/g$ (SIDM) and compare to a model without self-interactions (CDM). We run our simulations with and without a time-dependent embedded potential to capture effects of the baryonic disk and bulge contributions. The CDM and SIDM simulations with the embedded baryonic potential exhibit remarkably similar host halo profiles, subhalo abundances and radial distributions within the virial radius. The SIDM host halo is denser in the center than the CDM host and has no discernible core, in sharp contrast to the case without the baryonic potential (core size $\sim 7 \, \rm kpc$). The most massive subhalos (with $V_{\mathrm{peak}}> 20 \, \rm km/s$) in our SIDM simulations, expected to host the classical satellite galaxies, have density profiles that are less dense than their CDM analogs at radii less than 500 pc but the deviation diminishes for less massive subhalos. With the baryonic potential included in the CDM and SIDM simulations, the most massive subhalos do not display the too-big-to-fail problem. However, the least dense among the massive subhalos in both these simulations tend to have the smallest pericenter values, a trend that is not apparent among the bright MW satellite galaxies., Comment: 5 figures, submitted to MNRAS letters

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., Comment: 22 pages, 16 figures

We present a suite of FIRE-2 cosmological zoom-in simulations of isolated field dwarf galaxies, all with masses of $M_\mathrm{halo} \approx 10^{10}\,$M$_\odot$ at $z=0$, across a range of dark matter models. For the first time, we compare how both self-interacting dark matter (SIDM) and/or warm dark matter (WDM) models affect the assembly histories as well as the central density structure in fully hydrodynamical simulations of dwarfs. Dwarfs with smaller stellar half-mass radii (r$_{1/2}<500$ pc) have lower $\sigma_\star/V_\mathrm{max}$ ratios, reinforcing the idea that smaller dwarfs may reside in halos that are more massive than is naively expected. The majority of dwarfs simulated with self-interactions actually experience contraction of their inner density profiles with the addition of baryons relative to the cores produced in dark-matter-only runs, though the simulated dwarfs are always less centrally dense than in $\Lambda$CDM. The V$_{1/2}-$r$_{1/2}$ relation across all simulations is generally consistent with observations of Local Field dwarfs, though compact objects such as Tucana provide a unique challenge. Spatially-resolved rotation curves in the central regions ($<400$ pc) of small dwarfs could provide a way to distinguish between CDM, WDM, and SIDM, however: at the masses probed in this simulation suite, cored density profiles in dwarfs with small r$_{1/2}$ values can only originate from dark matter self-interactions., Comment: 16 pages, 12 figures. V2: matches version accepted by MNRAS

The quenching "maintenance'" and related "cooling flow" problems are important in galaxies from Milky Way mass through clusters. We investigate this in halos with masses $\sim 10^{12}-10^{14}\,{\rm M}_{\odot}$, using non-cosmological high-resolution hydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model. We specifically focus on physics present without AGN, and show that various proposed "non-AGN" solution mechanisms in the literature, including Type Ia supernovae, shocked AGB winds, other forms of stellar feedback (e.g. cosmic rays), magnetic fields, Spitzer-Braginskii conduction, or "morphological quenching" do not halt or substantially reduce cooling flows nor maintain "quenched" galaxies in this mass range. We show that stellar feedback (including cosmic rays from SNe) alters the balance of cold/warm gas and the rate at which the cooled gas within the galaxy turns into stars, but not the net baryonic inflow. If anything, outflowing metals and dense gas promote additional cooling. Conduction is important only in the most massive halos, as expected, but even at $\sim 10^{14}\,{\rm M}_{\odot}$ reduces inflow only by a factor $\sim 2$ (owing to saturation effects and anisotropic suppression). Changing the morphology of the galaxies only slightly alters their Toomre-$Q$ parameter, and has no effect on cooling (as expected), so has essentially no effect on cooling flows or maintaining quenching. This all supports the idea that additional physics, e.g., AGN feedback, must be important in massive galaxies., Comment: 16 pages, 12 figures

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

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

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 analysis of the rotation curves of a large sample of galaxies with very diverse stellar properties reveal a relation between the radial acceleration purely due to the baryonic matter and the one inferred directly from the observed rotation curves. Assuming the dark matter (DM) exists, this acceleration relation is tantamount to an acceleration relation between DM and baryons. This leads us to a universal maximum acceleration for all halos. Using the latter in DM profiles that predict inner cores implies that the central surface density $\mu_{DM} = \rho_s r_s$ must be a universal constant, as suggested by previous studies in selected galaxies, revealing a strong correlation between the density $\rho_s$ and scale $r_s$ parameters in each profile. We then explore the consequences of the constancy of $\mu_{DM}$ in the context of the ultra-light scalar field dark matter model (SFDM). We find that for this model $\mu_{DM} = 648 \, M_\odot {\rm pc}^{-2}$, and that the so-called WaveDM soliton profile should be an universal feature of the DM halos. Comparing with data from the Milky Way and Andromeda satellites, we find that they are consistent with a boson mass of the scalar field particle of the order of $10^{-21} \, {\rm eV}/c^2$, which puts the SFDM model in agreement with recent cosmological constraints., Comment: 6 pages, 1 figure

One alternative to the cold dark matter (CDM) paradigm is the scalar field dark matter (SFDM) model, which assumes dark matter is a spin-0 ultra-light scalar field (SF) with a typical mass $m\sim10^{-22}\mathrm{eV}/c^2$ and positive self-interactions. Due to the ultra-light boson mass, the SFDM could form Bose-Einstein condensates (BEC) in the very early Universe, which are interpreted as the dark matter haloes. Although cosmologically the model behaves as CDM, they differ at small scales: SFDM naturally predicts fewer satellite haloes, cores in dwarf galaxies and the formation of massive galaxies at high redshifts. The ground state (or BEC) solution at zero temperature suffices to describe low-mass galaxies but fails for larger systems. A possible solution is adding finite-temperature corrections to the SF potential which allows combinations of excited states. In this work, we test the finite-temperature multistate SFDM solution at galaxy cluster scales and compare our results with the Navarro-Frenk-White (NFW) and BEC profiles. We achieve this by fitting the mass distribution of 13 \textit{Chandra} X-ray clusters of galaxies, excluding the region of the brightest cluster galaxy. We show that the SFDM model accurately describes the clusters' DM mass distributions offering an equivalent or better agreement than the NFW profile. The complete disagreement of the BEC model with the data is also shown. We conclude that the theoretically motivated multistate SFDM profile is an interesting alternative to empirical profiles and ad hoc fitting-functions that attempt to couple the asymptotic NFW decline with the inner core in SFDM., Comment: 17 pages, 8 figures. This version matches the published version in MNRAS

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 modeled explicitly. The additional physics have no systematic effect on galactic star formation rates (SFRs) . In contrast, removing stellar feedback leads to SFRs being over-predicted by factors of $\sim 10 -100$. Without feedback, neither galactic winds nor volume filling hot-phase gas exist, and discs tend to runaway collapse to ultra-thin scale-heights with unphysically dense clumps congregating at the galactic center. 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_{\odot}$, the additional physics investigated here (MHD, conduction, viscosity, metal diffusion) have only weak ($\sim10\%$-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 AGN feedback is not included in the present work.

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

Recent observations of tidal debris around galaxies have revealed that the structural properties of the spheroidal components of tidally disturbed galaxies are similar to those found in non-interacting early-type galaxies(ETGs), likely due to minor merging events that do not strongly affect the bulge region or to major mergers that happened a long time ago. We show that independently of merger events, tidal features like shells or rings can also arise if the the dark matter is an ultra light scalar field of mass ~10$^{-22}$eV/c$^2$. In the scalar field dark matter (SFDM) model the small mass precludes halo formation below ~10$^8$ M$_{\odot}$ reducing the number of small galaxies today, it produces shallow density profiles due to the uncertainty principle in contrast to the steep profiles found in the standard cold dark matter (CDM) paradigm, in addition to the usual soliton solution there exists dark matter haloes in multistates, characterized by ripples in their density profiles, which are stable provided that most of the halo mass resides in the ground state. We use the hydrodynamics code ZEUS to track the gas evolution in a background potential given by a superposition of the ground and first excited state of the scalar field, we study this configuration when it is initially unstable (excited state more massive than ground state) but by a population inversion in the states it eventually becomes stable, this could happen when haloes decoupled from the expansion of the universe and collapse to reach a state of equilibrium. We found that tidal structures like rings are formed at a particular radii as a direct consequence of the wavelike structure of the dark matter halo(abridged), Comment: 11 pages, 4 figures. Submitted to MNRAS

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

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., Comment: To be published in the Proceedings of the IV International Meeting on Gravitation and Cosmology, 32 pages, 7 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. One mechanism to reconcile the simulations with the observed data is the feedback from star formation, this might be successful in isolated dwarf galaxies but its success in 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, the feedback recipes usually require to remove a large quantity of baryons from the center of galaxies, however, some 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 at 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 that it can satisfy the two desired requirements and that the rotation curve profile is not longer universal., Comment: 11 pages, 3 figures, this version matches the one accepted for publication in the Astrophysical Journal

In the last time the cold dark matter (CDM) model has suggested more and more that it is not able to describe all the properties of nearby galaxies that can be observed in great detail as well as that it has some problems in the mechanism by which matter is more rapidly gathered into large-scale structure such as galaxies and clusters of galaxies. In this work we revisit an alternative model, the scalar field dark matter (SFDM) model, which proposes that the galactic haloes form by condensation of a scalar field (SF) very early in the Universe, i.e., in this model the haloes of galaxies are astronomical Bose-Einstein Condensate drops of SF. On the other hand, large-scale structures like clusters or superclusters of galaxies form similar to the $\Lambda$CDM model, by hierarchy, thus all the predictions of the $\Lambda$CDM model at cosmological scales are reproduced by SFDM. This model predicts that all galaxy haloes must be very similar and exist for higher redshifts than in the $\Lambda$CDM model. In the first part of this review we revisit the cosmological evolution of SFDM model with a scalar potential $m^2\Phi^2/2+\lambda\Phi^4/4$ with two different frameworks: the field and fluid approach. The scalar fluctuations have an oscillating growing mode and therefore, this kind of dark matter could lead to the early formation of gravitational structures in the Universe. In the last part, we study the core central density profiles of BEC dark matter haloes and fit high-resolution rotation curves of low surface brightness galaxies. The mean value of the logarithmic inner density slopes is $\alpha = - $0.27 $\pm$ 0.18 and we show 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 $\Lambda$CDM model is currently facing the SFDM model can be a worthy alternative to keep exploring further., Comment: To be published in Proceedings of the XIII Mexican Workshop on Particles and Fields, 20 pages, 10 figures

The scalar field dark matter (SFDM) model proposes that galaxies form by condensation of a scalar field (SF) very early in the universe forming Bose-Einstein Condensates (BEC) drops, i.e., in this model haloes of galaxies are gigantic drops of SF. Here big structures form like in the LCDM model, by hierarchy, thus all the predictions of the LCDM model at big scales are reproduced by SFDM. This model predicts that all galaxies must be very similar and exist for bigger redshifts than in the LCDM model. In this work we show that BEC dark matter haloes fit high-resolution rotation curves of a sample of thirteen low surface brightness galaxies. We compare our fits to those obtained using a Navarro-Frenk-White and Pseudo-Isothermal (PI) profiles and found a better 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 can be a worthy alternative to keep exploring further., Comment: Submitted to MNRAS, 9 pages, 32 Figures, 2 Tables.The paper with better resolution figures can be downloaded at "http://estudiantes.fis.cinvestav.mx/vrobles/SFDMfile.pdf