Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of $H_0$ measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high $\ell$ polarisation data constrain the sterile neutrino mass to be less than approximately 1~eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in a narrow range around 1~eV remains consistent with both cosmology and SBL data.

Numerical simulations of massive neutrino cosmologies consistently find a spoon-like feature in the non-linear matter power spectrum ratios of cosmological models that differ only in the neutrino mass fraction f_N. Typically, the ratio approaches unity at low wave numbers k, decreases by ~ 10 f_N at k ~ 1 h/Mpc, and turns up again at large k. Using the halo model of large-scale structure, we show that this spoon feature originates in the transition from the two-halo power spectrum to the one-halo power spectrum. The former's sensitivity to f_N rises with k, while that of the latter decreases with k. The presence of this spoon feature is robust with respect to different choices of the halo mass function and the halo density profile, and does not require any parameter tuning within the halo model. We demonstrate that a standard halo model calculation is already able to predict the depth, width, and position of this spoon as well as its evolution with redshift z with remarkable accuracy. Predictions at z >= 1 can be further improved using non-linear perturbative inputs.

We investigate how sterile neutrinos with a range of masses influence cosmology in MeV-scale reheating temperature scenarios. By computing the production of sterile neutrinos through the combination of mixing and scattering in the early Universe, we find that light sterile neutrinos, with masses and mixings as inferred from short-baseline neutrino oscillation experiments, are consistent with big-bang nucleosynthesis (BBN) and cosmic microwave background (CMB) radiation for the reheating temperature of ${\cal O}(1)$ MeV if the parent particle responsible for reheating decays into electromagnetic components (radiative decay). In contrast, if the parent particle mainly decays into hadrons (hadronic decay), the bound from BBN becomes more stringent. In this case, the existence of the light sterile neutrinos can be cosmologically excluded, depending on the mass and the hadronic branching ratio of the parent particle.

We study the impact of assumptions made about the neutrino mass ordering on cosmological parameter estimation with the purpose of understanding whether in the future it will be possible to infer the specific neutrino mass distribution from cosmological data. We find that although the commonly used assumption of a degenerate neutrino hierarchy is manifestly wrong and leads to changes in cosmological observables such as the cosmic microwave background and large scale structure compared to the correct (normal or inverted) neutrino hierarchy, the induced changes are so small that even with extremely optimistic assumptions about future data they will remain undetectable. We are thus able to conclude that while cosmology can probe the neutrino contribution to the cosmic energy density extremely precisely (and hence provide a detection of a non-zero total neutrino mass at high significance), it will not be possible to directly measure the individual neutrino masses.

We propose an alternative approach to the construction of fitting functions to the nonlinear matter power spectrum extracted from $N$-body simulations based on the relative matter power spectrum $\delta(k,a)$, defined as the fractional deviation in the absolute matter power spectrum produced by a target cosmology away from a reference $\Lambda$CDM prediction. From the computational perspective, $\delta(k,a)$ is fairly insensitive to the specifics of the simulation settings, and numerical convergence at the 1%-level can be readily achieved without the need for huge computing capacity. Furthermore, $\delta(k,a)$ exhibits several interesting properties that enable a piece-wise construction of the full fitting function, whereby component fitting functions are sought for single-parameter variations and then multiplied together to form the final product. Then, to obtain 1%-accurate absolute power spectrum predictions for any target cosmology only requires that the community as a whole invests in producing one single ultra-precise reference $\Lambda$CDM absolute power spectrum, to be combined with the fitting function to produce the desired result. To illustrate the power of this approach, we have constructed the fitting function RelFit using only five relatively inexpensive $w$CDM simulations (box length $L=256 h^{-1}$Mpc, $N=1024^3$ particles, initialised at $z_i=49$). In a 6-parameter space spanning $\{\omega_m,A_s,n_s,w,\omega_b,h\}$, the output relative power spectra of RelFit are consistent with the predictions of the CosmicEmu emulator to 1% or better for a wide range of cosmologies up to $k\simeq 10$/Mpc. Thus, our approach could provide an inexpensive and democratically accessible route to fulfilling the 1%-level accuracy demands of the upcoming generation of large-scale structure probes, especially in the exploration of "non-standard" or "exotic" cosmologies on nonlinear scales.

From a theoretical point of view, there is a strong motivation to consider an MeV-scale reheating temperature induced by long-lived massive particles with masses around the weak scale, decaying only through gravitational interaction. In this study, we investigate lower limits on the reheating temperature imposed by big-bang nucleosynthesis assuming both radiative and hadronic decays of such massive particles. For the first time, effects of neutrino self-interactions and oscillations are taken into account in the neutrino thermalization calculations. By requiring consistency between theoretical and observational values of light element abundances, we find that the reheating temperature should conservatively be $T_{\rm RH} \gtrsim 1.8$ MeV in the case of the 100% radiative decay, and $T_{\rm RH} \gtrsim$ 4-5 MeV in the case of the 100% hadronic decays for particle masses in the range of 10 GeV to 100 TeV.

Current upper bounds on the sum of 3 active neutrino masses, $\sum m_{\nu}$ from analyses of cosmological data in the backdrop of $\Lambda\textrm{CDM}+\sum m_{\nu}$ model are close to the minimum sum of neutrino masses required by the inverted hierarchy, which is around 0.1 eV. However, these analyses are usually done with the assumption of degenerate masses, which is not a good approximation any more since the bounds are strong enough that the neutrino mass-squared splittings can no longer be considered negligible. In this work we update the bounds on $\sum m_{\nu}$ from latest publicly available cosmological data and likelihoods while explicitly considering particular neutrino mass hierarchies. In the minimal $\Lambda\textrm{CDM}+\sum m_{\nu}$ model with most recent CMB data from Planck 2018 TT,TE,EE, lowE, and lensing; and BAO data from BOSS DR12, MGS, and 6dFGS, we find that at 95% C.L. the bounds are: $\sum m_{\nu}<0.121$ eV (degenerate), $\sum m_{\nu}<0.146$ eV (normal), $\sum m_{\nu}<0.172$ eV (inverted); i.e., the bounds vary significantly across the different mass orderings. Also, we find that the normal hierarchy is very mildly preferred relative to the inverted: $\Delta \chi^2 \equiv \chi^2_{\textrm{NH}}- \chi^2_{\textrm{IH}} = -0.95$ (best-fit). In this paper we also provide bounds on $\sum m_{\nu}$ considering different hierarchies in various extended cosmological models: $\Lambda\textrm{CDM}+\sum m_{\nu}+r$, $w\textrm{CDM}+\sum m_{\nu}$, $w_0 w_a \textrm{CDM}+\sum m_{\nu}$, $w_0 w_a \textrm{CDM}+\sum m_{\nu}$ with $w(z)\geq -1$, $\Lambda \textrm{CDM} + \sum m_{\nu} + \Omega_k$, and $\Lambda \textrm{CDM} + \sum m_{\nu} + A_{\textrm{Lens}}$. We do not find any strong evidence of normal hierarchy over inverted hierarchy from looking at the $\chi^2$ values in the extended models either.(abstract abridged).

We present $N$-body simulations in which either all, or a fraction of, the cold dark matter decays non-relativistically to a relativistic, non-interacting dark radiation component. All effects from radiation and general relativity are self-consistently included at the level of linear perturbation theory, and our simulation results therefore match those from linear Einstein-Boltzmann codes such as CLASS in the appropriate large-scale limit. We also find that standard, Newtonian $N$-body simulations adequately describe the non-linear evolution at smaller scales ($k \gtrsim 0.1 \, h/{\mathrm{Mpc}}$) in this type of model, provided that the mass of the decaying component is modified correctly, and that the background evolution is correctly treated. That is, for studies of small scales, effects from general relativity and radiation can be safely neglected.

We present $N$-body simulations which are fully compatible with general relativity, with dark energy consistently included at both the background and perturbation level. We test our approach for dark energy parameterised as both a fluid, and using the parameterised post-Friedmann (PPF) formalism. In most cases, dark energy is very smooth relative to dark matter so that its leading effect on structure formation is the change to the background expansion rate. This can be easily incorporated into Newtonian $N$-body simulations by changing the Friedmann equation. However, dark energy perturbations and relativistic corrections can lead to differences relative to Newtonian $N$-body simulations at the tens of percent level for scales $k < (10^{-3} \unicode{x2013} 10^{-2})\,\mathrm{Mpc}^{-1}$, and given the accuracy of upcoming large scale structure surveys such effects must be included. In this paper we will study both effects in detail and highlight the conditions under which they are important. We also show that our $N$-body simulations exactly reproduce the results of the Boltzmann solver CLASS for all scales which remain linear.

Bayesian parameter inference depends on a choice of prior probability distribution for the parameters in question. The prior which makes the posterior distribution maximally sensitive to data is called the Jeffreys prior, and it is completely determined by the response of the likelihood to changes in parameters. Under the assumption that the likelihood is a Gaussian distribution, the Jeffreys prior is a constant, i.e. flat. However, if one parameter is constrained by physical considerations, the Gaussian approximation fails and the flat prior is no longer the Jeffreys prior. In this paper we compute the correct Jeffreys prior for a multivariate normal distribution constrained in one dimension, and we apply it to the sum of neutrino masses $\Sigma m_\nu$ and the tensor-to-scalar ratio $r$. We find that one-dimensional marginalised posteriors for these two parameters change considerably and that the 68% and 95% Bayesian upper limits increase by 9% and 4% respectively for $\Sigma m_\nu$ and 22% and 3% for $r$. Adding the prior to an existing chain can be done as a trivial importance sampling in the final step of the analysis proces.

In this paper we investigate a time-varying neutrino mass model, motivated by the mild tension between cosmic microwave background (CMB) measurements of the matter fluctuations and those obtained from low-redshift data. We modify the minimal case of the model proposed by Dvali and Funcke (2016) that predicts late neutrino mass generation in a post-recombination cosmic phase transition, by assuming that neutrino asymmetries allow for the presence of relic neutrinos in the late-time Universe. We show that, if the transition is supercooled, current cosmological data (including CMB temperature, polarization and lensing, baryon acoustic oscillations, and Type Ia supernovae) prefer the scale factor $a_s$ of the phase transition to be very large, peaking at $a_s\sim 1$, and therefore supporting a cosmological scenario in which neutrinos are almost massless until very recent times. We find that in this scenario the cosmological bound on the total sum of the neutrino masses today is significantly weakened compared to the standard case of constant-mass neutrinos, with $\sum m_\nu<4.8$~eV at 95\% confidence, and in agreement with the model predictions. The main reason for this weaker bound is a large correlation arising between the dark energy and neutrino components in the presence of false vacuum energy that converts into the non-zero neutrino masses after the transition. This result provides new targets for the coming KATRIN and PTOLEMY experiments. We also show that the time-varying neutrino mass model considered here does not provide a clear explanation to the existing cosmological $\Omega_m$-$\sigma_8$ discrepancies.

Cosmological $N$-body simulations are typically purely run with particles using Newtonian equations of motion. However, such simulations can be made fully consistent with general relativity using a well-defined prescription. Here, we extend the formalism previously developed for $\Lambda$CDM cosmologies with massless neutrinos to include the effects of massive, but light neutrinos. We have implemented the method in two different $N$-body codes, CONCEPT and PKDGRAV, and demonstrate that they produce consistent results. We furthermore show that we can recover all appropriate limits, including the full GR solution in linear perturbation theory at the per mille level of precision.

Including massive neutrinos in $N$-body simulations is a challenging task due to the large thermal velocities of the neutrinos. In particle based codes this leads to problems of shot-noise due to insufficient sampling of the neutrino momentum distribution function. In this paper we investigate the benefits and drawbacks of a scheme first suggested in a paper by Banerjee et al. in which the initial neutrino distribution is symmetrised in momentum space. We confirm that this method reduce shot-noise significantly, but we also find that it generates some spurious power in the neutrino power spectrum at intermediate and small scales. We speculate that this happens because many neutrinos in the simulation sample the same underlying dark matter structures while moving through the simulation. By carefully tuning the number of directions in momentum space of the initial neutrino distribution we show that some improvements can be made over the case where initial neutrino directions are purely random. At redshifts $z\gtrsim 3$ the method works very well, but at smaller redshifts significant improvements are not possible due to the spurious power generation.

In this paper the non-linear effect of massive neutrinos on cosmological structures is studied in a conceptually new way. We have solved the non-linear continuity and Euler equations for the neutrinos on a grid in real space in $N$-body simulations, and closed the Boltzmann hierarchy at the non-linear Euler equation using the stress and pressure perturbations from linear theory. By comparing with state-of-the art cosmological neutrino simulations, we are able to simulate the non-linear neutrino power spectrum very accurately. This translates into a negligible error in the matter power spectrum, and so our CONCEPT code is ideally suited for extracting the neutrino mass from future high precision non-linear observational probes such as Euclid.

We study models in which neutrino masses are generated dynamically at cosmologically late times. Our study is purely phenomenological and parameterized in terms of three effective parameters characterizing the redshift of mass generation, the width of the transition region, and the present day neutrino mass. We also study the possibility that neutrinos become strongly self-interacting at the time where the mass is generated. We find that in a number of cases, models with large present day neutrino masses are allowed by current CMB, BAO and supernova data. The increase in the allowed mass range makes it possible that a non-zero neutrino mass could be measured in direct detection experiments such as KATRIN. Intriguingly we also find that there are allowed models in which neutrinos become strongly self-interacting around the epoch of recombination.

Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and the impact of general relativity (GR) beyond the background expansion. This approximation can be relaxed and we discuss three different approaches that are accurate to leading order in GR. For simulations that start at redshift less than about 100 we find that the presence of early radiation typically leads to percent-level effects on the numerical power spectra at large scales. Our numerical results agree across the three methods, and we conclude that all of the three methods are suitable for simulations in a standard cosmology. Two of the methods modify the N-body evolution directly, while the third method can be applied as a post-processing prescription.

One of the open questions in modern cosmology is the small scale crisis of the cold dark matter paradigm. Increasing attention has recently been devoted to self-interacting dark matter models as a possible answer. However, solving the so-called "missing satellites" problem requires in addition the presence of an extra relativistic particle (dubbed dark radiation) scattering with dark matter in the early universe. Here we investigate the impact of different theoretical models devising dark matter dark radiation interactions on large scale cosmological observables. We use cosmic microwave background data to put constraints on the dark radiation component and its coupling to dark matter. We find that the values of the coupling allowed by the data imply a cut-off scale of the halo mass function consistent with the one required to match the observations of satellites in the Milky Way.

