Articles
https://doi.org/10.1038/s41567-019-0606-5
Electronic correlations in twisted bilayer graphene
near the magic angle
Youngjoon Choi
1,2,3
, Jeannette Kemmer
1,2
, Yang Peng
2,3,4
, Alex Thomson
2,3,4
, Harpreet Arora
1,2
,
Robert Polski
1,2
, Yiran Zhang
1,2,3
, Hechen Ren
1,2
, Jason Alicea
2,3,4
, Gil Refael
2,3,4
, Felix von Oppen
2,5
,
Kenji Watanabe
6
, Takashi Taniguchi
6
and Stevan Nadj-Perge
1,2
*
1
T. J. Watson Laboratory of Applied Physics, California Institute of Technology, Pasadena, CA, USA.
2
Institute for Quantum Information and Matter,
California Institute of Technology, Pasadena, CA, USA.
3
Department of Physics, California Institute of Technology, Pasadena, CA, USA.
4
Walter Burke
Institute for Theoretical Physics, California Institute of Technology, Pasadena, CA, USA.
5
Dahlem Center for Complex Quantum Systems and Fachbereich
Physik, Freie Universität Berlin, Berlin, Germany.
6
National Institute for Materials Science, Tsukuba, Ibaraki, Japan. *e-mail: s.nadj-perge@caltech.edu
SUPPLEMENTARY INFORMATION
In the format provided by the authors and unedited.
NATuRE PHYSiCS
|
www.nature.com/naturephysics
Supplementary information for “
Electronic
Correlations
in Twisted Bilayer Graphene
Near
the Magic Angle”
Authors
:
Youngjoon Choi, Jeannette Kemmer, Yang Peng, Alex Thomson, Harpreet Arora,
Robert Polski, Yiran Zhang, Hechen Ren, Jason Alicea,
Gil Refael, Felix von Oppen, Kenji
Watanabe, Takashi Taniguchi and Stevan Nadj
-
Perge
Table of c
ontents
Section 1: Sample fabrication and measurement details
................................
................................
.............
1
Section 2: Data on 4.8nm moir
é
pattern
................................
................................
................................
......
1
Section 3: Tip
-
related effects
................................
................................
................................
........................
2
Section
4: Tip
-
gating correction at the CNP and the enhancement of the width
................................
........
7
Section 5: Areas with different twist angle
θ
................................
................................
................................
9
Section 6: Additional data
................................
................................
................................
...........................
10
Section 7: Theoretical modeling
................................
................................
................................
.................
12
Refere
nces:
................................
................................
................................
................................
.................
29
Supplementary Text
Section 1
:
Sample fabrication and measurement details
Figure S1 outlines the fabrication steps for twisted bilayer
graphene (TBG). We
us
ed
the
tear
-
and
-
twist
technique
following similar procedures outlined in Refs
1
-
3
. After picking up a 30
-
50nm
-
thick
boron nitride
flake using
Polydimethylsiloxane
(PDMS) coated with a P
oly(bisphenol
A carbonate)
(PC) polymer, the graph
ene is torn into two parts which were subsequently picked
up while controlling the twist between the parts. The entire stack is then transferred onto a
separate PDMS film in order to flip the order of the layers (Fig. S1 left panel). In this step, PC is
di
ssolved in
N
-
Methyl
-
2
-
pyrrolidone
(NMP). Afterward the PDMS with the inverted stack structure
is transferred onto a prepared chip
(silicon chip covered with 300nm of silicon
-
oxide)
with pre
-
patterned electrodes and a 10nm thick metallic graphene multilayer
that is used as the back gate
(Fig. S1 right panel). The twisted bilayer graphene (TBG) is then contacted to the gold electrodes
using additional few
-
layer graphene contacts. During device fabrication special care is taken so
that the temperature of the s
ample never exceeds 150
̊
C to avoid untwisting of the TBG
4
.
Fig
ure
S1
. Critical steps of
the
device fabrication
sequence. After picking up the BN and graphene
flakes, the stack is transferred and flipped onto a
second PDMS stamp that is used for a second
transfer onto the graphite back gate. In a separate
step, few
-
layer graphene is transferred to contact
TBG with a bias electrode.
Section
2
:
Data on
4.8nm
moir
é
pattern
Fig
ure
S2 shows a dI/dV map for the L
M
= 4.8nm (device D2) moir
é
pattern. The charge neutrality
point and both vHs peaks move in parallel
following the quadratic relation between charge
density
n
and chemical potential
μ
, |
n
|
=
α
μ
2
(up to an offset in V
Bias
and V
Bg
). This suggests that
in this regime the density of states between the vHs peaks follows approximately a linear
dispersion as expected for large twist angles. The vertical offset between the two branches is a
1
consequence of tip
-
induced gating, which is ex
plained in the following section. We note that this
effect may explain the apparent small increase in the distance between the van Hove singularities
at the CNP observed in Ref.
5
. Importantly, the relative enhancement of the peak splitting for
angles clos
e to the magic angle value is too large to be explained by this effect
as
dis
cu
ssed
in
S
ection 4
.
Fig
ure
S2. Spectroscopy
data taken on
the AA site for
θ
≈
3
°
, device D2. Dashed lines correspond
to a fit
|
n
|
=
α
μ
2
where
μ =
-
eV
Bias
and
n
= C
Bg
V
Bg
.
The small vertical offset of the two branches
that give rise to a small enhancement in splitting in the V
Bias
direction
that
is due to tip gating as
discussed in
the
following section.
Section
3
:
Tip
-
related effects
Dynamic strain effects
O
ccasionally
we
observe
topographies with pronounced hysteresis, especially for parameters
when the STM tip
-
sample distance is reduced (for example V
Bias
=
±
200mV, set point current
100pA or more)
.
Similar STM
-
tip
-
induced strain effects have been previously
reported
for
moir
é
patterns in single
-
layer graphene/boron nitride heterostructures
6
as well as in graphene on
silicon oxide
7
. The intuitive understanding of this effect is that when the interaction between the
tip and the sample is strong, the tip can slightly displace the carbon atoms as it scans over the
surface and consequently shift the moir
é
pattern boundari
es. We have performed a similar
analysis as in Ref.
6
and obtained strain maps illustrating this effect (see Fig. S3).
In addition to
2
th
ese
strain
-
related effects,
we have also observed other signatures of tip
-
graphene interaction,
such as hysteresis in th
e
current vs height
relation
I(z).
Figure S3.
(a)
T
he t
opography of the strained region close to the magic angle (L
M
= 13nm) in device
D1 and the corresponding strain
map (b).
The map represents the local lattice constant obtained
by coarsening
the topography grid into discrete pixels and performing a Fourier analysis of the
3nmx3nm area around each pixel and selecting one particular direction. Other directions show a
similar dependence. Scale bar is 10nm for both plots.
The dynamic strain effects are not observed consistently throughout both samples. The strain is
typically observed in areas close to the magic angle in device D1 and in certain areas away from
the magic
angle in device D2.
Also
,
in regions close to the mag
ic angle in device D2, the moir
é
superlattice is more uniform and there is overall less strain. From these observations we conclude
that the presence of intrinsic external strain is needed in order to observe the dynamic ti
p
-
induced strain effect.
D
ifferent STM microtips may have different sensitivities to these effects.
I
mportantly
,
we d
o
not see
a
qualitative difference in
the
spectroscopy data between the two
devices indicating that th
ese
effects
do not play a significant role
in determining the l
ocal density
of states
. This is consistent with
the
previous
ly
reported observation
s in Ref.
6
.
While this effect
produces artifacts in the topography, it does not change it permanently so the distance between
AA sites and hence the twist angle
θ
can still
be obtained accurately.
3
E
ffects of tip screening
Due to the close proximity between the STM tip and the sample, one might expect that the
presence of metallic objects will facilitate screening of electronic correlations. In this section, we
estimate
the s
trength of the electron
-
electron interaction in this system. Indeed, in the presence
of the metal tip which provides additional screening the interaction energy scale can be greatly
suppressed compared to the pristine unscreened case. A naïve estimate for
the electrostatic
energy for two electrons placed L
M
=
13nm apart is given by e
2
/(4
πεε
BN
L
M
) = 36meV. This estimate
is enhanced further to approximately 50meV if we consider that we have the boron nitride
dielectric only on one side. However, in the presenc
e of a metallic tip this interaction is locally
screened considerably. For a tip
-
sample distance of 1nm (taken as an upper limit), the Coulomb
-
energy scale would scale as (
1
/