Velocity distribution of a uniformly heated hard sphere granular gas

This paper presents a molecular dynamics simulation of an inelastic gas system under the influence of a uniform thermostat. We investigate the nonequilibrium properties of a granular gas where collisions between molecules are characterized by a coefficient of restitution less than unity. The simulation employs an event-driven algorithm to efficiently propagate the system in time, tracking molecular positions and velocities. A thermostat mechanism is incorporated to maintain the system's temperature by applying Gaussian white noise to the molecular velocities. The system's kinetic energy evolves towards a nonequilibrium steady state, with the initial dynamics governed by the interplay between energy input from the thermostat and energy dissipation through inelastic collisions. This steady state emerges when the energy gain from the thermostat balances the energy loss due to inelastic collisions. We calculate the coefficients of the Sonine polynomial expansion of the velocity distribution function to demonstrate the system's departure from Maxwell-Boltzmann statistics in the steady state.