We present a comparative confrontation of both the Bose-Einstein Condensate (BEC) and the Navarro-Frenk-White (NFW) dark halo models with galactic rotation curves. We employ 6 High Surface Brightness (HSB), 6 Low Surface Brightness (LSB), and 7 dwarf galaxies with rotation curves falling into two classes. In the first class rotational velocities increase with radius over the observed range. The BEC and NFW models give comparable fits for HSB and LSB galaxies of this type, while for dwarf galaxies the fit is significantly better with the BEC model. In the second class the rotational velocity of HSB and LSB galaxies exhibits long flat plateaus, resulting in better fit with the NFW model for HSB galaxies and comparable fits for LSB galaxies. We conclude that due to its central density cusp avoidance the BEC model fits better dwarf galaxy dark matter distribution. Nevertheless it suffers from sharp cutoff in larger galaxies, where the NFW model performs better. The investigated galaxy sample obeys the Tully-Fisher relation, including the particular characteristics exhibited by dwarf galaxies. In both models the fitting enforces a relation between dark matter parameters: the characteristic density and the corresponding characteristic distance scale with an inverse power.
European UnionState of HungaryEuropean Social FundJapan Society for the Promotion of Science1. Introduction
The visible part of most galaxies is embedded in a dark matter (DM) halo of yet unknown composition, observable only through its gravitational interaction with the baryonic matter. Assuming the standard ΛCDM cosmological model, the Planck satellite measurements of the cosmic microwave background anisotropy power spectrum support 4.9% baryonic matter, 26.8% DM, and 68.3% dark energy in the Universe [1, 2].
Investigation of mass distribution of spiral galaxies is an essential tool in the research of DM. Beside the stellar disk and central bulge, most of the galaxies harbour a spherically symmetric, massive DM halo, which dominates the dynamics in the outer regions of the stellar disk. Nevertheless there are examples of galaxies which at larger radii are better described by a flattened baryonic mass distribution (global disk model) [3].
Several DM candidates and alternatives have been proposed, the latter assuming Einstein’s theory of gravity breaking down on the galactic scale and above ([4–12]). In brane-world and f(R)-gravity models, the galactic rotation curves could be explained without DM ([13–16]).
At this moment there are strong experimental constraints for all proposed dark matter candidates. Supersymmetric dark matter has been strongly constrained by LHC [17, 18], sterile neutrinos disruled with 99% confidence level by IceCube [19], Weekly Interacting Massive Particles (WIMPs) severely bounded by the LUX [20], PandaX-II [21], and Xenon100 [22] experiments. Extra dimensional effects as dark matter substitutes have been also contained by LHC [23]. Massive Compact Halo Objects (MACHOs) with masses below 20 solar masses have been shown to give at most 10% of dark matter by microlensing experiments on the Large Magellanic Cloud [24]. There is still hope for larger mass MACHOs as dark matter candidates, revived after the spectacular first direct detection of gravitational waves [25], sourced by black holes of approximately 30 solar masses.
It is well known that hot dark matter (HDM) consisting of light (m∝ eV) particles cannot reproduce the cosmological structure formation, as they imply that the superclusters of galaxies are the first structures to form contradicting CMB observations, according to which superclusters would form at the present epoch [26]. Warm dark matter (m∝ keV) models seem to be compatible with the astronomical observations on galactic and also cosmological scales [27, 28]. Leading candidates for warm dark matter are the right handed neutrinos, which in contrast with their left handed counterparts do not participate in the weak interaction. The decay of these sterile neutrinos produces high amount of X-rays, which can boost the star formation rate leading to an earlier reionization [29]. The existence of sterile neutrinos was however severely constrained by recent IceCube Neutrino Observatory experiments [30]. Cold dark matter (CDM) also shows remarkably good agreement with observations over kpc scales ([31, 32]). Particular CDM candidates, like neutralinos (which are stable and can be produced thermally in the early Universe) and other WIMPs originating in supersymmetric extensions of the Standard Model were severely constrained by recent LHC results, rendering them into the range 200 GeV ≲mn≲500 GeV [33]. In a Higgs-portal DM scenario the Higgs boson acts as the mediator particle between DM and Standard Model particles, and it can decay to a pair of DM particles. Very recent constraints established by the ATLAS Collaboration on DM-nucleon scattering cross-section impose upper limits of approximately 60 GeV for each of the scalar, fermion, and vector DM candidates (see Figure 4 of [34]), within the framework of this scenario. While MACHOs of masses less than 10 solar masses (like white dwarfs, neutron stars, brown dwarfs and unassociated planets, and primordial black holes in the astrophysical mass range) were disruled either by Big Bang Nucleosynthesis constraints or microlensing experiments as dominant DM candidates, primordial black holes with intermediate mass could still be viable candidates [35, 36].
Large N-body simulations (e.g., [39]) performed in the framework of the ΛCDM-model (Λ being the cosmological constant) were compatible with CDM halos with central density cusps [40]. They are modeled by the Navarro-Frenk-White (NFW) DM density profile ρNFW(r)=ρs/(r/rs)1+r/rs2, where rs is a scale radius and ρs is a characteristic density. Some observations support such a steep cuspy density profile [41, 42]; nevertheless certain high-resolution rotation curves instead indicate that the distribution of DM in the centres of DM dominated dwarf and Low Surface Brightness (LSB) galaxies is much shallower, exhibiting a core with nearly constant density [43]. In turn, the baryonic matter distribution may also affect the DM density profile. As shown in [44] a dark matter core within an isolated, initially cuspy dark matter halo may form due to strong supernova feedback. By contrast, adiabatic contraction of baryonic gas tends to produce even cuspier dark matter halos [45].
The surface number-density profiles of satellites decline with the projected distance as a power law with the slope -2÷-1.5, while the line-of-sight velocity dispersion declines gradually [46]. These observations support the NFW model on scales of 50÷500 kpc.
In a cosmological setup various scalar field DM models were also discussed ([47, 48] and references therein). A particular scalar field DM model describes light bosons in a dilute gas. The thermal de Broglie wavelength of the particles is λT∝1/mT, which can be large for light bosons (m<eV) and for low temperature. Below a critical temperature (Tc), the bosons’ wave packets, which are the order of λT overlap, result in correlated particles. Such bosons share the same quantum ground state, behaving as a Bose-Einstein Condensate (BEC), characterized by a single macroscopic wave function. It has been proposed that galactic DM halos could be gigantic BECs [49].
It has been shown that caustics of ring shape appear in rotating BEC models, which have an effect on rotation curves, by causing bumps [50, 51]. Such ring shaped caustics degenerate into the origin in the nonrotating BEC limit, adopted in this paper.
The self-gravitating condensate is described by the Gross-Pitaevskii-Poisson equation system in the mean-field approximation [10, 52–54]. In the Thomas-Fermi approximation, a 2-parameter (mass m and scattering length a) density distribution of the BEC halo is obtained [see (3) below], which is less concentrated towards the centre as compared to the NFW model, relaxing the cuspy halo problem.
In model [55] where a normal dark matter phase with an equation of state P=ρc2σtr2 condensed into a BEC with self-interaction (σtr=0.0017 being the one-dimensional velocity dispersion and c the speed of light), the stability of the BEC halo depends on the particle mass and scattering length. For a given mass the stability occurs for larger scattering length and for given scattering length the stability appears at smaller mass. For the following scattering lengths: a=103 fm, a=10-14 fm, and a=10-55 fm, the mass of the BEC particle arises as m>1 eV, m>2×10-6 eV, and m>4.57×10-20 eV, respectively. Galactic size stable halos can form with m>10-24 eV (Figure 3 in [56]).
A stable BEC halo can form as a result of gravitational collapse [57]. The model has been tested on kpc scales confronting it with galactic rotation curve observations [10]. It was pointed out by [58] that the effects of BEC DM should be seen in the matter power spectrum if the boson mass is in the range 15 meV <m<35 meV and 300 meV <m<700 meV for the scattering lengths a=106 fm and a=1010 fm, respectively. In [59] the authors showed that the observed collisional behaviour of DM in Abell 520 cluster can also be recovered within the framework of the BEC model. All of the mentioned BEC particle masses are consistent with the limit m<1.87 eV imposed from galaxy observations and N-body simulation [60]. A discrepancy was however pointed out between the best fit density profile parameters derived from the strong lensing and the galactic rotational curves data. In conclusion the BEC halo should be denser in lens galaxies than in dwarf spheroidals [61].
In this work we critically examine the BEC model as a possible DM candidate against rotation curve data, pointing out both advantages and disadvantages over the NFW model. Previous studies on the compatibility of the BEC model and galactic rotation curves were promising but relied on a less numerous and less diversified set of galaxies than employed here ([62, 63]). The paper has the following structure. The basic properties of the BEC DM model are reviewed in Section 2. In Section 3 a comparison is made between the theoretical predictions of the BEC model and the observed rotation curve data of three types of galaxies: the High Surface Brightness (HSB), LSB, and dwarf galaxies. The conclusions are presented in Section 4.
2. The Bose-Einstein Condensate Galactic Dark Matter Halo
An ideal, dilute Bose gas at very low temperature forms a Bose-Einstein Condensate in which all particles are in the same ground state. In the thermodynamic limit, the critical temperature for the condensation is Tc=2πℏ2n/ζ3/2/mkB [64]. Here n and m are the number density and the mass of the bosons, respectively, and ζ=2.612 is a constant, while ℏ and kB denote the reduced Planck and Boltzmann constants, respectively. Atoms can be regarded as quantum-mechanical wave packets of the order of their thermal de Broglie wavelength λT=2πℏ2/mkBT. The condition for the condensation T<Tc can be reformulated as l<λT/ζ-1/3, where l is the average distance between pairs of bosons, and it occurs when the temperature, hence the momentum of the bosons, decreases and as a consequence their de Broglie wavelengths overlap. The thermodynamic limit is only approximately realized, the finite size giving corrections to the critical temperature [65–68]. A dilute, nonideal Bose gas also displays BEC; on the other hand, the condensate fraction is smaller than unity at zero temperature and the critical temperature is also modified [69–72]. Experimentally, BEC (which could be formed by bosonic atoms but also form fermionic Cooper pairs) has been realized first in Rb87 [73–75], then in Na23 [76, 77], and in Li7 [78].
In a dilute gas, only two-particle interactions dominate. The repulsive, two-body interparticle potential is approximated as Vself=λδr-r′, with a self-coupling constant λ=4πℏ2a/m, where a is the scattering length. Then in the mean-field approximation (in case when we neglect the contribution of the excited states) the BEC is described by the Gross-Pitaevskii equation [52–54]:(1)iℏ∂∂tψr,t=-ℏ22mΔ+Vselfgravr+λρr,tψr,t,where ψ(r,t) is the wave function of the condensate and Δ is the 3-dimensional Laplacian. The probability density ρr,t=ψ(r,t)2 is normalized to(2)n0t=∫drρr,t,where n0t is the number of particles in the ground state and ρr,t is the number density of the condensate. The potential Vselfgrav(r)/m is the Newtonian gravitational potential produced by the Bose-Einstein Condensate.
Stationary solutions of the Gross-Pitaevskii equation can be found in a simple way by using the Madelung representation of complex wave functions [79, 80], then deriving the Madelung hydrodynamic equations [79]. Madelung’s equations can be interpreted as the continuity and Euler equations of fluid mechanics, with quantum corrections included. However, the quantum correction potential in the generalized Euler equation contributes significantly only close to the boundary of the system [81]. In the Thomas-Fermi approximation the quantum correction potential is neglected compared to the self-interaction term. This approximation becomes more accurate as the particle number increases [82].
Assuming a spherically symmetric distribution of the condensate the following solution was found [10, 81]:(3)ρBECr=ρBECcsinkrkr,where ρBEC=mρr and(4)k=Gm3ℏ2a.The central density ρBECc≡ρBEC0 is determined from normalization condition (2) as(5)ρBECc=n0mk34π2. The Thomas-Fermi approximation remains valid for n0≫1/ka [81].
The BEC galactic DM halo’s size is defined by ρ(RBEC)=0, giving k=π/RBEC, that is,(6)RBEC=πℏ2aGm3. The mass profile of the BEC halo is then given as(7)mBECr=4π∫0rρBECrr2dr=4πρBECck2rsinkrkr-coskr. The BEC halo contributes to the velocity profile of the particles which are moving on circular orbit as dictated by the Newtonian gravitational force [10]. This can be taken into account by the following equation:(8)v2r=4πGρBECck2sinkrkr-coskr,which needs to be added to the baryonic contribution, respectively.
3. Confronting the Model with Rotation Curve Data
The validity of our model was tested by confronting the rotation curve data of a sample of 6 HSB, 6 LSB, and 7 dwarf galaxies, with both the NFW DM and the BEC density profiles. For reasons to become obvious during our analysis, we split both the HSB and LSB data sets into two groups (type I. and II.), based on the shapes of the curves. In the first group the rotational velocities increase over the whole observed range, while in the second set the rotation curves exhibit long flat regions.
The commonly used NFW model has the mass density profile(9)ρNFWr=ρsr/rs1+r/rs2, where ρs and rs are a characteristic density and distance scale, to be determined from the fit.
The mass within a sphere with radius r=yrs is then given as(10)mNFWr=4πρsrs3ln1+y-y1+y, where y is a positive dimensionless radial coordinate.
3.1. HSB Galaxies
In this subsection we will follow the method described in [15]. In HSB galaxy the baryonic component was decomposed into a thin stellar disk and a spherically symmetric bulge. It was assumed that the mass distribution of bulge component follows the deprojected luminosity distribution with a factor known as the mass-to-light ratio. The bulge parameters were estimated from a Sérsic r1/n bulge model, which was obtained by the fitting of the optical I-band galaxy light profiles.
Each galaxy’s spheroidal bulge component has a surface brightness profile which is described by a generalized Sérsic function [83];(11)Ibr=I0,bexp-rr01/n,wherein I0,b is the central surface brightness of the bulge, r0 is the characteristic radius of the bulge, and the magnitude-radius curve’s shape parameter is denoted by n.
The mass-to-light ratio for the Sun is γ⊙=5133 kg W-1. The mass-to-light ratio of the bulge σ will be given in units of γ⊙ (solar units). We will also give the mass in units of the solar mass M⊙=1.98892×1030 kg. We assume that the radial distribution of visible mass follows the radial distribution of light derived from the bulge-disk decomposition. Accordingly the mass of the bulge inside the projected radius r can be derived from the surface brightness observed within this radius:(12)mbr=σNDF⊙2π∫0rIbrrdr, where F⊙D is the apparent flux density of the Sun at a distance D Mpc, F⊙D=2.635×106-0.4f⊙ mJy, with f⊙= 4.08 + 5lg(D/1Mpc) + 25 mag, and(13)ND=4.4684×10-35D-2m-2arcsec2. The rotational velocity related to the bulge is(14)vb2r=Gmbrr, where G is the gravitational constant.
In case of a spiral galaxy, the radial surface brightness profile of the disk decreases exponentially as a function of the radius [84];(15)Idr=I0,dHSBexp-rhHSB,where I0,dHSB is the central surface brightness of the disk and hHSB is a characteristic disk length scale. The disk contributes to the circular velocity as follows ([84]):(16)vd2x=GMDHSB2hHSBx2I0K0-I1K1, where x=r/hHSB and In and Kn are the modified Bessel functions evaluated at x/2, while MDHSB is the total mass of the disk.
Accordingly in HSB galaxy the rotational velocity adds up as(17)vtg2x=vb2x+vd2x+vDM2x.
In order to validate the BEC+baryonic model, we confront it with rotation curve data of 6 well-tested galaxies (which were already employed in [15] for testing a brane-world model). The data was obtained from a sample given in [37] and meets the following criteria: (i) it has to be among the best accuracies obtained from the sample and (ii) the bulge has to be spherically symmetric. As a check we also fitted the NFW + baryonic model with the same data set. The respective rotation curves are plotted for both models on Figures 1 and 2. The small humps on both figures are due to the baryonic component. From the available photometric data the best fitting values were derived for the baryonic model parameters I0,b, n, r0, rb, and I0,dHSBhHSB. By fitting BEC and NFW models to the investigated rotation curve data, the parameters for these models (as well as the corresponding baryonic parameters) were calculated. The parameter values are indicated in Tables 1 and 2.
The distances (D) and the photometric parameters of the 6 HSB galaxy sample as determined by the fit with available photometric data [37]. Bulge parameters: the central surface brightness (I0,b), the shape parameter (n), the characteristic radius (r0), and radius of the bulge (rb). Disk parameters: central surface brightness (I0,dHSB) and length scale (hHSB) of the disk.
Galaxy
D
I0,b
n
r0
rb
I0,dHSB
hHSB
Mpc
mJy/arcsec2
kpc
kpc
mJy/arcsec2
kpc
ESO215G39
61.29
0.1171
0.6609
0.78
2.58
0.0339
4.11
ESO322G77
38.19
0.1949
0.7552
0.33
1.37
0.0744
2.20
ESO509G80
92.86
0.2090
0.7621
1.10
4.69
0.0176
11.03
ESO323G25
59.76
0.1113
0.4626
0.43
0.99
0.0825
3.47
ESO383G02
85.40
0.6479
0.7408
0.42
1.94
0.5118
3.82
ESO446G01
98.34
0.2093
0.8427
1.28
6.33
0.0357
5.25
The best fit parameters and the minimum values (χmin2) of the χ2 statistics for the HSB I and II galaxies (the first and last three galaxies, respectively). Columns 2–5 give the BEC model parameters (radius RBEC and central density ρBEC(c) of the BEC halo) and the corresponding baryonic parameters (mass-to-light ratio σBEC of the bulge and total mass of the disk MDHSBBEC). Columns 7–10 give the NFW model parameters (scale radius rs and characteristic density ρs of the halo) and the corresponding baryonic parameters (mass-to-light ratio σNFW of the bulge and total mass of the disk MDHSBNFW). The 1σ confidence levels are shown in the last column (these are the same for both models). For HSB I galaxies the two models give similar χmin2 values (within 1σ confidence level); however in case of HSB II galaxies with extended flat regions, the NFW model fits better the rotation curves. The χmin2 values in the case of BEC model are outside the 1σ confidence level for HSB II galaxies.
Galaxy
σ (BEC)
MDHSB (BEC)
RBEC
ρBEC(c)
χmin2 (BEC)
σ (NFW)
MDHSB (NFW)
rs
ρs
χmin2 (NWF)
1σ
⊙
1010M⊙
kpc
10-21 kg/m3
⊙
1010M⊙
kpc
10-24 kg/m3
ESO215G39
0.3
5.61
3.8
2.0
23.07
0.6
3.84
187
14.7
22.22
34.18
ESO322G77
1.6
5.1
0.8
89.0
9.15
2.5
3.79
709
8
7.69
11.53
ESO509G80
1.4
48.74
9.7
1.2
12.52
0.9
11
22
800
33.48
36.3
ESO323G25
2.5
12.18
2.5
11.8
222.74
6
9.43
436
6
80.55
66.74
ESO383G02
0.13
8.77
3.0
5.7
48.83
1.7
6.32
459
4.2
23.3
47.9
ESO446G01
0.6
12.77
4.6
5.9
86.02
1.4
6.7
786
4.1
43.37
44.74
Best fit curves for the HSB I. galaxy sample where the solid black lines hold for the baryonic matter + BEC model, while the dashed red lines refer to the baryonic matter + NFW model.
Best fit curves for the HSB II. galaxy sample. The solid black lines hold for the baryonic matter + BEC model, while the dashed red lines hold for the baryonic matter + NFW model. The BEC model does not describe well the extended flat regions.
Both the BEC and NFW DM models give comparable χmin2 values (within 1σ confidence level) for HSB I galaxies. In case of galaxies with extended flat regions (HSB II), the NFW DM model fits better the rotation curves; nevertheless BEC model gives rotational curves which fall outside the 1σ confidence level.
3.2. LSB Galaxies
The surface brightness of LSB galaxies is substantially fainter than the brightness of the sky at night. They belong to an early stage class of galaxies [85]. LSB galaxies were found to be metal poor, which indicates a lower star formation rate than what is generally found in HSB galaxies [86]. Wide spectrum of colors can be measured in case of LSB galaxies ranging from red to blue [87] and they are diverse as regards morphologies and other properties. Most of the LSB galaxies that were observed are dwarf galaxies; however there is also a significant number of large spirals among LSB galaxies [88].
According to our model the LSB galaxy is made up of two main components, one being a thin stellar + gas disk and the other one being a CDM component which is assumed to be a BEC. We use the same model for the disk component as in the case of the HSB galaxies. The surface brightness profile can be described by the following equation [84]:(18)Idr=I0,dLSBexp-rhLSB, where I0,dLSB is the central surface brightness and hLSB is the disk length scale. The contribution of the disk to the circular velocity can be expressed as (19)vd2r=GMDLSB2hLSBq2I0K0-I1K1,where q=r/hLSB and MDLSB are the total mass of the disk while the modified Bessel functions In and Kn are evaluated at q/2.
Consequently, for an arbitrary projected radius r the rotational velocity can be calculated based on the combined model resulting in the following equation:(20)vtg2r=vd2r+vDM2.
A preliminary check confirmed that the BEC + baryonic model represents a better fit than the purely BEC model.
We confronted the BEC model with 6 LSB galaxies chosen from a larger sample [38]. The applied data were obtained from both HI and Hα measurements. From a χ2-test the parameters in both the BEC + baryonic and NFW + baryonic models were identified: these are shown in Table 3. The best fit rotation curves are represented on Figures 3 and 4.
The best fit BEC and NFW parameters of the LSB I and II type galaxies (the first and last three galaxies, respectively). D is taken from [38]. The rest of the parameters are rotation curve fits. For LSB I galaxies the BEC DM model gives significantly better fitting velocity curves (within 1σ confidence level) than the NFW model. However the velocity curves are outside the 1σ confidence level for LSB II galaxies.
Galaxy
D
hLSB
MDLSB (BEC)
ρBEC(c)
RBEC
χmin2 (BEC)
MDLSB (NFW)
ρs
rs
χmin2 (NFW)
1σ
Mpc
kpc
109M⊙
10-21 kg/m3
kpc
109M⊙
10-24 kg/m3
kpc
DDO 189
12.6
1.9
2.71
0.38
8.3
0.519
2.16
16
70
1.09
7.03
NGC 4455
6.8
2.3
0.231
1.44
5.5
9.29
0.11
25.04
66
5.39
18.11
UGC 10310
15.6
5.2
0.443
0.98
7.8
2.66
0.9
14.9
88
5.76
13.74
NGC 2366
3.4
1.5
2.43
0.22
5.3
110.73
2.5
0.2
1000
116.93
26.72
NGC 5023
4.8
0.8
0.894
2.45
5.6
53.2
0.0449
457
13
143.08
32.05
NGC 3274
6.7
0.5
1.1
1.69
6.4
269.8
0.252
2373
4
148.44
20.27
Best fit curves for the LSB I. galaxy sample. The solid black lines indicate the baryonic matter + BEC model, while the dashed red lines indicate the baryonic matter + NFW model.
Best fit curves for the LSB II. galaxy sample. The solid black lines refer to the baryonic matter + BEC model, while the dashed red lines to the baryonic matter + NFW model. As for HSB galaxies, the BEC model fails to explain the extended flat regions of the rotation curves.
For the LSB I galaxies the BEC DM model gives significantly better fitting velocity curves (all within the 1σ confidence level) compared to the NFW model (which in two cases out of the three gives fits falling outside 1σ). For LSB II galaxies the quality of the fits is comparable, but in both models they are beyond the 1σ confidence level.
3.3. Dwarf Galaxies
Approximately 85% of the explored galaxies in the Local Volume [89] are dwarf galaxies. The dwarfs are defined by having an absolute magnitude which is fainter than MB~-16 mag. On the other hand they are larger than globular clusters [90].
Although little is known about their formation, it is generally accepted that dwarfs are formed at the centres of subhalos. Dwarf galaxies can be categorised in five groups according to their optical appearance. The five groups being dwarf ellipticals, dwarf irregulars, dwarf spheroidals, blue compact dwarfs, and dwarf spirals. The dwarfs falling in the last group represent the very small ends of spirals [91]. Dwarf spheroidals are old systems and among the most DM dominated galaxies in the Universe.
The central velocity dispersion of most dwarf galaxies is in the range 6 ÷ 25 km/s [92]. In a typical dwarf galaxy, assuming dynamical equilibrium, the mass derived from the observed velocity dispersion is substantially greater than the observed total visible mass. This implies that the mass-to-light ratio is very high compared to other types of galaxies; hence they can greatly contribute to the understanding of DM distribution on small scales. Dwarf galaxies allow for proving or falsifying different alternative gravity theories [93].
We decided to use 7 dwarf galaxies for testing the BEC model. We have selected the sample dwarf galaxies such as to ensure that sufficient high-resolution rotation curve data would be available for our study. We fitted both the BEC + baryonic and the NFW + baryonic models, respectively, with similar baryonic components as for the LSB galaxies. As the length scales of the stellar disks were not available for the selected sample, they were calculated by χ2 minimization, too.
A preliminary check showed that the addition of the BEC dark matter halo to the baryonic model improved (giving lower χmin2 values) on the fit in all cases. By contrast, the NFW model was unable to improve on the purely baryonic fit in four out of seven cases. We note that since the data does not contain the error margins, the χmin2 values are relatively high (beyond the 1σ confidence level in most cases). The best fit BEC and NFW parameters are shown in Table 4 and the corresponding rotation curves are represented on Figure 5. The inclusion of the BEC DM model gives significantly (in some cases one order of magnitude in the value of χ2) better fits compared to the case of NFW model. This is due to the cusp avoidance in the central density profile of the BEC model and the fact that dwarf galaxies do not exhibit extended flat regions in their rotation curves.
The best rotation curve data fit BEC + baryonic and NFW + baryonic parameters for the dwarf galaxy sample.
Galaxy
hdwarf (BEC)
MDdwarf (BEC)
ρBEC(c)
RBEC
χmin2 (BEC)
hdwarf (NFW)
MDdwarf (NFW)
ρs
rs
χmin2 (NFW)
1σ
kpc
109M⊙
10-21 kg/m3
kpc
kpc
109M⊙
10-24 kg/m3
kpc
IC 2574
1.2
0.1122
0.4
13
68.47
7.9
28.44
0
0
714.73
44.74
HoI
0.2
0.0107
3.6
1.9
95.26
0.9
0.533
0
0
241.30
20.27
HoII
1.2
0.4431
0.2
7.69
33.33
1.7
0.642
4
92
43.86
26.72
DDO 39
1.3
1.1235
0.7
10.01
69.39
4.3
7.21
43
35
69.82
17.02
DDO 53
0.2
0.0061
1.8
2.5
20.05
1.6
0.976
1
24
51.53
10.42
DDO 154
3.1
3.3502
0.2
5.8
1.48
3.2
4.52
0
0
9.40
9.30
M81dwB nor
0.9
1.023
3.7
0.7
6.19
0.7
0.705
0
0
8.4
10.42
The best fit curves for the dwarf galaxy sample. The BEC + baryonic model (solid black curves) give a better fit in all cases than the NFW + baryonic model (dashed red lines). In both cases the fit was performed with the same baryonic model.
4. Discussions and Final Remarks
We have performed a χ2-test of the BEC and NFW DM models, with the rotation curves of 6 HSB, 6 LSB, and 7 dwarf galaxy samples. For improved accuracy we also included realistic baryonic models in every case. For the HSB galaxy sample, both the rotation curve and the surface photometry data were available. Most of the rotation curves were smooth, symmetric, and uniform in quality.
The circular velocity of the investigated galaxies was decomposed into its baryonic and DM contribution: vmodel2(r)=vbaryonic2+vDM2. For the BEC model the DM contribution to the rotational velocity can be described as (8). Then the rotation curves are fitted with the parameters of the baryonic and DM halo models (BEC and NFW) using χ2 minimization method.
The analysis of the HSB I galaxies showed a remarkably good agreement for both DM models with observations. The BEC and NFW models show similar fits. However, the rotation curves of the HSB II type galaxies are significantly better described by the NFW model.
It was previously known that for LSB galaxies and without including the baryonic sector, the BEC model gave a better fit than the NFW model [62]. We additionally found that including the baryonic component improves on the fit of [62]. Our detailed analysis showed a significantly better performance of the BEC model for LSB type I galaxies, while comparable fits for LSB type II galaxies were obtained. These latter fits were however outside the 2σ confidence level.
The unsatisfactory large distance behaviour of the BEC model for both the HSB and LSB galaxies of type II originates in the sharp cutoff of the BEC DM distribution and clearly indicates that it would be desirable to modify the BEC model on larger scale, also to comply with the behaviour of the universal rotation curves (URCs) at larger radii [94].
From the above analysis of HSB and LSB galaxies it is also obvious that (while on large distances the BEC model suffers from problems due to the sharp cutoff) close to the core it works overall better than the NFW model. This is also supported by our fit of both the BEC + baryonic and NFW + baryonic DM models with rotation curve data of a sample of 7 dwarf galaxies. Since dwarf galaxies are DM dominated, they allow for the best comparison between the various models. The results can be seen in Figure 5. We also note that the NFW DM improved over the pure baryonic fit in four cases out of seven, while including the BEC component improved over the fit with the baryonic component in all cases.
The BEC parameters were determined for all cases. The parameters ρBEC(c), RBEC are given in Tables 2, 3, and 4. The averages of the radii RDM of the BEC halos for the HSB, LSB, and dwarf galaxies are 〈RBECHSB〉≈4.06 kpc, 〈RDMLSB〉≈6.48 kpc, and 〈RDMdwarf〉≈5.94 kpc, respectively. The scatter however is large; there are no universal BEC parameters which globally fit all the galaxies, not even at 3σ confidence level. The closer to this goal were the HSB galaxies, where 3 out of 6 had overlapping 3σ domains. Nonetheless the given values of RDM are consistent within the order of magnitude with the halo radii of 59 other galaxies determined from weak lensing [95].
We represent the density parameter ρBEC(c) of the BEC model as function of RBEC in Figure 6(a) and the density parameter ρs of the NFW model as function of rs in Figure 6(b) (four dwarf galaxies are absent, as the NFW halo does not improve the fit over the pure baryonic case). The fitting enforces a relation between the dark matter parameters: the characteristic density scales with an inverse power with the corresponding characteristic distance.
The density parameter ρBEC(c) of the BEC model is shown as function of RBEC in (a), and the density parameter ρs of the NFW model is shown as function of rs in (b). The HSB, LSB, and dwarf galaxies are represented by filled circles, empty triangles, and filled triangles, respectively.
We verify the Tully-Fisher relation for the investigated galaxy sample and present the results on Figure 7. Apparent B magnitudes and galaxy distances were collected from the NASA/IPAC extragalactic database [96] and were corrected for extinction based on Landolt standard fields to calculate the absolute magnitudes. It is known that the Tully-Fisher relation holds for spiral and lenticular galaxies with the same slope (e.g., [97]). A larger slope and scatter characterize the Tully-Fisher relation for the dwarf galaxies (e.g., [97, 98]). The investigated sample exactly exhibits these features.
The baryonic Tully-Fisher relation of our galaxy sample. Absolute B magnitudes are presented as function of the logarithm of the maximal rotational velocity. The HSB, LSB, and dwarf galaxies are represented by filled circles, empty triangles, and filled triangles, respectively.
There is a relation among the mass m of the BEC particle, its coherent scattering length a, and the radius of the DM halo RDM [10]:(21)m=π2ℏ2aGRBEC21/3≈6.73×10-2afm1/3RBECkpc-2/3eV. Axions have been proposed as the Peccei-Quinn solution to the strong CP problem [99] and they are among the best dark matter candidates. Being bosons, they may also form BEC. The Axion Dark Matter Experiment has already established limits on the dark matter axions [100, 101].
Assuming the BEC is formed of axions with mass of 10-6 eV, the scattering lengths for the three types of galaxies emerge as aHSB≈5.4×10-14 fm, aLSB≈1.37×10-13 fm, and adwarf≈1.15×10-13 fm. These values are consistent with the results of [95], which are based on a statistical analysis of 61 DM dominated galaxies. The total energy of the BEC halo is negative with these scattering lengths and particle mass, meaning the halo is stable (see Figure 3 of [56]).
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
In the earlier stages of this work M. Dwornik, Z. Keresztes, and L. Á. Gergely were supported by the European Union and the State of Hungary, cofinanced by the European Social Fund in the framework of TÁMOP 4.2.4. A/2-11-/1-2012-0001 “National Excellence Program”. L. Á. Gergely was also supported by the Japan Society for the Promotion of Science.
AdeP. A. R.AghanimN.Armitage-CaplanC.Planck 2013 results. I. Overview of products and scientific results2014571, article 1FrancisM.First Planck results: the Universe is still weird and interesting2013Ja lochaJ.BratekL.KutscheraM.SkindzierP.Global disc models for galaxies NGC 1365, 6946, 7793 and UGC 6446201040642805281610.1111/j.1365-2966.2010.16887.xMilgromM.A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis198327036537010.1086/161130SandersR. H.Anti-gravity and galaxy rotation curves19841362L21L23MoffatJ.SokolovI.Galaxy dynamics predictions in the nonsymmetric gravitational theory19963781–4596710.1016/0370-2693(96)00366-8MannheimP. D.Are galactic rotation curves really flat?1997479265966410.1086/3039332-s2.0-33645696779RobertsM. D.Galactic metrics200436112423243110.1023/B:GERG.0000046830.85831.4aMR2102145ZBL1064.85001BöhmerC. G.HarkoT.On Einstein clusters as galactic dark matter haloes2007379139339810.1111/j.1365-2966.2007.11977.x2-s2.0-34547175149BöhmerC. G.HarkoT.Can dark matter be a Bose-Einstein condensate?20072007610.1088/1475-7516/2007/06/0252-s2.0-34548623343BertolamiO.B{\"o}hmerC. G.HarkoT.LoboF. S.Extra force in f(R) modified theories of gravity2007751010401610.1103/physrevd.75.104016MR2326814BöhmerC. G.HarkoT.LoboF. S. N.Dark matter as a geometric effect in f (R) gravity200829638639210.1016/j.astropartphys.2008.04.0032-s2.0-46749144491MakM. K.HarkoT.Can the galactic rotation curves be explained in brane world models?200470202401010.1103/physrevd.70.0240102-s2.0-19244365405RahamanF.KalamM.DeBenedictisA.UsmaniA. A.SaibalR.Galactic rotation curves and brane world models20083891273310.1111/j.1365-2966.2008.13559.xGergelyL. Á.HarkoT.DwornikM.KupiG.KeresztesZ.Galactic rotation curves in brane world models201141543275329010.1111/j.1365-2966.2011.18941.x2-s2.0-80051810414StabileA.CapozzielloS.Galaxy rotation curves in f(R,ϕ) gravity201387606400210.1103/physrevd.87.064002ConoverE.Physicists narrow in on electrical short in Large Hadron Collider201510.1126/science.aab0390Arkani-HamedN.2016Copenhagen, DenmarkCurrent Themes in High Energy Physics and CosmologyAartsenM. G.AbrahamK.AckermannM.Searches for sterile neutrinos with the IceCube detector20161177907180110.1103/PhysRevLett.117.071801AkeribD. S.AraujoH. M.BaiX.First results from the LUX dark matter experiment at the sanford underground research facility2013112909130310.1103/PhysRevLett.112.091303PandaX-II CollaborationDark matter results from First 98.7-day data of pandaX-II experiment201611712130310.1103/PhysRevLett.117.121303AprileE.AalbersJ.AgostiniF.XENON100 dark matter results from a combination of 477 live days2016941212200110.1103/PhysRevD.94.122001ChoudhuryD.GhoshK.Bounds on universal extra dimension from LHC run I and II data201676315516010.1016/j.physletb.2016.10.010MR3577296AlcockC.AllsmanR. A.AlvesD. R.AxelrodT. S.BeckerA. C.BennettD. P.CookK. H.DalalN.DrakeA. J.FreemanK. C.GehaM.GriestK.LehnerM. J.MarshallS. L.MinnitiD.NelsonC. A.PetersonB. A.PopowskiP.PrattM. R.QuinnP. J.StubbsC. W.SutherlandW.TomaneyA. B.VandeheiT.WelchD.The MACHO project: microlensing results from 5.7 years of large magellanic cloud observations2000542128130710.1086/3095122-s2.0-0034633466AbbottB. P.AbbottR.AbbottT. D.Observation of gravitational waves from a binary black hole merger2016116606110210.1103/PhysRevLett.116.061102CaldwellD. O.2001Berlin, HeidelbergSpringer10.1007/978-3-662-04597-8de VegaH. J.SanchezN. G.Warm dark matter in the galaxies: theoretical and observational progresses. Highlights and conclusions of the chalonge meudon workshop 2011https://arxiv.org/abs/1109.3187WeiH.ChenZ.-C.LiuJ.Cosmological constraints on variable warm dark matter20137204-527127610.1016/j.physletb.2013.02.0272-s2.0-84875553653BiermannP. L.KusenkoA.Relic keV sterile neutrinos and reionization200696909130110.1103/PhysRevLett.96.0913012-s2.0-33644872808AartsenM. G.AckermannM.AdamsJ.Searches for sterile neutrinos with the icecube detector2016117707180110.1103/physrevlett.117.071801PadmanabhanT.Cosmological constant—the weight of the vacuum20033805-623532010.1016/s0370-1573(03)00120-0MR19926802-s2.0-0038457945PeeblesP. J.RatraB.The cosmological constant and dark energy200375255960610.1103/RevModPhys.75.559MR1973358ZBL1205.830822-s2.0-0038579350FowlieA.KowalskaK.RoszkowskiL.SessoloE. M.TsaiY.-L. S.Dark matter and collider signatures of the MSSM201388505501210.1103/physrevd.88.0550122-s2.0-84884872217AadG.AbajyanT.AbbottB.Search for invisible decays of a higgs boson produced in association with a Z boson in ATLAS2014112201920180210.1103/PhysRevLett.112.201802FramptonP. H.Angular momentum of dark matter black holeshttps://arxiv.org/abs/1608.05009SasakiM.SuyamaT.TanakaT.YokoyamaS.Primordial black hole scenario for the gravitational-wave event GW1509142016117606110110.1103/physrevlett.117.0611012-s2.0-84982085327PalunasP.WilliamsT. B.Maximum disk mass models for spiral galaxies200012062884290310.1086/3168782-s2.0-0000562520de BlokW. J. G.BosmaA.High-resolution rotation curves of low surface brightness galaxies2002385381684610.1051/0004-6361:200200802-s2.0-0036530781SpringelV.WhiteS. D. M.JenkinsA.FrenkC. S.YoshidaN.GaoL.NavarroJ.ThackerR.CrotonD.HellyJ.PeacockJ. A.ColeS.ThomasP.CouchmanH.EvrardA.ColbergJ.PearceF.Simulations of the formation, evolution and clustering of galaxies and quasars2005435704262963610.1038/nature035972-s2.0-20444461186NavarroJ. F.FrenkC. S.WhiteS. D. M.The structure of cold dark matter halos1996462256357510.1086/1771732-s2.0-25044479731ValenzuelaO.RheeG.KlypinA.GovernatoF.StinsonG.QuinnT.WadsleyJ.Is there evidence for flat cores in the halos of dwarf galaxies? The case of NGC 3109 and NGC 682220076572 I77378910.1086/5086742-s2.0-33947612980JardelJ. R.GebhardtK.FabriciusM. H.DroryN.WilliamsM. J.Measuring dark matter profiles non-parametrically in dwarf spheroidals: an application to Draco20137632, article 9110.1088/0004-637x/763/2/912-s2.0-84872709611BurkertA.Klapdor-KleingrothausH. V.RamachersY.The structure of dark matter halos. Observation versus theory1997SingaporeWorld Scientific35TeyssierR.PontzenA.DuboisY.ReadJ. I.Cusp-core transformations in dwarf galaxies: observational predictions201342943068307810.1093/mnras/sts5632-s2.0-84874194664InoueS.SaitohT. R.Cores and revived cusps of dark matter haloes in disc galaxy formation through clump clusters201141842527253110.1111/j.1365-2966.2011.19873.x2-s2.0-83155181328KlypinA.PradaF.Testing gravity with motion of satellites around galaxies: Newtonian gravity against modified Newtonian dynamics200969021488149610.1088/0004-637x/690/2/14882-s2.0-60349088738GergelyL. Á.TsujikawaS.Effective field theory of modified gravity with two scalar fields: dark energy and dark matter201489606405910.1103/physrevd.89.064059Rodríguez-MontoyaI.MagañaJ.MatosT.Pérez-LorenzanaA.Ultra light bosonic dark matter and cosmic microwave background201072121509151410.1088/0004-637X/721/2/15092-s2.0-78149318327SinS.-J.Late-time phase transition and the galactic halo as a Bose liquid1994506, article 365010.1103/physrevd.50.36502-s2.0-0002577434SikivieP.Caustic rings of dark matter19984321-21391442-s2.0-0342291678SikivieP.Caustic ring singularity19996061162-s2.0-18044395933GrossE. P.Structure of a quantized vortex in boson systems1961203454477MR0128907GrossE. P.Hydrodynamics of a superfluid condensate19634219510.1063/1.1703944PitaevskiiL. P.Vortex lines in an imperfect bose gas1961402646HarkoT.Cosmological dynamics of dark matter Bose-Einstein condensation2011831212351510.1103/physrevd.83.1235152-s2.0-79960760918De SouzaJ. C. C.PiresM. O. C.Discussion on the energy content of the galactic dark matter Bose-Einstein condensate halo in the Thomas-Fermi approximation201420143, article no. 01010.1088/1475-7516/2014/03/0102-s2.0-84897491240HarkoT.Gravitational collapse of Bose-Einstein condensate dark matter halos201489808404010.1103/physrevd.89.0840402-s2.0-84898656239VeltenH.WambaE.Power spectrum for the Bose-Einstein condensate dark matter20127091-21510.1016/j.physletb.2012.01.0712-s2.0-84857059538LeeJ. W.LimS.ChoiD.BEC dark matter can explain collisions of galaxy clustershttps://arxiv.org/abs/0805.3827BoyanovskyD.De VegaH. J.SanchezN. G.Constraints on dark matter particles from theory, galaxy observations, and N-body simulations200877404351810.1103/physrevd.77.0435182-s2.0-41049109844González-MoralesA. X.Diez-TejedorA.Ureña-LópezL. A.ValenzuelaO.Hints on halo evolution in scalar field dark matter models with galaxy observations201387210.1103/physrevd.87.021301RoblesV. H.MatosT.Flat central density profile and constant dark matter surface density in galaxies from scalar field dark matter2012422128228910.1111/j.1365-2966.2012.20603.x2-s2.0-84859863274DwornikM.KeresztesZ.GergelyL. Á.KinjoN.NakajimaA.Rotation curves in bose-einstein condensate dark matter halos2014Nova Science195219PitaevskiiL. P.StringariS.2003New York, NY, USAOxford University PressGrossmannS.HolthausM.On Bose-Einstein condensation in harmonic traps1995208318819210.1016/0375-9601(95)00766-V2-s2.0-0006712406KetterleW.van DrutenN. J.Two-step condensation of the ideal bose gas in highly anisotropic traps199654656KirstenK.TomsD. J.Bose-Einstein condensation of atomic gases in a general harmonic-oscillator confining potential trap1996545, article 418810.1103/physreva.54.41882-s2.0-0000680986HaugerudH.HaugsetT.RavndalF.A more accurate analysis of Bose-Einstein condensation in harmonic traps19972251-318222-s2.0-0039384215GiorginiS.PitaevskiiL. P.StringariS.Theory of ultracold atomic Fermi gases200880410.1103/RevModPhys.80.12152-s2.0-54249144336GlaumK.PelsterA.KleinertH.PfauT.Critical temperature of weakly interacting dipolar condensates200798808040710.1103/PhysRevLett.98.0804072-s2.0-33847309628SchütteM.PelsterA.JankeW.PelsterA.Critical temperature of a bose-einstein condensate with 1/r interactionsProceedings of the 9th International Conference on Path Integrals: New Trends and Perspectives (PI '07)September 2007Dresden, GermanyWorld Scientific4174202-s2.0-84893126762DalfovoF.GiorginiS.PitaevskiiL. P.StringariS.Theory of Bose-Einstein condensation in trapped gases199971346351210.1103/RevModPhys.71.4632-s2.0-0033246313AndersonM. H.EnsherJ. R.MatthewsM. R.WiemanC. E.CornellE. A.Observation of Bose-Einstein condensation in a dilute atomic vapor1995269522119820110.1126/science.269.5221.1982-s2.0-11944274056HanD. J.WynarR. H.CourteilleP.HeinzenD. J.Bose-Einstein condensation of large numbers of atoms in a magnetic time-averaged orbiting potential trap1998576R4114R411710.1103/physreva.57.r41142-s2.0-0000661053ErnstU.MarteA.SchreckF.SchusterJ.RempeG.Bose-Einstein condensation in a pure Ioffe-Pritchard field configuration19984111610.1209/epl/i1998-00107-22-s2.0-0031629521DavisK. B.MewesM.-O.AndrewsM. R.Van DrutenN. J.DurfeeD. S.KurnD. M.KetterleW.Bose-Einstein condensation in a gas of sodium atoms19957522, article 396910.1103/physrevlett.75.39692-s2.0-4244115335HauL. V.BuschB. D.LiuC.DuttonZ.BurnsM. M.GolovchenkoJ. A.Near-resonant spatial images of confined Bose-Einstein condensates in a 4-Dee magnetic bottle1998581R54R5710.1103/physreva.58.r542-s2.0-1642337912BradleyC. C.SackettC. A.TollettJ. J.HuletR. G.Evidence of Bose-Einstein condensation in an atomic gas with attractive interactions19957591687169010.1103/PhysRevLett.75.16872-s2.0-4243132347MadelungE.Quantum theory in hydrodynamic form192640322326SonegoS.Interpretation of the hydrodynamical formalism of quantum mechanics199121101135118110.1007/BF00734264MR11381062-s2.0-0001055955WangX. Z.Cold Bose stars: self-gravitating Bose-Einstein condensates200164121240092-s2.0-33750542370LiebE. H.SeiringerR.YngvasonJ.A rigorous derivation of the Gross-Pitaevskii energy functional20006140436022-s2.0-0347107495SérsicJ. L.1968Cordoba, ArgentinaObservatorio AstronomicoFreemanK. C.On the disks of spiral and so galaxies197016081110.1086/150474ImpeyC.BothunG.Low surface brightness galaxies199735126730710.1146/annurev.astro.35.1.2672-s2.0-0031440504McgaughS. S.Oxygen abundances in low surface brightness disk galaxies1994426113514910.1086/1740492-s2.0-12044256936O'NeilK.BothunG. D.SchombertJ.CornellM. E.ImpeyC. D.A wide field ccd survey for low surface brightness galaxies. II. Color distributions, stellar populations, and missing baryons19971146244810.1086/1186592-s2.0-0037893793BeijersbergenM.De BlokW. J. G.Van Der HulstJ. M.Surface photometry of bulge dominated low surface brightness galaxies199935139039192-s2.0-0039826218KarachentsevI. D.KarachentsevaV. E.HuchtmeierW. K.MakarovD. I.A catalog of neighboring galaxies200412742031206810.1086/3829052-s2.0-2442431376TammannG. A.Dwarf galaxies in the past19943MatthewsL. D.GallagherJ. S.IIIB and V CCD photometry of southern, extreme late-type spiral galaxies199711451899191910.1086/1186132-s2.0-0000140588MateoM.Dwarf galaxies of the local group199836143550610.1146/annurev.astro.36.1.4352-s2.0-0032349773CapozzielloS.CardoneV. F.TroisiA.Low surface brightness galaxy rotation curves in the low energy limit of Rn gravity: no need for dark matter?200737541423144010.1111/j.1365-2966.2007.11401.x2-s2.0-33847336367PersicM.SalucciP.StelF.The universal rotation curve of spiral galaxies—I. The dark matter connection19962811274710.1093/mnras/278.1.272-s2.0-1542505414PiresM. O. C.De SouzaJ. C. C.Galactic cold dark matter as a Bose-Einstein condensate of WISPs2012201211, article 02410.1088/1475-7516/2012/11/0242-s2.0-84870278604HelouG.MadoreB. F.SchmitzM.BicayM. D.WuX.BennettJ.The NASA/IPAC extragalactic database1991171Springer Netherlands89106Astrophysics and Space Science Library10.1007/978-94-011-3250-3_10McgaughS. S.SchombertJ. M.BothunG. D.De BlokW. J. G.The baryonic tully-fisher relation20005332L99L10210.1086/3126282-s2.0-0034689818McGaughS. S.OlfJ.Local group dwarf spheroidals: correlated deviations from the Baryonic Tully-Fisher Relation2010722124826110.1088/0004-637x/722/1/2482-s2.0-78149294776PecceiR. D.QuinnH. R.CP conservation in the presence of pseudoparticles197738251440144310.1103/PhysRevLett.38.14402-s2.0-9344251553AsztalosS. J.CaosiG.HagmannC.The Axion Dark Matter eXperimentProceedings of the 31st International Symposium on Physics in CollisionAugust 2011Vancouver, CanadaRosenbergL. J.Dark-matter QCD-axion searches201511240122781228110.1073/pnas.13087881122-s2.0-84943302542