--- title: Machine-Learned Interatomic Potentials url: https://www.emergentmind.com/topics/machine-learned-interatomic-potential-ml-iap type: topic --- # Machine-Learned Interatomic Potentials Machine-learned interatomic potentials (ML-IAPs), often denoted MLIPs in the recent literature, are data-driven models of the potential energy surface used to provide computationally efficient predictions of energy and Newtonian forces for atomistic simulation. Their standard construction decomposes total energy into local atomic contributions, \(E=\sum_i E_i\), maps local atomic environments \(\mathbf{R}_i\) to descriptors \(\mathbf{G}_i\) that are invariant to translations, rotations, and permutations, and fits a regression model to a quantum-mechanical reference database, most often generated with density functional theory (DFT). Within the interpolation regime defined by that database, ML-IAPs can achieve near quantum accuracy at a computational cost orders of magnitude below DFT, while active-learning and physically informed variants are designed to reduce failures associated with extrapolation to previously unseen environments [2102.06163] [1611.09346]. ## 1. Formal structure and conceptual interpretation A canonical ML-IAP expresses the total energy as a sum of atomic energies and evaluates each atomic energy from a descriptor of the local environment. In the review literature, this is written schematically as \[ \mathbf{R}_i \rightarrow \mathbf{G}_i \overset{\mathcal{R}}{\rightarrow} E_i, \] with \(\mathcal{R}\) implemented by neural networks, Gaussian process regression, kernel ridge regression, moment tensor potentials, SNAP, or related regressors. The same literature emphasizes that the predictive ability and transferability of an ML-IAP are strictly limited to the atomic environments represented in the reference database, which is why database design, active learning, and uncertainty monitoring are central rather than auxiliary components of the field [2102.06163]. One influential conceptual interpretation is that high-accuracy ML-IAPs are not merely black-box surrogates. For elemental Ti, linearized MLIPs were presented as a generalization of the embedded atom method (EAM) and modified embedded atom method (MEAM), with the embedding energy treated beyond the uniform density approximation. In that treatment, pairwise descriptors take the form \[ b_{n0}^{(i)}=\sum_j f_n(r_{ij}), \] and angular-dependent descriptors are written as \[ b_{nl}^{(i)}=\sum_{j,k} f_n(r_{ij})f_n(r_{ik})\cos^l(\gamma_{jik}), \qquad l\ge 1. \] The central claim is that the high accuracy and transferability of MLIPs stem mainly from employing a large, systematic set of physically meaningful descriptors, especially angular-dependent ones, rather than from using sophisticated black-box functions alone [1708.02741]. A parallel line of work introduces hybrid or physically informed ML-IAPs in which the regression model predicts parameters of a physics-based potential rather than energies directly: \[ \mathbf{R}_i \rightarrow \mathbf{G}_i \overset{\mathcal{R}}{\rightarrow} \mathbf{p}_i \overset{\Phi}{\rightarrow} E_i. \] In the review formulation, this class is intended to improve transferability to unknown atomic environments while retaining the accuracy of data-driven fitting inside the database domain [2102.06163]. ## 2. Descriptor and model families The field contains several distinct but partially overlapping model classes. A large fraction of early methodological work concerns linearly parameterized or polynomial models. In the active-learning formulation for linearly parameterized potentials, the energy of a configuration \(x\) is written as \[ E(x)=\sum_{j=1}^{m}\theta_j b_j(x), \] which covers models such as Moment Tensor Potentials (MTPs) and, in the cited discussion, other linearized families including SNAP and GAP expansions [1611.09346]. For Ti, linearized MLIPs used polynomial functions of pairwise and angular descriptors; for Si–C, the UF3 formulation used explicit two-body and three-body terms expanded in cubic B-splines, \[ E(S)=\sum_{i\gamma_{\text{tsh}}\), the configuration is added to the training set. In that formulation, \(\gamma(x^*)\) is the extrapolation grade, and values with all \(|C_k|\le 1\) are interpreted as interpolation rather than extrapolation [1611.09346]. Other data-generation strategies target specific regions of the potential energy surface. For bulk crystalline matter near stationary points, an automated non-diagonal-supercell protocol was proposed to sample the small-displacement regime efficiently. Starting from a geometry-optimized primitive cell, one enumerates NDSCs for a target \(q\)-point grid, applies Gaussian random displacements with \(\sigma=0.10\) Å, isotropic volume changes, and random shear strains through \[ \mathbf{L}_{\text{rand}}=(\mathbf{I}+\boldsymbol{\epsilon})\mathbf{L}, \] with \(\boldsymbol{\epsilon}\) drawn componentwise from \(\mathcal{U}(-0.01,0.01)\). This procedure was shown for Al, W, Mg, and Si, yielding phonon RMSE below \(0.5\) THz and elastic properties often within \(1\)–\(6\%\) of DFT [2207.11828]. For chemically diverse systems, global structure search can replace manually curated enumeration. The Si–C UF3 model was trained predominantly on structures discovered by the GASP genetic algorithm over the entire silicon–carbon composition space, including both relaxed and unrelaxed structures, producing 41,145 GA structures plus 80 AIMD structures [2403.15897]. A different route is indirect learning: an accurate but slower teacher ML-IAP, such as GAP, generates both configurations and labels for a much larger secondary dataset, which is then used to train a faster student model such as MTP, without new quantum-mechanical calculations at the secondary stage [2111.11120]. Recent work also formalizes data selection as a cost–accuracy optimization problem. In an application-specific qSNAP study for Be, 20,000 configurations were generated by information entropy maximization, and leverage-score sampling based on the hat matrix \[ H=X(X^T X)^{-1}X^T \] was used to identify informative training subsets. That study found that leverage sampling can reduce the required training-set size by up to an order of magnitude compared with random sampling, and that the optimal DFT precision should be matched to the intrinsic precision limit of the target ML-IAP rather than maximized indiscriminately [2506.05646]. ## 4. Training objectives, benchmarking practice, and uncertainty quantification Most modern ML-IAPs are trained against energies, forces, and often stresses or virials simultaneously. The review literature gives a representative loss combining these terms with regularization, while application-specific studies often use weighted least squares. In the Be qSNAP analysis, the loss was written as \[ L=\sum_{m=1}^{M}\left\{\frac{w_E^2(\hat{E}_m-E_m)^2}{N_m^2} +\sum_{i=1}^{3N_m}w_F^2(\hat{F}_{mi}-F_{mi})^2\right\}, \] and one of the main findings was that, when lower-precision DFT data are used, increasing the force weight can lower both force and energy errors on high-precision test data because forces converge faster than total energies [2506.05646]. Model-specific fitting procedures differ substantially. Linearized Ti MLIPs used ridge regression on energies and forces [1708.02741]. Carbon SNAP at extreme conditions used weighted linear regression through FitSNAP, with group weights optimized by a genetic algorithm in DAKOTA and final hyperparameters \(r_{\text{cut}}=2.7\) Å and \(J=4\) [2205.01209]. The aluminum SNAP model used FitSNAP and DAKOTA to optimize both cutoff radius and training weights, yielding an optimal \(R_{\text{cut}}=5.323\) Å and \(2J_{\max}=6\) [2304.09703]. BCM training for CNH molecules used SGPR experts and combined committee predictions by uncertainty-weighted aggregation, while the covariance loss \[ s(\rho_i)=K(\rho_i,\rho_i)-\bm{k}_{\rho_i m}\bm{k}_{mm}^{-1}\bm{k}_{\rho_i m}^T \] served as an uncertainty indicator for active sampling [2402.17132]. Validation practice has broadened beyond train/test RMSE. In bcc Fe, a three-step protocol evaluated ML-IAPs by Pareto fronts of computational speed versus RMSE, a benchmark-property quality factor \[ Q=\frac{1}{N_p}\sum_{i=1}^{N_p}\left|\frac{f_i^{ML}-f_i^{DFT}}{f_i^{DFT}}\right|, \] and per-atom uncertainty during extended-defect simulations. That study found that \(Q\) correlates with RMSE, but reductions in RMSE below about \(5\) meV/atom do not necessarily improve benchmark-property accuracy, and it explicitly argued for indicators beyond RMSE when assessing transferability [2307.10072]. Uncertainty quantification is likewise model-family dependent. In D-optimal active learning, the extrapolation grade \(\gamma(x^*)\) can be computed without QM data and was reported to be strongly correlated with model error [1611.09346]. In Fe defect studies, safe regimes were defined by per-atom uncertainty below \(15\) meV/atom for GAP and by per-atom extrapolation grade below \(3\) for ACE-derived models [2307.10072]. This benchmarking practice has become particularly important when the intended applications involve dislocations, fracture, shock loading, or warm dense matter rather than equilibrium crystalline properties alone. ## 5. Accuracy, transferability, and representative domains of application The most direct demonstrations of ML-IAP capability come from property-specific case studies. For elemental Ti, a pairwise MLIP trained on 2,700 DFT configurations and about 430,000 data points achieved an energy RMSE of \(3.8\) meV/atom with 2,925 coefficients, whereas the angular-dependent MLIP reached \(0.5\) meV/atom with 35,245 coefficients, and the latter accurately reproduced elastic constants, bulk moduli, and phonon dispersion curves, unlike EAM/MEAM and unlike the pairwise model on some elastic properties [1708.02741]. For extreme-condition carbon, a quadratic SNAP model was trained on a QMD database covering pressures from \(0\) to \(5\) TPa, temperatures from \(0\) to \(40{,}000\) K, densities from \(2.9\) to \(13.6\) g/cm\(^3\), and 636 structures comprising 124,907 unique atomic environments. In two-phase MD, the resulting melting lines for diamond, BC8, and SC were within about \(200\) K, or \(\sim 3\%\), of QMD across \(0\)–\(5{,}000\) GPa, and the model enabled simulations of 1 million-atom polycrystalline diamond [2205.01209]. A different carbon study used GAP-20 with nested sampling up to \(1\) TPa and found remarkably good agreement with experiment up to \(200\) GPa, but also identified nonphysical graphite layer spacings and high-pressure phases not confirmed by DFT, making it a canonical example of both impressive transferability and clear extrapolation failure [2208.09692]. For aluminum, a single SNAP ML-IAP was trained on DFT-MD data between \(300\) and \(10{,}000\) K at ambient density \(2.7\) g/cm\(^3\), using 108- and 256-atom cells and about 100 snapshots per temperature. The reported prediction errors were \(6.27\) meV/atom MAE and \(8.7\) meV/atom RMSE for energies, and \(0.22\) eV/Å MAE and \(0.32\) eV/Å RMSE for forces. Despite the training range ending at \(10{,}000\) K, the model was used to examine thermal conductivity, diffusion coefficient, viscosity, sound velocity, and ion-ion structure factors up to about \(58{,}000\) K, and diffusion trends were checked up to \(116{,}000\) K [2304.09703]. Transferability across chemical families has likewise been demonstrated in molecular systems. The BCM potential for organic nitrogen compounds partitioned CNH data into chemical-family experts such as amines, azoles, alkaloids, cyanides, hydrazines, imidazoles, and nitriles; it then showed high consistency with DFT on long-chain amines such as \(\ce{C21NH45}\) even though training used only short-chain amines. Reported test-set energy MAEs ranged from \(0.036\) meV for azoles to \(0.242\) meV for alkaloids, and force MAEs ranged from \(0.091\) to \(0.120\) eV/Å across the listed groups [2402.17132]. Defect mechanics has been a particularly stringent proving ground. In bcc Fe, GAP, PACE, MTP, and neural-network models were compared for dislocations and cracks. All optimized models reached energy RMSE \(\le 5\) meV/atom; mean relative benchmark-property errors \(Q<3\%\) were reported; the screw-dislocation Peierls barrier was predicted in the range \(48\)–\(54\) meV per Burgers vector \(b\), compared with DFT values of \(48\)–\(58\) meV/\(b\); and cleavage on the pre-existing crack plane was found to be the zero-temperature atomistic fracture mechanism of pure bcc Fe under mode-I loading, independent of ML package and training database [2307.10072]. Large-scale materials applications extend beyond equilibrium thermodynamics. A SNAP ML-IAP for Ni–Mo alloys was used in LAMMPS to simulate polycrystalline Ni and Ni-10%Mo with average grain diameters from \(4.1\) to \(11.6\) nm and system sizes up to about 454,600 atoms. It reproduced Hall–Petch and inverse Hall–Petch behavior, linked strengthening at very small grain size to reductions in atomic von Mises strain by solute doping and annealing, and showed that annealing can weaken large-grain polycrystals because of large atomic movements in grain-boundary regions [2302.08698]. For Si–C, the UF3-based model was reported to capture peritectic decomposition of SiC and carbon dimer formation during SiC surface reconstruction, phenomena that the paper states cannot be captured with prior classical interatomic potentials for that material [2403.15897]. ## 6. Software ecosystems, scalability, and universal deployment Methodological progress has been accompanied by increasingly integrated software stacks. The active-learning paper that introduced D-optimal selection published software, test cases, and usage examples at the `mlip` repository associated with MTP fitting and on-the-fly training [1611.09346]. SNAP development for carbon and aluminum relied on FitSNAP, with DAKOTA used for hyperparameter or weight optimization [2205.01209] [2304.09703]. The BCM workflow for CNH compounds integrated SGPR models with molecular dynamics through the AutoForce Python package [2402.17132]. The NDSC phonon-data workflow was accompanied by scripts and datasets intended for reproducibility and insertion into broader MLIP pipelines [2207.11828]. A more explicit unification effort appears in the 2025 MLIP library white paper. That library is written entirely in JAX, exposes three architectures—MACE, NequIP, and ViSNet—supports ASE and JAX-MD wrappers, and provides pre-trained organics models. On SPICE2-derived data, models were trained for 220 epochs with reported training times of 245 h for NequIP, 158 h for ViSNet, and 266 h for MACE on NVIDIA H100 GPUs. Runtime benchmarks reported, for example, \(6.3\) ms/step for MACE-large on 1UAO in JAX+JAX-MD versus \(44.2\) ms in Torch+ASE, and the paper summarized this as an approximately \(5\)–\(10\times\) speed-up for JAX+JAX-MD over Torch+ASE [2505.22397]. AI2Pot extends this unification toward high-throughput training and large-scale MD deployment. It supports MTP and NEP within a PyTorch/PyTorch-Lightning frontend, but re-engineers descriptor construction, energy, force, and virial evaluation as hand-crafted C++/CUDA operators shared between training and inference. Models can be deployed in ASE and LAMMPS, and benchmarked throughput on a single NVIDIA RTX 4090 reached 923,500 atom-steps/s for NEP and 355,300 atom-steps/s for MTP at 3.4 million atoms. Reported full training times on a Ge–Sb–Te dataset were 1.3 h for NEP and 0.75 h for MTP with batch size 32 [2607.06969]. At the universal-model scale, MACE-Osaka26 merged HE26 heavy-element data with MPtrj and OFF23 using total-energy alignment and produced a 97-element model trained on 32 NVIDIA A100 GPUs across 4 nodes. Reported test errors included an MPtrj energy MAE improvement from 33.7 to 29.3 meV/atom and an OFF23 RMSE reduction from 25.2 to 12.3 meV/atom relative to the cited Osaka24 baseline, while heavy-element compounds in HE26 achieved \(44.7\) meV/atom energy MAE and \(26.3\) meV/Å force MAE [2603.03223]. This development places universal ML-IAPs alongside application-specific models as a distinct software and data-management paradigm. ## 7. Limitations, controversies, and evolving scope The central limitation of ML-IAPs remains extrapolation. The 2021 review states the problem directly: outside the database domain, extrapolation to atomic environments not present or underrepresented in the training set leads to unpredictable or even unphysical results, because no underlying physical model constrains the potential in the way traditional analytic forms do [2102.06163]. Carbon phase-diagram studies illustrate this concretely: GAP-20 reproduced graphite, diamond, and melting behavior well up to about \(200\) GPa despite being trained mainly at zero pressure, yet at higher pressure it produced multiple minima in graphite interlayer spacing and stable strained phases not confirmed by DFT [2208.09692]. A second controversy concerns what counts as sufficient validation. Several studies warn that low train/test errors are not enough. The Fe defect paper showed that significant quantitative predictions of dislocation glide paths and \(K_{Ic}\) can be sensitive to database, ML package, and cutoff radius, and argued that RMSE must be complemented by benchmark properties and uncertainty indicators [2307.10072]. The Be qSNAP study similarly argued that MLIP quality depends jointly on training-set selection, reference-data precision, energy-versus-force weighting, and model complexity, and used Pareto analysis rather than a single accuracy measure to identify efficient operating points [2506.05646]. A third controversy concerns whether ML-IAPs are restricted to interpolation-only use. Typical usage guidelines emphasize exhaustive training sets and caution against applying models to situations not represented in training. A 2025 study explicitly challenged that position by constructing SOAP/GAP and Allegro models from minimalist on-the-fly datasets—247 structures for BaTiO\(_3\), 488 for BiFeO\(_3\), 1300 for PbZrO\(_3\), and 463 for HfO\(_2\), all with default or near-default hyperparameters—and showing qualitatively and quasi-quantitatively correct predictions for vortex–antivortex electric textures, ferroelectric switching paths, competing polymorphs, and the oIII ferroelectric HfO\(_2\) polymorph. The same study also noted that accuracy drops for states very far from the training domain [2511.17449]. This suggests that the operative question is not simply whether extrapolation occurs, but how aggressively it is monitored, constrained, and validated. The likely long-term direction is therefore plural rather than singular. Active learning based on extrapolation grades, physically informed hybrids that couple ML with physics-based potentials, optimal data-generation schemes such as NDSC or leverage sampling, distilled teacher–student workflows, universal many-element foundation-style models, and high-performance deployment frameworks all address different points on the same design surface: accuracy, transferability, interpretability, and computational cost [1611.09346] [2102.06163] [2111.11120] [2603.03223]. The modern ML-IAP field is best understood not as a single method, but as a family of atomistic surrogate models whose success depends on matching representation, database, uncertainty control, and simulation target with unusual precision. Source: https://www.emergentmind.com/topics/machine-learned-interatomic-potential-ml-iap