Equilibrium and Kinetic Adsorption of Hydrogen Fluoride onto Zeolite 3A

Equilibrium and kinetic adsorption behaviours of hydrogen fluoride (HF) onto zeolite 3A were investigated under different temperatures (298–338 K) and pressures (0.15–1.1 bar). The HF adsorption isotherms on zeolite 3A were well described by the Dubinin–Astakhov (DA) model. The isosteric heat of HF adsorption was calculated based on the DA model and Clausius–Clapeyron equation, and different kinetic models were used to analyze the HF adsorption kinetic. The results showed HF adsorption isotherms are type I of the IUPAC classification. The maximum adsorption capacities by the DA model are 0.443–0.5631 mg/g. The H s t values decrease with increasing surface loading indicating the surface energy heterogeneity of zeolite 3A. The HF adsorption kinetics on zeolite 3A follows simple first-order kinetics. The calculated film diffusion parameter B t values show nonlinear characteristics with time. The order of magnitude of interparticle HF diffusion coefficient D c is 10−9 m2/s. The higher pressure and temperature are in favour of HF diffusion. The intraparticle diffusion curves for HF adsorption on zeolite 3A show quart-linear characteristics, indicating the presence of four consecutive HF adsorption steps including film diffusion, interparticle diffusion, intraparticle diffusion, and surface adsorption. The intraparticle diffusion is a rate-controlling stage.


Introduction
Fluorine-containing chemicals are widely used in atomic energy, metallurgy, machinery, chemical industry, and other industries because of their unique physical and chemical characteristics [1]. As a basic fluorine chemical product, hydrogen fluoride (HF) is widely used in producing refrigerants, fluorine resin, fluorine rubber, and fluorine-containing intermediate and fine chemicals [2]. HF is corrosive and toxic, and it poisons humans and animals through the food chain, destroys the normal metabolism of calcium and phosphorus, inhibits enzyme activity, and affects the nervous system [1]. In the use of fluorine (F 2 ), e.g., uranium enrichment, optics, and material fluorination, HF is often accompanied by F 2 , which reduces the purity of fluorine [3][4][5]. Therefore, HF purification treatment and recovery are indispensable to save resources and protect the environment.
The common methods of collecting and recovering HF are compression condensation, absorption, and adsorption [6]. The compression condensation method can purify HF gas well in the-ory, but this method needs low temperature and consumes more energy [7]. The HF amount treated by the absorption method is huge. However, absorption purification is wet with high requirements for equipment and chemical treatment problems [8]. Collecting HF by adsorption method cannot only remove HF gas but also reduce the emission of toxic gas HF and energy consumption. The adsorption method is mainly based on the HF adsorption performance. Thus, the basic HF adsorption data is required [9]. The frequently used adsorbents for HF adsorption are porous, e.g., sodium fluoride, activated alumina, activated carbon, and silica gel [2]. Zeolite is a commonly used adsorption material in gas separation and gas purification [6]. The application of zeolite can be extended to HF treatment.
The studies on HF treatment by adsorption have been investigated by some scholars. Sodium fluoride (NaF) is a commonly used adsorption material for HF [10][11][12]. NaF reacts with HF to produce NaHF 2 . When the temperature is higher than 250°C, NaHF 2 can be decomposed into NaF and HF. Afzal et al. conducted experiments on HF adsorption with NaF by the fixed bed, and the influence of concentration, temperature, and other parameters on the surface loading has been studied in detail [10][11][12]. The experimental results show that Langmuir and Temkin can better describe the HF adsorption isotherm than Freundlich [11]. Moreover, Liu et al. used NaF to adsorb and remove trace HF in the production of high-purity tungsten hexafluoride, and the HF adsorption by NaF under different temperatures, pressures, and gas flows was measured by orthogonal test [13]. The results showed the suitable adsorption temperature was 20-80°C, and the flow rate has a significant effect on adsorption efficiency. Activated alumina (Al 2 O 3 ) is also used for HF adsorption. Activated alumina reacts with HF to produce AlF 3 and H 2 O. In addition, Wu et al. studied the adsorption characteristics and influencing factors for HF adsorption by Al 2 O 3 , and the relationships between crystal structure, the specific surface area of Al 2 O 3 , and HF adsorption have been considered [14]. Sheng studied the HF adsorption behaviour on activated carbon, calcium oxide, activated alumina, iron oxide, and fly ash, and the Freundlich model described the adsorption isotherms well [15]. In addition to NaF and Al 2 O 3 , some other porous materials are also gradually applied to HF adsorption. Che et al. studied the HF adsorption behaviour on DB18C6/SiO2 composite [16]. Siahooe et al. studied the HF adsorption on carbon nanotubes and found that the HF adsorption density decreased with temperature increase [17]. Wang et al. investigated the HF adsorption behaviour on three metal-organic frameworks (MOFs) and found that MOFs had better HF adsorption performance than commercial Al 2 O 3 and 13X zeolite [18]. Bahrami et al. investigated the HF gas adsorption by three kinds of activated carbon under vacuum and described the experimental data by isotherm adsorption models [19].
Zeolite 3A is an excellent performance adsorbent for gas storage and separation [20]. However, the studies about HF adsorption on zeolite 3A are very lacking. Therefore, the study of the adsorption characteristics and influencing factors of HF on zeolite 3A is necessary. In this study, the equilibrium and kinetic adsorption curves of HF on zeolite 3A at different temperatures and pressures were measured. The HF adsorption isotherms on zeolite 3A were described by the three isotherms model. The DA model and the Clausius-Clapeyron equation were used to calculate the isosteric heat of HF adsorption, and different kinetic models were used to analyze the kinetic process for HF adsorption on zeolite 3A. The adsorption mechanism of HF onto zeolite 3A was evaluated based on kinetics and thermodynamics. 2.2. Characterization of the Zeolite 3A. The compositions were measured by X-ray diffraction (XRD) using a Rigaku Smartlab 9kW with Cu Kα radiation (40 kV and 40 mA). The Fourier transform-infrared (FT-IR) of zeolite 3A spectra was obtained from an FT-IR spectrometer (Thermo Nicolet iS5, USA) using pellets of zeolite/KBr. The surface morphology was obtained by SEM (Zeiss sigma500) at an accelerating voltage of 10 kV. The pore structure distribution (PSD) of zeolite 3A was measured by the low-temperature CO 2 adsorption method at 273 K using 3H-2000 PM.

