Proton decay is a smoking gun signature of Grand Unified Theories (GUTs). Searches by Super-Kamiokande have resulted in stringent limits on the GUT symmetry breaking scale. The large-scale multipurpose neutrino experiments DUNE, Hyper-Kamiokande and JUNO will either discover proton decay or further push the symmetry breaking scale above $10^{16}$ GeV. Another possible observational consequence of GUTs is the formation of a cosmic string network produced during the breaking of the GUT to the Standard Model gauge group. The evolution of such a string network in the expanding Universe produces a stochastic background of gravitational waves which will be tested by a number of gravitational wave detectors over a wide frequency range. We demonstrate the non-trivial complementary between the observation of proton decay and gravitational waves produced from cosmic strings in determining $SO(10)$ GUT breaking chains. We show that such observations could exclude $SO(10)$ breaking via flipped $SU(5)\times U(1)$ or standard $SU(5)$, while breaking via a Pati-Salam intermediate symmetry, or standard $SU(5)\times U(1)$, may be favoured if a large separation of energy scales associated with proton decay and cosmic strings is indicated.

The experiment to detect the global epoch of reionization signature (EDGES) collaboration reported the detection of a line at 78MHz in the sky-averaged spectrum due to neutral hydrogen (HI) 21-cm hyperfine absorption of cosmic microwave background photons at $z\sim 17$. This requires that the spin temperature of HI be coupled to the kinetic temperature of the gas at this $z$ through the scattering of Ly-$\alpha$ photons emitted by massive stars. To explain the experimental result, star formation needs to be sufficiently efficient at $z\sim 17$ and this can be used to constrain models in which small-scale structure formation is suppressed (DMF models), either due to dark matter free-streaming or non-standard inflationary dynamics. We combine simulations of structure formation with a simple recipe for star formation to investigate whether these models emit enough Ly-$\alpha$ photons to reproduce the experimental signal for reasonable values of the star formation efficiency, $f_\star$. We find that a thermal warm dark matter (WDM) model with mass $m_\mathrm{WDM}\sim 4.3\,\mathrm{keV}$ is consistent with the timing of the signal for $f_\star\lesssim 2\%$. The exponential growth of structure around $z\sim 17$ in such a model naturally generates a sharp onset of the absorption. A warmer model with $m_\mathrm{WDM}\sim 3\,\mathrm{keV}$ requires a higher star formation efficiency, $f_\star\sim 6\%$, which is a factor of few above predictions of current star formation models and observations of satellites in the Milky Way. However, uncertainties in the process of star formation at these $z$ do not allow to derive strong constrains on such models using 21-cm absorption line. The onset of the 21-cm absorption is generally faster in DMF compared to cold dark matter (CDM), unless some process significantly suppresses star formation in halos with masses below $\sim 10^8\,h^{-1}\,\mathrm{M}_\odot$.

Thermal inflation models (which feature two inflationary stages) can display damped primordial curvature power spectra on small scales which generate damped matter fluctuations. For a reasonable choice of parameters, thermal inflation models naturally predict a suppression of the matter power spectrum on galactic and sub-galactic scales, mimicking the effect of warm or interacting dark matter. Matter power spectra in these models are also characterised by an excess of power (w.r.t. the standard $\Lambda$CDM power spectrum) just below the suppression scale. By running a suite of N-body simulations we investigate the non-linear growth of structure in models of thermal inflation. We measure the non-linear matter power spectrum and extract halo statistics, such as the halo mass function, and compare these quantities with those predicted in the standard $\Lambda$CDM model and in other models with damped matter fluctuations. We find that the thermal inflation models considered here produce measurable differences in the matter power spectrum from $\Lambda$CDM at redshifts $z>5$, while the halo mass functions are appreciably different at all redshifts. The halo mass function at $z=0$ for thermal inflation displays an enhancement of around $\sim 20\%$ w.r.t. $\Lambda$CDM and a damping at lower halo masses, with the position of the enhancement depending on the value of the free parameter in the model. The enhancement in the halo mass function (w.r.t. $\Lambda$CDM ) increases with redshift, reaching $\sim 40\%$ at $z=5$. We also study the accuracy of the analytical Press-Schechter approach, using different filters to smooth the density field, to predict halo statistics for thermal inflation. We find that the predictions with the smooth-$k$ filter agree with the simulation results over a wider range of halo masses than is the case with other filters commonly used in the literature.

