We investigate generalized interacting dark matter-dark energy scenarios with a time-dependent coupling parameter, allowing also for freedom in the neutrino sector. The models are tested in the phantom and quintessence regimes, characterized by an equation of state $w_x<-1$ and $w_x>-1$, respectively. Our analyses show that for some of the scenarios the existing tensions on the Hubble constant $H_0$ and on the clustering parameter $S_8$ can be significantly alleviated. The relief is either due to \textit{(a)} a dark energy component which lies within the phantom region; or \textit{(b)} the presence of a dynamical coupling in quintessence scenarios. The inclusion of massive neutrinos into the interaction schemes does not affect neither the constraints on the cosmological parameters nor the bounds on the total number or relativistic degrees of freedom $N_{\rm eff}$, which are found to be extremely robust and, in general, strongly consistent with the canonical prediction $N_{\rm eff}=3.045$. The most stringent bound on the total neutrino mass $M_{\nu}$ is $M_{\nu}<0.116$ eV and it is obtained within a quintessence scenario in which the matter mass-energy density is only mildly affected by the presence of a dynamical dark sector coupling.

The Phenomenologically Emergent Dark Energy model, a dark energy model with the same number of free parameters as the flat $\Lambda$CDM, has been proposed as a working example of a minimal model which can avoid the current cosmological tensions. A straightforward question is whether or not the inclusion of massive neutrinos and extra relativistic species may spoil such an appealing phenomenological alternative. We present the bounds on $M_{\nu}$ and $N_{\rm eff}$ and comment on the long standing $H_0$ and $\sigma_8$ tensions within this cosmological framework with a wealth of cosmological observations. Interestingly, we find, at $95\%$ confidence level, and with the most complete set of cosmological observations, $M_{\nu}\sim 0.21^{+0.15}_{-0.14}$ eV and $N_{\rm eff}= 3.03\pm 0.32$ i.e. an indication for a non-zero neutrino mass with a significance above $2\sigma$. The well known Hubble constant tension is considerably easened, with a significance always below the $2\sigma$ level.

We examine interactions between dark matter and dark energy in light of the latest cosmological observations, focusing on a specific model with coupling proportional to the dark energy density. Our data includes Cosmic Microwave Background (CMB) measurements from the Planck 2018 legacy data release, late-time measurements of the expansion history from Baryon Acoustic Oscillations (BAO) and Supernovae Type Ia (SNeIa), galaxy clustering and cosmic shear measurements from the Dark Energy Survey Year 1 results, and the 2019 local distance ladder measurement of the Hubble constant $H_0$ from the Hubble Space Telescope. Considering Planck data both in combination with BAO or SNeIa data reduces the $H_0$ tension to a level which could possibly be compatible with a statistical fluctuation. The very same model also significantly reduces the $\Omega_{\rm m}-\sigma_8$ tension between CMB and cosmic shear measurements. Interactions between the dark sectors of our Universe remain therefore a promising joint solution to these persisting cosmological tensions.

Type Ia Supernovae (SNeIa) used as standardizable candles have been instrumental in the discovery of cosmic acceleration, usually attributed to some form of dark energy (DE). Recent studies have raised the issue of whether intrinsic SNeIa luminosities might evolve with redshift. While the evidence for cosmic acceleration is robust to this possible systematic, the question remains of how much the latter can affect the inferred properties of the DE component responsible for cosmic acceleration. This is the question we address in this work. We use SNeIa distance moduli measurements from the Pantheon and JLA samples. We consider models where the DE equation of state is a free parameter, either constant or time-varying, as well as models where DE and dark matter interact, and finally a model-agnostic parametrization of effects due to modified gravity (MG). When SNeIa data are combined with Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements, we find strong degeneracies between parameters governing the SNeIa systematics, the DE parameters, and the Hubble constant $H_0$. These degeneracies significantly broaden the DE parameter uncertainties, in some cases leading to ${\cal O}(\sigma)$ shifts in the central values. However, including low-redshift Baryon Acoustic Oscillation and Cosmic Chronometer measurements, as well as CMB lensing measurements, considerably improves the previous constraints, and the only remaining effect of the examined systematic is a $\lesssim 40\%$ broadening of the uncertainties on the DE parameters. The constraints we derive on the MG parameters are instead basically unaffected by the systematic in question. We therefore confirm the overall soundness of dark energy properties.

We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Namely, we have included all T2K measurements as of December 2019, the most recent NO$\nu$A antineutrino statistics, and data collected by the Daya Bay and RENO reactor experiments. Concerning the atmospheric and solar sectors, we have also updated our analyses of DeepCore and SNO data, respectively. All in all, these new analyses result in more accurate measurements of $\theta_{13}$, $\theta_{12}$, $\Delta m_{21}^2$ and $|\Delta m_{31}^2|$. The best fit value for the atmospheric angle $\theta_{23}$ lies in the second octant, but first octant solutions remain allowed at $\sim2\sigma$. Regarding CP violation measurements, the preferred value of $\delta$ we obtain is 1.20$\pi$ (1.54$\pi$) for normal (inverted) neutrino mass ordering. These new results should be regarded as extremely robust due to the excellent agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a preference for the normal neutrino mass ordering at the $2.7\sigma$ level. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the significance to $3.1\sigma$ within a conservative approach. A more aggressive data set combination of cosmological observations leads to a stronger preference for normal with respect to inverted mass ordering, at the $3.3\sigma$ level. This cosmological data set provides $2\sigma$ upper limits on the total neutrino mass corresponding to $\sum m_\nu<0.13$ ($0.15$)~eV in the normal (inverted) neutrino mass ordering scenario. These bounds are among the most complete ones in the literature, as they include all currently available neutrino physics inputs.

The observation of space-time variations in fundamental constants would provide strong evidence for the existence of new light degrees of freedom in the theory of Nature. Robustly constraining such scenarios requires exploiting observations that span different scales and probe the state of the Universe at different epochs. In the context of cosmology, both the cosmic microwave background and the Lyman-$\alpha$ forest have proven to be powerful tools capable of constraining variations in electromagnetism, however at the moment there do not exist cosmological probes capable of bridging the gap between recombination and reionization. In the near future, radio telescopes will attempt to measure the 21cm transition of neutral hydrogen during the epochs of reionization and the cosmic dawn (and potentially the tail end of the dark ages); being inherently sensitive to electromagnetic phenomena, these experiments will offer a unique perspective on space-time variations of the fine-structure constant and the electron mass. We show here that large variations in these fundamental constants would produce features on the 21cm power spectrum that may be distinguishable from astrophysical uncertainties. Furthermore, we forecast the sensitivity for the Square Kilometer Array, and show that the 21cm power spectrum may be able to constrain variations at the level of ${\cal O}(10^{-3})$.

In this paper we explore possible extensions of Interacting Dark Energy cosmologies, where Dark Energy and Dark Matter interact non-gravitationally with one another. In particular, we focus on the neutrino sector, analyzing the effect of both neutrino masses and the effective number of neutrino species. We consider the Planck 2018 legacy release data combined with several other cosmological probes, finding no evidence for new physics in the dark radiation sector. The current neutrino constraints from cosmology should be therefore regarded as robust, as they are not strongly dependent on the dark sector physics, once all the available observations are combined. Namely, we find a total neutrino mass $M_\nu<0.15$ eV and a number of effective relativistic degrees of freedom of $N_{\rm eff}=3.03^{+0.33}_{-0.33}$, both at 95\% CL, which are close to those obtained within the $\Lambda$CDM cosmology, $M_\nu<0.12$ eV and $N_{\rm eff}=3.00^{+0.36}_{-0.35}$ for the same data combination.

We explore whether non-standard dark sector physics might be required to solve the existing cosmological tensions. The properties we consider in combination are an interaction between the dark matter and dark energy components, and a dark energy equation of state $w$ different from that of the canonical cosmological constant $w=-1$. In principle, these two parameters are independent. In practice, to avoid early-time, superhorizon instabilities, their allowed parameter spaces are correlated. We analyze three classes of extended interacting dark energy models in light of the 2019 Planck CMB results and Cepheid-calibrated local distance ladder $H_0$ measurements of Riess et al. (R19), as well as recent BAO and SNeIa distance data. We find that in quintessence coupled dark energy models, where $w > -1$, the evidence for a non-zero coupling between the two dark sectors can surpass the $5\sigma$ significance. On the other hand, in phantom coupled dark energy models, there is no such preference for a non-zero dark sector coupling. All the models we consider significantly raise the value of the Hubble constant easing the $H_0$ tension. The addition of low-redshift BAO and SNeIa measurements leaves some residual tension with R19 but at a level that could be justified by a statistical fluctuation. Bayesian evidence considerations mildly disfavour both the coupled quintessence and phantom models, while mildly favouring a coupled vacuum scenario, even when late-time datasets are considered. We conclude that non-minimal dark energy cosmologies, such as coupled quintessence, phantom, or vacuum models, are still an interesting route towards softening existing cosmological tensions, even when low-redshift datasets and Bayesian evidence considerations are taken into account. (abstract severely abridged)

We consider the possibility that dark energy and baryons might scatter off each other. The type of interaction we consider leads to a pure momentum exchange, and does not affect the background evolution of the expansion history. We parametrize this interaction in an effective way at the level of Boltzmann equations. We compute the effect of dark energy-baryon scattering on cosmological observables, focusing on the Cosmic Microwave Background (CMB) temperature anisotropy power spectrum and the matter power spectrum. Surprisingly, we find that even huge dark energy-baryon cross-sections $\sigma_{xb} \sim {\cal O}({\rm b})$, which are generically excluded by non-cosmological probes such as collider searches or precision gravity tests, only leave an insignificant imprint on the observables considered. In the case of the CMB temperature power spectrum, the only imprint consists in a sub-percent enhancement or depletion of power (depending whether or not the dark energy equation of state lies above or below $-1$) at very low multipoles, which is thus swamped by cosmic variance. These effects are explained in terms of differences in how gravitational potentials decay in the presence of a dark energy-baryon scattering, which ultimately lead to an increase or decrease in the late-time integrated Sachs-Wolfe power. Even smaller related effects are imprinted on the matter power spectrum. The imprints on the CMB are not expected to be degenerate with the effects due to altering the dark energy sound speed. We conclude that, while strongly appealing, the prospects for a direct detection of dark energy through cosmology do not seem feasible when considering realistic dark energy-baryon cross-sections. As a caveat, our results hold to linear order in perturbation theory.

In this paper, we explore a scenario where the dark matter is a mixture of interacting and non interacting species. Assuming dark matter-photon interactions for the interacting species, we find that the suppression of the matter power spectrum in this scenario can mimic that expected in the case of massive neutrinos. Our numerical studies include present limits from Planck Cosmic Microwave Background data, which render the strength of the dark matter photon interaction unconstrained when the fraction of interacting dark matter is small. Despite the large entangling between mixed dark matter and neutrino masses, we show that future measurements from the Dark Energy Instrument (DESI) could help in establishing the dark matter and the neutrino properties simultaneously, provided that the interaction rate is very close to its current limits and the fraction of interacting dark matter is at least of $\mathcal{O}\left(10\%\right)$. However, for that region of parameter space where a small fraction of interacting DM coincides with a comparatively large interaction rate, our analysis highlights a considerable degeneracy between the mixed dark matter parameters and the neutrino mass scale.

We examine the maximum possible strength of the global 21-cm absorption dip on the Cosmic Background Radiation at high-redshift caused by the atomic intergalactic medium, when the Lyman-$\alpha$ coupling is maximum, assuming no exotic cooling mechanisms from interactions with dark matter. This maximum absorption is limited by three inevitable factors that need to be accounted for: $(a)$ heating by energy transferred from the Cosmic Background Radiation to the hydrogen atoms via 21-cm transitions, dubbed as 21-cm heating; $(b)$ Ly$\alpha$ heating by scatterings of Ly$\alpha$ photons from the first stars; $(c)$ the impact of the expected density fluctuations in the intergalactic gas in standard Cold Dark Matter theory, which reduces the mean 21-cm absorption signal. Inclusion of this third novel effect reduces the maximum global 21-cm absorption by $\sim 10\%$. Overall, the three effects studied here reduce the 21-cm global absorption by $\sim 20\%$ at $z \simeq 17$.

