Abstract: Nonthermal atmospheric pressure plasmas are gaining increasing attention for biomedical applications. However, very little fundamental information on the interaction mechanisms between the plasma species and biological cells is currently available. We investigate the interaction of important plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, with bacterial peptidoglycan by means of reactive molecular dynamics simulations, aiming for a better understanding of plasma disinfection. Our results show that OH, O, O3, and H2O2 can break structurally important bonds of peptidoglycan (i.e., CO, CN, or CC bonds), which consequently leads to the destruction of the bacterial cell wall. The mechanisms behind these breakups are, however, dependent on the impinging plasma species, and this also determines the effectiveness of the cell wall destruction.

Abstract: The application of nonthermal atmospheric pressure plasma is emerging as an alternative and efficient technique for the inactivation of bacterial biofilms. In this study, reactive molecular dynamics simulations were used to examine the reaction mechanisms of hydroxyl radicals, as key reactive oxygen plasma species in biological systems, with several organic molecules (i.e., alkane, alcohol, carboxylic acid, and amine), as prototypical components of biomolecules in the biofilm. Our results demonstrate that organic molecules containing hydroxyl and carboxyl groups may act as trapping agents for the OH radicals. Moreover, the impact of OH radicals on N-acetyl-glucosamine, as constituent component of staphylococcus epidermidis biofilms, was investigated. The results show how impacts of OH radicals lead to hydrogen abstraction and subsequent molecular damage. This study thus provides new data on the reaction mechanisms of plasma species, and particularly the OH radicals, with fundamental components of bacterial biofilms.

Abstract: The possibility of application of nanomaterials is determined by our ability to control the properties of the materials, which are ultimately determined by their structure and hence their growth processes. We employ hybrid molecular dynamics / Monte Carlo (MD/MC) simulations to explore the growth of SWNTs and graphene on nickel as a catalyst, with the specific goal of unraveling the growth mechanisms. While the general observations are in agreement with the literature, we find a number of interesting phenomena to be operative which are crucial for the growth, and which are not accessible by MD simulations alone due to the associated time scale. Specifically, we observe metal mediated healing and restructuring processes to take place, reorganizing the carbon network during the initial nucleation step. In the case of carbon nanotube growth, this leads to the growth of tubes with a determinable chirality. In the case of graphene formation, we find that graphene is only formed at temperatures above 700 K. These results are of importance for understanding the growth mechanisms of these carbon nanomaterials on the fundamental level.

Abstract: Nanocrystalline diamond films have generated much interested due to their diamond-like properties and low surface roughness. Several techniques have been used to obtain a high re-nucleation rate, such as hydrogen poor or high methane concentration plasmas. In this work, the properties of nano-diamond films grown on silicon substrates using a continuous DC bias voltage during the complete duration of growth are studied. Subsequently, the layers were characterised by several morphological, structural and optical techniques. Besides a thorough investigation of the surface structure, using SEM and AFM, special attention was paid to the bulk structure of the films. The application of FTIR, XRD, multi wavelength Raman spectroscopy, TEM and EELS yielded a detailed insight in important properties such as the amount of crystallinity, the hydrogen content and grain size. Although these films are smooth, they are under a considerable compressive stress. FTIR spectroscopy points to a high hydrogen content in the films, while Raman and EELS indicate a high concentration of sp2 carbon. TEM and EELS show that these films consist of diamond nano-grains mixed with an amorphous sp2 bonded carbon, these results are consistent with the XRD and UV Raman spectroscopy data.

Abstract: The recent advantages of the fabrication of artificial nanochannels enabled new research on the molecular transport, permeance, and selectivity of various gases and molecules. However, the physisorption/chemisorption of the unwanted molecules (usually hydrocarbons) inside nanochannels results in the alteration of the functionality of the nanochannels. We investigated contamination due to hydrocarbon molecules, nanochannels made of graphene, hexagonal boron nitride, BC2N, and molybdenum disulfide using molecular dynamics simulations. We found that for a certain size of nanochannel (i.e., h = 0.7 nm), as a result of the anomalous hydrophilic nature of nanochannels made of graphene, the hydrocarbons are fully adsorbed in the nanochannel, giving rise to full uptake. An increasing temperature plays an important role in unclogging, while pressure does not have a significant role. The results of our pioneering work contribute to a better understanding and highlight the important factors in alleviating the contamination and unclogging of nanochannels, which are in good agreement with the results of recent experiments.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.

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: The application of coreshell SiSiO2 nanowires as nanoelectronic devices strongly depends on their structure, which is difficult to tune precisely. In this work, we investigate the formation of the coreshell nanowires at the atomic scale, by reactive molecular dynamics simulations. The occurrence of two temperature-dependent oxidation mechanisms of ultra-small diameter Si-NWs is demonstrated. We found that control over the Si-core radius and the SiOx (x ≤ 2) oxide shell is possible by tuning the growth temperature and the initial Si-NW diameter. Two different structures were obtained, i.e., ultrathin SiO2 silica nanowires at high temperature and Si core|ultrathin SiO2 silica nanowires at low temperature. The transition temperature is found to linearly decrease with the nanowire curvature. Finally, the interfacial stress is found to be responsible for self-limiting oxidation, depending on both the initial Si-NW radius and the oxide growth temperature. These novel insights allow us to gain control over the exact morphology and structure of the wires, as is needed for their application in nanoelectronics.

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: We performed molecular dynamics simulations to investigate the effect of lipid peroxidation products on the structural and dynamic properties of the cell membrane. Our simulations predict that the lipid order in a phospholipid bilayer, as a model system for the cell membrane, decreases upon addition of lipid peroxidation products. Eventually, when all phospholipids are oxidized, pore formation can occur. This will allow reactive species, such as reactive oxygen and nitrogen species (RONS), to enter the cell and cause oxidative damage to intracellular macromolecules, such as DNA or proteins. On the other hand, upon increasing the cholesterol fraction of lipid bilayers, the cell membrane order increases, eventually reaching a certain threshold, from which cholesterol is able to protect the membrane against pore formation. This finding is crucial for cancer treatment by plasma technology, producing a large number of RONS, as well as for other cancer treatment methods that cause an increase in the concentration of extracellular RONS. Indeed, cancer cells contain less cholesterol than their healthy counterparts. Thus, they will be more vulnerable to the consequences of lipid peroxidation, eventually enabling the penetration of RONS into the interior of the cell, giving rise to oxidative stress, inducing pro-apoptotic factors. This provides, for the first time, molecular level insight why plasma can selectively treat cancer cells, while leaving their healthy counterparts undamaged, as is indeed experimentally demonstrated.

