The Baryon acoustic oscillations from Integrated Neutral Gas Observations (BINGO) telescope is a 40-m~class radio telescope under construction that has been designed to measure the large-angular-scale intensity of HI emission at 980--1260 MHz and hence to constrain dark energy parameters. A large focal plane array comprising of 1.7-metre diameter, 4.3-metre length corrugated feed horns is required in order to optimally illuminate the telescope. Additionally, very clean beams with low sidelobes across a broad frequency range are required, in order to facilitate the separation of the faint HI emission from bright Galactic foreground emission. Using novel construction methods, a full-sized prototype horn has been assembled. It has an average insertion loss of around 0.15 dB across the band, with a return loss around -25 dB. The main beam is Gaussian with the first sidelobe at around $-25 dB. A septum polariser to separate the signal into the two hands of circular polarization has also been designed, built and tested.

HI Intensity Mapping (IM) will be used to do precision cosmology using many existing and upcoming radio observatories. The signal will be contaminated due to absorption, the largest component of which will be the flux absorbed by the HI emitting sources themselves from the flux incident on them from background radio continuum sources. We, for the first time, provide a quantitative estimate of the magnitude of the absorbed flux compared to the emitted HI flux for various voxels placed at redshifts between 0.1 and 2.5. We use a cosmological sky simulation of the atomic HI emission line, and sum over the emitted and absorbed fluxes for all sources within voxels at different redshifts. For estimating the absorbed flux we use various relations based on existing observations as well as simulations. We find that for the same co-moving volume of sky, the HI emission falls off quickly with increasing redshift, while the absorption varies much less with redshift and follows the redshift distribution of faint sources that dominate the number counts of radio continuum sources. This results in the fraction of absorption compared to emission to be negligible in the nearby Universe (up to a redshift of ~0.5), increases to about 10% at a redshift of 1, and continues to increase to about 30% up to a redshift of 2.5. These numbers can vary significantly due to the uncertainties on the exact forms of the various relations used, the largest variation being driven by the uncertainty on the number counts of radio continuum sources at sub-mJy flux densities. Absorption of flux incident from background radio continuum sources might become an important contaminant to HI IM signals beyond redshifts of 0.5, and needs to be quantified more accurately using inputs from upcoming deep high resolution surveys of radio continuum sources, HI absorption, and HI emission with the SKA and its precursors.

The Baryon acoustic oscillations from Integrated Neutral Gas Observations (BINGO) telescope is a new 40-m class radio telescope to measure the large-angular-scale intensity of Hi emission at 980-1260 MHz to constrain dark energy parameters. As it needs to measure faint cosmological signals at the milliKelvin level, it requires a site that has very low radio frequency interference (RFI) at frequencies around 1 GHz. We report on measurement campaigns across Uruguay and Brazil to find a suitable site, which looked at the strength of the mobile phone signals and other radio transmissions, the location of wind turbines, and also included mapping airplane flight paths. The site chosen for the BINGO telescope is a valley at Serra do Urubu, a remote part of Paraiba in North-East Brazil, which has sheltering terrain. During our measurements with a portable receiver we did not detect any RFI in or near the BINGO band, given the sensitivity of the equipment. A radio quiet zone around the selected site has been requested to the Brazilian authorities ahead of the telescope construction.

The existence of spiral DRAGNs challenges standard galaxy formation theories. We present the first observation of 0313$-$192, the archetypal spiral DRAGN, at VLBI resolutions. Spiral DRAGNs are Double Radio Sources Associated with Galactic Nuclei (DRAGNs) that are hosted by spiral galaxies. 0313$-$192 is an edge-on spiral galaxy that appears to host a 360 kpc double-lobed radio source. The core of this galaxy is clearly detected at L, S, and X-bands using the VLBA, signifying an ongoing active nucleus in the galaxy. This rules out the possibility that the spiral DRAGN is merely a chance alignment. The radio core has L$_{1.4\,GHz} \sim 3.0 \times 10^{23}$W Hz$^{-1}$. Radio components are detected to the South-West of the core, but there are no detections of a counterjet. Assuming a symmetric, relativistic jet, we estimate an upper limit to the inclination angle of $\theta \lesssim 72$ degrees. The VLBI-detected radio jet components are extremely well-aligned with the larger-scale radio source suggesting little to no jet disruption or interaction with the ISM of the host galaxy.

We have reanalysed the 1996/1997 VLA monitoring data of the gravitational lens system JVAS B0218+357 to produce improved total flux density and polarization variability curves at 15, 8.4 and 5 GHz. This has been done using improved calibration techniques, accurate subtraction of the emission from the Einstein ring and careful correction of various systematic effects, especially an offset in polarization position angle that is hour-angle dependent. The variations in total and polarized flux density give the best constraints and we determine a combined delay estimate of $11.3 \pm 0.2$ d (1$\sigma$). This is consistent with the $\gamma$-ray value recently derived using the Fermi Gamma-ray Space Telescope and thus we find no evidence for a positional shift between the radio and $\gamma$-ray emitting regions. Combined with the previously published lens model found using LensClean, the new delay gives a value for the Hubble constant of $H_0 = 72.9 \pm 2.6$ km s$^{-1}$ Mpc$^{-1}$ (1$\sigma$).

Prompted by a recent claim by Barrows et al. that X-ray AGN are often found significantly offset from the centres of their host galaxies, we have looked for examples of compact radio sources which are offset from the optical centroids of nearby (z < 0.2) galaxies. We have selected a sample of 345 galaxies from the Sloan Digital Sky Survey (SDSS) galaxy catalog which have nearby compact radio sources listed in the Cosmic-Lens All Sky Survey (CLASS) catalog. We find only three matches (0.87 per cent of the sample) with offsets greater than 600 milliarcsec (mas), which is considerably fewer than we would have expected from the Barrows et al. X-ray survey. We fit our histogram of offsets with a Rayleigh distribution with {\sigma} = 60.5 mas, but find that there is an excess of objects with separations greater than approximately 150 mas. Assuming that this excess represents AGN with real offsets, we place an upper limit of approximately 17 per cent on the fraction of offset AGN in our radio-selected sample. We select 38 objects with offsets greater than 150 mas, and find they have some diverse properties: some are well known, such as Mrk 273 and Arp 220, and others have dust lanes which may have affected the optical astrometry, while a few are strong new candidates for offset AGN.

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

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

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

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

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

Gaia will open up a huge volume of new parameter space in which to explore the physics of AGN and black hole evolution. We address the question as to how far along the relativistic jets blazar radio, optical and gamma ray emission originated. In some models the optical centroid wander should be detectable as the relative contributions of thermal and non-thermal optical emission change. Black holes powering AGN do not necessarily reside at the centres of their host galaxies; they can be one member of a binary pair or they could have received a kick after binary coalescence. For radio-loud AGN comparison of astrometric radio and and optical positions can reveal such displacements. It is suggested that it would be feasible to do this using Gaia and e-MERLIN for a sample of thousands of elliptical radio galaxies.

Unified schemes of radio sources, which account for different types of radio AGN in terms of anisotropic radio and optical emission, together with different orientations of the ejection axis to the line of sight, have been invoked for many years. Recently, large samples of optical quasars, mainly from the Sloan Digital Sky Survey, together with large radio samples, such as FIRST, have become available. These hold the promise of providing more stringent tests of unified schemes but, compared to previous samples, lack high resolution radio maps. Nevertheless they have been used to investigate unified schemes, in some cases yielding results which appear inconsistent with such theories. Here we investigate using simulations how the selection effects to which such investigations are subject can influence the conclusions drawn. In particular, we find that the effects of limited resolution do not allow core-dominated radio sources to be fully represented in the samples, that the effects of limited sensitivity systematically exclude some classes of sources and the lack of deep radio data make it difficult to decide to what extent closely separated radio sources are associated. Nevertheless, we conclude that relativistic unified schemes are entirely compatible with the current observational data. For a sample selected from SDSS and FIRST which includes weak-cored triples we find that the equivalent width of the [OIII] emission line decreases as core-dominance increases, as expected, and also that core-dominated quasars are optically brighter than weak-cored quasars.

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

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

