Exact average many-body interatomic interaction model for random alloys

Machine learning (ML) interatomic potentials (ML-IAPs) are generated for alkane and polyene hydrocarbons using on-the-fly adaptive sampling and a sparse Gaussian process regression (SGPR) algorithm. The ML model is generated based on the PBE+D3 level of density functional theory (DFT) with molecular dynamics (MD) for small alkane and polyene molecules. Intermolecular interactions are also trained with clusters and condensed phases of small molecules. It shows excellent transferability to long alkanes and closely describes the ab inito potential energy surface for polyenes. Simulation of liquid ethane also shows reasonable agreement with experimental reports. This is a promising initiative toward a universal ab initio quality force-field for hydrocarbons and organic molecules.

Two-dimensional (2D) quantum materials are poised to transform conventional electronics for a wide spectrum of applications that will encompass chemical sciences. For the study of thermal transport in single-layer (1L) or multi-layer transition metal dichalcogenides (TMDs), this work explores the combination of density functional theory (DFT) and algorithmic training for the generation of a moment tensor potential (MTP) that models 1L-MoS2, 1L-WS2 and their alloys, and demonstrates a synergy of theoretical techniques that is anticipated to play an important role in the field. From a high-performance computing perspective, these yield very convenient inter-atomic (or inter-molecular in other contexts) potentials that are useful to predict the response of quantum materials to thermal perturbations, or other driving forces. We show that our trained MTP functions successfully describe vibrational properties of the systems, and their thermal conductivities. The trained potential displays consistent agreement with DFT calculations, as well as the Stillinger-Weber (SW) potential. We also find that the thermal conductivity of the 2D alloys is little affected by sulfur vacancies. This is a behavior that may aid the fine-tuning of material's thermal properties for heat management and energy storage and conversion applications.

Exploring phononic properties of two-dimensional materials using machine learning interatomic potentials

Dopants can tune the performance of MoS2 in various applications, but the use of molecular dynamics simulations for doped MoS2 materials discovery is limited by the lack of multidopant interatomic potentials. Universal machine learning interatomic potentials (MLIPs) could be a solution, but the accuracy of these potentials must first be evaluated. Here, we evaluate the accuracy of a recently developed MLIP, META's Universal Model for Atoms, for 25 different MoS2 dopants spanning metals, nonmetals, and transition metals in Mo-substitution, S-substitution, and intercalated positions by benchmarking the MLIP-predicted formation energy and the dopant-induced structural change against density functional theory (DFT) calculations. The computational framework for MLIP validation and simulations is described in detail, and the source code is made available online. The MLIP is then demonstrated by performing heating-cooling simulations of MoS2 supercells with all 25 dopants. These simulations capture complex phenomena including dopant clustering, MoS2 layer fracturing, interlayer diffusion, and chemical compound formation at orders-of-magnitude reduced computational cost compared to DFT. This work provides a computational workflow for the application-oriented design of doped-MoS2, enabling high-throughput screening of dopant candidates and optimization of compositions for targeted tribological, electronic, and optoelectronic performance.

Chemical heuristics have been fundamental to the advancement of chemistry and materials science. These heuristics are typically established by scientists using knowledge and creativity to extract patterns from limited datasets. Machine learning offers opportunities to perfect this approach using computers and larger datasets. Here, we discuss the relationships between traditional heuristics and machine learning approaches. We show how traditional rules can be challenged by large-scale statistical assessment and how traditional concepts commonly used as features are feeding the machine learning techniques. We stress the waste involved in relearning chemical rules and the challenges in terms of data size requirements for purely data-driven approaches. Our view is that heuristic and machine learning approaches are at their best when they work together.

