Research Article | Open Access

Bohua Sun, "Formulation of 2D Graphene Deformation Based on Chiral-Tube Base Vectors", *Journal of Nanomaterials*, vol. 2010, Article ID 402591, 7 pages, 2010. https://doi.org/10.1155/2010/402591

# Formulation of 2D Graphene Deformation Based on Chiral-Tube Base Vectors

**Academic Editor:**Rakesh Joshi

#### Abstract

The intrinsic feature of graphene honeycomb lattice is defined by its chiral index (), which can be taken into account when using molecular dynamics. However, how to introduce the index into the continuum model of graphene is still an open problem. The present manuscript adopts the continuum shell model with single director to describe the mechanical behaviors of graphene. In order to consider the intrinsic features of the graphene honeycomb latticeโchiral index (), the chiral-tube vectors of graphene in real space have been used for construction of reference unit base vectors of the shell model; therefore, the formulations will contain the chiral index automatically, or in an explicit form in physical components. The results are quite useful for future studies of graphene mechanics.

#### 1. Introduction

A one-atom-thick layer of graphite called graphene is the โmotherโ of all graphitic forms. Graphene is a 2D building material for carbon materials of all other dimensionalities. It can be wrapped up to 0D buckyballs, rolled into 1D nanotubes, or stacked into 3D graphite [1โ3]. Figure 1 is a famous picture take on from [3].

The graphene has attracted attention because of its unusual two-dimensional structure and potential for applications. Owing to its exceptional mechanical properties and low mass density, graphene is an ideal material for use in nanoelectromechanical systems (NEMS), which are of great interest both for fundamental studies of mechanics at the nanoscale and for a variety of applications, including force, position, mass, and gas sensing [4].

Graphene research has developed at a truly relentless pace. Several papers appear every day; comprehensive reviews have been done by the graphene pioneers [2, 3]. Geometrically, the graphene can be considered as an โalmostโ 2D material due to its very small thickness at only one-atom. The 2D feature of the graphene makes it possible to directly use well-formulated 2D continuum shell model to simulate its deformation, vibration, and buckling [1โ3].

Although many features of the graphene can be described by relativistic quantum mechanics and molecular dynamics, the continuum modeling is still useful and crucial for rational device design and interpretation of experimental results. Because experiments at the nanoscale are extremely difficult and atomistic modelling remains prohibitively expensive for large-sized atomic system, it has been accepted that continuum models will continue to play an essential role in the study of carbon nanotube [5]. The validity of using continuum models for nanotube and graphene has been supported by both experimental results and molecular-dynamics simulation, all previous investigations on nanotube have indicated that the laws of continuum mechanics are still valid to some extent even in nanoscale. The success of continuum models to nanotube gives us confidence to predict that contunuum modelling is going to be also valid for graphene mechanics analysis [6, 7].

Within the frame of classic continuum mechanics, there is a well-formulated branch-theory of shell [8โ14]. It deals with a special 3D continuum shell; its thickness is very much smaller than its width and length. The deformation of the thinner 3D shell can be considered as a 2D middle surface with a small thickness, and all off-middle surface quantities can be presented in terms of middle surface. In this way, the 3D thinner shell problem is being converted into a 2D surface problem with a small thickness. Based on this understanding, it is to believe that the shell theory should be the best continuum mechanics model for graphene mechanics.

Regarding the mechanics modeling of the graphene, most of literature considered the graphene as a 2D pre-stretched membrane with zero thickness [15โ18], and some considered the graphene bending properties [19โ22]. Since monolayer graphene is very thin, to some extent the membrane theory could be very simple and quite good estimate model, but for tiny and delicated NEMS devices, any small thing must be taken into account as it might has a significant impact on graphene, for instance, grapheneโs thickness [3, 6].

Although grapheneโs thickness is very small, one layer graphene thickness is 0.335โnm [5], but it cannot be ignored as zero, because the thickness plays a key role to supply a bending stiffness. The graphene should be modelled as a very thinner 3D object by shell theory [8โ14].

The current shell theory [8โ14] can be directly applied to the deformation analysis of nanostructures, for instance, [7, 22]. Reference [22] considered the effects of surface tension on the elastic properties of nanostructures, and [7] discussed a promised nanomaterials nanotorus. However, how to introduce the chiral index () into continuum model of the graphene is still an open problem. The present manuscript adopts the continuum shell model with single director to describe the mechanical behaviors of graphene. In order to consider the intrinsic features of the graphene honeycomb lattice-chiral index (), the chiral-tube vectors of graphene in real space have been used for construction of reference unit base vectors of the shell model; therefore, the formulations will contain the chiral index automatically, or in explicit form in physical components. The results are quite useful for future studies of graphene mechanics.