The Fermi-LAT survey provides a large sample of blazars selected on the strength of their inverse Compton emission. We cross-correlate the first Fermi-LAT catalogue with the CRATES radio catalogue and use this sample to investigate whether blazar gamma-ray luminosities are influenced by the availability of external photons to be up-scattered. Using the 8.4 GHz flux densities of their compact radio cores as a proxy for their jet power, we calculate their Compton Efficiency parameters, which measure the ability of jets to convert power in the form of ultra-relativistic electrons into Compton gamma-rays. We find no clear differences in Compton efficiencies between BL Lac objects and FSRQs and no anti-correlation between Compton efficiency and synchrotron peak frequency. This suggests that the scattering of external photons is energetically unimportant compared to the synchrotron self-Compton process. These results contradict the predictions of the blazar sequence.

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

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

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

The Hubble constant value is currently known to 10% accuracy unless assumptions are made for the cosmology (Sandage et al. 2006). Gravitational lens systems provide another probe of the Hubble constant using time delay measurements. However, current investigations of ~20 time delay lenses, albeit of varying levels of sophistication, have resulted in different values of the Hubble constant ranging from 50-80 km/s/Mpc. In order to reduce uncertainties, more time delay measurements are essential together with better determined mass models (Oguri 2007, Saha et al. 2006). We propose a more efficient technique for measuring time delays which does not require regular monitoring with a high-resolution interferometer array. The method uses double image and long-axis quadruple lens systems in which the brighter component varies first and dominates the total flux density. Monitoring the total flux density with low-resolution but high sensitivity radio telescopes provides the variation of the brighter image and is used to trigger high-resolution observations which can then be used to see the variation in the fainter image. We present simulations of this method together with a pilot project using the WSRT (Westerbork Radio Synthesis Telescope) to trigger VLA (Very Large Array) observations. This new method is promising for measuring time delays because it uses relatively small amounts of time on high-resolution telescopes. This will be important because many SKA pathfinder telescopes, such as MeerKAT (Karoo Array Telescope) and ASKAP (Australian Square Kilometre Array Pathfinder), have high sensitivity but limited resolution.

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

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

We present the first results from the new Survey of Extragalactic Nuclear Spectral Energies (SENSE) sample of "blazars". The sample has been chosen with minimal selection effects and is therefore ideal to probe the intrinsic properties of the blazar population. We report evidence for negative cosmological evolution in this radio selected sample and give an outline of future work related to the SENSE sample.

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

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

Aims. We investigate the properties of the host galaxy of the blazar J0324+3410 (B2 0321+33) by the analysis of B and R images obtained with the NOT under good photometric conditions. Methods: The galaxy was studied using different methods: Sersic model fitting, unsharp-masked images, B-R image and B-R profile analysis. Results: The images show that the host galaxy has a ring-like morphology. The B-R colour image reveals two bluish zones: one that coincides with the nuclear region, interpreted as the signature of emission related to the active nucleus, the other zone is extended and is located in the host ring-structure. We discuss the hypothesis that the later is thermal emission from a burst of star formation triggered by an interacting/merging process.

We report the results of a survey of 442 planetary nebulae at 30 GHz. The purpose of the survey is to develop a list of planetary nebulae as calibration sources which could be used for high frequency calibration in future. For 41 PNe with sufficient data, we test the emission mechanisms in order to evaluate whether or not spinning dust plays an important role in their spectra at 30 GHz. The 30-GHz data were obtained with a twin-beam differencing radiometer, OCRA-p, which is in operation on the Torun 32-m telescope. Sources were scanned both in right ascension and declination. We estimated flux densities at 30 GHz using a free-free emission model and compared it with our data. The primary result is a catalogue containing the flux densities of 93 planetary nebulae at 30 GHz. Sources with sufficient data were compared with a spectral model of free-free emission. The model shows that free-free emission can generally explain the observed flux densities at 30 GHz thus no other emission mechanism is needed to account for the high frequency spectra.

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

Context: We present the soft X-ray properties obtained in the ROSAT All-Sky survey and from pointed PSPC observations for the AGN in the complete flux-density limited Caltech-Jodrell Bank flat spectrum sample (hereafter CJF). CJF is a VLBI survey (VLBA observations at 5 GHz) of 293 AGN with detailed information on jet component motion. Aims: We investigate and discuss the soft X-ray properties of this AGN sample and examine the correlations between X-ray and VLBI properties, test beaming scenarios, and search for the discriminating properties between the sub-samples detected and not detected by ROSAT. Methods: Comparing the observed and the predicted X-ray fluxes by assuming an Inverse Compton (IC) origin for the observed X-rays, we compute the beaming or Doppler factor and contrast the Doppler factors with other beaming indicators derived from the VLBI observations, such as the value of the expansion velocity, and the observed and intrinsic brightness temperature. In addition, we investigate the large-scale radio structure of the AGN and the difference between the pc- and kpc-scale structure (misalignment) with regard to the X-ray observations.Results: We find a nearly linear relation between X-ray and radio luminosities, and a similar but less stringent behaviour for the relation between optical and X-ray luminosities. The quasars detected by ROSAT have a different $\beta_{\rm app}$-redshift relationship compared to the non-detected ones. ROSAT-detected sources tend to reveal extended large-scale radio structures more often. Conclusions: We conclude that beaming alone cannot explain the observed dichotomy of ROSAT detection or non-detection and assume that the large-scale jet structure plays a decisive role.

Context: This is the second of a series of papers presenting VLBI observations of the 293 Caltech-Jodrell Bank Flat-Spectrum (hereafter CJF) sources and their analysis. Aims: To obtain a consistent motion dataset large enough to allow the systematic properties of the population to be studied. Methods: We present the detailed kinematic analysis of the complete flux-density limited CJF survey. We computed 2-D kinematic models based on the optimal model-fitting parameters of multi-epoch VLBA observations. This allows us to calculate not only radial, but also orthogonal motions, and thus to study curvature and acceleration. Statistical tests of the motions measured and their reliability have been performed. A correlation analysis between the derived apparent motions, luminosities, spectral indices, and core dominance and the resulting consequences is described. Results: With at least one velocity in each of 237 sources, this sample is much larger than any available before and allows a meaningful statistical investigation of apparent motions and any possible correlations with other parameters in AGN jets. Conclusions: This AGN survey provides the basis for any statistical analysis of jet and jet-component properties.

Context: This is the first of a series of papers presenting VLBI observations of the 293 Caltech-Jodrell Bank Flat-Spectrum (hereafter CJF) sources and their analysis. Aims: One of the major goals of the CJF is to make a statistical study of the apparent velocities of the sources. Methods: We have conducted global VLBI and VLBA observations at 5 GHz since 1990, accumulating thirteen separate observing campaigns. Results: We present here an overview of the observations, give details of the data reduction and present the source parameters resulting from a model-fitting procedure. For every source at every observing epoch, an image is shown, built up by restoring the model-fitted components, convolved with the clean beam, into the residual image, which was made by Fourier transforming the visibility data after first subtracting the model-fitted components in the uv-plane. Overplotted we show symbols to represent the model components. Conclusions: We have produced VLBI images of all but 5 of the 293 sources in the complete CJF sample at several epochs and investigated the kinematics of 266 AGN.

To measure the 30-GHz flux densities of the 293 sources in the Caltech-Jodrell Bank flat-spectrum (CJF) sample. The measurements are part of an ongoing programme to measure the spectral energy distributions of flat spectrum radio sources and to correlate them with the milliarcsecond structures from VLBI and other measured astrophysical properties. The 30-GHz data were obtained with a twin-beam differencing radiometer system mounted on the Torun 32-m telescope. The system has an angular resolution of 1.2 arcmin. Together with radio spectral data obtained from the literature, the 30-GHz data have enabled us to identify 42 of the CJF sources as Giga-hertz Peaked Spectrum (GPS) sources. Seventeen percent of the sources have rising spectra (alpha > 0) between 5 and 30 GHz.

We report a probable gravitational lens J0316+4328, one of 19 candidate asymmetric double lenses (2 images at a high flux density ratio) from CLASS. Observations with the Very Large Array (VLA), MERLIN and the Very Long Baseline Array (VLBA) imply that J0316+4328 is a lens with high confidence. It has 2 images separated by 0.40", with 6 GHz flux densities of 62 mJy and 3.2 mJy. The flux density ratio of ~19 (constant over the frequency range 6-22 GHz) is the largest for any 2 image gravitational lens. High resolution optical imaging and deeper VLBI maps should confirm the lensing interpretation and provide inputs to detailed lens models. The unique configuration will give strong constraints on the lens galaxy's mass profile.