Adsorption Isotherm and Kinetics
Measurements. Equilibrium and kinetics measurements of HF adsorption on zeolite 3A were performed in a static volumetric instrument. The detailed composition of the experimental instrument is shown in Figure 1. The pressures of the experimental system were recorded using a precision pressure sensor (Keller, Switzerland, ±0.01% precision). An automatic thermostatic oil bath with ±0.01 K accuracy was used to keep the system temperature constant.
The main experimental procedures for equilibria and kinetics adsorption are: (1) Sample pretreatment and filling. The 3A sample used in the adsorption experiment was first dried in the drying oven at 385 K for 24 h, weighed, and then quickly transferred to the sample cell. The constant temperature of the oil bath was set to 358 K and vacuumed for 12 h to remove the water and other impurities adsorbed in the sample transfer process (2) Leakage test for adsorption system. After the sample pretreatment, the temperature of the oil bath was set to 338 K. After the system temperature is constant, helium (He) was fed into the reference and adsorption cells to a certain value (maximum pressure required for adsorption). The leakage testing time was 12 h. The system was heated and vacuumed for 4 h after the leakage test (3) Measurement of void volume. The void volume of the adsorption cell filled with adsorbent was measured by the helium expansion method. The temperature for the helium expansion test was set to the adsorption experimental temperature. The helium density used in void volume measurement is calculated by the BWR equation [21]. Helium expansion tests were repeated at least four times and the values were averaged to minimize error (4) Equilibrium adsorption experiment. A certain HF amount was injected into the reference cell when the temperature of the adsorption system was constant with the experimental setting value. The gas pressure for the reference cell was recorded until it is stable. The valve between the reference and sample cells was opened, the HF flowed into the sample cell from the reference cell and the adsorption began.
The system pressure was recorded when the system pressure was constant for 5 h. The above operations were repeated until the maximum pressure point of the equilibrium adsorption is reached 2 Adsorption Science & Technology (5) Adsorption kinetics experiment. HF was injected into the reference cell to a calculated pressure value when the temperature of the adsorption system was constant with the experimental set value. The valve between the reference and sample cells was then opened, the HF flowed into the adsorption cell, and the system pressure was recorded in real time by the high-precision pressure transmitter until the final system pressure was stabilized at the desired value. The kinetic experimental conditions were 0.15, 0.35, 0.5, 0.70, 0.95, and 1.10 bar at 298 K and 298, 308, 318, 328, and 338 K at 0.5 bar During the experiment, the calculation process for HF adsorption capacity was based on the mass conservation under the ideal gas state, so the adsorption capacity is [19]: where q is the adsorption capacity, p 0 is the gas pressure in the reference cell before adsorption and p 1 is the system gas pressure during the adsorption process. R and T are the ideal gas constant and the system temperature, respectively.  The relationship between adsorption capacity, temperature, and pressure in the adsorption process is usually described by the adsorption model. The Langmuir and Freundlich models had been successful in adapting the described HF adsorption process [11,12]. The DA model is often used to describe the gas adsorption in micropores [22]. The DA model is expressed as follows: where q and q m are the adsorption capacity and saturated adsorption capacity of gas, p 0 is the saturated vapor pressure of gas, and E is the adsorption characteristic energy.
The saturated vapor pressure p 0 of HF is calculated by the Antoine equation with extension [23]: where the values of A, B, C, D, and E are 23.7347, −1799.6, −6.1764, 5.0046 × 10 −10 , and 6.15 × 10 −6 , respectively. The unit of p 0 is mmHg in the calculation formula, which needs to be multiplied by 0.0013332 to convert to bar.   Adsorption Science & Technology than the unadsorbed state for adsorbate molecule [24]. The Clausius-Clapeyron equation was used to calculate the isosteric heat at constant adsorption capacity, and its expression is as follows [24]: where H st is the calculated isosteric heat at the adsorption capacity q.
The isosteric heat of adsorption (H st ) can be obtained by combining the DA model with the Clausius-Clapeyron equation [24]: where α is the thermal expansion coefficient of adsorbed phase and α is calculated by 1/T. h f g is the enthalpy of evaporation, which is expressed as follows:

