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 use the formalism of Gould and Raffelt [1] to compute the dimensionless thermal conduction coefficients for scattering of dark matter particles with standard model nucleons via cross-sections that depend on the relative velocity or momentum exchanged between particles. Motivated by models invoked to reconcile various recent results in direct detection, we explicitly compute the conduction coefficients alpha and kappa for cross-sections that go as v(rel)(2), v(rel)(4), v(rel)(-2), q(2), q(4) and q(-2), where v(rel) is the relative DM-nucleus velocity and q is the momentum transferred in the collision. We find that a v(rel)(-2) depend ence can significantly enhance energy transport from the inner solar core to the outer core. The same can true for any q-dependent coupling, if the dark matter mass lies within some specific range for each coupling. This effect can complement direct searches for dark matter; combining these results with state-of-the-art solar simulations should greatly increase sensitivity to certain DM models. It also seems possible that the so-called Solar Abundance Problem could be resolved by enhanced energy transport in the solar core due to such velocity-or momentum-dependent scatterings.

Abstract: Sterile neutrinos with mass in the eV-scale and large mixings of order theta(0) similar or equal to 0.1 could explain some anomalies found in short-baseline neutrino oscillation data. Here, we revisit a neutrino portal scenario in which eV-scale sterile neutrinos have self-interactions via a new gauge vector boson phi. Their production in the early Universe via mixing with active neutrinos can be suppressed by the induced effective potential in the sterile sector. We study how different cosmological observations can constrain this model, in terms of the mass of the new gauge boson, M-phi, and its coupling to sterile neutrinos, g(s). Then, we explore how to probe part of the allowed parameter space of this particular model with future observations of the diffuse supernova neutrino background by the Hyper-Kamiokande and DUNE detectors. For M-phi similar to 5 – 10 keV and g(s) similar to 10-(4) – 10(-2), as allowed by cosmological constraints, we find that interactions of diffuse supernova neutrinos with relic sterile neutrinos on their way to the Earth would result in significant dips in the neutrino spectrum which would produce unique features in the event spectra observed in these detectors.

Abstract: We show that an invariant an universal length scale can be consistently introduced in a generally covariant theory through the gravitational sector using the Palatini approach. The resulting theory is able to capture different aspects of quantum gravity phenomenology in a single framework. In particular, it is found that in this theory field excitations propagating with different energy-densities perceive different background metrics, which is a fundamental characteristic of the DSR and Rainbow Gravity approaches. We illustrate these properties with a particular gravitational model and explicitly show how the soccer ball problem is avoided in this framework. The isotropic and anisotropic cosmologies of this model also avoid the big bang singularity by means of a big bounce.

Abstract: The metric-affine bumblebee model in the presence of fermionic matter minimally coupled to the connection is studied. We show that the model admits an Einstein frame representation in which the matter sector is described by a non-minimal Dirac action without any analogy in the literature. Such non-minimal terms involve unconventional couplings between the bumblebee and the fermion field. We then rewrite the quadratic fermion action in the Einstein frame in the basis of 16 Dirac matrices in order to identify the coefficients for Lorentz/CPT violation in all orders of the non-minimal coupling xi. The exact result for the fermionic determinant in the Einstein frame, including all orders in xi, is also provided. We demonstrate that the axial contributions are at least of second order in the perturbative expansion of xi. Furthermore, we compute the one-loop effective potential within the weak field approximation.