Water on the MgO(001) Surface: Surface ... - ACS Publications

Jun 5, 2015 - Milan Ončák, Radosław Włodarczyk, and Joachim Sauer*. Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 100...
0 downloads 0 Views 1MB Size
Letter pubs.acs.org/JPCL

Water on the MgO(001) Surface: Surface Reconstruction and Ion Solvation Milan Ončaḱ , Radosław Włodarczyk, and Joachim Sauer* Institut für Chemie, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany S Supporting Information *

ABSTRACT: The interaction of water with the MgO(001) surface under ambient conditions is investigated by density functional theory combined with statistical thermodynamics. For water loadings of more than one monolayer, we show that the standard structure model, a fully hydroxylated surface, needs to be revised. Reconstructed surfaces, involving hydrated/hydroxylated Mg2+ ions above the surface, are more stable. These findings provide a consistent picture for surface hydroxylation between low and high water coverage that is in agreement with available XPS data.

K

We show that this model with the fully hydroxylated topmost MgO(001) layer is not a stable structure. Instead, we present alternative structural models that feature hydrated/hydroxylated Mg2+ ions above the surface. Using density functional theory (DFT) and statistical thermodynamics, we examine the stability of different surface structures as a function of temperature and pressure. We compare our results with the X-ray photoelectron spectra (XPS) that yield different O 1s signals for O2− ions of the MgO structure, for OH− ions formed on dissociation of H2O, and for molecularly adsorbed H2O.6,13−16 For both ultrahigh vacuum and ambient pressure conditions, the structure models we propose are consistent with the XPS data. Using semiempirical quantum-chemical calculations, Jug et al.24 showed that the hydroxylated monolayer structure4 shown in Figure 1a is not favorable at ambient temperature and suggested another OH monolayer model with OH groups located in bridging sites between two Mg2+ ions; see the Supporting Information (SI).25 We examined both structures by DFT and found that neither of them represents a minimum on the potential energy surface but rather high-order saddle points. Hydroxyl groups recombine easily when the symmetry constraints are lifted and small displacements are made along an imaginary frequency mode, resulting in a structure with only one-quarter of the water monolayer dissociated. (See Table 1.) Hence, a simple monolayer of OH− (Figure 1a) is not a realistic model. For the MgO(001) surfaces with at least one water monolayer, computational studies14,19−26 have found that water tends to partially dissociate on the MgO(001) surface.

nowing the oxide−water interface in atomic detail is prerequisite to understanding diverse and important phenomena such as nucleation, corrosion and dissolution of oxides, mineral weathering, and the surface chemistry of dust particles. Magnesium oxide with its simple rock-salt structure is a popular model system and the interaction of water with various surfaces has been studied,1−3 both experimentally4−18 and computationally,9,13,14,19−26 in a wide range of temperatures and for pressures reaching from high vacuum to ambient ones. For low water coverage on the MgO(001) surface,5,8−12,14,17,18 a clear picture emerges: At low temperatures (100−180 K), water forms a layer with c(4 × 2) symmetry, which transforms to p(3 × 2) symmetry at higher temperatures (185−221 K). In these structures, 1/5 and 1/3, respectively, of the water molecules are dissociated. At higher temperatures, water molecules desorb from the surface. For water coverages beyond a monolayer, the standard structure model, proposed 1965 by Anderson et al.,4 assumes a fully hydroxylated MgO(001) surface (Figure 1a). It is formed by dissociative adsorption of one water molecule per MgO surface ion pair H 2O + [Mg 2 +, O2 −]surf → [OH−Mg 2 +, H+O2 −]surf

(1)

Protons are placed on top of surface oxygen ions, OH− groups are positioned on top of the Mg2+ ions, and nondissociated water molecules are located above the hydroxylated surface layer. Although the multitude of surface OH groups is consistent with numerous experiments performed under ambient water conditions, in his 2006 review Ewing concludes, “What is not settled in the discussion of thin film water on MgO(001) is to what extent the oxide surface has been hydroxylated.”2 © 2015 American Chemical Society

Received: April 29, 2015 Accepted: June 4, 2015 Published: June 5, 2015 2310

DOI: 10.1021/acs.jpclett.5b00885 J. Phys. Chem. Lett. 2015, 6, 2310−2314

Letter

The Journal of Physical Chemistry Letters

