Abstract: A gliding arc plasmatron (GAP) is very promising for CO2 conversion into value-added chemicals, but to further improve this important application, a better understanding of the arc behavior is indispensable. Therefore, we study here for the first time the dynamic arc behavior of the GAP by means of a high-speed camera, for different reactor configurations and in a wide range of operating conditions. This allows us to provide a complete image of the behavior of the gliding arc. More specifically, the arc body shape, diameter, movement and rotation speed are analyzed and discussed. Clearly, the arc movement and shape relies on a number of factors, such as gas turbulence, outlet diameter, electrode surface, gas contraction and buoyance force. Furthermore, we also compare the experimentally measured arc movement to a state-of-the-art 3D-plasma model, which predicts the plasma movement and rotation speed with very good accuracy, to gain further insight in the underlying mechanisms. Finally, we correlate the arc dynamics with the CO2 conversion and energy efficiency, at exactly the same conditions, to explain the effect of these parameters on the CO2 conversion process. This work is important for understanding and optimizing the GAP for CO2 conversion.

Abstract: This work explains the need for plasma models, introduces arguments for choosing the type of model that better fits the purpose of each study, and presents the basics of the most common nonequilibrium low-temperature plasma models and the information available from each one, along with an extensive list of references for complementary in-depth reading. The paper presents the following models, organised according to the level of multi-dimensional description of the plasma: kinetic models, based on either a statistical particle-in-cell/Monte-Carlo approach or the solution to the Boltzmann equation (in the latter case, special focus is given to the description of the electron kinetics); multi-fluid models, based on the solution to the hydrodynamic equations; global (spatially-average) models, based on the solution to the particle and energy rate-balance equations for the main plasma species, usually including a very complete reaction chemistry; mesoscopic models for plasma–surface interaction, adopting either a deterministic approach or a stochastic dynamical Monte-Carlo approach. For each plasma model, the paper puts forward the physics context, introduces the fundamental equations, presents advantages and limitations, also from a numerical perspective, and illustrates its application with some examples. Whenever pertinent, the interconnection between models is also discussed, in view of multi-scale hybrid approaches.

Abstract: Although plasma catalysis is gaining increasing interest for various environmental applications, the underlying mechanisms are still far from understood. For instance, it is not yet clear whether and how plasma streamers can propagate in catalyst pores, and what is the minimum pore size to make this happen. As this is crucial information to ensure good plasma-catalyst interaction, we study here the mechanism of plasma streamer propagation in a catalyst pore, by means of a twodimensional particle-in-cell/Monte Carlo collision model, for various pore diameters in the nm range to μm-range. The so-called Debye length is an important criterion for plasma penetration into catalyst pores, i.e. a plasma streamer can penetrate into pores when their diameter is larger than the Debye length. The Debye length is typically in the order of a few 100 nm up to 1 μm at the conditions under study, depending on electron density and temperature in the plasma streamer. For pores in the range of ∼50 nm, plasma can thus only penetrate to some extent and at
very short times, i.e. at the beginning of a micro-discharge, before the actual plasma streamer reaches the catalyst surface and a sheath is formed in front of the surface. We can make plasma streamers penetrate into smaller pores (down to ca. 500 nm at the conditions under study) by increasing the applied voltage, which yields a higher plasma density, and thus reduces the Debye length. Our simulations also reveal that the plasma streamers induce surface charging of the catalyst pore sidewalls, causing discharge enhancement inside the pore, depending on pore diameter and depth.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

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: A gliding arc (GA) plasma has great potential for CO2 conversion into value-added chemicals, because of its high energy efficiency. To improve the application, a 2D/3D fluid model is needed to investigate the CO2 conversion mechanisms in the actual discharge geometry. Therefore, the complex CO2 chemical kinetics description must be reduced due to the huge computational cost associated with 2D/3D models. This paper presents a chemistry reduction method for CO2 plasmas, based on the so-called directed relation graph method. Depending on the defined threshold values, some marginal species are identified. By means of a sensitivity analysis, we can further reduce the chemistry set by removing one by one the marginal species. Based on the socalled flux-sensitivity coupling, we obtain a reduced CO2 kinetics model, consisting of 36 or 15 species (depending on whether the 21 asymmetric mode vibrational states of CO2 are explicitly included or lumped into one group), which is applied to a GA discharge. The results are compared with those predicted with the full chemistry set, and very good agreement is reached. Moreover, the range of validity of the reduced CO2 chemistry set is checked, telling us that this reduced set is suitable for low power GA discharges. Finally, the time and spatial evolution of the CO2 plasma characteristics are presented, based on a 2D model with the reduced kinetics.

