Cosmological constraints on neutrino plus axion hot dark matter

We use observations of the cosmological large- scale structure to derive limits on two- component hot dark matter consisting of mass- degenerate neutrinos and hadronic axions, both components having velocity dispersions corresponding to their respective decoupling temperatures. We restrict the data samples to the safely linear regime, in particular excluding the Lyman-alpha forest. Using standard Bayesian inference techniques we derive credible regions in the two- parameter space of m(a) and Sigma m(nu) Marginalizing over Sigma m(nu) provides ma < 1.2 eV ( 95% C. L.). In the absence of axions the same data and methods give Sigma m. < 0.65 eV ( 95% C. L.). We also derive limits on ma for a range of axion - pion couplings up to one order of magnitude larger or smaller than the hadronic value.