Leading-order anisotropic hydrodynamics for central collisions

We use leading-order anisotropic hydrodynamics to study an azimuthally-symmetric boost-invariant quark-gluon plasma. We impose a realistic lattice-based equation of state and perform self-consistent anisotropic freeze-out to hadronic degrees of freedom. We then compare our results for the full spatiotemporal evolution of the quark-gluon plasma and its subsequent freeze-out to results obtained using 1+1d Israel-Stewart second-order viscous hydrodynamics. We find that for small shear viscosities, 4 pi eta/s ~ 1, the two methods agree well for nucleus-nucleus collisions, however, for large shear viscosity to entropy density ratios or proton-nucleus collisions we find important corrections to the Israel-Stewart results for the final particle spectra and the total number of charged particles. Finally, we demonstrate that the total number of charged particles produced is a monotonically increasing function of 4 pi eta/s in Israel-Stewart viscous hydrodynamics whereas in anisotropic hydrodynamics it has a maximum at 4 pi eta/s ~ 10. For all 4 pi eta/s > 0, we find that for Pb-Pb collisions Israel-Stewart viscous hydrodynamics predicts more dissipative particle production than anisotropic hydrodynamics.