Abstract: We derive lower bounds on the density of sources of ultra-high energy cosmic rays from the lack of significant clustering in the arrival directions of the highest energy events detected at the Pierre Auger Observatory. The density of uniformly distributed sources of equal intrinsic intensity was found to be larger than similar to (0.06 – 5) x 10(-4) Mpc(-3) at 95% CL, depending on the magnitude of the magnetic defections. Similar bounds, in the range (0.2 – 7) x 10(-4) Mpc(-3), were obtained for sources following the local matter distribution.

Abstract: A new class of modified theories of gravity, consisting of the superposition of the metric Einstein-Hilbert Lagrangian with an f(R) term constructed a la Palatini was proposed recently. The dynamically equivalent scalar-tensor representation of the model was also formulated, and it was shown that even if the scalar field is very light, the theory passes the Solar System observational constraints. Therefore the model predicts the existence of a long-range scalar field, modifying the cosmological and galactic dynamics. An explicit model that passes the local tests and leads to cosmic acceleration was also obtained. In the present work, it is shown that the theory can be also formulated in terms of the quantity X equivalent to kappa T-2 + R, where T and R are the traces of the stress-energy and Ricci tensors, respectively. The variable X represents the deviation with respect to the field equation trace of general relativity. The cosmological applications of this hybrid metric-Palatini gravitational theory are also explored, and cosmological solutions coming from the scalar-tensor representation of f(X)-gravity are presented. Criteria to obtain cosmic acceleration are discussed and the field equations are analyzed as a dynamical system. Several classes of dynamical cosmological solutions, depending on the functional form of the effective scalar field potential, describing both accelerating and decelerating Universes are explicitly obtained. Furthermore, the cosmological perturbation equations are derived and applied to uncover the nature of the propagating scalar degree of freedom and the signatures these models predict in the large-scale structure.

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.