Abstract: Hydrogen–graphite interactions are relevant to a wide variety of applications, ranging from astrophysics to fusion devices and nano-electronics. In order to shed light on these interactions, atomistic simulation using Molecular Dynamics (MD) has been shown to be an invaluable tool. It suffers, however, from severe timescale
limitations. In this work we apply the recently developed Collective Variable-Driven Hyperdynamics (CVHD) method to hydrogen etching of graphite for varying inter-impact times up to a realistic value of 1 ms, which corresponds to a flux of
1020 m2 s1. The results show that the erosion yield, hydrogen surface coverage and species distribution are significantly affected by the time between impacts. This can be explained by the higher probability of C–C bond breaking due to the prolonged exposure to thermal stress and the subsequent transition from ion- to thermal-induced etching. This latter regime of thermal-induced etching – chemical erosion – is here accessed for the first time using atomistic simulations. In conclusion, this study demonstrates that accounting for long time-scales significantly affects ion bombardment simulations and should not be neglected in a wide range of conditions, in contrast to what is typically assumed.

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: Bacteria can be inactivated through various physical and chemical means, and these have always been the focus of extensive research. To further improve the methodology for these ends, two types of plasma systems were investigated: nano-second pulsed plasma (NPP) as liquid discharge plasma and an Argon gas-feeding dielectric barrier discharge (Ar-DBD) as a form of surface plasma. To understand the sterilizing action of these two different plasma sources, we performed experiments with Staphylococcus aureus (S. aureus) bacteria (wild type) and multidrug resistant bacteria (Penicillum-resistant, Methicillin-resistant and Gentamicin-resistant). We observed that both plasma sources can inactivate both the wild type and multidrug-resistant bacteria to a good extent. Moreover, we observed a change in the surface morphology, gene expression and β-lactamase activity. Furthermore, we used X-ray photoelectron spectroscopy to investigate the variation in functional groups (C-H/C-C, C-OH and C=O) of the peptidoglycan (PG) resulting from exposure to plasma species. To obtain atomic scale insight in the plasma-cell interactions and support our experimental observations, we have performed molecular dynamics simulations to study the effects of plasma species, such as OH, H2O2, O, O3, as well as O2 and H2O, on the dissociation/formation of above mentioned functional groups in PG.

Abstract: The deposition of ultra-thin tungsten films and their related 2D chalcogen compounds on large area dielectric substrates by gas phase reactions is challenging. The lack of nucleation sites complicates the adsorption of W-related precursors and subsequent sulfurization usually requires high temperatures. We propose here a technique in which a thin solid amorphous silicon film is used as reductant for the gas phase precursor WF6 leading to the conversion to metallic W. The selectivity of the W conversion towards the underlying dielectric surfaces is demonstrated. The role of the Si surface preparation, the conversion temperature, and Si thickness on the formation process is investigated. Further, the in situ conversion of the metallic tungsten into thin stoichiometric WS2 is achieved by a cyclic approach based on WF6 and H2S pulses at the moderate temperature of 450 1C, which is much lower than usual oxide sulfurization processes.

Abstract: Heterostructures of low-dimensional semiconducting materials, such as transition metal dichalcogenides (MX2), are promising building blocks for future electronic and optoelectronic devices. The patterning of one MX2 material on top of another one is challenging due to their structural similarity. This prevents an intrinsic etch stop when conventional anisotropic dry etching processes are used. An alternative approach consist in a two-step process, where a sacrificial silicon layer is pre-patterned with a low damage plasma process, stopping on the underlying MoS2 film. The pre-patterned layer is used as sacrificial template for the formation of the top WS2 film. This study describes the optimization of a cyclic Ar/Cl-2 atomic layer etch process applied to etch silicon on top of MoS2, with minimal damage, followed by a selective conversion of the patterned Si into WS2. The impact of the Si atomic layer etch towards the MoS2 is evaluated: in the ion energy range used for this study, MoS2 removal occurs in the over-etch step over 1-2 layers, leading to the appearance of MoOx but without significant lattice distortions to the remaining layers. The combination of Si atomic layer etch, on top of MoS2, and subsequent Si-to-WS2 selective conversion, allows to create a WS2/MoS2 heterostructure, with clear Raman signals and horizontal lattice alignment. These results demonstrate a scalable, transfer free method to achieve horizontally individually patterned heterostacks and open the route towards wafer-level processing of 2D materials.