Slowly rotating Bose Einstein Condensate galactic dark matter halos, and their rotation curves

If dark matter is composed of massive bosons, a Bose-Einstein Condensation process must have occurred during the cosmological evolution. Therefore galactic dark matter may be in a form of a condensate, characterized by a strong self-interaction. We consider the effects of rotation on the Bose-Einstein Condensate dark matter halos, and we investigate how rotation might influence their astrophysical properties. In order to describe the condensate we use the Gross-Pitaevskii equation, and the Thomas-Fermi approximation, which predicts a polytropic equation of state with polytropic index $n=1$. By assuming a rigid body rotation for the halo, with the use of the hydrodynamic representation of the Gross-Pitaevskii equation we obtain the basic equation describing the density distribution of the rotating condensate. We obtain the general solutions for the condensed dark matter density, and we derive the general representations for the mass distribution, boundary (radius), potential energy, velocity dispersion, tangential velocity and for the logarithmic density and velocity slopes, respectively. Explicit expressions for the radius, mass, and tangential velocity are obtained in the first order of approximation, under the assumption of slow rotation. In order to compare our results with the observations we fit the theoretical expressions of the tangential velocity of massive test particles moving in rotating Bose-Einstein Condensate dark halos with the data of 12 dwarf galaxies, and the Milky Way, respectively.


