{"id":"45bc85bf-006f-4ffe-a41d-cfaf8fbc6110","arxiv_id":"2505.05039","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For non-BPS Abelian-Higgs vortices, the normal modes are classified into Derrick-type, Type A, and Type B multipolar families, with k=1 Type A modes as translational zero modes and k>1 Type A modes as the instability modes for lambda>1.","lead":"Higher-charge Abelian-Higgs vortices have a richer internal vibration spectrum than previously catalogued, with modes that either preserve or break the vortex's rotational symmetry. The authors derive reduced one-dimensional equations for all these modes and compute them numerically for charges one through five.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The angular ansatz (31)-(32) is shown to be compatible with H+, not exhaustive; the claimed completeness of the internal-mode classification rests on an unproven #chθ≤2 reduction.","rationale":"The reader's weakest assumption correctly identifies the main load-bearing gap: the Lemma proves that the chosen two-angle ansatz can satisfy the spectral equation, not that all solutions are captured by it. This is precisely the point on which the Theorem's completeness claim depends. The concern is not a manufactured one: the paper's own language in the Lemma ('admissible eigenfunctions', 'compatible') concedes that exhaustiveness is not proven, and the numerical method solves only the reduced operators, so it cannot falsify the presence of other angular sectors. Direct 2D diagonalization of H+ is the natural arbiter, since it does not presuppose the angular reduction. The BPS-limit checks and the explicit translational modes are genuine supporting evidence, but they constrain only λ=1 and the k=1 sector. The minor contradiction about Type A eigenvalues being 'always positive' is real but secondary; it does not change the conditional verdict. Because the reader already assigned CONDITIONAL for the same reason, no verdict adjustment is needed.","tokens_in":26406,"tokens_out":9521,"duration_ms":115700,"concrete_test":"Compute the low-lying spectrum of the original four-component operator H+ of Eq. (15) directly in two dimensions, e.g. on a polar grid with N=2000 radial points, rmax=20, and a Fourier decomposition in θ retaining all angular modes |m|≤M for M=10 and M=20, for n=2 at λ=1.2 and n=3 at λ=0.8. Compare every eigenvalue below the continuum threshold min(1,λ) with the union of the spectra of H0 (39) and all Hk (41) from the paper. If any extra eigenvalue persists as M is increased, the angular ansatz is not exhaustive; if the full discrete spectrum coincides, the completeness concern is settled.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central Theorem claims a complete classification, but the Lemma in Sec. 4 only shows that the two-angle ansatz (31) with parameter choice (32) is compatible with H+; it does not show that every eigenfunction has angular characteristic number ≤2 or the specific coefficient structure g1 sin(aθ)+g2 sin(bθ), etc. H+ is a coupled second-order system on R2 with background winding n, so in the natural angular-momentum basis the coupling terms can connect harmonics differing by 1 and n; nothing in the Lemma rules out eigenfunctions with more than two (or infinitely many) angular harmonics. The Proposition and Theorem then analyze radial regularity only within the truncated sector (33). Consequently, the phrase 'the spectrum comprises' is not established: the reduction to H0 and Hk is necessary, but sufficiency is assumed. The numerical section solves the reduced 3×3 operator (41) directly, so it cannot detect omitted angular sectors; the agreement with known BPS results only validates λ=1, not the general case. A separate minor inconsistency is that the Sec. 4 summary says Type A modes 'are always associated with a positive eigenvalue', while the numerical tables show Type A k>1 modes with negative ω2 for λ>1; this does not affect the main argument but should be corrected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the linear fluctuation spectrum around rotationally invariant n-vortex solutions of the Abelian-Higgs model for general coupling λ, extending earlier BPS-only analyses. The authors derive a reduced spectral problem: after imposing an angular ansatz and the background gauge, they obtain a 2×2 radial operator H0 for Derrick-type modes and a 3×3 radial operator Hk for multipolar modes, and they classify the latter into Type A and Type B according to the radial behavior at the origin. Numerical eigenvalues and eigenfunction profiles are presented for n=1,...,5 over a range of λ, including negative squared frequencies for λ>1 that are identified with vortex-splitting instabilities. The paper claims to provide the complete internal mode structure of non-BPS Abelian-Higgs vortices.","tokens_in":26687,"tokens_out":2281,"duration_ms":26531,"significance":"If the completeness claim is established, this would be a valuable and long-sought result: a full classification of the internal excitations of higher-charge Abelian-Higgs vortices away from the BPS point, with direct implications for vortex scattering, decay channels, and cosmic-string phenomenology. The derivation of the ODE reductions is a genuine technical asset, the reduction from a 4×4 PDE system on R2 to one-dimensional matrix operators is explicit, and the reproduction of the BPS zero/shape-mode structure at λ=1 in the tables is a useful consistency check. The numerical tabulations provide concrete, falsifiable spectral predictions, and the identification of the k=1 sector with translational modes is exact. The main significance is conditional on the unproven exhaustiveness of the angular reduction and on the internal consistency of the classification statements.","major_comments":[{"comment":"The Lemma shows that the two-angle ansatz (31) with the parameter choice (32) is compatible with the spectral equation, but it does not prove that every eigenfunction of H+ has angular characteristic number at most two. As the proof itself exhibits, H+ couples harmonics with angular momenta differing by 1 and by n, so in the natural angular-momentum basis the coupling network can in principle connect more than two (or infinitely many) harmonics. The Theorem then asserts that the spectrum 'comprises' exactly the Derrick-type, Type A, and Type B modes, but this completeness rests entirely on the unproven exhaustiveness of the ansatz. The numerical section solves only the reduced operators H0 and Hk, so it cannot detect omitted angular sectors. This is a load-bearing gap in the central classification claim.","section":"Sec. 4, Lemma and Theorem"},{"comment":"The text states that Type A multipolar modes 'are always associated with a positive eigenvalue', immediately followed by the claim that the rest of the Type A modes are unstable for λ>1, zero for λ=1, and stable for λ<1. The numerical tables contradict the first statement: for n=2, k=2 at λ=1.2 and 1.4, the entries are ω2=−0.021487 and −0.043283, and similar negative Type A entries appear for n=3,4,5. Since the classification's physical interpretation (n−1 unstable modes for λ>1) depends on Type A modes carrying the negative eigenvalues, this internal inconsistency must be corrected and the wording reconciled with the tables.","section":"Sec. 4, final summary; Figs. 2-6"},{"comment":"The numerical eigenvalues are presented without error estimates, convergence checks, or a systematic study of the dependence on the discretization parameters N and rmax. The appendix fixes N=2000 and rmax=20 and describes the finite-difference scheme, but it does not report how many significant digits of the tabulated eigenvalues are stable under mesh refinement, nor how the results change with rmax. Given that the paper's instability statement ('n−1 negative eigenvalues for λ>1') is a quantitative spectral claim, the numerical section needs at least a convergence table or an estimated error for representative eigenvalues, especially near thresholds where eigenvalues approach zero or the continuum.","section":"Sec. 5 and Appendix A"},{"comment":"The manuscript invokes quasi-bound modes in the ranges λ<ω2<1 and 1<ω2<λ, and the conclusion discusses their existence and physical relevance, but the Theorem's classification and the phrase 'the spectrum comprises' do not explain how quasi-bound modes fit into the completeness claim. Since quasi-bound modes are not square-integrable eigenfunctions, either the Theorem must be restricted to L2 discrete modes and the quasi-bound sector explicitly excluded, or the classification must be extended to include them. As written, the scope of the central claim is ambiguous.","section":"Sec. 2, Sec. 5 final paragraph, Sec. 6"}],"minor_comments":[{"comment":"There is a typo in 'the condiguration is split' (should be 'configuration'), and similar minor typographical issues appear elsewhere (e.g., 'dissappear' in Sec. 6 and 'the review' in the Sec. 3 opening sentence of the Introduction's paper structure paragraph).","section":"Sec. 5.1"},{"comment":"The angular characteristic list for H+ξν is written with semicolons separating rows but is not explained in terms of the definition of chθ; explicitly stating which row corresponds to vector versus scalar components would improve readability.","section":"Sec. 4, proof of Lemma"},{"comment":"The figures plot ω2 versus λ, but it is not stated whether the curves include the quasi-bound regime or only discrete modes; adding a marker or shading for the continuum thresholds would help the reader interpret the graphs.","section":"Sec. 5, Figs. 2-6"},{"comment":"The Data Availability statement says 'no datasets were generated or analysed', yet the paper contains extensive numerical tables and figures; a statement clarifying that numerical data are reproducible from the described method, or providing code or data files, would be more accurate and useful.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a well-motivated problem and the reduction to H0/Hk is a useful contribution, but the central completeness claim currently rests on an unproven angular-sector truncation. This is fixable in principle, but it requires either a genuine proof of exhaustiveness or a careful reformulation of the claim as a classification of the modes within the constructed ansatz. The internal inconsistency about Type A eigenvalues should also be resolved before publication. I do not see evidence of circularity: the BPS results are used as consistency checks, and no fitted parameters enter the derivation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis paper does something real: it tackles the full fluctuation spectrum for rotationally invariant non-BPS Abelian-Higgs vortices, where the supersymmetric BPS toolkit does not apply, and it reduces the problem to coupled radial ODEs (a 2x2 operator for Derrick modes and a 3x3 operator for multipolar modes). The idea of classifying modes by angular characteristics and radial behavior near the origin is clean, and the Type A / Type B split with the explicit k-ranges (k=1,...,n for Type A, unrestricted for Type B) is a substantial unifying step. I checked the algebra where I could: the reduction from the four-field system to the 3x3 operator, the treatment of the background gauge condition, and the consistency argument using the vortex equations look correct. The identification of the k=1 multipolar modes as translational zero modes for every lambda is exact and well-verified.\n\nWhat is genuinely new is the extension beyond the BPS limit, including the prediction of n-1 negative squared frequencies for lambda>1, with the n=1 case correctly showing no instability. The paper also correctly recovers the BPS results at lambda=1, which serves as a useful consistency check. The numerical tables are plausible and match known limits, but they are not evidence of completeness, because they only solve the already-reduced operators. The abstract's claim of \"full internal mode structure\" is therefore a bit strong.\n\nThe soft spots are real but not fatal. The Lemma shows that the ansatz (31)-(32) is compatible with H+, not that every eigenfunction has this angular form. If modes with more angular harmonics exist, the classification misses them. I don't see an obvious way to rule that out from the paper alone. The numerical section gives no error bars or convergence study (mesh size 2000, rmax=20), and the data-availability statement says no datasets were generated, which is odd for a paper full of tables. The minor inconsistency in Section 4—Type A said to be always positive, while the tables show negative eigenvalues for lambda>1—is clearly a wording slip, since the text around it correctly describes instabilities.\n\nVerdict: this deserves a serious referee. The core reduction is a strong, citable result even if the completeness proof needs a sharper argument or a careful statement of assumptions. The paper is honest about prior partial results and does not overstate what is derived versus what is assumed. I would send it to a good referee with instructions to focus on the angular completeness question and to ask for convergence checks, at least for a few representative eigenvalues.\n\nBest,\n[Your name]","headline":"A serious first-principles attack on the non-BPS vortex fluctuation spectrum, with a genuinely new reduction to 2x2 and 3x3 radial operators, but the claimed completeness is not fully proven and the numerics lack error control.","tokens_in":27179,"tokens_out":663,"would_cite":true,"duration_ms":8731,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["11.15.Kc","11.27.+d","11.10.Gh"],"model":"deepseek-v4-flash","headline":"This paper establishes the complete internal-mode spectrum of rotationally invariant n-vortices in the Abelian-Higgs model for arbitrary coupling, reducing every fluctuation to two small radial operators.","keywords":["Abelian-Higgs model","vortex solitons","internal modes","fluctuation spectrum","Derrick modes","multipolar modes","vortex instability","self-dual limit"],"falsifier":"Diagonalize the full four-component fluctuation operator on a two-dimensional domain for a fixed $n$ and $\\lambda$, using a spectral basis that includes angular harmonics up to a high cutoff, and compare every discrete eigenvalue below the continuum thresholds with the union of spectra of $H_0$ and $H_k$; any discrete eigenvalue (or any eigenfunction with angular characteristic number larger than two) not reproduced would disprove the claimed completeness.","tokens_in":26230,"feed_emoji":"🌀","tokens_out":9029,"duration_ms":84527,"temperature":0.7,"pith_summary":"This paper asks what happens when a higher-charge Abelian-Higgs vortex is nudged: what are all the ways it can vibrate, and which of those vibrations make it unstable? Away from the special self-dual coupling $\\lambda=1$, it claims a complete answer for rotationally invariant $n$-vortices. Every normal mode is either a Derrick-type mode that keeps the vortex angular symmetry or one of two classes of multipolar mode, and all of them are governed by a $2\\times2$ radial operator $H_0$ and a $3\\times3$ radial operator $H_k$. If correct, this settles a long-standing problem and explains why type II $n$-vortices split: for $\\lambda>1$ the spectrum contains exactly $n-1$ negative squared frequencies, one for each way the vortex can break apart. The same classification also reproduces the self-dual BPS spectrum at $\\lambda=1$ as a special case.","feed_headline":"Abelian-Higgs vortex modes fully classified","feed_subtitle":"Every small vibration of a charge-n vortex reduces to two radial equations, pinning down type II instabilities.","key_machinery":"The load-bearing object is a two-harmonic angular ansatz for the four-component fluctuation. With the parameter choice $b=-2-a$, $c=n-1-a$, $d=1+n+a$, the full operator $H_+$ acting on this ansatz produces only two independent angular functions in each component, so the spectral problem closes on four radial ODEs. Imposing the background gauge reduces those to a $2\\times2$ operator $H_0$ for the symmetry-preserving $k=0$ modes and a $3\\times3$ operator $H_k$ for the $k\\ge1$ multipolar modes; a consistency check shows the gauge condition is compatible. Finally, power-series analysis at $r=0$ fixes the allowed radial behaviors and separates multipolar modes into Type A and Type B, which is what turns the reduction into a complete spectral classification.","core_discovery":"The central claim is the theorem in Section 4: for a rotationally invariant non-self-dual $n$-vortex, in the background gauge, the fluctuation spectrum of $H_+$ consists exactly of Derrick-type modes with radial characteristic $(1,n)$, Type A multipolar modes with radial characteristic $(k-1,n-k)$ for $k=1,\\dots,n$, and Type B multipolar modes with radial characteristic $(k-1,n+k)$ for $k\\ge1$. The Derrick mode is nondegenerate and has positive squared frequency; every multipolar eigenvalue is doubly degenerate. Type A modes with $k=1$ are the translational zero modes for every $\\lambda$, Type A modes with $k>1$ have positive squared frequency for $\\lambda<1$, vanish at $\\lambda=1$, and become negative for $\\lambda>1$; Type B modes are always positive. The theorem is backed by numerically tabulated eigenvalues for $n=1,\\dots,5$ and by the explicit recovery of the BPS zero modes and shape modes when $\\lambda=1$.","pith_inferences":["The angular-reduction argument uses only rotational invariance and the gauge-covariant structure, so the same two-harmonic ansatz is likely to classify fluctuations around rotationally symmetric vortices in closely related gauge theories; this is our inference, not a claim of the paper.","The linearized picture suggests a concrete dynamical test: excite a Type A mode with $k=m$ in a time-dependent simulation of an $n$-vortex and check whether the vortex splits into $m$ unit vortices arranged in the predicted symmetric pattern while an $(n-m)$-vortex remains at the center.","One sharp numerical check of completeness would be to search the full two-dimensional spectral problem for discrete eigenvalues whose eigenfunctions have three or more angular harmonics; the paper's theorem predicts none exist below the continuum, and finding one would falsify the classification."],"forward_implications":["For $n>1$ and $\\lambda>1$, every rotationally invariant $n$-vortex has exactly $n-1$ unstable directions, matching the ways it can split into smaller vortices under Type A excitations.","The two $k=1$ Type A modes are zero modes for every $\\lambda$, so rigid translations of the vortex persist as exact zero modes away from the BPS limit.","At $\\lambda=1$ the classification degenerates cleanly: Type A modes become the $2n$ BPS zero modes and Type B modes become the positive shape modes, unifying the previously separate self-dual analyses.","The full discrete spectrum can be obtained from the radial operators $H_0$ and $H_k$, reducing a huge two-dimensional eigenproblem to a small one-dimensional one and making high-vorticity spectra numerically accessible.","For large $\\lambda$, the scalar part of the fluctuation stays bound and approaches the modes of an ungauged global vortex, while the vector part becomes scattering, so the vortex develops quasibound scalar modes."],"supporting_citations":[{"why":"Supplies the BPS-limit spectral results (zero modes and shape modes) that the new classification must reproduce and whose radial functions it recovers for $\\lambda=1$.","marker":"[21]"},{"why":"Extends the self-dual spectral analysis and provides the zero-mode ODE and shape-mode solutions used as consistency checks.","marker":"[22]"},{"why":"Earlier partial computation of vortex bound states and instabilities that the paper's complete classification extends.","marker":"[24]"},{"why":"Standard reference for the rotationally invariant $n$-vortex ansatz, the radial ODEs, regularity expansions, and the explicit translational modes.","marker":"[5]"},{"why":"Establishes existence and asymptotic properties of the Abelian-Higgs vortex solutions whose fluctuation spectrum is analyzed.","marker":"[6]"},{"why":"Index-theorem proof of the $2|n|$ zero modes at self-dual coupling, which the Type A modes are claimed to reduce to when $\\lambda=1$.","marker":"[30]"},{"why":"Gives the gauge-invariant translational zero-mode expressions used to identify the $k=1$ Type A modes for arbitrary $\\lambda$.","marker":"[31]"},{"why":"Source for the claim that at large $\\lambda$ the vortex core behaves like a global vortex, supporting the quasibound-mode discussion.","marker":"[32]"}],"fun_headline_variants":["Abelian-Higgs n-vortex modes: three families exactly","Complete spectral classification for non-self-dual vortices","Derrick and multipolar modes: full spectrum announced","All charge-n vortex vibrations now categorized","Two radial equations give every vortex fluctuation mode"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classification assumes that every vibration's shape around the vortex is one of the specific angle patterns written down in ansatz (31); the paper proves these patterns solve the equations, but not that every solution must have that shape.","fun_headline_variants_meta":{"raw":{"variants":["Abelian-Higgs n-vortex modes: three families exactly","Complete spectral classification for non-self-dual vortices","Derrick and multipolar modes: full spectrum announced","All charge-n vortex vibrations now categorized","Two radial equations give every vortex fluctuation mode"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001081,"raw_usage":{"total_tokens":4456,"prompt_tokens":814,"completion_tokens":3642,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":430,"completion_tokens_details":{"reasoning_tokens":3568}},"tokens_in":430,"tokens_out":3642,"duration_ms":29018,"temperature":1.0,"reasoning_tokens":3568,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:14:22.003288+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Diagonalize the full four-component fluctuation operator on a two-dimensional domain for a fixed $n$ and $\\lambda$, using a spectral basis that includes angular harmonics up to a high cutoff, and compare every discrete eigenvalue below the continuum thresholds with the union of spectra of $H_0$ and $H_k$; any discrete eigenvalue (or any eigenfunction with angular characteristic number larger than two) not reproduced would disprove the claimed completeness.","supporting_citations":[{"cited_title":"Alonso-Izquierdo, W","cited_arxiv_id":null,"evidence_quote":"Supplies the BPS-limit spectral results (zero modes and shape modes) that the new classification must reproduce and whose radial functions it recovers for $\\lambda=1$."},{"cited_title":"Alonso-Izquierdo, W","cited_arxiv_id":null,"evidence_quote":"Extends the self-dual spectral analysis and provides the zero-mode ODE and shape-mode solutions used as consistency checks."},{"cited_title":"Goodband and M","cited_arxiv_id":null,"evidence_quote":"Earlier partial computation of vortex bound states and instabilities that the paper's complete classification extends."},{"cited_title":"Manton and P","cited_arxiv_id":null,"evidence_quote":"Standard reference for the rotationally invariant $n$-vortex ansatz, the radial ODEs, regularity expansions, and the explicit translational modes."},{"cited_title":"Jaffe and C.H","cited_arxiv_id":null,"evidence_quote":"Establishes existence and asymptotic properties of the Abelian-Higgs vortex solutions whose fluctuation spectrum is analyzed."},{"cited_title":"Weinberg, Multivortex solutions of the Ginzburg-Landau equations , Phys","cited_arxiv_id":null,"evidence_quote":"Index-theorem proof of the $2|n|$ zero modes at self-dual coupling, which the Type A modes are claimed to reduce to when $\\lambda=1$."},{"cited_title":"Tong and K","cited_arxiv_id":null,"evidence_quote":"Gives the gauge-invariant translational zero-mode expressions used to identify the $k=1$ Type A modes for arbitrary $\\lambda$."},{"cited_title":"Alonso-Izquierdo, J.J","cited_arxiv_id":null,"evidence_quote":"Source for the claim that at large $\\lambda$ the vortex core behaves like a global vortex, supporting the quasibound-mode discussion."}],"review_version":1}