# Simulated Nanoscale Peeling Process of Monolayer Graphene Sheet: Effect of Edge Structure and Lifting Position

Authors: Naruo Sasaki; Hideaki Okamoto; Shingen Masuda; Kouji Miura; Noriaki Itamura
Journal: Journal of Nanomaterials
 (2010)

---

## Abstract

The nanoscale peeling of the graphene sheet on the graphite surface is numerically studied by molecular mechanics simulation. For center-lifting case, the successive partial peelings of the graphene around the lifting center appear as discrete jumps in the force curve, which induce the arched deformation of the graphene sheet. For edge-lifting case, marked atomic-scale friction of the graphene sheet during the nanoscale peeling process is found. During the surface contact, the graphene sheet takes the atomic-scale sliding motion. The period of the peeling force curve during the surface contact decreases to the lattice period of the graphite. During the line contact, the graphene sheet also takes the stick-slip sliding motion. These findings indicate the possibility of not only the direct observation of the atomic-scale friction of the graphene sheet at the tip/surface interface but also the identification of the lattice orientation and the edge structure of the graphene sheet.


---

## Body

## 1. Introduction

Adhesion and peeling phenomena play important roles in connecting two objects regardless of whether they are inorganic, organic, or biological materials, which contributes to building up microscopic devices. The carbon nanostructures such as, carbon nanotube (CNT) and graphene have recently attracted great interests as the components of the electronic, magnetic, and optical devices. We have so far studied the peeling mechanics of the carbon nanotube (CNT) adsorbed onto the graphite surface both theoretically [1–3] and experimentally [4, 5]. It is clarified that the transition from the line- to the point-contact between the CNT and the graphite surface occurs during the peeling process [1–5]. The CNT on the sub-microscale has the same size as the spatulae of the microscopic hairs aligned on the gecko foot [6, 7]. Therefore the study of the peeling process of the nanoscale objects such as, CNT is useful for not only developing the gecko-foot-mimic adhesives [8] but also understanding the elementary process of adhesion.On the other hand, since the success of its experimental isolation [9], the potential of various applications of the graphene such as, the components of the electronic devices [10, 11] has been discussed by many researchers. There is also a possibility that its adhesion with the substrate is applied to the adhesive tape at nanoscale. Therefore the peeling mechanics of the graphene sheet is very important, which can be regarded as the elementary process of the macroscopic sticky tape such as, the gecko-foot-mimic adhesives [6–8] or that of the microscopic extension of the crack in the fracture process. In our preliminary experiments, we have already succeeded in peeling the multilayered graphene plate with a thickness of several μm by using atomic-force microscopy tip [12]. Here the two-component epoxy resin adhesive is used to bond the graphene plate to the AFM tip. Here the standard Si3N4 tip for the contact AFM experiment is used. The junction formed between the AFM tip and the graphene should be mechanically rigid enough to measure the elasticity of the graphene sheet during the peeling process. The two-component epoxy resin adhesive satisfies the above condition. If the thickness of the peeled graphene plate is reduced, the comparison between the present simulation and the experiment will become possible.Therefore, in this paper, ahead of experiment, we have theoretically reported the nanoscale peeling behaviors of the monolayer graphene sheet based on the molecular mechanics simulation [13, 14]. The peeling force curve exhibits the nanoscale change of the graphene shape from the surface to the line contact. The center position and the left edge are chosen as the lifting position. In Section 3, the peeling of the monolayer graphene sheet with the armchair edge for lifting the center position is discussed. In Sections 4 and 5, the peeling of the monolayer graphene sheet with the armchair- and zigzag-edge for lifting the edge position is discussed, respectively.


## 2. Model and Method of Simulation

In the simulation, a rectangular-shaped monolayer graphene sheet with each side of38 Å ×20-21 Å, comprised of 310 carbon atoms, is peeled from the rigid rectangular graphene sheet (which is called, the “graphite surface,” hereafter) with each side of 164-165 Å ×58 Å, comprised of 3536 carbon atoms [Figure 1(a)]. First, both the above graphene sheets are separately optimized by minimizing the covalent bonding energy described by the Tersoff potential energy [15], Vcov, using the Polak-Rebiere-type conjugate gradient (CG) method [16]. Here the convergence criterion is set so that the maximum of absolute value of all the forces acting on the movable atoms becomes lower than 10-5 eV/Å. Next, the graphene sheet is put and adsorbed onto the graphite surface, so that the AB stacking registry between the graphene sheet and the graphite surface is satisfied as shown in Figures 1(b) and 1(c). Here the green-colored six-membered ring at the center position or the outermost left edge of the graphene sheet is assumed to be attached to the AFM tip apex (Figure 1(a)), and then it is gradually moved upward along the z direction, parallel to the [0001] axis, by 0.1 Å. For each lifting position of the graphene sheet, z, the total energy Vtotal=Vcov+VvdW is minimized using the CG method, where VvdW is the nonbonding vdW interaction described by the modified Lennard-Jones (LJ) potential energy [17, 18], acting between the graphene sheet and the graphite surface. Thus the optimized positions of the movable carbon atoms of the graphene sheet, (x,y,z), the vertical peeling force Fz, and the lateral sliding forces Fx and Fy, acting on the lifting center, are calculated during the peeling process. In this paper, the graphene sheets with armchair-(Figure 1(b)) and zigzag-edges (Figure 1(c)) are discussed.(a) The schematic illustration of the model of the monolayer graphene sheet physically adsorbed onto the rigid graphite surface used in the simulation. The green-colored six-membered ring at the center position or left edge of the graphene sheet is assumed to be adsorbed onto the atomic force microscopy tip apex indicated by broken lines, and it is moved upward along thez (or [0001]) direction, by z=0.1 Å. Initial AB stacking registry of the red-colored graphene sheet with (b) armchair and (c) zigzag edge adsorbed onto the blue-colored graphite surface within the x-y plane.
(a)(b)(c)


