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.

VV 655, a dwarf irregular galaxy with HI tidal debris, is a companion to the lenticular luminous infrared galaxy (LIRG) NGC 4418. NGC 4418 stands out among nearby LIRGs due to its dense central concentration of molecular gas and the dusty, bi-polar structures along its minor axis suggestive of a wind driven by a central starburst and possible nuclear activity. We seek to understand the consequences of the ongoing minor interaction between VV 655 and NGC 4418 for the evolution of the LIRG, including the origin of the gas supply responsible for its unusual nuclear properties. We investigate the structural, kinematic, and chemical properties of VV 655 and NGC 4418 by analyzing archival imaging data and optical spectroscopic observations from the SDSS-III and new spectra from SALT-RSS. We characterize their gas-phase metal abundances and spatially resolved, ionized gas kinematics, and demonstrate that the gas-phase metallicity in NGC 4418 significantly exceeds that in VV 655. No kinematic disturbances in the ionized gas are observed along the minor axis of NGC 4418, but we see evidence for ionized gas outflows from VV 655 that may increase the cross-section for gas stripping in grazing collisions. A faint, asymmetric outer arm is detected in NGC 4418 of the type normally associated with galaxy-galaxy interactions. The simplest model suggests that the minor interaction between VV 655 and NGC 4418 produced the unusual nuclear properties of the LIRG via tidal torquing of the interstellar medium of NGC 4418 rather than through a significant gas transfer event. In addition to inducing a central concentration of gas in NGC 4418, this interaction also produced an enhanced star formation rate and an outer tidal arm in the LIRG. The VV 655-NGC 4418 system offers an example of the potential for minor collisions to alter the evolutionary pathways of giant galaxies.

High-resolution (0."03 - 0."09) (9 - 26 pc)) ALMA (100 - 350 GHz (3 - 0.8 mm)) and VLA 45 GHz measurements have been used to image continuum and spectral line emission from the inner region of the nearby infrared luminous galaxy IC 860. We detect compact (r=8-10 pc), luminous, mm continuum emission in the core of IC 860, with brightness temperatures $T_{\rm B}>$ 280 K. The 45 GHz continuum is equally compact, but fainter in flux. We suggest that the 3 to 0.8 mm continuum is opaque at mm-wavelengths, implying a very large H$_2$ column density $N$(H$_2$) of 1e26 cm-2, and that it emerges from hot dust with temperature $T_{\rm d} \sim 280$ K. Vibrationally excited lines of HCN $\nu_2$=1f J=4-3 and 3-2 (HCN-VIB) are seen in emission and resolved on scales of 40-50 pc. The HCN-VIB emission reveals a north-south nuclear velocity gradient with projected rotation velocities of v=100 kms$^{-1}$ at r=10 pc. The brightest HCN-VIB emission is oriented perpendicular to the velocity gradient. Ground state lines of HCN 3-2, 4-3, HC$^{15}$N 4-3, HCO$^+$ 3-2, 4-3 and CS 7-6, show complex line absorption and emission features. HCN and HCO$^+$ have red-shifted, reversed P-Cygni profiles consistent with gas inflows of $\sim$50 kms$^{-1}$. The absorption features can be traced from the north-east into the nucleus. In contrast, CS 7-6 shows blue-shifted line wings extending to -180 kms$^{-1}$. We suggest that a dense and slow outflow is hidden behind a foreground layer of inflowing gas. It appears that the centre of IC 860 is in a phase of rapid evolution where an inflow is building up the nuclear column density of gas. A slow, dense outflow may be signalling the onset of feedback. The inner, r=10 pc, IR luminosity can be powered by an accreting black hole and/or a compact starburst with a top-heavy initial mass function.

