Single-site Lennard-Jones models via polynomial chaos surrogate of Monte Carlo molecular simulation

by A. Kadoura, A. Sripatana, O. Knio, S. Sun, I. Hoteit
Year: 2016

Bibliography

Single-site Lennard-Jones models via polynomial chaos surrogate of Monte Carlo molecular simulation
A. Kadoura, A. Sripatana, O. Knio, S. Sun, and I. Hoteit
Journal of Chemical Physics, 144, 214301, 2016

Abstract

​In this work, two Polynomial Chaos (PC) surrogates were generated to reproduce Monte Carlo (MC) molecular simulation results of the canonical (single-phase) and the NVT-Gibbs (two-phase) ensembles for a system of normalized structureless Lennard-Jones (LJ) particles. The main advantage of such surrogates, once generated, is the capability of accurately computing the needed thermodynamic quantities in a few seconds, thus efficiently replacing the computationally expensive MC molecular simulations. Benefiting from the tremendous computational time reduction, the PC surrogates were used to conduct large-scale optimization in order to propose single-site LJ models for several simple molecules. Experimental data, a set of supercritical isotherms, and part of the two-phase envelope, of several pure components were used for tuning the LJ parameters (ϵ, σ). Based on the conducted optimization, excellent fit was obtained for different noble gases (Ar, Kr, and Xe) and other small molecules (CH4, N2, and CO). On the other hand, due to the simplicity of the LJ model used, dramatic deviations between simulation and experimental data were observed, especially in the two-phase region, for more complex molecules such as CO2 and C2 H6.

DOI: 10.1063/1.4952976

Keywords

Complex Molecules Computational Time Large-scale Optimization Lennard-Jones Model Molecular Simulations Polynomial Chaos Supercritical Isotherm