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.

Studies of fundamental parameters of very low-mass objects are indispensable to provide tests of stellar evolution models that are used to derive theoretical masses of brown dwarfs and planets. However, only objects with dynamically determined masses and precise photometry can effectively evaluate the predictions of stellar models. AB Dor C (0.090 solar masses) has become a prime benchmark for calibration of theoretical evolutionary models of low-mass young stars. One of the ambiguities remaining in AB Dor C is the possible binary nature of this star. We observed AB Dor C with the VLTI/AMBER instrument in low-resolution mode at the J, H and K bands. The interferometric observables at the K-band are compatible with a binary brown dwarf system with tentative components AB Dor Ca/Cb with a K-band flux ratio of 5$\pm$1% and a separation of 38$\pm$1 mas. This implies theoretical masses of 0.072$\pm$0.013 M$_{\rm \odot}$ and 0.013$\pm$0.001 M$_{\rm \odot}$ for each component, near the hydrogen-burning limit for AB Dor Ca, and near the deuterium-burning limit, straddling the boundary between brown dwarfs and giant planets, for AB Dor Cb. The possible binarity of AB Dor C alleviates the disagreement between observed magnitudes and theoretical mass-luminosity relationships.

(Abridged) The CASA software suite, can now reduce very long baseline interferometry (VLBI) data with the recent addition of a fringe fitter. Here, we present the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is an open-source VLBI calibration and imaging pipeline built on top of the CASA framework. The pipeline is capable of reducing data from different VLBI arrays. It can be run non-interactively after only a few non-default input parameters are set and delivers high-quality calibrated data. CPU scalability based on a message-passing interface (MPI) implementation ensures that large bandwidth data from future arrays can be processed within reasonable computing times. Phase calibration is done with a Schwab-Cotton fringe fit algorithm. For the calibration of residual atmospheric effects, optimal solution intervals are determined based on the signal-to-noise ratio (S/N) of the data for each scan. Different solution intervals can be set for different antennas in the same scan to increase the number of detections in the low S/N regime. These novel techniques allow rPICARD to calibrate data from different arrays, including high-frequency and low-sensitivity arrays. The amplitude calibration is based on standard telescope metadata, and a robust algorithm can solve for atmospheric opacity attenuation in the high-frequency regime. Standard CASA tasks are used for CLEAN imaging and self-calibration. In this work we demonstrate the capabilities of rPICARD by calibrating and imaging 7 mm VLBA data of the central radio source in the M87 galaxy. The reconstructed jet image reveals a complex collimation profile and edge-brightened structure. A potential counter-jet is detected that has 10 % of the brightness of the approaching jet. This constrains jet speeds close to the radio core to about half the speed of light for small inclination angles.

Currently, HOPS and AIPS are the primary choices for the time-consuming process of (millimeter) Very Long Baseline Interferometry (VLBI) data calibration. However, for a full end-to-end pipeline, they either lack the ability to perform easily scriptable incremental calibration or do not provide full control over the workflow with the ability to manipulate and edit calibration solutions directly. The Common Astronomy Software Application (CASA) offers all these abilities, together with a secure development future and an intuitive Python interface, which is very attractive for young radio astronomers. Inspired by the recent addition of a global fringe-fitter, the capability to convert FITS-IDI files to measurement sets, and amplitude calibration routines based on ANTAB metadata, we have developed the the CASA-based Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD). The pipeline will be able to handle data from multiple arrays: EHT, GMVA, VLBA and the EVN in the first release. Polarization and phase-referencing calibration are supported and a spectral line mode will be added in the future. The large bandwidths of future radio observatories ask for a scalable reduction software. Within CASA, a message passing interface (MPI) implementation is used for parallelization, reducing the total time needed for processing. The most significant gain is obtained for the time-consuming fringe-fitting task where each scan be processed in parallel.

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.

