Abstract: This paper is a continuation of an earlier work (arXiv:1210.7928) which computed analytically the tree-level annihilation rates of a collection of non-relativistic neutralino and chargino two-particle states in the general MSSM. Here we extend the results by providing the next-to-next-to-leading order corrections to the rates in the non-relativistic expansion in momenta and mass differences, which include leading P-wave effects, in analytic form. The results are a necessary input for the calculation of the Sommerfeld-enhanced dark matter annihilation rates including short-distance corrections at next-to-next-to-leading order in the non-relativistic expansion in the general MSSM with neutralino LSP.

Abstract: We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations which incorporate cold dark matter (CDM) and neutrinos as independent particle species. Our set of simulations explore the properties of neutrinos in a reference Lambda CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass 10(11) – 10(15) h(-1) M-circle dot, over a redshift range z = 0 – 2. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula, once the neutrino contribution to the total matter is removed. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than similar to 10(13.5) h(-1) M-circle dot. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. The results we obtain are numerically converged in terms of neutrino profiles at the 10% level for scales above similar to 200 h(-1) kpc at z = 0, and are stable with respect to box-size and starting redshift of the simulation. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.

Abstract: We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts measured from Sloan Digital Sky Survey III Data Release 8. The galaxies have photometric redshifts between z = 0.45 and z = 0.65 and cover 10,000 deg(2), thus probing a volume of 3 h(-3) Gpc(3) and enabling tight constraints to be derived on the amount of dark matter in the form of massive neutrinos. A new bound on the sum of neutrino masses Sigma m nu < 0.27 eV, at the 95% confidence level (CL), is obtained after combining our sample of galaxies, which we call “CMASS,” with Wilkinson Microwave Anisotropy Probe (WMAP) seven-year cosmic microwave background data and the most recent measurement of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a conservative multipole range of 30 < l < 200 in order to minimize nonlinearities, and a free bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy bias model using mock catalogs and find that this model causes a small (similar to 1 sigma-1.5 sigma) bias in Omega(DM)h(2). For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing additional, shot-noise-like free parameters. In this conservative case, the bounds are significantly weakened, e. g., Sigma m(nu) < 0.38 eV (95% CL) for WMAP+HST+CMASS (l(max) = 200). We also study the dependence of the neutrino bound on the multipole range (l(max) = 150 versus l(max) = 200) and on which combination of data sets is included as a prior. The addition of supernova and/or baryon acoustic oscillation data does not significantly improve the neutrino mass bound once the HST prior is included. A companion paper describes the construction of the angular power spectra in detail and derives constraints on a general cosmological model, including the dark energy equation of state w and the spatial curvature Omega(K), while a second companion paper presents a measurement of the scale of baryon acoustic oscillations from the same data set. All three works are based on the catalog by Ross et al.