The Langevin equation m·dv/dt = −γv + η(t) describes a particle buffeted by random thermal forces η(t).
In the overdamped limit, mean squared displacement ⟨r²⟩ = 2dDt where D = k_BT/γ (Einstein relation).
The displacement histogram should approach a Gaussian. Increase temperature for faster, wilder motion.