Dynamic evolution in coarse-grained molecular dynamics simulations of polyethylene melts

Roughness Volumes: An Improved RoughMob Concept for Predicting the Increase of Molecular Mobility upon Coarse-Graining

The reduced number of degrees of freedom in a coarse-grained molecular model compared to its parent atomistic model not only makes it possible to simulate larger systems for longer time scales but also results in an artificial mobility increase. The RoughMob method [Meinel, M. K. and Müller-Plathe, F. J. Chem. Theory Comput. 2020, 16, 1411.] linked the acceleration factor of the dynamics to the loss of geometric information upon coarse-graining. Our hypothesis is that coarse-graining a multiatom molecule or group into a single spherical bead smooths the molecular surface and, thus, leads to reduced intermolecular friction. A key parameter is the molecular roughness difference, which is calculated via a numerical comparison of the molecular surfaces of both the atomistic and coarse-grained models. Augmenting the RoughMob method, we add the concept of the region where the roughness acts. This information is contained in four so-called roughness volumes. For 17 systems of homogeneous hydrocarbon fluids, simple one-bead coarse-grained models are derived by the structure-based iterative Boltzmann inversion. They include 13 different homogeneous aliphatic and aromatic molecules and two different mapping schemes. We present a simple way to correlate the roughness volumes to the acceleration factor. The resulting relation is able to a priori predict the acceleration factors for an extended size and shape range of hydrocarbon molecules, with different mapping schemes and different densities.

Coarse-grained molecular dynamics simulations of poly(ethylene terephthalate)

We have constructed efficient coarse-grained (CG) models of poly(ethylene terephthalate) (PET), using three mapping schemes, in which a repeat unit is lumped into either three or four beads. The CG potentials are parameterized to reproduce target distributions of an underlying accurate atomistic model [H. Eslami and F. Müller-Plathe, Macromolecules 42, 8241-8250 (2009)]. The CG simulations allow equilibration of long PET chains at all length scales. The CG results on the density of PET in melt and glassy states, chain dimension, local packing, and structure factor are in good agreement with experiment. We have established a link between the glass transition temperature and the local movements including conformational transitions and mean-square displacements of chain segments. Temperature transferabilities of the three proposed models were studied by comparing CG results on the static and thermodynamic properties of a polymer with atomistic and experimental findings. One of the three CG models has a good degree of transferability, following all inter- and intra-structural rearrangements of the atomistic model, over a broad range of temperature. Furthermore, as a distinct point of strength of CG, over atomistic, simulations, we have examined the dynamics of PET long chains, consisting of 100 repeat units, over a regime where entanglements dominate the dynamics. Performing long-time (550 ns) CG simulations, we have noticed the signature of a crossover from Rouse to reptation dynamics. However, a clear separation between the Rouse and the reptation dynamics needs much longer time simulations, confirming the experimental findings that the crossover to full reptation dynamics is very protracted.

Bottom-up approach to represent dynamic properties in coarse-grained molecular simulations

Several molecular coarse-graining methods have been proposed in recent years to derive chemical- and state-point transferable force fields. While these force fields describe structural and thermodynamic properties in good agreement with fine-grained models and experiments, dynamic properties are usually overestimated. Herein, we examine if the long-time dynamic properties of molecular coarse-grained (CG) systems can be correctly represented by employing a dissipative particle dynamics (DPD) thermostat, which is "bottom-up informed" by means of a variant of the Markovian Mori-Zwanzig (MZ) DPD coarse-graining method. We report single-site and multiple-site CG models for a monomer, dimer, and 24mer based on 2,2-dimethyl propane as a chemical repeat unit and report data obtained from MZ-DPD simulations of liquids, polymer solutions, and polymer melts. We find that despite incomplete time scale separation of the molecular CG model, MZ-DPD achieves quantitative accuracy in predicting diffusive dynamics in single-component liquids and polymer solutions (24mers in a dimer solvent). We also find that MZ-DPD simulations of molecular penetrant diffusion in polymer networks do not reach quantitative agreement with the fine-grained model. Modeling diffusion governed by the activated barrier crossing of small molecular penetrants in these dense systems requires an accurate description of energy barriers, presumably combined with the treatment of memory effects. The use of a MZ-DPD thermostat extends the scope and applicability of molecular CG models for multicomponent systems where a correct description of the relative diffusion rates of the different components is important.

Systematic derivation of hybrid coarse-grained models

