Nonlinear multimode optical platforms can induce effective quantum interactions between photons, resulting in the thermalization of single- and two-photon states. This thermalization is manifested by the emergence of Rayleigh-Jeans and Boltzmann statistical distributions.
