# On the Vibration of Single-Walled Carbon Nanocones: Molecular Mechanics Approach versus Molecular Dynamics Simulations

Authors: R. Ansari; A. Momen; S. Rouhi; S. Ajori
Journal: Shock and Vibration
 (2014)

---

## Abstract

The vibrational behavior of single-walled carbon nanocones is studied using molecular structural method and molecular dynamics simulations. In molecular structural approach, point mass and beam elements are employed to model the carbon atoms and the connecting covalent bonds, respectively. Single-walled carbon nanocones with different apex angles are considered. Besides, the vibrational behavior of nanocones under various types of boundary conditions is studied. Predicted natural frequencies are compared with the existing results in the literature and also with the ones obtained by molecular dynamics simulations. It is found that decreasing apex angle and the length of carbon nanocone results in an increase in the natural frequency. Comparing the vibrational behavior of single-walled carbon nanocones under different boundary conditions shows that the effect of end condition on the natural frequency is more prominent for nanocones with smaller apex angles.


---

## Body

## 1. Introduction

Recently, different nanostructures such as carbon nanotubes (CNTs) [1], fullerenes [2], carbon nanorings [3], and carbon nanocones (CNCs) [4, 5] have attracted a great deal of researches. Due to diverse potential applications of single-walled carbon nanocones (SWCNCs) in different areas such as cold electron and field emitter [6], adsorbent [7], and mechanical sensors [8, 9], a comprehensive understanding of their mechanical, physical, and electronic properties is necessary.Researches on CNCs have been started at the same time as CNTs. However, the CNC and its properties attracted less attention. The cones with apex angles of 19.2°, 38.9°, 60°, 84.6°, and 112.9° were experimentally synthesized by Krishnan et al. [11]. Subsequently, Naess et al. [12] studied the morphologies of the CNCs with different apex angles by using transmission electron microscopy (TEM), synchrotron X-ray, and electron diffraction. Employing a quantum chemical optimization, the accurate geometries of CNCs were determined by Compernolle et al. [13]. It was shown that both armchair and zigzag edge atoms are present simultaneously in the same cone.Different theoretical approaches have been used to study the properties of nanostructures which can be divided as atomistic and continuum approaches. Molecular dynamics (MD) simulations [14, 15] and density functional theory (DFT) [16] are the most significant approaches of the former approaches, while the latter one can be subdivided to Bernoulli-Euler/Timoshenko beam models [17–19], the shell models [10, 20, 21], the molecular structural models [22–25], and meshless approaches [26–28]. Due to considerable time consumption of the atomistic approaches, continuum methods have attracted more attention. Due to similarity between the molecular structure of nanostructures and a space frame structure, the molecular structural methods have been used to study the mechanical behavior of nanostructures. In these methods, to obtain the physical properties of the model, the energy terms in molecular and structural models are equated.As in meshless methods, the problem domain is modeled by using a set of scattered nodes; the common difficulties of mesh usage including inaccurate and unstable solutions caused by distortion of meshes [26] and high costs of remeshing are omitted [27]. So, they can be usefully applied to study the nanostructures. Recently, meshless methods have been used extensively in the investigation of functionally graded panels reinforced by nanostructures [26–28].Some theoretical studies on mechanical behavior of CNCs can be found in the literature. MD simulation was used by Wei et al. [29] to study the mechanical behavior of CNCs. Liew et al. [30] analyzed the axial buckling of CNCs. Using the beam model and MD simulations, Hu et al. [8] investigated the free transverse vibration of SWCNCs. Using a nonlocal continuum shell model, Firouz-Abadi et al. [31] studied free vibration characteristics of nanocones. They used Galerkin technique to study the influences of geometrical parameters of the nanocone on these frequencies. MD simulations were also employed by Jordan and Crespi [32] to study the mechanical chiral inversion of closed CNCs. An atomically accurate chiral inversion was also obtained in their work, without breaking the chemical bonds. Utilizing the higher order continuum theory, Yan et al. [33] studied the structural parameters and elastic properties of SWCNCs. They employed all five types of CNCs to test the influence of the conical angle on the mechanical properties. A modal analysis of SWCNTs and SWCNCs was performed by J. H. Lee and B. S. Lee [34] employing finite element method (FEM). In their work, natural frequencies and corresponding mode shapes of these nanostructures were calculated. Using FEM, Fakhrabadi et al. [35, 36] predicted vibrational properties and the elastic and buckling behavior of CNCs with different dimensions. They obtained the natural frequencies and corresponding mode shapes as well as the elastic modulus and compressive forces of the axial buckling for CNCs with the various dimensions and boundary conditions.In this paper, a finite element (FE) model is applied to investigate the vibration behavior of SWCNCs. The effect of boundary condition, apex angle, and length on the fundamental natural frequency of SWCNCs are studied. Besides, MD simulations are used to obtain the natural frequency of SWCNCs. The results of FE model are verified by MD simulations.