We have calculated the non-linear effects of generic fermionic and bosonic hot dark matter components in cosmological N-body simulations. For sub-eV masses, the non-linear power spectrum suppression caused by thermal free-streaming resembles the one seen for massive neutrinos, whereas for masses larger than 1eV, the non-linear relative suppression of power is smaller than in linear theory. We furthermore find that in the non-linear regime, one can map fermionic to bosonic models by performing a simple transformation.

We perform a Fourier space decomposition of the dynamics of non-linear cosmological structure formation in LCDM models. From N-body simulations involving only cold dark matter we calculate 3-dimensional non-linear density, velocity divergence and vorticity Fourier realizations, and use these to calculate the fully non-linear mode coupling integrals in the corresponding fluid equations. Our approach allows for a reconstruction of the amount of mode coupling between any two wavenumbers as a function of redshift. With our Fourier decomposition method we identify the transfer of power from larger to smaller scales, the stable clustering regime, the scale where vorticity becomes important, and the suppression of the non-linear divergence power spectrum as compared to linear theory. Our results can be used to improve and calibrate semi-analytical structure formation models.

We study the expected variance of measurements of the Hubble constant, $H_0$, as calculated in either linear perturbation theory or using non-linear velocity power spectra derived from $N$-body simulations. We compare the variance with that obtained by carrying out mock observations in the N-body simulations, and show that the estimator typically used for the local Hubble constant in studies based on perturbation theory is different from the one used in studies based on N-body simulations. The latter gives larger weight to distant sources, which explains why studies based on N-body simulations tend to obtain a smaller variance than that found from studies based on the power spectrum. Although both approaches result in a variance too small to explain the discrepancy between the value of $H_0$ from CMB measurements and the value measured in the local universe, these considerations are important in light of the percent determination of the Hubble constant in the local universe.

In the preprint arxiv:1703.03425 "strong evidence" for the normal neutrino mass ordering is claimed. The authors obtain Bayesian odds of 42:1 in favour of the normal ordering. Their conclusion is based on adopting a flat logarithmic prior for the three neutrino masses. Such an assumption favours a hierarchical spectrum for the masses, which is much easier to accommodate for the normal mass ordering, and hence their prior assumption makes the inverted ordering much less likely a priori. We argue that the claimed "evidence" for normal ordering is almost entirely driven by the adopted prior and not due to the data itself.

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.

With the upcoming sky survey with the Large Synoptic Survey Telescope a great sample of type Ia supernovae will be observed, allowing for a precise mapping of the velocity structure of the universe. Since the source of peculiar velocities is variations in the density field, cosmological parameters related to the matter distribution can subsequently be extracted from the velocity power spectrum. One way to quantify this is through the angular power spectrum of radial peculiar velocities on spheres at different redshifts. We investigate how well this observable can be measured, despite the problems caused by areas with no information. To obtain a realistic distribution of supernovae, we create mock supernova catalogs by using a semi-analytical code for galaxy formation on the merger trees extracted from N-body simulations. We measure the cosmic variance in the velocity power spectrum by repeating the procedure many times for differently located observers, and vary different aspects of the analysis, such as the observer environment, to see how this affects the measurements. Our results confirm the findings from earlier studies regarding the precision with which the angular velocity power spectrum can be determined in the near future. This level of precision has been found to imply, that the angular velocity power spectrum from type Ia supernovae is competitive in its potential to measure parameters such as $\sigma_8$. This makes the peculiar velocity power spectrum from type Ia supernovae a promising new observable, which deserves further attention.

Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian approximation, thus discarding general relativistic effects such as the coupling between matter and radiation ($\equiv$ photons and neutrinos). In this paper we investigate novel hybrid simulations which incorporate interactions between radiation and matter to the leading order in General Relativity, whilst evolving the matter dynamics in full non-linearity according to Newtonian theory. Our hybrid simulations come with a relativistic space-time and make it possible to investigate structure formation in a unified framework. In the present work we focus on simulations initialized at $z=99$, and show that the extracted matter power spectrum receives up to $3\%$ corrections on very large scales through radiation. Our numerical findings compare favourably with linear analytical results from Fidler et al. (2016), from which we deduce that there cannot be any significant non-linear mode-coupling induced through linear radiation corrections.

We propose a simple method to quantify a possible exclusion of the inverted neutrino mass ordering from cosmological bounds on the sum of the neutrino masses. The method is based on Bayesian inference and allows for a calculation of the posterior odds of normal versus inverted ordering. We apply the method for a specific set of current data from Planck CMB data and large-scale structure surveys, providing an upper bound on the sum of neutrino masses of 0.14 eV at 95% CL. With this analysis we obtain posterior odds for normal versus inverted ordering of about 2:1. If cosmological data is combined with data from oscillation experiments the odds reduce to about 3:2. For an exclusion of the inverted ordering from cosmology at more than 95% CL, an accuracy of better than 0.02 eV is needed for the sum. We demonstrate that such a value could be reached with planned observations of large scale structure by analysing artificial mock data for a EUCLID-like survey.

The Short BaseLine (SBL) neutrino oscillation anomalies hint at the presence of a sterile neutrino with a mass of around 1 eV. However, such a neutrino is incompatible with cosmological data, in particular observations of the Cosmic Microwave Background (CMB) anisotropies. However, this conclusion can change by invoking new physics. One possibility is to introduce a secret interaction in the sterile neutrino sector mediated by a light pseudoscalar. In this pseudoscalar model, CMB data prefer a sterile neutrino mass that is fully compatible with the mass ranges suggested by SBL anomalies. In addition, this model predicts a value of the Hubble parameter which is completely consistent with local measurements.

I review the current status of structure formation bounds on neutrino properties such as mass and energy density. I also discuss future cosmological bounds as well as a variety of different scenarios for reconciling cosmology with the presence of light sterile neutrinos.

This paper explores the possibility of determining the spin of the supermassive black hole (SMBH) in Sgr A*, by using secondary images of stars orbiting the SMBH. The photons propagate close to the SMBH and their trajectories probe the space time in a region where the spin of the SMBH is important. We find the appearance of spikes in the secondary image, which depends on the angular momentum and spin axis of the SMBH and study the specific case of the star S2 in detail. The spikes has a magnitude of $\sim 29$ in the K-band and the required angular resolution is of order 15-20 $\mu$as. The combination of these two requirements poses an extreme observational challenge, but might be possible with interferometric observations in the sub-mm regime. The next possible time frame for observing this effect on the star S2 is in the late 2017 and then it repeats with the period of the star.

The effects of local inhomogeneities on low redshift $H_0$ determinations are studied by estimating the redshift-distance relation of mock sources in N-body simulations. The results are compared to those obtained using the standard approach based on Hubble's law. The comparison shows a clear tendency for the standard approach to yield lower values of $H_0$ than the approach based on the scheme using light rays. The difference is, however, small. More precisely, it is found that the overall effect of inhomogeneities on the determination of $H_0$ is a small increase in the local estimates of about $0.3\%$ compared to the results obtained with Hubble's law, when based on a typical distribution of supernovae in the redshift range $0.01 < z < 0.1$. The overall conclusion of the study is a verification of the results that have earlier been obtained by using Hubble's law: The effects of inhomogeneities on local $H_0$ estimates are not significant enough to make it plausible that differences in high- and low-redshift estimates of $H_0$ are due to small inhomogeneities within the setting of standard cosmology.

An expression for the average redshift drift in a statistically homogeneous and isotropic dust universe is given. The expression takes the same form as the expression for the redshift drift in FLRW models. It is used for a proof-of-principle study of the effects of backreaction on redshift drift measurements by combining the expression with two-region models. The study shows that backreaction can lead to positive redshift drift at low redshifts, exemplifying that a positive redshift drift at low redshifts does not require dark energy. Moreover, the study illustrates that models without a dark energy component can have an average redshift drift observationally indistinguishable from that of the standard model according to the currently expected precision of ELT measurements. In an appendix, spherically symmetric solutions to Einstein's equations with inhomogeneous dark energy and matter are used to study deviations from the average redshift drift and effects of local voids.

We study in detail how neutrino perturbations can be followed in linear theory by using only terms up to $l=2$ in the Boltzmann hierarchy. We provide a new approximation to the third moment and demonstrate that the neutrino power spectrum can be calculated to a precision of better than $\sim$ 5% for masses up to $\sim$ 1 eV. The matter and CMB power spectra can be calculated far more precisely and typically at least a factor of a few better than with existing approximations. We then proceed to study how the neutrino power spectrum can be reliably calculated even in the presence of non-linear gravitational clustering by using the full non-linear gravitational potential derived from semi-analytic methods based on $N$-body simulations in the Boltzmann evolution hierarchy. Our results agree extremely well with results derived from $N$-body simulations.

Within the framework of the Standard Model of particle physics and standard cosmology, observations of the Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations (BAO) set stringent bounds on the sum of the masses of neutrinos. If these bounds are satisfied, the upcoming KATRIN experiment which is designed to probe neutrino mass down to $\sim 0.2$ eV will observe only a null signal. We show that the bounds can be relaxed by introducing new interactions for the massive active neutrinos, making neutrino masses in the range observable by KATRIN compatible with cosmological bounds. Within this scenario, neutrinos convert to new stable light particles by resonant production of intermediate states around a temperature of $T\sim$ keV in the early Universe, leading to a much less pronounced suppression of density fluctuations compared to the standard model.

The simplest standard ray tracing scheme employing the Born and Limber approximations and neglecting lens-lens coupling is used for computing the convergence along individual rays in mock N-body data based on Szekeres swiss cheese and onion models. The results are compared with the exact convergence computed using the exact Szekeres metric combined with the Sachs formalism. A comparison is also made with an extension of the simple ray tracing scheme which includes the Doppler convergence. The exact convergence is reproduced very precisely as the sum of the gravitational and Doppler convergences along rays in Lemaitre-Tolman-Bondi models. This is not the case when the models are based on non-symmetric Szekeres models. For such models, there is a significant deviation between the exact and ray traced paths and hence also the corresponding convergences.

It is proposed to use exact, cosmologically relevant solutions to Einstein's equations to accurately quantify the precision of ray tracing techniques through Newtonian N-body simulations. As an initial example of such a study, the recipe in (Green & Wald, 2012) for going between N-body results and a perturbed FLRW metric in the Newtonian gauge is used to study light propagation through quasi-spherical Szekeres models. The study is conducted by deriving a set of ODEs giving an expression for the angular diameter distance in the Newtonian gauge metric. The accuracy of the results obtained from the ODEs is estimated by using the ODEs to determine the distance-redshift relation in mock N-body data based on quasi-spherical Szekeres models. The results are then compared to the exact relations. From this comparison it is seen that the obtained ODEs can accurately reproduce the distance-redshift relation along both radial and non-radial geodesics in spherically symmetric models. The reproduction of geodesics in non-symmetric Szekeres models is slightly less accurate, but still good. These results indicate that the employment of perturbed FLRW metrics for standard ray tracing techniques yields fairly accurate results, at least regarding distance-redshift relations. It is possible though, that this conclusion will be rendered invalid if other typical ray tracing approximations are included and if light is allowed to travel through several structures instead of just one.

Sterile neutrinos in the electronvolt mass range are hinted at by a number of terrestrial neutrino experiments. However, such neutrinos are highly incompatible with data from the Cosmic Microwave Background and large scale structure. This paper discusses how charging sterile neutrinos under a new pseudoscalar interaction can reconcile eV sterile neutrinos with terrestrial neutrino data. We show that this model can reconcile eV sterile neutrinos in cosmology, providing a fit to all available data which is way better than the standard $\Lambda$CDM model with one additional fully thermalized sterile neutrino. In particular it also prefers a value of the Hubble parameter much closer to the locally measured value.

Sterile neutrinos are thermalised in the early Universe via oscillations with the active neutrinos for certain mixing parameters. The most detailed calculation of this thermalisation process involves the solution of the momentum-dependent quantum kinetic equations, which track the evolution of the neutrino phase space distributions. Until now the collision terms in the quantum kinetic equations have always been approximated using equilibrium distributions, but this approximation has never been checked numerically. In this work we revisit the sterile neutrino thermalisation calculation using the full collision term, and compare the results with various existing approximations in the literature. We find a better agreement than would naively be expected, but also identify some issues with these approximations that have not been appreciated previously. These include an unphysical production of neutrinos via scattering and the importance of redistributing momentum through scattering, as well as details of Pauli blocking. Finally, we devise a new approximation scheme, which improves upon some of the shortcomings of previous schemes.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. Such sterile neutrinos are incompatible with cosmology because they suppress structure formation unless they can be prevented from thermalising in the early Universe or removed by subsequent decay or annihilation. Here we present a novel scenario in which both sterile neutrinos and dark matter are coupled to a new, light pseudoscalar. This can prevent thermalisation of sterile neutrinos and make dark matter sufficiently self-interacting to have an impact on galactic dynamics and possibly resolve some of the known problems with the standard cold dark matter scenario. Even more importantly it leads to a strongly self-interacting plasma of sterile neutrinos and pseudoscalars at late times and provides an excellent fit to CMB data. The usual cosmological neutrino mass problem is avoided by sterile neutrino annihilation to pseudoscalars. The preferred value of $H_0$ is substantially higher than in standard $\Lambda$CDM and in much better agreement with local measurements.

We study the potential of a future, large-volume photometric survey to constrain the axion mass $m_a$ in the hot dark matter limit. Future surveys such as Euclid will have significantly more constraining power than current observations for hot dark matter. Nonetheless, the lowest accessible axion masses are limited by the fact that axions lighter than $\sim 0.15$ eV decouple before the QCD epoch, assumed here to occur at a temperature $T_{\rm QCD} \sim 170$ MeV; this leaves an axion population of such low density that its late-time cosmological impact is negligible. For larger axion masses, $m_a \gtrsim 0.15$ eV, where axions remain in equilibrium until after the QCD phase transition, we find that a Euclid-like survey combined with Planck CMB data can detect $m_a$ at very high significance. Our conclusions are robust against assumptions about prior knowledge of the neutrino mass. Given that the proposed IAXO solar axion search is sensitive to $m_a\lesssim 0.2$ eV, the axion mass range probed by cosmology is nicely complementary.

