Modified embedded-atom method interatomic potentials for the Mg-Al alloy system

We developed new modified embedded-atom method (MEAM) interatomic potentials for the Mg-Al alloy system using a first-principles method based on density functional theory (DFT). The materials parameters, such as the cohesive energy, equilibrium atomic volume, and bulk modulus, were used to determine the MEAM parameters. Face-centered cubic, hexagonal close packed, and cubic rock salt structures were used as the reference structures for Al, Mg, and MgAl, respectively. The applicability of the new MEAM potentials to atomistic simulations for investigating Mg-Al alloys was demonstrated by performing simulations on Mg and Al atoms in a variety of geometries. The new MEAM potentials were used to calculate the adsorption energies of Al and Mg atoms on Al (111) and Mg (0001) surfaces. The formation energies and geometries of various point defects, such as vacancies, interstitial defects and substitutional defects, were also calculated. We found that the new MEAM potentials give a better overall agreement with DFT calculations and experiments when compared against the previously published MEAM potentials.Comment: Fixed a referenc

An interatomic potential for saturated hydrocarbons based on the modified embedded-atom method

In this work, we developed an interatomic potential for saturated hydrocarbons using the modified embedded-atom method (MEAM), a reactive semi-empirical many-body potential based on density functional theory and pair potentials. We parameterized the potential by fitting to a large experimental and first-principles (FP) database consisting of 1) bond distances, bond angles, and atomization energies at 0 K of a homologous series of alkanes and their select isomers from methane to n-octane, 2) the potential energy curves of H$_2$, CH, and C$_2$ diatomics, 3) the potential energy curves of hydrogen, methane, ethane, and propane dimers, i.e., (H$_2$)$_2$, (CH$_4$)$_2$, (C$_2$H$_6$)$_2$, and (C$_3$H$_8$)$_2$, respectively, and 5) pressure-volume-temperature (PVT) data of a dense high-pressure methane system with the density of 0.5534 g/cc. We compared the atomization energies and geometries of a range of linear alkanes, cycloalkanes, and free radicals calculated from the MEAM potential to those calculated by other commonly used reactive potentials for hydrocarbons, i.e., second-generation reactive empirical bond order (REBO) and reactive force field (ReaxFF). MEAM reproduced the experimental and/or FP data with accuracy comparable to or better than REBO or ReaxFF. The experimental PVT data for a relatively large series of methane, ethane, propane, and butane systems with different densities were predicted reasonably well by MEAM. Although the MEAM formalism has been applied to atomic systems with predominantly metallic bonding in the past, the current work demonstrates the promising extension of the MEAM potential to covalently bonded molecular systems, specifically saturated hydrocarbons and saturated hydrocarbon-based polymers

Thermal conductivity of strained silicon: molecular dynamics insight and kinetic theory approach

In this work, we investigated tensile and compression forces effect on the thermal conductivity of silicon. We used equilibrium molecular dynamics approach for the evaluation of thermal conductivity considering different interatomic potentials. More specifically, we tested Stillinger-Weber, Tersoff, Environment-Dependent Interatomic Potential and Modified Embedded Atom Method potentials for the description of silicon atom motion under different strain and temperature conditions. Additionally, we extracted phonon density of states and dispersion curves from molecular dynamics simulations. These data were used for direct calculations of thermal conductivity considering the kinetic theory approach. Comparison of molecular dynamics and kinetic theory simulations results as a function of strain and temperature allowed us to investigate the different factors affecting the thermal conductivity of strained silicon

Strain-induced structural transformation of a silver nanowire

We have investigated the structural characteristics of the experimentally observed phase transition of a silver nanowire into a tube under tensile strain. In the simulations, atoms are allowed to interact via a model potential extracted from the modified embedded atom method. Our calculations demonstrate that the formation of the hollow structure is governed by the nature of the applied strain, the length of the wire, and the initial cross-sectional shape. The results further offer insights into the atomistic nature of this specific structural transformation into a nanotube with the smallest possible cross-section

Structural, elastic and thermal properties of cementite (Fe3_3C) calculated using Modified Embedded Atom Method

Structural, elastic and thermal properties of cementite (Fe$_3$C) were studied using a Modified Embedded Atom Method (MEAM) potential for iron-carbon (Fe-C) alloys. Previously developed Fe and C single element potentials were used to develop an Fe-C alloy MEAM potential, using a statistically-based optimization scheme to reproduce structural and elastic properties of cementite, the interstitial energies of C in bcc Fe as well as heat of formation of Fe-C alloys in L$_{12}$ and B$_1$ structures. The stability of cementite was investigated by molecular dynamics simulations at high temperatures. The nine single crystal elastic constants for cementite were obtained by computing total energies for strained cells. Polycrystalline elastic moduli for cementite were calculated from the single crystal elastic constants of cementite. The formation energies of (001), (010), and (100) surfaces of cementite were also calculated. The melting temperature and the variation of specific heat and volume with respect to temperature were investigated by performing a two-phase (solid/liquid) molecular dynamics simulation of cementite. The predictions of the potential are in good agreement with first-principles calculations and experiments.Comment: 12 pages, 9 figure

Non-clasical Nucleation in Supercooled Nickel

The dynamics of homogeneous nucleation and growth of crystalline nickel from the super-cooled melt is examined during rapid quenching using molecular dynamics and a modified embedded atom method potential. The character of the critical nuclei of the crystallization transition is examined using common neighbor analysis and visualization. At nucleation the saddle point droplet consists of randomly stacked planar structures with an in plane triangular order. These results are consistent with previous theoretical results that predict that the nucleation process in some metals is non-classical due to the presence of long-range forces and a spinodal.Comment: 4 pages, 5 figure

Understanding Homogeneous Nucleation in Solidification of Aluminum by Molecular Dynamics Simulations

Homogeneous nucleation from aluminum (Al) melt was investigated by million-atom molecular dynamics (MD) simulations utilizing the second nearest neighbor modified embedded atom method (MEAM) potentials. The natural spontaneous homogenous nucleation from the Al melt was produced without any influence of pressure, free surface effects and impurities. Initially isothermal crystal nucleation from undercooled melt was studied at different constant temperatures, and later superheated Al melt was quenched with different cooling rates. The crystal structure of nuclei, critical nucleus size, critical temperature for homogenous nucleation, induction time, and nucleation rate were determined. The quenching simulations clearly revealed three temperature regimes: sub-critical nucleation, super-critical nucleation, and solid-state grain growth regimes. The main crystalline phase was identified as face-centered cubic (fcc), but a hexagonal close-packed (hcp) and an amorphous solid phase were also detected. The hcp phase was created due to the formation of stacking faults during solidification of Al melt. By slowing down the cooling rate, the volume fraction of hcp and amorphous phases decreased. After the box was completely solid, grain growth was simulated and the grain growth exponent was determined for different annealing temperatures.Comment: 41 page