Nobuyoshi Koga

Loïc Favergeon

Satoki Kodani

Loic Favergeon

Impact of atmospheric water vapor on the thermal decomposition of calcium hydroxide: a universal kinetic approach to a physico-geometrical consecutive reaction in solid-gas systems under different partial pressures of product gas

come

Impact of atmospheric water vapor on the thermal decomposition of calcium hydroxide: a universal kinetic approach to a physico-geometrical consecutive reaction in solid-gas systems under different partial pressures of product gas

Introduction

The impact of the partial pressure of a gaseous reactant or a product in the reaction atmosphere on the kinetics of the reversible reactions in solid-gas systems, which may be thermal decomposition of solids or solid-gas reactions, is one of the intrinsic problems for the kinetic description.

A(s) ⇄ B(s) + C(g) (1) 49 One possible way to formulate the kinetic equation with consideration for the impact of the partial pressure of C(g), p(C(g)), is to introduce an accommodation function, [50][51][52][53][START_REF] Koga | SIn Thermal Analysis of Micro, Nano-and Non-Crystalline Materials[END_REF][START_REF] Koga | Handbook of Thermal Analysis and Calorimetry[END_REF] a(p(C(g)), P eq (T)), into the fundamental kinetic equation that describes the reaction rate of a single-step reaction as a function of temperature (T) and the fractional reaction (α).

d d = exp - C g , (2)
where A, E a , and R are the Arrhenius preexponential factor, activation energy, and gas constant, respectively.

P eq (T) is the equilibrium pressure of the reversible reaction of eqn (1). The function f(α) is the kinetic model function that describes the changes in the reaction rate as the reaction progresses at a constant temperature. The equilibrium pressure, P eq (T), for the forward reaction of eqn ( 1) is expressed by the following equation:

= exp !- ∆ # $ ° & (3)
where ∆ r G° is the standard Gibbs energy for the reaction. At a constant temperature, eqn (2) is rewritten as follows:

d d = ' C g , (4)
where k is the rate constant. Various forms of a(p(C(g)), P eq (T)) have been derived and used for describing reaction kinetics in solid-gas systems. 1-10, 41, 56,57 Among these accommodation functions, a simple function has been widely used in the practical kinetic analysis for the thermal decomposition of solids with consideration for the influence of p(C(g)). [1][2][3]7 C g , = 1 -

) * + , -. / (5)
The accommodation function of eqn ( 5) is derived simply by considering the contribution of the reverse reaction in the overall reaction rate. 1 A comparable form of the accommodation function has also been derived by focusing on the component reaction steps of gas absorption on the intermediate layer, 42 gaseous diffusion via the solid-product layer, [START_REF] Searcy | [END_REF] and the formation of an active solid product. 3 The accommodation function of eqn ( 5) is also derived more generally based on the transition state theory. 48 However, application of eqn (2) with the accommodation function of eqn (5) for kinetic calculations sometimes results in the varying Arrhenius-type plots for the respective reaction processes under different p(C(g)) values; these plots shift systematically to higher temperatures with the increase of p(C(g)). 7 These results imply that the equation cannot universally describe the kinetic behaviors of the reaction under different p(C(g)) values. For a rigorous description based on the fundamental kinetic equation, which assumes a single-step reaction characterized by a constant E a value, alternative forms for the accommodation function should be utilized to formulate a kinetic equation universally applicable to the process in solid-gas systems under different p(C(g)) conditions. On the other hand, apparent variations of calculated E a with p(C(g)) have been subjected to correlation analysis, and arguments have been made for finding a theoretical background of possible variations in E a with p(C(g)). 47 As stated previously, thermal decomposition of solids and solid-gas reactions are heterogeneous processes that are constrained by the physico-geometrical properties of the reaction mechanism. [START_REF] Galwey | Thermal Decomposition of Ionic Solids[END_REF]52,58 In many cases, mass and heat transfer phenomena largely influence on the apparent kinetic behavior. [59][60][61][62] Since a surface of reactant solid is the most reactive site, the overall reaction proceeds via the physico-geometrical consecutive processes comprising the surface and subsequent internal reactions. In turn, these reactions are characterized by surface nucleation and growth processes, as well as advancement of the reaction interface toward the center of the original solid reactant. [START_REF] Galwey | Thermal Decomposition of Ionic Solids[END_REF]52,58 In some cases, a significant induction period (IP) is also observed prior to the surface reaction (SR). [63][64][65][66][67][68][69] Such heterogeneous characteristics of the reaction and any significant impact of the mass and heat transfer phenomena are responsible for the numerous observations induced by atmospheric gases. Therefore, in addition to having a universal kinetic description for these processes under various p(C(g)) conditions, investigation into the impact of p(C(g)) on each component physico-geometrical reaction step must be explored separately in order to gain further insights into the fundamental problems of reaction kinetics that are often encountered in solid-gas systems. To this end, in this study, we focus on the thermal decomposition of Ca(OH) 2 under various water vapor pressure values, p(H 2 O).

Ca(OH) 2 (s) ⇄ CaO (s) + H 2 O (g) (6) This reaction is a well-known reversible process, and thermal decomposition rate is significantly suppressed by the presence of atmospheric water vapor. 11,[70][71][72] The overall mass-loss behavior observed during thermal decomposition has been explained via the contracting-geometry-type reaction model. 72 At the same time, the two-step physico-geometrical consecutive processes of surface and internal reactions was predicted using kinetic analysis. 11,70 In addition, prereaction structural distortion is a reported trigger for thermal decomposition reaction, [73][74][75] in which the characteristic temperatures of structural distortion are strongly dependent on atmospheric p(H 2 O). Previously reported trends in the thermal decomposition of Ca(OH) 2 indicate the characteristics of the physico-geometrical consecutive process and the impact of atmospheric p(H 2 O) on the kinetics of each component reaction step, which are desirable features for the purpose of the present kinetic study.

There is a need to describe universal kinetics for the thermal decomposition processes under various p(H 2 O) conditions and to separately interpret the impact of p(H 2 O) on each physico-geometrical reaction step through logically designed kinetic analyses for the systematically recorded kinetic data under controlled p(H 2 O) values.

Experimental

Sample Characteristics

Reagent grade Ca(OH) 2 (special grade, 499.9% in metal base, Wako Pure Chem. Ind.) was used without any further purification or pulverization. In our previous study, 11 which was focused on the kinetic modeling for the thermal decomposition of Ca(OH) 2 in an inert gas atmosphere and thermally induced carbonation of Ca(OH) 2 in a CO 2 atmosphere, the sample was characterized using powder X-ray diffraction (XRD), Fourier-transform infrared (FTIR) spectroscopy, and simultaneous thermogravimetry and differential thermal analysis (TG-DTA). The sample exhibited partial carbonation, in which a mixed phase of Ca(OH) 2 and CaCO 3 was formed, and the purity of Ca(OH) 2 was determined to be 94.30 mass% from the mass-loss values during the thermal decomposition of Ca(OH) 2 and CaCO 3 . The morphological characteristics of the sample powders were also examined using a scanning electron microscopy (SEM). The specific surface area of the sample was determined using the single-point Brunauer-Emmett-Teller (BET) method. It was revealed that the sample was the agglomerate of submicron-sized particles, and the specific surface area was 3.41±0.03 m 2 g -1 . The same sample was used for this study, and all measurements were performed during the same period alongside the aforementioned kinetic studies for the thermal decomposition of Ca(OH) 2 in an inert atmosphere and the thermally induced carbonation of Ca(OH) 2 in a CO 2 atmosphere.

TG measurements under various water vapor pressures

Approximately 3.00±0.05 mg of the Ca(OH) 2 sample was weighed into a platinum pan that was 5 mm in diameter and 2.5 mm deep. The weighed sample was set in a TG-DTA instrument (TG8120; Rigaku) that had been constructed using a horizontal thermobalance with a differential mechanism. Water vapor pressure in the reaction atmosphere was controlled using an instrument (HUM-1, Rigaku) according to the following procedures. 76 Hot water in the controlled temperature range of 323-343 K was circulated around the furnace tube of the TG-DTA instrument prior to the measurement run in order to avoid possible water vapor condensation when the atmospheric gas with a higher p(H 2 O) is later introduced into the furnace tube. The sample was initially heated to a predetermined temperature within the range of 323-343 K at a heating rate (β) of 5 K min -1 . During the linear heating phase, dry N 2 was introduced at a rate of 400 cm 3 min -1 at the forefront of the furnace tube. At the same time, a purge gas consisting of dry N 2 was also introduced at the back of the balance system at a rate of 50 cm 3 min -1 . These gases were discharged from a gas line positioned at the end of the furnace tube, at a spot that was close to the joint of the furnace tube and the balance system. While maintaining the sample within the aforementioned temperature range, the gas introduced at the forefront of the furnace tube was switched to a wet

