High-Coverage H2 Adsorption on the ... - ACS Publications

Sep 26, 2017 - The probability distribution of the reconstructed Cu2O(111) surfaces as ... the lattice oxygen to the surface at higher coverage (five ...
0 downloads 0 Views 5MB Size
Subscriber access provided by Purdue University Libraries

Article 2

2

High Coverage H Adsorption on the Reconstructed CuO(111) Surface Xiaohu Yu, Xuemei Zhang, Hongtao Wang, Zhiyin Wang, and Gang Feng J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b06361 • Publication Date (Web): 26 Sep 2017 Downloaded from http://pubs.acs.org on September 27, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

High Coverage H2 Adsorption on the Reconstructed Cu2O(111) Surface Xiaohu Yua*, Xuemei Zhanga,b, Hongtao Wangb*, Zhiyin Wanga and Gang Fengc a) Institute of Theoretical and Computational Chemistry, Shaanxi Key Laboratory of Catalysis, School of Chemical & Environment Sciences, Shaanxi University of Technology, Hanzhong 723000, P.R. China; b) College of Environmental Science and Engineering, Taiyuan University of Technology,Taiyuan 030024, P.R. China; c) College of Chemistry, Nanchang University, Nanchang, Jiangxi 330031, P.R. China. Xiaohu Yu Tel: +86 18739700918; E-mail: [email protected] Hongtao Wang Tel: +86 13834607727; Email: [email protected] RECEIVED DATE (xxxx) Abstract The adsorption of H2 on the Cu2O(111) surface has been studied by spin-polarized density functional theory (DFT+U) calculations and atomic thermodynamics. It has been found that there exists reconstruction on stoichiometric Cu2O(111) surface. Probability distribution of the reconstructed Cu2O(111) surfaces as function of temperature has been analyzed using Boltzman statistics. It has been found that the molecular H2 prefers to adsorption on the CuCUS atom at low coverages (1/4 or 1/2 monolayer), while totally dissociative H2 is preferred on reconstructed Cu2O(111) surface at higher coverages (3/4 or 1 monolayer). For H2 splitting on Cu2O(111) surface, homolytical dissociative adsorption on two surface un-coordinated CuCUS atoms is preferred which is new mechanism for H2 on metal oxides. More interesting is that the surface reconstruction will be recovered for eight hydrogen atoms binding on four un-coordinated CuCUS and four un-coordinated OCUS atoms at saturation coverage. It has been found that the adsorbed H atoms will put out the lattice oxygen to the surface at higher coverage (five and six H2), which agrees well with the experimental findings. The phase diagrams of H2 binding on ideal and reconstructed Cu2O(111) surfaces were plotted and analyzed. In addition, we compared and analyzed the adsorption mechanisms of H2 splitting on different metal oxides.

ACS Paragon Plus Environment

1

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 28

1. Introduction The study of the dissociation of the hydrogen molecule at interfaces is an active field that is motivated by the desire to understand many processes, such as heterogeneous catalysis, hydrogen storage, fuels, sensors, doping of semiconducting oxides, etc.1 It is well known that there are two mechanisms of H2 splitting on metal oxide surfaces: heterolytic or homolytic mechanisms, depending on the nature of the metal oxide 2,3. Cuprous oxide (Cu2O)4 is an attractive p-type oxide for photoelectrochemical hydrogen production with a direct bandgap of 2 eV. However, interaction of H2 on Cu2O surfaces is still not clear theoretically. The interaction of H2 with copper oxides attracted many attentions. Althrough the reduction of metal oxide using hydrogen is a frequently used method to prepare active catalysts and electronic devices, the reduction process is generally not well understood at a molecular or atomic level. Using Time-resolved X-ray diffraction (XRD), X-ray absorption fine structure (XAFS), and first-principles density functional calculations, Kim et al.3,5,6 studied the reduction of CuO and Cu2O with H2 and they found that the mechanism for the reduction of CuO is complex, involving an induction period and the embedding of H into the bulk of the oxide. In addition, they found that the reduction of CuO is easier than the reduction of Cu2O and the apparent activation energy for reduction of CuO is 14.5 kcal/mol, while the value is 27.4 kcal/mol for Cu2O. Tabuchi et al. 7 studied control of carrier concentration in thin cuprous oxide Cu2O films by atomic hydrogen, and found that atomic H exposure is a useful tool for improving the properties of Cu2O. Poulston et al.8 studied surface oxidation and reduction of CuO and Cu2O using XPS and XAES, and found that hydrogen-reduced CuO and Cu2O were both reoxidized on vacuum annealing, demonstrating the diffusion of lattice oxygen to the surface. Using ultraviolet photoelectron spectroscopy (UPS) and thermal desorption spectroscopy (TDS) , Cox and Schulz9 studied water adsorption on the Cu-terminated Cu2O(100) surface and found that adsorption at 110 K is both dissociative and molecular, with only about 10% of a monolayer of water dissociated, while adsorption at 300 K is dissociative only. In addition, they10 studied the interaction of atomic hydrogen with Cu2O(100) and found that surface hydride formation is formed on the metal oxide surface. Nygren et al.11 studied H2O interaction with the polar Cu2O(100) surface and found that dissociative adsorption is favoured than molecular adsorption. Li et al.12 studied H2O adsorption and dissociation on Cu2O(100) and found that activation barrier for HxO(x=1, 2) dehydrogenation are 0.42 eV and 1.86 eV. Using ACS Paragon Plus Environment

2

Page 3 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DFT+U, Maimaiti et al.13 studied the reduction mechanisms of the CuO(111) surface through surface oxygen vacancy formation and hydrogen adsorption, and found that the CuO surface is reduced to Cu2O at approximately 360 K and that complete reduction from Cu2O to metallic Cu occurs at 780 K using H2 as the reducing agent. Recently, water adsorption on CuO surfaces at different coverage was studied by us14,15. Using high-resolution photoelectron spectroscopy (PES) and Hubbard U and dispersion corrected density functional theory (PBE-D3+U) calculations, Stenlid et al.16 studied water interactions and reactions on a common Cu2O(100):Cu surface and found that the results are consistent for ambient temperatures under wet/humic and oxygen lean conditions. Cu2O(111) surface attracted extensive researches experimentally and theoretically as main grown surface of Cu2O surfaces17 and having higher catalytic activity in many surface processes18,19. Using scanning tunneling microscopy (STM) and low-energy electron diffraction, Sträter et al.20 studied the growth and surface properties of cuprous oxide films on Au(111) surface and found that the structure of postannealed films is similar to the one of bulk-cut Cu2O(111). Using STM and hybrid density functional theory, Nilius et al.21 studied the stability of Cu2O(111) surface and found that the stoichiometric surface is the most stable at oxygen-lean conditions which is confirmed by an excellent matching between STM spectroscopy data and the calculated surface electronic structure. It is reported that the octahedral Cu2O with exposed (111) facets exhibited much higher photocatalytic activity than cubes

22,23

. Zheng et al.24 researched the surface stabilities and photocatalytic properties of Cu2O

microcrystals and found that Cu2O exposing (111) facets can be as a stable photocatalyst. Önsten et al. 25

studied role of defects in suface chemistry on Cu2O(111), and found that there are two different

defects on the surface. Using cluster model, Casarin et al.26 studied H2O and H2S adsorption on Cu2O(111) cluster and found that molecular H2O chemisorption is favoured than the dissociative adsorption on Cu2O cluster. Recently, surface properties of Cu2O surface27, high coverage water adsorption on Cu2O(111) surface28 and O2 adsorption Cu2O surfaces29 were studied by us. CO2 adsorption on Cu2O(111) surface was researched by different groups

30-33

and water interaction with

Cu2O(111) surface was also studied by Riplinger and Carter 34. To shed light on the interaction of H2 on Cu2O(111) surface and unraveling the role of different coverages of H2 adsorbates, systematically studying the different H2 adsorption is very important. To the best of our knowledge, there is only a few reports about the H2 adsorption on the Cu2O surfaces have ACS Paragon Plus Environment

