Pith. sign in

REVIEW 3 major objections 4 minor 43 references

The emergence of one-dimensional channels in marginal-angle twisted bilayer graphene

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read Applying an interlayer electric bias to relaxed twisted bilayer graphene at marginal angles frees a stable network of quasi-one-dimensional channel states, leaving a bias-independent low-energy spectrum.

desk verdict A useful extension of the continuum model to relaxed marginal-angle TBG, with credible bias-liberated channel states; the missing harmonic-convergence check at 0.221° is the main thing standing between this and a fully convincing claim. read the letter →

arxiv 1908.05941 v1 pith:WEYYKLOU submitted 2019-08-16 cond-mat.str-el cond-mat.mes-hall

Niels R. Walet , Francisco Guinea This is my paper · ORCID
classification cond-mat.str-elcond-mat.mes-hall
keywords twistedbilayergraphenemarginaltwistanglelatticerelaxationone-dimensionalchannelsinterlayerbiascontinuummodelchiraledgestatesmoiré
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

This paper argues that in relaxed twisted bilayer graphene at marginal twist angles near 0.2 degrees, an applied interlayer electric bias separates the electronic spectrum into two parts: states tied to AB-stacked regions, which shift in energy with the bias, and a stable network of quasi-one-dimensional states confined to the relaxation-induced channels between AB and BA regions, which stay put. The authors show that as the bias grows, the low-energy bands and density of states converge to a bias-independent limit, which they read as the signature of a standalone channel network. They establish this using a generalized continuum model that includes lattice relaxation and a large number of interlayer harmonics, making small-angle calculations tractable without a costly tight-binding treatment of hundreds of thousands of atoms. If correct, the work identifies a concrete experimental route—applying a bias to marginally twisted bilayers—for isolating and studying one-dimensional chiral channel states.

What carries the argument

The machinery is a generalized continuum model for twisted bilayer graphene that starts from a relaxed atomic structure (computed with the LCBOPI+KC potential), projects the interlayer hopping onto Bloch states, and includes 75 reciprocal-lattice harmonics rather than the usual three of the Bistritzer-MacDonald model. This makes small-angle calculations feasible with accuracy comparable to tight binding at larger angles, despite over 268,000 atoms in the moiré unit cell at 0.221 degrees. The paper uses an exponential Koster-Slater interlayer hopping, scaled so the first magic angle sits near 1.05 degrees, and applies an interlayer bias as a probe to distinguish bias-sensitive AB-region states from bias-insensitive channel states.

What would settle it

A direct tight-binding or other numerically exact calculation for a relaxed 0.2-degree bilayer that shows the low-energy spectrum remaining strongly bias-dependent at high bias—or a scanning probe experiment at 0.2 degrees with bias that finds bias-sensitive states coexisting with the channel network—would contradict the claim that the low-energy limit is bias independent and channel-only.

Watch

Extended reading notes

Core claim

The central discovery is that lattice relaxation at marginal twist angles creates a network of narrow channels separating AB and BA stacked domains, and that these channels support quasi-one-dimensional states that are insensitive to an interlayer bias. In the authors' continuum calculations for a relaxed bilayer at 0.221 degrees, the low-energy bands and density of states approach a bias-independent limit as the bias is increased from 0 to 200 meV, whereas unrelaxed layers show little spectral change and larger angles show only weak sensitivity. The bias liberates the channel states by gapping the AB-region states, leaving a stable network whose wave functions concentrate on the channels with chiral interlayer beats. The paper also finds that a single pair of nearly straight, linearly dispersing bands persists, and that crossings occur at K, K' and near M, differing from simpler one-channel models. The authors conclude that the occurrence of channel states is a robust consequence of relaxation in bilayer graphene, with details depending on the theoretical model and experimental environment.

Load-bearing premise

The calculations assume that the relaxed atomic structure from the LCBOPI+KC potential, together with the 75-vector harmonic truncation of interlayer coupling, correctly reproduces the low-energy electronic structure at angles near 0.2 degrees, without an independent check against a full tight-binding calculation at those angles.