Results and Discussion

Impact of the water vapor pressure on the thermal decomposition process

Irrespective of the atmospheric p(H 2 O) values applied for the thermal decomposition of the sample, two separate mass-loss processes were observed under linearly increasing temperature conditions, as in the dry N 2 gas flow conditions reported previously 11 (see Fig. S2 in the ESI †). These mass-loss steps originated from the thermal decomposition of Ca(OH) 2 and CaCO 3 . 11,77 The first mass-loss step attributed to the thermal decomposition of Ca(OH) 2 was subjected to kinetic studies focused on the impact of atmospheric water vapor. Fig. 1 illustrates the characteristics of the thermal decomposition of Ca(OH) 2 under isothermal conditions and the impact of atmospheric water vapor on the reaction. In the presence of water vapor, the IP was clearly observed as a certain time period without any mass changes after the sample has reached the pre-determined reaction temperature and before the start of the mass-loss process (Fig. 1(a)). This was not evident in our previous study of the same sample in a stream of dry N 2 gas. 11 Once mass loss was initiated, the mass-loss rate gradually increased, reaching its maximum midway through the reaction, followed by gradual deceleration. Generally, the mass-loss process for the thermal decomposition of Ca(OH) 2 in a stream of wet N 2 gas was characterized by a smooth curve with a sigmoidal shape and was vastly different in shape from that observed in a stream of dry N 2 gas. 11 In a stream of dry N 2 gas, the mass-loss process was characterized by two partially overlapping reaction steps, which were attributed to the thermal decomposition on Ca(OH) 2 surfaces that were activated via partial carbonation of the sample and of the main structure of the particles under investigation. The SR appeared as the first kinetic step in a stream of dry N 2 gas with the contribution of approximately 0.2 to the overall process was disappeared in a stream of wet N 2 gas. This phenomenon indicates a restraining influence of atmospheric water vapor on the SR process and explains the appearance of the IP. Isothermal mass-loss curves recorded at a fixed temperature under different p(H 2 O) values (Fig. 1(b) and (c)) clearly indicated that the IP and reaction time for the mass-loss process were prolonged under the higher p(H 2 O) values. Therefore, it is evident that atmospheric water vapor exerts a restraining influence on both the SR and the phase-boundary-controlled reaction (PBR) at the reaction interface that is advancing toward the center of reacting particles of Ca(OH) 2 .

Fig. 2 shows the mass-loss curves of the thermal decomposition of Ca(OH) 2 under different p(H 2 O) conditions, which were achieved by linearly heating the sample at β = 5 K min -1 . Nonisothermal mass-loss curves shifted systematically to higher temperatures with increasing p(H 2 O) (Fig. 2(a)), from which a qualitative relationship between atmospheric p(H 2 O) and its restraining influence on the overall mass-loss process during thermal decomposition was observed as a normal effect of p(H2O) for a reversible solid-gas reaction system in view of chemical equilibrium. Both the reaction initiation temperature expressed by the extrapolated onset temperature (T eo ), which is the temperature at the intersection point of baseline and the tangent line of the initial acceleration period, and the peak top temperature (T p ) of the DTG curve shifted to higher temperatures with increasing p(H 2 O) (Fig. 2(b)). However, the difference between T eo and T p decreased with increasing p(H 2 O), which was accompanied by a systematic rise in the height of the DTG peak. Systematic changes in the reaction temperature interval and the shape of the DTG peak under nonisothermal conditions at a selected b were indicative of changes in the overall kinetic behavior of the reaction; this could be linked to the effect of p(H 2 O) and its diverse impact on specific reaction stages, such as SR and the subsequent reaction at the reaction interface.

Kinetics of the induction period

Figure 3 shows the isothermal mass-loss curves for the thermal decomposition of Ca(OH) 2 at various temperatures in a stream of wet N 2 gas with three different p(H 2 O) values. Irrespective of the p(H 2 O) value, it was noted that the IP was systematically prolonged with a reduction in the reaction temperature. A systematic decrease in the reaction rate with a reduction in the reaction temperature was also observed for subsequent mass-loss processes.

In order to record the isothermal mass-loss curves for kinetic calculations, the reaction temperature range for the isothermal measurements was forcibly shifted to a higher temperature region with increasing p(H 2 O) values to compensate for the restraining influence of p(H 2 O).

The kinetic data for the IP was extracted from the isothermal mass-loss curves using an empirical procedure that defined the duration time of IP (t IP ) as the time interval from the moment when the sample achieved the desired measurement temperature to the time at which the sample indicated a mass-loss value of 0.23%. This value corresponded to the fractional reaction a of 0.01 for the thermal decomposition of Ca(OH) 2 . Fig. 4 8/ 45 (7) where A IP and E a,IP are the apparent preexponential factor and activation energy for the IP process, respectively.

f(α IP ) is a kinetic model function that describes the changes in the reaction rate at a constant temperature as a function of the fractional conversion noted during the IP process, α IP . Taking the logarithms of eqn (7), the conventional Arrhenius-type plot of ln(t IP -1 ) versus T -1 can be applied to the series of data (T, t IP ). 65,66,68,69,[78][79][80] Application of the Arrhenius-type plot based on eqn ( 8) is justified when zero-order kinetics with f (α IP ) = 1 is assumed for the IP process 78 or when the value of a IP with the average conversion rate (t IP -1 ) is assumed to be invariant for IP processes at the various temperatures being fulfilled by an isoconversional relationship. 65,66,68,69,80 In Fig Introducing an accommodation function in the form of eqn ( 5) is a possible primary approach. [1][2][3]7 H > O , = 1 - S1 in the ESI. † Changes in P eq for the reversible reaction of eqn ( 6) with temperature variations, as calculated using a thermodynamic database (MALT2, Kagaku Gijutsu-Sha), 81,82 are shown in Fig. S4 (ESI †), along with (T, p(H 2 O))

) @ A B , -. / (10)
values applied to mass-loss measurements that were carried out under isothermal conditions. For the data points recorded under specific p(H 2 O) conditions, the value of a(p(H 2 O), P eq (T)), as determined using eqn (10), tended to vary because of the temperature-dependent variations of P eq . Changes in the slope of the modified Arrhenius plot under specific p(H 2 O) conditions (Fig. S3(b), ESI †) compared to the observations made from the conventional Arrhenius plot (Fig. S3(a), ESI †) appeared to be due to corrections in the overall rate, which have been determined according to the accommodation function in eqn (10). However, the modified Arrhenius plots under various p(H 2 O) conditions indicated that there were different linear correlations with different slopes (Fig. S3(b), ESI †), as in the case of the original Arrhenius plots (Fig. S3(a), ESI †). The results indicate that the IP processes taking place within vastly different temperature ranges due to the effect of p(H 2 O) cannot be universally described by applying the conventional accommodation function in eqn (10). If we adhere to the fundamental assumption, which is a single-step reaction characterized by a constant E a value for the processes under different p(H 2 O) conditions, it becomes clear that an alternative accommodation function that can be used as a universal kinetic description is sorely needed. As it relates to finding alternative forms of the accommodation function, possible frameworks were sought by considering the elementary steps of the IP process. The IP process for reactions in the solid state is empirically defined as the time required for activating the reaction at a constant temperature. In the case of thermal decomposition of solids, which produces both solids and gases as products, the reaction is usually initiated on the reactant surfaces. Therefore, the physico-chemical processes of nucleation on the surfaces of the reactant solid should be considered when looking at the kinetics of the IP process. [63][64][65][66][67][68][69] According to the classical nucleation theory, [START_REF] Garner | Chemistry of the Solid State[END_REF][START_REF] Young | Decomposition of Solids[END_REF][START_REF] Tompkins | Treates on Solid State Chemistry[END_REF][START_REF] Brown | Reactions in the Solid State[END_REF][START_REF] Galwey | Thermal Decomposition of Ionic Solids[END_REF][START_REF] Delmon | Introduction a La Cinetique Heterogene[END_REF][START_REF] Barret | Cinetique Heterogene[END_REF] various heterogeneous processes may take part in the nucleation process on the surfaces of reactant crystals as elementary steps. For the sake of simplicity, one rate-limiting step was selected from these elementary steps, and the other steps were assumed to be at equilibrium by applying rate-limiting step approximation. For example, using Kröger's notation, [START_REF] Kröger | The Chemistry of Imperfect Crystals[END_REF] five successive steps, as listed in Table 1, could be considered as elementary steps that took part in surface nucleation processes for the thermal decomposition of divalent metal hydroxides to form a divalent metal oxide and water vapor.

