A DFT Study

temperature NO oxidation, one possible strategy is to increase the surface O2 coverage on Cr2O3(012), or to weaken the binding strength of NO2 on the ...
11 downloads 6 Views 1MB Size
Subscriber access provided by La Trobe University Library

Article

Insight into Room-temperature Catalytic Oxidation of Nitric oxide by Cr2O3: A DFT Study Jiamin Jin, Ningling Sun, Wende Hu, Haiyang Yuan, Haifeng Wang, and Peijun Hu ACS Catal., Just Accepted Manuscript • DOI: 10.1021/acscatal.8b00081 • Publication Date (Web): 03 May 2018 Downloaded from http://pubs.acs.org on May 3, 2018

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 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 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.

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 23 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

ACS Catalysis

Insight into Room-temperature Catalytic Oxidation of Nitric oxide by Cr2O3: A DFT Study Jiamin Jin1, Ningling Sun1, Wende Hu1, Haiyang Yuan1, Haifeng Wang1* and Peijun Hu1,2* 1

Key Laboratories 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, P. R. China. Email: [email protected] 2 School of Chemistry and Chemical Engineering, Queen’s University Belfast, Belfast, BT9 5AZ, UK. Email: [email protected]

ABSTRACT: Cr-based catalysts has drawn attentions as a promising roomtemperature NO oxidation catalyst. However, the intrinsic active component and reaction mechanism at the atomic level remain unclear. Here, taking the Cr2O3, one of the most stable chromium oxides, as object, we systematically investigated NO oxidation processes on Cr2O3(001) and (012) surface by virtue of DFT+U calculations, aiming to uncover the activity limiting factors and basic structure-activity relationship of Cr2O3 catalyst. It was revealed that NO oxidation could not proceed via Mars-vanKrevelen mechanism involving the lattice oxygen on both surfaces. For Cr2O3(001) surface exposing the isolated three-coordinated Cr3c, the reactions are inclined to occur through the Eley-Rideal route, in which the NO couples directly with the molecular O2* or atomic O* adsorbed at the Cr3c site to form two key intermediate species (ONOO* and NO2*) following a barrierless process. Nevertheless, the overall activity is limited by the irreversible adsorption of NO2 species on the highly unsaturated Cr3c. In contrast, on the (012) termination which exposes the five-coordinated Cr5c, the NO2* can be easily released, but the reactant O2 cannot be efficiently adsorbed and also results in a limited overall activity at the room temperature. To achieve a higher activity, a thermodynamically favored interface model of mono-chain CrO3 supported on Cr2O3(012) was proposed, which shows an improved O2 adsorption energy of -0.99 eV and thus an enhanced activity of Cr2O3(012), possibly accounting for the experimentally high activity of Cr-based catalysts usually involving the Cr3+/Cr6+ redox. This study demonstrated the catalytic ability of Cr2O3 for NO oxidation at room temperature, and the presented systematic picture may facilitate the further design of more active Cr-based catalyst. KEYWORDS: density functional theory, chromium oxides; NO oxidation; room temperature; catalytic mechanism; mono-chain CrO3

1

ACS Paragon Plus Environment

ACS Catalysis 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

1.

Introduction NO oxidation to NO2 is a key process in a number of NOx remediation technologies1,

including selective catalytic reduction (SCR)2-4, NOx storage and reduction (NSR)5-6. Platinum-based catalysts7-10 have been widely used for NO oxidation; however, the high costs restrict their applications. Transition metal oxides (TMOs) have recently attracted great attentions as promising catalysts for NO oxidation, and previous studies have accomplished much progress in finding highly efficient TMOs catalysts. The perovskites such as the La1-xSrxCoO3, La1-xSrxMnO3, and Mn-mullite (Sm, Gd)Mn2O5 are expected to be viable alternatives to precious metal catalysts.11-14 These observations/suggestions are particularly noteworthy because these catalysts are much less expensive than Pt group metal based catalysts. However, they are still not satisfactory for NO oxidation at room temperature. It has been reported that the activated carbon, especially modified carbon15 and activated carbon fibers16, are active for catalytic NO oxidation at room temperature, but the catalytic activity diminished quickly because of the irreversible adsorption of NO2. Thus, exploring non-precious and low-temperature catalysts for NO oxidation and understanding the catalytic mechanism are of great importance. Recently, there has been great interest in the development of Cr-based catalysts17-24 for the effective removal of NOx from stationary sources. Smirniotis et al.17, 25 found that Cr oxides supported on Hombikat TiO2 showed superior performance on lowtemperature selective catalytic reduction (SCR) of NO with NH3, and pointed out the existence of Cr2O3 crystal phase. Shiba et al.26 early reported that TiO2 supported chromium component had a high ability for the oxidation of NO to NO2. Cai et al.27 studied the effect of Cr on the NO oxidation over the ceria-zirconia solid solutions, proposing that the synergistic mixture of the well dispersed CrOx and Ce0.5Zr0.5O2 exhibited higher catalytic activity than the separate ones. More recently, Liu et al.22 found that chromium oxides prepared by ammonia precipitation method showed higher NO to NO2 conversions at ambient temperature, and suggested that the amorphous phase structure (neither Cr2O3 nor CrO3) of the prepared chromic oxide catalysts is critical in obtaining high performance and long-time stability for NO conversion. These attractive findings represent much progress because Cr-based oxides are significantly less expensive and more abundant, and more importantly, they were demonstrated to possess excellent performances on complete oxidation of NO at room temperature.22-24 However, several important issues on Cr-based catalysts remain to be addressed. Firstly, as a metal oxide with the multivalence, the real active component for NO oxidation is 2

ACS Paragon Plus Environment

Page 2 of 23

Page 3 of 23 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

ACS Catalysis