## 2. Atomic Structure of SWCNCs

To form a SWCNC, a sector of circular graphene sheet is removed; then the graphene sheet is rolled around the apex and the open edges are joined. Disclination angledθ, which is defined as the angle of the sector removed from a plain sheet, is utilized to characterize the nanocones. The value of dθ can only have certain values including 60°, 120°, 180°, 240°, and 300°. Figure 1 displays the sectors of graphene sheet with disclination angles of 60°, 120°, 180°, and 240°. The apex angle of cones is obtained as 2sin-1⁡(1-dθ/360).Figure 1
Cone sheets with disclination angles of 60°, 120°, 160°, and 240°.Transformation from a graphene sheet with disclination angle 180° in a polar coordinate system to a three-dimensional cone in Cartesian coordinate system is displayed in Figure2. The polar coordinates of point P in the graphene sector are indicated as (l,α); L and φ represent the parameters of the graphene sector in which the lengths of OA- and OB- are equal. Sector OAB is folded at apex O, and OA- is overlapped by OB- [20].Symbols and variables used in the transformation from a graphene sheet for a cone 180°.
(a)
(b)Considering the pointP(X,Y,Z), the following equations are obtained [34]:
(1)X=rcos⁡⁡β,Y=rsin⁡β,Z=-l2-r2,
in which β and r are the unknown variables. The transformation from the graphene sheet to the three-dimensional cone leads to a relation between the corresponding angles as follows:
(2)β2π=αφ.
Rearranging (2) leads to the dihedral angle β between O′2O′C′ and O′2O′A′ in the following form:
(3)β=α(2πφ).
By involving lengths l and L, the relation between radii r and R can also be specified as follows:
(4)rR=lL.
The base radius of the cone is obtained as follows:
(5)R=Lφ2π.
Substitution of (5) into (4) gives the radius r as follows:
(6)r=l2π.
Substituting (3) and (6) into (1) yields the coordinates of point P in the three-dimensional cone in the following form:
(7)X=(lφ2π)cos⁡⁡[α(2πφ)],Y=(lφ2π)sin⁡[α(2πφ)],Z=-l1-(φ2π),
where X, Y, and Z indicate the nanocone atom coordinates, and angles are given in radians.


## 3. Molecular Structural Modeling

