No Arabic abstract
Recent application of neural networks (NNs) to modeling interatomic interactions has shown the learning machines encouragingly accurate performance for select elemental and multicomponent systems. In this study, we explore the possibility of building a library of NN-based models by introducing a hierarchical NN training. In such a stratified procedure NNs for multicomponent systems are obtained by sequential training from the bottom up: first unaries, then binaries, and so on. Advantages of constructing NN sets with shared parameters include acceleration of the training process and intact description of the constituent systems. We use an automated generation of diverse structure sets for NN training on density functional theory-level reference energies. In the test case of Cu, Pd, Ag, Cu-Pd, Cu-Ag, Pd-Ag, and Cu-Pd-Ag systems, NNs trained in the traditional and stratified fashions are found to have essentially identical accuracy for defect energies, phonon dispersions, formation energies, etc. The models robustness is further illustrated via unconstrained evolutionary structure searches in which the NN is used for the local optimization of crystal unit cells.
Graph neural networks, trained on experimental or calculated data are becoming an increasingly important tool in computational materials science. Networks, once trained, are able to make highly accurate predictions at a fraction of the cost of experiments or first-principles calculations of comparable accuracy. However these networks typically rely on large databases of labelled experiments to train the model. In scenarios where data is scarce or expensive to obtain this can be prohibitive. By building a neural network that provides a confidence on the predicted properties, we are able to develop an active learning scheme that can reduce the amount of labelled data required, by identifying the areas of chemical space where the model is most uncertain. We present a scheme for coupling a graph neural network with a Gaussian process to featurise solid-state materials and predict properties textit{including} a measure of confidence in the prediction. We then demonstrate that this scheme can be used in an active learning context to speed up the training of the model, by selecting the optimal next experiment for obtaining a data label. Our active learning scheme can double the rate at which the performance of the model on a test data set improves with additional data compared to choosing the next sample at random. This type of uncertainty quantification and active learning has the potential to open up new areas of materials science, where data are scarce and expensive to obtain, to the transformative power of graph neural networks.
Module for ab initio structure evolution (MAISE) is an open-source package for materials modeling and prediction. The codes main feature is an automated generation of neural network (NN) interatomic potentials for use in global structure searches. The systematic construction of Behler-Parrinello-type NN models approximating ab initio energy and forces relies on two approaches introduced in our recent studies. An evolutionary sampling scheme for generating reference structures improves the NNs mapping of regions visited in unconstrained searches, while a stratified training approach enables the creation of standardized NN models for multiple elements. A more flexible NN architecture proposed here expands the applicability of the stratified scheme for an arbitrary number of elements. The full workflow in the NN development is managed with a customizable MAISE-NET wrapper written in Python. The global structure optimization capability in MAISE is based on an evolutionary algorithm applicable for nanoparticles, films, and bulk crystals. A multitribe extension of the algorithm allows for an efficient simultaneous optimization of nanoparticles in a given size range. Implemented structure analysis functions include fingerprinting with radial distribution functions and finding space groups with the SPGLIB tool. This work overviews MAISEs available features, constructed models, and confirmed predictions.
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.
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.
We propose a simple scheme to construct composition-dependent interatomic potentials for multicomponent systems that when superposed onto the potentials for the pure elements can reproduce not only the heat of mixing of the solid solution in the entire concentration range but also the energetics of a wider range of configurations including intermetallic phases. We show that an expansion in cluster interactions provides a way to systematically increase the accuracy of the model, and that it is straightforward to generalise this procedure to multicomponent systems. Concentration-dependent interatomic potentials can be built upon almost any type of potential for the pure elements including embedded atom method (EAM), modified EAM, bond-order, and Stillinger-Weber type potentials. In general, composition-dependent N-body terms in the total energy lead to explicit (N+1)-body forces, which potentially renders them computationally expensive. We present an algorithm that overcomes this problem and that can speed up the calculation of the forces for composition-dependent pair potentials in such a way as to make them computationally comparable in efficiency and scaling behaviour to standard EAM potentials. We also discuss the implementation in Monte-Carlo simulations. Finally, we exemplarily review the composition-dependent EAM model for the Fe-Cr system [PRL 95, 075702 (2005)].