## Dynamical friction in rotating ultralight dark matter galactic cores

V.M. Gorkavenko 1 , ∗ , O.V. Barabash , T.V. Gorkavenko , O.M. Teslyk , 1 1 1 A.O. Zaporozhchenko , Junji Jia , A.I. Yakimenko 1 3 1 4 5 , , , E.V. Gorbar 1 2 ,

1 Faculty of Physics, Taras Shevchenko National University of Kyiv,

64, Volodymyrs'ka str., Kyiv 01601, Ukraine

- 2
- Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 14-b Metrolohichna str., Kyiv 03143, Ukraine
- 3 School of Physics and Technology, Wuhan University, 299 Bayi Road, Wuhan, Hubei Prov., China 4 Dipartimento di Fisica e Astronomia 'Galileo Galilei', Universit'a di Padova,

via Marzolo 8, 35131 Padova, Italy

- 5 Istituto Nazionale di Fisica Nucleare, Sezione di Padova, via Marzolo 8, 35131 Padova, Italy

## Abstract

Dynamical friction and stellar orbital motion in spiral galaxies with dark matter composed of ultralight bosons in the state of rotating Bose-Einstein condensate (BEC) are studied. It is found that the dynamical friction force is significantly affected by the topological charge of the vortex structure of the BEC core with the strongest effect at distances near the galactic center. It is also shown that the ultralight dark matter self-interaction plays an important role in studying the dynamical friction.

Keywords: dark matter, ultralight dark matter bosons, vortex structure

## 1 Introduction

The problem of nature and composition of dark matter (DM) is one of the most pressing problems of modern physics and astrophysics. DM makes up 26.8% of the mass-energy composition of the Universe [1] and is assumed to consist of beyond the Standard Model particles which do not interact with Standard Model particles or interact with them very feebly. For a recent review of observational evidence for the distribution of dark matter in galaxies, see [2]. There are many candidates for the explanation of DM including primordial black holes [3, 4]. Among popular candidates for DM particles are weakly interacting massive particles (WIMPs) with mass of ∼ 100 GeV [5], sterile right-handed neutrinos with masses of several keV [6, 7], and axion-like (pseudoscalar) particles with masses of µ eV [8]. At the moment, none of these particles has been directly observed and the corresponding models are known as the cold dark matter (CDM) models [9].

∗ Corresponding author. Email address: gorkavol@gmail.com (Volodymyr Gorkavenko)

In this paper, we will consider another attractive possibility of ultralight DM composed of spinless bosons with masses ∼ 10 -23 -10 -21 eV whose enormous occupation number results in the formation of the Bose-Einstein Condensate (BEC) with superfluid properties (for a review and discussion of the observable restrictions on the bosonic particle mass, see, e.g., [10-12]). Ultralight dark matter (ULDM) models have very interesting phenomenology. Because of the very large de Broglie wavelength of order kpc, the formation of small-scale structures in ULDM models is suppressed, unlike CDM models which predict too many dwarf galaxies and too much dark matter in the innermost regions of galaxies.

The fascinating physics of superfluid DM leads to a number of interesting ULDM structures whose observational signatures might help to elucidate the basic properties of DM particles. Recently, it was found that wavelike ULDM may resolve lensing anomalies of Einstein rings in gravitationally lensed images [13]. Further, quantum wave interference of CDM without self-interaction (also known as fuzzy DM [14], ultralight pseudo NambuGoldstone boson appearing in the late time cosmological phase transition theories was proposed in [15] as a natural realization of fuzzy dark matter) produces stochastically fluctuating density granulation which could drive significant disc thickening, providing a natural explanation for galactic thick discs [16]. It was shown that repulsive selfinteractions, even weak, could be important for galactic halo composed of ultralight bosons [17]. A scalar field model of dark matter in spiral galaxies was proposed in [18] which produces the stellar orbital velocity profile in good agreement with observations. For purposes of this work, it is important that ULDM could form stable vortex structures [19] which survive even after head-on collision [20]. Repulsive self-interaction is critical for the existence of stable vortex solutions.

Given that direct observation of DM remains elusive, we focus on the influence of dynamic friction on baryonic matter. We consider how different soliton structures of bosonic ULDM affect the motion of stars in galaxies. We concentrate on stellar motion because it admits semi-analytic treatment of dynamical friction, unlike the case of motion of supermassive black holes in the ULDM environment near the galactic centers. Indeed, a recent analysis in [21] showed that the presence of a moving supermassive black hole significantly affects ULDM soliton, leading to dynamical friction oscillations, requiring numerical analysis. Further, it was shown that stochastic density fluctuations of ULDM may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred parsecs and can heat and expand the central regions of galaxies [22].

The classical work by Chandrasekhar [23] on dynamical friction experienced by a moving star due to the fluctuating gravitational force of neighbouring stars was generalized to the case of the motion of stars in ULDM in the absence of self-interaction (the case of fuzzy DM [14]) in Ref. [24] where analytic expression for dynamical friction was derived. Later the obtained results were numerically checked in [25,26]. The ULDM self-interaction modifies the dynamical friction force and its role for the star's motion in self-interacting ultralight dark matter was studied in [27].

