  In this paper, we develop and characterize the fully dissipative Lattice
Boltzmann method for ultra-relativistic fluids in two dimensions using three
equilibrium distribution functions: Maxwell-J\"uttner, Fermi-Dirac and
Bose-Einstein. Our results stem from the expansion of these distribution
functions up to fifth order in relativistic polynomials. We also obtain new
Gaussian quadratures for square lattices that preserve the spatial resolution.
Our models are validated with the Riemann problem and the limitations of lower
order expansions to calculate higher order moments are shown. The kinematic
viscosity and the thermal conductivity are numerically obtained using the
Taylor-Green vortex and the Fourier flow respectively and these transport
coefficients are compared with the theoretical prediction from Grad's theory.
In order to compare different expansion orders, we analyze the temperature and
heat flux fields on the time evolution of a hot spot.
