Shape Evolution of Metal Nanoparticles in Binary Gas Environment

Technology, Shanghai Institute of Applied Physics, Chinese Academy of Sciences,. Shanghai ... Introduction. Metal nanoparticles (NPs) have drawn lots ...
0 downloads 11 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Article

Shape Evolution of Metal Nanoparticles in Binary Gas Environment Jun Meng, Beien Zhu, and Yi Gao J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b00052 • Publication Date (Web): 28 Feb 2018 Downloaded from http://pubs.acs.org on March 2, 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.

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

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

The Journal of Physical Chemistry

Shape Evolution of Metal Nanoparticles in Binary Gas Environment Jun Meng1,2, Beien Zhu*1 and Yi Gao*1 1Division

of Interfacial Water and Key Laboratory of Interfacial Physics and

Technology, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai, 201800 China. 2University

of Chinese Academy of Sciences, Beijing 100049, China

*Email: [email protected]; [email protected].

1 ACS Paragon Plus Environment

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

Page 2 of 24

Abstract

To control the shape and structure of metal nanoparticle (NP) is a crucial strategy to improve its catalytic properties, but the understanding and quantitative description of the structure reconstruction of the catalysts under reaction conditions has not been achieved. Previous works are mostly focused on the single gas conditions, which is apparently not the case in the real catalytic reactions. In this work, a multi-scale structure reconstruction model is established to describe the equilibrium structures of metal NPs in mixed gas environment quantitatively. Taking NO and CO reaction as a model system, the structures of the Pd, Pt and Rh NPs in a large range of temperature and pressure are fully presented. Moreover, we show the variation of 𝑃 : 𝑃

plays

the critical role in determining the structures and therefore the number of active sites of the NPs at certain conditions. This work provides an efficient model to guide the future experiments in the real reactions.

2 ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

Introduction Metal nanoparticles (NPs) have drawn lots of research interests because of their remarkable catalytic performance in a large range of chemical reactions. One important issue is the shape control of metal NPs, since the physical properties and catalytic performance of metal NPs mainly depend on their sizes, shapes and structures. 1-2 To improve the activity of metal NPs, many experimental procedures have been developed to prepare metal NPs with specific shapes.3-6 However, well-defined NPs may change their shape in the reactive environment,7-17 resulting in dramatically altered catalytic performance.18-24

Using

time-resolved

high-resolution

transmission

electron

microscopy, a periodic refacetting of the Pt NPs was observed in the CO oxidation reaction.12 Eren et al. have reported that the Cu(111) surface decomposed into clusters when exposed to CO gas, and the cluster-covered Cu(111) surface was very active in dissociating water.21 Thus, it is a general idea now that the study of the morphology change of metal NPs in reactive environment is a key step to understand and control the real catalytic properties. However, up to date, most of the experiments and theoretical studies are limited in the condition of single gas, while in a real catalysis the NPs are mostly exposed to mixed gas conditions. NO and CO as the major pollutants of car exhausts, cause great damage to our living environment. One of the most promising way to eliminate the exhaust toxic gases is the reduction of NO by CO.25-26 The three-way catalysts including Pd, Pt, and Rh have been used in the automobile industry to oxidize CO to O2 and reduce NO to N2 simultaneously.27-29 Both NO and CO can influence the equilibrium shape of metal 3 ACS Paragon Plus Environment

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

NPs.10, 12, 19-23, 30-31 As a typical instance, Newton et al. have observed the size and shape changes of Pd NPs during NO/CO cycling.23 It is also known that the different facets of metal nanoparticles exhibit quite different catalytic activities. Goodman and coworkers have investigated the NO+CO reaction on Pd(111) and (100) facets. 32 At temperatures from 525K-650K, 𝑃 /𝑃

from 1:16 to 16:1, and total pressure of 2-