The ground and vortex states of the ULDM yield different DM density and velocity distributions. Therefore, the doughnut-like density distribution (vanishing at the vortex core) and vortex flows (rapidly increasing at the vortex axis) of the BEC core of vortex soliton should directly affect the dynamical friction and, hence, star motion and their relaxation time. Since stellar motion is directly observed, the problem of the impact of

the vortex structure of ULDM on the dynamical friction provides the main motivation for the study performed here, which to the best of our knowledge was not considered in the literature.

The paper is organized as follows. Dynamical friction acting on moving stars in ULDM with trivial and nontrivial vortex structures is considered in a simplified approach in the absence of the ULDM self-interaction in Sec.2. Taking the ULDM self-interaction into account, the stellar orbital motion is studied in Sec.3. Conclusions are drawn in Sec.4.

## 2 Model and dynamical friction in simplified model

When a star moves in dark matter in the form of the BEC, it produces perturbation in galactic BEC leading to an overdense wake. This wake induces a drag force on a moving star resulting in dynamical friction. Conceptually, the situation is quite similar to the motion of bodies in fluids. Since the de Broglie wavelength of the BEC is much larger than stellar radii, we could approximate stars as point objects. As stars move through the galactic BEC, a gravitational wake forms behind them and the associated overdensity causes resistance producing dynamic friction acting on moving stars.

Since we are interested in the impact of vortex states of ULDM on the motion of stars and the repulsive self-interaction of ULDM is crucial for the existence of stable vortex structures, we employ in our analysis of dynamical friction the analytic results obtained in [27,28], which determined the coefficient of dynamical friction in self-interacting ULDM.

Let us begin our analysis by reminding the system of the Gross-Pitaevskii-Poisson equations [11,29] which govern the dynamical evolution of self-gravitating and self-interacting BEC field ψ and the gravitational potential Φ

$$i \hbar { \frac { } { } } \partial \psi = \left ( - \frac { \hbar { ^ { 2 } } } { 2 m } \nabla ^ { 2 } + g N | \psi | ^ { 2 } + m \Phi _ { g } \right ) \psi,$$

$$\nabla ^ { 2 } \Phi _ { g } = 4 \pi G m N | \psi | ^ { 2 },$$

where g = 4 π ℏ 2 a /m s is the coupling strength of the self-interaction of DM particles with mass m a , s is the s -wave scattering length, N is the number of boson particles, ℏ is the Planck constant, and G is the gravitational constant. The first term on the right-hand side of the Gross-Pitaevskii equation (1) corresponds to the kinetic energy. Repulsive self-interaction is described by the nonlinear term in Eq. (1) with g &gt; 0. The last term contains the self-consistent gravitational potential Φ which supports the formation of a localized BEC core.

As for the DM density in the BEC state, it was investigated for the case of the Milky Way galaxy in Refs. [30-33]. Approximate variational and numerical solutions to the GPP equations were found in [19,34]. In present work, we use a variational approach for stationary states with the following ansatz for ψ in cylindrical coordinates z, r, ϕ :

$$\psi ( t, r, \phi, z ) = A \left ( \frac { r } { R } \right ) ^ { s } e ^ { - i \mu t - \frac { r ^ { 2 } } { 2 R ^ { 2 } } - \frac { z ^ { 2 } } { 2 ( R _ { \eta } ) ^ { 2 } } + i s \phi },$$

where R and η are variational parameters, A is the normalization constant, and µ is the chemical potential. The quantum number s is the topological charge of the galactic BEC condensate. While the case s = 0 corresponds to the ground BEC state, the s = 1 state

describes the BEC soliton rotating around the center of the galaxy. States with higher topological charge s &gt; 1 are unstable [19,35]. The gravitational potential Φ can be found by solving the Poisson equation (2). This equation allows for an analytical solution in the s = 0 [34] state and Φ can be determined numerically for s &gt; 0.

It is very important for the study of dynamical friction that vortex states of BEC affect the DM profile differently. Here we use the DM profile ρ DM ∼ | ψ | 2 obtained in the case of the Milky Way galaxy in [34] for distances from the galaxy center less than 1 kpc for the trivial ground s = 0 and nontrivial s = 1 vortex states given by (for r &gt; 1 kpc, the Navarro-Frenk-White profile [36] was used)

$$\rho _ { s = 0 } = 0. 3 4 \cdot 1 0 ^ { - 1 7 } e ^ { - \left ( \frac { r _ { s p h } } { 0. 6 1 } \right ) ^ { 2 } } \frac { \text{kg} } { \text{m} ^ { 3 } }, \quad r < 1 \, k p c \\.$$

$$\rho _ { s = 1 } & = 0. 8 8 \cdot 1 0 ^ { - 1 7 } \left ( \frac { r } { 0. 4 } \right ) ^ { 2 } e ^ { - \left ( \frac { r } { 0. 4 } \right ) ^ { 2 } - \left ( \frac { z } { 0. 5 8 } \right ) ^ { 2 } } \frac { \ k g } { m ^ { 3 } }, \quad r < 1 \ k p c, \quad \cdots \quad \cdots \quad \cdots \quad \cdots$$

