Pith. sign in

REVIEW 2 major objections 1 minor 58 references

Electric field-induced chiral d+id superconducting state in AA-stacked bilayer graphene: A quantum Monte Carlo study

T0 review · 2 major / 1 minor · reviewed 2026-05-24 · grok-4.3

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

Pith's one-line read An electric field induces dominant chiral d+id superconducting pairing in AA-stacked bilayer graphene at half filling.

desk verdict CPQMC on AA bilayer Hubbard model finds electric field favors d+id pairing at half filling, but the nodal constraint could be selecting the channel. read the letter →

arxiv 1907.10236 v1 pith:R2TTBGOU submitted 2019-07-24 cond-mat.supr-con

Shi-Chao Fang , Yan Zhang , Xiaojun Zheng , Guangkun Liu , Zhongbing Huang This is my paper · ORCID
classification cond-mat.supr-con
keywords electricfieldchirald+idsuperconductivityAA-stackedbilayergrapheneHubbardmodelquantumMonteCarlopairingcorrelationshalffillingsuperconductingstate
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

The paper applies the constrained-path quantum Monte Carlo method to the Hubbard model on AA-stacked honeycomb lattices in the presence of an electric field. It reports that the field stabilizes a chiral d+id wave pairing channel at half filling, with the pairing correlations growing as the on-site repulsion is increased. The authors link this to an electric-field-driven rise in the density of states near the Fermi energy together with a suppression of antiferromagnetic spin correlations. A reader would care because the result identifies a concrete, gate-tunable route to chiral superconductivity in a minimal graphene lattice without requiring doping.

What carries the argument

Constrained-path quantum Monte Carlo applied to the electric-field-tuned Hubbard model on the AA-stacked honeycomb lattice, used to extract the dominant pairing correlations in the d+id channel.

What would settle it

An unbiased calculation such as exact diagonalization on small clusters that shows no dominant d+id pairing signal once the electric field is applied would falsify the central claim.

Watch

Extended reading notes

Core claim

Our simulation demonstrates a dominant chiral d+id wave pairing induced by the electric field at half filling. In particular, as the on-site Coulomb interaction increases, the effective pairing correlation of chiral d+id superconducting state exhibits increasing behavior. We attribute the electric field induced d+id superconductivity to an increased density of states near the Fermi energy and a suppressed antiferromagnetic spin correlation after turning on the electric field. Our results strongly suggest the AA-stacked graphene system with electric field is a good candidate for chiral d+id superconductors.

Load-bearing premise

The constrained-path quantum Monte Carlo method accurately captures the pairing correlations without introducing significant bias from the sign or phase constraint in this field-tuned system.

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, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 1 minor

Summary. The manuscript applies the constrained-path quantum Monte Carlo method to the Hubbard model on AA-stacked bilayer graphene subject to a perpendicular electric field. It reports that the field induces dominant chiral d+id pairing correlations at half filling, with these correlations strengthening as the on-site repulsion U increases; the effect is attributed to an enhanced density of states at the Fermi level together with suppressed antiferromagnetic spin correlations.

Significance. If the numerical results prove robust against methodological bias, the work supplies non-perturbative evidence that an external electric field can stabilize chiral d+id superconductivity in a graphene multilayer, offering a concrete route to field-tunable topological pairing. The direct simulation of the microscopic Hamiltonian, free of fitted parameters beyond U and the field strength, constitutes a clear strength.

major comments (2)
  1. [Simulation method and pairing-correlation analysis] The central claim of d+id dominance rests on CPQMC pairing correlations. The manuscript does not report any diagnostic that varies the symmetry or nodal structure of the trial wavefunction used to impose the constrained-path approximation while the electric field is applied; because the field breaks layer equivalence and shifts the single-particle spectrum, a mismatch between trial and true nodal surface can preferentially weight one pairing channel over others (e.g., d+id versus s or p). This test is load-bearing for the reported channel selection.
  2. [Discussion of physical mechanism] The attribution of enhanced d+id correlations to increased DOS and suppressed AF order is stated qualitatively. No quantitative comparison (e.g., field-induced change in DOS extracted from the single-particle spectrum or AF structure factor versus field strength) is provided to establish that these mechanisms, rather than the constraint itself, drive the observed trend with U.