still controversial, despite that Cr6+ species with high dispersion was suggested to be more favorable than Cr2O3 to oxidize NO very recently.28-30 Secondly, the favorable reaction path behind the room-temperature NO oxidation is not yet clear. Thirdly, as the most common chromium oxide, what the role Cr2O3 might play in NO conversion remains elusive. Specially, one may ask the following questions: (i) How easily can the involved molecules (O2, NO and NO2 etc.) in NO oxidation chemisorb and convert on Cr2O3? (ii) Is the catalytic activity surface sensitive or insensitive to Cr2O3? (iii) What are the advantages and shortcomings of Cr2O3 for catalyzing room-temperature NO oxidation? Could Cr2O3 cooperate with CrO3 to avoid the shortcomings and obtain an enhanced activity? Inspired by these questions, therefore, it is highly desirable to investigate the nature of Cr2O3 for NO oxidation theoretically. Aiming at revealing the abovementioned issues, in this work a molecular-level investigation of catalytic oxidation of the nitric oxide on Cr2O3 and at the Cr2O3/CrO3 interface is studied: Two typically exposed surfaces of Cr2O3, (001) and (012) surface, are examined to address (i) the catalytic mechanism; (ii) the activity-limiting factor of Cr2O3; (iii) the basic structure-activity relationship of Cr2O3 on NO oxidation; and (iv) the possible Cr2O3/CrO3 interface model with improved activity. It was noted that the first principle studies have been used extensively to understand NO catalytic oxidation on metal,10, 31-32 metal oxides33-35 and other systems36-37. However, to the best of our knowledge, this is the first computational study focusing on the room-temperature NO oxidation catalyzed by chromium oxides. 2. Calculation Method and Model All the spin-polarized calculations in this work were performed within the framework of density functional theory (DFT) using the Vienna ab initio simulation package (VASP).38-39 Generalized gradient approximation (GGA) with Perdew-BurkeErnzerhof (PBE)40 functional was employed to describe the exchange and correlation energy. Electron-ion interaction were treated within the projected augmented wave (PAW) method.41 The valence electronic states were expanded in plane wave basis sets with a cutoff energy of 450 eV. The Cr[3d54s1], Cr[3p63d54s1], O[2s22p4] and N[2s22p3] were treated as valence states for describing Cr3+ in Cr2O3, Cr6+ in CrO3, O and N elements, respectively. We used a Monkhorst-Pack K-mesh of 8×8×3 for the bulk of Cr2O3. For both the Cr2O3(001) and (012) surfaces, a p(2×2) supercell with four Cr2O3unit layers was used, and the vacuum between slabs are 12 Å, in which a 2×2×1 Kpoint mesh was used. During structural optimization, the bottom two layers of the slab were fixed at bulk truncated position and the top two layers and the adsorbates were 3

ACS Paragon Plus Environment

ACS Catalysis 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

fully relaxed. All the adsorption geometries were optimized until the forces on all atoms were below 0.05 eV/Å, and the transition states were located with a constrained minimization technique.42-45 The DFT-D3(BJ) method was used to describe the dispersion effects in the system.46

Fig. 1 hexagonal unit cell of Cr2O3 (a) together with the rhombohedral primitive cell (b). Red spheres stand for O atoms; blue spheres represent Cr atoms. The Cr atoms denoting with ↑↓ stand for the spin arrangements. (c) and (d) are the top views of the optimized surface configuration of the Cr2O3(001) and (012) surface, respectively. This notation is used throughout the paper.

Bulk Cr2O3 is an antiferromagnetic (AFM) semiconductor, which crystallizes in the hexagonal corundum structure (space group R3� c, see Fig. 1a). Three different AFM

configurations: + + - -, + - - +, and + - + - could exist for the Cr cations in the Rhombohedral primitive cell (see Fig. 1b), where “+” and “-” denotes spin up and spin down, respectively. The tests showed that the equilibrium magnetic order of this material in its pure state is the AFM + - + - configuration (shown in Fig. 1), which is consistent with previous studies.47 In addition, to correctly describe the strong

correlation of Cr2O3 materials, the DFT+U method, proved to have a significantly improved description of the structure, magnetic, and electronic properties of material,4850

has been adopted. To obtain a proper U term for Cr-3d orbital, the electronic structure

(the magnetic moment of Cr3+ and the band gap of Cr2O3 bulk) and the chemical properties such as the formation energy of Cr2O3 (referring to 1/2Cr2O3 + 1/4O2 → CrO2) and NO adsorption energy on Cr2O3(001) surface, which reflect the Cr-O bond strength and the adsorption ability of Cr3+ site, respectively, are calculated under different U values. As shown in Table 1, compared with the experimental value for each term, a proper value of Ueff=3.0 eV can be obtained on comprehensive consideration, which basically gives a minor error for all the terms, especially the two 4

ACS Paragon Plus Environment

Page 4 of 23

Page 5 of 23 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

ACS Catalysis

chemical properties. For example, the calculated formation energy of Cr2O3 and the NO adsorption energy on Cr2O3(001) surface is -0.38 eV and -0.79 eV at Ueff=3.0 eV, respectively, being close to the respective experimental value of -0.29 eV and -0.80 eV. Also, the identified Ueff value here is close to the findings of Shi et al.50 and Mosey et al.51. Therefore, the parameter of Ueff=3.0 eV for Cr2O3 are utilized hereinafter. The adsorption energy, Ead, of an adsorbate on the surface is calculated as Ead = Esurf/adsorbate -Esurf -Eadsorbate, where Esurf/adsorbate and Esurf are the total energy of the surface slab with and without the adsorbate, respectively, and Eadsorbate is the energy of the NO or other small molecules in vacuum. Based on the above definition, a negative Ead indicates favorable (exothermic) adsorption. Notably, the Gibbs free energy change (∆G) of each elementary step was estimated according to ∆G =∆H -T∆S; the reaction enthalpy ΔH is approximated with the total energy difference (∆E) including the zeropoint energy correction (∆ZPE) and neglecting the small heat capacity correction and ∆(pV) term. Regarding the entropy effect term (T∆S), it is usually trivial for the surface reactions with no adsorption/desorption involved, owing to that the entropies of the surface species are typically small and can also be largely cancelled between the initial state and the transition state or the final state. For the adsorption/desorption processes, the large entropy contribution of gaseous molecules (T∆S), including the vibrational, rotational and translational entropies, have to be considered to estimate ∆G at a given temperature52-55. Table 1 Formation energies of Cr2O3 with respect to the reaction (1/2Cr2O3 + 1/4O2 → CrO2, ∆H/eV). The adsorption energies (Ead/eV) of NO on the Cr2O3 (001) surface, the magnetic moment of Cr3+ and band gap of Cr2O3 on a series of U values are also listed. For a comparison, the experimental values are also given. U

0

1

2

3

4

5

exp

ΔH /eV

-1.10

-0.87

-0.62

-0.38

-0.14

0.05

-0.29

Ead /eV

-1.79

-1.34

-0.99

-0.79

-0.56

-0.44

-0.80

Mag. Moment /μB

2.62

2.74

2.82

2.88

2.94

2.98

2.76

Band Gap /eV

1.59

2.03

2.46

2.84

3.13

3.32

3.30

3.

Results and discussion

3.1 NO oxidation on Cr2O3(001) surface. The optimized surface structure of Cr2O3(001) is shown in Fig. 1c, from which one can see that the surface is terminated by three-coordinated Cr3c and three-coordinated lattice oxygen (O3c), and the distance of the adjacent Cr3c-Cr3c on the (001) surface is very long (5.051 Å), implying that the Cr3c site is basically isolated. Notably, the 5

ACS Paragon Plus Environment

ACS Catalysis 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