17 Torr, they found out that Pd(111) always shows a higher activity than Pd(100) for both CO2 and N2O production, while Pd(100) gives a higher selectivity for N2 than Pd(111). Nakao et al. discovered that Pd(110) is more active than Pd(111) surface in the pressure range of 10-2~10-1 Torr.33 Thus, the understanding of the morphology change of the metal NPs under CO and NO would be critical for the tuning and controlling its catalytic activities for real applications. In this work we present a multi-scale structure reconstruction model (MSRM) to quantitatively describe the equilibrium structures of metal NPs larger than 3nm interacting with mixed gas. In our previous studies, the corresponding single gas model gave precise description of metal NPs under single gas environment.13-15 We have discovered that the introducing of water vapor change the shape of Cu NPs,13 which have been observed by the in situ TEM experiments.9 Besides, CO, NO, O2, or H2 also has influences on the equilibrium shape of metal NPs depending on the temperature and the pressure.14-15 In order to figuring out how the mixed gas could affect the shape of metal NPs, taking NO + CO reaction as a model system, the shape evolutions of Pd, Pt and Rh NPs under NO and CO environment are investigated. The results show that the structures of Pd, Pt and Rh NPs could be not only adjusted by the temperature, pressure 4 ACS Paragon Plus Environment

Page 4 of 24

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

The Journal of Physical Chemistry

but by the partial pressure ratio of NO to CO (𝑃 ⁄𝑃 ). Under certain conditions, the structures of the NPs are highly dependent on the 𝑃 ⁄𝑃

and largely differ from the

ones in pure NO or CO environment. The partial pressure of the mixed gas, together with the temperature and the total pressure, should be considered as crucial parameters in the rational design of nanocatalysts. Our model provides an effective and useful method for advancing the structure reconstruction study from single gas effects to mixed gas effects. Computational Details According to the Wulff rule, the equilibrium shape of a nanoparticle(NP) could be constructed if the surface tensions of the NP’s facets are known. 34 When the NP is exposed to the reaction environment, the surface tensions 𝛾 the interface tensions 𝛾

should be revised as

, considering the interactions between the NP surfaces and

the adsorbates.14 To obtain revised interface tension, Fowler-Guggenheim (F-G) adsorption isotherm35 is employed to deal with the localized monolayer adsorption with the lateral interactions = 𝑃𝐾 𝑒𝑥𝑝 (1) where 𝜃 is the coverage of the adsorbates on the NP surface, 𝑃 is the gas pressure, 𝐾 is the equilibrium constant, 𝑧 is the number of the nearest neighbor sites, 𝑤 represents the interaction energy between two adjacent adsorbates, 𝑅 is the gas constant, and 𝑇 is the temperature. When the NP interacts with binary gases, an extended competitive adsorption isotherm can be derived based on this isotherm. 36 5 ACS Paragon Plus Environment

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

(

)

(

)

= 𝑃 𝐾 𝑒𝑥𝑝

(

)

(

)

Page 6 of 24

(2) = 𝑃 𝐾 𝑒𝑥𝑝

(3) where the subscripts A and B correspond to the two component A and B, respectively. 𝑤 (𝑤 ) represents the lateral interaction between A(B) molecules. 𝑤

represents the

interaction energy between A molecule and B molecule.𝐾 and 𝐾 is the equilibrium constant of the adsorption process for gas A and B, respectively, which are defined as 𝐾 = 𝑒𝑥𝑝 −

_

_

_

𝐾 = 𝑒𝑥𝑝 −

_

_

_

(4)

(5) where the 𝐸 𝑆

_

/𝑆

_

_

/𝐸

and 𝑆

is the adsorption energy of a gas A/B molecule.

_

_

/𝑆

is the entropy of gas A/B before and after adsorption,

_

respectively. By solving the above equations, we can get the value of 𝜃 and 𝜃 . Then the interface tension 𝛾 𝛾

= 𝛾

+

could be obtained by the following equation (

_

( A

+ AB

))

(

_

( AB

+

B

))

(6) 𝐴

is the surface area per atom. In this work, the surface tension 𝛾

, the adsorption energy 𝐸

and the lateral

interaction 𝑤 are computed by spin-polarized density functional theory (DFT) calculations using the Vienna Ab-initio Simulation Package (VASP). Three low-index 6 ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

facets (100), (110) and (111) are considered in this study to build the atomic structures of the NPs. The Perdew, Burke, and Erzernhof functional (PBE) is used to calculate the surface tension 𝛾

. By optimizing a series of (1×1) slab models with increasing

layers, the accurate surface energy could be extracted. The surface energies per atom are presented in Table S1 in the supporting information. Revised gradient-corrected functional of Perdew, Burke and Erzernhof (RPBE) is used to get the adsorption energies 𝐸

and the lateral interactions 𝑤 of NO/CO on

Pd/Pt/Rh surfaces. The adsorption energy 𝐸

is obtained by a (4×4) slab model.

