Ganes Shukria,
Wilson Agerico Diño*ab,
Hermawan K. Dipojonoef,
Mohammad Kemal Agustaef and
Hideaki Kasaiacd
aDepartment of Applied Physics, Osaka University, Suita Campus, Osaka 565-0871, Japan. E-mail: wilson@dyn.ap.eng.osaka-u.ac.jp
bCenter for Atomic and Molecular Technologies, Osaka University, Suita Campus, Osaka 565-0871, Japan
cNational Institute of Technology, Akashi College, Akashi, Hyogo 974-8501, Japan
dInstitute of Industrial Science, The University of Tokyo, Komaba, Meguro, Tokyo 153-8505, Japan
eNational Research Center for Nanotechnology, Institute of Technology Bandung, Jl. Ganesha 10, Bandung 40132, Indonesia
fLaboratory of Computational Materials Design, Research Group of Engineering Physics, Institute of Technology Bandung, Jl. Ganesha 10, Bandung 40132, Indonesia
First published on 12th September 2016
A density functional theory (DFT)-based study of ethylene (C2H4) adsorption on a reduced anatase titanium dioxide (TiO2) (001) surface, i.e., with a surface oxygen vacancy (Ovac), is presented. It was found that C2H4 preferably adsorbs on the Ovac-site. The excess electrons originating from the removed oxygen weaken the C2H4 CC double bond by filling the lowest unoccupied molecular orbital (LUMO) or
of C2H4, and simultaneously enhance the binding of C2H4 to the reduced anatase TiO2 (001). The bonding between the two C2H4 C atoms and the two Ti atoms nearest to the Ovac-site produces two σ-type bonds, leading to the emergence of a new localized mid-gap state. This hybrid
3d defect state can account for the ∼0.7 eV decrease in the band-gap observed in previous optical measurements of carbon-coated TiO2, formed after TiO2 exposure to C2H4 gas (Mater. Lett. 108, 2013, 134). Subsequent calculations on two initial decomposition pathways of C2H4 adsorbed on the Ovac-site show that the C–H bond and C–C bond cleaving require high activation barriers where the C–H bond cleaving is slightly easier compared to the C–C bond cleaving, 2.94 eV and 3.01 eV, respectively, (as calculated using DFT-D2+Ud = 3 eV) and the final states of the two initial decomposition pathways show similar endothermic characteristics. This finding indicates that the surface Ovac-sites tend to favor the formation of molecularly adsorbed Ti-bound sp3-C2H4 (–Ti–CH2CH2–Ti–) as compared to the dissociative adsorption case.
In this paper particular attention is given to the ethylene–titanium dioxide (C2H4–TiO2) interaction. Previous studies of the photocatalytic decomposition of C2H4 using TiO2 found that in the absence of light irradiation, C2H4 did not adsorb nor dissociate with TiO2. Experimental measurement shows no gas pressure change occurred even after long exposure of TiO2 to C2H4. More recently Fourier-transform-infrared measurements on low temperature C2H4 adsorption using a Degussa P25-type of TiO2 showed that C2H4 interacts weakly with TiO2. The C2H4 vibrational frequencies did not undergo significant modification and C2H4 retains its planar configuration upon adsorption.11–14 These results agree with theoretical calculations on C2H4 adsorption on a pristine anatase TiO2 (001) surface.15 It was found that C2H4 adsorbs on anatase TiO2 (001) via a π-complex interaction with the unsaturated five-coordinated Ti (Ti5c). Furthermore, it was found that the dispersion (van der Waals) interaction further enhances the binding of C2H4.
A number of further studies on C2H4–TiO2 have since been conducted, most of them focusing on C2H4 decomposition using TiO2 as a photocatalyst.16–20 Some recent experimental work have tried to modify TiO2 surfaces with carbon by using a hydrocarbon such as C2H4 and acetylene as carbon sources.21,22 Some other researchers have also have attempted to grow carbonaceous structures on TiO2 for solar cell device application.23,24 Indeed, these studies have given many valuable insights to understand the nature of C2H4–TiO2 interaction. However, not much is known regarding the role of intrinsic defects to C2H4 interaction with TiO2. Intrinsic defects, in particular oxygen vacancy (Ovac), are expected to have important roles in determining the physics and chemistry of transition metal oxide surfaces. The excess electrons originating from the Ovac often become a critical factor that governs the interaction of oxide surfaces with various kind of electron-acceptor molecules.25–29 Therefore, it is necessary to study further how C2H4 adsorption – which can be regarded as an olefin prototype – behaves when an Ovac exists in TiO2.
In the present work, TiO2 anatase was chosen because of its relatively higher catalytic activity compared to its rutile counterpart.30–32 For the surface facet, research was concentrated on the bulk-terminated (1 × 1) anatase TiO2 (001) because of its superior catalytic activity, in particular when coupled with the presence of an Ovac.33,34 Indeed, bulk-terminated anatase TiO2 (001) is well known for its highly unstable characteristics, and thus, it usually undergoes surface reconstruction (the so-called (1 × 4)-reconstructed) to compensate for this instability. However, recent experimental work has succeeded in synthesizing anatase nanoparticles with the dominant bulk-terminated (001) facet. This is achieved by using a particular non-metallic atom (typically fluorine), which acts as a shape-controlling agent. Furthermore, this shape-controlling agent can be cleaned using heat treatment to produce a fluorine free surface without changing the basic crystal structure and morphology.35–37 Based on these useful results, an attempt has been made to further explore the interaction between C2H4 and reduced anatase TiO2 (001) surface. In this study, it will be shown that the presence of the Ovac promotes the formation of molecularly adsorbed C2H4. Furthermore, it will also be shown how the coupling between the C2H4-reduced anatase TiO2 (001) can possibly extend the photoabsorption of TiO2 to the visible light range. A similar study has been performed on the more stable (101) facet, as well. There are some differences in the energetics and configurations of the adsorbed C2H4, in particular, on the reduced anatase TiO2 (101) surface with that on the reduced (001) surface. This originates mainly from the difference of the preferable site of the Ovac. The (101) surface favors the formation of the Ovac in the subsurface, which is different to the (001) surface which favors the formation of the Ovac on its outer most layer (this will be explained in Section 3.1). However, it is planned to publish the results in a future work to ensure a more rigorous and concise report.
Approximately 16 Å of vacuum space was utilized and included a dipole correction to eliminate the artificial interaction between the periodically repeated slabs. A c(2 × 2) unit cell, was utilized containing four unsaturated Ti5c atoms, each bonding to two-raised 2-coordinated (O2c) and two-lowered 3-coordinated (O3c) oxygen atoms along the [100] and [010] directions, respectively. Brillouin zone integration was sampled with 6 × 6 × 1 k-points. Each slab consists of four O–Ti–O layers. This amount of layers is sufficient to produces good convergence of the (001) surface structure32,45 and has been widely used in previous studies.46–48 Two upper O–Ti–O layers of the slab were allowed to relax during geometry optimization, keeping the other remaining O–Ti–O layers at the bottom side fixed to its initial bulk position. The adsorption energy is defined as:
Eads = Eadsorbate/surface − (Eadsorbate + Esurface) | (1) |
![]() | (2) |
The damping factor, S6, is 0.75 as the GGA-PBE was utilized in all calculations. The dispersion coefficient (C6) and the vdW radius (R0) for O, C, and H are 0.7 J nm6 mol−1, 1.75 J nm6 mol−1, 0.14 J nm6 mol−1 and 1.342 Å, 1.452 Å, 1.001 Å, respectively.49 For Ti, C6 and R0 are set to be 10.8 J nm6 mol−1 and 1.562 Å, following previous theoretical studies.50
Previous theoretical studies on TiO2 with an Ovac have shown that a common GGA functional failed to describe the experimentally observed localized defect states, because of inherent spurious self-interaction within the DFT formalism. This results in unphysical delocalization of the electronic states.51–53 The utilization of DFT-D2+Ud, however, could overcome this problem. This method is implemented by introducing the on-site Coulomb interaction in the localized d-orbital and exchange interactions, by adding an effective Hubbard-U parameter to better describe the repulsion between the electrons on the same orbital.54,55 Here, Ud = 3 eV was mainly used to determine the bonding mechanism and energetics of the initial decomposition reactions of C2H4 on reduced anatase TiO2 (001). Indeed, the chosen Ud is not sufficient to fully open the band-gap of pristine anatase TiO2 to its experimentally measured value (∼2.4 eV versus ∼3.2 eV). However, there are some considerations that have been taken into account as to why the present Ud was chosen: (1) the present choice of Ud properly describes the energy level of the reduced anatase TiO2 (001) defect state with respect to the conduction band minimum (this will be elaborated further in the next section); (2) it gives a reasonable estimation for the reaction energies of TiO2 with other molecules [e.g., hydrogen (H2)];56 (3) if the Ud-parameter (typically >5 eV) is too high it may open the band-gap more, but simultaneously it may also result in the wrong description of the adsorbate-surface bond when it comes to adsorption properties and also the energetics of corresponding chemical reactions.56 In this work, the hybrid DFT calculation (HSE06)57 was also performed to check the consistency of the result calculated using DFT-D2+Ud, in particular for C2H4 adsorbing on the Ovac site. Finally, the climbing image nudged elastic band method (CI-NEB) was utilized with four to eight transition images to determine the activation energy barriers for the initial C–H and CC bond cleaving of adsorbed C2H4. Vibrational frequency calculations were then used for the confirmation of the transition state (TS).
Evac-formation = ETiO2−x − (ETiO2 + ½ EO2) | (3) |
Site | Evac-formation (eV) | Difference |
---|---|---|
1 | 3.50 | — |
2 | 4.35 | 0.85 |
3 | 4.47 | 0.97 |
4 | 4.80 | 1.30 |
Five different adsorption configurations of C2H4 on the reduced anatase TiO2 (001) with surface Ovac as shown in Fig. 3 were studied. Note that only C2H4 adsorption on the reduced anatase TiO2 (001) with Ovac created at site #1 which is the most feasible site for creating Ovac was considered. Single point calculations were first performed for each of the five configurations, where all degrees of freedom were constrained and only the vertical distance of C2H4 relative to the surface was varied (Fig. S3, ESI†). Based on the one-dimensional potential energy curves obtained from each of the configurations, the distance at which the C2H4-reduced TiO2 (001) surface system gives the lowest energy was further fully optimized. In previous work by Shukri and Kasai, it was found that the long range dispersion interaction shows a significant contribution to the adsorption energy of C2H4. This affects overall C2H4 stability on the pristine TiO2 surface.15 The contribution of vdW interaction was also checked by comparing the calculated adsorption energy of C2H4 on each adsorption site. Indeed, the vdW interaction still contributes a significant part to the overall calculated adsorption energies (Table S1, ESI†). Therefore, in the present study, the long range dispersion interaction (vdW interaction) was also included in the calculations.
The energetics and geometric properties of C2H4 adsorbed on a defective anatase TiO2 (001) surface calculated at DFT-D2+Ud (Ud = 3 eV) are shown in Table 2. C2H4 adsorption on Ti5c, O3c, hollow, and O3c-tilted sites resulted in relatively weak adsorption energies. These weak interactions were also indicated by the geometrical structure of C2H4, having almost no perturbation on its initial sp2-hybridization and adsorbed at a relatively far distance from the reduced TiO2 (001) surface. For molecularly adsorbed C2H4 on those four previously mentioned sites, it can be inferred that vdW is the dominant interaction that stabilizes the adsorption.
Configuration | ΔEads (eV) | dC–C (Å) | dC–H (Å) | dC–Ti (Å) |
---|---|---|---|---|
Ti4c | −0.10 | 1.33 | 1.10 | 3.18 |
O3c | −0.18 | 1.33 | 1.09 | 3.32 |
Ovac | −0.94 | 1.47 | 1.09 | 2.10 |
Hol | −0.32 | 1.33 | 1.09 | 3.87 |
O3c-tilted | −0.15 | 1.33 | 1.10 | 3.40 |
C2H4 adsorption on an O-vacancy site (O2c-site), however, has a quite distinct property when compared to the other four configurations. Upon adsorption, the initial (gas phase) sp2-hybridization of C2H4 changed into sp3-type structure, following the elongation of the C2H4 CC bond. The C
C bond was elongated to ∼0.16 Å compared to its gas phase bond length, indicating a weakened bonding interaction between the two carbon atoms. This weakening also corresponds to the decrement of the C2H4 bond order from the initial double bond (C
C) to become a single bond (C–C). From now on, the weakened carbon–carbon bond of adsorbed C2H4 with the C–C notation will be referred to. In order to further elucidate the mechanism of C2H4 adsorption on the Ovac-site, the orbital-resolved LDOS of the isolated C2H4 gas phase and C2H4 adsorbed on the Ovac-site [Fig. 4 and 5] was plotted. Here, the pz-orbitals of the two C atoms formed the highest occupied molecular orbital (HOMO; πcc)–LUMO
of C2H4. These two frontier molecular orbitals (MOs) often play a significant part in C2H4 interaction with various solid surfaces. By observing the plot of the orbital-resolved LDOS of C2H4 adsorbed on the Ovac site (Fig. 5), it can be seen that there is a strong hybridization between the p-states of C2H4 C–C with the Ti 3d-states, i.e., two adjacent Ti atoms in the vicinity of the Ovac.65 Another noticeable feature is the emergence of a localized occupied state below the Fermi level, which formed by the mixing of the
of C2H4 with the occupied Ti 3d defect states.66 This hybridization implies that charges from occupied defect states of Ti 3d are transferred to the
of C2H4, indicating a back donation from the defect states to C2H4 that subsequently weakened the adsorbed C2H4 C–C bond. The weakening because of the depletion of charges in the C–C region upon adsorption can also be observed, whereas accumulation of charges took place in the two C–Ti bonding and C–H bonding regions (Fig. 5d). The accumulation of charges in the C–H bonding regions subsequently tilted the four hydrogen atoms, and as a consequence the C2H4 minimized the repulsive interactions between its bonding electrons. Fig. 5 also shows the plot of the spin density profile within the energy range of Ti 3d defect states –
derived-bonding localized mid-gap state that shows two σ-type bonding profiles between both C atoms of C2H4 with the two Ti atoms near the vacancy site. These two σ-type interactions (or the so-called di-σ complex) characterize the rehybridization of C2H4 into the sp3-type upon its adsorption on the Ovac site. It is interesting to note that previous experimental work has shown that coating a anatase TiO2 surface using C2H4 and acetylene as carbon sources resulted in the decrease of the TiO2 band-gap by about ∼0.7 eV, subsequently increasing the photoactivity of the C-modified TiO2 under visible light.22 These present results clarified at least one possible origin of this band-gap decrement. The hybrid
– 3d defect state resulting from C2H4 – reduced anatase TiO2 (001) coupling (located at ∼1.0 eV below the conduction band minimum) can account for the increasing photoactivity of C-modified TiO2, thus extending the system photoactivity to the visible light range. Thus, exposing the system to ultraviolet-visible light can trigger electronic excitation from the hybrid
– 3d defect state to the conduction band that accounts for the experimentally observed decreasing band-gap. This is similar to other electron-acceptor molecules adsorbed on a reduced TiO2 surface such as oxygen (O2), where a localized mid-gap state was also observed upon adsorption.67 In addition, a hybrid-DFT calculation was also performed to check the consistency of the present result (C2H4 adsorbed on the Ovac-site) obtained by DFT+Ud (Ud = 3 eV). The optimized geometry and adsorption energy results obtained using hybrid-DFT are in good agreement with the results obtained from DFT+Ud (Ud = 3 eV) (Table 3). This consistency ensures the reliability of the chosen Ud value. To this end, it can be concluded that excess electrons because of the oxygen vacancy play a vital role to stabilize C2H4 adsorption on the reduced anatase TiO2 (001) surface. The present finding is in contrast to the C2H4 adsorption on clean/pristine anatase TiO2 (001), in which C2H4 adsorption is mainly driven by the dispersion (or vdW) interaction.15
Calculation method | Eads (eV) | C–C (Å) | C–H (Å) | C–Ti (Å) |
---|---|---|---|---|
DFT-D2+Ud (Ud = 4 eV) | −0.8 | 1.49 | 1.1 | 2.14 |
DFT-D2+Ud (Ud = 5 eV) | −0.32 | 1.35 | 1.1 | 2.64 |
Hybrid-DFT (HSE06) | −1.12 | 1.50 | 1.10 | 2.10 |
Previous research has shown that by controlling the work function (W), it may be possible to control the redox chemistry of oxide materials as well.25,28 In conjunction with the present work on the C2H4-reduced TiO2 system, the influence of increasing reduced anatase TiO2 (001) work function (WTiO2) on the adsorption of C2H4 was also checked. This idea is tested in an ad hoc manner, whereby the position of the defect state (which simultaneously determines the WTiO2) is varied upon altering the Ud values. Here, the C2H4 adsorption properties on Ovac for two other Ud values (Ud = 4 eV and 5 eV) have been calculated. The different behaviors of C2H4 adsorption calculated with Ud = 5 eV were observed, whereas the result for Ud = 4 eV showed overall similarity to the result calculated with Ud = 3 eV (Fig. S4, ESI†). C2H4 retained its planar geometry and its adsorption energy was also significantly lowered (Table 3). In order to further elucidate the decrease of the adsorption strength of C2H4 on the vacancy site calculated with Ud = 5 eV, the LDOS for the present case were plotted (Fig. 6a). Compared to the DFT-D2+Ud (Ud = 3 eV and 4 eV) cases, some less significant orbital hybridization can be observed, in particular between the Ti 3d defect states with C2H4 (LUMO). Indeed, some fractions of C2H4
can still be observed within the gap states and also near the valence band edge, showing some electronic filling of the
orbital. However, those hybridizations are not strong enough to produce significant stabilization for the C2H4-reduced TiO2 (001) coupling, because most of the broadening of
still occupies the region above the Fermi level, resulting in less significant interplay with the πcc (HOMO). This indicates that charge transfer from the surface to the C2H4 is significantly hindered because of the deeper emergence of the defect states in the energy level (or more tightly bound excess electrons) (Fig. 6b). Nevertheless, the tendency of the di-σ bonding formation from the charge density difference (CDD) plot can still be seen, as the charges accumulate in both the C–Ti regions (Fig. 6c). The results obtained from DFT-D2+Ud (Ud = 5 eV) imply another important point that the deeper the energy level of excess electrons (or the more the excess electrons are tightly bound), the more difficult the charge transfer to the C2H4 LUMO will be. This is similar to if the work function of a substrate is controlled, as has also been demonstrated by previous work of O2 adsorption on reduced rutile TiO2 (110) surface.29
![]() | (4) |
![]() | (5) |
The ‘*’-symbol indicates the adsorbed species on the anatase TiO2 (001) surface; O2c-surf corresponds to the nearest two-fold coordinated oxygen (from the Ovac-site) on the topmost layer. Note that adsorbed C2H4 on Ovac-site was adopted as the initial configuration for both decomposition pathways. The products for the two initial decomposition pathways are determined based on the two criteria: (1) the adsorbed C2H4 dissociated fragments that give the most energetically stable configuration (relative to isolated gas-phase C2H4 and reduced TiO2 (001) slab); (2) the shortest diffusion pathway for C2H4 H-atom upon dehydrogenation. The first pathway (dehydrogenation) results in the breaking of one H atom of C2H4, followed by the formation of a surface hydroxyl with the nearest O2c (two-fold coordination oxygen), whereas the remaining fraction (CHCH2 or vinyl) makes a bond with the two Ti atoms nearest to the Ovac-site. For the second pathway (C–C bond cleaving), C2H4 breaks into two CH2 fragments, where one fragment occupies the Ovac-site and another fragment forms adsorbed formaldehyde with the nearest neighbor O2c (two-fold coordination oxygen) as well. In addition, the consideration of formaldehyde as a C–C cleaving product is based on the previous experimental work in which formaldehyde is frequently observed during hydrocarbon–oxide interaction.68,69
As shown in Fig. 7, it can be see that both decomposition pathways are highly activated reactions with calculated activation barriers of 2.94 eV and 3.01 eV for C–H and C–C bond cleaving, respectively. The calculated reaction energies for both decomposition pathways show similar endothermic characteristics. It is interesting to see that both decomposition pathways show almost similar activation energies. We can attribute this similarity in high activation barrier, partly, to the fundamental characteristics of the C–H and C–C bonds. For dehydrogenation, note that the C–H bond is formed by a strong σ-type interaction, so that breaking this particular bond requires a large amount of energy. It is also important to note that diffusion energy needed by the H to move to the nearest O2c (to form the surface hydroxyl) is also accounted for in the present calculated barrier. This is also in agreement with results from previous work, where atomic H needs to overcome a certain amount of activation barriers to diffuse on the surface of TiO2, and this can range from 0.7 eV–2.06 eV, depending on the diffusion pathways and simulation environments.70,71
![]() | ||
Fig. 7 First two initial decomposition pathways of C2H4 adsorbed on Ovac-site. Reaction (4) corresponds to C2H4 dehydrogenation, whereas (5) corresponds to direct C![]() |
For C–C bond cleaving, a slightly higher activation barrier by about 0.07 eV can be observed, when compared to the C–H bond cleaving. However, similar to the C–H bond, the same argument also prevails in the C–C bond cleaving case, because, once again, the breaking of the σ-type interaction is being dealt with. The C–C bond cleaving consequently requires a certain amount of energy as well. Another important point is that the activation barrier for the C–C bond cleaving also involves the contribution from the energy needed to break the Ti4c–O2c bonding nearest to the Ovac-site, in order to facilitate the formation of adsorbed formaldehyde species. To this end, it can be inferred that this high activation barrier for C–C bond cleaving stems from two contributions: the energy needed to break the C–C σ-bond, and also, subsequently, the Ti4c–O2c bond in order to facilitate H2CO formation as the final product. To some extent, this similarity is analogous to C–H and C–C bond breaking of gas phase C2H6 in which all the CC and CH bonds are composed by σ-type bonding (the energies needed to break CC and CH in C2H6 are ∼3.95 eV and ∼4.36 eV, respectively). Indeed, to the best of our knowledge, there are direct experimental data available that can be referred to and compared with the present results. However, the present tendency of the formation of C–C coupling on the Ovac is similarly observed when reduced rutile TiO2 (110) with surface bridging Ovac is exposed to formaldehyde. By using scanning tunneling microscopy, Zhu et al. observed that, mediated by the Ovac, adsorbed formaldehyde and/or methylene can couple with other adsorbed formaldehyde/methylene species on the surface through several reaction channels, forming C2H4 as the final product as the temperature increased.68,69 The present results support these experimental findings, in which coupled C–C bond (or described here as molecularly adsorbed C2H4) is more favorable when compared to the dissociative adsorption cases when Ovac exists on the surface of TiO2 [as shown by the high dissociation barriers and the comparison of thermodynamic stability of the initial reactant (molecularly adsorbed C2H4) with the final product, namely, C–H and C–C bond cleaving]. An interesting recent work on acetylene (C2H2) adsorption on anatase TiO2 (101) and (001) surfaces by Chen et al. should also be mentioned briefly.72 The published report focuses on the study of C2H2 cyclo-oligomerisation mechanism on anatase TiO2 surfaces. In particular, it was reported that C2H2 spontaneously dissociates upon adsorption on the pristine anatase TiO2 (001) via H–CCH bond cleaving. The dissociated CCH fragment then binds to the surface of the Ti, forming a Ti–CCH surface complex where the H binds to the nearest O2c producing surface OH. Therefore, there seems to be a contradiction between the present results with those obtained by Chen et al., as a relatively high activation energy of C2H4 CH bond cleaving was reported in the current paper. However, it is also important to note some aspects of the two systems in order to come to terms with these seemingly different results. The present results deal with reduced anatase TiO2 (001) with surface Ovac whereas Chen et al. deal with the pristine surface of anatase TiO2 (001). This difference in the surface of interest may already lead to a different adsorption and reaction mechanism as well. A qualitative (and to some extent semi-quantitative) similarity that may be noticed is that in terms of the thermodynamics of the reaction, this result agrees with that obtained by Chen et al., in which reduced anatase TiO2 (001) favors the dehydrogenation or the CH bond cleaving. Indeed, it is more reasonable and interesting to further discuss the C2H2 adsorption on the pristine anatase TiO2 (001) results with the previous work of C2H4 on pristine anatase TiO2 (001).15 However, this is already beyond the scope of the present work. More rigorous work related to the decomposition pathways of unsaturated hydrocarbon (both C2H2 and C2H4) on the pristine (as well as hydroxylated) surface of anatase surfaces is currently being done and will soon be covered in a future report.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c6ra13633h |
This journal is © The Royal Society of Chemistry 2016 |