Quantized Vortex Dynamics of the Coupled Nonlinear Schr\"odinger Equation

Yongxing Zhu

REVIEW 2 major objections 5 minor 13 references

Quantized Vortex Dynamics of the Coupled Nonlinear Schr\"odinger Equation

T0 review · 2 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The coupled nonlinear Schrödinger equation's vortices decouple in the small-core limit: each component follows the single-component vortex motion law.

desk verdict New and likely correct result on CNLS vortex dynamics, but Lemma 3.1 has a real gap and Lemma 3.2 has a factor-2 error, so it needs revision before acceptance. read the letter →

arxiv 2411.14753 v1 pith:S57IFSHX submitted 2024-11-22 math.AP

Yongxing Zhu This is my paper · ORCID

classification math.AP MSC 35Q5535Q41

keywords couplednonlinearSchrödingerequationquantizedvortexcanonicalharmonicmapreduceddynamicallawrenormalizedenergyfractionalpathBose-Einsteincondensation

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper gives a rigorous derivation of the limiting vortex motion for the coupled nonlinear Schrödinger equation without Josephson junction, in two space dimensions, as the vortex core size $\varepsilon$ tends to zero. The main claim is that the vortices in the two components decouple: the motion of a vortex in $u$ is governed by the same renormalized-energy ODE as in the scalar nonlinear Schrödinger equation, and the $v$-component does not appear in it, and vice versa. A sympathetic reader should care because multi-component Bose-Einstein condensates are described by this system, and fractional vortices had previously been studied numerically or in two-vortex special cases, but no limiting ODE for the many-vortex, general case had been established. The paper also constructs the single-vortex core constant $\gamma_g$ through a radial ODE and reports a numerical estimate for one coupling value.

What carries the argument

The central objects are the renormalized energy $W_d(a)$ and the canonical harmonic map $H_d(x;a)$, the solution of $\nabla\cdot j(H)=0$ and $\nabla\cdot(Jj(H))=2\pi\sum_j d_j\delta_{a_j}$ with zero normal current on the boundary. $W_d$ encodes the logarithmic interactions of the vortex positions, and its gradient supplies the right-hand side of the vortex ODE. The argument works by comparing the true solution to the two-component comparison profile $(u^*,v^*)=((1+g)^{-1/2}H_{d1}(x;a),(1+g)^{-1/2}H_{d2}(x;b))$, using a refined lower bound that controls the excess energy in a ball around each vortex in terms of the distance of the vorticity from a single delta. A second ingredient is the radial ODE problem for the vortex profile, whose solution defines $\gamma_g$ and whose decay estimates make the energy expansion precise. These pieces are assembled through the derivative identity that converts the equations of motion into an estimate on the distance between the true and comparison vortex paths.

What would settle it

Numerically solve the coupled equations with decreasing core size $\varepsilon$ (for example $10^{-3},10^{-4},10^{-5}$) for a configuration with one $u$-vortex and one $v$-vortex separated by a fixed positive distance, starting from data that satisfies the theorem's assumptions. If the tracked $v$-vortex path does not converge to the solution of $\dot b = -\frac{d2}{\pi}J\nabla_b W_{d2}(b)$, with no dependence on the $u$-vortex position, then the decoupling claim would be falsified.

Watch

Extended reading notes

Core claim

Under Theorem 1.1, the paper proves the following. Assume the initial Jacobians of $u^\varepsilon$ and $v^\varepsilon$ concentrate as $\frac{\pi}{1+g}\sum_j d1_j\delta_{a^0_j}$ and $\frac{\pi}{1+g}\sum_k d2_k\delta_{b^0_k}$, and the initial energy matches the sum of two single-component renormalized energies with core constant $\gamma_g$. Then there is a time interval $[0,T)$ on which the vorticity measures converge in $W^{-1,1}$ to delta measures at Lipschitz paths $a_j(t)$, $b_k(t)$, and those paths satisfy the decoupled system $\dot a_j = -\frac{d1_j}{\pi}J\nabla_{a_j}W_{d1}(a)$, $\dot b_k = -\frac{d2_k}{\pi}J\nabla_{b_k}W_{d2}(b)$. In words: in the small-core limit each component's vortices move exactly as vortices of a single nonlinear Schrödinger equation, and the two components interact only through the common coupling constant $g$ that rescales the vortex strength. The proof also establishes existence of the core constant $\gamma_g$ via radial minimizers and gives a numerical value for $g=1/2$.

