We present Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimetre Array (ALMA) observations of SDSS J0924+0219, a z = 1.524 radio-quiet lensed quasar with an intrinsic radio flux density of about 3 micro-Jy. The four lensed images are clearly detected in the radio continuum and the CO(5-4) line, whose centroid is at z = 1.5254 +/- 0.0001, with a marginal detection in the submillimetre continuum. The molecular gas displays ordered motion, in a structure approximately 1--2.5 kpc in physical extent, with typical velocities of 50-100 km/s. Our results are consistent with the radio emission being emitted from the same region, but not with a point source of radio emission. SDSS J0924+0219 shows an extreme anomaly in the flux ratios of the two merging images in the optical continuum and broad emission lines, suggesting the influence of microlensing by stars in the lensing galaxy. We find the flux ratio in the radio, submillimetre continuum and CO lines to be slightly greater than 1 but much less than that in the optical, which can be reproduced with a smooth galaxy mass model and an extended source. Our results, supported by a microlensing simulation, suggest that the most likely explanation for the optical flux anomaly is indeed microlensing.

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

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

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

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged)

We present e-MERLIN and EVN observations which reveal unambiguous jet activity within radio quiet quasar HS~0810+2554. With an intrinsic flux density of 880~nJy, this is the faintest radio source ever imaged. The findings present new evidence against the idea that radio loud and radio quiet quasars are powered by different underlying radio emission mechanisms, showing instead that the same AGN mechanism can operate as the dominant source of radio emission even in the very lowest radio luminosity quasars. Thanks to strong gravitational lensing, our source is not only visible, but with VLBI is imaged to a scale of just 0.27~pc: the highest ever resolution image of a radio quiet quasar. Brightness temperatures of at least $8.4\times 10^6$~K are associated with two highly compact components, and subsequent modelling of the lensed system has revealed that the components are linearly aligned on opposing sides of the optical quasar core, with the typical morphology of a compact symmetric object (CSO). Given that this source has been found to fall on the radio--FIR correlation, we suggest that the radio--FIR correlation cannot always be used to rule out AGN activity in favour of star-formation activity. The correlation -- or at least its scatter -- may conceal the coexistence of kinetic and radiative feedback modes in AGN. Modelling of the lensing mass itself points to a non-smooth mass distribution, hinting at the presence of dark matter substructure which has manifested as astrometric perturbations of the VLBI lensed images, posing no threat to the CDM paradigm.

Determining the star-forming properties of radio-quiet quasars is important for understanding the co-evolution of star formation and black hole accretion. Here, we present the detection of the gravitationally-lensed radio-quiet quasars SDSS J1055+4628, SDSS J1313+5151 and SBS 1520+530 at 144 MHz that fall in the HETDEX Spring Field targeted in the LOFAR Two-metre Sky Survey (LoTSS) first full data release. We compare their radio and far-infrared luminosities relative to the radio-infrared correlation and find that their radio luminosities can be explained by star formation. The implied star formation rates derived from their radio and infrared luminosities are between 20 and 300 $\rm{M_{\odot}~yr^{-1}}$. These detections represent the first study of gravitationally lensed sources with LOFAR, opening a new frequency window for investigating the star-forming properties of high-redshift quasars at radio wavelengths. We consider the implications for future data releases and estimate that many of the objects in our parent sample will be detected during LoTSS, significantly increasing the fraction of gravitationally lensed radio-quiet quasars with radio detections.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168 MHz survey of the Northern sky with diverse and ambitious science goals. Many of the scientific objectives of LoTSS rely upon, or are enhanced by, the association or separation of the sometimes incorrectly catalogued radio components into distinct radio sources, and the identification and characterisation of the optical counterparts to these sources. Here we present the source associations and optical and/or IR identifications for sources in the first data release, which are made using a combination of statistical techniques and visual association and identification. We document in detail the colour- and magnitude-dependent likelihood ratio method used for statistical identification as well as the Zooniverse project, called LOFAR Galaxy Zoo, used for the visual classification. We describe the process used to select which of these two different methods is most appropriate for each LoTSS source. The final LoTSS-DR1-IDs value-added catalogue presented contains 318,520 radio sources, of which 231,716 (73%) have optical and/or IR identifications in Pan-STARRS and WISE. The value-added catalogue is available online at https://lofar-surveys.org/, as part of this data release.

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

The imminent advent of very large-scale optical sky surveys, such as Euclid and LSST, makes it important to find efficient ways of discovering rare objects such as strong gravitational lens systems, where a background object is multiply gravitationally imaged by a foreground mass. As well as finding the lens systems, it is important to reject false positives due to intrinsic structure in galaxies, and much work is in progress with machine learning algorithms such as neural networks in order to achieve both these aims. We present and discuss a Support Vector Machine (SVM) algorithm which makes use of a Gabor filterbank in order to provide learning criteria for separation of lenses and non-lenses, and demonstrate using blind challenges that under certain circumstances it is a particularly efficient algorithm for rejecting false positives. We compare the SVM engine with a large-scale human examination of 100000 simulated lenses in a challenge dataset, and also apply the SVM method to survey images from the Kilo-Degree Survey.