## 3. Center-Lifting Case of Armchair-Edge Graphene

When the six-membered ring located at the center position of the monolayer graphene sheet is lifted, the graphene sheet exhibits the characteristic transition of its shape during the peeling process within thex–z plane as illustrated in Figures 2(A)–2(J), corresponding to Figures 3(A)–3(J), the vertical force acting on the lifting center position Fz plotted as a function of the displacement from the initial position along z-direction, z. At first the monolayer graphene sheet takes an initial planar structure parallel to the rigid graphite surface (Figure 2(A): z=0 Å). Here the surface contact is formed between the graphene sheet and the graphite surface. The vertical force Fz is zero (Figure 3(a)). Just after the beginning of the peeling (Figure 2(B): z=2.0 Å), the attractive interaction force takes the minimum value, -3.1 eV/ Å (Figure 3(B)).The transition of the shape of the monolayer graphene sheet during the peeling process from A to J within thex-z plane. The red-colored graphene sheet and blue-colored graphite surface are shown. The displacement of the lifting center position from the initial position, z(Å), is indicated on the upper-right positions of each picture.Figure 3
The vertical force,Fz, acting on the center six-membered ring, plotted as a function of the lifting displacement z. The positions A–J correspond to those of Figure 2.Betweenz=2.0 Å and 2.1 Å, the first discrete partial peeling of the graphene occurs (Figures 2(B)→2(C)), which produces the 1st discontinuous jump in the force curve (Figures 3(B)→3(C)). The partial peeled area around the lifting center of the graphene is shown in Figures 4(B)→4(C). Then, between z=2.4 Å and 2.5 Å, the second discrete partial peeling of the graphene occurs (Figures 2(D)→2(E)), which produces the 2nd discontinuous jump in the force curve (Figures 3(D)→3(E)). The partial peeled area of the graphene is shown in Figures 4(D)→4(E). Which of these two areas is peeled first is expected to be actually the stochastic process under the room temperature condition. Now the surface contact region is split into the left and right sections (Figure 2(E)). After the two discrete jumps, Fz increases as the peeling proceeds, since the attractive surface contact region gradually decreases (Figure 2(F): z=5.0 Å). Then the surface contact continuously turns into the line contact at z=7.3 Å (Figure 2(G)). Here the 'line contact' is defined by the following two criteria: (1) The carbon atoms on the left and right outermost arrays of the graphene sheet (Figure 1(b)) receive the repulsive interaction force from the graphite surface. (2) The carbon atoms on the second arrays (Figure 1(b)) next to the outermost arrays receive the attractive interaction force. As illustrated in Figure 5, the average forces acting on one carbon atom on the outermost and the second arrays satisfy the above criteria at z=7.3 Å, which corresponds to Figure 2(G).Figure 4
The atomic structures of the graphene sheet just before and after the discrete change, B→C, D→E, and I→J. The regions surrounded by dotted ellipses show the partial peeled areas.Figure 5
The averaged forces acting on one atom on the left and right outermost arrays (red-colored) and those on the left and right second arrays (blue-colored), as a function of the displacement of the lifting center position from the initial position,z(Å).Once the line contacts are formed between the free edges (outermost arrays) of the peeled graphene sheet and the graphite surface, they clearly slide on the graphite surface as indicated by a circle in Figures2(G)→2(H)→2(I) with a rapid increase of the bending of the graphene sheet. Within x-y plane, the right outermost array of the graphene sheet slides nearly straightforward along -x direction, not so sensitive to the lattice structure of the surface as illustrated in Figures 6(a) and 6(b), which show the trajectories of the two carbon atoms on the right outermost array illustrated in Figure 1(b). The sliding of the outermost arrays during G and I appears much more clearly than that during A and G. During H and I, the decrease of Fz (Figures 3(H) and 3(I)) can be explained by the decrease of the repulsive force acting on the carbon atoms on the left and right edges of the graphene sheet as shown in Figure 5, that is to say, the relative increase of the effect of the attractive interaction force.The trajectory of the two carbon atoms on the right free edge (outermost array) indicated by white circles in Figure1(b). (a) The whole trajectory A→I, and (b) the part of the trajectory A→E, including the discrete jumps, B→C and D→E, are indicated. White circles mean carbon atoms of the graphite surface. The indices A–I correspond to those in Figure 2.
(a)(b)When the bending of the graphene sheet becomes larger than a certain range, both the left and right line contacts break and the graphene sheet is completely peeled from the surface (Figures2(I)→2(J): z=14.8Å→14.9 Å), which produces the 3rd discontinuous jump in the force curve (Figures 3(I)→3(J)). As illustrated in Figures 2(I)→2(J) and Figures 4(I)→4(J), the graphene sheet exhibits the transition from the arched shape to the planer shape.Thus the vertical peeling forceFz exhibits the characteristic shape as shown in Figure 3, which reflects the transition from the surface to the line contact between the graphene sheet and the graphite surface. On the other hand, the lateral sliding force Fx is zero due to the structural symmetry of the system. However, the lateral sliding force Fy shows a finite value with an oscillation whose period and amplitude decrease as z increases. This oscillation of Figure 7 reflects the trajectory of the graphene edges illustrated in Figure 6 at the graphene-substrate interface during the peeling process. The maximum lateral force Fy⋍0.1 eV/Å which is only about 3% of the absolute value of the maximum adhesion force |Fz|=3.1 eV/Å.Figure 7
The lateral force,Fy, acting on the center six-membered ring, plotted as a function of the displacement z. The positions A–J correspond to those of Figure 2.