Abstract: This study focuses on a magnetically stabilized gliding arc (MGA) plasma. Two fully coupled flow-plasma models (in 3D and 2D) are presented. The 3D model is applied to compare the arc dynamics of the MGA with a traditional gas-driven gliding arc. The 2D model is used for a detailed parametric study on the effect of the external magnetic field. The results show that the relative velocity between the plasma and feed gas is generated due to the Lorentz force, which can increase the plasma-treated gas fraction. The magnetic field also helps to decrease the gas temperature by enhancing heat transfer and to increase the electron number density. This work shows the potential of an external magnetic field to control the gliding arc behavior, for enhanced gas conversion at low gas flow rates.

Abstract: We studied the evolution of an isolated pulsed plasma in a vortex flow stabilised microwave (MW) discharge in N2 at 25 mbar via the combination of 0D kinetics modelling, iCCD imaging and laser scattering diagnostics. Quenching of electronically excited N2 results in fast gas heating and the onset of a thermal-ionisation instability, contracting the discharge volume. The onset of a thermal-ionisation instability driven by vibrational excitation pathways is found to facilitate significantly higher N2 conversion (i.e. dissociation to atomic N2 ) compared to pre-instability conditions, emphasizing the potential utility of this dynamic in future fixation applications. The instability onset is found to be instigated by super-elastic heating of the electron energy distribution tail via vibrationally excited N2 . Radial contraction of the discharge to the skin depth is found to occur post instability, while the axial elongation is found to be temporarily contracted during the thermal instability onset. An increase in power reflection during the thermal instability onset eventually limits the destabilising effects of exothermic electronically excited N2 quenching. Translational and vibrational temperature reach a quasi-non-equilibrium after the discharge contraction, with translational temperatures reaching ∼1200 K at the pulse end, while vibrational temperatures are found in near equilibrium with the electron energy (1 eV, or ∼11 600 K). This first description of the importance of electronically excited N2 quenching in thermal instabilities gives an additional fundamental understanding of N2 plasma behaviour in pulsed MW context, and thereby brings the eventual implementation of this novel N2 fixation method one step closer.

Abstract: Nanosecond pulsed discharge plasma shows a high degree of non-equilibrium, and exhibits relatively high conversions in the dry reforming of methane. To further improve the application, a good insight of the underlying mechanisms is desired. We developed a chemical kinetics model to explore the underlying plasma chemistry in nanosecond pulsed discharge. We compared the calculated conversions and product selectivities with experimental results, and found reasonable agreement in a wide range of specific energy input. Hence, the chemical kinetics model is able to provide insight in the underlying plasma chemistry. The modeling results predict that the most important dissociation reaction of CO₂and CH₄is electron impact dissociation. C₂H₂is the most abundant hydrocarbon product, and it is mainly formed upon reaction of two CH₂radicals. Furthermore, the vibrational excitation levels of CO₂contribute for 85% to the total dissociation of CO₂.

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: 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.

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: Background
During normal pregnancy, placental oxidative stress (OS) is present during all three trimesters and is necessary to obtain normal cell function. However, if OS reaches a certain level, pregnancy complications might arise. In preeclampsia (PE), a dangerous pregnancy specific hypertensive disorder, OS induced in the ischemic placenta causes a systemic inflammatory response and activates maternal endothelial cells. In this study, we aimed to quantify superoxide concentrations (as a measure of systemic OS) using electron paramagnetic resonance (EPR) and correlate them to markers of systemic inflammation, iron status and vascular function.
Methods
Fifty-nine women with a healthy pregnancy (HP), 10 non-pregnant controls (NP) and 28 PE patients (32±3.3weeks) were included. During HP, blood samples for superoxide, neutrophil to lymphocyte ratio (NLR), mean platelet volume (MPV) and iron status were taken at 10, 25 and 39 weeks. Vascular measurements for arterial stiffness (carotid-femoral pulse wave velocity (CF-PWV), augmentation index (AIx), augmentation Pressure (AP)) and microvascular endothelial function (reactive hyperemia index (RHI)) were performed at 35 weeks. In PE, all measurements were performed at diagnosis. CMH (1-hydroxy-3-methoxycarbonyl-2,2,5,5-tetramethylpyrrolidine) was used as spin probe for EPR, since the formed CM radical
corresponds to the amount of superoxide.
Results
Superoxide concentration remains stable during pregnancy (p = 0.92), but is significantly higher compared to the NP controls (p<0.0001). At 25 weeks, there is a significant positive correlation between superoxide and ferritin concentration. (p = 0.04) In PE, superoxide, systemic inflammation and iron status are much higher compared to HP (all p<0.001). During HP, superoxide concentrations correlate significantly with arterial stiffness (all p<0.04), while in PE superoxide is significantly correlated to microvascular endothelial function (p = 0.03).
Conclusions
During HP there is an increased but stable oxidative environment, which is correlated to ferritin concentration. If superoxide levels increase, there is an augmentation in arterial stiffness. In PE pregnancies, systemic inflammation and superoxide concentrations are higher and result in a deterioration of endothelial function. Together, these findings support the hypothesis that vascular function is directly linked to the amount of OS and that measurement of OS in combination with vascular function tests might be used in the prediction of PE.

