The Two Higgs Doublet Model predicts the emergence of 3 distinct domain wall solutions arising from the breaking of 3 accidental global symmetries, $Z_2$, CP1 and CP2, at the electroweak scale for specific choices of the model parameters. We present numerical kink solutions to the field equations in all three cases along with dynamical simulations of the models in (2+1) and (3+1) dimensions. For each kink solution we define an associated topological current. In all three cases simulations produce a network of domain walls which deviates from power law scaling in Minkowski and FRW simulations. This deviation is attributed to a winding of the electroweak group parameters around the domain walls in our simulations. We observe a local violation of the neutral vacuum condition on the domain walls in our simulations. This violation is attributed to relative electroweak transformations across the domain walls which is a general feature emerging from random initial conditions.

We discuss axion dark matter detection via two mechanisms: spontaneous decays and resonant conversion in neutron star magnetospheres. For decays, we show that the brightness temperature signal, rather than flux, is a less ambiguous measure for selecting candidate objects. This is owing principally to the finite beam width of telescopes which prevents one from being sensitive to the total flux from the object. With this in mind, we argue that the large surface-mass-density of the galactic centre or the Virgo cluster centre offers the best chance of improving current constraints on the axion-photon coupling via spontaneous decays. For the neutron star case, we first carry out a detailed study of mixing in magnetised plasmas. We derive transport equations for the axion-photon system via a controlled gradient expansion, allowing us to address inhomogeneous mass-shell constraints for arbitrary momenta. We then derive a non-perturbative Landau-Zener formula for the conversion probability valid across the range of relativistic and non-relativistic axions and show that the standard perturbative resonant conversion amplitude is a truncation of this result in the non-adiabatic limit. Our treatment reveals that that infalling dark matter axions typically convert non-adiabatically in magnetospheres. We describe the limitations of one-dimensional mixing equations and explain how three-dimensional effects activate new photon polarisations, including longitudinal modes and illustrate these arguments with numerical simulations in higher dimensions. We find that the bandwidth of the radio signal is dominated by Doppler broadening from the relative motion of the neutron star with respect to the observer. Therefore, we conclude that the radio signal from the resonant decay is weaker than previously thought, which means one relies on local density peaks to probe weaker axion-photon couplings.

Cosmological Birefringence (CB) is a phenomenon, caused by parity violating modifications to electrodynamics, whereby the linear polarisation angle of light changes as photons traverse a vacuum. It is possible to use a number of different analysis techniques to constrain this effect using Cosmic Microwave Background (CMB) polarisation observations. We investigate two different methods of constraining direction dependent birefringence for present and future CMB experiments including BICEP/Keck, Simons Observatory (SO), and LiteBIRD . Specifically we compare the constraints placed on anisotropic CB from a quadratic estimator technique to those derived from estimates of the $B$-mode power-spectrum for the three different experiments. The constraints derived from estimates of the $B$-mode power spectrum are found to be comparable to those derived from quadratic estimator for BICEP/Keck and SO, but not LiteBIRD due to its larger sky coverage. These forecasted upper bounds for CB are converted to constraints on primordial magnetic fields and the coupling between photons and pseudo Nambu-Goldstone bosons. Finally we show that even with the best constraints on CB, for the respective experiments, the potentially induced $B$-mode power can act as a significant contaminant in the prospective measurement of primordial $B$-modes.

We present optical and near-infrared imaging covering a $\sim$1.53 deg$^2$ region in the Super-Cluster Assisted Shear Survey (SuperCLASS) field, which aims to make the first robust weak lensing measurement at radio wavelengths. We derive photometric redshifts for $\approx$176,000 sources down to $i^\prime_{\rm AB}\sim24$ and present photometric redshifts for 1.4 GHz $e$-MERLIN and VLA detected radio sources found in the central 0.26 deg$^{2}$. We compile an initial catalog of 149 radio sources brighter than S$_{1.4}>75$ $\mu$Jy and find their photometric redshifts span $0<z_{\rm phot}<4$ with radio luminosities between $10^{21}-10^{25}$ W Hz$^{-1}$, with medians of $\langle z \rangle =0.55$ and $\langle L_{1.4}\rangle =1.9\times10^{23}$ W Hz$^{-1}$ respectively. We find 95% of the \uJy\ radio source sample (141/149) have SEDs best fit by star-forming templates while 5% (8/149) are better fit by AGN. Spectral indices are calculated for sources with radio observations from VLA and GMRT at 325 MHz, with an average spectral slope of $\alpha=0.59\pm0.04$. Using the full photometric redshift catalog we construct a density map at the redshift of the known galaxy clusters, $z=0.20\pm0.08$. Four of the five clusters are prominently detected at $>7 \sigma$ in the density map and we confirm the photometric redshifts are consistent with previously measured spectra from a few galaxies at the cluster centers.

We describe the first results on weak gravitational lensing from the SuperCLASS survey: the first survey specifically designed to measure the weak lensing effect in radio-wavelength data, both alone and in cross-correlation with optical data. We analyse 1.53 square degrees of optical data from the Subaru telescope and 0.26 square degrees of radio data from the e-MERLIN and VLA telescopes (the DR1 data set). Using standard methodologies on the optical data only we make a significant (10 sigma) detection of the weak lensing signal (a shear power spectrum) due to the massive supercluster of galaxies in the targeted region. For the radio data we develop a new method to measure the shapes of galaxies from the interferometric data, and we construct a simulation pipeline to validate this method. We then apply this analysis to our radio observations, treating the e-MERLIN and VLA data independently. We achieve source densities of 0.5 per square arcmin in the VLA data and 0.06 per square arcmin in the e-MERLIN data, numbers which prove too small to allow a detection of a weak lensing signal in either the radio data alone or in cross-correlation with the optical data. Finally, we show preliminary results from a visibility-plane combination of the data from e-MERLIN and VLA which will be used for the forthcoming full SuperCLASS data release. This approach to data combination is expected to enhance both the number density of weak lensing sources available and the fidelity with which their shapes can be measured.

The SuperCLuster Assisted Shear Survey (SuperCLASS) is a legacy programme using the e-MERLIN interferometric array. The aim is to observe the sky at L-band (1.4 GHz) to a r.m.s. of 7 uJy per beam over an area of ~1 square degree centred on the Abell 981 supercluster. The main scientific objectives of the project are: (i) to detect the effects of weak lensing in the radio in preparation for similar measurements with the Square Kilometre Array (SKA); (ii) an extinction free census of star formation and AGN activity out to z~1. In this paper we give an overview of the project including the science goals and multi-wavelength coverage before presenting the first data release. We have analysed around 400 hours of e-MERLIN data allowing us to create a Data Release 1 (DR1) mosaic of ~0.26 square degrees to the full depth. These observations have been supplemented with complementary radio observations from the Karl G. Jansky Very Large Array (VLA) and optical/near infra-red observations taken with the Subaru, Canada-France-Hawaii and Spitzer Telescopes. The main data product is a catalogue of 887 sources detected by the VLA, of which 395 are detected by e-MERLIN and 197 of these are resolved. We have investigated the size, flux and spectral index properties of these sources finding them compatible with previous studies. Preliminary photometric redshifts, and an assessment of galaxy shapes measured in the radio data, combined with a radio-optical cross-correlation technique probing cosmic shear in a supercluster environment, are presented in companion papers.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

We use the Equation of State (EoS) approach to study the evolution of the dark sector in Horndeski models, the most general scalar-tensor theories with second order equations of motion. By including the effects of the dark sector into our code EoS\_class, we demonstrate the numerical stability of the formalism and excellent agreement with results from other publicly available codes for a range of parameters describing the evolution of the function characterising the perturbations for Horndeski models, $\alpha_{\rm x}$, with ${\rm x}=\{{\rm K}, {\rm B}, {\rm M}, {\rm T}\}$. After demonstrating that on sub-horizon scales ($k\gtrsim 10^{-3}~{\rm Mpc}^{-1}$ at $z=0$) velocity perturbations in both the matter and the dark sector are typically subdominant with respect to density perturbations in the equation of state for perturbations, we find an attractor solution for the dark sector gauge-invariant density perturbation $\Delta_{\rm ds}$ by neglecting its time derivatives in the equation describing its time evolution, as commonly done in the well-known quasi-static approximation. Using this result, we provide simplified expressions for the equation-of-state functions: the dark sector entropy perturbations $w_{\rm ds}\Gamma_{\rm ds}$ and anisotropic stress $w_{\rm ds}\Pi_{\rm ds}$. From this we derive a growth factor-like equation for both matter and dark sector and are able to capture the relevant physics for several observables with great accuracy. We finally present new analytical expressions for the well-known modified gravity phenomenological functions $\mu$, $\eta$ and $\Sigma$ for a generic Horndeski model as functions of $\alpha_{\rm x}$. We show that on small scales they reproduce expressions presented in previous works, but on large scales, we find differences with respect to other works.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

We investigate the impact of $1/f$ noise on cosmology for an intensity mapping survey with SKA1-MID Band\,1 and Band\,2. We use a Fisher matrix approach to forecast constraints on cosmological parameters under the influence of $1/f$ noise, adopting a semi-empirical model from an earlier work, which results from the residual $1/f$ noise spectrum after applying a component separation algorithm to remove smooth spectral components. Without $1/f$ noise, the projected constraints are $4\%$ on $w_0$, $1\%$ on $h$, $2\%$ on $b_{\rm HI}$ using Band\,1+\emph{Planck}, and $3\%$ on $w_0$, $0.5\%$ on $h$, $2\%$ on $b_{\rm HI}$ using Band\,2+\emph{Planck}. A representative baseline $1/f$ noise degrades these constraints by a factor of $\sim1.5$ for Band\,1+\emph{Planck}, and $\sim1.2$ for Band\,2+\emph{Planck}. On the power spectrum measurement, higher redshift and smaller scales are more affected by $1/f$ noise, with minimal contamination comes from $z\lesssim1$ and $\ell\lesssim100$. Subject to the specific scan strategy of the adopted $1/f$ noise model, one prefers a correlated in frequency with minimised spectral slope, a low knee frequency, and a large telescope slew speed in order to reduce its impact.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