I. INTRODUCTION
The assumption of the existence of dark matter (DM) is one of the cornerstones of present day cosmology and astrophysics [1][2][3][4][5] . The first evidence for its presence in the Universe was provided by the study of the galactic rotation curves. More exactly, the idea of DM was first suggested to explain the rotation curves of spiral galaxies, whose rotation curves decay far more slowly than one would expect by taking into account the effects of baryonic matter (gas and stars) only. This behavior is considered as evidence for the existence of a supplementary (and exotic) missing mass component, most likely consisting of new particle(s) lying beyond the standard model of particle physics. The rotation curves still represent the most powerful and convincing evidence for DM [6][7][8][9]. But various other astrophysical and cosmological observations have also provided evidence for the existence of dark matter, like, for example, the recent determination of the cosmological parameters from the Planck satellite observations of the Cosmic Background Radiation [10]. These observations have also shown that dark matter cannot be explained by baryonic matter only, thus confirming the standard Λ Cold Dark Matter (ΛCDM) cosmological paradigm. Other types of observations, such as gravitational lensing also require the existence of dark matter for a consistent interpretation of the data [11][12][13]. A particularly strong evidence for the existence of dark matter is provided by the observations of a galaxy cluster called the Bullet Cluster. In this cluster the baryonic matter and the dark matter components are separated due of a collision of its two components that occurred in the past [14]. Measurements of the cosmological parameters by using the Planck data on the cosmic microwave background radiation indicate that the Universe is composed of 4% baryons, 22% non-baryonic dark matter and 74% dark energy [10].
Dark matter models can be divided into three types, cold, warm and hot dark matter models, respectively, by the energy of the particles composing them [15]. The main candidates for dark matter are the WIMPs (Weakly Interacting Massive Particles) and the axions [15]. WIMPs are heavy particles that interact via the weak force [16,17]. Axions are bosons that were first proposed to solve the strong CP problem [18,19]. If the axions form the dark matter, then at low temperature the axion gas will form a Bose-Einstein Condensate (BEC).
There are also other theories that try to explain the observations without introducing dark matter. These theories are based on a modification of the law of gravity at the galactic scales. The earliest one of them is the MOND theory (Modified Newtonian Dynamics) [20]. Modified gravity theories have also been used extensively as an alternative to dark matter [21][22][23][24][25][26][27][28][29][30][31]. An interesting possibility to detect the presence of dark matter is via its possible annihilation. If such a physical process does indeed occur, then a large numbers gamma ray photons and positrons could be produced. Observationally, some excess positron emission in our galaxy has been detected [32][33][34][35][36][37][38][39]. Therefore, it may be possible that the excess positron and gamma-ray emissions could be explained by the annihilation of dark matter with mass m ∼ 10 − 100 GeV [32][33][34]. For a detailed discussion of this problem, as well as of the possibility of alternative interpretations of the observational data see [35][36][37][38][39].
Even that dark matter models can give a good explanation of the qualitative behavior and constancy of the rotation curves, an important contradiction did arise as a result of the in depth comparison of the simulation results with the observations. Data on almost all observed rotation curves show that they rise less steeply than cosmological simulations of structure formation in the standard ΛCDM model in presence of a single pressureless dark matter component do predict. The simulations indicate a central dark matter density profile that behaves as ρ ∼ 1/r (a cusp) [40], while the observed rotation curves indicate the presence of constant density cores [41,42]. This is the so-called core-cusp problem in dark matter physics. Another important open question dark matter models have to face is the "too big to fail" problem [43,44]. By using the Aquarius simulations it was shown that the most massive subhalos in the dark matter halos predicted in the ΛCDM model are inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies [44]. While the best-fitting hosts of the dwarf spheroidals galaxies have 12 < V max < 25 km/s, the λCDM simulations all predict at least ten subhalos with V max > 25 km/s. These results cannot be explained in the framework of the ΛCDM-based models of the satellite population of the Milky Way. The main problem emerging here is related to the densities of the satellites, with the dwarf spheroidals required to have dark matter halos that are a factor of ∼5 more massive than it is observed.
These problems related to the physics of the dark matter may be solvable if one goes one step beyond the standard ΛCDM model, and assume that the dark matter particles may posses some forms of self-interaction. Such a possibility has gained some observational support after the study of the data provided by the observations of 72 cluster collisions, including both 'major' and 'minor' mergers, with the observations done by using the Chandra and Hubble Space Telescopes [45]. Collisions between galaxy clusters can provide an important test of the non-gravitational forces acting on dark matter, and the analysis done in [45] gives an upper limit of the ratio of the self-interaction cross-section σ DM and of the mass m of the dark matter particle as σ DM /m < 0.47 cm 2 /g (95% Confidence Level). A new upper limit on the self-interaction cross-section of dark matter of σ DM < 1.28 cm 2 /g (68% Confidence Level), was obtained in [46]. From a theoretical point of view different self-interacting dark matter models were investigated in [47][48][49][50]. The effects of self-interacting dark matter on the tidal stripping and evaporation of satellite galaxies in a Milky Way-like host were investigated in [51]. Velocity-independent self-interacting dark matter models show a modest increase in the stellar stripping effect with satellite mass, whereas velocity-dependent selfinteracting dark matter models show a large increase in this effect towards lower masses, making observations of ultra-faint dwarfs prime targets for distinguishing between and constraining self-interacting dark matter models. The response of self-interacting dark matter halos to the growth of galaxy potentials using idealized simulations, each run in tandem with standard collisionless Cold Dark Matter (CDM) was investigated in [52]. A greater diversity in the self-interacting dark matter halo profiles was found, as compared to the standard CDM halo profiles. A self-interacting dark matter halo simulated with cross section over mass σ DM /m = 0.1 cm 2 /g provides a good match to the measured dark matter density profile of A2667, while an adiabatically-contracted CDM halo is denser and cuspier. The cored profile of the same halo simulated with σ DM /m = 0.5 cm 2 /g is not dense enough to match A2667. These findings are in agreement with previous results [45] that σ DM /m ≥ 0.1 cm 2 /g is disfavored for dark matter collision velocities in excess of about 1500 km/s. Therefore the possibility that dark matter is a self-interacting component of the Universe cannot be rejected a priori, and physical models whose component particles naturally exhibit this property may provide valuable explanations and suggestions for the dark matter candidates, and for their properties. From both a fundamental theoretical point of view, as well as from a phenomenological perspective, the phys-ically best motivated self-interacting dark matter model can be obtained by assuming that dark matter is in a Bose-Einstein Condensate phase.
The idea that at very low temperatures all integer spin particles may occupy the lowest quantum state, at which point macroscopic quantum phenomena become apparent, was proposed, from a statistical physical point of view by Bose and Einstein in the 1920s [53][54][55]. The Bose-Einstein Condensation process is determined by the quantum mechanical correlation of the gas particles, which implies that the de Broglie thermal wavelength is greater than the mean interparticle distance. The transition to the condensate phase begins when the temperature T of the boson gas is lower than the critical one, T cr , given by [56][57][58][59] T cr = 2π 2 ρ 2/3 cr where m is the particle mass in the condensate, ρ cr is the critical transition density, k B is Boltzmann's constant, and ζ denotes the Riemmann zeta function. It took around seventy years to achieve the experimental realization of the Bose-Einstein Condensation, which was first observed in dilute alkali gases in 1995 [60][61][62]. From a physical point of view the presence of a BEC in an experimental framework is indicated by the appearance of sharp peaks in both coordinate and momentum space distributions of the particles.
Up to now, the main evidence for the existence of BECs comes from laboratory experiments, performed on a very small scale. However, the possibility of the presence of some forms of condensates in the cosmic environment cannot be rejected a priori, and the implications of the possible existence of a condensate state of matter in a astrophysical or cosmological background is certainly worth to investigate. It has been hypothesized that due to their superfluid properties in general relativistic compact objects, like neutron or quark stars, the neutrons may form Cooper pairs, which would condense eventually. Bose-Einstein Condensate stars could have maximum masses of the order of 2 M ⊙ , maximum central densities of the order of 0.1 − 0.3 × 10 16 g/cm 3 , and minimum radii in the range of 10-20 km, respectively. Their interesting physical and astrophysical properties were investigated in [63][64][65][66][67][68][69][70][71][72].
The idea that dark matter is in the form of a Bose-Einstein Condensate was proposed initially in [73], and then rediscovered/reinvestigated, in [74][75][76][77][78][79][80][81][82][83][84]. A systematic study of the properties of the BEC dark matter halos, based on the non-relativistic Gross-Pitaevskii (GP) equation in the presence of a confining gravitational potential, was initiated in [85]. A further simplification of the mathematical formalism of the gravitationally bounded BECs can be achieved by introducing the Madelung representation of the wave function, which allows the representation of the GP equation in the equivalent form of a continuity equation, and of a hydrodynamic Euler type equation. Hence with the use of the Madelung representation we obtain the fundamental result that dark matter can be described as a non-relativistic, Newtonian Bose-Einstein Condensed gas in a gravitational trapping potential, with the pressure and density obeying a polytropic equation of state, with polytropic index n = 1. The validity of the BEC dark matter model was tested by fitting the Newtonian tangential velocity equation to a sample of rotation curves of low surface brightness and dwarf galaxies, respectively.
The thermal correction to the dark matter density profile where obtained in [86]. In [87] it was shown that the density profiles of the Bose-Einstein Condensed dark matter generally show the presence of an extended core, whose formation is explained by the strong interaction between dark matter particles. A further observational test of the model can be obtained by computing the mean value of the logarithmic inner slope of the mass density profile of dwarf galaxies, and by comparing it with the observations. The study of the properties of the Bose -Einstein Condensate dark matter on a cosmological and astrophysical scales is presently a very active field of research . The properties of the Fuzzy Dark Matter, assumed to be formed of a light (m ∼ 10 −22 eV) boson having a de Broglie wavelength λ ∼ 1 kpc, were recently investigated in [126].
If the static properties of the BEC dark matter halos have been studied extensively, their rotational properties have attracted less attention. In [81] the presence of vortices in a self gravitating BEC dark halo, consisting of ultra-low mass scalar bosons was investigated, and it was pointed out that rotation of the dark matter imprints a background phase gradient on the condensate, which induces a harmonic trap potential for vortices. A detailed study of the vortices in rotating BEC dark matter halos was performed in [91], where strong bounds for the boson mass, interaction strength, the shape and quantity of vortices in the halo, and the critical rotational velocity for the nucleation of vortices were found. An exact solution for the mass density of a single, axisymmetric vortex was also found. The effects of rotation on a superfluid BEC dark matter halo were explored in [92], by assuming that a vortex lattice forms. With fine-tuning of the bosonic particle mass and the two-body repulsive interaction strength, it was found that sub-structures on rotation curves that resembles some observations in spiral galaxies could exist. The study of the equilibrium of self-gravitating, rotating BEC haloes, which satisfy the Gross-Pitaevskii-Poisson equations was performed in [103]. Vortices are expected to form for a wide range of BEC parameters. However, vortices cannot form for vanishing self-interaction. The question if and when vortices are energetically favored was also considered, and it was found that vortices form as long as self-interaction is strong enough.
It is the goal of the present paper to study the properties of the BEC dark matter halos in the presence of rotation. Rotation might be a general property of galaxies, whose origin may be traced back to some physical processes in the early Universe. In order to describe the Bose Einstein Condensate dark matter we adopt the Gross-Pitaevskii equation, which gives an effective meanfield description of the multi-particle bosonic system. The mathematical description of the condensate is significantly simplified after introducing the Thomas-Fermi approximation, which allows the description of the dark matter as a gas obeying a polytropic equation of state, with polytropic index n = 1. By assuming a rigid body rotation for the halo, with the use of the hydrodynamic representation of the Gross-Pitaevskii equation, we obtain the basic relation describing the density distribution of the rotating condensate, which naturally generalizes the previously obtained static profile. From the density distribution of the rotating condensate we derive the general representations for the mass distribution, boundary (radius), potential energy, velocity dispersion, tangential velocity, as well as for the logarithmic density and velocity slopes. From the general results we obtain explicit expressions for the radius, mass, and tangential velocity in the first order of approximation, under the assumption of slow rotation. A comparison of our results with the observations is performed by fitting the theoretical expressions of the tangential velocity of massive test particles moving in rotating Bose-Einstein Condensate dark halos with the data of 12 dwarf galaxies, and of the Milky Way galaxy, respectively.
The present paper is organized as follows. The mathematical and physical description of the Bose-Einstein Condensate dark matter is introduced in Section II, where the Gross-Pitaevskii equation, and the Thomas-Fermi approximation are presented. The rotating Bose-Einstein Condensate dark matter structures are investigated in Section III, by using the general approach that allows us to obtain the exact expression of density as expressed in terms of Legendre polynomials. The slowly rotating dark matter halo in the framework of the first order approximation is also considered, and the density profile, as well as the radius are also presented. The gravitational and astrophysical properties of the rotating Bose-Einstein Condensate dark matter halos are studied in Section IV. In this Section we derive the expressions of a number of important astrophysical quantities, like, for example, the mass distribution, potential energy, logarithmic slopes of the density and velocity, which could allow an in depth comparison between the theoretical model and observations. The fitting of the theoretical model with astronomical/astrophysical data is performed in Section V, where we compare the predicted Bose-Einstein Condensate galactic rotation curves with the observational data from 12 dwarf galaxies, and the Milky Way galaxy. We discuss and conclude our results in Section VI.