3

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 28

been reported35,36. Zhang et al.35 studied one single H2 adsorption on Cu2O(111) p(2×2) surface with six-layer model and found that surface un-coordinated CuCUS is the most advantageous position with the side-on type of H2. Li et al.

36

researched single H2 adsorption on the different Cu2O(111) p(1×1)

surfaces and found that parallelly H2 adsorbs on surface un-coordinated CuCUS atom. In the present paper, we have performed a systemtically DFT+U study of different coverages of H2 on reconstructed Cu2O(111) surface. The probability distribution of different reconstructed Cu2O(111) surfaces as function of temperature are calculated using Boltzmann statistics. Atomic thermodynamics is used to analyze the adsorption properties of H2 on ideal and reconstructed Cu2O(111) surfaces. The overall trends of the adsorption energy as a function of coverage are discussed and the adsorption properties of H2 on Cu2O(111) surfaces and other metal oxide surfaces are compared. 2. Methodology and Surface Model 2.1. Computational Methodology The exchange and correlation energies were calculated using the Perdew, Burke and Ernzerhof (PBE) functional.37 The projector augmented wave (PAW) method described the electron-ion interaction.38 The well-known GGA+U method, Ueff

= 6 eV for Cu2O which is same to former literature

34

, was

introduced for the Cu-3d electrons to describe the on-site Coulomb interaction. The Kohn-Sham oneelectron states were expanded in a plane wave basis set up to 500 eV for bulk, and smaller 350 eV for slab in order to decrease the computational cost. The integrations over the Brillouin zone were performed using Monkhorst-Pack grids39 and the conjugate gradient algorithm was applied to determine the electronic ground state. The results presented in this work were obtained using k-point meshes of 5×5×5 for the bulk and 3×3×1 for the surface calculations. A Gaussian broadening of the Fermi surface of 0.2 eV was applied to improve the convergence of the solutions. In all calculations the positions of all atoms were fully relaxed until the Hellmann-Feynman forces were smaller than 0.02 eV/Å and the energy was smaller than 10-4 eV. The dipole moment induced by adsorbed H molecules was corrected by dipole moment which is perpendicular to the surface.40,41 All calculations were performed using spin polarized plane-wave periodic density functional method as implemented in the Vienna ab initio simulation package (VASP).

40,41

The DFT-D3 method of Grimme42,43, which takes into account the

dispersive interactions, was used to test the adsorption of H2 on the Cu2O(111) surface. It was found that the adsorption energies of the most stable molecular (-0.54 eV), homolytical (-0.34 eV) and heterolytical ACS Paragon Plus Environment

4

Page 5 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(-0.31 eV) H2 adsorptions on the Cu2O(111) surface using (GGA+U)-D3 method are larger than corresponding molecular (-0.47 eV), homolytical (-0.30 eV) and heterolytical (-0.28 eV) adsorptions using GGA+U. However, the vdw force correction does not affect our conclusion. So, the results reported in this work use the GGA+U method. For the location of transition states, the climbing image nudged elastic band (CI-NEB) algorithm was employed44. The nature of the transition states was assessed by performing numerical frequency analysis. The reported Bader charges were calculated using the charge densities obtained with VASP and the program developed by Henkelman et al. 45 2.2. Model and Properties Calculated. It is well known that the Cu2O crystal has a simple cubic unit cell with (Pn3) structure and there are six atoms in Cu2O unit cell: 2O atoms form a bcc lattice and each O atom is at a center of a tetrahedron, whose four vertices are occupied by Cu atoms. Using Ueff = 6 eV, we derived of the unit cell parameters of 4.285 Å which is in reasonable agreement with the experimental value of about 4.270 Å and former theoretical result (4.270 27 with DFT and 4.317 Å46 with DFT+U). For the Cu2O(111) surface, there are three different terminations: O-terminated surface (2.526 J/m2), Cu-terminated surface (2.300 J/m2) and Cu/O terminated Cu2O(111) surface (0.784 J/m2)27,29. Among three terminations, Cu/O terminated Cu2O(111) surface was chosen in this work since it is most stable termination17 and has higher activity in many reactions47,48. In experiment, there are reported that two experimentally characterized structures of the (111) surface: the ideal Cu/O terminated surface with (1× 1) periodicity and the surface with (√3×√3)R30°periodicity and 2/3 of the terminating oxygen ML.49 The most stable perfect Cu/O terminated Cu2O(111) surface is non-polar surface. It is reported that there are six-layered slabs to simulate H235, O250, H2S51, N2O52 adsorption on Cu2O(111) surface, nine-layered slabs to simulate CO oxidation53 and CO2 electrochemical reduction32 on Cu2O(111) surface, twelve-layer slabs to simulate CO adsorption54 on Cu2O(111) surface and fifteen-layer slabs to simulate CO231 and water34 binding on Cu2O(111) surface. Considering computational accuracy and efficient, a nine-layered slabs made of (2×2) supercell was used to represent Cu2O(111) surface in this work (Figure 1). A 15 Å vacuum region was used to separate slabs in order to minimize the spurious interactions between the periodic images. In the whole calculating process, a (12.204 Å×12.204 Å× 21.228 Å) hexagonal supercell was used to simulate the Cu2O(111) surface. It exposed four different types of surface atoms in Cu2O(111) surface17, which were denoted as OCUS (the outmost surface oxygen) ACS Paragon Plus Environment

5

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

in the first layer, OCSA (the subsurface oxygen) in the third layer, CuCUS (coordinative unsaturated copper) in the second layer, and CuCSA (coordinative saturated copper) in the second layer. The most stable molecular H2 adsorption on two nine-layered models was tested: the adsorption energy on nine-layered slab with fixed bottom three-layer is larger than on nine-layered slab with fixed bottom six-layer about 0.03 eV. So in all following reported calculations, the atoms of the top three layers and adsorbed H atoms are fully relaxed, while the atoms of the bottom six layers are fixed to their bulk positions. All possible adsorption sites and structures are considerred for each adsorbate, and only those with the largest adsorption energies are used for discussion and comparison; the less stable adsorption sites and structures are given in the Supporting Information. The total adsorption energies were calculated using equation: ΔE(H2)n = E(nH2 /slab)- E(slab)-nE(H2) where E(nH2/slab) is total energy of the adsorbed system, E(slab) is the energy of slab, and E(H2) is the energy of H2 in the gas phase. Gas-phase H2 molecule was calculated in a box of (10×10×10) Å3 using the same cutoff parameters as for the slab calculation. The stepwise adsorption enery was defined as: ∆∆E(H2)n =∆E(H2)n -∆E(H2)n-1, and the average adsorption energy was defined as: ΔE =ΔE(H2)n/n. The bond length of molecular H2 calculated from our approach is 0.750 Å, the bond energy is about 6.09 eV, which is in well agreement with the experimental values of 0.740 Å 36

55

and other theoretical results

. Above results obtained in these tests made us confident in pursuing the next step of our researchs,

namely the interaction of H2 with Cu2O surfaces. (Table 1) and (Figure 1) 3. Results and Discussion 3.1. Reconstructed Cu2O(111) surface: First, we studied the reconstructed Cu2O(111) surface and the most stable reconstructed Cu2O(111) surface is shown in Figure 1a. The optimized structure of perfect Cu2O(111) surface is shown in Figure 1b and there are two different oxygen sites: four un-coordinated oxygen atoms (OCUS) in the first layer, four coordinated oxygen atoms (OCSA) in the third layer. The CuCUS-O bond is 1.848 Å, which is much shorter than the bulk Cu-O bond length of 1.868 Å. In addition, there are four surface un-coordinated CuCUS atoms in the second layer on the Cu2O(111) surface as shown in Figure 1b, A, B, C, D represent four surface un-coordinated CuCUS atoms. In order to shed light on the reconstructed Cu2O(111) surface, ACS Paragon Plus Environment

