Abstract: The electrical asymmetry effect is realized by applying multiple frequency power sources
(13.56 MHz and 27.12 MHz) to a capacitively biased substrate electrode in a specific inductively
coupled plasma reactor. On the one hand, by adjusting the phase angle θ between the multiple
frequency power sources, an almost linear self-bias develops on the substrate electrode, and
consequently the ion energy can be well modulated, while the ion flux stays constant within a
large range of θ. On the other hand, the plasma density and ion flux can be significantly
modulated by tuning the inductive power supply, while only inducing a small change in the self-
bias. Independent control of self-bias/ion energy and ion flux can thus be realized in this specific
inductively coupled plasma reactor.

Abstract: In this work we study energy transport in a gliding arc discharge with two diverging flat
electrodes in argon gas at atmospheric pressure. The discharge is ignited at the shortest electrode
gap and it is pushed downstream by a forced gas flow. The current values considered are
relatively low and therefore a non-equilibrium plasma is produced. We consider two cases, i.e.
with high and low discharge current (28 mA and 2.8mA), and a constant gas flow of 10 lmin −1 ,
with a significant turbulent component to the velocity. The study presents an analysis of the
various energy transport mechanisms responsible for the redistribution of Joule heating to the
plasma species and the moving background gas. The objective of this work is to provide a
general understanding of the role of the different energy transport mechanisms in arc formation
and sustainment, which can be used to improve existing or new discharge designs. The work is
based on a three-dimensional numerical model, combining a fluid plasma model, the shear stress
transport Reynolds averaged Navier–Stokes turbulent gas flow model, and a model for gas
thermal balance. The obtained results show that at higher current the discharge is constricted
within a thin plasma column several hundred kelvin above room temperature, while in the low-
current discharge the combination of intense convective cooling and low Joule heating prevents
discharge contraction and the plasma column evolves to a static non-moving diffusive plasma,
continuously cooled by the flowing gas. As a result, the energy transport in the two cases is
determined by different mechanisms. At higher current and a constricted plasma column, the
plasma column is cooled mainly by turbulent transport, while at low current and an unconstricted
plasma, the major cooling mechanism is energy transport due to non-turbulent gas convection. In
general, the study also demonstrates the importance of turbulent energy transport in
redistributing the Joule heating in the arc and its significant role in arc cooling and the formation
of the gas temperature profile. In general, the turbulent energy transport lowers the average gas
temperature in the arc, thus allowing additional control of thermal non-equilibrium in the
discharge.

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

Abstract: We present a zero-dimensional plasma kinetics model, including both surface and gas phase kinetics, to determine the role of vibrationally excited states in plasma-catalytic ammonia synthesis. We defined a new method to systematically capture the conditions of dielectric barrier discharges (DBDs), including those found in packed bed DBDs. We included the spatial and temporal nature of such discharges by special consideration of the number of micro-discharges in the model. We introduce a parameter that assigns only a part of the plasma power to the microdischarges, to scale the model conditions from filamentary to uniform plasma. Because of the spatial and temporal behaviour of the micro-discharges, not all micro-discharges occurring in the plasma reactor during a certain gas residence time are affecting the molecules. The fraction of power considered in the model ranges from 0.005 %, for filamentary plasma, to 100 %, for uniform plasma. If vibrational excitation is included in the plasma chemistry, these different conditions, however, yield an ammonia density that is only varying within one order of magnitude. At only 0.05 % of the power put into the uniform plasma component, a model neglecting vibrational excitation clearly does not result in adequate amounts of ammonia. Thus, our new model, which accounts for the concept in which not all the power is deposited by the micro-discharges, but some part may also be distributed in between them, suggests that vibrational kinetic processes are really important in (packed bed) DBDs. Indeed, vibrational excitation takes place in both the uniform plasma between the micro-discharges and in the strong micro-discharges, and is responsible for an increased N2 dissociation rate. This is shown here for plasma-catalytic ammonia synthesis, but might also be valid for other gas conversion applications.

Abstract: Plasma has been studied for several years to convert CO2 into value-added products. If CO2 could be converted in the presence of H2O as a cheap H-source for making syngas and oxygenates, it would mimic natural photosynthesis. However, CO2/H2O plasmas have not yet been extensively studied, not by experiments, and certainly not computationally. Therefore, we present here a kinetic modelling study to obtain a greater understanding of the vibrational kinetics of a CO2/H2O microwave plasma. For this purpose, we first created an electron impact cross section set for H2O, using a swarm-derived method. We added the new cross section set and CO2/H2O-related chemistry to a pure CO2 model. While it was expected that H2O addition mainly causes quenching of the CO2 asymmetric mode vibrational levels due to the additional CO2/H2O vibrational-translational relaxation, our model shows that the modifications in the vibrational kinetics are mainly induced by the strong electron dissociative attachment to H2O molecules, causing a reduction in electron density, and the corresponding changes in the input of energy into the CO2 vibrational levels by electron impact processes.

