No Arabic abstract
The analysis of the dynamics of tracer particles in a complex bath can provide valuable information about the microscopic behaviour of the bath. In this work, we study the dynamics of a forced tracer in a colloidal bath by means of Langevin dynamics simulations and a theory model within continuum mechanics. In the simulations, the bath is comprised by quasi-hard spheres with a volume fraction of 50% immersed in a featureless quiescent solvent, and the tracer is pulled with a constant small force (within the linear regime). The theoretical analysis is based on the Navier Stokes equation, where a term proportional to the velocity arises from coarse-graining the friction of the colloidal particles with the solvent. As a result, the final equation is similar to the Brinkman model, although the interpretation is different. A length scale appears in the model, 1/k_0, where the transverse momentum transport crosses over to friction with the solvent. The effective friction coefficient experienced by the tracer grows with the tracer size faster than the prediction from Stokes law. Additionally, the velocity profiles in the bath decay faster than in a Newtonian fluid. The comparison between simulations and theory points to a boundary condition of effective partial slip at the tracer surface. We also study the fluctuations in the tracer position, showing that it reaches diffusion at long times, with a subdiffusive regime at intermediate times. The diffusion coefficient, obtained from the long-time slope of the mean squared displacement, fulfills the Stokes-Einstein relation with the friction coefficient calculated from the steady tracer velocity, confirming the validity of the linear response formalism.
The transient response of model hard sphere glasses is examined during the application of steady rate start-up shear using Brownian Dynamics (BD) simulations, experimental rheology and confocal microscopy. With increasing strain the glass initially exhibits an almost linear elastic stress increase, a stress peak at the yield point and then reaches a constant steady state. The stress overshoot has a non-monotonic dependence with Peclet number, Pe, and volume fraction, {phi}, determined by the available free volume and a competition between structural relaxation and shear advection. Examination of the structural properties under shear revealed an increasing anisotropic radial distribution function, g(r), mostly in the velocity - gradient (xy) plane, which decreases after the stress peak with considerable anisotropy remaining in the steady-state. Low rates minimally distort the structure, while high rates show distortion with signatures of transient elongation. As a mechanism of storing energy, particles are trapped within a cage distorted more than Brownian relaxation allows, while at larger strains, stresses are relaxed as particles are forced out of the cage due to advection. Even in the steady state, intermediate super diffusion is observed at high rates and is a signature of the continuous breaking and reformation of cages under shear.
Normal dynamics in a quasi-one-dimensional channel of length L (toinfty) of N hard spheres are analyzed. The spheres are heterogeneous: each has a diffusion coefficient D that is drawn from a probability density function (PDF), W D^(-{gamma}), for small D, where 0leq{gamma}<1. The initial spheres density {rho} is non-uniform and scales with the distance (from the origin) l as, {rho} l^(-a), 0leqaleq1. An approximation for the N-particle PDF for this problem is derived. From this solution, scaling law analysis and numerical simulations, we show here that the mean square displacement for a particle in such a system obeys, <r^2>~t^(1-{gamma})/(2c-{gamma}), where c=1/(1+a). The PDF of the tagged particle is Gaussian in position. Generalizations of these results are considered.
We report results of dynamic light scattering measurements of the coherent intermediate scattering function (ISF) of glasses of hard spheres for several volume fractions and a range of scattering vectors around the primary maximum of the static structure factor. The ISF shows a clear crossover from an initial fast decay to a slower non-stationary decay. Ageing is quantified in several different ways. However, regardless of the method chosen, the perfect aged glass is approached in a power-law fashion. In particular, the coupling between the fast and slow decays, as measured by the degree of stretching of the ISF at the crossover, also decreases algebraically with waiting time. The non-stationarity of this coupling implies that even the fastest detectable processes are themselves non-stationary.
In this work we compare and characterize the behavior of Langevin and Dissipative Particle Dynamics (DPD) thermostats in a broad range of non-equilibrium simulations of polymeric systems. Polymer brushes in relative sliding motion, polymeric liquids in Poiseuille and Couette flows, and brush-melt interfaces are used as model systems to analyze the efficiency and limitations of different Langevin and DPD thermostat implementations. Widely used coarse-grained bead-spring models under good and poor solvent conditions are employed to assess the effects of the thermostats. We considered equilibrium, transient, and steady state examples for testing the ability of the thermostats to maintain constant temperature and to reproduce the underlying physical phenomena in non-equilibrium situations. The common practice of switching-off the Langevin thermostat in the flow direction is also critically revisited. The efficiency of different weight functions for the DPD thermostat is quantitatively analyzed as a function of the solvent quality and the non-equilibrium situation.
The nonlinear dynamics associated with sliding friction forms a broad interdisciplinary research field that involves complex dynamical processes and patterns covering a broad range of time and length scales. Progress in experimental techniques and computational resources has stimulated the development of more refined and accurate mathematical and numerical models, capable of capturing many of the essentially nonlinear phenomena involved in friction.