## 4. Edge-Lifting Case of Armchair-Edge Graphene

### 4.1. Nanoscale Peeling along Vertical Direction

When the left edge of the monolayer graphene sheet with armchair edge (Figure1(b)) is lifted, the shape of the graphene sheet markedly changes during the peeling process within the x–z plane as illustrated in Figures 8(A)–8(J), corresponding to Figures 9(a)A–9(a)J, the vertical force acting on the lifting edge Fz plotted as a function of the edge height z.The transition of the shape of the monolayer graphene sheet with free edge of armchair-type during the peeling process from A to J within thex–z plane. The red-colored graphene sheet and blue-colored graphite surface are shown. The height of the lifting left edge z(Å) is indicated on the upper-right positions of each picture.(a) The vertical force,Fz, acting on the lifting edge, plotted as a function of the lifting edge height z for the graphene with armchair-type free edge. The indices A–J correspond to those of Figure 8. (b) The red-colored averaged force per one atom acting on the outermost array, and the blue-colored one acting on the second array, as a function of the lifting edge height z(Å). The indices E, F, H, and I correspond to those of Figures 8 and 9(a).
(a)(b)


#### 4.1.1. Surface-Contact Region

At first the monolayer graphene sheet takes an initial planar structure parallel to the rigid graphite surface (Figure8(A): z=0 Å), and the vertical force Fz is zero (Figure 9(a)A), which means the graphene sheet completely takes the surface contact with the graphite surface. Just after the start of the peeling (Figure 8(B): z=1.1 Å), the attractive interaction force |Fz| becomes the maximum, 0.74 eV/Å (Figure 9(a)B). After that the surface contact area gradually decreases as the peeling proceeds (Figures 8(C)–8(E)), where Figure 9(a) exhibits the atomic-scale zigzag structures (Figures 9(a)C–E), which will be explained in the Section 4.2.


#### 4.1.2. Line-Contact Region

After the surface contact vanishes, the line contact appears (Figure8(F)). Here the “line contact” is defined by the following two criteria similar to the case of the center-lifting peeling. (1) The outermost array of the free edge of the graphene sheet (Figure 1(b)) receives the averaged repulsive interaction force per one carbon atom from the graphite surface. (2) The second array next to the outermost array (Figure 1(b)) receives the averaged attractive interaction force per one carbon atom. As illustrated in Figure 9(b), the edge height z=31.3 Å satisfies the above criteria (Figure 8(F)). Here the free edge adsorbed onto the graphite surface is nearly fixed even if the left edge is lifted (Figure 8(F)→Figure 8(G): z=33.9 Å). As a result, the in-plane bending of the graphene sheet markedly decreases, and |Fz| decreases to zero (Figures 9(a)F–G). After that the free edge slides on the surface until the graphene sheet takes planar structure within the y-z plane (Figure 8(H): z=37.9 Å), where |Fz| increases again (Figure 9(a)H).


#### 4.1.3. Toward Complete Peeling

