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

We performed simultaneous observations of the H2O 6(1,6) - 5(2,3) (22.235080 GHz) and SiO v= 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1, 3 - 2 (43.122080, 42.820587, 86.243442, and 129.363359 GHz) masers towards the suspected D-type symbiotic star, V627 Cas, using the Korean VLBI Network. Here, we present astrometrically registered maps of the H2O and SiO v = 1, 2, J = 1 - 0, SiO v = 1, J = 2 - 1 masers for five epochs from January 2016 to June 2018. Distributions of the SiO maser spots do not show clear ring-like structures, and those of the H2O maser are biased towards the north-north-west to west with respect to the SiO maser features according to observational epochs. These asymmetric distributions of H2O and SiO masers are discussed based on two scenarios of a bipolar outflow and the presence of the hot companion, a white dwarf, in V627 Cas. We carried out ring fitting of SiO v = 1, and v = 2 masers and estimated the expected position of the cool red giant. The ring radii of the SiO v = 1 maser are slightly larger than those of the SiO v = 2 maser, as previously known. Our assumption for the physical size of the SiO maser ring of V627 Cas to be the typical size of a SiO maser ring radius (\sim4 au) of red giants yields the distance of V627 Cas to be \sim1 kpc.

High precision astrometry provides the foundation to resolve many fundamental problems in astrophysics. The application of astrometric studies spans a wide range of fields, and has undergone enormous growth in recent years. This is as a consequence of the increasing measurement precision and wide applicability, which is due in turn to the development of new techniques. Forthcoming next generation observatories have the potential to further increase the astrometric precision, providing there is a matching improvement in the methods to correct for systematic errors. The EVN and other observatories are providing demonstrations of these and are acting as pathfinders for next-generation telescopes such as the SKA and ngVLA. We will review the perspectives for the coming facilities and examples of the current state-of-the-art for astrometry.

The SKA will deliver orders of magnitude increases in sensitivity, but most astrometric VLBI observations are limited by systematic errors. In these cases improved sensitivity offers no benefit. The best current solution for improving the accuracy of the VLBI calibration is \MV\ VLBI, where multiple simultaneous observations around the target are used to deduce the corrections required for the line of sight to the target. We have estimated and quantified the applicability of \MV\ from real-world ionospheric studies, making projections into achievable astrometric accuracies. These predict systematic measurement errors, with calibrators separated by several degrees, of $\sim$10\uas\ with current VLBI facilities. For closer calibrators, that are in-beam for single dish VLBI facilities, we predict systematic measurement errors of a few \uas. This is the ideal combination, where the sensitivity of the SKA will provide the precision and \MV\ will provide the accuracy. Based on these results we suggest that the SKA design should increase the number of VLBI beams it can form from four to as many as ten.

The KVN(Korean VLBI Network)-style simultaneous multi-frequency receiving mode is demonstrated to be promising for mm-VLBI observations. Recently, other Very long baseline interferometry (VLBI) facilities all over the globe start to implement compatible optics systems. Simultaneous dual/multi-frequency VLBI observations at mm wavelengths with international baselines are thus possible. In this paper, we present the results from the first successful simultaneous 22/43 GHz dual-frequency observation with KaVA(KVN and VERA array), including images and astrometric results. Our analysis shows that the newly implemented simultaneous receiving system has brought a significant extension of the coherence time of the 43 GHz visibility phases along the international baselines. The astrometric results obtained with KaVA are consistent with those obtained with the independent analysis of the KVN data. Our results thus confirm the good performance of the simultaneous receiving systems for the non-KVN stations. Future simultaneous observations with more global stations bring even higher sensitivity and micro-arcsecond level astrometric measurements of the targets.

We obtained, for the first time, astrometrically registered maps of the 22.2 GHz H2O and 42.8, 43.1, and 86.2 GHz SiO maser emission toward the semiregular b-type variable (SRb) R Crateris, at three epochs (2015 May 21, and 2016 January 7 and 26) using the Korean Very-long-baseline Interferometry Network. The SiO masers show a ring-like spatial structure, while the H2O maser shows a very asymmetric one-side outflow structure, which is located at the southern part of the ring-like SiO maser feature. We also found that the 86.2 GHz SiO maser spots are distributed in an inner region, compared to those of the 43.1 GHz SiO maser, which is different from all previously known distributions of the 86.2 GHz SiO masers in variable stars. The different distribution of the 86.2 GHz SiO maser seems to be related to the complex dynamics caused by the overtone pulsation mode of the SRb R Crateris. Furthermore, we estimated the position of the central star based on the ring fitting of the SiO masers, which is essential for interpreting the morphology and kinematics of a circumstellar envelope. The estimated stellar coordinate corresponds well to the position measured by Gaia.

The ambitious scientific goals of the SKA require a matching capability for calibration of atmospheric propagation errors, which contaminate the observed signals. We demonstrate a scheme for correcting the direction-dependent ionospheric and instrumental phase effects at the low frequencies and with the wide fields of view planned for SKA-Low. It leverages bandwidth smearing, to filter-out signals from off-axis directions, allowing the measurement of the direction-dependent antenna-based gains in the visibility domain; by doing this towards multiple directions it is possible to calibrate across wide fields of view. This strategy removes the need for a global sky model, therefore all directions are independent. We use MWA results at 88 and 154 MHz under various weather conditions to characterise the performance and applicability of the technique. We conclude that this method is suitable to measure and correct for temporal fluctuations and direction-dependent spatial ionospheric phase distortions on a wide range of scales: both larger and smaller than the array size. The latter are the most intractable and pose a major challenge for future instruments. Moreover this scheme is an embarrassingly parallel process, as multiple directions can be processed independently and simultaneously. This is an important consideration for the SKA, where the current planned architecture is one of compute-islands with limited interconnects. Current implementation of the algorithm and on-going developments are discussed.

