Surface Structure and Reactivity of Anatase TiO2 Crystals with

Riley E. Rex , Fritz J. Knorr , and Jeanne L. McHale .... Bart Kenens , Maha Chamtouri , Remko Aubert , Kana Miyakawa , Yasufumi Hayasaka , Hiroyuki N...
0 downloads 0 Views 2MB Size
Article pubs.acs.org/JPCC

Surface Structure and Reactivity of Anatase TiO2 Crystals with Dominant {001} Facets Sencer Selçuk and Annabella Selloni* Department of Chemistry, Princeton University, Frick Laboratory, Princeton, New Jersey 08544, United States S Supporting Information *

ABSTRACT: Hydrofluoric acid (HF)-assisted hydrothermal/solvothermal methods are widely used to synthesize anatase TiO2 single crystals with a high percentage of {001} facets, which are generally considered to be highly reactive. We have used Density Functional Theory calculations and first principles molecular dynamics simulations to investigate the structure of these facets, which is not yet well understood. Our results suggest that (001) surfaces exhibit the bulk-terminated structure when in contact with concentrated HF solutions. However, (1 × 4)-reconstructed surfaces, as observed in UHV, become always more stable at the typical temperatures, 400−600 °C, used to clean the as-prepared crystals in experiments. Since the (1 × 4)-reconstructed surfaces are only weakly reactive, our results predict that synthetic anatase crystals with dominant {001} facets should not exhibit enhanced photocatalytic activity, consistent with recent experimental observations.

1. INTRODUCTION Increasing the efficiency of photocatalytic materials is a major scientific and technological challenge.1−4 Among the available methods, intense research efforts have been recently focused on tailoring the crystal shape to expose more reactive facets.5−14 One of the most widely used photocatalytic materials is anatase TiO2, for which the equilibrium crystal shape is a truncated square bipyramid15−17 exposing majority {101} and minority {001} facets. Prompted by computational results suggesting that the minority (001) surface could play an important role in the reactivity of anatase,18,19 Yang et al.6 first showed that anatase single crystals with large amounts of (001) surface can be prepared by solvothermal synthesis using hydrofluoric acid (HF) as the controlling agent. (The samples prepared by this procedure will be called F-anatase in the following.) Their work was followed by many others,8−11 and crystals with up to almost 100% (001) surfaces have been synthesized.11 The structure of the resulting (001) surface is still unclear however. In particular, while anatase TiO2(001) is known to exhibit a 1 × 4 reconstruction under vacuum conditions,20 it is not yet established whether the (001) surface of F-anatase samples is also reconstructed or has instead the structure of the bulk terminated surface.21 Here, we use Density Functional Theory (DFT) calculations and first principles molecular dynamics (FPMD)22 simulations to investigate the structure and energetics of H2O and HF adsorption on the 1 × 1 bulk-terminated and 1 × 4 reconstructed anatase (001) surfaces, and determine the corresponding surface free energy diagrams in wide temperature and pressure ranges. We consider extended surface models, which are appropriate for crystals of μm size (as, e.g., in ref 6.), but are less well justified for crystals of nanometric dimensions (as, e.g., in ref 12.). We also assume anatase crystals © 2013 American Chemical Society

without defects and/or added metal impurities, as often are present in real catalytic materials,12,14 since these do not affect the equilibrium structure of extended surfaces.

2. METHODS AND MODELS Our DFT calculations are based on the Generalized Gradient Approximation (GGA) of Perdew−Burke−Ernzerhof (PBE),23 and the plane-wave-pseudopotential scheme as implemented in the Quantum Espresso package.24 Bulk-terminated slabs (Figure 1a) were generated by cleaving bulk anatase and

Figure 1. Side view of (a) bulk-terminated and (b) 1 × 4 reconstructed anatase (001) slabs. The bottom trilayer (shaded area) is fixed at bulk positions in the DFT optimizations. Red and gray spheres represent oxygen and Ti atoms, respectively.

inserting a vacuum layer (∼12 Å) between the two resulting surfaces. To describe the reconstructed surface, we used the “ad-molecule” model of ref 25. (Figure 1b), which agrees well with the available experimental information and is widely accepted. This model is characterized by added TiO2 rows forming “ridges” which run parallel to the [010] direction and Received: February 28, 2013 Published: March 6, 2013 6358

