Abstract: Reactive molecular dynamics simulations are performed to study the interaction of reactive oxygen species, such as OH, HO2 and H2O2, with the endotoxic biomolecule lipid A of the gram-negative bacterium Escherichia coli. It is found that the aforementioned plasma species can destroy the lipid A, which consequently results in reducing its toxic activity. All bond dissociation events are initiated by hydrogen-abstraction reactions. However, the mechanisms behind these dissociations are dependent on the impinging plasma species, i.e. a clear difference is observed in the mechanisms upon impact of HO2 radicals and H2O2 molecules on one hand and OH radicals on the other hand. Our simulation results are in good agreement with experimental observations.

Abstract: Aquaporins (AQPs) are transmembrane proteins that conduct not only water molecules across the cell membrane but also other solutes, such as reactive oxygen and nitrogen species (RONS), produced (among others) by cold atmospheric plasma (CAP). These RONS may induce oxidative stress in the cell interior, which plays a role in cancer treatment. The underlying mechanisms of the transport of RONS across AQPs, however, still remain obscure. We apply molecular dynamics simulations to investigate the permeation of both hydrophilic (H₂O₂and OH) and hydrophobic (NO₂and NO) RONS through AQP1. Our simulations show that these RONS can all penetrate across the pores of AQP1. The permeation free energy barrier of OH and NO is lower than that of H₂O₂and NO₂, indicating that these radicals may have easier access to the pore interior and interact with the amino acid residues of AQP1. We also study the effect of RONS-induced oxidation of both the phospholipids and AQP1 (i.e., sulfenylation of Cys₁₉₁) on the transport of the above-mentioned RONS across AQP1. Both lipid and protein oxidation seem to slightly increase the free energy barrier for H₂O₂and NO₂permeation, while for OH and NO, we do not observe a strong effect of oxidation. The simulation results help to gain insight in the underlying mechanisms of the noticeable rise of CAP-induced RONS in cancer cells, thereby improving our understanding on the role of AQPs in the selective anticancer capacity of CAP.

Abstract: Background: Strong electric fields are knownto affect cell membrane permeability,which can be applied for therapeutic purposes, e.g., in cancer therapy. A synergistic enhancement of this effect may be accomplished by the presence of reactive oxygen species (ROS), as generated in cold atmospheric plasmas. Little is known about the synergy between lipid oxidation by ROS and the electric field, nor on howthis affects the cell membrane permeability.
Method: We here conduct molecular dynamics simulations to elucidate the dynamics of the permeation process under the influence of combined lipid oxidation and electroporation. A phospholipid bilayer (PLB), consisting of di-oleoyl-phosphatidylcholine molecules covered with water layers, is used as a model system for the plasma membrane.
Results and conclusions:Weshow howoxidation of the lipids in the PLB leads to an increase of the permeability of the bilayer to ROS, although the permeation free energy barriers still remain relatively high. More importantly, oxidation of the lipids results in a drop of the electric field threshold needed for pore formation (i.e., electroporation) in the PLB. The created pores in the membrane facilitate the penetration of reactive plasma species deep into the cell interior, eventually causing oxidative damage.
General significance: This study is of particular interest for plasma medicine, as plasma generates both ROS and electric fields, but it is also of more general interest for applications where strong electric fields and ROS both come into play.

Abstract: Experiments have demonstrated the potential selective anticancer capacity of cold atmospheric plasmas (CAPs), but the underlying mechanisms remain unclear. Using computer simulations, we try to shed light on the mechanism of selectivity, based on aquaporins (AQPs), i.e. transmembrane protein channels transferring external H 2 O 2 and other reactive oxygen species, created e.g. by CAPs, to the cell interior. Specifically, we perform molecular dynamics simulations for the permeation of H 2 O 2 through AQP1 (one of the members of the AQP family) and the palmitoyl-oleoyl-phosphatidylcholine (POPC) phospholipid bilayer (PLB). The free energy barrier of H 2 O 2 across AQP1 is lower than for the POPC PLB, while the permeability coefficient, calculated using the free energy and diffusion rate profiles, is two orders of magnitude higher. This indicates that the delivery of H 2 O 2 into the cell interior should be through AQP. Our study gives a better insight into the role of AQPs in the selectivity of CAPs for treating cancer cells.

