We present the radio luminosity function (RLF) of optically-selected quasars below 1~mJy, constructed by applying a Bayesian-fitting stacking technique to objects well below the nominal radio flux-density limit. We test the technique using simulated data, confirming that we can reconstruct the RLF over three orders of magnitude below the typical $5\sigma$ detection threshold. We apply our method to 1.4-GHz flux-densities from the Faint Images of the Radio Sky at Twenty-cm survey (FIRST), extracted at the positions of optical quasars from the Sloan Digital Sky Survey (SDSS) over seven redshift bins up to $z=2.15$ {and measure the RLF down to two orders of magnitude below the FIRST detection threshold}. In the lowest redshift bin ($0.2<z<0.45$), we find that our measured RLF agrees well with deeper data from the literature. The RLF for the radio-loud quasars flattens below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 25.5$ and becomes steeper again below $\log_{10}[L_{1.4}/{\rm W\,Hz}^{-1}] \approx 24.8$, where radio-quiet quasars start to emerge. The radio luminosity where radio-quiet quasars emerge coincides with the luminosity where star-forming galaxies are expected to start to dominate the radio source counts. This implies that there could be a significant contribution from star formation in the host galaxies, but additional data is required to investigate this further. The higher-redshift bins show a similar behaviour as for the lowest-$z$ bin, implying that the same physical process may be responsible.

In this paper we present the enhanced X-ray Timing and Polarimetry mission - eXTP. eXTP is a space science mission designed to study fundamental physics under extreme conditions of density, gravity and magnetism. The mission aims at determining the equation of state of matter at supra-nuclear density, measuring effects of QED, and understanding the dynamics of matter in strong-field gravity. In addition to investigating fundamental physics, eXTP will be a very powerful observatory for astrophysics that will provide observations of unprecedented quality on a variety of galactic and extragalactic objects. In particular, its wide field monitoring capabilities will be highly instrumental to detect the electro-magnetic counterparts of gravitational wave sources. The paper provides a detailed description of: (1) the technological and technical aspects, and the expected performance of the instruments of the scientific payload; (2) the elements and functions of the mission, from the spacecraft to the ground segment.

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

