This white paper describes the science case for Very Long Baseline Interferometry (VLBI) and provides suggestions towards upgrade paths for the European VLBI Network (EVN). The EVN is a distributed long-baseline radio interferometric array, that operates at the very forefront of astronomical research. Recent results, together with the new science possibilities outlined in this vision document, demonstrate the EVN's potential to generate new and exciting results that will transform our view of the cosmos. Together with e-MERLIN, the EVN provides a range of baseline lengths that permit unique studies of faint radio sources to be made over a wide range of spatial scales. The science cases are reviewed in six chapters that cover the following broad areas: cosmology, galaxy formation and evolution, innermost regions of active galactic nuclei, explosive phenomena and transients, stars and stellar masers in the Milky Way, celestial reference frames and space applications. The document concludes with identifying the synergies with other radio, as well as multi-band/multi-messenger instruments, and provide the recommendations for future improvements. The appendices briefly describe other radio VLBI arrays, the technological framework for EVN developments, and a selection of spectral lines of astrophysical interest below 100 GHz. The document includes a glossary for non-specialists, and a list of acronyms at the end.

We propose a probabilistic framework for performing simultaneous estimation of source structure and fringe-fitting parameters in Very Long Baseline Interferometry (VLBI) observations. As a first step, we demonstrate this technique through the analysis of synthetic short-duration Event Horizon Telescope (EHT) observations of various geometric source models at 230 GHz, in the presence of baseline-dependent thermal noise. We perform Bayesian parameter estimation and model selection between the different source models to obtain reliable uncertainty estimates and correlations between various source and fringe-fitting related model parameters. We also compare the Bayesian posteriors with those obtained using widely-used VLBI data reduction packages such as CASA and AIPS, by fringe-fitting 200 Monte Carlo simulations of each source model with different noise realisations, to obtain distributions of the Maximum A Posteriori (MAP) estimates. We find that, in the presence of resolved asymmetric source structure and a given array geometry, the traditional practice of fringe-fitting with a point source model yields appreciable offsets in the estimated phase residuals, potentially biasing or limiting the dynamic range of the starting model used for self-calibration. Simultaneously estimating the source structure earlier in the calibration process with formal uncertainties improves the precision and accuracy of fringe-fitting and establishes the potential of the available data especially when there is little prior information. We also note the potential applications of this method to astrometry and geodesy for specific science cases and the planned improvements to the computational performance and analyses of more complex source distributions.

We present MeerKAT 1.28 GHz total-intensity, polarization, and spectral-index images covering the giant (projected length $l \approx 1.57$~Mpc) X-shaped radio source PKS~2014$-$55 with an unprecedented combination of brightness sensitivity and angular resolution. They show the clear "double boomerang" morphology of hydrodynamical backflows from the straight main jets deflected by the large and oblique hot-gas halo of the host galaxy PGC~064440. The magnetic field orientation in PKS~2014$-$55 follows the flow lines from the jets through the secondary wings. The radio source is embedded in faint ($T_\mathrm{b} \approx 0.5 \mathrm{\,K}$) cocoons having the uniform brightness temperature and sharp outer edges characteristic of subsonic expansion into the ambient intra-group medium. The position angle of the much smaller ($l \sim 25$~kpc) restarted central source is within $5^\circ$ of the main jets, ruling out models that invoke jet re-orientation or two independent jets. Compression and turbulence in the backflows probably produce the irregular and low polarization bright region behind the apex of each boomerang as well as several features in the flow with bright heads and dark tails.

