Abstract: Warm plasma offers a promising route for CO2 splitting into valuable CO, yet recombination reactions of CO with oxygen, forming again CO2, have recently emerged as critical limitation. This study combines experiments and fluid dynamics + chemical kinetics modelling to comprehensively analyse the recombination reactions upon CO2 splitting in an atmospheric plasmatron. We introduce an innovative in-situ gas sampling technique, enabling 2D spatial mapping of gas product compositions and temperatures, experimentally confirming for the first time the substantial limiting effect of CO recombination reactions in the afterglow region. Our results show that the CO mole fraction at a 5 L/min flow rate drops significantly from 11.9 % at a vertical distance of z = 20 mm in the afterglow region to 8.6 % at z = 40 mm. We constructed a comprehensive 2D model that allows for spatial reaction rates analysis incorporating crucial reactions, and we validated it to kinetically elucidate this phenomenon. CO2 +M⇌O+CO+M and CO2 +O⇌CO+O2 are the dominant reactions, with the forward reactions prevailing in the plasma region and the backward reactions becoming prominent in the afterglow region. These results allow us to propose an afterglow quenching strategy for performance enhancement, which is further demonstrated through a meticulously developed plasmatron reactor with two-stage cooling. Our approach substantially increases the CO2 conversion (e.g., from 6.6 % to 19.5 % at 3 L/min flow rate) and energy efficiency (from 13.5 % to 28.5 %, again at 3 L/min) and significantly shortens the startup time (from ~ 150 s to 25 s). Our study underscores the critical role of inhibiting recombination reactions in plasma-based CO2 conversion and offers new avenues for performance enhancement.

Abstract: The generation of high-density plasmas on the surface of porous catalysts is very important for plasma catalysis, as it determines the active surface of the catalyst that is available for the reaction. In this work, we investigate the mechanism of surface and volume plasma streamer formation and propagation near micro-sized pores in dielectric barrier discharges operating in air at atmospheric pressure. A two-dimensional particle-in-cell/ Monte Carlo collision model is used to model the individual kinetic behavior of plasma species. Our calculations indicate that the surface discharge is enhanced on the surface of the catalyst pores compared with the microdischarge inside the catalyst pores. The reason is that the surface ionization wave induces surface charging along the catalyst pore sidewalls, leading to a strong electric field along the pore sidewalls, which in turn further enhances the surface discharge. Therefore, highly concentrated reactive species occur on the surfaces of the catalyst pores, indicating high-density plasmas on the surface of porous catalysts. Indeed, the maximum electron impact excitation and ionization rates occur on the pore surface, indicating the more pronounced production of excited state and electron-ion pairs on the pore surface than inside the pore, which may profoundly affect the plasma catalytic process. Published under license by AIP Publishing.

Abstract: While plasma treatment of skin diseases and wound healing has been proven highly effective, the underlying mechanisms, and more generally the effect of plasma radicals on skin tissue, are not yet completely understood. In this paper, we perform ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of plasma generated OH radicals with a model system composed of free fatty acids, ceramides, and cholesterol molecules. This model system is an approximation of the upper layer of the skin (stratum corneum). All interaction mechanisms observed in our simulations are initiated by H-abstraction from one of the ceramides. This reaction, in turn, often starts a cascade of other reactions, which eventually lead to the formation of aldehydes, the dissociation of ceramides or the elimination of formaldehyde, and thus eventually to the degradation of the skin barrier function.

