Search for: All records

Coarse-grained (CG) modeling has gained significant attention in recent years due to its wide applicability in enhancing the spatiotemporal scales of molecular simulations. While CG simulations, often performed with Hamiltonian mechanics, faithfully recapitulate structural correlations at equilibrium, they lead to ambiguously accelerated dynamics. In Paper I [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034103 (2023)], we proposed the excess entropy scaling relationship to understand the CG dynamics. Then, in Paper II [J. Jin, K. S. Schweizer, and G. A. Voth, J. Chem. Phys. 158(3), 034104 (2023)], we developed a theory to map the CG system into a dynamically consistent hard sphere system to analytically derive an expression for fast CG dynamics. However, many chemical and physical systems do not exhibit hard sphere-like behavior, limiting the extensibility of the developed theory. In this paper, we aim to generalize the theory to the non-hard sphere system based on the Weeks–Chandler–Andersen perturbation theory. Since non-hard sphere-like CG interactions affect the excess entropy term as it deviates from the hard sphere description, we explicitly account for the extra entropy to correct the non-hard sphere nature of the system. This approach is demonstrated for two different types of interactions seen in liquids, and we further provide a generalized description for any CG models using the generalized Gaussian CG models using Gaussian basis sets. Altogether, this work allows for extending the range and applicability of the hard sphere CG dynamics theory to a myriad of CG liquids. 

Bottom-up coarse-grained (CG) modeling is an effective means of bypassing the limited spatiotemporal scales of conventional atomistic molecular dynamics while retaining essential information from the atomistic model. A central challenge in CG modeling is the trade-off between accuracy and efficiency, as the inclusion of often pivotal many-body interaction terms in the CG force-field renders simulation markedly slower than simple pairwise models. The Ultra Coarse-Graining (UCG) method incorporates many-body terms through discrete internal state variables that modulate the CG force-field according to, e.g., changes in local environment when substantial chemical heterogeneities exist. However, assigning optimal internal states systematically from atomistic simulation data, as well as the practical application of bottom-up UCG theory to biomolecular systems, remain open problems. We develop two synergistic methods to aid in the development of UCG models that can capture inhomogeneities in atomistic systems such as those induced by phase coexistence. The first method establishes the systematic construction of UCG force-fields from a relative entropy minimization principle, while the second method utilizes machine-learning to obtain optimal local order parameters for enhanced model efficiency and transferability. We apply these methods to a methanol liquid–vapor interface and the ripple phase of a 1,2-dipalmitoyl-sn-glycero-3-phosphocholine lipid bilayer and demonstrate that UCG modeling alone recapitulates aspects of phase coexistence that are otherwise not observed in CG modeling. 

Applying an excess entropy scaling formalism to the coarse-grained (CG) dynamics of liquids, we discovered that missing rotational motions during the CG process are responsible for artificially accelerated CG dynamics. In the context of the dynamic representability between the fine-grained (FG) and CG dynamics, this work introduces the well-known Stokes–Einstein and Stokes–Einstein–Debye relations to unravel the rotational dynamics underlying FG trajectories, thereby allowing for an indirect evaluation of the effective rotations based only on the translational information at the reduced CG resolution. Since the representability issue in CG modeling limits a direct evaluation of the shear stress appearing in the Stokes–Einstein and Stokes–Einstein–Debye relations, we introduce a translational relaxation time as a proxy to employ these relations, and we demonstrate that these relations hold for the ambient conditions studied in our series of work. Additional theoretical links to our previous work are also established. First, we demonstrate that the effective hard sphere radius determined by the classical perturbation theory can approximate the complex hydrodynamic radius value reasonably well. Furthermore, we present a simple derivation of an excess entropy scaling relationship for viscosity by estimating the elliptical integral of molecules. In turn, since the translational and rotational motions at the FG level are correlated to each other, we conclude that the “entropy-free” CG diffusion only depends on the shape of the reference molecule. Our results and analyses impart an alternative way of recovering the FG diffusion from the CG description by coupling the translational and rotational motions at the hydrodynamic level. 

Coarse-grained (CG) molecular dynamics can be a powerful method for probing complex processes. However, most CG force fields use pairwise nonbonded interaction potentials sets, which can limit their ability to capture complex multi-body phenomena such as the hydrophobic effect. As the hydrophobic effect primarily manifests itself due to the nonpolar solute affecting the nearby hydrogen bonding network in water, capturing such effects using a simple one CG site or “bead” water model is a challenge. In this work, we systematically test the ability of CG one site water models for capturing critical features of the solvent environment around a hydrophobe as well as the potential of mean force (PMF) of neopentane association. We study two bottom-up models: a simple pairwise (SP) force-matched water model constructed using the multiscale coarse-graining method and the Bottom-Up Many-Body Projected Water (BUMPer) model, which has implicit three-body correlations. We also test the top-down monatomic (mW) and the Machine Learned mW (ML-mW) water models. The mW models perform well in capturing structural correlations but not the energetics of the PMF. BUMPer outperforms SP in capturing structural correlations and also gives an accurate PMF in contrast to the two mW models. Our study highlights the importance of including three-body interactions in CG water models, either explicitly or implicitly, while in general highlighting the applicability of bottom-up CG water models for studying hydrophobic effects in a quantitative fashion. This assertion comes with a caveat, however, regarding the accuracy of the enthalpy–entropy decomposition of the PMF of hydrophobe association. 

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. 

This paper series aims to establish a complete correspondence between fine-grained (FG) and coarse-grained (CG) dynamics by way of excess entropy scaling (introduced in Paper I). While Paper II successfully captured translational motions in CG systems using a hard sphere mapping, the absence of rotational motions in single-site CG models introduces differences between FG and CG dynamics. In this third paper, our objective is to faithfully recover atomistic diffusion coefficients from CG dynamics by incorporating rotational dynamics. By extracting FG rotational diffusion, we unravel, for the first time reported to our knowledge, a universality in excess entropy scaling between the rotational and translational diffusion. Once the missing rotational dynamics are integrated into the CG translational dynamics, an effective translation-rotation coupling becomes essential. We propose two different approaches for estimating this coupling parameter: the rough hard sphere theory with acentric factor (temperature-independent) or the rough Lennard-Jones model with CG attractions (temperature-dependent). Altogether, we demonstrate that FG diffusion coefficients can be recovered from CG diffusion coefficients by (1) incorporating “entropy-free” rotational diffusion with translation-rotation coupling and (2) recapturing the missing entropy. Our findings shed light on the fundamental relationship between FG and CG dynamics in molecular fluids.