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: The interaction between cold atmospheric pressure plasma and liquids is receiving increasing attention for various applications. In particular, the use of plasma-treated liquids (PTL) for biomedical applications is of growing importance, in particular for sterilization and cancer treatment. However, insight into the
underlying mechanisms of plasma–liquid interactions is still scarce. Here, we present a 2D fluid dynamics model for the interaction between a plasma jet and liquid water. Our results indicate that the formed reactive species originate from either the gas phase (with further solvation) or are formed at the liquid interface. A clear increase in the aqueous density of H2O2, HNO2/NO2- and NO3-
is observed as a function of time, while the densities of O3, HO2/O2- and ONOOH/ONOO- are found to quickly reach a maximum due to chemical reactions in solution. The trends observed in our model correlate well with experimental observations from the literature.

Abstract: Amorphous semiconductors are usually characterized by a low charge carrier mobility, essentially related to their lack of long-range order. The development of such material with higher charge carrier mobility is hence challenging. Part of the issue comes from the difficulty encountered by first-principles simulations to evaluate concepts such as the electron effective mass for disordered systems since the absence of periodicity induced by the disorder precludes the use of common concepts derived from condensed matter physics. In this paper, we propose a methodology based on first-principles simulations that partially solves this problem, by quantifying the degree of delocalization of a wave function and of the connectivity between the atomic sites within this electronic state. We validate the robustness of the proposed formalism on crystalline and molecular systems and extend the insights gained to disordered/amorphous InGaZnO4 and Si. We also explore the properties of p-type oxide semiconductor candidates recently reported to have a low effective mass in their crystalline phases [G. Hautier et al., Nat. Commun. 4, 2292 (2013)]. Although in their amorphous phase none of the candidates present a valence band with delocalization properties matching those found in the conduction band of amorphous InGaZnO4, three of the seven analyzed materials show some potential. The most promising candidate, K2Sn2O3, is expected to possess in its amorphous phase a slightly higher hole mobility than the electron mobility in amorphous silicon.

Abstract: In this study, we performed density functional theory (DFT) calculations to investigate different reaction mechanisms of CO oxidation catalyzed by the Si atom embedded defective BC2N nanostructures as well as the analysis of the structural and electronic properties. The structures of all the complexes are optimized and characterized by frequency calculations at the M062X/6-31G* computational level. Also, The electronic structures and thermodynamic parameters of adsorbed CO and O-2 molecules over Si-doped BC2N nanostructures are examined in detail. Moreover, to investigate the curvature effect on the CO oxidation reaction, all the adsorption and CO oxidation reactions on a finite-sized armchair (6,6) Si-BC2NNT are also studied. Our results indicate that there can be two possible pathways for the CO oxidation with O-2 molecule: O-2(g) + CO(g) -> O-2(ads) + CO(ads) -> CO2(g) + O-(ads) and O-(ads) + CO(g) -> CO2(g). The first reaction proceeds via the Langmuir-Hinshelwood (LH) mechanism while the second goes through the Eley-Rideal (ER) mechanism. On the other hand, by increasing the tube diameter, the energy barrier increases due to the strong adsorption energy of the O-2 molecule which is related to its dissociation over the tube surface. Our calculations indicate that the two step energy barrier of the oxidation reaction over Si-BC2NNS is less than that over the Si-BC2NNT. Hence, Si-BC2NNS may serve as an efficient and highly activated substrate to CO oxidation rather than (4,4) Si-BC2NNT. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: A packed bed dielectric barrier discharge (DBD) is widely used for plasma catalysis applications, but the exact plasma characteristics in between the packing beads are far from understood. Therefore, we study here these plasma characteristics by means of fluid modelling and experimental observations using ICCD imaging, for packing materials with different dielectric constants. Our study reveals that a packed bed DBD reactor in dry air at atmospheric pressure may show three types of discharges, i.e. positive restrikes, filamentary microdischarges, which can also be localized between two packing beads, and surface discharges (so-called surface ionization
waves). Restrikes between the dielectric surfaces result in the formation of filamentary microdischarges, while surface charging creates electric field components parallel to the dielectric surfaces, leading to the formation of surface discharges. A transition in discharge mode occurs from surface discharges to local filamentary discharges between the packing beads when the dielectric constant of the packing rises from 5 to 1000. This may have implications for the efficiency of plasma catalytic gas treatment, because the catalyst activation may be limited by constraining the discharge to the contact points of the beads. The production of reactive species occurs most in the positive restrikes, the surface discharges and the local microdischarges in between the beads, and is less significant in the longer filamentary microdischarges. The faster streamer propagation and discharge development with higher dielectric constant of the packing beads leads to a faster production of reactive species. This study is of great interest for plasma catalysis, where packing beads with different dielectric constants are often used as supports for the catalytic materials. It allows us to better understand how different packing materials can influence the performance of packed bed plasma reactors for environmental applications.