Adsorption Kinetics
(1) Adsorption Rate Models. The pseudo first-order (PFO) model assumes that the adsorption rate is directly proportional to the number of active sites. The rate expression is given by [25,26]: where q 1 and q represent the adsorption capacity at equilibrium and time t, respectively, and k 1 is the first-order rate constant.
After integrating Eq. (7) with the boundary conditions (t =0, q =0 and t = t ∞ , q = q 1 ), The pseudo second-order model (PSO) assumes the adsorption rate is proportional to the square of number of active sites; therefore, the rate model is given by [25,26]: where q 2 represents the equilibrium adsorption capacity, and k 2 is the second-order rate constant. With the boundary conditions (t =0, q =0 and t = t ∞ , q = q 2 ), the PSO is given by [25,26]: (2) Diffusion Models.
(1) Film diffusion Boyd's film diffusion model assumes that the main diffusion resistance is boundary film around the adsorbent particles [26,27]. Therefore, the model was used to distinguish whether the film diffusion was a rate-controlling step in HF adsorption on 3A. The expression of this model is given by [26]: where F is the fraction at times t and can be expressed using adsorption capacity: B t can be expressed in F [28]: when F > 0:85, (2) Interparticle diffusion The model assumes interparticle diffusion control the adsorption rate, the model is given by the following expression [26,27]: where q t is adsorption amount at time t, q e is adsorption amount at infinite time, D c is interparticle diffusion coefficient, and r p is the radius of the particle. When q t /q e value is >70%, Eq. (15) can be simplified to [26]: 3400 3600 3800 4000 4200 4400 4600 4800 5000 5200   By plotting ln ½1 − q t /q e versus t, the linear slope is − π 2 ðD c /r 2 p Þ and the intercept is ln6/π 2 .
(3) Intraparticle diffusion If the HF adsorption process on zeolite 3A is controlled by intraparticle diffusion, the HF adsorption (q t ) should change linearly with time square root (t 1/2 ) [26,27]. The intraparticle diffusion model assumes three steps typically in adsorption: bulk or film diffusion, gradual adsorption, and the final adsorption equilibrium [26]. The intraparticle diffusion model can be used for determining the adsorption mechanism, and the model is given by [29]: where k i is the intraparticle diffusion constant and C is the intercept. The intercept value is related to the film

