Twisted Bilayer Graphene Magic Angle Calculation

Estimate the first magic angle of twisted bilayer graphene using the continuum-model relation between interlayer coupling, Fermi velocity, lattice constant, and the target dimensionless coupling parameter. The tool also calculates the moire wavelength, wavevector mismatch, and a comparison against any trial twist angle you enter.

Results

Expert Guide to Twisted Bilayer Graphene Magic Angle Calculation

Twisted bilayer graphene, often abbreviated TBG, is one of the most important model systems in modern condensed matter physics. It consists of two atomically thin graphene sheets placed on top of one another with a small rotational offset. That tiny twist changes the electronic structure dramatically. At ordinary angles the electronic bands resemble slightly perturbed graphene bands, but near a narrow window around roughly 1.1 degrees, the low-energy bands become exceptionally flat. Those flat bands suppress kinetic energy, amplify the importance of electron-electron interactions, and open the door to correlated insulating behavior, superconductivity, and a broad family of many-body quantum phases.

The phrase magic angle refers to the twist angle at which the low-energy bands become maximally flat in the continuum approximation first popularized by Bistritzer and MacDonald. Although experimental devices can differ in relaxation, strain, heterostrain, dielectric screening, and local tunneling amplitudes, a first-pass estimate of the magic angle is still extremely useful. It helps researchers decide which rotational alignment they are targeting during fabrication, how to interpret scanning probe and transport data, and how parameter choices such as interlayer tunneling strength influence the expected angle window.

What the calculator is actually computing

This calculator implements a standard continuum-model estimate. In this framework, the important dimensionless coupling parameter is

alpha = w / (hbar v_F k_theta)

where w is the effective interlayer coupling energy, v_F is the monolayer graphene Fermi velocity, and k_theta is the momentum separation between the two Dirac cones induced by twisting. That momentum mismatch is computed as

k_theta = 2K sin(theta / 2), with K = 4pi / 3a

where a is the graphene lattice constant. In the small-angle regime, the first magic angle is often estimated by setting alpha close to 0.605. Solving for theta gives an angle typically near 1 degree, depending on the assumed values of w and v_F. In other words, the magic angle is not a universal constant independent of the model. It emerges from a competition between interlayer hybridization and the monolayer kinetic scale.

The first magic angle commonly quoted in the literature is about 1.05 to 1.15 degrees, but the exact prediction shifts if you change interlayer coupling, band-structure renormalization, lattice relaxation, or tunneling anisotropy.

Why the magic angle matters physically

Flat bands mean electrons move sluggishly in momentum space. When the single-particle bandwidth narrows, interaction energy can dominate over kinetic energy. This is the regime in which strongly correlated phases become likely. In twisted bilayer graphene, correlated insulating states and superconducting behavior attracted global attention because they appear in a tunable and remarkably clean two-dimensional material platform. The same ideas now extend into a wider family of moire systems, including twisted transition-metal dichalcogenides and multilayer graphene-based heterostructures.

For TBG specifically, a magic-angle estimate is also useful because the moire length scale becomes large enough to be probed by microscopy and spectroscopy. The moire wavelength near 1.1 degrees is roughly a dozen nanometers. That enlarged superlattice creates a mini Brillouin zone and reconstructed bands with reduced bandwidth. Even a few tenths of a degree away from the optimal angle can noticeably change the moire period and therefore the low-energy physics.

Key parameters and realistic values

Several physical inputs control the estimate:

• Graphene lattice constant a: approximately 0.246 nm.
• Carbon-carbon bond length: about 0.142 nm, distinct from the lattice constant and often confused with it.
• Fermi velocity v_F: commonly around 0.8 x 10⁶ to 1.1 x 10⁶ m/s depending on model details and renormalization.
• Interlayer coupling w: often modeled in the neighborhood of 90 to 130 meV, though effective tunneling can be split into more than one parameter once relaxation effects are included.
• Magic-condition alpha: first-order estimate near 0.605 for the first magic angle in the isotropic continuum approximation.

Parameter	Typical value	Units	Why it matters
Graphene lattice constant a	0.246	nm	Sets the graphene K-point wavevector and therefore the twist-induced momentum mismatch.
Fermi velocity vF	0.8 to 1.1 x 10^6	m/s	Controls the monolayer kinetic energy scale.
Interlayer coupling w	90 to 130	meV	Controls how strongly the two graphene layers hybridize.
First magic-condition alpha	about 0.605	dimensionless	Approximate benchmark for maximal low-energy band flattening in the basic continuum model.
First magic angle	about 1.05 to 1.15	degrees	The angle region where flat-band physics is expected to emerge.

How to interpret the moire wavelength