The supergiant VX Sagittarii is a strong emitter of both H2O and SiO masers. However, previous VLBI observations have been performed separately, which makes it difficult to spatially trace the outward transfer of the material consecutively. Here we present the astrometrically registered, simultaneous maps of 22.2 GHz H2O and 43.1/42.8/86.2/129.3 GHz SiO masers toward VX Sagittarii. The H2O masers detected above the dust-forming layers have an asymmetric distribution. The multi-transition SiO masers are nearly circular ring, suggesting spherically symmetric wind within a few stellar radii. These results provide the clear evidence that the asymmetry in the outflow is enhanced after the smaller molecular gas clump transform into the inhomogeneous dust layers. The 129.3 GHz maser arises from the outermost region compared to that of 43.1/42.8/86.2 GHz SiO masers. The ring size of the 129.3 GHz maser is maximized around the optical maximum, suggesting that radiative pumping is dominant.

We report on the astrometric registration of VLBI images of the SiO and water masers in OH231.8+4.2, the iconic Proto-Planetary Nebula also known as the Calabash nebula, using the KVN and Source/Frequency Phase Referencing. This, for the first time, robustly confirms the alignment of the SiO masers, close to the AGB star, driving the bi-lobe structure with the water masers in the out-flow. We are able to trace the bulk motions for the water masers over the last few decades to be 19 km/s and deduce that the age of this expansion stage is 38$\pm$2 years. The combination of this result with the distance allows a full 3D reconstruction, and confirms that the water masers lie on and expand along the known large-scale symmetry axis and that the outflow is only a few decades old, so mass loss is almost certainly on-going. Therefore we conclude that the SiO emission marks the stellar core of the nebular, the water emission traces the expansion, and that there must be multiple epochs of ejection to drive the macro-scale structure.

We present a multi-calibrator solution, i.e. MultiView, to achieve accurate astrometry on the level of the thermal noise at low VLBI frequencies dominated by ionospheric residuals. We demonstrate on L-band VLBA observations how MultiView provides superior astrometry to conventional phase referencing techniques (Rioja et al. 2017). We also introduce a new trial method to detect antenna based systematic errors in the observations (Orosz et al. 2017}. All presented methods and results are based on our recent papers (Orosz et al. 2017; Rioja et al. 2017).

Atmospheric propagation effects at millimeter wavelengths can significantly alter the phases of radio signals and reduce the coherence time, putting tight constraints on high frequency Very Long Baseline Interferometry (VLBI) observations. In previous works, it has been shown that non-dispersive (e.g. tropospheric) effects can be calibrated with the frequency phase transfer (FPT) technique. The coherence time can thus be significantly extended. Ionospheric effects, which can still be significant, remain however uncalibrated after FPT as well as the instrumental effects. In this work, we implement a further phase transfer between two FPT residuals (i.e. so-called FPT-square) to calibrate the ionospheric effects based on their frequency dependence. We show that after FPT-square, the coherence time at 3 mm can be further extended beyond 8~hours, and the residual phase errors can be sufficiently canceled by applying the calibration of another source, which can have a large angular separation from the target (>20 deg) and significant temporal gaps. Calibrations for all-sky distributed sources with a few calibrators are also possible after FPT-square. One of the strengths and uniqueness of this calibration strategy is the suitability for high-frequency all-sky survey observations including very weak sources. We discuss the introduction of a pulse calibration system in the future to calibrate the remaining instrumental effects and allowing the possibility of imaging the source structure at high frequencies with FPT-square, where all phases are fully calibrated without involving any additional sources.

This review arose from the European Radio Astronomy Technical Forum (ERATec) meeting held in Firenze, October 2015, and aims to highlight the breadth and depth of the high-impact science that will be aided and assisted by the use of simultaneous mm-wavelength receivers. Recent results and opportunities are presented and discussed from the fields of: continuum VLBI (observations of weak sources, astrometry, observations of AGN cores in spectral index and Faraday rotation), spectral line VLBI (observations of evolved stars and massive star-forming regions) and time domain observations of the flux variations arising in the compact jets of X-ray binaries. Our survey brings together a large range of important science applications, which will greatly benefit from simultaneous observing at mm-wavelengths. Such facilities are essential to allow these applications to become more efficient, more sensitive and more scientifically robust. In some cases without simultaneous receivers the science goals are simply unachievable. Similar benefits would exist in many other high frequency astronomical fields of research.

In order to meet its cutting-edge scientific objectives, the Square Kilometre Array (SKA) telescope requires high-precision frequency references to be distributed to each of its antennas. The frequency references are distributed via fiber-optic links and must be actively stabilized to compensate for phase-noise imposed on the signals by environmental perturbations on the links. SKA engineering requirements demand that any proposed frequency reference distribution system be proved in "astronomical verification" tests. We present results of the astronomical verification of a stabilized frequency reference transfer system proposed for SKA-mid. The dual-receiver architecture of the Australia Telescope Compact Array was exploited to subtract the phase-noise of the sky signal from the data, allowing the phase-noise of observations performed using a standard frequency reference, as well as the stabilized frequency reference transfer system transmitting over 77 km of fiber-optic cable, to be directly compared. Results are presented for the fractional frequency stability and phase-drift of the stabilized frequency reference transfer system for celestial calibrator observations at 5 GHz and 25 GHz. These observations plus additional laboratory results for the transferred signal stability over a 166 km metropolitan fiber-optic link are used to show that the stabilized transfer system under test exceeds all SKA phase-stability requirements under a broad range of observing conditions. Furthermore, we have shown that alternative reference dissemination systems that use multiple synthesizers to supply reference signals to sub-sections of an array may limit the imaging capability of the telescope.