II. THE BOSE-EINSTEIN CONDENSATE DARK MATTER MODEL
In the present Section we briefly introduce the fundamental physical concepts related to the Bose-Einstein Condensation, as well as the basic equations describing the rotating condensate. It has been shown that if the dark matter is composed of ultralight boson particles with mass m ∼ 10 −22 eV and wavelength λ ∼ 1 kpc, then the transition temperature to a Bose-Einstein Condensate is of the order T c ∼ 10 9 K [83]. Hence, if dark matter is composed of Bose particles, like, for example, the axion, it is quite natural to assume that dark matter is in a Bose-Einstein Condensate state. For a recent discussion of the arguments from particle physics that may motivate the existence of the ultra-light dark matter, as well as of its properties and astrophysical signatures see [126].

A. The Gross-Pitaevskii Equation
A Bose-Einstein Condensate is a phase of matter in which all the particles are localized in the ground state. The Bose-Einstein Condensation occurs for particles that have integer spins, and obey the Bose-Einstein statistics. We will consider in the following that the bosons are weakly interacting, and that their interaction is described by a two-body interparticle potential. We start our analysis by writing down first the many-body Hamiltonian of the bosonic system, whereΨ(r) andΨ † (r ′ ) are the annihilation operator and the creation operator at the position r, respectively, V denotes the interparticle interaction potential, and m is the mass of the particle in the condensate. In the case of a rotating dark matter halo, the external potential V ext is the gravitational potential, and V rot is the effective centrifugal potential. In the following we will adopt the comoving frame, that is, the frame that's rotating with the same speed as the system.
To simplify the mathematical formalism, we introduce the mean field description, in which we decompose the field operator in the formΨ(r) = Ψ 0 +Ψ ′ (r), and treat the operatorΨ ′ (r) as a small perturbation. Then for the mean field component Ψ 0 we have Ψ 0 = N/V , where N is the total particle number, and V is the volume. Hence Ψ 0 is equal to the square root of the number density of the particles [56].
In the general time-dependent case, the field operator in the Heisenberg picture is given bŷ where ψ(r, t) = Ψ (r, t) is also called the condensate wave function. Then for the number density of the condensate we have ρ n (r, t) = |ψ(r, t)| 2 . The normalisation condition is N = ρ n (r, t)d 3 r.
In the Heisenberg representation the equation of motion of the field operator is (4) In the theory of the Bose-Einstein Condensation one usually assumes that the interparticle interaction is a short range interaction, and hence we can write the interaction potential as being proportional to a constant, which is related to the scattering length, times a Dirac delta function [127]: where a is the scattering length. To obtain a more general description, in the following we introduce the function g(|ψ(r, t)| 2 ) to describe the self-interaction term [85].
In the standard approach to Bose-Einstein Condensation the self-interaction is assumed to have a quadratic form, so that g(|ψ(r, t)| 2 ) = 1 2 λ|ψ(r, t)| 4 [56]. With this approximation of the potential, and with the use of the mean field approximation, by integrating over Eq.(4) we obtain the Gross-Pitaevskii equation, describing the main properties of a Bose-Einstein Condensate, as To give a more direct physical interpretation of the Gross-Pitaevskii equation, we introduce Madelung representation of the wave function, which separates the wave function into two components, its magnitude, and a phase factor, respectively. The function S(r, t) has the dimensions of an action. Then in the Madelung representation we have [59] i ψ Then the Gross-Pitaevskii equation is separated into two parts. From the imaginary part we obtain, where v = ∇S m is the velocity of the quantum fluid. This is the continuity equation. On the other hand, from the real part we obtain the equation which is the momentum conservation, or the hydrodynamic Euler equation [117]. In Eq. (11) P is the thermodynamic pressure, which is related to the mass density ρ = ρ n m of the condensate by a barotropic type equation of state [85,127] The term σ Q is given by √ ρ , and its divergence is called the quantum stress tensor.
Therefore the two equations describing the evolution of a Bose-Einstein Condensate are the continuity and the Eular equations of classical fluid dynamics, which describe viscosity free flows. Also, we can see from the definition of the velocity that the flow must be automatically irrotational. We will discuss this issue later.