Load-bearing premise

The load-bearing premise is that the initial data is well prepared: the vortex centers of the two components are separated by a distance much larger than the core size, and the leading-order energy is exactly the sum of two single-component energies with no comparable cross-component term; if a $u$-vortex and a $v$-vortex start within $O(\varepsilon)$ of each other, or if a cross energy of the same order is present, the decoupled ODE is not established and may fail.

Editorial extensions

If this is right

As $\varepsilon\to0$, each component's vortex paths solve the scalar nonlinear Schrödinger vortex ODE; the other component never appears in the limiting equations.
The fractional vortices of the coupled system can thus be tracked by integer-winding zeros of each component, with the coupling constant $g$ entering only through the prefactor $1/(1+g)$ multiplying vortex strengths and the core constant $\gamma_g$.
The total vortex dynamics is a superposition of two independent vortex systems, so for well-separated vortices there is no cross-component scattering at leading order.
The existence and asymptotic form of the core constant $\gamma_g$ follow from a radial ODE, making the core energy computable from a one-dimensional problem.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

A likely extension is to remove the separation assumption: if a $u$-vortex and a $v$-vortex start within $O(\varepsilon)$, the decoupled ODE is not proven, and a collision-scale cross term may enter at leading order.
A related question the paper leaves open is whether the next-order correction gives a $O(1/\log(1/\varepsilon))$ cross-component force; numerical experiments at moderate $\varepsilon$ could check when the decoupled law becomes reliable.
The same two-scaled Jacobian comparison might apply to three-component condensates with more general coupling matrices, provided a similar single-vortex core constant exists.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Bottom line: this is the first rigorous reduced vortex-dynamics law for CNLS without Josephson coupling, and I think the main theorem is correct. But the manuscript as written has a genuine gap in Lemma 3.1 and a wrong factor 2 in Lemma 3.2. It deserves a serious referee, and the referee should ask for a real revision, not just typo fixes.

What is actually new: Theorem 1.1 shows that, under well-prepared, pairwise-distinct initial data, each component's vortices converge to Lipschitz paths satisfying the single-component NLS ODEs with no cross term. That is a meaningful extension of the Colliander-Jerrard and Jerrard-Spirn machinery to the coupled system. The existence of the renormalized constant gamma_g via the radial ODE in Section 5 is also a genuine contribution. The energy lower bound Lemma 3.3 is the expected core estimate and is handled with real detail. The paper cites the relevant prior work, including the special two-vortex results, and the use of [13] for earlier identities is legitimate.

The soft spots, in decreasing order. The stress-test note is right about Lemma 3.1. The estimate (3.21) controls |a_j(t2)-a_j(t1)| by sup_t (||grad u_eps||^2 + ||grad v_eps||^2)|t2-t1|, and that sup is of order log(1/eps) by energy conservation. So the constant is not independent of eps and Arzela-Ascoli does not follow as written. This is load-bearing because it produces the limiting vortex paths. The likely repair is to use that phi is affine near the vortex, so Hess phi vanishes there and the log divergence is cut off by (3.19), but the paper does not make that argument. Second, (3.30) is off by a factor of 2: since J(u) = (1/2) div(J j(u)), the weak limit must satisfy div(J j_a) = 2pi/(1+g) sum d delta, not pi/(1+g) sum. Without the 2, Lemma 3.2's identification j_a = j(u*) does not close. This looks like a typo, but it is at the hinge of the convergence proof. Third, Section 5 has several typos, most visibly (5.23) and the integration limits in (5.24); these are minor once the ODE argument is checked.