dx.doi.org/10.1021/jp402100v | J. Phys. Chem. C 2013, 117, 6358−6362

The Journal of Physical Chemistry C

Article

Figure 2. Free energy diagram of the anatase (001) surface in the presence of water vapor (center), and selected water adsorption structures: (a) 0.25 ML, and (b) 1 ML water on the 1 × 1 surface; (c) 0.20 ML, and (d) 0.80 ML water on the 1 × 4 surface. Surface formation energies are plotted as a function of the water chemical potential (bottom) or water vapor pressure at different temperatures (top). The shaded area highlights the region of chemical potential corresponding to liquid water (298−500 K).

Figure 3. Free energy diagram of anatase (001) surface interacting with gaseous HF (center), and selected HF adsorption structures: (a) 0.25 ML, and (b) 1 ML HF on the 1 × 1 surface; (c) 0.20 ML, and (d) 0.80 ML HF on the 1 × 4 surface. Surface formation energies are plotted as a function of the HF chemical potential (bottom) or HF pressure at different temperatures (top). The shaded area highlights the region of chemical potential corresponding to HF aqueous solutions (2%-40%).

with cleavage of surface Ti−O bonds. At higher coverage, our results differ from those of previous studies, which were based on smaller unit cells.18,19 Here we find a mixed dissociatedmolecular adsorption structure where only 25% of the water molecules is dissociated, while the others form hydrogen bonds with the surface and between themselves (Figure 2b). This leads to a decrease of the adsorption energy per molecule which is reduced to 0.98 eV at full monolayer coverage. From Figure 2b, we can remark an interesting similarity between the structure of the fully hydrated bulk terminated surface and that of the clean reconstructed surface (Figure 1b). Calculations of water adsorption on the 1 × 4 reconstructed surface show that the latter is significantly less reactive than the 1 × 1 surface, consistent with the relative stabilities of the two surfaces.25,28 The site of highest reactivity for the 1 × 4 surface is the ridge, where water adsorbs dissociatively with an energy of 1.54 eV at low coverage. In the terrace region, water is physisorbed or not adsorbed at all (Figure 2d). The water adsorption energy is 0.63 eV/molecule at full coverage. The computed water adsorption structures and energies described above were used to construct the surface free energy diagram as a function of the water chemical potential (μH2O), or

expose 4-fold coordinated Ti atoms (Ti4c). In-between ridges, “terraces” expose 5-fold coordinated Ti (Ti5c) and 2-fold coordinated O (O2c) atoms, and have the structure of the unreconstructed surface with a lateral compression along [100]. For consistency, (2 × 4) surface supercells (twice the size of the primitive cell of the (1 × 4)-reconstructed surface) were used for all calculations. The approach of ref 26 was used to construct the surface free energy diagrams, with gas phase thermodynamic data taken from JANAF tables.27 Additional details on the methods and models used for the calculations are given in the Supporting Information (SI).

3. RESULTS AND DISCUSSION Anatase (001) in Humid Environment. A large number of structures were optimized in order to accurately determine the most stable configurations for water adsorbed on the bulk terminated and 1 × 4 reconstructed (001) surfaces at different coverages up to one monolayer (ML). In agreement with previous studies,18,19 we find that water readily dissociates on the 1 × 1 surface at low coverage. This adsorption is strongly exothermic, by 2.33 (2.38) eV/H2O in the case of one (two) adsorbed water molecule(s) per 2 × 4 supercell, and occurs 6359

dx.doi.org/10.1021/jp402100v | J. Phys. Chem. C 2013, 117, 6358−6362

The Journal of Physical Chemistry C

Article

Figure 4. Snapshots from FPMD simulations of the bulk terminated (a,b) and reconstructed (c,d) anatase (001) surface in 32% HF aqueous solution (16 HF + 34 H2O molecules).

