Scalability of Coarse-Grained Potentials Generated from Iterative Boltzmann Inversion for Polymers: Case Study on Polycarbonates

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.

Energy renormalization for coarse-graining polymers having different segmental structures

Multiscale coarse-grained (CG) modeling of soft materials, such as polymers, is currently an art form because CG models normally have significantly altered dynamics and thermodynamic properties compared to their atomistic counterparts. We address this problem by exploiting concepts derived from the generalized entropy theory (GET), emphasizing the central role of configurational entropy s c in the dynamics of complex fluids. Our energy renormalization (ER) method involves varying the cohesive interaction strength in the CG models in such a way that dynamic properties related to s c are preserved. We test this ER method by applying it to coarse-graining polymer melts (i.e., polybutadiene, polystyrene, and polycarbonate), representing polymer materials having a relatively low, intermediate, and high degree of glass "fragility". We find that the ER method allows the dynamics of the atomistic polymer models to be faithfully described to a good approximation by CG models over a wide temperature range.

One-bead coarse-grained model for RNA dynamics

We present a revised version of a coarse-grained model for RNA dynamics. In such approach, the description of nucleotides is reduced to single points that interact between them through a series of effective pair potentials that were obtained from an improved analysis of RNA structures from the Protein Data Bank. These interaction potentials are the main constituents of a Brownian dynamics simulation algorithm that allows to perform a variety of tasks by taking advantage of the reduced number of variables. Such tasks include the prediction of the three-dimensional configuration of a series of test molecules. Moreover, the model permits the inclusion of effective magnesium ions and the ends of the RNA chains can be pulled with an external force to study the process of unfolding. In spite of the simplicity of the model, we obtain a good agreement with the experimental results.

Coarse-Grained Molecular Dynamics Force-Field for Polyacrylamide in Infinite Dilution Derived from Iterative Boltzmann Inversion and MARTINI Force-Field

We present a mesoscale model of aqueous polyacrylamide in the infinitely dilute concentration regime, by combining an extant coarse-grained (CG) force-field, MARTINI, and the Iterative Boltzmann Inversion protocol (IBI). MARTINI force-field was used to retain the thermodynamics of solvation of the polymer in water, whereas the structural properties and intrapolymer interactions were optimized by IBI. Atomistic molecular dynamics simulations of polymer in water were performed to benchmark the mesoscale simulations. Our results from the CG model show excellent agreement in structure with the atomistic system. We also studied the dynamical behavior of our CG system by computing the shear viscosity and compared it with the standard IBI model. The viscosity trends of our model were similar to the atomistic system, whereas the standard IBI model was highly dissimilar as expected. In summary, our hybrid CG model sufficiently mimics an infinitely dilute system, and is superior to both MARTINI and IBI in representing the structure and thermodynamics of the atomistic system, respectively. Our hybrid coarse-graining strategy promises applicability in large-scale simulations of polymeric/biological systems where the structure needs to be replicated accurately while preserving the thermodynamics of a smoother surrounding.