6

Page 7 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

we move the surface un-coordinated CuCUS atom in 1% about 0.122 Å in x and y axis, respectively, before fully relaxation. So the AxyBCD means move the A atom in x and y about 0.122 Å, and AxyBxyCD means move the A and B atoms in x and y axis about 0.122 Å, respectively. In total, eight different reconstructed Cu2O(111) surface models were considered in this work. After fully relaxation, the AxyBxCD model was found to be the most stable reconstructed Cu2O(111) surfaces as shown in Figure 1a, which is more stable than ideal Cu2O(111) surface about 0.764 eV. It indicates that the fully relaxed ideal Cu2O(111) surface is metastable surface which is in well agreement with former theoretical results56,57. The four surface un-coordinated CuCUS atoms move to one side of reconstructed Cu2O(111) surface (Figure 1a), the distances of CuCUS-CuCSA change from 3.580 to 2.446, 2.547, 2.649, 2.446, 2.547, 2.649, 2.643, 2.450, 2.555, 2.643, 2.450, and 2.555 Å, respectively. The extra stability of reconstructed Cu2O(111) surface comes from the shortened CuCUS-CuCSA distances. The configurations of least stable reconstructed Cu2O(111) surfaces were shown in Supporting Information Figure S1. In addition, the probability distribution as function of temperature was analyzed according the Boltzmann statistics (Figure 2). Each reconstructed Cu2O(111) surface has a probability of occurrence according Boltzmann statistics, Pm, which is function of the temperature: Pm = 1/Zexp[-∆E/(KBT)], where Z is the canonical partition function, and ∆E is the free energy of reconstructed and ideal Cu2O(111) surfaces. The AxyBxCD model of the reconstructed Cu2O(111) (Figure 2a) is the most preferred at 0 K. Along the much higher temperature, probability distribution of different least stable reconstructed Cu2O(111) surfaces increase. One can clearly see that probability distribution of the different least stable reconstructed Cu2O(111) surfaces (Figure 2c) increase at 1000 K, while the ideal Cu2O(111) surface only have very small probability. It indicates that different reconstructed Cu2O(111) surfaces can exchangeable at higher temperature and the reconstructed Cu2O(111) surface may play an important role in the surface chemistry of Cu2O(111) surface. (Figure 2) 3.2. H2 molecules on reconstructed Cu2O surface: Same to the ideal Cu2O(111) surface, there exposed four surface un-coordinated OCUS atoms in first layer, four surface un-coordinated CuCUS atoms in the second layer, 12 surface coordinated CuCSA atoms in the second layer and four coordinated OCSA in the third layer on the reconstructed Cu2O(111) surface. Here the coverage is defined as the ratio of the number of adsorbed H2 molecules to the number of ACS Paragon Plus Environment

7

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 28

surface un-coordinated CuCUS atoms. In the following, if we do not particular stress, the total adsorption energy is corresponding to the reconstructed surface. (Figure 3) For a single H2 molecule at a 1/4 monolayer (ML) coverage, molecular, heterolytic and homolytic adsorptions were computed. In the most stable molecular adsorption (Figure 3a), molecular H2 parallelly adsorbs on surface un-coordinated CuCUS atom with H-CuCUS bond lengths of 1.662 and 1.651 Å, and the bond length of H-H changes from 0.750 (free H2 molecule) to 0.825 Å: the total adsorption energy is -0.47 eV. Our adsorption mode agrees well with the former theoretical result 36. The adsorption energy (-0.73 eV)36 of the former result is much larger than our result (-0.47 eV). There should be two reasons: one is that they did not consider the surface reconstruction in former work, the other is that the slab that they used is p(1×1) in former work and we used is slab p(2×2). In the most stable heterolytic adsorption (Figure 3b), one hydrogen atom adsorbs on surface un-coordinated CuCUS atom with H-CuCUS bond length of 1.488 Å, and another hydrogen atom binds on the surface un-coordinated OCUS with HOCUS bond length of 0.988 Å: the total adsorption energy is -0.28 eV. In the most stable homolytic adsorption (Figure 3c), one hydrogen atom adsorbs on one surface un-coordinated CuCUS atom with HCuCUS bond length of 1.490 Å, and another hydrogen goes to another surface un-coordinated CuCUS with H-CuCUS bond length of 1.483 Å; the total adsorption energy is -0.30 eV. Comparing with H2 adsorption on two surface CuCUS atoms via homolytical dissociation, H2 adsorbs on two surface OCUS atoms via homolytical dissociation with total adsorption energy of 0.38 eV (Figure S2f). One can clearly see that there exists an unique mechanism for H2 splitting on Cu2O(111) surface which is different with former mechanisms of H2 splitting on other metal oxides: H2 adsorbs on two surface un-coordinated CuCUS atoms via homolytical dissociation. It should be mentioned that at the initial stage put two hydrogen atoms on surface un-coordinated OCUS atom, the surface oxygen atom will be pulled out and forms a H2O molecule with total adsorption energy of -0.11 eV (Figure S2b). It indicates that the water molecule formation on reconstructed Cu2O(111) surface by two hydrogen atoms binding with surface uncoordinated OCUS atom is favorable thermodynamically. On the basis of the most stable molecular and heterolytic dissociative adsorption (Figure 3a and Figure 3b), we computated dissociation transition state44, in which the breaking H-H distance is elongated from 0.825 to 1.160 Å and formed H-O distance is 1.198 Å: the computated heterolytic dissociation barrier is 1.22 eV (Figure S2j). In the same time, the ACS Paragon Plus Environment

8

Page 9 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

H atom which binds on surface CuCUS derives 0.35 e by Bader analysis45, while another H atom which binds on surface OCUS loses 0.46 e. On the basis of the most stable heterolytic dissociative adsorption and homolytic dissociative adsorption (Figure 3b and Figure 3c), the H-OCUS distance is 1.216 Å and the H-CuCUS distance is 1.805 Å in the transition state (Figure S2k), and the diffusion barrier of H atom is 0.62 eV. Comparing three adsorption modes, one can easily see that molecular adsorption is the most preferred on reconstructed Cu2O(111) surface at 1/4 ML coverage and H2 adsorption on two CuCUS atoms via homolytical dissociation is new mechamism for H2 splitting. For two H2 molecules at a 2/4 ML coverage, molecular, mixed and dissociative adsorptions were computed. In the most stable molecular H2 adsorption configuration (Figure 3d), two H2 molecules adsorb on two surface un-coordinated CuCUS atoms with four H-CuCUS bond lengths of 1.658, 1.660, 1.659 and 1.650 Å, and two H-H bond lengths change from 0.750 Å (free H2 molecule) to 0.822 and 0.825 Å; the total adsorption energy is -0.87 eV. In the most stable mixed molecular and dissociative adsorption configuration (Figure 3e), one molecular H2 binds on one surface un-coordinated CuCUS atom with two H-CuCUS bond lengths of 1.646 and 1.656 Å, and the bond length of H-H changes from 0.750 to 0.825 Å, and another H2 heterolytically dissociates with one hydrogen atom binding on surface uncoordinated CuCUS atom with H-CuCUS bond length of 1.490 Å, and another H atom goes to surface uncoordinated OCUS atom with H-OCUS bond length of 0.990 Å; the total adsorption energy is -0.66 eV. In the most stable dissociative adsorption configuration (Figure 3f), two hydrogen atoms bind on two surface un-coordinated CuCUS atoms with both H-CuCUS bond lengths of 1.492 Å, other two hydrogen atoms go to two surface un-coordinated OCUS atoms with two H-OCUS bond lengths of 0.989 and 0.991 Å; the total adsorption energy is -0.84 eV. Comparing three adsorption modes, one can clearly see that the molecular adsorption is the most preferred for two H2 molecules at 2/4 ML coverage. For three H2 molecules at a 3/4 ML coverage, molecular, mixed and dissociative adsorptions were computed. In the most stable molecular adsorption configuration (Figure 3g), three molecular H2 parallelly bind on three surface un-coordinated CuCUS atoms with six H-CuCUS distances of 1.655, 1.645, 1.650, 1.644, 1.655 and 1.649 Å, and three H-H bond lengths change from 0.750 to 0.825, 0.827, and 0.826 Å: the total adsorption energy is -1.25 eV. In the most stable mixed molecular and dissociative adsorption configuration (Figure 3h), two molecular H2 parallelly adsorb on the surface un-coordinated CuCUS atoms with four H-CuCUS bond lengths of 1.649, 1.643, 1.654, and 1.646 Å, and the H-H bond ACS Paragon Plus Environment