We report on the measurement of the trigonometric parallaxes of 1612 MHz hydroxyl masers around two asymptotic giant branch stars, WX Psc and OH138.0+7.2, using the NRAO Very Long Baseline Array with in-beam phase referencing calibration. We obtained a 3-sigma upper limit of <=5.3 mas on the parallax of WX Psc, corresponding to a lower limit distance estimate of >~190 pc. The obtained parallax of OH138.0+7.2 is 0.52+/-0.09 mas (+/-18%), corresponding to a distance of 1.9(+0.4,-0.3) kpc, making this the first hydroxyl maser parallax below one milliarcsecond. We also introduce a new method of error analysis for detecting systematic errors in the astrometry. Finally, we compare our trigonometric distances to published phase-lag distances toward these stars and find a good agreement between the two methods.

In this paper we describe a new approach for mm-VLBI calibration that provides bona-fide astrometric alignment of the mm-wavelength images from a single source, for the measurement of frequency dependent effects, such as `core-shifts' near the black hole of AGN jets. We achieve our astrometric alignment by solving firstly for the ionospheric (dispersive) contributions using wide-band cm-wavelength observations. Secondly we solve for the tropospheric (non-dispersive) contributions by using fast frequency-switching at the target mm-wavelengths. These solutions can be scaled and transferred from the low frequency to the high frequency. To complete the calibration chain one additional step was required to remove a residual constant phase offset on each antenna. The result is an astrometric calibration and the measurement of the core-shift between 22 and 43 GHz for the jet in BL Lacertae to be -8$\pm$5, 20$\pm$6 $\mu$as, in RA and Declination, respectively. By comparison to conventional phase referencing at cm-wavelengths we are able to show that this core shift at mm-wavelengths is significantly less than what would be predicted by extrapolating the low frequency result, which closely followed the predictions of the Blandford \& K\"onigl conical jet model. As such it would be the first demonstration for the association of the VLBI core with a recollimation shock, normally hidden at low frequencies due to the optical depth, which could be responsible for the $\gamma$-ray production in blazar jets.

The arrival of the Square Kilometer Array (SKA) will revitalise all aspects of Very Long Baseline Interferometry (VLBI) astronomy at the lower frequencies. In the last decade there have been huge strides towards routinely achieving high precision VLBI astrometry at frequencies dominated by the tropospheric contributions, most notably at 22GHz, using advanced phase referencing techniques. Nevertheless to increase the capability for high precision astrometric measurements at low radio frequencies (<8GHz) an effective calibration strategy of the systematic ionospheric propagation effects that is widely applicable is required. Observations at low frequencies are dominated by distinct direction dependent ionospheric propagation errors, which place a very tight limit on the angular separation of a suitable phase referencing calibrator. The MultiView technique holds the key to the compensation of atmospheric spatial-structure errors, by using observations of multiple calibrators and 2-D interpolation. In this paper we present the first demonstration of the power of MultiView using three calibrators, several degrees from the target, along with a comparative study of the astrometric accuracy between MultiView and phase-referencing techniques. MultiView calibration provides an order of magnitude improvement in astrometry with respect to conventional phase referencing, achieving ca. 100 micro-arcseconds astrometry errors in a single epoch of observations, effectively reaching the thermal noise limit. MultiView will achieve its full potential with the enhanced sensitivity and multibeam capabilities of SKA and the pathfinders, which will enable simultaneous observations of the target and calibrators. Our demonstration indicates that the 10 micro-arcseconds goal of astrometry at ~1.6GHz using VLBI with SKA is feasible using the MultiView technique.

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

Simultaneous observations at multiple frequency bands have the potential to overcome the fundamental limitation imposed by the atmospheric propagation in mm-VLBI observations. The propagation effects place a severe limit in the sensitivity achievable in mm-VLBI, reducing the time over which the signals can be coherently combined, and preventing the use of phase referencing and astrometric measurements. We carried out simultaneous observations at 22, 43, 87 and 130 GHz of a group of five AGNs, the weakest of which is ca. 200 mJy at 130 GHz, with angular separations ranging from 3.6 to 11 degrees, using the KVN. We analysed this data using the Frequency Phase Transfer (FPT) and the Source Frequency Phase Referencing (SFPR) techniques, which use the observations at a lower frequency to correct those at a higher frequency. The results of the analysis provide an empirical demonstration of the increase in the coherence times at 130 GHz from a few tens of seconds to about twenty minutes, with FPT, and up to many hours with SFPR. Moreover the astrometric analysis provides high precision relative position measurements between two frequencies, including, for the first time, astrometry at 130 GHz. Finally we demonstrate a method for the generalised decomposition of the relative position measurements into absolute position shifts for bona fide astrometric registration of the maps of the individual sources at multiple frequencies, up to 130 GHz.

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

We have conducted the first parallax and proper motion measurements of 6.7 GHz methanol maser emission using the Australian Long Baseline Array (LBA). The parallax of G339.884$-$1.259 measured from five epochs of observations is 0.48$\pm $0.08 mas, corresponding to a distance of $2.1^{+0.4}_{-0.3}$ kpc, placing it in the Scutum spiral arm. This is consistent (within the combined uncertainty) with the kinematic distance estimate for this source at 2.5$\pm $0.5 kpc using the latest Solar and Galactic rotation parameters. We find from the Lyman continuum photon flux that the embedded core of the young star is of spectral type B1, demonstrating that luminous 6.7 GHz methanol masers can be associated with high-mass stars towards the lower end of the mass range.

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.

Oxygen-rich Asymptotic Giant Branch (AGB) stars can be intense emitters of SiO ($v$=1 and 2, J=1$\rightarrow$0) and H$_2$O maser lines at 43 and 22 GHz, respectively. VLBI observations of the maser emission provide a unique tool to probe the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned \water\ and $v$=1 and $v$=2 SiO maser maps have traditionally limited the physical constraints that can be placed on the SiO maser pumping mechanism. We present phase referenced simultaneous spectral-line VLBI images for the SiO $v$=1 and $v$=2, J=1$\rightarrow$0, and H$_2$O maser emission around the AGB star R\,LMi, obtained from the Korean VLBI Network (KVN). The simultaneous multi-channel receivers of the KVN offer great possibilities for astrometry in the frequency domain. With this facility we have produced images with bona-fide absolute astrometric registration between high frequency maser transitions of different species to provide the positions of the \water\ maser emission, and the centre of the SiO maser emission, and hence reducing the uncertainty in the proper motion for R\,LMi by an order of magnitude over that from Hipparcos. This is the first successful demonstration of source frequency phase referencing for mm-VLBI spectral-line observations and also where the ratio between the frequencies is not an integer.

We present our work towards using the Korean VLBI (Very Long Baseline Interferometer) Network (KVN) and VLBI Exploration of Radio Astronomy (VERA) arrays combined into the KVN and VERA Array (KaVA) for observations of radio pulsars at high frequencies ($\simeq$22-GHz). Pulsar astronomy is generally focused at frequencies approximately 0.3 to several GHz and pulsars are usually discovered and monitored with large, single-dish, radio telescopes. For most pulsars, reduced radio flux is expected at high frequencies due to their steep spectrum, but there are exceptions where high frequency observations can be useful. Moreover, some pulsars are observable at high frequencies only, such as those close to the Galactic Center. The discoveries of a radio-bright magnetar and a few dozen extended Chandra sources within 15 arc-minute of the Galactic Center provide strong motivations to make use of the KaVA frequency band for searching pulsars in this region. Here, we describe the science targets and report progresses made from the KVN test observations for known pulsars. We then discuss why KaVA pulsar observations are compelling.

The Korean VLBI Network (KVN) is a new mm-VLBI dedicated array with capability for simultaneous observations at multiple frequencies, up to 129 GHz. The innovative multi-channel receivers present significant benefits for astrometric measurements in the frequency domain. The aim of this work is to verify the astrometric performance of the KVN using a comparative study with the VLBA, a well established instrument. For that purpose, we carried out nearly contemporaneous observations with the KVN and the VLBA, at 14/7 mm, in April 2013. The KVN observations consisted of simultaneous dual frequency observations, while the VLBA used fast frequency switching observations. We used the Source Frequency Phase Referencing technique for the observational and analysis strategy. We find that having simultaneous observations results in a superior performance for compensation of all atmospheric terms in the observables, in addition to offering other significant benefits for astrometric analysis. We have compared the KVN astrometry measurements to those from the VLBA. We find that the structure blending effects introduce dominant systematic astrometric shifts and these need to be taken into account. We have tested multiple analytical routes to characterize the impact of the low resolution effects for extended sources in the astrometric measurements. The results from the analysis of KVN and full VLBA datasets agree within 2-$\sigma$ of the thermal error estimate. We interpret the discrepancy as arising from the different resolutions. We find that the KVN provides astrometric results with excellent agreement, within 1-$\sigma$, when compared to a VLBA configuration which has a similar resolution. Therefore this comparative study verifies the astrometric performance of KVN using SFPR at 14/7 mm, and validates the KVN as an astrometric instrument.

(Abridged) Aims. To investigate the first stages of the process of high-mass star formation, we selected a sample of massive clumps previously observed with the SEST at 1.2 mm and with the ATNF ATCA at 1.3 cm. We want to characterize the physical conditions in such sources, and test whether their properties depend on the evolutionary stage of the clump. Methods. With ATCA we observed the selected sources in the NH3(1,1) and (2,2) transitions and in the 22 GHz H2O maser line. Ammonia lines are a good temperature probe that allow us to accurately determine the mass and the column-, volume-, and surface densities of the clumps. We also collected all data available to construct the spectral energy distribution of the individual clumps and to determine if star formation is already occurring, through observations of its most common signposts, thus putting constraints on the evolutionary stage of the source. We fitted the spectral energy distribution between 1.2 mm and 70 microns with a modified black body to derive the dust temperature and independently determine the mass. Results. The clumps are cold (T~10-30 K), massive (M~10^2-10^3 Mo), and dense (n(H2)>~10^5 cm^-3) and they have high column densities (N(H2)~10^23 cm^-2). All clumps appear to be potentially able to form high-mass stars. The most massive clumps appear to be gravitationally unstable, if the only sources of support against collapse are turbulence and thermal pressure, which possibly indicates that the magnetic field is important in stabilizing them. Conclusions. After investigating how the average properties depend on the evolutionary phase of the source, we find that the temperature and central density progressively increase with time. Sources likely hosting a ZAMS star show a steeper radial dependence of the volume density and tend to be more compact than starless clumps.

We present new wideband radio observations with the Australia Telescope Compact Array of a sample of 12 giant radio galaxies. The radio observations are part of a larger radio-optical study aimed at relating the radio structures with the ambient medium on large scales. With projected linear sizes larger than 0.7 Mpc, these objects are ideal candidates for the study of the Warm-Hot Intergalactic Medium (WHIM). The sample includes sources with sizes spanning 0.8 to 3.2 Mpc and total powers of 1.2*10^24 to 4.0*10^26 W Hz^-1 at 2.1 GHz. Redshifts were limited to z<0.15 to permit spectroscopic observations of the hosts and neighbouring galaxies, which were obtained using the AAOmega spectrograph on the Anglo-Australian Telescope. We derive lobe energy densities from the radio observations via equipartition arguments. The inferred pressures in the lobes of the giant radio sources, which range from 1.1*10^-15 to 2.0*10^-14 Pa (80 to 1500 cm^-3 K), are lower than previously inferred from X-ray observations of dense filaments. Comparison with the OverWhelmingly Large Simulations (OWLS) suggests that the WHIM in pressure balance with the radio lobes has a temperature in excess of ~10^6.5 K or a particle overdensity in the range 50 to 500. This study highlights the capability of next generation surveys, such as the Evolutionary Map of the Universe (EMU) survey with the Australian Square Kilometre Array Pathfinder (ASKAP), to study populations of giant radio sources at lower surface brightness and thereby discriminate between models for the cosmological evolution of the intergalactic medium and examine the validity of cosmological hydrodynamical simulations.

High precision astrometric Space Very Long Baseline Interferometry (S-VLBI) at the low end of the conventional frequency range, i.e. 20cm, is a requirement for a number of high priority science goals. These are headlined by obtaining trigonometric parallax distances to pulsars in Pulsar--Black Hole pairs and OH masers anywhere in the Milky Way Galaxy and the Magellanic Clouds. We propose a solution for the most difficult technical problems in S-VLBI by the MultiView approach where multiple sources, separated by several degrees on the sky, are observed simultaneously. We simulated a number of challenging S-VLBI configurations, with orbit errors up to 8m in size and with ionospheric atmospheres consistant with poor conditions. In these simulations we performed MultiView analysis to achieve the required science goals. This approach removes the need for beam switching requiring a Control Moment Gyro, and the space and ground infrastructure required for high quality orbit reconstruction of a space-based radio telescope. This will dramatically reduce the complexity of S-VLBI missions which implement the phase-referencing technique.

We have carried out full imaging simulation studies to explore the impact of frequency standards in millimeter and sub-millimeter Very Long Baseline Interferometry (VLBI), focusing on the coherence time and sensitivity. In particular, we compare the performance of the H-maser, traditionally used in VLBI, to that of ultra-stable cryocooled sapphire oscillators over a range of observing frequencies, weather conditions and analysis strategies. Our simulations show that at the highest frequencies, the losses induced by H-maser instabilities are comparable to those from high quality tropospheric conditions. We find significant benefits in replacing H-masers with cryocooled sapphire oscillator based frequency references in VLBI observations at frequencies above 175 GHz in sites which have the best weather conditions; at 350 GHz we estimate a 20-40% increase in sensitivity, over that obtained when the sites have H-masers, for coherence losses of 20-10%, respectively. Maximum benefits are to be expected by using colocated Water Vapour Radiometers for atmospheric correction. In this case, we estimate a 60-120% increase in sensitivity over the H-maser at 350 GHz.

Infrared-dark high-mass clumps are among the most promising objects to study the initial conditions of the formation process of high-mass stars and rich stellar clusters. In this work, we have observed the (3-2) rotational transition of C18O with the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 1.2 mm continuum, with the aim of measuring basic chemical and physical parameters such as the CO depletion factor (fD), the gas kinetic temperature and the gas mass. In particular, the C18O (3-2) line allows us to derive fD in gas at densities higher than that traced by the (1-0) and (2-1) lines, typically used in previous works. We have detected NH3 and C18O in all targets. The clumps possess mass, H2 column and surface densities consistent with being potentially the birthplace of high-mass stars. We have measured fD in between 5 and 78, with a mean value of 32 and a median of 29. These values are, to our knowledge, larger than the typical CO depletion factors measured towards infrared-dark clouds and high-mass dense cores, and are comparable to or larger than the values measured in low-mass pre-stellar cores close to the onset of the gravitational collapse. This result suggests that the earliest phases of the high-mass star and stellar cluster formation process are characterised by fD larger than in low-mass pre-stellar cores. Thirteen out of 21 clumps are undetected in the 24 {\mu}m Spitzer images, and have slightly lower kinetic temperatures, masses and H2 column densities with respect to the eight Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less evolved or are going to form less massive objects.

We present the results from an ultra-high-resolution 7mm Very Long Baseline Array (VLBA) study of the relativistic jet in the BL Lacertae object OJ287 from 1995 to 2011 containing 136 total intensity images. Analysis of the image sequence reveals a sharp jet-position-angle swing by >100 deg. during [2004,2006], as viewed in the plane of the sky, that we interpret as the crossing of the jet from one side of the line of sight to the other during a softer and longer term swing of the inner jet. Modulating such long term swing, our images also show for the first time a prominent erratic wobbling behavior of the innermost ~0.4mas of the jet with fluctuations in position angle of up to ~40 deg. over time scales ~2yr. This is accompanied by highly superluminal motions along non-radial trajectories, which reflect the remarkable non-ballistic nature of the jet plasma on these scales. The erratic nature and short time scales of the observed behavior rules out scenarios such as binary black hole systems, accretion disk precession, and interaction with the ambient medium as possible origins of the phenomenon on the scales probed by our observations, although such processes may cause longer-term modulation of the jet direction. We propose that variable asymmetric injection of the jet flow; perhaps related to turbulence in the accretion disk; coupled with hydrodynamic instabilities, leads to the non-ballistic dynamics that cause the observed non-periodic changes in the direction of the inner jet.

Space Very-Long-Baseline-Interferometry (S-VLBI) observations at high frequencies hold the prospect of achieving the highest angular resolutions and astrometric accuracies, resulting from the long baselines between ground and satellite telescopes. Nevertheless, space-specific issues, such as limited accuracy in the satellite orbit reconstruction and constraints on the satellite antenna pointing operations, limit the application of conventional phase referencing. We investigate the feasibility of an alternative technique, source frequency phase referencing (SFPR), to the S-VLBI domain. With these investigations we aim to contribute to the design of the next-generation of S-VLBI missions. We have used both analytical and simulation studies to characterize the performance of SFPR in S-VLBI observations, applied to astrometry and increased coherence time, and compared these to results obtained using conventional phase referencing. The observing configurations use the specifications of the ASTRO-G mission for their starting point. Our results show that the SFPR technique enables astrometry at 43 GHz, using alternating observations with 22 GHz, regardless of the orbit errors, for most weathers and under a wide variety of conditions. The same applies to the increased coherence time for the detection of weak sources. Our studies show that the capability to carry out simultaneous dual frequency observations enables the application to higher frequencies, and a general improvement of the performance in all cases, hence we recommend its consideration for S-VLBI programs.

We describe a new method which achieves high precision Very Long Baseline Interferometry (VLBI) astrometry in observations at millimeter wavelengths. It combines fast frequency-switching observations, to correct for the dominant non-dispersive tropospheric fluctuations, with slow source-switching observations, for the remaining ionospheric dispersive terms. We call this method Source-Frequency Phase Referencing. Provided that the switching cycles match the properties of the propagation media, one can recover the source astrometry. We present an analytic description of the two-step calibration strategy, along with an error analysis to characterize its performance. Also, we provide observational demonstrations of a successful application with observations using the Very Long Baseline Array at 86 GHz of the pairs of sources 3C274 & 3C273 and 1308+326 & 1308+328, under various conditions. We conclude that this method is widely applicable to millimeter VLBI observations of many target sources, and unique in providing bona-fide astrometrically registered images and high precision relative astrometric measurements in mm-VLBI using existing and newly built instruments.

IRAS 20126+4104 is one of the best candidates for a high-mass (proto)star surrounded by an accretion disk. Such a fact may be used to set constraints on theories of high-mass star formation, but requires confirmation that the mass and luminosity of IRAS 20126+4104 are indeed typical of a B0.5 star, which in turn requires an accurate estimate of the distance. We used the Very Long Baseline Array and the European VLBI Network to observe the 22.2 GHz water and 6.7 GHz methanol masers in IRAS 20126+4104 at a number of epochs suitably distributed in time. The absolute positions of the maser features were established with respect to reference quasars, which allowed us to derive absolute proper motions. From the parallax of the water masers we obtain a distance of 1.64 \pm 0.05 kpc, which is very similar to the value adopted so far in the literature (1.7 kpc) and confirms that IRAS 20126+4104 is a high-mass (proto)star. From the methanol masers we derive the component in the plane of the sky of the systemic velocity of the disk+star system (-16 km/s in right-ascension and +7.6 km/s in declination). Accurate knowledge of the distance and systemic velocity allows us to improve on the model fit to the water maser jet presented in a previous study. Finally, we identify two groups of methanol maser features, one undergoing rotation in the disk and possibly distributed along a narrow ring centered on the star, the other characterised by relative proper motions indicating that the features are moving away from the disk, perpendicular to it. We speculate that the latter group might be tracing the disk material marginally entrained by the jet.

This document describes the advantages of applying "Source/Frequency Phase Referencing" (SFPR) techniques to the analysis of VLBI observations with VSOP-2, for high precision astrometric measurements and/or increased sensitivity. The SFPR calibration technique basics and a demonstration of the method applied to highest frequency VLBA observations are described in detail in VLBA Scientific Memo n. 31. Here we outline its importance in the context of space VLBI astrometry with VSOP-2, where errors in the satellite orbit determination and rapid tropospheric phase fluctuations set extreme challenges for the successful application of conventional phase referencing techniques, specially at the higher frequencies. SFPR is ideally suited for full calibration of those - regardless of the orbit determination accuracy - and, in general, of any non-dispersive terms. The requirements for application of SFPR techniques are fully compatible with current technical specifications of VSOP-2. Hence we foresee that SFPR will play an important role in helping expanding the scientific outcome of the space VLBI mission.

In this document we layout a new method to achieve "bona fide" high precision Very-Long-Baseline-Interferometry (VLBI) astrometric measurements of frequency-dependent positions of celestial sources (even) in the high (mm-wavelength) frequency range, where conventional phase referencing techniques fail. Our method, dubbed "Source/Frequency Phase Referencing" (SFPR) combines fast frequency-switching (or dual-frequency observations) with the source switching of conventional phase referencing techniques. The former is used to calibrate the dominant highly unpredictable rapid atmospheric fluctuations, which arise from variations of the water vapor content in the troposphere, and ultimately limit the application of conventional phase referencing techniques; the latter compensates the slower time scale remaining ionospheric/instrumental, non-negligible, phase variations. For cm-VLBI, the SFPR method is equivalent to conventional phase referencing applied to the measurement of frequency-dependent source positions changes ("core-shifts"). For mm-VLBI, the SFPR method stands as the only approach which will provide astrometry. In this memo we layout the scope and basis of our new method, along with a description of the strategy and a successful demonstration of the application of this new astrometric analysis technique to the highest frequency VLBA observations, at 86 GHz. Our previous comparative astrometric analysis of cm-VLBI observations, presented elsewhere, produced equivalent results using both methods.

We revisit the "Cluster-Cluster" or multi-view Very-Long-Baseline-Interferometry (VLBI) technique from the perspective of its synergy with the multi-beam features inherent in the Australian Square Kilometer Array Pathfinder (ASKAP)and its potential to improve the outcomes of VLBI observations with ASKAP. We include a list of candidate VLBI sites that already support or can be upgraded to support multi-view VLBI located in Australia and overseas, and which have common visibility with ASKAP. The results of our previous "cluster-cluster" observations at 1.6 GHz demonstrated the advantages of this configuration to calibrate the ionospheric distortions responsible for the loss of positional accuracy at low frequencies, using multiple calibrators in a range between 1 to 6 degrees away from the target. Therefore, we conclude that joint observations of ASKAP with other multi-view sites using "cluster-cluster" techniques would improve the outcomes of the high spatial resolution component of ASKAP applied to astrometric projects, achieving higher precision for many more targets, and with lower detection thresholds. Also very wide-field VLBI mapping becomes a possibility. Looking to the future, this would contribute to the development of new techniques that are relevant for future high resolution observations with the SKA.

Oxygen-rich Asymptotic Giant Branch (AGB) stars are intense emitters of SiO and H$_2$O maser lines at 43 (J=1-0, v=1 and 2) and 22 GHz, respectively. VLBI observations of the maser emission provides a unique tool to sample the innermost layers of the circumstellar envelopes in AGB stars. Nevertheless, the difficulties in achieving astrometrically aligned v=1 and v=2 SiO maser maps have traditionally prevented a unique interpretation of the observations in terms of physical underlying conditions, which depend on the nature of the SiO pumping mechanism. We have carried out observations of the SiO and H$_2$O maser emission towards RLMi, using the astrometric capabilities of VERA. Due to the too-weak emission of the reference calibrator we had to develop a special method to accurately relate the coordinates for both transitions. We present relative astrometrically aligned v=1 and v=2 J=1-0 SiO maser maps, at multiple epochs, and discuss the astrophysical results. The incorporation of astrometric information into the maps of SiO masers challenges the weak points in the current theoretical models, which will need further refinements to address the observations results.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. Approximately twenty-five percent of the mission time was dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz. This paper, the fifth in the series, presents images and models for the remaining 140 sources not included in Paper III, which contained 102 sources. For most sources, the plots of the uv-coverage, the visibility amplitude versus uv-distance, and the high resolution image are presented. Model fit parameters to the major radio components are determined, and the brightness temperature of the core component for each source is calculated. The brightness temperature distributions for all of the sources in the VSOP AGN survey are discussed.

VLBI is unique, among the space geodetic techniques, in its contribution to defining and maintaining the International Celestial Reference Frame, providing precise measurements of coordinates of extragalactic radiosources. The quest for increasing accuracy of VLBI geodetic products has lead to a deeper revision of all aspects that might introduce errors in the analysis. The departure of the observed sources from perfect, stable, compact and achromatic celestial targets falls within this category. This paper is concerned with the impact of unaccounted frequency-dependent position shifts of source cores in the analysis of dual-band S/X VLBI geodesy observations, and proposes a new method to measure them. The multi-frequency phase transfer technique developed and demonstrated by Middelberg et al. (2005) increases the high frequency coherence times of VLBI observations, using the observations at a lower frequency. Our proposed SOURCE/FREQUENCY PHASE REFERENCING method endows it with astrometric applications by adding a strategy to estimate the ionospheric contributions. Here we report on the first successful application to measure the core shift of the quasar 1038+528 A at S and X-bands, and validate the results by comparison with those from standard phase referencing techniques. In this particular case, and in general in the cm-wavelength regime, both methods are equivalent. Moreover the proposed method opens a new horizon with targets and fields suitable for high precision astrometric studies with VLBI, especially at high frequencies where severe limitations imposed by the rapid fluctuations in the troposphere prevent the use of standard phase referencing techniques.

We report on the multi-epoch observations of H2O maser emission in star forming region OH 43.8-0.1 carried out with VLBI Exploration of Radio Astrometry (VERA). The large-scale maser distributions obtained by single-beam VLBI mapping reveal new maser spots scattered in area of 0.7 x 1.0 arcsec, in addition to a `shell-like' structure with a scale of 0.3 x 0.5 arcsec which was previously mapped by Downes et al.(1979). Proper motions are also obtained for 43 spots based on 5-epoch monitoring with a time span of 281 days. The distributions of proper motions show a systematic outflow in the north-south direction with an expansion velocity of ~8 km/s, and overall distributions of maser spots as well as proper motions are better represented by a bipolar flow plus a central maser cluster with a complex structure, rather than a shell with a uniform expansion such as those found in Cep A R5 and W75N VLA2. The distance to OH 43.8-0.1 is also estimated based on the statistical parallax, yielding D = 2.8 +/- 0.5 kpc. This distance is consistent with a near kinematic distance and rules out a far kinematics distance (~9 kpc), and the LSR velocity of OH 43.8-0.1 combined with the distance provides a constraint on the flatness of the galactic rotation curve, that there is no systematic difference in rotation speeds at the Sun and at the position of OH 43.8-0.1, which is located at the galacto-centric radius of ~6.3 kpc.

