No Arabic abstract
Following the Flory ideality hypothesis intrachain and interchain excluded volume interactions are supposed to compensate each other in dense polymer systems. Multi-chain effects should thus be neglected and polymer conformations may be understood from simple phantom chain models. Here we provide evidence against this phantom chain, mean-field picture. We analyze numerically and theoretically the static correlation function of the Rouse modes. Our numerical results are obtained from computer simulations of two coarse-grained polymer models for which the strength of the monomer repulsion can be varied, from full excluded volume (`hard monomers) to no excluded volume (`phantom chains). For nonvanishing excluded volume we find the simulated correlation function of the Rouse modes to deviate markedly from the predictions of phantom chain models. This demonstrates that there are nonnegligible correlations along the chains in a melt. These correlations can be taken into account by perturbation theory. Our simulation results are in good agreement with these new theoretical predictions.
Nonlinear extensional flows are common in polymer processing but remain challenging theoretically because dramatic stretching of chains deforms the entanglement network far from equilibrium. Here, we present coarse-grained simulations of extensional flows in entangled polymer melts for Rouse-Weissenberg numbers $Wi_R=0.06$-$52$ and Hencky strains $epsilongeq6$. Simulations reproduce experimental trends in extensional viscosity with time, rate and molecular weight. Studies of molecular structure reveal an elongation and thinning of the confining tube with increasing $Wi_R$. The rising stress is quantitatively consistent with the decreasing entropy of chains at the equilibrium entanglement length. Molecular weight dependent trends in viscosity are related to a crossover from the Newtonian limit to a high rate limit that scales differently with chain length.
Conformation-dependent design of polymer sequences can be considered as a tool to control macromolecular self-assembly. We consider the monomer unit sequences created via the modification of polymers in a homogeneous melt in accordance with the spatial positions of the monomer units. The geometrical patterns of lamellae, hexagonally packed cylinders, and balls arranged in a body-centered cubic lattice are considered as typical microphase-separated morphologies of block copolymers. Random trajectories of polymer chains are described by the diffusion-type equations and, in parallel, simulated in the computer modeling. The probability distributions of block length $k$, which are analogous to the first-passage probabilities, are calculated analytically and determined from the computer simulations. In any domain, the probability distribution can be described by the asymptote $~k^{-3/2}$ at moderate values of $k$ if the spatial size of the block is less than the smallest characteristic size of the domain. For large blocks, the exponential asymptote $exp(-const , k a^2/d_{as}^2)$ is valid, $d_{as}$ being the asymptotic domain length (a is the monomer unit size). The number average block lengths and their dispersities change linearly with the block length for lamellae, cylinders, and balls, when the domain is characterized by a single characteristic size. If the domain is described by more than one size, the number average block length can grow nonlinearly with the domain sizes and the length das can depend on all of them.
Based on non-equilibrium molecular dynamics simulations of entangled polymer melts, a recent Letter [Phys. Rev. Lett. $textbf{121}$, 047801 (2018), arXiv:1806.09509] claims that the rising extensional stress is quantitatively consistent with the decreasing entropy of chains at the equilibrium entanglement length. We point out that exactly the opposite is true: the intrachain entropic stress arising from individual entanglement strands generally does not agree with the total macroscopic stress. The conclusion of the Letter is based on an incomplete and questionable analysis of a limited range of the simulation trajectory. The opposite conclusion should have been drawn from their data, had they examined the full simulation trajectory in a proper way.
We study the relaxation dynamics of a coarse-grained polymer chain at different degrees of stretching by both analytical means and numerical simulations. The macromolecule is modelled as a string of beads, connected by anharmonic springs, subject to a tensile force applied at the end monomer of the chain while the other end is fixed at the origin of coordinates. The impact of bond non-linearity on the relaxation dynamics of the polymer at different degrees of stretching is treated analytically within the Gaussian self-consistent approach (GSC) and then compared to simulation results derived from two different methods: Monte-Carlo (MC) and Molecular Dynamics (MD). At low and medium degrees of chain elongation we find good agreement between GSC predictions and the Monte-Carlo simulations. However, for strongly stretched chains the MD method, which takes into account inertial effects, reveals two important aspects of the nonlinear interaction between monomers: (i) a coupling and energy transfer between the damped, oscillatory normal modes of the chain, and (ii) the appearance of non-vanishing contributions of a continuum of frequencies around the characteristic modes in the power spectrum of the normal mode correlation functions.
We present an effective and simple multiscale method for equilibrating Kremer Grest model polymer melts of varying stiffness. In our approach, we progressively equilibrate the melt structure above the tube scale, inside the tube and finally at the monomeric scale. We make use of models designed to be computationally effective at each scale. Density fluctuations in the melt structure above the tube scale are minimized through a Monte Carlo simulated annealing of a lattice polymer model. Subsequently the melt structure below the tube scale is equilibrated via the Rouse dynamics of a force-capped Kremer-Grest model that allows chains to partially interpenetrate. Finally the Kremer-Grest force field is introduced to freeze the topological state and enforce correct monomer packing. We generate $15$ melts of $500$ chains of $10.000$ beads for varying chain stiffness as well as a number of melts with $1.000$ chains of $15.000$ monomers. To validate the equilibration process we study the time evolution of bulk, collective and single-chain observables at the monomeric, mesoscopic and macroscopic length scales. Extension of the present method to longer, branched or polydisperse chains and/or larger system sizes is straight forward.