The abstract also overstates the range of the decoupling by not stating the well-preparedness and separation assumptions, but the theorem itself is honest.

Who this is for: researchers in vortex dynamics for multi-component BEC and Ginzburg-Landau. I would not desk-reject it. Send it to a referee who knows the Colliander-Jerrard machinery, and require the Lemma 3.1 gap and the factor 2 to be fixed. If the announced repair works, which I expect, the theorem stands.

Referee Report

2 major / 5 minor

Summary. The paper studies the coupled nonlinear Schrödinger system (CNLS) without Josephson coupling in a bounded 2D domain with Neumann boundary conditions. Its main theorem states that for well-prepared initial data with finitely many degree ±1 vortices in each component, the Jacobians of u_ε and v_ε converge as ε→0 to sums of Dirac measures concentrated at Lipschitz vortex paths a_j(t), b_k(t), and these paths satisfy the decoupled single-component ODEs (1.18). The proof follows the standard vortex-method structure: energy conservation and lower bounds, extraction of vortex paths, convergence of currents, an energy-deficit estimate, and a Gronwall comparison with the limiting ODE. Section 5 establishes existence and asymptotics of the core-energy constant γ_g and reports a numerical estimate γ_{1/2}≈0.5377.

Significance. If the main theorem is correct, this is a substantial rigorous extension of single-component Ginzburg-Landau vortex dynamics to a coupled two-component system, and it gives a concrete, falsifiable prediction: in the small-core limit the vortex motion of one component does not affect the other component, and each component follows the classical NLS vortex law. The treatment of γ_g in Section 5, including the explicit two-term expansions in Lemma 5.6 and the numerical value, is a useful contribution. The paper is not machine-checked and relies on several nontrivial imported results from [3], [7], [8], and [13]; nevertheless, the overall strategy is coherent and the specific gaps identified below appear locally repairable.

major comments (2)

[§3.2, Lemma 3.1, Eq. (3.21)] The compactness step producing the limiting Lipschitz paths is not justified as written. In (3.21) the chain of inequalities ends with a constant C asserted to be independent of ε, but the previous line contains sup_t(‖∇u_ε(t)‖²_{L²(Ω)} + ‖∇v_ε(t)‖²_{L²(Ω)}), which by conservation of energy and assumption (1.16) is of order (M+N)π/(1+g) log(1/ε). Thus the displayed estimate does not give a uniform Lipschitz bound, and Arzelà-Ascoli cannot be applied from the written inequality. This is load-bearing because Lemma 3.1 supplies the paths ã_j, b̃_k on which Theorem 1.1 and all later arguments rest. The gap seems repairable: choose φ to be linear on a slightly enlarged fixed ball B_{3r0/4}(a0_j), so that Hess φ vanishes on the ε-core of the j-th vortex, and bound the remaining integral using the uniform energy bound (3.19) on the annulus away from all vortex cores. With that modification the relevant gradient L² norms are uniformly bounded in ε. The proof should also explicitly show that the maximal time T_ε in (3.9) is bounded below by a positive constant independent of ε; otherwise the subsequent limiting construction may be vacuous.
[§3.3, Lemma 3.2, Eq. (3.30)] Equation (3.30) is missing a factor 2. Since J(u) is defined in (1.14) by J(u)=1/2 ∇·(Jj(u)), the convergence J(u_ε)→π/(1+g)∑_j d1_j δ_{ã_j(t)} in (3.7) implies div(J j_a)=2π/(1+g)∑_j d1_j δ_{ã_j(t)}, not π/(1+g) as written. With the displayed (3.30), the subsequent identity div(J(j_a−j(u*)))=0 in (3.31) is false, because j(u*)=j(H_{d1})/(1+g) satisfies div(Jj(u*))=2π/(1+g)∑_j d1_j δ_{ã_j(t)} by (2.7). This is a local typo, but as written it invalidates the identification j_a=j(u*) in Lemma 3.2 and must be corrected before the proof proceeds.