Many low-luminosity active galactic nuclei (AGNs) contain a compact radio core which can be observed with high angular resolution using very long baseline interferometry (VLBI). Combining arcsec-scale structural information with milliarcsec-resolution VLBI imaging is a useful way to characterise the objects and to find compact cores on parsec scales. VLBI imaging could also be employed to look for dual AGNs when the sources show kpc-scale double symmetric structure with flat or inverted radio spectra. We observed five such sources at redshifts 0.36 < z < 0.58 taken from an optically selected sample of Type 2 quasars with the European VLBI Network (EVN) at 1.7 and 5 GHz. Out of the five sources, only one (SDSS J1026-0042) shows a confidently detected compact VLBI core at both frequencies. The other four sources are marginally detected at 1.7 GHz only, indicating resolved-out radio structure and steep spectra. Using first-epoch data from the ongoing Karl G. Jansky Very Large Array Sky Survey, we confirm that indeed all four of these sources have steep radio spectra on arcsec scale, contrary to the inverted spectra reported earlier in the literature. However, the VLBI-detected source, SDSS J1026-0042, has a flat integrated spectrum. Radio AGNs that show kpc-scale symmetric structures with truly flat or inverted spectra could still be promising candidates of dual AGNs, to be targeted with VLBI observations in the future.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.

In April 2019, the Event Horizon Telescope (EHT) collaboration revealed the first image of the candidate super-massive black hole (SMBH) at the centre of the giant elliptical galaxy Messier 87 (M87). This event-horizon-scale image shows a ring of glowing plasma with a dark patch at the centre, which is interpreted as the shadow of the black hole. This breakthrough result, which represents a powerful confirmation of Einstein's theory of gravity, or general relativity, was made possible by assembling a global network of radio telescopes operating at millimetre wavelengths that for the first time included the Atacama Large Millimeter/ submillimeter Array (ALMA). The addition of ALMA as an anchor station has enabled a giant leap forward by increasing the sensitivity limits of the EHT by an order of magnitude, effectively turning it into an imaging array. The published image demonstrates that it is now possible to directly study the event horizon shadows of SMBHs via electromagnetic radiation, thereby transforming this elusive frontier from a mathematical concept into an astrophysical reality. The expansion of the array over the next few years will include new stations on different continents - and eventually satellites in space. This will provide progressively sharper and higher-fidelity images of SMBH candidates, and potentially even movies of the hot plasma orbiting around SMBHs. These improvements will shed light on the processes of black hole accretion and jet formation on event-horizon scales, thereby enabling more precise tests of general relativity in the truly strong field regime.

With the publication of Gaia DR2, 1.3 billion stars now have public parallax and proper motion measurements. In this contribution, we compare the results for sources that have both optical and radio measurements, focusing on circumstellar masers. For these large, variable and bright AGB stars, the VLBI astrometry results can be more robust. Moreover, there are a number of applications where VLBI astrometry provides unique data for studying stellar populations and Galactic structure. The BeSSel project not only provides parallax and proper motions at much larger distances than Gaia can reach, but it also uniquely samples the spiral arms of the Galaxy. The evolved stars in the BAaDE sample can potentially constrain the dynamics and stellar content of the inner bulge and bar of the Milky Way, not reachable in the optical.

We report interferometric observations tuned to the redshifted neutral hydrogen (HI) 21cm emission line in three strongly lensed galaxies at $z \sim 0.4$ with the Giant Metrewave Radio Telescope (GMRT). One galaxy spectrum (J1106+5228 at z=0.407) shows evidence of a marginal detection with an integrated signal-to-noise ratio of 3.8, which, if confirmed by follow-up observations, would represent the first strongly lensed and most distant individual galaxy detected in HI emission. Two steps are performed to transcribe the lensed integrated flux measurements into HI mass measurements for all three target galaxies. First, we calculate the HI magnification factor $\mu$ by applying general relativistic ray-tracing to a physical model of the source-lens system. The HI magnification generally differs from the optical magnification and depends largely on the intrinsic HI mass $M_{\rm HI}$ due to the HI mass-size relation. Second, we employ a Bayesian formalism to convert the integrated flux, amplified by the $M_{\rm HI}$-dependent magnification factor $\mu$, into a probability density for $M_{\rm HI}$, accounting for the asymmetric uncertainty due to the declining HI mass function (Eddington bias). In this way, we determine a value of $\log_{\rm 10} (M_{\rm HI}/M_\odot) = 10.2^{+0.3}_{-0.7}$ for J1106+5228, consistent with the estimate of $9.4\pm0.3$ from the optical properties of this galaxy. The HI mass of the other two sources are consistent with zero within a 95 per cent confidence interval however we still provide upper limits for both sources and a $1\sigma$ lower limit for J1250-0135 using the same formalism.

