Tkachenko Mode in Rotating Superfluids

• Tkachenko Mode is the low-frequency transverse shear oscillation of the triangular vortex lattice in rotating superfluids, exhibiting soft, elastic behavior.
• The mode features elliptically polarized vortex trajectories with a dispersion relation that transitions from linear (incompressible limit) to quadratic (compressible or rapid rotation regimes).
• Experimental and theoretical studies use imaging, Bragg spectroscopy, and BdG analyses in Bose–Einstein condensates to explore its role in vortex lattice dynamics.

Tkachenko modes are the low-frequency transverse shear oscillations of the vortex lattice that forms in a rotating superfluid or Bose–Einstein condensate. When rotation generates a triangular array of quantized vortices, small collective distortions of that array support an unusually soft elastic wave in which vortex cores move around equilibrium sites with predominantly transverse, often slightly elliptical, polarization and only weak density modulation at low frequency. In an incompressible homogeneous fluid the mode is linear, $\omega=c_T k$, whereas compressibility and rapid rotation produce the characteristic long-wavelength quadratic behavior; in the lowest-Landau-level regime this quadratic law is the leading long-wavelength result (Sonin, 2013, Jeevanesan et al., 2022).

1. Physical character and distinction from related modes

When a neutral superfluid is set into steady rotation with angular velocity $\Omega$, it nucleates an ordered array of quantized vortices with areal density fixed by the Feynman relation $n_v=2\Omega/\kappa$, $\kappa=h/m$. The energetically preferred arrangement is triangular. Tkachenko waves are the long-wavelength shear deformations of this vortex crystal: the restoring force comes from vortex-lattice elasticity, while the inertia and coupling to density fluctuations come from the surrounding superfluid (Sonin, 2013).

These modes are distinct from ordinary sound and from Kelvin waves. Sound modes are longitudinal density waves with $\omega \simeq c_s k$, whereas Tkachenko modes are dominantly transverse lattice-shear oscillations with little density fluctuation at low $\omega$. Kelvin waves are helical bending excitations along individual vortex lines and do not require lattice shear. In trapped three-dimensional condensates the two can hybridize, producing Kelvin–Tkachenko branches whose even-parity members are largely planar and whose odd-parity members carry gyroscopic tilt of vortex lines (Simula et al., 2010).

The polarization is not exactly circular. In the hydrodynamic and trapped-BdG descriptions, vortices typically execute elliptic trajectories about equilibrium positions. In the incompressible limit the motion is nearly purely transverse; in compressible systems a small longitudinal component is induced, and in rapid-rotation or lowest-Landau-level settings time-reversal breaking fixes a definite chirality of the elliptical motion (Jeevanesan et al., 2022).

2. Continuum elasticity, compressibility, and dispersion

A coarse-grained description uses the vortex displacement field $u(x,y,t)$, superfluid velocity $v(x,y,t)$, and the elastic moduli of the triangular lattice. In an incompressible homogeneous fluid the transverse displacement obeys a wave equation with shear-wave speed

$c_T=\sqrt{\frac{C_{66}}{\rho}},$

so that

$\omega=c_T k.$

In the vortex-line-lattice regime one has

$\Omega$0

while in the lowest-Landau-level regime reviewed for rapidly rotating weakly interacting condensates,

$\Omega$1

(Sonin, 2013).

Compressibility softens the mode. In a homogeneous compressible superfluid the standard long-wavelength dispersion is

$\Omega$2

with two asymptotic limits: $\Omega$3 This interpolates between the “soft” quadratic regime and the “stiff” linear regime, and it is the central hydrodynamic expression used in both the review literature and microscopic trapped-BEC analyses (Sonin, 2013, Simula et al., 2010).

A dual-gauge-field effective theory makes the same structure explicit. In that formulation the low-energy spectrum contains a gapped Kohn mode with $\Omega$4 and a Tkachenko mode whose full long-wavelength dispersion is again

$\Omega$5

Because rotation breaks parity and time-reversal symmetry, the Tkachenko eigenmode is elliptically polarized; for a plane wave with $\Omega$6 along $\Omega$7, the small longitudinal component obeys

$\Omega$8

showing a $\Omega$9 phase shift between components (Moroz et al., 2018).

3. Microscopic formulations in Bose–Einstein condensates

Microscopically, Tkachenko modes are obtained by solving Gross–Pitaevskii ground states and then linearizing with Bogoliubov–de Gennes theory. In rotating trapped condensates this is done in the rotating frame,

