We derive a model for Sunyaev--Zel'dovich data from a galaxy cluster which uses an Einasto profile to model the cluster's dark matter component. This model is similar to the physical models for clusters previously used by the Arcminute Microkelvin Imager (AMI) consortium, which model the dark matter using a Navarro-Frenk-White (NFW) profile, but the Einasto profile provides an extra degree of freedom. We thus present a comparison between two physical models which differ only in the way they model dark matter: one which uses an NFW profile (PM I) and one that uses an Einasto profile (PM II). We illustrate the differences between the models by plotting physical properties of clusters as a function of cluster radius. We generate AMI simulations of clusters which are \textit{created} and \textit{analysed} with both models. From this we find that for 14 of the 16 simulations, the Bayesian evidence gives no preference to either of the models according to the Jeffreys scale, and for the other two simulations, weak preference in favour of the correct model. However, for the mass estimates obtained from the analyses, the values were within $1\sigma$ of the input values for 14 out of 16 of the clusters when using the correct model, but only in 6 out of 16 cases when the incorrect model was used to analyse the data. Finally we apply the models to real data from cluster A611 obtained with AMI, and find the mass estimates to be consistent with one another except in the case of when PM II is applied using an extreme value for the Einasto shape parameter.

We present a comparison between three cluster models applied to data obtained by the Arcminute Microkelvin Imager radio interferometer system. The physical model (PM) parameterises a cluster in terms of its physical quantities to model the dark matter and baryonic components of the cluster using NFW and GNFW profiles respectively. The observational models (OM I and OM II) model only the gas content of the cluster. The two OMs vary only in the priors they use in Bayesian inference: OM I has a joint prior on angular radius $\theta$ and integrated Comptonisation $Y$, derived from simulations, while OM II uses separable priors on $\theta$ and $Y$ which are based on calculations of the physical model. For the comparison we consider a sample of $54$ clusters which are a subsample of the second Planck catalogue of Sunyaev-Zel'dovich sources. We first compare the $Y$ estimates of the three models, and find that the PM generally yields lower estimates relative to the OMs. We then compute the Earth Mover's Distance between the $\theta$ - $Y$ posterior distributions obtained from each model for each cluster, and find that the two models which are most discrepant are PM and OM I. Finally, we compare the Bayesian evidence values obtained from each model for each cluster. OM I generally provides the best fit to the data but not at a statistically significant level, according to the Jeffreys scale. The highest evidence ratio obtained is actually in favour of the PM over OM I.

We present a comparison of mass estimates for $54$ galaxy cluster candidates from the second Planck catalogue (PSZ2) of Sunyaev-Zel'dovich sources. We compare the mass values obtained with data taken from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite. The former of these uses a Bayesian analysis pipeline that parameterises a cluster in terms of its physical quantities, and models the dark matter & baryonic components of a cluster using NFW and GNFW profiles respectively. Our mass estimates derived from Planck data are obtained from the results of the Bayesian detection algorithm PowellSnakes (PwS), are based on the methodology detailed in the PSZ2 paper, and produce two sets of mass estimates; one estimate is calculated directly from the angular radius $\theta$ - integrated Comptonisation parameter $Y$ posterior distributions, and the other uses a `slicing function' to provide information on $\theta$ based on X-ray measurements and previous Planck mission samples. We find that for $37$ of the clusters, the AMI mass estimates are lower than both values obtained from Planck data. However the AMI and slicing function estimates are within one combined standard deviation of each other for $31$ clusters. We also generate cluster simulations based on the slicing-function mass estimates, and analyse them in the same way as we did the real AMI data. We find that inclusion in the simulations of radio-source confusion & CMB noise and measurable radio-sources causes AMI mass estimates to be systematically low.

We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the analysis to simulated clusters and find that when the cluster pressure profile is known a-priori, the joint dataset provides precise and accurate constraints on the cluster parameters, removing the need for external information to reduce the parameter degeneracy. When the pressure profile deviates from that assumed for the fit, the constraints become biased. Allowing the pressure profile shape parameters to vary in the analysis allows an unbiased recovery of the integrated cluster signal and produces constraints on some shape parameters, depending on the angular size of the cluster. When applied to real data from Planck-detected cluster PSZ2 G063.80+11.42, our method resolves the discrepancy between the AMI and Planck $Y$-estimates and usefully constrains the gas pressure profile shape parameters at intermediate and large radii.

The Arcminute Microkelvin Imager (AMI) telescopes located at the Mullard Radio Astronomy Observatory near Cambridge have been significantly enhanced by the implementation of a new digital correlator with 1.2 MHz spectral resolution. This system has replaced a 750-MHz resolution analogue lag-based correlator, and was designed to mitigate the effects of radio frequency interference, particularly from geostationary satellites that contaminate observations at low declinations. The upgraded instrument consists of 18 ROACH2 Field Programmable Gate Array platforms used to implement a pair of real-time FX correlators -- one for each of AMI's two arrays. The new system separates the down-converted RF baseband signal from each AMI receiver into two 2.3 GHz-wide sub-bands which are each digitized at 5-Gsps with 8 bits of precision. These digital data streams are filtered into 2048 frequency channels and cross-correlated using FPGA hardware, with a commercial 10 Gb Ethernet switch providing high-speed data interconnect. Images formed using data from the new digital correlator show over an order of magnitude improvement in dynamic range over the previous system. The ability to observe at low declinations has also been significantly improved.

The Arcminute Microkelvin Imager (AMI) carried out a blind survey for galaxy clusters via their Sunyaev-Zel'dovich effect decrements between 2008 and 2011. The first detection, known as AMI-CL J0300+2613, has been reobserved with AMI equipped with a new digital correlator with high dynamic range. The combination of the new AMI data and more recent high-resolution sub-mm and infra-red maps now shows the feature in fact to be a ring of positive dust-correlated Galactic emission, which is likely to be anomalous microwave emission (AME). If so, this is the first completely blind detection of AME at arcminute scales.

AMI observations towards CIZA J2242+5301, in comparison with observations of weak gravitational lensing and X-ray emission from the literature, are used to investigate the behaviour of non-baryonic dark matter (NBDM) and gas during the merger. Analysis of the Sunyaev-Zel'dovich (SZ) signal indicates the presence of high pressure gas elongated perpendicularly to the X-ray and weak-lensing morphologies which, given the merger-axis constraints in the literature, implies that high pressure gas is pushed out into a linear structure during core passing. Simulations in the literature closely matching the inferred merger scenario show the formation of gas density and temperature structures perpendicular to the merger axis. These SZ observations are challenging for modified gravity theories in which NBDM is not the dominant contributor to galaxy-cluster gravity.

A new method is presented for modelling the physical properties of galaxy clusters. Our technique moves away from the traditional approach of assuming specific parameterised functional forms for the variation of physical quantities within the cluster, and instead allows for a 'free-form' reconstruction, but one for which the level of complexity is determined automatically by the observational data and may depend on position within the cluster. This is achieved by representing each independent cluster property as some interpolating or approximating function that is specified by a set of control points, or 'nodes', for which the number of nodes, together with their positions and amplitudes, are allowed to vary and are inferred in a Bayesian manner from the data. We illustrate our nodal approach in the case of a spherical cluster by modelling the electron pressure profile Pe(r) in analyses both of simulated Sunyaev-Zel'dovich (SZ) data from the Arcminute MicroKelvin Imager (AMI) and of real AMI observations of the cluster MACS J0744+3927 in the CLASH sample. We demonstrate that one may indeed determine the complexity supported by the data in the reconstructed Pe(r), and that one may constrain two very important quantities in such an analysis: the cluster total volume integrated Comptonisation parameter (Ytot) and the extent of the gas distribution in the cluster (rmax). The approach is also well-suited to detecting clusters in blind SZ surveys.

Using Arcminute Microkelvin Imager (AMI) SZ observations towards ten CLASH clusters we investigate the influence of cluster mergers on observational galaxy cluster studies. Although selected to be largely relaxed, there is disagreement in the literature on the dynamical states of CLASH sample members. We analyse our AMI data in a fully Bayesian way to produce estimated cluster parameters and consider the intrinsic correlations in our NFW/GNFW-based model. Varying pressure profile shape parameters, illustrating an influence of mergers on scaling relations, induces small deviations from the canonical self-similar predictions -- in agreement with simulations of Poole et al. 2007 who found that merger activity causes only small scatter perpendicular to the relations. We demonstrate this effect observationally using the different dependencies of SZ and X-ray signals to $n_{\rm e}$ that cause different sensitivities to the shocking and/or fractionation produced by mergers. Plotting $Y_{\rm X}$--$M_{\rm gas}$ relations (where $Y_{\rm X}=M_{\rm gas}T$) derived from AMI SZ and from $Chandra$ X-ray gives ratios of AMI and $Chandra$ $Y_{\rm X}$ and $M_{\rm gas}$ estimates that indicate movement of clusters \textit{along} the scaling relation, as predicted by Poole et al. 2007. Clusters that have moved most along the relation have the most discrepant $T_{\rm SZ}$ and $T_{\rm X}$ estimates: all the other clusters (apart from one) have SZ and X-ray estimates of $M_{\rm gas}$, $T$ and $Y_{\rm X}$ that agree within $r_{500}$. We use SZ vs X-ray discrepancies in conjunction with $Chandra$ maps and $T_{\rm X}$ profiles, making comparisons with simulated cluster merger maps in Poole et al. 2006, to identify disturbed members of our sample and estimate merger stages.