found to leave the surface. This was interpreted as a possible first step of brucite (Mg(OH)2) layer formation on the MgO surface. Here we build on these results and investigate two types of surface structures by DFT: (i) Mg(OH)2-like structures and (ii) hydrated MgO(001) surface structures. We apply the PBE functional27,28 augmented with the empirical D2 dispersion term of Grimme (PBE+D).29 Adding the D2 term to PBE is not always an improvement, but this combination has met with success for one water monolayer on MgO(001),14 and it also yields an improved binding energy for an isolated H2O molecule on MgO(001), 60 kJ/mol compared with the observed 63 kJ/mol (SI, Table SI1); see the SI also for further information on the performance of PBE+D compared with PBE. Figure 1 shows the structures obtained and Table 1 shows the water coverage, the number of dissociated water molecules, the energy of formation from the pristine MgO(001) surface and gas-phase water molecules, and the respective surface related Gibbs free energies, whereas the chemical composition and surface area are given in the SI. For the Mg2+(O2−)k(OH−)l (H2O)m surface ions, the coordination numbers k, l, and m, with k + l + m = 6, are given in Table SI5 of the SI. All surface structures include four MgO layers, which correspond to the thickness of thin films studied experimentally.13,15 Mg(OH)2-Like Structures. Because of the exothermicity of the MgO(s) + H2O(g) → Mg(OH)2(s) reaction (−81 kJ/mol, Table SI1 in the SI), the formation of a four-layered Mg(OH)2(0001) surface slab (brucite, Figure 1b) represents the most stable hydration structure of the four-layered MgO surface slab at 0 K (Table 1); however, the experiments performed under ambient conditions show that the presence of the OH groups is limited to the surface.13,15 Moreover, brucite formation would require a deep reconstruction of the MgO(001) surface and is known to be slow requiring hours to days. (See ref 15 and references therein.) To create a single brucite layer on the MgO(001) surface (Figure 1c), one water monolayer is required. To account for the fact that the Mg(OH)2(0001) surface possesses trigonal cell symmetry, while the MgO(001) structure is cubic, we calculated the whole structure in the cubic cell and added the relaxation energy of the brucite layer (see SI). This, however, leads only to a marginal stabilization (Table 1). Hydrated MgO(001) Surface Structures. For low water coverages, 1 ML). In summary, for low OH signal levels, we confirmed previous conclusions14 that a partially dissociated water layer on the MgO(001) surface with p(3 × 2) or c(4 × 2) symmetry is most probably formed.14 For high OH-signal levels, structures with hydrated/hydroxylated Mg2+ ions above the surface explain the experimental data. OH groups are formed both in contact with the Mg2+ ions above the surface and on the surface as a part of the defect healing process (eq 2). Further reconstruction of the hydrated MgO(001) surface can eventually lead to a STRIPE structure. We have presented structure models for hydrated MgO(001) surfaces beyond monolayer coverage that explain observed XPS spectra and, according to our DFT calculations, are stable phases (minima on the potential energy surface). The latter is not true for previous models that assume a full monolayer of hydroxyl groups. In contrast, our models feature reconstructed surfaces with hydrated/hydroxylated Mg2+ ions above the surface layer. The understanding gained for the MgO(001) model surface will be useful when studying the hydration of more complex surfaces, such as LaCoO3 (001);31 however, more work is required to elucidate the role of the proposed structures for corrosion and dissolution processes. Specifically, future work should address the kinetics of the formation of the proposed structures, for example, the fact that the ability of the MgO surfaces to dissociate water is inhibited after several water adsorption−desorption cycles.7



ASSOCIATED CONTENT

S Supporting Information *

Details of computational methods, benchmarking calculations, interpretation of experiments, and Cartesian coordinates of included structures. The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpclett.5b00885.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work has been supported by the German Research Foundation (DFG) within CRC 1109 “Metal Oxide − Water Interfaces” and the Cluster of Excellence “Unifying Concepts in Catalysis”. MO thanks for a fellowship provided by the Alexander von Humboldt Foundation. The authors gratefully acknowledge the computing time granted by the John von Neumann Institute for Computing (NIC) and provided on the supercomputer JUROPA at Jülich Supercomputing Centre (JSC).



REFERENCES

