# Analysis of equation of state for carbon nanotubes.

1. IntroductionAfter their discovery by Iijima [1] carbon nanotubes have got a lot of attention because of their unique electronic and mechanical properties, but the most important property of carbon nanotubes is their elastic response to the external force or stress [2,3]. As a one-dimensional structure, carbon nanotubes can be thought of as one sheet or multiple sheets of graphene rolled into a cylinder. Single-walled carbon nanotube bundles typically consist of several nested tubes, each like a graphene sheet bent into the cylindrical form with an overall diameter of a few nanometers. According to different chiral angles, single-walled carbon nanotubes can be classified into zigzag ([theta] = 0[degrees]), armchair ([theta] = 30[degrees]), and chiral tubes (0[degrees] < [theta] < 30[degrees]) [4]. The Young's modulus and Poisson's ratio of nanotube have been studied by the previous workers [5-9], while the hardness as one of the most important parameters characterizing the mechanical properties of single-walled carbon nanotubes has been intensively studied during the last decade [6,10].

High pressures which are encountered from deep down the earth to the astrophysical objects may cause many effects such as compression, pressure ionization, modification in electronic properties, phase changes, and several phenomenons in applied fields [11]. For this, pressure versus volume relations of condensed matter known as equation of state is a vital input. Many equations of state exist in the literature, but still there is a need to judge on their suitability under whole range of compressions as most of them give the same result under small compression.

In high pressure, generally used theory is the finite strain theory which means the theory due to Birch [12]. However, Birch's theory rates no more than a passing mention as discussed in detail by Stacy [13]. The attention has also been given to the theory based on atomic potential as presented by Rydberge [14]. The Mie Gruneisen theory was followed by Suzuki et al. [15,16] and reported what became known as Suzuki formulation [16-18]. The Gruneisen theory of thermal expansion as formulated by Born and Huang [19] has been used by Shanker et al. [20]. These authors included a higher order term for the change in the expansion of potential energy. Using some thermodynamical relation Tait formulated a relation [21-23] known as UTE.

In the present work, we have investigated the mechanical properties especially volume compression (V/[V.sub.0]) of nanotube bundles and individual carbon nanotubes using widely used Suzuki, Shanker, and usual Tait equation of state (EOS) [15-23]. We have then compared our calculated results with the experimental values as well as with the results obtained for graphite. For input parameters the values of bulk modulus and its first order pressure derivative of nanotube bundle ([K.sub.0] = 37 GPa, [K.sup.1.sub.0] = 11) and individual tube within the bundle ([K.sub.0] = 230 GPa, [K.sup.1.sub.0] = 4.5) have been taken from the work of Reich et al. [24] who have used local-density approximation of density-functional theory to calculate these values. For graphite the values of bulk modulus and its first order pressure derivative are taken to be [K.sub.0] = 33.8 GPa and [K.sup.1.sub.0] = 8.9, as measured by Hanfland et al. [25] from X-ray diffraction studies.

2. Method of Analysis

Based on Mie-Gruneisen Debye theory Suzuki et al. [15,16] have derived an EOS, known as the Suzuki equation [18,26], to study the change in volume V/[V.sub.0] with temperature T, along an isobar. This Suzuki equation has been widely used in the literature [15-18], and its mathematical form reads as follows [27]:

V/[V.sub.0] = [[1 + 2k - (1 - 4k[E.sub.Th]/Q)].sup.1/2]/2k, (1)

