Abstract
The kinetics of the thermal desorption of CO2 adsorbed on zeolite 13X were obtained using a differential thermogravimetric analyser under two different carrier gas conditions. The varying heating rates were set as 8, 12, 16, and 20 K min−1, respectively. The desorption activation energy of the physisorption sites for this experiment evaluated by an integral method without prediction of the reaction order ranged from 12.15 to 14.12 kJ mol−1 (CO2 as the carrier gas) and 43.32 to 50.42 kJ mol−1 (Ar as the carrier gas), respectively. The desorption activation energy of the chemisorption sites ranged from 57.95 to 58.53 kJ mol−1 (CO2 as the carrier gas) and 74.02 to 79.92 kJ mol−1 (Ar as the carrier gas), respectively.
Introduction
CO2 is considered to be the main greenhouse gas responsible for global warming, formed primarily through combustion processes in power plants and through burning of flammable gases, such as methane, in gas flares in crude oil extractions sites (Cavenati et al., 2005a, 2005b; Chen et al., 2014). Carbon capture and storage (CCS) technologies are increasingly viewed as a potential solution for carbon sequestration (Ordorica-Garcia et al., 2010). Liquid amine sorption of CO2 is one of the main sequestration methods but has drawbacks, such as high equipment corrosion rate, high energy consumption in regeneration, and a large absorber volume required. Adsorption processes with solid sorbents, such as zeolites, have been suggested and studied to overcome many of these drawbacks through their greater stability over a wide range of temperatures and pressures, lower energy requirements, and lower costs. Consequently, these types of materials are being increasingly reported and investigated as practical methods for CO2 capture (Cavenati et al., 2004; Huang et al., 2007; Ko et al., 2002; Liang et al., 2009; Walton et al., 2006). Present studies of physical adsorbents are listed in Table 1.
Present studies of physical adsorbents.
Studies on the adsorption of CO2 on zeolite 13X have been reported regarding the influence of sorbent characteristics (Shen et al., 2013), gas composition (Mulgundmath et al., 2012; Siriwardane et al., 2003), pressure (Siriwardane et al., 2003), and temperature (Siriwardane et al., 2005). However, the existing literature on the desorption of CO2 from 13X is scarce and limited. Desorption characteristics play a significant role in the selection of appropriate adsorbents for the CCS process and also provide guidance for the exploration of adsorption and desorption mechanisms. An adequate understanding of desorption kinetics is useful for further investigations of zeolites as practical CO2 adsorbents, offering reference for the CCS process at a low regeneration energy cost. Thermogravimetric analysis (TGA) and differential thermogravimetric analysis (DTG) have been widely used to study thermal desorption in a molecular sieve. The key kinetic parameters, such as desorption activation energy, reaction order, and pre-exponential factor, could be obtained by mathematical calculations under a non-isothermal condition (Pérez-Maqueda et al., 2002). Several authors have studied the activation energy by the integral method (El-Sayed and Mostafa, 2014) and differential method (Yang et al., 2001) (or ‘DTG curve fitting method’).
In this study, we conducted experiments under two different carrier gases, i.e. argon and CO2, respectively. The parameters of non-isothermal desorption of CO2 from zeolite 13X under two types of carrier gases were examined. Furthermore, the mechanism of adsorption and desorption performances over 13X, as affected by the pore size distribution (PSD) and adsorption sites of the sorbents, is discussed.
Materials and methods
Experimental materials
Adsorption experiments were carried out with a TA Q-500 chemisorption instrument. The commercial zeolite 13X-APG adsorbent used in this study was purchased from Shanghai Jiuzhou Chemicals Co., Ltd. The CO2 adsorbate (99.999%) was obtained from Beijing Qianxi Gas Co. Ltd, and the Ar shielding gas (99.999%) was acquired from Praxair (China) Investment Co. Ltd. The specific surface and porosity analyser Autosorb-1 manufactured by Quantachrome Inc. was used to obtain the PSD of zeolite 13X. Liquid nitrogen was used to maintain the temperature of the adsorption environment at 77.3 K.
Experimental method
Characterization of zeolite 13X
The characterization of the 13X sorbents was obtained by Autosorb-1. Prior to the determination, the samples were maintained at 393 K and 10−5 Torr for 12 h for complete desorption. The specific surface areas were determined by the Brunauer–Emmett–Teller (BET) method based on the adsorption isotherms of N2 at 77 K. The non-local density functional theory was used to obtain the PSD.
Thermogravimetric experiment
An isothermal process was set for 30 min, in which 45 ± 2 mg of the 13X sample was placed on a platinum pan and loaded with CO2 for saturation adsorption. After that step, for the purge stage, the sample was treated, respectively, with CO2 (condition1) or Ar (condition2) at a flow rate of 40 ml·min−1 at the adsorption temperature (293 K) for at least half an hour. Next, the thermogravimetric non-isothermal experiments were performed from 293 to 700 K at four heating rates (β = 8, 12, 16, and 20 K·min−1) purged with the same gas as the previous stage. The curves of mass loss versus temperature (TG) were subsequently recorded.
Summit-differential method
Generally, the desorption rate under non-isothermal conditions could be expressed by the following equation
The uniform kinetic reaction is obtained according to the Polanyi–Wigner equation (Cvetanović and Amenomiya, 1972), widely used in the temperature-programmed desorption process, and could be written as follows
In this experiment, the temperature T is a function of time t since heating occurs in a linear fashion (T0 = initial temperature)
Hence, t = (T − T0)/β, and combining equations (1) to (3), the desorption rate equation becomes
As the DTG curve is a differential of the TG curve, for any point on the curve in this experiment, θ can be expressed as follows
Rewriting equation (8) as
Combining equations (7) and (9), the reaction order is shown as
Taking the logarithm of equation (7), it can be expressed as follows
Equation (10) gives the expression of n, so equation (11) could be rewritten as
Results and discussion
Texture properties of 13X
The PSD of zeolite 13X, which was obtained by the BET method, is illustrated in Figure 2 and includes micropores (0.7–2 nm) that consist of 52% of the total pore volume, and mesopores (2–50 nm) that consist of 48% of the total pore volume. This finding confirmed that 13X is a microporous molecular sieve (Lillo, 2005; Rege et al., 2000). The adsorption capacity of CO2 on zeolite 13X is large as a result of the strong adsorption capacity of micropores and capillary condensation in the mesopores.
Desorption characteristic
Adsorption could be divided into physisorption and chemisorption on the basis of different forces between adsorbents and adsorbates, which have different thermal behaviours during the desorption process (Groß, 2009). Physisorption, based on Van der Waals’ (VDW) force, reaches adsorption equilibrium extremely quickly and desorbs at a lower temperature. Chemisorption, i.e. ‘activated adsorption’, is the result of chemical bond force and the adsorbates are difficult to escape. From the perspective of heat, the adsorption heat released by physisorption is usually similar to the liquefaction heat of gaseous adsorbate, while the adsorption heat of chemisorption, reaching the order of chemical reaction heat, is considerably larger than the former.
Figures 3 and 4 are the combined TG and DTG profiles showing the thermal desorption characteristics of CO2 on zeolite 13X at a heating rate of 8 K·min−1 under CO2 and Ar as carrier gases, respectively. The weight loss process of the CO2-loaded sample can be divided into two main regions (Klepel and Hunger, 2005): escape at low temperature (293–450K); outgas at high temperature (>450 K). As shown in Figures 3 and 4, the removal of very light components is followed by the first region of weight loss at ∼400 K. Much of the desorption occurred in the first region of weight loss, a result of the thermal breaking of weak bonds formed by VDW and the unstable adsorption equilibrium.
Meanwhile, although the amount of desorption of the first zone is significantly different under the two carrier gas conditions, the amount of the second zone remains essentially unchanged. In the DTG curve, the summit position of the curve describes the temperature at which the maximum weight loss rate occurs. The DTG curve of the samples during the desorption stage shows two distinct peaks (which are represented by a noticeable change in the slope of the TG curve). Based on the temperature range at which CO2 has been previously observed to escape, the DTG peaks can be assigned as follows: one between 293 and 450 K and one between 450 and 600 K, written as PL and PR. The highest desorption rate of PL is significantly greater than that of PR. Since the temperature range for the two peaks is comparable, this means that under the same β the total desorption taking place in the temperature range of PL is always greater than that of PR. Figure 3 shows that in the case of CO2 as the carrier gas, the separation of the two desorption peaks is not obvious, since the carrier gas introduced is the same as the loaded gas and is equal to an indirect partial pressure applied on the sample. Meanwhile, as seen from Figure 4, the demarcation of the two zones is trenchant.
Adsorption thermodynamics shows the relationship between the adsorption–desorption process and pressure, as well as temperature. As shown in Figure 5, the curve is the critical state of adsorption and desorption which is generally written as the adsorption–desorption equilibrium state. When the pressure and temperature are Pe and Te, respectively, point e is on the equilibrium curve. Point e′ is in the adsorption zone along the pressure rise to P′ without variation of temperature and cannot escape from the adsorbent. The temperature should be set slightly higher than T′ for the desorption of adsorbates.
The equilibrium adsorption amount has an enhancement as a result of the partial pressure increase of the component. By comparing Figures 3 and 4, a conclusion can be obtained in that elevating the pressure is a non-positive factor for desorption and leads to an increase of the peak height and end temperature of PL; meanwhile, it has an undetectable effect on PR. The surface science tells us that the adsorbates captured by chemisorption could only escape under a sufficiently high temperature and are particularly unaffected by pressure change (Groß, 2009).
The desorption rates versus temperature under different heating rates are shown in Figures 6 and 7. The desorption amounts of the first zone under CO2 atmosphere were higher than those of Ar, corresponding to larger mass losses. The desorption process on a heterogeneous surface during a TPD experiment has been described according to the multiadsorption centre model (Dumesic et al., 1993), which could be used to explain the formation of two peaks in the desorption DTG curves. Two hypotheses have been proposed which explain the presence of the two desorption peaks in this experiment (Yang, 2002): (1) there are two separate active sites for CO2 adsorption on the surface of zeolite 13X, with one clearly exhibiting a stronger adsorption than the other; (2) diffusion towards the positive and negative directions occurred between the surface-adsorbed molecules and sub-layer molecules. It has been found that both adsorption peaks remain when β is varied, so it is safe to ignore any possible diffusion between the adsorbed CO2 multilayers. We propose that the first hypothesis is more likely to occur under the aforementioned conditions. Therefore, for physisorption, we believe it occurs primarily in the micropores along with the amount of capillary condensation in the mesopores. In summary, it could be concluded that PR is chemisorbed and PL is physically adsorbed from another point of view (Siriwardane et al., 2005).
Kinetic parameters study
Summit-differential method
As Figure 8 shows, the slope of the fitted straight line of the plot represents –E/R, which can be used to evaluate the activation energy. Adding this E to equation (10) will derive the reaction order n. The pre-exponential factor A can be determined from the intercept of the line with the vertical axis.

