Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Some effective operators for graphene monolayer superlattices, from variational perturbation theory

Louis Garrigue

TL;DR

The paper develops accurate effective operators for graphene monolayer superlattices by coupling variational basis constructions with perturbation theory in a two-scale framework. It builds a reduced two-scale space using derivatives of Dirac Bloch modes, derives Schur-reduced 2×2 matrix operators, and provides explicit parameterizations and symmetry-based matrix computations. Numerical simulations demonstrate improved band diagrams and eigenvectors over the conventional massless Dirac model, with careful control of spectral pollution via a plane-wave cutoff and scale parameter. The approach offers a scalable, systematically improvable effective model for low-energy Dirac fermions in periodically modulated graphene, with potential relevance to moiré-like and other multiscale graphene systems.

Abstract

Our goal is to provide precise effective operators for monolayer graphene at Fermi energy. We consider the microscopic potential created by a lattice, and add a macroscopic potential with the same periodicity but varying at a scale $\varepsilon^{-1} \in \mathbb{N}$, creating a superlattice. Our approach consists in coupling the variational approximation, perturbation theory together with a multiscale method. At the effective level the usual massless Dirac operator is replaced by other operators, and we provide simulations in the case of graphene.

Some effective operators for graphene monolayer superlattices, from variational perturbation theory

TL;DR

The paper develops accurate effective operators for graphene monolayer superlattices by coupling variational basis constructions with perturbation theory in a two-scale framework. It builds a reduced two-scale space using derivatives of Dirac Bloch modes, derives Schur-reduced 2×2 matrix operators, and provides explicit parameterizations and symmetry-based matrix computations. Numerical simulations demonstrate improved band diagrams and eigenvectors over the conventional massless Dirac model, with careful control of spectral pollution via a plane-wave cutoff and scale parameter. The approach offers a scalable, systematically improvable effective model for low-energy Dirac fermions in periodically modulated graphene, with potential relevance to moiré-like and other multiscale graphene systems.

Abstract

Our goal is to provide precise effective operators for monolayer graphene at Fermi energy. We consider the microscopic potential created by a lattice, and add a macroscopic potential with the same periodicity but varying at a scale , creating a superlattice. Our approach consists in coupling the variational approximation, perturbation theory together with a multiscale method. At the effective level the usual massless Dirac operator is replaced by other operators, and we provide simulations in the case of graphene.
Paper Structure (57 sections, 9 theorems, 177 equations, 8 figures)

This paper contains 57 sections, 9 theorems, 177 equations, 8 figures.

Table of Contents

  1. Introduction
  2. Introduction of the setting
  3. The lattice
  4. Scaling
  5. Symmetries of the potentials
  6. The scaled exact operator
  7. The microscopic exact problem
  8. The exact operator
  9. Exact eigenvectors when $V=0$
  10. The general effective operators
  11. Two-scales reduced basis
  12. Effective operator
  13. First simplification
  14. Particular effective models
  15. Preliminary information about the unscaled system
  16. ...and 42 more sections

Key Result

Proposition 3.1

Take $\alpha, \beta \in \mathbb{C}^M \otimes \mathcal{C}^{\infty}_{\rm per}(\Omega)$, and $V,A \in \mathcal{C}^{\infty}_{\rm per}(\Omega)$. For any $N \in \mathbb{N}$, there exists $C_N > 0$ such that for any $\varepsilon \in ]0,1[ \cap (1/\mathbb{N})$,

Figures (8)

  • Figure 1: Band diagram when $V=0$ and $\varepsilon=1$.
  • Figure 2: Letting $\nu$ vary at fixed $\varepsilon = \frac{1}{7}$, and $\mathcal{F} = \mathcal{F}_0$, $V = 0$.
  • Figure 3: Letting $\varepsilon$ decrease at fixed but large $\nu = 21$, and where only bands of the effective model are plotted. We chose $\mathcal{F}_{1}$, $\nu = 21$, $V = 0$.
  • Figure 4: Letting $\ell$ vary, with $\varepsilon = \frac{1}{7}$, $\nu = 2$ and $V = 0$.
  • Figure 5: Letting $\ell$ vary, with $\varepsilon = \frac{1}{7}$, $\nu = 2$ and a non-zero $V = V_{\text{ng}}$ (with $\lambda = 5$) which does not have honeycomb symmetry, presented in \ref{['eq:V_non_graphene']}.
  • ...and 3 more figures

Theorems & Definitions (16)

  • Proposition 3.1: Weak convergence
  • Remark 3.2: Undo the Bloch transform
  • Remark 3.3: Effective operators are defined for any $\varepsilon >0$
  • Proposition 4.1: Effective matrices for $\mathcal{F}_{1,k}$
  • Proposition 4.2: Effective matrices for $\mathcal{F}_{1}$
  • Proposition 4.3: Effective matrices for $\mathcal{F}^3$
  • Lemma 6.1
  • proof
  • Lemma 7.1: One derivative
  • proof
  • ...and 6 more