Graphene
Vol.04 No.04(2015), Article ID:59477,8 pages
10.4236/graphene.2015.44008

The Surface Reactivity of Pure and Monohydrogenated Nanocones Formed from Graphene Sheets

Ahlam A. El-Barbary1,2*, Mohamed A. Kamel1, Khaled M. Eid1,3, Hayam O. Taha1, Rasha A. Mohamed1, Mohammed A. Al-Khateeb1

1Physics Department, Faculty of Education, Ain Shams University, Cairo, Egypt

2Physics Department, Faculty of Science, Jazan University, Jazan, KSA

3Bukairiayh for Science, Qassim University, Qassim, KSA

Email: *ahla_eg@yahoo.co.uk

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 29 July 2015; accepted 6 September 2015; published 9 September 2015

ABSTRACT

A systematic computational study of surface reactivity for pure and mono-hydrogenated carbon nanocoes (CNCs) formed from graphene sheets as functions of disclination angle, cone size and hydrogenation sites has been investigated through density functional (DFT) calculations and at the B3LYP/3-21G level of theory. Five disclination angles (60˚, 120˚, 180˚, 240˚ and 300˚) are applied and at any disclination angle four structures with different sizes are studied. For comparison, pure and mono-hydrogenated boron nitride nanocones (BNNCs) with disclination angles 60˚, 120˚, 180˚, 240˚ and 300˚ are also investigated. The hydrogenation is done on three different sites, HS1 (above the first neighbor atom of the apex atoms), HS2 (above one atom of the apex atoms) and HS3 (above one atom far from the apex atoms). Our calculations show that the highest surface reactivity for pure CNCs and BNNCs at disclination angles 60˚, 180˚ and 300˚ is 23.50 Debye for B41N49H10 cone and at disclination angles 120˚ and 240˚ is 15.30 Debye for C94H14 cone. For mono-hydroge- nated CNCs, the highest surface reactivity is 22.17 Debye for C90H10-HS3 cone at angle 300˚ and for mono-hydrogenated BNNCs the highest surface reactivity is 28.97 Debye for B41N49H10-HS1 cone when the hydrogen atom is adsorbed on boron atom at cone angle 240˚.

Keywords:

Carbon Nanocones, Boron Nitride Nanocones, DFT, Surface Reactivity, Hydrogenation

1. Introduction

The hallow shape and the curvature of CNTs that are not present in bulk graphite make them easier to storage hydrogen inside and outside the tube surface. At the same line, one of the most important applications of CNCs (representing the fifth allotropic form of carbon) is hydrogen storage. Hydrogen storage on CNCs has been explored theoretically [1] -[10] and experimentally [11] -[13] . In addition, the fabrication of CNCs has attracted attention because of their mechanical stability and sharp tip structures [14] [15] . Also, CNCs can be used as practical phononic devices, scanning probe tips and electron field emitters [16] -[18] . The field emission properties of cone-shaped CNTs whose tip are stacked by CNCs have been investigated [19] . Furthermore, decorating the network of carbon nanostructures by nitrogen and boron has been studied [20] -[22] . The main problem of hydrogen storage is that a hydrogen fuel-cell must be of comparable weight and size to current gas tanks for the technology to gain acceptance where the determined hydrogen storage density of 9 wt% is required. Therefore, the main idea behind this work is to investigate the ability of CNCs and BNNCs for hydrogen storage via surface reactivity study. By increasing the surface reactivity, the ability of hydrogen storage is expected to be increased. To provide a better understanding of surface reactivity, pure CNCs and BNNCs, in addition to mono- hydrogenated CNCs and BNNCs are studied. Five disclination angles 60˚, 120˚, 180˚, 240˚ and 300˚ are applied on graphene sheet and the cones 60˚, cones 120˚, cones 180˚, cones 240˚ and cones 300˚ are obtained. At any disclination angle for CNCs and BNNCs, four structures with different sizes are studied in order to understand the effect of the size on the surface reactivity of NCs, except for BNNCs with disclination angles 60˚, 180˚ and 300˚ eight structures for each disclination angle are studied. Also, the hydrogenation is done on three different sites, HS1 (above the first neighbor atom of the apex atoms), HS2 (above one atom of the apex atoms) and HS3 (above one atom far from the apex atoms). Hopefully, these results could be useful for testing the capability of CNCs and BNNs as hydrogen storage systems.

