Universal Machine-Learned Interatomic Potential (u-MLIP)

• Universal MLIPs are pre-trained, transferable models that decompose atomic contributions to predict energies and forces with near-DFT fidelity.
• They employ diverse architectures, including kernel regression and graph neural networks, to navigate broad chemical and configurational spaces.
• Fine-tuning on task-specific datasets enhances reliability, enabling applications from crystal structure prediction to high-pressure simulations.

Searching arXiv for recent u-MLIP literature to ground the article in current papers. Search query: universal machine learning interatomic potentials arXiv 2025 2026 Universal machine-learned interatomic potentials (u-MLIPs) are pre-trained, data-driven interatomic models intended to approximate the Born–Oppenheimer or density-functional-theory potential-energy surface across broad chemical and configurational space with a single transferable model. They differ from specialized MLIPs, which are trained on one material or a narrow class of systems, and from classical force fields, which rely on fixed analytic forms; in current practice, u-MLIPs are used both zero-shot and as foundation models for subsequent fine-tuning on smaller task-specific datasets. Across crystals, molecules, surfaces, defects, frameworks, and disordered phases, the central promise is near-DFT energies and forces at speedups reported from roughly $10^3$ to $10^6$ relative to direct electronic-structure calculations, although benchmark studies increasingly emphasize that this promise is strongly regime-dependent (Han et al., 2 Jun 2025, Liu et al., 9 Jun 2025, Shinagawa et al., 9 Mar 2026).

1. Concept and mathematical structure

A u-MLIP is commonly defined as a single surrogate model that predicts energies and forces across multiple elements, bonding motifs, crystal structures, and thermodynamic conditions without system-specific reparameterization. In one explicit benchmark formulation, the desired criteria of “universality” are transferability, near-equilibrium fidelity, global PES generalization, and numerical robustness (Tahmasbi et al., 23 Dec 2025).

The formal core of most u-MLIPs is an atomwise decomposition of the total energy. A representative expression used across multiple studies is

$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$

where $D_i$ is a descriptor of the local environment of atom $i$, and $\Theta$ denotes learnable parameters. Forces are obtained by differentiation,

$$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$

The same basic structure appears in graph models, tensor networks, and other local-environment formalisms (Han et al., 2 Jun 2025, Loew et al., 25 Aug 2025).

This decomposition encodes a specific compromise between locality and transferability. The model learns atomic energy contributions from neighborhoods defined by species, distances, angles, and related geometric information, then reconstructs the system energy by summation. A plausible implication is that most questions about universality ultimately reduce to how well the model’s local environment representation spans the relevant configuration manifold.

2. Architectural families and atomic representations

Current u-MLIPs span at least two broad implementation families. One family uses kernel regression over global or local descriptors, with a form written as

$$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$

where $K$ is a positive-definite kernel on structural descriptors. The second family uses message-passing or graph neural networks, including M3GNet, MACE, and eSEN, where node features are initialized from atomic species, edge features are built from distance basis functions and optional angular channels, and iterative message updates generate the final atomic embedding from which $E_i$ is read out (Loew et al., 25 Aug 2025).

Descriptor design follows a recurring pattern. Atomic species are mapped to learned embeddings; radial information is expanded with basis functions such as Gaussians or Chebyshev polynomials up to a cutoff; angular descriptors or three-body terms encode bonding geometry; and the structure is cast as a neighbor graph with edges inside a cutoff radius. In benchmarked public models, these design choices appear as SOAP in ORB, ACSF in SevenNet variants, message-passing features in MACE, GRACE, and MatterSim, local density channels in CHGNet, and graph-transformer or equivariant backbones in models such as EquiformerV2 (Han et al., 2 Jun 2025).