Combined with the calculation of a (1×1) slab model, the lateral interaction energy 𝑤 between the same kind of adsorbates can be derived in such way ( × )

𝑤=

( × )

(7) where 𝐸

( × )

represents the monomolecular adsorption energy, 𝐸

( × )

is

the adsorption energy at 1 monolayer adsorption. The adsorption sites are carefully tested in the (4×4) and (1×1) slab models and the most stable structures are used to calculate the lateral interaction. In order to obtain the lateral interaction between NO and CO molecules, denoted by 𝑤

, a (2×2) slab model with the coadsorption of NO and CO molecules is

established. NO and CO coadsorbed on the surfaces side by side and the total coverage is 1 monolayer. The 𝑤

can be derived by solving several equations. Details of

the calculation methods and the results of the adsorption energies 𝐸 interaction 𝑤

, 𝑤

and 𝑤

are displayed in the ESI file. 7

ACS Paragon Plus Environment

, lateral

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

Page 8 of 24

The slab models are separated with a 15 Å vacuum space along Z axis. Each model contains 6 atom layers and the bottom 2 layers are fixed. For all the DFT calculations, the cutoff energy of plane wave is 400 eV, the convergence of the electronic selfconsistent is 10-5 eV, and the force convergence criteria in a conjugate-gradient algorithm is 0.05 eV/Å respectively. The Kpoints are set as (2×2×1), (4×4×1) and (6×6×1) for the (4×4), (2×2) and (1×1) slab models, respectively. When the gas molecules adsorbed on the surfaces, the entropy 𝑆 as zero. The entropy of gas phase molecules 𝑆 𝑆 = 𝑆(𝑇, 𝑃 ) − 𝑅𝑙𝑛(

is assumed

is described by

)

(8) 𝑃 is 1 atm, and 𝑆(𝑇, 𝑃 ) is the entropy at 1 atm with different temperature. The 𝑆(𝑇, 𝑃 ) is fitted by data from the NIST-JANAF Thermochemical Tables37. Applying this multi-scale construction model, the equilibrium shape of Pd, Pt and Rh NPs under NO and CO environment are predicted. All the constructed NPs are around 10 nm, 28000 atoms. Results & Discussion The equilibrium structures of Pd, Pt and Rh NPs are first constructed with fixed partial pressure ratio ( 𝑃 /𝑃

= 1: 1) , which is a common 𝑃 /𝑃

in real

application due to the stoichiometry of the reaction 2NO + 2CO → N2 + 2CO2. The temperature is ranging from 200 K to 1200 K and the total pressure is ranging from 1 Pa to 10000 Pa. Fig. 1 shows the structure evolution of a Pd NP with respect to the temperature and the total pressure. Generally, at a constant total pressure, the shape of 8 ACS Paragon Plus Environment

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

The Journal of Physical Chemistry

the Pd NP changes from a rhombic dodecahedron (blue region) to a truncated octahedron (red region) with increasing temperature. At a constant temperature, the shape of the Pd NP changes from a truncated octahedron (red region) to rhombic dodecahedron (blue region) by increasing the total pressure. The dynamic response of Pd NP to temperature and pressure is attributed to the change of coverage of NO and CO molecules on each of the three low-index facets under different conditions.

Fig. 1

(a) Contour plot of the structures of a 10nm Pd NP changing with the

temperature and the total pressure (𝑷𝑵𝑶 /𝑷𝑪𝑶 = 𝟏: 𝟏). Different structures are distinguished by the structure index number from 1 to 6. The black area indicates the negative interface energies area. (b) The typical structures of the Pd NP labeled by the structure index number. (c) The fractions of (111), (100) and (110) facets of the typical Pd NPs. Taking 1000 Pa as an example, which is a normal condition for the NO+CO reaction to take place,29 the coverage of NO and CO on the Pd (100), (110) and (111) surfaces, and the corresponding interface tensions are shown in Fig. 2. When the temperature is higher than 850 K, the Pd NP almost exposes (111) and (100) facets 9 ACS Paragon Plus Environment

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