We present 30 GHz Sunyaev Zel'dovich (SZ) observations of a sample of four galaxy clusters with a prototype of the One Centimetre Receiver Array (OCRA-p) which is mounted on the Torun 32-m telescope. The clusters (Cl0016+16, MS0451.6-0305, MS1054.4-0321 and Abell 2218) are popular SZ targets and serve as commissioning observations. All four are detected with clear significance (4-6 sigma) and values for the central temperature decrements are in good agreement with measurements reported in the literature. We believe that systematic effects are successfully suppressed by our observing strategy. The relatively short integration times required to obtain these results demonstrate the power of OCRA-p and its successors for future SZ studies.

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

We have used the very large JVAS/CLASS 8.4-GHz surveys of flat-spectrum radio sources to obtain a large, uniformly observed and calibrated, sample of radio source polarizations. These are useful for many investigations of the properties of radio sources and the interstellar medium. We discuss comparisons with polarization measurements from this survey and from other large-scale surveys of polarization in flat-spectrum sources.

Central gravitational image detection is very important for the study of the mass distribution of the inner parts ($\sim 100$ pc) of lens galaxies. However, the detection of such images is extremely rare and difficult. We present a 1.7-GHz High Sensitivity Array (HSA) observation of the double-image radio lens system B1030+074. The data are combined with archive VLBA and global-VLBI observations, and careful consideration is given to the effects of noise, {\sc clean}ing and self-calibration. An upper limit is derived for the strength of the central image of 180 $\mu$Jy (90% confidence level), considerably greater than would have been expected on the basis of a simple analysis. This gives a lower limit of $\sim 10^3$ for the ratio of the brightest image to the central image. For cusped models of lens mass distributions, we have made use of this non-detection to constrain the relation between inner power-law slope $\beta$ of the lensing galaxy mass profile, and its break radius $r_b$. For $r_b>130$ pc the power-law slope is required to be close to isothermal ($\beta>1.8$). A flatter inner slope is allowed if a massive black hole is present at the centre of the lensing galaxy, but the effect of the black hole is small unless it is $\sim 10$ times more massive than that implied by the relation between black hole mass and stellar velocity dispersion. By comparing four epochs of VLBI observations, we also detected possible superluminal motion in the jet in the brighter A image. The B jet remains unresolved, as expected from a simple lens model of the system.

We present the selection of the Jodrell Bank Flat-spectrum (JBF) radio source sample, which is designed to reduce the uncertainties in the Cosmic Lens All-Sky Survey (CLASS) gravitational lensing statistics arising from the lack of knowledge about the parent population luminosity function. From observations at 4.86 GHz with the Very Large Array, we have selected a sample of 117 flat-spectrum radio sources with flux densities greater than 5 mJy. These sources were selected in a similar manner to the CLASS complete sample and are therefore representative of the parent population at low flux densities. The vast majority (~90 per cent) of the JBF sample are found to be compact on the arcsecond scales probed here and show little evidence of any extended radio jet emission. Using the JBF and CLASS complete samples we find the differential number counts slope of the parent population above and below the CLASS 30 mJy flux density limit to be -2.07+/-0.02 and -1.96+/-0.12, respectively.

Many lens surveys have hitherto used observations of large samples of background sources to select the small minority which are multiply imaged by lensing galaxies along the line of sight. Recently surveys such as SLACS and OLS have improved the efficiency of surveys by pre-selecting double-redshift systems from SDSS. We explore other ways to improve survey efficiency by optimum use of astrometric and morphological information in existing large-scale optical and radio surveys. The method exploits the small position differences between FIRST radio positions of lensed images and the SDSS lens galaxy positions, together with the marginal resolution of some larger gravitational lens systems by the FIRST beam. We present results of a small pilot study with the VLA and MERLIN, and discuss the desirable criteria for future surveys.

We present the results of phase-referenced VLBA+Effelsberg observations at five frequencies of the double-image gravitational lens JVAS B0218+357, made to establish the precise registration of the A and B lensed image positions. The motivation behind these observations is to investigate the anomalous variation of the image flux density ratio (A/B) with frequency - this ratio changes by almost a factor of two over a frequency range from 1.65 GHz to 15.35 GHz. We investigate whether frequency dependent image positions, combined with a magnification gradient across the image field, could give rise to the anomaly. Our observations confirm the variation of image flux ratio with frequency. The results from our phase-reference astrometry, taken together with the lens mass model of Wucknitz et al. (2004), show that shifts of the image peaks and centroids are too small to account for the observed frequency-dependent ratio.

We present deep optical observations of the gravitational lens system CLASS B0218+357, from which we derive an estimate of the Hubble constant. Extensive radio observations have reduced the degeneracies between Hubble's constant and the mass model in this lens to one involving only the position of the radio-quiet lensing galaxy relative to the lensed images. B0218+357 has an image separation of only 334 mas, so optical observations have previously been unable to resolve the lens galaxy from the bright lensed images. Using the new Advanced Camera for Surveys installed on the Hubble Space Telescope, we have measured the separation between the lens galaxy centre and the brightest image. The position found, and hence our estimate of Hubble's constant, depends on our approach to the spiral arms in B0218+357. If the most prominent arms are left unmasked, we find a value for Hubble's constant of 70+/-5 km/s /Mpc (95% confidence). If the spiral arms are masked out, we find a value of 61+/-7 km/s /Mpc (95% confidence).

We address the question of whether or not the properties of all low-luminosity flat spectrum radio sources, not just the obvious BL Lac objects, are consistent with them being the relativistically beamed counterparts of the low radio luminosity radio galaxies (the Fanaroff-Riley type 1 - FRI). We compare the emission-line luminosities of a well studied sample of core-dominated radio sources with a matched sample of FR I radio galaxies. The emission lines in the core-dominated objects are on average significantly more luminous than those in the comparison sample, inconsistent with the simplest unified models in which there is no orientation dependence of the line emission. However, the results are equally consistent with a scenario in which, for the majority of objects, beaming has minimal effect on the observed core luminosities of a large fraction of the FRI population and that intrinsically stronger cores simply give rise to stronger emission lines. We conclude that FRI unification is much more complex than usually portrayed, and models combining beaming with an intrinsic relationship between core and emission line strengths need to be explored.

We report the discovery of a new gravitational lens system from the CLASS survey, CLASS B0631+519. VLA, MERLIN and VLBA observations show a doubly-imaged radio core, a doubly-imaged lobe and a second lobe that is probably quadruply-imaged. The maximum image separation is 1.16 arcseconds. The VLBA resolves the most magnified image of the flat-spectrum radio core into a number of sub-components spread across approximately 20 milli-arcseconds. Optical and near-infrared imaging with the ACS and NICMOS cameras on the HST show that there are two galaxies along the line of sight to the lensed source, as was previously discovered by optical spectroscopy. The nearer galaxy at z=0.0896 is a small blue irregular, while the more distant galaxy at z=0.6196 is an elliptical type and appears to contribute most of the lensing effect. The host galaxy of the lensed source is visible in NICMOS imaging as a set of arcs that form an almost complete Einstein ring. Mass modelling using non-parametric techniques can reproduce the ring and indicates that the irregular galaxy has a (localised) effect on the flux density distribution in the Einstein ring at the 5-10% level.

We present observations of CLASS B2108+213, the widest separation gravitational lens system discovered by the Cosmic Lens All-Sky Survey. Radio imaging using the VLA at 8.46 GHz and MERLIN at 5 GHz shows two compact components separated by 4.56 arcsec with a faint third component in between which we believe is emission from a lensing galaxy. 5-GHz VLBA observations reveal milliarcsecond-scale structure in the two lensed images that is consistent with gravitational lensing. Optical emission from the two lensed images and two lensing galaxies within the Einstein radius is detected in Hubble Space Telescope imaging. Furthermore, an optical gravitational arc, associated with the strongest lensed component, has been detected. Surrounding the system are a number of faint galaxies which may help explain the wide image separation. A plausible mass distribution model for CLASS B2108+213 is also presented.

We analyse a group of radio sources, a subset of the 200 mJy sample, all of which have core-jet radio structures measured with VLBI and have flat spectra stretching from the radio to the millimetre/sub-millimetre band. Thus the objects have most of the properties expected of blazars. However, they display varied optical properties ranging from ``Seyfert-like'' objects, through BL Lac objects, to ``normal'' elliptical galaxies. We investigate the distribution of synchrotron peak frequencies in their Spectral Energy Distributions (SEDs) and find a broad distribution between 10^12 and 10^16 Hz. Our conclusion is that we should consider virtually all objects in the sample as blazars since much of the diversity in their classification based on traditional optical criteria arises from differences in the frequency at which the non-thermal emission begins to decline. Specifically, an object is only classified as BL Lac when its peak frequency falls in the near IR/optical range. We determine peak frequencies using the same method for objects from other blazar samples. An important result is that our objects do not follow the blazar spectral sequence proposed by Fossati et al. and Ghisellini et al. in which, on average, peak frequencies increase as the radio luminosity decreases. Most of our low radio-luminosity sources have peaks in their SEDs at low frequencies, not at the expected high frequencies. We suggest that at least part of the systematic trend seen by Fossati et al. and Ghisellini et al. results from selection effects.

