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 we find the cusp-core problem and the abundance of small structures of the standard Cold Dark Matter (CDM) model. In this letter we show that mixed state solutions of the low energy and weak gravitational field limit of field equations, interpreted as galactic halo density profiles, could explain the anisotropic distribution of satellite galaxies observed in the Milky Way, M31 and Centaurus A, where satellites seem to concentrate more towards 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 explain the anisotropy of satellite trajectories.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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)

We analyse the velocity dispersion for eight of the Milky Way dwarf spheroidal satellites in the context of finite temperature scalar field dark mater. In this model the finite temperature allows the scalar field to be in configurations that possess excited states, a feature that has proved to be necessary in order to explain the asymptotic rotational velocities found in low surface brightness (LSB) galaxies. In this work we show that excited states are not only important in large galaxies but also have visible effects in dwarf spheroidals. Additionally, we stress that contrary to previous works where the scalar field dark matter haloes are consider to be purely Bose-Einstein condensates, the inclusion of excited states in these halo configurations provides a consistent framework capable of describing LSBs and dwarf galaxies of different sizes without arriving to contradictions within the scalar field dark matter model. Using this new framework we find that the addition of excited states accounts very well for the raise in the velocity dispersion in Milky Way dwarf spheroidal galaxies improving the fit compared to the one obtained assuming all the DM to be in the form of a Bose Einstein Condensate.

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.

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.

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.

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.

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.