Does Hessian Data Improve the Performance of MachineLearning Potentials?

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.

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.

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.

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.

In nuclear and/or solar applications that involve molten salts, impurities frequently enter the salt as either fission products or via corrosion. Impurities can interact and make complexes, but the impact of such complexation on the properties of the salts and corrosion rates has not been understood. Common impurities in molten salts, such as FLiBe, include Cr, Ni, and Fe. Here, we investigate the complexation of Cr and Ni in FLiBe using molecular dynamics based on a machine learning interatomic potential (MLIP) fitted using the atomic cluster expansion (ACE) method. The MLIP allows us to overcome the challenges of simultaneously needing accurate energetics and long time scale to study complexation. We demonstrate that impurity behavior is more difficult to capture than that of concentrated elements with MLIPs due to less sampling in training data, but that this can be overcome by using active learning strategies to obtain a robust fit. Our findings suggest that there is a weak but potentially significant binding free energy between CrF2 and NiF2 in eutectic FLiBe of -0.112eV. Under certain conditions this binding creates a significant short-range order between the species and lowers the redox potential of NiF2 in the presence of CrF2 in FLiBe, making Ni dissolution more favorable in the presence of Cr as compared to its dissolution in pure FLiBe. However, we find little impact of this complexation on the diffusivity of Ni and Cr. Overall, the methodology presented here suggests an approach to modeling complexation with MLIPs and suggests that interactions between dissolved cations could be playing a significant role in some salt thermophysical properties.

Variational Monte Carlo (VMC) can be used to train accurate machine learning interatomic potentials (MLIPs), enabling molecular dynamics (MD) simulations of complex materials on time scales and system sizes previously unattainable. VMC training sets are often based on partially optimized wave functions (WFs) to circumvent expensive energy optimizations of the whole set of WF parameters. However, frozen variational parameters lead to VMC forces and pressures not consistent with the underlying potential energy surface, a bias called the self-consistency error (SCE). Here, we demonstrate how the SCE can spoil the accuracy of MLIPs trained on these data, taking high-pressure hydrogen as the test case. We then apply a recently introduced SCE correction [Phys. Rev. B 2024 109, 205151] to generate unbiased VMC training sets based on a Jastrow-correlated single determinant WF with frozen Kohn-Sham orbitals. The MLIPs generated within this framework are significantly improved and can approach in quality those trained on data sets built with fully optimized WFs. Our conclusions are further supported by MD simulations, which show how MLIPs trained on SCE-corrected data sets systematically yield more reliable physical observables. Our framework opens the possibility of constructing extended high-quality training sets with VMC.

Accurate modeling of aqueous metal salt solutions is essential for understanding processes relevant to environmental safety, energy storage, and separation technologies. Trace metals such as As3+ at low concentrations pose significant health and environmental risks. They are challenging to simulate due to limitations in both classical force fields (CFFs), which lack accuracy, and ab initio methods, which are restricted to short trajectories. In this study, we develop machine-learned interatomic potentials (MLIPs) to model aqueous AsCl3 and MgCl2 using the NequIP/Allegro equivariant graph neural network architecture trained on ab initio molecular dynamics (AIMD) and density functional theory data. Our MLIP models accurately reproduce ab initio energies and forces while capturing solvation structure, ion diffusion, and hydration dynamics more effectively than CFFs (AMBER and UFF models). Our MLIPs achieve energy MAEs < 1 meV/atom and force RMSEs < 40 meV/Å, while providing an O(104) speedup over AIMD. These MLIPs offer a reliable and efficient alternative for modeling trace metal speciation and transport, with implications for improved separation and environmental processes.