We investigate the nonlinear evolution of structure in variants of the standard cosmological model which display damped density fluctuations relative to cold dark matter (e.g. in which cold dark matter is replaced by warm or interacting DM). Using N-body simulations, we address the question of how much information is retained from different scales in the initial linear power spectrum following the nonlinear growth of structure. We run a suite of N-body simulations with different initial linear matter power spectra to show that, once the system undergoes nonlinear evolution, the shape of the linear power spectrum at high wavenumbers does not affect the non-linear power spectrum, while it still matters for the halo mass function. Indeed, we find that linear power spectra which differ from one another only at wavenumbers larger than their half-mode wavenumber give rise to (almost) identical nonlinear power spectra at late times, regardless of the fact that they originate from different models with damped fluctuations. On the other hand, the halo mass function is more sensitive to the form of the linear power spectrum. Exploiting this result, we propose a two parameter model of the transfer function in generic damped scenarios, and show that this parametrisation works as well as the standard three parameter models for the scales on which the linear spectrum is relevant.

We propose a new filter, a smooth-$k$ space filter, to use in the Press-Schechter approach to model the dark matter halo mass function which overcomes shortcomings of other filters. We test this against the mass function measured in N-body simulations. We find that the commonly used sharp-$k$ filter fails to reproduce the behaviour of the halo mass function at low masses measured from simulations of models with a sharp truncation in the linear power spectrum. We show that the predictions with our new filter agree with the simulation results over a wider range of halo masses for both damped and undamped power spectra than is the case with the sharp-$k$ and real-space top-hat filters.

We propose a new mechanism to generate a lepton asymmetry based on the vacuum CP-violating phase transition (CPPT). This approach differs from classical thermal leptogenesis as a specific seesaw model, and its UV completion, need not be specified. The lepton asymmetry is generated via the dynamically realised coupling of the Weinberg operator during the phase transition. This mechanism provides a connection with low-energy neutrino observables.

We investigate the impact of thermal velocities in N-body simulations of structure formation in warm dark matter models. Adopting the commonly used approach of adding thermal velocities, randomly selected from a Fermi-Dirac distribution, to the gravitationally-induced velocities of the simulation particles, we compare the matter and velocity power spectra measured from CDM and WDM simulations, in the latter case with and without thermal velocities. This prescription for adding thermal velocities introduces numerical noise into the initial conditions, which influences structure formation. At early times, the noise affects dramatically the power spectra measured from simulations with thermal velocities, with deviations of the order of $\sim \mathcal{O}(10)$ (in the matter power spectra) and of the order of $\sim \mathcal{O}(10^2)$ (in the velocity power spectra) compared to those extracted from simulations without thermal velocities. At late times, these effects are less pronounced with deviations of less than a few percent. Increasing the resolution of the N-body simulation shifts these discrepancies to higher wavenumbers. We also find that spurious haloes start to appear in simulations which include thermal velocities at a mass that is $\sim$3 times larger than in simulations without thermal velocities.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

In the thermal dark matter (DM) paradigm, primordial interactions between DM and Standard Model particles are responsible for the observed DM relic density. In Boehm et al. (2014), we showed that weak-strength interactions between DM and radiation (photons or neutrinos) can erase small-scale density fluctuations, leading to a suppression of the matter power spectrum compared to the collisionless cold DM (CDM) model. This results in fewer DM subhaloes within Milky Way-like DM haloes, implying a reduction in the abundance of satellite galaxies. Here we use very high resolution N-body simulations to measure the dynamics of these subhaloes. We find that when interactions are included, the largest subhaloes are less concentrated than their counterparts in the collisionless CDM model and have rotation curves that match observational data, providing a new solution to the "too big to fail" problem.

We study lensing by voids in Cubic Galileon and Nonlocal gravity cosmologies, which are examples of theories of gravity that modify the lensing potential. We find voids in the dark matter and halo density fields of N-body simulations and compute their lensing signal analytically from the void density profiles, which we show are well fit by a simple analytical formula. In the Cubic Galileon model, the modifications to gravity inside voids are not screened and they approximately double the size of the lensing effects compared to GR. The difference is largely determined by the direct effects of the fifth force on lensing and less so by the modified density profiles. For this model, we also discuss the subtle impact on the force and lensing calculations caused by the screening effects of haloes that exist in and around voids. In the Nonlocal model, the impact of the modified density profiles and the direct modifications to lensing are comparable, but they boost the lensing signal by only $\approx 10\%$, compared with that of GR. Overall, our results suggest that lensing by voids is a promising tool to test models of gravity that modify lensing.

