Constraining Bose-Einstein condensate dark matter models with the galaxy rotation curves of the SPARC sample

Abstract

If dark matter did exist in the form of a self-interacting boson gas, when its temperature became lower than the critical, density dependent temperature, a phase transition to a Bose-Einstein Condensate did occur in the early Universe. The presence of the condensate dark matter leads to specific signatures in the galaxy rotation, and in the dark matter density distribution. In particular Bose-Einstein Condensate (BEC) dark matter models predict a finite and well defined radius of the dark matter halo, whose numerical value depends on the mass of the dark matter particle, and of the scattering length.

We compare the theoretical predictions of the functional form of the rotation curves in the slowly rotating BEC models with the SPARC sample of measured rotation curves, by using genetic algorithms, to fit the observational data, and to obtain estimates of the relevant physical parameters of the BEC dark matter halos (central density, angular velocity and static radius).

The density profiles of the dark matter distribution are also considered, and it follows that the presence of the condensate dark matter could also provide an alternative solution for the core/cusp problem.

Authors

M. Crăciun
(Tiberiu Popoviciu Institute of Numerical Analysis)

T. Harko
(Department of Physics, Babes-Bolyai University, Cluj-Napoca, Romania and
School of Physics, Sun Yat-Sen University Guangzhou, People’s Republic of China)

Keywords

Cosmology; galactic astronomy; dark matter; Bose-Einstein condensation

References

See the expanding block below.

Paper coordinates

M. Crăciun, T. Harko, Constraining Bose-Einstein condensate dark matter models with the galaxy rotation curves of the SPARC sample, Romanian Astronomical Journal, 29 (2019) no.2, 101-118.

PDF

About this paper

Journal

Romanian Astronomical Journal

Publisher Name

Romanian Academy Publishing House (Editura Academiei Romane)

DOI

not available yet.

Print ISSN

1220-5168  

Online ISSN

2285-3758

Google Scholar Profile
Constraining Bose-Einstein Condensate Dark Matter Models with the Galaxy Rotation Curves of the SPARC Sample

Constraining Bose-Einstein Condensate Dark Matter Models with the Galaxy Rotation Curves of the SPARC Sample

Maria CRĂCIUN 1    Tiberiu HARKO 2,3

1 ‘Tiberiu Popoviciu’ Institute of Numerical Analysis, Romanian Academy Fântânele St. 57, 400320, Cluj-Napoca, Romania Email:craciun@ictp.acad.ro
2 Department of Physics, Babes-Bolyai University Kogalniceanu St., Cluj-Napoca 400084, Romania
3 School of Physics, Sun Yat-Sen University Guangzhou 510275, People’s Republic of China Email:t.harko@ucl.ac.uk

Abstract

If dark matter did exist in the form of a self-interacting boson gas, when its temperature became lower than the critical, density dependent temperature, a phase transition to a Bose-Einstein Condensate did occur in the early Universe. The presence of the condensate dark matter leads to specific signatures in the galaxy rotation, and in the dark matter density distribution. In particular Bose-Einstein Condensate (BEC) dark matter models predict a finite and well defined radius of the dark matter halo, whose numerical value depends on the mass of the dark matter particle, and of the scattering length. We compare the theoretical predictions of the functional form of the rotation curves in the slowly rotating BEC models with the SPARC sample of measured rotation curves, by using genetic algorithms, to fit the observational data, and to obtain estimates of the relevant physical parameters of the BEC dark matter halos (central density, angular velocity and static radius). The density profiles of the dark matter distribution are also considered, and it follows that the presence of the condensate dark matter could also provide an alternative solution for the core/cusp problem.

keywords:
Cosmology – Galactic astronomy – Dark matter – Bose-Einstein Condensation

1 Introduction

The standard ΛCDM cosmological paradigm has proven to be exceptionally successful in describing the major properties of the Universe, as well as its large scale structure (Ade et al., 2016). However, despite its success, two fundamental open questions remain unanswered by the ΛCDM model, namely, the precise nature of the dark matter, and of the dark energy, respectively. Moreover, there is a large discrepancy between the theoretically expected and the observed values of the cosmological constant Λ (Bull et al., 2016).

