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.

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.

The role of massive stars is central to an understanding of galactic ecology. It is important to establish the details of how massive stars provide radiative, chemical, and mechanical feedback in galaxies. Central to these issues is an understanding of the evolution of massive stars, and the critical role of mass loss via strongly structured winds and stellar binarity. Ultimately, massive stellar clusters shape the structure and energetics of galaxies. We aim to conduct high-resolution, deep field mapping at 21cm of the core of the massive Cygnus OB2 association and to characterise the properties of the massive stars and colliding winds at this waveband. We used seven stations of the e-MERLIN radio facility, with its upgraded bandwidth and enhanced sensitivity to conduct a 21cm census of Cygnus OB2. Based on 42 hours of observations, seven overlapping pointings were employed over multiple epochs during 2014 resulting in 1 sigma sensitivities down to ~21microJy and a resolution of ~180mas. A total of 61 sources are detected at 21cm over a ~0.48deg x 0.48deg region centred on the heart of the Cyg OB2 association. Of these 61 sources, 33 are detected for the first time. We detect a number of previously identified sources including four massive stellar binary systems, two YSOs, and several known X-ray and radio sources. We also detect the LBV candidate (possible binary system) and blue hypergiant (BHG) star of Cyg OB2 #12. The 21cm observations secured in the COBRaS Legacy project provide data to constrain conditions in the outer wind regions of massive stars; determine the non-thermal properties of massive interacting binaries; examine evidence for transient sources, including those associated with young stellar objects; and provide unidentified sources that merit follow-up observations. The 21cm data are of lasting value and will serve in combination with other key surveys of Cyg OB2.

Massive stars and their stellar winds are important for a number of feedback processes. The mass lost in the stellar wind can help determine the end-point of the star as a NS or a BH. However, the impact of mass-loss on the post-Main Sequence evolutionary stage of massive stars is not well understood. Westerlund 1 is an ideal astrophysical laboratory in which to study massive stars and their winds in great detail over a large range of different evolutionary phases. Aims: We aim to study the radio emission from Westerlund 1, in order to measure radio fluxes from the population of massive stars, and determine mass-loss rates and spectral indices where possible. Methods: Observations were carried out in 2015 and 2016 with the Australia telescope compact array (ATCA) at 5.5 and 9 GHz using multiple configurations, with maximum baselines ranging from 750m to 6km. Results: 30 stars were detected in the radio from the fully concatenated dataset, 10 of which were WRs (predominantly late type WN stars), 5 YHGs, 4 RSGs, 1 LBV star, the sgB[e] star W9, and several O and B supergiants. New source detections in the radio were found for 5 WR stars, and 5 OB supergiants. These detections have led to evidence for 3 new OB supergiant binary candidates, inferred from derived spectral index limits. Conclusions: Spectral indices and index limits were determined for massive stars in Westerlund 1. For cluster members found to have partially optically thick emission, mass-loss rates were calculated. Under the approximation of a thermally emitting stellar wind and a steady mass-loss rate, clumping ratios were then estimated for 8 WRs. Diffuse radio emission was detected throughout the cluster. Detections of knots of radio emission with no known stellar counterparts indicate the highly clumped structure of this intra-cluster medium, likely shaped by a dense cluster wind.

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.

We present a large-scale study of diffuse X-ray emission in the nearby massive stellar association Cygnus OB2 as part of the Chandra Cygnus OB2 Legacy Program. We used 40 Chandra X-ray ACIS-I observations covering $\sim$1.0 deg$^2$. After removing 7924 point-like sources detected in our survey, background-corrected X-ray emission, the adaptive smoothing reveals large-scale diffuse X-ray emission. Diffuse emission was detected in the sub-bands Soft [0.5 : 1.2] and Medium [1.2 : 2.5], and marginally in the Hard [2.5 : 7.0] keV band. From X-ray spectral analysis of stacked spectra we compute a total [0.5 : 7.0 keV] diffuse X-ray luminosity of L$_{\rm x}^{\rm diff}\approx$4.2$\times$10$^{\rm 34}$ erg s$^{-1}$, characterized with plasma temperature components at kT$\approx$ 0.11, 0.40 and 1.18 keV, respectively. The HI absorption column density corresponding to these temperatures has a distribution consistent with N$_{\rm H}$ = 0.43, 0.80 and 1.39 $\times$10$^{22}$ cm$^{-2}$. The extended medium band energy emission likely arises from O-type stellar winds thermalized by wind-wind collisions in the most populated regions of the association, while the soft band emission probably arises from less energetic termination shocks against the surrounding Interstellar-Medium. Super-soft and Soft diffuse emission appears more widely dispersed and intense than the medium band emission. The diffuse X-ray emission is generally spatially coincident with low-extinction regions that we attribute to the ubiquitous influence of powerful stellar winds from massive stars and their interaction with the local Interstellar-Medium. Diffuse X-ray emission is volume-filling, rather than edge-brightened, oppositely to other star-forming regions. We reveal the first observational evidence of X-ray haloes around some evolved massive stars.

