Breakdown of the Sonine expansion for the velocity distribution of granular gases

The velocity distribution of a granular gas is analyzed in terms of the Sonine polynomials expansion. We derive an analytical expression for the third Sonine coefficient a(3). In contrast to frequently used assumptions this coefficient is of the same order of magnitude as the second Sonine coefficient a(2). For small inelasticity the theoretical result is in good agreement with numerical simulations. The next-order Sonine coefficients a(4), a(5) and a(6) are determined numerically. While these coefficients are negligible for small dissipation, their magnitude grows rapidly with increasing inelasticity for 0 < epsilon less than or similar to 0.6. We conclude that this behavior of the Sonine coefficients manifests the breakdown of the Sonine polynomial expansion caused by the increasing impact of the overpopulated high-energy tail of the distribution function.