Search

[Abridged] The relation between nuclear star clusters (NSCs) and the growth of the central SMBHs, as well as their connection to the properties of the host galaxies, is crucial for understanding the evolution of galaxies. Recent observations have revealed that about 10 per cent of nucleated galaxies host hybrid nuclei, consisting of both NSCs and accreting SMBHs that power active galactic nuclei (AGN). Motivated by the potential of the recently published multi-wavelength data sets from LeMMINGs survey, here we present the most thorough investigation to date of the incidence of hybrid nuclei in a large sample of 100 nearby nucleated galaxies (10 E, 25 S0, 63 S, and 2 Irr), covering a wide range in stellar mass ($M_{*,\rm gal} \sim 10^{8.7}-10^{12}~\rm M_{sun}$). We identify the nuclei and derive their properties by performing detailed 1D and 2D multi-component decompositions of the optical and near-infrared $HST$ stellar light distributions of the galaxies using S\'ersic and core-S\'ersic models. Our AGN diagnostics are based on homogeneously derived nuclear 1.5 GHz $e$-MERLIN radio, $Chandra$ X-ray (0.3--10 keV) and optical emission-line data. We determine the nucleation fraction ($f_{\rm nuc} $) as the relative incidence of nuclei across the LeMMINGs $HST$ sample and find $f_{\rm nuc} =~ $100/149 (= 67 $\pm$ 7 per cent), confirming previous work, with a peak value of 49/56~(= $88 \pm 13$ per cent) at bulge masses $M_{*,\rm bulge} \sim 10^{9.4}$- $10^{10.8}~\rm M_{sun}$. We identify 30 nucleated LeMMINGs galaxies that are optically active, radio-detected and X-ray luminous ($L_{X} > 10^{39}$ erg s$^{-1}$). This indicates that our nucleated sample has a lower limit $\sim$ 30 per cent occupancy of hybrid nuclei, which is a function of $M_{*,\rm bulge}$ and $M_{*,\rm gal}$. We find that hybrid nuclei have a number density of $(1.5 \pm 0.4) \times 10^{-5}$ Mpc$^{-3}$., Comment: Accepted for publication in MNRAS (21 pages, 10 figures, 9 tables)

JWST/NIRCam obtained high angular-resolution (0.05-0.1''), deep near-infrared 1--5 micron imaging of Supernova (SN) 1987A taken 35 years after the explosion. In the NIRCam images, we identify: 1) faint H2 crescents, which are emissions located between the ejecta and the equatorial ring, 2) a bar, which is a substructure of the ejecta, and 3) the bright 3-5 micron continuum emission exterior to the equatorial ring. The emission of the remnant in the NIRCam 1-2.3 micron images is mostly due to line emission, which is mostly emitted in the ejecta and in the hot spots within the equatorial ring. In contrast, the NIRCam 3-5 micron images are dominated by continuum emission. In the ejecta, the continuum is due to dust, obscuring the centre of the ejecta. In contrast, in the ring and exterior to the ring, synchrotron emission contributes a substantial fraction to the continuum. Dust emission contributes to the continuum at outer spots and diffuse emission exterior to the ring, but little within the ring. This shows that dust cooling and destruction time scales are shorter than the synchrotron cooling time scale, and the time scale of hydrogen recombination in the ring is even longer than the synchrotron cooling time scale. With the advent of high sensitivity and high angular resolution images provided by JWST/NIRCam, our observations of SN 1987A demonstrate that NIRCam opens up a window to study particle-acceleration and shock physics in unprecedented details, probed by near-infrared synchrotron emission, building a precise picture of how a SN evolves., Comment: Accepted for publication in MNRAS; 18 pages

[Abridged] We use high-resolution $HST$ imaging and $e$-MERLIN 1.5-GHz observations of galaxy cores from the LeMMINGs survey to investigate the relation between optical structural properties and nuclear radio emission for a large sample of galaxies. We perform accurate, multi-component decompositions of new surface brightness profiles extracted from $HST$ images for 163 LeMMINGs galaxies and fit up to six galaxy components (e.g., bulges, discs, AGN, bars, rings, spiral arms, and nuclear star clusters) simultaneously with S\'ersic and/or core-S\'ersic models. By adding such decomposition data for 10 LeMMINGs galaxies from our past work, the final sample of 173 nearby galaxies (102 Ss, 42 S0s, 23 Es plus 6 Irr) with bulge stellar mass (typically) M_*, bulge ~ 10^6-10^12.5 M_sun, encompasses all optical spectral classes (LINER, Seyfert, ALG and H II). We show that the bulge mass can be significantly overestimated in many galaxies when components such as bars, rings and spirals are not included in the fits. We additionally implement a Monte Carlo method to determine errors on bulge, disc and other fitted structural parameters. Moving (in the opposite direction) across the Hubble sequence, i.e., from the irregular to elliptical galaxies, we confirm that bulges become larger, more prominent and round. Such bulge dominance is associated with a brighter radio core luminosity. We also find that the radio detection fraction increases with bulge mass. At M_*,bulge > 10^11 M_sun, the radio detection fraction is 77%, declining to 24% for M_bulge < 10^10 M_sun. Furthermore, we observe core-S\'ersic bulges tend to be systematically round and to possess high radio core luminosities and boxy-distorted or pure elliptical isophotes., Comment: 102 pages, 26 figures, 10 tables, accepted for publication in A&A

Radio observations of ultracool dwarfs, objects comprising brown dwarfs and the very lowest mass stars, have mainly focused on analyzing their light-curve and spectral energy distributions providing valuable insights into their magnetic fields. However, spatially-resolved studies of such magnetospheres have been elusive so far. Radio interferometric observations of the brown dwarf LSR J1835+3259 reveal an extended magnetosphere with a morphology compatible with the presence of a radiation belt, similar to that of Jupiter and Earth, consisting of energetic particles confined via magnetic mirroring. Our finding suggests that radio emitting ultracool dwarfs may behave as scaled up versions of Jupiter, validating the connection between dipole-ordered magnetic fields and the presence of belt-like morphologies and aurorae beyond our Solar System.