M(OH) 2 (s) ⇄ MO(s) + H 2 O(g) (12)
The equilibrium constant (K i ) for the elementary step i is expressed as listed in the third column of Table 1, whereas the fourth column shows kinetic equations for the overall formation rate (v i ) of the nuclei at a constant temperature, as expressed using the rate constant (k i ) of the rate-limiting step i. All kinetic equations contain the p(H 2 O) component. P eq (T) is also included in all kinetic equations because the overall equilibrium constant (K) for the reversible reaction in eqn ( 6) is expressed by the product of K i for all elementary steps i and is equivalent to the equilibrium pressure P eq (T).

G = G 0 G > G H G I G J 0/L = exp - ∆M °8/ = (13)
At a constant temperature, each K i has a certain value. Therefore, the transformation rate v IP of the overall IP process which is controlled by process i = 1, 2, or 4 can be reduced to the following equation:

N 45 = ' 45 O 1 - ) @ A B , -. / (14)
Similarly, for the overall process controlled by the process i = 5, the rate equation is expressed as

N 45 = ' 45 OO 0 ) @ A B P Q1 - ) @ A B , -. / L R T ≥ 1 (15)
where k IP ' and k IP '' are the overall rate constants. In eqn ( 13)-( 15), the unit of all pressure terms is in atm. The rate equation for the overall process controlled by the elementary step i = 3 is reduced to eqn ( 14) and ( 15) with n = 1, when K 4 is much larger and smaller than p(H 2 O), respectively. Reduction of the rate equation for the overall process, which is controlled by the elementary step i = 4, to eqn ( 14) requires the assumption that the value of K 1 K 2 K 3 K 5 1/n is much larger or smaller than unity. The overall rate expressed by eqn ( 14) is influenced by the relative value of p(H 2 O) with reference to P eq (T) in the form of the conventional accommodation function in eqn (10). An alternative form of the accommodation function may be derived from the rate equation in eqn (15), in which the overall rate depends on the product of two p(H 2 O) components, which are composed of the reciprocal of the n th powered p(H 2 O) and the n th powered relative value of p(H 2 O) with reference to P eq (T), respectively.

Based on eqn ( 14) and ( 15), a generalized form of the accommodation function was proposed, which would universally describe the kinetic behavior of IP processes under different p(H 2 O) conditions.

H > O , = 0 ) @ A B V C1 - ) @ A B , -. / W F ( 16
)
where a and b are constants. P 0 is the standard pressure, which is introduced to express all pressure terms in the unit of Pa. and vacancy (V B@ • ) at the interstitial and OH sites, respectively

OH B@ ⇄ OH XYZ O + V B@ • G 0 = ;OH XYZ O <;V B@ • < N 0 = ' 0 !1 - H > O &
(2) Formation of hydrogen ( H XYZ • ) and oxygen (O B@ O ) defects at the interstitial and OH sites, respectively

OH B@ ⇄ H XYZ • + O B@ O G > = ;H XYZ • <;O B@ O < N > = ' > !1 - H > O &
(3) Consumption of the interstitial defects and formation of H 2 O absorbed on the surface (S)

OH XYZ O + H XYZ • + S ⇄ H > O-S G H = ;H > O-S < ; S <;H XYZ • <;OH XYZ O < N H = ' H G 0 G > G I G J 0/L G I + `>a !1 - H > O & (4) Desorption of water molecule from the surface H > O-S ⇄ H > O g + S G I = ; b < H > O ;H > O-S < N I = ' I G 0 G > G H G J 0/L 1 + G 0 G > G H G J 0/L !1 - H > O &
(5) Consumption of the defects and formation of the O elements

TV B@ • + TO B@ O ⇄ TO B G J = 1 ;V B@ • < L ;O B@ O < L N J = ' J G 0 G > G H G I H > O L C1 -! H > O & L F
Having (a, b) = (0, 1), the accommodation function is comparable to the conventional accommodation function in eqn (10). The dependence of the overall rate on p(H 2 O), as expressed by eqn (11), can be modified by applying eqn ( 16) with a = b = n (n ≥ 1, integer). For a more empirical approach, the values of a and b can be determined numerically through kinetic calculations based on eqn (11) with the accommodation function in eqn ( 16) to obtain the most appropriate kinetic expression for IP processes under various p(H 2 O) conditions. Such an approach may lead to empirical (a, b) values in which it is difficult to interpret the physico-chemical significance of the results.

Even so, the empirical (a, b) values that were determined through the aforementioned kinetic calculations were useful for evaluating the relevance of the newly introduced accommodation function and for revealing the role of each p(H 2 O) component in eqn (16) in relation to its influence on the kinetic description for the reaction under investigation.

All the data of the IP process recorded at different temperatures under various p(H 2 O) conditions were universally analyzed using the modified Arrhenius plot based on eqn (11) with the accommodation function in eqn (16). 2 summarizes the kinetic parameters that were determined using a modified Arrhenius plot that was based on eqn (11) with the accommodation function of eqn (16).

Kinetics of the mass-loss process

Fig. 6 shows the mass-loss curves obtained for the thermal decomposition of Ca(OH) 2 at different β values under various p(H 2 O) conditions. For each of the selected p(H 2 O) values, the mass-loss curves shifted systematically to a higher temperature range with β, as was generally observed for kinetic processes. The systematic shift in the mass-loss curves to higher temperatures with increasing p(H 2 O), which was demonstrated previously in Fig. 2, was observed as a general trend irrespective of the selected β value. Kinetic analysis of the mass-loss process during the thermal decomposition of Ca(OH) 2 was performed using mass-loss curves recorded under isothermal (Fig. 3) and linear nonisothermal (Fig. 6) conditions. The fundamental kinetic equation in differential form, which assumed a(p(H 2 O),P eq (T)) = 1 or constant in the expression of eqn (2), was employed for the preliminary kinetic approach.

de d1 = exp - 6 7 8/ (17)
where A and E a are the apparent preexponential factor and activation energy for the mass-loss process, respectively. The fractional reaction (α) was defined as the relative mass-loss value at time t with reference to the overall mass-loss value. The kinetic model function f(α) describes a-dependent changes in the overall reaction rate.

For the mass-loss process in the thermal decomposition of solids, the functional form of f(α) was theoretically derived by considering the probability of the reaction initiation via surface nucleation, the geometry of the reaction interface advancement, and the rate-limiting step of this interface advancement, that is, chemical reaction or diffusion-controlled reaction. [START_REF] Galwey | Thermal Decomposition of Ionic Solids[END_REF]52,58,[START_REF] Khawam | [END_REF] In addition to obtaining an ideal f(α) based on the physico-geometry of the reaction, empirical functions f(α) have also been applied in order to accommodate the actual rate behavior of the reaction, 85 for example, the Sestak-Berggren model SB(m,n,p): 50,86,87 The perfect fit of the experimentally resolved kinetic curves that were determined using the empirical f(α) serves as a way to avoid any undesired distortions of the Arrhenius parameters that originated from an inappropriate f(α), [88][89][90] because the distortions usually result in an empirical linear correlation between E a and ln A known as the kinetic compensation effect. [91][92][93][94] When the reaction is well defined as the single-step reaction and the mechanistic feature does not change during the course of the reaction and within the temperature region for kinetic analysis, the problem of mutual dependence among f(α), E a , and ln A can be avoided by separately determining the E a value at a restricted a based on an isoconversional kinetic relationship. 95,96 Taking the logarithms in eqn (17), ln de d1 = ln; < - Among the data points (T, dα/dt) at a selected a, the plot of ln(dα/dt) versus T -1 represents a linear correlation in which the E a value at particular a is calculated from the slope of the linear plot. The differential isoconversional method based on eqn ( 18) is known as the Friedman plot 97 and is applicable to the kinetic data recorded under any temperature profiles. [98][99][100][101] Fig. S5 in the ESI † shows the results of the isoconversional analysis that was applied to the mass-loss process for the thermal decomposition of Ca(OH)