Architectural specialization is substantial even within the common atomwise-energy template. MACE constructs high-body-order equivariant features by tensor products of radial bases and spherical harmonics and uses only two message-passing layers because each layer already encodes up to 4-body interactions (Liu et al., 9 Jun 2025). PFP v8 is built on TeaNet, an equivariant graph neural network that propagates scalars, vectors, and rank-2 Euclidean tensors over a radius $10^6$0 nm and uses Bader-charge labels as auxiliary targets during training (Shinagawa et al., 9 Mar 2026). UMA uses an equivariant graph neural network with a Mixture-of-Linear-Experts architecture, in which a router activates only a sparse subset of experts at inference (Kwon et al., 19 Jun 2026). These differences are not merely cosmetic: they alter angular resolution, short-range behavior, memory use, and the balance between numerical smoothness and PES fidelity.

3. Data regimes, objectives, and transfer across fidelity levels

The training ecology of u-MLIPs is heterogeneous. Representative datasets include MPtrj with $10^6$1 million DFT calculations and 89 elements, a 2021.2.8 Materials Project release with $10^6$2 million relaxed structures, a pressure benchmark based on 162k relaxed Alexandria structures expanded to $10^6$3 million DFT single-point calculations from $10^6$4 to $10^6$5 GPa, and PFP v8 training corpora comprising 3 million r2SCAN structures across 70 elements, 60 million PBE/+U structures across 96 elements, 6 million $10^6$6B97X-D molecules, and $10^6$7 OC20 slab–adsorbate structures (Loew et al., 25 Aug 2025, Shinagawa et al., 9 Mar 2026). Expansion to heavier chemistry has also been demonstrated: MACE-Osaka26 merges HE26, OFF23, and MPtrj to reach 97 elements, with the heavy-element subset covering Am, Cm, Cf, Fr, At, Ra, Po, and Rn (Kuroda et al., 3 Mar 2026).

Despite architectural diversity, training objectives are remarkably standardized. Losses almost always combine energy, force, and often stress errors; some models additionally include magnetic moments or atomic charges. Typical forms are weighted mean-squared or Huber/L1 hybrids. CHGNet fine-tuning across functionals uses a weighted Huber-type combination of energy, force, stress, and magnetic moment with $10^6$8, while pressure fine-tuning of eSEN uses an $10^6$9 loss on energies, forces, and virial stresses (Huang et al., 7 Apr 2025, Loew et al., 25 Aug 2025).

Dataset composition has direct consequences for downstream reliability. In MOFSimBench, MPtraj-only training systematically softens the PES, whereas OMat24-based training with extensive out-of-equilibrium sampling reduces softening and improves structural optimization, MD stability, and bulk-property prediction (Kraß et al., 16 Jul 2025). A plausible implication is that universality is constrained at least as much by sampling policy as by network expressivity.

Transfer from one electronic-structure fidelity to another is nontrivial. Within CHGNet, raw single-point GGA and r$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$0SCAN energies show a Pearson correlation of only $$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$1, but after removing each functional’s fitted atomic reference the correlation increases to $$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$2. On the 0.24 M-structure MP-r$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$3SCAN set, replacing the source AtomRef with the target r$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$4SCAN AtomRef before transfer learning reduces test errors to 17 meV/atom for energy and 38 meV/$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$5, compared with 27 meV/atom and 45 meV/$$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$6 for scratch training (Huang et al., 7 Apr 2025). This identifies elemental energy referencing as a central technical issue in cross-functional u-MLIPs.

4. Benchmarked performance and scientific use-cases

Benchmark results show that u-MLIPs are already useful across multiple scientific workflows, but the relevant metric depends strongly on task. In phonon and neutron-scattering applications, the best models achieve sub-meV mean phonon errors; in metals and alloys, pretrained defect benchmarks reach low-meV/atom energy errors; in crystal-structure prediction, u-MLIPs materially shift the computational bottleneck from force evaluation to search; and in nanoporous materials they can exceed both classical force fields and fine-tuned MOF-specific baselines on several tasks (Han et al., 2 Jun 2025, Shuang et al., 5 Feb 2025, Kraß et al., 16 Jul 2025, An et al., 3 Feb 2026).