(1) Verdaguer, A.; Sacha, G. M.; Bluhm, H.; Salmeron, M. Molecular Structure of Water at Interfaces: Wetting at the Nanometer Scale. Chem. Rev. 2006, 106, 1478−1510. (2) Ewing, G. E. Ambient Thin Film Water on Insulator Surfaces. Chem. Rev. 2006, 106, 1511−1526. (3) Woodruff, D. P. Quantitative Structural Studies Of Corundum and Rocksalt Oxide Surfaces. Chem. Rev. 2013, 113, 3863−3886. (4) Anderson, P. J.; Horlock, R. F.; Oliver, J. F. Interaction of Water with Magnesium Oxide Surface. Trans. Faraday Soc. 1965, 61, 2754− 2762. (5) Ferry, D.; Glebov, A.; Senz, V.; Suzanne, J.; Toennies, J. P.; Weiss, H. Observation of the Second Ordered Phase of Water on the MgO(100) Surface: Low Energy Electron Diffraction and Helium Atom Scattering Studies. J. Chem. Phys. 1996, 105, 1697−1701. (6) Liu, P.; Kendelewicz, T.; Gordon, G. E.; Parks, G. A. Reaction of Water with MgO(100) Surfaces. Part I: Synchrotron X-ray Photoemission Studies of Low-Defect Surfaces. Surf. Sci. 1998, 412−13, 287−314. (7) Abriou, D.; Jupille, J. Self-Inhibition of Water Dissociation on Magnesium Oxide Surfaces. Surf. Sci. 1999, 430, L527−L532. (8) Kim, Y. D.; Stultz, J.; Kim, Y. D. Dissociation of Water on MgO(100). J. Phys. Chem. B 2002, 106, 1515−1517. (9) Kim, Y. D.; Lynden-Bell, R. M.; Alavi, A.; Stulz, J.; Goodman, D. W. Evidence for Partial Dissociation of Water on Flat MgO(100) Surfaces. Chem. Phys. Lett. 2002, 352, 318−322. (10) Yu, Y. H.; Guo, Q. L.; Liu, S.; Wang, E. G.; Møller, P. J. Partial Dissociation of Water on a MgO(100) Film. Phys. Rev. B 2003, 68, 115414. (11) Savio, L.; Celasco, E.; Vattuone, L.; Rocca, M. Enhanced Reactivity at Metal-Oxide Interface: Water Interaction with MgO Ultrathin Films. J. Phys. Chem. B 2004, 108, 7771−7778. (12) Altieri, S.; Contri, S. F.; Valeri, S. Hydrolysis at MgO(100)/ Ag(100) Oxide-Metal Interfaces Studied by O 1s X-Ray Photoelectron and MgKL23L23 Auger Spectroscopy. Phys. Rev. B 2007, 76, 205413. (13) Carrasco, E.; Brown, M. A.; Sterrer, M.; Freund, H.-J.; Kwapien, K.; Sierka, M.; Sauer, J. Thickness-Dependent Hydroxylation of MgO(001) Thin Films. J. Phys. Chem. C 2010, 114, 18207−18214.



COMPUTATIONAL METHODS The DFT calculations employed the Vienna ab initio simulation package (VASP),32,33 with the projector-augmented wave method (PAW).34 The cutoff energy for plane waves was 600 eV, and the 2 × 2 × 1 Monkhorst−Pack k-point mesh was applied. The slab model was employed for the surface systems with a cell size of 30 Å in the z direction. For the surface cells, we used the lattice constant calculated for the bulk phase for both MgO and Mg(OH)2 at the same level of calculation with the 2 × 2 × 2 Monkhorst−Pack k-point mesh (see SI). For structure optimizations and frequency calculations, an energy convergence criterion of 10−6 eV was applied for both electronic and ionic steps. The two bottom layers were kept fixed to mimic their bulk nature. Molecular dynamics runs were performed with a cutoff energy for plane waves of 300 eV and an energy convergence criterion of 10−4 eV. All reported energies include zero point vibrational energies (ZPE). Gibbs free energies are calculated per surface area for selected 2313

DOI: 10.1021/acs.jpclett.5b00885 J. Phys. Chem. Lett. 2015, 6, 2310−2314

The Journal of Physical Chemistry Letters



