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: 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: The unique plasma-specific features and physical phenomena in the organization of nanoscale soild-state systems in a broad range of elemental composition, structure, and dimensionality are critically reviewed. These effects lead to the possibility to localize and control energy and matter at nanoscales and to produce self-organized nano-solids with highly unusual and superior properties. A unifying conceptual framework based on the control of production, transport, and self-organization of precursor species is introduced and a variety of plasma-specific non-equilibrium and kinetics-driven phenomena across the many temporal and spatial scales is explained. When the plasma is localized to micrometer and nanometer dimensions, new emergent phenomena arise. The examples range from semiconducting quantum dots and nanowires, chirality control of single-walled carbon nanotubes, ultra-fine manipulation of graphenes, nano-diamond, and organic matter to nano-plasma effects and nano-plasmas of different states of matter.

Abstract: A seriesof [TiN/TiAlN]nmultilayer coatingswith different bilayer numbers n=5, 10, 25, 50, and 100 were deposited on stainless steel substrate AISI 304 by a lateral rotating cathode arc technique in a flowing nitrogen atmosphere. The composition and microstructure of the coatings have been analyzed by using energy dispersive X-ray spectroscopy, X-ray diffraction (XRD), and conventional and high-resolution transmission electron microscopy (HRTEM). XRD analysis shows that the preferential orientation growth along the (111) direction is reduced in the multilayer coatings. TEM analysis reveals that the grain size of the coatings decreases with increasing bilayer number. HRTEMimaging of the multilayer coatings shows a high density misfit dislocation between the TiN and TiAlN layers. The cross-plane thermal conductivity of the coatings was measured by a pulsed photothermal reflectance technique. With increasing bilayer number, the multilayer coatings' thermal conductivity decreases gradually. This reduction of thermal conductivity can be ascribed to increased phonon scattering due to the disruption of columnar structure, reduced preferential orientation, decreased grain size of the coatings and present misfit dislocations at the interfaces.

Abstract: Low-temperature growth of carbon nanotubes (CNTs) has been claimed to provide a route towards chiral-selective growth, enabling a host of applications. In this contribution, we employ reactive molecular dynamics simulations to demonstrate how plasma-based deposition allows such low-temperature growth. We first show how ion bombardment during the growth affects the carbon dissolution and precipitation process. We then continue to demonstrate how a narrow ion energy window allows CNT growth at 500 K. Finally, we also show how CNTs in contrast cannot be grown in thermal CVD at this low temperature, but only at high temperature, in agreement with experimental data. (C) 2013 Elsevier Ltd. All rights reserved.

Abstract: In this study, we investigated the diffusion of H-atoms to the subsurface and their further diffusion into the bulk of a Ni(111) crystal by means of density functional theory calculations in the context of thermal and plasma-assisted catalysis. The H-atoms at the surface can originate from the dissociative adsorption of H2 or CH4 molecules, determining the surface H-coverage. When a threshold H-coverage is passed, corresponding to 1.00 ML for the crystalline Ni(111) surface, the surface-bound H-atoms start to diffuse to the subsurface. A similar threshold coverage is observed for the interstitial H-coverage. Once the interstitial sites are filled up with a coverage above 1.00 ML of H, dissolution of interstitial H-atoms to the layer below the interstitial sites will be initiated. Hence, by applying a high pressure or inducing a reactive plasma and high temperature, increasing the H-flux to the surface, a large amount of hydrogen can diffuse in a crystalline metal like Ni and can be absorbed. The formation of metal hydride may modify the entire reaction kinetics of the system. Equivalently, the H-atoms in the bulk can easily go back to the surface and release a large amount of heat. In a plasma process, H-atoms are formed in the plasma, and therefore the energy barrier for dissociative adsorption is dismissed, thus allowing achievement of the threshold coverage without applying a high pressure as in a thermal process. As a result, depending on the crystal plane and type of metal, a large number of H-atoms can be dissolved (absorbed) in the metal catalyst, explaining the high efficiency of plasma-assisted catalytic reactions. Here, the mechanism of H-dissolution is established as a chemical identifier for the investigation of the reaction kinetics of a chemical process.

Abstract: In recent years, the potential use of hydrogen as a clean energy source has gained considerable attention. Especially H2 formation by Ni-catalyzed reforming of methane at elevated temperatures is an attractive process. However, a more fundamental knowledge at the atomic level is needed for a full comprehension of the reactions at the catalyst surface. In this contribution, we therefore investigate the H2 formation after CHx impacts on a Ni(1 1 1) surface in the temperature range 4001600 K, by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. While some H2 formation is already observed at the lower temperatures, substantial H2 formation is only obtained at elevated temperatures of 1400 K and above. At 1600 K, the H2 molecules are even the most frequently formed species. In direct correlation with the increasing dehydrogenation at elevated temperatures, an increased surface-to-subsurface C-diffusivity is observed as well. This study highlights the major importance of the temperature on the H2 formation.

