1 Introduction
To avoid the serious effects of global warming on both human society and the ecosystem, CO
2 reduction has become a consensus worldwide. CO
2 capture is a vital to achieving this goal, especially in the energy section. Toward the end of the 20th century, a supported potassium carbonate (K
2CO
3) solid sorbent was first proposed to successfully recovery carbon dioxide (CO
2) from flue gases, demonstrating a thermal-swing approach for CO
2 capture [
1,
2]. CO
2 is chemically absorbed by K
2CO
3, which is dispersedly attached on activated carbon, following the reaction of K
2CO
3 + CO
2 + H
2O = 2KHCO
3 to form potassium bicarbonate species. The entrapped CO
2 is released by the decomposition of KHCO
3 at the elevated temperature, and afterwards the solid sorbent is regenerated for the next cyclic operation. Before long, sodium carbonate (Na
2CO
3) was also proved to be another promising candidate for CO
2 capture from the humid flue gases [
3,
4]. The concept map for low-temperature CO
2 capture by alkali carbonate-based sorbents (ACSs), including K
2CO
3-based and Na
2CO
3-based sorbents, is illustrated in Fig.1.
ACSs are feasible for low-temperature CO2 reduction from stationary sources, especially from fossil−fuel fired power plants because they are:
1) Convenient to retrofit the existing power plants: The nature of the reversible carbonation reaction determines that the forward step prefers a high moist atmosphere with a mild temperature, which is attainable after the wet flue gas desulfurization process [
5]. In addition, the higher temperature required for the backward step, the regeneration of ACSs, is easily achieved by suction heating from steam turbine [
6]. Dual fluidized-bed reactor is demonstrated to be suitable for the dry regenerable CO
2 capture process using ACSs [
7,
8]. Therefore, the CO
2 capture unit is individually incorporated with the existing power plants without making many changes.
2) Cost-effective and energy-efficient to operate: First, the minor retrofit to the power plants costs less than other CO
2 capture technologies. In addition, the active components, K
2CO
3 and Na
2CO
3, are abundant and cheap, as well as various common supports, for instance, alumina, silica, and activated carbon. Moreover, the heat demand for the sorbent regeneration and the pressure drop for the sorbent fluidization are less than that of the liquid scrubber technologies, such as the monoethanolamine-based method [
9]. Furthermore, the efficiency penalty can be further reduced by incorporating with heat recovery facilities [
10], system integration optimization [
11], and the solar thermal power system [
12,
13].
3) Avoidable from the secondary pollution: According to the decomposition of the bicarbonate species, CO
2 and H
2O are the only two gaseous products which are both non-toxic, non-corrosive, and easy-to-collect. No secondary pollution, such as ammonia, VOC, and NO
x, is generated during the regeneration period compared to amine-based solvents or sorbents [
14].
However, some drawbacks of the alkali carbonate-based dry carbonate process should be overcome to facilitate its commercial applications:
1) Non-sufficient CO
2 capacity and sorption/desorption kinetics: Owing to the poor pore structure and scarce surface exposure, non-supported carbonates are commonly not feasible to be used for CO
2 capture, despite of their high theoretical CO
2 capacity [
15,
16]. Porous materials are applied to decrease the size of active components (K
2CO
3 or Na
2CO
3) and increase the CO
2 capacity and sorption/desorption kinetics. However, the dispersion of active components on the supports is influenced by many factors, such as the support structure, support chemical properties, and sorbent preparation methods, resulting in different CO
2 reduction performances [
17]. In general, the CO
2 capacity and sorption/desorption kinetics of ACSs are slightly worse than those of amine-based solvents or sorbents.
2) Hypersensitive to the operation parameters: Aside from the physical/chemical properties of ACSs, the sorbent performance is highly dependent on the operation parameters, which affects the real environment of the gas-solid reaction. For a dual fluidized-bed reactor system, gas inlet velocity, solid circulation rate, temperature and water vapor content are the four factors related to the CO
2 removal performance [
18].
Scientists made great efforts to improve the sorbent performance experimentally in the past two decades. Four key factors were mostly regulated in the design of sorbent, including active components, supports and additives, and preparation method (Fig.2(a)). Crystal structure and dosage of active components were commonly investigated. Zhao et al. reported that CO
2 capture on hexagonal K
2CO
3 was more kinetically favorable than that on monoclinic K
2CO
3 under the same carbonation condition [
19]. Ryu et al. manufactured different alumina-supported sorbents with various loadings of K
2CO
3 and found the CO
2 capacity of PAl40 was around 124.5 mgCO
2/g sorbent, which is much higher than that of other sorbents with lower K
2CO
3 loadings [
20]. The effects of types, pore structures, and the surface area of support were extensively studied as well [
21‒
23]. Bararpour et al. investigated the effect of textual properties, including surface area, pore size, and pore volume, on the capturing performance of mesoporous alumina supported K
2CO
3, and concluded that the larger pore size and pore volume facilitate K
2CO
3 penetration into the pores of supports during impregnation, providing more active sites for carbonation reactions and thereby achieving a better CO
2 capturing performance [
21]. Khuong et al. highlighted that the coexistence of micropores and macropores in activated carbons contributed to great CO
2 sorption on K
2CO
3-activated samples by forming solutions for carbon sequestration [
22]. Kazemi et al. successfully improved the CO
2 capacity of a
γ-Al
2O
3 supported K
2CO
3 sorbent by increasing its specific surface area, reaching 131.6 mgCO
2/g sorbent [
23]. Additives were also investigated for the specific aims, like enhancing the kinetics of carbonation [
24,
25], improving the regeneration kinetics, and decreasing the regeneration heat consumption [
26]. The wet impregnation method [
27,
28], the spray-dryer method [
29], and even direct mixing [
17] were studied in the ACSs preparation. Some theoretical modeling work were done as well to facilitate sorbent design from the atomic level by rationally simplifying the sorbent structure and the reaction conditions. The density functional theory (DFT) was first introduced to study the crystal structures of both alkali carbonate and bicarbonate species, including Li
2CO
3/LiHCO
3, Na
2CO
3/NaHCO
3, and K
2CO
3/KHCO
3. The thermodynamic properties of each crystal were calculated and the CO
2 capture by the thermal swing mode using Na- or K- carbonate/bicarbonate pair was predicted to be feasible under the post-combustion condition [
30,
31]. Potassium carbonate sesquihydrate (K
2CO
3·1.5H
2O) was investigated under the same computational condition and it was proved to be more energy-efficient using K
2CO
3·1.5H
2O to capture CO
2 than that of K
2CO
3, provided that the operation conditions are carefully maintained to reduce the dehydration of K
2CO
3·1.5H
2O [
32]. Besides the CO
2 sorption/desorption capacity and reaction rate, SO
2 and abrasion resistance are the other two key assessment criteria to evaluate sorbent (Fig.2(a)) [
24,
26,
29,
33]. Moreover, unraveling the reaction mechanism is also important to the sorbent development. At present, the intrinsic reaction mechanism is mainly validated by offline characterization methods, such as X-ray diffraction (XRD), N
2 adsorption/desorption curve, Fourier transform infrared spectroscopy (FTIR), X-ray photoelectron spectroscopy (XPS), and electron microscopy techniques (Fig.2(b)) [
28]. However, the
ex-
situ characterization results may not reveal the real species, especially the intermediates, owing to the inconsistency between the environment (such as temperature, pressure and the relative humidity) of the sorbent performance tests and that of the characterization process. For example, Luo et al. indicated that K
2CO
3·1.5H
2O was the intermediate during K
2CO
3 carbonation based on their XRD results [
34] while Jayakumar et al. suggested that K
2CO
3 hydration and carbonation occurred in parallel with no direct conversion from K
2CO
3·1.5H
2O to KHCO
3 via thermo-gravimetric analysis (TGA) [
35]. This inconsistency of the role of K
2CO
3·1.5H
2O may result from improper preservation before the
ex-
situ XRD test. Furthermore, the intrinsic kinetics of sorbents are commonly investigated in a particle-motionless reaction cell like thermal gravimetric analyzer, ignoring the effect of particle motion. However, the sorbent performance was proved to be notably influenced by flow patterns in a fluidized reactor, indicating that the reaction mechanism is dependent on momentum/heat/mass transfer [
36] (Fig.2(b)). Currently, intrusive pressure sensor is applied to reflect the average flow field in the reactor [
37], nevertheless the accurate particle behaviors are hard to analysis. DFT modeling were used to reveal the reaction mechanism as well. The surface carbonation steps on different facets of K
2CO
3 were simulated by DFT calculations [
38,
39]. The results indicate that at least two surficial carbonation pathways exist, depending on the surface textural structure [
39].
By far, neither the experimental characterization nor the theoretical simulations can unfold what exactly happen during the ACSs working period. In this perspective paper, multi-scale theoretical modeling and a compact operando measurement system are introduced to understand the operando behaviors on the sorbent chemically and physically, contributing to the design of sorbents and reactors.
2 Multi-scale theoretical modeling
CO
2 capture using ACSs belongs to non-catalytic gas-solid reaction. The operando modeling is very challengeable. First, the microstructure of solid sorbents is complex. Due to the heavy loading of alkali carbonates, the surface-active sites are diverse and the surface textual are complicated. Some active components assume to form clusters rather than be atomically dispersed on the supports [
28]. Moreover, the sorbent structure can be dynamic fluxional under an operando condition. The carbonation reaction happens at a relative low temperature (40‒70 °C for Na
2CO
3-based sorbents, and 50‒90 °C for K
2CO
3-based sorbents, respectively) with a high humidity (it is sometimes beyond a relative humidity of 90%). As a result, deliquescence and hydration, can happen, further complicating the reaction mechanism. Multi-scale theoretical simulations are a plausible understanding of a working heterogeneous catalyst system at the atomic scale [
40]. Here, DFT, ab initio molecular dynamics (AIMD), classic molecular dynamics (CMD), and machine learning (ML) technologies are proposed to realize the operando modeling of CO
2 capture on ACSs.
2.1 DFT calculations
The DFT method can be used to deal with some small systems (100‒300 atoms), which can accurately predict the surficial process (step 1 in Fig.3(a)) and bulk properties (steps 2/3 in Fig.3(a)) involving charge transfer in the thermodynamic steady-state and the quasi-steady-state (e.g., transition states).
Crystal structure refinement, surface bicarbonate formation, bulk defect pair formation, and bulk defect diffusion on Na
2CO
3 were conducted for the first time via DFT approaches, indicating that the surficial carbonation reaction was facile without the obstruction from the high surface humidity and that the generation of the ionic defect pairs (I
H+/V
Na−) in Na
2CO
3 was the rate-limiting step in the kinetics of bulk carbonation of Na
2CO
3, followed by the Na
+ cation diffusion in the NaHCO
3 layer [
41,
42]. Furthermore, the Li-doping effect on the enhancement of the I
H+/V
Na− defect pairs generation was predicted by DFT calculations. The foreseen acceleration of the initial CO
2 sorption on Na
2CO
3 was experimentally validated to be over 125% [
41]. The effect of the texture of support on the regeneration of bicarbonates was investigated at an atomic level as well. It was predicted that the surficial basic OH groups of ZrO
2 and TiO
2 were beneficial to the bicarbonate decomposition, causing the energy reduction of the sorbent regeneration. It was proved by TGA that the activation energy of KHCO
3 decomposition was decreased by 46.6% and 23.3% if adding 5% ZrO
2 or 5% TiO
2, respectively [
26,
43].
2.2 AIMD/CMD simulations
The AIMD method combines the advantages of both the first principle methods and the CMD methods. Finite-temperature dynamical trajectories were generated by using forces computed “on the fly” from electronic structure calculations, and hence, the subsequent molecular dynamics were much more “realistic” [
44]. The AIMD method was first introduced to solve the surface chemical bonding and bond-breaking between CO
2 and H
2O species on Na
2CO
3 (001) and Na
2CO
3 (−402) [
42]. It was revealed that a high coverage of H
2O molecules on Na
2CO
3 was unable to stop CO
2 binding excess H
2O or asymmetric OH groups (Fig.3(b)). To simulate the surface chemo-physical changes involving thousands of atoms is time-consuming and unaffordable by utilizing the DFT method. Moreover, AIMD simulations can give dependable dynamical trajectories in several or tens of picosecond, facilitating the analysis of some dynamic parameters. It facilitates the dynamic modeling of reaction environment and takes insights into the fluxionality of the gas-solid interface. However, AIMD methods are barely applied in the development of ACSs by far.
The CMD method based on molecular mechanics force fields was performed to demonstrate the different hydrophilicity and deliquescence between Na
2CO
3 and K
2CO
3 (Fig.3(c))). The simulation duration is longer than 500 picoseconds, which is enough to tell the discrepancy between Na
2CO
3 and K
2CO
3. Although Na
2CO
3 deformed earlier than K
2CO
3, Na
+ cations did not separate from the main cluster before 200 picoseconds. K
2CO
3 maintained a relatively regular crystal structure for about 50 picoseconds and then the crystal structure collapsed rapidly. At 500 picoseconds, most K
+ cations and
anions are separately distributed around by H
2O molecules. Based on the different deliquescence of Na
2CO
3 and K
2CO
3 predicted by CMD, a strategy was proposed to increase the CO
2 sorption of ACSs by tuning the hydrophilicity of sorbent. The experimental results confirmed that Na
2CO
3 mixed with
w(K
2CO
3) = 1.3% (
w is a mass fraction) could absorb 50% more CO
2 than pure Na
2CO
3, which was less hydrophilic [
45]. Similarly, Gayan and Vicki found that the adsorbed water could enhance surface reactivity by providing a medium for ionic dissociation [
46]. CMD methods has been used to solve many other problems, such as, particle sintering [
47], crystal thermal decomposition [
48], and molecular diffusion in pore channels [
49].
For the complicated sorbent system, which may involve active components, additives, supports and gas molecules, it is preferred to construct the more realistic particle model (beyond thousands of atoms) and perform the dynamic simulation under reasonable operation conditions (e.g., temperature and pressure) via AIMD and CMD methods.
2.3 ML approaches
As illustrated in Fig.3(d), high throughput screening (HTS) with ML approaches can contribute to figuring out candidate materials for CO
2 capture using ACSs. The accuracy of the ML algorithm relies heavily on the selection of gross descriptors. Here, the results from the theoretical simulations (DFT, AIMD and CMD) are preferred, considering the fact that uniform settings can be achieved easily on simulation while the experiments may induce a lot of anthropogenic errors. Due to the fact that the massive data are needed for ML approaches, the descriptors should be attained time-efficiently. HTS is widely implemented in the fields of battery electrolytes [
50], membranes [
51], and catalysts [
52]. However, it is still not induced into the development of ACSs. It is prospected that ML approaches can provide help in quickly finding the optimal materials and the optimal operation parameters.
3 Compact operando measurement system
The term “operando” was coined by Miguel A. Bañares, which means “operating” in Latin. Operando methods require the data from the real working condition, including temperature, pressure, and atmosphere. Besides, the data should be obtained inline or real-time. Hence, all operando methods belong to
in situ methods. The most commonly used operando method is operando spectroscopy, mainly composed of an infrared spectroscopy (IR), a Raman spectroscopy, and an X-ray absorption spectroscopy. An
in situ environmental transmission electron microscopy and an
in situ environmental scanning electron microscopy are included in operando methods though the operating pressure is currently lower than the atmospheric pressure. Operando characterization have been successfully applied for electrocatalysis [
53,
54], photocatalysis [
55], thermocatalysis [
56], and gas-solid interfacial adsorption [
57,
58]. Nevertheless, the reaction cell is designed for non-motion particles, gas, or liquids. Therefore, the existing operando methods are not suitable for gas-solid reactions in fluidized reactors. A compact operando measurement system is proposed to reveal the real kinetics of CO
2 sorption on ACSs (Fig.4).
3.1 Micro fluidized bed (MFB)
An MFB reactor is often less than 20 mm in diameter and filled with only a few to a few hundred milligrams of sample, ranging from several micrometers to several hundred micrometers in diameter. The operating gas velocity is often between 3 times to 7 times of the critical fluidization wind speed, providing an ideal mode of plug flow and reducing the effects of gas backmixing and diffusion. Thanks to the slight temperature difference between samples, the minimization of the gas backmixing, and the miniaturization of the reaction area, MFB is an ideal differential chemical reactor. The kinetics of redox performance on oxygen carrier in chemical looping combustion [
59,
60], and the kinetics of carbonation on calcium carbonate [
61] and calcium oxide [
62] are successfully obtained by MFB reactors. Two technical routes are alternatively used to determine the kinetics, including an online process mass spectrometer set at the gas outlet of the MFB reactor [
61] or a TGA holding the entire MFB reactor to monitor the weight change [
62]. The later technical route (MFB-TGA) is adopted in the proposed compact operando measurement prototype system for intrinsic kinetics analysis (Fig.4).
3.2 Multi-camera 3D imaging
Although the ideal plug flow should be kept during the experiment, the undesirable operation parameters (e.g. gas feed rate, gas composition, operating temperature) and the varying sorbent properties (e.g. particle components, size and number, moisture content, particle sphericity factor, abrasion rate) can alter the flow pattern, inducing diverse CO
2 sorption performance. A 3D particle tracking velocimetry by defocused imaging technique [
63] is proposed to online monitor particle fluidization behavior (Fig.4). Under the backlighting condition (Light source 1), the illuminated particles are imaged on two charge coupled device (CCD)-cameras (Camera 1 and Camera 2). The optical paths, which are adjusted by the light splitter, from particles to the two CCDs are unequal. Hence, the extent of defocus is dichotomous, giving rise to the determination of the particle positions. The particle velocity can be calculated from two sequential frames by cross-correlation algorithm. Besides, multi-camera tomography, can be alternatively used to obtain 3D particle tracking velocimetry as well. The particle concentration can be adjustable from dilute to dense by using particle tracking velocimetry or particle imaging velocimetry.
3.3 Operando infrared/Raman measurements
Operando spectrum technologies, including IR and Raman, are extensively implemented in food and medicine industry for process analysis. Direct spectrum detection for solids in flow is difficult since particles may not run through the detection area (the focused light spot), causing a bad signal intensity or a poor resolution. In the prototype, IR source (Light source 2 in Fig.4) is used as the backlighting for the IR Camera or IR spectrometer. Due to the strong infrared absorption effect of water and CO
2, the signals received on the IR CCD are diminished and the concentrations of gases are acquired by the inverse algorithm. Moreover, the spontaneous far IR radiation of solid sorbents and gas components (wavelength: 8‒10 μm) can be recorded by the IR camera and used to calculate temperature theoretically. However, the weak signal, the stray light, and the relatively poor photon absorption efficiency and the small number of pixels of the far IR CCD or the complementary metal oxide semiconductor (CMOS) may affect the results. Light source 3 at 785 nm is used to stimulate Raman scattering from the solid sorbents and the Raman signals are analyzed by the Raman spectrometer to monitor the generation of the solid products and the consumption of the solid reactants. The time resolution of the operando Raman measurement is supposed to be worse than that of the other measurements, because of the weak Raman effect. However, some enhancement means, including surface enhanced Raman scattering, stimulated Raman scattering, wide-angle Raman scattering collection, and optical path enlargement can be applied to improve the time resolution to seconds [
64]. In virtue of the combined Raman and IR sub-system, the new solid products (bicarbonate or carbonate hydrate) and the gas composition (mainly H
2O concentration) are monitored simultaneously. The mechanism of carbonation and hydration can be deduced, and these values can be validated by the weight fluctuation obtained by the MFB sub-system as well. By adjusting the flow pattern, the sorbent performance under the various operating condition is monitored. Moreover, the 3D particle imaging sub-system can capture some characteristic flow area, and the detailed local reaction intensity can be calculated by the combined Raman and IR sub-system. Hence, some hints for reactor design and operation condition optimization are obtained, facilitating the industrial scale-up of ACSs for CO
2 capture.
4 Conclusions and perspectives
CO
2 capture using dry regenerable ACSs are competitive and promising among various CO
2 reduction technologies. In this perspective, this paper first briefly introduced the challenges of developing more efficient ACSs: ① solid ACSs are component-complicated and sensitive to the operation parameters, resulting in massive factors to be considered in the development of sorbent, ② current techniques, including offline characterizations and particle-motionless reactor validation, are neither sufficient to reveal the intrinsic reaction mechanism nor to decouple the kinetics with the fluidization. Then, it proposed multi-scale theoretical simulations to understand the nature of reaction at the micro-level. DFT methods were demonstrated to reveal both the surface and bulk chemical behavior assuming that the sorbent model was simple. AIMD and CMD methods could be used to investigate the dynamics of sorbent in a realistic environment. ML approaches, involving DFT/AIMD/CMD results, were illustrated to facilitate material screening and ACSs development. Finally, it proposed a compact operando measurement system as well. The MFB-TGA technique allows the kinetic analysis to decouple from fluidization effects. 3D imaging method assists the understanding of the fluidization effects on chemical kinetics performance. Operando IR and Raman tools can reveal the gas composition, solid components and temperature, simultaneously. With the aid of these new research tools, it is believed that the complicated phenomena, for instance, the multi component synergistic effects, the steric-hinerance effects, and the effects of the intermediate products, the dynamic interface structure between sorbents and atmosphere can be systematically studied, thereby some inconsistency in literature can be well discussed, more convinced results can be concluded, and more efficient sorbents can be designed and manufactured. Furthermore, new ACSs are ready to be produced not only for CO
2 capture from flue gas, but direct from air. Combined Operando modeling and measurements are also very useful to investigate other heterogeneous reactions which are highly sensitive to atmosphere and flow conditions, such as the surface chemistry of environmental interfaces [
46].