The paper is organized as follows. Some elasticity concepts and kinematics of a deformable surface will be introduced in Section 2; the parameterization of single director shell and deformation mode will be demonstrated in Sections 3 and 4. The basis vectors will be established in Section 5. The deformation gradient for both with and without predeformation will be formulated in Sections 6 and 7. The explicit form of the chiral index () in the physical components will be shown in the Section 8.

#### 2. Some Elasticity Concepts

##### 2.1. Configuration, Deformation, and Parametrization

Let represent the stress-free or reference placement of the physical body with boundary . Points of *B* are indentified by . Let *B* be parametrized by the one to one point mapping function , such that where 2 is the parametric image of the body under the inverse mapping . Denote the deformed configuration of the body by with the boundary of given by . Points in S are identified by . The time-dependent configurations are parametrized by one to one point mapping function , such that where and is the prescribed placement of points on the boundary where the placement is given. The deformation mapping is given parametrically by . โNote that is the identity mapping.

##### 2.2. Tangent Basis Vectors and the Deformation Gradient

By differentiating the mapping function with respect to curvilinear coordinate (convected coordinate) function , the tangent basis vectors are obtained . The fact that the tangent basis vectors are linearly independent is expressed by the relationship . The reference or undeformed tangent basis vectors are obtained by particularizing the current basis vectors to , hence . The dual basis vectors and are defined by the orthogonality condition , . The deformation gradient is the mapping and is obtained by taking the material gradient of the deformation mapping as , where . In terms of the parametrization, **F** is given by . Using the above relations, the deformation gradient *F* can be given in terms of the current tangent basis vectors and the reference dual basis vectors by or in a coordinate format as , where , with relations , , and . Then, the determinant of the deformation gradient is given by , , . Since the current and reference tangent basis vectors are linearly independent, so have .

#### 3. Parameterization of Graphene Shell Model

The variable pair (**r, d**) is chosen normally as kinematic variable, in which **r** is the position vector and **d** is the director field of mid-surface. Since they are vectors, so the calculations corresponding to these quantities also operate in vector space.

The reasonable parameterization of director field plays a central role in development of shell models. The director field is decomposed into a positive magnitude parameter and the inextensible component or unit director field *t*, , where , , . The unit director field **t** is parameterized by two independent angles, that is, . The advantage of this kind of representation of **t** is that the variational calculation is very simple, for instance, .

#### 4. Single Director Kinematic Modelling of Graphene

The variable pair (**r, d**) is chosen normally as kinematic variable, in which **r** is the position vector and **d** is the director field of deforemed mid-surface and denoting **R** and **D** being the position vector and single director on the middle surface of undeformed configuration.

According to the long tradition of shell research [8โ14], the single director field in reference or undeformed configuration can be written as
The single director field in current configuration can be written as
Thus, the points in the body are indentified as the mid-surface mapping **r** or **R** plus the distance along the director **d** or **D**.

The deformation mode of shell is such that points intially defined along straight fibers, that is, indentified as the distance above or below the mid-surface along the line defined by the initial director **D**, remain a straight fibers. From the physical point of view, the solution to the shell problem, under the single director kinematic assumption, can be considered as the solution for the shell mid-surface plus the solution for the director which orients points along thickness fibers.

Let us now check, from geometric understanding, whether the single director kinematic assumption could capture the main characteristics of shell deformation. Finite membrane stretch is captured by the deformation of the mid-surface map **r**; finite bending strains are modeled by the spatial gradient of the director **d**; transverse shear strain is accounted for by measuring the relative rotation of the director **d** with respect to the normal to the mid-surface. The thickness change can be represented by the difference of magnitude of initial and current director.

#### 5. Establishment of Unit Base Vectors in the Reference State (Undeformed Graphene Sheet)

The way the graphene sheet is wrapped is represented by a pair of indices () called the chiral vector. The integers and denote the number of unit vectors along two directions in the honeycomb crystal lattice of graphene. If , the nanotubes are called โzigzagโ. If , the nanotubes are called โarmchairโ. Otherwise, they are called โchiralโ.

The orthogonal base vectors in the reference state of the grapheme can defined as follows: let **a **_{1} be a base vector **e **_{1}, normal base vector of the grapheme plane is , and then . The chiral vector can be expressed as follows in terms of orthogonal base vectors: . The tube axes vector **T** is orthogonal to the chiral vectors, that is, , then we have .The no-orthogonal base vectors can be represented by orthogonal base vectors .

Since **T** and are orthogonal each other, we can form another chiral-tube base vectors , where the normal vector or single director field in undeformed grapheme plane is .