B. Thomas-Fermi Approximation
When the number of particles in the condensate get large enough, the contribution to the energy of the quantum pressure term ∇ · σ Q can be neglected except near the boundary [58]. Then the equations describing the condensate become purely classical in their mathematical form, even that their physical interpretation must be given in the framework of quantum statistical physics.
If we consider the Bose-Einstein Condensate to be static, or have a rigid body rotation, and we work in the corotating frame, then v = 0. Thus from Eq. (11) we obtain where ∇h(ρ) = (1/ρ) ∇P (ρ). For the exterior potential we assume that it is the gravitational potential, V ext = V grav = V , and that it satisfies the Poisson equation, where G is the gravitational constant. For the rotational potential we adopt the expression where ω is the angular velocity of the dark matter halo. In the case of the quadratic nonlinearity we have g ′ (ρ) = λρ/m, and g(ρ) = λρ 2 /(2m 2 ), respectively. Thus for the equation of state of the condensate we obtain From the equation of state of the condensate we can see that P ∝ ρ 2 , and since for a general polytropic equation of state Γ = 1 + 1/n = 2, it follows that the polytropic index of the condensate is n = 1.
Hence with the use of the equation of state of the BEC dark matter from Eq. (13) we obtain the equation describing the variation of the density of the rotating dark matter halo as where Eq. (18) is the Helmholtz equation. If the system had a different polytropic index n = 1, instead, we will get a general Lane-Emden equation [85], which will be nonlinear.
If the halo is nonrotating, ω = 0, and, under the assumption of spherical symmetry, Eq. (18) has the solution [85] where A 0 is an integration constant. One can obtain the radius R of the static halo from the boundary condition ρ(R) = 0, which gives The central density is ρ c = ρ(0) = A 0 . We can see from the expression of the density that the radius only depends on the mass and scattering length of the particles. The size of the halo is independent of the central density. We can determine ρ c from the normalisation condition, d 3 rρ = M , thus obtaining Once the BEC dark matter density profile is known, all the global parameters of the BEC dark matter halo (mass, radius, central density), as well as the rotational speeds of particles in stable circular motion can be obtained in an exact form. These results open the possibility of the observational test of the BEC dark matter model.

C. Emergence of Vortices
We have already mentioned that from the definition of the velocity v = ∇S, the flow must automatically satisfy the condition ∇ × v = 0, and hence the quantum BEC fluid motion should be irrotational. However, if we expect the dark matter halo to be rotating, we must introduce singularities of the vorticity, and therefore the halo will contain a vortex lattice.
First, let us recall the concepts of vortex and vorticity. A vortex is a region of fluid in which the flow is rotating around an axis line. Its vorticity is the curl of the velocity, w = ∇ × v. If the fluid rotates like a rigid body with an angular velocity Ω = (0, 0, Ω), we have the velocity v = Ω × r = (−Ωy, Ωx, 0) and the vorticity w = ∇ × v = 2Ω. A vortex can also be irrotational, if it has angular velocity Ω = (0, 0, αr −2 ), its vorticity will be 0 except at the axis line, where the vorticity will be infinite. If a fluid is irrotational, then although the particles have an angular velocity, they will not rotate over themselves.
If we expect the dark matter halo to rotate like a rigid body, it will give rise to a vortex lattice [91,92]. It was already shown in laboratory experiments that such vortex lattices arise when an asymmetry is introduced [128,129]. It has been shown experimentally that vortices arise at a critical angular velocity, at which the energy of the system is lower if it generates vortices [91]. When the angular velocity gets higher, instead of generating a bigger vortex, a lattice of several vortices will be generated [128].
How a vortex will influence the properties of the dark matter halo was studied in [91]. Significantly, a core appears at the center of the vortex, and within the core the density is zero. In [92], this feature was used to explain the wiggles in the rotation curves of the galaxies.
For simplicity, we will not study the vortices in this paper. We will assume that the halo rotates like a rigid body, we will ignore the cores that generate inhomogeneities in the density profile, and we will concentrate our attention on how rotation causes the deformation of the halo, and influences its observable physical properties.

