Yuxuan
Xu
ab,
Yuan
Gao
a,
Liguang
Dou
a,
Dengke
Xi
a,
Chong
Qi
ab,
Baowang
Lu
c and
Tao
Shao
*ab
aBeijing International S&T Cooperation Base for Plasma Science and Energy Conversion, Institute of Electrical Engineering, Chinese Academy of Sciences, Beijing 100190, China. E-mail: st@mail.iee.ac.cn
bUniversity of Chinese Academy of Sciences, Beijing 100049, China
cNational Institute of Clean-and-Low-Carbon Energy, Beijing 102211, China
First published on 10th July 2025
CO2 dissociation, as one of the key pathways for carbon utilization, plays a critical role in sustainable carbon emission reduction. Low-temperature plasma (LTP) technology, with its highly reactive characteristics, can effectively lower the energy barrier for CO2 activation, thus facilitating efficient CO2 dissociation with reduced energy requirement. LTP-enabled CO2 dissociation has attracted widespread research interests aiming to enhance conversion performance and advance scale-up applications in recent years. Achieving such an objective necessitates dual-level considerations from both the plasma setups, such as power supplies, in situ and post-treatment control strategies, and the conversion mechanisms which encompass insights from diagnostic techniques, reaction kinetics simulations, and multi-physics modeling. Following an overview of the basic characteristics of different plasma systems and the distinctions between the CO2 dissociation processes driven by thermal catalysis and plasma catalysis, this review examines the current developments and existing limitations in LTP setups and conversion mechanisms through comprehensive analyses of experimental achievements and mechanistic investigations in LTP-enabled CO2 dissociation. Based on these findings, we propose a strategic outlook for the progression of LTP-enabled CO2 valorization from laboratory research toward scale-up implementations.
Green foundation1. This review explores recent developments in low-temperature plasma technology as an innovative green energy conversion method for CO2 dissociation, focusing on the status and limitations of plasma setups and conversion mechanisms.2. Low-temperature plasma-enabled CO2 dissociation represents a crucial sustainable carbon emission reduction strategy. Leveraging the capability of efficient CO2 valorization, it holds promising potential for broader applications across various carbon emission scenarios. 3. Future advancements require synergistic optimization of CO2 conversion performance and techno-economic feasibility for scale-up implementations. Through comparative analysis of various plasma technologies, this review provides critical guidance for facilitating holistic process optimization and maturation of plasma-enabled CO2 dissociation systems, ultimately contributing to the establishment of sustainable carbon cycling infrastructures aligned with green chemistry principles. |
CO2(g) → CO(g) + ½O2(g), ΔH = 283 kJ mol−1 | (1) |
Low-temperature plasma (LTP) technology, as a novel energy conversion pathway, represents an important alternative for CO2 conversion.41,42 The LTP is a multi-component system, with numerous charged and neutral species co-existing at different temperatures. The thermal equilibrium state of LTP can be assessed by examining the electron temperature (Te) and the gas temperature (Tg).43,44 Typically, LTP is in a non-thermal equilibrium state, where Te ≫ Tg. In terms of CO2 activation mechanisms, LTP can generate high-energy electrons (1–10 eV) and reactive species under limited energy input, thereby reducing the activation energy for CO2 dissociation without the need for catalysts.45,46 Regarding its application prospects, LTP-based CO2 dissociation, as an electric-driven reaction process with “instant-on/off” characteristics, is well suited for coupling with intermittent renewable energy sources, contributing to distributed and flexible CO2 utilization.47–51 Furthermore, considering that CO2 dissociation does not involve intricate reaction pathways in multi-reactant systems,52–54 the issue of product selectivity in plasma discharge processes is often minimal, making it highly valuable for both research and practical deployment in CO2 dissociation approaches.
Nevertheless, after examining a substantial amount of research on LTP-enabled CO2 dissociation, it has been observed that while many experimental results demonstrate high CO2 dissociation conversion rates, the corresponding energy efficiency still lacks sufficient competitiveness.55–57 Hence, balancing conversion with energy consumption is one of the key challenges in advancing LTP technology from the laboratory level to scale-up applications.58–60 Existing improvement strategies typically focus on reactor modifications,61,62 optimization of reaction parameters,63,64 and the design of plasma-specific catalysts,65,66 all of which are out of iterations within a specific plasma form. Yet, the absence of comparisons of different LTP-enabled CO2 dissociation pathways should provide new insights into addressing the aforementioned challenges from a holistic perspective.
Based on a comprehensive investigation of the current status, advantages, and drawbacks of various LTP-enabled CO2 dissociation approaches, this work provides a detailed analysis of the macroscopic plasma setups and the microscopic CO2 dissociation mechanisms. The former directly impacts CO2 conversion performance but is constrained by the core reaction steps within the latter, which suggests an interacting relationship. On top of that, we further discuss the technological bottlenecks that need to be overcome for the outlook toward industrial-scale applications of LTP-enabled CO2 dissociation.
From the perspective of discharge characteristics, DBD plasma exhibits non-equilibrium properties with significantly different vibrational temperature (Tv) and rotational temperature (Tr). Under the high reduced electric field (E/n) of several hundred Townsends (Td), a large amount of power is used to generate high-energy electrons, which trigger the electronic excitation and dissociative ionization of ground-state CO2 molecules. Such reaction pathways limit the energy efficiency of CO2 dissociation in DBD plasma. Moreover, in order to maintain stable and uniform discharge while ensuring sufficient interaction between CO2 molecules and the catalysts for high conversion rates, the throughput of the DBD plasma process is often limited.
In addition to the GAD reactor structure with quasi-two-dimensional blade electrodes, some studies have designed three-dimensional multi-blade electrodes68 and rotating GAD (RGAD) plasma reactors.69 Expanding the plasma region is considered beneficial for increasing the residence time of reactant gases, which can enhance the CO2 dissociation conversion in GAD plasma to some extent. Another approach to improve CO2 conversion is catalyst integration in the post-plasma region of the reactor, where the residual heat downstream of the discharge channel further ignites the catalysts to facilitate the conversion of unreacted CO2 molecules.
To suppress the post-plasma temperature of MWD systems, the addition of a downstream quenching structure (such as a converging-diverging nozzle, CDN) can improve heat distribution through altering the flow velocity, ensuring that the temperature downstream the core discharge region remains outside the temperature range of thermal equilibrium, thereby reducing the extent of reverse reactions. Additionally, besides coupling catalysts in the post-plasma region, similar to the approach used in GAD plasma, some studies have also explored the coupling of MWD plasma with electrocatalysis to enhance the overall conversion of CO2 dissociation.70
With regard to the CO2 dissociation reaction pathway, the vibrational excitation process remains dominant in spark discharge plasma, offering advantages in balancing the conversion rate and energy efficiency. However, due to the relatively small volume of the spark discharge channel within the reactor, enhancing the contact between CO2 reactants and the discharge region, as well as improving the conversion level within the discharge channel, represents an important research breakthrough. Some studies have pointed out that, compared to the quasi-uniform electric fields generated by plate and arc-shaped electrode structures, utilizing the highly non-uniform electric field created by the needle-plate electrode configuration, as shown in Fig. 1(d), can yield better CO2 dissociation conversion results.71
CO(g) → C(g) + ½O2(g), ΔH = 824 kJ mol−1 | (2) |
![]() | ||
Fig. 2 Thermal catalysis of CO2 dissociation. (a) Calculated mole fraction of the main molecules in CO2 dissociation with increase in temperature. (b) A typical thermal catalysis process of CO2 dissociation on the Cu(I)/TiO2−x surface. Reproduced from ref. 76 with permission from the American Chemical Society, copyright 2012. (c) Qualitative enthalpy and Gibbs energy plots presenting a basic two step mechanism of thermal-catalytic CO2 dissociation. Reproduced from ref. 77 with permission from John Wiley and Sons, copyright 2024. |
The addition of catalysts can significantly lower the operating temperature required for the thermal dissociation of CO2. Through the specific structural design of catalysts and the Langmuir–Hinshelwood (L–H) mechanism of reactant molecules on the catalyst surface, the activation energy for CO2 dissociation is substantially reduced, allowing the dissociation reaction to occur at temperatures well below those predicted by theoretical calculations. Liu et al.76 prepared a Cu(I)/TiO2−x catalyst with oxygen vacancies by pretreating it under a He atmosphere. At room temperature, CO2 molecules undergo a series of processes: surface adsorption, electron transfer, capture of oxygen atoms from CO2 by oxygen vacancies, formation and breaking of Cu+–CO bonds, and finally, the desorption of CO from the surface, completing the dissociation process. Fig. 2(b) illustrates a schematic of these key steps.
To better understand the CO2 thermal catalytic reaction process, Lefferts77 explored two key reaction steps after the introduction of catalysts:
CO2 + es → CO + Oads | (3) |
2Oads ↔ O2 + 2es | (4) |
In the above equation, “es” represents the vacancy, and “Oads” represents the adsorbed oxygen atom. Since the adsorption of CO on the metal surface is generally weaker than the dissociative adsorption of O2 molecules,78,79 the adsorption state of CO is typically not considered. The dissociative adsorption of CO2 on the catalyst surface in reaction R3 is the rate-limiting step. If this process is reversible, the composite reaction (R5) must be considered to obtain the actual reaction rate of R3. The desorption of dissociatively adsorbed oxygen atoms (R4) is always at dynamic equilibrium, with its reaction rate constant equal to the inverse of the dissociative adsorption rate constant of O2 on the catalyst surface.
CO + Oads → CO2 + es | (5) |
If the catalyst has a strong dissociative adsorption effect on O2, the reaction rate constant for R4 is much smaller than 1, meaning that the active sites of the catalyst surface are largely occupied, which suppresses the rate of R3. On the other hand, when the catalyst has a weak dissociative adsorption effect on O2, the reaction rate constant for R4 is much greater than 1, and the rate of R3 primarily depends on the concentration of CO2 molecules in the gas phase.
In terms of qualitative changes in reaction enthalpy and Gibbs free energy depicted in Fig. 2(c), the solid black line represents the CO2 dissociation process on the catalyst surface that has a dissociative adsorption effect on O2. After overcoming the activation energy barrier for CO2, it enters the Oads state, which requires further overcoming of the energy barrier from Oads to O2. The dashed black line, on the other hand, represents the reaction process on the catalyst surface with a weak dissociative adsorption effect on O2. After the CO2 dissociative adsorption, it directly transitions to O2 without going through the Oads state. The corresponding change in Gibbs free energy reveals that the transition from the dashed black line to O2 is spontaneous (ΔG < 0), while the transition from Oads to O2 on the solid black line has a positive Gibbs free energy change (ΔG > 0). This suggests that high temperatures favor the desorption of O2. As can be seen from the colored solid lines on the left side of Fig. 2(c), as the temperature increases, the enthalpy difference between Oads and O2 reduces the Gibbs energy value, gradually driving the total ΔG < 0 for the reaction to take place spontaneously.
![]() | ||
Fig. 3 Plasma catalysis of CO2 dissociation. (a) Qualitative enthalpy plot of plasma-activation of CO2. Reproduced from ref. 77 with permission from John Wiley and Sons, copyright 2024. (b) Schematic of part of CO2 electronic and vibrational levels. Reproduced from ref. 80 with permission from the Royal Society of Chemistry, copyright 2015. (c) Schematic of the basic process of plasma-activated CO2 dissociation. |
The effectiveness of plasma catalysis in reducing the reaction activation energy barrier can be attributed to its unique molecular activation pathway at the microscopic level. As shown in Fig. 3(b), compared to direct electron excitation that causes the dissociation of ground-state CO2 molecules, high-energy electrons in the plasma initiate a V–V collision process, in which CO2 molecules gradually transition from low-energy vibrational states to higher-energy states and eventually dissociate. This mechanism is vividly referred to as the “ladder-climbing process”, where lower energy input guides the CO2 molecules to step up their vibrational energy levels until they are activated and participate in dissociation.
Under the influence of these mechanisms, the basic reactions in plasma-activated CO2 dissociation include:
Electron impact ionization:
e + CO2 → e + e + CO2+ | (6) |
Electron impact excitation:
e + CO2 ↔ e + CO2(vi) | (7) |
e + CO2 ↔ e + CO2(ej) | (8) |
Dissociation upon collision with heavy particles:
M + CO2 ↔ M + CO + O | (9) |
O + CO2 ↔ O2 + CO | (10) |
The reactions R6 to R10 constitute a primary process, where “e” represents an electron, with the energy level of e on the left side of the reaction being higher than that on the right. The notation “vi” refers to the vibrational excitation energy levels of CO2, typically considering both symmetric (i = a, b, c, d) and asymmetric levels (1 ≤ i ≤ 21). “ej” indicates the electronic excitation energy levels of CO2, usually considering two levels (j = 1, 2). “M” represents any heavy particle, such as CO2, CO, O2, etc. In addition to these reactions, several important secondary processes also occur:
Recombination of CO and O/O2:
M + O + CO ↔ M + CO2 | (11) |
O2 + CO ↔ O + CO2 | (12) |
Vibrational–translational (V–T) relaxation:
M + CO2(vi) → M + CO2(vx) | (13) |
Electron detachment:
O− + CO ↔ CO2 + e | (14) |
Electron attachment:
e + CO ↔ O− + C | (15) |
In reactions R11 to R15, “vx” represents a vibrational energy level that is lower than vi.
Fig. 3(c) provides a detailed illustration of the interactions between the primary and secondary processes, as well as the key species and products involved in the reactions. Under the influence of high-energy electrons generated in the plasma, CO2 molecules in their ground state undergo electron impact ionization (R6), producing CO2+ ions. Simultaneously, electron collisions with CO2 molecules also lead to the formation of vibrationally excited CO2(vi) and electronically excited CO2(ej) species (R7 and R8). Electrons transitioning from higher to lower energy states can induce energy level upshift through V–V collisions between CO2(vi). In addition to the formation of vibrationally and electronically excited CO2 species, electronically excited CO and O2 species are also generated. Furthermore, electron loss and gain during collisions between different species result in the formation of various charged intermediates, such as CO+, O−, and others. CO2 molecules undergo dissociation reactions upon collisions with heavy particles, producing CO and O/O2 species. Apart from the primary processes, secondary processes include V–T relaxation between CO2(vi) and ground state molecules, leading to a decrease in vibrational energy levels. This process occurs at a rate lower than that of V–V collisions, with the lost energy primarily dissipated as heat into the system. The generated CO and O/O2 species can also undergo recombination processes, while electron adsorption and desorption between species like CO, O, O2, C, and O− are also significant secondary processes. Among the three products of CO2 dissociation, CO can be further utilized to synthesize a variety of chemicals, while carbon products with various morphologies find applications in membrane technologies, battery storage, drilling, and extraction industries. O2 holds significant value in specific scenarios, such as in situ resource utilization (ISRU) on Mars.
It is evident that CO2 dissociation follows different reaction pathways in thermal catalysis and plasma catalysis systems. Taking advantage of the ability to activate CO2 molecules efficiently at a lower energy input, LTP in combination with catalysts can significantly reduce the reaction difficulty of CO2 dissociation. This synergy enables CO2 conversion under milder conditions, thereby contributing to carbon reduction and the production of valuable chemicals.
Regarding the plasma power supplies, the outer ring chart of Fig. 4, marked with solid black lines, illustrates the types and proportions of plasma power supplies used to generate each discharge form. In the case of DBD plasma, alternating current (AC) power supplies are predominantly used, followed by pulsed power supplies. For GAD plasma, direct current (DC) power supplies are mainly employed, followed by AC power supplies, with some studies utilizing AC-pulsed power supplies for plasma excitation. For MWD plasma, the power supplies are primarily fixed-frequency solid-state microwave supplies, predominantly operating at 2.45 GHz, although 24 GHz and 915 MHz frequencies are also used.
In terms of in situ plasma regulation techniques, the second layer of the annular chart from the inside out in Fig. 4 lists the main in situ techniques corresponding to the three most used plasma forms. For DBD plasma, modulation towards plasma setups mainly relies on the addition of catalysts and the adjustment of dielectric properties of the barrier materials, with some studies exploring time-cycling strategies for in situ regulation. As for GAD plasma, magnetic fields can be applied to the discharge channel, or the vortex flow field can be created to alter the fluid dynamics distribution, thereby enabling in situ regulation of reaction parameters. In situ control strategies for MWD plasma typically focus on manipulating gas pressure and optimizing waveguide design.
As shown in the outermost ring chart of Fig. 4, post-plasma coupling strategies primarily concern reaction systems involving GAD and MWD plasma, where catalysts cannot be directly placed in the discharge region. For GAD plasma, there have been studies involving the use of catalysts or carbon beds for coupling downstream, while for MWD plasma, CO2 dissociation can be enhanced through integration with electrocatalysis or quenching structures.
Plasma power supply | Waveform | Discharge type | Output parameters | Ref. |
---|---|---|---|---|
DC |
![]() |
GAD | ▶ 5–15 kV | 122, 124, 126, 129, 130, 132, 134 and 137 |
Corona | ▶ 8–40 W | 120 | ||
Glow | ▶ 30 V | 138 and 139 | ||
Micro-slit | ▶ 400–850 V | 144 | ||
AC |
![]() |
DBD | ▶ 3–15 kV ▶ 4–100 kHz | 81–102 |
GAD | ▶ 0.25–15 kV ▶ 20–60 kHz | 123, 125, 127, 128, 131, 134 and 135 | ||
Glow | ▶ 2–15 kV ▶ 2–10 kHz | 72 and 140 | ||
Pulsed |
![]() |
Spark | ▶ 6–20 kV ▶ 30–100 kHz | 141–143 |
DBD | ▶ 5–30 W | 83 | ||
Corona | ▶ 35–51 kV ▶ 20–200 Hz | 119 | ||
Non-self-sustained | ▶ 10 kV ▶ 20 kHz | 121 | ||
AC-pulsed |
![]() |
GAD | ▶ 10 kV ▶ 0.05–40 kHz | 133 and 136 |
Spark | ▶ 10 kV ▶ 20 kHz | 71 | ||
MW power supply |
![]() |
MWD | ▶ 2.45 GHz ▶ 0.2–2 kW | 103, 105, 107, 109, 110 and 113–118 |
▶ 24 GHz ▶ 0.4–5 kW | 104, 106 and 111 | |||
▶ 915 MHz ▶ 1–10 kW | 108 and 112 |
DC power sources are predominantly used to drive GAD plasma. To prevent the discharge channel from transitioning into a high-current arc mode and to avoid overheating of the reactor, a current-limiting resistor is typically added.122 During operation, the current amplitude generally ranges from mA to A. Additionally, the GAD process involves a repetitive arc initiation–gliding–extinction cycle, and its frequently changing load characteristics place high demands on the stability of the power supply. When DC power supplies are used to drive low-intensity discharge forms such as corona discharge, glow discharge, and micro-slit discharge plasmas, the output voltage is typically low.
AC power supplies are mainly used in DBD plasma systems because the load of a DBD reactor can be approximated as a capacitor, which has a blocking effect on DC. The setting of the AC output frequency directly influences the level of power injection. Increasing the frequency while maintaining constant output voltage can increase the power, often resulting in enhanced CO2 conversion. When the output power is fixed, increasing the frequency can make the discharge more uniform, but it does not necessarily promote the CO2 conversion rate.84,89
Pulsed power supplies can facilitate a variety of plasma forms. Among them, nanosecond pulsed power supplies, with their short pulse width, low duty cycle, and high instantaneous power density brought by fast pulse rise times, offer unique advantages in terms of energy efficiency. For instance, Yong et al.142 used a nanosecond pulsed power supply in spark discharge plasma for CO2 dissociation, achieving a 16% CO2 conversion rate with 23% energy efficiency. AC-pulsed power supplies, which combine both AC and pulsed forms, can drive GAD and spark discharge plasmas and show promising energy efficiency.
Microwave power supplies typically have a fixed output frequency, but different waveguide and resonator configurations lead to various reflection processes. Since microwave energy directly interacts with the plasma region in the form of electromagnetic fields, power levels are typically used to assess the strength of the discharge. The output power of microwave power supplies generally ranges in the level of kilowatts (kW). Common 2.45 GHz microwave power supplies can achieve an output efficiency of around 90% (the ratio of plasma absorbed power to the total output power).104 Increasing the frequency (e.g., to 24 GHz) can enhance the specific energy input (SEI) and plasma density, while the output efficiency may decrease.145
LTP can be modeled as an electrical circuit load based on the discharge form, and the degree of matching between the power supplies and the plasma load determines its output efficiency, which is essentially the conversion efficiency of electrical energy. Lisi et al.82 considered an equivalent circuit model (ECM) that couples the internal circuit of the power supply with the load characteristics of the plasma before and after discharge. They highlighted the importance of frequency regulation in stabilizing the non-linear processes of plasma discharge.
However, current explorations on LTP-enabled CO2 dissociation using various power supplies demonstrate insufficient analysis of the matching relationships between plasma power supplies and plasma loads (typically the plasma reactors) based on ECMs. From the perspective of efficient energy conversion, two critical challenges persist. On the one hand, the transmission efficiency of electricity from power supplies to plasma reactors, as an important determinant of effective energy utilization in electrically driven plasma processes, is frequently overlooked. This efficiency metric is particularly critical for evaluating the whole energy consumption during industrial-scale process optimization. On the other hand, plasma discharge parameters are intrinsically governed by the output characteristics of power supplies. Establishing fundamental circuit-level correlations between power supplies and plasma loads could advance the optimization of critical parameters, including power supply configurations and reactor designs. Overall, the ECMs of different plasma loads have not yet been studied systematically and in depth. For DBD plasma with planar or coaxial reactor geometries, relatively well-developed resistor-capacitor models exist, with key circuit parameters that can be determined from the Lissajous characteristics of stable discharge.146,147 However, for plasma systems with transient and nonlinear nature (spark discharge and GAD) or complex and variable reactor geometries (MWD), it is very difficult to extract equivalent circuit parameters for ECM analysis. In these cases, computationally expensive numerical methods or electromagnetic simulations remain preferred over circuit approaches.148–150 Nowadays, as emerging artificial intelligence (AI) technologies have catalyzed the paradigm-shifting research perspectives, which may facilitate the elucidation of source-load matching mechanisms through iterative learning of multi-dimensional parameter datasets. AI models, especially machine-learning (ML) algorithms based on various neural network frameworks, can significantly streamline the process of extracting and optimizing equivalent circuit parameters that are critical for ensuring efficient power transfer and desired frequency responses.151,152 For future improvements toward efficient and in-depth ECM investigations, leveraging a hybrid approach that combines first-principles physics with ML algorithms to overcome analytical intractability in parameter extraction is a promising way. Such advancements could ultimately enable circuit-level optimization strategies for plasma-enabled CO2 dissociation systems.
Moreover, the miniaturization and cost control of the power supplies used in laboratory-scale systems often represent one of the major challenges in scaling up LTP-based CO2 dissociation systems for practical applications. Xu et al.71 attempted to use an improved cost-effective commercial neon power supply to drive spark discharge plasma for CO2 dissociation, achieving a good balance between conversion rate and energy efficiency. For other plasma power supplies, deriving key electrical parameters through ECM analysis may help in the development of compact plasma power supply systems.
As shown in the statistical results listed in Table 2, alkaline metal oxides and various forms of SiO2 materials (such as glass beads, quartz wool, and quartz sand) are the most commonly used packing materials in DBD reactors, with the dielectric constant being a key factor influencing the discharge characteristics after packing. A higher dielectric constant not only enhances the electric field strength within the discharge channel,80 but also leads to an increase in Tgvia Joule heating.153 Xu et al.86 directly filled BaTiO3 beads between the electrodes, which have a relative dielectric constant of around 10000. This allowed discharge in an originally non-breakdown gap, and after doping with Ar as the carrier gas, the conversion rate reached 36%. The addition of Ar effectively lowers the breakdown voltage, and its metastable states alter the electron energy distribution function (EEDF). The increased electron density is favorable for electron-excitation dissociation of CO2 in DBD plasma. However, Ar dilution also consumes part of the discharge energy, so the overall energy efficiency is reduced. Further attachment of catalyst materials to the surface of materials with a high dielectric constant is proven to further enhance the CO2 dissociation conversion effect in DBD plasma. Hosseini Rad et al.91 packed glass beads coated with CeO2 catalysts into the reactor, achieving a CO2 conversion rate of 12.3% under 2-bar pressure, compared to 9% for the blank control, with the corresponding energy efficiency increasing from 16% to 20%.
Dielectric | Catalyst | Conversion before (%) | Conversion after (%) | Ref. |
---|---|---|---|---|
BaTiO3 | — | 5 | 18 | 81 |
Ni/SiO2 | 18 | 24 | ||
— | — | 36 | 86 | |
— | 5.2 | 20 | 87 | |
— | 1.3 | 1.7 | 102 | |
Glass balls | — | — | 12 | 83 |
— | 9 | 8 | 91 | |
CeO2 | 9 | 12.3 | ||
— | 12.5 | 16 | 101 | |
SiO2 | — | — | 4 | 93 |
— | 5.2 | 5.5 | 87 | |
ZrO2 | — | 5.2 | 5.8 | |
— | 5.2 | 19 | ||
Al2O3 | — | 26 | 42 | 96 |
— | — | 4.9 | 93 | |
— | 0.6 | 1.4 | 102 | |
γ—Al2O3 | — | 14.5 | 16.3 | 97 |
— | 12.5 | 16 | 101 | |
α-Al2O3 | — | 12.5 | 15 | |
Foamed Ni | — | — | 48.6 | 92 |
Foamed Cu | — | — | 49.2 | |
Ca0.7Sr0.3TiO3 + Li2Si2O5 | — | — | 15.5 | 93 |
CaTiO3 | — | 12.5 | 20.5 | 101 |
CaO-B2O3-SiO2 | — | — | 48.7 | 98 |
CaO | — | 14.5 | 32.8 | 97 |
Quartz wool | — | 14.5 | 24.5 | |
Quartz sand | — | 14.5 | 15.2 | |
MgO | — | 14.5 | 23.7 | |
Silica gel | — | 12.5 | 14 | 101 |
In addition to spherical dielectric materials, the use of mesh-type porous materials in the gas gap has also shown good performance. Zhu et al.92 filled a DBD reactor with porous Cu foam, which serves as an effective channel for electron and energy transfer. The resulting micro-discharges were evenly distributed within the gas gap, achieving a CO2 conversion rate of up to 49.2%. Beyond enhancing the discharge, the packing materials can also act as sorbents for CO2 storage, enabling time-sequenced manipulation of the plasma CO2 dissociation process. Li et al.85 used Mg/Al hydrotalcite powder with particle sizes ranging from 250 to 355 μm as a CO2 sorbent, which was packed into two identical DBD reactors. After one reactor (A) completely absorbed CO2, discharge was initiated within, and the products along with unreacted CO2 were transferred to another reactor (B) for CO2 desorption. CO2 was then introduced into reactor B to saturate the sorbent, and after discharge and CO2 conversion, the mixed gas was passed into reactor A to absorb unreacted CO2. Followed by CO2 being fed into reactor A until saturation was reached, the cycle was repeated with the steps above. With Ar dilution, a single cycle achieved a conversion rate of up to 67.44%. However, the time–space distribution process for conversion is inefficient in terms of plasma energy utilization, resulting in an energy efficiency below 1%, which is insufficient for practical applications.
The in situ modulation strategies aimed at enhancing the CO2 conversion in GAD plasma primarily focus on increasing the coverage volume of the arc channel within the reactor, thereby increasing the contact probability with reactant molecules. One approach involves utilizing the electromagnetic field characteristics of the arc itself, where an external magnetic field is applied to guide the evolution of the discharge channel. Ivanov et al.128 investigated the effects of different directions of applied magnetic fields on the performance of a GAD configuration with blade electrodes. The results indicated that in the magnetic-accelerated GAD configuration, the direction of the Lorentz force acting on the arc aligns with the direction of the gas flow (i.e., the direction of arc propagation). This configuration shortened the time between arc ignition and extinction, resulting in higher CO2 conversion rates at high flow rates (8–14 L min−1), though with a decrease in discharge stability and energy efficiency. In contrast, in the magnetically stabilized GAD configuration, the direction of the electromagnetic force was opposite to that of the gas flow, stabilizing the discharge channel within a certain range. At low flow rates (1–7 L min−1), this setup achieved higher CO2 conversion rates and energy efficiency, both of which decreased significantly as the flow rate increased. Zhang et al.129 studied the discharge characteristics and CO2 dissociation conversion of a magnetically assisted RGAD plasma driven by a DC power supply. Under the influence of the applied magnetic field and tangential gas inlet, the discharge channel rotated at a frequency of approximately 200 Hz, achieving CO2 conversion rates of 4–4.4% and energy efficiencies of 16–17% at flow rates of 6–7 L min−1.
Another in situ modulation technique for GAD plasma is the utilization of the flow field distribution induced by a turbulent vortex to interfere with the movement of the discharge channel. Compared to the effects of an external magnetic field, the influence of the vortex extends beyond the discharge channel to affect the entire reactor space and various stages of the reaction process. Nunnally et al.134 employed a reactor structure that generated reverse vortex flow, achieving CO2 conversion rates of 2–9% at high flow rates (14–40 L min−1), with a peak energy efficiency of 43%. The key to achieving high energy efficiency was that the reverse vortex field improved the mixing of the discharge region with the reactant gas, increasing the residence time. Additionally, under the influence of the rotating gas flow, products were rapidly moved away from the core high-temperature region, thereby suppressing the extent of recombination reactions.154 Thus, designing adjustable external magnetic fields and coupling them with appropriate vortex-generating structures may significantly enhance both the CO2 dissociation conversion rates and energy efficiency in GAD plasma systems.
In MWD plasma systems, the in situ regulation methods are limited due to the fixed resonance cavity structure parameters at specific microwave frequencies. One of the key factors affecting MWD plasmas is the gas pressure, which typically needs to be low (around 0.15 atm) to maintain the plasma in a diffuse discharge form, ensuring that the system remains in a non-equilibrium state.155 In specific scenarios with naturally low atmospheric pressures, such as in Martian ISRU, the atmospheric pressure on Mars is around 0.1 atm, and the atmosphere is composed of 96% CO2,156,157 making it particularly suitable for MWD processes. Kelly et al.110 conducted MWD for CO2 dissociation experiments in a laboratory-simulated Martian atmospheric environment, achieving CO and O2 production rates of 76.1 g h−1 and 47 g h−1, respectively, which are about 30 times higher than the yields of NASA's current MOXIE experiment.
Although the resonator dimensions are relatively fixed, some studies have employed surface wave excitation structures, such as coaxial waveguide and surfaguide,158,159 to facilitate MWD.
Qin et al.105 utilized a coaxial resonator under atmospheric pressure to concentrate microwave energy within a central quartz tube to generate MWD plasma, achieving a CO2 conversion rate of 50% and an energy efficiency of 5.25%. Silva et al.113 used a surfaguide as the microwave energy transmission medium, achieving a high local CO2 conversion rate of 80%, though the CO2 conversion was unevenly distributed within the quartz tube.
To sum up, the in situ regulation techniques for different plasma types influence discharge parameters, heat and mass transfer processes, and energy distribution to varying degrees. DBD plasma systems demonstrate inherent suitability for in situ modulation strategies in CO2 dissociation systems, owing to their unique reactor configurations and discharge characteristics. The introduction of tailored dielectric packing materials enables versatile control over plasma-chemical reaction pathways, often yielding promising improvements in CO2 conversion performance. However, current methodologies for determining the types, compositions, packing densities, and spatial distributions of the packing materials primarily rely on empirical trial-and-error approaches across extensive material libraries, lacking systematic optimization aligned with specific plasma parameters. To address this limitation, integrating AI-assisted screening protocols with targeted design strategies emerges as a promising pathway to enhance optimization efficiency.160–163 The AI techniques for catalyst material selection have primarily relied on ML frameworks that integrate both supervised learning models, like convolutional neural networks trained on large-scale material property databases, and active learning algorithms that predict material properties. The efficacy of such AI-driven frameworks, however, is contingent upon the extraction of critical process descriptors and the establishment of sufficiently representative datasets including composition and structure–performance databases, and density functional theory-calculated data, helping identify promising materials efficiently. Similarly, AI models offer unique advantages in capturing the nonlinear dynamics of multi-component plasma systems during reactions and thus assist in determining plasma parameters that favor conversion performance.164,165 For GAD and MWD plasmas, the reactions occurring within the discharge channels are often insufficient to significantly alter the overall conversion performance. This is due to the decline in Tg outside the plasma core region and the energy decay of active species, which leads to changes in the dominant reaction pathways in those areas. Therefore, considering modulation strategies for the post-plasma region is of significant importance for enhancing CO2 dissociation.
![]() | ||
Fig. 5 Post-plasma techniques for CO2 dissociation in GAD and MWD plasma systems. (a) Downstream catalyst addition. Reproduced from ref. 137 with permission from Springer Nature, copyright 2020. (b) Carbon bed integration in GAD plasma systems. Reproduced from ref. 125 with permission from Elsevier, copyright 2022. (c) Converging-diverging nozzle facilitated quenching. Reproduced from ref. 109 with permission from Elsevier, copyright 2023. (d) The coupling with electrolysis in MWD plasma systems. Reproduced from ref. 70 with permission from Elsevier, copyright 2022. |
In addition to using the post-plasma energy to activate catalysts and enhance CO2 dissociation, another important post-plasma strategy is to remove products promptly in order to promote the forward reaction based on Le Chatelier's principle.168 As shown in Fig. 5(b), Girard-Sahun et al.125 coupled a carbon bed with a continuous feed silo in the post-plasma region of GAD plasma, and monitored the carbon bed temperature distribution using thermocouples. The results showed that the heat generated by GAD caused O/O2 to react with carbon (the so-called reverse Boudouard reaction) which allows the consumption of the undesirable O/O2 and to enrich the exhaust in CO. It not only increased CO2 conversion from 7.6% to 12.6%, but also led to a rise in CO concentration from 7.2% to 21.9%, as the removal of O/O2 further generated CO. Overall, the energy efficiency increased from 27.9% to 45.4%. Furthermore, some studies have explored the effects of oxygen carriers to achieve complete O2 removal in plasma chemical cycles, but a much lower throughput (0.025 L min−1) was required compared to coupling with carbon beds.169
In MWD plasma systems, one of the challenges is addressing the decline in energy efficiency due to the excessively high Tg value in the discharge region. Specifically, in atmospheric pressure MWD processes, the Tg in the core region of the plasma torch can reach 6000 K. In the nearby regions, temperatures range from 2000 to 3000 K, which promotes the recombination of CO and O/O2.106 While reducing the operating pressure can help lower the Tg in the core region and reduce the temperature gradient within the reactor, it is still necessary to regulate the temperature distribution to enhance CO2 dissociation performance. Mercer et al.109 utilized a CDN placed at a specific position downstream of the MWD region to generate a vortex flow geometry (Fig. 5(c)). This approach optimized the temperature field in the plasma region while enhancing the radial driving force for transporting products to the cooler downstream region, thus achieving a simultaneous improvement in both CO2 conversion and energy efficiency. MWD plasma can also be integrated with other CO2 conversion methods. For example, Pandiyan et al.70 integrated an electrocatalytic process into the downstream of the MWD plasma, where a solid oxide electrolysis cell was used to further convert the CO2, CO, and O2 mixed gases generated by the plasma, thereby separating CO and O2 (Fig. 5(d)). However, there is still some CO loss at present, and optimization of electrode materials is needed to improve the effective yield of CO2 dissociation in MWD plasma coupled with electrocatalysis.
CO2 conversion rate:
![]() | (16) |
CO production rate:
![]() | (17) |
SEI:
![]() | (18) |
![]() | (19) |
It should be noted that the molar value “24.5 (L mol−1)” in R17 and R19 is only valid at 298 K and atmospheric pressure; for other conditions, certain calibrations are required.
Energy efficiency of CO2 dissociation (R1):
![]() | (20) |
Here “ΔH298” (283 kJ mol−1 or 2.93 eV per molecule) represents the enthalpy change of the CO2 dissociation reaction (R1) at 298 K.
Energy cost of CO2 dissociation (R1):
![]() | (21) |
Based on equations R16 to R21, the CO2 conversion performance under different plasma setups can be evaluated, with the two most critical metrics being XCO2 and η. It should be noted that not all the parameters are listed by the authors, but they can be calculated using the experimental data provided in these articles, ultimately yielding all the evaluation metrics. As shown in Fig. 6(a), (c) and (e), the XCO2 and η achieved by the three types of plasma systems under different SEIs are statistically summarized.
For DBD plasma, XCO2 increases with the SEI, which typically exceeds 40% within the range of 3–110 eV per molecule. In terms of η, higher SEIs generally result in lower η for CO2 dissociation. However, at lower SEIs, some studies reported η reaching around 40%. The work conducted by Wang et al.98 accomplished the optimal balance between XCO2 and η. In their study, they packed the DBD reactor with microstructured Ca0.8Sr0.2TiO3 materials, and when the addition weight percentage was 5 wt%, XCO2 and η reached 48.71% and 44.62%, respectively.
For GAD plasma-activated CO2 dissociation, a notable difference compared to the DBD plasma system is its generally lower (less than 30%) and smaller SEI (less than 4 eV per molecule). However, it exhibits a significant advantage in terms of η, which can reach an average of 30% and a maximum of 61.55%.123 Furthermore, Zhang et al.122 incorporated biochar into the GAD plasma system, leveraging the synergistic effects of both plasma chemistry and thermochemistry to achieve a balance between XCO2 (27.1%) and η (39.59%).
MWD plasma for CO2 dissociation achieves the highest average XCO2, reaching approximately 34%. As the SEI increases from 0.5 to 100 eV per molecule, a significant enhancement in conversion performance is observed, with a maximum XCO2 of up to 92% (with Ar dilution).115 Furthermore, the MWD plasma also demonstrates relatively high η within the lower SEI range. Chen et al.112 utilized surface wave excitation to drive MWD, achieving an η of up to 56% with an SEI of around 2.9 eV per molecule after coupling with a downstream NiO/TiO2 catalyst chamber. The corresponding XCO2 reached 45%, resulting in a favorable balance of CO2 dissociation performance within this system.
Based on the estimation by Snoeckx and Bogaerts,170 LTP-based energy conversion technologies for syngas production need to achieve an η of at least 60% in order to be competitive with existing mature processes, such as thermal catalysis and electrolysis. This η target can also be applied to the CO2 dissociation field. While factors like renewable energy integration and modularization facilitate the competitiveness of LTP technology even below this target, it still provides a relevant benchmark. As indicated by the red dashed lines in Fig. 6(a), (c) and (e), there is still a gap between the η achieved by DBD plasma for CO2 dissociation and this target, while GAD and MWD plasmas have already met or are close to meeting this goal. However, for MWD plasma, although there are reported results exceeding the 60% efficiency target and even reaching 90%,155,171,172 these results were not reproducible in subsequent studies. In the recent work of van Rooij et al.,117 they did not obtain comparable energy efficiencies (the maximum value was near 50%) using conditions like those in high-efficiency studies. It is worth noting that, with the plasma setup tailored to achieve non-equilibrium conditions throughout the whole discharge region, theoretical models indicate that the energy efficiency could reach 90% via the vibrational excitation mechanism.173,174 However, because the aforementioned studies did not provide sufficient evidence to substantiate the underlying mechanism, they were not included in the statistical analysis presented in this work.
In addition to discussing the XCO2 and η associated with SEI, practical energy conversion applications often focus on the process throughput, which refers to the amount of CO2 processed per unit of time. This is closely related to the CO2 inlet flow rate. Therefore, we further analyzed the distribution patterns of XCO2, YCO, and ECCO2 of CO2 conversion for the three plasma forms, as shown in Fig. 6(b) and (d).
Regarding the CO2 flow rate, the DBD plasma typically operates with an inlet flow rate of less than 1 L min−1, while the flow rates in MWD and GAD plasmas can reach up to 75–80 L min−1, which are significantly higher than that of the DBD plasma. With regard to XCO2, the normalized normal distribution curves of the XCO2 for different plasma types are plotted on the right axis of Fig. 6(b). The peak height of these curves corresponds to the distribution density, while the position of the peak along the vertical axis indicates the mean XCO2 (expectation value). It can be observed that the XCO2 distribution for MWD plasma is quite scattered, with the highest mean XCO2, while the conversion rate for DBD plasma is more concentrated, with a lower mean XCO2 than that of the MWD plasma. The XCO2 for GAD plasma is the most concentrated but at a lower level. However, due to the larger throughput of GAD plasma setups, the normalized log–normal distribution curve of YCO shown on the right axis of Fig. 6(d) indicates that the average YCO in GAD plasma is higher than that in the other two plasma forms, with a more concentrated distribution. The average YCO in MWD plasma is lower than that in GAD plasma but still significantly higher than that in DBD plasma.
The ECCO2 metric which integrates energy consumption with XCO2, offers a new perspective for evaluating the CO2 conversion performance of the three plasma setups. As shown in Fig. 6(f), the normalized log–normal distribution curve of ECCO2 plotted on the right axis reveals that the GAD and MWD plasmas exhibit relatively similar average ECCO2. However, the ECCO2 distribution for GAD plasma is more concentrated across different flow rates. In contrast, DBD plasma shows a significantly higher ECCO2 compared to the other two plasma types.
In summary, GAD and MWD plasmas have demonstrated the potential for mass-scale applications in terms of energy efficiency and processing throughput in existing studies, yet there exist inherent limitations in the conversion rates of CO2 dissociation. To achieve complete feedstock utilization, the systematic implementation of post-reaction gas recycling protocols becomes imperative. In contrast, DBD plasma demonstrates a more balanced performance concerning conversion rates and energy efficiencies. However, scaling DBD reactors may require distributed cascades to increase the processing capacity, which introduces non-trivial trade-offs: (1) increased consumption of dielectric packing materials, raising operational costs and waste management challenges; (2) potential degradation of conversion rates and energy efficiencies due to inhomogeneous field distribution and transport limitations at elevated flow rates. These scaling-induced compromises underscore the critical need for multi-objective optimization frameworks that reconcile throughput enhancement with sustained CO2 conversion performance.
Through optical emission spectra (OES), it is not only feasible to effectively identify intermediate species involved in the CO2 dissociation process but also to calculate critical parameters such as electron density (ne), Tv, and Tr (which reflect local Tg) by fitting specific spectral lines corresponding to species like CO, C2, and those introduced by added carrier gases such as Ar and N2. These parameters are closely related to the microscopic mechanisms of the CO2 dissociation reaction.
In GAD plasma systems, due to the continuous movement of the discharge channel, OES diagnostics are generally not applicable.175 To illustrate the differences in the CO2 dissociation OES across different plasma systems, a comparison is made between the OES results of spark discharge plasma, which is similar to the GAD plasma form in terms of discharge parameters, and those of DBD and MWD plasmas.
As shown in Table 3, in the DBD plasma system, a significant presence of the CO2+ Fox–Duffendack–Barker (FDB) band system is observed, along with the Ångström band system of CO and atomic peaks of O. In the spark discharge plasma, the CO2+ species are nearly identical to those in DBD, but CO (A1Π → X1Π) and CO (c3Π → a3Π) species are also present. In the MWD plasma, the predominant species shift to the C2 Swan band system, specifically (0, 0), (1, 0), and (0, 1), with band heads appearing at wavelengths of 516 nm, 473 nm, and 562 nm, respectively. Additionally, the CO Ångström band system is still observed in the range of 400–700 nm, along with the atomic peak of C. The differences in the composition of these species can be attributed to the varying intensities of the discharges. For instance, the formation of C2 species in MWD plasma indicates a more intense CO2 dissociation process under the influence of the discharge.
Plasma type | Main species | Wavelength (nm) | Ref. |
---|---|---|---|
DBD | CO2+ | ▶ 288.6 ▶ 351 ▶ 367 ▶ 385 ▶ 434.2 | 82, 90, 99 and 100 |
CO | ▶ 451 ▶ 471.7 ▶ 483 ▶ 520 ▶ 561 ▶ 579 ▶ 610 | ||
O | ▶ 751 ▶ 777 ▶ 845 | ||
Spark discharge | CO2+ | ▶ 288.6 ▶ 404.7 ▶ 410.9 ▶ 434.2 | 71 |
CO | ▶ 191.4–271.1 ▶ 447.9 ▶ 471.7 | ||
▶ 483.5 ▶ 519.8 ▶ 561 ▶ 607.9 | |||
O | ▶ 777.5 ▶ 844.7 | ||
MWD | C2 | ▶ 420–570 | 105, 108, 112–116 and 118 |
CO | ▶ 400–700 | ||
C | ▶ 249 | ||
O | ▶ 777.4 ▶ 844.6 |
An in-depth review of studies on 0D reaction kinetics of CO2 dissociation in DBD plasma,84,86,91,173,176–179 GAD plasma,125,132,136,175,180,181 and MWD plasma103,107,110,117,182–184 has been conducted, summarizing the CO2 reaction pathways in these three plasma systems as shown in Fig. 7.
![]() | ||
Fig. 7 Schematic overview of CO2 dissociation mechanisms of different plasma systems. The width of the arrow line indicates the importance of the corresponding reaction path. |
As shown in the upper part of Fig. 7, in the DBD plasma system, electronic impact dissociation of CO2 molecules is the dominant dissociation pathway. The primary mechanism involves electron collision with ground-state CO2 molecules, leading to direct dissociation into CO and O. Additionally, CO2 may undergo ionization, producing CO2+ ions, or dissociate into O− ions. Furthermore, a small portion of CO2+ ions undergo dissociation through electron collisions, generating CO and O. The O/O− species can collide with other molecules, further generating O2 or O3, and O2 may also undergo collisions with electrons, leading to the formation of O2− ions. The primary recombination reaction pathway involves collisions between CO and O with other molecules. This reaction mechanism results in relatively low electron energy utilization efficiency. However, filling the discharge region with dielectric materials primarily aims to increase the number of high-energy electrons during discharge, thereby promoting the dissociation of more CO2 molecules. Although direct modeling of packed materials is not feasible within 0D models, attention to the related changes in simulation input parameters (particularly the reaction set, cross-section data, and reaction coefficients) is vital for accurately reflecting the conversion performance of DBD plasma-enabled CO2 dissociation. In the work by Ponduri et al.,185 the 0D simulation results were found to be in reasonable agreement with experimental observations after incorporating Itikawa's dissociation cross-section data186 alongside the electron excitation cross sections reported by Phelps with a threshold of 7 eV.187
The lower left side of Fig. 7 summarizes the CO2 dissociation mechanism in the GAD plasma, as revealed by the 0D reaction model. Unlike the DBD plasma system, the dissociation of CO2 in the GAD reactors depends on the electronic excitation of ground-state CO2 molecules, which undergo a V–V process to reach a vibrationally excited state (mainly in the asymmetric mode). These vibrationally excited molecules then undergo electron collisions, resulting in the formation of CO and O particles. O can further collide with CO2 at high-energy vibrational states to generate CO and O2. A small number of vibrationally excited CO2 molecules also participate in the generation of ions such as CO2+, CO3−, and O−. The formation of CO2 is mainly attributed to the V–T relaxation process of vibrationally excited molecules, followed by recombination through collisions between CO and O/O2 with other particles. Through the vibrational excitation mechanism, the CO2 dissociation path in the GAD plasma exhibits higher energy transfer efficiency. In practice, the 0D modeling results of Sun et al.136 emphasized the importance of the vibrational-excitation mechanism in GAD plasma. The simulation outcomes exhibited good agreement with experimental measurements, indicating that the identified key species and reaction pathways offer critical guidance for deriving metrics that govern CO2 dissociation performance. Specifically, constraining the gas temperature within the arc channel and increasing the power density serve to elevate the significance of non-equilibrium vibrational excitation in the CO2 dissociation pathways, while facilitating the prompt removal of O atoms from the system helps suppress the recombination of CO and O. Besides, the corresponding macroscopic control strategy is primarily focused on enhancing the interaction between the plasma region and the gas flow to maximize the utilization of this efficient dissociation mechanism.
The lower right side of Fig. 7 illustrates the CO2 dissociation path in the MWD plasma. The 0D simulation results for this system show variability under different experimental conditions. Some studies suggest that at lower pressures, CO2 molecules dissociate primarily through vibrational excitation, mainly producing CO2 at different vibrational modes in the symmetric form, with a small portion excited to higher energy levels in the asymmetric mode.107,110,183 Other research indicates that at low pressures, electron-collision dissociation dominates, and as pressure increases, the discharge transitions into a torch mode. The V–T relaxation process hinders the formation of vibrationally excited species, and the dissociation mechanism shifts to one predominantly governed by thermal effects, especially involving collisions between CO2 and O or other molecules.178,182,184 In the work of Kozák and Annimie,183 the 0D simulation results accurately captured the trends in CO2 conversion and energy efficiency as functions of SEI and ne, and the calculated Tg aligned well with experimental results. Further analysis revealed that increasing ne and reducing Tg are critical for limiting the V–T relaxation processes that deplete vibrational levels, thereby enhancing energy efficiency. However, once the gas temperature exceeds a certain threshold, further improvements in energy efficiency become marginal, as the non-equilibrium vibrational excitation mechanism is increasingly suppressed. Therefore, when adjusting microwave discharge parameters, it is essential to carefully control operating pressure and gas flow distribution to avoid excessive heating in the plasma region, which could lead to a decrease in the energy efficiency of CO2 dissociation.
Overall, a systematic analysis of the dominant reaction mechanisms governing CO2 dissociation across distinct plasma systems reveals critical modality-specific energy transfer pathways. In DBD plasma systems, CO2 dissociation predominantly proceeds via direct electron-impact dissociation mediated by high-energy electrons, whereas GAD and MWD plasma systems additionally exploit vibrational excitation mechanisms targeting ground-state CO2 molecules. This distinction in energy transfer strategies underpins the superior energy efficiency observed in GAD and MWD systems compared to DBD systems.
Combining the plasma setups and the corresponding conversion performances, it further elucidates the following mechanistic drivers: DBD reactors can achieve elevated CO2 conversion rates through high-energy electrons generated from synergistic interactions between plasma and the packing materials, where material defects provide abundant active sites to lower dissociation activation barriers. However, the energy efficiency of DBD systems remains constrained by the inherent limitations of electron-impact dissociation requiring excessive energy inputs (SEI > 100 eV per molecule) that poorly exploit injected power.
In GAD plasma systems, the gliding discharge channels allow sufficient contact with CO2 molecules, enabling effective dissociation at higher flow rates. Meanwhile, the increase of energy input becomes challenging (SEI < 10 eV per molecule) with elevated processing throughputs. Nonetheless, the energy efficiencies remain high due to the effective utilization of vibrational excitation mechanisms.
MWD plasma systems can effectively leverage electromagnetic energy from microwave power supplies to achieve high energy inputs (SEI > 100 eV per molecule), thus enabling promising CO2 conversion rates. But their energy efficiencies undergo a drastic decline under elevated pressure conditions (>100 Torr). This degradation stems from a mechanistic transition in CO2 dissociation pathways from vibrationally excited to electron-impact dissociation and thermally driven processes. Consequently, this pressure-dependent efficiency collapse confines their applicability to low-pressure operational regimes, limiting the practical applications.
The dimension of the simulation model often determines its proximity to the real multi-physics field. Commonly used simulation models in the literature can be classified into one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) models. As shown in Fig. 8, the 1D model simplifies the spatial distribution of multi-physics fields to a single degree of freedom, typically applied to uniform discharge modes or idealized simplified physical processes. For instance, Ponduri et al.185 equated the plasma region between symmetric parallel plate electrodes of DBD to a 1D model and studied the effect of flow and temperature distribution on CO2 dissociation conversion. In the case of GAD plasma, Wang et al.188 simplified the discharge channels within parallel plates and expanding electrodes, constructing two types of 1D models along the gas propagation direction to simulate the multi-physics field of GAD. MWD plasma can also be approximated by a 1D multi-physics model. Berthelot et al.189 simplified the energy injection in the MWD region of a cylindrical reactor to a power pulse and combined it with the temperature field along the axial direction to form a simplified 1D model.
The 2D multi-physics model is typically used in axisymmetric or rotationally symmetric discharge structures, as it can better fit the field distribution of actual discharge processes. For example, in a DBD reactor with dielectric material packing, a 1D model fails to embody the micro-discharges between dielectric particles. However, by modeling the basic unit consisting of electrodes, dielectric layers, filling particles, and gas gaps in 2D, the spatial distribution of the electric field can be well approximated.96 For the blade-electrode structure in GAD, a 2D model established for the plane where the GAD occurs, with appropriate initial and boundary conditions, can effectively represent both the discharge initiation area and the propagation region of the arc channel.190 In a spark discharge plasma reactor with rotating-symmetric needle-plate electrodes, modeling the vertical half cross-section in 2D allows for the physical field distribution of the entire space to be derived.71
When the model dimension is further increased to 3D, the accuracy of the multi-physics field evolution process in the spatial domain is further enhanced. Taking a packed bed DBD reactor as an example, Van Laer et al.174 established a 3D multi-physics field model of a reactor unit, which validated the electric field enhancement effect of the dielectric packing. For plasma reactors with more complex multi-physics processes, such as GAD137 and MWD,109 the 3D models developed are nearly identical to the actual reactor structures and can resolve details that are often overlooked by 1D or 2D models. Although the 3D modeling configurations of the GAD and MWD plasma systems bear certain resemblances to their respective reactor geometries, they differ significantly in the way their multi-physics parameters are prescribed. For the blade-structured GAD reactor illustrated, the electromagnetic boundary conditions closely mirror those used in the analogous 2D model. However, when coupling the electromagnetic field with the thermal and flow fields, the presence of the electrodes and the underlying tray, with each defined by its unique geometry and material properties, modulates heat and mass transfer. It further leads to spatial variations in temperature and gas velocity that, in turn, alter the evolution of the discharge channel. In contrast, the MWD plasma system, which operates in an electrodeless discharge mode, is typically simulated by focusing on the spatial propagation and distribution of the electromagnetic field. For a reactor of fixed dimensions, this yields a discharge region that remains essentially stationary. In the case of the converging-diverging nozzle shown, the tangential gas inlet induces a disturbed vortex flow, characterized by counter-rotating eddies that effectively mitigate thermal hot spots within the discharge zone.
Practically, the macroscopic physical laws fitted by the aforementioned multi-physics field models need to be further integrated with the microscopic reaction mechanisms calculated through 0D reaction kinetics, as discussed before. While the increased complexity of multi-physics field models brings them closer to actual physical processes, it also introduces challenges in coupling with 0D models. Currently, 1D models can be effectively coupled with 0D simulation calculations. However, in 2D models, the inclusion of reaction kinetics significantly increases the computational difficulty. As for 3D models, due to the excessive number of variables and constraints, it is still difficult to achieve a coupling with 0D models while maintaining the accuracy of the multi-physics field model.
In the development of plasma-assisted CO2 dissociation technologies for scale-up applications, multi-physics modeling serves as an indispensable analytical tool during the preliminary reactor design and operational monitoring to post-deployment maintenance phases. However, such modeling methodologies remain in their nascent stages of development.
The 0D kinetic simulations with temporal evolution of reaction mechanisms provide foundational insights into plasma-chemical pathways, which can be integrated with 1D multi-physics models for simplified reactor geometries under quasi-steady-state plasma discharges, enabling the incorporation of macroscopic transport phenomena (e.g., thermal gradients and species diffusion) into conversion performance analyses. Advancing to 2D modeling significantly enhances fidelity to physical realities through spatially resolved field distributions, albeit at the cost of exponentially increased computational complexity and convergence challenges. The coupling of 3D multi-physics simulations with 0D reaction kinetics currently faces prohibitive barriers due to incompatible temporal-spatial resolution requirements and unaccounted multi-scale plasma heterogeneities. The complexity of incorporating detailed reaction mechanisms into higher-dimensional models necessitates simplifications to make them computationally feasible.
Future progress demands dual advancements: (1) fundamental breakthroughs in plasma-specific reaction mechanisms through in situ diagnostics of critical parameters and (2) development of adaptive mesh algorithms capable of bridging disparate timescales. Such innovations will enable the transition from low-dimensional kinetic models to high-fidelity 3D simulations that holistically optimize CO2 dissociation systems across molecular to reactor scales.
The inherent small gas throughput presents a fundamental scalability challenge for current LTP reactors. Under laboratory conditions, DBD units process only 0.4–1.0 L min−1 CO2, GAD reactors handle around 15 L min−1, and MWD setups reach 15–30 L min−1. All of them are orders of magnitude below conventional chemical plant feed rates of hundreds of kg h−1, necessitating extensive parallelization.195 Modularization is thus the primary scale-up strategy, envisioning thousands of DBD tubes or hundreds of GAD units sharing manifolds, or parallel MWD reactors.191 Detz and van der Zwaan196 also concluded that no real economy of scale existed for MWD systems since the capital cost scales linearly with the total output capacity, favoring small to midsize modular plants. However, uniform power distribution, balancing between flow and temperature, and maintenance complexity remain significant hurdles.
From the perspective of economic viability, the cost per mole of CO produced is a key metric for energy and economic assessment. For DBD systems, a high SEI of about 30 MJ mol−1 translates to €0.4–0.7 per mol CO which counts for the operational expenditure (OPEX), with the levelized cost of electricity (LCOE) at €0.05–0.08 per kWh from wind and photovoltaic (PV) energy.197 Moreover, the capital expenditure (CAPEX) of a DBD system at a pilot scale of 144 kW was estimated at €1300 per kW, while the CAPEX estimated for a multi-tube DBD reactor at a flow rate of 100 m3 h−1 would increase linearly with tube numbers and the power consumed.198 Conversely, GAD systems achieve higher efficiencies. A recent TEA study by Osorio-Tejada et al.199 calculated a GAD system with a downstream carbon bed that consumed 19.5 GJ per ton (0.55 MJ mol−1) of the CO produced, corresponding to a production cost at $671 per ton, lower than the equivalent electrolysis of $962 per ton. Meanwhile, the CAPEX of pilot GAD systems was envisioned at about €2400 per kW, which can be applied to plants processing up to 1000 ton per year of CO2. As for MWD plasma systems, while they demonstrate significant advantages in energy efficiency (>50% at reduced pressure), the high CAPEX poses a potential threat to scalability. The work by O'Modhrain et al.126 concluded that a 20000 ton per year MWD plant would need around €7000 per kW with 3 MW generators.
On the other hand, comprehensive lifecycle assessment (LCA) addressing the holistic sustainability and environmental impacts of LTP technology provides essential metrics for evaluating the scale-up viability. A cradle-to-gate LCA comparing LTP-based and electrolytic CO production (each produced 100 ton per day of CO) found that the plasma route consumed 40% less total energy.200 While CO2 feedstock dominated impacts (60–80%), the plasma process showed environmental advantages over fossil-based CO in 7 out of 10 categories (e.g., 86% lower acidification, 91% lower freshwater ecotoxicity, and 83% lower fossil resource use) and a 7% lower global warming potential (GWP). Electrolysis, though offering moderate reductions in some categories, incurred an 8% higher GWP than the conventional processes.
The integration with renewable energy sources like solar and wind is a key potential advantage of LTP toward practical applications. LTP systems offer inherent modularity, enabling deployment matching local renewable capacity without massive single units. Utilizing low-marginal-cost renewable electricity significantly lowers OPEX. Against the backdrop of “net-zero” targets, anticipated rises in carbon taxes and higher carbon pricing in emissions trading systems (ETS) will further boost the competitiveness of renewable energy, thereby driving the deployment of LTP technology. For example, the economic competitiveness of MWD systems over conventional thermocatalytic processes using natural gas will occur by 2050 with a carbon price of €150 per t CO2.196 However, continuous CO2 feed requires buffering or storage if plasma operates only during surplus power, and downstream processing must accommodate intermittent input. While TEA suggests fuel costs could approach €2–3 per L under direct coupling with very cheap renewables, the low capacity factor risks capital underutilization. Realizing the economic potential hinges on continued research and development to improve efficiency, reduce capital and operational costs, develop robust catalysts and reactors, and lower renewable energy prices, alongside solutions for managing variable operation and product separation. Overall, while current costs are high, LTP-driven CO2 dissociation demonstrates potential techno-economic and environmental advantages over electrolytic or thermocatalytic routes, particularly when integrated with low-cost renewable electricity and targeted for specific modular applications.
At the macro level, we observed that the characteristics of different plasma power supplies lead to distinct plasma parameters. The control strategies for CO2 dissociation processes can mainly be divided into in situ and post-plasma approaches. For DBD plasma, in situ techniques focus on the use of various medium materials within the reactor. For GAD plasma, plasma regulation is achieved by introducing external magnetic fields or a vortex to alter the heat and mass transfer in the reaction zone. MWD plasma, on the other hand, can adjust the gas pressure and waveguide structure to modulate the distribution of microwave energy. In terms of post-plasma techniques, improving the conversion of high-throughput GAD plasma often requires the use of post-plasma catalysts or carbon beds. MWD plasma can optimize the temperature field by adding a nozzle and can also be coupled with electrocatalytic processes to enhance conversion. After comparing the CO2 conversion performance of different plasma forms, we observed that GAD and MWD plasmas have advantages in balancing processing throughput and energy consumption.
At the microscopic level, plasma parameters were obtained through in situ diagnostics using OES. Current research attempts to clarify the reaction mechanisms of CO2 dissociation in different types of plasmas from both 0D reaction kinetics and multi-physics simulation perspectives. For the representative plasma forms, namely DBD, GAD, and MWD plasma, the core CO2 dissociation pathways are centered around electronic excitation collision dissociation, vibrational excitation dissociation, and electronic or thermal effect-dominated dissociation, respectively. As for multi-physics simulations within LTP reactors, although 3D models offer the highest degree of accuracy in reflecting actual conditions, coupling with 0D models still requires simplification of the physical field spatial distribution using 1D or 2D models, through which the computational complexity can be reduced while maintaining reasonable accuracy.
From a detailed techno-economic perspective encompassing CO2 conversion performance, CO production costs, and benefits of renewable-energy integration, GAD and MWD plasma systems already demonstrate promising CO2 conversion rates and energy efficiencies. However, further reductions in capital and operational costs are essential during the modularization processes. Upon coupling with renewables, the economic competitiveness of LTP technology is projected to become increasingly prominent alongside global decarbonization efforts, with the potential to outperform mature CO2 valorization routes such as electrolysis and thermal catalysis.
Based on the key findings above, we believe that the LTP-based CO2 dissociation technology has the potential for broader adoption and mass-scale applications. As shown in Fig. 9, on the one hand, LTP technology needs specialized design in several areas, such as power supply improvements (compact modularization, adaptation to various load characteristics), catalytic coupling (development of plasma-specific catalysts, in situ/post-coupling strategies), and renewable energy systems (integrated system solutions). On the other hand, in the full-process LTP technology for valorization of CO2 dissociation, there are multiple layers of techno-economic constraints.201,202 The key breakthrough for scaling up the applications of plasma technology lies in the balance of various process indicators and improvement of its competitiveness. Upon considering the entire process flow of plasma-assisted CO2 dissociation technology, from the perspective of the upstream input, it necessitates the comprehensive evaluation of feed supply and separation, electricity, and construction/maintenance costs. At the midstream plasma conversion section, it is imperative to synergistically investigate discharge configurations and reactor geometries to elucidate reaction mechanisms while employing diverse characterization techniques to monitor real-time reaction states. This integrated approach aims to minimize energy loss during the conversion process. For downstream product management, it is vital to achieve optimal processing throughput while enhancing overall energy efficiency. Furthermore, the strategic development of application scenarios compatible with LTP technology should be prioritized. Concurrently, preparatory work for scale-up applications requires alignment with future policies that support CO2 valorization technologies and adaptation to renewable energy-driven electrification transition trends. This multi-scale optimization framework bridges fundamental plasma physics with practical engineering implementation requirements.
This journal is © The Royal Society of Chemistry 2025 |