催化学报  2020, Vol. 41 Issue (9): 1360-1368      DOI: 10.1016/S1872-2067(20)63564-1   PDF    
扩展功能
加入收藏夹
复制引文信息
加入引用管理器
Email Alert
RSS
本文作者相关文章
Zhu-Yuan Zheng
Dong Wang
Yi Zhang
Fan Yang
Xue-Qing Gong
Structures and reactivities of the CeO2/Pt(111) reverse catalyst: A DFT+U study
Zhu-Yuan Zhenga, Dong Wanga, Yi Zhangb, Fan Yangb, Xue-Qing Gonga     
a. Key Laboratory for Advanced Materials, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, Shanghai 200237, China;
b. State Key Laboratory of Catalysis, Collaborative Innovation Center of Chemistry for Energy Materials, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, Liaoning, China
* Corresponding author. Dong Wang, Tel: +86-21-64252623; Fax: +86-21-64253755; E-mail: wangd@ecust.edu.cn;
Xue-Qing Gong, E-mail: xgong@ecust.edu.cn
This work was supported by National Key R & D Program of China (2018YFA0208602), National Natural Science Foundation of China (21825301, 21573067, 21421004, 21903025) and Program of Shanghai Academic Research Leader (17XD1401400). The authors also thank the National Super Computing Center in Jinan for the computing time.
Abstract: For heterogeneous catalysts, the build-up of interface contacts can influence markedly their activities. Being different from the conventional supported metal/oxide catalysts, the reverse type of oxide/metal structures, e.g. the ceria/Pt composite, have emerged as novel catalytic materials in many fields. However, it remains challenging to determine the optimal interface structure and/or the metal-oxide synergistic effect that can boost catalytic activities. In this work, we conducted density functional theory calculations with on-site Coulomb interaction correction to determine the optimal structures and investigate the physical as well as catalytic properties of various CeO2/Pt(111) composites containing CeO2(111) monolayer, bilayer, and trilayer at Pt(111). We found that the interaction strength between CeO2(111) and Pt(111) substrate first reduces as the ceria slab grows from monolayer to bilayer, and then largely gets converged when the trilayer occurs. Such trend was well rationalized by analyzing the number and distances of O-Pt bonds at the interface. Calculated Bader charges uncovered the significant charge redistribution occurring around the interface, whereas the net electron transfer across the interface is non-significant and decreases as ceria thickness increases. Moreover, comparative calculations on oxygen vacancy formation energies clarified that oxygen removal can be promoted on the CeO2/Pt(111) composites, especially at the interface. We finally employed CO oxidation as a model reaction to probe the surface reactivity, and determined an intrinsic activity order of monolayer CeO2(111) > monolayer CeO2(111)/Pt(111) > regular CeO2(111). More importantly, we emphasized the significant role of the moderate ceria-Pt interaction at the interface that endows the CeO2/Pt reverse catalyst both good thermostability and high catalytic activity. The monolayer CeO2(111)/Pt(111) composite was theoretically predicted highly efficient for catalyzing CO oxidation.
© 2020, Dalian Institute of Chemical Physics, Chinese Academy of Sciences.
Published by Elsevier B.V. All rights reserved.
Key words: CeO2/Pt reverse catalyst    Interface structure    Oxygen vacancy    CO oxidation    Density functional theory    
CeO2/Pt(111)反向催化剂结构和活性的密度泛函理论研究
郑珠媛a, 王栋a, 张毅b, 杨帆b, 龚学庆a     
a. 华东理工大学化学与分子工程学院, 计算化学中心, 工业催化研究所, 结构可控先进功能材料及其制备教育部重点实验室, 上海 200237;
b. 中国科学院大连化学物理研究所, 能源材料化学协同创新中心, 催化基础国家重点实验室, 辽宁大连 116023
摘要:铈基材料因其独特的Ce4+/Ce3+转化性质而广泛运用于非均相催化反应中.尽管在实验和理论上对纯净二氧化铈表面的物理和/或化学性质进行了深入研究,但是与二氧化铈有关的界面结构和反应性能引起了人们的极大兴趣.其中,已有报道表明,氧化铈/金属反向催化剂相较于氧化铈、金属或者金属/氧化铈负载材料能明显提高CO催化氧化和水汽转化等反应活性.然而多数前期研究并没有从理论上给出合理解释,同时也并未说明反向催化剂中氧化铈结构(层数)和性质的关系.可以预见,因受到金属基板的影响,二氧化铈表面的物化性质,如氧空位形成能、电子分布、催化活性等必然会发生变化.本文通过库伦作用校正的密度泛函理论(DFT+U)计算,系统地研究了不同厚度的CeO2/Pt(111)反向催化剂几何结构和电子性质,催化CO氧化的性能.本文首先在Pt(111)载体上明确了单层CeO2(111)的最佳结构,然后研究随着二氧化铈厚度增加,各复合结构界面热力学稳定性、几何结构和电荷性质的变化.计算结果表明:首先,单层CeO2/Pt(111)比双层和三层CeO2/Pt(111)复合结构在界面处表现出更强的相互作用,并且其强度与界面结合结构密切相关,如界面O-Pt键的数量及其长度等;其次,氧化铈板层和Pt基板之间的接触会显著影响界面处一个氧化铈层和两个金属层内的电子分布,使氧化铈外暴露表面的氧空位形成能降低~0.3eV,而界面氧空位形成能则显著降低1.3-1.8eV,并且当表面上沉积≥ 2个氧化铈层时,氧化铈/铂复合材料的物理性能会趋向收敛;最后,通过计算单层CeO2/Pt(111),单层CeO2和模拟体相结构的三层CeO2(111)表面上的CO氧化过程,结果表明三者均遵循Mvk机理,并且关键步骤OCOs偶联的反应能垒分别是0.45,0.33和0.61eV,表明三者的活性趋势为ML CeO2 > ML CeO2/Pt(111)> TL CeO2(111).综合考虑到单层CeO2/Pt(111)界面处适度的二氧化铈-铂相互作用,一方面可以极大提高复合材料热力学稳定性,另一方面还成功保留了单层二氧化铈的优异催化活性,因此单层CeO2/Pt(111)复合材料从理论上认为是一种优异的CO氧化催化剂.
关键词CeO2/Pt(111)反向催化剂    界面结构    氧空位    CO氧化    密度泛函理论    

