Coarse-Grained Description of Polymer Liquids and their Mixtures as Interacting Soft-Colloidal Particles (original) (raw)
Related papers
Accuracy, Transferability, and Efficiency of Coarse-Grained Models of Molecular Liquids
The Journal of Physical Chemistry B
Coarse-graining (CG) approaches are becoming essential tools in the study of complex systems because they can considerably speed up computer simulations, with the promise of determining properties in a range of length scales and time scales never before possible. While much progress in this field has been achieved in recent years, application of CG methods is still inhibited by the limited understanding of a number of conceptual points that need to be resolved to open up the field of CG to a wide range of applications in material science and biology. In this paper, we present some of the key findings that emerged from the development of the integral equation theory of coarse-graining (IECG), which addresses some of the fundamental questions in coarse-graining. Although the IECG method pertains to the CG of polymer liquids, and specifically homopolymer melts are illustrated here, many of the results that emerge from the study of the IECG approach are general and apply to the CG of any molecular liquid. Through this method, we developed a formal relation between the statistical mechanics of CG and a number of predicted physical properties. On the basis of the theory of liquids, the IECG affords the analytical solution of the intermolecular potential for macromolecules represented by a Markov chain of CG sites, thus providing a transparent tool for analysis of the properties in coarse-graining. We identify three key requirements that render a CG model useful: accuracy, transferability, and computational efficiency. When these three requirements are fulfilled, the CG model becomes widely applicable and useful for studying regions in the phase space that are not covered by atomistic simulations. In the process, the IECG answers formally a number of relevant questions on how structural, thermodynamic, and dynamical properties are modified during coarse-graining. It sheds light upon how the level of CG affects the shape of the CG potential and how, in turn, the shape of the potential affects the physical properties. It tests the validity of selecting the potential-of-mean force as the effective pairwise CG potential and the role of higher-order many-body corrections to the pairwise potential to recover structural and thermodynamic consistency of the CG model. Because the IECG theory can be analytically formalized, it does not suffer from the problem of transferability and, in the canonical ensemble, leads to consistent pair distribution functions, pressure, isothermal compressibility, and excess free energy at variable levels of CG from the atomistic to the ultra-CG model, where macromolecules are represented as interpenetrable soft spheres.
Multiscale Modeling of Coarse-Grained Macromolecular Liquids
The Journal of Physical Chemistry B, 2009
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.
Coarse-grained description of polymer blends as interacting soft-colloidal particles
The Journal of Chemical Physics, 2005
We present a theoretical approach which maps polymer blends onto mixtures of soft-colloidal particles. The analytical mesoscale pair distribution functions reproduce well data from united atom molecular dynamics simulations of polyolefin mixtures without fitting parameters. The theory exactly recovers the analytical expressions for density and concentration fluctuation structure factors of soft colloidal mixtures (liquid alloys).
Effective Soft-Core Potentials and Mesoscopic Simulations of Binary Polymer Mixtures
Macromolecules, 2010
Mesoscopic molecular dynamics simulations are used to determine the large scale structure of several binary polymer mixtures of various chemical architecture, concentration, and thermodynamic conditions. By implementing an analytical formalism, which is based on the solution to the Ornstein-Zernike equation, each polymer chain is mapped onto the level of a single soft colloid. From the appropriate closure relation, the effective, soft-core potential between coarse-grained units is obtained and used as input to our mesoscale simulations. The potential derived in this manner is analytical and explicitly parameter dependent, making it general and transferable to numerous systems of interest. From computer simulations performed under various thermodynamic conditions the structure of the polymer mixture, through pair correlation functions, is determined over the entire miscible region of the phase diagram. In the athermal regime mesoscale simulations exhibit quantitative agreement with united atom simulations. Furthermore, they also provide information at larger scales than can be attained by united atom simulations and in the thermal regime approaching the phase transition.
Analytical Soft-Core Potentials for Macromolecular Fluids and Mixtures
Physical Review Letters, 2004
An analytical description of polymer melts and their mixtures as liquids of interacting soft colloidal particles is obtained from liquid-state theory. The derived center-of-mass pair correlation functions with no adjustable parameters reproduce those computed from united atom molecular dynamics simulations. The coarse-grained description correctly bridges micro-and mesoscopic fluid properties. Molecular dynamics simulations of soft colloidal particles interacting through the calculated effective pair potentials are consistent with data from microscale simulations and analytical formulas.
Theoretical models for bridging timescales in polymer dynamics
Journal of Physics: Condensed Matter, 2008
The dynamics of macromolecules are characterized by the presence of several length scales and related timescales in which relevant phenomena take place. This defines the complex nature of the liquid and renders its theoretical treatment a difficult matter. The necessity of developing theoretical approaches that can describe in a comprehensive manner properties observed at many different length scales is a fundamental challenge in polymer physics.
The Journal of Chemical Physics, 2013
The theory to reconstruct the atomistic-level chain diffusion from the accelerated dynamics that is measured in mesoscale simulations of the coarse-grained system, is applied here to the dynamics of cis-1,4-polybutadiene melts where each chain is described as a soft interacting colloidal particle. The rescaling formalism accounts for the corrections in the dynamics due to the change in entropy and the change in friction that are a consequence of the coarse-graining procedure. By including these two corrections the dynamics is rescaled to reproduce the realistic dynamics of the system described at the atomistic level. The rescaled diffusion coefficient obtained from mesoscale simulations of coarse-grained cis-1,4-polybutadiene melts shows good agreement with data from united atom simulations performed by Tsolou et al. [Macromolecules 38, 1478]. The derived monomer friction coefficient is used as an input to the theory for cooperative dynamics that describes the internal dynamics of a polymer moving in a transient regions of slow cooperative motion in a liquid of macromolecules. Theoretically predicted time correlation functions show good agreement with simulations in the whole range of length and time scales in which data are available.
Theoretical coarse-graining approach to bridge length scales in diblock copolymer liquids
Physical Review E, 2007
A microscopic theory for coarse graining diblock copolymers into dumbbells of interacting soft colloidal particles has been developed, based on the solution of liquid-state integral equations. The Ornstein-Zernike equation is solved to provide a mesoscopic description of the diblock copolymer system at the level of block centers of mass, and at the level of polymer centers of mass. Analytical forms of the total correlation functions for block-block, block-monomer, and center-of-mass pairs are obtained for a liquid of structurally symmetric diblock copolymers as a function of temperature, density, chain length, and chain composition. The theory correctly predicts thermodynamicallydriven segregation of diblocks into microdomains as a function of temperature (chi parameter). The coarse-grained description contains contributions from density and concentration fluctuations, with the latter becoming dominant as temperature decreases. Numerical calculations for the block coarse-grained total correlation functions, as a function of the proximity of the system to its phase transition, are presented. Comparison with united atom molecular dynamics simulations are carried out in the athermal regime, where simulations and theory quantitatively agree with no need of adjustable parameters. PACS numbers: 61.20.Gy/61.25.Em/61.25.Hq/83.80.Sg
Theoretical Reconstruction of Realistic Dynamics of Polymer Melts from Soft Sphere Representation
2011
The theory to reconstruct the atomistic-level chain diffusion from the accelerated dynamics that is measured in mesoscale simulations of the coarse-grained system, is applied here to the dynamics of cis-1,4-polybutadiene melts where each chain is described as a soft interacting colloidal particle. The rescaling formalism accounts for the corrections in the dynamics due to the change in entropy and the change in friction that are a consequence of the coarse-graining procedure. By including these two corrections the dynamics is rescaled to reproduce the realistic dynamics of the system described at the atomistic level. The rescaled diffusion coefficient obtained from mesoscale simulations of coarse-grained cis-1,4-polybutadiene melts shows good agreement with data from united atom simulations performed by Tsolou et al. [Macromolecules 38, 1478]. The derived monomer friction coefficient is used as an input to the theory for cooperative dynamics that describes the internal dynamics of a polymer moving in a transient regions of slow cooperative motion in a liquid of macromolecules. Theoretically predicted time correlation functions show good agreement with simulations in the whole range of length and time scales in which data are available.
Mapping of polymer melts onto liquids of soft-colloidal chains
The Journal of Chemical Physics, 2010
Microscopic computer simulations of fluids of long polymers are greatly restricted by the limits of current computational power, and so course-grained descriptions, accurate on molecular length scales, are essential to extending the range of accessible systems. For some phenomena, particularly dynamical entanglement, descriptions that eliminate all internal degrees of freedom from the polymers are too drastic, as intermediate wavelength degrees of freedom are essential to the effect. Employing first-principles liquid-state theory, we have developed a course-grained model for the intermolecular structure of melts of long homopolymer chains that maps each chain of hard-sphere monomers onto a chain of connected soft colloids. All dependence on system parameters is analytically expressed so the results may be immediately applied to melts with different polymer and thermodynamic properties to calculate effective potentials between the soft colloids on the chains, which can then be used to perform molecular dynamics simulations. These simulations will be able to capture the large wavelength structure of the system at greatly reduced computational cost, while still retaining enough internal degrees of freedom explicitly to describe the phenomena that occur on length scales much larger than the monomeric units that comprise the chain, but shorter than the size of the molecule.