{"id":"1f3ad5c5-1856-47fb-9934-460ee8436dce","arxiv_id":"2608.07990","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Quantum Monte Carlo simulations of an SO(3)-symmetric bilayer SLAC fermion model with six Dirac cones find a continuous Gross-Neveu critical point into a partially gapless ordered phase and a first-order transition into an inter-layer superconducting state.","lead":"This paper simulates a special two-layer electron model and finds two phase transitions: first into an ordered state where part of the electrons remain mobile, then into a superconducting state. The result gives a new platform for studying exotic quantum critical points and pairing.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The strong-coupling SC phase identification rests on a finite-size m^2_SC that is never extrapolated to the thermodynamic limit.","rationale":"The paper is a sign-problem-free projector DQMC study of a bilayer SLAC fermion model. The continuous DSM-to-SO(3)-broken transition is supported by correlation-ratio crossings, scaling collapse, and power-law fits of m^2_SO(3) and the fermion correlator. These elements give reasonable evidence for a continuous transition in the expected universality class, even though the extracted η exponents deviate from RG predictions; this is a known issue also seen for N=12, and the text acknowledges it. The most fragile part of the central claim is the strong-coupling superconducting phase. The evidence there is an abrupt onset of the finite-size structure factor m^2_SC (Eq. 8) and a simultaneous drop in m^2_SO(3); no thermodynamic-limit scaling is presented. This matters because m^2_SC is a volume-averaged correlation: in a disordered phase it generically decays as L^{-2}, while in a true SC phase it tends to a positive constant. The paper demonstrates that competing pairing channels decrease with L (Fig. S1), but does not show that the selected channel grows or saturates with L. Without such an extrapolation, the 'sudden onset' is consistent with a finite-size artifact, and the claim of SO(3)-symmetric superconductivity—a headline result—is unverified. Because the same finite-size data are used to argue for a fluctuation-induced first-order transition, both parts of the strong-coupling story hinge on the same missing scaling analysis. This is exactly the condition that must hold for the central claim, and it is the least secure. Therefore the reader's conditional verdict is appropriate; no change is needed.","tokens_in":13263,"tokens_out":16676,"duration_ms":167226,"concrete_test":"At a fixed J well above J_c2 (e.g., J=2.0), compute m^2_SC for L=9,11,13,15,17 (and larger if reachable) and plot m^2_SC against 1/L^2. If m^2_SC extrapolates to a positive constant as L→∞, superconducting order is established; if it vanishes, the strong-coupling phase is not SC and the first-order interpretation must be revisited.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In the section 'SC phase in the strong-coupling region', the identification of the strong-coupling phase as an SO(3)-symmetric superconductor rests on the abrupt onset of the finite-size structure factor m^2_SC defined in Eq. (8). For a local pairing operator in d=2, a disordered phase gives m^2_SC ~ O(L^-2) (short-range correlations summed over the volume and divided by L^4), whereas a phase with off-diagonal long-range order gives m^2_SC -> constant > 0 as L -> infinity. The paper does not report the L-dependence of m^2_SC at fixed J > J_c2, nor a thermodynamic-limit extrapolation. Fig. 4(a) shows m^2_SC only for L=9 and L=13 (per the legend); the competing channels in Fig. S1 are explicitly shown to decrease with L, but no analogous demonstration is given for m^2_SC. Thus the 'sudden onset' could be a finite-size pairing enhancement rather than genuine SC order. The claim of a simultaneous first-order transition at J_c2 inherits this fragility, since the discontinuity in m^2_SO(3) alone does not identify the ordered state. This is the least secure condition for the paper's central claim that the model realizes a superconducting phase.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a bilayer SLAC fermion model with SO(3) flavor symmetry using sign-problem-free projector determinant quantum Monte Carlo. It claims a continuous quantum phase transition from a Dirac semimetal to an SO(3)-broken Dirac semimetal, assigned to the Gross-Neveu-SO(3) universality class with N=6 irreducible Dirac cones, with extracted critical exponents ν=0.81(5), η_b=0.44(2), η_f=0.051(2). At stronger coupling, the model is claimed to enter an inter-layer SO(3)-symmetric superconducting phase via a first-order transition. The paper includes an unrestricted Hartree-Fock-Bogoliubov mean-field analysis to identify the pairing channel.","tokens_in":13551,"tokens_out":13037,"duration_ms":131211,"significance":"If the claims hold, the paper would provide the first numerical realization of Gross-Neveu-SO(3) criticality for N=6 and a new microscopic model for SO(3)-symmetric superconductivity, extending the earlier N=12 study. The QMC method is sign-problem-free and the model construction is physically motivated. The paper is transparent about the deviations of anomalous dimensions from RG predictions and about the mean-field/QMC discrepancy at the second transition. However, the strength of the superconductivity claim is limited by the absence of a thermodynamic-limit analysis of the pairing structure factor.","major_comments":[{"comment":"The identification of the strong-coupling phase as a superconductor rests on the finite-size structure factor m^2_SC, defined in Eq. (8), shown only for L=9 and L=13. In d=2, a disordered phase gives m^2_SC ~ O(L^{-2}), while ODLRO gives m^2_SC → const > 0. Without a thermodynamic-limit extrapolation of m^2_SC at fixed J > J_c2, the abrupt onset could be a finite-size pairing enhancement rather than true SC order. Please provide m^2_SC as a function of L for several J > J_c2 and an extrapolation, or otherwise demonstrate that the order parameter survives in the thermodynamic limit.","section":"SC phase in the strong-coupling region, Fig. 4(a)"},{"comment":"The first-order character of the transition at J_c2 is inferred from the sharp drop in m^2_SO(3) and the simultaneous onset of m^2_SC. For a first-order quantum phase transition, one expects a discontinuity that sharpens with system size, typically evidenced by a double-peaked histogram, a Binder cumulant that develops a minimum, or a clear L-dependence of the transition width. The present finite-size data (L up to 15 for SO(3) and L up to 13 for SC) do not rule out a very sharp continuous transition. Please provide additional diagnostics (e.g., energy histograms or L-dependence of the jumps) to support the claim of a first-order transition.","section":"SC phase in the strong-coupling region, Fig. 4(a)"},{"comment":"The assignment of the transition to the GN-SO(3) universality class relies on matching the correlation-length exponent ν and the anomalous dimensions η_b and η_f to RG predictions. While ν=0.81(5) is consistent with the 1/N expansion, the quoted η_b=0.44(2) and η_f=0.051(2) are both stated to deviate substantially from the RG results. To make the universality-class assignment convincing, please provide a quantitative table comparing all extracted exponents with the RG predictions, and discuss whether the deviations could be a consequence of limited system sizes or subleading corrections. In particular, consider fitting with correction-to-scaling terms before asserting a precise match.","section":"SO(3)-broken phase and GN-SO(3) universality class, after Eq. (5) and comparison with Ref. [60]"}],"minor_comments":[{"comment":"The title 'Metallic Gross-Neveu criticality and superconductivity on theSO(3)SLAC fermion' is missing spaces; please correct the typo.","section":"Title"},{"comment":"Δq is introduced as 'the minimum lattice momentum'; please specify the actual momentum components used in the simulation.","section":"Eq. (3)"},{"comment":"The fitting procedure for R_SO(3) selects the fifth-order polynomial based on the smallest root-mean-square residual; since higher-order polynomials always reduce residuals, please report the goodness of fit (e.g., chi-square per degree of freedom) or the spread of J_c1 and ν across different n_max.","section":"Fig. 2 and fitting procedure"},{"comment":"The abstract says the exponents are 'precisely extracted,' but the text notes substantial deviations from RG predictions; please consider a more cautious wording.","section":"Abstract and Introduction"},{"comment":"The operator ordering for m^2_SC1 is not immediately transparent; a brief explanation of the ordering would improve readability.","section":"Supplemental Material, Eq. (S12)"},{"comment":"The HFB calculation uses L=71 with a (π,π) twist, while QMC uses periodic boundary conditions up to L=17; please state whether boundary conditions could affect the comparison.","section":"Supplemental Material, Mean-field analysis"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is potentially interesting and the QMC methodology is sound. My main concern is that the superconducting phase and the first-order transition are not yet established at the thermodynamic level. The authors should be encouraged to perform the additional finite-size scaling analysis for m^2_SC. I would also suggest toning down the claim of 'precise' exponent extraction given the noted deviations from RG. The paper is, however, within the scope of the journal and the issues are addressable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The first thing you should know: this paper delivers the first numerical realization of metallic Gross-Neveu criticality at N=6, and it does so with a clean, sign-problem-free QMC setup. The model is new, the SO(3)-symmetric bilayer construction is sensible, and the transition out of the Dirac semimetal is supported by a convincing scaling collapse. The extracted exponents, nu=0.81(5), eta_b=0.44(2), eta_f=0.051(2), are honestly compared with RG predictions, and the authors acknowledge the deviations. That is more than many papers in this area do.\n\nWhere the paper is weaker is the superconducting phase. The identification of an inter-layer SO(3)-symmetric SC state rests on the abrupt onset of m_SC^2 in Fig. 4(a), but the figure shows only L=9 and L=13, and there is no finite-size scaling of m_SC^2 to show the order parameter survives in the thermodynamic limit. The competing channels in Fig. S1 do decay with L, but the main channel is not given the same test. Given that local pairing in d=2 gives m_SC^2 ~ O(L^-2) in a disordered phase, the sudden onset could be a finite-size pairing enhancement rather than true SC order. This is the least secure claim in the paper. The first-order transition at J_c2 inherits that fragility, since it is identified from the discontinuity in m_SO(3)^2 and the simultaneous onset of m_SC^2, both at finite size.\n\nOn the universality class assignment: the continuous transition and the SO(3) structure make the GN-SO(3) class plausible, but the exponent agreement is mixed. nu matches the 1/N expansion reasonably, but eta_b and eta_f are both smaller than RG predictions. The authors are upfront about this, and similar discrepancies exist for N=12, so the claim is not overreaching, but it is not definitive either.\n\nOverall, this is a well-executed numerical study of a new model, with a transparent discussion of its limitations. The SC phase needs a genuine thermodynamic-limit analysis—at minimum an L-scaling study of m_SC^2 at fixed J>J_c2 and a first-order diagnostic like a histogram or Binder cumulant. I would send it to peer review, because the N=6 GN-SO(3) result is new and the model will be useful to the community, but the SC claim should be substantially strengthened before publication.\n\nWho should read it? Anyone working on Dirac quantum criticality, sign-problem-free fermion QMC, or unconventional pairing in Dirac systems. It deserves a serious referee, and the referee should push on the SC phase.","headline":"A solid QMC study of a new N=6 SO(3) Dirac model with a clean continuous transition and an intriguing but unproven SC phase; the SC claim needs a thermodynamic-limit extrapolation before it can be accepted.","tokens_in":14094,"tokens_out":2083,"would_cite":true,"duration_ms":24467,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a bilayer SO(3)-symmetric SLAC fermion model hosts two successive transitions: a continuous Gross-Neveu-SO(3) quantum critical point with N=6 that gaps only a subset of Dirac cones, followed by a first-order…","keywords":["Gross-Neveu-SO(3) universality class","Dirac semimetal","SLAC fermion","quantum Monte Carlo","SO(3) symmetry breaking","inter-layer superconductivity","first-order transition","critical exponents"],"falsifier":"Compute $m_SC^{2}$ of Eq. (8) on a sequence of larger linear sizes L at fixed J above J_c2 and check whether it saturates to a nonzero value or decays with L; if it decays, the strong-coupling phase is not a true SO(3)-symmetric superconductor. Additionally, scan J back and forth across J_c2 to look for hysteresis or a latent-heat peak, which would confirm the claimed first-order transition.","tokens_in":13061,"feed_emoji":"⚛️","tokens_out":4664,"duration_ms":53955,"temperature":0.7,"pith_summary":"This paper introduces a bilayer model of SLAC fermions with SO(3)-symmetric interactions and uses quantum Monte Carlo to map its ground-state phase diagram. It claims that as inter-layer coupling increases, the six-cone Dirac semimetal first undergoes a continuous transition into an SO(3)-broken semimetal in which only two thirds of the Dirac cones become gapped, with two cones remaining gapless. This transition is identified as a realization of the Gross-Neveu-SO(3) universality class with N=6, with exponents nu=0.81(5), eta_b=0.44(2), and eta_f=0.051(2). At stronger coupling, the SO(3) order drops abruptly while an inter-layer, intra-flavor, SO(3)-symmetric superconducting structure factor rises, evidence for a direct first-order transition rather than a deconfined critical point. The work matters because it extends metallic Gross-Neveu criticality beyond the standard Gross-Neveu-Yukawa paradigm and provides a microscopic platform for SO(3)-symmetric superconductivity.","feed_headline":"SO(3) Dirac model shows new criticality, then superconducts","feed_subtitle":"Monte Carlo finds the N=6 Gross-Neveu-SO(3) critical point, then a first-order jump into inter-layer pairing.","key_machinery":"The central object is the bilayer SLAC fermion Hamiltonian of Eq. (1), whose layer-resolved hopping takes conjugate amplitudes t_R and t_R^* so the two layers carry opposite chiralities. The model is sign-problem-free under the antiunitary symmetry T=i tau_y K, which pairs the two layers and guarantees nonnegative fermion weights in the auxiliary-field quantum Monte Carlo. The order parameter for SO(3) breaking is the staggered structure factor S(0) built from the flavor generators K^$\\alpha$, and the transition is located through the dimensionless correlation ratio R_SO(3)=1-S($\\Delta$ q)/S(0), whose finite-size scaling collapse yields J_c1 and nu. The anomalous dimensions eta_b and eta_f follow from power-law scaling of the squared order parameter m_SO(3)^2 and the fermion correlation G_f at criticality. The superconducting channel is probed by the inter-layer pairing structure factor $m_SC^{2}$ of Eq. (8), and its simultaneous jump with the drop in m_SO(3)^2 is used to diagnose the first-order transition.","core_discovery":"The central claim is that an SO(3)-symmetric bilayer of SLAC fermions with N=6 irreducible Dirac cones exhibits a continuous Dirac-semimetal-to-SO(3)-broken transition governed by the Gross-Neveu-SO(3) universality class, and that at larger coupling the system enters an SO(3)-symmetric superconducting state through a fluctuation-induced first-order transition. Unlike conventional Gross-Neveu transitions where all Dirac cones gap out, here the SO(3) order gaps only two of the three flavor copies, leaving one flavor (two cones) gapless, a scenario the paper calls metallic Gross-Neveu criticality. The paper extracts the critical point J_c1=0.803(4) and critical exponents from finite-size scaling of the correlation ratio, the SO(3) order parameter, and the fermion correlation function. It further shows that at J_c2 the SO(3) and U(1) order parameters drop sharply while the intra-flavor inter-layer pairing structure factor $m_SC^{2}$ onsets abruptly, ruling out a deconfined quantum critical point and supporting a first-order boundary. The superconducting order preserves SO(3) symmetry and is consistent with an unrestricted Hartree-Fock-Bogoliubov analysis identifying same-flavor inter-layer pairing as the dominant channel.","pith_inferences":["A natural extension would be to test whether the N=6 Gross-Neveu-SO(3) fixed point and the N=12 one flow to the same large-N analytical limit or whether finite-N corrections are unusually large, given the substantial discrepancy with RG predictions for anomalous dimensions.","The paper's case for a true superconducting phase would be strengthened by an independent thermodynamic-limit scaling analysis of m_SC^2 at fixed J > J_c2; without it, the strong-coupling signal could in principle be a finite-size pairing enhancement rather than a genuine ordered state.","One could check the robustness of the first-order transition by varying the trial-state construction or boundary conditions, since the flavor-balanced closed-shell trial state explicitly preserves SO(3) but may affect finite-size estimates near the transition.","If the SO(3)-symmetric superconducting state is realized in cold-atom or synthetic-layer platforms, it would constitute a rare example of pairing that does not break the full spin-rotation symmetry, which may host distinct topological or thermodynamic signatures."],"forward_implications":["If the continuous transition is in the Gross-Neveu-SO(3) universality class with N=6, then metallic Gross-Neveu criticality is not restricted to N=12 and can occur when the ordered phase leaves a fraction of Dirac cones gapless.","The extracted exponents, especially the anomalous dimensions eta_b=0.44(2) and eta_f=0.051(2), provide a direct numerical benchmark for future analytical or conformal-bootstrap studies of the SO(3) Gross-Neveu fixed point.","The abrupt simultaneous drop in m_SO(3)^2 and rise in m_SC^2 at J_c2 implies a direct first-order transition from the SO(3)-broken semimetal to the SO(3)-symmetric superconductor, with no intervening deconfined critical point as seen in the N=12 model.","The superconducting state is channel-selective: intra-flavor inter-layer pairing dominates, while two competing inter-flavor channels show no comparable onset, so the strong-coupling phase is not a generic pairing enhancement.","The comparison with Hartree-Fock-Bogoliubov results shows that quantum fluctuations of the SO(3) order parameter preempt a mean-field coexistence region, so the first-order character is a quantum-fluctuation-driven effect rather than a mean-field feature."],"supporting_citations":[{"why":"Proposes the Gross-Neveu-SO(3) universality class that the paper tests and extends to N=6.","marker":"[14]"},{"why":"Provides the N=12 bilayer SLAC model, the sign-problem-free construction, and the DQCP scenario that this paper compares against.","marker":"[22]"},{"why":"Supplies the polynomial finite-size scaling fitting procedure used to extract J_c1 and nu from the correlation ratio.","marker":"[34]"},{"why":"The projector determinant quantum Monte Carlo method used for all ground-state simulations.","marker":"[58]"},{"why":"Renormalization-group and 1/N expansion predictions for the critical exponents that the numerical results are compared with.","marker":"[60]"},{"why":"Theoretical framework for quantum-fluctuation-induced first-order transitions that is invoked to interpret the preemption of the mean-field coexistence region.","marker":"[61–64]"}],"fun_headline_variants":["Metallic Gross-Neveu criticality then SO(3) superconductivity","N=6 Gross-Neveu-SO(3) critical point, first-order SC","Beyond GNY: metallic criticality and SO(3) pairing","SO(3) bilayer shows new criticality, then superconducts","Fluctuation-driven first-order transition into SO(3) SC"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the abrupt onset of the inter-layer pairing structure factor $m_SC^{2}$ at J > J_c2 reflects a true thermodynamic superconducting order; the paper does not show $m_SC^{2}$ surviving the thermodynamic limit, so the strong-coupling phase could instead be a finite-size pairing enhancement.","fun_headline_variants_meta":{"raw":{"variants":["Metallic Gross-Neveu criticality then SO(3) superconductivity","N=6 Gross-Neveu-SO(3) critical point, first-order SC","Beyond GNY: metallic criticality and SO(3) pairing","SO(3) bilayer shows new criticality, then superconducts","Fluctuation-driven first-order transition into SO(3) SC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000653,"raw_usage":{"total_tokens":3065,"prompt_tokens":1087,"completion_tokens":1978,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":703,"completion_tokens_details":{"reasoning_tokens":1881}},"tokens_in":703,"tokens_out":1978,"duration_ms":16353,"temperature":1.0,"reasoning_tokens":1881,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:35:18.178408+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute $m_SC^{2}$ of Eq. (8) on a sequence of larger linear sizes L at fixed J above J_c2 and check whether it saturates to a nonzero value or decays with L; if it decays, the strong-coupling phase is not a true SO(3)-symmetric superconductor. Additionally, scan J back and forth across J_c2 to look for hysteresis or a latent-heat peak, which would confirm the claimed first-order transition.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the N=12 bilayer SLAC model, the sign-problem-free construction, and the DQCP scenario that this paper compares against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the polynomial finite-size scaling fitting procedure used to extract J_c1 and nu from the correlation ratio."},{"cited_title":"Quantum monte carlo methods on lat- tices: The determinantal approach,","cited_arxiv_id":null,"evidence_quote":"The projector determinant quantum Monte Carlo method used for all ground-state simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Renormalization-group and 1/N expansion predictions for the critical exponents that the numerical results are compared with."}],"review_version":1}