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: A Langmuir probe and an ICCD are employed to study the discharge mode transition in Ar inductively coupled plasma. Electron density and plasma emission intensity are measured during the E (capacitive discharge) to H (inductive discharge) mode transitions at different pressures. It is found that plasma exists with a low electron density and a weak emission intensity in the E mode, while it has a high electron density and a strong emission intensity in the H mode. Meanwhile, the plasma emission intensity spatial (2D image) profile is symmetrical in the H mode, but the 2D image is an asymmetric profile in the E mode. Moreover, the electron density and emission intensity jump up discontinuously at high pressure, but increase almost continuously at the E to H mode transition under low pressure.

Abstract: The requirements for new technologies to serve as anticancer agents go far beyond their toxicity potential. Novel applications also need to be safe on a molecular and patient level. In a broader sense, this also relates to cancer metastasis and inflammation. In a previous study, the toxicity of an atmospheric pressure argon plasma jet in four human pancreatic cancer cell lines was confirmed and plasma treatment did not promote metastasis in vitro and in ovo. Here, these results are extended by additional types of analysis and new models to validate and define on a molecular level the changes related to metastatic processes in pancreatic cancer cells following plasma treatment in vitro and in ovo. In solid tumors that were grown on the chorion-allantois membrane of fertilized chicken eggs (TUM-CAM), plasma treatment induced modest to profound apoptosis in the tissues. This, however, was not associated with a change in the expression levels of adhesion molecules, as shown using immunofluorescence of ultrathin tissue sections. Culturing of the cells detached from these solid tumors for 6d revealed a similar or smaller total growth area and expression of ZEB1, a transcription factor associated with cancer metastasis, in the plasma-treated pancreatic cancer tissues. Analysis of in vitro and in ovo supernatants of 13 different cytokines and chemokines revealed cell line-specific effects of the plasma treatment but a noticeable increase of, e.g., growth-promoting interleukin 10 was not observed. Moreover, markers of epithelial-to-mesenchymal transition (EMT), a metastasis-promoting cellular program, were investigated. Plasma-treated pancreatic cancer cells did not present an EMT-profile. Finally, a realistic 3D tumor spheroid co-culture model with pancreatic stellate cells was employed, and the invasive properties in a gel-like cellular matrix were investigated. Tumor outgrowth and spread was similar or decreased in the plasma conditions. Altogether, these results provide valuable insights into the effect of plasma treatment on metastasis-related properties of cancer cells and did not suggest EMT-promoting effects of this novel cancer therapy.

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: We investigated the mode transition from volume to surface discharge in a packed bed dielectric barrier discharge reactor by a two-dimensional particle-in-cell/Monte Carlo collision method. The calculations are performed at atmospheric pressure for various driving voltages and for gas mixtures with different N2 and O2 compositions. Our results reveal that both a change of the driving voltage and gas mixture can induce mode transition. Upon increasing voltage, a mode transition from hybrid (volume+surface) discharge to pure surface discharge occurs, because the charged species can escape much more easily to the beads and charge the bead surface due to the strong electric field at high driving voltage. This significant surface charging will further enhance the tangential component of the electric field along the dielectric bead surface, yielding surface ionization waves (SIWs). The SIWs will give rise to a high concentration of reactive species on the surface, and thus possibly enhance the surface activity of the beads, which might be of interest for plasma catalysis. Indeed, electron impact excitation and ionization mainly take place near the bead surface. In addition, the propagation speed of SIWs becomes faster with increasing N2 content in the gas mixture, and slower with increasing O2 content, due to the loss of electrons by attachment to O2
molecules. Indeed, the negative O-2 ion density produced by electron impact attachment is much higher than the electron and positive O+2 ion density. The different ionization rates between N2 and O2 gases will create different amounts of electrons and ions on the dielectric bead surface, which might also have effects in plasma catalysis.

Abstract: Using a molecular dynamics model the crystallinity of MgxAlyOz thin films with a variation in the stoichiometry of the thin film is studied at operating conditions similar to the experimental operating conditions of a dual magnetron sputter deposition system. The films are deposited on a crystalline or amorphous substrate. The Mg metal content in the film ranged from 100% (i.e. MgO film) to 0% (i.e. Al2O3 film). The radial distribution function and density of the films are calculated. The results are compared with x-ray diffraction and transmission electron microscopy analyses of experimentally deposited thin films by the dual magnetron reactive sputtering process. Both simulation and experimental results show that the structure of the MgAlO film varies from crystalline to amorphous when the Mg concentration decreases. It seems that the crystalline MgAlO films have a MgO structure with Al atoms in between.

Abstract: It is known that film structure may change dramatically with the extent of surface diffusion during the film growth process. In the present work, surface diffusion, induced thermally or activated by energetic impacts, is investigated theoretically under conditions appropriate for magnetron sputter-deposition of MgAlO thin films with varying stoichiometry. The distribution of surface diffusion energy barriers available to the system was determined for each stoichiometry, which allowed assessing in a qualitative way how much surface diffusion will take place on the time scale available between deposition events. The activation energy barriers increase with the Al concentration in the film, and therefore, the surface diffusion rates in the time frame of typical deposition rates drop, which can explain the decrease in crystallinity in the film structure and the transition to amorphous structure. The deposition process and the immediate surface diffusion enhanced by the energetic adatoms are simulated by means of a molecular dynamics model. The longer-time thermal surface diffusion and the energy landscape are studied by the temperature accelerated dynamics method, applied in an approximate way. The surface diffusion enhanced by the energetic impacts appears to be very important for the film structure in the low-temperature deposition regime.

Abstract: Molecular dynamics simulations are essential for a better understanding of dissociative chemisorption on metal surfaces, which is often the rate-controlling step in heterogeneous and plasma catalysis. The workhorse quasi-classical trajectory approach ubiquitous in molecular dynamics is able to accurately predict reactivity only for high translational and low vibrational energies. In contrast, catalytically relevant conditions generally involve low translational and elevated vibrational energies. Existing quantum dynamics approaches are intractable or approximate as a result of the large number of degrees of freedom present in molecule−metal surface reactions. Here, we extend a ring polymer molecular dynamics approach to fully include, for the first time, the degrees of freedom of a moving metal surface. With this approach, experimental sticking probabilities for the dissociative chemisorption of methane on Pt(111) are reproduced for a large range of translational and vibrational energies by including nuclear quantum effects and employing full-dimensional simulations.