1 Introduction

Ceria-based materials have found wide applications in many catalytic processes [1], such as vehicle emission purification [2], water-gas shift reactions [3], CO oxidation [4] etc. Their excellent performance primarily stems from the unique Ce 4f orbital with strong electron localization characteristics, typically exhibiting the dynamic Ce4+/Ce3+ transition [5-8]. Accordingly, the process of oxygen vacancy (Ov) formation has aroused great interest, anticipating that the geometric/electronic structure and/or thermostability of Ov may directly affect the performance of CeO2 catalysts [9, 10]. For example, Li et al. [11] found that the subsurface Ov is 0.18 eV more stable than the surface one at CeO2(111), and they also proposed a preferable two-step Ov diffusion mechanism in contrast to the conventional surface hopping mechanism. Nolan et al. [12] and Fabris et al. [13] investigated Ov formation on the three low index CeO2 surfaces and found that the formation energies follow the order of (111) > (110) > (100). Nevertheless, most of previous theoretical studies employed pure ceria models for simulating physical and/or chemical properties, while the significant role of the interface structure in composite catalysts remains largely elusive.

Being different from the conventional metal/oxide catalysts, the reverse type of 'oxide supported on metal (i.e. oxide/metal)' composites have recently emerged as novel catalytic materials with distinct activities [14-16]. In particular, thin ceria films deposited on metal substrates, such as Pt(111) [17], Rh(111) [18, 19], Cu(111) [20], and Au(111) [21], were successfully synthesized and characterized in experiments. Using the low-energy electron diffraction (LEED) technique, Hardacre et al. [22] observed a (4 × 4) geometry pattern at ceria/Pt composites, which corresponds to a 3:4 lattice matching between CeO2(111) and Pt(111). Chan et al. [19] carefully characterized the structure of ceria films on the Rh(111) surface using the LEED, scanning tunneling microscopy (STM), and X-ray photoemission spectroscopy (XPS). They found that the ceria/Rh composites are composed of several atomic layers of CeO2(111) on the Rh(111) substrate, and the deposited CeO2 shows a O:Ce ratio of 1.96:1 with slightly contracted lattice parameter relative to the bulk ceria. There are also STM measurements showing that ultrathin ceria films at metal substrate often yield small islands with a thickness about one or two O-Ce-O layers [23, 24]. Theoretically, Spiel et al. [17] constructed three types of ceria/Pt composites composed of monolayer CeO2(111) on the Pt(111) substrate, and found that aligning surface oxygens of CeO2(111) right on the top sites of interface Pt atoms could enhance the stability of CeO2(111)/Pt(111) composites. Chan and coworkers [15] investigated the thermostabilities of defective CeO2/Rh(111) composites, and found that Ov prefers to form at the interface rather than on the ceria surface. Nevertheless, fine structural resolution of the reverse oxide/metal catalysts at the atomic level is still very rare, hindering the progress in the clarification and regulation of the structure-property relationship.

From the catalytic perspective, CO oxidation was widely considered as a model reaction to probe the reactivity of CeO2-based catalysts [25-27]. As both the geometric and electronic structures of ceria thin films are certainly affected by the metal support, it was anticipated that the reverse ceria/metal catalysts might show distinct catalytic properties compared with the pure ceria materials [26-29]. Recently, by growing well-defined monolayer ceria islands on the Pt(111) surface and monitoring Ce oxidation state, Suchorski et al. [27] observed enhanced CO oxidation activity on the ceria/Pt(111) composites with respect to the free-standing Pt(111) or CeO2(111), and they attributed such enhancement to the optimized electronic structures and accelerated oxygen supply brought by ceria deposition. Besides the ceria/Pt composites, Palomino et al. [30] and Mullins et al. [31] prepared a series of oxide/Au(111) reverse catalysts (e.g. MgO, CeO2, TiO2, Fe2O3), among which they consistently found enhanced oxide-metal interactions and CO oxidation activities relative to the conventional metal/oxide structures. Obviously, the unique catalytic performance of reverse catalysts has been evidenced in many experiments [32, 33], though the fundamental understandings on the origins of the promotion effect and catalytic activities are currently very limited, largely owing to the difficulty in determining reliable interface structures.

In this work, we conducted density functional theory (DFT) calculations corrected by on-site Coulomb interaction (for Ce) to systematically examine the structures and stabilities of various CeO2(111)/Pt(111) reverse catalysts. We first focused on the identification of the optimal structure of monolayer CeO2(111) on the Pt(111) support, and then moved to investigate the changes of the thermodynamic, geometric, and electronic properties at the interface as the ceria thickness increases. It was found that the significant impact of metal substrate on ceria slabs, particularly for interface oxygens, can be directly captured by a simple parameter of the Ov formation energy. After clarifying the physical properties of the CeO2/Pt(111) composites, we used the model reaction of CO oxidation to probe the unique catalytic features of the reverse catalysts, through comparative studies of the monolayer CeO2(111) with and without the Pt(111) support as well as the regular CeO2(111) surface. Our results underline the significant role of the moderate ceria-Pt interaction at the interface that endows the CeO2/Pt reverse catalysts both good thermostability and high catalytic activity.

2 Computational methods

