Neural Network-Based Interatomic Potential for the Study of Thermal and Mechanical Properties of Siliceous Zeolites

François-Xavier Coudert,Luca Brugnoli,Maxime Ducamp
DOI: https://doi.org/10.26434/chemrxiv-2024-np8wf
2024-06-10
Abstract:Molecular dynamics simulations of zeolites are commonly employed for the characterization of their framework dynamics and response to the application of temperature and pressure. While classical interatomic potentials are commonly used for this task, they offer a description of the interactions in the system with limited accuracy. Density Functional Theory, meanwhile, is accurate but its high computational expense limits its scalability for large systems or long dynamics. Recent advances in machine learning interatomic potentials, trained on computational data obtained at the quantum chemical level, offer a promising alternative combining high accuracy with computational efficiency. In this study, we developed an MLIP specifically for pure silica zeolites, trained on data from high- temperature ab initio MD simulations across various zeolitic topologies. This MLIP was then applied to predict structural properties, thermal expansion, and pressure response of different zeolites, demonstrating its potential for accurate and generalizable in simulations of topologies beyond its initial training set.
Chemistry
What problem does this paper attempt to address?
This paper explores how to use machine learning (ML) to develop a neural network-based interatomic potential (MLIP) for studying the thermodynamic and mechanical properties of silicate zeolites. Traditional molecular dynamics simulations typically use classical interatomic potentials, which have limited accuracy, while density functional theory (DFT) is accurate but computationally expensive. The researchers obtained data from ab initio molecular dynamics simulations at high temperatures to train an MLIP specifically for pure silica zeolites. The trained MLIP was then used to predict different zeolite structures, thermal expansion, and pressure response, demonstrating its generalization ability beyond the initial training set. Specifically, the researchers first used DFT to compute a set of training data, and then trained the MLIP using NequIP, a method based on equivariant graph neural networks, optimizing hyperparameters to improve the accuracy and efficiency of the model. Once training was completed, the MLIP was tested on unseen zeolite structures, and the results showed good performance in predicting lattice parameters, volume, and thermodynamic properties, with higher accuracy compared to DFT calculations, and the ability to effectively simulate zeolites with different topologies. This indicates that MLIP is an efficient and accurate tool for studying zeolite materials.