The discrete structures of CNCs, as presented in Figure3, with different lengths and apex angles are utilized in the present finite element modeling. A three-dimensional elastic beam element is employed to model the covalent bonds between the carbon atoms. Carbon atoms are considered as point mass elements. The total potential energy is the summation of energies that emerged from the different interatomic interactions and can be defined as follows [25]:
(8)E=∑Ur+∑Uθ+∑Uϕ+∑Uω+∑UvdW,
in which Ur, Uθ, Uϕ, Uω, and UvdW are energies corresponding to the bond stretching, angle variation, dihedral angle torsion, out-of-plane torsion, and van der Waals (vdW) interactions, respectively. The physical interpretation of these energy terms is given in Figure 4. Due to the small contribution of the vdW interaction in the SWCNC structure, it can be ignored in comparison with the other energies. Discarding the final term of (8), the other terms can be obtained in the forms of
(9)Ur=12kr(Δr)2,d2UrdΔr2=kr,Uθ=12kθ(Δθ)2,d2UθdΔθ2=kθ,Uτ=Uϕ+Uω=12kφ(Δφ)2,d2UτdΔφ2=kφ,
where kr, kθ, and kφ represent the bond stretching force constant, bond angle variation force constant, and out-of-plane angle variation resistance, respectively. In (9), the terms Δr, Δθ, and Δφ indicate the increase in bond length, bond angle change, and out-of-plane change, respectively. Equalization of the well-known structural energies with molecular energies in (9) leads to the following formulations [25]:
(10)EAL=kr,EIL=kθ,GJL=kφ,
where E and G are the elastic modulus and shear modulus, and A, L, I, and J are cross section area, length, moment of inertia, and polar moment of inertia of the beam element, respectively. After simplifications, the mechanical constants of the three-dimensional beam element with circular cross section having diameter d are obtained as follows [25]:
(11)d=4kθkr,E=kr2L4πkθ,G=kr2kφL8πkθ.
The values of the cross-sectional area, the area moment of inertia (Iyy=Izz), and the torsional moment of inertia (Ixx) are considered to be 1.68794 Å2, 0.453456 Å4, and 0.22682 Å4, respectively, and the covalent bond distance of the carbon atoms is 0.142 nm [25]. Substitution of these values into (11) results in the following:
(12)d=1.466Å,E=5.488×10-8N/Å2,G=8.701×10-9N/Å2.Cones with apex angles of (a) 38.9°, (b) 60°, (c) 86.6°, and (d) 123.6°, with cone heights of 20 Å.
(a)
(b)
(c)
(d)Figure 4
Various bonds between carbon atoms.


## 4. Application of FEM in the Vibrational Analysis of the CNC

In this section, the finite element model which is used to analyze the vibrational behavior of the CNCs is clarified. The FEM is performed using the commercial code ANSYS. The CNCs with different boundary conditions are modeled to obtain the natural frequencies and the corresponding mode shapes. As it was mentioned in Section3, the covalent bonds between two carbon atoms are modeled using the beam element and each atom is considered as point mass element. The free vibration of a nanostructure can be formulated in the form of
(13)[M]{q¨}+[K]{q}=0,
where [M], [K], {q}, and {q¨} indicate the global mass matrix, global stiffness matrix, nodal displacement vector, and nodal acceleration vector, respectively. The aforementioned matrices should be assembled for the whole of a structure. For a linear system, the nodal displacement vector is considered as follows:
(14){q}={φ}icos⁡⁡ωit,
in which ωi is the ith natural frequency, {φ}i is the vector that represents the mode shapes associated with ωi, and t denotes the time variable. Substituting (14) into (13) leads to the following:
(15)([K]-ωi2[M]){φ}i={0},
in which {φ}i is ith mode shape. To obtain the natural frequencies, the determinant of the term in the parentheses in (15) must be equal to zero. By substitution of the natural frequencies into (15), the corresponding mode shapes are determined. The above finite element formulation and simulation were carried out using shifted Block Lanczos technique in ANSYS which is applicable for symmetric eigenvalue problems [37].


## 5. Molecular Dynamics

In order to assess the results generated by FEM for vibrational behavior of SWCNC, MD simulation is utilized herein. In this regard, the Tersoff-Brenner [38, 39] potential energy function is used to calculate energetics of carbon bonds within the canonical ensemble known as NVT at room temperature (300 K). This potential function is given by
(16)V(rij)=(12)fc(rij)[VR(rij)+bijVA(rij)],
where rij expresses the distance between atoms i and j, VR(rij) and VA(rij) stand for the repulsive and attractive pairwise interaction terms, respectively, and bij is the strength of VA(rij).To integrate the equations of motion through time, Velocity-Verlet algorithm [40] is employed together with Nose-Hoover algorithm [41]. Moreover, to conserve energy and reducing temperature fluctuation and prevent system instability, the time step of 1 fs is chosen. Further, to compute the natural frequency of SWCNC, the system initially was relaxed for 3 ps to reach its optimum state of energy; then after applying initial proper deflection similar to first mode shape [42], it is allowed to vibrate freely. Finally, using the time history of center of mass position during simulation, natural frequencies are computed by using the Fast Fourier Transform (FFT) method.


## 6. Results and Discussions