We studied the polarization behavior of the quasar 3C 273 over the 1 mm wavelength band at ALMA with a total bandwidth of 7.5GHz across 223 to 243 GHz at 0.8 arcsec resolution, corresponding to 2.1 kpc at the distance of 3C 273. With these observations we were able to probe the optically thin polarized emission close to the jet base, and constrain the magnetic field structure. We computed the Faraday rotation measure using simple linear fitting and Faraday rotation measure synthesis. In addition, we modeled the broadband behavior of the fractional Stokes Q and U parameters (qu-fitting). The systematic uncertainties in the polarization observations at ALMA were assessed through Monte Carlo simulations. We find the unresolved core of 3C 273 to be 1.8% linearly polarized. We detect a very high rotation measure (RM) of ~5.0 x10^5 rad/m^2 over the 1 mm band when assuming a single polarized component and an external RM screen. This results in a rotation of >40 deg of the intrinsic electric vector position angle, which is significantly higher than typically assumed for millimeter wavelengths. The polarization fraction increases as a function of wavelength, which according to our qu-fitting could be due to multiple polarized components of different Faraday depth within our beam, or internal Faraday rotation. With our limited wavelength coverage, we cannot distinguish between the cases, and additional multifrequency and high angular resolution observations are needed to determine the location and structure of the magnetic field of the Faraday active region. Comparing our RM estimate with values obtained at lower frequencies, the RM increases as a function of observing frequency, following a power law with an index of ~2.0 consistent with a sheath surrounding a conically expanding jet. We also detect ~0.2% circular polarization, although further observations are needed to confirm this result.

The Galactic Center supermassive black hole Sagittarius A* (Sgr A*) is one of the most promising targets to study the dynamics of black hole accretion and outflow via direct imaging with very long baseline interferometry (VLBI). At 3.5 mm (86 GHz), the emission from Sgr A* is resolvable with the Global Millimeter VLBI Array (GMVA). We present the first observations of Sgr A* with the phased Atacama Large Millimeter/submillimeter Array (ALMA) joining the GMVA. Our observations achieve an angular resolution of ~87{\mu}as, improving upon previous experiments by a factor of two. We reconstruct a first image of the unscattered source structure of Sgr A* at 3.5 mm, mitigating effects of interstellar scattering. The unscattered source has a major axis size of 120 $\pm$ 34{\mu}as (12 $\pm$ 3.4 Schwarzschild radii), and a symmetrical morphology (axial ratio of 1.2$^{+0.3}_{-0.2}$), which is further supported by closure phases consistent with zero within 3{\sigma}. We show that multiple disk-dominated models of Sgr A* match our observational constraints, while the two jet-dominated models considered are constrained to small viewing angles. Our long-baseline detections to ALMA also provide new constraints on the scattering of Sgr A*, and we show that refractive scattering effects are likely to be weak for images of Sgr A* at 1.3 mm with the Event Horizon Telescope. Our results provide the most stringent constraints to date for the intrinsic morphology and refractive scattering of Sgr A*, demonstrating the exceptional contribution of ALMA to millimeter VLBI.

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

We report results from very long baseline interferometric (VLBI) observations of the supermassive black hole in the Galactic center, Sgr A*, at 1.3 mm (230 GHz). The observations were performed in 2013 March using six VLBI stations in Hawaii, California, Arizona, and Chile. Compared to earlier observations, the addition of the APEX telescope in Chile almost doubles the longest baseline length in the array, provides additional {\it uv} coverage in the N-S direction, and leads to a spatial resolution of $\sim$30 $\mu$as ($\sim$3 Schwarzschild radii) for Sgr A*. The source is detected even at the longest baselines with visibility amplitudes of $\sim$4-13% of the total flux density. We argue that such flux densities cannot result from interstellar refractive scattering alone, but indicate the presence of compact intrinsic source structure on scales of $\sim$3 Schwarzschild radii. The measured nonzero closure phases rule out point-symmetric emission. We discuss our results in the context of simple geometric models that capture the basic characteristics and brightness distributions of disk- and jet-dominated models and show that both can reproduce the observed data. Common to these models are the brightness asymmetry, the orientation, and characteristic sizes, which are comparable to the expected size of the black hole shadow. Future 1.3 mm VLBI observations with an expanded array and better sensitivity will allow a more detailed imaging of the horizon-scale structure and bear the potential for a deep insight into the physical processes at the black hole boundary.

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

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

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.