concentrated aqueous HF solutions, as used for the synthesis of the F-anatase crystals.9 More importantly, Figure 3 shows that the clean 1 × 4 surface is the stable structure at low μHF/ pressure, in particular at pressures below ∼0.01 atm at 900 K. Since, experimentally, heat treatment at 400−600 °C is required to clean the F-anatase samples,6 our calculations predict that the surface obtained after heat treatment is 1 × 4 reconstructed. Interaction of Anatase (001) with Gaseous HF−Water Mixtures and with a Liquid HF−Water Solution. We investigated the adsorption of mixed HF−water monolayers with varying amounts of the two molecules on the bulkterminated and reconstructed surfaces. In these mixtures HF is always found to be adsorbed dissociatively, whereas water forms H-bonds with the surface (SI, Figure S2). The results of these calculations were combined with those for the adsorption of a single species to build the surface phase diagram shown in Figure S3 of the SI. They predict that the reconstructed surface is always preferred over the bulk-terminated one in the presence of a gaseous HF mixture. The above static calculations still provide a very simplified description of the HF aqueous solution environment used for the synthesis of F-anatase. To further investigate the relative stabilities of the two surfaces in such an environment, we performed FPMD simulations at 350 K for these surfaces in contact with a 32% HF aqueous solution. The simulations were started using the most stable geometries established by total energy calculations, and, after 2 ps equilibration, were continued for additional 10 ps. The total potential energy evolutions during the simulations and their corresponding probability distributions are shown in SI, Figure S4 (see also Figure S5). Within the error bars of our simulations, we find no significant difference between the total energies of two surfaces in HF solution. However, analysis of the surface structures during the simulations provides some interesting clues, see Figure 4. While the unreconstructed surface is very stable, the ridges of the 1 × 4 reconstructed surface show important structural fluctuations in HF solution. A complete disruption of the ridges and transformation to the unreconstructed surface cannot be observed since the number of particles is fixed in our simulations. However these results suggest that the reconstructed surface is likely to be unstable in HF aqueous solution.

water pressure at different temperatures, see Figure 2. It is immediately evident from this diagram that the 1 × 4 surface is more stable than the 1 × 1 one in the whole range of water chemical potentials considered. In particular, the fully hydrated 1 × 4 surface is predicted to be the most stable (001) surface configuration at room temperature and a pressure of 0.023 atm, the partial pressure of water vapor in air. Moreover, a partially hydrated 1 × 4 surface is stable in most of the shaded region of the diagram, which corresponds to the range of μH2O values of liquid water. Even though the diagram refers to the surface in contact with water vapor, these results may be considered as an indication that the (001) surface in contact with liquid water is also 1 × 4 reconstructed. Anatase (001) in the Presence of HF(g). Calculations similar to those of water adsorption were performed to study the interaction of the bulk terminated and 1 × 4 reconstructed (001) surfaces with gaseous HF. On the bulk-terminated surface, the structure of adsorbed HF is similar to that of water at low coverage; that is, the adsorption is dissociative and is accompanied by cleavage of a surface Ti−O bond (Figure 3a). The adsorption energy is 2.98 eV per molecule at 0.25 ML, about 25% larger than that of water at the same coverage. In contrast to water, HF adsorption is dissociative also at high coverage, but in this case all Ti−O bonds remain intact (see Figure 3b). The adsorption energy decreases with increasing coverage, down to the value 1.22 eV per molecule at full coverage. On the 1 × 4 surface, HF adsorbs dissociatively on the Ti atoms of the ridge; the computed adsorption energy is 2.14 eV (1.54 eV/HF) in the case of a single (two) HF molecule(s) per cell. At higher coverage, terrace sites are occupied by a mixture of dissociatively adsorbed and H-bonded HF molecules. The adsorption energy is 0.81 eV/HF at full coverage (Figure 3c). Figure 3 shows the free energy diagram of the anatase (001) surface interacting with gaseous HF. In comparison to water (Figure 2), HF adsorption stabilizes the bulk-terminated surface more than the reconstructed one, due to a less favorable F bonding configuration on the latter, where the adsorbed fluorines cause significant distortions of the ridges. In addition, at variance with the bulk-terminated surface, some repulsion between the adsorbed fluorines and the surface oxygens appears to be present on the reconstructed surface. This changes the relative energies of the two surfaces making the 1 × 1 surface essentially degenerate with the reconstructed one in a small range of HF chemical potential (μHF) values characteristic of 6360

dx.doi.org/10.1021/jp402100v | J. Phys. Chem. C 2013, 117, 6358−6362

The Journal of Physical Chemistry C

Article

