Next Article in Journal
Correct Stability Condition and Fundamental Performance Analysis of the α-β-γ-δ Filter
Next Article in Special Issue
An Interdisciplinary Approach to the Nanomanipulation of SiO2 Nanoparticles: Design, Fabrication and Feasibility
Previous Article in Journal
Characterizing Situations of Dock Overload in Bicycle Sharing Stations
Previous Article in Special Issue
The Effect of Production Parameters of Oxide Layers on Their Nanostructure, Nanomorphology, and Surface Free Energy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Deeper Insight on the Stability of Water-Induced Reconstruction of Anatase (001) Surface

1
Dipartimento di Scienze di Base e Applicate per l’Ingegneria, University of Rome “La Sapienza”, Via A. Scarpa 14–16, 00161 Rome, Italy
2
Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Appl. Sci. 2018, 8(12), 2522; https://doi.org/10.3390/app8122522
Submission received: 28 November 2018 / Revised: 3 December 2018 / Accepted: 4 December 2018 / Published: 6 December 2018
(This article belongs to the Special Issue Nanocharacterization and Innovation at Nanoscale)

Abstract

:
TiO 2 anatase (001) surface has been indicated for many years as a potential system for water dissociation and hydrogen production. Surface reconstruction periodicity of TiO 2 anatase (001) in water is revised on the basis of the new water induced reconstruction model that accounts for dissociative water adsorption in the first monolayer and self-assembling of surface hydroxyls. The study has been performed in the context of first principles total energy calculations on the basis of state of the art Density Functional Theory. Different surface periodical structures have been studied and compared in terms of residual surface stress and surface reactivity. While a preference seems to emerge for the (2 × 3) surface reconstruction, there are indications that this configuration might not occur spontaneously in bulk water.

1. Introduction

Bulk TiO 2 and TiO 2 surfaces of the various crystal phases (anatase, rutile, and brookite) are being considered since many years for a variety of applications with special mention to environment-and energy-related applications, such as the photo-decomposition of organic pollutants, solar cells, and solar-hydrogen production [1,2,3]. Titanium oxides are also important in surgical implants where the implant surface is oxidized and in contact with the bio-moieties in the human body [4]. All the mentioned systems and problems require the knowledge of the surface-water interaction processes because they imply the usage of either nano-sized TiO 2 or TiO 2 surfaces in water. For this reason, the atomistic scale characterization and modeling of the TiO 2 -water interface is being studied intensively [5,6,7,8]. Among the crystal phases, the metastable anatase form is known to be photo-catalytically more active than the ground state rutile phase, and becomes the most stable phase at the nano-scale. This is why the anatase is increasingly considered. While the (101) is the most stable anatase surfaces, the minoritary (001) has been claimed to be more reactive in water [7,9,10,11,12,13,14] and this is the reason why many articles have focussed on the (001) anatase surface in water. Moreover, recent results have reported on the successful synthesis of TiO 2 particles with minority facets [15,16,17]. The reactivity of the (001) surface with respect to the water adsorption and dissociation has been debated for long time with no conclusive statements yet. While the bare (unreconstructed) TiO 2 (001) surface has been predicted as a dissociative reactive surface for water absorption and, therefore, a good candidate to obtain photo-reactive water splitting for energetic purposes, the experimental results are rather controversial and do not confirm such predictions [18,19,20]. The (001) unreconstructed (dry) surface is made of rows of under-coordinated atoms along the [100] direction formed by sequences of five-fold titanium atoms (Ti 5 c ) and two-fold oxygen atoms (O 2 c ) with alternating long and short Ti 5 c -O 2 c bonds and an inward relaxation of the topmost layer with respect to the unrelaxed slab (Figure 1). There is a general consensus that such a surface should be highly reactive for water dissociation [7,9,10,11,12,13,14] but this is in marked contrast with the experiments that show, quite surprisingly, that this is not the case [18,19,20,21].
Moreover, in-situ experimental findings have shown that stable (001) anatase surface shows a (1 × 4) reconstruction character in vacuum. Atomistic models, based on popular first principles techniques such as the Density Functional Theory (DFT), have been constructed starting from the well known ADM model [22] that is able to reproduce the (1 × 4) surface reconstruction observed in vacuum. However there have been increasing difficulties to explain the experimentally observed surface intertness for water dissociation using the ADM model [21,23]. The reactivity loss has been related to the decreasing of the surface stress with respect to the unreconstructed surface [11,12,13] but, while the ADM model predicts a reduced surface stress due to the addition of TiO 2 molecules on the surface, this appears to be not sufficient because some residual surface reactivity, that is not confirmed by experiments, has been found [23]. Moreover it has been shown also that point defects, such as vacancies, oxygen ad-atoms or interstitials, might even enhance the ADM reactivity [23]. For this reason the ADM model has been further modified proposing that the observed surface reactivity might be related to the oxidized or reduced state of the surface [21]. Regarding the behavior of TiO 2 (001) surface in water, some authors have suggested that water patterned TiO 2 (001) surface might compete with the ADM models from the stability point of view [20,24,25,26,27]. The present authors have recently revised the properties of water patterned (001) anatase surface in a large supercell [27] and two main new findings emerged: (1) Dissociative water adsorption forms a ridge-terrace geometry similar the one of the ADM model. This surface modification is accompanied by the formation of a stabilizing network of hydrogen bonds so that the water induced reconstruction (WIR) geometry shows an apparent (2 × 4) periodicity (see Figure 2); (2) such a water induced surface reconstruction reduces the surface stress at a larger amount than the ADM model and the surface reactivity disappears. However, the above mentioned WIR model was not studied and discussed in terms of the stability hierarchy of the possible various surface periodicities, even though previous studies have indicated the (1 × 3) reconstructed surface as the most stable one [24]. Then, in the present article, we try to go deeper into this aspect concerning the newly found stable ridge geometry of the WIR model.

