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: In metallized film capacitors, there exists an air gap of about 0.2 μm between the films, with a pressure ranging generally from 130 atm. Because of the created potential difference between the two films, a microdischarge is formed in this gap. In this paper, we use an implicit particle-in-cell Monte Carlo collision simulation method to study the discharge properties in this direct-current microdischarge with 0.2 μm gap in a range of different voltages and pressures. The discharge process is significantly different from a conventional high pressure discharge. Indeed, the high electric field due to the small gap sustains the discharge by field emission. At low applied voltage (~15 V), only the electrons are generated by field emission, while both electrons and ions are generated as a stable glow discharge at medium applied voltage (~50 V). At still higher applied voltage (~100 V), the number of electrons and ions rapidly multiplies, the electric field reverses, and the discharge changes from a glow to an arc regime.

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: A reactive ReaxFF force field is developed for studying the influence of intrinsic point defects on the chemistry with TiO2 condensed phases. The force field parameters are optimized to ab initio data for the equations of state, relative phase stabilities for titanium and titanium dioxide, potential energy differences for (TiO2)n-clusters (n = 116). Also data for intrinsic point defects in anatase were added. These data contain formation energies for interstitial titanium and oxygen vacancies, diffusion barriers of the oxygen vacancies and molecular oxygen adsorption on a reduced anatase (101) surface. Employing the resulting force field, we study the influence of concentration of oxygen vacancies and expansion or compression of an anatase surface on the diffusion of the oxygen vacancies. Also the barrier for oxygen diffusion in the subsurface region is evaluated using this force field. This diffusion barrier of 27.7 kcal/mol indicates that the lateral redistribution of oxygen vacancies on the surface and in the subsurface will be dominated by their diffusion in the subsurface, since both this barrier as well as the barriers for diffusion from the surface to the subsurface and vice versa (17.07 kcal/mol and 21.91 kcal/mol, respectively, as calculated with DFT), are significantly lower than for diffusion on the surface (61.12 kcal/mol as calculated with DFT).

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.

Abstract: Plasma medicine is a rapidly evolving multidisciplinary field at the intersection of chemistry, biochemistry, physics, biology, medicine and bioengineering. It holds great potential in medical, health care, dentistry, surgical, food treatment and other applications. This multidisciplinary nature and variety of possible applications come along with an inherent and intrinsic complexity. Advancing plasma medicine to the stage that it becomes an everyday tool in its respective fields requires a fundamental understanding of the basic processes, which is lacking so far. However, some major advances have already been made through detailed experiments over the last 15 years. Complementary, computer simulations may provide insight that is difficultif not impossibleto obtain through experiments. In this review, we aim to provide an overview of the various simulations that have been carried out in the context of plasma medicine so far, or that are relevant for plasma medicine. We focus our attention mostly on atomistic simulations dealing with plasmabiomolecule interactions. We also provide a perspective and tentative list of opportunities for future modelling studies that are likely to further advance the field.

Abstract: Particle-in-cell/Monte Carlo collision simulations, coupled with an external circuit, are used to investigate the behavior of pulsed dual-frequency (DF) capacitively coupled plasmas (CCPs). It is found that the phase shift θ between the high (or low) frequency source and the pulse modulation has a great influence on the ion density and the ionization rate. By pulsing the high frequency source, the time-averaged ion density shows a maximum when θ = 90∘. The time-averaged ion energy distribution functions (IEDFs) at the driven electrode, however, keep almost unchanged, illustrating the potential of pulsed DF-CCP for independent control of ion density (and flux) and ion energy. A detailed investigation of the temporal evolution of the plasma characteristics indicates that several high frequency harmonics can be excited at the initial stage of a pulse period by tuning the phase shift θ, and this gives rise to strong sheath oscillations, and therefore high ionization rates. For comparison, the pulsing of the low frequency source is also studied. In this case, the ion density changes slightly as a function of time, and the time-averaged ion density shows the same trend as in the HF modulation for different phase shifts θ. Moreover, the time-averaged IEDFs at the driven electrode can be modulated, showing the potential to reduce the maximum ion bombardment energy.

