An Aging Model of NH 3 Storage Sites for Predicting Kinetics of NH 3 Adsorption, Desorption and Oxidation over Hydrothermally Aged Cu-Chabazite

: A uniﬁed transient kinetic model which can predict the adsorption, desorption and oxidation kinetics of NH 3 over hydrothermally aged Cu-chabazite was developed. The model takes into account the variation of fractional coverages of NH 3 storage sites due to hydrothermal aging. In order to determine the fractional coverage of these sites, the catalyst was aged for various times at a certain temperature followed by NH 3 adsorption, desorption and temperature-programmed desorption (TPD) experiments. TPD proﬁles were deconvoluted mainly into three peaks with centres at 317, 456 and 526 ◦ C, respectively. Hydrothermal aging resulted in the progressive increase in the intensity of the peak at 317 ◦ C and decrease in the intensity of the peaks at 456 and 526 ◦ C, along with decreased NH 3 oxidation at high temperatures. A model for hydrothermal aging kinetics of the fractional coverage of storage sites was developed using three reactions with appropriate rate expressions with parameters regressed from experimental data. The model was then incorporated into a multi-site kinetic model for the degreened Cu-Chabazite by the addition of aging reactions on each storage site. The e ﬀ ects of both aging time and temperature on the kinetics NH 3 adsorption, desorption and oxidation were successfully predicted in the 155-540 ◦ C range. This study is the ﬁrst step towards the development of a hydrothermal aging-uniﬁed kinetic model of NH 3 -Selective Catalytic Reduction over Cu-chabazite.