Multi-wavelength astronomical studies brings a wealth of science within reach. One way to achieve a cross-wavelength analysis is via `stacking', i.e. combining precise positional information from an image at one wavelength with data from one at another wavelength in order to extract source-flux distributions and other derived quantities. For the first time we extend stacking to include the effects of confusion. We develop our algorithm in a fully Bayesian framework and apply it to the Square Kilometre Array Design Study (SKADS) simulation in order to extract galaxy number counts. Previous studies have shown that recovered source counts are highly biased high when confusion is non-negligible. With this new method, source counts are returned correctly. We also describe a novel estimator for quantifying the impact of confusion on stacking analyses. This method is an essential step in exploiting scientific return for upcoming deep radio surveys, e.g. MIGHTEE on MeerKAT.

The MeerKAT telescope will be one of the most sensitive radio arrays in the pre-SKA era. Here we discuss a low-frequency SZ-selected cluster survey with MeerKAT, the MeerKAT Extended Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey. The primary goal of this survey is to detect faint signatures of diffuse cluster emission, specifically radio halos and relics. SZ-selected cluster samples offer a homogeneous, mass-limited set of targets out to higher redshift than X-ray samples. MeerKAT is sensitive enough to detect diffuse radio emission at the faint levels expected in low-mass and high-redshift clusters, thereby enabling radio halo and relic formation theories to be tested with a larger statistical sample over a significantly expanded phase space. Complementary multiwavelength follow-up observations will provide a more complete picture of any clusters found to host diffuse emission, thereby enhancing the scientific return of the MERGHERS survey.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

Almost two decades of observations of radio emission in galaxy clusters have proven the existence of relativistic particles and magnetic fields that generate extended synchrotron emission in the form of radio halos. In the current scenario, radio halos are generated through re--acceleration of relativistic electrons by turbulence generated by cluster mergers. Although this theoretical framework has received increasingly supporting observational evidence over the last ten years, observations of statistically complete samples are needed in order to fundamentally test model predictions. In this paper we briefly review our 7--element Karoo Radio Telescope observations of a sample of nearby clusters aimed to test the predictions of the turbulent re--acceleration model in small systems ($M_{500} > 4 \times 10^{14}$ M$_{\odot}$). We conclude by presenting two galaxy cluster surveys to be carried out with MeerKAT in order to provide crucial test of models of radio halo formation in nearby ($z < 0.1$) and high redshift ($z > 0.4$) systems respectively.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

We analyse Very Long Baseline Interferometry (VLBI) observations of the blazar CGRaBS J0809+5341 using Bayesian inference methods. The observation was carried out at 5 GHz using 8 telescopes that form part of the European VLBI Network. Imaging and deconvolution using traditional methods imply that the blazar is unresolved. To search for source structure beyond the diffraction limit, we perform Bayesian model selection between three source models (point, elliptical Gaussian, and circular Gaussian). Our modelling jointly accounts for antenna-dependent gains and system equivalent flux densities. We obtain posterior distributions for the various source and instrumental parameters along with the corresponding uncertainties and correlations between them. We find that there is very strong evidence (>1e9 :1) in favour of elliptical Gaussian structure and using this model derive the apparent brightness temperature distribution of the blazar, accounting for uncertainties in the shape estimates. To test the integrity of our method, we also perform model selection on synthetic observations and use this to develop a Bayesian criterion for the minimum resolvable source size and consequently the maximum measurable brightness temperature for a given interferometer, dependent on the signal-to-noise ratio (SNR) of the data incorporating the aforementioned systematics. We find that calibration errors play an increasingly important role in determining the over-resolution limit for SNR>>100. We show that it is possible to exploit the resolving power of future VLBI arrays down to about 5 per cent of the size of the naturally-weighted restoring beam, if the gain calibration is precise to 1 per cent or less.

The birth of the first luminous sources and the ensuing epoch of reionization are best studied via the redshifted 21-cm emission line, the signature of the first two imprinting the last. In this work we present a fully-Bayesian method, \textsc{hibayes}, for extracting the faint, global (sky-averaged) 21-cm signal from the much brighter foreground emission. We show that a simplified (but plausible), Gaussian model of the 21-cm emission from the Cosmic Dawn epoch ($15 \lesssim z \lesssim 30$), parameterized by an amplitude $A_{\rm HI}$, a frequency peak $\nu_{\rm HI}$ and a width $\sigma_{\rm HI}$, can be extracted even in the presence of a structured foreground frequency spectrum (parameterized as a $7^{\rm th}$-order polynomial), provided sufficient signal-to-noise (400~hours of observation with a single dipole). We apply our method to an early, 19-minute long observation from the Large aperture Experiment to detect the Dark Ages, constraining the 21-cm signal amplitude and width to be $-890 < A_{\rm HI} < 0$ mK and $\sigma_{\rm HI} > 6.5$ MHz (corresponding to $\Delta z > 1.9$ at redshift $z \simeq 20$) respectively at the 95-per-cent confidence level in the range $13.2 < z < 27.4$ ($100 > \nu > 50$ MHz).

The construction of the KAT-7 array in the Karoo region of the Northern Cape in South Africa was intended primarily as an engineering prototype for technologies and techniques applicable to the MeerKAT telescope. This paper looks at the main engineering and scien- tific highlights from this effort, and discusses their applicability to both MeerKAT and other next-generation radio telescopes. In particular we found that the composite dish surface works well, but it becomes complicated to fabricate for a dish lacking circular symmetry; the Stir- ling cycle cryogenic system with ion pump to achieve vacuum works but demands much higher maintenance than an equivalent Gifford-McMahon cycle system; the ROACH (Recon- figurable Open Architecture Computing Hardware)-based correlator with SPEAD (Stream- ing Protocol for Exchanging Astronomical Data) protocol data transfer works very well and KATCP (Karoo Array Telescope Control Protocol) control protocol has proven very flexible and convenient. KAT-7 has also been used for scientific observations where it has a niche in mapping low surface-brightness continuum sources, some extended HI halos and OH masers in star-forming regions. It can also be used to monitor continuum source variability, observe pulsars, and make VLBI observations

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Bo\"otes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5\sigma$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

Radio interferometers suffer from the problem of missing information in their data, due to the gaps between the antennas. This results in artifacts, such as bright rings around sources, in the images obtained. Multiple deconvolution algorithms have been proposed to solve this problem and produce cleaner radio images. However, these algorithms are unable to correctly estimate uncertainties in derived scientific parameters or to always include the effects of instrumental errors. We propose an alternative technique called Bayesian Inference for Radio Observations (BIRO) which uses a Bayesian statistical framework to determine the scientific parameters and instrumental errors simultaneously directly from the raw data, without making an image. We use a simple simulation of Westerbork Synthesis Radio Telescope data including pointing errors and beam parameters as instrumental effects, to demonstrate the use of BIRO.

We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information.

We present Montblanc, a GPU implementation of the Radio interferometer measurement equation (RIME) in support of the Bayesian inference for radio observations (BIRO) technique. BIRO uses Bayesian inference to select sky models that best match the visibilities observed by a radio interferometer. To accomplish this, BIRO evaluates the RIME multiple times, varying sky model parameters to produce multiple model visibilities. Chi-squared values computed from the model and observed visibilities are used as likelihood values to drive the Bayesian sampling process and select the best sky model. As most of the elements of the RIME and chi-squared calculation are independent of one another, they are highly amenable to parallel computation. Additionally, Montblanc caters for iterative RIME evaluation to produce multiple chi-squared values. Modified model parameters are transferred to the GPU between each iteration. We implemented Montblanc as a Python package based upon NVIDIA's CUDA architecture. As such, it is easy to extend and implement different pipelines. At present, Montblanc supports point and Gaussian morphologies, but is designed for easy addition of new source profiles. Montblanc's RIME implementation is performant: On an NVIDIA K40, it is approximately 250 times faster than MeqTrees on a dual hexacore Intel E5-2620v2 CPU. Compared to the OSKAR simulator's GPU-implemented RIME components it is 7.7 and 12 times faster on the same K40 for single and double-precision floating point respectively. However, OSKAR's RIME implementation is more general than Montblanc's BIRO-tailored RIME. Theoretical analysis of Montblanc's dominant CUDA kernel suggests that it is memory bound. In practice, profiling shows that is balanced between compute and memory, as much of the data required by the problem is retained in L1 and L2 cache.

New telescopes like the Square Kilometre Array (SKA) will push into a new sensitivity regime and expose systematics, such as direction-dependent effects, that could previously be ignored. Current methods for handling such systematics rely on alternating best estimates of instrumental calibration and models of the underlying sky, which can lead to inadequate uncertainty estimates and biased results because any correlations between parameters are ignored. These deconvolution algorithms produce a single image that is assumed to be a true representation of the sky, when in fact it is just one realization of an infinite ensemble of images compatible with the noise in the data. In contrast, here we report a Bayesian formalism that simultaneously infers both systematics and science. Our technique, Bayesian Inference for Radio Observations (BIRO), determines all parameters directly from the raw data, bypassing image-making entirely, by sampling from the joint posterior probability distribution. This enables it to derive both correlations and accurate uncertainties, making use of the flexible software MEQTREES to model the sky and telescope simultaneously. We demonstrate BIRO with two simulated sets of Westerbork Synthesis Radio Telescope data sets. In the first, we perform joint estimates of 103 scientific (flux densities of sources) and instrumental (pointing errors, beamwidth and noise) parameters. In the second example, we perform source separation with BIRO. Using the Bayesian evidence, we can accurately select between a single point source, two point sources and an extended Gaussian source, allowing for 'super-resolution' on scales much smaller than the synthesized beam.

Measuring radio source counts is critical for characterizing new extragalactic populations, brings a wealth of science within reach and will inform forecasts for SKA and its pathfinders. Yet there is currently great debate (and few measurements) about the behaviour of the 1.4-GHz counts in the microJy regime. One way to push the counts to these levels is via 'stacking', the covariance of a map with a catalogue at higher resolution and (often) a different wavelength. For the first time, we cast stacking in a fully bayesian framework, applying it to (i) the SKADS simulation and (ii) VLA data stacked at the positions of sources from the VIDEO survey. In the former case, the algorithm recovers the counts correctly when applied to the catalogue, but is biased high when confusion comes into play. This needs to be accounted for in the analysis of data from any relatively-low-resolution SKA pathfinders. For the latter case, the observed radio source counts remain flat below the 5-sigma level of 85 microJy as far as 40 microJy, then fall off earlier than the flux hinted at by the SKADS simulations and a recent P(D) analysis (which is the only other measurement from the literature at these flux-density levels, itself extrapolated in frequency). Division into galaxy type via spectral-energy distribution reveals that normal spiral galaxies dominate the counts at these fluxes.

Astronomy at or below the 'survey threshold' has expanded significantly since the publication of the original 'Science with the Square Kilometer Array' in 1999 and its update in 2004. The techniques in this regime may be broadly (but far from exclusively) defined as 'confusion' or 'P(D)' analyses (analyses of one-point statistics), and 'stacking', accounting for the flux-density distribution of noise-limited images co-added at the positions of objects detected/isolated in a different waveband. Here we discuss the relevant issues, present some examples of recent analyses, and consider some of the consequences for the design and use of surveys with the SKA and its pathfinders.

We present polarization measurements of extragalactic radio sources observed during the Cosmic Microwave Background polarization survey of the Q/U Imaging Experiment (QUIET), operating at 43 GHz (Q-band) and 95 GHz (W-band). We examine sources selected at 20 GHz from the public, $>$40 mJy catalog of the Australia Telescope (AT20G) survey. There are $\sim$480 such sources within QUIET's four low-foreground survey patches, including the nearby radio galaxies Centaurus A and Pictor A. The median error on our polarized flux density measurements is 30--40 mJy per Stokes parameter. At S/N $> 3$ significance, we detect linear polarization for seven sources in Q-band and six in W-band; only $1.3 \pm 1.1$ detections per frequency band are expected by chance. For sources without a detection of polarized emission, we find that half of the sources have polarization amplitudes below 90 mJy (Q-band) and 106 mJy (W-band), at 95% confidence. Finally, we compare our polarization measurements to intensity and polarization measurements of the same sources from the literature. For the four sources with WMAP and Planck intensity measurements $>1$ Jy, the polarization fraction are above 1% in both QUIET bands. At high significance, we compute polarization fractions as much as 10--20% for some sources, but the effects of source variability may cut that level in half for contemporaneous comparisons. Our results indicate that simple models---ones that scale a fixed polarization fraction with frequency---are inadequate to model the behavior of these sources and their contributions to polarization maps.

We measure star-formation rates (SFRs) and specific SFRs (SSFRs) of Ks-selected galaxies from the VIDEO survey by stacking 1.4-GHz Very Large Array data. We split the sample, which spans 0 < z < 3 and stellar masses 10**8.0 < Mstellar/Msol < 10**11.5, into elliptical, irregular or starburst galaxies based on their spectral-energy distributions. We find that SSFR falls with stellar mass, in agreement with the `downsizing' paradigm. We consider the dependence of the SSFR-mass slope on redshift: for our full and elliptical samples the slope flattens, but for the irregular and starburst samples the slope is independent of redshift. The rate of SSFR evolution reduces slightly with stellar mass for ellipticals, but irregulars and starbursts co-evolve across stellar masses. Our results for SSFR as a function of stellar mass and redshift are in agreement with those derived from other radio-stacking measurements of mass-selected passive and star-forming galaxies, but inconsistent with those generated from semi-analytic models, which tend to underestimate SFRs and SSFRs. There is a need for deeper high-resolution radio surveys such as those from telescopes like the next-generation MeerKAT in order to probe lower masses at earlier times and to permit direct detections, i.e. to study individual galaxies in detail.