The comparison of the theoretical models with observations have also generated a number of still unresolved tensions on the galactic scales (Bullock and Boylan-Kolchin, 2017). There are (at least) three problems that have been intensively investigated. The first one is the missing satellite problem (Moore et al., 1999; Klypin et al., 1999), which follows from the significant difference between the observed number of satellite galaxies of the Milky Way and the number of substructures that appear in ΛCDM cosmological simulations. The second problem is the too big to fail problem, which is a consequence of the observations of the low central densities of the Milky Way satellites (Boylan-Kolchin, Bullock and Kaplinghat, 2011).

A third problem that confronts present day galactic astrophysics is the so-called ‘core/cusp’ problem. The dark matter-only simulations in the ΛCDM model predict cuspy density profiles for dark matter haloes, that is, profiles with steep central slopes, given by (Navarro, Frenk and White, 1997),

ρ=ρs(r/rs)(1+r/rs)2, (1)

which behave near the origin as ρr1. On the other hand, the observations of the galaxy rotation curves indicate the existence of flat or cored profiles in the inner regions of the halos (Oh et al., 2015). There are several solutions proposed to this problem, like, for example, taking into account in the numerical simulations the effects of the baryonic feedback (Di Cintio et al., 2014). Supernovae explosions may blow large amounts of gas in the central galactic region of a galaxy, leading to a fluctuation in the gravitational potential of the gas that may also influence the dark matter distribution.

However, from a physical point of view the core-cusp problem can be easily solved (Harko, 2011) by assuming that the dark matter is in the form of a BEC (Boehmer and Harko, 2007), that is, it consists of massive bosons that gain a repulsive interaction by occupying the same ground energy state. The Bose-Einstein Condensation process takes place when the temperature T of the boson gas becomes smaller than the critical temperature Tcr, given by (Pethick and Smith, 2008)

Tcr=2π2ρcr2/3ζ2/3(3/2)m5/3kB, (2)

where m is the mass of the particle in the condensate, ρcr is the critical transition density, kB is Boltzmann’s constant, and ζ is the Riemmann zeta function, respectively. The Bose-Einstein condensation process was observed experimentally in 1995 in dilute alkali gases, cooled to extremely low temperatures, and confined in a magnetic trap (Anderson et al., 1995; Bradley et al., 1995; Davis et al., 1995).

The Bose-Einstein condensate dark matter is described by the nonlinear Gross-Pitaevskii equation, with the gravitational potential playing the role of the confining potential. In the Thomas-Fermi approximation, condensate dark matter can be described by a polytropic equation of state with polytropic index n=1 (Boehmer and Harko, 2007). The properties of the galactic clusters in the BEC dark matter model were considered in (Harko et al., 2015). The study of the rotating BEC halos was performed in (Zhang et al., 2018) by assuming a rigid body rotation for the halo. By using the hydrodynamic representation of the Gross-Pitaevskii equation one can obtain the basic equation describing the density distribution of the rotating condensate. The general solutions for the condensed dark matter density can then be found, as well as the general representations for the mass distribution, boundary (radius), and tangential velocity, respectively. In the first order of approximation, and under the assumption of slow rotation, explicit representations for the radius, mass, and tangential velocity can be obtained, which generalizes to the rotating case the static dark matter distribution expressions.

From the BEC dark matter model one can obtain a systematic description of the rotation velocity of massive particles in condensate galactic halos, as well of the density distribution of the halo. It is the goal of the present study to compare the predictions of the BEC dark matter models with the SPARC sample of high-quality rotation curves (Lelli, McGaugh and Schombert, 2016), and to search for the specific signatures of the presence of condensate dark matter. After performing several cuts to the sample, we keep a number of 139 galaxies, each modeled by a gaseous disc, a stellar disc, a dark matter halo, and a stellar bulge, if appropriate. For these models we apply a fitting technique based on genetic algorithms, and we obtain the physical parameters for each galaxy (central density, rotational velocity, and halo radius).

The present paper is organized as follows. In Section 2 we briefly present the basic properties of the static and rotating Bose-Einstein Condensatedark matter halos. The comparison of the theoretical predictions of the BEC model with the high quality rotation curves data of the SPARC sample is performed in Section 3. We discuss and conclude our results in Section 4.

