Clustering of the large scale structure provides complementary information to the measurements of the cosmic microwave background anisotropies through power spectrum and bispectrum of density perturbations. Extracting the bispectrum information, however, is more challenging than it is from the power spectrum due to the complex models and the computational cost to measure the signal and its covariance. To overcome these problems, we adopt a proxy statistic, skew spectrum which is a cross-spectrum of the density field and its quadratic field. By applying a large smoothing filter to the density field, we show the theory fits the simulations very well. With the spectra and their full covariance estimated from $N$-body simulations as our "mock" Universe, we perform a global fits for the cosmological parameters. The results show that adding skew spectrum to power spectrum the $1\sigma$ marginalized errors for parameters $ b_1^2A_s, n_s$ and $f_{\rm NL}^{\rm loc}$ are reduced by $31\%, 22\%, 44\%$, respectively. This is the answer to the question posed in the title and indicates that the skew spectrum will be a fast and effective method to access complementary information to that enclosed in the power spectrum measurements, especially for the forthcoming generation of wide-field galaxy surveys.

We forecast the benefits induced by adding the bispectrum anisotropic signal to the standard, two- and three-point, clustering statistics analysis. In particular, we forecast cosmological parameter constraints including the bispectrum higher multipoles terms together with the galaxy power spectrum (monopole plus quadrupole) and isotropic bispectrum (monopole) data vectors. To do so, an analytical covariance matrix model is presented. This template is carefully calibrated on well-known terms of a numerical covariance matrix estimated from a set of simulated galaxy catalogues. After testing the calibration using the power spectrum and isotropic bispectrum measurements from the same set of simulations, we extend the covariance modelling to the galaxy bispectrum higher multipoles. Using this covariance matrix we proceed to perform cosmological parameter inference using a suitably generated mock data vector. Including the bispectrum mutipoles up to the hexadecapole, yields 1-D $68\%$ credible regions for the set of parameters $(b_1,b_2,f,\sigma_8,f_\mathrm{NL},\alpha_\perp, \alpha_\parallel)$ tighter by a factor of 30$\%$ on average for $k_\mathrm{max}=0.09\,h$/Mpc, significantly reducing at the same time the degeneracies present in the posterior distribution.

The concept of blind analysis, a key procedure to remove the human-based systematic error called confirmation bias, has long been an integral part of data analysis in many research areas. In cosmology, blind analysis is recently making its entrance, as the field progresses into a fully fledged high-precision science. The credibility, reliability and robustness of results from future sky-surveys will dramatically increase if the effect of confirmation bias is kept under control by using an appropriate blinding procedure. Here, we present a catalog-level blinding scheme for galaxy clustering data apt to be used in future spectroscopic galaxy surveys. We shift the individual galaxy positions along the line of sight based on 1) a geometric shift mimicking the Alcock-Paczynski effect and 2) a perturbative shift akin to redshift-space distortions. This procedure has several advantages. After combining the two steps above, it is almost impossible to accidentally unblind. The procedure induces a shift in cosmological parameters without changing the galaxies' angular positions, hence without interfering with the effects of angular systematics. Since the method is applied at catalog level, there is no need to adopt several blinding schemes tuned to different summary statistics, likelihood choices or types of analyses. By testing the method on mock catalogs and the BOSS DR12 catalog we demonstrate its performance in blinding galaxy clustering data for relevant cosmological parameters sensitive to the background expansion rate and the growth rate of structures.

Are the stellar-mass merging binary black holes, recently detected by their gravitational wave signal, of stellar or primordial origin? Answering this question will have profound implications for our understanding of the Universe, including the nature of dark matter, the early Universe and stellar evolution. We build on the idea that the clustering properties of merging binary black holes can provide information about binary formation mechanisms and origin. The cross-correlation of galaxy with gravitational wave catalogues carries information about whether black hole mergers trace more closely the distribution of dark matter -- indicative of primordial origin -- or that of stars harboured in luminous and massive galaxies -- indicative of a stellar origin. We forecast the detectability of such signal for several forthcoming and future gravitational wave interferometers and galaxy surveys, including, for the first time in such analyses, an accurate modelling for the different merger rates, lensing magnification and other general relativistic effects. Our results show that forthcoming experiments could allow us to test most of the parameter space of the still viable models investigated, and shed more light on the issue of binary black hole origin and evolution.

In the era of precision cosmology, establishing the correct magnitude of statistical errors in cosmological parameters is of crucial importance. However, widely used approximations in galaxy surveys analyses can lead to parameter uncertainties that are grossly mis-estimated. These approximations can be introduced at three different levels: in the form of the likelihood, in the theoretical modelling of the observable and in the numerical computation of the observable. Their consequences are important both in data analysis through e.g., Markov Chain Monte Carlo parameter inference, and when survey instrument and strategy are designed and their constraining power on cosmological parameters is forecasted, for instance using Fisher matrix analyses. In this work, considering the galaxy angular power spectrum as the target observable, we report one example for each of such three categories and derive their impact on the estimated parameters uncertainties and degeneracies. We show that there are cases where these commonly used approximations lead, perhaps counter-intuitively, to unacceptably large mis-estimates of parameters errors and correlations. Furthermore, we stress how these approximations might even spoil the benefits of the nascent multi-tracer and multi-messenger cosmology. Hence we recommend that the type of analysis presented here should be repeated for every approximation adopted in survey design or data analysis, to quantify how it may affect the results. To this aim, we have developed \texttt{Multi\_CLASS}, a new extension of \texttt{CLASS} that includes the angular power spectrum for multiple (galaxy and other tracers such as gravitational waves) populations. The public release of \texttt{Multi\_CLASS} is associated with this paper.

Cosmological parameter estimation from forthcoming experiments promise to reach much greater precision than current constraints. As statistical errors shrink, the required control over systematic errors increases. Therefore, models or approximations that were sufficiently accurate so far, may introduce significant systematic biases in the parameter best-fit values and jeopardize the robustness of cosmological analyses. We present a general expression to estimate a priori the systematic error introduced in parameter inference due to the use of insufficiently good approximations in the computation of the observable of interest or the assumption of an incorrect underlying model. Although this methodology can be applied to measurements of any scientific field, we illustrate its power by studying the effect of modeling the angular galaxy power spectrum incorrectly. We also introduce Multi_CLASS, a new, public modification of the Boltzmann code CLASS, which includes the possibility to compute angular cross-power spectra for two different tracers. We find that significant biases in most of the cosmological parameters are introduced if one assumes the Limber approximation or neglects lensing magnification in modern galaxy survey analyses, and the effect is in general larger for the multi-tracer case, especially for the parameter controlling primordial non-Gaussianity of the local type, $f_{\rm NL}$.

During the red giant phase, stars loose mass at the highest rate since birth. The mass-loss rate is not fixed, but varies from star-to-star by up to 5%, resulting in variations of the star's luminosity at the tip of the red giant branch (TRGB). Also, most stars, during this phase, engulf part of their planetary system, including their gas giant planets. Gas giant planet masses range between 0.1 to 2% of the host star mass. The engulfing of their gas giants planets can modify their luminosity at the TRGB, i.e. the point at which the He-core degeneracy is removed. We show that the increase in mass of the star by the engulfing of the gas giant planets only modifies the luminosity of a star at the TRGB by less than 0.1%, while metallicity can modify the luminosity of a star at the TRGB by up to 0.5%. However, the increase in turbulence of the convective envelope of the star, i.e., modification of the mixing length, has a more dramatic effect, on the star's luminosity, which we estimate could be as large as 5%. The effect is always in the direction to increase the turbulence and thus the mixing length which turns into a systematic decrease of the luminosity of the star at the TRGB. We find that the star-to-star variation of the mass-loss rate will dominate the variations in the luminosity of the TRGB with a contribution at the 5% level. If the star-to-star variation is driven by environmental effects --as it is reasonable to assume--, the same effects can potentially create an environmentally-driven mean effect on the luminosity of the tip of the red giant branch of a galaxy. Finally, we touch upon how to infer the frequency, and identify the engulfment, of exoplanets in low-metallicity RGB stars through high resolution spectroscopy as well as how to quantify mass loss rate distributions from the morphology of the horizontal branch.

The evolution of differential ages of passive galaxies at different redshifts (cosmic chronometers) has been proved to be a method potentially able to constrain the Hubble parameter in a cosmology-independent way, but the systematic uncertainties must be carefully evaluated. In this paper, we compute the contribution to the full covariance matrix of systematic uncertainties due to the choice of initial mass function, stellar library, and metallicity, exploring a variety of stellar population synthesis models; the effect of a residual young component and star formation history was discussed elsewhere. Through simulations in the redshift range 0<z<1.5 we find that the choice of the stellar population synthesis model dominates the total error budget on H(z), with contributions at ~4.5% level, discarding the most discordant model. The contribution due to the choice of initial mass function is <0.5%, while that due to the stellar library is ~6.6% on average. In parallel, we also assess the impact of an uncertainty in the determination of the stellar metallicity, finding that an error of ~10% (5%) on the stellar metallicity propagates to a 9% (4%) error on H(z). These results are used to provide the combined contribution of these systematic effects on the error budget. For current H(z) measurements, where the uncertainties due to metallicity and star formation history were already included, we show that, using the more modern stellar libraries, the additional systematic uncertainty is between 5.4% (at z=0.2) and 2.3% (at z=1.5). To reach the goal of keeping the systematic error budget below the 1% level we discuss that efforts needed to obtain higher resolution and signal-to-noise spectra and improvements in the modeling of stellar population synthesis.

Due to instrumental limitations, the nature of which vary from case to case, spectroscopic galaxy redshift surveys usually do not collect redshifts for all galaxies in the population of potential targets. Especially problematic is the entanglement between this incompleteness and the true cosmological signal, arising from the fact that the proportion of successful observations is typically lower in regions with higher density of galaxies. The result is a fictitious suppression of the galaxy clustering that, if not properly corrected, can impact severely on cosmological-parameter inference. Recent developments in the field have shown that an unbiased estimate of the 2-point correlation function in the presence of missing observations can be obtained by weighting each pair by its inverse probability of being targeted. In this work we expand on the concept of probability weights by developing a more mature statistical formalism, which provides us with a deeper understanding of their fundamental properties. We take advantage of this novel perspective to handle the problem of estimating the inverse probability, specifically, we discuss how to efficiently determine the weights from a finite set of realisations of the targeting and how to model exactly the resulting sampling effects. This allows us to derive an inverse-probability-based estimator for the power spectrum, which is the main result of this work, but also to improve robustness and computational efficiency of the already existing configuration-space estimator. Finally, we propose a strategy to further extend the concept of inverse probability, providing examples of how traditional missing-observation countermeasures can be included in this more general picture. The effectiveness of the different models and weighting schemes discussed in this work is demonstrated using realisations of an idealised simple survey strategy.

Constraints on dark matter annihilation or decay offer unique insights into the nature of dark matter. We illustrate how surveys dedicated to detect the highly redshifted 21 cm signal from the dark ages will offer a new window into properties of particle dark matter. The 21 cm intensity mapping signal and its fluctuations are sensitive to energy injection from annihilation or decay of long-lived particles in a way that is complementary to other probes. We present forecasted constraints from forthcoming and next-generation radio surveys. We show that, while SKA might be capable of a detection for some cases, the most promising opportunity to detect dark matter in the 21 cm intensity mapping signal is with a futuristic radio array on the lunar far-side, with the potential to detect a signal many orders of magnitude weaker than current or maximal constraints from other probes.

We study the clustering properties of dark matter halos in real- and redshift-space in cosmologies with massless and massive neutrinos through a large set of state-of-the-art N-body simulations. We provide quick and easy-to-use prescriptions for the halo bias on linear and mildly non-linear scales, both in real and redshift space, which are valid also for massive neutrinos cosmologies. Finally we present a halo bias emulator,$\textbf{BE-HaPPY}$, calibrated on the N-body simulations, which is fast enough to be used in the standard Markov Chain Monte Carlo approach to cosmological inference. For a fiducial standard $\Lambda$CDM cosmology $\textbf{BE-HaPPY}$ provides percent or sub-percent accuracy on the scales of interest (linear and well into the mildly non-linear regime), meeting therefore for the halo-bias the accuracy requirements for the analysis of next-generation large--scale structure surveys.

In the model where Primordial Black Holes (PBHs) form from large primordial curvature (C) perturbations, i.e., CPBHs, constraints on PBH abundance provide in principle constraints on the primordial curvature power spectrum. This connection however depends necessarily on the details of PBH formation mechanism. In this paper we provide, for the first time, constraints on the primordial curvature power spectrum from the latest limits on PBH abundance, taking into account all the steps from gravitational collapse in real space to PBH formation. In particular, we use results from numerical relativity simulations and peak theory to study the conditions for PBH formation for a range of perturbation shapes, including non-linearities, perturbation profile and a careful treatment of smoothing and filtering scales. We then obtain updated PBH formation conditions and translate that into primordial spectrum constraints for a wide range of shapes and abundances. These updated constraints cover a range of scales not probed by other cosmological observables. Our results show that the correct and accurate modelling of non-linearities, filtering and typical perturbation profile, is crucial for deriving meaningful cosmological implications.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Neutral hydrogen pervades the infant Universe, and its redshifted 21-cm signal allows one to chart the Universe. This signal allows one to probe astrophysical processes such as the formation of the first stars, galaxies, (super)massive black holes and enrichment of the pristine gas from z~6 to z~30, as well as fundamental physics related to gravity, dark matter, dark energy and particle physics at redshifts beyond that. As one enters the Dark Ages (z>30), the Universe becomes pristine. Ground-based low-frequency radio telescopes aim to detect the spatial fluctuations of the 21-cm signal. Complementary, global 21-cm experiments aim to measure the sky-averaged 21-cm signal. Escaping RFI and the ionosphere has motivated space-based missions, such as the Dutch-Chinese NCLE instrument (currently in lunar L2), the proposed US-driven lunar or space-based instruments DAPPER and FARSIDE, the lunar-orbit interferometer DSL (China), and PRATUSH (India). To push beyond the current z~25 frontier, though, and measure both the global and spatial fluctuations (power-spectra/tomography) of the 21-cm signal, low-frequency (1-100MHz; BW~50MHz; z>13) space-based interferometers with vast scalable collecting areas (1-10-100 km2), large filling factors (~1) and large fields-of-view (4pi sr.) are needed over a mission lifetime of >5 years. In this ESA White Paper, we argue for the development of new technologies enabling interferometers to be deployed, in space (e.g. Earth-Sun L2) or in the lunar vicinity (e.g. surface, orbit or Earth-Moon L2), to target this 21-cm signal. This places them in a stable environment beyond the reach of most RFI from Earth and its ionospheric corruptions, enabling them to probe the Dark Ages as well as the Cosmic Dawn, and allowing one to investigate new (astro)physics that is inaccessible in any other way in the coming decades. [Abridged]