We report on the first multi-epoch, phase referenced VLBI observations of the water maser emission in a high-mass protostar associated with a disk-jet system. The source under study, IRAS 20126+4104, has been extensively investigated in a large variety of tracers, including water maser VLBA data acquired by us three years before the present observations. The new findings fully confirm the interpretation proposed in our previous study, namely that the maser spots are expanding from a common origin coincident with the protostar. We also demonstrate that the observed 3-D velocities of the maser spots can be fitted with a model assuming that the spots are moving along the surface of a conical jet, with speed increasing for increasing distance from the cone vertex. We also present the results of single-dish monitoring of the water maser spectra in IRAS 20126+4104. These reveal that the peak velocity of some maser lines decreases linearly with time. We speculate that such a deceleration could be due to braking of the shocks from which the maser emission originates, due to mass loading at the shock front or dissipation of the shock energy.

We present VLBI observations of the H$_2$O maser emission towards a selection of massive young stellar objects (YSOs). In one of these sources, IRAS20126+4104, the measured proper motions confirm that the water masers spots are tracing the very root (around 100 AU) of a conical bipolar outflow emerging near the position of the embedded YSO, and are an excellent tool to investigate the structure and kinematics of the outflow/jet. We also present preliminary results of a first epoch EVN observations of a selection of 5 high-mass YSO to assess the precise association of maser spots and molecular outflows in these sources.

