On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

Previous Very Long Baseline Interferometry (VLBI) observations of the quasar B1152+199 at 5GHz has revealed two images of a strongly lensed jet with seemingly discordant morphologies. Whereas the jet appears straight in one of the images, the other exhibits slight curvature on milliarcsecond scales. This is unexpected from the lensing solution and has been interpreted as possible evidence for secondary, small-scale lensing (millilensing) by a compact object with a mass of $~10^5$-$10^7\ M_\odot$ located close to the curved image. The probability for such a superposition is extremely low unless the millilens population has very high surface number density. Here, we revisit the case for millilensing in B1152+199 by combining new global-VLBI data at 8.4GHz with two datasets from the European VLBI Network (EVN) at 5GHz (archival) and at 22GHz (new dataset), and the previously published 5GHz Very Long Baseline Array (VLBA) data. We find that the new data with a more circular synthesized beam, exhibits no apparent milliarcsecond-scale curvature in image B. Various observations of the object spanning $\sim$15 years apart enable us to improve the constraints on lens system (thanks also to the improved astrometry resulting from 22GHz observations) to the point that the only plausible explanation left for the apparent curvature is the artifact due to the shape of the synthesized beam.

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

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

The blazar 1156+295 was active at gamma-ray energies, exhibiting three prominent flares during the year 2010. Here, we present results using the combination of broadband (X-ray through mm single dish) monitoring data and radio band imaging data at 43 GHz on the connection of gamma-ray events to the ejections of superluminal components and other changes in the jet of 1156+295. The kinematics of the jet over the interval 2007.0-2012.5 using 43 GHz Very Long Baseline Array observations, reveal the presence of four moving and one stationary component in the inner region of the blazar jet. The propagation of the third and fourth components in the jet corresponds closely in time to the active phase of the source in gamma rays. We briefly discuss the implications of the structural changes in the jet for the mechanism of gamma-ray production during bright flares. To localise the gamma-ray emission site in the blazar, we performed the correlation analysis between the 43 GHz radio core and the gamma-ray light curve. The time lag obtained from the correlation constrains the gamma-ray emitting region in the parsec-scale jet.

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

We have developed a new Bayesian Markov Chain Monte Carlo method to deconvolve light curves of blazars into individual flares, including proper estimation of the fit errors. We use the method to fit 15\,GHz light curves obtained within the OVRO 40-m blazar monitoring program where a large number of AGN have been monitored since 2008 in support of the Fermi Gamma-Ray Space Telescope mission. The time scales obtained from the fitted models are used to calculate the variability brightness temperature of the sources. Additionally, we have calculated brightness temperatures of a sample of these objects using Very Long Baseline Array data from the MOJAVE survey. Combining these two data sets enables us to study the intrinsic brightness temperature distribution in these blazars at 15 GHz. Our preliminary results indicate that the mean intrinsic brightness temperature in a sample of 14 sources is near the equipartition brightness temperature of $\sim 10^{11}$K.

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

Dark halo substructure may reveal itself through secondary, small-scale gravitational lensing effects on light sources that are macrolensed by a foreground galaxy. Here, we explore the prospects of using Very Long Baseline Interferometry (VLBI) observations of multiply-imaged quasar jets to search for submilliarcsecond-scale image distortions produced by various forms of dark substructures in the 1e3-1e8 Msolar mass range. We present lensing simulations relevant for the angular resolutions attainable with the existing European VLBI Network (EVN), the global VLBI array, and an upcoming observing mode in which the Atacama Large Millimeter Array (ALMA) is connected to the global VLBI array. While observations of this type would not be sensitive to standard cold dark matter subhalos, they can be used to detect more compact forms of halo substructure predicted in alternative structure formation scenarios. By mapping ~5 strongly lensed systems, it should be possible to detect or robustly rule out primordial black holes in the 1e3-1e6 Msolar mass range if they constitute >1% percent of the dark matter in these lenses. Ultracompact minihalos are harder to detect using this technique, but 1e6-1e8 Msolar ultracompact minihalos could in principle be detected if they constitute >10% of the dark matter.

The blazar OJ287 is the most promising (and the only) case for an extragalactic binary black hole system inspiralling under the action of gravitational radiation reaction. At present, though it is not possible to directly observe the binary components, it is possible to observe the jet emanating form the primary black hole. We argue that the orbital motion of the secondary black hole is reflected in the wobble of the jet and demonstrate that the wobble is orbital position dependent. The erratic wobble of the jet, reported in Agudo et al. (2012), is analyzed by taking into account the binary nature of the system and we find that the erratic component of jet wobble is very small.