entirely (structure 6 in Fig. 1). The (111) facet, as the dominant surface, occupies 85% of the surface region. Since the temperature is so high that the gases could not adsorb on the NPs (Fig. 2), the corresponding morphology is thus the same as the structure in vacuum. There is a large morphology change of the Pd NP at the temperature ranging from 700K to 850 K. As the temperature decreases, more NO and CO molecules could adsorb steadily on Pd surfaces. NO is the dominant adsorbate compared to CO due to its stronger adsorption energy (Table S2), as shown in Fig. 2a. The coverage of NO increases rapidly on (110) facet, resulting in the rapid decline of the interface energy of (110) facet (Fig. 2b). The proportion of (110) facets keeps to enlarge and reaches 100% when the temperature is smaller than 700 K (from structure 6 to 1). It is very interesting to see that the interface tension becomes negative when the temperature is below 600 K. As discussed in our previous work and in the literatures,14, 31 when the interaction between the adsorbates and the surface is strong enough, the interface tension could be negative. This will break up the NP into smaller nanoclusters. This prediction is in consistent with the experimental observation that Pd atoms with lower coordination numbers induced by NO.23

10 ACS Paragon Plus Environment

Page 10 of 24

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

The Journal of Physical Chemistry

Fig. 2

(a) The coverages of NO and CO on Pd (100), (110) and (111) surfaces,

changing with the temperature. (b) The interface energies between Pd surfaces with NO and CO gases, changing with the temperature. The total pressure is 1000 Pa (𝑷𝑵𝑶 /𝑷𝑪𝑶 = 𝟏: 𝟏). Above we present the morphologies that the Pd NPs could have under different conditions at given 𝑃 /𝑃 . One interesting and important thing in this work is to see how the 𝑃 /𝑃

may alter the structure at given temperature and total pressure. The

experiments have shown the metal NPs display different catalysis selectivity under different partial pressure ratio of NO to CO.32 To investigate the effect of 𝑃 /𝑃

on

the structure of Pd NPs, the equilibrium structures of a 10 nm Pd NP are constructed with different 𝑃 /𝑃

at different conditions.

11 ACS Paragon Plus Environment

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

Fig. 3

The structures evolution of a 10nm Pd NP changing with the temperature

and 𝑷𝑵𝑶 /𝑷𝑪𝑶 ratio under a total pressure of (a) 10 Pa, (b) 100 Pa, (c) 1000 Pa and (d) 10000 Pa. The structure evolution of Pd NP under pure CO and pure NO environment are displayed at the bottom and the top. The structures are represented by the structure index number. The typical structures of the Pd NP labeled by the structure index number have been shown in Fig. 1(b). The results of pure CO and pure NO are also provided to make comparison. From Fig. 3 we can see an overall conclusion that whether the variation of the 𝑃 /𝑃 would induce the structure reconstruction depends on the structures of the NP in pure CO and pure NO environment. If the structure of the NP in pure NO is the same as the one in pure CO environment, then it is the same for all the 𝑃 /𝑃 . There is no difference between the structure of the NP in single gas and in mixed gas. As we can 12 ACS Paragon Plus Environment

Page 12 of 24

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

The Journal of Physical Chemistry

see there are “solid colour” regions in Fig. 3. On the other hand, the 𝑃 /𝑃

has

significant effect on the structure of the NP if the structure in the pure NO environment is different to that in the pure CO environment. The “gradient colour” regions can be found in the middle of each graph in Fig. 3. As an example, at 10 Pa and 580K, the Pd NP is enclosed by 15% (100) and 85% (111) facets labelled by red colour (structure 6 in Fig. 1-b) in the pure CO environment, while under the same condition, it is enclosed by 88% (110) and 12% (111) facets labelled by cyan colour (structure 2 in Fig. 1-b) in the pure NO environment. Since NO has greater adsorption energies on Pd surfaces than CO, under this condition it has greater effects to decrease the interface energies, especially for the (110) facet. That is why the Pd NP has a large (110) facet in pure NO environment. When the 𝑃 /𝑃

increases from 1:10 to 10:1, there is a transition

procedure that the structure of the Pd NP gradual changes from the shape in pure CO to the shape in pure NO. In this transition procedure, the proportion of (110) facet of the Pd NP keeps growing, meanwhile the (111) and (100) facet keeps declining. So it is very interested to see that the NP goes through a complete red-yellow-light greendark green-cyan structure evolution. Apparently, the structure of the NP could be very different at different 𝑃 /𝑃 . Since different facet of the NP has different activity and selectivity, the change of the facet fraction caused by the 𝑃 /𝑃

also influences the NP’s catalytic property.

Notice that the boundaries between the “solid color” regions and the “gradient color” regions shift with temperature and pressure. The 𝑃 /𝑃

, together with the