Schematic diagram of Si.

PSD for the sorbent of 13X.
Equation (12) shows that E, n, and A can be obtained by plotting
Table 2 shows the key influential parameters for different desorption conditions. With the same adsorption carrier gas, the α value of PL or PR has a regular variation with increasing heating rate whereas the α value of the total desorption process remains almost unchanged.
Key influential parameters.
From Table 2 and Figure 3, PL and PR will overlap when CO2 is the carrier gas. Yang (1983) presented the following points in the desorption spectrum analysis: (1) the contributions of the two peaks to the overlapping parts are generally unequal; (2) the non-overlapped side of the peak is less affected or even negligible by the other one. Therefore, the desorption activation energy should be calculated from the range close to the summit of each peak, which is distant from the overlapping region.

TG and DTG curves with CO2 as the carrier gas.

TG and DTG curves with Ar as the carrier gas.

Adsorption–desorption process.

Desorption rates with CO2 as the carrier gas under different β.
Results and reaction mechanism
The three-dimensional structure of zeolite 13X, a faujasite-type zeolite, is shown in Figure 9(a). It is composed of alumina tetrahedra and silicon oxygen tetrahedron consisting of oxygen bridge according to a certain spatial structure to form a uniform aperture crystal. The basic structural unit of 13X is a beta cage, and adjacent beta cages form the molecular sieve by the six-angle prism connection. The alumina tetrahedra show electronegativity, and Na+ neutralizes the negative charge. The Na+ in 13X occupies three positions (Feuerstein and Engelhardt, 1998) as shown in Figure 9(b). The SI and SI′ positions are not exposed and cannot be used to interact with the adsorbed molecules after the dehydration. The SII and SIII positions are exposed and can participate in adsorption.

