When neutral graphene is placed nearly in rotational alignment on a hexagonal boron nitride (hBN) substrate, a long-wavelength moiré superlattice forms and an energy gap of a few to a few tens of meV often opens at the charge-neutrality (Dirac) point, but its presence and magnitude vary widely between samples. The theory of Jung, DaSilva, MacDonald and Adam (Nature Communications 6, 6308 (2015)) attributed the gap to the interplay of moiré-scale atomic relaxation, which generates a spatially averaged sublattice-staggering (mass) term, and nonlocal electron-electron exchange, estimated with a static effective dielectric constant. That paper explicitly acknowledged five unresolved, quantitatively important limitations: the dependence of relaxation and of the gap on BN substrate thickness and encapsulation; whether the macroscopic lattice mismatch and twist angle are set by kinetics or thermodynamics; the crudeness of the static-screening Hartree-Fock treatment; the neglect of higher moiré reciprocal-lattice shells and van der Waals tail corrections; and sensitivity to the parametrizations of the interaction kernels and elastic constants. This problem asks for the completion of that theory: a quantitative, experiment-comparable account of the gap and its sample-to-sample variability that resolves these five items, most importantly the interaction-driven enhancement on realistic relaxed structures and the substrate-thickness/encapsulation dependence.

Background

Graphene is a single layer of carbon atoms on a honeycomb lattice with two inequivalent sublattices, conventionally labeled A and B. Because the two sublattices are related by symmetry, isolated graphene has massless Dirac-like low-energy electrons and no band gap. A term in the Hamiltonian that gives different site energies to the A and B sublattices (a 'sublattice mass' or 'staggered potential') breaks this symmetry and opens a gap at the Dirac point, the energy where conduction and valence bands touch in neutral graphene. Hexagonal boron nitride (hBN) is an insulator with the same honeycomb structure but about 1.8% larger lattice constant, so stacking graphene on hBN with a small relative twist angle θ\theta produces a moiré superlattice: a long-period triangular pattern of locally varying stacking registries (AA, AB, BA), with period set by the relative lattice-constant mismatch ε\varepsilon and θ\theta. Density-functional theory for a perfectly commensurate, aligned stacking predicts gaps of order 50 meV, but any incommensurability strongly suppresses the electronic coupling, so rigid-lattice continuum models predict nearly vanishing gaps, while experiments on nearly aligned devices report gaps anywhere from a few meV to roughly 30-40 meV (from thermal-activation transport) or 14-17 meV (from tunneling and photocurrent spectroscopy). The 2015 paper of Jung, DaSilva, MacDonald and Adam resolved this tension in three steps. First, it computed moiré-scale atomic relaxation by minimizing an elastic-plus-interlayer energy: in-plane and vertical displacements of the carbon and boron/nitrogen atoms, with the interlayer stacking energy parametrized from exact-exchange-plus-random-phase-approximation (EXX+RPA) density-functional energy curves and the elasticity described by Lamé constants (the two standard elastic moduli of an isotropic two-dimensional continuum); this is essentially a two-dimensional Frenkel-Kontorova problem (the classic model of a chain or sheet of atoms in a periodic substrate potential, which develops commensurate domains separated by dislocation-like domain walls when the natural periods mismatch). Second, it fed the relaxed local stacking field into a continuum moiré Hamiltonian for graphene's low-energy pi electrons. Third, it added electron-electron interactions at the self-consistent Hartree-Fock level (a mean-field treatment of the Coulomb interaction, with exchange included nonlocally), using an effective static dielectric constant bracketed between 2.5 and 4 to model screening. The resulting hierarchy is the key quantitative landscape of the field: a rigid-lattice gap of only about 1 meV, roughly 7 meV when relaxation of graphene and the top BN layer is allowed, and roughly 20 meV with exchange enhancement, in reasonable agreement with experiment. The paper treats the moiré period parameters ε\varepsilon and θ\theta as experimentally fixed observables (arguing they are set by sample-transfer kinetics rather than by thermodynamic equilibration, since changing them requires collective motion of many atoms), truncates its interaction-energy expansion at the first shell of moiré reciprocal-lattice vectors, and neglects van der Waals tail corrections from the bulk substrate. Its Discussion and Methods sections explicitly flag the five limitations listed in the abstract, which this problem takes as its scope.

Problem Statement

The problem is to complete the microscopic, quantitative theory of the charge-neutrality (primary Dirac-point) band gap of nearly aligned graphene on hexagonal boron nitride, by resolving the five limitations that Jung, DaSilva, MacDonald and Adam (Nature Communications 6, 6308 (2015)) explicitly attached to their relaxation-plus-exchange model, under the single objective of predicting the observed gap magnitudes and their sample-to-sample variability.

Concretely: (i) determine how the moiré-scale atomic relaxation of the carbon and boron/nitrogen atoms, the resulting spatially averaged sublattice mass term, and hence the primary gap depend on the thickness of the BN substrate (i.e. on the stiffness or relaxation of remote BN layers) and on the presence of additional encapsulating layers; (ii) establish whether the macroscopic parameters — the relative lattice-constant difference and the twist angle — are locked in by transfer/growth kinetics or by thermodynamic equilibration, and what values or distributions of them experiments should produce; (iii) compute the enhancement of the primary gap by electron-electron interactions beyond the static effective dielectric constant (2.5 to 4) self-consistent Hartree-Fock treatment used in 2015 — including the effect of dynamical screening and of correlation beyond Hartree-Fock — using realistic relaxed moiré Hamiltonian parameters, and state how much of the experimentally observed gap (roughly 14-17 meV in tunneling and photocurrent spectroscopy, up to 30-40 meV in thermal-activation transport) is single-particle relaxation physics versus many-body enhancement; (iv) assess the quantitative importance, in the parameter regimes where the moiré period is short or local strains are large, of higher moiré reciprocal-lattice shells beyond the first-shell expansion and of longer-range van der Waals tail corrections to the stacking energy that were neglected in 2015; and (v) quantify the sensitivity of the predicted strain fields, averaged mass term, and gap to the parametrization of the EXX+RPA-derived interaction kernels and to the uncertainties in the elastic (Lamé) constants.