Molecular dynamics represents a key enabling technology for applications ranging from biology to the development of new materials. However, many real-world applications remain inaccessible to fully resolved simulations due to their unsustainable computational costs and must therefore rely on semiempirical coarse-grained models. Significant efforts have been devoted in the last decade towards improving the predictivity of these coarse-grained models and providing a rigorous justification of their use, through a combination of theoretical studies and data-driven approaches. One of the most promising research efforts is the (re)discovery of the Mori-Zwanzig projection as a generic, yet systematic, theoretical tool for deriving coarse-grained models. Despite its clean mathematical formulation and generality, there are still many open questions about its applicability and assumptions. In this work, we propose a detailed derivation of a hybrid multiscale system, generalizing and further investigating the approach developed in Español [Europhys. Lett. 88, 40008 (2009)10.1209/0295-5075/88/40008]. Issues such as the general coexistence of atoms (fully resolved degrees of freedom) and beads (larger coarse-grained units), the role of the fine-to-coarse mapping chosen, and the approximation of effective potentials are discussed. The theoretical discussion is supported by numerical simulations of a monodimensional nonlinear periodic benchmark system with an open-source parallel Julia code, easily extensible to arbitrary potential models and fine-to-coarse mapping functions. The results presented highlight the importance of introducing, in the macroscopic model, nonconstant fluctuating and dissipative terms, given by the Mori-Zwanzig approach, to correctly reproduce the reference fine-grained results, without requiring ad hoc calibration of interaction potentials and thermostats.

Onset of static and dynamic universality among molecular models of polymers

A quantitatively accurate prediction of properties for entangled polymers is a long-standing challenge that must be addressed to enable efficient development of these materials. The complex nature of polymers is the fundamental origin of this challenge. Specifically, the chemistry, structure, and dynamics at the atomistic scale affect properties at the meso and macro scales. Therefore, quantitative predictions must start from atomistic molecular dynamics (AMD) simulations. Combined use of atomistic and coarse-grained (CG) models is a promising approach to estimate long-timescale behavior of entangled polymers. However, a systematic coarse-graining is still to be done for bridging the gap of length and time scales while retaining atomistic characteristics. Here we examine the gaps among models, using a generic mapping scheme based on power laws that are closely related to universality in polymer structure and dynamics. The scheme reveals the characteristic length and time for the onset of universality between the vastly different scales of an atomistic model of polyethylene and the bead-spring Kremer-Grest (KG) model. The mapping between CG model of polystyrene and the KG model demonstrates the fast onset of universality, and polymer dynamics up to the subsecond time scale are observed. Thus, quantitatively traceable timescales of polymer MD simulations can be significantly increased.

Directly Modifying the Nonbonded Potential Based on the Standard Iterative Boltzmann Inversion Method for Coarse-Grained Force Fields

Effective potentials are of great importance for coarse-grained (CG) simulations, which can be obtained by the structure-based iterative Boltzmann inversion (IBI) method. However, the standard IBI method is incapable of maintaining the mechanical and thermodynamic properties of the CG model in agreement with those of the all-atom model. Unlike the existing techniques, such as introducing friction force as the dissipative force to reduce the superatom motion while keeping the conservative force arising from the CG potential intact, we directly modified the standard IBI nonbonded potential by adding an empirical function. According to an analysis of the dissipative particle dynamics, the additional function did compensate for the friction reduction of the standard IBI CG model. In this work, the thermal fluctuation information from the nonbonded radial distribution function was incorporated into the additional empirical function. As an illustration of the new CG force fields, we presented simulations of the stress-strain relation and thermodynamic properties in terms of cis-polyisoprene and compared the statistical structure information of the superatoms with those of the IBI CG model and the all-atom model. It should be emphasized that the additional empirical function contributed to compensating for the friction reduction, irrespective of the functional form it took. In this sense, the proposed method was easily operable.

Multi-scale molecular dynamics study of cholera pentamer binding to a GM1-phospholipid membrane

The AB5 type toxin produced by the Vibrio cholerae bacterium is the causative agent of the cholera disease. The cholera toxin (CT) has been shown to bind specifically to GM1 glycolipids on the membrane surface. This binding of CT to the membrane is the initial step in its endocytosis and has been postulated to cause significant disruption to the membrane structure. In this work, we have carried out a combination of coarse-grain and atomistic simulations to study the binding of CT to a membrane modelled as an asymmetrical GM1-DPPC bilayer. Simulation results indicate that the toxin binds to the membrane through only three of its five B subunits, in effect resulting in a tilted bound configuration. Additionally, the binding of the CT can increase the area per lipid of GM1 leaflet, which in turn can cause the membrane regions interacting with the bound subunits to experience significant bilayer thinning and lipid tail disorder across both the leaflets.

