Ab Initio Thermodynamics of Surface Oxide Structures under

Jan 10, 2017 - U. S. A. 2011, 108, 917– 924 DOI: 10.1073/pnas.1006669107 ...... A.; Alloyeau , D. New Insights into the Mixing of Gold and Copper in...
0 downloads 0 Views 1MB Size
Article pubs.acs.org/JPCC

Ab Initio Thermodynamics of Surface Oxide Structures under Controlled Growth Conditions Taehun Lee,† Yonghyuk Lee,† Simone Piccinin,‡ and Aloysius Soon*,† †

Department of Materials Science and Engineering, Yonsei University, Seoul 03722, South Korea Consiglio Nazionale delle Ricerche, Istituto Officina dei Materiali (CNR-IOM), c/o SISSA, Via Bonomea 265, 34136 Trieste, Italy



S Supporting Information *

ABSTRACT: Having a robust and predictive ab initio thermodynamic model to examine and describe the interplay of the oxygen gas and evaporated metal atoms on another metal substrate may prove to be very helpful in understanding the surface phase diagrams of these oxygen/metal systems. In this work, we examine the O/Cu/Au(111) system and provide a refined atomistic thermodynamic model which takes different definitions of the chemical potential of the less abundant metal, Cu into account. We argue that the latter highly depends on the various surface structures (overlayers and alloys) that forms on the metal substrate under growth conditions. We demonstrate that our improved thermodynamic model rationalizes new experimentally observed oxide structures and may pave a systematic way to predict new surface structures of reduced stoichiometries, which would otherwise be missed by the common practice of taking only the bulk limits.



INTRODUCTION Ultrathin oxide films grown on metal substrates show remarkably different structural properties from those of bulk oxides in the monolayers and the submonolayer regime, depending on the growth conditions.1−4 The ultrathin oxides on metal supports exhibit distinct physical and chemical properties, which can increase the performance of nanoscale devices, as well as representing an alternative model to study heterogeneous catalysis.5,6 The development of experimental spectroscopic and scanning probe techniques has led to vast improvements in the characterization of the films at the atomic scale.7,8 In line with the progress of nanoscale measurements, theory also has successfully provided insights to complement the experimental observations via atomistic modeling coupled with advanced thermodynamic approaches, namely ab initio atomistic thermodynamics (aiAT).9−12 These approaches make it possible to predict accurate local atomic structures and the thermodynamic stability of thin film oxides on metal and oxide surfaces as a function of temperature and pressure of reactive molecules species, by combining general concepts in thermodynamics with total energy calculations based on density-functional theory (DFT).13−17 These thin oxide layers are commonly prepared by two different methods: Either depositing metal atoms in an O2 atmosphere or evaporating metal atoms in ultrahigh vacuum (UHV) followed by an oxidation step.2−4,18 In these processes, the interfacial coupling between the thin oxide layer and the relevant metal supports has been commonly exploited to obtain the required morphologies and stoichiometries of the ultrathin film oxides, which are highly linked to their physical and chemical properties.3,18−20 Interestingly, it has been reported that the amount of evaporated reactive metal source © XXXX American Chemical Society

significantly influences the morphology and the stoichiometry of the grown oxide layer,4 but to date a quantitative description of these correlations is still lacking. Nanostructures of Cu2O/Au have caught a lot of interest due to their extraordinary performance in photocatalytic and photochemical applications, due to the fact that they display photoluminescence properties that differ significantly from that of Cu2O bulk.21,22 Via X-ray photoelectron spectroscopy (XPS) and scanning tunneling microscopy (STM) experiments at 400 and 700 K, various Cu thin film oxides were observed and characterized.18 Low energy thin film oxides were modeled using DFT, displaying simulated STM images in good agreement with experiments.16 DFT calculations conventionally consider the chemical reservoir of Cu originated from its bulk phase. In our study, the deposited Cu atoms may first form only a few monolayers on the Au surface, deviating from our conventional assumption of using bulk Cu as the chemical reservoir. According to the experiments, the deposited copper indeed forms a few monolayers of metallic Cu on Au surface.18 In this step, the film thickness was initially adjusted to about 10 to 15 ML and it then decreased to about 5 ML upon annealing at 700 K. In the case of the CeO2/Cu(111) system,4 it was revealed that the final surface structure highly depends on the initial deposition conditions of cerium. Likewise in the Cu/ Au(111) system, nanocomposites such as overlayers or surface alloys can form depending the deposition conditions of Cu. However, it is still an open question which chemical potential Received: November 14, 2016 Revised: January 10, 2017 Published: January 10, 2017 A

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

