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.

A scalar field $\phi$ endowed with a trigonometric potential has been proposed to play the role of Dark Matter. A deep study of the cosmological evolution of linear perturbations, and its comparison to the Cold Dark Matter (CDM) and Fuzzy Dark Matter (FDM) cases (scalar field with quadratic potential), reveals an enhancement in the amplitude of the mass power spectrum for large wave numbers due to the nonlinearity of the axion-like potential. For the first time, we study the scale-dependence on physical quantities such as the growth factor $D_k$, the velocity growth factor $f_k$, and $f_k \sigma_8$. We found that for $z<10$, all these quantities recover the CDM evolution, whereas for high redshift there is a clear distinction between each model (FDM case, and axion-like potential) depending on the wavenumber $k$ and on the decay parameter of the axion-like potential as well. A semi-analytical Halo Mass Function is also revisted, finding a suppression of the number of low mass halos, as in the FDM case, but with a small increment in the amplitude of the variance and halo mass function due to the nonlinearity of the axion-like potential. Finally, we present constraints on the axion mass $m_{\phi}\geq 10^{-24}$eV and the axion decay parameter is not constrained within the prior $0\leq \lambda \leq 10^4$ by using data of the Planck Collaboration 2015.

We present a general strategy to solve the stationary Schr\"odinger-Poisson (SP) system of equations for multistates with axial symmetry. The approach allows us to obtain the well known single and multistate solutions with spherical symmetry, Newtonian multistate $\ell-$boson stars and axially symmetric multistate configurations. For each case we construct particular examples that illustrate the method, whose stability properties are studied by numerically solving the time-dependent SP system. Among the stable configurations there are the mixed-two-state configurations including spherical and dipolar components, which might have an important value as potential anisotropic dark matter halos in the context of ultralight bosonic dark matter scenarios. This is the reason why we also present a possible process of formation of these mixed-two-state configurations that could open the door to the exploration of more general multistate structure formation scenarios.

Dark matter models in which the constituent particle is an ultra-light boson have become part of the mainstream discussion in cosmology and astrophysics. At the classical level, the models are represented by the dynamics of a (real or complex) scalar field endowed with a potential that contains its self-interactions, and for this reason are generically known as scalar field dark matter models. Here, we revisit the properties of such a model with a cosh potential and compare it with other known examples in the literature. Within the cosmological context, the self-interaction in the potential induces a radiation-like behavior at early times of the scalar field density, which is followed by a proper matter-like behavior at the onset of rapid field oscillations around the minimum of the potential. The solutions are found by numerical means, and from them we obtain information about the cosmological observables up to the level of linear density perturbations. We also study the general properties of self-gravitating objects in the non-relativistic limit and determine the role in them of the self-interaction obtained from the cosh potential. An overall conclusion is that, for the range of values in its parameters allowed by different constraints, a cosh potential behaves almost indistinguishable from the simpler quadratic one, which also means that the two models suffer the same tight constraints from cosmological and astrophysical observations.

Wave Dark Matter (WaveDM) has recently gained attention as a viable candidate to account for the dark matter content of the Universe. In this paper we explore the extent to which dark matter halos in this model, and under what conditions, are able to reproduce strong lensing systems. First, we analytically explore the lensing properties of the model -- finding that a pure WaveDM density profile, a soliton profile, produces a weaker lensing effect than other similar cored profiles. Then we analyze models with a soliton embedded in an NFW profile, as has been found in numerical simulations of structure formation. We use a benchmark model with a boson mass of $m_a=10^{-22}{\rm eV}$, for which we see that there is a bi-modality in the contribution of the external NFW part of the profile, and actually some of the free parameters associated with it are not well constrained. We find that for configurations with boson masses $10^{-23}$ -- $10^{-22}{\rm eV}$, a range of masses preferred by dwarf galaxy kinematics, the soliton profile alone can fit the data but its size is incompatible with the luminous extent of the lens galaxies. Likewise, boson masses of the order of $10^{-21}{\rm eV}$, which would be consistent with Lyman-$\alpha$ constraints and consist of more compact soliton configurations, necessarily require the NFW part in order to reproduce the observed Einstein radii. We then conclude that lens systems impose a conservative lower bound $m_a > 10^{-24}$ and that the NFW envelope around the soliton must be present to satisfy the observational requirements.

Nonlinearities in a realistic axion field potential may play an important role in the cosmological dynamics. In this paper we use the Boltzmann code CLASS to solve the background and linear perturbations evolution of an axion field and contrast our results with those of CDM and the free axion case. We conclude that there is a slight delay in the onset of the axion field oscillations when nonlinearities in the axion potential are taken into account. Besides, we identify a tachyonic instability of linear modes resulting in the presence of a bump in the power spectrum at small scales. Some comments are in turn about the true source of the tachyonic instability, how the parameters of the axionlike potential can be constrained by Ly-$\alpha$ observations, and the consequences in the stability of self-gravitating objects made of axions.