2. Methodology

Density Functional Theory (DFT) calculations have been performed employing the B3LYP exchange-correla- tion functional [23] [24] and the 3-21G standard basis set as implemented in the Gaussian 03W program [25] [26] . The samples of study are pure and mono-hydrogenated CNCs and BNNCs with five disclination angles, 60˚ (a five-membered ring at the apex), 120˚ (a four-membered ring at the apex), 180˚ (a three-membered ring at the apex), 240˚ (a bicyclic system at the apex) and 300˚ (a complex with a three-membered and four-membered rings at the apex), see Figure 1. The hydrogenation is done on three different sites, HS1 (above the first neighbor

Figure 1. Optimized structures of CNCs and BNNCs: (a) and (f) cones 60˚, (b) and (g) cones 120˚, (c) and (h) cones 180˚, (d) and (i) cones 240˚, (e) and (j) cones 300˚. Grey atoms represent carbon atoms, blue atoms represent nitrogen atom and pink atoms represent boron atoms. Arrows refer to the bonds between connected atoms edges (boron atoms, model M2).

atom of the apex atoms), HS2 (above one atom of the apex atoms) and HS3 (above one atom far from the apex atoms). To avoid the dangling effects, the hydrogen atoms have been used to saturate the ending atoms for CNCs and BNNCs. All the atomic geometries of pure and mono-hydrogenated CNCs and BNNCs have been allowed to fully relaxation during the optimization processes.

3. Results

3.1. Geometric Structures of CNCs and BNNCs

The CNCs are consisting of curved graphite sheets formed as open cones and are constructed by cutting out sectors of n × 60˚ (n = 1 - 5) from the flat sheet of graphene and connecting the edges. In the present work, we have investigated cone 60˚ (a five-membered ring at the apex), cone 120˚ (a four-membered ring), cone 180˚ (a three-memberedring), cone 240˚ (a bicyclic system), and cone 300˚ (a complex with a three-membered and four- membered rings at the apex), see Figures 1(a)-(e). At the same line the BNNCs are constructed, see Figures 1(f)-(j). For all CNCs and BNNCs with disclination angles 120˚ and 240˚, there is only one model. However, for BNNCs with disclination angles 60˚, 180˚, and 300˚ there are two models (M1 and M2), resulting from the connected atoms edges can be either nitrogen atoms (named M1) or can be boron atoms (named M2).

3.2. Surface Reactivity of Pure CNCs and BNNCs

The surface reactivity of pure CNCs and BNNCs at different disclination angels n × 60˚ (n = 1 to 5) are studied. Also, at any disclination angle the effect of the cone size is tested through working on four different structures. By this way, one can test the effects of size, disclination angle and types of NCs on the surface reactivity of NCs.

The surface reactivity of CNCs and BNNCs at disclination angle 60˚ is investigated and is listed in Table 1. The dipole moments are calculated and are used as indicator for the surface reactivity [27] [28] where the high values of the dipole moments indicate the high surface reactivity and vice versa. As shown in Table 1, the surface reactivity of cone 60˚ is increased by increasing the cone size. Also, it found that the B82N87H30 (M1) cone 60˚ possesses the highest surface reactivity (13.1 Debye), followed by C169H30 cone 60˚ (12.1 Debye) and the smallest surface reactivity is for B87N82H30 (M2) cone 60˚ (7.1 Debye). In other words, by increasing the number of nitrogen atoms (in case of the type of connected atoms edges is nitrogen atoms) and decreasing the number of boron atoms, the surface reactivity of BNNCs is increased.

Table 2 shows the calculated surface reactivity of CNCs and BNNCs at disclination angle 120˚. It is found that the calculated surface reactivity of cones 120˚ is increased by increasing the cone size. Also, it is noticed that the surface reactivity of the C136H24 cone 120˚ (13.6 Debye) is higher than the surface reactivity of B68N68H24 cones 120˚ (9.2 Debye). This indicates that when the number of boron atoms is equal to number of nitrogen atoms (i.e. type of connected atoms edges is both of nitrogen and boron atoms), the surface reactivity of BNNCs is reduced. From Table 1 and Table 2, it is clear that by increasing the disclination angle (from cones 60˚ to cones 120˚), the surface reactivity is increased for both of CNCs and BNNCs.

