pith. sign in

arxiv: 1907.10236 · v1 · pith:R2TTBGOUnew · submitted 2019-07-24 · ❄️ cond-mat.supr-con

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

Shi-Chao Fang , Yan Zhang , Xiaojun Zheng , Guangkun Liu , Zhongbing Huang This is my paper

Pith reviewed 2026-05-24 16:56 UTC · model grok-4.3

classification ❄️ cond-mat.supr-con
keywords electric fieldchiral d+id superconductivityAA-stacked bilayer grapheneHubbard modelquantum Monte Carlopairing correlationshalf fillingsuperconducting state
0
0 comments X

The pith

An electric field induces dominant chiral d+id superconducting pairing in AA-stacked bilayer graphene at half filling.

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

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.

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.

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.

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.

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.

Figures

Figures reproduced from arXiv: 1907.10236 by Guangkun Liu, Shi-Chao Fang, Xiaojun Zheng, Yan Zhang, Zhongbing Huang.

Figure 1
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 … view at source ↗
Figure 2
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 ,… view at source ↗
Figure 3
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 … view at source ↗
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 = 4 U/t = 3.0 10-3 0.0 0.2 0.4 0.6 0.8 1.0 -0.03 0.00 0.03 0.06 V (R>3) /t NN-s NN-d+id NN-p+ip NNN-p+ip NNN-f (b) L = 5 U/t = 3.0 [PITH_FU… view at source ↗
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] view at source ↗
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 lines represent the position of 0.0. 3.2. Impact of electric field on electron correlations Firstly, the impact of potential difference ε on var… view at source ↗
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 possible to accurately determine the pairing form in the system from the point of view of the traditional pairing correlation function. The co… view at source ↗
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 of pairing range on the lattice size L = 6 in [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

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

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.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 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.
axioms (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.

pith-pipeline@v0.9.0 · 5651 in / 1192 out tokens · 32632 ms · 2026-05-24T16:56:30.868580+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

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

  9. [9]

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

    work page 2012

  10. [10]

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

    work page 2014

  11. [11]

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

    work page 2012

  12. [12]

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

    work page 2014

  13. [13]

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

    work page 2010

  14. [14]

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

    work page 2011

  15. [15]

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

    work page 2018

  16. [16]

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

    work page 2007

  17. [17]

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

    work page 2012

  18. [18]

    Honerkamp C 2008 Phys. Rev. Lett. 100 146404

    work page 2008

  19. [19]

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

    work page 2010

  20. [20]

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

    work page 2014

  21. [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

    work page 2017

  22. [22]

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

    work page 2017

  23. [23]

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

    work page 2015

  24. [24]

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

    work page 2016

  25. [25]

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

    work page 2016

  26. [26]

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

    work page 2000

  27. [27]

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

    work page 1991

  28. [28]

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

    work page 1991

  29. [29]

    Phys.: Condens

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

    work page 1989

  30. [30]

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

    work page 2008

  31. [31]

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

    work page 2018

  32. [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

    work page 2018

  33. [33]

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

    work page 2018

  34. [34]

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

    work page 2019

  35. [35]

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

    work page 2018

  36. [36]

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

    work page 2018

  37. [37]

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

    work page 2018

  38. [38]

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

    work page 2018

  39. [39]

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

    work page 2018

  40. [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

    work page arXiv 1903

  41. [41]

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

    work page 2009

  42. [42]

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

    work page 2016

  43. [43]

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

    work page 2012

  44. [44]

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

    work page 2008

  45. [45]

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

    work page 2014

  46. [46]

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

    work page 2013

  47. [47]

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

    work page 2013

  48. [48]

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

    work page 1997

  49. [49]

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

    work page 1997

  50. [50]

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

    work page 2001

  51. [51]

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

    work page 2001

  52. [52]

    Phys.: Condens

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

    work page 2014

  53. [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

    work page 2016

  54. [54]

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

    work page 1992

  55. [55]

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

    work page 2002

  56. [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

    work page 2009

  57. [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

    work page 2007

  58. [58]

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

    work page 1989