The Simons Observatory (SO) is a new cosmic microwave background experiment being built on Cerro Toco in Chile, due to begin observations in the early 2020s. We describe the scientific goals of the experiment, motivate the design, and forecast its performance. SO will measure the temperature and polarization anisotropy of the cosmic microwave background in six frequency bands: 27, 39, 93, 145, 225 and 280 GHz. The initial configuration of SO will have three small-aperture 0.5-m telescopes (SATs) and one large-aperture 6-m telescope (LAT), with a total of 60,000 cryogenic bolometers. Our key science goals are to characterize the primordial perturbations, measure the number of relativistic species and the mass of neutrinos, test for deviations from a cosmological constant, improve our understanding of galaxy evolution, and constrain the duration of reionization. The SATs will target the largest angular scales observable from Chile, mapping ~10% of the sky to a white noise level of 2 $\mu$K-arcmin in combined 93 and 145 GHz bands, to measure the primordial tensor-to-scalar ratio, $r$, at a target level of $\sigma(r)=0.003$. The LAT will map ~40% of the sky at arcminute angular resolution to an expected white noise level of 6 $\mu$K-arcmin in combined 93 and 145 GHz bands, overlapping with the majority of the LSST sky region and partially with DESI. With up to an order of magnitude lower polarization noise than maps from the Planck satellite, the high-resolution sky maps will constrain cosmological parameters derived from the damping tail, gravitational lensing of the microwave background, the primordial bispectrum, and the thermal and kinematic Sunyaev-Zel'dovich effects, and will aid in delensing the large-angle polarization signal to measure the tensor-to-scalar ratio. The survey will also provide a legacy catalog of 16,000 galaxy clusters and more than 20,000 extragalactic sources.

We investigate generalized Einstein-Aether theories that are compatible with the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation, and lensing data. For a given dark energy equation of state, $w_{\rm de}$, we formulate a designer approach and we investigate their impact on the CMB temperature anisotropy and matter power spectra. We use the Equation of State approach to parametrize the perturbations and find that this approach is particularly useful in identifying the most suitable and numerically efficient parameters to explore in a Markov chain Monte Carlo (MCMC) analysis. We find the data constrains models with $w_{\rm de} = -1$ to be compatible with $\Lambda$CDM. For $w_{\rm de} \not= -1$ models, which avoid the gravitational waves constraint through the entropy perturbation, we constrain $w_{\rm de}$ to be $w_{\rm de } = -1.06^{+0.08}_{-0.03}$ (CMB) and $w_{\rm de } =-1.04^{+0.05}_{-0.02}$ (CMB+Lensing) at $68\%$C.L., and find that these models can be different from $\Lambda$CDM and still be compatible with the data. We also find that these models can ameliorate some anomalies in $\Lambda$CDM when confronted with data, such as the low-$\ell$ and high-$k$ power in the CMB temperature anisotropy and matter power spectra respectively, but not simultaneously. We also investigate the anomalous lensing amplitude, quantified by $A_{\rm lens}$, and find that that for $w_{\rm de} = -1$ models, $A_{\rm lens} = 1.15^{+0.07}_{-0.08}$ (CMB) and $A_{\rm lens} = 1.12\pm0.05$ (CMB+Lensing) at $68\%$C.L. $\sim$ 2$\sigma$ larger than expected, similar to previous analyses of $\Lambda$CDM.

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.

We develop a new method of extracting simultaneous measurements of weak lensing shear and a local rotation of the plane of polarization using observations of resolved radio sources. We show that the direction of polarization is statistically linked with that of the gradient of the total intensity field, and this provides the basis of our method. Using a number of sources spread over the sky, this method allows constraints to be placed on cosmic shear and birefringence, and it can be applied to any resolved radio sources for which such a correlation exists. Assuming that the rotation and shear are constant across the source, we use this relationship to construct a quadratic estimator and investigate its properties using simulated observations. We develop a calibration scheme using simulations based on the observed images to mitigate a bias which occurs in the presence of measurement errors and an astrophysical scatter on the polarization. The method is applied directly to archival data of radio galaxies where we measure a mean rotation signal of $\left<\omega\right>=-2.03^{\circ}\pm0.75^{\circ}$ and an average shear compatible with zero using 30 reliable sources. This level of constraint on an overall rotation is comparable with current leading constraints from CMB experiments and is expected to increase by at least an order of magnitude with future high precision radio surveys, such as those performed by the SKA. We also measure the shear and rotation two-point correlation functions and estimate the number of sources required to detect shear and rotation correlations in future surveys.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

We explore the constraints on dark sector models imposed by the recent observation of coincident gravitational waves and gamma rays from a binary neutron star merger, GW170817. Rather than focusing on specific models as has been considered by other authors, we explore this in the context of the equation of state approach of which the specific models are special cases. After confirming the strong constraints found by others for Horndeski, Einstein-Aether and massive gravity models, we discuss how it is possible to construct models which might evade the constraints from GW170817 but still leading to cosmologically interesting modifications to gravity. Possible examples are ``miracle cancellations" such as in $f(R)$ models, nonlocal models and higher-order derivatives. The latter two rely on the dimensionless ratio of the wave number of the observed gravitational waves to the Hubble expansion rate being very large ($\sim10^{19}$) which is used to suppress modifications to the speed of gravitational waves.

We review the Equation of State (EoS) approach to dark sector perturbations and apply it to $f(\mathcal{R})$ gravity models of dark energy. We show that the EoS approach is numerically stable and use it to set observational constraints on designer models. Within the EoS approach we build an analytical understanding of the dynamics of cosmological perturbations for the designer class of $f(\mathcal{R})$ gravity models, characterised by the parameter $B_0$ and the background equation of state of dark energy $w$. When we use the Planck Cosmic Microwave Background (CMB) temperature anisotropy, polarisation and lensing data as well as the Baryonic Acoustic Oscillation (BAO) data from SDSS and WiggleZ, we find $B_0<0.006$ (95\%CL) for the designer models with $w=-1$. Furthermore, we find $B_0<0.0045$ and $|w+1|<0.002$ (95\%CL) for the designer models with $w\neq -1$. Previous analyses found similar results for designer and Hu-Sawicki $f(\mathcal{R})$ gravity models using the Effective Field Theory (EFT) approach (Raveri et al. 2014, Hu et al. 2016), therefore this hints for the fact that generic $f(\mathcal{R})$ models with $w\neq-1$ can be tightly constrained by current cosmological data, complementary to solar system tests (Brax et al. 2008, Faulkner et al. 2007). When compared to a $w$CDM fluid with the same sound speed, we find that the equation of state for $f(\mathcal{R})$ models is better constrained to be close to -1 by about an order of magnitude, due to the strong dependence of the perturbations on $w$.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

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 present a deep survey of the SuperCLASS super-cluster - a region of sky known to contain five Abell clusters at redshift $z\sim0.2$ - performed using the Arcminute Microkelvin Imager (AMI) Large Array (LA) at 15.5$~$GHz. Our survey covers an area of approximately 0.9 square degrees. We achieve a nominal sensitivity of $32.0~\mu$Jy beam$^{-1}$ toward the field centre, finding 80 sources above a $5\sigma$ threshold. We derive the radio colour-colour distribution for sources common to three surveys that cover the field and identify three sources with strongly curved spectra - a high-frequency-peaked source and two GHz-peaked-spectrum sources. The differential source count (i) agrees well with previous deep radio source count, (ii) exhibits no evidence of an emerging population of star-forming galaxies, down to a limit of 0.24$~$mJy, and (iii) disagrees with some models of the 15$~$GHz source population. However, our source count is in agreement with recent work that provides an analytical correction to the source count from the SKADS Simulated Sky, supporting the suggestion that this discrepancy is caused by an abundance of flat-spectrum galaxy cores as-yet not included in source population models.

Cosmology has entered an era where the experimental limitations are not due to instrumental sensitivity but instead due to inherent systematic uncertainties in the instrumentation and data analysis methods. The field of HI intensity mapping (IM) is still maturing, however early attempts are already systematics limited. One such systematic limitation is 1/f noise, which largely originates within the instrumentation and manifests as multiplicative gain fluctuations. To date there has been little discussion about the possible impact of 1/f noise on upcoming single-dish HI IM experiments such as BINGO, FAST or SKA. Presented in this work are Monte-Carlo end-to-end simulations of a 30 day HI IM survey using the SKA-MID array covering a bandwidth of 950 and 1410 MHz. These simulations extend 1/f noise models to include not just temporal fluctuations but also correlated gain fluctuations across the receiver bandpass. The power spectral density of the spectral gain fluctuations are modelled as a power-law, and characterised by a parameter $\beta$. It is found that the degree of 1/f noise frequency correlation will be critical to the success of HI IM experiments. Small values of $\beta$ ($\beta$ < 0.25) or high correlation is preferred as this is more easily removed using current component separation techniques. The spectral index of temporal fluctuations ($\alpha$) is also found to have a large impact on signal-to-noise. Telescope slew speed has a smaller impact, and a scan speed of 1 deg s$^{-1}$ should be sufficient for a HI IM survey with the SKA.

HI intensity mapping is a new observational technique to survey the large-scale structure of matter using the 21 cm emission line of atomic hydrogen (HI). In this work, we simulate BINGO (BAO from Integrated Neutral Gas Observations) and SKA (Square Kilometre Array) phase-1 dish array operating in auto-correlation mode. For the optimal case of BINGO with no foregrounds, the combination of the HI angular power spectra with Planck results allows $w$ to be measured with a precision of $4\%$, while the combination of the BAO acoustic scale with Planck gives a precision of $7\%$. We consider a number of potentially complicating effects, including foregrounds and redshift dependent bias, which increase the uncertainty on $w$ but not dramatically; in all cases the final uncertainty is found to be $\Delta w < 8\%$ for BINGO. For the combination of SKA-MID in auto-correlation mode with Planck, we find that, in ideal conditions, $w$ can be measured with a precision of $4\%$ for the redshift range $0.35 < z < 3$ (i.e., for the bandwidth of $\Delta \nu = [350, 1050]$ MHz) and $2\%$ for $0 < z < 0.49$ (i.e., $\Delta \nu = [950, 1421]$ MHz). Extending the model to include the sum of neutrino masses yields a $95\%$ upper limit of $\sum m_\nu < 0.24$ eV for BINGO and $\sum m_\nu < 0.08$ eV for SKA phase 1, competitive with the current best constraints in the case of BINGO and significantly better than them in the case of SKA.

We investigate the evolution of cosmological perturbations in models of dark energy described by a time-like unit normalized vector field specified by a general function $\mathcal{F}(\mathcal{K})$, so-called Generalized Einstein-Aether models. First we study the background dynamics of such models via a designer approach in an attempt to model this theory as dark energy. We find that only one specific form of this designer approach matches $\Lambda$CDM at background order and we also obtain a differential equation which $\mathcal{F}(\mathcal{K})$ must satisfy for general $w$CDM cosmologies. We also present the equations of state for perturbations in Generalized Einstein-Aether models, which completely parametrize these models at the level of linear perturbations. A generic feature of modified gravity models is that they introduce new degrees of freedom. By fully eliminating these we are able to express the gauge invariant entropy perturbation and the scalar, vector, and tensor anisotropic stresses in terms of the perturbed fluid variables and metric perturbations only. These can then be used to study the evolution of perturbations in the scalar, vector, and tensor sectors and we use these to evolve the Newtonian gravitational potentials.

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

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