2. Methods

All the calculations have been performed in the context of state of the art DFT with plane waves basis for both the electron density and wave-functions and periodic boundary conditions (PBCs). We have employed ultrasoft pseudopotentials [28] and the generalized gradient approximation (GGA) as implemented by the Perdew-Burke-Ernzerhof (PBE) scheme [29]. This computational approach has been already successfully employed in the previous literature revealing that PBE is the best GGA exchange-correlation scheme for the TiO 2 -water system [30]. Some technical aspects are worth to be mentioned such as the energy cutoff values employed for the wave functions and the electron density (60 Ry and 600 Ry respectively) and the threshold values employed for convergence of the atomistic optimization procedure, namely 10 4 Ry/Bohr for residual forces and 10 8 Ry for the total energy. It is important to evidence the necessity to add a dispersion correction term to the total energy in order to represent correctly the water-surface interaction. Generally speaking, this term is always expected to improve the energetics even though its inclusion causes a slight shortening of bond lengths and lattice constants [31,32] that, however, remain within the DFT-PBE error bar. In the particular case of the TiO 2 -water system, we have recently shown that the inclusion of the dispersion term is mandatory to get a reliable description of the adsorption and the surface energy values for the reconstructed surface. For all these reasons, we have included in the Hamiltonian a dispersion energy term by employing the well known Grimme-D2 correction [33]. Special mention deserves to be devoted to the well known self-interaction problem affecting standard DFT resulting in a wrong description of the electron localization of 3d orbitals. This inaccurate description affects the electronic properties of the system and should be taken into account especially for transition metals such as Ti. A common method to correct the self-interaction error is the so called DFT+U scheme where the orbital occupancy is forced to integer values by the inclusion of the Hubbard term in the Hamiltonian [20]. However we have shown that, concerning the present system, the inclusion of the Hubbard term does not affect the energetic properties we are interested in because these properties, depending on the difference between the energy values obtained from different configurations (see below), are rather insensitive to the inclusion of the Hubbard term [27]. For this reason the Hubbard term has not been included in our model.
Further details of the adopted computational scheme, together with the comparison of the bulk anatase properties obtained with the present scheme and other ones from the previous literature, can be found in a recent article [27].
The anatase (001) surface has been modeled through a slab geometry with the surface orthogonal to the z axis. The anatase bulk crystal was cleaved to expose the (001) surface and then relaxed. The slab contains four layers that provide an adequate geometry to describe in a correct way the relaxation pattern of the topmost layers in a real surface [23,25,32,34]. Moreover the supercell size along the z-axis has been chosen to have a ∼40 Å-thick vacuum region in order to avoid artifacts due to the PBCs.
An important choice is to keep fixed just the atoms belonging to bottom-most layer in order to reproduce the bulk atomic positions. Indeed, the less rigorous choice, adopted by many authors, of keeping fixed two bottom layers does not guarantee the correct relaxation patten in the slab [27].
Because the aim of the present article is to give a deeper insight to the water induced surface reconstruction periodicity and to the relevant water adsorption energy values, we have adopted different supercell sizes depending on the periodical ridge-terrace geometry investigated, namely (2 × 2), (2 × 3), (2 × 4), and (2 × 6) where the indexes refer respectively to the ridge direction [010] and to the orthogonal one [100]. For all the supercell sizes and surface reconstructions we have adopted a rigorous Monkhorst-Pack k-point sampling scheme (MP) [35] with denser and more demanding grids than the ones typically employed in the recent literature: our reference was to keep the grid at least as dense as the 2 × 2 × 1 MP grid for a (4 × 4) surface geometry. In other words we have employed a 4 × 4 × 1 MP grid for the (2 × 2) slab geometry, a 4 × 3 × 1 MP grid for a (2 × 3) surface, a 4 × 2 × 1 for (2 × 4) surface, and lastly 4 × 2 × 1 for (2 × 6). Lastly, the artificial electric field across the slab generated by periodic boundary conditions (PBC) along the z-axis is handled by using a self-consistent dipole correction scheme [36]. All the calculations have been performed with the Quantum ESPRESSO (QE) package [37].