We show that, in the presence of massive neutrinos, the Galileon gravity model provides a very good fit to the current CMB temperature, CMB lensing and BAO data. This model, which we dub ${\nu} \rm{Galileon}$, when assuming its stable attractor background solution, contains the same set of free parameters as $\Lambda\rm{CDM}$, although it leads to different expansion dynamics and nontrivial gravitational interactions. The data provide compelling evidence ($\gtrsim 6\sigma$) for nonzero neutrino masses, with $\Sigma m_\nu \gtrsim 0.4\ {\rm eV}$ at the $2\sigma$ level. Upcoming precision terrestrial measurements of the absolute neutrino mass scale therefore have the potential to test this model. We show that CMB lensing measurements at multipoles $l \lesssim 40$ will be able to discriminate between the ${\nu} \rm{Galileon}$ and $\Lambda\rm{CDM}$ models. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model is consistent with local determinations of the Hubble parameter. The presence of massive neutrinos lowers the value of $\sigma_8$ substantially, despite of the enhanced gravitational strength on large scales. Unlike $\Lambda\rm{CDM}$, the ${\nu} \rm{Galileon}$ model predicts a negative ISW effect, which is difficult to reconcile with current observational limits.

We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the clusters. We discuss a number of important issues that need to be taken into account, associated with the use of nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, the same as in $\Lambda$CDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the cluster radius. However, at distances $\sim \left[2-20\right] {\rm Mpc}/h$ from the cluster center, we find that the surrounding force profiles are enhanced by $\sim20-40\%$ in the Cubic Galileon model. This has an impact on dynamical mass estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can also be applied to the Cubic Galileon model.

Interactions between dark matter (DM) and radiation (photons or neutrinos) in the early Universe suppress density fluctuations on small mass scales. Here we perform a thorough analysis of structure formation in the fully non-linear regime using N-body simulations for models with DM-radiation interactions and compare the results to a traditional calculation in which DM only interacts gravitationally. Significant differences arise due to the presence of interactions, in terms of the number of low-mass DM haloes and their properties, such as their spin and density profile. These differences are clearly seen even for haloes more massive than the scale on which density fluctuations are suppressed. We also show that semi-analytical descriptions of the matter distribution in the non-linear regime fail to reproduce our numerical results, emphasizing the challenge of predicting structure formation in models with physics beyond collisionless DM.

We model the linear and nonlinear growth of large scale structure in the Cubic Galileon gravity model, by running a suite of N-body cosmological simulations using the {\tt ECOSMOG} code. Our simulations include the Vainshtein screening effect, which reconciles the Cubic Galileon model with local tests of gravity. In the linear regime, the amplitude of the matter power spectrum increases by $\sim 25%$ with respect to the standard $\Lambda$CDM model today. The modified expansion rate accounts for $\sim 20%$ of this enhancement, while the fifth force is responsible for only $\sim 5%$. This is because the effective unscreened gravitational strength deviates from standard gravity only at late times, even though it can be twice as large today. In the nonlinear regime ($k \gtrsim 0.1 h\rm{Mpc}^{-1}$), the fifth force leads to only a modest increase ($\lesssim 8%$) in the clustering power on all scales due to the very efficient operation of the Vainshtein mechanism. Such a strong effect is typically not seen in other models with the same screening mechanism. The screening also results in the fifth force increasing the number density of halos by less than 10%, on all mass scales. Our results show that the screening does not ruin the validity of linear theory on large scales which anticipates very strong constraints from galaxy clustering data. We also show that, whilst the model gives an excellent match to CMB data on small angular scales ($l \gtrsim 50$), the predicted integrated Sachs-Wolf effect is in tension with Planck/WMAP results.

The cold dark matter (CDM) model faces persistent challenges on small scales. In particular, taken at face value, the model significantly overestimates the number of satellite galaxies around the Milky Way. Attempts to solve this problem remain open to debate and have even led some to abandon CDM altogether. However, current simulations are limited by the assumption that dark matter feels only gravity. Here, we show that including interactions between CDM and radiation (photons or neutrinos) leads to a dramatic reduction in the number of satellite galaxies, alleviating the Milky Way satellite problem and indicating that physics beyond gravity may be essential to make accurate predictions of structure formation on small scales. The methodology introduced here gives constraints on dark matter interactions that are significantly improved over those from the cosmic microwave background.