Discordance in the $\Lambda$CDM cosmological model can be seen by comparing parameters constrained by CMB measurements to those inferred by probes of large scale structure. Recent improvements in observations, including final data releases from both Planck and SDSS-III BOSS, as well as improved astrophysical uncertainty analysis of CFHTLenS, allows for an update in the quantification of any tension between large and small scales. This paper is intended, primarily, as a discussion on the quantifications of discordance when comparing the parameter constraints of a model when given two different data sets. We consider KL-divergence, comparison of Bayesian evidences and other statistics which are sensitive to the mean, variance and shape of the distributions. However, as a by-product, we present an update to the similar analysis in (Battye, Charnock and Moss; 2015) where we find that, considering new data and treatment of priors, the constraints from the CMB and from a combination of LSS probes are in greater agreement and any tension only persists to a minor degree. In particular, we find the parameter constraints from the combination of LSS probes which are most discrepant with the Planck2015+Pol+BAO parameter distributions can be quantified at a 2.55$\sigma$ tension using the method introduced in (Battye, Charnock and Moss; 2015). If instead we use the distributions constrained by the combination of LSS probes which are in greatest agreement with those from Planck2015+Pol+BAO this tension is only 0.76$\sigma$.

21cm intensity mapping is a novel approach aimed at measuring the power spectrum of density fluctuations and deducing cosmological information, notably from the Baryonic Acoustic Oscillations (BAO). We give an update on the progress of BAO from Integrated Neutral Gas Observations (BINGO) which is a single dish intensity mapping project. First we explain the basic ideas behind intensity mapping concept before updating the instrument design for BINGO. We also outline the survey we plan to make and its projected science output including estimates of cosmological parameters.

We study the nature of potentials in scalar field based models for dark energy - with both canonical and noncanonical kinetic terms. We calculate numerically, and using an analytic approximation around $a\approx 1$, potentials for models with constant equation-of-state parameter, $w_{\phi}$. We find that for a wide range of models with canonical and noncanonical kinetic terms there is a simple approximation for the potential that holds when the scale factor is in the range $0.6\lesssim a\lesssim 1.4$. We discuss how this form of the potential can also be used to represent models with non-constant $w_{\phi}$ and, hence, how it could be used in reconstruction from cosmological data.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

We present the results of 325 MHz GMRT observations of a super-cluster field, known to contain five Abell clusters at redshift $z \sim 0.2$. We achieve a nominal sensitivity of $34\,\mu$Jy beam$^{-1}$ toward the phase centre. We compile a catalogue of 3257 sources with flux densities in the range $183\,\mu\rm{Jy}\,-\,1.5\,\rm{Jy}$ within the entire $\sim 6.5$ square degree field of view. Subsequently, we use available survey data at other frequencies to derive the spectral index distribution for a sub-sample of these sources, recovering two distinct populations -- a dominant population which exhibit spectral index trends typical of steep-spectrum synchrotron emission, and a smaller population of sources with typically flat or rising spectra. We identify a number of sources with ultra-steep spectra or rising spectra for further analysis, finding two candidate high-redshift radio galaxies and three gigahertz-peaked-spectrum radio sources. Finally, we derive the Euclidean-normalised differential source counts using the catalogue compiled in this work, for sources with flux densities in excess of $223 \, \mu$Jy. Our differential source counts are consistent with both previous observations at this frequency and models of the low-frequency source population. These represent the deepest source counts yet derived at 325 MHz. Our source counts exhibit the well-known flattening at mJy flux densities, consistent with an emerging population of star-forming galaxies; we also find marginal evidence of a downturn at flux densities below $308 \, \mu$Jy, a feature so far only seen at 1.4 GHz.

It has been proposed recently that the SKA1-MID could be used to conduct an HI intensity mapping survey that could rival upcoming Stage IV dark energy surveys. However, measuring the weak HI signal is expected to be very challenging due to contaminations such as residual Galactic emission, RFI, and instrumental 1/f noise. Modelling the effects of these contaminants on the cosmological HI signal requires numerical end-to-end simulations. Here we present how 1/f noise within the receiver can double the effective uncertainty of an SKA-like survey to HI on large angular scales (l < 50).

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts.

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $\tau$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $\tau$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $\tau$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $\Lambda$CDM cosmological model and recent models of reionization.

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

The B-mode polarisation power spectrum in the Cosmic Microwave Background (CMB) is about four orders of magnitude fainter than the CMB temperature power spectrum. Any instrumental imperfections that couple temperature fluctuations to B-mode polarisation must therefore be carefully controlled and/or removed. We investigate the role that a scan strategy can have in mitigating certain common systematics by averaging systematic errors down with many crossing angles. We present approximate analytic forms for the error on the recovered B-mode power spectrum that would result from differential gain, differential pointing and differential ellipticity for the case where two detector pairs are used in a polarisation experiment. We use these analytic predictions to search the parameter space of common satellite scan strategies in order to identify those features of a scan strategy that have most impact in mitigating systematic effects. As an example we go on to identify a scan strategy suitable for the CMB satellite proposed for the ESA M5 call. considering the practical considerations of fuel requirement, data rate and the relative orientation of the telescope to the earth. Having chosen a scan strategy we then go on to investigate the suitability of the scan strategy.

We discuss the detectability of large-scale HI intensity fluctuations using the FAST telescope. We present forecasts for the accuracy of measuring the Baryonic Acoustic Oscillations and constraining the properties of dark energy. The FAST $19$-beam L-band receivers ($1.05$--$1.45$ GHz) can provide constraints on the matter power spectrum and dark energy equation of state parameters ($w_{0},w_{a}$) that are comparable to the BINGO and CHIME experiments. For one year of integration time we find that the optimal survey area is $6000\,{\rm deg}^2$. However, observing with larger frequency coverage at higher redshift ($0.95$--$1.35$ GHz) improves the projected errorbars on the HI power spectrum by more than $2~\sigma$ confidence level. The combined constraints from FAST, CHIME, BINGO and Planck CMB observations can provide reliable, stringent constraints on the dark energy equation of state.

The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,\mu$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,\mu$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,\sigma$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy.

We study the equations for the evolution of cosmological perturbations in $f\left(\mathcal{R}\right)$ and conclude that this modified gravity model can be expressed as a dark energy fluid at background and linearised perturbation order. By eliminating the extra scalar degree of freedom known to be present in such theories, we are able to characterise the evolution of the perturbations in the scalar sector in terms of equations of state for the entropy perturbation and anisotropic stress which are written in terms of the density and velocity perturbations of the dark energy fluid and those in the matter, or the metric perturbations. We also do the same in the much simpler vector and tensor sectors. In order to illustrate the simplicity of this formulation, we numerically evolve perturbations in a small number of cases.

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

With the temperature power spectrum of the cosmic microwave background (CMB) at least four orders of magnitude larger than the B-mode polarisation power spectrum, any instrumental imperfections that couple temperature to polarisation must be carefully controlled and/or removed. Here we present two new map-making algorithms that can create polarisation maps that are clean of temperature-to-polarisation leakage systematics due to differential gain and pointing between a detector pair. Where a half wave plate is used, we show that the spin-2 systematic due to differential ellipticity can also by removed using our algorithms. The algorithms require no prior knowledge of the imperfections or temperature sky to remove the temperature leakage. Instead, they calculate the systematic and polarisation maps in one step directly from the time ordered data (TOD). The first algorithm is designed to work with scan strategies that have a good range of crossing angles for each map pixel and the second for scan strategies that have a limited range of crossing angles. The first algorithm can also be used to identify if systematic errors that have a particular spin are present in a TOD. We demonstrate the use of both algorithms and the ability to identify systematics with simulations of TOD with realistic scan strategies and instrumental noise.

HI intensity mapping is an emerging tool to probe dark energy. Observations of the redshifted HI signal will be contaminated by instrumental noise, atmospheric and Galactic foregrounds. The latter is expected to be four orders of magnitude brighter than the HI emission we wish to detect. We present a simulation of single-dish observations including an instrumental noise model with 1/f and white noise, and sky emission with a diffuse Galactic foreground and HI emission. We consider two foreground cleaning methods: spectral parametric fitting and principal component analysis. For a smooth frequency spectrum of the foreground and instrumental effects, we find that the parametric fitting method provides residuals that are still contaminated by foreground and 1/f noise, but the principal component analysis can remove this contamination down to the thermal noise level. This method is robust for a range of different models of foreground and noise, and so constitutes a promising way to recover the HI signal from the data. However, it induces a leakage of the cosmological signal into the subtracted foreground of around 5%. The efficiency of the component separation methods depends heavily on the smoothness of the frequency spectrum of the foreground and the 1/f noise. We find that as, long as the spectral variations over the band are slow compared to the channel width, the foreground cleaning method still works.

We present current and future constraints on equations of state for dark sector perturbations. The equations of state considered are those corresponding to a generalized scalar field model and time-diffeomorphism invariant $L(g)$ theories that are equivalent to models of a relativistic elastic medium and also Lorentz violating massive gravity. We develop a theoretical understanding of the observable impact of these models. In order to constrain these models we use CMB temperature data from Planck, BAO measurements, CMB lensing data from Planck and the South Pole Telescope, and weak galaxy lensing data from CFHTLenS. We find non-trivial exclusions on the range of parameters, although the data remains compatible with $w=-1$. We gauge how future experiments will help to constrain the parameters. This is done via a likelihood analysis for CMB experiments such as CoRE and PRISM, and tomographic galaxy weak lensing surveys, focussing in on the potential discriminatory power of Euclid on mildly non-linear scales.

We develop and apply the position angle-only shear estimator of Whittaker et al. (2014) to realistic galaxy images. This is done by demonstrating the method on the simulations of the GREAT3 challenge (Mandelbaum et al. 2014b), which include contributions from anisotropic PSFs. We measure the position angles of the galaxies using three distinct methods - the integrated light method, quadrupole moments of surface brightness, and using model-based ellipticity measurements provided by IM3SHAPE. A weighting scheme is adopted to address biases in the position angle measurements which arise in the presence of an anisotropic PSF. Biases on the shear estimates, due to measurement errors on the position angles and correlations between the measurement errors and the true position angles, are corrected for using simulated galaxy images and an iterative procedure. The properties of the simulations are estimated using the deep field images provided as part of the challenge. A method is developed to match the distributions of galaxy fluxes and half-light radii from the deep fields to the corresponding distributions in the field of interest. We recover angle-only shear estimates with a performance close to current well-established model and moments-based methods for all three angle measurement techniques. The Q-values for all three methods are found to be Q~400. The code is freely available online at http://www.jb.man.ac.uk/~mbrown/angle_only_shear/.

Asymmetric beams can create significant bias in estimates of the power spectra from CMB experiments. With the temperature power spectrum many orders of magnitude stronger than the B-mode power spectrum any systematic error that couples the two must be carefully controlled and/or removed. Here, we derive unbiased estimators for the CMB temperature and polarisation power spectra taking into account general beams and general scan strategies. A simple consequence of asymmetric beams is that, even with an ideal scan strategy where every sky pixel is seen at every orientation, there will be residual coupling from temperature power to B-mode power if the orientation of the beam asymmetry is not aligned with the orientation of the co-polarisation. We test our correction algorithm on simulations of two temperature-only experiments and demonstrate that it is unbiased. The simulated experiments use realistic scan strategies, noise levels and highly asymmetric beams. We also develop a map-making algorithm that is capable of removing beam asymmetry bias at the map level. We demonstrate its implementation using simulations and show that it is capable of accurately correcting both temperature and polarisation maps for all of the effects of beam asymmetry including the effects of temperature to polarisation leakage.