Then the graphene sheet is continuously moved upward (Figure8(I): z=38.0 Å → Figure 8(J): z=38.4 Å) and is completely peeled from the surface. It is noted that the line contact clearly vanishes at z=38.0 Å (Figure 9(b)I). The attractive interaction force |Fz| increases to take the maximum value (Figures 9(a)I→J), and then it gradually decreases to zero toward the complete peeling.


### 4.2. Atomic-Scale Sliding within Lateral Plane

Figure9(a) shows the atomic-scale zigzag structures within the surface- and line-contact regions, which can be explained by the following atomic-scale sliding motions of the graphene sheet within the x-y plane.


#### 4.2.1. Surface-Contact Region

During the surface contact region between C and E in Figure9(a), z-Fz curve takes the atomic-scale zigzag structures from I to VII. The zigzag behaviors exhibit the transition from the continuous (Figures 10(a)I-II) to the sawtooth shapes (Figures 10(a)III–VII). First Figures 10(a)1→2→3→4→5 correspond to Figures 10(b)1→2→3→4→5, which show that the graphene sheet continuously slides passing over the nearest neighboring AB stacking sites with the graphite surface. The trajectories of the graphene sheet exhibit the continuous zigzag paths as shown in Figure 10(b)5. Next Figures 10(a)6→7→8→9→10 correspond to Figures 10(c)6→7→8→9→10 which show that the graphene sheet takes the zigzag stick-slip motions between the nearest neighboring AB stacking sites. Just before the slip, the graphene sheet deviates quite a little from the AB-stacking site (Figure 10(c)6: z=17.3 Å). Then it discretely jumps or slips to the neighboring AB-stacking site (Figure 10(c)7: z=17.4 Å). As the peeling proceeds, the graphene sheet continuously slides quite a little (Figure 10(c)8: z=18.2 Å), then it discretely slips again to the neighboring AB stacking site (Figure 10(c)9: z=18.3 Å). After that the graphene sheet continuously slides quite a little again (Figure 10(c)10: z=20.1 Å) until the next slip toward the neighboring AB stacking site occurs. As a result the trajectories of the graphene sheet exhibit the discrete zigzag paths connecting the nearest neighboring AB-stacking sites as shown in Figure 10(c)10. The period of the zigzag behavior of the Fz curve decreases from 3.7 Å to 2.5 Å as shown in Figure 10(a) as the peeling proceeds. The lattice spacing of the graphite surface, 2.5 Å, appears in the peeling force curve particularly for the stick-slip region.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the continuous and stick-slip process during the surface contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 5 indicated in (a). (c) The trajectories of the two carbon atoms on the free edge from 6 to 10 indicated in (a).
(a)(b)(c)


#### 4.2.2. Line-Contact Region

During the line contact region between G and H in Figure9(a), z-Fz curve takes another atomic-scale zigzag structures as shown in Figure 11(a). One of the zigzag behaviors in the force curve (Figures 11(a)1→2→3) corresponds to the stick-slip sliding motion of the graphene sheet (Figures 11(b)1→2→3). Here the free edge of the graphene sheet slides with nearly the straight stick-slip motions. One of the carbon atoms on the free edge passes over the carbon-carbon bonds as shown in Figures 11(b)1 and 3.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the nearly straight stick-slip region during the line contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 3 indicated in (a).
(a)(b)


## 4.1. Nanoscale Peeling along Vertical Direction

When the left edge of the monolayer graphene sheet with armchair edge (Figure1(b)) is lifted, the shape of the graphene sheet markedly changes during the peeling process within the x–z plane as illustrated in Figures 8(A)–8(J), corresponding to Figures 9(a)A–9(a)J, the vertical force acting on the lifting edge Fz plotted as a function of the edge height z.The transition of the shape of the monolayer graphene sheet with free edge of armchair-type during the peeling process from A to J within thex–z plane. The red-colored graphene sheet and blue-colored graphite surface are shown. The height of the lifting left edge z(Å) is indicated on the upper-right positions of each picture.(a) The vertical force,Fz, acting on the lifting edge, plotted as a function of the lifting edge height z for the graphene with armchair-type free edge. The indices A–J correspond to those of Figure 8. (b) The red-colored averaged force per one atom acting on the outermost array, and the blue-colored one acting on the second array, as a function of the lifting edge height z(Å). The indices E, F, H, and I correspond to those of Figures 8 and 9(a).
(a)(b)


### 4.1.1. Surface-Contact Region

At first the monolayer graphene sheet takes an initial planar structure parallel to the rigid graphite surface (Figure8(A): z=0 Å), and the vertical force Fz is zero (Figure 9(a)A), which means the graphene sheet completely takes the surface contact with the graphite surface. Just after the start of the peeling (Figure 8(B): z=1.1 Å), the attractive interaction force |Fz| becomes the maximum, 0.74 eV/Å (Figure 9(a)B). After that the surface contact area gradually decreases as the peeling proceeds (Figures 8(C)–8(E)), where Figure 9(a) exhibits the atomic-scale zigzag structures (Figures 9(a)C–E), which will be explained in the Section 4.2.


