Abstract: A one-dimensional, implicit particle-in-cell Monte Carlo collision model is used to simulate the plasma kinetic properties at a steady state in a parallel-plate direct current argon glow microdischarge under various operating conditions, such as driving voltage (301000 V) and gap size (101000 µm) at atmospheric pressure. First, a comparison between rf and dc modes is shown for the same pressure, driving voltage and gap spacing. Furthermore, the effect of gap size scaling (in the range of 101000 µm) on the breakdown voltage, peak electron density and peak electron current density at the breakdown voltage is examined. The breakdown voltage is lower than 150 V in all gaps considered. The microdischarge is found to have a neutral bulk plasma region and a cathode sheath region with size varying with the applied voltage and the discharge gap. In our calculations, the electron and ion densities are of the order of 10181023 m−3, which is in the glow discharge limit, as the ionization degree is lower than 1% . The electron energy distribution function shows a two-energy group distribution at a gap of 10 µm and a three-energy group distribution at larger gaps such as 200 µm and 1000 µm, emphasizing the importance of the gap spacing in dc microdischarges.

Abstract: The global relaxation process after pulsed laser induced photodetachment in a rf electronegative SIH4 discharge is studied by a self-consistent kinetic one-dimensional particle-in-cell-Monte Carlo model. Our results reveal a comprehensive physical picture of the relaxation process, including the main plasma variables, after a perturbation up to the full recovery of the steady state. A strong influence of the photodetachment on the discharge is found, which results from an increase of the electron density, leading to a weaker bulk field, and hence to a drop in the high energy tail of the electron energy distribution function (EEDF), a reduction of the reaction rates of electron impact attachment and ionization, and a subsequent decrease of the positive and negative ion densities. All the plasma quantities related to electrons recover synchronously. The recovery time of the ion densities is about 1-2 orders of magnitude longer than that of the electrons due to different recovery mechanisms. The modeled behavior of all the charged particles agrees very well with experimental results from the literature. In addition, our work clarifies some unclear processes assumed in the literature, such as the relaxation of the EEDF, the evolution of the electric field, and the recovery of negative ions.

Abstract: We have computationally investigated the introduction of copper elemental particles in an inductively coupled plasma torch connected to a sampling cone, including for the first time the ionization of the sample. The sample is inserted as liquid particles, which are followed inside the entire torch, i.e., from the injector inlet up to the ionization and reaching the sampler. The spatial position of the ion clouds inside the torch as well as detailed information on the copper species fluxes at the position of the sampler orifice and the exhausts of the torch are provided. The effect of on- and off-axis injection is studied. We clearly show that the ion clouds of on-axis injected material are located closer to the sampler with less radial diffusion. This guarantees a higher transport efficiency through the sampler cone. Moreover, our model reveals the optimum ranges of applied power and flow rates, which ensure the proper position of ion clouds inside the torch, i.e., close enough to the sampler to increase the fraction that can enter the mass spectrometer and with minimum loss of material toward the exhausts as well as a sufficiently high plasma temperature for efficient ionization.

Abstract: We have developed a comprehensive kinetic model to study the O atom kinetics in an O₂plasma and its afterglow. By adopting a pseudo-1D plug-flow formalism within the kinetic model, our aim is to assess how far the O atoms travel in the plasma afterglow, evaluating its potential as a source of O atoms for post-plasma gas conversion applications. Since we could not find experimental data for pure O₂plasma at atmospheric pressure, we first validated our model at low pressure (1–10 Torr) where very good experimental data are available. Good agreement between our model and experiments was achieved for the reduced electric field, gas temperature and the densities of the dominant neutral species, i.e. O₂(a), O₂(b) and O. Subsequently, we confirmed that the chemistry set is consistent with thermodynamic equilibrium calculations at atmospheric pressure. Finally, we investigated the O atom densities in the O₂plasma and its afterglow, for which we considered a microwave O₂plasma torch, operating at a pressure between 0.1 and 1 atm, for a flow rate of 20 slm and an specific energy input of 1656 kJ mol⁻¹. Our results show that for both pressure conditions, a high dissociation degree of ca. 92% is reached within the discharge. However, the O atoms travel much further in the plasma afterglow for<italic>p</italic>= 0.1 atm (9.7 cm) than for<italic>p</italic>= 1 atm (1.4 cm), attributed to the longer lifetime (3.8 ms at 0.1 atm vs 1.8 ms at 1 atm) resulting from slower three-body recombination kinetics, as well as a higher volumetric flow rate.

Abstract: A parameter study is carried out for an inductively coupled plasma used for the synthesis of carbon nanotubes or carbon nanofibres (CNTs/CNFs), by means of the Hybrid Plasma Equipment Model. The influence of processing parameters including gas ratio for four different gas mixtures typically used for CNT/CNF growth (i.e. CH4/H2, CH4/NH3, C2H2/H2 and C2H2/NH3), inductively coupled plasma (ICP) power (501000 W), operating pressure (10 mTorr1 Torr), bias power (01000 W) and temperature of the substrate (01000 °C) on the plasma chemistry is investigated and the optimized conditions for CNT/CNF growth are analysed. Summarized, our calculations suggest that a lower fraction of hydrocarbon gases (CH4 or C2H2, i.e. below 20%) and hence a higher fraction of etchant gases (H2 or NH3) in the gas mixture result in more 'clean' conditions for controlled CNT/CNF growth. The same applies to a higher ICP power, a moderate ICP gas pressure above 100 mTorr (at least for single-walled carbon nanotubes), a high bias power (for aligned CNTs) and an intermediate substrate temperature.

