Pith. sign in

REVIEW 2 major objections 6 minor 42 references

Interacting fermions on a honeycomb bilayer show a continuous quantum phase transition with emergent relativistic symmetry, matching Gross-Neveu-Ising critical exponents within 5%.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · grok-4.5

2026-07-13 20:18 UTC pith:LXAG6VVF

load-bearing objection Clean large-scale QMC confirmation of emergent z=1 Gross-Neveu-Ising criticality in a spinless bilayer model; the <5% exponent claim is a bit tight for L≤27 but the assignment itself holds. the 2 major comments →

arxiv 2603.22259 v2 pith:LXAG6VVF submitted 2026-03-23 cond-mat.str-el

Emergent relativistic symmetry from interacting fermions on the honeycomb bilayer

Zi Hong Liu , Lukas Janssen This is my paper
classification cond-mat.str-el PACS 71.10.Fd71.30.+h05.30.Rt64.70.Tg
keywords honeycomb bilayerquadratic band touchingGross-Neveu-Isingemergent relativityquantum Monte Carlocharge density waveDirac fermionsdynamical critical exponent
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper maps the phase diagram of spinless fermions on the Bernal-stacked honeycomb bilayer using large-scale quantum Monte Carlo simulations free of the sign problem. Without interactions the bands touch quadratically, producing a finite density of states. Weak to intermediate repulsion splits each quadratic touching into four Dirac cones, while strong repulsion drives a layer-polarized charge-density-wave insulator that breaks Z2 layer inversion. The authors show that the intervening quantum phase transition is continuous and is controlled by emergent relativistic invariance: the measured correlation-length exponent, order-parameter anomalous dimension and fermion anomalous dimension all agree with the 2+1D Gross-Neveu-Ising universality class for eight two-component Dirac fermions to better than 5%. They also chart the finite-temperature crossover from a high-temperature non-relativistic (z=2) regime into an intermediate-temperature Dirac (z=1) regime, and the subsequent classical 2D Ising ordering transition. The result supplies clean numerical evidence that relativistic symmetry can emerge at a quantum critical point even when the microscopic spectrum is non-relativistic.

Core claim

The interaction-driven semimetal-to-insulator transition of spinless fermions on the Bernal-stacked honeycomb bilayer is continuous and belongs to the 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions; the extracted exponents 1/ν=1.065(48), η_φ+z=1.856(42) and η_ψ=0.0199(46) match theoretical predictions within less than 5%.

What carries the argument

Crossing-point finite-size scaling of the charge-correlation ratio Rc, the CDW order parameter and the quasiparticle weight, combined with Bayesian extrapolation of the form aL^{-p}, which yields the three independent critical exponents that diagnose Gross-Neveu-Ising criticality.

Load-bearing premise

The thermodynamic-limit exponents rest on power-law finite-size extrapolations performed on lattices no larger than L=27 and on the prior assignment that the dynamical exponent equals one.

What would settle it

A measurement of the same three critical exponents on substantially larger lattices (or with an independent method that does not assume z=1) that yields values differing from Gross-Neveu-Ising predictions by more than the quoted uncertainties.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Hall coefficient should cross over from RH∝T^{-1} (z=2) above the crossover scale to RH∝T^{-2} (z=1) below it.
  • Landau-level spectroscopy should switch from quadratic-band-touching scaling EN∝B√[N(N-1)] at high T to Dirac scaling EN∝√(BN) at intermediate T.
  • Transport and thermodynamic probes of Bernal bilayer graphene above a putative Mott transition can reveal the same intermediate-temperature relativistic window.
  • The same microscopic route—interaction-induced splitting of quadratic band touchings—can host Gross-Neveu-Ising criticality in other bilayer or multi-orbital systems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Because the microscopic model is sign-problem free, the same lattice can serve as a controlled testbed for higher-order corrections or multi-flavor Gross-Neveu fixed points that remain inaccessible in most continuum simulations.
  • The clean separation of the z=2, z=1 and Ising temperature scales suggests that finite-temperature transport or optical conductivity measurements could map the entire crossover diagram without needing ultra-low temperatures.
  • If residual corrections-to-scaling are larger than estimated, the apparent sub-5% agreement may shrink, providing a quantitative benchmark for future larger-scale simulations.

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 / 6 minor

