Titanium-capped carbon chains as promising new hydrogen storage media
Chun-Sheng Liu, Hui An, and Zhi Zeng
Received Xth XXXXXXXXXX 20XX,
Accepted Xth XXXXXXXXX 20XX
First published on the web Xth XXXXXXXXXX 200X
The capacity of Ti-capped sp carbon atomic chains for use as hydrogen storage media is studied using first-principles density functional theory. The Ti atom is strongly attached at one end of the carbon chains via d-p hybridization, forming stable TiC complexes. We demonstrate that the number of adsorbed H on Ti through Kubas interaction depends upon the chain types. For polyyne (n even) or cumulene (n odd) structures, each Ti atom can hold up to five or six H molecules, respectively. Furthermore, the TiC chain effectively terminated on a C fullerene can store hydrogen with optimal binding of 0.52 eV/H. Our results reveal a possible way to explore high-capacity hydrogen storage materials in truly one-dimensional carbon structures.
1 Introduction††footnotetext: AKey Laboratory of Materials Physics, Institute of Solid State Physics, Chinese Academy of Sciences, Hefei 230031, People’s Republic of China. Fax: +86-551-5591434; Tel: +86-551-5591407; E-mail: firstname.lastname@example.org††footnotetext: Graduate School of the Chinese Academy of Sciences, Beijing 100094, People’s Republic of China.
To develop economical hydrogen energy, carbon nanostructures with sp-like bonding functionalized by transition metal (TM),1, 2, 3, 4, 5, 6 alkali metal (AM) 7, 8 and alkaline-earth metal (AEM) atoms 9 have been expected to be promising storage materials due to their light weights and large surface areas. Recently, Iijima et al. 10 transformed graphene to single sp-hybridized carbon chains containing 16 carbon atoms. In addition, chains with metal atoms connected to the ends have been previously generated,11 and their magnetic, electronic, and transport properties have been studied extensively.12, 13 These results inspired us to consider whether metal-capped carbon atomic chains are efficient hydrogen storage media. In this study, we conduct theoretical studies of high-capacity hydrogen storage media consisting of a Ti atom capped on two kinds of atomic carbon chains, cumulene (with double C-C bonds) or polyyne (with alternating singlet and triplet bonds). The number of adsorbed H molecules is only determined by the type of chain. Each Ti atom in TiC and TiC can bind up to six and five H molecules, respectively, corresponding to storage capacities of 10 wt % and 6.5 wt %. The average binding energy of H molecules on TiC (TiC) is 0.59 eV/H (0.57 eV/H), which is between the physisorption and chemisorption energies.
2 Computational method
Numerical calculations were performed using spin-polarized density functional theory (DFT) with the Perdew-Wang (1991) exchange-correlation function,14 as implemented in the DMol package (Accelrys Inc.).15 A double numerical atomic orbital augmented by d-polarization functions (DNP) was employed as the basis set. In the self-consistent field calculations, the electronic-density convergence threshold was set to 110 electron/Å. Geometric optimization was performed with convergence thresholds of 10 Ha for the energy, 210 Ha/Å for the force, and 10 Å for the atomic displacements. We performed normal-mode analysis of the obtained structures to ensure that the structures optimized without any symmetry constraints were true minima of the potential-energy surface.
3 Results and discussion
We first consider the bonding energetics of a single Ti atom at one end of the C (n=5-16) chains. The binding energies of Ti, shown in Fig. 1(a), are much larger than those of Ti on sp- hybridized carbon materials,2, 3, 4, 5 and the adsorption strengths are thus strong enough to maintain stable Ti-C structures at room temperature. Moreover, the Ti binding energies exhibit an interesting odd-even oscillatory behavior, consistent with the trend of Mulliken charge [see Fig. 1(a)]. This behavior suggests that the bonding between the Ti atom and carbon chains has a significant ionic characteristic. In addition, the configurations of TiC complexes exhibit significant differences for odd or even values of n. For instance, the C-C distances in TiC, depicted in Fig .1(b), are rather uniform. In contrast, the bond lengths of TiC alternate, consistent with - bonding [Fig. 1(c)]. Contour plots of the total electronic charge density confirm these bonding patterns.
To elucidate the nature of the bonding, the projected density of states (DOS) shown in Fig. 1(d) illustrates the bonding orbitals of the C (C) and Ti atom result from the hybridization of the carbon chain orbitals and the Ti-d orbitals. The C-p and Ti-d hybridization is evidently stronger in TiC than in TiC, in accordance with the larger binding energy of TiC. Since odd-numbered and even-numbered free chains have full and half-occupied HOMOs (highest occupied molecular orbitals), respectively, different states are responsible for the bonding of the Ti atom to C or C. From the isosurface plots of these states [Fig. 1(d)], it is clearly that the LUMO (a linear combination of degenerate and antibonding orbitals) of C hybridizes with two Ti-d orbitals to create a new -symmetry bond. However, for TiC there is an overlap between one Ti-d orbital and the antibonding state of C. As mentioned above, the bonding mechanisms for TiC and TiC seem to be very similar to the Dewar coordination.16
The stability of the TiC and TiC complexes was further tested by normal mode analysis, and no negative mode was found. However, it remains to be shown whether it is possible to realize linear C-Ti chains experimentally. To address this issue, we calculated the minimum-energy path (MEP) for different reaction processes with the nudged elastic band method to determine the energy barrier.17 The image number considered is 16 to ensure that the obtained MEP is correct. In Figs. 2(a) and 2(b), the MEPs show that a single Ti atom is attached to the right end of C and C without any energy barrier. The other end is free and ready to be attached to any other structure with a high surface-to-volume ratio. The main Ti vibrational frequencies, 623 cm in TiC and 359 cm in TiC, correspond to the stretching modes of C-Ti. These characteristic modes will provide a reference for the Raman/IR spectra of the synthesized materials. To this end, a monatomic carbon chain, realized by removing the carbon atoms from the graphene through energetic electron irradiation inside a transmission electron microscope,10 can be used as the initial material in the fabrication process.
Now we investigate the interaction between these complexes and hydrogen molecules. The first H molecule is found to dissociate on both TiC and TiC complexes and bind atomically to Ti, with interatomic H distances of 3.003 Å [Fig. 3(a)] and 2.748 Å [Fig. 3(b)], respectively. Such a dissociation is observed in most TM-hydrogen binding [Ref. 2] due to excessive charge transfer from the TM to the antibonding state of H. Additional H molecules, however, bind molecularly around the Ti atom, since the charge transfer per H molecule is not enough to destabilize the dihydrogen state when more hydrogen molecules are added to the system.
For TiC, several initial configurations of the second H molecule were considered in the search for the lowest energy structure. The structure in which one of the H molecules binds as a dihydride is 0.09 eV higher in energy than the structure in which both H molecules bind in molecular form [Fig. 3(a)]. We continue to introduce successive hydrogen molecules near the Ti atom. Surprisingly, it is energetically favorable for the TiC complex to adsorb 6 H molecules, corresponding to a 10 wt % gravimetric density. A close examination of the geometry of the TiC(H) configuration shows that the adsorbed hydrogen molecules are divided into two sets [Fig. 2(a)]. First, the molecule on the top binding site, which is vertical to the carbon chain axis, has a bond length of 0.808 Å and the largest bonding energy, 0.82 eV. Second, the H-H bonds of the side H molecules lie almost parallel to the chain axis, with the lower H atoms of each H molecule tending to tilt toward the chain. The bond lengths between the Ti and the lower (upper) H atoms are 1.910 Å (1.970 Å). The side molecules are more weakly bound, with a binding energy of 0.54 eV/H and an elongated H-H bond distance of 0.810 Å, on average. All of these features suggest that different bonding mechanisms (to be discussed below) for the two sets of H molecules on TiC.
On the other hand, when the second H molecule is bound to TiC, the structural configuration shown in Fig. 3(b) is more energetically favorable than its isomer (not shown) by 1.02 eV. As we add a third H molecule close to the Ti atom, all of the hydrogen molecules are almost parallel to the chain and are bonded molecularly. Interestingly, the H-H bonds turn nearly perpendicular to the chain when the fourth H molecule is attached. The maximum number of H molecules that could be bound to the TiC complex is 5, which corresponds to a hydrogen wt % of 6.5. Note that the adsorption of the fifth H molecule does not affect the degree of the bond-length elongation of the side hydrogen molecules [Fig. 3(b)]. The H-H bond of the top H molecule is 0.756 Å and the bond length between this H molecule and the Ti atom is 2.334 Å, leading to an adsorption energy of 0.29 eV, which is only about half the value of the side molecules (0.66 eV).
Insight into the nature of Ti-H bonding and the orientation of H molecules can be gained from the DOS and molecular orbital (MO) analyses. The metal-dihydrogen binding of TiC(H), shown in Fig. 3(c), stems primarily from the Kubas interaction.18 In the energy range from -10 to -6 eV (region I), the d orbitals of Ti are hybridized with the orbitals of the hydrogen molecules, resulting in charge transfer from the hydrogen to the metal. Region II of Fig. 3(c) highlights that the hybridization of the Ti-d orbitals with the H -antibonding orbitals is responsible for keeping the side H molecules parallel and the top H molecule perpendicular to the TiC chain. The isosurface plots of these states, as shown in Fig. 3(d), clearly confirm that the bonding is dominated by the overlap between the Ti-d and H-s states. To better understand the orientation of the side H molecules, we also plotted the deformation charge in Fig. 3(a). The charge distribution around the lower H sites, marked in yellow, indicates that these H atoms lose a greater amount of charge. With the help of Mulliken charge population analysis, we identified that the upper H atom has a charge of 0.04 e, while the lower one carries 0.12 e due to charge transfer from H to the C atom. Thus, the distances between the Ti and the upper or lower H depend on the strength of the Coulomb interaction between the H and C atoms. The Kubas interaction is likewise found in TiC(H). Figure 3(e) shows that the Ti-d orbitals interact with the orbitals of the H molecules. The hybridization between the Ti-d orbitals and the side H molecules’ -antibonding orbitals plays a role in keeping them planar [Fig. 3(f)], in contrast to the situation in TiC(H). The top hydrogen in TiC(H) interacts with Ti rather weakly compared to the rest of the hydrogen molecules, as the relevant MO is more heavily polarized toward the C than toward the H [see Fig. 3(g)]. The Kubas interaction mentioned above leads to variations in the effective charge of the Ti atom. The Mulliken charge of the Ti atom varies from positive to negative as the number of hydrogen molecules increases. When the fifth (sixth) H is adsorbed onto TiC (TiC), the Ti atom carries -0.13 (-0.49) electrons.
The above results demonstrate that the interaction between one Ti atom and the carbon chain is very strong, and the binding energy of H on TiC is optimal for room temperature application. Now the question arises: can the storage performance of TiC be influenced by introducing another Ti? We first discuss the second Ti atom saturating the free end of TiC (n= 5, 8). The binding energy for Ti to the TiC (TiC) is 4.31 eV (4.39 eV), which is obviously much large than that of Ti on the sp carbon nanostructures.2, 3, 4, 5 However, further functionalization through Ti atoms bridging other C-C bonds makes the chain unstable, due to the fact that the small separation between Ti atoms leads to stronger Ti-Ti coupling. At the end, the linear chains are transformed to the clusters composed of C and Ti atoms, which is also confirmed by other theoretical studies.19 We therefore only consider two Ti atoms each capping one end of the C.
Next, we studied the H storage capacity of the TiCTi complex. The TiCTi and TiCTi can adsorb 12 and 10 H molecules, respectively, with average binding energies of 0.53 and 0.55 eV/H. As shown in Figs. 4(a) and 4(b), the gravimetric densities of stored H molecules for TiCTi becomes 13.3 and 9.4 wt % when n=5 and n=8, respectively. Analysis of the electronic structures of the TiCTi-(H) and TiCTi-(H) shows that the Kubas interaction between H and TiC mentioned above also applies to the TiCTi complex. Note that introduction of the second Ti atom to TiC not only keeps the number of adsorbed H per Ti on odd- or even-numbered chains, but also remains the binding energies of the H with metal atoms .
Because controlling the number of atoms per chain in the synthesis of carbon chains may be difficult with today’s technology, it is important to know whether the results reported above for TiC and TiC are applicable to other Ti-capped carbon chains and how they vary with chain length. To this end, we have also studied the potential of the TiC and the TiC complexes as storage media, and we found that the number of adsorbed H molecules is six and five, respectively. The corresponding bond lengths of H and Ti-H, as depicted in Fig. 5(a) and 5(b), are very similar to those in TiC(H) and TiC(H). Furthermore, the average binding energies of the H molecules to TiC and TiC are, respectively, 0.54 and 0.59 eV/H, which are very close to the values found in the case of TiC(H) and TiC(H). The bonding mechanisms between Ti and H presented above also hold for these chains. This suggests that the length of the chains does not affect the hydrogen adsorption performance. In summary, the finding that a single Ti atom-capped, odd-numbered (even-numbered) carbon chain can bind up to six (five) hydrogen molecules is very general.
While the above results are promising for isolated TiC systems, one can imagine terminating the other end of these complexes with suitable graphitic nano-objects. These new structures not only represent the typical interface in realistic nanostructures produced by cluster beam deposition.20 Here we choose C fullerene as the end-capping candidate, as it is the most frequently experimentally synthesized curved sp system. An increasing amount of experimental 21 evidence shows that these hybrid sp+sp carbon-based systems (with linear chains connecting sp-type fragments) exhibit unusual electronic and optical properties. Figure 5(c) illustrates that a TiC chain can be effectively stabilized by termination on C. The calculated binding energy (with respect to the isolated TiC plus the fully relaxed C) is 2.79 eV. As shown in Fig. 5(d), the TiC assembled on the C structure can hold 6 H molecules, with an average binding energy of 0.52 eV/H.
In conclusion, using all-electron DFT calculations, we have shown that each Ti atom adsorbed on even or odd-numbered carbon atomic chains can bind up to five or six hydrogen molecules, respectively. Note that the number of adsorbed H molecules depends only on the type of chain. We propose that the TiC chain terminated effectively on a C fullerene can also store 6 H molecules with an average binding energy of 0.52 eV/H. Recent experiments have produced junctions between a single carbon chain and two fullerenes,22 which has provided a way to synthesis the sp+sp system we proposed here. We hope that the theoretical results presented here will provide a useful reference for the design of high-capacity hydrogen storage materials in the laboratory.
This work was supported by the National Science Foundation of China under Grant no. 10774148, the special funds for the Major State Basic Research Project of China (973) under Grant no. 2007CB925004, and the Knowledge Innovation Program of the Chinese Academy of Sciences. Part of the calculations were performed at the Center for Computational Science of CASHIPS.
- 1 Y. H. Kim, Y. Y. Sun, W. I. Choj, J. G. Kang and S. B. Zhang, Phys. Chem. Chem. Phys., 2009, 11, 11400.
- 2 T. Yildirim and S. Ciraci, Phys. Rev. Lett., 2005, 94, 175501.
- 3 Y. F. Zhao, Y. H. Kim, A. C. Dillon, M. J. Heben and S. B. Zhang, Phys. Rev. Lett., 2005, 94, 155504.
- 4 H. Lee, W. I. Choi and J. Ihm, Phys. Rev. Lett., 2006, 97, 056104.
- 5 E. Durgun, S. Ciraci, W. Zhou and T. Yildirim, Phys. Rev. Lett., 2006, 97, 226102.
- 6 Q. Sun, Q. Wang, P. Jena and Y. Kawazoe, J. Am. Chem. Soc., 2005, 127, 14582.
- 7 Q. Sun, P. Jena, Q. Wang and M. Marquez, J. Am. Chem. Soc., 2006, 128, 9741.
- 8 C. S. Liu and Z. Zeng, Appl. Phys. Lett., 2010, 96, 123101.
- 9 M. Yoon, S. Yang, C. Hicke, E. Wang, D. Geohegan and Z. Zhang, Phys. Rev. Lett., 2008, 100, 206806.
- 10 C. Jin, H. Lan, L. Peng, K. Suenage and S. Iijima, Phys. Rev. Lett., 2009, 102, 205501.
- 11 G. Roth and H. Fischer, Organometallics, 1996, 15, 5766, and references therein.
- 12 N. D. Lang and Ph.Avouris, Phys. Rev. Lett., 1998, 81, 3515; Phys. Rev. Lett., 2000, 84, 358.
- 13 E. Durgun, R. T. Senger, H. Sevincli, H. Mehrez and S. Ciraci, Phys. Rev. B, 2006, 74, 235413; S. Dag, S. Tongay, T. Yildirim, E. Durgun, R. T. Senger, C. Y. Fong and S. Ciraci, Phys. Rev. B, 2005, 72, 155444.
- 14 Y. Wang and J. P. Perdew, Phys. Rev. B, 1991, 44, 13298; J. P. Perdew and Y. Wang, Phys. Rev. B, 1992, 45, 13244.
- 15 B. Delley, J. Chem. Phys., 1990, 92, 508.
- 16 R. H. Crabtree, The Organometallic Chemistry of the Transition Metals. 3rd ed. (Wiley, New York, 2001).
- 17 G. Henkelman and H. Jósson, J. Chem. Phys., 2000, 13, 9978.
- 18 G. J. Kubas, R. R. Ryan, B. I. Swanson, P. J. Vergamini and H. J. Wasserman, J. Am. Chem. Soc., 106, 451.
- 19 L. Largo, A. Cimas, P. Redondo, V. M. Rayón and C. Barrientos, Chem. Phys., 2006, 330, 431.
- 20 M. P. Bogana and L. Colombo, Appl. Phys. A, 2007, 86, 275.
- 21 L. Ravagnan, N. Manini, E. Cinquanta, G. Onida, D. Sangalli, C. Motta, M. Devetta, A. Bordoni, P. Piseri and P. Milani, Phys. Rev. Lett., 2009, 102, 245502; L. Ravagnan, F. Siviero, C. Lenardi, P. Piseri, E. Barborini, P. Milani, C. S. Casari, A. Li Bassi and C. E. Bottani, Phys. Rev. Lett., 2002, 89, 285506; L. Ravagnan, G. Bongiorno, D. Bandiera, E. Salis, P. Piseri, P. Milani, C. Lenardi, M. Coreno, M. de Simone and K. C. Prince, Carbon, 2006, 44, 1518.
- 22 F. Börrnert, C. Börrnert, S. Gorantla, X. Liu, A. Bachmatiuk, J. O. Joswig, F. R. Wagner, F. Schäffel, J. H. Warner, R. Schönfelder, B. Rellinghaus, T. Gemming, J. Thomas, M. Knupfer, B. Büchner and M. H. Rümmeli, Phys. Rev. B, 2010, 81, 085439.