The nearby ultra-luminous infrared galaxy (ULIRG) Arp 220 is an excellent laboratory for studies of extreme astrophysical environments. For 20 years, Very Long Baseline Interferometry (VLBI) has been used to monitor a population of compact sources thought to be supernovae (SNe), supernova remnants (SNRs) and possibly active galactic nuclei (AGNs). Using new and archival VLBI data spanning 20 years, we obtain 23 high-resolution radio images of Arp 220 at wavelengths from 18 cm to 2 cm. From model-fitting to the images we obtain estimates of flux densities and sizes of all detected sources. We detect radio continuum emission from 97 compact sources and present flux densities and sizes for all analysed observation epochs. We find evidence for a LD-relation within Arp 220, with larger sources being less luminous. We find a compact source LF $n(L)\propto L^\beta$ with $\beta=-2.19\pm0.15$, similar to SNRs in normal galaxies. Based on simulations we argue that there are many relatively large and weak sources below our detection threshold. The observations can be explained by a mixed population of SNe and SNRs, where the former expand in a dense circumstellar medium (CSM) and the latter interact with the surrounding interstellar medium (ISM). Nine sources are likely luminous, type IIn SNe. This number of luminous SNe correspond to few percent of the total number of SNe in Arp 220 which is consistent with a total SN-rate of 4 yr$^{-1}$ as inferred from the total radio emission given a normal stellar initial mass function (IMF). Based on the fitted luminosity function, we argue that emission from all compact sources, also below our detection threshold, make up at most 20\% of the total radio emission at GHz frequencies.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

Previous Very Long Baseline Interferometry (VLBI) observations of the quasar B1152+199 at 5GHz has revealed two images of a strongly lensed jet with seemingly discordant morphologies. Whereas the jet appears straight in one of the images, the other exhibits slight curvature on milliarcsecond scales. This is unexpected from the lensing solution and has been interpreted as possible evidence for secondary, small-scale lensing (millilensing) by a compact object with a mass of $~10^5$-$10^7\ M_\odot$ located close to the curved image. The probability for such a superposition is extremely low unless the millilens population has very high surface number density. Here, we revisit the case for millilensing in B1152+199 by combining new global-VLBI data at 8.4GHz with two datasets from the European VLBI Network (EVN) at 5GHz (archival) and at 22GHz (new dataset), and the previously published 5GHz Very Long Baseline Array (VLBA) data. We find that the new data with a more circular synthesized beam, exhibits no apparent milliarcsecond-scale curvature in image B. Various observations of the object spanning $\sim$15 years apart enable us to improve the constraints on lens system (thanks also to the improved astrometry resulting from 22GHz observations) to the point that the only plausible explanation left for the apparent curvature is the artifact due to the shape of the synthesized beam.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

Minor mergers are important processes contributing significantly to how galaxies evolve across the age of the Universe. Their impact on supermassive black hole growth and star formation is profound. The detailed study of dense molecular gas in galaxies provides an important test of the validity of the relation between star formation rate and HCN luminosity on different galactic scales. We use observations of HCN, HCO+1-0 and CO3-2 to study the dense gas properties in the Medusa merger. We calculate the brightness temperature ratios and use them in conjunction with a non-LTE radiative line transfer model. The HCN and HCO+1-0, and CO3-2 emission do not occupy the same structures as the less dense gas associated with the lower-J CO emission. The only emission from dense gas is detected in a 200pc region within the "Eye of the Medusa". No HCN or HCO+ is detected for the extended starburst. The CO3-2/2-1 brightness temperature ratio inside "the Eye" is ~2.5 - the highest ratio found so far. The line ratios reveal an extreme, fragmented molecular cloud population inside "the Eye" with large temperatures (>300K) and high gas densities (>10^4 cm^-3). "The Eye" is found at an interface between a large-scale minor axis inflow and the Medusa central region. The extreme conditions inside "the Eye" may be the result of the radiative and mechanical feedback from a deeply embedded, young, massive super star cluster, formed due to the gas pile-up at the intersection. Alternatively, shocks from the inflowing gas may be strong enough to shock and fragment the gas. For both scenarios, however, it appears that the HCN and HCO+ dense gas tracers are not probing star formation, but instead a post-starburst and/or shocked ISM that is too hot and fragmented to form new stars. Thus, caution is advised in linking the detection of emission from dense gas tracers to evidence of ongoing or imminent star formation.

