{"id":"59bd9f1b-8477-4354-ae89-3c9e6d5d4b4e","arxiv_id":"2511.11263","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A nominally single-component superconductor can exhibit type-1.5 vortex clustering because a suppressed competing pairing channel generates a second coherence length.","lead":"This paper shows that a superconductor with only one type of pairing active can still have two different coherence lengths — the length scales dictating how superconductivity recovers around defects — because a second, hidden pairing tendency that never appears in the ground state still affects fluctuations. The result predicts vortex clustering, a signature of type-1.5 superconductivity, in materials where only one order parameter is visible.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (C4) Hessian omits the subdominant-channel mass term; as written it predicts an infinite second coherence length, so the analytic two-length hierarchy is not self-consistent.","rationale":"The reader's weakest assumption identifies exactly the same technical flaw: Eq. (C4) as written gives a zero second-diagonal Hessian entry for a pure s-wave or d-wave ground state, leading to an infinite second coherence length. My reading of the manuscript confirms that no alpha-dependent curvature terms are included in the Hessian, and that the text attributes the second length to the mixed-gradient terms in Γ. A direct calculation shows that Γ off-diagonal terms alone cannot produce a second finite eigenvalue if H is diagonal with a zero entry; the mass of the suppressed component is required. This is not merely a cosmetic issue: the analytic estimates are used to draw the type-1.5 region in Fig. 1 and to motivate the central claim that a suppressed competing order parameter changes magnetic response. However, the numerical vortex-clustering results in Figs. 3 and 8 are obtained from the full GL functional, which does include α2 and β3 terms, so those results remain credible. Therefore the appropriate verdict remains CONDITIONAL, matching the reader's assessment: the claimed mechanism needs a corrected analytic derivation, but the numerical evidence for vortex clustering is not directly undermined.","tokens_in":20970,"tokens_out":3657,"duration_ms":32371,"concrete_test":"Recompute the coherence lengths in Fig. 2 for the pure s-wave parameter set with a corrected Hessian: replace H22 by 2(α2 + β3|Δs^GS|^2) (and analogously for pure d-wave), keep H12=0, and diagonalize Γ^{-1}H at the same temperatures. Check whether both eigenvalues are positive and whether λ lies between them over the temperature range shown. If the corrected hierarchy disappears, the estimated type-1.5 region and the claimed gradient-term origin of the second length would need revision; if it persists, the analytic mechanism is confirmed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing weak point is the analytic coherence-length calculation for pure s-wave/d-wave ground states. In Eq. (C4), the Hessian is written with diagonal entries 8β1|Δs^GS|^2 and 8β2|Δd^GS|^2 and off-diagonal entries proportional to |Δs^GS||Δd^GS|. For a pure s-wave state, Δd^GS=0, so H22=0 and H12=0. Then Γ^{-1}H has rank one, yielding a zero eigenvalue and hence an infinite coherence length from Eq. (6). This contradicts Fig. 2, which shows both ξ1 and ξ2 finite in the pure s-wave and pure d-wave regions. The finite second length must come from curvature terms in the potential of the uncondensed component, specifically ∂²Fp/∂|Δd|² = 2α2 + 2β3|Δs^GS|² (plus quartic terms that vanish at zero amplitude), none of which appear in Eq. (C4). The paper's statement that the off-diagonal terms in Γ are 'crucial' for producing two coherence lengths is also not correct when H has zero off-diagonal entries; those terms only hybridize modes once the mass term makes Γ^{-1}H full rank. The numerical vortex-cluster calculations solve the full GL equations and are not invalidated, but the analytic derivation of the two-length hierarchy in Fig. 2 and the type-1.5 boundary in Fig. 1 is not self-consistent as presented.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper considers a single-band square-lattice superconductor with nearest-neighbor pairing, whose link order parameter decomposes into extended s-wave and d-wave irreducible representations. The authors derive a Ginzburg–Landau functional from the microscopic model, compute the mean-field phase diagram (s-wave, d-wave, and s+id regions), estimate coherence lengths and magnetic penetration depth, and solve the resulting GL equations numerically. The central claim is that even when only one order-parameter component is nucleated in the ground state (pure s-wave or pure d-wave), the suppressed competing component generates a second finite correlation length, producing type-1.5 superconductivity with vortex clustering. The paper reports negative vortex binding energies for pure s-wave, pure d-wave, and s+id parameter sets, and additionally reports skyrmionic textures in the s+id case.","tokens_in":21313,"tokens_out":4668,"duration_ms":44433,"significance":"If the mechanism holds, this would substantially broaden the class of systems exhibiting multiple coherence lengths and type-1.5 behavior, showing that a nominally single-component U(1) superconductor can display vortex attraction and clustering. The numerical vortex-clustering results, obtained from microscopically derived GL parameters, are a concrete falsifiable prediction and are the strongest part of the paper. The microscopic derivation of the GL coefficients and the systematic phase diagram are also valuable. However, the analytic coherence-length derivation—which underpins the predicted type-1.5 region in Fig. 1—contains a gap that must be repaired before the paper's central claim is fully supported.","major_comments":[{"comment":"For a pure s-wave ground state, Δd^GS = 0, so the Hessian in Eq. (C4) becomes diag(8β1|Δs^GS|^2, 0) with zero off-diagonal entries. The matrix Γ^{-1}H then has rank one, giving a zero eigenvalue and hence an infinite coherence length from Eq. (C5). This contradicts Fig. 2, which shows both ξ1 and ξ2 finite in the pure s-wave and pure d-wave regions. The finite second length must come from curvature terms in the potential of the uncondensed component—specifically, from terms such as ∂²F_p/∂|Δd|² = 2α2 + 2β3|Δs^GS|² (plus β4-dependent contributions) that are not included in Eq. (C4). The authors should include these terms in the Hessian and recompute the coherence lengths, or explicitly explain why they are negligible. As written, the analytic two-length hierarchy in Fig. 2 is not self-consistent.","section":"Appendix C.1, Eq. (C4) and Eq. (C5)"},{"comment":"The statement that the off-diagonal terms in Γ are 'crucial' for obtaining two coherence lengths in a pure s-wave or d-wave ground state is not correct as it stands. When one order-parameter amplitude vanishes, the off-diagonal entries of H in Eq. (C4) also vanish. Off-diagonal gradient terms cannot create a second nonzero eigenvalue from a rank-one matrix Γ^{-1}H; they only rotate the single nonzero mode. Hybridization through Γ is relevant only when both amplitudes are nonzero. This explanatory claim should be rewritten once the Hessian is corrected.","section":"Sec. II/C.1, text near Eqs. (6)-(7)"},{"comment":"The estimated type-1.5 region in Fig. 1 and the displayed length-scale hierarchy in Fig. 2 rest on the approximate coherence lengths derived from Eq. (C4). Since that derivation yields an infinite coherence length for the pure states, the shaded type-1.5 region is not justified by the equations presented. The numerical vortex clustering in Figs. 3, 6–8 is independent and remains evidence, but the analytic estimate used to select the parameter points needs to be recomputed with the full Hessian of the potential including the subdominant-component mass term.","section":"Fig. 1 and Fig. 2"}],"minor_comments":[{"comment":"The notation γαβ_ij in Eq. (4) is not defined; clarify how it relates to the matrix Γ(φ) in Eq. (7), and specify the index conventions for the anisotropic gradient term.","section":"Eq. (4)"},{"comment":"E1 is used in the normalized binding energy but never defined in the main text. State explicitly that it is the energy of a single isolated vortex and describe how that energy was computed in the simulation.","section":"Fig. 8"},{"comment":"The claim that the two-vortex state has skyrmion charge Q=2 'up to numerical accuracy' should be accompanied by the computed numerical value and an estimate of the discretization error.","section":"Sec. C.3, Eq. (C26)"},{"comment":"The paper would benefit from a data/code availability statement, since the GL coefficients in Table I are provided but the numerical methods are not fully reproducible without additional details on grid sizes, boundary relaxation, and convergence criteria.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The numerical vortex-clustering results are strong and indicate that the central phenomenon may be real. The analytic inconsistency in Eq. (C4) is nevertheless load-bearing for the claimed mechanism and for the predicted type-1.5 region; it is fixable within the manuscript's scope by including the omitted subdominant-channel mass curvature. If, after correction, the second finite coherence length disappears in the pure states, the paper's main claim should be reformulated accordingly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper’s core idea is new and worth taking seriously: a suppressed subdominant pairing channel in a different irreducible representation can imprint an additional correlation length on a nominally single-component U(1) superconductor, and this can produce type-1.5 vortex clustering. The numerical evidence—vortex clusters with negative binding energy for pure s-wave, pure d-wave, and s+id ground states—is direct and looks solid. The GL coefficients are derived from a microscopic model, not fitted, which adds weight. The vortex morphologies, including the diagonal chains and the skyrmion charge in the s+id case, are nice extras. This is a real contribution to the type-1.5 literature.\n\nThe soft spot is real and located exactly where the stress-test puts it. In Eq. (C4), the Hessian for a pure s-wave ground state has H22 = 0 because Δd^GS = 0. Taken literally, that gives a zero eigenvalue and an infinite second coherence length from Eq. (6), contradicting Fig. 2. The finite second length must come from the curvature of the potential in the subdominant direction—the 2α2 + 2β3|Δs^GS|^2 term—which is absent from Eq. (C4). The statement that the off-diagonal gradient terms are “crucial” is also wrong: those terms cannot generate a finite eigenvalue when H has a zero row. This makes the analytic length-scale hierarchy and the type-1.5 boundary in Fig. 1 not self-consistent as presented. It is a fixable error, but the paper’s explanation of the mechanism is incorrect as written.\n\nThat said, the numerical vortex-cluster calculations solve the full GL equations and do not depend on the flawed approximation. They are the real evidence for the effect, and they appear convincing. So the qualitative conclusion may survive a correction, but the analytic derivation needs to be redone and the narrative about gradient terms reframed.\n\nMinor issues: no code or data are provided, and the numerical method is somewhat underspecified, but the physics is clear enough for a specialist to reproduce.\n\nThis paper deserves a serious referee. With a corrected Hessian and a revised explanation of the origin of the second length scale, it would be a solid contribution. I would bring it to a reading group and would cite it, with a caveat about the current form of the analytic derivation.","headline":"A genuinely new idea—multiple correlation lengths and type-1.5 vortex clustering in a nominally single-component U(1) superconductor—backed by convincing numerics, but the analytic coherence-length derivation has a missing Hessian term that needs correction before the stated mechanism is credible.","tokens_in":21836,"tokens_out":3746,"would_cite":true,"duration_ms":35129,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82D55"],"pacs":["74.20.De"],"model":"deepseek-v4-flash","headline":"A single-component U(1) superconductor can hide a second correlation length, and that changes its vortex physics.","keywords":["multiple correlation lengths","type-1.5 superconductivity","vortex clustering","s-wave/d-wave pairing","Ginzburg-Landau theory","irreducible representations","nonlocal pairing","time-reversal symmetry breaking"],"falsifier":"A calculation of the exact linearized fluctuation spectrum of the microscopic self-consistency equation at a point where the ground state is pure s-wave or pure d-wave would settle the question: if only one finite coherence length survives when the subdominant amplitude is exactly zero, the central claim collapses. Experimentally, measuring the vortex-vortex interaction energy as a function of separation in a material with a near-degenerate second pairing channel would distinguish type-1.5 clustering from a conventional type-II lattice.","tokens_in":20798,"feed_emoji":"🌀","tokens_out":5083,"duration_ms":44344,"temperature":0.7,"pith_summary":"The paper argues that even a superconductor whose ground state breaks only U(1) symmetry and contains a single order-parameter component (pure s-wave or pure d-wave) can be effectively multicomponent in its response to perturbations. A subdominant pairing channel belonging to a different irreducible representation of nonlocal pairing does not nucleate in the ground state, but its virtual presence in the gradient-energy structure imprints an additional coherence length. When this second length exceeds the magnetic field penetration depth, the system enters the type-1.5 regime: vortices repel at short range but attract at long range, producing vortex clusters instead of a vortex lattice. The claim matters because it challenges the usual identification of one order parameter with one coherence length and enlarges the class of materials that can show type-1.5 behavior.","feed_headline":"One condensate, two lengths: hidden channel makes type 1.5","feed_subtitle":"Even when only s-wave or d-wave order is nucleated, a suppressed partner imprints a second coherence length and clusters vortices.","key_machinery":"The central object is the two-component Ginzburg-Landau free energy derived from the microscopic lattice model, with order parameters Delta_s and Delta_d, anisotropic gradient terms gamma1, gamma2, a mixed gradient term gamma12, and interband couplings beta3, beta4. The carrier of the argument is the matrix Gamma(phi) describing the stiffness of amplitude fluctuations as a function of direction; its off-diagonal entries, proportional to gamma12 cos 2phi, hybridize the two channels. Through the eigenvalues of Gamma^{-1} H, one obtains two direction-dependent coherence lengths xi1, xi2 even when the subdominant amplitude is zero in the ground state, because the mixed gradient term couples fluc","core_discovery":"The central discovery is that multiple correlation lengths do not require multiple broken symmetries or multiple bands. In a model with nearest-neighbor pairing on a square lattice, the link order parameter decomposes into extended s-wave (A1g) and d-wave (B1g) irreducible representations. Even when only one of these is condensed, the mixed gradient coupling gamma12 between the two channels makes amplitude fluctuations hybridize, so the linearized Ginzburg-Landau spectrum has two finite coherence lengths. Near a second transition into an s+id state with time-reversal symmetry breaking, one of these lengths diverges; with the penetration depth lying between the two coherence lengths, vortex i","pith_inferences":["Editorial inference: the mechanism should generalize beyond s/d channels; any superconductor with a nearby instability in a different irreducible representation (for example, a d+ig state) could display a second correlation length even when only one component is condensed.","Editorial inference: a direct microscopic test would be to compute the exact fluctuation spectrum of the self-consistency equation (without the GL Hessian approximation) in the pure s-wave or d-wave state and check whether a finite second length survives at zero subdominant amplitude.","Editorial inference: scanning-tunneling measurements of the local density of states around a vortex could reveal the subdominant order parameter nucleating in the core, providing an experimental fingerprint of the hidden length scale.","Editorial inference: if the analytic estimate is corrected by including the omitted curvature terms, the type-1.5 region on the phase diagram may shift; the numerical vortex-clustering results stand independently of that estimate."],"forward_implications":["Ground states that look single-component in equilibrium can exhibit two correlation lengths; one coherence length diverges at both the U(1) transition and the transition to a time-reversal-symmetry-broken s+id state, even when the second component never condenses.","Within the estimated parameter region, the magnetic penetration depth can fall between the two coherence lengths, so vortices attract each other on the scale of the longer length and form clusters rather than a triangular lattice.","Vortex clusters have anisotropic morphology: in the s+id state they form diagonal chains for two, three and five vortices, whereas four-vortex clusters can form a diamond or a chain depending on the ground state, offering a diagnostic of the pairing symmetry.","Vortex clusters in the s+id state carry skyrmion charge Q=2 as bound pairs of Q=1 skyrmions, with an additional density-density attraction in the type-1.5 regime.","A nominally single-component U(1) superconductor cannot be type-II in proximity to a continuous transition to a different pairing state with additional broken symmetry."],"fun_headline_variants":["Two coherence lengths without two condensates","Hidden order channel creates type-1.5 superconductor","Suppressed pairing still sets two lengths","Single-band superconductor shows multi-length behavior","Vortex clustering from a ghost pairing channel"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The analytic claim that two finite coherence lengths exist for a pure s-wave or d-wave ground state rests on the Hessian in Eq. (C4) being supplemented by curvature terms (alpha2 and beta3|Delta|^2 terms) for the uncondensed component; as written, the Hessian has a zero diagonal entry, which would give an infinite second coherence length.","fun_headline_variants_meta":{"raw":{"variants":["Two coherence lengths without two condensates","Hidden order channel creates type-1.5 superconductor","Suppressed pairing still sets two lengths","Single-band superconductor shows multi-length behavior","Vortex clustering from a ghost pairing channel"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1077,"prompt_tokens":685,"completion_tokens":392,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":429,"completion_tokens_details":{"reasoning_tokens":323}},"tokens_in":429,"tokens_out":392,"duration_ms":4261,"temperature":1.0,"reasoning_tokens":323,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T22:14:05.874470+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A calculation of the exact linearized fluctuation spectrum of the microscopic self-consistency equation at a point where the ground state is pure s-wave or pure d-wave would settle the question: if only one finite coherence length survives when the subdominant amplitude is exactly zero, the central claim collapses. Experimentally, measuring the vortex-vortex interaction energy as a function of separation in a material with a near-degenerate second pairing channel would distinguish type-1.5 clustering from a conventional type-II lattice.","supporting_citations":[],"review_version":1}