Abstract: Methane reforming by plasma catalysis is a complex process that is far from understood. It requires a multidisciplinary approach which ideally takes into account all effects from the plasma on the catalyst, and vice versa. In this contribution, we focus on the interactions of CHx (x = {1,2,3}) radicals that are created in the plasma with several nickel catalyst surfaces. To this end, we perform reactive molecular dynamics simulations, based on the ReaxFF potential, in a wide temperature range of 4001600 K. First, we focus on the H2 formation as a function of temperature and surface structure. We observe that substantial H2 formation is obtained at 1400 K and above, while the role of the surface structure seems limited. Indeed, in the initial stage, the type of nickel surface influences the CH bond breaking efficiency of adsorbed radicals; however, the continuous carbon diffusion into the surface gradually diminishes the surface crystallinity and therefore reduces the effect of surface structure on the H2 formation probability. Furthermore, we have also investigated to what extent the species adsorbed on the catalyst surface can participate in surface reactions more in general, for the various surface structures and as a function of temperature. These results are part of the ongoing research on the methane reforming by plasma catalysis, a highly interesting yet complex alternative to conventional reforming processes.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

Abstract: Vibrational excitation represents an efficient channel to drive the dissociation of CO₂in a non-thermal plasma. Its viability is investigated in low-pressure pulsed discharges, with the intention of selectively exciting the asymmetric stretching mode, leading to stepwise excitation up to the dissociation limit of the molecule. Gas heating is crucial for the attainability of this process, since the efficiency of vibration–translation (V–T) relaxation strongly depends on temperature, creating a feedback mechanism that can ultimately thermalize the discharge. Indeed, recent experiments demonstrated that the timeframe of V–T non-equilibrium is limited to a few milliseconds at ca. 6 mbar, and shrinks to the<italic>μ</italic>s-scale at 100 mbar. With the aim of backtracking the origin of gas heating in pure CO₂plasma, we perform a kinetic study to describe the energy transfers under typical non-thermal plasma conditions. The validation of our kinetic scheme with pulsed glow discharge experiments enables to depict the gas heating dynamics. In particular, we pinpoint the role of vibration–vibration–translation relaxation in redistributing the energy from asymmetric to symmetric levels of CO₂, and the importance of collisional quenching of CO₂electronic states in triggering the heating feedback mechanism in the sub-millisecond scale. This latter finding represents a novelty for the modelling of low-pressure pulsed discharges and we suggest that more attention should be paid to it in future studies. Additionally, O atoms convert vibrational energy into heat, speeding up the feedback loop. The efficiency of these heating pathways, even at relatively low gas temperature and pressure, underpins the lifetime of V–T non-equilibrium and suggests a redefinition of the optimal conditions to exploit the ‘ladder-climbing’ mechanism in CO₂discharges.

Abstract: Carbon nanotubes (CNTs) are nowadays routinely grown in a thermal CVD setup. State-of-the-art plasma-enhanced CVD (PECVD) growth, however, offers advantages over thermal CVD. A lower growth temperature and the growth of aligned freestanding single-walled CNTs (SWNTs) makes the technique very attractive. The atomic scale growth mechanisms of PECVD CNT growth, however, remain currently entirely unexplored. In this contribution, we employed molecular dynamics simulations to focus on the effect of applying an electric field on the SWNT growth process, as one of the effects coming into play in PECVD. Using sufficiently strong fields results in (a) alignment of the growing SWNTs, (b) a better ordering of the carbon network, and (c) a higher growth rate relative to thermal growth rate. We suggest that these effects are due to the small charge transfer occurring in the Ni/C system. These simulations constitute the first study of PECVD growth of SWNTs on the atomic level.

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: Plasma catalysis is gaining increasing interest for various environmental applications, but the crucial question is whether plasma can be created inside catalyst pores and under which conditions. In practice, various catalytic support materials are used, with various dielectric constants. We investigate here the influence of the dielectric constant on the plasma properties inside catalyst pores and in the sheath in front of the pores, for various pore sizes. The calculations are performed by a two-dimensional fluid model for an atmospheric pressure dielectric barrier discharge in helium. The electron impact ionization rate, electron temperature, electron and ion density, as well as the potential distribution and surface charge density, are analyzed for a better understanding of the discharge behavior inside catalyst pores. The results indicate that, in a 100 μm pore, the electron impact ionization in the pore, which is characteristic for the plasma generation inside the pore, is greatly enhanced for dielectric constants below 300. Smaller pore sizes only yield enhanced ionization for smaller dielectric constants, i.e., up to εr = 200, 150, and 50 for pore sizes of 50, 30, and 10 μm. Thus, the most common catalyst supports, i.e., Al2O3 and SiO2, which have dielectric constants around εr = 8−11 and 4.2, respectively, should allow more easily that microdischarges can be formed inside catalyst pores, even for smaller pore sizes. On the other hand, ferroelectric materials with dielectric constants above 300 never seem to yield plasma enhancement inside catalyst pores, not even for 100 μm pore sizes. Furthermore, it is clear that the dielectric constant of the material has a large effect on the extent of plasma enhancement inside the catalyst pores, especially in the range between εr = 4 and εr = 200. The obtained results are explained in detail based on the surface charge density at the pore walls,
and the potential distribution and electron temperature inside and above the pores. The results obtained with this model are
important for plasma catalysis, as the production plasma species in catalyst pores might affect the catalyst properties, and thus
improve the applications of plasma catalysis.