We report on a strong outburst of the W49N H2O maser observed with VERA. Single-dish monitoring with VERA 20 m telescopes detected a strong outburst of the maser feature at V_LSR = -30.7 km/s in 2003 October. The outburst had a duration of ~100 days and a peak intensity of 7.9 x 10^4 Jy, being one of the strongest outbursts in W49N observed so far. VLBI observations with the VERA array were also carried out near to the maximum phase of the outburst, and the outburst spot was identified in the VLBI map. While the map was in good agreement with previous studies, showing three major concentrations of maser spots, we found a newly formed arc-like structure in the central maser concentration, which may be a shock front powered by a forming star or a star cluster. The outburst spot was found to be located on the arc-like structure, indicating a possible connection of the present outburst to a shock phenomenon.

We describe, and present preliminary results from the analysis of 1.6 GHz VLBI observations, made in Cluster-Cluster mode between subgroups of antennas from the WSRT, VLA and MERLIN arrays. We observed for two 6-hour periods on consecutive days, simultaneously monitoring 4 compact radio sources with angular separations ranging from ~1 to 9 degrees. The data were recorded in a non-standard MkIII mode using special set-ups, and correlated at the MPIfR MkIIIa correlator. The preliminary results from standard and phase reference analyses are very encouraging. The experiment has provided a useful database for investigating the potential of this new observing technique at low frequencies.

