We model a 21 cm intensity mapping survey in the redshift range 0.01<z<1.5 designed to simulate the skies as seen by future radio telescopes such as the Square Kilometre Array (SKA), including instrumental noise and Galactic foregrounds. In our pipeline, we remove the introduced Galactic foregrounds with a fast independent component analysis (fastica) technique. We present the power spectrum of the large-scale matter distribution, C(l), before and after the application of this foreground removal method and calculate the resulting systematic errors. We attempt to reduce systematics in the foreground subtraction by optimally masking the maps to remove high foregrounds in the Galactic plane. Our simulations show a certain level of bias remains in the power spectrum at all scales l<400. At large-scales l<30 this bias is particularly significant. We measure the impact of these systematic effects in two different ways: firstly we fit cosmological parameters to the broadband shape of the power spectrum and secondly we extract the position of the Baryon Acoustic Oscillations (BAO). In the first analysis, we find that the systematics introduce an significant shift in the best fit cosmological parameters at the 2 to 3 sigma level which depends on the masking and noise levels. However, cosmic distances can be recovered in an unbiased way after foreground removal at all simulated redshifts by fitting the BAOs in the power spectrum. We conclude that further advances in foreground removal are needed in order to recover unbiased information from the broadband shape of the power spectrum, however, intensity mapping experiments will be a powerful tool for mapping cosmic distances across a wide redshift range.

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of CO in galaxies at redshift z=0 to z=10. This prediction relies on a recently presented simulation of the molecular cold gas content in ~3e7 evolving galaxies based on the Millennium Simulation. We combine this simulation with a model for the conversion between molecular mass and CO-line intensities, which incorporates the following mechanisms: (i) molecular gas is heated by the CMB, starbursts (SBs), and active galactic nuclei (AGNs); (ii) molecular clouds in dense or inclined galaxies can overlap; (iii) compact gas can attain a smooth distribution in the densest part of disks; (iv) CO-luminosities scale with metallicity changes between galaxies; (v) CO-luminosities are always detected against the CMB. We analyze the relative importance of these effects and predict the cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies (i.e. no AGN, no massive SB) at high z>7 in CO-emission will be dramatically hindered by the weak contrast against the CMB, in contradiction to earlier claims that CMB-heating will ease the detection of high-redshift CO. The full simulation of extragalactic CO-lines and the predicted CO-LFs at any redshift can be accessed online, prior registration required} and they should be useful for the modeling of CO-line surveys with future telescopes, such as ALMA, the LMT, or the SKA.

This is the first paper in a series that present a multi-wavelength analysis of the archetype Ultra-Luminous InfraRed Galaxy (ULIRG) IRAS FSC10214+4724, a gravitationally lensed, starburst/AGN at z=2.3. Here we present a new lens model and spatially-resolved radio data, as well as a deep HST F160W map. The lens modelling employs a Bayesian Markov Chain Monte Carlo algorithm with extended-source, forward ray-tracing. Using these high resolution HST, MERLIN and VLA maps, the algorithm allows us to constrain the level of distortion to the continuum spectral energy distribution resulting from emission components with differing magnification factors, due to their size and proximity to the caustic. Our lens model finds the narrow line region (NLR), and by proxy the active nucleus, is preferentially magnified. This supports previous claims that preferential magnification could mask the expected polycyclic aromatic hydrocarbon spectral features in the Spitzer mid-infrared spectrum which roughly trace the star-forming regions. Furthermore, we show the arc-to-counter-image flux ratio is not a good estimate of the magnification in this system, despite its common use in the IRAS FSC10214+4724 literature. Our lens modelling suggests magnifications of \mu ~ 15-20+-2 for the HST F814W, MERLIN 1.7 GHz and VLA 8 GHz maps, significantly lower than the canonical values of \mu = 50-100 often used for this system. Systematic errors such as the dark matter density slope and co-location of stellar and dark matter centroids dominate the uncertainties in the lens model at the 40 percent level.