Editorial extensions

If this is right

  • At angles below about 0.5 degrees, applying an interlayer bias of order 100 meV should drive the system into a regime where only channel states remain near the Fermi energy, enabling transport through a quasi-1D network.
  • The generalized continuum model with many harmonics extends quantitative band-structure calculations to twist angles as small as 0.1 degrees, where the moiré unit cell contains about a million atoms, at no significant increase in computational cost.
  • The channel states show a chiral structure with interlayer beats, so their detection in scanning tunneling microscopy would provide a direct signature of the relaxation-induced channels.
  • Only one pair of nearly straight, linearly dispersing bands survives, with crossings at K, K' and near M, suggesting that simplified one-channel models miss important scattering and energy-dependent velocity effects.

Reading between the lines

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

  • The bias-independent limit implies that the channel network could serve as a platform for studying one-dimensional conductance quantization and interferometry, since the AB-region states can be gapped away without destroying the channels.
  • If the channel states survive at finite temperature and moderate disorder, the marginal-angle regime with bias might exhibit transport dominated by the channel network's topology, potentially mimicking features of electrical networks in other moiré systems.
  • The method's reliance on a relaxed structure means that the channel network's geometry—and thus the device behavior—could be tunable by strain or substrate encapsulation, a testable prediction beyond the free-standing planar case.
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

3 major / 4 minor

Summary. The paper generalizes the continuum model of twisted bilayer graphene to include lattice relaxation and a large number of interlayer harmonics, and applies it to marginal twist angles near 0.2 degrees. The authors compute relaxed atomic structures using the LCBOPI+KC potential, project the interlayer hopping onto a 75-G-vector basis, and study the effect of an interlayer bias on the low-energy spectrum. Their central claim is that at small angles an applied bias gaps out states localized in AB-stacked regions, liberating a set of bias-insensitive quasi-one-dimensional states confined to the relaxation-induced channels separating AB and BA domains. They compare their results with the Efimkin-MacDonald helical network model and the Zhang-MacDonald-Mele domain-wall model, and also provide wave-function plots showing channel localization and interlayer beats. The main conclusion is that the low-energy bands and density of states tend to a bias-independent limit, consistent with a situation in which only channel states remain.

Significance. If the central claim holds, the paper provides a concrete and testable prediction: at marginal twist angles around 0.2 degrees, an interlayer bias can isolate a robust network of quasi-one-dimensional channel states that are insensitive to the bias and confined to AB/BA domain walls. This is directly relevant to recent experiments on marginally twisted bilayers (Yoo et al., Xu et al., Rickhaus et al.) and goes beyond earlier model studies by including realistic relaxation and a more complete interlayer-coupling expansion. Strengths of the paper include the extensive numerical spectra, density-of-states plots, and wave-function visualizations, and the fact that the central bias-liberation prediction is not used to fit any constant, giving it independent predictive content. The method extends the authors' prior continuum-projection framework to smaller angles, and the robustness of the qualitative behavior is asserted across alternative relaxation and hopping models, though the supporting data for those alternatives is not shown.

major comments (3)
  1. [Methods and Sec. I (Fig. 3b, Fig. S5)]
  2. [Sec. I, Eq. (1)]
  3. [Sec. I, Fig. 5 and Conclusions]
minor comments (4)
  1. [Sec. I, after Eq. (1)]
  2. [Fig. 2a and Fig. S5]
  3. [Sec. I, first Results paragraph]
  4. [Sec. I, Fig. 4 discussion]

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the channel-state prediction is not fitted to the bias-response data it claims to explain.

full rationale