The Q/U Imaging ExperimenT (QUIET) has observed the cosmic microwave background (CMB) at 43 and 95GHz. The 43-GHz results have been published in QUIET Collaboration et al. (2011), and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers with a total array sensitivity of 87 uK sqrt(s). Four low-foreground fields were observed, covering a total of ~1000 square degrees with an effective angular resolution of 12.8', allowing for constraints on primordial gravitational waves and high-signal-to-noise measurements of the E-modes across three acoustic peaks. The data reduction was performed using two independent analysis pipelines, one based on a pseudo-Cl (PCL) cross-correlation approach, and the other on a maximum-likelihood (ML) approach. All data selection criteria and filters were modified until a predefined set of null tests had been satisfied before inspecting any non-null power spectrum. The results derived by the two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between l=25 and 975 and find that the EE spectrum is consistent with LCDM, while the BB power spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar ratio to r=1.1+0.9-0.8 (r<2.8 at 95% C.L.) as derived by the ML pipeline, and r=1.2+0.9-0.8 (r<2.7 at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the dust component of the Planck Sky Model, though the corresponding excess power is small compared to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that these correspond to a tensor-to-scalar ratio smaller than r=0.01, the lowest level yet reported in the literature.

In this paper we describe the first data release of the the Visible and Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations (VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and K_s bands, specifically designed to enable the evolution of galaxies and large structures to be traced as a function of both epoch and environment from the present day out to z=4, and active galactic nuclei (AGN) and the most massive galaxies up to and into the epoch of reionization. With its depth and area, VIDEO will be able to fully explore the period in the Universe where AGN and starburst activity were at their peak and the first galaxy clusters were beginning to virialize. VIDEO therefore offers a unique data set with which to investigate the interplay between AGN, starbursts and environment, and the role of feedback at a time when it was potentially most crucial. We provide data over the VIDEO-XMM3 tile, which also covers the Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5, J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of Y=24.6 will be reached within the full integration time in future releases). The data are compared to previous surveys over this field and we find good astrometric agreement with the Two-Micron All Sky Survey, and source counts in agreement with the recently released UltraVISTA survey data. The addition of the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of photometric redshifts and significantly reduces the fraction of catastrophic outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence of an i-band luminosity prior. (Truncated Abstract)

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.