9

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 28

lengths change from 0.750 Å to 0.827 and 0.828 Å, and another H2 heterolytically dissociates with one H atom binding surface un-coordinated CuCUS atom with H-CuCUS distance of 1.491 Å and another H atom binding on surface un-coordinated OCUS atom with H-OCUS distance of 0.990 Å: the total adsorption energy is -1.21 eV. In the most stable dissociative adsorption configuration (Figure 3i), three hydrogen atoms adsorb on three surface un-coordinated OCUS atoms and one dissociated hydrogen atom binds on one surface un-coordinated CuCUS atom, and other two dissociated hydrogen atoms go to bridge sites of surface un-coordinated CuCUS and surface coordinated CuCSA atoms; the total adsorption energy is -1.34 eV. Comparing three adsorption modes, it is clearly see that the totally dissociative adsorption is most preferred on reconstructed Cu2O(111) surface at 3/4 ML coverage. For four H2 molecules at a 1 ML coverage, molecular, mixed and dissociative adsorptions were computed. In the most stable molecular adsorption configuration (Figure 3j), four molecular H2 adsorb on four surface un-coordinated CuCUS atoms with total adsorption energy of -1.77 eV. In the most stable mixed molecular and dissociative adsorption configuration (Figure 3k), one molecular H2 adsorbs on one surface un-coordinated CuCUS atom, and other three H2 heterolytically dissociate to six hydrogen atoms: three hydrogen atoms bind on three un-coordinated CuCUS atoms and other three hydrogen atoms go to three surface un-coordinated OCUS atoms; the total adsorption energy is -1.73 eV. In the most stable dissociative adsorption configuration (Figure 3l), four hydrogen atoms adsorb on four surface uncoordinated CuCUS atoms with four H-CuCUS bond lengths of 1.490 Å and other four hydrogen atoms bind on four surface un-coordinated OCUS atoms with four H-Osurf distances of 0.898 Å: the total adsorption energy is -1.98 eV. More interesting is that four H2 adsorption on reconstructed Cu2O(111) surface via heterolytical dissociation can get rid of the surface reconstruction. One can see that totally dissociative adsorption is most preferable for four H2 molecules at 1 ML coverage in which the adsorption of four H2 molecules completely saturates the surface cationic and anionic sites. For five H2 molecules at a 5/4 ML coverage, the molecular, mixed and dissociative adsorptions were computed. No molecular configuration was located for five H2 molecules. In the most stable mixed molecular and dissociative coadsorption configuration (Figure 3m), one can easily see that there are four dissociated H atoms adsorption on surface un-coordinated CuCUS atoms and other four hydrogen atoms binding on four surface un-coordinated OCUS atoms, one molecular H2 physisorbs on surface uncoordinated CuCUS and the distance of H-CuCUS is 3.501 Å: the total adsorption energy is -1.94 eV. Here ACS Paragon Plus Environment

10

Page 11 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

it should be mentioned that in the most stable chemisorption configuration (Figure S5c), four molecular H2 adsorb on four surface un-coordinated CuCUS atoms and two dissociated H atoms adsorb on two surface un-coordinated OCUS atoms with the total adsorption energy of -0.61 eV. In the most stable dissociative adsorption configuration (Figure 3n), four H atoms adsorb on four surface un-coordinated OCUS atoms, and six H atoms bind on the top of surface un-coordinated CuCUS and bridge sites of CuCUS and CuCSA atoms with the total adsorption energy of -0.75 eV. More interesting is that at the beginning, put two H atoms on two surface coordinated OCSA atoms in the third layer (Figure S6a), the two OCSA atoms will move out to surface layer by formed two H-OCSA bonds and induce the surface strong reconstruction with the total adsorption energy of -3.16 eV, which is in well agreement with experimental report8: the diffusion of lattice oxygen atom to surface by H adsorption. Comparing the four H2 and five H2 adsorption on Cu2O(111) surface, one can clearly see that five H2 adsorption on Cu2O(111) is not favorable thermodynamically. For six H2 molecules at a 6/4 ML coverage, the molecular, mixed and dissociative adsorptions were computed. No configuration of six molecular H2 adsorption on the Cu2O(111) surface was located. In the most stable mixed molecular and dissociative coadsorption configuration (Figure 3o), four H atoms favor binding on the four surface un-coordinated CuCUS atoms and other four H atoms favor adsorption on the four surface un-coordinated OCUS atoms, and last two molecular H2 physisorb on the Cu2O(111) surface as indicated the long H2-CuCUS distance of more than 3.5 Å: the total adsorption energy is -1.92 eV. It should be mentioned that four molecular H2 chemisorb on surface un-coordinated CuCUS atoms (Figure S7b), and other four H atoms adsorb on surface un-coordinated OCUS atom with the total adsorption energy of 0.67 eV. In the most stable dissociative adsorption configuration (Figure 3p), four H atoms adsorb on four surface un-coordinated OCUS atoms and other dissociative H atoms bind at bridge sites of CuCUS and CuCSA atoms with the total adsorption energy of -0.85 eV. More interesting is that initially put four hydrogen atoms on the four coordinated OCSA atoms in the third layer (Figure S7c), and four OCSA atoms in the third layer will move out of surface and induce the strong reconstruction with the total adsorption energy of 0.22 eV, which agrees well with experimental report8: the diffusion of lattice oxygen atom to surface by H adsorption. One can easily see that six H2 adsorption on Cu2O(111) surface at 6/4 ML coverage is not favorable thermodynamically. The stepwise adsorption energy for most stable one to six H2 molecules is -0.47, -0.40, -0.37, -0.65, ACS Paragon Plus Environment

11

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 28

0.04, and 0.02 eV, and corresponding average adsorption energy for one to six H2 molecules is -0.47, 0.44, -0.45, -0.50, -0.35, -0.39, and -0.32 eV, respectively. This indicates that the adsorption of four H2 molecules on the reconstructed Cu2O(111) surface is favorable thermodynamically since the negative value of the stepwise adsorption energy. It is found that there is new mechanism for H2 splitting on Cu2O(111) surface: H2 adsorbs on two surface CuCSA atoms via homolytic dissociation. There is no clear trend about the coverage effect: along higher H2 coverage, both stewpwise and average adsorption energies first increase and then decrease, which is mainly due to the surface reconstruction and Coulomb repulsion between adsorbed H atoms. 3.3. Phase diagram of H2 adsorption at given conditions: To shed light on the stable adsorption structure of H2 on the Cu2O(111) surface under experimental conditions, the effect of different temperature and pressure was considerred using atomic thermodynamics 58,59. It is well known that atomic thermodynamics is very useful tool to study the effect of environment on the bulk60, surface61,62 and cluster63-65. According atomic thermodynamics, we chosed the Gibbs free energy (∆G) of nH2 adsorption on the Cu2O(111) surface as the criterion and defined it as: ads ∆GCu (T , P, nH 2 ) = G[nH 2 /Cu 2 O(111)] - G[Cu 2 O(111)] - nGH 2 (T , pH2 ) 2O

(1)