Abstract: The gliding arc plasmatron (GAP) is a highly efficient atmospheric plasma source, which is very promising for CO2 conversion applications. To understand its operation principles and to improve its application, we present here comprehensive modeling results, obtained by means of computational fluid dynamics simulations and plasma modeling. Because of the complexity of the CO2 plasma, a full 3D plasma model would be computationally impractical. Therefore, we combine a 3D turbulent gas flow model with a 2D plasma and gas heating model in order to calculate the plasma parameters and CO2 conversion characteristics. In addition, a complete 3D gas flow and plasma model with simplified argon chemistry is used to evaluate the gliding arc evolution in space and time. The calculated values are compared with experimental data from literature as much as possible in order to validate the model. The insights obtained in this study are very helpful for improving the application of CO2 conversion, as they allow us to identify the limiting factors in the performance, based on which solutions can be provided on how to further improve the capabilities of CO2 conversion in the GAP.

Abstract: Synthetic glucocorticoids (GC) are the mainstay therapy for treatment of acute and chronic inflammatory disorders. Due to the high adverse effects associated with long-term use, GC pharmacology has focused since the nineties on more selective GC ligand-binding strategies, classified as selective glucocorticoid receptor (GR) agonists (SEGRAs) or selective glucocorticoid receptor modulators (SEGRMs). In the current study, GSK866 analogs with electrophilic covalent-binding warheads were developed with potential SEGRA properties to improve their clinical safety profile for long-lasting topical skin disease applications. Since the off-rate of a covalently binding drug is negligible compared to that of a non-covalent drug, its therapeutic effects can be prolonged and typically, smaller doses of the drug are necessary to reach the same level of therapeutic efficacy, thereby potentially reducing systemic side effects. Different analogs of SEGRA GSK866 coupled to cysteine reactive warheads were characterized for GR potency and selectivity in various biochemical and cellular assays. GR- and NFκB-dependent reporter gene studies show favorable anti-inflammatory properties with reduced GR transactivation of two non-steroidal GSK866 analogs UAMC-1217 and UAMC-1218, whereas UAMC-1158 and UAMC-1159 compounds failed to modulate cellular GR activity. These results were further supported by GR immuno-localization and S211 phospho-GR western analysis, illustrating significant GR phosphoactivation and nuclear translocation upon treatment of GSK866, UAMC-1217, or UAMC-1218, but not in case of UAMC-1158 or UAMC-1159. Furthermore, mass spectrometry analysis of tryptic peptides of recombinant GR ligand-binding domain (LBD) bound to UAMC-1217 or UAMC-1218 confirmed covalent cysteine-dependent GR binding. Finally, molecular dynamics simulations, as well as glucocorticoid receptor ligand-binding domain (GR-LBD) coregulator interaction profiling of the GR-LBD bound to GSK866 or its covalently binding analogs UAMC-1217 or UAMC-1218 revealed subtle conformational differences that might underlie their SEGRA properties. Altogether, GSK866 analogs UAMC-1217 and UAMC-1218 hold promise as a novel class of covalent-binding SEGRA ligands for the treatment of topical inflammatory skin disorders.

Abstract: By means of chemical kinetics modeling, it is possible to elucidate the main dissociation mechanisms of CO2 in a gliding arc plasmatron (GAP). We obtain good agreement between the calculated and experimental conversions and energy efficiencies, indicating that the model can indeed be used to study the underlying mechanisms. The calculations predict that vibration-induced dissociation is the main dissociation mechanism of CO2, but it occurs mainly from the lowest vibrational levels because of fast thermalization of the vibrational distribution. Based on these findings, we propose ideas for improving the performance of the GAP, but testing of these ideas in the simulations reveals that they do not always lead to significant enhancement, because of other side effects, thus illustrating the complexity of the process. Nevertheless, the model allows more insight into the underlying mechanisms to be obtained and limitations to be identified.