The advent of new observational facilities in the last two decades has allowed the rapid discovery and high-resolution optical imaging of many strong lens systems from galaxy to cluster scales, as well as their spectroscopic follow-up. Radio telescopes have played the dominant role in the systematic detection of dozens of new arcsec-scale lens systems. For the future, we expect nothing less! The next major ground- and space-based facilities, especially the Square Kilometer Array can discover tens of thousands of new lens systems in large sky surveys. For optical imaging and spectroscopic follow-up a strong synergy with planned optical facilities is needed. Here, we discuss the field where strong gravitational lensing is expected to play the dominant role and where SKA can have a major impact: The study of the internal mass structure and evolution of galaxies and clusters to z~1. In addition, studies of more exotic phenomena are contemplated. For example, milli- and microlensing can provide a way to measure the mass-functions of stars and CDM substructure at cosmological distances. All-sky radio monitoring will also rapidly develop the field of time-domain lensing.

We address the question of why low-luminosity radio sources with similar flat radio spectra show a range of optical activity. The investigation is based on the spectral energy distributions (SEDs) of objects from the 200 mJy sample. We gathered new data from the VLA at 43 GHz, from SCUBA in the JCMT at 2000, 1350 and 850 micron and from the ISOPHOT instrument on ISO at 170, 90, 60 and 25 micron. There is considerable diversity amongst the SEDs of the objects: there are objects with steep broadband spectra between centimetre and millimetre bands (14% of the sample), there are those with flat broadband spectra over most of the spectral range (48% of the sample) and there are those which show pronounced sub-mm/infrared excesses (27% of the sample). Some objects of the first group have two sided radio morphology indicating that their pc-scale emission is not dominated by beamed jet emission. Amongst the objects that have smooth broadband spectra from the radio to the infrared there are passive elliptical galaxies as well as the expected BL Lacs objects. The most pronounced sub-mm/infrared excesses are shown by the broad emission line objects.

We present new observations of the gravitational lens system CLASS B0128+437. HST observations detect a very faint, extended object in I-band with no emission from the lensed images visible; no detection at all is made in V-band. The lens system is detected with much higher signal to noise with UKIRT in K-band, but the resolution is not sufficient to allow the lensed images and the lens galaxy to be separated. A careful astrometric calibration, however, suggests that the peak of the infrared emission corresponds to the two merging images A and B and therefore that the lensed images dominate at infrared wavelengths. The new radio data consist of VLBI radio images at three frequencies, 2.3, 5 and 8.4GHz, made with the VLBA and the 100-m Effelsberg telescope. The lensed source consists of three well-defined sub-components embedded in a more extended jet. Due to the fact that the sub-components have different spectral indices it is possible to determine which part of each image corresponds to the same source sub-component. Our main finding is that one of the images, B, looks very different to the others, there being no obvious division into separate sub-components and the image being apparently both broader and smoother. This is a consequence we believe of scatter-broadening in the ISM of the lensing galaxy. The large number of multiply-imaged source sub-components also provide an abundance of modelling constraints and we have attempted to fit an SIE+external shear model to the data, as well as utilising the novel method of Evans & Witt. It proves difficult in both cases, however, to obtain a satisfactory fit which strongly suggests the presence of sub-structure in the mass distribution of the lensing galaxy, perhaps of the kind that is predicted by CDM theories of structure formation.

We present the optical spectra of four newly discovered gravitational lenses from the Cosmic Lens All-Sky Survey (CLASS). These observations were carried out using the Low Resolution Imaging Spectrograph on the W. M. Keck-I Telescope as part of a program to study galaxy-scale gravitational lenses. From our spectra we found the redshift of the background source in CLASS B0128+437 (z_s=3.1240+-0.0042) and the lensing galaxy redshifts in CLASS B0445+123 (z_l=0.5583+-0.0003) and CLASS B0850+054 (z_l=0.5883+-0.0006). Intriguingly, we also discovered that CLASS B0631+519 may have two lensing galaxies (z_l,1=0.0896+-0.0001, z_l,2=0.6196+-0.0004). We also found a single unidentified emission line from the lensing galaxy in CLASS B0128+437 and the lensed source in CLASS B0850+054. We find the lensing galaxies in CLASS B0445+123 and CLASS B0631+519 (l,2) to be early-type galaxies with Einstein Radii of 2.8-3.0 h^{-1} kpc. The deflector in CLASS B0850+054 is a late-type galaxy with an Einstein Radius of 1.6 h^{-1} kpc.

B0218+357 is one of the most promising systems to determine the Hubble constant from gravitational lenses. Consisting of two bright resolved images plus an Einstein ring, it provides better constraints for the mass model than other systems. The main problem left until now was the poorly determined position of the lensing galaxy. After presenting results from classical lens modelling, we apply our improved version of LensClean which utilizes the Einstein ring for lens modelling purposes. The primary result using isothermal models is a well defined lens position which allows the first reliable measurement of the Hubble constant from this system. The result of H0=(78+-6) km/s/Mpc (2 sigma) is high compared with other lenses but compatible with the HST key project and WMAP results. We furthermore discuss effects of different radial mass profiles. The power-law exponent of the potential is constrained by VLBI data to be beta=1.04+-0.02, very close to isothermal. The effect on H0 is expected to be very small. We also present a composite map (lensed and unlensed) which shows the rich structure of B0218+357 on scales from milli-arcseconds to arcseconds. Finally we use a comparison of observations at different frequencies to investigate the question of possible weakening of one of the images by propagation effects and/or source shifts with frequency. The data clearly favour the model of significant extinction without noticeable source position shifts. The technical details of our variant of the LensClean method are presented in the accompanying Paper I.

We present flux-ratio curves of the fold and cusp (i.e. close multiple) images of six JVAS/CLASS gravitational lens systems. The data were obtained over a period of 8.5 months in 2001 with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5-GHz with 50 mas resolution, as part of a MERLIN Key-Project. Even though the time delays between the fold and cusp images are small (<~1 day) compared to the time-scale of intrinsic source variability, all six lens systems show evidence that suggests the presence of extrinsic variability. In particular, the cusp images of B2045+265 -- regarded as the strongest case of the violation of the cusp relation (i.e. the sum of the magnifications of the three cusp images add to zero) -- show extrinsic variations in their flux-ratios up to ~40 percent peak-to-peak on time scales of several months. Its low Galactic latitude of b=-10 degree and a line-of-sight toward the Cygnus superbubble region suggest that Galactic scintillation is the most likely cause. The cusp images of B1422+231 at b=+69 degree do not show strong extrinsic variability. Galactic scintillation can therefore cause significant scatter in the cusp and fold relations of some radio lens systems (up to 10 percent rms), even though these relations remain violated when averaged over a <~1 year time baseline.

In the Caltech-Jodrell Survey of bright flat-spectrum radio sources, 3-5 epochs have now been observed for nearly all 293 sources; a uniquely large sample. The derivation of component motions has not yet been completed; it is complicated by, for example, variability and jet curvature. Nevertheless, some basic results are clear. The average apparent velocity in CJF quasars is near 3c, while for radio galaxies and BL Lacs it is around 1c. The distribution of velocities is broad, indicating a broad distribution of jet Lorentz factors, with many low values, and/or a decoupling of the Lorentz factors between the cores and the moving jet components, due to bending, speed changes, or due to pattern motions (shocks). A hint has emerged at this workshop that the average apparent velocities might be lower at lower radio frequencies; this will need careful verification.