(10) Fang, W. Q.; Zhou, J. Z.; Liu, J.; Chen, Z. G.; Yang, C.; Sun, C. H.; Qian, G. R.; Zou, J.; Qiao, S. Z.; Yang, H. G. Hierarchical Structures of Single-Crystalline Anatase TiO2 Nanosheets Dominated by {001} Facets. Chem.Eur. J. 2011, 17, 1423−1427. (11) Chen, J. S.; Tan, Y. L.; Li, C. M.; Cheah, Y. L.; Luan, D.; Madhavi, S.; Boey, F. Y.; Archer, L. A.; Lou, X. W. Constructing Hierarchical Spheres from Large Ultrathin Anatase TiO2 Nanosheets with Nearly 100% Exposed (001) Facets for Fast Reversible Lithium Storage. J. Am. Chem. Soc. 2010, 132, 6124−6130. (12) Gordon, T. R.; Cargnello, M.; Paik, T.; Mangolini, F.; Weber, R. T.; Fornasiero, P.; Murray, C. B. Nonaqueous Synthesis of TiO2 Nanocrystals Using TiF4 to Engineer Morphology, Oxygen Vacancy Concentration, and Photocatalytic Activity. J. Am. Chem. Soc. 2012, 134, 6751−6761. (13) Fang, W. Q.; Gong, X.-Q.; Yang, H. G. On the Unusual Properties of Anatase TiO2 Exposed by Highly Reactive Facets. J. Phys. Chem. Lett. 2011, 2, 725−734. (14) Jiang, H. B.; Cuan, Q.; Wen, C. Z.; Xing, J.; Wu, D.; Gong, X.Q.; Li, C.; Yang, H. G. Anatase TiO2 Crystals with Exposed HighIndex Facets. Angew. Chem., Int. Ed. 2011, 50, 3764−3768. (15) Ohno, T.; Sarukawa, K.; Matsumura, M. Crystal Faces of Rutile and Anatase TiO2 Particles and Their Roles in Photocatalytic Reactions. New. J. Chem. 2002, 26, 1167−1170. (16) Shklover, V.; Nazeeruddin, M. K.; Zakeeruddin, S. M.; Barbe, C.; Kay, A.; Haibach, T.; Steurer, W.; Hermann, R.; Nissen, H. U.; Gratzel, M. Structure of Nanocrystalline TiO2 Powders and Precursor to Their Highly Efficient Photosensitizer. Chem. Mater. 1997, 9, 430− 439. (17) Lazzeri, M.; Vittadini, A.; Selloni, A. Structure and Energetics of Stoichiometric TiO2 Anatase Surfaces. Phys. Rev. B 2001, 63, 1−9. (18) Gong, X. Q.; Selloni, A. Reactivity of Anatase TiO 2 Nanoparticles: The Role of the Minority (001) Surface. J. Phys. Chem. B 2005, 109, 19560−19562. (19) Vittadini, A.; Selloni, A.; Rotzinger, F. P.; Gratzel, M. Structure and Energetics of Water Adsorbed at TiO2 Anatase (101) and (001) Surfaces. Phys. Rev. Lett. 1998, 81, 2954−2957. (20) Herman, G. S.; Sievers, M. R.; Gao, Y. Structure Determination of the Two-Domain (1 × 4) Anatase TiO2 (001) Surface. Phys. Rev. Lett. 2000, 84, 3354−3357. (21) Selloni, A. Crystal Growth - Anatase Shows Its Reactive Side. Nat. Mater. 2008, 7, 613−615. (22) Car, R.; Parrinello, M. Unified Approach for Molecular Dynamics and Density-Functional Theory. Phys. Rev. Lett. 1985, 55, 2471−2474. (23) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (24) 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. (25) Lazzeri, M.; Selloni, A. Stress-Driven Reconstruction of an Oxide Surface: The Anatase TiO2 (001)-(1 × 4) Surface. Phys. Rev. Lett. 2001, 87, 24−27. (26) Arrouvel, C.; Digne, M.; Breysse, M.; Toulhoat, H.; Raybaud, P. Effects of Morphology on Surface Hydroxyl Concentration: A DFT Comparison of Anatase- TiO2 and γ-Alumina Catalytic Supports. J. Catal. 2004, 222, 152−166. (27) Nist Chemistry Webbook. http://webbook.nist.gov/chemistry/ (accessed on 22 September 2013). (28) Gong, X.-Q.; Selloni, A.; Vittadini, A. Density Functional Theory Study of Formic Acid Adsorption on Anatase TiO2(001): Geometries, Energetics, and Effects of Coverage, Hydration, and Reconstruction. J Phys. Chem. B 2006, 110, 2804−2811. (29) Pan, J.; Liu, G.; Lu, G. Q.; Cheng, H. M. On the True Photoreactivity Order of {001}, {010}, and {101} Facets of Anatase TiO2 Crystals. Angew. Chem., Int. Ed. 2011, 50, 2133−2137.