We have carefully studied how local measurements of the Hubble constant, $H_0$, can be influenced by a variety of different parameters related to survey geometry, depth, and size, as well as observer position in space. Our study is based on N-body simulations of structure in the standard $\Lambda$CDM model and our conclusion is that the expected variance in measurements of $H_0$ is far too small to explain the current discrepancy between the low value of $H_0$ inferred from measurements of the cosmic microwave background (CMB) by the Planck collaboration and the value measured directly in the local universe by use of Type Ia supernovae. This conclusion is very robust and does not change with different assumptions about survey geometry and effective sky coverage or observer position in space.

We assess the potential of a future large-volume photometric redshift survey to constrain observational inflationary parameters using three large-scale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Planck-like CMB measurements, we find that the spectral index n_s can be constrained to a 1 sigma precision of up to 0.0025. The sensitivity to the running of the spectral index can potentially improve to 0.0017, roughly a factor of five better than the present 1 sigma~constraint from Planck and auxiliary CMB data, allowing us to test the assumptions of the slow-roll scenario with unprecedented accuracy. Interestingly, neither CMB+shear nor CMB+galaxy nor CMB+clusters alone can achieve this level of sensitivity; it is the combined power of all three probes that conspires to break the different parameter degeneracies inherent in each type of observations. We make our forecast software publicly available via download or upon request from the authors.

We provide updated bounds on non-standard neutrino interactions based on data from the Planck satellite as well as auxiliary cosmological measurements. Two types of models are studied - A Fermi-like 4-point interaction and an interaction mediated by a light pseudoscalar - and we show that these two models are representative of models in which neutrinos either decouple or recouple in the early Universe. Current cosmological data constrain the effective 4-point coupling to be $G_X \leq \left(0.06 \, {\rm GeV}\right)^{-2}$, corresponding to $G_X \leq 2.5 \times 10^7 G_F$. For non-standard pseudoscalar interactions we set a limit on the diagonal elements of the dimensionless coupling matrix, $g_{ij}$, of $g_{ii} \leq 1.2 \times 10^{-7}$. For the off-diagonal elements which induce neutrino decay the bound is significantly stronger, corresponding to $g_{ij} \leq 2.3 \times 10^{-11}(m/0.05 \, {\rm eV})^{-2}$, or a lifetime constraint of $\tau \geq 1.2 \times 10^{9} \, {\rm s} \, (m/0.05 \, {\rm eV})^{3} \,$. This is currently the strongest known bound on this particular type of neutrino decay. We finally note that extremely strong neutrino self-interactions which completely suppress anisotropic stress over all of cosmic history are very highly disfavored by current data ($\Delta \chi^2 \sim 10^4$).

The recent discovery of B-modes in the polarization pattern of the Cosmic Microwave Background by the BICEP2 experiment has important implications for neutrino physics. We revisit cosmological bounds on light sterile neutrinos and show that they are compatible with all current cosmological data provided that the mass is relatively low. Using CMB data, including BICEP-2, we find an upper bound of $m_s < 0.85$ eV ($2\sigma$ Confidence Level). This bound is strengthened to 0.48 eV when HST measurements of $H_0$ are included. However, the inclusion of SZ cluster data from the Planck mission and weak gravitational measurements from the CFHTLenS project favours a non-zero sterile neutrino mass of $0.44^{+0.11}_{-0.16}$ eV. Short baseline neutrino oscillations, on the other hand, indicate a new mass state around 1.2 eV. This mass is highly incompatible with cosmological data if the sterile neutrino is fully thermalised ($\Delta \chi^2>10$). However, if the sterile neutrino only partly thermalises it can be compatible with all current data, both cosmological and terrestrial.

We perform a detailed forecast on how well a {\sc Euclid}-like survey will be able to constrain dark energy and neutrino parameters from a combination of its cosmic shear power spectrum, galaxy power spectrum, and cluster mass function measurements. We find that the combination of these three probes vastly improves the survey's potential to measure the time evolution of dark energy. In terms of a dark energy figure-of-merit defined as $(\sigma(w_{\mathrm p}) \sigma(w_a))^{-1}$, we find a value of 690 for {\sc Euclid}-like data combined with {\sc Planck}-like measurements of the cosmic microwave background (CMB) anisotropies in a 10-dimensional cosmological parameter space, assuming a $\Lambda$CDM fiducial cosmology. For the more commonly used 7-parameter model, we find a figure-of-merit of 1900 for the same data combination. We consider also the survey's potential to measure dark energy perturbations in models wherein the dark energy is parameterised as a fluid with a nonstandard non-adiabatic sound speed, and find that in an \emph{optimistic} scenario in which $w_0$ deviates by as much as is currently observationally allowed from $-1$, models with $\hat{c}_\mathrm{s}^2 = 10^{-6}$ and $\hat{c}_\mathrm{s}^2 = 1$ can be distinguished at more than $2\sigma$ significance. We emphasise that constraints on the dark energy sound speed from cluster measurements are strongly dependent on the modelling of the cluster mass function; significantly weaker sensitivities ensue if we modify our model to include fewer features of nonlinear dark energy clustering. Finally, we find that the sum of neutrino masses can be measured with a $1 \sigma$ precision of 0.015~eV, (abridged)

Neutrino-neutrino refraction can lead to non-periodic flavor oscillations in dense neutrino gases, and it has been hypothesized that some solutions are chaotic in nature. This is of particular interest in the case of neutrino emission from core-collapse supernovae where the measurement of the spectral shape for different flavors can provide crucial information about both neutrino physics and the physical conditions close to the proto-neutron star. Whether a system is chaotic or not can be assessed by the Lyapunov exponents which quantify the rate of divergence of nearby trajectories in the system. We have done a numerical case study for a simple toy model of two neutrino flavors with two momentum states traveling against each other which is known to exhibit flavor transition instabilities. We find the leading Lyapunov exponent to be positive in all cases, confirming the chaoticity of the system for both the normal and the inverted neutrino mass hierarchy. However, more Lyapunov exponents were approximately zero in the inverted hierarchy compared to the normal which has implications for the stability of the system. To investigate this, we have calculated a generalized set of normal modes, the so-called covariant Lyapunov vectors. The covariant Lyapunov vectors associated with vanishing Lyapunov exponents showed the existence of marginally stable directions in phase space for some cases. While our analysis was done for a toy model example, it should work equally well for more realistic cases of neutrinos streaming from a proto-neutron star and provide valuable insight into the nature of the flavor instability. We finally stress that our approach captures many more properties of the physical system than the linear stability analyses which have previously been performed.

In recent years precision cosmology has become an increasingly powerful probe of particle physics. Perhaps the prime example of this is the very stringent cosmological upper bound on the neutrino mass. However, other aspects of neutrino physics, such as their decoupling history and possible non-standard interactions, can also be probed using observations of cosmic structure. Here, I review the current status of cosmological bounds on neutrino properties and discuss the potential of future observations, for example by the recently approved EUCLID mission, to precisely measure neutrino properties.

Short baseline neutrino oscillation experiments have shown hints of the existence of additional sterile neutrinos in the eV mass range. However, such neutrinos seem incompatible with cosmology because they have too large an impact on cosmic structure formation. Here we show that new interactions in the sterile neutrino sector can prevent their production in the early Universe and reconcile short baseline oscillation experiments with cosmology.

We use cosmological observations in the post-Planck era to derive limits on thermally produced cosmological axions. In the early universe such axions contribute to the radiation density and later to the hot dark matter fraction. We find an upper limit m_a < 0.67 eV at 95% C.L. after marginalising over the unknown neutrino masses, using CMB temperature and polarisation data from Planck and WMAP respectively, the halo matter power spectrum extracted from SDSS-DR7, and the local Hubble expansion rate H_0 released by the Carnegie Hubble Program based on a recalibration of the Hubble Space Telescope Key Project sample. Leaving out the local H_0 measurement relaxes the limit somewhat to 0.86 eV, while Planck+WMAP alone constrain the axion mass to 1.01 eV, the first time an upper limit on m_a has been obtained from CMB data alone. Our axion limit is therefore not very sensitive to the tension between the Planck-inferred H_0 and the locally measured value. This is in contrast with the upper limit on the neutrino mass sum, which we find here to range from 0.27 eV at 95% C.L. combining all of the aforementioned observations, to 0.84 eV from CMB data alone.

New measurements of the cosmic microwave background (CMB) by the Planck mission have greatly increased our knowledge about the Universe. Dark radiation, a weakly interacting component of radiation, is one of the important ingredients in our cosmological model which is testable by Planck and other observational probes. At the moment the possible existence of dark radiation is an unsolved question. For instance, the discrepancy between the value of the Hubble constant, H_0, inferred from the Planck data and local measurements of H_0 can to some extent be alleviated by enlarging the minimal LambdaCDM model to include additional relativistic degrees of freedom. From a fundamental physics point of view dark radiation is no less interesting. Indeed, it could well be one of the most accessible windows to physics beyond the standard model. An example of this is that sterile neutrinos, hinted at in terrestrial oscillation experiments, might also be a source of dark radiation, and cosmological observations can therefore be used to test specific particle physics models. Here we review the most recent cosmological results including a complete investigation of the dark radiation sector in order to provide an overview of models that are still compatible with new cosmological observations. Furthermore we update the cosmological constraints on neutrino physics and dark radiation properties focussing on tensions between data sets and degeneracies among parameters that can degrade our information or mimic the existence of extra species.

Cosmology and short baseline neutrino oscillation data both hint at the existence of light sterile neutrinos with masses in the 1 eV range. Here we perform a detailed analysis of the sterile neutrino scenario using both cosmological and SBL data. We have additionally considered the possibility that the extra neutrino degrees of freedom are not fully thermalised in the early universe. Even when analyzing only cosmological data we find a preference for the existence of massive sterile neutrinos in both (3+1) and (3+2) scenarios, and with the inclusion of SBL data the evidence is formally at the 3.3sigma level in the case of a (3+1) model. Interestingly, cosmological and SBL data both point to the same mass scale of approximately 1 eV. In the (3+1) framework WMAP9+SPT provide a value of the sterile mass eigenstate m_4 = (1.72 \pm 0.65) eV: this result is strenghtened by adding the prior from SBL posterior to m_4 = (1.27 \pm 0.12) eV (m_4 = (1.23 \pm 0.13) eV when SDSS is also considered in the cosmological analysis). In the (3+2) scheme, two additional, non--fully thermalized, neutrinos are compatible with the whole set of cosmological and SBL data, leading to mass values of m_4 = (0.95 \pm 0.30) eV and m_5 = (1.59 \pm 0.49) eV. The inclusion of Planck data does not change our considerations about the mass scale; concerning the extra neutrino degrees of freedom, invoking a partial thermalisation the 3+1 model is still consistent with the latest data.

While the cosmic baryon asymmetry has been measured at high accuracy to be 6.1*10^-10, a corresponding lepton asymmetry could be as large as 10^-2 if it hides in the neutrino sector. It has been known for some time that a large asymmetry can be generated from a small initial asymmetry given the existence of a sterile neutrino with a mass less than the mass of the active neutrino. While the magnitude of the final lepton asymmetry is deterministic, its sign has been conjectured to be chaotic in nature. This has been proven in the single momentum approximation, also known as the quantum rate equations, but has up to now not been established using the full momentum dependent quantum kinetic equations. Here we investigate this problem by solving the quantum kinetic equations for a system of 1 active and 1 sterile neutrino on an adaptive grid. We show that by increasing the resolution oscillations in the lepton asymmetry are eliminated and that the sign of the lepton asymmetry is in fact deterministic. This paper also serves as a launch paper for the adaptive solver LASAGNA which is available at http://users-phys.au.dk/steen.

We perform a detailed forecast on how well a Euclid-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of sum m_nu ~ 0.06 eV in the normal hierarchy can be detected at 1.5 sigma to 2.5 sigma significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from Planck. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4 sigma. Interestingly, neither Planck+shear nor Planck+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to sum m_nu, Euclid-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias (< 1 sigma) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Delta chi_eff ^2< 1).

Recent cosmological data favour additional relativistic degrees of freedom beyond the three active neutrinos and photons, often referred to as 'dark' radiation. Light sterile neutrinos is one of the prime candidates for such additional radiation. However, constraints on sterile neutrinos based on the current cosmological data have been derived using simplified assumptions about thermalisation of the sterile neutrino at the Big Bang Nucleosynthesis (BBN) epoch. These assumptions are not necessarily justified and here we solve the full quantum kinetic equations in the (1 active + 1 sterile) scenario and derive the number of thermalised species just before BBN begins (T~1MeV) for null (L=0) and large (L=0.01) initial lepton asymmetry and for a range of possible mass-mixing parameters. We find that the full thermalisation assumption during the BBN epoch is justified for initial small lepton asymmetry only. Partial or null thermalisation occurs when the initial lepton asymmetry is large.

Future cluster surveys will observe galaxy clusters numbering in the hundred thousands. We consider this work how these surveys can be used to constrain dark energy parameters: in particular, the equation of state parameter w and the non-adiabatic sound speed c_s^2. We demonstrate that, in combination with Cosmic Microwave Background (CMB) observations from Planck, cluster surveys such as that in the ESA Euclid project will be able to determine a time-independent w with subpercent precision. Likewise, if the dark energy sound horizon falls within the length scales probed by the cluster survey, then c_s^2 can be pinned down to within an order of magnitude. In the course of this work, we also investigate the process of dark energy virialisation in the presence of an arbitrary sound speed. We find that dark energy clustering and virialisation can lead to dark energy contributing to the total cluster mass at approximately the 0.1% level at maximum.

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

We present a novel method for implementing massive neutrinos in N-body simulations. Instead of sampling the neutrino velocity distribution by individual point particles we take neutrino free-streaming into account by treating it as an effective redshift dependent sound speed in a perfect isothermal fluid, and assume a relation between the sound speed and velocity dispersion of the neutrinos. Although the method fails to accurately model the true neutrino power spectrum, it is able to calculate the total matter power spectrum to the same accuracy as more complex hybrid neutrino methods, except on very small scales. We also present an easy way to update the publicly available Gadget-2 version with this neutrino approximation.