What determines the nuclear radio emission in local galaxies? We combine optical [O III] line emission, robust black hole (BH) mass estimates, and high-resolution e-MERLIN 1.5-GHz data, from the LeMMINGs survey, of a statistically-complete sample of 280 nearby, optically active (LINER and Seyfert) and inactive HII and Absorption line galaxies [ALG]) galaxies. Using [O III] luminosity ($L_{\rm [O~III]}$) as a proxy for the accretion power, local galaxies follow distinct sequences in the optical-radio planes of BH activity, which suggest different origins of the nuclear radio emission for the optical classes. The 1.5-GHz radio luminosity of their parsec-scale cores ($L_{\rm core}$) is found to scale with BH mass ($M_{\rm BH}$) and [O~III] luminosity. Below $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, stellar processes from non-jetted HII galaxies dominate with $L_{\rm core} \propto M_{\rm BH}^{0.61\pm0.33}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.79\pm0.30}$. Above $M_{\rm BH} \sim$10$^{6.5}$ M$_{\odot}$, accretion-driven processes dominate with $L_{\rm core} \propto M_{\rm BH}^{1.5-1.65}$ and $L_{\rm core} \propto L_{\rm [O~III]}^{0.99-1.31}$ for active galaxies: radio-quiet/loud LINERs, Seyferts and jetted HII galaxies always display (although low) signatures of radio-emitting BH activity, with $L_{\rm 1.5\, GHz}\gtrsim$10$^{19.8}$ W Hz$^{-1}$ and $M_{\rm BH}\gtrsim10^{7}$ M$_{\odot}$, on a broad range of Eddington-scaled accretion rates ($\dot{m}$). Radio-quiet and radio-loud LINERs are powered by low-$\dot{m}$ discs launching sub-relativistic and relativistic jets, respectively. Low-power slow jets and disc/corona winds from moderately high to high-$\dot{m}$ discs account for the compact and edge-brightened jets of Seyferts, respectively. Jetted HII galaxies may host weakly active BHs. Fuel-starved BHs and recurrent activity account for ALG properties. [abridged], Comment: Accpeted for publication on MNRAS (27 pages, 9 Figures, 3 Tables, Appendix A)

Context: The fast rotator, pre-main sequence star AB Dor A is a strong and persistent radio emitter. The extraordinary coronal flaring activity is thought to be the origin of compact radio emission and other associated phenomena as large slingshot prominences. Aim: We aim to investigate the radio emission mechanism and the milliarcsecond radio structure around AB Dor A. Methods: We performed phase-referenced VLBI observations at 22.3 GHz, 8.4 GHz, and 1.4 GHz over more than one decade using the Australian VLBI array. Results: Our 8.4 GHz images show a double core-halo morphology, similar at all epochs, with emission extending at heights between 5 and 18 stellar radii. Furthermore, the sequence of the 8.4 GHz maps shows a clear variation of the source structure within the observing time. However, images at 1.4 GHz and 22.3 GHz are compatible with a compact source. The phase-reference position at 8.4 GHz and 1.4 GHz are coincident with those expected from the well-known milliarcsecond-precise astrometry of this star, meanwhile the 22.3 GHz position is 4$\sigma$ off the prediction in the north-west direction. The origin of this offset is still unclear. Conclusions: We have considered several models to explain the morphology and evolution of the inner radio structure detected in AB Dor A which include emission from the stellar polar caps, a flaring, magnetically-driven loop structure, and the presence of helmet streamers. A possible close companion to AB Dor A has been also investigated. Our results confirm the extraordinary coronal magnetic activity of this star, able to produce compact radio structures at very large heights, so far only seen in binary interacting systems., Comment: 11 pages, 6 figures, accepted for publication in Astronomy and Astrophysics

Context. Red supergiant stars possess surface features and extended molecular atmospheres. Photospheric convection may be a crucial factor of the levitation of the outer atmospheric layers. However, the mechanism responsible is still poorly understood. Aims. We image the stellar surface of V602 Carinae (V602 Car) to constrain the morphology and contrast of the surface features and of the extended atmospheric layers. Methods. We observed V602 Car with the Very Large Telescope Interferometer (VLTI) PIONIER instrument (1.53-1.78 $\mathrm{\mu}$m) between May and July 2016, and April and July 2019 with different telescope configurations. We compared the image reconstructions with 81 temporal snapshots of 3D radiative-hydrodynamics (RHD) CO$^5$BOLD simulations in terms of contrast and morphology, using the Structural Similarity Index. Results. The interferometric data are compatible with an overall spherical disk of angular diameter 4.4$\pm$0.2 mas, and an extended molecular layer. In 2016, the reconstructed image reveals a bright arc-like feature toward the northern rim of the photospheric surface. In 2019, an arc-like feature is seen at a different orientation and a new peak of emission is detected on the opposite side. The contrasts of the reconstructed surface images are 11%$\pm$2% and 9%$\pm$2% for 2016 and 2019, respectively. The morphology and contrast of the two images are consistent with 3D RHD simulations, within our achieved spatial resolution and dynamic range. The extended molecular layer contributes 10--13% of the total flux with an angular diameter of 6--8 mas. It is present but not clearly visible in the reconstructed images because it is close to the limits of the achieved dynamic range. The presence of the molecular layer is not reproduced by the 3D RHD simulations. Conclusions..., Comment: 16 pages, 13 figures, accepted for publication in Astronomy and Astrophysics

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., Comment: 8 pages, 5 figures, accepted for publication in ApJ Letters

We present high angular resolution (~80 mas) ALMA continuum images of the SN 1987A system, together with CO $J$=2 $\!\rightarrow\!$ 1, $J$=6 $\!\rightarrow\!$ 5, and SiO $J$=5 $\!\rightarrow\!$ 4 to $J$=7 $\!\rightarrow\!$ 6 images, which clearly resolve the ejecta (dust continuum and molecules) and ring (synchrotron continuum) components. Dust in the ejecta is asymmetric and clumpy, and overall the dust fills the spatial void seen in H$\alpha$ images, filling that region with material from heavier elements. The dust clumps generally fill the space where CO $J$=6 $\!\rightarrow\!$ 5 is fainter, tentatively indicating that these dust clumps and CO are locationally and chemically linked. In these regions, carbonaceous dust grains might have formed after dissociation of CO. The dust grains would have cooled by radiation, and subsequent collisions of grains with gas would also cool the gas, suppressing the CO $J$=6 $\!\rightarrow\!$ 5 intensity. The data show a dust peak spatially coincident with the molecular hole seen in previous ALMA CO $J$=2 $\!\rightarrow\!$ 1 and SiO $J$=5 $\!\rightarrow\!$ 4 images. That dust peak, combined with CO and SiO line spectra, suggests that the dust and gas could be at higher temperatures than the surrounding material, though higher density cannot be totally excluded. One of the possibilities is that a compact source provides additional heat at that location. Fits to the far-infrared--millimeter spectral energy distribution give ejecta dust temperatures of 18--23K. We revise the ejecta dust mass to $\mathrm{M_{dust}} = 0.2-0.4$M$_\odot$ for carbon or silicate grains, or a maximum of $<0.7$M$_\odot$ for a mixture of grain species, using the predicted nucleosynthesis yields as an upper limit., Comment: 32 pages, containing 19 figures and three appendices