Context. Precise determination of stellar masses is necessary to test the validity of pre-main-sequence (PMS) stellar evolutionary models, whose predictions are in disagreement with measurements for masses below 1.2 Msun. To improve such a test, and based on our previous studies, we selected the AB Doradus moving group (AB Dor-MG) as the best-suited association on which to apply radio-based high-precision astrometric techniques to study binary systems. Aims. We seek to determine precise estimates of the masses of a set of stars belonging to the AB Dor-MG using radio and infrared observations. Methods. We observed in phase-reference mode with the Very Large Array (VLA) at 5 GHz and with the European VLBI Network (EVN) at 8.4 GHz the stars HD 160934, EK Dra, PW And, and LO Peg. We also observed some of these stars with the near-infrared CCD AstraLux camera at the Calar Alto observatory to complement the radio observations. Results. We determine model-independent dynamical masses of both components of the star HD 160934, A and c, which are 0.70+/-0.07 Msun and 0.45+/-0.04 Msun , respectively. We revised the orbital parameters of EK Dra and we determine a sum of the masses of the system of 1.38+/-0.08 Msun. We also explored the binarity of the stars LO Peg and PW And. Conclusions. We found observational evidence that PMS evolutionary models underpredict the mass of PMS stars by 10%-40%, as previously reported by other authors. We also inferred that the origin of the radio emission must be similar in all observed stars, that is, extreme magnetic activity of the stellar corona that triggers gyrosynchrotron emission from non-thermal, accelerated electrons.

The new generation of low-frequency radio telescopes, such as the Low Frequency Array (LOFAR: a Square Kilometre Array-low pathfinder), provides advancements in our capability of probing Galactic magnetism through low-frequency polarimetry. Maps of diffuse polarized radio emission and Faraday rotation can be used to infer properties of, and trace structure in, the magnetic fields in the ISM. However, to date very little of the sky has been probed at high angular and Faraday depth resolution. We observed a 5x5 degree region centred on the nearby galaxy IC342 using LOFAR in the frequency range 115-178 MHz at 4 arcmin resolution and performed Faraday tomography to detect foreground Galactic polarized synchrotron emission separated by Faraday depth (different amounts of Faraday rotation). Our Faraday depth cube shows rich polarized structure, with up to 30 K of polarized emission at 150 MHz. We detect two overlapping diffuse polarized features that are clearly separated in Faraday depth. Faraday-thick structures at such low frequencies would be too strongly depolarized to explain the observations and are therefore rejected. Only Faraday thin structures will not be strongly depolarized; producing such structures requires localized variations in the ratio of synchrotron emissivity to Faraday depth per unit distance, which can arise from several physical phenomena, such as a transition between regions of ionized and neutral gas. We conclude that the observed polarized emission is Faraday thin, and propose that the emission originates from two neutral clouds in the local ISM. We have modeled the Faraday rotation for this line of sight and estimated that the line of sight component of magnetic field of the local ISM for this direction varies between -0.86 and +0.12 uG. We propose that this may be a useful method for mapping magnetic fields within the local ISM.

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 report on multi-frequency Very Large Array (VLA) and Very Long Baseline Interferometry (VLBI) radio observations for a monitoring campaign of supernova SN$\,$2008iz in the nearby irregular galaxy M$\,$82. We fit two models to the data, a simple time power-law, S$\propto$ t$^\beta$, and a simplified Weiler model, yielding decline indices, $\beta$ of -1.22$\pm$0.07 (days 100-1500) and -1.41$\pm$0.02 (days 76-2167), respectively. The late-time radio light curve evolution shows flux-density flares at $\sim\,$970 and $\sim\,$1400 days which are a factor of $\sim\,$2 and $\sim\,$4 higher than the expected flux, respectively. The later flare, besides being brighter, does not show signs of decline at least from results examined so far (2014 January 23, day 2167). We derive the spectral index, $\alpha$, S$\propto$ $\nu^\alpha$ for frequencies 1.4 to 43$\,$GHz for SN$\,$2008iz during the period from $\sim\,$430 to 2167 days after the supernova explosion. The value of $\alpha$ shows no signs of evolution and remains steep $\approx$$-$1 throughout the period, unlike that of SN$\,$1993J which started flattening at $\sim\,$day 970. From the 4.8 and 8.4$\,$GHz VLBI images, the supernova expansion is seen to start with shell like structure that gets more and more asymmetric, then breaks up in the later epochs with bright structures dominating the southern part of the ring. This structural evolution differs significantly from SN$\,$1993J, which remains circularly symmetric over 4000 days after the explosion. The VLBI 4.8 and 8.4$\,$GHz images are used to derive a deceleration index, $m$, for SN$\,$2008iz, of 0.86$\,\pm\,$0.02, and the average expansion velocity between days 73 and 1400 as (12.1$\pm$0.2)$\times$10$^{3}\,$km$\,$s$^{-1}$.