We study the nonlinear growth of structure in nonlocal gravity models with the aid of N-body simulation and the spherical collapse and halo models. We focus on a model in which the inverse-squared of the d'Alembertian operator acts on the Ricci scalar in the action. For fixed cosmological parameters, this model differs from $\Lambda{\rm CDM}$ by having a lower late-time expansion rate and an enhanced and time-dependent gravitational strength ($\sim 6\%$ larger today). Compared to $\Lambda{\rm CDM}$ today, in the nonlocal model, massive haloes are slightly more abundant (by $\sim 10\%$ at $M \sim 10^{14} M_{\odot}/h$) and concentrated ($\approx 8\%$ enhancement over a range of mass scales), but their linear bias remains almost unchanged. We find that the Sheth-Tormen formalism describes the mass function and halo bias very well, with little need for recalibration of free parameters. The fitting of the halo concentrations is however essential to ensure the good performance of the halo model on small scales. For $k \gtrsim 1 h/{\rm Mpc}$, the amplitude of the nonlinear matter and velocity divergence power spectra exhibits a modest enhancement of $\sim 12\%$ to $15\%$, compared to $\Lambda{\rm CDM}$ today. This suggests that this model might only be distinguishable from $\Lambda{\rm CDM}$ by future observational missions. We point out that the absence of a screening mechanism may lead to tensions with Solar System tests due to local time variations of the gravitational strength, although this is subject to assumptions about the local time evolution of background averaged quantities.

We use the latest CMB data from Planck, together with BAO measurements, to constrain the full parameter space of Galileon gravity. We constrain separately the three main branches of the theory known as the Cubic, Quartic and Quintic models, and find that all yield a very good fit to these data. Unlike in $\Lambda{\rm CDM}$, the Galileon model constraints are compatible with local determinations of the Hubble parameter and predict nonzero neutrino masses at over $5\sigma$ significance. We also identify that the low-$l$ part of the CMB lensing spectrum may be able to distinguish between $\Lambda{\rm CDM}$ and Galileon models. In the Cubic model, the lensing potential deepens at late times on sub-horizon scales, which is at odds with the current observational suggestion of a positive ISW effect. Compared to $\Lambda$CDM, the Quartic and Quintic models predict less ISW power in the low-$l$ region of the CMB temperature spectrum, and as such are slightly preferred by the Planck data. We illustrate that residual local modifications to gravity in the Quartic and Quintic models may render the Cubic model as the only branch of Galileon gravity that passes Solar System tests.

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

We study spherical collapse in the Quartic and Quintic Covariant Galileon gravity models within the framework of the excursion set formalism. We derive the nonlinear spherically symmetric equations in the quasi-static and weak-field limits, focusing on model parameters that fit current CMB, SNIa and BAO data. We demonstrate that the equations of the Quintic model do not admit physical solutions of the fifth force in high density regions, which prevents the study of structure formation in this model. For the Quartic model, we show that the effective gravitational strength deviates from unity at late times ($z \lesssim 1$), becoming larger if the density is low, but smaller if the density is high. This shows that the Vainshtein mechanism at high densities is not enough to screen all of the modifications of gravity. This makes halos that collapse at $z \lesssim 1$ feel an overall weaker gravity, which suppresses halo formation. However, the matter density in the Quartic model is higher than in standard $\Lambda$CDM, which boosts structure formation and dominates over the effect of the weaker gravity. In the Quartic model there is a significant overabundance of high-mass halos relative to $\Lambda$CDM. Dark matter halos are also less biased than in $\Lambda$CDM, with the difference increasing appreciably with halo mass. However, our results suggest that the bias may not be small enough to fully reconcile the predicted matter power spectrum with LRG clustering data.

We develop a numerical algorithm to solve the high-order nonlinear derivative-coupling equation associated with the quartic Galileon model, and implement it in a modified version of the RAMSES N-body code to study the effect of the Galileon field on the large-scale matter clustering. The algorithm is tested for several matter field configurations with different symmetries, and works very well. This enables us to perform the first simulations for a quartic Galileon model which provides a good fit to the cosmic microwave background (CMB) anisotropy, supernovae and baryonic acoustic oscillations (BAO) data. Our result shows that the Vainshtein mechanism in this model is very efficient in suppressing the spatial variations of the scalar field. However, the time variation of the effective Newtonian constant caused by the curvature coupling of the Galileon field cannot be suppressed by the Vainshtein mechanism. This leads to a significant weakening of the strength of gravity in high-density regions at late times, and therefore a weaker matter clustering on small scales. We also find that without the Vainshtein mechanism the model would have behaved in a completely different way, which shows the crucial role played by nonlinearities in modified gravity theories and the importance of performing self-consistent N-body simulations for these theories.