Context. Colliding winds in massive binaries are able to accelerate particles up to relativistic speeds as the result of the interaction between the winds of the different stellar components. HD 167971 exhibits this phenomenology which makes it a strong radio source. Aims. We aim at characterizing the morphology of the radio emission and its dependence on the orbital motion, traced independently by NIR-interferometry, of the spectroscopic binary and the tertiary component that conforms HD 167971. Methods. We analyze 2006 and 2016 very long baseline interferometric data at C and X bands. We complement our analysis with a geometrical model of the wind-wind collision region, and with an astrometric description of the system. Results. We confirm that the detected non-thermal radio emission is associated with the wind-wind collision region of the spectroscopic binary and the tertiary component in HD 167971. The wind-wind collision region changes orientation in agreement with the orbital motion of the tertiary around the spectroscopic binary. The total intensity also changes between the two observing epochs in a way inversely proportional to the separation between the SB and T, with a negative-steep spectral index typical of an optically thin synchrotron emission possibly steepened by an inverse Compton cooling effect. The wind-wind collision bow-shock shape and its position with respect to the stars indicates that the wind momentum from the spectroscopic binary is stronger than that of the tertiary. Finally, the astrometric solution derived for the stellar system and the wind-wind collision region is consistent with independent Gaia data., Comment: Accepted for publication in Astronomy and Astrophysics, 7 pages, 6 figures

Despite more than 30 years of searches, the compact object in Supernova (SN) 1987A has not yet been detected. We present new limits on the compact object in SN 1987A using millimeter, near-infrared, optical, ultraviolet, and X-ray observations from ALMA, VLT, HST, and Chandra. The limits are approximately 0.1 mJy ($0.1\times 10^{-26}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) at 213 GHz, 1 Lsun ($6\times 10^{-29}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in optical if our line-of-sight is free of ejecta dust, and $10^{36}$ erg s$^{-1}$ ($2\times 10^{-30}$ erg s$^{-1}$ cm$^{-2}$ Hz$^{-1}$) in 2-10 keV X-rays. Our X-ray limits are an order of magnitude less constraining than previous limits because we use a more realistic ejecta absorption model based on three-dimensional neutrino-driven SN explosion models (presented in an accompanying article). The allowed bolometric luminosity of the compact object is 22 Lsun if our line-of-sight is free of ejecta dust, or 138 Lsun if dust-obscured. Depending on assumptions, these values limit the effective temperature of a neutron star to <4-8 MK and do not exclude models, which typically are in the range 3-4 MK. For the simplest accretion model, the accretion rate for an efficiency $\eta$ is limited to $< 10^{-11} \eta^{-1}$ Msun yr$^{-1}$, which excludes most predictions. For pulsar activity modeled by a rotating magnetic dipole in vacuum, the limit on the magnetic field strength ($B$) for a given spin period ($P$) is $B < 10^{14} P^2$ G s$^{-2}$. By combining information about radiation reprocessing and geometry, it is likely that the compact object is a dust-obscured thermally-emitting neutron star, which may appear as a region of higher-temperature ejecta dust emission., Comment: 35 pages, 10 figures. Accepted for publication in ApJ. Updated to match accepted version; slightly shortened and clarified the discussion

The star V766 Cen (=HR 5171A) was originally classified as a yellow hypergiant but lately found to more likely be a 27-36 Msun red supergiant (RSG). Recent observations indicated a close eclipsing companion in the contact or common-envelope phase. Here, we aim at imaging observations of V766 Cen to confirm the presence of the close companion. We used near-infrared H -band aperture synthesis imaging at three epochs in 2014, 2016, and 2017, employing the PIONIER instrument at the Very Large Telescope Interferometer (VLTI). The visibility data indicate a mean Rosseland angular diameter of 4.1+/-0.8 mas, corresponding to a radius of 1575+/-400 Rsun. The data show an extended shell (MOLsphere) of about 2.5 times the Rosseland diameter, which contributes about 30% of the H-band flux. The reconstructed images at the 2014 epoch show a complex elongated structure within the photospheric disk with a contrast of about 10%. The second and third epochs show qualitatively and quantitatively different structures with a single very bright and narrow feature and high contrasts of 20-30%. This feature is located toward the south-western limb of the photospheric stellar disk. We estimate an angular size of the feature of 1.7+/-0.3 mas, corresponding to a radius of 650+/-150 Rsun, and giving a radius ratio of 0.42+0.35/-0.10} compared to the primary stellar disk. We interpret the images at the 2016 and 2017 epochs as showing the close companion, or a common envelope toward the companion, in front of the primary. At the 2014 epoch, the close companion is behind the primary and not visible. Instead, the structure and contrast at the 2014 epoch are typical of a single RSG harboring giant photospheric convection cells. The companion is most likely a cool giant or supergiant star with a mass of 5+15/-3 Msun., Comment: 4 pages, 2 figures plus 4 pages appendix. Accepted for publication in Astronomy and Astrophysics (A&A) as a letter

Most massive stars end their lives in core-collapse supernova explosions and enrich the interstellar medium with explosively nucleosynthesized elements. Following core collapse, the explosion is subject to instabilities as the shock propagates outwards through the progenitor star. Observations of the composition and structure of the innermost regions of a core-collapse supernova provide a direct probe of the instabilities and nucleosynthetic products. SN 1987A in the Large Magellanic Cloud (LMC) is one of very few supernovae for which the inner ejecta can be spatially resolved but are not yet strongly affected by interaction with the surroundings. Our observations of SN 1987A with the Atacama Large Millimeter/submillimeter Array (ALMA) are of the highest resolution to date and reveal the detailed morphology of cold molecular gas in the innermost regions of the remnant. The 3D distributions of carbon and silicon monoxide (CO and SiO) emission differ, but both have a central deficit, or torus-like distribution, possibly a result of radioactive heating during the first weeks ("nickel heating"). The size scales of the clumpy distribution are compared quantitatively to models, demonstrating how progenitor and explosion physics can be constrained., Comment: ApJL accepted, 9 pages, 5 figures

