Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

We study the displacements between the centres of galaxies and their supermassive black holes (BHs) in the cosmological hydrodynamical simulation Horizon-AGN, and in a variety of observations from the literature. The BHs in Horizon-AGN feel a sub-grid dynamical friction force, sourced by the surrounding gas, which prevents recoiling BHs being ejected from the galaxy. We find that i) the fraction of spatially offset BHs increases with cosmic time, ii) BHs live on prograde orbits in the plane of the galaxy with an orbital radius that decays with time but stalls near $z=0$, and iii) the magnitudes of offsets from the galaxy centres are substantially larger in the simulation than in observations. We attribute the stalling of the infall and excessive offset magnitudes to the fact that dynamical friction from stars and dark matter is not modelled in the simulation, and hence provide a way to improve the black hole dynamics of future simulations.

We construct a model for the angular power spectrum of the instrumental noise in interferometer networks mapping gravitational wave backgrounds (GWBs) as a function of detector noise properties, network configuration and scan strategy. We use the model to calculate the noise power spectrum for current and future ground-based experiments, as well as for planned space missions. We present our results in a language similar to that used in cosmic microwave background and intensity mapping experiments, and connect the formalism with the sensitivity curves that are common lore in GWB analyses. Our formalism is implemented in a lightweight python module that we make publicly available at https://github.com/damonge/schNell.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

The decay timescales of the quasinormal modes of a massive scalar field have an intriguing behavior: they either grow or decay with increasing angular harmonic numbers $\ell$, depending on whether the mass of the scalar field is small or large. We identify the properties of the effective potential of the scalar field that leads to this behavior and characterize it in detail. If the scalar field is non-minimally coupled, considered here, the scalar quasinormal modes will leak into the gravitational wave signal and will have decaying times that are comparable or smaller than those typical in General Relativity. Hence, these modes could be detectable in the future. Finally, we find that the anomalous behavior in the decay timescales of quasinormal modes is present in a much larger class of models beyond a simple massive scalar field.

The late time acceleration of the Universe can be characterized in terms of an extra, time dependent, component of the universe -- dark energy. The simplest proposal for dark energy is a scalar-tensor theory -- quintessence -- which consists of a scalar field, $\phi$, whose dynamics is solely dictated by its potential, $V(\phi)$. Such a theory can be uniquely characterized by the equation of state of the scalar field energy momentum-tensor. We find the time dependence of the equation of state for a broad family of potentials and, using this information, we propose an analytic prior distribution for the most commonly used parametrization. We show that this analytic prior can be used to accurately predict the distribution of observables for the next generation of cosmological surveys. Including the theoretical priors in the comparison with observations considerably improves the constraints on the equation of state.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

The spatial and temporal discreteness of gravitational wave sources leads to shot noise that may, in some regimes, swamp any attempts at measuring the anisotropy of the gravitational wave background. Cross-correlating a gravitational wave background map with a sufficiently dense galaxy survey can alleviate this issue, and potentially recover some of the underlying properties of the gravitational wave background. We quantify the shot noise level and we explicitly show that cross-correlating the gravitational wave background and a galaxy catalog improves the chances of a first detection of the background anisotropy with a gravitational wave observatory operating in the frequency range (10Hz,100Hz), given sufficient sensitivity.

Self-interacting dark matter (SIDM) is an able alternative to collisionless dark matter. If dark matter does have self-interactions, we would expect this to cause a separation between the collisionless stars and the dark matter halo of a galaxy as it falls through a dark matter medium. For stars arranged in a disk, this would generate a U-shaped warp. The magnitude of this warping depends on the SIDM cross section, type of self-interaction, relative velocity of galaxy and background, halo structure, and density of the dark matter medium. In this paper, we set constraints on long-range (light mediator) dark matter self-interaction by means of this signal. We begin by measuring U-shaped warps in $3,213$ edge-on disk galaxies within the Sloan Digital Sky Survey. We then forward-model the expected warp from SIDM on a galaxy-by-galaxy basis by combining models of halo structure, density and velocity field reconstructions, and models for the dark matter interactions. We find no evidence for a contribution to the warps from SIDM. Our constraints are highly dependent on the uncertain velocities of our galaxies: for a normalized Rutherford-like cross section we find $\tilde{\sigma}/m_{\rm{DM}} \lesssim 3\times 10^{-13}~\rm{cm}^2/\rm{g}$ at fixed velocity $v = 300~\rm{km/s}$ -- a bound that scales roughly linearly with increasing $v$. In the appendix we translate these bounds into limits on the momentum transfer cross section, finding $\sigma_T(300~\rm{km/s})/m_{\rm{DM}} \lesssim 0.1~\rm{cm}^2/\rm{g}$. [abridged]

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

Scale invariant theories of gravity give a compelling explanation to the early and late time acceleration of the Universe. Unlike most scalar-tensor theories, fifth forces are absent and it would therefore seem impossible to distinguish scale invariant gravity from general relativity. We show that the ringdown of a Schwarschild-de Sitter black hole may have a set of massive modes which are characteristic of scale invariant gravity. In principle these new modes can be used to distinguish scale invariant gravity from general relativity. In practice, we discuss the obstacles to generating these new massive modes and their detectability with future gravitational wave experiments but also speculate on their role in Kerr black holes.

Through numerical simulations of a minimally coupled massive Klein-Gordon scalar field, we show that it is possible to grow hair on a Schwarzschild black hole if one assumes an initial periodically time-varying but spatially homogeneous scalar background. By "hair", we mean a non-trivial profile in the scalar field. We find that this profile emerges on a timescale related to the mass of the black hole, with features related to the mass of the scalar particle. We undertake simulations with and without backreaction on the metric and see that the essential, qualitative features remain consistent. We also contrast the results from higher mass scalars to the case of a low mass with a large Compton wavelength. The results are particularly relevant for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g. axions.

One of the main science goals of the Large Synoptic Survey Telescope (LSST) is to uncover the nature of cosmic acceleration. In the base analysis, possible deviations from the Lambda-Cold-Dark-Matter ($\Lambda$CDM) background evolution will be probed by fitting a $w(z)$CDM model, which allows for a redshift-dependent dark energy equation of state with $w(z)$, within general relativity (GR). A rich array of other phenomena can arise due to deviations from the standard $\Lambda$CDM+GR model though, including modifications to the growth rate of structure and lensing, and novel screening effects on non-linear scales. Concrete physical models are needed to provide consistent predictions for these (potentially small) effects, to give us the best chance of detecting them and separating them from astrophysical systematics. A complex plethora of possible models has been constructed over the past few decades, with none emerging as a particular favorite. This document prioritizes a subset of these models along with rationales for further study and inclusion into the LSST Dark Energy Science Collaboration (DESC) data analysis pipelines, based on their observational viability, theoretical plausibility, and level of theoretical development. We provide references and theoretical expressions to aid the integration of these models into DESC software and simulations, and give justifications for why other models were not prioritized. While DESC efforts are free to pursue other models, we provide here guidelines on which theories appear to have higher priority for collaboration efforts due to their perceived promise and greater instructional value.

We study the phenomenology of the beyond Horndeski class of scalar-tensor theories of gravity, which on cosmological scales can be characterised in terms of one extra function of time, $\alpha_{\rm H}$, as well as the usual four Horndeski set of free functions. We show that $\alpha_{\rm H}$ can be directly related to the the damping of the matter power spectrum on both large and small scales. We also find that the temperature power spectrum of the cosmic microwave background (CMB) is enhanced at low multipoles and the lensing potential is decreased, as a function of $\alpha_{\rm H}$. We find constraints on $\alpha_{\rm H}$ of order ${\cal O}(1)$ using measurements of the temperature and polarisation of the CMB, as well as the lensing potential derived from it, combined with large scale structure data. We find that redshift space distortion measurements can play a significant role in constraining these theories. Finally, we comment on the recent constraints from the observation of an electromagnetic counterpart to a gravitational wave signal; we find that these constraints reduce the number of free parameters of the model but do not significantly change the constraints on the remaining parameters.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

Weyl (scale) invariant theories of scalars and gravity can generate all mass scales spontaneously. In this paper we study a particularly simple version -- scale invariant $R^2$ gravity -- and show that, during an inflationary period, it leads to fluctuations which, for a particular parameter choice, are almost indistinguishable from normal $R^2$ inflation. Current observations place tight constraints on the primary coupling constant of this theory and predict a tensor to scalar ratio, $0.0033 > r > 0.0026$, which is testable with the next generation of cosmic microwave background experiments.

Observational cosmology in the next decade will rely on probes of the distribution of matter in the redshift range between $0<z<3$ to elucidate the nature of dark matter and dark energy. In this redshift range, galaxy formation is known to have a significant impact on observables such as two-point correlations of galaxy shapes and positions, altering their amplitude and scale dependence beyond the expected statistical uncertainty of upcoming experiments at separations under 10 Mpc. Successful extraction of information in such a regime thus requires, at the very least, unbiased models for the impact of galaxy formation on the matter distribution, and can benefit from complementary observational priors. This work reviews the current state of the art in the modelling of baryons for cosmology, from numerical methods to approximate analytical prescriptions, and makes recommendations for studies in the next decade, including a discussion of potential probe combinations that can help constrain the role of baryons in cosmological studies. We focus, in particular, on the modelling of the matter power spectrum, $P(k,z)$, as a function of scale and redshift, and of the observables derived from this quantity. This work is the result of a workshop held at the University of Oxford in November of 2018.