Abstract: Computer simulations based on the particle-in-cell/Monte Carlo collision method are performed to study the plasma characteristics and especially the transition in electron heating mechanisms in a hybrid direct current (dc)/dual-frequency (DF) capacitively coupled CF 4 discharge. When applying a superposed dc voltage, the plasma density first increases, then decreases, and finally increases again, which is in good agreement with experiments. This trend can be explained by the transition between the four main heating modes, i.e., DF coupling, dc and DF coupling, dc source dominant heating, and secondary electron dominant heating.

Abstract: By using first-principles simulations, we investigate the interaction of a 2D silicon layer with two classes of chalcogenide-layered compounds, namely MoX2 and GaX (X = S, Se, Te). A rather weak (van der Waals) interaction between the silicene layers and the chalcogenide layers is predicted. We found that the buckling of the silicene layer is correlated to the lattice mismatch between the silicene layer and the MoX2 or GaX template. The electronic properties of silicene on these different templates largely depend on the buckling of the silicene layer: highly buckled silicene on MoS2 is predicted to be metallic, while low buckled silicene on GaS and GaSe is predicted to be semi-metallic, with preserved Dirac cones at the K points. These results indicate new routes for artificially engineering silicene nanosheets, providing tailored electronic properties of this 2D layer on non-metallic substrates. These non-metallic templates also open the way to the possible integration of silicene in future nanoelectronic devices.

Abstract: In this paper we analyse the gas phase production and loss pathways for several biomedically active species, i.e. N2(A), O, O3, O2(a), N, H, HO2, OH, NO, NO2, N2O5, H2O2, HNO2 and HNO3, in an argon plasma jet flowing into an open humid air atmosphere. For this purpose, we employ a zero-dimensional reaction kinetics model to mimic the typical experimental conditions by fitting several parameters to experimentally measured values. These include ambient air diffusion, the gas temperature profile and power deposition along the jet effluent. We focus in detail on how the pathways of the biomedically active species change as a function of the position in the effluent, i.e. inside the discharge device, active plasma jet effluent and afterglow region far from the nozzle. Moreover, we demonstrate how the reaction kinetics and species production are affected by different ambient air humidities, total deposited power into the plasma and gas temperature along the jet. It is shown that the dominant pathways can drastically change as a function of the distance from the nozzle exit or experimental conditions.

Abstract: Plasma catalysis holds great promise for environmental applications, provided that the process viability can be maximized in terms of energy efficiency and product selectivity. This requires a fundamental understanding of the various processes taking place and especially the mutual interactions between plasma and catalyst. In this review, we therefore first examine the various effects of the plasma on the catalyst and of the catalyst on the plasma that have been described in the literature. Most of these studies are purely experimental. The urgently needed fundamental understanding of the mechanisms underpinning plasma catalysis, however, may also be obtained through modelling and simulation. Therefore, we also provide here an overview of the modelling efforts that have been developed already, on both the atomistic and the macroscale, and we identify the data that can be obtained with these models to illustrate how modelling and simulation may contribute to this field. Last but not least, we also identify future modelling opportunities to obtain a more complete understanding of the various underlying plasma catalytic effects, which is needed to provide a comprehensive picture of plasma catalysis.

Abstract: Ammonia is a crucial nutrient used for plant growth and as a building block in pharmaceutical and chemical industry, produced via nitrogen fixation of the ubiquitous atmospheric N2. Current industrial ammonia production relies heavily on fossil resources, but a lot of work is put into developing non-fossil based pathways. Among these is the use of nonequilibrium plasma. In this work, we investigated water vapor as H source for nitrogen fixation into NH3 by non-equilibrium plasma. The highest selectivity towards NH3 was observed with low amounts of added H2O vapor, but the highest production rate was reached at high H2O vapor.