$n_v=2\Omega/\kappa$0

followed by the BdG eigenproblem for quasiparticle amplitudes $n_v=2\Omega/\kappa$1. Full three-dimensional calculations in oblate traps show that a condensate with $n_v=2\Omega/\kappa$2 vortices has $n_v=2\Omega/\kappa$3 low-lying Kelvin–Tkachenko branches, with axial-node structure and mode patterns determined directly from the quasiparticle fields (Simula, 2013).

In the finite lowest-Landau-level disk problem, the condensate is expanded in lowest-Landau-level orbitals,

$n_v=2\Omega/\kappa$4

with the droplet radius set by $n_v=2\Omega/\kappa$5. The Gross–Pitaevskii energy becomes a quartic form in the coefficients $n_v=2\Omega/\kappa$6, the ground state is a vortex crystal with a triangular bulk lattice and a circularly ordered rim, and BdG linearization produces a paired spectrum $n_v=2\Omega/\kappa$7 together with two exact zero modes from global $n_v=2\Omega/\kappa$8 phase and global rotation (Jeevanesan et al., 2022).

In the thermodynamic lowest-Landau-level limit, the Tkachenko mode has an exact analytic BdG solution. The full dispersion is

$n_v=2\Omega/\kappa$9

where $\kappa=h/m$0 and $\kappa=h/m$1 are lattice sums that can be written with Jacobi theta functions. At low momentum this reduces to the isotropic quadratic law

$\kappa=h/m$2

with $\kappa=h/m$3, $\kappa=h/m$4, and $\kappa=h/m$5. This microscopic result coincides with the hydrodynamic expectation for a compressible rapidly rotating vortex lattice (Matveenko et al., 2010).

4. Finite geometry, branch structure, and special modes

In finite trapped systems the continuum mode becomes a discrete family. For few-vortex arrays, both microscopic BdG work and classical point-vortex reasoning yield $\kappa=h/m$6 distinct Kelvin–Tkachenko branches for a configuration with $\kappa=h/m$7 vortices. In three-dimensional oblate traps these branches include the universal common mode $\kappa=h/m$8, the lowest-energy Tkachenko shear mode $\kappa=h/m$9, quadratic modes $\omega \simeq c_s k$0 for multi-orbital arrays, and rational modes $\omega \simeq c_s k$1 in which the central vortex or cluster remains stationary while outer rings shear (Simula, 2013).

The detailed taxonomy depends strongly on geometry. In single-ring arrays the branches can be labeled by the phase increment

$\omega \simeq c_s k$2

between neighboring vortices; for two vortices this yields acoustic and optical branches, while for three and four vortices it produces discrete triad and “molecular” torsional patterns. The lowest-frequency modes are usually torsional, out-of-phase shear modes with strongly suppressed radial motion (Simula et al., 2010).

A particularly important finite-disk realization is the Ruderman mode, which is an axisymmetric torsional Tkachenko solution. In the lowest-Landau-level disk, imposing the stress-free boundary condition $\omega \simeq c_s k$3 quantizes the wave number through

$\omega \simeq c_s k$4

and the corresponding frequency is

$\omega \simeq c_s k$5

In that same finite disk, bulk Tkachenko modes form a low-frequency standing-wave manifold, while chiral surface waves localized on the outer vortex ring occur at higher frequencies, with $\omega \simeq c_s k$6 at large $\omega \simeq c_s k$7 (Jeevanesan et al., 2022).

Finite systems also produce special soft modes. In a seven-vortex three-dimensional trapped condensate, the lowest Tkachenko branch was found to reach zero excitation energy near $\omega \simeq c_s k$8. The identified mechanism is an angular-momentum matching condition,

$\omega \simeq c_s k$9

at which the Magnus-force-induced coupling that ordinarily enforces elliptic motion vanishes and the vortex motion becomes purely azimuthal. The analysis explicitly characterizes this as a discrete trapped-system soft mode rather than a symmetry-protected Goldstone mode (Simula, 2012).

5. Symmetry, damping, and nonlinear effective theory

The symmetry structure of the vortex lattice is subtler than the naive counting of broken generators suggests. An effective-field-theory analysis of the rotating superfluid with a vortex lattice argues that only one gapless Nambu–Goldstone mode remains: the Tkachenko mode. The general counting rule is

$\omega$0

and the key operator identity in the rotating frame,

$\omega$1

