Nonlinear diffusion in relativistic kinetic theory

A nonlinear Lorentz invariant kinetic diffusion equation is introduced, which is consistent with the conservation laws of particles number, energy and momentum. The equilibrium solution converges to the Maxwellian density in the Newtonian limit, but it is not given by the J\"uttner distribution commonly employed in relativistic kinetic theory. The nonlinear kinetic diffusion equation on a general Lorentzian manifold is consistent with the contracted Bianchi identities and therefore can be coupled to the Einstein equations of general relativity.