We report the final results of the search for gravitationally lensed flat-spectrum radio sources found in the combination of CLASS (Cosmic Lens All Sky Survey) and JVAS (Jodrell-Bank VLA Astrometric Survey). VLA observations of 16,503 sources have been made, resulting in the largest sample of arcsec-scale lens systems available. Contained within the 16,503 sources is a complete sample of 11,685 sources having two-point spectral indices between 1.4 and 5 GHz flatter than -0.5 and 5 GHz flux densities $\geq$30 mJy. A subset of 8,958 sources form a well-defined statistical sample suitable for analysis of the lens statistics. We describe the systematic process by which 149 candidate lensed sources were picked from the statistical sample on the basis of possessing multiple compact components in the 0.2 arcsec-resolution VLA maps. Candidates were followed up with 0.05 arcsec resolution MERLIN and 0.003 arcsec VLBA observations at 5 GHz and rejected as lens systems if they failed well-defined surface brightness and/or morphological tests. Maps for all the candidates can be found on the World Wide Web at http://www.jb.man.ac.uk/research/gravlens/index.html We summarize the properties of each of the 22 gravitational lens systems in JVAS/CLASS. Twelve are double-image systems, nine are four-image systems and one is a six-image system. Thirteen constitute a statistically well-defined sample giving a point-source lensing rate of 1:690$\pm$190. The interpretation of the results in terms of the properties of the lensing galaxy population and cosmological parameters will be published elsewhere. (Abridged)

The Cosmic Lens All-Sky Survey (CLASS) is an international collaborative program which has obtained high-resolution radio images of over 10000 flat-spectrum radio sources in order to create the largest and best studied statistical sample of radio-loud gravitationally lensed systems. With this survey, combined with detailed studies of the lenses found therein, constraints can be placed on the expansion rate, matter density, and dark energy (e.g. cosmological constant, quintessence) content of the Universe that are complementary to and independent of those obtained through other methods. CLASS is aimed at identifying lenses where multiple images are formed from compact flat-spectrum radio sources, which should be easily identifiable in the radio maps. Because CLASS is radio-based, dust obscuration in lensing galaxies is not a factor, and the relative insensitivity of the instrument to environmental conditions leads to nearly uniform sensitivity and resolution over the entire survey. In four observing seasons from 1994-1999, CLASS has observed 13783 radio sources with the VLA at 8.4 GHz at 0.2 arcsecond resolution. When combined with the JVAS survey, the CLASS sample contains over 16,000 images. A complete sample of 11685 flat-spectrum sources was observed, selected from GB6 catalogue at 4.85 GHz and the NVSS at 1.4 GHz. So far, CLASS has found 16 new gravitational lens systems, and the JVAS/CLASS survey contains a total of 22 lenses. (Abridged)

We report the discovery of a new gravitational lens system from the CLASS survey. Radio observations with the VLA, the WSRT and MERLIN show that the radio source B0850+054 is comprised of two compact components with identical spectra, a separation of 0.7 arcsec and a flux density ratio of 6:1. VLBA observations at 5 GHz reveal structures that are consistent with the gravitational lens hypothesis. The brighter of the two images is resolved into a linear string of at least six sub-components whilst the weaker image is radially stretched towards the lens galaxy. UKIRT K-band imaging detects an 18.7 mag extended object, but the resolution of the observations is not sufficient to resolve the lensed images and the lens galaxy. Mass modelling has not been possible with the present data and the acquisition of high-resolution optical data is a priority for this system.

A new two-image gravitational lens system has been discovered as a result of the Cosmic Lens All-Sky Survey (CLASS). Radio observations with the VLA, MERLIN and the VLBA at increasingly higher resolutions all show two components with a flux density ratio of ~7:1 and a separation of 1.34". Both components are compact and have the same spectral index. Followup observations made with the VLA at 8.4 GHz show evidence of a feature to the south-east of the brighter component and a corresponding extension of the weaker component to the north-west. Optical observations with the WHT show ~1.7" extended emission aligned in approximately the same direction as the separation between the radio components with an R-band magnitude of 21.8 +/- 0.4.

We derive constraints on cosmological parameters and the properties of the lensing galaxies from gravitational lens statistics based on the final Cosmic Lens All Sky Survey (CLASS) data. For a flat universe with a classical cosmological constant, we find that the present matter fraction of the critical density is $\Omega_{\rm m}=0.31^{+0.27}_{-0.14}$ (68%) $^{+0.12}_{-0.10}$ (systematic). For a flat universe with a constant equation of state for dark energy $w = p_x({pressure})/\rho_x({energy density})$, we find $w < -0.55^{+0.18}_{-0.11}$ (68%).

In this paper we present new observations of the gravitational lens system JVAS B0218+357 made with a global VLBI network at a frequency of 8.4 GHz. Our maps have an rms noise of 30 microJy/beam and with these we have been able to image much of the extended structure of the radio jet in both the A and B images at high resolution (~1 mas). The main use of these maps will be to enable us to further constrain the lens model for the purposes of H0 determination. We are able to identify several sub-components common to both images with the expected parity reversal, including one which we identify as a counter-jet. We have not been successful in detecting either the core of the lensing galaxy or a third image. Using a model of the lensing galaxy we have back-projected both of the images to the source plane and find that they agree well. However, there are small, but significant, differences which we suggest may arise from multi-path scattering in the ISM of the lensing galaxy. We also find an exponent of the radial mass distribution of approximately 1.04, in agreement with lens modelling of published 15-GHz VLBI data. Polarisation maps of each image are presented which show that the distributions of polarisation across images A and B are different. We suggest that this results from Faraday rotation and associated depolarisation in the lensing galaxy.

Multi-frequency observations of Mrk 348 are presented consisting of 2 epochs of MERLIN data at 5 GHz, ISOPHOT data at 170, 90, 60 and 25 micron, NOT images at U, B, R and I bands and data at 1350 micron from the archive for SCUBA. The new optical images reveal a disrupted arm that ends towards the eastern companion of Mrk 348, consistent with the hypothesis that Mrk 348 and its companion form an interacting system. 5 GHz MERLIN imaging shows that only one of the radio components of Mrk 348 is polarized (%P=5). The broadband spectrum of Mrk 348 is flat between the radio and millimetre bands, suggesting that synchrotron emission extends to high frequencies. Mrk 348 has many of the characteristics of a radio-loud object. We discuss the orientation of the radio axis of Mrk 348 with respect to the line of sight. We conclude that the evidence is conflicting, and the geometry in Mrk 348 is not well-described by a simple edge-on or face-on model.

This paper describes the automated optical APM identification of radio sources from the Jodrell Bank - VLA Astrometric Survey (JVAS), as used for the search for distant radio-loud quasars. The sample has been used to investigate possible relations between optical and radio properties of flat spectrum radio sources. From the 915 sources in the sample, 756 have an optical APM identification at a red (e) and/or blue (o) plate,resulting in an identification fraction of 83% with a completeness and reliability of 98% and 99% respectively. About 20% are optically identified with extended APM objects on the red plates, e.g. galaxies. However the distinction between galaxies and quasars can not be done properly near the magnitude limit of the POSS-I plates. The identification fraction appears to decrease from >90% for sources with a 5 GHz flux density of >1 Jy, to <80% for sources at 0.2 Jy. The identification fraction, in particular that for unresolved quasars, is found to be lower for sources with steeper radio spectra. In agreement with previous studies, we find that the quasars at low radio flux density levels also tend to have fainter optical magnitudes, although there is a large spread. In addition, objects with a steep radio-to-optical spectral index are found to be mainly highly polarised quasars, supporting the idea that in these objects the polarised synchrotron component is more prominent. It is shown that the large spread in radio-to-optical spectral index is possibly caused by source to source variations in the Doppler boosting of the synchrotron component [Abridged].

We present a series of high resolution radio and optical observations of the CLASS gravitational lens system B1152+199 obtained with the Multi-Element Radio-Linked Interferometer Network (MERLIN), Very Long Baseline Array (VLBA) and Hubble Space Telescope (HST). Based on the milliarcsecond-scale substructure of the lensed radio components and precise optical astrometry for the lensing galaxy, we construct models for the system and place constraints on the galaxy mass profile. For a single galaxy model with surface mass density Sigma(r) propto r^-beta, we find that 0.95 < beta < 1.21 at 2-sigma confidence. Including a second deflector to represent a possible satellite galaxy of the primary lens leads to slightly steeper mass profiles.

This paper presents a new complete and well defined sample of Flat Spectrum Radio Sources (FSRS) selected from the CLASS survey with the further constraint of a bright (mag<17.5) optical counterpart. The sample has been designed to produce a large number of low-luminosity blazars in order to test the current unifying models in the low luminosity regime. In this first paper the new sample is presented and the radio properties of the 325 sources contained therein discussed.