Abstract: Low-temperature plasma chemical kinetic models are particularly important to the plasma community. These models typically require dozens of inputs, especially rate coefficients. The latter are not always precisely known and it is not surprising that the error on the rate coefficient data can propagate to the model output. In this paper, we present a model that uses N = 400 different combinations of rate coefficients based on the uncertainty attributed to each rate coefficient, giving a good estimation of the uncertainty on the model output due to the rate coefficients. We demonstrate that the uncertainty varies a lot with the conditions and the type of output. Relatively low uncertainties (about 15%) are found for electron density and temperature, while the uncertainty can reach more than an order of magnitude for the population of the vibrational levels in some cases and it can rise up to 100% for the CO2 conversion. The reactions that are mostly responsible for the largest uncertainties are identified. We show that the conditions of pressure, gas temperature and power density have a great effect on the uncertainty and on which reactions lead to this uncertainty. In all the cases tested here, while the absolute values may suffer from large uncertainties, the trends observed in previous modeling work are still valid. Finally, in accordance with the work of Turner, a number of ‘good practices’ is recommended.

Abstract: Epitaxially grown ultra-thin Si layers are often used to passivate Ge surfaces in the high-k gate module of (strained) Ge FinFET devices. We use Si4H10 as Si precursor as it enables epitaxial Si growth at temperatures down to 330 degrees C. C-V characteristics of blanket capacitors made on Ge virtual substrates point to the presence of an optimal Si thickness. In case of compressively strained Ge fin structures, the Si growth results in non-uniform and high strain levels in the strained Ge fin. These strain levels have been calculated for different shapes of the Ge fin and in function of the grown Si thickness. The high strain is the driving force for potential (unwanted) Ge surface reflow during the Si deposition. The Ge surface reflow is strongly affected by the strength of the H-passivation during Si-capping and can be avoided by carefully selected process conditions.

Abstract: Low-temperature plasmas (LTPs) have a proven bactericidal activity governed by the generated reactive oxygen and nitrogen species (RONS) that target microbial cell components. However, RONS also interact with biomolecules in the environment. Here we assess the impact of these interactions upon exposure of liquid suspensions with variable organic content to an atmospheric-pressure dielectric barrier discharge plasma jet. Salmonella enterica serovar Typhimurium viability in the suspension was reduced in the absence [e. g., phosphate buffered saline (PBS)], but not in the presence of (high) organic content [Dulbecco's Modified Eagle's Medium (DMEM), DMEM supplemented with foetal calf serum, and Lysogeny Broth]. The reduced viability of LTP-treated bacteria in PBS correlated to a loss of membrane integrity, whereas double-strand DNA breaks could not be detected in treated single cells. The lack of bactericidal activity in solutions with high organic content correlated with a relative decrease of center dot OH and O-3/O-2(a(1)Delta g)/O, and an increase of H2O2 and NO2- in the plasma-treated solutions. These results indicate that the redox reactions of LTP-generated RONS with nontarget biomolecules resulted in a RONS composition with reduced bactericidal activity. Therefore, the chemical composition of the bacterial environment should be considered in the development of LTP for antimicrobial treatment, and may affect other biomedical applications as well.

Abstract: Dry reforming of methane (DRM) in a gliding arc plasmatron is studied for different CH4 fractions in the mixture. The CO2 and CH4 conversions reach their highest values of approximately 18 and 10%, respectively, at 25% CH4 in the gas mixture, corresponding to an overall energy cost of 10 kJ L@1 (or 2.5 eV per molecule) and an energy efficiency of 66%. CO and H2 are the major products, with the formation of smaller fractions of C2Hx (x=2, 4, or 6) compounds and H2O. A chemical kinetics model is used to investigate the underlying chemical processes. The calculated CO2 and CH4 conversion and the energy efficiency are in good agreement with the experimental data. The model calculations reveal that the reaction of CO2 (mainly at vibrationally excited levels) with H radicals is mainly responsible for
the CO2 conversion, especially at higher CH4 fractions in the mixture, which explains why the CO2 conversion increases with increasing CH4 fraction. The main process responsible for CH4 conversion is the reaction with OH radicals. The excellent energy efficiency can be explained by the non-equilibrium character of the plasma, in which the electrons mainly activate the gas molecules, and by the important role of the vibrational kinetics of CO2. The results demonstrate that a gliding arc plasmatron is very promising for DRM.