Abstract: By means of steady-state experiments and a global model, we studied the effects of plasma-generated reactive species on the surface chemistry and coking in plasma-catalytic CH4/CO2 reforming at reduced pressure (8–40 kPa). We used a hybrid ZDPlasKin-CHEMKIN model to predict the species densities over time. The detailed plasma-catalytic mechanism consists of the plasma discharge scheme, a gas-phase chemistry set and a surface mechanism. Our experimental results show that the coupling of Ni/SiO2 catalyst with plasma is more effective in CH4/CO2 activation and conversion than unpacked DBD plasma, with syngas being the main products. The
highest total conversion of 16 % was achieved at 8000 V and 473 K, with corresponding CO and H2 yields of 15 % and 12 %, respectively. The reactants conversion and product selectivity are well captured by the kinetic model. Our simulation results suggest that vibrational species and radicals can accelerate the dissociative adsorption and Eley-Rideal (E-R) reactions. Path flux analysis shows that E-R reactions dominate the surface reaction pathways, which differs from thermal catalysis, indicating that the coupling of non-equilibrium plasma and catalysis can effectively shift the formation and consumption pathways of important adsorbates. For instance, our model suggests that HCOO(s) is primarily generated through the E-R reaction CO2(v) + H(s) → HCOO(s), while the hydrogenation reaction HCOO(s) + H → HCOOH(s) is the main source of HCOOH(s). Carbon deposition on the
catalyst surface is primarily formed through the stepwise dehydrogenation of CH4, while the E-R reactions enhanced by plasma-generated H and O atoms dominate the consumption of carbon deposition. This work provides new insights into the effects of reactive species on the surface chemistry in plasma-catalytic CH4/CO2 reforming.

Abstract: renewable energy. Plasma technology is promising for this purpose, as it can crack NH3 without the need for a catalyst and is highly compatible with renewable electricity, reducing the environmental footprint of the cracking process. This work investigates the NH3 cracking performance of four different warm plasma reactors with different configurations and operating in a wide range of conditions. We show that the NH3 conversion in warm plasma reactors is primarily determined by the specific energy input, with the main difference observed in the energy cost (EC) of cracking. The lowest EC obtained is 146 kJ/mol but at a conversion of only 8 %. A more reasonable conversion of around 50 % yields an EC of around 200 kJ/mol in two of the reactors investigated. Plasma reactors operating at higher feed flow rates are more efficient and yield a higher H2 production rate. Our data indicate that NH3 cracking in these warm plasma reactors occurs mainly via thermal chemistry, with nonthermal plasma chemistry playing a less prominent role. NH3 decomposes not only inside the plasma core but also in a hot volume around it, which reduces the EC. Our study shows that warm plasmas are significantly more efficient for NH3 cracking than cold plasmas, even when the latter are combined with catalysts.

Abstract: Microdischarges are an important type of plasma discharge that possess several unique characteristics, such as the presence of a stable glow discharge, high plasma density and intense excimer radiation, leading to several potential applications. The intense and controllable gas heating within the extremely small dimensions of microdischarges has been exploited in microthruster technologies by incorporating a micro-nozzle to generate the thrust. This kind of microthruster has a significantly improved specific impulse performance compared to conventional cold gas thrusters, and can meet the requirements arising from the emerging development and application of micro-spacecraft. In this paper, we performed a self-consistent 2D particle-in-cell simulation, with a Monte Carlo collision model, of a microdischarge operating in a prototype micro-plasma thruster with a hollow cylinder geometry and a divergent micro-nozzle. The model takes into account the thermionic electron emission including the Schottky effect, the secondary electron emission due to cathode bombardment by the plasma ions, several different collision processes, and a non-uniform argon background gas density in the cathode-anode gap. Results in the high-pressure (several hundreds of Torr), high-current (mA) operating regime showing the behavior of the plasma density, potential distribution, and energy flux towards the hollow cathode and anode are presented and discussed. In addition, the results of simulations showing the effect of different argon gas pressures, cathode material work function and discharge voltage on the operation of the microdischarge thruster are presented. Our calculated properties are compared with experimental data under similar conditions and qualitative and quantitative agreements are reached.

Abstract: Dielectric barrier discharges (DBDs) typically operate in the filamentary regime and thus exhibit great spatial and temporal non-uniformity. In order to optimize DBDs for various applications, such as in plasma catalysis, more fundamental insight is needed. Here, we consider how the millions of microdischarges, characteristic for a DBD, influence individual gas molecules. We use a Monte Carlo approach to determine the number of microdischarges to which a single molecule would be exposed, by means of particle tracing simulations through a full-scale packed bed DBD reactor, as well as an empty DBD reactor. We find that the fraction of microdischarges to which the molecules are exposed can be approximated as the microdischarge volume over the entire reactor gas volume. The use of this concept provides good agreement between a plasma-catalytic kinetics model and experiments for plasma-catalytic NH3 synthesis. We also show that the concept of the fraction of microdischarges indicates the efficiency by which the plasma power is transferred to the gas molecules. This generalised concept is also applicable for other spatially and temporally non-uniform plasmas.