Abstract: Hoofd-halskanker vormt een aanzienlijke uitdaging met bijna 900.000 nieuwe diagnoses per jaar, waarbij de jaarlijkse incidentie blijft stijgen. Vaak wordt de diagnose pas in een laat stadium gesteld, wat complexe behandelingen noodzakelijk maakt. Terugval van patiënten is helaas een veelvoorkomend probleem. De gemiddelde overlevingsduur is beperkt tot enkele maanden. Daarom is er een dringende behoefte om nieuwe, veelbelovende behandelingen te ontwikkelen voor patiënten met hoofd-halskanker. Voor het bereiken van deze vooruitgang spelen innovatieve studiemodellen een cruciale rol. Het ontwikkelen van deze nieuwe behandelingen start met laboratoriumonderzoek, waarbij traditionele tweedimensionale celculturen hun beperkingen hebben. Daarom verschuiven onderzoekers hun aandacht meer en meer naar geavanceerdere driedimensionale modellen, met hoofd-halskankerorganoïden als beloftevol nieuw model. Dit model behoudt immers zowel het genetische profiel als de morfologische kenmerken van de originele tumor van de hoofd-halskankerpatiënt. Hoofdhalskankerorganoïden bieden daarom de mogelijkheid om innovatieve behandelingen te testen en kunnen mogelijk zelfs de respons van een patiënt op bepaalde therapieën voorspellen. Hoewel tumororganoïden als ‘patiënt-in-het-lab’ veelbelovend zijn, zijn er uitdagingen te overwinnen, zoals de ontwikkelingstijd en de toepasbaarheid bij alle tumortypes, evenals het ontbreken van immuuncellen en andere micro-omgevingscomponenten. Er is daarom een grote behoefte aan gestandaardiseerde protocollen voor de ontwikkeling van organoïden en verkorting van de ontwikkelingstijd. Concluderend bieden driedimensionale hoofd-halskankerorganoïden een veelbelovend perspectief voor de toekomst van kankerbehandelingen. Ze hebben het potentieel om bij te dragen aan de ontwikkeling van gepersonaliseerde behandelingen en zo de overlevingskansen van kankerpatiënten te verbeteren. Het is echter belangrijk om hun voorspellend vermogen en toepassingsmogelijkheden verder te onderzoeken, voordat ze op grote schaal worden geïmplementeerd.

Abstract: In this work, a chemical kinetics study on methane activation for hydrogen production in a warm plasma, i.e., N-2 rotating gliding arc (RGA), was performed for the first time to get new insights into the underlying reaction mechanisms and pathways. A zero-dimensional chemical kinetics model was developed, which showed a good agreement with the experimental results in terms of the conversion of CH4 and product selectivities, allowing us to get a better understanding of the relative significance of various important species and their related reactions to the formation and loss of CH4, H-2, and C2H2 etc. An overall reaction scheme was obtained to provide a realistic picture of the plasma chemistry. The results reveal that the electrons and excited nitrogen species (mainly N-2(A)) play a dominant role in the initial dissociation of CH4. However, the H atom induced reaction CH4+ H -> CH3+ H-2, which has an enhanced reaction rate due to the high gas temperature (over 1200 K), is the major contributor to both the conversion of CH4 and H-2 production, with its relative contributions of > 90% and > 85%, respectively, when only considering the forward reactions. The coexistence and interaction of thermochemical and plasma chemical processes in the rotating gliding arc warm plasma significantly enhance the process performance. The formation of C-2 hydrocarbons follows a nearly one-way path of C2H6 -> C2H4 -> C2H2, explaining why the selectivities of C-2 products decreased in the order of C2H2 > C2H4 > C2H6.