3. Results

The formation energy of the unreconstructed surface γ is calculated according to the following equation:
γ = E s l a b n T i O 2 E T i O 2 A γ f = 2 E s l a b E s l a b f n T i O 2 E T i O 2 2 A
where γ f , E s l a b , E s l a b f are the surface formation energy of the unrelaxed slab and the total energy values respectively of the relaxed and the unrelaxed slab, n T i O 2 is the number of TiO 2 units in the slab, E T i O 2 is the energy of the TiO 2 unit calculated in the anatase bulk structure, and, finally, A is the surface area [27]. The unreconstructed surface energy is, of course, independent on the surface periodicity employed, i.e., on the n value (with n = 2, 3, 4, 6) of the (2xn) supercell. At the present level of the theory and using the above specified MP sampling scheme in the reciprocal space, we get γ = 1.105 J/m 2 , which is quite a large value reflecting the reactivity of the ideal unreconstructed surface.
Concerning the (2xn) water reconstructed surface, the surface energy can be calculated provided the average adsorption energy of dissociated waters is known: Indeed, if this is calculated as
E ( 2 x n ) a = E ( 2 x n ) w e t E ( 2 x n ) d r y n H 2 O E H 2 O n H 2 O
where E ( 2 x n ) w e t and E ( 2 x n ) d r y are the total energy values of respectively the water reconstructed surface and the unreconstructed one, n H 2 O is the number of water molecules dissociatively adsorbed and E H 2 O is the energy of an isolated water molecule, the surface energy can be obtained through:
γ ( 2 x n ) = γ + n H 2 O E ( 2 x n ) a A
The resulting values of the average adsorption surface energy are reported in Table 1.
The above values, which have been obtained in vacuum, can be further corrected by taking into account the dispersion energy contribution that should be considered in water. In bulk water, indeed, the adsorption process and energy, and the corresponding surface energy, must contain also the contribution from the dispersion energy change due to the water molecules that move from the bulk water to the surface where they are dissociatively adsorbed. This aspect has been detailed in a recent article [27]: According to this model, an estimate of such correction can be done with the following assumptions:
  • Provided that the total energy of the real system is obtained by summing the DFT and the dispersion energy terms, E = E D F T + E d , the interaction between the water molecules in bulk water is accounted for by the dispersion energy only ( E n H 2 O D F T = n E H 2 O D F T );
  • the dispersion energy change occurring in a real system when one water molecule is dissociatively adsorbed from the bulk water to the surface is obtained in a Born Oppenheimer Molecular Dynamics (BOMD) run as detailed recently [27]; this contribution is supposed to be due only to the change of the dispersion energy interaction between water molecules and between the TiO 2 slab and water molecules (the dispersion energy change due to the TiO 2 re-arrangement is neglected); BOMD simulations were performed on on the anatase unreconstructed and fully hydrated (001) surface with the mixed Gaussian and plane wave (GPW) method as implemented in the CP2K package [38]. The details and technicalities of the simulations can be found in a recent article that was mostly focussed on the study of the water structure [39];
  • the contribution to the dispersion energy change due to the interaction between the adsorbed (and dissociated) water and the TiO 2 slab is measured from the ground states of the water reconstructed and the unreconstructed surfaces;
  • the dispersion energy of one water molecule in bulk water is measured and averaged from BOMD of bulk water at room temperature; for this term we have obtained E H 2 O d 0.0245 ± 2 10 3 R y ;
  • the contribution to the adsorption energy coming from the interaction between the slab and the bulk water above it is dominated by the medium-long range interactions and thus by the dispersion part of the energy.
