SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

A correction for this article has been published in Nanoscale Research Letters 2014, 9:228

Open Access Nano Express

Nanoindentation experiments for single-layer rectangular graphene films: a molecular dynamics study

Weidong Wang1*, Shuai Li1, Jiaojiao Min1, Chenglong Yi1, Yongjie Zhan2 and Minglin Li3

Author Affiliations

1 School of Electrical and Mechanical Engineering, Xidian University, Xi'an 710071, China

2 Physics Department, Northwest University, Xi'an 710069, China

3 School of Mechanical Engineering and Automation, Fuzhou University, Fuzhou 350108, China

For all author emails, please log on.

Nanoscale Research Letters 2014, 9:41  doi:10.1186/1556-276X-9-41

The electronic version of this article is the complete one and can be found online at: http://www.nanoscalereslett.com/content/9/1/41


Received:9 October 2013
Accepted:14 January 2014
Published:22 January 2014

© 2014 Wang et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

A molecular dynamics study on nanoindentation experiments is carried out for some single-layer rectangular graphene films with four edges clamped. Typical load–displacement curves are obtained, and the effects of various factors including indenter radii, loading speeds, and aspect ratios of the graphene film on the simulation results are discussed. A formula describing the relationship between the load and indentation depth is obtained according to the molecular dynamics simulation results. Young’s modulus and the strength of the single-layer graphene film are measured as about 1.0 TPa and 200 GPa, respectively. It is found that the graphene film ruptured in the central point at a critical indentation depth. The deformation mechanisms and dislocation activities are discussed in detail during the loading-unloading-reloading process. It is observed from the simulation results that once the loading speed is larger than the critical loading speed, the maximum force exerted on the graphene film increases and the critical indentation depth decreases with the increase of the loading speed.

Keywords:
Molecular dynamics simulation; Nanoindentation; Rectangular graphene film; Young’s modulus; Strength

Background

The perfect graphene is a kind of carbonaceous material which consists of two-dimensional honeycomb lattice structures of single layer of carbon atoms. It is the basic unit to build other dimensional carbonaceous materials, such as zero-dimensional fullerenes, one-dimensional carbon nanotubes, and three-dimensional graphite [1,2]. Graphene sheets/ribbons/films have attracted the interest of the scientific community because of recent exciting experimental results [3-6]. Their growth, atomic makeup, electronics, doping, and intercalation have attracted many investigations [7-10]. A suspended graphene sheet [1,11] can be used in a variety of ways, such as for pressure sensors or gas detectors [12] or mechanical resonators [13].

It is still debatable whether a graphene sheet is truly a two-dimensional structure or if it should be regarded as a three-dimensional structure since it exhibits a natural tendency to ripple, as observed in recent experiments [2,14-16]. Carlsson addressed that an understanding of the coupling behaviors between bending and stretching of graphene sheets is necessary to fully explain the intrinsic ripples in a graphene sheet [15]. In addition to theoretical investigations, recent research has been carried out to measure the mechanical properties of suspended graphene sheets by utilizing an atomic force microscope (AFM) [17]. Through weak van der Waals forces, graphene sheets were suspended over silicon dioxide cavities where an AFM tip was probed to test its mechanical properties. Their Young’s modulus differs from that of bulk graphite. Poot and van der Zan [18] measured the nanomechanical properties of graphene sheets suspended over circular holes by using an AFM and suggested that graphene sheets can sustain very large bending and stretching prior to the occurrence of fracture, which indicates that the classical Kirchhoff plate theory used in the bending and vibration analysis of graphene sheets may not be suitable since deflection and stretching are considerable [19]. Some researchers thought that the large deflection plate theory of von Kármán may be a better candidate to model the graphene sheet, and they have characterized its bending and stretching through that theory [20,21]. Lee et al. measured Young’s modulus and the maximum stress of graphene by using an AFM in the nanoindentation experiment [22] and reported the effect of grain boundaries on the measurement of chemical vapor-deposited graphene [23]. Fang et al. [24] has studied the mechanical behavior of a rectangular graphene film under various indentation depths, velocities, and temperatures using molecular dynamics (MD) simulations. The physical models of the rectangular graphene film established by Fang et al. are doubly clamped using a bridge-type support and are loaded by a flat-bottomed diamond tip.