minor comments (1)
  1. [Figure captions and methods summary] Lattice sizes, inverse temperatures, and statistical error bars on the pairing correlations should be stated explicitly in the figure captions or a methods table so that convergence can be assessed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address the two major comments below and will incorporate the suggested checks and quantitative analysis into a revised manuscript.

read point-by-point responses

  1. Referee: The central claim of d+id dominance rests on CPQMC pairing correlations. The manuscript does not report any diagnostic that varies the symmetry or nodal structure of the trial wavefunction used to impose the constrained-path approximation while the electric field is applied; because the field breaks layer equivalence and shifts the single-particle spectrum, a mismatch between trial and true nodal surface can preferentially weight one pairing channel over others (e.g., d+id versus s or p). This test is load-bearing for the reported channel selection.

    Authors: We agree that testing the sensitivity of the constrained-path approximation to the trial-wavefunction nodal structure is important when the electric field breaks layer symmetry. Our calculations employ a trial wave function obtained from the non-interacting Hamiltonian that already includes the perpendicular field, thereby incorporating the correct single-particle spectrum and layer asymmetry. Nevertheless, to rule out bias toward the d+id channel, we will perform additional runs with trial states that impose alternative pairing symmetries (e.g., s-wave or p-wave nodes) and report the resulting pairing correlations. These diagnostics will be added to the revised manuscript. revision: yes

  2. Referee: The attribution of enhanced d+id correlations to increased DOS and suppressed AF order is stated qualitatively. No quantitative comparison (e.g., field-induced change in DOS extracted from the single-particle spectrum or AF structure factor versus field strength) is provided to establish that these mechanisms, rather than the constraint itself, drive the observed trend with U.

    Authors: The referee correctly notes that the mechanistic discussion remains qualitative. Although the AF structure factor is computed in our simulations and shows suppression with increasing field, and the density of states can be extracted from the single-particle Green's function, we did not present explicit quantitative correlations between these quantities and the pairing strength. In the revision we will add plots of the field dependence of the DOS at the Fermi level (obtained from both non-interacting and QMC spectra) together with the AF structure factor, and we will overlay these against the d+id pairing correlations to provide a quantitative link. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: results are direct outputs of numerical simulation

full rationale

The paper's central claim of electric-field-induced dominant chiral d+id pairing at half filling is obtained by running constrained-path quantum Monte Carlo on the Hubbard Hamiltonian for the AA-stacked lattice; pairing correlations are measured directly from the sampled configurations. No parameter is fitted to the target pairing channel and then re-reported as a prediction, no self-citation supplies a uniqueness theorem that forces the result, and no ansatz or known empirical pattern is renamed as a derivation. The simulation outputs are independent of the interpretive statements about DOS and AF suppression.

Assumptions & free parameters 2 free parameters · 1 assumptions · 0 invented entities

The central claim rests on the output of constrained-path QMC applied to the Hubbard model with an added electric-field term; the model itself and the numerical method are the main inputs.

free parameters (2)
  • on-site Coulomb interaction U
    Varied across simulations; pairing correlation increases with U.
  • electric field strength
    Applied to induce the reported superconducting state.
assumptions (1)
  • domain assumption AA-stacked bilayer graphene is adequately described by the Hubbard model on a honeycomb lattice plus an electric-field term.
    This is the Hamiltonian used for all simulations.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Electric field-induced chiral d+id superconducting state in AA-stacked bilayer graphene: A quantum Monte Carlo study." pith.science (2026). https://pith.science/paper/R2TTBGOU

@misc{pith2026190710236,
  author       = {Pith},
  title        = {Pith review of: Electric field-induced chiral d+id superconducting state in AA-stacked bilayer graphene: A quantum Monte Carlo study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R2TTBGOU}},
  note         = {Machine review of arXiv:1907.10236}
}
read the original abstract