Spin-polarized density functional theory calculations were performed by using the Vienna Ab-initio Simulation Package (VASP) [34, 35]. The exchange-correlation interaction was treated by the Perdew-Wang 1991 (PW91) functional on the basis of the generalized gradient approximation (GGA) [36]. The Kohn-Sham wavefunctions were expanded in a plane wave basis with an energy cutoff of 400 eV, utilizing the projector augmented wave (PAW) method [37, 38]. The Pt 5p, 5d, and 6s, the Ce 4f, 5s, 5p, 5d, 6s, and the 2s, 2p orbitals of carbon and oxygen atoms were treated as valence states. Throughout the DFT calculations [39], we applied an on-site Coulomb interaction correction, i.e. the DFT + U method, with a U value of 5 eV for the Ce 4f orbital, consistent with the previous reports [40-42]. The transition states (TSs) of reactions were located using a constrained optimization scheme [43], and all the calculations for structural optimization and TS search were converged until the forces on each ion were less than 0.05 eV/Å. For the STM simulation, we adopted the Tersoff-Hamann approximation, using a constant-current mode with bias voltage between -2 V and +2 V, to simulate the tunneling current proportional to the local density of states [44]. In addition, we checked the influence of the Van der Waals [45] interaction on both surface reactions and the Ov formation (Table S1), and the negligible energy difference indicates that DFT + U is competent to capture the catalytic trends in CO oxidation on ceria/Pt composites.

To construct the composite models, we adopted the 3:4 lattice-matching pattern observed in experiments, and built a p(3 × 3) CeO2(111) surface cell on the top of a four-layer p(4 × 4) Pt(111) substrate. Considering the lattice parameters of 3.87 Å for ceria and 2.84 Å for Pt, such pattern can cause only small lattice contraction of 2.15% for CeO2(111). The Brillouin-zone integration was treated using a 1 × 1 × 1 k-points mesh. During all calculations, the bottom two layers of Pt substrate were kept fixed, while the other Pt and ceria layers were allowed to relax. We altogether considered three types of CeO2/Pt(111) composites with different thickness of ceria. As one can see from Fig. 1(a‒c), the monolayer (ML), bilayer (BL), and trilayer (TL) CeO2(111) slabs were deposited on the Pt(111) surface with a large vacuum gap (> 10 Å) separating periodic cells, denoted as the ML, BL, and TL CeO2/Pt(111) models, respectively. For calculations of the regular CeO2(111) surface, we used a TL CeO2(111) surface model with the bottom two layers being fixed during calculations, as frequently suggested by previous reports [46].

Fig. 1. Calculated structures of (a, d) ML CeO2/Pt(111), (b, e) BL CeO2/Pt(111) and (c, f) TL CeO2/Pt(111) composites from the side and top views, respectively. The cerium, oxygen and platinum atoms are presented as ivory, red and blue balls, respectively. This color notation is used throughout the manuscript.

In this work, the oxygen vacancy formation energy (EOv) was defined as

where E(O2) is the energy of a gaseous O2 molecule, and Eslab and Evac are energies of the catalyst without and with one oxygen vacancy, respectively.

The adsorption energy of species X (Ead(X)) was calculated according to

where Etotal represents the energy of the interacting system, and EX is the energy of the adsorbate X species in the gas phase. We can also calculate the adsorption energy of ceria at Pt substrate, Ead(ceria), as follows

where EPt and Eceria are energies of the separated Pt and ceria slabs, respectively.

Moreover, we also calculated the normalized adsorption strength of ceria on the substrate with interface adhesion energies (Eadh) according to Eadh = Ead(ceria)/A, where A is the interface area, and defined the deposition energy of ceria, Edepos(ceria), as the adsorption energy of the outmost ceria layer on the underneath CeO2 slabs in the BL and TL CeO2/Pt(111) composites. Obviously, positive Ead(X), Edepos(ceria) or Eadh corresponds to the exothermic adsorption/adhesion process, and the more positive the Ead(X), Edepos(ceria) or Eadh is, the stronger the adsorbate X or ceria slab binds to the surface.

3 Results and discussion
3.1 Structures of CeO2/Pt(111) composites

As the structure of the oxide/metal interface largely determines the thermostability of the corresponding composite catalysts, we first systematically calculated the various adhesion geometries of ML CeO2(111) on the Pt(111) support. In total, 9 types of ML CeO2/Pt(111) structures with different in-plane positions of ceria relative to the Pt substrate are presented in Fig. 2, from which we showed how these structures were constructed by moving the CeO2(111) layer along the x and y axes to the different positions. From the calculated energetics, we found that the interaction between ceria and Pt is generally not strong, and the relative shift of the two materials only causes small impact on the Ead(ceria), showing a narrow energy window of 3.20‒3.50 eV (Table 1). Even for the calculated ML CeO2/Pt(111) structure with the lowest energy (Fig. 2(a)), the calculated Eadh of 0.43 J﹞m‒2 is just comparable to that of the van der Waals junctions in MoS2 nanosheets (0.36‒0.47 J﹞m‒2) [47]. A simple reason for the weak ceria-Pt interaction is that Pt is generally not an oxophilic metal (for example, the Ead(O2) was calculated to be as low as 0.76 eV on the close-packed Pt(111) surface).

Fig. 2. Top views of the calculated structures of ML CeO2/Pt(111) composites (a‒i) and their relative positions (in Å) and stabilities with respect to the optimal one with the lowest energy (x0, y0; E0) in (a). Two lattice vectors are labelled as x and y, respectively.
Table 1
Calculated adsorption and adhesion energies of ceria slabs as well as the structural parameters of the interface structures (the distance and number of Oi-Pt bonds: lOi-Pt and NOi-Pt) of the ML, BL, and TL CeO2/Pt(111) composites. For each type of composite, only the results of the most and least stable structures are listed for comparison.

