{"id":"5cf4f7a3-6749-4091-b46c-616b90fc1347","arxiv_id":"2411.14753","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the coupled nonlinear Schrödinger equation without Josephson coupling, vortex motion in each component converges as the core size vanishes to the classical NLS vortex ODE, with no cross-component forcing.","lead":"This paper proves a mathematical law for how quantum vortices move in a two-component Bose-Einstein condensate. Under carefully prepared, well-separated starting conditions, each component's vortices follow the same simple point-vortex ODE as in a single-component condensate, and the two components do not push each other.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Lemma 3.1's vortex-path compactness is not established: (3.21) bounds |a_j(t2)-a_j(t1)| by an L^2 norm that grows like log(1/epsilon), so the asserted uniform Lipschitz constant is unjustified.","rationale":"The central claim is that, under the well-prepared initial data (1.15)-(1.16), each component's vortices obey the single-component NLS ODE (1.18) with no cross-coupling. The proof strategy is standard and the decoupling itself is plausible: for well-separated vortices the leading-order renormalized energy is the sum Wd1(a)+Wd2(b), with no cross term. However, the theorem's proof depends on Lemma 3.1 producing Lipschitz vortex paths with a uniform time-regularity constant. The only displayed estimate for this, (3.21), uses the full L^2 norm of the gradients, which by conservation of energy and (1.16) diverges logarithmically. This is not a cosmetic typo: without a replacement argument using Hess(phi)=0 near the cores and the local bound (3.19), the existence of the paths ~a_j, ~b_j is not proved. The reader's weakest-assumption concern about cross terms and well-preparedness is a valid scope limitation, and the abstract should indeed be qualified, but the path-compactness gap is more directly load-bearing for the submitted proof. The factor-2 error in (3.30) is also a key displayed-equation correction, as it affects the identification of the limiting current in Lemma 3.2. I do not claim the theorem is false; the necessary estimate is likely obtainable by a local-energy argument. For that reason the appropriate verdict remains conditional rather than rejection.","tokens_in":25678,"tokens_out":25856,"duration_ms":262586,"concrete_test":"Repair or refute the estimate in (3.21) by replacing the full L^2 gradient norm with the local energy on the support of Hess(phi). Check whether (3.19) gives a uniform bound on int_{supp Hess(phi)} e0(u_eps,v_eps); if yes, Lemma 3.1 can be repaired. If not, exhibit the logarithmic growth from (1.16) to show the displayed C|t2-t1| has no epsilon-independent constant. Separately, recompute (3.30) from (3.7) and the identity J = (1/2) div(Jj); the factor 2 should appear, and Lemma 3.2 should be re-verified with that correction.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing gap is in Lemma 3.1, the step that produces the Lipschitz vortex paths ~a_j, ~b_j and hence the convergence (3.7) on which Theorem 1.1 rests. In (3.21) the displayed chain bounds |a_j^eps(t2)-a_j^eps(t1)| by (1+g)/pi ||phi||_{C^2} sup_t (||grad u_eps(t)||^2_{L^2(Omega)} + ||grad v_eps(t)||^2_{L^2(Omega)}) |t2-t1| + o(1), and then asserts this is <= C|t2-t1| with C independent of epsilon. But conservation of energy and assumption (1.16) give ||grad u_eps||^2_{L^2} + ||grad v_eps||^2_{L^2} ~ 2(M+N) pi/(1+g) log(1/epsilon), so the displayed constant diverges as epsilon -> 0. Arzela-Ascoli therefore does not apply from the written estimate. The argument is likely repairable: phi is chosen linear in B_{r0/2}(a0_j), so Hess(phi)=0 inside the vortex core and the Hessian integrand is supported in an annulus where (3.19) supplies a uniform L^2 bound. But that replacement is not made, and as written the inequality is false with a uniform constant. A secondary but related typo is the missing factor 2 in (3.30): since J(u_eps) -> pi/(1+g) sum d_j delta and J = (1/2) div(Jj), the weak limit ja of j(u_eps) must satisfy div(J ja) = 2pi/(1+g) sum d_j delta, not pi/(1+g) sum d_j delta. As written this invalidates the identification ja = j(u*) in Lemma 3.2 unless corrected.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":26074,"tokens_out":16720,"duration_ms":156980,"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":[{"comment":"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 ã_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.","section":"§3.2, Lemma 3.1, Eq. (3.21)"},{"comment":"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 δ_{ã_j(t)} in (3.7) implies div(J j_a)=2π/(1+g)∑_j d1_j δ_{ã_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 δ_{ã_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.","section":"§3.3, Lemma 3.2, Eq. (3.30)"}],"minor_comments":[{"comment":"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.","section":"§3.4, Eq. (3.45)"},{"comment":"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.","section":"§3.3, Lemma 3.2 proof"},{"comment":"The displayed identity '(rf2)′ = rf2′ + f2 = rh(r)' is inconsistent with the subsequent integration step; the correct identity is (r f2′)′ = r h(r).","section":"§5, Lemma 5.3, Eq. (5.23)"},{"comment":"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.","section":"§5, Lemma 5.2 proof, around (5.17)"},{"comment":"The notation ‖φ‖_{C²(Ω×[t1,t2])} is misleading because φ is time-independent; this should be ‖φ‖_{C²(Ω)}.","section":"§3.2, Eq. (3.21)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of math.AP and the central claim is plausible. The main obstacle to acceptance is not the strategy but the writing-level gaps: the uniform Lipschitz bound in Lemma 3.1 and the factor 2 in Lemma 3.2 are both repairable, but they are load-bearing. The paper also relies on several nontrivial estimates imported from the author's earlier work [13]; if the editorial standard favors self-containedness, I recommend asking the author to state the identity (3.16) in [13] used in Lemma 3.5. There is no indication of circular reasoning, and the numerical estimate for γ_g is a useful complement."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","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.\n\nWhat 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.\n\nThe 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.\n\nThe abstract also overstates the range of the decoupling by not stating the well-preparedness and separation assumptions, but the theorem itself is honest.\n\nWho 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.","headline":"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.","tokens_in":26600,"tokens_out":4859,"would_cite":true,"duration_ms":69255,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35Q41"],"pacs":[],"model":"deepseek-v4-flash","headline":"The coupled nonlinear Schrödinger equation's vortices decouple in the small-core limit: each component follows the single-component vortex motion law.","keywords":["coupled nonlinear Schrödinger equation","quantized vortex","canonical harmonic map","reduced dynamical law","renormalized energy","fractional vortex","vortex path","Bose-Einstein condensation"],"falsifier":"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.","tokens_in":25453,"feed_emoji":"🌀","tokens_out":8631,"duration_ms":80228,"temperature":0.7,"pith_summary":"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.","feed_headline":"Two-component vortex motion decouples as core size shrinks","feed_subtitle":"Rigorous proof: each component's vortices follow the same renormalized-energy ODE as scalar NLS.","key_machinery":"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.","core_discovery":"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$.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the scalar Ginzburg-Landau vortex-path theory: vortex existence via Theorem 1.4.4, energy bounds via Theorem 1.4.3, and the derivative identity (2.11) that the proof adapts componentwise.","marker":"[3]"},{"why":"Introduces the canonical harmonic map H, refined Jacobian estimates, and the renormalized energy W_d; its asymptotic (2.8) gives the leading energy of the comparison configuration.","marker":"[8]"},{"why":"Provides the refined Jacobian estimates used to control winding-number degrees and locate vortex centers in Lemma 4.1.","marker":"[7]"},{"why":"Gives existence, uniqueness, and decay estimates for the radial ODE governing the vortex profile, used to define gamma_g in Section 5.","marker":"[6]"},{"why":"Supplies identities from scalar NLS vortex dynamics used to convert path-derivative and Hessian identities into the proof of the ODE.","marker":"[13]"}],"fun_headline_variants":["Vortex motion decouples in coupled NLS as core shrinks","Each component's vortices move independently in zero-core limit","Coupled NLS vortices decouple: single-component law emerges","Small-core limit yields independent vortex dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Vortex motion decouples in coupled NLS as core shrinks","Each component's vortices move independently in zero-core limit","Coupled NLS vortices decouple: single-component law emerges","Small-core limit yields independent vortex dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000338,"raw_usage":{"total_tokens":1842,"prompt_tokens":890,"completion_tokens":952,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":893}},"tokens_in":506,"tokens_out":952,"duration_ms":8640,"temperature":1.0,"reasoning_tokens":893,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:56:11.255047+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the scalar Ginzburg-Landau vortex-path theory: vortex existence via Theorem 1.4.4, energy bounds via Theorem 1.4.3, and the derivative identity (2.11) that the proof adapts componentwise."},{"cited_title":"Refined jacobian es timates and gross–pitaevsky vortex dynamics","cited_arxiv_id":null,"evidence_quote":"Introduces the canonical harmonic map H, refined Jacobian estimates, and the renormalized energy W_d; its asymptotic (2.8) gives the leading energy of the comparison configuration."},{"cited_title":"Jerrard and Daniel Spirn","cited_arxiv_id":null,"evidence_quote":"Provides the refined Jacobian estimates used to control winding-number degrees and locate vortex centers in Lemma 4.1."},{"cited_title":"Uniqueness results for an ode related to a generaliz ed ginzburg–landau model for liquid crystals","cited_arxiv_id":null,"evidence_quote":"Gives existence, uniqueness, and decay estimates for the radial ODE governing the vortex profile, used to define gamma_g in Section 5."},{"cited_title":"Quantized vo rtex dynamics of the nonlinear schr¨ odinger equation on tor us with non-vanishing momentum","cited_arxiv_id":null,"evidence_quote":"Supplies identities from scalar NLS vortex dynamics used to convert path-derivative and Hessian identities into the proof of the ODE."}],"review_version":1}