Abstract: Cold atmospheric plasma (CAP) and plasma-treated liquids (PTLs) have recently become a promising option for cancer treatment, but the underlying mechanisms of the anti-cancer effect are still to a large extent unknown. Although hydrogen peroxide () has been recognized as the major anti-cancer agent of PTL and may enable selectivity in a certain concentration regime, the co-existence of nitrite can create a synergistic effect. We develop a mathematical model to describe the key species and features of the cellular response toward PTL. From the numerical solutions, we define a number of dependent variables, which represent feasible measures to quantify cell susceptibility in terms of the membrane diffusion rate constant and the intracellular catalase concentration. For each of these dependent variables, we investigate the regimes of selective versus non-selective, and of synergistic versus non-synergistic effect to evaluate their potential role as a measure of cell susceptibility. Our results suggest that the maximal intracellular concentration, which in the selective regime is almost four times greater for the most susceptible cells compared to the most resistant cells, could be used to quantify the cell susceptibility toward exogenous . We believe our theoretical approach brings novelty to the field of plasma oncology, and more broadly, to the field of redox biology, by proposing new ways to quantify the selective and synergistic anti-cancer effect of PTL in terms of inherent cell features.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

Abstract: Ammonia is an industrial large-volume chemical, with its main application in fertilizer production. It also attracts increasing attention as a green-energy vector. Over the past century, ammonia production has been dominated by the Haber-Bosch process, in which a mixture of nitrogen and hydrogen gas is converted to ammonia at high temperatures and pressures. Haber-Bosch processes with natural gas as the source of hydrogen are responsible for a significant share of the global CO(2)emissions. Processes involving plasma are currently being investigated as an alternative for decentralized ammonia production powered by renewable energy sources. In this work, we present the PNOCRA process (plasma nitrogen oxidation and catalytic reduction to ammonia), combining plasma-assisted nitrogen oxidation and lean NO(x)trap technology, adopted from diesel-engine exhaust gas aftertreatment technology. PNOCRA achieves an energy requirement of 4.6 MJ mol(-1)NH(3), which is more than four times less than the state-of-the-art plasma-enabled ammonia synthesis from N(2)and H(2)with reasonable yield (>1 %).

Abstract: Numerous applications have required the study of CO2 plasmas since the 1960s, from CO2 lasers to spacecraft heat shields. However, in recent years, intense research activities on the subject have restarted because of environmental problems associated with CO2 emissions. The present review provides a synthesis of the current state of knowledge on the physical chemistry of cold CO2 plasmas. In particular, the different modeling approaches implemented to address specific aspects of CO2 plasmas are presented. Throughout the paper, the importance of conducting joint experimental, theoretical and modeling studies to elucidate the complex couplings at play in CO2 plasmas is emphasized. Therefore, the experimental data that are likely to bring relevant constraints to the different modeling approaches are first reviewed. Second, the calculation of some key elementary processes obtained with semi-empirical, classical and quantum methods is presented. In order to describe the electron kinetics, the latest coherent sets of cross section satisfying the constraints of “electron swarm” analyses are introduced, and the need for self-consistent calculations for determining accurate electron energy distribution function (EEDF) is evidenced. The main findings of the latest zero-dimensional (0D) global models about the complex chemistry of CO2 and its dissociation products in different plasma discharges are then given, and full state-to-state (STS) models of only the vibrational-dissociation kinetics developed for studies of spacecraft shields are described. Finally, two important points for all applications using CO2 containing plasma are discussed: the role of surfaces in contact with the plasma, and the need for 2D/3D models to capture the main features of complex reactor geometries including effects induced by fluid dynamics on the plasma properties. In addition to bringing together the latest advances in the description of CO2 non-equilibrium plasmas, the results presented here also highlight the fundamental data that are still missing and the possible routes that still need to be investigated.