We present the first constraints on the full parameter space of the Galileon modified gravity model, considering both the cosmological parameters and the coefficients which specify the additional terms in the Lagrangian due to the Galileon field, which we call the Galileon parameters. We use the latest cosmic microwave background measurements, along with distance measurements from supernovae and baryonic acoustic oscillations, performing a Monte Carlo Markov Chain exploration of the 9-dimensional parameter space. The integrated Sachs-Wolfe signal can be very different in Galileon models compared to standard gravity, making it essential to use the full CMB data rather than the CMB distance priors. We demonstrate that meaningful constraints are only possible in the Galileon parameter space after taking advantage of a scaling degeneracy. We find that the Galileon model can fit the WMAP 9-year results better than the standard \Lambda-Cold Dark Matter model, but gives a slightly worse fit overall once lower redshift distance measurements are included. The best-fitting cosmological parameters (e.g. matter density, scalar spectral index, fluctuation amplitude) can differ by more than 2\sigma\ in the Galileon model compared with \Lambda CDM. We highlight other potential constraints of the Galileon model using galaxy clustering and weak lensing measurements.

We study the cosmology of Galileon modified gravity models in the linear perturbation regime. We derive the fully covariant and gauge invariant perturbed field equations using two different methods, which give consistent results, and solve them using a modified version of the {\tt CAMB} code. We find that, in addition to modifying the background expansion history and therefore shifting the positions of the acoustic peaks in the cosmic microwave background (CMB) power spectrum, the Galileon field can cluster strongly from early times, and causes the Weyl gravitational potential to grow, rather than decay, at late times. This leaves clear signatures in the low-$l$ CMB power spectrum through the modified integrated Sachs-Wolfe effect, strongly enhances the linear growth of matter density perturbations and makes distinctive predictions for other cosmological signals such as weak lensing and the power spectrum of density fluctuations. The quasi-static approximation is shown to work quite well from small to the near-horizon scales. We demonstrate that Galileon models display a rich phenomenology due to the large parameter space and the sensitive dependence of the model predictions on the Galileon parameters. Our results show that some Galileon models are already ruled out by present data and that future higher significance galaxy clustering, ISW and lensing measurements will place strong constraints on Galileon gravity.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We propose the liquid-scintillator detector LENA (Low Energy Neutrino Astronomy) as a next-generation neutrino observatory on the scale of 50 kt. The outstanding successes of the Borexino and KamLAND experiments demonstrate the large potential of liquid-scintillator detectors in low-energy neutrino physics. LENA's physics objectives comprise the observation of astrophysical and terrestrial neutrino sources as well as the investigation of neutrino oscillations. In the GeV energy range, the search for proton decay and long-baseline neutrino oscillation experiments complement the low-energy program. Based on the considerable expertise present in European and international research groups, the technical design is sufficiently mature to allow for an early start of detector realization.

The distribution of angles subtended between pairs of galaxies and the line of sight,which is uniform in real space, is distorted by their peculiar motions, and has been proposed as a probe of cosmic expansion. We test this idea using N-body simulations of structure formation in a cold dark matter universe with a cosmological constant and in two variant cosmologies with different dark energy models. We find that the distortion of the distribution of angles is sensitive to the nature of dark energy. However, for the first time, our simulations also reveal dependences of the normalization of the distribution on both redshift and cosmology that have been neglected in previous work. This introduces systematics that severely limit the usefulness of the original method. Guided by our simulations, we devise a new, improved test of the nature of dark energy. We demonstrate that this test does not require prior knowledge of the background cosmology and that it can even distinguish between models that have the same baryonic acoustic oscillations and dark matter halo mass functions. Our technique could be applied to the completed BOSS galaxy redshift survey to constrain the expansion history of the Universe to better than 2%. The method will also produce different signals for dark energy and modified gravity cosmologies even when they have identical expansion histories, through the different peculiar velocities induced in these cases.

Cosmological and astrophysical observations provide increasing evidence of the existence of dark matter in our Universe. Dark matter particles with a mass above a few GeV can be captured by the Sun, accumulate in the core, annihilate, and produce high energy neutrinos either directly or by subsequent decays of Standard Model particles. We investigate the prospects for indirect dark matter detection in the IceCube/DeepCore neutrino telescope and its capabilities to determine the dark matter mass.