We investigate the prospect of using black hole spectroscopy to constrain the parameters of Horndeski gravity through observations of gravitational waves from perturbed black holes. We study the gravitational waves emitted during ringdown from black holes without hair in Horndeski gravity, demonstrating the qualitative differences between such emission in General Relativity and Horndeski theory. In particular, Quasi-Normal Mode frequencies associated with the scalar field spectrum can appear in the emitted gravitational radiation. Analytic expressions for error estimates for both the black hole and Horndeski parameters are calculated using a Fisher Matrix approach, with constraints on the `effective mass' of the Horndeski scalar field of order $\sim 10^{-17}$eV$c^{-2}$ or tighter being shown to be achievable in some scenarios. Estimates for the minimum signal-noise-ratio required to observe such a signal are also presented.

Cosmological observations are beginning to reach a level of precision that allow us to test some of the most fundamental assumptions in our working model of the Universe. One such an assumption is that gravity is governed by the General Theory of Relativity. In this review we discuss how one might go about extending General Relativity and how such extensions can be described in a unified way on large scales. This allows us to describe the impact of modified gravity on the growth and morphology of the large scale structure of the Universe. On smaller scales we explore the physics of gravitational screening and how it might manifest itself in galaxies, clusters and, more generally, in the cosmic web. We then analyze the current constraints from large scale structure and conclude by discussing the future prospects of the field in light of the plethora of surveys currently being planned.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

Cosmological large-scale structure analyses based on two-point correlation functions often assume a Gaussian likelihood function with a fixed covariance matrix. We study the impact on cosmological parameter estimation of ignoring the parameter dependence of this covariance matrix, focusing on the particular case of joint weak-lensing and galaxy clustering analyses. Using a Fisher matrix formalism (calibrated against exact likelihood evaluation in particular simple cases), we quantify the effect of using a parameter dependent covariance matrix on both the bias and variance of the parameters. We confirm that the approximation of a parameter-independent covariance matrix is exceptionally good in all realistic scenarios. The information content in the covariance matrix (in comparison with the two point functions themselves) does not change with the fractional sky coverage. Therefore the increase in information due to the parameter dependent covariance matrix becomes negligible as the number of modes increases. Even for surveys covering less than $1\%$ of the sky, this effect only causes a bias of up to ${\cal O}(10\%)$ of the statistical uncertainties, with a misestimation of the parameter uncertainties at the same level or lower. The effect will only be smaller with future large-area surveys. Thus for most analyses the effect of a parameter-dependent covariance matrix can be ignored both in terms of the accuracy and precision of the recovered cosmological constraints.

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

Extensions of the standard models of particle physics and cosmology often lead to long-range fifth forces with properties dependent on gravitational environment. Fifth forces on astrophysical scales are best studied in the cosmic web where perturbation theory breaks down. We present constraints on chameleon- and symmetron-screened fifth forces with Yukawa coupling and megaparsec range -- as well as unscreened fifth forces with differential coupling to galactic mass components -- by searching for the displacements they predict between galaxies' stars and gas. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we set upper bounds on fifth-force strength relative to Newtonian gravity from $\Delta G/G_N < \text{few} \: \times 10^{-4}$ for range $\lambda_C = 50$ Mpc, to $\Delta G/G_N \lesssim 0.1$ for $\lambda_C = 500$ kpc. In $f(R)$ gravity this requires $f_{R0} < \text{few} \: \times \: 10^{-8}$. The analogous bounds without screening are $\Delta G/G_N < \text{few} \: \times 10^{-4}$ and $\Delta G/G_N < \text{few} \times 10^{-3}$. These are the tightest and among the only fifth-force constraints on galaxy scales. We show how our results may be strengthened with future survey data and identify the key features of an observational programme for furthering fifth-force tests beyond the Solar System.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

Intra-galaxy signals contain a wealth of information on fundamental physics, both the dark sector and the nature of gravity. While so far largely unexplored, such probes are set to rise dramatically in importance as upcoming surveys provide data of unprecedented quantity and quality on galaxy structure and dynamics. In this paper, we use warping of stellar disks to test the chameleon- or symmetron-screened fifth forces which generically arise when new fields couple to matter. We take r-band images of mostly late-type galaxies from the Nasa Sloan Atlas and develop an automated algorithm to quantify the degree of U-shaped warping they exhibit. We then forward-model the warp signal as a function of fifth-force strength $\Delta G/G_N$ and range $\lambda_C$, and the gravitational environments and internal properties of the galaxies, including full propagation of the non-Gaussian uncertainties. Convolving this fifth-force likelihood function with a Gaussian describing astrophysical and observational noise and then constraining $\Delta G/G_N$ and $\lambda_C$ by Markov Chain Monte Carlo, we find the overall likelihood to be significant increased ($\Delta\log(\mathcal{L}) \simeq 20$) by adding a screened fifth force with $\lambda_C \simeq 2$ Mpc, $\Delta G/G_N \simeq 0.01$. The variation of $\Delta\log(\mathcal{L})$ with $\lambda_C$ is quantitatively as expected from the correlation of the magnitude of the fifth-force field with the force's range, and a similar model without screening achieves no increase in likelihood over the General Relativistic case $\Delta G=0$. Although these results are in good agreement with a previous analysis of the same model using offsets between galaxies' stellar and gas mass centroids (Desmond et al. 2018), we caution that the effects of confounding baryonic and dark matter physics must be thoroughly investigated for the results of the inference to be unambiguous.

The geometric optics approximation traditionally used to study the propagation of gravitational waves on a curved background, breaks down in the vicinity of compact and extended astrophysical objects, where wave-like effects like diffusion and generation of polarization occur. We provide a framework to study the generation of polarization of a stochastic background of gravitational waves propagating in an inhomogeneous universe. The framework is general and can be applied to both cosmological and astrophysical gravitational wave backgrounds in any frequency range. We derive an order of magnitude estimate of the amount of polarization generated for cosmological and astrophysical backgrounds, in the frequency range covered by present and planned gravitational wave experiments. For an astrophysical background in the PTA and LISA band, the amount of polarization generated is suppressed by a factor $10^{-4}$ $(10^{-5})$ with respect to anisotropies. For a cosmological background we get an additional $10^{-2}$ suppression. We speculate on using our approach to map the distribution of (unresolvable) structures in the Universe.

One of the most common consequences of extensions to the standard models of particle physics or cosmology is the emergence of a fifth force. While generic fifth forces are tightly constrained at Solar System scales and below, they may escape detection by means of a screening mechanism which effectively removes them in dense environments. We constrain the strength $\Delta G/G_N$ and range $\lambda_C$ of a chameleon- or symmetron-screened fifth force with Yukawa coupling -- as well as an unscreened fifth force with differential coupling to galactic mass components -- by searching for the displacement it predicts between galaxies' stellar and gas mass centroids. Taking data from the Alfalfa HI survey, identifying galaxies' gravitational environments with the maps of Desmond et al. (2018a) and forward-modelling with a Bayesian likelihood framework, we find $6.6\sigma$ evidence for $\Delta G>0$ at $\lambda_C \simeq 2$ Mpc, with $\Delta G/G_N = 0.025$ at maximum-likelihood. A similar fifth-force model without screening gives no increase in likelihood over the case $\Delta G = 0$ for any $\lambda_C$. Although we validate these results by several methods, we do not claim screened modified gravity to provide the only possible explanation for the signal: this would require knowing that "galaxy formation" physics could not be responsible. We show also the results of a more conservative -- though less well motivated -- noise model which yields only upper limits on $\Delta G/G_N$, ranging from $\sim10^{-1}$ for $\lambda_C \simeq 0.5$ Mpc to $\sim \: \text{few} \times 10^{-4}$ at $\lambda_C \simeq 50$ Mpc. We show how these constraints may be improved by future galaxy surveys and identify the key features of an observational programme for directly constraining fifth forces on galactic scales. This paper provides a complete description of the analysis summarised in Desmond et al. (2018b).

We consider the cosmological dynamics of a scalar field in a potential with multiple troughs and peaks. We show that the dynamics of the scalar field will evolve from light dark matter-like behaviour (such as that of a light axion) to a combination of heavy dark matter-like and dark energy-like behaviour. We discuss the phenomenology of such a model, explaining how it can give rise to the cosmological constant, as well as how it can decouple the dark sector densities between the time of recombination and today, for both the homogeneous background and perturbations. The final form of the dark matter is axion-like, but with abundance and primordial isocurvature modes taking very different values from traditional, axionic, dark matter.

A scale-invariant universe can have a period of accelerated expansion at early times: inflation. We use a frame-invariant approach to calculate inflationary observables in a scale invariant theory of gravity involving two scalar fields - the spectral indices, the tensor to scalar ratio, the level of isocurvature modes and non-Gaussianity. We show that scale symmetry leads to an exact cancellation of isocurvature modes and that, in the scale-symmetry broken phase, this theory is well described by a single scalar field theory. We find the predictions of this theory strongly compatible with current observations.

We study the perturbations to General Relativistic black holes (i.e. those without scalar hair) in Horndeski scalar-tensor gravity. First, we derive the equations of odd and even parity perturbations of both the metric and scalar field in the case of a Schwarzschild black hole, and show that the gravitational waves emitted from such a system contain a mixture of quasi-normal mode frequencies from the usual General Relativistic spectrum and those from the new scalar field spectrum, with the new scalar spectrum characterised by just two free parameters. We then specialise to the sub-family of Horndeski theories in which gravitational waves propagate at the speed of light $c$ on cosmological backgrounds; the scalar quasi-normal mode spectrum of such theories is characterised by just a single parameter $\mu$ acting as an effective mass of the scalar field. Analytical expressions for the quasi-normal mode frequencies of the scalar spectrum in this sub-family of theories are provided for both static and slowly rotating black holes. In both regimes comparisons to quasi-normal modes calculated numerically show good agreement with those calculated analytically in this work.

A new concept for the direct measurement of muons in air showers is presented. The concept is based on resistive plate chambers (RPCs), which can directly measure muons with very good space and time resolution. The muon detector is shielded by placing it under another detector able to absorb and measure the electromagnetic component of the showers such as a water-Cherenkov detector, commonly used in air shower arrays. The combination of the two detectors in a single, compact detector unit provides a unique measurement that opens rich possibilities in the study of air showers.

The recent detection of GRB 170817A and GW170817 constrains the speed of gravity waves $c_T$ to be that of light, which severely restricts the landscape of modified gravity theories that impact the cosmological evolution of the universe. In this work, we investigate the presence of black hole hair in the remaining viable cosmological theories of modified gravity that respect the constraint $c_T=1$. We focus mainly on scalar-tensor theories of gravity, analyzing static, asymptotically flat black holes in Horndeski, Beyond Horndeski, Einstein-Scalar-Gauss-Bonnet, and Chern-Simons theories. We find that in all of the cases considered here, theories that respect $c_T=1$ do not allow for hair, or have negligible hair. We further comment on vector-tensor theories including Einstein Yang-Mills, Einstein-Aether, and Generalised Proca theories, as well as bimetric theories.

We report on the measurement of the $^{7}$Be($n, p$)$^{7}$Li cross section from thermal to approximately 325 keV neutron energy, performed in the high-flux experimental area (EAR2) of the n_TOF facility at CERN. This reaction plays a key role in the lithium yield of the Big Bang Nucleosynthesis (BBN) for standard cosmology. The only two previous time-of-flight measurements performed on this reaction did not cover the energy window of interest for BBN, and showed a large discrepancy between each other. The measurement was performed with a Si-telescope, and a high-purity sample produced by implantation of a $^{7}$Be ion beam at the ISOLDE facility at CERN. While a significantly higher cross section is found at low-energy, relative to current evaluations, in the region of BBN interest the present results are consistent with the values inferred from the time-reversal $^{7}$Li($p, n$)$^{7}$Be reaction, thus yielding only a relatively minor improvement on the so-called Cosmological Lithium Problem (CLiP). The relevance of these results on the near-threshold neutron production in the p+$^{7}$Li reaction is also discussed.

We use the covariant formulation proposed in Tattersall et al (2017) to analyse the structure of linear perturbations about a spherically symmetric background in different families of gravity theories, and hence study how quasi-normal modes of perturbed black holes may be affected by modifications to General Relativity. We restrict ourselves to single-tensor, scalar-tensor and vector-tensor diffeomorphism-invariant gravity models in a Schwarzschild black hole background. We show explicitly the full covariant form of the quadratic actions in such cases, which allow us to then analyse odd parity (axial) and even parity (polar) perturbations simultaneously in a straightforward manner.

In this paper are discussed possible many body generalizations of the expanding locally anisotropic metric ansatz with respect to approximately Newtonian gravitational systems. This ansatz consistently describes local point-like matter distributions on the expanding Universe also allowing for a covariant parameterization of gravitational interactions at intermediate length scales. As an example of applicability it is modeled a disk galaxy model matching the physical parameters of the galaxy UGC2885 and it is shown that, by fine-tuning the metric functional parameter, the flattening of the galaxy rotation curve is fully parameterized by this metric. In addition it is numerically computed the mass-energy density corrections due to the expanding anisotropic background and explicitly shown that although there are negative contributions within the galaxy plane the total mass-energy density is strictly positive both at the galaxy plane and outside the galaxy plane. As the functional parameter for this metric is an exponential factor is required a floating point precision of $250$ significant digits for root finding routines and $200$ significant digits to evaluate the effective mass-energy density rendering a final precision of the results presented above double precision ($16$ significant digits). It is further shown that these results are consistent with the interpretation of the gravitational corrections as due to Dark Matter, in particular constituting a novel heuristic local parameterization for the Dark Matter distribution within the galaxy plane consistent with both local scale and cosmological scale physical laws which is useful to further investigate the local properties of Dark Matter.

Future surveys will access large volumes of space and hence very long wavelength fluctuations of the matter density and gravitational field. It has been argued that the set of secondary effects that affect the galaxy distribution, relativistic in nature, will bring new, complementary cosmological constraints. We study this claim in detail by focusing on a subset of wide-area future surveys: Stage-4 cosmic microwave background experiments and photometric redshift surveys. In particular, we look at the magnification lensing contribution to galaxy clustering and general relativistic corrections to all observables. We quantify the amount of information encoded in these effects in terms of the tightening of the final cosmological constraints as well as the potential bias in inferred parameters associated with neglecting them. We do so for a wide range of cosmological parameters, covering neutrino masses, standard dark-energy parametrizations and scalar-tensor gravity theories. Our results show that, while the effect of lensing magnification to number counts does not contain a significant amount of information when galaxy clustering is combined with cosmic shear measurements, this contribution does play a significant role in biasing estimates on a host of parameter families if unaccounted for. Since the amplitude of the magnification term is controlled by the slope of the source number counts with apparent magnitude, $s(z)$, we also estimate the accuracy to which this quantity must be known to avoid systematic parameter biases, finding that future surveys will need to determine $s(z)$ to the $\sim$5-10\% level. On the contrary, large-scale general-relativistic corrections are irrelevant both in terms of information content and parameter bias for most cosmological parameters, but significant for the level of primordial non-Gaussianity.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

We analyse cosmological perturbations around a homogeneous and isotropic background for scalar-tensor, vector-tensor and bimetric theories of gravity. Building on previous results, we propose a unified view of the effective parameters of all these theories. Based on this structure, we explore the viable space of parameters for each family of models by imposing the absence of ghosts and gradient instabilities. We then focus on the quasistatic regime and confirm that all these theories can be approximated by the phenomenological two-parameter model described by an effective Newton's constant and the gravitational slip. Within the quasistatic regime we pinpoint signatures which can distinguish between the broad classes of models (scalar-tensor, vector-tensor or bimetric). Finally, we present the equations of motion for our unified approach in such a way that they can be implemented in Einstein-Boltzmann solvers.

Tests of gravity at the galaxy scale are in their infancy. As a first step to systematically uncovering the gravitational significance of galaxies, we map three fundamental gravitational variables -- the Newtonian potential, acceleration and curvature -- over the galaxy environments of the local universe to a distance of approximately 200 Mpc. Our method combines the contributions from galaxies in an all-sky redshift survey, halos from an N-body simulation hosting low-luminosity objects, and linear and quasi-linear modes of the density field. We use the ranges of these variables to determine the extent to which galaxies expand the scope of generic tests of gravity and are capable of constraining specific classes of model for which they have special significance. Finally, we investigate the improvements afforded by upcoming galaxy surveys.

In this work are computed analytical solutions for orbital motion on a background described by an Expanding Locally Anisotropic (ELA) metric ansatz. This metric interpolates between the Schwarzschild metric near the central mass and the Robertson-Walker metric describing the expanding cosmological background far from the central mass allowing for a fine-tuneable covariant parameterization of gravitational interactions corrections in between these two asymptotic limits. In particular it is shown that the decrease of the Sun's mass by radiation emission plus the General Relativity corrections due to the ELA metric background with respect to Schwarzschild backgrounds can be mapped to the reported yearly variation of the gravitational constant $\dot{G}$ through Kepler's third law. Based on the value of the heuristic fit corresponding to the more recent heliocentric ephemerides of the Solar System are derived bounds for the value of a constant parameter $\alpha_0$ for the ELA metric as well as the maximal corrections to perihelion advance and orbital radii variation within this framework. Hence it is shown that employing the ELA metric as a functional covariant parameterization to model gravitational interactions corrections within the Solar System allows to maintain the measurement projection standards constant over time, specifically both the Astronomical Unit ($AU$) and the Gravitational constant ($G$). Also it is noted that the effect obtained is not homogeneous for all planetary orbits consistently with the diversity of estimates in the literature obtained assuming Schwarzschild backgrounds.

The detection of an electromagnetic counterpart (GRB 170817A) to the gravitational wave signal (GW170817) from the merger of two neutron stars opens a completely new arena for testing theories of gravity. We show that this measurement allows us to place stringent constraints on general scalar-tensor and vector-tensor theories, while allowing us to place an independent bound on the graviton mass in bimetric theories of gravity. These constraints severely reduce the viable range of cosmological models that have been proposed as alternatives to general relativistic cosmology.

We discuss the ground-breaking science that will be possible with a wide area survey, using the MeerKAT telescope, known as MeerKLASS (MeerKAT Large Area Synoptic Survey). The current specifications of MeerKAT make it a great fit for science applications that require large survey speeds but not necessarily high angular resolutions. In particular, for cosmology, a large survey over $\sim 4,000 \, {\rm deg}^2$ for $\sim 4,000$ hours will potentially provide the first ever measurements of the baryon acoustic oscillations using the 21cm intensity mapping technique, with enough accuracy to impose constraints on the nature of dark energy. The combination with multi-wavelength data will give unique additional information, such as exquisite constraints on primordial non-Gaussianity using the multi-tracer technique, as well as a better handle on foregrounds and systematics. Such a wide survey with MeerKAT is also a great match for HI galaxy studies, providing unrivalled statistics in the pre-SKA era for galaxies resolved in the HI emission line beyond local structures at z > 0.01. It will also produce a large continuum galaxy sample down to a depth of about 5\,$\mu$Jy in L-band, which is quite unique over such large areas and will allow studies of the large-scale structure of the Universe out to high redshifts, complementing the galaxy HI survey to form a transformational multi-wavelength approach to study galaxy dynamics and evolution. Finally, the same survey will supply unique information for a range of other science applications, including a large statistical investigation of galaxy clusters as well as produce a rotation measure map across a huge swathe of the sky. The MeerKLASS survey will be a crucial step on the road to using SKA1-MID for cosmological applications and other commensal surveys, as described in the top priority SKA key science projects (abridged).

We present a covariant formulation for constructing general quadratic actions for cosmological perturbations, invariant under a given set of gauge symmetries for a given field content. This approach allows us to analyse scalar, vector and tensor perturbations at the same time in a straightforward manner. We apply the procedure to diffeomorphism invariant single-tensor, scalar-tensor and vector-tensor theories and show explicitly the full covariant form of the quadratic actions in such cases, in addition to the actions determining the evolution of vector and tensor perturbations. We also discuss the role of the symmetry of the background in identifying the set of cosmologically relevant free parameters describing these classes of theories, including calculating the relevant free parameters for an axisymmetric Bianchi-I vacuum universe.

We consider a universe with an arbitrary number of extra dimensions, $N$. We present a new method for constructing the cosmological equations of motion and find analytic solutions with an explicit dependence on $N$. When we take the $N\rightarrow\infty$ limit we find novel, emergent behaviour which distinguishes it from normal Kaluza-Klein universes.

We present the public version of hi_class (www.hiclass-code.net), an extension of the Boltzmann code CLASS to a broad ensemble of modifications to general relativity. In particular, hi_class can calculate predictions for models based on Horndeski's theory, which is the most general scalar-tensor theory described by second-order equations of motion and encompasses any perfect-fluid dark energy, quintessence, Brans-Dicke, $f(R)$ and covariant Galileon models. hi_class has been thoroughly tested and can be readily used to understand the impact of alternative theories of gravity on linear structure formation as well as for cosmological parameter extraction.

Imaging surveys of galaxies will have a high number density and angular resolution yet a poor redshift precision. Intensity maps of neutral hydrogen (HI) will have accurate redshift resolution yet will not resolve individual sources. Using this complementarity, we show how the clustering redshifts approach, proposed for spectroscopic surveys can also be used in combination with intensity mapping observations to calibrate the redshift distribution of galaxies in an imaging survey and, as a result, reduce uncertainties in photometric redshift measurements. We show how the intensity mapping surveys to be carried out with the MeerKAT, HIRAX and SKA instruments can improve photometric redshift uncertainties to well below the requirements of DES and LSST. The effectiveness of this method as a function of instrumental parameters, foreground subtraction and other potential systematic errors is discussed in detail.

The next generation of surveys will greatly improve our knowledge of cosmological gravity. In this paper we focus on how Stage IV photometric redshift surveys, including weak lensing and multiple tracers of the matter distribution and radio experiments combined with measurements of the cosmic microwave background will lead to precision constraints on deviations from General Relativity. We use a broad subclass of Horndeski scalar-tensor theories to forecast the accuracy with which we will be able to determine these deviations and their degeneracies with other cosmological parameters. Our analysis includes relativistic effects, does not rely on the quasi-static evolution and makes conservative assumptions about the effect of screening on small scales. We define a figure of merit for cosmological tests of gravity and show how the combination of different types of surveys, probing different length scales and redshifts, can be used to pin down constraints on the gravitational physics to better than a few percent, roughly an order of magnitude better than present probes. Future cosmological experiments will be able to constrain the Brans-Dicke parameter at a level comparable to Solar System and astrophysical tests

We implement the method developed in [1] to construct the most general parametrised action for linear cosmological perturbations of bimetric theories of gravity. Specifically, we consider perturbations around a homogeneous and isotropic background, and identify the complete form of the action invariant under diffeomorphism transformations, as well as the number of free parameters characterising this cosmological class of theories. We discuss, in detail, the case without derivative interactions, and compare our results with those found in massive bigravity.

Scalar fields, $\phi_i$ can be coupled non-minimally to curvature and satisfy the general criteria: (i) the theory has no mass input parameters, including the Planck mass; (ii) the $\phi_i$ have arbitrary values and gradients, but undergo a general expansion and relaxation to constant values that satisfy a nontrivial constraint, $K(\phi_i) =$ constant; (iii) this constraint breaks scale symmetry spontaneously, and the Planck mass is dynamically generated; (iv) there can be adequate inflation associated with slow roll in a scale invariant potential subject to the constraint; (v) the final vacuum can have a small to vanishing cosmological constant (vi) large hierarchies in vacuum expectation values can naturally form; (vii) there is a harmless dilaton which naturally eludes the usual constraints on massless scalars. These models are governed by a global Weyl scale symmetry and its conserved current, $K_\mu$ . At the quantum level the Weyl scale symmetry can be maintained by an invariant specification of renormalized quantities.

We revisit the possibility that the Planck mass is spontaneously generated in scale invariant scalar-tensor theories of gravity, typically leading to a "dilaton." The fifth force, arising from the dilaton, is severely constrained by astrophysical measurements. We explore the possibility that nature is fundamentally Weyl-scale invariant and argue that, as a consequence, the fifth force effects are dramatically suppressed and such models are viable. We discuss possible obstructions to maintaining scale invariance and how these might be resolved.

We discuss models involving two scalar fields coupled to classical gravity that satisfy the general criteria: (i) the theory has no mass input parameters, (ii) classical scale symmetry is broken only through $-\frac{1}{12}\varsigma \phi^2 R$ couplings where $\varsigma$ departs from the special conformal value of $1$; (iii) the Planck mass is dynamically generated by the vacuum expectations values (VEVs) of the scalars (iv) there is a stage of viable inflation associated with slow roll in the two--scalar potential; (v) the final vacuum has a small to vanishing cosmological constant and an hierarchically small ratio of the VEVs and the ratio of the scalar masses to the Planck scale. This assumes the paradigm of classical scale symmetry as a custodial symmetry of large hierarchies.

We present a method for parametrizing linear cosmological perturbations of theories of gravity, around homogeneous and isotropic backgrounds. The method is sufficiently general and systematic that it can be applied to theories with any degrees of freedom (DoFs) and arbitrary gauge symmetries. In this paper, we focus on scalar-tensor and vector-tensor theories, invariant under linear coordinate transformations. In the case of scalar-tensor theories, we use our framework to recover the simple parametrizations of linearized Horndeski and "Beyond Horndeski" theories, and also find higher-derivative corrections. In the case of vector-tensor theories, we first construct the most general quadratic action for perturbations that leads to second-order equations of motion, which propagates two scalar DoFs. Then we specialize to the case in which the vector field is time-like (\`a la Einstein-Aether gravity), where the theory only propagates one scalar DoF. As a result, we identify the complete forms of the quadratic actions for perturbations, and the number of free parameters that need to be defined, to cosmologically characterize these two broad classes of theories.

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.

Future ground-based CMB experiments will generate competitive large-scale structure datasets by precisely characterizing CMB secondary anisotropies over a large fraction of the sky. We describe a method for constraining the growth rate of structure to sub-1% precision out to $z\approx 1$, using a combination of galaxy cluster peculiar velocities measured using the kinetic Sunyaev-Zel'dovich (kSZ) effect, and the velocity field reconstructed from galaxy redshift surveys. We consider only thermal SZ-selected cluster samples, which will consist of $\mathcal{O}(10^4-10^5)$ sources for Stage 3 and 4 CMB experiments respectively. Three different methods for separating the kSZ effect from the primary CMB are compared, including a novel blind "constrained realization" method that improves signal-to-noise by a factor of $\sim 2$ over a commonly-used aperture photometry technique. Measurements of the integrated tSZ $y$-parameter are used to break the kSZ velocity-optical depth degeneracy, and the effects of including CMB polarization and SZ profile uncertainties are also considered. A combination of future Stage 4 experiments should be able to measure the product of the growth and expansion rates, $\alpha\equiv f H$, to better than 1% in bins of $\Delta z = 0.1$ out to $z \approx 1$ -- competitive with contemporary redshift-space distortion constraints from galaxy surveys.

We measure the redshift-space correlation function from a spectroscopic sample of 2783 emission line galaxies from the FastSound survey. The survey, which uses the Subaru Telescope and covers the redshift ranges of $1.19<z<1.55$, is the first cosmological study at such high redshifts. We detect clear anisotropy due to redshift-space distortions (RSD) both in the correlation function as a function of separations parallel and perpendicular to the line of sight and its quadrupole moment. RSD has been extensively used to test general relativity on cosmological scales at $z<1$. Adopting a LCDM cosmology with the fixed expansion history and no velocity dispersion $\sigma_{\rm v}=0$, and using the RSD measurements on scales above 8Mpc/h, we obtain the first constraint on the growth rate at the redshift, $f(z)\sigma_8(z)=0.482\pm 0.116$ at $z\sim 1.4$ after marginalizing over the galaxy bias parameter $b(z)\sigma_8(z)$. This corresponds to $4.2\sigma$ detection of RSD. Our constraint is consistent with the prediction of general relativity $f\sigma_8\sim 0.392$ within the $1-\sigma$ confidence level. When we allow $\sigma_{\rm v}$ to vary and marginalize it over, the growth rate constraint becomes $f\sigma_8=0.494^{+0.126}_{-0.120}$. We also demonstrate that by combining with the low-z constraints on $f\sigma_8$, high-z galaxy surveys like the FastSound can be useful to distinguish modified gravity models without relying on CMB anisotropy experiments.

Homogeneous and isotropic, nonsingular, bouncing world models are designed to evade the initial singularity at the beginning of the cosmic expansion. Here, we study the thermodynamics of the subset of these models governed by general relativity. Considering the entropy of matter, radiation and that the entropy of the apparent horizon is proportional to its area, we argue that these models do not respect the generalised second law of thermodynamics, also away from the bounce.

We present a complete derivation of the observationally motivated definition of the modified gravity statistic $E_G$. Using this expression, we investigate how variations to theory and survey parameters may introduce uncertainty in the general relativistic prediction of $E_G$. We forecast errors on $E_G$ for measurements using two combinations of upcoming surveys, and find that theoretical uncertainties may dominate for a futuristic measurement. Finally, we compute predictions of $E_G$ under modifications to general relativity in the quasistatic regime, and comment on the pros and cons of using $E_G$ to test gravity with future surveys.

We introduce a designer approach for extended Brans-Dicke gravity that allows us to obtain the evolution of the scalar field by fixing the Hubble parameter to that of a $w$CDM model. We obtain analytical approximations for $\phi$ as a function of the scale factor and use these to build expressions for the effective Newton's constant at the background and at the linear level and the slip between the perturbed Newtonian potentials. By doing so, we are able to explore their dependence on the fundamental parameters of the theory.

Future surveys of large-scale structure will be able to measure perturbations on the scale of the cosmological horizon, and so could potentially probe a number of novel relativistic effects that are negligibly small on sub-horizon scales. These effects leave distinctive signatures in the power spectra of clustering observables and, if measurable, would open a new window on relativistic cosmology. We quantify the size and detectability of the effects for the most relevant future large-scale structure experiments: spectroscopic and photometric galaxy redshift surveys, intensity mapping surveys of neutral hydrogen, and radio continuum surveys. Our forecasts show that next-generation experiments, reaching out to redshifts $z\simeq 4$, will not be able to detect previously-undetected general-relativistic effects by using individual tracers of the density field, although the contribution of weak lensing magnification on large scales should be clearly detectable. We also perform a rigorous joint forecast for the detection of primordial non-Gaussianity through the excess power it produces in the clustering of biased tracers on large scales, finding that uncertainties of $\sigma(f_{\rm NL})\sim 1-2$ should be achievable. We study the level of degeneracy of these large-scale effects with several tracer-dependent nuisance parameters, quantifying the minimal priors on the latter that are needed for an optimal measurement of the former. Finally, we discuss the systematic effects that must be mitigated to achieve this level of sensitivity, and some alternative approaches that should help to improve the constraints. The computational tools developed to carry out this study, which requires the full-sky computation of the theoretical angular power spectra for $\mathcal{O}(100)$ redshift bins, as well as realistic models of the luminosity function, are publicly available.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

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.

As well as primary fluctuations, CMB temperature maps contain a wealth of additional information in the form of secondary anisotropies. Secondary effects that can be identified with individual objects, such as the thermal and kinetic Sunyaev-Zel'dovich (SZ) effects due to galaxy clusters, are difficult to unambiguously disentangle from foreground contamination and the primary CMB however. We develop a Bayesian formalism for rigorously characterising anisotropies that are localised on the sky, taking the TSZ and KSZ effects as an example. Using a Gibbs sampling scheme, we are able to efficiently sample from the joint posterior distribution for a multi-component model of the sky with many thousands of correlated physical parameters. The posterior can then be exactly marginalised to estimate properties of the secondary anisotropies, fully taking into account degeneracies with the other signals in the CMB map. We show that this method is computationally tractable using a simple implementation based on the existing Commander component separation code, and also discuss how other types of secondary anisotropy can be accommodated within our framework.

A number of alternatives to general relativity exhibit gravitational screening in the non-linear regime of structure formation. We describe a set of algorithms that can produce weak lensing maps of large scale structure in such theories and can be used to generate mock surveys for cosmological analysis. By analysing a few basic statistics we indicate how these alternatives can be distinguished from general relativity with future weak lensing surveys.

Multiple tracers of the cosmic density field, with different bias, number and luminosity evolution, can be used to measure the large-scale properties of the Universe. We show how an optimal combination of tracers can be used to detect general-relativistic effects in the observed density of sources. We forecast for the detectability of these effects, as well as measurements of primordial non-Gaussianity and large-scale lensing magnification with current and upcoming large-scale structure experiments. In particular we quantify the significance of these detections in the short term with experiments such as the Dark Energy Survey (DES), and in the long term with the Large Synoptic Survey Telescope (LSST) and the Square Kilometre Array (SKA). We review the main observational challenges that must be overcome to carry out these measurements.

We propose a simple and computationally fast method for performing N-body simulations for a large class of modified gravity theories with a screening mechanism such as chameleons, symmetrons and galileons. By combining the linear Klein-Gordon equation with a screening factor, calculated from analytical solutions of spherical symmetric configurations, we obtain a modified field equation whose solution is exact in the linear regime while at the same time takes screening into account on non-linear scales. The resulting modified field equation remains linear and can be solved just as quickly as the Poisson equation without any of the convergence problems that can arise when solving the full equation. We test our method with N-body simulations and find that it compares remarkably well with full simulations well into the non-linear regime.

Ultra-light axions (ULAs) with masses in the range 10^{-33} eV <m <10^{-20} eV are motivated by string theory and might contribute to either the dark-matter or dark-energy density of the Universe. ULAs could suppress the growth of structure on small scales, or lead to an enhanced integrated Sachs-Wolfe effect on large-scale cosmic microwave-background (CMB) anisotropies. In this work, cosmological observables over the full ULA mass range are computed, and then used to search for evidence of ULAs using CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP), Planck satellite, Atacama Cosmology Telescope, and South Pole Telescope, as well as galaxy clustering data from the WiggleZ galaxy-redshift survey. In the mass range 10^{-32} eV < m <10^{-25.5} eV, the axion relic-density \Omega_{a} (relative to the total dark-matter relic density \Omega_{d}) must obey the constraints \Omega_{a}/\Omega_{d} < 0.05 and \Omega_{a}h^{2} < 0.006 at 95%-confidence. For m> 10^{-24} eV, ULAs are indistinguishable from standard cold dark matter on the length scales probed, and are thus allowed by these data. For m < 10^{-32} eV, ULAs are allowed to compose a significant fraction of the dark energy.

Theories of modified gravity, in both the linear and fully non-linear regime, are often studied under the assumption that the evolution of the new (often scalar) degree of freedom present in the theory is quasi-static. This approximation significantly simplifies the study of the theory, and one often has good reason to believe that it should hold. Nevertheless it is a crucial assumption that should be explicitly checked whenever possible. In this paper we do so for the Vainshtein mechanism. By solving for the full spatial and time evolution of the Dvali-Gabadadze-Porrati and the Cubic Galileon model, in a spherical symmetric spacetime, we are able to demonstrate that the Vainshtein solution is a stable attractor and forms no matter what initial conditions we take for the scalar field. Furthermore,the quasi-static approximation is also found to be a very good approximation whenever it exists. For the best-fit Cubic Galileon model, however, we find that for deep voids at late times, the numerical solution blows up at the same time as the quasi-static solution ceases to exist. We argue that this phenomenon is a true instability of the model.

We present a framework for forecasting cosmological constraints from future neutral hydrogen intensity mapping experiments at low to intermediate redshifts. In the process, we establish a simple way of comparing such surveys with optical galaxy redshift surveys. We explore a wide range of experimental configurations and assess how well a number of cosmological observables (the expansion rate, growth rate, and angular diameter distance) and parameters (the densities of dark energy and dark matter, spatial curvature, the dark energy equation of state, etc.) will be measured by an extensive roster of upcoming experiments. A number of potential contaminants and systematic effects are also studied in detail. The overall picture is encouraging -- if autocorrelation calibration can be controlled to a sufficient level, Phase I of the SKA should be able to constrain the dark energy equation of state about as well as a DETF Stage IV galaxy redshift survey like Euclid, in roughly the same timeframe.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

By considering linear-order departures from general relativity, we compute a novel expression for the weak lensing convergence power spectrum under alternative theories of gravity. This comprises an integral over a 'kernel' of general relativistic quantities multiplied by a theory-dependent 'source' term. The clear separation between theory-independent and -dependent terms allows for an explicit understanding of each physical effect introduced by altering the theory of gravity. We take advantage of this to explore the degeneracies between gravitational parameters in weak lensing observations.

We construct the effective field theory of a perfect fluid in the early universe. Focusing on the case where the fluid has the equation of state of radiation, we show that it may lead to corrections to the background dynamics that can dominate over those of an effective field theory of gravity alone. We describe the periods of accelerated expansion, in the form of inflationary and bounce solutions, that arise in the background dynamics and discuss their regime of validity within EFT.

We present an analysis of the diffuse emission at 5 GHz in the first quadrant of the Galactic plane using two months of preliminary intensity data taken with the C-Band All Sky Survey (C-BASS) northern instrument at the Owens Valley Radio Observatory, California. Combining C-BASS maps with ancillary data to make temperature-temperature plots we find synchrotron spectral indices of $\beta = -2.65 \pm 0.05$ between 0.408 GHz and 5 GHz and $ \beta = -2.72 \pm 0.09$ between 1.420 GHz and 5 GHz for $-10^{\circ} < |b| < -4^{\circ}$, $20^{\circ} < l < 40^{\circ}$. Through the subtraction of a radio recombination line (RRL) free-free template we determine the synchrotron spectral index in the Galactic plane ($ |b| < 4^{\circ}$) to be $\beta = -2.56 \pm 0.07$ between 0.408 GHz and 5 GHz, with a contribution of $53 \pm 8$ per cent from free-free emission at 5\,GHz. These results are consistent with previous low frequency measurements in the Galactic plane. By including C-BASS data in spectral fits we demonstrate the presence of anomalous microwave emission (AME) associated with the HII complexes W43, W44 and W47 near 30 GHz, at 4.4 sigma, 3.1 sigma and 2.5 sigma respectively. The CORNISH VLA 5 GHz source catalogue rules out the possibility that the excess emission detected around 30\;GHz may be due to ultra-compact HII regions. Diffuse AME was also identified at a 4 sigma level within $30^{\circ} < l < 40^{\circ}$, $-2^{\circ} < b < 2^{\circ}$ between 5 GHz and 22.8 GHz.

The study of the Universe on ultra-large scales is one of the major science cases for the Square Kilometre Array (SKA). The SKA will be able to probe a vast volume of the cosmos, thus representing a unique instrument, amongst next-generation cosmological experiments, for scrutinising the Universe's properties on the largest cosmic scales. Probing cosmic structures on extremely large scales will have many advantages. For instance, the growth of perturbations is well understood for those modes, since it falls fully within the linear regime. Also, such scales are unaffected by the poorly understood feedback of baryonic physics. On ultra-large cosmic scales, two key effects become significant: primordial non-Gaussianity and relativistic corrections to cosmological observables. Moreover, if late-time acceleration is driven not by dark energy but by modifications to general relativity, then such modifications should become apparent near and above the horizon scale. As a result, the SKA is forecast to deliver transformational constraints on non-Gaussianity and to probe gravity on super-horizon scales for the first time.

HI intensity mapping (IM) is a novel technique capable of mapping the large-scale structure of the Universe in three dimensions and delivering exquisite constraints on cosmology, by using HI as a biased tracer of the dark matter density field. This is achieved by measuring the intensity of the redshifted 21cm line over the sky in a range of redshifts without the requirement to resolve individual galaxies. In this chapter, we investigate the potential of SKA1 to deliver HI intensity maps over a broad range of frequencies and a substantial fraction of the sky. By pinning down the baryon acoustic oscillation and redshift space distortion features in the matter power spectrum -- thus determining the expansion and growth history of the Universe -- these surveys can provide powerful tests of dark energy models and modifications to General Relativity. They can also be used to probe physics on extremely large scales, where precise measurements of spatial curvature and primordial non-Gaussianity can be used to test inflation; on small scales, by measuring the sum of neutrino masses; and at high redshifts where non-standard evolution models can be probed. We discuss the impact of foregrounds as well as various instrumental and survey design parameters on the achievable constraints. In particular we analyse the feasibility of using the SKA1 autocorrelations to probe the large-scale signal.

The imprint of baryon acoustic oscillations (BAO) in large-scale structure can be used as a standard ruler for mapping out the cosmic expansion history, and hence for testing cosmological models. In this article we briefly describe the scientific background to the BAO technique, and forecast the potential of the Phase 1 and 2 SKA telescopes to perform BAO surveys using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys. We find that a 25,000 sq. deg. intensity mapping survey on a Phase 1 array will preferentially constrain the radial BAO, providing a highly competitive 2% constraint on the expansion rate at z ~ 2. A 30,000 sq. deg. galaxy redshift survey on SKA2 will outperform all other planned experiments for z < 1.4.

The peculiar motion of galaxies can be a particularly sensitive probe of gravitational collapse. As such, it can be used to measure the dynamics of dark matter and dark energy as well the nature of the gravitational laws at play on cosmological scales. Peculiar motions manifest themselves as an overall anisotropy in the measured clustering signal as a function of the angle to the line-of-sight, known as redshift-space distortion (RSD). Limiting factors in this measurement include our ability to model non-linear galaxy motions on small scales and the complexities of galaxy bias. The anisotropy in the measured clustering pattern in redshift-space is also driven by the unknown distance factors at the redshift in question, the Alcock-Paczynski distortion. This weakens growth rate measurements, but permits an extra geometric probe of the Hubble expansion rate. In this chapter we will briefly describe the scientific background to the RSD technique, and forecast the potential of the SKA phase 1 and the SKA2 to measure the growth rate using both galaxy catalogues and intensity mapping, assessing their competitiveness with current and future optical galaxy surveys.

21cm intensity mapping experiments aim to observe the diffuse neutral hydrogen (HI) distribution on large scales which traces the Cosmic structure. The Square Kilometre Array (SKA) will have the capacity to measure the 21cm signal over a large fraction of the sky. However, the redshifted 21cm signal in the respective frequencies is faint compared to the Galactic foregrounds produced by synchrotron and free-free electron emission. In this article, we review selected foreground subtraction methods suggested to effectively separate the 21cm signal from the foregrounds with intensity mapping simulations or data. We simulate an intensity mapping experiment feasible with SKA phase 1 including extragalactic and Galactic foregrounds. We give an example of the residuals of the foreground subtraction with a independent component analysis and show that the angular power spectrum is recovered within the statistical errors on most scales. Additionally, the scale of the Baryon Acoustic Oscillations is shown to be unaffected by foreground subtraction.

The archetypal theory of dark energy is quintessence: a minimally coupled scalar field with a canonical kinetic energy and potential. By studying random potentials we show that quintessence imposes a restricted set of priors on the equation of state of dark energy. Focusing on the commonly-used parametrisation, $w(a)\approx w_0+w_a(1-a)$, we show that there is a natural scale and direction in the $(w_0, w_a)$ plane that distinguishes quintessence as a general framework. We calculate the expected information gain for a given survey and show that, because of the non-trivial prior information, it is a function of more than just the figure of merit. This allows us to make a quantitative case for novel survey strategies. We show that the scale of the prior sets target observational requirements for gaining significant information. This corresponds to a figure of merit FOM$\gtrsim 200$, a requirement that future galaxy redshift surveys will meet.

We present a comprehensive analysis of classical scalar, vector and tensor cosmological perturbations in ghost-free massive bigravity. In particular, we find the full evolution equations and analytical solutions in a wide range of regimes. We show that there are viable cosmological backgrounds but, as has been found in the literature, these models generally have exponential instabilities in linear perturbation theory. However, it is possible to find stable scalar cosmological perturbations for a very particular choice of parameters. For this stable subclass of models we find that vector and tensor perturbations have growing solutions. We argue that special initial conditions are needed for tensor modes in order to have a viable model.

In the quasistatic regime, generic modifications to gravity can give rise to novel scale-dependence of the gravitational field equations. Crucially, the detectability of the new scale-dependent terms hinges upon the existence of an effective mass scale or length scale at which corrections to General Relativity become relevant. Starting from only a few basic principles, we derive the general form of this scale-dependence. Our method recovers results previously known in the specific case of Horndeski gravity, but also shows that they are valid more generally, beyond the regime of scalar field theories. We forecast the constraints that upcoming experiments will place on the existence of a new fundamental mass scale or length scale in cosmology.

We revisit and extend the `Effective field theory for massive gravitons' constructed by Arkani-Hamed, Georgi and Schwartz in the light of recent progress in constructing ghost-free theories with multiple interacting spin-2 fields. We show that there exist several dual ways of restoring gauge invariance in such multi-gravity theories, find a generalised Fierz-Pauli tuning condition relevant in this context and highlight subtleties in demixing tensor and scalar modes. The generic multi-gravity feature of scalar mixing and its consequences for higher order interactions are discussed. In particular we show how the decoupling limit is qualitatively changed in theories of interacting spin-2 fields. We relate this to dRGT (de Rham, Gabadadze, Tolley) massive gravity, Hassan-Rosen bigravity and the multi-gravity constructions by Hinterbichler and Rosen. As an additional application we show that EBI (Eddington-Born-Infeld) bigravity and higher order generalisations thereof possess ghost-like instabilities.