The standard cosmological model successfully describes many observations from widely different epochs of the Universe, from primordial nucleosynthesis to the accelerating expansion of the present day. However, as the basic cosmological parameters of the model are being determined with increasing and unprecedented precision, it is not guaranteed that the same model will fit more precise observations from widely different cosmic epochs. Discrepancies developing between observations at early and late cosmological time may require an expansion of the standard model, and may lead to the discovery of new physics. The workshop "Tensions between the Early and the Late Universe" was held at the Kavli Institute for Theoretical Physics on July 15-17 2019 (More details of the workshop (including on-line presentations) are given at the website: https://www.kitp.ucsb.edu/activities/enervac-c19) to evaluate increasing evidence for these discrepancies, primarily in the value of the Hubble constant as well as ideas recently proposed to explain this tension. Multiple new observational results for the Hubble constant were presented in the time frame of the workshop using different probes: Cepheids, strong lensing time delays, tip of the red giant branch (TRGB), megamasers, Oxygen-rich Miras and surface brightness fluctuations (SBF) resulting in a set of six new ones in the last several months. Here we present the summary plot of the meeting that shows combining any three independent approaches to measure H$_0$ in the late universe yields tension with the early Universe values between 4.0$\sigma$ and 5.8$\sigma$. This shows that the discrepancy does not appear to be dependent on the use of any one method, team, or source. Theoretical ideas to explain the discrepancy focused on new physics in the decade of expansion preceding recombination as the most plausible. This is a brief summary of the workshop.

We propose a method to probe the homogeneity of a general universe, without assuming symmetry. We show that isotropy can be tested at remote locations on the past lightcone by comparing the line-of-sight and transverse expansion rates, using the time dependence of the polarization of Cosmic Microwave Background photons that have been inverse-Compton scattered by the hot gas in massive clusters of galaxies. This probes a combination of remote transverse and parallel components of the expansion rate of the metric, and we may use radial baryon acoustic oscillations or cosmic clocks to measure the parallel expansion rate. Thus we can test remote isotropy, which is a key requirement of a homogeneous universe. We provide explicit formulas that connect observables and properties of the metric.

The ages of the oldest stellar objects in our galaxy provide an independent test of the current cosmological model as they give a lower limit to the age of the Universe. Recent accurate parallaxes by the Gaia space mission, accurate measurements of the metallicity of stars, via individual elemental abundances, and advances in the modelling of stellar evolution, provide new, higher-precision age estimates of the oldest stellar populations in the galaxy: globular clusters and very-low-metallicity stars. The constraints on the age of the Universe, $t_U$, so obtained are determined from the local Universe and at late time. It is well known that local and early-Universe determinations of another cosmological parameter closely related to the age of the Universe, the Hubble constant $H_0$, show a $\gtrsim 3 \sigma$ tension. In the standard cosmological model, $\Lambda$CDM, $t_U$ and $H_0$ are related by the matter density parameter $\Omega_{m,0}$. We propose to combine local $t_U$ constraints with late-time $\Omega_{m,0}$ estimates in a $\Lambda$CDM framework, to obtain a low-redshift $H_0$ determination that does not rely on early Universe physics. A proof-of-principle of this approach with current data gives $H_0=71\pm2.8$ ($H_0= 69.3 \pm 2.7$) km s$^{-1}$ Mpc$^{-1}$ from globular clusters (very-low-metallicity stars) with excellent prospects for improved constraints in the near future.

We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$\alpha$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$\alpha$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and Gaussian process interpolation. We show that this emulator has a typical accuracy of 1.5% and a worst-case accuracy of 4%, which compares well to the current statistical error of 3 - 5% at $z < 3$ from BOSS DR9. We compare to the previous state of the art, quadratic polynomial interpolation. The Latin hypercube samples the entire volume of parameter space, while quadratic polynomial emulation samples only lower-dimensional subspaces. The Gaussian process provides an estimate of the emulation error and we show using test simulations that this estimate is reasonable. We construct a likelihood function and use it to show that the posterior constraints generated using the emulator are unbiased. We show that our Gaussian process emulator has lower emulation error than quadratic polynomial interpolation and thus produces tighter posterior confidence intervals, which will be essential for future Lyman-$\alpha$ surveys such as DESI.

The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood requires a robust and accurate mapping between the parameters and observables, in this case the 1D flux power spectrum. Cosmological simulations enable such a mapping, but due to computational time constraints can only be evaluated at a handful of sample points; "emulators" are designed to interpolate between these. The problem then reduces to placing the sample points such that an accurate mapping is obtained while minimising the number of expensive simulations required. To address this, we introduce an emulation procedure that employs Bayesian optimisation of the training set for a Gaussian process interpolation scheme. Starting with a Latin hypercube sampling (other schemes with good space-filling properties can be used), we iteratively augment the training set with extra simulations at new parameter positions which balance the need to reduce interpolation error while focussing on regions of high likelihood. We show that smaller emulator error from the Bayesian optimisation propagates to smaller widths on the posterior distribution. Even with fewer simulations than a Latin hypercube, Bayesian optimisation shrinks the 95% credible volume by 90% and, e.g., the 1 sigma error on the amplitude of small-scale primordial fluctuations by 38%. This is the first demonstration of Bayesian optimisation applied to large-scale structure emulation, and we anticipate the technique will generalise to many other probes such as galaxy clustering, weak lensing and 21cm.

The expansion history of the Universe can be constrained in a cosmology-independent way by measuring the differential age evolution of cosmic chronometers. This yields a measurement of the Hubble parameter $H(z)$ as a function of redshift. The most reliable cosmic chronometers known so far are extremely massive and passively evolving galaxies. Age-dating these galaxies is, however, a difficult task, and even a small contribution of an underlying young stellar population could, in principle, affect the age estimate and its cosmological interpretation. We present several spectral indicators to detect, quantify and constrain such contamination in old galaxies, and study how their combination can be used to maximize the purity of cosmic chronometers selection. In particular, we analyze the CaII H/K ratio, the presence (or absence) of H$\alpha$ and [OII] emission lines, higher order Balmer absorption lines, and UV flux; each indicator is especially sensitive to a particular age range, allowing us to detect young components ranging between 10 Myr and 1 Gyr. The combination of these indicators minimizes the contamination to a level below 1% in the case of ideal data. More importantly, it offers a way to control the systematic error on $H(z)$ as a function of the contamination by young stellar populations. We show that for our previous measurements of the Hubble parameter, the possible bias induced by the presence of a younger component is well below the current errors. We envision that these indicators will be instrumental in paving the road for a robust and reliable dating of the old population and its cosmological interpretation.

The determination of the inflationary energy scale represents one of the first step towards the understanding of the early Universe physics. The (very mild) non-Gaussian signals that arise from any inflation model carry information about the energy scale of inflation and may leave an imprint in some cosmological observables, for instance on the clustering of high-redshift, rare and massive collapsed structures. In particular, the graviton exchange contribution due to interactions between scalar and tensor fluctuations leaves a specific signature in the four-point function of curvature perturbations, thus on clustering properties of collapsed structures. We compute the contribution of graviton exchange on two- and three-point function of halos, showing that at large scales $k\sim 10^{-3}\ \mathrm{Mpc}^{-1}$ its magnitude is comparable or larger to that of other primordial non-Gaussian signals discussed in the literature. This provides a potential route to probe the existence of tensor fluctuations which is alternative and highly complementary to B-mode polarisation measurements of the cosmic microwave background radiation.

We compute in detail how deviations from Einstein gravity at the inflation energy scale could appear as non-Gaussian features in the sky. To illustrate this we use multi-field $\alpha-$attractor models in the framework of supergravity to realise inflation. We find no obvious obstacle for having choices of model's parameters that generate non-Gaussian features of the equilateral and local type at the $\cal O$(1) level in the $f_{\rm NL}$ non-Gaussianity parameter, thus being potentially detectable in future cosmological surveys. This non-Gaussianity has its origin in either the non-canonical kinetic term (which, in turn, is an immediate consequence of assuming an hyperbolic geometry of the moduli space), the interactions of the fields in the potential or the $\alpha-$parameter, or a combination of these three. This opens up the exciting possibility of constraining the law of gravity at energy scales close to the Planck one.

We develop a practical methodology to remove modes from a galaxy survey power spectrum that are associated with systematic errors. We apply this to the BOSS CMASS sample, to see if it removes the excess power previously observed beyond the best-fit $\Lambda$CDM model on very large scales. We consider several possible sources of data contamination, and check whether they affect the number of targets that can be observed and the power spectrum measurements. We describe a general framework for how such knowledge can be transformed into template fields. Mode subtraction can then be used to remove these systematic contaminants at least as well as applying corrective weighting to the observed galaxies, but benefits from giving an unbiased power. Even after applying templates for all known systematics, we find a large-scale power excess, but this is reduced compared with that observed using the weights provided by the BOSS team. This excess is at much larger scales than the BAO scale and does not affect the main results of BOSS. However, it will be important for the measurement of a scale-dependent bias due to primordial non-Gaussianity. The excess is beyond that allowed by any simple model of non-Gaussianity matching Planck data, and is not matched in other surveys. We show that this power excess can further be reduced but is still present using "phenomenological" templates, designed to consider further potentially unknown sources of systematic contamination. As all discrepant angular modes can be removed using "phenomenological" templates, the potentially remaining contaminant acts radially.

The number of observed dwarf galaxies, with dark matter mass $\lesssim 10^{11}$ M$_{\odot}$ in the Milky Way or the Andromeda galaxy does not agree with predictions from the successful $\Lambda$CDM paradigm. To alleviate this problem a suppression of dark matter clustering power on very small scales has been conjectured. However, the abundance of dark matter halos outside our immediate neighbourhood (the Local Group) seem to agree with the $\Lambda$CDM--expected abundance. Here we connect these problems to observations of weak lensing cosmic shear, pointing out that cosmic shear can make significant statements about the missing satellites problem in a statistical way. As an example and pedagogical application we use recent constraints on small-scales power suppression from measurements of the CFHTLenS data. We find that, on average, in a region of $\sim $Gpc$^3$ there is no significant small-scale power suppression. This implies that suppression of small-scale power is not a viable solution to the `missing satellites problem' or, alternatively, that on average in this volume there is no `missing satellites problem' for dark matter masses $\gtrsim 5 \times 10^9$ M$_{\odot}$. Further analysis of current and future weak lensing surveys will probe much smaller scales, $k > 10h$ Mpc$^{-1}$ corresponding roughly to masses $M < 10^9 M_{\odot}$.

It is broadly accepted that Supermassive Black Holes (SMBHs) are located in the centers of most massive galaxies, although there is still no convincing scenario for the origin of their massive seeds. It has been suggested that primordial black holes (PBHs) of masses $\gtrsim 10^{2} M_\odot$ may provide such seeds, which would grow to become SMBHs. We suggest an observational test to constrain this hypothesis: gas accretion around PBHs during the cosmic dark ages powers the emission of high energy photons which would modify the spin temperature as measured by 21cm Intensity Mapping (IM) observations. We model and compute their contribution to the standard sky-averaged signal and power spectrum of 21cm IM, accounting for its substructure and angular dependence for the first time. If PBHs exist, the sky-averaged 21cm IM signal in absorption would be higher, while we expect an increase in the power spectrum for $\ell~\gtrsim 10^2-10^3$. We also forecast PBH detectability and measurement errors in the abundance and Eddington ratios for different fiducial parameter configurations for various future experiments, ranging from SKA to a futuristic radio array on the dark side of the Moon. While the SKA could provide a detection, only a more ambitious experiment would provide accurate measurements.

We forecast constraints on primordial non-Gaussianity (PNG) and bias parameters from measurements of galaxy power spectrum and bispectrum in future radio continuum and optical surveys. In the galaxy bispectrum, we consider a comprehensive list of effects, including the bias expansion for non-Gaussian initial conditions up to second order, redshift space distortions, redshift uncertainties and theoretical errors. These effects are all combined in a single PNG forecast for the first time. Moreover, we improve the bispectrum modelling over previous forecasts, by accounting for trispectrum contributions. All effects have an impact on final predicted bounds, which varies with the type of survey. We find that the bispectrum can lead to improvements up to a factor $\sim 5$ over bounds based on the power spectrum alone, leading to significantly better constraints for local-type PNG, with respect to current limits from \textit{Planck}. Future radio and photometric surveys could obtain a measurement error of $\sigma(f_{\mathrm{NL}}^{\mathrm{loc}}) \approx 0.2$. In the case of equilateral PNG, galaxy bispectrum can improve upon present bounds only if significant improvements in the redshift determinations of future, large volume, photometric or radio surveys could be achieved. For orthogonal non-Gaussianity, expected constraints are generally comparable to current ones.

An interesting test on the nature of the Universe is to measure the global spatial curvature of the metric in a model independent way, at a level of $|\Omega_k|<10^{-4}$, or, if possible, at the cosmic variance level of the amplitude of the CMB fluctuations $|\Omega_k|\approx10^{-5}$. A limit of $|\Omega_k|<10^{-4}$ would yield stringent tests on several models of inflation. Further, improving the constraint by an order of magnitude would help in reducing "model confusion" in standard parameter estimation. Moreover, if the curvature is measured to be at the value of the amplitude of the CMB fluctuations, it would offer a powerful test on the inflationary paradigm and would indicate that our Universe must be significantly larger than the current horizon. On the contrary, in the context of standard inflation, measuring a value above CMB fluctuations will lead us to conclude that the Universe is not much larger than the current observed horizon, this can also be interpreted as the presence of large fluctuations outside the horizon. However, it has proven difficult, so far, to find observables that can achieve such level of accuracy, and, most of all, be model-independent. Here we propose a method that can in principle achieve that, this is done by making minimal assumptions and using distance probes that are cosmology-independent: gravitational waves, redshift drift and cosmic chronometers. We discuss what kind of observations are needed in principle to achieve the desired accuracy.

Given the important role that the galaxy bispectrum has recently acquired in cosmology and the scale and precision of forthcoming galaxy clustering observations, it is timely to derive the full expression of the large-scale bispectrum going beyond approximated treatments which neglect integrated terms or higher-order bias terms or use the Limber approximation. On cosmological scales, relativistic effects that arise from observing on the past light-cone alter the observed galaxy number counts, therefore leaving their imprints on N-point correlators at all orders. In this paper we compute for the first time the bispectrum including all general relativistic, local and integrated, effects at second order, the tracers' bias at second order, geometric effects as well as the primordial non-Gaussianity contribution. This is timely considering that future surveys will probe scales comparable to the horizon where approximations widely used currently may not hold; neglecting these effects may introduce biases in estimation of cosmological parameters as well as primordial non-Gaussianity.

The model in which Primordial Black Holes (PBHs) constitute a non-negligible fraction of the dark matter has (re)gained popularity after the first detections of binary black hole mergers. Most of the observational constraints to date have been derived assuming a single mass for all the PBHs, although some more recent works tried to generalize constraints to the case of extended mass functions. Here we derive a general methodology to obtain constraints for any PBH Extended Mass Distribution (EMD) and any observables in the desired mass range. Starting from those obtained for a monochromatic distribution, we convert them into constraints for EMDs by using an equivalent, effective mass $M_{\rm eq}$ that depends on the specific observable. We highlight how limits of validity of the PBH modelling affect the EMD parameter space. Finally, we present converted constraints on the total abundance of PBH from microlensing, stellar distribution in ultra-faint dwarf galaxies and CMB accretion for Lognormal and Power Law mass distributions, finding that EMD constraints are generally stronger than monochromatic ones.

The possibility that a relevant fraction of the dark matter might be comprised of Primordial Black Holes (PBHs) has been seriously reconsidered after LIGO's detection of a $\sim 30 M_{\odot}$ binary black holes merger. Despite the strong interest in the model, there is a lack of studies on possible cosmological implications and effects on cosmological parameters inference. We investigate correlations with the other standard cosmological parameters using cosmic microwave background observations, finding significant degeneracies, especially with the tilt of the primordial power spectrum and the sound horizon at radiation drag. However, these degeneracies can be greatly reduced with the inclusion of small scale polarization data. We also explore if PBHs as dark matter in simple extensions of the standard $\Lambda$CDM cosmological model induces extra degeneracies, especially between the additional parameters and the PBH's ones. Finally, we present cosmic microwave background constraints on the fraction of dark matter in PBHs, not only for monochromatic PBH mass distributions but also for popular extended mass distributions. Our results show that extended mass distribution's constraints are tighter, but also that a considerable amount of constraining power comes from the high-$\ell$ polarization data. Moreover, we constrain the shape of such mass distributions in terms of the correspondent constraints on the PBH mass fraction.

It has been recently suggested that small mass black holes (BHs) may become unstable due to quantum-gravitational effects and eventually decay, producing radiation, on a timescale shorter than the Hawking evaporation time. We argue that the existence of a population of low-mass Primordial Black Holes (PBHs) acting as a fraction of the Universe dark matter component can be used to test proposed models of quantum decay of BHs via their effect on galaxy number counts. We study what constraints future galaxy clustering measurements can set on quantum-gravity parameters governing the BH lifetime and PBH abundance. In case of no detection of such effects, this would rule out either the existence of a non-negligible number of small PBHs, or the BH quantum decay scenario (or both). In case of independent observations of PBHs, the observables discussed here could be used to study the quantum effects that modify the final fate of BHs.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

The configuration of the three neutrino masses can take two forms, known as the normal and inverted hierarchies. We compute the Bayesian evidence associated with these two hierarchies. Previous studies found a mild preference for the normal hierarchy, and this was driven by the asymmetric manner in which cosmological data has confined the available parameter space. Here we identify the presence of a second asymmetry, which is imposed by data from neutrino oscillations. By combining constraints on the squared-mass splittings with the limit on the sum of neutrino masses of $\Sigma m_\nu < 0.13$ eV, and using a minimally informative prior on the masses, we infer odds of 42:1 in favour of the normal hierarchy, which is classified as "strong" in the Jeffreys' scale. We explore how these odds may evolve in light of higher precision cosmological data, and discuss the implications of this finding with regards to the nature of neutrinos. Finally the individual masses are inferred to be $m_1 = 3.80^{+26.2}_{-3.73} \, \text{meV}, m_2 = 8.8^{+18}_{-1.2} \, \text{meV}, m_3 = 50.4^{+5.8}_{-1.2} \, \text{meV}$ ($95\%$ credible intervals).

In this paper we discuss the commonly-used limiting cases, or approximations, for two-point cosmic shear statistics. We discuss the most prominent assumptions in this statistic: the flat-sky (small angle limit), the Limber (Bessel-to-delta function limit) and the Hankel transform (large l-mode limit) approximations; that the vast majority of cosmic shear results to date have used simultaneously. We find that the combined effect of these approximations can suppress power by >1% on scales of l<40. A fully non-approximated cosmic shear study should use a spherical-sky, non-Limber-approximated power spectrum analysis; and a transform involving Wigner small-d matrices in place of the Hankel transform. These effects, unaccounted for, would constitute at least 11% of the total budget for systematic effects for a power spectrum analysis of a Euclid-like experiment; but they are unnecessary.

The relation between the halo field and the matter fluctuations (halo bias), in the presence of massive neutrinos depends on the total neutrino mass, massive neutrinos introduce an additional scale-dependence of the bias which is usually neglected in cosmological analyses. We investigate the magnitude of the systematic effect on interesting cosmological parameters induced by neglecting this scale dependence, finding that while it is not a problem for current surveys, it is non-negligible for future, denser or deeper ones depending on the neutrino mass, the maximum scale used for the analyses and the details of the nuisance parameters considered. However there is a simple recipe to account for the bulk of the effect as to make it fully negligible, which we illustrate and advocate should be included in analysis of forthcoming large-scale structure surveys.

Obtaining high-sensitivity measurements of degree-scale cosmic microwave background (CMB) polarization is the most direct path to detecting primordial gravitational waves. Robustly recovering any primordial signal from the dominant foreground emission will require high-fidelity observations at multiple frequencies, with excellent control of systematics. We explore the potential for a new platform for CMB observations, the Airlander 10 hybrid air vehicle, to perform this task. We show that the Airlander 10 platform, operating at commercial airline altitudes, is well-suited to mapping frequencies above 220 GHz, which are critical for cleaning CMB maps of dust emission. Optimizing the distribution of detectors across frequencies, we forecast the ability of Airlander 10 to clean foregrounds of varying complexity as a function of altitude, demonstrating its complementarity with both existing (Planck) and ongoing (C-BASS) foreground observations. This novel platform could play a key role in defining our ultimate view of the polarized microwave sky.

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

Primordial features are one of the most important extensions of the Standard Model of cosmology, providing a wealth of information on the primordial universe, ranging from discrimination between inflation and alternative scenarios, new particle detection, to fine structures in the inflationary potential. We study the prospects of future large-scale structure (LSS) surveys on the detection and constraints of these features. We classify primordial feature models into several classes, and for each class we present a simple template of power spectrum that encodes the essential physics. We study how well the most ambitious LSS surveys proposed to date, including both spectroscopic and photometric surveys, will be able to improve the constraints with respect to the current Planck data. We find that these LSS surveys will significantly improve the experimental sensitivity on features signals that are oscillatory in scales, due to the 3D information. For a broad range of models, these surveys will be able to reduce the errors of the amplitudes of the features by a factor of 5 or more, including several interesting candidates identified in the recent Planck data. Therefore, LSS surveys offer an impressive opportunity for primordial feature discovery in the next decade or two. We also compare the advantages of both types of surveys.

We investigate our knowledge of early universe cosmology by exploring how much additional energy density can be placed in different components beyond those in the $\Lambda$CDM model. To do this we use a method to separate early- and late-universe information enclosed in observational data, thus markedly reducing the model-dependency of the conclusions. We find that the 95\% credibility regions for extra energy components of the early universe at recombination are: non-accelerating additional fluid density parameter $\Omega_{\rm MR} < 0.006$ and extra radiation parameterised as extra effective neutrino species $2.3 < N_{\rm eff} < 3.2$ when imposing flatness. Our constraints thus show that even when analyzing the data in this largely model-independent way, the possibility of hiding extra energy components beyond $\Lambda$CDM in the early universe is seriously constrained by current observations. We also find that the standard ruler, the sound horizon at radiation drag, can be well determined in a way that does not depend on late-time Universe assumptions, but depends strongly on early-time physics and in particular on additional components that behave like radiation. We find that the standard ruler length determined in this way is $r_{\rm s} = 147.4 \pm 0.7$ Mpc if the radiation and neutrino components are standard, but the uncertainty increases by an order of magnitude when non-standard dark radiation components are allowed, to $r_{\rm s} = 150 \pm 5$ Mpc.

We measure and analyse the bispectrum of the final, Data Release 12, galaxy sample provided by the Baryon Oscillation Spectroscopic Survey, splitting by selection algorithm into LOWZ and CMASS galaxies. The LOWZ sample contains 361\,762 galaxies with an effective redshift of $z_{\rm LOWZ}=0.32$, and the CMASS sample 777\,202 galaxies with an effective redshift of $z_{\rm CMASS}=0.57$. Combining the power spectrum, measured relative to the line-of-sight, with the spherically averaged bispectrum, we are able to constrain the product of the growth of structure parameter, $f$, and the amplitude of dark matter density fluctuations, $\sigma_8$, along with the geometric Alcock-Paczynski parameters, the product of the Hubble constant and the comoving sound horizon at the baryon drag epoch, $H(z)r_s(z_d)$, and the angular distance parameter divided by the sound horizon, $D_A(z)/r_s(z_d)$. After combining pre-reconstruction RSD analyses of the power spectrum monopole, quadrupole and bispectrum monopole; with post-reconstruction analysis of the BAO power spectrum monopole and quadrupole, we find $f(z_{\rm LOWZ})\sigma_8(z_{\rm LOWZ})=0.427\pm 0.056$, $D_A(z_{\rm LOWZ})/r_s(z_d)=6.60 \pm 0.13$, $H(z_{\rm LOWZ})r_s(z_d)=(11.55\pm 0.38)10^3\,{\rm kms}^{-1}$ for the LOWZ sample, and $f(z_{\rm CMASS})\sigma_8(z_{\rm CMASS})=0.426\pm 0.029$, $D_A(z_{\rm CMASS})/r_s(z_d)=9.39 \pm 0.10$, $H(z_{\rm CMASS})r_s(z_d)=(14.02\pm 0.22)10^3\,{\rm kms}^{-1}$ for the CMASS sample. We find general agreement with previous BOSS DR11 and DR12 measurements. Combining our dataset with {\it Planck15} we perform a null test of General Relativity (GR) through the $\gamma$-parametrisation finding $\gamma=0.733^{+0.068}_{-0.069}$, which is $\sim2.7\sigma$ away from the GR predictions.

We perform a comprehensive cosmological study of the $H_0$ tension between the direct local measurement and the model-dependent value inferred from the Cosmic Microwave Background. With the recent measurement of $H_0$ this tension has raised to more than $3\sigma$. We consider changes in the early time physics without modifying the late time cosmology. We also reconstruct the late time expansion history in a model independent way with minimal assumptions using distances measures from Baryon Acoustic Oscillations and Type Ia Supernovae, finding that at $z<0.6$ the recovered shape of the expansion history is less than 5 % different than that of a standard LCDM model. These probes also provide a model insensitive constraint on the low-redshift standard ruler, measuring directly the combination $r_s h$ where $H_0=h \times 100$ km/s/Mpc and $r_s$ is the sound horizon at radiation drag (the standard ruler), traditionally constrained by CMB observations. Thus $r_s$ and $H_0$ provide absolute scales for distance measurements (anchors) at opposite ends of the observable Universe. We calibrate the cosmic distance ladder and obtain a model-independent determination of the standard ruler for acoustic scale, $r_s$. The tension in $H_0$ reflects a mismatch between our determination of $r_s$ and its standard, CMB-inferred value. Without including high-l Planck CMB polarization data (i.e., only considering the "recommended baseline" low-l polarisation and temperature and the high l temperature data), a modification of the early-time physics to include a component of dark radiation with an effective number of species around 0.4 would reconcile the CMB-inferred constraints, and the local $H_0$ and standard ruler determinations. The inclusion of the "preliminary" high-l Planck CMB polarisation data disfavours this solution.

Assuming the existence of standard rulers, standard candles and standard clocks, requiring only the cosmological principle, a metric theory of gravity, a smooth expansion history, and using state-of-the-art observations, we determine the length of the "low-redshift standard ruler". The data we use are a compilation of recent Baryon acoustic oscillation data (relying on the standard ruler), Type 1A supernov\ae\ (as standard candles), ages of early type galaxies (as standard clocks) and local determinations of the Hubble constant (as a local anchor of the cosmic distance scale). In a standard $\Lambda$CDM cosmology the "low-redshift standard ruler" coincides with the sound horizon at radiation drag, which can also be determined --in a model dependent way-- from CMB observations. However, in general, the two quantities need not coincide. We obtain constraints on the length of the low-redshift standard ruler: $r^h_{\rm s}=101.0 \pm 2.3 h^{-1}$ Mpc, when using only Type 1A supernov\ae\ and Baryon acoustic oscillations, and $r_{\rm s}=150.0\pm 4.7 $ Mpc when using clocks to set the Hubble normalisation, while $r_{\rm s}=141.0\pm 5.5 $ Mpc when using the local Hubble constant determination (using both yields $r_{\rm s}=143.9\pm 3.1 $ Mpc). The low-redshift determination of the standard ruler has an error which is competitive with the model-dependent determination from cosmic microwave background measurements made with the {\em Planck} satellite, which assumes it is the sound horizon at the end of baryon drag.

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum.

We demonstrate that a scalar field is unable to reverse its direction of motion while continuously exchanging energy with another fluid. If the rate of transfer is modulated by the scalar's acceleration, the field can undergo a rapid process of freezing, despite being displaced from the local minimum of its potential. This enables dark energy to form from any potential, regardless of its shape. The field's equation of state mimicks that of a cosmological constant. We present a physically motivated realisation in the form of a derivative neutrino-majoron coupling. Coherent motions, which form only once the neutrinos become non-relativistic, could be responsible for instigating the freezing process. This would provide a natural resolution to the dark energy coincidence problem, while avoiding the dynamical instabilities associated with mass-varying neutrino models. Finally we discuss possible means by which this model could be experimentally verified.

Recent anomalies found in cosmological datasets such as the low multipoles of the Cosmic Microwave Background or the low redshift amplitude and growth of clustering measured by e.g., abundance of galaxy clusters and redshift space distortions in galaxy surveys, have motivated explorations of models beyond standard $\Lambda$CDM. Of particular interest are models where general relativity (GR) is modified on large cosmological scales. Here we consider deviations from $\Lambda$CDM+GR within the context of Horndeski gravity, which is the most general theory of gravity with second derivatives in the equations of motion. We adopt a parametrization in which the four additional Horndeski functions of time $\alpha_i(t)$ are proportional to the cosmological density of dark energy $\Omega_{DE}(t)$. Constraints on this extended parameter space using a suite of state-of-the art cosmological observations are presented for the first time. Although the theory is able to accommodate the low multipoles of the Cosmic Microwave Background and the low amplitude of fluctuations from redshift space distortions, we find no significant tension with $\Lambda$CDM+GR when performing a global fit to recent cosmological data and thus there is no evidence against $\Lambda$CDM+GR from an analysis of the value of the Bayesian evidence ratio of the modified gravity models with respect to $\Lambda$CDM, despite introducing extra parameters. The posterior distribution of these extra parameters that we derive return strong constraints on any possible deviations from $\Lambda$CDM+GR in the context of Horndeski gravity. We illustrate how our results can be applied to a more general frameworks of modified gravity models.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We present a minimally parametric, model independent reconstruction of the shape of the primordial power spectrum. Our smoothing spline technique is well-suited to search for smooth features such as deviations from scale invariance, and deviations from a power law such as running of the spectral index or small-scale power suppression. We use a comprehensive set of the state-of the art cosmological data: {\it Planck} observations of the temperature and polarisation anisotropies of the cosmic microwave background, WiggleZ and Sloan Digital Sky Survey Data Release 7 galaxy power spectra and the Canada-France-Hawaii Lensing Survey correlation function. This reconstruction strongly supports the evidence for a power law primordial power spectrum with a red tilt and disfavours deviations from a power law power spectrum including small-scale power suppression such as that induced by significantly massive neutrinos. This offers a powerful confirmation of the inflationary paradigm, justifying the adoption of the inflationary prior in cosmological analyses.

Deriving the expansion history of the Universe is a major goal of modern cosmology. To date, the most accurate measurements have been obtained with Type Ia Supernovae and Baryon Acoustic Oscillations, providing evidence for the existence of a transition epoch at which the expansion rate changes from decelerated to accelerated. However, these results have been obtained within the framework of specific cosmological models that must be implicitly or explicitly assumed in the measurement. It is therefore crucial to obtain measurements of the accelerated expansion of the Universe independently of assumptions on cosmological models. Here we exploit the unprecedented statistics provided by the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 to provide new constraints on the Hubble parameter $H(z)$ using the em cosmic chronometers approach. We extract a sample of more than 130000 of the most massive and passively evolving galaxies, obtaining five new cosmology-independent $H(z)$ measurements in the redshift range $0.3<z<0.5$, with an accuracy of $\sim$11-16\% incorporating both statistical and systematic errors. Once combined, these measurements yield a 6\% accuracy constraint of $H(z=0.4293)=91.8\pm5.3$ km/s/Mpc. The new data are crucial to provide the first cosmology-independent determination of the transition redshift at high statistical significance, measuring $z_{t}=0.4\pm0.1$, and to significantly disfavor the null hypothesis of no transition between decelerated and accelerated expansion at 99.9\% confidence level. This analysis highlights the wide potential of the cosmic chronometers approach: it permits to derive constraints on the expansion history of the Universe with results competitive with standard probes, and most importantly, being the estimates independent of the cosmological model, it can constrain cosmologies beyond -and including- the $\Lambda$CDM model.

We present cosmological upper limits on the sum of active neutrino masses using large-scale power spectrum data from the WiggleZ Dark Energy Survey and from the Sloan Digital Sky Survey - Data Release 7 (SDSS-DR7) sample of Luminous Red Galaxies (LRG). Combining measurements on the Cosmic Microwave Background temperature and polarisation anisotropies by the Planck satellite together with WiggleZ power spectrum results in a neutrino mass bound of 0.37 eV at 95% C.L., while replacing WiggleZ by the SDSS-DR7 LRG power spectrum, the 95% C.L. bound on the sum of neutrino masses is 0.38 eV. Adding Baryon Acoustic Oscillation (BAO) distance scale measurements, the neutrino mass upper limits greatly improve, since BAO data break degeneracies in parameter space. Within a $\Lambda$CDM model, we find an upper limit of 0.13 eV (0.14 eV) at 95% C.L., when using SDSS-DR7 LRG (WiggleZ) together with BAO and Planck. The addition of BAO data makes the neutrino mass upper limit robust, showing only a weak dependence on the power spectrum used. We also quantify the dependence of neutrino mass limit reported here on the CMB lensing information. The tighter upper limit (0.13 eV) obtained with SDSS-DR7 LRG is very close to that recently obtained using Lyman-alpha clustering data, yet uses a completely different probe and redshift range, further supporting the robustness of the constraint. This constraint puts under some pressure the inverted mass hierarchy and favours the normal hierarchy.

In General Relativity, the constraint equation relating metric and density perturbations is inherently nonlinear, leading to an effective non-Gaussianity in the dark matter density field on large scales - even if the primordial metric perturbation is Gaussian. Intrinsic non-Gaussianity in the large-scale dark matter overdensity in GR is real and physical. However, the variance smoothed on a local physical scale is not correlated with the large-scale curvature perturbation, so that there is no relativistic signature in the galaxy bias when using the simplest model of bias. It is an open question whether the observable mass proxies such as luminosity or weak lensing correspond directly to the physical mass in the simple halo bias model. If not, there may be observables that encode this relativistic signature.

We use the latest compilation of observational H(z) measurements obtained with cosmic chronometers in the redshift range $0<z<2$ to place constraints on cosmological parameters. We consider the sample alone and in combination with other state-of-the art cosmological probes: CMB data from the latest Planck 2015 release, the most recent estimate of the Hubble constant $H_{0}$, a compilation of recent BAO data, and the latest SNe sample. Since cosmic chronometers are independent of the assumed cosmological model, we are able to provide constraints on the parameters that govern the expansion history of the Universe in a way that can be used to test cosmological models. We show that the H(z) measurements obtained with cosmic chronometer from the BOSS survey provide enough constraining power in combination with CMB data to constrain the time evolution of dark energy, yielding constraints competitive with those obtained using SNe and/or BAO. From late-Universe probes alone we find that $w_0=-0.9\pm0.18$ and $w_a=-0.5\pm1.7$, and when combining also CMB data we obtain $w_0=-0.98\pm0.11$and $w_a=-0.30\pm0.4$. These new constraints imply that nearly all quintessence models are disfavoured, only phantom models or a pure cosmological constant being allowed. For the curvature we find $\Omega_k=0.003\pm0.003$, including CMB data. Cosmic chronometers data are important also to constrain neutrino properties by breaking or reducing degeneracies with other parameters. We find that $N_{eff}=3.17\pm0.15$, thus excluding the possibility of an extra (sterile) neutrino at more than $5\sigma$, and put competitive limits on the sum of neutrino masses, $\Sigma m_{\nu}< 0.27$ eV at 95% confidence level. Finally, we constrain the redshift evolution of dark energy, and find w(z) consistent with the $\Lambda$CDM model at the 40% level over the entire redshift range $0<z<2$. [abridged]

There is currently a discrepancy in the measured value of the amplitude of matter clustering, parameterised using $\sigma_8$, inferred from galaxy weak lensing, and cosmic microwave background data, which could be an indication of new physics, such as massive neutrinos or a modification to the gravity law, or baryon feedback. In this paper we make the assumption that the cosmological parameters are well determined by Planck, and use weak lensing data to investigate the implications for baryon feedback and massive neutrinos, as well as possible contributions from intrinsic alignments and biases in photometric redshifts. We apply a non-parametric approach to model the baryonic feedback on the dark matter clustering, which is flexible enough to reproduce the OWLS and Illustris simulation results. The statistic we use, 3D cosmic shear, is a method that extracts cosmological information from weak lensing data using a spherical-Bessel function power spectrum approach. We analyse the CFHTLenS weak lensing data and, assuming best fit cosmological parameters from the Planck CMB experiment, find that there is no evidence for baryonic feedback on the dark matter power spectrum, but there is evidence for a bias in the photometric redshifts in the CFHTLenS data, consistent with a completely independent analysis by Choi et al. (2015), based on spectroscopic redshifts; and that these conclusions are robust to assumptions about the intrinsic alignment systematic. We also find an upper limit on the sum of neutrino masses conditional on other $\Lambda$CDM parameters being fixed, of $< 0.28$ eV ($1\sigma$).

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

Recent constrains on the sum of neutrino masses inferred by analyzing cosmological data, show that detecting a non-zero neutrino mass is within reach of forthcoming cosmological surveys, implying a direct determination of the absolute neutrino mass scale. The measurement relies on constraining the shape of the matter power spectrum below the neutrino free streaming scale: massive neutrinos erase power at these scales. Detection of a lack of small-scale power, however, could also be due to a host of other effects. It is therefore of paramount importance to validate neutrinos as the source of power suppression at small scales. We show that, independent on hierarchy, neutrinos always show a footprint on large, linear scales; the exact location and properties can be related to the measured power suppression (an astrophysical measurement) and atmospheric neutrinos mass splitting (a neutrino oscillation experiment measurement). This feature can not be easily mimicked by systematic uncertainties or modifications in the cosmological model. The measurement of such a feature, up to 1% relative change in the power spectrum, is a smoking gun for confirming the determination of the absolute neutrino mass scale from cosmological observations. It also demonstrates the synergy of astrophysics and particle physics experiments.

We perform an empirical consistency test of General Relativity/dark energy by disentangling expansion history and growth of structure constraints. We replace each late-universe parameter that describes the behavior of dark energy with two meta-parameters: one describing geometrical information in cosmological probes, and the other controlling the growth of structure. If the underlying model (a standard wCDM cosmology with General Relativity) is correct, that is under the null hypothesis, the two meta-parameters coincide. If they do not, it could indicate a failure of the model or systematics in the data. We present a global analysis using state-of-the-art cosmological data sets which points in the direction that cosmic structures prefer a weaker growth than that inferred by background probes. This result could signify inconsistencies of the model, the necessity of extensions to it or the presence of systematic errors in the data. We examine all these possibilities. The fact that the result is mostly driven by a specific sub-set of galaxy clusters abundance data, points to the need of a better understanding of this probe.

Galactic Gamma-Ray Bursts (GRBs) are copious sources of gamma-rays that can pose a threat to complex life. Using recent determinations of their rate and the probability of GRBs causing massive extinction, we explore what type of universes are most likely to harbour advanced forms of life. For this purpose we use cosmological N-body simulations to determine at what time and for what value of the cosmological constant ($\Lambda$) the chances of life being unaffected by cosmic explosions are maximised. We find that $\Lambda-$dominated universes favour the survival of life against GRBs. Within a $\Lambda$CDM cosmology, the parameters that govern the likelihood of life survival to GRBs are dictated by the value of $\Lambda$ and the age of the Universe. We find that we seem to live in a favorable point in this parameter phase space which minimises the exposure to cosmic explosions, yet maximises the number of main sequence (hydrogen-burning) stars around which advanced life forms can exist.

We present distance scale measurements from the baryon acoustic oscillation signal in the CMASS and LOWZ samples from the Data Release 12 of the Baryon Oscillation Spectroscopic Survey (BOSS). The total volume probed is 14.5 Gpc$^3$, a 10 per cent increment from Data Release 11. From an analysis of the spherically averaged correlation function, we infer a distance to $z=0.57$ of $D_V(z)r^{\rm fid}_{\rm d}/r_ {\rm d}=2028\pm21$ Mpc and a distance to $z=0.32$ of $D_V(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1264\pm22$ Mpc assuming a cosmology in which $r^{\rm fid}_{\rm d}=147.10$ Mpc. From the anisotropic analysis, we find an angular diameter distance to $z=0.57$ of $D_{\rm A}(z)r^{\rm fid}_{\rm d}/r_{\rm d}=1401\pm21$ Mpc and a distance to $z=0.32$ of $981\pm20$ Mpc, a 1.5 per cent and 2.0 per cent measurement respectively. The Hubble parameter at $z=0.57$ is $H(z)r_{\rm d}/r^{\rm fid}_{\rm d}=100.3\pm3.7$ km s$^{-1}$ Mpc$^{-1}$ and its value at $z=0.32$ is $79.2\pm5.6$ km s$^{-1}$ Mpc$^{-1}$, a 3.7 per cent and 7.1 per cent measurement respectively. These cosmic distance scale constraints are in excellent agreement with a $\Lambda$CDM model with cosmological parameters released by the recent Planck 2015 results.

Applying a transformation to a non-Gaussian field can enhance the information content of the resulting power spectrum, by reducing the correlations between Fourier modes. In the context of weak gravitational lensing, it has been shown that this gain in information content is significantly compromised by the presence of shape noise. We apply clipping to mock convergence fields, a technique which is known to be robust in the presence of noise and has been successfully applied to galaxy number density fields. When analysed in isolation the resulting convergence power spectrum returns degraded constraints on cosmological parameters. However substantial gains can be achieved by performing a combined analysis of the power spectra derived from both the original and transformed fields. Even in the presence of realistic levels of shape noise, we demonstrate that this approach is capable of reducing the area of likelihood contours within the $\Omega_m - \sigma_8$ plane by more than a factor of three.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We examine the cosmological implications of the measurements of the linear growth rate of cosmological structure obtained in a companion paper from the power spectrum and bispectrum monopoles of the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey Data, Release 11, CMASS galaxies. This measurement was of $f^{0.43}\sigma_8$, where $\sigma_8$ is the amplitude of dark matter density fluctuations, and $f$ is the linear growth rate, at the effective redshift of the survey, $z_{\rm eff}=0.57$. In conjunction with Cosmic Microwave Background (CMB) data, interesting constraints can be placed on models with non-standard neutrino properties and models where gravity deviates from general relativity on cosmological scales. In particular, the sum of the masses of the three species of the neutrinos is constrained to $m_\nu<0.49\,{\rm eV}$ (at 95\% confidence level) when the $f^{0.43}\sigma_8$ measurement is combined with state-of-the-art CMB measurements. Allowing the effective number of neutrinos to vary as a free parameter does not significantly change these results. When we combine the measurement of $f^{0.43}\sigma_8$ with the complementary measurement of $f\sigma_8$ from the monopole and quadrupole of the two-point correlation function we are able to obtain an independent measurements of $f$ and $\sigma_8$. We obtain $f=0.63\pm0.16$ and $\sigma_8=0.710\pm0.086$ (68\% confidence level). This is the first time when these parameters have been able to be measured independently using the redshift-space power spectrum and bispectrum measurements from galaxy clustering data only.

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

We present a detailed study of second-order matter perturbations for the general Horn- deski class of models. Being the most general scalar-tensor theory having second-order equations of motion, it includes many known gravity and dark energy theories and General Relativity with a cosmological constant as a specific case. This enables us to estimate the leading order dark matter bispectrum generated at late-times by gravitational instability. We parametrize the evolution of the first and second-order equations of motion as proposed by Bellini and Sawicki (2014), where the free functions of the theory are assumed to be proportional to the dark energy density. We show that it is unnatural to have large 10% ( 1%) deviations of the bispectrum introducing even larger ~ 30% (~ 5%) deviations in the linear growth rate. Considering that measurements of the linear growth rate have much higher signal-to-noise than bispectrum measurements, this indicates that for Horndeski models which reproduce the expansion history and the linear growth rate as predicted by GR the dark matter bispectrum kernel can be effectively modelled as the standard GR one. On the other hand, an observation of a large bispectrum deviation that can not be explained in terms of bias would imply either that the evolution of perturbations is strongly different than the evolution predicted by GR or that the theory of gravity is exotic (e.g., breaks the weak equivalence principle) and/or fine-tuned.

We analyse the anisotropic clustering of the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS Data Release 11 sample, which consists of $690 827$ galaxies in the redshift range $0.43 < z < 0.70$ and has a sky coverage of $8 498$ deg$^2$ corresponding to an effective volume of $\sim6\,\rm{Gpc}^3$. We fit the Fourier space statistics, the power spectrum and bispectrum monopoles to measure the linear and quadratic bias parameters, $b_1$ and $b_2$, for a non-linear non-local bias model, the growth of structure parameter $f$ and the amplitude of dark matter density fluctuations parametrised by $\sigma_8$. We obtain $b_1(z_{\rm eff})^{1.40}\sigma_8(z_{\rm eff})=1.672\pm 0.060$ and $b_2^{0.30}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.579\pm0.082$ at the effective redshift of the survey, $z_{\rm eff}=0.57$. The main cosmological result is the constraint on the combination $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.582\pm0.084$, which is complementary to $f\sigma_8$ constraints obtained from 2-point redshift space distortion analyses. A less conservative analysis yields $f^{0.43}(z_{\rm eff})\sigma_8(z_{\rm eff})=0.584\pm0.051$. We ensure that our result is robust by performing detailed systematic tests using a large suite of survey galaxy mock catalogs and N-body simulations. The constraints on $f^{0.43}\sigma_8$ are useful for setting additional constrains on neutrino mass, gravity, curvature as well as the number of neutrino species from galaxy surveys analyses (as presented in a companion paper).

This is the summary of two lectures that aim to give an overview of cosmology. I will not try to be too rigorous in derivations, nor to give a full historical overview. The idea is to provide a "taste" of cosmology and some of the interesting topics it covers. The standard cosmological model is presented and I highlight the successes of cosmology over the past decade or so. Keys to the development of the standard cosmological model are observations of the cosmic microwave background and of large-scale structure, which are introduced. Inflation and dark energy and the outlook for the future are also discussed. Slides from the lectures are available from the school website: physicschool.web.cern.ch/PhysicSchool/CLASHEP/CLASHEP2011/.

The existence of a cosmic neutrino background can be probed indirectly by CMB experiments, not only by measuring the background density of radiation in the universe, but also by searching for the typical signatures of the fluctuations of free-streaming species in the temperature and polarisation power spectrum. Previous studies have already proposed a rather generic parametrisation of these fluctuations, that could help to discriminate between the signature of ordinary free-streaming neutrinos, or of more exotic dark radiation models. Current data are compatible with standard values of these parameters, which seems to bring further evidence for the existence of a cosmic neutrino background. In this work, we investigate the robustness of this conclusion under various assumptions. We generalise the definition of an effective sound speed and viscosity speed to the case of massive neutrinos or other dark radiation components experiencing a non-relativistic transition. We show that current bounds on these effective parameters do not vary significantly when considering an arbitrary value of the particle mass, or extended cosmological models with a free effective neutrino number, dynamical dark energy or a running of the primordial spectrum tilt. We conclude that it is possible to make a robust statement about the detection of the cosmic neutrino background by CMB experiments.

We exploit cosmological-model independent measurements of the expansion history of the Universe to provide a cosmic distance ladder. These are supernovae type Ia used as standard candles (at redshift between 0.01 and 1.3) and baryon acoustic oscillations (at redshifts between 0.1 and 0.8) as standard rulers. We calibrate (anchor) the ladder in two ways: first using the local $H_0$ value as an anchor at $z$ = 0 (effectively calibrating the standard candles) and secondly using the cosmic microwave background-inferred sound-horizon scale as an anchor (giving the standard ruler length) as an inverse distance ladder. Both methods are consistent, but the uncertainty in the expansion history $H(z)$ is smaller if the sound horizon scale is used. We present inferred values for the sound horizon at radiation drag $r_d$ which do not rely on assumptions about the early expansion history nor on cosmic microwave background measurements but on the cosmic distance ladder and baryon acoustic oscillations measurements. We also present derived values of $H_0$ from the inverse distance ladder and we show that they are in very good agreement with the extrapolated value in a $\Lambda$CDM model from Planck cosmic microwave background data.

We present a phenomenological modification of the standard perturbation theory prediction for the bispectrum in redshift space that allows us to extend the model to mildly non-linear scales over a wide range of redshifts, $z\leq1.5$. We find that we can describe the bispectrum of dark matter particles with $\sim5%$ accuracy for $k_i\lesssim0.10\,h/{\rm Mpc}$ at $z=0$, for $k_i\lesssim0.15\,h/{\rm Mpc}$ at $z=0.5$, for $k_i\lesssim0.17\,h/{\rm Mpc}$ at $z=1.0$ and for $k_i\lesssim0.20\,h/{\rm Mpc}$ at $z=1.5$. We also test that the fitting formula is able to describe with similar accuracy the bispectrum of cosmologies with different $\Omega_m$, in the range $0.2\lesssim \Omega_m \lesssim 0.4$, and consequently with different values of the logarithmic grow rate $f$ at $z=0$, $0.4\lesssim f(z=0) \lesssim 0.6$. We apply this new formula to recover the bias parameters, $f$ and $\sigma_8$, by combining the redshift space power spectrum monopole and quadrupole with the bispectrum monopole for both dark matter particles and haloes. We find that the combination of these three statistics can break the degeneracy between $b_1$, $f$ and $\sigma_8$. For dark matter particles the new model can be used to recover $f$ and $\sigma_8$ with $\sim1%$ accuracy. For dark matter haloes we find that $f$ and $\sigma_8$ present larger systematic shifts, $\sim10%$. The systematic offsets arise because of limitations in the modelling of the interplay between bias and redshift space distortions, and represent a limitation as the statistical errors of forthcoming surveys reach this level. Conveniently, we find that these residual systematics are mitigated for combinations of parameters. The improvement on the modelling of the bispectrum presented in this paper will be useful for extracting information from current and future galaxy surveys. [abridged]

We measure the length of the Baryon Acoustic Oscillation (BAO) feature, and the expansion rate of the recent Universe, from low-redshift data only, almost model-independently. We make only the following minimal assumptions: homogeneity and isotropy; a metric theory of gravity; a smooth expansion history, and the existence of standard candles (supernov\ae) and a standard BAO ruler. The rest is determined by the data, which are compilations of recent BAO and Type IA supernova results. Making only these assumptions, we find for the first time that the standard ruler has length $103.9 \pm 2.3\, h^{-1}$ Mpc. The value is a measurement, in contrast to the model-dependent theoretical prediction determined with model parameters set by Planck data ($99.3 \pm 2.1 \, h^{-1}$ Mpc). The latter assumes $\Lambda$CDM, and that the ruler is the sound horizon at radiation drag. Adding passive galaxies as standard clocks or a local Hubble constant measurement allows the absolute BAO scale to be determined ($142.8\pm 3.7$ Mpc), and in the former case the additional information makes the BAO length determination more precise ($101.9\pm 1.9 \, h^{-1}\,$Mpc). The inverse curvature radius of the Universe is weakly constrained and consistent with zero, independently of the gravity model, provided it is metric. We find the effective number of relativistic species to be $N_{\rm eff} = 3.53\pm 0.32$, independent of late-time dark energy or gravity physics.

General relativistic cosmology cannot be reduced to linear relativistic perturbations superposed on an isotropic and homogeneous (Friedmann-Robertson-Walker) background, even though such a simple scheme has been successfully applied to analyse a large variety of phenomena (such as Cosmic Microwave Background primary anisotropies, matter clustering on large scales, weak gravitational lensing, etc.). The general idea of going beyond this simple paradigm is what characterises most of the efforts made in recent years: the study of second and higher-order cosmological perturbations including all general relativistic contributions -- also in connection with primordial non-Gaussianities -- the idea of defining large-scale structure observables directly from a general relativistic perspective, the various attempts to go beyond the Newtonian approximation in the study of non-linear gravitational dynamics, by using e.g., Post-Newtonian treatments, are all examples of this general trend. Here we summarise some of these directions of investigation, with the aim of emphasising future prospects in this area of cosmology, both from a theoretical and observational point of view.

Recent work has shown that the z~2.5 Lyman-alpha forest on large scales encodes information about the galaxy and quasar populations that keep the intergalactic medium photoionized. We present the first forecasts for constraining the populations with data from current and next-generation surveys. At a minimum the forest should tell us whether galaxies or, conversely, quasars dominate the photon production. The number density and clustering strength of the ionising sources might be estimated to sub-10% precision with a DESI-like survey if degeneracies (e.g., with the photon mean-free-path, small-scale clustering power normalization and potentially other astrophysical effects) can be broken by prior information. We demonstrate that, when inhomogeneous ionisation is correctly handled, constraints on dark energy do not degrade.

We show that modifications of Einstein gravity during inflation could leave potentially measurable imprints on cosmological observables in the form of non-Gaussian perturbations. This is due to the fact that these modifications appear in the form of an extra field that could have non-trivial interactions with the inflaton. We show it explicitly for the case $R+\alpha R^2$, where nearly scale-invariant non-Gaussianity at the level of $f_{\rm NL} \approx -(1$ to $30)$ can be obtained, in a quasi-local configuration.

It has been claimed recently that massive sterile neutrinos could bring about a new concordance between observations of the cosmic microwave background (CMB), the large-scale structure (LSS) of the Universe, and local measurements of the Hubble constant, $H_0$. We demonstrate that this apparent concordance results from combining datasets which are in significant tension, even within this extended model, possibly indicating remaining systematic biases in the measurements. We further show that this tension remains when the cosmological model is further extended to include significant tensor modes, as suggested by the recent BICEP2 results. Using the Bayesian evidence, we show that the minimal $\Lambda$CDM model is strongly favoured over its neutrino extensions by various combinations of datasets. Robust data combinations yield stringent limits of $\sum m_\nu\lesssim0.3$ eV and $m_{\nu,{\rm sterile}}^{\rm eff} \lesssim 0.3$ eV at $95\%$ CL for the sum of active and sterile neutrinos, respectively.

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.

We give an analytical form for the weighted correlation function of peaks in a Gaussian random field. In a cosmological context, this approach strictly describes the formation bias and is the main result here. Nevertheless, we show its validity and applicability to the evolved cosmological density field and halo field, using Gaussian random field realisations and dark matter N-body numerical simulations. Using this result from peak theory we compute the bias of peaks (and dark matter halos) and show that it reproduces results from the simulations at the ${\mathcal O}(10\%)$ level. Our analytical formula for the bias predicts a scale-dependent bias with two characteristics: a broad band shape which, however, is most affected by the choice of weighting scheme and evolution bias, and a more robust, narrow feature localised at the BAO scale, an effect that is confirmed in simulations. This scale-dependent bias smooths the BAO feature but, conveniently, does not move it. We provide a simple analytic formula to describe this effect. We envision that our analytic solution will be of use for galaxy surveys that exploit galaxy clustering.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

We use cosmology-independent measurements of the expansion history in the redshift range 0.1 < z <1.2 and compare them with the Cosmic Microwave Background-derived expansion history predictions. The motivation is to investigate if the tension between the local (cosmology independent) Hubble constant H0 value and the Planck-derived H0 is also present at other redshifts. We conclude that there is no tension between Planck and cosmology independent-measurements of the Hubble parameter H(z) at 0.1 < z < 1.2 for the LCDM model (odds of tension are only 1:15, statistically not significant). Considering extensions of the LCDM model does not improve these odds (actually makes them worse), thus favouring the simpler model over its extensions. On the other hand the H(z) data are also not in tension with the local H0 measurements but the combination of all three data-sets shows a highly significant tension (odds ~ 1:400). Thus the new data deepen the mystery of the mismatch between Planck and local H0 measurements, and cannot univocally determine wether it is an effect localised at a particular redshift. Having said this, we find that assuming the NGC4258 maser distance as the correct anchor for H0, brings the odds to comfortable values. Further, using only the expansion history measurements we constrain, within the LCDM model, H0 = 68.5 +- 3.5 and Omega_m = 0.32 +- 0.05 without relying on any CMB prior. We also address the question of how smooth the expansion history of the universe is given the cosmology independent data and conclude that there is no evidence for deviations from smoothness on the expansion history, neither variations with time in the value of the equation of state of dark energy.

We use Bayesian model comparison to determine whether extensions to Standard-Model neutrino physics -- primarily additional effective numbers of neutrinos and/or massive neutrinos -- are merited by the latest cosmological data. Given the significant advances in cosmic microwave background (CMB) observations represented by the Planck data, we examine whether Planck temperature and CMB lensing data, in combination with lower redshift data, have strengthened (or weakened) the previous findings. We conclude that the state-of-the-art cosmological data do not show evidence for deviations from the standard cosmological model (which has three massless neutrino families). This does not mean that the model is necessarily correct -- in fact we know it is incomplete as neutrinos are not massless -- but it does imply that deviations from the standard model (e.g., non-zero neutrino mass) are too small compared to the current experimental uncertainties to be inferred from cosmological data alone.

It has been shown by Shchekinov & Vasiliev2006 (SV06) that HD molecules can be an important cooling agent in high redshift z >10 haloes if they undergo mergers under specific conditions so suitable shocks are created. Here we build upon Prieto et al. (2012) who studied in detail the merger-generated shocks, and show that the conditions for HD cooling can be studied by combining these results with a suite of dark-matter only simulations. We have performed a number of dark matter only simulations from cosmological initial conditions inside boxes with sizes from 1 to 4 Mpc. We look for haloes with at least two progenitors of which at least one has mass M > M_cr (z), where M_cr (z) is the SV06 critical mass for HD over-cooling. We find that the fraction of over-cooled haloes with mass between M_cr (z) and 10^{0.2} M_cr (z), roughly below the atomic cooling limit, can be as high as ~ 0.6 at z ~ 10 depending on the merger mass ratio. This fraction decreases at higher redshift reaching a value ~0.2 at z ~ 15. For higher masses, i.e. above 10^{0.2} M_cr (z) up to 10^{0.6} M_cr (z), above the atomic cooling limit, this fraction rises to values ~ 0.8 until z ~ 12.5. As a consequence, a non negligible fraction of high redshift z > 10 mini-haloes can drop their gas temperature to the Cosmic Microwave Background temperature limit allowing the formation of low mass stars in primordial environments.

We use the latest Planck constraints, and in particular constraints on the derived parameters (Hubble constant and age of the Universe) for the local universe and compare them with local measurements of the same quantities. We propose a way to quantify whether cosmological parameters constraints from two different experiments are in tension or not. Our statistic, T, is an evidence ratio and therefore can be interpreted with the widely used Jeffrey's scale. We find that in the framework of the LCDM model, the Planck inferred two dimensional, joint, posterior distribution for the Hubble constant and age of the Universe is in "strong" tension with the local measurements; the odds being ~ 1:50. We explore several possibilities for explaining this tension and examine the consequences both in terms of unknown errors and deviations from the LCDM model. In some one-parameter LCDM model extensions, tension is reduced whereas in other extensions, tension is instead increased. In particular, small total neutrino masses are favored and a total neutrino mass above 0.15 eV makes the tension "highly significant" (odds ~ 1:150). A consequence of accepting this interpretation of the tension is that the degenerate neutrino hierarchy is highly disfavoured by cosmological data and the direct hierarchy is slightly favored over the inverse.

We provide an exact expression for the multi-variate joint probability distribution function of non-Gaussian fields primordially arising from local transformations of a Gaussian field. This kind of non-Gaussianity is generated in many models of inflation. We apply our expression to the non- Gaussianity estimation from Cosmic Microwave Background maps and the halo mass function where we obtain analytical expressions. We also provide analytic approximations and their range of validity. For the Cosmic Microwave Background we give a fast way to compute the PDF which is valid up to 7{\sigma} for fNL values (both true and sampled) not ruled out by current observations, which consists of expressing the PDF as a combination of bispectrum and trispectrum of the temperature maps. The resulting expression is valid for any kind of non-Gaussianity and is not limited to the local type. The above results may serve as the basis for a fully Bayesian analysis of the non-Gaussianity parameter.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

It has been suggested that recent cosmological and flavor-oscillation data favor the existence of additional neutrino species beyond the three predicted by the Standard Model of particle physics. We apply Bayesian model selection to determine whether there is indeed any evidence from current cosmological datasets for the standard cosmological model to be extended to include additional neutrino flavors. The datasets employed include cosmic microwave background temperature, polarization and lensing power spectra, and measurements of the baryon acoustic oscillation scale and the Hubble constant. We also consider other extensions to the standard neutrino model, such as massive neutrinos, and possible degeneracies with other cosmological parameters. The Bayesian evidence indicates that current cosmological data do not require any non-standard neutrino properties.

We explore how local, cosmology-independent measurements of the Hubble constant and the age of the Universe help to provide a powerful consistency check of the currently favored cosmological model (flat LambdaCDM) and model-independent constraints on cosmology. We use cosmic microwave background (CMB) data to define the model-dependent cosmological parameters, and add local measurements to assess consistency and determine whether extensions to the model are justified. At current precision, there is no significant tension between the locally measured Hubble constant and age of the Universe (with errors of 3% and 5% respectively) and the corresponding parameters derived from the CMB. However, if errors on the local measurements could be decreased by a factor of two, one could decisively conclude if there is tension or not. We also compare the local and CMB data assuming simple extensions of the flat, $\Lambda$CDM model (including curvature, dark energy with a constant equation of state parameter not equal to -1, non-zero neutrino masses and a non-standard number of neutrino families) and find no need for these extra parameters; in particular, we constrain the effective number of neutrino species to be Neff < 4 at 95% confidence. We show that local measurements provide constraints on the curvature and equation of state of dark energy nearly orthogonal to those of the CMB. We argue that cosmology-independent measurements of local quantities at the percent level would be very useful to explore cosmology in a model-independent way.

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.

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 {\AA} break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a \Lambda CDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2.

We have assembled a compilation of observational Hubble parameter measurements estimated with the differential evolution of cosmic chronometers, in the redshift range 0<z<1.75. This sample has been used, in combination with CMB data and with the most recent estimate of the Hubble constant H_0, to derive new constraints on several cosmological parameters. The new Hubble parameter data are very useful to break some of the parameter degeneracies present in CMB-only analysis, and to constrain possible deviations from the standard (minimal) flat \Lambda CDM model. The H(z) data are especially valuable in constraining \Omega_k and \Omega_DE in models that allow a variation of those parameters, yielding constraints that are competitive with those obtained using Supernovae and/or baryon acoustic oscillations. We also find that our H(z) data are important to constrain parameters that do no affect directly the expansion history, by breaking or reducing degeneracies with other parameters. We find that Nrel=3.45\pm0.33 using WMAP 7-years data in combination with South Pole Telescope data and our H(z) determinations (Nrel=3.71\pm0.45 using Atacama Cosmology Telescope data instead of South Pole Telescope). We exclude Nrel>4 at 95% CL (74% CL) using the same datasets combinations. We also put competitive limits on the sum of neutrino masses, \Sigma m_\nu<0.24 eV at 68% confidence level. These results have been proven to be extremely robust to many possible systematic effects, such as the initial choice of stellar population synthesis model adopted to estimate H(z) and the progenitor-bias.

We investigate the accuracy of Eulerian perturbation theory for describing the matter and galaxy power spectra in real and redshift space in light of future observational probes for precision cosmology. Comparing the analytical results with a large suite of N-body simulations (160 independent boxes of 13.8 (Gpc/h)^3 volume each, which are publicly available), we find that re-summing terms in the standard perturbative approach predicts the real-space matter power spectrum with an accuracy of < 2% for k < 0.20 h/Mpc at redshifts z < 1.5. This is obtained following the widespread technique of writing the resummed propagator in terms of 1-loop contributions. We show that the accuracy of this scheme increases by considering higher-order terms in the resummed propagator. By combining resummed perturbation theories with several models for the mappings from real to redshift space discussed in the literature, the multipoles of the dark-matter power spectrum can be described with sub-percent deviations from N-body results for k < 0.15h/Mpc at z < 1. As a consequence, the logarithmic growth rate, f, can be recovered with sub-percent accuracy on these scales. Extending the models to massive dark-matter haloes in redshift space, our results describe the monopole term from N-body data within 2% accuracy for scales k < 0.15 h/Mpc at z < 0.5; here f can be recovered within < 5% when the halo bias is known. We conclude that these techniques are suitable to extract cosmological information from future galaxy surveys.

We work out an effective theory of accelerated expansion to describe general phenomena of inflation and acceleration (dark energy) in the Universe. Our aim is to determine from theoretical grounds, in a physically-motivated and model independent way, which and how many (free) parameters are needed to broadly capture the physics of a theory describing cosmic acceleration. Our goal is to make as much as possible transparent the physical interpretation of the parameters describing the expansion. We show that, at leading order, there are five independent parameters, of which one can be constrained via general relativity tests. The other four parameters need to be determined by observing and measuring the cosmic expansion rate only, H(z). Therefore we suggest that future cosmology surveys focus on obtaining an accurate as possible measurement of $H(z)$ to constrain the nature of accelerated expansion (dark energy and/or inflation).

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.

We consider whether the non-Gaussian scale-dependent halo bias can be used not only to constrain the local form of non-Gaussianity but also to distinguish among different shapes. In particular, we ask whether it can constrain the behavior of the primordial three-point function in the squeezed limit where one of the momenta is much smaller than the other two. This is potentially interesting since the observation of a three-point function with a squeezed limit that does not go like the local nor equilateral templates would be a signal of non-trivial dynamics during inflation. To this end we use the quasi-single field inflation model of Chen and Wang as a representative two-parameter model, where one parameter governs the amplitude of non-Gaussianity and the other the shape. We also perform a model-independent analysis by parametrizing the scale-dependent bias as a power-law on large scales, where the power is to be constrained from observations. We find that proposed large-scale structure surveys (with characteristics similar to the dark energy task force stage IV surveys) have the potential to distinguish among the squeezed limit behavior of different bispectrum shapes for a wide range of fiducial model parameters. Thus the halo bias can help discriminate between different models of inflation.

We present new observational constraints on inhomogeneous models based on observables independent of the CMB and large-scale structure. Using Bayesian evidence we find very strong evidence for the homogeneous LCDM model, thus disfavouring inhomogeneous models. Our new constraints are based on quantities independent of the growth of perturbations and rely on cosmic clocks based on atomic physics and on the local density of matter.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We compare the stellar population properties in the central regions of visually classified non-starforming spiral and elliptical galaxies from Galaxy Zoo and SDSS DR7. The galaxies lie in the redshift range $0.04<z<0.1$ and have stellar masses larger than $logM_*=10.4$. We select only face-on spiral galaxies in order to avoid contamination by light from the disk in the SDSS fiber and enabling the robust visual identification of spiral structure. Overall, we find that galaxies with larger central stellar velocity dispersions, regardless of morphological type, have older ages, higher metallicities, and an increased overabundance of alpha-elements. Age and alpha-enhancement, at fixed velocity dispersion, do not depend on morphological type. The only parameter that, at a given velocity dispersion, correlates with morphological type is metallicity, where the metallicity of the bulges of spiral galaxies is 0.07 dex higher than that of the ellipticals. However, for galaxies with a given total stellar mass, this dependence on morphology disappears. Under the assumption that, for our sample, the velocity dispersion traces the mass of the bulge alone, as opposed to the total mass (bulge+disk) of the galaxy, our results imply that the formation epoch of galaxy and the duration of its star-forming period are linked to the mass of the bulge. The extent to which metals are retained within the galaxy, and not removed as a result of outflows, is determined by the total mass of the galaxy.

We explore the benefits of using a passively evolving population of galaxies to measure the evolution of the rate of structure growth between z=0.25 and z=0.65 by combining data from the SDSS-I/II and SDSS-III surveys. The large-scale linear bias of a population of dynamically passive galaxies, which we select from both surveys, is easily modeled. Knowing the bias evolution breaks degeneracies inherent to other methodologies, and decreases the uncertainty in measurements of the rate of structure growth and the normalization of the galaxy power-spectrum by up to a factor of two. If we translate our measurements into a constraint on sigma_8(z=0) assuming a concordance cosmological model and General Relativity (GR), we find that using a bias model improves our uncertainty by a factor of nearly 1.5. Our results are consistent with a flat Lambda Cold Dark Matter model and with GR.

We implement Slow Roll Reconstruction -- an optimal solution to the inverse problem for inflationary cosmology -- within ModeCode, a publicly available solver for the inflationary dynamics. We obtain up-to-date constraints on the reconstructed inflationary potential, derived from the WMAP 7-year dataset and South Pole Telescope observations, combined with large scale structure data derived from SDSS Data Release 7. Using ModeCode in conjunction with the MultiNest sampler, we compute Bayesian evidence for the reconstructed potential at each order in the truncated slow roll hierarchy. We find that the data are well-described by the first two slow roll parameters, \epsilon and \eta, and that there is no need to include a nontrivial \xi parameter.

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

With future galaxy surveys, a huge number of Fourier modes of the distribution of the large scale structures in the Universe will become available. These modes are complementary to those of the CMB and can be used to set constraints on models of the early universe, such as inflation. Using a MCMC analysis, we compare the power of the CMB with that of the combination of CMB and galaxy survey data, to constrain the power spectrum of primordial fluctuations generated during inflation. We base our analysis on the Planck satellite and a spectroscopic redshift survey with configuration parameters close to those of the Euclid mission as examples. We first consider models of slow-roll inflation, and show that the inclusion of large scale structure data improves the constraints by nearly halving the error bars on the scalar spectral index and its running. If we attempt to reconstruct the inflationary single-field potential, a similar conclusion can be reached on the parameters characterizing the potential. We then study models with features in the power spectrum. In particular, we consider ringing features produced by a break in the potential and oscillations such as in axion monodromy. Adding large scale structures improves the constraints on features by more than a factor of two. In axion monodromy we show that there are oscillations with small amplitude and frequency in momentum space that are undetected by CMB alone but can be measured by including galaxy surveys in the analysis.

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.

We analyze the anisotropic clustering of massive galaxies from the Sloan Digital Sky Survey III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) sample, which consists of 264,283 galaxies in the redshift range 0.43 < z < 0.7 spanning 3,275 square degrees. Both peculiar velocities and errors in the assumed redshift-distance relation ("Alcock-Paczynski effect") generate correlations between clustering amplitude and orientation with respect to the line-of-sight. Together with the sharp baryon acoustic oscillation (BAO) standard ruler, our measurements of the broadband shape of the monopole and quadrupole correlation functions simultaneously constrain the comoving angular diameter distance (2190 +/- 61 Mpc) to z=0.57, the Hubble expansion rate at z=0.57 (92.4 +/- 4.5 km/s/Mpc), and the growth rate of structure at that same redshift (d sigma8/d ln a = 0.43 +/- 0.069). Our analysis provides the best current direct determination of both DA and H in galaxy clustering data using this technique. If we further assume a LCDM expansion history, our growth constraint tightens to d sigma8/d ln a = 0.415 +/- 0.034. In combination with the cosmic microwave background, our measurements of DA, H, and growth all separately require dark energy at z > 0.57, and when combined imply \Omega_{\Lambda} = 0.74 +/- 0.016, independent of the Universe's evolution at z<0.57. In our companion paper (Samushia et al. prep), we explore further cosmological implications of these observations.

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 have performed cosmological N-body simulations which include the effect of the masses of the individual neutrino species. The simulations were aimed at studying the effect of different neutrino hierarchies on the matter power spectrum. Compared to the linear theory predictions, we find that non-linearities enhance the effect of hierarchy on the matter power spectrum at mildly non-linear scales. The difference between the different hierarchies is about 0.5% for a sum of neutrino masses of 0.1eV. Albeit this is a small effect, it is potentially measurable from upcoming surveys. In combination with neutrinoless double-beta decay experiments, this opens up the possibility of using the sky to determine if neutrinos are Majorana or Dirac fermions.

In this paper we present an improved fitting formula for the dark matter bispectrum motivated by the previous phenomenological approach of Scoccimarro & Couchman (2001). We use a set of LCDM simulations to calibrate the fitting parameters in the k-range of 0.03 h/Mpc<k<0.4 h/Mpc and in the redshift range of 0<z<1.5. This new proposed fit describes well the BAO-features although it was not designed to. The deviation between the simulations output and our analytic prediction is typically less than 5% and in the worst case is never above 10%. We envision that this new analytic fitting formula will be very useful in providing reliable predictions for the non-linear dark matter bispectrum for LCDM models.

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

Conversion of photons into axions under the presence of a strong magnetic field can dim the radiation from magnetized astrophysical objects. Here we perform a detailed calculation aimed at quantifying the signatures of photon-axion conversion in the spectra, light curves, and polarization of neutron stars (NSs). We take into account the energy and angle-dependence of the conversion probability and the surface thermal emission from NSs. The latter is computed from magnetized atmosphere models that include the effect of photon polarization mode conversion due to vacuum polarization. The resulting spectral models, inclusive of the general-relativistic effects of gravitational redshift and light deflection, allow us to make realistic predictions for the effects of photon to axion conversion on observed NS spectra, light curves, and polarization signals. We identify unique signatures of the conversion, such as an increase of the effective area of a hot spot as it rotates away from the observer line of sight. For a star emitting from the entire surface, the conversion produces apparent radii that are either larger or smaller (depending on axion mass and coupling strength) than the limits set by NS equations of state. For an emission region that is observed phase-on, photon-axion conversion results in an inversion of the plane of polarization with respect to the no-conversion case. While the quantitative details of the features that we identify depend on NS properties (magnetic field strength, temperature) and axion parameters, the spectral and polarization signatures induced by photon-axion conversion are distinctive enough to make NSs very interesting and promising probes of axion physics.

Using observational data on the expansion rate of the universe (H(z)) we constrain the effective Lagrangian of the current accelerated expansion. Our results show that the effective potential is consistent with being flat i.e., a cosmological constant; it is also consistent with the field moving along an almost flat potential like a pseudo-Goldstone boson. We show that the potential of dark energy does not deviate from a constant at more than 6% over the redshift range 0 < z < 1. The data can be described by just a constant term in the Lagrangian and do not require any extra parameters; therefore there is no evidence for augmenting the number of parameters of the LCDM paradigm. We also find that the data justify the effective theory approach to describe accelerated expansion and that the allowed parameters range satisfy the expected hierarchy. Future data, both from cosmic chronometers and baryonic acoustic oscillations, that can measure H(z) at the % level, could greatly improve constraints on the flatness of the potential or shed some light on possible mechanisms driving the accelerated expansion. Besides the above result, it is shown that the effective Lagrangian of accelerated expansion can be constrained from cosmological observations in a model-independent way and that direct measurements of the expansion rate H(z) are most useful to do so.

We critically investigate current statistical tests applied to high redshift clusters of galaxies in order to test the standard cosmological model and describe their range of validity. We carefully compare a sample of high-redshift, massive, galaxy clusters with realistic Poisson sample simulations of the theoretical mass function, which include the effect of Eddington bias. We compare the observations and simulations using the following statistical tests: the distributions of ensemble and individual existence probabilities (in the >M,>z sense), the redshift distributions, and the 2d Kolmogorov-Smirnov test. Using seemingly rare clusters from Hoyle et al. (2011), and Jee et al. (2011) and assuming the same survey geometry as in Jee et al. (2011, which is less conservative than Hoyle et al. 2011), we find that the (>M,>z) existence probabilities of all clusters are fully consistent with LCDM. However assuming the same survey geometry, we use the 2d K-S test probability to show that the observed clusters are not consistent with being the least probable clusters from simulations at >95% confidence, and are also not consistent with being a random selection of clusters, which may be caused by the non-trivial selection function and survey geometry. Tension can be removed if we examine only a X-ray selected sub sample, with simulations performed assuming a modified survey geometry.

We measure the acoustic scale from the angular power spectra of the Sloan Digital Sky Survey III (SDSS-III) Data Release 8 imaging catalog that includes 872,921 galaxies over ~ 10,000 deg^2 between 0.45<z<0.65. The extensive spectroscopic training set of the Baryon Oscillation Spectroscopic Survey (BOSS) luminous galaxies allows precise estimates of the true redshift distributions of galaxies in our imaging catalog. Utilizing the redshift distribution information, we build templates and fit to the power spectra of the data, which are measured in our companion paper, Ho et al. 2011, to derive the location of Baryon acoustic oscillations (BAO) while marginalizing over many free parameters to exclude nearly all of the non-BAO signal. We derive the ratio of the angular diameter distance to the sound horizon scale D_A/r_s= 9.212 + 0.416 -0.404 at z=0.54, and therefore, D_A= 1411+- 65 Mpc at z=0.54; the result is fairly independent of assumptions on the underlying cosmology. Our measurement of angular diameter distance D_A is 1.4 \sigma higher than what is expected for the concordance LCDM (Komatsu et al. 2011), in accordance to the trend of other spectroscopic BAO measurements for z >~ 0.35. We report constraints on cosmological parameters from our measurement in combination with the WMAP7 data and the previous spectroscopic BAO measurements of SDSS (Percival et al. 2010) and WiggleZ (Blake et al. 2011). We refer to our companion papers (Ho et al. 2011; de Putter et al. 2011) for investigations on information of the full power spectrum.

We present N-body simulations for generic non-Gaussian initial conditions with the aim of exploring and modelling the scale-dependent halo bias. This effect is evident at very large scales requiring large simulation boxes. In addition, the previously available prescription to implement generic non-Gaussian initial conditions has been improved to keep under control higher-order terms which were spoiling the power spectrum on large scales. We pay particular attention to the differences between physical, inflation-motivated primordial bispectra and their factorizable templates, and to the operational definition of the non-Gaussian halo bias (which has both a scale-dependent and an approximately scale-independent contributions). We find that analytic predictions for both the non-Gaussian halo mass function and halo bias work well once a calibration factor (which was introduced before) is calibrated on simulations. The halo bias remains therefore an extremely promising tool to probe primordial non-Gaussianity and thus to give insights into the physical mechanism that generated the primordial perturbations. The simulation outputs and tables of the analytic predictions will be made publicly available via the non-Gaussian comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html

In this paper we analyze a suite of cosmological simulations of modified gravitational action f(R) models, where cosmic acceleration is induced by a scalar field that acts as a fifth force on all forms of matter. In particular, we focus on the bispectrum of the dark matter density field on mildly non-linear scales. For models with the same initial power spectrum, the dark matter bispectrum shows significant differences for cases where the final dark matter power spectrum also differs. Given the different dependence on bias of the galaxy power spectrum and bispectrum, bispectrum measurements can close the loophole of galaxy bias hiding differences in the power spectrum. Alternatively, changes in the initial power spectrum can also hide differences. By constructing LCDM models with very similar final non-linear power spectra, we show that the differences in the bispectrum are reduced (<4%) and are comparable with differences in the imperfectly matched power spectra. These results indicate that the bispectrum depends mainly on the power spectrum and less sensitively on the gravitational signatures of the f(R) model. This weak dependence of the matter bispectrum on gravity makes it useful for breaking degeneracies associated with galaxy bias, even for models beyond general relativity.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

We consider cosmological parameters estimation in the presence of a non-zero isocurvature contribution in the primordial perturbations. A previous analysis showed that even a tiny amount of isocurvature perturbation, if not accounted for, could affect standard rulers calibration from Cosmic Microwave Background observations such as those provided by the Planck mission, affect Baryon Acoustic Oscillations interpretation and introduce biases in the recovered dark energy properties that are larger than forecasted statistical errors from future surveys. Extending on this work, here we adopt a general fiducial cosmology which includes a varying dark energy equation of state parameter and curvature. Beside Baryon Acoustic Oscillations measurements, we include the information from the shape of the galaxy power spectrum and consider a joint analysis of a Planck-like Cosmic Microwave Background probe and a future, space-based, Large Scale Structure probe not too dissimilar from recently proposed surveys. We find that this allows one to break the degeneracies that affect the Cosmic Microwave Background and Baryon Acoustic Oscillations combination. As a result, most of the cosmological parameter systematic biases arising from an incorrect assumption on the isocurvature fraction parameter $f_{iso}$, become negligible with respect to the statistical errors. We find that the Cosmic Microwave Background and Large Scale Structure combination gives a statistical error $\sigma(f_{iso}) \sim 0.008$, even when curvature and a varying dark energy equation of state are included, which is smaller than the error obtained from Cosmic Microwave Background alone when flatness and cosmological constant are assumed. These results confirm the synergy and complementarity between Cosmic Microwave Background and Large Scale Structure, and the great potential of future and planned galaxy surveys.

In this paper we measure the angular power spectra $C_\ell$ of three high-redshift large-scale structure probes: the radio sources from the NRAO VLA Sky Survey (NVSS), the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs) and the MegaZ-LRG (DR7), the final SDSS II Luminous Red Galaxy (LRG) photometric redshift survey. We perform a global analysis of the constraints on the amplitude of primordial non-Gaussianity from these angular power spectra, as well as from their cross-correlation power spectra with the cosmic microwave background (CMB) temperature map. In particular, we include non-Gaussianity of the type arising from single-field slow roll, multifields, curvaton (local type), and those which effects on the halo clustering can be described by the equilateral template (related to higher-order derivative type non-Gaussianity) and by the enfolded template (related to modified initial state or higher-derivative interactions). When combining all data sets, we obtain limits of $f_{\rm NL}=48\pm20$, $f_{\rm NL}=50\pm265$ and $f_{\rm NL}=183\pm95$ at 68% confidence level for local, equilateral and enfolded templates, respectively. Furthermore, we explore the constraint on the cubic correction $g_{\rm NL}\phi^3$ on the bias of dark matter haloes and obtain a limit of $-1.2\times10^5<g_{\rm NL}<11.3\times10^5$ at 95% confidence level.