Up to now, from unit vectors of grapheme in real space, we have formulated two base vectors and . In order to capture more information about the deformed grapheme sheet, we will try to formulate its shell model based on the chiral-tube base vectors; this work has not been found in the literature. For tensor calculation, we introduce to represent the chiral-tube base vectors as follows: Corresponding to the chiral-tube unit base vectors , we can define three axes (), for instance, the any point in the graphene plane can be expressed as follows: and any point out of the graphene plane is .

#### 6. Tangent Basis and Deformation Gradient of Graphene without Predeformation

If the grapheme sheet has no any predeformation, it will be like a plane. Since , we have the tangent basis vectors in the undeformed configuration as follows: After the deformation, the convected tangent basis vectors in current graphene shell configuration is defined by The dual basis vectors are also defined by the relationships and .

Then, tangent of deformation or deformation gradient **F** can be defined by
If we put in (5), then the basis vectors of the mid-surface in the current configuration can be obtained; that is,
and the corresponding dual basis vectors in the mid-surface are defined by .

#### 7. Deformation Gradient of Graphene with Pre-Deformation

If 2D graphene is being wrapped up to 0D nanodot, 1D nanotube, 3D stacked graphite, or been prestressed and then undergo deformation, it means that the graphene got pre-deformation.

After the pre-deformation, the convected tangent basis vectors in predeformation configuration are defined by The dual basis vectors are also defined by the relationships .

The predeformation gradient will be The after deformation gradient will be According to the superposition of deformation gradient, the total deformation gradient will be The associated Jacobian determinants are given by and from (8) and (9), and if we denote the director field gradient as then, we have where can be thought of as the mean and Gaussian curvature in the deformed configuration, and can be thought of as mean and Gaussian curvature in the pre-deformation configuration. Since no deformation in the undeformed graphene plane exists, so both its mean and Gaussian curvature are zero.

For instance, the cylindrical pre-deformation [21],, where *L* is the width of the grapheme sheet before rolling and *R* is the nanotube radius. The pre-deformation gradient in this case is

#### 8. Physical Components of Some Quantities in Terms of Chiral Index ()

The chiral index () is an intrinsic parameter of graphene. With our chiral-tube basis vectors, the influence of the chiral index () has been introduced into deformation gradient and can be shown explicitly in physical components.

Before going to discuss the physical components of tensor, let us here make a remark about the physical components of physical quantities. We should remember that the tensor (vector) components of a physical quantity which is referred to a particular curvilinear coordinate system may or may not have the same physical dimensions. In this case a great convenience arises because we would like to keep our freedom in choosing arbitrary curvilinear coordinates. For this reason, we must distinguish the tenos components from the โphysical componentsโ, which must have uniform physical dimensions, but they do not transfer conveniently under coordinate transformations.

First of all, let us consider the physical components of a vector . It is clearly that the base vectors are in general not unit vectors. In fact, their lengths are , i not in summation. Then are unit vectors. With the notation of the unit vectors, the vector can be rewriten as
where will have the same physical dimensions. It is seen that are the components of *v* resolved in the direction of unit vectors which are tangent to the coordinate lines, and that are the components of **u** resolved in the direction of unit vectors which are perpendicular to the coordinate planes. So the components are called the physical components of the vector **u**, and are physical base vectors. This means that that the physical components of a vector **u **can be defined as the components of **u** resolved in the directions of a set of unit vectors which are parallel either to the set of base vectors or to the set of reciprocal base vectors.

Similarly, for the reference or undeformed state, then . For tensor calculation, we define the chiral-tube unit physical base vectors as where , , and . It is indicated in (14) that chiral index has been included into the physical components of base vectors. Since all quantities and variables must be established on the base vectors, all of them will be represented by the chiral index. This is the beauty of this work. For โzigzagโ with , we have ; for โarchairโ with , we have

Using the above physical component in (14), for instance, the 2nd Piola-stress tensor of the shell model can be written in physical component form: Similarly, all other quantities can be explicitly expressed in terms of the chiral index ().

#### 9. Conclusion

This paper investigates the possibility of introducing the chiral index () into the deformation formulation of graphene continuum shell model. For this purpose, we established a chiral-tube basis vectors from graphene unit vectors in real space and formulated graphene deformation gradient for both with and without predeformation. The importance of the chiral index on the deformation and stress has been shown explicitly in the physical components.

#### Acknowledgments

This work was supported by South Africa National Research Foundation. This support is gratefully acknowledged. The author would like to express his thankfulness to the constructive comments and insightful suggestions given by three anonymous referees.

#### References