Results and Discussion
3.1. Characterizations of the Zeolite 3A. The powder XRD pattern of zeolite 3A is shown in Figure 2. The XRD characteristic diffraction peaks are in accordance with 3A crystals by Wang [30], indicating that the sample is a typical 3A zeolite. The FT-IR of zeolite 3A spectra is shown in Figure 3. The broad bands around 3,447.9 cm −1 belong to the stretching vibration of the O-H group. The bands at around 1,659.3 cm −1 belong to the bending and tensile vibrations of the O-H and Si-OH groups, respectively [31]. At about 1,005.5, 665.8, and 462.9 cm −1 , three typical absorption peaks formed by amorphous aluminosilicate structure were observed, which belong to the asymmetric tensile vibration of the tetrahedron in the sample, the symmetric tensile vibration of the Si-O-Si bond, and the bending vibration of T-O (T is Si or Al) bond, respectively [31]. The same absorption peaks have also been found by Cao et al. [31]. The absorption peak at 554.6 cm −1 belongs to the vibration of double binary ring D4R, the main subunit structure unit of zeolite 3A.
The surface morphology of zeolite 3A is shown in Figure 4, which indicates that the surfaces of the samples are not homogeneous. Zeolite 3A has a large grain size, good dispersibility, and many surface fractures, which has also been found in the literature [30]. The low-temperature CO 2 isotherm and PSD of zeolite 3A are shown in Figure 5. The pore size of zeolite 3A is mainly concentrated at~0.3 and 0.37 nm. The specific surface area of zeolite 3A and the pore volume are 402.81 m 2 /g and 0.1325 cm 3 /g, respectively.

Adsorption Equilibrium
3.2.1. Adsorption Isotherms. The HF adsorption isotherms on zeolite 3A at 298 to 338 K are shown in Figure 6. All adsorption isotherms are type I according to the IUPAC classification. The adsorption amount of HF on zeolite 3A is close to that on activated carbon by Shen et al. [15] but less than that on activated carbon by Bahrami et al. [19] and MOF [18], which is mainly due to the pore characteristics and surface physicochemical properties of the adsorbent. The adsorption capacity increases with increasing pressure at the same temperature. The adsorption capacity decreases with the temperature increase because the binding strength of the zeolite 3A surface decreases with temperature, indicating that the HF adsorption on the zeolite 3A surface is exothermic. The Langmuir, Freundlich, and DA models were used to correlate the HF adsorption data using regression analysis [19]. Three statistical functions, including SSE, RMSE, and χ 2 , were used to estimate the goodness of fit [24].
The fitting parameter values for the three models are listed in Table 1. The Langmuir parameters q 0 and b L all decreased with the increasing temperature, which implied that the bonding of HF to zeolite decreased. The Freundlich parameter k F represents the adsorption capacity, and n shows the adsorption driving force and indicates the favourability of adsorption [10]. The k F values showed a negative correlation with temperature. The n values (1.271-1.368), determined from the Freundlich model suggested the moderately difficult nature of HF adsorption on zeolite 3A [10]. The same results for the Freundlich parameter n have been found in HF adsorption on produced NaF pellets and activated carbon [10][11][12].
The values of SSE, RMSE, and χ 2 in Table 1 for the DA model are the smallest among the three adsorption models, 9 Adsorption Science & Technology which indicates that the DA model predicted the adsorption data better than the other models. Figure 6 also shows the predicted HF adsorption data by the DA model. The DA model can describe well the HF adsorption curve, indicating that micropore filling existed in HF adsorption on zeolite 3A. Both characteristic energy E and DA parameter n increased with the increase of temperature. The DA parameter n can be used to describe the surface heterogeneity of the adsorbent [13,24]. The study shows that the closer the n value is to 1, the higher the heterogeneity of the zeolite 3A surface, which is more related to the active adsorption site at low temperatures. The n value shows a positive correlation with the characteristics energy value E (Figure 7). The lower n value results in E value.
The regeneration of 3A zeolite affects the reuse of materials [18]. The adsorption capacity after regeneration was tested at 298 K. The zeolite 3A after HF adsorption was treated at 393 K in the muffle furnace for 12 h and then quickly transferred and sealed to the adsorption equipment (sample cell in Figure 5) and then vacuumed at 373 K for 10 h to reactivate its capacity. The HF repeatability experiments on zeolite 3A at 298 K are shown in Figure 8(a). The remaining HF capacity after regeneration was slightly reduced, which suggested the high stability of exposure to dry HF gas. The IR spectra of the zeolite 3A before and after HF adsorption are shown in Figure 8(b), suggesting no change in the surface functional group of the zeolite 3A molecular sieve after HF adsorption.