Future galaxy surveys hope to distinguish between the dark energy and modified gravity scenarios for the accelerating expansion of the Universe using the distortion of clustering in redshift space. The aim is to model the form and size of the distortion to infer the rate at which large scale structure grows. We test this hypothesis and assess the performance of current theoretical models for the redshift space distortion using large volume N-body simulations of the gravitational instability process. We simulate competing cosmological models which have identical expansion histories - one is a quintessence dark energy model with a scalar field and the other is a modified gravity model with a time varying gravitational constant - and demonstrate that they do indeed produce different redshift space distortions. This is the first time this approach has been verified using a technique that can follow the growth of structure at the required level of accuracy. Our comparisons show that theoretical models for the redshift space distortion based on linear perturbation theory give a surprisingly poor description of the simulation results. Furthermore, the application of such models can give rise to catastrophic systematic errors leading to incorrect interpretation of the observations. We show that an improved model is able to extract the correct growth rate. Further enhancements to theoretical models of redshift space distortions, calibrated against simulations, are needed to fully exploit the forthcoming high precision clustering measurements.

Despite being very successful in explaining the wide range of precision experimental results obtained so far, the Standard Model (SM) of elementary particles fails to address two of the greatest observations of the recent decades: tiny but nonzero neutrino masses and the well-known problem of missing mass in the Universe. Typically the new models beyond the SM explain only one of these observations. Instead, in the present article, we take the view that they both point towards the same new extension of the Standard Model. The new particles introduced are responsible simultaneously for neutrino masses and for the dark matter of the Universe. The stability of dark matter and the smallness of neutrino masses are guaranteed by a U(1) global symmetry, broken to a remnant Z_2. The canonical seesaw mechanism is forbidden and neutrino masses emerge at the loop level being further suppressed by the small explicit breaking of the U(1) symmetry. The new particles and interactions are invoked at the electroweak scale and lead to rich phenomenology in colliders, in lepton flavour violating rare decays and in direct and indirect dark matter searches, making the model testable in the coming future.

The anisotropy of clustering in redshift space provides a direct measure of the growth rate of large scale structure in the Universe. Future galaxy redshift surveys will make high precision measurements of these distortions, and will potentially allow us to distinguish between different scenarios for the accelerating expansion of the Universe. Accurate predictions are needed in order to distinguish between competing cosmological models. We study the distortions in the redshift space power spectrum in $\Lambda$CDM and quintessence dark energy models, using large volume N-body simulations, and test predictions for the form of the redshift space distortions. We find that the linear perturbation theory prediction by Kaiser (1987) is a poor fit to the measured distortions, even on surprisingly large scales $k \ge 0.05 h$Mpc$^{-1}$. An improved model for the redshift space power spectrum, including the non-linear velocity divergence power spectrum, is presented and agrees with the power spectra measured from the simulations up to $k \sim 0.2 h$Mpc$^{-1}$. We have found a density-velocity relation which is cosmology independent and which relates the non-linear velocity divergence spectrum to the non-linear matter power spectrum. We provide a formula which generates the non-linear velocity divergence $P(k)$ at any redshift, using only the non-linear matter power spectrum and the linear growth factor at the desired redshift. This formula is accurate to better than 5% on scales $k<0.2 h $Mpc$^{-1}$ for all the cosmological models discussed in this paper. Our results will extend the statistical power of future galaxy surveys.

We model the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models which feature late $(z<2)$, rapid transitions towards $w=-1$ in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM, despite being very different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$. We find that a second class of models which feature non-negligible amounts of dark energy at early times cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO. These results highlight the need to accurately model quintessence dark energy in N-body simulations when testing cosmological probes of dynamical dark energy.

The feasibility of a next generation neutrino observatory in Europe is being considered within the LAGUNA design study. To accommodate giant neutrino detectors and shield them from cosmic rays, a new very large underground infrastructure is required. Seven potential candidate sites in different parts of Europe and at several distances from CERN are being studied: Boulby (UK), Canfranc (Spain), Fr\'ejus (France/Italy), Pyh\"asalmi (Finland), Polkowice-Sieroszowice (Poland), Slanic (Romania) and Umbria (Italy). The design study aims at the comprehensive and coordinated technical assessment of each site, at a coherent cost estimation, and at a prioritization of the sites within the summer 2010.

We study the nonlinear growth of cosmic structure in different dark energy models, using large volume N-body simulations. We consider a range of quintessence models which feature both rapidly and slowly varying dark energy equations of state, and compare the growth of structure to that in a universe with a cosmological constant. The adoption of a quintessence model changes the expansion history of the universe, the form of the linear theory power spectrum and can alter key observables, such as the horizon scale and the distance to last scattering. We incorporate these effects into our simulations in stages to isolate the impact of each on the growth of structure. The difference in structure formation can be explained to first order by the difference in growth factor at a given epoch; this scaling also accounts for the nonlinear growth at the 15% level. We find that quintessence models that are different from $\Lambda$CDM both today and at high redshifts $(z \sim 1000)$ and which feature late $(z<2)$, rapid transitions in the equation of state, can have identical baryonic acoustic oscillation (BAO) peak positions to those in $\Lambda$CDM. We find that these models have higher abundances of dark matter haloes at $z>0$ compared to $\Lambda$CDM and so measurements of the mass function should allow us to distinguish these quintessence models from a cosmological constant. However, we find that a second class of quintessence models, whose equation of state makes an early $(z>2)$ rapid transition to $w=-1$, cannot be distinguished from $\Lambda$CDM using measurements of the mass function or the BAO, even if these models have non-negligible amounts of dark energy at early times.