Abstract: <script type='text/javascript'>document.write(unpmarked('The crystal and defect structure of SnS crystals grown using chemical vapour deposition for application in electronic devices are investigated. The structural analysis shows the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nanometer scale thickness, and much thicker but smaller crystallites. Both show similar Raman response associated with SnS. The structural analysis with transmission electron microscopy shows that the flakes are single crystals of -SnS with [010] normal to the substrate. Parallel with the surface of the flakes, lamellae with varying thickness of a new SnS phase are observed. High-resolution transmission electron microscopy (TEM), scanning transmission electron microscopy (STEM), first-principles simulations (DFT) and nanobeam diffraction (NBD) techniques are employed to characterise this phase in detail. DFT results suggest that the phase is a strain stabilised \u0027 one grown epitaxially on the -SnS crystals. TEM analysis shows that the crystallites are also -SnS with generally the [010] direction orthogonal to the substrate. Contrary to the flakes the crystallites consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. Under high-dose electron irradiation, the SnS structure is reduced and -Sn formed. It is shown that this damage only occurs for SnS in direct contact with SiO2. Lay description SnS is a p-type semiconductor, which has attracted significant interest for electronic devices due to its unique properties, low-toxicity and abundance of Sn in nature. Although in the past it has been most extensively studied as the absorber material in solar cells, it has recently garnered interest for application as a p-type two-dimensional semiconductor in nanoelectronic devices due to its anisotropic layered structure similar to the better known phosphorene. Tin sulphide can take the form of several phases and the electronic properties of the material depend strongly on its crystal structure. It is therefore crucial to study the crystal structure of the material in order to predict the electronic properties and gain insight into the growth mechanism. In this work, SnS crystals deposited using a chemical vapour deposition technique are investigated extensively for their crystal and defect structure using transmission electron microscopy (TEM) and related techniques. We find the presence of two distinct crystal morphologies, that is thin flakes with lateral sizes up to 50 m and nm scale thickness, and much thicker but smaller crystallites. The flakes are single crystals of -SnS and contain lamellae with varying thickness of a different phase which appear to be -SnS at first glance. High-resolution scanning transmission electron microscopy is used to characterise these lamellae where the annular bright field (ABF) mode better reveals the position of the sulphur columns. The sulphur columns in the lamellae are found to be shifted relative to the -SnS structure which indicates the formation of a new phase which is a distorted version of the phase which we tentatively refer to as \u0027-SnS. Simulations based on density functional theory (DFT) are used to model the interface and a similar shift of sulphur columns in the -SnS layer is observed which takes place as a result of strong interaction at the interface between the two phases resulting in strain transfer. Nanobeam electron diffraction (NBD) is used to map the lattice mismatch in the thickness of the flakes which reveals good in-plane matching and some expansion out-of-plane in the lamellae. Contrary to the flakes the crystallites are made solely of -SnS and consist of two to four grains which are tilted up to 15 degrees relative to the substrate. The various grain boundary structures and twin relations are discussed. At high electron doses, SnS is reduced to -Sn, however the damage occurs only for SnS in direct contact with SiO2.'));

Abstract: Studies using xenograft mouse models have shown that plasma applied to the skin overlying tumors results in tumor shrinkage. Plasma is considered a nonpenetrating treatment; however, these studies demonstrate plasma effects that occur beyond the postulated depth of physical penetration of plasma components. The present study examines the propagation of plasma effects through a tissue model using three-dimensional, cell-laden extracellular matrices (ECMs). These ECMs are used as barriers against direct plasma penetration. By placing them onto a monolayer of target cancer cells to create an in-vitro analog to in-vivo studies, we distinguished between cellular effects from direct plasma exposure and cellular effects due to cell-to-cell signaling stimulated by plasma. We show that nanosecond-pulsed dielectric barrier discharge plasma treatment applied atop an acellular barrier impedes the externalization of calreticulin (CRT) in the target cells. In contrast, when a barrier is populated with cells, CRT externalization is restored. Thus, we demonstrate that plasma components stimulate signaling among cells embedded in the barrier to transfer plasma effects to the target cells.