where G[nH 2 /Cu 2 O(111)] , G[Cu 2 O(111)] and GH 2 (T , pH2 ) are the Gibbs free energies of adsorbed system, support and adsorbates, respectively, n is number of adsorbed H2 molecules, T is the temperature, pH2 is the partial pressure of H2 in the gas atmosphere. Here a more negative ∆G indicates the more stable adsorption structure. We can write GH 2 (T , pH2 ) as:

GH2 (T, pH2 ) = EHtotal + µH2 (T, p0 ) + kΒT ln( p / p0 ) 2

(2)

where EHtotal is the total enery of H2 molecules including zero point vibration energy, µH2(T, p0) term 2 includes vibrational and rotational contributions for H2 gas, and can be taken from tables of thermodynamic data 66. kΒT ln( p / p 0 ) is the contribution of temperature and H2 partial pressure to the H2 chemical potential and kΒ is the Boltzmann constant. Since frequency motions of surface-slab have large entropy contribution to the surface, we included the vibration contribution of Cu2O(111) surface with adsorbed H2 in the DFT calculated total energy to substitute the Gibbs free energies of solid surface. Then, we can rewrite equation (1) as: ads ∆GCu (T , P, nH 2 ) = G[nH 2 /Cu 2 O(111)] - Ε[Cu 2 O(111)] - nGH2 (T , pH2 ) 2O

(3)

One can see that there are three stable surfaces: clean Cu2O(111), one H2 adsorption and four H2 ACS Paragon Plus Environment

12

Page 13 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

adsorption on ideal Cu2O(111) surface, while there are two stable surfaces: reconstructed Cu2O(111), four H2 adsorption on reconstructed Cu2O(111) surface, as shown in convex hulls of Figure 4a and 4c, One can recast Figure 4a and 4c in another insightful forms, as shown in the phase diagram (Figure 4b and 4d). Inserting equation (2) into (3), one can derive: ads ∆GCu (T , P, nH 2 ) = G[nH 2 /Cu 2 O(111)] - Ε[Cu 2 O(111)] - nEHtotal - nµH2 (T , p 0 ) − nkΒT ln( p / p 0 ) (4) 2O 2

On the basis of H2 adsorption, atomic thermodynamics is applied to study the effects of temperature and pressure on its coverage. The phase diagram of H2 on ideal and reconstructed Cu2O(111) surface at wide pressure and temperature ranges is shown in Figure 5a and 5b. One can see that there are three stable adsorption configurations on ideal Cu2O(111) surface: clean Cu2O(111) surface, one H2 molcules adsorption and four H2 molecules adsorption, while there are two stable adsorption configurations on reconstructed Cu2O(111) surface: reconstructed Cu2O(111) surface and four H2 molecules adsorption. One can see that H2 molecules will desorb on the Cu2O(111) surface at high temperature and low H2 pressure from the diagram which agrees well with physical intuition. Two and three H2 molecules binding on Cu2O(111) surface which is not shown on the two phase diagrams maybe that the process is controlled by kinetically not by thermodynamically. The stable single H2 binding on ideal Cu2O(111) surface shown on phase diagram can be understood from Cu2O(111) surface reconstruction which induced by the H2 adsorption. Under UHV conditions full desorption of the adsorbed H2 molecules takes place in the temperature range of 150-200 K from phase diagram (Figure 5a) of H2 binding on reconstructed Cu2O(111), while full desorption of the adsorbed H2 molecules takes place in the temperature range of 250-300 K from phase diagram (Figure 5b) of H2 binding on ideal Cu2O(111) surface. Overall, the 2H and 8H terminations on ideal Cu2O(111) surface are stabilized in the pressure ranges accessible in UHV experiments, while only 8H termination on reconstructed Cu2O(111) surface is stabilized in the pressure ranges accessible in UHV experiments. (Figure 4) and (Figure 5)

3.4. Discussion: To derive a full picture of the H2 interaction with reconstructed Cu2O(111) surface, we discuss here the energetic trends as a function of H2 coverage. The total adsorption energies for different coverages are plotted in Figure 4 and also displayed in Table 1. Comparing three adsorption modes, parallelly ACS Paragon Plus Environment

13

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 28

molecular H2 adsorption on surface un-coordinated CuCUS atom is preferred for an isolated H2 on reconstructed Cu2O(111) surface. It indicates that surface un-coordinated CuCUS is the most advantageous position with the parallelly adsorbed H2 which agrees well with former theoretical result35,36. From above data, one can easily see that molecular H2 parallelly adsorption on surface uncoordinated CuCUS atoms with side-on mode is most preferred for one and two H2 molecules, while totally heterolytic dissociative adsorption is most preferred for three and four H2 molecules on reconstructed Cu2O(111) surface. H2 adsorption on reconstructed Cu2O(111) surface is not favorable thermodynamically for five and six H2 molecules. It is interesting to compare the H2 splitting mechanisms on Cu2O(111) surface at different coverage. For one H2 molecule, homolytic (-0.30 eV) and heterolytic (-0.28 eV) dissociation is competitive on reconstructed Cu2O(111) surface; for two H2 molecules, heterolytic dissociation (-0.84 eV) is much more preferred than homolytic dissociation (0.54 eV Figure S3i); for three and four H2 molecules, heterolytic dissociation (-1.34 eV and -1.98 eV) is most preferred on reconstructed Cu2O(111) surface. For much higher H2 coverage (five or six H2 molecules), the coordinatively saturated OCSA atoms in the third layer (lattice oxygen) will diffuse to surface of Cu2O(111), which is in well agreement with experimental report8: the diffusion of lattice oxygen atom to surface by H adsorption. (Table 2) It is also interesting to compare H2 adsorption mechanisms on different metal oxide surfaces (Table 2). For the most stable H2 adsorption, there is no correlation between the surface stability and adsorption energy. For Fe3O4(110), the more stable A layer (0.947 J/m2)67 has homolytic dissociative adsorption energy of -1.73 eV68, and the less stable B layer (1.020 J/m2)67 has homolytic dissociative adsorption energy of -0.78 eV68. For Fe3O4(111), the most stable Fetet1 termination (0.944 J/m2)67 has homolytic dissociative adsorption energy of -1.62 eV69, and the second stable Feoct2 termination (1.132J/m2)67 has heterolytic dissociative adsorption energy of -1.21 eV69. On the different Fe3O4 directions, the H2 adsorption rule is reverse with CO adsorption70: much more stable terminations with large H2 adsorption energy, while much more stable terminations with small CO adsorption energy. For Fe3O4(001), H2 prefers homolytical dissociation with adsorption energy of -0.88 eV71. For CeO2, the more stable (111) surface has homolytic dissociative adsorption of -2.80 eV72, and less stable (110) surface has homolytic dissociative adsorption energy of -3.57 eV72. One can easily see that usually H2 prefers homolytic ACS Paragon Plus Environment

14

Page 15 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

