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: 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: Cold atmospheric pressure plasmas are gaining increased interest in the medical sector and clinical trials to treat skin diseases are underway. Plasmas are capable of producing several reactive oxygen and nitrogen species (RONS). However, there are open questions how plasma-generated RONS interact on a molecular level in a biological environment, e.g. cells or cell components. The redox pair glutathione (GSH) and glutathione disulphide (GSSG) forms the most important redox buffer in organisms responsible for detoxification of intracellular reactive species. We apply Raman spectroscopy, mass spectrometry, and molecular dynamics simulations to identify the time-dependent chemical modifications on GSH and GSSG that are caused by dielectric barrier discharge under ambient conditions. We find GSSG, S-oxidised glutathione species, and S-nitrosoglutathione as oxidation products with the latter two being the final products, while glutathione sulphenic acid, glutathione sulphinic acid, and GSSG are rather reaction intermediates. Experiments using stabilized pH conditions revealed the same main oxidation products as were found in unbuffered solution, indicating that the dominant oxidative or nitrosative reactions are not influenced by acidic pH. For more complex systems these results indicate that too long treatment times can cause difficult-to-handle modifications to the cellular redox buffer which can impair proper cellular function.

Abstract: The diverse pharmacological properties of the diaryltriazenes have sparked the interest to investigate their potential to be repurposed as antitubercular drug candidates. In an attempt to improve the antitubercular activity of a previously constructed diaryltriazene library, eight new halogenated nitroaromatic triazenides were synthesized and underwent biological evaluation. The potency of the series was confirmed against the Mycobacterium tuberculosis lab strain H37Ra, and for the most potent derivative, we observed a minimal inhibitory concentration of 0.85 μm. The potency of the triazenide derivatives against M. tuberculosis H37Ra was found to be highly dependent on the nature of the halogenated phenyl substituent and less dependent on cationic species used for the preparation of the salts. Although the inhibitory concentration against J774A.1 macrophages was observed at 3.08 μm, the cellular toxicity was not mediated by the generation of nitroxide intermediate as confirmed by electron paramagnetic resonance spectroscopy, whereas no in vitro mutagenicity could be observed for the new halogenated nitroaromatic triazenides when a trifluoromethyl substituent was present on both the aryl moieties.

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: Journal of Physics D: Applied Physics published the first Plasma Roadmap in 2012
consisting of the individual perspectives of 16 leading experts in the various sub-fields of low temperature plasma science and technology. The 2017 Plasma Roadmap is the first update of a planned series of periodic updates of the Plasma Roadmap. The continuously growing interdisciplinary nature of the low temperature plasma field and its equally broad range of applications are making it increasingly difficult to identify major challenges that encompass all of the many sub-fields and applications. This intellectual diversity is ultimately a strength of the field. The current state of the art for the 19 sub-fields addressed in this roadmap demonstrates the enviable track record of the low temperature plasma field in the development of plasmas as an enabling technology for a vast range of technologies that underpin our modern society. At the same time, the many important scientific and technological challenges shared in this roadmap show that the path forward is not only scientifically rich but has the potential to make wide and far reaching contributions to many societal challenges.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: above atmospheric) pressure regimes (1–3.5 bar). It was demonstrated that these operational conditions significantly influence both the discharge dynamics and the process efficiencies of O2 and CO2 discharges. For the case of the O2 DBD, the pressure rise results in the amplification of the discharge current, the appearance of emission lines of the metal electrode material (Fe, Cr, Ni) in the optical emission spectrum and the formation of a granular film of the erosion products (10–300 nm iron oxide nanoparticles) on the reactor walls. Somewhat similar behavior was observed also for the CO2 DBD. The discharge current, the relative intensity of the CO Angstrom band measured by Optical Emission Spectroscopy (OES) and the CO2 conversion rates could be stimulated to some extent by the rise in pressure. The optimal conditions for the O2 DBD (P = 2 bar) and the CO2 DBD (P = 1.5 bar) are demonstrated. It can be argued that the dynamics of the microdischarges (MD) define the underlying process of this behavior. It could be
demonstrated that the pressure increase stimulates the formation of more intensive but fewer MDs. In this way, the operating pressure can represent an additional tool to manipulate the properties of the MDs in a DBD, and as a result also the discharge performance.

Abstract: Background Excessive maternal pre-pregnancy and gestational weight gain are related to pregnancy- and birth outcomes. The interpregnancy time window offers a unique opportunity to intervene in order to acquire a healthy lifestyle before the start of a new pregnancy. Methods INTER-ACT is an e-health driven multicentre randomised controlled intervention trial targeting women at high risk of pregnancy- and birth related complications. Eligible women are recruited for the study at day 2 or 3 postpartum. At week 6 postpartum, participants are randomised into the intervention or control arm of the study. The intervention focuses on weight, diet, physical activity and mental well-being, and comprises face-to-face coaching, in which behavioural change techniques are central, and use of a mobile application, which is Bluetooth-connected to a weighing scale and activity tracker. The intervention is rolled out postpartum (4 coaching sessions between week 6 and month 6) and in a new pregnancy (3 coaching sessions, one in each trimester of pregnancy); the mobile app is used throughout the two intervention phases. Data collection includes data from the medical record of the participants (pregnancy outcomes and medical history), anthropometric data (height, weight, waist- and hip circumferences, skinfold thickness and body composition by bio-electrical impedance analysis), data from the mobile app (physical activity and weight; intervention group only) and questionnaires (socio-demographics, breastfeeding, food intake, physical activity, lifestyle, psychosocial factors and process evaluation). Medical record data are collected at inclusion and at delivery of the subsequent pregnancy. All other data are collected at week 6 and month 6 postpartum and every subsequent 6 months until a new pregnancy, and in every trimester in the new pregnancy. Primary outcome is the composite endpoint score of pregnancy-induced hypertension, gestational diabetes mellitus, caesarean section, and large-for-gestational-age infant in the subsequent pregnancy.

Abstract: In recent years there has been growing interest in the use of plasma technology for CO2 conversion. To improve this application, a good insight into the underlying mechanisms is of great importance. This can be obtained from modeling the detailed plasma chemistry in order to understand the chemical reaction pathways leading to CO2 conversion (either in pure form or mixed with another gas). Moreover, in practice, several plasma reactor types are being investigated for CO2 conversion, so in addition it is essential to be able to model these reactor geometries so that their design can be improved, and the most energy efficient CO2 conversion can be achieved. Modeling the detailed plasma chemistry of CO2 conversion in complex reactors is, however, very time-consuming. This problem can be overcome by using a combination of two different types of model: 0D chemical reaction kinetics models are very suitable for describing the detailed plasma chemistry, while the characteristic features of different reactor geometries can be studied by 2D or 3D fluid models. In the first instance the latter can be developed in argon or helium with a simple chemistry to limit the calculation time; however, the ultimate aim is to implement the more complex CO2 chemistry in these models. In the present paper, examples will be given of both the 0D plasma chemistry models and the 2D and 3D fluid models for the most common plasma reactors used for CO2 conversion in order to emphasize the complementarity of both approaches. Furthermore, based on the modeling insights, the paper discusses the possibilities and limitations of plasma-based CO2 conversion in different types of plasma reactors, as well as what is needed to make further progress in this field.