Coupled dark matter-dark energy scenarios are modeled via a dimensionless parameter $\xi$, which controls the strength of their interaction. While this coupling is commonly assumed to be constant, there is no underlying physical law or symmetry that forbids a time-dependent $\xi$ parameter. The most general and complete interacting scenarios between the two dark sectors should therefore allow for such a possibility, and it is the main purpose of this study to constrain two possible and well-motivated coupled cosmologies by means of the most recent and accurate early and late-time universe observations. We find that CMB data alone prefers $\xi(z) >0$ and therefore a smaller amount of dark matter, alleviating some crucial and well-known cosmological data tensions. An objective assessment of the Bayesian evidence for the coupled models explored here shows no particular preference for the presence of a dynamical dark sector coupling.

The discoveries of a number of binary black hole mergers by LIGO and VIRGO has reinvigorated the interest that primordial black holes (PBHs) of tens of solar masses could contribute non-negligibly to the dark matter energy density. Should even a small population of PBHs with masses $\gtrsim \mathcal{O}(M_\odot)$ exist, they could profoundly impact the properties of the intergalactic medium and provide insight into novel processes at work in the early Universe. We demonstrate here that observations of the 21cm transition in neutral hydrogen during the epochs of reionization and cosmic dawn will likely provide one of the most stringent tests of solar mass PBHs. In the context of 21cm cosmology, PBHs give rise to three distinct observable effects: ${\textit{(i)}}$ the modification to the primordial power spectrum (and thus also the halo mass function) induced by Poisson noise, ${\textit{(ii)}}$ a uniform heating and ionization of the intergalactic medium via X-rays produced during accretion, and ${\textit{(iii)}}$ a local modification to the temperature and density of the ambient medium surrounding isolated PBHs. Using a four-parameter astrophysical model, we show that experiments like SKA and HERA could potentially improve upon existing constraints derived using observations of the cosmic microwave background by more than one order of magnitude.

Mixed damping is a physical effect that occurs when a heavy species is coupled to a relativistic fluid which is itself free streaming. As a cross-case between collisional damping and free-streaming, it is crucial in the context of neutrino-dark matter interactions. In this work, we establish the parameter space relevant for mixed damping, and we derive an analytical approximation for the evolution of dark matter perturbations in the mixed damping regime to illustrate the physical processes responsible for the suppression of cosmological perturbations. Although extended Boltzmann codes implementing neutrino-dark matter scattering terms automatically include mixed damping, this effect has not been systematically studied. In order to obtain reliable numerical results, it is mandatory to reconsider several aspects of neutrino-dark matter interactions, such as the initial conditions, the ultra-relativistic fluid approximation and high order multiple moments in the neutrino distribution. Such a precise treatment ensures the correct assessment of the relevance of mixed damping in neutrino-dark matter interactions.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

We propose here a \emph{novel} method which singles out the \emph{a priori} unavoidable dependence on the underlying cosmological model when extracting parameter constraints, providing robust limits which only depend on the considered dataset. Interestingly, when dealing with several possible cosmologies and interpreting the Bayesian preference in terms of the Gaussian statistical evidence, the preferred model is much less favored than when only two cases are compared. As a working example, we apply our approach to the cosmological neutrino mass bounds, which play a fundamental role not only in establishing the contribution of relic neutrinos to the dark matter of the Universe, but also in the planning of future experimental searches of the neutrino character and of the neutrino mass ordering.

Dark matter interactions with massless or very light Standard Model particles, as photons or neutrinos, may lead to a suppression of the matter power spectrum at small scales and of the number of low mass haloes. Bounds on the dark matter scattering cross section with light degrees of freedom in such interacting dark matter (IDM) scenarios have been obtained from e.g. early time cosmic microwave background physics and large scale structure observations. Here we scrutinize dark matter microphysics in light of the claimed 21 cm EDGES 78 MHz absorption signal. IDM is expected to delay the 21 cm absorption features due to collisional damping effects. We identify the astrophysical conditions under which the existing constraints on the dark matter scattering cross section could be largely improved due to the IDM imprint on the 21 cm signal, providing also an explicit comparison to the WDM scenario.

The ordering of the neutrino masses is a crucial input for a deep understanding of flavor physics, and its determination may provide the key to establish the relationship among the lepton masses and mixings and their analogous properties in the quark sector. The extraction of the neutrino mass ordering is a data-driven field expected to evolve very rapidly in the next decade. In this review, we both analyze the present status and describe the physics of subsequent prospects. Firstly, the different current available tools to measure the neutrino mass ordering are described. Namely, reactor, long-baseline (accelerator and atmospheric) neutrino beams, laboratory searches for beta and neutrinoless double beta decays and observations of the cosmic background radiation and the large scale structure of the universe are carefully reviewed. Secondly, the results from an up-to-date comprehensive global fit are reported: the Bayesian analysis to the 2018 publicly available oscillation and cosmological data sets provides \emph{strong} evidence for the normal neutrino mass ordering versus the inverted scenario, with a significance of 3.5 standard deviations. This preference for the normal neutrino mass ordering is mostly due to neutrino oscillation measurements. Finally, we shall also emphasize the future perspectives for unveiling the neutrino mass ordering. In this regard, apart from describing the expectations from the aforementioned probes, we also focus on those arising from alternative and novel methods, as 21~cm cosmology, core-collapse supernova neutrinos and the direct detection of relic neutrinos.

We explore cosmological constraints on the sum of the three active neutrino masses $M_{\nu}$ in the context of dynamical dark energy (DDE) models with equation of state (EoS) parametrized as a function of redshift $z$ by $w(z)=w_0+w_a\,z/(1+z)$, and satisfying $w(z)\geq-1$ for all $z$. We perform a Bayesian analysis and show that, within these models, the bounds on $M_{\nu}$ \textit{do not degrade} with respect to those obtained in the $\Lambda$CDM case; in fact the bounds are slightly tighter, despite the enlarged parameter space. We explain our results based on the observation that, for fixed choices of $w_0\,,w_a$ such that $w(z)\geq-1$ (but not $w=-1$ for all $z$), the upper limit on $M_{\nu}$ is tighter than the $\Lambda$CDM limit because of the well-known degeneracy between $w$ and $M_{\nu}$. The Bayesian analysis we have carried out then integrates over the possible values of $w_0$-$w_a$ such that $w(z)\geq-1$, all of which correspond to tighter limits on $M_{\nu}$ than the $\Lambda$CDM limit. We find a 95\% confidence level (C.L.) upper bound of $M_{\nu}<0.13\,\mathrm{eV}$. This bound can be compared with $M_{\nu}<0.16\,\mathrm{eV}$ at 95\%~C.L., obtained within the $\Lambda$CDM model, and $M_{\nu}<0.41\,\mathrm{eV}$ at 95\%~C.L., obtained in a DDE model with arbitrary EoS (which allows values of $w < -1$). Contrary to the results derived for DDE models with arbitrary EoS, we find that a dark energy component with $w(z)\geq-1$ is unable to alleviate the tension between high-redshift observables and direct measurements of the Hubble constant $H_0$. Finally, in light of the results of this analysis, we also discuss the implications for DDE models of a possible determination of the neutrino mass hierarchy by laboratory searches. (abstract abridged)

The elastic scattering between dark matter particles and radiation represents an attractive possibility to solve a number of discrepancies between observations and standard cold dark matter predictions, as the induced collisional damping would imply a suppression of small-scale structures. We consider this scenario and confront it with measurements of the ionization history of the Universe at several redshifts and with recent estimates of the counts of Milky Way satellite galaxies. We derive a conservative upper bound on the dark matter-photon elastic scattering cross section of $\sigma_{\gamma \rm{DM}} < 8 \times 10^{-10} \, \sigma_T \, \left(m_{\rm DM}/{\rm GeV}\right)$ at $95\%$~CL, about one order of magnitude tighter than previous {constraints from satellite number counts}. Due to the strong degeneracies with astrophysical parameters, the bound on the dark matter-photon scattering cross section derived here is driven by the estimate of the number of Milky Way satellite galaxies. Finally, we also argue that future 21~cm probes could help in disentangling among possible non-cold dark matter candidates, such as interacting and warm dark matter scenarios. Let us emphasize that bounds of similar magnitude to the ones obtained here could be also derived for models with dark matter-neutrino interactions and would be as constraining as the tightest limits on such scenarios.

We examine the results from the Experiment to Detect the Global Epoch of Reionization Signature (EDGES), which has recently claimed the detection of a strong absorption in the 21 cm hyperfine transition line of neutral hydrogen, at redshifts demarcating the early stages of star formation. More concretely, we study the compatibility of the shape of the EDGES absorption profile, centered at a redshift of $z \sim 17.2$, with measurements of the reionization optical depth, the Gunn-Peterson optical depth, and Lyman-$\alpha$ emission from star-forming galaxies, for a variety of possible reionization models within the standard $\Lambda$CDM framework (that is, a Universe with a cosmological constant and cold dark matter). When, conservatively, we only try to accommodate the location of the absorption dip, we identify a region in the parameter space of the astrophysical parameters that successfully explains all of the aforementioned observations. However, one of the most abnormal features of the EDGES measurement is the absorption amplitude, which is roughly a factor of two larger than the maximum allowed value in the $\Lambda$CDM framework. We point out that the simple considered astrophysical models that produce the largest absorption amplitudes are unable to explain the depth of the dip and of reproducing the observed shape of the absorption profile.

In warm dark matter scenarios structure formation is suppressed on small scales with respect to the cold dark matter case, reducing the number of low-mass halos and the fraction of ionized gas at high redshifts and thus, delaying reionization. This has an impact on the ionization history of the Universe and measurements of the optical depth to reionization, of the evolution of the global fraction of ionized gas and of the thermal history of the intergalactic medium, can be used to set constraints on the mass of the dark matter particle. However, the suppression of the fraction of ionized medium in these scenarios can be partly compensated by varying other parameters, as the ionization efficiency or the minimum mass for which halos can host star-forming galaxies. Here we use different data sets regarding the ionization and thermal histories of the Universe and, taking into account the degeneracies from several astrophysical parameters, we obtain a lower bound on the mass of thermal warm dark matter candidates of $m_X > 1.3$ keV, or $m_s > 5.5$ keV for the case of sterile neutrinos non-resonantly produced in the early Universe, both at 90\% confidence level.

The simple $m^2\phi^2$ potential as an inflationary model is coming under increasing tension with limits on the tensor-to-scalar ratio $r$ and measurements of the scalar spectral index $n_s$. Cubic Galileon interactions in the context of the Horndeski action can potentially reconcile the observables. However, we show that this cannot be achieved with only a constant Galileon mass scale because the interactions turn off too slowly, leading also to gradient instabilities after inflation ends. Allowing for a more rapid transition can reconcile the observables but moderately breaks the slow-roll approximation leading to a relatively large and negative running of the tilt $\alpha_s$ that can be of order $n_s-1$. We show that the observables on CMB and large scale structure scales can be predicted accurately using the optimized slow-roll approach instead of the traditional slow-roll expansion. Upper limits on $|\alpha_s|$ place a lower bound of $r\gtrsim 0.005$ and conversely a given $r$ places a lower bound on $|\alpha_s|$, both of which are potentially observable with next generation CMB and large scale structure surveys.

A deep understanding of the Epoch of Reionization is still missing in our knowledge of the universe. While future probes will allow us to test the precise evolution of the free electron fraction from redshifts between $z\simeq 6$ and $z\simeq 20$, at present one could ask what kind of reionization processes are allowed by present Cosmic Microwave Background temperature and polarization measurements. An early contribution to reionization could imply a departure from the standard picture where star formation determines the reionization onset. BBy considering a broad class of possible reionization parameterizations, we find that current data do not require an early reionization component in our universe and that only one marginal class of models, based on a particular realization of reionization, may point to that. In addition, the frequentist Akaike Information Criterion (AIC) provides strong evidence against alternative reionization histories, favoring the most simple reionization scenario, which describes reionization by means of only one (constant) reionization optical depth $\tau$.