where r 2 sph = r 2 + z 2 is the radial spherical distance and r , z are coordinates in the cylindrical coordinate system measured in kpc, see Fig.1. This density function was obtained for the Milky Way halo mass M = 1 3 . × 10 12 M ⊙ , where M ⊙ is the solar mass, and radius R halo = 287 kpc . We use the model described in [37] and studied in [34] and choose the DM particle mass m = 2 92 . × 10 -22 eV/c 2 = 0 52 . × 10 -57 kg with the scattering length a s = 8 17 . × 10 -77 meter.

For the vortex state s = 1, the velocity of DM is given by [34]

$$\mathbf u = c \frac { r _ { u } } { r } \mathbf e _ { \phi },$$

where r u = 2 2 . · 10 -5 kpc. To determine dynamical friction, we use the rotation curve of the Milky Way galaxy observed by astrophysical means [38-40] in the interval 10 pc &lt; r &lt; 10 3 pc and discussed in detail in [40]. The first reason for considering distances larger 10 pc is connected with the presence of the supermassive black hole Sgr A* in the Milky Way center, which could severely affect the DM density profile for r ≤ 10 pc producing a dark matter spike with power law profile [41]. Another reason is that near the center of the galaxy, the contribution to the matter density from the baryonic component is significant [42]. Therefore, we consider larger distances where the impact of the supermassive black hole and baryonic component is not essential. We restrict our consideration to distances less 1 kpc where bosonic DM is in the BEC state. For larger distances, there is isothermal ULDM halo [37] where DM density is no longer given by Eqs.(4) and (5).

In this study, we quantify the effect of dynamic friction on moving stars. Since the self-interaction of ULDM particles is crucial for the existence of stable vortex states of BEC, we leave the analysis of dynamical friction of the s = 1 state for the next section. Here we will consider dynamical friction only for the ground s = 0 state which can be realized also in the fuzzy DM.

Although the DM profile of the s = 0 state in Eq.(4) was found in [34] in a model with ULDM self-interaction, we will assume that this profile is valid also in the case of fuzzy DM to compare the obtained results for dynamical friction with those found in the correct approach in the next section where self-interaction is taken into account.

In the case of fuzzy DM, dynamical friction was studied in [24]. For parameter β ≪ 1, where

$$\beta = \frac { 2 \pi G m _ { s t a r } m } { \hbar { v } _ { s t a r } },$$

Figure 1: Schematics of a stellar object in a circular orbit subject to effective dynamical friction caused by ULDM. (a) The left panel displays the density isosurface of ULDM solitonic core with zero angular momentum (state s = 0). The right panel depicts the Bose-Einstein condensate (BEC) core in the vortex state ( s = 1, the direction of circulation is shown by arrow). The dynamical friction force acting on a moving star, F fr , is shown by arrow. (b) Schematics of the condensate density distribution (shaded blue for the radial distribution and gray-scale density plot for the distribution in ( x, z )-plane) as a function of the radial coordinate r for the soliton (left panel) and vortex (right panel) states. The red dashed line shows the radial dependence of the ULDM velocity u in the vortex state.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00104v2-figure-a85e5769-ddac-49e2-9dbb-9bcfc7fda638.jpg)

the friction force equals [24]

$$F _ { f r } = 4 \pi \frac { G ^ { 2 } m _ { s t a r } ^ { 2 } \rho _ { D M } } { v _ { s t a r } ^ { 2 } } C ( \Lambda ),$$

where

$$\Lambda = \frac { k r } { \beta } = \frac { v _ { s t a r } ^ { 2 } r } { G m _ { s t a r } }$$

is the gravitational counterpart of the Coulomb logarithm , 1 m star is the star's mass, v star is the absolute value of the star's velocity, k = mv star / ℏ , and r is the radius of the stellar orbit. Since m star is much less than the host system in which it resides, we have the following simple relation for the function C [24]:

$$C ( \Lambda ) = \frac { 1 + \Lambda } { \Lambda } \arctan h \frac { \cdot \cdot \cdot } { 1 + \Lambda } - \sqrt { 1 + \frac { 2 } { \Lambda } }.$$

1 In plasma kinetic theory the Coulomb logarithm appears in the description of charged particle collisions with the long-range Coulomb interaction (clearly, a similar situation is realized in the case of the long-range Newtonian potential). The corresponding logarithmic divergence at large r can be eliminated by introducing a cut-off at small scattering angle. Physically, the elimination of divergences is explained by the fact that the Coulomb field of the charged particle at large distances is shielded by other charged particles [43,44].

In our case, parameter Λ ≫ 1 for distances 10 pc &lt; r &lt; 10 3 pc and it is convenient to use more simple asymptotic expression for the function C given by

$$C ( \Lambda ) = \ln ( 2 \Lambda ) ( 1 + 1 / \Lambda ) - 1.$$

Equation (7) implies that the friction force increases for stars with small velocities and in the region with high DM density, i.e., near the galactic center.