### 4.1.2. Line-Contact Region

After the surface contact vanishes, the line contact appears (Figure8(F)). Here the “line contact” is defined by the following two criteria similar to the case of the center-lifting peeling. (1) The outermost array of the free edge of the graphene sheet (Figure 1(b)) receives the averaged repulsive interaction force per one carbon atom from the graphite surface. (2) The second array next to the outermost array (Figure 1(b)) receives the averaged attractive interaction force per one carbon atom. As illustrated in Figure 9(b), the edge height z=31.3 Å satisfies the above criteria (Figure 8(F)). Here the free edge adsorbed onto the graphite surface is nearly fixed even if the left edge is lifted (Figure 8(F)→Figure 8(G): z=33.9 Å). As a result, the in-plane bending of the graphene sheet markedly decreases, and |Fz| decreases to zero (Figures 9(a)F–G). After that the free edge slides on the surface until the graphene sheet takes planar structure within the y-z plane (Figure 8(H): z=37.9 Å), where |Fz| increases again (Figure 9(a)H).


### 4.1.3. Toward Complete Peeling

Then the graphene sheet is continuously moved upward (Figure8(I): z=38.0 Å → Figure 8(J): z=38.4 Å) and is completely peeled from the surface. It is noted that the line contact clearly vanishes at z=38.0 Å (Figure 9(b)I). The attractive interaction force |Fz| increases to take the maximum value (Figures 9(a)I→J), and then it gradually decreases to zero toward the complete peeling.


## 4.1.1. Surface-Contact Region

At first the monolayer graphene sheet takes an initial planar structure parallel to the rigid graphite surface (Figure8(A): z=0 Å), and the vertical force Fz is zero (Figure 9(a)A), which means the graphene sheet completely takes the surface contact with the graphite surface. Just after the start of the peeling (Figure 8(B): z=1.1 Å), the attractive interaction force |Fz| becomes the maximum, 0.74 eV/Å (Figure 9(a)B). After that the surface contact area gradually decreases as the peeling proceeds (Figures 8(C)–8(E)), where Figure 9(a) exhibits the atomic-scale zigzag structures (Figures 9(a)C–E), which will be explained in the Section 4.2.


## 4.1.2. Line-Contact Region

After the surface contact vanishes, the line contact appears (Figure8(F)). Here the “line contact” is defined by the following two criteria similar to the case of the center-lifting peeling. (1) The outermost array of the free edge of the graphene sheet (Figure 1(b)) receives the averaged repulsive interaction force per one carbon atom from the graphite surface. (2) The second array next to the outermost array (Figure 1(b)) receives the averaged attractive interaction force per one carbon atom. As illustrated in Figure 9(b), the edge height z=31.3 Å satisfies the above criteria (Figure 8(F)). Here the free edge adsorbed onto the graphite surface is nearly fixed even if the left edge is lifted (Figure 8(F)→Figure 8(G): z=33.9 Å). As a result, the in-plane bending of the graphene sheet markedly decreases, and |Fz| decreases to zero (Figures 9(a)F–G). After that the free edge slides on the surface until the graphene sheet takes planar structure within the y-z plane (Figure 8(H): z=37.9 Å), where |Fz| increases again (Figure 9(a)H).


## 4.1.3. Toward Complete Peeling

Then the graphene sheet is continuously moved upward (Figure8(I): z=38.0 Å → Figure 8(J): z=38.4 Å) and is completely peeled from the surface. It is noted that the line contact clearly vanishes at z=38.0 Å (Figure 9(b)I). The attractive interaction force |Fz| increases to take the maximum value (Figures 9(a)I→J), and then it gradually decreases to zero toward the complete peeling.


## 4.2. Atomic-Scale Sliding within Lateral Plane

Figure9(a) shows the atomic-scale zigzag structures within the surface- and line-contact regions, which can be explained by the following atomic-scale sliding motions of the graphene sheet within the x-y plane.


### 4.2.1. Surface-Contact Region