temperature and pressure should be taken into consider as critical parameters in the 13 ACS Paragon Plus Environment

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

rational design of nanocatalysts. It actually provides more possible ways to tailor the desired structure and property of the catalyst. In the case of CO+NO reaction, the Pd(111) surface is the surface attracting lots of interests.28, 32, 38 If one wants to have more (111) facets when using Pd NPs as catalysts in this reaction, the structure 3 and 4 of the NP in Fig. 1-b are good candidates. The structure 5 and 6 are not good choice although they have the largest (111) facet fractions since the adsorption of CO and NO need to be weak enough to have those structures. According to our work, the structure 3 and 4 can be achieved by different ways. It could be done by controlling the temperature and total pressure at fixed 𝑃 /𝑃 tuning the 𝑃 /𝑃

(e.g. 1:1). It could also be done by

at the temperature and total pressure beneficial to the reaction.

Besides the structure, the 𝑃

/𝑃

could also influence the size of the NP. As we

discussed above, the black region in our contour plots representing the dispersion of the NP induced by the negative interface tension. In Fig. 3 it can be seen that under certain condition (e.g. 700 K, 10000 Pa) the NP transfers from a blue region to a black region with the increased 𝑃 /𝑃 . This indicates the NP disperses from bigger one to smaller ones. A similar phenomenon have been reported by Newton et al. in their experimental study of the size and shape changes of Pd NPs during NO/CO cycling. 23 At 673 K, the Pd NP was found to sinter in pure CO environment and to disperse in mixed CO+NO and pure NO environment. Their observations could be well understood by our results.

14 ACS Paragon Plus Environment

Page 14 of 24

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

The Journal of Physical Chemistry

Fig. 4

(a) Contour plot of the structures of a 10nm Pt NP changing with the

temperature and the total pressure (𝑷𝑵𝑶 /𝑷𝑪𝑶 = 𝟏: 𝟏). (b) The structures evolution of the Pt NP changing with the temperature and 𝑷𝑵𝑶 /𝑷𝑪𝑶 ratio under a total pressure of 10 Pa. Different structures are distinguished by the structure index number from 1 to 6 in (a) and (b). (c) Typical structures of Pt NP labeled by the structure index number. (d) Contour plot of the structures of a 10nm Rh NP changing with the temperature and the total pressure (𝑷𝑵𝑶 /𝑷𝑪𝑶 = 𝟏: 𝟏). (e) The structures evolution of the Rh NP changing with the temperature and 𝑷𝑵𝑶 /𝑷𝑪𝑶 ratio under a total pressure of 10 Pa. Different structures are distinguished by the structure index number from 1 to 6 in (d) and (e). (f) Typical structures of Rh NP labelled by the structure index number.

15 ACS Paragon Plus Environment

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

For Pt and Rh NPs, similar transformation trends have been observed. The structure of Pt and Rh NPs under NO and CO environment are shown in Fig. 4. Although the general evolution trends of Pt and Rh NP structure are consistent with Pd NP, there are still some differences. DFT calculations shows that the adsorption energies and the lateral interactions of NO and CO are different on Pd, Pt and Rh surfaces (Table S2). There is a great difference between the shape of Rh NPs in pure CO and mixed gas condition, because NO has much larger adsorption energies than CO on Rh (Table S2). The results show that the dispersion of Rh NPs from CO to mixed gas condition is easier than Pd and Pt. Meanwhile, there is a considerable region of negative surface energies for Rh NP under NO and CO environment, which means nanoclusters could be easily formed on Rh surfaces under NO and CO conditions. The changing of the facet ratio for Pd, Pt and Rh are consistent in general, except that the structure 2 of Pt NP is enclosed by (110) and (100) surfaces, Pd NP is enclosed by (110) and (111) surfaces and Rh NP is enclosed by (110), (111) and (100) facets. The effect of the 𝑃 /𝑃

on the structures are also investigated at 10 Pa shown in Fig (4-b,e)

(Results of other pressures are also presented in the ESI file). The results show that the partial pressure ratio of the mixed gas has a general effect on the structure of different metal NP. It should be mentioned the work of our MSRM model is not limited in NO and CO environment, but a general model suitable for other systems. Further application of this model to other binary or even more complex systems will be done in our future works. 16 ACS Paragon Plus Environment

Page 16 of 24

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

The Journal of Physical Chemistry

