The Fokker-Planck equation for multivibron solitons interacting with lattice vibrations in a molecular chain has been derived by means of the nonequilibrium statistical operator method. It was shown that a soliton undergoes diffusive motion characterized by two substantially different diffusion coefficients. The first one corresponds to the ordinary (Einsteinian or dissipative) diffusion and characterizes the soliton Brownian motion, while the second one corresponds to the anomalous diffusion connected with frictionless displacement of the soliton center of mass coordinate due to the interaction with phonons. Both processes are the consequence of the Cherenkov-like radiation of phonon quanta arising when soliton velocity approaches the phase speed of sound.