Transferability of Local Density-Assisted Implicit Solvation Models for Homogeneous Fluid Mixtures

Insight into the Density-Dependence of Pair Potentials for Predictive Coarse-Grained Models

We investigate the temperature- and density-dependence of effective pair potentials for 1-site coarse-grained (CG) models of two industrial solvents, 1,4-dioxane and tetrahydrofuran. We observe that the calculated pair potentials are much more sensitive to density than to temperature. The generalized-Yvon-Born-Green framework reveals that this striking density-dependence reflects corresponding variations in the many-body correlations that determine the environment-mediated indirect contribution to the pair mean force. Moreover, we demonstrate, perhaps surprisingly, that this density-dependence is not important for accurately modeling the intermolecular structure. Accordingly, we adopt a density-independent interaction potential and transfer the density-dependence of the calculated pair potentials into a configuration-independent volume potential. Furthermore, we develop a single global potential that accurately models the intermolecular structure and pressure-volume equation of state across a very wide range of liquid state points. Consequently, this work provides fundamental insight into the density-dependence of effective pair potentials and also provides a significant step toward developing predictive CG models for efficiently modeling industrial solvents.

Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability

Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.

Stability, Speed, and Constraints for Structural Coarse-Graining in VOTCA

Structural coarse-graining involves the inverse problem of deriving pair potentials that reproduce target radial distribution functions. Despite its clear mathematical formulation, there are open questions about the existing methods concerning speed, stability, and physical representability of the resulting potentials. In this work, we make progress on several aspects of iterative methods used to solve the inverse problem. Based on integral equation theory, we derive fast Gauss-Newton schemes applicable to very general systems, including molecules with bonds and mixtures. Our methods are similar to inverse Monte Carlo in terms of convergence speed and have a similar cost per iteration as iterative Boltzmann inversion. We investigate stability problems in our schemes and in the inverse Monte Carlo method and propose modifications to fix them. Furthermore, we establish how the pair potential can be constrained at each iteration to reproduce the pressure, Kirkwood-Buff integral, or the enthalpy of vaporization. We demonstrate the potential of our approach in deriving coarse-grained force fields for nine different solvents and their mixtures. All methods described are implemented in the free and open VOTCA software framework for systematic coarse-graining.

Consistent and Transferable Force Fields for Statistical Copolymer Systems at the Mesoscale

The statistical trajectory matching (STM) method was applied successfully to derive coarse grain (CG) models for bulk properties of homopolymers. The extension of the methodology for building CG models for statistical copolymer systems is much more challenging. We present here the strategy for developing CG models for styrene-butadiene-rubber, and we compare the quality of the resulting CG force fields on the structure and thermodynamics at different chemical compositions. The CG models are used through the use of a genuine mesoscopic method called the dissipative particle dynamics method and compared to high-resolution molecular dynamics simulations. We conclude that the STM method is able to produce coarse-grained potentials that are transferable in composition by using only a few reference systems. Additionally, this methodology can be applied on any copolymer system.

Broad chemical transferability in structure-based coarse-graining

Compared to top-down coarse-grained (CG) models, bottom-up approaches are capable of offering higher structural fidelity. This fidelity results from the tight link to a higher-resolution reference, making the CG model chemically specific. Unfortunately, chemical specificity can be at odds with compound-screening strategies, which call for transferable parametrizations. Here we present an approach to reconcile bottom-up, structure-preserving CG models with chemical transferability. We consider the bottom-up CG parametrization of 3,441 C7O2 small-molecule isomers. Our approach combines atomic representations, unsupervised learning, and a large-scale extended-ensemble force-matching parametrization. We first identify a subset of 19 representative molecules, which maximally encode the local environment of all gas-phase conformers. Reference interactions between the 19 representative molecules were obtained from both homogeneous bulk liquids and various binary mixtures. An extended-ensemble parametrization over all 703 state points leads to a CG model that is both structure-based and chemically transferable. Remarkably, the resulting force field is on average more structurally accurate than single-state-point equivalents. Averaging over the extended ensemble acts as a mean-force regularizer, smoothing out both force and structural correlations that are overly specific to a single state point. Our approach aims at transferability through a set of CG bead types that can be used to easily construct new molecules, while retaining the benefits of a structure-based parametrization.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

Investigating the energetic and entropic components of effective potentials across a glass transition

By eliminating unnecessary details, coarse-grained (CG) models provide the necessary efficiency for simulating scales that are inaccessible to higher resolution models. However, because they average over atomic details, the effective potentials governing CG degrees of freedom necessarily incorporate significant entropic contributions, which limit their transferability and complicate the treatment of thermodynamic properties. This work employs a dual-potential approach to consider the energetic and entropic contributions to effective interaction potentials for CG models. Specifically, we consider one- and three-site CG models for ortho-terphenyl (OTP) both above and below its glass transition. We employ the multiscale coarse-graining (MS-CG) variational principle to determine interaction potentials that accurately reproduce the structural properties of an all-atom (AA) model for OTP at each state point. We employ an energy-matching variational principle to determine an energy operator that accurately reproduces the intra- and inter-molecular energy of the AA model. While the MS-CG pair potentials are almost purely repulsive, the corresponding pair energy functions feature a pronounced minima that corresponds to contacting benzene rings. These energetic functions then determine an estimate for the entropic component of the MS-CG interaction potentials. These entropic functions accurately predict the MS-CG pair potentials across a wide range of liquid state points at constant density. Moreover, the entropic functions also predict pair potentials that quite accurately model the AA pair structure below the glass transition. Thus, the dual-potential approach appears a promising approach for modeling AA energetics, as well as for predicting the temperature-dependence of CG effective potentials.