Building on the legacy of the Sloan Digital Sky Survey (SDSS-I and II), SDSS-III is a program of four spectroscopic surveys on three scientific themes: dark energy and cosmological parameters, the history and structure of the Milky Way, and the population of giant planets around other stars. In keeping with SDSS tradition, SDSS-III will provide regular public releases of all its data, beginning with SDSS DR8 (which occurred in Jan 2011). This paper presents an overview of the four SDSS-III surveys. BOSS will measure redshifts of 1.5 million massive galaxies and Lya forest spectra of 150,000 quasars, using the BAO feature of large scale structure to obtain percent-level determinations of the distance scale and Hubble expansion rate at z<0.7 and at z~2.5. SEGUE-2, which is now completed, measured medium-resolution (R=1800) optical spectra of 118,000 stars in a variety of target categories, probing chemical evolution, stellar kinematics and substructure, and the mass profile of the dark matter halo from the solar neighborhood to distances of 100 kpc. APOGEE will obtain high-resolution (R~30,000), high signal-to-noise (S/N>100 per resolution element), H-band (1.51-1.70 micron) spectra of 10^5 evolved, late-type stars, measuring separate abundances for ~15 elements per star and creating the first high-precision spectroscopic survey of all Galactic stellar populations (bulge, bar, disks, halo) with a uniform set of stellar tracers and spectral diagnostics. MARVELS will monitor radial velocities of more than 8000 FGK stars with the sensitivity and cadence (10-40 m/s, ~24 visits per star) needed to detect giant planets with periods up to two years, providing an unprecedented data set for understanding the formation and dynamical evolution of giant planet systems. (Abridged)