We study cosmological models that contain sterile neutrinos with eV-range masses as suggested by reactor and short-baseline oscillation data. We confront these models with both precision cosmological data (probing the CMB decoupling epoch) and light-element abundances (probing the BBN epoch). In the minimal LambdaCDM model, such sterile neutrinos are strongly disfavoured by current data because they contribute too much hot dark matter. However, if the cosmological framework is extended to include also additional relativistic degrees of freedom -- beyond the three standard neutrinos and the putative sterile neutrinos, then the hot dark matter constraint on the sterile states is considerably relaxed. A further improvement is achieved by allowing a dark energy equation of state parameter w<-1. While BBN strongly disfavours extra radiation beyond the assumed eV-mass sterile neutrino, this constraint can be circumvented by a small nu_e degeneracy. Any model containing eV-mass sterile neutrinos implies also strong modifications of other cosmological parameters. Notably, the inferred cold dark matter density can shift up by 20 to 75% relative to the standard LambdaCDM value.

Recent data from the OPERA experiment seem to point to neutrinos propagating faster than light. One possible physics explanation for such a result is the existence of light sterile neutrinos which can propagate in a higher dimensional bulk and achieve apparent superluminal velocities when measured by an observer confined to the 4D brane of the standard model. Such a model has the advantage of easily being able to explain the non-observation of superluminal neutrinos from SN1987A. Here we discuss the phenomenological implications of such a model and show that it can provide an explanation for the observed faster than light propagation of neutrinos.

We perform a detailed study of the sensitivity to the anisotropies related to Dark Matter (DM) annihilation in the Isotropic Gamma-Ray Background (IGRB) as measured by the Fermi Large Area Telescope (Fermi-LAT). For the first time, we take into account the effects of the Galactic foregrounds and use a realistic representation of the Fermi-LAT. We consider DM anisotropies of extra-galactic origin and of Galactic origin (which can be generated through annihilation in the Milky Way sub-structures) as opposed to a background of anisotropies generated by sources of astrophysical origin, blazars for example. We find that with statistics from 5 years of observation Fermi is sensitive to a DM contribution at the level of the thermal-relic cross section depending on the DM mass and annihilation mode. The anisotropy method for DM searches has a sensitivity comparable to the usual methods based only on the energy spectrum and thus constitutes an independent and complementary piece of information in the DM puzzle. (abridged)

We propose to use the threshold-free process of neutrino capture on beta-decaying nuclei (NCB) using all available candidate nuclei in the Milky Way as target material in order to detect the presence of the Cosmic neutrino background. By integrating over the lifetime of the galaxy one might be able to see the effect of NCB processes as a slightly eschewed abundance ratio of selected beta-decaying nuclei. First, the candidates must be chosen so that both the mother and daughter nuclei have a lifetime comparable to that of the Milky Way or the signal could be easily washed out by additional decays. Secondly, relic neutrinos have so low energy that their de Broglie wavelengths are macroscopic and they may therefore scatter coherently on the electronic cloud of the candidate atoms. One must therefore compare the cross sections for the two processes (induced beta-decay by neutrino capture, and coherent scattering of the neutrinos on atomic nuclei) before drawing any conclusions. Finally, the density of target nuclei in the galaxy must be calculated. We assume supernovae as the only production source and approximate the neutrino density as a homogenous background. Here we perform the full calculation for 187-Re and 138-La and find that one needs abundance measurements with 24 digit precision in order to detect the effect of relic neutrinos. Or alternatively an enhancement of of the relic neutrino density by a factor of $10^{15}$ to produce an effect within the current abundance measurement precision.

The KATRIN (KArlsruhe TRItium Neutrino) experiment will be analyzing the tritium beta-spectrum to determine the mass of the neutrino with a sensitivity of 0.2 eV (90% C.L.). This approach to a measurement of the absolute value of the neutrino mass relies only on the principle of energy conservation and can in some sense be called model-independent as compared to cosmology and neutrino-less double beta decay. However by model independent we only mean in case of the minimal extension of the standard model. One should therefore also analyse the data for non-standard couplings to e.g. righthanded or sterile neutrinos. As an alternative to the frequentist minimization methods used in the analysis of the earlier experiments in Mainz and Troitsk we have been investigating Markov Chain Monte Carlo (MCMC) methods which are very well suited for probing multi-parameter spaces. We found that implementing the KATRIN chi squared function in the COSMOMC package - an MCMC code using Bayesian parameter inference - solved the task at hand very nicely.

Cosmological and astrophysical measurements provide powerful constraints on neutrino masses complementary to those from accelerators and reactors. Here we provide a guide to these different probes, for each explaining its physical basis, underlying assumptions, current and future reach.

A sterile neutrino with mass in the eV range, mixing with the electron antineutrino, is allowed and possibly even preferred by cosmology and oscillation experiments. If such eV-mass neutrinos exist they provide a much better target for direct detection in beta decay experiments than the active neutrinos which are expected to have sub-eV masses. Their relatively high mass would allow for an easy separation from the primary decay signal in experiments such as KATRIN.

Axions with mass greater than 0.7 eV are excluded by cosmological precision data because they provide too much hot dark matter. While for masses above 20 eV the axion lifetime drops below the age of the universe, we show that the cosmological exclusion range can be extended from 0.7 eV till 300 keV, primarily by the cosmic deuterium abundance: axion decays would strongly modify the baryon-to-photon ratio at BBN relative to the one at CMB decoupling. Additional arguments include neutrino dilution relative to photons by axion decays and spectral CMB distortions. Our new cosmological constraints complement stellar-evolution limits and laboratory bounds.

In the last few years there has been some interest in WIMP Dark Matter models featuring a velocity dependent cross section through the Sommerfeld enhancement mechanism, which is a nonrelativistic effect due to massive bosons in the dark sector. In the first part of this article, we find analytic expressions for the boost factor for three different model potentials, the Coulomb potential, the spherical well and the spherical cone well and compare with the numerical solution of the Yukawa potential. We find that the resonance pattern of all the potentials can be cast into the same universal form. In the second part of the article we perform a detailed computation of the Dark Matter relic density for models having Sommerfeld enhancement by solving the Boltzmann equation numerically. We calculate the expected distortions of the CMB blackbody spectrum from WIMP annihilations and compare these to the bounds set by FIRAS. We conclude that only a small part of the parameter space can be ruled out by the FIRAS observations.

Precision cosmology and big-bang nucleosynthesis mildly favor extra radiation in the universe beyond photons and ordinary neutrinos, lending support to the existence of low-mass sterile neutrinos. We use the WMAP 7-year data, small-scale CMB observations from ACBAR, BICEP and QuAD, the SDSS 7th data release, and measurement of the Hubble parameter from HST observations to derive credible regions for the assumed common mass scale m_s and effective number N_s of thermally excited sterile neutrino states. Our results are compatible with the existence of one or perhaps two sterile neutrinos, as suggested by LSND and MiniBooNE, if m_s is in the sub-eV range.

Cosmology provides an excellent laboratory for testing various aspects of neutrino physics. Here, I review the current status of cosmological searches for neutrino mass, as well as other properties of neutrinos. Future cosmological probes of neutrino properties are also discussed in detail.

The matter power spectrum as derived from large scale structure (LSS) surveys contains two important and distinct pieces of information: an overall smooth shape and the imprint of baryon acoustic oscillations (BAO). We investigate the separate impact of these two types of information on cosmological parameter estimation, and show that for the simplest cosmological models, the broad-band shape information currently contained in the SDSS DR7 halo power spectrum (HPS) is by far superseded by geometric information derived from the baryonic features. An immediate corollary is that contrary to popular beliefs, the upper limit on the neutrino mass m_\nu presently derived from LSS combined with cosmic microwave background (CMB) data does not in fact arise from the possible small-scale power suppression due to neutrino free-streaming, if we limit the model framework to minimal LambdaCDM+m_\nu. However, in more complicated models, such as those extended with extra light degrees of freedom and a dark energy equation of state parameter w differing from -1, shape information becomes crucial for the resolution of parameter degeneracies. This conclusion will remain true even when data from the Planck surveyor become available. In the course of our analysis, we introduce a new dewiggling procedure that allows us to extend consistently the use of the SDSS HPS to models with an arbitrary sound horizon at decoupling. All the cases considered here are compatible with the conservative 95%-bounds \sum m_\nu < 1.16 eV, N_eff = 4.8 \pm 2.0.

We update cosmological hot dark matter constraints on neutrinos and hadronic axions. Our most restrictive limits use 7-year data from the Wilkinson Microwave Anisotropy Probe for the cosmic microwave background anisotropies, the halo power spectrum (HPS) from the 7th data release of the Sloan Digital Sky Survey, and the Hubble constant from Hubble Space Telescope observations. We find 95% C.L. upper limits of \sum m_\nu<0.44 eV (no axions), m_a<0.91 eV (assuming \sum m_\nu=0), and \sum m_\nu<0.41 eV and m_a<0.72 eV for two hot dark matter components after marginalising over the respective other mass. CMB data alone yield \sum m_\nu<1.19 eV (no axions), while for axions the HPS is crucial for deriving m_a constraints. This difference can be traced to the fact that for a given hot dark matter fraction axions are much more massive than neutrinos.

The mass varying neutrino scenario is a model that successfully explains the origin of dark energy while at the same time solves the coincidence problem. The model is, however, heavily constrained by its stability towards the formation of neutrino bound states when the neutrinos become nonrelativistic. We discuss these constraints and find that natural, adiabatic, stable models with the right amount of dark energy today do not exist. Secondly, we explain why using the lightest neutrino, which is still relativistic, as an explanation for dark energy does not work because of a feedback mechanism from the heavier neutrinos.

We use N-body simulations to find the effect of neutrino masses on halo properties, and investigate how the density profiles of both the neutrino and the dark matter components change as a function of the neutrino mass. We compare our neutrino density profiles with results from the N-one-body method and find good agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos is not saturated. Finally we study how the halo mass function changes as a function of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic formulae. Our results are important for surveys which aim at probing cosmological parameters using clusters, as well as future experiments aiming at measuring the cosmic neutrino background directly.

Cosmology provides a unique and very powerful laboratory for testing neutrino physics. Here, I review the current status of cosmological neutrino measurements. Future prospects are also discussed, with particular emphasis on the interplay with experimental neutrino physics. Finally I discuss the possibility of a direct detection of the cosmic neutrino background and its associated anisotropy.

We study the primordial non-Gaussinity predicted from simple models of inflation with a linear potential and superimposed oscillations. This generic form of the potential is predicted by the axion monodromy inflation model, that has recently been proposed as a possible realization of chaotic inflation in string theory, where the monodromy from wrapped branes extends the range of the closed string axions to beyond the Planck scale. The superimposed oscillations in the potential can lead to new signatures in the CMB spectrum and bispectrum. In particular the bispectrum will have a new distinct shape. We calculate the power spectrum and bispectrum of curvature perturbations in the model, as well as make analytic estimates in various limiting cases. From the numerical analysis we find that for a wide range of allowed parameters the model produces a feature in the bispectrum with fnl ~ 50 or larger while the power spectrum is almost featureless. This model is therefore an example of a string-inspired inflationary model which is testable mainly through its non-Gaussian features. Finally we provide a simple analytic fitting formula for the bispectrum which is accurate to approximately 5% in all cases, and easily implementable in codes designed to provide non-Gaussian templates for CMB analyses.

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 discuss current cosmological constraints on axions, as well as future sensitivities. Bounds on axion hot dark matter are discussed first, and subsequently we discuss both current and future sensitivity to models in which axions play the role as cold dark matter, but where the Peccei-Quinn symmetry is not restored during reheating.

The cosmic neutrino background anisotropy is calculated for massive neutrino states by solving the full Boltzmann equation. The effect of weak gravitational lensing, including the Limber approximation, is also derived for massive particles, and subsequently applied to the case of massive neutrinos.

We present the first simulation capable of resolving the structure of neutrino clustering on Mpc scales. The method combines grid- and particle-based methods and achieves very good accuracy on both small and large scales, while keeping CPU consumption under control. Such simulations are not only ideal for calculating the non-linear matter power spectrum but also particularly relevant for studies of how neutrinos cluster in galaxy- or cluster-sized halos. We perform the largest neutrino N-body simulation to date, effectively containing 10 different neutrino hot dark matter components with different thermal properties.

We perform an autocorrelation study of the Auger data with the aim to constrain the number density n_s of ultrahigh energy cosmic ray (UHECR) sources, estimating at the same time the effect on $n_s$ of the systematic energy scale uncertainty and of the distribution of UHECR. The use of global analysis has the advantage that no biases are introduced, either in n_s or in the related error bar, by the a priori choice of a single angular scale. The case of continuous, uniformly distributed sources is nominally disfavored at 99% C.L. and the fit improves if the sources follow the large-scale structure of matter in the universe. The best fit values obtained for the number density of proton sources are within a factor ~2 around n_s ~ 10^{-4}/Mpc^3 and depend mainly on the Auger energy calibration scale, with lower densities being preferred if the current scale is correct. The data show no significant small-scale clustering on scales smaller than a few degrees. This might be interpreted as a signature of magnetic smearing of comparable size, comparable with the indication of a ~3 degrees magnetic deflection coming from cross-correlation results. The effects on the above results of some approximations done is also discussed.

We explore the cosmological sensitivity to the amplitude of isocurvature fluctuations that would be caused by axions in the "anthropic window" where the axion decay constant f_a >> 10^12 GeV and the initial misalignment angle Theta_i << 1. In a minimal Lambda-CDM cosmology extended with subdominant scale-invariant isocurvature fluctuations, existing data constrain the isocurvature fraction to alpha < 0.09 at 95% C.L. If no signal shows up, Planck can improve this constraint to 0.042 while an ultimate CMB probe limited only by cosmic variance in both temperature and E-polarisation can reach 0.017, about a factor of five better than the current limit. In the parameter space of f_a and H_I (Hubble parameter during inflation) we identify a small region where axion detection remains within the reach of realistic cosmological probes.

We present a novel, fast and precise method for including the effect of light neutrinos in cosmological N-body simulations. The effect of the neutrino component is included by using the linear theory neutrino perturbations in the calculation of the gravitational potential in the N-body simulation. By comparing this new method with the full non-linear evolution first presented in \cite{Brandbyge1}, where the neutrino component was treated as particles, we find that the new method calculates the matter power spectrum with an accuracy better than 1% for \sum m_\nu \lesssim 0.5 eV at z = 0. This error scales approximately as (\sum m_\nu)^2, making the new linear neutrino method extremely accurate for a total neutrino mass in the range 0.05 - 0.3 eV. At z = 1 the error is below 0.3% for \sum m_\nu \lesssim 0.5 eV and becomes negligible at higher redshifts. This new method is computationally much more efficient than representing the neutrino component by N-body particles.