III. DEFORMATION OF THE SLOWLY ROTATING BEC DARK MATTER HALOS
There have been many studies considering the problem of rotating polytropes, using different approximations and building different models, like, for example, in [130,131]. For a detailed discussion of the rotational properties of n = 1 polytropes see [132]. We will assume the halo to be rotating slowly, and thus a first order approximation is sufficient. Chandrasekhar has worked on this problem in 1933 [133,134], and we will mainly follow his method (for a comparative study of the different approaches to the rotation problem see [131]).
A. The energy density and the gravitational potential of the rotating BEC dark matter halos We have already obtained the Helmholtz equation (18) describing the matter distribution inside a rotating BEC dark matter halo. Its general solution is given by where j l (X) are spherical Bessel functions, which are the solutions to the radial part of the equation, while P l (cos θ) are the Legendre polynomials -they are the solutions to the angular part of the Helmholtz equation. In the solution, we have neglected all the terms of odd order, since it has been shown that a rotating mass must be symmetric about its equator. This result is called Lichtenstein's theorem [135]. To determine the coefficients A 2l in the first order of approximation, we will write down the solution for the gravitational potential, and use the fact that it is continuous at the boundary of the dark matter halo. The gravitational potential satisfies the equations where r 0 is the radius (boundary) of the halo. We can represent the potential in a general form as Then the gravitational potential equation inside the halo becomes In the following we denote ξ = kr, and µ = cos θ, respectively. Thus Eq. (27) takes the form Now we separate the terms in the above equation. For the 0th order term we obtain The solution of the above equation is For the higher order terms, since the P 2l 's satisfy the equation we obtain the equations Since ∇ 2 (j 2l P 2l ) = −j 2l P 2l , and by taking into account that the Bessel functions j 2l satisfy the equation it follows that By substituting the above relations into Eq. (32) we obtain A particular solution of the above equation is The regular solution of the equation which behaves well near ξ = 0 is where C 2l are arbitrary integration constants. Hence the general solution for the V 2l 's is giving for the gravitational potential V the general solution

B. The continuity conditions
To determine the B 2l 's, we consider the behavior of the pressure P , which is given by the static limit of the Euler equation (11) as or, equivalently, After substituting the pressure with the use of the Bose-Einstein Condensate equation of state, and equating the coefficients of the series expansions, we obtain , B 2l = 0, l = 1. Therefore To the accuracy we are working, this potential should be continuous with the external potential on the sphere of radius ξ 1 = kR, where R is the boundary of the nonrotating sphere, To determine the coefficient A 2 , we will equate the potentials and their first derivatives at ξ 1 , Hence we obtain and for all the l = 0, 1, A 2l = 0. A 0 is not determined by the boundary condition of potential at radius R, since an arbitrary constant can always be added to the gravitational potential. We can determine it by the boundary condition of the matter density at the center of the dark matter halo, Hence we find C. The first order corrections to density and radius With the use of the expressions for the coefficients obtained above, after substituting ξ 1 = π and ξ = kr, respectively, we find the first-order correction to the density ρ of the rotating Bose-Einstein Condensate dark matter halo as or, equivalently, where We will describe the effect of the rotation on the structure of the dark matter halo by using the dimensionless parameter Ω 2 . In particular, the case of slow rotation corresponds to values of Ω 2 so that Ω 2 << 1. In the limiting case Ω 2 → 0 we recover from Eq. (51) the static limit for the halo density, ρ(r)ρ c = j 0 (kr). The angular momentum is usually described by using the dimensionless spin parameter λ = J|E| 1/2 /GM 5/2 , where J is the angular momentum, and E is the gravitational energy of the halo [136]. From a physical point of view the spin parameter λ is the ratio of the actual angular momentum of the galaxy and of the maximum angular momentum value needed for rotational support. Simulations have shown that λ is in the range 0.02 − 0.10 [136].
The comparative variation of the density in the nonrotating and rotating cases at θ = π/2 (corresponding to the equatorial plane) is presented in Fig. 1.
In order to obtain the first order approximation we will further expand this expression to the first order of ω 2 /2πG, and obtain This has the same form as in Chandrasekhar's work [133]. Hence the boundary of the slowly rotating Bose-Einstein Condensate dark matter halo is located at For the equatorial radius of the dark matter halo we obtain In the the non-rotating case only the first term exists. Rotation adds an expansion and an ellipticity term to the halo radius.