Domain	Representative result	Source
Phonons and neutron scattering	ORB-v3 phonon frequency MAE $$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$7 meV; a 100-atom phonon bandstructure + DOS takes $$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$8 s on one GPU vs. $$E_{\mathrm{tot}}(\{R\})=\sum_{i=1}^{N} E_i(D_i;\Theta),$$9 min with plane-wave DFT	(Han et al., 2 Jun 2025)
Defects in metals and random alloys	EquiformerV2 (86 M, omat+mp+salex) reaches 3.0 meV/atom energy RMSE and 60 meV/$D_i$0 force RMSE	(Shuang et al., 5 Feb 2025)
Nanoporous materials	eSEN-OAM gives bulk-modulus MAE $D_i$1 GPa; ORB-v3 gives heat-capacity MAE $D_i$2	(Kraß et al., 16 Jul 2025)
Crystal structure prediction	In one case, DFT-CSP costs $D_i$3 CPU-h, whereas M3GNet-CSP takes 37 CPU-h; seven new compounds are reported with $D_i$4 meV/atom	(An et al., 3 Feb 2026)
r2SCAN-level experimental agreement	PFP–R2SCAN gives crystal formation-energy MAE $D_i$5 eV/atom and melting-point MAE $D_i$6 K	(Shinagawa et al., 9 Mar 2026)

For lattice dynamics, the phonon benchmark over nearly 5,000 inorganic crystals finds ORB-v3, SevenNet-MF-ompa, GRACE-2L-OAM, and eSEN-30M-OAM all below 1 meV mean phonon-frequency MAE, with ORB-v3 at 0.50 meV and phonon DOS Spearman correlation $D_i$7. The same study reports a real-time fitting loop for inelastic neutron scattering that completes in $D_i$8 min per dataset, rather than hours if phonons are recomputed by DFT (Han et al., 2 Jun 2025).

For defect-rich metallic systems, state-of-the-art pretrained u-MLIPs can match or surpass specialized MLIPs on broad benchmarks. In metals and random alloys, the latest EquiformerV2 models achieve root mean square errors below 5 meV/atom for energies and 100 meV/$D_i$9 for forces, and specifically 3.0–3.2 meV/atom and 60–65 meV/$i$0 on representative defect datasets. The tungsten case study is especially notable: only 514 W configurations were present in pretraining, yet EquiformerV2 still achieved RMSE(E) $i$1 meV/atom and RMSE(F) $i$2 meV/$i$3 on a complex W defect genome (Shuang et al., 5 Feb 2025).

For nanoporous materials, MOFSimBench evaluates more than 20 models on structural optimization, MD stability, bulk properties, and host–guest interactions. Top u-MLIPs such as eSEN-OAM, ORB-v3-omat, MatterSim, and SevenNet-ompa outperform classical force fields and match or exceed a MOF-specific fine-tuned baseline on geometry, dynamics, bulk moduli, heat capacities, and adsorption energetics. In particular, top models converge structures more than 90% of the time, predict bulk moduli within 3 GPa and heat capacities within $i$4, and achieve sub-1 kJ/mol interaction-energy error for CO$i$5/H$i$6O host–guest screening (Kraß et al., 16 Jul 2025).

For structure prediction, M3GNet-driven CSP rediscovered experimentally known quaternary oxides absent from training and identified seven new thermodynamically and dynamically stable compounds. The same study emphasizes that semilocal-PBE stability predictions require cross-validation with SCAN, R2SCAN, and RPA, because PBE misranks at least some close polymorphs, including a 71 meV/f.u. misordering for two $i$7 polymorphs (An et al., 3 Feb 2026).