We revisit the issue of ripples in the primordial power spectra caused by trans-Planckian physics, and the potential for their detection by future cosmological probes. We find that for reasonably large values of the first slow-roll parameter epsilon (> 0.001), a positive detection of trans-Planckian ripples can be made even if the amplitude is as low as 10^-4. Data from the Large Synoptic Survey Telescope (LSST) and the proposed future 21 cm survey with the Fast Fourier Transform Telescope (FFTT) will be particularly useful in this regard. If the scale of inflation is close to its present upper bound, a scale of new physics as high as 0.2 M_Planck could lead to observable signatures.

We explore and compare the performances of two nonlinear correction and scale-dependent biasing models for the extraction of cosmological information from galaxy power spectrum data, especially in the context of beyond-LCDM cosmologies. The first model is the well known Q model, first applied in the analysis of 2dFGRS data. The second, the P model, is inspired by the halo model, in which nonlinear evolution and scale-dependent biasing are encapsulated in a single non-Poisson shot noise term. We find that while both models perform equally well in providing adequate correction for a range of galaxy clustering data in standard LCDM cosmology and in extensions with massive neutrinos, the Q model can give unphysical results in cosmologies containing a subdominant free-streaming dark matter whose temperature depends on the particle mass, e.g., relic thermal axions, unless a suitable prior is imposed on the correction parameter. This last case also exposes the danger of analytic marginalisation, a technique sometimes used in the marginalisation of nuisance parameters. In contrast, the P model suffers no undesirable effects, and is the recommended nonlinear correction model also because of its physical transparency.

We have performed detailed studies of non-linear structure formation in cosmological models with light neutrinos. For the first time the effect of neutrino thermal velocities has been included in a consistent way, and the effect on the matter power spectrum is found to be significant. The effect is large enough to be measured in future, high precision surveys. Additionally, we provide a simple but accurate analytic expression for the suppression of fluctuation power due to massive neutrinos. Finally, we describe a simple and fast method for including the effect of massive neutrinos in large-scale N-body simulations which is accurate at the 1% level for \sum m_nu < 0.15 eV.

The Pierre Auger collaboration has reported a correlation between Ultra-High Energy Cosmic Rays (UHECR) and nearby Active Galactic Nuclei (AGNs) within 75 Mpc. Two of these events fall within 3 degrees from Centaurus A, the nearest AGN, clearly suggesting that this object is a strong UHECR emitter. Here we pursue this hypothesis and forecast the expected rate of ultra-high energy neutrinos in detectors like IceCube. In our baseline model we find a rate of 0.4--0.6 yr^-1 events above a threshold of 100 TeV, the uncertainty of which is mainly related to the poor knowledge of the physical parameters of the source and on the details of the model. This situation will improve with detailed high energy gamma ray measurements of Cen A by the upcoming GLAST satellite. This would make Cen A the first example where the potential of high energy multi-messenger astronomy is finally realized.

The extragalactic cosmic gamma-ray background (CGB) is an interesting channel to look for signatures of dark matter annihilation. In particular, besides the imprint in the energy spectrum, peculiar anisotropy patterns are expected compared to the case of a pure astrophysical origin of the CGB. We take into account the uncertainties in the dark matter clustering properties on sub-galactic scales, deriving two possible anisotropy scenarios. A clear dark matter angular signature is achieved when the annihilation signal receives only a moderate contribution from sub-galactic clumps and/or cuspy haloes. Experimentally, if galactic foregrounds systematics are efficiently kept under control, the angular differences are detectable with the forthcoming GLAST observatory, provided that the annihilation signal contributes to the CGB for a fraction >10-20%. If, instead, sub-galactic structures have a more prominent role, the astrophysical and dark matter anisotropies become degenerate, correspondingly diluting the DM signature. As complementary observables we also introduce the cross-correlation between surveys of galaxies and the CGB and the cross-correlation between different energy bands of the CGB and we find that they provide a further sensitive tool to detect the dark matter angular signatures.

Precision data from cosmology suggest neutrinos stream freely and hence interact very weakly around the epoch of recombination. We study this issue in a simple framework where neutrinos recouple instantaneously and stop streaming freely at a redshift z_i. The latest cosmological data imply z_i < 1500, the exact constraint depending somewhat on the assumed prior on z_i. This bound translates into a limit on the coupling strength between neutrinos and majoron-like particles phi, implying tau > 1 x 10^10 s (m_2/50 meV)^3 for the decay nu_2 -> nu_1+phi.

The arrival directions of ultrahigh energy cosmic rays (UHECRs) may show anisotropies on all scales, from just above the experimental angular resolution up to medium scales and dipole anisotropies. We find that a global comparison of the two-point auto-correlation function of the data with the one of catalogues of potential sources is a powerful diagnostic tool. In particular, this method is far less sensitive to unknown deflections in magnetic fields than cross-correlation studies while keeping a strong discrimination power among source candidates. We illustrate these advantages by considering ordinary galaxies, gamma ray bursts and active galactic nuclei as possible sources. Already the sparse publicly available data suggest that the sources of UHECRs may be a strongly clustered sub-sample of galaxies or of active galactic nuclei. We present forecasts for various cases of source distributions which can be checked soon by the Pierre Auger Observatory.

We update our previous constraints on two-component hot dark matter (axions and neutrinos), including the recent WMAP 5-year data release. Marginalising over sum m_nu provides m_a < 1.02 eV (95% C.L.) for the axion mass. In the absence of axions we find sum m_nu < 0.63 eV (95% C.L.).

A coupling between a light scalar field and neutrinos has been widely discussed as a mechanism for linking (time varying) neutrino masses and the present energy density and equation of state of dark energy. However, it has been pointed out that the viability of this scenario in the non-relativistic neutrino regime is threatened by the strong growth of hydrodynamic perturbations associated with a negative adiabatic sound speed squared. In this paper we revisit the stability issue in the framework of linear perturbation theory in a model independent way. The criterion for the stability of a model is translated into a constraint on the scalar-neutrino coupling, which depends on the ratio of the energy densities in neutrinos and cold dark matter. We illustrate our results by providing meaningful examples both for stable and unstable models.

We present a Markov Chain Monte Carlo global analysis of neutrino parameters using both cosmological and experimental data. Results are presented for the combination of all presently available data from oscillation experiments, cosmology, and neutrinoless double beta decay. In addition we explicitly study the interplay between cosmological, tritium decay and neutrinoless double beta decay data in determining the neutrino mass parameters. We furthermore discuss how the inference of non-neutrino cosmological parameters can benefit from future neutrino mass experiments such as the KATRIN tritium decay experiment or neutrinoless double beta decay experiments.

We consider the inference of the cosmic radiation density, traditionally parameterised as the effective number of neutrino species N_eff, from precision cosmological data. Paying particular attention to systematic effects, notably scale-dependent biasing in the galaxy power spectrum, we find no evidence for a significant deviation of N_eff from the standard value of N_eff^0=3.046 in any combination of cosmological data sets, in contrast to some recent conclusions of other authors. The combination of all available data in the linear regime prefers, in the context of a ``vanilla+N_eff'' cosmological model, 1.1<N_eff<4.8 (95% C.L.) with a best-fit value of 2.6. Adding data at smaller scales, notably the Lyman-alpha forest, we find 2.2<N_eff<5.8 (95% C.L.) with 3.8 as the best fit. Inclusion of the Lyman-alpha data shifts the preferred N_eff upwards because the sigma_8 value derived from the SDSS Lyman-alpha data is inconsistent with that inferred from CMB. In an extended cosmological model that includes a nonzero mass for N_eff neutrino flavours, a running scalar spectral index and a w parameter for the dark energy, we find 0.8<N_eff<6.1 (95% C.L.) with 3.0 as the best fit.

The existence of an unparticle sector, weakly coupled to the standard model, would have a profound impact on supernova (SN) physics. Emission of energy into the unparticle sector from the core of SN1987A would have significantly shortened the observed neutrino burst. The unparticle interaction with nucleons, neutrinos, electrons and muons is constrained to be so weak that it is unlikely to provide any missing-energy signature at colliders. One important exception are models where scale invariance in the hidden sector is broken by the Higgs vacuum expectation value. In this case the SN emission is suppressed by threshold effects.

We calculate the sensitivity of future cosmic microwave background probes and large scale structure measurements from galaxy redshift surveys to the neutrino mass. We find that, for minimal models with few parameters, a measurement of the matter power spectrum over a large range of redshifts has more constraining power than a single measurement at low redshifts. However, this improvement in sensitivity does not extend to larger models. We also quantify how the non-Gaussian nature of the posterior distribution function with respect to the individual cosmological parameter influences such quantities as the sensitivity and the detection threshold. For realistic assumptions about future large scale structure data, the minimum detectable neutrino mass at 95 % C.L. is about 0.05 eV in the context of a minimal 8-parameter cosmological model. In a more general model framework, however, the detection threshold can increase by as much as a factor of three.

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 sum(m_nu). Marginalising over sum(m_nu) provides m_a < 1.2 eV (95% C.L.). In the absence of axions the same data and methods give sum(m_nu) < 0.65 eV (95% C.L.). We also derive limits on m_a for a range of axion-pion couplings up to one order of magnitude larger or smaller than the hadronic value.

Neutrino-neutrino interactions can lead to collective flavour conversion effects in supernovae and in the early universe. We demonstrate that the case of "bipolar" oscillations, where a dense gas of neutrinos and antineutrinos in equal numbers completely converts from one flavour to another even if the mixing angle is small, is equivalent to a pendulum in flavour space. Bipolar flavour conversion corresponds to the swinging of the pendulum, which begins in an unstable upright position (the initial flavour), and passes through momentarily the vertically downward position (the other flavour) in the course of its motion. The time scale to complete one cycle of oscillation depends logarithmically on the vacuum mixing angle. Likewise, the presence of an ordinary medium can be shown analytically to contribute to a logarithmic increase in the bipolar conversion period. We further find that a more complex (and realistic) system of unequal numbers of neutrinos and antineutrinos is analogous to a spinning top subject to a torque. This analogy easily explains how such a system can oscillate in both the bipolar and the synchronised mode, depending on the neutrino density and the size of the neutrino-antineutrino asymmetry. Our simple model applies strictly only to isotropic neutrino gasses. In more general cases, and especially for neutrinos streaming from a supernova core, different modes couple to each other with unequal strength, an effect that can lead to kinematical decoherence in flavour space rather than collective oscillations. The exact circumstances under which collective oscillations occur in non-isotropic media remain to be understood.

Type Ia supernovae are currently the best known standard candles at cosmological distances. In addition to providing a powerful probe of dark energy they are an ideal source of information about the peculiar velocity field of the local universe. Even with the very small number of supernovae presently available it has been possible to measure the dipole and quadrupole of the local velocity field out to z~0.025. With future continuous all-sky surveys like the LSST project the luminosity distances of tens of thousands of nearby supernovae will be measured accurately. This will allow for a determination of the local velocity structure of the universe as a function of redshift with unprecedented accuracy, provided the redshifts of the host galaxies are known. Using catalogues of mock surveys we estimate that future low redshift supernova surveys will be able to probe sigma-8 to a precision of roughly 5% at 95% C.L. This is comparable to the precision in future galaxy and weak lensing surveys and with a relatively modest observational effort it will provide a crucial cross-check on future measurements of the matter power spectrum.

If the diffuse extragalactic gamma ray emission traces the large scale structures of the universe, peculiar anisotropy patterns are expected in the gamma ray sky. In particular, because of the cutoff distance introduced by the absorption of 0.1-10 TeV photons on the infrared/optical background, prominent correlations with the local structures within a range of few hundreds Mpc should be present. We provide detailed predictions of the signal based on the PSCz map of the local universe. We also use mock N-body catalogues complemented with the halo model of structures to study some statistical features of the expected signatures. The results are largely independent from cosmological details, and depend mostly on the index of correlation (or bias) of the sources with respect to the large scale distribution of galaxies. For instance, the predicted signal in the case of a quadratic correlation (as it may happen for a dark matter annihilation contribution to the diffuse gamma flux) differs substantially from a linear correlation case, providing a complementary tool to unveil the nature of the sources of the diffuse gamma ray emission. The chances of the present and future space and ground based observatories to measure these features are discussed.

We propose a novel mechanism for dark energy, based on an extended seesaw for scalar fields, which does not require any new physics at energies below the TeV scale. A very light quintessence mass is usually considered to be technically unnatural, unless it is protected by some symmetry broken at the new very light scale. We propose that one can use an extended seesaw mechanism to construct technically natural models for very light fields, protected by SUSY softly broken above a TeV.

We present a measurement of the velocity flow of the local universe relative to the CMB rest frame, based on the Jha, Riess & Kirshner (2007) sample of 133 low redshift type Ia supernovae. At a depth of 4500 km/s we find a dipole amplitude of 279+-68 km/s in the direction (l,b) = (285+-18,-10+-15), consistent with earlier measurements and with the assumption that the local velocity field is dominated by the Great Attractor region. At a larger depth of 5900 km/s we find a shift in the dipole direction towards the Shapley concentration. We also present the first measurement of the quadrupole term in the local velocity flow at these depths. Finally, we have performed detailed studies based on N-body simulations of the expected precision with which the lowest multipoles in the velocity field can be measured out to redshifts of order 0.1. Our mock catalogues are in good agreement with current observations, and demonstrate that our results are robust with respect to assumptions about the influence of local environment on the type Ia supernova rate.

High-precision data from observation of the cosmic microwave background and the large scale structure of the universe provide very tight constraints on the effective parameters that describe cosmological inflation. Indeed, within a constrained class of LambdaCDM models, the simple lambda phi^4 chaotic inflation model already appears to be ruled out by cosmological data. In this paper, we compute constraints on inflationary parameters within a more general framework that includes other physically motivated parameters such as a nonzero neutrino mass. We find that a strong degeneracy between the tensor-to-scalar ratio r and the neutrino mass prevents lambda phi^4 from being excluded by present data. Reversing the argument, if lambda phi^4 is the correct model of inflation, it predicts a sum of neutrino masses at 0.3-0.5 eV, a range compatible with present experimental limits and within the reach of the next generation of neutrino mass measurements. We also discuss the associated constraints on the dark matter density, the dark energy equation of state, and spatial curvature, and show that the allowed regions are significantly altered. Importantly, we find an allowed range of 0.094 < Omega_c h^2 < 0.136 for the dark matter density, a factor of two larger than that reported in previous studies. This expanded parameter space may have implications for constraints on SUSY dark matter models.