We make use of a large set of fast simulations of an intensity mapping experiment with characteristics similar to those expected of the Square Kilometre Array (SKA) in order to study the viability and limits of blind foreground subtraction techniques. In particular, we consider different approaches: polynomial fitting, principal component analysis (PCA) and independent component analysis (ICA). We review the motivations and algorithms for the three methods, and show that they can all be described, using the same mathematical framework, as different approaches to the blind source separation problem. We study the efficiency of foreground subtraction both in the angular and radial (frequency) directions, as well as the dependence of this efficiency on different instrumental and modelling parameters. For well-behaved foregrounds and instrumental effects we find that foreground subtraction can be successful to a reasonable level on most scales of interest. We also quantify the effect that the cleaning has on the recovered signal and power spectra. Interestingly, we find that the three methods yield quantitatively similar results, with PCA and ICA being almost equivalent.

We present a code to generate mock observations of 21 cm intensity mapping experiments. The emphasis of the code is on reducing the computational cost of running a full-blown simulation, trading computational time for accuracy. The code can be used to generate independent realizations of the cosmological signal and foregrounds, which are necessary, for instance, in order to obtain realistic forecasts for future intensity mapping experiments. The code is able to reproduce the correct angular and radial clustering pattern for the cosmological signal, including redshift-space distortions, lightcone evolution and bias. Furthermore, it is possible to simulate a variety of foregrounds, including the potentially problematic polarized synchrotron emission.

