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: A packed bed dielectric barrier discharge plasma reactor is computationally studied with a fluid model. Two different complementary axisymmetric 2D geometries are used to mimic the intrinsic 3D problem. It is found that a packing enhances the electric field strength and electron temperature at the contact points of the dielectric material due to polarization of the beads by the applied potential. As a result, these contact points prove to be of direct importance to initiate the plasma. At low applied potential, the discharge stays at the contact points, and shows the properties of a Townsend discharge. When a high enough potential is applied, the plasma will be able to travel through the gaps in between the beads from wall to wall, forming a kind of glow discharge. Therefore, the inclusion of a so-called ‘channel of voids’ is indispensable in any type of packed bed modelling.

Abstract: Packed bed plasma reactors (PBPRs) are gaining increasing interest for use in environmental applications, such as greenhouse gas conversion into value-added chemicals or renewable fuels and volatile pollutant removal (e.g. NOx, VOC, K), as they enhance the conversion and energy efficiency of the process compared to a non-packed reactor. However, the plasma behaviour in a PBPR is not well understood. In this paper we demonstrate, by means of a fluid model, that the discharge behaviour changes considerably when changing the size of the packing beads and their dielectric constant, while keeping the interelectrode spacing constant. At low dielectric constant, the plasma is spread out over the full discharge gap, showing significant density in the voids as well as in the connecting void channels. The electric current profile shows a strong peak during each half cycle. When the dielectric constant increases, the plasma becomes localised in the voids, with a current profile consisting of many smaller peaks during each half cycle. For large bead sizes, the shift from full gap discharge to localised discharges takes place at a higher dielectric constant than for smaller beads. Furthermore, smaller beads or beads with a lower dielectric constant require a higher breakdown voltage to cause plasma formation.

Abstract: In recent years, cold atmospheric plasma (CAP) has been proposed as a novel method to enhance transdermal drug delivery, while avoiding tissue damage. However, the underlying mechanism for the increasing skin permeability upon CAP treatment is still undefined. We propose a mechanism in which CAP-generated reactive species induce cross-linking of skin lipids, leading to the generation of nanopores, thereby facilitating the permeation of drug molecules. Molecular dynamics simulations support this proposed mechanism. Furthermore, our results indicate that to achieve maximum enhancement of the permeability, the optimal treatment will depend on the exact lipid composition of the skin, as well as on the CAP source used.

Abstract: Promising biomedical uses for nonthermal plasma (NTP) in the fields of regenerative medicine, cancer therapy, and vaccine delivery involve the noninvasive application of uniform nonequilibrium plasma (including dielectric barrier discharge plasma) to living skin. Whereas most investigations have focused on achieving desired therapeutic outcomes, fewer studies have examined the mechanisms and pathways by which epithelial cells respond to NTP exposure. Using a transwell apical-basolateral-chambered system to culture the human keratinocyte HaCaT cell line, in vitro experiments were performed to demonstrate the effects of nanosecond-pulsed dielectric barrier discharge (nsDBD) plasma on polarized epithelial cell viability, monolayer permeability, intracellular oxidative stress, and the release of adenosine triphosphate (ATP). Application of nsDBD plasma at 60 Hz or below had minimal or no effect on HaCaT monolayer viability or permeability. nsDBD plasma exposure did, however, result in frequency-dependent reductions in intracellular glutathione (indicating direct induction of oxidative stress by nsDBD plasma) and increased extracellular ATP concentrations in the ba-solateral (subepithelial) media, which are indicators of cellular stress and an NTP-induced inflammatory response. These studies provide new insights into nsDBD plasma-induced inflammation and local innate immune responses initiated by polarized epithelial tissues.

Abstract: In this computational study, a gliding arc plasma reactor with a reverse-vortex flow stabilization is modelled for the first time by a fluid plasma description. The plasma reactor operates with argon gas at atmospheric pressure. The gas flow is simulated using the k-ε Reynolds-averaged Navier–Stokes turbulent model. A quasi-neutral fluid plasma model is used for computing the plasma properties. The plasma arc movement in the reactor is observed, and the results for the gas flow, electrical characteristics, plasma density, electron temperature, and gas temperature are analyzed.

