Diffusion in a periodic Lorentz gas

Abstract

We use a constant "driving force"Fd together with a Gaussian thermostatting "constraint force"Fd to simulate a nonequilibrium steady-state current (particle velocity) in a periodic, two-dimensional, classical Lorentz gas. The ratio of the average particle velocity to the driving force (field strength) is the Lorentz-gas conductivity. A regular "Galton-board" lattice of fixed particles is arranged in a dense triangular-lattice structure. The moving scatterer particle travels through the lattice at constant kinetic energy, making elastic hard-disk collisions with the fixed particles. At low field strengths the nonequilibrium conductivity is statistically indistinguishable from the equilibrium Green-Kubo estimate of Machta and Zwanzig. The low-field conductivity varies smoothly, but in a complicated way, with field strength. For moderate fields the conductivity generally decreases nearly linearly with field, but is nearly discontinuous at certain values where interesting stable cycles of collisions occur. As the field is increased, the phase-space probability density drops in apparent fractal dimensionality from 3 to 1. We compare the nonlinear conductivity with similar zero-density results from the two-particle Boltzmann equation. We also tabulate the variation of the kinetic pressure as a function of the field strength, © 1987 Plenum Publishing Corporation.