Experimental and Computational Interrogation of Fast SCR

Jun 21, 2017 - For a more comprehensive list of citations to this article, users are encouraged to ... Local dynamics of copper active sites in zeolit...
0 downloads 0 Views 2MB Size
Subscriber access provided by Binghamton University | Libraries

Article

Experimental and Computational Interrogation of Fast SCR Mechanism and Active Sites on H-Form SSZ-13 Sichi Li, Yang Zheng, Feng Gao, János Szanyi, and William F. Schneider ACS Catal., Just Accepted Manuscript • Publication Date (Web): 21 Jun 2017 Downloaded from http://pubs.acs.org on June 21, 2017

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

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

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

Experimental and Computational Interrogation of Fast SCR Mechanism and Active Sites on H-Form SSZ-13 Sichi Li,† Yang Zheng,‡ Feng Gao,‡ Janos Szanyi,‡ and William F. Schneider∗,† †Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, IN 46556, United States ‡Institute for Integrated Catalysis, Pacific Northwest National Laboratory, P.O. Box 999, Richland, WA 99352, United States E-mail: [email protected]

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

Abstract Experiments and density functional theory (DFT) models are combined to develop a unified, quantitative model of the mechanism and kinetics of fast selective catalytic reduction (SCR) of NO/NO2 mixtures over H-SSZ-13 zeolite. Rates, rate orders, and apparent activation energies collected under differential conditions reveal two distinct kinetic regimes. First-principles thermodynamics simulations are used to determine the relative coverages of free Brønsted sites, chemisorbed NH4+ and physisorbed NH3 as a function of reaction conditions. First-principles metadynamics calculations show that all three sites can contribute to the rate-limiting N−N bond forming step in fast SCR. The results are used to parameterize a kinetic model that encompasses the full range of reaction conditions and recovers observed rate orders and apparent activation energies. Observed kinetic regimes are related to changes in most-abundant surface intermediates.

Keywords: SCR catalysis, H-SSZ-13, density functional theory, kinetic modeling, ab inito molecular dynamics

2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

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

Introduction Nitrogen oxides NOx produced by industrial processes, during the combustion of fossil fuels, and by operation of light and heavy duty vehicles contribute adversely to air quality. The selective catalytic reduction (SCR) of NOx with ammonia is a widely used approach for abating these emissions. 1–4 Cation-exchanged zeolites have long been known to be active for NOx SCR and are the preferred catalysts for applications that demand high thermal stability and wide temperature windows of operation. 1,5,6 Exchanged redox-active metal ions, such as copper or iron, are usually required for a material to be active for “standard” SCR: 7,8

4 NO + O2 + 4NH3 −−→ 4N2 + 6H2 O

(1)

Catalytic performance for “fast” SCR, in which equimolar amounts of NO and NO2 are consumed: 2NO + 2 NO2 + 4 NH3 −−→ 4N2 + 6H2 O

(2)

is less sensitive to the nature of exchanged ions, and even Brønsted acid, or H-form, zeolites have sufficient SCR activity for some stationary source applications. 3,9–11 Even metalexchanged zeolites typically contain significant densities of residual Brønsted sites due to incomplete metal exchange. 12 Thus, it is important to isolate the contributions of these Brønsted acid sites to fast SCR even in studies related to the performance of metal exchange sites. The precise role of Brønsted sites in fast SCR remains unclear. NH3 binds to Brønsted sites as NH4+ , and a variety of experiments and theoretical studies on metal oxides 13–15 and H-form zeolites 16–18 implicate NH4+ as a key fast SCR reaction intermediate. Br¨ uggemann et al. 19 used density functional theory (DFT) calculations and cluster models to explore the reaction of N2 O3 with a zeolitic NH4+ : ZNH4 + N2 O3 −→ ZH + H2 NNO + HONO 3

ACS Paragon Plus Environment

(3)

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 4 of 30

Here we adopt “Z” as a symbol for an Al-substituted tetrahedral site in a zeolite. In this reaction NH4+ is oxidized by N2 O3 postulated to form in the well-known equilibrium: 20,21 NO + NO2 − ↽− −⇀ − N2 O3

(4)

Following reaction 3, decomposition of H2 NNO and reaction between HONO and NH3 catalyzed by Brønsted sites take place, producing final products N2 and H2 O. Br¨ uggemann et al. subsequently proposed NO adsorbed at Brønsted sites as an alternative fast SCR intermediate that is quickly consumed by NH3 to recover active Brønsted sites: 22 ZH + N2 O3 −−→ ZNO + HONO

(5)

ZNO + NH3 −−→ ZH + H2 NNO

(6)