Abstract: The vast biomedical potential of cold atmospheric pressure plasmas (CAPs) is governed by the formation of reactive species. These biologically active species are formed upon the interaction of CAPs with the surroundings. In biological milieu, water plays an essential role. The development of biomedical CAPs thus requires understanding of the sources of the reactive species in aqueous media exposed to the plasma. This is especially important in case of the COST RF plasma jet, which is developed as a reference microplasma system. In this work, we investigated the formation of the OH radicals, H atoms and H2O2 in aqueous solutions exposed to the COST plasma jet. This was done by combining experimental and modelling approaches. The liquid phase species were analysed using UV-Vis spectroscopy and spin trapping with hydrogen isotopes and electron paramagnetic resonance (EPR) spectroscopy. The discrimination between the species formed from the liquid phase and the gas phase molecules was performed by EPR and 1H-NMR analyses of the liquid samples. The concentrations of the reactive species in the gas phase plasma were obtained using a zero-dimensional (0D) chemical kinetics computational model. A three-dimensional (3D) fluid dynamics model was developed to provide information on the induced humidity in the plasma effluent. The comparison of the experimentally obtained trends for the formation of the species as a function of the feed gas and effluent humidity with the modelling results suggest that all reactive species detected in our system are mostly formed in the gas phase plasma inside the COST jet, with minor amounts arising from the plasma effluent humidity.

Abstract: Cold atmospheric plasma in contact with solutions has many applications, but its chemistry contains many unknowns such as the undescribed reactions with solutes. By combining experiments and modelling, we report the first direct demonstration of the reaction of chloride with oxygen atoms in aqueous solutions exposed to cold plasma.

Abstract: Plasma-treated liquids have great potential for biomedical applications. However, insight into the underlying mechanisms and the exact chemistry is still scarce. In this study, we present the combination of a 0D chemical kinetics and a 2D fluid dynamics model to investigate the plasma treatment of a buffered water solution with the kINPen (R) plasma jet. Using this model, we calculated the gas and liquid flow profiles and the transport and chemistry of all species in the gas and the liquid phase. Moreover, we evaluated the stability of the reactive oxygen and nitrogen species after plasma treatment. We found that of all species, only H2O2, HNO2/NO2-, and HNO3/NO3- are stable in the buffered solution after plasma treatment. This is because both their production and loss processes in the liquid phase are dependent on short-lived radicals (e.g. OH, NO, and NO2). Apart from some discrepancy in the absolute values of the concentrations, which can be explained by the model, all general trends and observations in our model are in qualitative agreement with experimental data and literature.

Abstract: Synchrotron small angle X-ray scattering (SAXS) was used to characterize silicalite-1 zeolite crystallization from TEOS/TPAOH/water clear sol. SAXS patterns were recorded over a broad range of length scales, enabling the simultaneous monitoring of nanoparticles and crystals occurring at various stages of the synthesis. A simple two-population model accurately described the patterns. Nanoparticles were modeled by polydisperse coreshell spheres and crystals by monodisperse oblate ellipsoids. These models were consistent with TEM images. The SAXS results, in conjunction with in situ light scattering, showed that nucleation of crystals occurred in a short period of time. Crystals were uniform in size and shape and became increasingly anisotropic during growth. In the presence of nanoparticles, crystal growth was fast. During crystal growth, the number of nanoparticles decreased gradually but their size was constant. These observations suggested that the nanoparticles were growth units in an aggregative crystal growth mechanism. Crystals grown in the presence of nanoparticles developed a faceted habit and intergrowths. In the final stages of growth, nanoparticles were depleted. Concurrently, the crystal growth rate decreased significantly.

Abstract: Ammonia is one of the most produced chemicals, mainly synthesized from fossil fuels for fertilizer applications. Furthermore, ammonia may be one of the energy carriers of the future, when it is produced from renewable electricity. This has spurred research on alternative technologies for green ammonia production. Research on plasma-driven ammonia synthesis has recently gained traction in academic literature. In the current review, we summarize the literature on plasma-driven ammonia synthesis. We distinguish between mechanisms for ammonia synthesis in the presence of a plasma, with and without a catalyst, for different plasma conditions. Strategies for catalyst design are discussed, as well as the current understanding regarding the potential plasma-catalyst synergies as function of the plasma conditions and their implications on energy efficiency. Finally, we discuss the limitations in currently reported models and experiments, as an outlook for research opportunities for further unravelling the complexities of plasma-catalytic ammonia synthesis, in order to bridge the gap between the currently reported models and experimental results.