Using constrained-path quantum Monte Carlo method, we systematically study the Hubbard model on AA-stacked honeycomb lattices with electric field. Our simulation demonstrates a dominant chiral d+id wave pairing induced by the electric field at half filling. In particular, as the on-site Coulomb interaction increases, the effective pairing correlation of chiral d+id superconducting state exhibits increasing behavior. We attribute the electric field induced d+id superconductivity to an increased density of states near the Fermi energy and an suppressed antiferromagnetic spin correlation after turning on the electric field. Our results strongly suggest the AA-stacked graphene system with electric field is a good candidate for chiral d+id superconductors.

Figures

Figures reproduced from arXiv: 1907.10236 by the authors.

Figure 1
Figure 1. (color online) (a) Sketch of an AA-stacked honeycomb lattice. Blue (red) dots represent sublattice A (B). A1,B1 (A2,B2) represent sublattices in layer 1 (layer 2). t and t⊥ are the intralayer and interlayer hoppings, respectively. (b) Geometry of each graphene layer. Atom number on each layer is 2 × 3L 2 . The presenting lattice is corresponding to L = 4. Hubbard model can be described as follows, H = − t X hi,ji,m,… view at source ↗
Figure 2
Figure 2. (color online) Sketch of calculated intralayer pairing channels. (a) nearest￾neighbor (NN) s-wave (b) NN d+id-wave (c) NN p+ip-wave (d) next-nearest-neighbor (NNN) p + ip wave (e) NNN f-wave pairings. where ∆† α (i) (∆α(i)) is the electron pair creation (annihilation) operator with pairing symmetry α. Singlet or triplet pair creation operator can be written as, ∆ † α (i) = 1 √ Nα X l f † α (δl)(ci,↑ci+δl ,↓ ∓ ci,↓ci… view at source ↗
Figure 3
Figure 3. (color online) (a) Averaged electron density of each layer hnmi versus potential difference ε for different lattice sizes. The filled symbols and open symbols correspond to the lattice L = 4 and L = 5, respectively. (b) Effective electron doping density hn1(ε)i − 1 of layer 1 for various potential difference ε using 3rd order polynomials in 1/ √ N, where N is the each layer atom number of the system and N = 6L 2 (L … view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: (color online) Long-range averaged pairing correlations of different pairing channels versus potential difference ε on different Hubbard U for the lattice size L = 4 (a) U/t = 0.0 and (b) U/t = 3.0. 0.0 0.2 0.4 0.6 0.8 1.0 -0.10 -0.05 0.00 0.05 0.10 10-3 V (R>3) (a) L …
Figure 5
Figure 5. Figure 5: (color online) Long-distance averaged effective pairing correlations of different pairing channels versus potential difference ε at U/t = 3.0 on different lattice sizes (a) L = 4 and (b) L = 5 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: (color online) (a) Long-distance effective pairing correlation of d+id pairing channels versus distance R between pairs on the lattice size L = 6 for different potential difference ε at U/t = 3.0 and (b) different on-site interaction U at ε/t = 0.8. The green dash line…
Figure 7
Figure 7. Figure 7: (color online) Band structures, Fermi surfaces and density of states of the non-interacting Hubbard model on the AA-stacked bilayer honeycomb lattice on various potential difference ε. The green dash lines represent the position of Fermi levels. Therefore, it is not po…
Figure 8
Figure 8. Figure 8: (color online) The intralayer NN spin correlation function versus potential difference ε at U/t = 3.0 for different lattice sizes.(a) L = 4 and (b) L = 5. repulsive interaction U on the dominant d+id-wave, we also plotted the effective pairing correlation as a function…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

58 extracted references · 58 canonical work pages

  1. [1]

    Dagotto E 2013 Rev. Mod. Phys. 85 849

    work page 2013

  2. [2]

    Johnston D C 2010 Adv. Phys. 59 803

    work page 2010

  3. [3]

    Mackenzie A P and Maeno Y 2003 Rev. Mod. Phys. 75 657

    work page 2003

  4. [4]

    Scalapino D J 2010 Rev. Mod. Phys. 84 1383

    work page 2010

  5. [5]

    Dagotto E 1994 Rev. Mod. Phys. 66 763

    work page 1994

  6. [6]

    Qi X L and Zhang S C 2011 Rev. Mod. Phys. 83 057

    work page 2011

  7. [7]

    Mazin I I, Singh D J, Johannes M D, and Du M H 2008 Phys. Rev. Lett . 101 057003

    work page 2008

  8. [8]

    Castro Neto A H, Guinea F, Peres N M R, Novoselov K S, and Geim A K 2 009 Rev. Mod. Phys. 81 109

    work page

Show all 58 references
  1. [9]

    Black-Schaffer A M 2012 Phys. Rev. Lett. 109 197001

    2012

  2. [10]

    Black-Schaffer A M, Wu W, Hur K L 2014 Phys. Rev. B 90 054521

    2014

  3. [11]

    Nandkishore R, Levitov L S, Chubukov A V 2012 Nature Physics 8 158C163

    2012

  4. [12]

    Nandkishore R, Thomale Ronny, Chubukov A V 2014 Phys. Rev. B 89 144501

    2014

  5. [13]

    Pathak S, Shenoy V B, Baskaran G 2010 Phys. Rev. B 81 085431

    2010

  6. [14]

    Ma T, Huang Z, Hu F, Lin H Q 2011 Phys. Rev. B 84 121410(R) 12

    2011

  7. [15]

    Ying T, Wessel S 2018 Phys. Rev. B 97 075127

    2018

  8. [16]

    Black-Schaffer A M, Doniach Sebastian 2007 Phys. Rev. B 75 134512

    2007

  9. [17]

    Kiesel M L, Platt M,Hanke W, Abanin D A, Thomale R 2012 Phys. Rev. B 86 020507(R)

    2012

  10. [18]

    Honerkamp C 2008 Phys. Rev. Lett. 100 146404

    2008

  11. [19]

    Raghu S, Kivelson S A, Scalapino 2010 Phys. Rev. B 81 224505

    2010

  12. [20]

    Ma T X, Yang F, Yao H, and Lin H Q 2014 Phys. Rev. B 90 245114

    2014

  13. [21]

    Di Bernardo A, Millo O, Barbone M, Alpern H, Kalcheim Y, Sassi U, O tt A K, Fazio D D, Yoon D, Amado M, Ferrari A C, Linder J and Robinson J W A 2017 Nature Com munications 8 14024

    2017

  14. [22]

    Ma T X, Yang F, Huang Z B and Lin H Q 2017 Scientific Report 7 19

    2017

  15. [23]

    Faye J P L, Sahebsara P, and Snchal D 2015 Phys. Rev. B 92 085 121

    2015

  16. [24]

    Xiao L Y, Yu S L, Wang W, Yao Z J and Li J X 2016 EPL (Europhysics L etters) 115 2

    2016

  17. [25]

    Xu X Y, Wessel S, and Meng Z Y 2016 Phys. Rev. B 94 115105

    2016

  18. [26]

    Vojta M, Zhang Y, and Sachdev S 2000 Phys. Rev. Lett. 85 4940

    2000

  19. [27]

    Moore Gregory, Read Nicholas 1991 Nuclear Physics B 360 362-396

    1991

  20. [28]

    Kopnin N B, Salomaa M M 1991 Phys. Rev. B 44 9667

    1991

  21. [29]

    Phys.: Condens

    Volovik G E, Yakovenko V M 1989 J. Phys.: Condens. Matter, 1 5263

    1989

  22. [30]

    Gonz´alez J 2008 Phys. Rev. B 78 205431

    2008

  23. [31]

    Cao Y, Fatemi V, Fang S, Watanabe K, Taniguchi T, Kaxiras E an d Jarillo-Herrero P 2018 Nature 556 43C50

    2018

  24. [32]

    Cao Y, Fatemi V, Demir A, Fang S, Tomarken S L, Luo J Y, Sanche z-Yamagishi J D, Watanabe K, Taniguchi T, Kaxiras E, Ashoori R C and Jarillo-Herrero P 2018 Na ture 556 80C84

    2018

  25. [33]

    Isobe H, Yuan N F Q ,Fu L 2018 Phys. Rev. X 8 041041

    2018

  26. [34]

    Gonz´alez J, Stauber T 2019 Phys. Rev. Lett. 122 026801

    2019

  27. [35]

    Liu C C, Zhang L D, Chen W Q,Yang F 2018 Phys. Rev. Lett. 121 217001

    2018

  28. [36]

    Xu C, Balents L 2018 Phys. Rev. Lett. 121 087001

    2018

  29. [37]

    Kennes D M, Lischner Johannes, Karrasch C 2018 Phys. Rev. B 98 241407

    2018

  30. [38]

    Peltonen T J, Ojaj¨ arvi R, Heikkil¨a T 2018 Phys. Rev. B 98 220504(R)

    2018

  31. [39]

    Fidrysiak M, Zegrodnik M, Spalek J 2018 Phys. Rev. B 98 085436

    2018

  32. [40]

    Shen C, Li N, Wang S P, Zhao Y C, Tang J, Liu J Y, Tian J P, Chu Y B , W atanabe K, Taniguchi T, Yang R, Meng Z Y, Shi D X, Zhang G Y arXiv:1903.06952

    1903

  33. [41]

    Liu Z, Suenaga K, J.F.Harris P, and Iijima S 2009 Phys. Rev. Lett. 102 015501

    2009

  34. [42]

    Rozhkov A V, Sboychakov A O, Rakhmanov A L, Nori F 2016 Phys ics Reports 648 1-104

    2016

  35. [43]

    Rakhmanov A L, Rozhkov A V, Sboychakov A O, and Nori F 2012 P hys. Rev. Lett. 109 206801

    2012

  36. [44]

    P. L. de Andres, Ramrez R, and Vergs J A 2008 Phys. Rev. B 77 045403

    2008

  37. [45]

    Akzyanov R S, Sboychakov A O, Rozhkov A V, Rakhmanov A L, an d Nori F 2014 Phys. Rev. B 90 11415

    2014

  38. [46]

    Sboychakov A O, Rozhkov A V, Rakhmanov A L, and Nori F 2013 P hys. Rev. B 88 045409

    2013

  39. [47]

    Sboychakov A O, Rakhmanov A L, Rozhkov A V, and Nori F 2013 P hys. Rev. B 87 121401(R)

    2013

  40. [48]

    Zhang S, Carlson J, and Gubernatis J E 1997 Phys. Rev. B 55 7464

    1997

  41. [49]

    Zhang S, Carlson J, and Gubernatis J E 1997 Phys. Rev. Lett. 78 4486

    1997

  42. [50]

    Huang Z, Lin H Q and Gubernatis J E 2001 Phys. Rev. B 63 115112

    2001

  43. [51]

    Huang Z, Lin H Q and Gubernatis J E 2001 Phys. Rev. B 64 205101

    2001

  44. [52]

    Phys.: Condens

    Liu G K, Huang Z B and Wang Y J 2014 J. Phys.: Condens. Matter, 26 325601

    2014

  45. [53]

    Liu G K, Kaushal N, Li S L, Bishop C B, Wang Y, Johnston S, Alvare z G, Moreo A, and Dagotto E 2016 Phys. Rev. E 93 063313

    2016

  46. [54]

    Charlier J C, Michenaud J P and Gonze X 1992 Phys. Rev. B 46 4531

    1992

  47. [55]

    Dresselhaus M S and Dresselhaus G 2002 Adv. Phys. 51 1

    2002

  48. [56]

    Zhang Y, Tang T-T, Girit C, Hao Z, C.Martin M, Zettl A, F.Cormmie M , Shen Y R, Wang F 2009 Nature 459 820

    2009

  49. [57]

    V.Castro E, Novoselov K S, Morozov S V, Peres N M R, Lopes dos S antos J M R, Nilsson J, Guinea F, Geim A K, and Castro Neto A H 2007 Phys. Rev. Lett. 99 216802 13

    2007

  50. [58]

    White S R, Scalapino D J, Sugar R L, Bickers N E, Scalettar R T 1989 Phys. Rev. B 39 839

    1989

Pith tools

Reviewed May 24, 2026 · model on record in the stance chip above.