Using the above assumptions it can be shown [27] that the adsorption energy, appropriately corrected with the dispersion terms that should be present in bulk water, casts into:
E a = E a + Δ E s l a b b u l k w a t e r d 2 E H 2 O d
with
Δ E s l a b b u l k w a t e r d = E ( T i O 2 ) ( n 2 ) ( H 2 O ) m d + E 2 ( H 2 O ) d ( n 2 ) ( H 2 O ) m d E ( T i O 2 ) n ( H 2 O ) m d
where ( H 2 O ) m ( d ) stands for molecular (dissociated) water molecules and E ( T i O 2 ) q ( H 2 O ) m d , E 2 ( H 2 O ) d p ( H 2 O ) m d are respectively the dispersion energy of the TiO 2 slab interacting with q water molecules in bulk water and the dispersion energy term of the two adsorbed and dissociated water molecules interacting with p water molecules in bulk water. On the basis of this simple model, we have calculated the dispersion corrected values of the average adsorption and surface energy for (2xn) reconstructed surfaces with n = 2–4 as listed in Table 2:
The surface energy values corrected with the dispersion energy contribution from the bulk water show a new stability hierarchy where the second most stable surface reconstruction is the (2 × 4) one instead of the (2 × 2) one. Moreover, while the surface energy values increase with respect to the one in vacuum, the difference between (2 × 3) and (2 × 4) is reduced. At this point, we have employed the (2 × 6) surface reconstruction in order to reproduce in a doubled supercell the (2 × 3) surface reconstruction. We have attempted to obtain a spontaneous dissociative adsorption in the middle of the terrace so as to reproduce a (2 × 3) periodicity and we have checked several different initial configurations. However we have only obtained molecular water adsorption on the terrace, as reported in Figure 3.

4. Discussion

Two main results of the previous section are worth to be analyzed and discussed: first the stability hierarchy between the three WIR surface geometries that have been considered, namely (2 × 2), (2 × 3), and (2 × 4), is changed when the missing part of the dispersion energy, the one that is typically due to the bulk water above the surface, is included. Hence, the results obtained in vacuum might change in bulk water as evidenced by the change of the positions in the stability hierarchy between (2 × 2) and (2 × 4): In bulk water, indeed, the (2 × 4) surface is more stable than the (2 × 2) that, on the contrary, appears to be more stable in vacuum. This behavior is related to the different values measured for the Δ E s l a b b u l k w a t e r d energy data shown in Table 2: the dispersion energy change, due to the interaction with the bulk water, between the reconstructed and the unreconstructed surfaces increases with the surface area (also this term per unit area increases with the area as well). Then, according to the Equation (4), the dispersion energy of bulk water E H 2 O d contributes to lower the adsorption energy. The surface energy is obtained from the corrected values of the average adsorption energy by dividing for the surface area and this reduces the positive correction term (that lowers the absolute value of the adsorption energy) at larger extent in larger areas. In other words, the correction per unit area is greater for (2 × 2) than for (2 × 4) thus resulting in the inversion of stability. It is interesting to note that the correction due to the dispersion energy results in closer surface energy values between (2 × 3) and (2 × 4). Moreover, the evidence that the adsorption in the (2 × 6) reconstructed surface is molecular also in the middle of the terrace, evidence that we have got after various attempts to obtain a dissociative adsorption, indicates that the (2 × 3) reconstruction should be, most probably, thermally activated because we expects the existence of an energy barrier for a further dissociative adsorption along [100] as soon as the ridge is already formed.
Lastly, because the adsorption on the (2 × 3) and (2 × 4) water reconstructed surfaces is molecular [27], it is plausible that in bulk water the hierarchy could be further changed if one takes into account the molecular adsorption energy occurring on the available sites. Indeed also molecular adsorption events might release some energy per unit area due to the small, but not negligible, surface relaxation. For this reason we have sampled all the possible adsorption sites on the (2 × 3) and (2 × 4) surface to check whether such inversion could occur in bulk water.
The adsorption sites sampled are reported in Figure 4 and the molecular adsorption values obtained in vacuum are reported in Table 3.
The adsorption energy values shown in Table 3 are larger for (2 × 4) than for (2 × 3) and this is consistent with the lower residual stress found in (2 × 3). This circumstance, anyway, tells us that the surface relaxation induced by molecular adsorption onto (2 × 4) should be larger than the one occurring on (2 × 3). In addition, considering the (2 × 2) case, we see that although this surface evidences the lowest residual tensile stress, the surface energy is the largest one occurring in bulk water teaching us that the further inclusion of the dispersion energy correction makes the relationship between the residual surface stress and the surface energy in bulk water not that straightforward as in vacuum.
In addition, the inclusion of the dispersion energy correction would presumably play the same role also for molecular adsorption as the one found for dissociative adsorption. Considering the close values of the surface energy in bulk water for (2 × 3) and (2 × 4) water induced reconstructed surfaces, this circumstances, together with the presumable existence of an activation energy for dissociative water on a (2 × 6) surface, makes the question of finding the most stable WIR configuration in bulk water still open.