We present a Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and CMB observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and c) find that there is a weak-to-moderate preference for logarithmic priors. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass. We also show that the limits on the total neutrino mass $\sum m_\nu$ can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

We use present cosmological observations and forecasts of future experiments to illustrate the power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling potential deviations from the standard $\Lambda$CDM scenario. To quantify this statement, we focus on an extension of $\Lambda$CDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular and matter power spectra. After finding that future CMB experiments (such as COrE+) will not significantly improve the constraints set by the Planck satellite, we show that the next generation of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative cosmologies and even have the potential to make a discovery. Typically we find that DESI would be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that until now have only been accessible via $N$-body simulations.

Using some of the latest cosmological datasets publicly available, we derive the strongest bounds in the literature on the sum of the three active neutrino masses, $M_\nu$, within the assumption of a background flat $\Lambda$CDM cosmology. In the most conservative scheme, combining Planck cosmic microwave background (CMB) temperature anisotropies and baryon acoustic oscillations (BAO) data, as well as the up-to-date constraint on the optical depth to reionization ($\tau$), the tightest $95\%$ confidence level (C.L.) upper bound we find is $M_\nu<0.151$~eV. The addition of Planck high-$\ell$ polarization data, which however might still be contaminated by systematics, further tightens the bound to $M_\nu<0.118$~eV. A proper model comparison treatment shows that the two aforementioned combinations disfavor the IH at $\sim 64\%$~C.L. and $\sim 71\%$~C.L. respectively. In addition, we compare the constraining power of measurements of the full-shape galaxy power spectrum versus the BAO signature, from the BOSS survey. Even though the latest BOSS full shape measurements cover a larger volume and benefit from smaller error bars compared to previous similar measurements, the analysis method commonly adopted results in their constraining power still being less powerful than that of the extracted BAO signal. Our work uses only cosmological data; imposing the constraint $M_\nu>0.06\,{\rm eV}$ from oscillations data would raise the quoted upper bounds by ${\cal O}(0.1\sigma)$ and would not affect our conclusions.

We present a novel approach to derive constraints on neutrino masses from cosmological data, while taking into account our ignorance of the neutrino mass ordering. We derive constraints from a combination of current and future cosmological datasets on the total neutrino mass $M_\nu$ and on the mass fractions carried by each of the mass eigenstates, after marginalizing over the (unknown) neutrino mass ordering, either normal (NH) or inverted (IH). The bounds take therefore into account the uncertainty related to our ignorance of the mass hierarchy. This novel approach is carried out in the framework of Bayesian analysis of a typical hierarchical problem. In this context, the choice of the neutrino mass ordering is modeled via the discrete hyperparameter $h_{type}$. The preference for either the NH or the IH scenarios is then encoded in the posterior distribution of $h_{type}$ itself. Current CMB measurements assign equal odds to the two hierarchies, and are thus unable to distinguish between them. However, after the addition of BAO measurements, a weak preference for NH appears, with odds of 4:3 from Planck temperature and large-scale polarization in combination with BAO (3:2 if small-scale polarization is also included). Forecasts suggest that the combination of upcoming CMB (COrE) and BAO surveys (DESI) may determine the neutrino mass hierarchy at a high statistical significance if the mass is very close to the minimal value allowed by oscillations, as for NH and $M_\nu=0.06$ eV there is a 9:1 preference of NH vs IH. On the contrary, if $M_\nu$ is of the order of 0.1 eV or larger, even future cosmological observations will be inconclusive. The unbiased limit on $M_\nu$ we obtain with this innovative statistical strategy is crucial for ongoing and planned neutrinoless double beta decay searches.

In models with dark matter made of particles with keV masses, such as a sterile neutrino, small-scale density perturbations are suppressed, delaying the period at which the lowest mass galaxies are formed and therefore shifting the reionization processes to later epochs. In this study, focusing on Warm Dark Matter (WDM) with masses close to its present lower bound, i.e. around the $3$ keV region, we derive constraints from galaxy luminosy functions, the ionization history and the Gunn-Peterson effect. We show that even if star formation efficiency in the simulations is adjusted to match the observed UV galaxy luminosity functions in both CDM and WDM models, the full distribution of Gunn-Peterson optical depth retains the strong signature of delayed reionization in the WDM model. However, until the star formation and stellar feedback model used in modern galaxy formation simulations is constrained better, any conclusions on the nature of dark matter derived from reionization observables remain model-dependent.

Various particle physics models suggest that, besides the (nearly) cold dark matter that accounts for current observations, additional but sub-dominant dark relics might exist. These could be warm, hot, or even contribute as dark radiation. We present here a comprehensive study of two-component dark matter scenarios, where the first component is assumed to be cold, and the second is a non-cold thermal relic. Considering the cases where the non-cold dark matter species could be either a fermion or a boson, we derive consistent upper limits on the non-cold dark relic energy density for a very large range of velocity dispersions, covering the entire range from dark radiation to cold dark matter. To this end, we employ the latest Planck Cosmic Microwave Background data, the recent BOSS DR11 and other Baryon Acoustic Oscillation measurements, and also constraints on the number of Milky Way satellites, the latter of which provides a measure of the suppression of the matter power spectrum at the smallest scales due to the free-streaming of the non-cold dark matter component. We present the results on the fraction $f_{\rm ncdm}$ of non-cold dark matter with respect to the total dark matter for different ranges of the non-cold dark matter masses. We find that the 2$\sigma$ limits for non-cold dark matter particles with masses in the range 1--10 keV are $f_{\rm ncdm}\leq0.29$ (0.23) for fermions (bosons), and for masses in the 10--100 keV range they are $f_{\rm ncdm}\leq0.43$ (0.45), respectively.

The answer is Yes! We indeed find that interacting dark energy can alleviate the current tension on the value of the Hubble constant $H_0$ between the Cosmic Microwave Background anisotropies constraints obtained from the Planck satellite and the recent direct measurements reported by Riess et al. 2016. The combination of these two datasets points towards an evidence for a non-zero dark matter-dark energy coupling $\xi$ at more than two standard deviations, with $\xi=-0.26_{-0.12}^{+0.16}$ at $95\%$ CL. However the $H_0$ tension is better solved when the equation of state of the interacting dark energy component is allowed to freely vary, with a phantom-like equation of state $w=-1.184\pm0.064$ (at $68 \%$ CL), ruling out the pure cosmological constant case, $w=-1$, again at more than two standard deviations. When Planck data are combined with external datasets, as BAO, JLA Supernovae Ia luminosity distances, cosmic shear or lensing data, we find good consistency with the cosmological constant scenario and no compelling evidence for a dark matter-dark energy coupling.

We explore an extended cosmological scenario where the dark matter is an admixture of cold and additional non-cold species. The mass and temperature of the non-cold dark matter particles are extracted from a number of cosmological measurements. Among others, we consider tomographic weak lensing data and Milky Way dwarf satellite galaxy counts. We also study the potential of these scenarios in alleviating the existing tensions between local measurements and Cosmic Microwave Background (CMB) estimates of the $S_8$ parameter, with $S_8=\sigma_8\sqrt{\Omega_m}$, and of the Hubble constant $H_0$. In principle, a sub-dominant, non-cold dark matter particle with a mass $m_X\sim$~keV, could achieve the goals above. However, the preferred ranges for its temperature and its mass are different when extracted from weak lensing observations and from Milky Way dwarf satellite galaxy counts, since these two measurements require suppressions of the matter power spectrum at different scales. Therefore, solving simultaneously the CMB-weak lensing tensions and the small scale crisis in the standard cold dark matter picture via only one non-cold dark matter component seems to be challenging.

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.

Current measurements of the temperature and polarization anisotropy power spectra of the Cosmic Microwave Background (CMB) seem to indicate that the naive expectation for the slow-roll hierarchy within the most simple inflationary paradigm may not be respected in nature. We show that a primordial power spectra with localized features could in principle give rise to the observed slow-roll anarchy when fitted to a featureless power spectrum. Future CMB missions have the key to disentangle among the two possible paradigms and firmly establish the slow-roll mechanism as the responsible one for the inflationary period in the early universe. From a model comparison perspective, and assuming that nature has chosen a featureless primordial power spectrum, we find that, while with mock Planck data there is only weak evidence against a model with localized features, upcoming CMB measurements may provide strong evidence against such a non-standard primordial power spectrum.

After four years of data taking, the IceCube neutrino telescope has detected 54 high-energy starting events (HESE, or contained-vertex events) with deposited energies above 20TeV. They represent the first ever detection of high-energy extraterrestrial neutrinos and therefore, the first step in neutrino astronomy. In order to study the energy, flavor and isotropy of the astrophysical neutrino flux arriving at Earth, we perform different analyses of two different deposited energy intervals, [10 TeV $-$ 10 PeV] and [60 TeV $-$ 10 PeV]. We first consider an isotropic unbroken power-law spectrum and constrain its shape, normalization and flavor composition. Our results are in agreement with the preliminary IceCube results, although we obtain a slightly softer spectrum. We also find that current data are not sensitive to a possible neutrino-antineutrino asymmetry in the astrophysical flux. Then, we show that although a two-component power-law model leads to a slightly better fit, it does not represent a significant improvement with respect to a single power-law flux. Finally, we analyze the possible existence of a North-South asymmetry, hinted at by the combination of the HESE sample with the through-going muon data. If only using HESE data, the scarce statistics from the northern hemisphere does not allow us to reach any conclusive answer, which indicates that the HESE sample alone is not driving the potential North-South asymmetry.

We study the impact of assumptions about neutrino properties on the estimation of inflationary parameters from cosmological data, with a specific focus on the allowed contours in the $n_s/r$ plane. We study the following neutrino properties: (i) the total neutrino mass $ M_\nu =\sum_i m_i$; (ii) the number of relativistic degrees of freedom $N_{eff}$; and (iii) the neutrino hierarchy: whereas previous literature assumed 3 degenerate neutrino masses or two massless neutrino species (that do not match neutrino oscillation data), we study the cases of normal and inverted hierarchy. Our basic result is that these three neutrino properties induce $< 1 \sigma$ shift of the probability contours in the $n_s/r$ plane with both current or upcoming data. We find that the choice of neutrino hierarchy has a negligible impact. However, the minimal cutoff on the total neutrino mass $M_{\nu,{min}}=0 $ that accompanies previous works using the degenerate hierarchy does introduce biases in the $n_s/r$ plane and should be replaced by $M_{\nu,min}= 0.059$ eV as required by oscillation data. Using current CMB data from Planck and Bicep/Keck, marginalizing over $ M_\nu$ and over $r$ can lead to a shift in the mean value of $n_s$ of $\sim0.3\sigma$ towards lower values. However, once BAO measurements are included, the standard contours in the $n_s/r$ plane are basically reproduced. Larger shifts of the contours in the $n_s/r$ plane (up to 0.8$\sigma$) arise for nonstandard values of $N_{eff}$. We also provide forecasts for the future CMB experiments COrE and Stage-IV and show that the incomplete knowledge of neutrino properties, taken into account by a marginalization over $M_\nu$, could induce a shift of $\sim0.4\sigma$ towards lower values in the determination of $n_s$ (or a $\sim 0.8\sigma$ shift if one marginalizes over $N_{eff}$). Comparison to specific inflationary models is shown.

Recent Cosmic Microwave Background (CMB) temperature and polarization anisotropy measurements from the Planck mission have significantly improved previous constraints on the neutrino masses as well as the bounds on extended models with massless or massive sterile neutrino states. However, due to parameter degeneracies, additional low redshift priors are mandatory in order to sharpen the CMB neutrino bounds. We explore here the role of different priors on low redshift quantities, such as the Hubble constant, the cluster mass bias, and the reionization optical depth $\tau$. Concerning current priors on the Hubble constant and the cluster mass bias, the bounds on the neutrino parameters may differ appreciably depending on the choices adopted in the analyses. With regard to future improvements in the priors on the reionization optical depth, a value of $\tau=0.05\pm 0.01$, motivated by astrophysical estimates of the reionization redshift, would lead to $\sum m_\nu<0.0926$~eV at $90\%$~CL, when combining the full \textit{Planck} measurements, Baryon Acoustic Oscillation and Planck clusters data, thereby opening the window to unravel the neutrino mass hierarchy with existing cosmological probes.