We present results from 3.6 cm observations of the quasar pair 1038+528 A and B, made in 1995 using the VLBA together with the Effelsberg 100m telescope. We describe the use of a phase-referencing technique to measure the astrometric separation between the quasars. We also introduce a new data analysis method -"hybrid double mapping"- which preserves the relative astrometric information in a single VLBI hybrid map for close source pairs. We combine our measurements with those from three previous epochs, the earliest in 1981. Our new observations confirm the evolution within the structure of quasar B, previously proposed to explain the measured change in the relative separation of the pair. Our upper bound for any systematic proper motion between the mass centres of quasars A and B is 10 micro-arcsec/year. This is set by the limited precision in defining the reference points in the quasars at different epochs and by possible instabilities of the source "core" locations. A separate analysis enables us to put more stringent upper limits to any core motions along the two source axes.

We present here new VLBI observations of one FR-I radio galaxy (NGC2484) and two Broad Line FR-II radio galaxies (3C109 and 3C382). For 3C109 new VLA maps are also shown. These sources belong to a complete sample of radio galaxies under study for a better knowledge of their structures at parsec resolution. The parsec structure of these 3 objects is very similar: asymmetric emission which we interpret as the core plus a one-sided jet. The parsec scale jet is always on the same side of the main kpc-scale jet. The limit on the jet to counterjet brightness ratio, the ratio of the core radio power to the total radio power and the Synchrotron-self Compton model allow us to derive some constraints on the jet velocity and orientation with respect to the line of sight. From these data and from those published on 2 other sources of our sample, we suggest that parsec scale jets are relativistic in both FR-I and FR-II radio galaxies and that pc scale properties in FR-I and FR-II radio galaxies are very similar despite of the large difference between these two classes of radio galaxies on the kpc scale.