Introduction
An engine aftertreatment system (ATS) of a heavy-duty vehicle comprises an NH 3 -Selective Catalytic Reduction (NH 3 -SCR) reactor in which NOx and NH 3 react to form N 2 and H 2 O via a variety of reactions including Standard, Fast and NO 2 -SCR. Cu-Chabazite (CHA) is the catalyst of choice due to its good deNOx performance and better high temperature stability as compared to other Cu-zeolites [1]. Much-speculated upcoming Euro7 emission regulations will impose more stringent NOx emission standards than before and also restrict CO 2 emissions. The latter would indeed require increasing the combustion efficiency of the engine by lowering engine outlet temperatures, which would lower SCR efficiencies. Therefore, intense engineering efforts are directed towards the design, calibration and control of the SCR reactors in the future ATSs of heavy-duty vehicles with lean burn diesel engines. The design and calibration of SCR reactors is usually carried out using data from engine dynamometer and vehicle tests that scan a large operation region of the engine. However, such tests are very long, between the ammonia storage capacity loss and the corresponding changes in the pseudo-first-order rate constants of the two reactions most affected by aging, i.e., ammonia oxidation and standard SCR, during the hydrothermal aging process for Fe-zeolite [49]. The developed aging model was able to describe effectively all their experimental data, collected at different aging times and aging temperatures, on a predictive basis. The study concluded that such an approach may also be applicable to other types of catalyst formulations. A similar treatment cannot be easily carried out for Cu-CHA since the total NH 3 storage does not change up to very long aging times. There have also been studies on the effects of aging on the nature and concentration of active sites on Cu-CHA which do not involve kinetic models. It has been found that the resistance of Cu-CHA to HA is due to the exceptional hydrothermal stability of the Z 2 Cu II sites located on 6MR, which also act as structure stabilizers during HA [50][51][52] and prevent dealumination. HA mostly affects the ZCu II OH sites which in turn convert to CuO x clusters [9,50,53] and cause some dealumination [28,54], and clusters in turn grow and destroy the zeolite cage [55] resulting in activity loss with hard HA. For mild aging, these changes in catalyst performance were attributed to the changes in the concentration of the individual active Cu sites. In a pioneering study on aging, Luo et al. proposed that ZCu II OH reacted with ZH to produce Z 2 Cu II and water upon HA, based on the behaviour of NH 3 -TPD profiles. A consistent increase in the intensity of the lower temperature peak around 300 • C and a decrease in the intensity of the higher temperature peaks (>400 • C) were observed. The decrease in the intensity of higher temperature peaks was attributed to the decrease in the Brönsted sites, whereas the increase of the 300 • C peak was attributed to the formation of new Z 2 Cu II sites which were assumed to accommodate more NH 3 per Cu than ZCu II OH [56,57]. The changes in the active site concentrations upon HA were also supported by H 2 -TPR and DRIFTS studies [56]. Recently, the same group were able to link the HA kinetics of the active sites to a transient kinetic model of NH 3 adsorption and desorption of fresh Cu-CHA in order to simulate the NH 3 -TPD profiles of the aged catalysts. This was based on the mean-field approximation, implying constant turnover Arrhenius rate constants for the adsorption and desorption reactions, but with changing individual site densities directly affecting the rates associated with the aforementioned reactions [58]. However, it has also been shown that via H 2 -TPR, NH 3 -TPD, DRIFTS [29,30,50] and in situ time-resolved XAS [26] that the NH 3 -TPD profiles of Cu-CHA are much more complicated than proposed in [56]. Additionally, the multi-site kinetic models to predict deNOx performance [30,44,59] suggest that the highest temperature peak of NH 3 -TPD profile of Cu-CHA is associated with NH 3 desorption from ZCu II OH sites residing at the 8MR.
Despite the numerous spectroscopic and kinetic studies on the issue, it appears that there is no unambiguous agreement as to which species are responsible for which TPD peak and there is a need for more detailed analysis and kinetic modelling of the NH 3 -TPD profiles of degreened Cu-CHA. Moreover, there is no report in the literature on the kinetic modelling of HA effects coupled with NH 3 storage and oxidation models of the degreened catalysts to predict the performance over hydrothermally aged Cu-CHA. Here, a HA-unified transient kinetic model which can predict the NH 3 adsorption breakthrough, isothermal desorption of NH 3 and NH 3 -TPD behaviors upon HA between 650-800 • C range up to 46h, is presented for Cu-CHA. Using this model, the NH 3 oxidation performance in the 155-540 • C range, from its freshly degreened state up to the mild HA state (in the 650-800 • C range, up to 46h of isothermal HA), is predicted for the first time in the literature, based on the kinetic parameters of the degreened catalyst. For this purpose, an HA kinetics model consisting of four NH 3 storage sites, was developed by determining the changes in fractional coverage of these sites by deconvoluting experimental NH 3 -TPD data only. Then, the active site HA kinetics model was coupled with the transient kinetic model for the degreened Cu-CHA while keeping the turnover rate constants and activation energies associated with adsorption, desorption and oxidation reactions unchanged upon HA. Figure 1 shows a typical profile for the adsorption, isothermal desorption and TPD of NH 3 over degreened Cu-CHA. The TPD profile was fitted with four Gaussian peaks associated with different NH 3 storage sites, hereafter referred to as ZCu1, ZCu2 and ZB, with peak centers at 317, 456 and 526 • C, respectively, along with a small contribution from Z 2 W with a peak center at 126 • C. Z 2 W site is associated with sites on which either physically adsorbed NH 3 [30] or Cu(NH 3 ) 4 species with very low or no reactivity towards SCR related reactions reside [26]. Its concentration was obtained by adding the amount of NH 3 desorbed during the isothermal desorption period and the amount desorbed corresponding to the small peak with the center at 126 • C. ZCu1 and ZCu2 are associated with Cu sites responsible for SCR-related reactions commonly referred to as Z 2 Cu II and ZCu II OH, which reside at the 6 and 8 MRs, respectively [23]. ZB sites are associated with Brönsted sites that were not exchanged with Cu during catalyst synthesis.  Figure 2 shows the NH 3 -TPD profiles for Cu-CHA with increasing aging times at 650, 700 and 800 • C. The intensity of the peaks in the higher temperature region of the curves decreased and the intensity of the peaks in the lower temperature region increased with increasing time, in line with the observations in the literature [57]. Meanwhile, the total NH 3 storage stayed nearly constant with increasing aging time at all the temperatures investigated in this study ( Figure 2d). Figure 3 illustrates the peaks obtained from the deconvolution of the data in Figure 2. The peaks are associated with Z 2 W, ZCu1, ZCu2 and ZB sites. Data clearly show the rapid decrease in the intensity of the peak associated with ZB, a much slower decrease in the intensity of the peak associated with ZCu2 and an increase in the intensity of the peak associated with ZCu1 with increasing times of HA. Figure 4 shows the changes in the fractional coverages of sites with time ( Figure 4a,c,e) along with the sequential HA and NH 3 adsorption/TPD sets. HA sets given in Figure 4b,d,f are for durations of 109, 60 and 50 h (i.e., totals of 46, 25 and 21 h isothermal HA at 650, 700 and 800 • C, respectively, heating and cooling ramps, and also includes the time required for NH 3 adsorption and TPD experiments). After HA at 650 • C (Figure 4a), the fractional coverage of ZCu1 slightly decreased and then increased with time to reach a value 28% higher than the degreened one, whereas the coverage of ZCu2 slightly increased by 12%, and then decreased with time to reach a value 22% lower than the degreened one. Initial changes in fractional coverages upon HA indicates that the initial degreening treatment (550 • C, 2 h) was not enough to produce the maximum possible concentration of ZCu2. Cu mobility during HA is a very complex process where both Cu sites can interchange among themselves up to a saturation in Cu2 concentration. The rate of initial increase in fractional coverage of ZCu2 with aging time was higher than its rate of decrease due to HA. Meanwhile, the fractional coverage of ZB decreased much faster than the fractional coverage of ZCu2, indicating that the decay mechanism for these two species could be different during HA. The fractional coverage of Z 2 W initially increased slightly and then stayed constant throughout. For runs at 700 • C ( Figure 4c) and 800 • C (Figure 4e), the initial decrease and increase in the fractional coverages of ZCu1 and ZCu2 were not observed, most probably because this process is happening much faster than HA kinetics. Similar to the data at 650 • C, fractional coverage of ZCu1 increased with time whereas fractional coverage of ZCu2 and ZB decreased. Again, the decay rates of ZB were much faster than those of ZCu2. Based on the above observations, three reactions for describing the changes in fractional coverages with aging time were proposed. The reactions for aging and their rate expressions are shown in Table 1 (Reactions 1-3). The main reaction responsible for HA is Reaction 3, where ZCu2 sites convert to ZCu1. The reaction responsible for the decrease in coverage of ZB, and initial decrease and increase of coverages of ZCu1 and ZCu2, respectively, was Reaction 2. Since two empty sites (Z) were used as reactants in Reaction 2, an empty Z site was formed. It was assumed that this site was converted to Z 2 W via reacting with ZCu1, as given by Reaction 1, to describe the slight increase in coverage of Z 2 W upon HA. Equation (12) (given in Section 3.3.1) was written for each surface species (Z 2 W, ZCu1, ZCu2, ZB) using the rate expressions for Reactions 1-3. Rate equation parameters were then regressed from fractional coverage versus time data given in Figure 4 and are presented in Table 2. Model results are in excellent agreement with the measured values at all times and at each aging temperature. Thus, all of the changes in fractional coverages of Z 2 W, ZCu1, ZCu2 and ZB with aging time were predicted very well at each aging temperature. The HA reaction mechanism in [56] was also investigated in this study using elementary rate laws. However, the mechanism in [56] was not able to describe accurately the HA data obtained in this study. Particularly, that model could not describe the fast decrease in fractional coverage of ZB with aging time and rather slow decrease in fractional coverage of ZCu2 at longer times.

