Propagation of Acoustic Wave inside the Carbon Nanotube : Comparative Study with Other Hexagonal Material

Carbon nanotube is a novel and more explored material. In this paper, ultrasonic acoustic velocity of the carbon nanotube has been calculated along unique axis at room temperature. For the evaluations of ultrasonic properties, secondand third-order elastic constants have been computed from Lennard-Jones interaction potential. Attenuation of ultrasonic waves due to phonon-phonon interaction is predominant over thermoelastic loss. Carbon nanotube shows the unique behavior with the chiral number. Chiral number not only affect the band gap and tube radius of the carbon nanotube but also affect the mechanical properties like stiffness, bulk modulus, shear modulus of the tube. The peculiar behavior is obtained at 55 ̊ due to their least thermal relaxation time and highest Debye average velocity. Results are also compared with other hexagonal metallic materials which present in periods and group of the periodic table. They show the optimum behavior with other hexagonal materials.


Introduction
Carbon nanotube (CNT) are cylindrical carbon molecules witch possess novel outstanding better properties (like mechanical, electrical, thermal and chemical properties, 100 times stronger than steel, best field emission emitter, can maintain current density more than 10 −9 A/cm 2 , thermal conductivity comparable to that of the diamond) which make then potential useful in a wide variety of the applications (e.g., optical nano-electronics, composite material, conducting polymers, sensor) [1,2].Nanotubes are composed of sp 2 bonds, similar to those observed in graphite and they are naturally aligned themselves into ropes held together by Vander Waals force [3].Then nono tubes are usually described using the chiral vector, , which connects two crystallographically equivalent sites (A and A¢) on a graphene sheet (where a 1 and a 2 are unit vectors of the hexagonal honeycomb lattice and n and m are integers).This chiral vector C h also defines a chiral angle õ, which is the angle between C h and the zigzag direction of the graphene sheet .Each nanotube topology is usually characterized by these two integer numbers (n, m), thus defining some nxa mxa   C peculiar symmetries such as armchair (n, n) and zigzag (n, 0) classes.The circumference vector or chiral vector describes the structure of the CNTs.The structure of the single-wall carbon nanotube (except for cap region on both ends) is specified by a vector of original hexagonal (also called hexagonal/honey Comb) lattice called the chiral vector.The unusual high stiffness of CNT makes it suitable ingredients for the application of several composite materials.CNT shows the similar properties as like other hexagonal materials (Ti, Y, Be, Mg, lanthanides ceramic etc.).Wave propagation velocity inside the chosen materials is another key factor in ultrasonic characterization, which in combination with attenuation can provide important tools in understanding, the inspectability of materials; for example, it can provide the information about crystallographic texture, grain sizes, delectability of flaws, and microstructure.Ultrasonic velocity is directly calculated from the elastic constants i.e.V  c d , where C is the elastic constant and d, the density of that particular material.Elastic constants and their variations with temperature provide useful insight into the nature of inter-atomic forces and also give the information of the stability and stiffness of the materials [4,5].
Elastic properties of materials at high pressure and high temperature are of great interest to researchers in many fields, such as physical sciences, earth sciences, and material sciences.
In the present investigation, we have chosen CNTs under the study due to their important features.In the present study, an applied approach for the calculation of second and third order elastic constants of the CNTs validating the interaction potential model was calculated.These nonlinear elastic constants have been used to calculate the ultrasonic velocity and attenuation of the CNTs for their characterization.Their properties were compared with other hexagonal material and optimized results were produced with comparative study.

