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.

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

The population of intermediate-mass black holes (IMBHs) in nearby dwarf galaxies plays an important "ground truth" role in exploring black hole formation and growth in the early Universe. In the dwarf elliptical galaxy SDSS J090613.77+561015.2 (z=0.0465), an accreting IMBH has been revealed by optical and X-ray observations. Aiming to search for possible radio core and jet associated with the IMBH, we carried out very long baseline interferometry (VLBI) observations with the European VLBI Network (EVN) at 1.66 GHz. Our imaging results show that there are two 1-mJy components with a separation of about 52 mas (projected distance 47 pc) and the more compact component is located within the 1-sigma error circle of the optical centroid from available Gaia astrometry. Based on their positions, elongated structures and relatively high brightness temperatures, as well as the absence of star-forming activity in the host galaxy, we argue that the radio morphology originates from the jet activity powered by the central IMBH. The existence of the large-scale jet implies that violent jet activity might occur in the early epochs of black hole growth and thus help to regulate the co-evolution of black holes and galaxies.

The astrophysical agenda of the 21st century requires a very sharp view of celestial objects. High angular resolution studies are essential for fundamental studies of a broad variety of astrophysical phenomena ranging from relativistic physics of black holes, their gravitational and electromagnetic imprints, violent transient processes, including those producing detectable gravitational waves, birth and evolution of planetary systems. Over the past decades, radio astronomy made huge leap in achieving ground-breaking angular resolution measured in tens of microarcseconds (one tenth of nanoradian and better). Recently a global Event Horizon Telescope (EHT) collaboration obtained first direct images of the shadow of a super-massive black hole in the nucleus of the active galaxy M87. These observations were conducted at 230 GHz. The two first generation Space Very Long Baseline Interferometry (VLBI) missions, VSOP/HALCA led by the Japan Aerospace Exploration Agency (JAXA) and RadioAstron led by the Russia Roscosmos State Corporation and Russia Academy of Sciences, achieved the highest angular resolution at frequencies from 0.3 to 22 GHz in observations conducted in the period 1997 - 2019. The next step in advancing high angular resolution radio astronomy is in combining high frequency (millimeter and sub-millimeter wavelengths) and interferometric baselines exceeding the Earth diameter. The present THEZA White Paper describes a combination which would unify technology developments in giga-/tera-hertz instrumentation and space-borne radio astronomy. The current preprint version of the THEZA White Paper is slightly re-formatted and edited comparing to the official submitted version.

We present IRAM-30m Telescope $^{12}$CO and $^{13}$CO observations of a sample of 55 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in the local universe. This sample is a subset of the Great Observatory All-Sky LIRG Survey (GOALS), for which we use ancillary multi-wavelength data to better understand their interstellar medium and star formation properties. Fifty-three (96%) of the galaxies are detected in $^{12}$CO, and 29 (52%) are also detected in $^{13}$CO above a 3$\sigma$ level. The median full width at zero intensity (FWZI) velocity of the CO line emission is 661km s$^{-1}$, and $\sim$54% of the galaxies show a multi-peak CO profile. Herschel photometric data is used to construct the far-IR spectral energy distribution of each galaxy, which are fit with a modified blackbody model that allows us to derive dust temperatures and masses, and infrared luminosities. We make the assumption that the gas-to-dust mass ratio of (U)LIRGs is comparable to local spiral galaxies with a similar stellar mass (i.e., gas/dust of mergers is comparable to their progenitors) to derive a CO-to-H$_2$ conversion factor of $\langle\alpha\rangle=1.8^{+1.3}_{-0.8}M_\odot$(K km s$^{-1}$pc$^{2}$)$^{-1}$; such a value is comparable to that derived for (U)LIRGs based on dynamical mass arguments. We derive gas depletion times of $400-600$Myr for the (U)LIRGs, compared to the 1.3Gyr for local spiral galaxies. Finally, we re-examine the relationship between the $^{12}$CO/$^{13}$CO ratio and dust temperature, confirming a transition to elevated ratios in warmer systems.

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.

The shape of low-frequency radio continuum spectra of normal galaxies is not well understood, the key question being the role of physical processes such as thermal absorption in shaping them. In this work we take advantage of the LOFAR Multifrequency Snapshot Sky Survey (MSSS) to investigate such spectra for a large sample of nearby star-forming galaxies. Using the measured 150MHz flux densities from the LOFAR MSSS survey and literature flux densities at various frequencies we have obtained integrated radio spectra for 106 galaxies. The spectra are explained through the use of a three-dimensional model of galaxy radio emission, and radiation transfer dependent on the galaxy viewing angle and absorption processes. Spectra of our galaxies are generally flatter at lower compared to higher frequencies but as there is no tendency for the highly inclined galaxies to have more flattened low-frequency spectra, we argue that the observed flattening is not due to thermal absorption, contradicting the suggestion of Israel & Mahoney (1990). According to our modelled radio maps for M51-like galaxies, the free-free absorption effects can be seen only below 30MHz and in the global spectra just below 20MHz, while in the spectra of starburst galaxies, like M82, the flattening due to absorption is instead visible up to higher frequencies of about 150MHz. Locally, within galactic disks, the absorption effects are distinctly visible in M51-like galaxies as spectral flattening around 100-200MHz in the face-on objects, and as turnovers in the edge-on ones, while in M82-like galaxies there are strong turnovers at frequencies above 700MHz, regardless of viewing angle. Our modelling suggests that the weak spectral flattening observed in the nearby galaxies studied here results principally from synchrotron spectral curvature due to cosmic ray energy losses and propagation effects.

We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of intrinsic line strengths, adding to the difficulties in determining abundances from line observations. Our formalism is applicable to arbitrary distributions of cloud properties, provided the cloud volume filling factor is small; numerical simulations show it to hold up to filling factors of about 10%. We show that irrespective of the complexity of the cloud ensemble, the radiative effect of clumpiness can be parametrized at each frequency by a single multiplicative correction to the overall optical depth; this multiplier is derived from appropriate averaging over individual cloud properties. Our main finding is that cloud shapes have only a negligible effect on radiation propagation in clumpy media; the results of calculations employing point-like clouds are practically indistinguishable from those for finite-size clouds with arbitrary geometrical shapes.

The new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of the paper is to develop a computationally efficient and rigorously defined source-finding algorithm for linearly polarized sources. We used a calibrated data set from the LOw Frequency ARray (LOFAR) at 150 MHz centered on the nearby galaxy M51 to search for polarized background sources. With a new imaging software, we re-imaged the field at a resolution of 18''x15'' and cataloged a total of about 3000 continuum sources within 2.5 degrees of the center of M51. We made small Stokes Q and U images centered on each source brighter than 100 mJy in total intensity (201 sources) and used RM synthesis to create corresponding Faraday cubes that were analyzed individually. For each source, the noise distribution function was determined from a subset of measurements at high Faraday depths where no polarization is expected; the peaks in polarized intensity in the Faraday spectrum were identified and the p-value of each source was calculated. Finally, the false discovery rate method was applied to the list of p-values to produce a list of polarized sources and quantify the reliability of the detections. We also analyzed sources fainter than 100 mJy but that were reported as polarized in the literature at at least another radio frequency. Of the 201 sources that were searched for polarization, 6 polarized sources were detected confidently (with a false discovery rate of 5 percent). This corresponds to a number density of 1 polarized source per 3.3 square degrees, or 0.3 source/square degree. [Abridged]

Context: We describe the new SEPIA (Swedish-ESO PI Instrument for APEX) receiver, which was designed and built by the Group for Advanced Receiver Development (GARD), at Onsala Space Observatory (OSO) in collaboration with ESO. It was installed and commissioned at the APEX telescope during 2015 with an ALMA Band 5 receiver channel and updated with a new frequency channel (ALMA Band 9) in February 2016. Aims: This manuscript aims to provide, for observers who use the SEPIA receiver, a reference in terms of the hardware description, optics and performance as well as the commissioning results. Methods: Out of three available receiver cartridge positions in SEPIA, the two current frequency channels, corresponding to ALMA Band 5, the RF band 158--211 GHz, and Band 9, the RF band 600--722 GHz, provide state-of-the-art dual polarization receivers. The Band 5 frequency channel uses 2SB SIS mixers with an average SSB noise temperature around 45K with IF (intermediate frequency) band 4--8 GHz for each sideband providing total 4x4 GHz IF band. The Band 9 frequency channel uses DSB SIS mixers with a noise temperature of 75--125K with IF band 4--12 GHz for each polarization. Results: Both current SEPIA receiver channels are available to all APEX observers.