IV. GRAVITATIONAL AND ASTROPHYSICAL PROPERTIES OF ROTATING BOSE-EINSTEIN CONDENSATE DARK MATTER HALOS
In the present Section we will obtain some basic gravitational and astrophysical properties of the slowly rotating Bose-Einstein Condensate dark matter halos, which could allow an in depth comparison of the theoretical model with the astronomical observations. In particular, we will consider the mass distribution within the halo, as well as to its gravitational potential energy. Moreover, we will concentrate on astrophysical parameters like velocity dispersion, logarithmic density and velocity slopes, and the tangential velocity expression, which allows a detailed comparison of the model predictions with observational data.
A. Mass and gravitational potential of the slowly rotating BEC halo As we have already seen, the general solution for the matter density distribution ρ inside the rotating halo in spherical coordinates is given by The boundary (radius) r 0 of the halo is defined as the surface whose points satisfy the condition ρ(r 0 ) = 0. The equation for the density involves infinitely many terms in general.
The mass profile m(r, θ = π) within a radius r is given by which involves an integration of the spherical Bessel functions, which can be done by using the relation where 1F2 is the is the regularized generalized hypergeometric function p F q (a; b; z) / (Γ (b 1 ) ...Γ (b q )). Thus we obtain for the total mass the expression For dark matter halos located within a radius r ≤ π/k − 3ω 2 /4Gρ c k, we can calculate the mass distribution within radius r as being given by As compared to the non-rotating case, a second term appears, which is due to the presence of the rigid body type rotation. The mass profile within radius r is bigger than in the non-rotating case, and it depends on the central density ρ c . These results are consistent with the slower decay of the density profile for larger values of the radial coordinate r. The variation of the dimensionless ratio m(ξ)/M * , where M * = 4πρ c /k 3 is represented in Fig. 2. The total mass of the condensate in the first order of For a given galaxy, the total halo mass is fixed by the physical processes leading to its formation, and it is unchanged due to the presence of rotation. But, as one can see from Eq. (64), a rotating halo is able to hold more mass than a static one, and in this sense the rotation of the dark matter halo becomes a stabilizing factor against gravitational collapse.
To perform the integration in the accuracy of the first order of ω 2 /2πG, we first integrate over r, then expand the result in the first order, and then we perform the integration over µ. This procedure does not affect the final result obtained by doing the series expansion after performing the full integration, since ω 2 /2πG is independent of r and µ, and thus the power of it is not changed after each integration.
Then we can express the central density ρ c as a function of the total mass and angular velocity as The total volume of the halo becomes From the above expression of the volume it follows that due to rotation the halo has expanded. The added volume is proportional to ω 2 , and inversely proportional to ρ c . With the help of the total mass and of the volume we obtain the mean density of the BEC halo as which is smaller than in the non-rotating case.
In the previous Section, we have already calculated the gravitational potential inside the halo, We can determine the constant by using the continuity of the potential near radius ξ 1 = kR, Thus we obtain Hence the gravitational potential is given by Using the gravitational potential, we can calculate the gravitational binding energy U defined as Hence the total gravitational potential energy of the BEC dark matter halo is given by As compared to the non-rotating case, the total gravitational energy has a second negative term, and hence for certain values of ρ c it is lower than in the non-rotating case.
The centrifugal potential of the rotating dark matter halo is Hence we can calculate the centrifugal potential energy as The centrifugal potential energy is always lower than 0. Hence the effective potential energy in the corotating frame can be obtained as B. Velocity dispersion of particles in slowly rotating BEC halos In the following we consider the dynamics of a collection of particles (stars) in the gravitational field of a Bose-Einstein Condensate dark matter halo. The statistical properties of the motion are described by the Jeans equation, which for a system of particles with number density n = n(x i , t) is given by [137] where the notation < > means an average at a given point and time (x, t), and For a static spherical symmetrical system, one can further simplify the radial Jeans equation by adopting the assumptions: 1) steady-state hydrodynamic equilibrium, which implies ∂vj ∂t = 0, and < v r >= 0, respectively, 2) < v θ >=< v φ >= 0, and σ 2 rθ = σ 2 rφ = σ 2 θφ = 0, respectively, which follows from the spherical symmetry of the system, and 3) a single tangential velocity dispersion for all directions σ 2 tt = σ 2 θθ = σ 2 φφ . Hence when the tangential velocity dispersion tensor is isotropic, σ 2 rr = σ 2 tt = σ 2 , the Jeans equation reduces to 1 ρ n ∂ρ n σ 2 ∂r = − ∂V ∂r . (79) By assuming that all particle have the same mass, the particle velocity dispersion is obtained as For a rotating system, we can still adopt the following set of assumptions, 1) Thus in presence of rotation the radial Jeans equation becomes Therefore for the mean value of the square of the radial velocity we obtain v 2 r (r, θ) = 1 ρ(r, θ) r0 r ρ(r, θ) ∂V (r, θ) ∂r ′ − ρ(r, θ)ω 2 r ′ sin 2 θ dr ′ . (82) In order to find an explicit expression for v 2 r (r, θ) , we expand the integrand to the first order in ω 2 /2πG. Hence for r < r 0 we find v 2 r (r, θ) = 2ρ c Gπ k 3 r sin kr + ω 2 24k 5 r 3 cos(kr) − 12k 2 r 2 + 10P 2 (cos θ) 6k 4 r 4 + π 2 (6 − 7k 2 r 2 ) + 12k 2 r 2 + 10P 2 (cos θ)π 2 (−6 + k 2 r 2 ) cos(2kr) + 24k 3 r 3 sin(kr) + 30kP 2 (cos θ)π 2 r sin(2kr) . (83) We define the kinetic energy K of the halo in terms of the average velocity dispersion σ v as If the velocity dispersion is a constant, and at the boundary K takes the value, For the ratio of the kinetic and potential energy, after expanding to the first order of ω 2 /2πG, we obtain Comparing this expression to the non-rotating case, we can see that rotation generates a second (negative) term in the parentheses, and hence in the presence of rotation the ratio of the kinetic and potential energy of particles in motion in rotating BEC dark matter halo is lower than in the non-rotating case.
As compared to the non-rotating case, a third term is added to the expression of the logarithmic slope of the density. The third term is smaller than the the nonrotating value of α BEC , and varies with θ. This seems reasonable, since the rotation pushes the matter outward, and the ratio of the center density to the density at larger radii is smaller.
Since the logarithmic density slope varies with θ, it's no longer convenient to define the core radius R core as α BEC (R core ) = 1. If we simply define it as R core = nR = nπ/k, where n is a constant which must be determined from observations (n = 0.6 gives the value of the core radius in the static case [87]), we can define the mean value of the logarithmic density slope within the radius 0 ≤ r ≤ R core as The variation of the mean value of the logarithmic slope of the density α BEC is represented, as a function of Ω, and for different values of R core , in Fig. 3.
Also, we can calculate the density at the core radius, The core density is smaller as compared to the nonrotating case at θ = 0, and larger at θ = π/2. For a given ρ c , the quantity ρ c R core can be obtained as and its value depends on ω 2 .

D. Tangential velocity of test particles in slowly rotating BEC ark matter halos
In the Newtonian approximation the tangential velocity of a test particle moving in the Bose-Einstein condensed dark matter halo is given by In the slow rotation approximation, and for r ≤ π/k − 3ω 2 /4Gρ c k, we have We can see that due to rotation a second positive term is added to the expression of the tangential velocity, so that at a distance r, the tangential velocity of a test particle becomes higher as compared to the non-rotating case. This is because the mass profile within the radius is higher. The variation of the ratio V tg (ξ)/V * as a function of ξ = kr, where V * = 4πGρ c /k 2 , is represented in Fig. 4. At the equatorial boundary of the halo we have, We can see that the second term proportional to ω 2 is also positive. Hence in a rotating BEC dark matter halo, the tangential velocity of a test particle is larger than in the static one. Since the ratio of the density at larger and smaller radii is greater than in the non-rotating case, it follows that the tangential velocity is bigger at larger radii than in the non-rotating case.
An important observational quantity is the logarithmic slope of the tangential velocity β tg , defined by Generally, β tg can be obtained as At the center, where r = 0, we have β tg (0) = −1. This gives the same result as in the static case, since ω = 0 at the center. By expanding the logarithmic slope of the tangential velocity to the first order in ω 2 /2πG, we obtain β tg ≃ 1 2 1 + k 2 r 2 kr cot(kr) − 1 + k 3 r 3 3kr cos(kr) + (−3 + k 2 r 2 ) sin kr 6(−kr cos kr + sin kr) 2 Ω 2 .
This is smaller than the value without rotation. The variation of the logarithmic slope of the tangential velocity is shown, for different values of the dimensionless parameter Ω, in Fig. 5.