Lithium-based disordered rocksalts (LDRs), which are an important class of positive electrode materials that can increase the energy density of current Li-ion batteries, represent a significantly complex chemical and configurational space for conventional density functional theory (DFT)-based high-throughput screening approaches. Notably, atom-centered machine-learned interatomic potentials (MLIPs) are a promising pathway to accurately model the potential energy surface of highly disordered chemical spaces, such as LDRs, where the performance of such MLIPs has not been rigorously explored yet. Here, we represent a comprehensive evaluation of the accuracy, transferability, and ease of training of five atom-centered MLIPs, including the artificial neural network potentials developed by the atomic energy network (AENET), the Gaussian approximation potential (GAP), the spectral neighbor analysis potential (SNAP) and its quadratic extension (qSNAP), and the moment tensor potential (MTP), in modeling a 11-component LDR chemical space. Specifically, we generate a DFT-calculated data set of 10,842 configurations of disordered LiTMO2 and TMO2 compositions, where TM = Sc, Ti, V, Cr, Mn, Fe, Co, Ni, and/or Cu. To provide a point-of-comparison on the performance of atom-centered MLIPs, we also trained the neural equivariant interatomic potential (NequIP) on a subset of our data. Importantly, we find AENET to be the best potential in terms of accuracy and transferability for energy predictions, while MTP is the best for atomic forces. While AENET is the fastest to train among the MLIPs considered at low number of epochs (300), the training time increases significantly as epochs increase (3300), with a corresponding reduction in training errors (∼60%). Note that AENET and GAP tend to overfit in small data sets, with the extent of overfitting reducing with larger data sets. Finally, we observe AENET to provide reasonable predictions of average Li-intercalation voltages in layered, single-TM LiTMO2 frameworks, compared to DFT (∼10% error on average). Our study should pave the way both for discovering novel disordered rocksalt electrodes and for modeling other configurationally complex systems, such as high-entropy ceramics and alloys.

The suppression of heat transport in disordered crystals arises from a competition between mass fluctuations and bond disorder, but their relative contributions remain difficult to disentangle. We address this challenge using a machine-learned interatomic potential trained on ab initio data across the PbTe1-xSex alloy series. Molecular dynamics simulations with the trained machine-learned interatomic potential reproduce experimental lattice thermal conductivities and density-functional theory phonon dispersions while enabling frequency-resolved analysis of heat transport from 300 to 800 K. We find a narrow window near 1.7-2.0 THz dominates heat transport across all compositions, where heat is primarily carried by mixed longitudinal acoustic and optical modes. Alloying dramatically reduces spectral diffusivity near 2 THz leading to the deterioration of thermal conductivity. For the parent compounds both the spectral diffusivity and overall thermal conductivity decrease at elevated temperatures. However, the thermal degradation is weaker for the mid-range composition (x = 0.5) due to the greater thermal occupation of vibrational modes and increased heat capacity. Alchemical simulations show that force-constant disorder, not mass contrast, plays the dominant role in this suppression. These results highlight the microscopic mechanisms underlying thermal transport breakdown in alloys and demonstrate how machine-learned interatomic potentials now offer a tractable path toward predictive, physics-rich thermal transport modeling in complex disordered solids.

Low-dimensional hybrid organic-inorganic perovskites (HOIPs) are promising electronically active materials for light absorption and emission. The design space of HOIPs is extremely large, as a variety of organic cations can be combined with different inorganic frameworks. This not only allows for tunable electronic and mechanical properties but also necessitates the development of new tools for in silico high throughput analysis of candidate materials. In this work, we present an accurate, efficient, and widely applicable machine learning interatomic potential (MLIP) trained on 86 diverse experimentally reported HOIP materials. This MLIP was tested on 73 experimentally reported perovskite compositions and achieves a high accuracy, relative to density functional theory (DFT). We also introduce a novel random structure search algorithm designed for the crystal structure prediction of 2D HOIPs. The combination of MLIP and the structure search algorithm reliably recovers the crystal structure of 14 known 2D perovskites by specifying only the organic molecule and inorganic cation/halide. Performing this crystal structure search with ab initio methods would be computationally prohibitive but is relatively inexpensive with the MLIP. Finally, the developed procedure is used to predict the structure of a totally new HOIP with cation (cis-1,3-cyclohexanediamine). Subsequently, the new compound was synthesized and characterized, which matches the predicted structure, confirming the accuracy of our method. This capability will enable the efficient and accurate screening of thousands of combinations of organic cations and inorganic layers for further investigation.

