Search

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., 45 pages, 23 figures

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., Comment: 15 pages, 8 figures, submitted to MNRAS

Axions are well-motivated candidates for dark matter. Recently, much interest has focused on the detection of photons produced by the resonant conversion of axion dark matter in neutron star magnetospheres. Various groups have begun to obtain radio data to search for the signal, however, more work is needed to obtain a robust theory prediction for the corresponding radio lines. In this work we derive detailed properties for the signal, obtaining both the line shape and time-dependence. The principal physical effects are from refraction in the plasma as well as from gravitation which together lead to substantial lensing which varies over the pulse period. The time-dependence from the co-rotation of the plasma with the pulsar distorts the frequencies leading to a Doppler broadened signal whose width varies in time. For our predictions, we trace curvilinear rays to the line of sight using the full set of equations from Hamiltonian optics for a dispersive medium in curved spacetime. Thus, for the first time, we describe the detailed shape of the line signal as well as its time dependence, which is more pronounced compared to earlier results. Our prediction of the features of the signal will be essential for this kind of dark matter search., 16 pages, 10 figures plus table: updated journal version, including extra appendix

Matched filters are routinely used in cosmology in order to detect galaxy clusters from mm observations through their thermal Sunyaev–Zeldovich (tSZ) signature. In addition, they naturally provide an observable, the detection signal-to-noise or significance, which can be used as a mass proxy in number counts analyses of tSZ-selected cluster samples. In this work, we show that this observable is, in general, non-Gaussian, and that it suffers from a positive bias, which we refer to as optimization bias. Both aspects arise from the fact that the signal-to-noise is constructed through an optimization operation on noisy data, and hold even if the cluster signal is modelled perfectly well, no foregrounds are present, and the noise is Gaussian. After reviewing the general mathematical formalism underlying matched filters, we study the statistics of the signal-to-noise with a set Monte Carlo mock observations, finding it to be well-described by a unit-variance Gaussian for signal-to-noise values of 6 and above, and quantify the magnitude of the optimization bias, for which we give an approximate expression that may be used in practice. We also consider the impact of the bias on the cluster number counts of Planck and the Simons Observatory (SO), finding it to be negligible for the former and potentially significant for the latter.

Author(s): Aghanim, N; Akrami, Y; Ashdown, M; Aumont, J; Baccigalupi, C; Ballardini, M; Banday, AJ; Barreiro, RB; Bartolo, N; Basak, S; Battye, R; Benabed, K; Bernard, JP; Bersanelli, M; Bielewicz, P; Bock, JJ; Bond, JR; Borrill, J; Bouchet, FR; Boulanger, F; Bucher, M; Burigana, C; Butler, RC; Calabrese, E; Cardoso, JF; Carron, J; Challinor, A; Chiang, HC; Chluba, J; Colombo, LPL; Combet, C; Contreras, D; Crill, BP; Cuttaia, F; De Bernardis, P; De Zotti, G; Delabrouille, J; Delouis, JM; DI Valentino, E; DIego, JM; Dore, O; Douspis, M; Ducout, A; Dupac, X; Dusini, S; Efstathiou, G; Elsner, F; Enslin, TA; Eriksen, HK; Fantaye, Y; Farhang, M; Fergusson, J; Fernandez-Cobos, R; Finelli, F; Forastieri, F; Frailis, M; Fraisse, AA; Franceschi, E; Frolov, A; Galeotta, S; Galli, S; Ganga, K; Genova-Santos, RT; Gerbino, M; Ghosh, T; Gonzalez-Nuevo, J; Gorski, KM; Gratton, S; Gruppuso, A; Gudmundsson, JE; Hamann, J; Handley, W; Hansen, FK; Herranz, D; Hildebrandt, SR; Hivon, E; Huang, Z; Jaffe, AH; Jones, WC; Karakci, A; Keihanen, E; Keskitalo, R; Kiiveri, K; Kim, J; Kisner, TS | Abstract: In the original version, the bounds given in Eqs. (87a) and (87b) on the contribution to the early-time optical depth, (15,30), contained a numerical error in deriving the 95th percentile from the Monte Carlo samples. The corrected 95% upper bounds are: τ(15,30) l 0:018 (lowE, flat τ(15, 30), FlexKnot), (1) τ(15, 30) l 0:023 (lowE, flat knot, FlexKnot): (2) These bounds are a factor of 3 larger than the originally reported results. Consequently, the new bounds do not significantly improve upon previous results from Planck data presented in Millea a Bouchet (2018) as was stated, but are instead comparable. Equations (1) and (2) give results that are now similar to those of Heinrich a Hu (2021), who used the same Planck 2018 data to derive a 95% upper bound of 0.020 using the principal component analysis (PCA) model and uniform priors on the PCA mode amplitudes.

