Improved hot dark matter bound on the QCD axion

We obtain a reliable cosmological bound on the axion mass ma by (1) deriving the production rate directly from pion-pion scattering data, which overcomes the breakdown of chiral perturbation theory and results in & SIM;30% differences from previous estimates; (2) including momentum dependence in the Boltzmann equations for axion-pion scatterings, which enhances the relic abundance by & SIM;40%. Using present cosmological datasets we obtain ma < 0.24 eV, at 95% C.L. We also constrain the sum of neutrino masses, P m & nu; < 0.14 eV at 95% C.L., in the presence of relic axions and neutrinos. Finally, we show that reliable nonperturbative calculations above the QCD crossover are needed to exploit the reach of upcoming cosmological surveys for axion detection.