Cosmological parameter constraints from observations of time-delay lenses are becoming increasingly precise. However, there may be significant bias and scatter in these measurements due to, among other things, the so-called mass-sheet degeneracy. To estimate these uncertainties, we analyze strong lenses from the largest EAGLE hydrodynamical simulation. We apply a mass-sheet transformation to the radial density profiles of lenses, and by selecting lenses near isothermality, we find that the bias on H0 can be reduced to 5% with an intrinsic scatter of 10%, confirming previous results performed on a different simulation data set. We further investigate whether combining lensing observables with kinematic constraints helps to minimize this bias. We do not detect any significant dependence of the bias on lens model parameters or observational properties of the galaxy, but depending on the source--lens configuration, a bias may still exist. Cross lenses provide an accurate estimate of the Hubble constant, while fold (double) lenses tend to be biased low (high). With kinematic constraints, double lenses show bias and intrinsic scatter of 6% and 10%, respectively, while quad lenses show bias and intrinsic scatter of 0.5% and 10%, respectively. For lenses with a reduced $\chi^2 > 1$, a power-law dependence of the $\chi^2$ on the lens environment (number of nearby galaxies) is seen. Lastly, we model, in greater detail, the cases of two double lenses that are significantly biased. We are able to remove the bias, suggesting that the remaining biases could also be reduced by carefully taking into account additional sources of systematic uncertainty.

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

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

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

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($\alpha < -1$; $S_{\nu} \sim \nu^{\alpha}$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

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

We present new VLA 22-GHz and e-MERLIN 5-GHz observations of CLASS B1030+074, a two-image strong gravitational lens system whose background source is a compact flat-spectrum radio quasar. In such systems we expect a third image of the background source to form close to the centre of the lensing galaxy. The existence and brightness of such images is important for investigation of the central mass distributions of lensing galaxies, but only one secure detection has been made so far in a galaxy-scale lens system. The noise levels achieved in our new B1030+074 images reach 3 microJy/beam and represent an improvement in central image constraints of nearly an order of magnitude over previous work, with correspondingly better resulting limits on the shape of the central mass profile of the lensing galaxy. Simple models with an isothermal outer power law slope now require either the influence of a central supermassive black hole, or an inner power law slope very close to isothermal, in order to suppress the central image below our detection limit. Using the central mass profiles inferred from light distributions in Virgo galaxies, moved to z=0.5, and matching to the observed Einstein radius, we now find that 45% of such mass profiles should give observable central images, 10% should give central images with a flux density still below our limit, and the remaining systems have extreme demagnification produced by the central SMBH. Further observations of similar objects will therefore allow proper statistical constraints to be placed on the central properties of elliptical galaxies at high redshift.

We report the discovery of a gravitationally lensed hyperluminous infrared galaxy (L_IR~10^13 L_sun) with strong radio emission (L_1.4GHz~10^25 W/Hz) at z=2.553. The source was identified in the citizen science project SpaceWarps through the visual inspection of tens of thousands of iJKs colour composite images of Luminous Red Galaxies (LRGs), groups and clusters of galaxies and quasars. Appearing as a partial Einstein ring (r_e~3") around an LRG at z=0.2, the galaxy is extremely bright in the sub-millimetre for a cosmological source, with the thermal dust emission approaching 1 Jy at peak. The redshift of the lensed galaxy is determined through the detection of the CO(3-2) molecular emission line with the Large Millimetre Telescope's Redshift Search Receiver and through [OIII] and H-alpha line detections in the near-infrared from Subaru/IRCS. We have resolved the radio emission with high resolution (300-400 mas) eMERLIN L-band and JVLA C-band imaging. These observations are used in combination with the near-infrared imaging to construct a lens model, which indicates a lensing magnification of ~10x. The source reconstruction appears to support a radio morphology comprised of a compact (<250 pc) core and more extended component, perhaps indicative of an active nucleus and jet or lobe.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

We present VLA detections of radio emission in four four-image gravitational lens systems with quasar sources: HS0810+2554, RXJ0911+0511, HE0435$-$1223 and SDSSJ0924+0219, and e-MERLIN observations of two of the systems. The first three are detected at a high level of significance, and SDSS J0924+0219 is detected. HS0810+2554 is resolved, allowing us for the first time to achieve 10-mas resolution of the source frame in the structure of a radio quiet quasar. The others are unresolved or marginally resolved. All four objects are among the faintest radio sources yet detected, with intrinsic flux densities in the range 1-5$\mu$Jy; such radio objects, if unlensed, will only be observable routinely with the Square Kilometre Array. The observations of HS0810+2554, which is also detected with e-MERLIN, strongly suggest the presence of a mini-AGN, with a radio core and milliarcsecond scale jet. The flux densities of the lensed images in all but HE0435-1223 are consistent with smooth galaxy lens models without the requirement for smaller-scale substructure in the model, although some interesting anomalies are seen between optical and radio flux densities. These are probably due to microlensing effects in the optical.

Lens modeling of resolved image data has advanced rapidly over the past two decades. More recently pixel-based approaches, wherein the source is reconstructed on an irregular or adaptive grid, have become popular. Generally, the source reconstruction takes place in a Bayesian framework and is guided by a set of sensible priors. We discuss the integration of a shapelets-based method into a Bayesian framework and quantify the required regularization. In such approaches, the source is reconstructed analytically, using a subset of a complete and orthonormal set of basis functions, known as shapelets. To calculate the flux in an image plane pixel, the pixel is split into two or more triangles (depending on the local magnification), and each shapelet basis function is integrated over the source plane. Source regularization (enforcement of priors on the source) can also be performed analytically. This approach greatly reduces the number of source parameters from the thousands to hundreds and results in a posterior probability distribution that is much less noisy than pixel-based approaches.

Strong gravitational lenses provide an important tool to measure masses in the distant Universe, thus testing models for galaxy formation and dark matter; to investigate structure at the Epoch of Reionization; and to measure the Hubble constant and possibly w as a function of redshift. However, the limiting factor in all of these studies has been the currently small samples of known gravitational lenses (~10^2). The era of the SKA will transform our understanding of the Universe with gravitational lensing, particularly at radio wavelengths where the number of known gravitational lenses will increase to ~10^5. Here we discuss the technical requirements, expected outcomes and main scientific goals of a survey for strong gravitational lensing with the SKA. We find that an all-sky (3pi sr) survey carried out with the SKA1-MID array at an angular resolution of 0.25-0.5 arcsec and to a depth of 3 microJy / beam is required for studies of galaxy formation and cosmology with gravitational lensing. In addition, the capability to carryout VLBI with the SKA1 is required for tests of dark matter and studies of supermassive black holes at high redshift to be made using gravitational lensing.

We provide an overview of the science benefits of combining information from the Square Kilometre Array (SKA) and the Large Synoptic Survey Telescope (LSST). We first summarise the capabilities and timeline of the LSST and overview its science goals. We then discuss the science questions in common between the two projects, and how they can be best addressed by combining the data from both telescopes. We describe how weak gravitational lensing and galaxy clustering studies with LSST and SKA can provide improved constraints on the causes of the cosmological acceleration. We summarise the benefits to galaxy evolution studies of combining deep optical multi-band imaging with radio observations. Finally, we discuss the excellent match between one of the most unique features of the LSST, its temporal cadence in the optical waveband, and the time resolution of the SKA.

We present the results of the first strong lens time delay challenge. The motivation, experimental design, and entry level challenge are described in a companion paper. This paper presents the main challenge, TDC1, which consisted of analyzing thousands of simulated light curves blindly. The observational properties of the light curves cover the range in quality obtained for current targeted efforts (e.g.,~COSMOGRAIL) and expected from future synoptic surveys (e.g.,~LSST), and include simulated systematic errors. \nteamsA\ teams participated in TDC1, submitting results from \nmethods\ different method variants. After a describing each method, we compute and analyze basic statistics measuring accuracy (or bias) $A$, goodness of fit $\chi^2$, precision $P$, and success rate $f$. For some methods we identify outliers as an important issue. Other methods show that outliers can be controlled via visual inspection or conservative quality control. Several methods are competitive, i.e., give $|A|<0.03$, $P<0.03$, and $\chi^2<1.5$, with some of the methods already reaching sub-percent accuracy. The fraction of light curves yielding a time delay measurement is typically in the range $f = $20--40\%. It depends strongly on the quality of the data: COSMOGRAIL-quality cadence and light curve lengths yield significantly higher $f$ than does sparser sampling. Taking the results of TDC1 at face value, we estimate that LSST should provide around 400 robust time-delay measurements, each with $P<0.03$ and $|A|<0.01$, comparable to current lens modeling uncertainties. In terms of observing strategies, we find that $A$ and $f$ depend mostly on season length, while P depends mostly on cadence and campaign duration.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

I review the current state of determinations of the Hubble constant, which gives the length scale of the Universe by relating the expansion velocity of objects to their distance. There are two broad categories of measurements. The first uses individual astrophysical objects which have some property that allows their intrinsic luminosity or size to be determined, or allows the determination of their distance by geometric means. The second category comprises the use of all-sky cosmic microwave background, or correlations between large samples of galaxies, to determine information about the geometry of the Universe and hence the Hubble constant, typically in a combination with other cosmological parameters. Many, but not all, object-based measurements give $H_0$ values of around 72-74km/s/Mpc , with typical errors of 2-3km/s/Mpc. This is in mild discrepancy with CMB-based measurements, in particular those from the Planck satellite, which give values of 67-68km/s/Mpc and typical errors of 1-2km/s/Mpc. The size of the remaining systematics indicate that accuracy rather than precision is the remaining problem in a good determination of the Hubble constant. Whether a discrepancy exists, and whether new physics is needed to resolve it, depends on details of the systematics of the object-based methods, and also on the assumptions about other cosmological parameters and which datasets are combined in the case of the all-sky methods.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

We present LOFAR Low Band observations of the Bootes and 3C295 fields. Our images made at 34, 46, and 62 MHz reach noise levels of 12, 8, and 5 mJy beam$^{-1}$, making them the deepest images ever obtained in this frequency range. In total, we detect between 300 and 400 sources in each of these images, covering an area of 17 to 52 deg$^{2}$. From the observations we derive Euclidean-normalized differential source counts. The 62 MHz source counts agree with previous GMRT 153 MHz and VLA 74 MHz differential source counts, scaling with a spectral index of $-0.7$. We find that a spectral index scaling of $-0.5$ is required to match up the LOFAR 34 MHz source counts. This result is also in agreement with source counts from the 38 MHz 8C survey, indicating that the average spectral index of radio sources flattens towards lower frequencies. We also find evidence for spectral flattening using the individual flux measurements of sources between 34 and 1400 MHz and by calculating the spectral index averaged over the source population. To select ultra-steep spectrum ($\alpha < -1.1$) radio sources, that could be associated with massive high redshift radio galaxies, we compute spectral indices between 62 MHz, 153 MHz and 1.4 GHz for sources in the Bo\"otes field. We cross-correlate these radio sources with optical and infrared catalogues and fit the spectral energy distribution to obtain photometric redshifts. We find that most of these ultra-steep spectrum sources are located in the $ 0.7 \lesssim z \lesssim 2.5$ range.

Obtaining lensing time delay measurements requires long-term monitoring campaigns with a high enough resolution (< 1 arcsec) to separate the multiple images. In the radio, a limited number of high-resolution interferometer arrays make these observations difficult to schedule. To overcome this problem, we propose a technique for measuring gravitational time delays which relies on monitoring the total flux density with low-resolution but high-sensitivity radio telescopes to follow the variation of the brighter image. This is then used to trigger high-resolution observations in optimal numbers which then reveal the variation in the fainter image. We present simulations to assess the efficiency of this method together with a pilot project observing radio lens systems with the Westerbork Synthesis Radio Telescope (WSRT) to trigger Very Large Array (VLA) 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 instruments that have high sensitivity but limited resolution, together with an optimum usage of followup high-resolution observations from appropriate radio telescopes may in the future be useful for gravitational lensing time delay measurements by means of this new method.

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

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.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

I review the observations of gravitationally lensed quasars. These systems are important because they allow us to probe the properties of the lensing galaxies at various scales, and they also allow insights into the structures of the quasars themselves. Samples of quasar lenses also have the potential to act as cosmographic probes. These areas are described, together with observational and scientific prospects for the future.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

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.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

Deep radio observations of galaxy clusters have revealed the existence of diffuse radio sources related to the presence of relativistic electrons and weak magnetic fields in the intracluster volume. The role played by this non-thermal intracluster component on the thermodynamical evolution of galaxy clusters is debated, with important implications for cosmological and astrophysical studies of the largest gravitationally bound structures of the Universe. The low surface brightness and steep spectra of diffuse cluster radio sources make them more easily detectable at low-frequencies. LOFAR is the first instrument able to detect diffuse radio emission in hundreds of massive galaxy clusters up to their formation epoch. We present the first observations of clusters imaged by LOFAR and the huge perspectives opened by this instrument for non-thermal cluster studies.

We present results of spectral investigations of the Sco X-1 like Z-track sources Sco X-1, GX 349+2 and GX 17+2 based on Rossi-XTE observations using an extended accretion disk corona model. The results are compared with previous results for the Cyg X-2 like group: Cyg X-2, GX 340+0 and GX 5-1 and a general model for the Z-track sources proposed. On the normal branch, the Sco-like and Cyg-like sources are similar, the results indicating an increase of mass accretion rate Mdot between soft and hard apex, not as in the standard view that this increases around the Z. In the Cyg-like sources, increasing Mdot causes the neutron star temperature kT to increase from ~1 to ~2 keV. At the lower kT, the radiation pressure is small, but at the higher kT, the emitted flux of the neutron star is several times super-Eddington and the high radiation pressure disrupts the inner disk launching the relativistic jets observed on the upper normal and horizontal branches. In the Sco-like sources, the main physical difference is the high kT of more than 2 keV on all parts of the Z-track suggesting that jets are always possible, even on the flaring branch. The flaring branch in the Cyg-like sources is associated with release of energy on the neutron star consistent with unstable nuclear burning. The Sco-like sources are very different as flaring appears to be a combination of unstable burning and an increase of Mdot which makes flaring much stronger. Analysis of 15 years or RXTE ASM data on all 6 classic Z-track sources shows the high rate and strength of flaring in the Sco-like sources suggesting that continual release of energy heats the neutron star causing the high kT. A Sco X-1 observation with unusually little flaring supports this. GX 17+2 appears to be transitional between the Cyg and Sco-like types. Our results do not support the suggestion that Cyg or Sco-like nature is determined by luminosity.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

Measuring the relative mass contributions of luminous and dark matter in spiral galaxies is important for understanding their formation and evolution. The combination of a galaxy rotation curve and strong lensing is a powerful way to break the disk-halo degeneracy that is inherent in each of the methods individually. We present an analysis of the 10-image radio spiral lens B1933+503 at z_l=0.755, incorporating (1) new global VLBI observations, (2) new adaptive-optics assisted K-band imaging, (3) new spectroscopic observations for the lens galaxy rotation curve and the source redshift. We construct a three-dimensionally axisymmetric mass distribution with 3 components: an exponential profile for the disk, a point mass for the bulge, and an NFW profile for the halo. The mass model is simultaneously fitted to the kinematics and the lensing data. The NFW halo needs to be oblate with a flattening of a/c=0.33^{+0.07}_{-0.05} to be consistent with the radio data. This suggests that baryons are effective at making the halos oblate near the center. The lensing and kinematics analysis probe the inner ~10 kpc of the galaxy, and we obtain a lower limit on the halo scale radius of 16 kpc (95% CI). The dark matter mass fraction inside a sphere with a radius of 2.2 disk scale lengths is f_{DM,2.2}=0.43^{+0.10}_{-0.09}. The contribution of the disk to the total circular velocity at 2.2 disk scale lengths is 0.76^{+0.05}_{-0.06}, suggesting that the disk is marginally submaximal. The stellar mass of the disk from our modeling is log_{10}(M_{*}/M_{sun}) = 11.06^{+0.09}_{-0.11} assuming that the cold gas contributes ~20% to the total disk mass. In comparison to the stellar masses estimated from stellar population synthesis models, the stellar initial mass function of Chabrier is preferred to that of Salpeter by a probability factor of 7.2.

Gravitational lens systems containing lensed quasars are important as cosmological probes, as diagnostics of structural properties of the lensing galaxies and as tools to study the quasars themselves. The largest lensed quasar sample is the SDSS Quasar Lens Search, drawn from the Sloan Digital Sky Survey (SDSS). We are attempting to extend this survey using observations of lens candidates selected from a combination of the quasar sample from the SDSS and the UKIRT Infrared Deep Sky Survey (UKIDSS). This adds somewhat higher image quality together with a wider range of wavelength for the selection process. In previous pilot surveys we observed 5 objects, finding 2 lenses; here we present further observations of 20 objects in which we find 4 lenses, of which 2 are independently discovered in SQLS (in preparation). Following earlier work on the combination of these two surveys, we have refined our method and find that use of a colour-separation diagnostic, where we select for separations between components which appear to decrease in wavelength, is an efficient method to find lensed quasars and may be useful in ongoing and future large-scale strong lensing surveys with instruments such as Pan-STARRS and LSST. The new lenses have mostly high flux ratios, with faint secondaries buried in the lensing galaxy and typically 6-10 times less bright than the primary. Our survey brings the total number of lenses discovered in the SDSS quasar sample to 46, plus 13 lenses already known. This is likely to be up to 60-70% of the total number of lensed quasars; we briefly discuss strategies by which the rest might be found.

At very low frequencies, the new pan-European radio telescope LOFAR is opening the last unexplored window of the electromagnetic spectrum for astrophysical studies. The revolutionary APERTIF phased arrays that are about to be installed on the Westerbork radio telescope (WSRT) will dramatically increase the survey speed for the WSRT. Combined surveys with these two facilities will deeply chart the northern sky over almost two decades in radio frequency from \sim 15 up to 1400 MHz. Here we briefly describe some of the capabilities of these new facilities and what radio surveys are planned to study fundamental issues related the formation and evolution of galaxies and clusters of galaxies. In the second part we briefly review some recent observational results directly showing that diffuse radio emission in clusters traces shocks due to cluster mergers. As these diffuse radio sources are relatively bright at low frequencies, LOFAR should be able to detect thousands of such sources up to the epoch of cluster formation. This will allow addressing many question about the origin and evolution of shocks and magnetic fields in clusters. At the end we briefly review some of the first and very preliminary LOFAR results on clusters.

The Low Frequency Array (LOFAR) is under construction in the Netherlands and in several surrounding European countries. In this contribution, we describe the layout and design of the telescope, with a particular emphasis on the imaging characteristics of the array when used in its "standard imaging" mode. After briefly reviewing the calibration and imaging software used for LOFAR image processing, we show some recent results from the ongoing imaging commissioning efforts. We conclude by summarizing future prospects for the use of LOFAR in observing the little-explored low frequency Universe.

The Low Frequency Array (LOFAR) will operate between 10 and 250 MHz, and will observe the low frequency Universe to an unprecedented sensitivity and angular resolution. The construction and commissioning of LOFAR is well underway, with over 27 of the Dutch stations and five International stations routinely performing both single-station and interferometric observations over the frequency range that LOFAR is anticipated to operate at. Here, we summarize the capabilities of LOFAR and report on some of the early commissioning imaging of Cygnus A.

We present new radio observations of the large-separation gravitationally-lensed quasar SDSS J1004+4112, taken in a total of 6 hours of observations with the Extended Very Large Array (EVLA). The maps reach a thermal noise level of approximately 4microJy. We detect four of the five lensed images at the 15-35microJy level, representing a source of intrinsic flux density, after allowing for lensing magnification, of about 1microJy, intrinsically probably the faintest radio source yet detected. This reinforces the utility of gravitational lensing in potentially allowing us to study nanoJy-level sources before the advent of the SKA. In an optical observation taken three months after the radio observation, image C is the brightest image, whereas the radio map shows flux density ratios consistent with previous optical observations. Future observations separated by a time delay will give the intrinsic flux ratios of the images in this source.

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.

We have carried out a pilot study for the SCUBA-2 'All-Sky' Survey, SASSy, a wide and shallow mapping project at 850 microns, designed to find rare objects, both Galactic and extragalactic. Two distinct sets of exploratory observations were undertaken, and used to test the SASSy approach and data reduction pipeline. The first was a 0.5 by 0.5 degrees map around the nearby galaxy NGC 2559. The galaxy was easily detected at 156 mJy, but no other convincing sources are present in the map. Comparison with other galaxies with similar wavelength coverage indicates that NGC 2559 has relatively warm dust. The second observations cover 1 square degree around the W5-E HII region. As well as diffuse structure in the map, a filtering approach was able to extract 27 compact sources with signal-to-noise greater than 6. By matching with data at other wavelengths we can see that the SCUBA-2 data can be used to discriminate the colder cores. Together these observations show that the SASSy project will be able to meet its original goals of detecting new bright sources which will be ideal for follow-up observations with other facilities.

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.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

Based on the results of applying the extended ADC emission model for low mass X-ray binaries to three Z-track sources: GX340+0, GX5-1 and CygX-2, we propose an explanation of the CygnusX-2 like Z-track sources. The Normal Branch is dominated by the increasing radiation pressure of the neutron star caused by a mass accretion rate that increases between the soft apex and the hard apex. The radiation pressure continues to increase on the Horizontal Branch becoming several times super-Eddington. We suggest that this disrupts the inner accretion disk and that part of the accretion flow is diverted vertically forming jets which are detected by their radio emission on this part of the Z-track. We thus propose that high radiation pressure is the necessary condition for the launching of jets. On the Flaring Branch there is a large increase in the neutron star blackbody luminosity at constant mass accretion rate indicating an additional energy source on the neutron star. We find that there is good agreement between the mass accretion rate per unit emitting area of the neutron star mdot at the onset of flaring and the theoretical critical value at which burning becomes unstable. We thus propose that flaring in the CygnusX-2 like sources consists of unstable nuclear burning. Correlation of measurements of kilohertz QPO frequencies in all three sources with spectral fitting results leads to the proposal that the upper kHz QPO is an oscillation always taking place at the inner accretion disk edge, the radius of which increases due to disruption of the disk by the high radiation pressure of the neutron star.

The incidence of sub-galactic level substructures is an important quantity, as it is a generic prediction of high-resolution Cold Dark Matter (CDM) models which is susceptible to observational test. Confrontation of theory with observations is currently in an uncertain state. In particular, gravitational lens systems appear to show evidence for flux ratio anomalies, which are expected from CDM substructures although not necessarily in the same range of radius as observed. However, the current small samples of lenses suggest that the lens galaxies in these systems are unusually often accompanied by luminous galaxies. Here we investigate a large sample of unlensed elliptical galaxies from the COSMOS survey, and determine the fraction of objects with satellites, in excess of background counts, as a function of satellite brightness and separation from the primary object. We find that the incidence of luminous satellites within 20 kpc is typically a few tenths of one percent for satellites of a few tenths of the primary flux, comparable to what is observed for the wider but shallower SDSS survey. Although the environments of lenses in the SLACS survey are compatible with this observation, the CLASS radio survey lenses are significantly in excess of this.

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 observations of a new double-image gravitational lens system, ULAS J082016.1+081216, of image separation 2.3" and high (~6) flux ratio. The system is selected from the Sloan Digital Sky Survey spectroscopic quasar list using new high-quality images from the UKIRT Deep Sky Survey (UKIDSS). The lensed quasar has a source redshift of 2.024, and we identify the lens galaxy as a faint red object of redshift 0.803+/-0.001. Three other objects from the UKIDSS survey, selected in the same way, were found not to be lens systems. Together with the earlier lens found using this method, the SDSS-UKIDSS lenses have the potential to significantly increase the number of quasar lenses found in SDSS, to extend the survey to higher flux ratios and lower separations, and to give greater completeness which is important for statistical purposes.

Whereas considerable effort has been afforded in understanding the properties of galaxies, a full physical picture, connecting their baryonic and dark-matter content, super-massive black holes, and (metric) theories of gravity, is still ill-defined. Strong gravitational lensing furnishes a powerful method to probe gravity in the central regions of galaxies. It can (1) provide a unique detection-channel of dark-matter substructure beyond the local galaxy group, (2) constrain dark-matter physics, complementary to direct-detection experiments, as well as metric theories of gravity, (3) probe central super-massive black holes, and (4) provide crucial insight into galaxy formation processes from the dark matter point of view, independently of the nature and state of dark matter. To seriously address the above questions, a considerable increase in the number of strong gravitational-lens systems is required. In the timeframe 2010-2020, a staged approach with radio (e.g. EVLA, e-MERLIN, LOFAR, SKA phase-I) and optical (e.g. LSST and JDEM) instruments can provide 10^(2-4) new lenses, and up to 10^(4-6) new lens systems from SKA/LSST/JDEM all-sky surveys around ~2020. Follow-up imaging of (radio) lenses is necessary with moderate ground/space-based optical-IR telescopes and with 30-50m telescopes for spectroscopy (e.g. TMT, GMT, ELT). To answer these fundamental questions through strong gravitational lensing, a strong investment in large radio and optical-IR facilities is therefore critical in the coming decade. In particular, only large-scale radio lens surveys (e.g. with SKA) provide the large numbers of high-resolution and high-fidelity images of lenses needed for SMBH and flux-ratio anomaly studies.

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

We present results of a combined investigation of the spectral and kHz QPO evolution around the Z-track in GX 5-1 based on high-quality RXTE data. The Extended ADC emission model provides very good fits to the spectra, the results pointing clearly to a model for the nature of the Z-track, in agreement with previous results for the similar source GX 340+0. In this model, at the soft apex of the Z-track, the mass accretion rate Mdot is minimum and the neutron star has its lowest temperature; but as the source moves along the normal branch, the luminosity of the Comptonized emission increases, indicating that Mdot increases and the neutron star gets hotter. The measured flux f of the neutron star emission increases by a factor of ten becoming super-Eddington, and we propose that this disrupts the inner disk so forming jets. In flaring, the luminosity of the dominant Comptonized emission from the ADC is constant, while the neutron star emission increases, and we propose for the first time that flaring consists of unstable nuclear burning on the neutron star, and the measured mass accretion rate per unit area mdot at the onset of flaring agrees well with the theoretical critical value at which burning becomes unstable. There is a striking correlation between the frequencies of the kHz QPO and the ratio of the flux to the Eddington value: f/f_Edd, suggesting an explanation of the higher frequency QPO and of its variation along the Z-track. It is well known that a Keplerian orbit in the disk at this frequency corresponds to a position some distance from the neutron star; we propose that the oscillation always occurs at the inner disk edge, which moves radially outwards on the upper normal and horizontal branches as the measured increasing radiation pressure increasingly disrupts the inner disk.

We present a new model for the Z-track phenomenon, based on analysis of the spectral evolution around the Z-track in several Z-track sources, in which radiation pressure plays a major role. Increasing mass accretion rate on the normal branch causes heating of the neutron star with the emissive flux from the surface increasing by an order of magnitude to become super-Eddington at the horizontal branch where radio detection shows the presence of jets. We propose that the radiation pressure disrupts the inner disk leading to the launching of the jets. Secondly, by timing analysis of the same data we find a correlation of the frequency of kHz QPO with the emissive flux and propose that the higher frequency QPO is an oscillation at the inner disk edge which progressively moves to larger radial positions as the disk is disrupted by radiation pressure.

A complete manual search has been carried out of the list of 285423 objects, nearly all of them galaxies, identified in the COSMOS field that are brighter than I=25. Two certain and one highly probable new gravitational lenses are found, in addition to the lenses and candidate lens systems previously found by Faure et al. (2008). A further list of 112 candidate lens systems is presented. Few of these are likely to be true gravitational lens systems, most being star-forming rings or pairs of companion galaxies. It is possible to examine of order 10^6 objects by eye in a reasonable time, although reliable detection of lenses by such methods is likely to be possible only with high-resolution data. The loss of completeness involved in a rapid search is estimated as up to a factor of 2, depending on the morphology of the lens candidate.

We report the discovery of a new gravitational lens system. This object, ULAS J234311.93-005034.0, is the first to be selected by using the new UKIRT Infrared Deep Sky Survey (UKIDSS), together with the Sloan Digital Sky Survey (SDSS). The ULAS J234311.93-005034.0 system contains a quasar at redshift 0.788 which is doubly imaged, with separation 1.4". The two quasar images have the same redshift and similar, though not identical, spectra. The lensing galaxy is detected by subtracting point-spread functions from R-band images taken with the Keck telescope. The lensing galaxy can also be detected by subtracting the spectra of the A and B images, since more of the galaxy light is likely to be present in the latter. No redshift is determined from the galaxy, although the shape of its spectrum suggests a redshift of about 0.3. The object's lens status is secure, due to the identification of two objects with the same redshift together with a lensing galaxy. Our imaging suggests that the lens is found in a cluster environment, in which candidate arc-like structures, that require confirmation, are visible in the vicinity. Further discoveries of lenses from the UKIDSS survey are likely as part of this programme, due to the depth of UKIDSS and its generally good seeing conditions.

Based on the results of applying the extended ADC emission model to three Z-track sources: GX340+0, GX5-1 and CygX-2, we propose an explanation of the Z-track sources in which the Normal and Horizontal Branches are dominated by the increasing radiation pressure of the neutron star. The emitted flux becomes several times super-Eddington at the Hard Apex and Horizontal Branch and we suggest that the inner accretion disk is disrupted by this and that part of the accretion flow is diverted vertically. This position on the Z-track is exactly the position where radio emission is detected showing the presence of jets. We thus propose that high radiation pressure is a necessary condition for the launching of jets. We also show that flaring must consist of unstable nuclear burning and that the mass accretion rate per unit emitting area of the neutron star mdot at the onset of flaring agrees well with the critical theoretical value at which burning becomes unstable.

We report the discovery of a new gravitationally lensed quasar, SDSS J131339.98+515128.3, at a redshift of 1.875 with an image separation of 1.24". The lensing galaxy is clearly detected in visible-light follow-up observations. We also identify three absorption-line doublets in the spectra of the lensed quasar images, from which we measure the lens redshift to be 0.194. Like several other known lenses, the lensed quasar images have different continuum slopes. This difference is probably the result of reddening and microlensing in the lensing galaxy. The lensed quasar was selected by correlating Sloan Digital Sky Survey (SDSS) spectroscopic quasars with Two Micron All Sky Survey (2MASS) sources and choosing quasars that show near-infrared (IR) excess. The near-IR excess can originate, for example, from the contribution of the lensing galaxy at near-IR wavelengths. We show that the near-IR excess technique is indeed an efficient method to identify lensed systems from a large sample of quasars.

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

New observations are reported of J-band spectra (1.04micron -- 1.4 micron) of three Seyfert 2 galaxies, Mkn 34, Mkn 78 and NGC 5929. In each case the spectral range includes the near-infrared lines of [FeII], [PII], HeI and Pa beta. Each Seyfert galaxy has a known radio jet, and we investigate the infrared line ratios of the nuclear and extended regions of each galaxy compared to the radio structure. In Mkn 34 there is a clear indication of an extranuclear region, probably coincident with a shock induced by the radio jet, in which [FeII] is considerably enhanced, although the nuclear emission is almost certainly the result of photoionization by the continuum of the active nucleus. Similar effects in extranuclear regions are seen in the other objects, in the case of Mkn 78 confirming recent studies by Ramos Almeida et al. A possible detection of extranuclear [PII] emission suggests, if real, that photoionization by the active nucleus is the dominant line excitation mechanism over the whole source, including the regions coincident with the radio jet.

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 new HST near-infrared polarimetry, broad and narrow-band imaging, and MERLIN 4.5GHz Multi-Frequency Synthesis radio imaging of 3C 293, a unique radio galaxy whose host is an obvious merger remnant, in an exceptionally under-dense region of space. We have discovered near-infrared, optical, and ultra-violet synchrotron emission from the jet. In the optical, the jet is mostly obscured by a dust lane, but three knots are clear in our HST NICMOS images at 1.6 and 2.0 microns, clearly aligning with features in the radio. The outer jet knot is highly polarized (~15%) at 2 microns, confirming the synchrotron emission mechanism. The radio-IR spectral index steepens significantly with distance from the nucleus, as in 3C 273 and in contrast to M 87. The inner knot is visible (with hindsight) on the WFPC2 and STIS images obtained for the earlier 3CR HST snapshot surveys. There is no [Fe II] emission seen associated with the jet, constraining the role of shock-induced ionisation by the jet. Overall there is a strong implication that the NIR jet emission is indeed synchrotron. From our NIR images, the core of the galaxy is clearly identifiable with the main feature in the western extension of the radio ``jet'' image, although no unresolved AGN component is identifiable even at K-band, consistent with an FRII-like nucleus obscured by an optically thick torus. The galaxy appears to have a single nucleus, with any multiple nuclei falling within the central </~100 pc.

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

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 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 have obtained near-infrared spectroscopy of a small sample of powerful Compact Steep-Spectrum (CSS) radio sources mainly, but not exclusively, from the 3CR sample. We find no differences between the distributions in the equivalent width and luminosity of the [OIII] 5007A line for our sample and other larger, presumably older, high-redshift 3C objects, suggesting that the underlying quasar luminosity remains roughly constant as quasars age. We also find a possible broad line in 3C241, adding to recent evidence for broad lines in some radio galaxies.

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.

We use a combination of high spatial resolution optical and near-IR spectroscopic data to make a detailed study of the kinematics of the NLR gas in the near-nuclear regions of the powerful, FRII radio galaxy Cygnus A (z=0.0560), with the overall goal of placing limits on the mass of any supermassive black hole in the core. Our K-band infrared observations (0.75 arcsec seeing) -- taken with NIRSPEC on the Keck II telescope -- show a smooth rotation pattern across the nucleus in the Paschen alpha and H_2 emission lines along a slit position (PA180) close to perpendicular to the radio axis, however, there is no evidence for such rotation along the radio axis (PA105). Higher spatial resolution observations of the [OIII]5007 emission line -- taken with STIS on the Hubble Space Telescope (HST) -- confirm the general rotation pattern of the gas in the direction perpendicular to the radio axis, and provide evidence for steep velocity gradients within a radius of 0.1 arcsec of the core. The circular velocities measured from both the Keck and HST data lead to an estimate of the mass of the supermassive black hole of 2.5+/-0.7x10^9 solar masses. For the host galaxy properties of Cygnus A, this mass is consistent with the global correlations between black hole mass and host galaxy properties deduced for non-active galaxies. Therefore, despite the extreme power of its radio source and the quasar-like luminosity of its AGN, the black hole in Cygnus A is not unusually massive considering theluminosity of its host galaxy. Indeed, the estimated mass of the black hole in Cygnus A is similar to that inferred for the supermassive black hole in the FRI radio galaxy M87, despite the fact that the AGN and radio jets of Cygnus A are 2 -- 3 orders of magnitude more powerful.

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.

We present MERLIN and Hubble Space Telescope (HST) observations of the central region of the nearby radio galaxy 3C305 and use them to study the gas and dust in this object. The MERLIN observations are of neutral hydrogen (HI) absorption against the strong non-thermal 20cm continuum seen towards the central 4kpc of 3C305. Our 0.2-arcsec (160 pc) resolution observations show that the HI absorption is highly localised against the south-western radio-emission with column densities 1.9x10^21/cm2. The absorption is broad and red-shifted by 130 km/s relative to the systemic velocity. The HST images in 430nm, 702nm, [OIII] 500.7nm, [FeII] 1.64um and K-band polarization are presented. Evidence is seen for coincidence of the [FeII] emission with the knot at the end of the radio jet, suggesting the presence of shocks. We compare the optical and radio images in order to investigate the relationship between the dust and neutral gas. An unresolved (0.07-arcsec) nucleus is detected in H and K and its properties are consistent with a quasar reddened by A(V)>4. We propose that the absorption arises in a region of neutral gas and dust whose structure is complex but is broadly consistent with an inclined disc encircling, but not covering, the active galactic nucleus. A comparison of the neutral gas observations and previous emission-line observations suggests that both the neutral and ionised gas are undergoing galactic rotation towards the observer in the north-east and away from the observer in the south-west. We propose that the outflow giving rise to the radio emission has a component towards the observer in the north-east and away from the observer in the south-west.

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.

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.

We report on the search for distant radio-loud quasars in the Cosmic Lens All Sky Survey (CLASS) of flat spectrum radio sources with S_5GHz>30 mJy. Unresolved optical counterparts were selected from APM scans of POSS-I plates, with e<19.0 and red o-e>2.0 colours, in an effective area of ~6400 deg^2. Four sources were found to be quasars with z>4, of which one was previously known. This sample bridges the gap between the strong radio surveys with S_5GHz>200 mJy and the samples of radio-weak quasars that can be generated via radio observations of optically selected quasars. In addition, 4 new quasars at z>3 have been found. The selection criteria result in a success-rate of ~1:7 for radio-loud quasars at z>4, which is a significant improvement over previous studies. This search yields a surface density of 1 per 1600 deg^2, which is about a factor of ~15 lower than that found in a similar search for radio-quiet quasars at z>4. The study presented here is strongly biased against quasars beyond z>4.5, since the e-passband of the POSS-I only samples the spectra shortward of 1200 Angstrom at these redshifts.

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.

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.

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 near-IR imaging polarimetry observations of the nucleus of Cygnus A, taken with the NICMOS camera of the HST at a wavelength of 2.0 microns.These maps reveal a highly collimated region of polarized emission straddling the nucleus and extending to a radius of 1.2 arcseconds. Remarkably, this feature coincides with one, but only one, limb of the edge-brightened bicone structure seen in the total intensity image. The high degree (P_k ~ 25%) and orientation of the extended polarization feature are consistent with a scattering origin. Most plausibly, the detection of polarization along only one limb of the bicone is a consequence of intrinsic anisotropy of the near-IR continuum within the radiation cones, with the direction of maximum intensity of the near-IR radiation field significantly displaced from the direction of the radio axis. The unresolved nuclear core source is also highly polarized (P_k > 28%), with a position angle close to the perpendicular to the radio axis. Given that this high degree of nuclear polarization can only be explained in terms of dichroic extinction if the dichroic mechanism is unusually efficient in Cygnus A, it is more likely that the nuclear polarization is due to the scattering of nuclear light in an unresolved scattering region close to the AGN. In this case, the flux of the core source in the K-band is dominated by scattered rather than transmitted quasar light, and previous extinction estimates based on K-band photometry of the core substantially underestimate the true nuclear extinction.

The CERES Astronomical Archive was created initially in order to make the large amount of data that were collected from the two surveys, the Jodrell-VLA Astrometric Survey (JVAS) and the Cosmic Lens All-Sky Survey (CLASS), easily accessible to the partner members of the Consortium for European Research on Extragalactic Surveys (CERES) through the web. Here we describe the web based archive of the 15,000 flat spectrum radio sources. There is a wealth of information that one can access including downloading the raw and processed data for each of the sources. We also describe the search engines that were developed so that a) one can search through the archive for an object or a sample of objects by setting specific criteria and b) select new samples. Since new data are continually gathered on the same sources from follow up programs, an automatic update engine was created so that the new information and data can be added easily in the archive.

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

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.

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 present new observations of the gravitational lens (GL) system B1600+434, strongly suggesting that the lens is an edge-on spiral galaxy. These observations are used to constrain the mass model of the system, in particular the oblateness and velocity dispersion of the dark matter halo around the lensing galaxy. From an analytical model we find a lower limit on the halo oblateness $q_\rmn{halo}=(c/a)_{\rho}\ga0.4$; more detailed numerical models give a lower limit $q_\rmn{h}\ga0.5$. We determine an average halo velocity dispersion of $\sigma_{halo}=190\pm 15$ km/s over all non-singular isothermal elliptical (NIE) halo models. Constraining the models to larger and more massive disks, decreases this average by only 10 km/s. A lower limit of $\sigma_{halo}\ga 150$ km/s is found, even for disk masses larger then the mass inside the Einstein radius. This lower limit indicates the need for a massive dark matter halo, contributing at least half of the mass inside the Einstein radius. Time delay calculations give $(54\pm3)/h_{50}$ days for the NIE halo model and $(70\pm4)/h_{50}$ days for the modified Hubble profile (MHP) halo model. Although the time delay for both NIE and MHP halo models is well constrained on our parameter grid, it strongly depends on the halo surface density profile. We furthermore find that the presence of a flat luminous mass distribution can severely alter the statistical properties of the lens.

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.

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.