Supermassive black holes (SMBH) are essential for the production of jets in radio-loud active galactic nuclei (AGN). Theoretical models based on Blandford & Znajek extract the rotational energy from a Kerr black hole, which could be the case for NGC1052, to launch these jets. This requires magnetic fields of the order of $10^3\,$G to $10^4\,$G. We imaged the vicinity of the SMBH of the AGN NGC1052 with the Global Millimetre VLBI Array and found a bright and compact central feature, smaller than 1.9 light days (100 Schwarzschild radii) in radius. Interpreting this as a blend of the unresolved jet bases, we derive the magnetic field at 1 Schwarzschild radius to lie between 200 G and ~80000 G consistent with Blandford & Znajek models.

Our aim is to characterize the polarized continuum emission properties including intensity, polarization position angle, and polarization percentage of Sgr A* at $\sim$100 (3.0 mm), $\sim$230 (1.3 mm), $\sim$345 (0.87 mm), $\sim$500 (0.6 mm), and $\sim$700 GHz (0.43 mm). We report continuum emission properties of Sgr A* at the above frequency bands, based on the Atacama Large Millimeter Array (ALMA) observations. We measured flux densities of Sgr A* from ALMA single pointing and mosaic observations. We performed sinusoidal fittings to the observed (XX-YY)/I intensity ratios, to derive the polarization position angles and polarization percentages. We successfully detect polarized continuum emission from all observed frequency bands. We observed lower Stokes I intensity at $\sim$700 GHz than that at $\sim$500 GHz, which suggests that emission at $\gtrsim$500 GHz is from optically thin part of a synchrotron emission spectrum. Both the Stokes I intensity and the polarization position angle at our highest observing frequency of $\sim$700 GHz, may be varying with time. However, we do not yet detect variation in the polarization percentage at $>$500 GHz. The polarization percentage at $\sim$700 GHz is likely lower than that at $\sim$500 GHz. By comparing the $\sim$500 GHz and $\sim$700 GHz observations with the observations at lower frequency bands, we suggest that the intrinsic polarization position angle of Sgr A* is varying with time. This paper also reports the measurable polarization properties from the observed calibration quasars. The future simultaneous multi-frequency polarization observations are required for clarifying the time and frequency variation of polarization position angle and polarization percentage.