Heat of Adsorption.
The isosteric heat of adsorption (H st ) calculated using Eq. (5) is plotted with the HF surface load (q/q m ), as shown in Figure 9. H st values of HF on zeolite 3A decrease with the increase of surface loading, which indicates the surface energy heterogeneity of zeolite 3A. A large number of vacant adsorption sites exist on the zeolite 3A surface at the initial stage, which results in strong interaction between the HF molecules and zeolite 3A. Therefore, the H st value is high at the initial adsorption stage. With further adsorption capacity increase, the adsorption sites are occupied, resulting in the weakening of the interaction between the HF molecules and zeolite 3A. Therefore, H st decreases with the surface loading.

Effect of Pressure on Adsorption
Kinetic. The adsorption kinetic curve of HF on zeolite 3A at 298 K is shown in Figure 10. The adsorption capacity continuously increases with the adsorption time. The adsorption capacity rapidly increases in the initial stage, which is mainly due to the high adsorption energy sites on the adsorbent surface and the rapid occupation of HF molecules onto the adsorption sites, so the adsorption rate is high. When large numbers of active adsorption sites are occupied, the adsorption sites available for HF molecules become less; hence, the adsorption rate becomes slow. When the active sites are almost occupied,  10 Adsorption Science & Technology HF molecules reach dynamic equilibrium on the surface of zeolite 3A. At this time, the number of HF molecules on the adsorbent surface will not increase; therefore, the total adsorption rate tends to zero and the adsorption capacity will not increase.
(1) Adsorption Rate. The HF adsorption kinetics curves and the corresponding curves predicted by the PFO model are shown in Figure 10. The kinetic constants calculated with the PFO and PSO models and the characteristic parameters of the kinetic model are listed in Table 2. The larger R 2 of the  Figure 11: Adsorption mechanism of HF on 3A zeolite. PFO model indicates a better predictivity of the adsorption kinetic process over the whole pressure range well, so the adsorption of HF by zeolite 3A follows simple first-order kinetics. The calculated q 1 values from the PFO model are more in accordance with the experimental data than the calculated q 2 values by the PSO model.
Both q 1 and q 2 increase with pressure, which is consistent with the high pressure is conducive to HF adsorption. Meanwhile, the PFO and PSO model parameters (k 1 and k 2 ) increase with pressure because the velocity of the HF molecule will increase and the migration speed of HF molecules inside the pores will be faster at higher pressure, leading to an increase in the diffusion rate. The increased mass transfer coefficient makes the kinetic curve steeper. The k 1 values (0.0395-0.0743 min −1 ) for HF on zeolite 3A are larger than that of H 2 S (0.00177min −1 ) on zeolite 3A, which suggested that the HF adsorption rate is faster than that of H 2 S [6,32].
(2) Diffusion Mechanism. The HF adsorption process on zeolite 3A can be divided into five continuous processes, as shown in Figure 11: (a) HF molecules diffuse from free state to adsorbent gas film (bulk diffusion), (b) HF molecules diffuse in the film (film diffusion), (c) HF molecules diffuse in pores between particles (interparticle diffusion), (d) HF molecules diffuse in 3A crystal (intraparticle diffusion), and (e) surface adsorption. Generally, surface adsorption and bulk diffusion are very fast, so they are ignored in the analysis of the HF adsorption kinetic process, and the role of the other three processes in the adsorption process is mainly considered.
(1) Film diffusion Equations (11)- (14) can be used to predict HF adsorption control by particle or film diffusion [26,27]. By plotting B t vs. time, the mass transfer rate is controlled by pore diffusion if the plot is a straight line passing through the origin [26][27][28]. If the curve is nonlinear or does not pass through the origin, film diffusion can also be considered to affect the adsorption rate [26][27][28]. The calculated B t values as a function of time are shown in Figure 12, and the curves    Table 3. The existence of intercept values suggests that the ratecontrolling process for HF adsorption on zeolite 3A is mainly controlled by external mass transfer [26]. The intercept increases close to zero with increasing pressure, which indicates that pore diffusion appears to be more important at slightly higher pressures.
(2) Interparticle diffusion The fitting results for the adsorption kinetics data with the interparticle diffusion linear form are shown in Figure 13 and Table 4. A good correlation between ln ½1 − q t /q e with t indicates that the experimental data can be well described by the model. The D c /r 2 p and D c values calculated 13 Adsorption Science & Technology by the slope of the linear fit increase continuously with increasing pressure, indicating that high pressure favours the HF gas diffusion in zeolite 3A, which is consistent with the previous calculation that high pressure favours the adsorption rate. The order of magnitude of interparticle diffusion coefficient D c of HF is 10 −9 m 2 /s. The intercept values obtained by linear fitting the experimental data are not close to −0.4977 (ln6/π 2 ), which indicates that the HF adsorption rate on zeolite 3A is not completely controlled by the interparticle diffusion. The calculated interparticle diffusion coefficient D c /r 2 p versus first-order rate constant k 1 does not satisfy the 1/15 relationship and also reflects that the HF diffusion over zeolite 3A is not completely controlled by interparticle diffusion [26].
(3) Intraparticle diffusion The intraparticle diffusion curves for HF adsorption on zeolite 3A are shown in Figure 14. The curves are observed to be nonlinear. According to the intraparticle diffusion model, the adsorption amount q t plot vs. t 1/2 should be a straight line and pass through the origin of the mass transfer rate solely controlled by intraparticle diffusion [26]. Multilinearity can be observed when different processes are related to the adsorption mechanism. The slope of each linear portion represents the corresponding adsorption rate [27]. The curves in Figure 14 show quart-linear characteristics, indicating the presence of four consecutive HF adsorption steps on zeolite 3A. The initial part, region I, was ascribed to the HF diffusion from the free state to and through the gas film of zeolite 3A, and the initial stage was associated with the HF film diffusion. The second and third parts, regions II and III describe the adsorption-progressive stage where diffusion occurs within  Figure 16: Plots of Boyd's film model for HF adsorption on 3A.  Figure 14 shows that the extended line for straight lines corresponding to regions II and III does not pass through the origin, which may result from the different diffusion rates for the first and other stages of HF adsorption. The adsorption rate of each HF molecule is controlled by four processes at a certain pressure. In the first stage, the adsorption rate is controlled by the film diffusion, as shown in segment I of Figure 14. The interparticle diffusion would be the rate-controlling process when the HF molecules diffuse into the particle pore. The intraparticle diffusion then controls the rate. Finally, the rate is governed by surface adsorption. The slope of each linear segment before region IV can be used to deduce which stage controls the whole HF adsorption rate [27]. Moreover, Figure 14 shows that the linear slope in region III is the smallest. The linear portion with the lowest slope is the rate-controlling step [27,28]. This suggests that HF diffusion into the 3A particle is slower and the intraparticle diffusion stage is the ratecontrolling stage. The linear slopes for regions II and III gradually increase with increasing pressure, which may be because the film thickness decreases as the pressure rises and the film diffusion rate will increase.