The moire wavelength is a geometric consequence of the twist. A common expression is

L_m = a / [2 sin(theta / 2)]

for identical layers twisted by an angle theta. As theta decreases, the moire pattern becomes larger. Near 1 degree, this superlattice length becomes much larger than the atomic spacing, which is why continuum models work so well and why moire minibands can be treated as emergent low-energy electronic structures.

Twist angle	Approximate moire wavelength	Interpretation
2.0 degrees	about 7.0 nm	Significant moire reconstruction, but generally less flat-band enhancement than near the first magic angle.
1.5 degrees	about 9.4 nm	Intermediate regime with larger moire unit cell and stronger interlayer band reconstruction.
1.10 degrees	about 12.8 nm	Classic flat-band regime associated with correlated transport phenomena.
1.00 degrees	about 14.1 nm	Very large moire cell, sensitive to relaxation and local disorder.
0.80 degrees	about 17.6 nm	Even larger supercell, with stronger sensitivity to strain and inhomogeneity.

Step-by-step approach to calculating the magic angle

1. Choose a graphene lattice constant. For most practical calculations, use 0.246 nm.
2. Choose the monolayer Fermi velocity. A default of 1.0 x 10⁶ m/s is a common starting point.
3. Choose an effective interlayer coupling in meV. A value near 110 to 120 meV often reproduces a first magic angle around 1 degree in simple estimates.
4. Set the target alpha. For the first magic angle in the isotropic continuum model, 0.605 is widely used.
5. Convert the coupling from meV to joules and the lattice constant from nm to meters.
6. Compute K = 4pi / 3a.
7. Solve k_theta = w / (alpha hbar v_F).
8. Recover theta from k_theta = 2K sin(theta / 2). The calculator uses the exact inverse sine form, not just the small-angle approximation.
9. Compute the moire wavelength at the resulting angle to estimate the moire superlattice size.

Important limitations of any simple magic-angle calculator

A calculator like this one is best viewed as a high-quality first estimate, not a complete many-body simulation. Real TBG devices show several complications:

• Lattice relaxation: Atomic reconstruction changes local stacking areas and effectively modifies tunneling amplitudes.
• Tunneling anisotropy: Some advanced continuum models distinguish AA and AB tunneling terms rather than using one single coupling value.
• Heterostrain and twist-angle disorder: Spatial variation of angle or strain broadens or shifts the correlated regime.
• Dielectric environment: Encapsulation in hBN and gate screening can modify interaction scales.
• Interaction effects: The experimentally observed phase diagram is not determined by single-particle band flattening alone.

Even with those caveats, the first magic-angle estimate is extremely valuable because it captures the central kinematic balance in the problem: stronger interlayer coupling pushes the magic angle upward, while a larger Fermi velocity pushes it downward. This tradeoff explains why different papers can report slightly different preferred angles while still describing the same essential physics.

How experimentalists use this information

In fabrication, the twist angle is often targeted through tear-and-stack assembly. Researchers attempt to align flakes so the final rotational mismatch lands close to the desired angle. Small angle errors matter. A device intended for 1.10 degrees but fabricated at 1.30 degrees may still show moire effects, yet the band flatness and correlated behavior can be significantly altered. During characterization, the moire period can be cross-checked by scanning probe techniques, diffraction, or transport signatures tied to moire filling. A quick calculator helps translate among these different observables.

In theory, the same estimate guides parameter scans. If you vary w from 100 meV to 125 meV while holding the Fermi velocity fixed, the predicted first magic angle shifts upward. If instead you increase the Fermi velocity, the twist angle needed to reach the same alpha decreases. That sensitivity is why quoted values should always be tied to a specific model definition rather than treated as a universal exact number.

Authoritative scientific and educational resources

For readers who want to dig deeper into the underlying physics and experimental background, these sources are useful starting points:

• NIST overview of magic-angle graphene and superconductivity
• Moire bands in twisted double-layer graphene hosted through Princeton-related publication access
• MIT Physics explanation of correlated behavior in magic-angle graphene

Bottom line

Twisted bilayer graphene magic-angle calculation is the process of identifying the twist where interlayer hybridization and monolayer kinetic energy balance in a way that flattens the low-energy bands. In the simplest and most widely used estimate, you supply the graphene lattice constant, the Fermi velocity, the effective interlayer coupling, and a target alpha near 0.605. From those inputs, you obtain a predicted first magic angle and the corresponding moire wavelength. For many realistic parameter sets, the result falls close to 1.1 degrees, which is why this number has become central to modern moire materials research.

Scientific note: this calculator is intended for educational and preliminary design use. Precise device analysis may require full continuum-model band calculations including relaxation, anisotropic tunneling, strain, and interaction corrections.