Reaction Number Reaction Rate Expression
Reactions related to HA Reactions related to NH 3 adsorption and desorption Reactions related to NH 3 oxidation

4-Site Kinetic Model for Degreened Cu-CHA
In order to model the adsorption, desorption and oxidation of NH 3 over Cu-CHA from a degreened state up to a hydrothermally aged state, firstly, a transient reaction kinetics model for the degreened Cu-CHA was developed. For the adsorption and desorption of NH 3 , the reactions and the associated rate expressions for each reaction are given in Table 1 (Reactions 4-11). The site densities of Z 2 W, ZCu1, ZCu2 and ZB of the degreened Cu-CHA were 64.0, 45.1, 24.6 and 14.6 mol site /m 3 , respectively, which were obtained from NH 3 adsorption/TPD profile of the degreened catalyst. The site Z 2 W occcupied two Z sites which resulted in a total Z concentration of 212.3 mol site /m 3 . The fractional coverages of Z 2 W, ZCu1, ZCu2 and ZB for the degreened sample were, therefore, 0.60, 0.21, 0.12 and 0.07, respectively. The stoichiometry between a certain active site and NH 3 [18,47,56] was lumped into to the reaction rate parameters, similar to the prior literature [30,44]. Figure 5 shows the NH 3 concentration at the outlet of the reactor during NH 3 adsorption, isothermal desorption and TPD along with the fractional coverages obtained using the 4-site transient kinetic model for the degreened Cu-CHA. The parameters for the rate expressions of Reactions 4-11 were obtained by fitting the experimental data shown in Figure 5a to the 4-site transient kinetic model (given in Section 3.3.1) and are presented in Table 2. The model captured very well all of the vital features of the experimental data, which are the adsorption breakthrough, isothermal desorption and the TPD profiles. The peaks at 317 and 456 • C in the TPD profile, along with the shoulder at 526 • C, were also very well predicted by the model. Here, in order for our approach to predict NH 3 storage and oxidation after HA, the kinetic model for the adsorption and desorption should be consistent with the positions of the fitted TPD peaks for the degreened catalyst ( Figure 1). To elaborate, the centers of the fitted peaks (in terms of temperature) should match the temperatures where the coverages of the associated individual sites in the kinetic model are halved during TPD. Figure 5b shows the changes in the coverages of ZCu1NH 3 , ZCu2NH 3 , and ZBNH 3 species during TPD of NH 3 of the degreened Cu-CHA. During TPD, as the temperature was raised from 110 to 600 • C, NH 3 started to desorb from the aforementioned sites according to Reactions 5,7,9 and 11. The predicted temperatures where the coverages were halved were at 304, 451 and 513 • C for ZCu1NH 3 , ZCu2NH 3 and ZBNH 3 (Figure 5b), respectively, which very closely match the fitted values of 317, 456 and 526 • C (Figures 1 and 2a). This shows that the kinetic model for the adsorption and desorption of NH 3 is consistent with the deconvoluted peaks of the experimental TPD curve of the degreened Cu-CHA. For the development of a model for NH 3 oxidation on degreened Cu-CHA, two reactions (Reactions 12-13) were considered. The oxidation reactions were assumed to occur on ZCu1 and ZCu2 sites. The parameters for rates of NH 3 oxidation reactions were obtained by fitting transient experimental NH 3 oxidation data of the degreened Cu-CHA to the model while keeping the parameters of the rates of Reactions 4-11 constant. Parameters are presented in Table 2. The model described the NH 3 oxidation performance of the degreened catalyst very well.

