Thermalization in classical systems with discrete phase space

We study the emergence of statistical mechanics in isolated classical systems with local interactions and discrete phase spaces. We establish that thermalization in such systems does not require global ergodicity; instead, it arises from effective local ergodicity, where dynamics in a subsystem may appear pseudorandom. To corroborate that, we analyze the spectrum of the unitary evolution operator and propose an ansatz to describe statistical properties of local observables expanded in the eigenfunction basis - the classical counterpart of the Eigenstate Thermalization Hypothesis. Our framework provides a unified perspective on thermalization in classical and quantum systems with discrete spectra.