Zgryźliwość kojarzy mi się z radością, która źle skończyła.
PHYSICAL REVIEW B
77
, 233406
2008
Plasmon spectroscopy of free-standing graphene films
T. Eberlein,
1
U. Bangert,
2
R. R. Nair,
2,4
R. Jones,
1
M. Gass,
3
A. L. Bleloch,
3
K. S. Novoselov,
4
A. Geim,
4
and P. R. Briddon
5
1
School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL, United Kingdom
2
School of Materials, The University of Manchester, Manchester M1 7HS, United Kingdom
3
SuperSTEM Laboratories, CCLRC Daresbury Laboratory, Warrington WA4 4AD, United Kingdom
4
School of Physics and Astronomy, The University of Manchester, Manchester M1 7HS, United Kingdom
5
School of Natural Sciences, University of Newcastle upon Tyne, Newcastle upon Tyne NE1 7RU, United Kingdom
Received 5 November 2007; revised manuscript received 25 January 2008; published 18 June 2008
Plasmon spectroscopy of the thinnest possible membrane, a single layer of carbon atoms: graphene, has been
carried out in conjunction with
ab initio
calculations of the low loss function. We observe
and
+
-surface
plasmon modes in free-standing single sheets at 4.7 and 14.6 eV, which are substantially redshifted from their
values in graphite. These modes are in very good agreement with the theoretical spectra, which find the
- and
+
in-plane modes of graphene at 4.8 and 14.5 eV. We also find that there is little loss caused by out-of-plane
modes for energies less than about 10 eV.
DOI:
PACS number
s
: 73.20.
r, 81.05.Uw, 71.45.Gm, 79.20.Uv
There has recently been intense interest in the properties
of graphene
1
,
2
and, in particular, those properties that distin-
guish it from graphite. Among the many suggested areas
2
in
which graphene may excel is, e.g., its application for sensors
due to the sensitivity of its electronic structure to
adsorbates.
3
Low loss energy electron spectroscopy provides
a way of detecting changes in the electronic structure, which
are highly spatially resolved. To be able to employ this tech-
nique, the loss spectra of graphene and graphite must be
clearly understood.
From an experimental viewpoint, an essential task in the
research into two-dimensional
butions of the latter split in energy when the product of
thickness
d
and momentum
q
,
dq
0; in graphite the maxi-
→
mum and/or minimum energy of the
→
transition tends
to roughly 7 and 5 eV and of the
→
transition to 20 and
15 eV, respectively.
5
The
E
field of a fast moving particle is elongated along its
direction of travel, therefore, when passing perpendicularly
through a graphene foil, mainly, the out-of plane mode with
momentum
q
parallel to
E
should be excited. However, as
shown below, these modes are forbidden in a single layer and
they have weak intensity in graphite. In an EELS experiment
carried out in a scanning transmission electron microscope
structures is to provide
evidence that they do indeed exist. This is especially impor-
tant because theory does not allow the existence of perfect
crystals in 2D space. The most conclusive evidence for the
existence of free-standing graphene has so far been obtained
from electron diffraction experiments
4
and all previously
published high resolution electron microscopy
2D
, although the momentum transfer is close to zero,
nonetheless,
q
has a considerable in-plane component be-
cause the collection angle is several millirad. For this reason,
we will observe those
STEM
surface and bulk plasmons excited
with
q
parallel to
a
.
Large graphene membranes were prepared by a microme-
chanical cleavage
1
of natural graphite on top of an oxidized
Si wafer. This deposition technique has the benefit of allow-
ing a quick and easy identification of mono- and multiple
layers by the additional optical contrast with respect to the
oxidized wafer. By using photolithography, a perforated
copper-gold film, serving as scaffold that could be used in
standard transmission electron microcopy holders, was de-
posited on top of the graphene crystallites. This scaffold was
lifted off the Si wafer, leaving the graphene attached to it.
15
Prior to electron microscopy, optical microscopy was used to
identify regions of one to few layers. HAADF images of
one-, two-, and five-layer regions are shown in Fig.
1
; the
rectangular boxes show the regions in which the intensity
im-
ages, to our knowledge, are of bilayers. By carrying out
highly spatially resolved electron energy loss spectroscopy
HREM
, we observe specific redshifts in the frequency of
plasmons in sample positions concomitant with single
graphene sheets, very similar to those observed in single-
wall carbon nanotubes
EELS
.
5
We also find further evi-
dence for existence of a single sheet by combining optical
and high angle annular dark field
SWCNTs
imaging.
Surface plasmon behavior in thin metal sheets is well
documented experimentally and explained using dielectric
theory.
6
,
7
There are also numerous reports on surface plas-
mons in graphite and carbon nanotubes.
5
,
8
–
11
Plasmon behav-
ior of truly 2D graphite, i.e., of monolayer graphite foils has
been theoretically suggested in a number of papers
12
,
13
and
has been experimentally studied for flat monolayers grown
on TiC.
14
However, there has been, so far, no experimental
study for free-standing sheets. There have, however, been
studies of SWCNTs and much of the interpretation of plas-
mon behavior for SWCNTs for radius
r
HAADF
, integrated over the width of the box, is traced.
The heights of the contrast profiles are multiples of the
smallest height measured in Fig.
1
Fig.
1
d
a
, confirming that the
contrast in the dark region
arises from one single
graphene sheet. Close inspection of the dark areas in each of
the images reveals uniform contrast upon which the hexago-
nal atomic lattice can be seen. Filtering of the HAADF im-
ages enhances the visibility of the lattice
Fig.
1
a
can be applied
to free-standing single graphene sheets,
5
particularly as the
tubes are free standing also. A characteristic of thin foils is
the vanishing of the bulk plasmon mode, leaving only the
surface plasmon mode; the out-of-plane and in-plane contri-
→
white-framed
boxes in Figs.
1
. The images are proof that the dark
areas are pristinely clean, i.e., free from contaminants.
a
–
1
c
1098-0121/2008/77
23
/233406
4
233406-1
©2008 The American Physical Society
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
FIG. 1.
Color online
a
–
c
HAADF
STEM images of ‘clean’ patches revealing
one, two, and five layers of graphene. The
area in the white-framed boxes has been sub-
jected to a low pass filter to disclose the atoms
more clearly: In the single layer the atoms in
the six rings are white and the hexagon center
is black.
d
Intensity traces have been taken
along the long dimensions of the rectangular
cyan-framed boxes in
a
–
c
; these traces
were then averaged over 40 pixels
width of
boxes
. They show step-wise increase of the
HAADF intensity corresponding to the layer
number. The vacuum intensity is at
268.
e
EEL spectra of one, two, five and several lay-
ers of graphene showing the
-and the
+
plasmon. The spectra are extracted from spec-
trum images; they are background subtracted
and each summed over 25 pixels.
(c)
(a)
(b)
(e)
(d)
The Daresbury aberration corrected SuperSTEM with an
Enfina EELS spectrometer was used for HREM bright field
and HAADF imaging in conjunction with EELS. EELS of
the plasmon region was carried out with a dispersion of 0.05
eV per channel at an energy resolution
approximately doubles. A plateau between 15 and 20 eV
starts to appear, which becomes pronounced for increasing
numbers of sheets: In the five-layer patch features above 15
eV start to appear, the plasmon maximum keeps moving to
higher energies, accompanied by further broadening, and the
integrated intensity increases to
full width at half
maximum of the zero loss peak
of 0.3–0.4 eV. The collec-
tion semiangle was approximately 19 mrad. The operating
voltage was 100 kV. In order to minimize beam damage, the
acquisition time per EEL spectrum was kept short
5 times the value of the
single sheet. The plasmon structure for more than 10 sheets
strongly resembles that of the graphite. The plasmon charac-
teristics are, thus, supreme indicators for the presence of
single layers. However, the smallest amount of contamina-
tion, e.g., caused by individual molecular adsorbates, intro-
duces a “three-dimensional” component and leads to a break-
down of the 2D behavior, as is frequently evidenced.
17
Contaminants can easily be discerned in HREM and
HAADF images; in fact the whitish contrast bordering the
clean patches in the HAADF images Figs.
1
0.1 s
;
the dose each pixel
in a spectrum image received was
10
7
electrons. When the acquisition time was in-
creased or when repeated scans were performed on the same
area, electron beam damage became manifest in formation of
point defects
1.6
.
16
single and multiple vacancies
, we show raw EEL data obtained by subtract-
ing the zero loss peak under identical conditions. The spectra
were all taken under the same acquisition conditions, hence,
direct comparison between them is possible.
The plasmon excitations in graphitic structures consist of
In Fig.
1
e
repre-
sent exactly this. We note that unlike low loss spectra, core
loss spectra were not found appropriate to reveal character-
istics that could help distinguish one layer from few layers
by their shape. In the spectrum of a larger number of
graphene sheets in Fig.
1
a
–
1
c
the
plasmons both exhibiting bulk and surface
modes. In Fig.
1
- and
+
, the plasmon spectra taken in clean
patches of one, two, five and multiple sheets shows that the
e
e
, resembling much that of a thin
graphitic slab,
the well-known graphitic bulk
- and
-mode, at 7 eV in graphite, has shifted to 4.7 eV in the
single layer. Furthermore, the spectrum here exhibits only
the
-modes can be observed at 7 and 26 eV, respectively. We
have consistently observed
-plasmon behavior with distinct
occurrence of higher energy components, accompanied by a
step-wise increase in the HAADF intensity when going from
one to two or more sheets. Values of the HAADF intensity
for one, two, and more sheets
-surface mode at 14.6 eV and the 26 eV plasmon in
graphite is not present. The shape and intensity of the one-
layer plasmon structure was repeatedly measured in different
places on the same sample, in different samples, and even in
different experimental sessions. Given the same acquisition
conditions
+
were measured
against the vacuum HAADF intensity and displayed incre-
ments of approximately integer multiples of the first incre-
ment above the vacuum intensity, i.e., the intensity of one
sheet. Data acquired from a number of experiments, carried
out on different samples and different days, yielded the same
results for the HAADF and the integrated EEL intensity.
The graphene plasmon behavior can be compared with
see Fig.
1
d
e.g., energy dispersion, electron beam current,
and dwell
, nearly identical spectra were obtained,
which, at the same time, constituted the lowest plasmon sig-
nal measured overall. For two sheets, the triangular shape of
the
time
resonance shifts toward higher energies and the
integrated intensity under the peak
+
between 10 and 40 eV
233406-2
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
that of thin carbon nanotubes, where the bulk plasmon sub-
sides and only tangential and radial surface modes prevail. A
trend has been observed for the
3
q||c 20(c/a)
q||c 10(c/a)
q||c 1(c/a)
2.5
plasmon in tubes, where
1,
5
namely, that the tangential
r
/
R
−
or in plane
mode
2
increases on the expense of the radial
or out of plane
and,
1.5
looking at Fig. 3
in Ref.
5
, where plasmon spectra for
tubes with increasing
r
are presented, it furthermore appears
that for large
r
, only the tangential mode
a
1
re-
mains, very similar to the graphene case. However, although
the resemblance with single-wall nanotube EEL data is large,
it cannot, from the start, be assumed that the graphene low
loss function reduces to the same theoretical approximation,
namely, to that derived for an electron passing in aloof
geometry.
18
For large tube diameters, in the work by Taverna
et al.
,
18
the dielectric response tends to that of an anisotropic
thin foil in the weak coupling regime, however, in their case,
the electron passes parallel at a distance, whereas for a pen-
etrating geometry, as in our case, the potential would have to
be modified for electron transit within the foil. There has
been no expression of the low loss function derived using
dielectric theory and continuum models for the latter case.
Here, we derive the microscopic loss function for this sce-
nario via
ab initio
calculations.
These deal with supercells composed of carbon layers. We
imagine that the EELS experiment can be modeled by a
monochromatic beam of electrons whose momentum transfer
q
to the sample is along either the
c
or
a
axis of the unit cell.
The unit cell contains planes of carbon atoms, which are
separated by multiples of the interlayer separation found in
graphite. Thus, to model graphene, layers of carbon atoms
are removed from the supercell leaving single layers well
separated from each other. Bilayers and trilayers can be mod-
eled by removing layers as before but leaving pairs of planes
with the separation and stacking sequence found in graphite
but each bilayer or trilayer is well separated from its neigh-
bor.
The rate at which energy and momentum
h
at
15 eV
0.5
0
0
5
10
15
20
25
30
35
40
3.5
q||a 20(c/a)
q||a 10(c/a)
q||a 1(c/a)
3
2.5
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
The loss function for graphite with a
separation between
c
-planes equal to 20 and 10 times its separation
in graphite and for graphite itself. Figure
2
top
is case where
q
c
and Fig.
2
bottom
when
q
a
.
Y
axis: arbitrary units,
x
axis: energy
loss in eV.
FIG. 2.
Color online
We first consider graphite. The top and bottom panels of
Fig.
2
show the loss function for
q
parallel to
c
and
a
blue
dotted-line spectra
, respectively. We find that the former has
peaks around 4, 12, 14, and 18 eV while the latter has a peak
at 7 eV and a broad peak at 27 eV. The 7 and 27 eV peaks are
due to
-electron plasmons
one electron per atom
and
plus
plasmons reflecting the van-
ishing of the real part of the dielectric constant.
The
four electrons per atom
q
is lost from a
charged particle moving through a homogeneous dielectric
with speed
v
due
to ionization is proportional to the loss
function, which is the imaginary part of the inverse longitu-
dinal dielectric constant
4 eV peak is due to a
to
transition. It is
noteworthy that the height of the
4 eV peak is about 10%
of the height of the 7 eV peak and, hence, the main plasmon
loss below 10 eV is due to the plasmon moving along
a
.
These results are in reasonable agreement with studies of
the
c
and
a
plasmons in graphite:
23
Here, for
q
parallel to
c
,
peaks are observed at 4.6, 13, 15, and 19 eV but for
q
par-
allel to
a
, only peaks at 6.8 eV and a broad peak at 23 eV are
seen. Moreover, the intensity of the 4.6 eV peak is about
20% that of the 6.8 eV peak. This experiment shows that
energy loss in a transmission experiment, such as used here,
will be due to both
q
parallel to
a
and
q
parallel to
c
plas-
mons. The close agreement with the calculations suggests
that the loss function Im
ij
.
19
Within the random phase
approximation,
20
,
21
the expression for
ij
is
k
e
2
0
r
·
i
k
+
qc
,
kv
E
k
+
q
,
c
−
E
k
,
v
−
E
kv
,
k
+
qc
r
·
j
ij
q
,
E
=
ij
+
.
Here,
E
k,c
and
E
k,
v
denote the empty and filled bands and
Ω
is the volume of the unit cell. The sum is over a special set of
k
vectors. For graphene and graphite, the principal values of
the loss tensor lie along
c
and
a
. This formulation includes
the effect of multiple inelastic scattering but assumes a ho-
mogeneous system. We use the AIMPRO local density func-
tional code
22
to evaluate the imaginary part of the dielectric
tensor for
q
=0 and use the Kramers–Kronig relations to de-
rive its real part. We use 43 000
k
-points and a broadening of
0.5 eV to fully converge the results. The loss function is then
found for graphite, graphene, and bimultilayers and trimulti-
layers.
is an appropriate way of mod-
eling the loss for planar sheets of graphite.
We now turn to graphene. This was modeled by expand-
ing the lattice parameter along
c
. The loss functions for
q
parallel to
c
and
a
for 10 and 20 fold expansions are shown
in Fig.
2
as green dotted-line and red full-line spectra, re-
spectively. In the latter case where
q
is parallel to
c
1
/
top
panel
and after
this the onset occurs at a similar energy to graphite. How-
ever, the peak heights are very different from graphite. The
, the loss is almost zero up to 12 eV
red line
233406-3
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
disappearance of the 4 eV peak can be explained by a selec-
tion rule. As stated above, the peak is due to a transition
between occupied
3
q||c single graphene
q||c bilayer graphene
q||c trilayer graphene
q||c graphite
2.5
bands at the
M
point
of the Brillouin zone. Inspection of the wave functions for
these two states shows that both transform as
p
z
and are odd
under reflection symmetry present in the basal plane of
graphene but not AB graphite and, hence, the dipole matrix
element between them vanishes for transitions for which
q
is
parallel to
c
although the transition is allowed for
q
parallel
to
a
. This shows that the dielectric constant of graphene is
not the same as graphite.
Figure
2
, bottom panel, compares the loss function for
q
parallel to
a
for graphene and graphite. We note there are
substantial redshifts of the peaks found in graphite.
13
The 7
eV
and unoccupied
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
3
q||a single graphene
q||a bilayer graphene
q||a trilayer graphene
q||a graphite
2.5
plasmon peak has shifted downward to about 4.8 eV
while the broad peak around 27 eV has sharpened and shifted
to below 15 eV. Such shifts are seen in the experimental
spectra shown in Fig.
1
.
We now investigate the loss functions for bilayers and
trilayers stacked as in graphite. The precise peak positions
depend on the separation of planes and to compare the loss
functions for graphite and multiple layers, we choose the
separation between periodically repeated multilayers to be 5
times the separation in graphite. We also use a broadening of
1.5 eV. The supercell containing the trilayer for example, has
three layers of graphene separated by the interlayer separa-
tion found in graphite, but the separation of these planes
from similar planes in adjacent unit cells along the
c
axis is
now five times the separation in graphite. The resulting loss
function is shown in Fig.
3
for single, double, and triple
layers, as well as for graphite. It is clear that there is an
increasing redshift of the peaks above, about 10 eV as the
number of layers decreases. The relative increase in ampli-
tudes of peaks for the different layers
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
FIG. 3.
Color online
Comparison of the loss function for a
single and multilayers of graphene for
q
c
top
and
q
a
bottom
.
Note, the loss around 4 eV in the
q
c
case for the multilayer case
and its absence for single layer graphene
shown as red line
. Note
also the increasing redshift of the main peaks above about 10 eV as
the number of layers decreases.
Y
axis: arbitrary units,
x
axis: en-
ergy loss in eV.
for a bilayer with AA stacking. It implies that any observa-
tions of a loss below 10 eV due to these plasmons must be
due to adsorbates lying on graphene and makes graphene
peculiarly sensitive to such adsorbates. We note that these
plasmons could be excited by light of grazing incidence and
polarized along
c
. A further feature, which is unique to the
plasmon behavior of graphene, is the shift of the 7 eV in-
plane plasmon seen in graphite to about 4.7 eV, as well as a
substantial redshift of the broad plasmon peak at 25 eV, seen
in graphite to about 14.6 eV in graphene accompanied by a
distinct shape change to skewed triangular. The redshifts de-
crease as the number of close-by layers increases.
seem roughly consistent with experimental spectra in Fig.
1
for one, two, and five layers but the observed spectra are
broader. All the layers except the single one show the out-
of-plane plasmon peak around 4 eV.
In conclusion, the 4.6 eV out-of-plane loss peak found in
graphite disappears for a single layer graphene. This is re-
lated to a selection rule operating in graphene but not in
graphite. This, however, does not occur for in-plane modes
or a graphene bilayer with AB stacking but does also occur
Fig.
3
, bottom panel
1
K. S. Novoselov
et al.
, Science
306
, 666
.
2
A. K. Geim and K. S. Novoselov, Nat. Mater.
6
, 183
2004
.
13
A. G. Marinopoulos
et al.
, Appl. Phys. A: Mater. Sci. Process.
78
, 1157
2007
2007
.
3
F. Schedin
et al.
, Nat. Mater.
6
, 652
2007
.
.
14
A. Nagashima
et al.
, Solid State Commun.
83
, 581
2004
4
C. J. Meyer
et al.
, Nature
London
446
,60
2007
.
1992
.
5
O. Stéphan
et al.
, Phys. Rev. B
66
, 155422
2002
.
15
R. R. Nair
et al.
, Science
to be published
.
6
H. Raether, Phys. Thin Films
9
, 149
1977
.
7
A. Otto, Phys. Status Solidi
22
, 401
1967
.
8
J. F. Annett
et al.
, Phys. Rev. B
37
, 2408
1988
.
9
P. Laitenberger and R. E. Palmer, Phys. Rev. Lett.
76
, 1952
1996
.
10
R. Perez and W. Que, J. Phys.: Condens. Matter
18
, 3197
2006
.
11
M. Kociak
et al.
, Phys. Rev. Lett.
87
, 075501
2001
.
12
E. H. Hwang and S. Das Sarma, Phys. Rev. B
75
, 205418
16
U. Bangert
et al.
unpublished
.
17
M. H. Gass
et al.
unpublished
.
18
D. Taverna
et al.
, Phys. Rev. B
66
, 235419
2002
.
19
P. Nozieres and D. Pines, Phys. Rev.
113
, 1254
1959
.
20
H. Ehrenreich and M. H. Cohen, Phys. Rev.
115
, 786
1959
.
21
S. L. Adler, Phys. Rev.
126
, 413
1962
.
22
U. Bangert
et al.
, Philos. Mag.
86
, 4757
2006
.
23
K. Zeppenfeld, Phys. Lett.
25A
, 335
1967
.
233406-4
zanotowane.pl doc.pisz.pl pdf.pisz.pl hannaeva.xlx.pl
77
, 233406
2008
Plasmon spectroscopy of free-standing graphene films
T. Eberlein,
1
U. Bangert,
2
R. R. Nair,
2,4
R. Jones,
1
M. Gass,
3
A. L. Bleloch,
3
K. S. Novoselov,
4
A. Geim,
4
and P. R. Briddon
5
1
School of Physics, University of Exeter, Stocker Road, Exeter, EX4 4QL, United Kingdom
2
School of Materials, The University of Manchester, Manchester M1 7HS, United Kingdom
3
SuperSTEM Laboratories, CCLRC Daresbury Laboratory, Warrington WA4 4AD, United Kingdom
4
School of Physics and Astronomy, The University of Manchester, Manchester M1 7HS, United Kingdom
5
School of Natural Sciences, University of Newcastle upon Tyne, Newcastle upon Tyne NE1 7RU, United Kingdom
Received 5 November 2007; revised manuscript received 25 January 2008; published 18 June 2008
Plasmon spectroscopy of the thinnest possible membrane, a single layer of carbon atoms: graphene, has been
carried out in conjunction with
ab initio
calculations of the low loss function. We observe
and
+
-surface
plasmon modes in free-standing single sheets at 4.7 and 14.6 eV, which are substantially redshifted from their
values in graphite. These modes are in very good agreement with the theoretical spectra, which find the
- and
+
in-plane modes of graphene at 4.8 and 14.5 eV. We also find that there is little loss caused by out-of-plane
modes for energies less than about 10 eV.
DOI:
PACS number
s
: 73.20.
r, 81.05.Uw, 71.45.Gm, 79.20.Uv
There has recently been intense interest in the properties
of graphene
1
,
2
and, in particular, those properties that distin-
guish it from graphite. Among the many suggested areas
2
in
which graphene may excel is, e.g., its application for sensors
due to the sensitivity of its electronic structure to
adsorbates.
3
Low loss energy electron spectroscopy provides
a way of detecting changes in the electronic structure, which
are highly spatially resolved. To be able to employ this tech-
nique, the loss spectra of graphene and graphite must be
clearly understood.
From an experimental viewpoint, an essential task in the
research into two-dimensional
butions of the latter split in energy when the product of
thickness
d
and momentum
q
,
dq
0; in graphite the maxi-
→
mum and/or minimum energy of the
→
transition tends
to roughly 7 and 5 eV and of the
→
transition to 20 and
15 eV, respectively.
5
The
E
field of a fast moving particle is elongated along its
direction of travel, therefore, when passing perpendicularly
through a graphene foil, mainly, the out-of plane mode with
momentum
q
parallel to
E
should be excited. However, as
shown below, these modes are forbidden in a single layer and
they have weak intensity in graphite. In an EELS experiment
carried out in a scanning transmission electron microscope
structures is to provide
evidence that they do indeed exist. This is especially impor-
tant because theory does not allow the existence of perfect
crystals in 2D space. The most conclusive evidence for the
existence of free-standing graphene has so far been obtained
from electron diffraction experiments
4
and all previously
published high resolution electron microscopy
2D
, although the momentum transfer is close to zero,
nonetheless,
q
has a considerable in-plane component be-
cause the collection angle is several millirad. For this reason,
we will observe those
STEM
surface and bulk plasmons excited
with
q
parallel to
a
.
Large graphene membranes were prepared by a microme-
chanical cleavage
1
of natural graphite on top of an oxidized
Si wafer. This deposition technique has the benefit of allow-
ing a quick and easy identification of mono- and multiple
layers by the additional optical contrast with respect to the
oxidized wafer. By using photolithography, a perforated
copper-gold film, serving as scaffold that could be used in
standard transmission electron microcopy holders, was de-
posited on top of the graphene crystallites. This scaffold was
lifted off the Si wafer, leaving the graphene attached to it.
15
Prior to electron microscopy, optical microscopy was used to
identify regions of one to few layers. HAADF images of
one-, two-, and five-layer regions are shown in Fig.
1
; the
rectangular boxes show the regions in which the intensity
im-
ages, to our knowledge, are of bilayers. By carrying out
highly spatially resolved electron energy loss spectroscopy
HREM
, we observe specific redshifts in the frequency of
plasmons in sample positions concomitant with single
graphene sheets, very similar to those observed in single-
wall carbon nanotubes
EELS
.
5
We also find further evi-
dence for existence of a single sheet by combining optical
and high angle annular dark field
SWCNTs
imaging.
Surface plasmon behavior in thin metal sheets is well
documented experimentally and explained using dielectric
theory.
6
,
7
There are also numerous reports on surface plas-
mons in graphite and carbon nanotubes.
5
,
8
–
11
Plasmon behav-
ior of truly 2D graphite, i.e., of monolayer graphite foils has
been theoretically suggested in a number of papers
12
,
13
and
has been experimentally studied for flat monolayers grown
on TiC.
14
However, there has been, so far, no experimental
study for free-standing sheets. There have, however, been
studies of SWCNTs and much of the interpretation of plas-
mon behavior for SWCNTs for radius
r
HAADF
, integrated over the width of the box, is traced.
The heights of the contrast profiles are multiples of the
smallest height measured in Fig.
1
Fig.
1
d
a
, confirming that the
contrast in the dark region
arises from one single
graphene sheet. Close inspection of the dark areas in each of
the images reveals uniform contrast upon which the hexago-
nal atomic lattice can be seen. Filtering of the HAADF im-
ages enhances the visibility of the lattice
Fig.
1
a
can be applied
to free-standing single graphene sheets,
5
particularly as the
tubes are free standing also. A characteristic of thin foils is
the vanishing of the bulk plasmon mode, leaving only the
surface plasmon mode; the out-of-plane and in-plane contri-
→
white-framed
boxes in Figs.
1
. The images are proof that the dark
areas are pristinely clean, i.e., free from contaminants.
a
–
1
c
1098-0121/2008/77
23
/233406
4
233406-1
©2008 The American Physical Society
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
FIG. 1.
Color online
a
–
c
HAADF
STEM images of ‘clean’ patches revealing
one, two, and five layers of graphene. The
area in the white-framed boxes has been sub-
jected to a low pass filter to disclose the atoms
more clearly: In the single layer the atoms in
the six rings are white and the hexagon center
is black.
d
Intensity traces have been taken
along the long dimensions of the rectangular
cyan-framed boxes in
a
–
c
; these traces
were then averaged over 40 pixels
width of
boxes
. They show step-wise increase of the
HAADF intensity corresponding to the layer
number. The vacuum intensity is at
268.
e
EEL spectra of one, two, five and several lay-
ers of graphene showing the
-and the
+
plasmon. The spectra are extracted from spec-
trum images; they are background subtracted
and each summed over 25 pixels.
(c)
(a)
(b)
(e)
(d)
The Daresbury aberration corrected SuperSTEM with an
Enfina EELS spectrometer was used for HREM bright field
and HAADF imaging in conjunction with EELS. EELS of
the plasmon region was carried out with a dispersion of 0.05
eV per channel at an energy resolution
approximately doubles. A plateau between 15 and 20 eV
starts to appear, which becomes pronounced for increasing
numbers of sheets: In the five-layer patch features above 15
eV start to appear, the plasmon maximum keeps moving to
higher energies, accompanied by further broadening, and the
integrated intensity increases to
full width at half
maximum of the zero loss peak
of 0.3–0.4 eV. The collec-
tion semiangle was approximately 19 mrad. The operating
voltage was 100 kV. In order to minimize beam damage, the
acquisition time per EEL spectrum was kept short
5 times the value of the
single sheet. The plasmon structure for more than 10 sheets
strongly resembles that of the graphite. The plasmon charac-
teristics are, thus, supreme indicators for the presence of
single layers. However, the smallest amount of contamina-
tion, e.g., caused by individual molecular adsorbates, intro-
duces a “three-dimensional” component and leads to a break-
down of the 2D behavior, as is frequently evidenced.
17
Contaminants can easily be discerned in HREM and
HAADF images; in fact the whitish contrast bordering the
clean patches in the HAADF images Figs.
1
0.1 s
;
the dose each pixel
in a spectrum image received was
10
7
electrons. When the acquisition time was in-
creased or when repeated scans were performed on the same
area, electron beam damage became manifest in formation of
point defects
1.6
.
16
single and multiple vacancies
, we show raw EEL data obtained by subtract-
ing the zero loss peak under identical conditions. The spectra
were all taken under the same acquisition conditions, hence,
direct comparison between them is possible.
The plasmon excitations in graphitic structures consist of
In Fig.
1
e
repre-
sent exactly this. We note that unlike low loss spectra, core
loss spectra were not found appropriate to reveal character-
istics that could help distinguish one layer from few layers
by their shape. In the spectrum of a larger number of
graphene sheets in Fig.
1
a
–
1
c
the
plasmons both exhibiting bulk and surface
modes. In Fig.
1
- and
+
, the plasmon spectra taken in clean
patches of one, two, five and multiple sheets shows that the
e
e
, resembling much that of a thin
graphitic slab,
the well-known graphitic bulk
- and
-mode, at 7 eV in graphite, has shifted to 4.7 eV in the
single layer. Furthermore, the spectrum here exhibits only
the
-modes can be observed at 7 and 26 eV, respectively. We
have consistently observed
-plasmon behavior with distinct
occurrence of higher energy components, accompanied by a
step-wise increase in the HAADF intensity when going from
one to two or more sheets. Values of the HAADF intensity
for one, two, and more sheets
-surface mode at 14.6 eV and the 26 eV plasmon in
graphite is not present. The shape and intensity of the one-
layer plasmon structure was repeatedly measured in different
places on the same sample, in different samples, and even in
different experimental sessions. Given the same acquisition
conditions
+
were measured
against the vacuum HAADF intensity and displayed incre-
ments of approximately integer multiples of the first incre-
ment above the vacuum intensity, i.e., the intensity of one
sheet. Data acquired from a number of experiments, carried
out on different samples and different days, yielded the same
results for the HAADF and the integrated EEL intensity.
The graphene plasmon behavior can be compared with
see Fig.
1
d
e.g., energy dispersion, electron beam current,
and dwell
, nearly identical spectra were obtained,
which, at the same time, constituted the lowest plasmon sig-
nal measured overall. For two sheets, the triangular shape of
the
time
resonance shifts toward higher energies and the
integrated intensity under the peak
+
between 10 and 40 eV
233406-2
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
that of thin carbon nanotubes, where the bulk plasmon sub-
sides and only tangential and radial surface modes prevail. A
trend has been observed for the
3
q||c 20(c/a)
q||c 10(c/a)
q||c 1(c/a)
2.5
plasmon in tubes, where
1,
5
namely, that the tangential
r
/
R
−
or in plane
mode
2
increases on the expense of the radial
or out of plane
and,
1.5
looking at Fig. 3
in Ref.
5
, where plasmon spectra for
tubes with increasing
r
are presented, it furthermore appears
that for large
r
, only the tangential mode
a
1
re-
mains, very similar to the graphene case. However, although
the resemblance with single-wall nanotube EEL data is large,
it cannot, from the start, be assumed that the graphene low
loss function reduces to the same theoretical approximation,
namely, to that derived for an electron passing in aloof
geometry.
18
For large tube diameters, in the work by Taverna
et al.
,
18
the dielectric response tends to that of an anisotropic
thin foil in the weak coupling regime, however, in their case,
the electron passes parallel at a distance, whereas for a pen-
etrating geometry, as in our case, the potential would have to
be modified for electron transit within the foil. There has
been no expression of the low loss function derived using
dielectric theory and continuum models for the latter case.
Here, we derive the microscopic loss function for this sce-
nario via
ab initio
calculations.
These deal with supercells composed of carbon layers. We
imagine that the EELS experiment can be modeled by a
monochromatic beam of electrons whose momentum transfer
q
to the sample is along either the
c
or
a
axis of the unit cell.
The unit cell contains planes of carbon atoms, which are
separated by multiples of the interlayer separation found in
graphite. Thus, to model graphene, layers of carbon atoms
are removed from the supercell leaving single layers well
separated from each other. Bilayers and trilayers can be mod-
eled by removing layers as before but leaving pairs of planes
with the separation and stacking sequence found in graphite
but each bilayer or trilayer is well separated from its neigh-
bor.
The rate at which energy and momentum
h
at
15 eV
0.5
0
0
5
10
15
20
25
30
35
40
3.5
q||a 20(c/a)
q||a 10(c/a)
q||a 1(c/a)
3
2.5
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
The loss function for graphite with a
separation between
c
-planes equal to 20 and 10 times its separation
in graphite and for graphite itself. Figure
2
top
is case where
q
c
and Fig.
2
bottom
when
q
a
.
Y
axis: arbitrary units,
x
axis: energy
loss in eV.
FIG. 2.
Color online
We first consider graphite. The top and bottom panels of
Fig.
2
show the loss function for
q
parallel to
c
and
a
blue
dotted-line spectra
, respectively. We find that the former has
peaks around 4, 12, 14, and 18 eV while the latter has a peak
at 7 eV and a broad peak at 27 eV. The 7 and 27 eV peaks are
due to
-electron plasmons
one electron per atom
and
plus
plasmons reflecting the van-
ishing of the real part of the dielectric constant.
The
four electrons per atom
q
is lost from a
charged particle moving through a homogeneous dielectric
with speed
v
due
to ionization is proportional to the loss
function, which is the imaginary part of the inverse longitu-
dinal dielectric constant
4 eV peak is due to a
to
transition. It is
noteworthy that the height of the
4 eV peak is about 10%
of the height of the 7 eV peak and, hence, the main plasmon
loss below 10 eV is due to the plasmon moving along
a
.
These results are in reasonable agreement with studies of
the
c
and
a
plasmons in graphite:
23
Here, for
q
parallel to
c
,
peaks are observed at 4.6, 13, 15, and 19 eV but for
q
par-
allel to
a
, only peaks at 6.8 eV and a broad peak at 23 eV are
seen. Moreover, the intensity of the 4.6 eV peak is about
20% that of the 6.8 eV peak. This experiment shows that
energy loss in a transmission experiment, such as used here,
will be due to both
q
parallel to
a
and
q
parallel to
c
plas-
mons. The close agreement with the calculations suggests
that the loss function Im
ij
.
19
Within the random phase
approximation,
20
,
21
the expression for
ij
is
k
e
2
0
r
·
i
k
+
qc
,
kv
E
k
+
q
,
c
−
E
k
,
v
−
E
kv
,
k
+
qc
r
·
j
ij
q
,
E
=
ij
+
.
Here,
E
k,c
and
E
k,
v
denote the empty and filled bands and
Ω
is the volume of the unit cell. The sum is over a special set of
k
vectors. For graphene and graphite, the principal values of
the loss tensor lie along
c
and
a
. This formulation includes
the effect of multiple inelastic scattering but assumes a ho-
mogeneous system. We use the AIMPRO local density func-
tional code
22
to evaluate the imaginary part of the dielectric
tensor for
q
=0 and use the Kramers–Kronig relations to de-
rive its real part. We use 43 000
k
-points and a broadening of
0.5 eV to fully converge the results. The loss function is then
found for graphite, graphene, and bimultilayers and trimulti-
layers.
is an appropriate way of mod-
eling the loss for planar sheets of graphite.
We now turn to graphene. This was modeled by expand-
ing the lattice parameter along
c
. The loss functions for
q
parallel to
c
and
a
for 10 and 20 fold expansions are shown
in Fig.
2
as green dotted-line and red full-line spectra, re-
spectively. In the latter case where
q
is parallel to
c
1
/
top
panel
and after
this the onset occurs at a similar energy to graphite. How-
ever, the peak heights are very different from graphite. The
, the loss is almost zero up to 12 eV
red line
233406-3
BRIEF REPORTS
PHYSICAL REVIEW B
77
, 233406
2008
disappearance of the 4 eV peak can be explained by a selec-
tion rule. As stated above, the peak is due to a transition
between occupied
3
q||c single graphene
q||c bilayer graphene
q||c trilayer graphene
q||c graphite
2.5
bands at the
M
point
of the Brillouin zone. Inspection of the wave functions for
these two states shows that both transform as
p
z
and are odd
under reflection symmetry present in the basal plane of
graphene but not AB graphite and, hence, the dipole matrix
element between them vanishes for transitions for which
q
is
parallel to
c
although the transition is allowed for
q
parallel
to
a
. This shows that the dielectric constant of graphene is
not the same as graphite.
Figure
2
, bottom panel, compares the loss function for
q
parallel to
a
for graphene and graphite. We note there are
substantial redshifts of the peaks found in graphite.
13
The 7
eV
and unoccupied
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
3
q||a single graphene
q||a bilayer graphene
q||a trilayer graphene
q||a graphite
2.5
plasmon peak has shifted downward to about 4.8 eV
while the broad peak around 27 eV has sharpened and shifted
to below 15 eV. Such shifts are seen in the experimental
spectra shown in Fig.
1
.
We now investigate the loss functions for bilayers and
trilayers stacked as in graphite. The precise peak positions
depend on the separation of planes and to compare the loss
functions for graphite and multiple layers, we choose the
separation between periodically repeated multilayers to be 5
times the separation in graphite. We also use a broadening of
1.5 eV. The supercell containing the trilayer for example, has
three layers of graphene separated by the interlayer separa-
tion found in graphite, but the separation of these planes
from similar planes in adjacent unit cells along the
c
axis is
now five times the separation in graphite. The resulting loss
function is shown in Fig.
3
for single, double, and triple
layers, as well as for graphite. It is clear that there is an
increasing redshift of the peaks above, about 10 eV as the
number of layers decreases. The relative increase in ampli-
tudes of peaks for the different layers
2
1.5
1
0.5
0
0
5
10
15
20
25
30
35
40
FIG. 3.
Color online
Comparison of the loss function for a
single and multilayers of graphene for
q
c
top
and
q
a
bottom
.
Note, the loss around 4 eV in the
q
c
case for the multilayer case
and its absence for single layer graphene
shown as red line
. Note
also the increasing redshift of the main peaks above about 10 eV as
the number of layers decreases.
Y
axis: arbitrary units,
x
axis: en-
ergy loss in eV.
for a bilayer with AA stacking. It implies that any observa-
tions of a loss below 10 eV due to these plasmons must be
due to adsorbates lying on graphene and makes graphene
peculiarly sensitive to such adsorbates. We note that these
plasmons could be excited by light of grazing incidence and
polarized along
c
. A further feature, which is unique to the
plasmon behavior of graphene, is the shift of the 7 eV in-
plane plasmon seen in graphite to about 4.7 eV, as well as a
substantial redshift of the broad plasmon peak at 25 eV, seen
in graphite to about 14.6 eV in graphene accompanied by a
distinct shape change to skewed triangular. The redshifts de-
crease as the number of close-by layers increases.
seem roughly consistent with experimental spectra in Fig.
1
for one, two, and five layers but the observed spectra are
broader. All the layers except the single one show the out-
of-plane plasmon peak around 4 eV.
In conclusion, the 4.6 eV out-of-plane loss peak found in
graphite disappears for a single layer graphene. This is re-
lated to a selection rule operating in graphene but not in
graphite. This, however, does not occur for in-plane modes
or a graphene bilayer with AB stacking but does also occur
Fig.
3
, bottom panel
1
K. S. Novoselov
et al.
, Science
306
, 666
.
2
A. K. Geim and K. S. Novoselov, Nat. Mater.
6
, 183
2004
.
13
A. G. Marinopoulos
et al.
, Appl. Phys. A: Mater. Sci. Process.
78
, 1157
2007
2007
.
3
F. Schedin
et al.
, Nat. Mater.
6
, 652
2007
.
.
14
A. Nagashima
et al.
, Solid State Commun.
83
, 581
2004
4
C. J. Meyer
et al.
, Nature
London
446
,60
2007
.
1992
.
5
O. Stéphan
et al.
, Phys. Rev. B
66
, 155422
2002
.
15
R. R. Nair
et al.
, Science
to be published
.
6
H. Raether, Phys. Thin Films
9
, 149
1977
.
7
A. Otto, Phys. Status Solidi
22
, 401
1967
.
8
J. F. Annett
et al.
, Phys. Rev. B
37
, 2408
1988
.
9
P. Laitenberger and R. E. Palmer, Phys. Rev. Lett.
76
, 1952
1996
.
10
R. Perez and W. Que, J. Phys.: Condens. Matter
18
, 3197
2006
.
11
M. Kociak
et al.
, Phys. Rev. Lett.
87
, 075501
2001
.
12
E. H. Hwang and S. Das Sarma, Phys. Rev. B
75
, 205418
16
U. Bangert
et al.
unpublished
.
17
M. H. Gass
et al.
unpublished
.
18
D. Taverna
et al.
, Phys. Rev. B
66
, 235419
2002
.
19
P. Nozieres and D. Pines, Phys. Rev.
113
, 1254
1959
.
20
H. Ehrenreich and M. H. Cohen, Phys. Rev.
115
, 786
1959
.
21
S. L. Adler, Phys. Rev.
126
, 413
1962
.
22
U. Bangert
et al.
, Philos. Mag.
86
, 4757
2006
.
23
K. Zeppenfeld, Phys. Lett.
25A
, 335
1967
.
233406-4