In this paper, we extend an alternative, phenomenological approach to inflation by means of an equation of state and a sound speed, both of them functions of the number of $e$-folds and three phenomenological parameters. This approach captures a number of possible inflationary models, including those with non-canonical kinetic terms or scale-dependent non-gaussianities. We perform Markov Chain Monte Carlo analyses using the latest cosmological publicly available measurements, which include Cosmic Microwave Background (CMB) data from the Planck satellite. In the context of the phenomenological description studied here, the sound speed of the primordial curvature perturbation is found to be $c_s> 0.2$ at $95\%$~CL for a number of $e$-folds $N=60$. Within this parametrization, we discard scale invariance with a significance of about $10\sigma$, and the running of the spectral index is constrained as $\alpha_s=-0.62\,^{+0.08}_{-0.09} \times 10^{-3}$. The limit on the tensor-to-scalar ratio is $r<0.003$ at $95\%$~CL from CMB data alone. The maximum amplitude of the equilateral non-gaussianity that we obtain, $|f^{\text{equil}}_{\text{NL}}|\sim 2$, is much smaller than the current Planck mission errors. Future high-redshift, all-sky surveys could reach the required accuracy on equilateral non-gaussianities to allow for additional and independent tests of the parametrization explored here.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The simplest inflationary models predict a primordial power spectrum (PPS) of the curvature fluctuations that can be described by a power-law function that is nearly scale-invariant. It has been shown, however, that the low-multipole spectrum of the CMB anisotropies may hint the presence of some features in the shape of the scalar PPS, which could deviate from its canonical power-law form. We study the possible degeneracies of this non-standard PPS with the neutrino anisotropies, the neutrino masses, the effective number of relativistic species and a sterile neutrino or a thermal axion mass. The limits on these additional parameters are less constraining in a model with a non-standard PPS when only including the temperature auto-correlation spectrum measurements in the data analyses. The inclusion of the polarization spectra noticeably helps in reducing the degeneracies, leading to results that typically show no deviation from the $\Lambda$CDM model with a standard power-law PPS.

The most recent measurements of the temperature and low-multipole polarization anisotropies of the Cosmic Microwave Background (CMB) from the Planck satellite, when combined with galaxy clustering data from the Baryon Oscillation Spectroscopic Survey (BOSS) in the form of the full shape of the power spectrum, and with Baryon Acoustic Oscillation measurements, provide a $95\%$ confidence level (CL) upper bound on the sum of the three active neutrinos $\sum m _\nu< 0.183$ eV, among the tightest neutrino mass bounds in the literature, to date, when the same datasets are taken into account. This very same data combination is able to set, at $\sim70\%$ CL, an upper limit on $\sum m _\nu$ of $0.0968$ eV, a value that approximately corresponds to the minimal mass expected in the inverted neutrino mass hierarchy scenario. If high-multipole polarization data from Planck is also considered, the $95\%$ CL upper bound is tightened to $\sum m _\nu< 0.176$ eV. Further improvements are obtained by considering recent measurements of the Hubble parameter. These limits are obtained assuming a specific non-degenerate neutrino mass spectrum; they slightly worsen when considering other degenerate neutrino mass schemes. Current cosmological data, therefore, start to be mildly sensitive to the neutrino mass ordering. Low-redshift quantities, such as the Hubble constant or the reionization optical depth, play a very important role when setting the neutrino mass constraints. We also comment on the eventual shifts in the cosmological bounds on $\sum m_\nu$ when possible variations in the former two quantities are addressed.

Future dedicated radio interferometers, including HERA and SKA, are very promising tools that aim to study the epoch of reionization and beyond via measurements of the 21 cm signal from neutral hydrogen. Dark matter (DM) annihilations into charged particles change the thermal history of the Universe and, as a consequence, affect the 21 cm signal. Accurately predicting the effect of DM strongly relies on the modeling of annihilations inside halos. In this work, we use up-to-date computations of the energy deposition rates by the products from DM annihilations, a proper treatment of the contribution from DM annihilations in halos, as well as values of the annihilation cross section allowed by the most recent cosmological measurements from the Planck satellite. Given current uncertainties on the description of the astrophysical processes driving the epochs of reionization, X-ray heating and Lyman-$\alpha$ pumping, we find that disentangling DM signatures from purely astrophysical effects, related to early-time star formation processes or late-time galaxy X-ray emissions, will be a challenging task. We conclude that only annihilations of DM particles with masses of $\sim100$ MeV, could leave an unambiguous imprint on the 21 cm signal and, in particular, on the 21 cm power spectrum. This is in contrast to previous, more optimistic results in the literature, which have claimed that strong signatures might also be present even for much higher DM masses. Additional measurements of the 21 cm signal at different cosmic epochs will be crucial in order to break the strong parameter degeneracies between DM annihilations and astrophysical effects and undoubtedly single out a DM imprint for masses different from $\sim100$ MeV.

The value of the tensor-to-scalar ratio $r$ in the region allowed by the latest $Planck$ 2015 measurements can be associated to a large variety of inflationary models. We discuss here the potential of future Cosmic Microwave Background cosmological observations in disentangling among the possible theoretical scenarios allowed by our analyses of current $Planck$ temperature and polarization data. Rather than focusing only on $r$, we focus as well on the running of the primordial power spectrum, $\alpha_s$ and the running of thereof, $\beta_s$. Our Fisher matrix method benefits from a detailed and realistic appraisal of the expected foregrounds. Future cosmological probes, as the COrE mission, may be able to reach an unprecedented accuracy in the extraction of $\beta_s$ and rule out the most favoured inflationary models.

The simplest models of inflation predict small non-gaussianities and a featureless power spectrum. However, there exist a large number of well-motivated theoretical scenarios in which large non-gaussianties could be generated. In general, in these scenarios the primordial power spectrum will deviate from its standard power law shape. We study, in a model-independent manner, the constraints from future large scale structure surveys on the local non-gaussianity parameter $f_{\rm NL}$ when the standard power law assumption for the primordial power spectrum is relaxed. If the analyses are restricted to the large scale-dependent bias induced in the linear matter power spectrum by non-gaussianites, the errors on the $f_{\rm NL}$ parameter could be increased by $60\%$ when exploiting data from the future DESI survey, if dealing with only one possible dark matter tracer. In the same context, a nontrivial bias $|\delta f_{\rm NL}| \sim 2.5$ could be induced if future data are fitted to the wrong primordial power spectrum. Combining all the possible DESI objects slightly ameliorates the problem, as the forecasted errors on $f_{\rm NL}$ would be degraded by $40\%$ when relaxing the assumptions concerning the primordial power spectrum shape. Also the shift on the non-gaussianity parameter is reduced in this case, $|\delta f_{\rm NL}| \sim 1.6$. The addition of Cosmic Microwave Background priors ensure robust future $f_{\rm NL}$ bounds, as the forecasted errors obtained including these measurements are almost independent on the primordial power spectrum features, and $|\delta f_{\rm NL}| \sim 0.2$, close to the standard single-field slow-roll paradigm prediction.

Axions currently provide the most compelling solution to the strong CP problem. These particles may be copiously produced in the early universe, including via thermal processes. Therefore, relic axions constitute a hot dark matter component and their masses are strongly degenerate with those of the three active neutrinos, as they leave identical signatures in the different cosmological observables. In addition, thermal axions, while still relativistic states, also contribute to the relativistic degrees of freedom, parameterised via $N_{eff}$. We present the cosmological bounds on the relic axion and neutrino masses, exploiting the full Planck mission data, which include polarization measurements. In the mixed hot dark matter scenario explored here, we find the tightest and more robust constraint to date on the sum of the three active neutrino masses, $\sum m_\nu <0.136$ eV at $95\%$ CL, obtained in the well-known linear perturbation regime. The Planck Sunyaev-Zeldovich cluster number count data further tightens this bound, providing a $95\%$ CL upper limit of $\sum m_\nu <0.126$ eV in this very same mixed hot dark matter model, a value which is very close to the expectations in the inverted hierarchical neutrino mass scenario. Using this same combination of data sets we find the most stringent bound to date on the thermal axion mass, $m_a<0.529$ eV at $95\%$ CL.

A full energy and flavor-dependent analysis of the three-year high-energy IceCube neutrino events is presented. By means of multidimensional fits, we derive the current preferred values of the high-energy neutrino flavor ratios, the normalization and spectral index of the astrophysical fluxes, and the expected atmospheric background events, including a prompt component. A crucial assumption resides on the choice of the energy interval used for the analyses, which significantly biases the results. When restricting ourselves to the ~30 TeV - 3 PeV energy range, which contains all the observed IceCube events, we find that the inclusion of the spectral information improves the fit to the canonical flavor composition at Earth, (1:1:1), with respect to a single-energy bin analysis. Increasing both the minimum and the maximum deposited energies has dramatic effects on the reconstructed flavor ratios as well as on the spectral index. Imposing a higher threshold of 60 TeV yields a slightly harder spectrum by allowing a larger muon neutrino component, since above this energy most atmospheric tracklike events are effectively removed. Extending the high-energy cutoff to fully cover the Glashow resonance region leads to a softer spectrum and a preference for tau neutrino dominance, as none of the expected electron antineutrino induced showers have been observed so far. The lack of showers at energies above 2 PeV may point to a broken power-law neutrino spectrum. Future data may confirm or falsify whether or not the recently discovered high-energy neutrino fluxes and the long-standing detected cosmic rays have a common origin.

We present cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a "piecewise cubic Hermite interpolating polynomial" (PCHIP). When using Cosmic Microwave Background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of $\sigma_8$ and $\Omega_m$ from the 2013 Planck cluster catalogue on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favouring a $\sim 1$~eV axion mass, regardless of the model adopted for the primordial power spectrum.However, in general, such a preference disappears if the sum of the three active neutrino masses is also considered as a free parameter in our numerical analyses, due to the strong correlation between the masses of these two hot thermal relics.

We present an in-depth analysis of the flavour and spectral composition of the 36 high-energy neutrino events observed after three years of observation by the IceCube neutrino telescope. While known astrophysical sources of HE neutrinos are expected to produce a nearly $(1:1:1)$ flavour ratio (electron : muon : tau) of neutrinos at earth, we show that the best fits based on the events detected above $E_\nu \ge 28$ TeV do not necessarily support this hypothesis. Crucially, the energy range that is considered when analysing the HE neutrino data can have a profound impact on the conclusions. We highlight two intriguing puzzles: an apparent deficit of muon neutrinos, seen via a deficit of track-like events; and an absence of $\bar \nu_e$'s at high energy, seen as an absence of events near the Glashow resonance. We discuss possible explanations, including the misidentification of tracks as showers, and a broken power law, in analogy to the observed HE cosmic ray spectrum.

We examine the impact of a non-minimal coupling of the inflaton to the Ricci scalar, $\frac12 \xi R\phi^2$, on the inflationary predictions. Such a non-minimal coupling is expected to be present in the inflaton Lagrangian on fairly general grounds. As a case study, we focus on the simplest inflationary model governed by the potential $V\propto \phi^2$, using the latest combined 2015 analysis of Planck and BICEP2/Keck Array. We find that the presence of a coupling $\xi$ is favoured at a significance of $99\%$ CL, assuming that nature has chosen the potential $V\propto \phi^2$ to generate the primordial perturbations and a number of e-foldings $N=60$. Within the context of the same scenario, we find that the value of $\xi$ is different from zero at the $2\sigma$ level. When considering the cross-correlation polarization spectra from BICEP2/Keck Array and Planck, a value of $r=0.038_{-0.030}^{+0.039}$ is predicted in this particular non-minimally coupled scenario. Future cosmological observations may therefore test these values of $r$ and verify or falsify the non-minimally coupled model explored here.

