{"id":"fd3b3cf8-9273-4349-8936-a52aaf329f5c","arxiv_id":"2507.10438","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The short-range n-vortex interaction energy in the abelian Higgs model is derived from Jacobi-operator spectral data, with the pair potential scaling as separation to the fourth power.","lead":"Two nearby vortices in the Ginzburg-Landau (abelian Higgs) model feel an interaction energy that scales as the fourth power of their separation, not the square, according to this paper's explicit formula. The result gives a quantitative handle on short-range vortex forces, a regime previously accessible only through expensive numerical simulation.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The posited equality behind (5.11) — Hessian of constrained E_int equals unconstrained Hessian restricted to V(λ) — is verified at only two λ values, with one key curve missing data; if the inequality is strict, the spectral coefficients are not the true interaction Hessian.","rationale":"The paper's headline consequence, the R^4 scaling of the two-vortex interaction energy, is robust: it does not depend on the spectral identification, only on the critical-point argument and the U(1) action on polynomial coefficients. The quantitative formula (5.11), however, reaches beyond that scaling and encodes the curvature of E_int at the coincident point. The authors explicitly flag the Hessian identification as a posit and test it numerically, which is honest and provides real support. But the test is limited: two couplings, two vortex numbers, and the collinear n=3 curve at λ=2 has no data in the small-separation regime due to a disclosed lattice artifact. Because E_int is an infimum, the second derivative of the constrained problem is bounded above by the second derivative along any feasible curve; equality requires the eigenmode curves to be exactly the constrained minimizers to second order. This equality is not derived, so the coefficients in (5.11) are conditional on it. The reader's weakest_assumption identifies exactly this issue, and my independent reading agrees. The verdict should remain CONDITIONAL: no change from the reader's assessment, since the concern is real but the numerical agreement and the robustness of the scaling law justify a conditional rather than a reject decision.","tokens_in":12830,"tokens_out":6335,"duration_ms":79376,"concrete_test":"Run constrained minimization for the n=3 collinear curve p(z)=z^3-R^2z at λ=2 (where Fig. 4 omits R<1) using a finer lattice (h1=0.025) and a zero-pinning constraint that prevents the central zero from spreading; for R=0.2,...,0.9 fit E_int(R)-(E3-3E1) to (c2/2)R^4. If the fitted c2 differs from (f0/b2)^2 Λ2 = -0.0238484 by more than ~5% after finite-size extrapolation, the posited Hessian equality in §5 is false. A confirmatory run at λ=0.8 along the same curve would test a parameter value not used in the paper's calibration.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 5 defines the central coefficient formula by 'it is natural to posit that the Hessian of E_int at 0 coincides with the restriction to V(λ) of the Hessian of the Ginzburg-Landau energy functional.' The derivation up to this point is sound: for the eigenmode curve (5.1), the unconstrained energy has second derivative Λ_k, while E_int along the corresponding polynomial curve has second derivative c_k b_k^2/f_0^2, giving (5.10). The missing link is the equality. For any feasible curve of fields realizing the prescribed zero set, E_int(a(t)) ≤ E(Φ_t), with equality at t=0 and both first derivatives zero, so the true Hessian of the constrained infimum is at most the restriction of the unconstrained Hessian. Thus Λ_k is an upper bound for the coefficient unless the optimal fields are exactly the eigenmode linearization. The numerical checks in §6 cover n=2,3 at λ=0.5 and λ=2 only; for the n=3 collinear case at λ=2, the paper explicitly removes all R<1 data because of a lattice artifact, so the most unstable short-range curve is not tested. No convergence or error-bar data accompany the comparisons. The positivity/negativity pattern and R^4 scaling are robust, but the quantitative coefficients (5.11) rest entirely on this unproved equality.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper derives a short-range expansion for the interaction energy of n nearly coincident vortices in the abelian Higgs (Ginzburg-Landau) model. Using U(1) invariance of the centred configuration space, the authors show in Section 2 that the Hessian of the interaction energy at the cocentred n-vortex must be diagonal in the elementary symmetric polynomial coefficients a_k, yielding E_int(a) = E_n - nE_1 + (1/2) Σ c_k |a_k|^2 + O(|a|^3). The main proposal, developed in Sections 3-5, is to compute the coefficients c_k from the spectrum of the Jacobi operator of the n-vortex: for each k = 2,...,n, c_k = (f_0/b_k)^2 Λ_k, where Λ_k is the eigenvalue of the relevant splitting mode and b_k is its leading coefficient at the origin. This gives formula (5.11), and in particular for n=2 the interaction energy varies as the fourth power of the half-separation R, E_int(R) = E_2 - 2E_1 + (c_2/2) R^4. The coefficients are computed numerically for n=2,3 for couplings 0.1 ≤ λ ≤ 2.5 and compared with direct constrained field-theory simulations at λ = 0.5 and λ = 2, with good visual agreement at small to moderate separations.","tokens_in":12962,"tokens_out":8181,"duration_ms":97113,"significance":"If the central formula is correct, the paper provides a valuable quantitative tool for vortex interactions, reducing a difficult constrained field-theory problem to a sequence of linear ODE spectral problems. The rigorous part of the argument — that by U(1) invariance and criticality the leading short-range interaction is quadratic in the polynomial coefficients, giving R^4 for a vortex pair — is clean and independent of the spectral machinery. The paper also contributes a reproducible shooting scheme with a translation-mode consistency check (|Λ_1| < 10^{-5}), tabulated spectral coefficients, and open data. The potential extension to multicomponent Ginzburg-Landau models for type-1.5 superconductivity and core splitting is an interesting outlook. The main weakness is that the quantitative coefficients rely on an unproved identification between the Hessian of the constrained infimum and the restriction of the unconstrained Hessian, which is only partially tested numerically.","major_comments":[{"comment":"The coefficient formula (5.11) rests on the unproved identification, stated as 'it is natural to posit' in Section 5, that the Hessian of the constrained infimum E_int at the cocentred point coincides with the restriction to V(λ) of the Hessian of the full Ginzburg-Landau energy. The preceding derivation only computes the second derivative of the unconstrained energy along the eigenmode curve (5.1), and the matching to the E_int curve uses that this eigenmode realizes the polynomial deformation (5.8) to first order. For a constrained optimization problem, however, the Hessian of the reduced (infimum) function is generically a Schur complement of the full Hessian, not simply its restriction to the constraint-transverse directions; the equality requires a separate argument (for example, an envelope-theorem calculation showing that Lagrange multipliers vanish at the cocentred solution, or an explicit cancellation of cross-terms between V(λ) and the level-set directions of the divisor map). The paper acknowledges this is a posit and tests it numerically in Section 6, but only for n=2,3 at λ=0.5 and λ=2. For the n=3 collinear case at λ=2, the data with R<1 are removed because of a lattice artifact (Section 6, discussion of Figure 4), so the most short-range regime for the most unstable mode is not tested at all. The authors should either prove the identification or substantially extend the numerical verification, including convergence studies, error bars, and an attempt to access the missing short-range collinear data (e.g., with a different discretization or constraint implementation). Without this, the quantitative coefficients in (5.11) are not established beyond the specific tested cases.","section":"Section 5, Eqs. (5.4)-(5.11)"}],"minor_comments":[{"comment":"The matrix H is used both for the Hermitian quadratic-form matrix and for the Higgs field elsewhere; consider renaming one of them to avoid confusion.","section":"Section 2, Eq. (2.8)"},{"comment":"The sign convention in the small-r expansion of α3, with a minus sign before the e4 term, should be checked against the smoothness conditions (4.13)-(4.14) to ensure no typographical error.","section":"Section 4, Eq. (4.16)"},{"comment":"The backward shooting map S2 is defined on R^3; it would be helpful to state explicitly why the decaying boundary conditions at r2 leave exactly three free parameters (b1,b2,b3), given that the system is first-order in seven variables.","section":"Section 4, after Eq. (4.21)"},{"comment":"The column headers 'cn=2 2', 'cn=3 2', and 'cn=3 2' are garbled; use notation such as c_2^{(2)}, c_2^{(3)}, c_3^{(3)} to distinguish the coefficients for n=2 and n=3.","section":"Section 6, Table 1"},{"comment":"Minor typos: 'coeffcients' should be 'coefficients', and 'colinear' should be 'collinear' in the discussion of the vortex line.","section":"Section 6"},{"comment":"The phrase 'Eint vanishes indentically' contains a typo; it should read 'identically'.","section":"Section 7"},{"comment":"The notation 'the zeros of ϕt(z)' should be typeset with the subscript on ϕ_t to match the surrounding text.","section":"Section 5, Eq. (5.5)"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and the R^4 scaling result is elegant and rigorous. The main concern is the unproved identification in Section 5 that converts spectral data into interaction coefficients; this is a genuinely load-bearing step, and the numerical evidence is limited to two couplings with one important short-range curve missing. If the authors can supply a proof (or a much more extensive numerical verification), the paper would be a strong candidate for acceptance. The current manuscript, while interesting, leaves the central quantitative claim insufficiently supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I'll get straight to it: the paper earns its headline. The observation that the centred n-vortex configuration space has U(1) symmetry acting on the polynomial coefficients a_k with weight k, and that the coincident point is a critical point of E_int, forces the pair interaction to be quartic in R rather than quadratic. That argument is rigorous and does not depend on the spectral machinery. The paper is also genuinely useful: it gives the first quantitative short-range theory in the abelian Higgs model, and the coefficient data are deposited with the DOI, so the numbers are checkable. The shooting method for the Jacobi operator looks well posed, and the translation-mode consistency check |Lambda_1| < 10^{-5} is a good sign.\n\nThe soft spot is exactly where the stress-test placed it: equation (5.11) identifies the Hessian of the constrained infimum with the restriction of the unconstrained Hessian to the near-zero eigenspace V(lambda). That is a real assumption, not a theorem. For any feasible curve realizing the prescribed zeros, E_int(a(t)) <= E(Phi_t), so the true Hessian is at most the restricted Hessian; the coefficients could in principle be smaller in magnitude. The authors flag this, and they test it in section 6 at lambda=0.5 and lambda=2 for n=2,3, finding good agreement. That is reasonable support. But the n=3 collinear curve at lambda=2, which is exactly the most unstable short-range case, has no data for R<1 because of a lattice artifact. So for the case where the assumption would be most severely tested, the paper cannot test it. I would not call that fatal — the agreement elsewhere is encouraging — but the authors should either prove the saturation condition or extend the numerics before I'd call the coefficients definitive.\n\nThe paper is a serious piece of work: careful symmetry arguments, explicit ODE reduction, independent verification against full field theory with no fitted parameters. The citation pattern is fine; the new result is clearly distinguished from prior work on nonpairwise interactions and long-range forces. It deserves a full peer review. I would suggest the referee ask for convergence data for the spectral and lattice computations, more lambda values for the identification check, and ideally some discussion of when the variational inequality is saturated. But this is a strong candidate for acceptance, not a desk reject.","headline":"The R^4 pair-force scaling is rigorous and elegant, but the quantitative coefficient formula rests on an unproved variational identification that the authors verify at only two couplings.","tokens_in":13640,"tokens_out":1830,"would_cite":true,"duration_ms":20188,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T13","35Q51"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper derives an explicit short-range formula for the interaction energy of $n$ coincident abelian Higgs vortices, showing that a vortex pair's energy varies as the fourth power of half-separation, with coefficients computed from…","keywords":["abelian Higgs vortices","Ginzburg-Landau model","interaction energy","Jacobi operator","short-range forces","spectral methods","topological solitons","type-1.5 superconductivity"],"falsifier":"Compute $E_{\\mathrm{int}}$ for two vortices at very small half-separations $R$ (for example $0.02$ to $0.3$) using high-resolution constrained minimisation, and fit $E_{\\mathrm{int}}(R)-E_2+2E_1$ to a power law; if the leading exponent is not $4$ within discretisation error, or if the fitted coefficient disagrees with $(f_0/b_2)^2\\Lambda_2$ beyond numerical tolerance, the Section 5 identification is wrong.","tokens_in":2219,"feed_emoji":"🌀","tokens_out":2408,"duration_ms":94800,"temperature":0.7,"pith_summary":"In the abelian Higgs (Ginzburg-Landau) model, vortices exert forces on each other, but until now only the long-range part of those forces was quantitatively understood. This paper aims to establish the short-range counterpart: when $n$ vortices nearly coincide, the interaction energy is a quadratic form in the coefficients of the polynomial whose roots are the vortex positions, with coefficients fixed by spectral data of the Jacobi operator of the cocentred $n$-vortex. The immediate payoff is that a vortex pair separated by distance $2R$ has interaction energy $(c_2/2)R^4$ at short range, quartic rather than quadratic in $R$. The coefficients are computed for $n=2,3$ over couplings $0.1\\le\\lambda\\le2.5$ and matched against full field theory simulations at $\\lambda=0.5$ and $\\lambda=2$, giving quantitative short-range potentials that were previously missing.","feed_headline":"Vortex pairs attract or repel with the fourth power of separation","feed_subtitle":"A single spectral computation reproduces two- and three-vortex attraction and repulsion in full field simulations.","key_machinery":"The central object is the Jacobi operator $J$ of the cocentred $n$-vortex, the Hessian of the Ginzburg-Landau energy at that solution. Rotational equivariance splits $J$ into ODE operators $J_k^+$ on subspaces $C_k^+$ labelled by integer $k$; the paper needs the splitting modes, the $2n$ eigenmodes whose eigenvalues $\\Lambda_k$ pass through $0$ at critical coupling $\\lambda=1$ (one for each $k=1,\\dots,n$ plus a degenerate partner). For each mode, $b_k$ is the leading coefficient in the small-radius expansion of the Higgs-field perturbation, and $f_0$ is the leading coefficient of the unperturbed $n$-vortex profile. The load-bearing identification, posited in Section 5, is that the Hessian of the constrained interaction energy $E_{\\mathrm{int}}$ at the coincident point equals the restriction of the unconstrained energy Hessian to the space $V(\\lambda)$ of splitting modes; this converts the spectral pair $(\\Lambda_k,b_k)$ into the interaction coefficient $c_k$. The polynomial coefficients $a_k$ provide the natural coordinates on the centred $n$-vortex configuration space, and make the $R^4$ law transparent because $a_2=R^2$ for a pair.","core_discovery":"For $a=(a_2,\\dots,a_n)$ the complex coefficients of the centred monic polynomial $p(z)=z^n+a_2z^{n-2}+\\cdots+a_n$ whose roots are the $n$ vortex positions, the paper's short-range formula is $E_{\\mathrm{int}}^{(0)}(a)=E_n-nE_1+\\frac{f_0^2}{2}\\sum_{k=2}^n \\frac{\\Lambda_k}{b_k^2}|a_k|^2$, where $E_n$ is the cocentred $n$-vortex energy, $f_0$ is the leading coefficient of the vortex profile near the origin, and $(\\Lambda_k,b_k)$ are spectral data of the Jacobi operator: $\\Lambda_k$ is the eigenvalue of the splitting mode in the $k$-th rotational symmetry class and $b_k$ is the leading coefficient of its $L^2$-normalised eigenmode. A $U(1)$ rotation symmetry forces the Hessian to be diagonal, eliminating linear terms and cross terms; each coefficient $c_k=(f_0/b_k)^2\\Lambda_k$ is separately computable from a linear ODE problem. The paper computes these coefficients numerically for $n=2,3$ and $\\lambda\\in[0.1,2.5]$, and verifies the resulting potentials against constrained energy minimisation in the full field theory at $\\lambda=0.5$ and $\\lambda=2$, finding close agreement from coincidence out to separations where the potential crosses the well-known long-range formula. For two vortices, $a_2=R^2$, so $E_{\\mathrm{int}}(R)=E_2-2E_1+\\frac{c_2}{2}R^4$.","pith_inferences":["Editorial inference: if the Section 5 identification holds generally, then any soliton model with a critical coupling where splitting modes become zero modes should show a quartic (or higher even) leading interaction law, making the $R^4$ behaviour a generic degeneracy phenomenon rather than a peculiarity of abelian Higgs vortices.","Editorial inference: because $c_k(\\lambda)$ passes through zero at $\\lambda=1$, the short-range potential is extremely flat near coincidence near critical coupling, so in the crossover region the next Taylor order may dominate; measuring that order would test whether the splitting-mode subspace captures the full Hessian.","Editorial inference: the same spectral construction could be applied to $n=4$; the one-parameter-per-mode structure gives an independent cross-check, since a two-vortex measurement fixes $c_2$ and a collinear three-vortex measurement fixes $c_3$ without extra fitting parameters.","Editorial inference: the static potentials computed here provide a baseline for interpreting recent results on interaction forces between vortices whose shape modes are excited, where fluctuation-induced attractive forces appear even at critical coupling."],"forward_implications":["For two vortices at half-separation $R$, the short-range interaction energy is $E_2-2E_1+(c_2/2)R^4$, so the leading force is quartic, not quadratic, in separation.","Short-range interactions are not a sum of pairwise terms for $n\\ge3$: three collinear vortices at spacing $R$ give quartic dependence, while an equilateral triangle of side $R$ gives sextic dependence.","The coefficients $c_k(\\lambda)$ are tabulated for $0.1\\le\\lambda\\le2.5$, and the short-range formula overlaps the long-range point-source formula near $R\\approx3$, allowing the two approximations to be spliced into a global explicit potential.","Computing the short-range interaction energy reduces from constrained energy minimisation in the full field theory to solving linear ODE eigenvalue problems for the cocentred $n$-vortex.","In multicomponent Ginzburg-Landau models, the same spectral criterion (a negative $c_2$ together with a magnetic long-range length scale) flags type-1.5 superconductivity, and a negative core-splitting mode of the one-vortex Jacobi operator flags vortex core splitting."],"supporting_citations":[{"why":"Supplies the symmetry reduction of vortex normal modes that the paper rederives and extends, and identifies the shape mode used in later discussion.","marker":"[2]"},{"why":"Provides the previous spectral decomposition of fluctuations around $n$-vortices that motivates the space $V(\\lambda)$ of splitting modes.","marker":"[3]"},{"why":"Original study of the vortex Jacobi operator's spectrum, giving the foundational eigenvalue framework the paper builds on.","marker":"[8]"},{"why":"Supplies the long-range point-source formula and the vortex monopole and dipole charges used for comparison and for splicing the approximations.","marker":"[13]"},{"why":"Provides the constrained energy-minimisation algorithm used for the full field theory numerical interaction energies.","marker":"[15]"},{"why":"Makes the computed spectral coefficients available for practical use and for the proposed splined global approximation.","marker":"[16]"},{"why":"Demonstrates nonpairwise intervortex interactions in type-1.5 superconductors, motivating the need for a systematic short-range formula rather than a pairwise sum.","marker":"[5]"}],"fun_headline_variants":["Vortex interaction energy scales as d^4 at short range","Close vortices: fourth power force law emerges","Short-range vortex forces: d^4, not d^2","Spectral method predicts vortex interaction at close distances","Fourth power controls short-range vortex pairs"],"cache_read_input_tokens":15616,"weakest_assumption_plain":"The load-bearing premise is that the Hessian of the constrained interaction energy at the coincident-vortex point equals the second variation of the unconstrained Ginzburg-Landau energy along the near-zero splitting modes of the Jacobi operator, an identification the paper posits and tests numerically but does not derive from the constrained variational problem.","fun_headline_variants_meta":{"raw":{"variants":["Vortex interaction energy scales as d^4 at short range","Close vortices: fourth power force law emerges","Short-range vortex forces: d^4, not d^2","Spectral method predicts vortex interaction at close distances","Fourth power controls short-range vortex pairs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000349,"raw_usage":{"total_tokens":1963,"prompt_tokens":1057,"completion_tokens":906,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":673,"completion_tokens_details":{"reasoning_tokens":830}},"tokens_in":673,"tokens_out":906,"duration_ms":9018,"temperature":1.0,"reasoning_tokens":830,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T17:34:02.633315+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $E_{\\mathrm{int}}$ for two vortices at very small half-separations $R$ (for example $0.02$ to $0.3$) using high-resolution constrained minimisation, and fit $E_{\\mathrm{int}}(R)-E_2+2E_1$ to a power law; if the leading exponent is not $4$ within discretisation error, or if the fitted coefficient disagrees with $(f_0/b_2)^2\\Lambda_2$ beyond numerical tolerance, the Section 5 identification is wrong.","supporting_citations":[{"cited_title":"Dissecting normal modes of vibration on vortices in Ginzburg-Landau superconductors","cited_arxiv_id":null,"evidence_quote":"Supplies the symmetry reduction of vortex normal modes that the paper rederives and extends, and identifies the shape mode used in later discussion."},{"cited_title":"Spectral structure of fluc- tuations around n-vortices in the abelian-higgs model","cited_arxiv_id":null,"evidence_quote":"Provides the previous spectral decomposition of fluctuations around $n$-vortices that motivates the space $V(\\lambda)$ of splitting modes."},{"cited_title":"Bound states and instabilities of vortices","cited_arxiv_id":null,"evidence_quote":"Original study of the vortex Jacobi operator's spectrum, giving the foundational eigenvalue framework the paper builds on."},{"cited_title":"Static intervortex forces","cited_arxiv_id":null,"evidence_quote":"Supplies the long-range point-source formula and the vortex monopole and dipole charges used for comparison and for splicing the approximations."},{"cited_title":"Intervortex forces in competing-order superconductors","cited_arxiv_id":null,"evidence_quote":"Provides the constrained energy-minimisation algorithm used for the full field theory numerical interaction energies."},{"cited_title":"Short range intervortex force coefficients","cited_arxiv_id":null,"evidence_quote":"Makes the computed spectral coefficients available for practical use and for the proposed splined global approximation."},{"cited_title":"Semi-meissner state and nonpairwise intervortex interactions in type-1.5 superconductors","cited_arxiv_id":null,"evidence_quote":"Demonstrates nonpairwise intervortex interactions in type-1.5 superconductors, motivating the need for a systematic short-range formula rather than a pairwise sum."}],"review_version":1}