We report on the first sub-arcsecond (0.44 $\times$ 0.41 arcsec$\rm ^2$) angular resolution image at 150 MHz of the A-nucleus in the Luminous Infrared Galaxy Arp$\,$299, from International Low Frequency Array (LOFAR) Telescope observations. The most remarkable finding is that of an intriguing two-sided, filamentary structure emanating from A-nucleus, which we interpret as an outflow that extends up to at least 14 arcseconds from the A-nucleus in the N-S direction ($\approx$ 5 kpc deprojected size) and accounts for almost 40% of the extended emission of the entire galaxy system. We also discuss HST/NICMOS [FeII] 1.64 $\rm \mu m$ and H$\rm_2$ 2.12 $\rm \mu m$ images of Arp$\,$299-A, which show similar features to those unveiled by our 150 MHz LOFAR observations, thus giving string morphological support for the outflow scenario. Finally, we discuss unpublished NaI D spectra that confirm the outflow nature of this structure. From energetic arguments, we rule out the low-luminosity active galactic nucleus in Arp$\,$299-A as a driver for the outflow. On the contrary, the powerful, compact starburst in the central regions of Arp$\,$299-A provides plenty of mechanical energy to sustain an outflow, and we conclude that the intense supernova (SN) activity in the nuclear region of Arp299-A is driving the observed outflow. We estimate that the starburst wind can support a mass-outflow rate in the range (11-63) $\rm M_{\odot} yr^{-1}$ at speeds of up to (370 - 890) $\rm km \, s^{-1}$, and is relatively young, with an estimated kinematic age of (3 - 7) Myr. Those results open an avenue to the use of low-frequency (150 MHz), sub-arcsecond imaging with LOFAR to detect outflows in the central regions of local luminous infrared galaxies.

The majority of searches for radio emission from exoplanets have to date focused on short period planets, i.e., the so-called hot Jupiter type planets. However, these planets are likely to be tidally locked to their host stars and may not generate sufficiently strong magnetic fields to emit electron cyclotron maser emission at the low frequencies used in observations (typically >150 MHz). In comparison, the large mass-loss rates of evolved stars could enable exoplanets at larger orbital distances to emit detectable radio emission. Here, we first show that the large ionized mass-loss rates of certain evolved stars relative to the solar value could make them detectable with the Low Frequency Array (LOFAR) at 150 MHz ($\lambda$ = 2 m), provided they have surface magnetic field strengths >50 G. We then report radio observations of three long period (>1 au) planets that orbit the evolved stars $\beta$ Gem, $\iota$ Dra, and $\beta$ UMi using LOFAR at 150 MHz. We do not detect radio emission from any system but place tight 3$\sigma$ upper limits of 0.98, 0.87, and 0.57 mJy on the flux density at 150 MHz for $\beta$ Gem, $\iota$ Dra, and $\beta$ UMi, respectively. Despite our non-detections these stringent upper limits highlight the potential of LOFAR as a tool to search for exoplanetary radio emission at meter wavelengths.

Context. The importance of the magnetic field in high-mass-star formation is not yet fully clear and there are still many open questions concerning its role in the accretion processes and generation of jets and outflows. In the past few years, masers have been successfully used to probe the magnetic field morphology and strength at scales of a few au around massive protostars, by measuring linear polarisation angles and Zeeman splitting. The massive protostar IRAS 18089-1732 is a well studied high-mass-star forming region, showing a hot core chemistry and a disc-outflow system. Previous SMA observations of polarised dust revealed an ordered magnetic field oriented around the disc of IRAS 18089-1732. Aims. We want to determine the magnetic field in the dense region probed by 6.7 GHz methanol maser observations and compare it with observations in dust continuum polarisation, to investigate how the magnetic field in the compact maser region relates to the large-scale field around massive protostars. Methods. We reduced MERLIN observations at 6.7 GHz of IRAS 18089-1732 and we analysed the polarised emission by methanol masers. Results. Our MERLIN observations show that the magnetic field in the 6.7 GHz methanol maser region is consistent with the magnetic field constrained by the SMA dust polarisation observations. A tentative detection of circularly polarised line emission is also presented. Conclusions. We found that the magnetic field in the maser region has the same orientation as in the disk. Thus the large-scale field component, even at the au scale of the masers, dominates over any small-scale field fluctuations. We obtained, from the circular polarisation tentative detection, a field strength along the line of sight of 5.5 mG which appeared to be consistent with the previous estimates.

