Pith. sign in

REVIEW 4 major objections 6 minor 57 references

Variational Monte Carlo Optimization of Topological Chiral Superconductors

T0 review · 4 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Pure Coulomb repulsion can stabilize topological chiral superconducting states without a Fermi-surface pairing instability.

desk verdict A credible, well-benchmarked VMC comparison that overreaches from 'beats the QFL' to 'relevant ground state above the Wigner crystal.' read the letter →

arxiv 2507.18582 v3 pith:XGKHI4L5 submitted 2025-07-24 cond-mat.str-el cond-mat.supr-con

Minho Luke Kim , Abigail Timmel , Xiao-Gang Wen This is my paper · ORCID
classification cond-mat.str-elcond-mat.supr-con
keywords variationalMonteCarlochiralsuperconductivitytopologicalorderPfaffianwavefunctionK-matrixLaughlinstatesquarterFermiliquidrhombohedralgrapheneCoulombinteraction
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Using variational Monte Carlo, this paper compares the ground-state energies of three topological chiral superconducting trial states (the single-species Pfaffian and the two-species K2a and K2b Laughlin-type states) against a fully optimized Slater-Jastrow description of the spin-valley polarized quarter Fermi liquid. For the dispersion $E_k = c_2 k^2 + c_4 k^4$, tuned to rhombohedral graphene parameters, it finds that the Pfaffian and K2a states have lower energy per electron than the Fermi liquid at densities up to about $0.5\times10^{12}\,\mathrm{cm}^{-2}$, with a condensation energy near 5% of the Coulomb energy scale, roughly 1 meV per electron. The preference is strongest when $c_2$ is between zero and a negative value, i.e. just before a hole pocket forms at $k=0$. The paper reads this as evidence that superconductivity can emerge from purely repulsive Coulomb interactions in a nearly flat band, without a Fermi-surface pairing instability.

What carries the argument

The carrying object is the family of generalized Laughlin-type chiral superconducting wavefunctions built from holomorphic and antiholomorphic factors, with a K-matrix $K = K^+ - K^-$ whose diagonal entries are odd integers and whose null vector fixes the species density ratios through $\sum_J K_{IJ} f_J = 0$, cancelling the macroscopic angular momentum that would otherwise make the kinetic energy diverge. For the one-species case the paper uses an improved Pfaffian ansatz $\psi_{\mathrm{Pf}} = \mathrm{Pf}\left(\frac{1}{z_i-z_j}\right)$ times a rotationally symmetric two-body product with four variational parameters, replacing the Gaussian decay of the Laughlin-type factor so that algebraic long-range density correlations are allowed. The comparison is made through the energy formula $E_{\mathrm{tot}} = (e^2\sqrt{n_e}/\epsilon) V + c_2\langle k^2\rangle + c_4\langle k^4\rangle$, where $V$ is obtained from Monte Carlo pair-distribution functions and the kinetic averages from numerical derivatives of the wavefunction.

What would settle it

Compute the Wigner crystal energy for the same parameters and densities; if it lies more than roughly 1 meV per electron below the optimized Pfaffian and K2a energies, the claimed ordering fails. Alternatively, run the same variational comparison with wavefunctions that include the compressible superfluid density fluctuations and Berry curvature; a shift of more than the roughly 1 meV condensation energy in either state's favor would overturn the conclusion.

Watch

Extended reading notes

Core claim

The central claim is that strong repulsive Coulomb interactions alone can stabilize topological chiral superconducting phases over the spin-valley polarized quarter Fermi liquid in the density window relevant to the rhombohedral graphene experiments. After optimizing the variational parameters in the Pfaffian and two-species K-matrix wavefunctions, and in the Slater-Jastrow Fermi liquid, the Pfaffian and K2a states win with a condensation energy of about 1 meV per electron. The win survives up to densities around $0.5\times10^{12}\,\mathrm{cm}^{-2}$, and it is largest for $c_2$ slightly negative, corresponding to a Fermi sea on the verge of developing a hole pocket at $k=0$. The paper argues these states are not BCS superconductors: they are driven by flux attachment and have a gapless superfluid density mode, short coherence length, broken time-reversal symmetry, and nontrivial topological order.

Load-bearing premise

The argument stands on the assumption, flagged in Section I, that the trial wavefunction family—the K-matrix Laughlin-type chiral states and the Pfaffian, plus the Slater-Jastrow Fermi liquid—contains the true ground-state correlations of the repulsive Coulomb Hamiltonian, so that the winner of the variational competition is the real phase. The paper itself notes in Section I that the superconducting trial states do not incorporate density fluctuations of the compressible superfluid mode, omit Berry curvature, and do not compute the Wigner crystal.

