Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userMultiscale approach to equilibrating model polymer melts
58
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
It is commonly accepted that in concentrated solutions or melts high-molecular weight polymers display random-walk conformational properties without long-range correlations between subsequent bonds. This absence of memory means, for instance, that the bond-bond correlation function, $P(s)$, of two bonds separated by $s$ monomers along the chain should exponentially decay with $s$. Presenting numerical results and theoretical arguments for both monodisperse chains and self-assembled (essentially Flory size-distributed) equilibrium polymers we demonstrate that some long-range correlations remain due to self-interactions of the chains caused by the chain connectivity and the incompressibility of the melt. Suggesting a profound analogy with the well-known long-range velocity correlations in liquids we find, for instance, $P(s)$ to decay algebraically as $s^{-3/2}$. Our study suggests a precise method for obtaining the statistical segment length bstar in a computer experiment.
We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, $G(t)$, into the plateau regime for chains with $Z=40$ entanglements and into the terminal relaxation regime for $Z=10$. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter -free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.
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.
We present a numerical study of the slip link model introduced by Likhtman for describing the dy- namics of dense polymer melts. After reviewing the technical aspects associated with the implemen- tation of the model, we extend previous work in several directions. The dependence of the relaxation modulus with the slip link density and the slip link stiffness is reported. Then the nonlinear rheolog- ical properties of the model, for a particular set of parameters, are explored. Finally, we introduce excluded volume interactions in a mean field such as manner in order to describe inhomogeneous systems, and we apply this description to a simple nanocomposite model. With this extension, the slip link model appears as a simple and generic model of a polymer melt, that can be used as an alternative to molecular dynamics for coarse grained simulations of complex polymeric systems.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