where k = ([K'.sub.0] - 1)/2, [E.sub.Th] is the thermal energy, Q = [K.sub.0] [V.sub.0]/[[gamma].sub.0], and [[gamma].sub.0] is the Gruneisen ratio. [K.sub.0], [K'.sub.0] are the isothermal bulk modulus and its first order pressure derivative, respectively. The subscript o refers to their value at P = 0 and T = 300 K. Using the well known relation [27] [[gamma].sub.0] [E.sub.Th]/[V.sub.0] = [P.sub.Th], (1) reads as follows [28]:

V/[V.sub.0] - 1 = 1 - [[1 - 2(([K'.sub.0] - 1)/[K.sub.0])[P.sub.Th]].sup.1/2]/[K'.sub.0] - 1, (2)

where [P.sub.Th] is the thermal pressure.

Shanker et al. [20] have argued that if P is not equal to zero, (2) may be rewritten as follows:

V/[V.sub.0] - 1 = 1 - [[1 - 2(([K'.sub.0] - 1)/[K.sub.0])([P.sub.Th] - P)].sup.1/2]/[K'.sub.0] - 1. (3)

Now, when thermal pressure is zero ([P.sub.Th] = 0), (3) gives the following simple relation:

V/[V.sub.0] - 1 = 1 - [[1 + 2(([K'.sub.0] - 1)/[K.sub.0])P].sup.1/2]/[K'.sub.0] - 1. (4)

This is Suzuki EOS.

Shanker et al. [20] have reported an isobaric EOS, using the Gruneisen theory of thermal expansion as formulated by Born and Huang [19] whose mathematical form reads as follows:

V/[V.sub.0] - 1 = 1 - [[1 - 2(([K'.sub.0] + 1)/[K.sub.0])[P.sub.Th]].sup.1/2]/[K'.sub.0] + 1. (5)

It has been argued by Shanker et al. [20] that when P is not equal to zero, (5) may be rewritten as follows:

V/[V.sub.0] - 1 = 1 - [[1 - 2(([K'.sub.0] - 1)/[K.sub.0])[P.sub.Th - P]].sup.1/2]/[K'.sub.0] - 1. (6)

When thermal pressure is zero ([P.sub.Th] = 0), (6) gives

V/[V.sub.0] - 1 = 1 - [1 + 2(([K'.sub.0] + 1)/[K.sub.0])P].sup.1/2]/[[K'.sub.0] + 1. (7)

This is Shanker EOS.

The usual Tait equation of state is obtained by assuming the fact that the product of the thermal expansion coefficient (a) and the bulk modulus ([K.sub.T]) is constant under the effect of pressure [29], that is,

[alpha]K = constant. (8)

Differentiation of (8) with respect to volume at constant temperature gives

[alpha][(dK/dV).sub.T] + K[(d[alpha]/dV).sub.T] = 0. (9)

Anderson-Gruneisen parameter is defined as

[[delta].sub.T] = V/[alpha][(d[alpha]/dV).sub.T], (10)

where [[delta].sub.T] is Anderson-Gruneisen parameter at constant temperature.

From (9)and (10), we get

[[delta].sub.T] = V/[alpha][(d[alpha]/dV).sub.T] = - V/K[(dK/dV).sub.T]. (11)

Assuming [[delta].sub.T] to be independent of V,

[[delta].sub.T] = [(dK/dP).sub.T] = [K.sup.1.sub.0]. (12)

Anderson-Gruneisen parameter [[delta].sub.T] and [eta] = V/[V.sub.0] are related by the following relation [28]:

([[delta].sub.T] + 1)/[eta] = A, (13)

where A is a constant for a given solid.

In view of (13), (11) can be written as

dK/K = [- A/[V.sub.0] + 1/V]dV. (14)

Integrating the previous equation, we get

K/[K.sub.0] = V/[V.sub.0] exp A [1 - V/[V.sub.0]], (15)

where

K = -V[(dP/dV).sub.T]. (16)

In view of (16), (15) is written as

K/[K.sub.0] exp A [1 - V/[V.sub.0]]dV = -dP. (17)

The integration of (17) gives

p = [K.sub.0]/A [exp A(1 - V/[V.sub.0]) - 1]. (18)

Here, [K.sub.0] is the bulk modulus at the zero pressure, and the constant A is determined from the initial conditions, namely, at V = [V.sub.0], A = [[delta].sup.0.sub.T] + 1.

Substituting the value of A in (18) and taking the natural log, we get the following final form of usual Tait equation of state:

V/[V.sub.0] - 1 = - ln [1 + (([K.sup.1.sub.0] + 1)/[K.sub.0])P]/[K.sup.1.sub.0] + 1, (19)

here the approximation [[delta].sup.0.sub.T] [approximately equal to] [K.sup.1.sub.0] is used.

To test the validity of this usual Tait's equation of state we have employed it to predict the high pressure compression behavior of carbon nanotubes and graphite along with the Suzuki and Shanker EOSs.

Within the framework of UTE, the expression for the isothermal bulk modulus can be written as [30,31]

[K.sub.T]/[K.sub.0] = V/[V.sub.0] {1 + [K.sup.1.sub.0] + 1/[K.sub.0]P}. (20)

3. Results and Discussion

Present analysis is based on the three widely used formulations, namely, Suzuki, Shanker, and usual Tait EOSs. We have calculated the values of compression (V/[V.sub.0]) at different pressures by using (4), (7), and (19)for individual CNT as well as nanotube bundles. Table 1 shows the variation of (V/[V.sub.0]) with pressure for individual as well as nanotube bundles. In Table 2, we have reported the values of the compression for graphite which may be treated as a bulk material. Thus we have analyzed the validity of equation of state for nanomaterials as well as bulk materials. For high pressure compression behavior, a graphical presentation between the calculated and experimental values has been shown in Figures 1 and 2. From tables and figures, it is clear that the results obtained by Suzuki formulation do not agree with the experimental data [24,25] while Shanker formulation improves the results slightly. Moreover, it seems that at low compressions all equation give the similar results. On the other hand usual Tait EOS (19) gives the results which are in close agreement to the experimental data [24,25] for the entire range of pressure. For usual Tait EOS the agreement between the theoretically calculated and experimental values demonstrates that the product of bulk modules and the coefficient of volume thermal expansion remain constant for carbon nanotubes. While the failure of Suzuki and Shanker formulation may be attributed to the fact that their analysis is based on approximation. The Suzuki model is based on the Taylor series expansion of potential energy, and during this expansion higher order terms beyond second order have been neglected. This approximation must introduce serious errors as discussed in detail by Wang and Reeber [26]. Inspired with the Suzuki model, Shanker et al. [20] have also neglected higher order terms in deriving (7).

In the present study from Tables 1 and 2 it is found that the graphite is more compressible than the nanotube bundle which is formed by rolling graphene sheets. We further extended usual Tait's EOS to calculate bulk modulus and the variation of isothermal bulk modulus (KT) with pressure for nanotube bundle and individual tube within the bundle which are shown in Figures 3 and 4, while calculated values are given in Table 3. From Figures 3 and 4, it is clear that the isothermal bulk modulus increases with pressure and the compression decreases with pressure in these nanotubes and thus follows the same trend as in the bulk solids. The variation of bulk modulus with pressure could not be compared with the experimental values, as to the best of our knowledge the experimental data about the pressure dependence of bulk modulus for carbon nanotubes is not available in the literature. However, in future our calculated values may stimulate some experimental investigations regarding high pressure behavior of bulk modulus for carbon nanotubes.

4. Conclusions

On the basis of overall descriptions, it may thus be concluded that the Suzuki formulation fails to explain the high pressure compression behavior of carbon nanotubes Shanker formulation improves the result of Suzuki formulation, but among the three formulations, the usual Tait EOS is more suitable to explain the elastic properties of nanomaterials. The reason for its success lies in the fact that the product of bulk modules and the coefficient of volume thermal expansion remain constant for carbon nanotubes.

http://dx.doi.org/10.1155/2013/639068

References

[1] S. Iijima, "Helical microtubules of graphitic carbon," Nature, vol. 354, no. 6348, pp. 56-58, 1991.

[2] M. S. Dresselhaus, G. Dresselhaus, and P. C. Eklund, Science of Fullerenes and Carbon Nanotubes, Academic, New York, NY, USA, 1996.

[3] B. Vigolo, A. Penicaud, C. Coulon et al., "Macroscopic fibers and ribbons of oriented carbon nanotubes," Science, vol. 290, no. 5495, pp. 1331-1334, 2000.

[4] M. S. Dresselhaus, G. Dresselhaus, and R. Saito, "Physics of carbon nanotubes," Carbon, vol. 33, no. 7, pp. 883-891, 1995.

[5] J. P. Lu, "Elastic properties of single and multilayered nanotubes," Journal of Physics and Chemistry of Solids, vol. 58, no. 11, pp. 1649-1652, 1997.

[6] J. P. Lu, "Elastic properties of carbon nanotubes and nanoropes," Physical Review Letters, vol. 79, no. 7, pp. 1297-1300, 1997.

[7] E. Hernandez, C. Goze, and A. Rubio, "Elastic properties of single-wall nanotubes," Applied Physics A, vol. 68, pp. 287-292, 1999.

[8] D. Sanchez-Portal, E. Artacho, J. M. Soler, A. Rubio, and P. Ordejon, "Ab initio structural, elastic, and vibrational properties of carbon nanotubes," Physical Review B, vol. 59, no. 19, pp. 12678-12688, 1999.

[9] G. van Lier, C. van Alsenoy, V. van Doren, and P. Geerlings, "Ab initio study of the elastic properties of single-walled carbon nanotubes and graphene," Chemical Physics Letters, vol. 326, no. 1-2, pp. 181-185, 2000.

[10] B. I. Yakobson, C. J. Brabec, and J. Bernholc, "Nanomechanics of carbon tubes: instabilities beyond linear response," Physical Review Letters, vol. 76, no. 14, pp. 2511-2514, 1996.

[11] J. S. Schilling, "The use of high pressure in basic and materials science," Journal of Physics and Chemistry of Solids, vol. 59, no. 4, pp. 553-568, 1998.

[12] F. Birch, "Elasticity and constitution of the Earth's interior," Journal of Geophysical Research, vol. 57, no. 2, pp. 227-286, 1952.

[13] F. D. Stacy, "Finite strain, thermodynamics and the earth's core," Physics of the Earth and Planetary Interiors, vol. 128, pp. 179-193, 2001.

[14] R. Rydberge, "Graphical representation of results of band spectroscopy," Zeitschrift fur Physik, vol. 73, pp. 376-385, 1931.

[15] I. Suzuki, "Thermal expansion of periclase and olivine, and their anharmonic properties," Journal of Physics of the Earth, vol. 23, no. 2, pp. 145-159, 1975.

[16] I. Suzuki, S. Okajima, and K. Seya, "Thermal expansion of single-crystal manganosite," Journal of Physics of the Earth, vol. 27, no. 1, pp. 63-69, 1979.

[17] O. L. Anderson, Equation of State for Geophysics and Ceramic Science, Oxford University Press, Oxford, UK, 1995.

[18] G. Helfrich, "Practical use of Suzuki's thermal expansivity formulation," Physics of the Earth and Planetary Interiors, vol. 116, pp. 133-136, 1999.

[19] M. Born and K. Huang, Dynamical Theory of Crystal Lattice, Oxford University Press, Oxford, UK, 1995.

[20] J. Shanker, S. S. Kushwah, and P. Kumar, "Theory of thermal expansivity and bulk modulus for MgO and other minerals at high temperatures," Physica B, vol. 233, no. 1, pp. 78-83, 1997.

[21] M. Taravillo, V. G. Baonza, J. E. F. Rubio, J. Nuez, and M. Caceres, "The temperature dependence of the equation of state at high pressures revisited: a universal model for solids," Journal of Physics and Chemistry of Solids, vol. 63, no. 9, pp. 1705-1715, 2002.

[22] M. Singh, P. P. Singh, B. R. K. Gupta, and M. Kumar, "Temperature and pressure dependence of elastic constants," High Temperatures, vol. 33, no. 2, pp. 199-206, 2001.

[23] D. Kandpal, K. Y. Singh, and B. R. K. Gupta, "On pressure dependence of the relative compression (V/V0) at room temperature for the solids: copper and lead as prototypes," Indian Journal of Physics A, vol. 78, no. 3, pp. 393-395, 2004.

[24] S. Reich, C. Thomsen, and P. Ordejon, "Elastic properties of carbon nanotubes under hydrostatic pressure," Physical Review B, vol. 65, no. 15, Article ID 153407, 4 pages, 2002.

[25] M. Hanfland, H. Beister, and K. Syassen, "Graphite under pressure: equation of state and first-order Raman modes," Physical Review B, vol. 39, no. 17, pp. 12598-12603, 1989.

[26] K. Wang and R. R. Reeber, "The role of defects on thermophysical properties: thermal expansion of V, Nb, Ta, Mo and W," Materials Science and Engineering R, vol. 23, no. 3, pp. 101-137, 1998.

[27] M. Kumar, "A comparative study of Suzuki and Kumar formulations for thermal expansivity," Indian Journal of Pure and Applied Physics, vol. 42, no. 1, pp. 67-70, 2004.

[28] O. L. Anderson, D. G. Isaak, and H. Oda, "High-temperature elastic constant data on minerals relevant to geophysics," Reviews of Geophysics, vol. 30, no. 1, pp. 37-90, 1992.

[29] J. Shanker and M. Kumar, "Thermodynamic approximations in high-pressure and high-temperature physics of solids," Physica Status Solidi B, vol. 179, no. 2, pp. 351-356, 1993.

[30] H. Schlosser, J. Ferrante, and P. Vinet, "Pressure dependence of the melting temperature of metals," Physical Review B, vol. 40, no. 9, pp. 5929-5935, 1989.

[31] J. R. MacDonald, "Some simple isothermal equations of state," Reviews of Modern Physics, vol. 38, no. 4, pp. 669-679, 1966.

Jeewan Chandra, (1) Pooja Kapri Bhatt, (2) and Kuldeep Kholiya (3)

(1) Department of Physics, National Institute of Technology Uttarakhand, Srinagar 247174, India

(2) Kumaon University, Nainital, Uttarakhand 263002, India

(3) Department of Applied Science, B. T. Kumaon Institute of Technology Dwarahat263653, India

Correspondence should be addressed to Jeewan Chandra; jeewan25may@rediffmail.com

Received 9 May 2013; Revised 6 July 2013; Accepted 14 July 2013

Academic Editor: Bobby G. Sumpter

TABLE 1: Volume compression with pressure for carbon nanotube bundles and individual carbon-nanotubes using different EOSs. P (GPa) v/[v.sub.0] v/[v.sub.0] v/[v.sub.0] v/[v.sub.0] (19) (7) (4) (Exp [24]) bundle bundle bundle bundle 1.0 0.9766 0.9763 0.9759 0.9769 1.5 0.9670 0.9663 0.9654 0.9638 2.0 0.9583 0.9570 0.9557 0.9513 2.5 0.9505 0.9484 0.9467 0.9425 3.0 0.9434 0.9403 0.9381 0.9400 3.5 0.9368 0.9326 0.9299 0.9350 4.0 0.9307 0.9253 0.9222 0.9238 4.5 0.9250 0.9184 0.9147 0.9213 5.0 0.9197 0.9117 0.9076 0.9188 5.5 0.9147 0.9052 0.9007 0.9088 6.0 0.9100 0.8990 0.8940 0.9063 6.5 0.9055 0.8930 0.8875 0.8981 7.0 0.9013 0.8872 0.8813 0.8963 7.5 0.8972 0.8815 0.8752 0.8900 8.0 0.8934 0.8760 0.8693 0.8875 P (GPa) v/[v.sub.0] v/[v.sub.0] v/[v.sub.0] v/[v.sub.0] (19) (7) (4) (Exp [24]) individ. individ. individ. individ. 1.0 0.9957 0.9957 0.9957 0.9944 1.5 0.9936 0.9936 0.9936 0.9925 2.0 0.9915 0.9915 0.9914 0.9900 2.5 0.9894 0.9894 0.9893 0.9888 3.0 0.9874 0.9874 0.9872 0.9863 3.5 0.9854 0.9854 0.9852 0.9838 4.0 0.9834 0.9834 0.9831 0.9825 4.5 0.9814 0.9814 0.9811 0.9800 5.0 0.9795 0.9794 0.9790 0.9781 5.5 0.9776 0.9775 0.9770 0.9776 6.0 0.9756 0.9754 0.9750 0.9757 6.5 0.9737 0.9736 0.9730 0.9738 7.0 0.9719 0.9718 0.9710 0.9720 7.5 0.9700 0.9699 0.9691 0.9702 8.0 0.9682 0.9680 0.9671 0.9683 TABLE 2: Volume compression with pressure for graphite. P (GPa) V/[V.sub.0] V/[V.sub.0] V/[V.sub.0] V/[V.sub.0] (Exp [25]) (19) (7) (4) 0 1 1 1 1 2.0870 0.9618 0.9518 0.9504 0.9487 3.0434 0.9507 0.9356 0.9325 0.9296 4.8696 0.9204 0.9105 0.9027 0.8975 6.9565 0.8961 0.8878 0.8735 0.8656 8.3478 0.8809 0.8751 0.8559 0.8463 9.0434 0.8737 0.8693 0.8475 0.8372 10.434 0.8665 0.8586 0.8316 0.8197 11.913 0.8566 0.8483 0.8157 0.8022 13.913 0.8500 0.8359 0.7955 0.7798 TABLE 3: Variation of bulk modulus with pressure. P (GPa) [K.sub.T] bundle KT individual 1 47.853 23.487 1.5 53.183 236.725 2 58.459 238.952 2.5 63.685 241.166 3 68.866 243.394 3.5 74.007 245.611 4 79.109 24.817 4.5 84.176 250.012 5 89.209 25.21 5.5 94.212 254.394 6.0 99.186 256.583 6.5 104.132 258.761 7 109.053 260.955 7.5 113.948 263.112 8 118.820 265.287

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Chandra, Jeewan; Bhatt, Pooja Kapri; Kholiya, Kuldeep |

Publication: | Journal of Nanotechnology |

Article Type: | Report |

Date: | Jan 1, 2013 |

Words: | 3476 |

Previous Article: | Optimizing the processing conditions for the reinforcement of epoxy resin by multiwalled carbon nanotubes. |

Next Article: | Some fundamental aspects of mechanics of nanocomposite materials and structural members. |

Topics: |