Abstract: Low-temperature plasma has received extensive attention due to its promising application prospects in the field of air pollutants degradation and energy conversion. To fulfill the need for particular applications, constructing stable plasma sources and investigating the interaction mechanisms between plasma and substances have been hot research topics. This thesis reports the diagnosis and improvement of plasma sources, diagnosis of the active species in plasma and a modeling study of chemical kinetics processes. The main research contents are as follows: In Chapter 3, a diffuse sine AC dielectric barrier discharge (DBD) is successfully obtained by optimizing the electrode structure. It is found that using double-layer dielectric plates can limit the discharge current intensity and significantly improve the discharge uniformity. The electrical characteristics and gas temperature with different operating time show that the discharge stability is also improved by using double-layer dielectric plates. In Chapter 4, nanosecond pulses are employed to generate diffuse DBD plasmas. Three main discharge stages are distinguished by ICCD images, i.e., the streamer breakdown from the needle tip to the plate electrode, the regime transition from streamer to diffuse plasma, and the propagation of surface discharge on the plate electrode surface. The chapter reveales that in nanosecond pulsed discharges the vibrational temperature of N2 increases with the discharge duration, while the rotational temperature mainly stays constant, which means electron energy is transferred into the vibrational levels, but gas heating is not obvious during the discharge pulse. In Chapter 5, both sine AC DBD and nanosecond pulsed DBD, studied in Chapter 2 and 3, are used for formaldehyde degradation. It is found that nanosecond pulsed DBD has more homogenous characteristics, better stability, and lower plasma gas temperature. Moreover, the energy consumption of nanosecond pulsed DBD is much lower than that of AC DBD. In Chapter 6, a 0D chemical kinetics model is developed to investigate the underlying plasma chemistry of methane dry reforming in a nanosecond pulsed discharge. An overview of the dominant reaction pathways of CO2 and CH4 conversion into the major products is given. Furthermore, most of the CO2 molecules are populated into vibrational states during the pulse. Hence, the vibrational states of CO2 play an important role in its dissociation process. In general, this PhD thesis contributes to a better insight in the mechanisms of sinusoidal AC DBD and nanosecond pulsed DBD plasmas and their applications, i.e., decomposition of formaldehyde and dry reforming of methane.

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 catalysis is gaining increasing interest for various environmental applications. Catalytic
material can be inserted in different shapes in the plasma, e.g., as pellets, (coated) beads, but also
as honeycomb monolith and 3DFD structures, also called ‘structured catalysts’, which have high
mass and heat transfer properties. In this work, we examine the streamer discharge propagation
and the interaction between plasma and catalysts, inside the channels of such structured catalysts,
by means of a two-dimensional particle-in-cell/Monte Carlo collision model. Our results reveal
that plasma streamers behave differently in various structured catalysts. In case of a honeycomb
structure, the streamers are limited to only one channel, with low or high plasma density when
the channels are parallel or perpendicular to the electrodes, respectively. In contrast, in case of a
3DFD structure, the streamers can distribute to different channels, causing discharge
enhancement due to surface charging on the dielectric walls of the structured catalyst, and
especially giving rise to a broader plasma distribution. The latter should be beneficial for plasma
catalysis applications, as it allows a larger catalyst surface area to be exposed to the plasma.

Abstract: Porous materials are commonly used in microelectronics, as they can meet the demand for continuously shrinking electronic feature dimensions. However, they are facing severe challenges in plasma etching, due to plasma induced damage. In this paper, we present both the plasma characteristics and surface processing during the etching of porous materials. We explain how the damage occurs in the porous material during plasma etching for a wide range of chuck temperatures and the responsible mechanism for plasma damage-free etching at cryogenic temperature, by a combination of experiments and numerical modeling.