The Q/U Imaging ExperimenT (QUIET) is designed to measure polarization in the Cosmic Microwave Background, targeting the imprint of inflationary gravitational waves at large angular scales (~ 1 degree). Between 2008 October and 2010 December, two independent receiver arrays were deployed sequentially on a 1.4 m side-fed Dragonian telescope. The polarimeters which form the focal planes use a highly compact design based on High Electron Mobility Transistors (HEMTs) that provides simultaneous measurements of the Stokes parameters Q, U, and I in a single module. The 17-element Q-band polarimeter array, with a central frequency of 43.1 GHz, has the best sensitivity (69 uK sqrt(s)) and the lowest instrumental systematic errors ever achieved in this band, contributing to the tensor-to-scalar ratio at r < 0.1. The 84-element W-band polarimeter array has a sensitivity of 87 uK sqrt(s) at a central frequency of 94.5 GHz. It has the lowest systematic errors to date, contributing at r < 0.01. The two arrays together cover multipoles in the range l= 25-975. These are the largest HEMT-based arrays deployed to date. This article describes the design, calibration, performance of, and sources of systematic error for the instrument.

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

We present an interesting Sunyaev-Zel'dovich (SZ) detection in the first of the Arcminute Microkelvin Imager (AMI) 'blind', degree-square fields to have been observed down to our target sensitivity of 100{\mu}Jy/beam. In follow-up deep pointed observations the SZ effect is detected with a maximum peak decrement greater than 8 \times the thermal noise. No corresponding emission is visible in the ROSAT all-sky X-ray survey and no cluster is evident in the Palomar all-sky optical survey. Compared with existing SZ images of distant clusters, the extent is large (\approx 10') and complex; our analysis favours a model containing two clusters rather than a single cluster. Our Bayesian analysis is currently limited to modelling each cluster with an ellipsoidal or spherical beta-model, which do not do justice to this decrement. Fitting an ellipsoid to the deeper candidate we find the following. (a) Assuming that the Evrard et al. (2002) approximation to Press & Schechter (1974) correctly gives the number density of clusters as a function of mass and redshift, then, in the search area, the formal Bayesian probability ratio of the AMI detection of this cluster is 7.9 \times 10^4:1; alternatively assuming Jenkins et al. (2001) as the true prior, the formal Bayesian probability ratio of detection is 2.1 \times 10^5:1. (b) The cluster mass is MT,200 = 5.5+1.2\times 10^14h-1M\odot. (c) Abandoning a physical model with num- -1.3 70 ber density prior and instead simply modelling the SZ decrement using a phenomenological {\beta}-model of temperature decrement as a function of angular distance, we find a central SZ temperature decrement of -295+36 {\mu}K - this allows for CMB primary anisotropies, receiver -15 noise and radio sources. We are unsure if the cluster system we observe is a merging system or two separate clusters.

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the Canada-France-Hawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.16--0.41. The cluster gas is modelled using the Sunyaev--Zel'dovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we: i) find very good agreement between SZ measurements (assuming large-scale virialisation and a gas-fraction prior) and lensing measurements of the total cluster masses out to r_200; ii) perform the first multiple-component weak-lensing analysis of A115; iii) confirm the unusual separation between the gas and mass components in A1914; iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived mass-temperature relation for parameterizing measurements of the SZ effect.