We present a method to minimize, or even cancel out, the nuisance parameters affecting a measurement. Our approach is general and can be applied to any experiment or observation. We compare it with the bayesian technique used to deal with nuisance parameters: marginalization, and show how the method compares and improves by avoiding biases. We illustrate the method with several examples taken from the astrophysics and cosmology world: baryonic acoustic oscillations, cosmic clocks, Supernova Type Ia luminosity distance, neutrino oscillations and dark matter detection. By applying the method we recover some known results but also find some interesting new ones. For baryonic acoustic oscillation (BAO) experiments we show how to combine radial and angular BAO measurements in order to completely eliminate the dependence on the sound horizon at radiation drag. In the case of exploiting SN1a as standard candles we show how the uncertainty in the luminosity distance by a second parameter modeled as a metallicity dependence can be eliminated or greatly reduced. When using cosmic clocks to measure the expansion rate of the universe, we demonstrate how a particular combination of observables nearly removes the metallicity dependence of the galaxy on determining differential ages, thus removing the age-metallicity degeneracy in stellar populations. We hope that these findings will be useful in future surveys to obtain robust constraints on the dark energy equation of state.

It has been suggested that both recent cosmological data and the results of flavor oscillation experiments (MiniBooNE and LSND) lend support to the existence of low-mass sterile neutrinos. The cosmological data appear to weakly favor additional forms of radiation in the Universe, beyond photons and three standard neutrino families. We reconsider the cosmological evidence by making the resulting confidence intervals on the additional effective neutrino species as prior-independent as possible. We find that, once the prior-dependence is removed, the latest cosmological data show no evidence for deviations from the standard number of neutrino species.