f(α) = α m (1-α) n [-ln(1-α)] p .
It was also expected that the introduction of an appropriate accommodation function greatly improve the kinetic description for the reactions under different p(H 2 O) conditions, because similar trends in the variation of the E a value as the reaction progressed were noted in Fig. S5 Therefore, possible alternative forms of the accommodation function of a(p(H 2 O), P eq (T)) are again required for the mass-loss processes of the thermal decomposition of Ca(OH) 2 .

The overall mass-loss process of the thermal decomposition of solids can be interpreted as a physico-geometrical consecutive process of the SR and PBR, which are characterized by nucleation and growth of the solid product on the surface and subsequent movement of the as-produced reaction interface inward the reacting particle, respectively. [START_REF] Garner | Chemistry of the Solid State[END_REF][START_REF] Young | Decomposition of Solids[END_REF][START_REF] Tompkins | Treates on Solid State Chemistry[END_REF][START_REF] Brown | Reactions in the Solid State[END_REF][START_REF] Galwey | Thermal Decomposition of Ionic Solids[END_REF][START_REF] Delmon | Introduction a La Cinetique Heterogene[END_REF][START_REF] Barret | Cinetique Heterogene[END_REF] PBR is accompanied by the diffusional removal of the gaseous product via the surface product layer. Thus, the influence of p(H 2 O) on the interfacial processes of the PBR was further considered for the overall mass-loss process in addition to nucleation processes as illustrated for the IP process. Similar to nucleation during the IP process, PBR during the mass-loss process can also be described by successive elementary steps that occur at the reaction interface with an appreciable thickness. In addition to the formation/consumption of defects and the adsorption/desorption of gaseous products as assumed for surface nucleation, diffusion of the as-generated defect across the phase boundary (reaction interface) should be considered in the interfacial processes. Successive elementary steps as listed in Table 3 may describe the overall process of PBR. The rate equations listed in the fourth column of Table 3 are obtained when one of the elementary steps is assumed to be the rate-limiting step and the other elementary steps are set at equilibrium. These can be reduced to two types of rate expressions in relation to the dependence of the overall rate of PBR (υ PBR ) on p(H 2 O).

N 5fg = ' 5fg O Q1 - ) @ A B , -. / L R T = 0 > hi 1 (21)
N 5fg = ' 5fg OO 0

) @ A B P Q1 -

) @ A B , -. / L R T = 0 > hi ≥ 1 ( 22
)
where the unit of all pressure terms in eqn ( 21) and ( 22) is in atm. Eqn ( 21) is derived when the elementary step of i = 1, 2, 3a, or 3b is assumed as the rate-limiting step, which is equivalent to eqn ( 14) when n = 1 and can have n = 1/2 when the reaction rate is controlled by diffusion of the as-formed defects via the reaction interface. Eqn ( 21) is also derived from the rate equations for steps i = 4 and 5 when p(H 2 O) is much smaller than K 5 and K 1 K 2 K 4 K 6 1/n is much smaller or larger than unity, respectively. Eqn (22), which is expressed by the product of two p(H 2 O) components with exponent n in each component, is derived from the rate equation of i = 5 if K 5 is much smaller than p(H 2 O) and from i = 6. The frameworks for eqn ( 21) and ( 22) are comparable with eqn ( 14) and ( 15) derived for the IP process, respectively.

Subsequently, it was found that H 2 O(g) produced at the reaction interface diffused to the surface of the reacting particle through the surface product layer and was discharged into the reaction atmosphere via the desorption from the outer surface of the product layer. For processes controlled by the diffusion of gaseous product through the surface product layer, an accommodation function comparable to that mentioned in eqn ( 5) had already been derived by Searcy and Beruto. [START_REF] Searcy | [END_REF] Therefore, all possible rate expressions concerning the effect of p(H 2 O) were related to the analytical form of the accommodation function described in eqn (16). Fig. 7 shows the results of the kinetic analysis using the modified

Friedman plot with the accommodation function of eqn (16). Application of eqn ( 20) combined with eqn (16) provided the universal linear correlation of the Friedman plot for mass-loss processes under various p(H 2 O) conditions (Fig. 7 16) did not significantly contribute to the universal kinetic description of the mass-loss process at the selected a value under various p(H 2 O) conditions and that the contribution did not change over the course of the reaction. Thus, the exponent a for the first p(H 2 O) component, (P 0 /p(H 2 O)) a , was again responsible for the universal kinetic description for the isoconversional relationship of the mass-loss processes. However, the exponent a changed systematically as the reaction advanced, exhibiting an initial decrease from approximately 3.2 to 2.0 in the first half of the process (approximately 0.01 ≤ α ≤ 0.50), a nearly constant value of approximately 2.0 in the second half of the process (approximately 0.50 ≤ α ≤ 0.85), and a gradual increase from approximately 2.0 to 2.5 at the end (approximately 0.85 ≤ α ≤ 0.99). Variations in the a value were indicative of the changing impact of p(H 2 O) over the course of the reaction, which was accompanied by the variations in the contribution of the diffusional removal process of the as-produced water vapor. It was noted that the change in the a value was synchronized with the variation of the apparent E a value over the course of the reaction (Fig. 7(c)). Therefore, the kinetic approach that was developed on the basis of the modified

Friedman plot with the accommodation function of eqn ( 16) revealed that there were changes in the physico-geometrical reaction mechanism and the impact of p(H 2 O) over the course of the reaction.

Phenomenologically, the existence of the IP and the isothermal mass-loss process characterized by the sigmoidal mass-loss curve were positive evidence of physico-geometrical consecutive processes comprising the IP, SR, and PBR steps. The initial decrease in the E a value during the first half of the reaction and the constant E a value noted during the second half, as revealed through the isoconversional kinetic approach, have also been observed for this type of reaction. 65,66,68,69 It is expected that a rigorous kinetic solution for the physico-geometrical consecutive process provides further insights into the impact of p(H 2 O) on each component reaction step.

OH B@ ⇄ OH XYZ O + V B@ • G 0 = ;OH XYZ O <;V B@ • < N 0 = ' 0 !1 - H > O &
(2) Formation of hydrogen (H XYZ • ) and oxygen (O B@ O ) defects at the interstitial and OH sites, respectively

(internal) OH B@ ⇄ H XYZ • + O B@ O G > = ;H XYZ • <;O B@ O < N > = ' > !1 - H > O & (3a) Diffusion of (OH XYZ O ) from internal to external of the reaction interface (interface) OH XYZ@XYZ # O ⇄ OH XYZ@ kZ # O G H = lOH XYZ@ kZ # O m lOH XYZ@XYZ # O m = 1 N H = ' nH o p G 0 G > G q 0/L 0/> r1 -! H > O & 0/> s (3b) Diffusion of (H XYZ • ) from internal to external of the reaction interface (interface) H XYZ@XYZ # • ⇄ H XYZ@ kZ # • G Ht = lH XYZ@ kZ # • m lH XYZ@XYZ # • m = 1 N Ht = ' nHt o p G 0 G > G q 0/L 0/> r1 -! H > O & 0/> s (
OH XYZ O + H XYZ • + S ⇄ H > O-S G I = ;H > O-S < ; S <;H XYZ • <;OH XYZ O < N I = ' I G 0 G > G J G q 0/L G J + `>a !1 - H > O &
(5) Desorption of water molecule from the external surface (S) of the reaction interface

(external) H > O-S ⇄ H > O g + S G J = ; S < H > O ;H > O-S < N J = ' J G 0 G > G I G q 0/L 1 + G 0 G > G I G q 0/L !1 - H > O & (6)
Consumption of the defects and formation of the O elements

