We report observations with the Atacama Large Millimetre Array (ALMA) of six submillimetre galaxies (SMGs) within 3 arcmin of the Distant Red Core (DRC) at $z=4.0$, a site of intense cluster-scale star formation, first reported by Oteo et al. (2018). We find new members of DRC in three SMG fields; in two fields, the SMGs are shown to lie along the line of sight towards DRC; one SMG is spurious. Although at first sight this rate of association is consistent with earlier predictions, associations with the bright SMGs are rarer than expected, which suggests caution when interpreting continuum over-densities. We consider the implications of all 14 confirmed DRC components passing simultaneously through an active phase of star formation. In the simplest explanation, we see only the tip of the iceberg in terms of star formation and gas available for future star formation, consistent with our remarkable finding that the majority of newly confirmed DRC galaxies are not the brightest continuum emitters in their immediate vicinity. Thus while ALMA continuum follow-up of SMGs identifies the brightest continuum emitters in each field, it does not necessarily reveal all the gas-rich galaxies. To hunt effectively for protocluster members requires wide and deep spectral-line imaging to uncover any relatively continuum-faint galaxies that are rich in atomic or molecular gas. Searching with short-baseline arrays or single-dish facilities, the true scale of the underlying gas reservoirs may be revealed.

Measurements of the cosmic far-infrared background (CIB) indicate that emission from many extragalactic phenomena, including star formation and black hole accretion, in the Universe can be obscured by dust. Resolving the CIB to study the population of galaxies in which this activity takes place is a major goal of submillimetre astronomy. Here, we present interferometric 650$\mu$m submillimetre number counts. Using the Band 8 data from the ALMACAL survey, we have analysed 81 ALMA calibrator fields together covering a total area of 5.5~arcmin$^2$. The typical central rms in these fields is $\sim 100 \mu$Jy~beam$^{-1}$ with the deepest maps reaching $\sigma = 47 \mu$Jy~beam$^{-1}$ at sub-arcsec resolution. Multi-wavelength coverage from ALMACAL allows us to exclude contamination from jets associated with the calibrators. However, residual contamination by jets and lensing remain a possibility. Using a signal-to-noise threshold of $4.5\sigma$, we find 21 dusty, star-forming galaxies with 650$\mu$m flux densities of $\geq 0.7 $mJy. At the detection limit we resolve $\simeq 100$ per cent of the CIB at 650$\mu$m, a significant improvement compared to low resolution studies at similar wavelength. We have therefore identified all the sources contributing to the EBL at 650 microns and predict that the contribution from objects with flux 0.7<mJy will be small.

We present an overview and description of the eMERLIN Galaxy Evolution survey (eMERGE) Data Release 1 (DR1), a large program of high-resolution 1.5 GHz radio observations of the GOODS-N field comprising $\sim140$ hours of observations with eMERLIN and $\sim40$ hours with the Very Large Array (VLA). We combine the long baselines of eMERLIN (providing high angular resolution) with the relatively closely-packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5 GHz radio survey with the sensitivity ($\sim 1.5\mu$Jy beam$^{-1}$), angular resolution ($0.2"$--$0.7"$) and field-of-view ($\sim15' \times 15'$) to detect and spatially resolve star-forming galaxies and AGN at $z\gtrsim 1$. The goal of eMERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over a $\sim20$ year baseline and the development of new point spread function/primary beam models to seamlessly merge eMERLIN and VLA data in the $uv$ plane. We present early science results, including measurements of the luminosities and/or linear sizes of $\sim500$ galaxes selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of $r_{\rm eMERGE}/r_{\rm HST}\sim1.02\pm0.03$, suggesting that in most high-redshift galaxies, the $\sim$GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers which will explore the $\sim$kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by eMERGE DR1.

With the spatial resolution of the Atacama Large Millimetre Array (ALMA), dusty galaxies in the distant Universe typically appear as single, compact blobs of dust emission, with a median half-light radius, $\approx$ 1 kpc. Occasionally, strong gravitational lensing by foreground galaxies or galaxy clusters has probed spatial scales 1-2 orders of magnitude smaller, often revealing late-stage mergers, sometimes with tantalising hints of sub-structure. One lensed galaxy in particular, the Cosmic Eyelash at $z=$ 2.3, has been cited extensively as an example of where the interstellar medium exhibits obvious, pronounced clumps, on a spatial scale of $\approx$ 100 pc. Seven orders of magnitude more luminous than giant molecular clouds in the local Universe, these features are presented as circumstantial evidence that the blue clumps observed in many $z\sim$ 2-3 galaxies are important sites of ongoing star formation, with significant masses of gas and stars. Here, we present data from ALMA which reveal that the dust continuum of the Cosmic Eyelash is in fact smooth and can be reproduced using two S\'ersic profiles with effective radii, 1.2 and 4.4 kpc, with no evidence of significant star-forming clumps down to a spatial scale of $\approx$ 80 pc and a star-formation rate of $<$ 3 M$_\odot$ yr$^{-1}$.

We present results of the MUSE-ALMA Halos, an ongoing study of the Circum-Galactic Medium (CGM) of low redshift galaxies (z < 1.4), currently comprising 14 strong HI absorbers in five quasar fields. We detect 43 galaxies associated with absorbers down to star formation rate (SFR) limits of 0.01-0.1 solar masses/yr, found within impact parameters (b) of 250 kpc from the quasar sightline. Excluding the targeted absorbers, we report a high detection rate of 89 per cent and find that most absorption systems are associated with pairs or groups of galaxies (three to eleven members). We note that galaxies with the smallest impact parameters are not necessarily the closest to the absorbing gas in velocity space. Using a multi-wavelength dataset (UVES/HIRES, HST, MUSE), we combine metal and HI column densities, allowing for derivation of the lower limits of neutral gas metallicity as well as emission line diagnostics (SFR, metallicities) of the ionised gas in the galaxies. We find that groups of associated galaxies follow the canonical relations of N(HI) -- b and W_r(2796) -- b, defining a region in parameter space below which no absorbers are detected. The metallicity of the ISM of associated galaxies, when measured, is higher than the metallicity limits of the absorber. In summary, our findings suggest that the physical properties of the CGM of complex group environments would benefit from associating the kinematics of individual absorbing components with each galaxy member.

We are just starting to understand the physical processes driving the dramatic change in cosmic star-formation rate between $z\sim 2$ and the present day. A quantity directly linked to star formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, thanks to the large redshift path surveyed ($\Delta z=182$), we provide constraints on the molecular column density distribution function beyond $z\sim 0$. In addition, we probe column densities of N(H$_2$) > 10$^{16}$ atoms~cm$^{-2}$, five orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3} \text{M}_{\odot} \text{Mpc}^{-3}$ at $0 < z \leq 1.7$ already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star-formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to $z\sim 2$.

The matched filter (MF) represents one of the main tools to detect signals from known sources embedded in the noise. In the Gaussian case the noise is assumed to be the realization of a Gaussian random field (GRF). The most important property of the MF, the maximization of the probability of detection subject to a constant probability of false detection or false alarm (PFA), makes it one of the most popular techniques. However, the MF technique relies upon the a priori knowledge of the number and the position of the searched signals in the GRF which usually are not available. A typical way out is to assume that the position of a signal coincides with one of the peaks in the matched filtered data. A detection is claimed when the probability that a given peak is due only to the noise (i.e. the PFA) is smaller than a prefixed threshold. In this case the probability density function (PDF) of the amplitudes has to be used for the computation of the PFA, which is different from the Gaussian. Moreover, the probability that a detection is false depends on the number of peaks present in the filtered GRF, the greater the number of peaks in a GRF, the higher the probability of peaks due to the noise that exceed the detection threshold. If not taken into account, the PFA can be severely underestimated. Many solutions proposed to this problem are non-parametric hence not able to exploit all the available information. This limitation has been overcome by means of two efficient parametric approaches, one based on the PDF of the peak amplitudes of a smooth and isotropic GRF whereas the other uses the Gumbel distribution (the asymptotic PDF of the corresponding extreme). Simulations and ALMA maps show that, although the two methods produce almost identical results, the first is more flexible and allows us to check the reliability of the detection procedure.