Abstract: Cold atmospheric plasma (CAP) and plasma-treated liquids (PTLs) have recently become a promising option for cancer treatment, but the underlying mechanisms of the anti-cancer effect are still to a large extent unknown. Although hydrogen peroxide () has been recognized as the major anti-cancer agent of PTL and may enable selectivity in a certain concentration regime, the co-existence of nitrite can create a synergistic effect. We develop a mathematical model to describe the key species and features of the cellular response toward PTL. From the numerical solutions, we define a number of dependent variables, which represent feasible measures to quantify cell susceptibility in terms of the membrane diffusion rate constant and the intracellular catalase concentration. For each of these dependent variables, we investigate the regimes of selective versus non-selective, and of synergistic versus non-synergistic effect to evaluate their potential role as a measure of cell susceptibility. Our results suggest that the maximal intracellular concentration, which in the selective regime is almost four times greater for the most susceptible cells compared to the most resistant cells, could be used to quantify the cell susceptibility toward exogenous . We believe our theoretical approach brings novelty to the field of plasma oncology, and more broadly, to the field of redox biology, by proposing new ways to quantify the selective and synergistic anti-cancer effect of PTL in terms of inherent cell features.

Abstract: A pulsed microwave surfaguide discharge operating at 2.45 GHz was used for the conversion of molecular nitrogen into valuable compounds in several gas mixtures: N2 :O2 , N2 :O2 :CO2 and N2 :CO2 . The ro-vibrational absorption bands of the molecular species were monitored by a Fourier transform infrared apparatus in the post-discharge region in order to evaluate the relative number density of species, specifically NO production. The effects of specific energy input, pulse frequency, gas flow fraction, gas admixture and gas flow rate were studied for better understanding and optimization of the NO production yield and the corresponding energy cost (EC). By both the experiment and modelling, a highest NO yield is obtained at N2 :O2 (1:1) gas ratio in N2 :O2 mixture. The NO yield reveals a small growth followed by saturation when pulse repetition frequency increases. The energy efficiency start decreasing after the energy input reaches about 5 eV/molec, whereas the NO yield rises steadily at the same time. The lowest EC of about 8 MJ mol−1 corresponding to the yield and the energy efficiency of about 7% and 1% are found, respectively, in an optimum discharge condition in our case.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is a challenging neoplastic disease, mainly due to the development of resistance to radio- and chemotherapy. Cold atmospheric plasma (CAP) is an alternative technology that can eliminate cancer cells through oxidative damage, as shown in vitro, in ovo, and in vivo. However, how CAP affects the pancreatic stellate cells (PSCs), key players in the invasion and metastasis of PDAC, is poorly understood. This study aims to determine the effect of an anti-PDAC CAP treatment on PSCs tissue developed in ovo using mono- and co-cultures of RLT-PSC (PSCs) and Mia PaCa-2 cells (PDAC). We measured tissue reduction upon CAP treatment and mRNA expression of PSC activation markers and extracellular matrix (ECM) remodelling factors via qRT-PCR. Protein expression of selected markers was confirmed via immunohistochemistry. CAP inhibited growth in Mia PaCa-2 and co-cultured tissue, but its effectiveness was reduced in the latter, which correlates with reduced ki67 levels. CAP did not alter the mRNA expression of PSC activation and ECM remodelling markers. No changes in MMP2 and MMP9 expression were observed in RLT-PSCs, but small changes were observed in Mia PaCa-2 cells. Our findings support the ability of CAP to eliminate PDAC cells, without altering the PSCs.

Abstract: Plasma catalysis is gaining increasing interest for various applications, but the underlying mechanisms are still far from understood. Hence, more fundamental research is needed to understand these mechanisms. This can be obtained by both modelling and experiments. This foundations paper describes the fundamental insights in plasma catalysis, as well as efforts to gain more insights by modelling and experiments. Furthermore, it discusses the state-of-the-art of the major plasma catalysis applications, as well as successes and challenges of technology transfer of these applications.

Abstract: The Quantemol database (QDB) provides cross sections and rates of processes important for plasma models; heavy particle collisions (chemical reactions) and electron collision processes are considered. The current version of QDB has data on 28 917 processes between 2485 distinct species plus data for surface processes. These data are available via a web interface or can be delivered directly to plasma models using an application program interface; data are available in formats suitable for direct input into a variety of popular plasma modeling codes including HPEM, COMSOL, ChemKIN, CFD-ACE+, and VisGlow. QDB provides ready assembled plasma chemistries plus the ability to build bespoke chemistries. The database also provides a Boltzmann solver for electron dynamics and a zero-dimensional model. Thesedevelopments, use cases involving O₂, Ar/NF₃, Ar/NF₃/O₂, and He/H₂O/O₂chemistries, and plans for the future are presented.

