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: Structural, mechanical and electronic properties of two-dimensional single-layer hexagonal structures in the (111) crystal plane of IIIAs-ZnS systems (III = B, Ga and In) are studied by first-principles calculations based on density functional theory (DFT). Elastic and phonon dispersion relation display that 2D h-IIIAs systems (III = B, Ga and In) are both mechanical and dynamically stable. Electronic structures analysis show that the semiconducting nature of the 3D-IIIAs compounds is retained by their 2D single layer counterpart. Furthermore, density of states reveals the influence of sigma and pi bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Calculations of elastic constants show that the Young's modulus, bulk modulus and shear modulus decrease for 2D h-IIIAs binary compounds as we move down on the group of elements of the periodic table. In addition, as the bond length between the neighboring cation-anion atoms increases, the 2D h-IIIAs binary compounds display less stiffness and more plasticity. Our findings can be used to understand the contribution of the r and p bonding in the most stable geometry (planar or buckled) for 2D h-IIIAs systems. Structural and electronic properties of h-IIIAs systems as a function of the number of layers have been also studied. It is shown that h-BAs keeps its planar geometry while both h-GAs and h-InAs retained their buckled ones obtained by their single layers. Bilayer h-IIIAs present the same bandgap nature of their counterpart in 3D. As the number of layers increase from 2 to 4, the bandgap width for layered h-IIIAs decreases until they become semimetal or metal. Interestingly, these results are different to those found for layered h-GaN. The results presented in this study for single and few-layer h-IIIAs structures could give some physical insights for further theoretical and experimental studies of 2D h-IIIV-like systems. (C) 2017 Elsevier B.V. All rights reserved.

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).