dissociative adsorption on reducible transition metal oxide surface, which agrees well with experimental results. In addition, H2 adsorption mechanism on reducible metal oxide surfaces is in fact often controlled by specific properties of metal oxide surfaces, including electronic and structural surface defects. For example, H2 prefers to heterolytical dissociation on Feoct2 termination of Fe3O4(111) 69. On the contrary, H2 prefers to heterolytical dissociation on nonreducible oxides. For MgO(001), H2 prefers to heterolytical dissociation on terrace with adsorption energy of 1.80 eV73, and H2 prefers to heterolytical dissociation on step with adsorption energy of -0.43 eV73. H2 prefers to hereolytical dissociation on γ-Al2O3(110) surface with adsorption energy of -0.97 eV74. ZrO2 is a poor reducible oxide, H2 prefers to heterolytical dissociation on ZrO2(111) with adsorption energy of 0.24 eV75, and H2 prefers to heterolytical dissociation on ZrO2(101) with adsorption energy of -0.08 eV75. On Cu2O(111) surface, molecular H2 adsorption on surface un-coordinated Cuus atom from side-on mode is most preferred for one or two H2 molecules, while totally heterolytically dissociative adsorption on Cu2O(111) surface is most preferred for three and four H2 molecules. Comparing H2 dissociative adsorption on Cu2O(111) surface for one H2 molecule, H2 adsorption on two un-coordinated CuCUS atoms via homolytical dissociation is much more favorable. It indicates that there is unique mechanism for H2 dissociative adsorption on Cu2O(111) surface which is different with former H2 splitting mechanisms. Table 2 also shows the adsorption energies of H2 dissociative adsorption on different metal oxides, but these results are not comparable. 4. Conclusion Using spin-polarized density functional theory calculations (DFT+U) and atomic thermodynamics we systematically study H2 adsorption on the Cu2O(111) surface at different coverage. It has been found that there exists surface reconstruction of Cu/O terminated Cu2O(111) surface. Boltzmann statistics was used to analyze the probability distribution of different reconstructed Cu2O(111) surfaces as function of temperature. It has been found that molecular H2 adsorption on reconstructed Cu2O(111) surface is favorable at low coverage (1/4 and 2/4 ML), while totally dissociative adsorption is preferred at high coverage (3/4 and 1 ML). At higher H2 coverage (5/4 and 6/4 ML), the lattice oxygen atoms will diffuse to surface induced by H adsorption, which is in well agreement with experimental results. Four H2 molecules adsorption on reconstructed Cu2O(111) surface is favorable thermodynamically according to the stepwise adsorption energy. In addition, four totally dissociative adsorbed H2 molecules will get rid ACS Paragon Plus Environment

15

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 28

of the surface reconstruction. The phase diagram of H2 molecules on ideal and reconstructed Cu2O(111) surfaces were plotted to analyze the different H2 adsorption mechanisms. A deeper understanding of the reaction mechanism of H2 adsorption on Cu2O(111) surfaces has been obtained. More interesting is that four H2 molecules adsorption on four surface un-coordinated CuCUS atoms and four surface un-coordinated OCUS atoms via heterolytical dissociation is saturation adsorption. The adsorption mechanisms of H2 on different metal oxide surfaces are compared. It has been found that there exists an unique mechanism for H2 splitting on Cu2O(111) surface which is different with former mechanisms. The different H2 adsorption mechanisms on different metal oxides may benefit for new H2splitting metal oxides catalyst designing. The carefully designed experiments are also suggested to be conducted to verify the phenomenon.

Supporting Information: The computed less stable reconstructed Cu2O(111) surfaces, less stable configurations of one, two, three, four, five and six H2 adsorption on Cu2O(111) surface. This material is available free of charge via the Internet at http://xxx/xxx

ACKNOWLEDGMENTS This work was supported by National Natural Science Foundation of China (No. 21673270) and the National Natural Science Foundation of Henan (No. 162300410001).

ACS Paragon Plus Environment

16

Page 17 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 1. Top and side views of (a) the reconstructed and (b) the ideal Cu2O(111) surface. The energy of the ground state is set as the reference: A more positive value denotes a less stable ideal surface. (first layer O atoms in bigger red ball; second layer coordinatively un-saturated Cu atoms in green ball; other Cu and O atoms in blue and red ball)

ACS Paragon Plus Environment

17

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 28

Figure 2. Probability distribution of different reconstructed Cu2O(111) surfaces as function of temperature.

ACS Paragon Plus Environment

18

Page 19 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3. Configurations of the most preferred (a)-(c) 1H2; (d)-(f) 2H2; (g)-(i) 3H2; (j)-(l) 4H2; (m)-(n) 5H2; (o)-(p) 6H2 adsorption on the Cu2O(111) surface at different coverage (first layer O atoms in bigger red ball; second layer coordinatively un-saturated Cu atoms in green ball; other Cu and O atoms in blue and red ball , adsorbed H atoms in smaller white ball)

ACS Paragon Plus Environment

19

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 28

Figure 4. Phase diagram of stable H2 adsorption as a function of (a) number of H2 molecules and (b) H2 chemical potential on reconstructed Cu2O(111) surface; (c) number of H2 molecules and (d) H2 chemical potential on ideal Cu2O(111) surface. (Solid lines show the most stable adsorption configurations at different coverage)

ACS Paragon Plus Environment

20

Page 21 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. Phase diagram of stable H2 adsorption on (a) reconstructed and (b) ideal Cu2O(111) surfaces under different conditions.

ACS Paragon Plus Environment

21

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 28

Table 1. Total adsorption energy (eV) of H2 on reconstructed Cu2O(111) surface at different coverage. 1 H2

2H2

3H2

4H2

molecular

-0.47

-0.87

-1.25

-1.77

mixed

-0.26

-0.66

-1.21

dissociative

-0.30

-0.84

-1.34

Adsorption

5H2

6H2

-1.73

-0.75

-0.85

-1.98

-0.61

0.67

states Cu2O(111)

ACS Paragon Plus Environment

22

Page 23 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 2. Comparison of the largest adsorption energies (eV) for H2 on different surfaces of different metal oxides at low coverage

a

adsorption states

site

PBE

PW91

GGA+U

Fe3O4(110) A-layer

homolytic

O-top

-1.7368

Fe3O4(110) B layer

homolytic

O-top

-0.7868

Fe3O4(001)

homolytic

O-top

-0.8871

Fe3O4(111) Fetet1

homolytic

O-top

-1.6269

Fe3O4(111) Feoct2

heterolytic

Fe, O

-1.2169

CeO2(111)

homolytic

O-top

-0.1872

-2.4076, -2.8072

CeO2(110)

homolytic

O-top

-1.1672

-3.5772

ZrO2(111)

heterolytic

Zr, O

0.2475

ZrO2(101)

heterolytic

Zr, O

-0.0875

MgO(001) terrace

heterolytic

Mg, O

1.8073

MgO(001) step

heterolytic

Mg, O

-0.4373

γ-Al2O3(110)

heterolytic

Al,O

Cu2O(111)

molecular

CuCUS

-0.1977

-0.9774 -0.7336

-0.47a

this work

ACS Paragon Plus Environment

23

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 28

REFERENCES (1)

Christmann, K. Interaction of Hydrogen with Solid Surfaces. Surf. Sci. Rep. 1988, 9, 1-