The past decade has seen significant advances in cm-wave VLBI extragalactic observations due to a wide range of technical successes, including the increase in processed field-of-view and bandwidth. The future inclusion of MeerKAT into global VLBI networks would provide further enhancement, particularly the dramatic sensitivity boost to >7000 km baselines. This will not be without its limitations, however, considering incomplete MeerKAT band overlap with current VLBI arrays and the small (real-time) field-of-view afforded by the phased up MeerKAT array. We provide a brief overview of the significant contributions MeerKAT-VLBI could make, with an emphasis on the scientific output of several MeerKAT extragalactic Large Survey Projects.

The MeerKAT telescope will be one of the most sensitive radio arrays in the pre-SKA era. Here we discuss a low-frequency SZ-selected cluster survey with MeerKAT, the MeerKAT Extended Relics, Giant Halos, and Extragalactic Radio Sources (MERGHERS) survey. The primary goal of this survey is to detect faint signatures of diffuse cluster emission, specifically radio halos and relics. SZ-selected cluster samples offer a homogeneous, mass-limited set of targets out to higher redshift than X-ray samples. MeerKAT is sensitive enough to detect diffuse radio emission at the faint levels expected in low-mass and high-redshift clusters, thereby enabling radio halo and relic formation theories to be tested with a larger statistical sample over a significantly expanded phase space. Complementary multiwavelength follow-up observations will provide a more complete picture of any clusters found to host diffuse emission, thereby enhancing the scientific return of the MERGHERS survey.

The MIGHTEE large survey project will survey four of the most well-studied extragalactic deep fields, totalling 20 square degrees to $\mu$Jy sensitivity at Giga-Hertz frequencies, as well as an ultra-deep image of a single ~1 square degree MeerKAT pointing. The observations will provide radio continuum, spectral line and polarisation information. As such, MIGHTEE, along with the excellent multi-wavelength data already available in these deep fields, will allow a range of science to be achieved. Specifically, MIGHTEE is designed to significantly enhance our understanding of, (i) the evolution of AGN and star-formation activity over cosmic time, as a function of stellar mass and environment, free of dust obscuration; (ii) the evolution of neutral hydrogen in the Universe and how this neutral gas eventually turns into stars after moving through the molecular phase, and how efficiently this can fuel AGN activity; (iii) the properties of cosmic magnetic fields and how they evolve in clusters, filaments and galaxies. MIGHTEE will reach similar depth to the planned SKA all-sky survey, and thus will provide a pilot to the cosmology experiments that will be carried out by the SKA over a much larger survey volume.

The SKA era is set to revolutionize our understanding of neutral hydrogen (HI) in individual galaxies out to redshifts of z~0.8; and in the z > 6 intergalactic medium through the detection and imaging of cosmic reionization. Direct HI number density constraints will, nonetheless, remain relatively weak out to cosmic noon (z~2) - the epoch of peak star formation and black hole accretion - and beyond. However, as was demonstrated from the 1990s with molecular line observations, this can be overcome by utilising the natural amplification afforded by strong gravitational lensing, which results in an effective increase in integration time by the square of the total magnification (\mu^2) for an unresolved source. Here we outline how a dedicated lensed HI survey will leverage MeerKAT's high sensitivity, frequency coverage, large instantaneous bandwidth, and high dynamic range imaging to enable a lasting legacy of high-redshift HI emission detections well into the SKA era. This survey will not only provide high-impact, rapid-turnaround MeerKAT science commissioning results, but also unveil Milky Way-like systems towards cosmic noon which is not possible with any other SKA precursors/pathfinders. An ambitious lensed HI survey will therefore make a significant impact from MeerKAT commissioning all the way through to the full SKA era, and provide a more complete picture of the HI history of the Universe.