Summary. The manuscript presents large-scale sign-problem-free determinantal quantum Monte Carlo simulations of spinless fermions on the Bernal-stacked honeycomb bilayer at half filling. In the noninteracting limit the spectrum has quadratic band touchings; weak-to-intermediate nearest-neighbor repulsion V splits each touching into four Dirac cones, while strong V drives a continuous transition into a layer-polarized Z2 charge-density-wave insulator. Crossing-point finite-size scaling on lattices up to 4 imes27² yields critical exponents 1/ν=1.065(48), η_φ+z=1.856(42) and η_ψ=0.0199(46) that the authors compare with the 2+1D Gross-Neveu-Ising class for eight two-component Dirac flavors, claiming agreement within <5%. Finite-temperature data map a 2D Ising melting line of the CDW and a higher crossover scale separating a high-T quadratic-band-touching regime (z=2) from an intermediate-T Dirac regime (z=1), supported by uniform charge susceptibility and single-particle spectral functions.

Significance. If the identification holds, the work supplies one of the cleanest numerical demonstrations of emergent relativistic (Lorentz) symmetry at a fermionic quantum critical point whose microscopic spectrum is nonrelativistic. The model is free of the fermion sign problem by Majorana reflection positivity, reaches comparatively large volumes, and reports three independent exponents together with a finite-T phase diagram that includes both the classical Ising line and a DSM–QBT crossover. The quantitative comparison with Gross-Neveu-Ising theory (Table I) and the explicit spectral evidence for interaction-generated Dirac cones are falsifiable and of direct interest for bilayer graphene and related Dirac materials. These strengths make the manuscript a substantial contribution to the literature on fermionic quantum criticality.

major comments (2)
  1. Abstract and Table I claim agreement of 1/ν, η_φ+z and η_ψ with Gross-Neveu-Ising values “within less than 5%.” The main-text numbers come from Bayesian power-law extrapolations of crossing-point estimators on L≤27 (Figs. 2b–d). Supplemental data-collapse analyses (Table S1) with L_min=21 produce systematically shifted central values (e.g. η_φ+z=1.773(7) versus main-text 1.856(42) and theory 1.868(4); V_c=0.919(1) versus 0.900(5)). The SM itself notes that corrections-to-scaling remain visible and prefers the crossing-point route, yet the abstract’s quantitative “<5%” statement does not reflect this procedure dependence. The Gross-Neveu-Ising assignment remains plausible within enlarged uncertainties, but the manuscript should either (i) qualify the precision claim to match the residual finite-size uncertainty or (ii) present a joint analysis that folds both estimators and their discrepa
  2. In the zero-temperature scaling analysis (text surrounding Figs. 2c–d), the order-parameter estimator is written η_φ(1/L)=−z−ln[m²_CDW(L+3)/m²_CDW(L)]/ln r with z=1 inserted by hand before reporting η_φ=0.856(42) and the combination η_φ+z. The subsequent claim that the QCP has z=1 therefore partly rests on the same universality-class assignment that the exponents are meant to test. Finite-T susceptibility scaling (χ_uni∝T in the intermediate regime) and spectral functions support z=1 away from the QCP, but an independent zero-temperature estimate of z (e.g. from imaginary-time correlation lengths or dynamical structure-factor collapse at V_c) would make the emergent-relativistic conclusion less circular. At minimum the manuscript should state clearly which observables fix z without presupposing Gross-Neveu-Ising.
minor comments (6)
  1. Abstract: “intermediates temperatures” → “intermediate temperatures.”
  2. Supplemental Material title uses “interacting electrons” while the main text uses “fermions”; unify terminology.
  3. Fig. 1(b) caption and main text refer to red squares/dots for T_cross inconsistently; align symbol description with the plotted markers.
  4. The equal-residual criterion used to locate T_cross (Figs. 4a,b and S6) is pragmatic but under-specified (choice of fit windows, sensitivity). A short sentence on robustness under window variation would help reproducibility.
  5. Table I reprints theoretical values from Ref. [24] (co-authored by one of the present authors). A brief note that those numbers are independent field-theoretic estimates (not fits to the present QMC data) would forestall any appearance of circularity.
  6. Units and conventions: t=1, t_⊥=t are fixed early; a single sentence reminding the reader that all energies are in units of t would aid readers skimming the finite-T figures.