We report the discovery of a new two-image gravitational lens system from the Cosmic Lens All-Sky Survey, CLASS B2319+051. Radio imaging with the Very Large Array (VLA) and Multi-Element Radio-Linked Interferometer Network (MERLIN) shows two compact components with a flux density ratio of 5:1, separated by 1.36 arcsec. Observations with the Very Long Baseline Array (VLBA) resolve each of the radio components into a pair of parity-reversed subcomponents. Hubble Space Telescope (HST) observations with the Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) show a bright elliptical galaxy (G1) coincident with the radio position, and a second irregular galaxy (G2) 3.4 arcsec to the northwest. Previous spectroscopic studies have indicated that these galaxies are at different redshifts: z(G1) = 0.624, z(G2) = 0.588. Infrared counterparts to the lensed radio components are not detected in the NICMOS image, and the source redshift has not yet been determined. Preliminary mass modeling based on the VLBA subcomponent data indicates that the lensing potential includes a strong external shear contribution. A VLA monitoring program is currently being undertaken to measure the differential time delay.

Using Very Long Baseline Interferometry we have searched a sample of 300 compact radio sources for examples of multiple imaging produced by gravitational lensing; no multiple images were found with separations in the angular range 1.5--50 milliarcsec. This null result allows us to place a limit on the cosmological abundance of intergalactic supermassive compact objects in the mass range $\sim 10^{6}$ to $\sim 10^{8}$M$_{\odot}$; such objects cannot make up more than $\sim 1%$ of the closure density (95% confidence). A uniformly distributed population of supermassive black holes forming soon after the Big Bang do not, therefore, contribute significantly to the dark matter content of the Universe.

Radio monitoring of the gravitational lens system B0218+357 reveals it to be a highly variable source with variations on timescales of a few days correlated in both images. This shows that the variability is intrinsic to the background lensed source and suggests that similar variations in other intraday variable sources can also be intrinsic in origin.

We present VLA and MERLIN monitoring data of the JVAS gravitational lens system B1030+074. The system was monitored with the VLA from February 1998 to October 1998 at 8.4-GHz during which the VLA was at its A, BnA and B configuration. The 47 epochs of observations have an average spacing of approximately 5 days. Ten MERLIN snapshots were obtained in the L-band (1.7 GHz) during the months of April, May and June 1998. Preliminary light curves of the two components of the lens system obtained from the VLA data indicate that during the period of the monitoring the A flux density showed a steady decrease. No changes are observed in the B light curve.

Gravitational lenses offer the possibility of accurately determining the Hubble parameter (H_0) over cosmological distances, and B0218+357 is one of the most promising systems for an application of this technique. In particular this system has an accurately measured time delay (10.5+/-0.4 d; Biggs et al. 1999) and preliminary mass modelling has given a value for H_0 of 69 +13/-19 km/s/Mpc. The error on this estimate is now dominated by the uncertainty in the mass modelling. As this system contains an Einstein ring it should be possible to constrain the model better by imaging the ring at high resolution. To achieve this we have combined data from MERLIN and the VLA at a frequency of 5 GHz. In particular MERLIN has been used in multi-frequency mode in order to improve substantially the aperture coverage of the combined data set. The resulting map is the best that has been made of the ring and contains many new and interesting features. Efforts are currently underway to exploit the new data for lensing constraints using the LensClean algorithm (Kochanek & Narayan 1992).

We present the strategy and status of a gravitational lens search for multiple imaging with angular separations between 6'' and 15'' with the combined JVAS and CLASS data set of ~15,000 flat-spectrum radio sources. Currently all but one lens candidate have been rejected.

We present the results of an unbiased radio search for gravitational lensing events with image separations between 15 and 60 arcsec, which would be associated with clusters of galaxies with masses >10^{13-14}M_{\sun}. A parent population of 1023 extended radio sources stronger than 35 mJy with stellar optical identifications was selected using the FIRST radio catalogue at 1.4 GHz and the APM optical catalogue. The FIRST catalogue was then searched for companions to the parent sources stronger than 7 mJy and with separation in the range 15 to 60 arcsec. Higher resolution observations of the resulting 38 lens candidates were made with the VLA at 1.4 GHz and 5 GHz, and with MERLIN at 5 GHz in order to test the lens hypothesis in each case. None of our targets was found to be a gravitational lens system. These results provide the best current constraint on the lensing rate for this angular scale, but improved calculations of lensing rates from realistic simulations of the clustering of matter on the relevant scales are required before cosmologically significant constraints can be derived from this null result. We now have an efficient, tested observational strategy with which it will be possible to make an order-of-magnitude larger unbiased search in the near future.

High resolution MERLIN observations of a newly-discovered four-image gravitational lens system, B0128+437, are presented. The system was found after a careful re-analysis of the entire CLASS dataset. The MERLIN observations resolve four components in a characteristic quadruple-image configuration; the maximum image separation is 542 mas and the total flux density is 48 mJy at 5 GHz. A best-fit lens model with a singular isothermal ellipsoid results in large errors in the image positions. A significantly improved fit is obtained after the addition of a shear component, suggesting that the lensing system is more complex and may consist of multiple deflectors. The integrated radio spectrum of the background source indicates that it is a GigaHertz-Peaked Spectrum (GPS) source. It may therefore be possible to resolve structure within the radio images with deep VLBI observations and thus better constrain the lensing mass distribution.

We present results from VLBA+Effelsberg 1.7 GHz, VLBA 8.4 GHz and VLBA 15 GHz observations of the JVAS gravitational lens system B1030+074. The VLBA+Effelsberg 1.7 GHz data with 3 mas resolution reveal a detailed 60 mas jet structure in the strong A component of the lens system and maybe a first hint of the corresponding jet structure in the faint component B. As we go to the higher resolution 8.4 GHz and 15 GHz data we can see more of the inner structure of the A core, which is resolved, while the B core still remains unresolved.

We present the discovery of CLASS B0739+366, a new gravitational lens system from the Cosmic Lens All-Sky Survey. Radio imaging of the source with the Very Large Array (VLA) shows two compact components separated by $0\farcs54$, with a flux density ratio of $\sim$ 6:1. High-resolution follow-up observations using the Very Long Baseline Array (VLBA) at 1.7 GHz detect weak, parity-reversed jet emission from each of the radio components. Hubble Space Telescope NICMOS F160W observations detect infrared counterparts to the lensed images, as well as an extended object between them which we identify as the lensing galaxy. Redshifts for the galaxy and lensed source have not yet been obtained. For typical lens and source redshifts of $z=0.5$ and $z=1.5$, respectively, preliminary mass modeling predicts a time delay of $\sim7h^{-1}$ days in a flat $\Omega_{M}=1.0$ universe. The small predicted time delay and weak radio components will make CLASS B0739+366 a challenging target for Hubble constant determination.

We present radio, optical, near-infrared and spectroscopic observations of the source B0827+525. We consider this source as the best candidate from the Cosmic Lens All-Sky Survey (CLASS) for a `dark lens' system or binary radio-loud quasar. The system consists of two radio components with somewhat different spectral indices, separated by 2.815 arcsec. VLBA observations show that each component has substructure on a scale of a few mas. A deep K-band exposure with the W.M.Keck-II Telescope reveals emission near both radio components. The K-band emission of the weaker radio component appears extended, whereas the emission from the brighter radio component is consistent with a point source. Hubble Space Telescope F160W-band observations with the NICMOS instrument confirms this. A redshift of 2.064 is found for the brighter component, using the LRIS instrument on the W.M.Keck-II Telescope. The probability that B0827+525 consists of two unrelated compact flat-spectrum radio sources is ~3%, although the presence of similar substructure in both component might reduce this. We discuss two scenarios to explain this system: (i) CLASS B0827+525 is a `dark lens' system or (ii) B0827+525 is a binary radio-loud quasar. B0827+525 has met all criteria that thus far have in 100% of the cases confirmed a source as an indisputable gravitational lens system. Despite this, no lens galaxy has been detected with m_F160W<=23 mag. Hence, we might have found the first binary radio-loud quasar. At this moment, however, we feel that the `dark lens' hypothesis cannot yet be fully excluded.