Abstract: A study of ns-laser ablation is presented, which focuses on the transient behavior of the physical processes that act in and above a copper sample. A dimensionless multiphase collisional radiative model describes the interplay between the ablation, collisional, and radiative mechanisms. Calculations are done for a 6 ns-Nd:YAG laser pulse operating at 532 nm and fluences up to 15 J/cm2. Temporal intensity profiles as well as transmissivities are in good agreement with experimental results. It is found that volumetric ablation mechanisms and photo-processes both play an essential role in the onset of ns-laser induced breakdown.

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: Plasma modelling is an exciting subject in which virtually all physical disciplines are represented. Plasma models combine the electromagnetic, statistical and fluid dynamical theories that have their roots in the 19th century with the modern insights concerning the structure of matter that were developed throughout the 20th century. The present cluster issue consists of 20 invited contributions, which are representative of the state of the art in plasma modelling and numerical simulation. These contributions provide an in-depth discussion of the major theories and modelling and simulation strategies, and their applications to contemporary plasma-based technologies. In this editorial review, we introduce and complement those papers by providing a bird's eye perspective on plasma modelling and discussing the historical context in which it has surfaced.

Abstract: A modelling study is performed of high-density low-pressure inductively coupled CF4/CHF3/H2/Cl2/O2/HBr plasmas under different gas mixing ratios. A reaction set describing the complete plasma chemistry is presented and discussed. The gas fraction of each component in this mixture is varied to investigate the sensitivity of the plasma properties, like electron density, plasma potential and species densities, towards the gas mixing ratios. This research is of great interest for microelectronics applications because these gases are often combined in two (or more)-component mixtures, and mixing gases or changing the fraction of a gas can sometimes yield unwanted reaction products or unexpected changes in the overall plasma properties due to the increased chemical complexity of the system. Increasing the CF4 fraction produces more F atoms for chemical etching as expected, but also more prominently lowers the density of Cl atoms, resulting in an actual drop in the etch rate under certain conditions. Furthermore, CF4 decreases the free electron density when mixed with Cl2. However, depending on the other gas components, CF4 gas can also sometimes enhance free electron density. This is the case when HBr is added to the mixture. The addition of H2 to the gas mixture will lower the sputtering process, not only due to the lower overall positive ion density at higher H2 fractions, but also because more H+, H2 + and H3 + are present and they have very low sputter yields. In contrast, a larger Cl2 fraction results in more chemical etching but also in less physical sputtering due to a smaller abundance of positive ions. Increasing the O2 fraction in the plasma will always lower the etch rate due to more oxidation of the wafer surface and due to a lower plasma density. However, it is also observed that the density of F atoms can actually increase with rising O2 gas fraction. This is relevant to note because the exact balance between fluorination and oxidation is important for fine-tuning the overall etch rate and for control of the sidewall profile. Finally, HBr is often used as a chemical etcher, but when mixed with F- or Cl-containing gases, HBr creates the same diluting effects as Ar or He, because a
higher fraction results in less chemical etching but more (physical) sputtering.

Abstract: The lowest-energy configurations of Cn(n≤55) clusters are obtained using the energy-minimization technique with the conjugate gradient method where a modified Brenner potential is invoked to describe the carbon and hydrocarbon interaction. We found that the ground-state configuration consists of a single ring for small number of C atoms and multiring structures are found with increasing n, which can be in planar, bowl-like or caplike form. Contrary to previous predictions, the binding energy Eb does not show even-odd oscillations and only small jumps are found in the Eb(n) curve as a consequence of specific types of edges or equivalently the number of secondary atoms. We found that hydrogenation of the edge atoms may change the ground-state configuration of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to trigonal and hexagonal shaped carbon clusters and to clusters having a graphenelike configuration. Trigonal clusters are never the ground state while hexagonal-shaped clusters are only the ground state when they have zigzag edges.