We discuss methods for performing weak lensing using radio observations to recover information about the intrinsic structural properties of the source galaxies. Radio surveys provide unique information that can benefit weak lensing studies, such as HI emission, which may be used to construct galaxy velocity maps, and polarized synchrotron radiation; both of which provide information about the unlensed galaxy and can be used to reduce galaxy shape noise and the contribution of intrinsic alignments. Using a proxy for the intrinsic position angle of an observed galaxy, we develop techniques for cleanly separating weak gravitational lensing signals from intrinsic alignment contamination in forthcoming radio surveys. Random errors on the intrinsic orientation estimates introduce biases into the shear and intrinsic alignment estimates. However, we show that these biases can be corrected for if the error distribution is accurately known. We demonstrate our methods using simulations, where we reconstruct the shear and intrinsic alignment auto and cross-power spectra in three overlapping redshift bins. We find that the intrinsic position angle information can be used to successfully reconstruct both the lensing and intrinsic alignment power spectra with negligible residual bias.

Any Dark Energy (DE) or Modified Gravity (MG) model that deviates from a cosmological constant requires a consistent treatment of its perturbations, which can be described in terms of an effective entropy perturbation and an anisotropic stress. We have considered a recently proposed generic parameterisation of DE/MG perturbations and compared it to data from the Planck satellite and six galaxy catalogues, including temperature-galaxy (Tg), CMB lensing-galaxy and galaxy-galaxy (gg) correlations. Combining these observables of structure formation with tests of the background expansion allows us to investigate the properties of DE/MG both at the background and the perturbative level. Our constraints on DE/MG are mostly in agreement with the cosmological constant paradigm, while we also find that the constraint on the equation of state w (assumed to be constant) depends on the model assumed for the perturbation evolution. We obtain $w=-0.92^{+0.20}_{-0.16}$ (95% CL; CMB+gg+Tg) in the entropy perturbation scenario; in the anisotropic stress case the result is $w=-0.86^{+0.17}_{-0.16}$. Including the lensing correlations shifts the results towards higher values of w. If we include a prior on the expansion history from recent Baryon Acoustic Oscillations (BAO) measurements, we find that the constraints tighten closely around $w=-1$, making it impossible to measure any DE/MG perturbation evolution parameters. If, however, upcoming observations from surveys like DES, Euclid or LSST show indications for a deviation from a cosmological constant, our formalism will be a useful tool towards model selection in the dark sector.

We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

We develop a method for performing a weak lensing analysis using only measurements of galaxy position angles. By analysing the statistical properties of the galaxy orientations given a known intrinsic ellipticity distribution, we show that it is possible to obtain estimates of the shear by minimizing a $\chi^2$ statistic. The method is demonstrated using simulations where the components of the intrinsic ellipticity are taken to be Gaussian distributed. Uncertainties in the position angle measurements introduce a bias into the shear estimates which can be reduced to negligible levels by introducing a correction term into the formalism. We generalize our approach by developing an algorithm to obtain direct shear estimators given any azimuthally symmetric intrinsic ellipticity distribution. We introduce a method of measuring the position angles of the galaxies from noisy pixelized images, and propose a method to correct for biases which arise due to pixelization and correlations between measurement errors and galaxy ellipticities. We also develop a method to constrain the sample of galaxies used to obtain an estimate of the intrinsic ellipticity distribution such that fractional biases in the resulting shear estimates are below a given threshold value. We demonstrate the angle only method by applying it to simulations where the ellipticities are taken to follow a log-normal distribution. We compare the performance of the position angle only method with the standard method based on full ellipticity measurements by reconstructing lensing convergence maps from both numerical simulations and from the CFHTLenS data. We find that the difference between the convergence maps reconstructed using the two methods is consistent with noise.

There is a tension between measurements of the amplitude of the power spectrum of density perturbations inferred using the Cosmic Microwave Background (CMB) and directly measured by Large-Scale Structure (LSS) on smaller scales. We show that this tension exists, and is robust, for a range of LSS indicators including clusters, lensing and redshift space distortions and using CMB data from either $Planck$ or WMAP+SPT/ACT. One obvious way to try to reconcile this is the inclusion of a massive neutrino which could be either active or sterile. Using $Planck$ and a combination of all the LSS data we find that (i) for an active neutrino $\sum m_{\nu}= (0.357\pm0.099)\,{\rm eV}$ and (ii) for a sterile neutrino $m_{\rm sterile}^{\rm eff}= (0.67\pm0.18)\,{\rm eV}$ and $\Delta N_{\rm eff}= 0.32\pm0.20$. This is, however, at the expense of a degraded fit to $Planck$ temperature data, and we quantify the residual tension at $2.5\sigma$ and $1.6 \sigma$ for massive and sterile neutrinos respectively. We also consider alternative explanations including a lower redshift for reionization that would be in conflict with polarisation measurements made by WMAP and $ad$-$hoc$ modifications to primordial power spectrum.

We consider a novel null test for contamination which can be applied to CMB polarization data that involves analysis of the statistics of the polarization position angles. Specifically, we will concentrate on using histograms of the measured position angles to illustrate the idea. Such a test has been used to identify systematics in the NVSS point source catalogue with an amplitude well below the noise level. We explore the statistical properties of polarization angles in CMB maps. If the polarization angle is not correlated between pixels, then the errors follow a simple $\sqrt{N_{pix}}$ law. However this is typically not the case for CMB maps since these have correlations which result in an increase in the variance since the effective number of independent pixels is reduced. Then we illustrate how certain classes of systematic errors can result in very obvious patterns in these histograms, and thus that these errors could possibly be identified using this method. We discuss how this idea might be applied in a realistic context, and make a preliminary analysis of the WMAP7 data, finding evidence of a systematic error in the Q and W band data, consistent with a constant offset in Q and U.

The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

The evolution of perturbations is a crucial part of the phenomenology of the dark sector cosmology. We advocate parameterizing these perturbations using equations of state for the entropy perturbation and the anisotropic stress. For small perturbations, these equations of state will be linear in the density, velocity and metric perturbations, and in principle these can be related back to the field content of the underlying model allowing for confrontation with observations. We illustrate our point by constructing gauge invariant entropy perturbations for theories where the dark sector Lagrangian is a general function of a scalar field, its first and second derivatives, and the metric and its first derivative, ${\cal L}={\cal L}(\phi,\partial_\mu\phi,\partial_\mu\partial_\nu\phi,g_{\mu\nu},\partial_{\alpha}g_{\mu\nu})$. As an example, we show how to apply this approach to the case of models of Kinetic Gravity Braiding.

We present a methodology for computing model independent perturbations in dark energy and modified gravity. This is done from the Lagrangian for perturbations, by showing how field content, symmetries, and physical principles are often sufficient ingredients for closing the set of perturbed fluid equations. The fluid equations close once "equations of state for perturbations" are identified: these are linear combinations of fluid and metric perturbations which construct gauge invariant entropy and anisotropic stress perturbations for broad classes of theories. Our main results are the proof of the equation of state for perturbations presented in a previous paper, and the development of the required calculational tools.

We present constraints on cosmological parameters using number counts as a function of redshift for a sub-sample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the Sunyaev--Zeldovich (SZ) effect, and the sub-sample used here has a signal-to-noise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass $M$ and SZ signal $Y$ calibrated to X-ray measurements, we derive constraints on the power spectrum amplitude $\sigma_8$ and matter density parameter $\Omega_{\mathrm{m}}$ in a flat $\Lambda$CDM model. We test the robustness of our estimates and find that possible biases in the $Y$--$M$ relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the X-ray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that $\sigma_8=0.75\pm 0.03$, $\Omega_{\mathrm{m}}=0.29\pm 0.02$, and $\sigma_8(\Omega_{\mathrm{m}}/0.27)^{0.3} = 0.764 \pm 0.025$. The value of $\sigma_8$ is degenerate with the mass bias; if the latter is fixed to a value of 20% we find $\sigma_8(\Omega_{\mathrm{m}}/0.27)^{0.3}=0.78\pm 0.01$ and a tighter one-dimensional range $\sigma_8=0.77\pm 0.02$. We find that the larger values of $\sigma_8$ and $\Omega_{\mathrm{m}}$ preferred by Planck's measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the $\Lambda$CDM model, such as a component of massive neutrinos (abridged).

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

The accumulation of redshifts provides a significant observational bottleneck when using galaxy cluster surveys to constrain cosmological parameters. We propose a simple method to allow the use of samples where there is a fraction of the redshifts that are not known. The simplest assumption is that the missing redshifts are randomly extracted from the catalogue, but the method also allows one to take into account known selection effects in the accumulation of redshifts. We quantify the reduction in statistical precision of cosmological parameter constraints as a function of the fraction of missing redshifts for simulated surveys, and also investigate the impact of making an incorrect assumption for the distribution of missing redshifts.

The QUIJOTE (Q-U-I JOint Tenerife) CMB Experiment is designed to observe the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range of 10-40 GHz. The first of the two QUIJOTE telescopes and the multi-frequency (10-20 GHz) instrument have been in operation since November 2012. In 2014 a second telescope and a new instrument at 30GHz will be ready for commissioning, and an additional instrument at 40 GHz is in its final design stages. After three years of effective observations, the data obtained by these telescopes and instruments will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r = 0.05. At the moment, we have completed half of the wide Galactic survey with the multi-frequency instrument covering 18 000 square degrees of the Northern hemisphere. When we finish this survey in early 2014, we shall have reached approximately 14{\mu}K per one degree beam at 11, 13, 17 and 19 GHz, in both Q and U.

The two fundamental assumptions of the standard cosmological model - that the initial fluctuations are statistically isotropic and Gaussian - are rigorously tested using maps of the cosmic microwave background (CMB) anisotropy from the Planck satellite. Deviations from isotropy have been found and demonstrated to be robust against component separation algorithm, mask choice and frequency dependence. Many of these anomalies were previously observed in the WMAP data, and are now confirmed at similar levels of significance (about 3 sigma). However, we find little evidence for non-Gaussianity, with the exception of a few statistical signatures that seem to be associated with specific anomalies. In particular, we find that the quadrupole-octopole alignment is also connected to a low observed variance of the CMB signal. A power asymmetry is now found to persist to scales corresponding to about l=600, and can be described in the low-l regime by a phenomenological dipole modulation model. However, any primordial power asymmetry is strongly scale-dependent and does not extend to arbitrarily small angular scales. Finally, it is plausible that some of these features may be reflected in the angular power spectrum of the data, which shows a deficit of power on similar scales. Indeed, when the power spectra of two hemispheres defined by a preferred direction are considered separately, one shows evidence for a deficit in power, while its opposite contains oscillations between odd and even modes that may be related to the parity violation and phase correlations also detected in the data. Although these analyses represent a step forward in building an understanding of the anomalies, a satisfactory explanation based on physically motivated models is still lacking.

