Co-rotating Vortices on Surfaces of Variable Negative Curvature: Hamiltonian Structure and Curvature-Induced Drift
Pith reviewed 2026-05-21 08:16 UTC · model grok-4.3
The pith
On a catenoid, identical co-rotating vortices rotate rigidly at fixed latitude with speed set by the Gaussian curvature gradient.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For two identical vortices the authors obtain an exact analytic solution in which the pair rotates rigidly at fixed latitude with angular velocity Ω = (Γ/16π) K'(V)/√(-K(V)), where K(V) is the Gaussian curvature. This establishes that vortex motion is governed by the curvature gradient. Conservation of the Hamiltonian and rotational momentum reduces the generic co-rotating case to a quadrature, producing bounded relative oscillations plus a secular azimuthal drift. The rigid state is linearly unstable with growth rate λ = √3 |Ω|, and the same curvature-induced drift appears in localized many-vortex clusters.
What carries the argument
The Hamiltonian structure of point-vortex motion on the catenoid together with conservation of rotational momentum, which reduces the relative dynamics to a quadrature whose effective velocity contains the curvature gradient explicitly.
If this is right
- The rigid-rotation state is linearly unstable with growth rate √3 |Ω|.
- Generic co-rotating pairs exhibit bounded relative oscillations together with a secular azimuthal drift.
- Localized many-vortex clusters display the same curvature-induced azimuthal drift.
- The reduction to quadrature supplies an exact benchmark for testing full numerical simulations of the vortex equations.
Where Pith is reading between the lines
- The explicit dependence on the curvature gradient suggests a similar drift mechanism may operate on other surfaces whose Gaussian curvature varies spatially.
- The secular drift identified in clusters provides a starting point for a collective description of many-vortex systems on curved manifolds.
- The exact rigid-rotation solution offers a clean test case for validating numerical integrators in non-flat geometries.
Load-bearing premise
Vortex motion on the surface obeys a Hamiltonian structure that conserves rotational momentum and thereby permits reduction of the relative dynamics to a quadrature.
What would settle it
A direct numerical integration or laboratory measurement of the angular velocity of two identical co-rotating vortices held at fixed latitude on a catenoid that fails to match the formula Ω = (Γ/16π) K'(V)/√(-K(V)).
Figures
read the original abstract
Vortices in fluids and superfluids are fundamental to phenomena ranging from Bose-Einstein condensates and superfluid films to neutron stars and hydrodynamic micro-rotors, where background geometry often plays an important role. Curvature can induce vortex motion distinct from planar domains. We study Hamiltonian vortex motion on a catenoid, a minimal surface of variable negative curvature, and derive explicit equations of motion and conserved quantities for co-rotating vortex pairs. For two identical vortices we find an exact analytic solution in which the pair rotates rigidly at fixed latitude, with angular velocity $\Omega=(\Gamma/16\pi)\,K'(V)/\sqrt{-K(V)}$, where $K(V)$ is the Gaussian curvature. Thus the motion is governed by the curvature gradient rather than the curvature itself. This state is linearly unstable, with growth rate $\lambda=\sqrt{3}|\Omega|$, in agreement with numerical simulations. For generic co-rotating pairs, conservation of the Hamiltonian and rotational momentum reduces the nonlinear dynamics to a single quadrature, yielding bounded relative oscillations together with a secular azimuthal drift. Simulations of the full equations confirm this and reveal the same curvature-induced azimuthal drift in a localized many-vortex cluster, motivating a broader theory of collective vortex drift on curved surfaces.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript derives the Hamiltonian structure and equations of motion for co-rotating vortex pairs on a catenoid (a minimal surface with variable negative Gaussian curvature K). For two identical vortices it presents an exact analytic solution in which the pair rotates rigidly at fixed latitude V with angular velocity Ω = (Γ/16π) K'(V)/√(-K(V)), governed by the curvature gradient. This equilibrium is shown to be linearly unstable with growth rate λ = √3 |Ω|, in agreement with direct numerical simulations. Conservation of the Hamiltonian and rotational momentum reduces the generic two-vortex problem to a single quadrature, yielding bounded relative oscillations accompanied by secular azimuthal drift. Simulations of the full system confirm the drift and extend the observation to a localized many-vortex cluster.
Significance. If the derivations are correct, the work supplies an explicit, parameter-free illustration that vortex motion on surfaces of non-constant negative curvature is driven by the gradient of K rather than by K itself. The reduction to quadrature, the closed-form instability rate, and the numerical confirmation for both pairs and clusters constitute concrete, falsifiable predictions that could seed a broader theory of collective curvature-induced drift in superfluids and fluids on curved geometries.
major comments (2)
- [§3] §3 (Hamiltonian construction) and the paragraph leading to Eq. (Ω): the exact rigid-rotation solution at constant V requires that the reduced two-vortex Hamiltonian be strictly of the form H = f(V, Δφ) with the precise prefactor arising from the Laplace-Beltrami Green's function on the catenoid metric. Any omitted curvature-dependent correction to the stream function would shift the equilibrium condition away from the stated Ω. Please exhibit the explicit Green's function (or stream function) used for the catenoid and confirm that no additional metric terms enter the effective potential.
- [Linear stability section] Linear-stability calculation (near Eq. (λ)): the growth rate λ = √3 |Ω| is obtained by linearizing the reduced equations about the fixed-latitude equilibrium. It is not immediately clear whether this factor of √3 is universal for any K(V) or whether it relies on the specific catenoid profile; an explicit Jacobian or characteristic equation should be shown to establish independence from higher-order curvature corrections.
minor comments (3)
- [Introduction / coordinate setup] The coordinate V (latitude) and the range over which -K(V) > 0 should be defined explicitly at first use, together with the relation between the catenoid metric and the Gaussian curvature expression.
- [Numerical results] Figure captions for the simulation panels should state the number of vortices, the value of Γ, and the integration time step or method used, to allow direct reproduction of the reported drift.
- [Many-vortex simulations] A brief remark on how the many-vortex cluster is initialized (e.g., random or lattice placement) would clarify the robustness of the observed collective drift.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of the work and for the constructive comments on the Hamiltonian derivation and linear stability analysis. We address each point below and have revised the manuscript to incorporate the requested explicit details.
read point-by-point responses
-
Referee: [§3] §3 (Hamiltonian construction) and the paragraph leading to Eq. (Ω): the exact rigid-rotation solution at constant V requires that the reduced two-vortex Hamiltonian be strictly of the form H = f(V, Δφ) with the precise prefactor arising from the Laplace-Beltrami Green's function on the catenoid metric. Any omitted curvature-dependent correction to the stream function would shift the equilibrium condition away from the stated Ω. Please exhibit the explicit Green's function (or stream function) used for the catenoid and confirm that no additional metric terms enter the effective potential.
Authors: We have added the explicit Green's function for the catenoid to the revised Section 3. On the metric ds² = du² + cosh²(u) dφ² the Green's function satisfying ΔG = δ(u-u₀,φ-φ₀) (with appropriate normalization) takes the integral form G(u,φ;u₀,φ₀) = (1/2π) ∫ dk e^{ik(φ-φ₀)} [cosh(k(u_<)) / sinh(k(u_>)) ] or the equivalent closed expression involving the hyperbolic distance on the surface. The resulting stream function for a vortex of strength Γ yields the two-vortex Hamiltonian H = -(Γ²/4π) log|sinh((u-u₀)/2)| plus the azimuthal term, with no further curvature corrections beyond those already encoded in the Laplace-Beltrami operator. Substituting into the reduced equations reproduces exactly the stated equilibrium Ω = (Γ/16π) K'(V)/√(-K(V)) at constant latitude, confirming that the rigid-rotation solution is free of omitted metric contributions. revision: yes
-
Referee: [Linear stability section] Linear-stability calculation (near Eq. (λ)): the growth rate λ = √3 |Ω| is obtained by linearizing the reduced equations about the fixed-latitude equilibrium. It is not immediately clear whether this factor of √3 is universal for any K(V) or whether it relies on the specific catenoid profile; an explicit Jacobian or characteristic equation should be shown to establish independence from higher-order curvature corrections.
Authors: The revised linear-stability section now displays the explicit 2×2 Jacobian matrix obtained by differentiating the reduced Hamiltonian equations with respect to (V,Δφ) and evaluating at the equilibrium (V,Δφ=π). The characteristic equation det(J-λI)=0 reduces to λ² - 3Ω² = 0, giving the unstable eigenvalue λ = √3 |Ω|. The factor √3 is a direct algebraic consequence of the trigonometric identities that appear when the catenoid curvature K(u) = -sech²(u) and its first derivative are inserted into the second derivatives of H; it is therefore specific to this surface rather than universal for arbitrary K(V). Because the linearization uses only the local value of K' at equilibrium and does not invoke higher-order curvature terms, the result is exact within the reduced two-vortex model and independent of any truncation of the metric expansion. revision: yes
Circularity Check
No significant circularity; derivation is self-contained via standard Hamiltonian reduction
full rationale
The paper derives the Hamiltonian for point vortices on the catenoid from the Laplace-Beltrami Green's function on the surface metric, identifies the ignorable azimuthal coordinate, and invokes conservation of rotational momentum to reduce the two-vortex dynamics to a single quadrature in the relative coordinate. The rigid-rotation solution at constant latitude follows by setting the derivative of the resulting effective potential to zero, producing the stated Ω = (Γ/16π) K'(V)/√(-K(V)) directly from the curvature gradient term. This is a standard integrable reduction with no parameter fitting, no self-referential definition of the target quantity, and no load-bearing self-citation chain. The analytic result is cross-checked against direct numerical integration of the unreduced equations, confirming the derivation remains independent of its own output.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The vortex system on the curved surface admits a Hamiltonian formulation with conserved quantities including rotational momentum.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
angular velocity Ω=(Γ/16π) K'(V)/√(-K(V)) ... motion is governed by the curvature gradient rather than the curvature itself
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
conservation of the Hamiltonian and rotational momentum reduces the nonlinear dynamics to a single quadrature
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Aref,Integrable, chaotic, and turbulent vortex motion in two-dimensional flows,Annu
H. Aref,Integrable, chaotic, and turbulent vortex motion in two-dimensional flows,Annu. Rev. Fluid Mech.15(1983) 345–389
work page 1983
-
[2]
P. G. Saffman,Vortex Dynamics(Cambridge Univ. Press, 1993)
work page 1993
-
[3]
C. C. Lin,On the motion of vortices in two dimensions: I. Existence of the Kirchhoff–Routh function,Proc. Natl. Acad. Sci. USA27(1941) 570–575. 25
work page 1941
-
[4]
C. C. Lin,On the motion of vortices in two dimensions: II. Some further investigations on the Kirchhoff–Routh function,Proc. Natl. Acad. Sci. USA27(1941) 575–577
work page 1941
-
[5]
A. A. Tchieu, E. Kanso, and P. K. Newton,The finite-dipole dynamical system, Proc. R. Soc. A468, 3006–3026 (2012)
work page 2012
- [6]
-
[7]
R. Samanta and N. Oppenheimer,Vortex flows and streamline topology in curved biological membranes,Phys. Fluids33(2021) 092111
work page 2021
- [8]
-
[9]
Aswathy K. R. , U. Maurya, S. T. Gavva and R. Samanta,Dynamics of vortex clusters on a torus,Phys. Fluids37(2025) 093324
work page 2025
-
[10]
K. Banthia and R. Samanta,A self propelled vortex dipole model on a surface of variable negative curvature, J. Phys. A: Math. Theor.59, 145701 (2026)
work page 2026
-
[11]
Aswathy K. R. and R. Samanta,Collective dynamics of vortex clusters in compact fluid do- mains: From pair interactions to a quadrupole description, arXiv:2604.07373 [physics.flu-dyn] (2026)
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[12]
V. A. Bogomolov,Dynamics of vorticity on a sphere,Fluid Dyn.12(1977) 863–870
work page 1977
-
[13]
Y. Kimura and H. Okamoto,Vortex motion on a sphere,J. Phys. Soc. Jpn.56(1987) 4203– 4206
work page 1987
-
[14]
Kimura,Vortex motion on surfaces with constant curvature,Proc
Y. Kimura,Vortex motion on surfaces with constant curvature,Proc. R. Soc. A455(1999) 245–259
work page 1999
-
[15]
R. Kidambi and P. K. Newton,Streamline topologies for integrable vortex motion on a sphere, Physica D140(2000) 95–125
work page 2000
-
[16]
R. Kidambi and P. K. Newton,Point vortex motion on a sphere with solid boundaries,Phys. Fluids12(2000) 581–588
work page 2000
-
[17]
Hally,Stability of streets of vortices on surfaces of revolution with a reflection symmetry, J
D. Hally,Stability of streets of vortices on surfaces of revolution with a reflection symmetry, J. Math. Phys.21(1980) 211–217
work page 1980
-
[18]
D.G. Crowdy and J. Marshall,Analytical formulae for the Kirchhoff–Routh path function in multiply connected domains,Proc. R. Soc. A461(2005) 2477–2501. 26
work page 2005
-
[19]
S. Boatto and J. Koiller,Vortices on closed surfaces, inGeometry, Mechanics, and Dynamics: The Legacy of Jerry Marsden, edited by D. E. Chang, D. D. Holm, G. Patrick, and T. Ratiu, Fields Institute Communications, Vol. 73 (Springer, New York, 2015), pp. 185–237
work page 2015
-
[20]
J. Koiller and S. Boatto,Vortex pairs on surfaces, AIP Conf. Proc.,1130, 77–88 (2009)
work page 2009
-
[21]
D. G. Dritschel and S. Boatto,The motion of point vortices on closed surfaces,Proc. R. Soc. A471(2015) 20140890
work page 2015
-
[22]
A. M. Turner, V. Vitelli and D. R. Nelson,Vortices on curved surfaces,Rev. Mod. Phys.82 (2010) 1301–1348
work page 2010
-
[23]
S. Reuther and A. Voigt,The interplay of curvature and vortices in flow on curved surfaces, Multiscale Model. Simul.13, 632–643 (2015)
work page 2015
-
[24]
Gustafsson,Vortex pairs and dipoles on closed surfaces,J
B. Gustafsson,Vortex pairs and dipoles on closed surfaces,J. Nonlinear Sci.32(2022) 62
work page 2022
-
[25]
T. D. Drivas, D. Glukhovskiy, and B. Khesin,Singular Vortex Pairs Follow Magnetic Geodesics, Int. Math. Res. Not.2024(14), 10880–10894 (2024)
work page 2024
-
[26]
M. A. Caracanhas, P. Massignan, and A. L. Fetter,Superfluid vortex dynamics on an ellipsoid and other surfaces of revolution, Phys. Rev. A105, 023307 (2022)
work page 2022
-
[27]
T. W. Neely, E. C. Samson, A. S. Bradley, M. J. Davis and B. P. Anderson,Observation of vortex dipoles in an oblate Bose–Einstein condensate,Phys. Rev. Lett.104(2010) 160401
work page 2010
-
[28]
D. V. Freilich, D. M. Bianchi, A. M. Kaufman, T. K. Langin and D. S. Hall,Real-time dynamics of single vortex lines and vortex dipoles in a Bose–Einstein condensate,Science329(2010) 1182–1185
work page 2010
-
[29]
G. Gauthier, M. T. Reeves, X. Yu, A. S. Bradley, M. Baker, T. A. Bell, H. Rubinsztein- Dunlop, M. J. Davis, and T. W. Neely,Giant vortex clusters in a two-dimensional quantum fluid, Science364, 1264–1267 (2019)
work page 2019
-
[30]
S. J. Rooney, P. B. Blakie, B. P. Anderson and A. S. Bradley,Suppression of Kelvon-induced decay of quantized vortices in oblate Bose-Einstein condensates,Phys. Rev. A84, 023637 (2011)
work page 2011
-
[31]
R. H. Goodman, P. G. Kevrekidis and R. Carretero-Gonz´ alez,Dynamics of Vortex Dipoles in Anisotropic Bose–Einstein Condensates,SIAM J. Appl. Dyn. Syst.14, no. 2, 699–729 (2015)
work page 2015
-
[32]
A. C. White, C. F. Barenghi and N. P. Proukakis,Creation and Characterization of Vortex Clusters in Atomic Bose-Einstein Condensates,Physical Review A86, 013635, (2012). 27
work page 2012
-
[33]
A. C. White, B. P. Anderson, and V. S. Bagnato,Vortices and turbulence in trapped atomic condensates, Proc. Natl. Acad. Sci. U.S.A.111, 4719–4726 (2014)
work page 2014
-
[34]
G. W. Stagg, N. G. Parker, and C. F. Barenghi,Ultraquantum turbulence in a quenched homogeneous Bose gas, Phys. Rev. A94, 053632 (2016)
work page 2016
-
[35]
R. Samanta,Co-rotating Vortex Clusters on Negatively Curved Geometries, manuscript in preparation (2026)
work page 2026
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.