We present JVLA observations of the cold (CO (1-0)) molecular gas in IRAS F10214+4724, a lensed ULIRG at z=2.3 with an obscured active nucleus. The galaxy is spatially and spectrally well-resolved in the CO (1-0) emission line. A CO (1-0) counter-image is detected at the 3-sigma level. Five of the 42 km/s channels (with >5-sigma detections) are mapped back into the source plane and their total magnification posterior PDFs sampled. This reveals a roughly linear arrangement, tentatively a rotating disk. We derive a molecular gas mass of M_gas = 1.2 +- 0.2 x 10^10 M_sun, assuming a ULIRG L_{CO}-to-M_{gas} conversion ratio of \alpha = 0.8 M_sun / (K km/s pc^2) that agrees well with the derived range of \alpha = 0.3 - 1.3 for separate dynamical mass estimates at assumed inclinations of i = 90 - 30 degrees. Based on the AGN and CO (1-0) peak emission positions and the lens model, we predict a distortion of the CO Spectral Line Energy Distribution (SLED) where higher order J lines that may be partially excited by AGN heating will be preferentially lensed owing to their smaller solid angles and closer proximity to the AGN and therefore the cusp of the caustic. Comparison with other lensing inversion results shows that the narrow line region and AGN radio core in IRAS F10214+4724 are preferentially lensed by a factor >~ 3 and 11 respectively, relative to the molecular gas emission. This distorts the global continuum emission Spectral Energy Distribution (SED) and suggests caution in unsophisticated uses of IRAS F10214+4724 as an archetype high-redshift ULIRG. We explore two Large Velocity Gradient (LVG) models, incorporating spatial CO (1-0) and (3-2) information and present tentative evidence for an extended, low excitation cold gas component that implies that the total molecular gas mass in IRAS F10214+4724 is a factor >~2 greater than that calculated using spatially unresolved CO observations.

We report 1.7 GHz Very Long Baseline Interferometry (VLBI) observations of IRAS F10214+4724, a lensed z=2.3 obscured quasar with prodigious star formation. We detect what we argue to be the obscured active nucleus with an effective angular resolution of < 50 pc at z = 2.3 . The S_{1.7} = 210 micro-Jy (9-\sigma) detection of this unresolved source is located within the HST rest-frame ultraviolet/optical arc, however, >~100 mas northward of the arc centre of curvature. This leads to a source plane inversion that places the European VLBI Network detection to within milli-arcseconds of the modelled cusp caustic, resulting in a very large magnification (\mu ~70), over an order of magnitude larger than the CO (1-0) derived magnification of a spatially resolved JVLA map, using the same lens model. We estimate the quasar bolometric luminosity from a number of independent techniques and with our X-ray modelling find evidence that the AGN may be close to Compton-thick, with an intrinsic bolometric luminosity log(L_{bol,QSO} / L_sun) = 11.34 +- 0.27 dex. We make the first black hole mass estimate of IRAS F10214+4724 and find log(M_{BH}/M_sun) = 8.36 +- 0.56 which suggests a low black hole accretion rate (\lambda = \dot{M} / \dot{M}_{Edd} ~ 3\pm^7_2 percent). We find evidence for a M_{BH}/M_{spheroid} ratio that is 1-2 orders of magnitude larger than that of submillimetre galaxies (SMGs) at z~2. At face value, this suggests IRAS F10214+4724 has undergone a different evolutionary path compared to SMGs at the same epoch. A primary result of this work is the demonstration that emission regions of differing size and position can undergo significantly different magnification boosts (> 1 dex) and therefore distort our view of high-redshift, gravitationally lensed galaxies.

We describe a 325-MHz survey, undertaken with the Giant Metrewave Radio Telescope (GMRT), which covers a large part of the three equatorial fields at 9, 12 and 14.5 h of right ascension from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) in the area also covered by the Galaxy And Mass Assembly survey (GAMA). The full dataset, after some observed pointings were removed during the data reduction process, comprises 212 GMRT pointings covering ~90 deg^2 of sky. We have imaged and catalogued the data using a pipeline that automates the process of flagging, calibration, self-calibration and source detection for each of the survey pointings. The resulting images have resolutions of between 14 and 24 arcsec and minimum rms noise (away from bright sources) of ~1 mJy/beam, and the catalogue contains 5263 sources brighter than 5 sigma. We investigate the spectral indices of GMRT sources which are also detected at 1.4 GHz and find them to agree broadly with previously published results; there is no evidence for any flattening of the radio spectral index below S_1.4=10 mJy. This work adds to the large amount of available optical and infrared data in the H-ATLAS equatorial fields and will facilitate further study of the low-frequency radio properties of star formation and AGN activity in galaxies out to z~1.

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

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

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