chosen to be Cu bulk phase.13,16,17 In this assumption, it is regarded that the amount of evaporated Cu atoms is enough to reach equilibrium condition to the bulk phase (μCu ≃ Ebulk Cu ). However, under the experimental conditions where the amount of atomic source of Cu is minute compared to that of Au and controlled as a certain amount to prepare the oxide layer in the monolayer regime,18 the reservoir of Cu cannot be simply assumed as the bulk Cu (μCu ≠ Ebulk Cu ). When a minute amount of Cu is deposited in an O2 ambience, the evaporated Cu atoms can react with substrate and locally form nanocomposite alloys, CuxAuy/Au(111) (where x ≪ y). Experiments monitoring the formation of thin Cu oxide layers after Cu evaporation on Au(111) have detected the appearance of few monolayers of Cu on Au(111) as well as AuxCuy alloys prior to the formation of the oxide.18,31 In modeling these systems, we therefore need to consider that the sources of Cu atoms forming the oxide are these nanocomposite structures (overlayers and alloys) formed on Au(111). Since the chemical potential of Cu atoms in these structures will be different depending on their structure and composition, we need to investigate the impact this has on the thermodynamic stability of the surface oxides. Our approach implicitly assumes that the amount of Cu deposited on the Au substrate is very large compared to the amount of Cu that ends up forming the thin oxide layers, so that the chemical potential of Cu is determined by structure adopted by Au/Cu system after the deposition of Cu. In particular, in previous theoretical and experimental studies, various nanocomposite (CuxAuy) structures have been reported as thermodynamic products under different conditions, showing structural diversities from overlayer to heterostructures.32−35 In this spirit, we examine and consider various Cu−Au bimetallic structures in this work. In a schematic view, we propose the various chemical origins of Cu from different AuxCuy/Au(111) nanocomposite systems to form thin oxide layers in Figure 1. Specifically, we present a scheme displaying various AuxCuy/Au(111) nanocomposite systems that can act as a potential source of Cu atoms for the oxide layer. The thin oxide formation reactions could be examined and approximated as a two-step process, i.e. the initial formation of these mildly doped bimetallic alloy surfaces (cf. eq 2), and then

values do nanocomposites have and how they affect the energetic stability in the thin film oxides. Here, we use a refined aiAT approach to search for the most favorable atomic configurations adopted by Cu−O thin oxides on Au(111). We assume that Cu atoms deposited on Au(111) act as a reservoir of Cu for the formation of the surface thin oxides. Our results show that different definitions of the Cu chemical potential can have a dramatic effect on the stability of metal supported oxides. We therefore investigate the dependence of the Cu chemical potential on the Cu coverage as well as on the atomic structure of the Cu/Au system. Furthermore, we propose the ideal amount of Cu to deposit for obtaining a specific oxidic structure.



