Abstract: Modelling at the molecular or atomic scale can be very useful for obtaining a better insight in plasma medicine. This paper gives an overview of different atomic/molecular scale modelling approaches that can be used to study the direct interaction of plasma species with biomolecules or the consequences of these interactions for the biomolecules on a somewhat longer time-scale. These approaches include density functional theory (DFT), density functional based tight binding (DFTB), classical reactive and non-reactive molecular dynamics (MD) and united-atom or coarse-grained MD, as well as hybrid quantum mechanics/molecular mechanics (QM/MM) methods. Specific examples will be given for three important types of biomolecules, present in human cells, i.e. proteins, DNA and phospholipids found in the cell membrane. The results show that each of these modelling approaches has its specific strengths and limitations, and is particularly useful for certain applications. A multi-level approach is therefore most suitable for obtaining a global picture of the plasma–biomolecule interactions.

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: 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: 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: Given the ecological problems associated to the CO2 emissions of fossil fuels, plasma technology has gained
interest for conversion of CO2 into value-added products. Microwave plasmas operating at atmospheric pressure
have proven to be especially interesting, due to the high gas temperatures inside the reactor (i.e. up to 6000 K)
allowing for efficient thermal dissociation of CO2 into CO and O2. However, the performance of these high
temperature plasmas is limited by recombination of CO back into CO2 once the gas cools down in the afterglow.
In this work, we computationally investigated several quenching nozzles, developed and experimentally tested
by Hecimovic et al., [1] for their ability to quickly cool the gas after the plasma, thereby quenching the CO
recombination reactions. Using a 3D computational fluid dynamics model and a quasi-1D chemical kinetics
model, we reveal that a reactor without nozzle lacks gas mixing between hot gas in the center and cold gas near
the reactor walls. Especially at low flow rates, where there is an inherent lack of convective cooling due to the
low gas flow velocity, the temperature in the afterglow remains high (between 2000 and 3000 K) for a relatively
long time (in the 0.1 s range). As shown by our quasi-1D chemical kinetics model, this results in a important loss
of CO due to recombination reactions. Attaching a nozzle in the effluent of the reactor induces fast gas quenching
right after the plasma. Indeed, it introduces (i) more convective cooling by forcing cool gas near the walls to mix
with hot gas in the center of the reactor, as well as (ii) more conductive cooling through the water-cooled walls of
the nozzle. Our model shows that gas quenching and the suppression of recombination reactions have more
impact at low flow rates, where recombination is the most limiting factor in the conversion process.

Abstract: Plasma catalysis is emerging for plasma-assisted gas conversion
processes. However, the underlying mechanisms of plasma catalysis are poorly
understood. In this work, we present a 1D heterogeneous catalysis model with axial
dispersion (i.e., accounting for back-mixing and molecular diffusion of fluid elements in
the process stream in the axial direction), for plasma-catalytic NO production from
N2/O2 mixtures. We investigate the concentration and reaction rates of each species
formed as a function of time and position across the catalyst, in order to determine the
underlying mechanisms. To obtain insights into how the performance of the process
can be further improved, we also study how changes in the postplasma gas flow
composition entering the catalyst bed and in the operation conditions of the catalytic
stage affect the performance of NO production.

Abstract: Plasma-based CO2 and CH4 conversion is gaining increasing interest, and a great portion of research is dedicated to adapting the process to actual industrial conditions. In an industrial context, the process needs to be able to process N2 admixtures, since most industrial gas emissions contain significant amounts of N2, and gas separations are financially costly. In this paper we therefore investigate the effect of N2 on the CO2 and CH4 conversion in a gliding arc plasmatron reactor. The addition of 20 % N2 reduces the energy cost of the conversion process by 21 % compared to a pure CO2/CH4 mixture, from 2.9 down to 2.2 eV/molec (or from 11.5 to 8.7 kJ/L), yielding a CO2 and CH4 (absolute) conversion of 28.6 and 35.9 % and an energy efficiency of 58 %. These results are among the best reported in literature for plasma-based DRM, demonstrating the benefits of N2 present in the mix. Compared to DRM results in different plasma reactor types, a low energy cost was achieved. To understand the underlying mechanisms of N2 addition, we developed a combination of four different computational models, which reveal that the beneficial effect of N2 addition is attributed to (i) a rise in the electron density (increasing the plasma conductivity, and therefore reducing the plasma power needed to sustain the plasma, which reduces the energy cost), as well as (ii) a rise in the gas temperature, which accelerates the CO2 and CH4 conversion reactions.

Abstract: Transforming CO2 into value-added chemicals is crucial to realizing a carbon–neutral economy, and plasma-based conversion, a Power-2-X technology, offers a promising route to realizing an efficient and scalable process. This paper investigates the effects of post-plasma placement of a converging–diverging nozzle in a vortex-stabilized 2.45 GHz CO2 microwave plasma reactor to increase energy efficiency and conversion. The CDN leads to a 21 % relative increase in energy efficiency (31 %) and CO2 conversion (13 %) at high flow rates and near-atmospheric conditions. The most significant performance improvement was seen at low flow rates and sub-atmospheric pressure (300 mbar), where energy efficiency was 23 % and conversion was 28 %, a 71 % relative increase over conditions without the CDN. Using CFD simulations, we found that the CDN produces a change in the flow geometry, leading to a confined temperature profile at the height of the plasma, and forced extraction of CO to the post-CDN region.