Abstract: Many current anti-cancer therapies rely on increasing the intracellular reactive oxygen and nitrogen species (RONS) contents with the aim to induce irreparable damage, which subsequently results in tumor cell death. A novel tool in cancer therapy is the use of cold atmospheric plasma (CAP), which has been found to be very effective in the treatment of many different cancer cell types in vitro as well as in vivo, mainly through the vast generation of RONS. One of the key determinants of the cell's fate will be the interaction of RONS, generated by CAP, with important proteins, i.e. redox-regulatory proteins. One such protein is cytoglobin (CYGB), a recently discovered globin proposed to be involved in the protection of the cell against oxidative stress. In this study, the effect of plasma-produced RONS on CYGB was investigated through the treatment of CYGB with CAP for different treatment times. Spectroscopic analysis of CYGB showed that although chemical modifications occur, its secondary structure remains intact. Mass spectrometry experiments identified these modifications as oxidations of mainly sulfur-containing and aromatic amino acids. With longer treatment time, the treatment was also found to induce nitration of the heme. Furthermore, the two surface-exposed cysteine residues of CYGB were oxidized upon treatment, leading to the formation of intermolecular disulfide bridges, and potentially also intramolecular disulfide bridges. In addition, molecular dynamics and docking simulations confirmed, and further show, that the formation of an intramolecular disulfide bond, due to oxidative conditions, affects the CYGB 3D structure, thereby opening the access to the heme group, through gate functioning of His117. Altogether, the results obtained in this study (1) show that plasma-produced RONS can extensively oxidize proteins and (2) that the oxidation status of two redox-active cysteines lead to different conformations of CYGB.

Abstract: Melanoma arises from pigment-producing cells called melanocytes located in the basal layers of the epidermis of the skin. Cytoglobin (CYGB) is a ubiquitously expressed hexacoordinated globin that is highly enriched in me­lanocytes and frequently downregulated during melanomagenesis. Previously, we showed that non-thermal plasma (NTP)-produced reactive oxygen and nitrogen species (RONS) lead to the formation of an intra­ molecular disulfide bridge that would allow CYGB to function as a redox-sensitive protein. Here, we investigate the cytotoxic effect of indirect NTP treatment in two melanoma cell lines with divergent endogenous CYGB expression levels, and we explore the role of CYGB in determining treatment outcome. Our findings are consistent with previous studies supporting that NTP cytotoxicity is mediated through the production of RONS and leads to apoptotic cell death in melanoma cells. Furthermore, we show that NTP-treated solutions elicit an antioxidant response through the activation of nuclear factor erythroid 2–related factor 2 (NRF2). The knock­ down and overexpression of CYGB respectively sensitizes and protects melanoma cells from RONS-induced apoptotic cell death. The presence of CYGB enhances heme-oxygenase 1 (HO-1) and NRF2 protein expression levels, whereas the absence impairs their expression. Moreover, analysis of the CYGB-dependent transcriptome demonstrates the tumor suppressor long non-coding RNA maternally expressed 3 (MEG3) as a hitherto unde­ scribed link between CYGB and NRF2. Thus, the presence of CYGB, at least in melanoma cells, seems to play a central role in determining the therapeutic outcome of RONS-inducing anticancer therapies, like NTP-treated solutions, possessing both tumor-suppressive and oncogenic features. Hence, CYGB expression could be of in­ terest either as a biomarker or as a candidate for future targeted therapies in melanoma.