Radio relics are patches of diffuse synchrotron radio emission that trace shock waves. Relics are thought to form when intra-cluster medium electrons are accelerated by cluster merger induced shock waves through the diffusive shock acceleration mechanism. In this paper, we present observations spanning 150 MHz to 30 GHz of the `Sausage' and `Toothbrush' relics from the Giant Metrewave and Westerbork telescopes, the Karl G. Jansky Very Large Array, the Effelsberg telescope, the Arcminute Microkelvin Imager and Combined Array for Research in Millimeter-wave Astronomy. We detect both relics at 30 GHz, where the previous highest frequency detection was at 16 GHz. The integrated radio spectra of both sources clearly steepen above 2 GHz, at the >6$\sigma$ significance level, supports the spectral steepening previously found in the `Sausage' and the Abell 2256 relic. Our results challenge the widely adopted simple formation mechanism of radio relics and suggest more complicated models have to be developed that, for example, involve re-acceleration of aged seed electrons.

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

We consider the high radio frequency (15 GHz - 353 GHz) properties and variability of 35 Brightest Cluster Galaxies (BCGs). These are the most core-dominated sources drawn from a parent sample of more than 700 X-ray selected clusters, thus allowing us to relate our results to the general population. We find that >6.0% of our parent sample (>15.1% if only cool-core clusters are considered) contain a radio-source at 150 GHz of at least 3mJy (~1x10^23 W/Hz at our median redshift of z~0.13). Furthermore, >3.4% of the BCGs in our parent sample contain a peaked component (Gigahertz Peaked Spectrum, GPS) in their spectra that peaks above 2 GHz, increasing to >8.5% if only cool-core clusters are considered. We see little evidence for strong variability at 15 GHz on short (week-month) time-scales although we see variations greater than 20% at 150 GHz over 6-month times-frames for 4 of the 23 sources with multi-epoch observations. Much more prevalent is long-term (year-decade time-scale) variability, with average annual amplitude variations greater than 1% at 15 GHz being commonplace. There is a weak trend towards higher variability as the peak of the GPS-like component occurs at higher frequency. We demonstrate the complexity that is seen in the radio spectra of BCGs and discuss the potentially significant implications of these high-peaking components for Sunyaev-Zel'dovich cluster searches.

We present observations and analysis of a sample of 123 galaxy clusters from the 2013 Planck catalogue of Sunyaev-Zel'dovich sources with the Arcminute Microkelvin Imager (AMI), a ground-based radio interferometer. AMI provides an independent measurement with higher angular resolution, 3 arcmin compared to the Planck beams of 5-10 arcmin. The AMI observations thus provide validation of the cluster detections, improved positional estimates, and a consistency check on the fitted 'size' ($\theta_{s}$) and 'flux' ($Y_{\rm tot}$) parameters in the Generalised Navarro, Frenk and White (GNFW) model. We detect 99 of the clusters. We use the AMI positional estimates to check the positional estimates and error-bars produced by the Planck algorithms PowellSnakes and MMF3. We find that $Y_{\rm tot}$ values as measured by AMI are biased downwards with respect to the Planck constraints, especially for high Planck-SNR clusters. We perform simulations to show that this can be explained by deviation from the 'universal' pressure profile shape used to model the clusters. We show that AMI data can constrain the $\alpha$ and $\beta$ parameters describing the shape of the profile in the GNFW model for individual clusters provided careful attention is paid to the degeneracies between parameters, but one requires information on a wider range of angular scales than are present in AMI data alone to correctly constrain all parameters simultaneously.

QUIJOTE (Q-U-I JOint TEnerife) is a new polarimeter aimed to characterize the polarization of the Cosmic Microwave Background and other Galactic and extragalactic signals at medium and large angular scales in the frequency range 10-40 GHz. The multi-frequency (10-20~GHz) instrument, mounted on the first QUIJOTE telescope, saw first light on November 2012 from the Teide Observatory (2400~m a.s.l). During 2014 the second telescope has been installed at this observatory. A second instrument at 30~GHz will be ready for commissioning at this telescope during summer 2015, and a third additional instrument at 40~GHz is now being developed. These instruments will have nominal sensitivities to detect the B-mode polarization due to the primordial gravitational-wave component if the tensor-to-scalar ratio is larger than r=0.05.

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

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

We present the first public release of our Bayesian inference tool, Bayes-X, for the analysis of X-ray observations of galaxy clusters. We illustrate the use of Bayes-X by analysing a set of four simulated clusters at z=0.2-0.9 as they would be observed by a Chandra-like X-ray observatory. In both the simulations and the analysis pipeline we assume that the dark matter density follows a spherically-symmetric Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. We then perform four sets of analyses. By numerically exploring the joint probability distribution of the cluster parameters given simulated Chandra-like data, we show that the model and analysis technique can robustly return the simulated cluster input quantities, constrain the cluster physical parameters and reveal the degeneracies among the model parameters and cluster physical parameters. We then analyse Chandra data on the nearby cluster, A262, and derive the cluster physical profiles. To illustrate the performance of the Bayesian model selection, we also carried out analyses assuming an Einasto profile for the matter density and calculated the Bayes factor. The results of the model selection analyses for the simulated data favour the NFW model as expected. However, we find that the Einasto profile is preferred in the analysis of A262. The Bayes-X software, which is implemented in Fortran 90, is available at http://www.mrao.cam.ac.uk/facilities/software/bayesx/.

On 2014 April 23, the Swift satellite detected a gamma-ray superflare from the nearby star system DG CVn. This system comprises a M-dwarf binary with extreme properties: it is very young and at least one of the components is a very rapid rotator. The gamma-ray superflare is one of only a handful detected by Swift in a decade. As part of our AMI-LA Rapid Response Mode, ALARRM, we automatically slewed to this target, were taking data at 15 GHz within six minutes of the burst, and detected a bright (~100 mJy) radio flare. This is the earliest detection of bright, prompt, radio emission from a high energy transient ever made with a radio telescope, and is possibly the most luminous incoherent radio flare ever observed from a red dwarf star. An additional bright radio flare, peaking at around 90 mJy, occurred around one day later, and there may have been further events between 0.1-1 days when we had no radio coverage. The source subsequently returned to a quiescent level of 2-3 mJy on a timescale of about 4 days. Although radio emission is known to be associated with active stars, this is the first detection of large radio flares associated with a gamma ray superflare, and demonstrates both feasibility and scientific importance of rapid response modes on radio telescopes.

The CIZA J2242.8+5301 cluster hosts an extraordinary, narrow, Mpc-wide, diffuse patch of radio emission tracing travelling shock waves, called a relic. We perform a multi-wavelength of observations of the cluster. We discover radio spectral steepening and increasing spectral curvature in the shock downstream area, indicative of emission produced by spectrally-aged electrons, following a JP ageing model. The curved integrated spectrum towards $16$ GHz challenges the favoured relic formation model. Further, we find a boosting the number density of H$\alpha$ emitters by a factor of $10$, possibly caused by interactions of the galaxies with the shock wave.

We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

We observed the cluster CIZA J2242.8+5301 with the Arcminute Microkelvin Imager at $16$ GHz and present the first high radio-frequency detection of diffuse, non-thermal cluster emission. This cluster hosts a variety of bright, extended, steep-spectrum synchrotron-emitting radio sources, associated with the intra-cluster medium, called radio relics. Most notably, the northern, Mpc-wide, narrow relic provides strong evidence for diffusive shock acceleration in clusters. We detect a puzzling, flat-spectrum, diffuse extension of the southern relic, which is not visible in the lower radio-frequency maps. The northern radio relic is unequivocally detected and measures an integrated flux of $1.2\pm0.3$ mJy. While the low-frequency ($<2$ GHz) spectrum of the northern relic is well represented by a power-law, it clearly steepens towards $16$ GHz. This result is inconsistent with diffusive shock acceleration predictions of ageing plasma behind a uniform shock front. The steepening could be caused by an inhomogeneous medium with temperature/density gradients or by lower acceleration efficiencies of high energy electrons. Further modelling is necessary to explain the observed spectrum.

We present one of the best sampled early time light curves of a gamma-ray burst (GRB) at radio wavelengths. Using the Arcminute Mircrokelvin Imager (AMI) we observed GRB 130427A at the central frequency of 15.7 GHz between 0.36 and 59.32 days post-burst. These results yield one of the earliest radio detections of a GRB and demonstrate a clear rise in flux less than one day after the gamma-ray trigger followed by a rapid decline. This early time radio emission probably originates in the GRB reverse shock so our AMI light curve reveals the first ever confirmed detection of a reverse shock peak in the radio domain. At later times (about 3.2 days post-burst) the rate of decline decreases, indicating that the forward shock component has begun to dominate the light-curve. Comparisons of the AMI light curve with modelling conducted by Perley et al. show that the most likely explanation of the early time 15.7 GHz peak is caused by the self-absorption turn-over frequency, rather than the peak frequency, of the reverse shock moving through the observing bands.

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