Einstein's General Theory of Relativity (GR) successfully describes gravity. The most fundamental predictions of GR are black holes (BHs), but in spite of many convincing BH candidates in the Universe, there is no conclusive experimental proof of their existence using astronomical observations in the electromagnetic spectrum. Are BHs real astrophysical objects? Does GR hold in its most extreme limit or are alternatives needed? The prime target to address these fundamental questions is in the center of our own Galaxy, which hosts the closest and best-constrained supermassive BH candidate in the Universe, Sagittarius A* (Sgr A*). Three different types of experiments hold the promise to test GR in a strong-field regime using observations of Sgr A* with new-generation instruments. The first experiment aims to image the relativistic plasma emission which surrounds the event horizon and forms a "shadow" cast against the background, whose predicted size (~50 microarcseconds) can now be resolved by upcoming VLBI experiments at mm-waves such as the Event Horizon Telescope (EHT). The second experiment aims to monitor stars orbiting Sgr A* with the upcoming near-infrared interferometer GRAVITY at the Very Large Telescope (VLT). The third experiment aims to time a radio pulsar in tight orbit about Sgr A* using radio telescopes (including the Atacama Large Millimeter Array or ALMA). The BlackHoleCam project exploits the synergy between these three different techniques and aims to measure the main BH parameters with sufficient precision to provide fundamental tests of GR and probe the spacetime around a BH in any metric theory of gravity. Here, we review our current knowledge of the physical properties of Sgr A* as well as the current status of such experimental efforts towards imaging the event horizon, measuring stellar orbits, and timing pulsars around Sgr A*.

The Event Horizon Telescope (EHT) aims to spatially resolve the silhouette (or shadow) of the supermassive black holes in the Galactic Centre (Sgr A$^\star$) and M87. The primary scientific objectives are to test general relativity in the strong-field regime and to probe accretion and jet-launch physics at event-horizon scales. This is made possible by the technique of Very Long Baseline Interferometry (VLBI) at (sub)millimetre wavelengths, which can achieve angular resolutions of order $\sim10~\mu$-arcsec. However, this approach suffers from unique observational challenges, including scattering in the troposphere and interstellar medium; rapidly time-variable source structure in both polarized and total intensity; as well as non-negligible antenna pointing errors. In this, the first paper in a series, we present the MeqSilhouette software package which is specifically designed to accurately simulate EHT observations. It includes realistic descriptions of a number of signal corruptions that can limit the ability to measure the key science parameters. This can be used to quantify calibration requirements, test parameter estimation and imaging strategies, and investigate systematic uncertainties that may be present. In doing so, a stronger link can be made between observational capabilities and theoretical predictions, with the goal of maximising scientific output from the upcoming order of magnitude increase in EHT sensitivity.

Strong gravitational lensing provides some of the deepest views of the Universe, enabling studies of high-redshift galaxies only possible with next-generation facilities without the lensing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected directly out to z~0.2, limited by the sensitivity and instantaneous bandwidth of current radio telescopes. We discuss how current and future radio interferometers such as the Square Kilometre Array (SKA) will detect lensed HI emission in individual galaxies at high redshift. Our calculations rely on a semi-analytic galaxy simulation with realistic HI disks (by size, density profile and rotation), in a cosmological context, combined with general relativistic ray tracing. Wide-field, blind HI surveys with the SKA are predicted to be efficient at discovering lensed HI systems, increasingly so at z > 2. This will be enabled by the combination of the magnification boosts, the steepness of the HI luminosity function at the high-mass end, and the fact that the HI spectral line is relatively isolated in frequency. These surveys will simultaneously provide a new technique for foreground lens selection and yield the highest redshift HI emission detections. More near term (and existing) cm-wave facilities will push the high redshift HI envelope through targeted surveys of known lenses.