The galaxy NGC 4418 harbours a compact ($<20$ pc) core with a very high bolometric luminosity ($\sim10^{11}$L$_\odot$). As most of the galaxy's energy output comes from this small region, it is of interest to determine what fuels this intense activity. An interaction with VV 655 has been proposed, where gas aquired by NGC 4418 could trigger intense star formation and/or black hole accretion in the centre. We aim to constrain the interaction hypothesis by studying neutral hydrogen structures around the two galaxies. We present observations at 1.4 GHz with the Very Large Array of radio continuum as well as emission and absorption from atomic hydrogen. Gaussian distributions are fitted to observed HI emission and absorption spectra. An atomic HI bridge is seen in emission, connecting NGC 4418 to VV 655. While NGC 4418 is bright in continuum emission and seen in HI absorption, VV 655 is barely detected in the continuum but show bright HI emission (M$_\mathrm{HI}\sim10^9$ M$_\odot$). We estimate SFRs from 1.4 GHz of 3.2 M$_\odot$ yr$^{-1}$ and 0.13 M$_\odot$ yr$^{-1}$ for NGC 4418 and VV 655 respectively. Systemic HI velocities of 2202$\pm$20 km s$^{-1}$ (emission) and 2105.4$\pm$10 km s$^{-1}$ (absorption) are measured for VV 655 and NGC 4418 respectively. Redshifted HI absorption is seen towards NGC 4418, suggesting gas infall. Blueshifted HI-emission is seen north-west of NGC 4418, which we interpret as a continuation of the outflow previously discussed by Sakamoto et al. (2013). The morphology and velocity structure seen in HI is consistent with an interaction scenario, where gas was transferred from VV 655 to NGC 4418, and may fuel the activity in the centre.

Extragalactic cosmic ray populations are important diagnostic tools for tracking the distribution of energy in nuclei and for distinguishing between activity powered by star formation versus active galactic nuclei (AGNs). Here, we compare different diagnostics of the cosmic ray populations of the nuclei of Arp 220 based on radio synchrotron observations and the recent gamma-ray detection. We find the gamma-ray and radio emission to be incompatible; a joint solution requires at minimum a factor of 4 - 8 times more energy coming from supernovae and a factor of 40 - 70 more mass in molecular gas than is observed. We conclude that this excess of gamma-ray flux in comparison to all other diagnostics of star-forming activity indicates that there is an AGN present that is providing the extra cosmic rays, likely in the western nucleus.

The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s).

(abridged). We outline LBCS (the LOFAR Long-Baseline Calibrator Survey), whose aim is to identify sources suitable for calibrating the highest-resolution observations made with the International LOFAR Telescope, which include baselines >1000 km. Suitable sources must contain significant correlated flux density (50-100mJy) at frequencies around 110--190~MHz on scales of a few hundred mas. At least for the 200--300-km international baselines, we find around 1 suitable calibrator source per square degree over a large part of the northern sky, in agreement with previous work. This should allow a randomly selected target to be successfully phase calibrated on the international baselines in over 50% of cases. Products of the survey include calibrator source lists and fringe-rate and delay maps of wide areas -- typically a few degrees -- around each source. The density of sources with significant correlated flux declines noticeably with baseline length over the range 200--600~km, with good calibrators on the longest baselines appearing only at the rate of 0.5 per square degree. Coherence times decrease from 1--3 minutes on 200-km baselines to about 1 minute on 600-km baselines, suggesting that ionospheric phase variations contain components with scales of a few hundred kilometres. The longest median coherence time, at just over 3 minutes, is seen on the DE609 baseline, which at 227km is close to being the shortest. We see median coherence times of between 80 and 110 seconds on the four longest baselines (580--600~km), and about 2 minutes for the other baselines. The success of phase transfer from calibrator to target is shown to be influenced by distance, in a manner that suggests a coherence patch at 150-MHz of the order of 1 degree.

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

The correlation between radio spectral index and redshift has been exploited to discover high redshift radio galaxies, but its underlying cause is unclear. It is crucial to characterise the particle acceleration and loss mechanisms in high redshift radio galaxies to understand why their radio spectral indices are steeper than their local counterparts. Low frequency information on scales of $\sim$1 arcsec are necessary to determine the internal spectral index variation. In this paper we present the first spatially resolved studies at frequencies below 100 MHz of the $z = 2.4$ radio galaxy 4C 43.15 which was selected based on its ultra-steep spectral index ($\alpha < -1$; $S_{\nu} \sim \nu^{\alpha}$ ) between 365 MHz and 1.4 GHz. Using the International Low Frequency Array (LOFAR) Low Band Antenna we achieve sub-arcsecond imaging resolution at 55 MHz with VLBI techniques. Our study reveals low-frequency radio emission extended along the jet axis, which connects the two lobes. The integrated spectral index for frequencies $<$ 500 MHz is -0.83. The lobes have integrated spectral indices of -1.31$\pm$0.03 and -1.75$\pm$0.01 for frequencies $\geq$1.4 GHz, implying a break frequency between 500 MHz and 1.4 GHz. These spectral properties are similar to those of local radio galaxies. We conclude that the initially measured ultra-steep spectral index is due to a combination of the steepening spectrum at high frequencies with a break at intermediate frequencies.