Abstract: The adsorption probability and reaction behavior of CHx plasma species on various nickel catalyst surfaces is investigated by means of reactive molecular dynamics (MD) simulations using the ReaxFF potential. Such catalysts are used in the reforming of hydrocarbons and in the growth of carbon nanotubes, and further insight in the underlying mechanisms of these processes is needed to increase their applicability. Single and consecutive impacts of CHx radicals (x={1,2,3}) were performed on four different Ni surfaces, at a temperature of 400 K. The adsorption probability is shown to be related to the number of free electrons, i.e. a higher number leads to more adsorptions, and the steric hindrance caused by the hydrogen atoms bonded to the impacting CHx species. Furthermore, some of the CH bonds break after adsorption, which generally leads to diffusion of the hydrogen atom over the surface. Additionally, these adsorbed H-atoms can be used in reactions to form new molecules, such as CH4 and C2Hx, although this is dependent on the precise morphology of the surface. New molecules are also formed by subtraction of H-atoms from adsorbed radicals, leading to occasional formation of H2 and C2Hx molecules.

Abstract: A crucial property for implants is their biocompatibility. To ensure biocompatibility, thin coatings of hydroxyapatite (HA) are deposited on the actual implant. In this study, we investigate the effects of the addition of silicate anions to the structure of hydroxyapatite coatings on their adhesion strength via a scratch test and ab initio calculations. We find that both the grain size and adhesion strength decrease with the increase in the silicon content in the HA coating (SiHA). The increase in the silicon content to 1.2 % in the HA coating leads to a decrease in the average crystallite size from 28 to 21 nm, and in the case of 4.6 %, it leads to the formation of an amorphous or nanocrystalline film. The decreases in the grain and crystallite sizes lead to peeling and destruction of the coating from the titanium substrate at lower loads. Further, our ab initio simulations demonstrate an increased number of molecular bonds at the amorphous SiHA-TiO2 interface. However, the experimental results revealed that the structure and grain size have more pronounced effects on the adhesion strength of the coatings. In conclusion, based on the results of the ab initio simulations and the experimental results, we suggest that the presence of Si in the form of silicate ions in the HA coating has a significant impact on the structure, grain size, and number of molecular bonds at the interface and on the adhesion strength of the SiHA coating to the titanium substrate.

Abstract: Experiments show that the etch rate of Si with SF6 inductively coupled plasma (ICP) is significantly influenced by the absolute gas flow rate in the range of 50–600 sccm, with a maximum at around 200 sccm. Therefore, we numerically investigate the effects of the gas flow rate on the bulk plasma properties and on the etch rate, to obtain more insight in the underlying reasons of this effect. A hybrid Monte Carlo—fluid model is applied to simulate an SF6 ICP. It is found that the etch rate is influenced by two simultaneous effects: (i) the residence time of the gas and (ii) the temperature profile of the plasma in the ICP volume, resulting indeed in a maximum etch rate at 200 sccm.

Abstract: Plasma medicine has been claimed to provide a novel route to heal wounds and regenerate skin, although very little is currently known about the elementary processes taking place. We carried out a series of ReaxFF-based reactive molecular dynamics simulations to investigate the interaction of O and OH radicals with lipids, more specifically with α-linolenic acid as a model for the free fatty acids present in the upper skin layer. Our calculations predict that the O and OH radicals most typically abstract a H atom from the fatty acids, which can lead to the formation of a conjugated double bond, but also to the incorporation of alcohol or aldehyde groups, thereby increasing the hydrophilic character of the fatty acids and changing the general lipid composition of the skin. Within the limitations of the investigated model, no formation of possibly toxic products was observed.

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.

Abstract: Cold atmospheric pressure plasmas in and in contact with liquids represent a growing field of research for various applications. Understanding the interactions between the plasma generated species and the liquid is crucial. In this work we perform molecular dynamics (MD) simulations based on a quantum mechanical method, i.e. density-functional based tight-binding (DFTB), to examine the interactions of OH radicals and O atoms in bulk water. Our calculations reveal that the transport of OH radicals through water is not only governed by diffusion, but also by an equilibrium reaction of H-abstraction with water molecules. Furthermore, when two OH radicals encounter each other, they either form a stable cluster, or react, resulting in the formation of a new water molecule and an O atom. In addition,
the O atoms form either oxywater (when in singlet configuration) or they remain stable in solution (when in triplet configuration), stressing the important role that O atoms can play in aqueous solution, and in contact with biomolecules. Our observations are in line with both experimental and ab initio results from the literature.

Abstract: Molecular dynamics (MD) simulations were performed to provide atomic scale insight in the initial interaction between hydroxyl radicals (OH) and peptide systems in solution. These OH radicals are representative reactive oxygen species produced by cold atmospheric plasmas. The use of plasma for biomedical applications is gaining increasing interest, but the fundamental mechanisms behind the plasma modifications still remain largely elusive. This study helps to gain more insight in the underlying mechanisms of plasma medicine but is also more generally applicable to peptide oxidation, of interest for other applications. Combining both reactive and nonreactive MD simulations, we are able to elucidate the reactivity of the amino acids inside the peptide systems and their effect on their structure up to 1 μs. Additionally, experiments were performed, treating the simulated peptides with a plasma jet. The computational results presented here correlate well with the obtained experimental data and highlight the importance of the chemical environment for the reactivity of the individual amino acids, so that specific amino acids are attacked in higher numbers than expected. Furthermore, the long time scale simulations suggest that a single oxidation has an effect on the 3D conformation due to an increase in hydrophilicity and intra- and intermolecular interactions.