nV • OH + nO' OH ⇄ nO O TV B@ • + TO B@ O ⇄ TO B G q = 1 ;V B@ • < L ;O B@ O < L N q = ' q G 0 G > G I G J H > O L C1 -! H > O & L F

Kinetic modeling of the consecutive IP-SR-PBR under various water vapor pressures

Kinetic modeling of the consecutive physico-geometrical processes for the surface nucleation and the subsequent reaction interface shrinkage was first proposed by Mampel for solid-state reactions under isothermal conditions. 102 Practical kinetic analysis based on the integral form of Mampel's model has been studied by several researchers in order to gain further insights into the real kinetic feature of various solid-state reactions. 103-106 Favergeon et al.

discussed in detail the constitutional equations for this type of physico-geometrical consecutive process within the context of nucleation and isotropic or anisotropic growth models. 107 Several examples of the physico-geometrical consecutive process were analyzed kinetically using the principles from the nucleation and anisotropic growth models. [108][109][110][111] The integral kinetic behavior, which was characterized by a sigmoidal mass-loss curve under isothermal conditions, is superficially similar to those obtained for the reactions comprising the bulk nucleation and growth of nuclei, as formalized by the Johnson-Mehl-Avrami (JMA) model, [112][113][114][115] and the homogeneous-like autocatalytic model. 50,86,87 For the kinetic analysis conducted under isothermal conditions, reactions regulated by JMA and autocatalytic models could be treated as the pseudo-single-step reaction described by one rate constant;

on the other hand, in reactions regulated by Mampel's model, two rate constants for SR and PBR should be evaluated from a single kinetic curve at a constant temperature. Because of the difficulty of conducting kinetic calculations for the experimental kinetic curves that were obtained under isothermal conditions based on Mampel's model, in many cases, the SR-PBR processes were analyzed using the JMA or auto-catalytic models with an acknowledgment as the empirical approach, regardless of the vastly different physico-geometrical features assumed while formalizing the respective kinetic models. However, difficulties in performing kinetic calculations based on Mampel's model have been largely alleviated by the increased quality of the kinetic data recorded through the use of the modern thermoanalytical instruments and advancements in software technology. For practical kinetic calculations based on Mampel's model, Ogasawara and Koga recently formalized the differential kinetic equations for physico-geometrical consecutive processes comprising IP-SR-PBR(n), which occur under isothermal conditions: this was done by assuming a zero-order rate process for IP, a first-order rate process for SR, and a PBR with different geometrical dimensions (n). 68 The differential kinetic equations obtained for the consecutive IP-SR-PBR(n) process under isothermal conditions are listed in Table S2 in the ESI. † Application of the kinetic calculation generates the chance to separately discuss the impact of p(H 2 O) on each component step of the physico-geometrical reaction, as well as gain further insights into the mechanistic feature of the reaction.

The isothermal kinetic data recorded at numerous temperature settings under different p(H 2 O) conditions (Fig. 3)

were used to conduct kinetic analysis based on the IP-SR-PBR(n) model. The rate constants, k IP (= 1/t IP ), k SR , and k PBR(n) for IP, SR, and PBR with the interface shrinkage dimension of n, respectively, were simultaneously optimized so as to minimize the square sum F of the residue when fitting the experimental (dα/dt) exp versus t by the calculated (d a /dt) cal versus t based on the equations listed in Table S2 (ESI †) through nonlinear least-squares analysis according to the Levenberg-Marquardt algorithm.

u = ∑ Q de d1 kw,x - de d1 y z,x R > { x|0 (23)
where N is the number of data points in each kinetic curve. For each optimization run, the initial value for k IP was determined using the apparent Arrhenius parameters obtained via the conventional Arrhenius plot for IP processes under each p(H 2 O) condition (Table S1 in the ESI †). For k PBR(n) , the initial value was derived with reference to the average E a value that was evaluated using the conventional Friedman plot within the range of 0.5 ≤ α ≤ 0.8 and analyzing an experimental master plot drawn by assuming a constant E a value within the corresponding a range. 95,96,[99][100][101] After setting the initial values of k IP and k PBR(n) , the order of the initial k SR value was determined by graphically comparing the experimental and calculated kinetic curves. The temperature dependence of the rate constants for each reaction step was separately analyzed on the basis of the conventional Arrhenius plot. ln' = ln -6 7 8/ (24) where k can be k IP , k SR , or k PBR (2) . In Fig. S7 (ESI †), there are conventional Arrhenius plots applied to the k IP (=1/t IP ), k SR , and k PBR (2) values determined for the overall process of the thermal decomposition of Ca(OH) 2 at different temperatures under various p(H 2 O) conditions. For each reaction step, that is, IP, SR, and PBR(2), the conventional Arrhenius plot detailed the different linear relationships that existed under different p(H 2 O) conditions; this was already observed when the conventional Arrhenius and Friedman plots were applied to the IP process (Fig. S3(a) in the ESI †) and to the overall mass-loss process (Fig. S5(a) in the ESI †), respectively. For all reaction steps, the slope of the Arrhenius plot systematically rose with increasing p(H 2 O). Table S4 (ESI †) lists the apparent Arrhenius parameters that were determined for each reaction step under a particular p(H 2 O) condition.

Thus, introducing an accommodation function, a(p(H 2 O),P eq (T)), in the Arrhenius plot is a possible means of establishing the universal kinetic description for each physico-geometrical reaction step that occurs at different temperatures and under various p(H 2 O) conditions. ln C } E ) @ A B ,, -. / F = ln - conditions. The apparent Arrhenius parameters that were determined using the modified Arrhenius plot with the accommodation function of eqn (10) are also summarized in Table S4 (ESI †). Although the accommodation function of eqn (10) could not be used to universally describe the rate data obtained under different p(H 2 O) values, the effect of the correction appeared more extensively for reactions under the higher p(H 2 O) value, as changes in the slope of the Arrhenius plot, and thus the apparent E a value. The largest a value was observed for the IP step. This was expected because the IP step was triggered via exposure to the reaction atmosphere and was not influenced by any self-generated water vapor pressure. Smaller a values observed for the SR and PBR(2) steps were indicative of apparently the rate behaviors that were relatively less sensitive to p(H 2 O) in comparison to those found in the IP step. However, when interpreting the apparent rate behavior, it must be considered that the actual reaction sites of the SR and PBR(2) steps were not directly exposed to the reaction atmosphere. This situation was more distinguished for the PBR(2) step, since the reaction interface was covered by a surface layer of the solid product. In addition, these reaction steps were accompanied by the evolution of water vapor; therefore, the actual water vapor pressure at the reaction sites positioned within the reacting particles could not be precisely determined. Thus, the impact of atmospheric p(H 2 O) on the SR and PBR steps was the result of changes in the actual partial pressure of water vapor at the reaction interface, which were induced by changes in the atmospheric p(H 2 O), even if the evolution rate of water vapor during the reaction under isothermal conditions was much lower than that observed under the nonisothermal conditions. Since there were challenges in directly measuring or precisely calculating the partial pressure of water vapor at the reaction interface, systematic changes in the a value as the reaction step advanced could be used as a semi-quantitative measure that would empirically explain the observed variations in the apparent rate behavior via the effect of atmospheric p(H 2 O). The relatively constant b values for all reaction steps may have resulted from the same p(H 2 O)/P eq (T) ratio among the corresponding data points recorded at specific (T, p(H 2 O)) conditions in all the reaction steps. From the aforementioned results of the modified Arrhenius plots with the accommodation function of eqn (10) (Table S4 and Fig. S8 in the ESI †), it was expected that the second p(H 2 O) component, 1 -(p(H 2 O)/P eq (T)) b , would make slight corrections to the slope of the individual Arrhenius plot under a specific p(H 2 O) value. However, the effect of the second component in eqn (16) was smaller than the value that determined in the case of the conventional accommodation function of eqn (10) because the b value was found to be larger than unity. Therefore, the single Arrhenius plot for the rate constants of each reaction step, which had been calculated for the kinetic data recorded under different p(H 2 O) conditions, was again recognized via the contributions of the first p(H 2 O) component in eqn (16).

