No Arabic abstract
We develop an Fe-C-H interatomic potential based on the modified embedded-atom method (MEAM) formalism based on density functional theory to enable large-scale modular dynamics simulations of carbon steel and hydrogen.
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.
We have developed a multi-objective optimization (MOO) procedure to construct modified-embedded-atom-method (MEAM) potentials with minimal manual fitting. This procedure has been applied successfully to develop a new MEAM potential for magnesium. The MOO procedure is designed to optimally reproduce multiple target values that consist of important materials properties obtained from experiments and first-principles calculations based on density-functional theory (DFT). The optimized target quantities include elastic constants, cohesive energies, surface energies, vacancy formation energies, and the forces on atoms in a variety of structures. The accuracy of the new potential is assessed by computing several material properties of Mg and comparing them with those obtained from other potentials previously published. We found that the present MEAM potential yields a significantly better overall agreement with DFT calculations and experiments.
We propose a modification of the embedded-atom method-type potential aiming at reconciling simulated melting and ground-state properties of metals by means of classical molecular dynamics. Considering titanium, magnesium, gold, and platinum as case studies, we demonstrate that simulations performed with the modified force field yield quantitatively correctly both the melting temperature of the metals and their ground-state properties. It is shown that the accounting for the long-range interatomic interactions noticeably affect the melting point assessment. The introduced modification weakens the interaction at interatomic distances exceeding the equilibrium one by a characteristic vibration amplitude defined by the Lindemann criterion, thus allowing for the correct simulation of melting, while keeping its behavior in the vicinity of the ground state minimum. The modification of the many-body potential has a general nature and can be applicable to metals with different characteristics of the electron structure as well as for many different molecular and solid state systems experiencing phase transitions.
Molecular dynamics (MD) simulation of dislocation migration requires semi-empirical potentials of the interatomic interaction. While there are many reliable semi-empirical potentials for the bcc Fe, the number of the available potentials for the fcc is very limited. In the present study we tested three EAM potentials for the fcc Fe (ABCH97 [Phil. Mag. A, 75, 713-732 (1997)], BCT13 [MSMSE 21, 085004 (2013)] and ZFS18 [J. Comp. Chem. 39, 2420-2431 (2018)]) from literature. It was found that the ABCH97 potential does not provide that the fcc phase is the most stable at any temperature. On the other hand, the fcc phase is always more stable than the bcc phase for the BCT13, ZFS18 potentials. The hcp phase is the most stable phase for the BCT13 potential at any temperature. In order to fix these problems we developed two new EAM potentials (MB1 and MB2). The fcc phase is still more stable than the bcc phase for the MB1 potential but the MB2 potential provides that the bcc phase is the most stable phase from the upper fcc-bcc transformation temperature, T_gamma-delta, to the melting temperature, Tm, and the fcc phase is the most stable phase below T_gamma-delta. This potential also leads to an excellent agreement with the experimental data on the fcc elastic constants and reasonable stacking fault energy which makes it the best potential for the simulation of the dislocation migration in the fcc Fe among all semi-empirical potentials considered in the present study. The MD simulation demonstrated that only the ZFS18, MB1 and MB2 potentials are actually suitable for the simulation of the dislocation migration in the fcc Fe. They lead to the same orders of magnitude for the dislocation velocities and all of them show that the edge dislocation is faster than the screw dislocation. However, the actual values of the dislocation velocities do depend on the employed semi-empirical potential.
Thermal management materials are of critical importance for engineering miniaturized electronic devices, where theoretical design of such materials demands the evaluation of thermal conductivities which are numerically expensive. In this work, we applied the recently developed machine learning interatomic potential (MLIP) to evaluate the thermal conductivity of hexagonal boron nitride monolayers. The MLIP is obtained using the Gaussian approximation potential (GAP) method, and the resulting lattice dynamical properties and thermal conductivity are compared with those obtained from explicit frozen phonon calculations. It is observed that accurate thermal conductivity can be obtained based on MLIP constructed with about 30% representative configurations, and the high-order force constants provide a more reliable benchmark on the quality of MLIP than the harmonic approximation.