5. Conclusions

In summary we have revised the stability hierarchy of water induced reconstruction periodicities of anatase (001) surfaces on the basis of the newly found (2xn) geometry. Using the dispersion energy correction, we have shown that in bulk water the stability hierarchy found in vacuum is changed with the (2 × 4) surface reconstruction being more stable than the (2 × 2) one. The (2 × 3) surface appears to be the most stable one also in bulk water but the surface energy difference between (2 × 3) and (2 × 4) decreases with respect to the vacuum case. The molecular adsorption energy values on both the terrace and the ridge sites are larger (in absolute values) on the (2 × 4) than on the (2 × 3) in accordance with the surface energy values and the residual surface stress. This means that the molecular adsorption induces a larger relaxation pattern onto the (2 × 4) than onto the (2 × 3) surfaces. Thus we can argue that, if one considers the dispersion energy involved in the molecular adsorption, the same surface energy correction behavior could occur for molecular adsorption (to a lesser extent perhaps) as the one found for dissociative adsorption. Another point is that the water adsorption onto the (2 × 6) surface is not spontaneously dissociative, meaning that some activation energy must exist. For all these reasons, we can state that the real stability hierarchy of WIR anatase (001) surface in bulk water should not be considered as a closed question yet. Free energy measurement in the context of ab-initio molecular dynamics simulations could help to clarify this important issue.

Author Contributions

G.Z. has contributed with supervision, conceptualization, investigation, data curation, writing and editing. E.V. has contributed with conceptualization and data curation.

Funding

This research received no external funding.

Acknowledgments

We kindly acknowledge the CRESCO/ENEAGRID High Performance Computing infrastructure and its staff [40]. CRESCO/ENEAGRID High Performance Computing infrastructure is funded by ENEA, the Italian National Agency for New Technologies, Energy and Sustainable Economic Development and by Italian and European research programmes, see http://www.cresco.enea.it/english for information. We also acknowledge Lorenzo Agosta for providing the BOMD outputs.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DFTDensity Functional Theory
GGAGeneralized Gradient Approximation
PBEPerdew-Burke-Ernzerhof
WIRWater Induced Reconstruction
BOMDBorn-Oppenheimer Molecular Dynamics