[Abridged] Sub-mm observations of the William Herschel Deep Field using LABOCA revealed possible counterparts for 2 X-ray absorbed QSOs. The aim here is to exploit EVLA imaging at 8.4 GHz to establish the QSOs as radio/sub-mm sources. The challenge in reducing the EVLA data was the presence of a strong 4C source in the field. A new calibration algorithm was applied to the data to subtract it. The resulting thermal noise limited radio map covers the 16'x16' Extended WHDF. It contains 41 sources above a 4-sigma limit, 17 of which have primary beam corrected flux. The radio observations show that the absorbed AGN with LABOCA detections are coincident with radio sources, confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These sources show strong ultraviolet excess (UVX) suggesting the nuclear sightline is gas- but not dust-absorbed. Of the 3 remaining LABOCA sources within the ~5' half-power beam width, 1 is identified with a faint nuclear X-ray/radio source in a nearby galaxy, 1 with a faint radio source and 1 is unidentified in any other band. More generally, differential radio source counts are in good agreement with previous observations, showing at S<50 micro-Jy a significant excess over a pure AGN model. In the full area, of 10 sources fainter than this limit, 6 have optical counterparts of which 3 are UVX (i.e. likely QSOs) including the 2 absorbed quasar LABOCA sources. The other faint radio counterparts are not UVX but are only slightly less blue and likely to be star-forming/merging galaxies, predominantly at lower luminosities and redshifts. The 4 faint, optically unidentified radio sources may be either dust obscured QSOs or galaxies. These high-z obscured AGN and lower-z star-forming populations are thus the main candidates to explain the observed excess in faint source counts and hence the excess radio background found previously by the ARCADE2 experiment.

We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.

The benefits for the e-MERLIN and EVN arrays of using antennae at the satellite communication station at Goonhilly in Cornwall are discussed. The location of this site - new to astronomy - will provide an almost equal distribution of long baselines in the east-west- and north-south directions, and opens up the possibility to get significantly improved observations of equatorial fields with e-MERLIN. These additional baselines will improve the sensitivity on a set of critical spatial scales and will increase the angular resolution of e-MERLIN by a factor of two. e-MERLIN observations, including many allocated under the e-MERLIN Legacy programme, will benefit from the enhanced angular resolution and imaging capability especially for sources close to or below the celestial equator (where ESO facilities such as ALMA will operate) of including the Goonhilly telescopes. Furthermore, the baselines formed between Goonhilly and the existing stations will close the gap between the baselines of e-MERLIN and those of the European VLBI Network (EVN) and therefore enhance the legacy value of e-MERLIN datasets.

A consortium of universities has recently been formed with the goal of using the decommissioned telecommunications infrastructure at the Goonhilly Earth Station in Cornwall, UK, for astronomical purposes. One particular goal is the introduction of one or more of the ~30-metre parabolic antennas into the existing e-MERLIN radio interferometer. This article introduces this scheme and presents some simulations which quantify the improvements that would be brought to the e-MERLIN system. These include an approximate doubling of the spatial resolution of the array, an increase in its N-S extent with strong implications for imaging the most well-studied equatorial fields, accessible to ESO facilities including ALMA. It also increases the overlap between the e-MERLIN array and the European VLBI Network. We also discuss briefly some niche science areas in which an e-MERLIN array which included a receptor at Goonhilly would be potentially world-leading, in addition to enhancing the existing potential of e-MERLIN in its role as a Square Kilometer Array pathfinder instrument.

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

Gravitational lensing is a powerful astrophysical and cosmological probe and is particularly valuable at submillimeter wavelengths for the study of the statistical and individual properties of dusty starforming galaxies. However the identification of gravitational lenses is often time-intensive, involving the sifting of large volumes of imaging or spectroscopic data to find few candidates. We used early data from the Herschel Astrophysical Terahertz Large Area Survey to demonstrate that wide-area submillimeter surveys can simply and easily detect strong gravitational lensing events, with close to 100% efficiency.

We use the Herschel-ATLAS science demonstration data to investigate the star-formation properties of radio-selected galaxies in the GAMA-9h field as a function of radio luminosity and redshift. Radio selection at the lowest radio luminosities, as expected, selects mostly starburst galaxies. At higher radio luminosities, where the population is dominated by AGN, we find that some individual objects are associated with high far-infrared luminosities. However, the far-infrared properties of the radio-loud population are statistically indistinguishable from those of a comparison population of radio-quiet galaxies matched in redshift and K-band absolute magnitude. There is thus no evidence that the host galaxies of these largely low-luminosity (Fanaroff-Riley class I), and presumably low-excitation, AGN, as a population, have particularly unusual star-formation histories. Models in which the AGN activity in higher-luminosity, high-excitation radio galaxies is triggered by major mergers would predict a luminosity-dependent effect that is not seen in our data (which only span a limited range in radio luminosity) but which may well be detectable with the full Herschel-ATLAS dataset.