Aims: Our goal is to study the chemical composition of the outflows of active galactic nuclei and starburst galaxies. Methods: We obtained high-resolution interferometric observations of HCN and HCO$^+$ $J=1\rightarrow0$ and $J=2\rightarrow1$ of the ultraluminous infrared galaxy Mrk~231 with the IRAM Plateau de Bure Interferometer. We also use previously published observations of HCN and HCO$^+$ $J=1\rightarrow0$ and $J=3\rightarrow2$, and HNC $J=1\rightarrow0$ in the same source. Results: In the line wings of the HCN, HCO$^+$, and HNC emission, we find that these three molecular species exhibit features at distinct velocities which differ between the species. The features are not consistent with emission lines of other molecular species. Through radiative transfer modelling of the HCN and HCO$^+$ outflow emission we find an average abundance ratio $X(\mathrm{HCN})/X(\mathrm{HCO}^+)\gtrsim1000$. Assuming a clumpy outflow, modelling of the HCN and HCO$^+$ emission produces strongly inconsistent outflow masses. Conclusions: Both the anti-correlated outflow features of HCN and HCO$^+$ and the different outflow masses calculated from the radiative transfer models of the HCN and HCO$^+$ emission suggest that the outflow is chemically differentiated. The separation between HCN and HCO$^+$ could be an indicator of shock fronts present in the outflow, since the HCN/HCO$^+$ ratio is expected to be elevated in shocked regions. Our result shows that studies of the chemistry in large-scale galactic outflows can be used to better understand the physical properties of these outflows and their effects on the interstellar medium (ISM) in the galaxy.

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.

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

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

The source NRAO150 is a very prominent millimeter to radio emitting quasar at redshift z=1.52 for which previous millimeter VLBI observations revealed a fast counterclockwise rotation of the innermost regions of the jet. Here we present new polarimetric multi-epoch VLBI-imaging observations of NRAO150 performed at 8, 15, 22, 43, and 86GHz with the VLBA and GMVA between 2006 and 2010. All new and previous observational evidence are consistent with an interpretation of the source behavior where the jet is seen at an extremely small angle to the line of sight, and the high frequency emitting regions in NRAO150 rotate at high speeds on the plane of the sky with respect to a reference point that does not need to be related to any particularly prominent jet feature. The observed polarization angle distribution at 22, 43, and 86 GHz during observing epochs with high polarization degree suggests that we have detected the toroidal component of the magnetic field threading the innermost jet plasma regions. This is also consistent with the lower degree of polarization detected at progressively poorer angular resolutions, where the integrated polarization intensity produced by the toroidal field is explained by polarization cancellation inside the observing beam. All this evidence is fully consistent with a kinematic scenario where the main kinematic and polarization properties of the 43 GHz emitting structure of NRAO150 are explained by the internal rotation of such emission regions around the jet axis when the jet is seen almost face on. A simplified model developed to fit helical trajectories to the observed kinematics of the 43 GHz features fully supports this hypothesis. This explains the kinematics of the innermost regions of the jet in NRAO150 in terms of internal jet rotation.

The main goal of this research is to determine the angular size and the atmospheric structures of cool giant stars and to compare them with hydrostatic stellar model atmospheres, to estimate the fundamental parameters, and to obtain a better understanding of the circumstellar environment. We conducted spectro-interferometric observations of epsilon Oct, beta Peg, NU Pav, and psi Peg in the near-infrared K band (2.13-2.47 microm), and gamma Hya (1.9-2.47 microm) with the VLTI/AMBER instrument at medium spectral resolution. To obtain the fundamental parameters, we compared our data with hydrostatic atmosphere models (PHOENIX). We estimated the Rosseland angular diameters of epsilon Oct, beta Peg, NU Pav, psi Peg, and gamma Hya. Together with distances and bolometric fluxes, we estimated radii, effective temperatures, and luminosities of our targets. In the beta Peg visibility, we observed a molecular layer of CO with a size similar to that modeled with PHOENIX. However, there is an additional slope in absorption starting around 2.3 microm. This slope is possibly due to a shell of water that is not modeled with PHOENIX. The visibility of psi Peg shows a low increase in the CO bands, compatible with the modeling of the PHOENIX model. The visibility data of the other sources show no increase in molecular bands. The spectra and visibilities predicted by the PHOENIX atmospheres agree with the spectra and the visibilities observed in our stars (except for beta Peg). This indicates that the opacity of the molecular bands is adequately included in the model, and the atmospheres of our targets have an extension similar to the modeled atmospheres. The atmosphere of beta Peg is more extended than that predicted by the model. The targets are located close to the red limits of the evolutionary tracks of the STAREVOL model.

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 present very long baseline interferometry (VLBI) observations of the milliarcsecond-scale radio structure of the supernova remnant SNR 4449$-$1 in the galaxy NGC 4449. This young and superluminous remnant was observed at 1.6 GHz ($\lambda = 18$\,cm) with the European VLBI Network. The observations confirm earlier identifications of this object with a supernova remnant (SNR) while revealing a somewhat different morphology compared with the structure reported by Bietenholz et al. from VLBI observations at 1.4 GHz. This difference is discussed here in the context of structural sensitivity of both observations. The 1.6 GHz image yields accurate estimates of the size (0.0422 arcsec $\times$ 0.0285 arcsec and 0.8 $\times$ 0.5 pc) and age ($\sim$55 yr) of SNR 4449$-$1. With a total flux of 6.1 $\pm$ 0.6 mJy measured in the VLBI image, the historical lightcurve of the source can be well represented by a power-law decay with a power index of $-$1.19 $\pm$ 0.07. The SNR exhibits a decline rate of the radio emission of 2.2$%$ $\pm$ 0.1$%$ yr$^{-1}$ and a radio luminosity of 1.74 $\times$ 10$^{35}$ erg s$^{-1}$.

