No Arabic abstract
We present a theoretical approach to scale the artificially fast dynamics of simulated coarse-grained polymer liquids down to its realistic value. As coarse-graining affects entropy and dissipation, two factors enter the rescaling: inclusion of intramolecular vibrational degrees of freedom, and rescaling of the friction coefficient. Because our approach is analytical, it is general and transferable. Translational and rotational diffusion of unentangled and entangled polyethylene melts, predicted from mesoscale simulations of coarse-grained polymer melts using our rescaling procedure, are in quantitative agreement with united atom simulations and with experiments.
Structural and thermodynamic consistency of coarse-graining models across multiple length scales is essential for the predictive role of multi-scale modeling and molecular dynamic simulations that use mesoscale descriptions. Our approach is a coarse-grained model based on integral equation theory, which can represent polymer chains at variable levels of chemical details. The model is analytical and depends on molecular and thermodynamic parameters of the system under study, as well as on the direct correlation function in the k --> 0 limit, c0. A numerical solution to the PRISM integral equations is used to determine c0, by adjusting the value of the effective hard sphere diameter, d, to agree with the predicted equation of state. This single quantity parameterizes the coarse-grained potential, which is used to perform mesoscale simulations that are directly compared with atomistic-level simulations of the same system. We test our coarse-graining formalism by comparing structural correlations, isothermal compressibility, equation of state, Helmholtz and Gibbs free energies, and potential energy and entropy using both united atom and coarse-grained descriptions. We find quantitative agreement between the analytical formalism for the thermodynamic properties, and the results of Molecular Dynamics simulations, independent of the chosen level of representation. In the mesoscale description, the potential energy of the soft-particle interaction becomes a free energy in the coarse-grained coordinates which preserves the excess free energy from an ideal gas across all levels of description. The structural consistency between the united-atom and mesoscale descriptions means the relative entropy between descriptions has been minimized without any variational optimization parameters. The approach is general and applicable to any polymeric system in different thermodynamic conditions.
A lattice model is presented for the simulation of dynamics in polymeric systems. Each polymer is represented as a chain of monomers, residing on a sequence of nearest-neighbor sites of a face-centered-cubic lattice. The polymers are self- and mutually avoiding walks: no lattice site is visited by more than one polymer, nor revisited by the same polymer after leaving it. The dynamics occurs through single-monomer displacements over one lattice spacing. To demonstrate the high computational efficiency of the model, we simulate a dense binary polymer mixture with repelling nearest-neighbor interactions between the two types of polymers, and observe the phase separation over a long period of time. The simulations consist of a total of 46,080 polymers, 100 monomers each, on a lattice with 13,824,000 sites, and an interaction strength of 0.1 kT. In the final two decades of time, the domain-growth is found to be d(t) ~ t^1/3, as expected for a so-called Model B system.
Significant progress was made in recent years in the understanding of the proton spin kinetics in polymer melts. Generally, the proton spin kinetics is determined by intramolecular and intermolecular magnetic dipole-dipole contributions of proton spins. During many decades it was postulated that the main contribution is a result of intramolecular magnetic dipole-dipole interactions of protons belonging to the same polymer segment. It appears that this postulate is far from reality. The relative weights of intra- and intermolecular contributions are time dependent and sensitive to details of polymer chain dynamics. It is shown that for isotropic models of polymer dynamics the influence of the intermolecular magnetic dipole-dipole interactions increases faster with increasing evolution time (i.e. decreasing frequency) than the corresponding influence of the intramolecular counterpart. On the other hand, an inverted situation is predicted by the tube-reptation model: here the influence of the intramolecular magnetic dipole-dipole interactions increases faster with time than the contribution from intermolecular interactions. The intermolecular contribution in the proton relaxation of polymer melts can experimentally be isolated using the isotope dilution technique and this opens a new perspective for experimental investigations of polymer dynamics by proton NMR.
The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.
We study the conformational dynamics within homo-polymer globules by solvent-implicit Brownian dynamics simulations. A strong dependence of the internal chain dynamics on the Lennard-Jones cohesion strength {epsilon} and the globule size NG is observed. We find two distinct dynamical regimes: a liquid- like regime (for {epsilon} < {epsilon}s) with fast internal dynamics and a solid-like regime (for {epsilon} > {epsilon}s) with slow internal dynamics. The cohesion strength {epsilon}s of this freezing transition depends on NG. Equilibrium simulations, where we investigate the diffusional chain dynamics within the globule, are compared with non-equilibrium simulations, where we unfold the globule by pulling the chain ends with prescribed velocity (encompassing low enough velocities so that the linear-response, viscous regime is reached). From both simulation protocols we derive the internal viscosity within the globule. In the liquid-like regime the internal friction increases continuously with {epsilon} and scales extensive in NG. This suggests an internal friction scenario where the entire chain (or an extensive fraction thereof) takes part in conformational reorganization of the globular structure.