Abstract: The fast growing world population demands food to survive, and nitrogen-based fertilizers are essential to ensure sufficient food production. Today, fertilizers are mainly produced from non-sustainable fossil fuels<italic>via</italic>the Haber–Bosch process, leading to serious environmental problems. We propose here a novel rotating gliding arc plasma, operating in air, for direct NO_xproduction, which can yield high nitrogen content organic fertilizers without pollution associated with ammonia emission. We explored the efficiency of NO_xproduction in a wide range of feed gas ratios, and for two arc modes: rotating and steady. When the arc is in steady mode, record-value NO_xconcentrations up to 5.5% are achieved which are 1.7 times higher than the maximum concentration obtained by the rotating arc mode, and with an energy consumption of 2.5 MJ mol⁻¹(or<italic>ca.</italic>50 kW h kN⁻¹);<italic>i.e.</italic>the lowest value so far achieved by atmospheric pressure plasma reactors. Computer modelling, using a combination of five different complementary approaches, provides a comprehensive picture of NO_xformation in both arc modes; in particular, the higher NO_xproduction in the steady arc mode is due to the combined thermal and vibrationally-promoted Zeldovich mechanisms.

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: Nitrogen oxyanions and oxyacids are important agents in atmospheric chemistry and medical biology. Although their chemical behavior in solution is relatively well understood, they may behave very differently at the water/air interface of atmospheric aerosols or at the membrane/water interface of cells. Here, we developed a fully classical model for molecular dynamics simulations of NO3−, NO2−, HNO3, and HNO2 in the framework of the GROMOS 53A6 and 54A7 force field versions. The model successfully accounted for the poorly structured solvation shell and ion pairing tendency of NO3−. Accurate pure-liquid properties and hydration free energies were obtained for the oxyacids. Simulations at the water/air interface showed a local enrichment of HNO3 and depletion of NO3−. The effect was discussed in light of earlier spectroscopic data and ab initio calculations, suggesting that HNO3 behaves as a weaker acid at the surface of water. Our model will hopefully allow for efficient and accurate simulations of nitrogen oxyanions and oxyacids in solution and at microheterogeneous interface environments.

Abstract: In this paper, we present the results of combined molecular dynamics−Metropolis Monte Carlo (MD-MMC) simulations of hydrocarbon species at flat diamond (100)2 × 1 and (111)1 × 1 surfaces. The investigated species are considered to be the most important growth species for (ultra)nanocrystalline diamond ((U)NCD) growth. When applying the MMC algorithm to stuck species at monoradical sites, bonding changes are only seen for CH2. The sequence of the bond breaking and formation as put forward by the MMC simulations mimics the insertion of CH2 into a surface dimer as proposed in the standard growth model of diamond. For hydrocarbon species attached to two adjacent radical (biradical) sites, the MMC simulations give rise to significant changes in the bonding structure. For UNCD, the combinations of C3 and C3H2, and C3 and C4H2 (at diamond (100)2 × 1) and C and C2H2 (at diamond (111)1 × 1) are the most successful in nucleating new crystal layers. For NCD, the following combinations pursue the diamond structure the best: C2H2 and C3H2 (at diamond (100)2 × 1) and CH2 and C2H2 (at diamond (111)1 × 1). The different behaviors of the hydrocarbon species at the two diamond surfaces are related to the different sterical hindrances at the diamond surfaces.

Abstract: The initial stage of nanoparticle formation and growth in radiofrequency acetylene (C2H2) plasmas is investigated by means of a self-consistent one-dimensional fluid model. A detailed chemical kinetic scheme, containing electron impact, ion-neutral, and neutral-neutral reactions, has been developed in order to predict the underlying dust growth mechanisms and the most important dust precursors. The model considers 41 different species (neutrals, radicals, ions, and electrons) describing hydrocarbons (CnHm) containing up to 12 carbon atoms. Possible routes for particle growth are discussed. Both positive and negative ion reaction pathways are considered, as consecutive anion- and cation-molecule reactions seem to lead to a fast build up of the carbon skeleton.

Abstract: The behavior of nanoparticles in a low-pressure silane discharge is studied with the use of a self-consistent one-dimensional fluid model. Nanoparticles of a given (prescribed) radius are formed in the discharge by the incorporation of a dust growth mechanism, i.e., by including a step in which large anions (typically Si12H−25), produced in successive chemical reactions of anions with silane molecules, are transformed into particles. Typically a few thousand anions are used for one nanoparticle. The resulting particle density and the charge on the particles are calculated with an iterative method. While the spatial distribution and the charge of the particles are influenced by the plasma, the presence of the nanoparticles will in turn influence the plasma properties. Several simulations with different particle radii are performed. The resulting density profile of the dust will greatly depend on the particle size, as it reacts to the shift of the balance of the different forces acting on the particles.