Galaxy halos appear to be missing a large fraction of their baryons, most probably hiding in the circumgalactic medium (CGM), a diffuse component within the dark matter halo that extends far from the inner regions of the galaxies. A powerful tool to study the CGM gas is offered by absorption lines in the spectra of background quasars. Here, we present optical (MUSE) and mm (ALMA) observations of the field of the quasar Q1130-1449 which includes a log [N(H I)/cm^-2]=21.71+/-0.07 absorber at z=0.313. Ground-based VLT/MUSE 3D spectroscopy shows 11 galaxies at the redshift of the absorber down to a limiting SFR>0.01 M_sun yr^-1 (covering emission lines of [OII], Hbeta, [OIII], [NII] and Halpha), 7 of which are new discoveries. In particular, we report a new emitter with a smaller impact parameter to the quasar line-of-sight (b=10.6 kpc) than the galaxies detected so far. Three of the objects are also detected in CO(1-0) in our ALMA observations indicating long depletion timescales for the molecular gas and kinematics consistent with the ionised gas. We infer from dedicated numerical cosmological RAMSES zoom-in simulations that the physical properties of these objects qualitatively resemble a small group environment, possibly part of a filamentary structure. Based on metallicity and velocity arguments, we conclude that the neutral gas traced in absorption is only partly related to these emitting galaxies while a larger fraction is likely the signature of gas with surface brightness almost four orders of magnitude fainter that current detection limits. Together, these findings challenge a picture where strong-HI quasar absorbers are associated with a single bright galaxy and favour a scenario where the HI gas probed in absorption is related to far more complex galaxy structures.