We report the first molecular line survey of Supernova 1987A in the millimetre wavelength range. In the ALMA 210--300 and 340--360 GHz spectra, we detected cold (20--170 K) CO, 28SiO, HCO+ and SO, with weaker lines of 29SiO from ejecta. This is the first identification of HCO+ and SO in a young supernova remnant. We find a dip in the J=6--5 and 5--4 SiO line profiles, suggesting that the ejecta morphology is likely elongated. The difference of the CO and SiO line profiles is consistent with hydrodynamic simulations, which show that Rayleigh-Taylor instabilities cause mixing of gas, with heavier elements much more disturbed, making more elongated structure. We obtained isotopologue ratios of 28SiO/29SiO>13, 28SiO/30SiO>14, and 12CO/13CO>21, with the most likely limits of 28SiO/29SiO>128, 28SiO/30SiO>189. Low 29Si and 30Si abundances in SN 1987A are consistent with nucleosynthesis models that show inefficient formation of neutron-rich isotopes in a low metallicity environment, such as the Large Magellanic Cloud. The deduced large mass of HCO+ (~5x10^-6 Msun) and small SiS mass (<6x10^-5 Msun) might be explained by some mixing of elements immediately after the explosion. The mixing might have caused some hydrogen from the envelope to sink into carbon and oxygen-rich zones after the explosion, enabling the formation of a substantial mass of HCO+. Oxygen atoms may have penetrated into silicon and sulphur zones, suppressing formation of SiS. Our ALMA observations open up a new window to investigate chemistry, dynamics and explosive-nucleosynthesis in supernovae., Comment: Accepted by MNRAS, 16 pages

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., Comment: 15 pages, 19 figures, accepted for publication in A&A

We add four warmer late-type supergiants to our previous spectro-interferometric studies of red giants and supergiants. V766 Cen (=HR 5171 A) is found to be a high-luminosity log(L/L_sun)=5.8+-0.4 source of Teff 4290+-760 K and radius 1490+-540 Rsun located close to both the Hayashi and Eddington limits; this source is consistent with a 40 Msun evolutionary track without rotation and current mass 27-36 Msun. It exhibits NaI in emission arising from a shell of radius 1.5 Rphot and a photocenter displacement of about 0.1 Rphot. V766 Cen shows strong extended molecular (CO) layers and a dusty circumstellar background component. This suggest an optically thick pseudo-photosphere at about 1.5 Rphot at the onset of the wind. V766 Cen is a red supergiant located close to the Hayashi limit instead of a yellow hypergiant already evolving back toward warmer Teff as previously discussed. The stars sigma Oph, BM Sco, and HD 206859 are found to have lower luminosities of about log(L/Lsun)=3.4-3.5 and Teff of 3900-5300 K, corresponding to 5-9 Msun tracks. They do not show extended molecular layers as observed for higher luminosity red supergiants of our sample. BM Sco shows an unusually strong contribution by an over-resolved circumstellar dust component. These stars are more likely high-mass red giants instead of red supergiants. This leaves us with an unsampled locus in the HR diagram corresponding to luminosities log(L/Lsun)~3.8-4.8 or masses 10-13 Msun, possibly corresponding to the mass region where stars explode as type II-P supernovae during the RSG stage. Our previously found relation of increasing strength of extended molecular layers with increasing luminosities is now confirmed to extend to double our previous luminosities and up to the Eddington limit. This might further point to steadily increasing radiative winds with increasing luminosity. [Abridged], Comment: 16 pages, 14 figures, accepted for publication in Astronomy and Astrophysics (A&A)

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., Comment: Final version. In press

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., Comment: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Espa\~nola de Astronom\'ia

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., Comment: 9 pages, 4 tables, 5 figures. Accepted for publication by Astronomy & Astrophysics

We present the atmospheric structure and the fundamental parameters of three red supergiants, increasing the sample of RSGs observed by near-infrared spectro-interferometry. Additionally, we test possible mechanisms that may explain the large observed atmospheric extensions of RSGs. We carried out spectro-interferometric observations of 3 RSGs in the near-infrared K-band with the VLTI/AMBER instrument at medium spectral resolution. To comprehend the extended atmospheres, we compared our observational results to predictions by available hydrostatic PHOENIX, available 3-D convection, and new 1-D self-excited pulsation models of RSGs. Our near-infrared flux spectra are well reproduced by the PHOENIX model atmospheres. The continuum visibility values are consistent with a limb-darkened disk as predicted by the PHOENIX models, allowing us to determine the angular diameter and the fundamental parameters of our sources. Nonetheless, in the case of V602 Car and HD 95686, the PHOENIX model visibilities do not predict the large observed extensions of molecular layers, most remarkably in the CO bands. Likewise, the 3-D convection models and the 1-D pulsation models with typical parameters of RSGs lead to compact atmospheric structures as well, which are similar to the structure of the hydrostatic PHOENIX models. They can also not explain the observed decreases in the visibilities and thus the large atmospheric molecular extensions. The full sample of our RSGs indicates increasing observed atmospheric extensions with increasing luminosity and decreasing surface gravity, and no correlation with effective temperature or variability amplitude, which supports a scenario of radiative acceleration on Doppler-shifted molecular lines., Comment: Accepted for publication in A&A

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

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., Comment: 10 pages, 5 figures. Accepted for publication on The Astrophysical Journal

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., Comment: Accepted for publication in A&A

Precise determination of dynamical masses of pre-main-sequence (PMS) stars is essential to calibrate stellar evolution models that are widely used to derive theoretical masses of young low-mass objects. Binary stars in young, nearby loose associations are particularly good candidates for this calibration since all members share a common age. Interestingly, some of these young binaries present a persistent and compact radio emission, which makes them excellent targets for astrometric VLBI studies. We aim to monitor the orbital motion of the binary system HD 160934, a member of the AB Doradus moving group. We observed HD 160934 with the Very Large Array and the European VLBI Network at 8.4 and 5 GHz, respectively. The orbital information derived from these observations was analyzed along with previously reported orbital measurements. We show that the two components of the binary, HD 160934 A and HD 160934 c, display compact radio emission at VLBI scales, providing precise information on the relative orbit. Revised orbital elements were estimated. Future VLBI monitoring of this pair should determine precise model-independent mass estimates for the A and c components, which will serve as calibration tests for PMS evolutionary models., Comment: 5 pages, 5 figures, accepted for publication in A&A

