We present the VLA-COSMOS 3 GHz Large Project based on 384 hours of observations with the Karl G. Jansky Very Large Array (VLA) at 3 GHz (10 cm) toward the two square degree Cosmic Evolution Survey (COSMOS) field. The final mosaic reaches a median rms of 2.3 uJy/beam over the two square degrees at an angular resolution of 0.75". To fully account for the spectral shape and resolution variations across the broad (2 GHz) band, we image all data with a multiscale, multifrequency synthesis algorithm. We present a catalog of 10,830 radio sources down to 5 sigma, out of which 67 are combined from multiple components. Comparing the positions of our 3 GHz sources with those from the Very Long Baseline Array (VLBA)-COSMOS survey, we estimate that the astrometry is accurate to 0.01" at the bright end (signal-to-noise ratio, S/N_3GHz > 20). Survival analysis on our data combined with the VLA-COSMOS 1.4~GHz Joint Project catalog yields an expected median radio spectral index of alpha=-0.7. We compute completeness corrections via Monte Carlo simulations to derive the corrected 3 GHz source counts. Our counts are in agreement with previously derived 3 GHz counts based on single-pointing (0.087 square degrees) VLA data. In summary, the VLA-COSMOS 3 GHz Large Project simultaneously provides the largest and deepest radio continuum survey at high (0.75") angular resolution to date, bridging the gap between last-generation and next-generation surveys.

We study the physical properties of a sample of 6 SMGs in the COSMOS field, spectroscopically confirmed to lie at z>4. We use new GMRT 325 MHz and 3 GHz JVLA data to probe the rest-frame 1.4 GHz emission at z=4, and to estimate the sizes of the star-forming (SF) regions of these sources, resp. Combining our size estimates with those available in the literature for AzTEC1 and AzTEC3 we infer a median radio-emitting size for our z>4 SMGs of (0.63"+/-0.12")x(0.35"+/-0.05") or 4.1x2.3 kpc^2 (major times minor axis; assuming z=4.5) or lower if we take the two marginally resolved SMGs as unresolved. This is consistent with the sizes of SF regions in lower-redshift SMGs, and local normal galaxies, yet higher than the sizes of SF regions of local ULIRGs. Our SMG sample consists of a fair mix of compact and more clumpy systems with multiple, perhaps merging, components. With an average formation time of ~280 Myr, derived through modeling of the UV-IR SEDs, the studied SMGs are young systems. The average stellar mass, dust temperature, and IR luminosity we derive are M*~1.4x10^11 M_sun, T_dust~43 K, and L_IR~1.3x10^13L_sun, resp. The average L_IR is up to an order of magnitude higher than for SMGs at lower redshifts. Our SMGs follow the correlation between dust temperature and IR luminosity as derived for Herschel-selected 0.1<z<2 galaxies. We study the IR-radio correlation for our sources and find a deviation from that derived for z<3 ULIRGs (<q_IR>=1.95+/-0.26 for our sample, compared to q~2.6 for IR luminous galaxies at z<2). In summary, we find that the physical properties derived for our z>4 SMGs put them at the high end of the L_IR-T_dust distribution of SMGs, and that our SMGs form a morphologically heterogeneous sample. Thus, further in-depth analyses of large, statistical samples of high-redshift SMGs are needed to fully understand their role in galaxy formation and evolution.

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.

The Square Kilometre Array (SKA) will have a low frequency component (SKA-low) which has as one of its main science goals the study of the redshifted 21cm line from the earliest phases of star and galaxy formation in the Universe. This 21cm signal provides a new and unique window on both the formation of the first stars and accreting black holes and the later period of substantial ionization of the intergalactic medium. The signal will teach us fundamental new things about the earliest phases of structure formation, cosmology and even has the potential to lead to the discovery of new physical phenomena. Here we present a white paper with an overview of the science questions that SKA-low can address, how we plan to tackle these questions and what this implies for the basic design of the telescope.

The European VLBI Network (EVN) has been used at two epochs in 2003 and 2009 to obtain multi-frequency high-resolution images of the merger galaxy NGC6240 in order to study the radio properties of all compact high-brightness components in the galaxy. Our observations at milli-arcsecond resolution detected the northern and southern nuclei and two radio components, which we interpret as long-lived luminous supernovae associated with the circum-nuclear starburst activity at the southern nucleus. The new VLBI data support the presence of an active galactic nucleus (AGN) together with starburst activity at the southern nucleus and provides some evidence for an AGN at the northern nucleus. The two nuclei both display an inverted spectrum at lower GHz frequencies. The spectrum of the southern nucleus indicates thermal free-free absorption on parsec scales, consistent with the presence of an AGN.

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.

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

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.

We present a sky simulation of the atomic HI emission line and the first ten CO rotational emission lines of molecular gas in galaxies beyond the Milky Way. The simulated sky field has a comoving diameter of 500/h Mpc, hence the actual field-of-view depends on the (user-defined) maximal redshift zmax; e.g. for zmax=10, the field of view yields ~4x4 sqdeg. For all galaxies, we estimate the line fluxes, line profiles, and angular sizes of the HI and CO emission lines. The galaxy sample is complete for galaxies with cold hydrogen masses above 10^8 Msun. This sky simulation builds on a semi-analytic model of the cosmic evolution of galaxies in a Lambda-cold dark matter (LCDM) cosmology. The evolving CDM-distribution was adopted from the Millennium Simulation, an N-body CDM-simulation in a cubic box with a side length of 500/h Mpc. This side length limits the coherence scale of our sky simulation: it is long enough to allow the extraction of the baryon acoustic oscillations (BAOs) in the galaxy power spectrum, yet the position and amplitude of the first acoustic peak will be imperfectly defined. This sky simulation is a tangible aid to the design and operation of future telescopes, such the SKA, the LMT, and ALMA. The results presented in this paper have been restricted to a graphical representation of the simulated sky and fundamental dN/dz-analyzes for peak flux density limited and total flux limited surveys of HI and CO. A key prediction is that HI will be harder to detect at redshifts z>2 than predicted by a no-evolution model. The future verification or falsification of this prediction will allow us to qualify the semi-analytic models.

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.

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.

We present observations at 1.2 mm with MAMBO-II of a sample of z>~2 radio-intermediate obscured quasars, as well as CO observations of two sources with the Plateau de Bure Interferometer. Five out of 21 sources (24%) are detected at a significance of >=3sigma. Stacking all sources leads to a statistical detection of <S_1.2mm>= 0.96+-0.11 mJy and stacking only the non-detections also yields a statistical detection, with <S_1.2mm>= 0.51+-0.13 mJy. This corresponds to a typical far-infrared luminosity L_FIR~4x10^12 Lsol. If the far-infrared luminosity is powered entirely by star-formation, and not by AGN-heated dust, then the characteristic inferred star-formation rate is ~700 Msol yr-1. This far-infrared luminosity implies a dust mass of M_dust~3x10^8 Msol. We estimate that such large dust masses on kpc scales can plausibly cause the obscuration of the quasars. We present dust SEDs for our sample and derive a mean SED for our sample. This mean SED is not well fitted by clumpy torus models, unless additional extinction and far-infrared re-emission due to cool dust are included. There is a hint that the host galaxies of obscured quasars must have higher far-infrared luminosities and cool-dust masses and are therefore often found at an earlier evolutionary phase than those of unobscured quasars. For one source at z=2.767, we detect the CO(3-2) transition, with S_CO Delta nu=630+-50 mJy km s-1, corresponding to L_CO(3-2)= 3.2x10^7 Lsol, or L'_CO(3-2)=2.4x10^10 K km s-1 pc2. For another source at z=4.17, the lack of detection of the CO(4-3) line yields a limit of L'_CO(4-3)<1x10^10 K km s-1 pc2. Molecular gas masses, gas depletion timescales and gas-to-dust ratios are estimated (Abridged).

The radio properties of 11 obscured `radio-intermediate' quasars at redshifts z>~2 have been investigated using the European Very-Long-Baseline-Interferometry Network (EVN) at 1.66 GHz. A sensitivity of ~25 micro Jy per 14x17 mas2 beam was achieved, and in 7 out of 11 sources unresolved radio emission was securely detected. The detected radio emission of each source accounts for ~30-100 % of the total source flux density. The physical extent of this emission is ~<150 pc, and the derived properties indicate that this emission originates from an active galactic nucleus (AGN). The missing flux density is difficult to account for by star-formation alone, so radio components associated with jets of physical size >~150 pc, and ~< 40 kpc are likely to be present in most of the sources. Amongst the observed sample steep, flat, gigahertz-peaked and compact-steep spectrum sources are all present. Hence, as well as extended and compact jets, examples of beamed jets are also inferred, suggesting that in these sources, the obscuration must be due to dust in the host galaxy, rather than the torus invoked by the unified schemes. Comparing the total to core (~< 150 pc) radio luminosities of this sample with different types of AGN suggests that this sample of z >~2 radio-intermediate obscured quasars shows radio properties that are more similar to those of the high-radio-luminosity end of the low-redshift radio-quiet quasar population than those of FR I radio galaxies. This conclusion may reflect intrinsic differences, but could be strongly influenced by the increasing effect of inverse-Compton cooling of extended radio jets at high redshift.

We present simulated images of Supernova 1993J at 8.4 GHz using Very Long Baseline Interferometry (VLBI) techniques. A spherically symmetric source model is convolved with realistic uv-plane distributions, together with standard imaging procedures, to assess the extent of instrumental effects on the recovered brightness distribution. In order to facilitate direct comparisons between the simulations and published VLBI images of SN1993J, the observed uv-coverage is determined from actual VLBI observations made in the years following its discovery. The underlying source model only exhibits radial variation in its density profile, with no azimuthal dependence and, even though this model is morphologically simple, the simulated VLBI observations qualitatively reproduce many of the azimuthal features of the reported VLBI observations, such as appearance and evolution of complex azimuthal structure and apparent rotation of the shell. We demonstrate that such features are inexorably coupled to the uv-plane sampling. The brightness contrast between the peaks and the surrounding shell material are not as prominent in the simulations (which of course assume no antenna- or baseline-based amplitude or phase errors, meaning no self-calibration procedures will have incorporated any such features in models). It is conclusive that incomplete uv-plane sampling has a drastic effect on the final images for observations of this nature. Difference imaging reveals residual emission up to the 8 sigma level. Extreme care should be taken when using interferometric observations to directly infer the structure of objects such as supernovae.

We have developed a semi-empirical simulation of the extragalactic radio continuum sky suitable for aiding the design of next generation radio interferometers such as the Square Kilometre Array (SKA). The emphasis is on modelling the large-scale cosmological distribution of radio sources rather than the internal details of individual galaxies. Here we provide a description of the simulation to accompany the online release of a catalogue of 320 million simulated radio sources. The simulation covers 20x20 deg^2 - a plausible upper limit to the instantaneous field of view attainable with future (e.g. SKA) aperture array technologies - out to redshift z=20, and down to flux density limits of 10 nJy at 151, 610 MHz, 1.4, 4.86 and 18 GHz. Five distinct source types are included: radio-quiet AGN, radio-loud AGN of the FRI and FRII structural classes, and star-forming galaxies, the latter split into populations of quiescent and starbursting galaxies. In our semi-empirical approach, the simulated sources are drawn from observed (or extrapolated) luminosity functions and grafted onto an underlying dark matter density field with biases which reflect their measured large-scale clustering. A numerical Press-Schechter-style filtering of the density field is used to identify and populate clusters of galaxies. Radio source structures are built from point source and elliptical sub-components, and for FRI and FRII sources an orientation-based unification and beaming model is used to partition flux between the core and extended lobes and hotspots. The simulation output can be post-processed to achieve more complete agreement with observational data in the years ahead, with the aim of using these 'idealised skies' in telescope simulators to optimise the design of the SKA itself (abridged).

We discuss the equatorial imaging benefits that arise from the addition of the 25-metre dish at Chilbolton to the e-MERLIN array. Its inclusion considerably enhances the capabilities of e-MERLIN on and below the equator. This will become particularly important in the era of ALMA and other upcoming southern hemisphere facilities. We present simulated observations of point sources in the equatorial region of the sky which is the target area for many existing sky surveys. We find that the additional baselines created by the inclusion of the Chilbolton dish favourably adjust the beam shape of e-MERLIN to a more compact and circular shape, with significantly reduced sidelobe structure. Putting aside the benefits of increased collecting area, the modified beam shape has implications for more rapidly reaching a given completeness limit for equatorial surveys.

We present basic observational data on the 6C** sample. This is a new sample of radio sources drawn from the 151 MHz 6C survey, which was filtered with radio criteria chosen to optimize the chances of finding radio galaxies at z > 4. The filtering criteria are a steep-spectral index and a small angular size. The final sample consists of 68 sources from a region of sky covering 0.421 sr. We present VLA radio maps, and the results of K-band imaging and optical spectroscopy. Near-infrared counterparts are identified for 66 of the 68 sources, down to a 3-sigma limiting magnitude of K ~ 22 mag in a 3-arcsec aperture. Eight of these identifications are spatially compact, implying an unresolved nuclear source. The K-magnitude distribution peaks at a median K=18.7 mag, and is found to be statistically indistinguishable from that of the similarly selected 6C* sample, implying that the redshift distribution could extend to z > 4. Redshifts determined from spectroscopy are available for 22 (32 per cent) of the sources, over the range of 0.2 < z < 3.3 . We measure 15 of these, whereas the other 7 were previously known. Six sources are at z > 2.5. Four sources show broad emission lines in their spectra and are classified as quasars. Three of these show also an unresolved K-band identification. Eleven sources fail to show any distinctive emission and/or absorption features in their spectra. We suggest that these could be (i) in the so-called `redshift desert' region of 1.2 < z < 1.8, or (ii) at a greater redshift, but feature weak emission line spectra.

We present multi-frequency radio observations of a sample of z~2 obscured (type-2) quasars in the Spitzer extragalactic First Look Survey area. We combine the public data at 1.4 GHz, used in the selection of these sources, with new observations at 610 MHz (GMRT) and at 4.9 GHz (VLA). We find the sample includes sources with steep, flat and gigahertz-peaked spectra. There are no strong correlations between the presence or absence of emission lines in the optical spectra and the radio spectral properties of the sample. However, there are no secure flat-spectrum type-2 quasars with narrow emission lines which would be problematic for unified schemes. Most of the population have straight radio spectra with spectral index alpha~1 as is expected for developed, potentially FRI-like, jets in which continous injection of relativistic electrons is accompanied by inverse-Compton losses against the cosmic microwave background.

Broad hydroxyl (OH) absorption-lines in the 1667 MHz and 1665 MHz transition towards the central region of NGC 3079 have been observed at high resolution with the European VLBI Network (EVN). Velocity fields of two OH absorption components were resolved across the unresolved nuclear radio continuum of ~10 parsecs. The velocity field of the OH absorption close to the systemic velocity shows rotation in nearly the same sense as the edge-on galactic-scale molecular disk probed by CO(1-0) emission. The velocity field of the blue-shifted OH absorption displays a gradient in almost the opposite direction. The blue-shifted velocity field represents a non-rotational component, which may trace an outflow from the nucleus, or material driven and shocked by the kiloparsec-scale superbubble. This OH absorption component traces a structure that does not support a counter-rotating disk suggested on the basis of the neutral hydrogen absorption.