Abstract:Key to being able to accurately model the properties of realistic materials is being able to predict their properties in the thermodynamic limit. Nevertheless, because most many-body electronic structure methods scale as a high-order polynomial, or even exponentially, with system size, directly simulating large systems in their thermodynamic limit rapidly becomes computationally intractable. As a result, researchers typically estimate the properties of large systems that approach the thermodynamic limit by extrapolating the properties of smaller, computationally-accessible systems based on relatively simple scaling expressions. In this work, we employ Gaussian processes to more accurately and efficiently extrapolate many-body simulations to their thermodynamic limit. We train our Gaussian processes on Smooth Overlap of Atomic Positions (SOAP) descriptors to extrapolate the energies of one-dimensional hydrogen chains obtained using two, high-accuracy many-body methods: Coupled Cluster theory and Auxiliary Field Quantum Monte Carlo (AFQMC). In so doing, we show that Gaussian processes trained on relatively short, 10-30-atom chains can predict the energies of both homogeneous and inhomogeneous hydrogen chains in their thermodynamic limit with sub-milliHartree accuracy. Unlike standard scaling expressions, our GPR-based approach is highly generalizable given representative training data and is not dependent on systems' geometries or dimensionality. This work highlights the potential for machine learning to correct for the finite size effects that routinely complicate the interpretation of finite size many-body simulations.

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

Accelerated Bayesian Inference for Molecular Simulations using Local Gaussian Process Surrogate Models

An Adaptive Gaussian Process‐based Method for Efficient Bayesian Experimental Design in Groundwater Contaminant Source Identification Problems

Fast, accurate and lightweight sequential simulation-based inference using Gaussian locally linear mappings

Gaussian Processes for Surrogate Modeling of Discharged Fuel Nuclide Compositions

Iterative Construction of Gaussian Process Surrogate Models for Bayesian Inference

A Bayesian approach for parameter estimation and prediction using a computationally intensive model

Fast Bayesian Force Fields from Active Learning: Study of Inter-Dimensional Transformation of Stanene

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Additive Multi-Index Gaussian process modeling, with application to multi-physics surrogate modeling of the quark-gluon plasma

Accelerating a hybrid continuum-atomistic fluidic model with on-the-fly machine learning

Gaussian Processes for Finite Size Extrapolation of Many-Body Simulations

Accelerating inference for stochastic kinetic models

Parameter optimisation using Bayesian inference for spallation models

Langevin dynamics simulations of biomolecules on graphics processors

NNP/MM: Accelerating molecular dynamics simulations with machine learning potentials and molecular mechanic

An adaptive Gaussian process method for multi-modal Bayesian inverse problems

Computing linear optical spectra in the presence of nonadiabatic effects on Graphics Processing Units using molecular dynamics and tensor-network approaches