The complex ten-component gravitational lens system B1933+503 has been monitored with the VLA during the period February to June 1998 with a view to measuring the time delay between the four compact components and hence to determine the Hubble parameter. Here we present the results of an `A' configuration 8.4-GHz monitoring campaign which consists of 37 epochs with an average spacing of 2.8 days. The data have yielded light curves for the four flat-spectrum radio components (components 1, 3, 4 and 6). We observe only small flux density changes in the four flat-spectrum components which we do not believe are predominantly intrinsic to the source. Therefore the variations do not allow us to determine the independent time delays in this system. However, the data do allow us to accurately determine the flux density ratios between the four flat-spectrum components. These will prove important as modelling constraints and could prove crucial in future monitoring observations should these data show only a monotonic increase or decrease in the flux densities of the flat-spectrum components.

Spectroscopic observations of a sample of 42 flat-spectrum radio sources from the Cosmic Lens All-Sky Survey (CLASS) have yielded a mean redshift of $<z> = 1.27$ with an RMS spread of 0.95, at a completeness level of 64%. The sample consists of sources with a 5-GHz flux density of 25-50 mJy, making it the faintest flat-spectrum radio sample for which the redshift distribution has been studied. The spectra, obtained with the Willam Herschel Telescope (WHT), consist mainly of broad-line quasars at $z>1$ and narrow-line galaxies at $z<0.5$. Though the mean redshift of flat-spectrum radio sources exhibits little variation over more than two orders of magnitude in radio flux density, there is evidence for a decreasing fraction of quasars at weaker flux levels. In this paper we present the results of our spectroscopic observations, and discuss the implications for constraining cosmological parameters with statistical analyses of the CLASS survey.

We present an overview of all the observations (radio - VLA, MERLIN, VLBA,EVN - and optical - WFPC2 and NICMOS -) that were initially used to confirm the gravitational lens nature of the double JVAS system B1030+074. Since the 1.56 arcsec system showed some first indication of variability it has been monitored with the VLA and MERLIN to confirm its variable nature. We also present new VLBA observations of the lens system at 1.7 GHz that have unveiled detailed structure of the jet in the strong component and first detection of the jet in the faint component.

We present Hubble Space Telescope (HST) infrared images of four gravitational lens systems from the JVAS/CLASS gravitational lens survey and compare the new infrared HST pictures with previously published WFPC2 HST optical images and radio maps. Apart from the wealth of information that we get from the flux ratios and accurate positions and separations of the components of the lens systems that we can use as inputs for better constraints on the lens models we are able to discriminate between reddening and optical/radio microlensing as the possible cause of differences observed in the flux ratios of the components across the three wavelength bands. Substantial reddening has been known to be present in the lens system B1600+434 and has been further confirmed by the present infrared data. In the two systems B0712+472 and B1030+074 microlensing has been pinpointed down as the main cause of the flux ratio discrepancy both in the optical/infrared and in the radio, the radio possibly caused by the substructure revealed in the lensing galaxies. In B0218+357 however the results are still not conclusive. If we are actually seeing the two "true" components of the lens system then the flux ratio differences are attributed to a combination of microlensing and reddening or alternatively due to some variability in at least one of the images. Otherwise the second "true" component of B0218+357 maybe completely absorbed by a molecular cloud and the anomalous flux density ratios and large difference in separation between the optical/infrared and radio that we see can be explained by emission from either a foreground object or from part of the lensing galaxy.

The complex 10-component gravitational lens B1933+503 has been monitored with the VLA during the period February to June 1998. Here we present the results of this A-configuration 8.4 GHz monitoring campaign. The 37 epochs of observations have an average spacing of 3 days. The data have yielded light curves for the four flat-spectrum radio components (components 1, 3, 4 and 6). We find no significant flux density changes in any of the four flat spectrum components.

This poster presents a new 5 GHz combined MERLIN/VLA map of B0218+357 which shows for the first time believable substructure in the Einstein ring. This will now be exploited for further constraints on the model which presently dominates the error on the estimate of H_0 derived from the time delay (10.5+/-0.4 days) measured for this system.

We present the results of a program to monitor the four-image gravitational lens B1608+656 with the VLA. The system was observed over a seven month period from 1996 October to 1997 May. The 64 epochs of observation have an average spacing of 3.6~d. The light curves of the four images of the background source show that the flux density of the background source has varied at the ~5% level. We measure time delays in the system based on common features that are seen in all four light curves. The three independent time delays in the system are found to be Delta t_{BA} = 31 +/- 7~d, Delta t_{BC} = 36 +/- 7~d, and Delta t_{BD} = 76^{+9}_{-10}~d at 95% confidence. This is the first gravitational lens system for which three independent time delays have been measured. A companion paper presents a mass model for the lensing galaxy which correctly reproduces the observed image positions, flux density ratios, and time delay ratios. The last condition is crucial for determining H_0 with a four-image lens. We combine the time delays with the model to obtain a value for the Hubble constant of H_0 = 59^{+8}_{-7} km/s/Mpc at 95% confidence (statistical) for (Omega_M, Omega_{Lambda}) = (1,0). In addition, there is an estimated systematic uncertainty of +/- 15 km/s/Mpc from uncertainties in modeling the radial mass profiles of the lensing galaxies. The value of H_0 presented in this paper is comparable to recent measurements of H_0 from the gravitational lenses 0957+561, PG1115+080, B0218+357, and PKS1830-211.

The third phase of the Cosmic Lens All-Sky Survey (CLASS) has recently been completed, bringing the total number of sources imaged to over 15000 in the CLASS and JVAS combined survey. In the VLA observations carried out in March and April of 1998, two new candidate lensed systems were discovered: CLASS B1152+199 and B1359+154. B1152+199 is a 1.6 arcsecond double, with a background quasar at z=1.019 lensed by a foreground galaxy at z=0.439. The relatively flat radio spectra of the lensed images, combined with a previous ROSAT detection of the source, make B1152+199 a strong candidate for time delay studies at both radio and X-ray wavelengths. B1359+154 is a quadruply lensed quasar at z=3.235, with a maximum image separation of 1.7 arcseconds. As yet, the redshift of the lensing object in this system is undetermined. The steep spectral index of the source suggests that B1359+154 will not exhibit strong variability, and is therefore unlikely to be useful for determining the Hubble constant from measured time delays.

We present constraints on the cosmological constant $\lambda_{0}$ from gravitational lensing statistics of the Jodrell Bank-VLA Astrometric Survey (JVAS). Although this is the largest gravitational lens survey which has been analysed, cosmological constraints are only comparable to those from optical surveys. This is due to the fact that the median source redshifts of JVAS are lower, which leads to both relatively fewer lenses in the survey and a weaker dependence on the cosmological parameters. Although more approximations have to be made than is the case for optical surveys, the consistency of the results with those from optical gravitational lens surveys and other cosmological tests indicate that this is not a major source of uncertainty in the results. However, joint constraints from a combination of radio and optical data are much tighter. Thus, a similar analysis of the much larger Cosmic Lens All-Sky Survey should provide even tighter constraints on the cosmological constant, especially when combined with data from optical lens surveys. At 95% confidence, our lower and upper limits on $\lambda_{0}-\Omega_{0}$, using the JVAS lensing statistics information alone, are respectively -2.69 and 0.68. For a flat universe, these correspond to lower and upper limits on \lambda_{0} of respectively -0.85 and 0.84. Using the combination of JVAS lensing statistics and lensing statistics from the literature as discussed in Quast & Helbig (Paper I) the corresponding $\lambda_{0}-\Omega_{0}$ values are -1.78 and 0.27. For a flat universe, these correspond to lower and upper limits on $\lambda_{0}$ of respectively -0.39 and 0.64.

The JVAS (Jodrell Bank-VLA Astrometric Survey) and CLASS (Cosmic Lens All-Sky Survey) are well-defined surveys containing about ten thousand flat-spectrum radio sources. For many reasons, flat-spectrum radio sources are particularly well-suited as a population from which one can obtain unbiased samples of gravitational lenses. These are by far the largest gravitational (macro)lens surveys, and particular attention was paid to constructing a cleanly-defined sample for the survey itself and for the underlying luminosity function. Here we present the constraints on cosmological parameters, particularly the cosmological constant, derived from JVAS and combine them with constraints from optical gravitational lens surveys, `direct' measurements of $\Omega_{0}$, $H_{0}$ and the age of the universe, and constraints derived from CMB anisotropies, before putting this final result into the context of the latest results from other, independent cosmological tests.