It was observed that the apparent E a values for each reaction step evaluated from the modified Arrhenius plot with the accommodation function of eqn (16) systematically decreased as the reaction step advanced from IP to PBR (2) via SR. The trend agreed with the aforementioned behavior patterns that were observed via formal kinetic analyses for the IP and the overall mass-loss processes. It was also noted that the preexponential factors that were determined for each reaction step decreased as the reaction step advanced in accompanying with the decrease in the E a values, which followed a generally observed trend known as the kinetic compensation effect. [91][92][93][94] For the thermal decomposition of the current Ca(OH) 2 sample, the results of the kinetic analysis in relation to the impact the Arrhenius parameters for a specific reaction step and in a different variation trend of these correlations as the reaction step progresses. Classifications of the correlation patterns and these variation trends may provide information necessary for establishing a rigorous theoretical and practical basis for the comprehensive understanding of reaction kinetics in solid-gas systems under various gaseous atmospheres.

Conclusions

In the presence of atmospheric water vapor, the thermal decomposition of Ca(OH) 2 exhibits vastly different kinetic behavior patterns from those obtained in a stream of dry N 2 gas; this is characterized by the appearance of the IP and the sigmoidal mass-loss processes at constant temperatures. Reaction rates of the IP and mass-loss processes at a selected temperature are inhibited as the p(H 2 O) value rises. Under a linearly increasing temperature condition at a specific β value, the mass-loss curve systematically shifts to higher temperatures with rising p(H 2 O) values. Introducing an accommodation function related to p(H 2 O) into the conventional kinetic equation is one possible way to universally describe the apparent kinetic behavior that is dependent on α, T, and p(H 2 O). The accommodation function of eqn (16), which was derived from principles of classical nucleation and interface reaction theories, allows for a universal kinetic description for the IP and the mass-loss processes when applying formal kinetic calculations that use the Arrhenius plot and Friedman plot, respectively.

In our study for thermal decomposition of Ca(OH) (16). Similarly, the mass-loss process recorded under isothermal and linear nonisothermal conditions in various p(H 2 O) conditions was analyzed using a modified Friedman plot with the accommodation function of eqn (16), resulting in a single Friedman plot at each α. The a value in eqn (16), which was optimized at different α, decreased systematically from approximately 3.2 to 2.0 in the first half of the mass-loss process; it was subsequently maintained at approximately 2.0, followed by a detectable increase to approximately 2.5 in the final stage of the mass-loss process. On the other hand, the optimized b value was relatively constant within the wide range of α (0.04 ≤ α ≤ 0.97), with the average value of 1.28±0.07. The apparent E a value for the mass-loss process varied and followed the similar trend to the a values, exhibiting a systematic decrease from approximately 540 kJ mol -1 to 375 kJ mol -1 in the first half of the mass-loss process, plateauing at approximately 375 kJ mol -1 during the second phase, and then experiencing a detectable increase to approximately 420 kJ mol -1 at the end.

Therefore, changes in the physico-geometric reaction mechanism and the impact of p(H 2 O) over the course of the reaction could be deduced from the results of these modified kinetic plots.

The appearance of the IP and subsequent sigmoidal mass-loss behavior were interpreted as outcomes of the mechanistic feature of the reaction, which is characterized by physico-geometrical consecutive processes comprising the IP, SR, and PBR steps. Variations in the apparent E a value observed for the mass-loss process were also indicative of consecutive SR and PBR steps. Nonlinear least-squares analysis for overall thermal decomposition process under isothermal conditions with the IP-SR-PBR models indicated that the two-dimensional interface shrinkage occurred during the PBR step with the most statistically significant coefficient. Subsequently, establishing a universal kinetic description for each physico-geometrical reaction step at different T under various p(H 2 O) conditions was made possible by analyzing the optimized rate constants for each reaction step using the modified Arrhenius plot with the accommodation function of eqn (16). The kinetic parameters comprising the Arrhenius parameters and the optimized exponents (a, b) in eqn (16) for each physico-geometrical reaction step could be determined using the kinetic analysis based on the IP-SR-PBR models and subsequent application of the modified Arrhenius plot with the accommodation function of eqn (16). Kinetic results based on the physico-geometrical reaction mechanism provided the opportunity to discuss the impact of S2-S4).

Corresponding Author *Tel./fax: +81-82-424-7092. E-mail: nkoga@hiroshima-u.ac.jp

N 2

 2 gas with a controlled p(H 2 O), flowing at a rate of approximately 400 cm 3 min -1 . After the measurement system was stabilized in a stream of wet N 2 gas, mass-change measurements were performed by heating the sample to 873 K at different β values (1 ≤ β ≤ 10 K min -1 ). Similarly, the sample was heated to different pre-determined temperatures (T) at β = 10 K min -1 and held at these temperatures while the mass-change trace was recorded. In order to obtain the mass-loss traces for the kinetic calculations, three different p(H 2 O) values were selected at approximately 0.65, 1.85, and 5.65 kPa. The details of the experimental setup involving the instrument calibration and the control of p(H 2 O) values in the reaction atmosphere are described in Section S1 in the ESI, † together with the examples of the TG measurement procedures performed under controlled p(H 2 O) conditions (Fig. S1, ESI †).

Figure 1 .

 1 Figure 1. Mass-loss characteristics for the thermal decomposition of Ca(OH) 2 under isothermal conditions under a stream of wet N 2 gas (400 cm 3 min -1 ) and the p(H 2 O)-dependent changes: (a) TG-DTG curves at 663 K under p(H 2 O) = 1.85 kPa, (b) TG curves at 658 K under p(H 2 O) = 0.65 and 1.85 kPa, and (c) TG curves at 694 K under p(H 2 O) = 1.85 and 5.65 kPa.

Figure 2 .

 2 Figure 2. Mass-loss curves of the thermal decomposition of Ca(OH) 2 under nonisothermal conditions at β = 5 K min -1 under a stream of wet N 2 gas (400 cm 3 min -1 ) and the p(H 2 O)-dependent changes: (a) TG-DTG curves at various p(H 2 O) values and (b) changes in the extrapolated onset temperature (T eo ), peak top temperature (T p ), and peak height at T p of the DTG curve depending on p(H 2 O).

Figure 3 .

 3 Figure 3. TG-DTG curves for the thermal decomposition of Ca(OH) 2 under isothermal conditions at various temperature settings under a stream of wet N 2 gas (400 cm 3 min -1 ) as characterized by different p(H 2 O) values: (a) p(H 2 O) = 0.65 ± 0.01 kPa, (b) p(H 2 O) = 1.85 ± 0.01 kPa, and (c) p(H 2 O) = 5.65 ± 0.04 kPa.

  . S3 (ESI †), conventional Arrhenius plots that were applied to the IP process are shown for the thermal decomposition of Ca(OH) 2 in a stream of wet N 2 gas as characterized by different p(H 2 O) conditions. The IP processes under the different p(H 2 O) values were described by individual Arrhenius plots, while each plot exhibits statistically significant straight line (γ 2 > 0.98) (Fig.S3(a), ESI †). The slope of the plots tended to vary depending on p(H 2 O), which was indicative of a systematic rise in the slope with increasing p(H 2 O) values; thus, the increases in apparent E a,IP and ln A IP values were noted. The apparent kinetic parameters determined using the conventional Arrhenius plots are summarized in TableS1in the ESI. †The basic kinetic equation of eqn(7) does not consider the impact of p(H 2 O) on the kinetic behavior of the IP. Thus, kinetic descriptions for the IP processes under different p(H 2 O) conditions by the individual Arrhenius plots are inevitable and results in a distorted set of apparent kinetic parameters. Introducing an accommodation function a(p(H 2 O), P eq (T)), composed of p(H 2 O) and P eq (T) for the reversible reaction of eqn(6), into eqn (7) is one possible way of obtaining the kinetic description considering the p(H 2 O) as one of the parameters controlling the kinetic behavior of the IP process.