We have obtained deep SZ observations towards 15 of the apparently hottest XMM Cluster Survey (XCS) clusters that can be observed with the Arcminute Microkelvin Imager (AMI). We use a Bayesian analysis to quantify the significance of our SZ detections. We detect the SZ effect at high significance towards three of the clusters and at lower significance for a further two clusters. Towards the remaining ten clusters, no clear SZ signal was measured. We derive cluster parameters using the XCS mass estimates as a prior in our Bayesian analysis. For all AMI-detected clusters, we calculate large-scale mass and temperature estimates while for all undetected clusters we determine upper limits on these parameters. We find that the large- scale mean temperatures derived from our AMI SZ measurements (and the upper limits from null detections) are substantially lower than the XCS-based core-temperature estimates. For clusters detected in the SZ, the mean temperature is, on average, a factor of 1.4 lower than temperatures from the XCS. For clusters undetected in SZ, the average 68% upper limit on the mean temperature is a factor of 1.9 below the XCS temperature.

We present CARMA observations of a massive galaxy cluster discovered in the AMI blind SZ survey. Without knowledge of the cluster redshift a Bayesian analysis of the AMI, CARMA and joint AMI & CARMA uv-data is used to quantify the detection significance and parameterise both the physical and observational properties of the cluster whilst accounting for the statistics of primary CMB anisotropies, receiver noise and radio sources. The joint analysis of the AMI & CARMA uv-data was performed with two parametric physical cluster models: the {\beta}-model; and the model described in Olamaie et al. 2012 with the pressure profile fixed according to Arnaud et al. 2010. The cluster mass derived from these different models is comparable but our Bayesian evidences indicate a preference for the {\beta}-profile which we, therefore, use throughout our analysis. From the CARMA data alone we obtain a Bayesian probability of detection ratio of 12.8:1 when assuming that a cluster exists within our search area; alternatively assuming that Jenkins et al. 2001 accurately predicts the number of clusters as a function of mass and redshift, the Bayesian probability of detection is 0.29:1. From the analysis of the AMI or AMI & CARMA data the probability of detection ratio exceeds 4.5x10^3:1. Performing a joint analysis of the AMI & CARMA data with a physical cluster model we derive the total mass internal to r200 as MT,200 = 4.1x10^14Msun. Using a phenomenological {\beta}-model to quantify the temperature decrement as a function of angular distance we find a central SZ temperature decrement of 170{\mu}K in the AMI & CARMA data. The SZ decrement in the CARMA data is weaker than expected and we speculate that this is a consequence of the cluster morphology. In a forthcoming study we will assess the impact of cluster morphology on the SZ decrements that are observed with interferometers such as AMI and CARMA.

We present 16-GHz Sunyaev-Zel'dovich observations using the Arcminute Microkelvin Imager (AMI) and subsequent Bayesian analysis of six galaxy clusters at redshift ($z \approx 1$) chosen from an X-ray and Infrared selected sample from Culverhouse et al. (2010). In the subsequent analysis we use two cluster models, an isothermal \beta-model and a Dark Matter GNFW (DM-GNFW) model in order to derive a formal detection probability and the cluster parameters. We detect two clusters (CLJ1415+3612 & XMJ0830+5241) and measure their total masses out to a radius of 200 $\times$ the critical density at the respective cluster's redshift. For CLJ1415+3612 and XMJ0830+5241, we find M_{\mathrm{T},200} for each model, which agree with each other for each cluster. We also present maps before and after source subtraction of the entire sample and provide 1D and 2D posterior marginalised probability distributions for each fitted cluster profile parameter of the detected clusters. Using simulations which take into account the measured source environment from the AMI Large Array (LA), source confusion noise, CMB primordials, instrument noise, we estimate from low-radius X-ray data from Culverhouse et al. (2010), the detectability of each cluster in the sample and compare it with the result from the Small Array (SA) data. Furthermore, we discuss the validity of the assumptions of isothermality and constant gas mass fraction. We comment on the bias that these small-radius estimates introduce to large-radius SZ predictions. In addition, we follow-up the two detections with deep, single-pointed LA observations. We find a 3 sigma tentative decrement toward CLJ1415+3612 at high-resolution and a 5 sigma high-resolution decrement towards XMJ0830+5241.

The AMI Galactic Plane Survey (AMIGPS) is a large area survey of the outer Galactic plane to provide arcminute resolution images at milli-Jansky sensitivity in the centimetre-wave band. Here we present the first data release of the survey, consisting of 868 deg^2 of the Galactic plane, covering the area 76 deg \lessapprox l \lessapprox 170 deg between latitudes of |b| \lessapprox 5 deg, at a central frequency of 15.75 GHz (1.9 cm). We describe in detail the drift scan observations which have been used to construct the maps, including the techniques used for observing, mapping and source extraction, and summarise the properties of the finalized datasets. These observations constitute the most sensitive Galactic plane survey of large extent at centimetre-wave frequencies greater than 1.4 GHz.

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.

We present 16-GHz AMI SZ observations of 19 clusters with L_X >7x10^37 W (h50=1) selected from the LoCuS survey (0.142<z<0.295) and of A1758b, in the FoV of A1758a. We detect 17 clusters with 5-23sigma peak surface brightnesses. Cluster parameters are obtained using a Bayesian cluster analysis. We fit isothermal beta-models to our data and assume the clusters are virialized (with all the kinetic energy in gas internal energy). Our gas temperature, T_AMI, is derived from AMI SZ data, not from X-ray spectroscopy. Cluster parameters internal to r500 are derived assuming HSE. We find: (i) Different gNFW parameterizations yield significantly different parameter degeneracies. (ii) For h70 = 1, we find the virial radius r200 to be typically 1.6+/-0.1 Mpc and the total mass M_T(r200) typically to be 2.0-2.5xM_T(r500).(iii) Where we have found M_T X-ray (X) and weak-lensing (WL) values in the literature, there is good agreement between WL and AMI estimates (with M_{T,AMI}/M_{T,WL} =1.2^{+0.2}_{-0.3} and =1.0+/-0.1 for r500 and r200, respectively). In comparison, most Suzaku/Chandra estimates are higher than for AMI (with M_{T,X}/M_{T,AMI}=1.7+/-0.2 within r500), particularly for the stronger mergers.(iv) Comparison of T_AMI to T_X sheds light on high X-ray masses: even at large r, T_X can substantially exceed T_AMI in mergers. The use of these higher T_X values will give higher X-ray masses. We stress that large-r T_SZ and T_X data are scarce and must be increased. (v) Despite the paucity of data, there is an indication of a relation between merger activity and SZ ellipticity. (vi) At small radius (but away from any cooling flow) the SZ signal (and T_AMI) is less sensitive to ICM disturbance than the X-ray signal (and T_X) and, even at high r, mergers affect n^2-weighted X-ray data more than n-weighted SZ, implying significant shocking or clumping or both occur even in the outer parts of mergers.

We present 16-GHz observations using the Arcminute Microkelvin Imager (AMI) of 11 clusters with 7 x 10^{37}W < L_X < 11 x 10^{37}W (h_{50}=1.0) selected from the Local Cluster Substructure Survey (LoCuSS) and compare them to X-ray data. We use a fast, Bayesian cluster analysis to explore the high-dimensional parameter space of the cluster-plus-sources model and obtain robust cluster parameter estimates in the presence of radio point sources, receiver noise and primordial CMB anisotropy. Our analysis fits a spherical, isothermal beta-model to our data and assumes the cluster follows the theoretical mass-temperature relation. Large-scale cluster parameters internal to r_{500} are derived under the assumption of hydrostatic equilibrium. Posterior distributions for the large-scale parameters of 8 of our clusters are given; SZ effects towards Abell 1704 and Zw0857.9+2107 were not detected and our spherical beta-profile was found to be an inadequate fit to the decrement on our map for Abell 2409.

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

We 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 16 GHz (1.9 cm) deep radio continuum observations made with the Arcminute Microkelvin Imager (AMI) of a sample of low-mass young stars driving jets. We combine these new data with archival information from an extensive literature search to examine spectral energy distributions (SEDs) for each source and calculate both the radio and sub-mm spectral indices in two different scenarios: (1) fixing the dust temperature (Td) according to evolutionary class; (2) allowing Td to vary. We use the results of this analysis to place constraints on the physical mechanisms responsible for the radio emission. From AMI data alone, as well as from model fitting to the full SED in both scenarios, we find that 80 per cent of the objects in this sample have spectral indices consistent with free-free emission. We find an average spectral index in both Td scenarios consistent with free-free emission. We examine correlations of the radio luminosity with bolometric luminosity, envelope mass, and outflow force and find that these data are consistent with the strong correlation with envelope mass seen in lower luminosity samples. We examine the errors associated with determining the radio luminosity and find that the dominant source of error is the uncertainty on the opacity index, beta. We examine the SEDs for variability in these young objects, and find evidence for possible radio flare events in the histories of L1551 IRS 5 and Serpens SMM 1.

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.

