Abstract: In this review paper, we present several examples of reactive molecular dynamics simulations, which contribute to a better understanding of the underlying mechanisms in plasma medicine on the atomic scale. This includes the interaction of important reactive oxygen plasma species with the outer cell wall of both gram-positive and gram-negative bacteria, and with lipids present in human skin. Moreover, as most biomolecules are surrounded by a liquid biofilm, the behavior of these plasma species in a liquid (water) layer is presented as well. Finally, a perspective for future atomic scale modeling studies is given, in the field of plasma medicine in general, and for cancer treatment in particular.

Abstract: Melittin (MEL), a small peptide component of bee venom, has been reported to exhibit anti-cancer effects in vitro and in vivo. However, its clinical applicability is disputed because of its non-specific cytotoxicity and haemolytic activity in high treatment doses. Plasma-treated phosphate buffered saline solution (PT-PBS), a solution rich in reactive oxygen and nitrogen species (RONS) can disrupt the cell membrane integrity and induce cancer cell death through oxidative stress-mediated pathways. Thus, PT-PBS could be used in combination with MEL to facilitate its access into cancer cells and to reduce the required therapeutic dose. The aim of our study is to determine the reduction of the effective dose of MEL required to eliminate cancer cells by its combination with PT-PBS. For this purpose, we have optimised the MEL threshold concentration and tested the combined treatment of MEL and PT-PBS on A375 melanoma and MCF7 breast cancer cells, using in vitro, in ovo and in silico approaches. We investigated the cytotoxic effect of MEL and PT-PBS alone and in combination to reveal their synergistic cytological effects. To support the in vitro and in ovo experiments, we showed by computer simulations that plasma-induced oxidation of the phospholipid bilayer leads to a decrease of the free energy barrier for translocation of MEL in comparison with the non-oxidized bilayer, which also suggests a synergistic effect of MEL with plasma induced oxidation. Overall, our findings suggest that MEL in combination with PT-PBS can be a promising combinational therapy to circumvent the non-specific toxicity of MEL, which may help for clinical applicability in the future.

Abstract: Pancreatic ductal adenocarcinoma (PDAC) is one of the most aggressive cancers with a low response to treatment and a five-year survival rate below 5%. The ineffectiveness of treatment is partly because of an immunosuppressive tumor microenvironment, which comprises tumor-supportive pancreatic stellate cells (PSCs). Therefore, new therapeutic strategies are needed to tackle both the immunosuppressive PSC and pancreatic cancer cells (PCCs). Recently, physical cold atmospheric plasma consisting of reactive oxygen and nitrogen species has emerged as a novel treatment option for cancer. In this study, we investigated the cytotoxicity of plasma-treated phosphate-buffered saline (pPBS) using three PSC lines and four PCC lines and examined the immunogenicity of the induced cell death. We observed a decrease in the viability of PSC and PCC after pPBS treatment, with a higher efficacy in the latter. Two PCC lines expressed and released damage-associated molecular patterns characteristic of the induction of immunogenic cell death (ICD). In addition, pPBS-treated PCC were highly phagocytosed by dendritic cells (DCs), resulting in the maturation of DC. This indicates the high potential of pPBS to trigger ICD. In contrast, pPBS induced no ICD in PSC. In general, pPBS treatment of PCCs and PSCs created a more immunostimulatory secretion profile (higher TNF-α and IFN-γ, lower TGF-β) in coculture with DC. Altogether, these data show that plasma treatment via pPBS has the potential to induce ICD in PCCs and to reduce the immunosuppressive tumor microenvironment created by PSCs. Therefore, these data provide a strong experimental basis for further in vivo validation, which might potentially open the way for more successful combination strategies with immunotherapy for PDAC.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Direct oxidation of methane to methanol (DOMTM) is attractive for the increasing industrial demand of feedstock. In this review, the latest advances in heterogeneous catalysis and plasma catalysis for DOMTM are summarized, with the aim to pinpoint the differences between both, and to provide some insights into their reaction mechanisms, as well as the implications for future development of highly selective catalysts for DOMTM.