Abstract: Utilizing ammonia as a storage medium for hydrogen is currently receiving increased attention. A possible method to retrieve the hydrogen is by plasma-catalytic decomposition. In this work, we combined an experimental study, using a dielectric barrier discharge plasma reactor, with a plasma kinetic model, to get insights into the decomposition mechanism. The experimental results revealed a similar effect on the ammonia conversion when changing the flow rate and power, where increasing the specific energy input (higher power or lower flow rate) gave an increased conversion. A conversion as high as 82 % was achieved at a specific energy input of 18 kJ/Nl. Furthermore, when changing the discharge volume from 31 to 10 cm3, a change in the plasma distribution factor from 0.2 to 0.1 was needed in the model to best describe the conversions of the experiments. This means that a smaller plasma volume caused a higher transfer of energy through micro-discharges (non-uniform plasma), which was found to promote the decomposition of ammonia. These results indicate that it is the collisions between NH3 and the high-energy electrons that initiate the decomposition. Moreover, the rate of ammonia destruction was found by the model to be in the order of 1022 molecules/(cm3 s) during the micro-discharges, which is 5 to 6 orders of magnitude higher than in the afterglows. A considerable re-formation of ammonia was found to take place in the afterglows, limiting the overall conversion. In addition, the model revealed that implementation of packing material in the plasma introduced high concentrations of surface-bound hydrogen atoms, which introduced an additional ammonia re-formation pathway through an Eley-Rideal reaction with gas phase NH2. Furthermore, a more uniform plasma is predicted in the presence of MgAl2O4, which leads to a lower average electron energy during micro-discharges and a lower conversion (37 %) at a comparable residence time for the plasma alone (51 %).

Abstract: Experimental evidence in the literature has shown that low-current direct current nitrogen discharges can exist in both glow and arc regimes at atmospheric pressure. However, modelling investigations of the positive column that include the influence of the cathode phenomena are scarce. In this work we developed a 2D axisymmetric model of a plasma discharge in flowing nitrogen gas, studying the influence of the two cathode emission mechanisms—thermionic field emission and secondary electron emission—on the cathode region and the positive column. We show for an inlet gas flow velocity of 1 m s⁻¹in the current range of 80–160 mA, that the electron emission mechanism from the cathode greatly affects the size and temperature of the cathode region, but does not significantly influence the discharge column at atmospheric pressure. We also demonstrate that in the discharge column the electron density balance is local and the electron production and destruction is dominated by volume processes. With increasing flow velocity, the discharge contraction is enhanced due to the increased convective heat loss. The cross sectional area of the conductive region is strongly dependent on the gas velocity and heat conductivity of the gas.

Abstract: The recent energy crisis revealed that there is a strong need to replace hydrocarbon-fueled industrial nitrogen fixation processes by alternative, more sustainable methods. In light of this, plasma-based nitrogen fixation remains one of the most promising options, considering both theoretical and experimental aspects. Lately, plasma interacting with water has received considerable attention in nitrogen fixation applications as it can trigger a unique gas- and liquid-phase chemistry. Within this context, a critical exploration of plasma-assisted nitrogen fixation with or without water presence is of great interest with an emphasis on energy costs, particularly in plasma reactors which have potential for large-scale industrial application. In this work, the presence of water in a multi-pin plasma system on nitrogen oxidation is experimentally investigated by comparing two pulsed negative DC voltage plasmas in metal–metal and metal–liquid electrode configurations. The plasma setups are designed to create similar plasma properties, including plasma power and discharge regime in both configurations. The system energy cost is calculated, considering nitrogen-containing species generated in gas and liquid phases as measured by a gas analyzer, nitrate sensor, and a colorimetry method. The energy cost profile as a function of specific energy input showed a strong dependency on the plasma operational frequency and the gas flow rate, as a result of different plasma operation regimes and initiated reverse processes. More importantly, the presence of the plasma/liquid interface increased the energy cost up to 14 ± 8%. Overall, the results showed that the presence of water in the reaction zone has a negative impact on the nitrogen fixation process.

Abstract: The use of microfluidic devices in the field of plasma-liquid interaction can unlock unique possibilities to investigate the effects of plasma-generated reactive species for environmental and biomedical applications. So far, very little simulation work has been performed on microfluidic devices in contact with a plasma source. We report on the modelling and computational simulation of physical and chemical processes taking place in a novel plasma-microfluidic platform. The main production and transport pathways of reactive species both in plasma and liquid are modelled by a novel modelling approach that combines 0D chemical kinetics and 2D transport mechanisms. This combined approach, applicable to systems where the transport of chemical species occurs in unidirectional flows at high Péclet numbers, decreases calculation times considerably compared to regular 2D simulations. It takes advantage of the low computational time of the 0D reaction models while providing spatial information through multiple plug-flow simulations to yield a quasi-2D model. The gas and liquid flow profiles are simulated entirely in 2D, together with the chemical reactions and transport of key chemical species. The model correctly predicts increased transport of hydrogen peroxide into the liquid when the microfluidic opening is placed inside the plasma effluent region, as opposed to inside the plasma region itself. Furthermore, the modelled hydrogen peroxide production and transport in the microfluidic liquid differs by less than 50% compared with experimental results. To explain this discrepancy, the limits of the 0D–2D combined approach are discussed.