Ovonic threshold switching (OTS) selectors are pivotal in nonvolatile memory devices due to their nonlinear electrical characteristics and polarity-dependent threshold voltages. However, the atomic-scale origins of the defect states responsible for these behaviors remain unclear. In this study, we systematically investigate defects in amorphous GeSe using molecular dynamics simulations accelerated by machine learning interatomic potentials (MLIPs). We first benchmark several MLIP architectures, including descriptor-based potentials and graph neural network (GNN)-based potentials. Our results demonstrate that capturing higher-order interactions, at least four-body correlations, and medium-range structural order is essential for accurately representing amorphous GeSe structures. Our analysis indicates that GNN architectures with multiple interaction layers effectively capture higher-order correlations and medium-range order, thereby preventing spurious defects easily introduced by less descriptive MLIPs. Utilizing the optimized GNN model, we identify two distinct defect motifs across 20 independent 960-atom amorphous GeSe structures: aligned Ge chains associated with defect states near the conduction band, and overcoordinated Ge chains linked to defect states near the valence band. Moreover, we establish correlations between electronic defect levels and specific structural features─namely, the average alignment of bond angles in aligned chains and the degree of local Peierls distortion around overcoordinated Ge atoms. These insights provide a theoretical framework for interpreting experimental observations and deepening the understanding of defect-driven OTS phenomena in amorphous GeSe.

Machine learning interatomic potentials (MLIPs) enables molecular dynamics (MD) simulations with ab initio accuracy and has been applied to various fields of physical science. However, the performance and transferability of MLIPs are limited by insufficient labeled training data, which require expensive ab initio calculations to obtain the labels, especially for complex molecular systems. To address this challenge, we design a novel geometric structure learning paradigm that consists of two stages. We first generate a large quantity of 3D configurations of target molecular system with classical molecular dynamics simulations. Then, we propose geometry-enhanced self-supervised learning consisting of masking, denoising, and contrastive learning to better capture the topology and 3D 1 geometric information from the unlabeled 3D configurations. We evaluate our method on various benchmarks ranging from small molecule datasets to complex periodic molecular systems with more types of elements. The experimental results show that the proposed pre-training method can greatly enhance the accuracy of MLIPs with few extra computational costs and works well with different invariant or equivariant graph neural network architectures. Our method improves the generalization capability of MLIPs and helps to realize accurate MD simulations for complex molecular systems.

Machine learning interatomic potentials (MLIPs) enable more efficient molecular dynamics (MD) simulations with ab initio accuracy, which have been used in various domains of physical science. However, distribution shift between training and test data causes deterioration of the test performance of MLIPs, and even leads to collapse of MD simulations. In this work, we propose an online Test-time Adaptation Interatomic Potential (TAIP) framework to improve the generalization on test data. Specifically, we design a dual-level self-supervised learning approach that leverages global structure and atomic local environment information to align the model with the test data. Extensive experiments demonstrate TAIP’s capability to bridge the domain gap between training and test dataset without additional data. TAIP enhances the test performance on various benchmarks, from small molecule datasets to complex periodic molecular systems with various types of elements. TAIP also enables stable MD simulations where the corresponding baseline models collapse.

Machine learning interatomic potential (MLIP) is an emerging technique that has helped achieve molecular dynamics simulations with unprecedented balance between efficiency and accuracy. Recently, the body of MLIP literature has been growing rapidly, which propels the need to automatically process relevant information for researchers to understand and utilize. Named entity recognition (NER), a natural language processing technique that identifies and categorizes information from texts, may help summarize key approaches and findings of relevant papers. In this work, we develop an NER model for MLIP literature by fine‐tuning a pre‐trained language model. To streamline text annotation, we build a user‐friendly web application for annotation and proofreading, which is seamlessly integrated into the training procedure. Our model can identify technical entities with an F1 score of 0.8 for new MLIP paper abstracts using only 60 training paper abstracts and up to 0.75 for scientific texts on different topics. Notably, some “errors” in predictions are actually reasonable decisions, showcasing the model's ability beyond what the performance metrics indicate. This work demonstrates the linguistic capabilities of the NER approach in processing textual information of a specific scientific domain and has the potential to accelerate materials research using language models and contribute to a user‐centric workflow.