The launching mechanism of the jets of active galactic nuclei is observationally poorly constrained, due to the large distances to these objects and the very small scales (sub-parsec) involved. In order to better constrain theoretical models, it is especially important to get information from the region close to the physical base of the jet, where the plasma acceleration takes place. In this paper, we report multi-epoch and multi-frequency continuum observations of the z=2.5 blazar PKS1830-211 with ALMA, serendipitously coincident with a strong $\gamma$-ray flare reported by Fermi-LAT. The blazar is lensed by a foreground z=0.89 galaxy, with two bright images of the compact core separated by 1". Our ALMA observations individually resolve these two images (although not any of their substructures), and we study the change of their relative flux ratio with time (four epochs spread over nearly three times the time delay between the two lensed images) and frequency (between 350 and 1050 GHz, rest-frame of the blazar), during the $\gamma$-ray flare. In particular, we detect a remarkable frequency-dependent behaviour of the flux ratio, which implies the presence of a chromatic structure in the blazar (i.e., a core-shift effect). We rule out the possiblity of micro- and milli-lensing effects and propose instead a simple model of plasmon ejection in the blazar's jet to explain the time and frequency variability of the flux ratio. We suggest that PKS1830-211 is likely one of the best sources to probe the activity at the base of a blazar's jet at submillimeter wavelengths, thanks to the peculiar geometry of the system. The implications of the core-shift in absorption studies of the foreground z=0.89 galaxy (e.g., constraints on the cosmological variations of fundamental constants) are discussed., Comment: Accepted for publication in A&A

SN1993J in M81 is the best studied young radio-luminous supernova in the Northern Hemisphere. We recently reported results from the analysis of a complete set of VLBI observations of this supernova at 1.7, 2.3, 5.0, and 8.4 GHz, covering a time baseline of more than one decade. Those reported results were focused on the kinematics of the expanding shock, the particulars of its evolving non-thermal emission, the density profile of the circumstellar medium, and the evolving free-free opacity by the supernova ejecta. In the present paper, we complete our analysis by performing a search for any possible signal from a compact source (i.e., a stellar-mass black hole or a young pulsar nebula) at the center of the expanding shell. We have performed a stacking of all our VLBI images at each frequency, after subtraction of our best-fit shell model at each epoch, and measured the peak intensity in the stacked residual image. Given the large amount of available global VLBI observations, the stacking of all the residual images allows us to put upper limits to the eventual emission of a putative compact central source at the level of $\sim102$ $\mu$Jy at 5 GHz (or, more conservatively, $\sim192$ $\mu$Jy, if we make a further correction for the ejecta opacity) and somewhat larger at other wavelengths., Comment: 4 pages, 3 figures. Accepted for publication in A&A

We present the atmospheric structure and the fundamental properties of the red supergiants (RSGs) AH Sco, UY Sct, and KW Sgr based on VLTI/AMBER observations. We carried out spectro-interferometric observations of AH Sco, UY Sct, and KW Sgr in the near-infrared K band with the VLTI/AMBER instrument, and compared the data to a new grid of hydrostatic PHOENIX model atmospheres. In our visibility data, we observe molecular layers of water and CO in extended atmospheres. For a uniform disk modeling, we observe size increases at the water band of 10% to 25% and at the CO bandheads of 20%-35% with respect to the near continuum bandpass. The PHOENIX atmosphere models predict the spectra and the continuum visibility values, but cannot reproduce the large extensions of the molecular layers. This indicates that the opacities of the molecular bands are included, but that the model atmospheres are too compact compared to the observations. The observed extended layers may be levitated by processes such as pulsation or convection, which are not included in the hydrostatic atmospheric models. Comparing the continuum visibility values to PHOENIX models, we estimate the Rosseland-mean photospheric angular diameters. Together with the distance and the spectro-photometry, we calculate radii and effective temperatures. The location of the targets in the HR-diagram is confirmed to be close to, and possibly slightly to the right of, the Hayashi limit of recent evolutionary tracks corresponding to masses between about 20 Msun and 40 Msun, Comment: Accepted for publication in A&A

In a recent publication, we reported results of a multi-frequency VLBI campaign of observations of the Active Galactic Nucleus (AGN) in galaxy M\,81, phase-referenced to the supernova SN\,1993J. We were able to extract precise information on the relative astrometry of the AGN radio emission at different epochs and frequencies. We found strong evidence of precession in the AGN jet (i.e., a systematic evolution in the jet inclination at each frequency) coupled to changes in the overall flux density at the different frequencies. In these proceedings, we summarise the main contents of our previous publication and we report on (preliminary) new results from our follow-up VLBI observations, now phase-referenced to the young supernova SN2008iz. We also briefly discuss how these results match the picture of our previously-reported precession model., Comment: 6 pages, 3 figures. Proceedings of EVN meeting 2012

We investigate the atmospheric structure and fundamental properties of the red supergiant VY CMa. We obtained near-infrared spectro-interferometric observations of VY CMa with spectral resolutions of 35 and 1500 using the AMBER instrument at the VLTI. The visibility data indicate the presence of molecular layers of water vapor and CO in the extended atmosphere with an asymmetric morphology. The uniform disk diameter in the water band around 2.0 mu is increased by \sim20% compared to the near-continuum bandpass at 2.20-2.25 mu and in the CO band at 2.3-2.5 mu it is increased by up to \sim50%. The closure phases indicate relatively small deviations from point symmetry close to the photospheric layer, and stronger deviations in the extended H2O and CO layers. Making use of the high spatial and spectral resolution, a near-continuum bandpass can be isolated from contamination by molecular and dusty layers, and the Rosseland-mean photospheric angular diameter is estimated to 11.3 +/- 0.3 mas based on a PHOENIX atmosphere model. Together with recent high-precision estimates of the distance and spectro-photometry, this estimate corresponds to a radius of 1420 +/- 120 Rsun and an effective temperature of 3490 +/- 90 K. VY CMa exhibits asymmetric, possibly clumpy, atmospheric layers of H2O and CO, which are not co-spatial, within a larger elongated dusty envelope. Our revised fundamental parameters put VY CMa close to the Hayashi limit of recent evolutionary tracks of initial mass 25 Msun with rotation or 32 Msun without rotation, shortly before evolving blueward in the HR-diagram., Comment: 5 pages, 5 figures, accepted for publication in Astronomy and Astrophysics (A&A) as a Letter