minor comments (5)

[§3.4, Eq. (3.45)] In the sums involving b_j and b̃_j, the symbols d1_j and W_{d1} should be d2_j and W_{d2}; the displayed inequality otherwise confuses the two components in the Gronwall comparison.
[§3.3, Lemma 3.2 proof] The sentence 'Combing (3.17), (1.17) and Theorem 1.4.3 in [3]' should refer to (3.7), not (1.17), since (1.17) is the conclusion of Theorem 1.1 and is not yet available at that point.
[§5, Lemma 5.3, Eq. (5.23)] The displayed identity '(rf2)′ = rf2′ + f2 = rh(r)' is inconsistent with the subsequent integration step; the correct identity is (r f2′)′ = r h(r).
[§5, Lemma 5.2 proof, around (5.17)] The line 'h^k_2(r)=√r H3' is a typo: H3 was defined as the weak limit of h^k_1/√r, so the displayed relation should refer to h^k_1.
[§3.2, Eq. (3.21)] The notation ‖φ‖_{C²(Ω×[t1,t2])} is misleading because φ is time-independent; this should be ‖φ‖_{C²(Ω)}.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the vortex ODE is derived from the PDE and independently defined constants; self-citations are not load-bearing.

full rationale

The central claim of Theorem 1.1 is not an input to the argument. The initial-data assumptions (1.15)-(1.16) prescribe only the vortex locations and the additive energy expansion; they do not prescribe the velocities or the final ODE (1.18). The constant gamma_g is defined as a variational limit in (1.11) and its existence is proved independently in Section 5, not fitted to the target dynamics. The renormalized energy W_d(a) is defined through the canonical harmonic map in (2.7), with the needed identities (2.8) and (2.11) quoted from prior external work [8,3]. The only self-citation is to [13], used for algebraic identities in the energy decomposition and in the derivation of the path acceleration; these are published identities from a prior paper about single-component NLS and do not assume the CNLS decoupling law. Thus even the self-citation does not force the conclusion. The skeptical concerns about Lemma 3.1, such as the uniform Lipschitz constant in (3.21) and the apparent missing factor 2 in (3.30), concern correctness or rigor of the proof, not circularity: they do not make the target ODE equivalent to an assumption. No equation in the manuscript is shown to reduce to the theorem's conclusion by construction, so the paper receives a score of 0.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

No free parameters are fitted: gamma_g is defined as a variational limit and its existence is proved; the numerical value 0.5377 in Remark 5.8 is an illustration and is not used in the theorem. The central claim depends on standard imported vortex theorems and on the well-prepared-data hypothesis, listed as axioms. No new particles, forces, or entities are introduced.

assumptions (4)

standard math Known NLS vortex dynamics theorems from Colliander-Jerrard [3], including Theorems 1.4.3, 1.4.4 and Lemmas 3.1.3, 3.1.4
Used in Section 3.2 to locate vortex centers and to obtain energy lower bounds for sqrt(1+g)u_epsilon and sqrt(1+g)v_epsilon via the comparison (2.6).
standard math Canonical harmonic map estimates from Jerrard-Spirn [8], especially equation (2.8)
Defines H_d(x;a) and the renormalized energy expansion for the comparison maps u*,v* in Section 2 and throughout the proof.
standard math Uniqueness and asymptotic results for the radial ODE from Ignat-Nguyen-Slastikov-Zarnescu [6], Theorem 1.1 and Proposition 3.5
Used in Section 5 to establish existence, convergence, and decay of radial minimizers defining gamma_g.
domain assumption Well-prepared initial data assumption (1.15)-(1.16)
The theorem is conditional on exact delta-concentration of both Jacobians at pairwise distinct points and additivity of the renormalized energy at t=0. The proof uses this to control vortex paths and energy on positive time intervals.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Quantized Vortex Dynamics of the Coupled Nonlinear Schr\"odinger Equation." pith.science (2026). https://pith.science/paper/S57IFSHX