Effect of Temperature.
Kinetic experiments at 298-338 K and 0.50 bar were conducted to investigate the effect of temperature on the HF kinetic adsorption on zeolites 3A. The HF adsorption kinetics curves at five temperatures are shown in Figure 15. The HF kinetic curves at five temperatures can be divided into three different stages: rapid adsorption, slow adsorption, and adsorption equilibrium stages. The HF adsorption capacity gradually decreases with increasing temperature, which is consistent with the discovery in literature [19]. High temperature enables the HF molecules to obtain more energy to overcome the intermolecular force and electrostatic interaction; thus, HF molecules remain in the free state. Simultaneously, the affinity between HF molecules and zeolite 3A surface will be weakened by the increasing temperature and will result in the decrease in adsorption capacity. In addition, high temperature is not conducive to the formation of a more orderly and stable configuration of adsorbed molecules, resulting in the reduction of adsorption capacity. Higher temperature is not conducive to improving HF adsorption recovery. It can also be found that the time to reach equilibrium decreases with the increasing temperature, which also shows that the temperature is conducive to the gas adsorption rate.
(1) Adsorption Rate. The PFO and PSO models were used to predict the HF adsorption kinetic curves. The fitting results are shown in Figure 15 and Table 5. The adsorption rate constants k 1 and k 2 of HF increase with temperature,   indicating that the HF adsorption rate shows a positive correlation with temperature, which is due to the acceleration of molecular motion caused by temperature increase. q 1 and q 2 decrease with increasing temperature, suggesting that higher temperature is conducive to gas diffusion; however, it will reduce the gas adsorption capacity.
(1) Film diffusion The calculated B t values are shown in Figure 16. Moreover, B t is not completely linear. The slope and intercept of Boyd's film diffusion model are listed in Table 6. The existence of intercept value shows that the HF adsorption rate on zeolite 3A was mainly controlled by external mass transfer. The intercept values decrease slightly with increasing temperature, suggesting that pore diffusion become more important at high temperature. The B t values show a positive correlation with temperature, indicating that the film diffusion rate increases with temperature increase.
(2) Interparticle diffusion The linear fitting results for the HF adsorption kinetic curve using the interparticle diffusion model at five temperatures and 0.5 bar are shown in Figure 17 and Table 7. A good correlation between ln ½1 − q t /q e shows that t indicates that the experimental data can be well described by the model. In addition, Table 7 shows that D c /r 2 p and D c values increase with increasing temperature. Based on molecular dynamics theory, the essence of HF diffusion in zeolite 3A is due to the molecular irregular thermal movement. The average kinetic energy of HF molecules increases with temperature, which leads to an increase in amplitude and frequency of molecular vibration. The higher molecular irregular thermal movement intensity occurs at a higher temperature. Consequently, the HF diffusion rate in zeolite 3A becomes faster.