We present deep 1.8 cm (16 GHz) radio continuum imaging of seven young stellar objects in the Taurus molecular cloud. These objects have previously been extensively studied in the sub-mm to NIR range and their SEDs modelled to provide reliable physical and geometrical parametres.We use this new data to constrain the properties of the long-wavelength tail of the greybody spectrum, which is expected to be dominated by emission from large dust grains in the protostellar disk. We find spectra consistent with the opacity indices expected for such a population, with an average opacity index of beta = 0.26+/-0.22 indicating grain growth within the disks. We use spectra fitted jointly to radio and sub-mm data to separate the contributions from thermal dust and radio emission at 1.8 cm and derive disk masses directly from the cm-wave dust contribution. We find that disk masses derived from these flux densities under assumptions consistent with the literature are systematically higher than those calculated from sub-mm data, and meet the criteria for giant planet formation in a number of cases.

The Planck Early Release Compact Source Catalogue includes 26 sources with no obvious matches in other radio catalogues (of primarily extragalactic sources). Here we present observations made with the Arcminute Microkelvin Imager Small Array (AMI SA) at 15.75 GHz of the eight of the unmatched sources at declination > +10 degrees. Of the eight, four are detected and are associated with known objects. The other four are not detected with the AMI SA, and are thought to be spurious.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Serpens molecular cloud by the Spitzer c2d programme at a wavelength of 1.8cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). These observations have a resolution of ~30arcsec and an average sensitivity of 19microJy/beam. The targets are predominantly Class I sources, and we find the detection rate for Class I objects in this sample to be low (18%) compared to that of Class 0 objects (67%), consistent with previous works. For detected objects we examine correlations of radio luminosity with bolometric luminosity and envelope mass and find that these data support correlations found by previous samples, but do not show any indiction of the evolutionary divide hinted at by similar data from the Perseus molecular cloud when comparing radio luminosity with envelope mass. We conclude that envelope mass provides a better indicator for radio luminosity than bolometric luminosity, based on the distribution of deviations from the two correlations. Combining these new data with archival 3.6cm flux densities we also examine the spectral indices of these objects and find an average spectral index of 0.53+/-1.14, consistent with the canonical value for a partially optically thick spherical or collimated stellar wind. However, we caution that possible inter-epoch variability limits the usefulness of this value, and such variability is supported by our identification of a possible flare in the radio history of Serpens SMM 1.

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

We present follow-up observations of two galaxy clusters detected blindly via the Sunyaev-Zel'dovich (SZ) effect and released in the Planck Early Release Compact Source Catalogue. We use the Arcminute Microkelvin Imager, a dual-array 14-18 GHz radio interferometer. After radio source subtraction, we find a SZ decrement of integrated flux density -1.08+/-0.10 mJy toward PLCKESZ G121.11+57.01, and improve the position measurement of the cluster, finding the centre to be RA 12 59 36.4, Dec +60 04 46.8, to an accuracy of 20 arcseconds. The region of PLCKESZ G115.71+17.52 contains strong extended emission, so we are unable to confirm the presence of this cluster via the SZ effect.

We present deep radio continuum observations of the cores identified as deeply embedded young stellar objects in the Perseus molecular cloud by the Spitzer c2d programme at a wavelength of 1.8 cm with the Arcminute Microkelvin Imager Large Array (AMI-LA). We detect 72% of Class 0 objects from this sample and 31% of Class I objects. No starless cores are detected. We use the flux densities measured from these data to improve constraints on the correlations between radio luminosity and bolometric luminosity, infrared luminosity and, where measured, outflow force. We discuss the differing behaviour of these objects as a function of protostellar class and investigate the differences in radio emission as a function of core mass. Two of four possible very low luminosity objects (VeLLOs) are detected at 1.8 cm.

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 present observations from the Small Array of the Arcminute Microkelvin Imager (AMI) of eight high X-ray luminosity galaxy cluster systems selected from the Local Cluster Substructure Survey (LoCuSS) sample.We detect the Sunyaev-Zel'dovich (SZ) effect in seven of these clusters. With the assumptions that galaxy clusters are isothermal, have a density profile described by a spherical b -model and obey the theoretical M-T relation, we are able to derive cluster parameters at r200 from our SZ data. With the additional assumption of hydrostatic equilibrium we are able to derive parameters at r500. We present posterior probability distributions for cluster parameters such as mass, radius and temperature (TSZ, MT). Combining our sample with that of AMI Consortium: Rodr'\iguez-Gonz\'alvez et al. (2011) and using large-radius X-ray temperature estimates (TX) from Chandra and Suzaku observations, we find that there is reasonable correspondence between TX and TSZ,MT values at low TX, but that for clusters with TX above around 6keV the correspondence breaks down with TX exceeding TSZ, MT; we stress that this finding is based on just ten clusters.

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.