Massive stars play an important role in both cluster and galactic evolution and the rate at which they lose mass is a key driver of both their own evolution and their interaction with the environment up to and including their SNe explosions. Young massive clusters provide an ideal opportunity to study a co-eval population of massive stars. We performed 3mm continuum observations with the Atacama Large Millimetre/submillimetre Array of the Galactic cluster Westerlund 1, to study the constituent massive stars and determine mass-loss rates for the diverse post-main sequence population. We detected emission from 50 stars in Westerlund 1, comprising all 21 Wolf-Rayets within the field of view, eight cool and 21 OB super-/hypergiants. Emission nebulae were associated with a number of the cool hypergiants while, unexpectedly, a number of hot stars also appear spatially resolved. We measured the mass-loss rates for a unique population of massive post-main sequence stars at every stage of evolution, confirming a significant increase as stars transition from OB supergiant to WR states. The range of spectral types exhibited provides a critical test of radiatively driven wind theory and the reality of the bi-stability jump. The extreme mass-loss rate inferred for the interacting binary Wd1-9 in comparison to other cluster members confirmed the key role binarity plays in massive stellar evolution. The presence of compact nebulae around a number of OB and WR stars is unexpected; by analogy to the cool super-/hypergiants we attribute this to confinement and sculpting of the stellar wind via interaction with the intra-cluster medium/wind. Given the morphology of core collapse SNe depend on the nature of the pre-explosion circumstellar environment, if this hypothesis is correct then the properties of the explosion depend not just on the progenitor, but also the environment in which it is located.

We report new 5.5 GHz radio observations of the massive star cluster Westerlund 1, taken by the Australia Telescope Compact Array, detecting nine of the ten yellow hypergiants (YHGs) and red supergiants (RSGs) within the cluster. Eight of nine sources are spatially resolved. The nebulae associated with the YHGs Wd1-4a, -12a and -265 demonstrate a cometary morphology - the first time this phenomenon has been observed for such stars. This structure is also echoed in the ejecta of the RSGs Wd1-20 and -26; in each case the cometary tails are directed away from the cluster core. The nebular emission around the RSG Wd1-237 is less collimated than these systems but once again appears more prominent in the hemisphere facing the cluster. Considered as a whole, the nebular morphologies provide compelling evidence for sculpting via a physical agent associated with Westerlund 1, such as a cluster wind.

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

In the colliding-wind region of massive binaries, non-thermal radio emission occurs. This non-thermal radio emission (due to synchrotron radiation) has so far been observed at centimetre wavelengths. At millimetre wavelengths, the stellar winds and the colliding-wind region emit more thermal free-free radiation, and it is expected that any non-thermal contribution will be difficult or impossible to detect. We aim to determine if the material in the colliding-wind region contributes substantially to the observed millimetre fluxes of a colliding-wind binary. We also try to distinguish the synchrotron emission from the free-free emission. We monitored the massive binary Cyg OB2 #8A at 3 mm with the NOrthern Extended Millimeter Array (NOEMA) interferometer of the Institut de Radioastronomie Millimetrique (IRAM). The data were collected in 14 separate observing runs (in 2014 and 2016), and provide good coverage of the orbital period. The observed millimetre fluxes range between 1.1 and 2.3 mJy, and show phase-locked variability, clearly indicating that a large part of the emission is due to the colliding-wind region. A simple synchrotron model gives fluxes with the correct order of magnitude, but with a maximum that is phase-shifted with respect to the observations. Qualitatively this phase shift can be explained by our neglect of orbital motion on the shape of the colliding-wind region. A model using only free-free emission results in only a slightly worse explanation of the observations. Additionally, on the map of our observations we also detect the O6.5 III star Cyg OB2 #8B, for which we determine a 3 mm flux of 0.21 +- 0.033 mJy. The question of whether synchrotron radiation or free-free emission dominates the millimetre fluxes of Cyg OB2 #8A remains open. More detailed modelling of this system, based on solving the hydrodynamical equations, is required to give a definite answer.

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.

Mass-loss in massive stars plays a critical role in their evolution, although the precise mechanism(s) responsible - radiatively driven winds, impulsive ejection and/or binary interaction -remain uncertain. In this paper we present ALMA line and continuum observations of the supergiant B[e] star Wd1-9, a massive post-Main Sequence object located within the starburst cluster Westerlund 1. We find it to be one of the brightest stellar point sources in the sky at millimetre wavelengths, with (serendipitously identified) emission in the H41alpha radio recombination line. We attribute these properties to a low velocity (~100 km/s) ionised wind, with an extreme mass-loss rate 6.4x10^-5(d/5kpc)^1.5 Msol/yr. External to this is an extended aspherical ejection nebula indicative of a prior phase of significant mass-loss. Taken together, the millimetre properties of Wd1-9 show a remarkable similarity to those of the highly luminous stellar source MWC349A.We conclude that these objects are interacting binaries evolving away from the main sequence and undergoing rapid case-A mass transfer. As such they - and by extension the wider class of supergiant B[e] stars - may provide a unique window into the physics of a process that shapes the life-cycle of ~70% of massive stars found in binary systems.