Most Sunyaev--Zel'dovich (SZ) and X-ray analyses of galaxy clusters try to constrain the cluster total mass and/or gas mass using parameterised models and assumptions of spherical symmetry and hydrostatic equilibrium. By numerically exploring the probability distributions of the cluster parameters given the simulated interferometric SZ data in the context of Bayesian methods, and assuming a beta-model for the electron number density we investigate the capability of this model and analysis to return the simulated cluster input quantities via three rameterisations. In parameterisation I we assume that the T is an input parameter. We find that parameterisation I can hardly constrain the cluster parameters. We then investigate parameterisations II and III in which fg(r200) replaces temperature as a main variable. In parameterisation II we relate M_T(r200) and T assuming hydrostatic equilibrium. We find that parameterisation II can constrain the cluster physical parameters but the temperature estimate is biased low. In parameterisation III, the virial theorem replaces the hydrostatic equilibrium assumption. We find that parameterisation III results in unbiased estimates of the cluster properties. We generate a second simulated cluster using a generalised NFW (GNFW) pressure profile and analyse it with an entropy based model to take into account the temperature gradient in our analysis and improve the cluster gas density distribution. This model also constrains the cluster physical parameters and the results show a radial decline in the gas temperature as expected. The mean cluster total mass estimates are also within 1 sigma from the simulated cluster true values. However, we find that for at least interferometric SZ analysis in practice at the present time, there is no differences in the AMI visibilities between the two models. This may of course change as the instruments improve.

Statistically anomalous signals in the microwave background have been extensively studied in general in multipole space, and in real space mainly for circular and other simple patterns. In this paper we search for a range of non-trivial patterns in the temperature data from WMAP 7-year observations. We find a very significant detection of a number of such features and discuss their consequences for the essential character of the cosmos.

We have observed an area of approximatley 27 deg^2 to an rms noise level of less than 0.2 mJy at 15.7 GHz, using the Arcminute Microkelvin Imager Large Array. These observations constitute the most sensitive radio-source survey of any extent (greater than approximately 0.2 deg^2) above 1.4 GHz. This paper presents the techniques employed for observing, mapping and source extraction. We have used a systematic procedure for extracting information and producing source catalogues, from maps with varying noise and uv-coverage. We have performed simulations to test our mapping and source-extraction procedures, and developed methods for identifying extended, overlapping and spurious sources in noisy images. In an accompanying paper, AMI Consortium: Davies et al. 2010, the first results from the 10C survey, including the deep 15.7-GHz source count, are presented.

