• 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 have studied the polarization behavior of the quasar 3C 273 over the 1 mm wavelength band at ALMA with a total bandwidth of 7.5 GHz across 223 to 243 GHz at 0.8" resolution, corresponding to 2.1 kpc at the distance of 3C 273. With these observations we are able to probe the optically thin polarized emission close to the jet base, and constrain the magnetic field structure at the jet launching region. We compute the Faraday rotation measure using simple linear fitting and Faraday rotation measure synthesis. In addition, we model the broadband behavior of the fractional Stokes Q and U parameters. The systematic uncertainties in the polarization observations at ALMA are assessed through Monte Carlo simulations. We find the unresolved core of 3C 273 to be 3.3% linearly polarized. We detect a very high rotation measure (RM) of ~ 3.6x10$^5$rad/m$^2$ over the 1 mm band with all the methods we have used. This results in a rotation of >30 degrees of the intrinsic electric vector position angle, which is significantly higher than typically assumed for millimeter wavelengths. Our results strongly favor depolarization models over an external Faraday screen without depolarization. For a single region dominating the polarized emission, the RM and depolarization we observe can be explained by external Faraday dispersion or an ordered RM gradient within the synthesized beam. Further 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 ~1.9 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. (Abridged.)
  • 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 Atacama Millimeter/submillimeter Array (ALMA) Phasing Project (APP) has developed and deployed the hardware and software necessary to coherently sum the signals of individual ALMA antennas and record the aggregate sum in Very Long Baseline Interferometry (VLBI) Data Exchange Format. These beamforming capabilities allow the ALMA array to collectively function as the equivalent of a single large aperture and participate in global VLBI arrays. The inclusion of phased ALMA in current VLBI networks operating at (sub)millimeter wavelengths provides an order of magnitude improvement in sensitivity, as well as enhancements in u-v coverage and north-south angular resolution. The availability of a phased ALMA enables a wide range of new ultra-high angular resolution science applications, including the resolution of supermassive black holes on event horizon scales and studies of the launch and collimation of astrophysical jets. It also provides a high-sensitivity aperture that may be used for investigations such as pulsar searches at high frequencies. This paper provides an overview of the ALMA Phasing System design, implementation, and performance characteristics.
  • 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.
  • We use Atacama Large Millimeter/submillimeter Array Band 5 science verification observations of the red supergiant VY CMa to study the polarization of SiO thermal/masers lines and dust continuum at ~1.7 mm wavelength. We analyse both linear and circular polarization and derive the magnetic field strength and structure, assuming the polarization of the lines originates from the Zeeman effect, and that of the dust originates from aligned dust grains. We also discuss other effects that could give rise to the observed polarization. We detect, for the first time, significant polarization (~3%) of the circumstellar dust emission at millimeter wavelengths. The polarization is uniform with an electric vector position angle of $\sim8^\circ$. Varying levels of linear polarization are detected for the J=4-3 28SiO v=0, 1, 2, and 29SiO v=0, 1 lines, with the strongest polarization fraction of ~30% found for the 29SiO v=1 maser. The linear polarization vectors rotate with velocity, consistent with earlier observations. We also find significant (up to ~1%) circular polarization in several lines, consistent with previous measurements. We conclude that the detection is robust against calibration and regular instrumental errors, although we cannot yet fully rule out non-standard instrumental effects. Emission from magnetically aligned grains is the most likely origin of the observed continuum polarization. This implies that the dust is embedded in a magnetic field >13 mG. The maser line polarization traces the magnetic field structure. The magnetic field in the gas and dust is consistent with an approximately toroidal field configuration, but only higher angular resolution observations will be able to reveal more detailed field structure. If the circular polarization is due to Zeeman splitting, it indicates a magnetic field strength of ~1-3 Gauss, consistent with previous maser observations.
  • 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 have carried out a high-precision astrometric analysis of two very-long-baseline-interferometry (VLBI) epochs of observation of the 13 extragalactic radio sources in the complete S5 polar cap sample. The VLBI epochs span a time baseline of 10 years and enable us to achieve precisions in the proper motions of the source cores up to a few micro-arcseconds per year. The observations were performed at 14.4 GHz and 43.1 GHz, and enable us to estimate the frequency core-shifts in a subset of sources, for which the spectral-index distributions can be computed. We study the source-position stability by analysing the changes in the relative positions of fiducial source points (the jet cores) over a decade. We find motions of 0.1-0.9 mas among close-by sources between the two epochs, which imply drifts in the jet cores of approximately a few tens of micro-as per year. These results have implications for the standard Active Galactic Nucleus (AGN) jet model (where the core locations are supposed to be stable in time). For one of our sources, 0615+820, the morphological and spectral properties in year 2010, as well as the relative astrometry between years 2000 and 2010, suggest the possibility of either a strong parsec-scale interaction of the AGN jet with the ISM, a gravitational lens with ~1 mas diameter, or a resolved massive binary black hole.
  • Current mm/submm interferometers, like the Atacama Large mm/submm Array (ALMA), use receivers that register the sky signal in a linear polarization basis. In the case of observations performed in full-polarization mode (where the cross-correlations are computed among all the polarization channels) it is possible to reconstruct the full-polarization brightness distribution of the observed sources, as long as a proper calibration of delay offsets and leakage among polarization channels can be performed. Observations of calibrators, preferably with some linear polarization, with a good parallactic angle coverage are usually needed for such a calibration. In principle, dual-polarization observations only allow us to recover the Stokes $I$ intensity distribution of the sources, regardless of the parallactic angle coverage of the observations. In this paper, we present a novel technique of dual differential polarimetry that makes it possible to obtain information related to the full-polarization brightness distribution of the observed sources from dual-polarization observations. This technique is inspired in the Earth-rotation polarization synthesis and can be applied even to sources with complex structures.
  • We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($\alpha=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.
  • The Gamma-ray BL Lac object OJ 287 is known to exhibit inner-parsec "jet-wobbling", high degrees of variability at all wavelengths and quasi-stationary features including an apparent (~100 deg) position angle change in projection on the sky plane. Sub-50 micro-arcsecond resolution 86 GHz observations with the global mm-VLBI array (GMVA) supplement ongoing multi-frequency VLBI blazar monitoring at lower frequencies. Using these maps together with cm/mm total intensity and Gamma-ray observations from Fermi/LAT from 2008-2014, we aimed to determine the location of Gamma-ray emission and to explain the inner-mas structural changes. Observations with the GMVA offer approximately double the angular resolution compared with 43 GHz VLBA observations and allow us to observe above the synchrotron self-absorption peak frequency. The jet was spectrally decomposed at multiple locations along the jet. From this we derived estimates of the magnetic field. How the field decreases down the jet allowed an estimate of the distance to the jet apex and an estimate of the magnetic field strength at the jet apex and in the broad line region. Combined with accurate kinematics we attempt to locate the site of Gamma-ray activity, radio flares and spectral changes. Strong Gamma-ray flares appeared to originate from either the "core" region, a downstream stationary feature, or both, with Gamma-ray activity significantly correlated with radio flaring in the downstream quasi-stationary feature. Magnetic field estimates were determined at multiple locations along the jet, with the magnetic field found to be >1.6 G in the "core" and >0.4 G in the downstream quasi-stationary feature. We therefore found upper limits on the location of the "core" as >6.0 pc from the jet apex and determined an upper limit on the magnetic field near the jet base of the order of thousands of Gauss.
  • 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.
  • We discuss about how the relative flux densities among the images of gravitationally-lensed active galactic nuclei, AGN, can be used to study the intrinsic AGN variability with high accuracy. Multi-frequency monitoring observations of resolved gravitational lenses can allow us to detect signals of very weak variability and also provide information about the jet opacity and structure. As an example, we investigate the variability of the flux-density ratio between the two lensed images of the blazar B0218+357, using dual-frequency cm-wave observations. Similar to our previously reported submm-wave observations of the lensed blazar PKS1830-211, we observe a clear chromatic variability, starting short before an increase in the flux-density of the blazar. The evolution of the flux-density ratios between the blazar images shows a more clear and rich structure than that of the mere lightcurves of each individual image. The accuracy in the ratio measurements is allowing us to see variability episodes in the blazar that are weaker than the natural scatter in the absolute flux-density measurements. A simple opacity model in the jet is used to consistently explain the difference between the flux-density-ratio evolution at the two frequencies.
  • 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.
  • J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo, A. del Olmo, J.-F. Desmurs, J. M. Diego, R. Dominguez-Tenreiro, R. Estalella, A. Fernandez-Soto, E. Florido, J. Font, J. A. Font, A. Fuente, R. Garcia-Benito, S. Garcia-Burillo, B. Garcia-Lorenzo, A. Gil de Paz, J. M. Girart, J. R. Goicoechea, J. F. Gomez, M. Gonzalez-Garcia, O. Gonzalez-Martin, J. I. Gonzalez-Serrano, J. Gorgas, J. Gorosabel, A. Guijarro, J. C. Guirado, L. Hernandez-Garcia, C. Hernandez-Monteagudo, D. Herranz, R. Herrero-Illana, Y.-D. Hu, N. Huelamo, M. Huertas-Company, J. Iglesias-Paramo, S. Jeong, I. Jimenez-Serra, J. H. Knapen, R. A. Lineros, U. Lisenfeld, J. M. Marcaide, I. Marquez, J. Marti, J. M. Marti, I. Marti-Vidal, E. Martinez-Gonzalez, J. Martin-Pintado, J. Masegosa, J. M. Mayen-Gijon, M. Mezcua, S. Migliari, P. Mimica, J. Moldon, O. Morata, I. Negueruela, S. R. Oates, M. Osorio, A. Palau, J. M. Paredes, J. Perea, P. G. Perez-Gonzalez, E. Perez-Montero, M. A. Perez-Torres, M. Perucho, S. Planelles, J. A. Pons, A. Prieto, V. Quilis, P. Ramirez-Moreta, C. Ramos Almeida, N. Rea, M. Ribo, M. J. Rioja, J. M. Rodriguez Espinosa, E. Ros, J. A. Rubiño-Martin, B. Ruiz-Granados, J. Sabater, S. Sanchez, C. Sanchez-Contreras, A. Sanchez-Monge, R. Sanchez-Ramirez, A. M. Sintes, J. M. Solanes, C. F. Sopuerta, M. Tafalla, J. C. Tello, B. Tercero, M. C. Toribio, J. M. Torrelles, M. A. P. Torres, A. Usero, L. Verdes-Montenegro, A. Vidal-Garcia, P. Vielva, J. Vilchez, B.-B. Zhang
    June 17, 2015 astro-ph.IM
    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.
  • The Atacama Large mm-submm Array (ALMA) is, by far, the most sensitive mm/submm telescope in the World. The ALMA Phasing Project (APP) will allow us to phase-up all the ALMA antennas and use them as one single VLBI station. This will be a key component of the Event Horizon Telescope (EHT), a Global VLBI array at millimeter wavelengths. A problem in the APP is the calibration and conversion of the polarization channels. Most VLBI stations record their signals in a circular basis, but the ALMA receivers record in a linear basis. The strategy that will be followed in the phased-ALMA VLBI observations will be to correlate in "mixed" basis (i.e., linear versus circular) and convert the visibilities to a pure circular basis after the correlation. We have developed an algorithm to perform such a polarization conversion of the VLBI visibilities. In these proceedings, we present the basics of the PolConvert algorithm and discuss on the polarization conversion in the general case were single dishes (besides phased arrays) record with linear receivers in VLBI observations. We show some results of PolConvert applied to realistic simulations, as well as a test with real VLBI observations at 86\,GHz between the Onsala radiotelescope (recording in linear basis) and the Effelsberg radiotelescope (recording in circular basis).
  • The Global mm-VLBI Array (GMVA) is a network of 14 3\,mm and 7\,mm capable telescopes spanning Europe and the United States, with planned extensions to Asia. The array is capable of sensitive maps with angular resolution often exceeding 50\,$\mu$as. Using the GMVA, a large sample of prominent $\gamma$-ray blazars have been observed approximately 6 monthly from later 2008 until now. Combining 3\,mm maps from the GMVA with near-in-time 7\,mm maps from the VLBA-BU-BLAZAR program and 2\,cm maps from the MOJAVE program, we determine the sub-pc morphology and high frequency spectral structure of $\gamma$-ray blazars. The magnetic field strength can be estimated at different locations along the jet under the assumption of equipartition between magnetic field and relativistic particle energies. Making assumptions on the jet magnetic field configuration (e.g. poloidal or toroidal), we can estimate the separation of the mm-wave "core" and the jet base, and estimate the strength of the magnetic field there. The results of this analysis show that on average, the magnetic field strength decreases with a power-law $B \propto r^{-n}$, $n=0.3 \pm 0.2$. This suggests that on average, the mm-wave "core" is $\sim 1-3$\,pc downstream of the de-projected jet apex and that the magnetic field strength is of the order $B_{\rm{apex}} \sim 5-20$\,kG, broadly consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks (MAD)).
  • Context. ABDoradus is the main system of the ABDoradus moving group. It is a quadruple system formed by two widely separated binaries of pre-main-sequence (PMS) stars: ABDor A/C and ABDor Ba/Bb. The pair ABDor A/C has been extensively studied and its dynamical masses have been determined with high precision, thus making of ABDor C a benchmark for calibrating PMS stellar models. If the orbit and dynamical masses of the pair ABDor Ba/Bb can be determined, they could not only play a similar role to that of ABDor C in calibrating PMS models, but would also help to better understand the dynamics of the whole ABDoradus system. Aims. We aim to determine the individual masses of the pair ABDor Ba/Bb using VLBI observations and archive infrared data, as part of a larger program directed to monitor binary systems in the ABDoradus moving group. Methods. We observed the system ABDor B between 2007 and 2013 with the Australian Long Baseline Array (LBA), at a frequency of 8.4 GHz in phase-reference mode. Results. We detected, for the first time, compact radio emission from both stars in the binary, ABDor Ba and ABDor Bb. This result allowed us to determine the orbital parameters of both the relative and absolute orbits and, consequently, their individual dynamical masses: 0.28+/-0.05 Msun and 0.25+/-0.05 Msun, respectively. Conclusions. Comparisons of the dynamical masses with the prediction of PMS evolutionary models show that the models underpredict the dynamical masses of the binary components Ba and Bb by 10-30% and 10-40%, respectively, although they all still agree at the 2-sigma level. Some of stellar models considered favour an age between 50 and 100 Myr for this system, meanwhile others predict older ages. We also discuss the evolutionary status of ABDor Ba/Bb in terms of an earlier double-double star scenario that might explain the strong radio emission detected in both components.
  • Stacking as a tool for studying objects that are not individually detected is becoming popular even for radio interferometric data, and will be widely used in the SKA era. Stacking is typically done using imaged data rather than directly using the visibilities (the uv-data). We have investigated and developed a novel algorithm to do stacking using the uv-data. We have performed exten- sive simulations comparing to image-stacking, and summarize the results of these simulations. Furthermore, we disuss the implications in light of the vast data volume produced by the SKA. Having access to the uv-stacked data provides a great advantage, as it allows the possibility to properly analyse the result with respect to calibration artifacts as well as source properties such as size. For SKA the main challenge lies in archiving the uv-data. For purposes of robust stacking analysis, it would be strongly desirable to either keep the calibrated uv-data at least in an aver- age form, or implement a stacking queue where stacking positions could be provided prior to the observations and the uv-stacking is done almost in real time.
  • 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).
  • We report deep EVN and eMERLIN observations of the Type Ia SN 2014J in the nearby galaxy M 82. Our observations represent, together with JVLA observations of SNe 2011fe and 2014J, the most sensitive radio studies of Type Ia SNe ever. By combining data and a proper modeling of the radio emission, we constrain the mass-loss rate from the progenitor system of SN 2014J to $\dot{M} \lesssim 7.0\times 10^{-10}\, {\rm M_{\odot}\, yr^{-1}}$ (3-$\sigma$; for a wind speed of $100\, {\rm km s^{-1}}$). If the medium around the supernova is uniform, then $n_{\rm ISM} \lesssim 1.3 {\rm cm^3}$ (3-$\sigma$), which is the most stringent limit for the (uniform) density around a Type Ia SN. Our deep upper limits favor a double-degenerate (DD) scenario--involving two WD stars--for the progenitor system of SN 2014J, as such systems have less circumstellar gas than our upper limits. By contrast, most single-degenerate (SD) scenarios, i.e., the wide family of progenitor systems where a red giant, main-sequence, or sub-giant star donates mass to a exploding WD, are ruled out by our observations. Our estimates on the limits to the gas density surrounding SN 2011fe, using the flux density limits from Chomiuk et al. (2012), agree well with their results. Although we discuss possibilities for a SD scenario to pass observational tests, as well as uncertainties in the modeling of the radio emission, the evidence from SNe 2011fe and 2014J points in the direction of a DD scenario for both.
  • R.P.J. Tilanus, T.P. Krichbaum, J.A. Zensus, A. Baudry, M. Bremer, H. Falcke, G. Giovannini, R. Laing, H. J. van Langevelde, W. Vlemmings, Z. Abraham, J. Afonso, I. Agudo, A. Alberdi, J. Alcolea, D. Altamirano, S. Asadi, K. Assaf, P. Augusto, A-K. Baczko, M. Boeck, T. Boller, M. Bondi, F. Boone, G. Bourda, R. Brajsa, J. Brand, S. Britzen, V. Bujarrabal, S. Cales, C. Casadio, V. Casasola, P. Castangia, J. Cernicharo, P. Charlot, L. Chemin, Y. Clenet, F. Colomer, F. Combes, J. Cordes, M. Coriat, N. Cross, F. D'Ammando, D. Dallacasa, J-F. Desmurs, R. Eatough, A. Eckart, D. Eisenacher, S. Etoka, M. Felix, R. Fender, M. Ferreira, E. Freeland, S. Frey, C. Fromm, L. Fuhrmann, K. Gabanyi, R. Galvan-Madrid, M. Giroletti, C. Goddi, J. Gomez, E. Gourgoulhon, M. Gray, I. di Gregorio, R. Greimel, N. Grosso, J. Guirado, K. Hada, A. Hanslmeier, C. Henkel, F. Herpin, P. Hess, J. Hodgson, D. Horns, E. Humphreys, B. Hutawarakorn Kramer, V. Ilyushin, V. Impellizzeri, V. Ivanov, M. Julião, M. Kadler, E. Kerins, P. Klaassen, K. van 't Klooster, E. Kording, M. Kozlov, M. Kramer, A. Kreikenbohm, O. Kurtanidze, J. Lazio, A. Leite, M. Leitzinger, J. Lepine, S. Levshakov, R. Lico, M. Lindqvist, E. Liuzzo, A. Lobanov, P. Lucas, K. Mannheim, J. Marcaide, S. Markoff, I. Martí-Vidal, C. Martins, N. Masetti, M. Massardi, K. Menten, H. Messias, S. Migliari, A. Mignano, J. Miller-Jones, D. Minniti, P. Molaro, S. Molina, A. Monteiro, L. Moscadelli, C. Mueller, A. Müller, S. Muller, F. Niederhofer, P. Odert, H. Olofsson, M. Orienti, R. Paladino, F. Panessa, Z. Paragi, T. Paumard, P. Pedrosa, M. Pérez-Torres, G. Perrin, M. Perucho, D. Porquet, I. Prandoni, S. Ransom, D. Reimers, M. Rejkuba, L. Rezzolla, A. Richards, E. Ros, A. Roy, A. Rushton, T. Savolainen, R. Schulz, M. Silva, G. Sivakoff, R. Soria-Ruiz, R. Soria, M. Spaans, R. Spencer, B. Stappers, G. Surcis, A. Tarchi, M. Temmer, M. Thompson, J. Torrelles, J. Truestedt, V. Tudose, T. Venturi, J. Verbiest, J. Vieira, P. Vielzeuf, F. Vincent, N. Wex, K. Wiik, T. Wiklind, J. Wilms, E. Zackrisson, H. Zechlin
    July 1, 2014 astro-ph.IM
    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