Abstract: In this paper, an overview is given of different modeling approaches used for describing gas discharge plasmas, as well as plasma-surface interactions. A fluid model is illustrated for describing the detailed plasma chemistry in capacitively coupled rf discharges. The strengths and limitations of Monte Carlo simulations and of a particle-in-cell-Monte Carlo collisions model are explained for a magnetron discharge, whereas the capabilities of a hybrid Monte Carlo-fluid approach are illustrated for a direct current glow discharge used for spectrochemical analysis of materials. Finally, some examples of molecular dynamics simulations, for the purpose of plasma-deposition, are given.

Abstract: Reactive molecular dynamics simulations are performed to study the interaction of reactive oxygen species, such as OH, HO2 and H2O2, with the endotoxic biomolecule lipid A of the gram-negative bacterium Escherichia coli. It is found that the aforementioned plasma species can destroy the lipid A, which consequently results in reducing its toxic activity. All bond dissociation events are initiated by hydrogen-abstraction reactions. However, the mechanisms behind these dissociations are dependent on the impinging plasma species, i.e. a clear difference is observed in the mechanisms upon impact of HO2 radicals and H2O2 molecules on one hand and OH radicals on the other hand. Our simulation results are in good agreement with experimental observations.

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: The catalytic oxidative dehydrogenation of hydrocarbons is an industrially important process, in which selectivity is a key issue. We here investigate the conversion of methanol to formaldehyde on a vanadia surface employing long timescale simulations, reaching a time scale of seconds. In particular, we compare the thermal process to the case where an additional external electric field is applied, as would be the case in a direct plasma-catalysis setup. We find that the electric field influences the retention time of the molecules at the catalyst surface. These simulations provide an atomic scale insight in the thermal catalytic oxidative dehydrogenation process, and in how an external electric field may affect this process.

Abstract: We perform molecular dynamics simulations to study the flip-flop motion of phosphatidylserine (PS) across the plasma membrane upon increasing oxidation degree of the membrane. Our computational results show that an increase of the oxidation degree in the lipids leads to a decrease of the free energy barrier for translocation of PS through the membrane. In other words, oxidation of the lipids facilitates PS flip-flop motion across the membrane, because in native phospholipid bilayers this is only a “rare event” due to the high energy barriers for the translocation of PS. The present study provides an atomic-scale insight into the mechanisms of the PS flip-flop upon oxidation of lipids, as produced for example by cold atmospheric plasma, in living cells.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: Purpose Natural deep eutectic solvents (NADES) represent a green alternative to conventional organic solvents as reaction medium, offering more benign properties. To efficiently design NADES for biocatalysis, a better understanding of their effect on these reactions is needed. We hypothesize that this effect can be described by separately considering (1) the solvent interactions with the substrates, (2) the solvent viscosities and (3) the enzyme stability in NADES. Methods We investigated the effect of substrate solvation and viscosity on the reaction rate; and the stability of the enzyme in NADES. To this end, we monitored the conversion over time of the transesterification of vinyl laurate with 1- butanol by the lipase enzyme Candida antarctica B in NADES of different compounds and molar ratios. Results The initial reaction rate is higher in most NADES ( varying between 1.14 and 15.07 mu mol min(-1) mg(-1)) than in the reference n-hexane (4.0 mu mol min(-1) mg(-1))), but no clear relationship between viscosity and initial reaction rate was found. The increased reaction rate is most likely related to the solvation of the substrate due to a change in the activation energy of the reaction or a change in the conformation of the substrate. The enzyme retained part of its activity after the first 2 h of reaction (on average 20 % of the substrate reacted in the 2-24 h period). Enzyme incubation in ethylene glycol-based NADES resulted in a reduced reaction rate ( 15.07 vs. 3.34 mu mol min(-1) mg(-1)), but this may also be due to slow dissolution of the substrate. Conclusions The effect of viscosity seems to be marginal next to the effect of solvation and possible enzyme-NADES interaction. The enzyme retains some of its activity during the 24-hour measurements, but the enzyme incubation experiments did not yield accurate, comparable values. [GRAPHICS] .