Efficiently creating a concise but comprehensive data set for training machine-learned interatomic potentials (MLIPs) is an under-explored problem. Active learning, which uses biased or unbiased molecular dynamics (MD) to generate candidate pools, aims to address this objective. Existing biased and unbiased MD-simulation methods, however, are prone to miss either rare events or extrapolative regions—areas of the configurational space where unreliable predictions are made. This work demonstrates that MD, when biased by the MLIP’s energy uncertainty, simultaneously captures extrapolative regions and rare events, which is crucial for developing uniformly accurate MLIPs. Furthermore, exploiting automatic differentiation, we enhance bias-forces-driven MD with the concept of bias stress. We employ calibrated gradient-based uncertainties to yield MLIPs with similar or, sometimes, better accuracy than ensemble-based methods at a lower computational cost. Finally, we apply uncertainty-biased MD to alanine dipeptide and MIL-53(Al), generating MLIPs that represent both configurational spaces more accurately than models trained with conventional MD.

We present an automated procedure for computing stacking fault energies in random alloys from large-scale simulations using moment tensor potentials (MTPs) with the accuracy of density functional theory (DFT). To that end, we develop an algorithm for training MTPs on random alloys. In the first step, our algorithm constructs a set of $\ensuremath{\sim}10$ 000 or more training candidate configurations with 50--100 atoms that are representative for the atomic neighborhoods occurring in the large-scale simulation. In the second step, we use active learning to reduce this set to $\ensuremath{\sim}100$ most distinct configurations---for which DFT energies and forces are computed and on which the potential is ultimately trained. We validate our algorithm for the MoNbTa medium-entropy alloy by showing that the MTP reproduces the DFT $\frac{1}{4}[111]$ unstable stacking fault energy over the entire compositional space up to a few percent. Contrary to state-of-the-art methods, e.g., the coherent potential approximation (CPA) or special quasi-random structures (SQSs), our algorithm naturally accounts for relaxation, is not limited by DFT cell sizes, and opens opportunities to efficiently investigate follow-up problems, such as chemical ordering. In a broader sense, our algorithm can be easily modified to compute related properties of random alloys, for instance, misfit volumes, or grain boundary energies. Moreover, it forms the basis for an efficient construction of MTPs to be used in large-scale simulations of multicomponent systems.

Machine-learning interatomic potentials (MLIPs) enable large-scale atomistic simulations at moderate computational cost while retaining ab initio accuracy. In recent years, MLIPs trained on coupled-cluster data─particularly CCSD(T), which includes single, double, and perturbative triple excitations─have emerged as a promising route to achieve chemical accuracy (1 kcal/mol) beyond the limits of density functional theory (DFT) and to incorporate nonempirical van der Waals (vdW) interactions. Most existing approaches are, however, still not straightforwardly applicable for systems with extended covalent networks such as covalent organic frameworks (COFs) due to the limited availability of CCSD(T) under periodic boundary conditions. Here we present a methodology to train MLIPs with CCSD(T) accuracy for systems with extended covalent networks. The approach is based on the Δ-learning method with a dispersion-corrected tight-binding baseline and an MLIP trained on the differences of the target CCSD(T) energies from the baseline. This Δ-learning strategy enables training on compact molecular fragments while preserving transferability toward the periodic systems. Dispersion interactions are accounted for by including vdW-bound multimers in the training set, and the combination with a vdW-aware tight-binding baseline allows the formally local MLIP to attain CCSD(T)-level accuracy even for systems dominated by long-range vdW forces. The resulting potential yields root-mean-square energy errors below 0.4 meV/atom on both training and test sets and reproduces electronic total atomization energies, bond lengths, harmonic vibrational frequencies, and intermolecular interaction energies for benchmark molecular systems. We apply the method to a prototypical quasi-two-dimensional covalent organic framework (COF) composed of carbon and hydrogen. The COF structure, interlayer binding energies, and hydrogen absorption are analyzed at CCSD(T) accuracy. Overall, the developed methodology opens a practical route to large-scale atomistic simulations for systems with extended covalent networks and vdW interactions with chemical accuracy.