A quasi-static approximation (QSA) for modified gravity can be applied in a number of ways. We consider three different analytical formulations based on applying this approximation to: (1) the field equations; (2) the equations for the two metric potentials; (3) the use of the attractor solution derived within the Equation of State (EoS) approach. We assess the veracity of these implementations on the effective gravitational constant ($\mu$) and the slip parameter ($\eta$), within the framework of Horndeski models. In particular, for a set of models we compare cosmological observables, i.e., the matter power spectrum and the CMB temperature and lensing angular power spectra, computed using the QSA, with exact numerical solutions. To do that, we use a newly developed branch of the CLASS code: QSA_class. All three approaches agree exactly on very small scales. Typically, we find that, except for $f(R)$ models where all the three approaches lead to the same result, the quasi-static approximations differ from the numerical calculations on large scales ($k \lesssim 3 - 4 \times 10^{-3}\,h\,{\rm Mpc}^{-1}$). Cosmological observables are reproduced to within 1% up to scales ${\rm K} = k/H_0$ of the order of a few and $\ell>5$ for the approaches based on the field equations and on the EoS, and we also do not find any appreciable difference if we use the scale-dependent expressions for $\mu$ and $\eta$ with respect to the value on small scales, showing that the formalism and the conclusions are reliable and robust, fixing the range of applicability of the formalism. We discuss why the expressions derived from the equations for the potentials have limited applicability. Our results are in agreement with previous analytical estimates and show that the QSA is a reliable tool and can be used for comparison with current and future observations to constrain models beyond $\Lambda$CDM., Comment: The manuscript is divided into two parts: the main paper (37 pages, 3 figures, 1 table) and the Supplementary Materials (24 pages, 1 table). Improved discussion on the validity of the quasi-static approximation. Matches the published version on JCAP

Spectral distortions (SDs) of the cosmic microwave background (CMB) provide a powerful tool for studying particle physics. Here we compute the distortion signals from decaying particles that convert directly into photons at different epochs during cosmic history, focusing on injection energies $E_\mathrm{inj}\lesssim 20\,\mathrm{keV}$. We deliver a comprehensive library of SD solutions that can be used to study a wide range of particle physics scenarios. We use {\tt CosmoTherm} to compute the SD signals, including effects on the ionization history and opacities of the Universe. We also consider the effect of blackbody-induced stimulated decay, which can modify the injection history significantly. Then, we use data from COBE/FIRAS and EDGES to constrain the properties of the decaying particles. We explore scenarios where these provide a dark matter (DM) candidate or constitute only a small fraction of DM. We complement the SD constraints with CMB anisotropy constraints, highlighting new effects from injections at very-low photon energies ($h\nu\lesssim 10^{-4}\,\mathrm{eV}$). Our model-independent constraints exhibit rich structures in the lifetime-energy domain, covering injection energies $E_\mathrm{inj}\simeq 10^{-10}\mathrm{eV}-10\mathrm{keV}$ and lifetimes $\tau_X\simeq 10^5\,\mathrm{s}-10^{33}\mathrm{s}$. We discuss the constraints on axions and axion-like particles that convert directly into two photons, revising existing SD constraints in the literature. Our limits are competitive with other constraints for axion masses $m_a c^2\gtrsim 27\,\mathrm{eV}$ and we find that simple estimates based on the overall energetics are generally inaccurate. Future CMB spectrometers could significantly improve the obtained constraints, thus providing an important complementary probe of early-universe particle physics., Comment: match published version

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., Comment: 31 pages, 14 figures and 2 tables. Accepted by Phys. Rev. D. Comments welcome

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 $07 \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., Comment: 14 pages, 13 figures, accepted for publication in MNRAS on March 3 2020

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., 20 pages, 14 figures, 4 tables. Accepted in MNRAS. Links to Data Release catalogues will be updated when available upon publication. Before this, catalogues are available on request from the authors

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., Comment: 17 pages, 2 figures, prepared for submission to JCAP

