{"id":"c2a5ce16-7410-4bb7-91bb-2a6bf85098c5","arxiv_id":"2510.04830","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Vortices in the Chern–Simons–Ginzburg–Landau model can form stable multi-vortex bound states via short-range repulsion and long-range attraction, realizing hybrid type-I/II superconductivity in a single component.","lead":"This paper studies superconducting vortices in a 2D model with a Chern–Simons term and finds their interactions can be repulsive at short range yet attractive at long range. This suggests a single-component superconductor can realize 'hybrid' vortex bound states, not just the usual type-I or type-II behavior.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central bound-state claim rests on uncomputed asymptotic coefficients c_H, c_B, γ; Eq. (53) cannot be evaluated and is used beyond its large-R validity to infer short-range repulsion.","rationale":"The reader's verdict of CONDITIONAL identifies exactly the same load-bearing weakness: the interaction-energy formula (Eq. 53) is written with coefficients that are never computed, and the short-range repulsion/long-range attraction conclusion depends on these coefficients and on the gauge term dominating at short range. My stress-test concurs. This is not an internal inconsistency in the derivation of the complex screening masses (Eq. 25) or the oscillatory gauge tails, which appear coherent. It is a gap between the analysis and the central physical claim. The numerical figures referenced in the text are not reproducible from the manuscript, and no error bars or convergence tests are given. However, the mechanism is plausible and the missing pieces could plausibly be supplied by additional computation. Therefore, the verdict should remain CONDITIONAL: the paper should not be accepted as fully establishing hybrid superconductivity without either analytical derivations of c_H, c_B, γ or release of the numerical data/code. I do not see grounds to reject outright, since the core mechanism may be correct and the asymptotic screening result is a meaningful contribution. The minor sign typo in Eq. (32) does not affect the final operator (33) when corrected, so it is not the primary concern.","tokens_in":13642,"tokens_out":7670,"duration_ms":62240,"concrete_test":"For a representative parameter set (e.g., m=1, q=1, λ=0.5, κ=0.5), compute the single-vortex profile numerically using the described constrained Newton flow, extract c_H from ϕ(r) ≈ c_H K0(m_H r) and c_B, γ from B(r) ≈ 2Re[c_B K0(m_+ r)] in the tail region, and evaluate Vint(R) from Eq. (53). Then directly compute the two-vortex interaction energy by solving the full static equations for two vortices at separation R (e.g., using the same flow algorithm) and compare the location and depth of the predicted minimum with the direct result. If Eq. (53) does not reproduce the direct Vint(R) at the minimum, the bound-state claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that CSLG vortices realize hybrid superconductivity with short-range repulsion and long-range attraction, forming stable multi-vortex bound states—depends on the interaction energy Vint(R) in Eq. (53). This formula contains three undetermined parameters: the single-vortex tail amplitudes c_H and |c_B|, and the phase γ (equivalently ϑ = arg c_B). These are not computed from the nonlinear vortex profiles nor fixed by the linearized analysis. Consequently, the sign and magnitude of the gauge channel at any separation R cannot be evaluated analytically. Furthermore, Eq. (51) is a large-R asymptotic expansion of K0(m+R); the hypothesized short-range repulsion would occur at R comparable to the core size, where this expansion is not valid. The text itself acknowledges the conditionality: 'If the gauge term dominates at a shorter range than the Higgs term, it can provide a repulsive force initially...'. The only evidence for actual bound states is Fig. 1 and Fig. 2, but the numerical implementation is described without releasing code or data, and no convergence or error analysis is provided. Thus, the existence of stable multi-vortex bound states is asserted rather than demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Chern–Simons extension of the single-component Ginzburg–Landau model (CSLG) and claims that its vortices realize 'hybrid superconductivity': short-range repulsion and long-range attraction, leading to stable multi-vortex anyon bound states. The paper derives a positive-definite static energy using Gauss' law, obtains the charge–flux relation, introduces a constrained Newton-flow method, and linearizes the field equations to find complex-conjugate gauge screening masses. It then uses a point-particle source method to write an asymptotic interaction energy, Eq. (53), and combines this with numerical binding-energy plots to conclude that the type-I/type-II dichotomy is broken in a single-component condensate.","tokens_in":13943,"tokens_out":6397,"duration_ms":64552,"significance":"If established, this result would be notable: a topological gauge term in a single-component superconducting model would generate non-monotonic intervortex forces and molecular-like vortex bound states, extending the phenomenology of type-1.5 superconductivity. The linearized analysis is largely coherent, the positivity of the static energy is a useful step, and the source-method reconstruction of the asymptotic tails is a principled approach. The κ→0 consistency check is a genuine strength. However, the central bound-state claim is not yet supported: the asymptotic interaction energy contains undetermined amplitudes and phases, the asymptotic formula is used beyond its validity, and the numerical evidence is not reproducible from the information given.","major_comments":[{"comment":"The central prediction of hybrid superconductivity rests on Vint(R) in Eq. (53), yet the constants c_H, |c_B|, and γ are never determined. They are single-vortex far-field amplitudes that must be obtained from the nonlinear vortex profiles, but no calculation or numerical value is given. Consequently, the sign and magnitude of the gauge term at any separation cannot be evaluated; the claimed short-range repulsion is explicitly conditional in the text ('If the gauge term dominates...'). Without these coefficients, or an alternative computation of Vint, the analytic argument does not establish bound states.","section":"Hybrid superconductivity, Eq. (53)"},{"comment":"Eq. (51) is the large-R asymptotic expansion of K0(m+R), with the exponentially small prefactor 1/sqrt(R). Eq. (53) is then used to infer a repulsive region at R < π/κ + 2γ/κ, which is precisely the regime where this asymptotic expansion may not be valid. The short-range behavior should be obtained from the exact expression in Eq. (50) with determined c_B, or from direct numerical evaluation of the interaction energy. As written, the argument uses an asymptotic formula outside its domain of validity.","section":"Appendix 2, Eq. (51)"},{"comment":"The existence of stable multi-vortex anyon bound states is the central claim, but the numerical evidence is not reproducible. The paper does not release code or data, state the lattice size h, the number of grid points, the boundary conditions, or provide convergence tests with respect to grid resolution and flow tolerance. Figure 2 shows a single N=4 configuration, and Fig. 1 shows binding energies without error estimates. The reader cannot verify that the minima are not numerical artifacts.","section":"Multi-vortex anyon bound states; Figs. 1 and 2"},{"comment":"The damped oscillatory form and the common penetration depth require α = sqrt(m_A^2 - κ^2/4) to be real, i.e. κ < 2m_A. The manuscript never states this condition or discusses the regime κ ≥ 2m_A, where the asymptotic tails (38) and interaction energy (53) are not of the claimed form. Either the analysis should be explicitly restricted to κ < 2m_A, or the κ ≥ 2m_A case should be addressed.","section":"Hybrid superconductivity, Eqs. (24)–(26)"}],"minor_comments":[{"comment":"Applying ∇² to Eq. (30c) gives (−∇² + m_A²)∇² a0 = −κ∇² B, not (∇² + m_A²)∇² a0 as printed. The sign error appears to be a typo, since the final operator in Eq. (33) is correct, but the printed relation should be fixed.","section":"Appendix 1, Eq. (32)"},{"comment":"Several references are incomplete: [22] lacks year/article number, [27] lacks volume/page, and [42] lacks volume/page. These should be completed before publication.","section":"References"},{"comment":"The notation is confusing: λ denotes the Ginzburg–Landau parameter in Eq. (23), while the magnetic penetration depth is also denoted λ_H in the same paragraph. Please use distinct symbols for the GL parameter and the penetration depth.","section":"Multi-vortex anyon bound states"},{"comment":"The caption and axis labels should specify the exact parameters used (m, q, and the meaning of the curves labelled 'κ = 0, 0.5, 1'), and ideally include error bars or convergence indicators.","section":"Figure 1"}],"recommendation":"major_revision","confidential_remarks":"The paper's central idea is interesting and the analytic framework is mostly coherent, but the central claim is not yet supported. I recommend major revision rather than rejection because the missing coefficients, the validity-regime discussion, and the numerical documentation could in principle be supplied. The self-citations to the numerical method are not load-bearing. The paper would be acceptable only after the interaction-energy prediction is either derived with determined coefficients or replaced by robust numerical evidence with reproducibility details."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Paul,\n\nThe thing to know: the paper splits into a solid analytic core and a headline claim that is not yet proven. The analytic core—the linearized screening-mass calculation for the CSLG model—is worth taking seriously. The static masses are m± = sqrt(m_A^2 - kappa^2/4) ± i kappa/2, and the author correctly distinguishes them from the dynamical pole masses of Pisarski–Rao and corrects the sign in Paul–Khare. The common penetration depth plus oscillatory phase for both B and E is a genuine observation, and the kappa to 0 limit recovers GL.\n\nThe problem is the central bound-state claim. Eq. (53) depends on the single-vortex amplitudes c_H, c_B and phase gamma, which are never computed. Without their signs and magnitudes, the gauge channel in (53) cannot be evaluated at any R, and the advertised short-range repulsion/long-range attraction remains a possibility, not a result. The author does phrase it conditionally at one point, but the abstract and introduction do not. Also, the numerics in Figs. 1–2 are the only evidence for actual bound states, and no code, data, or convergence study is provided. The use of the large-R Bessel expansion to conclude something about short-range repulsion is questionable; the first zero can sit at separations where m_A R is of order one.\n\nThere are minor sign slips in intermediate equations (32) and (40c) that do not appear to propagate, but they add to the sense of a paper that has not been polished for the central claim. The positive-definite energy (7) is clean, and the point-particle source method is recognizably the Speight/Manton framework; the derivation is mostly consistent.\n\nBottom line: send to peer review. A referee can verify the analytic core and push the author for the missing tail coefficients and for the numerics. The paper is for people working on CS vortices and intervortex forces; the mass analysis alone is citable. The hybrid bound-state conclusion should be softened until the coefficients are computed.","headline":"Static-mass analysis is a real contribution; the hybrid bound-state claim is not yet supported by the analytic result or by reproducible numerics.","tokens_in":14404,"tokens_out":11589,"would_cite":true,"duration_ms":85666,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single-component superconductor with a Chern–Simons term can host stable multi-vortex anyon bound states, breaking the type-I/type-II dichotomy.","keywords":["anyons","Chern–Simons theory","Ginzburg–Landau model","superconductivity type","vortex interactions","bound states","flux-charge relation","type-1.5 superconductivity"],"falsifier":"Compute the single-vortex amplitudes c_H, c_B, γ from numerically relaxed vortex profiles, insert them into Eq. (53) for λ < 1, κ > 0, and check whether the pair interaction Vint(R) actually rises at short range and then falls; if Vint(R) is monotone, the hybrid bound-state claim fails. Alternatively, relax two well-separated vortices in the full field equations and measure whether the equilibrium separation is finite and nonzero.","tokens_in":13533,"feed_emoji":"🌀","tokens_out":3930,"duration_ms":31860,"temperature":0.7,"pith_summary":"The paper argues that adding a Chern–Simons term to the standard Ginzburg–Landau model turns each vortex into an anyon carrying both a flux quantum and a proportional electric charge, and rewires how vortices interact. Linearizing about the vacuum shows the gauge fields acquire complex-conjugate screening masses, so magnetic and electric fields share one penetration depth but decay with an oscillatory phase. The resulting pair interaction has a damped oscillatory gauge part competing with the usual Higgs attraction: repulsive at short range, attractive at longer range. If the gauge term dominates at the right scale, vortices settle into separated, molecular-like bound clusters rather than collapsing into one core or forming a type-II lattice. This would realize hybrid type-I/type-II behavior in a single-component condensate, without any second band.","feed_headline":"Chern–Simons vortices bind into stable anyon clusters","feed_subtitle":"A single-component superconductor gains hybrid type-I/type-II behavior from oscillatory gauge tails and short-range repulsion.","key_machinery":"The central object is the Chern–Simons–Landau–Ginzburg static energy with its Gauss-law constraint. It produces (i) the charge-flux relation Q_m = −κΦ that makes each vortex anyonic, and (ii) a fourth-order linearized operator Δκ = (−∇² + m_A²)² + κ²∇² for the gauge fields. Factorizing Δκ gives complex-conjugate screening masses m± = α ± iβ, α = sqrt(m_A² − κ²/4), β = κ/2, so the gauge channel oscillates as it decays; this is the mechanism that converts the usual type-I/type-II force dichotomy into a non-monotonic force law.","core_discovery":"The central claim is that the Chern–Simons extension of the single-component Ginzburg–Landau model breaks the conventional type-I/type-II dichotomy by giving vortices a dual interaction. Through Gauss' law, each flux quantum binds a Noether charge proportional to the Chern–Simons level, so vortices are anyons. The Chern–Simons term also changes the screening spectrum from two real gauge masses to one complex-conjugate pair m± = sqrt(m_A² − κ²/4) ± iκ/2, meaning the magnetic and electric tails have a common decay length but oscillate in sign. The asymptotic interaction energy between two vortex anyons is then Vint(R) = 2π|c_B|²√(2π/(m_A R)) e^{−αR} cos(βR−γ) − 2π c_H² K0(m_H R), a damped osci","pith_inferences":["Because the interaction formula depends on amplitudes c_H and c_B that must come from the single-vortex profile, the paper's bound-state claim is only as strong as those coefficients; a natural next step is to extract them numerically and check the sign of the short-range gauge force directly.","The same mechanism — a Chern–Simons term generating oscillatory gauge tails and non-monotonic intervortex forces — suggests that any parity-violating term producing complex screening masses could create type-1.5-like behavior in other single-component systems; the author draws the analogy to non-centrosymmetric superconductors, but does not test it.","If the bound clusters are stable and anyonic, braiding them may give rise to fractional or non-Abelian statistics at the cluster level, a consequence not explored in the paper.","The damped oscillatory interaction also implies preferred inter-vortex spacings, which could manifest as commensurate vortex cluster phases or novel lattices, a direction the paper leaves open."],"forward_implications":["Even at the Bogomolny point λ = 1, where ordinary Ginzburg–Landau vortices are neutrally stable, a nonzero Chern–Simons level makes vortex anyons repel.","In the nominal type-I regime (λ < 1), an intermediate Chern–Simons coupling frustrates collapse and stabilizes separated multi-vortex clusters with negative binding energy.","For sufficiently large Chern–Simons coupling, repulsion dominates at all distances and the system mimics type-II behavior even when λ < 1.","The resulting molecular-like vortex clusters are neither a giant type-I core nor an Abrikosov lattice, so anyon superconductors can support collective states outside the standard classification."],"fun_headline_variants":["Chern-Simons anyons break superconductivity's type-I/II rule","Vortex anyons attract at long range via oscillatory tails","Hybrid superconductivity from Chern-Simons vortex binding","Chern-Simons anyons: superconductor type-I/II hybrid"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The existence of bound states hinges on the gauge term in the asymptotic interaction energy being repulsive and dominant at short range; the paper never fixes the amplitudes c_B, c_H and phase γ from the vortex profiles, so the sign and range of that term are not established by the calculation itself.","fun_headline_variants_meta":{"raw":{"variants":["Chern-Simons anyons break superconductivity's type-I/II rule","Vortex anyons attract at long range via oscillatory tails","Hybrid superconductivity from Chern-Simons vortex binding","Chern-Simons anyons: superconductor type-I/II hybrid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001167,"raw_usage":{"total_tokens":4650,"prompt_tokens":713,"completion_tokens":3937,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":457,"completion_tokens_details":{"reasoning_tokens":3859}},"tokens_in":457,"tokens_out":3937,"duration_ms":24952,"temperature":1.0,"reasoning_tokens":3859,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T11:23:01.850051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the single-vortex amplitudes c_H, c_B, γ from numerically relaxed vortex profiles, insert them into Eq. (53) for λ < 1, κ > 0, and check whether the pair interaction Vint(R) actually rises at short range and then falls; if Vint(R) is monotone, the hybrid bound-state claim fails. Alternatively, relax two well-separated vortices in the full field equations and measure whether the equilibrium separation is finite and nonzero.","supporting_citations":[],"review_version":1}