To date, 14 high-redshift (z>1.0) galaxy clusters with mass measurements have been observed, spectroscopically confirmed and are reported in the literature. These objects should be exceedingly rare in the standard LCDM model. We conservatively approximate the selection functions of these clusters' parent surveys, and quantify the tension between the abundances of massive clusters as predicted by the standard LCDM model and the observed ones. We alleviate the tension considering non-Gaussian primordial perturbations of the local type, characterized by the parameter fnl and derive constraints on fnl arising from the mere existence of these clusters. At the 95% confidence level, fnl>467 with cosmological parameters fixed to their most likely WMAP5 values, or fnl > 123 (at 95% confidence) if we marginalize over WMAP5 parameters priors. In combination with fnl constraints from Cosmic Microwave Background and halo bias, this determination implies a scale-dependence of fnl at approx. 3 sigma. Given the assumptions made in the analysis, we expect any future improvements to the modeling of the non-Gaussian mass function, survey volumes, or selection functions to increase the significance of fnl>0 found here. In order to reconcile these massive, high-z clusters with an fnl=0, their masses would need to be systematically lowered by 1.5 sigma or the sigma8 parameter should be approx. 3 sigma higher than CMB (and large-scale structure) constraints. The existence of these objects is a puzzle: it either represents a challenge to the LCDM paradigme or it is an indication that the mass estimates of clusters is dramatically more uncertain than we think.

COrE (Cosmic Origins Explorer) is a fourth-generation full-sky, microwave-band satellite recently proposed to ESA within Cosmic Vision 2015-2025. COrE will provide maps of the microwave sky in polarization and temperature in 15 frequency bands, ranging from 45 GHz to 795 GHz, with an angular resolution ranging from 23 arcmin (45 GHz) and 1.3 arcmin (795 GHz) and sensitivities roughly 10 to 30 times better than PLANCK (depending on the frequency channel). The COrE mission will lead to breakthrough science in a wide range of areas, ranging from primordial cosmology to galactic and extragalactic science. COrE is designed to detect the primordial gravitational waves generated during the epoch of cosmic inflation at more than $3\sigma $ for $r=(T/S)>=10^{-3}$. It will also measure the CMB gravitational lensing deflection power spectrum to the cosmic variance limit on all linear scales, allowing us to probe absolute neutrino masses better than laboratory experiments and down to plausible values suggested by the neutrino oscillation data. COrE will also search for primordial non-Gaussianity with significant improvements over Planck in its ability to constrain the shape (and amplitude) of non-Gaussianity. In the areas of galactic and extragalactic science, in its highest frequency channels COrE will provide maps of the galactic polarized dust emission allowing us to map the galactic magnetic field in areas of diffuse emission not otherwise accessible to probe the initial conditions for star formation. COrE will also map the galactic synchrotron emission thirty times better than PLANCK. This White Paper reviews the COrE science program, our simulations on foreground subtraction, and the proposed instrumental configuration.

Current results from the Lyman alpha forest assume that the primordial power spectrum of density perturbations follows a simple power law form, with running. We present the first analysis of Lyman alpha data to study the effect of relaxing this strong assumption on primordial and astrophysical constraints. We perform a large suite of numerical simulations, using them to calibrate a minimally parametric framework for describing the power spectrum. Combined with cross-validation, a statistical technique which prevents over-fitting of the data, this framework allows us to reconstruct the power spectrum shape without strong prior assumptions. We find no evidence for deviation from scale-invariance; our analysis also shows that current Lyman alpha data do not have sufficient statistical power to robustly probe the shape of the power spectrum at these scales. In contrast, the ongoing Baryon Oscillation Sky Survey (BOSS) will be able to do so with high precision. Furthermore, this near-future data will be able to break degeneracies between the power spectrum shape and astrophysical parameters.

We examine whether future, nearly all-sky galaxy redshift surveys, in combination with CMB priors, will be able to detect the signature of the cosmic neutrino background and determine the absolute neutrino mass scale. We also consider what constraints can be imposed on the effective number of neutrino species. In particular we consider two spectroscopic strategies in the near-IR, the so-called "slitless" and "multi-slit" approaches, whose examples are given by future space-based galaxy surveys, as EUCLID for the slitless case, or SPACE, JEDI, and possibly WFIRST in the future, for the multi-slit case. We find that, in combination with Planck, these galaxy probes will be able to detect at better than 3--sigma level and measure the mass of cosmic neutrinos: a) in a cosmology-independent way, if the sum of neutrino masses is above 0.1 eV; b) assuming spatial flatness and that dark energy is a cosmological constant, otherwise. We find that the sensitivity of such surveys is well suited to span the entire range of neutrino masses allowed by neutrino oscillation experiments, and to yield a clear detection of non-zero neutrino mass. The detection of the cosmic relic neutrino background with cosmological experiments will be a spectacular confirmation of our model for the early Universe and a window into one of the oldest relic components of our Universe.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

We explore how the Peccei-Quinn (PQ) axion parameter space can be constrained by the frequency-dependence dimming of radiation from astrophysical objects. To do so we perform accurate calculations of photon-axion conversion in the presence of a variable magnetic field. We propose several tests where the PQ axion parameter space can be explored with current and future astronomical surveys: the observed spectra of isolated neutron stars, occultations of background objects by white dwarfs and neutron stars, the light curves of eclipsing binaries containing a white dwarf. We find that the lack of dimming of the light curve of a detached eclipsing white dwarf binary recently observed, leads to relevant constraints on the photon-axion conversion. Current surveys designed for Earth-like planet searches are well matched to strengthen and improve the constraints on the PQ axion using astrophysical objects radiation dimming.

We present a modification of the standard halo model with the goal of providing an improved description of galaxy clustering. Recent surveys, like the Sloan Digital Sky Survey (SDSS) and the Anglo-Australian Two-degree survey (2dF), have shown that there seems to be a correlation between the clustering of galaxies and their properties such as metallicity and star formation rate, which are believed to be environment-dependent. This environmental dependence is not included in the standard halo model where the host halo mass is the only variable specifying galaxy properties. In our approach, the halo properties i.e., the concentration, and the Halo Occupation Distribution --HOD-- prescription, will not only depend on the halo mass (like in the standard halo model) but also on the halo environment. We examine how different environmental dependence of halo concentration and HOD prescription affect the correlation function. We see that at the level of dark matter, the concentration of haloes affects moderately the dark matter correlation function only at small scales. However the galaxy correlation function is extremely sensitive to the HOD details, even when only the HOD of a small fraction of haloes is modified.

Using the publicly available VESPA database of SDSS Data Release 7 spectra, we calculate the stellar Mass Weighted Age (hereafter MWA) as a function of local galaxy density and dark matter halo mass. We compare our results with semi-analytic models from the public Millennium Simulation. We find that the stellar MWA has a large scatter which is inherent in the data and consistent with that seen in semi-analytic models. The stellar MWA is consistent with being independent (to first order) with local galaxy density, which is also seen in semi-analytic models. As a function of increasing dark matter halo mass (using the SDSS New York Value Added Group catalogues), we find that the average stellar MWA for member galaxies increases, which is again found in semi-analytic models. Furthermore we use public dark matter Mass Accretion History (MAH) code calibrated on simulations, to calculate the dark matter Mass Weighted Age as a function of dark matter halo mass. In agreement with earlier analyses, we find that the stellar MWA and the dark matter MWA are anti correlated for large mass halos, i.e, dark matter accretion does not seem to be the primary factor in determining when stellar mass was compiled. This effect can be described by down-sizing.

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.

The strong dependence of the large-scale dark matter halo bias on the (local) non-Gaussianity parameter, f_NL, offers a promising avenue towards constraining primordial non-Gaussianity with large-scale structure surveys. In this paper, we present the first detection of the dependence of the non-Gaussian halo bias on halo formation history using N-body simulations. We also present an analytic derivation of the expected signal based on the extended Press-Schechter formalism. In excellent agreement with our analytic prediction, we find that the halo formation history-dependent contribution to the non-Gaussian halo bias (which we call non-Gaussian halo assembly bias) can be factorized in a form approximately independent of redshift and halo mass. The correction to the non-Gaussian halo bias due to the halo formation history can be as large as 100%, with a suppression of the signal for recently formed halos and enhancement for old halos. This could in principle be a problem for realistic galaxy surveys if observational selection effects were to pick galaxies occupying only recently formed halos. Current semi-analytic galaxy formation models, for example, imply an enhancement in the expected signal of ~23% and ~48% for galaxies at z=1 selected by stellar mass and star formation rate, respectively.

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.

Comparing clustering of differently biased tracers of the dark matter distribution offers the opportunity to reduce the cosmic variance error in the measurement of certain cosmological parameters. We develop a formalism that includes bias non-linearities and stochasticity. Our formalism is general enough that can be used to optimise survey design and tracers selection and optimally split (or combine) tracers to minimise the error on the cosmologically interesting quantities. Our approach generalises the one presented by McDonald & Seljak (2009) of circumventing sample variance in the measurement of $f\equiv d \ln D/d\ln a$. We analyse how the bias, the noise, the non-linearity and stochasticity affect the measurements of $Df$ and explore in which signal-to-noise regime it is significantly advantageous to split a galaxy sample in two differently-biased tracers. We use N-body simulations to find realistic values for the parameters describing the bias properties of dark matter haloes of different masses and their number density. We find that, even if dark matter haloes could be used as tracers and selected in an idealised way, for realistic haloes, the sample variance limit can be reduced only by up to a factor $\sigma_{2tr}/\sigma_{1tr}\simeq 0.6$. This would still correspond to the gain from a three times larger survey volume if the two tracers were not to be split. Before any practical application one should bear in mind that these findings apply to dark matter haloes as tracers, while realistic surveys would select galaxies: the galaxy-host halo relation is likely to introduce extra stochasticity, which may reduce the gain further.

The Atacama Cosmology Telescope (ACT) is currently observing the cosmic microwave background with arcminute resolution at 148 GHz, 218 GHz, and 277 GHz. In this paper, we present ACT's first results. Data have been analyzed using a maximum-likelihood map-making method which uses B-splines to model and remove the atmospheric signal. It has been used to make high-precision beam maps from which we determine the experiment's window functions. This beam information directly impacts all subsequent analyses of the data. We also used the method to map a sample of galaxy clusters via the Sunyaev-Zel'dovich (SZ) effect, and show five clusters previously detected with X-ray or SZ observations. We provide integrated Compton-y measurements for each cluster. Of particular interest is our detection of the z = 0.44 component of A3128 and our current non-detection of the low-redshift part, providing strong evidence that the further cluster is more massive as suggested by X-ray measurements. This is a compelling example of the redshift-independent mass selection of the SZ effect.

We present a catalog of 105 rich and massive ($M>3\times10^{14}M_{\sun}$) optically-selected clusters of galaxies extracted from 70 square-degrees of public archival griz imaging from the Blanco 4-m telescope acquired over 45 nights between 2005 and 2007. We use the clusters' optically-derived properties to estimate photometric redshifts, optical luminosities, richness, and masses. We complement the optical measurements with archival XMM-Newton and ROSAT X-ray data which provide additional luminosity and mass constraints on a modest fraction of the cluster sample. Two of our clusters show clear evidence for central lensing arcs; one of these has a spectacular large-diameter, nearly-complete Einstein Ring surrounding the brightest cluster galaxy. A strong motivation for this study is to identify the massive clusters that are expected to display prominent signals from the Sunyaev-Zeldovich Effect (SZE) and therefore be detected in the wide-area mm-band surveys being conducted by both the Atacama Cosmology Telescope and the South Pole Telescope. The optical sample presented here will be useful for verifying new SZE cluster candidates from these surveys, for testing the cluster selection function, and for stacking analyzes of the SZE data.

We present a measurement of the angular power spectrum of the cosmic microwave background (CMB) radiation observed at 148 GHz. The measurement uses maps with 1.4' angular resolution made with data from the Atacama Cosmology Telescope (ACT). The observations cover 228 square degrees of the southern sky, in a 4.2-degree-wide strip centered on declination 53 degrees South. The CMB at arcminute angular scales is particularly sensitive to the Silk damping scale, to the Sunyaev-Zel'dovich (SZ) effect from galaxy clusters, and to emission by radio sources and dusty galaxies. After masking the 108 brightest point sources in our maps, we estimate the power spectrum between 600 < \ell < 8000 using the adaptive multi-taper method to minimize spectral leakage and maximize use of the full data set. Our absolute calibration is based on observations of Uranus. To verify the calibration and test the fidelity of our map at large angular scales, we cross-correlate the ACT map to the WMAP map and recover the WMAP power spectrum from 250 < ell < 1150. The power beyond the Silk damping tail of the CMB is consistent with models of the emission from point sources. We quantify the contribution of SZ clusters to the power spectrum by fitting to a model normalized at sigma8 = 0.8. We constrain the model's amplitude ASZ < 1.63 (95% CL). If interpreted as a measurement of sigma8, this implies sigma8^SZ < 0.86 (95% CL) given our SZ model. A fit of ACT and WMAP five-year data jointly to a 6-parameter LCDM model plus terms for point sources and the SZ effect is consistent with these results.

We present an analysis of the constraints on the amplitude of primordial non-Gaussianity of local type described by the dimensionless parameter $f_{\rm NL}$. These constraints are set by the auto-correlation functions (ACFs) of two large scale structure probes, the radio sources from NRAO VLA Sky Survey (NVSS) and the quasar catalogue of Sloan Digital Sky Survey Release Six (SDSS DR6 QSOs), as well as by their cross-correlation functions (CCFs) with the cosmic microwave background (CMB) temperature map (Integrated Sachs-Wolfe effect). Several systematic effects that may affect the observational estimates of the ACFs and of the CCFs are investigated and conservatively accounted for. Our approach exploits the large-scale scale-dependence of the non-Gaussian halo bias. The derived constraints on {$f_{\rm NL}$} coming from the NVSS CCF and from the QSO ACF and CCF are weaker than those previously obtained from the NVSS ACF, but still consistent with them. Finally, we obtain the constraints on $f_{\rm NL}=53\pm25$ ($1\,\sigma$) and $f_{\rm NL}=58\pm24$ ($1\,\sigma$) from NVSS data and SDSS DR6 QSO data, respectively.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

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.

We address the issue of setting up generic non-Gaussian initial conditions for N-body simulations. We consider inflationary-motivated primordial non-Gaussianity where the perturbations in the Bardeen potential are given by a dominant Gaussian part plus a non-Gaussian part specified by its bispectrum. The approach we explore here is suitable for any bispectrum, i.e. it does not have to be of the so-called separable or factorizable form. The procedure of generating a non-Gaussian field with a given bispectrum (and a given power spectrum for the Gaussian component) is not univocal, and care must be taken so that higher-order corrections do not leave a too large signature on the power spectrum. This is so far a limiting factor of our approach. We then run N-body simulations for the most popular inflationary-motivated non-Gaussian shapes. The halo mass function and the non-linear power spectrum agree with theoretical analytical approximations proposed in the literature, even if they were so far developed and tested only for a particular shape (the local one). We plan to make the simulations outputs available to the community via the non-Gaussian simulations comparison project web site http://icc.ub.edu/~liciaverde/NGSCP.html.

The measurement of Baryonic Acoustic Oscillations from galaxy surveys is well known to be a robust and powerful tool to constrain dark energy. This method relies on the knowledge of the size of the acoustic horizon at radiation drag derived from Cosmic Microwave Background Anisotropy measurements. In this paper we quantify the effect of non-standard initial conditions in the form of an isocurvature component on the determination of dark energy parameters from future BAO surveys. In particular, if there is an isocurvature component (at a level still allowed by present data) but it is ignored in the CMB analysis, the sound horizon and cosmological parameters determination is biased, and, as a consequence, future surveys may incorrectly suggest deviations from a cosmological constant. In order to recover an unbiased determination of the sound horizon and dark energy parameters, a component of isocurvature perturbations must be included in the model when analyzing CMB data. Fortunately, doing so does not increase parameter errors significantly.

The NRAO VLA Sky Survey (NVSS) is the only dataset that allows an accurate determination of the auto-correlation function (ACF) on angular scales of several degrees for Active Galactic Nuclei (AGNs) at typical redshifts $z \simeq 1$. Surprisingly, the ACF is found to be positive on such large scales while, in the framework of the standard hierarchical clustering scenario with Gaussian primordial perturbations it should be negative for a redshift-independent effective halo mass of order of that found for optically-selected quasars. We show that a small primordial non-Gaussianity can add sufficient power on very large scales to account for the observed NVSS ACF. The best-fit value of the parameter $f_{\rm NL}$, quantifying the amplitude of primordial non-Gaussianity of local type is $f_{\rm NL}=62 \pm 27$ ($1\,\sigma$ error bar) and $25<f_{\rm NL}<117$ ($2\,\sigma$ confidence level), corresponding to a detection of non-Gaussianity significant at the $\sim 3\,\sigma$ confidence level. The minimal halo mass of NVSS sources is found to be $M_{\rm min}=10^{12.47\pm0.26}h^{-1}M_{\odot}$ ($1\,\sigma$) strikingly close to that found for optically selected quasars. We discuss caveats and possible physical and systematic effects that can impact on the results.

We update constraints on cosmic opacity by combining recent SN Type Ia data compilation with the latest measurements of the Hubble expansion at redshifts between 0 and 2. The new constraint on the parameter $\epsilon$ parametrising deviations from the luminosity-angular diameter distance relation ($d_L=d_A(1+z)^{2+\epsilon}$), is $\epsilon=-0.04_{-0.07}^{+0.08}$ (2-$\sigma$). For the redshift range between 0.2 and 0.35 this corresponds to an opacity $\Delta\tau<0.012$ (95% C.L.), a factor of 2 stronger than the previous constraint. Various models of beyond the standard model physics that predict violation of photon number conservation contribute to the opacity and can be equally constrained. In this paper we put new limits on axion-like particles, including chameleons, and mini-charged particles.