- A. K. Geim and K. S. Novoselov, โThe rise of graphene,โ
*Nature Materials*, vol. 6, no. 3, pp. 183โ191, 2007. View at: Publisher Site | Google Scholar - A. K. Geim, โGraphene: status and prospects,โ
*Science*, vol. 324, no. 5934, pp. 1530โ1534, 2009. View at: Publisher Site | Google Scholar - A. K. Geim and K. S. Novoselov, โThe rise of graphene,โ
*Nature Materials*, vol. 6, no. 3, pp. 183โ191, 2007. View at: Publisher Site | Google Scholar - B. Sun and X. M. Henry Huang, โMechanical nano-resonators at ultra-high frequency and their potential applications,โ
*South Africa Journal of Science*, vol. 104, pp. 169โ171, 2008. View at: Google Scholar - C. Q. Ru, โAxially compressed buckling of a doublewalled carbon nanotube embedded in an elastic medium,โ
*Journal of the Mechanics and Physics of Solids*, vol. 49, no. 6, pp. 1265โ1279, 2001. View at: Publisher Site | Google Scholar - Y. Wu, M. Huang, F. Wang et al., โDetermination of the young's modulus of structurally defined carbon nanotubes,โ
*Nano Letters*, vol. 8, no. 12, pp. 4158โ4161, 2008. View at: Publisher Site | Google Scholar - B. Sun, โDeformation, vibration, and buckling of nanotorus,โ
*Journal of Nanomaterials*. In press. View at: Google Scholar - W. Z. Chien, โThe instrisic theory of thin shells and plates, I,โ
*Quarterly of Applied Mathematics*, vol. 1, pp. 297โ327, 1943. View at: Google Scholar - W. Z. Chien, โThe instrisic theory of thin shells and plates, II,โ
*Quarterly of Applied Mathematics*, vol. 2, pp. 120โ135, 1944. View at: Google Scholar - W. T. Koiter, โOn the nonlinear theory of thin elastic shells,โ
*Proceedings of the Koninklijke Nederlandse Akademie van Wetenschappen. Series B*, vol. 69, pp. 1โ54, 1966. View at: Google Scholar - P. M. Naghdi,
*The Theory of Shells, Handbuch der Physik*, vol. 2, Springer, Berlin, Germany, 1972, edited by C. Truesdell. - J. L. Sanders, โNon-linear theories for thin shells,โ
*Quarterly of Applied Mathematics*, vol. 21, pp. 21โ36, 1963. View at: Google Scholar - J. C. Simo and D. D. Fox, โOn a stress resultant geometrically exact shell model. Part I: formulation and optimal parametrization,โ
*Computer Methods in Applied Mechanics and Engineering*, vol. 72, no. 3, pp. 267โ304, 1989. View at: Google Scholar - B. Sun and R. Liu, โSingle director shell model without complex geometric concepts,โ
*Advance in Mechanics*, vol. 35, no. 2, pp. 181โ194, 2003. View at: Google Scholar - C. Lee, X. Wei, J. W. Kysar, and J. Hone, โMeasurement of the elastic properties and intrinsic strength of monolayer graphene,โ
*Science*, vol. 321, no. 5887, pp. 385โ388, 2008. View at: Publisher Site | Google Scholar - J. Atalaya, A. Isacsson, and J. M. Kinaret, โContinuum elastic modeling of graphene resonators,โ
*Nano Letters*, vol. 8, no. 12, pp. 4196โ4200, 2008. View at: Publisher Site | Google Scholar - J. S. Bunch, S. S. Verbridge, J. S. Alden et al., โImpermeable atomic membranes from graphene sheets,โ
*Nano Letters*, vol. 8, no. 8, pp. 2458โ2462, 2008. View at: Publisher Site | Google Scholar - B. I. Yakobson, C. J. Brabec, and J. Bernholc, โNanomechanics of carbon tubes: instabilities beyond linear response,โ
*Physical Review Letters*, vol. 76, pp. 2511โ2514, 1996. View at: Google Scholar - M. Arroyo and T. Belytschko, โFinite crystal elasticity of carbon nanotubes based on the exponential Cauchy-Born rule,โ
*Physical Review B*, vol. 69, no. 11, Article ID 115415, 11 pages, 2004. View at: Google Scholar - Y. Huang, J. Wu, and K. C. Hwang, โThickness of graphene and single-wall carbon nanotubes,โ
*Physical Review B*, vol. 74, no. 24, Article ID 245413, 2006. View at: Publisher Site | Google Scholar - Q. Lu and R. Huang, โNonlinear mechanics of single-atom-layer graphene sheets,โ
*International Journal of Applied Mechanics*, vol. 1, no. 3, pp. 443โ467, 2009. View at: Google Scholar - Z.-Q. Wang, Y.-P. Zhao, and Z.-P. Huang, โThe effects of surface tension on the elastic properties of nano structures,โ
*International Journal of Engineering Science*, vol. 48, no. 2, pp. 140โ150, 2010. View at: Google Scholar

#### Copyright

Copyright ยฉ 2010 Bohua Sun. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.