Abstract:A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

Ab initio machine learning simulation of calcium carbonate from aqueous solutions to the solid state

Machine learning potentials for complex aqueous systems made simple

An ab-initio deep neural network potential for accurate large-scale simulations of the muscovite mica-water interface

Molecular Dynamics Simulation on Aqueous Solution of Calcium Carbonate System

In-situ liquid SAXS studies on the early stage of calcium carbonate formation

Liquid to Crystal Si Growth Simulation Using Machine Learning Force Field

Machine Learning a Simple Interpretable Short-Range Potential for Silica

Phase Equilibrium of Water with Hexagonal and Cubic Ice Using the SCAN Functional

The nano- and meso-scale structure of amorphous calcium carbonate

Unravelling the dissolution dynamics of silicate minerals by deep learning molecular dynamics simulation: A case of dicalcium silicate

A deep learning potential applied in tobermorite phases and extended to calcium silicate hydrates

A deep learning-based potential developed for calcium silicate hydrates with both high accuracy and efficiency

Recalibrating the calcium trap in amino acid carboxyl groups via classical molecular dynamics simulations

Weaving the next-level structure of calcium silicate hydrate at the submicron scale via a remapping algorithm from coarse-grained to all-atom model

Efficient Machine Learning Force Field for Large-Scale Molecular Simulations of Organic Systems

Machine Learning in QM/MM Molecular Dynamics Simulations of Condensed-Phase Systems

A Coarse-Grained SPICA Makeover for Solvated and Bare Sodium and Chloride Ions

Machine learning for predictive condensed-phase simulation

Simulating the binding of key organic functional groups to aqueous calcium carbonate species

A first-principles machine-learning force field for heterogeneous ice nucleation on microcline feldspar