2 Bose-Einstein Condensate Dark Matter Halos in Slow Rotation

In the following we assume that the dark matter halos consist of a strongly - coupled, but dilute bosonic Bose-Einstein Condensate at absolute zero temperature (Membrado, Pacheco and Sanudo, 1989). Therefore almost all dark matter particles must be in the condensate. In a cold and dilute fluid, only low energy binary collisions are important. These collisions can be described, independently of the details of the two-body potential, with the help of a single parameter, the s-wave scattering length a. Hence, the interaction potential can be replaced with an effective interaction VI(rr)=λδ(rr), with the coupling constant λ related to the scattering length a through the relation λ=4π2a/m. The ground state properties of the dark matter are described by the mean-field Gross-Pitaevskii (GP) equation, which for the case of rotating dark matter halos can be derived from the GP energy functional, which takes the form

E[ψ]=[22m|ψ(r)|2+mVrot+U02|ψ(r)|4]𝑑r (3)
12Gm2|ψ(r)|2|ψ(r)|2|rr|𝑑r𝑑r, (4)

where ψ(r) is the wave function of the condensate, and U0=4π2a/m (Pethick and Smith, 2008). The first term in the GP energy functional is the quantum pressure, the second is the rotational energy Vrot=ω2r2/2, where ω is the angular velocity of the condensate, the third is the interaction energy, while the fourth term gives the gravitational potential energy. We also introduce the mass density of the condensate, defined as ρ(r)=m|ψ(r)|2. The wave function obeys the normalization condition is N=|ψ(r)|2𝑑r, where N is the total number of condensate dark matter particles. The variational procedure

δE[ψ]μδ|ψ(r)|2𝑑r=0, (5)

gives the GP equation as (Boehmer and Harko, 2007)

22mΔψ(r)+mVrot(r)ψ(r)+mV(r)ψ(r)+U0|ψ(r)|2ψ(r)=μψ(r), (6)

where μ=constant is the chemical potential. The gravitational potential V(r) satisfies the Poisson equation

ΔV=4πGρ. (7)

In the case of the galactic dark matter halos the number of particles in the condensate is very large, and therefore the quantum pressure term makes a significant contribution only near the boundary of the galaxy. Moreover, it is much smaller than the interaction energy term. Thus for condensate galactic halos, the quantum pressure term can be neglected, and this approximation is called the Thomas-Fermi approximation. For a very large number of particles, when N, the Thomas-Fermi approximation becomes exact (Pethick and Smith, 2008). Hence in this approximation we obtain

ρ(r)=mU0[μmVrot(r)mV(r)]. (8)

Then the Poisson equation gives the density profile of the dark matter halo as (Zhang et al., 2018)

Δρ+k2(ρω28πG)=0, (9)

where we have denoted

k=4πGm2U0. (10)

From a mathematical point of view Eq. (9) represents a Helmholtz type equation.

2.1 Static Bose-Einstein Condensate dark matter halos

In the case when the rotation of the condensate can be neglected, ω=0, the density distribution ρ of the static gravitationally bounded dark matter Bose-Einstein condensate is obtained as (Boehmer and Harko, 2007)

ρ(r)=ρcsinkrkr, (11)

where ρc is the central density of the condensate, ρc=ρ(0). The density profile has the important property of being non-singular at the galactic center r=0. The mass profile m(r)=4π0rρ(r)r2𝑑r of the galactic dark matter halo can be obtained as

m(r)=4πρrk2(1krcotkr), (12)

with a boundary radius R. At the boundary of the dark matter distribution we have ρ(R)=0, a condition that gives kR=π. Hence for the radius of the condensate dark matter halo we find the expression

R=π2aGm3. (13)

The tangential velocity vtg of a massive test particle in rotational notion in the condensed dark halo is obtained, in the Newtonian approximation, as

vtg2(r)=Gm(r)r=4πGρk2(1krcotkr), (14)

or, equivalently,

vtg2(r)=4GρcR2π[sin(πr/R)πr/Rcos(πrR)]. (15)

The total mass M of the galactic dark matter halo M is found to be

M=4π2(2aGm3)3/2ρc=4πR3ρc, (16)

giving for the mean value ρ of the condensate density the expression ρ=3ρc/π2. Therefore the tangential velocity can also be represented as

