Analytical Rheology of Polymer Melts: State of the Art

The extreme sensitivity of rheology to the microstructure of polymer melts has prompted the development of “analytical rheology,” which seeks inferring the structure and composition of an unknown sample based on rheological measurements. Typically, this involves the inversion of a model, which may be mathematical, computational, or completely empirical. Despite the imperfect state of existing models, analytical rheology remains a practically useful enterprise. I review its successes and failures in inferring the molecular weight distribution of linear polymers and the branching content in branched polymers.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.

Topological analysis of polymeric melts: chain-length effects and fast-converging estimators for entanglement length

Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead spring and atomistic polyethylene polymer melts. Estimators for the entanglement length N_{e} which operate on results for a single chain length N are shown to produce systematic O(1/N) errors. The mathematical roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian corrections to chain and primitive path dimensions. The prefactors for the O(1/N) errors may be large; in general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We propose, derive, and test new estimators which eliminate these systematic errors using information obtainable from the variation in entanglement characteristics with chain length. The new estimators produce accurate results for N_{e} from marginally entangled systems. Formulas based on direct enumeration of entanglements appear to converge faster and are simpler to apply.

Screening of hydrodynamic interactions in Brownian rod suspensions

We present the details and results of a simulation study addressing the dynamics and rheology of rod suspensions over a wide regime of concentrations ranging from dilute to concentrated systems. Our study compares the results of two complementary simulation methods. The first method adapts a recently proposed explicit solvent simulation strategy and incorporates both hydrodynamical effects and steric interactions between the rod units. We compare the results of such a method with those obtained from a Brownian dynamics simulation approach which retains the steric interactions but neglects the effects of hydrodynamic interactions. Overall, our results in the context of the translational and rotational diffusivities are in agreement with the hydrodynamical predictions in the dilute regime and the corresponding results of the tube model and its extensions thereof in the semidilute regimes. The latter results suggest that effects of hydrodynamic interactions on the translational and rotational diffusivities are secondary relative to the steric interactions and at best lead only to a small correction to the results of the classical tube model. Our results in the context of linear viscoelasticity also broadly confirms the predictions of the tube model for the storage and loss moduli and allows us to extract for the first time the independent hydrodynamic and Brownian contributions to the zero shear viscosity. While the relative magnitudes of these contributions are consistent with the theoretical predictions, the quantitative magnitudes are quite different from the theoretical predictions. Overall, these results confirm the validity of the hydrodynamic "screening" hypothesis and ratify the neglect of hydrodynamical stresses in quantifying the linear rheology of Brownian rod suspensions.

Analytic expressions for the statistics of the primitive-path length in entangled polymers

An analytic expression is proposed for the primitive-path length of entangled polymer chains. The expression is derived from statistical mechanics of a chain that is a random walk with randomly scattered entanglements. The only parameters are the number of Kuhn steps in the chain and a dimensionless parameter beta that contains information about the entanglement density and Kuhn step size. The expression is found to compare very favorably with numerical results recently found from examining topological constraints in microscopic simulations. The comparison also predicts well the plateau modulus of polyethylene, suggesting that the slip-link model is a viable intermediate in the search for true ab initio rheology predictions. Since the expression is analytic, it can be used to make predictions where the simulations cannot reach, and hence is applicable for coarse graining.

Conformational properties of blends of cyclic and linear polymer melts

An adapted version of the annealing algorithm to identify primitive paths of a melt of ring polymers is presented. This algorithm ensures that the primitive path length becomes zero for nonconcatenated rings, and that no entanglements are observed. The bond-fluctuation model was used to simulate ring-linear blends with N=150 and 300 monomers. The primitive path length and the average number of entanglements of the linear component were found to be independent of the blend composition. In contrast, the primitive path length and the average number of entanglements on a ring molecule increased approximately linearly with the fraction of linear chains, and for large N , they approached values comparable with linear chains. Threading of ring molecules by linear chains, and ring-ring interactions were observed only in the presence of linear chains. It is conjectured that for large N , these latter interactions facilitate the formation of a percolating entangled network, thereby resulting in a disproportionate retardation of the dynamical processes.