Although the above research has been carried out for mechanical properties of graphene, the variation of the values of elastic modulus has not yet reached a consensus. Especially for rectangular graphene films, the relationship between the load and the indentation depth is not clear. Furthermore, there are few papers available which describe the deformation mechanisms and dislocation activities of graphene film during the nanoindentation processes in detail. These investigations are concentrated on tension deformation [25-28] and shear deformation [29]. Almost all of the available literatures on dislocation activities in graphene focus on theoretical studies and numerical simulations, including density functional theory (DFT) [26], tight-binding molecular dynamics (TBMD) [30], ab initio total energy calculation [30], and quantum mechanical computations [31]. Researchers always artificially applied defects or dislocations and then studied their effects on the properties and activities in graphene. However, due to the bottleneck of experimental study at nanoscale, a very few experimental observations of dislocation activities are available at present. Warner et al. [32] also reported the observation of dislocation pairs through HRTEM experiments and gave five possible mechanisms that describe how these dislocation pairs could have formed, namely, during the CVD growth, electron beam sputtering of carbon dimers along a zigzag lattice direction, from surface adatom incorporation, from a monovacancy, and from a Stone-Wales defect. They then concluded that edge dislocations result in substantial deformation of the atomic structure of graphene, with bond compression or elongation of ±27%, plus shear strain and lattice rotations.

In this article, some MD simulations of nanoindentation experiments are performed on a set of single-layer rectangular graphene films with four clamped edges. The dislocation activities and the deformation mechanism are discussed, and a formula is introduced in order to describe the relationship of load and indentation depth and to measure the mechanical properties of graphene.

Methods

In order to carry out the nanoindentation experiments, one diamond sphere was introduced to simulate the indenter. Figure  1a shows the origin model for the nanoindentation experiment. Here, the upper ball is the indenter and constructed by diamond, which is considered as a rigid object so that the atomic configuration of the diamond indenter had no changes during MD simulations. The lower plane is a single-layer rectangular graphene film with different aspect ratios. For the inner atoms of the indenter and the graphene film, the energy function was described by adaptive intermolecular reactive empirical bond order (AIREBO) potential. Compared with the Tersoff-Brenner potential, the AIREBO potential not only introduces multi-body potential effects and the local atomic circumstance effect but also adds long-range interactions and a torsion term [33], and it has been found be more suitable than the Tersoff-type potential for accurately capturing Young’s modulus of graphene as well as bond breaking and reforming between carbon atoms [28]. In order to avoid the influence of nonphysical explanations with improper cutoff functions on the fracture process, the cutoff parameter of the AIREBO potential is set to be 2.0 Å. As for the interaction between the indenter and the graphene film, van der Waals forces were simulated based on the Lennard-Jones potential.

thumbnailFigure 1. Atomic configuration of the system model during the nanoindentation experiment. (a) The origin model, (b) the state during the loading process, and (c) at rupture state.

When performing MD simulations, we use the canonical (i.e., NVT) ensemble and control the temperatures at an ideal temperature of 0.01 K. In order to avoid the complex effects of the atomic thermal fluctuations, the temperature is regulated with the Nosé-Hoover method and the time step was set to 1 fs. During the simulation, one key step, named energy minimization and relaxation, should be carried out to make the system remain in the equilibrium state with lowest energy. Then, the indentation experiment was executed and the simulation results were output for further research.

Results and discussion

Loading and unloading properties

We take the case of the graphene film with an aspect ratio of 1.2 and the diamond indenter with a radius of 2 nm as an example to describe the indentation experiment in the following. The indenter was placed over the geometric center of the graphene film and forced to move in the direction perpendicular to the original graphene surface. Figure  1 gives the atomic configurations of the system model during the indentation experiment at a speed of 0.20 Å/ps. The atoms on the edge of the graphene film remained in a static state due to fixed boundary conditions. After enough loading time, the graphene film is eventually pierced through by the indenter, appearing some fractured graphene lattices. The load–displacement curves can be attained from the data of intender load (F) and indentation depth (d) calculated in MD simulations. The moment the load–displacement curve drops suddenly is considered to be a critical moment. In our simulations, the load suddenly decreased once the indentation depth exceeded 5.595 nm, defined as the critical indentation depth (dc), and the corresponding maximum load (Fmax) is 655.08 nN.