vtg2(r)=GMR[sin(πr/R)πr/Rcos(πrR)]. (17)

For r>R the rotation curves are given by the standard Keplerian law, with vtg2=GM/r. The mass of the dark matter particle in the dark matter halo is given by

m=(π22aGR2)1/36.73×102[a(fm)]1/3[R(kpc)]2/3eV. (18)

For a1 fm and R10 kpc, we obtain for the typical mass of the condensate particle a value of the order of m14 meV. For a106 fm, a value of the scattering length that corresponds to values of a specific for laboratory experiments, m1.44 eV. These values are consistent with the limit m<1.87 eV obtained for the mass of the condensate particle from cosmological considerations (Boyanovsky, de Vega and Sanchez, 2008).

2.2 Slowly rotating Bose-Einstein Condensate dark matter halos

In the analysis of the slowly rotating Bose-Einstein Condensate dark matter halos we follow the approach developed in (Zhang et al., 2018). The general solution of Eq. (9) can be obtained as (Zhang et al., 2018)

ρ(r,θ)=ω22πG+l=0A2lj2l(kr)P2l(cosθ), (19)

where by jl(X) we have denoted the spherical Bessel functions, giving the solutions of the radial part of the equation. The solutions of the angular part of the Helmholtz equation are given by the Legendre polynomials Pl(cosθ). The integration constants A2l must be determined from the condition of the continuity of the gravitational potential across the boundary of the galactic condensate dark matter halo. By assuming the condition

ω2G<<Mk3, (20)

which corresponds to the physical case of a slow rotation of the dark matter halo, it turns out that the deviation from a spherically symmetric matter distribution can be considered small. Therefore in the slow rotation approximation we can use with a good approximation perturbative methods, by introducing the perturbation parameter ω2/GMk3.

The boundary of the dark matter halo r0=r0(θ) is defined by the condition ρ(r,θ)=0. In the first order of approximation we obtain for the density distribution the expression Zhang et al. (2018)

ρ(r,θ)=ω22πG+(ρcω22πG)j0(kr)5ω2π12Gj2(kr)P2(cosθ). (21)

Alternatively, the density distribution of the dark matter in the BEC halo can be written as,

ρ(r,θ)ρc=(1Ω2)j0(kr)+Ω2[15π26j2(kr)P2(cosθ)], (22)

where we have denoted

Ω2 = ω22πGρc=0.02386×(ω1016s1)2×(ρc1024g/cm3)1. (23)

The radius of the slowly rotating Bose-Einstein Condensate dark matter halo, defined as the radius where the dark matter density vanishes, can be computed according to the relation (Zhang et al., 2018)

r0(θ) = πk+ω24Gρck[25P2(cosθ)]=πk×{1+Ω2[152P2(cosθ)]}. (24)

The equatorial radius of the galactic dark matter halo can be obtained as

r0(π2)=πk(1+94Ω2). (25)

If Ω=0 we have r0(π/2)=R=π/k, and hence we recover the static case. For slowly rotating dark matter halos the mass distribution inside radius r is given by (Zhang et al., 2018)

m(r) = 4πρck2r×[(1Ω2)sinkrkr(1Ω2)coskr+Ω2(kr)23]. (26)

Due to the presence of the rigid body type rotation, a second term appears in the mass profile, which is not present in the static case. In the first order of approximation the tangential velocity vtg2 of massive test particles rotating in the Bose-Einstein Condensate galactic dark matter is given by

vtg2(r)=4GρcR2π×[(1Ω2)sin(πr/R)πr/R(1Ω2)cosπrR+Ω23(πrR)2],

or, in an alternative representation, in terms of the radius R of the static configuration, by 111This equation corrects Eq. (103) of reference (Zhang et al., 2018), where, due to a typing mistake, the term proportional to (ρc/1024g/cm3)1 is missing in the definition of Ω2.,

vtg2(km2/s2)=80.861×(ρc1024g/cm3)×(Rkpc)2×
[(1Ω2)sin(πr/R)πr/R(1Ω2)cosπrR+Ω23(πrR)2]. (28)

The total mass of the galactic dark matter halo M is defined according to (Zhang et al., 2018)