The surface reactivity of CNCs and BNNCs at disclination angle 180˚ is listed in Table 3. As shown in Table 3, the surface reactivity of cones180˚ is increased by increasing the cone size. Also, it is found that the

Table 1. The configuration structures and the dipole moments of the CNCs 60˚ and BNNCs 60˚. The dipole moment is given by Debye.

M1 refers to the type of connected atoms edges is nitrogen atoms, M2 refers to the type of connected atoms edges is boron atoms.

Table 2. The configuration structures and the dipole moments of the CNCs 120˚ and BNNCs 120˚. The dipole moment is given by Debye.

Table 3. The configuration structures and the dipole moments of the CNCs 180˚ and BNNCs 180˚. The dipole moment is given by Debye.

B68N73H21 (M1) cone 180˚ possesses the highest surface reactivity (15.9 Debye), followed by C141H21 cone 180˚ (11.5 Debye) and the smallest surface reactivity is for B73N68H21 (M2) cone 180˚ (11.2 Debye). In addition, the surface reactivity of cones 180˚ is found to be higher than the surface reactivity of the cones 120˚ and the latter is higher than the cones 60˚.

From Table 4, the calculated surface reactivity of CNCs and BNNCs at disclination angle 240˚ is increased by increasing the cone size. Also, it is found that the surface reactivity of the C94H14 cone 240˚ (15.3 Debye) is higher than the surface reactivity of B47N47H14 cone 240˚ (11.9 Debye). From Tables 1-4, it is clear that by increasing the disclination angle (from cones 60˚ to cones 240˚), the surface reactivity is increased.

The surface reactivity of CNCs and BNNCs at disclination angle 300˚ is shown in Table 5. It is found that the surface reactivity of cones 300˚ is increased by increasing the cone size. Also, it is found that the B41N49H10 (M1) cone 300˚ possesses the highest surface reactivity (23.5 Debye), followed by C90H10 cone 300˚ (23.5 Debye) and the smallest surface reactivity is for B49N41H10 (M2) cone 300˚ (12.5 Debye).

From Tables 1-5, one can report that the surface reactivity is increased by increasing the cone size and the disclination angle. In addition, the highest surface reactivity is for BnNmHy (M1) cones when the type of connected atoms edges is nitrogen atoms (at disclination angles 60˚, 180˚ and 300˚), otherwise the surface reactivity of C2mHy cones is always higher than BmNmHy cones (at disclination angles 120˚ and 240˚).

3.3. Surface Reactivity of Mono-Hydrogenated CNCs and BNNCs

The surface reactivity of mono-hydrogenated CNCs at three different hydrogenation sites HS1 (above the first neighbor atom of the apex atoms), HS2 (above one atom of the apex atoms) and HS3 (above one atom far from the apex atoms) for each disclination angels n × 60˚ (n = 1 to 5) are studied, see Figure 2. Also, at any disclination angle the effect of the cone size is tested through working on four different structures. By this way, one can test the effects of size, disclination angle and hydrogenation site of CNCs on the surface reactivity of mono- hydrogenated CNCs. From Table 6, it is found that the surface reactivity of mono-hydrogenated CNCs is increased by increasing the size and dsiclination angles of CNCs. Also, the surface reactivity at hydrogenation site HS3 always possesses the highest dipole moment, followed by HS1 and HS2 sites. The highest surface reactivities at disclination angles 60˚, 120˚, 180˚, 240˚ and 300˚ are found to be 10.94 Debye for C170H30-HS1, 12.58 Debye for C136H24-HS3, 18.70 Debye for C141H21-HS1, 16.98 Debye for C94H14-HS3 and 22.17 Debye for C90H10-HS3, respectively. Finally, on can conclude that the best cone size, the best hydrogenation site and best disclination angle for hydrogen adsorption is for CNC C90H10-HS3 at 300˚ declination angle.

