Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation

Silica precipitation is a subject of big interest since it occurs in a wide variety of environmental and industrial processes. Even though there are many advances in atomistic simulation research of different forms of silica, the mechanism of silica precipitation has not been fully understood. We pr...

Full description

Bibliographic Details
Main Author: Orlova, Yuliia
Other Authors: Gómez-Bombarelli, Rafael
Format: Thesis
Published: Massachusetts Institute of Technology 2024
Online Access:https://hdl.handle.net/1721.1/154194
_version_ 1826193445984665600
author Orlova, Yuliia
author2 Gómez-Bombarelli, Rafael
author_facet Gómez-Bombarelli, Rafael
Orlova, Yuliia
author_sort Orlova, Yuliia
collection MIT
description Silica precipitation is a subject of big interest since it occurs in a wide variety of environmental and industrial processes. Even though there are many advances in atomistic simulation research of different forms of silica, the mechanism of silica precipitation has not been fully understood. We propose to study the following process using reactive force-field method (ReaxFF). Despite being a classical force field, ReaxFF can achieve quantum chemical accuracy once the optimal potential coefficients are found. However, the fitting of ReaxFF parameters is a challenge due to the complex functional form of the potential. Several techniques have been proposed to solve this problem, such as evolutionary algorithms, Monte Carlo methods, and simulated annealing. The stochastic nature of these methods requires millions of error evaluations to fit the parameters, which results in excessive optimization times. Recent advances in machine learning made it possible to drastically speed up the process by utilizing the gradient of the potential. In this work, the gradient-based optimization of reactive force-field parameters using Pytorch was performed. We have implemented ReaxFF potential as a Pytorch model. The model’s performance was validated against existing ReaxFF implementations. ReaxFF parameters were fitted to the dataset, which comprised 15345 geometries calculated using a long-range corrected hybrid functional 𝜔B97XD3.
first_indexed 2024-09-23T09:39:21Z
format Thesis
id mit-1721.1/154194
institution Massachusetts Institute of Technology
last_indexed 2024-09-23T09:39:21Z
publishDate 2024
publisher Massachusetts Institute of Technology
record_format dspace
spelling mit-1721.1/1541942024-04-18T03:55:25Z Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation Orlova, Yuliia Gómez-Bombarelli, Rafael Massachusetts Institute of Technology. Center for Computational Science and Engineering Silica precipitation is a subject of big interest since it occurs in a wide variety of environmental and industrial processes. Even though there are many advances in atomistic simulation research of different forms of silica, the mechanism of silica precipitation has not been fully understood. We propose to study the following process using reactive force-field method (ReaxFF). Despite being a classical force field, ReaxFF can achieve quantum chemical accuracy once the optimal potential coefficients are found. However, the fitting of ReaxFF parameters is a challenge due to the complex functional form of the potential. Several techniques have been proposed to solve this problem, such as evolutionary algorithms, Monte Carlo methods, and simulated annealing. The stochastic nature of these methods requires millions of error evaluations to fit the parameters, which results in excessive optimization times. Recent advances in machine learning made it possible to drastically speed up the process by utilizing the gradient of the potential. In this work, the gradient-based optimization of reactive force-field parameters using Pytorch was performed. We have implemented ReaxFF potential as a Pytorch model. The model’s performance was validated against existing ReaxFF implementations. ReaxFF parameters were fitted to the dataset, which comprised 15345 geometries calculated using a long-range corrected hybrid functional 𝜔B97XD3. S.M. 2024-04-17T21:10:34Z 2024-04-17T21:10:34Z 2023-09 2023-10-11T18:01:51.643Z Thesis https://hdl.handle.net/1721.1/154194 In Copyright - Educational Use Permitted Copyright MIT http://rightsstatements.org/page/InC-EDU/1.0/ application/pdf Massachusetts Institute of Technology
spellingShingle Orlova, Yuliia
Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title_full Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title_fullStr Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title_full_unstemmed Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title_short Gradient-Based Optimization of ReaxFF Parameters Using Pytorch for the Study of Silica Precipitation
title_sort gradient based optimization of reaxff parameters using pytorch for the study of silica precipitation
url https://hdl.handle.net/1721.1/154194
work_keys_str_mv AT orlovayuliia gradientbasedoptimizationofreaxffparametersusingpytorchforthestudyofsilicaprecipitation