Modeling the Effects of HA on NH 3 Adsorption, Desorption and Oxidation
HA resulted in various changes in the fractional coverages of the NH 3 storage sites over Cu-CHA, as previously demonstrated. These changes needed to be coupled with the 4-site transient kinetic model in order to predict the NH 3 adsorption, desorption and oxidation performance of the catalyst upon HA. Thus, the effects of HA on the adsorption and desorption was modeled using Reactions 1-11 and the effects of HA on the oxidation behavior was modeled using Reactions 1-13. The rate constants for Reactions 4-13 obtained using data for the degreened Cu-CHA were used during the simulations.
For adsorption and desorption of NH 3 , Figure 6a,c,e illustrates the experimental and predicted NH 3 concentration at the outlet of the reactor, together with temperature profiles (HA treatment at 650 • C (Figure 6a), 700 • C (Figure 6c) and 800 • C (Figure 6e)). Generally, there is a very good agreement between experimental data and model predictions upon sequential HA. Figure 6b shows a comparison of experimental and predicted NH 3 concentration after 46 h of isothermal HA at 650 • C (last NH 3 adsorption/desorption experiment in Figure 6a). Model predictions agreed very well with experimental data. The increase in the intensity of the peak centered at 317 • C and the decreases in the intensity of the peaks centered at 456 and 526 • C upon HA were captured very well by the model. Figure 6d,f show a comparison of experimental and predicted NH 3 concentrations after 25 h of HA at 700 • C and 21 h of HA at 800 • C, respectively. Here, the model again described the important features of the experiment very well, including the NH 3 breakthrough, isothermal desorption and bi-modal TPD behavior of NH 3 .  Figure 7 shows a comparison of the measured and predicted NH 3 conversion with temperature for the degreened and aged samples for the oxidation of NH 3 . Three NH 3 oxidation experiments were performed subsequent to the completion of the three HA sets carried out at 650, 700 and 800 • C, shown in Figure 6a,c,e (the conditions of aging also shown in Figure 4b,d,f). Figure 7a shows the measured and simulated data for NH 3 oxidation experiment after 46 h of isothermal HA at 650 • C.
Here, the experimental data show decreased NH 3 oxidation activity at temperatures above 330 • C, as compared to the degreened Cu-CHA. The model is in excellent agreement with the measured values. Thus, the phenomena associated with HA causing a decrease in NH 3 oxidation activity could be modeled via just accounting for the changes in the fractional coverage of NH 3 storage sites of the Cu-CHA, while keeping the turnover rate constants and the activation energies of the Reactions 12-13 obtained from the NH 3 oxidation experiment of degreened Cu-CHA unchanged. The predicted changes in the fractional coverage of storage sites upon HA are affecting the values required for Equations (3) and (12), thereby enabling the prediction of NH 3 oxidation performance after HA. Similarly, good agreement between experimental data and model predictions was obtained for NH 3 oxidation experiments after 25 h of isothermal HA at 700 • C (Figure 7b). After 21 h of isothermal HA at 800 • C (Figure 7c), the model could capture the decrease in conversion of NH 3 upon HA as compared to the degreened sample. However, the model underestimated the conversion of NH 3 at temperatures higher than 350 • C. The fact that the model could capture the decreases in conversion upon HA at 650 and 700 • C has several important implications. Firstly, this shows that ZCu2 was mostly responsible for the high-temperature NH 3 oxidation activity, since the decay in its coverage upon HA was well correlated with both the behavior of the TPD profile and NH 3 oxidation. Secondly, it illustrates that the HA reactions and rate expressions proposed here are valid and the aforementioned active species survive the HA treatment without forming CuO x clusters, or if these CuO x clusters were formed they had the same turnover rate constants. The under-estimation of the NH 3 conversion after very a high temperature HA at 800 • C could therefore be attributed to the formation of the aforementioned CuO x clusters with different reactivites than ZCu2.