We resurrect Eddington's proposal for the gravitational action in the presence of a cosmological constant and extend it to include matter fields. We show that the Newton-Poisson equation is modified in the presence of sources and that charged black holes show great similarities with those arising in Born-Infeld electrodynamics coupled to gravity. When we consider homogeneous and isotropic space-times we find that there is a minimum length (and maximum density) at early times, clearly pointing to an alternative theory of the Big Bang. We thus argue that the modern formulation of Eddington's theory, Born-Infeld gravity, presents us with a novel, non-singular description of the Universe.

The recent detection of B-modes by BICEP2 has non-trivial implications for axion dark matter implied by combining the tensor interpretation with isocurvature constraints from Planck. In this paper the measurement is taken as fact, and its implications considered, though further experimental verification is required. In the simplest inflation models $r=0.2$ implies $H_I=1.1\times 10^{14}\text{ GeV}$. If the axion decay constant $f_a<H_I/2\pi$ constraints on the dark matter (DM) abundance alone rule out the QCD axion as DM for $m_a \lesssim 52\chi^{6/7}\,\mu\text{eV}$ (where $\chi>1$ accounts for theoretical uncertainty). If $f_a>H_I/2\pi$ then vacuum fluctuations of the axion field place conflicting demands on axion DM: isocurvature constraints require a DM abundance which is too small to be reached when the back reaction of fluctuations is included. High $f_a$ QCD axions are thus ruled out. Constraints on axion-like particles, as a function of their mass and DM fraction, are also considered. For heavy axions with $m_a\gtrsim 10^{-22}\text{ eV}$ we find $\Omega_a/\Omega_d\lesssim 10^{-3}$, with stronger constraints on heavier axions. Lighter axions, however, are allowed and (inflationary) model-independent constraints from the CMB temperature power spectrum and large scale structure are stronger than those implied by tensor modes.

The growth rate of the large-scale structure of the universe has been advocated as the observable par excellence for testing gravity on cosmological scales. By considering linear-order deviations from General Relativity, we show that corrections to the growth rate, f, can be expressed as an integral over a `source' term, weighted by a theory-independent `response kernel'. This leads to an efficient and accurate `plug-and-play' expression for generating growth rates in alternative gravity theories, bypassing lengthy theory-specific computations. We use this approach to explicitly show that f is sensitive to a degenerate combination of modified expansion and modified clustering effects. Hence the growth rate, when used in isolation, is not a straightforward diagnostic of modified gravity.

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

There is evidence that Newton and Einstein's theories of gravity cannot explain the dynamics of a universe made up solely of baryons and radiation. To be able to understand the properties of galaxies, clusters of galaxies and the universe on the whole it has become commonplace to invoke the presence of dark matter. An alternative approach is to modify the gravitational field equations to accommodate observations. We propose a new class of gravitational theories in which we add a new degree of freedom, the Aether, in the form of a vector field that is coupled covariantly, but non-minimally, with the space-time metric. We explore the Newtonian and non-Newtonian limits, discuss the conditions for these theories to be consistent and explore their effect on cosmology.

Bekenstein's theory of relativistic gravity is conventionally written as a bi-metric theory. The two metrics are related by a disformal transformation defined by a dynamical vector field and a scalar field. In this comment we show that the theory can be re-written as Vector-Tensor theory akin to Einstein-Aether theories with non-canonical kinetic terms. We discuss some of the implications of this equivalence.

We propose and develop a general algorithm for finding the action for cosmological perturbations which rivals the conventional, gauge-invariant approach and can be applied to theories with more than one metric. We then apply it to a particular case of bigravity, focusing on the Eddington- inspired Born-Infeld theory, and show that we can obtain a nearly scale-invariant power spectrum for both scalar and tensor primordial quantum perturbations. Unfortunately, in the case of the minimal Eddington-inspired Born-Infeld theory, we find that the tensor-to-scalar ratio of perturbations is unacceptably large. We discuss the applicability of our general method and the possibility of resurrecting the specific theory we have looked at.

We outline the prospects for performing pioneering radio weak gravitational lensing analyses using observations from a potential forthcoming JVLA Sky Survey program. A large-scale survey with the JVLA can offer interesting and unique opportunities for performing weak lensing studies in the radio band, a field which has until now been the preserve of optical telescopes. In particular, the JVLA has the capacity for large, deep radio surveys with relatively high angular resolution, which are the key characteristics required for a successful weak lensing study. We highlight the potential advantages and unique aspects of performing weak lensing in the radio band. In particular, the inclusion of continuum polarisation information can greatly reduce noise in weak lensing reconstructions and can also remove the effects of intrinsic galaxy alignments, the key astrophysical systematic effect that limits weak lensing at all wavelengths. We identify a VLASS "deep fields" program (total area ~10-20 square degs), to be conducted at L-band and with high-resolution (A-array configuration), as the optimal survey strategy from the point of view of weak lensing science. Such a survey will build on the unique strengths of the JVLA and will remain unsurpassed in terms of its combination of resolution and sensitivity until the advent of the Square Kilometre Array. We identify the best fields on the JVLA-accessible sky from the point of view of overlapping with existing deep optical and near infra-red data which will provide crucial redshift information and facilitate a host of additional compelling multi-wavelength science.

We consider nonlocal gravity theories that include tensor nonlocalities. We show that in the cosmological context, the tensor nonlocalities, unlike scalar ones, generically give rise to growing modes. An explicit example with quadratic curvature terms is studied in detail. Possible consequences for recent nonlocal cosmological models proposed in the literature are also discussed.

The large-scale structure of the Universe supplies crucial information about the physical processes at play at early times. Unresolved maps of the intensity of 21 cm emission from neutral hydrogen HI at redshifts z~1-5 are the best hope of accessing the ultralarge-scale information, directly related to the early Universe. A purpose-built HI intensity experiment may be used to detect the large scale effects of primordial non-Gaussianity, placing stringent bounds on different models of inflation. We argue that it may be possible to place tight constraints on the non-Gaussianity parameter f_NL, with an error close to ~1.

Ultra-light axions ($m_a\lesssim 10^{-18}$eV), motivated by string theory, can be a powerful probe of the energy scale of inflation. In contrast to heavier axions the isocurvature modes in the ultra-light axions can coexist with observable gravitational waves. Here it is shown that large scale structure constraints severely limit the parameter space for axion mass, density fraction and isocurvature amplitude. It is also shown that radically different CMB observables for the ultra-light axion isocurvature mode additionally reduce this space. The results of a new, accurate and efficient method to calculate this isocurvature power spectrum are presented, and can be used to constrain ultra-light axions and inflation.

Relativistic scalar fields are ubiquitous in modified theories of gravity. An important tool in understanding their impact on structure formation, especially in the context of N-body simulations, is the quasi-static approximation in which the time evolution of perturbations in the scalar fields is discarded. We show that this approximation must be used with some care by studying linearly perturbed scalar field cosmologies and quantifying the errors that arise from taking the quasi-static limit. We focus on f(R) and chameleon models and link the accuracy of the quasi-static approximation to the fast/slow-roll behaviour of the background and its proximity to {\Lambda}CDM. Investigating a large range of scales, from super- to sub-horizon, we find that slow-rolling ({\Lambda}CDM-like) backgrounds generically result in good quasi-static behaviour, even on (super-)horizon scales. We also discuss how the approximation might affect studying the non-linear growth of structure in numerical N-body simulations.

An interesting approach to the cosmological coincidence problem is to allow dark matter and dark energy interact with each other also nongravitationally. We consider two general Ans\"{a}tze for such an interaction and appraise their ability to address the coincidence problem. We determine the average accuracy required on the cosmic expansion rate data to distinguish interacting cosmological models from the conventional $\Lambda$CDM scenario. We find that among the planned surveys the Wide Field Infrared Survey Telescope has the best chance to detect an interaction, though at a low significance level. To unambiguously determine the existence of an interaction one must, therefore, combine the said expansion data with other probes.

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.

This paper is concerned with string cosmology and the dynamics of multiple scalar fields in potentials that can become negative, and their features as (Early) Dark Energy models. Our point of departure is the "String Axiverse", a scenario that motivates the existence of cosmologically light axion fields as a generic consequence of string theory. We couple such an axion to its corresponding modulus. We give a detailed presentation of the rich cosmology of such a model, ranging from the setting of initial conditions on the fields during inflation, to the asymptotic future. We present some simplifying assumptions based on the fixing of the axion decay constant $f_a$, and on the effective field theory when the modulus trajectory is adiabatic, and find the conditions under which these assumptions break down. As a by-product of our analysis, we find that relaxing the assumption of fixed $f_a$ leads to the appearance of a new meta-stable de-Sitter region for the modulus without the need for uplifting by an additional constant. A dynamical systems analysis reveals the existence of many fixed point attractors, repellers and saddle points, which we analyse in detail. We also provide geometric interpretations of the phase space. The fixed points can be used to bound the couplings in the model. A systematic scan of certain regions of parameter space reveals that the future evolution of the universe in this model can be rich, containing multiple epochs of accelerated expansion.

A generic prediction of string theory is the existence of many axion fields. It has recently been argued that many of these fields should be light and, like the well known QCD axion, lead to observable cosmological consequences. In this paper we study in detail the effect of the so-called string axiverse on large scale structure, focusing on the morphology and evolution of density perturbations, anisotropies in the cosmic microwave background and weak gravitational lensing of distant galaxies. We quantify specific effects that will arise from the presence of the axionic fields and highlight possible degeneracies that may arise in the presence of massive neutrinos. We take particular care understanding the different physical effects and scales that come into play. We then forecast how the string axiverse may be constrained and show that with a combination of different observations, it should be possible to detect a fraction of ultralight axions to dark matter of a few percent.

A unified framework for theories of modified gravity will be an essential tool for interpreting the forthcoming deluge of cosmological data. We present such a formalism, the Parameterized Post-Friedmann framework (PPF), which parameterizes the cosmological perturbation theory of a wide variety of modified gravity models. PPF is able to handle spin-0 degrees of freedom from new scalar, vector and tensor fields, meaning that it is not restricted to simple models based solely on cosmological scalar fields. A direct correspondence is maintained between the parameterization and the underlying space of theories, which allows us to build up a `dictionary' of modified gravity theories and their PPF correspondences. In this paper we describe the construction of the parameterization and demonstrate its use through a number of worked examples relevant to the current literature. We indicate how the formalism will be implemented numerically, so that the dictionary of modified gravity can be pitted against forthcoming observations.

Simple models of gravitational arcs are crucial to simulate large samples of these objects with full control of the input parameters. These models also provide crude and automated estimates of the shape and structure of the arcs, which are necessary when trying to detect and characterize these objects on massive wide area imaging surveys. We here present and explore the ArcEllipse, a simple prescription to create objects with shape similar to gravitational arcs. We also present PaintArcs, which is a code that couples this geometrical form with a brightness distribution and adds the resulting object to images. Finally, we introduce ArcFitting, which is a tool that fits ArcEllipses to images of real gravitational arcs. We validate this fitting technique using simulated arcs and apply it to CFHTLS and HST images of tangential arcs around clusters of galaxies. Our simple ArcEllipse model for the arc, associated to a S\'ersic profile for the source, recovers the total signal in real images typically within 10%-30%. The ArcEllipse+S\'ersic models also automatically recover visual estimates of length-to-width ratios of real arcs. Residual maps between data and model images reveal the incidence of arc substructure. They may thus be used as a diagnostic for arcs formed by the merging of multiple images. The incidence of these substructures is the main factor preventing ArcEllipse models from accurately describing real lensed systems.

We study the dynamics of homogeneous, isotropic universes which are governed by the Eddington-inspired alternative theory of gravity which has a single extra parameter, $\kappa$. Previous results showing singularity-avoiding behaviour for $\kappa > 0$ are found to be upheld in the case of domination by a perfect fluid with equation of state parameter $w > 0$. The range $-1/3 < w < 0$ is found to lead to universes which experience unbounded expansion rate whilst still at a finite density. In the case $\kappa < 0$ the addition of spatial curvature is shown to lead to the possibility of oscillation between two finite densities. Domination by a scalar field with an exponential potential is found to also lead to singularity-avoiding behaviour when $\kappa > 0$. Certain values of the parameters governing the potential lead to behaviour in which the expansion rate of the universe changes sign several times before transitioning to regular GR-like behaviour.

A phenomenological attempt at alleviating the so-called coincidence problem is to allow the dark matter and dark energy to interact. By assuming a coupled quintessence scenario characterized by an interaction parameter $\epsilon$, we investigate the precision in the measurements of the expansion rate $H(z)$ required by future experiments in order to detect a possible deviation from the standard $\Lambda$CDM model ($\epsilon = 0$). We perform our analyses at two levels, namely: through Monte Carlo simulations based on $\epsilon$CDM models, in which $H(z)$ samples with different accuracies are generated and through an analytic method that calculates the error propagation of $\epsilon$ as a function of the error in $H(z)$. We show that our analytical approach traces simulations accurately and find that to detect an interaction {using $H(z)$ data only, these must reach an accuracy better than 1%.

The peculiar velocities of galaxies are an inherently valuable cosmological probe, providing an unbiased estimate of the distribution of matter on scales much larger than the depth of the survey. Much research interest has been motivated by the high dipole moment of our local peculiar velocity field, which suggests a large scale excess in the matter power spectrum, and can appear to be in some tension with the LCDM model. We use a composite catalogue of 4,537 peculiar velocity measurements with a characteristic depth of 33 h-1 Mpc to estimate the matter power spectrum. We compare the constraints with this method, directly studying the full peculiar velocity catalogue, to results from Macaulay et al. (2011), studying minimum variance moments of the velocity field, as calculated by Watkins, Feldman & Hudson (2009) and Feldman, Watkins & Hudson (2010). We find good agreement with the LCDM model on scales of k > 0.01 h Mpc-1. We find an excess of power on scales of k < 0.01 h Mpc-1, although with a 1 sigma uncertainty which includes the LCDM model. We find that the uncertainty in the excess at these scales is larger than an alternative result studying only moments of the velocity field, which is due to the minimum variance weights used to calculate the moments. At small scales, we are able to clearly discriminate between linear and nonlinear clustering in simulated peculiar velocity catalogues, and find some evidence (although less clear) for linear clustering in the real peculiar velocity data.

The Q/U Imaging ExperimenT (QUIET) employs coherent receivers at 43GHz and 95GHz, operating on the Chajnantor plateau in the Atacama Desert in Chile, to measure the anisotropy in the polarization of the CMB. QUIET primarily targets the B modes from primordial gravitational waves. The combination of these frequencies gives sensitivity to foreground contributions from diffuse Galactic synchrotron radiation. Between 2008 October and 2010 December, >10,000hours of data were collected, first with the 19-element 43GHz array (3458hours) and then with the 90-element 95GHz array. Each array observes the same four fields, selected for low foregrounds, together covering ~1000deg^2. This paper reports initial results from the 43GHz receiver which has an array sensitivity to CMB fluctuations of 69uK sqrt(s). The data were extensively studied with a large suite of null tests before the power spectra, determined with two independent pipelines, were examined. Analysis choices, including data selection, were modified until the null tests passed. Cross correlating maps with different telescope pointings is used to eliminate a bias. This paper reports the EE, BB and EB power spectra in the multipole range ell=25-475. With the exception of the lowest multipole bin for one of the fields, where a polarized foreground, consistent with Galactic synchrotron radiation, is detected with 3sigma significance, the E-mode spectrum is consistent with the LCDM model, confirming the only previous detection of the first acoustic peak. The B-mode spectrum is consistent with zero, leading to a measurement of the tensor-to-scalar ratio of r=0.35+1.06-0.87. The combination of a new time-stream double-demodulation technique, Mizuguchi-Dragone optics, natural sky rotation, and frequent boresight rotation leads to the lowest level of systematic contamination in the B-mode power so far reported, below the level of r=0.1

The Universe is inhomogeneous, and yet it seems to be incredibly well-characterised by a homogeneous relativistic model. One of the current challenges is to accurately characterise the properties of such a model. In this paper we explore how inhomogeneities may affect the overall optical properties of the Universe by quantifying how they can bias the redshift-distance relation in a number of toy models that mimic the real Universe. The models that we explore are statistically homogeneous on large scales. We find that the effect of inhomogeneities is of order of a few percent, which can be quite important in precise estimation of cosmological parameters. We discuss what lessons can be learned to help us tackle a more realistic inhomogeneous universe.

In this paper we consider an extension to Eddington's proposal for the gravitational action. We study tensor perturbations of a homogeneous and isotropic space-time in the Eddington regime, where modifications to Einstein gravity are strong. We find that the tensor mode is linearly unstable deep in the Eddington regime and discuss its cosmological implications.

There are a number of approaches to testing General Relativity (GR) on linear scales using parameterized frameworks for modifying cosmological perturbation theory. It is sometimes assumed that the details of any given parameterization are unimportant if one uses it as a diagnostic for deviations from GR. In this brief report we argue that this is not necessarily so. First we show that adopting alternative combinations of modifications to the field equations significantly changes the constraints that one obtains. In addition, we show that using a parameterization with insufficient freedom significantly tightens the apparent theoretical constraints. Fundamentally we argue that it is almost never appropriate to consider modifications to the perturbed Einstein equations as being constraints on the effective gravitational constant, for example, in the same sense that solar system constraints are. The only consistent modifications are either those that grant near-total freedom, as in decomposition methods, or ones which map directly to a particular part of theory space.

In this review we present a thoroughly comprehensive survey of recent work on modified theories of gravity and their cosmological consequences. Amongst other things, we cover General Relativity, Scalar-Tensor, Einstein-Aether, and Bimetric theories, as well as TeVeS, f(R), general higher-order theories, Horava-Lifschitz gravity, Galileons, Ghost Condensates, and models of extra dimensions including Kaluza-Klein, Randall-Sundrum, DGP, and higher co-dimension braneworlds. We also review attempts to construct a Parameterised Post-Friedmannian formalism, that can be used to constrain deviations from General Relativity in cosmology, and that is suitable for comparison with data on the largest scales. These subjects have been intensively studied over the past decade, largely motivated by rapid progress in the field of observational cosmology that now allows, for the first time, precision tests of fundamental physics on the scale of the observable Universe. The purpose of this review is to provide a reference tool for researchers and students in cosmology and gravitational physics, as well as a self-contained, comprehensive and up-to-date introduction to the subject as a whole.

It is discussed the possibility of a fine-tuneable contribution to the two way Doppler acceleration either towards, either outwards the Sun for heliocentric distances above 20 AU by considering a background described by an Expanding Locally Anisotropic (ELA) metric. This metric encodes both the standard local Schwarzschild gravitational effects and the cosmological Universe expansion effects allowing simultaneously to fine-tune other gravitational effects at intermediate scales, which may be tentatively interpreted as a covariant parameterization of either cold dark matter either gravitational interaction corrections. Are derived bounds for the ELA metric functional parameter by considering the bounds on the deviation from standard General Relativity imposed by the current updated limits for the Pioneer anomaly, taking in consideration both the natural outgassing and on-board radiation pressure, resulting in an average Doppler acceleration outwards the Sun of a_p = +0.4^{+2.1}_{-2.0} x 10^{-10} (m/s^2). It is also computed the mass-energy density for the ELA metric within the bounds obtained and are discussed the respective contributions to the cosmological mass-energy density which, for compatibility with the Lambda-CDM model, are included in Omega_{CDM}.

There is a distinct possibility that current and future cosmological data can be used to constrain Einstein's theory of gravity on the very largest scales. To be able to do this in a model-independent way, it makes sense to work with a general parameterization of modified gravity. Such an approach would be analogous to the Parameterized Post-Newtonian (PPN) approach which is used on the scale of the Solar System. A few such parameterizations have been proposed and preliminary constraints have been obtained. We show that the majority of such parameterizations are only exactly applicable in the quasistatic regime. On larger scales they fail to encapsulate the full behaviour of typical models currently under consideration. We suggest that it may be possible to capture the additions to the `Parameterized Post-Friedmann' (PPF) formalism by treating them akin to fluid perturbations.