We study the variation of the optical polarization angle in the blazar OJ287 and compare it with the precessing binary black hole model with a 'live' accretion disk. First, a model of the variation of the jet direction is calculated, and the main parameters of the model are fixed by the long term optical brightness evolution. Then this model is compared with the variation of the parsec scale radio jet position angle in the sky. Finally, the variation of the polarization angle is calculated using the same model, but using a magnetic field configuration which is at a constant angle relative to the optical jet. It is found that the model fits the data reasonably well if the field is almost parallel to the jet axis. This may imply a steady magnetic field geometry, such as a large-scale helical field.

We study the jet kinematics of the blazar S5 0716+714 during its active state in 2003-2004 with multi-epoch VLBI observations. Aims. We present a kinematic analysis of the large-scale (0-12 mas) jet of 0716+714, based on the results of six epochs of VLBA monitoring at 5 GHz. Additionally, we compare kinematic results obtained with two imaging methods based on different deconvolution algorithms. The blazar 0716+714 has a diffuse large-scale jet and a very faint bright compact core. Experiments with simulated data showed that the conventional data reduction procedure based on the CLEAN deconvolution algorithm does not perform well in restoring this type of structure. This might be the reason why previous kinematic studies of this source yielded ambiguous results. In order to obtain accurate kinematics of this source, we independently applied two imaging techniques to the raw data: the conventional method, based on difference mapping, which uses CLEAN deconvolution, and the generalized maximum entropy method (GMEM) realized in the VLBImager package developed at the Pulkovo Observatory in Russia. The results of both methods give us a consistent kinematic scenario: the large-scale jet of 0716+714 is diffuse and stationary. Differences between the inner (0-1 mas) and outer (1-12 mas) regions of the jet in brightness and velocity of the components could be explained by the bending of the jet, which causes the angle between the jet direction and the line of sight to change from ~5 deg to ~11 deg. For the source 0716+714 both methods worked at the limit of their capability.

The dark halo substructures predicted by current cold dark matter simulations may in principle be detectable through strong-lensing image splitting of quasars on small angular scales (0.01 arcseconds or below). Here, we estimate the overall probabilities for lensing by substructures in a host halo closely aligned to the line of sight to a background quasar. Under the assumption that the quasar can be approximated as a point source, the optical depth for strong gravitational lensing by subhalos typically turns out to be very small (tau < 0.01), contrary to previous claims. We therefore conclude that it is currently not feasible to use this strategy to put the simulation predictions for the dark matter subhalo population to the test. However, if one assumes the source to be spatially extended, as is the case for a quasar observed at radio wavelengths, there is a reasonable probability for witnessing substructure lensing effects even at rather large projected distances from the host galaxy, provided that the angular resolution is sufficient. While multiply-imaged, radio-loud quasars would be the best targets for unambiguously detecting dark matter subhalos, even singly-imaged radio quasars might be useful for setting upper limits on the abundance and central surface mass density of subhalos.

We observed the blazar \object{0716+714} with the VLBA during its active state in 2003-2004. In this paper we discuss multi-frequency analysis of the inner jet (first 1 mas) kinematics. The unprecedentedly dense time sampling allows us to trace jet components without misidentification and to calculate the component speeds with good accuracy. In the smooth superluminal jet we were able to identify and track three components over time moving outwards with relatively high apparent superluminal speeds (8.5-19.4 $c$), which contradicts the hypothesis of a stationary oscillating jet in this source. Component ejections occur at a relatively high rate (once in two months), and they are accompanied by mm-continuum outbursts. Superluminal jet components move along wiggling trajectories, which is an indication of actual helical motion. Fast proper motion and rapid decay of the components suggest that this source should be observed with the VLBI at a rate of at least once in one or two months in order to trace superluminal jet components without confusion.

In February 1997, the Japanese radio astronomy satellite HALCA was launched to provide the space-bourne element for the VLBI Space Observatory Programme (VSOP) mission. A significant fraction of the mission time was to be dedicated to the VSOP Survey of bright compact Active Galactic Nuclei (AGN) at 5 GHz, which was lead by ISAS. The VSOP Survey Sources are an unbiased dataset of 294 targets, of which 82% were successfully observed. These are now undergoing statistical analysis to tease out the characteristics of typical AGN sources. We present here the summary of the imaging and conclusions we have reached.