We determine the range of neutrino masses and cosmic radiation content allowed by the most recent CMB and large-scale structure data. In contrast to other recent works, we vary these parameters simultaneously and provide likelihood contours in the two-dimensional parameter space of N_eff}, the usual effective number of neutrino species measuring the radiation density, and \sum m_nu. The allowed range of \sum m_nu and N_eff has shrunk significantly compared to previous studies. The previous degeneracy between these parameters has disappeared, largely thanks to the baryon acoustic oscillation data. The likelihood contours differ significantly if \sum m_nu resides in a single species instead of the standard case of being equally distributed among all flavors. For \sum m_nu=0 we find 2.7 < N_eff < 4.6 at 95% CL while \sum m_nu < 0.62 eV at 95% CL for the standard radiation content.

Surveys of weak gravitational lensing of distant galaxies will be one of the key cosmological probes in the future. We study the ability of such surveys to constrain neutrino masses and the equation of state parameter of the dark energy, focussing on how tomographic information can improve the sensitivity to these parameters. We also provide a detailed discussion of systematic effects pertinent to weak lensing surveys, and the possible degradation of sensitivity to cosmological parameters due to these effects. For future probes such as the Large Synoptic Survey Telescope survey, we find that, when combined with cosmic microwave background data from the Planck satellite, a sensitivity to neutrino masses of sigma(sum m_nu) < 0.05 eV can be reached. This results is robust against variations in the running of the scalar spectral index, the time-dependence of dark energy equation of state, and/or the number of relativistic degrees of freedom.

The Fisher matrix formalism has in recent times become the standard method for predicting the precision with which various cosmological parameters can be extracted from future data. This approach is fast, and generally returns accurate estimates for the parameter errors when the individual parameter likelihoods approximate a Gaussian distribution. However, where Gaussianity is not respected (due, for instance, to strong parameter degeneracies), the Fisher matrix formalism loses its reliability. In this paper, we compare the results of the Fisher matrix approach with those from Monte Carlo simulations. The latter method is based on the publicly available CosmoMC code, but uses synthetic realisations of data sets anticipated for future experiments. We focus on prospective cosmic microwave background (CMB) data from the Planck satellite, with or without CMB lensing information, and its implications for a minimal cosmological scenario with eight parameters and an extended model with eleven parameters. We show that in many cases, the projected sensitivities from the Fisher matrix and the Monte Carlo methods differ significantly, particularly in models with many parameters. Sensitivities to the neutrino mass and the dark matter fraction are especially susceptible to change.

We have studied bounds on the neutrino mass using new data from the WMAP 3 year data, the Sloan Digital Sky Survey measurement of the baryon acoustic peak, the Type Ia supernovae from SNLS, and the Lyman-alpha forest. We find that even in the most general models with a running spectral index where the number of neutrinos and the dark energy equation of state are allowed to vary, the 95% C.L. bound on the sum of neutrino masses is sum m_nu < 0.62 eV (95% C.L.), a bound which we believe to be robust. In the more often used constrained analysis with N_nu =3, w = -1, and alpha_s = 0, we find a bound of 0.48 eV without using the Lyman-alpha data. If Lyman-alpha data is used, the bound shrinks to \sum m_nu < 0.2-0.4 eV (95% C.L.), depending strongly on the Lyman-alpha analysis used.

I review the basics of neutrino cosmology, from the question of neutrino decoupling and the presence of sterile neutrinos to the effects of neutrinos on the cosmic microwave background and large scale structure. Particular emphasis is put on cosmological neutrino mass measurements, both the present bounds and the future prospects.

We have derived new bounds on the relativistic energy density in the Universe from cosmic microwave background (CMB), large scale structure (LSS), and type Ia supernova (SNI-a) observations. In terms of the effective number of neutrino species a bound of N_\nu = 4.2^{+1.2}_{-1.7} is derived at 95% confidence. This bound is significantly stronger than previous determinations, mainly due to inclusion of new CMB and SNI-a observations. The absence of a cosmological neutrino background (N_\nu = 0) is now excluded at 5.4 \sigma. The value of N_\nu is compatible with the value derived from big bang nucleosynthesis considerations, marking one of the most remarkable successes of the standard cosmological model. In terms of the cosmological helium abundance, the CMB, LSS, and SNI-a observations predict a value of 0.240 < Y < 0.281.

Cosmology at present provides the nominally strongest constraint on the masses of standard model neutrinos. However, this constraint extremely dependent on the nature of the dark energy component of the Universe. When the dark energy equation of state parameter is taken as a free (but constant) parameter, the neutrino mass bound is sum m_nu < 1.48 eV (95% C.L.), compared with sum m_nu < 0.65 eV (95% C.L.) in the standard model where the dark energy is in the form of a cosmological constant. This has important consequences for future experiments aimed at the direct measurement of neutrino masses. We also discuss prospects for future cosmological measurements of neutrino masses.

Precision measurements of the acoustic peaks of the cosmic microwave background indicate that neutrinos must be freely streaming at the photon decoupling epoch when T ~ 0.3 eV. This requirement implies restrictive limits on ``secret neutrino interactions,'' notably on neutrino Yukawa couplings with hypothetical low-mass (pseudo)scalars \phi. For diagonal couplings in the neutrino mass basis we find g < 1 x 10^-7, comparable to limits from supernova 1987A. For the off-diagonal couplings and assuming hierarchical neutrino masses we find g < 1 x 10^-11 (0.05 eV/m)^2 where m is the heavier mass of a given neutrino pair connected by g. This stringent limit excludes that the flavor content of high-energy neutrinos from cosmic-ray sources is modified by \nu -> \nu' + \phi decays on their way to Earth.

The present status of our knowledge about the dark matter and dark energy is reviewed. Bounds on the content of cold and hot dark matter from cosmological observations are discussed in some detail. I also review current bounds on the physical properties of dark energy, mainly its equation of state and effective speed of sound.

We study the difference between thermally produced fermionic and bosonic hot dark matter in detail. In the linear regime of structure formation, their distinct free-streaming behaviours can lead to pronounced differences in the matter power spectrum. While not detectable with current cosmological data, such differences will be clearly observable with upcoming large scale weak lensing surveys for particles as light as m_HDM ~ 0.2 eV. In the nonlinear regime, bosonic hot dark matter is not subject to the same phase space constraints that severely limit the amount of fermionic hot dark matter infall into cold dark matter halos. Consequently, the overdensities in fermionic and bosonic hot dark matter of equal particle mass can differ by more than a factor of five in the central part of a halo. However, this unique manifestation of quantum statistics may prove very difficult to detect unless the mass of the hot dark matter particle and its decoupling temperature fall within a very narrow window, 1 < m_HDM/eV < 4 and g* < 30. In this case, hot dark matter infall may have some observable consequences for the nonlinear power spectrum and hence the weak lensing convergence power spectrum at l ~ 10^3 - 10^4 at the percent level.

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.

Observations of the cosmological large-scale structure provide well-established neutrino mass limits. We extend this argument to thermal relic axions. We calculate the axion thermal freeze-out temperature and thus their cosmological abundance on the basis of their interaction with pions. For hadronic axions we find a new mass limit $m_a<1.05$ eV (95% CL), corresponding to a limit on the axion decay constant of $f_a>5.7\times 10^6$ GeV. For other models this constraint is significantly weakened only if the axion-pion coupling is strongly suppressed. For comparison we note that the same approach leads to $\sum m_\nu<0.65$ eV (95% CL) for neutrinos.

We have studied constraints on the equation of state, $w$, and speed of sound, c_s, of the dark energy from a joint analysis of data from the cosmic microwave background, large scale structure and type-Ia supernovae. We find that current observations have no significant sensitivity to c_s. However, there is a slight difference between models in which there are no dark energy perturbations and models in which dark energy behaves as a fluid. Assuming that there are no dark energy perturbations shifts the allowed region for $w$ to slightly higher values. At present models with and without dark energy perturbations provide roughly equally good fits to observations, but the difference is potentially important for future parameter estimations. Finally, we have also performed error forecasts for future measurements of c_s.

We consider the angular power spectrum in a finite universe with different boundary conditions and perform a fit to the CMB, LSS and supernova data. A finite universe could be the consequence of a holographic constraint, giving rise to an effective IR cutoff at the future event horizon. In such a model there is a cosmic duality relating the dark energy equation of state and the power spectrum, which shows a suppression and oscillatory behaviour that is found to describe the low l features extremely well. However, much of the discussion here will also apply if we actually live inside an expanding bubble that describes our universe. The best fit to the CMB and LSS data turns out to be better than in the standard Lambda-CDM model, but when combined with the supernova data, the holographic model becomes disfavored. We speculate on the possible implications.

We investigate a model where neutrinos are strongly coupled to a new, light scalar field. In this model neutrinos annihilate as soon as they become non-relativistic in the early universe, and a non-zero neutrino mass has a marginal effect on the matter power spectrum. However, the angular power spectrum of the cosmic microwave background (CMB) is changed significantly because the strongly interacting fluid of neutrinos and scalars does not experience any anisotropic stress. Such models are strongly disfavoured by current observations. Interestingly, this leads to the conclusion that the relativistic energy density around the epoch of recombination must be in the form of very weakly interacting particles. This conclusion is independent of the specific interaction model.

In this paper we discuss the possibility to measure the Hubble parameter and the slope of galaxy density profiles using future supernova data. With future supernova surveys such as SNAP, large numbers of core collapse supernovae will be discovered, a small fraction of which will be multiply imaged. Measurements of the image separation, flux-ratio, time-delay and lensing foreground galaxy for these systems will provide tight constraints on the slope of galactic halos as well as providing complementary and independent information to other cosmological tests with respect to the Hubble parameter.

We have calculated constraints on the evolution of the equation of state of the dark energy, w(z), from a joint analysis of data from the cosmic microwave background, large scale structure and type-Ia supernovae. In order to probe the time-evolution of w we propose a new, simple parametrization of w, which has the advantage of being transparent and simple to extend to more parameters as better data becomes available. Furthermore it is well behaved in all asymptotic limits. Based on this parametrization we find that w(z=0)=-1.43^{+0.16}_{-0.38} and dw/dz(z=0) = 1.0^{+1.0}_{-0.8}. For a constant w we find that -1.34 < w < -0.79 at 95% C.L. Thus, allowing for a time-varying w shifts the best fit present day value of w down. However, even though models with time variation in w yield a lower chi^2 than pure LambdaCDM models, they do not have a better goodness-of-fit. Rank correlation tests on SNI-a data also do not show any need for a time-varying w.

We propose a new method for identifying new physics imprints on present observational data in cosmology whereby signatures of string theory are clearly distinguished from imprints of possible features on the inflaton potential. Our method relies on the cross-correlations spectra of cosmic shear from large scale structure (LSS) with the CMB temperature anisotropies and E-mode polarization, by using the following properties: inflationary cosmology provides only one source term for all CMB spectra and LSS which highly constrains any deviations from the standard predictions; string theory can add new non-inflationary channels to the source of perturbations as well as modify clustering properties at large scales. Discrepancies in the source terms of correlations and clustering properties provide the evidence for new physics. Models of single-field inflation with a feature are disfavored even with present data. Upcoming WMAP results and future data from weak lensing of LSS will further improve our ability to probe new physics in this manner and could open the first direct window to string theory.

The current status of neutrino cosmology is reviewed, from the question of neutrino decoupling and the presence of sterile neutrinos to the effects of neutrinos on the cosmic microwave background and large scale structure. Particular emphasis is put on cosmological neutrino mass measurements.

We propose an efficient and model independent method for reconstructing the primordial power spectrum from Cosmic Microwave Background (CMB) and large scale structure observations. The algorithm is based on a Monte Carlo principle and therefore very simple to incorporate into existing codes such as Markov Chain Monte Carlo. The algorithm has been used on present cosmological data to test for features in the primordial power spectrum. No significant evidence for features is found, although there is a slight preference for an overall bending of the spectrum, as well as a decrease in power at very large scales. We have also tested the algorithm on mock high precision CMB data, calculated from models with non-scale invariant primordial spectra. The algorithm efficiently extracts the underlying spectrum, as well as the other cosmological parameters in each case. Finally we have used the algorithm on a model where an artificial glitch in the CMB spectrum has been imposed, like the ones seen in the WMAP data. In this case it is found that, although the underlying cosmological parameters can be extracted, the recovered power spectrum can show significant spurious features, such as bending, even if the true spectrum is scale invariant.

We study models in which the universe exits reheating at temperatures in the MeV regime. By combining light element abundance measurements with cosmic microwave background and large scale structure data we find a fairly robust lower limit on the reheating temperature of T_RH > 4 MeV at 95% C.L. However, if the heavy particle whose decay reheats the universe has a direct decay mode to neutrinos, there are some small islands left in parameter space where a reheating temperature as low as 1 MeV is allowed. The derived lower bound on the reheating temperature also leads to very stringent bounds on models with $n$ large extra dimensions. For n=2 the bound on the compactification scale is M > 2000 TeV, and for n=3 it is 100 TeV. These are currently the strongest available bounds on such models.

In theories with large extra dimensions, supernova (SN) cores are powerful sources of Kaluza-Klein (KK) gravitons. A large fraction of these massive particles are gravitationally retained by the newly born neutron star (NS). The subsequent slow KK decays produce potentially observable gamma rays and heat the NS. We here show that the back-absorption of the gravitationally trapped KK gravitons does not significantly change our previous limits. We calculate the graviton emission rate in a nuclear medium by combining the low-energy classical bremsstrahlung rate with detailed-balancing arguments. This approach reproduces the previous thermal emission rate, but it is much simpler and allows for a calculation of the absorption rate by a trivial phase-space transformation. We derive systematically the dependence of the SN and NS limits on the number of extra dimensions.

Resonant particle production, along with many other physical processes which change the effective equation of state (EOS) during inflation, introduces a step-like feature in the primordial power spectrum. We calculate observational constraints on resonant particle production, parameterised in form of an effective step height, N_eff and location in k-space, k_break. Combining data from the cosmic microwave background and the 2dF Galaxy Redshift Survey yields strong constraints in some regions of parameter space, although the range in k-space which can be probed is restricted to k ~ 0.001 - 0.1 h Mpc^-1. We also discuss the implications of our findings for general models which change the effective EOS during inflation.

