Abstract: Plasma catalysis continues to gain increasing scientific interest, both in established fields like toxic waste abatement and emerging fields like greenhouse gas conversion into value-added chemicals. Attention is typically focused on the obtained conversion process selectivity, rates and energy efficiency. Much less attention is usually paid to the underlying mechanistic aspects of the processes that occur. In this contribution, we critically examine a number of fundamentally important nanoscale thermodynamic aspects of plasma catalysis, which are very relevant to these processes but so far have been overlooked or insufficiently covered in the plasma catalysis literature.

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: Recently, Sr-doped LaCrO3 has been experimentally introduced as a new p-type transparent conducting oxide. It is demonstrated that substituting Sr for La results in inducing p-type conductivity in LaCrO3. Performing first principles calculations we study the electronic structure and formation energy of various point defects in LaCrO3. Our results for the formation energies show that in addition to Sr, two more divalent defects, Ca and Ba, substituting for La in LaCrO3, behave as shallow acceptors in line with previous experimental reports. We further demonstrate that under oxygen-poor growth conditions, these shallow acceptors will be compensated by intrinsic donor-like defects (an oxygen vacancy and Cr on an oxygen site), but in the oxygen-rich growth regime the shallow acceptors have the lowest formation energies between all considered defects and will lead to p-type conductivity.

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: Cold atmospheric pressure plasmas have proven to provide an alternative treatment of cancer by targeting tumorous cells while leaving their healthy counterparts unharmed. However, the underlying mechanisms of the plasma–cell interactions are not yet fully understood. Reactive oxygen species, and in particular hydroxyl radicals (OH), are known to play a crucial role in plasma driven apoptosis of
malignant cells. In this paper we investigate the interaction of OH radicals, as well as H2O2 molecules and HO2 radicals, with DNA by means of reactive molecular dynamics simulations using the ReaxFF force field. Our results provide atomic-scale insight into the dynamics of oxidative stress on DNA caused by the OH radicals, while H2O2 molecules appear not reactive within the considered timescale. Among the observed processes are the formation of 8-OH-adduct radicals, forming the first stages towards the formation of 8-oxoGua and 8-oxoAde, H-abstraction reactions of the amines, and the partial opening of loose DNA ends in aqueous solution.

Abstract: Understanding the nature and effect of the multitude of plasma–surface interactions in plasma catalysis is a crucial requirement for further process development and improvement. A particularly intriguing and rather unique property of a plasma-catalytic setup is the ability of the plasma to modify the electronic structure, and hence chemical properties, of the catalyst through charging, i.e. the absorption of excess electrons. In this work, we develop a quantum chemical model based on density functional theory to study excess negative surface charges in a heterogeneous catalyst exposed to a plasma. This method is specifically applied to investigate plasma-catalytic CO2 activation on supported M/Al2O3 (M=Ti, Ni, Cu) single atom catalysts. We find that (1) the presence of a negative surface charge dramatically improves the reductive power of the catalyst, strongly promoting the splitting of CO2 to CO and oxygen, and (2) the relative activity of the investigated transition metals is also changed upon charging, suggesting that controlled surface charging is a powerful additional parameter to tune catalyst activity and selectivity. These results strongly point to plasma-induced surface charging of the catalyst as an important factor contributing to the plasma-catalyst synergistic effects frequently reported for plasma catalysis.

Abstract: Uniform acceptance force biased Monte Carlo (UFMC) simulations have previously been shown to be a powerful tool to simulate atomic scale processes, enabling one to follow the dynamical path during the simulation. In this contribution, we present a simple proof to demonstrate that this uniform acceptance still complies with the condition of detailed balance, on the condition that the characteristic parameter lambda = 1/2 and that the maximum allowed step size is chosen to be sufficiently small. Furthermore, the relation to Metropolis Monte Carlo (MMC) is also established, and it is shown that UFMC reduces to MMC by choosing the characteristic parameter lambda = 0 [Rao, M. et al. Mol. Phys. 1979, 37, 1773]. Finally, a simple example compares the UFMC and MMC methods.