Table 4. The configuration structures and the dipole moments of the CNCs 240˚ and BNNCs 240˚. The dipole moment is given by Debye.

Table 5. The configuration structures and the dipole moments of the CNCs 300˚ and BNNCs 300˚. The dipole moment is given by Debye.

Table 6. The configuration structures and the dipole moments of mono-hydrogenated CNCs. The dipole moment is given by Debye.

HS1 refers to the hydrogenation site is above the first neighbor atom of the apex atoms. HS2 refers to the hydrogenation site is above one atom of the apex atoms. HS3 refers to the hydrogenation site is above one atom far from the apex atoms.

Also, the surface reactivity of mono-hydrogenated BNNCs at three different hydrogenation sites HS1, HS2 and HS3 for each disclination angels n × 60˚ (n = 1 to 5), are studied. The hydrogen atom can be adsorbed on boron atom (named Type1) or can be adsorbed on nitrogen atom (named Type2). As we mention above, for disclination angles 60˚, 180˚ and 300˚ there are two models of BNNCs, BNNCs-M1 (the connected edges atoms are nitrogen atoms) and BNNCs-M2 (the connected edges atoms are boron atoms). Therefore, for the surface reactivity of BNNCs at disclination angles 120˚ and 240˚ for there are two systems of BNNCs, BNNCs-Type1 and BNNCs-Type2, see Table 6. However, for disclination angles 60˚, 180˚ and 300˚ there are four systems of BNNCs, BNNCs-M1-Type1, BNNCs-M1-Type2, BNNCs-M2-type1 and BNNCs-M2-Type2, see Table 7 and Table 8.

From Table 7, it is clear that the surface reactivities for monohydrogenated BNNCs-Type1 and BNNCs-Type2

Figure 2. Schematic representation for hydrogenation sites 1-HS1, 2-HS2, 3-HS3 of CNCs with disclination angle 120˚ for structures (a) C36H12, (b) C56H16, (c) C92H20 and (d) C136H24. Circles refer to the hydrogenation sites.

are increased by increasing the cone size and cone angle. For disclination angle 120˚, the highest surface reactivates for Type1 and Type2 are found to be 10.90 Debye for B68N68H24-HS2 and 11.75 Debye for B68N68H24-HS1, respectively. For disclination angle 240˚, the highest surface reactivates for Type1 and Type2 are found to be 12.85 Debye for B47N47H14-HS2 and 25.20 Debye for B34N34H12-HS3, respectively.

Table 8 shows that the surface reactivities for mono-hydrogenated BNNCs-M1-Type1 and BNNCs-M1- Type2 are increased by increasing the cone size and cone angle. For disclination angle 60˚, the highest surface reactivates for Type1 and Type2 are found to be 13.96 Debye for B38N42H20-HS2 and 10.77 Debye for B83N87H30-HS1, respectively. For disclination angle 180˚, the highest surface reactivates for Type1 and Type2 are found to be 20.81 Debye for B68N73H21-HS1 and 28.50 Debye for B35N40H15-HS3, respectively. For disclination angle 300˚, the highest surface reactivates for Type1 and Type2 are found to be 28.97 Debye for B41N49H10-HS1 and 27.55 Debye for B26N32H8-HS3, respectively.

From Table 9, it is found that the surface reactivities for mono-hydrogenated BNNCs-M2-Type1 and BNNCs-M2-Type2 are increased by increasing the cone size and cone angle. For disclination angle 60˚, the highest surface reactivates for Type1 and Type2 are found to be 7.44 Debye for B87N83H30-HS2 and 20.77 Debye for B42N38H20-HS1, respectively. For disclination angle 180˚, the highest surface reactivates for Type1 and Type2 are found to be 11.83 Debye and 17.28 Debye B73N68H21-HS3, respectively. For disclination angle 300˚, the highest surface reactivates for Type1 and Type2 are found to be 7.26 Debye and 25.94 Debye for B32N26H8-HS3, respectively. We can conclude that the best cone size, the best hydrogenation site and best disclination angle for mono-hydrogenated Type1 and Type2 of BNNCs is 28.97 Debye for B41N49H10-HS1 and 27.55 Debye for B26N32H8-HS3 at 300˚ declination angle.