coordinatively unsaturated Cr3c cations are generally considered to be vital for the catalytic activity, providing basic adsorption sites for NO and O2. As a start, the adsorption energetics and geometries of NO on Cr2O3(001) surface were explored. As shown in Fig. 2(a, b), NO on Cr2O3(001) prefers to adsorb via its N atom on the Cr3c in a bending configuration, corresponding to an adsorption energy of -0.79 eV. Firstly, we considered whether the adsorbed NO (NO*) could react with the lattice oxygen (the Mars-van Krevelen path). It was found that the reaction NO* +O3c → NO2* requires a barrier as high as 1.81 eV, and the desorption energy of NO2 is very high (1.96 eV), which is too high for the occurrence of reaction at room temperature, indicating that NO oxidation cannot proceed via the MvK path on Cr2O3(001) surface. Secondly, we resorted to the co-adsorption and reaction of NO and O2 on Cr2O3(001) via the Langmuir-Hinshelwood mechanism (LH), in which an important issue is to determine the adsorption manner of O2 and to identify whether the adsorbed O2 could dissociate into O atoms that then oxidize NO. As illustrated in Fig. 2(c, d), O2 can molecularly adsorb on the unsaturated Cr3c with the O-O bond paralleled to the surface (Fig. 2d, π-Cr3c-adsorption) or alternatively in a tilted manner (Fig. 2c, σ-Cr3cadsorption), giving an adsorption energy of -0.90 eV and -0.95 eV, respectively, which shows that O2 can bind more strongly than NO molecule on the Cr2O3(001) surface. Moreover, the reactant O2 is typically richer than NO under the practical condition of NO oxidation, and therefore it can be expected that O2 could occupy the Cr3c site preferentially. It should be noted that the O2 dissociation is difficult from both the two adsorption configurations, with the dissociation barriers being above 1.53 eV in which the distance of the O-O bond is activated to 1.90 Å at the TS, implying the infeasibility of obtaining O atoms via oxygen dissociation at room temperature. In other words, the reactant O2 has to participate in the subsequent reactions in the molecular state.

6

ACS Paragon Plus Environment

Page 6 of 23

Page 7 of 23 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

ACS Catalysis

Fig. 2 Optimized adsorption structures of NO (a, b), molecular O2 (c, d), dissociated O2 (e), and NO2 (f-j) on the Cr2O3(001) surface. The value on the top right of each configuration is the corresponding adsorption energy. All energies are in eV.

As depicted in SCHEME 1, the proposed reaction paths for NO oxidation via NO and O2 co-adsorption on the Cr2O3(001) surface comprise six consecutive elementary steps. (i) the co-adsorption of NO and O2: NO + O2* → NO* + O2*; (ii) the reaction between the adsorbed NO* and O2*: NO* + O2* → NO2* + O*; (iii) the release of the first NO2: NO2* + O* → NO2 + O*; (iv) the co-adsorption of O atom with another NO: NO + O* → NO* + O*; (v) the formation of the second NO2: NO* + O* → NO2* + *; and (vi) the desorption of the second NO2: NO2* → NO2 + *, completing the catalytic cycle. It was worth noting that, with respect to the initial co-adsorption of NO and O2 molecule, they could adsorb at two adjacent Cr3c sites or alternatively at the same single-atom Cr3c site, owing to the high unsaturation of three-coordinated Cr3c. These two adsorption configurations lead to two different reaction routes, named as LH-1 and LH-2, respectively, and Fig. 3 illustrates their energy profiles (blue and green curves, respectively). SCHEME 1. Elementary steps of NO Oxidation based on the Langmuir-Hinshelwood path (blue for LH-1: co-adsorption of NO and O2 at adjacent reaction sites, green for LH-2: co-binding at the single Cr atom), and red for Eley-Rideal path on the Cr2O3(001) surface.

7

ACS Paragon Plus Environment

ACS Catalysis 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 23

O *O Cr

NO

Cr

LH-2

O2 NO Cr

Cr

NO

LH-1

ER N

O

O

O

Cr

O

*O

*N

Cr

Cr

O N*

*O

Cr

O O*

NO2

Cr

Cr O

O N

NO2

O

O*

N*

Cr

Cr

O*

O O N* Cr

Cr

O* Cr

Cr

NO2 O* NO

Cr

Cr

NO2

For LH-1 path, the adsorption energy of NO at Cr3c with O2* sitting at the nearestneighbored Cr3c is -0.78 eV (see IS1 in Fig. 3), and the as-formed co-adsorption species can overcome a low barrier of 0.34 eV to form NO2* and atomic O* (see FS1 in Fig. 3), which is exothermic by -2.21 eV. The optimized transition sate structure (TS1) was shown in Fig. 3, and one can see that the forming N-O bond is 2.287 Å, and the O-O bond is elongated to 1.338 Å from 1.305 Å. Afterwards, NO2* could desorb from Cr3c but is limited by large energy requirement as high as 1.30 eV. By comparison, LH-2 path starting with the NO and O2 co-adsorption at the single Cr3c atom seems more favorable. Firstly, as NO co-adsorbs with O2* at the same Cr3c, they would spontaneously transform to a four-membered-ring ON*OO* species (see IS2 in Fig. 3) without an apparent barrier; more importantly, this process is strongly exothermic by 1.19 eV, indicating the more stability of IS2 compared with IS1 (see Fig. 3). Secondly, the ON*OO* intermediate can readily break its O-O bond, resulting in the co-adsorbed NO2* and O*, which is slightly exothermic by -0.08 eV and has a low barrier of 0.39 eV (see TS2 in Fig. 3 for the transition state, showing an elongated O-O bond of 1.825 Å). Moreover, the release of NO2* becomes easier at an energy release of 0.42 eV owing to the competing effect of atomic O* (see FS2 in Fig. 3), comparing to the strongly endothermic progress via LH-1 path (1.30 eV). By comparing the energy profiles of such two LH pathways for NO oxidation (NO* +O2* → NO2 +O*; see Fig. 3), one can clearly see that LH-2 is more favorable energetically. However, both the LH paths cannot remove the left O* through reacting 8

ACS Paragon Plus Environment

Page 9 of 23 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

ACS Catalysis

with NO. On one hand, O* cannot react with another NO* via the LH-1 path as a result of the long distance between the adjacent Cr cations; on the other hand, it is not plausible that NO co-binds with O* at the same Cr3c to follow LH-2 path owning to the strong binding of O* (-1.56 eV). Therefore, an alternative pathway in which NO directly couples with the surface O*, as well as the O2*, from the gas phase (EleyRideal, ER path) was examined below.

Fig. 3 Energy profiles of NO oxidation on Cr2O3(001) from three pathways: LH-1 (blue), LH-2 (green) and ER route (red). The solid and dashed lines represent the total energy profiles including the ZPE corrections and the Gibbs free energy profiles (including mainly the entropy contribution of gaseous molecule) at 300 K, respectively. The optimized structures of some important intermediates and transition states of each path are illustrated, and the key structural parameters are labeled. All the length values are in Å.

As illustrated in SCHEME1 (see red line), the gaseous NO can spontaneously couple with the surface adsorbed O2* to form a new monodentate cis-ONOO* configuration (see IS3 in Fig. 3) when the NO molecule is arbitrarily put around the O2* at a distance of ~3.0 Å or longer; the whole process is exothermic by -1.29 eV. To confirm this barrierless process, we examined the energy profile of ONOO* formation by putting the gaseous NO at a series of heights from the adsorbed O2*. As depicted in Fig. 4 (blue line), one can see that the energy monotonously declines as the ON…OO bond distance reduces, implying that there is indeed no total energy barrier. Then, the intermediate cis-ONOO* can readily break its ONO-O bond and release the NO2 into the gas phase, with a rather low barrier of 0.08 eV and a reaction energy change of -0.40 eV. The corresponding transition state (TS3) is shown in Fig. 3, in which the ONO-O bond 9