Cosmological probes are steadily reducing the total neutrino mass window, resulting in constraints on the neutrino-mass degeneracy as the most significant outcome. In this work we explore the discovery potential of cosmological probes to constrain the neutrino hierarchy, and point out some subtleties that could yield spurious claims of detection. This has an important implication for next generation of double beta decay experiments, that will be able to achieve a positive signal in the case of degenerate or inverted hierarchy of Majorana neutrinos. We find that cosmological experiments that nearly cover the whole sky could in principle distinguish the neutrino hierarchy by yielding 'substantial' evidence for one scenario over the another, via precise measurements of the shape of the matter power spectrum from large scale structure and weak gravitational lensing.

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.

With the advent of galaxy surveys which provide large samples of galaxies or galaxy clusters over a volume comparable to the horizon size (SDSS-III, HETDEX, Euclid, JDEM, LSST, Pan-STARRS, CIP etc.) or mass-selected large cluster samples over a large fraction of the extra-galactic sky (Planck, SPT, ACT, CMBPol, B-Pol), it is timely to investigate what constraints these surveys can impose on primordial non-Gaussianity. I illustrate here three different approaches: higher-order correlations of the three dimensional galaxy distribution, abundance of rare objects (extrema of the density distribution), and the large-scale clustering of halos (peaks of the density distribution). Each of these avenues has its own advantages, but, more importantly, these approaches are highly complementary under many respects.

We report on the finding of a correlation between galaxies' past star formation activity and the degree to which neighbouring galaxies rotation axes are aligned. This is obtained by cross-correlating star formation histories, derived with MOPED, and spin direction (chirality), as determined by the Galaxy Zoo project, for a sample of SDSS galaxies. Our findings suggest that spiral galaxies which formed the majority of their stars early (z > 2) tend to display coherent rotation over scales of ~10 Mpc/h. The correlation is weaker for galaxies with significant recent star formation. We find evidence for this alignment at more than the 5-sigma level, but no correlation with other galaxy stellar properties. This finding can be explained within the context of hierarchical tidal-torque theory if the SDSS galaxies harboring the majority of the old stellar population where formed in the past, in the same filament and at about the same time. Galaxies with significant recent star formation instead are in the field, thus influenced by the general tidal field that will align them in random directions or had a recent merger which would promote star formation, but deviate the spin direction.

We revisit the predictions for the expected Cosmic Microwave Background bispectrum signal from the primary-lensing-Rees-Sciama correlation; we point out that it can be a significant contaminant to the bispectrum signal from primordial non-Gaussianity of the local type. This non-Gaussianity, usually parameterized by the non-Gaussian parameter f_NL, arises, for example, in multi-field inflation. In particular both signals are frequency independent, and are maximized for nearly squeezed configurations. While their detailed scale-dependence and harmonic imprints are different for generic bispectrum shapes, we show that, if not included in the modeling, the primary-lensing-Rees-Sciama contribution yields an effective f_{NL} of 10 when using a bispectrum estimator optimized for local non-Gaussianity. Considering that expected 1-sigma errors on f_{NL} are < 10 from forthcoming experiments, we conclude that the contribution from this signal must be included in future constraints on f_{NL} from the Cosmic Microwave Background bispectrum.

We present a minimally-parametric reconstruction of the primordial power spectrum using the most recent cosmic microwave background and large scale structure data sets. Our goal is to constrain the shape of the power spectrum while simultaneously avoiding strong theoretical priors and over-fitting of the data. We find no evidence for any departure from a power law spectral index. We also find that an exact scale-invariant power spectrum is disfavored by the data, but this conclusion is weaker than the corresponding result assuming a theoretically-motivated power law spectral index prior. The reconstruction shows that better data are crucial to justify the adoption of such a strong theoretical prior observationally. These results can be used to determine the robustness of our present knowledge when compared with forthcoming precision data from Planck.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

The advent of large data-set in cosmology has meant that in the past 10 or 20 years our knowledge and understanding of the Universe has changed not only quantitatively but also, and most importantly, qualitatively. Cosmologists rely on data where a host of useful information is enclosed, but is encoded in a non-trivial way. The challenges in extracting this information must be overcome to make the most of a large experimental effort. Even after having converged to a standard cosmological model (the LCDM model) we should keep in mind that this model is described by 10 or more physical parameters and if we want to study deviations from it, the number of parameters is even larger. Dealing with such a high dimensional parameter space and finding parameters constraints is a challenge on itself. Cosmologists want to be able to compare and combine different data sets both for testing for possible disagreements (which could indicate new physics) and for improving parameter determinations. Finally, cosmologists in many cases want to find out, before actually doing the experiment, how much one would be able to learn from it. For all these reasons, sophisiticated statistical techniques are being employed in cosmology, and it has become crucial to know some statistical background to understand recent literature in the field. I will introduce some statistical tools that any cosmologist should know about in order to be able to understand recently published results from the analysis of cosmological data sets. I will not present a complete and rigorous introduction to statistics as there are several good books which are reported in the references. The reader should refer to those.

We consider the description of the clustering of halos for physically-motivated types of non-Gaussian initial conditions. In particular we include non-Gaussianity of the type arising from single field slow-roll, multi fields, curvaton (local type), higher-order derivative-type (equilateral), vacuum-state modifications (enfolded-type) and horizon-scale GR corrections type. We show that large-scale halo bias is a very sensitive tool to probe non-Gaussianity, potentially leading, for some planned surveys, to a detection of non-Gaussianity arising from horizon-scale GR corrections. In tandem with cosmic microwave background constraints, the halo-bias approach can help enormously to discriminate among different shapes of non-Gaussianity and thus among models for the origin of cosmological perturbations.

In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground minimum of ~100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization signal with tensor-to-scalar ratio r=0.01, in a small patch (~1%) of the sky known to have low Galactic emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects for detecting an r=0.01 signal including degree-scale measurements appear promising, with 5 sigma_r ~0.003 forecast from multiple methods. A mission that observes a range of scales offers better prospects from the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how optimizing the composition of frequency channels in the focal plane can maximize our ability to perform component separation, with a range of typically 40 < nu < 300 GHz preferred for ten channels. Foreground cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the optimization and detectability of the primordial signal will be useful for mission design.

The recent weak lensing measurement of the dark matter mass of the high-redshift galaxy cluster XMMUJ2235.3-2557 of (8.5 +- 1.7) x 10^{14} Msun at z=1.4, indicates that, if the cluster is assumed to be the result of the collapse of dark matter in a primordial gaussian field in the standard LCDM model, then its abundance should be < 0.002 clusters in the observed area. Here we investigate how to boost the probability of XMMUJ2235.3-2557 in particular resorting to deviations from Gaussian initial conditions. We show that this abundance can be boosted by factors > 3-10 if the non-Gaussianity parameter f^local_NL is in the range 150-200. This value is comparable to the limit for f_NL obtained by current constraints from the CMB. We conclude that mass determination of high-redshift, massive clusters can offer a complementary probe of primordial non-gaussianity.

Scaling relations among galaxy cluster observables, which will become available in large future samples of galaxy clusters, could be used to constrain not only cluster structure, but also cosmology. We study the utility of this approach, employing a physically motivated parametric model to describe cluster structure, and applying it to the expected relation between the Sunyaev-Zel'dovich decrement (S_\nu) and the emission-weighted X-ray temperature (T_ew). The slope and normalization of the entropy profile, the concentration of the dark matter potential, the pressure at the virial radius, and the level of non-thermal pressure support, as well as the mass and redshift-dependence of these quantities are described by free parameters. With a suitable choice of fiducial parameter values, the cluster model satisfies several existing observational constraints. We employ a Fisher matrix approach to estimate the joint errors on cosmological and cluster structure parameters from a measurement of S_\nu vs. T_ew in a future survey. We find that different cosmological parameters affect the scaling relation differently: predominantly through the baryon fraction (\Omega_m and \Omega_b), the virial overdensity (w_0 and w_a for low-z clusters) and the angular diameter distance (w_0, w_a for high-z clusters; \Omega_DE and h). We find that the cosmology constraints from the scaling relation are comparable to those expected from the number counts (dN/dz) of the same clusters. The scaling relation approach is relatively insensitive to selection effects and it offers a valuable consistency check; combining the information from the scaling relation and dN/dz is also useful to break parameter degeneracies and help disentangle cluster physics from cosmology.

We present new determinations of the cosmic expansion history from red-envelope galaxies. We have obtained for this purpose high-quality spectra with the Keck-LRIS spectrograph of red-envelope galaxies in 24 galaxy clusters in the redshift range 0.2 < z < 1.0. We complement these Keck spectra with high-quality, publicly available archival spectra from the SPICES and VVDS surveys. We improve over our previous expansion history measurements in Simon et al. (2005) by providing two new determinations of the expansion history: H(z) = 97 +- 62 km/sec/Mpc at z = 0.5 and H(z) = 90 +- 40 km/sec/Mpc at z = 0.8. We discuss the uncertainty in the expansion history determination that arises from uncertainties in the synthetic stellar-population models. We then use these new measurements in concert with cosmic-microwave-background (CMB) measurements to constrain cosmological parameters, with a special emphasis on dark-energy parameters and constraints to the curvature. In particular, we demonstrate the usefulness of direct H(z) measurements by constraining the dark- energy equation of state parameterized by w0 and wa and allowing for arbitrary curvature. Further, we also constrain, using only CMB and H(z) data, the number of relativistic degrees of freedom to be 4 +- 0.5 and their total mass to be < 0.2 eV, both at 1-sigma.

We present a spectroscopic catalog of (mostly) red galaxies in 24 galaxy clusters in the redshift range 0.17 < z < 0.92 obtained with the LRIS instrument on the Keck I telescope. Here we describe the observations and the galaxy spectra, including the discovery of three cD galaxies with LINER emission spectra, and the spectroscopic discovery of four new galaxy-galaxy lenses in cluster environments.

We present the power spectrum of the reconstructed halo density field derived from a sample of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey Seventh Data Release (DR7). The halo power spectrum has a direct connection to the underlying dark matter power for k <= 0.2 h/Mpc, well into the quasi-linear regime. This enables us to use a factor of ~8 more modes in the cosmological analysis than an analysis with kmax = 0.1 h/Mpc, as was adopted in the SDSS team analysis of the DR4 LRG sample (Tegmark et al. 2006). The observed halo power spectrum for 0.02 < k < 0.2 h/Mpc is well-fit by our model: chi^2 = 39.6 for 40 degrees of freedom for the best fit LCDM model. We find \Omega_m h^2 * (n_s/0.96)^0.13 = 0.141^{+0.009}_{-0.012} for a power law primordial power spectrum with spectral index n_s and \Omega_b h^2 = 0.02265 fixed, consistent with CMB measurements. The halo power spectrum also constrains the ratio of the comoving sound horizon at the baryon-drag epoch to an effective distance to z=0.35: r_s/D_V(0.35) = 0.1097^{+0.0039}_{-0.0042}. Combining the halo power spectrum measurement with the WMAP 5 year results, for the flat LCDM model we find \Omega_m = 0.289 +/- 0.019 and H_0 = 69.4 +/- 1.6 km/s/Mpc. Allowing for massive neutrinos in LCDM, we find \sum m_{\nu} < 0.62 eV at the 95% confidence level. If we instead consider the effective number of relativistic species Neff as a free parameter, we find Neff = 4.8^{+1.8}_{-1.7}. Combining also with the Kowalski et al. (2008) supernova sample, we find \Omega_{tot} = 1.011 +/- 0.009 and w = -0.99 +/- 0.11 for an open cosmology with constant dark energy equation of state w.

We explore consistency among different distance measures, including Supernovae Type Ia data, measurements of the Hubble parameter, and determination of the Baryon acoustic oscillation scale. We present new constraints on the cosmic transparency combining $H(z)$ data together with the latest Supernova Type Ia data compilation. This combination, in the context of a flat $\Lambda$CDM model, improves current constraints by nearly an order of magnitude. We re-examine the recently reported tension between the Baryon acoustic oscillation scale and Supernovae data in light of possible deviations from transparency, concluding that the source of the discrepancy may most likely be found among systematic effects of the modelling of the low redshift data or a simple $\sim 2-\sigma$ statistical fluke, rather than in exotic physics. Finally, we attempt to draw model-independent conclusions about the recent accelerated expansion, determining the acceleration redshift to be $z_{acc}=0.35^{+0.20}_{-0.13}$ (1-$\sigma$).

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We present first results from the Southern Cosmology Survey, a new multiwavelength survey of the southern sky coordinated with the Atacama Cosmology Telescope (ACT), a recently commissioned ground-based mm-band Cosmic Microwave Background experiment. This article presents a full analysis of archival optical multi-band imaging data covering an 8 square degree region near right ascension 23 hours and declination -55 degrees, obtained by the Blanco 4-m telescope and Mosaic-II camera in late 2005. We describe the pipeline we have developed to process this large data volume, obtain accurate photometric redshifts, and detect optical clusters. Our cluster finding process uses the combination of a matched spatial filter, photometric redshift probability distributions and richness estimation. We present photometric redshifts, richness estimates, luminosities, and masses for 8 new optically-selected clusters with mass greater than $3\times10^{14}M_{\sun}$ at redshifts out to 0.7. We also present estimates for the expected Sunyaev-Zel'dovich effect (SZE) signal from these clusters as specific predictions for upcoming observations by ACT, the South Pole Telescope and Atacama Pathfinder Experiment.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

The description of the abundance and clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys, which can potentially yield constraints of order unity on the non-Gaussianity parameter f_{NL}. We present tests on N-body simulations of analytical formulae describing the halo abundance and clustering for non-Gaussian initial conditions. We calibrate the analytic non-Gaussian mass function of Matarrese et al.(2000) and LoVerde et al.(2008) and the analytic description of clustering of halos for non-Gaussian initial conditions on N-body simulations. We find excellent agreement between the simulations and the analytic predictions if we make the corrections delta_c --> delta_c X sqrt{q} and delta_c --> \delta_c X q where q ~ 0.75, in the density threshold for gravitational collapse and in the non-Gaussian fractional correction to the halo bias, respectively. We discuss the implications of this correction on present and forecasted primordial non-Gaussianity constraints. We confirm that the non-Gaussian halo bias offers a robust and highly competitive test of primordial non-Gaussianity.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

Much of the science case for the next generation of deep, wide-field optical/infrared surveys has been driven by the further study of dark energy. This is a laudable goal (and the subject of a companion white paper by Zhan et al.). However, one of the most important lessons of the current generation of surveys is that the interesting science questions at the end of the survey are quite different than they were when the surveys were being planned. The current surveys succeeded in this evolving terrain by being very general tools that could be applied to a number of very fundamental measurements. Likewise, the accessibility of the data enabled the broader cosmological and astronomical community to generate more science than the survey collaborations could alone. With that in mind, we should consider some of the basic physical and cosmological questions that surveys like LSST and JDEM-Wide will be able to address.

The coming decade will be an exciting period for dark energy research, during which astronomers will address the question of what drives the accelerated cosmic expansion as first revealed by type Ia supernova (SN) distances, and confirmed by later observations. The mystery of dark energy poses a challenge of such magnitude that, as stated by the Dark Energy Task Force (DETF), "nothing short of a revolution in our understanding of fundamental physics will be required to achieve a full understanding of the cosmic acceleration." The lack of multiple complementary precision observations is a major obstacle in developing lines of attack for dark energy theory. This lack is precisely what next-generation surveys will address via the powerful techniques of weak lensing (WL) and baryon acoustic oscillations (BAO) -- galaxy correlations more generally -- in addition to SNe, cluster counts, and other probes of geometry and growth of structure. Because of their unprecedented statistical power, these surveys demand an accurate understanding of the observables and tight control of systematics. This white paper highlights the opportunities, approaches, prospects, and challenges relevant to dark energy studies with wide-deep multiwavelength photometric redshift surveys. Quantitative predictions are presented for a 20000 sq. deg. ground-based 6-band (ugrizy) survey with 5-sigma depth of r~27.5, i.e., a Stage 4 survey as defined by the DETF.

Baryon Acoustic Oscillations (BAO) in the radial direction offer a method to directly measure the Universe expansion history, and to set limits to space curvature when combined to the angular BAO signal. In addition to spectroscopic surveys, radial BAO might be measured from accurate enough photometric redshifts obtained with narrow-band filters. We explore the requirements for a photometric survey using Luminous Red Galaxies (LRG) to competitively measure the radial BAO signal and discuss the possible systematic errors of this approach. If LRG were a highly homogeneous population, we show that the photo-z accuracy would not substantially improve by increasing the number of filters beyond $\sim 10$, except for a small fraction of the sources detected at high signal-to-noise, and broad-band filters would suffice to achieve the target $\sigma_z = 0.003 (1+z)$ for measuring radial BAO. Using the LRG spectra obtained from SDSS, we find that the spectral variability of LRG substantially worsens the achievable photometric redshift errors, and that the optimal system consists of $\sim$ 30 filters of width $\Delta \lambda / \lambda \sim 0.02$. A $S/N > 20$ is generally necessary at the filters on the red side of the $H\alpha$ break to reach the target photometric accuracy. We estimate that a 5-year survey in a dedicated telescope with etendue in excess of 60 ${\rm m}^2 {\rm deg}^2$ would be necessary to obtain a high enough density of galaxies to measure radial BAO with sufficiently low shot noise up to $z= 0.85$. We conclude that spectroscopic surveys have a superior performance than photometric ones for measuring BAO in the radial direction.

Measurements of Baryonic Acoustic Oscillations in galaxy surveys have been recognized as a powerful tool for constraining dark energy. However, this method relies on the knowledge of the size of the acoustic horizon at recombination derived from Cosmic Microwave Background Anisotropy measurements. This estimate is typically derived assuming a standard recombination scheme; additional radiation sources can delay recombination altering the cosmic ionization history and the cosmological inferences drawn from CMB and BAO data. In this paper we quantify the effect of delayed recombination on the determination of dark energy parameters from future BAO surveys such as BOSS and WFMOS. We find the impact to be small but still not negligible. In particular, if recombination is non-standard (to a level still allowed by CMB data), but this is ignored, future surveys may incorrectly suggest the presence of a redshift dependent dark energy component. On the other hand, in the case of delayed recombination, adding to the analysis one extra parameter describing deviations from standard recombination, does not significantly degrade the error-bars on dark energy parameters and yields unbiased estimates.

We present catalogs of QSO candidates selected using photometry from GALEX combined with SDSS in the Stripe 82 region and Blanco Cosmology Survey (BCS) near declination -55 degrees. The SDSS region contains ~700 objects with magnitude i < 20 and ~3600 objects with i < 21.5 in a ~60 square degree sky region, while the BCS region contains ~280 objects with magnitude i < 20 and ~2000 objects with i < 21.5 for a 11 square degree sky region that is being observed by three current microwave Sunyaev-Zeldovich surveys. Our QSO catalog is the first one in the BCS region. Deep GALEX exposures (~2000 seconds in FUV and NUV, except in three fields) provide high signal-to-noise photometry in the GALEX bands (FUV, NUV < 24.5 mag). From this data, we select QSO candidates using only GALEX and optical r-band photometry, using the method given by Atlee and Gould (2008). In the Stripe 82 field, 60% (30%) of the GALEX selected QSO's with optical magnitude i<20 (i<21.5) also appear in the Richards et al. (2008) QSO catalog constructed using 5-band optical SDSS photometry. Comparison with the same catalog by Richards et al. shows that the completeness of the sample is approximately 40%(25%). However, for regions of the sky with very low dust extinction, like the BCS 23hr field and the Stripe 82 between 0 and 10 degrees in RA, our completeness is close to 95%, demonstrating that deep GALEX observations are almost as efficient as multi-wavelength observations at finding QSO's. GALEX observations thus provide a viable alternate route to QSO catalogs in sky regions where u-band optical photometry is not available. The full catalog is available at http://www.ice.csic.es/personal/jimenez/PHOTOZ

We summarize existing constraints on the epoch of reionization and discuss the observational probes that are sensitive to the process. We focus on the role large scale polarization can play. Polarization probes the integrated optical depth across the entire epoch of reionization. Future missions such as Planck and CMBPol will greatly enhance our knowledge of the reionization history, allowing us to measure the time evolution of the ionization fraction. As large scale polarization probes high redshift activity, it can best constrain models where the Universe was fully or partially ionized at early times. In fact, large scale polarization could be our only probe of the highest redshifts.

We discuss the cosmological degeneracy between the Hubble parameter H(z), the age of the universe and cosmological parameters describing simple variations from the minimal LCDM model. We show that independent determinations of the Hubble parameter H(z) such as those recently provided by Simon et al., Phys. Rev. D 71 (2005) 123001, combined with Cosmic Microwave Background data, can provide stringent constraints on possible deviations from the LCDM model. In particular we find that this data combination constrains at the 68% (95%) confidence level the following parameters: sum of the neutrino masses < 0.5 (1.0) eV, number of relativistic neutrino species N_rel = 4.1^{+0.4}_{-0.9} (^{+1.1}_{-1.5}), dark energy equation of state parameter w = -0.95 +- 0.17 (+- 0.32), and curvature \Omega_k = 0.002 +- 0.006 (+- 0.014), in excellent agreement with dataset combinations involving Cosmic Microwave Background, Supernovae and Baryon Acoustic Oscillations. This offers a valuable consistency check for systematic errors.

We use a Press-Schechter-like calculation to study how the abundance of voids changes in models with non-Gaussian initial conditions. While a positive skewness increases the cluster abundance, a negative skewness does the same for the void abundance. We determine the dependence of the void abundance on the non-Gaussianity parameter fnl for the local-model bispectrum-which approximates the bispectrum in some multi-field inflation models-and for the equilateral bispectrum, which approximates the bispectrum in e.g. string-inspired DBI models of inflation. We show that the void abundance in large-scale-structure surveys currently being considered should probe values as small as fnl < 10 and fnl^eq < 30, over distance scales ~10 Mpc.

We generalize to non-flat geometries the formalism of Simon et al. (2005) to reconstruct the dark energy potential. This formalism makes use of quantities similar to the Horizon-flow parameters in inflation, can, in principle, be made non-parametric and is general enough to be applied outside the simple, single scalar field quintessence. Since presently available and forthcoming data do not allow a non-parametric and exact reconstruction of the potential, we consider a general parametric description in term of Chebyshev polynomials. We then consider present and future measurements of H(z), Baryon Acoustic Oscillations surveys and Supernovae type 1A surveys, and investigate their constraints on the dark energy potential. We find that, relaxing the flatness assumption increases the errors on the reconstructed dark energy evolution but does not open up significant degeneracies, provided that a modest prior on geometry is imposed. Direct measurements of H(z), such as those provided by BAO surveys, are crucially important to constrain the evolution of the dark energy potential and the dark energy equation of state, especially for non-trivial deviations from the standard LambdaCDM model.

The latest cosmological data seem to indicate a significant deviation from scale invariance of the primordial power spectrum when parameterized either by a power law or by a spectral index with non-zero "running". This deviation, by itself, serves as a powerful tool to discriminate among theories for the origin of cosmological structures such as inflationary models. Here, we use a minimally-parametric smoothing spline technique to reconstruct the shape of the primordial power spectrum. This technique is well-suited to search for smooth features in the primordial power spectrum such as deviations from scale invariance or a running spectral index, although it would recover sharp features of high statistical significance. We use the WMAP 3 year results in combination with data from a suite of higher resolution CMB experiments (including the latest ACBAR 2008 release), as well as large-scale structure data from SDSS and 2dFGRS. We employ cross-validation to assess, using the data themselves, the optimal amount of smoothness in the primordial power spectrum consistent with the data. This minimally-parametric reconstruction supports the evidence for a power law primordial power spectrum with a red tilt, but not for deviations from a power law power spectrum. Smooth variations in the primordial power spectrum are not significantly degenerate with the other cosmological parameters.

