Physically Interpretable Interatomic Potentials via Symbolic Regression and Reinforcement Learning

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Varughese, Bilvin, Loeffler, Troy D., Banik, Suvo, Koneru, Aditya, Manna, Sukriti, Balasubramanian, Karthik, Batra, Rohit, Cherukara, Mathew J., Yildiz, Orcun, Peterka, Tom, Sumpter, Bobby G., Sankaranarayanan, Subramanian K. R. S.
Format: Preprint
Published: 2025
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909924716445696
author Varughese, Bilvin
Loeffler, Troy D.
Banik, Suvo
Koneru, Aditya
Manna, Sukriti
Balasubramanian, Karthik
Batra, Rohit
Cherukara, Mathew J.
Yildiz, Orcun
Peterka, Tom
Sumpter, Bobby G.
Sankaranarayanan, Subramanian K. R. S.
author_facet Varughese, Bilvin
Loeffler, Troy D.
Banik, Suvo
Koneru, Aditya
Manna, Sukriti
Balasubramanian, Karthik
Batra, Rohit
Cherukara, Mathew J.
Yildiz, Orcun
Peterka, Tom
Sumpter, Bobby G.
Sankaranarayanan, Subramanian K. R. S.
contents The development of next-generation molecular simulation models requires moving beyond pre-defined functional forms toward machine learning (ML) techniques that directly capture multiscale physics. Here, we demonstrate such an approach using symbolic regression (SR) with equation learner networks and a reinforcement learning search engine to derive interpretable equations for interatomic interactions. Training data were generated through nested ensemble sampling with density functional theory (DFT) energetics, spanning crystalline to highly disordered states. The optimization of the learner network employed continuous-action Monte Carlo Tree Search (MCTS) combined with gradient descent, enabling efficient exploration of function space. For copper as a representative transition metal, an unconstrained search produced models that outperformed fixed-form Sutton-Chen EAM potentials. The SR-derived models (SR1 and SR2) reproduced key material properties - lattice constants, cohesive energies, equations of state, elastic constants, phonon dispersion, defect formation energies, surface/bulk energetics, and phase transformation with significantly improved accuracy. Furthermore, stringent melting simulations using two-phase solid-amorphous interfaces confirmed that SR models accurately capture the interplay of vibrational entropy, cohesive energy, and structural dynamics, surpassing SC-EAM in both qualitative and quantitative predictions. This highlights the potential of SR to deliver fast, accurate, flexible, and physically meaningful potentials, advancing predictive modeling across scales.
format Preprint
id arxiv_https___arxiv_org_abs_2511_20506
institution arXiv
publishDate 2025
record_format arxiv
spellingShingle Physically Interpretable Interatomic Potentials via Symbolic Regression and Reinforcement Learning
Varughese, Bilvin
Loeffler, Troy D.
Banik, Suvo
Koneru, Aditya
Manna, Sukriti
Balasubramanian, Karthik
Batra, Rohit
Cherukara, Mathew J.
Yildiz, Orcun
Peterka, Tom
Sumpter, Bobby G.
Sankaranarayanan, Subramanian K. R. S.
Materials Science
The development of next-generation molecular simulation models requires moving beyond pre-defined functional forms toward machine learning (ML) techniques that directly capture multiscale physics. Here, we demonstrate such an approach using symbolic regression (SR) with equation learner networks and a reinforcement learning search engine to derive interpretable equations for interatomic interactions. Training data were generated through nested ensemble sampling with density functional theory (DFT) energetics, spanning crystalline to highly disordered states. The optimization of the learner network employed continuous-action Monte Carlo Tree Search (MCTS) combined with gradient descent, enabling efficient exploration of function space. For copper as a representative transition metal, an unconstrained search produced models that outperformed fixed-form Sutton-Chen EAM potentials. The SR-derived models (SR1 and SR2) reproduced key material properties - lattice constants, cohesive energies, equations of state, elastic constants, phonon dispersion, defect formation energies, surface/bulk energetics, and phase transformation with significantly improved accuracy. Furthermore, stringent melting simulations using two-phase solid-amorphous interfaces confirmed that SR models accurately capture the interplay of vibrational entropy, cohesive energy, and structural dynamics, surpassing SC-EAM in both qualitative and quantitative predictions. This highlights the potential of SR to deliver fast, accurate, flexible, and physically meaningful potentials, advancing predictive modeling across scales.
title Physically Interpretable Interatomic Potentials via Symbolic Regression and Reinforcement Learning
topic Materials Science
url https://arxiv.org/abs/2511.20506