Higher Order Elastic Constants
A carbon atom in graphene will assemble in a singlesheet hexagonal lattice resembling the surface of a honeycomb.This is also known as a trigonal-planar σ-bond framework, with an inter-atomic spacing, a cc , of approximately 1.41Å along the bonds that are separated by 120 degrees.The 2p electrons from all the atoms on the sheet constitute a cloud of delocalized π-orbital's surrounding the carbon cores, and these valence electrons, once excited, are responsible for conduction in graphene.Furthermore, multiple sheets of graphene may assemble in stacks, whereby two adjacent sheets are held together weakly by dispersion forces and have an inter-layer spacing of about 3.35 Å.The basis vectors generate the graphene lattice, where a = 0.142 nm is the carbon-carbon bond length.A and B are the two atoms in the unit cell of graphene.The chiral vector begins and ends at equivalent lattice points, so that the particular (n 1 ; n 2 ) tube is formed by rolling up the vector so that its head and tail join, forming a ring around the tube.The length of C h is thus the circumference of the tube, and the radius is given by the formula Many of the hexagonal crystals have axial ratio p = c/a different from the ideal value Based on the knowledge of higher order elastic constants of the material is essential for the study of the anharmonic properties of solids.Elastic constants also pro-vide insight the nature of binding forces between the atoms since they are represented by the derivatives of the internal energy.A complete set of elastic constants for materials (like second and third order elastic constant) is essential to estimate physical parameters such as Debye temperature, compressibility and acoustic anisotropy.The elastic energy density for a crystal of a cubic symmetry can be expanded up to quartic terms as shown below [6]: where C ijk and C ijklmn are the second and third order elastic constants in tensor form respectively.The second (C ij ) and third (C ijk ) order elastic constants of material are defined by following expressions.
2 , 1, , where, U is elastic energy density, e i = e ij (i or j = x, y, z, 1, , 6 i   ) is component of strain tensor.Equations ( 1) and ( 2) leads six second and ten third order elastic constants (SOEC and TOEC) for the hexagonal closed packed structured materials [7].In present approach, Lennard-Jones interaction potential {      where a 0 , b 0 are constants} was utilized due to many body interaction potential.If it is assumed that the basic structure of carbon nanotube is a hexagon then Equation (1) leads following six second-and ten third-order elastic constants.
m, n = integer quantity; b 0 = Lennard Jones parameter = 6.5 × 10 -79 Jm 7 .p = c/a: axial ratio; c is the height of the elastic constants are related rough Gruneisen parameters unit cell and "a" be the basal plane distance.The symbol "p" is the proposed parameter of axial ratio.

Ultrasonic Velocity
The second-and third-order to ultrasonic attenuation th and acoustic coupling constant.On the basis of mode of lattice vibration, there are three types of velocities (longitudinal, quasi shear and shear) in acoustical region [8].The ultrasonic velocities can be computed using calculated values of second-order elastic constants.These velocities depend on the direction of propagation of wave from the unique axis of hexagonal crystal [9].The ultrasonic velocities as a function of angle between direction of propagation and unique axis for hexagonal structured materials are given as follows [10]:

 
a parallellopiped by height c with a base which is a rhombus of side "a".A hexagonal crystal has one sixfold axis of symmetry.All belong to the crystal class P63/mmc.The ideal c/a ratio required for close packing of spheres to form the HCP structure is 1.633 (i.e.24 3 ).The tube radius is calculated from the Equation (1). Figure 1 shows the variation of tube radius with chiral number.In curve n 1 = varied but n 2 is assumed to be zero, n 1 -n 2 is the divisible of 3 (where n 1 and n 2 are the chiral number of the CNT either in zigzag structure or armchair structute).Hence, the tube is zigzag in nature and they show metallic in character.In zigzag structure, the chiral vector is oriented along the axis, hence in this case the chiral angle i.e., ζ = 0˚.In case of arm chair structure, the chiral number is equal, but it is not zero.The behaviour is again found in metallic nature, as shown in Figure 1(b).In Figure 1(c), the difference is not the divisible of 3; the found structure may be either zigzag or armchair.However the behaviour of the CNT is to be changed from metallic to semiconducting in nature, while the bond angle ζ is 30˚.The variation of the band bap with tube radius and chiral number in semiconducting CNT is shown in Figure 2.These band gaps are calculated from the following formula: cos sin where V L , V S1 and V S2 are longitudinal, qua pure shear wave ultrasonic velocities.Variab si shear and les ρ and θ represent the density of the material and angle with the unique axis of the crystal respectively.The Debye average velocity (V D ) is useful to study the Debye temperature and thermal relaxation time of the materials.The Debye temperature (T D ) is an important physical parameter of solids, which differentiate the classical and quantum behavior of phonon [11].The following expressions have been used for the evaluation of Debye average velocity and Debye temperature

 
The band gap of semiconducting CNT decreases with increase the tube radius as well as chiral number.Knowledge of the carbon nanotube dispersion relation allows one to identify whether a tube produces metallic or semiconducting behaviour, and in the latter case, determine the conduction and valence band edges and thus the band gap [12,13].Typical values for the band gap are in the range of tenths of an electron-volt to a few electron-volts, and for the previously mentioned (16, 0) tube it is ~0:62 eV.The parameters "R t " for SWCNT (10, 10) is 6.7834 Å. Hence the lattice parameter "a" becomes equal to 16.9664 Å.By applying the symmetric operation, it is also found that and C 66 can be obtained by measuring the velocity in c-direction (six fold symmetry axis).The constants C 11 , C 33 and C 44 have been determined from the direct relation, 2 , where V is the relevant sound velocity.They are related with Equation (7).Seven velocity measurements enable us to determine all the seven second-order

