Statistical mechanics far from equilibrium: prediction and test for a sheared system

We report the complete statistical treatment of a system of particles interacting via Newtonian forces in continuous boundary-driven flow, far from equilibrium. By numerically time-stepping the force-balance equations of a model fluid we measure occupancies and transition rates in simulation. The high-shear-rate simulation data verify the invariant quantities predicted by our statistical theory, thus demonstrating that a class of non-equilibrium steady states of matter, namely sheared complex fluids, is amenable to statistical treatment from first principles.Comment: 4 pages plus a 3-page pdf supplemen

Coarse-grained simulations of flow-induced nucleation in semi-crystalline polymers

We perform kinetic Monte Carlo simulations of flow-induced nucleation in polymer melts with an algorithm that is tractable even at low undercooling. The configuration of the non-crystallized chains under flow is computed with a recent non-linear tube model. Our simulations predict both enhanced nucleation and the growth of shish-like elongated nuclei for sufficiently fast flows. The simulations predict several experimental phenomena and theoretically justify a previously empirical result for the flow-enhanced nucleation rate. The simulations are highly pertinent to both the fundamental understanding and process modeling of flow-induced crystallization in polymer melts.Comment: 17 pages, 6 eps figure

Validation of the Jarzynski relation for a system with strong thermal coupling: an isothermal ideal gas model

We revisit the paradigm of an ideal gas under isothermal conditions. A moving piston performs work on an ideal gas in a container that is strongly coupled to a heat reservoir. The thermal coupling is modeled by stochastic scattering at the boundaries. In contrast to recent studies of an adiabatic ideal gas with a piston [R.C. Lua and A.Y. Grosberg, J. Phys. Chem. B 109, 6805 (2005); I. Bena et al., Europhys. Lett. 71, 879 (2005)], the container and piston stay in contact with the heat bath during the work process. Under this condition the heat reservoir as well as the system depend on the work parameter lambda and microscopic reversibility is broken for a moving piston. Our model is thus not included in the class of systems for which the nonequilibrium work theorem has been derived rigorously either by Hamiltonian [C. Jarzynski, J. Stat. Mech. (2004) P09005] or stochastic methods [G.E. Crooks, J. Stat. Phys. 90, 1481 (1998)]. Nevertheless the validity of the nonequilibrium work theorem is confirmed both numerically for a wide range of parameter values and analytically in the limit of a very fast moving piston, i.e., in the far nonequilibrium regime

Loss of solutions in shear banding fluids in shear banding fluids driven by second normal stress differences

Edge fracture occurs frequently in non-Newtonian fluids. A similar instability has often been reported at the free surface of fluids undergoing shear banding, and leads to expulsion of the sample. In this paper the distortion of the free surface of such a shear banding fluid is calculated by balancing the surface tension against the second normal stresses induced in the two shear bands, and simultaneously requiring a continuous and smooth meniscus. We show that wormlike micelles typically retain meniscus integrity when shear banding, but in some cases can lose integrity for a range of average applied shear rates during which one expects shear banding. This meniscus fracture would lead to ejection of the sample as the shear banding region is swept through. We further show that entangled polymer solutions are expected to display a propensity for fracture, because of their much larger second normal stresses. These calculations are consistent with available data in the literature. We also estimate the meniscus distortion of a three band configuration, as has been observed in some wormlike micellar solutions in a cone and plate geometry.Comment: 23 pages, to be published in Journal of Rheolog

Inhomogeneous High Frequency Expansion-Free Gravitational Waves

We describe a natural inhomogeneous generalization of high frequency plane gravitational waves. The waves are high frequency waves of the Kundt type whose null propagation direction in space-time has vanishing expansion, twist and shear but is not covariantly constant. The introduction of a cosmological constant is discussed in some detail and a comparison is made with high frequency gravity waves having wave fronts homeomorphic to 2-spheres.Comment: 18 pages, Latex file, accepted for publication in Physical Review

Nonequilibrium statistical mechanics of shear flow: invariant quantities and current relations

In modeling nonequilibrium systems one usually starts with a definition of the microscopic dynamics, e.g., in terms of transition rates, and then derives the resulting macroscopic behavior. We address the inverse question for a class of steady state systems, namely complex fluids under continuous shear flow: how does an externally imposed shear current affect the microscopic dynamics of the fluid? The answer can be formulated in the form of invariant quantities, exact relations for the transition rates in the nonequilibrium steady state, as discussed in a recent letter [A. Baule and R. M. L. Evans, Phys. Rev. Lett. 101, 240601 (2008)]. Here, we present a more pedagogical account of the invariant quantities and the theory underlying them, known as the nonequilibrium counterpart to detailed balance (NCDB). Furthermore, we investigate the relationship between the transition rates and the shear current in the steady state. We show that a fluctuation relation of the Gallavotti-Cohen type holds for systems satisfying NCDB.Comment: 24 pages, 11 figure

Lattice Boltzmann Simulations of Liquid Crystal Hydrodynamics

We describe a lattice Boltzmann algorithm to simulate liquid crystal hydrodynamics. The equations of motion are written in terms of a tensor order parameter. This allows both the isotropic and the nematic phases to be considered. Backflow effects and the hydrodynamics of topological defects are naturally included in the simulations, as are viscoelastic properties such as shear-thinning and shear-banding.Comment: 14 pages, 5 figures, Revte

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

Lateral phase separation in mixtures of lipids and cholesterol

In an effort to understand "rafts" in biological membranes, we propose phenomenological models for saturated and unsaturated lipid mixtures, and lipid-cholesterol mixtures. We consider simple couplings between the local composition and internal membrane structure, and their influence on transitions between liquid and gel membrane phases. Assuming that the gel transition temperature of the saturated lipid is shifted by the presence of the unsaturated lipid, and that cholesterol acts as an external field on the chain melting transition, a variety of phase diagrams are obtained. The phase diagrams for binary mixtures of saturated/unsaturated lipids and lipid/cholesterol are in semi-quantitative agreement with the experiments. Our results also apply to regions in the ternary phase diagram of lipid/lipid/cholesterol systems

Assessment of interatomic potentials for atomistic analysis of static and dynamic properties of screw dislocations in W

Screw dislocations in bcc metals display non-planar cores at zero temperature which result in high lattice friction and thermally activated strain rate behavior. In bcc W, electronic structure molecular statics calculations reveal a compact, non-degenerate core with an associated Peierls stress between 1.7 and 2.8 GPa. However, a full picture of the dynamic behavior of dislocations can only be gained by using more efficient atomistic simulations based on semiempirical interatomic potentials. In this paper we assess the suitability of five different potentials in terms of static properties relevant to screw dislocations in pure W. As well, we perform molecular dynamics simulations of stress-assisted glide using all five potentials to study the dynamic behavior of screw dislocations under shear stress. Dislocations are seen to display thermally-activated motion in most of the applied stress range, with a gradual transition to a viscous damping regime at high stresses. We find that one potential predicts a core transformation from compact to dissociated at finite temperature that affects the energetics of kink-pair production and impacts the mechanism of motion. We conclude that a modified embedded-atom potential achieves the best compromise in terms of static and dynamic screw dislocation properties, although at an expense of about ten-fold compared to central potentials