In Fig. 3(a), we present the interface structure of the most stable ML CeO2/Pt(111) composite from the side view. One can see that the interface oxygens (Oi) located on (or near) the top of substrate Pt atoms can slightly pull the corresponding Pt out of the plane, forming in total 7 Oi-Pt bonds. However, for the less stable structures (Fig. 2(b)‒(i)), the mismatched bonding pattern with slightly longer distances of Oi-Pt bonds (Table 1) can be easily recognized. For example, we determined that the average Oi-Pt distance decreases by 0.02 Å for the least stable composite (Fig. 2(i)) compared to the most stable one (Fig. 2(a)), though they involve the same number of Oi-Pt bond.

Fig. 3. (a) Enlarged interface structure of the ML CeO2/Pt(111) and (b) simulated STM image (bias voltage: +1 V) of the BL CeO2/Pt(111) composite.

With the same approach, we also systematically calculated the possible structures of BL and TL CeO2/Pt(111) composites, and the determined most stable ones are shown in Fig. 1. It needs to mention that Chan et al. [19] have conducted STM measurement of the BL CeO2/Rh(111) composite, which helped them resolve the atomic configurations of the exposed CeO2(111) surface. In the current work, we also performed STM simulation on the most stable BL CeO2/Pt(111) structure (Fig. 3(b)). The ordered bright spots representing surface (or near-surface) Ce cations can be easily recognized, and each surface cell is marked by the dashed line in red. It should be note that, Chan et al. [19] characterized a 5:7 lattice-matching pattern between CeO2(111) and Rh(111) from their experiments. In the current work, we adopted the 3:4 pattern observed in experiments for the CeO2/Pt(111) composites (Fig. 3(b)), and four bright spots along each side of the surface cell can be recognized. Besides, relative lower resolution contrast and absence of one bright spot within each surface cell in the simulated STM image might be largely caused by the rugged surface configuration of the supported CeO2(111).

To determine the effects of ceria thickness on the interface properties, we compared the interface structures of the ML, BL, and TL CeO2/Pt(111) composites with the most and least stable structures (see Table 1). It is evidenced that, for each type of CeO2/Pt(111) composites, the Eadh decreases as the lOi-Pt increases and/or the NOi-Pt decreases (see more samples in Table S2). From their most stable interface structures, we determined the same number of 7 Oi-Pt bonds with the average distances of 2.13, 2.17, 2.16 Å for the ML, BL and TL CeO2/Pt(111) composites, respectively. On the other hand, their least stable interface structures exhibit either identical NOi-Pt but elongated lOi-Pt (e.g. ML and TL composites), or identical lOi-Pt but fewer NOi-Pt (e.g. the BL composite). Most importantly, we determined that the ceria-Pt interaction strength (Ead(ceria) or Eadh) largely decreases as the thickness of the ceria slab increases. Accordingly, the BL and TL CeO2/Pt(111) structures both give longer lOi-Pt distances than the ML CeO2/Pt(111). Such changes can be attributed to the gradually improved bonding strength within ceria slabs as evidenced by the calculated Edepos(ceria) in Table 1. In general, our results regarding the interaction energies and interface structures show that the physical properties of the reverse CeO2/Pt(111) composites basically get converged when ≥ 2 ceria layers are deposited on the Pt substrate.

3.2 Bader charge analysis

Having obtained the optimal composite structures, we then performed electronic structure and Bader charge analyses to investigate the electronic properties of the ML, BL and TL CeO2/Pt(111) composites. From the density plots of calculated charge differences (Fig. 4), one can see that drastic charge redistributions occur at the interface regions for all the three types of composites. Specifically, one oxide layer and two metal layers around the interface were found to be significantly affected, in similar manners to the previous results regarding the metal/TiO2 systems [48]. For the thicker ceria slabs, the top CeO2 layers gradually become less disturbed as they are located away from the interface (Fig. 4(b), and (c)). Accordingly, calculated Bader charges in Table 2 show that each Oi atom is positively charged by ~0.07 |e| and interface Pt is positively charged by ~0.04 |e| for the CeO2/Pt (111) composites with respect to the free-standing ceria slabs and Pt(111) surface; while the charges of the surface oxygens (Os) on the top layer of CeO2/Pt (111) remain basically unchanged (≤ 0.01 |e|). Overall, we determined that the ceria slab only transfers 0.138, 0.031 |e| to the Pt substrate for the ML, BL CeO2/Pt (111) composites, respectively, whereas Pt donates as low as 0.003 |e| to ceria for the TL composite (last row in Table 2). These results show that the interfacial interactions are mainly contributed by the electrostatic interactions between the top Pt layer and the bottom ceria layer at the interface rather than the covalent hybridization. Such feature was also evidenced by the significant electronic polarization at the interface (Fig. 4) and redistribution within the individual ceria or Pt slabs (Table S3). In addition, considering the accumulation of positive charges at Oi and the negligible electronic changes at Os, distinct differences in these two types of O can be further expected in catalytic reactions.

Fig. 4. Calculated charge differential densities (iso-value of 0.003 |e|/Bohr3) of ML (a), BL (b) and TL (c) CeO2/Pt(111) composites caused by the ceria deposition, where the yellow and blue densities represent the electron depletion and accumulation, respectively.
Table 2
Calculated Bader charges (amount of charge per atom) of ML, BL and TL CeO2 slabs before (separated) and after forming contact with the Pt(111) substrate. For the charge transfer results in the last row, positive values indicate that the whole CeO2 slabs donate certain number of electrons to the Pt substrates and vice versa.
3.3 Ov formation

Considering that the excellent oxygen storage/release capacity of ceria is closely related to the Ov formation, [5, 6] we then systematically calculated and compared the EOv at different sites (Oi and Os) of CeO2/Pt(111) composites and on the CeO2(111) surface. Table 3 shows that, for the free-standing CeO2(111) slabs, the Ov formation energy at surface Os site (EOv-s) was computed to be 2.99 and 3.02 eV on the BL and TL surfaces respectively, being ~0.8 eV larger than that on the ML slab (EOv-s = 2.21 eV). These results again indicate that the ML CeO2(111) is quite active and the BL and TL CeO2(111) are less active and may give reduced adhesion strength with the Pt substrate (Table 1).