The first results from the Tenth Cambridge (10C) Survey of Radio Sources, carried out using the AMI Large Array (LA) at an observing frequency of 15.7 GHz, are presented. The survey fields cover an area of approximately 27 sq. degrees to a flux-density completeness of 1 mJy. Results for some deeper areas, covering approximately 12 sq. degrees, wholly contained within the total areas and complete to 0.5 mJy, are also presented. The completeness for both areas is estimated to be at least 93 per cent. The source catalogue contains 1897 entries and is available at www.mrao.cam.ac.uk/surveys/10C. It has been combined with that of the 9C Survey to calculate the 15.7-GHz source counts. A broken power law is found to provide a good parameterisation of the differential count between 0.5 mJy and 1 Jy. The measured count has been compared to that predicted by de Zotti et al. (2005). The model displays good agreement with the data at the highest flux densities but under-predicts the integrated count between 0.5 mJy and 1 Jy by about 30 per cent. Entries from the source catalogue have been matched to those contained in the catalogues of NVSS and FIRST (both of which have observing frequencies of 1.4 GHz). This matching provides evidence for a shift in the typical 1.4-to-15.7-GHz spectral index of the 15.7-GHz-selected source population with decreasing flux density towards sub-mJy levels - the spectra tend to become less steep. Automated methods for detecting extended sources have been applied to the data; approximately 5 per cent of the sources are found to be extended relative to the LA synthesised beam of approximately 30 arcsec. Investigations using higher-resolution data showed that most of the genuinely extended sources at 16 GHz are classical doubles, although some nearby galaxies and twin-jet sources were also identified.

We present 13.9-18.2 GHz observations of the Sunyaev-Zel'dovich (SZ) effect towards Abell 2146 using the Arcminute Microkelvin Imager (AMI). The cluster is detected with a peak SNR ratio of 13 sigma in the radio source subtracted map. Comparison of the SZ and X-ray images suggests that they both have extended regions which lie approximately perpendicular to one another, with their emission peaks significantly displaced. These features indicate non-uniformities in the distributions of the gas temperature and pressure, indicative of a cluster merger. We use a Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model to obtain cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy; the probability of SZ + CMB primordial structure + radio sources + receiver noise to CMB + radio sources + receiver noise is 3 x 10^{6}:1. We compare the results from three different cluster models. Our preferred model exploits the observation that the gas fractions do not appear to vary greatly between clusters. Given the relative masses of the two merging systems in Abell 2146, the mean gas temperature can be deduced from the virial theorem (assuming all of the kinetic energy is in the form of internal gas energy) without being affected significantly by the merger event, provided the primary cluster was virialized before the merger. In this model we fit a simple spherical isothermal beta-model, despite the inadequacy of this model for a merging system like Abell 2146, and assume the cluster follows the mass-temperature relation of a virialized, singular, isothermal sphere. We note that this model avoids inferring large-scale cluster parameters internal to r_200 under the widely used assumption of hydrostatic equilibrium. We find that at r_200 M_T= 4.1 \pm 0.5 x 10^{14} h^{-1}M_sun and T=4.5 \pm 0.5 keV.