The central claim—bias-liberated quasi-1D channel states at 0.221°—is obtained from a Hamiltonian constructed from relaxed geometries, projected interlayer couplings, and a 75-harmonic basis. The only fitted constant, the 1.3 enhancement of t(r12) in Eq. (1), is calibrated to place the magic angle at 1.05°, an external benchmark, and is not tuned to the 0.221° channel-state data. The conclusion that low-energy bands and the density of states tend to a bias-independent limit is supported by computed DOS curves and by the wavefunction localization shown in Fig. 6, not by reinserting the claim as an input. Although the colored-band identification of 'channel states' in Fig. 5 is guided by bias insensitivity, the paper independently demonstrates spatial confinement to relaxation-induced channels, so the interpretation is not purely definitional. Reliance on the authors' prior continuum-projection method [33] is a standard method citation with independent validation in that earlier work; no uniqueness theorem is invoked to forbid alternatives, and no prediction reduces to a fit by construction. The absence of an explicit convergence check at 0.221° is a correctness or robustness concern, not a circularity. The derivation is therefore self-contained with respect to its central claim.

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

The central claim rests on the calibrated interlayer hopping and the assumed relaxation potential, plus the validity of the continuum projection. No new physical entities are introduced.

free parameters (1)
  • interlayer hopping enhancement factor = 1.3
    The Koster-Slater interlayer hopping in Eq. (1) is multiplied by 1.3 relative to previous work so that the first magic angle occurs at 1.05 degrees for the relaxed LKC lattice, calibrating the model to a known result.
assumptions (5)
  • domain assumption A finite projection of Bloch states onto superlattice harmonics reproduces the tight-binding low-energy spectrum.
    Inherited from ref. [33]; required for the entire method, and not re-validated at the smallest angles in this paper.
  • domain assumption In-plane hopping parameters are unchanged by relaxation since bond stretching is smaller than 1 part in 10^4.
    Stated in the introduction; allows use of the unperturbed in-layer bands.
  • domain assumption The interlayer distance is fixed at d = 3.46 Å with no corrugation.
    Stated in Results; real samples may have distance variations that affect channel confinement.
  • domain assumption Intervalley (K-K') coupling is negligible at the small angles considered.
    The authors note large-angle results are incorrect due to neglect of intervalley coupling; it is assumed small at marginal angles.
  • domain assumption The LCBOPI+KC potential describes the relaxed channel structure.
    Used to generate relaxed positions; no ab initio cross-check is provided in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The emergence of one-dimensional channels in marginal-angle twisted bilayer graphene." pith.science (2026). https://pith.science/paper/WEYYKLOU

@misc{pith2026190805941,
  author       = {Pith},
  title        = {Pith review of: The emergence of one-dimensional channels in marginal-angle twisted bilayer graphene},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/WEYYKLOU}},
  note         = {Machine review of arXiv:1908.05941}
}
read the original abstract

We generalize the continuum model for Moir\'e structures made from twisted graphene layers, in order to include lattice relaxation and the formation of channels at very small (marginal) twist angles. We show that a precise description of the electronic structure at such small angles can be achieved by i) calculating first the relaxed atomic structure, ii) projecting the interlayer electronic hopping parameters using a suitable basis of Bloch states, and iii) increasing the number of harmonics in the continuum approximation to interlayer hopping. The results show a complex structure of quasi one dimensional states when a finite bias is applied.

Figures

Figures reproduced from arXiv: 1908.05941 by the authors.