Table 3
Calculated EOv (in eV) at various oxygen sites of free-standing CeO2(111) and CeO2/Pt(111) composites containing ML, BL, and TL ceria slabs.

In addition, the calculated EOv-s for all the three types of CeO2/Pt(111) composites were found to be ~0.3 eV lower than those in the corresponding free-standing ceria slabs, and such small energy differences are also consistent with the similar calculated Bader charges on Os (Table 2). Interestingly, we found that Ov formation turns to be markedly promoted by 1.0‒1.5 eV at the Oi site (EOv-i) at the interface with the Pt substrate, being 0.92, 1.28 and 1.26 eV for the ML, BL, and TL composites, respectively. The easy removal of the Oi relative to the Os is probably due to the proper structural accommodation to the as-formed Ov-i as provided by the Pt substrate. Besides, these Eov results again demonstrated the converged physicochemical properties of the reverse CeO2/Pt(111) composites when more than two layers of ceria slabs are deposited.

It needs to be particularly noted for the ML CeO2/Pt(111) composite that both the EOv-i and EOv-s are much smaller than that of bulk oxygens in ceria (EOv = 3.05 eV) [49] or other oxides such as TiO2 (EOv = 3.16 eV) [50]. After carefully examining the composite structures, we found that ceria slabs in the ML CeO2/Pt(111) is more flexible than that in the BL and TL composite. As illustrated in Fig. 5, the vertical Ce-Oi and Ce-Os distances (along the z direction; Table S4) in the ML CeO2/Pt(111) can change dynamically to better stabilize the composite, forming different extents of compressed Ce-Os interactions depending on the absence or presence of Ov at the interface (see vertical distances in Fig. 5(b)‒(c) and Ce-O bond length in Table S5). The dynamic distance changes were also observed in the BL and TL CeO2/Pt(111) composites, but their corresponding extents are relatively suppressed (Tables S4 and S5), owing to the increased structural rigidness in ceria as the slab thickness increases.

Fig. 5. Illustrations of variations of vertical distances between Ce and O atomic layers (along the z direction) in (a) free-standing ML CeO2 and (b, c) ML CeO2/Pt(111) composites without and with a Ov-i at the interface. Unit of bond distance in Å.

Particularly for the Ce-Os bond compress induced by the Ov-i at the ML CeO2/Pt (111) composite, we calculated the charge differential density of ML CeO2/Pt(111) composite without (Fig. S1(a)) and with (Fig. S1(b)) a Ov-i. It was found that the presence of Ov-i causes slightly stronger electronic interaction as evidenced by the broader differential density, which accordingly induces structural changes such as bond compression or extension. After comparing the calculated bader charges of the two structures (Table S3), we revealed the following features: (i) the interface first oxygen (O1) layer generally holds similar charge quantity before and after the Ov-i formation, whereas the first cerium (Ce1) layer becomes less positive in the presence of Ov-i, which therefore leads to the elongated O1-Ce1 distance (0.926 Å vs. 0.856 Å; Fig. 5) owing to the reduced electrostatic interaction; (ii) the upper oxygen layer (O2) carries more negative charges (1.222 vs. 1.189 |e|) after Ov-i formation, and taking into account the already reduced O1-Ce1 interaction, it is reasonable that the Ce1 would bind the O2 layer stronger with compressed Ce‒O distance; (iii) the overall charge transfer between ceria and Pt substrate remains fairly small (0.273 |e|) after Ov-i formation, which again signifies the intensive electrostatic interactions rather than the covalent hybridization between the ceria slab (containing Ov at the interface) and Pt substrate.

Overall, these geometric and electronic features regarding the compressed Ce‒Os (or extended Ce‒Oi) distances and stronger structural relaxation in the ML CeO2/Pt(111) could reasonably account for the lower Ov formation energies. Given that Ov formation is widely regarded as an important step in ceria catalysis (e.g. CO oxidation [51]), it is therefore expected that the ML CeO2/Pt(111) composite might show excellent catalytic activities.

3.4 CO oxidation

Regarding the catalytic performance, we took CO oxidation as a model reaction to study the reactivities of the reverse CeO2/Pt(111) catalysts. We first calculated CO adsorption at the ML, BL, TL composites and compared to that on the regular CeO2(111) surface (using a TL model). It was found that ML CeO2/Pt(111) gives the largest Ead(CO) of 0.20 eV, being twice the value on the CeO2(111), and the adsorption strength decreases as the ceria slab becomes thicker. These results agree well with experimental observations that the CO adsorption capacity (including both the amount and the strength of CO adsorption) decreases as more ceria are deposited on the metal substrate [26], and also indicate that the ML CeO2/Pt(111) might be more active than other composites. In addition, considering the weak adsorption of CO, one may also anticipate that CO oxidation via the Langmuir-Hinshelwood (LH) mechanism [52] is not likely to occur on the CeO2/Pt(111) composites. Therefore, we mainly focused on the Mars-van-Krevelen (MvK) mechanism [53] as widely validated for ceria based catalysts [54, 55], and compared the calculated energetics of CO oxidation at ML CeO2(111) and ML CeO2/Pt(111) with those on the regular CeO2(111) surface.