Figure  2 gives some detailed views on the graphene lattice fracture process starting from the critical moment. It is shown in Figure  2a that the carbon network was expanded largely, but there is no broken carbon-carbon (C-C) bond at the critical moment. Figure  2b represents the moment the bond-broken phenomenon emerged for the first time, with a pore appearing. The bond-broken process is irreversible and the load exerted on the graphene firstly declines. The first appearance of the pentagonal-heptagonal (5–7) and trilateral structures is shown in Figure  2c. Unlike with other Stone-Wales 5-7-7-5 (SW) defects mentioned in [34-36], in our indentation simulations, the loads are applied in the direction perpendicular to the surface of the graphene film and the 5–7 structure seems not conjugate but singular. In addition, it is found that the trilateral structure is an interim state in the evolution process from a pristine hexagonal structure to the 5–7 structure. A 5-3-6 structure including this trilateral structure and its adjacent structures would evolve into another 5–7 structure, the right one in Figure  2d, through bond breaking and bond reforming. Furthermore, a single-chain structure, shown in Figure  2e, can be observed during the fracture process, which can also be found in [26]. Afterwards, the single chain was broken and the indenter totally pierced through the graphene film.

thumbnailFigure 2. Evolution of graphene lattice fracture at different indentation depths. This group of figures shows the process from the state at which the indentation depth reaches the critical depth to the state the graphene film is totally ruptured with an indenter radius of 2 nm, loading speed of 0.20 Å/ps, and aspect ratio of 1.2. (a) At critical moment: indentation depth 55.95 Å, load 655.08 nN; (b) first broken bond emerged: indentation depth 55.97 Å, load 635.60 nN; (c) pentagonal-heptagonal (5–7) and trilateral structures emerged: indentation depth 55.99 Å, load 426.04 nN; (d) three 5–7 structures: indentation depth 56.01 Å, load 310.45 nN; (e) single-chain structure emerged: indentation depth 56.51 Å, load 112.03 nN; (f) fracture of the chain: indentation depth 56.61 Å, load 93.70 nN.

Generally speaking, elastic deformation which is reversible and plastic deformation which is irreversible are two typical kinds of deformation of an object or material in the view of engineering. In order to determine whether the deformation of the graphene film is elastic or plastic, a set of experiments of loading-unloading-reloading processes are conducted. As shown in Figure  3, during the continuous loading process of the indenter on the graphene film, it can be found that the graphene film mainly takes on two stages in sequence:

thumbnailFigure 3. Load–displacement curve of loading-unloading-reloading process with maximum indentation depth smaller than the critical indentation depth.

• Stage I. The unloading process is done before the indentation depth reaches the critical depth, dc. The graphene sheet almost can make a complete recovery, i.e., restore its initial structures, and the curves of reloading processes almost perfectly match the initial loading curve while the unloading curve shows very small deviations from the initial one, as shown in the inset of Figure  3. In general, the almost-perfect coincidence is due to the fact that the carbon covalent bonds and the graphene lattice structure are not destroyed. It can be concluded that there is no plastic deformation in this stage, i.e., the graphene undergoes elastic deformation.

• Stage II, i.e., the yellow region in Figure  3. In Figure  4a,b,c, it is after the indentation depth exceeded dc that the unloading process begins. The graphene sheet cannot make a complete recovery, and there exited broken covalent bonds after the unloading process. In the reloading process, the maximum force exerting on graphene is much smaller than that in Figure  3, which denotes the fracture of graphene lattices. Figure  4b describes the state where the unloading process begins, and Figure  4c describes the state where the unloading process ends. After the loading process, there exited broken bonds and fractured lattices in the middle of the graphene film and these defective structures did not recover during the unloading process. Therefore, the deformation of the graphene described in this figure can be considered as a plastic type.

thumbnailFigure 4. Loading-unloading-reloading process with the maximum indentation depth smaller than the critical indentation depth. (a) Load–displacement curve, (b) local atom configuration when the loading process is finished, and (c) local atom configuration when the unloading process is finished.

Young’s modulus and strength of the graphene film

According to the available correlation for the indentation experiments of a circular single-layer graphene film in [18,22,37], one new formula is constructed to describe the relationship between indention depth and load,

F = F σ d + F E d (1)