Light sterile neutrinos might mix with the active ones and be copiously produced in the early Universe. In the present paper, a detailed multi-flavor analysis of sterile neutrino production is performed. Making some justified approximations allows us to consider not only neutrino interactions with the primeval medium and neutrino coherence breaking effects, but also oscillation effects arising from the presence of three light (mostly-active) neutrino states mixed with two heavier (mostly-sterile) states. First, we emphasize the underlying physics via an analytical description of sterile neutrino abundances that is valid for cases with small mixing between active and sterile neutrinos. Then, we study in detail the phenomenology of (3+2) sterile neutrino models in light of short-baseline oscillation data, including the LSND and MiniBooNE results. Finally, by using the information provided by this analysis, we obtain the expected sterile neutrino cosmological abundances and then contrast them with the most recent available data from Cosmic Microwave Background and Large Scale Structure observations. We conclude that (3+2) models are significantly more disfavored by the internal inconsistencies between sterile neutrino interpretations of appearance and disappearance short-baseline data themselves, rather than by the used cosmological data.

If the dark matter of the Universe is constituted by weakly interacting massive particles (WIMP), they would accumulate in the core of astrophysical objects as the Sun and annihilate into particles of the Standard Model. High-energy neutrinos would be produced in the annihilations, both directly and via the subsequent decay of leptons, quarks and bosons. While Cherenkov neutrino detectors/telescopes can only count the number of neutrinos above some threshold energy, we study how, by exploiting their energy resolution, large magnetized iron calorimeter and, possibly, liquid argon and totally active scintillator detectors, planned for future long baseline neutrino experiments, have the capability of reconstructing the neutrino spectrum and might provide information on the dark matter properties. In particular, for a given value of the WIMP mass, we show that a future iron calorimeter could break the degeneracy between the WIMP-proton cross section and the annihilation branching ratios, present for Cherenkov detectors, and constrain their values with good accuracy.

It is commonly assumed that the cosmological and astrophysical bounds on the mixings of sterile with active neutrinos are much more stringent than those obtained from laboratory measurements. We point out that in scenarios with a very low reheating temperature T_RH << 100 MeV at the end of (the last episode of) inflation or entropy creation, the abundance of sterile neutrinos becomes largely suppressed with respect to that obtained within the standard framework. Thus, in this case cosmological bounds become much less stringent than usually assumed, allowing sterile neutrinos to be ``visible'' in future experiments. Here, we concentrate on massive (mostly sterile) neutrinos heavier than 1 MeV.

MeV particles have been advocated as Dark Matter (DM) candidates in different contexts. This hypothesis can be tested indirectly by searching for the Standard Model (SM) products of DM self-annihilations. As the signal from DM self-annihilations depends on the square of the DM density, we might expect a sizable flux of annihilation products from our galaxy. Neutrinos are the least detectable particles in the SM and a null signal in this channel would allow to set the most conservative bound on the total annihilation cross section. Here, we show that neutrino detectors with good energy resolution and low energy thresholds can not only set bounds on the annihilation cross section but actually test the hypothesis of the possible existence of MeV DM, i.e. test the values of the cross section required to explain the observed DM density. At present, the data in the (positron) energy interval [18-82] MeV of the Super-Kamiokande experiment is already able to put a very stringent bound on the annihilation cross section for masses between ~15-130 MeV. Future large experiments, like megaton water-Cherenkov or large scintillator detectors, will improve the present limits and, if MeV DM exists, would be able to detect it.

We present a scenario in which a remarkably simple relation linking dark matter properties and neutrino masses naturally emerges. This framework points towards a low energy theory where the neutrino mass originates from the existence of a light scalar dark matter particle in the MeV mass range. A very surprising aspect of this scenario is that the required MeV dark matter is one of the favoured candidates to explain the mysterious emission of 511 keV photons in the centre of our galaxy. A possible interpretation of these findings is that dark matter is the stepping stone of a theory beyond the standard model instead of being an embarrassing relic whose energy density must be accounted for in any successful model building.