The calculated reaction pathways of CO oxidation (toward CO2 production) via the MvK mechanism on the unsupported ML and regular CeO2(111) surfaces are shown in Fig. 6. This process starts with a gas-phase CO molecule (IS) which then forms weak interaction with the surface (IM1) and combines with an Os to form the OC-Os bond. At the coupling TS, the OC-Os distance was calculated to be ~1.50 Å while the O-C-Os angle is ~120º. We determined a coupling barrier of 0.61 eV on the regular CeO2(111) surface, but a low barrier of 0.33 eV on the ML CeO2(111). After the coupling, the linear CO2 (IM2) occurs with an exothermic energy of 0.73 and 1.38 eV on the regular and ML CeO2(111), respectively. Obviously, the lower barrier and larger energy release on the ML CeO2(111) relative to the regular surface can be rationalized by their different EOv-s values (2.21 vs. 3.02 eV; Table 3). For the last step of CO2 desorption, the as-formed CO2 can leave both surfaces quite easily by overcoming its adsorption energies of ~0.38 eV, leaving an Ov on the surface. Note that previous reports have suggested that the surface Ov at ceria surfaces can be readily recovered [56, 57] and consequently the overall reactivity of CO oxidation can be largely determined by the OC···O coupling. Therefore, judging from the above calculated energetics, we expect that the ultrathin ML CeO2(111) surface can show much higher CO oxidation activity than bulk ceria materials.

Fig. 6. Calculated energy profile and key intermediate structures of CO oxidation via the MvK mechanism on the ML (black) and regular (red) CeO2(111) surfaces. C atoms are in grey and the top-surface O are in green.