We report on the VLBI detection of supernova SN2011dh at 22GHz using a subset of the EVN array. The observations took place 14 days after the discovery of the supernova, thus resulting in a VLBI image of the youngest radio-loud supernova ever. We provide revised coordinates for the supernova with milli-arcsecond precision, linked to the ICRF. The recovered flux density is a factor 2 below the EVLA flux density reported by other authors at the same frequency and epoch of our observations. This discrepancy could be due to extended emission detected with the EVLA or to calibration problems in the VLBI and/or EVLA observations., Comment: Letter. Accepted in A&A

The pre-main sequence (PMS) star ABDorA is the main component of the quadruple system ABDoradus. The precise determination of the mass and photometry of the close companion to ABDorA, ABDorC, has provided an important benchmark for calibration of theoretical evolutionary models of low-mass stars. The limiting factor to the precision of this calibration is the age of the system, as both the mass and luminosity of ABDorA and C are well monitored by other ongoing programs. In this paper we present VLTI/AMBER observations of ABDorA which provide a direct measurement of the size of this star, 0.96+/-0.06 Rsun. The latter estimate, combined with other fundamental parameters also measured for this star, allows a precise test of PMS evolutionary models using both H-R diagrams and mass-radius relationships. We have found that our radius measurement is larger than that predicted by the models, which we interpret as an evidence of the oversizing produced by the strong magnetic activity of ABDorA. Considering, at least partially, this magnetic effect, theoretical isochrones have been used to derive constraints to the age of ABDorA, favouring an age about 40-50 Myr for this system. Older ages are not completely excluded by our data., Comment: 7 pages, 5 figures. Accepted for publication in A&A

(See the complete and formatted abstract in the paper). We report on VLBI monitoring of the low-luminosity AGN (LLAGN) in M81 at 1.7, 2.3, 5, and 8.4GHz. These observations are phase-referenced to the supernova SN1993J (located in the same galaxy) and cover from late 1993 to late 2005. The source consists at all frequencies of a slightly resolved core and a small jet extension towards the north-east direction (position angle of ~65 degrees) in agreement with previous publications. We find that the position of the intensity peak in the images at 8.4GHz is very stable in the galactic frame of M81 (proper motion upper limit about 0.010 mas per year). We confirm previous reports that the peaks at all frequencies are systematically shifted among them, possibly due to opacity effects in the jet as predicted by the standard relativistic jet model. We use this model to estimate the magnetic field in the jet and the mass of the central black hole. We obtain a black-hole mass of ~2.e+7 solar masses, comparable to estimates previously reported, but the magnetic fields obtained are thousands of times lower than previous estimates. We find that the positions of the cores at 1.7, 2.3, and 5GHz are less stable than that at 8.4GHz and evolve systematically, shifting southward at a rate of several tens of micro-arcsec per year. The evolution in the jet orientation seems to be related to changes in the inclination of the cores at all frequencies. These results can be interpreted as due to a precessing jet. The jet precession also seems to be related to a flare in the flux densities at 5.0 and 8.4GHz, which lasts ~4 years. A continued monitoring of the flux density and the jet structure evolution in this LLAGN will be necessary to further confirm our jet precession model., Comment: 16 pages; 13 figures; 8 tables. A&A (in press)

We report on K-band VLTI/AMBER observations at medium spectral resolution ($\sim$1500) of RS Capricorni, an M6/M7III semi-regular AGB star. From the spectrally-dispersed visibilities, we measure the star diameter as a function of observing wavelength from 2.13 to 2.47 microns. We derive a Rosseland angular diameter of $7.95 \pm 0.07$ mas, which corresponds to an effective temperature of $3160 \pm 160$ K. We detect size variations of around 10% in the CO band heads, indicating strong opacity effects of CO in the stellar photosphere. We also detect a linear increase of the size as a function of wavelength, beginning at 2.29 microns. Models of the stellar atmosphere, based on the mass of the star as estimated from stellar-evolution models, predict CO-size effects about half of those observed, and cannot reproduce the linear size increase with wavelength, redward of 2.29 microns. We are able to model this linear size increase with the addition of an extended water-vapor envelope around the star. However, we are not able to fit the data in the CO bandheads. Either the mass of the star is overestimated by the stellar-evolution models and/or there is an additional extended CO envelope in the outer part of the atmosphere. In any case, neither the water-vapor envelope, nor the CO envelope, can be explained using the current models., Comment: 8 pages, 6 figures. Accepted for publication in A&A

We have performed a complete re-calibration and re-analysis of all the available VLBI observations of supernova SN1993J, following an homogeneous and well-defined methodology. Observations of SN1993J at 69 epochs, spanning 13 years, were performed by two teams, which used different strategies and analysis tools. The results obtained by each group are similar, but their conclusions on the supernova expansion and the shape and evolution of the emitting region differ significantly. From our analysis of the combined set of observations, we have obtained an expansion curve with unprecedented time resolution and coverage. We find that the data from both teams are compatible when analyzed with the same methodology. One expansion index ($m_3 = 0.87 \pm 0.02$) is enough to model the expansion observed at 1.7\,GHz, while two expansion indices ($m_1 = 0.933\pm0.010$ and $m_2 = 0.796\pm0.005$), separated by a break time, $t_{br} = 390\pm30$ days, are needed to model the data, at frequencies higher than 1.7\,GHz, up to day 4000 after explosion. We thus confirm the wavelength dependence of the size of the emitting region reported by one of the groups. We also find that all sizes measured at epochs later than day 4000 after explosion are systematically smaller than our model predictions. We estimate the fractional shell width ($0.31 \pm 0.02$, average of all epochs and frequencies) and the level of opacity to the radio emission by the ejecta. We find evidence of a spectral-index radial gradient in the supernova shell, which is indicative of a frequency-dependent ejecta opacity. Finally, we study the distribution and evolution of the azimuthal anisotropies (hot spots) found around the radio shell during the expansion. These anisotropies have intensities of $\sim 20$% of the mean flux density of the shell, and appear to systematically evolve during the expansion., Comment: 13 pages, 9 figures, accepted for publication in A&A