Editorial extensions

If this is right

  • If the variational ordering is correct, the observed chiral superconductivity in rhombohedral multilayers should be understood as a strongly correlated topological superconductor rather than a BCS condensate of a Fermi surface.
  • A spin-unpolarized two-species state (K2a, in the same phase as a spin-triplet p+ip superconductor) remains competitive with the spin-polarized Pfaffian, so in-plane magnetic field robustness does not by itself rule out spin-unpolarized superconductivity.
  • The superconducting window is tied to a nearly flat band bottom: the energy advantage is largest when $c_2$ is slightly negative, just before a hole pocket forms, and shrinks as $c_2$ becomes positive.
  • The same pure-repulsion mechanism should be sought in other two-dimensional systems with an almost flat band bottom and strong Coulomb interactions, not only rhombohedral graphene.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable extension is to track superconductivity against displacement field: this mechanism predicts the strongest chiral superconducting order immediately below the Lifshitz transition where $c_2 = -4\pi n_e c_4$, and a rapid suppression once $c_2$ is clearly positive.
  • The paper did not include Berry curvature; since it expects that effect to lower the chiral states more than the Fermi liquid, including it would likely widen the superconducting region, with the K2b state benefiting most.
  • The trial chiral states omit density fluctuations of the compressible superfluid mode; a family that includes those fluctuations could lower the superconducting energies further, and comparing its effect state-by-state is a direct way to test whether the variational ordering is stable.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. This manuscript uses variational Monte Carlo (VMC) to compare the energies of topological chiral superconducting trial states—a single-species Pfaffian state and two-species K2a/K2b Laughlin-type states—against an optimized Slater-Jastrow quarter Fermi liquid (QFL), for the dispersion E_k = c2 k^2 + c4 k^4 motivated by rhombohedral graphene. The authors introduce a modified Pfaffian ansatz (Eq. 7) and a two-parameter Slater-Jastrow QFL wavefunction (Eq. 12), optimize the variational parameters, and construct phase diagrams in density versus -c2 and in density versus magnetic field. They conclude that the Pfaffian and K2a chiral superconducting states can be energetically favored over the QFL at densities relevant to experiment (0.2–0.7 x 10^12 cm^-2), with a condensation energy of about 1 meV per electron, and they propose this as a mechanism for superconductivity from pure repulsive Coulomb interactions without a Fermi-surface pairing instability.

Significance. If the result is correct, it provides numerical support for a qualitatively new route to superconductivity in flat-band two-dimensional systems, driven by Coulomb repulsion through flux attachment rather than through a weak-pairing BCS instability, and it connects directly to the chiral superconducting phase observed in rhombohedral multilayer graphene. The paper's internal checks are a real strength: the QFL energies reproduce the Tanatar-Ceperley quadratic-dispersion result within 0.2% of the Coulomb energy, the alternative kinetic-energy estimators agree within 0.5%, the Hartree-Fock limit is recovered, and the new Pfaffian ansatz is shown to beat the previously used Laughlin-type ansatz. The important limitations—missing Wigner crystal energetics, missing error bars, and incompletely optimized chiral wavefunctions—are acknowledged by the authors, but they are precisely the quantities needed to make the phase-diagram claim quantitative rather than conditional.