We present e-MERLIN 21 cm (L-band) observations of single luminous OB stars in the Cygnus OB2 association, from the COBRaS Legacy programme. The radio observations potentially offer the most straightforward, least model-dependent, determinations of mass-loss rates, and can be used to help resolve current discrepancies in mass-loss rates via clumped and structured hot star winds. We report here that the 21 cm flux densities of O3 to O6 supergiant and giant stars are less than ~ 70 microJy. These fluxes may be translated to `smooth' wind mass-loss upper limits of ~ 4.4 - 4.8 x 10^(-6) M_sol/yr for O3 supergiants and < 2.9 x 10^(-6) M_sol/yr for B0 to B1 supergiants. The first ever resolved 21 cm detections of the hypergiant (and LBV candidate) Cyg OB2 #12 are discussed; for multiple observations separated by 14 days, we detect a ~ 69% increase in its flux density. Our constraints on the upper limits for the mass-loss rates of evolved OB stars in Cyg OB2 support the model that the inner wind region close to the stellar surface (where Halpha forms) is more clumped than the very extended geometric region sampled by our radio observations.

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.

The Scripted E-merlin Rfi-mitigation PipelinE for iNTerferometry (SERPent) is an automated reduction and RFI-mitigation procedure utilising the SumThreshold methodology (Offringa et al. 2010b), originally developed for the LOFAR pipeline. SERPent is written in the Parseltongue language enabling interaction with the Astronomical Image Processing Software (AIPS) program. Moreover, SERPent is a simple "out of the box" Python script, which is easy to set up and is free of compilers. In addition to the flagging of RFI affected visibilities, the script also flags antenna zero-amplitude dropouts and Lovell telescope phase calibrator stationary scans inherent to the e-MERLIN system. Both the flagging and computational performances of SERPent are presented here, for e-MERLIN commissioning datasets for both L-band (1.3 - 1.8 GHz) and C-band (4 - 8 GHz) observations. RFI typically amounts to < 20 - 25% for the more problematic L-band observations and < 5% for the generally RFI quieter C-band. The level of RFI detection and flagging is more accurate and delicate than visual manual flagging, with the output immediately ready for AIPS calibration. SERPent is fully parallelised and has been tested on a range of computing systems. The current flagging rate is at 110 GB/day on a "high-end" computer (16 CPUs, 100 GB memory) which amounts to ~ 6.9 GB/CPU/day, with an expected increase in performance when e-MERLIN has completed its commissioning. The refining of automated reduction and calibration procedures is essential for the e-MERLIN legacy projects and future interferometers such as the SKA and the associated pathfinders (MeerKAT and ASKAP), where the vast data sizes (> TB) make traditional astronomer interactions unfeasible.

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 CLEAN algorithm, widely used in radio interferometry for the deconvolution of radio images, performs well only if the raw radio image (dirty image) is, to good approximation, a simple convolution between the instrumental point-spread function (dirty beam) and the true distribution of emission across the sky. An important case in which this approximation breaks down is during frequency synthesis if the observing bandwidth is wide enough for variations in the spectrum of the sky to become significant. The convolution assumption also breaks down, in any situation but snapshot observations, if sources in the field vary significantly in flux density over the duration of the observation. Such time-variation can even be instrumental in nature, for example due to jitter or rotation of the primary beam pattern on the sky during an observation. An algorithm already exists for dealing with the spectral variation encountered in wide-band frequency synthesis interferometry. This algorithm is an extension of CLEAN in which, at each iteration, a set of N `dirty beams' are fitted and subtracted in parallel, instead of just a single dirty beam as in standard CLEAN. In the wide-band algorithm the beams are obtained by expanding a nominal source spectrum in a Taylor series, each term of the series generating one of the beams. In the present paper this algorithm is extended to images which contain sources which vary over both frequency and time. Different expansion schemes (or bases) on the time and frequency axes are compared, and issues such as Gibbs ringing and non-orthogonality are discussed. It is shown that practical considerations make it often desirable to orthogonalize the set of beams before commencing the cleaning. This is easily accomplished via a Gram-Schmidt technique.

We present first results from the ongoing radio monitoring of SN 2008iz in M82. The VLBI images reveal a shell-like structure with circular symmetry, which expands in a self-similar way. There is strong evidence of a compact component with a steep spectrum at the center of the shell. The expansion curve obtained from our VLBI observations is marginally decelerated (m = 0.89) and can be modeled simultaneously with the available radio light curves. While the results of this simultaneous fitting are not conclusive (i.e. different combinations of values of the magnetic field, CSM density profile, and electron energy distribution, provide fits to the available data with similar quality), additional observations should allow a more robust and detailed modeling.

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.

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.

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.

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.