Abstract: In this paper, a new implementation of the Metropolis Monte Carlo (MMC) algorithm is presented. When combining the MMC model with a molecular dynamics (MD) code, crystal growth by plasma-enhanced chemical vapor deposition can be simulated. As the MD part simulates impacts of growth species onto the surface on a time scale of picoseconds, the MMC algorithm simulates the slower adatom surface processes. The implementation includes a criterion for the selection of atoms that are allowed to be displaced during the simulation, and a criterion of after how many MMC cycles the simulation is stopped. We performed combined MD-MMC simulations for hydrocarbon species that are important for the growth of ultrananocrystalline diamond (UNCD) films at partially hydrogenated diamond surfaces, since this implementation is part of a study of the growth mechanisms of (ultra)nanocrystalline diamond films. Exemplary for adatom arrangements during the growth of UNCD, the adatom surface behavior of C and C2H2 at diamond (111)1 × 1, C and C4H2 at diamond (111)1 × 1 and C3 at diamond (100)2 × 1 has been investigated. For all cases, the diamond crystal structure is pursued under the influence of MMC simulation. Additional longer time-scale MD simulations put forward very similar structures, verifying the MMC algorithm. Nevertheless, the MMC simulation time is typically one order of magnitude shorter than the MD simulation time.

Abstract: technologies have been expanding, and one of the most exciting and rapidly growing
applications is in biology and medicine. Most biomedical studies with DBD plasma systems are performed in vitro, which include cells grown on the surface of plastic well plates, or in vivo, which include animal research models (e.g. mice, pigs). Since many DBD systems use the biological target as the secondary electrode for direct plasma generation and treatment, they are sensitive to the surface properties of the target, and thus can be altered based on the in vitro or in vivo system used. This could consequently affect biological response from plasma treatment. Therefore, in this study, we investigated the DBD plasma behavior both in vitro (i.e. 96-well flat bottom plates, 96-well U-bottom plates, and 24-well flat bottom plates), and in vivo (i.e. mouse skin). Intensified charge coupled device (ICCD) imaging was performed and the plasma discharges were visually distinguishable between the different systems. The geometry of the wells did not affect DBD plasma generation for low application distances (≤ 2 mm), but differentially affected plasma uniformity on the bottom of the well at greater distances. Since DBD plasma treatment in vitro is rarely performed in dry wells for plasma medicine experiments, the effect of well wetness was also investigated. In all in vitro cases, the uniformity of the DBD plasma was affected when comparing wet versus dry wells, with the plasma in the wide-bottom wells appearing the most similar to plasma generated on mouse skin. Interestingly, based on quantification of ICCD images, the DBD plasma intensity per surface area demonstrated an exponential one-phase decay with increasing application distance, regardless of the in vitro or in vivo system. This trend is similar to that of the energy per pulse of plasma, which is used to determine the total plasma treatment energy for biological systems. Optical emission spectroscopy performed on the plasma revealed similar trends in radical species generation between the plastic well plates and mouse skin. Therefore, taken together, DBD plasma intensity per surface area may be a valuable parameter to be used as a simple method for in situ monitoring during biological treatment and active plasma treatment control, which can be applied for in vitro and in vivo systems.

Abstract: Despite recent advances in melanoma treatment, there are still patients that either do not respond or develop resistance. This unresponsiveness and/or acquired resistance to therapy could be explained by the fact that some melanoma cells reside in a dedifferentiated state. Interestingly, this dedifferentiated state is associated with greater sensitivity to ferroptosis, a lipid peroxidation-reliant, iron-dependent form of cell death. Cytoglobin (CYGB) is an iron hexacoordinated globin that is highly enriched in melanocytes and frequently downregulated during melanomagenesis. In this study, we investigated the potential effect of CYGB on the cellular sensitivity towards (1S, 3R)-RAS-selective lethal small molecule (RSL3)-mediated ferroptosis in the G361 melanoma cells with abundant endogenous expression. Our findings show that an increased basal ROS level and higher degree of lipid peroxidation upon RSL3 treatment contribute to the increased sensitivity of CYGB knockdown G361 cells to ferroptosis. Furthermore, transcriptome analysis demonstrates the enrichment of multiple cancer malignancy pathways upon CYGB knockdown, supporting a tumor-suppressive role for CYGB. Remarkably, CYGB knockdown also triggers activation of the NOD-, LRR- and pyrin domain-containing protein 3 (NLRP3) inflammasome and subsequent induction of pyroptosis target genes. Altogether, we show that silencing of CYGB expression modulates cancer therapy sensitivity via regulation of ferroptosis and pyroptosis cell death signaling pathways.