We report on a simultaneous modelling of the expansion and radio light curves of SN1993J. We have developed a simulation code capable of generating synthetic expansion and radio light curves of supernovae by taking into consideration the evolution of the expanding shock, magnetic fields, and relativistic electrons, as well as the finite sensitivity of the interferometric arrays used in the observations. Our software successfully fits all the available radio data of SN 1993J with an standard emission model for supernovae extended with some physical considerations, as an evolution in the opacity of the ejecta material, a radial drop of the magnetic fields inside the radiating region, and a changing radial density profile of the circumstellar medium beyond day 3100 after explosion., Comment: 12 pages, 12 figures, accepted for publication in A&A

Phase referencing is a standard calibration procedure in radio interferometry. It allows to detect weak sources by using quasi-simultaneous observations of closeby sources acting as calibrators. Therefore, it is assumed that, for each antenna, the optical paths of the signals from both sources are similar. However, atmospheric turbulence may introduce strong differences in the optical paths of the signals and affect, or even waste, phase referencing for cases of relatively large calibrator-to-target separations and/or bad weather. The situation is similar in wide-field observations, since the random deformations of the images, mostly caused by atmospheric turbulence, have essentially the same origin as the random astrometric variations of phase-referenced sources with respect to the phase center of their calibrators. In this paper, we present the results of a Monte Carlo study of the astrometric precision and sensitivity of an interferometric array (a realization of the Square Kilometre Array, SKA) in phase-referenced and wide-field observations. These simulations can be extrapolated to other arrays by applying the corresponding corrections. We consider several effects from the turbulent atmosphere (i.e., ionosphere and wet component of the troposphere) and also from the antenna receivers. We study the changes in dynamic range and astrometric precision as a function of observing frequency, source separation, and strength of the turbulence. We find that, for frequencies between 1 and 10 GHz, it is possible to obtain images with high fidelity, although the atmosphere strongly limits the sensitivity of the instrument compared to the case with no atmosphere. Outside this frequency window, the dynamic range of the images and the accuracy of the source positions decrease. [...] (Incomplete abstract. Please read manuscript.), Comment: 9 pages, 11 figures. Accepted for publication in A&A.

Context: Phase referencing is a standard calibration technique in radio interferometry, particularly suited for the detection of weak sources close to the sensitivity limits of the interferometers. However, effects from a changing atmosphere and inaccuracies in the correlator model may affect the phase-referenced images, leading to wrong estimates of source flux densities and positions. A systematic observational study of signal decoherence in phase referencing, and its effects in the image plane, has not been performed yet. Aims: We systematically studied how the signal coherence in Very-Long-Baseline-Interferometry (VLBI) observations is affected by a phase-reference calibration at different frequencies and for different calibrator-to-target separations. The results obtained should be of interest for a correct interpretation of many phase-referenced observations with VLBI. Methods: We observed a set of 13 strong sources (the S5 polar cap sample) at 8.4 and 15 GHz in phase-reference mode, with 32 different calibrator/target combinations spanning angular separations between 1.5 and 20.5 degrees. We obtained phase-referenced images and studied how the dynamic range and peak flux density depend on observing frequency and source separation. Results: We obtained dynamic ranges and peak flux densities of the phase-referenced images as a function of frequency and separation from the calibrator. We compared our results with models and phenomenological equations previously reported. Conclusions: The dynamic range of the phase-referenced images is strongly limited by the atmosphere at all frequencies and for all source separations. The limiting dynamic range is inversely proportional to the sine of the calibrator-to-target separation. We also find that the peak flux densities, relative to those obtained with the self-calibrated images, decrease with source separation., Comment: 5 pages, 3 figures. Submitted to A&A on 5.02.2010; accepted on 11.03.2010.