CRL 618 is an object that exhibits characteristics of both AGB and post-AGB star. It also displays a spectacular array of bipolar lobes with a dense equatorial region, which makes it an excellent object to study the development of asymmetries in evolved stars. In the recent decades, an elliptical compact HII region located in the center of the nebula has been seen to be increasing in size and flux. This seems to be due to the ionization of the circumstellar envelope by the central star, and it would be indicating the beginning of the planetary nebula phase for CRL 618. We analyzed interferometric radio continuum data at ~5 and 22 GHz from observations carried out at seven epochs with the VLA. We traced the increase of the flux of the ionized region over a period of ~26 years. We measured the dimensions of the HII region directly from the brightness distribution images to determine the increase of its size over time. For one of the epochs we analyzed observations at six frequencies from which we estimated the electron density distribution. We carried out model calculations of the spectral energy distribution at two different epochs to corroborate our observational results. We found that the radio continuum flux and the size of the ionized region have been increasing monotonically in the last three decades. The size of the major axis of the HII region shows a dependance with frequency, which has been interpreted as a result of a gradient of the electron density in this direction. The growth of the HII region is due to the expansion of an ionized wind whose mass-loss rate increased continuously for a period of ~100 years until a few decades ago, when the mass-loss rate experienced a sudden decline. Our results indicate that the beginning of the ionization of the circumstellar envelope began around 1971, which marks the start of the planetary nebula phase of CRL 618.

The remarkable stability of extragalactic jets is surprising, given the reasonable possibility of the growth of instabilities. In addition, much work in the literature has invoked this possibility in order to explain observed jet structures and obtain information from these structures. For example, it was recently shown that the observed helical structures in the jet in S5 0836+710 could be associated with helical pressure waves generated by Kelvin-Helmholtz instability. Our aim is to resolve the arc-second structure of the jet in the quasar S5 0836+710 and confirm the lack of a hot-spot (reverse jet-shock) found by present observing arrays, as this lack implies a loss of jet collimation before interaction with the intergalactic medium. In this work, we use an observation performed in 2008 using EVN and MERLIN. The combined data reduction has provided a complete image of the object at arc-second scales. The lack of a hot-spot in the arc-second radio structure is taken as evidence that the jet losses its collimation between the VLBI region and the region of interaction with the ambient medium. This result, together with the previous identification of the helical structures in the jet with helical pressure waves that grow in amplitude with distance, allow us to conclude that the jet is probably disrupted by the growth of Kelvin-Helmholtz instability. This observational evidence confirms that the physical parameters of jets can be extracted using the assumption that instability is present in jets and can be the reason for many observed structures. Interestingly, the observed jet is classified as a FRII object in terms of its luminosity, but its large-scale morphology does not correspond to this classification. The implications of this fact are discussed.