M(r0)=2π02πsinθdθ0r0ρ(r,θ)r2𝑑r, (29)

where r0 is the radius of the galaxy.

In the first order of ω2/2πG the total mass of the galactic condensate can be obtained as

M(r0)4π2ρck3+2πω2Gk3(π231)=4π2k3ρc[1+(π231)Ω2]. (30)

The total halo mass of a given galaxy is fully determined by the evolutionary processes that led to its formation. The halo mass is not modified by the presence of rotation. But, as it follows from Eq. (30), the rotation of the halo allows the galaxy to have more mass than a static one. Thus the rotation of the dark matter halo becomes a stabilizing factor that prevents gravitational collapse.

The central density ρc of the galaxy is determined by the total mass and angular velocity of the halo according to

ρc=π4R3[M2πω2k3(π231)]. (31)

The mean density of the rotating BEC halo is given by

ρ=MVBE3ρcπ2+(112π2)ω22πG=3ρcπ2[1+π23(112π2)Ω2]<3ρcπ2, (32)

where VBE is the galactic volume. The mean density of the rotating dark matter distribution is smaller than the mean density of the static configuration.

3 Testing the Bose-Einstein Condensate dark matter model with the SPARC high resolution rotation curves

In the present Section we consider a comparison of the theoretical predictions of the Bose-Einstein Condensate dark energy model with the observations of the galactic rotation curves.

3.1 The SPARC sample and the fitting procedure

For the present study we use the rotation curves from the SPARC database (Lelli, McGaugh and Schombert, 2016). This sample consists of 175 high-quality HI/Hα rotation curves, together with accompanying 3.6μm photometry from Spitzer. The SPARC sample contains a large variety of data, extending over 5 orders of magnitude in mass. It also encompasses a variety of morphologies. The considered 3.6 μm surface brightness is a good tracer of the stellar mass, and hence its combination with rotation curve data allows detailed modelling of the different components of a given galaxy.

The SPARC database includes a large variety of information for each galaxy. From the point of view of the test of the BEC dark matter models the most important information is contained in vobs(r), giving the rotation speed as a function of radial coordinate r, together with the corresponding error bars. From the photometrical information one can also obtain the contributions to the rotation speed of the stellar disc and of the stellar bulge, vdisc and vbulge, respectively. The bulge component is included only for galaxies for which the photometry profile significantly differs from an exponential disc profile in the central regions. vdisc and vbulge are obtained by taking a mass-to-light ratio of 1 M/L, and therefore these quantities must be rescaled for other ratios. For each galaxy, the information on the gas surface brightness was converted into a gas contribution to the total velocity vgas.

The distribution of the number of galaxies according to their type in the SPARC sample is represented in Fig. 1.

Refer to caption
Figure 1: The number of galaxies of each type in the SPARC sample.

The histogram of lengths for the galactic velocities in the SPARC sample is represented in Fig. 2.

Refer to caption
Figure 2: Histogram of lengths of rotation velocity curves in the SPARC sample.

We have taken into account for the fits of the theoretical model only the data with length greater than or equal to 7.

In order to test the velocity of the massive particles in the Bose-Einstein Condensate dark matter halo defined by Eq. (2.2) we have fitted the total velocity from SPARC database with a total velocity vtot, which also includes the baryonic matter contribution, and which is given by the relation

vtot=vgas|vgas|+Υd×vdisc|vdisc|+Υb×vbulge|vbulge|+vtg2, (33)

where Υd and Υb are the mass-to-light ratios for the stellar disc and the bulge, respectively. In order to determine the free parameters of the BEC model we have looked for the minimum of the function

χ2(R,ρc,ω,Υd,Υb)=1nki=1n(vtot,ivobs,i)2σi2, (34)

where n is the length of data and k is the number of parameters that must be estimated, i.e, n=4 for bulgeless galaxies and n=5 for galaxies with bulge.

We have used the so-called genetic algorithms to find the parameters that minimizes the objective function (34) for a given set of constraints. The genetic algorithms that belongs to the larger class of evolutionary algorithms are inspired by the Darwinian process of natural selection (‘survival of the fittest’), and they are used to generate solutions to optimization problems by relying on biologically inspired operators such as mutation, crossover and selection.

