are not weak, it merely holds that momentum relaxing collisions are weak while momentum-conserving collisions are not. A formulation in terms of non-local transport coefficients allow for a microscopic description of hydrodynamic flow pattern and goes beyond the usual description in terms of the linear Navier-Stokes equation. The latter corresponds to the leading gradient expansion of the theory. In addition, the inclusion of dynamical phenomena - here expressed in terms of the dependency on the time difference $t-t'$ between force and current – allows to determine the system’s collective modes.

In this paper we develop the theory of non-local transport in Dirac systems at charge neutrality in the collision-dominated hydrodynamic regime and find closed expressions for the frequency and wave-vector-dependent, charge and thermal conductivities as well as the non-local viscosity. Remarkably, the calculations of this paper are exact in the limit of a small graphene fine structure constant $\alpha$ in the regime of linear response. This is made possible by the peculiar $\propto\left|m\right|$ dependence of the scattering rates of collinear zero modes in higher angular momentum channels $m>2$ - a behavior that was shown to lead to a super-diffusive Lévy-flight-like phase space dynamics in an earlier work [@Kiselev2019b]. Collinear zero 