Abstract: Hydroxyapatite (HAP) is frequently used as biocompatible coating on Ti-based implants. In this context, the HAP-Ti adhesion is of crucial importance. Here, we report ab initio calculations to investigate the influence of Si incorporation into the amorphous calcium-phosphate (a-HAP) structure on the interfacial bonding mechanism between the a-HAP coating and an amorphous titanium dioxide (a-TiO2) substrate, contrasting two different density functionals: PBE-GGA, and DFT-D3, which are capable of describing the influence of the van der Waals (vdW) interactions. In particular, we discuss the effect of dispersion on the work of adhesion (W-ad), equilibrium geometries, and charge density difference (CDD). We find that replacement of P by Si in a-HAP (a-Si-HAP) with the creation of OH vacancies as charge compensation results in a significant increase in the bond strength between the coating and substrate in the case of using the PBE-GGA functional. However, including the vdW interactions shows that these forces considerably contribute to the W-ad. We show that the difference (W-ad – W-ad(vdW)) is on average more than 1.1 J m(-2) and 0.5 J m(-2) for a-HAP/a-TiO2 and a-Si-HAP/a-TiO2, respectively. These results reveal that including vdW interactions is essential for accurately describing the chemical bonding at the a-HAP/a-TiO2 interface.

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: The adsorption of C and CHx radicals on anatase (001) was studied using DFT within the generalized gradient approximation using the Perde-Burke-Ernzerhof (PBE) functional. We have studied the influence of oxygen vacancies in and at the surface on the adsorption properties of the radicals. For the oxygen vacancies in anatase (001), the most stable vacancy is located at the surface. For this vacancy, the maximal adsorption strength of C and CH decreases compared to the adsorption on the stoichiometric surface, but it increases for CH2 and CH3. If an oxygen vacancy is present in the first subsurface layer, the maximal adsorption strength increases for C, CH, CH2, and CH3. When the vacancy is present in the next subsurface layer, we find that only the CH3 adsorption is enhanced, while the maximal adsorption energies for the other radical species decrease. Not only does the precise location of the oxygen vacancy determine the maximal adsorption interaction, it also influences the adsorption strengths of the radicals at different surface configurations. This determines the probability of finding a certain adsorption configuration at the surface, which in turn influences the possible surface reactions. We find that C preferentially adsorbs far away from the oxygen vacancy, while CH2 and CH3 adsorb preferentially at the oxygen vacancy site. A fraction of CH partially adsorbs at the oxygen vacancy, and another fraction adsorbs further away from the vacancy.

Abstract: We discuss the synthesis of carbon nanotubes (CNTs) and graphene by catalytic chemical vapour deposition (CCVD) and plasma-enhanced CVD (PECVD), summarising the state-of-the-art understanding of mechanisms controlling their growth rate, chiral angle, number of layers (walls), diameter, length and quality (defects), before presenting a new model for 2D nucleation of a graphene sheet from amorphous carbon on a nickel surface. Although many groups have modelled this process using a variety of techniques, we ask whether there are any complementary ideas emerging from the different proposed growth mechanisms, and whether different modelling techniques can give the same answers for a given mechanism. Subsequently, by comparing the results of tight-binding, semi-empirical molecular orbital theory and reactive bond order force field calculations, we demonstrate that graphene on crystalline Ni(111) is thermodynamically stable with respect to the corresponding amorphous metal and carbon structures. Finally, we show in principle how a complementary heterogeneous nucleation step may play a key role in the transformation from amorphous carbon to graphene on the metal surface. We conclude that achieving the conditions under which this complementary crystallisation process can occur may be a promising method to gain better control over the growth processes of both graphene from flat metal surfaces and CNTs from catalyst nanoparticles.

Abstract: Metal-catalyzed growth mechanisms of carbon nanotubes (CNTs) were studied by hybrid molecular dynamics−Monte Carlo simulations using a recently developed ReaxFF reactive force field. Using this novel approach, including relaxation effects, a CNT with definable chirality is obtained, and a step-by-step atomistic description of the nucleation process is presented. Both root and tip growth mechanisms are observed. The importance of the relaxation of the network is highlighted by the observed healing of defects.

Abstract: In this contribution, we present an overview of the various techniques for combining atomistic molecular dynamics with Monte Carlo simulations, mainly in the context of condensed matter systems, as well as a brief summary of the main accelerated dynamics techniques. Special attention is given to the force bias Monte Carlo technique and its combination with molecular dynamics, in view of promising recent developments, including a definable timescale. Various examples of the application of combined molecular dynamics / Monte Carlo simulations are given, in order to demonstrate the enhanced simulation efficiency with respect to either pure molecular dynamics or Monte Carlo.

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: 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: 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: The interaction of thermal and hyperthermal Ni ions with gas-phase C60 fullerene was investigated at two temperatures with classical molecular dynamics simulations using a recently developed interatomic many-body potential. The interaction between Ni and C60 is characterized in terms of the NiC60 binding sites, complex formation, and the collision and temperature induced deformation of the C60 cage structure. The simulations show how ion implantation theoretically allows the synthesis of both endohedral Ni@C60 and exohedral NiC60 metallofullerene complexes.

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: Using reactive molecular dynamics simulations, we have investigated the effect of single-impact, low-energy (thermal-100 eV) bombardment of a Si(100){2 × 1} surface by atomic and molecular oxygen. Penetration probability distributions, as well as defect formation distributions, are presented as a function of the impact energy for both species. It is found that at low impact energy, defects are created chemically due to the chemisorption process in the top layers of the surface, while at high impact energy, additional defects are created by a knock-on displacement of Si. These results are of particular importance for understanding device performances of silica-based microelectronic and photovoltaic devices.

Abstract: This paper reports on the incorporation of three commercial fluorescent dyes, i.e., rhodamine 6G, fluorescein, and fluorescent brightener 184, in plasma coatings, by utilizing a dielectric barrier discharge (DBD) reactor, and the subsequent monitoring of the coatings homogeneity based on the emitted fluorescent light. The plasma coatings are qualitatively characterized with fluorescence microscopy, UVvis spectroscopy and profilometry for the determination of the coating thickness. The emitted fluorescent light of the coating correlates to the amount of dye per area, and deviations of these factors can hence be observed by monitoring the intensity of this light. This allows monitoring the homogeneity of the plasma coatings in a fast and simple way, without making major adjustments to the process.

Abstract: In this paper, we quantitatively investigate with atom probe tomography, the effect of temperature on the interfacial transition layer suboxide species due to the thermal oxidation of silicon. The chemistry at the interface was measured with atomic scale resolution, and the changes in chemistry and intermixing at the interface were identified on a nanometer scale. We find an increase of suboxide (SiOx) concentration relative to SiO2 and increased oxygen ingress with elevated temperatures. Our experimental findings are in agreement with reactive force field molecular dynamics simulations. This work demonstrates the direct comparison between atom probe derived chemical profiles and atomistic-scale simulations for transitional interfacial layer of suboxides as a function of temperature.