The search for non-toxic and low-cost perovskites with high thermoelectric performance is still a challenge despite low thermal conductivity. The thermoelectric properties of nitride anti-perovskites X3BN (B = Bi, Sb, X = Mg, Ca, Sr) with a cubic structure were investigated using density functional theory and machine learning interatomic potential. The low Debye temperature and thermal conductivity were obtained due to strong lattice anharmonicity, and the phonon vibration modes were also analyzed. The high band degeneracy and suitable bandgap lead to a large power factor. The maximum power factor is 7.54 mW/mk2 for Mg3BiN at 900 K. We obtained a maximum ZT of 1.49 for p-type Sr3BiN, and it is 1.22 for n-type doping at 900 K. The ZT for Mg3BiN is 1.18 and 1.19 for p-type and n-type doping, respectively. Our calculations reveal that these anti-perovskites are excellent materials for non-toxic and low-cost thermoelectric applications.

The accurate description of phonon dispersion of two-dimensional (2D) materials demonstrates significance in many research fields of condensed matter physics. In this paper, we systematically calculate the phonon spectra and transport properties of six representative 2D materials (encompassing single-element and binary compounds with flat, buckled, and puckered backbone geometries) by means of density functional theory (DFT) and two machine learning interatomic potentials [MLIPs, on-the-fly machine learning potential (FMLP), and moment tensor potential (MTP)]. The results show that the acoustic out-of-plane flexural (ZA) dispersion of the 2D materials are always and easily exhibiting non-quadratic dispersion phenomena near the center of the Brillouin zone by using the pure DFT calculation method. This phenomenon contradicts physics and reflects intuitively from the non-zero group velocity at Γ point. However, no matter which MLIP (FMLP/MTP) the calculation is based on, it could solve such behavior perfectly, where the ZA mode conforms to the quadratic dispersion relationship in the long-wavelength limit. Our results further demonstrate that compared to the pure DFT calculation, the FMLP and MTP method could quickly and relatively accurately obtain the lattice thermal conductivities of graphene, silicene, phosphorene, SiC, MoS2, and GeS. The findings presented in this work provide a solution about the pseudophysical phenomenon of ZA dispersions in 2D materials with the pure DFT calculation, which will greatly facilitate research areas such as phonon thermal transport, flexural mechanics, and electron–acoustic coupling.

Recent developments in machine learning interatomic potentials (MLIPs) have empowered even nonexperts in machine learning to train MLIPs for accelerating materials simulations. However, reproducibility and independent evaluation of presented MLIP results is hindered by a lack of clear standards in current literature. In this Perspective, we aim to provide guidance on best practices for documenting MLIP use while walking the reader through the development and deployment of MLIPs including hardware and software requirements, generating training data, training models, validating predictions, and MLIP inference. We also suggest useful plotting practices and analyses to validate and boost confidence in the deployed models. Finally, we provide a step-by-step checklist for practitioners to use directly before publication to standardize the information to be reported. Altogether, we hope that our work will encourage the reliable and reproducible use of these MLIPs, which will accelerate their ability to make a positive impact in various disciplines including materials science, chemistry, and biology, among others.

The development of differentiable invariant descriptors for accurate representations of atomic environments plays a central role in the success of interatomic potentials for chemistry and materials science. We introduce a method to generate fast proper orthogonal descriptors for the construction of many-body interatomic potentials and discuss its relation to exising empirical and machine learning interatomic potentials. A traditional way of implementing the proper orthogonal descriptors has a computational complexity that scales exponentially with the body order in terms of the number of neighbors. We present an algorithm to compute the proper orthogonal descriptors with a computational complexity that scales linearly with the number of neighbors irrespective of the body order. We show that our method can enable a more efficient implementation for a number of existing potentials and provide a scalable systematic framework to construct new many-body potentials. The new potentials are demonstrated on a data set of density functional theory calculations for Tantalum and compared with other interatomic potentials.