{"id":"8ead3933-e9c3-4754-8dc6-846c5ed6c92b","arxiv_id":"2607.25767","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Flat membrane-Dirac systems asymptotically decouple at charge neutrality; finite-momentum rippling realizes either Wilson-Fisher criticality with spectator fermions or, when the ripple carries Kekulé mass quantum numbers, chiral-XY Gross-Neveu-Yukawa criticality.","lead":"Crystalline membranes like graphene combine slow bending vibrations with fast relativistic electrons. This paper shows the flat sheet decouples from the electrons at long wavelengths, and sorts ripple phase transitions into two classes: an electron-blind elastic transition, or a hybrid electronic-structural Gross-Neveu-Yukawa transition when symmetry allows.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Flat-phase decoupling rests on an asserted power-counting bound in Appendix D; an explicit two-loop check would settle the central classification.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the flat-phase decoupling rests on the Appendix D power-counting bound that no s^0 logarithm is generated in the fast region, a point asserted rather than computed. My independent read of the paper confirms this is the most fragile link in the central classification. The finite-Q Wilson–Fisher spectator result is robust given the exact scaling dimension y_λ = 1/ν − D and the accepted XY exponent; the GNY route is explicitly conditional on broken mirror symmetry or a mirror-even bond component, and the paper discloses that limitation. The model restriction to the local scalar channel, while important, narrows the claim rather than threatening its internal validity. The Appendix D assertion, by contrast, is internal to the argument: if wrong, the flat-phase asymptotic decoupling — the first half of the classification — is false. No error has been demonstrated, and the paper's one-loop calculations are internally consistent, so the appropriate disposition remains CONDITIONAL: acceptance should require an explicit check of the mixed vertex correction or a more rigorous proof of the s^0 suppression. The reader's verdict already reflects this, so I recommend no change.","tokens_in":21480,"tokens_out":5617,"duration_ms":65618,"concrete_test":"Compute explicitly the order-g^3 (two-loop) correction to the 1PI vertex Γ^1PI_{hh\\barψψ} in eq. (D.13) for external flexural kinematics Ω ∼ q^2, using the full Dirac propagator without the slow-flexural s expansion. Project onto the tree tensor structure (p1·p2)γ0 at zero external fermion momentum and extract the coefficient of ln(Λ/µ). If the coefficient vanishes as s→0 (i.e., is O(s)), eq. (D.15) holds. If it is O(s^0), the flat-phase coupling receives a marginal fermionic correction, and the decoupling claim in eq. (3.17) fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The first half of the central classification — asymptotic decoupling of Dirac fermions in the flat membrane phase — depends on the claim in Appendix D that the projected local scalar vertex receives no infrared s^0 correction from the fermions. The order-g^2 diagram vanishes by oddness, so the first candidate correction is order g^3 (i.e., O(ḡ^2) after extracting the tree factor), and the paper bounds it as δg_R/g_R = O(ḡ^2 s) with s = ω/(v_F k), concluding γ_g^(ψ) = O(α_e) and hence irrelevance. This is a power-counting assertion, not a computed leading correction: the slow-flexural expansion used to obtain the s suppression breaks down in the fast region ω ∼ v_F k, where s is order one and the flexural propagator changes from k^-4 to k^-2. The statement that this fast region 'cannot generate an s^0 logarithm' is plausible but not demonstrated. If an unsuppressed logarithm did exist, it would contribute a nonvanishing term to d ln ḡ/dℓ, potentially halting the logarithmic decrease in eq. (3.17) and invalidating the asymptotic decoupling. The argument also defers derivative-operator mixing (Γ_der) to a separate composite-operator analysis, so the bound applies only to the zero-momentum projection. Thus the classification's flat-phase leg rests on an unverified technical estimate.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a low-energy field theory for a crystalline membrane coupled to massless Dirac fermions and proposes a classification of membrane–Dirac criticality. In the flat phase, the local scalar strain–density coupling is argued to be marginally irrelevant: rotational invariance fixes its pure-membrane renormalization, while fermionic feedback is suppressed by the running coupling α_e, so the two sectors asymptotically decouple. At finite momentum, two scenarios are distinguished: if the ripple order parameter does not share quantum numbers with any low-energy Dirac bilinear, the transition is governed by the Wilson–Fisher fixed point with spectator Dirac fermions; if it does share them (including horizontal-reflection parity), the transition belongs to the chiral-XY Gross–Neveu–Yukawa universality class, with one-loop exponents ν=4/5, η_Φ=2/3, η_ψ=1/6 and velocity locking with exponent ω_ζ=7ε/36. The GNY counterterms are re-derived in Appendix F rather than merely quoted.","tokens_in":21791,"tokens_out":17854,"duration_ms":178782,"significance":"If the central claim holds, the paper provides a useful organizing principle for a class of problems that mixes z=1 Dirac fermions with z=2 flexural phonons, and it sharpens the distinction between asymptotic flat-phase stability and finite-momentum order. The re-derivation of the one-loop membrane and GNY flows, the explicit symmetry bookkeeping for the Kekulé mass matrices, and the concrete predictions for velocities and exponents are strengths. The classification is, however, stated more broadly than the actual model, and the flat-phase leg relies on a power-counting estimate in Appendix D that is not backed by an explicit computation. These issues do not invalidate the overall picture, but they need to be addressed before the classification can be taken as established.","major_comments":[{"comment":"The flat-phase decoupling rests on the estimate δg_R/g_R = O(ḡ²s), with the statement that the fast region ω∼v_F k “cannot generate an s⁰ logarithm” in the projected local vertex. This is asserted but not demonstrated: the order-g³ vertex correction is not computed, and the slow-flexural expansion (D.11) breaks down precisely in the fast region, where s∼1 and the flexural propagator changes from k⁻⁴ to k⁻². Since eq. (3.24) and the claim γ_g^(ψ)=O(α_e) are load-bearing for the first half of the classification, the authors should either supply an explicit leading calculation of the projected local vertex or clearly state the weaker, robust statement (e.g., what happens if only O(ḡ²) is known; my estimate is that the marginal-log stiffening would still dominate, but this should be said in the paper).","section":"§3.2 and Appendix D, Eqs. (D.13)–(D.16)"},{"comment":"The model includes only the local trace coupling H_T ρψ; the momentum-projected deformation potential and the shear/pseudogauge valley coupling are explicitly excluded (“not included here”, §2.1), and Γ_der is deferred to a separate composite-operator analysis. Nevertheless the abstract and introduction present a “field-theoretic classification of membrane–Dirac criticality” with no such qualifier. The finite-momentum classification in §§4–5 is vulnerable to the omitted terms: at finite Q the pseudogauge coupling is linear in the ripple amplitude and could generate a symmetry-allowed Yukawa-type coupling not captured by the density coupling in eq. (4.7). The central claim should be reformulated as a classification within the local-scalar sector, or the omitted interactions need to be analyzed, at least by power counting.","section":"§2.1 and §6"}],"minor_comments":[{"comment":"For N_f=N_b=2 and ε=1, Eq. (5.16) gives 1/ν=4/5, i.e. ν=5/4, whereas Eq. (5.17) gives ν=4/5 by substituting ε=1 into the expanded series. This is the usual ε-expansion extrapolation ambiguity, but the text should state explicitly that 4/5 is obtained from the truncated series, not from the fixed-point relation (5.16).","section":"§5.2, Eqs. (5.16)–(5.17)"},{"comment":"The full functional derivation of the rotational Ward identity is said to be “not reproduced here.” Since the zero-momentum tension-insertion result is central to §3, it would help to give the complete argument or at least a clearer reference to the original derivation, rather than only a one-loop verification in Appendix D.","section":"Appendix B"},{"comment":"The velocity-locking analysis is linearized about ζ=1 and shows stability of the equal-velocity point. The statement that a generic anisotropic state is “expected” to flow to a common terminal velocity is stronger than what is derived; the text should mark this as an expectation based on the isotropic subspace.","section":"§5.3"},{"comment":"There are several typographical issues, including missing spaces (“Hereκ b” in §3.2, “Themixingtherefore” after eq. (5.9), “correspoding” in §5) and inconsistent spacing around cross-references. These should be cleaned up.","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely publishable after revision. The main substantive work is to close the Appendix D gap and to align the scope of the classification with the model actually analyzed. The GNY part is standard but competently re-derived; the membrane–Dirac combination is the novel element."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a serious paper and the central classification is worth taking seriously. The flat phase is claimed to decouple asymptotically because the local scalar coupling flows logarithmically to zero via the Ward identity, while the fermionic feedback alpha_e dies as a power law. For finite-momentum rippling, the paper separates Wilson-Fisher with spectator Dirac fermions from a chiral-XY GNY transition, depending on whether a mass-type bilinear shares the ripple's quantum numbers. That classification is new--earlier work on electron-induced rippling did not have this controlled RG structure.\n\nThe one-loop machinery is solid: the RG equations are internally consistent, the fixed points and exponents match earlier GNY results, and Appendix F re-derives the counterterms rather than quoting them. The citation pattern is clean--the membrane flows come from Mauri-Katsnelson and Kats-Lebedev, and the GNY exponents are independently derived or re-derived here. The authors also disclose the relevant limitations: the mirror-parity obstruction in graphene, the isolated-+/-Q condition, and the omission of the momentum-projected deformation potential and the shear/pseudogauge channel.\n\nThe soft spot is where the reader's take and stress-test land: the flat-phase decoupling depends on Appendix D's assertion that the fast region (omega ~ v_F k) cannot generate an s^0 logarithm in the projected local vertex. That is a power-counting claim, not a computed leading correction. The order-g^2 diagram vanishes by oddness, so the first candidate is order g^3; the slow-flexural expansion gives a factor s, but in the fast region s is O(1) and the k^-2 propagator argument is plausible without being verified. If an unsuppressed logarithm existed, the dimensionless coupling would stop its decrease and asymptotic decoupling would fail. This does not sink the paper--the bound is likely right--but it is precisely the kind of estimate that should be checked explicitly before the classification is taken as definitive.\n\nThe model restriction also matters. The flat-phase result applies only to the leading local scalar coupling; the full deformation potential and the valley-current channel are outside the analysis. That is honestly flagged, but it strictly limits what the 'stable flat phase' claim covers.\n\nIn short: the paper advances the field, the classification is useful, and the one-loop content is trustworthy. The weak spot is narrow and the authors have already named it; now it needs a strong justification or an explicit two-loop computation. Send it to a serious referee. Condensed matter theorists working on graphene, membranes, or Dirac criticality should read it.","headline":"Conditional accept: the classification is new and the one-loop RG holds up, but the flat-phase decoupling rests on an unverified power-counting bound that deserves a two-loop check.","tokens_in":22342,"tokens_out":4555,"would_cite":true,"duration_ms":48552,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["64.60.Fr","71.10.-w","81.05.ue"],"model":"deepseek-v4-flash","headline":"Flat graphene membranes decouple from electrons; ripples re-couple them","keywords":["crystalline membranes","Dirac fermions","renormalization group","flat phase","graphene","Kekulé ordering","Wilson–Fisher fixed point","Gross–Neveu–Yukawa"],"falsifier":"Compute the projected local vertex at P=0 to next order in the slow-flexural expansion for D=2: if delta g_R/g_R contains a term constant in the RG scale (an s^0 logarithm), the asymptotic decoupling flow fails. Alternatively, a careful numerical simulation of a graphene-like membrane–Dirac model measuring the effective strain–density coupling at long wavelengths could detect deviation from the predicted logarithmic decrease.","tokens_in":1344,"feed_emoji":"🌀","tokens_out":2390,"duration_ms":58458,"temperature":0.7,"pith_summary":"The paper establishes a classification of what happens at low energies when massless Dirac electrons live on a crystalline membrane whose bending modes are slow (dynamical exponent z=2) compared with the fermions (z=1). It argues that in the flat phase the local strain–density interaction becomes logarithmically weaker and electronic feedback dies as a power law, so the flat phase is stable. If the membrane ripples at a finite wavevector, the transition is ordinarily the purely bosonic Wilson–Fisher fixed point, with the Dirac fermions as spectators. Only when the ripple order parameter carries the same momentum and quantum numbers as a mass-type Dirac bilinear does the transition become a hybrid chiral-XY Gross–Neveu–Yukawa transition, with emergent velocity locking. A sympathetic reader cares because this organizes the possible electronic–structural critical phenomena in graphene-like membranes.","feed_headline":"Flat graphene membranes decouple from electrons; ripples re-couple them","feed_subtitle":"A symmetry rule decides when rippling is purely elastic and when electrons join the critical mode — and velocities lock.","key_machinery":"Three mechanisms carry the argument. First, the rotational/tension Ward identity of the pure membrane fixes the renormalization of the local trace-strain operator H_T, so the mixed coupling gbar=g/kappa_b inherits the membrane's marginal-irrelevant flow. Second, an expansion in s=omega_h/(v_F k), the ratio of flexural to Dirac energy scales, suppresses fermionic corrections by alpha_e=(N/16)gbar^2 s, making electronic feedback power-law irrelevant. Third, in the finite-Q sector, the scaling dimension y_lambda=1/nu-D and the Clifford-algebra matching of the ripple field with Kekulé mass matrices decide which fixed point controls the transition.","core_discovery":"The central claim is a field-theoretic classification: at charge neutrality, dynamical scale separation forces asymptotic decoupling in the flat phase, with the dimensionless coupling flowing as gbar(l)=gbar0[1+a y0 l]^{-theta}, fermionic feedback alpha_e(l) ~ e^{-l} l^{-3theta/2} -> 0, and no Landau damping. When the membrane is destabilized at finite momentum ±Q, the transition is either Wilson–Fisher with spectator Dirac fermions (since the fermion–boson coupling has scaling dimension y_lambda=1/nu-D < 0 in D=2) or, when symmetry permits a mass-type Dirac bilinear sharing the ripple's momentum and horizontal-reflection parity, a hybrid chiral-XY Gross–Neveu–Yukawa transition with one-loop","pith_inferences":["The flat-phase decoupling rests on the power-counting bound that no s^0 logarithm renormalizes the projected local vertex; a direct computation of the next correction to delta g_R/g_R would test this, and a hidden s^0 logarithm would open the possibility of weak-coupling instability at longer scales.","The paper explicitly omits the momentum-projected deformation potential and the shear/pseudogauge valley coupling; extending the classification to those channels could yield additional GNY universality classes, since they couple to different Dirac bilinears.","Asymptotic irrelevance does not rule out electron-induced rippling at intermediate scales, so the classification is compatible with experiments that observe ripples before the long-wavelength flat regime is reached.","For a nearly continuous ring of minima, the zero-temperature problem is marginal by power counting; working out the angularly resolved quartic flow would decide whether that transition is continuous or weakly first order."],"forward_implications":["The flat phase is stable against the local scalar deformation potential: the membrane stiffens logarithmically, flexural dynamics stay undamped, and the Dirac kinetic term retains canonical scaling.","At a finite-momentum ripple with no symmetry-matched Dirac bilinear, the transition is the ordinary bosonic Wilson–Fisher fixed point; the Dirac fermions are spectators at criticality.","In the ordered ripple phase the spectator coupling still shifts the Dirac-point energy relative to the chemical potential, producing a finite carrier density or chemical-potential shift.","When symmetry permits a mass-type Yukawa coupling, the critical point is chiral-XY Gross–Neveu–Yukawa; an electronic instability then induces a secondary structural distortion through Landau mixing.","At that hybrid critical point, fermionic and bosonic velocities lock to a common terminal value in the isotropic continuum limit, restoring emergent Lorentz symmetry."],"fun_headline_variants":["Symmetry decides if ripples drag electrons into graphene's critical mode","Flat graphene: electrons decouple; ripples re-couple them if symmetry allows","Ripples can re-couple graphene's electrons when mass bilinear matches","Hybrid critical point in graphene ripples hinges on reflection parity"],"cache_read_input_tokens":23552,"weakest_assumption_plain":"The flat-phase half of the classification assumes both that fermions produce no infrared s^0 logarithm in the projected local strain–density vertex (a power-counting bound in Appendix D) and that the model retains only the local scalar coupling, excluding momentum-projected deformation potential and shear/pseudogauge valley couplings.","fun_headline_variants_meta":{"raw":{"variants":["Symmetry decides if ripples drag electrons into graphene's critical mode","Flat graphene: electrons decouple; ripples re-couple them if symmetry allows","Ripples can re-couple graphene's electrons when mass bilinear matches","Hybrid critical point in graphene ripples hinges on reflection parity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000176,"raw_usage":{"total_tokens":1176,"prompt_tokens":845,"completion_tokens":331,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":252}},"tokens_in":589,"tokens_out":331,"duration_ms":4513,"temperature":1.0,"reasoning_tokens":252,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T01:32:30.238782+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the projected local vertex at P=0 to next order in the slow-flexural expansion for D=2: if delta g_R/g_R contains a term constant in the RG scale (an s^0 logarithm), the asymptotic decoupling flow fails. Alternatively, a careful numerical simulation of a graphene-like membrane–Dirac model measuring the effective strain–density coupling at long wavelengths could detect deviation from the predicted logarithmic decrease.","supporting_citations":[],"review_version":1}