In this work, we analyze two possible alternative and model-independent approaches to describe the inflationary period. The first one assumes a general equation of state during inflation due to Mukhanov, while the second one is based on the slow-roll hierarchy suggested by Hoffman and Turner. We find that, remarkably, the two approaches are equivalent from the observational viewpoint, as they single out the same areas in the parameter space, and agree with the inflationary attractors where successful inflation occurs. Rephrased in terms of the familiar picture of a slowly rolling, canonically normalized scalar field, the resulting inflaton excursions in these two approaches are almost identical. Furthermore, once the Galactic dust polarization data from Planck are included in the numerical fits, inflaton excursions can safely take sub-Planckian values.

We employ state-of-the art cosmological observables including supernova surveys and BAO information to provide constraints on the mass and mixing angle of a non-resonantly produced sterile neutrino species, showing that cosmology can effectively rule out sterile neutrinos which decay between BBN and the present day. The decoupling of an additional heavy neutrino species can modify the time dependence of the Universe's expansion between BBN and recombination and, in extreme cases, lead to an additional matter-dominated period; while this could naively lead to a younger Universe with a larger Hubble parameter, it could later be compensated by the extra radiation expected in the form of neutrinos from sterile decay. However, recombination-era observables including the Cosmic Microwave Background (CMB), the shift parameter $R_{CMB}$ and the sound horizon $r_s$ from Baryon Acoustic Oscillations (BAO) severely constrain this scenario. We self-consistently include the full time-evolution of the coupled sterile neutrino and standard model sectors in an MCMC, showing that if decay occurs after BBN, the sterile neutrino is essentially bounded by the constraint $\sin^2\theta \lesssim 0.026 (m_s/\mathrm{eV})^{-2}$.

The IceCube experiment has recently released 3 years of data of the first ever detected high-energy (>30 TeV) neutrinos, which are consistent with an extraterrestrial origin. In this talk, we compute the compatibility of the observed track-to-shower ratio with possible combinations of neutrino flavors with relative proportion (alpha_e:alpha_mu:alpha_tau). Although this observation is naively favored for the canonical (1:1:1) at Earth, once we consider the IceCube expectations for the atmospheric muon and neutrino backgrounds, this flavor combination presents some tension with data. We find that, for an astrophysical neutrino E_nu^{-2} energy spectrum, (1:1:1) at Earth is currently disfavored at 92% C.L. We discuss the trend of this result by comparing the results with the 2-year and 3-year data. We obtain the best-fit for (1:0:0) at Earth, which cannot be achieved from any flavor ratio at sources with averaged oscillations during propagation. Although it is not statistically significant at present, if confirmed, this result would suggest either a misunderstanding of the expected background events, or a misidentification of tracks as showers, or even more compellingly, some exotic physics which deviates from the standard scenario.

We revisit the cosmological viability of the Hu $\&$ Sawicki modified gravity scenario. The impact of such a modification on the different cosmological observables, including gravitational waves, is carefully described. The most recent cosmological data, as well as constraints on the relationship between the clustering parameter $\sigma_8$ and the current matter mass-energy density $\Omega_m$ from cluster number counts and weak lensing tomography, are considered in our numerical calculations. The strongest bound we find is $|f_{R0}| < 3.7 \times 10^{-6}$ at $95\%$~CL. Forthcoming cluster surveys covering 10,000 deg$^2$ in the sky, with galaxy surface densities of $\mathcal{O}(10)$~arcmin$^{-2}$ could improve the precision in the $\sigma_8$-$\Omega_m$ relationship, tightening the above constraint.

The IceCube experiment has recently reported the observation of 28 high-energy (> 30 TeV) neutrino events, separated into 21 showers and 7 muon tracks, consistent with an extraterrestrial origin. In this letter we compute the compatibility of such an observation with possible combinations of neutrino flavors with relative proportion (alpha_e:alpha_mu:alpha_tau). Although the 7:21 track-to-shower ratio is naively favored for the canonical (1:1:1) at Earth, this is not true once the atmospheric muon and neutrino backgrounds are properly accounted for. We find that, for an astrophysical neutrino E^(-2) energy spectrum, (1:1:1) at Earth is disfavored at 81% C.L. If this proportion does not change, 6 more years of data would be needed to exclude (1:1:1) at Earth at 3 sigma C.L. Indeed, with the recently-released 3-year data, that flavor composition is excluded at 92% C.L. The best-fit is obtained for (1:0:0) at Earth, which cannot be achieved from any flavor ratio at sources with averaged oscillations during propagation. If confirmed, this result would suggest either a misunderstanding of the expected background events, or a misidentification of tracks as showers, or even more compellingly, some exotic physics which deviates from the standard scenario.

Despite the great progress of current cosmological measurements, the nature of the dominant component of the universe, coined dark energy, is still an open question. Early Dark Energy is a possible candidate which may also alleviate some fine tuning issues of the standard paradigm. Using the latest available cosmological data, we find that the 95% CL upper bound on the early dark energy density parameter is $\Omega_{\textrm{eDE}}$. On the other hand, the dark energy component may be a stressed and inhomogeneous fluid. If this is the case, the effective sound speed and the viscosity parameters are unconstrained by current data. Future omniscope-like $21$cm surveys, combined with present CMB data, could be able to distinguish between standard quintessence scenarios from other possible models with $2\sigma$ significance, assuming a non-negligible early dark energy contribution. The precision achieved on the $\Omega_{\textrm{eDE}}$ parameter from these $21$ cm probes could be below $\mathcal{O} (10\%)$.