MHONGOOSE is a deep survey of the neutral hydrogen distribution in a representative sample of 30 nearby disk and dwarf galaxies with HI masses from 10^6 to ~10^{11} M_sun, and luminosities from M_R ~ -12 to M_R ~ -22. The sample is selected to uniformly cover the available range in log(M_HI). Our extremely deep observations, down to HI column density limits of well below 10^{18} cm^{-2} - or a few hundred times fainter than the typical HI disks in galaxies - will directly detect the effects of cold accretion from the intergalactic medium and the links with the cosmic web. These observations will be the first ever to probe the very low-column density neutral gas in galaxies at these high resolutions. Combination with data at other wavelengths, most of it already available, will enable accurate modelling of the properties and evolution of the mass components in these galaxies and link these with the effects of environment, dark matter distribution, and other fundamental properties such as halo mass and angular momentum. MHONGOOSE can already start addressing some of the SKA-1 science goals and will provide a comprehensive inventory of the processes driving the transformation and evolution of galaxies in the nearby universe at high resolution and over 5 orders of magnitude in column density. It will be a Nearby Galaxies Legacy Survey that will be unsurpassed until the advent of the SKA, and can serve as a highly visible, lasting statement of MeerKAT's capabilities.

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.

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 new IRAM 30m spectroscopic observations of the $\sim88$ GHz band, including emission from the CCH (n=1-0) multiplet, HCN (1-0), HCO+ (1-0), and HNC (1-0), for a sample of 58 local luminous and ultraluminous infrared galaxies from the Great Observatories All-sky LIRG Survey (GOALS). By combining our new IRAM data with literature data and Spitzer/IRS spectroscopy, we study the correspondence between these putative tracers of dense gas and the relative contribution of active galactic nuclei (AGN) and star formation to the mid-infrared luminosity of each system. We find the HCN (1-0) emission to be enhanced in AGN-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.84$), compared to composite and starburst-dominated systems ($\langle$L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}\rangle=1.14$, and 0.88, respectively). However, some composite and starburst systems have L'$_{HCN (1-0)}$/L'$_{HCO^+ (1-0)}$ ratios comparable to those of AGN, indicating that enhanced HCN emission is not uniquely associated with energetically dominant AGN. After removing AGN-dominated systems from the sample, we find a linear relationship (within the uncertainties) between $\log_{10}$(L'$_{HCN (1-0)}$) and $\log_{10}$(L$_{IR}$), consistent with most previous findings. L'$_{HCN (1-0)}$/L$_{IR}$, typically interpreted as the dense gas depletion time, appears to have no systematic trend with L$_{IR}$ for our sample of luminous and ultraluminous infrared galaxies, and has significant scatter. The galaxy-integrated HCN (1-0) and HCO+ (1-0) emission do not appear to have a simple interpretation, in terms of the AGN dominance or the star formation rate, and are likely determined by multiple processes, including density and radiative effects.

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 measure the sizes of redshift ~2 star-forming galaxies by stacking data from the Atacama Large Millimeter/submillimeter Array (ALMA). We use a uv-stacking algorithm in combination with model fitting in the uv-domain and show that this allows for robust measures of the sizes of marginally resolved sources. The analysis is primarily based on the 344 GHz ALMA continuum observations centred on 88 sub-millimeter galaxies in the LABOCA ECDFS Submillimeter Survey (ALESS). We study several samples of galaxies at z~2 with $M_* \sim{} 5\times{}10^{10} $ M$_\odot$ , selected using near-infrared photometry (distant red galaxies, extremely red objects, sBzK-galaxies, and galaxies selected on photometric redshift). We find that the typical sizes of these galaxies are ~0.6 arcsec which corresponds to ~5 kpc at z~2, this agrees well with the median sizes measured in the near-infrared z-band (~0.6 arcsec). We find errors on our size estimates of ~0.1-0.2 arcsec, which agree well with the expected errors for model fitting at the given signal-to-noise ratio. With the uv-coverage of our observations (18-160 m), the size and flux density measurements are sensitive to scales out to 2 . We compare this to a simulated ALMA Cycle 3 dataset with intermediate length baseline coverage, and we find that, using only these baselines, the measured stacked flux density would be an order of magnitude fainter. This highlights the importance of short baselines to recover the full flux density of high-redshift galaxies.

Aim: The aim of this work is to search Seyfert 2 galaxy NGC 4945, a well-known 22 GHz water megamaser galaxy, for water (mega)maser emission at 183 GHz. Method: We used APEX SEPIA Band 5 to perform the observations. Results: We detected 183 GHz water maser emission towards NGC 4945 with a peak flux density of ~3 Jy near the galactic systemic velocity. The emission spans a velocity range of several hundred km/s. We estimate an isotropic luminosity of > 1000 Lsun, classifying the emission as a megamaser. A comparison of the 183 GHz spectrum with that observed at 22 GHz suggests that 183 GHz emission also arises from the active galactic nucleus (AGN) central engine. If the 183 GHz emission originates from the circumnuclear disk, then we estimate that a redshifted feature at 1084 km/s in the spectrum should arise from a distance of 0.022 pc from the supermassive black hole (1.6 x 10(5) Schwarzschild radii), i.e. closer than the water maser emission previously detected at 22 GHz. This is only the second time 183 GHz maser emission has been detected towards an AGN central engine (the other galaxy being NGC 3079). It is also the strongest extragalactic millimetre/submillimetre water maser detected to date. Conclusions: Strong millimetre 183 GHz water maser emission has now been shown to occur in an external galaxy. For NGC 4945, we believe that the maser emission arises, or is dominated by, emission from the AGN central engine. Emission at higher velocity, i.e. for a Keplerian disk closer to the black hole, has been detected at 183 GHz compared with that for the 22 GHz megamaser. This indicates that millimetre/submillimetre water masers can indeed be useful probes for tracing out more of AGN central engine structures and dynamics than previously probed. Future observations using ALMA Band 5 should unequivocally determine the origin of the emission in this and other galaxies.

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.

We present LOFAR High-Band Array (HBA) observations of the Herschel-ATLAS North Galactic Pole survey area. The survey we have carried out, consisting of four pointings covering around 142 square degrees of sky in the frequency range 126--173 MHz, does not provide uniform noise coverage but otherwise is representative of the quality of data to be expected in the planned LOFAR wide-area surveys, and has been reduced using recently developed `facet calibration' methods at a resolution approaching the full resolution of the datasets ($\sim 10 \times 6$ arcsec) and an rms off-source noise that ranges from 100 $\mu$Jy beam$^{-1}$ in the centre of the best fields to around 2 mJy beam$^{-1}$ at the furthest extent of our imaging. We describe the imaging, cataloguing and source identification processes, and present some initial science results based on a 5-$\sigma$ source catalogue. These include (i) an initial look at the radio/far-infrared correlation at 150 MHz, showing that many Herschel sources are not yet detected by LOFAR; (ii) number counts at 150 MHz, including, for the first time, observational constraints on the numbers of star-forming galaxies; (iii) the 150-MHz luminosity functions for active and star-forming galaxies, which agree well with determinations at higher frequencies at low redshift, and show strong redshift evolution of the star-forming population; and (iv) some discussion of the implications of our observations for studies of radio galaxy life cycles.

We present the first wide area (19 deg$^2$), deep ($\approx120-150$ {\mu}Jy beam$^{-1}$), high resolution ($5.6 \times 7.4$ arcsec) LOFAR High Band Antenna image of the Bo\"otes field made at 130-169 MHz. This image is at least an order of magnitude deeper and 3-5 times higher in angular resolution than previously achieved for this field at low frequencies. The observations and data reduction, which includes full direction-dependent calibration, are described here. We present a radio source catalogue containing 6276 sources detected over an area of $19$\,deg$^2$, with a peak flux density threshold of $5\sigma$. As the first thorough test of the facet calibration strategy, introduced by van Weeren et al., we investigate the flux and positional accuracy of the catalogue. We present differential source counts that reach an order of magnitude deeper in flux density than previously achieved at these low frequencies, and show flattening at 150 MHz flux densities below 10 mJy associated with the rise of the low flux density star-forming galaxies and radio-quiet AGN.

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range.

Heterodyne receivers register the sky signal on either a circular polarization basis (where it is split into left-hand and right-hand circular polarization) or a linear polarization basis (where it is split into horizontal and vertical linear polarization). We study the problem of interferometric observations performed with telescopes that observe on different polarization bases, hence producing visibilities that we call "mixed basis" (i.e., linear in one telescope and circular in the other). We present novel algorithms for the proper calibration and treatment of such interferometric observations and test our algorithms with both simulations and real data. The use of our algorithms will be important for the optimum calibration of forthcoming observations with the Atacama Large mm/submm Array (ALMA) in very-long-baseline interferometry (VLBI) mode. Our algorithms will also allow us to optimally calibrate future VLBI observations at very high data rates (i.e., wide bandwidths), where linear-polarization feeds will be preferable at some stations, to overcome the polarimetric limitations due to the use of quarter-wave plates.

We present the results of a four-month campaign searching for low-frequency radio transients near the North Celestial Pole with the Low-Frequency Array (LOFAR), as part of the Multifrequency Snapshot Sky Survey (MSSS). The data were recorded between 2011 December and 2012 April and comprised 2149 11-minute snapshots, each covering 175 deg^2. We have found one convincing candidate astrophysical transient, with a duration of a few minutes and a flux density at 60 MHz of 15-25 Jy. The transient does not repeat and has no obvious optical or high-energy counterpart, as a result of which its nature is unclear. The detection of this event implies a transient rate at 60 MHz of 3.9 (+14.7, -3.7) x 10^-4 day^-1 deg^-2, and a transient surface density of 1.5 x 10^-5 deg^-2, at a 7.9-Jy limiting flux density and ~10-minute time-scale. The campaign data were also searched for transients at a range of other time-scales, from 0.5 to 297 min, which allowed us to place a range of limits on transient rates at 60 MHz as a function of observation duration.

Steep spectrum radio sources associated with active galactic nuclei (AGN) may contain remnants of past AGN activity episodes. Novel instruments like the LOw Frequency ARray (LOFAR) are enabling studies of these fascinating structures to be made at tens to hundreds of MHz with sufficient resolution to analyse their complex morphology. Our goal is to characterize the integrated and resolved spectral properties of VLSS J1431+1331 and estimate source ages based on synchrotron radio emission models, thus putting constraints on the AGN duty cycle. Using a broad spectral coverage, we have derived spectral and curvature maps, and used synchrotron ageing models to determine the time elapsed from the last time the source plasma was energized. We used LOFAR, Giant Metrewave Radio Telescope (GMRT) and Jansky Very Large Array (VLA) data. Based on our ageing analysis, we infer that the AGN that created this source currently has very low levels of activity or that it is switched off. The derived ages for the larger source component range from around 60 to 130 Myr, hinting that the AGN activity decreased or stopped around 60 Myr ago. Our analysis suggests that VLSS J1431.8+1331 is an intriguing, two-component source. The larger component seems to host a faint radio core, suggesting that the source may be an AGN radio relic. The spectral index we observe from the smaller component is distinctly flatter at lower frequencies than the spectral index of the larger component, suggesting the possibility that the smaller component may be a shocked plasma bubble. From the integrated source spectrum, we deduce that its shape and slope can be used as tracers of the activity history of this type of steep spectrum radio source.

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 report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

The Atacama Large mm-submm Array (ALMA) is, by far, the most sensitive mm/submm telescope in the World. The ALMA Phasing Project (APP) will allow us to phase-up all the ALMA antennas and use them as one single VLBI station. This will be a key component of the Event Horizon Telescope (EHT), a Global VLBI array at millimeter wavelengths. A problem in the APP is the calibration and conversion of the polarization channels. Most VLBI stations record their signals in a circular basis, but the ALMA receivers record in a linear basis. The strategy that will be followed in the phased-ALMA VLBI observations will be to correlate in "mixed" basis (i.e., linear versus circular) and convert the visibilities to a pure circular basis after the correlation. We have developed an algorithm to perform such a polarization conversion of the VLBI visibilities. In these proceedings, we present the basics of the PolConvert algorithm and discuss on the polarization conversion in the general case were single dishes (besides phased arrays) record with linear receivers in VLBI observations. We show some results of PolConvert applied to realistic simulations, as well as a test with real VLBI observations at 86\,GHz between the Onsala radiotelescope (recording in linear basis) and the Effelsberg radiotelescope (recording in circular basis).

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Core-collapse SNe (CCSNe): Systematic searches of radio emission from CCSNe are still lacking, and only targeted searches of radio emission from just some of the optically discovered CCSNe in the local universe have been carried out. Optical searches miss a significant fraction of CCSNe due to dust obscuration; CCSN radio searches are thus more promising for yielding the complete, unobscured star-formation rates in the local universe. The SKA yields the possibility to piggyback for free in this area of research by carrying out commensal, wide-field, blind transient survey observations. SKA1-SUR should be able to discover several hundreds of CCSNe in just one year, compared to about a dozen CCSNe that the VLASS would be able to detect in one year, at most. SKA, with an expected sensitivity ten times that of SKA1, is expected to detect CCSNe in the local Universe by the thousands. Therefore, commensal SKA observations could easily result in an essentially complete census of all CCSNe in the local universe, thus yielding an accurate determination of the volumetric CCSN rate. Type Ia SNe: We advocate for the use of the SKA to search for the putative prompt (~first few days after the explosion) radio emission of any nearby type Ia SN, via target-of-opportunity observations. The huge improvement in sensitivity of the SKA with respect to its predecessors will allow to unambiguously discern which progenitor scenario (single-degenerate vs. double-degenerate) applies to them.

Stacking as a tool for studying objects that are not individually detected is becoming popular even for radio interferometric data, and will be widely used in the SKA era. Stacking is typically done using imaged data rather than directly using the visibilities (the uv-data). We have investigated and developed a novel algorithm to do stacking using the uv-data. We have performed exten- sive simulations comparing to image-stacking, and summarize the results of these simulations. Furthermore, we disuss the implications in light of the vast data volume produced by the SKA. Having access to the uv-stacked data provides a great advantage, as it allows the possibility to properly analyse the result with respect to calibration artifacts as well as source properties such as size. For SKA the main challenge lies in archiving the uv-data. For purposes of robust stacking analysis, it would be strongly desirable to either keep the calibrated uv-data at least in an aver- age form, or implement a stacking queue where stacking positions could be provided prior to the observations and the uv-stacking is done almost in real time.

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.

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

We present a new algorithm for stacking radio interferometric data in the uv-domain. The performance of uv-stacking is compared to the stacking of fully imaged data using simulated Atacama Large Millimeter/sub-millimeter Array (ALMA) and the Karl G. Jansky Very Large Array (VLA) deep extragalactic surveys. We find that image- and uv-stacking produce similar results, however, uv-stacking is typically the more robust method. An advantage of the uv-stacking algorithm is the availability of uv-data post stacking, which makes it possible to identify and remove problematic baselines. For deep VLA surveys uv-stacking yields a signal-to-noise ratio that is up to 20 per cent higher than image-stacking. Furthermore, we have investigated stacking of resolved sources with a simulated VLA data set where 1.5" (10-12 kpc at z ~ 1-4) sources are stacked. We find that uv-stacking, where a model is fitted directly to the visibilities, significantly improves the accuracy and robustness of the size estimates. While scientific motivation for this work is studying faint, high-z galaxies, the algorithm analysed here would also be applicable in other fields of astronomy. Stacking of radio interferometric data is also expected to play a big role for future surveys with telescopes such as LOFAR and Square Kilometre Array (SKA).

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 grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates.

The Sun is an active source of radio emission which is often associated with energetic phenomena such as solar flares and coronal mass ejections (CMEs). At low radio frequencies (<100 MHz), the Sun has not been imaged extensively because of the instrumental limitations of previous radio telescopes. Here, the combined high spatial, spectral and temporal resolution of the Low Frequency Array (LOFAR) was used to study solar Type III radio bursts at 30-90 MHz and their association with CMEs. The Sun was imaged with 126 simultaneous tied-array beams within 5 solar radii of the solar centre. This method offers benefits over standard interferometric imaging since each beam produces high temporal (83 ms) and spectral resolution (12.5 kHz) dynamic spectra at an array of spatial locations centred on the Sun. LOFAR's standard interferometric output is currently limited to one image per second. Over a period of 30 minutes, multiple Type III radio bursts were observed, a number of which were found to be located at high altitudes (4 solar radii from the solar center at 30 MHz) and to have non-radial trajectories. These bursts occurred at altitudes in excess of values predicted by 1D radial electron density models. The non-radial high altitude Type III bursts were found to be associated with the expanding flank of a CME. The CME may have compressed neighbouring streamer plasma producing larger electron densities at high altitudes, while the non-radial burst trajectories can be explained by the deflection of radial magnetic fields as the CME expanded in the low corona.

This study aims to characterise the polarized foreground emission in the ELAIS-N1 field and to address its possible implications for the extraction of the cosmological 21-cm signal from the Low-Frequency Array - Epoch of Reionization (LOFAR-EoR) data. We use the high band antennas of LOFAR to image this region and RM-synthesis to unravel structures of polarized emission at high Galactic latitudes. The brightness temperature of the detected Galactic emission is on average 4 K in polarized intensity and covers the range from -10 to +13rad m^-2 in Faraday depth. The total polarized intensity and polarization angle show a wide range of morphological features. We have also used the Westerbork Synthesis Radio Telescope (WSRT) at 350 MHz to image the same region. The LOFAR and WSRT images show a similar complex morphology, at comparable brightness levels, but their spatial correlation is very low. The fractional polarization at 150 MHz, expressed as a percentage of the total intensity, amounts to 1.5%. There is no indication of diffuse emission in total intensity in the interferometric data, in line with results at higher frequencies. The wide frequency range, good angular resolution and good sensitivity make LOFAR an exquisite instrument for studying Galactic polarized emission at a resolution of 1-2 rad m^-2 in Faraday depth. The different polarised patterns observed at 150 MHz and 350 MHz are consistent with different source distributions along the line of sight wring in a variety of Faraday thin regions of emission. The presence of polarised foregrounds is a serious complication for Epoch of Reionization experiments. To avoid the leakage of polarized emission into total intensity, which can depend on frequency, we need to calibrate the instrumental polarization across the field of view to a small fraction of 1%.

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.

We report radio observations of two stripped-envelope supernovae (SNe), 2010O and 2010P, which exploded within a few days of each other in the luminous infrared galaxy Arp 299. Whilst SN 2010O remains undetected at radio frequencies, SN 2010P was detected (with an astrometric accuracy better than 1 milli arcsec in position) in its optically thin phase in epochs ranging from ~1 to ~3yr after its explosion date, indicating a very slow radio evolution and a strong interaction of the SN ejecta with the circumstellar medium. Our late-time radio observations toward SN 2010P probe the dense circumstellar envelope of this SN, and imply a mass-loss rate (Msun/yr) to wind velocity (in units of 10 km/s) ratio of (3.0-5.1)E-05, with a 5 GHz peak luminosity of ~1.2E+27 erg/s/Hz on day ~464 after explosion. This is consistent with a Type IIb classification for SN 2010P, making it the most distant and most slowly evolving Type IIb radio SN detected to date.

We present the first detection of carbon radio recombination line absorption along the line of sight to Cygnus A. The observations were carried out with the LOw Frequency ARray in the 33 to 57 MHz range. These low frequency radio observations provide us with a new line of sight to study the diffuse, neutral gas in our Galaxy. To our knowledge this is the first time that foreground Milky Way recombination line absorption has been observed against a bright extragalactic background source. By stacking 48 carbon $\alpha$ lines in the observed frequency range we detect carbon absorption with a signal-to-noise ratio of about 5. The average carbon absorption has a peak optical depth of 2$\times$10$^{-4}$, a line width of 10 km s$^{-1}$ and a velocity of +4 km s$^{-1}$ with respect to the local standard of rest. The associated gas is found to have an electron temperature $T_{e}\sim$ 110 K and density $n_{e}\sim$ 0.06 cm$^{-3}$. These properties imply that the observed carbon $\alpha$ absorption likely arises in the cold neutral medium of the Orion arm of the Milky Way. Hydrogen and helium lines were not detected to a 3$\sigma$ peak optical depth limit of 1.5$\times$10$^{-4}$ for a 4 km s$^{-1}$ channel width. Radio recombination lines associated with Cygnus A itself were also searched for, but are not detected. We set a 3$\sigma$ upper limit of 1.5$\times$10$^{-4}$ for the peak optical depth of these lines for a 4 km s$^{-1}$ channel width.

An international consortium is presently constructing a beamformer for the Atacama Large Millimeter/submillimeter Array (ALMA) in Chile that will be available as a facility instrument. The beamformer will aggregate the entire collecting area of the array into a single, very large aperture. The extraordinary sensitivity of phased ALMA, combined with the extremely fine angular resolution available on baselines to the Northern Hemisphere, will enable transformational new very long baseline interferometry (VLBI) observations in Bands 6 and 7 (1.3 and 0.8 mm) and provide substantial improvements to existing VLBI arrays in Bands 1 and 3 (7 and 3 mm). The ALMA beamformer will have impact on a variety of scientific topics, including accretion and outflow processes around black holes in active galactic nuclei (AGN), tests of general relativity near black holes, jet launch and collimation from AGN and microquasars, pulsar and magnetar emission processes, the chemical history of the universe and the evolution of fundamental constants across cosmic time, maser science, and astrometry.

The characteristic outer scale of turbulence and the ratio of the random to ordered components of the magnetic field are key parameters to characterise magnetic turbulence in the interstellar gas, which affects the propagation of cosmic rays within the Galaxy. We provide new constraints to those two parameters. We use the LOw Frequency ARray (LOFAR) to image the diffuse continuum emission in the Fan region at (l,b) (137.0,+7.0) at 80"x70" resolution in the range [146,174] MHz. We detect multi-scale fluctuations in the Galactic synchrotron emission and compute their power spectrum. Applying theoretical estimates and derivations from the literature for the first time, we derive the outer scale of turbulence and the ratio of random to ordered magnetic field from the characteristics of these fluctuations . We obtain the deepest image of the Fan region to date and find diffuse continuum emission within the primary beam. The power spectrum of the foreground synchrotron fluctuations displays a power law behaviour for scales between 100 and 8 arcmin with a slope of (-1.84+/-0.19). We find an upper limit of about 20 pc for the outer scale of the magnetic interstellar turbulence toward the Fan region. We also find a variation of the ratio of random to ordered field as a function of Galactic coordinates, supporting different turbulent regimes. We use power spectra fluctuations from LOFAR as well as earlier GMRT and WSRT observations to constrain the outer scale of turbulence of the Galactic synchrotron foreground, finding a range of plausible values of 10-20 pc. Then, we use this information to deduce lower limits of the ratio of ordered to random magnetic field strength. These are found to be 0.3, 0.3, and 0.5 for the LOFAR, WSRT and GMRT fields considered respectively. Both these constraints are in agreement with previous estimates.

LOFAR, the LOw-Frequency ARray, is a new-generation radio interferometer constructed in the north of the Netherlands and across europe. Utilizing a novel phased-array design, LOFAR covers the largely unexplored low-frequency range from 10-240 MHz and provides a number of unique observing capabilities. Spreading out from a core located near the village of Exloo in the northeast of the Netherlands, a total of 40 LOFAR stations are nearing completion. A further five stations have been deployed throughout Germany, and one station has been built in each of France, Sweden, and the UK. Digital beam-forming techniques make the LOFAR system agile and allow for rapid repointing of the telescope as well as the potential for multiple simultaneous observations. With its dense core array and long interferometric baselines, LOFAR achieves unparalleled sensitivity and angular resolution in the low-frequency radio regime. The LOFAR facilities are jointly operated by the International LOFAR Telescope (ILT) foundation, as an observatory open to the global astronomical community. LOFAR is one of the first radio observatories to feature automated processing pipelines to deliver fully calibrated science products to its user community. LOFAR's new capabilities, techniques and modus operandi make it an important pathfinder for the Square Kilometre Array (SKA). We give an overview of the LOFAR instrument, its major hardware and software components, and the core science objectives that have driven its design. In addition, we present a selection of new results from the commissioning phase of this new radio observatory.

Cassiopeia A was observed using the Low-Band Antennas of the LOw Frequency ARray (LOFAR) with high spectral resolution. This allowed a search for radio recombination lines (RRLs) along the line-of-sight to this source. Five carbon-alpha RRLs were detected in absorption between 40 and 50 MHz with a signal-to-noise ratio of > 5 from two independent LOFAR datasets. The derived line velocities (v_LSR ~ -50 km/s) and integrated optical depths (~ 13 s^-1) of the RRLs in our spectra, extracted over the whole supernova remnant, are consistent within each LOFAR dataset and with those previously reported. For the first time, we are able to extract spectra against the brightest hotspot of the remnant at frequencies below 330 MHz. These spectra show significantly higher (15-80 %) integrated optical depths, indicating that there is small-scale angular structure on the order of ~1 pc in the absorbing gas distribution over the face of the remnant. We also place an upper limit of 3 x 10^-4 on the peak optical depths of hydrogen and helium RRLs. These results demonstrate that LOFAR has the desired spectral stability and sensitivity to study faint recombination lines in the decameter band.

The aim of the LOFAR Epoch of Reionization (EoR) project is to detect the spectral fluctuations of the redshifted HI 21cm signal. This signal is weaker by several orders of magnitude than the astrophysical foreground signals and hence, in order to achieve this, very long integrations, accurate calibration for stations and ionosphere and reliable foreground removal are essential. One of the prospective observing windows for the LOFAR EoR project will be centered at the North Celestial Pole (NCP). We present results from observations of the NCP window using the LOFAR highband antenna (HBA) array in the frequency range 115 MHz to 163 MHz. The data were obtained in April 2011 during the commissioning phase of LOFAR. We used baselines up to about 30 km. With about 3 nights, of 6 hours each, effective integration we have achieved a noise level of about 100 microJy/PSF in the NCP window. Close to the NCP, the noise level increases to about 180 microJy/PSF, mainly due to additional contamination from unsubtracted nearby sources. We estimate that in our best night, we have reached a noise level only a factor of 1.4 above the thermal limit set by the noise from our Galaxy and the receivers. Our continuum images are several times deeper than have been achieved previously using the WSRT and GMRT arrays. We derive an analytical explanation for the excess noise that we believe to be mainly due to sources at large angular separation from the NCP.

M87 is a giant elliptical galaxy located in the centre of the Virgo cluster, which harbours a supermassive black hole of mass 6.4x10^9 M_sun, whose activity is responsible for the extended (80 kpc) radio lobes that surround the galaxy. The energy generated by matter falling onto the central black hole is ejected and transferred to the intra-cluster medium via a relativistic jet and morphologically complex systems of buoyant bubbles, which rise towards the edges of the extended halo. Here we present the first observations made with the new Low-Frequency Array (LOFAR) of M87 at frequencies down to 20 MHz. Images of M87 were produced at low radio frequencies never explored before at these high spatial resolution and dynamic range. To disentangle different synchrotron models and place constraints on source magnetic field, age and energetics, we also performed a detailed spectral analysis of M87 extended radio-halo using these observations together with archival data. We do not find any sign of new extended emissions; on the contrary the source appears well confined by the high pressure of the intra-cluster medium. A continuous injection of relativistic electrons is the model that best fits our data, and provides a scenario in which the lobes are still supplied by fresh relativistic particles from the active galactic nuclei. We suggest that the discrepancy between the low-frequency radio-spectral slope in the core and in the halo implies a strong adiabatic expansion of the plasma as soon as it leaves the core area. The extended halo has an equipartition magnetic field strength of ~10 uG, which increases to ~13 uG in the zones where the particle flows are more active. The continuous injection model for synchrotron ageing provides an age for the halo of ~40 Myr, which in turn provides a jet kinetic power of 6-10x10^44 erg/s.

Deep radio observations of galaxy clusters have revealed the existence of diffuse radio sources related to the presence of relativistic electrons and weak magnetic fields in the intracluster volume. The role played by this non-thermal intracluster component on the thermodynamical evolution of galaxy clusters is debated, with important implications for cosmological and astrophysical studies of the largest gravitationally bound structures of the Universe. The low surface brightness and steep spectra of diffuse cluster radio sources make them more easily detectable at low-frequencies. LOFAR is the first instrument able to detect diffuse radio emission in hundreds of massive galaxy clusters up to their formation epoch. We present the first observations of clusters imaged by LOFAR and the huge perspectives opened by this instrument for non-thermal cluster studies.

Aims: This paper discusses the spectral occupancy for performing radio astronomy with the Low-Frequency Array (LOFAR), with a focus on imaging observations. Methods: We have analysed the radio-frequency interference (RFI) situation in two 24-h surveys with Dutch LOFAR stations, covering 30-78 MHz with low-band antennas and 115-163 MHz with high-band antennas. This is a subset of the full frequency range of LOFAR. The surveys have been observed with a 0.76 kHz / 1 s resolution. Results: We measured the RFI occupancy in the low and high frequency sets to be 1.8% and 3.2% respectively. These values are found to be representative values for the LOFAR radio environment. Between day and night, there is no significant difference in the radio environment. We find that lowering the current observational time and frequency resolutions of LOFAR results in a slight loss of flagging accuracy. At LOFAR's nominal resolution of 0.76 kHz and 1 s, the false-positives rate is about 0.5%. This rate increases approximately linearly when decreasing the data frequency resolution. Conclusions: Currently, by using an automated RFI detection strategy, the LOFAR radio environment poses no perceivable problems for sensitive observing. It remains to be seen if this is still true for very deep observations that integrate over tens of nights, but the situation looks promising. Reasons for the low impact of RFI are the high spectral and time resolution of LOFAR; accurate detection methods; strong filters and high receiver linearity; and the proximity of the antennas to the ground. We discuss some strategies that can be used once low-level RFI starts to become apparent. It is important that the frequency range of LOFAR remains free of broadband interference, such as DAB stations and windmills.

Abell 2256 is one of the best known examples of a galaxy cluster hosting large-scale diffuse radio emission that is unrelated to individual galaxies. It contains both a giant radio halo and a relic, as well as a number of head-tail sources and smaller diffuse steep-spectrum radio sources. The origin of radio halos and relics is still being debated, but over the last years it has become clear that the presence of these radio sources is closely related to galaxy cluster merger events. Here we present the results from the first LOFAR Low band antenna (LBA) observations of Abell 2256 between 18 and 67 MHz. To our knowledge, the image presented in this paper at 63 MHz is the deepest ever obtained at frequencies below 100 MHz in general. Both the radio halo and the giant relic are detected in the image at 63 MHz, and the diffuse radio emission remains visible at frequencies as low as 20 MHz. The observations confirm the presence of a previously claimed ultra-steep spectrum source to the west of the cluster center with a spectral index of -2.3 \pm 0.4 between 63 and 153 MHz. The steep spectrum suggests that this source is an old part of a head-tail radio source in the cluster. For the radio relic we find an integrated spectral index of -0.81 \pm 0.03, after removing the flux contribution from the other sources. This is relatively flat which could indicate that the efficiency of particle acceleration at the shock substantially changed in the last \sim 0.1 Gyr due to an increase of the shock Mach number. In an alternative scenario, particles are re-accelerated by some mechanism in the downstream region of the shock, resulting in the relatively flat integrated radio spectrum. In the radio halo region we find indications of low-frequency spectral steepening which may suggest that relativistic particles are accelerated in a rather inhomogeneous turbulent region.

Context. The two nuclei of the starburst galaxy Arp220 contain multiple compact radio sources previously identified as radio supernovae or supernova remnants. Aims. In order to search for an embedded radio AGN, or other possible exotic objects, we have carried out a program of VLBI monitoring at 6 cm over three epochs each separated by four months. Methods. Combining the new data with existing data at 6 cm and 18 cm (spanning 4 and 12 years respectively) we are able to characterise source flux density variability on a range of time-scales. Additionally we analyse the variability of sources in shape and position. Results. We detect rapid (< 4 months) variability in three sources (W7, W26 and W29). These sources show possible superluminal motion (> 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources might be associated with an AGN or a highly beamed microquasar (i.e. microblazar). Other hypotheses include that the apparent variability is intrinsic and is produced by neutron star powered central components within a supernova remnant, by a sequence of several supernovae within super star clusters, or is extrinsic and is produced by Galactic interstellar scintillation of very compact non-varying objects. Conclusions. A microquasar/microblazar origin seems to be the best explanation for the nature of the variable sources in Arp220.

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 present 2 cm and 3.6 cm wavelength very long baseline interferometry images of the compact radio continuum sources in the nearby ultra-luminous infrared galaxy Arp220. Based on their radio spectra and variability properties, we confirm these sources to be a mixture of supernovae (SNe) and supernova remnants (SNRs). Of the 17 detected sources we resolve 7 at both wavelengths. The SNe generally only have upper size limits. In contrast all the SNRs are resolved with diameters {\geq} 0.27 pc. This size limit is consistent with them having just entered their Sedov phase while embedded in an interstellar medium (ISM) of density 10^4 cm^{-3} . These objects lie on the diameter-luminosity correlation for SNRs (and so also on the diameter-surface brightness relation) and extend these correlations to very small sources. The data are consistent with the relation L {\propto} D^{-9/4}. Revised equipartition arguments adjusted to a magnetic field to relativistic particle energy density ratio of 1% combined with a reasonable synchrotron-emitting volume filling factor of 10% give estimated magnetic field strengths in the SNR shells of ~ 15-50 mG. The SNR shell magnetic fields are unlikely to come from compression of ambient ISM fields and must instead be internally generated. We set an upper limit of 7 mG for the ISM magnetic field. The estimated energy in relativistic particles, 2%-20% of the explosion kinetic energy, is consistent with estimates from models that fit the IR-radio correlation in compact starburst galaxies.

At very low frequencies, the new pan-European radio telescope LOFAR is opening the last unexplored window of the electromagnetic spectrum for astrophysical studies. The revolutionary APERTIF phased arrays that are about to be installed on the Westerbork radio telescope (WSRT) will dramatically increase the survey speed for the WSRT. Combined surveys with these two facilities will deeply chart the northern sky over almost two decades in radio frequency from \sim 15 up to 1400 MHz. Here we briefly describe some of the capabilities of these new facilities and what radio surveys are planned to study fundamental issues related the formation and evolution of galaxies and clusters of galaxies. In the second part we briefly review some recent observational results directly showing that diffuse radio emission in clusters traces shocks due to cluster mergers. As these diffuse radio sources are relatively bright at low frequencies, LOFAR should be able to detect thousands of such sources up to the epoch of cluster formation. This will allow addressing many question about the origin and evolution of shocks and magnetic fields in clusters. At the end we briefly review some of the first and very preliminary LOFAR results on clusters.

The Low Frequency Array (LOFAR) is under construction in the Netherlands and in several surrounding European countries. In this contribution, we describe the layout and design of the telescope, with a particular emphasis on the imaging characteristics of the array when used in its "standard imaging" mode. After briefly reviewing the calibration and imaging software used for LOFAR image processing, we show some recent results from the ongoing imaging commissioning efforts. We conclude by summarizing future prospects for the use of LOFAR in observing the little-explored low frequency Universe.

The Low Frequency Array (LOFAR) will operate between 10 and 250 MHz, and will observe the low frequency Universe to an unprecedented sensitivity and angular resolution. The construction and commissioning of LOFAR is well underway, with over 27 of the Dutch stations and five International stations routinely performing both single-station and interferometric observations over the frequency range that LOFAR is anticipated to operate at. Here, we summarize the capabilities of LOFAR and report on some of the early commissioning imaging of Cygnus A.

Low frequency radio waves, while challenging to observe, are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric "radio window": 10-240MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth. LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals. We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR.

We present the most recent VLBI images of SN 1993J, taken at 1.7 GHz on 2010 March 5-6, along with a discussion of its evolution with time. The new image is the latest in a sequence covering almost the entire lifetime of the supernova. For these latest observations we used an "in beam calibrator" technique, and obtained a background rms brightness of 3.7 micro-Jy/beam. The supernova shell remains quite circular in outline. Modulations in brightness are seen around the rim which evolve relatively slowly, having remained generally similar over the last several years of observation. We determine the outer radius of the supernova using visibility-plane model-fitting. The supernova has slowed down to around 30% of its original expansion velocity, and continues to expand with radius approximately proportional to t^0.8, however, deviations from a strict power-law evolution are seen. We do not find any clear-cut evidence for systematically frequency-dependent evolution, suggesting that the radii as determined from visibility-plane model-fitting continue to provide reasonable estimates of the physical outer shock-front radius.

We present multi-epoch global VLBI observations at 2 cm, 3.6 cm and 6 cm of the compact radio sources in Arp220. We resolve many sources and estimate sizes, expansion velocities and source classes. We find most source properties are consistent with them being either radio supernovae (SNe) or supernova remnants (SNRs). We extend the luminosity-diameter relation for SNRs to very small sources and argue this supports models where shell magnetic fields are internally amplified. We also detect one probable SN/SNR transition object candidate and three highly variable sources with possible superluminal motion (approximately 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources, which show similarities to the recently discovered superluminal source in M82, might be associated with an AGN or a new radio source class (e.g. intermediate-mass black holes or micro-blazars).

We present results from 4.8 GHz VLA and Global-VLBI observations of the northern half of the moderate FIR luminosity (median L_IR = 10^11.01 L_Sol) COLA sample of star-forming galaxies. VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosities of these cores (~10^21 W/Hz) are too large to be explained by radio supernovae or supernova remnants and we argue that they are instead powered by AGN. These sub-parsec scale radio cores are preferentially detected toward galaxies whose VLA maps show bright 100-500 parsec scale nuclear radio components. Since these latter structures tightly follow the FIR to radio-continuum correlation for star-formation we conclude that the AGN powered VLBI sources are associated with compact nuclear starburst environments. The implications for possible starburst-AGN connections are discussed. The detected VLBI sources have a relatively narrow range of radio luminosity consistent with models in which intense compact Eddington-limited starbursts regulate the gas supply onto a central super-massive black hole. The high incidence of AGN radio cores in compact starbursts suggests little or no delay between the starburst phase and the onset of AGN activity.

We present Very Long Baseline Array (VLBA) HI absorption observations of the core region of the powerful radio galaxy Cygnus A. These data show both broad (FWHM = 231 \pm 21 km/s) and narrow (FWHM <30 km/s) velocity width absorption components. The broad velocity absorption shows high opacity on the counter jet, low opacity against the core and no absorption on the jet side. We argue that these results are most naturally explained by a circumnuclear HI absorbing disk orientated roughly perpendicular to the jet axis. We estimate that the HI absorbing gas lies at a radius of ~80 pc has a scale height of about 20 pc, density n > 10^{4} cm^{-3} and total column density in the range 10^{23}-10^{24} cm^{-2}. Models in which the HI absorption is primarily from an atomic or a molecular gas phase can both fit our data. Modelling taking into account the effective beam shows that the broad HI absorbing gas component does not cover the radio core in Cygnus A and therefore does not contribute to the gas column that blocks our view of the hidden quasar nucleus. If however Cygnus A were observed from a different direction, disk gas on ~100 pc radius scales would contribute significantly to the nuclear column density, implying that in some radio galaxies gas on these scales may contribute to the obscuration of the central engine. We argue that the circumnuclear torus in Cygnus A contains too little mass to power the AGN over > 10^{7} yr but that material in the outer HI absorbing gas disk can provide a reservoir to fuel the AGN and replenish torus clouds. The second narrow HI absorption component is significantly redshifted (by 186km/s) with respect to the systemic velocity and probably traces infalling gas which will ultimately fuel the source. [abridged]

Radio astronomy is entering the era of large surveys. This paper describes the plans for wide surveys with the LOw Frequency ARray (LOFAR) and their synergy with large surveys at higher frequencies (in particular in the 1-2 GHz band) that will be possible using future facilities like Apertif or ASKAP. The LOFAR Survey Key Science Project aims at conducting large-sky surveys at 15, 30, 60, 120 and 200 MHz taking advantage of the wide instantaneous field of view and of the unprecedented sensitivity of this instrument. Four topics have been identified as drivers for these surveys covering the formation of massive galaxies, clusters and black holes using z>6 radio galaxies as probes, the study of the intercluster magnetic fields using diffuse radio emission and Faraday rotation measures in galaxy clusters as probes and the study of star formation processes in the early Universe using starburst galaxies as probes. The fourth topic is the exploration of new parameter space for serendipitous discovery taking advantage of the new observational spectral window open up by LOFAR. Here, we briefly discuss the requirements of the proposed surveys to address these (and many others!) topics as well as the synergy with other wide area surveys planned at higher frequencies (and in particular in the 1-2 GHz band) with new radio facilities like ASKAP and Apertif. The complementary information provided by these surveys will be crucial for detailed studies of the spectral shape of a variety of radio sources (down to sub-mJy sources) and for studies of the ISM (in particular HI and OH) in nearby galaxies.

We report ~100 pc (0.3") resolution observations of (sub)millimeter HCO+ and CO lines in the ultraluminous infrared galaxy Arp 220. The lines peak at two merger nuclei, with HCO+ being more spatially concentrated than CO. Asymmetric line profiles with blueshifted absorption and redshifted emission are discovered in HCO+(3-2) and (4-3) toward the two nuclei and in CO(3-2) toward one nucleus. We suggest that these P-Cygni profiles are due to ~100 km/s outward motion of molecular gas from the nuclei. This gas is most likely outflowing from the inner regions of the two nuclear disks rotating around individual nuclei, clearing the shroud around the luminosity sources there.

A strong signature of a circumstellar disc around a high-mass protostar has been inferred from high resolution methanol maser observations in NGC7538-IRS1 N. This interpretation has however been challenged with a bipolar outflow proposed as an alternative explanation. We compare the two proposed scenarios for best consistency with the observations. Using a newly developed formalism we model the optical depth of the maser emission at each observed point in the map and LOS velocity for the two scenarios. We find that if the emission is symmetric around a central peak in both space and LOS velocity then it has to arise from an edge-on disc in sufficiently fast differential rotation. Disc models successfully fit ~100 independent measurement points in position-velocity space with 4 free parameters to an overall accuracy of 3-4%. Solutions for Keplerian rotation require a central mass of at least 4 solar masses. Close to best-fitting models are obtained if Keplerian motion is assumed around a central mass equaling ~30 solar masses as inferred from other observations. In contrast we find that classical bipolar outflow models cannot fit the data, although could be applicable in other sources. Our results strongly favour the differentially rotating disc hypothesis to describe the main feature of the 12.2 (and 6.7) GHz methanol maser emission in NGC7538 IRS1 N. Furthermore, for Keplerian rotation around a ~30 solar masses protostar we predict the position and velocity at which tangentially amplified masers should be detected in high dynamic range observations. [abridged]

We report the first detection at multiple radio wavelengths (13, 6 and 3.6 cm) of the compact sources within both nuclei of the Ultra Luminous Infra-Red Galaxy Arp 220. We present the radio spectra of the 18 detected sources. In just over half of the sources we find that these spectra and other properties are consistent with the standard model of powerful Type IIn supernovae interacting with their pre-explosion stellar wind. The rate of appearance of new radio sources identified with these supernova events suggests that an unusually large fraction of core collapse supernovae in Arp 220 are highly luminous; possibly implying a radically different stellar initial mass function or stellar evolution compared to galactic disks. Another possible explanation invokes very short (~3 x 10^5 year) intense (~10^3 M_Sol year^-1) star formation episodes with a duty cycle of ~10%. A second group of our detected sources, consisting of the brightest and longest monitored sources at 18 cm do not easily fit the radio supernova model. These sources show a range of spectral indexes from -0.2 to -1.9. We propose that these are young supernova remnants which have just begun interacting with a surrounding ISM with a density between 10^4 and 10^5 cm^-3. One of these sources is probably resolved at 3.6 cm wavelength with a diameter 0.9 pc. In the western nucleus we estimate that the ionized component of the ISM gives rise to foreground free-free absorption with opacity at 18 cm of <0.6 along the majority of lines of sight. Other sources may be affected by absorption with opacity in the range 1 to 2. These values are consistent with previous models as fitted to the radio recombination lines and the continuum spectrum.

Context: Methanol masers at 6.7GHz are recognised markers of high-mass star formation regions. The study of their distribution in the Galaxy gives important insights into the star formation activity of the Milky Way. We present a statistical analysis on the General Catalogue of 6.7GHz methanol masers in the Galaxy with the aim of extracting global properties of the masers. Aims: We provide constraints on the luminosity function of 6.7GHz methanol masers and on their total number in the Galaxy. Methods: We model the spatial distribution of the masers in the Milky Way by using their distribution in galactocentric distance which is unambiguous once a rotation curve for the Galaxy is assumed. This is the starting point for determining the luminosity function of the masers. Results: The luminosity function of 6.7GHz methanol masers is modelled as a power-law with sharp cutoffs and having an index lying between -1.5 and -2. We also predict the number of detections of methanol masers assuming different sensitivity limits in the observations.

Context: NGC7538 is known to host a 6.7 and 12.2GHz methanol maser cospatial with a Ultra Compact (UC) HII region, IRS1. Aims: We report on the serendipitous discovery of two additional 6.7GHz methanol masers in the same region, not associated with IRS1. Methods: Interferometry maser positions are compared with recent single-dish and interferometry continuum observations. Results: The positions of the masers agree to high accuracy with the 1.2mm continuum peak emission in NGC7538 IRS9 and NGC7538 S. This clear association is also confirmed by the positional agreement of the masers with existing high resolution continuum observations at cm and/or mm wavelengths. Conclusions: Making use of the established strong relation between methanol masers and high-mass star formation, we claim that we have accurately positioned the high-mass protostars within the regions where they are detected. The variety of objects hosting a 6.7GHz methanol maser in NGC7538 shows that this emission probably traces different evolutionary stages within the protostellar phase.

We model the OH megamaser emission from the luminous infrared galaxy IIIZw35 as arising from a narrow rotating starburst ring of radius 22 pc enclosing a mass of 7 x 10^6 Msun. We show how both the compact and apparently diffuse maser emission from this ring can arise from a single phase of unsaturated maser clouds amplifying background radio continuum. The masing clouds are estimated to have a diameter of <0.7 pc and internal velocity dispersion of about 20 km/s. We find that the clouds are neither self-gravitating nor pressure confined but are freely expanding. Their dispersal lifetimes may set the vertical thickness of the ring. For an estimated internal density of 3 x 10^3 cm^-3, cloud masses are of order 24 Msun. The observed spectral features and velocity gradients indicate that the clouds must be outflowing and escaping the nucleus. The cloud mass outflow rate is estimated to be 0.8 Msun/yr, while the star formation rate is about 19 Msun/yr. Associated ionised gas, possibly generated from dissipated clouds, provides free-free absorption along the source axis, explaining the observed East-West asymmetries. We show that the clumpiness of a maser medium can have a dramatic effect on what is observed even in a relatively low gain OH megamaser. Specifically, in IIIZw35 our clumpy maser model naturally explains the large line to continuum ratios, the large 1667MHz:1665MHz line ratios and the wide velocity dispersions seen in the compact maser spots. Other astrophysical masers showing both compact and apparently diffuse emission might be explained by similar clumpy structures.

Radio imaging of ULIR galaxies is ideal to explore the connection between the starburst and the AGN phenomenon since it is unaffected by dust obscuration, and provides the required high angular resolution to distinguish between an AGN and starburst emission. We have made combined 18cm radio continuum, EVN and MERLIN observations of 13 ULIRGs that have the parsec and deci-parsec scale resolution necessary to distinguish between an AGN and supernovae remnants at the centres of these galaxies, and assess the contribution of each to the total energy distribution. Images of four galaxies are presented here.

We present the results of simultaneous high resolution observations of v=1 and v=2, J=1-0 SiO masers toward TX Cam at four epochs covering a stellar cycle. Near maser maximum (Epochs III and IV), the individual components of both masers are distributed in ring-like structures but the ring is severely disrupted near stellar maser minimum (Epochs I and II). In Epochs III and IV there is a large overlap between the radii at which the two maser transitions occur. However in both epochs the average radius of the v=2 maser ring is smaller than for the v=1 maser ring, the difference being larger for Epoch IV. The observed relative ring radii in the two transitions, and the trends on the ring thickness, are close to those predicted by the model of Humphreys et al. (\cite{humphreys02}). In many individual features there is an almost exact overlap in space and velocity of emission from the two transitions, arguing against pure radiative pumping. At both Epochs III and IV in many spectral features only 50% of the flux density is recovered in our images, implying significant smooth maser structure. For both transitions we find that red- and blue-shifted masers occur in all parts of the rings, with relatively few masers at the systemic velocity. Thus there is no evidence for rotation, although the blue-shifted masers are somewhat more prominent to the west. At all four epochs red-shifted components are generally brighter than blue-shifted ones. At Epochs III and IV, we see many filamentary or spoke-like features in both v=1 and v=2 masers, especially in the red-shifted gas. These spokes show systematic velocity gradients consistent with a decelerating outward flow with increasing radius. We outline a possible model to explain why, given the presence of these spokes, there is a deficit of maser features at the systemic velocity.

We summarize the present status of VLBI experiments at 3 mm (86 GHz), 2 mm (129-150 GHz) and 1.3 mm (215-230 GHz). We present and discuss a new 3 mm VLBI map of M87 (Virgo A), which has a spatial resolution of only approx. 20 Schwarzschild radii. We discuss recent VLBI results for SgrA* and argue in favor of new observations within an extended European mm-VLBI network, in order to search for variability. We discuss the possibilities to image the `event horizon' of a super-massive black hole at wavelengths < 2mm, and conclude that the addition of large and sensitive millimetre telescopes such as CARMA, the SMA, the LMT and ALMA will be crucial for this.

Global millimetre VLBI allows detailed studies of the most central jet regions of AGN with unprecedent spatial resolution of a few 100-1000 Schwartzschild radii to be made. Study of these regions will help to answer the question how the highly relativistic AGN jets are launched and collimated. Since the early 1990s, bright mm-sources have been observed with global 3 mm VLBI. Here we present new images from an ongoing systematic analysis of the available observations. In particular, we focus on the structure and structural evolution of the best observed AGN jets, taking 3C 454.3 as a characteristic example. This core-dominated and highly variable quasar shows a complex morphology with individual jet components accelerating superluminally towards the outer structure. We briefly discuss the X-ray properties of 3C 454.3 and present its radio- to X-ray large-scale brightness distribution.

We present an edge-on Keplerian disk model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disk. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating disks. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disk to have inner and outer radii of about 270 AU and 750 AU.

We present a model for the OH megamaser emission in the starburst galaxy IIIZw35. The observed diffuse and compact OH maser components in this source are explained by a single phase of unsaturated clumpy gas distributed in a thin ring structure and amplifying background continuum. We emphasize the importance of clumpiness in the OH masing medium, an effect that has not been fully appreciated previously. The model explains why multiple bright spots are seen only at the ring tangents while smoother emission is found elsewhere. Both the observed velocity gradients and the line to continuum ratios around the ring enquire a geometry where most of the seed photons come from a continuum emission which lies outside the OH ring. To explain both the OH and continuum brightness, free-free absorbing gas is required along the ring axis to partially absorb the far side of the ring. It is proposed that the required geometry arises from an inwardly propagating ring of starburst activity.

We discuss the propagation of spectral line and continuum radiation in a clumpy medium and give general expressions for the observed absorption or emission from a cloud population. We show that the affect of the medium clumpiness can usually be characterised by a single number multiplying the mean column opacity. Our result provides a simpler proof and generalisation of the result of Martin et al (1984). The formalism provides a simple way to understand the effects of clumping on molecular line profiles and ratios, for example how clumping effects the interpretation of 13CO(1-0) to 12CO(1-0) line ratios. It can also be used as a propagation operator in physical models of clumpy media where the incident radiation effects the spectral line emissivity. We are working to extend the formalism to the propagation of masers in a clumpy medium, but in this case there are special difficulties because formal expectation values are not characteristic of observations because they are biased by rare events.

We present an edge-on Keplerian disc model to explain the main component of the 12.2 and 6.7 GHz methanol maser emission detected toward NGC7538-IRS1 N. The brightness distribution and spectrum of the line of bright masers are successfully modeled with high amplification of background radio continuum emission along velocity coherent paths through a maser disc. The bend seen in the position-velocity diagram is a characteristic signature of differentially rotating discs. For a central mass of 30 solar masses, suggested by other observations, our model fixes the masing disc to have inner and outer radii of about 350 AU and 1000 AU.

We report observations of the T Tauri system at 8.4 GHz with a VLBI array comprising the VLBA, VLA and Effelsberg 100m telescopes. We detected a compact source offset approximately 40 mas from the best infrared position of the T Tau Sb component. This source was unresolved, and constrained to be less than 0.5 mas in size, corresponding to 0.07 AU or 15 R$_{\odot}$ at a distance of 140 pc. The other system components (T Tau Sa, T Tau N) were not detected in the VLBI data. The separate VLA map contains extended flux not accounted for by the compact VLBI source, indicating the presence of extended emission on arcsecond scales. The compact source shows rapid variability, which together with circular polarization and its compact nature indicate that the observed flux arises from a magnetically-dominated region. Brightness temperatures in the MK range point to gyrosynchrotron as the emission mechanism for the steady component. The rapid variations are accompanied by dramatic changes in polarization, and we record an at times 100% polarized component during outbursts. This strongly suggests a coherent emission process, most probably an electron cyclotron maser. With this assumption it is possible to estimate the strength of the local magnetic field to be 1.5-3 kilogauss.

We consider the incidence of HI absorption in intrinsically small sub-galactic sized extragalactic sources selected from sources classified as Gigahertz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources. We find that the smaller sources (<0.5 kpc) have larger HI column densities than the larger sources (>0.5 kpc). Both a spherical and an axi-symmetric gas distribution, with a radial power law density profile, can be used to explain this anti-correlation between projected linear size and HI column density. Since most detections occur in objects classified as galaxies, we argue that if the unified schemes apply to GPS/CSSs a disk distribution for the HI is more likely. The most favoured explanation for the compact sizes of the GPS/CSSs is that they are young sources evolving in a power law density medium. For the GPSs with measured expansion velocities, our derived densities are within an order of magnitude of those estimated from ram-pressure confinement of the lobes assuming equipartition. Our results therefore support the youth model.

We present an overview of the occurrence and properties of atomic gas associated with compact radio sources at redshifts up to z=0.85. Searches for HI 21cm absorption were made with the Westerbork Synthesis Radio Telescope at UHF-high frequencies (725-1200 MHz). Detections were obtained for 19 of the 57 sources with usable spectra (33%). We have found a large range in line depths, from tau=0.16 to tau<=0.001. There is a substantial variety of line profiles, including Gaussians of less than 10km/s, to more typically 150km/s, as well as irregular and multi-peaked absorption profiles, sometimes spanning several hundred km/s. Assuming uniform coverage of the entire radio source, we obtain column depths of atomic gas between 1e19 and 3.3e21(Tsp/100K)(1/f)cm^(-2). There is evidence for significant gas motions, but in contrast to earlier results at low redshift, there are many sources in which the HI velocity is substantially negative (up to v=-1420km/s) with respect to the optical redshift, suggesting that in these sources the atomic gas, rather than falling into the centre, may be be flowing out, interacting with the jets, or rotating around the nucleus.

Combining our own observations with data from the literature, we consider the incidence of HI absorption in Gigahertz Peaked Spectrum (GPS) and Compact Steep Spectrum (CSS) sources. Here we present our preliminary results, where we find that the smaller GPS sources (<1 kpc) on average have larger HI column densities than the larger CSS sources (>1 kpc). Both a spherical and an axi-symmetric gas distribution, with a radial power law density profile, can be used to explain this anti-correlation between projected linear size and HI column density. Since most detections occur in galaxy classified objects, we argue that if the unified schemes apply to the GPS/CSSs, a disk distribution for the HI is more likely.

We report on the results of a two epoch study of the low power Compact Symmetric Object 4C 31.04. Observations performed with the VLBA at 5 GHz in 1995 and 2000 have yielded images of this source at milliarcsecond angular resolution. A central core is detected, with bright compact hot spots and extended lobes on both sides. Model-fitting and other analysis of the data (brightness profile, difference map) clearly indicate that the source is expanding. We estimate the velocity of this expansion to be (0.085 +/- 0.016) mas/yr, i.e. (0.33 +/- 0.06) c in both hot spots. Assuming a constant expansion velocity, we estimate the kinematic age of the source at 550 yrs. We also study the spectral index using VLBA observations at 1.3 GHz and MERLIN at 22 GHz. The derived spectral age is 3000-5000 years in equipartition conditions. The two estimates are discussed and found to be in agreement, given present uncertainties.

We discuss recent measurements of proper motions of the hotspots of Compact Symmetric Objects. Source expansion has been detected in ten CSOs so far and all these objects are very young (<3000 years). In a few sources ages have also been estimated from energy supply and spectral ageing arguments and these estimates are comparable. This argues that these sources are close to equipartition and that standard spectral ageing models apply. Proper motions studies are now constraining hotspot accelerations, side-to-side motions and differences in hotspot advance speeds between the two hotspots within sources. Although most CSOs are young sources their evolution is unclear. There is increasing evidence that in some objects the CSO structure represents a new phase of activity within a recurrent source.

PKS 1413+135 (z=0.24671) is one of very few radio-loud AGN with an apparent spiral host galaxy. Previous authors have attributed its nearly exponential infrared cutoff to heavy absorption but have been unable to place tight limits on the absorber or its location in the optical galaxy. In addition, doubts remain about the relationship of the AGN to the optical galaxy given the observed lack of re-emitted radiation. We present new HST, ASCA and VLBA observations which throw significant new light on these issues. The HST observations reveal an extrremely red color (V-H = 6.9 mag) for the active nucleus of PKS 1413+135, requiring both a spectral turnover at a few microns due to synchrotron aging and a GMC-sized absorber. We derive an intrinsic column N_H = 4.6^{+2.1}_{-1.6} times 10^{22}cm^{-2} and covering fraction f = 0.12^{+0.07}_{-0.05}. As the GMC is likely in the disk of the optical galaxy, our sightline is rather unlikely (P ~ 2 times 10^{-4}). The properties of the GMC typical of GMCs in our own galaxy. The HI absorber appears centered 25 milliarcseconds away from the nucleus, while the X-ray and nearly all of the molecular absorbers must cover the nucleus, implying a complicated geometry and cloud structure, with a molecular core along our line of sight to the nucleus. Interestingly, the HST/NICMOS data require the AGN to be decentered relative to the optical galaxy by 13 +/- 4 milliarcseconds. This could be interpreted as suggestive of an AGN location far in the background compared to the optical galaxy, but it can also be explained by obscuration and/or nuclear structure, which is more consistent with the observed lack of multiple images.

At 147 GHz (2mm wavelength), we detected three prominent AGN (NRAO150, 3C279, 1633+382) with Very Long Baseline Interferometry (VLBI) with an angular resolution of only 18 micro-arcseconds on the baseline between two antennas in Arizona (10m HHT and 12m KittPeak) and the IRAM 30m antenna on Pico Veleta in Spain. This is a new world record in radio interferometry and astronomical imaging and opens fascinating future possibilities to directly image and study the innermost regions in Quasars and other Active Galactic Nuclei.

We review the state of the Onsala blind survey of the galactic plane, searching for new 6.7 GHz methanol masers. We also describe preliminary results of millimeter follow-up observations of the new detections and high resolution observations using the EVN. We conclude that blind surveys are important to complement targeted searches done until now and give the possibility to detect new classes of objects.

We present 3-band HST imaging of three z</=0.1 compact symmetric objects (CSOs): 4C 31.04, 1946+708 and 1146+596. These objects were chosen on the basis of proximity to Earth as well as HI 21 cm absorption. The inner H-band isophotes of these galaxies are well fit by "nuker" models, typical of nearby ellipticals. Each shows a significant flattening in the isophotal profile at radii ~0.5", as well as variations in ellipticity and PA. However, neither is uncommon for elliptical galaxies. All three show modest departures from nuker models at radii of 1-5 $h_{60}^{-1}$ kpc. Each shows large, well distributed dust features, which are somewhat concentrated in the nuclear regions in features which resemble disks or tori. The amount of dust in these galaxies is about 10 times higher than normal for ellipticals and radio galaxy hosts. The nuclear dust disks tend to be oriented roughly perpendicular to the radio axis. One galaxy, 4C31.04, exhibits bright nuclear regions well-aligned with the radio axis, while another, 1146+596, shows a significant near-IR excess resembling a stellar bar along its dust disk. The combination of outwardly normal isophotal profiles with significant variations in PA and ellipticity is consistent with the host galaxies being relatively recent merger remnants, with the mergers having occurred >/= 10^8 years ago. Such a merger could have "triggered" the current activity in these objects, but our data require a significant time delay between the merger and the onset of nuclear activity. However, these data are also consistent with the hypothesis that the onset of nuclear activity in radio galaxies is due to relatively minor "feeding" events and/or the formation of "bars within bars", which would disturb the internal kinematics only slightly.

We present measurements of absorption by the 21cm hyperfine transition of neutral hydrogen toward radio sources at substantial look-back times. These data are used in combination with observations of rotational transitions of common interstellar molecules, to set limits on the evolution of the fine structure constant: (dot alpha)/alpha < 3.5e-15/yr to a look-back time of 4.8 Gyr. The neutral hydrogen observations employed Very Long Baseline Interferometry in order to mitigate the substantial uncertainty arising from the fact that observations at very different wavelengths may probe different lines-of-site due to frequency dependent structure of the background source. We discuss the implication of these results on theories unifying natural forces based on compact quantum dimensions. In the context of string theory, the limit on the secular evolution of the scale factor of these compact dimensions, R, is (dot R)/R < 1e-15/yr. Including terrestrial and other astronomical measurements places limits (2sigma) on slow oscillations of R from the present to the epoch of cosmic nucleosynthesis, just seconds after the big bang, of (Delta R)/R < 1e-5.

We present VLBI and MERLIN data at 1.66 and 4.99 GHz on the central component coincident with the nucleus of the giant radio galaxy, 3C236. The nuclear radio structure is composed of two complexes of emission which are resolved on scales from 1 milli-arcsec (mas) to 1 arcsec. Oscillations with an amplitude of about 5 degrees can be seen in the compact radio structure. Spectral index distributions are plotted at angular resolutions of 10 and 25 mas and allow us to identify the core component in the south-east emission complex. Re-examination of the HST WFPC-2 image of 3C236 by de Koff et al. (2000), shows that the normal to the dust disk in the nucleus is ~30 degrees from the plane of the sky and within 12 degrees of parallel to the overall orientation of the radio source. We suggest that the radio axis is also at an angle of ~30 degrees to the plane of the sky and that the north-west jet is on the approaching side. This orientation implies an overall size of 4.5 Mpc (H_o = 75 km/s/Mpc, q_o = 0.5) for 3C236. The coincidence of a dust feature and the south-east compact jet, within the astrometric errors, leads us to suggest that the dust may be in the form of a cloud encountered by the jet in the first ~400 pc of its journey out from the nucleus. One-sided emission at 5 GHz on 1 mas scales would suggest that the jets are ejected initially at <= 35 degrees to the line of sight, but this is difficult to reconcile with the obvious orientation stability of the jet system as a whole. Free-free absorption of the counter-jet may be an alternative explanation for the one-sideness.

We obtained EVN observations of the X-ray transient source GRS1915+105 in its radio-loud state. In 2 images obtained before a major flare, the source is resolved, showing two opposite asymmetric jets. This indicates that steady jets exist also before a major flare, with an estimated velocity of 0.2-0.6c.

We report on high sensitivity, spectral line VLBI observations of the HI absorption feature in the radio galaxy NGC4261. Although absorption is only detectable on the most sensitive baseline, it can be unambiguously associated with the counterjet and is interpreted to originate in a thin atomic circumnuclear disk. This structure is probably a continuation of the dusty accretion disk inferred from HST imaging, which could be feeding the massive black hole. HI column densities in front of the counterjet of the order of 10^{21}(T_sp/100 K) cm^{-2} are derived, consistent with X-ray data and VLBI scale free-free absorption. The data presented here are the result of the first scientific project processed on the new EVN MkIV data processor.

We report on 3.6 cm VLA and VLBA observations of YZ CMi and AD Leo, two nearby dMe stars. We resolve YZ CMi and can fit a circular symmetrical gaussian component of FWHP of 0.98 +/-0.2 mas, corresponding to an extent of the corona above the photosphere of 1.77 x 10^{10} +/-8.8 x 10^{9} cm or 0.7 +/-0.3 Rstar (Rstar refers to the photospheric radius). We obtain an estimate of the brightness temperature of 7.3 x 10^{7} K, which is consistent with that expected from gyrosynchrotron emission. For AD Leo the emitting region is unresolved. We therefore set a conservative upper limit to its diameter of 1.8 times the photosphere diameter, which leads to an extent of the corona above the photosphere of <2.8 x 10^{10} cm or <0.8 Rstar. We compare the radio emitting dMe stars with measured sizes with the Sun and conclude that these active stars have much more extended coronal radio emission than the Sun. The VLBA position of YZ CMi has been found to differ by 32 mas from the positions calculated from the Hipparcos catalogue. The discrepancy is caused by large errors in the listed proper motion. An improved value is given.

We present an analysis of multi-epoch global VLBI observations of the Compact Symmetric Objects: 2352+495 and 0710+439 at 5 GHz. Analysis of data spread over almost two decades shows strong evidence for an increase in separation of the outer components of both sources at a rate of ~0.2 h^{-1} c (for q_{0}=0.5 and H_{0}=100 h kms^{-1} Mpc^{-1}). Dividing the overall sizes of the sources by their separation rates implies that these Compact Symmetric Objects have a kinematic age < 10^{4} years. These results (and those for other CSOs) strongly argue that CSOs are indeed very young sources and that they are probably evolve into the much larger classical doubles.

We present results of VLA and VLBA observations of the 1.420 GHz neutral hydrogen absorption associated with the Compact Symmetric Object 1946+708 (z=0.101). We find significant structure in the gas on parsec scales. The peak column density in the HI (N_HI~2.2x10^23 cm^-2 (T_s/8000K)) occurs toward the center of activity of the source, as does the highest velocity dispersion (FWHM~350 \kms). In addition, we find that the continuum spectra of the various radio components associated with these jets strongly indicate free-free absorption. This effect is particularly pronounced toward the core and inner components of the receding jet, suggesting the presence of a screen local to the source, perhaps part of an obscuring torus.

We present the results of multi-epoch global VLBI observations of the Compact Symmetric Object (CSO), 0710+439 at 5 GHz. Analysis of data spread over 13 years shows strong evidence for an increase in the separation of the outer components at a rate of $0.251 \pm 0.029 h^{-1}c$. Dividing the overall size of 86.8 $h^{-1}$pc by this separation rate implies an estimated kinematic age of only $1100 \pm 100$ yrs. After taking into account possible temporal variations in hotspot advance speeds due to cloud collisions or hydrodynamic instabilities we argue that the upper limit to the age of 0710+439 is most likely within a factor of 2 of this estimate and certainly within a factor of 10 (i.e. $<11 000$ yrs). This result therefore strongly supports the idea that Compact Symmetric Objects are very young radio-loud sources. Furthermore the large radiative efficiency we calculate for 0710+439 is consistent with strong negative luminosity evolution as CSOs grow in size and with them evolving into classical double sources.

We present the results of multi-epoch global VLBI observations of the Compact Symmetric Object (CSO) 0108+388 at 5 GHz. Analysis of data spread over 12 years shows strong evidence for an increase in the separation of the outer components at a rate of $0.197\pm0.026 h^{-1}c$. Given an overall size of 22.2 $h^{-1}$pc this implies a kinematic age of only $367\pm48$ yrs. This result strongly supports the idea that radio emission in Compact Symmetric Objects arises from recently activated radio sources. The presence of weak radio emission on kpc-scales in 0108+388 suggests recurrent activity in this source, and that we are observing it just as a new period of activity is beginning.