Abstract: We derive dispersion relations for a system of identical particles confined in a two-dimensional annular harmonic well and which interact through a Yukawa potential, e.g., a dusty plasma ring. When the particles are in a single chain (i.e., a one-dimensional ring), we find a longitudinal acoustic mode and a transverse optical mode which show approximate agreement with the dispersion relation for a straight configuration for large radii of the ring. When the radius decreases, the dispersion relations modify: there appears an anticrossing of the modes near the crossing point resulting in a frequency gap between the lower and upper branches of the modified dispersion relations. For the double chain (i.e., a two-dimensional zigzag configuration), the dispersion relation has four branches: longitudinal acoustic and optical and transverse acoustic and optical.

Abstract: Single-file diffusion (SFD) of an infinite one-dimensional chain of interacting particles has a long-time mean-square displacement ∝t1/2, independent of the type of interparticle repulsive interaction. This behavior is also observed in finite-size chains, although only for certain intervals of time t depending on the chain length L, followed by the ∝t for t→∞, as we demonstrate for a closed circular chain of diffusing interacting particles. Here, we show that spatial correlation of noise slows down SFD and can result, depending on the amount of correlated noise, in either subdiffusive behavior ∝tα, where 0<α<1/2, or even in a total suppression of diffusion (in the limit N→∞). Spatial correlation can explain the subdiffusive behavior in recent SFD experiments in circular channels.

Abstract: Simulations and experiments are performed to obtain a better insight in the plasma enhanced chemical vapor deposition process of SiO2 by SiCl4/O2/Ar plasmas for introducing a SiO2-like coating in wafer processing reactors. Reaction sets describing the plasma and surface chemistry of the SiCl4/O2/Ar mixture are presented. Typical calculation results include the bulk plasma characteristics, i.e., electrical properties, species densities, and information on important production and loss processes, as well as the chemical composition of the deposited coating, and the thickness uniformity of the film on all reactor surfaces. The film deposition characteristics, and the trends for varying discharge conditions, are explained based on the plasma behavior, as calculated by the model.

Abstract: In this article, surface processes occurring during the etching of Si with a Cl2/O2/Ar plasma are investigated by means of experiments and modeling. Cl2-based plasmas are commonly used to etch silicon, while a small fraction of O2 is added to protect the sidewalls from lateral etching during the shallow trench isolation process. When the oxygen fraction exceeds a critical value, the wafer surface process changes from an etching regime to a deposition regime, drastically reducing the etch rate. This effect is commonly referred to as the etch stop phenomenon. To gain better understanding of this mechanism, the oxygen fraction is varied in the gas mixture and special attention is paid to the effects of oxygen and of the redeposition of non-volatile etched species on the overall etch/deposition process. It is found that, when the O2 flow is increased, the etch process changes from successful etching to the formation of a rough surface, and eventually to the actual growth of an oxide layer which completely blocks the etching of the underlying Si. The size of this etch stop island was found to increase as a function of oxygen flow, while its thickness was dependent on the amount of Si etched. This suggests that the growth of the oxide layer mainly depends on the redeposition of non-volatile etch products. The abrupt change in the etch rate as a function of oxygen fraction was not found back in the oxygen content of the plasma, suggesting the competitive nature between oxidation and chlorination at the wafer. Finally, the wafer and reactor wall compositions were investigated by modeling and it was found that the surface rapidly consisted mainly of SiO2 when the O2 flow was increased above about 15 sccm.