where d is the indentation depth and F denotes the concentrated force gotten by the graphene film. In Equation 1, the load F consists of two parts: the first part, Fσ(d), represents the term due to the axial tension of the two-dimensional (2-D) film,

F σ d = σ 0 2 D πR β 1 / 2 r R equ 3 / 4 d R equ (2)

where σ02D is the pre-tension of the single-layer graphene film, r is the indenter radius, β denotes the aspect ratio and is equal to L/b, and Requ represents the equivalent radius of the rectangular graphene sheet, (Lb/π)1/2. The second one, FE(d), represents the large deformation term,

F E d = E 2 D q 3 R equ β 1 / 2 r R equ 1 / 4 d R equ 3 (3)

where E2D is the 2-D elastic modulus, i.e., Young’s modulus, of the single layer graphene film. The strain energy density of graphene, as a standard 2-D material, can be represented by the energy of per unit area. Then, the corresponding pre-tension and elastic modulus can be expressed as σ02D and E2D, respectively, with the unit N/m. The common pre-tension and elastic modulus of a 3-D bulk material can be obtained through these 2-D values divided respectively by the effective thickness which is always treated as the layer spacing of the graphite crystal, i.e., 3.35 Å. q is an nondimensional value, q = 1/(1.05 - 0.15ν - 0.16ν2) = 0.9795, where ν denotes Poisson’s ratio, ν = 0.165 [3,18,21]. It is reported that when r/R > 0.1, the indenter radius has a significant influence on the load–displacement properties [38,39]. In our simulations, r/R > 0.1; thus, Equations 2 and 3 are corrected by a factor of (r/R)3/4 and (r/R)1/4, respectively.

According to the expression of maximum stress in the plane of a circular elastic film in [40], the maximum stress for a rectangular elastic film can be expressed by the following formula:

σ m 2 D = F max E 2 D 8 πr 1 / 2 (4)

Figure  5 presents the curves of indentation depth versus load for the nanoindentation experiment. The red solid curve is from the MD simulation results. According to Equation 1, nonlinear least squares method was used to fit the simulation results, and then the black curve in Figure  5 can be obtained. It is noted that when the indentation depth is about 5.597 nm, the load received by the graphene film suddenly drops from approximately 655.08 to approximately 522.172 nN. Corresponding to Figure  2b,c, the lengths of C-C bonds under the indenter quickly become larger than before, which indicates that the bonds were broken.

thumbnailFigure 5. Curves of indentation depth versus load for the nanoindentation experiment.

Table  1 gives the mechanical properties calculated from the MD simulation results. Young’s modulus and the maximum stress of the graphene are obtained as 1.0539 TPa and 205.1328 GPa, respectively. Young’s modulus obtained in this paper is in good agreement with those obtained by both experimental and numerical methods. Kudin et al. has predicted a Young’s modulus of 1.02 TPa using ab initio methods [41]. Lee et al. obtained a Young’s modulus of 1 ± 0.1 TPa by nanoindentation in an AFM of freestanding monolayer circular graphene membranes [22]. Neek-Amal and Peeters studied the nanoindentation of a bilayer graphene using molecular dynamics simulations and estimated a Young’s modulus of 0.8 TPa [42]. In addition, the maximum stress ranges from 130 to 240 GPa by means of both experiments and numerical simulations reported in other literatures [21,22,43,44]. The maximum stress obtained in this paper can also be included in the above range, which verified our simulation results. The changing trend of 2-D pre-tension demonstrates that the pre-tension of the rectangular graphene film is positively correlated with the loading speed of the indenter. The indenter size also affects the pre-tension, which, to some extent, explains why the correction factors were introduced in Equations 2 and 3.

Table 1. Mechanical properties of the single-layer graphene film from nanoindentation experiments

Other parameters’ influences on nanoindentation experiments