In this section, the numerical results are presented. To this end, SWCNCs are modeled using molecular structural model in which the C–C bonds are simulated using beam elements. Besides, mass elements are used to model the carbon atoms located at the sides of beams. The MASS21 element with the concentrated mass of 1.994 3 × 10−23 g is placed at the corners of beam elements.The usage of the present analysis for the vibration behavior of SWCNCs has been validated in Figure5. This figure demonstrates the variation of the natural frequency of SWCNC with apex angle of 38.9° versus the length of SWCNC. The three-dimensional elastic BEAM4 element is used to model the bonds which has six degrees of freedom (DOFs) in each node, including three translations in x, y, and z directions and three rotations about these axes. As it is observed, the present numerical results meet an excellent agreement with those reported in [35]. It is clear that by increasing the length of the nanocone, natural frequency decreases and this fact can be generalized to all of the possible apex angles.Comparison of the obtained natural frequency using BEAM4 and BEAM188 elements with the results of MD simulations and [10] (apex angle = 38.9°).
(a)
(b)To make a comparison between the results of the present study and those predicted by MD simulation, the three-dimensional elastic BEAM188 element is also utilized for the modeling of the bonds. The utilized beam theory for BEAM188 element is based on the Timoshenko beam model, which is a first-order shear deformation theory; but for BEAM4 element, the shear deformation is neglected [43]. The values of Young’s modulus, Poisson’s ratio, and cross section radius are given as the input of this element. These values are given in (12). Figure 5 also shows that the obtained results by using BEAM188 element are in excellent agreement with those generated from MD simulation and it implies the accuracy of the present method. Comparing the given curves in Figure 5, one can conclude that, to have rational results which can trace result of MD simulations, the element BEAM188 should be utilized. So, in all of the subsequent results, the BEAM188 element is used to simulate the C–C bonds.Figure6 gives the variation of natural frequency of the SWCNCs with different apex angles versus the SWCNC length. It can be found that, by increasing the apex angle of SWCNCs, the value of natural frequency decreases. It is also seen that increasing apex angle results in decreasing the sensitivity of vibrational behavior to SWCNC length. Besides, increasing side length results in decreasing the sensitivity to length change.Figure 6
Variation of natural frequencies of the SWCNCs versus the length for different apex angles.The variations of the natural frequencies of the SWCNCs with apex angles of 38.9°, 60°, 86.6°, and 123.6° for different boundary conditions are illustrated in Figure7 versus the SWCNC length. Three different boundary conditions are considered including clamped-clamped (CC), clamped-free (CF), and simply supported-simply supported (SS). It should be noted that all of the DOFs of the nodes located at the boundary are constrained for clamped boundary conditions. For simply supported ones, only translational DOFs are constrained and the boundary nodes can revolve in rotational DOFs. The natural frequencies of the SWCNCs with CC boundary conditions are higher than those with the two other end supports. However, in higher apex angles, the dependency of natural frequencies on the boundary conditions decreases. Finally, the first five mode shapes of SWCNCs with various apex angles are shown in Figure 8. It is seen that SWCNCs with different apex angles experience different mode shapes. So, one can conclude that the mode shapes of SWCNCs depend on their apex angles.Variation of natural frequencies of the SWCNCs with the length for different apex angles and various boundary conditions: (a) apex angle = 38.9°, (b) 60°, (c) 86.6°, and (d) 123.6°.
(a)
(b)
(c)
(d)First five mode shapes of the SWCNC with various apex angles: (a) apex angle = 38.9°, (b) 60°, (c) 86.6°, and (d) 123.6°.
(a)
(b)
(c)
(d)


## 7. Conclusion

In this paper, the vibrational behavior of SWCNCs with different lengths and apex angles was studied. Besides, the SWCNCs under various boundary conditions were examined. To this end, a molecular structural model was used. Comparing the results of the model with those of MD simulations, it was shown that modeling C–C bonds by BEAM188 element can accurately predict the fundamental natural frequencies of SWCNCs. It was shown that increasing side length of a SWCNC with a constant apex angle results in decreasing the frequency. Moreover, it was seen that the SWCNCs with larger apex angles have smaller frequencies. Finally, comparing the frequencies of SWCNCs under different boundary conditions, it was observed that increasing apex angle leads to decreasing the sensitivity of vibrational behavior of SWCNCs to their end conditions.



---

*Source: 410783-2014-07-14.xml*