Constructing higher-order hydrodynamics: The third order

Hydrodynamics can be formulated as the gradient expansion of conserved currents in terms of the fundamental fields describing the near-equilibrium fluid flow. In the relativistic case, the Navier-Stokes equations follow from the conservation of the stress-energy tensor to first order in derivatives. In this paper, we go beyond the presently understood second-order hydrodynamics and discuss the systematisation of obtaining the hydrodynamic expansion to an arbitrarily high order. As an example of the algorithm that we present, we fully classify the gradient expansion at third order for neutral fluids in four dimensions, thus finding the most general next-to-leading-order corrections to the relativistic Navier-Stokes equations in curved space-time. In doing so, we list $20$ new transport coefficient candidates in the conformal and $68$ in the non-conformal case. As we do not consider any constraints that could potentially arise from the local entropy current analysis, this is the maximal possible set of neutral third-order transport coefficients. To investigate the physical implications of these new transport coefficients, we obtain the third-order corrections to the linear dispersion relations that describe the propagation of diffusion and sound waves in relativistic fluids. We also compute the corrections to the scalar (spin-$2$) two-point correlation function of the third-order stress-energy tensor. Furthermore, as an example of a non-linear hydrodynamic flow, we calculate the third-order corrections to the energy density of a boost-invariant Bjorken flow. Finally, we apply our field theoretic results to the $\mathcal{N}=4$ supersymmetric Yang-Mills fluid at infinite 't Hooft coupling and infinite number of colours to find the values of five new linear combinations of the conformal transport coefficients.