Abstract: CO 2 conversion by a gliding arc plasma is gaining increasing interest, but the underlying mechanisms for an energy-efficient process are still far from understood. Indeed, the chemical complexity of the non-equilibrium plasma poses a challenge for plasma modeling due to the huge computational load. In this paper, a one-dimensional (1D) gliding arc model is developed in a cylindrical frame, with a detailed non-equilibrium CO 2 plasma chemistry set, including the CO 2 vibrational kinetics up to the dissociation limit. The model solves a set of time- dependent continuity equations based on the chemical reactions, as well as the electron energy balance equation, and it assumes quasi-neutrality in the plasma. The loss of plasma species and heat due to convection by the transverse gas flow is accounted for by using a characteristic frequency of convective cooling, which depends on the gliding arc radius, the relative velocity of the gas flow with respect to the arc and on the arc elongation rate. The calculated values for plasma density and plasma temperature within this work are comparable with experimental data on gliding arc plasma reactors in the literature. Our calculation results indicate that excitation to the vibrational levels promotes efficient dissociation in the gliding arc, and this is consistent with experimental investigations of the gliding arc based CO 2 conversion in the literature. Additionally, the dissociation of CO 2 through collisions with O atoms has the largest contribution to CO 2 splitting under the conditions studied. In addition to the above results, we also demonstrate that lumping the CO 2 vibrational states can bring a significant reduction of the computational load. The latter opens up the way for 2D or 3D models with an accurate description of the CO 2 vibrational kinetics.

Abstract: A gliding arc plasma is a potential way to convert CO2 into CO and O2, due to its non-equilibrium character, but little is known about the underlying mechanisms. In this paper, a self-consistent two-dimensional (2D) gliding arc model is developed, with a detailed non-equilibrium CO2 plasma chemistry, and validated with experiments. Our calculated values of the electron number density in the plasma, the CO2 conversion and energy efficiency show reasonable agreement with the experiments, indicating that the model can provide a realistic picture of the plasma chemistry. Comparison of the results with classical thermal conversion, as well as other plasma-based technologies for CO2 conversion reported in literature, demonstrates the non-equilibrium character of the gliding arc, and indicates that the gliding arc is a promising plasma reactor for CO2 conversion. However, some process modifications should be exploited to further improve its performance. As the model provides a realistic picture of the plasma behaviour, we use it first to investigate the plasma characteristics in a whole gliding arc cycle, which is necessary to understand the underlying mechanisms. Subsequently, we perform a chemical kinetics analysis, to investigate the different pathways for CO2 loss and formation. Based on the revealed discharge properties and the underlying CO2 plasma chemistry, the model allows us to propose solutions on how to further improve the
CO2 conversion and energy efficiency by a gliding arc plasma.

Abstract: Cancer cells are characterized by higher levels of reactive oxygen species (ROS) compared to normal cells as a result of an imbalance between oxidants and antioxidants. However, cancer cells maintain their redox balance due to their high antioxidant capacity. Recently, a high level of oxidative stress is considered a novel target for anticancer therapy. This can be induced by increasing exogenous ROS and/or inhibiting the endogenous protective antioxidant system. Additionally, the immune system has been shown to be a significant ally in the fight against cancer. Since ROS levels are important to modulate the antitumor immune response, it is essential to consider the effects of oxidative stress-inducing treatments on this response. In this review, we provide an overview of the mechanistic cellular responses of cancer cells towards exogenous and endogenous ROS-inducing treatments, as well as the indirect and direct antitumoral immune effects, which can be both immunostimulatory and/or immunosuppressive. For future perspectives, there is a clear need for comprehensive investigations of different oxidative stress-inducing treatment strategies and their specific immunomodulating effects, since the effects cannot be generalized over different treatment modalities. It is essential to elucidate all these underlying immune effects to make oxidative stress-inducing treatments effective anticancer therapy.

Abstract: Plasma catalysis is gaining increasing interest for CO2 conversion, but the interaction between the plasma and catalyst is still poorly understood. This is caused by limited systematic materials research, since most works combine a plasma with commercial supported catalysts and packings. In the present paper, we study the influence of specific material and reactor properties, as well as reactor/bead configuration, on the conversion and energy efficiency of CO2 dissociation in a packed bed dielectric barrier discharge (DBD) reactor. Of the various packing materials investigated, BaTiO3 yields the highest conversion and energy efficiency, i.e., 25% and 4.5%.
Our results show that, when evaluating the influence of catalysts, the impact of the packing (support) material itself cannot be neglected, since it can largely affect the conversion and energy efficiency. This shows the large potential for further improvement of packed bed plasma reactors for CO2 conversion and other chemical conversion reactions by adjusting both packing (support) properties and catalytically active sites. Moreover, we clearly prove that comparison of results obtained in different reactor setups should be done with care, since there is a large effect of the reactor setup and reactor/bead configuration.

Abstract: The use of plasma technology for CO2 splitting is gaining increasing interest, but one of the major obstacles to date for industrial implementation is the considerable energy cost. We demonstrate that the introduction of a packing of dielectric zirconia (ZrO2) beads into a dielectric barrier discharge (DBD) plasma reactor can enhance the CO2 conversion and energy efficiency up to a factor 1.9 and 2.2, respectively, compared to that in a normal (unpacked) DBD reactor. We obtained a maximum conversion of 42 % and a maximum energy efficiency of 9.6 %. However, it is the ability of the packing to almost double both the conversion and the energy efficiency simultaneously at certain input parameters that makes it very promising. The improved conversion and energy efficiency can be explained by the higher values of the local electric field and electron energy near the contact points of the beads and the lower breakdown voltage, demonstrated by 2 D fluid modeling.

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.