Cosmological observations have in the past few years become an increasingly powerful method for determining parameters within the neutrino sector, such as the presence of sterile states and the mass of neutrinos. I review the current status of the field in light of recent measurements of the cosmic microwave background by the WMAP collaboration, as well as current large scale galaxy surveys.

Entropy production from scalar decay in the era of low temperatures after neutrino decoupling will change the ratio of the relic neutrino temperature to the CMB temperature, and, hence, the value of N_eff, the effective number of neutrino species. Such scalar decay is relevant to reheating after thermal inflation, proposed to dilute massive particles, like the moduli and the gravitino, featuring in supersymmetric and string theories. The effect of such entropy production on the relic neutrino temperature ratio is calculated in a semi-analytic manner, and a recent lower bound on this ratio, obtained from the WMAP satellite and 2dF galaxy data, is used to set a lower bound of ~ 1.5 x 10^-23 Gev on the scalar decay constant, corresponding to a reheating temperature of about 3.3 Mev.

We investigate the potential of observations of anisotropies in the cosmic microwave background (CMB) and large-scale structure in the Universe to detect possible modifications of standard inflationary models by physics beyond the Planck scale. A generic model of primordial density fluctuations is investigated, and we derive constraints on its parameters from current data. We conclude that the currently available data do not put very stringent constraints on this model. Furthermore, we use simulated power spectra from the Sloan Digital Sky Survey (SDSS) and the Planck satellite to show that it is unlikely that a trans-Planckian signature of this type can be detected in CMB and large-scale structure data.

We measure the average gravitational shear profile of 6 massive clusters (M_vir ~ 10^15 M_sun) at z=0.3 out to a radius ~2h^-1 Mpc. The measurements are fitted to a generalized NFW-like halo model \rho(r) with an arbitrary r -> 0 slope \alpha. The data are well fitted by such a model with a central cusp with \alpha ~ 0.9 - 1.6 (68% confidence interval). For the standard-NFW case \alpha = 1.0, we find a concentration parameter c_vir that is consistent with recent predictions from high-resolution CDM N-body simulations. Our data are also well fitted by an isothermal sphere model with a softened core. For this model, our 1\sigma upper limit for the core radius corresponds to a limit \sigma_star \leq 0.1 cm^2 g^-1 on the elastic collision cross-section in a self-interacting dark matter model.

The paper has been withdrawn

We have performed a thorough analysis of the constraints which can be put on neutrino parameters from cosmological observations, most notably those from the WMAP satellite and the 2dF galaxy survey. For this data we find an upper limit on the sum of active neutrino mass eigenstates of \sum m_nu < 1.0 eV (95% conf.), but this limit is dependent on priors. We find that the WMAP and 2dF data alone cannot rule out the evidence from neutrinoless double beta decay reported by the Heidelberg-Moscow experiment. In terms of the relativistic energy density in neutrinos or other weakly interacting species we find, in units of the equivalent number of neutrino species, N_nu, that N_nu = 4.0+3.0-2.1 (95% conf.). When BBN constraints are added, the bound on N_\nu is 2.6+0.4-0.3 (95% conf.), suggesting that N_nu could possibly be lower than the standard model value of 3. This can for instance be the case in models with very low reheating temperature and incomplete neutrino thermalization. Conversely, if N_nu is fixed to 3 then the data from WMAP and 2dFGRS predicts that 0.2458 < Y_P < 0.2471, which is significantly higher than the observationally measured value. The limit on relativistic energy density changes when a small $\nu_e$ chemical potential is present during BBN. In this case the upper bound on N_nu from WMAP, 2dFGRS and BBN is N_nu < 6.5. Finally, we find that a non-zero \sum m_nu can be compensated by an increase in N_nu. One result of this is that the LSND result is not yet ruled out by cosmological observations.

We have carefully analysed the potential of future Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) measurements to probe neutrino masses. We perform a Fisher matrix analysis on a 9-dimensional cosmological parameter space and find that data from the Planck CMB experiment combined with the Sloan Digital Sky Survey (SDSS) can measure a neutrino mass of 0.12 eV at 95% conf. This is almost at the level of the 0.06 eV mass suggested by current neutrino oscillation data. A future galaxy survey with an order of magnitude larger survey volume than the SDSS would allow for a neutrino mass determination of 0.03-0.05 eV (95% conf.).

We have reanalysed constraints on the equation of state parameter, w_Q = P/rho, of the dark energy, using several cosmological data sets and relaxing the usual constraint w_Q > -1. We find that combining Cosmic Microwave Background, large scale structure and Type Ia supernova data yields a non-trivial lower bound on w_Q. At 95.4% confidence we find, assuming a flat geometry of the universe, a bound of -2.68 < w_Q < -0.78 if w_Q is taken to be a completely free parameter. Reassuringly we also find that the constraint w_Q > -1 does not significantly bias the overall allowed region for w_Q. When this constraint is imposed the 95.4% confidence bound is -1 < w_Q < -0.71. Also, a pure cosmological constant (w = -1) is an excellent fit to all available data. Based on simulations of future data from the Planck CMB experiment and the SNAP and SNfactory supernova experiments we estimate that it will be possible to constrain w_Q at the 5% level in the future.

We review the current status of neutrino cosmology, focusing mainly on the question of the absolute values of neutrino masses and the possibility of a cosmological neutrino lepton asymmetry.

In models with large extra dimensions all gauge singlet fields can in principle propagate in the extra dimensional space. We have investigated possible constraints on majoron models of neutrino masses in which the majorons propagate in extra dimensions. It is found that astrophysical constraints from supernovae are many orders of magnitude stronger than previous accelerator bounds. Our findings suggest that unnatural types of the "see-saw" mechanism for neutrino masses are unlikely to occur in nature, even in the presence of extra dimensions.

We have performed a careful analysis of constraints on the neutrino mass from current cosmological data. Combining data from the cosmic microwave background and the 2dF galaxy survey yields an upper limit on the sum of the three neutrino mass eigenstates of \sum m_nu < 3 eV (95% conf.), without including additional priors. Including data from SNIa observations, Big Bang nucleosynthesis, and HST Hubble key project data on H_0 tightens the limit to \sum m_nu < 2.5 eV (95% conf.). We also perform a Fisher matrix analysis which illustrates the cosmological parameter degeneracies affecting the determination of \sum m_nu.

We study the effect on the cosmic microwave background (CMB) anisotropy and large scale structure (LSS) power spectrum of a scattering interaction between cold dark matter and baryons. This scattering alters the CMB anisotropy and LSS spectrum through momentum transfer between the cold dark matter particles and the baryons. We find that current CMB observations can put an upper limit on the scattering cross section which is comparable with or slightly stronger than previous disk heating constraints at masses greater than 1 GeV, and much stronger at smaller masses. When large-scale structure constraints are added to the CMB limits, our constraint is more stringent than this previous limit at all masses. In particular, a dark matter-baryon scattering cross section comparable to the ``Spergel-Steinhardt'' cross section is ruled out for dark matter mass greater than 1 GeV.

Supernovae (SNe) are copious sources for Kaluza-Klein gravitons which are generic for theories with large extra dimensions. These massive particles are produced with average velocities ~0.5 c so that many of them are gravitationally retained by the SN core. Every neutron star thus has a halo of KK gravitons which decay into nu bar-nu, e^+e^- and gamma gamma on time scales \~10^9 years. The EGRET gamma-flux limits (E_gamma ~ 100 MeV) for nearby neutron stars constrain the fundamental scale for n=2 extra dimensions to M >500 TeV, and M>30 TeV for n=3. The upcoming GLAST satellite is a factor ~30 more sensitive and thus may detect KK decays, for example at the nearby neutron star RX J185635--3754. The requirement that neutron stars are not excessively heated by KK decays implies M>1700 TeV for n=2, and M>60 TeV for n=3.

In the early universe, neutrinos decouple from equilibrium with the electromagnetic plasma at a temperature which is only slightly higher than the temperature where electrons and positrons annihilate. Therefore neutrinos to some extent share in the entropy transfer from e^+e^- to other species, and their final temperature is slightly higher than the canonical value T_nu = (4/11)^{1/3} T_gamma. We study neutrino decoupling in the early universe with effects of neutrino oscillations included, and find that the change in neutrino energy density from e^+ e^- annihilations can be about 2-3% higher if oscillation are included. The primordial helium abundance can be changed by as much as 1.5 x 10^-4 by neutrino oscillations.

We constrain the spectrum of primordial curvature perturbations by using recent Cosmic Microwave Background (CMB) and Large Scale Structure (LSS) data. Specifically, we consider CMB data from the COBE, Boomerang and Maxima experiments, the real space galaxy power spectrum from the IRAS PSCz survey, and the linear matter power spectrum inferred from Ly-alpha forest spectra. We study the case of single field slow roll inflationary models, and we extract bounds on the scalar spectral index, n, the tensor to scalar ratio, r, and the running of the scalar spectral index, dlnk, for various combinations of the observational data. We find that CMB data, when combined with data from Lyman-alpha forest, place strong constraints on the inflationary parameters. Specifically, we obtain n \approx 0.9, r < 0.3 and dlnk \approx 0, indicating that single field hybrid models are ruled out.

If large extra dimensions exist in nature, supernova (SN) cores will emit large fluxes of Kaluza-Klein gravitons, producing a cosmic background of these particles with energies and masses up to about 100 MeV. Radiative decays then give rise to a diffuse cosmic gamma-ray background with E_gamma < 100 MeV which is well in excess of the observations if more than 0.5-1% of the SN energy is emitted into the new channel. This argument complements and tightens the well-known cooling limit from the observed duration of the SN1987A neutrino burst. For two extra dimensions we derive a conservative bound on their radius of R < 0.9 x 10^-4 mm, for three extra dimensions it is R < 1.9 x 10^-7 mm.

New precision Cosmic Microwave Background Radiation (CMBR) anisotropy data are beginning to constrain physics beyond the standard model, for example in the form of additional light particle species. These constraints are complementary to what can be obtained from big bang nucleosynthesis (BBN) considerations because they apply to much later times. We derive a constraint on the equivalent number of neutrino species, N_\nu, from the presently available data. Specifically we analyse two different CMBR data sets to test the robustness of our results. Analyzing only CMBR data yields an upper bound of N_\nu < 17 (95% confidence). Adding large scale structure (LSS) data from the PSC-z survey tightens the upper bound slightly. However, the addition of LSS data gives a non-trivial {\it lower} bound of N_\nu > 1.5/2.5 (95% confidence) for the two data sets. This is the first independent indication of the presence of the cosmological neutrino background which is predicted by the standard model, and seen in big bang nucleosynthesis. The value $N_\nu = 0$ is disfavoured at 3\sigma and 4\sigma for the two data sets respectively.

We have studied cosmological constraints on the number and radii of possible large extra dimensions. If such dimensions exist, Kaluza-Klein (KK) modes are copiously produced at high temperatures in the early universe, and can potentially lead to unacceptable cosmological effects. We show that during reheating, large numbers of KK modes are produced. These modes are not diluted completely by the entropy production during reheating because they are produced non-relativistically. This means that the modes produced during reheating can easily be the dominant component. For instance, for two extra dimensions the bound on their radii from considering only the thermally produced KK modes is R < 1.1 x 10^-4 mm. If the modes produced during reheating are also accounted for, the bound is strengthened to R < 2.2 x 10^-5 mm. This bound is stronger than all other known astrophysical or laboratory limits.

Interest in non-standard recombination scenarios has been spurred by recent cosmic microwave background (CMB) results from BOOMERANG and MAXIMA, which show an unexpectedly low second acoustic peak, resulting in a best-fit baryon density that is 50% larger than the prediction of big-bang nucleosynthesis (BBN). This apparent discrepancy can be avoided if the universe has a non-standard ionization history in which the recombination of hydrogen is significantly delayed relative to the standard model. While future CMB observations may eliminate this discrepancy, it is useful to develop a general framework for analyzing non-standard ionization histories. We develop such a framework, examining non-standard models in which the hydrogen binding energy E_b and the overall expression for the time rate of change of the ionized fraction of electrons are multiplied by arbitrary factors. This set of models includes a number of previously-proposed models as special cases. We find a wide range of models with delayed recombination that are able to fit the CMB data with a baryon density in accordance with BBN, but there are even allowed models with earlier recombination than in the standard model. A generic prediction of these models is that the third acoustic CMB peak should be very low relative to what is found in the standard model. This is the case even for the models with earlier recombination than in the standard model, because here the third peak is lowered by an increased diffusion damping at recombination relative to the standard model. Interestingly, the specific height of the third peak depends sensitively on the model parameters, so that future CMB measurements will be able to distinguish between different non-standard recombination scenarios.

Measurements of fluctuations in the Cosmic Microwave Background Radiation (CMBR) is one of the most promising methods for measuring the fundamental cosmological parameters. However, in order to infer parameters from precision measurements it is necessary to calculate the theoretical fluctuation spectrum to at least the measurement accuracy. Standard treatments assume that electrons, ions and neutral hydrogen are very tightly coupled during the entire recombination history, and that the baryon-photon plasma can be treated as a two-fluid system consisting of baryons and photons interacting via Thomson scattering. We investigate the validity of this approximation by explicitly writing down and solving the full set of Boltzmann equations for electrons, ions, neutral hydrogen and photons. The main correction to the standard treatment is from including Rayleigh scattering between photons and neutral hydrogen, a change of less than 0.1% in the CMBR power spectrum. Our conclusion is thus that the standard treatment of the baryon-photon system is a very good approximation, better than any possible measurement accuracy.

We constrain the primordial fluctuation spectrum P(k) by using the new data on the Cosmic Microwave Background (CMB) from the Boomerang and MAXIMA experiments. Our study is based on slow-roll inflationary models, and we consider the possibility of a running spectral index. Specifically, we expand the power spectrum P(k) to second order in ln(k), thus allowing the power spectrum to ``bend'' in k-space. We show that allowing the power spectrum to bend erases the ability of the present data to measure the tensor to scalar perturbation ratio. Moreover, if the primordial baryon density \Omega_b h^2 is as low as found from Big Bang nucleosynthesis (BBN), the data favor a negative bending of the power spectrum, corresponding to a bump-like feature in the power spectrum around a scale of k=0.004 Mpc^-1.

