1Key Laboratory of Hunan Province on Information Photonics and Freespace Optical Communications, School of Physics and Electronic Science, Hunan Institute of Science and Technology, Yueyang 414006, China
2Graduate School of China Academy of Engineering Physics, Beijing 100193, China
We aim to find one highly nontrivial example of the solutions to the vortex fluid dynamical equation on the unit sphere (S2) and compare it with the numerical simulation. Since the rigid rotating steady solution for vortex fluids on S2 is already known to us, we consider the perturbations above it. After decomposing the perturbation of the vortex number density and vortex charge density into spherical harmonics, we find that the perturbations are propagating waves. To be precise, the velocities for different single-mode vortex number density waves are all the same, while the velocities for single-mode vortex charge density waves depend on the degree of the spherical harmonics l, which is a signal of the existence of dispersion. Meanwhile, we find that there is a beat phenomenon for the positive (or negative) vortex density wave. Numerical simulation based on the canonical equations for the point vortex model agrees perfectly with our theoretical calculations.
Yanqi Xiong, Zhijun Zou, Liang Luo. Density waves of vortex fluids on a sphere[J]. Communications in Theoretical Physics, 2025, 77(4): 045702. DOI: 10.1088/1572-9494/ad7e97
1. Introduction
Fluids on curved surfaces exhibit abundant phenomena that are absent in planar systems. Significant studies have been conducted to explore the interplay between geometry, topology, and fluid dynamics in various systems, including active matter [1–3], quantum Hall liquids [4–9], and classical fluids [10–13].
Due to the presence of couplings to curvature, quantum vortices have a fundamental role in determining the properties of superfluids on curved surfaces [14, 15]. In gravity environments, superfluids on curved surfaces tend to settle at the bottom of the surface. However, recent advances in experiments with Bose–Einstein condensates on the International Space Station [16] have given rise to a promising possibility to investigate bubble-trapped superfluids in ultracold atomic bubbles [17]. Intrigued by this experimental progress, the research interest in few-body vortex dynamics on curved surfaces has been revitalized [18–20].
To investigate the dynamics of systems comprised of numerous quantum vortices, a vortex fluid [21, 22] is a suitable model. Prior studies have uncovered novel properties of vortex fluids on planar surfaces. For instance, binary vortex fluids are compressible [22] while chiral vortex fluids are incompressible [21]; Additionally, there exists an odd viscous tensor, and the circulation quantum acts as a non-dissipative coefficient for this odd viscosity; Vortex fluids are closely linked to quantum Hall liquids [23] and fractons [24, 25].
Recently, researchers have begun to explore the effects of curvature and topology on the collective dynamics of quantum vortices. Xiong and Yu have obtained a steady-state solution of a rigid rotation of vortex fluids on a spherical surface [26]. In this paper, we will mainly focus on the perturbation properties of this steady-state solution. We consider the problem in the dynamical regime where the number of point vortices is sufficiently large, and the evolution time is not too long such that non-equilibrium processes induced by stochasticity can be neglected. Our theoretical results demonstrate that despite the background being a steady-state flow of rigid rotation, perturbations imposed on it have their own propagation velocities. Moreover, the propagation velocities of the vortex number density wave and the vortex charge density wave differ. This phenomenon cannot be described by continuous vorticity models [27], as such models only have the vorticity ω that corresponds to the vortex charge density and lack a corresponding quantity for the vortex number density. Additionally, our numerical simulation results match well with the theory in a quantitative manner.
2. Vortex fluids on S2 and rigid rotating steady solutions
2.1. Vortex fluids on a sphere
2.1.1. Point vortex model
Generally, Bose–Einstein condensates (BECs) with interaction are a type of superfluid characterized by a macroscopic wave function $\Psi$ that satisfies the Gross-Pitaevski (G-P) equation [28]. The single valuedness of $\Psi$ requires that vortex excitations in BECs have quantized circulations:
Here, m is the mass of the atoms, $\vec{u}$ is the velocity field inside the superfluid, n is an integer called winding number (in this article, only ±1 are considered), and h is the Planck constant.
In many cases, the distance between vortices in a superfluid is much larger than their core size (the healing length [29]). In these cases, the vortex system can be treated as a system of point vortices, whose vorticity is given by a sum of delta functions
where θ and φ are the polar and azimuthal angles, respectively. However, the counter terms proportional to $-\tfrac{1}{4\pi }$ in the definition of vorticity are unphysical in real superfluids since there is no continuous vorticity in superfluids. The only way out is to ask the total circulation to be equal to 0, i.e. ∑iΓi = 0. In that case, all the counter terms cancel neatly among themselves.
2.1.2. Vortex fluids on a sphere
In this article, we only consider two types of point vortices, for which the winding number n = ± 1, namely the positive vortex and negative vortex, or vortex and anti-vortex for short. After coarse graining [22], the system is characterized by the vortex number density and vortex charge density
where $({\theta }_{i}^{+},{\phi }_{i}^{+})$ is the coordinate of the ith positive vortex and $({\theta }_{i}^{-},{\phi }_{i}^{-})$ is the coordinate of the ith negative vortex. Meanwhile, the number of positive vortices N+ must be equal to that of negative vortices N−, to ensure the cancelation of total vorticity. By continuity relation the vortex number current and vortex charge current are induced:
where ${v}_{i}^{\mu }$ is the velocity of the ith point vortex and its expression is given by Kimura [30], and the coarse-grained vortex number velocity field vμ and vortex charge velocity field wμ for vortex fluids are defined. The continuity equations read
where $\gamma =\tfrac{{\hslash }}{m}$ and N is the total number of vortices. Obviously, the system is rotating about the axis like a rigid body.
3. Density waves on the sphere
3.1. Theory
In this section, we will add perturbations to the rigid rotating background introduced at the end of the last section. To this purpose we write every physical quantity as a sum of two parts, namely the unperturbed part and the perturbation part
Besides, using spherical harmonics as the expansion basis for the perturbation part has the following advantages: The single-mode spherical harmonic perturbation above the rigid rotating background evolves in time as a traveling wave, while other types of perturbation generally do not have such a simple temporal behavior because of dispersion, which we will shortly see. It is worth mentioning that Shankar et al [2] also used the method of linear perturbation about the steady-state solution to explore the dynamic mode in an article on active matters on curved surfaces, although their work did not involve spherical harmonics. The absence of the Y00 component for δσ is the direct result of the requirement that the total vorticity must be equal to 0. A lengthy calculation (see appendix) leads to the following result
Obviously, this is a wave solution, and the waves are propagating along the equator of the unit sphere. What’s more, the wave velocities for δρlm and δσlm are different and the latter depends on l generally. In the large N limit, the wave velocity for δρlm agrees with the rotating velocity of the rigid background, while the wave velocity for δσlm does not. Therefore, we have found that there exists dispersion for propagating vortex charge density waves even in dynamical regimes.
There is another important result, namely the beat phenomenon. If we start from the following initial condition
The above equations clearly show that the amplitudes of the density waves of vortex and anti-vortex themselves oscillate with a frequency of $\tfrac{{mN}\gamma }{4l(l+1)}$ separately, and they differ from each other by a phase difference of π/2. This is indeed a beat phenomenon.
3.2. Numerical results
To check our theory, we use numerical methods to simulate the dynamical evolution of point vortex systems on the unit sphere based on the canonical equations [30]. Specifically, the fourth-order Runge-Kutta method is employed to solve the following equations for N vortices
Figure 1 shows the propagation of the vortex number density waves and vortex charge density waves. To obtain a satisfying collective behavior, the total number of vortices should not be too small, since our theory works best in the dynamical regime. When N is relatively small, the stochasticity will play an important role and we should develop a kinetic theory [31]. Meanwhile, to the rigid rotating background, we have to add perturbations, whose stochastic fluctuations also should not be too significant. Through several trials, we found that N = 80000 is a decent choice. We have chosen $\delta \rho =\delta \sigma \propto {\sin }^{2}\theta \cos 2\phi $ (proportional to Y22 + Y2,−2) and $\delta \rho =\delta \sigma \propto {\sin }^{3}\theta \cos 3\phi $ (proportional to Y33 + Y3,−3) instead of $\delta \rho =\delta \sigma \propto \sin \theta \cos \phi $ (proportional to Y11 + Y1,−1) as our examples for the following reasons. When θ → π, $\sin \theta \cos \phi \sim (\pi -\theta )$ . If we have chosen $\delta \rho =\delta \sigma \propto \sin \theta \cos \phi $ , then $\delta {\rho }_{+}\equiv \tfrac{1}{2}(\delta \rho +\delta \sigma )\,\propto \sin \theta \cos \phi \sim (\pi -\theta )$ . Remembering that for our background ${\sigma }_{0}={\rho }_{0}\cos \theta $ , therefore the unperturbed density of positive vortices ${\rho }_{0}^{+}=\tfrac{1}{2}{\rho }_{0}(1+\cos \theta )\sim {\left(\pi -\theta \right)}^{2}$ , which can be much smaller than (π − θ). Thus in some areas δρ+ is even greater than ${\rho }_{0}^{+}$ , which will make the perturbation theory invalid. What is even worse is that ρ0 + δρ < 0 in certain areas, which is absolutely not allowed. Figures 1(a) (b) show the propagation of the vortex number density wave and vortex charge density wave, respectively, with the initial perturbation (at t = t0) sampled according to the distribution $\delta \rho =\delta \sigma \propto {\sin }^{2}\theta \cos 2\phi $ which is the real part of Y22(θ, φ) up to a coefficient. Similarly, figures 1(c) (d) are according to the distribution $\delta \rho =\delta \sigma \propto {\sin }^{3}\theta \cos 3\phi $ which is the real part of Y33(θ, φ) up to a coefficient. t0 through t5 are six evenly-spaced time points. The φ axis is evenly partitioned into 40 intervals and the vortex number and vortex charge within each interval are counted at t = t0, resulting in the blue-colored histograms. The coordinates of positive and negative vortices at t1, t2, t3, t4 and t5 are calculated according to the canonical equations for point vortex systems on the unit sphere given the aforementioned initial perturbations. After the above procedure, sinusoidal curve fittings are used, resulting in the six colored curves in each subfigure. The theory and numerical results match very well. Through the root-finding algorithm, the displacements Δφs from t0 to t5 are found to be 0.398, 0.269, 0.405 and 0.318 for figures 1(a), (b), (c) and (d), respectively. The theoretical predictions are 2/5, 4/15, 2/5 and 1/3, respectively.
Figures 2(a)–(e) show the propagation of the density wave of positive vortices with the initial perturbation(at t = t0) sampled according to the distribution $\delta {\rho }^{+}\propto {\sin }^{2}\theta \cos 2\phi $ and δρ− = 0. Similarly, figures 2(f)–(j) show the propagation of the density wave of negative vortices in the same process. The total number of vortices N = 30000, with the numbers of positive and negative vortices being equal. t0 through t4 are five evenly-spaced time points. The φ axis is evenly partitioned into 40 intervals and the positive and negative vortex numbers within each interval are counted at t = t0 s. The coordinates of positive vortices and negative vortices at t1, t2, t3 and t4 are calculated according to the canonical equations for point vortex systems on the unit sphere given the aforementioned initial perturbations. Through the above procedure, 10 histograms are produced and then sinusoidal curve fittings are used, resulting in 10 colored curves. We have carefully chosen t1 such that t1 − t0 is one-eighth of the oscillation period of the amplitudes obtained from equations (44), (45) so that we expect the amplitude of the positive vortex density wave to reduce to 0 at t2. The theory and numerical results match very well. δρ+ and δρ− go to 0 at t2 and t4, respectively.
Figure 2. Propagating vortex density waves with beat phenomenon.
4. Conclusions
In this paper, we have investigated the behavior of the perturbations above the axisymmetric rigid rotating state. It is found that the perturbations propagate generally. The propagation velocities for single-mode vortex number density waves are constantly Nγ/4. However, the propagation velocities for single-mode vortex charge density waves depend on the degree of the spherical harmonics l. Thus a general wavelet of perturbation will disperse while propagating. This is a bit counterintuitive since one may naively anticipate that the rigid rotating background will drive everything attached to it to move with the same constant speed as itself.
It is worth noting that we have obtained the velocities for vortex number density waves. If we consider the continuous vortex rather than the point vortex system on S2, then it’s nonsense to talk about the vortex number density waves, since in that case we only have one density, namely the vortex charge density σ or vorticity ω to describe our system. Therefore, the motion of a vortex fluid is different from the motion of the classical vortex of an ideal fluid. The former is a collective description of a large number of point vortices based on the binary point vortex model. It is an emergent phenomenon after taking the positive and negative point vortices as elementary constituents, so it contains more physics than the classical vortex. Our study also provides an example of the interaction between topological defects and curvature, which can be seen when we insert the radius R of the sphere in equations (36), (37) according to dimensional requirement, and is expected to be useful for investigating rich phenomena involving a large number of quantum vortices [32–34] in bubble-trapped Bose–Einstein condensates [35, 36].
Meanwhile, we find that the amplitudes of the positive and negative vortex density waves oscillate themselves. This phenomenon is a mathematical consequence of the fact that the propagation velocity of the vortex number density wave is different from that of the vortex charge density wave.
Acknowledgments
This work was supported by the Scientific research projects of Hunan Provincial Department of Education (Grant Nos. 22A0477 and 20B273). We also acknowledge L Dai for useful discussions.
Thus once we know the vorticity ω or charge density σ, by solving the Poisson equation we can obtain the stream function $\Psi$ and the superfluid velocity field u. First, we need to find the Green function for S2 which satisfies the following equation [30],
GromovA, ChoG Y, YouY, AbanovA G, FradkinE2015 Framing anomaly in the effective theory of the fractional quantum Hall effect Phys. Rev. Lett.114 016805