Our main quantity of interest is a characteristic time when a moving star changes notably its velocity [24]

$$T = \frac { v _ { s t a r } m _ { s t a r } } { F _ { f r } }.$$

Using Eq.(7), we find

$$T = \frac { v _ { s t a r } ^ { 3 } } { G ^ { 2 } m _ { s t a r } \rho _ { D M } } \frac { 1 } { 4 \pi C ( \Lambda ) }.$$

Our results for T in ULDM without self-interaction are presented by dash-dotted line in Fig.5, where we used m star = M ⊙ . As one can see, T monotonously grows with r . The reason is that while the DM density is almost constant at small distances 0.1 kpc &lt; r &lt; 0.34 kpc, v star increases linearly. Then at larger distances r &gt; 0 34 kpc stellar velocity . v star does not change much but DM density decreases notably. Hence T should grow in view of Eq.(12). In the next subsection, we will compare these results for the s = 0 state with those obtained in a more accurate approach where the ULDM self-interaction is taken into account.

## 3 Motion of stars in ULDM with vortex structure

To determine dynamical friction of stars moving in ULDM with different vortex structures, we assume that stars move on circular orbits with constant velocity and constant angular velocity Ω in the self-interacting ULDM and use the expression for dynamical friction derived in Ref. [28].

As we mentioned in the Introduction, a stable DM structure for the case of rotating BEC with the quantum number s = 1 is possible only for self-interacting ULDM particles with coupling constant of repulsive self-interaction exceeding some critical value which is fixed by the soliton core's observable radius and the DM particle mass [11,19]. Following Ref. [34], we take a s = 8 17 . · 10 -77 meter and m = 0 52 . · 10 -57 kg. Then we find g/ ℏ 2 = 1 97 . · 10 -18 meter/kg.

For our purposes, it is sufficient to consider dynamical friction in the direction tangential to a circular orbit. The dynamical friction force ⃗ F fr = -m star ⃗ ϕ ⃗, t ∇ ( r ) is determined by the gravitational potential ϕ induced by perturbation of the DM density by a moving body, for example, a star with mass m star . It can be found by solving the system of the Euler and continuity equations for the DM density and velocity as well as taking into account the Poisson equation for the gravitational potential (see, for details [24,45]). For a star with mass m star moving at constant velocity v star on a circular orbit of radius r , the tangential dynamical friction equals [28,45,46]

$$F _ { f r } = \frac { 4 \pi G ^ { 2 } m _ { s t a r } ^ { 2 } \rho _ { D M } } { c _ { s } ^ { 2 } } \mathcal { F },$$

where function F is given in the form of the sum over the magnetic quantum number m l and the azimuthal quantum number ℓ (note that the derivation of this dynamical friction force is quite similar to that given by Ostriker's formula [47] for the motion of a perturber in gaseous medium)

$$\mathcal { F } = \sum _ { \ell = 1 } ^ { \ell _ { \max } } \sum _ { m _ { l } = - \ell } ^ { \ell - 2 } \gamma _ { \ell m _ { l } } \text{Im} \left ( S _ { \ell, \ell - 1 } ^ { m _ { l } } - S _ { \ell, \ell - 1 } ^ { m _ { l } + 1 ^ { * } } \right ),$$

$$\gamma _ { \ell m _ { l } } & = ( - 1 ) ^ { m _ { l } } \frac { ( \ell - m _ { l } )! } { ( \ell - m _ { l } - 2 )! } \\ & \times \left \{ \Gamma \left ( \frac { 1 - \ell - m _ { l } } { 2 } \right ) \Gamma \left ( 1 + \frac { \ell - m _ { l } } { 2 } \right ) \Gamma \left ( \frac { 3 - \ell + m _ { l } } { 2 } \right ) \Gamma \left ( 1 + \frac { \ell + m _ { l } } { 2 } \right ) \right \} ^ { - 1 },$$

$$S _ { \ell, \ell - 1 } ^ { m _ { l } } = \frac { \pi i } { 2 \sqrt { 1 + \frac { m _ { f } ^ { \prime } } { \ell _ { q } ^ { \prime } } } } \left [ ( - 1 ) ^ { 1 + \theta ( m _ { l } ) } j _ { \ell } ( \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { - } ) j _ { \ell - 1 } ( \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { - } ) \\ + \, & i j _ { \ell } ( \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { - } ) y _ { \ell - 1 } ( \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { - } ) - j _ { \ell } ( i \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { + } ) h _ { \ell - 1 } ^ { ( 1 ) } ( i \ell _ { q } \mathcal { M } f _ { m _ { l } } ^ { + } ) \right ],$$

j l ( x ) and y l ( x ) are the spherical Bessel functions of the first and second kind, h (1) l ( x ) is the spherical Hankel function of the first kind, θ x ( ) is the Heaviside step-function, and

$$f _ { m _ { l } } ^ { \pm } = \sqrt { 2 } \left ( \sqrt { 1 + ( m _ { l } / \ell _ { \mathbf q } ) ^ { 2 } } \pm 1 \right ) ^ { 1 / 2 },$$

where parameter ℓ q is defined as

$$\mathbf e _ { q } = \frac { r m c _ { s } } { \hbar { \mathcal { M } } }.$$

The important characteristic of the problem of dynamical friction is the Mach number

$$\mathcal { M } & = \frac { v } { c _ { s } },$$

where v is the absolute value of the relative velocity of the star and ULDM velocity u and c s is the adiabatic sound speed of DM superfluid. The latter is given by the first derivative of pressure concerning density [29]

$$c _ { s } = \frac { \cdot \cdot \cdot } { m }.$$

Note that pressure is proportional to the coupling constant of self-interaction g , therefore, c s vanishes in fuzzy dark matter.

For the ground BEC state s = 0, where u = 0, the relative velocity v coincides with the star's velocity v = v star . In the vortex state s = 1, there are two possiblities of counter-rotating ( v = v star + u ) and corotating star motion ( v = | v star -u | ). Since the BEC velocity rapidly decreases with distance from the galactic center, the difference between these two cases is significant only at a small distance from the galactic center.

Figure 2: The Mach number as a function of distance r to the galactic center for different vortex states of ULDM. We consider the ground BEC state s = 0 (solid black line) and vortex BEC state s = 1 for the cases where a probe star corotates (dotted red line) and counter-rotates (dashed green line) with respect to DM superfluid. The horizontal black solid line corresponds to the Mach number M = 1.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00104v2-figure-e29c32b3-beab-4325-aa44-d1d0b18b1a8e.jpg)

