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

A key characteristic of many active galactic nuclei (AGN) is their variability, but its origin is poorly understood, especially in the radio domain. Williams et al. (2017) reported a ~50 per cent increase in peak flux density of the AGN in the Seyfert galaxy NGC 4151 at 1.5 GHz with the e-MERLIN array. We present new high resolution e-MERLIN observations at 5 GHz and compare these to archival MERLIN observations to investigate the reported variability. Our new observations allow us to probe the nuclear region at a factor three times higher-resolution than the previous e-MERLIN study. We separate the core component, C4, into three separate components: C4W, C4E and X. The AGN is thought to reside in component C4W, but this component has remained constant between epochs within uncertainties. However, we find that the Eastern-most component, C4E, has increased in peak flux density from 19.35$\pm$1.10 to 37.09$\pm$1.86 mJy/beam, representing a 8.2 sigma increase on the MERLIN observations. We attribute this peak flux density increase to continued interaction between the jet and the emission line region (ELR), observed for the first time in a low-luminosity AGN such as NGC 4151. We identify discrete resolved components at 5 GHz along the jet axis, which we interpret as areas of jet-ELR interaction.

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.

We report deep radio observations of nearby Type Ia Supernovae (SNe Ia) with the electronic Multi-Element Radio Linked Interferometer Net-work (e-MERLIN), and the Australia Telescope Compact Array (ATCA). No detections were made. With standard assumptions for the energy densities of relativistic electrons going into a power-law energy distribution, and the magnetic field strength (epsilon_e = epsilon_B = 0.1), we arrive at the upper limits on mass-loss rate for the progenitor system of SN 2013dy (2016coj, 2018gv, 2018pv, 2019np), to be less than 12 (2.8,1.3, 2.1, 1.7)EE(-8) solar masses per year (for a wind speed of 100 km/s). To SNe 2016coj, 2018gv, 2018pv and 2019np we add radio data for 17 other nearby SNe Ia, and model their non-detections. With the same model as described, all 21 SNe Ia have mass-loss rates less than 4EE(-8) solar masses per year (for a wind speed of 100 km/s). We compare those limits with the expected mass loss rates in different single-degenerate progenitor scenarios. We also discuss how information on epsilon_e and epsilon_B can be obtained from late observations of SNe Ia and the youngest SN Ia remnant detected in radio, G1.9+0.3, as well as stripped-envelope core-collapse SNe. We highlight SN 2011dh, and argue for epsilon_e approximately equal to 0.1 and epsilon_B approximately equal to 0.0033. Finally, we discuss strategies to observe at radio frequencies to maximize the chance of detection, given the time since explosion, the distance to the supernova and the telescope sensitivity.

We report on VLBI observations of the fast and blue optical transient (FBOT), AT 2018cow. At ~62 Mpc, AT 2018cow is the first relatively nearby FBOT. The nature of AT 2018cow is not clear, although various hypotheses from a tidal disruption event to different kinds of supernovae have been suggested. It had a very fast rise time (3.5 d) and an almost featureless blue spectrum although high photospheric velocities (40,000 km s$^{-1}$) were suggested early on. The X-ray luminosity was very high, ~$1.4 \times 10^{43}$ erg s$^{-1}$, larger than those of ordinary SNe, and more consistent with those of SNe associated with gamma-ray bursts. Variable hard X-ray emission hints at a long-lived "central engine." It was also fairly radio luminous, with a peak 8.4-GHz spectral luminosity of ~$4 \times 10^{28}$ erg s$^{-1}$ Hz$^{-1}$, allowing us to make VLBI observations at ages between 22 and 287 d. We do not resolve AT 2018cow. Assuming a circularly symmetric source, our observations constrain the average apparent expansion velocity to be <0.49$c$ by t = 98 d (3$\sigma$ limit). We also constrain the proper motion of AT 2018cow to be <0.51$c$. Since the radio emission generally traces the fastest ejecta, our observations make the presence of a long-lived relativistic jet with a lifetime of more than one month very unlikely.