METHODOLOGY All calculations are performed using Kohn−Sham DFT and the generalized gradient approximation (GGA) to the exchangecorrelation functional due to Perdew, Burke, and Ernzerhof (PBE),23 as implemented in the Vienna ab initio simulation package (VASP).24,25 The Kohn−Sham orbitals are expanded in a planewave basis set with a kinetic energy cutoff of up to 500 eV and the electron−ion interactions are described using the projector augmented wave (PAW) method.26 All atomic calculations have been performed with spin polarization. For the surface structures, we have also tested and found that they all possess a nonmagnetic ground-state solution. Ultrathin films of O/Cu on Au(111) and bimetallic system of Cu on Au(111) are modeled using asymmetric periodic supercell models. In addition, following ref 15, to obtain a more accurate description of EO2, the experimental binding energy of an O2 molecule (5.22 eV)27 is used to correct for the large known error in the PBE binding energy of O2. As estimated in our earlier work,28,29 the vibrational contributions to the Gibbs free energy of adsorption in O/Cu systems are found to be negligible and is thus not included in this work. In addition, via the Stirling’s approximation,30 we have estimated the contribution of configurational entropy to our calculated Gibbs free energies of adsorption, and have found it to be less than 25 meV/Å2 for T < 700 K. According to this estimate, the inclusion of configurational entropy will not change the conclusions drawn here, and thus neglected in this work.



RESULTS AND DISCUSSION When a metal (Cu) source is deposited on different metal surface (Au surface) in contact with an oxygen atmosphere, the gas phase acts as a reservoir of oxygen, interchanging oxygen atoms with the metal surface. The Gibbs free energy of adsorption can be defined by ΔG(T , p) =

1 (GO/Cu/Au − GAu(111) − ΔNAuμAu A − ΔNCuμCu − NOμO)

(1)

where GO/Cu/Au and GAu(111) are the Gibbs free energies of the O/Cu/Au system and pristine Au(111), respectively. A is the surface area of unit cell and, μAu, μO, and μCu are the chemical potentials of the Au, O and Cu atoms, accordingly. ΔNAu and ΔNCu are the change of the number of species to involve in forming oxide layer and substrate. NO is the number of adsorbed O atoms. When ΔNAu is zero, the Gibbs free energy can be interpreted as a function of μCu and μO. The chemical reservoir of Cu, i.e. the source of Cu atoms forming the ultrathin oxide layer, has been conventionally

Figure 1. Schematic view of the formation of CuxAuy nanocomposites and thin Cu oxide layers on Au(111), reflecting the various values possible for the Cu chemical potential. The O, Cu, and Au atoms are depicted as red, blue, and gray circles, respectively. B

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Having established how the value of the chemical potential of Cu changes depending on the structures formed by the Cu/Au system upon evaporation of Cu on Au(111), we can now explore the stability of the thin oxide structures formed once the system is exposed to an oxygen atmosphere. In the case of oxygen, the chemical potential of oxygen is defined with respect to the half of the total energy of the O2 molecule, i.e. 1 ΔμO = μO − 2 EO2 . Notably, the thin oxide layer is not only in equilibrium with the O2 gas phase, but also in equilibrium with the underlying specific CuxAuy nanocomposite on Au(111). Here, we can write the surface free energy of adsorption as

the oxidation of these surfaces (cf. eq 3). First, the evaporated Cu atoms lead to the formation of Cu−Au nanocomposites: xCu atom + y Au bulk + Au(111) → CuxAu y /Au(111)

(2)

This Cu−Au nanocomposite system acts as the thermodynamic reservoir which provides Cu atoms to form the Cu oxide layer on Au as shown in eq 3: CuxAu y /Au(111) + O2 (g) ⇌ O/Cux /Au(111)

(3)

where x ≪ y. The Gibbs free energy of CuxAuy/Au(111) system can be written as GCuxAu y = NxμCu + NyμAu

(4)

ΔG(ΔμO , ΔμCu ) =

where Nx and Ny are the number of atoms of Cu and Au in the CuxAuy nanocomposite, respectively. Assuming thermodynamic equilibrium between the CuxAuy nanoalloy and the underlying Au(111) surface, the chemical potential of Au, μAu, can be taken equal to the chemical potential of bulk Au, Ebulk Au . We now define that the change of Cu chemical potential with respect to the atom energy of Cu, ΔμCu = μCu − Eatom Cu . Then eq 4 can be reformulated as 1 atom ΔμCu = (GCuxAu y − NyμAu − NxECu ) Nx (5)