@misc{pith2026241114753,
  author       = {Pith},
  title        = {Pith review of: Quantized Vortex Dynamics of the Coupled Nonlinear Schr\"odinger Equation},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/S57IFSHX}},
  note         = {Machine review of arXiv:2411.14753}
}

abstract

We derive rigorously the reduced dynamical law for quantized vortex dynamics of the coupled nonlinear Schr\"odinger equation without Josephson junction (CNLS) when the core size of vortex $\varepsilon\to 0$. It is proved that when $\varepsilon\to 0$, the vortex motion of one component won't affect the vortex motion on the other component. Moreover, the motion of vortices of each component follows the vortex motion law for the nonlinear Schr\"odinger.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

[13]

Quantized vo rtex dynamics of the nonlinear schr¨ odinger equation on tor us with non-vanishing momentum

Yongxing Zhu, W eizhu Bao, and Huaiyu Jian. Quantized vo rtex dynamics of the nonlinear schr¨ odinger equation on tor us with non-vanishing momentum. Physica D: Nonlinear Phenomena , 453:133812, 2023. 25

work page 2023
[3]

J. E. Colliander and R. L. Jerrard. Ginzburg-landau vort ices: weak stability and schr¨ odinger equation dynamics. Journal d’Analyse Math´ ematique, 77(1):129–205, 1999

work page 1999
[7]

Jerrard and Daniel Spirn

Robert L. Jerrard and Daniel Spirn. Reﬁned jacobian esti mates for ginzburg-landau functionals. Indiana University Mathematics Journal , 56:135–186, 2007

work page 2007
[8]

Reﬁned jacobian es timates and gross–pitaevsky vortex dynamics

Robert Leon Jerrard and Daniel Spirn. Reﬁned jacobian es timates and gross–pitaevsky vortex dynamics. Archive for Rational Mechanics and Analysis , 190(3):425–475, 2008

work page 2008
[1]

Mathematical models and num erical methods for spinor bose-einstein condensates

W eizhu Bao and Yongyong Cai. Mathematical models and num erical methods for spinor bose-einstein condensates. Commun. Comput. Phys. , 24(4):899–965, 2018

work page 2018
[2]

Numerical study of quantize d vortex interactions in the nonlinear Schr¨ odinger equati on on bounded domains

W eizhu Bao and Qinglin Tang. Numerical study of quantize d vortex interactions in the nonlinear Schr¨ odinger equati on on bounded domains. Multiscale Modeling & Simulation , 12(2):411–439, 2014

work page 2014
[4]

Collision dy namics and reactions of fractional vortex molecules in cohe rently coupled bose-einstein condensates

Minoru Eto, Kazuki Ikeno, and Muneto Nitta. Collision dy namics and reactions of fractional vortex molecules in cohe rently coupled bose-einstein condensates. Physical Review Research, 2(3):033373, 2020. PRRESEARCH

work page 2020
[5]

Many-Particle Entanglement, Einstein-Podolsky-Rosen St eering and Bell Correlations in Bose-Einstein Condensates

Matteo Fadel. Many-Particle Entanglement, Einstein-Podolsky-Rosen St eering and Bell Correlations in Bose-Einstein Condensates. Quantum Science and Technology. Springer Cham, 1 edition, 2021

work page 2021

Show all 13 references

[6]

Uniqueness results for an ode related to a generaliz ed ginzburg–landau model for liquid crystals

Radu Ignat, Luc Nguyen, Valeriy Slastikov, and Arghir Za rnescu. Uniqueness results for an ode related to a generaliz ed ginzburg–landau model for liquid crystals. SIAM Journal on Mathematical Analysis , 46(5):3390–3425, 2014

2014
[9]

F. H. Lin and J. X. Xin. On the incompressible ﬂuid limit an d the vortex motion law of the nonlinear schr¨ odinger equation. Communications in Mathematical Physics , 200(2):249–274, 1999

1999
[10]

Complex ginzburg-landau equations and dy namics of vortices, ﬁlaments, and codimension-2 submanifo lds

Fang Hua Lin. Complex ginzburg-landau equations and dy namics of vortices, ﬁlaments, and codimension-2 submanifo lds. Communications on Pure and Applied Mathematics , 51(4):385–441, 1998

1998
[11]

Multiple time scale dynam ics in coupled ginzburg-landau equations

Fanghua Lin and Tai-Chia Lin. Multiple time scale dynam ics in coupled ginzburg-landau equations. Communications in Mathematical Sciences, 1(4):671–695, 2003

2003
[12]

Skyrmions i n gross-pitaevskii functionals

Fanghua Lin, Taichia Lin, and Juncheng W ei. Skyrmions i n gross-pitaevskii functionals. Acta Mathematica Scientia , 29(3):751–776, 2009

2009

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.

[1] [13]

Quantized vo rtex dynamics of the nonlinear schr¨ odinger equation on tor us with non-vanishing momentum

Yongxing Zhu, W eizhu Bao, and Huaiyu Jian. Quantized vo rtex dynamics of the nonlinear schr¨ odinger equation on tor us with non-vanishing momentum. Physica D: Nonlinear Phenomena , 453:133812, 2023. 25

work page 2023

[2] [3]

J. E. Colliander and R. L. Jerrard. Ginzburg-landau vort ices: weak stability and schr¨ odinger equation dynamics. Journal d’Analyse Math´ ematique, 77(1):129–205, 1999

work page 1999

[3] [7]

Jerrard and Daniel Spirn

Robert L. Jerrard and Daniel Spirn. Reﬁned jacobian esti mates for ginzburg-landau functionals. Indiana University Mathematics Journal , 56:135–186, 2007

work page 2007

[4] [8]

Reﬁned jacobian es timates and gross–pitaevsky vortex dynamics

Robert Leon Jerrard and Daniel Spirn. Reﬁned jacobian es timates and gross–pitaevsky vortex dynamics. Archive for Rational Mechanics and Analysis , 190(3):425–475, 2008

work page 2008

[5] [1]

Mathematical models and num erical methods for spinor bose-einstein condensates

W eizhu Bao and Yongyong Cai. Mathematical models and num erical methods for spinor bose-einstein condensates. Commun. Comput. Phys. , 24(4):899–965, 2018

work page 2018

[6] [2]

Numerical study of quantize d vortex interactions in the nonlinear Schr¨ odinger equati on on bounded domains

W eizhu Bao and Qinglin Tang. Numerical study of quantize d vortex interactions in the nonlinear Schr¨ odinger equati on on bounded domains. Multiscale Modeling & Simulation , 12(2):411–439, 2014

work page 2014

[7] [4]

Collision dy namics and reactions of fractional vortex molecules in cohe rently coupled bose-einstein condensates

Minoru Eto, Kazuki Ikeno, and Muneto Nitta. Collision dy namics and reactions of fractional vortex molecules in cohe rently coupled bose-einstein condensates. Physical Review Research, 2(3):033373, 2020. PRRESEARCH

work page 2020

[8] [5]

Many-Particle Entanglement, Einstein-Podolsky-Rosen St eering and Bell Correlations in Bose-Einstein Condensates

Matteo Fadel. Many-Particle Entanglement, Einstein-Podolsky-Rosen St eering and Bell Correlations in Bose-Einstein Condensates. Quantum Science and Technology. Springer Cham, 1 edition, 2021

work page 2021