major comments (4)
  1. [Section VI and Fig. 1] The central claim in the abstract and conclusion that chiral superconductivity wins 'above the density of Wigner crystal phase' is not established by the calculations, because the Wigner crystal is never computed from the same Hamiltonian. Section VI states that 'including the energetics of the Wigner crystal also will be crucial,' and the Wigner crystal boundary in Fig. 1 is taken from experiment rather than from a variational energy of the model. Given that the reported SC-QFL margins are only about 1 meV, a Wigner crystal trial state could interpose or shift the phase boundary without contradicting any calculation in the paper. The phase diagram should either include a Wigner crystal energy calculation or the claims should be restricted to a comparison between the chiral SC states and the QFL.
  2. [Fig. 4 and Section V] No statistical or finite-size error bars are reported for the central energy differences shown in Fig. 4. The text describes substantial Monte Carlo samples (for example, 5 x 10^5 samples for the Pfaffian kinetic energy with N=70 and 3 x 10^6 samples for the two-species potential energy with N=200), but the condensation energy of approximately 1 meV is presented without an uncertainty estimate and without a thermodynamic-limit extrapolation for the chiral droplet states. Since the phase boundaries in Fig. 1 and the magnetic-field diagram in Fig. 5 are separated by energy differences of this size, the statistical significance of the ordering cannot be assessed from the data as presented.
  3. [Section I and Eq. (7)] The authors explicitly state that the optimized Pfaffian, K2a, and K2b wavefunctions 'do not incorporate the density fluctuations of the compressible superfluid mode' and that Berry curvature is omitted. These omissions are expected to lower the energies of the chiral states, but because the trial states are variational, the comparison to the QFL is not a complete energy competition. The likely improvement of the chiral states relative to the QFL is plausible but unquantified, and it could be comparable to the 1 meV differences reported in Fig. 4. A quantitative estimate of the missing chiral-state correlations, or a clear statement that the phase boundaries may shift when they are included, is needed.
  4. [Appendix A and Fig. 7] The QFL side of the comparison also has unquantified variational bias. Appendix A reports that the trial Jastrow factor differs from the Gaskell form by less than 10% of the correlation energy; near the transition density the correlation energy is several meV, so this residual can be a substantial fraction of the claimed 1 meV SC condensation energy. The Gaskell comparison is made only for the quadratic-dispersion parameter set c2 = 91 meV nm^2, whereas the phase boundary in Fig. 4 uses c2 = -4 pi c4 n_e. The conclusion that an 'ideal' QFL would still lie above the chiral states therefore needs a quantitative propagation of the 10% residual into the phase diagram, not just the statement that the residual is less than 10%.
minor comments (6)
  1. [Eq. (10)] There is a typo in the displayed K matrix: the lower-right entry reads 'm.' and should be 'm'.
  2. [Section II.A] The parameters labeled 'xi_1' and 'xi_2' should be typeset as xi_1 and xi_2, and the sampling ranges written as '0.2 intervals' do not specify whether the grid is inclusive of the upper bound.
  3. [Appendix B.3] In the first sentence, 'naley' should be 'namely'.
  4. [Section VI] In the paragraph discussing future directions, 'supercondcutor' should be 'superconductor'.
  5. [Section V] The sentence beginning 'The indeed shows that as a positive quadratic term...' has an incomplete subject and should be rephrased.
  6. [Appendix C] The allowed occupation numbers n_k < 2 for two-species states and n_k < 1 for single-species states need a brief justification in terms of species degeneracy and spin/valley labels; as written, the reader may misread these bounds as a violation of Pauli exclusion.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the chiral-SC vs QFL energy competition is a self-contained variational Monte Carlo calculation, benchmarked externally; the self-citations to prior chiral-SC theory are not load-bearing.

full rationale

The energy-ordering claim does not reduce to its inputs. The chiral superconducting energies are obtained by direct Monte Carlo evaluation of explicit trial wavefunctions (Eqs. 1 and 7) with independently optimized variational parameters, and the QFL energy is computed from a separate Slater-Jastrow ansatz and validated against the known quadratic-dispersion result ('The two results agree within 0.2% of the Coulomb energy'). Neither the SC nor the QFL parameters are fitted to the final phase boundary, so the roughly 1 meV condensation energy is a genuine variational outcome rather than a relabeled fit. The identification of the trial states as chiral superconductors is imported from Ref. 17 and the anyon-superfluid literature (Refs. 29, 40), which are authored or co-authored by X.-G. Wen; however, these citations are parameter-free theoretical classifications that do not contain the present energetics result and are externally checkable via edge-mode and Hall-response predictions. The acknowledged omissions (Wigner crystal energetics, Berry curvature, hole pockets, compressible-mode density fluctuations) are explicitly flagged as future work, so they are incompleteness rather than circularity. No equation in the paper is equivalent by construction to its own input, and no fitted parameter is renamed as a prediction; hence no circular step is exhibited. The score of 2 reflects only the presence of minor self-citations, not a circular derivation.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central result rests on the credibility of the trial wavefunction family from the same group's prior work, on several modeling simplifications for rhombohedral graphene, and on accepting finite-size VMC energies as thermodynamically representative. The variational parameters listed above are optimized, not fitted to experimental data, but they are numerous relative to the roughly 1 meV energy differences being resolved. No new physical entities are introduced.

