No Arabic abstract
A first-principle multiscale modeling approach is presented, which is derived from the solution of the Ornstein-Zernike equation for the coarse-grained representation of polymer liquids. The approach is analytical, and for this reason is transferable. It is here applied to determine the structure of several polymeric systems, which have different parameter values, such as molecular length, monomeric structure, local flexibility, and thermodynamic conditions. When the pair distribution function obtained from this procedure is compared with the results from a full atomistic simulation, it shows quantitative agreement. Moreover, the multiscale procedure accurately captures both large and local scale properties while remaining computationally advantageous.
We present a detailed derivation and testing of our approach to rescale the dynamics of mesoscale simulations of coarse-grained polymer melts (I. Y. Lyubimov et al. J. Chem. Phys. textbf{132}, 11876, 2010). Starting from the first-principle Liouville equation and applying the Mori-Zwanzig projection operator technique, we derive the Generalized Langevin Equations (GLE) for the coarse-grained representations of the liquid. The chosen slow variables in the projection operators define the length scale of coarse graining. Each polymer is represented at two levels of coarse-graining: monomeric as a bead-and-spring model and molecular as a soft-colloid. In the long-time regime where the center-of-mass follows Brownian motion and the internal dynamics is completely relaxed, the two descriptions must be equivalent. By enforcing this formal relation we derive from the GLEs the analytical rescaling factors to be applied to dynamical data in the coarse-grained representation to recover the monomeric description. Change in entropy and change in friction are the two corrections to be accounted for to compensate the effects of coarse-graining on the polymer dynamics. The solution of the memory functions in the coarse-grained representations provides the dynamical rescaling of the friction coefficient. The calculation of the internal degrees of freedom provides the correction of the change in entropy due to coarse-graining. The resulting rescaling formalism is a function of the coarse-grained model and thermodynamic parameters of the system simulated. The rescaled dynamics obtained from mesoscale simulations of polyethylene, represented as soft colloidal particles, by applying our rescaling approach shows a good agreement with data of translational diffusion measured experimentally and from simulations. The proposed method is used to predict self-diffusion coefficients of new polyethylene samples.
Using concepts from integral geometry, we propose a definition for a local coarse-grained curvature tensor that is well-defined on polygonal surfaces. This coarse-grained curvature tensor shows fast convergence to the curvature tensor of smooth surfaces, capturing with accuracy not only the principal curvatures but also the principal directions of curvature. Thanks to the additivity of the integrated curvature tensor, coarse-graining procedures can be implemented to compute it over arbitrary patches of polygons. When computed for a closed surface, the integrated curvature tensor is identical to a rank-2 Minkowski tensor. We also provide an algorithm to extend an existing C++ package, that can be used to compute efficiently local curvature tensors on triangulated surfaces.
Integral equation theory is applied to a coarse-grained model of water to study potential of mean force between hydrophobic solutes. Theory is shown to be in good agreement with the available simulation data for methane-methane and fullerene-fullerene potential of mean force in water; the potential of mean force is also decomposed into its entropic and enthalpic contributions. Mode coupling theory is employed to compute self-diffusion coefficient of water, as well as diffusion coefficient of a dilute hydrophobic solute; good agreement with molecular dynamics simulation results is found.
Room-temperature ionic liquids (RTILs) stand out among molecular liquids for their rich physicochemical characteristics, including structural and dynamic heterogeneity. The significance of electrostatic interactions in RTILs results in long characteristic length- and timescales, and has motivated the development of a number of coarse-grained (CG) simulation models. In this study, we aim to better understand the connection between certain CG parametrization strategies and the dynamical properties and transferability of the resulting models. We systematically compare five CG models: a model largely parametrized from experimental thermodynamic observables; a refinement of this model to increase its structural accuracy; and three models that reproduce a given set of structural distribution functions by construction, with varying intramolecular parametrizations and reference temperatures. All five CG models display limited structural transferability over temperature, and also result in various effective dynamical speedup factors, relative to a reference atomistic model. On the other hand, the structure-based CG models tend to result in more consistent cation-anion relative diffusion than the thermodynamic-based models, for a single thermodynamic state point. By linking short- and long-timescale dynamical behaviors, we demonstrate that the varying dynamical properties of the different coarse-grained models can be largely collapsed onto a single curve, which provides evidence for a route to constructing dynamically-consistent CG models of RTILs.
We study DNA self-assembly and DNA computation using a coarse-grained DNA model within the directional dynamic bonding framework {[}C. Svaneborg, Comp. Phys. Comm. 183, 1793 (2012){]}. In our model, a single nucleotide or domain is represented by a single interaction site. Complementary sites can reversibly hybridize and dehybridize during a simulation. This bond dynamics induces a dynamics of the angular and dihedral bonds, that model the collective effects of chemical structure on the hybridization dynamics. We use the DNA model to perform simulations of the self-assembly kinetics of DNA tetrahedra, an icosahedron, as well as strand displacement operations used in DNA computation.