Abstract: Computer simulations have been performed to describe the effect of small admixtures of hydrogen to an argon glow discharge in the Grimm-type configuration. The two-dimensional density profiles of the various plasma species (i.e., electrons, Ar+, ArH+, H+, H2+ and H3+ ions, H atoms and H2 molecules, Ar metastable atoms and sputtered Cu atoms) are presented for 1% H2 added to the argon glow discharge, and the effect of different H2 additions (varying between 0.1 and 10%) on the species densities, the hydrogen dissociation degree, and the sputtering process, are investigated. Finally, the relative contributions of various production and loss processes for the different plasma species are calculated.

Abstract: In this paper, the effects of prolonged bias-enhanced nucleation (prolonged BEN) on the growth mechanisms of diamond are investigated by molecular dynamics (MD) and combined MD-Metropolis Monte Carlo (MD-MMC) simulations. First, cumulative impacts of CxHy+ and Hx+ on an a-C:H/nanodiamond composite were simulated; second, nonconsecutive impacts of the dominant ions were simulated in order to understand the observed phenomena in more detail. As stated in the existing literature, the growth of diamond structures during prolonged BEN is a process that takes place below the surface of the growing film. The investigation of the penetration behavior of CxHy+ and Hx+ species shows that the carbon-containing ions remain trapped within this amorphous phase where they dominate mechanisms like precipitation of sp3 carbon clusters. The H+ ions, however, penetrate into the crystalline phase at high bias voltages (>100 V), destroying the perfect diamond structure. The experimentally measured reduction of grain sizes at high bias voltage, reported in the literature, might thus be related to penetrating H+ ions. Furthermore, the CxHy+ ions are found to be the most efficient sputtering agents, preventing the build up of defective material.

Abstract: In this paper we study dry reforming of methane (DRM) in a gliding arc plasmatron (GAP) in the presence of N₂and O₂. N₂is added to create a stable plasma at equal fractions of CO₂and CH₄, and because emissions from industrial plants typically contain N₂, while O₂is added to enhance the process. We test different gas mixing ratios to evaluate the conversion and energy cost. We obtain conversions between 31 and 52% for CO₂and between 55 and 99% for CH₄, with total energy costs between 3.4 and 5.0 eV per molecule, depending on the gas mixture. This is very competitive when benchmarked with the literature. In addition, we present a chemical kinetics model to obtain deeper insight in the underlying plasma chemistry. This allows determination of the major reaction pathways to convert CO₂and CH₄, in the presence of O₂and N₂, into CO and H₂. We show that N₂assists in the CO₂conversion, but part of the applied energy is also wasted in N₂excitation. Adding O₂enhances the CH₄conversion, and lowers the energy cost, while the CO₂conversion remains constant, and only slightly drops at the highest O₂fractions studied, when CH₄is fully oxidized into CO₂.

Abstract: We have developed a sustainable method to produce NH3 directly from air using a plasma tandem-electrocatalysis system that operates via the N2−NOx−NH3 pathway. To efficiently reduce NO2− to NH3, we propose a novel electrocatalyst consisting of defective N-doped molybdenum sulfide nanosheets on vertical graphene arrays (N-MoS2/VGs). We used a plasma engraving process to form the metallic 1T phase, N doping, and S vacancies in the electrocatalyst simultaneously. Our system exhibited a remarkable NH3 production rate of 7.3 mg h−1 cm−2 at −0.53 V vs RHE, which is almost 100 times higher than the state-of-the-art electrochemical nitrogen reduction reaction and more than double that of other hybrid systems. Moreover, a low energy consumption of only 2.4 MJ molNH3−1 was achieved in this study. Density functional theory calculations revealed that S vacancies and doped N atoms play a dominant role in the selective reduction of NO2− to NH3. This study opens up new avenues for efficient NH3 production using cascade systems.

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: In this paper, an overview is given of modeling activities going on in our research group, for describing the plasma chemistry and plasmasurface interactions in reactive plasmas. The plasma chemistry is calculated by a fluid approach or by hybrid Monte Carlo (MC)fluid modeling. An example of both is illustrated in the first part of the paper. The example of fluid modeling is given for a dielectric barrier discharge (DBD) in CH4/O2, to describe the partial oxidation of CH4 into value-added chemicals. The example of hybrid MCfluid modeling concerns an inductively coupled plasma (ICP) etch reactor in Ar/Cl2/O2, including also the description of the etch process. The second part of the paper deals with the treatment of plasmasurface interactions on the atomic level, with molecular dynamics (MD) simulations or a combination of MD and MC simulations.