Recent results by the Planck collaboration have shown that cosmological parameters derived from the cosmic microwave background anisotropies and cluster number counts are in tension, with the latter preferring lower values of the matter density parameter, $\Omega_\mathrm{m}$, and power spectrum amplitude, $\sigma_8$. Motivated by this, we investigate the extent to which the tension may be ameliorated once the effect of baryonic depletion on the cluster mass function is taken into account. We use the large-volume Millennium Gas simulations in our study, including one where the gas is pre-heated at high redshift and one where the gas is heated by stars and active galactic nuclei (in the latter, the self-gravity of the baryons and radiative cooling are omitted). In both cases, the cluster baryon fractions are in reasonably good agreement with the data at low redshift, showing significant depletion of baryons with respect to the cosmic mean. As a result, it is found that the cluster abundance in these simulations is around 15 per cent lower than the commonly-adopted fit to dark matter simulations by Tinker et al (2008) for the mass range $10^{14}-10^{14.5}h^{-1} \mathrm{M}_\odot$. Ignoring this effect produces a significant artificial shift in cosmological parameters which can be expressed as $\Delta[\sigma_8(\Omega_\mathrm{m}/0.27)^{0.38}]\simeq -0.03$ at $z=0.17$ (the median redshift of the $\mathit{Planck}$ cluster sample) for the feedback model. While this shift is not sufficient to fully explain the $\mathit{Planck}$ discrepancy, it is clear that such an effect cannot be ignored in future precision measurements of cosmological parameters with clusters. Finally, we outline a simple, model-independent procedure that attempts to correct for the effect of baryonic depletion and show that it works if the baryon-dark matter back-reaction is negligible.

We discuss whether massive neutrinos (either active or sterile) can reconcile some of the tensions within cosmological data that have been brought into focus by the recently released {\it Planck} data. We point out that a discrepancy is present when comparing the primary CMB and lensing measurements both from the CMB and galaxy lensing data using CFHTLenS, similar to that which arises when comparing CMB measurements and SZ cluster counts. A consistent picture emerges and including a prior for the cluster constraints and BAOs we find that: for an active neutrino model with 3 degenerate neutrinos, $\sum m_{\nu}= (0.320 \pm 0.081)\,{\rm eV}$, whereas for a sterile neutrino, in addition to 3 neutrinos with a standard hierarchy and $\sum m_{\nu}= 0.06\,{\rm eV}$, $m_{\nu, \, \rm sterile}^{\rm eff}= (0.450 \pm 0.124)\,{\rm eV}$ and $\Delta N_{\rm eff} = 0.45 \pm 0.23$. In both cases there is a significant detection of modification to the neutrino sector from the standard model and in the case of the sterile neutrino it is possible to reconcile the BAO and local $H_0$ measurements. However, a caveat to our result is some internal tension between the CMB and lensing/cluster observations, and the masses are in excess of those estimated from the shape of the matter power spectrum from galaxy surveys.

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.

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.

Planck data have been used to provide stringent new constraints on cosmic strings and other defects. We describe forecasts of the CMB power spectrum induced by cosmic strings, calculating these from network models and simulations using line-of-sight Boltzmann solvers. We have studied Nambu-Goto cosmic strings, as well as field theory strings for which radiative effects are important, thus spanning the range of theoretical uncertainty in strings models. We have added the angular power spectrum from strings to that for a simple adiabatic model, with the extra fraction defined as $f_{10}$ at multipole $\ell=10$. This parameter has been added to the standard six parameter fit using COSMOMC with flat priors. For the Nambu-Goto string model, we have obtained a constraint on the string tension of $G\mu/c^2 < 1.5 x 10^{-7}$ and $f_{10} < 0.015$ at 95% confidence that can be improved to $G\mu/c^2 < 1.3 x 10^{-7}$ and $f_{10} < 0.010$ on inclusion of high-$\ell$ CMB data. For the abelian-Higgs field theory model we find, $G\mu_{AH}/c^2 < 3.2 x 10^{-7}$ and $f_{10} < 0.028$. The marginalized likelihoods for $f_{10}$ and in the $f_{10}$--$\Omega_b h^2$ plane are also presented. We have also obtained constraints on $f_{10}$ for models with semi-local strings and global textures for which $G\mu/c^2 < 1.1 x 10^{-6}$. We have made complementarity searches for the specific non-Gaussian signatures of cosmic strings, calibrating with all-sky Planck resolution CMB maps generated from networks of post-recombination strings. We have obtained upper limits on the string tension at 95% confidence of $G\mu/c^2 < 8.8 x 10^{-7}$ using modal bispectrum estimation and $G\mu/c^2 < 7.8 x 10^{-7}$ for real space searches with Minkowski functionals. These are conservative upper bounds because only post-recombination string contributions have been included in the non-Gaussian analysis.

We present projected constraints on the cosmic string tension, $G\mu/c^2$, that could be achieved by future gravitational wave detection experiments and express our results as semi-analytic relations of the form $G\mu(\Omega_{\rm gw}h^2)/c^2$, to allow for direct computation of the tension constraints for future experiments. These results can be applied to new constraints on $\ogwh$ as they are imposed. Experiments operating in different frequency bands probe different parts of the gravitational wave spectrum of a cosmic string network and are sensitive to different uncertainties in the underlying cosmic string model parameters. We compute the gravitational wave spectra of cosmic string networks based on the one-scale model, covering all the parameter space accessed by each experiment which is strongly dependent on the birth scale of loops relative to the horizon, $\alpha$. The upper limits on the string tension avoid any assumptions on the model parameters. We perform this investigation for Pulsar Timing Array experiments of different durations as well as ground-based and space-borne interferometric detectors.

We consider a class of modified gravity models where the terms added to the standard Einstein-Hilbert Lagrangian are just a function of the metric only. For linearized perturbations around an isotropic space-time, this class of models is entirely specified by a rank 4 tensor that encodes possibly time-dependent masses for the gravitons. This tensor has the same symmetries as an elasticity tensor, suggesting an interpretation of massive gravity as an effective rigidity of space-time. If we choose a form for this tensor which is compatible with the symmetries of FRW and enforce full reparameterization invariance, then the only theory possible is a cosmological constant. However, in the case where the theory is only time translation invariant, the ghost-free massive gravity theory is equivalent to the elastic dark energy scenario with the extra Lorentz violating vector giving rise to 2 transverse and 1 longitudinal degrees of freedom, whereas when one demands spatial translation invariance one is left with scalar field theory with a non-standard kinetic term.

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y for a sample of 19 objects for which weak lensing (WL) mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from XMM-Newton archive data and the SZ effect signal is measured from Planck all-sky survey data. We find an M_WL-D_A^2 Y relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our SZ effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods, and, for the present sample, the mass discrepancy and difference in mass concentration is especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations.

BINGO is a concept for performing a 21cm intensity mapping survey using a single dish telescope. We briefly discuss the idea of intensity mapping and go on to define our single dish concept. This involves a \sim 40 m dish with an array of \sim 50 feed horns placed \sim 90 m above the dish using a pseudo-correlation detection system based on room temperature LNAs and one of the celestial poles as references. We discuss how such an array operating between 960 and 1260 MHz could be used to measure the acoustic scale to 2.4% over the redshift range 0.13<z<0.48 in around 1 year of on-source integration time by performing a 10 deg \times 200 deg drift scan survey with a resolution of \sim 2/3 deg.

We discuss the detection of large scale HI intensity fluctuations using a single dish approach with the ultimate objective of measuring the Baryonic Acoustic Oscillations and constraining the properties of dark energy. We present 3D power spectra, 2D angular power spectra for individual redshift slices, and also individual line-of-sight spectra, computed using the S^3 simulated HI catalogue which is based on the Millennium Simulation. We consider optimal instrument design and survey strategies for a single dish observation at low and high redshift for a fixed sensitivity. For a survey corresponding to an instrument with T_sys=50 K, 50 feed horns and 1 year of observations, we find that at low redshift (z \approx 0.3), a resolution of 40 arc min and a survey of 5000 deg^2 is close to optimal, whereas at higher redshift (z \approx 0.9) a resolution of 10 arcmin and 500 deg^2 would be necessary. Continuum foreground emission from the Galaxy and extragalactic radio sources are potentially a problem. We suggest that it could be that the dominant extragalactic foreground comes from the clustering of very weak sources. We assess its amplitude and discuss ways by which it might be mitigated. We then introduce our concept for a single dish telescope designed to detect BAO at low redshifts. It involves an under-illumintated static 40 m dish and a 60 element receiver array held 90 m above the under-illuminated dish. Correlation receivers will be used with each main science beam referenced against an antenna pointing at one of the Celestial Poles for stability and control of systematics. We make sensitivity estimates for our proposed system and projections for the uncertainties on the power spectrum after 1 year of observations. We find that it is possible to measure the acoustic scale at z\approx 0.3 with an accuracy 2.4% and that w can be measured to an accuracy of 16%.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

In light of upcoming observations modelling perturbations in dark energy and modified gravity models has become an important topic of research. We develop an effective action to construct the components of the perturbed dark energy momentum tensor which appears in the perturbed generalized gravitational field equations, {\delta}G_{\mu\nu} = 8{\pi}G{\delta}T_{\mu\nu} + {\delta}U_{\mu\nu} for linearized perturbations. Our method does not require knowledge of the Lagrangian density of the dark sector to be provided, only its field content. The method is based on the fact that it is only necessary to specify the perturbed Lagrangian to quadratic order and couples this with the assumption of global statistical isotropy of spatial sections to show that the model can be specified completely in terms of a finite number of background dependent functions. We present our formalism in a coordinate independent fashion and provide explicit formulae for the perturbed conservation equation and the components of {\delta}U_{\mu\nu} for two explicit generic examples: (i) the dark sector does not contain extra fields, L = L(g_{\mu\nu}) and (ii) the dark sector contains a scalar field and its first derivative L = L(g_{\mu\nu}, {\phi}, \nabla_{\mu}{\phi}). We discuss how the formalism can be applied to modified gravity models containing derivatives of the metric, curvature tensors, higher derivatives of the scalar fields and vector fields.