During the surface contact region between C and E in Figure9(a), z-Fz curve takes the atomic-scale zigzag structures from I to VII. The zigzag behaviors exhibit the transition from the continuous (Figures 10(a)I-II) to the sawtooth shapes (Figures 10(a)III–VII). First Figures 10(a)1→2→3→4→5 correspond to Figures 10(b)1→2→3→4→5, which show that the graphene sheet continuously slides passing over the nearest neighboring AB stacking sites with the graphite surface. The trajectories of the graphene sheet exhibit the continuous zigzag paths as shown in Figure 10(b)5. Next Figures 10(a)6→7→8→9→10 correspond to Figures 10(c)6→7→8→9→10 which show that the graphene sheet takes the zigzag stick-slip motions between the nearest neighboring AB stacking sites. Just before the slip, the graphene sheet deviates quite a little from the AB-stacking site (Figure 10(c)6: z=17.3 Å). Then it discretely jumps or slips to the neighboring AB-stacking site (Figure 10(c)7: z=17.4 Å). As the peeling proceeds, the graphene sheet continuously slides quite a little (Figure 10(c)8: z=18.2 Å), then it discretely slips again to the neighboring AB stacking site (Figure 10(c)9: z=18.3 Å). After that the graphene sheet continuously slides quite a little again (Figure 10(c)10: z=20.1 Å) until the next slip toward the neighboring AB stacking site occurs. As a result the trajectories of the graphene sheet exhibit the discrete zigzag paths connecting the nearest neighboring AB-stacking sites as shown in Figure 10(c)10. The period of the zigzag behavior of the Fz curve decreases from 3.7 Å to 2.5 Å as shown in Figure 10(a) as the peeling proceeds. The lattice spacing of the graphite surface, 2.5 Å, appears in the peeling force curve particularly for the stick-slip region.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the continuous and stick-slip process during the surface contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 5 indicated in (a). (c) The trajectories of the two carbon atoms on the free edge from 6 to 10 indicated in (a).
(a)(b)(c)


### 4.2.2. Line-Contact Region

During the line contact region between G and H in Figure9(a), z-Fz curve takes another atomic-scale zigzag structures as shown in Figure 11(a). One of the zigzag behaviors in the force curve (Figures 11(a)1→2→3) corresponds to the stick-slip sliding motion of the graphene sheet (Figures 11(b)1→2→3). Here the free edge of the graphene sheet slides with nearly the straight stick-slip motions. One of the carbon atoms on the free edge passes over the carbon-carbon bonds as shown in Figures 11(b)1 and 3.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the nearly straight stick-slip region during the line contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 3 indicated in (a).
(a)(b)


## 4.2.1. Surface-Contact Region

During the surface contact region between C and E in Figure9(a), z-Fz curve takes the atomic-scale zigzag structures from I to VII. The zigzag behaviors exhibit the transition from the continuous (Figures 10(a)I-II) to the sawtooth shapes (Figures 10(a)III–VII). First Figures 10(a)1→2→3→4→5 correspond to Figures 10(b)1→2→3→4→5, which show that the graphene sheet continuously slides passing over the nearest neighboring AB stacking sites with the graphite surface. The trajectories of the graphene sheet exhibit the continuous zigzag paths as shown in Figure 10(b)5. Next Figures 10(a)6→7→8→9→10 correspond to Figures 10(c)6→7→8→9→10 which show that the graphene sheet takes the zigzag stick-slip motions between the nearest neighboring AB stacking sites. Just before the slip, the graphene sheet deviates quite a little from the AB-stacking site (Figure 10(c)6: z=17.3 Å). Then it discretely jumps or slips to the neighboring AB-stacking site (Figure 10(c)7: z=17.4 Å). As the peeling proceeds, the graphene sheet continuously slides quite a little (Figure 10(c)8: z=18.2 Å), then it discretely slips again to the neighboring AB stacking site (Figure 10(c)9: z=18.3 Å). After that the graphene sheet continuously slides quite a little again (Figure 10(c)10: z=20.1 Å) until the next slip toward the neighboring AB stacking site occurs. As a result the trajectories of the graphene sheet exhibit the discrete zigzag paths connecting the nearest neighboring AB-stacking sites as shown in Figure 10(c)10. The period of the zigzag behavior of the Fz curve decreases from 3.7 Å to 2.5 Å as shown in Figure 10(a) as the peeling proceeds. The lattice spacing of the graphite surface, 2.5 Å, appears in the peeling force curve particularly for the stick-slip region.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the continuous and stick-slip process during the surface contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 5 indicated in (a). (c) The trajectories of the two carbon atoms on the free edge from 6 to 10 indicated in (a).
(a)(b)(c)


## 4.2.2. Line-Contact Region

During the line contact region between G and H in Figure9(a), z-Fz curve takes another atomic-scale zigzag structures as shown in Figure 11(a). One of the zigzag behaviors in the force curve (Figures 11(a)1→2→3) corresponds to the stick-slip sliding motion of the graphene sheet (Figures 11(b)1→2→3). Here the free edge of the graphene sheet slides with nearly the straight stick-slip motions. One of the carbon atoms on the free edge passes over the carbon-carbon bonds as shown in Figures 11(b)1 and 3.(a) Enlargement of part of thez-Fz curve (Figure 9(a)) corresponding to the nearly straight stick-slip region during the line contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 3 indicated in (a).
(a)(b)


## 5. Edge-Lifting Case of Zigzag-Edge Graphene