The apparent accelerating expansion of the Universe, determined from observations of distant supernovae, and often taken to imply the existence of dark energy, may alternatively be explained by the effects of a giant underdense void if we relax the assumption of homogeneity on large scales. Recent studies have made use of the spherically-symmetric, radially-inhomogeneous Lemaitre-Tolman-Bondi (LTB) models to derive strong constraints on this scenario, particularly from observations of the kinematic Sunyaev-Zel'dovich (kSZ) effect which is sensitive to large scale inhomogeneity. However, most of these previous studies explicitly set the LTB 'bang time' function to be constant, neglecting an important freedom of the general solutions. Here we examine these models in full generality by relaxing this assumption. We find that although the extra freedom allowed by varying the bang time is sufficient to account for some observables individually, it is not enough to simultaneously explain the supernovae observations, the small-angle CMB, the local Hubble rate, and the kSZ effect. This set of observables is strongly constraining, and effectively rules out simple LTB models as an explanation of dark energy.

We consider the optical properties of Lindquist-Wheeler (LW) models of the Universe. These models consist of lattices constructed from regularly arranged discrete masses. They are akin to the Wigner-Seitz construction of solid state physics, and result in a dynamical description of the large-scale Universe in which the global expansion is given by a Friedmann-like equation. We show that if these models are constructed in a particular way then the redshifts of distant objects, as well as the dynamics of the global space-time, can be made to be in good agreement with the homogeneous and isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions of Einstein's equations, at the level of <3% out to z~2. Angular diameter and luminosity distances, on the other hand, differ from those found in the corresponding FLRW models, while being consistent with the 'empty beam' approximation, together with the shearing effects due to the nearest masses. This can be compared with the large deviations found from the corresponding FLRW values obtained in a previous study that considered LW models constructed in a different way. We therefore advocate the improved LW models we consider here as useful constructions that appear to faithfully reproduce both the dynamical and observational properties of space-times containing discrete masses.

We consider a model of the Universe in which the matter content is in the form of discrete islands, rather than a continuous fluid. In the appropriate limits the resulting large-scale dynamics approach those of a Friedmann-Robertson-Walker (FRW) universe. The optical properties of such a space-time, however, do not. This illustrates the fact that the optical and `average' dynamical properties of a relativistic universe are not equivalent, and do not specify each other uniquely. We find the angular diameter distance, luminosity distance and redshifts that would be measured by observers in these space-times, using both analytic approximations and numerical simulations. While different from their counterparts in FRW, the effects found do not look like promising candidates to explain the observations usually attributed to the existence of Dark Energy. This incongruity with standard FRW cosmology is not due to the existence of any unexpectedly large structures or voids in the Universe, but only to the fact that the matter content of the Universe is not a continuous fluid.

It is well known that observations of the cosmic microwave background (CMB) are highly sensitive to the spatial curvature of the Universe, k. Here we find that what is in fact being tightly constrained by small angle fluctuations is spatial curvature near the surface of last scattering, and that if we allow k to be a function of position, rather than taking a constant value everywhere, then considerable spatial curvature is permissible within our own locale. This result is of interest for the giant void models that attempt to explain the supernovae observations without Dark Energy. We find voids models with a homogeneous big bang can be compatible with the observed small angle CMB, but only if they exist in a positively curved universe. To be compatible with local measurements of H_0, however, we find that a radially varying bang time is required.

We propose a new statistic that has been designed to be used in situations where the intrinsic dispersion of a data set is not well known: The Crossing Statistic. This statistic is in general less sensitive than `chi^2' to the intrinsic dispersion of the data, and hence allows us to make progress in distinguishing between different models using goodness of fit to the data even when the errors involved are poorly understood. The proposed statistic makes use of the shape and trends of a model's predictions in a quantifiable manner. It is applicable to a variety of circumstances, although we consider it to be especially well suited to the task of distinguishing between different cosmological models using type Ia supernovae. We show that this statistic can easily distinguish between different models in cases where the `chi^2' statistic fails. We also show that the last mode of the Crossing Statistic is identical to `chi^2', so that it can be considered as a generalization of `chi^2'.

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.

The peculiar motions of galaxies can be used to infer the distribution of matter in the Universe. It has recently been shown that measurements of the peculiar velocity field indicates an anomalously high bulk flow of galaxies in our local volume. In this paper we find the implications of the high bulk flow for the power spectrum of density fluctuations. We find that analyzing only the dipole moment of the velocity field yields an average power spectrum amplitude which is indeed much higher than the LCDM value. However, by also including shear and octupole moments of the velocity field, and marginalizing over possible values for the growth rate, an average power spectrum amplitude which is consistent with LCDM is recovered. We attempt to infer the shape of the matter power spectrum from moments of the velocity field, and find a slight excess of power on scales ~ h-1 Gpc.

Ultra-light scalar fields, with masses of between m=10^{-33} eV and m=10^{-22} eV, can affect the growth of structure in the Universe. We identify the different regimes in the evolution of ultra-light scalar fields, how they affect the expansion rate of the universe and how they affect the growth rate of cosmological perturbations. We find a number of interesting effects, discuss how they might arise in realistic scenarios of the early universe and comment on how they might be observed.

It is suggested an expanding locally anisotropic metric (ELA) ansatz describing matter in a flat expanding universe which interpolates between the Schwarzschild (SC) metric near point-like central bodies of mass 'M' and the Robertson-Walker (RW) metric for large radial coordinate: 'ds^2=Z(cdt)2 - 1/Z (dr1-(Hr1/c) Z^(alpha/2+1/2)(cdt))^2-r1^2 dOmega', where 'Z=1-U' with 'U=2GM/(c^2r1)', 'G' is the Newton constant, 'c' is the speed of light, 'H=H(t)=\dot(a)/a' is the time-dependent Hubble rate, 'dOmega=dtheta^2+sin^2(theta) dvarphi^2' is the solid angle element, 'a' is the universe scale factor and we are employing the coordinates 'r1=ar', being 'r' the radial coordinate for which the RW metric is diagonal. For constant exponent 'alpha=alpha0=0' it is retrieved the isotropic McVittie (McV) metric and for 'alpha=alpha0=1' it is retrieved the locally anisotropic Cosmological-Schwarzschild (SCS) metric, both already discussed in the literature. However it is shown that only for constant exponent 'alpha=alpha0> 1' exists an event horizon at the SC radius 'r1=2GM/c^2' and only for 'alpha=alpha0>= 3' space-time is singularity free for this value of the radius. These bounds exclude the previous existing metrics, for which the SC radius is a naked extended singularity. In addition it is shown that for 'alpha=alpha0>5' space-time is approximately Ricci flat in a neighborhood of the event horizon such that the SC metric is a good approximation in this neighborhood. It is further shown that to strictly maintain the SC mass pole at the origin 'r1=0' without the presence of more severe singularities it is required a radial coordinate dependent correction to the exponent 'alpha(r1)=alpha0+alpha1 '2GM/(c^2 r1)' with a negative coefficient 'alpha1<0'. The energy-momentum density, pressures and equation of state are discussed.

A possible solution to the dark energy problem is that Einstein's theory of general relativity is modified. A suite of models have been proposed that, in general, are unable to predict the correct amount of large scale structure in the distribution of galaxies or anisotropies in the Cosmic Microwave Background. It has been argued, however, that it should be possible to constrain a general class of theories of modified gravity by focusing on properties such as the growing mode, gravitational slip and the effective, time varying Newton's constant. We show that assuming certain physical requirements such as stability, metricity and gauge invariance, it is possible to come up with consistency conditions between these various parameters. In this paper we focus on theories which have, at most, 2nd derivatives in the metric variables and find restrictions that shed light on current and future experimental constraints without having to resort to a (as yet unknown) complete theory of modified gravity. We claim that future measurements of the growth of structure on small scales (i.e. from 1-200 h^{-1} Mpc) may lead to tight constraints on both dark energy and modified theories of gravity.

Future cosmic microwave background (CMB) polarisation experiments aim to measure an unprecedentedly small signal - the primordial gravity wave component of the polarisation field B-mode. To achieve this, they will analyse huge datasets, involving years worth of time-ordered data (TOD) from massively multi-detector focal planes. This creates the need for fast and precise methods to complement the M-L approach in analysis pipelines. In this paper, we investigate fast map-making methods as applied to long duration, massively multi-detector, ground-based experiments, in the context of the search for B-modes. We focus on two alternative map-making approaches: destriping and TOD filtering, comparing their performance on simulated multi-detector polarisation data. We have written an optimised, parallel destriping code, the DEStriping CARTographer DESCART, that is generalised for massive focal planes, including the potential effect of cross-correlated TOD 1/f noise. We also determine the scaling of computing time for destriping as applied to a simulated full-season data-set for a realistic experiment. We find that destriping can out-perform filtering in estimating both the large-scale E and B-mode angular power spectra. In particular, filtering can produce significant spurious B-mode power via EB mixing. Whilst this can be removed, it contributes to the variance of B-mode bandpower estimates at scales near the primordial B-mode peak. For the experimental configuration we simulate, this has an effect on the possible detection significance for primordial B-modes. Destriping is a viable alternative fast method to the full M-L approach that does not cause the problems associated with filtering, and is flexible enough to fit into both M-L and Monte-Carlo pseudo-Cl pipelines.

We present a comprehensive investigation of cosmological constraints on the class of vector field formulations of modified gravity called Generalized Einstein-Aether models. Using linear perturbation theory we generate cosmic microwave background and large-scale structure spectra for general parameters of the theory, and then constrain them in various ways. We investigate two parameter regimes: a dark-matter candidate where the vector field sources structure formation, and a dark-energy candidate where it causes late-time acceleration. We find that the dark matter candidate does not fit the data, and identify five physical problems that can restrict this and other theories of dark matter. The dark energy candidate does fit the data, and we constrain its fundamental parameters; most notably we find that the theory's kinetic index parameter $n_{\mathrm{ae}}$ can differ significantly from its $\Lambda$CDM value.

The observed matter in the universe accounts for just 5 percent of the observed gravity. A possible explanation is that Newton's and Einstein's theories of gravity fail where gravity is either weak or enhanced. The modified theory of Newtonian dynamics (MOND) reproduces, without dark matter, spiral-galaxy orbital motions and the relation between luminosity and rotation in galaxies, although not in clusters. Recent extensions of Einstein's theory are theoretically more complete. They inevitably include dark fields that seed structure growth, and they may explain recent weak lensing data. However, the presence of dark fields reduces calculability and comes at the expense of the original MOND premise -- that the matter we see is the sole source of gravity. Observational tests of the relic radiation, weak lensing, and the growth of structure may distinguish modified gravity from dark matter.

The matter content of the Universe is strongly inhomogeneous on small scales. Motivated by this fact, we consider a model of the Universe that has regularly spaced discrete masses, rather than a continuous fluid. The optical properties of such space-times can differ considerably from the continuous fluid case, even if the 'average' dynamics are the same. We show that these differences have consequences for cosmological parameter estimation, and that fitting to recent supernovae observations gives a correction to the inferred value of Omega_Lambda of ~10%.

It is suggested a metric ansatz to describe local matter in an expanding universe, hence interpolating between the Schwarzschild metric at small spatial scales and the FLRW metric at large spatial scales. This is acomplished maintaining space-time free of singularities except for the Schwarzschild mass pole at the origin as opposed to metrics already considered in the literature with the same purpose, namely the McVittie metric. The modified Newton law is analyzed and the static orbit solutions computed. It is concluded that the effects of expantion in the solar system are negligible, however depending on the metric parameter value, at galactic scales there is a significant deviation from the General Relativity Newton law which may contribute to dark matter effects allowing for a flattening of galaxy rotation curves.

It has been argued that a Universe governed by Eddington-Born-Infeld gravity can be compatible with current cosmological constraints. The extra fields introduced in this theory can behave both as dark matter and dark energy, unifying the dark sector in one coherent framework. We show the various roles the extra fields can play in the expansion of the Universe and study the evolution of linear perturbations in the various regimes. We find that, as a unified theory of the dark sector, Eddington-Born-Infeld gravity will lead to excessive fluctuations in the Cosmic Microwave Background on large scales. In the presence of a cosmological constant, however, the extra fields can behave as a form of non-particulate dark matter and can lead to a cosmology which is entirely compatible with current observations of large scale structure. We discuss the interpretation of this form of dark matter and how it can differ from standard, particulate dark matter.

Map-making presents a significant computational challenge to the next generation of kilopixel CMB polarisation experiments. Years worth of time ordered data (TOD) from thousands of detectors will need to be compressed into maps of the T, Q and U Stokes parameters. Fundamental to the science goal of these experiments, the observation of B-modes, is the ability to control noise and systematics. In this paper, we consider an alternative to the maximum-likelihood method, called destriping, where the noise is modelled as a set of discrete offset functions and then subtracted from the time-stream. We compare our destriping code (Descart: the DEStriping CARTographer) to a full maximum-likelihood map-maker, applying them to 200 Monte-Carlo simulations of time-ordered data from a ground based, partial-sky polarisation modulation experiment. In these simulations, the noise is dominated by either detector or atmospheric 1/f noise. Using prior information of the power spectrum of this noise, we produce destriped maps of T, Q and U which are negligibly different from optimal. The method does not filter the signal or bias the E or B-mode power spectra. Depending on the length of the destriping baseline, the method delivers between 5 and 22 times improvement in computation time over the maximum-likelihood algorithm. We find that, for the specific case of single detector maps, it is essential to destripe the atmospheric 1/f in order to detect B-modes, even though the Q and U signals are modulated by a half-wave plate spinning at 5-Hz.

