Self-avoiding polymers in two-dimensional ($d=2$) melts are known to adopt compact configurations of typical size $R(N) sim N^{1/d}$ with $N$ being the chain length. Using molecular dynamics simulations we show that the irregular shapes of these chains are characterized by a perimeter length $L(N) sim R(N)^{dpm}$ of fractal dimension $dpm = d-Theta_2 =5/4$ with $Theta_2=3/4$ being a well-known contact exponent. Due to the self-similar structure of the chains, compactness and perimeter fractality repeat for subchains of all arc-lengths $s$ down to a few monomers. The Kratky representation of the intramolecular form factor $F(q)$ reveals a strong non-monotonous behavior with $q^2F(q) sim 1/(qN^{1/d})^{Theta_2}$ in the intermediate regime of the wavevector $q$. Measuring the scattering of labeled subchains %($s F(q) sim L(s)$) the form factor may allow to test our predictions in real experiments.
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.
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.
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.