No Arabic abstract
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.
The temperature dependence of the hydrodynamic boundary condition between a PDMS melt and two different non-attractive surfaces made of either an OTS (octadecyltrichlorosilane) self-assembled monolayer (SAM) or a grafted layer of short PDMS chains has been characterized. A slip length proportional to the fluid viscosity is observed on both surfaces. The slip temperature dependence is deeply influenced by the surfaces. The viscous stress exerted by the polymer liquid on the surface is observed to follow exactly the same temperature dependences as the friction stress of a cross-linked elastomer sliding on the same surfaces. Far above the glass transition temperature, these observations are rationalized in the framework of a molecular model based on activation energies: increase or decrease of the slip length with increasing temperatures can be observed depending on how the activation energy of the bulk viscosity compares to that of the interfacial Naviers friction coefficient.
Hydrodynamic slip of a liquid at a solid surface represents a fundamental phenomenon in fluid dynamics that governs liquid transport at small scales. For polymeric liquids, de Gennes predicted that the Navier boundary condition together with the theory of polymer dynamics imply extraordinarily large interfacial slip for entangled polymer melts on ideal surfaces; this Navier-de Gennes model was confirmed using dewetting experiments on ultra-smooth, low-energy substrates. Here, we use capillary leveling - surface tension driven flow of films with initially non-uniform thickness - of polymeric films on these same substrates. Measurement of the slip length from a robust one-parameter fit to a lubrication model is achieved. We show that at the lower shear rates involved in leveling experiments as compared to dewetting ones, the employed substrates can no longer be considered ideal. The data is instead consistent with physical adsorption of polymer chains at the solid/liquid interface. We extend the Navier-de Gennes description using one additional parameter, namely the density of physically adsorbed chains per unit surface. The resulting formulation is found to be in excellent agreement with the experimental observations.
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.
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.
Following Florys ideality hypothesis the chemical potential of a test chain of length $n$ immersed into a dense solution of chemically identical polymers of length distribution P(N) is extensive in $n$. We argue that an additional contribution $delta mu_c(n) sim +1/rhosqrt{n}$ arises ($rho$ being the monomer density) for all $P(N)$ if $n ll <N>$ which can be traced back to the overall incompressibility of the solution leading to a long-range repulsion between monomers. Focusing on Flory distributed melts we obtain $delta mu_c(n) approx (1- 2 n/<N>) / rho sqrt{n}$ for $n ll <N>^2$, hence, $delta mu_c(n) approx - 1/rho sqrt{n}$ if $n$ is similar to the typical length of the bath $<N>$. Similar results are obtained for monodisperse solutions. Our perturbation calculations are checked numerically by analyzing the annealed length distribution P(N) of linear equilibrium polymers generated by Monte Carlo simulation of the bond-fluctuation model. As predicted we find, e.g., the non-exponentiality parameter $K_p equiv 1 - <N^>/p!<N>^p$ to decay as $K_p approx 1 / sqrt{<N>}$ for all moments $p$ of the distribution.