Pairing of Lyapunov Exponents for a Hard-Sphere Gas under Shear in the Thermodynamic Limit

The phase space trajectories of many body systems charateristic of simple fluids are highly unstable. We quantify this instability by a set of Lyapunov exponents, which are the rates of exponential divergence, or convergence, of initial (infinitesimal) perturbations along carefully selected directions in phase space. It is demonstrated that the perturbation associated with the maximum Lyapunov exponent is localized in space. This localization persists in the large-particle li...

A dilute gas of particles with short range interactions is considered in a shearing stationary state. A Gaussian thermostat keeps the total kinetic energy constant. For infinitely many particles it is shown that the thermostat becomes a friction force with constant friction coefficient. For finite number of particles N, the fluctuations around this constant are of order one over the square root of N, and distributed approximately Gaussian with deviations for large fluctuation...

The fundamental question addressed in this paper is whether the partial Chapman--Enskog expansion $P_{xy}=-\sum_{k=0}^\infty \eta_k ({\partial u_x}/{\partial y})^{2k+1}$ of the shear stress converges or not for a gas of inelastic hard spheres. By using a simple kinetic model it is shown that, in contrast to the elastic case, the above series does converge, the radius of convergence increasing with inelasticity. It is argued that this paradoxical conclusion is not an artifact ...

We compute the Lyapunov spectrum and the Kolmogorov-Sinai entropy for a moving particle placed in a dilute, random array of hard disk or hard sphere scatterers - i.e. the dilute Lorentz gas model. This is carried out in two ways: First we use simple kinetic theory arguments to compute the Lyapunov spectrum for both two and three dimensional systems. In order to provide a method that can easily be generalized to non-uniform systems we then use a method based upon extensions of...

We use molecular dynamics simulations to compute the Lyapunov spectra of many-particle systems resembling simple fluids in thermal equilibrium and in non-equilibrium stationary states. Here we review some of the most interesting results and point to open questions.

The velocity fluctuations present in macroscopically homogeneous suspensions of neutrally buoyant, non-Brownian spheres undergoing simple shear flow, and their dependence on the microstructure developed by the suspensions, are investigated in the limit of vanishingly small Reynolds numbers using Stokesian dynamics simulations. We show that, in the dilute limit, the standard deviation of the velocity fluctuations is proportional to the volume fraction, in both the transverse a...

The evolution of a gas can be described by different models depending on the observation scale. A natural question, raised by Hilbert in his sixth problem, is whether these models provide consistent predictions. In particular, for rarefied gases, it is expected that continuum laws of kinetic theory can be obtained directly from molecular dynamics governed by the fundamental principles of mechanics. In the case of hard sphere gases, Lanford showed that the Boltzmann equation e...

The dynamics of a system consisting of many spherical hard particles can be described as a single point particle moving in a high-dimensional space with fixed hypercylindrical scatterers with specific orientations and positions. In this paper, the similarities in the Lyapunov exponents are investigated between systems of many particles and high-dimensional billiards with cylindrical scatterers which have isotropically distributed orientations and homogeneously distributed pos...

Via event driven molecular dynamics simulations and experiments, we study the packing fraction and shear-rate dependence of single particle fluctuations and dynamic correlations in hard sphere glasses under shear. At packing fractions above the glass transition, correlations increase as shear rate decreases: the exponential tail in the distribution of single particle jumps broadens and dynamic four-point correlations increase. Interestingly, however, upon decreasing the packi...

We study the Lyapunov exponents for a moving, charged particle in a two-dimensional Lorentz gas with randomly placed, non-overlapping hard disk scatterers placed in a thermostatted electric field, $\vec{E}$. The low density values of the Lyapunov exponents have been calculated with the use of an extended Lorentz-Boltzmann equation. In this paper we develop a method to extend these results to higher density, using the BBGKY hierarchy equations and extending them to include the...