The Mach number as a function of distance to the galactic center r is plotted in Fig.2 for the ground ( s = 0) and vortex ( s = 1) for the cases where a probe star corotates and counter-rotates with respect to DM superfluid.

Let us discuss now parameter ℓ q . This parameter as a function of distance r to the center of the galaxy is presented in Fig.3 for the ground s = 0 BEC state and s = 1 vortex BEC state for corotating and counter-rotating stellar orbits. As one can see, parameter ℓ q is typically much larger than the unity. Further, the maximum azimuthal quantum number ℓ max up to which the sum over l is performed is determined by the size R probe of probe body and orbit's radius r

$$\text{is $r$} \\ \ell _ { \max } = \frac { \pi r } { R _ { \text{probe} } } \,.$$

Since the radius of a star is usually much less than its orbit's radius, we have ℓ max ≫ 1. Since series (14) is convergent, we could introduce instead of ℓ max a cut-off whose value is chosen by the required accuracy of calculations.

Figure 3: Parameter ℓ q as a function of distance r to the galactic center for the s = 0 ground BEC state (solid black line) and the cases of counter-rotating (dashed green line) and corotating (dotted red line) star's orbits for the s = 1 vortex BEC state.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00104v2-figure-c88395b2-cf59-47f7-9916-dbfbc8905a73.jpg)

It is worth mentioning that there is a simple analytic expression which approximates the tangential component of the friction force for M &lt; 1 and ℓ q →∞ [28]

$$F _ { f r } = \frac { 4 \pi G ^ { 2 } m _ { s t a r } \rho _ { D M } } { c _ { s } ^ { 2 } } \frac { \arctanh \mathcal { M } - \mathcal { M } } { M ^ { 2 } }.$$

Comparing this relation with (13) one can see that function F in this case has a simple form

$$\mathcal { F } = \frac { \arctanh \mathcal { M } - \mathcal { M } } { M ^ { 2 } }.$$

As one can see from Fig.2, there is a region where the Mach number is less the unity, i.e., M &lt; 1, but parameter ℓ q , in view of Fig.3, is large in this region but not infinite. This interval defines the region where the above approximate formula for the dynamical friction force could be applied.

For stars in their rotation motion in the Milky Way, our results for function F defined by Eq.(14) and the approximate analytic formula (23) are given in Fig.4. It can be seen that the results agree quite well for the region where the Mach number is M ≲ 1, otherwise, the friction force defined by the approximate expression (22) blows up.

Having defined the dynamical friction force for ULDM with self-interaction, we can easily determine the characteristic time for a star to decrease its velocity using Eq.(11). The corresponding results are plotted by solid, dashed, and dotted lines in Fig.5 for distances r &gt; 1 pc in the z = 0 plane and and can be summarized as follows. (i) The characteristic time needed for a star to decrease its velocity is practically the same for the ground s = 0 and vortex s = 1 states for r ≳ 0 3 kpc, where the ULDM velocity . u is small and can be neglected. (ii) For r ≲ 50 pc, T is approximately the same for the counter-rotating and corotating motion in the s = 1 state for which the star's velocity can be neglected compared to the ULDM velocity at small r . Clearly, the decrease of ULDM density in the s = 1 state compared to the s = 0 state, see Fig.1, strongly diminishes the friction force and makes T very large at a small distance. (iii) T strongly increases for the co-rotating case in the s = 1 state at r ≈ 0 06 kpc, where the relative velocity of . the star's velocity and the ULDM velocity vanishes. This means that the ULDM rotation makes stellar orbital motion more stable at a certain distance from the galactic center.