4. CONCLUSIONS We have studied the (001) surface structure of anatase crystals prepared by HF-assisted hydrothermal synthesis using first principles calculations and extended surface models, which are well justified for crystals of micrometer size. Our results indicate that the {001} facets of these crystals exhibit a 1 × 1 bulk-terminated structure in concentrated HF solutions. However, (1 × 4) reconstructed surfaces should form after cleaning the samples at 400−600 °C to remove the adsorbed fluorine. Since the 1 × 4 reconstructed surface is only weakly reactive, our results predict that synthetic anatase samples with a large percentage of {001} facets should not exhibit enhanced photocatalytic activity, consistent with recent experimental observations.12,29,30



ASSOCIATED CONTENT

S Supporting Information *

Detailed information about the methodology including convergence criteria, structural data, snapshots from the MD simulation, and the two-component phase diagram. This material is available free of charge via the Internet at http:// pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by DoE-BES, Division of Chemical Sciences, Geosciences and Biosciences under Award DE-FG0212ER16286. We used resources of the Center for Functional Nanomaterials, Brookhaven National Laboratory. We also acknowledge use of the TIGRESS high performance computer center at Princeton University.



REFERENCES

(1) Diebold, U. The Surface Science of Titanium Dioxide. Surf. Sci. Rep. 2003, 48, 53−229. (2) Henderson, M. A. A Surface Science Perspective on TiO2 Photocatalysis. Surf. Sci. Rep. 2011, 66, 185−297. (3) Thompson, T. L.; Yates, J. T., Jr. Surface Science Studies of the Photoactivation of TiO2New Photochemical Processes. Chem. Rev. 2006, 106, 4428−4453. (4) Fujishima, A.; Zhang, X. T.; Tryk, D. A. TiO2 Photocatalysis and Related Surface Phenomena. Surf. Sci. Rep. 2008, 63, 515−582. (5) Jun, Y. W.; Casula, M. F.; Sim, J. H.; Kim, S. Y.; Cheon, J.; Alivisatos, A. P. Surfactant-Assisted Elimination of a High Energy Facet as a Means of Controlling the Shapes of TiO2 Nanocrystals. J. Am. Chem. Soc. 2003, 125, 15981−15985. (6) Yang, H. G.; Sun, C. H.; Qiao, S. Z.; Zou, J.; Liu, G.; Smith, S. C.; Cheng, H. M.; Lu, G. Q. Anatase TiO2 Single Crystals with a Large Percentage of Reactive Facets. Nature 2008, 453, 638−641. (7) Dinh, C. T.; Nguyen, T. D.; Kleitz, F.; Do, T. O. ShapeControlled Synthesis of Highly Crystalline Titania Nanocrystals. ACS Nano 2009, 3, 3737−3743. (8) Yang, X. H.; Li, Z.; Liu, G.; Xing, J.; Sun, C. H.; Yang, H. G.; Li, C. Z. Ultra-Thin Anatase TiO2 Nanosheets Dominated with {001} Facets: Thickness-Controlled Synthesis, Growth Mechanism and Water-Splitting Properties. CrystEngComm 2011, 13, 1378−1383. (9) Wen, C. Z.; Jiang, H. B.; Qiao, S. Z.; Yang, H. G.; Lu, G. Q. Synthesis of High-Reactive Facets Dominated Anatase TiO2. J. Mater. Chem. 2011, 21, 7052−7061. 6361

dx.doi.org/10.1021/jp402100v | J. Phys. Chem. C 2013, 117, 6358−6362

The Journal of Physical Chemistry C

Article

(30) Tachikawa, T.; Yamashita, S.; Majima, T. Evidence for CrystalFace-Dependent TiO2 Photocatalysis from Single-Molecule Imaging and Kinetic Analysis. J. Am. Chem. Soc. 2011, 133, 7197−7204.

6362

dx.doi.org/10.1021/jp402100v | J. Phys. Chem. C 2013, 117, 6358−6362