A fundamental presupposition of modern cosmology is the Copernican Principle; that we are not in a central, or otherwise special region of the Universe. Studies of Type Ia supernovae, together with the Copernican Principle, have led to the inference that the Universe is accelerating in its expansion. The usual explanation for this is that there must exist a `Dark Energy', to drive the acceleration. Alternatively, it could be the case that the Copernican Principle is invalid, and that the data has been interpreted within an inappropriate theoretical frame-work. If we were to live in a special place in the Universe, near the centre of a void where the local matter density is low, then the supernovae observations could be accounted for without the addition of dark energy. We show that the local redshift dependence of the luminosity distance can be used as a clear discriminant between these two paradigms. Future surveys of Type Ia supernovae that focus on a redshift range of ~0.1-0.4 will be ideally suited to test this hypothesis, and hence to observationally determine the validity of the Copernican Principle on new scales, as well as probing the degree to which dark energy must be considered a necessary ingredient in the Universe.

The Einstein-Aether theory provides a simple, dynamical mechanism for breaking Lorentz invariance. It does so within a generally covariant context and may emerge from quantum effects in more fundamental theories. The theory leads to a preferred frame and can have distinct experimental signatures. In this letter, we perform a comprehensive study of the cosmological effects of the Einstein-Aether theory and use observational data to constrain it. Allied to previously determined consistency and experimental constraints, we find that an Einstein-Aether universe can fit experimental data over a wide range of its parameter space, but requires a specific rescaling of the other cosmological densities.

We describe the objectives, design and predicted performance of Clover, which is a ground-based experiment to measure the faint ``B-mode'' polarisation pattern in the cosmic microwave background (CMB). To achieve this goal, clover will make polarimetric observations of approximately 1000 deg^2 of the sky in spectral bands centred on 97, 150 and 225 GHz. The observations will be made with a two-mirror compact range antenna fed by profiled corrugated horns. The telescope beam sizes for each band are 7.5, 5.5 and 5.5 arcmin, respectively. The polarisation of the sky will be measured with a rotating half-wave plate and stationary analyser, which will be an orthomode transducer. The sky coverage combined with the angular resolution will allow us to measure the angular power spectra between 20 < l < 1000. Each frequency band will employ 192 single polarisation, photon noise limited TES bolometers cooled to 100 mK. The background-limited sensitivity of these detector arrays will allow us to constrain the tensor-to-scalar ratio to 0.026 at 3sigma, assuming any polarised foreground signals can be subtracted with minimal degradation to the 150 GHz sensitivity. Systematic errors will be mitigated by modulating the polarisation of the sky signals with the rotating half-wave plate, fast azimuth scans and periodic telescope rotations about its boresight. The three spectral bands will be divided into two separate but nearly identical instruments - one for 97 GHz and another for 150 and 225 GHz. The two instruments will be sited on identical three-axis mounts in the Atacama Desert in Chile near Pampa la Bola. Observations are expected to begin in late 2009.

It has been suggested that Einstein's theory of General Relativity can be modified to accomodate mismatches between the gravitational field and luminous matter on a wide range of scales. Covariant theories of modified gravity generically predict the existence of extra degrees of freedom which may be interpreted as dark matter. We study a subclass of these theories where the overall energy density in these extra degrees of freedom is subdominant relative to the baryon density and show that they favour the presence of massive neutrinos. For some specific cases (such as a flat Universes with a cosmological constant) one finds a conservative lower bound on the neutrinos mass of $m_\nu>0.31$ eV.

We calculate the effect that a primordial homogeneous magnetic field, $\B_0$, will have on the different CMB power spectra due to Faraday rotation. Concentrating on the $TB$, $EB$ and $BB$ correlations, we forecast the ability for future CMB polarization experiments to constrain $\B_0$. Our results depend on how well the foregrounds can be subtracted from the CMB maps, but we find a predicted error between $\sigma_{\B_0} = 4 \times 10^{-11}$Gauss (for the QUIET experiment with foregrounds perfectly subtracted) and $3 \times 10^{-10}$Gauss (with the Clover experiment with no foreground subtraction). These constraints are two orders of magnitudes better than the present limits on $\B_0$.

It has been shown that generalized Einstein-Aether theories may lead to significant modifications to the non-relativistic limit of the Einstein equations. In this paper we study the effect of a general class of such theories on the Solar System. We consider corrections to the gravitational potential in negative and positive powers of distance from the source. Using measurements of the perihelion shift of Mercury and time delay of radar signals to Cassini, we place constraints on these corrections. We find that a subclass of generalized Einstein-Aether theories are compatible with these constraints.

We study the cosmological stability of a class of theories with a dynamical preferred frame. For a range of actions, we find cosmological solutions which are compatible with observations of the recent history of the Universe: a matter dominated era followed by accelerated expansion. We then study the evolution of linear perturbations on these backgrounds and find conditions on the parameters of the theory which allow for the growth of structure sourced by the new degrees of freedom.

We analyze the BOOMERanG 2003 (B03) 145 GHz temperature map to constrain the amplitude of a non Gaussian, primordial contribution to CMB fluctuations. We perform a pixel space analysis restricted to a portion of the map chosen in view of high sensitivity, very low foreground contamination and tight control of systematic effects. We set up an estimator based on the three Minkowski functionals which relies on high quality simulated data, including non Gaussian CMB maps. We find good agreement with the Gaussian hypothesis and derive the first limits based on BOOMERanG data for the non linear coupling parameter f_NL as -300<f_NL<650 at 68% CL and -800<f_NL<1050 at 95% CL.

It is widely accepted that the Universe underwent a period of thermal equilibrium at very early times. One expects a residue of this primordial state to be imprinted on the large scale structure of space time. In this paper we study the morphology of this thermal residue in a universe whose early dynamics is governed by a scalar field. We calculate the amplitude of fluctuations on large scales and compare it to the imprint of vacuum fluctuations. We then use the observed power spectrum of fluctuations on the cosmic microwave background to place a constraint on the temperature of the Universe before and during inflation. We also present an alternative scenario where the fluctuations are predominantly thermal and near scale-invariant.

Widrawn: The solutions presented in this work are not compatible with the equation of motion for g_{00}, which we did not properly verify. In the framework of this work, the only bulk solutions compatible with periodicity of the fifth dimension are constant fields. Also for a scalar field only constant solutions are obtained.

A range of experimental results point to the existence of a massive neutrino. The recent high precision measurements of the cosmic microwave background and the large scale surveys of galaxies can be used to place an upper bound on this mass. In this paper we perform a thorough analysis of all assumptions that go into obtaining a credible limit on $\sum m_{\nu}$. In particular we explore the impact of extending parameter space beyond the current standard cosmological model, the importance of priors and the uncertainties due to biasing in large scale structure. We find that the mass constraints are independent of the choice of parameterization as well as the inclusion of spatial curvature. The results of including the possibility of dark energy and tensors perturbations are shown to depend critically on the data sets used. The difference between an upper bound of 2.2 eV, assuming generic initial conditions, and an upper bound of 0.63 eV, assuming adiabaticity and a galaxy bias of 1, demonstrate the dependence of such a constraint on the assumptions in the analysis.

We study the background cosmology governed by the Tensor-Vector-Scalar theory of gravity proposed by Bekenstein. We consider a broad family of potentials that lead to modified gravity and calculate the evolution of the field variables both numerically and analytically. We find a range of possible behaviors, from scaling to the late time domination of either the additional gravitational degrees of freedom or the background fluid.

We present new measurements of the power spectra of the E-mode of CMB polarization, the temperature T, the cross-correlation of E and T, and upper limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager (CBI) observations. Both raw maps and optimal signal images in the uv-plane and real space show strong detections of the E-mode (11.7 sigma for the EE power spectrum overall) and no detection of the B-mode. The power spectra are used to constrain parameters of the flat tilted adiabatic Lambda-CDM models: those determined from EE and TE bandpowers agree with those from TT, a powerful consistency check. There is little tolerance for shifting polarization peaks from the TT-forecast locations, as measured by the angular sound crossing scale theta = 100 ell_s = 1.03 +/- 0.02 from EE and TE cf. 1.044 +/- 0.005 with the TT data included. The scope for extra out-of-phase peaks from subdominant isocurvature modes is also curtailed. The EE and TE measurements of CBI, DASI and BOOMERANG are mutually consistent, and, taken together rather than singly, give enhanced leverage for these tests.

A self consistent effective field theory of modified gravity has recently been proposed with spontaneous breaking of local Lorentz invariance. The symmetry is broken by a vector field with the wrong-sign mass term and it has been shown to have additional graviton modes and modified dispersion relations. In this paper we study the evolution of a homogeneous and isotropic universe in the presence of such a vector field with a minimum lying along the time-like direction. A plethora of different regimes is identified, such as accelerated expansion, loitering, collapse and tracking.

The trispectrum of the cosmic microwave background can be used to assess the level of non-Gaussianity on cosmological scales. It probes the fourth order moment, as a function of angular scale, of the probability distribution function of fluctuations and has been shown to be sensitive to primordial non-gaussianity, secondary anisotropies (such as the Ostriker-Vishniac effect) and systematic effects (such as astrophysical foregrounds). In this paper we develop a formalism for estimating the trispectrum from high resolution sky maps which incorporates the impact of finite sky coverage. This leads to a series of operations applied to the data set to minimize the effects of contamination due to the Gaussian component and correlations between estimates at different scales. To illustrate the effect of the estimation process, we apply our procedure to the BOOMERanG data set and show that it is consistent with Gaussianity. This work presents the first estimation of the CMB trispectrum on sub-degree scales.

A relativistic theory of modified gravity has been recently proposed by Bekenstein. The tensor field in Einstein's theory of gravity is replaced by a scalar, a vector, and a tensor field which interact in such a way to give Modified Newtonian Dynamics (MOND) in the weak-field non-relativistic limit. We study the evolution of the universe in such a theory, identifying its key properties and comparing it with the standard cosmology obtained in Einstein gravity. The evolution of the scalar field is akin to that of tracker quintessence fields. We expand the theory to linear order to find the evolution of perturbations on large scales. The impact on galaxy distributions and the cosmic microwave background is calculated in detail. We show that it may be possible to reproduce observations of the cosmic microwave background and galaxy distributions with Bekenstein's theory of MOND.

The Cosmic Microwave Background (CMB) anisotropy constrains the geometry of the Universe because the positions of the acoustic peaks of the angular power spectrum depend strongly on the curvature of underlying three-dimensional space. In this Letter we exploit current observations to determine the spatial geometry of the Universe in the presence of isocurvature modes. Previous analyses have always assumed that the cosmological perturbations were initially adiabatic. A priori one might expect that allowing additional isocurvature modes would substantially degrade the constraints on the curvature of the Universe. We find, however, that if one considers additional data sets, the geometry remains well constrained. When the most general isocurvature perturbation is allowed, the CMB alone can only poorly constrain the geometry to Omega_0=1.6+-0.3. Including large-scale structure (LSS) data one obtains Omega_0=1.07+-0.03, and Omega_0=1.06+-0.02 when supplemented by the Hubble Space Telescope (HST) Key Project determination of H_0 and SNIa data.

We use the isotropy of the Cosmic Microwave Background to place stringent constraints on a possible electrical charge asymmetry of the universe. We find the excess charge per baryon to be $q_{e-p}<10^{-26}e$ in the case of a uniform distribution of charge, where $e$ is the charge of the electron. If the charge asymmetry is inhomogeneous, the constraints will depend on the spectral index, $n$, of the induced magnetic field and range from $q_{e-p}<5\times 10^{-20}e$ ($n=-2$) to $q_{e-p}<2\times 10^{-26}e$ ($n\geq 2$). If one could further assume that the charge asymmetries of individual particle species are not anti-correlated so as to cancel, this would imply, for photons, $q_\gamma< 10^{-35}e$; for neutrinos, $q_\nu<4\times10^{-35}e$; and for heavy (light) dark matter particles $q_{\rm dm}<4\times10^{-24}e$ ($q_{\rm dm}<4\times10^{-30}e$).

A simple method is presented for the rapid simulation of statistically-isotropic non-Gaussian maps of CMB temperature fluctuations with a given power spectrum and analytically-calculable bispectrum and higher-order polyspectra. The $n$th-order correlators of the pixel values may also be calculated analytically. The cumulants of the simulated map may be used to obtain an expression for the probability density function of the pixel temperatures. The statistical properties of the simulated map are determined by the univariate non-Gaussian distribution from which pixel values are drawn independently in the first stage of the simulation process. We illustrate the method using a non-Gaussian distribution derived from the wavefunctions of the harmonic oscillator. The basic simulation method is easily extended to produce non-Gaussian maps with a given power spectrum and diagonal bispectrum.

Observations of the CMB can be contaminated by diffuse foreground emission from sources such as Galactic dust and synchrotron radiation. In these cases, the morphology of the contaminating source is known from observations at different frequencies, but not its amplitude at the frequency of interest for the CMB. We develop a technique for accounting for the effects of such emission in this case, and for simultaneously estimating the foreground amplitude in the CMB observations. We apply the technique to CMB data from the MAXIMA-1 experiment, using maps of Galactic dust emission from combinations of IRAS and DIRBE observations, as well as compilations of Galactic synchrotron emission observations. The spectrum of the dust emission over the 150--450 GHz observed by MAXIMA is consistent with preferred models but the effect on CMB power spectrum observations is negligible.

We investigate the constraints imposed by the first-year WMAP CMB data extended to higher multipole by data from ACBAR, BOOMERANG, CBI and the VSA and by the LSS data from the 2dF galaxy redshift survey on the possible amplitude of primordial isocurvature modes. A flat universe with CDM and Lambda is assumed, and the baryon, CDM (CI), and neutrino density (NID) and velocity (NIV) isocurvature modes are considered. Constraints on the allowed isocurvature contributions are established from the data for various combinations of the adiabatic mode and one, two, and three isocurvature modes, with intermode cross-correlations allowed. Since baryon and CDM isocurvature are observationally virtually indistinguishable, these modes are not considered separately. We find that when just a single isocurvature mode is added, the present data allows an isocurvature fraction as large as 13+-6, 7+-4, and 13+-7 percent for adiabatic plus the CI, NID, and NIV modes, respectively. When two isocurvature modes plus the adiabatic mode and cross-correlations are allowed, these percentages rise to 47+-16, 34+-12, and 44+-12 for the combinations CI+NID, CI+NIV, and NID+NIV, respectively. Finally, when all three isocurvature modes and cross-correlations are allowed, the admissible isocurvature fraction rises to 57+-9 per cent. The sensitivity of the results to the choice of prior probability distribution is examined.

Markov Chain Monte Carlo (MCMC) techniques are now widely used for cosmological parameter estimation. Chains are generated to sample the posterior probability distribution obtained following the Bayesian approach. An important issue is how to optimize the efficiency of such sampling and how to diagnose whether a finite-length chain has adequately sampled the underlying posterior probability distribution. We show how the power spectrum of a single such finite chain may be used as a convergence diagnostic by means of a fitting function, and discuss strategies for optimizing the distribution for the proposed steps. The methods developed are applied to current CMB and LSS data interpreted using both a pure adiabatic cosmological model and a mixed adiabatic/isocurvature cosmological model including possible correlations between modes. For the latter application, because of the increased dimensionality and the presence of degeneracies, the need for tuning MCMC methods for maximum efficiency becomes particularly acute.

We measure the angular power spectrum of radio galaxies in the NRAO VLA Sky Survey (NVSS) using two independent methods: direct spherical harmonic analysis and maximum likelihood estimation. The results are consistent and can be understood using models for the spatial matter power spectrum and for the redshift distribution of radio galaxies at mJy flux-density levels. A good fit to the angular power spectrum can only be achieved if radio galaxies possess high bias with respect to mass fluctuations; by marginalizing over the other parameters of the model we derive a 68% confidence interval 1.53 < b_0 sigma_8 < 1.87, where b_0 is the linear bias factor for radio galaxies and sigma_8 describes the normalization of the matter power spectrum. Our models indicate that the majority of the signal in the NVSS angular power spectrum is generated at low redshifts. Individual redshifts for the NVSS sources are thus required to alleviate projection effects and probe directly the matter power spectrum on large scales.

We investigate the constraints imposed by the current data on correlated mixtures of adiabatic and non-adiabatic primordial perturbations. We discover subtle flat directions in parameter space that tolerate large (~60%) contributions of non-adiabatic fluctuations. In particular, larger values of the baryon density and a spectral tilt are allowed. The cancellations in the degenerate directions are explored and the role of priors elucidated.

The temporal behaviour of the radial abundance gradients has important consequences on models for the chemical evolution of the Galaxy. We present a comparison of the time variation of the abundance gradients in the Milky Way disk as determined from a sample of planetary nebulae and open clusters. We conclude that the [Fe/H] gradients as measured in open cluster stars strongly support the time flattening of the abundance gradient as recently determined from O/H measurements in planetary nebulae. We estimate the average flattening rate in the time interval 10 Myr < t < 5 Gyr as roughly 0.005 dex kpc-1 Gyr-1.

We cross-correlate the cosmic microwave background temperature anisotropy maps from the WMAP, MAXIMA-I, and MAXIMA-II experiments. We use the cross-spectrum, which is the spherical harmonic transform of the angular two-point correlation function, to quantify the correlation as a function of angular scale. We find that the three possible pairs of cross-spectra are in close agreement with each other and with the power spectra of the individual maps. The probability that there is no correlation between the maps is smaller than 1 * 10^(-8). We also calculate power spectra for maps made of differences between pairs of maps, and show that they are consistent with no signal. The results conclusively show that the three experiments not only display the same statistical properties of the CMB anisotropy, but also detect the same features wherever the observed sky areas overlap. We conclude that the contribution of systematic errors to these maps is negligible and that MAXIMA and WMAP have accurately mapped the cosmic microwave background anisotropy.

Using mean relative peculiar velocity measurements for pairs of galaxies, we estimate the cosmological density parameter $\Omega_m$ and the amplitude of density fluctuations $\sigma_8$. Our results suggest that our statistic is a robust and reproducible measure of the mean pairwise velocity and thereby the $\Omega_m$ parameter. We get $\Omega_m = 0.30^{+0.17}_{-0.07}$ and $\sigma_8 = 1.13^{+0.22}_{-0.23}$. These estimates do not depend on prior assumptions on the adiabaticity of the initial density fluctuations, the ionization history, or the values of other cosmological parameters.

We describe the Millimeter wave Anisotropy eXperiment IMaging Array (MAXIMA), a balloon-borne experiment designed to measure the temperature anisotropy of the Cosmic Microwave Background (CMB) on angular scales of 10' to 5 degrees . MAXIMA mapped the CMB using 16 bolometric detectors in spectral bands centered at 150 GHz, 240 GHz, and 410 GHz, with 10' resolution at all frequencies. The combined receiver sensitivity to CMB anisotropy was ~40 microK/rt(sec). Systematic parasitic contributions were minimized by using four uncorrelated spatial modulations, thorough crosslinking, multiple independent CMB observations, heavily baffled optics, and strong spectral discrimination. Pointing reconstruction was accurate to 1', and absolute calibration was better than 4%. Two MAXIMA flights with more than 8.5 hours of CMB observations have mapped a total of 300 deg^2 of the sky in regions of negligible known foreground emission. MAXIMA results have been released in previous publications. MAXIMA maps, power spectra and correlation matrices are publicly available at http://cosmology.berkeley.edu/maxima

The MAXIMA cosmic microwave background anisotropy experiment had a significant impact on cosmology. Results from the program have played a significant role in determining the geometry of the universe, given strong supporting evidence to inflation, and, in combination with other astrophysical data, showed that the universe is filled with dark matter and energy. We present a subset of the internal consistency checks that were carried out on the MAXIMA-1 data prior to their release, which demonstrate that systematics errors were much smaller than statistical errors. We also discuss the MAXIMA-2 flight and data, compare the maps of MAXIMA-1 and -2 in areas where they overlap and show that the two independent experiments confirm each other. All of these results demonstrate that MAXIMA mapped the cosmic microwave background anisotropy with high accuracy

We review the physics of the Grishchuck-Zel'dovich effect which describes the impact of large amplitude, super-horizon gravitational field fluctuations on the Cosmic Microwave Background anisotropy power spectrum. Using the latest determination of the spectrum by WMAP, we infer a lower limit on the present length-scale of such fluctuations of 3927 times the cosmological particle horizon (at the 95% confidence level).

We discuss MAXIPOL, a bolometric balloon-borne experiment designed to measure the E-mode polarization anisotropy of the cosmic microwave background radiation (CMB) on angular scales of 10 arcmin to 2 degrees. MAXIPOL is the first CMB experiment to collect data with a polarimeter that utilizes a rotating half-wave plate and fixed wire-grid polarizer. We present the instrument design, elaborate on the polarimeter strategy and show the instrument performance during flight with some time domain data. Our primary data set was collected during a 26 hour turnaround flight that was launched from the National Scientific Ballooning Facility in Ft. Sumner, New Mexico in May 2003. During this flight five regions of the sky were mapped. Data analysis is in progress.

MAXIMA is a balloon-borne platform for measuring the anisotropy of the Cosmic Microwave Background (CMB). It has measured the CMB power spectrum with a ten-arcminute FWHM beam, corresponding to a detection of the power spectrum out to spherical harmonic multipole l~1000. The spectrum is consistent with a flat Universe with a nearly scale-invariant initial spectrum of adiabatic density fluctuations. Moreover, the MAXIMA data are free from any notable non-Gaussian contamination and from foreground dust emission. In the same region, the WMAP experiment observes the same structure as that observed by MAXIMA, as evinced by analysis of both maps and power spectra. The next step in the evolution of the MAXIMA program is MAXIPOL, which will observe the polarization of the CMB with comparable resolution and high sensitivity over a small patch of the sky.

Although the broad outlines of the appropriate pipeline for cosmological likelihood analysis with CMB data has been known for several years, only recently have we had to contend with the full, large-scale, computationally challenging problem involving both highly-correlated noise and extremely large datasets ($N > 1000$). In this talk we concentrate on the beginning and end of this process. First, we discuss estimating the noise covariance from the data itself in a rigorous and unbiased way; this is essentially an iterated minimum-variance mapmaking approach. We also discuss the unbiased determination of cosmological parameters from estimates of the power spectrum or experimental bandpowers.

We describe different methods for estimating the bispectrum of Cosmic Microwave Background data. In particular we construct a minimum variance estimator for the flat-sky limit and compare results with previously-studied frequentist methods. Application to the MAXIMA dataset shows consistency with primordial Gaussianity. Weak quadratic non-Gaussianity is characterised by a tunable parameter $f_{NL}$, corresponding to non-Gaussianity at a level $\sim 10^{-5}f_{NL}$ (ratio of non-Gaussian to Gaussian terms), and we find limits of $|f_{NL}|<950$ for the minimum-variance estimator and $|f_{NL}|<1650$ for the usual frequentist estimator. These are the tightest limits on primordial non-Gaussianity which include the full effects of the radiation transfer function.

We consider the most general parametrization of flat topologically compact universes, complementing the work of Scannapieco, Levin and Silk to include non-trivial shapes. We find that modifications in shape of the fundamental domain will lead to distinct signatures in the anisotropy of the cosmic microwave radiation. We make a preliminary assessment of the effect on three statistics: the angular power spectrum, the distribution of identified ``circles'' on the surface of last scattering and the correlation function of antipodal points.

We use a frequentist statistical approach to set confidence intervals on the values of cosmological parameters using the MAXIMA-1 and COBE measurements of the angular power spectrum of the cosmic microwave background. We define a $\Delta \chi^{2}$ statistic, simulate the measurements of MAXIMA-1 and COBE, determine the probability distribution of the statistic, and use it and the data to set confidence intervals on several cosmological parameters. We compare the frequentist confidence intervals to Bayesian credible regions. The frequentist and Bayesian approaches give best estimates for the parameters that agree within 15%, and confidence interval-widths that agree within 30%. The results also suggest that a frequentist analysis gives slightly broader confidence intervals than a Bayesian analysis. The frequentist analysis gives values of \Omega=0.89{+0.26\atop -0.19}, \Omega_{\rm B}h^2=0.026{+0.020\atop -0.011} and n=1.02{+0.31\atop -0.10}, and the Bayesian analysis gives values of \Omega=0.98{+0.14\atop -0.19}, \Omega_{\rm B}h^2=0.0.029{+0.015\atop-0.010}, and $n=1.18{+0.10\atop -0.23}$, all at the 95% confidence level.

We use the measurement of the cosmic microwave background taken during the MAXIMA-1 flight to estimate the bispectrum of cosmological perturbations. We propose an estimator for the bispectrum that is appropriate in the flat sky approximation, apply it to the MAXIMA-1 data and evaluate errors using bootstrap methods. We compare the estimated value with what would be expected if the sky signal were Gaussian and find that it is indeed consistent, with a $\chi^2$ per degree of freedom of approximately unity. This measurement places constraints on models of inflation.

We propose an estimator for the trispectrum of a scalar random field on a sphere, discuss its geometrical and statistical properties, and outline its implementation. By estimating the trispectrum of the 4 year COBE-DMR data (in HEALPix pixelization) we find new evidence of a non-Gaussian signal associated with a known systematic effect. We find that by removing data from the sky maps for those periods of time perturbed by this effect, the amplitudes of the trispectrum coefficients become completely consistent with predictions for a Gaussian sky. These results reinforce the importance of statistical methods based in harmonic space for quantifying non-Gaussianity.

This work describes Cosmic Microwave Background (CMB) data analysis algorithms and their implementations, developed to produce a pixelized map of the sky and a corresponding pixel-pixel noise correlation matrix from time ordered data for a CMB mapping experiment. We discuss in turn algorithms for estimating noise properties from the time ordered data, techniques for manipulating the time ordered data, and a number of variants of the maximum likelihood map-making procedure. We pay particular attention to issues pertinent to real CMB data, and present ways of incorporating them within the framework of maximum likelihood map-making. Making a map of the sky is shown to be not only an intermediate step rendering an image of the sky, but also an important diagnostic stage, when tests for and/or removal of systematic effects can efficiently be performed. The case under study is the MAXIMA data set. However, the methods discussed are expected to be applicable to the analysis of other current and forthcoming CMB experiments.

We discuss the cosmological implications of the new constraints on the power spectrum of the Cosmic Microwave Background Anisotropy derived from a new high resolution analysis of the MAXIMA-1 measurement (Lee et al. 2001). The power spectrum shows excess power at $\ell \sim 860$ over the average level of power at $411 \le\ell \le 785.$ This excess is statistically significant on the 95% confidence level. Such a feature is consistent with the presence of a third acoustic peak, which is a generic prediction of inflation-based models. The height and the position of the excess power match the predictions of a family of inflationary models with cosmological parameters that are fixed to fit the CMB data previously provided by BOOMERANG-LDB and MAXIMA-1 experiments (e.g., Jaffe et al.2001). Our results, therefore, lend support for inflationary models and more generally for the dominance of coherent perturbations in the structure formation of the Universe. At the same time, they seem to disfavor a large variety of the non-standard (but still inflation-based) models that have been proposed to improve the quality of fits to the CMB data and consistency with other cosmological observables. Within standard inflationary models, our results combined with the COBE-DMR data give best fit values and 95% confidence limits for the baryon density, $\Omega_b h^2\simeq 0.033{\pm 0.013}$, and the total density, $\Omega=0.9{+0.18\atop -0.16}$. The primordial spectrum slope ($n_s$) and the optical depth to the last scattering surface ($\tau_c$) are found to be degenerate and to obey the relation $n_s \simeq 0.46 \tau_c + (0.99 \pm 0.14)$, for $\tau_c \le 0.5$ (all 95% c.l.).

The angular power spectrum of the Cosmic Microwave Background has been measured out to sufficiently small angular scale to encompass a few acoustic oscillations. We use a phenomenological fit to the angular power spectrum to quantify the statistical significance of these oscillations and discuss the cosmological implications of such a finding.

We set new constraints on a seven-dimensional space of cosmological parameters within the class of inflationary adiabatic models. We use the angular power spectrum of the cosmic microwave background measured over a wide range of \ell in the first flight of the MAXIMA balloon-borne experiment (MAXIMA-1) and the low \ell results from COBE/DMR. We find constraints on the total energy density of the universe, \Omega=1.0^{+0.15}_{-0.30}, the physical density of baryons, \Omega_{b}h^2=0.03 +/- 0.01, the physical density of cold dark matter, \Omega_{cdm}h^2=0.2^{+0.2}_{-0.1}$, and the spectral index of primordial scalar fluctuations, n_s=1.08+/-0.1, all at the 95% confidence level. By combining our results with measurements of high-redshift supernovae we constrain the value of the cosmological constant and the fractional amount of pressureless matter in the universe to 0.45<\Omega_\Lambda<0.75 and 0.25<\Omega_{m}<0.50, at the 95% confidence level. Our results are consistent with a flat universe and the shape parameter deduced from large scale structure, and in marginal agreement with the baryon density from big bang nucleosynthesis.

We reexamine Wald's no-hair theorem for global anisotropy in the brane world scenarios. We derive a set of sufficient conditions which must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane asymptotically evolves to a de Sitter spacetime in the presence of a positive cosmological constant on the brane. We discuss the violations of these sufficient conditions and we show that a negative nonlocal energy density or the presence of strong anisotropic stress (i.e., a magnetic field) may lead the brane to collapse. We discuss the generality of these conditions.

We extend the analysis of the MAXIMA-1 cosmic microwave background (CMB) data to smaller angular scales. MAXIMA, a bolometric balloon experiment, mapped a 124 deg$^2$ region of the sky with 10\arcmin resolution at frequencies of 150, 240 and 410 GHz during its first flight. The original analysis, which covered the multipole range $36 \leq \ell \leq 785$, is extended to $\ell = 1235$ using data from three 150 GHz photometers in the fully cross-linked central 60 deg$^2$ of the map. The main improvement over the original analysis is the use of 3\arcmin square pixels in the calculation of the map. The new analysis is consistent with the original for $\ell < 785$. For $\ell > 785$, where inflationary models predict a third acoustic peak, the new analysis shows power with an amplitude of $56 \pm 7$ \microk at $\ell \simeq 850$ in excess to the average power of $42 \pm 3$ \microk in the range $441 < \ell < 785$.

Recent results from BOOMERANG-98 and MAXIMA-1, taken together with COBE-DMR, provide consistent and high signal-to-noise measurements of the CMB power spectrum at spherical harmonic multipole bands over $2<\ell\lta800$. Analysis of the combined data yields 68% (95%) confidence limits on the total density, $\Omega_{\rm {tot}}\simeq 1.11 \pm 0.07 (^{+0.13}_{-0.12})$, the baryon density, $\Omega_b h^2\simeq 0.032^{+0.005}_{-0.004} (^{+0.009}_{-0.008})$, and the scalar spectral tilt, $n_s\simeq1.01^{+0.09}_{-0.07} (^{+0.17}_{-0.14})$. These data are consistent with inflationary initial conditions for structure formation. Taken together with other cosmological observations, they imply the existence of both non-baryonic dark matter and dark energy in the universe.

Gaussianity of the cosmological perturbations in the universe is one of the key predictions of standard inflation, but it is violated by other models of structure formation such as topological defects. We present the first test of the Gaussianity of the cosmic microwave background (CMB) on sub-degree angular scales, where deviations from Gaussianity are most likely to occur. We apply the methods of moments, cumulants, the Kolmogorov test, the $\chi^2$ test, and Minkowski functionals in eigen, real, Wiener-filtered and signal-whitened spaces, to the MAXIMA-1 CMB anisotropy data. We conclude that the data, which probe angular scales between 10 arcminutes and 5 degrees, are consistent with Gaussianity.

We describe here an iterative method for jointly estimating the noise power spectrum from a CMB experiment's time-ordered data, together with the maximum-likelihood map. We test the robustness of this method on simulated Boomerang datasets with realistic noise.

We review selected methods of the Cosmic Microwave Background data analysis appropriate for the analysis of the largest currently available data sets. We focus on techniques of the time-ordered data manipulation and map making algorithms based on the maximum-likelihood approach. The presented methods have been applied to the MAXIMA data analysis (Hanany et al 2000) and the description of the algorithms is illustrated with the examples drawn from that experience. The more extensive presentation of the here-mentioned issues will be given in the forthcoming paper (Stompor et al 2001).

BOOMERanG has recently resolved structures on the last scattering surface at redshift $\sim$ 1100 with high signal to noise ratio. We review the technical advances which made this possible, and we focus on the current results for maps and power spectra, with special attention to the determination of the total mass-energy density in the Universe and of other cosmological parameters.

We report the first results from the BOOMERanG experiment, which mapped at 90, 150, 240 and 410 GHz a wide (3%) region of the microwave sky with minimal local contamination. From the data of the best 150 GHz detector we find evidence for a well defined peak in the power spectrum of temperature fluctuations of the Cosmic Microwave Background, localized at $\ell = 197 \pm 6$, with an amplitude of $(68 \pm 8) \mu K_{CMB}$. The location, width and amplitude of the peak is suggestive of acoustic oscillations in the primeval plasma. In the framework of inflationary adiabatic cosmological models the measured spectrum allows a Bayesian estimate of the curvature of the Universe and of other cosmological parameters. With reasonable priors we find $\Omega = (1.07 \pm 0.06)$ and $n_s = (1.00 \pm 0.08)$ (68%C.L.) in excellent agreement with the expectations from the simplest inflationary theories. We also discuss the limits on the density of baryons, of cold dark matter and on the cosmological constant.

We show how estimates of parameters characterizing inflation-based theories of structure formation localized over the past year when large scale structure (LSS) information from galaxy and cluster surveys was combined with the rapidly developing cosmic microwave background (CMB) data, especially from the recent Boomerang and Maxima balloon experiments. All current CMB data plus a relatively weak prior probability on the Hubble constant, age and LSS points to little mean curvature (Omega_{tot} = 1.08\pm 0.06) and nearly scale invariant initial fluctuations (n_s =1.03\pm 0.08), both predictions of (non-baroque) inflation theory. We emphasize the role that degeneracy among parameters in the L_{pk} = 212\pm 7 position of the (first acoustic) peak plays in defining the $\Omega_{tot}$ range upon marginalization over other variables. Though the CDM density is in the expected range (\Omega_{cdm}h^2=0.17\pm 0.02), the baryon density Omega_bh^2=0.030\pm 0.005 is somewhat above the independent 0.019\pm 0.002 nucleosynthesis estimate. CMB+LSS gives independent evidence for dark energy (Omega_\Lambda=0.66\pm 0.06) at the same level as from supernova (SN1) observations, with a phenomenological quintessence equation of state limited by SN1+CMB+LSS to w_Q<-0.7 cf. the w_Q=-1 cosmological constant case.

The past, present and future of cosmic microwave background (CMB) anisotropy research is discussed, with emphasis on the Boomerang and Maxima balloon experiments. These data are combined with large scale structure (LSS) information and high redshift supernova (SN1) observations to explore the inflation-based cosmic structure formation paradigm. Here we primarily focus on a simplified inflation parameter set, {omega_b,omega_{cdm},Omega_{tot}, Omega_Q,w_Q, n_s,tau_C, sigma_8}. After marginalizing over the other cosmic and experimental variables, we find the current CMB+LSS+SN1 data gives Omega_{tot}=1.04\pm 0.05, consistent with (non-baroque) inflation theory. Restricting to Omega_{tot}=1, we find a nearly scale invariant spectrum, n_s =1.03 \pm 0.07. The CDM density, omega_{cdm}=0.17\pm 0.02, is in the expected range, but the baryon density, omega_b=0.030\pm 0.004, is slightly larger than the current nucleosynthesis estimate. Substantial dark energy is inferred, Omega_Q\approx 0.68\pm 0.05, and CMB+LSS Omega_Q values are compatible with the independent SN1 estimates. The dark energy equation of state, parameterized by a quintessence-field pressure-to-density ratio w_Q, is not well determined by CMB+LSS (w_Q<-0.3 at 95%CL), but when combined with SN1 the resulting w_Q<-0.7 limit is quite consistent with the w_Q=-1 cosmological constant case. Though forecasts of statistical errors on parameters for current and future experiments are rosy, rooting out systematic errors will define the true progress.

We describe the implications of cosmic microwave background (CMB) observations and galaxy and cluster surveys of large scale structure (LSS) for theories of cosmic structure formation, especially emphasizing the recent Boomerang and Maxima CMB balloon experiments. The inflation-based cosmic structure formation paradigm we have been operating with for two decades has never been in better shape. Here we primarily focus on a simplified inflation parameter set, {omega_b,omega_{cdm},Omega_{tot}, Omega_\Lambda,n_s,\tau_C, \sigma_8}. Combining all of the current CMB+LSS data points to the remarkable conclusion that the local Hubble patch we can access has little mean curvature (Omega_{tot}=1.08\pm 0.06) and the initial fluctuations were nearly scale invariant (n_s=1.03\pm 0.08), both predictions of (non-baroque) inflation theory. The baryon density is found to be slightly larger than that preferred by independent Big Bang Nucleosynthesis estimates (omega_b=0.030\pm 0.005 cf. 0.019\pm 0.002). The CDM density is in the expected range (omega_{cdm}=0.17 \pm 0.02). Even stranger is the CMB+LSS evidence that the density of the universe is dominated by unclustered energy akin to the cosmological constant (Omega_\Lambda=0.66\pm 0.06), at the same level as that inferred from high redshift supernova observations. We also sketch the CMB+LSS implications for massive neutrinos.

We present a map and an angular power spectrum of the anisotropy of the cosmic microwave background (CMB) from the first flight of MAXIMA. MAXIMA is a balloon-borne experiment with an array of 16 bolometric photometers operated at 100 mK. MAXIMA observed a 124 square degrees region of the sky with 10 arcminute resolution at frequencies of 150, 240 and 410 GHz. The data were calibrated using in-flight measurements of the CMB dipole anisotropy. A map of the CMB anisotropy was produced from three 150 and one 240 GHz photometer without need for foreground subtractions. Analysis of this CMB map yields a power spectrum for the CMB anisotropy over the range 36 < l < 785. The spectrum shows a peak with an amplitude of 78 +/- 6 micro-Kelvin at l ~ 220 and an amplitude varying between ~40 micro-Kelvin and ~50 micro-Kelvin for 400 < l < 785.

We propose a new formalism to handle asymmetric beams in the data analysis of cosmic microwave background anisotropy experiments. For any beam shape, the formalism finds the optimal circularly symmetric equivalent and is thus easily adaptable to existing data analysis methods. We demonstrate certain key points by using a simulated highly elliptic beam, and the beams and data of the MAXIMA-1 experiment, where the asymmetry is mild. In particular, we show that in both cases the formalism does not bias the angular power spectrum estimates. We analyze the limitations of the formalism and find that it is well suited for most practical situations.

The anisotropy of the cosmic microwave background radiation contains information about the contents and history of the universe. We report new limits on cosmological parameters derived from the angular power spectrum measured in the first Antarctic flight of the BOOMERANG experiment. Within the framework of inflation-motivated adiabatic cold dark matter models, and using only weakly restrictive prior probabilites on the age of the universe and the Hubble expansion parameter $h$, we find that the curvature is consistent with flat and that the primordial fluctuation spectrum is consistent with scale invariant, in agreement with the basic inflation paradigm. We find that the data prefer a baryon density $\Omega_b h^2$ above, though similar to, the estimates from light element abundances and big bang nucleosynthesis. When combined with large scale structure observations, the BOOMERANG data provide clear detections of both dark matter and dark energy contributions to the total energy density $\Omega_{\rm {tot}}$, independent of data from high redshift supernovae.

The blackbody radiation left over from the Big Bang has been transformed by the expansion of the Universe into the nearly isotropic 2.73K Cosmic Microwave Background. Tiny inhomogeneities in the early Universe left their imprint on the microwave background in the form of small anisotropies in its temperature. These anisotropies contain information about basic cosmological parameters, particularly the total energy density and curvature of the universe. Here we report the first images of resolved structure in the microwave background anisotropies over a significant part of the sky. Maps at four frequencies clearly distinguish the microwave background from foreground emission. We compute the angular power spectrum of the microwave background, and find a peak at Legendre multipole $\ell_{peak}=(197 \pm 6)$, with an amplitude $DT_{200}=(69 \pm 8)\mu K$. This is consistent with that expected for cold dark matter models in a flat (euclidean) Universe, as favoured by standard inflationary scenarios.

The motions of galaxies can be used to constrain the cosmological density parameter Omega and the clustering amplitude of matter on large scales. The mean relative velocity of galaxy pairs, estimated from the Mark III survey, indicates that Omega = 0.35 +0.35/-0.25. If the clustering of galaxies is unbiased on large scales, Omega = 0.35 +/- 0.15, so that an unbiased Einstein-de Sitter model (Omega = 1) is inconsistent with the data.

It is well known that estimating the pairwise velocity of galaxies, v_{12}, from the redshift space galaxy correlation function is difficult because this method is highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate v_{12} directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. In contrast to other dynamical measures which determine beta = sigma_8 x Omega^{0.6}, our method can provide an estimate of (sigma_8)^2 x Omega^{0.6} for a range of sigma_8 (here Omega is the cosmological mass density parameter while sigma_8 is the standard normalization parameter for the spectrum of matter density fluctuations). We demonstrate how to measure this quantity from realistic catalogues.

We use the angular power spectrum of the Cosmic Microwave Background, measured during the North American test flight of the BOOMERANG experiment, to constrain the geometry of the universe. Within the class of Cold Dark Matter models, we find that the overall fractional energy density of the universe, Omega, is constrained to be 0.85 < Omega < 1.25 at the 68% confidence level. Combined with the COBE measurement and the high redshift supernovae data we obtain new constraints on the fractional matter density and the cosmological constant.

We describe a measurement of the angular power spectrum of anisotropies in the Cosmic Microwave Background (CMB) from 0.3 degrees to ~10 degrees from the North American test flight of the BOOMERANG experiment. BOOMERANG is a balloon-borne telescope with a bolometric receiver designed to map CMB anisotropies on a Long Duration Balloon flight. During a 6-hour test flight of a prototype system in 1997, we mapped > 200 square degrees at high galactic latitudes in two bands centered at 90 and 150 GHz with a resolution of 26 and 16.6 arcmin FWHM respectively. Analysis of the maps gives a power spectrum with a peak at angular scales of ~1 degree with an amplitude ~70 uK.

We derive an expression for the angular power spectrum of cosmic microwave background anisotropies due to gravity waves generated by a stochastic magnetic field and compare the result with current observations; we take into account the non-linear nature of the stress energy tensor of the magnetic field. For almost scale invariant spectra, the amplitude of the magnetic field at galactic scales is constrained to be of order 10^{-9} Gauss. If we assume that the magnetic field is damped below the Alfven damping scale, we find that its amplitude at 0.1 h^{-1}Mpc, B_\lambda, is constrained to be B_\lambda<7.9 x10^{-6} e^{3n} Gauss, for n<-3/2, and B_\lambda<9.5x10^{-8} e^{0.37n} Gauss, for n>-3/2, where n is the spectral index of the magnetic field and H_0=100h km s^{-1}Mpc^{-1} is the Hubble constant today.

We propose a formalism for estimating the skewness and angular power spectrum of a general Cosmic Microwave Background data set. We use the Edgeworth Expansion to define a non-Gaussian likelihood function that takes into account the anisotropic nature of the noise and the incompleteness of the sky coverage. The formalism is then applied to estimate the skewness of the publicly available 4 year Cosmic Background Explorer (COBE) Differential Microwave Radiometer data. We find that the data is consistent with a Gaussian skewness, and with isotropy. Inclusion of non Gaussian degrees of freedom has essentially no effect on estimates of the power spectrum, if each $C_\ell$ is regarded as a separate parameter or if the angular power spectrum is parametrized in terms of an amplitude (Q) and spectral index (n). Fixing the value of the angular power spectrum at its maxiumum likelihood estimate, the best fit skewness is $S=6.5\pm6.0\times10^4(\muK)^3$; marginalizing over Q the estimate of the skewness is $S=6.5\pm8.4\times10^4(\muK)^3$ and marginalizing over n one has $S=6.5\pm8.5\times10^4(\muK)^3$.

To correctly analyse data sets from current microwave detection technology, one is forced to estimate the sky signal and experimental noise simultaneously. Given a time-ordered data set we propose a formalism and method for estimating the signal and associated errors without prior knowledge of the noise power spectrum. We derive the method using a Bayesian formalism and relate it to the standard methods; in particular we show how this leads to a change in the estimate of the noise covariance matrix of the sky signal. We study the convergence and accuracy of the method on two mock observational strategies and discuss its application to a currently-favoured calibration procedure.

The mean pairwise velocity of galaxies has traditionally been estimated from the redshift space galaxy correlation function. This method is notorious for being highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate the streaming velocity directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. This method can provide an estimate of $\Omega^{0.6}\sigma_8^2$ for a range of $\sigma_8$ where $\Omega$ is the cosmological density parameter, while $\sigma_8$ is the standard normalization for the power spectrum of density fluctuations. We demonstrate how to measure this quantity from realistic catalogues and identify the main sources of bias and errors

I review our recent claim that there is evidence of non-Gaussianity in the 4 Year COBE DMR data. I describe the statistic we apply, the result we obtain and make a detailed list of the systematics we have analysed. I finish with a qualitative understanding of what it might be and its implications.

We describe the MAXIMA experiment, a balloon-borne measurement designed to map temperature anisotropy in the Cosmic Microwave Background (CMB) over a wide range of angular scales (multipole range 80 < l < 800). The experiment consists of a 1.3 m diameter off-axis Gregorian telescope and a receiver with a 16 element array of bolometers cooled to 100 mK. The frequency bands are centered at 150, 240, and 410 GHz. The 10' FWHM beam sizes are well matched to the scale of acoustic peaks expected in the angular power spectrum of the CMB. The first flight of the experiment in its full configuration was launched in August 1998. A 122 sq-deg map of the sky was made near the Draco constellation during the 7 hour flight in a region of extremely low galactic dust contamination. This map covers 0.3% of the sky and has 3200 independent beamsize pixels. We describe the MAXIMA instrument and its performance during the recent flight.

We review our recent claim that there is evidence of non-Gaussianity in the 4 Year COBE DMR data. We present some new results concerning the effect of the galactic cut upon the non-Gaussian signal. These findings imply a localization of the non-Gaussian signal on the Northern galactic hemisphere.

In a recent Letter we have shown how COBE-DMR maps may be used to disprove Gaussianity at a high confidence level. In this report we digress on a few issues closely related to this Letter. We present the general formalism for surveying non-Gaussianity employed. We present a few more tests for systematics. We wonder about the theoretical implications of our result.

A scalar field with an exponential potential has the particular property that it is attracted into a solution in which its energy scales as the dominant component (radiation or matter) of the Universe, contributing a fixed fraction of the total energy density. We briefly discuss the dynamics of such a scalar field and its impact on Big Bang nucleosynthesis, the growth of large scale structure and abundance of damped Lyman$-\alpha$ systems at high redshift. Given the simplicity of the model, its theoretical motivation, and its success in matching observations, we argue that it should be taken on par with other currently viable models of structure formation.

A weakly coupled scalar field $\Phi$ with a simple exponential potential $V=M_P^4\exp(-\lambda\Phi/M_P)$ where $M_P$ is the reduced Planck mass, and $\lambda > 2$, has an attractor solution in a radiation or matter dominated universe in which it mimics the scaling of the dominant component, contributing a fixed fraction $\Omega_\phi$ (determined by $\lambda$) to the energy density. Such fields arise generically in particle physics theories involving compactified dimensions, with values of $\lambda$ which give a cosmologically relevant $\Omega_\phi$. For natural initial conditions on the scalar field in the early universe the attractor solution is established long before the epoch of structure formation, and in contrast to the solutions used in other scalar field cosmologies, it is one which does not involve an energy scale for the scalar field characteristic of late times . We study in some detail the evolution of matter and radiation perturbations in a standard inflation-motivated $\Omega=1$ dark-matter dominated cosmology with this extra field. Using a full Einstein-Boltzmann calculation we compare observable quantities with current data. We find that, for $\Omega_\phi\simeq 0.08-0.12$, these models are consistent with large angle cosmic microwave background anisotropies as detected by COBE, the linear mass variance as compiled from galaxy surveys, big bang nucleosynthesis, the abundance of rich clusters and constraints from the Lyman-$\alpha$ systems at high redshift. Given the simplicity of the model, its theoretical motivation and its success in matching observations, we argue that it should be taken on a par with other currently viable models of structure formation.

We introduce and study the distribution of an estimator for the normalized bispectrum of the Cosmic Microwave Background (CMB) anisotropy. We use it to construct a goodness of fit statistic to test the coadded 53 and 90 GHz COBE-DMR 4 year maps for non-Gaussianity. Our results indicate that Gaussianity is ruled out at the confidence level in excess of 98%. This value is a lower bound, given all the investigated systematics. The dominant non-Gaussian contribution is found near the multipole of order $\ell=16$. Our attempts to explain this effect as caused by the diffuse foreground emission from the Galaxy have failed. We conclude that unless there exists a microwave foreground emission which spatially correlates neither with the DIRBE nor Haslam maps, the cosmological CMB anisotropy is genuinely non-Gaussian.

A scalar field with an exponential potential has the particular property that it is attracted into a solution in which its energy scales as the dominant component (radiation or matter) of the Universe, contributing a fixed fraction of the total energy density. We study the growth of perturbations in a CDM dominated $\Omega=1$ universe with this extra field, with an initial flat spectrum of adiabatic fluctuations. The observational constraints from structure formation are satisfied as well, or better, than in other models, with a contribution to the energy density from the scalar field $\Omega_\phi \sim 0.1$ which is small enough to be consistent with entry into the attractor prior to nucleosynthesis.

We propose a new method for constraining a primordial homogeneous magnetic field with the cosmic microwave background. Such a field will induce an observable parity odd cross correlation between the polarization anisotropy and the temperature anisotropy by Faraday rotation. We analyze the necessary experimental features to match, and improve, current constraints of such a field by measuring this correlation.

We discuss the requirements of good statistics for quantifying non-Gaussianity in the Cosmic Microwave Background. The importance of rotational invariance and statistical independence is stressed, but we show that these are sometimes incompatible. It is shown that the first of these requirements prefers a real space (or wavelet) formulation, whereas the latter favours quantities defined in Fourier space. Bearing this in mind we decide to be eclectic and define two new sets of statistics to quantify the level of non-Gaussianity. Both sets make use of the concept of cumulants of a distribution. However, one set is defined in real space, with reference to the wavelet transform, whereas the other is defined in Fourier space. We derive a series of properties concerning these statistics for a Gaussian random field and show how one can relate these quantities to the higher order moments of temperature maps. Although our frameworks lead to an infinite hierarchy of quantities we show how cosmic variance and experimental constraints give a natural truncation of this hierarchy. We then focus on the real space statistics and analyse the non-Gaussian signal generated by points sources obscured by large scale Gaussian fluctuations. We conclude by discussing the practical implementations of these techniques.

In this paper we explore the connection between anisotropic Gaussian fluctuations and isotropic non-Gaussian fluctuations. We first set up a large angle framework for characterizing non-Gaussian fluctuations: large angle non-Gaussian spectra. We then consider anisotropic Gaussian fluctuations in two different situations. Firstly we look at anisotropic space-times and propose a prescription for superimposed Gaussian fluctuations; we argue against accidental symmetry in the fluctuations and that therefore the fluctuations should be anisotropic. We show how these fluctuations display previously known non-Gaussian effects both in the angular power spectrum and in non-Gaussian spectra. Secondly we consider the anisotropic Grischuk-Zel'dovich effect. We construct a flat space time with anisotropic, non-trivial topology and show how Gaussian fluctuations in such a space-time look non-Gaussian. In particular we show how non-Gaussian spectra may probe superhorizon anisotropy.

We propose a set of new statistics which can be extracted out of the angular distribution of the Fourier transform of the temperature anisotropies in the small field limit. They quantify generic non-Gaussian structure and complement the power spectrum in characterizing the sampled distribution function of a data set (contribution to the 18th Texas Symposium Proceedings)

We derive an upper limit of $B_0<3.4\times 10^{-9}(\Omega_0h_{50}^2)^{1/2}$ Gauss on the present strength of any primordial homogeneous magnetic field. The microwave background anisotropy created by cosmological magnetic fields is calculated in the most general flat and open anisotropic cosmologies containing expansion-rate and 3-curvature anisotropies. Our limit is derived from a statistical analysis of the 4-year Cosmic Background Explorer data for anisotropy patterns characteristic of homogeneous anisotropy averaged over all possible sky orientations with respect to the COBE receiver. The limits we obtain are considerably stronger than those imposed by primordial nucleosynthesis and ensure that other magnetic field effects on the microwave background structure are unobservably small.

Gaussian cosmic microwave background skies are fully specified by the power spectrum. The conventional method of characterizing non-Gaussian skies is to evaluate higher order moments, the n-point functions and their Fourier transforms. We argue that this method is inefficient, due to the redundancy of information existing in the complete set of moments. In this paper we propose a set of new statistics or non-Gaussian spectra to be extracted out of the angular distribution of the Fourier transform of the temperature anisotropies in the small field limit. These statistics complement the power spectrum and act as localization, shape, and connectedness statistics. They quantify generic non-Gaussian structure, and may be used in more general image processing tasks. We concentrate on a subset of these statistics and argue that while they carry no information in Gaussian theories they may be the best arena for making predictions in some non-Gaussian theories. As examples of applications we consider superposed Gaussian and non-Gaussian signals, such as point sources in Gaussian theories or the realistic Kaiser-Stebbins effect. We show that in these theories non-Gaussianity is only present in a ring in Fourier space, which is best isolated in our formalism. Subtle but strongly non-Gaussian theories are also written down for which only non-Gaussian spectra may accuse non-Gaussianity.

I report on recent developments in the theory of cosmic background radiation perturbations. I describe ways of modeling alternatives to the canonical Gaussian theories within the standard framework of cosmological perturbation theory. Some comments are made on using these techniques to resolve the uncertainties in theories of structure formation with topological defects. (To appear in the proceedings of the XXXIth Moriond meeting, ``Microwave Background Anisotropies'')

Large-scale cosmic microwave background anisotropies in homogeneous, globally anisotropic cosmologies are investigated. We perform a statistical analysis in which the four-year data from the Cosmic Background Explorer satellite is searched for the specific anisotropy patterns predicted by these models and thereby set definitive upper limits on the amount of shear, $(\sigma/H)_0$ and vorticity, $(\omega/H)_0$, which are orders of magnitude stronger than previous constraints. We comment on how these results might impact our understanding of primordial global anisotropy.

We investigate how the qualitative structure of Doppler peaks in the angular power spectrum of the cosmic microwave anisotropy is affected by basic assumptions going into theories of structure formation. We define the concepts of ``coherent'' and ``incoherent'' fluctuations, and also of ``active'' and ``passive'' fluctuations. In these terms inflationary fluctuations are passive and coherent while topological defects are active incoherent fluctuations. Causality and scale invariance are shown to have different implementations in theories differing in the above senses. We then extend the formalism of Hu and Sugiyama to treat models with cosmic defects. Using this formalism we show that the existence or absence of secondary Doppler peaks and the rough placing of the primary peak are very sensitive to the fundamental properties defined. We claim therefore that even a rough measurement of the angular power spectrum $C_l$ shape at $100<l<1500$ ought to tell us which are the basic ingredients to be used in the right structure formation theory. We also apply our formalism to cosmic string theories. These are shown to fall into the class of active incoherent theories for which one can robustly predict the absence of secondary Doppler peaks. The placing of the cosmic strings' primary peak is more uncertain, but should fall in $l\approx400-600$.

Fluctuations in the cosmic microwave background (CMB) temperature are being studied with ever increasing precision. Two competing types of theories might describe the origins of these fluctuations: ``inflation'' and ``defects''. Here we show how the differences between these two scenarios can give rise to striking signatures in the microwave fluctuations on small scales, assuming a standard recombination history. These should enable high resolution measurements of CMB anisotropies to distinguish between these two broad classes of theories, independent of the precise details of each.

We examine how the qualitative structure of the Doppler peaks in the angular power spectrum of the cosmic microwave anisotropy depends on the fundamental nature of the perturbations which produced them. The formalism of Hu and Sugiyama is extended to treat models with cosmic defects. We discuss how perturbations can be ``active'' or ``passive'' and ``incoherent'' or ``coherent'', and show how causality and scale invariance play rather different roles in these various cases. We find that the existence of secondary Doppler peaks and the rough placing of the primary peak unambiguously reflect these basic properties.

We study the effects of cosmic strings on structure formation in open universes. We calculate the power spectrum of density perturbations for two class of models: one in which all the dark matter is non baryonic (CDM) and one in which it is all baryonic (BDM). Our results are compared to the 1 in 6 IRAS QDOT power spectrum. The best candidates are then used to estimate $\mu$, the energy per unit length of the string network. Some comments are made on mechanisms by which structures are formed in the two theories.

High resolution maps of the anisotropy of the microwave sky will yield invaluable clues as to the mechanisms involved in cosmic structure formation. One fundamental question they should answer is whether the fluctuations were Gaussian random noise, as predicted in inflationary models, or were nonGaussian as in theories based on symmetry breaking and cosmic defects. In the latter case there is the prospect of obtaining information regarding symmetry breaking at high energy scales, specifically the homotopy classes $\Pi_n$ of the vacuum manifold. In this paper we report on detailed calculations of the degree scale anisotropies predicted in the cosmic string, monopole texture and nontopological texture theories of structure formation, emphasising their distinct character from those predicted by inflation, and the bright prospects for experimental tests.

The inflationary cosmology is analyzed from the point of view of squeezed quantum states. As noted by Grishchuk and Sidorov, the amplification of quantum fluctuations into macroscopic perturbations which occurs during cosmic inflation is a process of quantum squeezing. We carefully develop the squeezed state formalism and derive the equations that govern the evolution of a gaussian initial state. We derive the power spectrum of density perturbations for a simple inflationary model and discuss its features. We conclude that the squeezed state formalism provides an interesting framework within which to study the amplification process, but,in disagreement with the claims of Grishchuk and Sidorov, that it does {\em not} provide us with any new physical results.