Gas-rich galaxies are selected efficiently via quasar absorption lines. Recently, a new perspective on such absorption-selected systems has opened up by studying the molecular gas content of absorber host galaxies using ALMA CO emission line observations. Here, we present an analysis of multiple CO transitions ($L'_{\rm CO} \sim 10^9$ K km s$^{-1}$) in two $z \sim 0.5$ galaxies associated with one Ly$\alpha$ absorber towards J0238+1636. The CO spectral line energy distribution (CO SLED) of these galaxies appear distinct from that of typical star-forming galaxies at similar redshifts and is comparable with that of luminous infrared galaxies or AGN. Indeed, these galaxies are associated with optically identified AGN activity. We infer that the CO line ratios and the $\alpha_{\rm CO}$ conversion factor differ from the Galactic values. Our findings suggest that at least a fraction of absorption selected systems shows ISM conditions deviating from those of normal star-forming galaxies. For a robust molecular gas mass calculation, it is therefore important to construct the CO SLED. Absorption-line-selection identifies systems with widely distributed gas, which may preferentially select interacting galaxies, which in turn will have more excited CO SLEDs than isolated galaxies. Furthermore, we raise the question whether quasar absorbers preferentially trace galaxy overdensities.

We have analysed archival VLA 8.4-GHz monitoring data of the gravitational lens system JVAS B1030+074 with the goal of determining the time delay between the two lensed images via the polarization variability. In contrast to the previously published total intensity variations, we detect correlated variability in polarized flux density, percentage polarization and polarization position angle. The latter includes a fast ($<$5d) 90-degree rotation event. Our best estimate of the time delay is $146\pm6$d (1$\sigma$), considerably longer than that predicted by the lens model presented in the discovery paper. Additional model constraints will be needed before this system can be used to measure $H_0$, for example through a detection of the lensed source's VLBI jet in image B. No time delay is visible in total flux density and this is partially due to much greater scatter in the image B measurements. This must be due to a propagation effect as the radio waves pass through the ISM of the lensing galaxies or the Galaxy.

We present deep high resolution (0.03", 200pc) ALMA Band 7 observations covering the dust continuum and [CII] $\lambda157.7\mu$m emission in four $z\sim4.4-4.8$ sub-millimeter galaxies (SMGs) selected from the ALESS and AS2UDS surveys. The data show that the rest-frame 160$\mu$m (observed 345 GHz) dust emission is consistent with smooth morphologies on kpc scales for three of the sources. One source, UDS47.0, displays apparent substructure but this is also consistent with a smooth morphology, as indicated by simulations showing that smooth exponential disks can appear clumpy when observed at high angular resolution (0.03") and depth of these observations ($\sigma_{345\text{GHz}} \sim27-47\mu$Jy beam$^{-1}$). The four SMGs are bright [CII] emitters, and we extract [CII] spectra from the high resolution data, and recover $\sim20-100$% of the [CII] flux and $\sim40-80$% of the dust continuum emission, compared to the previous lower resolution observations. When tapered to 0.2" resolution our maps recover $\sim80-100$% of the continuum emission, indicating that $\sim60$% of the emission is resolved out on $\sim200$pc scales. We find that the [CII] emission in high-redshift galaxies is more spatially extended than the rest-frame 160$\mu$m dust continuum by a factor of $1.6\pm0.4$. By considering the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of the star-formation rate surface density ($\Sigma_{\text{SFR}}$) we revisit the [CII] deficit, and suggest that the decline in the $L_{\text{[CII]}}$/$L_{\text{FIR}}$ ratio as a function of $\Sigma_{\text{SFR}}$ is consistent with local processes. We also explore the physical drivers that may be responsible for these trends and can give rise to the properties found in the densest regions of SMGs.

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

We present the results of a new study of the relationship between infrared excess (IRX), UV spectral slope (beta) and stellar mass at redshifts 2<z<3, based on a deep Atacama Large Millimeter Array (ALMA) 1.3-mm continuum mosaic of the Hubble Ultra Deep Field (HUDF). Excluding the most heavily-obscured sources, we use a stacking analysis to show that z~2.5 star-forming galaxies in the mass range 9.25 <= log(M/Msun) <= 10.75 are fully consistent with the IRX-beta relation expected for a relatively grey attenuation curve, similar to the commonly adopted Calzetti law. Based on a large, mass complete, sample of 2 <= z <= 3 star-forming galaxies drawn from multiple surveys, we proceed to derive a new empirical relationship between beta and stellar mass, making it possible to predict UV attenuation (A_1600) and IRX as a function of stellar mass, for any assumed attenuation law. Once again, we find that z~2.5 star-forming galaxies follow A_1600-mass and IRX-mass relations consistent with a relatively grey attenuation law, and find no compelling evidence that star-forming galaxies at this epoch follow a reddening law as steep as the Small Magellanic Cloud (SMC) extinction curve. In fact, we use a simple simulation to demonstrate that previous determinations of the IRX-beta relation may have been biased toward low values of IRX at red values of beta, mimicking the signature expected for an SMC-like dust law. We show that this provides a plausible mechanism for reconciling apparently contradictory results in the literature and that, based on typical measurement uncertainties, stellar mass provides a cleaner prediction of UV attenuation than beta. Although the situation at lower stellar masses remains uncertain, we conclude that for 2<z<3 star-forming galaxies with log(M/Msun) >= 9.75, both the IRX-beta and IRX-mass relations are well described by a Calzetti-like attenuation law.

Studying the flow of baryons into and out of galaxies is an important part of understanding the evolution of galaxies over time. We present a detailed case study of the environment around an intervening Ly $\alpha$ absorption line system at $z_{\rm abs} = 0.633$, seen towards the quasar J0423$-$0130 ($z_{\rm QSO} = 0.915$). We detect with ALMA the $^{12}$CO(2--1), $^{12}$CO(3--2) and $1.2$~mm continuum emission from a galaxy at the redshift of the Ly $\alpha$ absorber at a projected distance of $135$ kpc. From the ALMA detections, we infer ISM conditions similar to those in low redshift Luminous Infrared Galaxies. DDT MUSE integral field unit observations reveal the optical counterpart of the $^{12}$CO emission line source and three additional emission line galaxies at the absorber redshift, which together form a galaxy group. The $^{12}$CO emission line detections originate from the most massive galaxy in this group. While we cannot exclude that we miss a fainter host, we reach a dust-uncorrected star-formation rate (SFR) limit of > $0.3 \text{M}_{\odot} \text{ yr}^{-1}$ within $100$ kpc from the sightline to the background quasar. We measure the dust-corrected SFR (ranging from $3$ to $50$ M$_{\odot}$ yr$^{-1}$), the morpho-kinematics and the metallicities of the four group galaxies to understand the relation between the group and the neutral gas probed in absorption. We find that the Ly $\alpha$ absorber traces either an outflow from the most massive galaxy or intra-group gas. This case study illustrates the power of combining ALMA and MUSE to obtain a census of the cool baryons in a bounded structure at intermediate redshift.

We present the results of the first, deep ALMA imaging covering the full 4.5 sq arcmin of the Hubble Ultra Deep Field (HUDF) as previously imaged with WFC3/IR on HST. Using a mosaic of 45 pointings, we have obtained a homogeneous 1.3mm image of the HUDF, achieving an rms sensitivity of 35 microJy, at a resolution of 0.7 arcsec. From an initial list of ~50 >3.5sigma peaks, a rigorous analysis confirms 16 sources with flux densities S(1.3) > 120 microJy. All of these have secure galaxy counterparts with robust redshifts (<z> = 2.15), and 12 are also detected at 6GHz in new deep JVLA imaging. Due to the wealth of supporting data in this unique field, the physical properties of the ALMA sources are well constrained, including their stellar masses (M*) and UV+FIR star-formation rates (SFR). Our results show that stellar mass is the best predictor of SFR in the high-z Universe; indeed at z > 2 our ALMA sample contains 7 of the 9 galaxies in the HUDF with M* > 2 x 10^10 Msun and we detect only one galaxy at z > 3.5, reflecting the rapid drop-off of high-mass galaxies with increasing redshift. The detections, coupled with stacking, allow us to probe the redshift/mass distribution of the 1.3-mm background down to S(1.3) ~ 10 micro-Jy. We find strong evidence for a steep `main sequence' for star-forming galaxies at z ~ 2, with SFR \propto M* and a mean specific SFR = 2.2 /Gyr. Moreover, we find that ~85% of total star formation at z ~ 2 is enshrouded in dust, with ~65% of all star formation at this epoch occurring in high-mass galaxies (M* > 2 x 10^10 Msun), for which the average obscured:unobscured SF ratio is ~200. Finally, we combine our new ALMA results with the existing HST data to revisit the cosmic evolution of star-formation rate density; we find that this peaks at z ~ 2.5, and that the star-forming Universe transits from primarily unobscured to primarily obscured thereafter at z ~ 4.

We present $\simeq$0$.\!\!^{\prime\prime}4$-resolution extinction-independent distributions of star formation and dust in 11 star-forming galaxies (SFGs) at $z = 1.3-3.0$. These galaxies are selected from sensitive, blank-field surveys of the $2' \times 2'$ Hubble Ultra-Deep Field at $\lambda = 5$ cm and 1.3 mm using the Karl G. Jansky Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA). They have star-formation rates (SFRs), stellar masses, and dust properties representative of massive main-sequence SFGs at $z \sim 2$. Morphological classification performed on spatially-resolved stellar mass maps indicates a mixture of disk and morphologically disturbed systems; half of the sample harbor X-ray active galactic nuclei (AGN), thereby representing a diversity of $z \sim 2$ SFGs undergoing vigorous mass assembly. We find that their intense star formation most frequently occurs at the location of stellar-mass concentration and extends over an area comparable to their stellar-mass distribution, with a median diameter of $4.2 \pm 1.8$ kpc. This provides direct evidence for galaxy-wide star formation in distant, blank-field-selected main-sequence SFGs. The typical galactic-average SFR surface density is 2.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$, sufficiently high to drive outflows. In X-ray-selected AGN where radio emission is enhanced over the level associated with star formation, the radio excess pinpoints the AGN, which are found to be co-spatial with star formation. The median extinction-independent size of main-sequence SFGs is two times larger than those of bright submillimeter galaxies whose SFRs are $3-8$ times larger, providing a constraint on the characteristic SFR ($\sim300$ M$_{\odot}$yr$^{-1}$) above which a significant population of more compact star-forming galaxies appears to emerge.

We present high-resolution (0.16$"$) 870um Atacama Large Millimeter/submillimeter Array (ALMA) imaging of 16 luminous (L_IR ~ 4 x 10^12 L_sun) submillimeter galaxies (SMGs) from the ALESS survey of the Extended Chandra Deep Field South. This dust imaging traces the dust-obscured star formation in these z~2.5 galaxies on ~1.3 kpc scales. The emission has a median effective radius of $R_e=0.24" \pm 0.02"$, corresponding to a typical physical size of $R_{e}=1.8\pm$0.2 kpc. We derive a median S\'ersic index of $n=0.9\pm0.2$, implying that the dust emission is remarkably disk-like at the current resolution and sensitivity. We use different weighting schemes with the visibilities to search for clumps on 0.12$"$ (~1.0 kpc) scales, but we find no significant evidence for clumping in the majority of cases. Indeed, we demonstrate using simulations that the observed morphologies are generally consistent with smooth exponential disks, suggesting that caution should be exercised when identifying candidate clumps in even moderate S/N interferometric data. We compare our maps to comparable-resolution HST H$_{160}$-band images, finding that the stellar morphologies appear significantly more extended and disturbed, and suggesting that major mergers may be responsible for driving the formation of the compact dust disks we observe. The stark contrast between the obscured and unobscured morphologies may also have implications for SED fitting routines that assume the dust is co-located with the optical/near-IR continuum emission. Finally, we discuss the potential of the current bursts of star formation to transform the observed galaxy sizes and light profiles, showing that the $z\sim0$ descendants of these SMGs are expected to have stellar masses, effective radii, and gas surface densities consistent with the most compact massive (M_* ~ 1-2 x 10^11 M_sun) early-type galaxies observed locally.

We present ALMA ultra-high-spatial resolution ($\sim 20 \, {\rm mas}$) observations of dust continuum at $920 \, {\rm \mu m}$ and $1.2 \, {\rm mm}$ in a pair of submm galaxies (SMGs) at $z = 3.442$, ALMACAL-1 (A-1: $S_{\rm 870 \mu m} = 6.5 \pm 0.2 \, {\rm mJy}$) and ALMACAL-2 (A-2: $S_{\rm 870 \mu m} = 4.4 \pm 0.2 \, {\rm mJy}$). The spectroscopic redshifts of A-1 and A-2 have been confirmed via serendipitous detection of up to nine emission lines. Our ultra-high-spatial resolution data reveal that about half of the star formation in each of these starbursts is dominated by a single compact clump (FWHM size of $\sim 350 \, {\rm pc}$). This structure is confirmed by independent datasets at $920 \, {\rm \mu m}$ and $1.2 \, {\rm mm}$. The star-formation rate (SFR) surface densities of all these clumps are extremely high, $\Sigma_{\rm SFR} \sim 1200$ to $\sim 3000 \, {M_\odot \, {\rm yr}^{-1} \, {\rm kpc}^{-2}}$, the highest found in high-redshift galaxies. There is a small probability that A-1 and A-2 are the lensed components of a background source gravitationally amplified by the blazar host. If this was the case, the effective radius of the source would be $R_{\rm eff} \sim 40 \, {\rm pc}$, and the de-magnified SFR surface density would be $\Sigma_{\rm SFR} \sim 10000 \, {M_\odot \, {\rm yr}^{-1} \, {\rm kpc}^{-2}}$, comparable with the eastern nucleus of Arp 220. Despite being unable to rule out an AGN contribution, our results suggest that a significant percentage of the enormous far-IR luminosity in some dusty starbursts is concentrated in very small star-forming regions. The high $\Sigma_{\rm SFR}$ in our pair of SMGs could only be measured thanks to the ultra-high-resolution ALMA observations used in this work, demonstrating that long-baseline observations are essential to study and interpret the properties of dusty starbursts in the early Universe.

We present global VLBI observations of the 21-cm transition of atomic hydrogen seen in absorption against the radio source J0855+5751. The foreground absorber (SDSS~J085519.05+575140.7) is a dwarf galaxy at $z$ = 0.026. As the background source is heavily resolved by VLBI, the data allow us to map the properties of the foreground HI gas with a spatial resolution of 2pc. The absorbing gas corresponds to a single coherent structure with an extent $>$35pc, but we also detect significant and coherent variations, including a change in the HI optical depth by a factor of five across a distance of $\leq$6pc. The large size of the structure provides support for the Heiles & Troland model of the ISM, as well as its applicability to external galaxies. The large variations in HI optical depth also suggest that caution should be applied when interpreting $T_S$ measurements from radio-detected DLAs. In addition, the distorted appearance of the background radio source is indicative of a strong jet-cloud interaction in its host galaxy. We have measured its redshift ($z$ = 0.54186) using optical spectroscopy on the William Herschel Telescope and this confirms that J0855+5751 is a FRII radio source with a physical extent of $<$1kpc and supports the previous identification of this source as a Compact Symmetric Object. These sources often show absorption associated with the host galaxy and we suggest that both HI and OH should be searched for in J0855+5751.

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 have exploited ALMA calibration observations to carry out a novel, wide and deep submm survey, ALMACAL. These calibration data comprise a large number of observations of calibrator fields in a variety of frequency bands and array configurations. Gathering together data acquired during multiple visits to many ALMA calibrators, it is possible to reach noise levels which allow the detection of faint dusty, star-forming galaxies (DSFGs) over a significant area. In this paper we outline our survey strategy and report the first results. We have analysed data for 69 calibrators, reaching depths of $\sim 25 \, {\rm \mu Jy \, beam^{-1}}$ at sub-arcsec resolution. Adopting a conservative approach based on $\geq 5 \sigma$ detections, we have found eight and 11 DSFGs in ALMA bands 6 and 7, respectively, with flux densities $S_{\rm 1.2 mm} \geq 0.2 \, {\rm mJy}$. The faintest galaxies would have been missed by even the deepest \emph{Herschel} surveys. Our cumulative number counts have been determined independently at 870 $\mu$m and 1.2 mm, from a sparse sampling of the astronomical sky, and are thus relatively free of cosmic variance. The counts are lower than reported previously by a factor of at least $2\times$. Future analyses will yield large, secure samples of DSFGs, with redshifts determined via detection of submm spectral lines. Uniquely, our strategy then allows morphological studies of very faint DSFGs - representative of more normal star-forming galaxies than conventional submm galaxies (SMGs) - in fields where self-calibration is feasible, yielding milliarcsecond spatial resolution.

We present the results of a survey for intervening HI 21-cm absorbers at intermediate and low redshift (0<z<1.2). For our total sample of 24 systems, we obtained high quality data for 17 systems, the other seven being severely affected by radio frequency interference (RFI). Five of our targets are low redshift (z<0.17) optical galaxies with small impact parameters (<20 kpc) toward radio-bright background sources. Two of these were detected in 21-cm absorption, showing narrow, high optical depth absorption profiles, the narrowest having a velocity dispersion of only 1.5 km/s, which puts an upper limit on the kinetic temperature of T_k<270 K. Combining our observations with results from the literature, we measure a weak anti-correlation between impact parameter and integral optical depth in local (z<0.5) 21-cm absorbers. Of eleven CaII and MgII systems searched, two were detected in 21-cm absorption, and six were affected by RFI to a level that precludes a detection. For these two systems at z~0.6 we measure spin temperatures of T_s=(65+/-17) K and T_s>180 K. A subset of our systems were also searched for OH absorption, but no detections were made.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Adding VLBI capability to the SKA arrays will greatly broaden the science of the SKA, and is feasible within the current specifications. SKA-VLBI can be initially implemented by providing phased-array outputs for SKA1-MID and SKA1-SUR and using these extremely sensitive stations with other radio telescopes, and in SKA2 by realising a distributed configuration providing baselines up to thousands of km, merging it with existing VLBI networks. The motivation for and the possible realization of SKA-VLBI is described in this paper.

We analyse HST WFC3/$H_{160}$-band observations of a sample of 48 ALMA-detected submillimeter galaxies (SMGs) in the Extended Chandra Deep Field South field, to study their stellar morphologies and sizes. We detect 79$\pm$17% of the SMGs in the $H_{160}$-band imaging with a median sensitivity of 27.8 mag, and most (80%) of the non-detections are SMGs with 870$\mu$m fluxes of $S_{870} < $3 mJy. With a surface brightness limit of $\mu_H \sim $26 mag arcsec$^{-2}$, we find that 82$\pm$9% of the $H_{160}$-band detected SMGs at $z =$ 1-3 appear to have disturbed morphologies, meaning they are visually classified as either irregulars or interacting systems, or both. By determining a S\'ersic fit to the $H_{160}$ surface-brightness profiles we derive a median S\'ersic index of $n = $1.2$\pm$0.3 and a median half-light radius of $r_e = $4.4$^{+1.1}_{-0.5}$ kpc for our SMGs at $z = $1-3. We also find significant displacements between the positions of the $H_{160}$-component and 870$\mu$m emission in these systems, suggesting that the dusty star-burst regions and less-obscured stellar distribution are not co-located. We find significant differences in the sizes and the S\'ersic index between our $z = $2-3 SMGs and $z \sim $2 quiescent galaxies, suggesting a major transformation of the stellar light profile is needed in the quenching processes if SMGs are progenitors of the red-and-dead $z\sim$2 galaxies. Given the short-lived nature of SMGs, we postulate that the majority of the $z = $2-3 SMGs with $S_{870} \gtrsim $2 mJy are early/mid-stage major mergers.

We performed a kinematical analysis of the [CII] line emission of the BR 1202-0725 system at z~4,7 using ALMA observations. The most prominent sources of this system are a quasar and a submillimeter galaxy, separated by a projected distance of about 24 kpc and characterized by very high SFR, higher than 1000 Msun/yr. However, the ALMA observations reveal that these galaxies apparently have undisturbed rotating disks, which is at variance with the commonly accepted scenario in which strong star formation activity is induced by a major merger. We also detected faint components which, after spectral deblending, were spatially resolved from the main QSO and SMG emissions. The relative velocities and positions of these components are compatible with orbital motions within the gravitational potentials generated by the QSO host galaxy and the SMG, suggesting that they are smaller galaxies in interaction or gas clouds in accretion flows of tidal streams. We did not find any clear spectral evidence for outflows caused by AGN or stellar feedback. This suggests that the high star formation rates might be induced by interactions or minor mergers with these companions, which do not affect the large-scale kinematics of the disks, however. Our kinematical analysis also indicates that the QSO and the SMG have similar Mdyn, mostly in the form of molecular gas, and that the QSO host galaxy and the SMG are seen close to face-on with slightly different disk inclinations: the QSO host galaxy is seen almost face-on (i~15), while the SMG is seen at higher inclinations (i~25). Finally, the ratio between the black hole mass of the QSO, obtained from XShooter spectroscopy, and the Mdyn of the host galaxy is similar to value found in very massive local galaxies, suggesting that the evolution of black hole galaxy relations is probably better studied with dynamical than with stellar host galaxy masses.

We present an Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 survey of 126 submillimeter sources from the LABOCA ECDFS Submillimeter Survey (LESS). Our 870 micron survey with ALMA (ALESS) has produced maps ~3X deeper and with a beam area ~200X smaller than the original LESS observations, doubling the current number of interferometrically-observed submillimeter sources. The high resolution of these maps allows us to resolve sources that were previously blended and accurately identify the origin of the submillimeter emission. We discuss the creation of the ALESS submillimeter galaxy (SMG) catalog, including the main sample of 99 SMGs and a supplementary sample of 32 SMGs. We find that at least 35% (possibly up to 50%) of the detected LABOCA sources have been resolved into multiple SMGs, and that the average number of SMGs per LESS source increases with LESS flux density. Using the (now precisely known) SMG positions, we empirically test the theoretical expectation for the uncertainty in the single-dish source positions. We also compare our catalog to the previously predicted radio/mid-infrared counterparts, finding that 45% of the ALESS SMGs were missed by this method. Our ~1.6" resolution allows us to measure a size of ~9 kpc X 5 kpc for the rest-frame ~300 um emission region in one resolved SMG, implying a star formation rate surface density of 80 M_sol yr^-1 kpc^-2, and we constrain the emission regions in the remaining SMGs to be <10 kpc. As the first statistically reliable survey of SMGs, this will provide the basis for an unbiased multiwavelength study of SMG properties.

In the past decade, our understanding of galaxy evolution has been revolutionized by the discovery that luminous, dusty, starburst galaxies were 1,000 times more abundant in the early Universe than at present. It has, however, been difficult to measure the complete redshift 2 distribution of these objects, especially at the highest redshifts (z > 4). Here we report a redshift survey at a wavelength of three millimeters, targeting carbon monoxide line emission from the star-forming molecular gas in the direction of extraordinarily bright millimetrewave-selected sources. High-resolution imaging demonstrates that these sources are strongly gravitationally lensed by foreground galaxies. We detect spectral lines in 23 out of 26 sources and multiple lines in 12 of those 23 sources, from which we obtain robust, unambiguous redshifts. At least 10 of the sources are found to lie at z > 4, indicating that the fraction of dusty starburst galaxies at high redshifts is greater than previously thought. Models of lens geometries in the sample indicate that the background objects are ultra-luminous infrared galaxies, powered by extreme bursts of star formation.

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we have conducted a blind redshift survey in the 3 mm atmospheric transmission window for 26 strongly lensd dusty star-forming galaxies (DSFGs) selected with the South Pole Telescope (SPT). The sources were selected to have S_1.4mm>20 mJy and a dust-like spectrum and, to remove low-z sources, not have bright radio (S_843MHz<6mJy) or far-infrared counterparts (S_100um<1 Jy, S_60um<200mJy). We robustly detect 44 line features in our survey, which we identify as redshifted emission lines of 12CO, 13CO, [CI], H2O, and H2O+. We find one or more spectral features in 23 sources yielding a ~90% detection rate for this survey; in 12 of these sources we detect multiple lines, while in 11 sources we detect only a single line. For the sources with only one detected line, we break the redshift degeneracy with additional spectroscopic observations if available, or infer the most likely line identification based on photometric data. This yields secure redshifts for ~70% of the sample. The three sources with no lines detected are tentatively placed in the redshift desert between 1.7<z<2.0. The resulting mean redshift of our sample is <z>=3.5. This finding is in contrast to the redshift distribution of radio-identified DSFGs, which have a significantly lower mean redshift of <z>=2.3 and for which only 10-15% of the population is expected to be at z>3. We discuss the effect of gravitational lensing on the redshift distribution and compare our measured redshift distribution to that of models in the literature.

We report the first counts of faint submillimetre galaxies (SMG) in the 870-um band derived from arcsecond resolution observations with the Atacama Large Millimeter Array (ALMA). We have used ALMA to map a sample of 122 870-um-selected submillimetre sources drawn from the (0.5x0.5)deg^2 LABOCA Extended Chandra Deep Field South Submillimetre Survey (LESS). These ALMA maps have an average depth of sigma(870um)~0.4mJy, some ~3x deeper than the original LABOCA survey and critically the angular resolution is more than an order of magnitude higher, FWHM of ~1.5" compared to ~19" for the LABOCA discovery map. This combination of sensitivity and resolution allows us to precisely pin-point the SMGs contributing to the submillimetre sources from the LABOCA map, free from the effects of confusion. We show that our ALMA-derived SMG counts broadly agree with the submillimetre source counts from previous, lower-resolution single-dish surveys, demonstrating that the bulk of the submillimetre sources are not caused by blending of unresolved SMGs. The difficulty which well-constrained theoretical models have in reproducing the high-surface densities of SMGs, thus remains. However, our observations do show that all of the very brightest sources in the LESS sample, S(870um)>12mJy, comprise emission from multiple, fainter SMGs, each with 870-um fluxes of <9mJy. This implies a natural limit to the star-formation rate in SMGs of <10^3 M_Sun/yr, which in turn suggests that the space densities of z>1 galaxies with gas masses in excess of ~5x10^10 M_Sun is <10^-5 Mpc^-3. We also discuss the influence of this blending on the identification and characterisation of the SMG counterparts to these bright submillimetre sources and suggest that it may be responsible for previous claims that they lie at higher redshifts than fainter SMGs.

We present ALMA 870-um (345GHz) observations of two sub-millimetre galaxies (SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data identify the counterparts to these previously unidentified sub-millimeter sources and serendipitously detect bright emission lines in their spectra which we show are most likely to be [C II]157.74um emission yielding redshifts of z=4.42 and z=4.44. This blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the previously derived photometric redshift distribution of SMGs and suggests a modest, but not dominant (<25%), tail of 870-um selected SMGs at z>4. We find that the ratio of L_CII/L_FIR in these SMGs is much higher than seen for similarly far-infrared-luminous galaxies at z~0, which is attributed to the more extended gas reservoirs in these high-redshift ULIRGs. Indeed, in one system we show that the [C II] emission shows hints of extended emission on >3kpc scales. Finally, we use the volume probed by our ALMA survey to show that the bright end of the [C II] luminosity function evolves strongly between z=0 and z~4.4, reflecting the increased ISM cooling in galaxies as a result of their higher star-formation rates. These observations demonstrate that even with short integrations, ALMA is able to detect the dominant fine structure cooling lines from high-redshift ULIRGs, measure their energetics and trace their evolution with redshift.

We present Herschel-PACS spectroscopy of the [OI]63um far-infrared cooling line from a sample of six unlensed and spectroscopically-confirmed 870um-selected submillimetre (submm) galaxies (SMGs) at 1.1<z<1.6 from the LABOCA Extended Chandra Deep Field South (ECDFS) Submm Survey (LESS). This is the first survey of [OI]63um, one of the main photodissociation region (PDR) cooling lines, in SMGs. New high-resolution ALMA interferometric 870um continuum imaging confirms that these six Herschel-targeted SMG counterparts are bona fide sources of submm emission. We detect [OI]63um in two SMGs with a SNR >3, tentatively detect [OI]63um in one SMG, and constrain the line flux for the non-detections. We also exploit the combination of submm continuum photometry from 250-870um and our new PACS continuum measurements to constrain the far-infrared (FIR) luminosity, L_FIR, in these SMGs to < 30%. We find that SMGs do not show a deficit in their [OI]63um-to-far-infrared continuum luminosity ratios (with ratios ranging from ~0.5-1.5%), similar to what was seen previously for the [CII]158um-to-FIR ratios in SMGs. These observed ratios are about an order of magnitude higher than what is seen typically for local ultra luminous infrared galaxies (ULIRGs), which adds to the growing body of evidence that SMGs are not simply `scaled up' versions of local ULIRGs. Rather, the PDR line-to-L_FIR ratios suggest that the star formation modes of SMGs are likely more akin to that of local normal (lower-luminosity) star-forming galaxies, with the bulk of the star formation occurring in extended regions, galaxy-scale (~kpc) in size. These observations represent the first step towards a census of the major PDR cooling lines in typical SMGs that will be attainable with ALMA, enabling detailed modelling to probe the global properties of the star formation and the evolutionary status of SMGs.

We present radio and infrared (3.6-24um) counterparts to submillimetre galaxies (SMGs) detected in the Extended Chandra Deep Field South with the LABOCA 870-um bolometer camera on the 12-m Atacama Pathfinder Experiment. Using the Very Large Array at 1.4GHz and Spitzer we have identified secure counterparts to 79 of the 126 SMGs (SNR>3.7, S_870>4.4mJy) in the field, 62 via their radio and/or 24-um emission, the remainder using a colour-flux cut on IRAC 3.6- and 5.8-um sources chosen to maximise the number of secure, coincident radio and 24-um counterparts. In constructing our radio catalogue, we have corrected for the effects of `flux boosting', then used the corrected flux densities to estimate the redshifts of the SMGs based on the radio/submm spectral indices. The effect of the boosting correction is to increase the median redshift by 0.2 resulting in a value of z=2.2 (+0.7-0.8) (1-sigma errors) for the secure radio counterparts, in agreement with other studies, both spectroscopic and photometric.

The periodogram is a popular tool that tests whether a signal consists only of noise or if it also includes other components. The main issue of this method is to define a critical detection threshold that allows identification of a component other than noise, when a peak in the periodogram exceeds it. In the case of signals sampled on a regular time grid, determination of such a threshold is relatively simple. When the sampling is uneven, however, things are more complicated. The most popular solution in this case is to use the "Lomb-Scargle" periodogram, but this method can be used only when the noise is the realization of a zero-mean, white (i.e. flat-spectrum) random process. In this paper, we present a general formalism based on matrix algebra, which permits analysis of the statistical properties of a periodogram independently of the characteristics of noise (e.g. colored and/or non-stationary), as well as the characteristics of sampling.

[abridged] We derive photometric redshifts from 17-band optical to mid-IR photometry of 74 robust counterparts to 68 of the 126 submillimetre galaxies (SMGs) selected at 870um by LABOCA observations in the ECDFS. The median photometric redshift of identified SMGs is z=2.2\pm0.1, the interquartile range is z=1.8-2.7 and we identify 10 (~15%) high-redshift (z>3) SMGs. We derive a simple redshift estimator for SMGs based on the 3.6 and 8um fluxes, which is accurate to Delta_z~0.4 for SMGs at z<4. A statistical analysis of sources around unidentified SMGs identifies a population of likely counterparts with a redshift distribution peaking at z=2.5\pm0.3, which likely comprises ~60% of the unidentified SMGs. This confirms that the bulk of the undetected SMGs are coeval with those detected in the radio/mid-IR. We conclude that ~15% of all the SMGs are below the flux limits of our survey and lie at z>3 and hence ~30% of all SMGs have z>3. We estimate that the full S_870um>4mJy SMG population has a median redshift of 2.5\pm0.6. In contrast to previous suggestions we find no significant correlation between S_870um and redshift. The median stellar mass of the SMGs derived from SED fitting is (9.2\pm0.9)x10^10Msun and the interquartile range is (4.7-14)x10^10Msun, although we caution that uncertainty in the star-formation histories results in a factor of ~5 uncertainty in these stellar masses. The median characteristic dust temperature of SMGs is 35.9\pm1.4K and the interquartile range is 28.5-43.3K. The infrared luminosity function shows that SMGs at z=2-3 typically have higher far-IR luminosities and luminosity density than those at z=1-2. This is mirrored in the evolution of the star-formation rate density (SFRD) for SMGs which peaks at z~2. The maximum contribution of bright SMGs to the global SFRD (~5% for SMGs with S_870um>4mJy; ~50% for SMGs with S_870um>1mJy) also occurs at z~2.

We present continent-scale VLBI - obtained with the European VLBI Network (EVN) at a wavelength of 18cm - of six distant, luminous submm-selected galaxies (SMGs). Our images have a synthesized beam width of ~30 milliarcsec FWHM - three orders of magnitude smaller in area than the highest resolution VLA imaging at this wavelength - and are capable of separating radio emission from ultra-compact radio cores (associated with active super-massive black holes - SMBHs) from that due to starburst activity. Despite targeting compact sources - as judged by earlier observations with the VLA and MERLIN - we identify ultra-compact cores in only two of our targets. This suggests that the radio emission from SMGs is produced primarily on larger scales than those probed by the EVN, and therefore is generated by star formation rather than an AGN - a result consistent with other methods used to identify the presence of SMBHs in these systems.

We present results from the earliest observations of DEBRIS, a Herschel Key Programme to conduct a volume- and flux-limited survey for debris discs in A-type through M-type stars. PACS images (from chop/nod or scan-mode observations) at 100 and 160 micron are presented toward two A-type stars and one F-type star: beta Leo, beta UMa and eta Corvi. All three stars are known disc hosts. Herschel spatially resolves the dust emission around all three stars (marginally, in the case of beta UMa), providing new information about discs as close as 11 pc with sizes comparable to that of the Solar System. We have combined these data with existing flux density measurements of the discs to refine the SEDs and derive estimates of the fractional luminosities, temperatures and radii of the discs.

We present high spatial resolution (0.4", ~3.5 kpc) PdBI interferometric data on three ultra-luminous infrared galaxies (ULIRGs) at z~2: two submillimetre galaxies and one submillimetre faint star forming radio galaxy. The three galaxies have been ro- bustly detected in CO rotational transitions, either 12CO(J=4-3) or 12CO(J=3-2), allowing their sizes and gas masses to be accurately constrained. These are the highest spatial resolution observations observed to date (by a factor of ~2) for intermediate-excitation CO emission in z~2 ULIRGs. The galaxies appear extended over several resolution elements, having a mean radius of 3.7 kpc. High-resolution (0.3") combined MERLIN-VLA observations of their radio continua allow an analysis of the star formation behaviour of these galaxies, on comparable spatial scales to that of the CO observations. This 'matched beam' approach sheds light on the spatial distribution of both molecular gas and star formation, and we can therefore calculate accurate star formation rates and gas surface densities: this allows us to place the three systems in the context of a Kennicutt-Schmidt (KS)-style star formation law. We find a difference in size between the CO and radio emission regions, and as such we suggest that using the spatial extent of the CO emission region to estimate the surface density of star formation may lead to error. This size difference also causes the star formation efficiencies within systems to vary by up to a factor of 5. We also find, with our new accurate sizes, that SMGs lie significantly above the KS relation, indicating that stars are formed more efficiently in these extreme systems than in other high-z star forming galaxies.

We present an 880 micron Submillimeter Array (SMA) detection of the submillimeter galaxy SXDF850.6. SXDF850.6 is a bright source (S(850 micron) = 8 mJy) detected in the SCUBA Half Degree Extragalactic Survey (SHADES), and has multiple possible radio counterparts in its deep radio image obtained at the VLA. Our new SMA detection finds that the submm emission coincides with the brightest radio emission that is found ~8" north of the coordinates determined from SCUBA. Despite the lack of detectable counterparts in deep UV/optical images, we find a source at the SMA position in near-infrared and longer wavelength images. We perform SED model fits to UV-optical-IR photometry (u, B, V, R, i', z', J, H, K, 3.6 micron, 4.5 micron, 5.8 micron, and 8.0 micron) and to submm-radio photometry (850 micron, 880 micron, 1100 micron, and 21 cm) independently, and we find both are well described by starburst templates at a redshift of z ~= 2.2 (+/- 0.3). The best-fit parameters from the UV-optical-IR SED fit are a redshift of z = 1.87 (+0.15/-0.07), a stellar mass of M_star = 2.5 +2.2/-0.3 x 10^11 M_sun, an extinction of A_V = 3.0 (+0.3/-1.0) mag, and an age of 720 (+1880/-210) Myr. The submm-radio SED fit provides a consistent redshift of z ~ 1.8-2.5, an IR luminosity of L_IR = (7-26) x 10^12 L_sun, and a star formation rate of 1300-4500 M_sun/yr. These results suggest that SXDF850.6 is a mature system already having a massive amount of old stellar population constructed before its submm bright phase and is experiencing a dusty starburst, possibly induced by major mergers.

We investigate the correlation between far-infrared (FIR) and radio luminosities in distant galaxies, a lynchpin of modern astronomy. We use data from the Balloon-borne Large Aperture Submillimetre Telescope (BLAST), Spitzer, the Large Apex BOlometer CamerA (LABOCA), the Very Large Array (VLA) and the Giant Metre-wave Radio Telescope (GMRT) in the Extended Chandra Deep Field South (ECDFS). For a catalogue of BLAST 250-micron-selected galaxies, we re-measure the 70--870-micron flux densities at the positions of their most likely 24-micron counterparts, which have a median [interquartile] redshift of 0.74 [0.25, 1.57]. From these, we determine the monochromatic flux density ratio, q_250 = log_10 (S_250micron / S_1400MHz), and the bolometric equivalent, q_IR. At z ~= 0.6, where our 250-micron filter probes rest-frame 160-micron emission, we find no evolution relative to q_160 for local galaxies. We also stack the FIR and submm images at the positions of 24-micron- and radio-selected galaxies. The difference between q_IR seen for 250-micron- and radio-selected galaxies suggests star formation provides most of the IR luminosity in ~< 100-uJy radio galaxies, but rather less for those in the mJy regime. For the 24-micron sample, the radio spectral index is constant across 0 < z < 3, but q_IR exhibits tentative evidence of a steady decline such that q_IR is proportional to (1+z)^(-0.15 +/- 0.03) - significant evolution, spanning the epoch of galaxy formation, with major implications for techniques that rely on the FIR/radio correlation. We compare with model predictions and speculate that we may be seeing the increase in radio activity that gives rise to the radio background.

A complete sample of 33 sources believed to be starbursts ("5.8micron-peakers") was selected in the (0.5 sq. deg.) J1046+56 field with selection criteria F_(24micron)>400muJy, the presence of a redshifted stellar emission peak at 5.8um, and r'(Vega)>23. The field, part of the SWIRE Lockman Hole field, benefits from very deep VLA/GMRT 20cm, 50cm, and 90cm radio data (all 33 sources are detected at 50cm), and deep 160micron and 70micron Spitzer data. The 33 sources, with photometric redshifts ~1.5-2.5, were observed at 1.2mm with IRAM-30m/MAMBO to an rms ~0.7-0.8mJy in most cases. Their millimeter, radio, 7-band Spitzer, and near-IR properties were jointly analyzed. The entire sample of 33 sources has an average 1.2mm flux density of 1.56+/-0.22mJy and a median of 1.61mJy, so the majority of the sources can be considered SMGs. Four sources have confirmed 4sigma detections, and nine were tentatively detected at the 3sigma level. Because of its 24micron selection, our sample shows systematically lower F_(1.2mm)/F_(24micron) flux ratios than classical SMGs, probably because of enhanced PAH emission. A median FIR SED was built by stacking images at the positions of 21 sources in the region of deepest Spitzer coverage. Its parameters are T_(dust)=37+/-8K, L_(FIR)=2.5x10^{12}Lo, and SFR=450Mo/yr. The FIR-radio correlation provides another estimate of L_(FIR) for each source, with an average value of 4.1x10^{12}Lo; however, this value may be overestimated because of some AGN contribution. Most of our targets are also luminous star-forming BzK galaxies which constitute a significant fraction of weak SMGs at 1.7<z<2.3.

We identify eight z>1 radio sources undetected at 850um but robustly detected at 70um, confirming that they represent ultraluminous infrared galaxies (ULIRGs) with hotter dust temperatures (T_dust=52+-10 K) than Submillimetre Galaxies (SMGs) at similar luminosities and redshifts. These galaxies share many properties with SMGs: ultra-violet (UV) spectra consistent with starbursts, high stellar masses and radio luminosities. We can attribute their radio emission to star formation since high-resolution MERLIN radio maps show extended emission regions (with half light radii of 2-3kpc), which are unlikely to be generated by AGN activity. These observations provide the first direct confirmation of hot, dusty ULIRGs which are missed by current submillimetre surveys. They have significant implications for future observations from the Herschel Space Observatory and SCUBA2, which will select high redshift luminous galaxies with less selection biases.

We describe multi-frequency data designed to probe the emission mechanism that dominates in faint radio sources. Our analysis is based on observations of the Lockman Hole using the GMRT - the deepest 610-MHz imaging yet reported - together with 1.4-GHz imaging from the VLA, well matched in resolution and sensitivity to the GMRT data. The GMRT and VLA data are cross-matched to obtain the radio spectral indices for the faint radio emitters. Statistical analyses show no clear evolution for the median spectral index, alpha, as a function of flux density. alpha is found to be -0.6 to -0.7, based on an almost unbiased 10-sigma criterion, down to a flux level of S(1.4GHz) >~ 100uJy. The fraction of inverted spectrum sources (alpha > 0) is less than 10%. The results suggest that the most prevalent emission mechanism in the sub-mJy regime is optically-thin synchrotron, ruling out a dominant flat spectrum or ultra-steep spectrum radio population. The spectral index distribution has a significant scatter, 0.4-0.5, which suggests a mixture of different populations at all flux levels. Spectroscopic classification of radio sources with X-ray emission has allowed us to estimate that the fraction of radio-quiet AGN at 30 <~ S(1.4GHz) < 300uJy is roughly 25 +/- 10%, suggesting that star-forming galaxies dominate the sub-mJy regime.

High-resolution submm imaging of the HzRG, 4C60.07, at z=3.8, has revealed two dusty components. Spitzer imaging shows that one of these components (B) is coincident with an extremely red AGN, offset by ~4" (~30 kpc) from the HzRG core. The other submm component (A) - resolved by our beam and devoid of emission at 3.6-8.0um - lies between B and the HzRG core. Since the HzRG was discovered via its young, steep-spectrum lobes and their creation was likely triggered by the interaction, we argue that we are witnessing an early-stage merger, prior to its eventual equilibrium state. The interaction is between the host galaxy of an actively-fueled BH, and a gas-rich starburst/AGN (B) marked by the compact submm component and coincident with broad CO emission. `A' is a plume of cold, dusty gas, associated with a narrow (~150 km/s) CO feature, and may represent a short-lived tidal structure. It has been claimed that HzRGs and SMGs differ only in the activity of their AGNs, but such complex submm morphologies are seen only rarely amongst SMGs. Our study has important implications: where a galaxy's gas is not aligned with its central BH, CO may be an unreliable probe of dynamical mass, affecting work on the co-assembly of BHs and spheroids. Our data support the picture wherein close binary AGN are induced by mergers. They also raise the possibility that some supposedly jet-induced starbursts may have formed co-evally with (yet independently of) the radio jets, both triggered by the same interaction. We note that the HzRG host would have gone unnoticed without its jets/companion, so there may be many other unseen BHs at high redshift, lost in the sea of ~5 x 10^8 similarly bright IRAC sources - sufficiently massive to drive a >10^27-W radio source, yet practically invisible unless actively fueled (abridged).

We present high-resolution interferometric imaging of LH850.02, the brightest 850- and 1200-micron submillimetre (submm) galaxy in the Lockman Hole. Our observations were made at 890 micron with the Submillimetre Array (SMA). Our high-resolution submm imaging detects LH850.02 at >6-sigma as a single compact (size < 1 arcsec or < 8 kpc) point source and yields its absolute position to ~0.2-arcsec accuracy. LH850.02 has two alternative radio counterparts within the SCUBA beam (LH850.02N & S), both of which are statistically very unlikely to be so close to the SCUBA source position by chance. However, the precise astrometry from the SMA shows that the submm emission arises entirely from LH850.02N, and is not associated with LH850.02S (by far the brighter of the two alternative identifications at 24-micron). Fits to the optical-infrared multi-colour photometry of LH850.02N & S indicate that both lie at z~3.3, and are therefore likely to be physically associated. At these redshifts, the 24 micron--to--submm flux density ratios suggest that LH850.02N has an Arp220-type starburst-dominated far-IR SED, while LH850.02S is more similar to Mrk231, with less dust-enshrouded star-formation activity, but a significant contribution at 24-micron (rest-frame ~5-6 micron) from an active nucleus. This complex mix of star-formation and AGN activity in multi-component sources may be common in the high redshift ultraluminous galaxy population, and highlights the need for precise astrometry from high resolution interferometric imaging for a more complete understanding.

We have analysed the 24um properties of a radio-selected sample in the Subaru-XMM/Newton Deep Field in order to explore the behaviour of the FIR/radio relation at high redshifts. Statistically, the correlation is described by q24, the ratio between the observed flux densities at 24um and 1.4GHz, respectively. Using 24um data results in considerably more scatter in the correlation than previous work using data at 60-70um. Nevertheless, we do observe a steady correlation as a function of redshift, up to z~3.5, suggesting its validity back to primeval times. We find q24 = 0.30 +/- 0.56 for the observed and q24 = 0.71 +/- 0.47 for the k-corrected radio sample, based on sources with 300uJy < S(1.4GHz) < 3.2mJy and 24um detections. A suitable k-correction given by a M82-like mid-IR template suggests no extreme silicate absorption in the bulk of our radio sample. Using thresholds in q24 to identify radio-excess sources, we have been able to characterise the transition from radio-loud AGN to star-forming galaxies and radio-quiet AGN at faint (<1mJy) radio flux densities. Our results are in broad agreement with previous studies which show a dominant radio-loud AGN population at >1mJy. The rest-frame U-B colours of the expected radio-excess population have redder distribution than those that follow the correlation. This is therefore a promising way to select obscured Type-2 AGN, with a radio loud nature, missed by deep X-ray observations. Spectroscopic follow-up of these sources is required to fully test this method.

We have produced sensitive, high-resolution radio maps of 12 SMGs in the Lockman Hole using combined MERLIN and VLA data at a frequency of 1.4 GHz. Integrating for 350hr yielded an r.m.s. noise of 6.0 uJy/beam and a resolution of 0.2-0.5". For the first time, wide-field data from the two arrays have been combined in the (u,v) plane and the bandwidth smearing response of the VLA data has been removed. All of the SMGs are detected in our maps as well as sources comprising a non-submm luminous control sample. We find evidence that SMGs are more extended than the general uJy radio population and that therefore, unlike in local ULIRGs, the starburst component of the radio emission is extended and not confined to the galactic nucleus. For the eight sources with redshifts we measure linear sizes between 1 and 8 kpc with a median of 5 kpc. Therefore, they are in general larger than local ULIRGs which may support an early-stage merger scenario for the starburst trigger. X-rays betray AGN in six of the 33 sources in the combined sample. All but one of these are in the control sample, suggesting a lower incidence of AGN amongst the submm-luminous galaxies which is, in turn, consistent with increased X-ray absorption in these dust-obscured starbursts. Only one of our sources is resolved into multiple, distinct components with our high-resolution data. Finally, compared to a previous study of faint radio sources in the GOODS-N field we find systematically smaller source sizes and no evidence for a tail extending to ~4". Possible reasons for this are discussed.

We present the scientific motivation and observing plan for an upcoming detection survey for debris disks using the James Clerk Maxwell Telescope. The SCUBA-2 Unbiased Nearby Stars (SUNS) Survey will observe 500 nearby main sequence and sub-giant stars (100 of each of the A, F, G, K and M spectral classes) to the 850 micron extragalactic confusion limit to search for evidence of submillimeter excess, an indication of circumstellar material. The survey distance boundaries are 8.6, 16.5, 22, 25 and 45 pc for M, K, G, F and A stars, respectively, and all targets lie between the declinations of -40 deg to 80 deg. In this survey, no star will be rejected based on its inherent properties: binarity, presence of planetary companions, spectral type or age. This will be the first unbiased survey for debris disks since IRAS. We expect to detect ~125 debris disks, including ~50 cold disks not detectable in current shorter wavelength surveys. A substantial amount of complementary data will be required to constrain the temperatures and masses of discovered disks. High resolution studies will likely be required to resolve many of the disks. Therefore, these systems will be the focus of future observational studies using a variety of observatories to characterize their physical properties. For non-detected systems, this survey will set constraints (upper limits) on the amount of circumstellar dust, of typically 200 times the Kuiper Belt mass, but as low as 10 times the Kuiper Belt mass for the nearest stars in the sample (approximately 2 pc).

Determining an accurate position for a submm galaxy (SMG) is the crucial step that enables us to move from the basic properties of an SMG sample - source counts and 2-D clustering - to an assessment of their detailed, multi-wavelength properties, their contribution to the history of cosmic star formation and their links with present-day galaxy populations. In this paper, we identify robust radio and/or IR counterparts, and hence accurate positions, for over two thirds of the SCUBA HAlf-Degree Extragalactic Survey (SHADES) Source Catalogue, presenting optical, 24-um and radio images of each SMG. Observed trends in identification rate have given no strong rationale for pruning the sample. Uncertainties in submm position are found to be consistent with theoretical expectations, with no evidence for significant additional sources of error. Employing the submm/radio redshift indicator, via a parameterisation appropriate for radio-identified SMGs with spectroscopic redshifts, yields a median redshift of 2.8 for the radio-identified subset of SHADES, somewhat higher than the median spectroscopic redshift. We present a diagnostic colour-colour plot, exploiting Spitzer photometry, in which we identify regions commensurate with SMGs at very high redshift. Finally, we find that significantly more SMGs have multiple robust counterparts than would be expected by chance, indicative of physical associations. These multiple systems are most common amongst the brightest SMGs and are typically separated by 2-6", or 15-50/(sin i) kpc at z ~ 2, consistent with early bursts seen in merger simulations.

We describe AEGIS20 - a radio survey of the Extended Groth Strip (EGS) conducted with the Very Large Array (VLA) at 1.4GHz. The resulting catalog contains 1,123 emitters and is sensitive to ultraluminous starbursts to z<=1.3, well matched to the redshift range of the DEEP2 spectroscopic survey in this region. We use stacking techniques to explore the uJy-level emission from a variety of galaxy populations selected via conventional criteria - Lyman-break galaxies (LBGs), distant red galaxies (DRGs), UV-selected galaxies and extremely red objects (EROs) - determining their properties as a function of color, magnitude and redshift and their extinction-free contributions to the history of star formation. We confirm the familiar pattern that the star-formation-rate (SFR) density, SFRD, rises by at least ~5x from z=0-1, although we note highly discrepant UV- and radio-based SFR estimates. Our radio-based SFRs become more difficult to interpret at z>1 where correcting for contamination by radio-loud active galactic nuclei (AGN) comes at the price of rejecting luminous starbursts. Whilst stacking radio images is a useful technique, accurate radio-based SFRs for z>>1 galaxies require precise redshifts and extraordinarily high-fidelity radio data to identify and remove accretion-related emission.

In this the first of a series of Letters, we present a description of the panchromatic data sets that have been acquired in the Extended Groth Strip region of the sky. Our survey, the All-wavelength Extended Groth strip International Survey (AEGIS), is intended to study the physical properties and evolutionary processes of galaxies at z ~ 1. It includes the following deep, wide-field imaging data sets: Chandra/ACIS X-ray (0.5 - 10 keV), GALEX ultraviolet (1200 - 2500 Angstrom), CFHT/MegaCam Legacy Survey optical (3600 - 9000 Angstroms), CFHT/CFH12K optical (4500 - 9000 Angstroms), Hubble Space Telescope/ACS optical (4400 - 8500 Angstroms), Palomar/WIRC near-infrared (1.2 - 2.2 microns), Spitzer/IRAC mid-infrared (3.6 - 8.0 microns), Spitzer/MIPS far-infrared (24 - 70 microns), and VLA radio continuum (6 - 20 cm). In addition, this region of the sky has been targeted for extensive spectroscopy using the DEIMOS spectrograph on the Keck II 10 m telescope. Our survey is compared to other large multiwavelength surveys in terms of depth and sky coverage.

We present catalogues of faint 1.4-GHz radio sources from extremely-deep Very Large Array pointings in the Lockman Hole, the Hubble Deep Field North (HDF-N) and ELAIS N2. Our analysis of the HDF-N data has produced maps that are significantly deeper than those previously published and we have used these to search for counterparts to submm sources. For each of the fields we have derived normalised differential source counts and in the case of the HDF-N find no evidence for the previously-reported under-density of sources; our counts are entirely consistent with those found for the majority of other fields. The catalogues are available as an online supplement to this paper and the maps are also available for download.

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

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

Previous VLBI observations of the lens system CLASS B0128+437 showed that the lensed source had a triple sub-component structure. However, one of the images looked very different to the others and it has been suggested that this is due to scatter-broadening in the ISM of the lensing galaxy. We present a new 21-cm EVN map of this system, the low frequency compared to previous epochs chosen as the scattering should be more pronounced at this frequency. Analysis of the new map is at an early stage, but the area of scattering in the affected image is extremely obvious. Other points of interest that we seek to explain are the very large sizes of the images (in conflict with the observed spectral turnover at about 1 GHz) and the different position angle of one of the images at different frequencies.

Spurred by the advent of disk-based recording systems and the nearly explosive increase of internet bandwidth, eVLBI (Parsley et al. 2004) has undergone a remarkable development over the past two years. From ftp-based transfers of small amounts of astronomical data, through near real-time correlation (disk-buffered at the correlator), it has culminated this spring in the first three telescope real-time correlation at JIVE (Onsala, Westerbork and Jodrell Bank). In this paper we will give a review of this development and the current state of affairs. We will also address the current limitations and the way we may improve both bandwidth and reliability and finally we will discuss the opportunities a true high-bandwidth real-time VLBI correlator will provide.

We present the results of phase-referenced VLBA+Effelsberg observations at five frequencies of the gravitational lens B0218+357 to establish the precise registration of the A and B lensed image positions.

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.

Some microquasars are permanently bright radio sources while others are faint but produce powerful radio outbursts. Most of the X-ray binaries (XRB) however are very faint or undetected in the radio regime. The European VLBI Network (EVN) recently introduced the Mark5 recording system which allows data rates of up to 1 Gbit/s. This increases the sensitivity of the array significantly. We briefly describe recent developments in the EVN in terms of reliability of the network and also data quality. We demonstrate the power of the EVN in detecting sub-mJy radio sources with modest integration times (order of hours). This high sensitivity capability will permit the study of variable Galactic sources at milliarcsecond resolution. An estimate is given for the lowest detectable mass of hypothetical intermediate-mass black holes (IMBH) in nearby galaxies, provided these are located in radio-jet systems analogous to microquasars and active galactic nuclei.

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.

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

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

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

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

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

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

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

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

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.

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

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

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.

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

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

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

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