To investigate the poorly constrained sub-mm counts and spectral properties of blazars we searched for these in the Herschel-ATLAS (H-ATLAS) science demostration phase (SDP) survey catalog. We cross-matched 500$\mu$m sources brighter than 50 mJy with the FIRST radio catalogue. We found two blazars, both previously known. Our study is among the first blind blazar searches at sub-mm wavelengths, i.e., in the spectral regime where little is still known about the blazar SEDs, but where the synchrotron peak of the most luminous blazars is expected to occur. Our early results are consistent with educated extrapolations of lower frequency counts and question indications of substantial spectral curvature downwards and of spectral upturns at mm wavelengths. One of the two blazars is identified with a Fermi/LAT $\gamma$-ray source and a WMAP source. The physical parameters of the two blazars are briefly discussed.These observations demonstrate that the H-ATLAS survey will provide key information about the physics of blazars and their contribution to sub-mm counts.

To exploit synergies between the Herschel Space Observatory and next generation radio facilities, we have extended the semi-empirical extragalactic radio continuum simulation of Wilman et al. (2008) to the mid- and far-infrared. Here we describe the assignment of infrared spectral energy distributions (SEDs) to the star-forming galaxies and active galactic nuclei, using Spitzer 24, 70 and 160 micron and SCUBA 850 micron survey results as the main constraints. Star-forming galaxies dominate the source counts, and a model in which their far-infrared-radio correlation and infrared SED assignment procedure are invariant with redshift underpredicts the observed 24 and 70 micron source counts. The 70 micron deficit can be eliminated if the star-forming galaxies undergo stronger luminosity evolution than originally assumed for the radio simulation, a requirement which may be partially ascribed to known non-linearity in the far-infrared--radio correlation at low luminosity if it evolves with redshift. At 24 micron, the shortfall is reduced if the star-forming galaxies develop SEDs with cooler dust and correspondingly stronger Polycyclic Aromatic Hydrocarbon (PAH) emission features with increasing redshift at a given far-infrared luminosity, but this trend may reverse at z>1 in order not to overproduce the sub-mm source counts. The resulting model compares favourably with recent BLAST results and we have extended the simulation database to aid the interpretation of Herschel surveys. Such comparisons may also facilitate further model refinement and revised predictions for the SKA and its precursors.

The current status of the HI simulation efforts is presented, in which a self consistent simulation path is described and basic equations to calculate array sensitivities are given. There is a summary of the SKA Design Study (SKADS) sky simulation and a method for implementing it into the array simulator is presented. A short overview of HI sensitivity requirements is discussed and expected results for a simulated HI survey are presented.

We present a range of new observations of the "broad absorption line radio galaxy" 4C +72.26 (z = 3.5), including sensitive rest-frame ultraviolet integral field spectroscopy using the Gemini/GMOS-N instrument and Subaru/CISCO K-band imaging and spectroscopy. We show that 4C +72.26 is a system of two vigorously star-forming galaxies superimposed along the line of sight separated by ~1300 +/- 200 km/s in velocity, with each demonstrating spectroscopically resolved absorption lines. The most active star-forming galaxy also hosts the accreting supermassive black hole which powers the extended radio source. We conclude that the star formation is unlikely to have been induced by a shock caused by the passage of the radio jet, and instead propose that a collision is a more probable trigger for the star formation. Despite the massive starburst, the UV-mid-infrared spectral energy distribution suggests that the pre-existing stellar population comprises ~10^12 Msolar of stellar mass, with the current burst only contributing a further ~2%, suggesting that 4C +72.26 has already assembled most of its final stellar mass.