Figure 1
Figure 1. A representation of the relaxation of a bilayer for a [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. a) The alignment measure from Ref. [33] for a flat relaxed bilayer. The green areas are [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 5
Figure 5. Bands for the continuum projected models with a [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figures from the paper (2 more)
Figure 4
Figure 4. Figure 4: a) Density of states for the continuum projected [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 6
Figure 6. Figure 6: The probability density for states near zero energy [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 33 canonical work pages

  1. [1]

    Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, R. C. Ashoori, and P. Jarillo- Herrero, Correlated insulator behaviour at half-filling in magic-angle graphene superlattices, Nature 10.1038/na- ture26154 (2018)

    doi:10.1038/na- 2018

  2. [2]

    Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniuchi, E. Kaxiras, and P. Jarillo-Herrero, Unconventional super- conductivity in magic-angle graphene superlattices, Na- ture 10.1038/nature26160 (2018)

    doi:10.1038/nature26160 2018

  3. [3]

    H. Yoo, R. Engelke, S. Carr, S. Fang, K. Zhang, P. Cazeaux, S. H. Sung, R. Hovden, A. W. Tsen, T. Taniguchi, K. Watanabe, G.-C. Yi, M. Kim, M. Luskin, E. B. Tadmor, E. Kaxiras, and P. Kim, Atomic and electronic reconstruction at the van der Waals interface in twisted bilayer graphene, Nature Ma- terials 18, 448 (2019)

    work page 2019

  4. [4]

    S. G. Xu, A. I. Berdyugin, P. Kumaravadivel, F. Guinea, R. K. Kumar, D. A. Bandurin, S. V. Morozov, W. Kuang, B. Tsim, S. Liu, J. H. Edgar, I. V. Grigorieva, V. I. Fal’ko, M. Kim, and A. K. Geim, Giant oscillations in a tri- angular network of one-dimensional states in marginally twisted graphene, arXiv:1905.12984 [cond-mat] (2019), arXiv: 1905.12984

    arXiv 2019

  5. [5]

    Chung, R

    T.-F. Chung, R. He, T.-L. Wu, and Y. P. Chen, Optical Phonons in Twisted Bilayer Graphene with Gate-Induced Asymmetric Doping, Nano Lett. 15, 1203 (2015)

    work page 2015

  6. [6]

    L. A. Gonzalez-Arraga, J. L. Lado, F. Guinea, and P. San-Jose, Electrically controllable magnetism in twisted bilayer graphene, Phys. Rev. Lett. 119, 107201 (2017)

    work page 2017

  7. [7]

    Martin, Y

    I. Martin, Y. M. Blanter, and A. F. Morpurgo, Topolog- ical Confinement in Bilayer Graphene, Phys. Rev. Lett. 100, 036804 (2008)

    work page 2008

  8. [8]

    J. Jung, F. Zhang, Z. Qiao, and A. H. MacDonald, Valley- Hall kink and edge states in multilayer graphene, Phys. Rev. B 84, 075418 (2011)

    work page 2011

Show all 43 references
  1. [9]

    Zarenia, J

    M. Zarenia, J. M. Pereira, G. A. Farias, and F. M. Peeters, Chiral states in bilayer graphene: Magnetic field dependence and gap opening, Phys. Rev. B 84, 125451 (2011)

    2011

  2. [10]

    Vaezi, Y

    A. Vaezi, Y. Liang, D. H. Ngai, L. Yang, and E.-A. Kim, Topological Edge States at a Tilt Boundary in Gated Multilayer Graphene, Phys. Rev. X 3, 021018 (2013)

    2013

  3. [11]

    Zhang, A

    F. Zhang, A. H. MacDonald, and E. J. Mele, Valley Chern numbers and boundary modes in gapped bilayer graphene, PNAS 110, 10546 (2013)

    2013

  4. [12]

    X. Li, F. Zhang, Q. Niu, and A. MacDonald, Spontaneous Layer-Pseudospin Domain Walls in Bilayer Graphene, Phys. Rev. Lett. 113, 116803 (2014)

    2014

  5. [13]

    San-Jose, R

    P. San-Jose, R. V. Gorbachev, A. K. Geim, K. S. Novoselov, and F. Guinea, Stacking boundaries and transport in bilayer graphene, Nano Lett. 14, 2052 (2014)

    2014

  6. [14]

    Yeh, Y.-C

    C.-H. Yeh, Y.-C. Lin, Y.-C. Chen, C.-C. Lu, Z. Liu, K. Suenaga, and P.-W. Chiu, Gating Electron–Hole Asymmetry in Twisted Bilayer Graphene, ACS Nano 8, 6962 (2014)

    2014

  7. [15]

    D. A. Cosma and V. I. Fal’ko, Trigonal warping effect on velocity and transverse confinement length of topologi- cally confined states in bilayer graphene, Phys. Rev. B 92, 165412 (2015)

    2015

  8. [16]

    L. Ju, Z. Shi, N. Nair, Y. Lv, C. Jin, J. Velasco Jr, C. Ojeda-Aristizabal, H. A. Bechtel, M. C. Martin, A. Zettl, J. Analytis, and F. Wang, Topological valley transport at bilayer graphene domain walls, Nature 520, 650 (2015)

    2015

  9. [17]

    M. Pelc, W. Jask´ olski, A. Ayuela, and L. Chico, Topo- logically confined states at corrugations of gated bilayer graphene, Phys. Rev. B 92, 085433 (2015)

    2015

  10. [18]

    L.-J. Yin, H. Jiang, J.-B. Qiao, and L. He, Direct imaging of topological edge states at a bilayer graphene domain wall, Nature Communications 7, 11760 (2016). 7

    2016

  11. [19]

    C. Lee, G. Kim, J. Jung, and H. Min, Zero-line modes at stacking faulted domain walls in multilayer graphene, Phys. Rev. B 94, 125438 (2016)

    2016

  12. [20]

    J. Li, K. Wang, K. J. McFaul, Z. Zern, Y. Ren, K. Watan- abe, T. Taniguchi, Z. Qiao, and J. Zhu, Gate-controlled topological conducting channels in bilayer graphene, Na- ture Nanotechnology 11, 1060 (2016)

    2016

  13. [21]

    T. L. M. Lane, M. An delkovi´ c, J. R. Wallbank, L. Co- vaci, F. M. Peeters, and V. I. Fal’ko, Ballistic electron channels including weakly protected topological states in delaminated bilayer graphene, Phys. Rev. B 97, 045301 (2018)

    2018

  14. [22]

    Rickhaus, J

    P. Rickhaus, J. Wallbank, S. Slizovskiy, R. Pisoni, H. Overweg, Y. Lee, M. Eich, M.-H. Liu, K. Watan- abe, T. Taniguchi, T. Ihn, and K. Ensslin, Transport through a network of topological channels in twisted bi- layer graphene, Nano Lett. 18, 6725 (2018)

    2018

  15. [23]

    Cheng, H

    S.-G. Cheng, H. Liu, H. Jiang, Q.-F. Sun, and X. Xie, Manipulation and Characterization of the Valley- Polarized Topological Kink States in Graphene-Based In- terferometers, Phys. Rev. Lett. 121, 156801 (2018)

    2018

  16. [24]

    Y. Sun, H. Zhao, Z.-M. Yu, and H. Pan, Valley current and spin-valley filter in topological domain wall, Journal of Applied Physics 125, 123904 (2019)

    2019

  17. [25]

    Y. Yang, Z. Jia, Y. Wu, Z.-H. Hang, H. Jiang, and X. C. Xie, Gapped Topological Kink States and Topological Corner States in Graphene, arXiv:1903.01816 [cond-mat] (2019), arXiv: 1903.01816

    2019 arXiv

  18. [26]

    Wu, C.-M

    X.-C. Wu, C.-M. Jian, and C. Xu, Coupled-wire descrip- tion of the correlated physics in twisted bilayer graphene, Phys. Rev. B 99, 161405 (2019)

    2019

  19. [27]

    J. Liu, J. Liu, and X. Dai, Pseudo landau level representa- tion of twisted bilayer graphene: Band topology and im- plications on the correlated insulating phase, Phys. Rev. B 99, 155415 (2019)

    2019

  20. [28]

    X.-C. Wu, A. Keselman, C.-M. Jian, K. A. Pawlak, and C. Xu, Ferromagnetism and Spin-Valley liquid states in Moir´ e Correlated Insulators, arXiv:1905.00033 [cond- mat] (2019), arXiv: 1905.00033

    2019 arXiv

  21. [29]

    San-Jose and E

    P. San-Jose and E. Prada, Helical networks in twisted bilayer graphene under interlayer bias, Phys. Rev. B 88, 121408 (2013)

    2013

  22. [30]

    D. K. Efimkin and A. H. MacDonald, Helical network model for twisted bilayer graphene, Phys. Rev. B 98, 035404 (2018)

    2018

  23. [31]

    T. Hou, Y. Ren, Y. Quan, J. Jung, W. Ren, and Z. Qiao, Current partition and conducting topological networks in twisted graphene bilayer, 1904.12826 (2019), arXiv:1904.12826 [cond-mat]

    2019 arXiv

  24. [32]

    N. N. T. Nam and M. Koshino, Lattice relaxation and en- ergy band modulation in twisted bilayer graphene, Phys. Rev. B 96, 075311 (2017)

    2017

  25. [33]

    Guinea and N

    F. Guinea and N. R. Walet, Continuum models for twisted bilayer graphene: Effect of lattice deformation and hopping parameters, Phys. Rev. B99, 205134 (2019)

    2019

  26. [34]

    Koshino and Y.-W

    M. Koshino and Y.-W. Son, Moir´ e phonons in twisted bilayer graphene, Phys. Rev. B 100, 075416 (2019)

    2019

  27. [35]

    Ochoa, Moir´ e-pattern fluctuations and electron- phason coupling in twisted bilayer graphene, arXiv:1905.10850 [cond-mat] (2019), arXiv:1905.10850

    H. Ochoa, Moir´ e-pattern fluctuations and electron- phason coupling in twisted bilayer graphene, arXiv:1905.10850 [cond-mat] (2019), arXiv:1905.10850

    2019 arXiv

  28. [36]

    Bistritzer and A

    R. Bistritzer and A. H. MacDonald, Moir´ e bands in twisted double-layer graphene, PNAS 108, 12233 (2011)

    2011

  29. [37]

    S. Carr, S. Fang, Z. Zhu, and E. Kaxiras, Exact con- tinuum model for low-energy electronic states of twisted bilayer graphene, Phys. Rev. Research 1, 013001 (2019)

    2019

  30. [38]

    S. Fang, S. Carr, Z. Zhu, D. Massatt, and E. Kaxi- ras, Angle-dependent ab initio low-energy hamiltonians for a relaxed twisted bilayer graphene heterostructure, arXiv:1908.00058 [cond-mat] (2019), arXiv:1908.00058

    2019 arXiv

  31. [39]

    J. M. B. Lopes dos Santos, N. M. R. Peres, and A. H. Castro Neto, Continuum model of the twisted graphene bilayer, Phys. Rev. B 86, 155449 (2012)

    2012

  32. [40]

    San-Jose, J

    P. San-Jose, J. Gonz´ alez, and F. Guinea, Non-abelian gauge potentials in graphene bilayers, Phys. Rev. Lett. 108, 216802 (2012)

    2012

  33. [41]

    8 Appendix A: Models of chiral bands

    In order to be able to overlay the probabilities, which are each defined on different lattice points, we have interpo- lated separately the wave functions for each of the four sublattices (top A and top B, bottom A and B) before adding the relevant wave functions in quadrature. ...

  34. [42]

    Efimkin-McDonald model In the Efimkin-McDonald model [30], the spectrum is (apart from a constant shift in energy) determined by the ballistic propagation of electrons along a channel, with scattering in all directions at the AA joints. Apart from some phases that only lead to a...

  35. [43]

    It is essentially a contin- uum model of aligned bilayer graphene, H =ντ0qxσx +τ0qyσy + 1 2γ (σxτx −µσyτy) + ∆σ0τz, (S1) expanded about the K (K′) valley for µ = 1 (µ = −1)

    Zhang-McDonald-Mele model The ZMM model [11] was instrumental in showing the fact that we have two chiral states at the interface be- tweenAB and BA alignment. It is essentially a contin- uum model of aligned bilayer graphene, H =ντ0qxσx +τ0qyσy + 1 2γ (σxτx −µσyτy) + ∆σ0τz, (...

Pith tools

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