We study numerically time evolution in classical lattices with weak or
moderate nonlinearity which leads to interactions between linear modes. Our
results show that in a certain strength range a moderate nonlinearity generates
a dynamical thermalization process which drives the system to the quantum Gibbs
distribution of probabilities, or average oscillation amplitudes. The effective
dynamical temperature of the lattice varies from large positive to large
negative values depending on energy of initially excited modes. This quantum
Gibbs distribution is drastically different from usually expected energy
equipartition over linear modes corresponding to a regime of classical
thermalization. Possible experimental observations of this dynamical
thermalization are discussed for cold atoms in optical lattices, nonlinear
photonic lattices and optical fiber arrays.Comment: 15 pages, 12 figures. Small modifs., video abstract 107MB at
http://www.quantware.ups-tlse.fr/dima/video/gibbs2013.mp