We investigate the constraints that can be placed on the cosmic string tension by using the current Pulsar Timing Array limits on the stochastic gravitational wave background (SGWB). We have developed a code to compute the spectrum of gravitational waves (GWs) based on the widely accepted one-scale model. In its simplest form the one-scale model allows one to vary: (i) the string tension, G\mu/c^2; (ii) the size of cosmic string loops relative to the horizon at birth, \alpha; (iii) the spectral index of the emission spectrum, q; (iv) the cut-off in the emission spectrum, n_*; and (v) the intercommutation probability, p. The amplitude and slope of the spectrum in the nHz frequency range is very sensitive to these unknown parameters. We have also investigated the impact of more complicated scenarios with multiple initial loop sizes, in particular the 2-\alpha models proposed in the literature and a log-normal distribution for \alpha. We have computed the constraint on G\mu/c^2 due to the limit on a SGWB imposed by data from the European Pulsar Timing Array. Taking into account all the possible uncertainties in the parameters we find a conservative upper limit of G\mu/c^2<5.3x 10^{-7} which typically occurs when the loop production scale is close to the gravitational backreaction scale, \alpha\approx\Gamma G\mu/c^2. Stronger limits are possible for specific values of the parameters which typically correspond to the extremal cases \alpha\ll \Gamma G\mu/c^2 and \alpha\gg \Gamma G\mu/c^2. This limit is less stringent than the previously published limits which are based on cusp emission, an approach which does not necessarily model all the possible uncertainties. We discuss the prospects for lowering this limit by two orders of magnitude, or even a detection of the SGWB, in the very near future in the context of the Large European Array for Pulsars and the Square Kilometre Array.

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

We have exploited the large-volume Millennium Gas cosmological N-body hydrodynamics simulations to study the SZ cluster population at low and high redshift, for three models with varying gas physics. We confirm previous results using smaller samples that the intrinsic (spherical) Y_{500}-M_{500} relation has very little scatter (sigma_{log_{10}Y}~0.04), is insensitive to cluster gas physics and evolves to redshift one in accord with self-similar expectations. Our pre-heating and feedback models predict scaling relations that are in excellent agreement with the recent analysis from combined Planck and XMM-Newton data by the Planck Collaboration. This agreement is largely preserved when r_{500} and M_{500} are derived using the hydrostatic mass proxy, Y_{X,500}, albeit with significantly reduced scatter (sigma_{log_{10}Y}~0.02), a result that is due to the tight correlation between Y_{500} and Y_{X,500}. Interestingly, this assumption also hides any bias in the relation due to dynamical activity. We also assess the importance of projection effects from large-scale structure along the line-of-sight, by extracting cluster Y_{500} values from fifty simulated 5x5 square degree sky maps. Once the (model-dependent) mean signal is subtracted from the maps we find that the integrated SZ signal is unbiased with respect to the underlying clusters, although the scatter in the (cylindrical) Y_{500}-M_{500} relation increases in the pre-heating case, where a significant amount of energy was injected into the intergalactic medium at high redshift. Finally, we study the hot gas pressure profiles to investigate the origin of the SZ signal and find that the largest contribution comes from radii close to r_{500} in all cases. The profiles themselves are well described by generalised Navarro, Frenk & White profiles but there is significant cluster-to-cluster scatter.

We develop an analytic formalism that allows one to quantify the stability properties of X-type and Y-type junctions in domain wall networks in two dimensions. A similar approach might be applicable to more general defect systems involving junctions that appear in a range of physical situations, for example, in the context of F- and D-type strings in string theory. We apply this formalism to a particular field theory, Carter's pentavac model, where the strength of the symmetry breaking is governed by the parameter $|\epsilon|< 1$. We find that for low values of the symmetry breaking parameter X-type junctions will be stable, whereas for higher values an X-type junction will separate into two Y-type junctions. The critical angle separating the two regimes is given by \alpha_c = 293^{\circ}\sqrt{|\epsilon|}$ and this is confirmed using simple numerical experiments. We go on to simulate the pentavac model from random initial conditions and we find that the dominant junction is of \ytype for $|\epsilon| \geq 0.02$ and is of \xtype for $|\epsilon| \leq 0.02$. We also find that for small $\epsilon$ the evolution of the number of domain walls $\qsubrm{N}{dw}$ in Minkowski space does not follow the standard $\propto t^{-1}$ scaling law with the deviation from the standard lore being more pronounced as $\epsilon$ is decreased. The presence of dissipation appears to restore the $t^{-1}$ lore.

The characterisation of the atomic and molecular hydrogen content of high-redshift galaxies is a major observational challenge that will be addressed over the coming years with a new generation of radio telescopes. We investigate this important issue by considering the states of hydrogen across a range of structures within high-resolution cosmological hydrodynamical simulations. Additionally, our simulations allow us to investigate the sensitivity of our results to numerical resolution and to sub-grid baryonic physics (especially feedback from supernovae and active galactic nuclei). We find that the most significant uncertainty in modelling the neutral hydrogen distribution arises from our need to model a self-shielding correction in moderate density regions. Future simulations incorporating radiative transfer schemes will be vital to improve on our empirical self-shielding threshold. Irrespective of the exact nature of the threshold we find that while the atomic hydrogen mass function evolves only mildly from redshift two to zero, the molecular hydrogen mass function increases with increasing redshift, especially at the high-mass end. Interestingly, the weak evolution of the neutral hydrogen mass function is insensitive to the feedback scheme utilised, but the opposite is true for the molecular gas, which is more closely associated with the star formation in the simulations.

We perform a systematic study of generic accidental Higgs-family and CP symmetries that could occur in the two-Higgs-doublet-model potential, based on a Majorana scalar-field formalism which realizes a subgroup of GL(8,C). We derive the general conditions of convexity and stability of the scalar potential and present analytical solutions for two non-zero neutral vacuum expectation values of the Higgs doublets for a typical set of six symmetries, in terms of the gauge-invariant parameters of the theory. By means of a homotopy-group analysis, we identify the topological defects associated with the spontaneous symmetry breaking of each symmetry, as well as the massless Goldstone bosons emerging from the breaking of the continuous symmetries. We find the existence of domain walls from the breaking of Z_2, CP1 and CP2 discrete symmetries, vortices in models with broken U(1)_PQ and CP3 symmetries and a global monopole in the SO(3)_HF-broken model. The spatial profile of the topological defect solutions is studied in detail, as functions of the potential parameters of the two-Higgs doublet model. The application of our Majorana scalar-field formalism in studying more general scalar potentials that are not constrained by the U(1)_Y hypercharge symmetry is discussed. In particular, the same formalism may be used to properly identify seven additional symmetries that may take place in a U(1)_Y-invariant scalar potential.

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

We present 30 GHz Sunyaev Zel'dovich observations of a statistically complete sample of galaxy clusters with OCRA-p. The clusters are the 18 most X-ray luminous clusters at z > 0.2 in the ROSAT Brightest Cluster Sample. We correct for contaminant radio sources via supplementary observations with the Green Bank Telescope, also at 30 GHz, and remove a cluster that is contaminated by an unresolved X-ray source. All 17 remaining clusters have central SZ effects with Comptonisation parameter y_0 exceeding 1.9x10^-4, and 13 are detected at significance > 3 sigma. We use our data to examine scalings between y_0 and X-ray temperature, X-ray luminosity, and the X-ray mass proxy Y_X, and find good agreement with predictions from self-similar models of cluster formation, with an intrinsic scatter in y_0 of about 25%. We also comment on the success of the observations in the face of the contaminant source population, and the implications for upcoming cm-wave surveys.

Knowledge of the population of radio sources in the range ~2-200 GHz is important for understanding their effects on measurements of the Cosmic Microwave Background power spectrum. We report measurements of the 30 GHz flux densities of 605 radio sources from the Combined Radio All-sky Targeted Eight-GHz Survey (CRATES), which have been made with the One Centimetre Receiver Array prototype (OCRA-p) on the Torun 32-m telescope. The flux densities of sources that were also observed by WMAP and previous OCRA surveys are in broad agreement with those reported here, however a number of sources display intrinsic variability. We find a good correlation between the 30 GHz and Fermi gamma-ray flux densities for common sources. We examine the radio spectra of all observed sources and report a number of Gigahertz-peaked and inverted spectrum sources. These measurements will be useful for comparison to those from the Low Frequency Instrument of the Planck satellite, which will make some of its most sensitive observations in the region covered here.

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.

We study $F(R)$ modified gravity models which are capable of driving the accelerating epoch of the Universe at the present time whilst not destroying the standard Big Bang and inflationary cosmology. Recent studies have shown that a weak curvature singularity with $|R|\to\infty$ can arise generically in viable $F(R)$ models of present dark energy (DE) signaling an internal incompleteness of these models. In this work we study how this problem is cured by adding a quadratic correction with a sufficiently small coefficient to the $F(R)$ function at large curvatures. At the same time, this correction eliminates two more serious problems of previously constructed viable $F(R)$ DE models: unboundedness of the mass of a scalar particle (scalaron) arising in $F(R)$ gravity and the scalaron overabundance problem. Such carefully constructed models can also yield both an early time inflationary epoch and a late time de Sitter phase with vastly different values of $R$. The reheating epoch in these combined models of primordial and present dark energy is completely different from that of the old $R + R^{2}/6M^{2}$ inflationary model, mainly due to the fact that values of the effective gravitational constant at low and intermediate curvatures are different for positive and negative $R$. This changes the number of e-folds during the observable part of inflation that results in a different value of the primordial power spectrum index.

In a recent paper (Brown & Battye 2011), we proposed the use of integrated polarization measurements of background galaxies in radio weak gravitational lensing surveys and investigated the potential impact on the statistical measurement of cosmic shear. Here we extend this idea to reconstruct maps of the projected dark matter distribution, or lensing convergence field. The addition of polarization can, in principle, greatly reduce shape noise due to the intrinsic dispersion in galaxy ellipticities. We show that maps reconstructed using this technique in the radio band can be competitive with those derived using standard lensing techniques which make use of many more galaxies. In addition, since the reconstruction noise is uncorrelated between these standard techniques and the polarization technique, their comparison can serve as a powerful check for systematics and their combination can reduce noise further. We examine the convergence reconstruction which could be achieved with two forthcoming facilities: (i) a deep survey, covering 1.75 square degrees using the e-MERLIN instrument currently being commissioned in the UK and (ii) the high resolution, deep wide field surveys which will eventually be conducted with the Square Kilometre Array.

We have studied the implications of high sensitivity polarization measurements of objects from the WMAP point source catalogue made using the VLA at 8.4, 22 and 43 GHz. The fractional polarization of sources is almost independent of frequency with a median of ~2 per cent and an average, for detected sources, of ~3.5 per cent. These values are also independent of the total intensity over the narrow range of intensity we sample. Using a contemporaneous sample of 105 sources detected at all 3 VLA frequencies, we have investigated the spectral behaviour as a function of frequency by means of a 2-colour diagram. Most sources have power-law spectra in total intensity, as expected. On the other hand they appear to be almost randomly distributed in the polarized intensity 2-colour diagram. This is compatible with the polarized spectra being much less smooth than those in intensity and we speculate on the physical origins of this. We have performed an analysis of the correlations between the fractional polarization and spectral indices including computation of the principal components. We find that there is little correlation between the fractional polarization and the intensity spectral indices. This is also the case when we include polarization measurements at 1.4 GHz from the NVSS. In addition we compute 45 rotation measures from polarization position angles which are compatible with a \lambda^2 law. We use our results to predict the level of point source confusion noise that contaminates CMB polarization measurements aimed at detecting primordial gravitational waves from inflation. We conclude that some level of source subtraction will be necessary to detect r~0.1 below 100 GHz and at all frequencies to detect r~0.01. We present estimates of the level of contamination expected and the number of sources which need to be subtracted as a function of the imposed cut flux density and frequency.