implies that the would-be phonons of broken translations are redundant with the $\omega$2 phase mode. In this formulation the Tkachenko mode is identified with the Bogoliubov phase mode, while the displacement field contains a gapped cyclotron-like excitation with gap $\omega$3 (Watanabe et al., 2013).

Microscopic damping theory shows that the mode is underdamped at zero temperature in the mean-field lowest-Landau-level regime but can become strongly damped at finite temperature. For the exact lowest-Landau-level Tkachenko mode, Beliaev decay gives

$\omega$4

and

$\omega$5

At finite temperature and low energy, however, the damping becomes momentum-independent,

$\omega$6

so that modes with

$\omega$7

are overdamped (Matveenko et al., 2010).

More recent lowest-Landau-level field theories reformulate the Tkachenko mode as a compact $\omega$8 Lifshitz scalar or as a noncommutative Goldstone field. In the noncommutative formulation, magnetic translations act as a noncommutative dipole symmetry, the density operators satisfy the Girvin–MacDonald–Platzman algebra, and the zero-temperature decay width scales as

$\omega$9

in the low-energy limit (Du et al., 2022). In the Lifshitz-scalar and tensor-gauge description of vortex-crystal melting, the leading quadratic action is

$u(x,y,t)$0

with

$u(x,y,t)$1

and the scaling dimension of vacancy/interstitial operators becomes marginal at

$u(x,y,t)$2

providing a concrete quantum-melting criterion in the fast-rotation regime (Nguyen et al., 2023).

6. Experimental realizations, extensions, and broader settings

In ultracold-atom systems, Tkachenko modes can be probed by direct imaging of vortex positions, by trap modulation, and by Bragg or related spectroscopy. In finite lowest-Landau-level droplets, imaging the standing-wave pattern $u(x,y,t)$3, varying $u(x,y,t)$4 and $u(x,y,t)$5, and separating the low-frequency bulk branch from the higher-frequency chiral edge branch are explicit diagnostic strategies. The robust hierarchy $u(x,y,t)$6 helps isolate the bulk torsional excitations (Jeevanesan et al., 2022).

The mode is also sensitive to synthetic gauge engineering. In a quasi-two-dimensional condensate with density-dependent nonlinear rotation $u(x,y,t)$7, BdG, hydrodynamic, and time-domain vortex-tracking analyses all identify Tkachenko, circular, quadratic, and rational vortex-displacement modes. The measured Tkachenko frequencies extracted from the transverse and longitudinal vortex motion are

$u(x,y,t)$8

and the surface-mode frequency shifts from $u(x,y,t)$9 at $v(x,y,t)$0 to $v(x,y,t)$1 at $v(x,y,t)$2. In the same setting the dipole frequency

$v(x,y,t)$3

depends explicitly on $v(x,y,t)$4, signaling a violation of Kohn’s theorem by the density–angular-momentum coupling (Boral et al., 2 May 2025).

Quadratic Tkachenko excitations also modify impurity physics. In the “Tkachenko polaron” problem, an impurity moving in a vortex lattice couples to phonons with $v(x,y,t)$5 and vertex $v(x,y,t)$6. Because the Tkachenko continuum touches the impurity parabola at every $v(x,y,t)$7, the impurity spectral function develops a Lorentzian peak with finite width at arbitrarily small momentum even at $v(x,y,t)$8, with

$v(x,y,t)$9

This contrasts directly with both optical and acoustic phonon problems (Caracanhas et al., 2013).

Beyond cold atoms, Tkachenko waves have long been discussed in neutron-star interiors. In a two-fluid neutron-star model, compressibility, chemical coupling, mutual friction, and pinning can drastically modify or even eliminate the oscillatory branch. For slow pulsars rotating at $c_T=\sqrt{\frac{C_{66}}{\rho}},$0 Hz, many parameter choices still produce Tkachenko periods in the $c_T=\sqrt{\frac{C_{66}}{\rho}},$1 day range relevant to timing noise; for faster pulsars above $c_T=\sqrt{\frac{C_{66}}{\rho}},$2 Hz, a large portion of parameter space contains no propagating Tkachenko modes and only modified sound waves at much higher frequencies (Haskell, 2010).

The cumulative picture is therefore regime-dependent but coherent. Across continuum hydrodynamics, microscopic BdG theory, lowest-Landau-level effective field theory, and modern dual formulations, the Tkachenko mode remains the defining low-energy shear excitation of an ordered vortex crystal: soft, predominantly transverse, highly sensitive to compressibility and geometry, and central to the spectroscopy of rotating bosonic superfluids.