No Arabic abstract
An interatomic potential for Al-Tb alloy around the composition of Al90Tb10 was developed using the deep neural network (DNN) learning method. The atomic configurations and the corresponding total potential energies and forces on each atom obtained from ab initio molecular dynamics (AIMD) simulations are collected to train a DNN model to construct the interatomic potential for Al-Tb alloy. We show the obtained DNN model can well reproduce the energies and forces calculated by AIMD. Molecular dynamics (MD) simulations using the DNN interatomic potential also accurately describe the structural properties of Al90Tb10 liquid, such as the partial pair correlation functions (PPCFs) and the bond angle distributions, in comparison with the results from AIMD. Furthermore, the developed DNN interatomic potential predicts the formation energies of crystalline phases of Al-Tb system with the accuracy comparable to ab initio calculations. The structure factor of Al90Tb10 metallic glass obtained by MD simulation using the developed DNN interatomic potential is also in good agreement with the experimental X-ray diffraction data.
GeTe is a prototypical phase change material of high interest for applications in optical and electronic non-volatile memories. We present an interatomic potential for the bulk phases of GeTe, which is created using a neural network (NN) representation of the potential-energy surface obtained from reference calculations based on density functional theory. It is demonstrated that the NN potential provides a close to ab initio quality description of a number of properties of liquid, crystalline and amorphous GeTe. The availability of a reliable classical potential allows addressing a number of issues of interest for the technological applications of phase change materials, which are presently beyond the capability of first principles molecular dynamics simulations.
Potentials that could accurately describe the irradiation damage processes are highly desired to figure out the atomic-level response of various newly-discovered materials under irradiation environments. In this work, we introduce a deep-learning interatomic potential for monolayer MoS2 by combining all-electron calculations, an active-learning sampling method and a hybrid deep-learning model. This potential could not only give an overall good performance on the predictions of near-equilibrium material properties including lattice constants, elastic coefficients, energy stress curves, phonon spectra, defect formation energy and displacement threshold, but also reproduce the ab initial irradiation damage processes with high quality. Further irradiation simulations indicate that one single highenergy ion could generate a large nanopore with a diameter of more than 2 nm, or a series of multiple nanopores, which is qualitatively verified by the subsequent 500 keV Au+ ion irradiation experiments. This work provides a promising and feasible approach to simulate irradiation effects in enormous newly-discovered materials with unprecedented accuracy.
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 compute the thermal conductivity of water within linear response theory from equilibrium molecular dynamics simulations, by adopting two different approaches. In one, the potential energy surface (PES) is derived on the fly from the electronic ground state of density functional theory (DFT) and the corresponding analytical expression is used for the energy flux. In the other, the PES is represented by a deep neural network (DNN) trained on DFT data, whereby the PES has an explicit local decomposition and the energy flux takes a particularly simple expression. By virtue of a gauge invariance principle, established by Marcolongo, Umari, and Baroni, the two approaches should be equivalent if the PES were reproduced accurately by the DNN model. We test this hypothesis by calculating the thermal conductivity, at the GGA (PBE) level of theory, using the direct formulation and its DNN proxy, finding that both approaches yield the same conductivity, in excess of the experimental value by approximately 60%. Besides being numerically much more efficient than its direct DFT counterpart, the DNN scheme has the advantage of being easily applicable to more sophisticated DFT approximations, such as meta-GGA and hybrid functionals, for which it would be hard to derive analytically the expression of the energy flux. We find in this way, that a DNN model, trained on meta-GGA (SCAN) data, reduce the deviation from experiment of the predicted thermal conductivity by about 50%, leaving the question open as to whether the residual error is due to deficiencies of the functional, to a neglect of nuclear quantum effects in the atomic dynamics, or, likely, to a combination of the two.
Accuracy of molecular dynamics simulations depends crucially on the interatomic potential used to generate forces. The gold standard would be first-principles quantum mechanics (QM) calculations, but these become prohibitively expensive at large simulation scales. Machine learning (ML) based potentials aim for faithful emulation of QM at drastically reduced computational cost. The accuracy and robustness of an ML potential is primarily limited by the quality and diversity of the training dataset. Using the principles of active learning (AL), we present a highly automated approach to dataset construction. The strategy is to use the ML potential under development to sample new atomic configurations and, whenever a configuration is reached for which the ML uncertainty is sufficiently large, collect new QM data. Here, we seek to push the limits of automation, removing as much expert knowledge from the AL process as possible. All sampling is performed using MD simulations starting from an initially disordered configuration, and undergoing non-equilibrium dynamics as driven by time-varying applied temperatures. We demonstrate this approach by building an ML potential for aluminum (ANI-Al). After many AL iterations, ANI-Al teaches itself to predict properties like the radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. To demonstrate transferability, we perform a 1.3M atom shock simulation, and show that ANI-Al predictions agree very well with DFT calculations on local atomic environments sampled from the nonequilibrium dynamics. Interestingly, the configurations appearing in shock appear to have been well sampled in the AL training dataset, in a way that we illustrate visually.