Physical constants and cosmological parameters could vary with position. On the largest scales such variations would manifest themselves as gradients across our Hubble volume, leading to dipole-modulation of the cosmic microwave anisotropies. This generically leads to a correlation between adjacent multipoles in the spherical harmonics expansion of the sky, a distinctive signal which should be searched for in future data sets.

It has been shown that superconducting domain walls in a model with U(1) x Z2 symmetry can form long-lived loops called kinky vortons from random initial conditions in the broken field and a uniform charged background in (2+1) dimensions. In this paper we investigate a similar model with a hyper-cubic symmetry coupled to an unbroken U(1) in which the domain walls can form junctions and hence a lattice. We call this model the charge-coupled cubic-anisotropy (CCCA) model. First, we present a detailed parametric study of the U(1) x Z2 model; features which we vary include the nature of the initial conditions and the coupling constants. This allows us to identify interesting parameters to vary in the more complicated, and hence more computationally intensive, CCCA models. In particular we find that the coefficient of the interaction term can be used to engineer three separate regimes: phase mixing, condensation and phase separation with the condensation regime corresponding to a single value of the coupling constant defined by the determinant of the quartic interaction terms being zero. We then identify the condensation regime in the CCCA model and show that, in this regime, the number of domain walls does not scale in the standard way if the initial conditions have a sufficiently high background charge. Instead of forming loops of domain wall, we find that, within the constraints of dynamic range, the network appears to be moving toward a glass-like configuration. We find that the results are independent of the dimension of the hyper-cube.

We propose a new technique for weak gravitational lensing in the radio band making use of polarization information. Since the orientation of a galaxy's polarized emission is both unaffected by lensing and is related to the galaxy's intrinsic orientation, it effectively provides information on the unlensed galaxy position angle. We derive a new weak lensing estimator which exploits this effect and makes full use of both the observed galaxy shapes and the estimates of the intrinsic position angles as provided by polarization. Our method has the potential to both reduce the effects of shot noise, and to reduce to negligible levels, in a model-independent way, all effects of intrinsic galaxy alignments. We test our technique on simulated weak lensing skies, including an intrinsic alignment contaminant consistent with recent observations, in three overlapping redshift bins. Adopting a standard weak lensing analysis and ignoring intrinsic alignments results in biases of 5-10% in the recovered power spectra and cosmological parameters. Applying our new estimator to one tenth the number of galaxies used for the standard case, we recover both power spectra and the input cosmology with similar precision as compared to the standard case and with negligible residual bias, even in the presence of a substantial (astrophysical) scatter in the relationship between the observed orientation of the polarized emission and the intrinsic orientation. Assuming a reasonable polarization fraction for star-forming galaxies, and no cosmological conspiracy in the relationship between polarization direction and intrinsic morphology, our estimator should prove a valuable tool for weak lensing analyses of forthcoming radio surveys, in particular, deep wide field surveys with e-MERLIN, MeerKAT and ASKAP and ultimately, definitive radio lensing surveys with the SKA.

We re-examine the constraints on the cosmic string tension from Cosmic Microwave Background (CMB) and matter power spectra, and also from limits on a stochastic background of gravitational waves provided by pulsar timing. We discuss the different approaches to modeling string evolution and radiation. In particular, we show that the unconnected segment model can describe CMB spectra expected from thin string (Nambu) and field theory (Abelian-Higgs) simulations using the computed values for the correlation length, rms string velocity and small-scale structure relevant to each variety of simulation. Applying the computed spectra in a fit to CMB and SDSS data we find that $G\mu/c^2< 2.6\times 10^{-7}$ ($2 \sigma$) if the Nambu simulations are correct and $G\mu /c^2< 6.4\times 10^{-7}$ in the Abelian-Higgs case. The degeneracy between $G\mu/c^2$ and the power spectrum slope $n_{\rm S}$ is substantially reduced from previous work. Inclusion of constraints on the baryon density from Big Bang Nucleosynthesis (BBN) imply that $n_{\rm S} <1$ at around the $4\sigma$ level for both the Nambu and Abelian-Higgs cases. As a by-product of our results, we find there is "moderate-to-strong" Bayesian evidence that the Harrison-Zel'dovich spectrum is excluded (odds ratio of $\sim 100:1$) by the combination of CMB, SDSS and BBN when compared to the standard 6 parameter fit. Using the contribution to the gravitational wave background from radiation era loops as a conservative lower bound on the signal for specific values of $G\mu/c^2$ and loop production size, $\alpha$, we find that $G\mu /c^2< 7\times 10^{-7} $ for $\alpha c^2/(\Gamma G\mu)\ll1$ and $G\mu/c^2 < 5\times 10^{-11}/\alpha$ for $\alpha c^2/(\Gamma G\mu) \gg1$.

The back-reaction of baryons on the dark matter halo density profile is of great interest, not least because it is an important systematic uncertainty when attempting to detect the dark matter. Here, we draw on a large suite of high resolution cosmological hydrodynamical simulations, to systematically investigate this process and its dependence on the baryonic physics associated with galaxy formation. The inclusion of baryons results in significantly more concentrated density profiles if radiative cooling is efficient and feedback is weak. The dark matter halo concentration can in that case increase by as much as 30 (10) per cent on galaxy (cluster) scales. The most significant effects occur in galaxies at high redshift, where there is a strong anti-correlation between the baryon fraction in the halo centre and the inner slope of both the total and the dark matter density profiles. If feedback is weak, isothermal inner profiles form, in agreement with observations of massive, early-type galaxies. However, we find that AGN feedback, or extremely efficient feedback from massive stars, is necessary to match observed stellar fractions in groups and clusters, as well as to keep the maximum circular velocity similar to the virial velocity as observed for disk galaxies. These strong feedback models reduce the baryon fraction in galaxies by a factor of 3 relative to the case with no feedback. The AGN is even capable of reducing the baryon fraction by a factor of 2 in the inner region of group and cluster haloes. This in turn results in inner density profiles which are typically shallower than isothermal and the halo concentrations tend to be lower than in the absence of baryons.

We have investigated the correlations between optical and radio isophotal position angles for 14302 SDSS galaxies with r magnitudes brighter than 18. All the galaxies are identified with extended FIRST radio sources. For passive early-type galaxies, which we distinguish from the others by using the colour, concentration and their principal components, we find a strong statistical alignment of the radio axes with the optical minor axes. Since the radio emission is driven by accretion on to a nuclear black hole we argue that the observed correlation gives new insight into the structure of elliptical galaxies, for example, whether or not the nuclear kinematics are decoupled from the rest of the galaxy. Our results imply that a significant fraction of the galaxies are oblate spheroids, perhaps rotationally supported, with their radio emission aligned with the stellar minor axis. Remarkably, the strength of the correlation of the radio major axis with the optical minor axis depends on radio loudness. Those objects with a low ratio of FIRST radio flux density to total stellar light show a strong minor axis correlation while the stronger radio sources do not. This split may reflect different formation histories and we suggest this may be a new manifestation of the better known dichotomy between slow rotating and fast rotating ellipticals.

We use present cosmological data from the cosmic microwave background, large-scale structure and deuterium at high redshifts to constrain supersymmetric F- and D-term hybrid inflation scenarios including possible contributions to the CMB anisotropies from cosmic strings. Using two different realizations of the cosmic string spectrum, we find that the minimal version of the D-term model is ruled out at high significance. F-term models are also in tension with the data. We also discuss possible non-minimal variants of the models.

We present polarization measurements at 8.4, 22, and 43 GHz made with the VLA of a complete sample of extragalactic sources stronger than 1 Jy in the 5-year WMAP catalogue and with declinations north of -34 degrees. The observations were motivated by the need to know the polarization properties of radio sources at frequencies of tens of GHz in order to subtract polarized foregrounds for future sensitive Cosmic Microwave Background (CMB) experiments. The total intensity and polarization measurements are generally consistent with comparable VLA calibration measurements for less-variable sources, and within a similar range to WMAP fluxes for unresolved sources. A further paper will present correlations between measured parameters and derive implications for CMB measurements.

If the equation of state of dark energy is anisotropic there will be additional quadrupole anisotropy in the cosmic microwave background induced by the time dependent anisotropic stress quantified in terms of $\Delta w$. Assuming that the entire amplitude of the observed quadrupole is due to this anisotropy, we conservatively impose a limit of $|\Delta w| < 2.1\times 10^{-4}$ for any value of $w\ge -1$ assuming that $\Omega_{\rm m}<0.5$. This is considerably tighter than that which comes from SNe. Stronger limits, upto a factor of 10, are possible for specific values of $\Omega_{\rm m}$ and $w$. Since we assume this component is uncorrelated with the stochastic component from inflation, we find that both the expectation value and the sample variance are increased. There no improvement in the likelihood of an anomalously low quadrupole as suggested by previous work on an elliptical universe.

Small angular scale (high l) studies of cosmic microwave background anisotropies require accurate knowledge of the statistical properties of extragalactic sources at cm-mm wavelengths. We have used a 30 GHz dual-beam receiver (OCRA-p) on the Torun 32-m telescope to measure the flux densities of 121 sources in VSA fields selected at 15 GHz with the Ryle Telescope. We have detected 57 sources above a limiting flux density of 5mJy, of which 31 sources have a flux density greater than 10mJy, which is our effective completeness limit. From these measurements we derive a surface density of sources above 10mJy at 30 GHz of 2.0+/-0.4 per square degree. This is consistent with the surface density obtained by Mason et al. (2009) who observed a large sample of sources selected at a much lower frequency (1.4 GHz). We have also investigated the dependence of the spectral index distribution on flux density by comparing our results with those for sources above 1 Jy selected from the WMAP 22 GHz catalogue. We conclude that the proportion of steep spectrum sources increases with decreasing flux density, qualitatively consistent with the predictions of deZotti et al. (2005). We find no evidence for an unexpected population of sources whose spectra rise towards high frequencies, which would affect our ability to interpret current high resolution CMB observations at 30 GHz and above.