We present the Multifrequency Snapshot Sky Survey (MSSS), the first northern-sky LOFAR imaging survey. In this introductory paper, we first describe in detail the motivation and design of the survey. Compared to previous radio surveys, MSSS is exceptional due to its intrinsic multifrequency nature providing information about the spectral properties of the detected sources over more than two octaves (from 30 to 160 MHz). The broadband frequency coverage, together with the fast survey speed generated by LOFAR's multibeaming capabilities, make MSSS the first survey of the sort anticipated to be carried out with the forthcoming Square Kilometre Array (SKA). Two of the sixteen frequency bands included in the survey were chosen to exactly overlap the frequency coverage of large-area Very Large Array (VLA) and Giant Metrewave Radio Telescope (GMRT) surveys at 74 MHz and 151 MHz respectively. The survey performance is illustrated within the "MSSS Verification Field" (MVF), a region of 100 square degrees centered at J2000 (RA,Dec)=(15h,69deg). The MSSS results from the MVF are compared with previous radio survey catalogs. We assess the flux and astrometric uncertainties in the catalog, as well as the completeness and reliability considering our source finding strategy. We determine the 90% completeness levels within the MVF to be 100 mJy at 135 MHz with 108" resolution, and 550 mJy at 50 MHz with 166" resolution. Images and catalogs for the full survey, expected to contain 150,000-200,000 sources, will be released to a public web server. We outline the plans for the ongoing production of the final survey products, and the ultimate public release of images and source catalogs.

We present ultraviolet (UV) observations of six nearby Type~Ia supernovae (SNe~Ia) obtained with the Hubble Space Telescope, three of which were also observed in the near-IR (NIR) with Wide-Field Camera~3. UV observations with the Swift satellite, as well as ground-based optical and near-infrared data provide complementary information. The combined data-set covers the wavelength range $0.2$--$2~\mu$m. By also including archival data of SN 2014J, we analyse a sample spanning observed colour excesses up to $E(B-V)=1.4~$mag. We study the wavelength dependent extinction of each individual SN and find a diversity of reddening laws when characterised by the total-to-selective extinction $R_V$. In particular, we note that for the two SNe with $E(B-V)\gtrsim1~$mag, for which the colour excess is dominated by dust extinction, we find $R_V=1.4\pm0.1$ and $R_V=2.8\pm0.1$. Adding UV photometry reduces the uncertainty of fitted $R_V$ by $\sim50\,$% allowing us to also measure $R_V$ of individual low-extinction objects which point to a similar diversity, currently not accounted for in the analyses when SNe~Ia are used for studying the expansion history of the universe.