Figure 4 .

 4 Figure 4. Duration of the IP, t IP , for the thermal decomposition of Ca(OH) 2 at various temperature settings under different p(H 2 O) values.

  Practically speaking, the accommodation function corrects the apparent overall reaction rate by simply assuming the linear relationship between the restraining influence of p(H 2 O) and the relative magnitude of p(H 2 O) with reference to P eq (T). By introducing the accommodation function, the Arrhenius plot for the IP process is modified as the plot of ln[t IP -1 •a(p(H 2 O),P eq (T)) -1 ] versus T -1 . Arrhenius plot based on eqn (11) with the accommodation function of eqn (10) cannot universally describe the influence of p(H 2 O) representing individual Arrhenius plots for the IP processes under different p(H 2 O) values (Fig. S3(b), ESI †). Similar results of the application of the modified Arrhenius plot were reported for the thermal decomposition of CaCO 3 under various partial pressures of CO 2 . 7 For the modified Arrhenius plot for each IP process under specific p(H 2 O) conditions, it was noted that the slope of each plot decreased slightly from the corresponding original Arrhenius plot (Fig. S3(a), ESI †); however, there was a general trend toward a rising slope with p(H 2 O), which did not deviate from the results observed for conventional Arrhenius plots. The apparent values of (E a,IP , ln[A IP f(α IP )]) determined by the modified Arrhenius plots for each process under different p(H 2 O) conditions are also summarized in Table

  During the kinetic calculations, the (a, b) values in the accommodation function of eqn(16) were optimized in accordance with the Levenberg-Marquardt algorithm to give the best linearity of the modified Arrhenius plot by conducting a linear regression analysis. Fig.5shows the results of the said modified Arrhenius plot in two cases: one with the restriction of a = b and the other without it. Irrespective of the restriction of (a, b) values, all data points of the IP processes could be universally described by a single Arrhenius plot that exhibits statistically significant linearity (γ 2 > 0.99). Under the restriction of a = b, exponents in the accommodation function of eqn (16) were determined to be a = b = 4.74. When the restriction was removed, thus allowing the different values of a and b, the linearly of the modified Arrhenius plot improved slightly without any practical change in the slope, and the values of a and b values were determined to be (a, b) = (4.53, 1.57). The a value was not significantly influenced by restrictions of a = b, whereas the b value decreased dramatically and approached unity. The second p(H 2 O) component in eqn (16), expressed by (1 -(p(H 2 O)/P eq (T)) b ), was equivalent to the conventional accommodation function of eqn(10), when b was equal to unity. As previously mentioned, introducing the conventional accommodation function of eqn(10) did not significantly improve the kinetic description of the IP process under different p(H 2 O) conditions; it only resulted in the slight changes in the slope of the Arrhenius plot for IP processes under specific p(H 2 O) conditions. Therefore, it was empirically understood that the first p(H 2 O) component in eqn(16), as expressed by (P 0 /p(H 2 O)) a , played a predominant role in accommodating the influence of p(H 2 O) on the IP process in the equation needed to determine the universal kinetic description of the process under various p(H 2 O) conditions. Superficially, the (P 0 /p(H 2 O)) a component corresponds to the expression for the contribution of a negative catalyst in the homogeneous-like kinetic equation. Table

Figure 5 .

 5 Figure 5. Modified Arrhenius plots with the accommodation function of eqn (16) applied to the IP of the thermal decomposition of Ca(OH) 2 at different temperatures under various p(H 2 O) values.

Figure 6 .

 6 Figure 6. TG-DTG curves for the thermal decomposition of Ca(OH) 2 under linear nonisothermal conditions at different β values in a stream of wet N 2 gas (400 cm 3 min -1 ) characterized by different p(H 2 O) conditions: (a) p(H 2 O) = 0.66 ± 0.01 kPa, (b) p(H 2 O) = 1.80 ± 0.02 kPa, and (c) p(H 2 O) = 5.52 ± 0.07 kPa.

  (b) (ESI †) among the reactions under different p(H 2 O) conditions. By introducing the accommodation function, a(p(H 2 O), P eq (T)), the modified form of the Friedman plot the conventional accommodation function of eqn (10) did not significantly improve the results of the conventional Friedman plot as shown in Fig. S6 in the ESI. † The modified Friedman plot represented individual linear correlations with different slopes and intercepts in ln[(dα/dt)/(1-p(H 2 O)/P eq (T))]versus T -1 coordinate for the mass-loss processes under each p(H 2 O) value (Fig.S6(a), ESI †). The apparent E a values calculated from the modified Friedman plots decreased slightly in comparison to those determined using conventional Friedman plots, where a more impactful effect of the accommodation function was observed for mass-loss processes recorded at higher p(H 2 O) values. However, introducing the accommodation function of eqn(10) did not achieve a universal kinetic description for the mass-loss processes under different p(H 2 O) conditions and did not affect the trends in the variation of the apparent E a values as the reaction progressed (Fig.S6(b), ESI †).

  (a)), when the appropriate values of a and b were evaluated for each Friedman plot at different α by allowing the different values for a and b. Variations of the optimized exponents (a, b) in the accommodation function of eqn (16) as the reaction progresses (Fig. 7(b)) tell us that there are the changes in the impact of p(H 2 O) during the course of the reaction. The exponent b was nearly constant within the range of a (0.04 ≤ α ≤ 0.97), with the average value of 1.28 ± 0.07. Constancy in the exponent b, which was close to unity, indicated that the second p(H 2 O) component in the accommodation function of eqn (

4 )

 4 Consumption of the interstitial defects and formation of H 2 O absorbed on the external surface (S) of the reaction interface (external)

Figure 7 .

 7 Figure 7. Results of the modified Friedman plots with the accommodation function of eqn (16) as applied to the mass-loss process of the thermal decomposition of Ca(OH) 2 under isothermal and linear nonisothermal conditions at different p(H 2 O) values: (a) modified Friedman plot with (a, b) = (1.95, 1.31) at α = 0.5, (b) a and b values at different α values, and (c) E a values at different α values.

Fig. 8

 8 Fig.8shows typical results of the optimization run as exemplified by the kinetic curve recorded at T = 668 K under p(H 2 O) = 1.86 kPa; this was indicative of the most satisfactory fit obtainable when two-dimensional shrinkage of the reaction interface (n = 2) was assumed for the PBR step (Fig.8(b)). Irrespective of p(H 2 O), the most statistically relevant fit was observed when assuming the IP-SR-PBR(2) model. Therefore, the k IP (= 1/t IP ),

  shown in Fig. S8 (ESI †), introducing the most widely used accommodation function of eqn (10) into eqn (25) did not dramatically improve the individual Arrhenius plots for the processes under different p(H 2 O)

Figure 8 .

 8 Figure 8. Typical results of nonlinear least-squares analysis for the appropriate differential and integral kinetic curves obtained for the thermal decomposition of Ca(OH) 2 (T = 668 K, p(H 2 O) = 1.86 kPa) on the basis of IP-SR-PBR(n) model: (a) n = 1, (b) n = 2, and (c) n = 3.

Figure 9 .

 9 Figure 9. Modified Arrhenius plots with the accommodation function of eqn (16) which describe the temperature dependence of the corrected rate constants for the IP, SR, and PBR(2) steps of the thermal decomposition of Ca(OH) 2 under various p(H 2 O) values: (a) IP, (b) SR, and (c) PBR(2).