Galaxies and supermassive black holes (SMBHs) are believed to evolve through a process of hierarchical merging and accretion. Through this paradigm, multiple SMBH systems are expected to be relatively common in the Universe. However, to date there are poor observational constraints on multiple SMBHs systems with separations comparable to a SMBH gravitational sphere of influence (<< 1 kpc). In this chapter, we discuss how deep continuum observations with the SKA will make leading contributions towards understanding how multiple black hole systems impact galaxy evolution. In addition, these observations will provide constraints on and an understanding of stochastic gravitational wave background detections in the pulsar timing array sensitivity band (nanoHz -microHz). We also discuss how targets for pointed gravitational wave experiments (that cannot be resolved by VLBI) could potentially be found using the large-scale radio-jet morphology, which can be modulated by the presence of a close-pair binary SMBH system. The combination of direct imaging at high angular resolution; low-surface brightness radio-jet tracers; and pulsar timing arrays will allow the SKA to trace black hole binary evolution from separations of a galaxy virial radius down to the sub-parsec level. This large dynamic range in binary SMBH separation will ensure that the SKA plays a leading role in this observational frontier.

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.

Galaxies are believed to evolve through merging, which should lead to multiple supermassive black holes in some. There are four known triple black hole systems, with the closest pair being 2.4 kiloparsecs apart (the third component is more distant at 3 kiloparsecs), which is far from the gravitational sphere of influence of a black hole with mass $\sim$10$^9$ M$_\odot$ (about 100 parsecs). Previous searches for compact black hole systems concluded that they were rare, with the tightest binary system having a separation of 7 parsecs. Here we report observations of a triple black hole system at redshift z=0.39, with the closest pair separated by $\sim$140 parsecs. The presence of the tight pair is imprinted onto the properties of the large-scale radio jets, as a rotationally-symmetric helical modulation, which provides a useful way to search for other tight pairs without needing extremely high resolution observations. As we found this tight pair after searching only six galaxies, we conclude that tight pairs are more common than hitherto believed, which is an important observational constraint for low-frequency gravitational wave experiments.

We measure star-formation rates (SFRs) and specific SFRs (SSFRs) of Ks-selected galaxies from the VIDEO survey by stacking 1.4-GHz Very Large Array data. We split the sample, which spans 0 < z < 3 and stellar masses 10**8.0 < Mstellar/Msol < 10**11.5, into elliptical, irregular or starburst galaxies based on their spectral-energy distributions. We find that SSFR falls with stellar mass, in agreement with the `downsizing' paradigm. We consider the dependence of the SSFR-mass slope on redshift: for our full and elliptical samples the slope flattens, but for the irregular and starburst samples the slope is independent of redshift. The rate of SSFR evolution reduces slightly with stellar mass for ellipticals, but irregulars and starbursts co-evolve across stellar masses. Our results for SSFR as a function of stellar mass and redshift are in agreement with those derived from other radio-stacking measurements of mass-selected passive and star-forming galaxies, but inconsistent with those generated from semi-analytic models, which tend to underestimate SFRs and SSFRs. There is a need for deeper high-resolution radio surveys such as those from telescopes like the next-generation MeerKAT in order to probe lower masses at earlier times and to permit direct detections, i.e. to study individual galaxies in detail.

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.

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.

The Arcminute Microkelvin Imager (AMI) is a telescope specifically designed for high sensitivity measurements of low-surface-brightness features at cm-wavelength and has unique, important capabilities. It consists of two interferometer arrays operating over 13.5-18 GHz that image structures on scales of 0.5-10 arcmin with very low systematics. The Small Array (AMI-SA; ten 3.7-m antennas) couples very well to Sunyaev-Zel'dovich features from galaxy clusters and to many Galactic features. The Large Array (AMI-LA; eight 13-m antennas) has a collecting area ten times that of the AMI-SA and longer baselines, crucially allowing the removal of the effects of confusing radio point sources from regions of low surface-brightness, extended emission. Moreover AMI provides fast, deep object surveying and allows monitoring of large numbers of objects. In this White Paper we review the new science - both Galactic and extragalactic - already achieved with AMI and outline the prospects for much more.

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