We present the discovery of a new gravitational lens system with two compact radio images separated by 0.701+-0.001 arcsec. The lens system was discovered in the Cosmic Lens All Sky Survey (CLASS) as a flat spectrum radio source. Both radio components show structure in a VLBA 8.4 GHz radio image. No further extended structure is seen in either the VLA, MERLIN or VLBA images. Hubble Space Telescope (HST) WFPC2 images in F555W and F814W show two extended objects close to the radio components, which we identify as two lens galaxies. Their colours and mass-to-light ratios seem to favour two late-type spiral galaxies at relatively high redshifts (z_d>~0.5). Faint emission is also detected at positions corresponding to the radio images. A two-lens mass model can explain the observed VLBA structure. The best fit model has a reduced chi^2 of 1.1. The relative positions of the VLBA subcomponents are reproduced within 0.08 mas, the flux density ratios within 0.19. We also reproduce the position angle and separation of the two VLBA subcomponents in A and B within the observational errors, which we consider strong evidence for the validity of the lens model. Moreover, we find a surface density axis ratio of 0.74^{+0.10}_{-0.12} for the primary lens (G1), consistent with the surface brightness axis ratio of 0.69+-0.15. Also, the surface density position angle of 4.9^{+28.2}_{-22.4} degree of G1 compares well with the -6+-13 degree position angle of the surface brightness distribution. The errors indicate the 99 per cent confidence interval.

Measurement of the time delay between multiple images of a gravitational lens system is potentially an accurate method of determining the Hubble constant over cosmological distances. One of the most promising candidates for an application of this technique is the system B0218+357 which was found in the Jodrell Bank/VLA Astrometric Survey (JVAS). This system consists of two images of a compact radio source, separated by 335 milliarcsec, and an Einstein ring which can provide a strong constraint on the mass distribution in the lens. We present here the results of a three-month VLA monitoring campaign at two frequencies. The data are of high quality and both images show clear variations in total flux density, percentage polarization and polarization position angle at both frequencies. The time delay between the variations in the two images has been calculated using a chi-squared minimization to be 10.5+/-0.4 days at 95 per cent confidence, with the error being derived from Monte-Carlo simulations of the light curves. Although mass modelling of the system is at a preliminary stage, taking the lensing galaxy to be a singular isothermal ellipsoid and using the new value for the time delay gives a value for the Hubble constant of 69+13/-19 km/s/Mpc, again at 95 per cent confidence.

We have discovered a new gravitational lens in the Cosmic Lens All-Sky Survey (CLASS). B2045+265 is a four-image system with a maximum separation of 1.9 arcsec. A fifth radio component is detected, but its radio spectrum and its positional coincidence with infrared emission from the lensing galaxy strongly suggests that it is the radio core of the lensing galaxy. This implies that the B2045+265 system consists of a flat-spectrum radio source being lensed by another flat-spectrum radio source. Infrared images detect the lensed images of the background source and the lensing galaxy. The lensed images have relative positions and flux densities that are consistent with those seen at radio wavelengths. Spectra of the system reveal a lens redshift of z_l=0.8673 and a source redshift of z_s=1.28. The image splitting and system redshifts imply that the projected mass inside the Einstein radius of the lensing galaxy is M_E=4.7x10^{11} h^{-1} M_sun. An estimate of the light emitted inside the Einstein radius gives a mass-to-light ratio in the rest frame B band of (M/L_B)_E = 20 h (M/L_B)_sun. Both the mass and mass-to-light ratio are higher than what is seen in nearby Sa galaxies. In fact, the implied rotation velocity for the lensing galaxy is two to three times higher than what is seen in nearby spirals. The large projected mass inside the Einstein ring radius may be the result of a significant amount of dark matter in the system, perhaps from a compact group of galaxies; however, it may also arise from a misidentification of the source redshift. A simple model of the gravitational potential of the lens reproduces the image positions well, but further modeling is required to satisfy the constraints from the image flux density ratios. With further observations and modeling, this lens may yield an estimate of H_0.

NICMOS observations of the complex gravitational lens system B1933+503 reveal infrared counterparts to two of the inverted spectrum radio images. The infrared images have arc-like structures. The corresponding radio images are also detected in a VLBA map made at 1.7 GHz with a resolution of 6 mas. We fail to detect two of the four inverted radio spectrum components with the VLBA even though they are clearly visible in a MERLIN map at the same frequency at a different epoch. The absence of these two components could be due to rapid variability on a time-scale less than the time delay, or to broadening of the images during propagation of the radio waves through the ISM of the lensing galaxy to an extent that they fall below the surface brightness detectability threshold of the VLBA observations. The failure to detect the same two images with NICMOS is probably due to extinction in the ISM of the lensing galaxy.

In a recent paper, Hawkins (1997) argues on the basis of statistical studies of double-image gravitational lenses and lens candidates that a large population of dark lenses exists and that these outnumber galaxies with more normal mass-to-light ratios by a factor of 3:1. If correct, this is a very important result for many areas of astronomy including galaxy formation and cosmology. In this paper we discuss our new radio-selected gravitational lens sample, JVAS/CLASS, in order to test and constrain this proposition. We have obtained ground-based and HST images of all multiple-image lens systems in our sample and in 12 cases out of 12 we find the lensing galaxies in the optical and/or near infrared. Our success in finding lensing galaxies creates problems for the dark lens hypothesis. If it is to survive, ad hoc modifications seem to be necessary: only very massive galaxies (more than about one trillion solar masses) can be dark, and the cutoff in mass must be sharp. Our finding of lens galaxies in all the JVAS/CLASS systems is complementary evidence which supports the conclusion of Kochanek et al. (1997) that many of the wide-separation optically-selected pairs are physically distinct quasars rather than gravitational lens systems.

We report the discovery of a new double image gravitational lens system B1030+074 which was found during the Jodrell Bank - VLA Astrometric Survey (JVAS). We have collected extensive radio data on the system using the VLA, MERLIN, the EVN and the VLBA and optical observations using WFPC2 on the HST. The lensed images are separated by 1.56 arcseconds and their flux density ratio at centimetric wavelengths is approximately 14:1 although the ratio is slightly frequency dependent and the images appear to be time variable. The HST pictures show both the lensed images and the lensing galaxy close to the weaker image. The lensing galaxy has substructure which could be a spiral arm or an interacting galaxy.

We report the discovery of the most complex arcsec-scale radio gravitational lens system yet known. B1933+503 was found during the course of the CLASS survey and MERLIN and VLA radio maps reveal up to 10 components. Four of these are compact and have flat spectra; the rest are more extended and have steep spectra. The background lensed object appears to consist of a flat spectrum core (quadruply imaged) and two compact "lobes" symmetrically disposed relative to the core. One of the lobes is quadruply imaged while the other is doubly imaged. An HST observation of the system with the WFPC2 shows a galaxy with an axial ratio of 0.5, but none of the images of the background object are detected. A redshift of 0.755 has been measured for the lens galaxy.

The Caltech-Jodrell Bank VLBI surveys of bright extragalactic radio sources north of declination 35 deg were carried out between 1990 and 1995 using the Mark-II system, achieving images with a resolution of about 1 mas at 5 GHz. The CJ1 survey (together with the older `PR' survey) includes 200 objects with 5 GHz flux density greater than 0.7 Jy; the CJ2 survey includes 193 flat-spectrum sources with 5 GHz flux density greater than 0.35 Jy; and we have defined a complete flux-density limited sample, CJF, of 293 flat-spectrum sources stronger than 0.35 Jy. We summarize the definition of the samples and the VLBI, VLA, MERLIN, and optical observations, and present some highlights of the astrophysical results. These include: (1) superluminal motion and cosmology; (2) morphology and evolution of the `compact symmetric objects' (CSOs); (3) two-sided motion in some CSOs; (4) the angular-size - redshift diagram; (5) misalignment of parsec-scale and kiloparsec-scale jets.

We report the discovery, using NICMOS on the Hubble Space Telescope, of an arcsecond-diameter Einstein ring in the gravitational lens system B1938+666. The lensing galaxy is also detected, and is most likely an early-type. Modelling of the ring is presented and compared with the radio structure from MERLIN maps. We show that the Einstein ring is consistent with the gravitational lensing of an extended infrared component, centred between the two radio components.

We report the discovery of a radio loud quasar with a redshift of 4.30. This object, which is the first radio selected quasar with a redshift greater than four, was discovered during an observational investigation into the evolution of the luminosity function of radio loud quasars. Here we describe results based on a sample of ~300, flat spectrum radio sources with S(5GHz) > 200mJy. In this study, rather than carry out an indiscriminate redshift campaign on all the radio sources, we have used the APM POSS catalogue to preselect a subset of stellar-like optical counterparts with red optical colours. Such a subsample is expected to contain a high fraction of high redshift quasars. 10 of the ~300 sources were selected for follow-up optical spectroscopy and three of these are identified as quasars with redshifts greater than 3.0. One of the radio sources, GB1508+5714, is a quasar with z=4.30 and m(R)~19.