Dynamics of short polymer chains in solution

We present numerical and analytical results describing the effect of hydrodynamic interactions on the dynamics of a short polymer chain in solution. A molecular dynamics algorithm for the polymer is coupled to a direct simulation Monte Carlo algorithm for the solvent. We give an explicit expression for the velocity autocorrelation function of the centre of mass of the polymer which agrees well with numerical results if Brownian dynamics, hydrodynamic correlations and sound wave scattering are included

Systematic coarse-graining of the dynamics of entangled polymer melts: the road from chemistry to rheology

For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together, leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and actions must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.Comment: Topical review, 31 pages, 10 figure

Coexistence and Phase Separation in Sheared Complex Fluids

We demonstrate how to construct dynamic phase diagrams for complex fluids that undergo transitions under flow, in which the conserved composition variable and the broken-symmetry order parameter (nematic, smectic, crystalline, etc.) are coupled to shear rate. Our construction relies on a selection criterion, the existence of a steady interface connecting two stable homogeneous states. We use the (generalized) Doi model of lyotropic nematic liquid crystals as a model system, but the method can be easily applied to other systems, provided non-local effects are included.Comment: 4 pages REVTEX, 5 figures using epsf macros. To appear in Physical Review E (Rapid Communications

Linear and nonlinear rheology of wormlike micelles

Several surfactant molecules self-assemble in solution to form long, cylindrical, flexible wormlike micelles. These micelles can be entangled with each other leading to viscoelastic phases. The rheological properties of such phases are very interesting and have been the subject of a large number of experimental and theoretical studies in recent years. We shall report on our recent work on the macrorheology, microrheology and nonlinear flow behaviour of dilute aqueous solutions of a surfactant CTAT (Cetyltrimethylammonium Tosilate). This system forms elongated micelles and exhibits strong viscoelasticity at low concentrations ($\sim$ 0.9 wt%) without the addition of electrolytes. Microrheology measurements of $G(\omega)$ have been done using diffusing wave spectroscopy which will be compared with the conventional frequency sweep measurements done using a cone and plate rheometer. The second part of the paper deals with the nonlinear rheology where the measured shear stress $\sigma$ is a nonmonotonic function of the shear rate $\dot{\gamma}$. In stress-controlled experiments, the shear stress shows a plateau for $\dot{\gamma}$ larger than some critical strain rate, similar to the earlier reports on CPyCl/NaSal system. Cates et al have proposed that the plateau is a signature of mechanical instability in the form of shear bands. We have carried out extensive experiments under controlled strain rate conditions, to study the time-dependence of shear stress. The measured time series of shear stress has been analysed in terms of correlation integrals and Lyapunov exponents to show unambiguously that the behaviour is typical of low dimensional dynamical systems.Comment: 15 pages, 10 eps figure

The Johnson-Segalman model with a diffusion term in Couette flow

We study the Johnson-Segalman (JS) model as a paradigm for some complex fluids which are observed to phase separate, or ``shear-band'' in flow. We analyze the behavior of this model in cylindrical Couette flow and demonstrate the history dependence inherent in the local JS model. We add a simple gradient term to the stress dynamics and demonstrate how this term breaks the degeneracy of the local model and prescribes a much smaller (discrete, rather than continuous) set of banded steady state solutions. We investigate some of the effects of the curvature of Couette flow on the observable steady state behavior and kinetics, and discuss some of the implications for metastability.Comment: 14 pp, to be published in Journal of Rheolog

Scale-free static and dynamical correlations in melts of monodisperse and Flory-distributed homopolymers: A review of recent bond-fluctuation model studies

It has been assumed until very recently that all long-range correlations are screened in three-dimensional melts of linear homopolymers on distances beyond the correlation length $\xi$ characterizing the decay of the density fluctuations. Summarizing simulation results obtained by means of a variant of the bond-fluctuation model with finite monomer excluded volume interactions and topology violating local and global Monte Carlo moves, we show that due to an interplay of the chain connectivity and the incompressibility constraint, both static and dynamical correlations arise on distances $r \gg \xi$. These correlations are scale-free and, surprisingly, do not depend explicitly on the compressibility of the solution. Both monodisperse and (essentially) Flory-distributed equilibrium polymers are considered.Comment: 60 pages, 49 figure

Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids

In this review, we describe and analyze a mesoscale simulation method for fluid flow, which was introduced by Malevanets and Kapral in 1999, and is now called multi-particle collision dynamics (MPC) or stochastic rotation dynamics (SRD). The method consists of alternating streaming and collision steps in an ensemble of point particles. The multi-particle collisions are performed by grouping particles in collision cells, and mass, momentum, and energy are locally conserved. This simulation technique captures both full hydrodynamic interactions and thermal fluctuations. The first part of the review begins with a description of several widely used MPC algorithms and then discusses important features of the original SRD algorithm and frequently used variations. Two complementary approaches for deriving the hydrodynamic equations and evaluating the transport coefficients are reviewed. It is then shown how MPC algorithms can be generalized to model non-ideal fluids, and binary mixtures with a consolute point. The importance of angular-momentum conservation for systems like phase-separated liquids with different viscosities is discussed. The second part of the review describes a number of recent applications of MPC algorithms to study colloid and polymer dynamics, the behavior of vesicles and cells in hydrodynamic flows, and the dynamics of viscoelastic fluids

Lattice Boltzmann simulations of soft matter systems

This article concerns numerical simulations of the dynamics of particles immersed in a continuum solvent. As prototypical systems, we consider colloidal dispersions of spherical particles and solutions of uncharged polymers. After a brief explanation of the concept of hydrodynamic interactions, we give a general overview over the various simulation methods that have been developed to cope with the resulting computational problems. We then focus on the approach we have developed, which couples a system of particles to a lattice Boltzmann model representing the solvent degrees of freedom. The standard D3Q19 lattice Boltzmann model is derived and explained in depth, followed by a detailed discussion of complementary methods for the coupling of solvent and solute. Colloidal dispersions are best described in terms of extended particles with appropriate boundary conditions at the surfaces, while particles with internal degrees of freedom are easier to simulate as an arrangement of mass points with frictional coupling to the solvent. In both cases, particular care has been taken to simulate thermal fluctuations in a consistent way. The usefulness of this methodology is illustrated by studies from our own research, where the dynamics of colloidal and polymeric systems has been investigated in both equilibrium and nonequilibrium situations.Comment: Review article, submitted to Advances in Polymer Science. 16 figures, 76 page

Recent experimental probes of shear banding

Recent experimental techniques used to investigate shear banding are reviewed. After recalling the rheological signature of shear-banded flows, we summarize the various tools for measuring locally the microstructure and the velocity field under shear. Local velocity measurements using dynamic light scattering and ultrasound are emphasized. A few results are extracted from current works to illustrate open questions and directions for future research.Comment: Review paper, 23 pages, 11 figures, 204 reference

Scaling laws for polymers in dissipative particle dynamics

A new technique, dissipative particle dynamics (DPD), appears promising as a means of studying the dynamical behaviour of polymers. Real polymers are known to obey a number of scaling laws, and a simulation method should reproduce these if it is to be relied on. The present work is a study of the properties of polymers, as simulated by DPD. Two cases are of interest: polymers in dilute solution, and polymers in the melt. The polymer in a good solvent shows satisfactory agreement with scaling and Kirkwood theory, and the polymer melt is in excellent agreement with the predictions of Rouse theory