Abstract: Using reactive molecular dynamics simulations, the melting behavior of nickelcarbon 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 nickelcarbon 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: Controlling and steering the growth of single walled carbon nanotubes is often believed to require controlling of the nucleation stage. Yet, little is known about the microscopic mechanisms governing the nucleation from hydrocarbon molecules. Specifically, we address here the dehydrogenation of hydrocarbon molecules and the formation of all-carbon graphitic islands on metallic nanoclusters from hydrocarbon molecules under conditions typical for carbon nanotube growth. Employing reactive molecular dynamics simulations, we demonstrate for the first time that the formation of a graphitic network occurs through the intermediate formation of vertically oriented, not fully dehydrogenated graphitic islands. Upon dehydrogenation of these vertical graphenes, the islands curve over the surface, thereby forming a carbon network covering the nanoparticle. The results indicate that controlling the extent of dehydrogenation offers an additional parameter to control the nucleation of carbon nanotubes.

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: Methane reforming by plasma catalysis is a complex process that is far from understood. It requires a multidisciplinary approach which ideally takes into account all effects from the plasma on the catalyst, and vice versa. In this contribution, we focus on the interactions of CHx (x = {1,2,3}) radicals that are created in the plasma with several nickel catalyst surfaces. To this end, we perform reactive molecular dynamics simulations, based on the ReaxFF potential, in a wide temperature range of 4001600 K. First, we focus on the H2 formation as a function of temperature and surface structure. We observe that substantial H2 formation is obtained at 1400 K and above, while the role of the surface structure seems limited. Indeed, in the initial stage, the type of nickel surface influences the CH bond breaking efficiency of adsorbed radicals; however, the continuous carbon diffusion into the surface gradually diminishes the surface crystallinity and therefore reduces the effect of surface structure on the H2 formation probability. Furthermore, we have also investigated to what extent the species adsorbed on the catalyst surface can participate in surface reactions more in general, for the various surface structures and as a function of temperature. These results are part of the ongoing research on the methane reforming by plasma catalysis, a highly interesting yet complex alternative to conventional reforming processes.

Abstract: We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined depositionsegregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.

Abstract: tWe investigated the plasma-assisted catalytic reactions for the production of value-added chemicalsfrom Ni-catalyzed plasma dry reforming of methane by means of density functional theory (DFT). Weinspected many activation barriers, from the early stage of adsorption of the major chemical fragmentsderived fromCH4andCO2molecules up to the formation of value-added chemicals at the surface, focusingon the formation of methanol, as well as the hydrogenation of C1and C2hydrocarbon fragments. Theactivation barrier calculations show that the presence of surface-bound H atoms and in some cases alsoremaining chemical fragments at the surface facilitates the formation of products. This implies that thehydrogenation of a chemical fragment on the hydrogenated crystalline surface is energetically favouredcompared to the simple hydrogenation of the chemical fragment at the bare Ni(111) surface. Indeed, thepresence of hydrogen modifies the electronic structure of the surface and the course of the reactions.We therefore conclude that surface-bound H atoms, and to some extent also the remaining chemicalfragments at the crystalline surface, induce the following effects: they facilitate associative desorption ofmethanol and ethane by increasing the rate of H-transfer to the adsorbed fragments while they impedehydrogenation of ethylene to ethane, thus promoting again the desorption of ethylene. Overall, they thusfacilitate the catalytic conversion of the formed fragments from CH4and CO2, into value-added chemicals.Finally, we believe that the retention of methane fragments, especially CH3, in the presence of surface-boundHatoms (as observed here for Ni) can be regarded as an identifier for the proper choice of a catalystfor the production of value-added chemicals.

Abstract: The exact structure and properties of the Si vertical bar SiO2 interface are very important in microelectronics and photovoltaic devices such as metal-oxide-semiconductor field-effect transistors (MOSFETs) and solar cells. Whereas Si vertical bar SiO2 structures are traditionally produced by thermal oxidation, hyperthermal oxidation shows a number of promising advantages. However, the Si vertical bar SiO2 interface induced in hyperthermal Si oxidation has not been properly investigated yet. Therefore, in this work, the interface morphology and interfacial stresses during hyperthermal oxidation at room temperature are studied using reactive molecular dynamics simulations based on the ReaxFF potential. Interface thickness and roughness, as well as the bond length and bond angle distributions in the interface are discussed and compared with other models developed for the interfaces induced by traditional thermal oxidation. The formation of a compressive stress is observed. This compressive stress, which at the interface amounts about 2 GPa, significantly slows down the inward silica growth. This value is close to the experimental value in the Si vertical bar SiO2 interface obtained in traditional thermal oxidation.

Abstract: In recent years, the ability of cold atmospheric pressure plasmas (CAPS) to selectively induce cell death in cancer cells has been widely established. This selectivity has been assigned to the reactive oxygen and nitrogen species (RONS) created in CAPs. To provide new insights in the search for an explanation
for the observed selectivity, we calculate the transfer free energy of multiple ROS across membranes containing a varying amount of cholesterol. The cholesterol fraction is investigated as a selectivity parameter because membranes of cancer cells are known to contain lower fractions of cholesterol compared to healthy cells. We find that cholesterol has a significant effect on the permeation of
reactive species across a membrane. Indeed, depending on the specific reactive species, an increasing cholesterol fraction can lead to (i) an increase of the transfer free energy barrier height and width, (ii) the formation of a local free energy minimum in the center of the membrane and (iii) the creation of extra free energy barriers due to the bulky sterol rings. In the context of plasma oncology, these observations suggest that the increased ingress of RONS in cancer cells can be explained by the decreased cholesterol fraction of their cell membrane.