For further study of nanoindentation properties, a series of simulations have been carried out with different loading speeds, indenter radii, and aspect ratios of graphene film. It is indicated that the speed of 0.20 Å/ps can be regarded as the critical loading speed, as shown in Figure  6a. When the loading speed is higher than the critical value, with the increase of speed, the maximum load increases rapidly; simultaneously, the critical indentation depth decreases rapidly. However, when the loading speed is lower than the critical value, the changes of Fmax and dc are not that obvious. When the loading speed decreases from 1.00 to 0.50 Å/ps, dropping by 50%, the value of dc increases by 33.35%, and the value of Fmax decreases by 8.43% correspondingly. Nevertheless, when the loading speed decreases from 0.20 to 0.10 Å/ps, dropping by 50%, the changes of Fmax and dc are only 1.68% and 0.21%, respectively. The results may be attributed to the fact that the higher the loading speed of the indenter, the less time it takes to go through the graphene sheet, resulting in a higher load and lower indentation depth than those at a lower loading speed, in which situation the load process is much slower. Secondarily, the spherical indenter’s influences on results are observed by changing the indenter radius. The simulations of various indenter radii (1, 2, and 3 nm) are carried out at the speed of 0.20 Å/ps. The results of the load–displacement curve are shown in Figure  6b. The stress is more uniform in the middle of the graphene, so the Fmax increases obviously and the critical indentation depth also becomes greater with the increase of the indenter radius. Finally, after changing the aspect ratio of the graphene film to 1.1 and 1.5, Young’s modulus and the maximum stress of the graphene are obtained using the methods mentioned above. It can be deduced from Figure  7 that Young’s modulus and the maximum stress are the inherent properties of graphene and irrelevant to its size, which also verifies the formula obtained above.

thumbnailFigure 6. Comparison of load versus indentation depth for different parameters. (a) The indenter is loaded at different loading speeds between 0.10 and 2 Å/ps. (b) The indenter is loaded with different indenter radii of 1, 2, and 3 nm.

thumbnailFigure 7. Comparison of Young’s modulus and maximum stress versus the aspect ratio of the graphene film.

Conclusions

Some MD simulations of nanoindentation experiments on single-layer rectangular graphene sheets have been carried out in order to obtain the mechanical properties of graphene. A correlation between the load and the indentation depth is constructed, and Young’s modulus and the strength of graphene are obtained in the end. The simulation results show that the unloaded graphene film could make a complete recovery if the maximum indentation depth is less than the critical indentation depth, and the graphene film undergoes elastic deformation during the whole loading-unloading-reloading process. However, if the maximum indentation depth is larger than the critical indentation depth, the graphene sheet could not restore its original structures after unloading and the graphene deforms plastically. Based on the simulations of nanoindentation at different loading speeds and indenter radii, it can be observed that the maximum load increased and the critical indentation depth decreased with the increase of loading speed. In addition, the indenter radius has a remarkable influence on the force-displacement curve. As the indenter radius increases, the critical load and the critical indentation depth also increase.

Abbreviations

2-D: two-dimensional; AFM: atomic force microscope; AIREBO: adaptive intermolecular reactive empirical bond order; GPa: gigapascals; MD: molecular dynamics; NVT: constant number of particles, volume, and temperature; SW: Stone-Wales 5-7-7-5; TPa: terapascals.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The analysis of the simulation results was mainly carried out by WDW. The simulation processes were mainly conducted by SL, JJM, and CLY. Some fairly helpful proposals about the construction of models were made by YJZ and MLL. All authors read and approved the final manuscript.

Acknowledgements

We acknowledge the financial support provided by the Fundamental Research Funds for the Natural Science Basic Research Plan in Shaanxi Province of China (grant no. 2013JM7017), subsequently by the National Natural Science Foundations of China (grant no. 51205302 and no. 50903017) and the Central Universities in Xidian University (grant no. K5051304006). We also would like to thank all the reviewers for their comments and kind suggestions to our manuscript and all the editors for their careful corrections on the final version of the article.