References

  1. Chen, X.; Mao, S.S. Titanium Dioxide Nanomaterials: Synthesis, Properties, Modifications and Applications. Chem. Rev. 2007, 107, 2891–2959. [Google Scholar] [CrossRef] [PubMed]
  2. Fujishima, A.; Zhang, X.; Trykc, D.A. TiO2 photocatalysis and related surface phenomena. Surf. Sci. Rep. 2008, 63, 515–582. [Google Scholar] [CrossRef]
  3. Yang, Z.P.; Cheng, H.E.; Chang, I.H.; Yu, I.S. Atomic Layer Deposition TiO2 Films and TiO2/SiNx Stacks Applied for Silicon Solar Cells. Appl. Sci. 2016, 6, 233. [Google Scholar] [CrossRef]
  4. Ivanoff, C.J.; Widmark, G.; Hallgren, C.; Sennerby, L.; Wennerberg, A. Histologic evaluation of the bone integration of TiO2 blasted and turned titanium microimplants in humans. Clin. Oral Implant. Res. 2002, 12, 128–134. [Google Scholar] [CrossRef]
  5. Vittadini, A.; Casarin, M.; Selloni, A. Chemistry of and on TiO2-anatase surfaces by DFT calculations: A partial review. Theor. Chem. Acc. 2007, 117, 663–671. [Google Scholar] [CrossRef]
  6. Sun, C.; Liu, L.M.; Selloni, A.; Lu, G.Q.M.; Smith, S.C. Titania-water interactions: A review of theoretical studies. J. Mater. Chem. 2010, 20, 10319–10334. [Google Scholar] [CrossRef]
  7. Gala, F.; Agosta, L.; Zollo, G. Water Kinetics and Clustering on the (101) TiO2 Anatase Surface. J. Phys. Chem. C 2016, 120, 450–456. [Google Scholar] [CrossRef]
  8. Agosta, L.; Zollo, G.; Arcangeli, C.; Buonocore, F.; Gala, F.; Celino, M. Water driven adsorption of amino acids on the (101) anatase TiO2 surface: An ab initio study. PCCP 2015, 17, 1556–1561. [Google Scholar] [CrossRef] [PubMed]
  9. He, Y.; Tilocca, A.; Dulub, O.; Selloni, A.; Diebold, U. Local ordering and electronic signatures of submonolayer water on anatase TiO2 (101). Nat. Mater. 2009, 8, 585–589. [Google Scholar] [CrossRef] [PubMed]
  10. Agosta, L.; Gala, F.; Zollo, G. Water Diffusion on TiO2 Anatase Surface. AIP Conf. Proc. 2015, 1667, 020006-1–020006-9. [Google Scholar]
  11. Gong, X.Q.; Selloni, A. Reactivity of Anatase TiO2 Nanoparticles: The Role of the Minority (001) Surface. J. Phys. Chem. B Lett. 2005, 109, 19560–19562. [Google Scholar] [CrossRef] [PubMed]
  12. Fazio, G.; Ferrighi, L.; Valentin, C.D. Spherical versus faceted anatase TiO2 nanoparticles: A model study of structural and electronic properties. J. Phys. Chem. C 2015, 119, 20735–20746. [Google Scholar] [CrossRef]
  13. Lazzeri, M.; Vittadini, A.; Selloni, A. Structure and energetics of stoichiometric TiO2 anatase surfaces. Phys. Rev. B 2001, 15, 155409. [Google Scholar] [CrossRef]
  14. Vittadini, A.; Selloni, A.; Rotzinger, F.P.; Grätzel, M. Structure and Energetics of Water Adsorbed at TiO2 Anatase (101) and (001) Surfaces. Phys. Rev. Lett. 1998, 81, 2954. [Google Scholar] [CrossRef]
  15. Hengerer, R.; Bolliger, B.; Erbudak, M.; Gratzel, M. Structure and stability of the anatase TiO2 (101) and (001) surfaces. Surf. Sci. 2000, 460, 162–169. [Google Scholar] [CrossRef]
  16. Han, X.; Kuang, Q.; Jin, M.; Xie, Z.; Zheng, L. Synthesis of Anatase Nanosheets with Exposed (001) Facets via Chemical Vapor Deposition. J. Am. Chem. Soc. 2009, 131, 3152–3153. [Google Scholar] [CrossRef]
  17. Lee, W.J.; Sung, Y.M. Synthesis of Anatase Nanosheets with Exposed (001) Facets via Chemical Vapor Deposition. Cryst. Growth Des. 2012, 12, 5792–5795. [Google Scholar] [CrossRef]
  18. Yu, J.; Low, J.; Xiao, W.; Zhou, P.; Jaroniec, M. Enhanced Photocatalytic CO2-Reduction Activity of Anatase TiO2 by Coexposed {001} and {101} Facets. J. Am. Chem. Soc. 2014, 136, 8839–8842. [Google Scholar] [CrossRef]
  19. Tachikawa, T.; Yamashita, S.; Majima, T. Evidence for Crystal-Face-Dependent TiO2 Photocatalysis from Single-Molecule Imaging and Kinetic Analysis. J. Am. Chem. Soc. 2011, 133, 7197–7204. [Google Scholar] [CrossRef]
  20. Angelis, F.D.; Valentin, C.D.; Fantacci, S.; Vittadini, A.; Selloni, A. Theoretical Studies on Anatase and Less Common TiO2 Phases: Bulk, Surfaces, and Nanomaterials. Chem. Rev. 2014, 114, 9708–9753. [Google Scholar] [CrossRef]
  21. Wang, Y.; Sun, H.; Tan, S.; Feng, H.; Cheng, Z.; Zhao, J.; Zhao, A.; Wang, B.; Luo, Y.; Yang, J.; Hou, J.G. Role of point defects on the reactivity of reconstructed anatase titanium dioxide (001) surface. Nat. Commun. 2013, 4, 2214. [Google Scholar] [CrossRef] [PubMed]
  22. Lazzeri, M.; Selloni, A. Stress-Driven Reconstruction of an Oxide Surface: The Anatase TiO2(001) (1 × 4) Surface. Phys. Rev. Lett. 2001, 87, 266105. [Google Scholar] [CrossRef] [PubMed]
  23. Shi, Y.; Sun, H.; Saidi, W.A.; Nguyen, M.C.; Wang, C.Z.; Ho, K.; Yang, J.; Zhao, J. Role of Surface Stress on the Reactivity of Anatase TiO2 (001). J. Phys. Chem. Lett. 2017, 8, 1764–1771. [Google Scholar] [CrossRef] [PubMed]
  24. Ignatchenko, A.; Nealon, D.G.; Dushane, R.; Humphries, K. Interaction of water with titania and zirconia surfaces. J. Mol. Catal. A Chem. 2006, 256, 57–74. [Google Scholar] [CrossRef]
  25. Selcuk, S.; Selloni, A. Facet-dependent trapping and dynamics of excess electrons at anatase TiO2 surfaces and aqueous interfaces. Nat. Mater. 2016, 15, 1107–1112. [Google Scholar] [CrossRef] [PubMed]
  26. Selçuk, S.; Selloni, A. Surface Structure and Reactivity of Anatase TiO2 Crystals with Dominant {001} Facets. J. Phys. Chem. C 2013, 117, 6358–6362. [Google Scholar] [CrossRef]
  27. Vitale, E.; Zollo, G.; Agosta, L.; Gala, F.; Brandt, E.G.; Lyubartsev, A. Stress Relief and Reactivity Loss of Hydrated Anatase (001) Surface. J. Phys. Chem. C 2018, 122, 22407–22417. [Google Scholar] [CrossRef]
  28. Vanderbilt, D. Soft self-consistent pseudopotentials in a generalized eigenvalue formalism. Phys. Rev. B 1990, 41, 7892. [Google Scholar] [CrossRef]
  29. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef]
  30. Santra, B.; Michaelides, A.; Scheffler, M. On the accuracy of density-functional theory exchange-correlation functionals for H bonds in small water clusters: Benchmarks approaching the complete basis set limit. J. Chem. Phys. 2007, 127, 184104. [Google Scholar] [CrossRef] [Green Version]
  31. Moellmann, J.; Ehrlich, S.; Tonner, R.; Grimme, S. A DFT-D study of structural and energetic properties of TiO2 modifications. J. Phys. Condens. Matter 2012, 24, 424206. [Google Scholar] [CrossRef] [PubMed]
  32. Araujo-Lopez, E.; Varilla, L.A.; Seriani, N.; Montoya, J.A. TiO2 anatase’s bulk and (001) surface, structural and electronic properties: A DFT study on the importance of Hubbard and van der Waals contributions. Surf. Sci. 2016, 653, 187–196. [Google Scholar] [CrossRef]
  33. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H–Pu. J. Chem. Phys. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  34. Esch, T.; Gadaczek, I.; Bredow, T. Surface structures and thermodynamics of low-index of rutile, brookite and anatase - a comparative DFT study. Appl. Surf. Sci. 2014, 288, 275–287. [Google Scholar] [CrossRef]
  35. Monkhorst, H.J.; Pack, J.D. Special points for Brillouin-zone integrations. Phys. Rev. B 1976, 13, 5188. [Google Scholar] [CrossRef]
  36. Bengtsson, L. Dipole correction for surface supercell calculations. Phys. Rev. B 1999, 59, 12301. [Google Scholar] [CrossRef]
  37. Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Dabo, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 2009, 21, 395502. [Google Scholar] [CrossRef]
  38. Hutter, J.; Iannuzzi, M.; Schiffmann, F.; VandeVondele, J. CP2K: Atomistic simulations of condensed matter systems. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2014, 4, 15. [Google Scholar] [CrossRef]
  39. Agosta, L.; Brandt, E.G.; Lyubartsev, A.P. Diffusion and reaction pathways of water near fully hydrated TiO2 surfaces from ab initio molecular dynamics. J. Chem. Phys. 2017, 147, 024704. [Google Scholar] [CrossRef]
  40. Ponti, G.; Palombi, F.; Abate, D.; Ambrosino, F.; Aprea, G.; Bastianelli, T.; Beone, F.; Bertini, R.; Bracco, G.; Caporicci, M.; et al. The role of medium size facilities in the HPC ecosystem: The case of the new CRESCO4 cluster integrated in the ENEAGRID infrastructure. IEEE HPCS 2014, 6903807, 1030–1033. [Google Scholar]