Abstract: The limitation in time and length scale is a major issue of molecular dynamics (MD) simulation. Although several methods have been developed to extend the MD time scale, their performance usually deteriorates with increasing system size. Therefore, an acceleration method which is applicable to large systems is required to bridge the gap between the MD simulations and target phenomena. In this study, an accelerated MD method for large system is developed based on the collective variable-driven hyperdynamics (CVHD) method [K.M. Bal and E.C. Neyts, 2015]. The key idea is to run CVHD in parallel with rate control and accelerate multiple possible events simultaneously. Using this novel method, carbon diffusion in bcc-iron bicrystal with grain boundary is examined as an application for practical materials. Carbon atoms reaching at the grain boundary are trapped whereas carbon atoms in the bulk region diffuse randomly, and both dynamic regimes can be simultaneously accelerated with the parallel CVHD technique.

Abstract: Plasma–liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. The review focusses on nonequilibrium plasmas.

Abstract: Plasma generation inside catalyst pores is of utmost importance for plasma catalysis, as the existence of plasma species inside the pores affects the active surface area of the catalyst available to the plasma species for catalytic reactions. In this paper, the electric field enhancement, and thus the plasma production inside catalyst pores with different pore shapes is studied with a two-dimensional fluid model. The results indicate that the electric field will be significantly enhanced near tip-like structures. In a conical pore with small opening, the strongest electric field appears at the opening and bottom corners of the pore, giving rise to a prominent ionization rate throughout the pore. For a cylindrical pore, the electric field is only enhanced at the bottom corners of the pore, with lower absolute value, and thus the ionization rate inside the pore is only slightly enhanced. Finally, in a conical pore with large opening, the electric field is characterized by a maximum at the bottom of the pore, yielding a similar behavior for the ionization rate. These results demonstrate that the shape of the pore has a significantly influence on the electric field enhancement, and thus modifies the plasma properties.

Abstract: The ever increasing global production and dispersion of methane requires novel chemistry to transform it into easily condensable energy carriers that can be integrated into the chemical infrastructure. In this context, single atom catalysts have attracted considerable interest due to their outstanding catalytic activity. We here use density functional theory (DFT) computations to compare the reaction and activation energies of M and MN4 embedded graphene (M = Ni and Si) on the methane-to-methanol conversion near room temperature. Thermodynamically, conversion of methane to methanol is energetically favorable at ambient conditions. Both singlet and triplet spin state of the studied systems are considered in all of the calculations. The DFT results show that the barriers are significantly lower when the complexes are in the triplet state than in the singlet state. In particular, Si-G with the preferred spin multiplicity of triplet seems to be viable catalysts for methane oxidation thanks to the corresponding lower energy barriers and higher stability of the obtained configurations. Our results provide insights into the nature of methane conversion and may serve as guidance for fabricating cost-effective graphene-based single atom catalysts.

Abstract: In this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO2 and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.

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: In this work, we use density functional theory calculations to study the combined effect of external electric fields, surface morphology, and surface charge on CO2 activation over Cu(111), Cu(211), Cu(110), and Cu(001) surfaces. We observe that the binding energy of the CO2 molecule on Cu surfaces increases significantly upon increasing the applied electric field strength. In addition, rougher surfaces respond more effectively to the presence of the external electric field toward facilitating the formation of a carbonate-like CO2 structure and the transformation of the most stable adsorption mode from physisorption to chemisorption. The presence of surface charges further strengthens the electric field effect and consequently causes an improved bending of the CO2 molecule and C−O bond length elongation. On the other hand, a net charge in the absence of an externally applied electric field shows only a marginal effect on CO2 binding. The chemisorbed CO2 is more stable and further activated when the effects of an external electric field, rough surface, and surface charge are combined. These results can help to elucidate the underlying factors that control CO2 activation in heterogeneous and plasma catalysis, as well as in electrochemical processes.