Abstract: N2 fixation into NH3 is one of the main processes in the chemical industry. Plasma catalysis is among the environmentally friendly alternatives to the industrial energy-intensive Haber-Bosch process. However, many questions remain open, such as the applicability of the conventional catalytic knowledge to plasma. In this work, we studied the performance of Al2O3-supported Fe, Ru, Co and Cu catalysts in plasma-catalytic NH3 synthesis in a DBD reactor. We investigated the effects of different active metals, and different ratios of the feed gas components, on the concentration and production rate of NH3, and the energy consumption of the plasma system. The results show that the trend of the metal activity (common for thermal catalysis) does not appear in the case of plasma catalysis: here, all metals exhibited similar performance. These findings are in good agreement with our recently published microkinetic model. This highlights the virtual independence of NH3 production on the metal catalyst material, thus validating the model and indicating the potential contribution of radical adsorption and Eley-Rideal reactions to the plasma-catalytic mechanism of NH3 synthesis.

Abstract: Plasma-based nitrogen fixation for fertilizer production is an attractive alternative to the fossil fuel-based industrial processes. However, many factors hinder its applicability, e.g., the commonly observed inverse correlation between energy consumption and production rates or the necessity to enhance the selectivity toward NO2, the desired product for a more facile formation of nitrate-based fertilizers. In this work, we investigated the use of a rotating gliding arc plasma for nitrogen fixation at elevated pressures (up to 3 barg), at different feed gas flow rates and composition. Our results demonstrate a dramatic increase in the amount of NOx produced as a function of increasing pressure, with a record-low EC of 1.8 MJ/(mol N) while yielding a high production rate of 69 g/h and a high selectivity (94%) of NO2. We ascribe this improvement to the enhanced thermal Zeldovich mechanism and an increased rate of NO oxidation compared to the back reaction of NO with atomic oxygen, due to the elevated pressure.

Abstract: The current global energy crisis indicated that increasing our
insight into nonfossil fuel nitrogen fixation pathways for synthetic fertilizer
production is more crucial than ever. Nonequilibrium plasma is a good candidate
because it can use N2 or air as a N source and water directly as a H source, instead
of H2 or fossil fuel (CH4). In this work, we investigate NH3 gas phase formation
pathways from humid N2 and especially humid air up to 2.4 mol % H2O (100%
relative humidity at 20 °C) by optical emission spectroscopy and Fouriertransform
infrared spectroscopy. We demonstrate that the nitrogen fixation
capacity is increased when water vapor is added, as this enables HNO2 and NH3
production in both N2 and air. However, we identified a significant loss
mechanism for NH3 and HNO2 that occurs in systems where these species are
synthesized simultaneously; i.e., downstream from the plasma, HNO2 reacts with NH3 to form NH4NO2, which rapidly decomposes
into N2 and H2O. We also discuss approaches to prevent this loss mechanism, as it reduces the effective nitrogen fixation when not
properly addressed and therefore should be considered in future works aimed at optimizing plasma-based N2 fixation. In-line removal
of HNO2 or direct solvation in liquid are two proposed strategies to suppress this loss mechanism. Indeed, using liquid H2O is
beneficial for accumulation of the N2 fixation products. Finally, in humid air, we also produce NH4NO3, from the reaction of HNO3
with NH3, which is of direct interest for fertilizer application.

Abstract: Plasma catalysis is emerging for plasma-assisted gas conversion
processes. However, the underlying mechanisms of plasma catalysis are poorly
understood. In this work, we present a 1D heterogeneous catalysis model with axial
dispersion (i.e., accounting for back-mixing and molecular diffusion of fluid elements in
the process stream in the axial direction), for plasma-catalytic NO production from
N2/O2 mixtures. We investigate the concentration and reaction rates of each species
formed as a function of time and position across the catalyst, in order to determine the
underlying mechanisms. To obtain insights into how the performance of the process
can be further improved, we also study how changes in the postplasma gas flow
composition entering the catalyst bed and in the operation conditions of the catalytic
stage affect the performance of NO production.