Laboratory Tests
The catalyst used in this study was a commercial Cu-chabazite washcoated to a cordierite monolith (300 cpsi -5 mils). Cylindrical cores with a length of 2.2 cm and a diameter of 1.8 cm were used in the runs. The experiments were carried out in a synthetic gas bench described in [60]. In a typical run, the monolith was first wrapped with a ceramic fibre paper and then loaded in a quartz reactor and placed in the isothermal zone of an electric tubular furnace (Thermo Scientific Lindberg Blue M). The reaction temperature was monitored using a J-type thermocouple placed 0.5 cm upstream of the catalyst. NH 3 , O 2 , CO 2 and N 2 were delivered using calibrated mass flow controllers (Brooks Instruments) whereas H 2 O was delivered using a peristaltic pump (Gilson Minipuls 3). All lines before and after the reactor were heated to 190 • C. The concentrations of the species at the outlet of the reactor were continuously monitored using a FTIR spectrometer (MKS Multigas 2030). Initially, the catalyst was degreened at 550 • C using a stream consisting of 5% H 2 O and 8% O 2 in N 2 for 2h. After all of the experiments, the catalyst was hydrothermally aged, as explained below. Runs were performed with a space velocity of 40,000 h −1 (STP). Gases used in this study were 5.0 grade or above. NH 3 adsorption/TPD experiment included three steps, which are the adsorption, isothermal desorption, and TPD. In these experiments, NH 3 was continuously fed to the reactor for 1h at 110 • C, which resulted in an adsorption breakthrough curve. Subsequently, NH 3 feed was stopped and physisorbed NH 3 desorbed during the isothermal desorption step which continued for 1h. This was followed by a temperature ramp (10 • C/min) from 110 to 600 • C. The feed during the adsorption step was 570 ppm NH 3 , 5% H 2 O, 10% CO 2 in an N 2 balance. During the isothermal desorption and TPD steps, no NH 3 was fed. Fractional coverages of the NH 3 storage sites were obtained using the NH 3 adsorption, desorption and TPD experiments. TPD profiles were deconvoluted into four Gaussian peaks using Fityk software (version 1.3.1). Initially, peaks were added manually according to the experimental data. The best fits were then obtained via minimizing the weighted sum of squared residuals between the experimental data and model using Levenberg-Marquardt algorithm. For NH 3 oxidation experiments, the catalyst was exposed to a stream containing 500 ppm NH 3 , 10% CO 2 , 8% O 2 , 5% H 2 O in N 2 , and temperature was increased from 155 to 540 • C in a stepwise manner. Downstream NH 3 concentration was measured. NH 3 conversion (%) at steady-state was given by [(F NH3,z0 -F NH3,zL )/F NH3,z0 ]x100. The NH 3 oxidation performance of the degreened catalyst was measured for a fresh core after the degreening.