Items (i), (iii), (iv), and (v) jointly determine the magnitude of the same object — the primary gap and its variation across devices — and are therefore one problem; item (ii) concerns the physical origin of the input parameters of that same theory and is included as the source paper framed it. The task therefore spans all five items (i)-(v) and asks for a quantitative account, comparable with the available spectroscopic and transport gap measurements, that identifies which mechanisms set the gap and its variability and determines whether the twist angle and lattice-constant mismatch are set by transfer/growth kinetics or by thermodynamic equilibration.

References

  1. Jeil Jung, Ashley M. DaSilva, Allan H. MacDonald, Shaffique Adam, 'Origin of band gaps in graphene on hexagonal boron nitride', Nature Communications 6, 6308 (2015), DOI 10.1038/ncomms7308, https://doi.org/10.1038/ncomms7308
  2. Jeil Jung, Evan Laksono, Ashley M. DaSilva, Allan H. MacDonald, Marcin Mucha-Kruczyński, Shaffique Adam, 'Moiré band model and band gaps of graphene on hexagonal boron nitride', Physical Review B 96, 085442 (2017), DOI 10.1103/PhysRevB.96.085442, arXiv:1706.06016, https://doi.org/10.1103/PhysRevB.96.085442
  3. G. J. Slotman, M. M. van Wijk, Pei-Liang Zhao, A. Fasolino, M. I. Katsnelson, Shengjun Yuan, 'Effect of Structural Relaxation on the Electronic Structure of Graphene on Hexagonal Boron Nitride', Physical Review Letters 115, 186801 (2015), DOI 10.1103/PhysRevLett.115.186801, https://doi.org/10.1103/PhysRevLett.115.186801
  4. Hakseong Kim, Nicolas Leconte, Bheema L. Chittari, Kenji Watanabe, Takashi Taniguchi, Allan H. MacDonald, Jeil Jung, Suyong Jung, 'Accurate Gap Determination in Monolayer and Bilayer Graphene/h-BN Moiré Superlattices', Nano Letters 18, 7732-7741 (2018), DOI 10.1021/acs.nanolett.8b03423, arXiv:1808.06633, https://doi.org/10.1021/acs.nanolett.8b03423
  5. Tianyi Han, Jixiang Yang, Qihang Zhang, Lei Wang, Kenji Watanabe, Takashi Taniguchi, Paul L. McEuen, Long Ju, 'Accurate Measurement of the Gap of Graphene/hBN Moiré Superlattice through Photocurrent Spectroscopy', Physical Review Letters 126, 146402 (2021), DOI 10.1103/PhysRevLett.126.146402, arXiv:2103.08120, https://doi.org/10.1103/PhysRevLett.126.146402
  6. Jiaqi An, Nicolas Leconte, Srivani Javvaji, Youngju Park, Jeil Jung, 'Atomistic substrate relaxation effects in the band gaps of graphene on hexagonal boron nitride', arXiv:2602.18893 (2026), https://arxiv.org/abs/2602.18893
  7. C. R. Woods, L. Britnell, A. Eckmann, R. S. Ma, J. C. Lu, H. M. Guo, X. Lin, G. L. Yu, Y. Cao, R. V. Gorbachev, A. V. Kretinin, J. Park, L. A. Ponomarenko, M. I. Katsnelson, Yu. N. Gornostyrev, K. Watanabe, T. Taniguchi, C. Casiraghi, H.-J. Gao, A. K. Geim, K. S. Novoselov, 'Commensurate-incommensurate transition in graphene on hexagonal boron nitride', Nature Physics 10, 451-456 (2014), DOI 10.1038/nphys2954, https://doi.org/10.1038/nphys2954
  8. Nicolas Leconte, Jeil Jung, 'Commensurate and incommensurate double moiré interference in graphene encapsulated by hexagonal boron nitride', 2D Materials 7, 031005 (2020), DOI 10.1088/2053-1583/ab891a, https://doi.org/10.1088/2053-1583/ab891a
  9. M. Anđelković, S. P. Milovanović, L. Covaci, F. M. Peeters, 'Double moiré with a twist: super-moiré in encapsulated graphene', arXiv:1910.00345 (2019), https://arxiv.org/abs/1910.00345
  10. Jonah Herzog-Arbeitman, Yuzhi Wang, Jiaxuan Liu, Pok Man Tam, Ziyue Qi, Yujin Jia, Dmitri K. Efetov, Oskar Vafek, Nicolas Regnault, Hongming Weng, Quansheng Wu, B. Andrei Bernevig, Jiabin Yu, 'Moiré fractional Chern insulators. II. First-principles calculations and continuum models of rhombohedral graphene superlattices', Physical Review B 109, 205122 (2024), DOI 10.1103/PhysRevB.109.205122, arXiv:2311.12920, https://doi.org/10.1103/PhysRevB.109.205122
  11. Luca Nashabeh, Héctor Ochoa, 'Lattice Relaxation Flattens Chern Bands in Rhombohedral Graphene Stacks', arXiv:2605.16218 (2026), https://arxiv.org/abs/2605.16218