A further development is the shift from PBE-trained universality toward higher-fidelity experimental alignment. PFP v8 argues that better zero-shot agreement with experiment must be an explicit design target and reports systematic gains from r2SCAN-level training, including formation enthalpy, molecular GMTKN55, surface-energy, and melting-point improvements over PBE-based variants (Shinagawa et al., 9 Mar 2026).

5. Reliability limits, out-of-distribution behavior, and fine-tuning practice

A common misconception is that “universal” implies uniformly reliable across all thermodynamic regimes. The pressure benchmark makes the opposite point with unusual clarity. Across $i$8–$i$9 GPa, most ambient-trained u-MLIPs degrade substantially: volume MAE at $\Theta$0 GPa ranges from 0.05 to 0.42 $\Theta$1/atom, but at $\Theta$2 GPa it ranges from 0.08 to 1.56 $\Theta$3/atom; the best untuned models at $\Theta$4 GPa give $\Theta$5–6 meV/atom energy MAE, whereas at $\Theta$6 GPa the untuned ensemble degrades to 41–347 meV/atom. The study attributes this primarily to a data limitation: first-neighbor distances shift from $\Theta$7 $\Theta$8 at ambient conditions to $\Theta$9 $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$0 at 150 GPa, moving the task into an out-of-distribution compressed regime (Loew et al., 25 Aug 2025).

Analogous limits appear in other domains. For surface energies, out-of-the-box MACE, CHGNet, and M3GNet underpredict $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$1 with errors clustered around RMSE $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$2–0.6 J m$$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$3, and error magnitude correlates with descriptor-space distance from the bulk-training cluster (Focassio et al., 2024). For high-temperature MOF chemistry, ORB-v3 and fairchem OMAT yield the lowest static errors among five tested uMLIPs, yet long-timescale MD reveals a generative weighted loss at 2000 K that is 3–4$$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$4 larger than the AIMD validation loss, with bond breaking near 1000 K triggering a sharper error jump than static validation suggests (Edwards et al., 28 Apr 2026).

Reactive and transport observables remain especially demanding. Across seven chemically diverse systems, one benchmark finds that strong zero-shot performance on standard energy and force metrics does not guarantee accurate reactive, transport, or high-barrier observables. In the sulfur-vacancy jump of MoS$$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$5, the best zero-shot model, PET-OAM-XL, predicts a barrier of $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$6 eV instead of the DFT reference 0.8 eV (Hänseroth et al., 22 Jun 2026).

These failures have made fine-tuning a central component of u-MLIP practice. For MACE-based foundation models, fine-tuning on task-specific data often converges 2–5$$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$7 faster than scratch training and can outperform models trained from scratch when dataset selection is careful. In silicon, filtering a 500k-point dataset to $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$8k relevant points reduced energy MAE from a plateau of $$F_i=-\frac{\partial E_{\mathrm{tot}}}{\partial R_i}.$$9–30 meV/atom to $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$0 meV/atom and force MAE from $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$1 meV/$$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$2 to $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$3 meV/$$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$4; in a high-entropy alloy, fine-tuning reduced energy RMSE from $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$5 to $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$6 meV/atom and force RMSE from $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$7–130 to $$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$8 meV/$$E_{\mathrm{tot}}(R)=\sum_{j\in \mathrm{train}} \alpha_j K\bigl(\mathcal{S}(R),\mathcal{S}(R_j^{\mathrm{train}})\bigr),$$9 in under 50 epochs (Liu et al., 9 Jun 2025).

The general practical workflow is now standardized. A step-by-step tutorial based on MACE-MP-0 recommends dataset construction by rattle, strain, and short ab initio MD or Monte Carlo; uncertainty filtering if a large MD pool is available; universal losses on energy, force, and optional stress; Adam or AMSGrad with weight decay $K$0, gradient clipping near 10, learning rates in the $K$1–$K$2 range, and EMA/SWA for stabilization. The same tutorial reports representative fine-tuning results for a solid-state electrolyte, stacking-fault defects in Mo, and a graphene–water interface (Liu et al., 27 Jun 2025).