Figure 4: The function F as a function of distance r to the galactic center for the s = 0 ground BEC state (panel a ) and corotating (panel b ) and counter-rotating (panel c ) stellar orbits for the s = 1 vortex BEC state. Solid black lines correspond to the function F given by Eq.(14) and dashed red lines correspond to F given by the approximate formula (23). Two vertical solid lines restrict the region where M &lt; 1 and the approximate analytic formula (23) can be applied.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00104v2-figure-cee6b71d-88a0-4a6d-860d-c60839535013.jpg)

Finally, let us compare the results for the ground s = 0 state found in the simplified approach using Eq.(7) for the dynamical friction force where the ULDM self-interaction was not taken into account (dash-dotted black line in Fig.5) with the results obtained in the case where the ULDM self-interaction is taken into account (solid black line). Clearly, we see that the ULDM self-interaction is very important for the relaxation time.

Figure 5: The characteristic time in years needed to decrease the stellar velocity in selfinteracting DM as a function of distance to the galactic center for the DM particle mass m = 2 92 . · 10 -22 eV/c 2 for the s = 0 ground BEC state (solid black line) and the cases of counter-rotating (dashed green line) and corotating (dotted red line) stellar orbits for the s = 1 vortex BEC state. The dash-dotted black line corresponds to fuzzy DM without self-interaction plotted using Eq.(12).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00104v2-figure-56b37d3f-1acf-4a35-aec4-b1563aa23f21.jpg)

## 4 Conclusions

For ULDM in the state of BEC, we considered two different configurations of BEC: ground (not rotating) and vortex states. We studied the dynamical friction force acting on stars moving in ultralight dark matter with different soliton vortex structures and investigated the characteristic time T for a star to decrease its velocity. The main results of our investigation are presented in Fig.5, where the characteristic relaxation time for stars in their rotation in the Milky Way is displayed and are shown by solid black ( s = 0 state), dashed green (the counter-rotating case of s = 1 vortex DM state) and dotted red (the corotating case of s = 1 vortex DM state) lines. According to this figure, the relaxation time always exceeds the Hubble time.

We found that for large distances r ≳ 0 3 . kpc, the characteristic time needed to decrease the stellar velocity is practically the same for the ground s = 0 and vortex s = 1 states because the ULDM velocity is small. At small distances, r ≲ 50 pc, T is approximately the same for the counter-rotating and corotating stellar motion in the s = 1 state for which the star's velocity can be neglected compared to the ULDM velocity at small r . T strongly increases for the corotating case in the s = 1 state at r ≈ 0 06 . kpc, where the relative velocity of the star's velocity and the ULDM velocity vanishes. This means the ULDM rotation makes the stellar orbital motion more stable at a certain distance from the galactic center.

Overall, our analysis shows that the orbital motion of stars of the solar mass is stable in the Milky Way galaxy at least for distances from the galactic center r &gt; 10 pc and

the dynamic friction induced by BEC is not sufficient to affect strongly star's motion. Still, we found that the BEC structure and ULDM self-interaction are very important for the analysis of dynamical friction on star's motion in galaxies and much work remains to be done in this direction especially for objects much more massive compared to the Sun because the acceleration due to the friction force grows with mass.

In addition, we confirmed the conclusion made in previous studies that the ULDM selfinteraction essentially affects the dynamical friction. A similar conclusion was reached in recent studies [48-51] of the dynamical friction force for inspiraling black holes in the case of cold dark matter, where self-interaction may solve the final parsec problem of supermassive black holes merger.

## Acknowledgements

The authors thank K. Korshynska for useful discussions. V.M.G. is grateful to Wuhan University for its hospitality during his stay there. A.I.Y. acknowledge support from the projects 'Ultracold atoms in curved geometries', 'Theoretical analysis of quantum atomic mixtures' of the University of Padova, and from INFN.

## References

- [1] G. Bertone, D. Hooper and J. Silk, Particle dark matter: Evidence, candidates and constraints , Phys. Rept. 405 (2005) 279 [ hep-ph/0404175 ].
- [2] P. Salucci, The distribution of dark matter in galaxies , Astron. Astrophys. Rev. 27 (2019) 2 [ 1811.08843 ].
- [3] M. Y. Khlopov, B. A. Malomed, I. B. Zeldovich and Y. B. Zeldovich, Gravitational instability of scalar fields and formation of primordial black holes , Mon. Not. Roy. Astron. Soc. 215 (1985) 575.
- [4] B. Carr, F. Kuhnel and M. Sandstad, Primordial Black Holes as Dark Matter , Phys. Rev. D 94 (2016) 083504 [ 1607.06077 ].
- [5] G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, M. Pierre et al., The waning of the WIMP? A review of models, searches, and constraints , Eur. Phys. J. C 78 (2018) 203 [ 1703.07364 ].
- [6] T. Asaka and M. Shaposhnikov, The ν MSM, dark matter and baryon asymmetry of the universe , Physics Letters B 620 (2005) 17-26 [ hep-ph/0505013 ].
- [7] T. Asaka, S. Blanchet and M. Shaposhnikov, The nuMSM, dark matter and neutrino masses , Phys. Lett. B 631 (2005) 151 [ hep-ph/0503065 ].
- [8] C. B. Adams et al., Axion Dark Matter , in Snowmass 2021 , 3, 2022, 2203.14923 .
- [9] G. Bertone, Particle Dark Matter: Observations, Models and Searches . Cambridge University Press, 2010.
- [10] E. G. M. Ferreira, Ultra-light dark matter , Astron. Astrophys. Rev. 29 (2021) 7 [ 2005.03254 ].