V. GALACTIC ROTATION CURVES IN THE ROTATING BOSE-EINSTEIN CONDENSATE DARK MATTER MODEL
As a next step in our analysis we compare the predictions of the slowly rotating Bose-Einstein Condensate dark matter model with the observational data obtained for a sample of HSB, LSB and dwarf galaxies. From a realistic astrophysical point of view, the matter content in a galaxy consists of a distribution of baryonic (normal) matter, obtained as the algebraic sum of the masses M star of the stars, of the ionized gas with mass M gas , of the neutral hydrogen of mass M HI etc., as well as of dark matter of mass M DM . In the following we assume that dark matter is in the form of a slowly rotating Bose-Einstein Condensate. Therefore the total mass of the galaxy is can be obtained as . is the total baryonic mass in the galaxy. The rotation velocity of a test particle v rot is given by the sum of the different matter contributions, as where v 2 gas and v 2 stars are the contributions of the baryonic gas and stars, respectively, while v 2 halo is the dark matter contribution, which we assume to be given by Eq. (97). Hence the contribution of the rotating BEC dark matter halo to the rotational velocity can be represented as In Eq. (103), R is the radius of the static BEC dark matter halo, which is fixed by the numerical values of the scattering length a and the mass m of the dark matter particle.

A. Fitting results
In order to test our model, we apply it to small nearby dwarf galaxies (size ≤ 12 kpc) and our Milky Way galaxy. We use the data of the Spitzer Photometry and Accurate Rotation Curves (SPARC) obtained in [138] for investigation. The baryonic components (bulge, disk and gas) included in the data are obtained by observations using the homogeneous surface photometry at 3.6 µm [138]. The rotation curve contributions of the baryonic components are only determined by the mass-to-luminosity ratios Υ ∼ 1M ⊙ /L ⊙ of the disk and of the bulge. Nevertheless, the surface brightness and the resultant rotation curves obtained may also suffer from some systematic uncertainties, such as the irregularities in brightness profiles, uncertainties in inclination and patchy distribution of gas [138]. We choose the candidates which are Hubble stage T = 0 − 6 galaxies, bulgeless galaxies, and with the distance to the galaxy D ≤ 20 Mpc. Generally speaking, these galaxies have less diffuse features and smaller uncertainties in the observational data. Based on these criteria, we consider 12 dwarf galaxies for testing the present Bose-Einstein Condensate dark matter model.
Based on Eq. (103), we have altogether four free parameters for fitting (R, ρ c , ω, Υ). The mass-toluminosity ratio of the disk Υ mainly affects the central rotation curve, while the other three control the entire shape of rotation curve. Here, we define the reduced χ 2 as where N dof is the number of degrees of freedom, v rot,i is the calculated rotation velocity, v obs,i is the observed rotation velocity, and σ i is the observational uncertainties of the rotation velocity. By minimizing the reduced χ 2 value, we can obtain the best-fit values of R, ρ c , ω and Υ for each of the considered dwarf galaxies (see Table I).  In Fig. 6, we show the best-fit rotation curves for the 12 dwarf galaxies.
We can see that most of the best-fit values of R fall into a small range R = 5 − 8 kpc. Although the bestfit R for the UGC07151 galaxy is quite small (R = 2.9 kpc), the acceptable range of R is R = 2.5 − 6.5 kpc for χ 2 red ≤ 4. Therefore, the results for the 12 dwarf galaxies are still consistent with a fixed value of R ≈ 6.5.
With the help of Eq. (64), which gives the mass of the Bose-Einstein Condensate dark halo, we can predict the mass of the considered sample of galaxies. In the limit of small rotational values we obtain for the total mass the expression The predicted masses of the 12 dwarf galaxies are presented in Table II. We also apply our model to our Milky Way galaxy. By using the data obtained in [139] and the baryonic model used in [140], we perform a similar fit for the Milky Way data (see Fig. 7). Since the mass-to-luminosity ratios of the bulge and disk are included in the baryonic model, we only have three free parameters to fit. The best-fit values are R = 6.9 kpc, ρ c = 3.0 × 10 −24 g cm −3 and ω = 2.6 × 10 −16 s −1 (χ 2 red = 0.54). We also show the fit for R = 10 kpc (χ 2 red = 1.02) for reference in Fig. 7. Generally speaking, the best-fit values for the Milky Way galaxy and the dwarf galaxies are consistent with each other.