Conclusions In summary, a multi-scale structure reconstruction model (MSRM) is proposed to predict the equilibrium structures of metal NPs under binary gas environment. Taking NO and CO reaction as a model system, the shape evolution of Pd, Pt and Rh NPs under NO and CO binary gas environment is investigated. With a fixed partial pressure ratio of 𝑃

/𝑃

= 1 , the shape evolution of Pd, Pt and Rh NPs in a large range of

temperature and pressure is fully presented, which changes between a rhombic dodecahedron and a truncated octahedron. More important, we show the equilibrium structures could be largely alerted not only by the temperature and total pressure but also by tuning the partial pressure ratio of NO to CO. It is very necessary to take into consider the effect of the partial pressure ratio of the mixed gas when tuning the structure and catalytic property of the NPs. On the other hand, it offers more possible ways to design the catalysts. Our study shows the importance to study the structure reconstruction under mixed gas condition and also paves the way for the future studies in this area by offering a quantitative model.

Supporting Information. The following file is available free of charge. Computational details of the lateral interaction energies. Data of surface tensions, surface area, adsorption energies, and lateral interactions. Structure evolution of 10nm Pt (Rh) NPs change with the temperature and the 𝑃 /𝑃

ratio with a total pressure of 100Pa, 1000Pa and

10000Pa. 17 ACS Paragon Plus Environment

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

Notes The authors declare no competing financial interests. Acknowledgments This work is supported by National Natural Science Foundation of China (11604357, 21773287, 11574340). B.Z. also thanks for the development fund for Shanghai talents (Y439011011), and Natural Science Foundation of Shanghai (16ZR1443200). The authors also thank the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase) under Grant No. U1501501. The computational resources utilized in this research were provided by Shanghai Supercomputer Center, National Supercomputer Centers in Tianjin, Shenzhen.

18 ACS Paragon Plus Environment

Page 18 of 24

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

The Journal of Physical Chemistry

Reference 1.

Zalineeva, A.; Baranton, S.; Coutanceau, C.; Jerkiewicz, G. Octahedral Palladium

Nanoparticles as Excellent Hosts for Electrochemically Adsorbed and Absorbed Hydrogen. Sci. Adv. 2017, 3, e1600542. 2.

Cuenya, B. R. Synthesis and Catalytic Properties of Metal Nanoparticles: Size,

Shape, Support, Composition, and Oxidation State Effects. Thin Solid Films 2010, 518, 3127-3150. 3.

Zhang, L.; Roling, L. T.; Wang, X.; Vara, M.; Chi, M.; Liu, J.; Choi, S.-I.; Park, J.;

Herron, J. A.; Xie, Z., et al. Platinum-Based Nanocages with Subnanometer-Thick Walls and Well-Defined, Controllable Facets. Science 2015, 349, 412-416. 4.

Sun, Y.; Xia, Y. Shape-Controlled Synthesis of Gold and Silver Nanoparticles.

Science 2002, 298, 2176-2179. 5.

Ahmadi, T. S.; Wang, Z. L.; Green, T. C.; Henglein, A.; El-Sayed, M. A. Shape-

Controlled Synthesis of Colloidal Platinum Nanoparticles. Science 1996, 272, 19241925. 6.

Lee, K.; Kim, M.; Kim, H. Catalytic Nanoparticles Being Facet-Controlled. J.

Mater. Chem. 2010, 20, 3791-3798. 7.

Narayan, T. C.; Baldi, A.; Koh, A. L.; Sinclair, R.; Dionne, J. A. Reconstructing

Solute-Induced Phase Transformations within Individual Nanocrystals. Nat. Mater. 2016, 15, 768-774. 8.

Baldi, A.; Narayan, T. C.; Koh, A. L.; Dionne, J. A. In Situ Detection of Hydrogen-

Induced Phase Transitions in Individual Palladium Nanocrystals. Nat. Mater. 2014, 13, 19 ACS Paragon Plus Environment

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

1143-1148. 9.

Hansen, P. L.; Wagner, J. B.; Helveg, S.; Rostrup-Nielsen, J. R.; Clausen, B. S.;