Ion-mediated charge transport in ionomeric electrolytes

Ionomers, or single-ion conductors, serve as a model system to study ion transport in polymeric systems. Conductivity is a system property that depends on the net charge transport in the system. The mechanism through which ions are transported can dramatically change the contribution of an ion's self-motion (i.e. diffusion coefficient) to the conductivity of the system. For example, positive and negative ions diffusing as a pair have no net contribution to conductivity. In a coarse-grained molecular dynamics simulation of sodium-neutralized poly(PEO-co-sulfoisophthalate), we show that ion transport is mediated through consecutive coordination with ion pairs and higher order clusters due to the high density of ions. This transport mechanism is highly efficient and shows evidence of cation relaying. We show that larger ion aggregates can serve as ion-conducting paths for positive charges, and demonstrate how a highly ordered ion aggregate network can improve conductivity by enhancing correlated ion transport.

A multiscale approach to model hydrogen bonding: The case of polyamide

We present a simple multiscale model for polymer chains in which it is possible to selectively remove degrees of freedom. The model integrates all-atom and coarse-grained potentials in a simple and systematic way and allows a fast sampling of the complex conformational energy surface typical of polymers whilst maintaining a realistic description of selected atomistic interactions. In particular, we show that it is possible to simultaneously reproduce the structure of highly directional non-bonded interactions such as hydrogen bonds and efficiently explore the large number of conformations accessible to the polymer chain. We apply the method to a melt of polyamide removing from the model only the degrees of freedom associated to the aliphatic segments and keeping at atomistic resolution the amide groups involved in the formation of the hydrogen bonds. The results show that the multiscale model produces structural properties that are comparable with the fully atomistic model despite being five times faster to simulate.

Phase segregation in bio-inspired multi-component vesicles encompassing double tail phospholipid species

Our aim is to investigate the phase segregation and the structure of multi-component bio-inspired phospholipid vesicles via dissipative particle dynamics. The chemical distinction in the phospholipid species arises due to different head and tail group moieties, and molecular stiffness of the hydrocarbon tails. The individual amphiphilic phospholipid molecular species are represented by a hydrophilic head group and two hydrophobic tails. The distinct chemical nature of the moieties is modeled effectively via soft repulsive interaction parameters, and the molecular rigidity is tuned via suitable three-body potential constants. We demonstrate the formation of a stable hybrid vesicle through the self-assembly of the amphiphilic phospholipid molecules in the presence of a hydrophilic solvent. We investigate and characterize the phase segregation and the structure of the binary vesicles for different phospholipid mixtures. Our results demonstrate macroscopic phase separation for phospholipid mixtures composed of species with different hydrocarbon tail groups. We also investigate the relationship between the phase segregation and thermodynamic variables such as interfacial line tension and surface tension, and obtain correspondence between existing theory and experiments, and our simulation results. We report variations in the molecular chain stiffness to have negligible contributions to the phase segregation in the mixed bilayer, and to demonstrate shape transformations of the hybrid vesicle. Our results can be used to design novel bio-inspired hybrid vehicles for potential applications in biomedicine, sensing, imaging and sustainability.

Perspective: Coarse-grained models for biomolecular systems

By focusing on essential features, while averaging over less important details, coarse-grained (CG) models provide significant computational and conceptual advantages with respect to more detailed models. Consequently, despite dramatic advances in computational methodologies and resources, CG models enjoy surging popularity and are becoming increasingly equal partners to atomically detailed models. This perspective surveys the rapidly developing landscape of CG models for biomolecular systems. In particular, this review seeks to provide a balanced, coherent, and unified presentation of several distinct approaches for developing CG models, including top-down, network-based, native-centric, knowledge-based, and bottom-up modeling strategies. The review summarizes their basic philosophies, theoretical foundations, typical applications, and recent developments. Additionally, the review identifies fundamental inter-relationships among the diverse approaches and discusses outstanding challenges in the field. When carefully applied and assessed, current CG models provide highly efficient means for investigating the biological consequences of basic physicochemical principles. Moreover, rigorous bottom-up approaches hold great promise for further improving the accuracy and scope of CG models for biomolecular systems.