Results and Discussion
ure, conductivity and e, it is necessary to In order to examine the band struct wave propagation inside the nanotub derive its ε-k relation.The hexagonal unit cell is charac-   elastic constants.This is metallic in character and ultra ed by sing lattice parameter and set of Equations ( 2) for (10, sonic attenuation is due to phonon-phonon interaction.The SOECs and TOECs of SWCNT are calculat u 10) chiral system, their results are tabulated in Table 1.Furthermore, the following physical parameters have been also determined from set of equation.
using s ond order elastic constants, which are tabulated in Table 1.The values of six independent elastic const s along with dependent variable were evaluated by calculating on different deformation All these elastic constants are positive and satisfy the well-known Born's criteria for mechanically stable hexagonal crystals: C > 0, (C 11 − C 12 ) > 0, C 44 > 0, and ec ant the stress tensors s applied to the equilibrium lattice of the hexagonal unit cell.The beautiful results a in constant with chiral number is shown in Figure 3.With increase the chiral number, the elastic constant decrease from higher to sudden lower values.Beyond (6,10), the values becomes constant with increase the chiral num For zigzag structure, the maximum values are at higher chiral number, but the optimum and useful values are found from n 1 = 12 to n 1 = 21.re tabulated in Table 2 for different chiral number.The ultrasonic velocities are evaluated using elastic constants and density which in turn provides Debye average velocity in all chosen materials.The variations of all the three types of velocity are shown in Figure 4.With increase the angle from, the longitudinal velocity increases from 0˚ to 60˚ but again decrease with increase the angle of propagation.But shear velocity direction-1 show the same effect in the opposite direction.is the p-block elements.For investigating the wave propagation in the material, the density decides its nature.The unit cell parameters "a" and "p" for Gd, Tb, Dy, Ho, Er and Tm are 3.63Å, 3.60Å, 3.59Å, 3.58 Å, 3.56Å, 3.54Å and 1.592, 1.583, 1.574, 1.570, 1.570, 1.571 respectively.The value of m, n and b 0 for these lanthanides metals are 6, 7 and 2.3 × 10 −64 erg cm 7 correspondingly.Three metals La, Pr and Nd have the DHCP (double HCP) structure where the stacking of close packed planes follows the order ABACABAC instead of the usual sequence ABAB [14].Consequently, close packing under these conditions leads to an ideal c/a ratio of 3.3256.Overall, most of the metals were lower and within −0.1 of their ideal c/a ratio, except for Cd and Zn which exceeded the ratio by 0.223 and 0.194, respectively.In case of transition metal like Ti, Y, Zr, and Hf, the SOEC shows the irregular behavior due involvement of 4f orbital.In case of heavy rare earth elements, the SOEC are found to be increasing from metal Gd to Tm.Thus the thermal softening behaviors are obtained as we move from Tm to Gd.These metals have low bulk modulus and higher order elastic constants in comparison to other hexagonal structured material [7][8][9][10]12,13].This reveals that large strain will develop for a given deforming force in these materials.Figure 6 shows the variation of average Debye velocity vs. angle with unique axis of crystal and compare with alkaline earth metal like Be and Mg (because they possess HCP crystals structure).The velocity of propagation in medium depends upon the elastic Average Velocity (m/s) constant, and density of the materials [15].
The density of the carbon nanotube is much similar to these metals but their second order elastic constant four oc same crystals structure but they show dramatic effect on the propagation of acoustic wave inside the materials.In this material, the average Debye velocity decreases with the angle and attain the minimum value around 45˚ -55˚ and then increase with the angle, as shown in Figure 8. Hf also possesses higher density (9 -10 times heavy from the carbon nanotube) as compared to carbon nanotube.Due to higher density, they show lower propagating velocity inside the materials.Due to involvement of 4f orbital in Hf and also expected double crystal structure the propagating velocity along axis, they show negative effect from the Ti, and Y. Re is again the member of the 5d series and exhibits the same crystal structure, but thei city first decreases and attains the minimum values and and very high elastic constant as compared to the known to five times less than "Be" but similar to Mg. it is evident from the Figure 6 that the propagating velocity inside the Carbon nanotude is two times less than the Be crystal structure.Ti, Co, and Zn belong in 3d series and Y and Zr in 4d series but they possess HCP crystal structure.The variation of the average velocity with propagating angle along the unique axis is presented in the Figure 7.The density as well as second order elastic constant of the Ti and Y is high as carbon nanotube but the velocity of propagation along the crystals axis shows negative effect.To compare with carbon nanotube the average velocity in Ti decrease with increasing angle, as evident from Figure 7(a).Y show similar effect as Ti, compared with the Figure 7(b).Hf is the member of d-bl k (5d) element like Ti, Y and also possesses the r trends for propagating velocity inside the materials show positive effect from the Hf , Ti and Zr. Figure 9 shows the variation of average Debye velocity with the propagating axis.It is evident from the Figure 9 that the velo then increase with increase the angle.Due to their high density (19 -20 times heavy from the carbon nanotube) In periodic table and known HCP crystal structure, Os is the heaviest element among all the elements.Figure 10 shows the variation of average Debye velocity with propagating angle along the crystal axis.It can be seen from the Figure 10 that the velocity decrease with increase the angle, beyond attending the minimum velocity they again increase.Out of the various materials studied in this research work, they possess the higher values because of the higher elastic constant along with the maximum density (18 -20 times more than the carbon nanotube).
The group Os, Ru, Re, Be and Co exhibit the highest moduli of all the HCP metals.All show a pronounced maximum value of E on the basal plane, for which θ is zero, and a tendency by Ru, Re and Co to exhibit a min cating G is highest on the basal plane.The elastic constants are important since they are related to hardness and are used for the determination of the ultrasonic velocity.The plot of V D (Figure 11(a)) for the propagation of wave at different angles with unique axis implies that these metals have slight variation in velocity for each direction of propagation and have maxima at θ = 55˚.Fig. 11b shows the variation of the average velocity with crystals axis or the rare earth elements along with carbon nanotube.Lanthanides possess higher velocity with carbon nanotube because of similar elastic constant and heavier density.Carbon nanotube is the allotrope of the carbon and metal in HCP crystals structure expect Os, they exhibit lower velocity because of density effect, evaluated from Equation (7).
imum between 0 < θ < 90.With the exception of Be, the G-behaviour exhibits a minimum when θ is zero and 90˚, and a maximum for 0o < θ < 90˚.In the case of Be, a maximum G occurs when θ is zero, indi Tl belongs to the same family.The pπ-pπ interacti n provides the propagation of the velocity inside the metal.
In other case, sσ-sσ, dπ-sπ dπ-dπ, dπ-fπ and fπ-fπ are the responsible factor for the propagating the acoustic wave inside the material.In case of metal either from the s, d or f-block, they bear the maximum lattice vibration from the initial stage due to interaction with the acoustic wave.The dominant mechanism for ultrasonic attenuation in these metals is phonon-phonon interaction.But the carbon nanotube shows the adverse effect because of third order elastic constan ebye average velocity.5) Therefore, in carbon nanotube, the bonding between the atoms along the a-b plane is stronger than that between the at meter, but also second and third order parameter gives the same indication.The anisotropic properties of SWCNT astic constant are well related with the ultrasonic velocities as well as with the higher order elastic constants.
6) The highest change in Debye average was observed at 50˚ -55˚along the unique axis.
7) Due to change in bonding from metal to SWCNT, the observed changes in velocity with different type of the metal along with CNT are different at the different angle.Their values depend upon the phonon-phonon vibration in the lattice.
8) In metallic material, phonon-phonon vibration is easily achieved due acoustic wave interaction as compared to CNT because of less tightly packing.

1 =6, n 2 =varied=not equal to n 1 Figure 1 .Figure 2 .
Figure 1.Variation of Tube radius with chiral number for different values of n 1 and n 2.

Figure 5 .
Figure 5. Variation average Debye velocity vs. angle with unique axis of crystal.
The other parameter was calculated from following formula by the evaluation of the C 11 , C 12 , C 13 , C 33 , C 44 , and C 66 :