Topsøe, H. Atom-Resolved Imaging of Dynamic Shape Changes in Supported Copper Nanocrystals. Science 2002, 295, 2053-2055. 10. Tao, F.; Dag, S.; Wang, L.-W.; Liu, Z.; Butcher, D. R.; Bluhm, H.; Salmeron, M.; Somorjai, G. A. Break-up of Stepped Platinum Catalyst Surfaces by High CO Coverage. Science 2010, 327, 850-853. 11. Yoshida, H.; Kuwauchi, Y.; Jinschek, J. R.; Sun, K.; Tanaka, S.; Kohyama, M.; Shimada, S.; Haruta, M.; Takeda, S. Visualizing Gas Molecules Interacting with Supported Nanoparticulate Catalysts at Reaction Conditions. Science 2012, 335, 317319. 12. Vendelbo, S. B.; Elkjær, C. F.; Falsig, H.; Puspitasari, I.; Dona, P.; Mele, L.; Morana, B.; Nelissen, B. J.; van Rijn, R.; Creemer, J. F., et al. Visualization of Oscillatory Behaviour of Pt Nanoparticles Catalysing CO Oxidation. Nat. Mater. 2014, 13, 884890. 13. Zhu, B.; Xu, Z.; Wang, C.; Gao, Y. Shape Evolution of Metal Nanoparticles in Water Vapor Environment. Nano Lett. 2016, 16, 2628-2632. 14. Zhu, B.; Meng, J.; Gao, Y. Equilibrium Shape of Metal Nanoparticles under Reactive Gas Conditions. J. Phys. Chem. C 2017, 121, 5629-5634. 15. Zhang, X.; Meng, J.; Zhu, B.; Yu, J.; Zou, S.; Zhang, Z.; Gao, Y.; Wang, Y. In Situ TEM Studies of Shape Evolution of Pd Nanocrystals under Oxygen and Hydrogen Environment at Atmospheric Pressure. Chem. Commun. 2017, 53, 13213-13216. 20 ACS Paragon Plus Environment

Page 20 of 24

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

The Journal of Physical Chemistry

16. Jiang, Y.; Li, H.; Wu, Z.; Ye, W.; Zhang, H.; Wang, Y.; Sun, C.; Zhang, Z. In Situ Observation

of

Hydrogen-Induced

Surface

Faceting

for

Palladium–Copper

Nanocrystals at Atmospheric Pressure. Angew. Chem. Int. Ed. 2016, 55, 12427-12430. 17. Jiang, Y.; Zhang, Z.; Yuan, W.; Zhang, X.; Wang, Y.; Zhang, Z. Recent Advances in Gas-Involved in Situ Studies via Transmission Electron Microscopy. Nano Res. 2018, 11, 42-67. 18. Thostrup, P.; Christoffersen, E.; Lorensen, H. T.; Jacobsen, K. W.; Besenbacher, F.; Norskov, J. K. Adsorption-Induced Step Formation. Phys. Rev. Lett. 2002, 88, 126102. 19. Nguyen, L.; Liu, L.; Assefa, S.; Wolverton, C.; Schneider, W. F.; Tao, F. F. AtomicScale Structural Evolution of Rh(110) During Catalysis. ACS Catal. 2017, 7, 664-674. 20. Tao, F. F.; Nguyen, L.; Zhang, S.; Li, Y.; Tang, Y.; Zhang, L.; Frenkel, A. I.; Xia, Y.; Salmeron, M. Formation of Second-Generation Nanoclusters on Metal Nanoparticles Driven by Reactant Gases. Nano Lett. 2016, 16, 5001-5009. 21. Eren, B.; Zherebetskyy, D.; Patera, L. L.; Wu, C. H.; Bluhm, H.; Africh, C.; Wang, L.-W.; Somorjai, G. A.; Salmeron, M. Activation of Cu(111) Surface by Decomposition into Nanoclusters Driven by CO Adsorption. Science 2016, 351, 475-478. 22. Avanesian, T.; Dai, S.; Kale, M. J.; Graham, G. W.; Pan, X.; Christopher, P. Quantitative and Atomic-Scale View of CO-Induced Pt Nanoparticle Surface Reconstruction at Saturation Coverage via DFT Calculations Coupled with in Situ TEM and IR. J. Am. Chem. Soc. 2017, 139, 4551-4558. 23. Newton, M. A.; Belver-Coldeira, C.; Martinez-Arias, A.; Fernandez-Garcia, M. Dynamic in Situ Observation of Rapid Size and Shape Change of Supported Pd 21 ACS Paragon Plus Environment

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

