Resonating valence bond pairing energy in graphene by quantum Monte Carlo
Pith reviewed 2026-05-17 23:18 UTC · model grok-4.3
The pith
A finite gap near the Fermi level stabilizes resonating valence bond pairing in graphene with a binding energy of 0.48 millihartree per atom.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Variational and diffusion quantum Monte Carlo calculations demonstrate that the resonating valence bond state carries no binding energy in the zero-gap limit of rectangular graphene but acquires a finite pairing energy of approximately 0.48(1) millihartree per atom once a size-induced gap appears near the Fermi level, with the gap vanishing only when the x-length satisfies Lx = 3n sqrt(3) d.
What carries the argument
The geometry-dependent single-particle gap near the Fermi level that appears in rectangular samples lacking pi/3 rotational symmetry, which controls whether diffusion Monte Carlo finds a stable resonating valence bond paired state.
If this is right
- No stable resonating valence bond pairing occurs in gapless graphene under the conditions studied.
- A finite gap near the Fermi level produces a pairing energy of 0.48 millihartree per atom that survives extrapolation to the thermodynamic limit.
- The pairing is driven by the opening of a gap through choice of rectangular dimensions rather than by intrinsic interaction strength alone.
- The result applies specifically to confined graphene nanostructures whose boundaries break the usual rotational symmetry.
Where Pith is reading between the lines
- If the gap-stabilized pairing survives in other geometries, then deliberately patterning graphene ribbons or flakes to open a gap could induce similar pairing without doping.
- The same size-dependent gap mechanism may appear in related two-dimensional carbon systems such as graphene nanoribbons or bilayer structures under confinement.
- Testing whether the pairing energy remains finite when the gap is opened by strain or substrates rather than by finite-size effects would extend the result beyond the rectangular samples examined here.
Load-bearing premise
The chosen Jastrow-Slater-determinant and Jastrow-antisymmetrized-geminal-power wave functions must capture the essential correlations of the true resonating valence bond state, and finite-size scaling from rectangular samples must reach the thermodynamic limit without dominant boundary-condition biases.
What would settle it
A calculation or measurement on successively larger rectangular samples that finds the extrapolated pairing energy approaching zero even in the finite-gap configurations would falsify the claim of stable gap-stabilized pairing.
Figures
read the original abstract
We determine the resonating-valence-bond (RVB) state in graphene using real-space quantum Monte Carlo with correlated variational wave functions. Variational and diffusion quantum Monte Carlo (DMC) calculations with Jastrow-Slater-determinant and Jastrow-antisymmetrized-geminal-power ansatze are employed to evaluate the RVB pairing energy. Using a rectangular graphene sample that lacks $\pi/3$ rotational symmetry, we found that the single-particle energy gap near the Fermi level depends on the system size along the $x$-direction. The gap vanishes when the length satisfies $L_x=3n\sqrt{3}d$, where $n$ is an integer and $d$ is the carbon-carbon bond length, otherwise, the system, exhibits a finite gap. Our DMC results show no stable RVB pairing in the zero-gap case, whereas the opening of a finite gap near the Fermi level stabilizes the electron pairing. The DMC predicted absolute value of pairing energy at the thermodynamic limit for a finite-gap system is $\sim 0.48(1)$ mHa/atom. Our results reveal a feometry-driven electron pairing mechanism in the confined graphene nanostructure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses variational and diffusion Monte Carlo (DMC) with Jastrow-Slater-determinant and Jastrow-antisymmetrized-geminal-power ansatze to compute the resonating-valence-bond (RVB) pairing energy in rectangular graphene samples. It reports that the single-particle gap near the Fermi level vanishes for L_x = 3n√3 d but is finite otherwise; DMC finds zero pairing energy in the zero-gap cases but a finite extrapolated value of ~0.48(1) mHa/atom in the thermodynamic limit for the gapped sequence, interpreted as evidence for a geometry-driven pairing mechanism.
Significance. If the central numerical results hold after addressing finite-size consistency, the work supplies a concrete DMC benchmark for gap-stabilized RVB pairing in graphene nanostructures, including explicit statistical error bars and a size-dependent gap condition that is independent of the final energy value. The dual use of JSD and JAGP wave functions and the focus on rectangular geometries without π/3 symmetry constitute clear methodological strengths.
major comments (2)
- [Abstract and §4] Abstract and §4 (Results): the extrapolation of pairing energy to the thermodynamic limit for the finite-gap sequence yields a non-zero value (~0.48(1) mHa/atom), yet the manuscript states that the gap magnitude depends on L_x and vanishes for specific sizes; because the gapped sequence has gap ~1/L_x (or faster) that closes in the L_x → ∞ limit, the extrapolated pairing energy should approach the zero value reported for the gapless sequence unless an additional mechanism keeps an effective gap open under the chosen boundary conditions.
- [§3 and §5] §3 (Methods) and §5 (Discussion): the claim that the Jastrow-AGP ansatz captures the true RVB state sufficiently to demonstrate gap stabilization rests on the variational parameters and rectangular geometry; however, the geminal structure of the AGP may preferentially stabilize pairing precisely when a single-particle gap is artificially imposed by the boundary conditions, and no explicit test (e.g., comparison of gap closure rate versus pairing energy across multiple L_x in the gapped family) is provided to rule out this bias.
minor comments (2)
- [Figure 2] Figure 2 caption: the labeling of zero-gap versus finite-gap data points should explicitly reference the L_x = 3n√3 d condition to aid readability.
- The manuscript would benefit from a brief comparison to prior variational Monte Carlo studies of RVB states in graphene (e.g., those using different boundary conditions) to contextualize the reported pairing energy scale.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment point by point below, indicating where revisions will be made to improve clarity and address the concerns raised.
read point-by-point responses
-
Referee: [Abstract and §4] Abstract and §4 (Results): the extrapolation of pairing energy to the thermodynamic limit for the finite-gap sequence yields a non-zero value (~0.48(1) mHa/atom), yet the manuscript states that the gap magnitude depends on L_x and vanishes for specific sizes; because the gapped sequence has gap ~1/L_x (or faster) that closes in the L_x → ∞ limit, the extrapolated pairing energy should approach the zero value reported for the gapless sequence unless an additional mechanism keeps an effective gap open under the chosen boundary conditions.
Authors: We agree that the single-particle gap for the gapped sequence scales to zero in the L_x → ∞ limit, as stated in the manuscript. Our extrapolation of the pairing energy is performed separately on the sequence of finite-size systems that exhibit a non-zero gap due to the rectangular boundary conditions. The resulting non-zero value of 0.48(1) mHa/atom reflects the pairing energy associated with this geometry-driven gap stabilization at finite sizes. In the true thermodynamic limit under periodic boundary conditions the gap vanishes and pairing is absent, consistent with the zero result for the gapless sequence. We do not invoke an additional mechanism to keep a gap open at infinite size. We have revised §4 to explicitly discuss the gap closure rate, the separate extrapolations for each sequence, and the physical interpretation of the finite extrapolated value as a benchmark for gap-stabilized RVB pairing. revision: yes
-
Referee: [§3 and §5] §3 (Methods) and §5 (Discussion): the claim that the Jastrow-AGP ansatz captures the true RVB state sufficiently to demonstrate gap stabilization rests on the variational parameters and rectangular geometry; however, the geminal structure of the AGP may preferentially stabilize pairing precisely when a single-particle gap is artificially imposed by the boundary conditions, and no explicit test (e.g., comparison of gap closure rate versus pairing energy across multiple L_x in the gapped family) is provided to rule out this bias.
Authors: The referee correctly identifies a potential source of bias in the AGP ansatz. We have therefore performed calculations with both the Jastrow-Slater-determinant and Jastrow-antisymmetrized-geminal-power wave functions. The DMC results are consistent between the two ansatze, with pairing energy vanishing in the gapless sequence and remaining finite in the gapped sequence. This agreement indicates that the stabilization is driven by the presence of the gap rather than by a specific feature of the geminal structure. The rectangular geometry without π/3 symmetry was deliberately chosen to allow the gap to be tuned by system size. While we did not originally include a direct comparison of gap magnitude versus pairing energy across the gapped family, we will add this analysis together with the relevant variational parameters in the revised §5 to further rule out ansatz bias. revision: yes
Circularity Check
No significant circularity; pairing energy from explicit DMC projection
full rationale
The derivation computes RVB pairing energy via variational and diffusion Monte Carlo on Jastrow-Slater and Jastrow-AGP wave functions for rectangular graphene supercells whose single-particle gap is determined by geometry (Lx = 3n√3 d). Zero-gap sizes yield zero pairing energy while gapped sizes yield finite values; the thermodynamic-limit value is obtained by extrapolation of these independent numerical results. No step defines the output energy in terms of itself, renames a fitted parameter as a prediction, or relies on a self-citation chain that is itself unverified. The ansatz and boundary conditions are chosen independently of the final energy number, satisfying the criteria for a self-contained computation.
Axiom & Free-Parameter Ledger
free parameters (1)
- Jastrow and geminal variational parameters
axioms (1)
- domain assumption The RVB state is adequately represented by Jastrow-Slater-determinant and Jastrow-antisymmetrized-geminal-power forms
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
DMC predicted absolute value of pairing energy at the thermodynamic limit for a finite-gap system is ∼0.48(1) mHa/atom... gap vanishes when L_x = 3n√3 d
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Jastrow-antisymmetrized-geminal-power ansatz... pairing energy δ_P
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
An arXiv:2511.06506v1 [cond-mat.str-el] 9 Nov 2025 2 FIG. 1. A spin-singlet dimer configuration (blue color) on graphene lattice. The resonating valence bond (RVB) state is composed of a superposition of such configurations. Threesp 2 orbitals form theσband containing three localized electrons, while the bonding among thep z orbitals of different lattice ...
work page internal anchor Pith review Pith/arXiv arXiv 2025
-
[2]
K. Novoselov, A. Geim, S. Morozov, D. Jiang, Y. Zhang, S. Dubonos, I. Grigorieva, and A. Firsov, Electric field effect in atomically thin carbon films, Science306, 666 (2004)
work page 2004
-
[3]
A. Geim and K. Novoselov, The rise of graphene, Nat. Mterials6, 183 (2007)
work page 2007
-
[4]
A. C. Neto, F. Guinea, N. M. Peres, K. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys.81, 109 (2009)
work page 2009
-
[5]
K. Novoselov, A. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. Grigorieva, S. Dubonos, and A. Firsov, Two-dimensional gas of massless dirac fermions in graphene, Nature438, 197 (2005)
work page 2005
-
[6]
V. N. Kotov, B. Uchoa, V. M. Pereira, F. Guinea, and A. H. C. Neto, Electron-electron interactions in graphene: Current status and perspectives, Rev. Mod. Phys.84, 1067 (2012)
work page 2012
-
[7]
S. D. Sarma, S. Adam, E. Hwang, and E. Rossi, Elec- tronic transport in two-dimensional graphene, Rev. Mod. Phys.83, 407 (2011)
work page 2011
-
[8]
Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, Unconventional super- conductivity in magic-angle graphene superlattices, Na- ture566, 43 (2018)
work page 2018
- [9]
- [10]
-
[11]
K. Novoselov, V. Falko, L. Colombo, P. Gellert, M. Schwab, and K. Kim, A roadmap for graphene, Na- ture490, 192 (2012)
work page 2012
- [12]
-
[13]
K. Novoselov, Z. Jiang, Y. Zhang, S. Morozov, H. Stormer, U. Zeitler, J. Maan, G. Boebinger, P. Kim, and A. Geim, Room-temperature quantum hall effect in graphene, Science3151379(2007)
work page 2007
-
[14]
J. Drut and T. L¨ ahde, Is graphene in vacuum an insula- tor?, Phys. Rev. Lett.102, 026802 (2009)
work page 2009
-
[15]
E. Rotenberg, A. Bostwick, T. Ohta, J. L. McChesney, T. Seyller, and K. Horn, Origin of the energy bandgap in epitaxial graphene, Nat. Materials7, 258 (2008)
work page 2008
-
[16]
S. Sorella and E. Tosatti, Semi-metal-insulator transition of the hubbard model in the honeycomb lattice, Euro- phys. Lett.19, 699 (1992)
work page 1992
-
[17]
Sachdev,Low Dimensional Quantum Field Theories for Condensed Matter Physicists, edited by L
S. Sachdev,Low Dimensional Quantum Field Theories for Condensed Matter Physicists, edited by L. Yu, S. Lundqvist, and G. Morandi(World Scientific, Singapore, 1995)
work page 1995
-
[18]
Z. Meng, T. Lang, S. Wessel, F. Assaad, and A. Mura- matsu, Quantum spin liquid emerging in two-dimensional correlated dirac fermions, Nature464, 847 (2010)
work page 2010
-
[19]
M.-T. Tran and K.-S. Kim, Spin liquids in graphene, Phys. Rev. B83, 2011 (2011)
work page 2011
-
[20]
P. W. Anderson, The resonating valence bond state in La 2CuO4 and superconductivity, Science235, 1196 (1987)
work page 1987
-
[21]
G. Baskaran and P. W. Anderson, Gauge theory of high-temperature superconductors and strongly corre- lated fermi systems, Phys. Rev. B37, 580 (1988)
work page 1988
-
[22]
Pauling,The Nature of the Chemical Bond(Cornell University Press, Ithaca, NY, 1960)
L. Pauling,The Nature of the Chemical Bond(Cornell University Press, Ithaca, NY, 1960)
work page 1960
-
[23]
C. Glittum, A. ˇStrkalj, D. Prabhakaran, P. Goddard, C. D. Batista, and C. Castelnovo, A resonant valence bond spin liquid in the dilute limit of doped frustrated mott insulators, Nat. Phys.21, 1211 (2024)
work page 2024
-
[24]
P. W. Anderson, G. Baskaran, Z. Zou, and T. Hsu, Resonating–valence-bond theory of phase transitions and superconductivity in La 2CuO4-based compounds, Phys. Rev. Lett.58, 2790 (1987)
work page 1987
-
[25]
Y. Zhou, K. Kanoda, and T.-K. Ng, Quantum spin liquid states, Rev. Mod. Phys.89, 025003 (2017)
work page 2017
-
[26]
S. Pathak and V. B. G. Baskaran, Possible high- temperature superconducting state with a d+id pairing symmetry in doped graphene, Phys. Rev. B81, 085431 (2010)
work page 2010
-
[27]
Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and G. Saito, Spin liquid state in an organic mott insulator with a triangular lattice, Phys. Rev. Lett.91, 107001 (2003)
work page 2003
-
[28]
M. R. Norman, Colloquium: Herbertsmithite and the search for the quantum spin liquid, Rev. Mod. Phys.88, 041002 (2016)
work page 2016
-
[29]
M. Casula and S. Sorella, Improper s-wave symmetry of the electronic pairing in iron-based superconductors by first-principles calculations, Phys. Rev. B88, 155125 (2013)
work page 2013
-
[30]
F. Becca and S. Sorella,Quantum Monte Carlo ap- proches for correlated systems(Cambridge University, Cambridge, UK, 2017)
work page 2017
-
[31]
S. Sorella, N. Devaux, M. Dagrada, G. Mazzola, and M. Casula, Geminal embedding scheme for optimal atomic basis set construction in correlated calculations, J. Chem. Phys.143, 244112 (2015)
work page 2015
- [32]
- [33]
- [34]
-
[35]
M. Casula and S. Sorella, Geminal wave functions with jastrow correlation: A first application to atoms, J. Chem. Phys.119, 6500 (2003)
work page 2003
- [36]
-
[37]
A. Annaberdiyev, G. Wang, C. A. Melton, M. C. Ben- nett, L. Shulenburger, and L. Mitas, A new generation of effective core potentials from correlated calculations: 3d transition metal series, J. Chem. Phys.149, 134108 (2018)
work page 2018
-
[38]
M. C. Bennett, C. A. Melton, A. Annaberdiyev, G. Wang, L. Shulenburger, and L. Mitas, A new genera- tion of effective core potentials for correlated calculations, J. Chem. Phys.147, 224106 (2017)
work page 2017
- [39]
-
[40]
J. P. Perdew and A. Zunger, Self-interaction correction to density-functional approximations for many-electron systems, Phys. Rev. B23, 5048 (1981)
work page 1981
-
[41]
S. Fahy, X. W. Wang, and S. G. Louie, Variational quantum monte carlo nonlocal pseudopotential approach to solids: Formulation and application to diamond, graphite, and silicon, Phys. Rev. B42, 3503 (1990)
work page 1990
-
[42]
C. J. Umrigar, J. Toulouse, C. Filippi, S. Sorella, and R. G. Hennig, Alleviation of the fermion-sign problem by optimization of many-body wave functions, Phys. Rev. Lett.98, 110201 (2007)
work page 2007
-
[43]
Sorella, Green function monte carlo with stochastic reconfiguration, Phys
S. Sorella, Green function monte carlo with stochastic reconfiguration, Phys. Rev. Lett.80, 4558 (1998)
work page 1998
-
[44]
W. M. C. Foulkes, L. Mitas, R. J. Needs, and G. Ra- jagopal, Quantum monte carlo simulations of solids, Rev. Mod. Phys.73, 33 (2001)
work page 2001
-
[45]
D. M. Ceperley, Ground state of the fermion one- component plasma: A monte carlo study in two and three dimensions, Phys. Rev. B18, 3126 (1978)
work page 1978
-
[46]
D. M. Ceperley and B. J. Alder, Ground state of the electron gas by a stochastic method, Phys. Rev. Lett. 45, 566 (1980)
work page 1980
-
[47]
C. Melton and L. Mitas, Quantum monte carlo with vari- able spins: Fixed-phase and fixed-node approximations, Phys. Rev. E96, 043305 (2017)
work page 2017
- [48]
- [49]
-
[50]
See Supplemental Material for more details of our calcu- lations and additional figures and data
-
[51]
A. Black-Schaffer and S. Doniach, Resonating va- lence bonds and mean-field d-wave superconductivity in graphite, Phys. Rev. B75, 134512 (2007)
work page 2007
-
[52]
Mitas, Quantum monte carlo pair orbital wave func- tions for periodic systems, Phys
L. Mitas, Quantum monte carlo pair orbital wave func- tions for periodic systems, Phys. Rev. Lett.135, 126401 (2025)
work page 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.