The nuclear starburst in the nearby galaxy M82 provides an excellent laboratory for understanding the physics of star formation. This galaxy has been extensively observed in the past, revealing tens of radio-bright compact objects embedded in a diffuse free-free absorbing medium. Our understanding of the structure and physics of this medium in M82 can be greatly improved by high-resolution images at low frequencies where the effects of free-free absorption are most prominent. The aims of this study are, firstly, to demonstrate imaging using international baselines of the Low Frequency Array (LOFAR), and secondly, to constrain low-frequency spectra of compact and diffuse emission in the central starburst region of M82 via high-resolution radio imaging at low frequencies. The international LOFAR telescope was used to observe M82 at 110-126MHz and 146-162MHz. Images were obtained using standard techniques from very long baseline interferometry. images were obtained at each frequency range: one only using international baselines, and one only using the longest Dutch (remote) baselines. The 154MHz image obtained using international baselines is a new imaging record in terms of combined image resolution (0.3$"$) and sensitivity ($\sigma$=0.15mJy/beam) at low frequencies ($<327$MHz). We detected 16 objects at 154MHz, six of these also at 118MHz. Four weaker but resolved features are also found: a linear (50pc) filament and three other resolved objects, of which two show a clear shell structure. We do not detect any emission from either supernova 2008iz or from the radio transient source 43.78+59.3. The images obtained using remote baselines show diffuse emission, associated with the outflow in M82, with reduced brightness in the region of the edge-on star-forming disk.

Aims. An efficient means of locating calibrator sources for International LOFAR is developed and used to determine the average density of usable calibrator sources on the sky for subarcsecond observations at 140 MHz. Methods. We used the multi-beaming capability of LOFAR to conduct a fast and computationally inexpensive survey with the full International LOFAR array. Sources were pre-selected on the basis of 325 MHz arcminute-scale flux density using existing catalogues. By observing 30 different sources in each of the 12 sets of pointings per hour, we were able to inspect 630 sources in two hours to determine if they possess a sufficiently bright compact component to be usable as LOFAR delay calibrators. Results. Over 40% of the observed sources are detected on multiple baselines between international stations and 86 are classified as satisfactory calibrators. We show that a flat low-frequency spectrum (from 74 to 325 MHz) is the best predictor of compactness at 140 MHz. We extrapolate from our sample to show that the density of calibrators on the sky that are sufficiently bright to calibrate dispersive and non-dispersive delays for the International LOFAR using existing methods is 1.0 per square degree. Conclusions. The observed density of satisfactory delay calibrator sources means that observations with International LOFAR should be possible at virtually any point in the sky, provided that a fast and efficient search using the methodology described here is conducted prior to the observation to identify the best calibrator.

Carbon radio recombination lines (RRLs) at low frequencies (<=500 MHz) trace the cold, diffuse phase of the interstellar medium, which is otherwise difficult to observe. We present the detection of carbon RRLs in absorption in M82 with LOFAR in the frequency range of 48-64 MHz. This is the first extragalactic detection of RRLs from a species other than hydrogen, and below 1 GHz. Since the carbon RRLs are not detected individually, we cross-correlated the observed spectrum with a template spectrum of carbon RRLs to determine a radial velocity of 219 +- 9 km/s . Using this radial velocity, we stack 22 carbon-{\alpha} transitions from quantum levels n = 468-508 to achieve an 8.5 sigma detection. The absorption line profile exhibits a narrow feature with peak optical depth of 0.003 and FWHM of 31 km/s. Closer inspection suggests that the narrow feature is superimposed on a broad, shallow component. The total line profile appears to be correlated with the 21 cm H I line profile reconstructed from H I absorption in the direction of supernova remnants in the nucleus. The narrow width and centroid velocity of the feature suggests that it is associated with the nuclear starburst region. It is therefore likely that the carbon RRLs are associated with cold atomic gas in the direction of the nucleus of M82.

The galaxy NGC 4418 contains one of the most compact obscured nuclei within a luminous infrared galaxy (LIRG) in the nearby Universe. This nucleus contains a rich molecular gas environment and an unusually high ratio of infrared to radio luminosity (q-factor). The compact nucleus is powered by either a compact starburst or an active galactic nucleus (AGN). The aim of this study is to constrain the nature of the nuclear region (starburst or AGN) within NGC 4418 via very-high-resolution radio imaging. Archival data from radio observations using the EVN and MERLIN interferometers are imaged. Sizes and flux densities are obtained by fitting Gaussian intensity distributions to the image. The average spectral index of the compact radio emission is estimated from measurements at 1.4 GHz and 5.0 GHz. The nuclear structure of NGC 4418 visible with EVN and MERLIN consists of eight compact (<49 mas i.e. <8 pc) features spread within a region of 250 mas, i.e. 41 pc. We derive an inverted spectral index $\alpha\ge0.7$ ($S_\nu\propto\nu^{\alpha}$) for the compact radio emission. Brightness temperatures $>10^{4.8}$ K indicate that these compact features cannot be HII-regions. The complex morphology and inverted spectrum of the eight detected compact features is evidence against the hypothesis that an AGN alone is powering the nucleus of NGC 4418. The compact features could be super star clusters (SSCs) with intense star formation, and their associated free-free absorption could then naturally explain both their inverted radio spectrum and the low radio to IR ratio of the nucleus. The required star formation area density is extreme, however, and close to the limit of what can be observed in a well-mixed thermal/non-thermal plasma produced by star-formation, and is also close to the limit of what can be physically sustained.