Theoretical study of the interaction between molecular hydrogen and [MC60]+ complexes

Maitreyi Robledo*a, Sergio Díaz-Tenderoac, Fernando Martínabc and Manuel Alcamíab
aDepartamento de Química, Módulo 13, Universidad Autónoma de Madrid, 28049 Madrid, Spain. E-mail: Maitreyi.robledo@uam.es
bInstituto Madrileño de Estudios Avanzados en Nanociencias (IMDEA_Nanociencia), Cantoblanco, 28049, Spain
cCondensed Matter Physics Center (IFIMAC), Universidad Autónoma de Madrid, 28049 Madrid, Spain

Received 7th January 2016 , Accepted 8th March 2016

First published on 10th March 2016


Abstract

In this work we present a density functional theory study of the interaction between a positively charged exohedral metallofullerene and several hydrogen molecules. For this purpose we have chosen Li+, Ti+, V+ and Cu+ as the metal cations, since they represent a good sample of ionic or covalent interactions with the carbon cage. We have found that the interaction between the hydrogen molecules and the metal cation strongly depends on the type of interaction between the metal cation and the fullerene. Furthermore, the hydrogen saturation around the metal is also determined by the nature of the metal.


Introduction

Hydrogen is both by mass and number the most abundant element in the known universe, and has become in the last years a candidate to replace fossil fuels. Hydrogen presents significant advantages over petroleum-derivative fuels, the most important being that hydrogen does not contribute to global warming, does not produce toxic pollutants, and has a very high energy content. One of the greatest challenges is to find out materials that can effectively store this promising non-contaminant fuel in a reversible way. According to the US Department of Environment (DOE) the target is to obtain a hydrogen storage capacity of 6.5 wt%. Many storage methods have been already proposed, for instance liquid hydrogen,1 cycloalkanes as chemical hydrides,2 Metal-Organic Frameworks (MOF)3,4 or carbon nanotubes.5–9 Though carbon nanotubes were at first the best option, due to their high surface-to-volume ratio, other studies10–12 have shown that the adsorption potential for hydrogen on this kind of materials is less than for an array of idealized slit pores geometries.

On the other hand, following this line of carbon structures, fullerenes have also been studied as promising hydrogen storage devices. A first approach mainly considered their capacity to store H2 inside the cage. For instance, C110 fullerene was found to encapsulate up to 19 H2 molecules;13 these hydrogen atoms cannot be attached to the inner wall of the structure and therefore they exist in their molecular form. Other studies14 have investigated, by molecular dynamics simulations, the possibility of H2 storage in carbon clusters (as well as in boron nitride clusters) finding out that H2 molecules are able to abide stable in a C60 fullerene at 298 K and 0.1 MPa, however it has been confirmed that pressures up to 5 MPa are needed to store H2 molecules in the C60 cage.

Recent studies5,15 propose to decorate carbon structures with metals, for instance titanium-decorated fullerenes5 or carbon nanotubes,16 as efficient materials to store H2. It was shown how a titanium atom coated on a single-walled nanotube or to a fullerene binds up to four hydrogen molecules, thus offering an efficient storage material for practical applications. Y. Zhao et al.16 showed by means of Density Functional Theory calculations that C60Sc12 could bind up to 5 hydrogen molecules per transition metal; with binding energy between 0.26 and 0.45 eV/H2 and storage density of 7 wt%. The main inconvenience using transition metals is that the former tend to cluster on the fullerene outer shell. Thus, an alkaline-earth metal as calcium,15,17 which binds directly to the hexagonal and pentagonal faces, without forming droplets, would be more suitable for the final goal of hydrogen storage. C60Ca32 can adsorb up to 62 H2 in two different layers. The first of them with an average binding energy of 0.45 eV/H and the second one, 0.11 eV/H2, which makes Ca coated C60 a practical hydrogen storage material at near ambient temperatures. Other promising elements to store hydrogen are alkaline metals such as lithium,18 since they enhance the H2 adsorption capacity of fullerenes; this effect is more pronounced for Na, with a hydrogen adsorption capacity of 9.5 wt%,19 than for Li. Other experimental studies20 have shown how a Li doped fullerene can reversibly desorb up to 5 wt% (at 270 °C), which is significantly less than the desorption temperature of an hydrogenated fullerene (C60Hx) (500–600 °C).