- [11] X. Calmet, ed., Quantum Aspects of Black Holes , vol. 178 of Fundamental Theories of Physics . Springer, 2015, 10.1007/978-3-319-10852-0.
- [12] T. Matos, L. A. Ure˜a L´pez and J.-W. Lee, n o Short review of the main achievements of the scalar field, fuzzy, ultralight, wave, BEC dark matter model , Front. Astron. Space Sci. 11 (2024) 1347518 [ 2312.00254 ].
- [13] A. Amruth, T. Broadhurst, J. Lim, M. Oguri, G. Smoot, J. Diego et al., Einstein rings modulated by wavelike dark matter from anomalies in gravitationally lensed images , Nature Astronomy 7 (2023) 736.
- [14] W. Hu, R. Barkana and A. Gruzinov, Cold and fuzzy dark matter , Phys. Rev. Lett. 85 (2000) 1158 [ astro-ph/0003365 ].
- [15] S.-J. Sin, Late time cosmological phase transition and galactic halo as Bose liquid , Phys. Rev. D 50 (1994) 3650 [ hep-ph/9205208 ].
- [16] H.-Y. Yang, B. T. Chiang, G.-M. Su, H.-Y. Schive, T. Chiueh and J. P. Ostriker, Galactic disc heating by density granulation in fuzzy dark matter simulations , Mon. Not. Roy. Astron. Soc. 530 (2024) 129 [ 2403.09845 ].
- [17] J.-w. Lee and I.-g. Koh, Galactic halos as boson stars , Phys. Rev. D 53 (1996) 2236 [ hep-ph/9507385 ].
- [18] T. Matos and F. S. Guzman, Scalar fields as dark matter in spiral galaxies , Class. Quant. Grav. 17 (2000) L9 [ gr-qc/9810028 ].
- [19] Y. O. Nikolaieva, A. O. Olashyn, Y. I. Kuriatnikov, S. I. Vilchynskii and A. I. Yakimenko, Stable vortex in Bose-Einstein condensate dark matter , Low Temp. Phys. 47 (2021) 684 [ 2103.07856 ].
- [20] Y. O. Nikolaieva, Y. M. Bidasyuk, K. Korshynska, E. V. Gorbar, J. Jia and A. I. Yakimenko, Stable vortex structures in colliding self-gravitating bose-einstein condensates , Phys. Rev. D 108 (2023) 023503.
- [21] R. Boey, Y. Wang, E. Kendall and R. Easther, Dynamical friction and black holes in ultralight dark matter solitons , Phys. Rev. D 109 (2024) 103526 [ 2403.09038 ].
- [22] B. Bar-Or, J.-B. Fouvry and S. Tremaine, Relaxation in a fuzzy dark matter halo , Astrophysical Journal 871 (2019) 28.
- [23] S. Chandrasekhar, Dynamical Friction. I. General Considerations: the Coefficient of Dynamical Friction , Astrophys. J. 97 (1943) 255.
- [24] L. Hui, J. P. Ostriker, S. Tremaine and E. Witten, Ultralight scalars as cosmological dark matter , Phys. Rev. D 95 (2017) 043541 [ 1610.08297 ].
- [25] L. Lancaster, C. Giovanetti, P. Mocz, Y. Kahn, M. Lisanti and D. N. Spergel, Dynamical Friction in a Fuzzy Dark Matter Universe , JCAP 01 (2020) 001 [ 1909.06381 ].