We forecast constraints on primordial non-Gaussianity achievable from forthcoming surveys by exploiting the scale-dependent halo bias introduced on large scales by non-Gaussian initial conditions. We explore the performance of exploiting both the shape of the galaxy power-spectrum on large scales and the cross-correlation of galaxies with Cosmic Microwave Background maps through the Integrated Sachs Wolfe effect. We find that future surveys can detect primordial non-Gaussianity of the local form with a non-Gaussianity parameter $|f_{\rm NL}|$ of order unity. This is particularly exciting because, while the simplest single-field slow-roll models of inflation predict a primordial $|f_{\rm NL}|\ll 1$, this signal sources extra contributions to the effective $f_{\rm NL}$ of large-scale structures that are expected to be above our predicted detection threshold.

Four large-area Sunyaev-Zeldovich (SZ) experiments -- APEX-SZ, SPT, ACT, and Planck -- promise to detect clusters of galaxies through the distortion of Cosmic Microwave Background photons by hot (> 10^6 K) cluster gas (the SZ effect) over thousands of square degrees. A large observational follow-up effort to obtain redshifts for these SZ-detected clusters is under way. Given the large area covered by these surveys, most of the redshifts will be obtained via the photometric redshift (photo-z) technique. Here we demonstrate, in an application using ~3000 SDSS stripe 82 galaxies with r<20, how the addition of GALEX photometry (FUV, NUV) greatly improves the photometric redshifts of galaxies obtained with optical griz or ugriz photometry. In the case where large spectroscopic training sets are available, empirical neural-network-based techniques (e.g., ANNz) can yield a photo-z scatter of $\sigma_z = 0.018 (1+z)$. If large spectroscopic training sets are not available, the addition of GALEX data makes possible the use simple maximum likelihood techniques, without resorting to Bayesian priors, and obtains $\sigma_z=0.04(1+z)$, accuracy that approaches the accuracy obtained using spectroscopic training of neural networks on ugriz observations. This improvement is especially notable for blue galaxies. To achieve these results, we have developed a new set of high resolution spectral templates based on physical information about the star formation history of galaxies. We envision these templates to be useful for the next generation of photo-z applications. We make our spectral templates and new photo-z catalogs available to the community at http://www.ice.csic.es/personal/jimenez/PHOTOZ .

The detection of primordial non-Gaussianity could provide a powerful means to test various inflationary scenarios. Although scale-invariant non-Gaussianity (often described by the $f_{NL}$ formalism) is currently best constrained by the CMB, single-field models with changing sound speed can have strongly scale-dependent non-Gaussianity. Such models could evade the CMB constraints but still have important effects at scales responsible for the formation of cosmological objects such as clusters and galaxies. We compute the effect of scale-dependent primordial non-Gaussianity on cluster number counts as a function of redshift, using a simple ansatz to model scale-dependent features. We forecast constraints on these models achievable with forthcoming data sets. We also examine consequences for the galaxy bispectrum. Our results are relevant for the Dirac-Born-Infeld model of brane inflation, where the scale-dependence of the non-Gaussianity is directly related to the geometry of the extra dimensions.

We discuss the cosmological degeneracy between the age of the Universe, the Hubble parameter and the effective number of relativistic particles N_eff. We show that independent determinations of the Hubble parameter H(z) as those recently provided by Simon,Verde, Jimenez (2006), combined with other cosmological data sets can provide the most stringent constraint on N_eff, yielding N_eff=3.7 (-1.2) (+1.1) at 95% confidence level. A neutrino background is detected with high significance: N_eff >1.8 at better than 99% confidence level. Constraints on the age of the universe in the framework of an extra background of relativistic particles are improved by a factor 3.

We show that cosmological observables can constrain the topology of the compact additional dimensions predicted by string theory. To do this, we develop a general strategy for relating cosmological observables to the microscopic parameters of the potentials and field-dependent kinetic terms of the multiple scalar fields that arise in the low-energy limit of string theory. We apply this formalism to the Large Volume Scenarios in Type IIB flux compactifications where analytical calculations are possible. Our methods generalize to other settings.

It has long been known how to analytically relate the clustering properties of the collapsed structures (halos) to those of the underlying dark matter distribution for Gaussian initial conditions. Here we apply the same approach to physically motivated non-Gaussian models. The techniques we use were developed in the 1980s to deal with the clustering of peaks of non-Gaussian density fields. The description of the clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys. For inflationary-motivated non-Gaussianites, we find an analytic expression for the halo bias as a function of scale, mass and redshift, employing only the approximations of high-peaks and large separations.

In this paper we investigate the potential of 3D cosmic shear to constrain massive neutrino parameters. We find that if the total mass is substantial (near the upper limits from LSS, but setting aside the Ly alpha limit for now), then 3D cosmic shear + Planck is very sensitive to neutrino mass and one may expect that a next generation photometric redshift survey could constrain the number of neutrinos N_nu and the sum of their masses m_nu to an accuracy of dN_nu ~ 0.08 and dm_nu ~ 0.03 eV respectively. If in fact the masses are close to zero, then the errors weaken to dN_nu ~ 0.10 and dm_nu~0.07 eV. In either case there is a factor 4 improvement over Planck alone. We use a Bayesian evidence method to predict joint expected evidence for N_nu and m_nu. We find that 3D cosmic shear combined with a Planck prior could provide `substantial' evidence for massive neutrinos and be able to distinguish `decisively' between many competing massive neutrino models. This technique should `decisively' distinguish between models in which there are no massive neutrinos and models in which there are massive neutrinos with |N_nu-3| > 0.35 and m_nu > 0.25 eV. We introduce the notion of marginalised and conditional evidence when considering evidence for individual parameter values within a multi-parameter model.

This is the summary of 4 Lectures given at the XIX Canary islands winter school of Astrophysics ''The Cosmic Microwave Background, from Quantum Fluctuations to the present Universe''. Lectures were intended for master/beginning PhD students in cosmology.

We analyze the Jordan-Brans-Dicke model (JBD) of gravity, where deviations from General Relativity (GR) are described by a scalar field non-minimally coupled to gravity. The theory is characterized by a constant coupling parameter, $\omega_{\rm JBD}$; GR is recovered in the limit $\omega_{\rm JBD} \to \infty$. In such theories, gravity modifications manifest at early times, so that one cannot rely on the usual approach of looking for inconsistencies in the expansion history and perturbations growth in order to discriminate between JBD and GR. However, we show that a similar technique can be successfully applied to early and late times observables instead. Cosmological parameters inferred extrapolating early-time observations to the present will match those recovered from direct late-time observations only if the correct gravity theory is used. We use the primary CMB, as will be seen by the Planck satellite, as the early-time observable; and forthcoming and planned Supernov{\ae}, Baryonic Acoustic Oscillations and Weak Lensing experiments as late-time observables. We find that detection of values of $\omega_{\rm JBD}$ as large as 500 and 1000 is within reach of the upcoming (2010) and next-generation (2020) experiments, respectively.

At the epoch of reionization, when the high-redshift inter-galactic medium (IGM) is being enriched with metals, the 63.2 micron fine structure line of OI is pumped by the ~ 1300 AA soft UV background and introduces a spectral distortion in the Cosmic Microwave Background (CMB). Here we use a toy model for the spatial distribution of neutral oxygen, assuming metal bubbles surround dark matter halos, and compute the fluctuations of this distortion, and the angular power spectrum it imprints on the CMB. We discuss the dependence of the power spectrum on the velocity of the winds polluting the IGM with metals, the minimum mass of the halos producing these winds, and on the cosmic epoch when the OI pumping occurs. We find that, although the clustering signal of the CMB distortion is weak \delta y_{rms} ~ 10^{-7} (roughly corresponding to a temperature anisotropy of few nK), it may be reachable in deep integrations with high-sensitivity infrared detectors. Even without a detection, these instruments should be able to useful constraints on the heavy element enrichment history of the IGM.