(3) Intraparticle diffusion
The quart-linear characteristics for intraparticle diffusion curves can be observed in Figure 18, which indicated that three different processes are involved in the adsorption mechanism. The curves of the five temperatures in the first stage I are relatively close, indicating that the influence of temperature on the HF diffusion velocity from the free state to and through the gas film of zeolite 3A is little. In the second and third stages, the difference between the five curves increases with time. The linear slope for each stage decreases with increasing temperature, suggesting that the pore diffusion rate in the particles decreases slightly at this stage, which may be due to the impediment of diffusion by the interaction of HF molecules under high pressure. The same trend is also found in CO 2 in porous MgO [27]. The extended line for straight lines corresponding to regions II and III does not pass through the origin. The linear portion for stage III with the lowest slope indicates that the diffusion rate is controlled by intraparticle diffusion.

Conclusion
HF adsorption isotherms data on zeolite 3A at 298-338 K were type I. The adsorption capacity increased with  The HF adsorption by zeolite 3A at different adsorption conditions followed simple first-order kinetics. The firstorder rate increased with pressure and temperature. The calculated B t values showed nonlinear characteristics with time. B t values increased with increasing pressure and temperature, indicating that high pressure and temperature favour film diffusion. The order of magnitude of the interparticle diffusion coefficient D c of HF on zeolite 3A was 10 −9 m 2 /s. The higher pressure and higher temperature favour the HF diffusion process. The intraparticle diffusion curves for HF adsorption on zeolite 3A showed quart-linear characteristics, indicating the presence of three consecutive HF adsorption steps on zeolite 3A. In the initial adsorption stage, the film diffusion controlled the adsorption rate. The intraparticle diffusion would be the rate-controlling process after HF molecules diffuse through the film. The intraparticle diffusion then controls the rate. Finally, the rate was governed by surface adsorption. HF diffusion into the particle of zeolite 3A was the slowest and therefore, the intraparticle diffusion is the rate-controlling stage.
In this study, the HF adsorption mechanism onto zeolite 3A was evaluated based on kinetics and thermodynamics. However, only one adsorbent zeolite 3A was studied due to limited conditions, and other porous materials have not been analyzed in this paper. The adsorption kinetics and thermodynamics of HF on different materials will be compared in subsequent studies.

Data Availability
All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest
The authors declare that they have no conflicts of interest.