We perform deep 1.8 cm radio continuum imaging towards thirteen protostellar regions selected from the Spitzer c2d small clouds and cores programme at high resolution (25") in order to detect and quantify the cm-wave emission from deeply embedded young protostars. Within these regions we detect fifteen compact radio sources which we identify as radio protostars including two probable new detections. The sample is in general of low bolometric luminosity and contains several of the newly detected VeLLO sources. We determine the 1.8 cm radio luminosity to bolometric luminosity correlation, L_rad -L_bol, for the sample and discuss the nature of the radio emission in terms of the available sources of ionized gas. We also investigate the L_rad-L_IR correlation and suggest that radio flux density may be used as a proxy for the internal luminosity of low luminosity protostars.

We present observations using the Small Array of the Arcminute Microkelvin Imager (AMI; 14-18 GHz) of four Abell and three MACS clusters spanning 0.171-0.686 in redshift. We detect Sunyaev-Zel'dovich (SZ) signals in five of these without any attempt at source subtraction, although strong source contamination is present. With radio-source measurements from high-resolution observations, and under the assumptions of spherical $\beta$-model, isothermality and hydrostatic equilibrium, a Bayesian analysis of the data in the visibility plane detects extended SZ decrements in all seven clusters over and above receiver noise, radio sources and primary CMB imprints. Bayesian evidence ratios range from 10^{11}:1 to 10^{43}:1 for six of the clusters and 3000:1 for one with substantially less data than the others. We present posterior probability distributions for, e.g., total mass and gas fraction averaged over radii internal to which the mean overdensity is 1000, 500 and 200, r_200 being the virial radius. Reaching r_200 involves some extrapolation for the nearer clusters but not for the more-distant ones. We find that our estimates of gas fraction are low (compared with most in the literature) and decrease with increasing radius. These results appear to be consistent with the notion that gas temperature in fact falls with distance (away from near the cluster centre) out to the virial radius.

We report new cm-wave measurements at five frequencies between 15 and 18GHz of the continuum emission from the reportedly anomalous "region 4" of the nearby galaxy NGC6946. We find that the emission in this frequency range is significantly in excess of that measured at 8.5GHz, but has a spectrum from 15-18GHz consistent with optically thin free-free emission from a compact HII region. In combination with previously published data we fit four emission models containing different continuum components using the Bayesian spectrum analysis package radiospec. These fits show that, in combination with data at other frequencies, a model with a spinning dust component is slightly preferred to those that possess better-established emission mechanisms.

We analyze correlations between the first letter of the name of an author and the number of citations their papers receive. We look at simple mean counts, numbers of highly-cited papers, and normalized h-indices, by letter. To our surprise, we conclude that orthographically senior authors produce a better body of work than their colleagues, despite some evidence of discrimination against them.

We present 25 arcsecond resolution radio images of five Lynds Dark Nebulae (L675, L944, L1103, L1111 & L1246) at 16 GHz made with the Arcminute Microkelvin Imager (AMI) Large Array. These objects were previously observed with the AMI Small Array to have an excess of emission at microwave frequencies relative to lower frequency radio data. In L675 we find a flat spectrum compact radio counterpart to the 850 micron emission seen with SCUBA and suggest that it is cm-wave emission from a previously unknown deeply embedded young protostar. In the case of L1246 the cm-wave emission is spatially correlated with 8 micron emission seen with Spitzer. Since the MIR emission is present only in Spitzer band 4 we suggest that it arises from a population of PAH molecules, which also give rise to the cm-wave emission through spinning dust emission.

Observations at 14.2 to 17.9 GHz made with the AMI Small Array towards fourteen Lynds Dark Nebulae with a resolution of 2' are reported. These sources are selected from the SCUBA observations of Visser et al. (2001) as small angular diameter clouds well matched to the synthesized beam of the AMI Small Array. Comparison of the AMI observations with radio observations at lower frequencies with matched uv-plane coverage is made, in order to search for any anomalous excess emission which can be attributed to spinning dust. Possible emission from spinning dust is identified as a source within a 2' radius of the Scuba position of the Lynds dark nebula, exhibiting an excess with respect to lower frequency radio emission. We find five sources which show a possible spinning dust component in their spectra. These sources have rising spectral indices in the frequency range 14.2--17.9 GHz. Of these five one has already been reported, L1111, we report one new definite detection, L675, and three new probable detections (L944, L1103 and L1246). The relative certainty of these detections is assessed on the basis of three criteria: the extent of the emission, the coincidence of the emission with the Scuba position and the likelihood of alternative explanations for the excess. Extended microwave emission makes the likelihood of the anomalous emission arising as a consequence of a radio counterpart to a protostar or a proto-planetary disk unlikely. We use a 2' radius in order to be consistent with the IRAS identifications of dark nebulae (Parker 1988), and our third criterion is used in the case of L1103 where a high flux density at 850 microns relative to the FIR data suggests a more complicated emission spectrum.

We present follow-up observations of 97 point sources from the Wilkinson Microwave Anisotropy Probe (WMAP) 3-year data, contained within the New Extragalactic WMAP Point Source (NEWPS) catalogue between declinations of -4 and +60 degrees; the sources form a flux-density-limited sample complete to 1.1 Jy (approximately 5 sigma) at 33 GHz. Our observations were made at 16 GHz using the Arcminute Microkelvin Imager (AMI) and at 33 GHz with the Very Small Array (VSA). 94 of the sources have reliable, simultaneous -- typically a few minutes apart -- observations with both telescopes. The spectra between 13.9 and 33.75 GHz are very different from those of bright sources at low frequency: 44 per cent have rising spectra (alpha < 0.0), where flux density is proportional to frequency^-alpha, and 93 per cent have spectra with alpha < 0.5; the median spectral index is 0.04. For the brighter sources, the agreement between VSA and WMAP 33-GHz flux densities averaged over sources is very good. However, for the fainter sources, the VSA tends to measure lower values for the flux densities than WMAP. We suggest that the main cause of this effect is Eddington bias arising from variability.

Using the Arcminute Microkelvin Imager (AMI) at 16 GHz and the Very Small Array (VSA) at 33 GHz to make follow-up observations of sources in the New Extragalactic WMAP Point Source Catalogue, we have investigated the flux-density variability in a complete sample of 97 sources over timescales of a few months to approximately 1.5 years. We find that 53 per cent of the 93 sources, for which we have multiple observations, are variable, at the 99 per cent confidence level, above the flux density calibration uncertainties of approximately 4 per cent at 16 GHz; the fraction of sources having varied by more than 20 per cent is 15 per cent at 16 GHz and 20 per cent at 33 GHz. Not only is this common occurrence of variability at high frequency of interest for source physics, but strategies for coping with source contamination in CMB work must take this variability into account. There is no strong evidence of a correlation between variability and flux density for the sample as a whole. Using a maximum-likelihood method, we calculate the correlation in the variability at the two frequencies in a subset of sources classified as variable from both the AMI and VSA data and find the Pearson product-moment correlation coefficient to be very high (0.955 +/- 0.034). We also find the degree of variability at 16 GHz (0.202 +/- 0.028) to be very similar to that at 33 GHz (0.224 +/- 0.039). Finally, we have investigated the relationship between variability and spectral index, alpha_13.9^33.75 (where S proportional to nu^-alpha), and find a significant difference in the spectral indices of the variable sources (-0.06 +/- 0.05) and non-variable sources (0.13 +/- 0.04).

We present observations between 1.4 and 18 GHz confirming that G64.5+0.9 is new Galactic shell supernova remnant, using the Very Large Array and the Arcminute Microkelvin Imager. The remnant is a shell ~8 arcmin in diameter with a spectral index of alpha = 0.47 +/- 0.03 (with alpha defined such that flux density S varies with frequency nu as S proportional to nu to the power of -alpha). There is also emission near the centre of the shell, ~1 arcmin in extent, with a spectral index of alpha = 0.81 +/- 0.02. We do not find any evidence for spectral breaks for either source within our frequency range. The nature of the central object is unclear and requires further investigation, but we argue that is most unlikely to be extragalactic. It is difficult to avoid the conclusion that it is associated with the shell, although its spectrum is very unlike that of known pulsar wind nebulae.

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

We present observations between 14.2 and 17.9 GHz of 12 reported supernova remnants (SNRs) made with the Arcminute Microkelvin Imager Small Array (AMI SA). In conjunction with data from the literature at lower radio frequencies, we determine spectra of these objects. For well-studied SNRs (Cas A, Tycho's SNR, 3C58 and the Crab Nebula), the results are in good agreement with spectra based on previous results. For the less well-studied remnants the AMI SA observations provide higher-frequency radio observations than previously available, and better constrain their radio spectra. The AMI SA results confirm a spectral turnover at ~11 GHz for the filled-centre remnant G74.9+1.2. We also see a possible steepening of the spectrum of the filled-centre remnant G54.1+0.3 within the AMI SA frequency band compared with lower frequencies. We confirm that G84.9+0.5, which had previously been identified as a SNR, is rather an HII region and has a flat radio spectrum.

We present observations of the Lynds' dark nebula LDN 1111 made at microwave frequencies between 14.6 and 17.2 GHz with the Arcminute Microkelvin Imager (AMI). We find emission in this frequency band in excess of a thermal free--free spectrum extrapolated from data at 1.4 GHz with matched uv-coverage. This excess is > 15 sigma above the predicted emission. We fit the measured spectrum using the spinning dust model of Drain & Lazarian (1998a) and find the best fitting model parameters agree well with those derived from Scuba data for this object by Visser et al. (2001).

We present a Bayesian approach to modelling galaxy clusters using multi-frequency pointed observations from telescopes that exploit the Sunyaev--Zel'dovich effect. We use the recently developed MultiNest technique (Feroz, Hobson & Bridges, 2008) to explore the high-dimensional parameter spaces and also to calculate the Bayesian evidence. This permits robust parameter estimation as well as model comparison. Tests on simulated Arcminute Microkelvin Imager observations of a cluster, in the presence of primary CMB signal, radio point sources (detected as well as an unresolved background) and receiver noise, show that our algorithm is able to analyse jointly the data from six frequency channels, sample the posterior space of the model and calculate the Bayesian evidence very efficiently on a single processor. We also illustrate the robustness of our detection process by applying it to a field with radio sources and primordial CMB but no cluster, and show that indeed no cluster is identified. The extension of our methodology to the detection and modelling of multiple clusters in multi-frequency SZ survey data will be described in a future work.

The Arcminute Microkelvin Imager is a pair of interferometer arrays operating with six frequency channels spanning 13.9-18.2 GHz, with very high sensitivity to angular scales 30''-10'. The telescope is aimed principally at Sunyaev-Zel'dovich imaging of clusters of galaxies. We discuss the design of the telescope and describe and explain its electronic and mechanical systems.

We present observations between 14.2 and 17.9 GHz of sixteen Galactic HII regions made with the Arcminute Microkelvin Imager (AMI). In conjunction with data from the literature at lower radio frequencies we investigate the possibility of a spinning dust component in the spectra of these objects. We conclude that there is no significant evidence for spinning dust towards these sources and measure an average spectral index of 0.15+/-0.07 between 1.4 and 17.9 GHz for the sample.

We report the first detection of a Sunyaev-Zel'dovich (S-Z) decrement with the Arcminute Microkelvin Imager (AMI). We have made commissioning observations towards the cluster A1914 and have measured an integrated flux density of -8.61 mJy in a uv-tapered map with noise level 0.19 mJy/beam. We find that the spectrum of the decrement, measured in the six channels between 13.5-18GHz, is consistent with that expected for a S-Z effect. The sensitivity of the telescope is consistent with the figures used in our simulations of cluster surveys with AMI.