In order to apply these algorithms one must define first a fitness function, which in the present case is the objective function defined by Eq.  (34). Then we initialize the genetic algorithms by considering an initial population of random candidates for obtaining a solution, and we compute the fitness value for each candidate. Then new populations are iteratively generated by selection, crossover and mutation until the entire process is stopped when the average relative change in the fitness function value is less than 1010.

The total population size values were considered in the range [4000,5000], while for the crossover fraction (the ratio from the population size that will be recombined for the next population) we have taken values in the interval [0.6,0.9]. For each galaxy we have adapted the initial population range. For the maximum number of generations we have adopted the value 1000. In our investigation we have limited the values of Υd and Υb to the range [0.1,5]. We also require that the density of the halo must be a monotonically decreasing function of the radial coordinate in the range r<rmax. On the rotation parameter Ω we have imposed the constraint

Ω2<0.109. (35)

For the lower bound for the central density of the dark matter halo ρc we have adopted the value 1027g/cm3, while the lower bound for the static radius of the dark matter halo was set to rmax/2.

3.2 Fitting results

Refer to caption
Figure 3: The rotation velocity for some galaxies. The observed total velocity is plotted with error bars, the estimated baryonic contributions (disk and gas contributions) with black dashed line, the tangential velocity with blue line and the estimated total velocity with magenta thick line.

The results of the fitting of the rotation curves of a set of 9 galaxies from the SPARC sample with the predictions of the slowly rotating BEC dark matter model are presented in Fig. 3. As one can see, the theoretical model fits well the observational data, despite the apparent complexity of the total rotation curves. For the considered galaxies only NGC1705 and NGC4088 do indicate the presence of flat rotation curves, extending for a large radius. They also show a clear evidence for the Keplerian profile of the baryonic matter distribution. For the other galaxies the baryonic contribution to the total velocity distribution is either increasing, or a constant, indicating a strong effect of the luminous matter on the rotation curves.

In Fig. 4 we have represented the density distribution of the dark matter halo for the galaxies whose rotation velocity curves are represented in Fig. 3. With vertical dotted, respectively dashed line we pointed the values of the estimated static and rotating radius of the dark matter halo R, respectively Rh, and with continuous line the maximum radius of observed data. For all considered cases the density of the BEC halo is a monotonically decreasing function of the radial coordinate r, and it becomes zero for a fixed value of r, which formally defines the radius of the halo. The presence of a well-defined boundary is one of the important characteristics of the BEC dark matter models, and of their density profiles, clearly indicating a finite physical size of the dark matter halo.

Refer to caption
Figure 4: The density of the dark matter halo for the galaxies whose velocity curves are plotted in Fig. 3. With red dottted line and blue dashed line are represented the values of the estimated R, respectively Rh, while with thin black line is represented the maximum value of r (for which we have observed data).
Galaxy R Rh ρc ω Υd χ2
(kpc) (kpc) (1024g/cm3) (1016s1) (M/L)
DDO168 4.281 4.281 1.608 0.000 0.100 2.610
F563-1 16.902 26.298 0.315 1.203 4.446 1.151
F563-V2 7.619 11.854 1.002 2.147 4.047 0.111
F565-V2 8.567 13.329 0.715 1.814 1.575 0.194
F568-1 8.691 13.522 1.323 2.467 2.637 0.540
F571-V1 12.521 19.481 0.319 1.211 1.561 0.176
NGC1705 4.577 7.120 1.885 2.945 2.416 0.455
NGC3741 5.720 8.899 0.703 1.798 2.109 1.310
NGC4088 42.376 65.932 0.174 0.894 0.445 0.563
Table 1: The optimal parameters for the galaxies whose velocity curves are plotted in Fig. 3.

The physical parameters of the BEC dark matter halos as obtained from the fitting of the galaxies presented in Fig. 3 are represented in Table 1. The obtained results predict that the central density of the dark matter distribution is of the order of 10241025 g/cm3, while their rotational velocities is of the order of 1016 s-1. One should also point out that some fitting results indicate the presence of static (nonrotating) galaxies, like, for example, DDO168. This raises the question if indeed all BEC galaxies are rotating, and, if not, of the physical/astrophysical processes that could lead to a global galactic dynamics.