163. (2) Copéret, C.; Estes, D. P.; Larmier, K.; Searles, K. Isolated Surface Hydrides: Formation, Structure, and Reactivity. Chem. Rev. 2016, 116, 8463-8505. (3) Kim, J. Y.; Hanson, J. C.; Frenkel, A. I.; Lee, P. L.; Rodriguez, J. A. Reaction of CuO with Hydrogen Studied by Using Synchrotron-based X-ray Diffraction. J. Phys.: Condens. Matter 2004, 16, S3479. (4) Hara, M.; Kondo, T.; Komoda, M.; Ikeda, S.; N. Kondo, J.; Domen, K.; Shinohara, K.; Tanaka, A. Cu2O as a Photocatalyst for Overall Water Splitting under Visible Light Irradiation. Chem. Commun. 1998, 357-358. (5) Kim, J. Y.; Rodriguez, J. A.; Hanson, J. C.; Frenkel, A. I.; Lee, P. L. Reduction of CuO and Cu2O with H2: H Embedding and Kinetic Effects in the Formation of Suboxides. J. Am. Chem. Soc. 2003, 125, 10684-10692. (6) Rodriguez, J. A.; Kim, J. Y.; Hanson, J. C.; Pérez, M.; Frenkel, A. I. Reduction of CuO in H2;In Situ Time-Resolved XRD Studies. Catal. Lett. 2003, 85, 247-254. (7) Norikazu, T.; Hideki, M. Control of Carrier Concentration in Thin Cuprous Oxide Cu2O Films by Atomic Hydrogen. Jpn. J. Appl. Phys. 2002, 41, 5060. (8) Poulston, S.; Parlett, P. M.; Stone, P.; Bowker, M. Surface Oxidation and Reduction of CuO and Cu2O Studied Using XPS and XAES. Surf. Interface. Anal. 1996, 24, 811-820. (9) Cox, D. F.; Schulz, K. H. H2O Adsorption on Cu2O(100). Surf. Sci. 1991, 256, 67-76. (10) Schulz, K. H.; Cox, D. F. Surface Hydride Formation on a Metal Oxide Surface: the Interaction of Atomic Hydrogen with Cu2O(100). Surf. Sci. 1992, 278, 9-18. (11) Nygren, M. A.; Pettersson, L. G. M. H2O Interaction with the Polar Cu2O(100) Surface: A Theoretical Study. J. Phys. Chem. 1996, 100, 1874-1878. (12) Li, Y.; Yan, L.; Wang, G. Adsorption and Dissociation of H2O on Cu2O(100): A Computational Study. J. Nat. Gas. Chem. 2011, 20, 155-161. (13) Maimaiti, Y.; Nolan, M.; Elliott, S. D. Reduction Mechanisms of the CuO(111) Surface through Surface Oxygen Vacancy Formation and Hydrogen Adsorption. Phys. Chem. Chem. Phys. 2014, 16, 3036-3046. (14) Yu, X.; Zhang, X.; Wang, H.; Feng, G. High Coverage Water Adsorption on the CuO(111) Surface. Appl. Surf. Sci. 2017, 425, 803-810. (15) Yu, X.; Zhang, X. High Coverage Water Adsorption on CuO(011) Surface. Phys. Chem. Chem. Phys. 2017, 19, 18652-18659. (16) Stenlid, J. H.; Soldemo, M.; Johansson, A. J.; Leygraf, C.; Gothelid, M.; Weissenrieder, J.; Brinck, T. Reactivity at the Cu2O(100):Cu-H2O Interface: a Combined DFT and PES Study. Phys. Chem. Chem. Phys. 2016, 18, 30570-30584. (17) Soon, A.; Todorova, M.; Delley, B.; Stampfl, C. Thermodynamic Stability and Structure of Copper Oxide Surfaces: A First-principles Investigation. Phys. Rev. B 2007, 75, 125420. (18) Xu, Y.; Wang, H.; Yu, Y.; Tian, L.; Zhao, W.; Zhang, B. Cu2O Nanocrystals: SurfactantFree Room-Temperature Morphology-Modulated Synthesis and Shape-Dependent Heterogeneous Organic Catalytic Activities. J. Phys. Chem. C 2011, 115, 15288-15296. (19) Su, Y.; Li, H.; Ma, H.; Robertson, J.; Nathan, A. Controlling Surface Termination and Facet Orientation in Cu2O Nanoparticles for High Photocatalytic Activity: A Combined Experimental and Density Functional Theory Study. ACS Appl. Mater. Interfaces 2017, 9, 8100-8106. (20) 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.

ACS Paragon Plus Environment

24

Page 25 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(21) Nilius, N.; Fedderwitz, H.; Gro; Noguera, C.; Goniakowski, J. Incorrect DFT-GGA Predictions of the Stability of Non-stoichiometric/polar Dielectric Surfaces: the Case of Cu2O(111). Phys. Chem. Chem. Phys. 2016, 18, 6729-6733. (22) Xu, H.; Wang, W.; Zhu, W. Shape Evolution and Size-Controllable Synthesis of Cu2O Octahedra and Their Morphology-Dependent Photocatalytic Properties. J. Phys. Chem. B 2006, 110, 13829-13834. (23) Kuo, C.-H.; Huang, M. H. Facile Synthesis of Cu2O Nanocrystals with Systematic Shape Evolution from Cubic to Octahedral Structures. J. Phys. Chem. C 2008, 112, 18355-18360. (24) Zheng, Z.; Huang, B.; Wang, Z.; Guo, M.; Qin, X.; Zhang, X.; Wang, P.; Dai, Y. Crystal Faces of Cu2O and Their Stabilities in Photocatalytic Reactions. J. Phys. Chem. C 2009, 113, 1444814453. (25) Önsten, A.; Weissenrieder, J.; Stoltz, D.; Yu, S.; Göthelid, M.; Karlsson, U. O. Role of Defects in Surface Chemistry on Cu2O(111). J. Phys. Chem. C 2013, 117, 19357-19364. (26) Casarin, M.; Maccato, C.; Vigato, N.; Vittadini, A. A Theoretical Study of the H2O and H2S Chemisorption on Cu2O(111). Appl. Surf. Sci. 1999, 142, 164-168. (27) Yu, X.; Zhang, X.; Wang, S.; Feng, G. A First Principle Study on the Magnetic Properties of Cu2O Surfaces. Curr. Appl. Phys. 2015, 15, 1303-1311. (28) Yu, X.; Zhang, X.; Wang, S.; Feng, G. A Computational Study on Water Adsorption on Cu2O(111) Surfaces: The Effects of Coverage and Oxygen Defect. Appl. Surf. Sci. 2015, 343, 33-40. (29) Yu, X.; Zhang, X.; Tian, X.; Wang, S.; Feng, G. Density Functional Theory Calculations on Oxygen Adsorption on the Cu2O Surfaces. Appl. Surf. Sci. 2015, 324, 53-60. (30) Wu, H.; Zhang, N.; Wang, H.; Hong, S. Adsorption of CO2 on Cu2O(111) OxygenVacancy Surface: First-principles Study. Chem. Phys. Lett. 2013, 568–569, 84-89. (31) Bendavid, L. I.; Carter, E. A. CO2 Adsorption on Cu2O(111): A DFT+U and DFT-D Study. J. Phys. Chem. C 2013, 117, 26048-26059. (32) Nie, X.; Griffin, G. L.; Janik, M. J.; Asthagiri, A. Surface Phases of Cu2O(111) under CO2 Electrochemical Reduction Conditions. Cata. Commun. 2014, 52, 88-91. (33) Mishra, A. K.; Roldan, A.; Leeuw, N. H. d. A Density Functional Theory Study of the Adsorption Behaviour of CO2 on Cu2O Surfaces. J. Chem. Phys. 2016, 145, 044709. (34) Riplinger, C.; Carter, E. A. Cooperative Effects in Water Binding to Cuprous Oxide Surfaces. J. Phys. Chem. C 2015, 119, 9311-9323. (35) Zhang, R.; Wang, B.; Ling, L.; Huang, W. Adsorption and Dissociation of H2 on the Cu2O(111) Surface: A Density Functional Theory Study. Appl. Surf. Sci. 2010, 257, 1175-1180. (36) Li, M.; Zhang, J.-y.; Zhang, Y.; Zhang, G.-f.; Wang, T.-m. Density Functional Theory Calculations of Surface Properties and H2 Adsorption on the Cu2O(111) Surface. Appl. Surf. Sci. 2011, 257, 10710-10714. (37) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. (38) Blöchl, P. E. Projector Augmented-wave Method. Phys. Rev. B 1994, 50, 17953-17979. (39) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-zone Integrations. Phys. Rev. B 1976, 13, 5188-5192. (40) Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for Ab Initio Total-energy Calculations Using a Plane-wave Basis Set. Phys. Rev. B 1996, 54, 11169-11186. (41) Kresse, G.; Furthmüller, J. Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comp. Mater. Sci. 1996, 6, 15-50. (42) Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456-1465. (43) 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. 2010, 132, 154104. ACS Paragon Plus Environment

25

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 28