p(H 2 O

 2 ) on each reaction step and enabled us to compare its effect on different reaction steps in connection with the mechanistic features for each. The kinetics of the thermal decomposition of Ca(OH) 2 under atmospheric water vapor were characterized by systematic reductions in the Arrhenius parameters. The exponent a was shown to decrease as the physico-geometrical reaction advanced from IP to PBR(2) via SR, whereas the exponent b was shown to remain relatively constant. Introducing the accommodation function of eqn (16) into formal kinetic analysis and kinetic analysis based on the physico-geometrical consecutive reaction model allowed for the universal kinetic description of the thermal decomposition of solids under different partial pressures of the gas that contribute to the reaction. Kinetic parameters, comprising the Arrhenius parameters and the optimized exponents (a, b), and those variation behavior patterns with the advancement of the overall reaction or of physico-geometrical reaction step provide useful information of characterizing the impact of the partial pressure of the gas on the kinetics of the reaction. It also hints at the possibility that there are differences among the thermal decomposition reactions for various solids. The kinetic approach to the thermal decomposition of solids under various gaseous atmospheres demonstrated by our study is expected to contribute toward providing further comprehensive theoretical foundation for the kinetic understanding of reactions in solid-gas systems.FootnotesElectronic Supplementary Information Electronic supplementary information (ESI) is available: S1: Experimental setup and calibration of instrument (FigureS1); S2: Sample characterization (FigureS2); S3: Kinetics of induction period (FiguresS3 and S4, TableS1); S4: Kinetics of mass-loss process (Figures S5 and S6); S5: Kinetic modelling of the consecutive IP-SR-PBR process under different water vapor pressure (Figures S7 and S8, Tables

Table 1 .

 1 Elementary steps that take part in the nucleation processes in the crystal surfaces of metal hydroxides during the IP and the possible kinetic equation for the chemical process of the IP derived by assuming the selected elementary step as the rate-limiting step

	Elementary step i (1) Formation of hydroxide defect (OH XYZ O )	Reaction	Equilibrium constant, K i	Overall rate v for the process controlled by the elementary step i

Table 2 .

 2 Apparent kinetic parameters for the IP process, as determined according to the modified Arrhenius plot

	a = b	4.74	4.74	731.8 ± 15.9	105.8 ± 2.9	0.9960
	a ≠ b	4.53	1.57	711.2 ± 15.0	103.1 ± 2.7	0.9962

based on eqn

(11)

with the accommodation function of eqn (

16

)

Condition a b E a,IP / kJ mol -1 ln[A IP f(α IP ) / s -1 ] -γ a

a Correlation coefficient of the linear regression analysis.

  2 under isothermal and linear nonisothermal conditions at various p(H 2 O) values. In each p(H 2 O) value, the Friedman plot at a constant a exhibited a linear correlation, including the data points recorded under isothermal and linear nonisothermal conditions (Fig. S5(a), ESI †). However, Friedman plots for the mass-loss processes under various p(H 2 O) values appeared at quite different positions in the ln(dα/dt) versus T -1 coordinate with different slopes. The slope of the Friedman plot at

	each a systematically increased with a rise in the p(H 2 O) values, and thus, a larger apparent E a value for the
	reaction under a higher p(H 2 O) value was obtained. However, trends in the variations of the E a value as the
	reaction progressed were comparable among the reactions under various p(H 2 O) values (Fig. S5(b), ESI †),
	thereby representing systematic decreases in the initial part of the reaction (approximately α ≤ 0.3) followed by a
	nearly constant E a value during the main part of the reaction (approximately 0.3 ≤ α ≤ 0.9) and a subsequent rise
	at the end of the reaction (approximately α ≥ 0.9).
	de d1	= exp -	6 7 8/	H > O ,

The results of the conventional Friedman plot clearly indicated that an accommodation function a(p(H 2 O), P eq (T)) should be introduced into the kinetic equation in the form of eqn (2) for describing the influence of p(H 2 O) on the reaction rate of the mass-loss process for the thermal decomposition of Ca(OH) 2 , as was the case for the IP process.

Table 3 .

 3 Physico-chemical processes that participate in the phase boundary reaction during the thermal decomposition of metal hydroxides and the kinetic equation for the overall reaction controlled by the selected elementary step

	Elementary step i	Reaction	Equilibrium constant, K i	Linear advancement rate v for the process controlled by the elementary step i
	(1) Formation of hydroxide defect ( OH XYZ O ) and vacancy ( V B@ interstitial and OH sites, respectively • ) at the	(internal)

Table 4 .

 4 Optimized exponents (a, b) in the accommodation function of eqn (16) and the apparent Arrhenius parameters for each reaction step as determined by applying the modified Arrhenius plot with the accommodation

	function of eqn (16)
	Reaction step	Exponents (a, b) in eqn (16) a b	E a / kJ mol -1	ln (A / s -1 )	-γ a
	IP	4.75	1.86	736.8 ± 15.4	74.8 ± 2.6	0.9965
	SR	3.79	1.65	610.4 ± 11.4	66.5 ± 2.0	0.9972
	PBR(2)	3.36	1.91	539.7 ± 14.0	52.7 ± 2.5	0.9947

a Correlation coefficient of the linear regression analysis.

  2 , the kinetic behavior of the IP process under different T and p(H 2 O) values was described by a single Arrhenius plot, resulting in (E a,IP /kJ mol -1 , ln[A IP f(α IP )/s -1 ]) values of (711.2±15.0, 103.1±2.7), and was accompanied by optimized (a, b) values of (4.53, 1.57) in the accommodation function of eqn

7,23 8/(8)

Acknowledgements

The present work was supported by JSPS KAKENHI Grant Numbers 17H00820 and 16K00966.

of p(H 2 O) on the physico-geometrical consecutive processes comprising IP, SR, and PBR (2) were difficult to correlate with previously reported kinetic results of the reaction in a stream of dry N 2 gas. This was due to the presence of vastly different apparent reaction behavior patterns. As mentioned, the IP step was not experimentally detectable for the reaction in a stream of dry N 2 gas. 11 Mass-loss processes were indicative of two distinguishable peak maxima corresponding to the SR and PBR. Kinetic deconvolution analysis [START_REF] Koga | Handbook of Thermal Analysis and Calorimetry[END_REF]116 based on the cumulative kinetic equation, which was comprised of the conventional kinetic equation of eqn (17), revealed that the contribution of SR and PBR steps was 0.20 and 0.80, respectively. The apparent Arrhenius parameters (E a /kJ mol -1 , A/s -1 ), which were determined for the reaction steps of SR and PBR under isothermal conditions in a stream of dry N 2 gas were (160.3±0.5, (3.02±0.01)×10 10 ) and (181.6± 0.2, (8.79±0.01)×10 11 ), respectively. Comparable values for these contributions and the apparent Arrhenius parameters for each reaction step were also evaluated for thermal decomposition reactions in a stream of dry N 2 gas under isothermal and constant mass-loss rate conditions.

The presence of water vapor in the reaction atmosphere hindered the observation of multiple peak maxima in the mass-loss rate under both isothermal and non-isothermal conditions, and the sigmoidal mass-loss curve that was obtained was the characteristic for general rate behavior patterns under isothermal conditions. The apparent Arrhenius parameters evaluated by the modified Arrhenius plot with the accommodation function of eqn ( 16) for the SR and PBR(2) steps under different p(H 2 O) values were much larger than those calculated for the corresponding reaction steps in a stream of dry N 2 . In addition, the relative magnitudes of the apparent Arrhenius parameters for the SR and PBR steps were opposite to each other for the reactions in a stream of dry N 2 gas and the reactions in a stream of N 2 -H 2 O mixed gas with controlled p(H 2 O) conditions. Differences in the kinetic results obtained for the thermal decomposition of Ca(OH)2 in a stream of dry N 2 gas and wet N 2 gas were the evidence of the impact of atmospheric p(H 2 O) on the mechanistic feature of the reaction.

At the same time, the physico-chemical significance of the Arrhenius parameters that were determined by the conventional Arrhenius-type plots and the modified Arrhenius-type plots with the accommodation function a(p(H 2 O), P eq (T)) were varied. Arrhenius parameters for IP, SR, and PBR(2) reaction steps under various p(H 2 O) conditions, which were established through the introduction of the accommodation function of eqn ( 16), could only be of practical importance when exponents (a, b) in the accommodation function were accompanied.

Although the framework for the accommodation function of eqn ( 16) was derived on the basis of classical nucleation and interface reaction theories, the values of (a, b) could be analytically applied to a universal kinetic description for reactions under different p(H 2 O) conditions by allowing the values obtained to deviate from theoretically meaningful values. In this sense, it was difficult to rigorously discuss the physico-chemical significance of the empirically determined values of (a, b) and, thus, the Arrhenius parameters. Even so, a universal kinetic description of the physico-geometrical consecutive processes of IP, SR, and PBR under different p(H 2 O) conditions, as exemplified by the thermal decomposition of Ca(OH) 2 , provided further insights into the impact of atmospheric p(H 2 O) on each reaction step and the variation of the said impact over the course of physico-geometrical reaction. This universal kinetic description allows us to empirically simulate reaction processes under various temperature and p(H 2 O) conditions and provides a chance to interpret the significance of p(H 2 O) in relation to a physico-geometrical reaction mechanism. Additionally, the application of this empirical kinetic approach to other thermal decomposition of solids, which can be described by the physico-geometrical consecutive process of IP-SR-PBR, may result in a different set of correlations between the exponents (a, b) and