4 Discussions and final remarks

In the present work, we have considered the potential observational signatures of the Bose-Einstein Condensate dark matter model in the high-quality HI/Hα measurements of the rotation curves of the SPARC data. After reconsidering the data of the sample, and excluding all galaxies with a low number of measurements, we have considered in our analysis only 139 of the original 175 galaxies. The galaxies in the SPARC database are modelled as containing a stellar disc, a gaseous disc, a dark matter halo, and (where appropriate) a stellar bulge. Then, we have fitted the data on the rotational velocities with Eq. (2.2), giving the expression of the tangential velocity of massive test particles in the Bose-Einstein Condensate dark matter model. The tangential velocity is fixed in terms of three quantities only: the central density of the halo, the rotation angular velocity of the dark matter distribution, and the radius R of the static configuration.

The main findings of our analysis are as follows:

  • The slowly rotating Bose-Einstein Condensate dark matter model obtained gives good results in the fitting of the data of the SPARC sample. The values of the objective function are smaller than 1 for 61 galaxies with data length greater than or equal to 7 and reasonable values of the optimal parameters. The distribution of the values of the objective function χ2 for χ2<1 is represented in Fig. 5

    Refer to caption
    Figure 5: The histogram of the values of the objective function χ2
  • The density distribution of the BEC dark matter is nonsingular at the galactic center, and it shows a clear cored central profile. It has been already suggested (Harko, 2011) that the condensation of dark matter can solve the core-cusp problem, and hence remove the unphysical features that appear in the Navarro-Frenk-White, ΛCDM type, density profile. Our results based on the SPARC sample clearly indicate the presence of a core at the galactic center, a result which is supported by the astronomical observations.

  • Our results indicate a wide range of variations of the angular velocity of the galactic halo. Fig. 6 shows the distribution of the estimated values of the angular velocity ω. However, for a large number of galaxies the angular velocity turns out to be very close to zero, indicating the static nature of the condensate dark matter halo. This raises the interesting question of the nature, origin and dynamics of the galactic rotation of dark matter, and of its dissipation for certain structures. Is the galactic rotation the result of some perturbations in the late evolutionary periods, or it is an intrinsic property that could be related to the initial collapsing phases?

    Refer to caption
    Figure 6: The distribution of the estimated values of the angular velocity ω.
  • Another fitting parameter, the central density of the condensate dark halo, has also been obtained from our analysis. The distribution of the estimated values of central density ρc is represented in Fig. 7. Even that obtaining the central density of dark matter at the galactic center may prove to be an extremely difficult observational task, the central density distribution of dark matter may provide, at least in principle, a specific signature that may help in discriminating BEC dark matter models from other models.

    Refer to caption
    Figure 7: The distribution of the estimated values of central density ρc.
  • The third fundamental parameter determining the structure and properties of the condensate dark matter halos is the radius R of the static configuration, which is fully determined by the two fundamental physical parameters of the dark matter, the mass of the dark matter particle, and by its scattering length. Hence R should have a universal value in all condensate dark matter models. From our investigation it follows that R is well constrained by the SPARC observational data, but the spread of the obtained values of R across the sample is very broad, spanning roughly one order of magnitude. This result seems to be inconsistent with the existence of a single universal value for R. The BEC dark matter models make a very definite prediction about the existence of a finite, and well defined radius of the dark matter distribution. The determination of R in our investigation is determined by the spatial extension of the analyzed sample, and by the fact that the value of R is strongly correlated by the maximum value of r from the observational data. In Fig. 8 we have represented the estimated values of the radius R of the static dark matter configuration with respect to rmax, for the galaxies for which we obtained χ2<1. These two quantities are strongly correlated, with the Pearson correlation coefficient given by 0.91.

    Refer to caption
    Figure 8: The values of the estimated static radius of the condensate dark matter halo R versus rmax.

    Presently there are no convincing data indicating the existence of a given radius of the galactic dark matter halo, and of its universal nature. In several dark matter models the dark matter distribution extends to infinity. However, if a future detection of the radii of a large sample of galactic objects would be possible, this would allow a definite test of the BEC dark matter model. An observational proof of the infinite extension of the galactic halos would certainly contradict one of the basic predictions of the BEC model.