The Very Small Array (VSA) has been used to survey the l = 27 to 46 deg, |b|<4 deg region of the Galactic plane at a resolution of 13 arcmin. The survey consists of 44 pointings of the VSA, each with a r.m.s. sensitivity of ~90 mJy/beam. These data are combined in a mosaic to produce a map of the area. The majority of the sources within the map are HII regions. We investigated anomalous radio emission from the warm dust in 9 HII regions of the survey by making spectra extending from GHz frequencies to the FIR IRAS frequencies. Acillary radio data at 1.4, 2.7, 4.85, 8.35, 10.55, 14.35 and 94 GHz in addition to the 100, 60, 25 and 12 micron IRAS bands were used to construct the spectra. From each spectrum the free-free, thermal dust and anomalous dust emission were determined for each HII region. The mean ratio of 33 GHz anomalous flux density to FIR 100 micron flux density for the 9 selected HII regions was 1.10 +/-0.21x10^(-4). When combined with 6 HII regions previously observed with the VSA and the CBI, the anomalous emission from warm dust in HII regions is detected with a 33 GHz emissivity of 4.65 +/- 0.4 micro K/ (MJy/sr) at 11.5{\sigma}. The anomalous radio emission in HII regions is on average 41+/-10 per cent of the radio continuum at 33 GHz.

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

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of $1.84\cdot 10^{10} M_\odot$, the secondary mass $1.46\cdot 10^{8} M_\odot$, major axis precession rate $39^\circ.1$ per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is $0.28\:\pm\:0.01$ (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the Post Newtonian equations of motion does not exactly follow Einstein's theory: a parameter $q$ is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein's value $q = 1$. This supports the $no-hair theorem$ of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein's value with the accuracy of 2% (1 sigma).

The compact binary system in OJ287 is modelled to contain a spinning primary black hole with an accretion disk and a non-spinning secondary black hole. Using Post Newtonian (PN) accurate equations that include 2.5PN accurate non-spinning contributions, the leading order general relativistic and classical spin-orbit terms, the orbit of the binary black hole in OJ287 is calculated and as expected it depends on the spin of the primary black hole. Using the orbital solution, the specific times when the orbit of the secondary crosses the accretion disk of the primary are evaluated such that the record of observed outbursts from 1913 up to 2007 is reproduced. The timings of the outbursts are quite sensitive to the spin value. In order to reproduce all the known outbursts, including a newly discovered one in 1957, the Kerr parameter of the primary has to be $0.28 \pm 0.08$. The quadrupole-moment contributions to the equations of motion allow us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ where 0.3 is the one sigma error. This supports the `black hole no-hair theorem' within the achievable precision. It should be possible to test the present estimate in 2015 when the next outburst is due. The timing of the 2015 outburst is a strong function of the spin: if the spin is 0.36 of the maximal value allowed in general relativity, the outburst begins in early November 2015, while the same event starts in the end of January 2016 if the spin is 0.2

The dust feature G159.6--18.5 in the Perseus region has previously been observed with the COSMOSOMAS experiment \citep{Watson:05} on angular scales of $\approx$ 1$^{\circ}$, and was found to exhibit anomalous microwave emission. We present new observations of this dust feature, performed with the Very Small Array (VSA) at 33 GHz, to help increase the understanding of the nature of this anomalous emission. On the angular scales observed with the VSA ($\approx$ 10 -- 40$^{\prime}$), G159.6--18.5 consists of five distinct components, each of which have been individually analysed. All five of these components are found to exhibit an excess of emission at 33 GHz, and are found to be highly correlated with far-infrared emission. We provide evidence that each of these compact components have anomalous emission that is consistent with electric dipole emission from very small, rapidly rotating dust grains. These components contribute $\approx$ 10 % to the flux density of the diffuse extended emission detected by COSMOSOMAS, and are found to have a similar radio emissivity.

Supermassive black holes are common in centers of galaxies. Among the active galaxies, quasars are the most extreme, and their black hole masses range as high as to $6\cdot 10^{10} M_\odot$. Binary black holes are of special interest but so far OJ287 is the only confirmed case with known orbital elements. In OJ287, the binary nature is confirmed by periodic radiation pulses. The period is twelve years with two pulses per period. The last four pulses have been correctly predicted with the accuracy of few weeks, the latest in 2007 with the accuracy of one day. This accuracy is high enough that one may test the higher order terms in the Post Newtonian approximation to General Relativity. The precession rate per period is $39^\circ.1 \pm 0^\circ.1$, by far the largest rate in any known binary, and the $(1.83\pm 0.01)\cdot 10^{10} M_\odot$ primary is among the dozen biggest black holes known. We will discuss the various Post Newtonian terms and their effect on the orbit solution. The over 100 year data base of optical variations in OJ287 puts limits on these terms and thus tests the ability of Einstein's General Relativity to describe, for the first time, dynamic binary black hole spacetime in the strong field regime. The quadrupole-moment contributions to the equations of motion allows us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ which supports the black hole no-hair theorem within the achievable precision.

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.

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.

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

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.

Accurate calibration of data is essential for the current generation of CMB experiments. Using data from the Very Small Array (VSA), we describe procedures which will lead to an accuracy of 1 percent or better for experiments such as the VSA and CBI. Particular attention is paid to the stability of the receiver systems, the quality of the site and frequent observations of reference sources. At 30 GHz the careful correction for atmospheric emission and absorption is shown to be essential for achieving 1 percent precision. The sources for which a 1 percent relative flux density calibration was achieved included Cas A, Cyg A, Tau A and NGC7027 and the planets Venus, Jupiter and Saturn. A flux density, or brightness temperature in the case of the planets, was derived at 33 GHz relative to Jupiter which was adopted as the fundamental calibrator. A spectral index at ~30 GHz is given for each. Cas A,Tau A, NGC7027 and Venus were examined for variability. Cas A was found to be decreasing at $0.394 \pm 0.019$ percent per year over the period March 2001 to August 2004. In the same period Tau A was decreasing at $0.22\pm 0.07$ percent per year. A survey of the published data showed that the planetary nebula NGC7027 decreased at $0.16\pm 0.04$ percent per year over the period 1967 to 2003. Venus showed an insignificant ($1.5 \pm 1.3$ percent) variation with Venusian illumination. The integrated polarization of Tau A at 33 GHz was found to be $7.8\pm 0.6$ percent at pa $ = 148^\circ \pm 3^\circ$.}

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

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

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 present results from observations made at 33 GHz with the Very Small Array (VSA) telescope towards potential candidates in the Galactic plane for spinning dust emission. In the cases of the diffuse HII regions LPH96 and NRAO591 we find no evidence for anomalous emission and, in combination with Effelsberg data at 1.4 and 2.7 GHz, confirm that their spectra are consistent with optically thin free--free emission. In the case of the infra-red bright SNR 3C396 we find emission inconsistent with a purely non-thermal spectrum and discuss the possibility of this excess arising from either a spinning dust component or a shallow spectrum PWN, although we conclude that the second case is unlikely given the strong constraints available from lower frequency radio images.

We have observed the Corona Borealis Supercluster with the Millimeter and Infrared Testa grigia Observatory (MITO), located in the Italian Alps, at 143, 214, 272, and 353 GHz. We present a description of the measurements, data analysis, and results of the observations together with a comparison with observations performed at 33 GHz with the Very Small Array (VSA) interferometer situated at the Teide Observatory (Tenerife-Spain). Observations have been made in the direction of the supercluster towards one Cosmic Microwave Background (CMB) cold spot previously detected in a VSA temperature map. Observational strategy and data analysis are described in detail, explaining the procedures used to disentangle primary and secondary anisotropies in the resulting maps. From a first level of data analysis we find evidence in MITO data of primary anisotropy, however still with room for the presence of secondary anisotropy, especially when VSA results are included. With a second level of data analysis using map-making and the maximum entropy method we claim a weak detection of a faint signal compatible with a SZ effect characterized at most by a Comptonization parameter $y=(7.8^{+5.3}_{-4.4})\times10^{-6}$ 68% CL. The low level of confidence in the presence of a SZ signal invite us to study this sky region with higher sensitivity and angular resolution experiments like the already planned upgraded versions of VSA and MITO.

(Abridged) We have used the Rayner & Best (1989) smooth tests of goodness-of-fit to study the Gaussianity of the Very Small Array (VSA) data. Out of the 41 published VSA individual pointings dedicated to cosmological observations, 37 are found to be consistent with Gaussianity, whereas four pointings show deviations from Gaussianity. In two of them, these deviations can be explained as residual systematic effects of a few visibility points which, when corrected, have a negligible impact on the angular power spectrum. The non-Gaussianity found in the other two (adjacent) pointings seems to be associated to a local deviation of the power spectrum of these fields with respect to the common power spectrum of the complete data set, at angular scales of the third acoustic peak (l = 700-900). No evidence of residual systematics is found in this case, and unsubstracted point sources are not a plausible explanation either. If those visibilities are removed, a cosmological analysis based on this new VSA power spectrum alone shows no differences in the parameter constraints with respect to our published results, except for the physical baryon density, which decreases by 10 percent. Finally, the method has been also used to analyse the VSA observations in the Corona Borealis supercluster region (Genova-Santos et al. 2005), which show a strong decrement which cannot be explained as primordial CMB. Our method finds a clear deviation (99.82%) with respect to Gaussianity in the second-order moment of the distribution, and which can not be explained as systematic effects. A detailed study shows that the non-Gaussianity is produced in scales of l~500, and that this deviation is intrinsic to the data (in the sense that can not be explained in terms of a Gaussian field with a different power spectrum).

We present coincident observations of the Cosmic Microwave Background (CMB) from the Very Small Array (VSA) and Cosmic Background Imager (CBI) telescopes. The consistency of the full datasets is tested in the map plane and the Fourier plane, prior to the usual compression of CMB data into flat bandpowers. Of the three mosaics observed by each group, two are found to be in excellent agreement. In the third mosaic, there is a 2 sigma discrepancy between the correlation of the data and the level expected from Monte Carlo simulations. This is shown to be consistent with increased phase calibration errors on VSA data during summer observations. We also consider the parameter estimation method of each group. The key difference is the use of the variance window function in place of the bandpower window function, an approximation used by the VSA group. A re-evaluation of the VSA parameter estimates, using bandpower windows, shows that the two methods yield consistent results.

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.

We present interferometric imaging at 33 GHz of the Corona Borealis supercluster, using the extended configuration of the Very Small Array. A total area of 24 deg^2 has been imaged, with an angular resolution of 11 arcmin and a sensitivity of 12 mJy/beam. The aim of these observations is to search for Sunyaev-Zel'dovich (SZ) detections from known clusters of galaxies in this supercluster and for a possible extended SZ decrement due to diffuse warm/hot gas in the intercluster medium. We measure negative flux values in the positions of the ten richest clusters in the region. Collectively, this implies a 3.0-sigma detection of the SZ effect. In the clusters A2061 and A2065 we find decrements of approximately 2-sigma. Our main result is the detection of two strong and resolved negative features at -70+-12 mJy/beam (-157+-27 microK) and -103+-10 mJy/beam (-230+-23 microK), respectively, located in a region with no known clusters, near the centre of the supercluster. We discuss their possible origins in terms of primordial CMB anisotropies and/or SZ signals related to either unknown clusters or to a diffuse extended warm/hot gas distribution. Our analyses have revealed that a primordial CMB fluctuation is a plausible explanation for the weaker feature (probability of 37.82%). For the stronger one, neither primordial CMB (probability of 0.33%) nor SZ can account alone for its size and total intensity. The most reasonable explanation, then, is a combination of both primordial CMB and SZ signal. Finally, we explore what characteristics would be required for a filamentary structure consisting of warm/hot diffuse gas in order to produce a significant contribution to such a spot taking into account the constraints set by X-ray data.

We describe the source subtraction strategy and observations for the extended Very Small Array, a CMB interferometer operating at 33 GHz. A total of 453 sources were monitored at 33 GHz using a dedicated source subtraction baseline. 131 sources brighter than 20 mJy were directly subtracted from the VSA visibility data. Some characteristics of the subtracted sources, such as spectra and variability, are discussed. The 33-GHz source counts are estimated from a sample selected at 15 GHz. The selection of VSA fields in order to avoid bright sources introduces a bias into the observed counts. This bias is corrected and the resulting source count is estimated to be complete in the flux-density range 20-114 mJy. The 33-GHz source counts are used to calculate a correction to the VSA power spectrum for sources below the subtraction limit.

We present VSA observations (~34GHz) on scales ~20 arcmin towards a complete, X-ray-flux-limited sample of seven clusters at redshift z<0.1. Four have significant SZ detections in the presence of CMB primordial anisotropy. We use a bayesian MCMC method for inference from the VSA data, with X-ray priors on cluster positions and temperatures, and radio priors on sources. We make assumptions of beta-model gas distributions and of hydrostatic equilibrium, to evaluate probability densities for the gas mass and total mass out to r_200. Our combined estimate of the gas fraction is 0.08^{+0.06}_{-0.04}h^{-1} The random errors are poor (note that the errors are higher than would have been obtained with the usual chi-squared method) but the control of bias is good. We have described the MCMC analysis method specifically in terms of SZ but hope the description will be of more general use. We find that the effects of primordial CMB contamination tend to be similar in the estimates of both the gas mass and total mass over our narrow range of angular scales, so that there is little effect of primordials on the gas fraction determination. Using our total mass estimates we find a normalisation of the mass-temperature relation based on the profiles from the VSA cluster pressure maps that is in good agreement with recent M-T determinations from X-ray cluster measurements.

We present deep Ka-band ($\nu \approx 33$ GHz) observations of the CMB made with the extended Very Small Array (VSA). This configuration produces a naturally weighted synthesized FWHM beamwidth of $\sim 11$ arcmin which covers an $\ell$-range of 300 to 1500. On these scales, foreground extragalactic sources can be a major source of contamination to the CMB anisotropy. This problem has been alleviated by identifying sources at 15 GHz with the Ryle Telescope and then monitoring these sources at 33 GHz using a single baseline interferometer co-located with the VSA. Sources with flux densities $\gtsim 20$ mJy at 33 GHz are subtracted from the data. In addition, we calculate a statistical correction for the small residual contribution from weaker sources that are below the detection limit of the survey. The CMB power spectrum corrected for Galactic foregrounds and extragalactic point sources is presented. A total $\ell$-range of 150-1500 is achieved by combining the complete extended array data with earlier VSA data in a compact configuration. Our resolution of $\Delta \ell \approx 60$ allows the first 3 acoustic peaks to be clearly delineated. The is achieved by using mosaiced observations in 7 regions covering a total area of 82 sq. degrees. There is good agreement with WMAP data up to $\ell=700$ where WMAP data run out of resolution. For higher $\ell$-values out to $\ell = 1500$, the agreement in power spectrum amplitudes with other experiments is also very good despite differences in frequency and observing technique.

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.

We estimate cosmological parameters using data obtained by the Very Small Array (VSA) in its extended configuration, in conjunction with a variety of other CMB data and external priors. Within the flat $\Lambda$CDM model, we find that the inclusion of high resolution data from the VSA modifies the limits on the cosmological parameters as compared to those suggested by WMAP alone, while still remaining compatible with their estimates. We find that $\Omega_{\rm b}h^2=0.0234^{+0.0012}_{-0.0014}$, $\Omega_{\rm dm}h^2=0.111^{+0.014}_{-0.016}$, $h=0.73^{+0.09}_{-0.05}$, $n_{\rm S}=0.97^{+0.06}_{-0.03}$, $10^{10}A_{\rm S}=23^{+7}_{-3}$ and $\tau=0.14^{+0.14}_{-0.07}$ for WMAP and VSA when no external prior is included.On extending the model to include a running spectral index of density fluctuations, we find that the inclusion of VSA data leads to a negative running at a level of more than 95% confidence ($n_{\rm run}=-0.069\pm 0.032$), something which is not significantly changed by the inclusion of a stringent prior on the Hubble constant. Inclusion of prior information from the 2dF galaxy redshift survey reduces the significance of the result by constraining the value of $\Omega_{\rm m}$. We discuss the veracity of this result in the context of various systematic effects and also a broken spectral index model. We also constrain the fraction of neutrinos and find that $f_{\nu}< 0.087$ at 95% confidence which corresponds to $m_\nu<0.32{\rm eV}$ when all neutrino masses are the equal. Finally, we consider the global best fit within a general cosmological model with 12 parameters and find consistency with other analyses available in the literature. The evidence for $n_{\rm run}<0$ is only marginal within this model.

We have tested Very Small Array (VSA) observations of three regions of sky for the presence of non-Gaussianity, using high-order cumulants, Minkowski functionals, a wavelet-based test and a Bayesian joint power spectrum/non-Gaussianity analysis. We find the data from two regions to be consistent with Gaussianity. In the third region, we obtain a 96.7% detection of non-Gaussianity using the wavelet test. We perform simulations to characterise the tests, and conclude that this is consistent with expected residual point source contamination. There is therefore no evidence that this detection is of cosmological origin. Our simulations show that the tests would be sensitive to any residual point sources above the data's source subtraction level of 20 mJy. The tests are also sensitive to cosmic string networks at an rms fluctuation level of $105 \mu K$ (i.e. equivalent to the best-fit observed value). They are not sensitive to string-induced fluctuations if an equal rms of Gaussian CDM fluctuations is added, thereby reducing the fluctuations due to the strings network to $74 \mu K$ rms . We especially highlight the usefulness of non-Gaussianity testing in eliminating systematic effects from our data.

We investigate the constraints on basic cosmological parameters set by the first compact-configuration observations of the Very Small Array (VSA), and other cosmological data sets, in the standard inflationary LambdaCDM model. Using a weak prior 40 < H_0 < 90 km/s/Mpc and 0 < tau < 0.5 we find that the VSA and COBE_DMR data alone produce the constraints Omega_tot = 1.03^{+0.12}_{-0.12}, Omega_bh^2 = 0.029^{+0.009}_{-0.009}, Omega_cdm h^2 = 0.13^{+0.08}_{-0.05} and n_s = 1.04^{+0.11}_{-0.08} at the 68 per cent confidence level. Adding in the type Ia supernovae constraints, we additionally find Omega_m = 0.32^{+0.09}_{-0.06} and Omega_Lambda = 0.71^{+0.07}_{-0.07}. These constraints are consistent with those found by the BOOMERanG, DASI and MAXIMA experiments. We also find that, by combining all the recent CMB experiments and assuming the HST key project limits for H_0 (for which the X-ray plus Sunyaev--Zel'dovich route gives a similar result), we obtain the tight constraints Omega_m=0.28^{+0.14}_{-0.07} and Omega_Lambda= 0.72^{+0.07}_{-0.13}, which are consistent with, but independent of, those obtained using the supernovae data.

The fields chosen for the first observations of the cosmic microwave background with the Very Small Array have been surveyed with the Ryle Telescope at 15 GHz. We have covered three regions around RA 00h20m Dec +30deg, RA 09h40m Dec +32deg and RA 15h40m Dec +43deg (J2000), an area of 520 deg^2. There are 465 sources above the current completeness limit of approximately 25 mJy, although a total of ~760 sources have been detected, some as faint as 10 mJy. This paper describes our techniques for observation and data analysis; it also includes source counts and some discussion of spectra and variability. Preliminary source lists are presented.

We present new Ryle Telescope (RT) observations of the Sunyaev Zel'dovich (SZ) decrement from the cluster Abell 773. The field contains a number of faint radio sources that required careful subtraction. We use ASCA observations to measure the gas temperature and a ROSAT HRI image to model the gas distribution. Normalising the gas distribution to fit the RT visibilities returns a value of H_0 of 77 (+19,-15) km/s/Mpc (1-sigma errors) for an Einstein-de-Sitter universe, or 85 (+20,-17) km/s/Mpc for a flat model with Omega_Lambda = 0.7. The errors quoted include estimates of the effects of the principal errors: noise in the SZ measurement, gas temperature uncertainty, and line-of sight depth uncertainty.

We have observed the cosmic microwave background (CMB) in three regions of sky using the Very Small Array (VSA) in an extended configuration with antennas of beamwidth 2 degrees at 34 GHz. Combined with data from previous VSA observations using a more compact array with larger beamwidth, we measure the power spectrum of the primordial CMB anisotropies between angular multipoles l = 160 - 1400. Such measurements at high l are vital for breaking degeneracies in parameter estimation from the CMB power spectrum and other cosmological data. The power spectrum clearly resolves the first three acoustic peaks, shows the expected fall off in power at high l and starts to constrain the position and height of a fourth peak.

The Very Small Array (VSA) is a synthesis telescope designed to image faint structures in the cosmic microwave background on degree and sub-degree angular scales. The VSA has key differences from other CMB interferometers with the result that different systematic errors are expected. We have tested the operation of the VSA with a variety of blank-field and calibrator observations and cross-checked its calibration scale against independent measurements. We find that systematic effects can be suppressed below the thermal noise level in long observations; the overall calibration accuracy of the flux density scale is 3.5 percent and is limited by the external absolute calibration scale.

We present the power spectrum of the fluctuations in the cosmic microwave background detected by the Very Small Array (VSA) in its first season of observations in its compact configuration. We find clear detections of first and second acoustic peaks at l~200 and l~550, plus detection of power on scales up to l=800. The VSA power spectrum is in very good agreement with the results of the Boomerang, Dasi and Maxima telescopes despite the differing potential systematic errors.

We constrain the basic comological parameters using the first observations by the Very Small Array (VSA) in its extended configuration, together with existing cosmic microwave background data and other cosmological observations. We estimate cosmological parameters for four different models of increasing complexity. In each case, careful consideration is given to implied priors and the Bayesian evidence is calculated in order to perform model selection. We find that the data are most convincingly explained by a simple flat Lambda-CDM cosmology without tensor modes. In this case, combining just the VSA and COBE data sets yields the 68 per cent confidence intervals Omega_b h^2=0.034 (+0.007, -0.007), Omega_dm h^2 = 0.18 (+0.06, -0.04), h=0.72 (+0.15,-0.13), n_s=1.07 (+0.06,-0.06) and sigma_8=1.17 (+0.25, -0.20). The most general model considered includes spatial curvature, tensor modes, massive neutrinos and a parameterised equation of state for the dark energy. In this case, by combining all recent cosmological data, we find, in particular, 95 percent limit on the tensor-to-scalar ratio R < 0.63 and on the fraction of massive neutrinos f_nu < 0.11; we also obtain the 68 per cent confidence interval w=-1.06 (+0.20, -0.25) on the equation of state of dark energy.

We have observed the cosmic microwave background temperature fluctuations in eight fields covering three separated areas of sky with the Very Small Array at 34 GHz. A total area of 101 square degrees has been imaged, with sensitivity on angular scales 3.6 - 0.4 degrees (equivalent to angular multipoles l=150-900). We describe the field selection and observing strategy for these observations. In the full-resolution images (with synthesised beam of FWHM ~ 17 arcmin) the thermal noise is typically 45 microK and the CMB signal typically 55 microK. The noise levels in each field agree well with the expected thermal noise level of the telescope, and there is no evidence of any residual systematic features. The same CMB features are detected in separate, overlapping observations. Discrete radio sources have been detected using a separate 15 GHz survey and their effects removed using pointed follow-up observations at 34 GHz. We estimate that the residual confusion noise due to unsubtracted radio sources is less than 14 mJy/beam (15 microK in the full-resolution images), which added in quadrature to the thermal noise increases the noise level by 6 %. We estimate that the rms contribution to the images from diffuse Galactic emission is less than 6 microK. We also present images which are convolved to maximise the signal-to-noise of the CMB features and are co-added in overlapping areas, in which the signal-to-noise of some individual CMB features exceeds 8.

Interferometers require accurate determination of the array configuration in order to produce reliable observations. A method is presented for finding the maximum-likelihood estimate of the telescope geometry, and of other instrumental parameters, astrometrically from the visibility timelines obtained from observations of celestial calibrator sources. The method copes systematically with complicated and unconventional antenna and array geometries, with electronic bandpasses that are different for each antenna radiometer, and with low signal-to-noise ratios for the calibrators. The technique automatically focusses on the geometry errors that are most significant for astronomical observation. We apply this method to observations made with the Very Small Array and constrain some 450 telescope parameters, such as the antenna positions, effective observing frequencies and correlator amplitudes and phase shifts; this requires only ~ 1 h of CPU time on a typical workstation.

We assess the global properties of associated C IV 1548,1550 absorption lines measured in the spectra of radio-loud quasars drawn from a near-complete, low-frequency selected sample. The observations span restframe C IV in two redshift ranges -- 0.7<z<1.0 and 1.5<z<3.0 -- which were targetted in the UV with HST/STIS and in the optical with ground-based telescopes, respectively. First, we corroborate trends for C IV associated absorption to be found preferentially in steep-spectrum and lobe-dominated quasars, implying the absorbing material tends to lie away from the radio-jet axis. Furthermore, we find a clear anticorrelation between C IV absorption strength and the projected linear size of steep-spectrum quasars, indicative of an evolutionary sequence. We also find that heavily-absorbed quasars are systematically redder, implying dust is mixed in with the C IV-absorbing gas. No redshift dependence was found in any of the trends considered. These new results show that radio sources are triggered in galaxies which are exceptionally rich in gas and dust, which then dissipates on a timescale comparable with, but less than, that of the radio source. This observational sequence, together with the lack of redshift-dependence, points to a direct causal link between the event which triggered the radio source and the build-up of absorbing gas and dust, whose make-up is tantalisingly similar to the products of a post-merger starburst. Thus, these new results provide direct evidence for the clearing of absorbing material around quasars with time, as well as the probable association of starburst activity with the onset of the radio activity in AGN.

We present the results of radio, optical and near-infrared observations of the field of TOC J0233.3+3021, a cluster of milliJansky radiosources from the TexOx Cluster survey. In an observation of this field with the Ryle Telescope (RT) at 15 GHz, we measure a decrement in the cosmic microwave background (CMB) of $-675 \pm 95 \mu$Jy on the RT's $\approx$ 0.65 k$\lambda$ baseline. Using optical and infrared imaging with the McDonald 2.7-m Smith Reflector, Calar Alto 3.5-m telescope and UKIRT, we identify the host galaxies of five of the radiosources and measure magnitudes of $R \approx 24$, $J \approx 20$, $K \approx 18$. The CMB decrement is consistent with the Sunyaev-Zel'dovich (SZ) effect of a massive cluster of galaxies, which if modelled as a spherical King profile of core radius $\theta_C = 20^{\prime\prime}$ has a central temperature decrement of $900 \mu$K. The magnitudes and colours of the galaxies are consistent with those of old ellipticals at $z \sim 1$. We therefore conclude that TOC J0233.3+3021 is a massive, high redshift cluster. These observations add to the growing evidence for a significant population of massive clusters at high redshift, and demonstrate the effectiveness of combining searches for AGN `signposts' to clusters with the redshift-independence of the SZ effect.

(Abridged) This is the second in a series of three papers which present and interpret basic observational data on the 6C* 151-MHz radio sample: a low-frequency selected sample which exploits filtering criteria based on radio properties (steep spectral index and small angular size) to find radio sources at redshift z > 4 within a 0.133sr patch of sky. We present results of a programme of optical spectroscopy which has yielded redshifts in the range 0.5 < z < 4.4 for the 29 sources in the sample, all but six of which are secure. We find that the fil tering criteria used for 6C* are very effective in excluding the low-redshift, low-luminosity radio sources: the median redshift of 6C* is z~1.9 compared to z~1.1 for a complete sample matched in 151-MHz flux density. By combining the emission-line dataset for the 6C* radio sources with those for the 3CRR, 6CE and 7CRS samples we establish that z > 1.75 radio galaxies follow a rough proportionality between Lyalpha- and 151 MHz-luminosity which, like similar correlations seen in samples of lower-redshift radio sources, are indicative of a primary link between the power in the source of the photoionising photons (most likely a hidden quasar nucleus) and the power carried by the radio jets. We argue that radio sources modify their environments and that the range of emission-line properties seen is determined more by the range of source age than by the range in ambient environment. This is in accord with the idea that all high-redshift, high-luminosity radio sources are triggered in similar environments, presumably recently collapsed massive structures.

We describe our algorithm for measuring the Hubble constant from Ryle Telescope (RT) interferometric observations of the Sunyaev-Zel'dovich (SZ) effect from a galaxy cluster and observation of the cluster X-ray emission. We analyse the error budget in this method: as well as radio and X-ray random errors, we consider the effects of clumping and temperature differences in the cluster gas, of the kinetic SZ effect, of bremsstrahlung emission at radio wavelengths, of the gravitational lensing of background radio sources and of primary calibration error. Using RT, ASCA and ROSAT observations of the cluster Abell 1413, we find that random errors dominate over systematic ones, and estimate H_0 = 57^{+23}_{-16} km/s/Mpc (1-sigma errors).

We discuss prospects for cluster detection via the Sunyaev-Zel'dovich (SZ) effect in a blank field survey with the proposed interferometer array, the Arcminute MicroKelvin Imager (AMI). Clusters of galaxies selected in the SZ effect probe cosmology and structure formation with little observational bias, because the effect directly measures integrated gas pressure, and does so independently of cluster redshift. We use hydrodynamical simulations in combination with the Press-Schechter expression to simulate SZ cluster sky maps. These are used with simulations of the observation process to gauge the expected SZ cluster counts. Even with a very conservative choice of parameters we find that AMI will discover at least several tens of clusters every year with total mass down to 10^14 solar masses; the numbers depends on factors such as the mean matter density, the density fluctuation power spectrum and cluster gas evolution. The AMI survey itself can distinguish between these to some degree, and parameter degeneracies are largely eliminated given optical and X-ray follow-up of these clusters; this will also permit direct investigation of cluster physics and what drives the evolution.

We have observed the Sunyaev-Zel'dovich effect in a sample of five moderate-redshift clusters with the Ryle Telescope, and used them in conjunction with X-ray imaging and spectral data from ROSAT and ASCA to measure the Hubble constant. This sample was chosen with a strict X-ray flux limit using both the BCS and NORAS cluster catalogues to be well above the surface-brightness limit of the ROSAT All-Sky Survey, and hence to be unbiased with respect to the orientation of the cluster. This controls the major potential systematic effect in the SZ/X-ray method of measureing H0. Taking the weighted geometric mean of the results and including the main sources of random error, namely the noise in the SZ measurement, the uncertainty in the X-ray temperatures and the unknown ellipticity of the clusters, we find H0 = 59 +8/-7 km/s/Mpc assuming a standard CDM model with Omega_M = 1.0, Omega_Lambda = 0.0, or H0 = 65 +8/-7 km/s/Mpc if Omega_M = 0.3, Omega_Lambda = 0.7.

We have mapped the high-redshift (z = 0.546) cluster Cl 0016+16 with the Ryle Telescope at 15 GHz. The Sunyaev-Zel'dovich decrement is clearly detected, and resolved. We combine our data with an X-ray image from ROSAT, and a gas temperature from ASCA to estimate the Hubble Constant H0 = 69 +21/-16 km/s/Mpc for an Omega_M=1.0 cosmology or H0 = 84 +25/-19 km/s/Mpc for Omega_M=0.3 and Omega_Lambda=0.7.

We describe a maximum-likelihood technique for the removal of contaminating radio sources from interferometric observations of the Sunyaev-Zel'dovich (SZ) effect. This technique, based on a simultaneous fit for the radio sources and extended SZ emission, is also compared to techniques previously applied to Ryle Telescope observations and is found to be robust. The technique is then applied to new observations of the cluster Abell 611, and a decrement of -540 +/- 125 microJy/beam is found. This is combined with a ROSAT HRI image and a published ASCA temperature to give an Hubble constant estimate of 52+24-16 km/s/Mpc.

We present the preliminary results of a new survey of radio sources using the Ryle telescope at 15.2 GHz. This is the highest frequency at which a survey has been done that is relevant to the issue of radio source contamination in CMB experiments. The differential source count of the 66 sources found in 63 sqdeg is 80(S/Jy)^-2 /Jy/sr from about 20 to 500 mJy. Extrapolating this to 34 GHz (where many cm-wave CMB experiments operate) gives an estimated temperature contribution from sources of 9 microK in a CMB image, with a beam corresponding to multipole l=500. A means of source subtraction is evidently necessary, otherwise the signal-to-noise ratio in CMB images will be limited to 4 or 5, becoming worse at higher resolution. We compare the population of sources observed in this new survey to that predicted by extrapolation from lower frequency surveys, finding that source fluxes, and indeed the existence of many sources, cannot be determined by extrapolation.

We demonstrate that the maximum-entropy method for gravitational lens reconstruction presented in Bridle et al. (1998) may be applied even when only shear \emph{or} magnification information is present. We also demonstrate that the method can easily handle irregularly shaped observing fields and, because shear is a non-local function of the lensing mass, reconstructions that use shear information can successfully bridge small gaps in observations. For our simulations we use a mass density distribution that is realistic for a z=0.4 cluster of total mass around 10^15 h^-1 M_solar. Using HST-quality shear data alone, covering the area of four WFPC2 observations, we detect 60 per cent of the mass of the cluster within the area observed, despite the mass sheet degeneracy. This is qualitatively because the shear provides information about the variations in the mass distribution, and our prior includes a positivity constraint. We investigate the effect of using various sizes of observing field and find that 50 to 100 per cent of the cluster mass is detected, depending on the observing strategy and cluster shape. Finally we demonstrate how this method can cope with strong lensing regions of a mass distribution.

We have carried out a study to simulate distant clusters of galaxies in deep ground-based optical images. We find that when model galaxies are added to deep images obtained with the William Herschel Telescope, there is considerable scatter of the recovered galaxy colours away from the model values; this scatter is larger than that expected from photometric errors and is significantly affected by confusion, due to ground-based seeing, between objects in the field. In typical conditions of $\approx$ 1-arcsec seeing, the combination of confusion and incompleteness causes a considerable underestimation of the true surface density of $z \approx 3$ galaxies. We argue that the actual surface density of $z \approx 3$ galaxies may be several times greater than that estimated by previous ground-based studies, consistent with the surface density of such objects found in the HDF.

We describe observations at frequencies near 15 GHz of the second 2x2 degree field imaged with the Cambridge Cosmic Anisotropy Telescope (CAT). After the removal of discrete radio sources, structure is detected in the images on characteristic scales of about half a degree, corresponding to spherical harmonic multipoles in the approximate range l= 330--680. A Bayesian analysis confirms that the signal arises predominantly from the cosmic microwave background (CMB) radiation for multipoles in the lower half of this range; the average broad-band power in a bin with centroid l=422 (theta = 51') is estimated to be Delta_T/T=2.1^{+0.4}_{-0.5} x 10^{-5}. For multipoles centred on l=615 (theta =35'), we find contamination from Galactic emission is significant, and constrain the CMB contribution to the measured power in this bin to be Delta_T/T <2.0 x 10^{-5} (1-sigma upper limit). These new results are consistent with the first detection made by CAT in a completely different area of sky. Together with data from other experiments, this new CAT detection adds weight to earlier evidence from CAT for a downturn in the CMB power spectrum on scales smaller than 1 degree. Improved limits on the values of H_0 and Omega are determined using the new CAT data.

We present deep optical images of the PC1643+4631 field obtained at the WHT. This field contains two quasars at redshifts z=3.79 & 3.83 and a cosmic microwave background (CMB) decrement detected with the Ryle Telescope. The images are in U,G,V,R and I filters, and are complete to 25th magnitude in R and G and to 25.5 in U. The isophotal galaxy counts are consistent with the results of Metcalde et al. (1996), Hogg et al. (1997), and others. We find an excess of robust high-redshift Ly-break galaxy candidates with 25.0<R<25.5 compared with the mean number found in the fields studied by Steidel et al. -we expect 7 but find 16 - but we do not find that the galaxies are concentrated in the direction of the CMB decrement. However, we are still not sure of the distance to the system causing the CMB decrement. We have also used our images to compare the commonly used object-finding algorithms of FOCAS and SExtractor: we find FOCAS the more efficient at detecting faint objects and the better at dealing with composite objects, whereas SExtractor's morphological classification is more reliable, especially for faint objects near the resolution limit. More generally, we have also compared the flux lost using isophotal apertures on a real image with that on a noise-only image: recovery of artificial galaxies from the noise-only image significantly overestimates the flux lost from the galaxies, and we find that the corrections made using this technique suffer a systematic error of some 0.4 magnitudes.

In an investigation of the cause of the cosmic microwave background decrement in the field of the z = 3.8 quasar pair PC1643+4631, we have carried out a study to photometrically estimate the redshifts of galaxies in deep multi-colour optical images of the field taken with the WHT. To examine the possibility that a massive cluster of galaxies lies in the field, we have attempted to recover simulated galaxies with intrinsic colours matching those of the model galaxies used in the photometric redshift estimation. We find that when such model galaxies are added to our images, there is considerable scatter of the recovered galaxy redshifts away from the model value; this scatter is larger than that expected from photometric errors and is the result of confusion, simply due to ground-based seeing, between objects in the field. We have also compared the likely efficiency of the photometric redshift technique against the colour criteria used to select z>3 galaxies via the strong colour signature of the Lyman-limit break. We find that these techniques may significantly underestimate the true surface density of z>3, due to confusion between the high-redshift galaxies and other objects near the line of sight. We argue that the actual surface density of z=3 galaxies may be as much as 6 times greater than that estimated by previous ground-based studies, and note that this conclusion is consistent with the surface density of high-redshift objects found in the HDF. Finally, we conclude that all ground-based deep field surveys are inevitably affected by confusion, and note that reducing the effective seeing in ground-based images will be of paramount importance in observing the distant universe.

The maximum-entropy method is applied to the problem of reconstructing the projected mass density of a galaxy cluster using its gravitational lensing effects on background galaxies. We demonstrate the method by reconstructing the mass distribution in a model cluster using simulated shear and magnification data to which Gaussian noise is added. The mass distribution is reconstructed directly and the inversion is regularised using the entropic prior for this positive additive distribution. For realistic noise levels, we find that the method faithfully reproduces the main features of the cluster mass distribution not only within the observed field but also slightly beyond it. We estimate the uncertainties on the reconstruction by calculating an analytic approximation to the covariance matrix of the reconstruction values of each pixel. This result is compared with error estimates derived from Monte-Carlo simulations for different noise realisations and found to be in good agreement.

In a 15-GHz Ryle-Telescope observation of PC1643+4631 A & B, a pair of quasars at redshifts z = 3.79 and 3.83 separated by 198 arcsec on the sky, we find a decrement in the cosmic microwave background (CMB) of -380 +/- 64 uJy in a 110 arcsec x 175 arcsec beam. Assuming this to be a Sunyaev-Zel'dovich effect due to an intervening cluster, the minimum magnitude of the central temperature decrement is 560 uK. A serendipitous ROSAT observation shows that there is no X-ray-luminous cluster in the direction of the decrement at z < 1. The implied gas mass is >~ 2 x 10^{14} solar masses (assuming a temperature of ~ 5 keV), indicating a total mass of > 10^{15} solar masses. This result demonstrates the existence of a massive system too distant to be detected by its emission, but which can be found via its imprint on the CMB.

In a companion paper (Paper I: Jones et al. 1996) we report the discovery of a cosmic microwave background decrement, indicative of a distant cluster with mass ~10^{15} solar masses, towards the quasar pair PC1643+4631A&B (z=3.79,3.83, separation 198''). To search for the cluster responsible, we have obtained R-, J- and K-band images of the field and have also carried out optical spectroscopy of selected objects in it. No such cluster is evident in these images. Assuming the cluster causing the decrement is similar to massive clusters already known, our magnitude limits imply that it must lie about or beyond z=1. This provides independent support for the X-ray-based distance argument of Paper I. The cluster must gravitationally lens objects behind it; for a cluster z around 1-2, the Einstein-ring radius for sources at z ~= 3.8 is ~ 100''. Simple modelling, producing simultaneously the S-Z effect and the lensing, shows that the source positions of quasars A and B lie within ~10'' of each other and may indeed be coincident. The two quasar spectra are found to be remarkably similar apart from their one-percent redshift difference. Assuming A and B are images of a single quasar, we present a possible explanation of this difference.

The Ryle Telescope has been used to provide images at 15 GHz of the Sunyaev-Zel'dovich decrements towards a dozen clusters in an X-ray luminosity-limited sample. So far, X-ray data have allowed H_0 estimates towards two of these, both giving ``low'' values that are self-consistent, though it is essential to obtain many more estimates to reduce the effects of cluster projection. We have also discovered a decrement towards the z=3.8 quasar pair PC1643+4631A&B (198'' separation); substantial X-ray, optical and infrared follow-up show that the cluster responsible is either at z>1 or very underluminous and that very significant gravitational lensing must be occuring.