- [26] Y. Wang and R. Easther, Dynamical friction from ultralight dark matter , Phys. Rev. D 105 (2022) 063523.
- [27] N. Glennon, N. Musoke, E. O. Nadler, C. Prescod-Weinstein and R. H. Wechsler, Dynamical friction in self-interacting ultralight dark matter , Phys. Rev. D 109 (2024) 063501.
- [28] L. Berezhiani, G. Cintia, V. De Luca and J. Khoury, Dynamical friction in dark matter superfluids: The evolution of black hole binaries , JCAP 06 (2024) 024 [ 2311.07672 ].
- [29] E. G. M. Ferreira, Ultra-light dark matter , Astron. Astrophys. Rev. 29 (2021) 7 [ 2005.03254 ].
- [30] G. L. Bryan and M. L. Norman, Statistical properties of x-ray clusters: Analytic and numerical comparisons , Astrophys. J. 495 (1998) 80 [ astro-ph/9710107 ].
- [31] H.-Y. Schive, M.-H. Liao, T.-P. Woo, S.-K. Wong, T. Chiueh, T. Broadhurst et al., Understanding the Core-Halo Relation of Quantum Wave Dark Matter from 3D Simulations , Phys. Rev. Lett. 113 (2014) 261302 [ 1407.7762 ].
- [32] H.-N. Lin and X. Li, The Dark Matter Profiles in the Milky Way , Mon. Not. Roy. Astron. Soc. 487 (2019) 5679 [ 1906.08419 ].
- [33] A. Pozo et al., Galaxy formation with Wave/Fuzzy Dark Matter: The core-halo structure , 2310.12217 .
- [34] K. Korshynska, Y. M. Bidasyuk, E. V. Gorbar, J. Jia and A. I. Yakimenko, Dynamical galactic effects induced by solitonic vortex structure in bosonic dark matter , Eur. Phys. J. C 83 (2023) 451 [ 2301.13110 ].
- [35] A. S. Dmitriev, D. G. Levkov, A. G. Panin, E. K. Pushnaya and I. I. Tkachev, Instability of rotating Bose stars , Phys. Rev. D 104 (2021) 023504 [ 2104.00962 ].
- [36] J. F. Navarro, C. S. Frenk and S. D. M. White, The Structure of cold dark matter halos , Astrophys. J. 462 (1996) 563 [ astro-ph/9508025 ].
- [37] P.-H. Chavanis, Predictive model of BEC dark matter halos with a solitonic core and an isothermal atmosphere , Phys. Rev. D 100 (2019) 083022 [ 1810.08948 ].
- [38] R. J. A. Lambourne, S. Serjeant and M. H. Jones, eds., An Introduction to Galaxies and Cosmology . Cambridge University Press, 2015.
- [39] Y. Sofue, M. Honma and T. Omodaka, Unified Rotation Curve of the Galaxy Decomposition into de Vaucouleurs Bulge, Disk, Dark Halo, and the 9-kpc Rotation Dip -, Publ. Astron. Soc. Jap. 61 (2009) 227 [ 0811.0859 ].
- [40] Y. Sofue, Face-on map of the molecular disc and 3-kpc expanding ring of the galaxy based on a high-accuracy rotation curve , Astrophys. Space Sci. 368 (2023) 74 [ 2309.01135 ].

- [41] S. Nampalliwar, K. Saurabh, K. Jusufi, Q. Wu, M. Jamil and P. Salucci, Modelling the Sgr A* Black Hole Immersed in a Dark Matter Spike , The Astrophysical Journal 916 (2021) 116.
- [42] K. Jusufi, M. Jamil and T. Zhu, Shadows of Sgr A ∗ black hole surrounded by superfluid dark matter halo , Eur. Phys. J. C 80 (2020) 354 [ 2005.05299 ].
- [43] J. Krommes, An introduction to the physics of the Coulomb logarithm, with emphasis on quantum-mechanical effects , J. Plasma Phys. 85 (2019) 925850101 [ 1806.04990 ].
- [44] L. Pitaevskii and E. Lifshitz, Physical Kinetics: Volume 10 . Elsevier Science, 2012.
- [45] V. Desjacques, A. Nusser and R. Buehler, Analytic solution to the dynamical friction acting on circularly moving perturbers , The Astrophysical Journal 928 (2022) 64.
- [46] R. Buehler and V. Desjacques, Dynamical friction in fuzzy dark matter: Circular orbits , Phys. Rev. D 107 (2023) 023516 [ 2207.13740 ].
- [47] E. C. Ostriker, Dynamical friction in a gaseous medium , Astrophys. J. 513 (1999) 252 [ astro-ph/9810324 ].
- [48] G. Alonso- ´ lvarez, J. M. Cline and C. Dewar, A Self-Interacting Dark Matter Solves the Final Parsec Problem of Supermassive Black Hole Mergers , Phys. Rev. Lett. 133 (2024) 021401 [ 2401.14450 ].
- [49] M. S. Fischer and L. Sagunski, Dynamical friction from self-interacting dark matter , Astron. Astrophys. 690 (2024) A299 [ 2405.19392 ].
- [50] H. Koo, D. Bak, I. Park, S. E. Hong and J.-W. Lee, Final parsec problem of black hole mergers and ultralight dark matter , Phys. Lett. B 856 (2024) 138908 [ 2311.03412 ].
- [51] B. C. Bromley, P. Sandick and B. Shams Es Haghi, Supermassive black hole binaries in ultralight dark matter , Phys. Rev. D 110 (2024) 023517 [ 2311.18013 ].
V. M. Gorkavenko O. V. Barabash T. V. Gorkavenko O. M. Teslyk A. O. Zaporozhchenko Junji Jia A. I. Yakimenko E. V. Gorbar
Dynamical friction in rotating ultralight dark matter galactic cores
Dynamical friction and stellar orbital motion in spiral galaxies with dark
matter composed of ultralight bosons in the state of rotating Bose-Einstein
condensate (BEC) are studied. It is found that the dynamical friction force is
significantly affected by the topological charge of the vortex structure of the
BEC core with the strongest effect at distances near the galactic center. It is
also shown that the ultralight dark matter self-interaction plays an important
role in studying the dynamical friction.