Scaling behavior and local structure of ion aggregates in single-ion conductors

Single-ion conductors are attractive electrolyte materials because of their inherent safety and ease of processing. Most ions in a sodium-neutralized PEO sulfonated-isophthalate ionomer electrolyte exist as one dimensional chains, restricted in dimensionality by the steric hindrance of the attached polymer. Because the ions are slow to reconfigure, atomistic MD simulations of this material are unable to adequately sample equilibrium ion structures. We apply a novel coarse-graining scheme using a generalized-YBG procedure in which the polymer backbone is completely removed and implicitly represented by the effective potentials of the remaining ions. The ion-only coarse-grained simulation allows for substantial sampling of equilibrium aggregate configurations. We extend the wormlike micelle theory to model ion chain equilibrium. Our aggregates are random walks which become more positively charged with increasing size. Defects occur on the string-like structure in the form of “dust” and “knots,” which form due to cation coordination with open sites along the string. The presence of these defects suggest that cation hopping along open third-coordination sites could be an important mechanism of charge transport using ion aggregates.

Systematic coarse-graining of the dynamics of entangled polymer melts: the road from chemistry to rheology

For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together (up to five backbone carbon atoms) and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together (more than ten backbone carbon atoms), leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and action must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.

Why are coarse-grained force fields too fast? A look at dynamics of four coarse-grained polymers

Coarse-grained models decrease the number of force sites and thus reduce computational requirements for molecular simulation. While these models are successful in describing structural properties, dynamic evolution is faster than the corresponding atomistic simulations or experiments. We consider coarse-grained models for four polymers and one polymer mixture, where accurate dynamics are obtained by scaling to match the mean-squared displacements of the corresponding atomistic descriptions. We show that the required scaling is dictated by local friction and that this scaling is only valid after the onset of continuous motion.

Coarse-grained simulations of rapid assembly kinetics for polystyrene-b-poly(ethylene oxide) copolymers in aqueous solutions

We present a coarse-grained, implicit solvent model for polystyrene-b-poly(ethylene oxide) in aqueous solution and study its assembly kinetics using Brownian dynamics simulations. The polymer is modeled as a chain of freely jointed beads interacting through effective potentials. Coarse-grained force field parameters are determined by matching experimental thermodynamic quantities including radius of gyration, second virial coefficient, aggregation number, and critical micelle concentration. We investigate the influence of cooling rate (analogous to the rate of solvent quality change in rapid precipitations), polymer concentration, and friction coefficient on the assembly kinetics and compare simulation results to flash nanoprecipitation experiments. We find that assembly kinetics show a linear scaling relation with inverse friction coefficient when the friction coefficient is larger than 1. When the cooling time is less than the characteristic micellization time, stable kinetically arrested clusters are obtained; otherwise, close-to-equilibrium micelles are formed. The characteristic micellization time is estimated to be only 3-6 ms, in contrast to 30-40 ms previously determined in experiments. We suggest that previous experiments probed the formation of micellar clusters while simulations in this work studied the kinetics of a single micelle assembled from free polymer chains.

The phase behavior of polyethylene ring chains

The equilibrium properties of an isolated polyethylene ring chain are studied by using molecular dynamics (MD) simulations. The results of an 80-bond linear chain are also presented, which are in agreement with previous studies of square-well chains and Lennard-Jones (LJ) homopolymers. Mainly, we focus on the collapse of polyethylene ring chains. At high temperatures, a fully oblate structure is observed for the ring chains with different chain lengths. For such an oblate structure, a shape factor of delta(*)=0.25 and a rodlike scaling relation between the radius of gyration and chain lengths could be deduced easily in theory, and the same results are obtained by our MD simulations. Such an oblate structure can be obtained by Monte Carlo simulation only for sufficient stiff ring chains. When the temperature decreases, an internal energy barrier is observed. This induces a strong peak in the heat capacity, denoting a gas-liquid-like transition. This energy barrier comes mainly from the local monomer-monomer interactions, i.e., the bond-stretching, the bond-bending, and the torsion potentials. A low temperature peak is also observed in the same heat capacity curve, representing a liquid-solid-like transition. These numerical simulation results support a two-stage collapse of polyethylene ring chains; however, the nature should be different from the square-well and LJ ring chains.