Stone-Wales–type transformations in carbon nanostructures
driven by electron irradiation
arXiv:1105.1617v1 [cond-mat.mtrl-sci] 9 May 2011
J. Kotakoski1,∗ , J. C. Meyer2,† , S. Kurasch2 , D.
Santos-Cottin1 , U. Kaiser2 and A. V. Krasheninnikov1,3
1
Department of Physics,
University of Helsinki, P.O. Box 43,
00014 Helsinki, Finland
2
Electron microscopy of materials science,
University of Ulm, Germany
3
Department of Applied Physics,
Aalto University, P.O. Box 1100,
00076 Aalto, Finland
(Dated: May 10, 2011)
Abstract
Observations of topological defects associated with Stone-Wales–type transformations (i.e., bond
rotations) in high resolution transmission electron microscopy (HRTEM) images of carbon nanostructures are at odds with the equilibrium thermodynamics of these systems. Here, by combining
aberration-corrected HRTEM experiments and atomistic simulations, we show that such defects
can be formed by single electron impacts, and remarkably, at electron energies below the threshold
for atomic displacements. We further study the mechanisms of irradiation-driven bond rotations,
and explain why electron irradiation at moderate electron energies (∼100 keV) tends to amorphize
rather than perforate graphene. We also show via simulations that Stone-Wales defects can appear in curved graphitic structures due to incomplete recombination of irradiation-induced Frenkel
defects, similar to formation of Wigner-type defects in silicon.
PACS numbers: 68.37.Og, 81.05.ue, 64.70.Nd, 31.15.es
∗
†
Corresponding author, email: jani.kotakoski@iki.fi
Present address: University of Vienna, Department of Physics, 1090 Wien, Austria
1
INTRODUCTION
Stone-Wales defect [1, 2] – SW (55-77) – is the simplest example of topological disorder
in graphene and other sp2 -hybridized carbon systems. It can be formed by rotating a C-C
bond by 90◦ with regard to the midpoint of the bond – referred to as the SW transformation
– so that four hexagons are turned into two pentagons and two heptagons. This defect
has received considerable amount of attention [3–8], because it has the lowest formation
energy among all intrinsic defects in graphenic systems, and because it presumably plays
an important role in plastic deformation of carbon nanotubes (CNT) under tension [9] by
relieving strain [10]. It can also act as a source for dislocation dipoles [4, 11].
Regardless of being the lowest energy defect in graphene [12] and other sp2 -hybridized
carbon nanostructures, the SW (55-77) needs about 5 eV to appear in graphene [3, 6],
and 3–5 eV in CNTs with a diameter above 1 nm [4, 5], which should lead to a negligible
equilibrium concentration of such defects at room temperature. However, recent advances
in HRTEM have allowed the identification of such defects in graphene [13, 14] and CNTs [8].
Moreover, SW transformations play an important role in the response of graphene to electron
irradiation [14, 15], leading to changes in the morphology of vacancy-type defects [16] and to
their migration. Such changes are equally surprising, because the barrier for bond rotation is
about 5 eV [6, 17], which should exclude thermal activation as a cause for SW transformation
at room temperature during experimentally relevant time scales. Regarding irradiation
effects, previous simulations [18] showed that an energy of ∼ 30 eV must be transferred
to a C atom in graphene in the in-plane direction for a bond rotation to occur. Also this
cannot explain the frequently observed SW transformations under the usual TEM imaging
conditions, since with typical acceleration voltages (. 300 kV) the transferred kinetic energy
in the direction almost perpendicular to the electron beam will remain significantly below
10 eV.
Here, by combining aberration-corrected (AC-) HRTEM with atomistic computer simulations, we show that topological defects associated with the SW transformation can be formed
in sp2 -hybridized carbon nanostructures by impacts of individual electrons at energies even
below the threshold for a carbon atom displacement. We further study in detail the mechanisms of irradiation-driven bond rotations for pre-existing vacancy-type defect structures
and how they transform and migrate due to electron impacts. At the same time we explain
2
why electron irradiation at moderate energies (∼ 100 keV) tends to rather amorphize [14]
than perforate graphene. We also show via simulations that the SW (55-77) can appear
in curved graphitic structures due to “incomplete” recombination of irradiation-induced
Frenkel defects, reminiscent of the formation of Wigner-type defects in silicon [19].
METHODS
Experimental Methods
Graphene membranes used in our experiments were prepared by mechanical exfoliation of
graphite on Si/SiO2 substrates and transfer of the resulting graphene sheets onto TEM grids
as described previously [20]. For TEM experiments we used an FEI Titan 80−300 equipped
with an image-side aberration corrector, operated at 80 kV. The spherical aberration was
set to 15 µm and images were recorded at Scherzer defocus. The extraction voltage of the
source was reduced to 2 kV and the condensor lens C2 was switched off in order to minimize
the energy spread. Under these conditions, dark contrast in the images can be directly
interpreted in terms of the atomic structure. Image sequences were recorded on a CCD
camera with exposure times of 1 s and intervals of approximately 2 s.
Computational Methods
We carried out atomistic computer simulations based on the spin-polarized density functional theory (DFT) implemented in the plane wave basis set vasp code [21]. The projector
augmented wave potentials [22] were used to describe the core electrons, and the generalized
gradient approximation of Perdew, Burke and Ernzernhof [23] for exchange and correlation. We included plane waves up to a kinetic energy of 300 eV. The k-point sampling of
the Brillouin zone was performed using the scheme of Monkhorst-Pack [24] for the periodic dimensions. Structure relaxation calculations were combined with molecular dynamics
(DFT-MD) simulations with a lower kinetic energy threshold and fewer k-points.
Due to the high computational cost of the DFT-MD method, only a few simulations were
carried out at this level. Whenever statistics needed to be gathered, we calculated the forces
using the non-orthogonal DFT-based tight binding (DFTB) method [25]. The main results
were checked against DFT-MD. In total, we carried out ∼ 27.000 dynamical DFTB-MD
3
simulations. The simulated structures consisted of 120–200 atoms and were fully optimized.
For the displacement threshold simulations, one of the atoms was assigned a kinetic energy
T with the initial velocity vector pointing to a preselected direction. The initial temperature
of the system was set to 5 K, although we observed no differences when carrying out the
simulations for initially stationary atoms. Displacement threshold Td (minimum kinetic
energy required to eject the atom) was found to be 22.50 eV, in a good agreement with
earlier DFTB results [26, 27]. It is also close to the DFT value (22.03 eV) [28]. For the
annihilation simulations, various system temperatures were studied (500–1200 K) both to
fasten the migration of the adatoms and to understand the effect of an elevated temperature
(as will be mentioned below).
RESULTS AND DISCUSSION
Stone-Wales Defects due to Single Electron Impacts
We begin the presentation of our results with the description of experimental observations
of SW (55-77) in HRTEM images. Several long image sequences, typically containing
hundreds of images from clean and initially defect-free graphene membranes, were recorded
at 80 kV. Occasionally, SW (55-77) defects appear in individual exposures, as in the example
shown in Figure 1b (Figure 1c with structure overlay). Remarkably, in most of the observed
cases, isolated SW (55-77) appeared in pristine graphene for one 1 s exposure, only to
disappear in the following frame. Hence, the lifetime of this defect under the 80 kV electron
beam in terms of irradiation dose is of the order of 107 e− /nm2 , the dose used for a single
exposure.
To understand the appearance and disappearance of SW (55-77) , we carried out atomistic
simulations for individual displacement events under the electron beam. After calculating
Td , we extended this calculation to all in-plane (φ) and out-of-plane angles in the range
θ ∈ [0◦ , 25◦], Figure 2. Displacements with θ > 25◦ would result in transferred kinetic
energies of more than 2 eV below Td for an electron beam even at 120 keV. Since the
displacement threshold increases for increasing θ, it is unlikely that this restriction would
lead us to miss any significant electron beam-induced structural changes, especially for
electron energies similar to those used in this study (80 keV). The calculated displacement
4
a
b
c
d
Figure 1: HRTEM images of an initially pristine graphene sample (a), followed by a SW (55-77)
(b). Frame (c) shows the same image with a structure overlay.
SW (55-77) disappears in the
following frame (d). Scale bar is 1 nm. (See also video S1 in Ref. [29].)
Relative threshold energy
60°
1.10
ϕ
1.08
45°
1.06
1.04
30°
1.02
ϕ
ϑ
SW(55-77) created
0°
0°
5°
10°
15°
20°
25°
ϑ
Figure 2: Relative displacement threshold Tdφ,θ /Tdθ=0 as a function of the displacement space angle
(θ, φ). Crosses mark the angles for which we observed impact-induced SW transformations.
thresholds are shown in Figure 2 in a relative scale along with the space angles for which
we observed the formation of SW (55-77) .
It is evident from Figure 2 that the SW transformation is a very likely event at displacement angles slightly away from the graphene plane normal (θ ≥ 2.5◦ ). The transferred
kinetic energies (T ) required for this process are below the displacement threshold for the
corresponding φ and θ (Tdφ,θ ) since no actual removal of the recoil atom is required for the
5
Displacement
ϕ
ϕ
9 5
5 7
7 5
ϑ: 12.5°
ϕ: 30°
T: 21.2 eV
9 5
5 7
7 5
ϑ: 6°
ϕ: 30°
T: 21.6 eV
Figure 3: Typical SW transformation processes for low out-of-plane displacement angles (θ ≤ 25◦ ):
(a) The “circle” process, and (b) the “nudge” process. The black spheres correspond to the recoil
(displaced) atoms, and the mono-vacancy structure is highlighted with numbered carbon rings in
the second panels. The last panel shows schematically the route of the displaced atom. (See also
videos S8 and S9 in Ref. [29].)
bond rotation to occur. Typically, T ≈ Tdφ,θ − 1 eV resulted in the SW (55-77) formation, although for some space angles even T ≈ Tdφ,θ − 2 eV was enough. The probability
for SW (55-77) formation is particularly high for certain space angles, which is related to
different mechanisms of SW transformation, as described below.
The above-presented result is in clear contrast with the earlier simulation results for
graphite [18] where no SW (55-77) formation was observed for low θ. This discrepancy is
caused by the neighboring graphene planes in the case of graphite: The displaced atom gets
attached to the adjacent layer and does not therefore initiate a bond rotation. In Figure 3 we
show the two processes which account for the majority of the SW transformations observed
during our simulations. In the “circle” process (Figure 3a), the displaced atom circles around
its neighbor, whereas in the “nudge” process (Figure 3b) it nudges the neighbor to cause the
bond rotation. Note that the example cases are for the same φ and almost same T , but for
different θ. The resulting process for each displacement is an interplay of all three variables
(T , θ, φ).
Similar mechanisms also exist for CNTs. However, as two new parameters (tube diameter
and chirality) should be introduced for quantitative analysis of SW transformation, we did
not study this process in nanotubes at length due to unreasonably high computational cost.
6
a
b
Figure 4: Experimental image sequence of a vacancy annihilation. Panel (a) shows a V1 (5-9) monovacancy, with an overlay in the inset, and (b) shows the same region in a later exposure, with no
defect visible. Scale bar is 1 nm. (See also video S7 in Ref. [29].)
Annihilation of Vacancies and Adatoms with Possible Formation of SW (55-77)
Since we frequently observed formation of vacancy–adatom pairs (adatoms play the role of
interstitials in graphene and CNTs) in our simulations of electron impacts onto graphene and
earlier in CNTs [26], we also explored another possible mechanism of SW (55-77) formation,
which is based on “incomplete” annihilation of a Frenkel defect. This study was motivated
by the peculiarities of the recombination of such a defect in bulk silicon. In that covalentlybonded material the recombination can give rise to either annihilation of the defect and
restoration of the perfect crystal lattice, or to a Wigner-type defect [19]. Such topological
defects are imperfections in the crystal lattice with the locally “correct” number of atoms
(as opposed to vacancies and interstitials), with the atomic configuration separated from
the perfect structure by a finite potential barrier. Such defects are deemed to also exist in
graphite [30, 31], and SW (55-77) can clearly be classified into this group.
While carbon adatoms on graphene [32, 33] and CNTs [34, 35] (especially those inside
nanotubes) are mobile at room temperature, they can easily find vacancies in the system
and annihilate. Indeed, we occasionally observed disappearance of vacancies in HRTEM
image sequences. Figure 4 shows an example of a mono-vacancy that disappears during
observation. This proves that mobile carbon atoms are present under our experimental
conditions, and may recombine with vacancy-type defects. However, we never noticed the
creation of a SW (55-77) after an observed mono-vacancy.
7
Annihilation
7
9
5
5
7
5
5
9
8
5
Figure 5: (a) SW (55-77) creation process as an outcome of a vacancy–adatom annihilation in a
(6,6) armchair nanotube. (b) Carbon dimer sputtering process in a (10,0) zigzag nanotube. The
black spheres stand for adatoms and the gray ones denote the dimer to be sputtered [marked in
the first panel of (b)]. The nanotube structures are shown from inside the tube. The tube axis is
in the horizontal direction. (See also videos S10 and S11 in Ref. [29].)
To simulate the annihilation process, we created nearby Frenkel defects (separated by
a few Ångströms) in a graphene layer and small (6,6) zigzag and (10,0) armchair CNTs
(radii r ≈ 4.1 Å and ∼ 3.9 Å, respectively). We then heated the structures and collected
statistics on the evolution of each system by running dynamical atomistic simulations at
various temperatures (500–1200 K). Three possible outcomes emerged from the simulations:
(1) perfect annihilation to the pristine structure (similar to the experimental images in
Figure 4), (2) formation of a SW (55-77) (Figure 5a) and, surprisingly, (3) sputtering of a
C2 dimer with a remaining reconstructed di-vacancy V2 (5-8-5) (Figure 5b).
For graphene, we always observed perfect annihilation in accordance with the experiments. However, if the adatom in graphene was placed on top of one of those two bonds
in the nine-membered carbon ring which are right next to the pentagon, a SW (55-77)
was spontaneously formed without an energy barrier. Thus SW (55-77) may also form in
graphene due to recombination of Frenkel defects, but the probability for this process must
be much lower than in CNTs.
For (10,0) CNT, we obtained perfect annihilation in approximately 54% of cases, SW (5577) was formed in approximately 34% of cases and dimer sputtering occurred in approximately 12% of cases. For (6,6) CNT, the values were 53%, 42% and 4%, respectively.
With increasing temperature the probability to sputter a dimer showed a slight tendency
to increase. We also ran the calculations for a (8,8) CNT (r ≈ 5.4 Å) at 800 K in order to
8
estimate the curvature dependency of the results. The values did not significantly differ from
those for the (6,6) CNT, except for a somewhat increased tendency to perfect annihilation
and decreased sputtering (with probabilities of 56% and 2%, respectively).
The formation energy of a spatially separated Frenkel defect in graphene within the
DFTB model is approximately 11.1 eV. Transformation to form SW defect from this initial
setup leads to an energy gain of 5.4 eV. By sputtering a dimer, graphene would instead gain
0.4 eV. Hence, all three observed outcomes are energetically reasonable also for graphene.
For nanotubes, the formation energies of both SW (55-77) and a di-vacancy are lowered due
to the curvature and stronger C-C bonds at pentagons [35]. The corresponding energy gains
are also higher, which can explain why the probability for SW (55-77) defect formation
and dimer evolution is higher in curved carbon nanostructures. The actual energies depend
on the local curvature. It is also plausible that the dimer evaporation process plays a role
in shrinking fullerenes under electron irradiation [36].
Stone-Wales Transformations in Vacancy-Type Defects
SW (55-77) defect represents the elementary case of a topological change in the graphene
structure, i.e., a single bond rotation in the otherwise perfect structure. More abundant,
however, are changes in the atomic configuration through bond rotations in the reconstructed
vacancy defects, as recent experiments indicate [14, 15]. In presence of a (multi-)vacancy, the
atomic configuration of a defect can be transformed between different metastable structures
via bond rotations. It was presumed [14] that such SW transformations of vacancy-type
defects would be stimulated by electron impacts, but the actual atomistic mechanism has
not been hitherto unraveled.
In order to get microscopic insight into irradiation-stimulated bond rotations near
vacancy-type defects in graphene, we carried out a set of experiments and dedicated simulations aimed at assessing the probability of SW transformations in the defect structures.
In our experiments, we initially generated “defective” graphene by brief 150 kV electron
irradiation, and then recorded image sequences of di-vacancies using 80 kV AC-HRTEM. As
observed previously [14], the vacancies can transform between different configurations under
the influence of the 80 kV electron beam. Moreover, the di-vacancies migrate and transition
between the different reconstructed configurations via SW transformations, typically until
9
a
b
c
d
Figure 6: Experimental images illustrating SW transformations in the atomic structure of divacancies and migration of these defects. Panels (a-d) show sequential HRTEM images of the same
defect, in the lower row with structure overlay. (a) V2 (5555-6-7777) transforming into V2 (555-777)
(b), and V2 (5-8-5) (c). Each of these transitions can be explained by a single bond rotation. In a
later frame (d), the defect is again a V2 (5555-6-7777), but shifted by one lattice parameter. Scale
bar is 1 nm. (See also videos S2–4 in Ref. [29].)
they cluster into larger defects (see videos S5 and S6 in Ref. [29]).
In Figure 6, we show an example of a di-vacancy defect that transforms between different
reconstructed shapes [V2 (5-8-5), V2 (555-777), V2 (5555-6-7777)] under the electron beam. The
changes in the atomic structure of these di-vacancy configurations can be described by SW
transformations at the defect. Moreover, multiple transformations allow migration of the
di-vacancy. Similar to SW (55-77) defect formation discussed above, the activation energy
for these transitions is far too high to allow a thermally activated process with an observable
rate at room temperature. Hence, on the basis of the observations, the activation energy for
the transition must be provided by the electron beam.
In order to confirm that the SW transformations at these defects are caused by single
electron impacts, we carried out atomistic simulations of such impacts onto atoms near divacancies. Due to the computational cost related to many non-equivalent atoms present
in the system and a large number of possible atomic configurations, we could not repeat
the detailed analysis of the role of initial space angle of the displacement similar to pristine
graphene, and therefore limited our simulations to the θ = 0◦ case for all non-equivalent
atoms at reconstructed di-vacancy structures. However, because the defects break the sym10
5
8
7 5
5
7
7 5
5
7 5
5
7
7 5
7 5 7
5
7 5 7
5
Figure 7: Changes in the atomic structures of reconstructed di-vacancy defects through SW transformations upon atom displacements perpendicular to the graphene sheet (θ = 0◦ ); (a) V2 (5-85)→ V2 (555-777) and (b) V2 (555-777)→ V2 (5555-6-7777). The black spheres indicate the recoil
(displaced) atoms. Structurally equivalent atoms to the displaced ones are marked with circles in
first panels. In the last panels, the circled atoms are those which could cause a backward transformation upon displacement in addition to the recoil atom. (See also videos S12–14 in Ref. [29].)
metry of the lattice, a directional preference for the displacements arises (atom displaced
in the perpendicular direction will change its direction due to local strain). This effect is
strong enough to facilitate bond rotations in reconstructed di-vacancy structures. In Figure 7 we present examples of such processes for both V2 (5-8-5)→ V2 (555-777) and V2 (555777)→ V2 (5555-6-7777) transformations.
Within these simulations, we never observed V2 (555-777)→ V2 (5-8-5) transformations.
Because the symmetry of the V2 (555-777) defect around the middle atom is the same as
that of pristine graphene (the middle atom is represented as a black sphere in the last
panel of Figure 7a), one would also expect that a θ 6= 0◦ displacement is required for this
transformation, similar to the SW (55-77) case, although in principle the surrounding
atoms could also cause this transformation. We noticed that the most likely di-vacancy
transformation, at least for θ = 0◦ , is the V2 (5555-6-7777)→ V2 (555-777). The V2 (5-85)→ V2 (555-777) was the least likely one of those observed. We never observed a V2 (5-8-5)→
2×(5-7) transformation [14] during our simulations, which we also attribute to the limited
simulated conditions (θ = 0◦ ). Curiously, however, we did observe one transformation in
which a V2 (5-8-5) di-vacancy directly migrated one step in the zigzag lattice direction.
11
Another interesting observation originating from these simulations is the fact that the
displacement threshold for atoms in the central part of the reconstructed defects [V2 (555777) and V2 (5555-6-7777)] are higher than that for pristine graphene (by as much as 5%).
This may explain why defect structures tend to grow into larger and larger amorphous
patches instead of collapsing into holes under continuous electron irradiation at low voltages
(≤ 100 keV) [14]: Even when atoms are removed from the defected area, the displacements
occur at the edges of the existing defects rather than at the central part where the local
atomic density is already lower. Clearly, since the core structure of these defects consists of
carbon hexagons, there must exist a limiting size above which the displacement threshold
becomes similar to that of ideal graphene.
CONCLUSIONS
To conclude, by combining AC-HRTEM experiments and atomistic simulations, we have
shown that the bond rotations which lead to creation of topological defects in carbon nanostructures are caused by single electron impacts or incomplete annihilation of Frenkel defects.
This explains the discrepancy between experimental observations of Stone-Wales defects and
their relatively high formation energy and even higher energy barrier for bond rotation.
The SW transformation in graphene can be initiated at least in two different ways upon
electron impact (involving a “circling” or “nudging” motion), and for almost any space angle, provided that enough energy is transferred from the electron to the target atom. Our
simulations indicate that SW (55-77) can appear as a result of “incomplete” recombination of a Frenkel defect reminiscent of the formation of Wigner-type defects in silicon [19].
However, this is much more likely in the case of local curvature, as in nanotubes. More
surprisingly, we also observed sputtering of C2 dimer as a result of annihilation of a Frenkel
defect in carbon structures with high curvature.
Moreover, we showed that the displacement threshold of atoms in the central area of
reconstructed defects is higher than that of pristine graphene, by as much as 5%, which
explains why defected graphene under low-energy electron irradiation (. 100 keV) tends
to turn graphene into a two-dimensional amorphous structure [14] instead of a perforated
membrane. For different di-vacancy structures, even displacements in direction perpendicular to the graphene layer can initiate SW transformation and thus local structural changes
12
and defect migration.
Our results provide microscopic insight into the irradiation-induced changes in the atomic
structure of carbon nanosystems under electron irradiation, and taking into account the
interesting electronic properties of defects associated with SW transformations [14, 37, 38],
may open new avenues for irradiation-mediated engineering [39] of carbon nanostructures
with next-generation electron microscopes.
ACKNOWLEDGMENTS
We acknowledge financial support by the German Research Foundation (DFG), the German Ministry of Science, Research and Arts (MWK) of the state Baden-Wuerttemberg
within the SALVE (sub angstrom low voltage electron microscopy) project and Academy of
Finland through several projects. We also thank CSC, Espoo, Finland, for generous grants
on computer time.
[1] Thrower, P. The study of defects in graphite by transmission electron microscopy. In Chemistry
and physics of carbon, vol. 5 ; Edited by P. L. Walker Jr., Marcel Dekker: New York, 1969.
[2] A. J. Stone and D. J. Wales, Chem. Phys. Lett. 128, 501 (1986).
[3] J. Ma, D. Alfé, A. Michaelides, E. Wang, Phys. Rev. B 80, 033407 (2009).
[4] E. Ertekin, D. C. Chrzan, and M. S. Daw, Phys. Rev. B 79, 155421 (2009).
[5] L. G. Zhou and S. Shi, Appl. Phys. Lett. 83, 1222 (2003).
[6] L. Li, S. Reich, and J. Robertson, Phys. Rev. B 72, 184109 (2005).
[7] Y. Miyamoto, A. Rubio, S. Berber, M. Yoon, and D. Tománek, Phys. Rev. B 69, 121413+
(2004).
[8] K. Suenaga, H. Wakabayashi, M. Koshino, Y. Sato, K. Urita, and S. Iijima, Nat. Nanotech.
2, 358 (2007).
[9] J. Y. Huang, S. Chen, Z. Q. Wang, K. Kempa, Y. M. Wang, S. H. Jo, G. Chen, M. S.
Dresselhaus, and Z. F. Ren, Nature 439, 281 (2006).
[10] M. B. Buogiorno Nardelli, B. I. Yakobson, and J. Bernholc, Phys. Rev. B 57, R4277 (1998).
[11] F. Ding, K. Jiao, M. Wu, and B. I. Yakobson, Phys. Rev. Lett. 98, 075503+ (2007).
13
[12] F. Banhart, J. Kotakoski, and A. V. Krasheninnikov, ACS Nano 5, 26 (2011).
[13] J. C. Meyer, C. Kisielowski, R. Erni, M. D. Rossell, M. F. Crommie, and A. Zettl, Nano
Letters 8, 3582 (2008).
[14] J. Kotakoski, A. V. Krasheninnikov, U. Kaiser, and J. C. Meyer, Phys. Rev. Lett. 106,
105505+ (2011).
[15] B. Song, G. F. Schneider, Q. Xu, G. Pandraud, C. Dekker, and H. Zandbergen, (2011).
[16] G. D. Lee, C. Z. Wang, E. Yoon, N. M. Hwang, D. Y. Kim, and K. M. Ho, Phys. Rev. Lett.
95, 205501 (2005).
[17] O. Cretu, A. V. Krasheninnikov, J. A. Rodrı́guez-Manzo, L. Sun, R. M. Nieminen, and F.
Banhart, Phys. Rev. Lett. 105, 196102+ (2010).
[18] O. V. Yazyev, I. Tavernelli, U. Rothlisberger, and L. Helm, Phys. Rev. B 75, 115418 (2007).
[19] M. Tang, L. Colombo, J. Zhu, and T. Diaz de la Rubia, Phys. Rev. B 55, 14279 (1997).
[20] J. C. Meyer, C. O. Girit, M. F. Crommie, and A. Zettl, Appl. Phys. Lett. 92, 123110 (2008).
[21] G. Kresse and J. Furthmüller, Phys. Rev. B 54, 11169 (1996).
[22] P. E. Blöchl, Phys. Rev. B 50, 17953 (1994).
[23] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996).
[24] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 (1976).
[25] T. Frauenheim, G. Seifert, M. Elstner, T. Nichaus, C. Köhler, M. Amkreutz, M. Sternberg,
Z. Hajnal, A. Di Carlo, and S. Suhai, J. Phys.: Condens. Matter 14, 3015 (2002).
[26] A. V. Krasheninnikov, F. Banhart, J. X. Li, A. S. Foster, and R. M. Nieminen, Phys. Rev. B
72, 125428 (2005).
[27] A. Zobelli, A. Gloter, C. P. Ewels, G. Seifert, and C. Colliex, Phys. Rev. B 75, 245402 (2007).
[28] J. Kotakoski, C. H. Jin, O. Lehtinen, K. Suenaga, and A. V. Krasheninnikov, Phys. Rev. B
82, 113404 (2010).
[29] See EPAPS Document No. to-be-inserted for experimental videos and additional images. For
more information on EPAPS, see http://www.aip.org/pubservs/epaps.html..
[30] C. P. Ewels, R. H. Telling, A. A. El-Barbary, M. I. Heggie, and P. R. Briddon, Phys. Rev.
Lett. 91, 025505 (2003).
[31] R. H. Telling, C. P. Ewels, A. A. El-Barbary, and M. I. Heggie, Nat. Mater. 2, 333 (2003).
[32] P. O. Lehtinen, A. S. Foster, A. Ayuela, A. V. Krasheninnikov, K. Nordlund, and R. M.
Nieminen, Phys. Rev. Lett. 91, 017202 (2003).
14
[33] R. Erni, M. D. Rossell, M. T. Nguyen, S. Blankenburg, D. Passerone, P. Hartel, N. Alem, K.
Erickson, W. Gannett, and A. Zettl, Phys. Rev. B 82, 165443+ (2010).
[34] Y. Gan, J. Kotakoski, A. V. Krasheninnikov, K. Nordlund, and F. Banhart, New J. Phys. 10,
023022 (2008).
[35] A. V. Krasheninnikov, P. O. Lehtinen, A. S. Foster, and R. M. Nieminen, Chem. Phys. Lett.
418, 132 (2006).
[36] J. Y. Huang, F. Ding, K. Jiao, and B. I. Yakobson, Phys. Rev. Lett. 99, 175503 (2007).
[37] A. Lherbier, Simon M.-M. Dubois, X. Declerck, S. Roche, Y. M. Niquet, and J. C. Charlier,
Phys. Rev. Lett. 106, 046803 (2011).
[38] D. J. Appelhans, L. D. Carr, and M. T. Lusk, New Journal of Physics 12, 125006 (2010).
[39] A. V. Krasheninnikov and K. Nordlund, J. Appl. Phys. 107, 071301 (2010).
15