Random-walk derivation of stellar relaxation ($t_{\mathrm{rel}}\sim 0.1(N/\ln N)t_{\mathrm{cross}}$)

Small-angle gravitational encounters occur at rate $2\pi bdbnu$ and supply kicks $2Gm/(bu)$. Independent kick variances add to $d⟨\Delta v^2⟩/dt=8\pi G^2{m}^{2}n\ln \Lambda /u$. Setting the accumulated variance to $u^2$ and using virial scaling gives a stellar relaxation time of order $0.1N/\ln N$ stellar crossing times, with the coefficient dependent on profile and crossing convention.