In this work we present a thorough study of the binding energies and geometries of a series of complexes of the form [C60M(H2)n]+, where M is Li, Ti, V or Cu and n = 1–5 hydrogen molecules. The choice of charged clusters is relevant since recent mass spectrometry experiments have allowed characterization of the sequential absorption of H2 molecules on C60 fullerenes21 as well as other molecules.22 Furthermore, Rao and Jena have shown23–25 that positively charged atoms would be able to bind a large amount of H2 molecules due to the charge polarization effect, which is in fact the main interaction mechanism we found between a metal cation and the fullerene cage in a previous work.26 All [M–C60]+ interactions are mainly governed by an ion-induced dipole interaction due to the positively charged metal and the large polarizability of the C60 cage. Concerning Li+ and Cu+ ions, our conclusion was that they present a larger ionic and covalent character, respectively; while for Ti+ and V+, the interaction between the metal and the carbon cage presented a shared character between ionic and covalent. Thus, in contrast with previous works we present four different metals (i.e. Li+, Cu+, Ti+ and V+), which interact distinctly with the C60 carbon cage (i.e. covalently, ionically or with a mixed character) and that present different abilities to storage H2 molecules, according to their type of interaction. We thus, provide a simple guide to predict the behavior of other metals coating the fullerene for their ability to storage H2.

Computational details

We have performed self-consistent field calculations for all the geometry optimizations in the framework of Density Functional Theory (DFT). In particular we have used the hybrid B3LYP functional, which combines the B3 non-local hybrid exchange potential proposed by Becke27 with the non-local gradient corrected correlation functional of Lee, Yang and Parr28 and the M06-2X hybrid meta-exchange correlation functional,29 which includes dispersion forces, that are in principle expected to play some role in the interaction of H2 with the fullerene cage and with neighboring H2 molecules. In both cases the 6-31G* Pople basis set has been used. In the case of the M06-2X functional, further optimizations have been performed using the 6-311++G(d,p) basis. Previous studies have compared the performance of the M06-2X functional with wavefunction methods, as the perturbative MP2 one, showing not only comparable but identical results.30 In addition, we have performed a zero point vibrational energy correction (ZPE) for all calculations. The optimization of the fullerene cage leads to a length of 1.39 Å for the bond between two hexagons (hh bond) and 1.45 Å for the one between a pentagon and a hexagon (hp bond). These results are in good agreement with previous theoretical works, e.g. Rhh = 1.38 Å, Rhp = 1.43 Å,31 Rhh = 1.40 Å, Rhp = 1.46 Å (ref. 32) and with electron diffraction measurements, Rhh = 1.40 Å, Rhp = 1.45 Å.33 Concerning the basis set, the 6-31G* in combination with the B3LYP functional has shown to give accurate results in comparison with those obtained with perturbational theory (MP2) in fullerenes.34 Comparison between 6-31G* and 6-311++G(d,p) allows to check the importance of including diffuse functions for a correct description of the weak H2⋯M+ interaction.

Calculations were carried out with the Gaussian 09 program package.35 In order to determine the most stable isomer, we considered in our previous work26 5 types of initial geometries with the metal atom interacting with the fullerene external surface: the hexagonal (h) and pentagonal (p) rings, the double (hh) and single (hp) bonds and on top of a carbon atom (C). Furthermore, several spin multiplicities were also computed for the ground state. According to these results, we present here only the interaction between the most stable isomer for each [MC60]+ complex and the H2 molecules. This is for [C60Li]+ the h isomer with singlet spin multiplicity, for [C60Ti]+ the h isomer with quadruplet spin multiplicity, for [C60V]+ the h isomer with quintuplet spin multiplicity and the hh isomer with singlet spin multiplicity for [C60Cu]+.

The hydrogenation energy has been determined as follows:

EHyd = E[C60M(H2)n+] − (E[C60M(H2)n−1+] + E[H2])
where E[C60M(H2)n+] is the total energy of the hydrogenated complex, E[C60M(H2)n−1+] is the energy of the hydrogenated complex with one less H2 molecule and E[H2] is the energy of one single hydrogen molecule. The hydrogenated energy is thus the energy needed to add one more hydrogen molecule to the previous hydrated complex.

We have as well calculated the adsorption energy per H2, defined as

image file: c6ra00501b-t1.tif
where E[C60M(H2)n+] is the total energy of the hydrogenated complex, E[C60M+] is the energy of the metallofullerene complex and E[H2] is the energy of a single hydrogen molecule. Thus, negative values of the EAds will show stable adsorbed complexes and therefore the hydrogenation process would be an exothermic course.

Results