free parameters (5)
  • QFL Jastrow amplitude A and range B = A scanned 0.5 to 15; B scanned 0.1 to 1.0
    Two-parameter optimization of the Slater-Jastrow correlation; chosen by energy minimization at each density, not from experimental data.
  • Pfaffian ansatz parameters m, p, xi1, xi2 = m in 2..5, p in 0..5-m, xi1 <= xi2 in 0..1.6, step 0.2
    Variational parameters of the modified Pfaffian trial wavefunction; optimized by VMC energy comparison.
  • Kbar matrix elements a, b, c for two-species states = 5 >= a = b >= c >= 0, 0.2 or 0.1 intervals
    Continuous deformations of the K-matrix wavefunctions K2a and K2b; optimized for kinetic and potential energy separately.
  • Edge decay parameters sigma, s = sigma = 3, s = 2.5
    Chosen by hand for the droplet decay function; authors state they do not affect bulk physics.
  • K-matrix integer m for two-species states = m = 1 for K2a and m = 3 for K2b
    Discrete choice of topological states from Ref. 17; higher odd m states were not tested because authors judge them less likely experimentally.
assumptions (6)
  • domain assumption The generalized Laughlin-type wavefunction (1) with a K matrix satisfying Eqs. (2) and (3) is a valid chiral superconducting state of the many-electron system.
    Inherited from Ref. 17; this paper optimizes parameters inside this family but does not re-derive the topological order or the superfluid interpretation.
  • domain assumption A gapless charged density mode implies superconductivity, and the K-matrix states host such a mode because cofluctuations of the species densities are unconstrained.
    Section II invokes Refs. 17, 29, and 40 for this chain; it connects the trial wavefunctions to the claim of superconductivity without a BCS pairing instability.
  • domain assumption The rhombohedral graphene superconductor is captured by one or two electron species with dispersion E_k = c2 k^2 + c4 k^4, with c4 fixed to 549 or 366 meV nm^4 and c2 tunable, with no Berry curvature or screening corrections.
    Section V sets these parameters from Ref. 12; the authors explicitly leave Berry curvature effects for future work in Sections I and VI.
  • domain assumption The quarter Fermi liquid (spin-valley polarized Fermi liquid) is the correct competing phase in the density window, and the Wigner crystal can be positioned by experimental density rather than by computation.
    The paper computes only QFL versus chiral states; Section VI says the Wigner crystal energetics are not determined here, yet the abstract places the result above the Wigner crystal density.
  • domain assumption Variational Monte Carlo estimates on finite droplets (200 electrons for two-species states, 70 for the Pfaffian) and on a torus (69 electrons for the QFL), with restricted bulk sampling, represent thermodynamic-limit bulk energies.
    Finite-size and geometry effects are controlled only by consistency checks, such as the 0.2% quadratic QFL benchmark and 0.5% kinetic cross-checks, not by systematic extrapolation for the chiral states.
  • ad hoc to paper The Jastrow factor ansatz f(r) = A exp(-r/B)(1 + r/B + r^2/(2B^2)) with 9 supercell images is sufficiently expressive for the quarter Fermi liquid.
    Introduced in Appendix A; benchmarked against the Gaskell correlation energy within 10% and against Tanatar-Ceperley within 0.2% for the quadratic case, but tested only for this ansatz family.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Variational Monte Carlo Optimization of Topological Chiral Superconductors." pith.science (2026). https://pith.science/paper/XGKHI4L5

@misc{pith2026250718582,
  author       = {Pith},
  title        = {Pith review of: Variational Monte Carlo Optimization of Topological Chiral Superconductors},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/XGKHI4L5}},
  note         = {Machine review of arXiv:2507.18582}
}
abstract

We perform the variational Monte Carlo calculation for recently proposed chiral superconducting states driven by strong Coulomb interactions. We compare the resulting energetics of these electronic phases for the electron dispersion relation $E_k = c_2 k^2+c_4 k^4$. Motivated by the recent discovery of chiral superconductivity in rhombohedral graphene systems, we apply our analysis to relevant parameter regimes. We demonstrate that topological chiral superconducting phases (including a spin-unpolarized state) can be energetically favored over the spin-valley polarized Fermi liquid above the density of Wigner crystal phase. Our results show that the preference for chiral superconductivity is strongest when $c_2$ lies between zero and a negative value corresponding to a Fermi sea on the verge of forming a hole pocket around $k=0$. This finding suggests that superconductivity can arise from pure repulsive Coulomb interactions in systems with an almost flat band bottom, without relying on the pairing instability of a Fermi surface. This mechanism opens a new pathway to superconductivity beyond the conventional BCS mechanism.

Figures

Figures reproduced from arXiv: 2507.18582 by the authors.

