Abstract: CO2 conversion into value-added chemicals is gaining increasing interest in recent years, and a gliding arc plasma has great potential for this purpose, because of its high energy efficiency. In this study, a chemical reaction kinetics model is presented to study the CO2 splitting in a gliding arc discharge. The calculated
conversion and energy efficiency are in good agreement with experimental data in a range of different operating conditions. Therefore, this reaction kinetics model can be used to elucidate the dominant chemical reactions contributing to CO2 destruction and formation. Based on this reaction pathway analysis, the restricting factors for CO2 conversion are figured out, i.e., the reverse reactions and the small treated gas fraction. This allows us to propose some solutions in order to improve the CO2 conversion, such as decreasing the gas temperature, by using a high frequency discharge, or increasing the power
density, by using a micro-scale gliding arc reactor, or by removing the reverse reactions, which could be realized in practice by adding possible scavengers for O atoms, such as CH4. Finally, we compare our results with other types of plasmas in terms of conversion and energy efficiency, and the results illustrate that gliding arc discharges are indeed quite promising for CO2 conversion, certainly when keeping in mind the possible solutions for further performance improvement.

Abstract: Microwave (MW) plasmas represent a promising solution for efficient CO2 dissociation. MW discharges are also very versatile and can be sustained at various pressure and gas flow regimes. To identify the most favorable conditions for the further scale-up of the CO2 decomposition reaction, a MW plasma reactor operating in pure CO2 in a wide pressure range (200 mbar–1 bar) is studied. Three different gas flow configurations are explored: a direct, reverse and a vortex regime. The CO2 conversion and energy efficiency drop almost linearly with increasing pressure, regardless of the gas flow regime. The results obtained in the direct flow configuration underline the importance of post-discharge cooling, as the exhaust of the MW plasma reactor in this regime expanded into the vacuum chamber without additional quenching. As a result, this system yields exhaust temperatures of up to 1000 K, which explains the lowest conversion (∼3.5% at 200 mbar and 2% at 1 bar). A post-discharge cooling step is introduced for the reverse gas inlet regime and allows the highest conversion to be achieved (∼38% at 200 mbar and 6.2% at 1 bar, with energy efficiencies of 23% and 3.7%). Finally, a tangential gas inlet is utilized in the vortex configuration to generate a swirl flow pattern. This results in the generation of a stable discharge in a broader range of CO2 flows (15–30 SLM) and the highest energy efficiencies obtained in this study (∼25% at 300 mbar and ∼13% at 1 bar, at conversions of 21% and 12%). The experimental results are complemented with computational fluid dynamics simulations and with the analysis of the latest literature to identify the further research directions.

Abstract: Atmospheric pressure gliding arc (GA) discharges are gaining increasing interest for CO2 conversion and other gas conversion applications, due to their simplicity and high energy efficiency. However, they are characterized by some drawbacks, such as non-uniform gas treatment, limiting the conversion, as well as the development of a hot cathode spot, resulting in severe electrode degradation. In this work, we built a dual-vortex plasmatron, which is a GA plasma reactor with innovative electrode configuration, to solve the above problems. The design aims to improve the CO2 conversion capability of the GA reactor by elongating the arc in two directions, to increase the residence time of the gas inside the arc, and to actively cool the cathode spot by rotation of the arc and gas convection. The measured CO2 conversion and corresponding energy efficiency indeed look very promising. In addition, we developed a fluid dynamics non-thermal plasma model with argon chemistry, to study the arc behavior in the reactor and to explain the experimental results.

Abstract: Plasma catalysis, a promising technology for conversion of CO2 into value-added chemicals near room temperature, is gaining increasing interest. A dielectric barrier discharge (DBD) plasma has attracted attention due to its simple design and operation at near ambient conditions, ease to implement catalysts in the plasma zone and upscaling ability to industrial applications. To improve its main drawbacks, being relatively low conversion and energy efficiency, a packing material is used in the plasma discharge zone of the reactor, sometimes decorated by a catalytic material. Nevertheless, the extent to which different properties of the packing material influence plasma performance is still largely unexplored and unknown. In this study, the particular effect of synthesis induced differences in the morphology of a TiO2 shell covering a SiO2 core packing material on the plasma conversion of CO2 is studied. TiO2 has been successfully deposited around 1.6–1.8 mm sized SiO2 spheres by means of spray coating, starting from aqueous citratoperoxotitanate(IV) precursors. Parameters such as concentration of the Ti(IV) precursor solutions and addition of a binder were found to affect the shells’ properties and surface morphology and to have a major impact on the CO2 conversion in a packed bed DBD plasma reactor. Core-shell SiO2@TiO2 obtained from 0.25 M citratoperoxotitante(IV) precursors with the addition of a LUDOX binder showed the highest CO2 conversion 37.7% (at a space time of 70 s corresponding to an energy efficiency of 2%) and the highest energy efficiency of 4.8% (at a space time of 2.5 s corresponding to a conversion of 3%).

Abstract: We present a confined atmospheric pressure glow discharge plasma reactor, with very good performance towards dry reforming of methane, i.e., CO2 and CH4 conversion of 64 % and 94 %, respectively, at an energy cost of 3.5–4 eV/molecule (or 14–16 kJ/L). This excellent performance is among the best reported up to now for all types of plasma reactors in literature, and is due to the confinement of the plasma, which maximizes the fraction of gas passing through the active plasma region. The main product formed is syngas, with H2O and C2H2 as byproducts. We developed a quasi-1D chemical kinetics model, showing good agreement with the experimental results, which provides a thorough insight in the reaction pathways underlying the conversion of CO2 and CH4 and the formation of the different products.

Abstract: In this work, we take a crucial step towards the industrial readiness of plasma-based CO2 conversion. We present a stepwise method to study plasma reactors in series as a first approach to a recycle flow. By means of this procedure, the CO2 conversion is enhanced by a factor of 3, demonstrating that a single-pass plasma treatment performs far below the optimal capacity of the reactor. Furthermore, we explore the effect of O2 in the mixture with our flexible procedure. Addition of O2 in the mixture has a clear detrimental effect on the conversion, in agreement with other experiments in atmospheric pressure plasmas. O2 removal is however highly beneficial, demonstrating a conversion per pass that is 1.6 times higher than the standard procedure. Indeed, extracting one of the products prevents recombination reactions. Based on these insights, we discuss opportunities for further improvements, especially in the field of specialised separation techniques.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

Abstract: A hybrid model is used to investigate the fragmentation of C4F8 inductive discharges. Indeed, the resulting reactive species are crucial for the optimization of the Si-based etching process, since they determine the mechanisms of fluorination, polymerization, and sputtering. In this paper, we present the dissociation degree, the density ratio of F vs. CxFy (i.e., fluorocarbon (fc) neutrals), the neutral vs. positive ion density ratio, details on the neutral and ion components, and fractions of various fc neutrals (or ions) in the total fc neutral (or ion) density in a C4F8 inductively coupled plasma source, as well as the effect of pressure and power on these results. To analyze the fragmentation behavior, the electron density and temperature and electron energy probability function (EEPF) are investigated. Moreover, the main electron-impact generation sources for all considered neutrals and ions are determined from the complicated C4F8 reaction set used in the model. The C4F8 plasma fragmentation is explained, taking into account many factors, such as the EEPF characteristics, the dominance of primary and secondary processes, and the thresholds of dissociation and ionization. The simulation results are compared with experiments from literature, and reasonable agreement is obtained. Some discrepancies are observed, which can probably be attributed to the simplified polymer surface kinetics assumed in the model.