HA Methodology
HA was performed on the test rig described above. HA procedure included three sets of aging at temperatures of 650, 700 and 800 • C. All sets started with a fresh and degreened Cu-CHA sample. Thus, three fresh Cu-CHA washcoated cores were used. Once the sample was degreened, a certain HA set was started with an NH 3 adsorption/TPD experiment, as previously described in Section 3.1, to measure fractional coverages, which was followed by heating to the desired temperature (either 650, 700 or 800 • C) for an isothermal HA step for a desired time. Subsequently, the temperature was lowered to 110 • C to perform NH 3 adsorption/TPD experiment, which was again followed by heating and an isothermal HA step at the desired temperature (either 650, 700 or 800 • C) of that HA set. The above-described sequential procedure continued for up to total isothermal HA times 46, 25 and 21 h, for the isothermal HA temperatures of 650, 700 and 800 • C, respectively. The isothermal HA times accounted for the times passed at the desired isothermal HA temperature. During the sets at 650, 700 and 800 • C, 9, 7 and 5 NH 3 adsorption/TPD experiments were performed, respectively, between the HA steps, and the fractional coverages of the NH 3 storage sites were obtained accordingly after each HA step with respect to time and temperature. Heating and cooling ramps along with the isothermal HA step were carried out with a feed of 5% H 2 O, 8% O 2 in N 2 balance with a space velocity of 40,000 h −1 (stp). NH 3 oxidation performance was measured as previously described in Section 3.1 after the completion of each set.

The Reactor Model
The kinetic modelling was performed using GT-POWER, version 2018 [61]. Material balances for the modelling of NH 3 adsorption, desorption and oxidation over Cu-CHA washcoated monolith are provided below. All of the models included necessary momentum and energy balances. The length of the monolith was discretized to ten parts with uniform channels. This was sufficient to obtain data that do not vary with increasing mesh size. The unsteady-state species term for the bulk gas phase was neglected due to a very small residence time as compared to the time steps of interest.
The material balance for the species in the gas phase is The material balance for the external mass transfer film is The material balance for the washcoat gives with the flux balance at the fluid-washcoat interface in the radial direction as boundary condition and with δ = f wc S GT-POWER (version 2018) addresses the washcoat pore diffusion in a manner similar to the published literature [62]. In the model, the following expressions were utilized to obtain the effective diffusivity, Knudsen diffusivity and the binary diffusion coefficients (Fuller equation [63]).
The fractional coverages are given by A k dθ k dt = j σ k j a j r j (12) During the modelling of HA kinetics, the changes in the fractional coverages of the four storage sites were obtained using Equation (12). This is based on the assumption that there is no transport of sites in and out of the reactor, and that the sites are homogeneously dispersed, are in close-proximity to react with each other, but immobile throughout the washcoat in all directions.