Nanoparticles During CO/NO Cycling. Nat. Mater. 2007, 6, 528-532. 24. Tang, M.; Li, H.; Yuan, W.; Zou, S.; Sun, C.; Wang, Y. First-Principles Study of the Interactions of Hydrogen with Low-Index Surfaces of PdCu Ordered Alloy. Prog. Nat. Sci.: Mater. Int. 2017, 27, 709-713. 25. Ueda, K.; Suzuki, K.; Toyoshima, R.; Monya, Y.; Yoshida, M.; Isegawa, K.; Amemiya, K.; Mase, K.; Mun, B. S.; Arman, M. A., et al. Adsorption and Reaction of CO and NO on Ir(111) under near Ambient Pressure Conditions. Top. Catal. 2016, 59, 487-496. 26. Sreekanth, P. M.; Smirniotis, P. G. Selective Reduction of NO with CO over Titania Supported Transition Metal Oxide Catalysts. Catal. Lett. 2008, 122, 37-42. 27. Roobol, S. B.; Onderwaater, W. G.; van Spronsen, M. A.; Carla, F.; Balmes, O.; Navarro, V.; Vendelbo, S.; Kooyman, P. J.; Elkjaer, C. F.; Helveg, S., et al. In Situ Studies of NO Reduction by H2 over Pt Using Surface X-Ray Diffraction and Transmission Electron Microscopy. Phys. Chem. Chem. Phys. 2017, 19, 8485-8495. 28. Rainer, D. R.; Vesecky, S. M.; Koranne, M.; Oh, W. S.; Goodman, D. W. The CO+NO Reaction over Pd: A Combined Study Using Single-Crystal, Planar-ModelSupported, and High-Surface-Area Pd/Al2O3 Catalysts. J. Catal. 1997, 167, 234-241. 29. Hu, Y. H.; Griffiths, K.; Norton, P. R. Surface Science Studies of Selective Catalytic Reduction of NO: Progress in the Last Ten Years. Surf. Sci. 2009, 603, 1740-1750. 30. Tao, F.; Grass, M. E.; Zhang, Y.; Butcher, D. R.; Renzas, J. R.; Liu, Z.; Chung, J. Y.; Mun, B. S.; Salmeron, M.; Somorjai, G. A. Reaction-Driven Restructuring of RhPd and Pt-Pd Core-Shell Nanoparticles. Science 2008, 322, 932-934. 22 ACS Paragon Plus Environment

Page 22 of 24

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

The Journal of Physical Chemistry

31. Goldsmith, B. R.; Sanderson, E. D.; Ouyang, R.; Li, W.-X. CO-and NO-Induced Disintegration and Redispersion of Three-Way Catalysts Rhodium, Palladium, and Platinum: An Ab Initio Thermodynamics Study. J. Phys. Chem. C 2014, 118, 95889597. 32. Vesecky, S. M.; Chen, P.; Xu, X.; Goodman, D. W. Evidence for Structure Sensitivity in the High Pressure CO+NO Reaction over Pd(111) and Pd(100). J. Vac. Sci. Technol. A 1995, 13, 1539-1543. 33. Nakao, K.; Ito, S.-i.; Tomishige, K.; Kunimori, K. Infrared Chemiluminescence Study of CO2 Formation in CO + NO Reaction on Pd(110) and Pd(111) Surfaces. J. Phys. Chem. B 2005, 109, 17579-17586. 34. Wulff, G. Xxv. Zur Frage Der Geschwindigkeit Des Wachsthums Und Der Auflösung Der Krystallflächen. zkri 1901, 34, 449. 35. Fowler, R. H.; Guggenheim, E. A. Statistical Thermodynamics; Cambridge University Press: Cambridge, 1939, pp 431-450. 36. Jaroniec, M.; Borowko, M.; Patrykiejew, A. Consequences of Assuming the BraggWilliams Approximation in Mixed-Gas Adsorption. Surf. Sci. 1978, 78, L501-L503. 37. M.W. Chase, J., C.A. Davies, J.R. Downey, Jr. D.J. Frurip, R.A. McDonald, and A.N. Syverud Nist-Janaf Thermochemical Tables. http://kinetics.nist.gov/janaf/ (accessed 09,2016). 38. Ozensoy, E.; Goodman, D. W. Vibrational Spectroscopic Studies on CO Adsorption, NO Adsorption CO+ NO Reaction on Pd Model Catalysts. Phys. Chem. Chem. Phys. 2004, 6, 3765-3778. 23 ACS Paragon Plus Environment

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

TOC

24 ACS Paragon Plus Environment

Page 24 of 24