We present the results of a search for galaxy clusters in Subaru-XMM Deep Field. We reach a depth for a total cluster flux in the 0.5-2 keV band of 2x10^{-15} ergs cm^{-2} s^{-1} over one of the widest XMM-Newton contiguous raster surveys, covering an area of 1.3 square degrees. Cluster candidates are identified through a wavelet detection of extended X-ray emission. The red sequence technique allows us to identify 57 cluster candidates. We report on the progress with the cluster spectroscopic follow-up and derive their properties based on the X-ray luminosity and cluster scaling relations. In addition, 3 sources are identified as X-ray counterparts of radio lobes, and in 3 further sources, X-ray counterpart of radio lobes provides a significant fraction of the total flux of the source. In the area covered by NIR data, our identification success rate achieves 86%. We detect a number of radio galaxies within our groups and for a luminosity-limited sample of radio galaxies we compute halo occupation statistics using a marked cluster mass function. We compare the cluster detection statistics in the SXDF with the predictions of concordance cosmology and current knowledge of the X-ray cluster properties, concluding that a reduction of concordance sigma_8 value by 5% is required in order to match the prediction of the model and the data. This conclusion still needs verification through the completion of cluster follow-up.

The S3-Tools are a set of Python-based routines and interfaces whose purpose is to provide user-friendly access to the SKA Simulated Skies (S3) set of simulations, an effort led by the University of Oxford in the framework of the European Union's SKADS program (http://www.skads-eu.org). The databases built from the S3 simulations are hosted by the Oxford e-Research Center (OeRC), and can be accessed through a web portal at http://s-cubed.physics.ox.ac.uk. This paper focuses on the practical steps involved to make radio images from the S3-SEX and S3-SAX simulations using the S3-Map tool and should be taken as a broad overview. For a more complete description, the interested reader should look up the user's guide. The output images can then be used as input to instrument simulators, e.g. to assess technical designs and observational strategies for the SKA and SKA pathfinders.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

Infrared spectroscopy of the H-alpha emission lines of a sub-sample of 19 high-redshift (0.8 < z < 2.3) Molonglo quasars, selected at 408 MHz, is presented. These emission lines are fitted with composite models of broad and narrow emission, which include combinations of classical broad-line regions of fast-moving gas clouds lying outside the quasar nucleus, and/or a theoretical model of emission from an optically-thick, flattened, rotating accretion disk. All bar one of the nineteen sources are found to have emission consistent with the presence of an optically-emitting accretion disk, with the exception appearing to display complex emission including at least three broad components. Ten of the quasars have strong Bayesian evidence for broad-line emission arising from an accretion disk together with a standard broad-line region, selected in preference to a model with two simple broad lines. Thus the best explanation for the complexity required to fit the broad H-alpha lines in this sample is optical emission from an accretion disk in addition to a region of fast-moving clouds. We derive estimates of the angle between the rotation axis of the accretion disk and the line of sight. A weak correlation is found between the accretion disk angle and the logarithm of the low-frequency radio luminosity. This is direct, albeit tenuous, evidence for the receding torus model. Velocity shifts of the broad H-alpha components are analysed and the results found to be consistent with a two-component model comprising one single-peaked broad line emitted at the same redshift as the narrow lines, and emission from an accretion disk which appears to be preferentially redshifted with respect to the narrow lines for high-redshift sources and blueshifted relative to the narrow lines for low-redshift sources.

We present the first results from the largest deep extragalactic millimetre-wavelength survey undertaken to date. These results are derived from maps covering over 0.7 deg^2, made at 1.1mm, using the AzTEC continuum camera mounted on the James Clerk Maxwell Telescope. The maps were made in the two fields originally targeted at 0.85mm with SCUBA in the SHADES project, namely the Lockman Hole East (mapped to a depth of 0.9-1.3 mJy rms) and the Subaru XMM Deep Field (1.0-1.7 mJy rms). The wealth of existing and forthcoming deep multi-frequency data in these two fields will allow the bright mm source population revealed by these images to be explored in detail in subsequent papers. Here we present the maps themselves, a catalogue of 114 high-significance sub-millimetre galaxy detections, and a thorough statistical analysis leading to the most robust determination to date of the 1.1mm source number counts. Through careful comparison, we find that both the COSMOS and GOODS North fields, also imaged with AzTEC, contain an excess of mm sources over the new 1.1mm source-count baseline established here. In particular, our new AzTEC/SHADES results indicate that very luminous high-redshift dust enshrouded starbursts (S_{1.1} > 3 mJy) are 25-50% less common than would have been inferred from these smaller surveys, thus highlighting the potential roles of cosmic variance and clustering in such measurements. We compare number count predictions from recent models of the evolving mm/sub-mm source population to these SMG surveys, which provide important constraints for the ongoing refinement of semi-analytic and hydrodynamical models of galaxy formation, and find that all recent models over-predict the number of bright sub-millimetre galaxies found in this survey.