For CO oxidation on the ML CeO2/Pt(111) composite (Fig. 7), the initial reaction steps are very similar to those on CeO2(111) surfaces, leading to the formation of surface CO2 (IM2), though lower coupling barrier (0.45 vs. 0.61 eV) and larger energy release (1.69 vs. 0.73 eV) were also determined. These results are again consistent with the vast Eov difference between the two catalysts (1.92 vs. 3.02 eV; Table 3). In contrast, CO oxidation on the ML composite could form two types of Ov defects at the Os or Oi site (FS and FS' in Fig. 7) after CO2 desorption. The former Ov-s pathway shows no obvious difference to that on free-standing ceria surfaces, releasing gaseous CO2 directly with a small endothermic energy of 0.36 eV; while the latter Ov-i pathway involves the migration of an interface Oi to fill the surface Ov-s (IM3), which, however, turns to be thermodynamically more preferable by ~1.0 eV than IM2.

Fig. 7. Calculated energy profile and key intermediate structures of CO oxidation via the MvK mechanism on the ML CeO2/Pt(111) composite. Two possible final states exposing Ov at the Os and Oi site are labelled as FS and FS', respectively.

In addition, though the result of much lower EOv-i than EOv-s was also evidenced in the BL and TL composites (Table 3), Oi migration requires multi-steps of oxygen relay within such thicker ceria slabs, thus being a kinetically much slower process. In this regard, we calculated the energy barrier of Ov migration in the ML CeO2/Pt(111) composite using the climbing image nudged elastic band (CI-NEB) method [58]. We determined an energy barrier of 0.92 eV, and the dynamic picture of Ov migration was presented in Fig. S2. Therefore, our results indicate that the additional CO oxidation pathway involving Ov-i formation is an unique feature of the ML CeO2/Pt(111) or perhaps other ML ceria/metal composites. In fact, this speculation was supported by experimental observation of Ov-i formation in the ML CeO2/Rh(111) catalyst under flowing CO atmosphere [18].

Finally, though an higher OC···O coupling barrier of 0.45 eV was determined on the ML CeO2/Pt(111) than that on the ML CeO2(111) film (0.33 eV), the calculated results still indicate that the moderate ceria-Pt interaction at the interface can not only greatly stabilize the composite structure, realizing the synthetization of ultrathin ceria/metal reverse catalysts in experiments [59, 60], and it may also largely preserve the superior reactivity of the ML ceria slab. Therefore, the ML ceria/Pt composite reaches a good balance between the thermostability and the reactivity, manifesting itself superior catalysts for catalyzing CO oxidation.

4 Conclusions

In this work, we performed DFT + U calculations to systematically investigate the geometric, electronic, and catalytic properties of the reverse type of CeO2/Pt(111) composites. The following main conclusions can be obtained.

The ceria-Pt interaction strength is closely relevant with the distance and number of Oi-Pt bonds formed at the interface, and ML CeO2/Pt(111) shows larger adhesion interaction (0.43 J﹞m‒2) than the BL or TL composites. Increased thickness of the ceria slab in the CeO2/Pt(111) composites is highly favored for their thermostability, though their physical properties basically get converged when ≥ 2 ceria layers being deposited on the Pt substrate.

The contact between the ceria slab and the Pt substrate can significantly affect the electron distribution of the one ceria layer and two metal layers around the interface, and lowers the Ov formation energy by ~0.3 eV at the exposed ceria surface whereas as high as 1.3‒1.8 eV for the interface Oi as compared with the free-standing CeO2(111) slabs.

The key step of OC···O coupling in CO oxidation gives the reaction barriers of 0.61, 0.33, and 0.45 eV on the regular and ML CeO2(111) surfaces and ML CeO2/Pt(111) composite, respectively, indicating an intrinsic activity trend of ML CeO2(111) > ML CeO2/Pt(111) > regular CeO2(111). Moderate ceria-Pt interaction at the interface not only greatly stabilizes the composite structure, but also successfully preserves the superior reactivity of the ML ceria slab. Our results demonstrate the ML CeO2/Pt(111) composite as a promising type of catalyst with high activity for CO oxidation.

Conflicts of interest

There are no conflicts of interest to declare.

References
[1]
H. F. Wang, X. Q. Gong, Y. L. Guo, Y. Guo, G. Z. Lu, P. Hu, J. Phys. Chem. C, 2009, 113, 10229-10232. DOI:10.1021/jp900942a
[2]
X. Zhang, N. Rui, X. Jia, X. Hu, C.-J. Liu, Chin. J. Catal., 2019, 40, 495-503. DOI:10.1016/S1872-2067(19)63289-4
[3]
X. P. Fu, L. W. Guo, W. W. Wang, C. Ma, C. J. Jia, K. Wu, R. Si, L. D. Sun, C. H. Yan, J. Am. Chem. Soc., 2019, 141, 4613-4623. DOI:10.1021/jacs.8b09306
[4]
Z. Guo, L. Song, T. Xu, D. Gao, C. Li, X. Hu, G. Chen, Mater. Chem. Phys., 2019, 226, 338-343. DOI:10.1016/j.matchemphys.2019.01.057
[5]
B. Liu, J. Liu, T. Li, Z. Zhao, X. Q. Gong, Y. Chen, A. Duan, G. Y. Jiang, Y. C. Wei, J. Phys. Chem. C, 2015, 119, 12923-12934. DOI:10.1021/acs.jpcc.5b00267
[6]
T. Montini, M. Melchionna, M. Monai, P. Fornasiero, Chem. Rev., 2016, 116, 5987-6041. DOI:10.1021/acs.chemrev.5b00603
[7]
K. Jia, H. Zhang, W. Li, Chin. J. Catal., 2008, 29, 1089-1092. DOI:10.1016/S1872-2067(09)60006-1
[8]
B. L. Liu, Y. Li, Y. L. Cao, L. Wang, S. J. Qing, K. Wang, D. Z. Jia, Eur. J. Inorg. Chem., 2019, 2019, 1714-1723. DOI:10.1002/ejic.201801451
[9]
L. Chen, J. H. Li, M. F. Ge, L. Ma, H. Chang, Chin. J. Catal, 2011, 32, 836-841. DOI:10.1016/S1872-2067(10)60195-7
[10]
M. Wang, M. Shen, X. Jin, J. Tian, M. Li, Y. Zhou, L. X. Zhang, J. L. Shi, ACS Catal., 2019, 9, 4573-4581. DOI:10.1021/acscatal.8b03975
[11]
H. Y. Li, H. F. Wang, Y. L. Guo, G. Z. Lu, P. Hu, Chem. Commun., 2011, 47, 6105-6107. DOI:10.1039/c1cc11226k
[12]
M. Nolan, S. C. Parker, G. W. Watson, Surf. Sci., 2005, 595, 223-232. DOI:10.1016/j.susc.2005.08.015
[13]
S. Fabris, G. Vicario, G. Balducci, S. de Gironcoli, S. Baroni, J. Phys. Chem. B, 2005, 109, 22860-22867. DOI:10.1021/jp0511698
[14]
H. Yan, C. Yang, W. P. Shao, L. H. Cai, W. W. Wang, Z. Jin, C. J. Jia, Nat. Commun., 2019, 10, 3470. DOI:10.1038/s41467-019-11407-2
[15]
G. H. Li, S. W. Hu, Q. Xu, J. F. Zhu, J. Phys. Chem. C, 2019, 123, 23563-26571. DOI:10.1021/acs.jpcc.9b07146
[16]
Y. Zhang, W. Feng, F. Yang, X. Bao, Chin. J. Catal., 2019, 40, 204-213. DOI:10.1016/S1872-2067(18)63171-7
[17]
C. Spiel, P. Blaha, Y. Suchorski, K.Schwarz, G. Rupprechter, Phys. Rev. B, 2011, 84, 045412/1-045412/11.
[18]
L. H. Chan, J. Yuhara, Surf. Sci., 2017, 661, 69-76. DOI:10.1016/j.susc.2017.02.008
[19]
L. H. Chan, J. Yuhara, J. Phys. Chem. C, 2015, 143, 074708/1-074708/8.
[20]
S. D. Senanayake, J. T. Sadowski, J. Evans, S. Kundu, S. Agnoli, F. Yang, D. Stacchiola, J. I. Flege, J. Hrbek, J. A. Rodriguez, J. Phys. Chem. Lett., 2012, 3, 839-843. DOI:10.1021/jz300159p
[21]
S. Ma, J. Rodriguez, J. Hrbek, Surf. Sci., 2008, 602, 3272-3278. DOI:10.1016/j.susc.2008.08.021
[22]
C. Hardacre, G. M. Roe, R. M. Lambert, Surf. Sci., 1995, 326, 1-10. DOI:10.1016/0039-6028(94)00783-7
[23]
D. C. Grinter, R. Ithnin, C. L. Pang, G. Thornton, J. Phys. Chem. C, 2010, 114, 17036-17041. DOI:10.1021/jp102895k
[24]
P. Luches, F. Pagliuca, S. Valeri, F. Boscherini, J. Phys. Chem. C, 2013, 117, 1030-1036.
[25]
C. C. Zong, C. L. Wang, L. Hu, R. R. Zhang, J. Chen, L. Z. Wei, Q. W. Chen, Inorg. Chem., 2019, 58, 14238-14243. DOI:10.1021/acs.inorgchem.9b02356
[26]
J. A. Rodriguez, P. Liu, J. Graciani, S. D. Senanayake, D. C. Grinter, D. Stacchiola, J. Hrbek, J. Fernandez-Sanz, J. Phys. Chem. Lett., 2016, 7, 2627-2639. DOI:10.1021/acs.jpclett.6b00499
[27]
Y. Suchorski, R. Wrobel, S. Becker, H. Weiss, J. Phys. Chem. C, 2008, 112, 20012-20017. DOI:10.1021/jp806033v
[28]
L. Barrio, J. M. Campos-Martín, J. L. G. Fierro, J. Phys. Chem. A, 2007, 111, 2166-2171.
[29]
J. A. Rodriguez, J. Graciani, J. Evans, J. B. Park, F. Yang, D. Stacchiola, S. D. Senanayake, S. Ma, M. Perez, P. Liu, F. Sanz, F. Javier, J. Hrbek, Angew. Chem. Int. Ed., 2009, 48, 8047-8050. DOI:10.1002/anie.200903918
[30]
R. M. Palomino, R. A. Gutiérrez, Z. Y. Liu, S. Tenney, D. C. Grinter, E. Crumlin, I. Waluyo, P. J. Ramirez, J. A. Rodriguez, S. D. Senanayake, ACS Sustain. Chem. Eng., 2017, 5, 10783-10791. DOI:10.1021/acssuschemeng.7b02744
[31]
T. Yan, D. W. Redman, W. Y. Yu, D. W. Flaherty, J. A. Rodriguez, C. B. Mullins, J. Catal., 2012, 294, 216-222. DOI:10.1016/j.jcat.2012.07.024
[32]
T. R. Reina, S. Ivanova, M. A. Centeno, J. A. Odriozola, Appl. Catal. B, 2016, 187, 98-107. DOI:10.1016/j.apcatb.2016.01.031
[33]
B. Jenewein, M. Fuchs, K. Hayek, Surf. Sci., 2003, 532-535.
[34]
G. Kresse, J. Hafner, Phys. Rev. B, 1993, 48, 13115-13118. DOI:10.1103/PhysRevB.48.13115
[35]
B. Song, L. M. Liu, C. Y. Yam, J. Phys. Chem. Lett., 2019, 10, 5564-5570. DOI:10.1021/acs.jpclett.9b02048
[36]
J. P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett., 1996, 77, 3865-3868. DOI:10.1103/PhysRevLett.77.3865
[37]
P. E. Blochl, Phys. Rev. B, 1994, 50, 17953-17979. DOI:10.1103/PhysRevB.50.17953
[38]
H. Levämäki, M. Kuisma, K. Kokko, J. Chem. Phys., 2019, 150, 54101/1-054101/9..
[39]
G. Kresse, J. Furthmuller, Phys. Rev. B, 1996, 54, 11169-1186. DOI:10.1103/PhysRevB.54.11169
[40]
D. O. Scanlon, N. M. Galea, B. J. Morgan, G. W. Watson, J. Phys. Chem. C, 2009, 11311095-11103.
[41]
S. D. Huang, C. Sheng, Z. P. Liu, J. Chem. Phys., 2019, 151, 17472..
[42]
M. Nolan, S. Grigoleit, D. C. Sayle, S. C. Parker, G. W. Watson, Surf. Sci., 2005, 576, 217-229. DOI:10.1016/j.susc.2004.12.016
[43]
H. Y. Li, H. F. Wang, X. Q. Gong, Y. L. Guo, Y. Guo, G. Z. Lu, P. Hu, Phys. Rev. B, 2009, 102, 193401/1-193401/4..
[44]
J. Tersoff, D. R. Hamann, Phys. Rev. Lett., 1983, 50, 1998-2001. DOI:10.1103/PhysRevLett.50.1998
[45]
J. W. Cao, F. Y. Li, W. S. Xia, W. S. Bian, Chin. J. Chem. Phys., 2019, 32, 157-166. DOI:10.1063/1674-0068/cjcp1901007
[46]
F. Chen, D. Liu, J. Zhang, P. Hu, X. Q. Gong, G. Lu, Phys. Chem. Chem. Phys., 2014, 14, 16573-16580.
[47]
Z. Y. Zhao, Q. L. Liu, Catal. Sci. Technol., 2018, 8, 1867-1879. DOI:10.1039/C7CY02252B
[48]
H. R. Chen, P. Li, N. Umezawa, H. Abe, J. Ye, K. Shiraishi, A. Ohta, S. Miyazaki, J. Phys. Chem. C, 2016, 120, 5549-5556.
[49]
R. S. Mulliken, J. Chem. Phys, 1933, 1, 492-503. DOI:10.1063/1.1749322
[50]
L. Wang, F. Meng, Mater. Res. Bull., 2013, 48, 3492-3498. DOI:10.1016/j.materresbull.2013.05.036
[51]
D. Wang, T. Sheng, J. F. Chen, H. F. Wang, P. Hu, Nat. Catal., 2018, 1, 291-299. DOI:10.1038/s41929-018-0055-z
[52]
Z. L. Wu, M. J. Li, S. H. Overbury, J. Catal., 2012, 285, 61-73. DOI:10.1016/j.jcat.2011.09.011
[53]
C. N. Costa, S. Y. Christou, G. Georgiou, A. M. Efstathiou, J. Catal., 2013, 219, 259-272.
[54]
N. J. Lawrence, J. R. Brewer, L. Wang, T. S. Wu, J. Wells-Kingsbury, M. M. Ihrig, G. H. Wang, Y. L. Soo, W. N. Mei, C. L. Cheung, Nano. Lett., 2011, 11, 2666-2671. DOI:10.1021/nl200722z
[55]
Y. L. Song, L. L. Yin, J. Zhang, X. Q. Gong, G. Z. Lu, Surf. Sci., 2013, 618, 140-47. DOI:10.1016/j.susc.2013.09.001
[56]
W. Zhang, Z. Y. Li, Y. Luo, J. L. Yang, J. Phys. Chem. C, 2009, 113, 8302-8305. DOI:10.1021/jp810751j
[57]
B. Liu, J. Liu, T. Li, Z. Zhao, X. Q. Gong, Y. Chen, A. J. Duan, G. Y. Jiang, Y. C. Wei, J. Phys. Chem. C, 2015, 119, 12923-12934. DOI:10.1021/acs.jpcc.5b00267
[58]
H. Aljam, F. Abild-Pedersen. Phys. Chem. Chem. Phys., 2019, 21, 25328-25333. DOI:10.1039/C9CP04897A
[59]
X. W. Liu, K. B. Zhou, L. Wang, B. Y. Wang, Y. D. Li, J. Am. Chem. Soc., 2009, 113, 3140-3141.
[60]
C. Breinlich, J. M. Essen, E. Barletta, K. Wandelt, Thin Solid Films, 2011, 519, 3752-3755. DOI:10.1016/j.tsf.2011.01.353