It has been suggested that the internal dynamics of dwarf spheroidal galaxies (dSphs) can be used to test whether or not ultralight axions with $m_a\sim 10^{-22}\text{eV}$ are a preferred dark matter candidate. However, comparisons to theoretical predictions tend to be inconclusive for the simple reason that while most cosmological models consider only dark matter, one observes only baryons. Here we use realistic kinematic mock data catalogs of Milky Way dSph's to show that the "mass-anisotropy degeneracy" in the Jeans equations leads to biased bounds on the axion mass in galaxies with unknown dark matter halo profiles. In galaxies with multiple chemodynamical components this bias can be partly removed by modelling the mass enclosed within each subpopulation. However, analysis of the mock data reveals that the least-biased constraints on the axion mass result from fitting the luminosity-averaged velocity dispersion of the individual chemodynamical components directly. Applying our analysis to two dSph's with reported stellar subcomponents, Fornax and Sculptor, and assuming that the halo profile has not been acted on by baryons, yields core radii $r_{c}>1.5$ kpc and $r_c> 1.2$ kpc respectively, and $m_a<0.4\times 10^{-22}\text{eV}$ at 97.5\% confidence. These bounds are in tension with the number of observed satellites derived from simple (but conservative) estimates of the subhalo mass function in Milky Way-like galaxies. We discuss how baryonic feedback might affect our results, and the impact of such a small axion mass on the growth of structures in the Universe.

As we are entering the era of precision cosmology, it is necessary to count on accurate cosmological predictions from any proposed model of dark matter. In this paper we present a novel approach to the cosmological evolution of scalar fields that eases their analytic and numerical analysis at the background and at the linear order of perturbations. We apply the method to a scalar field endowed with a quadratic potential and revisit its properties as dark matter. Some of the results known in the literature are recovered, and a better understanding of the physical properties of the model is provided. It is shown that the Jeans wavenumber defined as $k_J = a \sqrt{mH}$ is directly related to the suppression of linear perturbations at wavenumbers $k>k_J$. We also discuss some semi-analytical results that are well satisfied by the full numerical solutions obtained from an amended version of the CMB code CLASS. Finally we draw some of the implications that this new treatment of the equations of motion may have in the prediction for cosmological observables.

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.

Compact objects, like neutron stars and white dwarfs, may accrete dark matter, and then be sensitive probes of its presence. These compact stars with a dark matter component can be modeled by a perfect fluid minimally coupled to a complex scalar field (representing a bosonic dark matter component), resulting in objects known as fermion-boson stars. We have performed the dynamical evolution of these stars in order to analyze their stability, and to study their spectrum of normal modes, which may reveal the amount of dark matter in the system. Their stability analysis shows a structure similar to that of an isolated (fermion or boson) star, with equilibrium configurations either laying on the stable or on the unstable branch. The analysis of the spectrum of normal modes indicates the presence of new oscillation modes in the fermionic part of the star, which result from the coupling to the bosonic component through the gravity.

We study the production of primordial black holes (PBHs) during the preheating stage that follows a chaotic inflationary phase. The scalar fields present in the process are evolved numerically using a modified version of the HLATTICE code. From the output of the numerical simulation, we compute the probability distribution of curvature fluctuations, paying particular attention to sub-horizon scales. We find that in some specific models these modes grow to large amplitudes developing highly non-Gaussian probability distributions. We then calculate PBH abundances using the standard Press-Schechter criterion and find that overproduction of PBHs is likely in some regions of the chaotic preheating parameter space.

We study in detail the properties of gravitationally-bounded multi-state configurations, made of spin-zero bosons, in the Newtonian regime. We show that the properties of such configurations, in particular their stability, depend upon how the particles are distributed in the different states they are composed of. Numerical techniques are used to distinguish between stable and unstable solutions, and to determine the final configurations they evolve towards to. Multi-state equilibrium configurations can be used as models of galactic halos made of scalar field dark matter, whose rotation curves appear more realistic than in the case of single-state configurations.

Standard thermodynamical results of ideal Bose gases are used to study the possible formation of a cosmological Bose-Einstein condensate in Scalar Field Dark Matter models; the main hypothesis is that the boson particles were in thermal equilibrium in the early Universe. It is then shown that the only relevant case needs the presence of both particles and anti-particles, and that it corresponds to models in which the bosonic particle is very light. Contrary to common wisdom, the condensate should be a relativistic phenomenon. Some cosmological implications are discussed in turn.

We construct an explicit scenario whereby the same material driving inflation in the early Universe can comprise dark matter in the present Universe, using a simple quadratic potential. Following inflation and preheating, the density of inflaton/dark matter particles is reduced to the observed level by a period of thermal inflation, of a duration already invoked in the literature for other reasons. Within the context of the string landscape, one can further argue for a non-zero vacuum energy of this field, thus unifying inflation, dark matter and dark energy into a single fundamental field.

We consider the conditions needed to unify the description of dark matter, dark energy and inflation in the context of the string landscape. We find that incomplete decay of the inflaton field gives the possibility that a single field is responsible for all three phenomena. By contrast, unifying dark matter and dark energy into a single field, separate from the inflaton, appears rather difficult.