(44) Henkelman, G.; Uberuaga, B. P.; Jonsson, H. A Climbing Image Nudged Elastic Band Method for Finding Saddle Points and Minimum Energy Paths. J. Chem. Phys. 2000, 113, 9901-9904. (45) Henkelman, G.; Arnaldsson, A.; Jonsson, H. A Fast and Robust Algorithm for Bader Decomposition of Charge Density. Comp. Mater. Sci. 2006, 36, 354-360. (46) Le, D.; Stolbov, S.; Rahman, T. S. Reactivity of the Cu2O(100) Surface: Insights from First Principles Calculations. Surf. Sci. 2009, 603, 1637-1645. (47) Hua, Q.; Cao, T.; Gu, X.-K.; Lu, J.; Jiang, Z.; Pan, X.; Luo, L.; Li, W.-X.; Huang, W. Crystal-Plane-Controlled Selectivity of Cu2O Catalysts in Propylene Oxidation with Molecular Oxygen. Angew. Chem. 2014, 126, 4956-4961. (48) Hua, Q.; Cao, T.; Bao, H.; Jiang, Z.; Huang, W. Crystal-Plane-Controlled Surface Chemistry and Catalytic Performance of Surfactant-Free Cu2O Nanocrystals. ChemSusChem 2013, 6, 1966-1972. (49) Schulz, K. H.; Cox, D. F. Photoemission and Low-energy-electron-diffraction Study of Clean and Oxygen-dosed Cu2O(111) and (100) Surfaces. Phys. Rev. B 1991, 43, 1610-1621. (50) Zhang, R.; Liu, H.; Zheng, H.; Ling, L.; Li, Z.; Wang, B. Adsorption and Dissociation of O2 on the Cu2O(111) Surface: Thermochemistry, Reaction Barrier. Appl. Surf. Sci. 2011, 257, 47874794. (51) Zhang, R.; Liu, H.; Li, J.; Ling, L.; Wang, B. A Mechanistic Study of H2S Adsorption and Dissociation on Cu2O(111) Surfaces: Thermochemistry, Reaction Barrier. Appl. Surf. Sci. 2012, 258, 9932-9943. (52) Sun, B.-Z.; Chen, W.-K.; Wang, X.; Lu, C.-H. A Density Functional Theory Study on the Adsorption and Dissociation of N2O on Cu2O(111) Surface. Appl. Surf. Sci. 2007, 253, 7501-7505. (53) Sun, B.-Z.; Chen, W.-K.; Xu, Y.-J. Reaction Mechanism of CO Oxidation on Cu2O(111): A Density Functional Study. J. Chem. Phys. 2010, 133, 154502-154507. (54) Soon, A.; Sohnel, T.; Idriss, H. Plane-wave Pseudopotential Density Functional Theory Periodic Slab Calculations of CO Adsorption on Cu2O(111) Surface. Surf. Sci. 2005, 579, 131-140. (55) Wu, G.; Zhang, J.; Wu, Y.; Li, Q.; Chou, K.; Bao, X. Adsorption and Dissociation of Hydrogen on MgO Surface: A first-principles Study. J. Alloy. Compd. 2009, 480, 788-793. (56) Shen, Y.; Tian, F. H.; Chen, S.; Ma, Z.; Zhao, L.; Jia, X. Density Functional Theory Study on the Mechanism of CO Sensing on Cu2O(111) surface: Influence of the Pre-adsorbed Oxygen Atom. Appl. Surf. Sci. 2014, 288, 452-457. (57) Kokalj, A.; Peljhan, S. Density Functional Theory Study of Adsorption of Benzotriazole on Cu2O Surfaces. J. Phys. Chem. C 2015, 119, 11625-11635. (58) Reuter, K.; Scheffler, M. First-principles Atomistic Thermodynamics for Oxidation Catalysis: Surface Phase Diagrams and Catalytically Interesting Regions. Phys. Rev. Lett. 2003, 90, 46103. (59) Reuter, K.; Scheffler, M. Composition, Structure, and Stability of RuO2(110) as a Function of Oxygen Pressure. Phys. Rev. B 2001, 65, 35406. (60) Li, W.-X.; Stampfl, C.; Scheffler, M. Why is a Noble Metal Catalytically Active? The Role of the O-Ag Interaction in the Function of Silver as an Oxidation Catalyst. Phys. Rev. Lett. 2003, 90, 256102. (61) Sun, Q.; Reuter, K.; Scheffler, M. Effect of a Humid Environment on the Surface Structure of RuO2(110). Phys. Rev. B 2003, 67, 205424. (62) 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 2003, 68, 45407. (63) Bhattacharya, S.; Levchenko, S. V.; Ghiringhelli, L. M.; Scheffler, M. Stability and Metastability of Clusters in a Reactive Atmosphere: Theoretical Evidence for Unexpected Stoichiometries of MgMOx. Phys. Rev. Lett. 2013, 111, 135501. (64) Yu, X.; Oganov, A. R.; Popov, I. A.; Boldyrev, A. I. d-AO Spherical Aromaticity in Ce6O8. J. Comput. Chem. 2016, 37, 103-109. ACS Paragon Plus Environment

26

Page 27 of 28

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(65) Yu, X.; Oganov, A. R.; Popov, I. A.; Qian, G.; Boldyrev, A. I. Antiferromagnetic Stabilization in the Ti8O12 Cluster. Angew. Chem. Int. Edit. 2016, 55, 1699-1703. (66) Stull, D. R.; Prophet, H. JANAF thermochemical tables, DTIC Document, 1971. (67) Yu, X.; Huo, C.-F.; Li, Y.-W.; Wang, J.; Jiao, H. Fe3O4 Surface Electronic Structures and Stability from GGA+U. Surf. Sci. 2012, 606, 872-879. (68) Yu, X.; Zhang, X.; Wang, S. High Coverage Hydrogen Adsorption on the Fe3O4(110) Surface. Appl. Surf. Sci. 2015, 353, 973-978. (69) Yang, T.; Wen, X.-D.; Huo, C.-F.; Li, Y.-W.; Wang, J.; Jiao, H. Structure and Energetics of Hydrogen Adsorption on Fe3O4(111). J. Mol. Cata. A: Chem 2009, 302, 129-136. (70) Yu, X.; Zhang, X.; Jin, L.; Feng, G. CO Adsorption, Oxidation and Carbonate Formation Mechanisms on Fe3O4 Surfaces. Phys. Chem. Chem. Phys. 2017, 19, 17287-17299. (71) Mulakaluri, N.; Pentcheva, R. Hydrogen Adsorption and Site-Selective Reduction of the Fe3O4(001) Surface: Insights From First Principles. J. Phys. Chem. C 2012, 116, 16447-16453. (72) Chen, H.-T.; Choi, Y. M.; Liu, M.; Lin, M. C. A Theoretical Study of Surface Reduction Mechanisms of CeO2(111) and (110) by H2. Chemphyschem 2007, 8, 849-855. (73) Chen, H.-Y. T.; Giordano, L.; Pacchioni, G. From Heterolytic to Homolytic H2 Dissociation on Nanostructured MgO(001) Films As a Function of the Metal Support. J. Phys. Chem. C 2013, 117, 10623-10629. (74) Wischert, R.; Laurent, P.; Copéret, C.; Delbecq, F.; Sautet, P. γ-Alumina: The Essential and Unexpected Role of Water for the Structure, Stability, and Reactivity of “Defect” Sites. J. Am. Chem. Soc. 2012, 134, 14430-14449. (75) Syzgantseva, O.; Calatayud, M.; Minot, C. Hydrogen Adsorption on Monoclinic (1̅11) and (1̅01) ZrO2 Surfaces: A Periodic Ab Initio Study. J. Phys. Chem. C 2010, 114, 11918-11923. (76) García-Melchor, M.; López, N. Homolytic Products from Heterolytic Paths in H2 Dissociation on Metal Oxides: The Example of CeO2. J. Phys. Chem. C 2014, 118, 10921-10926. (77) Hofmann, A.; Clark, S. J.; Oppel, M.; Hahndorf, I. Hydrogen Adsorption on the Tetragonal ZrO2(101) Surface: A Theoretical Study of an Important Catalytic Reactant. Phys. Chem. Chem. Phys. 2002, 4, 3500-3508.

ACS Paragon Plus Environment

27

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 28

TOC Graphic

ACS Paragon Plus Environment

28