Recently it has been reported that the edge structure of the graphene sheet plays quite an important role in electronic, magnetic, and optical properties of graphene, which can be also expected to give influences on the mechanical properties such as, the peeling process. Therefore, in this section, the peeling process of the graphene sheet with zigzag edge is discussed. In the simulation, the model obtained by rotating Figure1(b) by 30∘ is used (Figure 1(c)), and the left zigzag edge is lifted to simulate the peeling process, while the right free edge is zigzag type. As a result, the nanoscale peeling process within the x-z plane and the global shape of the force curve (Figure 12) is similar to that of Figures 8 and 9, respectively. The qualitative tendency of the decrease of the period and amplitude of the force curve (Figure 12(b)I-IV) is similar to that for Figure 9(a)I–VII. However, the details of the atomic-scale mechanics of the zigzag edge are clearly different from those of the armchair edge as follows.Figure 12
The vertical force,Fz, acting on the lifting edge, plotted as a function of the lifting edge height z for the graphene with zigzag-type free edge.During the surface contact, the graphene sheet first takes zigzag (Figures13(b)1–6) and then straight stick-slip motions (Figures 13(c)7–11), passing over the nearest neighboring AB-stacking site along [101¯0] direction. It is noted to avoid AA-stacking registry, the graphene sheet takes zigzag slip toward the nearest neighboring AB-stacking site as shown in Figure 13(b)1→2, although it takes straight slip as shown in Figure 13(c)7→8. The minimum period of the force curve of 4.4 Å (Figure 13(a)IV) reflects the lattice period of the graphite surface along the [101¯0] direction, while 2.5 Å for the armchair-type edge (Figure 10(a)VII) reflects that along the [123¯0] direction. Thus the edge structure gives the marked effects on the atomic-scale dynamics depending on the lattice orientation of the surface.(a) Enlargement of part of thez-Fz curve (Figure 12(b)) corresponding to the zigzag and straight stick-slip process during the surface contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 6 indicated in (a). (c) The trajectories of the two carbon atoms on the free edge from 7 to 11 indicated in (a).
(a)(b)(c)During the line contact, the difference between the armchair- and zigzag-type edge is enhanced. Figure14(a) reflects the zigzag stick-slip motion of the graphene sheet (Figures 14(b)1–9) unlike nearly the straight stick-slip motion (Figures 10(b)1–3). Important point of the line-contact sliding is that each carbon atom on the free edge takes stick-slip motion between the nearest neighboring six-membered rings. When each atom is located on the hollow site of the six-membered ring, the graphene sheet does not deform along the y direction (Figures 14(b)1, 4-5, and 8-9). However, when each atom is located a little far from the hollow site or near the carbon bond, the graphene sheet bends toward the y direction to decrease the total interaction energy (Figures 14(b)2-3 and 6-7). Thus, in the case of the zigzag-type edge, collective motion of the single carbon “atom” on the free edge nearly dominates the graphene mechanics together with its deformation. On the other hand, for the armchair-type edge, collective motion of the single carbon “bond” is dominant.(a) Enlargement of part of thez-Fz curve (Figure 12(b)) corresponding to the zigzag stick-slip process during the line contact. (b) The trajectories of the two carbon atoms on the free edge from 1 to 9 indicated in (a).
(a)(b)


## 6. Discussions and Conclusions