ACS Paragon Plus Environment

ACS Catalysis 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

increases to 1.711 Å from the original 1.442 Å. These results demonstrate that the ER path could exist in the system. Specifically, comparing the ER path with LH-2 path from the energy profiles in Fig. 3, one can see that the identified cis-ONOO* species (IS3) is 0.10 eV more stable than the four-membered-ring one (IS2) in LH-2 path, and the cleavage barrier of the ONO-O bond in cis-ONOO* is lower than the latter one (0.08 eV vs 0.39 eV), which results in a smaller effective barrier for ONOO conversion ((IS3 → FS3) versus (IS2 → FS2)). Moreover, another merit of ER pathway at this stage is that NO2 can be directly released to the gas phase instead of the indirect energyconsuming NO2* desorption process in LH-2 path. In other word, we can speculate that the ER path could be a preferable route for the stage (NO +O2 +* → NO2 + O*, i.e. the formation process of the first NO2) in NO oxidation. Also, it was found that the gaseous NO can follow the ER route to readily couple with the left atomic O* from the above process, forming the second NO2* (a monodentate-ONO configuration, see Fig. 2i) exothermically by -1.11 eV. Similar to the case of the surface O2* with ER path, this is also a barrierless process, which is confirmed by scanning the potential energy profile of NO2* formation (see red line in Fig. 4). Furthermore, the natural population analysis has also been performed using the Periodic NBO software56 (see details in the Table S1-2) to understand these two barrierless processes of NO coupling with O* or O2*. It was found that either the surface O* or O2* on Cr2O3(001) carries an unpaired electron, being the radical-like species, which could readily couple with the open-shell NO molecule. In other word, it could be the radical-like chemistry that accounts for the origin of such two barrierless processes. To better understand the overall reaction behavior at the realistic temperature (T=300 K), the entropy effects of the reactant NO/O2/NO2 molecule (T∆SNO, T∆SO2 and T∆SNO2 are 0.66 eV, 0.64 eV and 0.75 eV, respectively) have been taken into consideration to estimate the free energy profile (see details in Calculation Method and Model). From the obtained free energy profiles (see dashed line in Fig. 3), one can see that the whole free energy is downhill for the formation of the first gaseous NO2 and the second NO2* intermediate in ER path, indicating that it is an easy process on the Cr2O3(001) surface. However, the adsorbed monodentate NO2* is so strongly bound to the Cr3c site (-1.56 eV), resulting in its difficult removal and eventually limiting the overall activity of Cr2O3(001) towards NO oxidation at room temperature.

10

ACS Paragon Plus Environment

Page 10 of 23

Page 11 of 23

0.00

-0.25

-0.50

Energy /eV

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

ACS Catalysis

-0.75

-1.00

-1.25

TS_ONO TS_ONOO

-1.50 4.0

3.5

3.0

2.5

2.0

1.5

1.0

Bond Distance

Fig. 4 Illustration of NO chemisorption processes on molecular (blue) and atomic (red) oxygen adsorbed on Cr2O3(001). The zero point (dashed) is aligned to the system energy when the NO molecule in the gas phase.

According to the above discussions, an optimal pathway of NO oxidation on Cr2O3(001) can be summarized as follows. Firstly, the NO reacts with the surface adsorbed O2* from the gas phase (ER path), leading to the formation of cis-ONOO* species without an evident barrier. Secondly, the cleavage of cis-ONO-O* bond results in the O* formation and the first NO2 desorbing into the gas phase, with a dissociation barrier of 0.08 eV. Thirdly, NO couples with the adsorbed atomic O* at the isolated Cr3c site spontaneously, resulting in the formation of the second NO2*. Finally, the NO2* desorbs, completing the catalytic cycle. The Cr2O3(001) surface can efficiently catalyze the first three steps, but the adsorbed NO2* binds with the highly unsaturated Cr3c sites so strongly that it becomes a poisoning species and limits the overall activity of Cr2O3(001) for NO oxidation at room temperature. 3.2 NO oxidation on the Cr2O3(012) surface. Considering that the (012) facet is another typically exposed surface,57 we further investigated the NO oxidation on the Cr2O3(012) surface to explore the structural effect on the activity of Cr2O3. As depicted in Fig. 1d, Cr2O3(012) exposes the threecoordinated lattice O3c and five-coordinated Cr5c, which constitute basic active sites; the distance of adjacent Cr5c-Cr5c on Cr2O3(012) is 3.559 Å, being relatively shorter than that of Cr3c-Cr3c (5.051 Å) on Cr2O3(001). By viewing the coordination environment of surface Cr site, it could be inferred that the NO2 adsorption could decrease at Cr5c, possibly giving rise in an improved activity for NO oxidation in 11

ACS Paragon Plus Environment

ACS Catalysis 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

comparison with Cr2O3(001).

Fig. 5 Optimized adsorption structures of NO (a), molecular O2 (b, e, f), dissociated O2 (g) and NO2 (c, d, h) on the Cr2O3(012) surface. The value on the top right of each configuration is the corresponding adsorption energy. All the energies are in eV.

As shown in Fig. 5, Cr2O3(012) has indeed a weaker adsorption ability on NO2 as expected, due to the increased coordination number of the surface Cr5c site. The most stable configuration of NO2* on the (012) surface is μ(N, O) nitrito (Fig. 5h) with an adsorption energy of -1.07 eV, which is much weaker than that on the (001) surface (1.56 eV). Especially, for the monodentate-ONO, which is the most stable case on the (001) surface, the corresponding binding energy on the (012) surface is only -0.77 eV (Fig. 5d). Thus, it is reasonable to suggest that NO2 adsorption might not be the limiting factor on the (012) surface. To evaluate the overall activity of Cr2O3(012), the adsorption of reactant NO and O2 and their possible interactions on Cr2O3(012) were investigated. It is found that NO tends to adsorb at Cr5c with the N-end with an adsorption energy of -0.75 eV (vs -0.32 eV with the O-end, see Fig. 5a). Then we explored the MvK path, and the formation of ONO3c was identified to be endothermic by 0.98 eV. Moreover, the desorption of ONO3c requires an energy as high as 2.47 eV, revealing that the NO oxidation cannot proceed via the MvK route on the (012) surface either. With respect to the adsorption of O2 molecule, it is found that O2 would preferably adsorb with a double σ-bond on Cr2O3(012) (Fig. 5f) at an adsorption energy of -0.65 eV, compared with the other two adsorption configurations, the single σ-bond (Fig. 5b, -0.41 eV) and π-bond (Fig. 5e, 0.31 eV), on the (001) surface. Besides, the dissociation of double σ-bond adsorbed O2 into O atoms is exothermic (-0.58 eV), but gives a barrier as high as 1.31 eV, which 12