Kinetic Model
All reactions and associated rate equations are given in Table 1. Activation energy for the adsorption reactions (Reactions 4, 6, 8 and 10) was assumed to be zero. The temperature dependence of the turnover rate constant, k j was given by the Arrhenius equation where A j is the turnover rate constant and E A,j is the activation energy in reaction j.
A coverage-dependent activation energy function was used to describe the desorption of NH 3 in Reaction 7: HA kinetics model parameters were optimized in GT-POWER to minimize the sum of squared error between the experimental data for fractional coverages obtained via deconvolution of the TPD profiles and the fractional coverages calculated by the model via a genetic algorithm. Equation (12) and rate expressions of the Reactions 1-3 given in Table 1 were used. The effects of gaseous species (O 2 and H 2 O) concentrations during HA were neglected since the HA experiments were carried out with a feed which includes O 2 and H 2 O to represent the general cases of HA in an engine ATS, and also not to further complicate the model. For this purpose, the HA aging sets given in Figure 4b,d,f, with durations of 109, 60 and 50 h (i.e., totals of 46, 25 and 21 h isothermal HA, heating and cooling ramps and for NH 3 adsorption and TPD experiments) for 650, 700 and 800 • C aging sets, respectively, were utilized. Thus, not only the isothermal parts of HA but also the associated heating and cooling ramps were considered during the development of the HA kinetics model and its kinetic parameter optimization process.
The parameters of rate expressions of Reactions 4-11 were also optimized in GT-POWER using the experimental NH 3 adsorption/TPD data of the degreened catalyst via a genetic algorithm. Parameters of the rate expressions of Reactions 12 and 13 were obtained using the experimental data of NH 3 oxidation for the degreened catalyst while keeping the parameters of Reactions 4-11 constant. Then, the HA kinetics model and the model for the adsorption and desorption of NH 3 was linked through the utilization of Equations (1)- (14), and along with Reactions 1-11, to describe the NH 3 adsorption and desorption after a certain HA treatment. The HA kinetics model was linked to the models of NH 3 adsorption, desorption and oxidation through the utilization of Equations (1)- (14) and Reactions 1-13.

Conclusions
An accurate HA kinetics model of the changes in the fractional coverages of the four NH 3 storage sites of the Cu-CHA with respect to HA time and temperature was developed based on NH 3 -TPD profile deconvolutions. The main conclusions are as the following: • HA kinetics model successfully captured the changes in the fractional coverages of the NH 3 storage sites, which could be identified with their NH 3 -TPD peak centres observed at 317, 456 and 526 • C, namely ZCu1, ZCu2 and ZB, respectively, for the degreened catalyst. The fractional coverages of ZCu2 and ZCu1 initially increased and decreased, respectively, with HA time after isothermal HA at 650 • C, which was well represented by the HA kinetics model. The major effects of HA were the decrease in the fractional coverages of ZCu2 and ZB and the increase in the fractional coverage of ZCu1 with further increases in both HA time and temperature. These effects were also well captured by the model; • HA kinetics model was then linked with a 4-site transient kinetic model of adsorption, desorption and oxidation of NH 3 developed for degreened Cu-CHA. This linkage was possible via incorporating the changes in the fractional coverages of the NH 3 storage sites to the coverages of the active sites within the reactor model through the utilization of HA kinetics, while keeping the turnover rate constants and the activation energies associated with reactions occurring on the degreened catalysts unchanged;

•
The kinetic model is able to describe the kinetics of NH 3 adsorption and isothermal desorption of NH 3 , TPD of NH 3 and NH 3 oxidation of the degreened and hydrothermally aged Cu-CHA up to 46, 24 and 21 h at 650, 700 and 800 • C, respectively. The decreases in high temperature NH 3 oxidation performance upon HA, at 650 • C for 46 h and at 700 • C for 24 h were very well predicted by the model; • Both the HA methodology and the method used to create a link between the HA kinetics and the reaction models developed in this study could be extended to a variety of other catalytic systems for the prediction catalytic activity after HA. Along this line, the development of an HA model capable of describing the NH 3 -SCR performance after mild HA is on-going in our laboratory.