Abstract: The behavior of hydrocarbon species at step edges of diamond terraces is investigated by means of combined molecular dynamics−Metropolis Monte Carlo simulations. The results show that the formation of ballas-like diamond films (like UNCD) and well-faceted diamond films (like NCD) can be related to the gas phase concentrations of CxHy in a new manner: Species that have high concentrations above the growing UNCD films suppress the extension of step edges through defect formation. The species that are present above the growing NCD film, however, enhance the extension of diamond terraces, which is believed to result in well-faceted diamond films. Furthermore, it is shown that, during UNCD growth, CxHy species with x ≥ 2 play an important role, in contrast to the currently adopted CVD diamond growth mechanism. Finally, the probabilities for the extension of the diamond (100) terrace are much higher than those for the diamond (111) terrace, which is in full agreement with the experimental observation that diamond (100) facets are more favored than diamond (111) facets during CVD diamond growth.

Abstract: The mechanism of interaction of cold nonequilibrium plasma jets with mammalian cells in physiologic liquid is reported. The major biological active species produced by an argon RF plasma jet responsible for cell viability reduction are analyzed by experimental results obtained through physical, biological, and chemical diagnostics. This is complemented with chemical kinetics modeling of the plasma source to assess the dominant reactive gas phase species. Different plasma chemistries are obtained by changing the feed gas composition of the cold argon based RF plasma jet from argon, humidified argon (0.27%), to argon/oxygen (1%) and argon/air (1%) at constant power. A minimal consensus physiologic liquid was used, providing isotonic and isohydric conditions and nutrients but is devoid of scavengers or serum constituents. While argon and humidified argon plasma led to the creation of hydrogen peroxide dominated action on the mammalian cells, argonoxygen and argonair plasma created a very different biological action and was characterized by trace amounts of hydrogen peroxide only. In particular, for the argonoxygen (1%), the authors observed a strong negative effect on mammalian cell proliferation and metabolism. This effect was distance dependent and showed a half life time of 30 min in a scavenger free physiologic buffer. Neither catalase and mannitol nor superoxide dismutase could rescue the cell proliferation rate. The strong distance dependency of the effect as well as the low water solubility rules out a major role for ozone and singlet oxygen but suggests a dominant role of atomic oxygen. Experimental results suggest that O reacts with chloride, yielding Cl2 − or ClO−. These chlorine species have a limited lifetime under physiologic conditions and therefore show a strong time dependent biological activity. The outcomes are compared with an argon MHz plasma jet (kinpen) to assess the differences between these (at least seemingly) similar plasma sources.

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: Recycling of carbon dioxide by its conversion into value-added products has gained significant interest owing to the role it can play for use in an anthropogenic carbon cycle. The combined conversion with H2O could even mimic the natural photosynthesis process. An interesting gas conversion technique currently being considered in the field of CO2 conversion is plasma technology. To investigate whether it is also promising for this combined conversion, we performed a series of experiments and developed a chemical kinetics plasma chemistry model for a deeper understanding of the process. The main products formed were the syngas components CO and H2, as well as O2 and H2O2, whereas methanol formation was only observed in the parts-per-billion to parts-per-million range. The syngas ratio, on the other hand, could easily be controlled by varying both the water content and/or energy input. On the basis of the model, which was validated with experimental results, a chemical kinetics analysis was performed, which allowed the construction and investigation of the different pathways leading to the observed experimental results and which helped to clarify these results. This approach allowed us to evaluate this technology on the basis of its underlying chemistry and to propose solutions on how to further improve the formation of value-added products by using plasma technology.