Circularity Check

1 steps flagged

Mild self-citation of target GN-Ising exponents (and initial z=1) from co-authored Ref. [24]; QMC data and comparison remain independent.

specific steps
  1. self citation load bearing [Abstract; Table I caption and surrounding text (p. 3); finite-size estimator for η_φ (p. 3)]
    "Our results for the values of the correlation-length exponent 𝜈, the order-parameter anomalous dimension 𝜂_𝜙, and the fermion anomalous dimension 𝜂_𝜓 agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5% deviation. ... (bottom row, reprinted from Ref. [24]). ... η_𝜙(1/L)=−z−ln[m^{2}_CDW(...)]/ln r, where z corresponds to the dynamical critical exponent. We show below that the quantum critical point is characterized by z=1."

    The numerical target values (1/ν≈1.018, η_φ+z≈1.868, η_ψ≈0.0195) and the identification of the universality class (including the insertion of z=1 into the order-parameter estimator) are taken from the co-authored field-theory paper [24]. The QMC exponents are measured independently and then declared to match, so the 'confirmation of emergent relativistic symmetry' partially rests on that self-citation for its benchmark; the match is not forced by construction, but the comparison is not fully external.

full rationale

The paper's central numerical claim is an independent large-scale QMC measurement of the continuous SM-insulator transition, its critical exponents via crossing-point and data-collapse analyses, and the finite-T crossover scales. These are not defined in terms of the target values, nor are any fitted lattice parameters re-labeled as predictions. The only mild circularity is that the benchmark exponents and the assignment to the 8-flavor 2+1D Gross-Neveu-Ising class (including the a-priori z=1 used when converting m^{2}_CDW into η_φ) are taken from Ref. [24] (Ray-Vojta-Janssen), which shares an author. That citation supplies the comparison numbers and the expected universality class, but the QMC results themselves constitute a genuine, falsifiable test rather than a tautology. No uniqueness theorem is imported, no ansatz is smuggled, and no quantity is predicted from a fit to a closely related subset of the same data. Score 2 reflects this single non-load-bearing self-reference; the derivation chain is otherwise self-contained against the external field-theoretic benchmark.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claim rests on standard QMC technology, finite-size scaling hypotheses, and the continuum identification of the critical theory with eight Dirac flavors previously derived by RG. No new particles or forces are postulated; the only free parameters are the usual simulation cut-offs and the fitted location of the critical point itself.

free parameters (3)
  • V_c (critical interaction) = 0.900(5)
    Thermodynamic-limit critical coupling extracted from finite-size crossings of the charge correlation ratio; value V_c=0.900(5) is used to evaluate all exponents.
  • Trotter step Δτ = 0.1
    Fixed by hand at 0.1; claimed to produce negligible discretization error but not systematically extrapolated to zero.
  • projection length Θ = 2L
    Set to 2L for ground-state projection; sufficiency is asserted but not proven for every observable.
axioms (4)
  • domain assumption Majorana reflection positivity guarantees a non-negative fermion determinant for the chosen Hubbard-Stratonovich channel, rendering the model sign-problem free.
    Invoked in the Model section and SM Sec. I to justify the applicability of large-scale DQMC.
  • domain assumption Finite-size estimates of critical exponents approach their thermodynamic values as power laws aL^{-p} (or aL^{-ω}), allowing Bayesian extrapolation.
    Used throughout the zero-temperature scaling analysis (Figs. 2b–d and SM Sec. II).
  • domain assumption After interaction-induced splitting of each quadratic band-touching point, the low-energy theory consists of eight two-component Dirac fermions coupled to an Ising order parameter (Gross-Neveu-Ising).
    Taken from prior RG analyses (Refs. [24,25]) and used as the theoretical benchmark in Table I.
  • standard math Standard 2D Ising exponents (ν=1, η=1/4) govern the finite-temperature melting of the layer-polarized CDW.
    Assumed when collapsing the charge susceptibility data in Fig. 3.