We present high-sensitivity 1.51 GHz e-MERLIN radio images of the nearby galaxy NGC 6217. We resolve the compact radio source at the centre of NGC 6217 for the first time, finding a twin-lobed structure with a total linear size of ~4 arcsec (~400pc). The radio source does not have a compact central core but there is an unresolved hot spot near the outer edge of the southern lobe. Combining our e-MERLIN data with new VLA A-Array data and archival multi-wavelength data, we explore possible scenarios which might explain this complex radio morphology. We conclude that this radio source is most likely powered by a low-luminosity AGN (LLAGN) but with a possible important contribution from nuclear star formation. We also investigate the origin of a putative X-ray jet in NGC 6217, previously suggested in the literature. This `jet' consists of three X-ray `knots' in a line, pointed away from the nucleus, with a total size of ~3 arcmin (~15 kpc). We find no extended radio emission coincident with the `jet'. An additional compact radio source, aligned with the knots but without an X-ray counterpart is detected. We detect IR/optical sources falling within the X-ray extraction regions of the `knots', but note that these sources could be chance associations due to high source density around the target field and we estimate the probability of three randomly aligned X-ray sources to be 0.3 per cent in this field.

Intermediate-mass black holes are theoretically predicted but observationally elusive, and evidence for them is often indirect. The nearby face-on spiral galaxy NGC3310 has hosted many supernovae in recent history, and recent Chandra observations have shown a group of strong off-nuclear X-ray sources that are coincident with radio emission seen in archival VLA and MERLIN observations. Their luminosity, spectrum and off-nuclear location make these sources excellent IMBH candidates. To investigate this possibility, we used combined EVN/e-MERLIN observations at both 1.4 and 5 GHz to look for compact radio emission and evidence of jet activity. We detect a compact radio source within one arcsecond of a Chandra source with an estimated mass ${\rm M}_{\rm BH}\sim3\times10^4 {\rm M}_{\odot}$.

We studied the nearby, interacting galaxy NGC 5195 (M51b) in the radio, optical and X-ray bands. We mapped the extended, low-surface-brightness features of its radio-continuum emission; determined the energy content of its complex structure of shock-ionized gas; constrained the current activity level of its supermassive nuclear black hole. In particular, we combined data from the European Very Long Baseline Interferometry Network (~1-pc scale), from our new e-MERLIN observations (~10-pc scale), and from the Very Large Array (~100-1000-pc scale), to obtain a global picture of energy injection in this galaxy. We put an upper limit to the luminosity of the (undetected) flat-spectrum radio core. We find steep-spectrum, extended emission within 10 pc of the nuclear position, consistent with optically-thin synchrotron emission from nuclear star formation or from an outflow powered by an active galactic nucleus (AGN). A linear spur of radio emission juts out of the nuclear source towards the kpc-scale arcs (detected in radio, Halpha and X-ray bands). From the size, shock velocity, and Balmer line luminosity of the kpc-scale bubble, we estimate that it was inflated by a long-term-average mechanical power ~3-6 x 10^{41} erg/s over the last 3-6 Myr. This is an order of magnitude more power than can be provided by the current level of star formation, and by the current accretion power of the supermassive black hole. We argue that a jet-inflated bubble scenario associated with previous episodes of AGN activity is the most likely explanation for the kpc-scale structures.

We present the first data release of high-resolution ($\leq0.2$ arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the LeMMINGs survey. This sample includes galaxies which are active (LINER and Seyfert) and quiescent (HII galaxies and Absorption line galaxies, ALG), which are reclassified based upon revised emission-line diagrams. We detect radio emission $\gtrsim$ 0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for HII galaxies and 5/14 for ALGs) with radio sizes typically of $\lesssim$100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 10$^{32}$ to 10$^{40}$ erg s$^{-1}$: LINERs and HII galaxies show the highest and the lowest radio powers respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to $\sim$10$^{7}$ M$_{\odot}$, but a break emerges at lower masses. Using [O III] line luminosity as a proxy for the accretion luminosity, active nuclei and jetted HII galaxies follow an optical fundamental plane of BH activity, suggesting a common disc-jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; HII galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties.

We have analysed a new high-resolution e-MERLIN 1.5 GHz radio continuum map together with $HST$ and SDSS imaging of NGC 5322, an elliptical galaxy hosting radio jets, aiming to understand the galaxy's central structure and its connection to the nuclear activity. We decomposed the composite $HST$ + SDSS surface brightness profile of the galaxy into an inner stellar disc, a spheroid, and an outer stellar halo. Past works showed that this embedded disc counter-rotates rapidly with respect to the spheroid. The $HST$ images reveal an edge-on nuclear dust disc across the centre, aligned along the major-axis of the galaxy and nearly perpendicular to the radio jets. After careful masking of this dust disc, we find a central stellar mass deficit $M_{\rm def}$ in the spheroid, scoured by SMBH binaries with final mass $M_{\rm BH}$ such that $M_{\rm def}/M_{\rm BH} \sim 1.3 - 3.4$. We propose a three-phase formation scenario for NGC 5322 where a few ($2-7$) "dry" major mergers involving SMBHs built the spheroid with a depleted core. The cannibalism of a gas-rich satellite subsequently creates the faint counter-rotating disc and funnels gaseous material directly onto the AGN, powering the radio core with a brightness temperature of $T_{\rm B,core} \sim 4.5 \times 10^{7}$ K and the low-power radio jets ($P_{\rm jets}\sim 7.04 \times 10^{20}$ W Hz$^{-1}$) which extend $\sim 1.6$ kpc. The outer halo can later grow via minor mergers and the accretion of tidal debris. The low-luminosity AGN/jet-driven feedback may have quenched the late-time nuclear star formation promptly, which could otherwise have replenished the depleted core.

V959 Mon (Nova Mon 2012) was first detected by the Fermi Large Area Telescope in June 2012, as a transient gamma-ray source. Subsequent optical observations showed that this gamma-ray emission was due to a classical nova explosion. Multi-frequency observations of V959 Mon with the VLA between June and September 2012 revealed dramatic brightening, and a spectrum that steepened with increasing frequency. High resolution radio images of V959 Mon using e-MERLIN are presented here, at six epochs between September 2012 and February 2014 which show morphological evolution of the source. While early e-MERLIN observations of V959 Mon show an east-west elonga- tion in the ejecta morphology, subsequent observations suggest that the ejecta become elongated in the north-south direction. Our high-resolution observations of this sur- prising evolution in the structure of V959 Mon can assist us in further understanding the behaviour and morphology of nova ejecta.

We present high-sensitivity eMERLIN radio images of the Seyfert galaxy NGC 4151 at 1.5 GHz. We compare the new eMERLIN images to those from archival MERLIN observations in 1993 to determine the change in jet morphology in the 22 years between observations. We report an increase by almost a factor of 2 in the peak flux density of the central core component, C4, thought to host the black hole, but a probable decrease in some other components, possibly due to adiabatic expansion. The core flux increase indicates an AGN which is currently active and feeding the jet. We detect no significant motion in 22 years between C4 and the component C3, which is unresolved in the eMERLIN image. We present a spectral index image made within the 512 MHz band of the 1.5 GHz observations. The spectrum of the core, C4, is flatter than that of other components further out in the jet. We use HST emission line images (H$\alpha$, [O III] and [O II]) to study the connection between the jet and the emission line region. Based on the changing emission line ratios away from the core and comparison with the eMERLIN radio jet, we conclude that photoionisation from the central AGN is responsible for the observed emission line properties further than 4" (360 pc) from the core, C4. Within this region, several evidences (radio-line co-spatiality, low [O III]/H$\alpha$ and estimated fast shocks) suggest additional ionisation from the jet.

We present new observations of the GOODS-N field obtained at 5.5 GHz with the Karl G. Jansky Very Large Array (VLA). The central region of the field was imaged to a median r.m.s. of 3 microJy/beam with a resolution of 0.5 arcsec. From a 14-arcmin diameter region we extracted a sample of 94 radio sources with signal-to-noise ratio greater than 5. Near-IR identifications are available for about 88 percent of the radio sources. We used different multi-band diagnostics to separate active galactic nuclei (AGN), both radiatively efficient and inefficient, from star-forming galaxies. From our analysis, we find that about 80 percent of our radio-selected sample is AGN-dominated, with the fraction raising to 92 percent when considering only the radio sources with redshift >1.5. This large fraction of AGN-dominated radio sources at very low flux densities (the median flux density at 5.5 GHz is 42 microJy), where star-forming galaxies are expected to dominate, is somewhat surprising and at odds with other results. Our interpretation is that both the frequency and angular resolution of our radio observations strongly select against radio sources whose brightness distribution is diffuse on scale of several kpc. Indeed, we find that the median angular sizes of the AGN-dominated sources is around 0.2-0.3 arcsec against 0.8 arcsec for star-forming galaxies. This highlights the key role that high frequency radio observations can play in pinpointing AGN-driven radio emission at microJy levels. This work is part of the eMERGE legacy project.

IC883 is a luminous infrared galaxy (LIRG) classified as a starburst-active galactic nucleus (AGN) composite. In a previous study we detected a low-luminosity AGN (LLAGN) radio candidate. Here we report on our radio follow-up at three frequencies which provides direct and unequivocal evidence of the AGN activity in IC883. Our analysis of archival X-ray data, together with the detection of a transient radio source with luminosity typical of bright supernovae, give further evidence of the ongoing star formation activity, which dominates the energetics of the system. At sub-parsec scales, the radio nucleus has a core-jet morphology with the jet being a newly ejected component showing a subluminal proper motion of 0.6c-1c. The AGN contributes less than two per cent of the total IR luminosity of the system. The corresponding Eddington factor is ~1E-3, suggesting this is a low-accretion rate engine, as often found in LLAGNs. However, its high bolometric luminosity (~10E44erg/s) agrees better with a normal AGN. This apparent discrepancy may just be an indication of the transition nature of the nucleus from a system dominated by star-formation, to an AGN-dominated system. The nucleus has a strongly inverted spectrum and a turnover at ~4.4GHz, thus qualifying as a candidate for the least luminous (L_5.0GHz ~ 6.3E28erg/s/Hz) and one of the youngest (~3000yr) gigahertz-peaked spectrum (GPS) sources. If the GPS origin for the IC883 nucleus is confirmed, then advanced mergers in the LIRG category are potentially key environments to unveil the evolution of GPS sources into more powerful radio galaxies.

We present high-resolution e-MERLIN radio continuum maps of the Dwarf Irregular galaxy IC10 at 1.5GHz and 5GHz. We detect 11 compact sources at 1.5GHz, 5 of which have complementary detections at 5GHz. We classify 3 extended sources as compact H{\sc ii} regions within IC10, 5 sources as contaminating background galaxies and identify 3 sources which require additional observations to classify. We do not expect that any of these 3 sources are Supernova Remnants as they will likely be resolved out at the assumed distance of IC10 (0.7Mpc). We correct integrated flux densities of IC10 from the literature for contamination by unrelated background sources and obtain updated flux density measurements of $354\pm11$\,mJy at 1.5GHz and $199\pm9$\,mJy at 4.85GHz. The background contamination does not contribute significantly to the overall radio emission from IC10, so previous analysis concerning its integrated radio properties remain valid.

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN.

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary.

The core of the nearby galaxy NGC660 has recently undergone a spectacular radio outburst; using a combination of archival radio and Chandra X-ray data, together with new observations, the nature of this event is investigated. Radio observations made using e-MERLIN in mid-2013 show a new compact and extremely bright continuum source at the centre of the galaxy. High angular resolution observations carried out with the European VLBI Network show an obvious jet-like feature to the north east and evidence of a weak extension to the west, possibly a counter-jet. We also examine high angular resolution HI spectra of these new sources, and the radio spectral energy distribution using the new wide-band capabilities of e-MERLIN. We compare the properties of the new object with possible explanations, concluding that we are seeing a period of new AGN activity in the core of this polar ring galaxy.

We present the highest resolution, wide-field radio survey of a nearby face-on star-forming galaxy to date. The multi-phase centre technique is used to survey the entire disk of M51 (77 square arc minutes) at a maximum resolution of 5 milli-arcseconds on a single 8 hr pointing with the European VLBI Network at 18 cm. In total, 7 billion pixels were imaged using 192 phase centres that resulted in the detection of six sources: the Seyfert nucleus, the supernova SN 2011dh, and four background AGNs. Using the wealth of archival data available in the radio (MERLIN and the VLA), optical (Hubble Space Telescope) and X-rays (Chandra) the properties of the individual sources were investigated in detail. The combined multi-wavelength observations reveal a very complex and puzzling core region that includes a low-luminosity parsec scale core-jet structure typical of AGNs, with evidence for a lateral shift corresponding to 0.27c. Furthermore, there is evidence for a fossil radio hotspot located 1.44 kpc from the Seyfert nucleus that may have resulted from a previous ejection cycle. Our study provides measures of the supernova and star-formation rates that are comparable to independent studies at other wavelengths, and places further limits on the radio and X-ray luminosity evolution of the supernovae SN 1994I, SN 2005cs and SN 2011dh. The radio images of background AGN reveal complex morphologies that are indicative of powerful radio galaxies, and confirmed via the X-ray and optical properties.

Written in Python and utilising ParselTongue to interface with the Astronomical Image Processing System (AIPS), the e-MERLIN data reduction pipeline is intended to automate the procedures required in processing and calibrating radio astronomy data from the e-MERLIN correlator. Driven by a plain text file of input parameters, the pipeline is modular and can be run in stages by the user, depending on requirements. The software includes options to load raw data, average in time and/or frequency, flag known sources of interference, flag more comprehensively with SERPent, carry out some or all of the calibration procedures including self-calibration), and image in either normal or wide-field mode. It also optionally produces a number of useful diagnostic plots at various stages so that the quality of the data can be assessed. The software is available for download from the e-MERLIN website or via Github.

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $\mu$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

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

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR.

Prior to the launch of NuSTAR, it was not feasible to spatially resolve the hard (E > 10 keV) emission from galaxies beyond the Local Group. The combined NuSTAR dataset, comprised of three ~165 ks observations, allows spatial characterization of the hard X-ray emission in the galaxy NGC 253 for the first time. As a follow up to our initial study of its nuclear region, we present the first results concerning the full galaxy from simultaneous NuSTAR, Chandra, and VLBA monitoring of the local starburst galaxy NGC 253. Above ~10 keV, nearly all the emission is concentrated within 100" of the galactic center, produced almost exclusively by three nuclear sources, an off-nuclear ultraluminous X-ray source (ULX), and a pulsar candidate that we identify for the first time in these observations. We detect 21 distinct sources in energy bands up to 25 keV, mostly consisting of intermediate state black hole X-ray binaries. The global X-ray emission of the galaxy - dominated by the off-nuclear ULX and nuclear sources, which are also likely ULXs - falls steeply (photon index >~ 3) above 10 keV, consistent with other NuSTAR-observed ULXs, and no significant excess above the background is detected at E > 40 keV. We report upper limits on diffuse inverse Compton emission for a range of spatial models. For the most extended morphologies considered, these hard X-ray constraints disfavor a dominant inverse Compton component to explain the {\gamma}-ray emission detected with Fermi and H.E.S.S. If NGC 253 is typical of starburst galaxies at higher redshift, their contribution to the E > 10 keV cosmic X-ray background is < 1%.

We report deep EVN and eMERLIN observations of the Type Ia SN 2014J in the nearby galaxy M 82. Our observations represent, together with JVLA observations of SNe 2011fe and 2014J, the most sensitive radio studies of Type Ia SNe ever. By combining data and a proper modeling of the radio emission, we constrain the mass-loss rate from the progenitor system of SN 2014J to $\dot{M} \lesssim 7.0\times 10^{-10}\, {\rm M_{\odot}\, yr^{-1}}$ (3-$\sigma$; for a wind speed of $100\, {\rm km s^{-1}}$). If the medium around the supernova is uniform, then $n_{\rm ISM} \lesssim 1.3 {\rm cm^3}$ (3-$\sigma$), which is the most stringent limit for the (uniform) density around a Type Ia SN. Our deep upper limits favor a double-degenerate (DD) scenario--involving two WD stars--for the progenitor system of SN 2014J, as such systems have less circumstellar gas than our upper limits. By contrast, most single-degenerate (SD) scenarios, i.e., the wide family of progenitor systems where a red giant, main-sequence, or sub-giant star donates mass to a exploding WD, are ruled out by our observations. Our estimates on the limits to the gas density surrounding SN 2011fe, using the flux density limits from Chomiuk et al. (2012), agree well with their results. Although we discuss possibilities for a SD scenario to pass observational tests, as well as uncertainties in the modeling of the radio emission, the evidence from SNe 2011fe and 2014J points in the direction of a DD scenario for both.

We present results from three nearly simultaneous NuSTAR and Chandra monitoring observations between 2012 Sep 2 and 2012 Nov 16 of local star-forming galaxy NGC 253. The 3-40 keV NuSTAR intensity of the inner 20 arcsec (~400 pc) nuclear region varied by a factor of ~2 across the three monitoring observations. The Chandra data reveal that the nuclear region contains three bright X-ray sources, including a luminous (L2-10 keV ~ few x 10^39 erg/s) point source ~1 arcsec from the dynamical center of the galaxy (within the 3sigma positional uncertainty of the dynamical center); this source drives the overall variability of the nuclear region at energies >3 keV. We make use of the variability to measure the spectra of this single hard X-ray source when it was in bright states. The spectra are well described by an absorbed (NH ~ 1.6 x 10^23 cm^-2) broken power-law model with spectral slopes and break energies that are typical of ultraluminous X-ray sources (ULXs), but not AGN. A previous Chandra observation in 2003 showed a hard X-ray point source of similar luminosity to the 2012 source that was also near the dynamical center (~0.4 arcsec); however, this source was offset from the 2012 source position by ~1 arcsec. We show that the probability of the 2003 and 2012 hard X-ray sources being unrelated is >>99.99% based on the Chandra spatial localizations. Interestingly, the Chandra spectrum of the 2003 source (3-8 keV) is shallower in slope than that of the 2012 hard X-ray source. Its proximity to the dynamical center and harder Chandra spectrum indicate that the 2003 source is a better AGN candidate than any of the sources detected in our 2012 campaign; however, we were unable to rule out a ULX nature for this source. Future NuSTAR and Chandra monitoring would be well equipped to break the degeneracy between the AGN and ULX nature of the 2003 source, if again caught in a high state.

The VLBA was used at 1.6 GHz to observe a target field 50 arcminutes in diameter including the core of M31. Novel VLBI correlation techniques were used to observe 200 sources simultaneously, of which 16 were detected. We classify all 16 as background AGN based on their X-ray properties and arcsecond- and mas-scale morphology. The detected sources were then analyzed for evidence of scatter-broadening due to the ionized ISM of M31. The detection of a compact background source only 0.25 kpc projected distance from M31 places a constraint on the extent of any extreme scattering region associated with center of M31. However, the two sources closest to the core show evidence of scatter broadening consistent with that which would be observed if a compact source were observed through the inner disk of our Galaxy at the inclination of M31. We interpret this as a detection of the ionized ISM of M31 along two lines of sight. With the increases in bandwidth and sensitivity envisaged for future long-baseline interferometers, this should prove to be a remarkably powerful technique for understanding the ionized ISM in external galaxies.

We present first results from electronic Multi-Element Remotely Linked Interferometer Network (e-MERLIN) and electronic European VLBI Network (e-EVN) observations of a small sample of ultra-steep spectrum radio sources, defined as those sources with a spectral index alpha < -1.4 between 74 MHz and 325 MHz, which are unresolved on arcsecond scales. Such sources are currently poorly understood and a number of theories as to their origin have been proposed in the literature. The new observations described here have resulted in the first detection of two of these sources at milliarcsecond scales and show that a significant fraction of ultra-steep spectrum sources may have compact structures which can only be studied at the high resolution available with very long baseline interferometry (VLBI).

Despite first being detected in the 1970s, surprisingly little is known about the OH main line maser population in the nearby starburst galaxy M82. Sometimes referred to as 'kilomasers', they have isotropic luminosities intermediate between Galactic masers and those found in more distant megamasers. Several observations have been carried by this group over the last ten years in an attempt to get a better handle on their nature. High velocity resolution VLA observations in 2006 showed that almost all of the maser spots, distributed across the central arcminute of the galaxy, were apparently coincident with background continuum features, and a handful displayed multiple velocity components. The majority of those with velocity structure are located on a blue-shifted arc in the pv-plane, spatially located on an arc northward of the peculiar source known as B41.95+57.5. Now, new results from high spatial and spectral resolution observations with the EVN have resolved several of these masers into multiple spatial components for the first time. The maser emission is compared with known continuum sources in the galaxy, and we conclude that at least some of the maser emission is from high-gain maser action.

This paper presents the results of the 2009-2010 monitoring sessions of the starburst galaxy M82, obtained with the Multi-Element Radio-Linked Interferometer Network (MERLIN) at 5GHz and e-MERLIN at 6GHz. Combining several 5GHz MERLIN epochs to form a map with 33.0 uJy/bm noise level, 52 discrete sources, mostly supernova remnants and HII regions, are identified. These include three objects which were not detected in the 2002 5GHz MERLIN monitoring session: supernova SN2008iz, the transient source 43.78+59.3, and a new supernova remnant shell. Flux density variations, in the long (1981 to 2010), medium (2002 to 2010) and short (2009 to 2010) term, are investigated. We find that flux densities of SNRs in M82 stay constant in most of the sample (~95%). In addition, aside from SN2008iz and the well-known variable source 41.95+57.5, two sources display short and medium term variations over the period 2009-2010. These sources being among the most compact SNR in M82, these flux density variations could be due to changes in the circumstellar and interstellar medium in which the shocks travel.

The high star formation rates of luminous infrared galaxies (LIRGs) make them ideal places for core-collapse supernova (CCSN) searches. At radio frequencies, free from dust extinction, it is possible to detect compact components within the innermost LIRG nuclear regions, such as SNe and SN remnants, as well as AGN buried deep in the LIRG nuclei. We studied the LIRG IC883 aiming at: (i) investigating its (circum-)nuclear regions using the e-EVN at 5GHz, and e-MERLIN at 6.9GHz, complemented by archival VLBI data; (ii) detecting at radio frequencies the two recently reported circumnuclear SNe 2010cu and 2011hi, which were discovered by near-IR (NIR) adaptive optics observations of IC883; and (iii) further investigating the nature of SN2011hi at NIR by means of observations with Gemini-North. The circumnuclear regions traced by e-MERLIN at 6.9GHz have an extension of ~1kpc, and show a striking double-sided structure, which very likely corresponds to a warped rotating ring, in agreement with previous studies. Our e-EVN observations at 5GHz and complementary archival VLBI data at 5GHz and 8.4GHz, reveal the presence of various milliarcsec compact components in the nucleus of IC883. A single compact source, an AGN candidate, dominates the emission at both nuclear and circumnuclear scales, as imaged with the e-EVN and e-MERLIN, respectively. The other milliarcsec components are very suggestive of ongoing nuclear CCSN activity. Our e-EVN observations also resulted in upper limits to the radio luminosity of the two SNe in IC883 recently discovered at NIR. We refine the classification of SN2011hi as a Type IIP SN according to our latest Gemini-North epoch from 2012, in agreement with a low-luminosity radio SN nature. We estimate a CCSN rate lower limit of 1.1_{-0.6}^{+1.3} yr^{-1} for the entire galaxy, based on three nuclear radio SNe and the circumnuclear SNe 2010cu and 2011hi. (abridged)

We report on the VLBI detection of supernova SN2011dh at 22GHz using a subset of the EVN array. The observations took place 14 days after the discovery of the supernova, thus resulting in a VLBI image of the youngest radio-loud supernova ever. We provide revised coordinates for the supernova with milli-arcsecond precision, linked to the ICRF. The recovered flux density is a factor 2 below the EVLA flux density reported by other authors at the same frequency and epoch of our observations. This discrepancy could be due to extended emission detected with the EVLA or to calibration problems in the VLBI and/or EVLA observations.

One of our closest neighbours, the Andromeda Galaxy (M31) has been the subject of numerous large area studies across the entire spectrum, but so far full-disk radio surveys have been conducted only at low resolution. The new wide-field capabilities of the DiFX software correlator present the possibility of imaging the entire primary beam of a VLBI array, thus enabling a high resolution wide-field study of the entire galaxy. Using the VLBA and EVN, pilot observations of M31 have been carried out with the aim of using these new wide-field techniques to characterise the population of compact components at VLBI resolution both within and behind one of our nearest neighbours. This contribution describes the observations carried out, the preliminary processing and first results.

From a combination of MERLIN (Multi-Element Radio-Linked Interferometer Network) and global VLBI (Very Long Baseline Interferometry) observations of the starburst galaxy M82, images of 36 discrete sources at resolutions ranging from ~3 to ~80 mas at 1.7 GHz are presented. Of these 36 sources, 32 are identified as supernova remnants, 2 are HII regions, and 3 remain unclassified. Sizes, flux densities and radio brightnesses are given for all of the detected sources. Additionally, global VLBI only data from this project are used to image four of the most compact radio sources. These data provide a fifth epoch of VLBI observations of these sources, covering a 19-yr time-line. In particular, the continued expansion of one of the youngest supernova remnants, 43.31+59.3 is discussed. The deceleration parameter is a power-law index used to represent the time evolution of the size of a supernova remnant. For the source 43.31+59.3, a lower limit to the deceleration parameter is calculated to be 0.53+/-0.06, based on a lower limit of the age of this source.

A faint new radio source has been detected in the nuclear region of the starburst galaxy M82 using MERLIN radio observations designed to monitor the flux density evolution of the recent bright supernova SN2008iz. This new source was initially identified in observations made between 1-5th May 2009 but had not been present in observations made one week earlier, or in any previous observations of M82. In this paper we report the discovery of this new source and monitoring of its evolution over its first 9 months of existence. The true nature of this new source remains unclear, and we discuss whether this source may be an unusual and faint supernova, a supermassive blackhole associated with the nucleus of M82, or intriguingly the first detection of radio emission from an extragalactic microquasar.

The class of type Ic supernovae have drawn increasing attention since 1998 owing to their sparse association (only four so far) with long duration gamma-ray bursts. Although both phenomena originate from the core collapse of a massive star, supernovae emit mostly at optical wavelengths, whereas GRBs emit mostly in soft gamma-rays or hard X-rays. Though the GRB central engine generates ultra-relativistic jets, which beam the early emission into a narrow cone, no relativistic outflows have hitherto been found in type Ib/c supernovae explosions, despite theoretical expectations and searches. Here we report radio (interferometric) observations that reveal a mildly relativistic expansion in a nearby type Ic supernova, SN 2007gr. Using two observational epochs 60 days apart, we detect expansion of the source and establish a conservative lower limit for the average apparent expansion velocity of 0.6c. Independently, a second mildly relativistic supernova has been reported. Contrary to the radio data, optical observations of SN 2007gr indicate a typical type Ic supernova with ejecta velocities ~6000 km/s, much lower than in GRB-associated supernovae. We conclude that in SN 2007gr a small fraction of the ejecta produced a low-energy mildly relativistic bipolar radio jet, while the bulk of the ejecta were slower and, as shown by optical spectro-polarimetry, mildly aspherical.

Using the VLA, a series of high velocity resolution observations have been made of the M82 starburst at 1.6 GHz. These observations follow up on previous studies of the main line OH maser emission in the central kiloparsec of this starburst region, but with far greater velocity resolution, showing significant velocity structure in some of the maser spots for the first time. A total of thirteen masers were detected, including all but one of the previously known sources. While some of these masers are still unresolved in velocity, these new results clearly show velocity structure in spectra from several of the maser regions. Position-velocity plots show good agreement with the distribution of H{\sc i} including interesting velocity structure on the blueward feature in the west of the starburst which traces the velocity distribution seen in the ionised gas.

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type IIb supernova SN 2008ax, made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal detection of the supernova. The total flux density recovered from our VLBI image is 0.8$\pm$0.3 mJy (one standard deviation). As it appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48$\pm$0.12 mJy, compatible with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion velocities derived from our observations are $(1.90 \pm 0.30) \times 10^5$ km s$^{-1}$ (case of a double source) and $(5.2 \pm 1.3) \times 10^4$ km s$^{-1}$ (taking the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations.

The results of an extremely deep, 8-day long observation of the central kpc of the nearby starburst galaxy M82 using MERLIN (Multi-Element Radio Linked Interferometer Network) at 5 GHz are presented. The 17E-06 Jy/beam, rms noise level in the naturally weighted image make it the most sensitive high resolution radio image of M82 made to date. Over 50 discrete sources are detected, the majority of which are supernova remnants, but with 13 identified as HII regions. Sizes, flux densities and radio brightnesses are given for all of the detected sources, which are all well resolved with a majority showing shell or partial shell structures. Those sources within the sample which are supernova remnants have diameters ranging from 0.3 to 6.7 pc, with a mean size of 2.9 pc. From a comparison with previous MERLIN 5 GHz observations made in July 1992, which gives a 9.75 year timeline, it has been possible to measure the expansion velocities of ten of the more compact sources, eight of which have not been measured before. These derived expansion velocities range between 2200 and 10500 km/s.

A study of the distribution of OH gas in the central region of the nearby active starburst galaxy M82 has confirmed two previously known bright masers and revealed several new main line masers. Three of these are seen only at 1665 MHz, one is detected only at 1667 MHz, while the rest are detected in both lines. Observations covering both the 1665 and 1667 MHz lines, conducted with both the Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (MERLIN), have been used to accurately measure the positions and velocities of these features. This has allowed a comparison with catalogued continuum features in the starburst such as HII regions and supernova remnants, as well as known water and satellite line OH masers. Most of the main line masers appear to be associated with known HII regions although the two detected only at 1665 MHz are seen along the same line of sight as known supernova remnants.

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J.

We present the results of a second epoch of 18cm global Very Long-Baseline Interferometry (VLBI) observations, taken on 23 February 2001, of the central kiloparsec of the nearby starburst galaxy Messier 82. These observations further investigate the structural and flux evolution of the most compact radio sources in the central region of M82. The two most compact radio objects in M82 have been investigated (41.95+575 and 43.31+592). Using this recent epoch of data in comparison with our previous global VLBI observations and two earlier epochs of European VLBI Network observations we measure expansion velocities in the range of 1500-2000km/s for 41.95+575, and 9000-11000km/s for 43.31+592 using various independent methods. In each case the measured remnant expansion velocities are significantly larger than the canonical expansion velocity (500km/s) of supernova remnants within M82 predicted from theoretical models. In this paper we discuss the implications of these measured expansion velocities with respect to the high density environment that the SNR are expected to reside in within the centre of the M82 starburst.

The star formation rate (SFR) in starburst galaxies can be measured by many methods, one of which is through the supernova rate. Due to the heavy dust obscuration in these galaxies, searches for new supernovae in the optical or infra-red can easily miss events occurring in the central starburst regions. As part of a long term program to estimate the SFR in a sample of nearby starbursts we are using MERLIN and the VLA to regularly observe the galaxies for new radio supernovae. As part of this project, regular MERLIN observations have been made of two recent optically bright supernovae: 2004dj and 2004et. Both supernovae are of Type II and have been monitored frequently over periods of a few months, resulting in well sampled radio "light" curves for both objects.

Supernova 2004dj in the nearby spiral galaxy NGC2403 was detected optically in July 2004. Peaking at a magnitude of 11.2, this is the brightest supernova detected for several years. Here we present Multi-Element Radio Linked Interferometer Network (MERLIN) observations of this source, made over a four month period, which give a position of R.A. = 07h37m17.044s, Dec =+65deg35'57.84" (J2000.0). We also present a well-sampled 5 GHz light curve covering the period from 5 August to 2 December 2004. With the exception of the unusual and very close SN 1987A, these observations represent the first detailed radio light curve for the prompt emission from a Type II-P supernova.

Several new main line OH masers have been detected in the nearby starburst galaxy M82. Eight masers have been detected to 5 sigma, six of which are new detections. Observations covering both the 1665 and 1667 MHz lines with both the Very Large Array (VLA) and the Multi-Element Radio Linked Interferometer Network (MERLIN) have been used to accurately measure the positions and velocities of these features. Following analysis of the data, another six objects below 5 sigma, but with velocities consistent with the distribution inferred from the more certain detections, have been detected. These are classified as possible detections.

As part of an on-going radio supernova monitoring program, we have discovered a variable, compact steep spectrum radio source ~65 arcsec (~4 kpc) from the centre of the starburst galaxy NGC 3310. If the source is at the distance of NGC 3310, then its 5 GHz luminosity is ~3 x 10^{19} WHz^-1. The source luminosity, together with its variability characteristics, compact structure (<17 mas) and its association with a group of HII regions, leads us to propose that it is a previously uncatalogued type II radio supernova. A search of archival data also shows an associated X-ray source with a luminosity similar to known radio supernova.

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