{"id":"6b225843-2d21-4fc1-abaf-0df29e716af1","arxiv_id":"1908.10847","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For large gauge coupling, SU(N) Ginzburg-Landau models with N > 2 show two phase transitions, with an intermediate CP^{N-1}-neutral phase that has composite order, no Meissner effect, and no O(M) description.","lead":"Simulations of multicomponent superconducting models with SU(N) symmetry and N > 2 reveal an intermediate phase where composite neutral order survives after the Meissner effect disappears. This CP^{N-1}-neutral phase does not reduce to the simpler O(M) models used for two-component systems, and it classifies new physics for theories relevant to quantum criticality and ultracold atom simulators.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The hard total-density constraint after Eq. (1) is what makes the neutral sector CP^{N-1}; without it the claimed intermediate phase may not survive in a soft-density SU(N) Ginzburg-Landau model.","rationale":"Reader's weakest_assumption is the hard constraint, and I agree. The paper's two transitions for N=3, q=6 are supported by L*Upsilon and L*rho crossings with heat-capacity peaks, and the derivation of helicity-modulus proportionality in Sec. IV.B is internally consistent: for zero-sum twists all moduli are proportional because of SU(N) symmetry and the zero phase-sum stiffness. So I do not see a fatal flaw. The concern is scope: the CP^{N-1} target space arises only under the fixed-density constraint; the radial mode is absent by construction. Since the abstract and title do not qualify 'SU(N) Ginzburg-Landau models' as fixed-density models, the generality of the main claim is unproven. A direct soft-density simulation would settle it. This does not change the reader's conditional verdict: the paper is publishable with the caveat, but the central claim should be framed as applying to the fixed-density model or tested against soft densities.","tokens_in":11057,"tokens_out":12987,"duration_ms":141747,"concrete_test":"Run the N=3, q=6 system with the constraint relaxed to a potential V = lambda * (sum_i |psi_i|^2 - 1)^2, and measure L*Upsilon and L*rho crossings for lambda = infinity, 10, 1, 0.1 at L = 16, 24, 32, 40. If the charged and neutral crossings approach each other or the L*Upsilon crossing disappears as lambda decreases, the intermediate phase is an artifact of the hard constraint; if the crossings remain split for finite lambda, the abstract's generalization is justified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The model is defined by Eq. (1) with the hard constraint sum_i |psi_i|^2 = 1 imposed immediately afterward, and the rewriting in Eq. (2) shows that the neutral sector is the CP^{N-1} model of Eq. (13). This is the load-bearing step for the central claim: the existence of a composite-neutral intermediate phase depends on freezing out the radial total-density mode. A standard Ginzburg-Landau expansion would keep the |psi_i| amplitudes soft and include a potential in |psi_i|^2, so the total-density mode can fluctuate and couple to the charged current j in Eq. (3). The paper never relaxes this constraint, so the abstract's statement about 'SU(N)-symmetric Ginzburg-Landau models' is, strictly, supported only for the fixed-density model. If soft density fluctuations disorder the relative phase/density order or merge the two transitions, the headline claim fails beyond a special limit. This is a limitation rather than an internal contradiction, but it is the single most load-bearing assumption because the whole CP^{N-1}-neutral concept is constructed on it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies SU(N)-symmetric Ginzburg-Landau models coupled to a non-compact Abelian gauge field in three dimensions, with a hard constraint on the total matter density. For N=3 and N=4 it presents Monte Carlo evidence that, at sufficiently large gauge coupling q, the single superconducting transition splits into two: a lower-temperature transition where the Meissner effect disappears and an upper transition where the neutral CP^{N-1} sector disorders. The intermediate phase is termed the CP^{N-1}-neutral phase, and the paper argues it cannot be mapped onto an O(M) model for N>2. The paper also studies vortex configurations in an external magnetic field, finding regular vortex lattices in both the q=0 and finite-q cases despite the absence of a conserved U(1) topological invariant.","tokens_in":11179,"tokens_out":4122,"duration_ms":43995,"significance":"If the central claim holds, the paper identifies a new type of composite-order phase in multicomponent gauge theories with N>2, going beyond the well-studied U(1)xU(1) and SU(2) cases. The numerical evidence for N=3, q=6 is a strength: the authors use finite-size scaling of both the dual stiffness and the helicity modulus, parallel tempering, reweighting, and bootstrap errors, with system sizes up to L=40-48. The vortex-lattice results in external fields are also concrete and falsifiable. However, the significance is moderated by the fact that the central claim is established only for the fixed-total-density model, and the paper does not quantify the extrapolation of crossing temperatures or rule out weakly first-order transitions.","major_comments":[{"comment":"The hard constraint sum_i |psi_i|^2 = 1 is load-bearing for the central claim, because it is what projects the matter fields onto S^{2N-1}/U(1) = CP^{N-1} and hence what makes Eq. (13) the CP^{N-1}-neutral sector. The abstract and conclusion state results for 'SU(N)-symmetric Ginzburg-Landau models' without this qualifier, yet the paper never studies the soft-density version in which the total-density mode fluctuates and couples to the charged current in Eq. (3). Since the existence and the non-O(M) character of the intermediate phase depend on this restriction, the claim as stated is overbroad; the authors should either restrict the abstract and conclusions to fixed-total-density models or provide evidence that the split transition survives when the total-density mode is allowed to fluctuate.","section":"Section II, after Eq. (1)"},{"comment":"The finite-size crossing temperatures for L*Upsilon and L*rho are not tabulated or extrapolated, and no criterion is given for distinguishing continuous from weakly first-order transitions. Since the separation of the two transitions is the quantitative basis for the split-transition claim, the paper should report the crossing values, their extrapolation to L -> infinity, and an order-of-transition analysis (e.g., histogram or Binder cumulant) for at least the N=3, q=6 case.","section":"Section IV.A-IV.C"},{"comment":"The phase diagram for N=4 is presented without any finite-size crossing data, while the detailed evidence is given only for N=3 at q=3 and q=6. Since the abstract emphasizes N>2, the N=4 branch of the central claim needs either corresponding data or an explicit statement that it is an extrapolation from the N=3 case.","section":"Section IV.C, Fig. 1"},{"comment":"The statement that for N>2 the neutral phase 'cannot be mapped onto an O(M) model' is asserted rather than demonstrated. The symmetry argument identifies CP^{N-1} as the target manifold of the neutral sector, but the paper does not rule out an effective description in terms of O(M) variables after the charged sector is disordered; the authors should specify what 'mapped onto' means (target-space topology, critical exponents, or both) and provide the corresponding argument.","section":"Section IV.C, after Eq. (13)"}],"minor_comments":[{"comment":"There are several typos, including 'quantium' (Introduction), 'gauage' (Introduction), 'disaplayed' (figure captions), 'apear' (Section IV.B), and 'reweigting' (Section IV.B).","section":"Throughout"},{"comment":"The phrase 'Form this is follows that' should be 'From this it follows that'.","section":"Section IV.B"},{"comment":"The text uses both 'q -> 0' and 'q = 0' for the zero-charge case; please unify the notation and clarify how the q=0 limit is implemented given that the vector potential still appears in the Hamiltonian.","section":"Section V.B and Figs. 5-10"},{"comment":"The wave vector q in the exponential of the dual stiffness conflicts notationally with the electric charge q used elsewhere; consider renaming one of them.","section":"Eq. (8)"},{"comment":"The statement that 'the phase-sum helicity modulus being the special case where the constant of proportionality is zero' is confusing; clarify that the proportionality constant for the phase-sum combination is zero, rather than implying the modulus itself is zero in a special case.","section":"Section IV.B"},{"comment":"The phase diagram would benefit from a table of the numerical crossing temperatures and their errors, since the lines are described only as a guide to the eye.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The hard-density constraint is the main risk to the central claim. If the authors can show that the split transition persists when the total-density mode is made soft, or alternatively restrict the claims to the fixed-density model, the paper would be publishable. The lack of tabulated crossing values and the absence of an order-of-transition analysis for the charged transition should be addressed in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this one. First, it is a genuinely new numerical result: for SU(N) models with N>2 and large gauge coupling, the Monte Carlo data show a clean split transition into an intermediate phase with composite neutral order and no Meissner effect. That phase does not map onto an O(M) model, which is a real departure from the SU(2) and U(1)^N cases. Second, the whole CP^{N-1} structure is built into the model by the hard constraint sum_i |psi_i|^2 = 1 imposed right after Eq. (1). The paper never tests the soft-density limit, so the abstract's phrase \"SU(N)-symmetric Ginzburg-Landau models\" overstates the generality of what is demonstrated.\n\nThe strengths are real. The N=3, q=6 case shows two well-separated transitions with growing heat-capacity peaks and clean finite-size crossings for L=12 to 40. The phase diagram for N=4 is a useful addition. The vortex simulations in external field, especially the half-quantum-flux hexagonal pattern in the q=0 SU(3) case, are a nice bonus and go beyond the earlier SU(2) work. Error handling looks careful: bootstrapping, parallel tempering, reweighting, and consistency checks between halves of the data.\n\nSoft spots, in proportion. The hard-density constraint is the load-bearing assumption, and the stress-test note is right: the neutral sector is CP^{N-1} by construction, so the intermediate phase may not survive in a more physical soft-density Ginzburg-Landau model. That is a limitation, not an internal contradiction, but it should have been flagged more clearly in the abstract. The transition order is not firmly established; if the charged transition is weakly first order, the L*rho crossings could be biased. The helicity-modulus proportionality is asserted rather than derived, though symmetry makes it plausible. No code or data are released, so independent checks are limited.\n\nWho should read it: people working on multicomponent superconductors, deconfined criticality, or synthetic gauge fields in cold atoms. It deserves a serious referee. I would send it to review and ask for a broader abstract, a comment on the soft-density limit, and perhaps a test of transition order at larger L. None of these issues is fatal. It is a solid, careful study of a well-defined model, and the N>2 phase diagram is worth having.","headline":"Solid Monte Carlo evidence for a new CP^{N-1}-neutral phase in SU(N) gauge models with N>2, under a fixed-density constraint the paper never relaxes.","tokens_in":11875,"tokens_out":1915,"would_cite":true,"duration_ms":21029,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["74.20.De","74.25.-q","05.70.Fh"],"model":"deepseek-v4-flash","headline":"At sufficiently large gauge-field coupling, SU(N) Ginzburg-Landau models with N>2 have two phase transitions, with a CP^{N-1}-neutral intermediate phase that cannot be mapped onto an O(M) model.","keywords":["SU(N) Ginzburg-Landau theory","non-compact Abelian gauge field","CP^{N-1} neutral phase","composite order","vortex lattice","Meissner effect","finite-temperature phase transitions","Monte Carlo simulation"],"falsifier":"Simulate the N=3, q=6 model with a soft total-density potential such as $\\lambda(\\sum_i|\\psi_i|^2-1)^2$ at finite $\\lambda$ and check whether the $L\\Upsilon$ crossing for the neutral transition persists as $\\lambda\\to\\infty$; separately, measure the latent heat or Binder cumulant at the charged transition for q=7, because a clear first-order signature would invalidate the continuous-transition crossings used to locate the split.","tokens_in":10704,"feed_emoji":"🌀","tokens_out":12042,"duration_ms":107374,"temperature":0.7,"pith_summary":"What the paper tries to establish: for SU(N)-symmetric Ginzburg-Landau models with three or more components coupled to an Abelian gauge field, raising the gauge-field coupling high enough splits the single ordering transition into two. On cooling through the upper transition the Meissner effect disappears while order remains in the phase differences and density differences between components; the paper calls this a $CP^{{N-1}}$-neutral phase and shows that for N>2 it cannot be mapped onto an O(M) model, unlike the SU(2) and U(1)^N cases. The same systems, although not superconductors or superfluids in the usual sense, respond to an external magnetic field at low temperature by forming lattices of composite, non-topological vortices. The interest is that a gauge-theoretic model can pass through a genuine intermediate state with composite neutral order, which is a new finite-temperature phase-structure scenario for N>2.","feed_headline":"SU(N) models split into a neutral phase with no Meissner effect","feed_subtitle":"For N>2, the intermediate CP^{N-1} state breaks only phase-difference and density-difference symmetries.","key_machinery":"The load-bearing machinery is the rewriting of the Hamiltonian into charged, magnetic, and neutral terms, $$h = \\tfrac12 $j^{2}$ + \\tfrac12(\\nabla\\times A)^2 + \\sum_{i,j>i}|\\psi_i|^2|\\psi_j|^2(\\nabla\\varphi_{ij})^2 + \\tfrac12\\sum_i(\\nabla|\\psi_i|)^2,$$ together with the hard constraint $\\sum_i|\\psi_i|^2=1$. The constraint turns the neutral sector into a $CP^{N-1}$ target space, the space of N complex amplitudes of fixed total density modulo a common phase, whose degrees of freedom are the phase differences $\\varphi_{ij}$ and relative-density gradients. The paper's scenario is that integer-flux composite vortices, objects with winding in every component, are the cheapest excitations when all densities are nonzero; at large q they are tightly bound composites of fractional-flux pieces, so their proliferation removes the Meissner effect while leaving the phase-difference and relative-density sector ordered. The measured quantities are the dual stiffness $\\rho$, which vanishes in the Meissner state, and the helicity modulus $\\Upsilon$, which detects neutral order; their finite-size crossings locate the two transitions.","core_discovery":"The discovery claimed is that at sufficiently large gauge-field coupling q the SU(N) Ginzburg-Landau model in three dimensions, for N=3 and N=4 studied here (with N=2 for comparison), exhibits two separate phase transitions as temperature is lowered. Between the symmetric phase and the fully ordered low-temperature phase sits a $CP^{{N-1}}$-neutral phase: a state with no Meissner effect but a nonzero helicity modulus for phase-difference combinations, meaning spontaneous breaking only of relative-phase and relative-density symmetries. Because the hard constraint $\\sum_i|\\psi_i|^2=1$ makes the neutral sector a $CP^{N-1}$ target space, and because for N>2 the residual symmetry group is not the $S^{2}$ or $S^{1}$ type that appears in SU(2) models, the authors argue this intermediate phase cannot be represented as an O(M) model. The Monte Carlo study of magnetic response shows that in an external field the low-temperature state is a vortex lattice; for q=0 and N=3 the vortices group into triplets, and for q=1 in the two-component case a hexagonal lattice of half-quantum-flux objects appears, interpreted as a lattice of composite integer-flux objects with split cores.","pith_inferences":["If the CP^{N-1}-neutral phase survives when the hard density constraint is relaxed to a soft potential, it would be a finite-temperature example of composite order outside the O(M) duality description, offering a concrete test case for deconfined-criticality scenarios in lattice gauge theories.","Varying q continuously in the two-component model and tracking the six-peak magnetic structure factor would show whether the half-quantum vortex lattice persists beyond q=1 or crosses over to a conventional Abrikosov lattice as composite vortices merge.","The N=4 response to an external field is not reported; if the N=3 triplet grouping becomes a quadruplet grouping, that would confirm the composite-vortex interpretation is generic in N."],"forward_implications":["For the studied N=3 and N=4 cases at large q, the phase diagram has two transitions, and the paper's mechanism implies the same split should occur for all N>2 at sufficiently strong coupling.","The intermediate CP^{N-1}-neutral phase has no Meissner effect but retains phase-difference and relative-density order, so it is a distinct thermodynamic state with two heat-capacity peaks.","Because the neutral sector is CP^{N-1} and not O(M) for N>2, this composite phase falls outside the SU(2) and U(1)^N paired-phase classifications.","In an external magnetic field at low temperature the systems form vortex lattices even though individual integer-flux vortices are not energetically stable; at q=1 the two-component lattice is a hexagonal array of half-quantum-flux objects."],"supporting_citations":[{"why":"Establishes the SU(2) strong-coupling two-transition case with a neutral intermediate state, the baseline that the N>2 phase diagrams extend.","marker":"[25]"},{"why":"Provides the U(1)^N analysis in which integer-flux composite vortices dominate, giving the mechanism for why the charged transition shifts little at large charge.","marker":"[13]"},{"why":"States the general principle that proliferation of composite vortices can disorder individual phases while preserving products, on which the intermediate-phase scenario rests.","marker":"[22]"},{"why":"Documents the instability of single integer-flux vortices in type-2 SU(2) gauge theories, defining why the vortices invoked here are non-topological.","marker":"[28]"},{"why":"Shows that an SU(2) gauge system still forms a vortex lattice in an external field, the direct precedent for the vortex-lattice result.","marker":"[29]"},{"why":"Supplies the neutral SU(2) vortex patterns under rotation used as the q→0 comparison for the two- and three-component patterns.","marker":"[31]"},{"why":"Introduces the dual-stiffness observable whose finite-size crossings locate the Meissner transition in the phase diagrams.","marker":"[40]"}],"fun_headline_variants":["SU(N) models split into neutral phase without Meissner effect","For N>2, SU(N) gauge theory hosts CP^{N-1} neutral phase","Vortex lattice emerges in low-T SU(N) gauge theory","Two phase transitions found in SU(N) gauge models","SU(N) neutral phase escapes O(M) model mapping"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model fixes the total density $\\sum_i|\\psi_i|^2=1$ exactly at every lattice site, and the transition temperatures are read from finite-size crossings; if that hard constraint is relaxed, or if the charged transition is weakly first order, the intermediate $CP^{N-1}$-neutral phase and its claimed non-O(M) character could shift or disappear.","fun_headline_variants_meta":{"raw":{"variants":["SU(N) models split into neutral phase without Meissner effect","For N>2, SU(N) gauge theory hosts CP^{N-1} neutral phase","Vortex lattice emerges in low-T SU(N) gauge theory","Two phase transitions found in SU(N) gauge models","SU(N) neutral phase escapes O(M) model mapping"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001252,"raw_usage":{"total_tokens":5159,"prompt_tokens":997,"completion_tokens":4162,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":4071}},"tokens_in":613,"tokens_out":4162,"duration_ms":29157,"temperature":1.0,"reasoning_tokens":4071,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:33:21.320749+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Simulate the N=3, q=6 model with a soft total-density potential such as $\\lambda(\\sum_i|\\psi_i|^2-1)^2$ at finite $\\lambda$ and check whether the $L\\Upsilon$ crossing for the neutral transition persists as $\\lambda\\to\\infty$; separately, measure the latent heat or Binder cumulant at the charged transition for q=7, because a clear first-order signature would invalidate the continuous-transition crossings used to locate the split.","supporting_citations":[{"cited_title":"Smiseth , author E","cited_arxiv_id":null,"evidence_quote":"Provides the U(1)^N analysis in which integer-flux composite vortices dominate, giving the mechanism for why the charged transition shifts little at large charge."},{"cited_title":"Achucarro \\ and\\ author T","cited_arxiv_id":null,"evidence_quote":"Documents the instability of single integer-flux vortices in type-2 SU(2) gauge theories, defining why the vortices invoked here are non-topological."},{"cited_title":"Garaud \\ and\\ author E","cited_arxiv_id":null,"evidence_quote":"Shows that an SU(2) gauge system still forms a vortex lattice in an external field, the direct precedent for the vortex-lattice result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the neutral SU(2) vortex patterns under rotation used as the q→0 comparison for the two- and three-component patterns."}],"review_version":1}