1⎡ ⎢GO/Cu/Au − GAu(111) A⎣ atom ) − NCu(ΔμCu + ECu

⎛ 1 ⎞⎤ − NO⎜ΔμO + EO2⎟⎥ ⎝ 2 ⎠⎦

where ΔμCu =

1 (GCuxAu y Nx

(6)

bulk atom − NyEAu − NxECu ). Via eq 6, we

plot a 2D surface phase diagram as a continuous function as shown Figure 3, illustrating accessible surface compositions in Figure 3. The surface structures we investigate in this work have been proposed in our previous publication,16 matching with experiment.18 Here, we briefly outline their structural character-

We computed ΔμCu for the various Cu−Au nanocomposite systems as shown in Figure 2a (and listed in Table S1 of the

Figure 2. Illustrations of various CuxAuy nanocomposites on Au(111) are shown in part a, where the Cu and Au atoms are depicted as blue and gray circles accordingly. The corresponding changes in the Cu chemical potential (ΔμCu) are plotted in part b, with respect to the Cu surface coverage (ΘCu). Figure 3. PBE-calculated surface phase diagram for thin Cu surface oxide layers on Au(111): The change in the Gibbs free energy of adsorption of oxygen plotted as a function of both the change in the chemical potentials of oxygen (ΔμO) and copper (ΔμCu). The topand side-views of surface atomic structures, p2, p2 + CuF, p2s, p2s + CuF, and p2s + CuF−OL on Au(111) are shown. The Au, Cu, and O atoms are depicted as white, blue, and red circles, respectively. ΔμO = 1 0 corresponds to μO = 2 EO2 , as determined from the experimentally corrected DFT binding energy of the O2 molecule.15,27 The pressure scales for 400 and 700 K are included at the bottom of the figure. To calculate the ΔG in Figure 4, ΔμCu is defined as (a) bulk Cu, (b) ΘCu = 2.63 ML, and (c) ΘCu = 0.33 ML (as indicated by the labeled horizontal dotted lines).

Supporting Information). In here, μAu is taken as Ebulk Au , since Au(111) is close to the infinite amount. We have tested various nanocomposite models on Au(111), controlling the surface coverage of Cu (ΘCu) on Au(111) substrate (see Figure 2a). The computed values of ΔμCu have similar tendencies for various structural models. First, ΔμCu values saturate to that of bulk Cu when ΘCu is larger than about 1.00 ML. Second, if ΘCu is less than 1.00 ML, ΔμCu values become positive, compared to that of bulk Cu. Later, we discuss more about the impact of these values of ΔμCu on the stability of the thin oxide films. C

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C

Figure 4. Calculated Gibbs free energy of adsorption, ΔG for various thin film oxides of Cu on Au(111) as a function of the change in oxygen chemical potential at specific ΔμCu values: (a) − 3.52 eV which corresponds to bulk Cu; (b) − 3.12 eV which corresponds to ΘCu = 2.63 ML, as an example of monolayer CuxAuy structures; (c) −0.53 eV which then corresponds to ΘCu = 0.33 ML, as an example of overlayer CuxAuy structures. At the bottom of the figure, surface phase transition points are shown to reflect the most stable surface structures at the corresponding pO2 for T = 700 K.

that the thin oxide film on substrate cannot exist when the values of ΔμO are larger than the bulk oxide’s heat of formation per O atom, due to the bulk oxide being thermodynamically more stable than thin oxide films.13 Of the possible bulk oxides that can form, the enthalpy of formation per O atom of Cu2O is f = −1.67 eV. However, here we the most negative, ΔHCu 2O argue that the upper limit of the oxygen chemical potential may be better described by extending the physically meaningful range for the O chemical potential to the formation of bulk Au oxide. Thus, the oxygen rich limit of adsorption energy can be 1 f , where ΔHfAu2O3 (−1.77 eV) defined as ΔμO − limit ≃ 3 ΔH Au 2O3 is the enthalpy of formation of Au2O3. These thermodynamic limits are noted in Figure 4. Notwithstanding the importance of thermodynamic analysis in this work, we stress that kinetic effects may also play a role and could alter the thermodynamically controlled surface phase diagram of O/Cu/Au(111). This has been discussed in a previous study where bulk oxide phases were kinetically hindered, leading to the formation of thin metal oxide layers under oxygen-rich conditions, even up to ambient pressures.36 Now, to further assess the thermodynamic stability of the thin Cu oxide layers, ΔG is computed when chemical reservoir of Cu is defined as (a) bulk, (b) ΘCu = 2.63 ML, and (c) ΘCu = 0.33 ML, as shown in Figure 4. In the case where ΔμCu is set to that of bulk Cu, the pristine Au(111) surface is thermodynamically most stable surface phase for the values of ΔμO less than −1.85 eV. In addition, p2 and p2s surface structures are found to be the most stable across the range of ΔμO from −1.85 to −0.59 eV. After that, bulk Au2O3 may start to form as the thermodynamic product. When ΔμCu is extracted from Cu monolayers model (ΘCu = 2.63 ML) on Au(111), both the p2s + CuF − OL (−2.54 to −2.33 eV) and p2s (−2.33 to −0.59 eV) are found to be stable

istics: The p2 surface structures have a distorted p(1 × 1) Cu2O(111)-like trilayer on p(2 × 2) Au(111). p2 + CuF builds on this p2 motif by adding extra Cu adatoms at the fcc-hollow site on Au(111) in the center of the Cu−O hexagonal network. These p2-family structures have been proposed as the thermodynamic product in previous experimental reports.18 Next, the p2s structure adds an extra row/stripe of CuL in p2 by artificially breaking the surface symmetry, and p2s + CuF has additional Cu adatoms at the center of Cu−O hexagonal network. In addition, the p2s + CuF − OL structure is formed by extracting the nearest O atoms along the Cu adatoms line (CuL) from p2s + CuF. For a more detailed atomic analysis of Cu oxide layer on Au(111), we refer the reader to our previous publication.16 When ΔμCu is set to the bulk Cu value (marked in Figure 3c), p2s structure is the most stable over a range of ΔμO. As ΔμCu (above −3.52 eV) becomes more positive, the stable surfaces contain an increasing amount of Cu atoms, transitioning from p2s to p2s + CuF and p2s + CuF − OL via the adsorption of Cu atoms at the center of the Cu−O hexagonal ring. As discussed earlier, we can link the specific ΔμCu values with corresponding Cu−Au nanocomposites, which is directly connected to the ΘCu, as shown in Figure 2. For further analysis of the relationship between ΔμCu and the stability of oxide layer, we choose two representative conditions for the surface underneath the thin oxide layers, namely Cu monolayers (ΘCu = 2.63 ML, ΔμCu = −3.12 eV) and Cu overlayers (ΘCu = 0.33 ML, ΔμCu = −0.53 eV) (marked using horizontal dotted lines in Figure 3 in parts b and c), and plot ΔG for various Cu oxide structures on Au(111) (see Figure 4) comparing with that when ΔμCu assumes the bulk value. To establish a physically meaningful upper limit to the range of values of the oxygen chemical potential, we need to consider D

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C Notes

in series, shown in Figure 4b. In this case, the difference of ΔμCu with respect to that of bulk is 0.40 eV. This small change of ΔμCu results in the appearance of new stable structures such as p2s + CuF − OL under oxygen lean condition. As the small amount of Cu is adsorbed on Au(111) in the overlayer form (ΘCu = 0.33 ML), it brings significant variations to the stable thin oxide layer configurations as shown in Figure 4c. The p2s + CuF − OL and p2s + CuF oxide structures are now found to be the most stable, which may rationalize the possible observation of the 3-fold CuF atom in the center of the hexagonal ring as suggested in ref 18. These oxide layers show a structural transition at ΔμO = −1.97 eV. These results demonstrate the clear correlation between evaporated amounts of Cu and stable oxide layer structures. Referring to Figure 3, for a fixed value of ΔμO, larger values of ΔμCu will stabilize thin oxide structures with a higher Cu content. This concept can be applied to actual Cu oxide layer preparation, providing the tunability of stoichiometries in these oxide layers. Particularly based on our results, if the deposited amount of evaporated atoms forms Cu−Au nanocomposites below about 1 ML, then the Cu atoms on hexagonal networks can be stabilized, leading reduced stoichiometries of thin CuxO layers.

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We gratefully acknowledge support from the National Research and Development Program by the NRF (Grant No. 2016M3A7B4024140). Computational resources have been provided by the KISTI supercomputing center (KSC-2016-C30038) and the Australian National Computational Infrastructure (NCI).





CONCLUSION In conclusion, we have refined the atomistic thermodynamic model for describing the interaction of oxygen with a bimetallic system where one metal is significantly minute as compared to the other metal (here, Cu ≪ Au). We have argued that deviating from the bulk limit of the deposited metal (i.e., Cu in this study), its chemical potential can take discrete values highly depending on the various surface structures (overlayers and alloys) formed on the metal substrate (i.e., Au(111) in this work). In an attempt to rationalize the experimentally observed O/Cu/Au(111) structures, we have then demonstrated that taking a different definition of the Cu chemical potential dramatically affects the stability of thin Cu surface oxide structures on Au(111). In particular, with respect to the controversy of the existence of the 3-fold Cu adatoms in the center of the hexagonal network, the p2s + CuF oxide structure can indeed be stabilized if the amount of deposited Cu atoms forms Cu−Au nanocomposites below a Cu surface coverage of about 1 ML. This approach of analyzing the thermodynamic surface phase diagram may pave a new way to discover novel surface oxide structures of reduced stoichiometries, which would otherwise be missed by commonly taking only the bulk limits.



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.6b11445. Computed Cu chemical potential (ΔμCu) values for various models of CuxAuy/Au(111) nanocomposite system, Gibbs free energy of adsorption, and ΔG for various thin film oxides of Cu on Au(111) (PDF)



REFERENCES

(1) Freysoldt, C.; Rinke, P.; Scheffler, M. Ultrathin Oxides: BulkOxide-Like Model Surfaces or Unique Films? Phys. Rev. Lett. 2007, 99, 086101. (2) Staudt, T.; Lykhach, Y.; Hammer, L.; Schneider, M. A.; Matolín, V.; Libuda, J. A Route to Continuous Ultra-thin Cerium Oxide Films on Cu(111). Surf. Sci. 2009, 603, 3382−3388. (3) Szabová, L.; Stetsovych, O.; Dvořaḱ , F.; Farnesi Camellone, M.; Fabris, S.; Mysliveček, J.; Matolín, V. Distinct Physicochemical Properties of the First Ceria Monolayer on Cu(111). J. Phys. Chem. C 2012, 116, 6677−6684. (4) Vári, G.; Ó vári, L.; Kiss, J.; Kónya, Z. LEIS and XPS Investigation into the Growth of Cerium and Cerium Dioxide on Cu(111). Phys. Chem. Chem. Phys. 2015, 17, 5124−5132. (5) Rodriguez, J. A.; Liu, P.; Graciani, J.; Senanayake, S. D.; Grinter, D. C.; Stacchiola, D.; Hrbek, J.; Fernández-Sanz, J. Inverse Oxide/ Metal Catalysts in Fundamental Studies and Practical Applications: A Perspective of Recent Developments. J. Phys. Chem. Lett. 2016, 7, 2627−2639. (6) Freund, H.-J.; Pacchioni, G. Oxide Ultra-thin Films on Metals: New Materials for the Design of Supported Metal Catalysts. Chem. Soc. Rev. 2008, 37, 2224−2242. (7) Nilius, N. Exploring Routes to Tailor the Physical and Chemical Properties of Oxides via Doping: an STM Study. J. Phys.: Condens. Matter 2015, 27, 303001. (8) Somorjai, G. A.; Li, Y. Impact of Surface Chemistry. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 917−924. (9) Reuter, K.; Scheffler, M. Composition, Structure, and Stability of RuO2 (110) as a Function of Oxygen Pressure. Phys. Rev. B: Condens. Matter Mater. Phys. 2001, 65, 035406. (10) Reuter, K.; Scheffler, M. First-principles Atomistic Thermodynamics for Oxidation Catalysis: Surface Phase Diagrams and Catalytically Interesting Regions. Phys. Rev. Lett. 2003, 90, 046103. (11) Stampfl, C. Surface Processes and Phase Transitions from Ab Initio Atomistic Thermodynamics and Statistical Mechanics. Catal. Today 2005, 105, 17−35. (12) Stampfl, C.; Soon, A.; Piccinin, S.; Shi, H.; Zhang, H. Bridging the Temperature and Pressure Gaps: Close-packed Transition Metal Surfaces in an Oxygen Environment. J. Phys.: Condens. Matter 2008, 20, 184021. (13) Soon, A.; Todorova, M.; Delley, B.; Stampfl, C. Oxygen Adsorption and Stability of Surface Oxides on Cu(111): A Firstprinciples Investigation. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 73, 165424. (14) Shi, H.; Stampfl, C. First-principles Investigations of the Structure and Stability of Oxygen Adsorption and Surface Oxide Formation at Au(111). Phys. Rev. B: Condens. Matter Mater. Phys. 2007, 76, 075327. (15) Richter, N. A.; Kim, C.-E.; Stampfl, C.; Soon, A. Re-visiting the O/Cu(111) SystemĐ-When Metastable Surface Oxides Could Become an Issue! Phys. Chem. Chem. Phys. 2014, 16, 26735−26740. (16) Lee, T.; Lee, Y.; Kang, K.; Soon, A. In Search of Nonconventional Surface Oxidic Motifs of Cu on Au(111). Phys. Chem. Chem. Phys. 2016, 18, 7349−7358. (17) Piccinin, S.; Stampfl, C.; Scheffler, M. AgĐCu Alloy Surfaces in an Oxidizing Environment: A First-principles Study. Surf. Sci. 2009, 603, 1467−1475.

AUTHOR INFORMATION

Corresponding Author

*(A.S.) E-mail: [email protected]. Telephone: +82-22123-5839. ORCID

Aloysius Soon: 0000-0002-6273-9324 E

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX

Article

The Journal of Physical Chemistry C (18) Sträter, H.; Fedderwitz, H.; Groß, B.; Nilius, N. Growth and Surface Properties of Cuprous Oxide Films on Au(111). J. Phys. Chem. C 2015, 119, 5975−5981. (19) Lunkenbein, T.; Schumann, J.; Behrens, M.; Schlögl, R.; Willinger, M. G. Formation of a ZnO Overlayer in Industrial Cu/ ZnO/Al?2O3 Catalysts Induced by Strong Metal-Support Interactions. Angew. Chem. 2015, 127, 4627−4631. (20) Fedderwitz, H.; Groß, B.; Sträter, H.; Nilius, N. Reactive Copper Deposition on Au(111) and Mo(001): Role of the Support in the Oxidation Process. J. Phys. Chem. C 2016, 120, 7591−7596. (21) Liu, Y. L.; Liu, Y. C.; Mu, R.; Yang, H.; Shao, C. L.; Zhang, J. Y.; Lu, Y. M.; Shen, D. Z.; Fan, X. W. The Structural and Optical Properties of Cu2O Films Electrodeposited on Different Substrates. Semicond. Sci. Technol. 2005, 20, 44−49. (22) Cushing, S. K.; Li, J.; Meng, F.; Senty, T. R.; Suri, S.; Zhi, M.; Li, M.; Bristow, A. D.; Wu, N. Photocatalytic Activity Enhanced by Plasmonic Resonant Energy Transfer from Metal to Semiconductor. J. Am. Chem. Soc. 2012, 134, 15033. (23) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865. (24) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Liquid Metals. Phys. Rev. B: Condens. Matter Mater. Phys. 1993, 47, 558. (25) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for ab initio Total-energy Calculations Using a Plane-wave Basis Set. Phys. Rev. B: Condens. Matter Mater. Phys. 1996, 54, 11169. (26) Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials to the Projector Augmented-wave Method. Phys. Rev. B: Condens. Matter Mater. Phys. 1999, 59, 1758. (27) Feller, D.; Peterson, K. A. Re-examination of Atomization Energies for the Gaussian-2 Set of Molecules. J. Chem. Phys. 1999, 110, 8384−8396. (28) Soon, A.; Todorova, M.; Delley, B.; Stampfl, C. Surface Oxides of the Oxygen-copper System: Precursors to the Bulk Oxide Phase? Surf. Sci. 2007, 601, 5809−5813. (29) Bohnen, K.-P.; Heid, R.; Pintschovius, L.; Soon, A.; Stampfl, C. Ab Initio Lattice Dynamics and Thermal Expansion of Cu2O. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 134304. (30) Reuter, K.; Scheffler, M. Composition and Structure of the RuO2(110) Surface in An O2 and CO Environment: Implications for The Catalytic Formation of CO2. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 68, 045407. (31) Tsuda, Y.; Oka, K.; Makino, T.; Okada, M.; Dino, W. A.; Hashinokuchi, M.; Yoshigoe, A.; Teraoka, Y.; Kasai, H. Initial Stages of Cu3Au(111) Oxidation: Oxygen Induced Cu Segregation and the Protective Au Layer Profile. Phys. Chem. Chem. Phys. 2014, 16, 3815− 3822. (32) Magnussen, O. M.; Hotlos, J.; Nichols, R. J.; Kolb, D. M.; Behm, R. J. Atomic Structure of Cu Adlayers on Au(100) and Au(111) Electrodes Observed by In Situ Scanning Tunneling Microscopy. Phys. Rev. Lett. 1990, 64, 2929. (33) Macur, J. E.; Vook, R. W. Initial Epitaxial Growth of (111)Au/ (111)Cu and (111)Cu/(111)Au. Thin Solid Films 1980, 66, 371−379. (34) Prunier, H.; Nelayah, J.; Ricolleau, C.; Wang, G.; Nowak, S.; Lamic-Humblot, A.; Alloyeau, D. New Insights into the Mixing of Gold and Copper in a Nanoparticle from a Structural Study of AuĐCu Nanoalloys Synthesized via a Wet Chemistry Method and Pulsed Laser Deposition. Phys. Chem. Chem. Phys. 2015, 17, 28339−28346. (35) Mettela, G.; Kulkarni, G. U. Site Selective Cu Deposition on Au Microcrystallites: Corners, Edges versus, Planar surfaces. CrystEngComm 2015, 17, 9459−9465. (36) Lundgren, E.; Gustafson, J.; Mikkelsen, A.; Andersen, J. N.; Stierle, A.; Dosch, H.; Todorova, M.; Rogal, J.; Reuter, K.; Scheffler, M. Kinetic Hindrance during the Initial Oxidation of Pd(100) at Ambient Pressures. Phys. Rev. Lett. 2004, 92, 046101.

F

DOI: 10.1021/acs.jpcc.6b11445 J. Phys. Chem. C XXXX, XXX, XXX−XXX