The Fisher matrix approach (Fisher 1935) allows one to calculate in advance how well a given experiment will be able to estimate model parameters, and has been an invaluable tool in experimental design. In the same spirit, we present here a method to predict how well a given experiment can distinguish between different models, regardless of their parameters. From a Bayesian viewpoint, this involves computation of the Bayesian evidence. In this paper, we generalise the Fisher matrix approach from the context of parameter fitting to that of model testing, and show how the expected evidence can be computed under the same simplifying assumption of a gaussian likelihood as the Fisher matrix approach for parameter estimation. With this `Laplace approximation' all that is needed to compute the expected evidence is the Fisher matrix itself. We illustrate the method with a study of how well upcoming and planned experiments should perform at distinguishing between Dark Energy models and modified gravity theories. In particular we consider the combination of 3D weak lensing, for which planned and proposed wide-field multi-band imaging surveys will provide suitable data, and probes of the expansion history of the Universe, such as proposed supernova and baryonic acoustic oscillations surveys. We find that proposed large-scale weak lensing surveys from space should be able readily to distinguish General Relativity from modified gravity models.

We compute the effect of local electrons on the CMB temperature anisotropies. The number density and distribution of free electrons in our Galaxy has been accurately measured from pulsar dispersion measurements. Because of their distribution, the dynamics of our Galaxy and the Galaxy peculiar velocity with respect to the Hubble flow, these free electrons leave a frequency-independent imprint on the cosmic microwave background (CMB). In particular, the coherent motion of the free electrons respect to us and to the CMB rest frame produce a kinetic Sunyaev-Zeldovich signal. We compute this effect and we note that the large-scale antisymmetry of the signal gives it an angular power spectrum with a sawtooth pattern where even multipoles are suppressed with respect to the odd ones. We find the signal to be small (~ 2 uK) and sub-dominant compared to the primary CMB and other foreground signals. However, since there are no free parameters in the modeling of this signal, it can be taken into account if more precise measurements of the primordial signal are required.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

We consider the pumping of the $63.2 \mu$m fine structure line of neutral OI in the high--redshift intergalactic medium (IGM), in analogy with the Wouthuysen--Field effect for the 21cm line of cosmic HI. We show that the soft UV background at $\sim 1300$\AA can affect the population levels, and if a significant fraction of the IGM volume is filled with ``fossil HII regions'' containing neutral OI, then this can produce a non--negligible spectral distortion in the cosmic microwave background (CMB). OI from redshift $z$ is seen in emission at $(1+z)63.2\mu$m, and between $7<z<10$ produces a mean spectral distortion of the CMB with a $y$--parameter of $y=(10^{-9} - 3\times10^{-8}) (Z/10^{-3}{\rm Z_{\odot}}) (I_{UV})$, where $Z$ is the mean metallicity of the IGM and $I_{UV}$ is the UV background at 1300\AA in units of $10^{-20}$ erg/s/Hz/cm$^2$/sr. Because O is in charge exchange equilibrium with H, a measurement of this signature can trace the metallicity at the end of the dark ages, prior to the completion of cosmic reionization and is complementary to cosmological 21cm studies. While future CMB experiments, such as Planck could constrain the metallicity to the $10^{-2} Z_{\odot}$ level, specifically designed experiments could potentially achieve a detection. Fluctuations of the distortion on small angular scale may also be detectable.

Using cosmological hydrodynamical simulations, we investigate the prospects of the thermal Sunyaev-Zel'dovich (tSZ) effect to detect the missing baryons in the local universe. We find that at least 80% of the tSZ luminosity is generated in collapsed structures, and that $\sim$ 70% of the remaining diffuse tSZ luminosity (i.e., $\sim 15$% of the total) comes from overdense regions with $\delta_{gas}>$10, such as filaments and superclusters. The gas present in slightly overdense and underdense regions with $\delta_{gas} < 10$, despite making up 50% of the total baryon budget, leaves very little tSZ signature: it gives rise to only $\sim$ 5% of the total tSZ luminosity. Thus, future Cosmic Microwave Background (CMB) observations will be sensitive to, at best, one half of the missing baryons, improving the current observational status, but still leaving one half unobserved. Since most of the tSZ is generated in haloes, we find a tight correlation between gas pressure and galaxy number density. This allows us to predict the CMB Comptonization from existing galaxy surveys and to forecast the tSZ effect from the local structures probed by the Two Micron All Sky Survey (2MASS) galaxy catalog.

A four-dimensional universe, arising from a flux compactification of Type IIB string theory, contains scalar fields with a potential determined by topological and geometric parameters of the internal -hidden- dimensions. We show that inflation can be realized via rolling towards the large internal volume minima that are generic in these scenarios, and we give explicit formulae relating the microscopic parameters (e.g., the Euler number of the internal space) to the cosmological observables (e.g., the spectral index). We find that the tensor-to-scalar ratio, the running of the spectral index, and the potential energy density at the minimum are related by consistency relations and are exponentially small in the number of e-foldings. Further, requiring that these models arise as low-energy limits of string theory eliminates most of them, even if they are phenomenologically valid. In this context, this approach provides a strategy for systematically falsifying stringy inflation models.

We study the constraints that future multifrequency Cosmic Microwave Background (CMB) experiments will be able to set on the metal enrichment history of the Inter Galactic Medium at the epoch of reionisation. We forecast the signal to noise ratio for the detection of the signal introduced in the CMB by resonant scattering off metals at the end of the Dark Ages. We take into account systematics associated to inter-channel calibration, PSF reconstruction errors and innacurate foreground removal. We develop an algorithm to optimally extract the signal generated by metals during reionisation and to remove accurately the contamination due to the thermal Sunyaev-Zel'dovich effect. Although demanding levels of foreground characterisation and control of systematics are required, they are very distinct from those encountered in HI-21cm studies and CMB polarization, and this fact encourages the study of resonant scattering off metals as an alternative way of conducting tomography of the reionisation epoch. An ACT-like experiment with optimistic assumtions on systematic effects, and looking at clean regions of the sky, can detect changes of 3%-12% (95% c.l.) of the OIII abundance (with respect its solar value) in the redshift range $z\in$ [12,22], for reionization redshift $z_{\rm re}>10$. However, for $z_{\rm re} <10$, it can only set upper limits on NII abundance increments of $\sim$ 60% its solar value in the redshift range $z\in$ [5.5,9], (95% c.l.). These constraints assume that inter-channel calibration is accurate down to one part in $10^{4}$, which constitutes the most critical technical requirement of this method, but still achievable with current technology.

We quantify the limiting factors in optimizing current-technology cosmic microwave background (CMB) polarization experiments in order to learn about inflationary physics. We consider space-based, balloon-borne and ground-based experiments. We find that foreground contamination and residuals from foreground subtraction are ultimately the limiting factors in detecting a primordial gravity wave signal. For full-sky space-based experiments, these factors hinder the detection of tensor-to-scalar ratios of r < 10^{-3} on large scales, while for ground-based experiments these factors impede the ability to apply delensing techniques. We consider ground-based/balloon-borne experiments of currently planned or proposed designs and find that it is possible for a value of r=0.01 to be measured at ~ 3-sigma level. A small space-based CMB polarization experiment, with current detector technology and full sky coverage, can detect r ~ 10^{-3} at the ~ 3-sigma level, but a markedly improved knowledge of polarized foregrounds is needed. We advocate using as wide a frequency coverage as possible in order to carry out foreground subtraction at the percent level, which is necessary to measure such a small primordial tensor amplitude. To produce a clearly detectable (>3-sigma) tensor component in a realistic CMB experiment, inflation must either involve large-field variations, \Delta\phi >~ 1 or multi-field/hybrid models. Hybrid models can be easily distinguished from large-field models due to their blue scalar spectral index. Therefore, an observation of a tensor/scalar ratio and n < 1 in future experiments with the characteristics considered here may be an indication that inflation is being driven by some physics in which the inflaton cannot be described as a fundamental field.

We show how to place constraints on cluster physics by stacking the weak lensing signals from multiple clusters found through the Sunyaev-Zeldovich (SZ) effect. For a survey that covers about 200 sq. deg. both in SZ and weak lensing observations, the slope and amplitude of the mass vs. SZ luminosity relation can be measured with few percent error for clusters at z~0.5. This can be used to constrain cluster physics, such as the nature of feedback. For example, we can distinguish a pre-heated model from a model with a decreased accretion rate at more than 5sigma. The power to discriminate among different non-gravitational processes in the ICM becomes even stronger if we use the central Compton parameter y_0, which could allow one to distinguish between models with pre-heating, SN feedback and AGN feedback, for example, at more than 5sigma. Measurement of these scaling relations as a function of redshift makes it possible to directly observe e.g., the evolution of the hot gas in clusters. With this approach the mass-L_SZ relation can be calibrated and its uncertainties can be quantified, leading to a more robust determination of cosmological parameters from clusters surveys. The mass-L_SZ relation calibrated in this way from a small area of the sky can be used to determine masses of SZ clusters from very large SZ-only surveys and is nicely complementary to other techniques proposed in the literature.

In the context of a cosmological study of the bulk flows in the Universe, we present a detailed study of the statistical properties of the kinematic Sunyaev-Zel'dovich (kSZ) effect. We first compute analytically the correlation function and the power spectrum of the projected peculiar velocities of galaxy clusters. By taking into account the spatial clustering properties of these sources, we perform a line-of-sight computation of the {\em all-sky} kSZ power spectrum and find that at large angular scales ($l<10$), the local bulk flow should leave a visible signature above the Poisson-like fluctuations dominant at smaller scales, while the coupling of density and velocity fluctuations should give much smaller contribution. We conduct an analysis of the prospects of future high resolution CMB experiments (such as ACT and SPT) to detect the kSZ signal and to extract cosmological information and dark energy constraints from it. We present two complementary methods, one suitable for ``deep and narrow'' surveys such as ACT and one suitable for ``wide and shallow'' surveys such as SPT. Both methods can constraint the equation of state of dark energy $w$ to about 5-10% when applied to forthcoming and future surveys, and probe $w$ in complementary redshift ranges, which could shed some light on its time evolution. These determinations of $w$ do not rely on the knowledge of cluster masses, although they make minimal assumptions on cluster physics.

We reconstruct the shape of the primordial power spectrum (PPS) using a smoothing spline. Our adapted smoothing spline technique provides a complementary method to existing efforts to search for smooth features in the PPS, such as a running spectral index. With this technique we find no significant indication with WMAP first-year data that the PPS deviates from Harrison-Zeldovich and no evidence for loss of power on large scales. We also examine the effect on the cosmological parameters of the additional PPS freedom. Smooth variations in the PPS are not significantly degenerate with other cosmological parameters, but the spline reconstruction greatly increases the errors on the optical depth and baryon fraction.

We present an attempt to constrain deviations from the gravitational inverse-square law on large-scale structure scales. A perturbed law modifies the Poisson equation, which implies a scale-dependent growth of overdensities in the linear regime and thus modifies the power spectrum shape. We use two large-scale structure surveys (the Sloan Digital Sky survey and the Anglo-Australian Two-degree field galaxy redshift survey) to constrain the parameters of two simple modifications of the inverse-square law. We find no evidence for deviations from normal gravity on the scales probed by these surveys (~ 10^(23) m.)

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

We develop a formalism to characterize the redshift evolution of the dark energy potential. Our formalism makes use of quantities similar to the Horizon-flow parameters in inflation and is general enough that can deal with multiscalar quintessence scenarios, exotic matter components, and higher order curvature corrections to General Relativity. We show how the shape of the dark energy potential can be recovered non parametrically using this formalism and we present approximations analogous to the ones relevant to slow-roll inflation. Since presently available data do not allow a non-parametric and exact reconstruction of the potential, we consider a general parametric description. This reconstruction can also be used in other approaches followed in the literature (e.g., the reconstruction of the redshift evolution of the dark energy equation of state w(z)). Using observations of passively evolving galaxies and supernova data we derive constraints on the dark energy potential shape in the redshift range 0.1 < z < 1.8. Our findings show that at the 1sigma level the potential is consistent with being constant, although at the same level of confidence variations cannot be excluded with current data. We forecast constraints achievable with future data from the Atacama Cosmology Telescope.

Using the spectroscopic sample of the SDSS DR1 we measure how gas was transformed into stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a clear correlation between early star formation activity and present-day stellar mass: the more massive galaxies have formed about 80% of their stars at $z>1$, while for the less massive ones the value is only about 20%. By comparing the BCT to the dark matter merger tree, we find indications that star formation efficiency at $z>1$ had to be about a factor of two higher than today ($\sim 10%$) in galaxies with present-day stellar mass larger than $2 \times 10^{11}M_\odot$, if this early star formation occurred in the main progenitor. Therefore, the LCDM paradigm can accommodate a large number of red objects. On the other hand, in galaxies with present-day stellar mass less than $10^{11}$ M$_{\odot}$, efficient star formation seems to have been triggered at $z \sim 0.2$. We show that there is a characteristic mass (M$_* \sim 10^{10}$ M$_{\odot}$) for feedback efficiency (or lack of star formation). For galaxies with masses lower than this, feedback (or star formation suppression) is very efficient while for higher masses it is not. The BCT, determined here for the first time, should be an important observable with which to confront theoretical

We improve the algorithm of Kosowsky, Milosavljevic, and Jimenez (2002) for computing power spectra of the cosmic microwave background. The present algorithm computes not only the temperature power spectrum but also the E-mode polarization and the temperature-polarization cross power spectra, providing the accuracy required for current cosmological parameter estimation. We refine the optimum set of cosmological parameters for computing the power spectra as perturbations around a fiducial model, leading to an accuracy better than 0.5% for the temperature power spectrum throughout the region of parameter space within WMAP's first-year 3sigma confidence region. This accuracy is comparable to the difference between the widely-used CMBFAST code (Seljak and Zaldarriaga 1996) and Boltzmann codes. Our algorithm (CMBwarp) makes possible a full exploration of the likelihood region for eight cosmological parameters in about one hour on a laptop computer. We provide the code to compute power spectra as well as the Markov Chain Monte Carlo algorithm for cosmological parameters estimation at http://www.physics.upenn.edu/~raulj/CMBwarp

We measure the large-scale real-space power spectrum P(k) using a sample of 205,443 galaxies from the Sloan Digital Sky Survey, covering 2417 square degrees with mean redshift z~0.1. We employ a matrix-based method using pseudo-Karhunen-Loeve eigenmodes, producing uncorrelated minimum-variance measurements in 22 k-bands of both the clustering power and its anisotropy due to redshift-space distortions, with narrow and well-behaved window functions in the range 0.02 h/Mpc < k < 0.3h/Mpc. We pay particular attention to modeling, quantifying and correcting for potential systematic errors, nonlinear redshift distortions and the artificial red-tilt caused by luminosity-dependent bias. Our final result is a measurement of the real-space matter power spectrum P(k) up to an unknown overall multiplicative bias factor. Our calculations suggest that this bias factor is independent of scale to better than a few percent for k<0.1h/Mpc, thereby making our results useful for precision measurements of cosmological parameters in conjunction with data from other experiments such as the WMAP satellite. As a simple characterization of the data, our measurements are well fit by a flat scale-invariant adiabatic cosmological model with h Omega_m =0.201+/- 0.017 and L* galaxy sigma_8=0.89 +/- 0.02 when fixing the baryon fraction Omega_b/Omega_m=0.17 and the Hubble parameter h=0.72; cosmological interpretation is given in a companion paper.

We extend the Fourier transform based method for the analysis of galaxy redshift surveys of Feldman, Kaiser & Peacock (1994: FKP) to model luminosity-dependent clustering. In a magnitude limited survey, galaxies at high redshift are more luminous on average than galaxies at low redshift. Galaxy clustering is observed to increase with luminosity, so the inferred density field is effectively multiplied by an increasing function of radius. This has the potential to distort the shape of the recovered power spectrum. In this paper we present an extension of the FKP analysis method to incorporate this effect, and present revised optimal weights to maximize the precision of such an analysis. The method is tested and its accuracy assessed using mock catalogues of the 2-degree field galaxy redshift survey (2dFGRS). We also show that the systematic effect caused by ignoring luminosity-dependent bias was negligible for the initial analysis of the 2dFGRS of Percival et al. (2001). However, future surveys, sensitive to larger scales, or covering a wider range of galaxy luminosities will benefit from this refined method.

We describe our methodology for comparing the WMAP measurements of the cosmic microwave background (CMB) and other complementary data sets to theoretical models. The unprecedented quality of the WMAP data, and the tight constraints on cosmological parameters that are derived, require a rigorous analysis so that the approximations made in the modeling do not lead to significant biases. We describe our use of the likelihood function to characterize the statistical properties of the microwave background sky. We outline the use of the Monte Carlo Markov Chains to explore the likelihood of the data given a model to determine the best fit cosmological parameters and their uncertainties. We add to the WMAP data the l>~700 CBI and ACBAR measurements of the CMB, the galaxy power spectrum at z~0 obtained from the 2dF galaxy redshift survey (2dFGRS), and the matter power spectrum at z~3 as measured with the Ly alpha forest. These last two data sets complement the CMB measurements by probing the matter power spectrum of the nearby universe. Combining CMB and 2dFGRS requires that we include in our analysis a model for galaxy bias, redshift distortions, and the non-linear growth of structure. We show how the statistical and systematic uncertainties in the model and the data are propagated through the full analysis.

We present limits to the amplitude of non-Gaussian primordial fluctuations in the WMAP 1-year cosmic microwave background sky maps. A non-linear coupling parameter, f_NL, characterizes the amplitude of a quadratic term in the primordial potential. We use two statistics: one is a cubic statistic which measures phase correlations of temperature fluctuations after combining all configurations of the angular bispectrum. The other uses the Minkowski functionals to measure the morphology of the sky maps. Both methods find the WMAP data consistent with Gaussian primordial fluctuations and establish limits, -58<f_NL<134, at 95% confidence. There is no significant frequency or scale dependence of f_NL. The WMAP limit is 30 times better than COBE, and validates that the power spectrum can fully characterize statistical properties of CMB anisotropy in the WMAP data to high degree of accuracy. Our results also validate the use of a Gaussian theory for predicting the abundance of clusters in the local universe. We detect a point-source contribution to the bispectrum at 41 GHz, b_src = (9.5+-4.4) X 1e-5 uK^3 sr^2, which gives a power spectrum from point sources of c_src = (15+-6) X 1e-3 uK^2 sr in thermodynamic temperature units. This value agrees well with independent estimates of source number counts and the power spectrum at 41 GHz, indicating that b_src directly measures residual source contributions.

WMAP precision data enables accurate testing of cosmological models. We find that the emerging standard model of cosmology, a flat Lambda-dominated universe seeded by nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. With parameters fixed only by WMAP data, we can fit finer scale CMB measurements and measurements of large scle structure (galaxy surveys and the Lyman alpha forest). This simple model is also consistent with a host of other astronomical measurements. We then fit the model parameters to a combination of WMAP data with other finer scale CMB experiments (ACBAR and CBI), 2dFGRS measurements and Lyman alpha forest data to find the model's best fit cosmological parameters: h=0.71+0.04-0.03, Omega_b h^2=0.0224+-0.0009, Omega_m h^2=0.135+0.008-0.009, tau=0.17+-0.06, n_s(0.05/Mpc)=0.93+-0.03, and sigma_8=0.84+-0.04. WMAP's best determination of tau=0.17+-0.04 arises directly from the TE data and not from this model fit, but they are consistent. These parameters imply that the age of the universe is 13.7+-0.2 Gyr. The data favors but does not require a slowly varying spectral index. By combining WMAP data with other astronomical data sets, we constrain the geometry of the universe, Omega_tot = 1.02 +- 0.02, the equation of state of the dark energy w < -0.78 (95% confidence limit assuming w >= -1), and the energy density in stable neutrinos, Omega_nu h^2 < 0.0076 (95% confidence limit). For 3 degenerate neutrino species, this limit implies that their mass is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization rules out warm dark matter.

Since this is a meeting on inflation the title of this talk should probably have been ``The cosmos as a lab for inflation''. I will illustrate how, by combining WMAP observations of the cosmic microwave background (CMB) with smaller scales CMB experiments and large-scale structure surveys, it is possible to constrain inflationary models. WMAP data in combination with these external data sets, offers a window to the very early universe and the physics of inflation: enables us to confirm the basic tenets of the inflationary paradigm and begin to quantitatively test inflationary models. Understanding the complex physics and complicated systematics of the low redshift observations is crucial to make progress in this direction.

We present full sky microwave maps in five bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 2<l<900 anisotropy power spectrum is cosmic variance limited for l<354 with a signal-to-noise ratio >1 per mode to l=658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is 0.17 +/- 0.04, which implies a reionization epoch of 180+220-80 Myr (95% CL) after the Big Bang at a redshift of 20+10-9 (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. The age of the best-fit universe is 13.7 +/- 0.2 Gyr old. Decoupling was 379+8-7 kyr after the Big Bang at a redshift of 1089 +/- 1. The thickness of the decoupling surface was dz=195 +/- 2. The matter density is Omega_m h^2 = 0.135 +0.008 -0.009, the baryon density is Omega_b h^2 = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega_tot = 1.02 +/- 0.02. The spectral index of scalar fluctuations is fit as n_s = 0.93 +/- 0.03 at wavenumber k_0 = 0.05 Mpc^-1, with a running index slope of dn_s/d ln k = -0.031 +0.016 -0.018 in the best-fit model. This flat universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy equation of state is limited to w<-0.78 (95% CL). Inflation theory is supported with n_s~1, Omega_tot~1, Gaussian random phases of the CMB anisotropy, and superhorizon fluctuations. An admixture of isocurvature modes does not improve the fit. The tensor-to-scalar ratio is r(k_0=0.002 Mpc^-1)<0.90 (95% CL).

We confront predictions of inflationary scenarios with the WMAP data, in combination with complementary small-scale CMB measurements and large-scale structure data. The WMAP detection of a large-angle anti-correlation in the temperature--polarization cross-power spectrum is the signature of adiabatic superhorizon fluctuations at the time of decoupling. The WMAP data are described by pure adiabatic fluctuations: we place an upper limit on a correlated CDM isocurvature component. Using WMAP constraints on the shape of the scalar power spectrum and the amplitude of gravity waves, we explore the parameter space of inflationary models that is consistent with the data. We place limits on inflationary models; for example, a minimally-coupled lambda phi^4 is disfavored at more than 3-sigma using WMAP data in combination with smaller scale CMB and large scale structure survey data. The limits on the primordial parameters using WMAP data alone are: n_s(k_0=0.002 Mpc^{-1})=1.20_{-0.11}^{+0.12}, dn/dlnk=-0.077^{+0.050}_{- 0.052}, A(k_0=0.002 Mpc}^{-1})=0.71^{+0.10}_{-0.11} (68% CL), and r(k_0=0.002 Mpc^{-1})<1.28 (95% CL).

We place tight constraints on the redshift-averaged, effective value of the equation of state of dark energy, w, using only the absolute ages of Galactic stars and the observed position of the first peak in the angular power spectrum of the CMB. We find w<-0.8 at the 68% confidence level. If we further consider that w > -1, this finding suggests that within our uncertainties, dark energy is indistinguishable from a classical vacuum energy term. We detect a correlation between the ages of the oldest galaxies and their redshift. This opens up the possibility of measuring w(z) by computing the relative ages of the oldest galaxies in the universe as a function of redshift, dz/dt. We show that this is a realistic possibility by computing dz/dt at z~0 from SDSS galaxies and obtain an independent estimate for the Hubble constant, H_0 = 69 \pm 12 km s-1 Mpc-1. The small number of galaxies considered at z>0.2 does not yield, currently, a precise determination of w(z), but shows that the age--redshift relation is consistent with a Standard LCDM universe with $w=-1$.

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking ~700 known bright sources from the maps, we estimate residual sources contribute ~3500 uK^2 at 41 GHz, and ~130 uK^2 at 94 GHz, to the power spectrum l*(l+1)*C_l/(2*pi) at l=1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to l~350. The spectrum clearly exhibits a first acoustic peak at l=220 and a second acoustic peak at l~540 and it provides strong support for adiabatic initial conditions. Kogut et al. (2003) analyze the C_l^TE power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by \~30% on degree angular scales, due to secondary scattering.

By studying the bispectrum of the galaxy distribution in the 2dF galaxy redhsift survey (2dFGRS), we have shown that 2dFGRS galaxies are unbiased tracers of the mass distribution. This allows us the break the degeneracy intrinsic to power spectrum studies, between the matter density parameter Omega_m and the bias parameter b, and to obtain an accurate measurement of Omega_m: Omega_{m,z_{eff}}=0.27 \pm 0.06, a measurement obtained from the 2dFGRS alone, independently from other data sets. This result has to be interpreted as Omega_m at the effective redshift of the survey z_{eff}=0.17. Extrapolated at z=0 we obtain Omega_m=0.2 \pm 0.06 for the LCDM model. This constraint on the matter density parameter when combined with cosmic microwave background constraints on the flatness of the Universe, show firm evidence that the Universe is dominated by a vacuum energy component.

The annihilation rate of weakly interacting cold dark matter particles at the galactic center could be greatly enhanced by the growth of a density spike around the central supermassive black hole (SBH). Here we discuss the effects of hierarchical mergers on the central spike. Mergers between halos containing SBHs lead to the formation of SBH binaries which transfer energy to the dark matter particles, lowering their density. The predicted flux of annihiliation photons from the galactic center is several orders of magnitude smaller than in models that ignore the effects of SBHs and mergers. Measurement of the annihilation radiation could in principle be used to constrain the merger history of the galaxy.

We compute the bispectrum of the 2dF Galaxy Redshift Survey (2dFGRS) and use it to measure the bias parameter of the galaxies. This parameter quantifies the strength of clustering of the galaxies relative to the mass in the Universe. By analysing 80 million triangle configurations in the wavenumber range 0.1 < k < 0.5 h/Mpc (i.e. on scales roughly between 5 and 30 Mpc/h) we find that the linear bias parameter is consistent with unity: b_1=1.04 pm 0.11, and the quadratic (nonlinear) bias is consistent with zero: b_2=-0.054 pm 0.08. Thus, at least on large scales, optically-selected galaxies do indeed trace the underlying mass distribution. The bias parameter can be combined with the 2dFGRS measurement of the redshift distortion parameter beta = Omega_m^{0.6}/b_1, to yield Omega_m = 0.27 pm 0.06 for the matter density of the Universe, a result which is determined entirely from this survey, independently of other datasets. Our measurement of the matter density of the Universe should be interpreted as Omega_m at the effective redshift of the survey (z=0.17).

We show that gas in a large fraction of low mass dark matter halos may form Toomre stable disks, if angular momentum is conserved when the gas contracts. Such halos would be stable to star formation and therefore remain dark. This may potentially explain the discrepancy between the predicted and observed number of dwarf satellites in the Local Group, as well as the deviation between the predicted and the observed faint end slope of the luminosity function. The above mechanism does not require a strong variation of the baryon fraction with the virial mass of the dark halo. We show that model fits to rotation curves are also consistent with this hypothesis: none of the observed galaxies lie in the region of parameter space forbidden by the Toomre stability criterion.

We compute the cosmic microwave background bispectrum arising from the cross-correlation of primordial, lensing and Rees-Sciama signals. The amplitude of the bispectrum signal is sensitive to the matter density parameter, Omega_0, and the equation of state of the dark energy, which we parameterize by w_Q. We conclude that the dataset of the Atacama Cosmology Telescope, combined with MAP 2-year data or the Planck data set alone will allow us to break the degeneracy between Omega_0 and w_Q that arises from the analysis of CMB power spectrum. In particular a joint measurement of Omega_0 and w_Q with 10% and 30% error on the two parameters respectively, at the 90% confidence level can realistically be achieved.

We study a large set of high spatial resolution optical rotation curves of galaxies with the goal of determining the model parameters for a disk embedded within a CDM halo that we model either with a NFW profile or pseudo-isothermal profile. We show that parameter degeneracies present in lower resolution data are lifted at these higher resolutions. 34% of the galaxies do not have a meaningful fit when using the NFW profile and 32% when using the pseudoisothermal profile, however only 14% do not have a meaningful fit in either model. In both models we find correlations between the disk baryon fraction $f_d$ and the spin parameter of the halo $\lambda'$, between $f_d$ and the dark halo mass $M_{200}$, and between $M_{200}$ and the concentration parameter $c$. We show that the distribution of the concentration parameters $c$, for a NFW halo, is in good agreement with CDM predictions; no significant galaxy population is found with very low values of $c$. The overall distribution of $\lambda'$ is in good agreement with theoretical predictions from hierarchical tidal torque theory. The whole sample is also well fitted by a pseudo-isothermal dark halo with a core, but the size of the core is rather small (for 70% of the sample the core size is less than 2 kpc). Thus we conclude that the profile of dark matter is steep ($r^{-1}$ or steeper) down to this radius; large dark matter cores (and therefore very low dark matter central densities) seem to be excluded. LSBs tend to have higher values of $\lambda'$ for a given $f_d$ and lower values of $c$ for a given mass than HSBs. In an appendix we give some useful formula for pseudo-isothermal profile halos and discuss in detail the issue of parameter degeneracies.

We show that the distribution of the sizes and temperatures of clusters can be used to constrain cosmological models. The size-temperature (ST) distribution predicted in a flat Gaussian cluster-abundance-normalized Omega_0=0.3 model agrees well with the fairly tight ST relation observed. A larger power-spectrum amplitude sigma_8 would give rise to a larger scatter about the ST relation as would a larger value of Omega_0 and/or long non-Gaussian high-density tails in the probability density function. For Gaussian initial conditions, the ST distribution suggests a constraint sigma_8 Omega_0^{0.26} \simeq 0.76. The ST relation is expected to get tighter at high redshifts. In the process, we derive a simple formula for the halo formation-redshift distribution for non-Gaussian models. We also suggest that the discrepancy between the naive zero-redshift ST relation and that observed may be due, at least in part, to the fact that lower-mass clusters form over a wider range of redshifts. An Appendix derives an equation for the formation-redshift distribution of halos.

An extensive sample of galaxy clusters will be available in the coming years, detected through their Sunyaev-Zeldovich effect (SZE). We use a semi-analytic model to study the scientific yield of combining SZE data with X-ray and optical follow-up observations. If clusters at a given redshift z_o can be identified with virialized, spherical halos, they populate a well-defined "fundamental plane" (FP) in the parameter space of the three observables: virial temperature (T), total Sunyaev-Zeldovich flux decrement (DS), and angular size (theta). The location and orientation of the FP, as well as its redshift-evolution, are sensitive to both the internal evolution of clusters, and to the underlying cosmological parameters. We show that if clusters are not standard candles (e.g. due to feedback, or energy injection), then this can be inferred from the FP. Likewise, we study the dependence of the FP on the cosmological parameters h, sigma_8, and Omega_0, and quantify constraints on these parameters. We also show that in the absence of any non-gravitational effects, the scatter in the the (DS-T) plane is significantly smaller than in either the (theta-T) or the (theta-DS) planes. As a result, the (DS-T) relation can be an exceptionally sensitive probe of both cluster physics and cosmological parameters, and a comparison of the amount of scatter in these three scaling relations will serve as a test of the origin (cosmological vs. stochastic) of the scatter.

We investigate the relative sensitivities of several tests for deviations from Gaussianity in the primordial distribution of density perturbations. We consider models for non-Gaussianity that mimic that which comes from inflation as well as that which comes from topological defects. The tests we consider involve the cosmic microwave background (CMB), large-scale structure (LSS), high-redshift galaxies, and the abundances and properties of clusters. We find that the CMB is superior at finding non-Gaussianity in the primordial gravitational potential (as inflation would produce), while observations of high-redshift galaxies are much better suited to find non-Gaussianity that resembles that expected from topological defects. We derive a simple expression that relates the abundance of high-redshift objects in non-Gaussian models to the primordial skewness.

In the standard model for structure formation, bound objects originate from the gravitational collapse of small perturbations arising from quantum fluctuations with random phases. In other scenarios, based on defects, structures are seeded by localized energy density. In principle, it is possible to differentiate between these models on the basis of their statistical properties; only in the former case is the initial density field an almost-perfect random gaussian field. In this paper, we investigate the use of the trispectrum of the galaxy density field, which is the connected four-point function in Fourier space, as a discriminant between gaussian and non-gaussian models. It has the advantage of having only weak non-linear growth. We define a related statistic $\tau$ which, as a test of the gaussian hypothesis, is independent of cosmology, the power spectrum and biasing, in real space, and which is, in principle, a measure of the departure from gaussian statistics. For galaxy redshift surveys, the statistic depends on cosmology and bias only through the potentially observable parameter $\beta$. We compute the expected errors on the estimate of $\tau$, and demonstrate with numerical simulations that it can be a useful discriminant of models, with the important proviso that any bias is linear on large scales. Whether it is the most effective method is uncertain and depends on the nature of the departure from gaussianity.

We present a theoretical and exact analysis of the bispectrum of projected galaxy catalogues. The result can be generalized to evaluate the projection in spherical harmonics of any 3D bispectrum and therefore has applications to cosmic microwave background and gravitational lensing studies. By expanding the 2D distribution of galaxies on the sky in spherical harmonics, we show how the 3-point function of the coefficients can be used in principle to determine the bias parameter of the galaxy sample. If this can be achieved, it would allow a lifting of the degeneracy between the bias and the matter density parameter of the Universe which occurs in linear analysis of 3D galaxy catalogues. In previous papers we have shown how a similar analysis can be done in three dimensions, and we show here, through an error analysis and by implementing the method on a simulated projected catalogue, that ongoing three-dimensional galaxy redshift surveys (even with all the additional uncertainties introduced by partial sky coverage, redshift-space distortions and smaller numbers) will do far better than all-sky projected catalogues with similar selection function.

Following the discovery of the CMB, the hot big-bang model has become the standard cosmological model. In this theory, small primordial fluctuations are subsequently amplified by gravity to form the large-scale structure seen today. Different theories for unified models of particle physics, lead to different predictions for the statistical properties of the primordial fluctuations, that can be divided in two classes: gaussian and non-gaussian. Convincing evidence against or for gaussian initial conditions would rule out many scenarios and point us towards a physical theory for the origin of structures. The statistical distribution of cosmological perturbations, as we observe them, can deviate from the gaussian distribution in several different ways. Even if perturbations start off gaussian, non-linear gravitational evolution can introduce non-gaussian features. Additionally, our knowledge of the Universe comes principally from the study of luminous material such as galaxies, but these might not be faithful tracers of the underlying mass distribution. The relationship between fluctuations in the mass and in the galaxies distribution (bias), is often assumed to be local, but could well be non-linear. Moreover, galaxy catalogues use the redshift as third spatial coordinate: the resulting redshift-space map is non-linearly distorted by peculiar velocities. Non-linear gravitational evolution, biasing, and redshift-space distortion introduce non-gaussianity, even in an initially gaussian fluctuation field. I will investigate the statistical tools that allow us, in principle, to disentangle the above different effects, and the observational datasets we require to do so in practice.

The observed abundance of high-redshift galaxies and clusters contains precious information about the properties of the initial perturbations. We present a method to compute analytically the number density of objects as a function of mass and redshift for a range of physically motivated non-Gaussian models. In these models the non-Gaussianity can be dialed from zero and is assumed to be small. We compute the probability density function for the smoothed dark matter density field and we extend the Press and Schechter approach to mildly non-Gaussian density fields. The abundance of high-redshift objects can be directly related to the non-Gaussianity parameter and thus to the physical processes that generated deviations from the Gaussian behaviour. Even a skewness parameter of order 0.1 implies a dramatic change in the predicted abundance of $z\gap 1$ objects. Observations from NGST and X-ray satellites (XMM) can be used to accurately measure the amount of non-Gaussianity in the primordial density field.

Since cosmic-microwave-background (CMB) and large-scale-structure (LSS) data will shortly improve dramatically with the Microwave Anisotropy Probe (MAP) and Planck Surveyor, and the Anglo-Australian 2-Degree Field (2dF) and Sloan Digital Sky Survey (SDSS), respectively, it is timely to ask which of the CMB or LSS will provide a better probe of primordial non-gaussianity. In this paper we consider this question, using the bispectrum as a discriminating statistic. We consider several non-gaussian models and find that in each case the CMB will provide a better probe of primordial non-gaussianity. Since the bispectrum is the lowest-order statistic expected to arise in a generic non-gaussian model, our results suggest that if CMB maps appear gaussian, then apparent deviations from gaussian initial conditions in galaxy surveys can be attributed with confidence to the effects of biasing. We demonstrate this precisely for the spatial bispectrum induced by local non-linear biasing.

We study the power spectrum of galaxies in redshift space, with third order perturbation theory to include corrections that are absent in linear theory. We assume a local bias for the galaxies: i.e. the galaxy density is sampled from some local function of the underlying mass distribution. We find that the effect of the nonlinear bias in real space is to introduce two new features: first, there is a contribution to the power which is constant with wavenumber, whose nature we reveal as essentially a shot-noise term. In principle this contribution can mask the primordial power spectrum, and could limit the accuracy with which the latter might be measured on very large scales. Secondly, the effect of second- and third-order bias is to modify the effective bias (defined as the square root of the ratio of galaxy power spectrum to matter power spectrum). The effective bias is almost scale-independent over a wide range of scales. These general conclusions also hold in redshift space. In addition, we have investigated the distortion of the power spectrum by peculiar velocities, which may be used to constrain the density of the Universe. We look at the quadrupole-to-monopole ratio, and find that higher-order terms can mimic linear theory bias, but the bias implied is neither the linear bias, nor the effective bias referred to above. We test the theory with biased N-body simulations, and find excellent agreement in both real and redshift space, providing the local biasing is applied on a scale whose fractional r.m.s. density fluctuations are $< 0.5$.

The determination of the density parameter $\Omega_0$ from the large-scale distribution of galaxies is one of the major goals of modern cosmology. However, if galaxies are biased tracers of the underlying mass distribution, linear perturbation theory leads to a degeneracy between $\Omega_0$ and the linear bias parameter $b$, and the density parameter cannot be estimated. In Matarrese, Verde & Heavens (1997) we developed a method based on second-order perturbation theory to use the bispectrum to lift this degeneracy by measuring the bias parameter in an $\Omega_0$-independent way. The formalism was developed assuming that one has perfect information on the positions of galaxies in three dimensions. In galaxy redshift surveys, the three-dimensional information is imperfect, because of the contaminating effects of peculiar velocities, and the resulting clustering pattern in redshift space is distorted. In this paper, we combine second-order perturbation theory with a model for collapsed, virialised structures, to extend the method to redshift space, and demonstrate that the method should be successful in determining with reasonable accuracy the bias parameter from state-of-the-art surveys such as the Anglo-Australian 2 degree field survey and the Sloan digital sky survey.