Molecular dynamics simulation of collision operator eigenvalues

We simulate numerically the time evolution of 1000 interacting hard spheres in a finite box with periodic boundary conditions, and repeat the simulations many times for a selected random distribution of initial conditions. We use the resulting data to compute, directly, the smallest nonzero eigenvalue of the collision operator for this gas. We also give exact expressions for the transport coefficients and compare them to approximate expressions commonly seen in the literature.