Abstract
We study the adsorption of hydrogen molecules on a titanium atom supported by a benzene molecule using generalized gradient corrected Density Functional Theory (DFT). This simple system is found to bear important analogies with titanium adsorption sites in (8, 0) titanium-coated single-walled carbon nanotubes (SWNTs) (T. Yildirim and S. Ciraci, 2005) In particular, we show that up to four molecules can coordinate to the metal ion center, with adsorption patterns similar to those observed in Ti-SWNTs and no more than one molecule dissociating in the process. We analyze in detail the orbital interactions responsible for Ti-benzene binding and for the electron transfer responsible for the dissociation. We find the latter to involve a transition from a triplet to a singlet ground state as the hydrogen molecule approaches the adsorption site, similar to what has been observed in Ti-SWNTs. The total Ti--binding energy for the first dissociative addition is somewhat inferior (~0.4βeV) to the value estimated for adsorption on Ti-SWNTs. We analyze in detail the orbital interactions responsible for the binding.
1. Introduction
Developing safe, cost-effective, and practical means of storing hydrogen is crucial for the advancement of hydrogen and fuel-cell technologies [1β3]. In this context, titanium-decorated organic compounds have received attention for their potential use as high-capacity hydrogen storage materials [4β9]. More in specific, studies based on generalized gradient DFT and first-principle molecular dynamics simulations have indicated that a single-titanium atom supported on an insulating () SWNT can bind up to four hydrogen molecules [10]. Since high titanium coverages on SWNTs are achievable in experimental conditions, these materials could provide means to reach, or even exceed, the minimum hydrogen storage capacity required for technological applications.
Adsorption of one H molecule to a titanium site was found to proceed with a vanishing activation energy and to involve interactions between carbon atoms, titanium, and the hydrogen molecule which are essentially short ranged and localized at the absorption site. It was proposed that the mechanism responsible for the H bonding involves d orbitals of the Ti atom, the antibonding of a H molecule, and p orbital at carbon sites.
In this paper we present a DFT study that elucidates the catalytic role played by the Ti atom in the dissociation of the H molecule and in the bonding mechanism of the Ti atom and the four hydrogen molecules, characterized as βan unusual combination of chemisorption and physisorptionβ [10]. According to the results presented here, the same number of hydrogen molecules can be stored in essentially the same configuration introducing a different support for the Ti catalyst, the benzene molecule. The benzene molecule thus offers a hexagonal ring as support for Ti which is electronically very similar to the () SWNT, so that bonding of H molecules can be studied in detail using well-established molecular energy decomposition techniques and programs [11].
The paper is organized as follows. Details about the calculations are given in Section 2. The adsorption and dissociation of a H molecule on benzene-supported titanium is analyzed in Section 3.1, and the possible addition patterns for two and four molecule additions are studied in Section 3.3. Our results are briefly summarized in Section 4.
2. Computational Details
Calculations have been performed using the Amsterdam Density Functional program (ADF) [12β14]. The molecular orbitals (MOs) were expanded in large uncontracted sets of Slater-type orbitals containing diffuse functions of TZ2P quality. A small frozen core (1s for C, 1sβ2p for Ti) was used in the present calculations. Equilibrium structures were optimized using analytical gradient techniques, with geometries and energies calculated at the OLYP [15, 16] level of theory. Our choice for the functional is motivated by its superior performance for the relative energies of different spin states of various transition-metal complexes in the gas phase, in particular when high-spin states are involved [17β19]. Interactions between H molecules and the CH-Ti complexes were analyzed using the molecular fragment method that is a standard tool in energy decomposition methods for interactions between molecules, as detailed in Reference [11]. We use here primarily the description of the orbitals of the overall system expressed as linear combinations of fragment orbitals. The obvious fragments here are Ti, benzene and H molecules. The orbitals of these fragments are obtained (in the exact geometry they have in the total system) in separate calculations. These orbitals are then used as basis functions in the overall calculation. This affords an analysis of the overall orbitals in terms of percentages of the contributing fragment orbitals.
3. Results and Discussion
3.1. Orbital Interactions and Bonding in CH-Ti
The structure of the CH-Ti complex has been theoretically [20β27] and experimentally [28] studied on the basis of DFT calculations and laser evaporation methods. Experimental evidence indicates the ground state to be a triplet in a planar C configuration. We show in Figure 1 the orbital interaction diagram of the benzene molecule and the Ti atom obtained at the OLYP level of theory for the optimized triplet C configuration. In this analysis, the CH-Ti (C) compound was divided in two fragments: the isolated Ti atom and the benzene support, their geometries fixed to the relaxed interacting configuration. From now on, we will take the -plane as the plane parallel to CH. When interacting with the planar CH, the 3d set of the metal atomic orbitals (AOs) are split into three levels of E, A, and E symmetry, the 3da (d) orbital, two 3e (d and Ti-3d), and two 3e (d and d) orbitals, respectively. The important bonding interaction is between the benzene orbitals of E symmetry and the Ti-3d orbitals, which is understandable considering the relatively large overlap between these favorably oriented orbitals. This leads to significant stabilization of the bonding combination and destabilization of the antibonding combination. There is a weaker interaction between the orbitals of E symmetry and the Ti-3d and Ti-3d orbitals, which are parallel to the benzene plane and little interaction of the 3d with the much lower lying CH- orbital of A symmetry. The 3da (Ti-3d) AO is somewhat stabilized by the small hybridization with the 4s and 4p of the same symmetry. The orbital pattern for the d orbitals is reminiscent of the one for sandwich complexes, but a peculiarity of the CH-Ti system is the presence of a low-lying 4s orbital, hardly displaced from its atomic position [26, 29]. This is related to Ti being coordinatively very unsaturated; in highly coordinated (hexa- or penta-coordinated) metal complexes the is invariably very much destabilized by the orthogonality requirement on the occupied ligand orbitals.
The four electrons of Ti are available to go into the CHTi-de and -da orbitals. Because there is only a small gap between these orbitals, the βhigh-spinβ configuration (e)3 (da)1 results, leading to a triplet ground state of E symmetry. This is Jahn-Teller distorted, so that the degeneracy of the e MO is lifted. Our OLYP calculations predict a triplet-C configuration as ground state of the CH-Ti compound, with the quintet-C (e)2 (da)1 (4s)1 and triplet-C (e)3 (da)1 (4s)0 configurations lying higher in energy by 0.372 and 0.242βeV, respectively. Previously, both triplet-C [20, 24, 26] and quintet-C [21β23, 25, 27, 28] states have been theoretically predicted as possible ground states. However, a quintet-C ground state was predicted only by the B3LYP method, and previous studies on the reaction products formed between group 5 transition metals (V, Nb, Ta) and benzene revealed the B3LYP method to calculate the correct ground state only for the CH-Nb complex [29]. It was already suggested that triplet-C is the real ground state of the CH-Ti complex [26] at DFT-PBW91 level of theory, and our OLYP evidence leads us to agree with such suggestions. Conclusive evidence from experimental data on the ground state of the CH-Ti complex is, however, still lacking.
As a consequence of the Jahn-Teller distortion, the planar C CH structure undergoes a deformation, characterized by a ruffling of the C ring (up or down displacement of the C atoms along the -axis), Figure 1. Also, an elongation of all the CβC bonds in the CH-Ti compound of up to 4% is observed, because of the charge donation from the CH bonding e states to the metal atom. The C planar configuration may distort in two different ways towards a C configuration, one characterized by a displacement along the -axis of two para carbon atoms towards Ti, labeled as [4 + 2], the other one by displacement of the other four carbon atoms towards Ti, labeled as [2 + 4] configuration. In both cases, the distance along the -axis between the Ti and the displaced C atoms is reduced from the C value of 1.66βΓ to 1.56 and 1.59βΓ values of the [2 + 4] and [4 + 2] C configurations, respectively. Energies reported in Figure 2 are relative to the C structure with planar CH.
Figure 3 shows the electronic structure of the CHTi (Cβ [2 + 4]) complex. The Jahn-Teller distortion from the planar C to the [2 + 4] distorted C structures causes a splitting of the degenerate CHTi de levels (Figure 1) into the 8a and 3a states as shown in Figure 3. These now lie lower in energy than the Ti-3d (he 9a in Figure 3) with 3da (Figure 1) parentage. In the corresponding [2 + 4] configuration of the benzene, the e LUMO splits into an a and an a orbital, which stabilize the 8a (Ti-3d) and 3a (Ti-3d) MOs by admixing in a bonding fashion. The -axis is perpendicular to the direction of the 2βC atoms of the [2 + 4] structure (and again the 2βC atoms of the [4 + 2] structure). The triplet C β [2 + 4] configuration is thus characterized by a doubly occupied 8a-Ti-3d state and singly occupied 9a-Ti-3d and 3a-Ti-3d states. In the alternative C β [4 + 2] Jahn-Teller configuration, the splitting of the e is reversed, the 3a-Ti-3d state lies lower in energy, while the Ti-3d lies in higher position. Thus, the electron configuration in the distorted [] geometry is characterized by a doubly occupied 3a-Ti-3d state of the metal atom and a singly occupied Ti-3d. However, in both configurations, the 3d state is singly occupied (9a state in both [2 + 4] and [4 + 2] configurations). In all cases the empty Ti-4s (10a) is present, which potentially could act as acceptor orbital for electrons from the approaching H.
3.2. H2 Dissociative Adsorption on CH-Ti
Taking the triplet Cβ [2 + 4] as the representative electronic ground state of the CH-Ti compound, we now turn to investigate the role played by the Ti addition in the adsorption of the four H molecules on the CH-Ti compound reported by Yildirim and Ciraci [10]. The first H addition proved to lead to dissociation of the corresponding hydrogen molecule, while an additional three molecules could be weakly bonded to the Ti atom on the SWNT surface. We therefore obtain here, for the CH-Ti compound, the same H adsorption pattern observed in the () Ti-SWNT studied in Reference [10].
Figure 4 shows the total energies of the CHTiH complex as a function of the distance between the metal atom and the center of mass of the approaching H. All energies are relative to the value of the triplet curve at large distance, that is, the chosen energy zero is the sum of the energies of isolated ground state H and ground state triplet CHTi[]. For each value of the constrained Ti-H distance a full geometry optimization was carried out in both triplet and singlet states by relaxing to equilibrium all unconstrained degrees of freedom. The H approach turns out to be with the center of mass on the -axis and the H bond axis perpendicular to the -axis and in the direction of the two equivalent C atoms of the [2 + 4] configuration, which is the -axis. Similar to the H adsorption on a SWNT [10] we find the CHTi-H complex to be in a triplet ground state at large Ti-H distances and to undergo a transition to the singlet at 1.72βΓ . At variance with the barrierless dissociative adsorption of H on the () Ti-SWNT studied in Reference [10], a finite barrier (0.17βeV) is computed here for the H addition to CHTi, given by the difference between the energy at the spin transition and the energy at the minimum of the triplet spin surface. The absence of a barrier for H dissociation in the () Ti-SWNT is likely to arise from a selective stabilization of the singlet state relative to the triplet in the Ti-SWNT compared to CHTi. In Figure 4 this would result in the singlet curve joining smoothly the triplet one at Ti-H distances lower than ~2βΓ .
The H lowest unoccupied orbital (LUMO) 1 has B symmetry in C, so it cannot mix with the de orbitals (8a1 and 3a2 in Figure 3) or with the or 4s (both A), but it mixes very strongly with the -derived Ti-3d (the 6b of the CHTi fragment). This orbital is a typical frontier orbital which is hybridized with the 4p and thereby acquires a large amplitude towards the incoming H. The strong interaction of 1 and results in a low-lying stabilized strongly mixed CHTiH-6b level (Tables 1 and 2 and Figure 5). As the H molecule approaches the CHTi compound, this level stabilizes so much that it becomes successively occupied with two (spin-paired) electrons, which may be pictured as coming from the originally singly occupied 9a and next 3a of CHTi[2 + 4] of Figure 3. This is a typical example of a bond-breaking reaction, where the initially empty strongly HβH-antibonding orbital becomes occupied. The stabilization of the 1 which makes the mixing with 3d very strong is of course concomitant with the stretching of the H molecule (the bond breaking). The whole process only requires a very low barrier of 0.17βeV.
The interaction of H with a transition metal is often weak, because occupied metal orbitals have Pauli repulsion with the occupied orbital of H, preventing H from coming so close that the interaction of a metal d orbital with the Hββ orbital can become strong. Typically another occupied metal d orbital then interacts only weakly with the high-lying H, resulting in weak backdonation. In the present case the coordinatively unsaturated Ti atom, with relatively few d electrons (as an early transition metal in the 3d series), offers a route to strong interaction and H bond breaking. The 1 mixes (allowed by symmetry) with the 4s and the 3d states (9a and 10a of the CHTi fragment) but not with the 3d (8a), see Table 2. The antibonding combination of 3d and , embodying the Paul repulsion, is destabilized in the process and loses its electron to the downwards moving 6b. The loss of the electron from the antibonding orbital implies disappearance of the Pauli repulsion, and the H can approach unhindered to optimize the interaction of with in the now fully occupied 6b, which leads to the bond breaking. At the same time, there is a low-lying bonding combination of the (mostly) orbital with the Ti-3d, resulting in the orbital 8a, which remains fully occupied (with the original two electrons of ). This orbital is further stabilized by favorable mixing with the Ti 4s (10a) and the Ti-3d (note the percentages of 12% and 15%, resp. in the 8a, Table 2). This contributes to the force pulling H to Ti. So ultimately the two electrons of H end up in the Ti-(H) bonding orbital 8a and out of the four Ti electrons two go into the Ti-(H) bonding orbital 6b (HβH bond breaking), while the two remaining ones reside in the largely nonbonding Ti-3d (9a). The final resulting electron configuration has only doubly occupied orbitals, that is, it is a singlet state. The transition from the triplet to the singlet state occurs at a distance of 1.72βΓ between the Ti atom and the incoming H molecule.
Figure 5 shows the orbital interactions between the CHTi complex and the adsorbed H molecule in the singlet equilibrium geometry. For analysis the systems was divided into two fragments: the CHTi C β [2 + 4] complex described before and the hydrogen molecule with the HβH distance lengthened to the dissociated value (2.61βΓ ). At such a long HβH distance the 1 (A symmetry in C) and 1 (B symmetry) orbitals have little bonding/antibonding character left and are almost degenerate. As described above, these two H-based orbitals enter the orbital energy diagram as the states 8a and 6b, below the 3d orbitals, with a good deal of mixing of these H orbitals and the metal orbitals. The internal caption in the figure shows the composition of the CHTiH C β [2 + 4] 8a and 6b MOs.
3.3. Adsorption of Two and Four H2 Molecules
Turning now to the addition of two H molecules, we observe that Figure 6 shows two adsorption patterns, one in which one of the H molecules is dissociated (Figure 6(a)) and one in which purely molecular adsorptions are observed (Figure 6(b)). In the dissociative configuration (a) the distance between the centers of mass of the first and the second H molecule and the Ti is predicted to be 0.76 and 1.81 Γ , respectively, with a lengthened 0.84βΓ bond for the second nondissociated molecule and a HβH distance of 3.22βΓ for the dissociated one. The HβH distance in a free H molecule is 0.74βΓ . In the configuration (b), the two H molecules were added in a position other than the on-top adsorption site of the Ti. The optimized geometry predicted a displacement between the center of mass of the two H molecules and the Ti atom of 0.98βΓ and 1.77βΓ along the -axis and the -plane, respectively. The adsorbed H bond length was calculated to be 0.78βΓ . The first configuration is energetically favored by 0.38βeV over the second.
(a)
(b)
Table 3 shows the composition of the MOs involved in the donation/backdonation bonding mechanisms. In the configuration (a), the system was divided again in two fragments: CHTiH and the second undissociated H molecule. In the configuration (b), the system was divided into a CHTi and two H fragments. In both configurations, the added H molecules bind weakly to the Ti and do not dissociate with bonding energies of 0.10 and 0.20βeV per H molecule, respectively.
In the configuration (a), the dissociation of the second H molecule is hindered by the double occupancy of the CHTiH-8a orbital, which prevents the second H molecule from approaching the adsorption site by the Pauli repulsion (occupied-occupied orbital repulsion of 1 with 8a). Strong mixing of the high-lying H2ββ with a suitable occupied metal fragment orbital is therefore blocked, and significant stabilizing mixing between the 1 and the empty 10a-3d and the 11a-4s is also prevented. Charge donation from the CHTi complex to the 1 of the second H molecule is thus avoided, (which could take place in the first dissociative H addition and caused the HβH bond breaking). The empty 10a and 11a orbitals, with large 3d and 4s character, do act to some extent as acceptor orbitals (cf. the 7% 10a mixing with the 1 in Table 3), but bond breaking would require the antibonding 1 orbital to be filled.
In the configuration (b), the dissociation of the two H molecules does not occur since the orbital overlaps between the CHTi and H remain too weak in this configuration. For instance, the calculated overlap between 9a and 1 was 0.114, while a 0.765 value was obtained for the CHTi-H dissociative adsorption described before. The two H molecules should approach the CHTi fragment more closely in order to build better overlaps but would then experience substantial repulsion. The first dissociation in configuration (a) occurs because the antibonding H orbital 1 gets filled through its strong interaction with the exposed Ti-3d. In general, dissociation requires that electrons are dumped into the antibonding (bond-breaking) orbital. This does not occur in configuration (b) because of the repulsion between CHTi-8a and 1's. In addition, bond breaking requires occupation of the 1 orbitals. In the depicted configuration the 1 orbitals cannot, however, establish strong interactions with unoccupied CHTiH orbitals and gain sufficient stabilization to bring about its occupation through the mechanism described above. A weak donation/backdonation mechanism is therefore responsible for the bonding between the Ti atom and the two H molecules also in this second configuration. Note that the CHTi-6b, 6a and 3a states indicated in the right panel of Table 3 consist mainly of the Ti-3d, a lower lying CH-e and the Ti-3d/C-p states, respectively.
We next consider the adsorption of four H molecules on the CHTi complex. In this case, two stable four-H adsorption patterns were computed (Figure 7). We denote these two different configurations as TiH + 3H (Figure 7(a)) and Ti + 4H (Figure 7(b)). These equilibrium structures are clearly related to the corresponding geometries in the two-H addition patterns (Figure 6). In both cases, no more than one H molecule undergoes dissociative adsorption, and when dissociation occurs an additional undissociated molecule resides on top of the adsorption site. In the configurations of type (b), the on-top site is left vacant. Similar configurations were found to be stable when using SWNT support in place of CH for the Ti catalyst [10]. The final optimized TiH + 3H structure is 0.32βeV lower than Ti + 4H. In the very symmetric Ti + 4H configuration, all the four molecules stay intact, with an average HβH bond distance and bond energy of 0.82βΓ and 0.23βeV, respectively.
(a)
(b)
Table 4 shows a comparison between the benzene and SWNT supports. It is important to notice that the results presented in this paper agree with the hypothesis of Yildirim and Ciraci [10] about the characteristic properties of a Ti atom supported by a hexagonal-carbon-based framework and its hydrogen adsorption properties. Also, our orbital analysis substantiates the picture of the TiβCβH interaction proposed by the authors for the SWNT support. There are nonetheless important differences between the two supports. In the first instance the adsorption energy of the first dissociative H addition is much reduced in the benzene support (from 0.83βeV to 0.37βeV). In addition, the dissociative addition is an activated process which suggests, according to the argument of Section 3.1, that the singlet state is stabilized in the SWNT compared to the benzene support. Furthermore, when four molecules are adsorbed, in the benzene-Ti system, one is dissociated, the others are not. In the SWNT none is dissociated. The very symmetric SWNT-Ti-4H configuration was found to lie lower in energy by 0.10βeV than the SWNT-TiH3H [10]. It may therefore be argued that the yield of molecular H adsorption in the CH-Ti system would be lower than in the SWNT (one need to recombine one molecule when releasing hydrogen). Finally, the average Ti-H adsorption energy decreases from 0.54βeV for the SWNT support to 0.23βeV for the CH. In conclusion, although the CH support for the Ti addition qualitatively reproduces the dominant local contribution of the large gap () SWNT, the different boundary conditions of the two materials, finite system in the case of CH, infinite periodic in the case of the SWNT, lead to sizable differences in the adsorption energies. The result is a favored Ti + 4H configuration for the periodic SWNT support of the Ti atom and the more localized TiH + 3H as the favored pattern for CH. We stress that, in the TiH + 3H configuration, the first H adsorption is dissociative, while the other three are not. From a material modelling viewpoint, the CH-Ti complex can be seen as an βalmost convergedβ model system for the SWNT-Ti with respect to cluster size. Nevertheless, CH-Ti can store the same number of H molecules per Ti atom, achieving a 6% weight percentage of stored hydrogen. Indeed, some electronic relationship between SWNT and CH may be expected from the modern theory of the macroscopic polarization [30]. According to such theory, the quadratic spread of the manybody electronic wavefunction in condensed phase has an upper bound, in the presence of a finite gap : where is the so called localization length [30], defined by a unitary operator based on a Berry phase in place of the position operator, and is the electron mass. Qualitatively, (1) can be read as βthe larger the gap, the more localized the electrons areβ. Indeed, () SWNT is a large gap insulator nanotube, thus limiting the quadratic spread of the manybody electronic wavefunction. This means that, as normal for insulators, the interaction between Ti adatom and the SWNT may be expected to consist of a dominant local contribution from the C hexagon surrounding the metal atom. When we isolate the carbon hexagon support of the Ti atom, the CH is obtained by hydrogen passivation of the dangling bonds. Benzene can be therefore provide the simplest model system in which the local effects occurring in the large gap insulator () SWNT may be represented.
4. Summary and Conclusions
We studied the adsorption of one, two, and four H molecules on a single titanium atom supported on a benzene molecule. In all cases we found one H molecule to adsorb dissociatively and the remaining one(s) to bind to the metal ion center through weak-charge transfer interaction. The dissociative addition was found to involve a transition from the triplet to the singlet energy surface, similar to what has been found in the adsorption of on a titanium atom supported on a () SWNT. At variance with the latter, this was, however, found to be an activated process, with a barrier of 0.17βeV. The -binding energy was estimated to be ~0.4βeV lower than in the SWNT. The main emphasis of this work has been on a detailed analysis of the orbital interactions responsible for the H binding to Ti and for the dissociation of one H molecule. Typically, the coordination of H to a transition metal atom (or ion) is very weak, because the occupied orbital will have Pauli repulsion with occupied metal orbitals. This prevents H from coming close and interacting strongly by way of its high-lying , which could lead to breaking of the HβH bond. In the Ti-benzene system, however, the situation is different due to the small number of d electrons and due to fact that the interaction with the benzene βpreparesβ the Ti 3d orbitals for the interaction with H. Notably, the Pauli repulsion, which is embodied in the antibonding interaction of the occupied of with an occupied mostly orbital, is relieved by the loss of the two electrons from the rising antibonding combination of the and . These electrons are transferred to an in-phase (bonding) combination of the (originally empty) Ti and ββ orbitals, which comes down due to the stabilizing interaction. We also find that the empty Ti 4s orbital can interact favorably with , which adds to the force pulling the in. In essence, Ti is specially favorable because the combination of the small number of four d electrons and the special orbital interactions with benzene lead to just the right electronic structure for easy dissociation of . Our results elucidate the electronic structure reasons for the high potential of Ti on carbon supports as new and efficient materials for hydrogen storage.
Acknowledgments
This work was supported by the Netherlands Organization for Scientific Research through the ACTS Programme for Sustainable Hydrogen. Computer resources were provided by the Netherlands' Scientific Research Council (NWO) through a grant from Stichting Nationale Computerfaciliteiten (NCF). Support was also given by the WCU (World Class University) program through the Korea Science and Engineering Foundation funded by the Ministry of Education, Science and Technology of the Republic of Korea (Project no. R32-2008-000-10180-0).