A recent observation of the 511 keV electron-positron annihilation line from the Galactic bulge has prompted a debate on the origin of the galactic positrons responsible for this emission. Assuming equilibrium between annihilation and injection of positrons in the Galaxy, we investigate the possibility that positrons were produced by past gamma ray bursts (GRBs). We compare the positron annihilation rate inferred from the observed 511 keV line with the predicted injection rate due to electron-positron pairs leaking out of individual GRBs and those pairs produced by GRB photons in the dense molecular clouds at the Galactic center. We conclude that the proposed scenario is marginally consistent with observations, and can reproduce the observed normalization of the annihilation line only for higher-than-average values of the GRB rate in the Galaxy and/or the molecular clouds optical depth for pair production.

We propose an explanation of the LSND evidence for electron antineutrino appearance based on neutrino decay. We introduce a heavy neutrino, which is produced in pion and muon decays because of a small mixing with muon neutrinos, and then decays into a scalar particle and a light neutrino, predominantly of the electron type. We require values of $g m_4\sim$ few eV, $g$ being the neutrino--scalar coupling and $m_4$ the heavy neutrino mass, e.g. $m_4$ in the range from 1 keV to 1 MeV and $g \sim 10^{-6} - 10^{-3}$. Performing a fit to the LSND data as well as all relevant null-result experiments, we show that all data can be explained within this decay scenario. In the minimal version of the decay model, we predict a signal in the upcoming MiniBooNE experiment corresponding to a transition probability of the same order as seen in LSND. In addition, we show that extending our model to two nearly degenerate heavy neutrinos it is possible to introduce CP violation in the decay, which can lead to a suppression of the signal in MiniBooNE running in the neutrino mode. We briefly discuss signals in future neutrino oscillation experiments, we show that our scenario is compatible with bounds from laboratory experiments, and we comment on implications in astrophysics and cosmology.

In general the gravitino mass and/or the soft supersymmetry breaking masses in the observable sector can be much larger than the TeV scale. Depending on the relation between the masses, new important channels for gravitino production in the early Universe can arise. Gravitinos with a mass above 50 TeV decay before big bang nucleosynthesis, which leads to relaxation of the well known bound on the reheating temperature $T_{\rm R} \leq 10^{10}$ GeV. However, if the heavy gravitinos are produced abundantly in the early Universe, their decay can alter the abundance of the lightest supersymmetric particle. Moreover, they may dominate the energy density of the Universe. Their decay will in this case increase entropy and dilute already created baryon asymmetry and dark matter. Such considerations put new constraints on gravitino and sfermion masses, and the reheating temperature. In this paper we examine various cosmological consequences of supermassive gravitinos. We discuss advnatges and disadvantages of a large reheating temperature in connection with thermal leptogenesis, and find that large parts of the parameter space are opened up for the lightest right-handed (s)neutrino mass. We also discuss the viability of Affleck-Dine baryogenesis under the constraints from gravitino decay, and gravitino production from the decay of Q-balls.

We present here a scenario, based on a low reheating temperature T_R << 100 MeV at the end of (the last episode of) inflation, in which the coupling of sterile neutrinos to active neutrinos can be as large as experimental bounds permit (thus making this neutrino ``visible'' in future experiments). In previous models this coupling was forced to be very small to prevent a cosmological overabundance of sterile neutrinos. Here the abundance depends on how low the reheating temperature is. For example, the sterile neutrino required by the LSND result does not have any cosmological problem within our scenario.

We show that a sterile neutrino with mass in the 1-20 keV range and a small mixing with the electron neutrino can simultaneously explain the origin of the pulsar motions and the dark matter in the universe. An asymmetric neutrino emission from a hot nascent neutron star can be the explanation of the observed pulsar velocities. In addition to the pulsar kick mechanism based on resonant neutrino transitions, we point out a new possibility: an asymmetric off-resonant emission of sterile neutrinos. The two cases correspond to different values of the masses and mixing angles. In both cases we identify the ranges of parameters consistent with the pulsar kick, as well as cosmological constraints.

We study the instability of a scalar field at the end of hybrid inflation, using both analytical techniques and numerical simulations. We improve previous studies by taking the inflaton field fully into account, and show that the range of unstable modes depends sensitively on the velocity of the inflaton field, and thereby on the Hubble rate, at the end of inflation. If topological defects are formed, their number density is determined by the shortest unstable wavelength. Finally, we show that the oscillations of the inflaton field amplify the inhomogeneities in the energy density, leading to local symmetry restoration and faster thermalization. We believe this explains why tachyonic preheating is so effective in transferring energy away from the inflaton zero mode.