The existence of BEC dark matter in nature depends essentially on the significant increase in the number of observational data, and of their precision. A universal length scale, determined by the physical properties of the dark matter particle, is one of the main predictions of the theory. On the other hand the outcome of the core/cusp debate may also provide evidence for BEC dark matter. If simulations done in the framework of the ΛCDM model with baryonic effects included could predict cored profiles, then the presence of galactic halos made of BEC may not be necessary. If, on the other hand, the results of the simulations would still indicate cuspy NFW-like galactic halo profiles, then the BEC dark matter model may provide a realistic explanation for the cored profiles observed in a large number of galaxies.

Acknowledgments

We would like to thank Prof. László Gergely for comments and suggestions that helped us to improve our manuscript. T. H. thanks the Yat Sen School of the Sun Yat-Sen University in Guangzhou, P. R. China, for the kind hospitality offered during the preparation of this work.

References

  • Ade et al. (2016) Ade, P. A. R. et al., Planck Collaboration: 2016, Astron. Astrophys. 594, A1.
  • Anderson et al. (1995) Anderson, M. H., Ensher, J. R., Matthews, M. R., Wieman, C. E., Cornell, E. A.: 1995, Science 269, 198.
  • Boehmer and Harko (2007) Boehmer, C. G., Harko, T.: 2007, JCAP 06, 025.
  • Boyanovsky, de Vega and Sanchez (2008) Boyanovsky, D., de Vega H. J., Sanchez, N.: 2008, Phys. Rev. D 77, 043518.
  • Boylan-Kolchin, Bullock and Kaplinghat (2011) Boylan-Kolchin, M., Bullock, J. S., Kaplinghat, M.: 2011, Mon. Not. Roy. Astron. Soc. 415, L40.
  • Bradley et al. (1995) Bradley, C. C., Sackett, C. A., Tollett, J. J., Hulet, R. G.: 1995, Phys. Rev. Lett. 75, 1687.
  • Bull et al. (2016) Bull, P., et al.: 2016, Physics of the Dark Universe, 12, 56.
  • Bullock and Boylan-Kolchin (2017) Bullock, J. S., Boylan-Kolchin, M.: 2017, Annual Review of Astronomy and Astrophysics 55, 343.
  • Davis et al. (1995) Davis, K. B., Mewes, M. O., Andrews, M. R., van Drutten, N. J., Durfee, D. S., Kurn, D. M., Ketterle, W.: 1995, Phys. Rev. Lett. 75, 3969.
  • Di Cintio et al. (2014) Di Cintio, A., Brook, C. B., Dutton, A. A., Maccio, A. V., Stinson, G. S., Knebe, A.: 2014, Mon. Not. Roy. Astron. Soc. 441, 2986.
  • Harko (2011) Harko, T.: 2011, JCAP 1105, 022.
  • Harko et al. (2015) Harko T., Liang P.-X, Liang S.-D., Mocanu G.: 2015, JCAP 11, 027.
  • Klypin et al. (1999) Klypin, A., Kravtsov, A. V., Valenzuela, O., Prada, F.: 1999, Astrophys. J. 522, 82.
  • Lelli, McGaugh and Schombert (2016) Lelli, F., McGaugh, S. S., Schombert, J. M.: 2016, Astronom. J. 152, 157.
  • Membrado, Pacheco and Sanudo (1989) Membrado, M., Pacheco, A. F. and Sanudo, J. S.: 1989, Astron. Astrophys. 217, 92.
  • Moore et al. (1999) Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, J., Tozzi, P.: 1999, Astrophys. J. 524, L19.
  • Navarro, Frenk and White (1997) Navarro, J. F., Frenk, C. S., White, S. D. M.: 1997, Astrophys. J. 490, 493
  • Oh et al. (2015) Oh, S.-H., et al.: 2015, Astronom. J. 149, 180.
  • Pethick and Smith (2008) Pethick, C. J. and Smith, H.: 2008, Bose-Einstein condensation in dilute gases, Cambridge, Cambridge University Press.
  • Zhang et al. (2018) Zhang, X., Chan, M. H., Harko, T., Liang, S.-D., Leung, C. S.: 2018, EPJC 78, 346.

not available yet

2019

Related Posts