Computed activation energies of reaction 3, reactions 5 and 6 are ≈ 100, 22 and 12 kJ mol−1 , respectively. Microkinetic modeling based on these reaction pathways can reproduce NOx conversion at high temperature reasonably close to experimental results but significantly underestimate low temperature ( 500 ppm at 633 K, while at 663 K total rate is constant up to 800 ppm NH3 .

Figure 11: Computed natural logarithm of total reaction rate versus logarithm of PNH3 at (a) 633 K and (b) 663 K. P ◦ = 1 atm. This behavior reflects the competing contributions of the three underlying routes to the total rate. The NH4+ route is least sensitive to PNH3 . It decreases slightly with PNH3 (< 0.5 in ln(rate) between 100 and 800 ppm NH3 ). The rate of the ZH route (reaction 5) also decreases 23

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

with PNH3 but with a much larger relative amplitude (≈ 2.5 in ln(rate) in same PNH3 range). The physisorbed NH3 route, conversely, increases with PNH3 (≈ 2.0 in ln(rate) in same PNH3 range). These rate-PNH3 trends of different routes can be rationalized by different dependence of site coverages on PNH3 as each individual rate is proportional to the corresponding site coverage. Within the simulated PNH3 range, coverage of ZNH4 (θZNH4 ) stays high and decreases slightly with PNH3 . In contrast, coverages of both ZNH4 · NH3 and ZH (θZNH4 · NH3 and θZH ) are relatively low and very sensitive to PNH3 that increase significantly with PNH3 . Under low PNH3 , NH4+ and ZH routes contribute more than physisorbed NH3 route, so the total rate follows the decreasing trend of NH4+ and ZH routes. As PNH3 increases further, rate of physisorbed NH3 route becomes comparable to the other two. When increase of rate of physisorbed NH3 route balances the decrease of the other two, total rate reaches the flat region. When rate of physisorbed NH3 route surpasses the other two, the total rate begins to increase with PNH3 . These simulated results are qualitatively consistent with experimental measurements, rationalizing the experimentally observed negative-zero(-positive) transition of NH3 rate order (Fig. 2c).

Conclusions The role of zeolitic Brønsted acid sites in fast SCR has been debated for over a decade. What remains unclear is the mechanistic origin of high fast SCR activity observed at low temperature. We show here through theory that governed by NH3 adsorption thermodynamics, the state of Brønsted sites in H-SSZ-13 changes dynamically as a function of temperature, from NH3 -solvated NH4+ to free NH4+ under reaction conditions. Abundance of physisorbed NH3 hydrogen bonded to NH4+ ions and low energy barrier of reaction between N2 O3 and such physisorbed NH3 species rationalize low temperature catalytic performance. Further, we show through experiment that fast SCR exhibits two distinct kinetics regimes in the temperature window between 473 and 673 K that is rationalized by change of active sites for

24

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30

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

N2 O3 consumption revealed by first-principle kinetic modeling. Since this mechanism has no dependence on zeolitic topology, we expect the same kinetics can be observed in H-form zeolites with different framework topologies, such as H-ZSM-5. These findings demonstrate the importance of isolating reaction conditions in comprehensive investigation of kinetics of a complex reaction network. The DFT approach based on non-equilibrium AIMD and first-principles thermodynamics is particularly powerful to study mechanisms of catalytic systems, especially for those containing multiple competing reactive intermediates, adsorbates and/or sites, such as SCR on metal exchanged zeolites.

Supporting Information Available XRD pattern of the catalyst, NOx conversion plots under different space velocities of gas feed, rate orders of NOx , computational details are included in the associated PDF file.

Acknowledgement Financial support was provided by the National Science Foundation GAOLI program under award number 1258690-CBET. We thank the Center for Research Computing at Notre Dame, and EMSL, a DOE Office of Science User Facility sponsored by the Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory, for support of computational resources. Authors from PNNL (YZ, FG and JZ) gratefully acknowledge supports from the United States Department of Energy (DOE), Energy Efficiency and Renewable Energy, Vehicle Technologies Office.

References (1) Paolucci, C.; Iorio, J. D.; Ribeiro, F.; Gounder, R.; Schneider, W. Adv. Catal. 2016, 59, 1–107. 25

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

(2) Zhang, R.; Liu, N.; Lei, Z.; Chen, B. Chem. Rev. 2016, 116, 3658–3721. (3) Beale, A. M.; Gao, F.; Lezcano-Gonzalez, I.; Peden, C. H. F.; Szanyi, J. Chem. Soc. Rev. 2015, 44, 7371–7405. (4) Kaiser, R.; Rusch, K. MTZ worldwide 2007, 68, 17–20. (5) Chen, H.-Y. In Urea-SCR Technology for deNOx After Treatment of Diesel Exhausts; Nova, I., Tronconi, E., Eds.; Springer New York: New York, NY, 2014; pp 123–147. (6) Centi, G.; Perathoner, S. Appl. Catal., A 1995, 132, 179 – 259. (7) Lezcano-Gonzalez, I.; Deka, U.; Arstad, B.; Van Yperen-De Deyne, A.; Hemelsoet, K.; Waroquier, M.; Van Speybroeck, V.; Weckhuysen, B. M.; Beale, A. M. Phys. Chem. Chem. Phys. 2014, 16, 1639–1650. (8) Gao, F.; Kwak, J. H.; Szanyi, J.; Peden, C. H. F. Toc. Catal. 2013, 56, 1441–1459. (9) Devadas, M.; Kr¨ ocher, O.; Elsener, M.; Wokaun, A.; S¨ oger, N.; Pfeifer, M.; Demel, Y.; Mussmann, L. Appl. Catal., B 2006, 67, 187–196. (10) Schwidder, M.; Heikens, S.; Toni, A. D.; Geisler, S.; Berndt, M.; Br¨ uckner, A.; Gr¨ unert, W. J. Catal. 2008, 259, 96–103. (11) Pence, D.; Thomas, T. Reduction of nitrogen oxides with catalytic acid resistant aluminosilicate molecular sieves and ammonia. 1980; US Patent 4,220,632. (12) Brandenberger, S.; Kr¨ ocher, O.; Tissler, A.; Althoff, R. Catal. Rev. 2008, 50, 492–531. (13) Takagi, M.; Kawai, T.; Soma, M.; Onishi, T.; Tamaru, K. J. Catal. 1977, 50, 441–446. (14) Takagi, M.; Kawai, T.; Soma, M.; Onishi, T.; Tamaru, K. J. Phys. Chem. 1976, 80, 430–431.

26

ACS Paragon Plus Environment

Page 26 of 30

Page 27 of 30

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

(15) Arnarson, L.; Falsig, H.; Rasmussen, S. B.; Lauritsen, J. V.; Moses, P. G. J. Catal. 2017, 346, 188–197. (16) Long, R.; Yang, R. J. Catal. 2002, 207, 224–231. (17) Sanchez-Escribano, V.; Montanari, T.; Busca, G. Appl. Catal., B 2005, 58, 19–23. (18) Wallin, M.; Karlsson, C.-J.; Palmqvist, A.; Skoglundh, M. Top. Catal. 2004, 30, 107– 113. (19) Br¨ uggemann, T. C.; Przybylski, M.-D.; Balaji, S. P.; Keil, F. J. J. Phys. Chem. C 2010, 114, 6567–6587. (20) Smith, I. W.; Yarwood, G. Chem. Phys. Lett. 1986, 130, 24–28. (21) Glendening, E. D.; Halpern, A. M. J. Chem. Phys. 2007, 127, 164307. (22) Br¨ uuggemann, T. C.; Vlachos, D. G.; Keil, F. J. J. Catal. 2011, 283, 178–191. (23) Li, J.; Li, S. Phys. Chem. Chem. Phys. 2007, 9, 3304–3311. (24) Hunger, B.; Hoffmann, J. Thermochim. Acta 1986, 106, 133–140. (25) Parrillo, D. J.; Gorte, R. J.; Farneth, W. E. J. Am. Chem. Soc. 1993, 115, 12441–12445. (26) Dragoi, B.; Gervasini, A.; Dumitriu, E.; Auroux, A. Thermochim. Acta 2004, 420, 127–134. (27) Armandi, M.; Bonelli, B.; Bottero, I.; Are´an, C. O.; Garrone, E. J. Phys. Chem. C 2010, 114, 6658–6662. (28) Jones, A. J.; Iglesia, E. ACS Catal. 2015, 5, 5741–5755. (29) Paolucci, C.; Parekh, A. A.; Khurana, I.; Iorio, J. R. D.; Li, H.; Caballero, J. D. A.; Shih, A. J.; Anggara, T.; Delgass, W. N.; Miller, J. T.; Ribeiro, F. H.; Gounder, R.; Schneider, W. F. J. Am. Chem. Soc. 2016, 138, 6028–6048. 27

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

(30) Giordanino, F.; Borfecchia, E.; Lomachenko, K. A.; Lazzarini, A.; Agostini, G.; Gallo, E.; Soldatov, A. V.; Beato, P.; Bordiga, S.; Lamberti, C. J. Phys. Chem. Lett. 2014, 5, 1552–1559. (31) L´ onyi, F.; Valyon, J. Microporous Mesoporous Mater. 2001, 47, 293–301. (32) Fickel, D. W.; Lobo, R. F. J. Phys. Chem. C 2010, 114, 1633–1640. (33) Fickel, D. W.; D’Addio, E.; Lauterbach, J. A.; Lobo, R. F. Appl. Catal., B 2011, 102, 441–448. (34) Perdew, J. P.; Wang, Y. Phys. Rev. B 1992, 45, 13244–13249. (35) CPMD, http://www.cpmd.org/, Copyright IBM Corp 1990-2016, Copyright MPI f¨ urFestk¨ orperforschung Stuttgart 1997-2001. (36) Vanderbilt, D. Phys. Rev. B 1990, 41, 7892–7895. (37) Laasonen, K.; Car, R.; Lee, C.; Vanderbilt, D. Phys. Rev. B 1991, 43, 6796–6799. (38) Laasonen, K.; Pasquarello, A.; Car, R.; Lee, C.; Vanderbilt, D. Phys. Rev. B 1993, 47, 10142–10153. (39) Wong, K.-Y.; York, D. M. J. Chem. Theory Comput. 2012, 8, 3998–4003. (40) den Otter, W. K. J. Chem. Theory Comput. 2013, 9, 3861–3865. (41) Grimme, S. J. Comput. Chem. 2006, 27, 1787–1799. (42) Iannuzzi, M.; Laio, A.; Parrinello, M. Phys. Rev. Lett. 2003, 90, 238302. (43) Laio, A.; Rodriguez-Fortea, A.; Gervasio, F. L.; Ceccarelli, M.; Parrinello, M. J. Phys. Chem. B 2005, 109, 6714–6721. (44) Car, R.; Parrinello, M. Phys. Rev. Lett. 1985, 55, 2471–2474.

28

ACS Paragon Plus Environment

Page 28 of 30

Page 29 of 30

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

(45) Grossale, A.; Nova, I.; Tronconi, E. Catal. Lett. 2009, 130, 525–531. (46) Mihai, O.; Tamm, S.; Stenfeldt, M.; Wang-Hansen, C.; Olsson, L. Ind. Eng. Chem. Res. 2015, 54, 11779–11791. (47) Tronconi, E.; Nova, I.; Ciardelli, C.; Chatterjee, D.; Weibel, M. J Catal. 2007, 245, 1–10. (48) Nova, I.; Ciardelli, C.; Tronconi, E.; Chatterjee, D.; Weibel, M. AIChE J. 2009, 55, 1514–1529. (49) Forzatti, P.; Nova, I.; Tronconi, E.; Kustov, A.; Thøgersen, J. R. Catal. Today 2012, 184, 153–159. (50) Grossale, A.; Nova, I.; Tronconi, E.; Chatterjee, D.; Weibel, M. J. Catal. 2008, 256, 312–322. (51) Nova, I.; Ciardelli, C.; Tronconi, E.; Chatterjee, D.; Weibel, M. Top. Catal. 2007, 42, 43–46. (52) Nova, I.; Ciardelli, C.; Tronconi, E.; Chatterjee, D.; Bandl-Konrad, B. Catal. Today 2006, 114, 3–12. (53) Campbell, C. T.; Sellers, J. R. V. J. Am. Chem. Soc. 2012, 134, 18109–18115. (54) Piccini, G.; Sauer, J. J. Chem. Theory Comput. 2014, 10, 2479–2487. (55) Piccini, G.; Alessio, M.; Sauer, J.; Zhi, Y.; Liu, Y.; Kolvenbach, R.; Jentys, A.; Lercher, J. A. J. Phys. Chem. C 2015, 119, 6128–6137. (56) Paolucci, C.; Verma, A. A.; Bates, S. A.; Kispersky, V. F.; Miller, J. T.; Gounder, R.; Delgass, W. N.; Ribeiro, F. H.; Schneider, W. F. Angew. Chem., Int. Ed. 2014, 53, 11828–11833. (57) Richter, M.; Eckelt, R.; Parlitz, B.; Fricke, R. Appl. Catal., B 1998, 15, 129–146. 29

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

(58) Sun, D.; Schneider, W. F.; Adams, J. B.; Sengupta, D. J. Phys. Chem. A 2004, 108, 9365–9374. (59) Tiwary, A. S.; Mukherjee, A. K. J. Mol. Struct.: THEOCHEM 2009, 909, 57–65. (60) Chen, H.-Y.; Sun, Q.; Wen, B.; Yeom, Y.-H.; Weitz, E.; Sachtler, W. M. Catal. Today 2004, 96, 1–10. (61) Dusenbury, J. H.; Powell, R. E. J. Am. Chem. Soc. 1951, 73, 3266–3268. (62) Veley, V. H. J. Chem. Soc., Trans. 1903, 83, 736–749. (63) Lu, X.; Musin, R. N.; Lin, M. C. J. Phys. Chem. A 2000, 104, 5141–5148. (64) Chase, M.W., Jr., NIST-JANAF Themochemical Tables, Fourth Edition, J. Phys. Chem. Ref. Data, Monograph 9, 1998, 1-1951.

TOC

30

ACS Paragon Plus Environment

Page 30 of 30