Comparison of hydrodynamic and "hybrid" hydrodynamics+transport calculations to heavy-ion data inevitably requires the conversion of the fluid to particles. For dissipative fluids the conversion is ambiguous without additional theory input complementing hydrodynamics. We obtain self-consistent shear viscous phase space corrections from linearized Boltzmann transport theory for a gas of hadrons. These corrections depend on the particle species, and incorporating them in Cooper-Frye freezeout affects identified particle observables. For example, with additive quark model cross sections, proton elliptic flow is larger than pion elliptic flow at moderately high pT in Au + Au collisions at RHIC. This is in contrast to Cooper-Frye freezeout with the commonly used "democratic Grad" ansatz that assumes no species dependence. Various analytic and numerical results are also presented for massless and massive two-component mixtures to better elucidate how species dependence arises. For convenient inclusion in pure hydrodynamic and hybrid calculations, Appendix G contains self-consistent viscous corrections for each species both in tabulated and parameterized form.

The quantitative extraction of quark-gluon plasma (QGP) properties from heavy-ion data, such as its specific shear viscosity η/s, typically requires comparison to viscous hydrodynamic or "hybrid" hydrodynamics+transport simulations. In either case, one has to convert the fluid to hadrons, yet without additional theory input the conversion is ambiguous for dissipative fluids. Here, shear viscous phase-space corrections calculated using linearized transport theory are applied in CooperFrye freezeout to quantify the effects on anisotropic flow coefficients vn(pT ) at both RHIC and LHC energies. Expanding upon our previous flow harmonics studies [1,2], we calculate pion and proton v2(pT ), v4(pT ), and v6(pT ), but here we incorporate a hadron gas that is chemically frozen below a temperature of 175 MeV and use hypersurfaces from realistic viscous hydrodynamic simulations. For additive quark model cross sections and relative phase-space corrections with p 3/2 momentum dependence rather than the quadratic Grad form, we find at moderately high transverse momentum noticeably higher v4(pT ) and v6(pT ) for protons than for pions. In addition, the value of η/s deduced from elliptic flow data differs by nearly 50% from the value extracted using the naive "democratic Grad" form of freeze-out distributions. To facilitate the use of the self-consistent viscous corrections calculated here in hydrodynamic and hybrid calculations, we also present convenient parameterizations of the corrections for the various hadron species (cf. Table I).

Comparing hydrodynamic simulations to heavy-ion data inevitably requires the conversion of the fluid to particles. This conversion, typically done in the Cooper-Frye formalism, is ambiguous for viscous fluids. We compute self-consistent phase space corrections by solving the linearized Boltzmann equation and contrast the solutions to those obtained using the ad-hoc "democratic Grad" ansatz typically employed in the literature where coefficients are independent of particle dynamics. Solutions are calculated analytically for a massless gas and numerically for both a pion-nucleon gas and for the general case of a hadron resonance gas. We find that the momentum dependence of the corrections in all systems investigated is best fit by a power close to 3/2 rather than the typically used quadratic ansatz. The effects on harmonic flow coefficients v2 and v4 are substantial, and should be taken into account when extracting medium properties from experimental data.