We report on 1.6 GHz Very-Long-Baseline-Interferometry (VLBI) observations of supernova SN 1979C made on 18 November 2002. We derive a model-dependent supernova size. We also present a reanalysis of VLBI observations made by us on June 1999 and by other authors on February 2005. We conclude that, contrary to our earlier claim of strong deceleration in the expansion, SN 1979C has been undergoing almost-free expansion ($m = 0.91\pm0.09$; $R \propto t^m$) for over 25 years., Comment: 4 pages, 4 figures; submitted to A&A on 14 May 2009. Accepted on 7 Jul 2009

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type IIb supernova SN 2008ax, made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal detection of the supernova. The total flux density recovered from our VLBI image is 0.8$\pm$0.3 mJy (one standard deviation). As it appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48$\pm$0.12 mJy, compatible with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion velocities derived from our observations are $(1.90 \pm 0.30) \times 10^5$ km s$^{-1}$ (case of a double source) and $(5.2 \pm 1.3) \times 10^4$ km s$^{-1}$ (taking the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations., Comment: 4 pages, 3 figures. Accepted in A&A on 24/03/2009

We have studied the growth of the shell-like radio structure of supernova SN1993J in M81 from September 1993 through October 2003 with very-long-baseline interferometry (VLBI) observations at the wavelengths of 3.6, 6, and 18cm. For this purpose, we have developed a method to accurately determine the outer radius (R) of any circularly symmetric compact radio structure like SN1993J. The source structure of SN1993J remains circularly symmetric (with deviations from circularity under 2%) over almost 4000 days. We characterize the decelerated expansion of SN 1993J through approximately day 1500 after explosion with an expansion parameter $m= 0.845\pm0.005$ ($R \propto t^{m}$). However, from that day onwards the expansion is different when observed at 6 and 18cm. Indeed, at 18cm, the expansion can be well characterized by the same $m$ as before day 1500, while at 6cm the expansion appears more decelerated, and is characterized by another expansion parameter, $m_{6}= 0.788\pm0.015$. Therefore, since about day 1500 on, the radio source size has been progressively smaller at 6cm than at 18cm. These findings are in stark contrast to previous reports by other authors on the details of the expansion. In our interpretation the supernova expands with a single expansion parameter, $m= 0.845\pm0.005$, and the 6cm results beyond day 1500 are due to physical effects, perhaps also coupled to instrumental limitations. Two physical effects may be involved: (a) a changing opacity of the ejecta to the 6cm radiation, and (b) a radial decrease of the magnetic field in the emitting region. (Long abstract cut. Please, read full abstract in manuscript)., Comment: 21 pages, 19 figures, accepted in A&A

Prior to explosion, a supernova progenitor slowly loses significant amounts of its hydrogen envelope in a stellar wind. After the explosion, the blastwave interacts with this wind producing synchrotron emission. A year of radio observations allows us to probe the progenitor evolution for a thousand years. The EVLA and SKA would represent more than an order of magnitude improvement in our ability to explore the pre-explosion lives of a significantly large population of supernova progenitor stars. It will allow us to move beyond the crude optical classifications and develop a deeper physical understanding of how massive stars live and die., Comment: 8 pages, 4 figures. White Paper submitted to the Astro2010 SSE panel

Phase self-calibration (or selfcal) is an algorithm often used in the calibration of interferometric observations in astronomy. Although a powerful tool, this algorithm presents strong limitations when applied to data with a low signal-to-noise ratio. We analyze the artifacts that the phase selfcal algorithm produces when applied to extremely noisy data. We show how the phase selfcal may generate a spurious source in the sky from a distribution of completely random visibilities. This spurious source is indistinguishable from a real one. We numerically and analytically compute the relationship between the maximum spurious flux density generated by selfcal from noise and the particulars of the interferometric observations. Finally, we present two simple tests that can be applied to interferometric data for checking whether a source detection is real or whether the source is an artifact of the phase self-calibration algorithm., Comment: 7 pages, 2 figures. Accepted in A&A on 12 December 2007

We report on the first wide-field, high-precision astrometric analysis of the 13 extragalactic radio sources of the complete S5 polar cap sample at 15.4 GHz. We describe new algorithms developed to enable the use of differenced phase delays in wide-field astrometric observations and discuss the impact of using differenced phase delays on the precision of the wide-field astrometric analysis. From this global fit, we obtained estimates of the relative source positions with precisions ranging from 14 to 200 $\mu$as at 15.4 GHz, depending on the angular separation of the sources (from $\sim$1.6 to $\sim$20.8 degrees). These precisions are $\sim$10 times higher than the achievable precisions using the phase-reference mapping technique., Comment: 9 pages, 7 figures

We present our observations of the radio emission from supernova (SN) 1993J, in M 81 (NGC 3031), made with the VLA, from 90 to 0.7 cm, as well as numerous measurements from other telescopes. The combined data set constitutes probably the most detailed set of measurements ever established for any SN outside of the Local Group in any wavelength range. Only SN 1987A in the LMC has been the subject of such an intensive observational program. The radio emission evolves regularly in both time and frequency, and the usual interpretation in terms of shock interaction with a circumstellar medium (CSM) formed by a pre-SN stellar wind describes the observations rather well considering the complexity of the phenomenon. However: 1) The 85 - 110 GHz measurements at early times are not well fitted by the parameterization, unlike the cm wavelength measurements. 2) At mid-cm wavelengths there is some deviation from the fitted radio light curves. 3) At a time ~3100 days after shock breakout, the decline rate of the radio emission steepens without change in the spectral index. This decline is best described as an exponential decay starting at day 3100 with an e-folding time of ~1100 days. 4) The best overall fit to all of the data is a model including both non-thermal synchrotron self-absorption (SSA) and a thermal free-free absorbing (FFA) components at early times, evolving to a constant spectral index, optically thin decline rate, until the break in that decline rate. Moreover, neither a purely SSA nor a purely FFA absorbing models can provide a fit that simultaneously reproduces the light curves, the spectral index evolution, and the brightness temperature evolution. 5) The radio and X-ray light curves exhibit similar behavior and suggest a sudden drop in the SN progenitor mass-loss rate at ~8000 years prior to shock breakout., Comment: 45 pages, 13 figures, accepted for ApJ

We present the results from the on-going radio monitoring of recent type II supernovae (SNe), including SNe 2004et, 2004dj, 2002hh, 2001em, and 2001gd. Using the Very Large Array to monitor these supernovae, we present their radio light-curves. From these data we are able to discuss parameterizations and modeling and make predictions of the nature of the progenitors based on previous research. Derived mass loss rates assume wind-established circumstellar medium, shock velocity ~10,000 km/s, wind velocity ~10 km/s, and CSM Temperature ~10,000 K., Comment: 3 pages, 6 figures Conference Proceedings: "Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters" AIP, New York, eds. S. Immler, K.W. Weiler, and R. McCray

We present the results of radio observations, taken primarily with the Very Large Array, of Supernovae 1993J, 2001gd, 2001em, 2002hh, 2004dj, and 2004et. We have fit a parameterized model to the multi-frequency observations of each supernova. We compare the observed and derived radio properties of these supernovae by optical classification and discuss the implications., Comment: 5 pages, 2 figures, 1 table Conference Proceedings: "Supernova 1987A: 20 Years After: Supernovae and Gamma-Ray Bursters" AIP, New York, eds. S. Immler, K.W. Weiler, and R. McCray

We report on 8.4GHz Very Long Baseline Interferometry (VLBI) observations of the type II-P supernova SN2004et in the spiral galaxy NGC 6946, made on 20 February 2005 (151 days after explosion). The Very Large Array (VLA) flux density was 1.23$\pm$0.07 mJy, corresponding to an isotropic luminosity at 8.4GHz of (4.45$\pm$0.3)$\times10^{25}$ erg s$^{-1}$ Hz$^{-1}$ and a brightness temperature of (1.3$\pm$0.3)$\times10^{8}$ K. We also provide an improved source position, accurate to about 0.5 mas in each coordinate. The VLBI image shows a clear asymmetry. From model fitting of the size of the radio emission, we estimate a minimum expansion velocity of 15,700$\pm$2,000 km s$^{-1}$. This velocity is more than twice the expected mean expansion velocity estimated from a synchrotron self-absorbed emission model, thus suggesting that synchrotron self-absorption is not relevant for this supernova. With the benefit of an optical spectrum obtained 12 days after explosion, we favor an emission model which consists of two hot spots on an underlying expanding shell of width comparable to that of SN 1993J., Comment: Accepted for publication in A&A (22/05/07)

We present an astrometric analysis of the binary systems ABDorA /ABDorC and ABDorBa / ABDorBb. These two systems of well-known late-type stars are gravitationally associated and they constitute the quadruple ABDoradus system. From the astrometric data available at different wavelengths, we report: (i) a determination of the orbit of ABDorC, the very low mass companion to ABDorA, which confirms the mass estimate of 0.090Msun reported in previous works; (ii) a measurement of the parallax of ABDorBa, which unambiguously confirms the long-suspected physical association between this star and ABDorA; and (iii) evidence of orbital motion of ABDorBa around ABDorA, which places an upper bound of 0.4Msun on the mass of the pair ABDorBa / ABDorBb (50% probability). Further astrometric monitoring of the system at all possible wavelengths would determine with extraordinary precision the dynamical mass of its four components., Comment: 7 pages, 4 figures