(14) Wlodarczyk, R.; Sierka, M.; Kwapien, K.; Sauer, J.; Carrasco, E.; Aumer, A.; Gomes, J. F.; Sterrer, M.; Freund, H.-J. Structures of the Ordered Water Monolayer on MgO(001). J. Phys. Chem. C 2011, 115, 6764−6774. (15) Newberg, J. T.; Starr, D. E.; Yamamoto, S.; Kaya, S.; Kendelewicz, T.; Mysak, E. R.; Porsgaard, S.; Salmeron, M. B.; Brown, G. E., Jr.; Nisson, A.; et al. Autocatalytic Surface Hydroxylation of MgO(100) Terrace Sites Observed under Ambient Conditions. J. Phys. Chem. C 2011, 115, 12864−12872. (16) Newberg, J. T.; Starr, D. E.; Yamamoto, S.; Kaya, S.; Kendelewicz, T.; Mysak, E. R.; Porsgaard, S.; Salmeron, M. B.; Brown, G. E., Jr.; Nilsson, A.; et al. Formation of Hydroxyl and Water Layers on MgO Films Studied with Ambient Pressure XPS. Surf. Sci. 2011, 605, 89−94. (17) Carrasco, E.; Aumer, A.; Gomes, J. F.; Fujimori, Y.; Sterrer, M. Vibrational Spectroscopic Observation of Ice Dewetting on MgO(001). Chem. Commun. 2013, 49, 4355−4357. (18) Coulomb, J. P.; Demirdjian, B.; Ferry, D.; Trabelsi, M. Thermodynamic and Structural Properties of Water Adsorbed Film on MgO (100) Ionic Surface. Adsorption 2013, 19, 861−867. (19) Giordano, L.; Goniakowski, J.; Suzanne, J. Partial Dissociation of Water Molecules in the (3 × 2) Water Monolayer Deposited on the MgO (100) Surface. Phys. Rev. Lett. 1998, 81, 1271−1273. (20) Odelius, M. Mixed Molecular and Dissociative Water Adsorption on MgO 100. Phys. Rev. Lett. 1999, 82, 3919−3922. (21) Delle Site, L.; Alavi, A.; Lynden-Bell, R. M. The Structure and Spectroscopy of Monolayers of Water on MgO: An Ab Initio Study. J. Chem. Phys. 2000, 113, 3344−3350. (22) Cho, J. H.; Park, J. M.; Kim, K. S. Influence of Intermolecular Hydrogen Bonding on Water Dissociation at the MgO(001) Surface. Phys. Rev. B 2000, 62, 9981−9984. (23) Lynden-Bell, R. M.; Delle Site, L.; Alavi, A. Structures of Adsorbed Water Layers on MgO: an ab Initio Study. Surf. Sci. 2002, 496, L1−L6. (24) Jug, K.; Heidberg, B.; Bredow, T. Cyclic Cluster Study on the Formation of Brucite from Periclase and Water. J. Phys. Chem. C 2007, 111, 13103−13108. (25) Jug, K.; Heidberg, B.; Bredow, T. Cyclic Cluster Study of Water Adsorption Structures on the MgO(100) Surface. Surf. Sci. 2007, 601, 1529−1535. (26) Jug, K.; Heidberg, B.; Bredow, T. Molecular dynamics study of water adsorption structures on the MgO(100) surface. J. Phys. Chem. C 2007, 111, 6846−6851. (27) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (28) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple (vol 77, pg 3865, 1996). Phys. Rev. Lett. 1997, 78, 1396−1396. (29) Grimme, S. Semiempirical GGA-Type Density Functional Constructed with a Long-Range Dispersion Correction. J. Comput. Chem. 2006, 27, 1787−1799. (30) Robie, R. A.; Hemingway, B. S.; Fisher, J. R. Thermodynamic Properties of Minerals and Related Substances at 298.15 K and 1 bar Pressure and at Higher Temperatures; United States Government Printing Office: Washington, DC, 1978. (31) Stoerzinger, K. A.; Hong, W. T.; Crumlin, E. J.; Bluhm, H.; Biegalski, M. D.; Shao-Horn, Y. Water Reactivity on the LaCoO3 (001) Surface: An Ambient Pressure X-ray Photoelectron Spectroscopy Study. J. Phys. Chem. C 2014, 118, 19733−19741. (32) Kresse, G.; Furthmuller, J. Efficient Iterative Schemes for ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Phys. Rev. B 1996, 54, 11169−11186. (33) Kresse, G.; Furthmuller, J. Efficiency of ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15−50. (34) Blöchl, P. E. Projector Augmented-Wave Method. Phys. Rev. B 1994, 50, 17953−17979.

Letter

NOTE ADDED AFTER ASAP PUBLICATION This paper was published ASAP on June 5, 2015. Eq 1 and main text were updated. The revised paper was reposted on June 10, 2015.

2314

DOI: 10.1021/acs.jpclett.5b00885 J. Phys. Chem. Lett. 2015, 6, 2310−2314