A stronger reinterpretation treats the u-MLIP not as the final predictive model but as a configuration-space generator. On this view, long u-MLIP MD is used to generate trajectories, a few thousand frames are recalculated with DFT, and a compact material-specific MLIP is then trained or fine-tuned. Across several systems, $K$3 DFT single points are often sufficient; in the MoS$K$4 vacancy-jump case, an iterative self-training loop recovers the DFT barrier with only 600 first-principles calculations in total, and the full workflow can produce 1 ns ab initio-quality trajectories within three days (Hänseroth et al., 22 Jun 2026).

6. Emerging directions in universality

One major direction is movement beyond PBE-level targets. Cross-functional transfer work within CHGNet shows that naive transfer from GGA/GGA+$K$5 to r$K$6SCAN is impaired by energy-scale shifts, whereas proper atomic referencing restores high correlation and substantial data efficiency (Huang et al., 7 Apr 2025). PFP v8 generalizes this logic into a full uMLIP program trained directly on r2SCAN and reports systematically improved agreement with experimental data or high-accuracy references for crystals, molecules, and surfaces (Shinagawa et al., 9 Mar 2026).

A second direction is simultaneous expansion of scale and scope. MACE-Osaka26 extends open u-MLIP coverage to 97 elements, including minor actinides and related heavy elements crucial for nuclear applications, and demonstrates promising accuracy on HE26 alongside solid performance on MPtrj and OFF23 (Kuroda et al., 3 Mar 2026). At the opposite end of the efficiency spectrum, SevenNet-Nano shows that universality can be distilled into a much smaller student model: it uses $K$7 trainable parameters versus $K$8 in SevenNet-Omni, attains validation MAEs of 23–27 meV/atom for energy and 0.078–0.083 eV/$K$9 for forces, and yields speedups of $E_i$0–3.0 below 1,000 atoms and up to $E_i$1 at $E_i$2 atoms (Oh et al., 13 Apr 2026).

A third direction concerns representation learning itself. Cross-model latent-feature analysis finds that different u-MLIPs encode chemical space in markedly distinct ways: off-diagonal average GFRE is $E_i$3 and average LFRE is $E_i$4 across MACE-MP-03b, PET-MAD, DPA-3.1, and UMA-S-1P1 on sampled MAD environments. Fine-tuned variants remain almost perfectly reconstructible from the pre-trained latent space, showing that fine-tuning retains a strong pre-training bias. The same work proposes progressive cumulants for compressing atom-level features into structure-level descriptors, with higher-order cumulants contributing strictly additional information up to at least order eight (Chorna et al., 5 Dec 2025).

A fourth direction is better benchmark design. On elemental systems, a dual-axis benchmark combining equation-of-state tests with minima hopping shows that search efficiency and structural fidelity can decouple: smoother learned PESs do not necessarily imply more accurate energetic landscapes. The same benchmark identifies persistent weaknesses for alkali and alkaline-earth metals even when transition-metal performance is strong (Tahmasbi et al., 23 Dec 2025). In nanoporous materials, MOFSimBench further argues that data quality—especially diversity and the inclusion of out-of-equilibrium conformations—matters more than architecture across the evaluated models, with a reported Pearson correlation of $E_i$5 between OMat24 inclusion and bulk-modulus error versus $E_i$6 between architecture class and the same error (Kraß et al., 16 Jul 2025).

Taken together, these developments suggest that the term “universal” is best understood as an evolving research target rather than a completed property. Present u-MLIPs already function as transferable, high-throughput surrogates for many atomistic tasks, but the strongest evidence in current literature indicates that robust universality depends on three coupled conditions: broader coverage of non-equilibrium and extreme-state data, careful treatment of target fidelity and energy referencing, and benchmark protocols that test not only near-equilibrium errors but also transport, reactivity, stress, high-pressure compression, and long-time generative stability.