Figure 1
Figure 1. FIG. 1: The phase diagram for electron dispersion [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2: Electron dispersion [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3: Comparison of two ground state energies (per [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: FIG. 4: Ground state energy per electron for [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5: The phase diagram for electron dispersion [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6: The phase diagram proposed in Ref. [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7: Comparison of the correlation energy for our [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8: A model distribution for fermion occupation [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]

Review history (2 revisions) →

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

57 extracted references · 29 canonical work pages

  1. [1]

    If we add a negative quadratic term (c 2 <0), the chiral superconducting states become more favorable

    The indeed shows that as a positive quadratic term (c2 >0) is added, the chiral superconducting states be- come gradually unfavorable, appearing at lower electron densities. If we add a negative quadratic term (c 2 <0), the chiral superconducting states become more favorable. However,c 2 cannot be too negative, since a very nega- tivec 2 will cause a Ferm...

    work page

  2. [2]

    The Potential Energy 10

    work page

  3. [3]

    A model of electron distribution ink-space 12 References 12 I

    Alternative Method for Kinetic Energy Sampling 11 C. A model of electron distribution ink-space 12 References 12 I. INTRODUCTION Tunable 2D electron systems have been demonstrated to be rich systems harboring various phases of matter derived from strong correlations and topology. One of the most striking features of physics discovered in these systems is ...

    work page

  4. [4]

    For ring-like Fermi surface above the dashed-line, the superconductor is non-topological with no chiral edge state. For disk-like Fermi surface below the dashed-line, the superconductor is topological with chiral edge state of a single Majorana fermion, as well as a Majorana zero mode in a magnetic vertex. Wigner crystal appears below electron density∼0.2...

    work page

  5. [5]

    As we did for the quarter Fermi liquid in the previous sec- tion, we can sample⟨k 2⟩and⟨k 4⟩from sampling the local energy

    Kinetic Energy The kinetic energy for the two-species chiral supercon- ductor can be directly sampled using the local energy. As we did for the quarter Fermi liquid in the previous sec- tion, we can sample⟨k 2⟩and⟨k 4⟩from sampling the local energy. For⟨k 2⟩, we use− ∇2ψ ψ . While we are running our Monte Carlo simulations, the integration is still over t...

    work page

  6. [6]

    The Potential Energy The interaction energy between a species-I 0 particle and a species-J 0 particle is given as: UI0J0 ∫︂ ∏︂ I,i d2zI iV(z I0−z J0)|Ψ(zI0,zJ0,{zI i})|2.(B1) The pair correlation function of the two electrons is de- fined as gI0J0(zI0−z J0) (πR2)2 = ∫︂ ∏︂ I,i d2zI i|Ψ(zI0,zJ0,{zI i})|2 (B2) as for the Laughlin-type wavefunction at the the...

    work page

  7. [7]

    For Laughlin-type wavefunction, we know that this quantity should be a uniform throughout the droplet which is a disk

    Alternative Method for Kinetic Energy Sampling The kinetic energy of the ansatz wavefunction can be calculated in a different way, naley via computing its Green’s functions, which we will quote the result [17]: GI0(z,z∗,˜︁z,˜︁z∗) (B10) =C (︁ 1 +g 2(˜︁z∗−z∗)(z−˜︁z) +g 4(˜︁z∗−z∗)2(z−˜︁z)2)︁ e ∑︁ IJ πnJK+ I0 Jz˜︁z∗−|z|2 4l2 I0 e ∑︁ IJ πnJK− I0 J ˜︁zz∗−|˜︁z|2...

    work page

  8. [8]

    For Pfaffian and QFL states,n k must be less then 1

    = (k2 2−k 2 1)2 12 (C2) We find nk = 4πne k2 2−k 2 1 = 2πne√︁ 3(⟨k4⟩−⟨k 2⟩2) (C3) ForK 2a andK 2b states,n k must be less then 2. For Pfaffian and QFL states,n k must be less then 1. Indeed, our numerical results satisfy this condition. From k2 2−k 2 1 = 4πne nk , k 2 2 +k 2 1 = 2⟨k2⟩,(C4) we find k2 2 4πne = ⟨k2⟩ 4πne + 1 2nk , k2 1 4πne = ⟨k2⟩ 4πne − 1 ...

    work page

Show all 57 references
  1. [9]

    Balents, C

    L. Balents, C. R. Dean, D. K. Efetov, and A. F. Young, Superconductivity and strong correlations in moir´ e flat bands, Nature Physics16, 725 (2020)

    2020

  2. [10]

    K. P. Nuckolls and A. Yazdani, A microscopic perspec- tive on moir´ e materials, Nature Reviews Materials9, 460 (2024)

    2024

  3. [11]

    E. Y. Andrei and A. H. MacDonald, Graphene bilayers with a twist, Nature Materials19, 1265 (2020)

    2020

  4. [12]

    H. Zhou, L. Holleis, Y. Saito, L. Cohen, W. Huynh, C. L. Patterson, F. Yang, T. Taniguchi, K. Watanabe, and A. F. Young, Isospin mag- netism and spin-polarized superconductivity in bernal bilayer graphene, Science375, 774 (2022), https://www.science.org/doi/pdf/10.1126/science.abm8386

    2022 doi

  5. [13]

    Zhang, R

    Y. Zhang, R. Polski, A. Thomson, ´E. Lantagne- Hurtubise, C. Lewandowski, H. Zhou, K. Watan- abe, T. Taniguchi, J. Alicea, and S. Nadj-Perge, En- hanced superconductivity in spin–orbit proximitized bi- layer graphene, Nature613, 268 (2023)

    2023

  6. [14]

    Holleis, C

    L. Holleis, C. L. Patterson, Y. Zhang, Y. Vituri, H. M. Yoo, H. Zhou, T. Taniguchi, K. Watanabe, E. Berg, S. Nadj-Perge, and A. F. Young, Nematicity and orbital depairing in superconducting bernal bilayer graphene, Nature Physics21, 444 (2025). 13

    2025

  7. [15]

    C. Li, F. Xu, B. Li, J. Li, G. Li, K. Watanabe, T. Taniguchi, B. Tong, J. Shen, L. Lu, J. Jia, F. Wu, X. Liu, and T. Li, Tunable superconductivity in electron- and hole-doped bernal bilayer graphene, Nature631, 300 (2024)

    2024

  8. [16]

    Zhang, G

    Y. Zhang, G. Shavit, H. Ma, Y. Han, C. W. Siu, A. Mukherjee, K. Watanabe, T. Taniguchi, D. Hsieh, C. Lewandowski, F. von Oppen, Y. Oreg, and S. Nadj- Perge, Twist-programmable superconductivity in spin– orbit-coupled bilayer graphene, Nature641, 625 (2025)

    2025

  9. [17]

    H. Zhou, T. Xie, T. Taniguchi, K. Watanabe, and A. F. Young, Superconductivity in rhombohedral trilayer graphene, Nature598, 434 (2021)

    2021

  10. [18]

    J. Yang, X. Shi, S. Ye, C. Yoon, Z. Lu, V. Kakani, T. Han, J. Seo, L. Shi, K. Watanabe, T. Taniguchi, F. Zhang, and L. Ju, Impact of spin-orbit coupling on superconductiv- ity in rhombohedral graphene (2025), arXiv:2408.09906 [cond-mat.supr-con]

    2025 arXiv

  11. [19]

    C. L. Patterson, O. I. Sheekey, T. B. Arp, L. F. W. Holleis, J. M. Koh, Y. Choi, T. Xie, S. Xu, E. Re- dekop, G. Babikyan, H. Zhou, X. Cheng, T. Taniguchi, K. Watanabe, C. Jin, E. Lantagne-Hurtubise, J. Alicea, and A. F. Young, Superconductivity and spin canting in spin-orbit ...

    2024 doi arXiv

  12. [20]

    T. Han, Z. Lu, Y. Yao, L. Shi, J. Yang, J. Seo, S. Ye, Z. Wu, M. Zhou, H. Liu, G. Shi, Z. Hua, K. Watanabe, T. Taniguchi, P. Xiong, L. Fu, and L. Ju, Signatures of Chiral Superconductivity in Rhombohedral Graphene 10.48550/arXiv.2408.15233 (2024), arXiv:2408.15233

    2024 doi arXiv

  13. [21]

    C. L. Patterson, O. I. Sheekey, T. B. Arp, L. F. W. Holleis, J. M. Koh, Y. Choi, T. Xie, S. Xu, Y. Guo, H. Stoyanov, E. Redekop, C. Zhang, G. Babikyan, D. Gong, H. Zhou, X. Cheng, T. Taniguchi, K. Watan- abe, M. E. Huber, C. Jin, ´E. Lantagne-Hurtubise, J. Al- icea, and A. F. ...

    2025

  14. [22]

    Morissette, P

    E. Morissette, P. Qin, H.-T. Wu, N. J. Zhang, K. Watan- abe, T. Taniguchi, and J. I. A. Li, Superconductivity, anomalous hall effect, and stripe order in rhombohe- dral hexalayer graphene (2025), arXiv:2504.05129 [cond- mat.mes-hall]

    2025 arXiv

  15. [23]

    Y.-Z. Chou, J. Zhu, and S. D. Sarma, Intravalley spin- polarized superconductivity in rhombohedral tetralayer graphene (2024), arXiv:2409.06701

    2024 arXiv

  16. [24]

    Geier, M

    M. Geier, M. Davydova, and L. Fu, Chiral and topo- logical superconductivity in isospin polarized multilayer graphene (2024), arXiv:2409.13829

    2024 arXiv

  17. [25]

    M. Kim, A. Timmel, L. Ju, and X.-G. Wen, Topological chiral superconductivity beyond pairing in a Fermi liquid, Phys. Rev. B (2025), arXiv:2409.18067

    2025 arXiv

  18. [26]

    Jahin and S.-Z

    A. Jahin and S.-Z. Lin, Enhanced Kohn-Luttinger topological superconductivity in bands with nontriv- ial geometry, arXiv e-prints , arXiv:2411.09664 (2024), arXiv:2411.09664

    2024

  19. [27]

    J. D. Sau and S. Wang, Theory of anomalous hall effect from screened vortex charge in a phase disordered super- conductor (2024), arXiv:2411.08969 [cond-mat.supr-con]

    2024 arXiv

  20. [28]

    Yang and Y.-H

    H. Yang and Y.-H. Zhang, Topological incommensu- rate fulde-ferrell-larkin-ovchinnikov superconductor and bogoliubov fermi surface in rhombohedral tetra-layer graphene (2024), arXiv:2411.02503 [cond-mat.supr-con]

    2024 arXiv

  21. [29]

    Qin and C

    Q. Qin and C. Wu, Chiral finite-momentum su- perconductivity in the tetralayer graphene (2024), arXiv:2412.07145 [cond-mat.supr-con]

    2024 arXiv

  22. [30]

    C. Yoon, T. Xu, Y. Barlas, and F. Zhang, Quarter metal superconductivity (2025), arXiv:2502.17555 [cond- mat.mes-hall]

    2025

  23. [31]

    Parra-Martinez, A

    G. Parra-Martinez, A. Jimeno-Pozo, V. T. Phong, H. Sainz-Cruz, D. Kaplan, P. Emanuel, Y. Oreg, P. A. Pantaleon, J. A. Silva-Guillen, and F. Guinea, Band renormalization, quarter metals, and chiral supercon- ductivity in rhombohedral tetralayer graphene (2025), arXiv:2502.19474...

    2025

  24. [32]

    Sedov and M

    D. Sedov and M. S. Scheurer, Probing superconductivity with tunneling spectroscopy in rhombohedral graphene (2025), arXiv:2503.12650 [cond-mat.supr-con]

    2025 arXiv

  25. [33]

    Y. Chen, M. S. Scheurer, and C. Schrade, Intrinsic su- perconducting diode effect and nonreciprocal supercon- ductivity in rhombohedral graphene multilayers (2025), arXiv:2503.16391 [cond-mat.supr-con]

    2025 arXiv

  26. [34]

    Christos, P

    M. Christos, P. M. Bonetti, and M. S. Scheurer, Finite- momentum pairing and superlattice superconductivity in valley-imbalanced rhombohedral graphene (2025), arXiv:2503.15471 [cond-mat.str-el]

    2025

  27. [35]

    Y. H. Chen, F. Wilczek, E. Witten, and B. Halperin, J. Mod. Phys. B3, 1001 (1989)

    1989

  28. [36]

    Lee, Anyon superconductivity and the fractional quantum hall effect, Physica B: Condensed Matter169, 37 (1991)

    D.-H. Lee, Anyon superconductivity and the fractional quantum hall effect, Physica B: Condensed Matter169, 37 (1991)

    1991

  29. [37]

    Wen and A

    X.-G. Wen and A. Zee, Topological structures, universal- ity classes, and statistics screening in the anyon super- fluid, Phys. Rev. B44, 274 (1991)

    1991

  30. [38]

    Wiegmann, Topological Superconductivity, Progress of Theoretical Physics Supplement107, 243 (1992)

    P. Wiegmann, Topological Superconductivity, Progress of Theoretical Physics Supplement107, 243 (1992)

    1992

  31. [39]

    Lee, Pairing via index theorem, Phys

    D.-H. Lee, Pairing via index theorem, Phys. Rev. B60, 12429 (1999), arXiv:cond-mat/9902287

    1999 arXiv

  32. [40]

    W.-H. Ko, P. A. Lee, and X.-G. Wen, Doped Kagome sys- tem as exotic superconductor, Phys. Rev. B79, 214502 (2009), arXiv:0804.1359

    2009 arXiv

  33. [41]

    Tang and X.-G

    E. Tang and X.-G. Wen, Superconductivity with intrinsic topological order induced by pure coulomb interaction and time-reversal symmetry breaking, Phys. Rev. B88, 195117 (2013), arXiv:1306.1528

    2013 arXiv

  34. [42]

    Wen, Vacuum degeneracy of chiral spin states in compactified space, Phys

    X.-G. Wen, Vacuum degeneracy of chiral spin states in compactified space, Phys. Rev. B40, 7387 (1989)

    1989

  35. [43]

    Wen, Topological orders in rigid states, Int

    X.-G. Wen, Topological orders in rigid states, Int. J. Mod. Phys. B04, 239 (1990)

    1990

  36. [44]

    Read and D

    N. Read and D. Green, Paired states of fermions in two dimensions with breaking of parity and time- reversal symmetries and the fractional quantum Hall ef- fect, Physical Review B61, 10267 (2000), arXiv:cond- mat/9906453

    2000

  37. [45]

    Ceperley, Ground state of the fermion one-component plasma: A monte carlo study in two and three dimen- sions, Phys

    D. Ceperley, Ground state of the fermion one-component plasma: A monte carlo study in two and three dimen- sions, Phys. Rev. B18, 3126 (1978)

    1978

  38. [46]

    Tanatar and D

    B. Tanatar and D. M. Ceperley, Ground state of the two- dimensional electron gas, Phys. Rev. B39, 5005 (1989)

    1989

  39. [47]

    May-Mann, T

    J. May-Mann, T. Helbig, and T. Devakul, How pairing mechanism dictates topology in valley- polarized superconductors with berry curvature (2025), arXiv:2503.05697 [cond-mat.supr-con]

    2025 arXiv

  40. [48]

    Wen and A

    X.-G. Wen and A. Zee, Compressibility and superfluidity in the fractional-statistics liquid, Phys. Rev. B41, 240 (1990)

    1990

  41. [49]

    Wen, Topological orders and edge excitations in fractional quantum Hall states, Adv

    X.-G. Wen, Topological orders and edge excitations in fractional quantum Hall states, Adv. Phys.44, 405 14 (1995), arXiv:cond-mat/9506066

    1995 arXiv

  42. [50]

    Z. Lu, T. Han, Y. Yao, A. P. Reddy, J. Yang, J. Seo, K. Watanabe, T. Taniguchi, L. Fu, and L. Ju, Frac- tional Quantum Anomalous Hall Effect in a Graphene Moire Superlattice 10.48550/arXiv.2309.17436 (2023), arXiv:2309.17436

    2023 doi arXiv

  43. [51]

    Nijboer and F

    B. Nijboer and F. De Wette, On the calculation of lattice sums, Physica23, 309 (1957)

    1957

  44. [52]

    S. G. Brush, H. L. Sahlin, and E. Teller, Monte carlo study of a one-component plasma. i, The Journal of Chemical Physics45, 2102 (1966), https://pubs.aip.org/aip/jcp/article- pdf/45/6/2102/18845826/2102 1 online.pdf

    1966

  45. [53]

    Geier, K

    M. Geier, K. Nazaryan, T. Zaklama, and L. Fu, Is atten- tion all you need to solve the correlated electron problem? (2025), arXiv:2502.05383 [cond-mat.str-el]

    2025 arXiv

  46. [54]

    Gaskell, The collective treatment of a fermi gas: Ii, Proceedings of the Physical Society77, 1182 (1961)

    T. Gaskell, The collective treatment of a fermi gas: Ii, Proceedings of the Physical Society77, 1182 (1961)

    1961

  47. [55]

    Gaskell, The collective treatment of many-body sys- tems: Iii, Proceedings of the Physical Society80, 1091 (1962)

    T. Gaskell, The collective treatment of many-body sys- tems: Iii, Proceedings of the Physical Society80, 1091 (1962)

    1962

  48. [56]

    R. B. Laughlin, Anomalous quantum Hall effect: An in- compressible quantum fluid with fractionally charged ex- citations, Physical Review Letter50, 1395 (1983)

    1983

  49. [57]

    Morf and B

    R. Morf and B. I. Halperin, Monte carlo evaluation of trial wave functions for the fractional quantized hall ef- fect: Disk geometry, Phys. Rev. B33, 2221 (1986)

    1986

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.