A combination of machine learned interatomic potentials (MLIPs) and enhanced sampling simulations is used to investigate the activation of methane on a Ni(111) surface. The work entails the development and iterative refinement of MLIPs, initially trained on a data set constructed via ab initio molecular dynamics simulations, supplemented by adaptive biasing forces, to enrich the sampling of catalytically relevant configurations. Our results reveal that upon incorporation of collective variables that capture the behavior of the reactant molecule, as well as additional frames that describe the dynamic response of the catalytic surface, it is possible to enhance considerably the accuracy of predicted energies and forces. By employing enhanced sampling schemes in the refinement of the MLIP, we systematically explore the potential energy surface, leading to a refined MLIP capable of predicting density functional theory-level energies and forces and replicating key geometric characteristics of the catalytic system. The resulting free energy landscapes at several temperatures provide a detailed view of the thermodynamics and dynamics of methane activation. Specifically, as methane approaches and dissociates on the catalytic surface, the process involves the dynamic interplay of CH4 and the Ni catalyst that includes both enthalpic and entropic contributions. The progression toward the transition state involves a CH4 moiety that is increasingly restrained in its ability to rotate or translate, while the stage following the transition state is characterized by a notable rise of the Ni atom that interacts with the cleaved C-H bond. This leads to an increase in the mobility of the adsorbed species, a feature that becomes more pronounced at higher temperatures.

Using machine learning (ML) to construct interatomic interactions and thus potential energy surface (PES) has become a common strategy for materials design and simulations. However, those current models of machine-learning interatomic potential (MLIP) consider no relevant physical constraints or global scaling and thus may owe intrinsic out-of-domain difficulty which underlies the challenges of model generalizability and physical scalability. Here, by incorporating the global universal scaling law, we develop an ultrasmall parameterized MLIP with superlinear expressive capability, named SUS2-MLIP. Due to the global scaling derived from the universal equation of state (UEOS), SUS2-MLIP not only has significantly reduced parameters by decoupling the element space from coordinate space but also naturally outcomes the out-of-domain difficulty and endows the model with inherent generalizability and scalability even with relatively small training dataset. The non-linearity-embedding transformation in radial function endows the model with superlinear expressive capability. SUS2-MLIP outperforms the state-of-the-art MLIP models with its exceptional computational efficiency, especially for multiple-element materials and physical scalability in property prediction. This work not only presents a highly efficient universal MLIP model but also sheds light on incorporating physical constraints into AI-aided materials simulation.

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.

The chemical potential (μ) of species in solution is essential for understanding various chemical processes at interfaces. Molecular dynamics (MD) simulations, constrained by fixed compositions, cannot maintain constant chemical potential with reference to a targeted concentration or chemical potential under nonequilibrium or dynamic conditions, as solute species can migrate to the interface and deplete (or enrich) the bulk due to solute-interface interactions. In this study, we introduce a simple and computationally efficient approach named iterative quasi-constant chemical potential molecular dynamics (iqCμMD) simulation, which helps simulate targeted molar concentrations of species in solution. iqCμMD overcomes the limitations of conventional MD by adjusting the number of species in the solution to reach a target bulk concentration (chemical potential), which allows simulation of the interface under the bulk conditions comparable to experiment. We demonstrate our approach using machine learning interatomic potential (MLIP)-based MD simulations of the Na2SO4,aq-graphene interface, and to show the transferability of our approach, we also perform classical force field-based MD simulations of NaClaq-air and NaClaq-graphite interfaces, which produce comparable results to previous CμMD simulations. Our results also show that the iqCμMD approach efficiently achieves the desired bulk ion concentration within two iterations, and by utilizing MLIPs, we can achieve converged results using relatively small-scale simulations compared to previous CμMD simulations. By combining iqCμMD with MLIP-driven simulations, solid-liquid interfaces can be modeled under an effective constant chemical potential with DFT-level accuracy. We show that iqCμMD offers a robust and simple computational framework for constant chemical potential simulations, as its only requirement is to be able to converge interfacial simulations with a measurable bulk region.