We present three 43 GHz images and a single 86 GHz image of Markarian 501 from VLBA observations in 2005. The 86 GHz image shows a partially resolved core with a flux density of about 200 mJy and a size of about 300 Schwarzschild radii, similar to recent results by Giroletti et al. Extreme limb-brightening is found in the inner parsec of the jet in the 43 GHz images, providing strong observational support for a `spine-layer' structure at this distance from the core. The jet is well resolved transverse to its axis, allowing Gaussian model components to be fit to each limb of the jet. The spine-layer brightness ratio and relative sizes, the jet opening angle, and a tentative detection of superluminal motion in the layer are all discussed.

Multi-frequency VLBI observations allow studies of the continuum spectrum in the different parts of the parsec scale jets of AGN, providing information on the physical properties of the plasma and magnetic fields in them. Since VLBI networks cannot be scaled, the range of spatial frequencies observed differs significantly between the different observing frequencies, which makes it difficult to obtain a broadband spectrum of the individual emission features in the jet. In this paper we discuss a model-fitting based spectral extraction method, which can significantly relieve this problem. The method uses a priori knowledge of the source structure, measured at high frequencies, to allow at lower frequencies the derivation of the sizes and flux densities of even those emission features that have mutual separations significantly less than the Rayleigh limit at the given frequency. We have successfully used this method in the analysis of 5-86 GHz VLBA data of 3C273. The spectra and sizes of several individual jet features were measured, thus allowing derivation of the magnetic flux density and the energy density of the relativistic electrons in the different parts of the jet. We discuss the results, which include e.g. a detection of a strong gradient in the magnetic field across the jet of 3C273.

In the CDM scenario, dark matter halos are assembled hierarchically from smaller subunits. A long-standing problem with this picture is that the number of sub-halos predicted by CDM simulations is orders of magnitudes higher than the known number of satellite galaxies in the vicinity of the Milky Way. A plausible way out of this problem could be that the majority of these sub-halos somehow have so far evaded detection. If such "dark galaxies" do indeed exist, gravitational lensing may offer one of the most promising ways to detect them. Dark matter sub-halos in the 1e6 - 1e10 solar mass range should cause strong gravitational lensing on (sub)milliarcsecond scales. We study the feasibility of a strong lensing detection of dark sub-halos by deriving the image separations expected for density profiles favoured by recent simulations and comparing these to the angular resolution of both existing and upcoming observational facilities. We find that there is a reasonable probability to detect sub-halo lensing effects in high resolution observations at radio wavelengths, such as produced by the upcoming VSOP-2 satellite, and thereby test the existence of dark galaxies.

The cold dark matter scenario predicts that a large number of dark subhalos should be located within the halo of each Milky-way sized galaxy. One tell-tale signature of such dark subhalos could be additional milliarcsecond-scale image splitting of quasars previously known to be multiply-imaged on arcsecond scales. Here, we estimate the image separations for the subhalo density profiles favoured by recent N-body simulations, and compare these to the angular resolution of both existing and upcoming observational facilities. We find, that the image separations produced are very sensitive to the exact subhalo density profile assumed, but in all cases considerably smaller than previous estimates based on the premise that subhalos can be approximated by singular isothermal spheres. Only the most optimistic subhalo models produce image separations that would be detectable with current technology, and many models produce image separations that will remain unresolved with all telescopes expected to become available in the foreseeable future. Detections of dark subhalos through image-splitting effects will therefore be far more challenging than currently believed, albeit not necessarily impossible.

According to the cold dark matter scenario, a large number of dark subhalos should be located within the halo of each Milky-way sized galaxy. One promising possibility for detecting such subhalos is to try to observe their gravitational lensing effects on background sources. Dark matter subhalos in the 10^6 - 10^10 M_solar mass range should cause strong gravitational lensing on the (sub)milliarcsecond scales, which can be observed only using space VLBI. We study the feasibility of a strong-lensing detection of dark subhalos by deriving the image separations expected for density profiles favoured by current simulations and comparing it to the angular resolution of both existing and upcoming observational facilities. We show that the detection of subhalos is likely much more difficult than suggested in previous studies, due to the smaller image separations predicted for subhalo density profiles more realistic than the singular isothermal sphere models often adopted.

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

We present first results from a multifrequency VLBA observations of 3C273 in 2003. The source was observed simultaneously at 5.0, 8.4, 15.3, 22.2, 43.2 and 86.2 GHz, and from this multifrequency data set, spectra of 16 emission features in the parsec scale jet were carefully constructed by using a new model-fitting based method. The measured spectra and sizes of the emission features were used to calculate the magnetic field density and the energy density of the relativistic electrons in the different parts of the parsec scale jet, independent of any equipartition assumption. We measure magnetic field density of an order of 1 Gauss in the core. The magnetic energy density in the core dominates over that of the relativistic electrons, while in the downstream region our data are roughly consistent with an equipartition. A strong gradient in the magnetic field density across the jet width, coincident with a transverse velocity structure at about 1.5 mas from the core, was found: the slower superluminal component B2 on the northern side of the jet has a magnetic field density two orders of magnitude lower than the faster southern components B3 and B4.

In February 1997 the Japanese radio astronomy satellite HALCA was launched to provide the space-borne element for the VSOP mission. HALCA provided linear baselines three-times greater than that of ground arrays, thus providing higher resolution and higher AGN brightness temperature measurements and limits. Twenty-five percent of the scientific time of the mission was devoted to the ``VSOP survey'' of bright, compact, extra-galactic radio sources at 5 GHz. A complete list of 294 survey targets were selected from pre-launch surveys, 91% of which were observed during the satellite's lifetime. The major goals of the VSOP Survey are statistical in nature: to determine the brightness temperature and approximate structure, to provide a source list for use with future space VLBI missions, and to compare radio properties with other data throughout the electromagnetic spectrum. All the data collected have now been analysed and is being prepared for the final image Survey paper. In these papers we present details of the mission, images of the sources, and some statistics of the contributions and completeness plus some statistics of the images and brightness temperatures. Full resolution images are to be found on http://www.vsop.isas.jaxa.jp/survey/publications/evn.ps.gz

We report the discovery of an extremely curved jet in the radio-loud quasar PKS2136+141. Multi-frequency Very Long Baseline Array (VLBA) images show a bending jet making a turn-around of 210 degrees in the plane of the sky, which is, to our knowledge, the largest ever observed change in the position angle of an astrophysical jet. Images taken at six different frequencies, from 2.3 to 43 GHz, reveal a spiral-like trajectory, which is likely a sign of an intrinsic helical geometry. A space-VLBI image, taken with the HALCA satellite at 5 GHz and having comparable resolution to our ground-based 15 GHz data, confirms that the bend is a frequency-independent structure. VLBA monitoring data at 15 GHz, covering eight years of observations, show knots in the jet clearly deviating from ballistic motion, which suggests that the bending may be caused by a growing helical Kelvin-Helmholtz normal mode. The jet appearance suggests a helical wave at a frequency well below the "resonant" frequency of the jet, which indicates that the wave is driven by a periodic perturbation at the base of the jet. We fit the observed structure in the source with a helical twist, and we find that a simple isothermal model with a constant wave speed and wavelength gives a good fit. The measured apparent velocities indicate some degree of acceleration along the jet, which together with an observed change in the apparent half-opening angle of the jet allow us to estimate the changes in the angle between the local jet direction and our line of sight. We suggest that the jet in PKS2136+141 is distorted by a helical Kelvin-Helmholtz normal mode externally driven into the jet (e.g. by precession), and that our line of sight falls within the opening angle of the helix cone.

In this first of a series of papers describing polarimetric multifrequency Very Long Baseline Array (VLBA) monitoring of 3C 273 during a simultaneous campaign with the INTEGRAL gamma-ray satellite in 2003, we present 5 Stokes I images and source models at 7 mm. We show that a part of the inner jet (1-2 milliarcseconds from the core) is resolved in a direction transverse to the flow, and we analyse the kinematics of the jet within the first 10 mas. Based on the VLBA data and simultaneous single-dish flux density monitoring, we determine an accurate value for the Doppler factor of the parsec scale jet, and using this value with observed proper motions, we calculate the Lorentz factors and the viewing angles for the emission components in the jet. Our data indicates a significant velocity gradient across the jet with the components travelling near the southern edge being faster than the components with more northern path. We discuss our observations in the light of jet precession model and growing plasma instabilities.

We report the discovery of an extremely curved jet in the high frequency peaking GPS quasar PKS 2136+141. Our multi-frequency VLBA images show a bending jet making a turn-around of 210 degrees on the plane of the sky, which is, to our knowledge, the largest ever observed change in a structural position angle of an extragalactic jet. Images taken at six different frequencies, from 2 to 43 GHz, beautifully reveal a spiral-like trajectory. We discuss possibilities to constrain the 3-D geometry of the source and suggest that it could be used as a testbed for models describing the bending of the relativistic jets.

We report first results of a multifrequency campaign from radio to hard X-ray energies of the prominent gamma-ray blazar 3C 279 during the first year of the INTEGRAL mission. The variable blazar was found at a low activity level, but was detected by all participating instruments. Subsequently a multifrequency spectrum could be compiled. The individual measurements as well as the compiled multifrequency spectrum are presented. In addition, this 2003 broadband spectrum is compared to one measured in 1999 during a high activity period of 3C 279.

Given the scientific goals of VSOP-2, including the possibility of observations of the shadows of black holes, we have investigated the fidelity of the recovered images given a typical uv-coverage. We find that we can achieve a dynamic range of better than 1000:1.

In February 1997 the Japanese radio astronomy satellite Halca was launched to provide the space-bourne element for the VSOP mission. Approximately twenty-five percent of the mission time has been dedicated to the VSOP Survey, a 5 GHz survey of bright, compact AGN. Both the final, calibrated, high resolution images and plots of visibility amplitude versus uv distance for the first 102 of the sources have been prepared and has been submitted. Papers on the methods and the models from fitting the cumulative {\em uv} amplitudes will also be submitted. The analysis of the second half is well underway.