In this work molecular mechanics study of the nanoscale peeling of the monolayer graphene sheet has been performed. The peeling force curve clearly exhibits the change of the graphene shape from the surface- to the line-contact.In Section3, the peeling of the monolayer graphene sheet with the armchair edge for lifting the center position is discussed. It is noted that the maximum lateral sliding force Fy⋍ 0.1 eV/Å is only about 3% of the absolute value of the maximum adhesion force |Fz|⋍3.1 eV/Å. This small sliding force Fy is derived from the superlubricity at the interface between the graphene sheet and the graphite surface [13] and atomic-scale wear [19]. There is a possibility that such anisotropy between the vertical force Fz and the lateral sliding force Fy can be applied to the adhesives, which can be strongly adhered to the substrate but can easily slide on it. Our AFM measurement exhibits that the maximum pull-off force is about several hundreds of nN, which is clearly much larger than the binding force, 3 eV/Å ⋍4.8 nN, assumed in our simulation.In Sections4 and 5, the peeling of the monolayer graphene sheet with the armchair- and zigzag-edge for lifting the edge position is discussed, respectively. The atomic-scale sliding motion of the monolayer graphene sheet during the peeling process is found. For the graphene sheet with armchair edge, the transition from the continuous to the stick-slip motion of the graphene sheet is found, which can be explained as follows. The peeling process induces the increase of the peeled area of the graphene sheet and the decrease of the surface contact area. Considering that the peeled area of the graphene sheet acts as an effective spring as shown in Figure 15, the increase of the peeled area makes the effective spring softer, and the decrease of the surface contact area decreases the energy barrier to slide the graphene sheet. Finally the peeling process induces the transition from the continuous to the stick-slip sliding motion of the graphene sheet, together with the decrease of the period and amplitude of the z-Fz curve. An important point is that the period of the peeling force curve for the armchair-edge graphene for the surface contact region corresponds to the lattice spacing of the graphite surface along [123¯0] direction, 2.5 Å. On the other hand, for the zigzag-edge graphene, the period becomes the lattice spacing along [101¯0] direction, 4.4 Å. This means the sliding length of the graphene sheet along x direction becomes nearly equal to the peeled length along z direction. The zigzag structures of the peeling force curve with the same period of about several Å have been also observed by our preliminary experiments using the multilayered graphene, which will be reported elsewhere [12]. Of course, if the number of the peeled graphene sheets is reduced, the direct comparison between the present simulation and the experiment will become possible.(a) Schematic illustration of the increase of the peeled area and the decrease of the surface contact area from C (z=6.1 Å) to D (z=18.3 Å) for the graphene sheet with armchair-type free edge. (b) -Fz and Fx plotted as a function of the edge height z, show qualitatively the same behavior to each other for the graphene sheet with armchair-type free edge.
(a)(b)As a result, the center-lifting case requires the largest attractive peeling force,-3.1 eV/Å, in order to peel the graphene sheet as shown in Figure 16(a). On the other hand, the edge-lifting case requires only -0.74 eV/Å, about 20% of that for the center-lifting case as shown in Figures 16(b) and 16(c). The edge structures give little influences on the basic features of the force curve. However, the sliding direction and the edge structure clearly give marked influences on the surface- and line-contact regions, respectively.Figure 16
Comparison of the vertical forces,Fz, among (a) center-lifting case of the armchair-edge graphene sheet, (b) edge-lifting case of the armchair-edge graphene sheet, and (c) edge-lifting case of the zigzag-edge graphene sheet.Another important point is that the behavior of the lateral force curve (Fx(z)) is qualitatively the same as that of the vertical force curve (Fz(z)) during the surface contact as shown in Figure 9(b). Therefore, it can be said that the peeling force curve, Fz(z), directly reflects the atomic-scale friction force, Fx(z), which decreases to 0.019 eV/Å ⋍30 pN for z=27.8 Å (Figure 9(b)). This ultralow friction force, Fx, is derived from the superlubricity at the interface between the graphene sheet and the graphite surface [19–21]. Furthermore, effect of the edge structure on the peeling process is clarified by comparison of the free edge between the armchair- and zigzag-types. As mentioned above, the atomic-scale structure of the force curve during the surface contact reflects the lattice spacing of the graphite surface. So the minimum period of the atomic-scale structure of the force curve can tell us the atomic-scale lattice orientation and structure of the free edge of graphene. Such information can be used for the control of the electronic properties of the graphene sheet adsorbed onto the substrate. Therefore, this paper indicates the possibility of the identification of the lattice orientation and the edge structure of the graphene sheet.In this paper, we discussed the importance of the dynamics of the free edge during the peeling process. On the other hand, we also found the importance of the shape of the graphene sheet. Additional simulated model and results are shown in Figure17. As shown in Figure 17(a), rectangular graphene sheet whose aspect ratio is different from that of the graphene sheet of Figure 1(b) is used. The basic shape of the vertical force curve, Figure 17(b), is similar to that of Figure 3. However, the armchair-type edge is peeled first for Figure 10(a), although the zigzag-type edge is peeled first for Figure 1(b) as discussed in Section III. This means that the shape of the graphene sheet plays an important role for deciding which edge is peeled first. Effect of the graphene shape on the peeling process will be discussed in detail somewhere in the near future.(a) The model of the red-colored monolayer graphene sheet physically adsorbed onto the blue-colored rigid graphite surface within thex-y plane. The green-colored six-membered ring at the center position is moved upward along the z (or [0001]) direction, by z=0.1 Å. Initial AB stacking registry of the red-colored graphene sheet with the blue-colored graphite surface is assumed. (b) The vertical force, Fz, acting on the center six-membered ring, plotted as a function of the lifting displacement z.
(a)(b)Lastly it should be noted that the peeling process discussed in this paper is closely related to the atomic-scale wear of the graphite and the graphene tip formation in the friction force microscopy [22]. When the tip is pushed onto the surface for less than the critical tip height, the outermost graphene layer is attached to the FFM tip, which results in the formation of the graphene tip. In that case, the graphene sheet takes the surface contact with the second layer graphene, and it takes the two-dimensional stick-slip motion. However, it is difficult to observe directly the stick-slip motion during the scan process, due to the very small gap between the FFM tip and the graphite surface. On the other hand, if the peeling process is used, it can be expected that the contact at the AFM tip/graphite interface has a wider space to be observed directly by ex Transmission Electron Microscopy (TEM). This paper indicates the possibility of a direct observation of the stick-slip motion of the graphene sheet, that is to say, the elementary process of the atomic-scale friction or superlubricity which occurs at the tip/graphite surface interface.



---

*Source: 742127-2010-09-27.xml*