Figure 1. Side view of the unreconstructed anatase (001) surface. Along the [100] direction, long and short Ti 5 c -O 2 c bonds are evidenced. Red: Oxygen; Grey: Titanium.
Figure 1. Side view of the unreconstructed anatase (001) surface. Along the [100] direction, long and short Ti 5 c -O 2 c bonds are evidenced. Red: Oxygen; Grey: Titanium.
Applsci 08 02522 g001
Figure 2. The water induced (2 × 4) reconstruction occurring in a 4 × 4 supercell: (a) Side view, (b) top view. In the top view, the periodical arrangement of the hydrogen bonds along the ridge are evidenced with the arrow indicating the OH direction of the hydroxyl aligned along [100]. Red: Oxygen; Grey: Titanium, white (a) or cyan (b): Hydrogen.
Figure 2. The water induced (2 × 4) reconstruction occurring in a 4 × 4 supercell: (a) Side view, (b) top view. In the top view, the periodical arrangement of the hydrogen bonds along the ridge are evidenced with the arrow indicating the OH direction of the hydroxyl aligned along [100]. Red: Oxygen; Grey: Titanium, white (a) or cyan (b): Hydrogen.
Applsci 08 02522 g002
Figure 3. Initial (a) and final (b) adsorption configurations of one water molecule on adsorption sites at the middle of the terrace in a (2 × 6) reconstructed surface. Red: Oxygen; Grey: Titanium, white: Hydrogen. The adsorption is molecular.
Figure 3. Initial (a) and final (b) adsorption configurations of one water molecule on adsorption sites at the middle of the terrace in a (2 × 6) reconstructed surface. Red: Oxygen; Grey: Titanium, white: Hydrogen. The adsorption is molecular.
Applsci 08 02522 g003
Figure 4. Molecular water adsorption sites onto the (2 × 4) and (2 × 3) water reconstructed surfaces. Terrace and ridge sites are indicated as t and r respectively.
Figure 4. Molecular water adsorption sites onto the (2 × 4) and (2 × 3) water reconstructed surfaces. Terrace and ridge sites are indicated as t and r respectively.
Applsci 08 02522 g004
Table 1. Average adsorption and surface energy values for the three periodical reconstructions. The residual stress values along [100] and [010] are also reported.
Table 1. Average adsorption and surface energy values for the three periodical reconstructions. The residual stress values along [100] and [010] are also reported.
Surface Reconstruction E ( 2 xn ) a (eV) γ ( 2 xn ) (J/m 2 ) σ xx (kbar) σ yy (kbar)
(2 × 2)−1.380.32−12.38−0.76
(2 × 3)−2.160.30−17.061.71
(2 × 4)−2.610.37−20.502.7
(2 × 6)−2.980.55−25.292.93
Table 2. Average adsorption and surface energy values corrected with the dispersion energy contribution from bulk water.
Table 2. Average adsorption and surface energy values corrected with the dispersion energy contribution from bulk water.
Surface Reconstruction Δ E slab bulk water d (eV) E a (eV) γ (J/m 2 )
(2 × 2)−0.0302−1.060.51
(2 × 3)0.0057−1.820.42
(2 × 4)0.0476−2.250.47
(2 × 6)0.1355−2.580.62
Table 3. Molecular water adsorption energy values onto the (2 × 3) and (2 × 4) water reconstructed surface.
Table 3. Molecular water adsorption energy values onto the (2 × 3) and (2 × 4) water reconstructed surface.
Site(2 × 3) eV(2 × 4) eV
t 1 −0.455−0.542
t 2 −0.329−0.367
t 3 −0.644−0.672
t 4 −0.486−0.528
t 5 -−0.542
t 6 -−0.528
r 1 −0.344−0.411
r 2 −0.689−0.705

Share and Cite

MDPI and ACS Style

Zollo, G.; Vitale, E. A Deeper Insight on the Stability of Water-Induced Reconstruction of Anatase (001) Surface. Appl. Sci. 2018, 8, 2522. https://doi.org/10.3390/app8122522

AMA Style

Zollo G, Vitale E. A Deeper Insight on the Stability of Water-Induced Reconstruction of Anatase (001) Surface. Applied Sciences. 2018; 8(12):2522. https://doi.org/10.3390/app8122522

Chicago/Turabian Style

Zollo, Giuseppe, and Eugenio Vitale. 2018. "A Deeper Insight on the Stability of Water-Induced Reconstruction of Anatase (001) Surface" Applied Sciences 8, no. 12: 2522. https://doi.org/10.3390/app8122522

APA Style

Zollo, G., & Vitale, E. (2018). A Deeper Insight on the Stability of Water-Induced Reconstruction of Anatase (001) Surface. Applied Sciences, 8(12), 2522. https://doi.org/10.3390/app8122522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop