{"id":"f0e5df6e-cf9c-4286-be15-420a6dc2569c","arxiv_id":"2504.19321","paper_version":1,"verdict":"REJECT","confidence":"LOW","novelty_score":7.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"Weak attractive interactions at a spin-polarized Van Hove singularity drive pair density wave or charge density wave order depending on the mass anisotropy, while repulsive interactions leave the system stable.","lead":"This paper predicts that in a spin and valley polarized two-dimensional metal tuned to a Van Hove singularity, weak attractive interactions between electrons drive either a pair density wave or a charge density wave, depending on the shape of the energy dispersion. The result offers a weak-coupling route to exotic superconductivity in rhombohedral multilayer graphene.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The printed one-loop coefficients in Appendix B have the wrong sign for 0<η<1/3 and are non-real for η>1/3, so the beta function does not support the claimed stability/instability dichotomy.","rationale":"The reader's stated weakest assumption was the existence of the fully spin- and valley-polarized normal state, but the reader's rationale identified the Appendix B sign problem as the load-bearing error. I agree with that rationale: the printed formulas give a(η)<0 where they are real and are non-real elsewhere, so Eq. (5) does not yield the claimed stable-repulsive/unstable-attractive dichotomy. The polarization assumption is explicitly acknowledged in Sec. IV B and is a limitation, not an internal inconsistency; the beta-function sign is more directly fatal to the central claim. Since this confirms the reader's verdict of REJECT with low confidence, the verdict should remain unchanged. A simple sign correction could in principle rescue the paper, but the written manuscript as submitted does not support its main conclusion.","tokens_in":18494,"tokens_out":6722,"duration_ms":72990,"concrete_test":"Evaluate the angular integrals in Eqs. (B8) and (B10) numerically at η=0.1, with their stated prefactors (2m√η/16π^2 and 2m√η/32π^2, respectively), to obtain a_{ZS'} and a_{BCS} directly from the loop calculation; compare these numbers with Eqs. (B1) and (B2). Also evaluate Eqs. (B1) and (B2) at η=0.5 to check whether they are real. If the numerically integrated a is positive while the closed forms are negative or non-real, the manuscript's sign claim fails; if the numerical integration confirms a<0, then the stability/instability conclusion is reversed.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central result rests on Eq. (5), which states dot g = -a(η)g^2 with a(η)>0. Substituting the printed Appendix B expressions: for η<1/3, the logarithm in Eq. (B1) is positive (its argument is greater than 1), so a_{ZS'}(η) = -[positive constant]×log(...) is negative; the square root in Eq. (B2) is real and the arctangent is positive, so a_{BCS}(η) is also negative. Hence a(η) < 0 throughout 0<η<1/3. With a<0, Eq. (5) makes repulsive g>0 grow rather than flow to zero, and makes attractive g<0 flow to zero rather than diverge, which is the opposite of the paper's stated conclusion. For η>1/3, the logarithm in Eq. (B1) has a negative argument and the square root in Eq. (B2) becomes imaginary, so the printed beta function is not real-valued there. The paper's claim that 'a(η) is positive for all η' is therefore not supported by the written calculation; the angular integrands in Eqs. (B8) and (B10) are positive where nonzero, so the minus signs in Eqs. (B1)-(B2) appear to be the point of failure. A sign typo is possible, but as printed the derivation inverts the central weak-coupling conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a single-flavor (spin- and valley-polarized) two-dimensional Fermi liquid with C3v symmetry tuned to a Van-Hove singularity, using a frequency-cutoff patch renormalization group. For three VH points it derives a one-loop beta function for the single inter-patch coupling and claims that repulsive interactions are marginally irrelevant, while attractive interactions produce either CDW or PDW order depending on mass anisotropy and patch orientation, with a 3q PDW supporting h/6e vortices. A six-VH-point extension and an RPA estimate of the initial coupling for rhombohedral tetralayer graphene are also presented. The paper is clearly written and the RG setup is well motivated, but the central one-loop coefficients as printed have the wrong sign in 0<η<1/3 and are non-real for η>1/3, so the main stability/instability claim is not supported by the calculation as written.","tokens_in":18849,"tokens_out":10411,"duration_ms":101939,"significance":"If the RG calculation is corrected and the claimed sign of a(η) holds, the paper would provide a rare weak-coupling mechanism for PDW and CDW order in a spin/valley-polarized system, directly relevant to recent experiments on rhombohedral multilayer graphene. The manuscript contains explicit derivations of the beta functions, susceptibility exponents, GL coefficients, and an RPA estimate of the initial g(0), which are valuable for reproducibility. The h/6e vortex result is a standard re-derivation and not new. At present, however, the sign/domain inconsistency in the central beta function prevents the main physical conclusion from being accepted.","major_comments":[{"comment":"The central claim that weak repulsive interactions are marginally irrelevant and weak attractive interactions drive CDW/PDW order rests entirely on a(η)>0 in Eq. (5). The printed coefficients in Eqs. (B1)–(B2) do not satisfy this. For 0<η<1/3, the logarithm in Eq. (B1) is positive (its argument exceeds 1) and the arctangent in Eq. (B2) is positive, while both prefactors carry an overall minus sign, so a_ZS'(η), a_BCS(η), and hence a(η) are negative. For 1/3<η<1, the logarithm in Eq. (B1) has a negative argument and the square root in Eq. (B2) is imaginary, so the beta function is not real. With a(η)<0, Eq. (5) makes positive g grow and negative g flow to zero, which is the opposite of the claimed behavior. Since Eqs. (8) and (9) inherit these signs, the predicted stability/instability dichotomy is not established by the manuscript as written. A similar domain problem appears in Eq. (E1), where the argument of arcsin can exceed unity for allowed values of η and φ, so the six-VH coefficient A(η,φ) is not real over the full domain of Fig. 6. The authors should recompute or correct these coefficients and verify the sign of a(η) on the entire interval 0<η<1.","section":"Eq. (5), Appendix B, and Appendix E"},{"comment":"The selection between single-q and multi-q order is not derived from the microscopic model. The quartic GL functional in Eq. (10) is unbounded below (Appendix D), and the authors stabilize it by adding a phenomenological sixth-order term w[Σ(|Δ|²+ζ|ρ|²)]³ with arbitrary positive coefficients. The 3q PDW region in Fig. 4 therefore depends both on this uncalculated sixth-order term and on the assumed T_PDW(η) and T_CDW(η) profiles stated in the caption; the computed GL coefficients are evaluated only at η=ηc. The one-loop RG instability calculation is not affected, but the 'multiple wavevector' part of the central claim is illustrative rather than controlled. I recommend either computing the sixth-order coefficients from the patch model or explicitly presenting the multi-q phase selection as model-dependent.","section":"Sec. II.D and Appendix D"}],"minor_comments":[{"comment":"There are several typos: 'This a consequence' in the Introduction, 'the the CDW phase' in the Introduction, and repeated misspellings of 'anisotropy' as 'anistropy'.","section":"Throughout"},{"comment":"The lower-left block of the 6x6 matrix in Eq. (D7) appears to be written incorrectly: rows 4–6 begin with cP, bP, bP rather than the corresponding u-type coefficients, so the matrix does not have the expected block structure. Please check this expression.","section":"Eq. (D7)"},{"comment":"The figure caption states that the phase diagram is independent of w as long as w>0, but the boundary between 3q CDW and 1q PDW depends on the phenomenological parameter ζ; this dependence should be stated in the main text as well.","section":"Fig. 4 caption"},{"comment":"Reference [48] is formatted as 'Jiang-Xiazi et al. Lin'; the author name should be cleaned up to a standard format.","section":"Reference [48]"}],"recommendation":"major_revision","confidential_remarks":"The sign/domain error in Appendix B is the main blocker. If the corrected calculation yields a(η)<0, the paper's principal conclusion would be reversed, so the revision must include a full recomputation of the one-loop coefficients. The six-VH formula in Appendix E should also be checked over the whole (η,φ) domain. The GL sixth-order stabilization is a separate concern but is less central to the rejection risk."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe stress-test note holds up. The one-loop coefficients in Eqs. (B1) and (B2) as printed give a(η) < 0 for 0<η<1/3 and non-real values for η>1/3. That is the opposite of what the text needs for the central claim, and it reverses the stability/instability dichotomy in Eq. (5). This is almost certainly a sign typo—one stray minus sign in each coefficient—but it is load-bearing. The paper as submitted cannot be accepted without the algebra being redone and the conclusion re-derived.\n\nThat said, the idea itself is worth taking seriously. The authors isolate a single-flavor (spin and valley polarized) C3v VH system, where both particle-particle and particle-hole susceptibilities diverge only logarithmically, so the weak-coupling RG is on better footing than the spinful VH studies that have to contend with log-squared terms. The three-patch and six-patch analyses are clear, the GL phase diagrams (though qualitative) organize the CDW/PDW competition nicely, and the RPA estimate for an attractive inter-patch interaction in R4G is a concrete attempt to connect to the recent experiments. The h/6e vortex discussion is a clean re-derivation of known PDW physics, correctly attributed.\n\nThe other soft spots are secondary. The GL free energy is unbounded at quartic order, so Fig. 4 relies on a phenomenological sixth-order term with undetermined ζ and assumed T(η) profiles—they are honest about this, but it limits the phase diagram's predictive power. And the fully polarized normal state is assumed, not derived. They say so in Sec. IV B, but it means the model applies only if the experiment truly is in a fully isospin-polarized regime. For R4G/R5G that is plausible but not established.\n\nWho gets value from this? People working on rhombohedral graphene superconductivity and on weak-coupling mechanisms for PDW. The paper deserves a serious referee—not a desk reject—because the core model is novel and the error looks fixable. I would not cite the central result until the sign is confirmed. If the authors confirm the flip, the qualitative conclusions likely survive.\n\nRecommendation: send it to review, with instructions to verify the beta function algebra and the sign of a(η) before judging the physics.","headline":"A promising weak-coupling PDW mechanism that is currently undermined by a sign error in the printed one-loop beta function—likely a typo, but load-bearing.","tokens_in":19360,"tokens_out":6366,"would_cite":false,"duration_ms":58965,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that a fully spin- and valley-polarized electron gas at a Van-Hove singularity develops pair density wave order when its effective interactions are attractive, while repulsive interactions leave it a stable metal.","keywords":["pair density wave","charge density wave","Van Hove singularity","spin and valley polarization","renormalization group","rhombohedral multilayer graphene","C3v symmetry","fractional vortices"],"falsifier":"Measure the spin and valley polarization of the normal state in rhombohedral multilayer graphene at the Van-Hove density, for example by quantum oscillation or compressibility probes: seeing Fermi surfaces from both valleys or both spin species would invalidate the single-flavor model's phase diagram. Alternatively, in a candidate 3q PDW state, magnetometry that counts only $h/2e$ vortices per applied flux quantum would contradict the predicted $h/6e$ vortices.","tokens_in":18274,"feed_emoji":"🌀","tokens_out":6985,"duration_ms":67972,"temperature":0.7,"pith_summary":"The paper studies a two-dimensional electron gas that is fully spin- and valley-polarized and tuned so that the Fermi level sits at a Van-Hove singularity, where the density of states diverges logarithmically. Its central claim is that the fate of the system is governed by a single effective interaction between symmetry-related saddle points, together with the mass anisotropy of the dispersion. For repulsive interpatch interactions the system remains a stable metal even at the singularity. For attractive interactions it develops either a charge density wave or a pair density wave, a superconducting state whose Cooper pairs carry finite momentum. The three-wavevector pair density wave supports vortices carrying one-sixth of the superconducting flux quantum, a concrete signature that experiments could look for.","feed_headline":"Attractive electrons at a Van-Hove point form pair density waves","feed_subtitle":"Weak-coupling RG shows CDW or PDW wins by saddle-point anisotropy; a three-wavevector PDW can host h/6e vortices.","key_machinery":"The central object is the Van-Hove patch model: a low-energy theory keeping only fermion fields near the three or six symmetry-related saddle points, where the density of states diverges logarithmically. The heavy lifting is done by a one-loop frequency-shell renormalization group, in which both the BCS (particle-particle) and ZS' (particle-hole) diagrams contribute to a flow coefficient $a(\\eta)$ or $A(\\eta,\\phi)$ that is positive everywhere and diverges at $\\eta=1/3$. The mass-anisotropy parameter $\\eta$ and the principal-axis angle $\\phi$ enter only through these flow coefficients, and the test-vertex flows then determine whether CDW or PDW is the leading instability.","core_discovery":"The paper's central claim is that in a fully spin- and valley-polarized $C_{3v}$-symmetric two-dimensional metal at a Van-Hove singularity, the weak-coupling physics is governed by the sign and geometry of the interpatch interaction. In the three-VH-point model the one-loop flow of the single coupling is $\\dot g=-a(\\eta)g^2$ with $a(\\eta)>0$, so repulsive $g$ flows to zero and the metal is stable, while attractive $g$ diverges at the energy scale $E_c=\\Lambda_0\\exp(-1/t_c)$. Tracking test vertices gives susceptibility exponents $\\alpha_{\\text{CDW}}=-a_{ZS'}(\\eta)/a(\\eta)$ and $\\alpha_{\\text{PDW}}=-2a_{BCS}(\\eta)/a(\\eta)$, producing a crossover at $\\eta_c\\approx 0.157$ between PDW-dominated ($\\eta<\\eta_c$) and CDW-dominated ($\\eta>\\eta_c$) regions. The six-VH-point model adds a second geometric parameter, the angle $\\phi$ between principal axes of mirror-related saddle points, and yields a phase diagram in $(\\eta,\\phi)$ in which CDW fills most of the plane and PDW appears at large mass anisotropy and small angles around $\\phi=n\\pi/3$. A three-component pair density wave, selected near the transition, breaks translation symmetry and supports $h/6e$ vortices bound to lattice dislocations.","pith_inferences":["If the polarized normal state is confirmed in rhombohedral multilayer graphene, this mechanism suggests the observed superconductivity is a PDW selected by the specific $\\eta$ and $\\phi$ of the band structure; tuning the displacement field across the VH point should move the system along the CDW/PDW boundary.","Because the RG treats particle-particle and particle-hole channels on equal footing here, the same single-flavor VH setup could serve as a controlled theoretical laboratory for competing PDW and CDW order in other fully valley-polarized moiré systems, not only rhombohedral graphene.","The predicted $h/6e$ vortices provide a direct test: counting vortices versus applied magnetic field in a SQUID scan, or searching for zero-field dislocation-bound fractional vortices, could distinguish a $3q$ PDW from a conventional $h/2e$ superconductor.","The RPA estimate that an effective interpatch attraction can emerge from repulsive Coulomb interactions depends on wavefunction overlap factors near the VH points; a more quantitative calculation of that overlap would sharpen the experimental regime where the predicted instability applies."],"forward_implications":["A fully polarized normal state with net attractive inter-VH interactions will not remain a Fermi liquid; it develops an ordered state at a parametrically low energy scale $E_c$.","Which order wins is controlled by saddle-point geometry: the mass-anisotropy ratio $\\eta$, and for six VH points also the angle $\\phi$, select between CDW and PDW, so tuning a displacement field or strain could switch the ground state.","The $3q$ PDW state is a translation-symmetry-breaking superconductor with an accompanying $3q$ CDW, and its elementary defects carry flux $h/6e$ rather than $h/2e$, giving a sharp experimental fingerprint.","In the weak-coupling limit, PDW and CDW orders leave gapless fermionic excitations; PDW states exhibit Bogoliubov Fermi surfaces, so the low-energy spectrum is not fully gapped.","Repulsive interactions are marginally irrelevant even at the diverging density of states, so no Stoner-type or nematic instability appears in this single-flavor model."],"supporting_citations":[{"why":"Provides the motivating experimental state: superconductivity in a fully spin- and valley-polarized rhombohedral multilayer graphene.","marker":"[27]"},{"why":"Supplies the R4G band structure used to locate three VH points and as input to the RPA estimate of the initial interaction.","marker":"[28]"},{"why":"Previous spinful RG analysis whose log-squared particle-hole divergences the single-flavor model avoids; serves as a comparison target for PDW at VH points.","marker":"[42]"},{"why":"Fixes the $\\eta=1/3$ honeycomb-tight-binding limit where the particle-hole channel acquires a log-squared divergence.","marker":"[40]"},{"why":"Provides the interacting-fermion renormalization-group framework and the BCS, ZS, and ZS' diagram classification used in the beta function.","marker":"[45]"},{"why":"Established the analysis of defects of three-component PDW states used to derive the $h/6e$ fractional vortices.","marker":"[49]"},{"why":"Gives the RPA screened-interaction formula used in Appendix A to estimate the sign of the effective interpatch coupling.","marker":"[53]"}],"fun_headline_variants":["Van-Hove attraction yields PDW with h/6e vortices","Spin-polarized 2DEG: PDW and CDW from Van-Hove geometry","Attractive Van-Hove electrons: from PDW to h/6e vortices","Van-Hove singularity drives pair and charge density waves","Fractional vortices from three-wavevector pair density wave"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole calculation assumes that a fully spin- and valley-polarized normal state already exists; if the actual normal state is not fully polarized, the single-flavor patch model and its predicted phase boundaries do not describe the system.","fun_headline_variants_meta":{"raw":{"variants":["Van-Hove attraction yields PDW with h/6e vortices","Spin-polarized 2DEG: PDW and CDW from Van-Hove geometry","Attractive Van-Hove electrons: from PDW to h/6e vortices","Van-Hove singularity drives pair and charge density waves","Fractional vortices from three-wavevector pair density wave"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000607,"raw_usage":{"total_tokens":2878,"prompt_tokens":1044,"completion_tokens":1834,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":1738}},"tokens_in":660,"tokens_out":1834,"duration_ms":12082,"temperature":1.0,"reasoning_tokens":1738,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:56:23.790645+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the spin and valley polarization of the normal state in rhombohedral multilayer graphene at the Van-Hove density, for example by quantum oscillation or compressibility probes: seeing Fermi surfaces from both valleys or both spin species would invalidate the single-flavor model's phase diagram. Alternatively, in a candidate 3q PDW state, magnetometry that counts only $h/2e$ vortices per applied flux quantum would contradict the predicted $h/6e$ vortices.","supporting_citations":[{"cited_title":"Pair-Density-Wave Superconductivity: A Microscopic Model on the 2D Honeycomb Lattice,","cited_arxiv_id":null,"evidence_quote":"Provides the motivating experimental state: superconductivity in a fully spin- and valley-polarized rhombohedral multilayer graphene."},{"cited_title":"Pair density wave and loop current promoted by Van Hove singularities in moir´ e systems,","cited_arxiv_id":null,"evidence_quote":"Supplies the R4G band structure used to locate three VH points and as input to the RPA estimate of the initial interaction."},{"cited_title":"Chiral superconductivity from repulsive interactions in doped graphene,","cited_arxiv_id":null,"evidence_quote":"Previous spinful RG analysis whose log-squared particle-hole divergences the single-flavor model avoids; serves as a comparison target for PDW at VH points."},{"cited_title":"Antifer- romagnetism and superconductivity in a quasi two- dimensional electron gas. Scaling theory of a generic Hubbard model,","cited_arxiv_id":null,"evidence_quote":"Fixes the $\\eta=1/3$ honeycomb-tight-binding limit where the particle-hole channel acquires a log-squared divergence."},{"cited_title":"Half- and quarter- metals in rhombohedral trilayer graphene,","cited_arxiv_id":null,"evidence_quote":"Provides the interacting-fermion renormalization-group framework and the BCS, ZS, and ZS' diagram classification used in the beta function."},{"cited_title":"d-Wave Superconductivity and Pomeranchuk Instability in the Two-Dimensional Hubbard Model,","cited_arxiv_id":null,"evidence_quote":"Established the analysis of defects of three-component PDW states used to derive the $h/6e$ fractional vortices."},{"cited_title":"Spectral signatures of modulated d-wave superconducting phases,","cited_arxiv_id":null,"evidence_quote":"Gives the RPA screened-interaction formula used in Appendix A to estimate the sign of the effective interpatch coupling."}],"review_version":1}