Taking the most stable isomer of each [C60M]+ species (M = Li, Ti, V, Cu),26 hydrogen molecules have been then placed on a top configuration, i.e. forming an angle of 180° between the center of the fullerene, the metal and the center of the H2 molecule, since it is well known from previous works15 that this is the most favorable arrangement if one wants to maintain the hydrogen atoms in their molecular form. Thus, one single H2 molecule is initially placed on top of a metal cation (Li+, Ti+, V+ and Cu+) with H–H distance set to 0.743 Å, which was found to be the equilibrium distance for an isolated H2 molecule after optimization. In order to study the interaction with more than one H2 molecule we have placed sequentially 2, 3, 4 and 5 H2 molecules following the most symmetric configuration with respect to the metal atoms. Starting from these initial configurations the geometry of the complexes were fully optimized. Fig. 1 shows the optimized structures for [C60Li(H2)n]+ and [C60Cu(H2)n]+.
image file: c6ra00501b-f1.tif
Fig. 1 Optimized structures for [C60Li(H2)n]+ and [C60Cu(H2)n]+, n = 1–5, with M06-2X/6-311++G(d,p).

The result shows that the four systems present a certain limit number of hydrogen molecules that can be attached to the metal, i.e. there is a saturation limit. This limit is four hydrogen molecules for [C60Cu(H2)n]+, three H2 molecules for the [C60Li(H2)n]+ and [C60Ti(H2)n]+ systems and two H2 molecules for the [C60V(H2)n]+ (see Fig. 2). This behavior is observed at the three different levels of theory considered in this work. In fact the differences between the H2-cage distances, when H2 molecules are adsorbed to the metallofullerene system (i.e. when the saturation limit is not reached yet) is in the order of 0.5–2%. Although the three theory levels show similar trends, differences increase significantly (up to ∼30%) for those distances corresponding to H2 molecules weakly bound by the system. On the other hand, the H–H bond distances of H2 molecules in [C60Li]+ are 0.74 Å until saturation is reached (i.e. up to three hydrogen molecules), then it decreases (0.73 Å for [C60Li(H2)4]+ and [C60Li(H2)5]+). For [C60Ti]+, it changes from 0.76 Å to 0.73 Å in [C60Ti(H2)4]+ and [C60Ti(H2)5]+. For [C60V]+ it decreases from 0.76 Å when the complex is not saturated to 0.73 Å from the third H2 molecule onwards. And finally, for [C60Cu]+ it changes from 0.76 Å up to saturation of the complex to 0.75 Å in [C60Cu(H2)3]+ and 0.73 Å in [C60Cu(H2)4]+ and [C60Cu(H2)5. Hence, the H–H distance increases when the H2 molecule is interacting with the metal cation.


image file: c6ra00501b-f2.tif
Fig. 2 Optimized M–H2 distances (M = Li, Ti, V, Cu) at the three levels of theory.

This is a consequence of the H–H bond disruption by the metal; similar results were obtained by previous works5 on C60Ti12 where it was described that each Ti atom was able to bind 4 H2 molecules, but the first one was attached as a single H atom to Ti. Thus, the metal cation prefers to attach one single H atom, and form the hydride.

Fig. 3 shows the hydrogenation energy and the total adsorption energy per H2, as defined before. Among the four studied metals, Li shows the lowest adsorption energies and its saturation limit is reached almost at 3 H2 molecules. The interaction is caused by a polarization of the hydrogen molecule by the Li+ cation, which results in an electrostatic interaction between the metallofullerene and the H2 molecule.


image file: c6ra00501b-f3.tif
Fig. 3 Hydrogenation energy, EHyd., and adsorption energy, EAds., in kcal mol−1, for [C60Li(H2)n]+, where M = Li, Ti, V, Cu and n = 1–5.

Ti hydrated complexes also show a saturation limit of 3 H2 molecules, though adsorption energies are larger (∼5 kcal mol−1) compared to [C60Li(H2)n]+. [C60V(H2)n]+ shows a saturation limit of two hydrogen molecules, and the third H2 molecule releases ∼7 kcal mol−1 to be adsorbed. Among all the studied metals, Cu has the lowest H2 adsorption capability; [C60Cu(H2)]+ releases almost 10 kcal mol−1 to attached a second H2 molecule. All the studied complexes show an electrostatic interaction due to polarization. However, the adsorption energies are significantly larger for Ti, V and Cu, than for Li, especially for the adsorption of the first H2 molecule (5.56 kcal for TiC60+, 8.77 kcal for VC60+ and 10.10 kcal for CuC60+). Furthermore, if we take into account our previous conclusions on the type of interaction between these cations and the fullerene cage,26 we can establish a relation between the pure interaction M+–C60 and the adsorption energies H2–MC60+: Cu interacts with the carbon cage with larger covalent character and therefore presents a stronger M–C binding. This feature involves stable adsorption of the hydrogen molecules as well as high adsorption energy. On the other hand, Ti+ and V+ present a shared character, in between ionic and covalent, in their interaction with the fullerene. The double character of the M+–C60 binding will enhance this polarization and consequently the adsorption energies will be higher when more than one H2 molecule is adsorbed by the metallofullerene. The enhancement of the adsorption energy is due to the polarization that the positive metal exerts on the H2 molecule.

In order to give a deeper insight on the interaction, we have also evaluated the charge, in the framework of the Mulliken population analysis, on the metal cation after the interaction, as well as the charge on each H2 molecule after adsorption to the hydrated complex (see Fig. 4). The metal cation is initially charged with one single positive charge. The interaction between the metal cation and the fullerene cage involves the polarization of the latter by the former and, thus a gain of electron density (reducing the charge on the metal). This leads to a very localized reactivity, which results in the involvement of the closest carbon atoms in the interaction with the metal cation.


image file: c6ra00501b-f4.tif
Fig. 4 Charge localized on the metal cation and the surrounding hydrogen molecules, with the Mulliken approach.

When hydrogen molecules are added to the MC60+ complex, positive charge decreases on the metal and increases in some of the hydrogen molecules. This loss of neutrality by the H2 molecules, once attached to the complex, is the signature of the induced polarizability and thus, of induced dipole type of interaction between the metal and the H2 molecule. The charge in the metal is reduced and also shows a limit similar to the saturation limit.

Conclusions

In this work we have presented a DFT study on the interaction between cationic metallofullerenes and several hydrogen molecules (from one up to five). For this purpose we have chosen to study Li+, Ti+, V+ and Cu+. Each of these metal cations presents a different kind of interaction with the carbon cage (i.e. ionic or covalent-like). Despite the different interactions, all metallofullerenes exhibit a rather similar saturation limit for the adsorption of H2 molecules, which is around three H2 molecules. However, adsorption energies differ significantly from one system to the other (for instance, [CuC60]+ shows the largest ones) due to the different nature of the interaction between the metal atoms and the fullerene cage. These differences are also important for the reversibility of the H2 adsorption. We can thus conclude that in the case of M+…C60 interactions with ionic/dipolar character a larger charge remains in the metal atom, giving rise to a stronger adsorption of the H2 molecules. Where as for M+…C60 interactions with covalent character, the opposite behavior is observed.

This study provides important new insight into the adsorption of hydrogen molecules in exohedral metalofullerenes cationic complexes, showing the importance of the charge and the correct choice of a metal for hydrogen storage at the nanoscale.

Acknowledgements

We acknowledge the generous allocation of computer time at the Centro de Computación Científica at the Universidad Autónoma de Madrid (CCC-UAM). This work was partially supported by the projects CTQ2013-43698 P and FIS2013-42002 R (MINECO), and NANOFRONTMAG (CAM), and the European COST Action XLIC. M. R. acknowledges the FPI grant associated with the project CTQ2010-17006 of the Spanish Ministerio de Economía y Competitividad. S. D.-T. gratefully acknowledges the “Ramón y Cajal” program of the Spanish Ministerio de Educación y Ciencia.

References

  1. S. M. Aceves, J. Martinez and O. Gracía-Villazana, Int. J. Hydrogen Energy, 2000, 25, 1075 CrossRef CAS.
  2. R. B. Biniwale, S. Rayalu, S. Devotta and M. Ichikawa, Int. J. Hydrogen Energy, 2008, 33, 360 CrossRef CAS.
  3. H. Furukawa, K. E. Cordova, M. O'Keeffe and O. M. Yaghi, Science, 2013, 341, 974 CrossRef CAS PubMed.
  4. J.-R. Li, R. J. Kuppler and H.-C. Zhou, Chem. Soc. Rev., 2009, 38, 1477–1504 RSC.
  5. T. Yildrim and S. Ciraci, Phys. Rev. Lett., 2005, 94, 175501 CrossRef PubMed.
  6. H. M. Cheng, Q. H. Yang and C. Liu, Carbon, 2002, 40, 81–85 CrossRef.
  7. F. L. Darkrim, P. Malbrunot and G. P. Tartaglia, Int. J. Hydrogen Energy, 2002, 27, 193–202 CrossRef CAS.
  8. A. C. Dillon, K. M. Jones, T. A. Bekkedahl, C. H. Kiang, D. S. Bethune and M. J. Heben, Nature, 1997, 386, 377 CrossRef CAS.
  9. C. Liu, Y. Y. Fan, M. Liu, H. T. Cong, H. M. Cheng and M. S. Dresselhaus, Science, 1999, 286, 1127 CrossRef CAS PubMed.
  10. Q. Wang and J. Karl Johnson, J. Chem. Phys., 1999, 110, 577 CrossRef CAS.
  11. S. Jalili and R. Majidi, Phys. E, 2007, 39(1), 166–170 CrossRef CAS.
  12. S. Jalili, A. Jaberi, M. G. Mahjani and M. Jafarian, Int. J. Nanosci., 2009, 8, 425 CrossRef CAS.
  13. E. Barajas-Barraza and R. A. Guirado, Diamond Relat. Mater., 2002, 11, 945 CrossRef.
  14. I. Narita and T. Oku, Diamond Relat. Mater., 2002, 11, 945 CrossRef CAS.
  15. Q. Yang, Q. Sung, P. Jena and Y. Kawazoe, J. Chem. Theory Comput., 2009, 5, 374 CrossRef PubMed.
  16. Y. Zhao, Y.-H. Kim, A. C. Dillon, M. J. Heben and S. B. Zhang, Phys. Rev. Lett., 2005, 94, 155504 CrossRef PubMed.
  17. M. Yoon, S. Yang, C. Hicke, E. Wang, D. Geohegan and Z. Zhang, Phys. Rev. Lett., 2008, 100, 206806 CrossRef PubMed.
  18. Q. Sung, P. Jena, Q. Wang and M. Marquez, J. Am. Chem. Soc., 2006, 128, 9741 CrossRef PubMed.
  19. K. R. S. Chandrakumar and S. K. Ghosh, Nano Lett., 2008, 8, 13–19 CrossRef CAS PubMed.
  20. J. A. Teprovich, M. S. Wellons, R. Lascola, S. J. Hwang, P. A. Ward, R. N. Compton and R. Zidan, Nano Lett., 2012, 12, 582–589 CrossRef CAS PubMed.
  21. A. Kaiser, C. Leidlmair, P. Bartl, S. Zottl, S. Denifl, A. Mauracher, M. Probst, P. Scheier and O. Echt, J. Chem. Phys., 2013, 138, 074311 CrossRef CAS PubMed.
  22. O. Echt, A. Kaiser, S. Zottl, A. Mauracher, S. Denifl and P. Scheier, ChemPlusChem, 2013, 78, 910–920 CrossRef CAS.
  23. B. K. Rao and P. Jena, Europhys. Lett., 1992, 20, 307 CrossRef CAS.
  24. J. Niu, B. K. Rao, P. Jena and M. Manninen, Phys. Rev. B: Condens. Matter Mater. Phys., 1995, 51, 4475 CrossRef CAS.
  25. J. Niu, B. K. Rao and P. Jena, Phys. Rev. Lett., 1992, 68, 2277 CrossRef CAS PubMed.
  26. M. Robledo, N. Aguirre, S. Díaz-Tendero, F. Martín and M. Alcamí, RSC Adv., 2014, 4, 53010 RSC.
  27. A. D. Becke, Phys. Rev. A: At., Mol., Opt. Phys., 1988, 38, 3098 CrossRef CAS.
  28. C. Lee, W. Yang and R. G. Parr, Phys. Rev. B: Condens. Matter Mater. Phys., 1988, 33, 785 CrossRef.
  29. Y. Zahao and D. G. Thrular, Theor. Chem. Acc., 2008, 120, 215 CrossRef.
  30. M. Walker, A. J. A. Harvey, A. Sen and C. E. H. Dessent, J. Phys. Chem. A, 2013, 117(47), 12590–12600 CrossRef CAS PubMed.
  31. N. Hamamoto, J. Jitsukawa and C. Satoko, Eur. Phys. J. D, 2002, 19, 211 CAS.
  32. F. J. Rabilloud, J. Phys. Chem. A, 2010, 114, 7241–7247 CrossRef CAS PubMed.
  33. K. Hedberg, L. Hedberg, D. S. Bethune, C. A. Brown, H. D. Dorn, R. D. Johnson and M. de Vries, Science, 1991, 254, 410 CAS.
  34. M. Haser, J. Almof and G. E. Scuseria, Chem. Phys. Lett., 1991, 181(6), 497 CrossRef CAS.
  35. M. J. Frisch, et al., Gaussian 09 Revised B01, Gaussian Inc, Wallington CT, 2010 Search PubMed.

This journal is © The Royal Society of Chemistry 2016
Click here to see how this site uses Cookies. View our privacy policy here.