VI. DISCUSSIONS AND FINAL REMARKS
In the present paper we have investigated the possibility that the BEC dark matter halos could be in fact rotating. Presently, there is a common paradigm in the study of superfluidity according to which rotational motion in a Bose-Einstein condensate can exist only in the presence of quantized vortices. However, in a recent numerical study [141] it was shown that the merging of two twodimensional concentric Bose-Einstein condensates with axial symmetry may lead to the formation of a spiral dark soliton. This happens if one of the two condensates has a non-zero initial angular momentum. The spiral dark soliton makes possible the transfer of angular momentum between the two condensates, and allows the merged condensate to rotate, even in the absence of quantized vortices. A similar physical process could have acted on a galactic scale, since galactic collisions and merging, which occurred frequently in the early stages of the evolution of the Universe, may have favored angular momentum transfer between galaxies. Therefore the rotation of BEC dark matter halos may be possible even in the absence of quantized vortices.
By using the Thomas-Fermi approximation, the basic evolution equations of a rotating BEC halo reduce to the case of the rotating polytropic spheres with index n = 1. Many mathematical methods have been developed for the investigation of such systems. For exam-ple, in [131] the first exact analytic solution for an arbitrarily rotating gaseous polytrope under the assumption of the oblate spheroidal shape was derived. To obtain the solution the authors have adopted oblate spheroidal coordinates (ξ, η, φ), which are related to the Cartesian coordinates by means of the transformation where R e and R p are the equatorial and polar radii, respectively. In [130] index n = 1 polytropes were studied under the assumption that the shape of such a rotating fluid sphere is spheroidal. By introducing the spheroidal coordinates (x, η), defined as κ 2 r 2 = x 2 + c 2 1 − η 2 , and κz = xη, where κ and c are constants, one can express the gravitational potential by using spheroidal wave functions.
The effect of the rotation of the dark matter halo in the framework of the BEC dark matter model has also been previously investigated. Two classes of models for rotating haloes were investigated in [103], in order to analyze their stability with respect to vortex formation. In the first model haloes were modeled as homogeneous Maclaurin spheroids, while in the second one an n = 1 polytropic Riemann-S ellipsoid was considered. Generally, it was shown that BEC haloes in the polytropic Thomas-Fermi regime typically form vortices. The dynamics of the rotating Bose Condensate galactic dark matter halos, made of an ultralight spinless boson gas was investigated in [104]. The basic approach in this study was to obtain the numerical solution of the system of the coupled Gross-Pitaevskii and Poisson equations, i ∂Ψ/∂t = − 2 /2m ∆Ψ + V Ψ + 2π 2 a/m 2 |Ψ| 2 Ψ and ∆V = 4πGm|Ψ| 2 , respectively, in Cartesian coordinates. It was found that ultralight spinless boson dark matter candidates can describe well the galactic astrophysical properties at local scales with the addition of angular momentum to halos.
In our study we have considered the case in which the halo has an overall rigid body rotation, and we have studied the astrophysically relevant properties that such a rotating halo may have. As a first (and basic) result we have obtained the rotational corrections to the halo density profile, due to rotation. In our investigations we have followed the approach initiated in [133], and we have considered the rotation problem in spherical coor-dinates. While [133] considers polytropic systems with arbitrary n, in the present paper we have systematically investigated the n = 1 configurations. In this case the general expression of the density involves an infinite summation over a set of radial Bessel functions, and angular Legendre polynomials. We have restricted our investigations of the astrophysical properties of the BEC halos to the slow rotation case, when the deformation of the halo is small, and the rotation parameter Ω 2 satisfies the condition Ω 2 << 1. In this case explicit, and simple expressions of the density of the rotating halo can be obtained. The knowledge of the density distribution is the first step in the investigation of the physical properties of BEC dark matter halos. The mass distribution inside the halo can then be obtained, and the knowledge of the mass profile leads to the expression of the tangential velocity of test particles, following circular orbits around the galactic center. In our approach we have obtained the rotational velocity in the first approximation of Ω 2 , an approximation that may allow the investigation of the effects of the slow rotation on the halo.
The tangential velocity, as well as the density profile, essentially depend on three physical parameters: the central density of the dark matter, the radius of the static density profile, and the angular velocity of the galaxy, respectively. In order to obtain these parameter, we have compared our theoretical results with a (small) set of observational data, and we have fitted our BEC rotating model with 13 observed rotation curves. The sample consisted of 12 dwarf galaxies, and the Milky Way galaxy. The fittings of the 12 dwarf galaxies provided a range of central densities of ρ c = (0.7 − 6) × 10 −24 g/cm 3 , indicating a relatively inhomogeneous distribution of the central densities. The angular velocities also present a relatively large spread, ranging from 0.1 × 10 −16 s −1 to 5 × 10 −16 s −1 , presenting an order of magnitude variation. The predicted masses of the halos also did present a large variation, from 0.17 × 10 10 M ⊙ to 4.12 × 10 10 M ⊙ , implying a difference by a factor of about 24 between the smallest and the highest galactic mass.
One of the important parameters in the theoretical models of BEC dark matter is the radius R of the static (nonrotating) halo. The value of R is determined by the scattering length a and the mass m of the dark matter particle. As such, R must be a universal constant, and its constancy as proven by the observations, could represent a strong argument in the favor of the condensate dark matter models. In the considered sample of 12 dwarf galaxies R did vary between 2.9 and 8.2 kpc, respectively, by a factor of around 3. However, for this sample of considered galaxies, a fixed value of R ≈ 6.5 kpc can give a good fit to all considered observational data. Moreover, such a value of R can give a good fit even to the Milky Way rotation curve, which extends up to 30 kpc, a result which, taking into account the numerous uncertainties in the data and in the (baryonic) fitting model, suggests that the assumption of an universal value of R cannot be ruled out by present day observations. Our theoretical results have also shown that introducing an angular velocity gives a smaller central density ρ c , and a larger radius at the boundary of the halo. Also, slow rotation slightly increases the total mass of the angular halo. On the other hand, a good knowledge of the baryonic mass distribution in the galaxies, and of the total central mass density can give, via the fitting of the galactic rotation curves, a good indication on whether the halo is rotating, or not, and the value of its angular velocity.
In the present analysis we did not include the contribution to the gravitational potential of the baryonic galactic matter when we derived the mass distribution of the dark matter halo. This require to modify the Poisson equation to ∆V = 4πG (ρ + ρ b ), where ρ b is the density of the baryonic matter, and to systematically include the effect of the baryonic matter in the model. In particular, the Helmholtz type equation describing the density distribution of the dark matter would also depend on the density of the baryonic matter. For more luminous galaxies, it may be necessary to take this effect into account.
The analysis of the galactic rotation curves alone cannot determine the basic physical properties of the condensate dark matter. Alternative physical effects must be taken into account to fully determine the properties of the dark matter particle. One is such physical effects would be the study of the gravitational lensing by BEC halos. It was already shown in [85] that the BEC dark matter gives a very different prediction for gravitational lensing as compared to other models of dark matter. One might also include the effects of the vortices in the theoretical analysis, and such an inclusion might explain the wiggles in the rotation curves [92].