ACS Paragon Plus Environment

Page 12 of 23

Page 13 of 23 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

ACS Catalysis

implies that the dissociation is difficult at ambient temperature. In light of the fact that the NO binding energy at one Cr5c site is better than O2 occupying two Cr5c sites (-0.75 eV vs -0.65 eV), it is worth noting that in reality O2 could exist in minority and prefer to be with the single σ-bond configuration to participate in the catalytic oxidation, owing to the competitive adsorption effects between NO* and O2*. At the adjacent Cr5c site around NO*, O2 can weakly co-adsorb with an adsorption energy of -0.41 eV, named as the IS1 state (see Fig. 6); IS1 can quickly transform into a more stable ON*OO* species (IS1’) without an apparent barrier. Subsequently, the formed ON*OO* can break its O-O bond by overcoming a barrier of 0.69 eV, yielding the intermediates ON*O and O*, which corresponds to a reaction enthalpy of -0.59 eV. The ON*O can readily desorb at an energy cost of 0.54 eV and form the first NO2 molecule. With respect to the remaining O* atom, it can readily react with another NO* adsorbed nearby to form the μ-(N, O)-nitrito NO2* (see Fig. 5h) because of the short Cr-Cr distance on the (012) surface. This progress is almost spontaneous; the barrier is rather low (0.10 eV, the forming O-N bond is 2.069 Å at the TS2), and the reaction is much exothermic by -0.76 eV. Finally, the second NO2 desorbs into the gas phase, completing the catalytic cycle. The free energy profile in Fig. 6 reveals that the successive conversion from the coadsorbed NO* and O2* into two NO2 molecules can proceed smoothly without high barriers. However, O2 adsorption is very weak (+0.23 eV for the adsorption free energy (∆Gads) including the entropy effect at room temperature), specially due to the competitive effect resulting from the stronger adsorption of NO (∆Gads= -0.09 eV) to repel O2* occupation on Cr2O3(012). Accordingly, the low coverage of O2* can be expected, which would lead to a limitation on the overall activity of Cr2O3(012). From the free energy profile (see dashed lines in Fig. 6), one could further ascertain that the O2 activation would be the rate determining step in the whole progress of NO oxidation on the Cr2O3(012) surface. In addition, it is worth mentioning that, due to only a small amount of O2 could adsorb on the Cr2O3(012) surface, the subsequent Eley-Rideal process through the gaseous NO random collision with the surface O2* can be expected to hardly occur. Also, on the (012) surface, it is unlikely that NO and O2 co-adsorb at the single Cr5c atom, which rules out the possibility of LH-2 process as does on the (001) surface.

13

ACS Paragon Plus Environment

ACS Catalysis 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

Fig. 6 Reaction profile of NO oxidation by O2 on the Cr2O3(012) surface, following the LangmiurHinshelwood mechanism. The solid and dashed lines represent the total energy profiles including the ZPE corrections and the Gibbs free energy profiles (including mainly the entropy contribution of gaseous molecule) at 300 K, respectively. The optimized structures of some key intermediates and transition states of each path are illustrated, with the important bond length information. All the values are in Å.

Now, we are in a position to comparatively discuss the NO oxidation activity on the (001) and (012) surfaces. In brief, the (001) surface exposing the isolated Cr3c site can efficiently adsorb and activate the reactant O2, and NO can favorably attack O2* and trigger the subsequent reactions through the ER path. However, the overall reaction would be depressed by the strongly adsorbed NO2* species on Cr3c. In contrast, as a result of the surface exposed Cr5c site, the (012) surface has a weaker adsorption ability, which can guarantee the efficient desorption of the formed NO2 species; however, in the LH path the reactant O2 binds too weak at Cr5c and thus results in a low coverage, limiting the overall conversion rate of NO oxidation. Therefore, it is the structural coordination feature of Cr3c versus C5c that determines their different activity constraints. Also, the different adsorption ability of Cr3c and Cr5c site can be rationalized from their energy levels of 3d electron. As the projected density of states show in Fig. 7, the Cr 3d-electrons dominate the valence electron maximum (VBM) for both the (001) and (012) surfaces, and the 3d electrons of Cr3c on (001) are more localized near the Fermi level (and thus a higher reactivity) compared with that of Cr5c on (012), verified by a higher d-band center of -1.23 eV and -2.20 eV for the 3d majority electron, 14

ACS Paragon Plus Environment

Page 14 of 23

Page 15 of 23

respectively. 4

Projected Density of State

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

ACS Catalysis

2

0

-2

Cr3c@Cr2O3(001) Cr5c@Cr2O3(012)

-4

-4

-2

0

2

4

E-Ef

Fig. 7 Projected density of states (DOS) for the Cr 3d on the Cr2O3(001) (blue) and (012) surfaces (red) with respect to the Fermi levels.

3.3 Synergetic mechanism between Cr2O3 and CrO3. Based on the above results, to boost the catalytic activity of Cr2O3 for lowtemperature NO oxidation, one possible strategy is to increase the surface O2 coverage on Cr2O3(012), or to weaken the binding strength of NO2 on the Cr2O3(001) surface. Noteworthily, under realistic condition of NO oxidation, the presence of the excess oxygen could potentially oxidize Cr2O3 into CrO3 species thermodynamically. Also, some evidence that the Cr6+ species could promote or even dominate the overall activity of Cr2O3 catalysts has been given.22,23,28-30 Here, we computationally explore the possible synergetic mechanism between Cr2O3(012) and CrO3, aiming to enhance the O2 adsorption on Cr2O3(012) and thus improve the whole catalytic performance.

Fig. 8 (a) Bulk structure of CrO3 oxide, which shows a chain-stacked configuration. (b) The proposed model structure of mono-chain-CrO3/Cr2O3. (c) The O2 adsorption on the interface site of Cr6+/Cr3+. Green, blue, red (pink) indicate the Cr6+, Cr3+ and O, respectively, while yellow gives the adsorbed O2.

15

ACS Paragon Plus Environment

ACS Catalysis 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 23

In order to obtain a quantitative estimation of the presence of CrO3, a simple phase diagram was calculated by considering the thermodynamic equilibrium between the two states, Cr2O3 and CrO3. The crystal structure of CrO3 was illustrated in Figure 8a, which is a stacked configuration of one-dimensional CrO3-chains through the van der Waals weak interaction.58-59 At a given temperature and O2 partial pressure, the Gibbs free energy change from the Cr2O3 to CrO3 can be written as:

1 3 bulk ) GCrO (T , pO2 ) − GCrbulk (T , pO2 ) − µO2 (T , pO2 ) ∆G (T , pO2= 3 2 O3 2 4

(1)

bulk bulk where GCrO (T , pO2 ) are the Gibbs free energies for CrO3 and (T , pO2 ) and GCr 3 2 O3

Cr2O3, respectively. µO2 (T , pO2 ) is the chemical potential of gas-phase O2 at the given temperature T and partial pressure pO2 . The chemical potential µO2 (T , pO2 ) was calculated based on the thermodynamic equation:

p

O µO (T= , pO ) µO (T , p°) + RT ln = µO (0K, p°) − TSO (T , p°) + RT ln 2

2

2

2



2

2

pO2 p°

(2)

Where µO2 (0K, p°) was calculated from DFT including the zero-point energy, and the entropies were taken from the experimental data.60 When ∆G (T , pO2 ) < 0 , the oxidation formation is favored. Therefore, we can obtain the phase diagram of Cr2O3/CrO3, illustrated in Fig. 9 (the black solid line). We can see that under the typical experimental condition (T=300 K, pO2 p° =0.2 ), the Cr2O3 are not inclined to form CrO3 thermodynamically. Nevertheless, owing to the weakly stacked configuration of CrO3, we further examined the possibility of the CrO3 mono-chain (maintaining the local structure of CrO359 without destroying the evident chemical bond) supported by Cr2O3(012) (denoted as mc-CrO3/Cr2O3(012)), which could be regarded as the initially evolved CrO3 species on Cr2O3. The optimized mc-CrO3/Cr2O3(012) interface structure is shown in Fig. 8b. Thermodynamically, it is interestingly found that such an interface structure can be favorably formed under the practical O2-rich condition (Cr2O3 + O2 → mc-CrO3/Cr2O3(012)), as shown by the phase diagram in Fig. 9, implying that mcCrO3/Cr2O3(012) could be thermodynamically self-produced on Cr2O3(012). At the mc-CrO3/Cr2O3(012) interface, our calculations show that the reactant NO can readily couple with the terminal oxygen of CrO3 and leave a vacant Cr6+ site, which 16

ACS Paragon Plus Environment

Page 17 of 23

can, together with the original Cr5c on Cr2O3(012), constitute a new active site Cr6+@Cr3+. Significantly, this interface Cr6+@Cr3+ site can adsorb O2 molecule with an improved adsorption energy of -0.99 eV (see the adsorption configuration in Fig. 8c), which may provide an alternative pathway to supply O2 species on Cr2O3(012). Also, at this interface site, the product NO2 would adsorb moderately with an adsorption energy of -1.09 eV, implying the relatively easy desorption once generated, i.e. no evident poisoning effect. Therefore, an enhanced activity for low-temperature NO oxidation catalyzed by Cr2O3(012) could be expected at the mc-CrO3/Cr2O3(012) interface. On the other hand, our model could also provide a possible explanation for the experimental XPS observation that Cr6+ could exist and improve the catalytic activity of Cr2O3. Additionally, such a supported mono-chain CrO3 itself may also facilitate the NO oxidation process as an individual site, which is under way in our follow-up work. 10

5

CrO3 log(pO2/p0)

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

ACS Catalysis

mc-CrO3/Cr2O3 0

-5

Cr2O3 -10

0

100

200

300

400

500

600

Temperature (K)

700

800

900

1000

Fig. 9 Phase diagram of Cr oxides against partial pressures of O2 ( pO2 ) and temperatures T. The black line represents thermodynamic boundaries between Cr2O3 and CrO3, and the red one for that between Cr2O3 and mc-CrO3/Cr2O3. The dashed lines show the typical experimental condition.

4.

Conclusion In summary, the reaction mechanism of NO oxidation catalyzed by two common

Cr2O3 facets, i.e. (001) and (012), at room temperature has been studied, aiming at shedding light on the catalytic ability of Cr2O3 and the inherent activity-limiting factors. It was found that the (001) termination with the highly unsaturated Cr3c exposed exhibits strong adsorption ability and the overall activity would be limited by the strong 17

ACS Paragon Plus Environment

ACS Catalysis 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

adsorption of NO2 species (poisoning effect), and the whole reaction tends to follow the Eley-Rideal mechanism. In contrast, on the (012) termination which exposes the five-coordinated Cr5c, the reactant O2 cannot be efficiently adsorbed and activated, implying the limited catalytic performance for NO oxidation, which are in good accordance with the experimental findings. To boost the catalytic activity of Cr2O3 catalyst, we proposed an interface structure by depositing a mono-chain CrO3 on Cr2O3(012), which could be self-evolved on Cr2O3(012) thermodynamically under the typical reaction condition. Significantly, it has been found that the interface Cr6+@Cr3+ sites could show a much improved adsorption ability of O2 and can enhance the activity of Cr2O3(012). This study presents a systematic picture on the catalytic activities of Cr2O3 on NO oxidation at room temperature, and the provided insight could facilitate the design of more efficient Cr-based catalysts. ASSOCIATED CONTENT *Supporting Information The Supporting Information is available free of charge via the Internet at http://pubs.acs.org. The natural population analysis of the NO coupling with O* or O2* on Cr2O3(001) surface; summary table of the optimized structures along with their total energy and the imaginary frequency for each TS. AUTHOR INFORMATION Corresponding Authors *H.W.: e-mail, [email protected]; tel, +86-21-64253453. *P.H.: e-mail, [email protected] Notes The authors declare no competing financial interest. ACKNOWLEDGEMENT This project was supported by the National Natural Science Foundation of China (21333003, 21622305), Young Elite Scientist Sponsorship Program by CAST (YESS20150131), The Shanghai Shuguang scholar program (17SG30), and the Fundamental Research Funds for the Central Universities (WJ1616007).

REFERENCES 1.

Roy, S.; Hegde, M. S.; Madras, G., Catalysis for NOx Abatement. Appl. Energy 2009, 86, 2283-

2297. 18

ACS Paragon Plus Environment

Page 18 of 23

Page 19 of 23 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

ACS Catalysis

2.

Amiridis, M. D.; Zhang, T.; Farrauto, R. J., Selective Catalytic Reduction of Nitric Oxide by

Hydrocarbons. Appl. Catal., B 1996, 10, 203-227. 3.

Ciardelli, C.; Nova, I.; Tronconi, E.; Chatterjee, D.; Bandl-Konrad, B., A "Nitrate Route" for the

Low Temperature "Fast SCR" Reaction Over a V2O5-WO3/TiO2 Commercial Catalyst. Chem. Commun. 2004, 0, 2718-2719. 4.

Mao, Y.; Wang, Z.; Wang, H.-F.; Hu, P., Understanding Catalytic Reactions over Zeolites: A Density

Functional Theory Study of Selective Catalytic Reduction of NOx by NH3 over Cu-SAPO-34. ACS Catal. 2016, 6, 7882-7891. 5.

Fridell, E.; Skoglundh, M.; Westerberg, B.; Johansson, S.; Smedler, G., NOx Storage in Barium-

Containing Catalysts. J. Catal. 1999, 183, 196-209. 6.

Sakamoto, Y.; Okumura, K.; Kizaki, Y.; Matsunaga, S.; Takahashi, N.; Shinjoh, H., Adsorption and

Desorption Analysis of NOx and SOx on a Pt/Ba Thin Film Model Catalyst. J. Catal. 2006, 238, 361-368. 7.

Després, J.; Elsener, M.; Koebel, M.; Kröcher, O.; Schnyder, B.; Wokaun, A., Catalytic Oxidation

of Nitrogen Monoxide over Pt/SiO2. Appl. Catal., B 2004, 50, 73-82. 8.

Bhatia, D.; McCabe, R. W.; Harold, M. P.; Balakotaiah, V., Experimental and Kinetic Study of NO

Oxidation on Model Pt Catalysts. J. Catal. 2009, 266, 106-119. 9.

Wu, C.; Schmidt, D. J.; Wolverton, C.; Schneider, W. F., Accurate Coverage-Dependence

Incorporated into First-Principles Kinetic Models: Catalytic NO Oxidation on Pt (111). J. Catal. 2012, 286, 88-94. 10. Bray, J. M.; Schneider, W. F., First-Principles Analysis of Structure Sensitivity in NO Oxidation on Pt. ACS Catal. 2015, 5, 1087-1099. 11. Kim, C. H.; Qi, G.; Dahlberg, K.; Li, W., Strontium-Doped Perovskites Rival Platinum Catalysts for Treating NOx in Simulated Diesel Exhaust. Science 2010, 327, 1624-1627. 12. Parks, J. E., Less Costly Catalysts for Controlling Engine Emissions. Science 2010, 327, 1584-1585. 13. Wang, W.; McCool, G.; Kapur, N.; Yuan, G.; Shan, B.; Nguyen, M.; Graham, U. M.; Davis, B. H.; Jacobs, G.; Cho, K.; Hao, X., Mixed-phase oxide catalyst based on Mn-mullite (Sm, Gd)Mn2O5 for NO oxidation in diesel exhaust. Science 2012, 337, 832-5. 14. Choi, S. O.; Penninger, M.; Kim, C. H.; Schneider, W. F.; Thompson, L. T., Experimental and Computational Investigation of Effect of Sr on NO Oxidation and Oxygen Exchange for La1–xSrxCoO3 Perovskite Catalysts. ACS Catal. 2013, 3, 2719-2728. 15. Sousa, J. P. S.; Pereira, M. F. R.; Figueiredo, J. L., Catalytic Oxidation of NO to NO2 on N-doped Activated Carbons. Catal. Today 2011, 176, 383-387. 16. Zhang, Z.; Atkinson, J. D.; Jiang, B.; Rood, M. J.; Yan, Z., Nitric Oxide Oxidation Catalyzed by Microporous Activated Carbon Fiber Cloth: An Updated Reaction Mechanism. Appl. Catal., B 2014, 148-149, 573-581.

19

ACS Paragon Plus Environment

ACS Catalysis 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

17. Thirupathi, B.; Smirniotis, P. G., Co-doping a Metal (Cr, Fe, Co, Ni, Cu, Zn, Ce, and Zr) on Mn/TiO2 Catalyst and Its Effect on the Selective Reduction of NO with NH3 at Low-temperatures. Appl. Catal., B 2011, 110, 195-206. 18. Chen, Z.; Yang, Q.; Li, H.; Li, X.; Wang, L.; Chi Tsang, S., Cr–MnOx Mixed-oxide Catalysts for Selective Catalytic Reduction of NOx with NH3 at Low Temperature. J. Catal. 2010, 276, 56-65. 19. Liu, H.; Wei, L.; Yue, R.; Chen, Y., CrOx–CeO2 Binary Oxide as a Superior Catalyst for NO Reduction with NH3 at Low Temperature in Presence of CO. Catal. Commun. 2010, 11, 829-833. 20. Cai, W.; Zhong, Q.; Zhang, S.; Zhao, W., Fractional-hydrolysis-driven formation of nonuniform dopant concentration catalyst nanoparticles of Cr/CexZr1−xO2 and their catalysis in oxidation of NO. Chem. Eng. J. 2014, 236, 223-232. 21. Cai, W.; Zhong, Q.; Zhao, W., Solvent Effects on Formation of Cr-Doped Ce0.2Zr0.8O2 Synthesized with Cinnamic Acid and Their Catalysis in Oxidation of NO. Chem. Eng. J. 2014, 246, 328-336. 22. Liu, S.; Zhang, M.; Huang, Y.; Zhao, K.; Gao, Z.; Wu, M.; Dong, Y.; Wang, T.; Shi, J.; He, D., A Novel Chromic Oxide Catalyst for NO Oxidation at Ambient Temperature. RSC Adv. 2014, 4, 2918029186. 23. Zhong, L.; Cai, W.; Yu, Y.; Zhong, Q., Insights into Synergistic Effect of Chromium Oxides and Ceria Supported on Ti-PILC for NO Oxidation and Their Surface Species Study. Appl. Surf. Sci. 2015, 325, 52-63. 24. Zhong, L.; Cai, W.; Zhong, Q., Evaluation of Cerium Modification over Cr/Ti-PILC for NO Catalytic Oxidation and Their Mechanism Study. RSC Adv. 2014, 4, 43529-43537. 25. Peña, D. A.; Uphade, B. S.; Smirniotis, P. G., TiO2-supported metal oxide catalysts for lowtemperature selective catalytic reduction of NO with NH3I. Evaluation and characterization of first row transition metals. J. Catal. 2004, 221, 421-431. 26. Shiba, K.; Hinode, H.; Wakihara, M., Catalytic Oxidation of NO to NO2 Over Cr/TiO2 and Cu/TiO2 Under Oxidizing Atmosphere. React. Kinet. Catal. Lett. 1996, 58, 133-137. 27. Cai, W.; Zhong, Q.; Zhang, S.; Zhang, J., Effects of Cr on the NO Oxidation over the Ceria–Zirconia Solid Solution. RSC Adv. 2013, 3, 7009-7015. 28. Zhong, L.; Yu, Y.; Cai, W.; Geng, X.; Zhong, Q., Structure-Activity Relationship of Cr/Ti-PILC Catalysts Using a Pre-Modification Method for NO Oxidation and Their Surface Species Study. Phys. Chem. Chem. Phys. 2015, 17, 15036-45. 29. Wang, A.; Guo, Y.; Gao, F.; Peden, C. H. F., Ambient-temperature NO oxidation over amorphous CrOx-ZrO2 mixed oxide catalysts: Significant promoting effect of ZrO2. Appl. Catal., B 2017, 202, 706714. 30. Wang, A.; Lin, B.; Zhang, H.; Engelhard, M. H.; Guo, Y.; Lu, G.; Peden, C. H. F.; Gao, F., Ambient temperature NO oxidation over Cr-based amorphous mixed oxide catalysts: effects from the second oxide

20

ACS Paragon Plus Environment

Page 20 of 23

Page 21 of 23 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

ACS Catalysis

components. Catal. Sci. Technol. 2017, 7, 2362-2370. 31. Bhattacharjee, D.; Mishra, B. K.; Chakrabartty, A. K.; Deka, R. C., Catalytic Activity of Anionic Au–Ag Dimer for Nitric Oxide Oxidation: A DFT Study. New J. Chem. 2015, 39, 2209-2216. 32. Wang, H.; Guo, Y.; Lu, G.; Hu, P., An understanding and implications of the coverage of surface free sites in heterogeneous catalysis. J. Chem. Phys. 2009, 130, 224701. 33. Wang, H.-F.; Guo, Y.-L.; Lu, G.; Hu, P., NO Oxidation on Platinum Group Metals Oxides: First Principles Calculations Combined with Microkinetic Analysis. J. Phys. Chem. C 2009, 113, 18746-18752. 34. Sun, B.-Z.; Xu, X.-L.; Chen, W.-K.; Dong, L.-H., Theoretical Insights into the Reaction Mechanisms of NO Oxidation Catalyzed by Cu2O(111). Appl. Surf. Sci. 2014, 316, 416-423. 35. Xu, G.; Zhang, L.; He, C.; Ma, D.; Lu, Z., Adsorption and Oxidation of NO on Various SnO2(110) Surfaces: A Density Functional Theory Study. Sens. Actuators, B 2015, 221, 717-722. 36. Hou, M.; Cen, W.; Zhang, H.; Liu, J.; Yin, H.; Wei, F., Adsorption and Oxidation of NO on Graphene Oxides: A Dispersion Corrected Density Functional Theory Investigation. Appl. Surf. Sci. 2015, 339, 5561. 37. Verma, A. A.; Bates, S. A.; Anggara, T.; Paolucci, C.; Parekh, A. A.; Kamasamudram, K.; Yezerets, A.; Miller, J. T.; Delgass, W. N.; Schneider, W. F.; Ribeiro, F. H., NO oxidation: A probe reaction on CuSSZ-13. J. Catal. 2014, 312, 179-190. 38. 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. 39. Kresse, G.; Furthmüller, J., Efficiency of Ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15-50. 40. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865-3868. 41. Kresse, G.; Joubert, D., From Ultrasoft Pseudopotentials to the Projector Augmented-wave Method. Phys. Rev. B 1999, 59, 1758-1775. 42. Wang, H.-F.; Kavanagh, R.; Guo, Y.-L.; Guo, Y.; Lu, G.-Z.; Hu, P., Structural Origin: Water Deactivates Metal Oxides to CO Oxidation and Promotes Low-Temperature CO Oxidation with Metals. Angew. Chem., Int. Ed. 2012, 51, 6657-6661. 43. Wang, H.-F.; Kavanagh, R.; Guo, Y.-L.; Guo, Y.; Lu, G.; Hu, P., Origin of extraordinarily high catalytic activity of Co3O4 and its morphological chemistry for CO oxidation at low temperature. J. Catal. 2012, 296, 110-119. 44. Xia, Q.; Chen, Z.; Shao, Y.; Gong, X.; Wang, H.; Liu, X.; Parker, S. F.; Han, X.; Yang, S.; Wang, Y., Direct hydrodeoxygenation of raw woody biomass into liquid alkanes. Nat. Commun. 2016, 7, 11162. 45. Wang, H.-F.; Wang, D.; Liu, X.; Guo, Y.-L.; Lu, G.-Z.; Hu, P., Unexpected C–C Bond Cleavage Mechanism in Ethylene Combustion at Low Temperature: Origin and Implications. ACS Catal. 2016, 6,

21

ACS Paragon Plus Environment

ACS Catalysis 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 23

5393-5398. 46. Grimme, S.; Ehrlich, S.; Goerigk, L., Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456-1465. 47. Maldonado, F.; Novillo, C.; Stashans, A., Ab Initio Calculation of Chromium Oxide Containing Ti Dopant. Chem. Phys. 2012, 393, 148-152. 48. Lebreau, F.; Islam, M. M.; Diawara, B.; Marcus, P., Structural, Magnetic, Electronic, Defect, and Diffusion Properties of Cr2O3: A DFT+U Study. J. Phys. Chem. C 2014, 118, 18133-18145. 49. Rohrbach, A.; Hafner, J.; Kresse, G., Ab Initio Study of the (0001) Surfaces of Hematite and Chromia: Influence of Strong Electronic Correlations. Phys. Rev. B 2004, 70, 125426-1-17. 50. Shi, S.; Wysocki, A.; Belashchenko, K., Magnetism of Chromia from First-principles Calculations. Phys. Rev. B 2009, 79, 104404-1-7. 51. Mosey, N. J.; Liao, P.; Carter, E. A., Rotationally Invariant Ab Initio Evaluation of Coulomb and Exchange Parameters for DFT+U Calculations. J. Chem. Phys. 2008, 129, 014103-1-13. 52. Hu, W.; Lan, J.; Guo, Y.; Cao, X.-M.; Hu, P., Origin of Efficient Catalytic Combustion of Methane over Co3O4(110): Active Low-Coordination Lattice Oxygen and Cooperation of Multiple Active Sites. ACS Catal. 2016, 6, 5508-5519. 53. Wang, D.; Jiang, J.; Wang, H. F.; Hu, P.,Revealing the volcano-Shaped Acitivity Trend of Triiodide Reduction Reaction: A DFT Study Coupled with Microkinetic Analysis. ACS Catal. 2016, 6, 733-741. 54. Nørskov, J. K.; Bligaard, T.; Logadottir, A.; Kitchin, J. R.; Chen, J. G.; Pandelov, S.; Stimming, U., Trends in the Exchange Current for Hydrogen Evolution. J. Electrochem. Soc. 2005, 152, J23-J26. 55. Peng, C.; Wang, H. F.; Hu, P., Theoretical Insights into How the First C-C Bond Forms in the Methanol-to-Olefin Process Catalysed by HSAPO-34. Phys. Chem. Chem. Phys. 2016, 18, 14495-14502. 56. Dunnington, B. D.; Schmidt, J. R., Generalization of Natural Bond Orbital Analysis to Periodic Systems: Applications to Solids and Surfaces via Plane-Wave Density Functional Theory. J. Chem. Theory Comput. 2012, 8, 1902-1911. 57. York, S. C.; Abee, M. W.; Cox, D. F., α-Cr2O3(10 1� 2): Surface Characterization and Oxygen Adsorption. Surf. Sci. 1999, 437, 386-396.

58. Wooster., W. A.; Wooster., N., Crystal Structure of Chromium Trioxide. Nature 1931, 127, 782-782. ∞

59. Stephens, J. S.; Cruickshank, D. W. J., The Crystal Structure of (CrO3) . Acta Crystallogr. Sect. BStruct. Sci.Cryst. Eng. Mat. 1970, 26, 222-226. 60. Haynes, W. M., CRC Handbook of Chemistry and Physics, 93rd Edition. CRC Press, Boca Raton, FL: 2012.

TOC

22

ACS Paragon Plus Environment

Page 23 of 23

ACS Catalysis

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 23

ACS Paragon Plus Environment