A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). I. General theory and model

Water is undoubtedly one of the most important molecules for a variety of chemical and physical systems, and constructing precise yet effective coarse-grained (CG) water models has been a high priority for computer simulations. To recapitulate important local correlations in the CG water model, explicit higher-order interactions are often included. However, the advantages of coarse-graining may then be offset by the larger computational cost in the model parameterization and simulation execution. To leverage both the computational efficiency of the CG simulation and the inclusion of higher-order interactions, we propose a new statistical mechanical theory that effectively projects many-body interactions onto pairwise basis sets. The many-body projection theory presented in this work shares similar physics from liquid state theory, providing an efficient approach to account for higher-order interactions within the reduced model. We apply this theory to project the widely used Stillinger-Weber three-body interaction onto a pairwise (two-body) interaction for water. Based on the projected interaction with the correct long-range behavior, we denote the new CG water model as the Bottom-Up Many-Body Projected Water (BUMPer) model, where the resultant CG interaction corresponds to a prior model, the iteratively force-matched model. Unlike other pairwise CG models, BUMPer provides high-fidelity recapitulation of pair correlation functions and three-body distributions, as well as N-body correlation functions. BUMPer extensively improves upon the existing bottom-up CG water models by extending the accuracy and applicability of such models while maintaining a reduced computational cost.

Temperature and Phase Transferable Bottom-up Coarse-Grained Models

Despite the high fidelity of bottom-up coarse-grained (CG) approaches to recapitulate the structural correlations in atomistic simulations, the general use of bottom-up CG methods is limited because of the nontransferable nature of these CG models under different thermodynamic conditions. Because bottom-up CG potentials usually correspond to configuration-dependent free energies of the system, recent studies have focused on adjusting enthalpic or entropic contributions to account for issues with transferability. However, these approaches can require a manual adjustment of the CG interaction a priori and are usually limited to constant volume ensembles. To overcome these limitations, we construct temperature and phase transferable CG models under constant pressure by developing the ultra-coarse-graining (UCG) methodology in the mean-field limit. In the mean-field ansatz, an embedded semi-global order parameter recapitulates global changes to the system by automatically adjusting the effective CG interactions, thus bridging free energy decompositions with UCG theory. The method presented is designed to faithfully capture structural correlations under different thermodynamic conditions, using a single UCG model. Specifically, we test the applicability of the developed theory in three distinct cases: (1) different temperatures at constant pressure in liquids, (2) different temperatures across thermodynamic phases, and (3) liquid/vapor interfaces. We demonstrate that the systematic construction of both temperature and phase transferable bottom-up CG models is possible using this generalized UCG theory. Based on our findings, this approach significantly extends the transferability and applicability of the bottom-up CG theory and method.

Coarse-graining of polyisoprene melts using inverse Monte Carlo and local density potentials

Bottom-up coarse-graining of polymers is commonly performed by matching structural order parameters such as distribution of bond lengths, bending and dihedral angles, and pair distribution functions. In this study, we introduce the distribution of nearest-neighbors as an additional order parameter in the concept of local density potentials. We describe how the inverse-Monte Carlo method provides a framework for forcefield development that is capable of overcoming challenges associated with the parameterization of interaction terms in polymer systems. The technique is applied on polyisoprene melts as a prototype system. We demonstrate that while different forcefields can be developed that perform equally in terms of matching target distributions, the inclusion of nearest-neighbors provides a straightforward route to match both thermodynamic and conformational properties. We find that several temperature state points can also be addressed, provided that the forcefield is refined accordingly. Finally, we examine both the single-particle and the collective dynamics of the coarse-grain models, demonstrating that all forcefields present a similar acceleration relative to the atomistic systems.

Analysis of local density potentials

Low resolution coarse-grained (CG) models are widely adopted for investigating phenomena that cannot be effectively simulated with all-atom (AA) models. Since the development of the many-body dissipative particle dynamics method, CG models have increasingly supplemented conventional pair potentials with one-body potentials of the local density (LD) around each site. These LD potentials appear to significantly extend the transferability of CG models, while also enabling more accurate descriptions of thermodynamic properties, interfacial phenomena, and many-body correlations. In this work, we systematically examine the properties of LD potentials. We first derive and numerically demonstrate a nontrivial transformation of pair and LD potentials that leaves the total forces and equilibrium distribution invariant. Consequently, the pair and LD potentials determined via bottom-up methods are not unique. We then investigate the sensitivity of CG models for glycerol to the weighting function employed for defining the local density. We employ the multiscale coarse-graining (MS-CG) method to simultaneously parameterize both pair and LD potentials. When employing a short-ranged Lucy function that defines the local density from the first solvation shell, the MS-CG model accurately reproduces the pair structure, pressure-density equation of state, and liquid-vapor interfacial profile of the AA model. The accuracy of the model generally decreases as the range of the Lucy function increases further. The MS-CG model provides similar accuracy when a smoothed Heaviside function is employed to define the local density from the first solvation shell. However, the model performs less well when this function acts on either longer or shorter length scales.