Abstract: We studied hydrogen sulfide (H2S) decomposition into hydrogen (H2) and sulfur (S2) in a gliding arc plasmatron (GAP) and microwave (MW) plasma by a combination of 0D and 2D models. The conversion, energy efficiency, and plasma distribution are examined for different discharge conditions, and validated with available experiments from literature. Furthermore, a comparison is made between GAP and MW plasma. The GAP operates at atmospheric pressure, while the MW plasma experiments to which comparison is made were performed at reduced pressure. Indeed, the MW discharge region becomes very much contracted near atmospheric pressure, at the conditions under study, as revealed by our 2D model. The models predict that thermal reactions play the most important role in H2S decomposition in both plasma types. The GAP has a higher energy efficiency but lower conversion than the MW plasma at their typical conditions. When compared at the same conversion, the GAP exhibits a higher energy efficiency and lower energy cost than the MW plasma.

Abstract: Both geometrical and electrical asymmetry effects in capacitive argon discharges are investigated using a two-dimensional particle-in-cell coupled with Monte Carlo collision model. When changing the ratio of the top and bottom electrode surface areas and the phase shift between the two applied harmonics, the induced self-bias was found to develop separately. By adjusting the ratio between the high and low harmonic amplitudes, the electrical asymmetry effect at a fixed phase shift can be substantially optimized. However, the self-bias caused by the geometrical asymmetry hardly changed. Moreover, the separate control of these two asymmetry effects can also be demonstrated from their power absorption profiles. Both the axial and radial plasma density distributions can be modulated by the electrical asymmetry effect.

Abstract: In this work, a time modulated RF driven DBD-like atmospheric pressure plasma jet in Ar + 2%O2, operating at a time averaged power of 6.5 W is investigated. Spatially resolved ozone densities and gas temperatures are obtained by UV absorption and Rayleigh scattering, respectively. Significant gas heating in the core of the plasma up to 700 K is found and at the position of this increased gas temperature a depletion of the ozone density is found. The production and destruction reactions of O3 in the jet effluent as a function of the distance from the nozzle are obtained from a zero-dimensional chemical kinetics model in plug flow mode which considers relevant air chemistry due to air entrainment in the jet fluent. A comparison of the measurements and the models show that the depletion of O3 in the core of the plasma is mainly caused by an enhanced destruction of O3 due to a large atomic oxygen density.

Abstract: Computer simulations are performed for an argon inductively coupled plasma (ICP) with a capacitive radio-frequency bias power, to investigate the bias effect on the discharge mode transition and on the plasma characteristics at various ICP currents, bias voltages, and bias frequencies. When the bias frequency is fixed at 13.56 MHz and the ICP current is low, e.g., 6A, the spatiotemporal averaged plasma density increases monotonically with bias voltage, and the bias effect is already prominent at a bias voltage of 90 V. The maximum of the ionization rate moves toward the bottom electrode, which indicates clearly the discharge mode transition in inductive/capacitive discharges. At higher ICP currents, i.e., 11 and 13 A, the plasma density decreases first and then increases with bias voltage, due to the competing mechanisms between the ion acceleration power dissipation and the capacitive power deposition. At 11 A, the bias effect is still important, but it is noticeable only at higher bias voltages. At 13 A, the ionization rate is characterized by a maximum at the reactor center near the dielectric window at all selected bias voltages, which indicates that the ICP power, instead of the bias power, plays a dominant role under this condition, and no mode transition is observed. Indeed, the ratio of the bias power to the total power is lower than 0.4 over a wide range of bias voltages, i.e., 0300V. Besides the effect of ICP current, also the effect of various bias frequencies is investigated. It is found that the modulation of the bias power to the spatiotemporal distributions of the ionization rate at 2MHz is strikingly different from the behavior observed at higher bias frequencies. Furthermore, the minimum of the plasma density appears at different bias voltages, i.e., 120V at 2MHz and 90V at 27.12 MHz.

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.