The Two Higgs Doublet Model (2HDM) with spontaneously broken $Z_2$ symmetry predicts a production of domain walls at the electroweak scale. We derive cosmological constraints on model parameters for both Type-I and Type-II 2HDMs from the requirement that domain walls do not dominate the Universe by the present day. For Type-I 2HDMs, we deduce the lower bound on the key parameter $\tan\beta > 10^5$ for a wide range of Higgs-boson masses $\sim$ 100 GeV or greater close to the Standard Model alignment limit. In addition, we perform numerical simulations of the 2HDM with an approximate as well as an exact $Z_2$ symmetry but biased initial conditions. In both cases, we find that domain wall networks are unstable and, hence, do not survive at late times. The domain walls experience an exponential suppression of scaling in these models which can help ameliorate the stringent constraints found in the case of an exact discrete symmetry. For a 2HDM with softly-broken $Z_2$ symmetry, we relate the size of this exponential suppression to the soft-breaking bilinear parameter $m_{12}$ allowing limits to be placed on this parameter of order $\mu$eV, such that domain wall domination can be avoided. In particular, for Type-II 2HDMs, we obtain a corresponding lower limit on the CP-odd phase $\theta$ generated by QCD instantons, $\theta \ \stackrel{>}{{}_\sim}\ 10^{-11}/(\sin\beta \cos\beta)$, which is in some tension with the upper limit of $\theta \ \stackrel{, Comment: 18 pages, 6 figures

Planck Collaboration: et al., arXiv:1807.06209v3, We present cosmological parameter results from the final full-mission Planck measurements of the cosmic microwave background (CMB) anisotropies, combining information from the temperature and polarization maps and the lensing reconstruction. Compared to the 2015 results, improved measurements of large-scale polarization allow the reionization optical depth to be measured with higher precision, leading to significant gains in the precision of other correlated parameters. Improved modelling of the small-scale polarization leads to more robust constraints on many parameters, with residual modelling uncertainties estimated to affect them only at the 0.5σ level. We find good consistency with the standard spatially-flat 6-parameter ΛCDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density Ωch2 = 0.120 ± 0.001, baryon density Ωbh2 = 0.0224 ± 0.0001, scalar spectral index ns = 0.965 ± 0.004, and optical depth τ = 0.054 ± 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θ* = 1.0411 ± 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-ΛCDM cosmology, the inferred (model-dependent) late-Universe parameters are: Hubble constant H0 = (67.4 ± 0.5) km s−1 Mpc−1; matter density parameter Ωm = 0.315 ± 0.007; and matter fluctuation amplitude σ8 = 0.811 ± 0.006. We find no compelling evidence for extensions to the base-ΛCDM model. Combining with baryon acoustic oscillation (BAO) measurements (and considering single-parameter extensions) we constrain the effective extra relativistic degrees of freedom to be Neff = 2.99 ± 0.17, in agreement with the Standard Model prediction Neff = 3.046, and find that the neutrino mass is tightly constrained to ∑mν < 0.12 eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base ΛCDM at over 2σ, which pulls some parameters that affect the lensing amplitude away from the ΛCDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. The joint constraint with BAO measurements on spatial curvature is consistent with a flat universe, ΩK = 0.001 ± 0.002. Also combining with Type Ia supernovae (SNe), the dark-energy equation of state parameter is measured to be w0 = −1.03 ± 0.03, consistent with a cosmological constant. We find no evidence for deviations from a purely power-law primordial spectrum, and combining with data from BAO, BICEP2, and Keck Array data, we place a limit on the tensor-to-scalar ratio r0.002 < 0.06. Standard big-bang nucleosynthesis predictions for the helium and deuterium abundances for the base-ΛCDM cosmology are in excellent agreement with observations. The Planck base-ΛCDM results are in good agreement with BAO, SNe, and some galaxy lensing observations, but in slight tension with the Dark Energy Survey’s combined-probe results including galaxy clustering (which prefers lower fluctuation amplitudes or matter density parameters), and in significant, 3.6σ, tension with local measurements of the Hubble constant (which prefer a higher value). Simple model extensions that can partially resolve these tensions are not favoured by the Planck data., The Planck Collaboration acknowledges the support of: ESA; CNES, and CNRS/INSU-IN2P3-INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, JA, and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC and PRACE (EU).

We present the first detailed case study using quadratic estimators (QE) to diagnose and remove systematics present in observed Cosmic Microwave Background (CMB) maps. In this work we focus on the temperature to polarization leakage. We use an iterative QE analysis to remove systematics, in analogy to de-lensing, recovering the primordial B-mode signal and the systematic maps. We introduce a new Gaussian filtering scheme crucial to stable convergence of the iterative cleaning procedure and validate with comparisons to semi-analytical forecasts. We study the limitations of this method by examining its performance both on idealized simulations and on more realistic, non-ideal simulations, where we assume varying de-lensing efficiencies. Finally, we quantify the systematic cleaning efficiency by presenting a likelihood analysis on the tensor to scalar ratio, $r$, and demonstrate that the blind cleaning results in an un-biased measurement of $r$, reducing the systematic induced B-mode power by nearly two orders of magnitude., 30 pages, 12 figures, prepared for submission to JCAP

Model-independent constraints on modified gravity models hitherto exist mainly on linear scales. A recently developed formalism presented a consistent parameterisation that is valid on all scales. Using this approach, we perform model-independent modified gravity $N$-body simulations on all cosmological scales with a time-dependent $\mu$. We present convergence tests of our simulations, and we examine how well existing fitting functions reproduce the non-linear matter power spectrum of the simulations. We find that although there is a significant variation in the accuracy of all of the fitting functions over the parameter space of our simulations, the ReACT framework delivers the most consistent performance for the matter power spectrum. We comment on how this might be improved to the level required for future surveys such as Euclid and the Vera Rubin Telescope (LSST). We also show how to compute weak-lensing observables consistently from the simulated matter power spectra in our approach, and show that ReACT also performs best when fitting the weak-lensing observables. This paves the way for a full model-independent test of modified gravity using all of the data from such upcoming surveys., Comment: 32 pages, 14 Figures; Comments welcome

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=-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., Comment: 19 pages, 11 figures, 3 tables. Accepted by MNRAS

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., Comment: 18 pages, 7 figures, 6 tables. Updated to match version accepted by MNRAS

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., Comment: 43 pages, 13 figures. It matches the accepted version on JCAP

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., Comment: 20 pages, 9 figures, 4 tables. Public version of the code will be released upon acceptance by PRD

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., Comment: 17 pages, 14 figures, 3 tables. Accepted for publication in MNRAS

The CORE collaboration thanks CNES, Thales Alenia Space, and Air Liquide Advanced Technologies for advice and technical support during the preparation of the CORE proposal. We also thank the ESA CDF team for the CMB Polarisation CDF study performed in March 2016, the results of which were extensively used to define the mission concept presented in this paper. J.G.N. acknowledges financial support from the Spanish MINECO for a Ramon y Cajal fellowship (RYC-2013-13256) and the I+D 2015 project AYA2015-65887-P (MINECO/FEDER). CJM is supported by an FCT Research Professorship, contract reference IF/00064/2012, funded by FCT/MCTES (Portugal) and POPH/FSE. F.J.C., R.F.-C., E.M.-G. and P.V. acknowledge support from the Spanish Ministerio de Econom´ıa y Competitividad project ESP2015-70646-C2-1-R (cofinanced with EU FEDER funds), ConsoliderIngenio 2010 project CSD2010-00064 and from the CSIC “Proyecto Intramural Especial” project 201550E091. FA is supported by the National Taiwan University (NTU) under Project No. 103R4000 and by the NTU Leung Center for Cosmology and Particle Astrophysics (LeCosPA) under Project No. FI121. BFR acknowledges support from the National Science Centre, Poland, under grant 2014/13/B/ST9/00845., Macìas-Pérez, J.F., Delabrouille, J., De Bernardis, P., Bouchet, F.R., Achúcarro, A., Ade, P.A.R., Allison, R., Arroja, F., Artal, E., Ashdown, M., Baccigalupi, C., Ballardini, M., Banday, A.J., Banerji, R., Barbosa, D., Bartlett, J., Bartolo, N., Basak, S., Baselmans, J.J.A., Basu, K., Battistelli, E.S., Battye, R., Baumann, D., Benoít, A., Bersanelli, M., Bideaud, A., Biesiada, M., Bilicki, M., Bonaldi, A., Bonato, M., Borrill, J., Boulanger, F., Brinckmann, T., Brown, M.L., Bucher, M., Burigana, C., Buzzelli, A., Cabass, G., Cai, Z.-Y., Calvo, M., Caputo, A., Carvalho, C.-S., Casas, F.J., Castellano, G., Catalano, A., Challinor, A., Charles, I., Chluba, J., Clements, D.L., Clesse, S., Colafrancesco, S., Colantoni, I., Contreras, D., Coppolecchia, A., Crook, M., D'Alessandro, G., D'Amico, G., Silva, A.D., De Avillez, M., De Gasperis, G., Petris, M.D., De Zotti, G., Danese, L., Désert, F.-X., Desjacques, V., Valentino, E.D., Dickinson, C., Diego, J.M., Doyle, S., Durrer, R., Dvorkin, C., Eriksen, H.K., Errard, J., Feeney, S., Fernández-Cobos, R., Finelli, F., Forastieri, F., Franceschet, C., Fuskeland, U., Galli, S., Génova-Santos, R.T., Gerbino, M., Giusarma, E., Gomez, A., González-Nuevo, J., Grandis, S., Greenslade, J., Goupy, J., Hagstotz, S., Hanany, S., Handley, W., Henrot-Versillé, S., Hernández-Monteagudo, C., Hervias-Caimapo, C., Hills, M., Hindmarsh, M., Hivon, E., Hoang, D.T., Hooper, D.C., Hu, B., Keihänen, E., Keskitalo, R., Kiiveri, K., Kisner, T., Kitching, T., Kunz, M., Kurki-Suonio, H., Lagache, G., Lamagna, L., Lapi, A., Lasenby, A., Lattanzi, M., Brun, A.M.C.L., Lesgourgues, J., Liguori, M., Lindholm, V., Lizarraga, J., Luzzi, G., Macìas-Pérez, J.F., Maffei, B., Mandolesi, N., Martin, S., Martinez-Gonzalez, E., Martins, C.J.A.P., Masi, S., Massardi, M., Matarrese, S., Mazzotta, P., McCarthy, D., Melchiorri, A., Melin, J.-B., Mennella, A., Mohr, J., Molinari, D., Monfardini, A., Montier, L., Natoli, P., Negrello, M., Notari, A., Noviello, F., Oppizzi, F., O'Sullivan, C., Pagano, L., Paiella, A., Pajer, E., Paoletti, D., Paradiso, S., Partridge, R.B., Patanchon, G., Patil, S.P., Perdereau, O., Piacentini, F., Piat, M., Pisano, G., Polastri, L., Polenta, G., Pollo, A., Ponthieu, N., Poulin, V., Prêle, D., Quartin, M., Ravenni, A., Remazeilles, M., Renzi, A., Ringeval, C., Roest, D., Roman, M., Roukema, B.F., Rubiño-Martin, J.-A., Salvati, L., Scott, D., Serjeant, S., Signorelli, G., Starobinsky, A.A., Sunyaev, R., Tan, C.Y., Tartari, A., Tasinato, G., Toffolatti, L., Tomasi, M., Torrado, J., Tramonte, D., Trappe, N., Triqueneaux, S., Tristram, M., Trombetti, T., Tucci, M., Tucker, C., Urrestilla, J., Väliviita, J., De Weygaert, R.V., Tent, B.V., Vennin, V., Verde, L., Vermeulen, G., Vielva, P., Vittorio, N., Voisin, F., Wallis, C., Wandelt, B., Wehus, I.K., Weller, J., Young, K., Zannoni, M.

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., Comment: 23 pages; 11 figures. Matches version accepted in PRD

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., Comment: 10 pages

We thank the anonymous referee for useful comments, which helped to clarify the paper. Some of the results in this paper have been derived using the HEALPix package. Parts of the cosmological analysis was made using Cosmo MC, CLASS and MontePython. C.H.-M. acknowledges financial support of the Spanish Ministry of Economy and Competitiveness via I+D project AYA-2015-66211-C2-2-P. CJM is supported by an FCT Research Professorship, contract reference IF/00064/2012, funded by FCT/MCTES (Portugal) and POPH/FSE. J.G.N. acknowledges financial support from the Spanish MINECO for a ‘Ramon y Cajal’ fellowship (RYC-2013-13256) and the I+D 2015 project AYA2015-65887-P (MINECO/FEDER)., Melin, J.-B., Bonaldi, A., Remazeilles, M., Hagstotz, S., Diego, J.M., Hernández-Monteagudo, C., Génova-Santos, R.T., Luzzi, G., Martins, C.J.A.P., Grandis, S., Mohr, J.J., Bartlett, J.G., Delabrouille, J., Ferraro, S., Tramonte, D., Rubiño-Martín, J.A., Macìas-Pérez, J.F., Achúcarro, A., Ade, P., Allison, R., Ashdown, M., Ballardini, M., Banday, A.J., Banerji, R., Bartolo, N., Basak, S., Basu, K., Battye, R.A., Baumann, D., Bersanelli, M., Bonato, M., Borrill, J., Bouchet, F., Boulanger, F., Brinckmann, T., Bucher, M., Burigana, C., Buzzelli, A., Cai, Z.-Y., Calvo, M., Carvalho, C.S., Castellano, M.G., Challinor, A., Chluba, J., Clesse, S., Colafrancesco, S., Colantoni, I., Coppolecchia, A., Crook, M., D'Alessandro, G., De Bernardis, P., De Gasperis, G., Petris, M.D., Zotti, G.D., Valentino, E.D., Errard, J., Feeney, S.M., Fernández-Cobos, R., Finelli, F., Forastieri, F., Galli, S., Gerbino, M., González-Nuevo, J., Greenslade, J., Hanany, S., Handley, W., Hervias-Caimapo, C., Hills, M., Hivon, E., Kiiveri, K., Kisner, T., Kitching, T., Kunz, M., Kurki-Suonio, H., Lamagna, L., Lasenby, A., Lattanzi, M., Brun, A.M.C.L., Lesgourgues, J., Lewis, A., Liguori, M., Lindholm, V., Lopez-Caniego, M., Maffei, B., Martinez-Gonzalez, E., Masi, S., Mazzotta, P., McCarthy, D., Melchiorri, A., Molinari, D., Monfardini, A., Natoli, P., Negrello, M., Notari, A., Paiella, A., Paoletti, D., Patanchon, G., Piat, M., Pisano, G., Polastri, L., Polenta, G., Pollo, A., Poulin, V., Quartin, M., Roman, M., Salvati, L., Tartari, A., Tomasi, M., Trappe, N., Triqueneaux, S., Trombetti, T., Tucker, C., Väliviita, J., De Weygaert, R.V., Tent, B.V., Vennin, V., Vielva, P., Vittorio, N., Weller, J., Young, K., Zannoni, M.

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., Comment: 11 pages, 6 figures

The discovery of apparent cosmological acceleration has spawned a huge number of dark energy and modified gravity theories. The f(R) models of gravity are obtained when one replaces the Ricci scalar in the Einstein-Hilbert action by an arbitrary function of the Ricci scalar. In this work we provide expressions for the equations of state (EoS) for perturbations which completely characterize the linearized perturbations in f(R) gravity, including the scalar, vector, and tensor modes. The EoS formalism is a powerful and elegant parametrization where the modification to General Relativity are treated as a dark-energy-fluid. The perturbed dark-energy-fluid variables such as the anisotropic stress or the entropy perturbation are explicitly given in terms of parameters of the model of interest.

Discordance in the $\mathrm{\ensuremath{\Lambda}}$ cold dark matter cosmological model can be seen by comparing parameters constrained by cosmic microwave background (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 Kullback-Leibler divergence, comparison of Bayesian evidences and other statistics which are sensitive to the mean, variance and shape of the distributions. However, as a byproduct, we present an update to the similar analysis in [R. A. Battye, T. Charnock, and A. Moss, Phys. Rev. D 91, 103508 (2015)], where we find that, considering new data and treatment of priors, the constraints from the CMB and from a combination of large scale structure (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 $\mathrm{Planck}2015+\mathrm{Pol}+\mathrm{BAO}$ parameter distributions can be quantified at a $\ensuremath{\sim}2.55\ensuremath{\sigma}$ tension using the method introduced in [R. A. Battye, T. Charnock, and A. Moss, Phys. Rev. D 91, 103508 (2015)]. If instead we use the distributions constrained by the combination of LSS probes which are in greatest agreement with those from $\mathrm{Planck}2015+\mathrm{Pol}+\mathrm{BAO}$ this tension is only $0.76\ensuremath{\sigma}$.

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., Comment: 22 pages, 15 figures, 2 tables

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., 21 pages, 11 figures, Comments welcome

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., Comment: 26 pages, 4 figures, 3 tables, submitted to PRD

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.

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., Comment: 17 pages, 6 figures. Accepted by Phys. Rev. D

We estimate the bispectrum of the Very Small Array data from the compact and extended configuration observations released in December 2002, and compare our results to those obtained from Gaussian simulations. There is a slight excess of large bispectrum values for two individual fields, but this does not appear when the fields are combined. Given our expected level of residual point sources, we do not expect these to be the source of the discrepancy. Using the compact configuration data, we put an upper limit of 5400 on the value of f_NL, the non-linear coupling parameter, at 95 per cent confidence. We test our bispectrum estimator using non-Gaussian simulations with a known bispectrum, and recover the input values., Comment: 17 pages, 16 figures, replaced with version accepted by MNRAS. Primordial bispectrum recalculated and figure 11 changed

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., Comment: 25 pages, 18 figures, 10 tables. Accepted for publication in MNRAS

We study the equations for the evolution of cosmological perturbations in $f(\mathcal{R})$ and conclude that this modified gravity model can be expressed as a dark energy fluid at background and linearized perturbation order. By eliminating the extra scalar degree of freedom known to be present in such theories, we are able to characterize 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.

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 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.

We have studied the implications of high-sensitivity polarization measurements of objects from the Wilkinson Microwave Anisotropy Probe (WMAP) point source catalogue made using the Very Large Array (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 three VLA frequencies, we have investigated the spectral behaviour as a function of frequency by means of a two-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 two-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 NRAO VLA Sky Survey (NVSS). In addition we compute 45 rotation measures from polarization position angles which are compatible with a ?� 2 law. We use our results to predict the level of point source confusion noise that contaminates cosmic microwave background 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.

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.

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., Comment: 16 pages, 5 figures. Accepted by MNRAS. Two large figures are available at http://www.jb.man.ac.uk/~njj/pol_maps and http://www.jb.man.ac.uk/~njj/pol_spec

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 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 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.