Finally, it is found that the surface reactivity for pure and mono-hydrogenated CNCs and BNNCs is increased by increasing the cone angle and the cone size. Also, it is found that the surface reactivity is increased by hydrogenation and the highest surface reactivity is found to be 28.97 Debye for B41N49H10-HS1 (M1) when the hydrogenation is done on the boron atom (Type1).

Table 7. The configuration structures and the dipole moments of mono-hydrogenated BNNCs-Type1 and BNNCs-Type2 for disclination angles 120˚ and 240˚. The dipole moment is given by Debye.

Table 8. The configuration structures and the dipole moments of mono-hydrogenated BNNCs-M1-Type1 and BNNCs-M1- Type2, for disclination angles 60˚, 180˚ and 300˚. The dipole moment is given by Debye.

Table 9. The configuration structures and the dipole moments of mono-hydrogenated BNNCs-M2-Type1 and BNNCs-M2- Type2, for disclination angles 60˚, 180˚ and 300˚. The dipole moment is given by Debye.

4. Conclusion

The surface reactivity of fifty-two structures of pure CNCs and BNNCs and two hundred and fifty-two structures for mono-hydrogenated CNCs and BNNCs is calculated using density functional (DFT) calculations at the B3LYP/3-21G level of theory. Five disclination angles (60˚, 120˚, 180˚, 240˚ and 300˚), four different nanocone sizes and three different hydrogenation sites are applied. The calculations show that the dipole moments are always increased by increasing the nanocone sizes and the highest surface reactivity for pure CNCs and BNNCs at disclination angles 60˚, 180˚ and 300˚ is 23.50 Debye for B41N49H10 cone and at disclination angles 120˚ and 240˚ is 15.30 Debye for C94H14 cone. For mono-hydrogenated CNCs, the highest surface reactivity is found 22.17 Debye C90H10-HS3 at CNC angle 300˚ and for mono-hydrogenated BNNCs the highest surface reactivity is 28.97 Debye for B41N49H10-HS1 when the hydrogen atom is adsorbed on boron atom at BNNC angle 240˚.

Cite this paper

Ahlam A.El-Barbary,Mohamed A.Kamel,Khaled M.Eid,Hayam O.Taha,Rasha A.Mohamed,Mohammed A.Al-Khateeb,11, (2015) The Surface Reactivity of Pure and Monohydrogenated Nanocones Formed from Graphene Sheets. Graphene,04,75-83. doi: 10.4236/graphene.2015.44008