Desorption rates with Ar as the carrier gas under different β.

Schematic diagram of the calculation method of kinetic parameters.
Siriwardane found that the adsorption mode of CO2 on 13X can be categorized into two types: (a) physisorption of CO2 occurs with CO2 in linear orientation by ion–dipole interaction

Structure diagram of zeolite 13X.

Chemisorption of CO2 on 13X.
Table 3 shows the key kinetic triplets obtained from the summit-differential method. Eave is the average value of calculated values of E at different heating rates. As the sample size (the number of experiments is 4, and therefore the degree of freedom is 4–1 = 3) is small, T distribution should be used to avoid the inaccuracies caused by the size of the sample in the error analysis. The significance level is selected as 0.05 (confidence probability is 95%) to calculate confidence intervals. Referring to the T distribution list,
Key kinetic triplets by the summit-differential method.
For the first zone, there is a perfect linearity of
For the second zone, the value of E under condition 1 is significantly lower than the value under condition 2. One of the possible reasons is that the demarcation of the first zone and the second zone is indistinct, causing difficulty in determining the exact starting point of the calculation and the result is presented as an average of desorption activation energies of both physisorption and chemisorption. Meanwhile, regardless of which gas is used as the carrier gas, the calculation results of E in the second region fluctuate within a certain range. The reaction order n can be considered as 2 in the second zone and 1 in the first zone, which means that the desorption mechanisms corresponding to the two zones are different. Although the value of A in the second zone is higher than that in the first zone, the desorption rate in the second zone being faster than the first zone cannot be explained due to the different desorption amount between the two zones, and the value of A used to determine the value of the desorption rate assumes the same desorption amount. Taking the heating rate as the only variable, the value of A increases along with the growing heating rate which represents an increase in the desorption rate. This is because the same number of CO2 molecules needs to escape from the surface of molecular sieves.
It can be seen that the unconsidered selection of the existing data of E in other articles is not advisable because the atmosphere conditions in actual engineering are not completely consistent with the experimental conditions. For different engineering applications, there are great discrepancies between the atmosphere and adsorption sites where the main adsorption occurs.
Conclusions
A pressure difference can be formed between the adsorption–desorption process by the method using two different gases as the carrier gas (one gas should be the adsorbate), which plays a supplementary role in the determination of adsorption and desorption mechanisms. A conclusion could be drawn that the physisorption and chemisorption of CO2 on zeolite 13X exist simultaneously.
The summit-differential method without a need of estimating n could be used to obtain the key kinetic triplet, i.e. activation energy, pre-exponential factor, and reaction order, accurately. For the calculations of the desorption activation energy of different engineering applications, it is necessary to adjust the carrier gas atmosphere on the basis of various actual conditions and adsorption mechanisms to obtain the proper desorption activation energies.
Footnotes
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The project was supported by the Fundamental Research Funds for the Central Universities (FRF-BD-16-009A).
