A machine learning closure for polymer integral equation theory

Polymer reference interaction site model (PRISM) theory, a descendant of Ornstein-Zernike liquid state theory, is a powerful tool to predict the structure and thermodynamics of equilibrium polymer systems, but its accuracy and applicability can be limited in some important cases. Typically, these shortcomings are traced to the analytical closure relationships used to solve the integral equations. Here, we propose a machine learning (ML)-based closure relation trained on a dataset of coarse-grained molecular dynamics simulations of homopolymer melts and solutions. PRISM theory with the ML closure outperforms traditional atomic closures (e.g., Percus-Yevick) in predicting the structure of typical coarse-grained model systems. We also use the ML closure to accurately model the results of small-angle neutron scattering experiments. This ML-enhanced PRISM theory can, therefore, enable rapid soft materials discovery and design efforts.