{"id":"3dccb7d2-3e46-4909-9ce1-606c6daa7011","arxiv_id":"2507.18590","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the planar Gross-Pitaevskii equation, explicit n-vortex solutions are constructed whose dynamics follow the Helmholtz-Kirchhoff system at leading order, with the first correction given by a linear wave equation.","lead":"This paper proves that a fundamental quantum physics equation has solutions made of several tiny spinning vortices that move along the classical point-vortex paths predicted by Neu in 1990. It also derives, rigorously for the first time, the small wave-induced correction to those paths, confirming a later formal prediction of Ovchinnikov and Sigal.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The construction's inner-correction stage depends on the companion-paper solvability theory for the vortex linearized operator, and the in-paper text delegates the general Fourier-mode extension to [17]; if that theory fails or does not cover the epsilon-cutoff RHS used here, Theorems 1-3 collapse.","rationale":"The reader's weakest-assumption analysis correctly identifies the companion-paper solvability theory as the single most load-bearing unresolved input. Good-faith reading: the paper is a serious constructive PDE argument; the approximation and energy-estimate architecture are coherent, and the wave-equation estimates in Section 6 appear internally consistent. The central claim, however, requires solving L_j[phi_j]=H_j for general right-hand sides with quantitative decay/growth and derivative bounds. That theory is not proved in the manuscript: Section 3.2 gives finite-mode representation formulae and then delegates the general case to [17] via a contradiction/maximum-principle argument. Every subsequent conclusion--the O(epsilon^3|log epsilon|^2) error after first improvement, the arbitrary-order approximation in Proposition 5.1, the energy coercivity in Proposition 7.9, and finally Theorems 1-3--inherits this dependency. The concrete check proposed would settle whether [17]'s published hypotheses cover the epsilon-dependent RHS's used here, especially the cut-off dependence on delta epsilon^{-1} and the logarithmic growth estimates. If the check passes, the conditional verdict can be upgraded; if it fails, the construction has a real gap. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":64557,"tokens_out":37448,"duration_ms":373065,"concrete_test":"Check the general-mode extension in [17, Section 5] against the RHS appearing in Lemma 4.11 (equations (4.30)-(4.31)) and Lemma 4.13. Specifically, write the full RHS (after the cut-off eta~^{(9/4)}_j) in the Fourier decomposition (3.7)-(3.8), apply the representation formula (3.18) mode by mode, and verify that the infinite sum converges and yields the uniform bounds (4.32) and (5.21) with constants independent of the cut-off radius delta epsilon^{-1}. In particular, confirm that the case Re(h)=O(r^{-2}), Im(h)=O(1) stated in Remark 3.4 accommodates the compactly supported h with |h_2| <= C epsilon^3 r log r inside B_{C epsilon^{-1}}, and that the rescaled Schauder argument in Remark 3.5 produces the D^ell_y estimates in Remark 4.12 without losing powers of epsilon.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The approximate solution u* is produced in Proposition 4.10 and iterated in Proposition 5.1 by solving, near each vortex, elliptic equations of the form L_j[phi_j]=H_j (Sections 4-5). The paper only proves a solvability statement for right-hand sides with finitely many Fourier modes (the representation formulae in Proposition 3.3 and Remark 3.4); the extension to general right-hand sides, including the r log r growth bounds in Lemma 4.11 and Lemma 4.13, is explicitly delegated: 'Existence and estimates extend to the general case ... using a contradiction argument and maximum principle type estimates. See [17, Section 5]' (Remark 3.4). This is the load-bearing step: the inner corrections phi_j determine the zero locations and feed the outer wave equation for psi^{out,1}, which in turn produces the O(epsilon^2|log epsilon|^2) trajectory correction and the modified dynamics (1.11). If [17]'s hypotheses (stated there for Re(h)=O(r^{-2}), Im(h)=O(1)) do not cover the RHS's used here--supported on balls of radius C delta epsilon^{-1} with imaginary parts of size epsilon^3 r log r--or if the derivative bounds asserted in Remark 3.5 do not follow with constants independent of epsilon, then the estimates (4.32), (5.21), and the higher-derivative versions used in Proposition 5.1 fail. The paper itself flags the omission, so this is not a manufactured concern. A secondary, related gap is that Proposition 5.1's induction and the coercivity argument in Proposition 7.9 are only summarized; however, the companion-paper dependency is the most direct single point of failure.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs, for any n >= 2, smooth solutions of the planar Gross-Pitaevskii equation that look like products of n degree +/-1 vortices with trajectories close to a given collisionless solution of the Helmholtz-Kirchhoff system. The main results are Theorem 1 (existence with O(epsilon^2 |log epsilon|^2) trajectory and remainder error), Theorem 2 (an asymptotic expansion of the phase in the outer region whose first correction is the solution of a linear wave equation), and Theorem 3 (locations of the zeroes up to O(epsilon^{3-sigma}), governed by a modified system with a forcing term equal to 2 grad psi_1^{out,1}). The proof is a gluing construction: elliptic corrections are built in the vortex cores using the linearized Ginzburg-Landau operator, wave corrections are built in the outer region, the parameters xi_j(t) are adjusted to kill resonant modes, and a final energy/quadratic-form argument controls the remainder. The paper is a 65-page manuscript with explicit error bookkeeping and an induction producing an approximation accurate to arbitrary algebraic order in epsilon.","tokens_in":64897,"tokens_out":9999,"duration_ms":102689,"significance":"If the proofs are complete, this is a major advance: it supplies the first rigorous multi-vortex solutions whose vortex trajectories are tracked beyond the leading-order Helmholtz-Kirchhoff law, and it gives the first rigorous justification of the Ovchinnikov-Sigal radiation correction. The paper is genuinely constructive, contains no fitted parameters, and makes precise falsifiable predictions about the profile and the first-order dynamics. The introduction of a tailored quadratic form for the linearized evolution is a promising methodological contribution. However, the strength of the paper rests on a small number of load-bearing black boxes, especially the solvability theory for the linearized vortex operator, which is delegated to the companion paper [17]. The manuscript also contains a possible coefficient error in the derivation of the wave equation for the outer phase correction that affects the statements of Theorems 2 and 3.","major_comments":[{"comment":"The general solvability theory for the linearized Ginzburg-Landau operator L_j is delegated to the companion paper [17, Section 5]. The manuscript proves sharp representation formulas only for right-hand sides with finitely many Fourier modes, but the applications in Lemmas 4.11, 4.13 and the recursive problems (5.19)-(5.20) require existence and weighted estimates for right-hand sides of epsilon-dependent support and with imaginary parts of size epsilon^3 r log r, together with derivative bounds uniform in epsilon. Remark 3.4 states the hypotheses of [17] as Re(h)=O(r^{-2}) and Im(h)=O(1), but the manuscript does not verify that these hypotheses cover the present right-hand sides. This black box underpins estimates (4.32), (5.21) and hence Proposition 5.1; if it fails or does not apply, the inner corrections, the zero locations, and the outer wave equation all collapse. Please either state and prove the required extension in this paper, or give a precise verifiable correspondence between the hypotheses of [17] and every right-hand side used here.","section":"Section 3.2, Remark 3.4"},{"comment":"The induction proving Proposition 5.1 is presented as an outline. In particular, the key derivative-loss estimate (5.45) is asserted rather than proved, the recursive definition of M_k is only sketched, and the constants c_k are not tracked explicitly. Since Proposition 5.1 produces the arbitrary-order approximate solution that is fed into the final error bound (8.1), this summarized induction is load-bearing. The paper would need either a complete induction proof or a precise induction statement with all constants and a detailed verification of the derivative bookkeeping before the claim can be accepted.","section":"Section 5.3, Proof of Proposition 5.1"},{"comment":"There appears to be a factor error in the definition of the wave forcing. With tau = sqrt(2) epsilon^{-1} t, substituting psi_2^{out} = (1/2)(E_2 + epsilon^2 partial_t psi_1^{out}) into the real part of (4.39) gives a term -epsilon^{-1}/sqrt(2) partial_tau E_2 in the rescaled equation, not -epsilon^{-1} sqrt(2) partial_tau E_2 as written in (4.42) and (4.45). As a consequence, the forcing F^{out,1} used in the wave equation, in Theorem 2, and in the modified dynamics (1.11) may be off by a factor of two. Please check the algebra carefully; if the displayed factor is intentional, explain the rescaled variables; otherwise correct the coefficient and re-verify the estimates in Lemma 4.15 that depend on the size of the solution of (4.46).","section":"Section 4.3, equations (4.42)-(4.45)"},{"comment":"The coercivity proof of the quadratic form is a contradiction/compactness argument and is only partially written out. In Step 3, the passage from the limit profile, the use of [16, Lemma 3.1] to obtain (7.49), and the treatment of boundary terms at |y_j|=R are summarized rather than demonstrated, and the uniformity in epsilon of the compactness is not shown. Since the coercivity estimate (7.55) is the basis of Proposition 7.13 and of the final linear estimate (7.2), this is another load-bearing point that needs to be either fully proved or replaced by a precise reference with verified hypotheses.","section":"Section 7.3, Proposition 7.9"}],"minor_comments":[{"comment":"The phrase 'arbitrarily large, finite time interval' could be misread as uniform in T; the theorems fix T and all constants depend on T. Please clarify that the interval is arbitrary but fixed before epsilon is chosen.","section":"Abstract and Section 1.1"},{"comment":"The computation of d/dt B[phi,phi] is stated with 'the details left to the reader.' Since this identity is central to the energy argument, include at least the main integration-by-parts steps or a supplementary calculation in an appendix.","section":"Section 7.1, Lemma 7.2"},{"comment":"The constants c_k=c_{k,m} in Proposition 5.1 are never assigned explicit values; give the recursive definition or at least state that they are finite and independent of epsilon, with the dependence on m and k made explicit.","section":"Section 5.3"},{"comment":"The notation in equation (1.10) and (4.45) uses O_c(epsilon^2 |log epsilon|), 'lower order terms', and a cut-off chi that are introduced informally in Theorem 2. Please give the exact definition of chi and the compactly supported convention in the theorem statement, not only in Section 4.","section":"Section 4.3"}],"recommendation":"major_revision","confidential_remarks":"This is a serious and potentially landmark paper, but its acceptance depends on content that is not fully contained in the manuscript: the elliptic solvability theory of [17], the complete induction in Proposition 5.1, and the correctness of the factor in (4.42)-(4.45). The factor discrepancy is concrete enough that it should be resolved before publication. If the authors can supply the missing verification and correct the factor, the paper is likely publishable in a top analysis journal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this is the first construction of genuine time-dependent multi-vortex solutions for the planar Gross-Pitaevskii equation with the field profile near the cores and the vortex trajectories tracked through the first correction. If the companion solvability paper [17] holds up, it answers the Pacard–Riviere question and makes the formal Neu and Ovchinnikov–Sigal expansions rigorous.\n\nWhat is new: earlier rigorous results only showed that the vortex Jacobian converges to Dirac masses moving by (1.5); they described neither the solution near the cores nor any correction to the motion. The elliptic constructions were stationary. This paper builds an n-vortex ansatz and, in Theorem 3, obtains the corrected trajectory from the gradient of a linear wave solution, with parameters frozen before solving the wave equation. There is no circular fitting and no parameter fitting.\n\nWhat is done well: the paper is honest about its structure. The error bookkeeping is detailed, the wave estimates in Section 6 are substantive, and the text explicitly flags where arguments are delegated. The self-citation to [17] is appropriate: the linearized Ginzburg-Landau theory really is the load-bearing black box, and the authors say so in Remark 3.4.\n\nThe soft spots are real but proportional. Remark 3.4 sends the extension from finitely many Fourier modes to general right-hand sides to [17]. Lemma 4.11, Lemma 4.13, and Proposition 5.1 all use that extension, and if [17] does not cover the epsilon-dependent cutoffs and r log r growth used here, Theorems 1–3 collapse. That is not a manufactured concern; the paper itself says the proof is in [17]. Lemma 4.13 and the induction behind Proposition 5.1 are also sketched rather than fully written, and the coercivity proof in Proposition 7.9 is a contradiction argument with several steps compressed. These are normal for a 65-page gluing paper, but they make independent verification slow.\n\nOn balance, I would not bet against the construction, but I also would not sign off without seeing [17] or having the general solvability argument written into this paper. The architecture is coherent, the novelty is real, and the citation pattern is clean.\n\nWho it is for: researchers in Gross-Pitaevskii vortex dynamics and singular perturbation/gluing methods. It deserves a serious referee; desk rejection would be wrong. Recommended course: send to peer review, ask for an independent check of the companion solvability theory, and request expansion of Lemma 4.13 and the induction in Proposition 5.1. If those check out, this is a significant advance.","headline":"This is the first construction of true time-dependent multi-vortex solutions for planar Gross-Pitaevskii with tracked profiles and a first-order correction to Helmholtz-Kirchhoff dynamics; the main risk is the load-bearing elliptic solvability theorem delegated to the companion paper [17].","tokens_in":65461,"tokens_out":2058,"would_cite":true,"duration_ms":24425,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q55","35B25","35C20","76B47","35L05"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper constructs smooth multi-vortex solutions of the planar Gross-Pitaevskii equation whose trajectories follow the Helmholtz-Kirchhoff point-vortex law to leading order, with a first correction fixed by a linear wave equation.","keywords":["Gross-Pitaevskii equation","vortex dynamics","Helmholtz-Kirchhoff system","Ginzburg-Landau vortex","matched asymptotics","radiation correction","linear wave equation","singular perturbation"],"falsifier":"Solve the modified system (1.11) for a symmetric pair of equal-degree vortices, compute the resulting first-order correction to the angular frequency and separation, and compare these predictions to direct numerical evolution of (1.1) at several small $\\epsilon$ (for example $\\epsilon=0.02,0.01,0.005$) over one rotation period. If the measured displacement of the zeroes from $\\xi^0(t)$ is not $O(\\epsilon^2|\\log\\epsilon|^2)$ with the wave-determined coefficient, or if the $\\epsilon^{3-\\sigma}$ zero-location statement of Theorem 3 fails, the central claim collapses.","tokens_in":64306,"feed_emoji":"🌀","tokens_out":7669,"duration_ms":80201,"temperature":0.7,"pith_summary":"This paper proves that the formal vortex-dynamics asymptotics of Neu are real: for any collisionless trajectory of the Helmholtz-Kirchhoff point-vortex system and any sufficiently small core size $\\epsilon$, there is a smooth solution of the planar Gross-Pitaevskii equation made of sharply rescaled degree-one vortices whose positions track that trajectory to within $O(\\epsilon^2|\\log \\epsilon|^2)$. The construction works for any number $n\\geq 2$ of vortices of degree $\\pm 1$ on an arbitrarily long finite time interval, and it is the first to give the solution profile near the cores rather than only a measure-theoretic limit of vorticity. The paper also identifies the first correction to the leading-order motion: a small forcing term obtained from a linear wave equation, which the vortices themselves emit as they move, validating the formal expansion of Ovchinnikov and Sigal. A sympathetic reader should care because the result turns a 35-year-old formal picture into theorems with explicit descriptions of both the zero set and the next-order dynamics.","feed_headline":"Vortices follow point-vortex law, and the fix is a wave","feed_subtitle":"First explicit multi-vortex solutions show the leading motion and its radiation correction.","key_machinery":"The carrying object is a matched two-scale ansatz, $u_*(y,t)=e^{i\\psi^{\\mathrm{out}}}\\prod_j(\\eta_j(W_j+\\varphi_j)+(1-\\eta_j)W_j e^{\\varphi_j/W_j})$, in which inner corrections $\\varphi_j$ near each vortex solve elliptic equations $L_j[\\varphi_j]=H_j$ for the linearized Ginzburg-Landau operator, and the outer correction $\\psi^{\\mathrm{out}}_1$ solves the linear wave equation (1.8). The elliptic step fixes the vortex parameters through orthogonality conditions; the wave step turns the far-field error into radiation. A separate linear stability analysis controls the remainder by a modified energy: a quadratic form built from the second variation of a conserved-type functional, with a vector field $A$ and scalar $B$ chosen so that the formal continuity and phase equations (2.20)-(2.21) hold, yielding coercivity on the complement of the approximate kernel.","core_discovery":"The central claim is Theorem 1 together with its refinements in Theorems 2 and 3. Given a smooth collisionless solution $\\xi^0(t)$ of (1.5) and $\\epsilon$ sufficiently small, there exists a smooth solution $u_\\epsilon$ of (1.1) of the form (1.7), a product of rescaled copies of the degree-one vortex and its conjugate at positions $\\xi_j(t)=\\xi^0_j(t)+O(\\epsilon^2|\\log\\epsilon|^2)$, with a remainder of size $O(\\epsilon^2|\\log\\epsilon|^2)$. Away from the cores, the phase admits the expansion (1.9) in which the first non-trivial term $\\psi^{\\mathrm{out},1}_1$ solves the linear wave equation (1.8) with forcing (1.10); and the zeroes of $u_\\epsilon$ lie at points $\\xi^*_j(t)+O(\\epsilon^{3-\\sigma})$ where $\\xi^*$ solves the modified system (1.11). In short, the paper establishes both Neu's leading-order dynamics and the Ovchinnikov-Sigal first-order radiation correction by constructing solutions rather than by passing to limits.","pith_inferences":["Editorial inference: a similar inner-outer decomposition with a radiation wave should apply to other soliton equations with topological charges, such as wave-map or complex Ginzburg-Landau models; the obstacle will be finding the analogue of the coercive quadratic form.","Editorial inference: the explicit forcing in (1.11) gives a concrete route to test the predicted slow separation of two same-sign vortices; one can compute the radiation-induced drift for the two-vortex configuration and compare it with the $O(t^{1/6})$ spiral predicted by formal asymptotics.","Editorial inference: a direct numerical experiment is now well posed: evolve (1.1) from initial data built from the ansatz (2.7), measure the zeroes at several $\\epsilon$, and check that the displacement from $\\xi^0(t)$ follows the $O(\\epsilon^2|\\log\\epsilon|^2)$ scale with the wave-determined coefficient; agreement would confirm that the detected correction is the dominant next-order effect."],"forward_implications":["For every $n\\geq 2$, every smooth collisionless point-vortex trajectory is realized by an actual family of Gross-Pitaevskii solutions, so the Helmholtz-Kirchhoff law is a theorem about smooth solutions, not only about limiting vorticity measures.","The vortex positions of the constructed solutions are known pointwise to order $\\epsilon^{3-\\sigma}$, which makes the zero set a controlled object rather than a weak limit.","The first correction to vortex motion is explicitly computable from a linear wave equation, so the Ovchinnikov-Sigal radiation term becomes a defined, checkable quantity rather than a formal one.","The same construction yields refined remainder estimates in weighted and unweighted $H^1$, giving quantitative control on how close the true solution remains to a product of rescaled vortices over the whole time interval."],"supporting_citations":[{"why":"supplies the formal matched-asymptotics vortex dynamics that the paper turns into theorems.","marker":"[38]"},{"why":"gives the formal first-order radiation correction to vortex dynamics that Theorems 2 and 3 rigorously detect.","marker":"[40]"},{"why":"provides the solvability theory for the linearized Ginzburg-Landau operator around the degree-one vortex, the main external input used in the inner corrections.","marker":"[17]"},{"why":"supplies the stationary multi-vortex construction and kernel classification that the dynamic ansatz extends.","marker":"[42]"},{"why":"proves minimality and non-degeneracy of the degree-one vortex, used for coercivity of the quadratic form.","marker":"[16]"},{"why":"states the earlier measure-theoretic vortex-motion limit that the constructed solutions sharpen into pointwise trajectory statements.","marker":"[3]"},{"why":"provides quantitative Jacobian estimates and the hydrodynamic-limit context that motivate the leading-order law the paper realizes.","marker":"[26]"}],"fun_headline_variants":["Vortex dynamics: Kirchhoff law plus wave correction","Rigorous vortex motion: Kirchhoff plus wave correction","Vortex cores follow Kirchhoff, with a wave correction","Rigorous multi-vortex dynamics: Kirchhoff and wave"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The construction relies on a technical result it does not prove in this paper — the solvability of the linearized vortex equation for general right-hand sides, stated in Section 3.2 and attributed to the companion paper [17] — and if that result is incomplete, the approximate profile, phase expansion, and corrected dynamics all fail.","fun_headline_variants_meta":{"raw":{"variants":["Vortex dynamics: Kirchhoff law plus wave correction","Rigorous vortex motion: Kirchhoff plus wave correction","Vortex cores follow Kirchhoff, with a wave correction","Rigorous multi-vortex dynamics: Kirchhoff and wave"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000827,"raw_usage":{"total_tokens":3600,"prompt_tokens":917,"completion_tokens":2683,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":533,"completion_tokens_details":{"reasoning_tokens":2616}},"tokens_in":533,"tokens_out":2683,"duration_ms":18376,"temperature":1.0,"reasoning_tokens":2616,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:10:52.549994+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Solve the modified system (1.11) for a symmetric pair of equal-degree vortices, compute the resulting first-order correction to the angular frequency and separation, and compare these predictions to direct numerical evolution of (1.1) at several small $\\epsilon$ (for example $\\epsilon=0.02,0.01,0.005$) over one rotation period. If the measured displacement of the zeroes from $\\xi^0(t)$ is not $O(\\epsilon^2|\\log\\epsilon|^2)$ with the wave-determined coefficient, or if the $\\epsilon^{3-\\sigma}$ zero-location statement of Theorem 3 fails, the central claim collapses.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"supplies the formal matched-asymptotics vortex dynamics that the paper turns into theorems."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"gives the formal first-order radiation correction to vortex dynamics that Theorems 2 and 3 rigorously detect."},{"cited_title":"del Pino, R","cited_arxiv_id":null,"evidence_quote":"provides the solvability theory for the linearized Ginzburg-Landau operator around the degree-one vortex, the main external input used in the inner corrections."},{"cited_title":"Pacard and T","cited_arxiv_id":null,"evidence_quote":"supplies the stationary multi-vortex construction and kernel classification that the dynamic ansatz extends."},{"cited_title":"del Pino, P","cited_arxiv_id":null,"evidence_quote":"proves minimality and non-degeneracy of the degree-one vortex, used for coercivity of the quadratic form."},{"cited_title":"Bethuel, R","cited_arxiv_id":null,"evidence_quote":"states the earlier measure-theoretic vortex-motion limit that the constructed solutions sharpen into pointwise trajectory statements."},{"cited_title":"Partial Differential Equations 40 (2015), 135–190","cited_arxiv_id":null,"evidence_quote":"provides quantitative Jacobian estimates and the hydrodynamic-limit context that motivate the leading-order law the paper realizes."}],"review_version":2}