3.1. Sorbents Characterization
The pure Li
4SiO
4 powders obtained in each run of the experimental design (both for the preliminary screening and the RSM design) were analyzed by XRD to verify the formation of lithium orthosilicate phase. In
Figure 1 are reported, as an example, the diffractograms of the sorbents obtained by the RSM experimentations calcined for 10 h at different temperatures (600-900 °C). As shown, the sorbents exhibit a pattern with major peaks of the Li
4SiO
4 crystalline phase, confirming that the selected temperatures are sufficient for assessing the sorbent synthesis. The sorbent calcined at 600 °C presents small peaks attributable to the reagent lithium carbonate. The presence of Li
2CO
3 is probably attributed to an incomplete conversion of SiO
2 and LiOH due to the low synthesis temperature; LiOH is then carbonated owed to contact with air during calcining [
27,
40].
The Li
4SiO
4 powders were also characterized by SEM to assess the sorbents morphology and particle dimensions.
Figure 2 displays SEM images of the pure sorbents calcined for 10 h at different temperatures (relative to sample 10, 11, and 19 of the experimental runs of
Table 2). The sorbent calcined at 600 °C (
Figure 2a,d) shows dense and non-porous particles with sizes of below 5 µm agglomerated in clusters of very large dimensions (greater than 100 µm). By increasing the calcining temperature to 750 and 900 °C, the Li
4SiO
4 sorbents present a significant sintering of the particles, showing essentially a non-porous particle morphology and greater diameters of about 30-50 µm. Besides,
Figure 3 reports images of the morphology of the doped sorbents calcined at different temperatures obtained by SEM analysis. As example, the sorbent with 10 vol% of potassium carbonate was analyzed. All the doped sorbents show the same morphology and particle size evidenced for the pure Li
4SiO
4, demonstrating that the solid-state doping method does not affect the overall structure of the sorbent. K
2CO
3 is constituted by porous particles that results well distributed among the dense Li
4SiO
4 particles, as can be seen in
Figure 3b,c.
3.2. Preliminary Variables Selection for Parametric Modeling
To fully comprehend the adsorption process of doped Li
4SiO
4 sorbents and how it is influenced by the sorbent’s synthesis and process conditions, a DoE methodology was applied. Initially, five variables were selected as potentially significant for the adsorption process. The variables concerned both the sorbent preparation conditions and the adsorption operating conditions: calcining temperature and time, K
2CO
3 content, adsorption temperature, and CO
2 concentration. For the synthesis temperature, the low level was set to 600 °C since it represents the minimum temperature necessary for the Li
4SiO
4 synthesis from the precursors [
41], while the high level (900 °C) is the typical synthesis temperature for a solid-state method. A higher value could lead to Li
4SiO
4 decomposition to Li
2SiO
3 by lithium sublimation [
42]. The calcining time was varied from 4 h, which is a typical synthesis time used in literature for a complete conversion of reagents to Li
4SiO
4 [
43], to 10 h. A longer synthesis time was not considered since, at 900 °C, could promote lithium sublimation [
42]. The K
2CO
3 content was chosen from a previous study on Li
4SiO
4-based sorbents [
44]. A promoter amount exceeding 40 wt% has not been considered appropriate since reduces the amount of active sorbent.
The CO
2 concentration in vol% refers to the volume fraction of CO
2 in the TGA feed gas, consisting of a CO
2/N
2 mixture flow. The low level was set to 4 vol% (corresponding to a CO
2 partial pressure of 0.04 atm), which is typical of exhaust hot gases from gas turbines [
45,
46,
47,
48], while the high level (50 vol%) depended on the thermogravimetric analyzer operating limit. Then, the adsorption temperature range (500-700 °C) was chosen considering the curve of the equilibrium partial pressure of CO
2 at different temperatures, which has been evaluated on the basis of the Gibbs free energy changes of the adsorption/desorption reaction (Equation (1)) [
49]. The heterogeneous reaction that occurs between Li
4SiO
4 and CO
2 is an equilibrium reaction; therefore, for each CO
2 partial pressure there is a thermodynamic equilibrium temperature below which the sorption reaction could proceed (ΔG < 0), otherwise the desorption process takes place (ΔG > 0). At equilibrium of Equation (1), the following equation can be valid since ΔG = 0 is satisfied:
where
,
and
represent the variation of standard Gibbs energy, standard enthalpy and standard entropy of the carbonation reaction, respectively, R is the ideal gas constant, and T is the temperature.
Keq is the equilibrium constant for the reaction (Equation (1)) and is defined as in Equation (6):
where
pCO2 is the CO
2 partial pressure in the treated gas.
The
pCO2–temperature equilibrium curve can be evaluated from Equation (5) and Equation (6), as found in a previous work of the authors [
49], and it is reported in
Figure 4. As shown, the two-dimensional design space (for the two factors of adsorption temperature and CO
2 concentration) covers the entire region where the adsorption reaction is activated both at the low and high CO
2 concentration levels (4 and 50 vol%).
Preliminary calcining and adsorption tests were conducted with experimental tests using TGA to assess the effect of calcining time on the synthesis and the adsorption performance of the sorbent. Calcining in TGA was performed by placing the mixed reagent powders in the platinum crucible and heating them to the selected synthesis temperature (600 or 900 °C) under air flow. The temperature was maintained for several hours thus to observe the weight loss during the synthesis reaction. The results are reported in
Figure 5.
As shown, the main difference among the two thermograms (
Figure 5a,b) is related to the synthesis temperature. A weight loss is observed in the range between 400-600 °C and this is attributed to the LiOH decomposition to lithium oxide, Li
2O [
40]. Therefore, when the Li
4SiO
4 synthesis takes place at 600 °C, the time needed to reach the complete LiOH decomposition and, supposedly, the complete conversion of reagents is higher than for the synthesis at 900 °C, where the sample weight becomes stable even before reaching the isotherm section. Whereas, for the same synthesis temperature, the synthesis time (4 or 10 h of isotherm) does not affect markedly the conversion of reagents. A variation in the weight loss less than 0.5 % is observed from 4 to 10 h.
Moreover, preliminary adsorption tests were carried out on 10 wt% K
2CO
3 doped-Li
4SiO
4 sorbents prepared with different synthesis times and temperatures. The results of the adsorption experiments in terms of both adsorption capacity and Li
4SiO
4 conversion are reported in
Figure 6, considering an adsorption temperature of 500 °C and a CO
2 concentration of 4 vol%.
As illustrated, the adsorption performance is almost unchanged by varying the calcining time from 4 to 10 h. The main difference in the CO2 uptake could be ascribed to the calcining temperature effect. An increase in the calcining temperature from 600 to 900 °C leads to a decrease in the adsorption capacity and sorbent conversion, which change from about 240 to 200 mg/g and 73 to 62 %, respectively. Further considerations on the calcining temperature effect will be dealt with in the modeling section. These experimental considerations evidence that the synthesis time is poorly significant in affecting the adsorption performance of the sorbent. Therefore, it was not considered in the following RSM analysis and set to a value that allow to obtain a complete conversion of the reactants into Li4SiO4 for all the synthesis temperatures of the design space. Thus, a synthesis time of 10 h was selected and the five variables were consequently reduced to four.
3.3. Modeling of the CO2 Adsorption Process and Statistical Analysis
A Response Surface Method DoE was conducted to investigate the combined effect of the variables on the adsorption process and to develop non-linear models of the selected responses Li
4SiO
4 conversion and sorbent’s adsorption capacity. The obtained experimental design matrix along with the responses results are reported in
Table 2. By using the software Design Expert 11, regression calculations were executed to find the polynomial model that best fits the experimental data for each response. The analysis of variance (ANOVA) was conducted to confirm the statistical significance of all terms in the regression models, by means of statistical parameters. Each term of the resulting models corresponds to a main or interaction effect of a variable (or combination of variables) on the response. ANOVA splits data variability into two sources of variation, the model and the experimental error, and it uses F-values and p-values calculated for each term to verify their significance. Significant model terms were considered the ones having a p-value < 0.05. The goodness of models also depends on the determination coefficient R
2 and on its modification, adjusted R
2, which only increases when terms really affecting the model are selected. Predicted R
2 was also checked as a measure of the predictivity of the model.
ANOVA results and modeling statistics related to the RSM design for both responses are shown in
Table 3 and
Table 4. Both for Li
4SiO
4 conversion and adsorption capacity a quadratic model with additional mid-cubic terms was found to best fit the experimental data. The ANOVA results show that the models are both statistically significant (p-value < 0.05). Moreover, not significant lack of fit denotes the reliability of their predictive quality (p-value > 0.1). An estimation of pure error is also reported, evaluated through the replicates on the central points. The sum of squares of pure error is found to be at least two orders of magnitude lower than that of the models, representing their goodness in fitting the experimental data. High values of the determination coefficients R
2 (0.9869 and 0.9710 for the Li
4SiO
4 conversion and adsorption capacity models, respectively) and adjusted R
2 > 0.85 also denote the goodness of the models selected.
As shown in
Table 3, the effects that result highly significant for the Li
4SiO
4 conversion (high values of sum of squares and p-value < 0.01) are the main effect of all factors, their quadratic terms and the two-factor interactions AC (synthesis temperature-adsorption temperature) and CD (adsorption temperature-CO
2 concentration). In particular, the main effect of factor A - synthesis temperature has a sum of square two times greater than other main effects, indicating that the synthesis temperature can affect markedly the sorbent morphology (as can be seen by SEM analysis reported in
Figure 2) and thus the conversion. Moreover, the interaction CD is found to be the most significant effect showing the highest sum of squares among all the effects. This is related to the thermodynamic of the adsorption reaction on Li
4SiO
4, for which at 700 °C the adsorption performance is high for 50 vol% of CO
2, whereas is very low at 4 vol% due to the activation of desorption process (
Figure 4). Likewise, for the adsorption capacity (
Table 4) the main effect of A - synthesis temperature and the interactions AC (synthesis temperature-adsorption temperature) and CD (adsorption temperature-CO
2 concentration) result the most significant effects, since they have the highest sum of squares. Whereas the main effect of factor D – CO
2 concentration is almost not significant (p-value ≈ 0.05), and the factor B – K
2CO
3 content results not significant. However, the factor B is considered in the model to support hierarchy, since the interaction effect BD (K
2CO
3 content-CO
2 concentration) results significant. BD significance indicates that the addition of K
2CO
3 to Li
4SiO
4 affects the adsorption capacity mostly in relation to the CO
2 content in the gas flow, and this was also noted in our previous study [
27]. The addition of K
2CO
3 to Li
4SiO
4 improves more markedly the CO
2 uptake at low CO
2 concentrations (4 vol%) than at higher concentrations (50 vol%), since a high
pCO2 seems to compensate for the lower diffusivity of the CO
2 in the product layer.
The adequacy of the obtained models was verified by diagnostic plots. For each response, the accordance between experimental data and predicted values is confirmed by the close distribution of run data to the line of the correspondent plot (
Figure 7). The relationship between the factors, single or in interaction, can be visualized through the responses model equations, which are reported in terms of coded factors in Equation (7) and (8).
Three-dimensional response surface plots of the models for both responses are reported in
Figure 8 and
Figure 9, and the relationships between factors and responses are discussed in the sections below. The 3D model graphs displayed the trend of both responses in the design space as a function of two significant variables, setting the other two variables at the center or factorial values.
3.3.1. Effect of Synthesis Temperature
According to the ANOVA results presented in
Table 3 and
Table 4, synthesis temperature is an important variable that affects both the Li
4SiO
4 conversion and the adsorption capacity, and it demonstrate a quadratic effect on the responses (due to A
2 term) visible by a curvature in the model graphs of
Figure 8a,b and
Figure 9a,b. However, its effect on the responses cannot be considered individually as this factor presents interactions with adsorption temperature (AC, A
2C) and CO
2 concentration (AD, A
2D).
Figure 8a and
Figure 9a displays higher values of
XLi4SiO4 and adsorption capacity when the synthesis temperature is 600 °C, reaching about 92 % and 270 mg/g, respectively, when the K
2CO
3 content is 25 wt%, the CO
2 in the treated gas is 27 vol% and the adsorption temperature is 600 °C. Moreover, these values increase as CO
2 concentration increase from 4 to 50 vol% (
Figure 8b and
Figure 9b), raising from near 40 to 100 % and 100 to 280 mg/g for
XLi4SiO4 and adsorption capacity, respectively (maintaining 25 wt% of K
2CO
3 and 600 °C of sorption temperature). These results are related both to sorbent morphology and adsorption reaction kinetics. As shown by SEM analysis reported in
Figure 2, a lower synthesis temperature leads to a less sintered sorbent with smaller particle size and thus to higher adsorption capacities and conversions [
40,
50,
51,
52]. Moreover, when the CO
2 concentration increase from 4 to 50 vol%, the CO
2 concentration gradient between the bulk and the particle surface increase, so that increasing the amount of CO
2 adsorbed [
17,
26,
53].
3.3.2. Effect of K2CO3 Content
As evident by ANOVA (
Table 3), the variable K
2CO
3 content affects Li
4SiO
4 conversion, and it demonstrate also a quadratic effect on the response (due to B
2 term) visible by a slight curvature in the model graph of
Figure 8c. On the contrary, for adsorption capacity, the factor K
2CO
3 content results not significant (
Table 4), and this is shown by a linear contour plot in
Figure 9c. This different result can be explained by considering that an increase in K
2CO
3 content helps the reaction kinetics by reducing diffusion resistance through the product layer that covers the sorbent particles, leading to higher Li
4SiO
4 conversion. Nevertheless, higher values of K
2CO
3 content reduce the amount of active sorbent (Li
4SiO
4, which actually reacts with CO
2), leading to a lower adsorption capacity. However, the interaction effect BD (K
2CO
3 content-CO
2 concentration) results significant, and this is represented on 3D response surface plots by a slight curvature of the surface when CO
2 concentration varies from 4 to 50 vol%. The addition of K
2CO
3 to Li
4SiO
4 improves more markedly the adsorption performance at low CO
2 concentrations, where the diffusion resistance is high due to lower gradient for the CO
2 mass transfer into the product layer. Li
4SiO
4 conversion and adsorption capacity increase from 40 to 55 % and from 135 to 160 mg/g, respectively, when the K
2CO
3 content increased from 10 to 40 wt % for 4 vol% of CO
2 in the treated gas (considering synthesis and adsorption temperature of 600 °C). Whereas at higher
pCO2, the increasing concentration gradient compensates for the lower diffusivity of the CO
2 in the product layer [
27,
30]. When the CO
2 concentration is 50 vol% and K
2CO
3 content increased from 10 to 40 wt %, Li
4SiO
4 conversion ranged 90-95 %, while adsorption capacity remains almost constant to 260 mg/g (synthesis and adsorption temperature set to 600 °C).
3.3.3. Effect of Adsorption Temperature and CO2 Concentration
Based on the ANOVA results (
Table 3 and
Table 4) and on the model’s equations (Equations (7) and (8)), the term mostly affecting either Li
4SiO
4 conversion and adsorption capacity is the two-factor interaction CD (adsorption temperature-CO
2 concentration). It shows the highest sum of squares and a positive effect on the responses (the coefficient in the model’s equations is the highest in term of coded factors and has a positive sign). The positive effect can be noted in response surface graphs of
Figure 7d and
Figure 8d, where both Li
4SiO
4 conversion and adsorption capacity increase with increasing the adsorption temperature and CO
2 concentration. This term is strictly correlated to the adsorption reaction thermodynamic, since the first condition for obtaining high CO
2 removal performance is to work in a region where adsorption is favored [
17,
26,
53]. It is important to highlight that, at 4 vol% CO
2, both the Li
4SiO
4 conversion and adsorption capacity show a decrease by increasing the adsorption temperature from 500 to 700 °C (
Figure 7d and
Figure 8d). This is associated to the adsorption reaction equilibrium, for which at 500 °C the adsorption process is favored, whereas at 700 °C the desorption process is activated (
Figure 4) [
17,
27,
54].
3.4. Optimization and Validation Tests
The process of CO
2 capture using K
2CO
3-doped Li
4SiO
4 sorbents requires high CO
2 adsorption capacities to be efficient for practical applications. Moreover, high Li
4SiO
4 conversions are necessary for the process to be economically feasible. Therefore, the parametric models obtained by RSM analysis were used for a multi-objective optimization that allow to predict the optimal combination of variables for maximizing both the Li
4SiO
4 conversion and the adsorption capacity. As evidenced by response surface plots, by varying input parameters the two responses do not behave the same way. As previously observed, K
2CO
3 kinetically improves the sorbent performance in terms of conversion but, on the other hand, it represents an inert that reduces the amount of active sorbent (Li
4SiO
4), with a consequent reduction of adsorption capacity for high K
2CO
3 content. The optimization was carried out for three different CO
2 concentrations (4, 27, 50 vol%), letting the other factors (synthesis temperature, K
2CO
3 content and adsorption temperature) vary in their range and maximizing the two responses. The optimum conditions for the process variables and the predicted values for the responses, obtained for the three different CO
2 concentrations by numerical optimization, are reported in
Table 5 and
Table 6. For each different optimization (at different CO
2 concentration), the models predict an optimum working point that is different from all the runs (factorial, center and axial points) already performed during experimentation, and with predicted values of Li
4SiO
4 conversion and adsorption capacity higher than the values already observed in the experimental runs. This suggests a good models predictivity.
All the optimum values obtained for the process variables (
Table 5) show that a synthesis temperature of 600 °C would lead to the highest CO
2 removal sorbent performance, whatever CO
2 concentration. This is probably due to the less sintered structure with smaller particle size of the sorbent (
Figure 2), as explained in section 3.3.1. The amount of K
2CO
3 required to maximize both Y
1 and Y
2 decrease as the CO
2 concentration increase. This result could be explained as a consequence of Fick’s law for diffusive mass transfer; the higher the CO
2 concentration the higher the driving force (concentration gradient between the bulk and the particle surface) for diffusion. Thus, the amount of K
2CO
3 required to improve CO
2 diffusivity trough the product layer is lower. Then, for low CO
2 concentration it is necessary to work with high K
2CO
3 content to contrast the low driving force with an improved diffusivity coefficient [
27,
30]. On the contrary, the optimum adsorption temperature required for maximize both responses increase with increasing CO
2 concentration. This result is associated to adsorption reaction equilibrium, so the optimum temperature would tend to be as close as possible to the thermodynamic equilibrium temperature, which increase as CO
2 concentration increases (
Figure 4), as widely discussed in section 3.3.3.
To evaluate the model adequacy and the validity of the optimization procedure, multiple confirmation runs were performed. Three additional experiments for each optimization at different CO
2 concentrations were carried out under the proposed optimal conditions reported in
Table 5. The average values of the three repeated experiments are presented in
Table 6, together with the calculated average error between the predicted values and the experimental validation tests. As shown, the experimental and predicted responses values are in agreement, with an error range of about 1-5 %. Therefore, the obtained RSM models can correlate the process variables to the Li4SiO4 conversion and to the adsorption capacity with high accuracy. These results demonstrate the suitability of RSM DoE methodology for the modeling and optimization of the CO
2 capture process on K
2CO
3-doped lithium orthosilicate sorbents.
The optimized sorbents, for the three CO
2 concentrations investigated, display excellent adsorption capacity and conversion in each sorption condition. Notably, at 4 vol% of CO
2, the optimized sorbent can capture 196.4 mg CO
2/g sorbent (corresponding to a conversion of almost 73 %) after 120 min of adsorption; the adsorption capacity increases to 239.2 and 295.6 mg/g (about 84 and 94 % of Li
4SiO
4 conversion) with increasing the CO
2 concentration to 27 and 50 vol%, respectively. These CO
2 uptakes are comparable with the adsorption capacities of alkali doped-Li
4SiO
4 sorbents reported in recent literature [
28,
29,
31,
32,
33,
34], especially considering the low CO
2 concentration (4 vol%) in the flue gas used in this work for the experimentations, while, all the other sorbent were tested at 15 or 100 vol% of CO
2.
The doped-Li
4SiO
4 sorbent optimized for maximum conversion and adsorption capacity, considering a CO
2 concentration in the treated gas of 4 vol%, was subjected to cyclic CO
2 adsorption/desorption tests to evaluate its regenerability and stability. These properties are essential for a good CO
2 sorbent in view of its use in industrial processes, such as in a fixed bed for capturing CO
2 emissions from exhaust flue gas from gas turbines. Multiple adsorption (CO
2/N
2 mixture of 4/96 vol%) and desorption (100 vol% N
2) cycles were performed in the TGA system at 1 bar. The adsorption temperature was set to 500 °C (as obtained by the optimization results) and maintained for 60 min, then the desorption process was carried out elevating the temperature to 700 °C for 30 min. The results are shown in
Figure 10, evidencing an excellent regeneration capacity for 20 adsorption/desorption cycles without activity decay.
The CO2 uptake of the first cycle is about 203 mg/g and remains almost constant for the subsequent cycles. The high adsorption capacities and the good regeneration and stability obtained at very low CO2 content suggest that K2CO3-doped Li4SiO4 sorbents could be suitable for CO2 capture in industrial applications, such as in gas-fired power plants. Furthermore, the RSM methodology could be effectively used for tailoring the sorbent synthesis conditions thus optimizing its adsorption performance for different adsorption conditions.