We investigate the correlations between optical and radio isophotal position angles for 14302 SDSS galaxies with $r$ magnitudes brighter than 18 and which have been associated with extended FIRST radio sources. We identify two separate populations of galaxies using the colour, concentration and their principal components. Surprisingly strong statistical alignments are found: late-type galaxies are overwhelmingly biased towards a position angle differences of $0^{\circ}$ and early-type galaxies to $90^{\circ}$. The late-type alignment can be easily understood in terms of the standard picture in which the radio emission is intimately related to areas of recent star-formation. In early-type galaxies the radio emission is expected to be driven by accretion on to a nuclear black hole. We argue that the observed correlation of the radio axis with the minor axis of the large-scale stellar distribution gives a fundamental insight into the structure of elliptical galaxies, for example, whether or not the nuclear kinematics are decoupled form the rest of the galaxy. Our results imply that the galaxies are oblate spheroids with their radio emission aligned with the minor axis. Remarkably the strength of the correlation of the radio major axis with the optical minor axis depends on radio loudness. Those objects with a low ratio of FIRST radio flux density to total stellar light show a strong minor axis correlation while the stronger radio sources do not. This may reflect different formation histories for the different objects and we suggest we may be seeing the different behaviour of rationally supported and non-rotationally supported ellipticals.

Using large numbers of simulations of the microwave sky, incorporating the Cosmic Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect due to clusters, we investigate the statistics of the power spectrum at microwave frequencies between spherical multipoles of 1000 and 10000. From these virtual sky maps, we find that the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than would be expected from purely Gaussian realizations, and has a distribution that is significantly skewed towards higher values, especially when small map sizes are used. The standard deviation is also increased by around 10 percent compared to the trispectrum calculation due to the clustering of galaxy clusters. We also consider the effects of including residual point sources and uncertainties in the gas physics. This has implications for the excess power measured in the CMB power spectrum by the Cosmic Background Imager and BIMA experiments. Our results indicate that the observed excess could be explained using a lower value of $\sigma_8$ than previously suggested, however the effect is not enough to match $\sigma_8=0.825$. The uncertainties in the gas physics could also play a substantial role. We have made our maps of the SZ effect available online.

We investigate the breaking of global statistical isotropy caused by a dark energy component with an energy-momentum tensor which has point symmetry, that could represent a cubic or hexagonal crystalline lattice. In such models Gaussian, adiabatic initial conditions created during inflation can lead to anisotropies in the cosmic microwave background whose spherical harmonic coefficients are correlated, contrary to the standard assumption. We develop an adaptation of the line of sight integration method that can be applied to models where the background energy-momentum tensor is isotropic, but whose linearized perturbations are anisotropic. We then show how this can be applied to the cases of cubic and hexagonal symmetry. We compute quantities which show that such models are indistinguishable from isotropic models even in the most extreme parameter choices, in stark contrast to models with anisotropic initial conditions based on inflation. The reason for this is that the dark energy based models contribute to the CMB anistropy via the inegrated Sachs-Wolfe effect, which is only relevent when the dark energy is dominant, that is, on the very largest scales. For inflationary models, however, the anisotropy is present on all scales.

We have examined the statistics of the polarization position angles determined for point sources in the NRAO-VLA sky survey (NVSS) and find that there is a statistically significant bias toward angles which are multiples of 45 degrees. The formal probability that the polarization angles are drawn from a uniform distribution is exponentially small. When the sample of those NVSS sources with polarizations detected with a signal to noise $\geq$3 is split either around the median polarized flux density or the median fractional polarization, the effect appears to be stronger for the more highly polarized sources. Regions containing strong sources and regions at low galactic latitudes are not responsible for the non-uniform distribution of position angles. We identify CLEAN bias as the probable cause of the dominant effect, coupled with small multiplicative and additive offsets on each of the Stokes parameters. Our findings have implications for the extraction of science, such as information concerning galactic magnetic fields, from large scale polarization surveys.

We present the current status of the QUIJOTE (Q-U-I JOint TEnerife) CMB Experiment, a new instrument which will start operations early 2009 at Teide Observatory, with the aim of characterizing the polarization of the CMB and other processes of galactic and extragalactic emission in the frequency range 10-30 GHz and at large angular scales. QUIJOTE will be a valuable complement at low frequencies for the PLANCK mission, and will have the required sensitivity to detect a primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

Global topological defects may account for the large cold spot observed in the Cosmic Microwave Background. We explore possibilities of constructing models of supersymmetric F-term hybrid inflation, where the waterfall fields are globally SU(2)-symmetric. In contrast to the case where SU(2) is gauged, there arise Goldstone bosons and additional moduli, which are lifted only by masses of soft-supersymmetry breaking scale. The model predicts the existence of global textures, which can become semi-local strings if the waterfall fields are gauged under U(1)_X. Gravitino overproduction can be avoided if reheating proceeds via the light SU(2)-modes or right-handed sneutrinos. For values of the inflaton- waterfall coupling >=10^-4, the symmetry breaking scale imposed by normalisation of the power spectrum generated from inflation coincides with the energy scale required to explain the most prominent of the cold spots. In this case, the spectrum of density fluctuations is close to scale-invariant which can be reconciled with measurements of the power spectrum by the inclusion of the sub-dominant component due to the topological defects.

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.

We study cosmological expansion in F(R) gravity using the trace of the field equations. High frequency oscillations in the Ricci scalar, whose amplitude increase as one evolves backward in time, have been predicted in recent works. We show that the approximations used to derive this result very quickly breakdown in any realistic model due to the non-linear nature of the underlying problem. Using a combination of numerical and semi-analytic techniques, we study a range of models which are otherwise devoid of known pathologies. We find that high frequency asymmetric oscillations and a singularity at finite time appear to be present for a wide range of initial conditions. We show that this singularity can be avoided with a certain range of initial conditions, which we find by evolving the models forwards in time. In addition we show that the oscillations in the Ricci scalar are highly suppressed in the Hubble parameter and scale factor.

We discuss the possibility of performing a substantial spectroscopic galaxy redshift survey selected via the 21cm emission from neutral hydrogen using the Five-hundred metre Aperture Spherical Telescope (FAST) to be built in China. We consider issues related to the estimation of the source counts and optimizations of the survey, and discuss the constraints on cosmological models that such a survey could provide. We find that a survey taking around two years could detect ~10^7 galaxies with an average redshift of ~0.15 making the survey complementary to those already carried out at optical wavelengths. These conservative estimates have used the z=0 HI mass function and have ignored the possibility of evolution. The results could be used to constrain Gamma = (Omega_m h) to 5 per cent and the spectral index, n_s, to 7 per cent independent of cosmic microwave background data. If we also use simulated power spectra from the Planck satellite, we can constrain w to be within 5 per cent of -1.

We present interferometric imaging at 33 GHz, with the new superextended configuration of the Very Small Array (VSA), of a very deep decrement in the cosmic microwave background (CMB) temperature. This decrement is located in the direction of the Corona Borealis supercluster, at a position with no known galaxy clusters, and was discovered by a previous VSA survey (Genova-Santos et al.). A total area of 3 sq.deg. has now been imaged, with an angular resolution of 7 arcmin and a flux sensitivity of 5 mJy/beam. These observations confirm the presence of this strong and resolved negative spot at -41+/-5 mJy/beam (-258+/-29 muK), with a signal to noise level of 8. This structure is also present in the WMAP 5-year data. The temperature of the W-band (94 GHz) data at the position of the decrement agrees within 1.2-sigma with that observed by the VSA at 33 GHz, and within 0.2-sigma with the Sunyaev-Zel'dovich (SZ) spectrum. Our analyses show that it is a non-Gaussian feature in the CMB at a level of 4.8-sigma. The probability of finding such a deviation or larger in CMB Gaussian simulations is only 0.19 per cent. Therefore, an explanation other than primordial CMB is required. We have considered the possibility of an SZ effect generated in a diffuse, extended warm/hot gas distribution. This hypothesis is especially relevant, as the presence of such structures, if confirmed, could provide the location for a significant fraction of the missing baryons in the Local Universe.

We present an analysis of the spectral properties of the extragalactic radio sources in the nearly-complete VSA sample at 33GHz. Data from different surveys are used to study source spectra between 1.4 and 33GHz. We find that, in general, spectra can not be well described by a single power law in the range of frequencies considered. In particular, most of the VSA sources that are steep between 1.4 and 5GHz, show a spectral flattening at \nu>5GHz. We identify 20 objects (19% of the sample) clearly characterized by an upturn spectrum, i.e., a spectrum falling at low frequencies (\nu<5GHz) and inverted at higher frequencies. Spectra with high-frequency flattening or upturn shape are supposed to occur when the emission from the AGN compact core begins to dominate over the component from extended lobes. This picture fits well with the AGN unified scheme, for objects observed at intermediate viewing angles of the AGN jet. Finally, we discuss implications that this class of sources can have on future CMB observations at high resolution.

By considering simple, but representative, models of brane inflation from a single brane-antibrane pair in the slow roll regime, we provide constraints on the parameters of the theory imposed by measurements of the CMB anisotropies by WMAP including a cosmic string component. We find that inclusion of the string component is critical in constraining parameters. In the most general model studied, which includes an inflaton mass term, as well as the brane-antibrane attraction, values n_s < 1.02 are compatible with the data at 95 % confidence level. We are also able to constrain the volume of internal manifold (modulo factors dependent on the warp factor) and the value of the inflaton field to be less than 0.66M_P at horizon exit. We also investigate models with a mass term. These observational considerations suggest that such models have r < 2*10^-5, which can only be circumvented in the fast roll regime, or by increasing the number of antibranes. Such a value of r would not be detectable in CMB polarization experiment likely in the near future, but the B-mode signal from the cosmic strings could be detectable. We present forecasts of what a similar analysis using PLANCK data would yield and find that it should be possible to rule out G\mu > 6.5*10^-8 using just the TT, TE and EE power spectra.

Modified gravity models are subject to a number of consistency requirements which restrict the form that the function $F(R)$ can take. We study a particular class of $F(R)$ functions which satisfy various constraints that have been found in the literature. These models have a late time accelerating epoch, and an acceptable matter era. We calculate the Friedmann equation for our models, and show that in order to satisfy the constraints we impose, they must mimic General Relativity plus $\Lambda$ throughout the cosmic history, with exponentially suppressed corrections. We also find that the free parameters in our model must be fine tuned to obtain an acceptable late time accelerating phase. We discuss the generality of this conclusion.

We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio sources to test the hypothesis that there is a systematic alignment of polarization position angle vectors on cosmological scales of the type claimed by Hutsemekers et al (2005). The polarization position angles of 4290 sources with polarized flux density >=1 mJy have been examined. They do not reveal large-scale alignments either as a whole or when split in half into high-redshift (z >= 1.24) and low-redshift sub-samples. Nor do the radio sources which lie in the specific areas covered by Hutsemekers et al (2005). show any significant effect. We have also looked at the position angles of parsec-scale jets derived from VLBI observations and again find no evidence for systematic alignments. Finally, we have investigated the correlation between the polarization position angle and those of the parsec-scale jets. As expected, we find that there is a tendency for the polarization angles to be perpendicular to the jet angles. However, the difference in jet and polarization position angles does not show any systematic trend in different parts of the sky.