References

  1. 1. Liao, M.L. (2012) Preparation of Ni/Cu Composite Nanowires. Journal of Nanoparticle Research, 14, 837. http://dx.doi.org/10.1007/s11051-012-0837-1

  2. 2. Gotzias, A., Heiberg-Andersen, H., Kainourgiakis, Th. and Steriotis, M. (2011) Study of Hydrogen Adsorption in Nano-Structured Carbon Materials, with a Combination of Experimental Methods and Monte Carlo Simulations. Carbon, 49, 2715-2724. http://dx.doi.org/10.1016/j.carbon.2011.02.062

  3. 3. Mirzaei, M. (2013) Investigating Pristine and Carbon-Decorated Silicon Nanocones: DFT Studies. Superlattices and Microstructures, 58, 130-134. http://dx.doi.org/10.1016/j.spmi.2013.03.007

  4. 4. Esrafili, M.D. and Mahdavinia, G. (2013) Nitrogen-Doped (6,0) and (4,4) Single-Walled SiC Nanotubes: A DFT Study on Surface Reactivity and NMR Parameters. Superlattices and Microstructures, 62, 140-148. http://dx.doi.org/10.1016/j.spmi.2013.07.015

  5. 5. Yu, X. and Raaen, S. (2013) The Influence of Potassium Doping on Hydrogen Adsorption on Carbon Nanocone Material Studied by Thermal Desorption and Photoemission. Applied Surface Science, 270, 364-369. http://dx.doi.org/10.1016/j.apsusc.2013.01.031

  6. 6. Nouri, A. and Mirzaei, M. (2009) A Detailed Theoretical Study of the Interaction of Thiourea with Cis-Diaqua (Ethylenediamine) Platinum(II). Journal of Molecular Structure: THEOCHEM, 913, 207-209. http://dx.doi.org/10.1016/j.theochem.2009.07.044

  7. 7. Qu, C.Q., Qiao, L., Wang, C., Yu, S.S., Jiang, Q. and Zheng, W.T. (2010) Experimental and Theoretical Studies on the Magnetic Property of Carbon-Doped ZnO. Physics Letters A, 374, 782-787. http://dx.doi.org/10.1016/j.physleta.2009.11.066

  8. 8. Yu, X., Tverdal, M., Raaen, S., Helgesen, G. and Knudsen, K.D. (2008) Hydrogen Storage in Carbon Cones. Applied Surface Science, 255, 1906-1910. http://dx.doi.org/10.1016/j.apsusc.2008.06.120

  9. 9. Azevedo, S. (2004) Spin Polarization in Carbon Nanostructures with Disclinations. Physics Letters A, 325, 283-286. http://dx.doi.org/10.1016/j.physleta.2004.03.065

  10. 10. Baei, M.T., Peyghan, A.A., Bagheri, Z. and Tabar, M.B. (2012) DFT Study of NH3 Adsorption on Pristine, Ni- and Si- Doped Graphynes. Physics Letters A, 377, 107-111. http://dx.doi.org/10.1016/j.physleta.2012.11.006

  11. 11. Mirzaei, M. and Meskinfam, M. (2011) Computational Studies of Effects of Tubular Lengths on the NMR Properties of Pristine and Carbon Decorated Boron Phosphide Nanotubes. Solid State Sciences, 13, 1926-1930. http://dx.doi.org/10.1016/j.solidstatesciences.2011.08.018

  12. 12. Ge, M. and Sattler, K. (1994) Bundles of Carbon Nanotubes Generated by Vapor-Phase Growth. Applied Physics Letters, 64, 710. http://dx.doi.org/10.1063/1.111042

  13. 13. Ge, M. and Sattler, K. (1994) Observation of Fullerene Cones. Chemical Physics Letters, 220, 192-196. http://dx.doi.org/10.1016/0009-2614(94)00167-7

  14. 14. Merkulov, V.I., Melechko, A.V., Guillom, M.A., Lowndes, D.H. and Simpson, M.L. (2011) Sharpening of Carbon Nanocone Tips during Plasma-Enhanced Chemical Vapor Growth. Chemical Physics Letters, 350, 381-385. http://dx.doi.org/10.1016/S0009-2614(01)01312-4

  15. 15. Shenderova, O.A., Lawson, B.L., Areshkin, D. and Brenner, D.W. (2001) Predicted Structure and Electronic Properties of Individual Carbon Nanocones and Nanostructures Assembled from Nanocones. Nanotechnology, 12, 191. http://dx.doi.org/10.1088/0957-4484/12/3/302

  16. 16. Baylor, L.R., Merkulov, V.I., Ellis, E.D., Guillorn, M.A., Lowndes, D.H., Melechko, A.V., Simpson, M.L. and Whealton, J.H. (2002) Field Emission from Isolated Individual Vertically Aligned Carbon Nanocones. Journal of Applied Physics, 91, 4602-4606. http://dx.doi.org/10.1063/1.1455138

  17. 17. Lu, X., Yang, Q., Xiao, C. and Hirose, A. (2006) Field Electron Emission of Carbon-Based Nanocone Films. Applied Physics A, 82, 293-296. http://dx.doi.org/10.1007/s00339-005-3410-2

  18. 18. Bonard, J.M., Gaal, R., Garaj, S., Thien-Nga, L., Forro, L., Takahashi, K., Kokai, F., Yudasaka, M. and Iijima, S. (2002) Field Emission Properties of Carbon Nanohorn Films. Journal of Applied Physics, 91, Article ID: 10107. http://dx.doi.org/10.1063/1.1481200

  19. 19. Saito, Y., Tsujimoto, Y., Koshio, A. and Kokai, F. (2007) Field Emission Patterns from Multiwall Carbon Nanotubes with a Cone-Shaped Tip. Applied Physics Letters, 90, Article ID: 213108. http://dx.doi.org/10.1063/1.2742637

  20. 20. McGuire, K., Gothard, N., Gai, P.L., Dresselhaus, M.S., Sumanasekera, G. and Rao, A.M. (2005) Synthesis and Raman Characterization of Boron-Doped Single-Walled Carbon Nanotubes. Carbon, 43, 219-227. http://dx.doi.org/10.1016/j.carbon.2004.11.001

  21. 21. Cruz-Silva, E., Cullen, D.A., Gu, L., Romo-Herrera, J.M., Munoz-Sandoval, E., Lopez-Urias, F., Sumpter, B.G., Meunier, V., Charlier, J.C., Smith, D.J., Terrrones, H. and Terrones, M. (2008) Heterodoped Nanotubes: Theory, Synthesis, and Characterization of Phosphorus-Nitrogen Doped Multiwalled Carbon Nanotubes. ACS Nano, 2, 441-448. http://dx.doi.org/10.1021/nn700330w

  22. 22. Terrones, M., Souza, A.G. and Rao, A.M. (2008) Nanostructured Photoelectrodes for Solar Powered Applications. In: Jorio, A., Dresselhaus, G. and Dresselhaus, M.S., Eds., Carbon Nanotubes, Topics in Applied Physics, Vol. 111, Springer-Verlag, Berlin, 531-566. http://dx.doi.org/10.1007/978-3-540-72865-8_17

  23. 23. Becke, A.D. (1993) Density Functional Thermochemistry. III. The Role of Exact Exchange. The Journal of Chemical Physics, 98, 5648. http://dx.doi.org/10.1063/1.464913

  24. 24. Becke, A.D. (1988) Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Physical Review A, 38, 3098. http://dx.doi.org/10.1103/PhysRevA.38.3098

  25. 25. Frisch, M.J., Trucks, G.W., Schlegel, H.B., Scuseria, G.E., Robb, M.A., Cheeseman, J.R., Zakrzewski, V.G., Montgomery, J.A., Stratmann, R.E., Burant, J.C., Dapprich, S., Millam, J.M., Daniels, A.D., Kudin, K.N., Strain, M.C., Farkas, O., Tomasi, J., Barone, V., Cossi, M., Cammi, R., Mennucci, B., Pomelli, C., Adamo, C., Clifford, S., Ochterski, J., Petersson, G.A., Ayala, P.Y., Cui, Q., Morokuma, K., Malick, D.K., Rabuck, A.D., Raghavachari, K., Foresman, J.B., Cioslowski, J., Ortiz, J.V., Stefanov, B.B., Liu, G., Liashenko, A., Piskorz, P., Komaromi, I., Gomperts, R., Martin, R.L., Fox, D.J., Keith, T., Al-Lamham, M.A., Peng, C.Y., Nanayakkara, A., Gonzalez, C., Challacombe, M., Gill, P.M.W., Johnson, B.G., Chen, W., Wong, M.W., Andres, J.L., Head-Gordon, M., Replogle, E.S. and Pople, J.A. (2004) Gaussian 2004. Gaussian Inc., Wallingford.

  26. 26. Hindi, A. and EL-Barbary, A.A. (2015) Hydrogen Binding Energy of Halogenated C40 Cage: An Intermediate between Physisorption and Chemisorption. Journal of Molecular Structure, 1080, 169-175. http://dx.doi.org/10.1016/j.molstruc.2014.09.034

  27. 27. El-Nahass, M.M., Kamel, M.A., El-Barbary, A.A., El-Mansy, M.A.M. and Ibrahim, M. (2013) FT-IR Spectroscopic Analyses of 3-Methyl-5-Pyrazolone (MP). Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 111, 37-41. http://dx.doi.org/10.1016/j.saa.2013.03.072

  28. 28. Kotz, J.C., Treichel, P. and Weaver, G.C. (2006) Chemistry and Chemical Reactivity. Thomson Brooks/Cole, Pacific Grove.

NOTES

*Corresponding author.