We present a one per cent measurement of the cosmic distance scale from the detections of the baryon acoustic oscillations in the clustering of galaxies from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). Our results come from the Data Release 11 (DR11) sample, containing nearly one million galaxies and covering approximately $8\,500$ square degrees and the redshift range $0.2<z<0.7$. We also compare these results with those from the publicly released DR9 and DR10 samples. Assuming a concordance $\Lambda$CDM cosmological model, the DR11 sample covers a volume of 13\,Gpc${}^3$ and is the largest region of the Universe ever surveyed at this density. We measure the correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of over $7\,\sigma$ in both the correlation function and power spectrum. Fitting for the position of the acoustic features measures the distance relative to the sound horizon at the drag epoch, $r_d$, which has a value of $r_{d,{\rm fid}}=149.28\,$Mpc in our fiducial cosmology. We find $D_V=(1264\pm25\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.32$ and $D_V=(2056\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ at $z=0.57$. At 1.0 per cent, this latter measure is the most precise distance constraint ever obtained from a galaxy survey. Separating the clustering along and transverse to the line-of-sight yields measurements at $z=0.57$ of $D_A=(1421\pm20\,{\rm Mpc})(r_d/r_{d,{\rm fid}})$ and $H=(96.8\pm3.4\,{\rm km/s/Mpc})(r_{d,{\rm fid}}/r_d)$. Our measurements of the distance scale are in good agreement with previous BAO measurements and with the predictions from cosmic microwave background data for a spatially flat cold dark matter model with a cosmological constant.

Inflation is the leading theory to describe elegantly the initial conditions that led to structure formation in our universe. In this paper, we present a novel phenomenological fit to the Planck, WMAP polarisation (WP) and the BICEP2 datasets using an alternative parameterisation. Instead of starting from inflationary potentials and computing the inflationary observables, we use a phenomenological parameterisation due to Mukhanov, describing inflation by an effective equation-of-state, in terms of the number of e-folds and two phenomenological parameters $\alpha$ and $\beta$. Within such a parametrisation, which captures the different inflationary models in a model-independent way, the values of the scalar spectral index $n_s$, its running and the tensor-to-scalar ratio $r$ are predicted, given a set of parameters $(\alpha,\beta)$. We perform a Markov Chain Monte Carlo analysis of these parameters, and we show that the combined analysis of Planck and WP data favours the Starobinsky and Higgs inflation scenarios. Assuming that the BICEP2 signal is not entirely due to foregrounds, the addition of this last data set prefers instead the $\phi^2$ chaotic models. The constraint we get from Planck and WP data alone on the derived tensor-to-scalar ratio is $r<0.18$ at $95\%$~CL, value which is consistent with the one quoted from the BICEP2 collaboration analysis, $r = 0.16^{+0-06}_{-0.05}$, after foreground subtraction. This is not necessarily at odds with the $2\sigma$ tension found between Planck and BICEP2 measurements when analysing data in terms of the usual $n_s$ and $r$ parameters, given that the parameterisation used here includes, implicitly, a running spectral index.

We investigate the axion dark matter scenario (ADM), in which axions account for all of the dark matter in the Universe, in light of the most recent cosmological data. In particular, we use the Planck temperature data, complemented by WMAP E-polarization measurements, as well as the recent BICEP2 observations of B-modes. Baryon Acoustic Oscillation data, including those from the Baryon Oscillation Spectroscopic Survey, are also considered in the numerical analyses. We find that, in the minimal ADM scenario, the full dataset implies that the axion mass m_a = 82.2 pm 1.1 {\mu}eV (corresponding to the Peccei-Quinn symmetry being broken at a scale f_a = (7.54 pm 0.10)*10^10 GeV), or m_a = 76.6 pm 2.6 {\mu}eV (f_a = (8.08 pm 0.27)*10^10 GeV) when we allow for a non- standard effective number of relativistic species Neff . We also find a 2{\sigma} preference for Neff > 3.046. The limit on the sum of neutrino masses is Sum m_{\nu} < 0.25eV at 95% CL for Neff = 3.046, or Sum m_{\nu} < 0.47 eV when Neff is a free parameter. Considering extended scenarios where either the dark energy equation-of-state parameter w, the tensor spectral index nt or the running of the scalar index dns/dlnk are allowed to vary does not change significantly the axion mass-energy density constraints. However, in the case of the full dataset exploited here, there is a preference for a non-zero tensor index or scalar running, driven by the different tensor amplitudes implied by the Planck and BICEP2 observations.

We present up to date cosmological bounds on the sum of active neutrino masses as well as on extended cosmological scenarios with additional thermal relics, as thermal axions or sterile neutrino species. Our analyses consider all the current available cosmological data in the beginning of year 2014, including the very recent and most precise Baryon Acoustic Oscillation (BAO) measurements from the Baryon Oscillation Spectroscopic Survey. In the minimal three active neutrino scenario, we find Sum m_nu < 0.22 eV at 95% CL from the combination of CMB, BAO and Hubble Space Telescope measurements of the Hubble constant. A non zero value for the sum of the three active neutrino masses of about 0.3 eV is significantly favoured at more than 3 standard deviations when adding the constraints on sigma_8 and Omega_m from the Planck Cluster catalog on galaxy number counts. This preference for non zero thermal relic masses disappears almost completely in both the thermal axion and massive sterile neutrino schemes. Extra light species contribute to the effective number of relativistic degrees of freedom, parameterised via Neff. We found that when the recent detection of B mode polarization from the BICEP2 experiment is considered, an analysis of the combined CMB data in the framework of LCDM+r models gives Neff=4.00pm0.41, suggesting the presence of an extra relativistic relic at more than 95 % c.l. from CMB-only data.

We use the latest cosmic microwave background (CMB) observations to provide updated constraints on the dark matter lifetime as well as on p-wave suppressed annihilation cross sections in the 1 MeV to 1 TeV mass range. In contrast to scenarios with an s-wave dominated annihilation cross section, which mainly affect the CMB close to the last scattering surface, signatures associated with these scenarios essentially appear at low redshifts ($z \lesssim 50$) when structure began to form, and thus manifest at lower multipoles in the CMB power spectrum. We use data from Planck, WMAP9, SPT and ACT, as well as Lyman-$\alpha$ measurements of the matter temperature at $z \sim 4$ to set a 95 % confidence level lower bound on the dark matter lifetime of $\sim 4 \times 10^{25}$ s for $m_\chi = 100$ MeV. This bound becomes lower by an order of magnitude at $m_\chi = 1$ TeV due to inefficient energy deposition into the intergalactic medium. We also show that structure formation can enhance the effect of p--wave suppressed annihilation cross sections by many orders of magnitude with respect to the background cosmological rate, although even with this enhancement, CMB constraints are not yet strong enough to reach the thermal relic value of the cross section.

Axion-Like Particles (ALPs), if exist in nature, are expected to mix with photons in the presence of an external magnetic field. The energy range of photons which undergo strong mixing with ALPs depends on the ALP mass, on its coupling with photons as well as on the external magnetic field and particle density configurations. Recent observations of blazars by the Fermi Gamma-Ray Space Telescope in the 0.1-300 GeV energy range show a break in their spectra in the 1-10 GeV range. We have modeled this spectral feature for the flat-spectrum radio quasar 3C454.3 during its November 2010 outburst, assuming that a significant fraction of the gamma rays convert to ALPs in the large scale jet of this blazar. Using theoretically motivated models for the magnetic field and particle density configurations in the kiloparsec scale jet, outside the broad-line region, we find an ALP mass $m_a\sim (1-3)\cdot 10^{-7}$ eV and coupling $g_{a\gamma} \sim (1-3)\cdot 10^{-10}$ GeV$^{-1}$ after performing an illustrative statistical analysis of spectral data in four different epochs of emission. The precise values of $m_a$ and $g_{a\gamma}$ depend weakly on the assumed particle density configuration and are consistent with the current experimental bounds on these quantities. We apply this method and ALP parameters found from fitting 3C454.3 data to another flat-spectrum radio quasar PKS1222+216 (4C+21.35) data up to 400 GeV, as a consistency check, and found good fit. We find that the ALP-photon mixing effect on the GeV spectra may not be washed out for any reasonable estimate of the magnetic field in the intergalactic media.

We compute the bounds on the dark matter (DM) annihilation cross section using the most recent Cosmic Microwave Background measurements from WMAP9, SPT'11 and ACT'10. We consider DM with mass in the MeV-TeV range annihilating 100% into either an e+e- or a mu+mu- pair. We consider a realistic energy deposition model, which includes the dependence on the redshift, DM mass and annihilation channel. We exclude the canonical thermal relic abundance cross section (<sigma v> = 3 x 10^{-26} cm^3 s^{-1}) for DM masses below 30 GeV and 15 GeV for the e+e- and mu+mu- channels, respectively. A priori, DM annihilating in halos could also modify the reionization history of the Universe at late times. We implement a realistic halo model taken from results of state-of-the-art N-body simulations and consider a mixed reionization mechanism, consisting on reionization from DM as well as from first stars. We find that the constraints on DM annihilation remain unchanged, even when large uncertainties on the halo model parameters are considered.

We present new constraints on coupled dark energy from the recent measurements of the Cosmic Microwave Background Anisotropies from the Planck satellite mission. We found that a coupled dark energy model is fully compatible with the Planck measurements, deriving a weak bound on the dark matter-dark energy coupling parameter \xi=-0.49^{+0.19}_{-0.31} at 68% c.l.. Moreover if Planck data are fitted to a coupled dark energy scenario, the constraint on the Hubble constant is relaxed to H_0=72.1^{+3.2}_{-2.3} km/s/Mpc, solving the tension with the Hubble Space Telescope value. We show that a combined Planck+HST analysis provides significant evidence for coupled dark energy finding a non-zero value for the coupling parameter \xi, with -0.90< \xi <-0.22 at 95% c.l.. We also consider the combined constraints from the Planck data plus the BAO measurements of the 6dF Galaxy Survey, the Sloan Digital Sky Survey and the Baron Oscillation Spectroscopic Survey.

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.

We present here bounds on neutrino masses from the combination of recent Planck Cosmic Microwave Background measurements and galaxy clustering information from the Baryon Oscillation Spectroscopic Survey (BOSS), part of the Sloan Digital Sky Survey-III. We use the full shape of either the photometric angular clustering (Data Release 8) or the 3D spectroscopic clustering (Data Release 9) power spectrum in different cosmological scenarios. In the Lambda$CDM scenario, spectroscopic galaxy clustering measurements improve significantly the existing neutrino mass bounds from Planck data. We find sum m_nu< 0.39 eV at 95% confidence level for the combination of the 3D power spectrum with Planck CMB data (with lensing included) and Wilkinson Microwave Anisoptropy Probe 9-year polarization measurements. Therefore, robust neutrino mass constraints can be obtained without the addition of the prior on the Hubble constant from HST. In extended cosmological scenarios with a dark energy fluid or with non flat geometries, galaxy clustering measurements are essential to pin down the neutrino mass bounds, providing in the majority of cases better results than those obtained from the associated measurement of the Baryon Acoustic Oscillation scale only. In the presence of a freely varying (constant) dark energy equation of state, we find sum m_nu<0.49 eV at 95% confidence level for the combination of the 3D power spectrum with Planck CMB data (with lensing included) and Wilkinson Microwave Anisoptropy Probe 9-year polarization measurements. This same data combination in non flat geometries provides the neutrino mass bound sum m_nu<0.35 eV at 95% confidence level.

Recent Cosmic Microwave Background (CMB) results from the Planck satellite, combined with previous CMB data and Hubble constant measurements from the Hubble Space Telescope, provide a constraint on the effective number of relativistic degrees of freedom of Neff=3.62^{+0.50}_{-0.48} at 95% CL. These new measurements provide a unique opportunity to place limits on models containing relativistic species at the decoupling epoch. Here we review the bounds or the allowed parameter regions in sterile neutrino models, hadronic axion models as well as on extended dark sectors with additional light species based on the latest Planck CMB observations.

We present measurements of the angular diameter distance to and Hubble parameter at z=0.57 from the measurement of the baryon acoustic peak in the correlation of galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey. Our analysis is based on a sample from Data Release 9 of 264,283 galaxies over 3275 square degrees in the redshift range 0.43<z<0.70. We use two different methods to provide robust measurement of the acoustic peak position across and along the line of sight in order to measure the cosmological distance scale. We find D_A(0.57) = 1408 +/- 45 Mpc and H(0.57) = 92.9 +/- 7.8 km/s/Mpc for our fiducial value of the sound horizon. These results from the anisotropic fitting are fully consistent with the analysis of the spherically averaged acoustic peak position presented in Anderson et al, 2012. Our distance measurements are a close match to the predictions of the standard cosmological model featuring a cosmological constant and zero spatial curvature.

Recent Cosmic Microwave Background measurements at high multipoles from the South Pole Telescope and from the Atacama Cosmology Telescope seem to disagree in their conclusions for the neutrino and dark radiation properties. In this paper we set new bounds on the dark radiation and neutrino properties in different cosmological scenarios combining the ACT and SPT data with the nine year data release of the Wilkinson Microwave Anisotropy Probe (WMAP9), Baryon Acoustic Oscillation data, Hubble Telescope measurements of the Hubble constant, and Supernovae Ia luminosity distance data. In the standard three massive neutrino case, the two high multipole probes give similar results if Baryon Acoustic Oscillation data are removed from the analyses and Hubble Telescope measurements are also exploited. A similar result is obtained within a standard cosmology with Neff massless neutrinos, although in this case the agreement between these two measurements is also improved when considering simultaneously Baryon Acoustic Oscillation data and Hubble Space Telescope measurements. In the Neff massive neutrino case the two high multipole probes give very different results regardless of the external data sets used in the combined analyses. In the case in which a dark radiation background with unknown clustering properties is also considered, SPT data seem to exclude the standard value for the dark radiation viscosity cvis=1/3 at the 2 sigma CL, finding evidence for massive neutrinos only when combining SPT data with BAO measurements.

Dark matter particles gravitationally trapped inside the Sun may annihilate into Standard Model particles, producing a flux of neutrinos. The prospects of detecting these neutrinos in future multi-\kton{} neutrino detectors designed for other physics searches are explored here. We study the capabilities of a 34/100 \kton{} liquid argon detector and a 100 \kton{} magnetized iron calorimeter detector. These detectors are expected to determine the energy and the direction of the incoming neutrino with unprecedented precision allowing for tests of the dark matter nature at very low dark matter masses, in the range of 5-50 GeV. By suppressing the atmospheric background with angular cuts, these techniques would be sensitive to dark matter - nucleon spin dependent cross sections at the fb level, reaching down to a few ab for the most favorable annihilation channels and detector technology.

An extra dark radiation component can be present in the universe in the form of sterile neutrinos, axions or other very light degrees of freedom which may interact with the dark matter sector. We derive here the cosmological constraints on the dark radiation abundance, on its effective velocity and on its viscosity parameter from current data in dark radiation-dark matter coupled models. The cosmological bounds on the number of extra dark radiation species do not change significantly when considering interacting schemes. We also find that the constraints on the dark radiation effective velocity are degraded by an order of magnitude while the errors on the viscosity parameter are a factor of two larger when considering interacting scenarios. If future Cosmic Microwave Background data are analysed assuming a non interacting model but the dark radiation and the dark matter sectors interact in nature, the reconstructed values for the effective velocity and for the viscosity parameter will be shifted from their standard 1/3 expectation, namely ceff=0.34 (+0.006 -0.003) and cvis=0.29 (+0.002 -0.001) at 95% CL for the future COrE mission data.

Cosmological constraints on the sum of neutrino masses and on the effective number of neutrino species in standard and non-standard scenarios are computed using the most recent available cosmological data. Our cosmological data sets include the measurement of the Baryonic Acoustic Oscillation (BAO) feature in the Data Release 9 CMASS sample of the Baryon Oscillation Spectroscopic Survey (BOSS). We study in detail the different degeneracies among the parameters, as well as the impact of the different data sets used in the analyses. When considering bounds on the sum of the three active neutrino masses, the information in the BAO signal from galaxy clustering measurements is approximately equally powerful as the shape information from the matter power spectrum. The most stringent bound we find is sum m_nu<0.32 eV at 95 % CL. When non-standard neutrino scenarios with neff massless or massive neutrino species are examined, power spectrum shape measurements provide slightly better bounds than the BAO signal only, due to the breaking of parameter degeneracies. Recent BOSS data combined with CMB and Hubble Space Telescope measurements give neff=3.66^{+0.20 +0.73}_{-0.21 -0.69} in the massless neutrino scenario, and similar results are obtained in the massive case. The evidence for extra radiation neff>3 often claimed in the literature therefore remains at the 2 sigma level when considering up-to-date cosmological data sets. Measurements from the Wilkinson Microwave Anisotropy Probe combined with a prior on the Hubble parameter from the Hubble Space Telescope are very powerful in constraining either the sum of the three active neutrino masses or the number of massless neutrino species. If the former two parameters are allowed to freely vary, however, the bounds from the combination of these two cosmological probes get worse by an order of magnitude.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

Asymmetric Dark Matter (ADM) models invoke a particle-antiparticle asymmetry, similar to the one observed in the Baryon sector, to account for the Dark Matter (DM) abundance. Both asymmetries are usually generated by the same mechanism and generally related, thus predicting DM masses around 5 GeV in order to obtain the correct density. The main challenge for successful models is to ensure efficient annihilation of the thermally produced symmetric component of such a light DM candidate without violating constraints from collider or direct searches. A common way to overcome this involves a light mediator, into which DM can efficiently annihilate and which subsequently decays into Standard Model particles. Here we explore the scenario where the light mediator decays instead into lighter degrees of freedom in the dark sector that act as radiation in the early Universe. While this assumption makes indirect DM searches challenging, it leads to signals of extra radiation at BBN and CMB. Under certain conditions, precise measurements of the number of relativistic species, such as those expected from the Planck satellite, can provide information on the structure of the dark sector. We also discuss the constraints of the interactions between DM and Dark Radiation from their imprint in the matter power spectrum.

Recent cosmological data have provided evidence for a "dark" relativistic background at high statistical significance. Parameterized in terms of the number of relativistic degrees of freedom Neff, however, the current data seems to indicate a higher value than the one expected in the standard scenario based on three active neutrinos. This dark radiation component can be characterized not only by its abundance but also by its clustering properties, as its effective sound speed and its viscosity parameter. It is therefore crucial to study the correlations among the dark radiation properties and key cosmological parameters, as the dark energy equation of state or the running of the scalar spectral index, with current and future CMB data. We find that dark radiation with viscosity parameters different from their standard values may be misinterpreted as an evolving dark energy component or as a running spectral index in the power spectrum of primordial fluctuations.

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 measurements of galaxy clustering from the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey III (SDSS-III). These use the Data Release 9 (DR9) CMASS sample, which contains 264,283 massive galaxies covering 3275 square degrees with an effective redshift z=0.57 and redshift range 0.43 < z < 0.7. Assuming a concordance Lambda-CDM cosmological model, this sample covers an effective volume of 2.2 Gpc^3, and represents the largest sample of the Universe ever surveyed at this density, n = 3 x 10^-4 h^-3 Mpc^3. We measure the angle-averaged galaxy correlation function and power spectrum, including density-field reconstruction of the baryon acoustic oscillation (BAO) feature. The acoustic features are detected at a significance of 5\sigma in both the correlation function and power spectrum. Combining with the SDSS-II Luminous Red Galaxy Sample, the detection significance increases to 6.7\sigma. Fitting for the position of the acoustic features measures the distance to z=0.57 relative to the sound horizon DV /rs = 13.67 +/- 0.22 at z=0.57. Assuming a fiducial sound horizon of 153.19 Mpc, which matches cosmic microwave background constraints, this corresponds to a distance DV(z=0.57) = 2094 +/- 34 Mpc. At 1.7 per cent, this is the most precise distance constraint ever obtained from a galaxy survey. We place this result alongside previous BAO measurements in a cosmological distance ladder and find excellent agreement with the current supernova measurements. We use these distance measurements to constrain various cosmological models, finding continuing support for a flat Universe with a cosmological constant.

Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This {\it beat coupling} effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated {\it average} density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this {\it local average effect} and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only $\sim 10 %$ of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations. For the variances, we find excellent agreement with the analytic expressions for $k < 0.2 h$Mpc$^{-1}$ at $z=0.5$, while the correlation coefficients agree to beyond $k=0.4 h$Mpc$^{-1}$. As expected, the range of agreement increases towards higher redshift and decreases slightly towards $z=0$. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix. [abridged]

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

The Sloan Digital Sky Survey (SDSS) surveyed 14,555 square degrees, and delivered over a trillion pixels of imaging data. We present a study of galaxy clustering using 900,000 luminous galaxies with photometric redshifts, spanning between $z=0.45$ and $z=0.65$, constructed from the SDSS using methods described in Ross et al. (2011). This data-set spans 11,000 square degrees and probes a volume of $3h^{-3} \rm{Gpc}^3$, making it the largest volume ever used for galaxy clustering measurements. We present a novel treatment of the observational systematics and its applications to the clustering signals from the data set. In this paper, we measure the angular clustering using an optimal quadratic estimator at 4 redshift slices with an accuracy of ~15% with bin size of delta_l = 10 on scales of the Baryon Acoustic Oscillations (BAO) (at l~40-400). We derive cosmological constraints using the full-shape of the power-spectra. For a flat Lambda CDM model, when combined with Cosmic Microwave Background Wilkinson Microwave Anisotropy Probe 7 (WMAP7) and H_0 constraints from 600 Cepheids observed by HST, we find \Omega_\Lambda = 0.73 +/- 0.019 and H_0 to be 70.5 +/- 1.6 km/s/Mpc. For an open Lambda CDM model, when combined with WMAP7 + HST, we find $\Omega_K = 0.0035 +/- 0.0054, improved over WMAP7+HST alone by 40%. For a wCDM model, when combined with WMAP7+HST+SN, we find w = -1.071 +/- 0.078, and H_0 to be 71.3 +/- 1.7 km/s/Mpc, which is competitive with the latest large scale structure constraints from large spectroscopic surveys such as SDSS Data Release 7 (DR7) (Reid et al. 2010, Percival et al. 2010, Montesano et al. 2011) and WiggleZ (Blake et al. 2011). The SDSS-III Data Release 8 (SDSS-III DR8) Angular Clustering Data allows a wide range of investigations into the cosmological model, cosmic expansion (via BAO), Gaussianity of initial conditions and neutrino masses. (abridged)

We present neutrino mass bounds using 900,000 luminous galaxies with photometric redshifts measured from Sloan Digital Sky Survey III Data Release Eight (SDSS DR8). The galaxies have photometric redshifts between $z = 0.45$ and $z = 0.65$, and cover 10,000 square degrees and thus probe a volume of 3$h^{-3}$Gpc$^3$, enabling tight constraints to be derived on the amount of dark matter in the form of massive neutrinos. A new bound on the sum of neutrino masses $\sum m_\nu < 0.26$ eV, at 95% confidence level (CL), is obtained after combining our sample of galaxies, which we call "CMASS", with WMAP 7 year Cosmic Microwave Background (CMB) data and the most recent measurement of the Hubble parameter from the Hubble Space Telescope (HST). This constraint is obtained with a conservative multipole range choice of $30 < \ell < 200$ in order to minimize non-linearities, and a free bias parameter in each of the four redshift bins. We study the impact of assuming this linear galaxy bias model using mock catalogs, and find that this model causes a small ($\sim 1-1.5 \sigma$) bias in $\Omega_{\rm DM} h^2$. For this reason, we also quote neutrino bounds based on a conservative galaxy bias model containing additional, shot noise-like free parameters. In this conservative case, the bounds are significantly weakened, e.g. $\sum m_\nu < 0.36$ eV (95% confidence level) for WMAP+HST+CMASS ($\ell_{\rm max}=200$). We also study the dependence of the neutrino bound on multipole range ($\ell_{\rm max}=150$ vs $\ell_{\rm max}=200$) and on which combination of data sets is included as a prior. The addition of supernova and/or Baryon Acoustic Oscillation data does not significantly improve the neutrino mass bound once the HST prior is included. [abridged]

Sterile massive neutrinos are a natural extension of the Standard Model of elementary particles. The energy density of the extra sterile massive states affects cosmological measurements in an analogous way to that of active neutrino species. We perform here an analysis of current cosmological data and derive bounds on the masses of the active and the sterile neutrino states as well as on the number of sterile states. The so-called (3+2) models with three sub-eV active massive neutrinos plus two sub-eV massive sterile species is well within the 95% CL allowed regions when considering cosmological data only. If the two extra sterile states have thermal abundances at decoupling, Big Bang Nucleosynthesis bounds compromise the viability of (3+2) models. Forecasts from future cosmological data on the active and sterile neutrino parameters are also presented. Independent measurements of the neutrino mass from tritium beta decay experiments and of the Hubble constant could shed light on sub-eV massive sterile neutrino scenarios.

Cosmological measurements are affected by the energy density of both active and sterile massive neutrinos. We extend here a recent analysis of current cosmological data to non minimal cosmologies. Several possible scenarios are examined: a constant w \neq -1 dark energy equation of state, a non flat universe, a time varying dark energy component and coupled dark matter dark energy universes or modified gravity scenarios. When considering cosmological data only, (3+2) massive neutrino models with ~0.5 eV sterile species are allowed at 95% CL. This scenario has been shown to reconcile reactor, LSND and MiniBooNE positive signals with null results from other searches. Big Bang Nucleosynthesis bounds could compromise the viability of (3+2) models if the two sterile species are fully thermalized states at decoupling.

General relativistic corrections to the galaxy power spectrum appearing at the horizon scale, if neglected, may induce biases on the measured values of the cosmological parameters. In this paper, we study the impact of general relativistic effects on non standard cosmologies such as scenarios with a time dependent dark energy equation of state, with a coupling between the dark energy and the dark matter fluids or with non-Gaussianities. We then explore whether general relativistic corrections affect future constraints on cosmological parameters in the case of a constant dark energy equation of state and of non-Gaussianities. We find that relativistic corrections on the power spectrum are not expected to affect the foreseen errors on the cosmological parameters nor to induce large biases on them.

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

We present current and future constraints on the Hu and Sawicki modified gravity scenario. This model can reproduce a late time accelerated universe and evade solar system constraints. While current cosmological data still allows for distinctive deviations from the cosmological constant picture, future measurements of the growth of structure combined with Supernova Ia luminosity distance data will greatly improve present constraints.

Coupled cosmologies can predict values for the cosmological parameters at low redshifts which may differ substantially from the parameters values within non-interacting cosmologies. Therefore, low redshift probes, as the growth of structure and the dark matter distribution via galaxy and weak lensing surveys constitute a unique tool to constrain interacting dark sector models. We focus here on weak lensing forecasts from future Euclid and LSST-like surveys combined with the ongoing Planck cosmic microwave background experiment. We find that these future data could constrain the dimensionless coupling to be smaller than a few $\times 10^{-2}$. The coupling parameter $\xi$ is strongly degenerate with the cold dark matter energy density $\Omega_{c}h^2$ and the Hubble constant $H_0$.These degeneracies may cause important biases in the cosmological parameter values if in the universe there exists an interaction among the dark matter and dark energy sectors.

Photons couple to Axion-Like Particles (ALPs) or more generally to any pseudo Nambu-Goldstone boson in the presence of an external electromagnetic field. Mixing between photons and ALPs in the strong magnetic field of a Gamma-Ray Burst (GRB) jet during the prompt emission phase can leave observable imprints on the gamma-ray polarization and spectrum. Mixing in the intergalactic medium is not expected to modify these signatures for ALP mass > 10^(-14) eV and/or for < nG magnetic field. We show that the depletion of photons due to conversion to ALPs changes the linear degree of polarization from the values predicted by the synchrotron model of gamma ray emission. We also show that when the magnetic field orientation in the propagation region is perpendicular to the field orientation in the production region, the observed synchrotron spectrum becomes steeper than the theoretical prediction and as detected in a sizable fraction of GRB sample. Detection of the correlated polarization and spectral signatures from these steep-spectrum GRBs by gamma-ray polarimeters can be a very powerful probe to discover ALPs. Measurement of gamma-ray polarization from GRBs in general, with high statistics, can also be useful to search for ALPs.

We study a coupled dark energy-dark matter model in which the energy-momentum exchange is proportional to the Hubble expansion rate. The inclusion of its perturbation is required by gauge invariance. We derive the linear perturbation equations for the gauge invariant energy density contrast and velocity of the coupled fluids, and we determine the initial conditions. The latter turn out to be adiabatic for dark energy, when assuming adiabatic initial conditions for all the standard fluids. We perform a full Monte Carlo Markov Chain likelihood analysis of the model, using WMAP 7-year data.

We illustrate how recently improved low-redshift cosmological measurements can tighten constraints on neutrino properties. In particular we examine the impact of the assumed cosmological model on the constraints. We first consider the new HST H0 = 74.2 +/- 3.6 measurement by Riess et al. (2009) and the sigma8*(Omegam/0.25)^0.41 = 0.832 +/- 0.033 constraint from Rozo et al. (2009) derived from the SDSS maxBCG Cluster Catalog. In a Lambda CDM model and when combined with WMAP5 constraints, these low-redshift measurements constrain sum mnu<0.4 eV at the 95% confidence level. This bound does not relax when allowing for the running of the spectral index or for primordial tensor perturbations. When adding also Supernovae and BAO constraints, we obtain a 95% upper limit of sum mnu<0.3 eV. We test the sensitivity of the neutrino mass constraint to the assumed expansion history by both allowing a dark energy equation of state parameter w to vary, and by studying a model with coupling between dark energy and dark matter, which allows for variation in w, Omegak, and dark coupling strength xi. When combining CMB, H0, and the SDSS LRG halo power spectrum from Reid et al. 2009, we find that in this very general model, sum mnu < 0.51 eV with 95% confidence. If we allow the number of relativistic species Nrel to vary in a Lambda CDM model with sum mnu = 0, we find Nrel = 3.76^{+0.63}_{-0.68} (^{+1.38}_{-1.21}) for the 68% and 95% confidence intervals. We also report prior-independent constraints, which are in excellent agreement with the Bayesian constraints.

Cosmological analysis based on currently available observations are unable to rule out a sizeable coupling among the dark energy and dark matter fluids. We explore a variety of coupled dark matter-dark energy models, which satisfy cosmic microwave background constraints, in light of low redshift and near universe observations. We illustrate the phenomenology of different classes of dark coupling models, paying particular attention in distinguishing between effects that appear only on the expansion history and those that appear in the growth of structure. We find that while a broad class of dark coupling models are effectively models where general relativity (GR) is modified --and thus can be probed by a combination of tests for the expansion history and the growth of structure--, there is a class of dark coupling models where gravity is still GR, but the growth of perturbations is, in principle modified. While this effect is small in the specific models we have considered, one should bear in mind that an inconsistency between reconstructed expansion history and growth may not uniquely indicate deviations from GR. Our low redshift constraints arise from cosmic velocities, redshift space distortions and dark matter abundance in galaxy voids. We find that current data constrain the dimensionless coupling to be |xi|<0.2, but prospects from forthcoming data are for a significant improvement. Future, precise measurements of the Hubble constant, combined with high-precision constraints on the growth of structure, could provide the key to rule out dark coupling models which survive other tests. We shall exploit as well weak equivalence principle violation arguments, which have the potential to highly disfavour a broad family of coupled models.

Detection of high-energy (~> 100 MeV) gamma rays by the Fermi Large Area Telescope (LAT) from a nova in the symbiotic binary system V407 Cygni has opened possibility of high-energy neutrino detection from this type of sources. Thermonuclear explosion on the white dwarf surface sets off a nova shell in motion that expands and slows down in a dense surrounding medium provided by the red giant companion. Particles are accelerated in the shocks of the shell, and interact with surrounding medium to produce observed gamma rays. We show that proton-proton interaction, which is most likely responsible for producing gamma rays via neutral pion decay, produces ~> 0.1 GeV neutrinos that can be detected by the current and future experiments at ~> 10 GeV.

We study a coupled quintessence model in which the interaction with the dark matter sector is a function of the quintessence potential. Such a coupling can arise from a field dependent mass term for the dark matter field. The dynamical analysis of a standard quintessence potential coupled with the interaction explored here shows that the system possesses a late time accelerated attractor. In light of these results, we perform a fit to the most recent Supernovae Ia, Cosmic Microwave Background and Baryon Acoustic Oscillation data sets. Constraints arising from weak equivalence principle violation arguments are also discussed.

Determination of whether the Harrison--Zel'dovich spectrum for primordial scalar perturbations is consistent with observations is sensitive to assumptions about the reionization scenario. In light of this result, we revisit constraints on inflationary models using more general reionization scenarios. While the bounds on the tensor-to-scalar ratio are largely unmodified, when different reionization schemes are addressed, hybrid models are back into the inflationary game. In the general reionization picture, we reconstruct both the shape and amplitude of the inflaton potential. We find a broader spectrum of potential shapes when relaxing the simple reionization restriction. An upper limit of $10^{16}$ GeV to the amplitude of the potential is found, regardless of the assumptions on the reionization history.

We study the impact of the cosmological parameters uncertainties on the measurements of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While this is not expected to be an issue for the standard LCDM model, it may not be the case for more general models that modify the large-scale shape of the power spectrum. We consider the so-called local non-Gaussianity model and forecasts from planned surveys, alone and combined with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast the correlations among $f_{\rm NL}$ and the running of the spectral index $\alpha_s$, the dark energy equation of state $w$, the effective sound speed of dark energy perturbations $c^2_s$, the total mass of massive neutrinos $M_\nu=\sum m_\nu$, and the number of extra relativistic degrees of freedom $N_\nu^{rel}$. Neglecting CMB information on $f_{\rm NL}$ and scales $k > 0.03 h$/Mpc, we find that, if $N_\nu^{\rm rel}$ is assumed to be known, the uncertainty on cosmological parameters increases the error on $f_{\rm NL}$ by 10 to 30% depending on the survey. Thus the $f_{\rm NL}$ constraint is remarkable robust to cosmological model uncertainties. On the other hand, if $N_\nu^{\rm rel}$ is simultaneously constrained from the data, the $f_{\rm NL}$ error increases by $\sim 80%$. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local form with a marginalized 1--$\sigma$ error of the order $\Delta f_{\rm NL} \sim 2-5$, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment.

Cosmic Microwave Background satellite missions as the on-going Planck experiment are expected to provide the strongest constraints on a wide set of cosmological parameters. Those constraints, however, could be weakened when the assumption of a cosmological constant as the dark energy component is removed. Here we show that it will indeed be the case when there exists a coupling among the dark energy and the dark matter fluids. In particular, the expected errors on key parameters as the cold dark matter density and the angular diameter distance at decoupling are significantly larger when a dark coupling is introduced. We show that it will be the case also for future satellite missions as EPIC, unless CMB lensing extraction is performed.

Inflation predicts primordial scalar perturbations with a nearly scale-invariant spectrum and a spectral index approximately unity (the Harrison--Zel'dovich (HZ) spectrum). The first important step for inflationary cosmology is to check the consistency of the HZ primordial spectrum with current observations. Recent analyses have claimed that a HZ primordial spectrum is excluded at more than 99% c.l.. Here we show that the HZ spectrum is only marginally disfavored if one considers a more general reionization scenario. Data from the Planck mission will settle the issue.

Geometrical tests such as the combination of the Hubble parameter H(z) and the angular diameter distance d_A(z) can, in principle, break the degeneracy between the dark energy equation of state parameter w(z), and the spatial curvature Omega_k in a direct, model-independent way. In practice, constraints on these quantities achievable from realistic experiments, such as those to be provided by Baryon Acoustic Oscillation (BAO) galaxy surveys in combination with CMB data, can resolve the cosmic confusion between the dark energy equation of state parameter and curvature only statistically and within a parameterized model for w(z). Combining measurements of both H(z) and d_A(z) up to sufficiently high redshifts around z = 2 and employing a parameterization of the redshift evolution of the dark energy equation of state are the keys to resolve the w(z)-Omega_k degeneracy.

A class of well-behaved modified gravity models with long enough matter domination epoch and a late-time accelerated expansion is confronted with SNIa, CMB, SDSS, BAO and H(z) galaxy ages data, as well as current measurements of the linear growth of structure. We show that the combination of geometrical probes and growth data exploited here allows to rule out f(R) gravity models, in particular, the logarithmic of curvature model. We also apply solar system tests to the models in agreement with the cosmological data. We find that the exponential of the inverse of the curvature model satisfies all the observational tests considered and we derive the allowed range of parameters. Current data still allows for small deviations of Einstein gravity. Future, high precision growth data, in combination with expansion history data, will be able to distinguish tiny modifications of standard gravity from the LambdaCDM model.

The two dark sectors of the universe - dark matter and dark energy - may interact with each other. Background and linear density perturbation evolution equations are developed for a generic coupling. We then establish the general conditions necessary to obtain models free from early time non-adiabatic instabilities. As an application, we consider a viable universe in which the interaction strength is proportional to the dark energy density. The scenario does not exhibit "phantom crossing" and is free from instabilities, including early ones. A sizeable interaction strength is compatible with combined WMAP, HST, SN, LSS and H(z) data. Neutrino mass and/or cosmic curvature are allowed to be larger than in non-interacting models. Our analysis sheds light as well on unstable scenarios previously proposed.

Protons accelerated in the internal shocks of a long duration gamma ray burst can escape the fireball as cosmic rays by converting to neutrons. Hadronic interactions of these neutrons inside a stellar wind bubble created by the progenitor star will produce TeV gamma rays via neutral meson decay and synchrotron radiation by charged pion-decay electrons in the wind magnetic field. Such gamma rays should be observable from nearby gamma ray bursts by currently running and upcoming ground-based detectors.

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

We show that the measurements of 10 GeV atmospheric neutrinos by an upcoming array of densely packed phototubes buried deep inside the IceCube detector at the South Pole can be used to determine the neutrino mass hierarchy for values of sin^2(2theta13) close to the present bound, if the hierarchy is normal. These results are obtained for an exposure of 100 Mton years and systematic uncertainties up to 10%.

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

A suppression in the spectrum of ultrahigh-energy (UHE, >= 10^{18} eV) neutrinos will be present in extra-dimensional scenarios, due to enhanced neutrino-antineutrino annihilation processes with the supernova relic neutrinos. In the n>4 scenario, being n the number of extra dimensions, neutrinos can not be responsible for the highest energy events observed in the UHE cosmic ray spectrum. A direct implication of these extra-dimensional interactions would be the absence of UHE neutrinos in ongoing and future neutrino telescopes.

Relic thermal axions could play the role of an extra hot dark matter component in cosmological structure formation theories. By combining the most recent observational data we improve previous cosmological bounds on the axion mass m_a in the so-called hadronic axion window. We obtain a limit on the axion mass m_a < 0.42eV at the 95% c.l. (m_a < 0.72eV at the 99% c.l.). A novel aspect of the analysis presented here is the inclusion of massive neutrinos and how they may affect the bound on the axion mass. If neutrino masses belong to an inverted hierarchy scheme, for example, the above constraint is improved to m_a < 0.38eV at the 95% c.l. (m_a < 0.67eV at the 99% c.l.). Future data from experiments as CAST will provide a direct test of the cosmological bound.

High-energy neutrinos are expected to be produced in a vareity of astrophysical sources as well as in optically thick hidden sources. We explore the matter-induced oscillation effects on emitted neutrino fluxes of three different flavors from the latter class. We use the ratio of electron and tau induced showers to muon tracks, in upcoming neutrino telescopes, as the principal observable in our analysis. This ratio depends on the neutrino energy, density profile of the sources and on the oscillation parameters. The largely unknown flux normalization drops out of our calculation and only affects the statistics. For the current knowledge of the oscillation parameters we find that the matter-induced effects are non-negligible and the enhancement of the ratio from its vacuum value takes place in an energy range where the neutrino telescopes are the most sensitive. Quantifying the effect would be useful to learn about the astrophysics of the sources as well as the oscillation parameters. If the neutrino telescopes mostly detect diffuse neutrinos without identifying their sources, then any deviation of the measured flux ratios from the vacuum expectation values would be most naturally explained by a large population of hidden sources for which matter-induced neutrino oscillation effects are important.

We discuss the present bounds on the relativistic energy density in the Universe parameterized in terms of the effective number of neutrinos N using the most recent cosmological data on Cosmic Microwave Background (CMB) temperature anisotropies and polarization, Large Scale galaxy clustering from the Sloan Digital Sky Survey (SDSS) and 2dF, luminosity distances of type Ia Supernovae, Lyman-alpha absorption clouds (Ly-alpha), the Baryonic Acoustic Oscillations (BAO) detected in the Luminous Red Galaxies of the SDSS and finally, Big Bang Nucleosynthesis (BBN) predictions for 4He and Deuterium abundances. We find N= 5.2+2.7-2.2 from CMB and Large Scale Structure data, while adding Ly-alpha and BAO we obtain N= 4.6+1.6-1.5 at 95 % c.l.. These results show some tension with the standard value N=3.046 as well as with the BBN range N= 3.1+1.4-1.2 at 95 % c.l., though the discrepancy is slightly below the 2-sigma level. In general, considering a smaller set of data weakens the constraints on N. We emphasize the impact of an improved upper limit (or measurement) of the primordial value of 3He abundance in clarifying the issue of whether the value of N at early (BBN) and more recent epochs coincide.

The black hole at the center of the galaxy is a powerful lens for supernova neutrinos. In the very special circumstance of a supernova near the extended line of sight from Earth to the galactic center, lensing could dramatically enhance the neutrino flux at Earth and stretch the neutrino pulse.

We show that the current accelerated expansion of the Universe can be explained without resorting to dark energy. Models of generalized modified gravity, with inverse powers of the curvature can have late time accelerating attractors without conflicting with solar system experiments. We have solved the Friedman equations for the full dynamical range of the evolution of the Universe. This allows us to perform a detailed analysis of Supernovae data in the context of such models that results in an excellent fit. Hence, inverse curvature gravity models represent an example of phenomenologically viable models in which the current acceleration of the Universe is driven by curvature instead of dark energy. If we further include constraints on the current expansion rate of the Universe from the Hubble Space Telescope and on the age of the Universe from globular clusters, we obtain that the matter content of the Universe is 0.07 <= omega_m <= 0.21 (95% Confidence). Hence the inverse curvature gravity models considered can not explain the dynamics of the Universe just with a baryonic matter component.

If the equation of state for "dark energy" varies periodically, the expansion of the Universe may have undergone alternating eras of acceleration and deceleration. We examine a specific form that survives existing observational tests, does not single out the present state of the Universe as exceptional, and suggests a future much like the matter-dominated past: a smooth expansion without a final inflationary epoch.