References

  1. Geim A, Novoselov K: The rise of graphene.

    Nat Mater 2007, 6:183-191. PubMed Abstract | Publisher Full Text OpenURL

  2. Wang W, Hao Y, Yi C, Ji X, Niu X: Relaxation properties of graphene nanoribbons at different ambient temperatures: a molecular dynamics.

    Acta Phys Sin 2012, 61:200207. OpenURL

  3. Novoselov K, Geim A, Morozov S, Jiang D, Zhang Y, Dubonos S, Grigorieva I, Firsov A: Electric field effect in atomically thin carbon films.

    Science 2004, 306:666-669. PubMed Abstract | Publisher Full Text OpenURL

  4. Novoselov K, Jiang Z, Zhang Y, Morozov S, Stormer H, Zeitler U, Maan J, Boebinger G, Kim P, Geim A: Room-temperature quantum hall effect in graphene.

    Science 2007, 315:1397. OpenURL

  5. Barzola-Quiquia J, Esquinazi P, Rothermel M, Spemann D, Butz T, Garcia N: Experimental evidence for two-dimensional magnetic order in proton bombarded graphite.

    Phys Rev B 2007, 76:1403. OpenURL

  6. Peter W, Jan-Ingo F, Eli A: Epitaxial graphene on ruthenium.

    Nat Mater 2008, 7:406-411. PubMed Abstract | Publisher Full Text OpenURL

  7. Chiu Y, Lai Y, Ho J, Chuu D, Lin M: Electronic structure of a two-dimensional graphene monolayer in a spatially modulated magnetic field: peierls tight-binding model.

    Phys Rev B 2008, 77:045407. OpenURL

  8. Dutta S, Lakshmi S, Pati S: Electron–electron interactions on the edge states of graphene: a many-body configuration interaction study.

    Phys Rev B 2008, 77:073412. OpenURL

  9. Lai Y, Ho J, Chang C, Lin M: Magnetoelectronic properties of bilayer Bernal graphene.

    Phys Rev B 2008, 77:085426. OpenURL

  10. Lherbier A, Biel B, Niquet Y, Roche S: Transport length scales in disordered graphene-based materials: strong localization regimes and dimensionality effects.

    Phys Rev Lett 2008, 100:036803. PubMed Abstract | Publisher Full Text OpenURL

  11. Meyer J, Geim A, Katsnelson M, Novoselov K, Booth T, Roth S: The structure of suspended graphene sheets.

    Nature 2007, 446:60-63. PubMed Abstract | Publisher Full Text OpenURL

  12. Schedin F, Geim A, Morozov S, Hill E, Blake P, Katsnelson M, Novoselov K: Detection of individual gas molecules adsorbed on graphene.

    Nat Mater 2007, 6:652-655. PubMed Abstract | Publisher Full Text OpenURL

  13. Bunch J, van der Zande A, Verbridge S, Frank I, Tanenbaum D, Parpia J, Craighead H, McEuen P: Electromechanical resonators from graphene sheets.

    Science 2007, 315:490-493. PubMed Abstract | Publisher Full Text OpenURL

  14. Fasolino A, Los J, Katsnelson M: Intrinsic ripples in graphene.

    Nat Mater 2007, 6:858-861. PubMed Abstract | Publisher Full Text OpenURL

  15. Carlsson J: Graphene: buckle or break.

    Nat Mater 2007, 6:801-802. OpenURL

  16. Zhou J, Huang R: Internal lattice relaxation of single-layer graphene under in-plane deformation.

    J Mech Phys Solids 2008, 56:1609-1623. Publisher Full Text OpenURL

  17. Frank I, Tanenbaum D, van der Zande A, McEuen P: Mechanical properties of suspended graphene sheets.

    J Vac Sci Technol B 2007, 25:2558-2561. Publisher Full Text OpenURL

  18. Poot M, van der Zan H: Nanomechanical properties of few-layer graphene membranes.

    Appl Phys Lett 2008, 92:063111. Publisher Full Text OpenURL

  19. Duan W, Wang C: Nonlinear bending and stretching of a circular graphene sheet under a central point load.

    Nanotechnology 2009, 20:077702. OpenURL

  20. Landau L, Lifshits E: Theory of Elasticity. New York: Pergamon; 1970. OpenURL

  21. Yang X, He P, Wu A, Zheng B: Molecular dynamics simulation of nanoindentation for graphene.

    Scientia Sinica: Phys, Mech, Astron 2010, 40:353-360. OpenURL

  22. Lee C, Wei X, Kysar J, Hone J: Measurement of the elastic properties and intrinsic strength of monolayer graphene.

    Science 2008, 321:385-388. PubMed Abstract | Publisher Full Text OpenURL

  23. Lee G, Cooper R, An S, Lee S, Zande A, Petrone N, Hammerberg A, Lee C, Crawford B, Oliver W, Kysar J, Hone J: High-strength chemical-vapor–deposited graphene and grain boundaries.

    Science 2013, 340:1073-1076. PubMed Abstract | Publisher Full Text OpenURL

  24. Fang T, Wang T, Yang J, Hsiao Y: Mechanical characterization of nanoindented graphene via molecular dynamics simulations.

    Nanoscale Res Lett 2011, 6:481. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Kiselev S, Zhirov E: Molecular dynamic simulation of deformation and fracture of graphene under uniaxial tension.

    Phys Mesomech 2013, 16:125-132. Publisher Full Text OpenURL

  26. Topsakal M, Ciraci S: Elastic and plastic deformation of graphene, silicene, and boron nitride honeycomb nanoribbons under uniaxial tension: a first-principles density-functional theory study.

    Phys Rev B 2010, 81:024107. OpenURL

  27. Xu Z: Graphene nano-ribbons under tension.

    J Comput Theor Nanos 2009, 6:1-3. Publisher Full Text OpenURL

  28. Zhao H, Min K, Aluru N: Size and chirality dependent elastic properties of graphene nanoribbons under uniaxial tension.

    Nano Lett 2009, 9:3012. PubMed Abstract | Publisher Full Text OpenURL

  29. Carpio A, Bonilla L: Periodized discrete elasticity models for defects in graphene.

    Phys Rev B 2008, 78:085406. OpenURL

  30. Lee G, Yoon E, Hwang N, Wang C, Ho K: Formation and development of dislocation in graphene.

    Appl Phys Lett 2013, 102:021603. Publisher Full Text OpenURL

  31. Dumitrica T, Hua M, Yakobson B: Symmetry-, time-, and temperature-dependent strength of carbon nanotubes.

    Proc Natl Acad Sci U S A 2006, 103:6105-6109. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Warner J, Margine E, Mukai M, Robertson A, Giustino F, Kirkland A: Dislocation-driven deformations in graphene.

    Science 2012, 337:209. PubMed Abstract | Publisher Full Text OpenURL

  33. Wang W, Yi C, Ji X, Niu X: Molecular dynamics study on relaxation characteristics of graphene nanoribbons at room temperature.

    Nanosci Nanotechnol Lett 2012, 4:1188-1193. Publisher Full Text OpenURL

  34. Zhao H, Aluru N: Temperature and strain-rate dependent fracture strength of graphene.

    J Appl Phys 2010, 108:064321. Publisher Full Text OpenURL

  35. Xu L, Wei N, Zheng Y: Mechanical properties of highly defective graphene: from brittle rupture to ductile fracture.

    Nanotechnology 2013, 24:505703. PubMed Abstract | Publisher Full Text OpenURL

  36. Xiao J, Staniszewski J, Gillespie J Jr: Fracture and progressive failure of defective graphene sheets and carbon nanotubes.

    Compos Struct 2009, 88:602-609. Publisher Full Text OpenURL

  37. Komaragiri U, Begley M, Simmonds J: The mechanical response of freestanding circular elastic films on the point and pressure loads.

    J Apple Mech 2005, 72:203-212. Publisher Full Text OpenURL

  38. Begley M, Mackin T: Spherical indentation of freestanding circular thin films in the membrane regime.

    J Mech Phys Solid 2004, 52:2005-2032. Publisher Full Text OpenURL

  39. Scott O, Begley M, Komaragiri U, Mackin T: Indentation of freestanding circular elastomer films using spherical indenters.

    Acta Mater 2004, 52:4877-4885. Publisher Full Text OpenURL

  40. Bhatia N, Nachbar W: Finite indentation of elastic-perfectly plastic membranes by a spherical indenter.

    AIAA J 1968, 6:1050-1057. Publisher Full Text OpenURL

  41. Kudin K, Scuseria G, Yakobson B: C2F, BN, and C nanoshell elasticity from ab initio computations.

    Phys Rev B 2001, 64:235406. OpenURL

  42. Neek-Amal M, Peeters FM: Nanoindentation of a circular sheet of bilayer graphene.

    Phys Rev B 2010, 81:235421. OpenURL

  43. Wu J, Hwang C, Huang Y: An atomistic based finite deformation shell theory for single-wall carbon nanotubes.

    J Mech Phys Solid 2008, 56:279-292. Publisher Full Text OpenURL

  44. Lu Q, Arroyo M, Huang R: Elastic bending modulus of monolayer graphene.

    J Phys D Appl Phys 2009, 42:245413. OpenURL