Abstract: Modeling results are presented to gain a better insight in the properties of a SiH4/O2/Ar inductively coupled plasma (ICP) and how it interacts with a silicon substrate (wafer), as applied in the microelectronics industry for the fabrication of electronic devices. The SiH4/O2/Ar ICP is used for the filling of microtrenches with isolating material (SiO2), as applied in shallow trench isolation (STI). In this article, a detailed reaction set that describes the plasma chemistry of SiH4/O2/Ar discharges as well as surface processes, such as sputtering, oxidation, and deposition, is presented. Results are presented on the plasma properties during the plasma enhanced chemical vapor deposition process (PECVD) for different gas ratios, as well as on the shape of the filled trenches and the surface compositions of the deposited layers. For the operating conditions under study it is found that the most important species accounting for deposition are SiH2, SiH3O, SiH3 and SiH2O, while SiH+2, SiH+3, O+2 and Ar+ are the dominant species for sputtering of the surface. By diluting the precursor gas (SiH4) in the mixture, the deposition rate versus sputtering rate can be controlled for a desired trench filling process. From the calculation results it is clear that a high deposition rate will result in undesired void formation during the trench filling, while a small deposition rate will result in undesired trench bottom and mask damage by sputtering. By varying the SiH4/O2 ratio, the chemical composition of the deposited layer will be influenced. However, even at the highest SiH4/O2 ratio investigated (i.e., 3.2:1; low oxygen content), the bulk deposited layer consists mainly of SiO2, suggesting that low-volatile silane species deposit first and subsequently become oxidized instead of being oxidized first in the plasma before deposition. Finally, it was found that the top surface of the deposited layer contained less oxygen due to preferential sputtering of O atoms, making the top layer more Si-rich. However, this effect is negligible at a SiH4/O2 ratio of 2:1 or lower.

Abstract: A numerical model for NxOy remediation in humid air plasma produced with a dielectric barrier discharge at atmospheric pressure is presented. Special emphasis is given to NO2 and N2O reduction with the decrease of O2 content in the feedstock gas. A detailed reaction mechanism including electronic and ionic processes, as well as the contribution of radicals and excited atomic/molecular species is proposed. The temporal evolution of the densities of NO, NO2 and N2O species, and some other by-products, is analyzed, and the major pathways for the NxOy remediation are discussed for one pulse. Subsequently, simulations are presented for a multi-pulses case, where three O2 contents are tested for optimization of the remediation process. It is found that when the gas mixture O2/N2/H2O/NOx has no initial O2 content, the best NOx and N2O remediation is achieved.

Abstract: One of the most challenging and recurring problems when modeling plasmas is the lack of data on the key atomic and molecular reactions that drive plasma processes. Even when there are data for some reactions, complete and validated datasets of chemistries are rarely available. This hinders research on plasma processes and curbs development of industrial applications. The QDB project aims to address this problem by providing a platform for provision, exchange, and validation of chemistry datasets. A new data model developed for QDB is presented. QDB collates published data on both electron scattering and heavy-particle reactions. These data are formed into reaction sets, which are then validated against experimental data where possible. This process produces both complete chemistry sets and identifies key reactions that are currently unreported in the literature. Gaps in the datasets can be filled using established theoretical methods. Initial validated chemistry sets for SF6/CF4/O2 and SF6/CF4/N2/H2 are presented as examples.

Abstract: Fundamental and applied research on plasma‐treated liquids for biomedical applications was boosted in the last few years, dictated by their advantages with respect to direct treatments. However, often, the lack of consistent analysis at a molecular level of these liquids, and of the processes used to produce them, have raised doubts of their usefulness in the clinic. The aim of this article is to critically discuss some basic aspects related to the use of plasma‐treated liquids in medicine, with a focus on their chemical composition. We analyze the main liquids used in the field, how they are affected by non‐thermal plasmas, and the possibility to replicate them without plasma treatment.

Abstract: We present a 3D and 2D Cartesian quasi-neutral plasma model for a low current argon gliding arc discharge, including strong interactions between the gas flow and arc plasma column.
The 3D model is applied only for a short time of 0.2 ms due to its huge computational cost. It mainly serves to verify the reliability of the 2D model. As the results in 2D compare well with those in 3D, they can be used for a better understanding of the gliding arc basic characteristics. More specifically, we investigate the back-breakdown phenomenon induced by an artificially controlled plasma channel, and we discuss its effect on the gliding arc characteristics. The
back-breakdown phenomenon, or backward-jump motion of the arc, as observed in the experiments, results in a drop of the gas temperature, as well as in a delay of the arc velocity with respect to the gas flow velocity, allowing more gas to pass through the arc, and thus increasing the efficiency of the gliding arc for gas treatment applications.