We have investigated the evolution of the neutrino asymmetry in active-sterile neutrino oscillations in the early universe. We find that there are large regions of parameter space where the asymmetry is extremely sensitive to variations in the initial asymmetry as well as the external parameters (the mass difference and the mixing angle). In these regions the system undergoes chaotic transitions; however, the system is never truly chaotic in the sense that all information about initial conditions is lost. In some cases though, enough information is lost that the final sign of the neutrino asymmetry is stochastic. We discuss the implications of our findings for Big Bang nucleosynthesis (BBN) and the cosmic microwave background (CMB).

The Cosmic Microwave Background Radiation (CMBR) holds information about almost all the fundamental cosmological parameters, and by performing a likelihood analysis of high precision CMBR fluctuation data, these parameters can be inferred. However, this analysis relies on assumptions about the initial power spectrum, which is usually taken to be a featureless power-law, P(k) ~ k^{n_s-1}. Many inflationary models predict power spectra with non-power law features. We discuss the possibility for detecting such features by describing the power spectrum as bins in k-space. This method for power spectrum reconstruction is demonstrated in practise by performing likelihood optimization on synthetic spectra, and the difficulties arising from reconstructing smooth features using discontinuous bins are discussed in detail.

We have performed a likelihood analysis of the recent data on the Cosmic Microwave Background Radiation (CMBR) anisotropy taken by the Boomerang experiment. We find that this data places a strong upper bound on the radiation density present at recombination. Expressed in terms of the equivalent number of neutrino species the $2\sigma$ bound is N_nu < 13, and the standard model prediction, N_nu = 3.04, is completely consistent the the data. This bound is complementary to the one found from Big Bang nucleosynthesis considerations in that it applies to any type of radiation, i.e. it is not flavour sensitive. It also applies to the universe at a much later epoch, and as such places severe limits on scenarios with decaying neutrinos. The bound also yields a firm upper limit on the lepton asymmetry in the universe.

Many independent high resolution simulations have indicated that the standard collisionless cold dark matter model does not reproduce the structure of observed present day galaxies well. Several possible solutions in the form of modifications to the physics of the dark matter particles have been proposed. One of the most promising is warm dark matter (WDM), particles with significant thermal motion in the early universe. It is usually assumed that such particles are relativistically decoupled particles with a mass of approximately 1 keV. However, here we have investigated the possibility that much more massive particles with highly non-thermal spectra could make up warm dark matter. Several possible production mechanisms are reviewed and the only one found to be viable is that the WDM is produced by the non-relativistic decay of some massive species in the early universe. Such very massive warm dark matter could possibly be detected in direct detection experiments, as opposed to standard thermal warm dark matter.

Many independent high-resolution simulations of structure formation in cold dark matter models show that galactic halos should have singular core profiles. This is in stark contrast with observations of both low- and high-surface brightness galaxies, which indicate that the dark matter has almost constant density in the central parts of halos. Basically there are three possible avenues to a solution to the problem, which we discuss in turn. Observations of halo profiles could be more uncertain than previously thought, and higher resolution observations could reveal that spirals do have a singular core feature. The highest resolution simulations do not include a baryonic component, and it is conceivable that violent star formation processes and similar phenomena can destroy the singular dark matter core and lead to an almost constant density core profile. Finally, we discuss in more detail the intriguing possibility that the discrepancy hints at some new exotic physics of the dark matter. Warm dark matter and self-interacting dark matter are two of the most promising candidates.

We study if the neutrino mixing parameters suggested by the atmospheric neutrino anomaly imply chemical equilibrium between mu- and tau-flavored leptons in a supernova (SN) core. The initial flavor-conversion rate would indeed be fast if the nu_mu-nu_tau-mixing angle were not suppressed by second-order refractive effects. The neutrino diffusion coefficients are different for nu_mu, anti-nu_mu, nu_tau and anti-nu_tau so that neutrino transport will create a net mu and tau lepton number density. This will typically lead to a situation where the usual first-order refractive effects dominate, further suppressing the rate of flavor conversion. Altogether, neutrino refraction has the nontrivial consequence of guaranteeing the separate conservation of e, mu, and tau lepton number in a SN core on the infall and cooling time scales, even when neutrino mixing angles are large.

It has been shown by many independent studies that the cold dark matter scenario produces singular galactic dark halos, in strong contrast with observations. Possible remedies are that either the dark matter is warm so that it has significant thermal motion or that the dark matter has strong self interactions. We have combined these ideas to calculate the linear mass power spectrum and the spectrum of cosmic microwave background (CMB) fluctuations for self-interacting warm dark matter. Our results indicate that such models have more power on small scales than is the case for the standard warm dark matter model, with a CMB fluctuation spectrum which is nearly indistinguishable from standard cold dark matter. This enhanced small-scale power may provide better agreement with the observations than does standard warm dark matter.

Recent, very accurate simulations of galaxy formation have revealed that the standard cold dark matter model has great difficulty in explaining the detailed structure of galaxies. One of the major problems is that galactic halos are too centrally concentrated. Dark matter self-interactions have been proposed as a possible means of resolving this inconsistency. Here, we investigate quantitatively the effect of dark matter self interactions on formation of galactic halos. Our numerical framework is extremely simple, while still keeping the essential physics. We confirm that strongly self-interacting dark matter leads to less centrally concentrated structures. Interestingly, we find that for a range of different interaction strengths, the dark matter halos are unstable to particle ejection on a timescale comparable to the Hubble time.

The reconstruction of the CMBR power spectrum from a map represents a major computational challenge to which much effort has been applied. However, once the power spectrum has been recovered there still remains the problem of extracting cosmological parameters from it. Doing this involves optimizing a complicated function in a many dimensional parameter space. Therefore efficient algorithms are necessary in order to make this feasible. We have tested several different types of algorithms and found that the technique known as simulated annealing is very effective for this purpose. It is shown that simulated annealing is able to extract the correct cosmological parameters from a set of simulated power spectra, but even with such fast optimization algorithms, a substantial computational effort is needed.

Recent measurements of type Ia supernovae indicate that distant supernovae are substantially fainter than expected from the standard flat cold dark matter model. One possible explanation is that the energy density in our universe is in fact dominated by a cosmological constant. Another possible solution is that there are large amounts of grey dust in the intergalactic medium. Dust grains can be grey either because they are non-spherical or very large. We have numerically investigated whether grey dust can be emitted from high redshift galaxies without also emitting standard, reddening dust, which would have been visible in the spectra of high redshift objects. Our finding is that grain velocities are almost independent of ellipticity so that if greyness are due to the grains being elongated, it will not be possible to separate grey dust from ordinary dust. We also find that velocities are fairly independent of grain size, but we cannot rule out possible sputtering of small grains, so that large, grey dust grains could be preferentially emitted. Therefore, our conclusion is that grey dust is an unlikely explanation of the data, but we cannot rule it out if the grey dust consists of large, spherical grains.

Measurements of the Cosmic Microwave Background Radiation (CMBR) provide a powerful tool for measuring the primary cosmological parameters. However, there is a large degree of parameter degeneracy in simultaneous measurements of the matter density, Omega_m, and the Hubble parameter, H_0. In the present paper we use the presently available CMBR data together with measurements of the cosmological baryon-to-photon ratio, eta, from Big Bang nucleosynthesis, and the relative mass fraction of baryons in clusters to break the parameter degeneracy in measuring Omega_m and H_0. We find that present data is inconsistent with the standard Omega=1, matter dominated model. Our analysis favours a medium density universe with a rather low Hubble parameter. This is compatible with new measurements of type Ia supernovae, and the joint estimate of the two parameters is Omega_m = 0.45^{+0.07}_{-0.07} and H_0 = 39^{+14}_{-13} km s^-1 Mpc^-1. We stress that the upper bound on the Hubble parameter is likely to be much more uncertain than indicated here, because of the limited number of free parameters in our analysis.

The formation of the cosmic microwave background radiation (CMBR) provides a very powerful probe of the early universe at the epoch of recombination. Specifically, it is possible to constrain the variation of fundamental physical constants in the early universe. We have calculated the effect of a varying electromagnetic coupling constant (\alpha) on the CMBR and find that new satellite experiments should provide a tight constraint on the value of \alpha at recombination which is complementary to existing constraints. An estimate of the obtainable precision is |\dot{\alpha}/\alpha| \leq 7 x 10^{-13} y^{-1} in a realistic experiment.

We investigate in detail the possibility of constraining neutrino decays with data from the cosmic microwave background radiation (CMBR). Two generic decays are considered \nu_H -> \nu_L \phi and \nu_H -> \nu_L \nu_L_bar \nu_L. We have solved the momentum dependent Boltzmann equation in order to account for possible relativistic decays. Doing this we estimate that any neutrino with mass m > 1 eV decaying before the present should be detectable with future CMBR data. Combining this result with other results on stable neutrinos, any neutrino mass of the order 1 eV should be detectable.

We investigate the possible effects that deviations from kinetic equilibrium can have on massive particles as they decouple from chemical equilibrium. Different methods of solving the Boltzmann equation yield significantly different relic number densities of such particles. General considerations concerning the Dirac or Majorana structure of the particles are discussed. It is shown that non-equilibrium effects are small for particles decoupling while strongly non-relativistic, as will be the case for most cold dark matter candidates.

We investigate the possible effect of the residual annihilations of heavy particles after freeze-out from equilibrium in the early universe. An error in previous work on this subject is pointed out and the correct method of solving the Boltzmann equation for this case is developed. For Majorana particles there are significant differences relative to previous work.

The existence of low-mass sterile neutrinos is suggested by the current status of solar and atmospheric neutrinos together with the LSND experiment. In typical four-flavor scenarios, neutrinos would contribute to a cosmic hot dark matter component and to an increased radiation content at the epoch of matter-radiation equality. These effects leave their imprint in sky maps of the cosmic microwave background radiation (CMBR) and may thus be detectable with the precision measurements of the upcoming MAP and PLANCK missions.

Elastic scattering on nucleons, \nu N -> N \nu, is the dominant supernova (SN) opacity source for \mu and \tau neutrinos. The dominant energy- and number-changing processes were thought to be \nu e^- -> e^- \nu and \nu\bar \nu <-> e^+ e^- until Suzuki (1993) showed that the bremsstrahlung process \nu\bar \nu NN <-> NN was actually more important. We find that for energy exchange, the related ``inelastic scattering process'' \nu NN <-> NN \nu is even more effective by about a factor of 10. A simple estimate implies that the \nu_\mu and \nu_\tau spectra emitted during the Kelvin-Helmholtz cooling phase are much closer to that of \nu\bar_e than had been thought previously. To facilitate a numerical study of the spectra formation we derive a scattering kernel which governs both bremsstrahlung and inelastic scattering and give an analytic approximation formula. We consider only neutron-neutron interactions, we use a one-pion exchange potential in Born approximation, nonrelativistic neutrons, and the long-wavelength limit, simplifications which appear justified for the surface layers of a SN core. We include the pion mass in the potential and we allow for an arbitrary degree of neutron degeneracy. Our treatment does not include the neutron-proton process and does not include nucleon-nucleon correlations. Our perturbative approach applies only to the SN surface layers, i.e. to densities below about 10^{14} g cm^{-3}.

The decay of massive neutrinos to final states containing only invisible particles is poorly constrained experimentally. In this letter we describe the constraints that can be put on neutrino mass and lifetime using CMBR measurements. We find that very tight lifetime limits on neutrinos in the mass range 10 eV - 100 keV can be derived using CMBR data from upcoming satellite measurements.

We investigate different types of neutrino hot dark matter with respect to structure formation and anisotropies in the cosmic microwave background radiation (CMBR). The possibility of neutrino hot dark matter produced through the decay of a heavier neutrino by the process nu_H -> nu_L + phi, where phi is a scalar particle, is discussed in detail. This type of dark matter can possibly be distinguished observationally from the standard neutrino dark matter by using new CMBR data from the upcoming satellite missions MAP and PLANCK.

We investigate the non-radiative decay during nucleosynthesis of a massive tau neutrino with mass 0.1 - 1 MeV into an electron neutrino and a scalar or pseudoscalar particle, $\phi$. The full Boltzmann equation is used and shown to give markedly different results than the usual non-relativistic formalism for relativistic or semi-relativistic neutrino decays. Indeed, the region we investigate is where the formalism that has previously been applied to solving this problem is expected to break down. We also compare the nucleosynthesis predictions from this scenario with results from the standard model and with some of the available observational determinations of the primordial abundances. It is found that for relativistic or semi-relativistic decays the helium abundance can be significantly lowered without changing other light element abundances. Since a problem with the standard model of Big Bang nucleosynthesis is that helium appears to be overproduced, a decay of the type we discuss can be a possible solution.

We present a calculation of a neutrino decay scenario in the early Universe. The specific decay is \nu_{2} \to \nu_{1} + \phi, where \phi is a boson. If there is a neutrino mass hierarchy, m_{\nu_{e}} < m_{\nu_{\mu}} < m_{\nu_{\tau}}, we show that it is possible to generate stimulated decay and effects similar to atomic lasing without invoking new neutrinos, even starting from identical neutrino distributions. Under the right circumstances the decay can be to very low momentum boson states thereby producing something similar to a Bose condensate, with possible consequences for structure formation. Finally, we argue that this type of decay may also be important other places in early Universe physics.

Our results for the primordial abundance of a massive Majorana tau neutrino (hep-ph/9603201) published in PRL76, 2848 (1996) were incorrect due to a trivial (but annoying) error in the numerical code. In this ERRATUM we present the corrected results.

The primordial abundance of long-lived heavy Majorana neutrinos is calculated from the full Boltzmann equation. Inclusion of scattering reactions drastically change the predicted abundance of a heavy neutrino species. This loosens the well known mass constraint on MeV neutrinos from Big Bang nucleosynthesis, and allows for the existence of a Majorana tau neutrino with mass m > 11 MeV. Further experimental efforts are therefore needed to investigate the range 11 MeV < m < 24 MeV. Some interesting cosmological consequences of an MeV nu_tau are also pointed out.

A calculation of neutrino decoupling in the early Universe, including full Fermi-Dirac statistics and electron mass dependence in the weak reaction rates, is presented. We find that after decoupling, the electron neutrinos contribute 0.83\% more to the relativistic energy density than in the standard scenario, where neutrinos are assumed not to share the heating from e$^\pm$ annihilation. The corresponding number for muon and tau neutrinos is 0.41\% . This has the consequence of modifying the primordial $^4$He abundance by $\Delta Y=+1.0\times 10^{-4}$, and the cosmological mass limit on light neutrinos by 0.2--0.5 eV.