pith-pipeline@v1.1.0-grok45 · 18849 in / 2886 out tokens · 36871 ms · 2026-07-13T20:18:52.711511+00:00 · methodology

0 comments
read the original abstract

We study the phase diagram of interacting spinless fermions on the honeycomb bilayer at charge neutrality using large-scale quantum Monte Carlo simulations. In the noninteracting limit, the low-energy spectrum features quadratically dispersing bands that touch at the corners of the hexagonal Brillouin zone. Weak to intermediate interactions induce a splitting of each of the quadratic band touching points into four Dirac points, located along high-symmetry directions of the reciprocal lattice. Strong interactions lead to the formation of a layer-polarized charge density wave, which spontaneously breaks the $\mathbb Z_2$ layer inversion symmetry and opens an insulating gap in the spectrum. We show that the semimetal-to-insulator quantum phase transition as a function of interaction is continuous and characterized by emergent relativistic symmetry. Our results for the values of the correlation-length exponent $\nu$, the order-parameter anomalous dimension $\eta_\phi$, and the fermion anomalous dimension $\eta_\psi$ agree with those of the theoretically predicted 2+1D Gross-Neveu-Ising universality class with eight two-component Dirac fermions within less than 5\%\ deviation. We also determine the crossover scale as a function of interaction strength between the nonrelativistic semimetal state at high temperatures, characterized by dynamical critical exponent $z = 2$, and the Dirac semimetal state at intermediates temperatures, characterized by $z=1$. Further reducing the temperature below the crossover scale at a fixed value of the interaction strength above the quantum critical point results in a classical ordering transition in the 2D Ising universality class.

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

42 extracted references · 2 linked inside Pith

  1. [1]

    Sachdev,Quantum Phase Transitions, 2nd ed

    S. Sachdev,Quantum Phase Transitions, 2nd ed. (Cambridge University Press, Cambridge, UK, 2011)

    2011

  2. [2]

    Pelissetto and E

    A. Pelissetto and E. Vicari, Critical phenomena and renormalization-group theory, Phys. Rep.368, 549 (2002)

    2002

  3. [3]

    J. Lou, A. W. Sandvik, and L. Balents, Emergence of U(1) Symmetry in the 3D𝑋𝑌Model with𝑍 𝑞 Anisotropy, Phys. Rev. Lett.99, 207203 (2007)

    2007

  4. [4]

    Patil, H

    P. Patil, H. Shao, and A. W. Sandvik, Unconventional U(1) to𝑍𝑞 crossover in quantum and classical𝑞-state clock models, Phys. Rev. B103, 054418 (2021)

    2021

  5. [5]

    Senthil, Deconfined Quantum Critical Points: A Review, in 50 Years of the Renormalization Group(World Scientific, 2023) pp

    T. Senthil, Deconfined Quantum Critical Points: A Review, in 50 Years of the Renormalization Group(World Scientific, 2023) pp. 169–195, arXiv:2306.12638

    Pith/arXiv arXiv 2023

  6. [6]

    B. Roy, V. Juri ˇci´c, and I. F. Herbut, Emergent Lorentz sym- metry near fermionic quantum critical points in two and three dimensions, J. High Energy Phys. 4 (2016) 18

    2016

  7. [7]

    Biedermann and L

    J. Biedermann and L. Janssen, Dirac quantum criticality in twisted double bilayer transition metal dichalcogenides, arXiv:2509.04561

    Pith/arXiv arXiv

  8. [8]

    Lee, Emergence of supersymmetry at a critical point of a lattice model, Phys

    S.-S. Lee, Emergence of supersymmetry at a critical point of a lattice model, Phys. Rev. B76, 075103 (2007)

    2007

  9. [9]

    Grover, D

    T. Grover, D. N. Sheng, and A. Vishwanath, Emergent Space- Time Supersymmetry at the Boundary of a Topological Phase, Science344, 280 (2014)

    2014

  10. [10]

    Schuler, L.-P

    M. Schuler, L.-P. Henry, Y.-M. Lu, and A. M. L¨auchli, Emergent XY* transition driven by symmetry fractionalization and anyon condensation, SciPost Phys.14, 001 (2023)

    2023

  11. [11]

    Nahum, P

    A. Nahum, P. Serna, J. T. Chalker, M. Ortu ˜no, and A. M. So- moza, Emergent SO(5) Symmetry at the N´eel to Valence-Bond- Solid Transition, Phys. Rev. Lett.115, 267203 (2015)

    2015

  12. [12]

    Nakayama and T

    Y. Nakayama and T. Ohtsuki, Necessary Condition for Emergent Symmetry from the Conformal Bootstrap, Phys. Rev. Lett.117, 131601 (2016)

    2016

  13. [13]

    Poland, S

    D. Poland, S. Rychkov, and A. Vichi, The conformal boot- strap: Theory, numerical techniques, and applications, Rev. Mod. Phys.91, 015002 (2019)

    2019

  14. [14]

    F. F. Assaad and I. F. Herbut, Pinning the Order: The Nature of Quantum Criticality in the Hubbard Model on Honeycomb Lattice, Phys. Rev. X3, 031010 (2013)

    2013

  15. [15]

    Parisen Toldin, M

    F. Parisen Toldin, M. Hohenadler, F. F. Assaad, and I. F. Her- but, Fermionic quantum criticality in honeycomb and𝜋-flux Hubbard models: Finite-size scaling of renormalization-group- invariant observables from quantum Monte Carlo, Phys. Rev. B 91, 165108 (2015)

    2015

  16. [16]

    Otsuka, S

    Y. Otsuka, S. Yunoki, and S. Sorella, Universal Quantum Crit- icality in the Metal-Insulator Transition of Two-Dimensional Interacting Dirac Electrons, Phys. Rev. X6, 011029 (2016)

    2016

  17. [17]

    I. F. Herbut, Interactions and Phase Transitions on Graphene’s Honeycomb Lattice, Phys. Rev. Lett.97, 146401 (2006)

    2006

  18. [18]

    I. F. Herbut, V. Juriˇci´c, and O. Vafek, Relativistic Mott criticality in graphene, Phys. Rev. B80, 075432 (2009)

    2009

  19. [19]

    Janssen and I

    L. Janssen and I. F. Herbut, Antiferromagnetic critical point on graphene’s honeycomb lattice: A functional renormalization group approach, Phys. Rev. B89, 205403 (2014)

    2014

  20. [20]

    Ladovrechis, S

    K. Ladovrechis, S. Ray, T. Meng, and L. Janssen, Gross-Neveu- Heisenberg criticality from 2+𝜖expansion, Phys. Rev. B107, 035151 (2023)

    2023

  21. [21]

    T. C. Lang and A. M. L¨auchli, Chiral Heisenberg Gross-Neveu- Yukawa criticality: Honeycomb versus SLAC fermions, Phys. Rev. B112, 245121 (2025)

    2025

  22. [22]

    F.-H. Wang, F. Sun, C. He, and X. Y. Xu, Resolving Quantum Criticality in the Honeycomb Hubbard Model, arXiv:2602.03656

    arXiv

  23. [23]

    Vafek, Interacting fermions on the honeycomb bilayer: From weak to strong coupling, Phys

    O. Vafek, Interacting fermions on the honeycomb bilayer: From weak to strong coupling, Phys. Rev. B82, 205106 (2010)

    2010

  24. [24]

    S. Ray, M. Vojta, and L. Janssen, Quantum critical behavior of two-dimensional Fermi systems with quadratic band touching, Phys. Rev. B98, 245128 (2018)

    2018

  25. [25]

    Pujari, T

    S. Pujari, T. C. Lang, G. Murthy, and R. K. Kaul, Interaction- Induced Dirac Fermions from Quadratic Band Touching in Bi- layer Graphene, Phys. Rev. Lett.117, 086404 (2016)

    2016

  26. [26]

    Z. C. Wei, C. Wu, Y. Li, S. Zhang, and T. Xiang, Majorana Positivity and the Fermion Sign Problem of Quantum Monte Carlo Simulations, Phys. Rev. Lett.116, 250601 (2016)

    2016

  27. [27]

    Sun and X

    F. Sun and X. Y. Xu, Boosting determinant quantum Monte Carlo with submatrix updates: Unveiling the phase diagram of the 3D Hubbard model, SciPost Phys.18, 055 (2025)

    2025

  28. [28]

    Z. H. Liu, M. Vojta, F. F. Assaad, and L. Janssen, Metallic and Deconfined Quantum Criticality in Dirac Systems, Phys. Rev. Lett.128, 087201 (2022)

    2022

  29. [29]

    Z. H. Liu, M. Vojta, F. F. Assaad, and L. Janssen, Critical properties of metallic and deconfined quantum phase transitions in Dirac systems, Phys. Rev. B110, 125123 (2024)

    2024

  30. [30]

    Campostrini, A

    M. Campostrini, A. Pelissetto, and E. Vicari, Finite-size scaling at quantum transitions, Phys. Rev. B89, 094516 (2014)

    2014

  31. [31]

    Harada, Bayesian inference in the scaling analysis of critical phenomena, Phys

    K. Harada, Bayesian inference in the scaling analysis of critical phenomena, Phys. Rev. E84, 056704 (2011)

    2011

  32. [32]

    T. C. Lang and A. M. L¨auchli, Quantum Monte Carlo Simulation of the Chiral Heisenberg Gross-Neveu-Yukawa Phase Transition with a Single Dirac Cone, Phys. Rev. Lett.123, 137602 (2019)

    2019

  33. [33]

    Hohenadler, Y

    M. Hohenadler, Y. Liu, T. Sato, Z. Wang, W. Guo, and F. F. As- saad, Thermodynamic and dynamical signatures of a quantum spin Hall insulator to superconductor transition, Phys. Rev. B 106, 024509 (2022)

    2022

  34. [34]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. 6 Mod. Phys.81, 109 (2009)

    2009

  35. [35]

    See Supplemental Material, which additionally includes Refs. [41, 42], for (i) technical details of the quantum Monte Carlo simulations, (ii) additional finite-size scaling analyses, (iii) a discussion of the temperature dependence of the charge susceptibility in the QBT regime, and (iv) additional finite- temperature Monte Carlo data for the charge susce...

  36. [36]

    F. F. Assaad, M. Bercx, F. Goth, A. G¨otz, J. S. Hofmann, E. Huff- man, Z. Liu, F. P. Toldin, J. S. E. Portela, and J. Schwab, The ALF (Algorithms for Lattice Fermions) project release 2.4. Doc- umentation for the auxiliary-field quantum Monte Carlo code, SciPost Phys. Codebases1-v2.4(2025)

    2025

  37. [37]

    K. S. D. Beach, Identifying the maximum entropy method as a special limit of stochastic analytic continuation, arXiv:cond- mat/0403055

    arXiv

  38. [38]

    G. Li, A. Luican, and E. Y. Andrei, Scanning Tunneling Spec- troscopy of Graphene on Graphite, Phys. Rev. Lett.102, 176804 (2009)

    2009

  39. [39]

    A. K. Geim and K. S. Novoselov, The rise of graphene, Nat. Mater.6, 183 (2007)

    2007

  40. [40]

    http://www.nhr-verein.de/en/our-partners

  41. [41]

    K. Sun, H. Yao, E. Fradkin, and S. A. Kivelson, Topological Insulators and Nematic Phases from Spontaneous Symmetry Breaking in 2D Fermi Systems with a Quadratic Band Crossing, Phys. Rev. Lett.103, 046811 (2009)

    2009

  42. [42]

    Emergent relativistic symmetry from interacting electrons on the honeycomb bilayer

    S. Ray, M. Vojta, and L. Janssen, Soluble fermionic quantum critical point in two dimensions, Phys. Rev. B102, 081112(R) (2020). Supplemental Material for “Emergent relativistic symmetry from interacting electrons on the honeycomb bilayer” Zi Hong Liu 1, 2 and Lukas Janssen 2 1State Key Laboratory of Surface Physics and Department of Physics, Fudan Univer...

    2020