Gravitational and plasma lensing share the same mathematical formalism in the limit of geometrical optics. Both phenomena can be effectively described by a projected, two-dimensional deflection potential whose gradient causes an instantaneous light deflection in a single, thin lens plane. We highlight the differences in the time-delay and lensing equations that occur because plasma lensing is caused by a potential directly proportional to the deflecting electron number density and gravitational lensing is caused by a potential related to the deflecting mass density by a Poisson equation. Since we treat plasma and gravitational lensing as thin-screen effective theories, their degeneracies are both caused by the unknown distribution of deflecting objects. Deriving the formalism-intrinsic degeneracies for plasma lensing, we find that they are analogous to those occurring in gravitational lensing. To break the degeneracies, galaxies and galaxy-cluster scale strong gravitational lenses must rely on additional assumptions or complementary observations. Physically realistic assumptions to arrive at self-consistent lens and source reconstructions can be provided by simulations and analytical effective theories. In plasma lensing, a deeper understanding of the deflecting electron density distributions is still under development, so that a model-based comprehensive lens reconstruction is not yet possible. However, we show that transient lenses and multi-wavelength observations help to break the arising degeneracies. We conclude that the development of an observation-based inference of local lens properties seems currently the best way to further probe the morphologies of plasma electron densities. Due to the simpler evidence-based breaking of the lensing degeneracies, we expect to obtain tighter constraints on the local plasma electron densities than on the gravitationally deflecting masses.

The shapes of individual self-gravitating structures of an ensemble of identical, collisionless particles have remained elusive for decades. In particular, a reason why mass density profiles like the Navarro-Frenk-White or the Einasto profile are good fits to simulation- and observation-based dark matter halos has not been found. Given the class of three dimensional, spherically symmetric power-law probability density distributions to locate individual particles in the ensemble mentioned above, we derive the constraining equation for the power-law index for the most and least likely joint ensemble configurations. We find that any dark matter halo can be partitioned into three regions: a core, an intermediate part, and an outskirts part up to boundary radius $r_\mathrm{max}$. The power-law index of the core is determined by the mean radius of the particle distribution within the core. The intermediate region becomes isothermal in the limit of infinitely many particles. The slope of the mass density profile far from the centre is determined by the choice of $r_\mathrm{max}$ with respect to the outmost halo particle, such that two typical limiting cases arise that explain the $r^{-3}$-slope for galaxy-cluster outskirts and the $r^{-4}$-slope for galactic outskirts. Hence, we succeed in deriving the mass density profiles of common fitting functions from a general viewpoint. These results also allow to find a simple explanation for the cusp-core-problem and to separate the halo description from its dynamics.

In our current best cosmological model, the vast majority of matter in the Universe is dark, consisting of yet undetected, non-baryonic particles that do not interact electro-magnetically. So far, the only significant evidence for dark matter has been found in its gravitational interaction, as observed in galaxy rotation curves or gravitational lensing effects. The inferred dark matter agglomerations follow almost universal mass density profiles that can be reproduced well in simulations, but have eluded an explanation from a theoretical viewpoint. Forgoing standard (astro-)physical methods, I show that it is possible to derive these profiles from an intriguingly simple mathematical approach that directly determines the most likely spatial configuration of a self-gravitating ensemble of collisionless dark matter particles.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The information about the mass density of galaxy clusters provided by the gravitational lens effect has inspired many inversion techniques. In this article, updates to the previously introduced method in Grale are described, and explored in a number of examples. The first looks into a different way of incorporating time delay information, not requiring the unknown source position. It is found that this avoids a possible bias that leads to "over-focusing" the images, i.e. providing source position estimates that lie in a considerably smaller region than the true positions. The second is inspired by previous reconstructions of the cluster of galaxies MACS J1149.6+2223, where a multiply-imaged background galaxy contained a supernova, SN Refsdal, of which four additional images were produced by the presence of a smaller cluster galaxy. The inversion for the cluster as a whole, was not able to recover sufficient detail interior to this quad. We show how constraints on such different scales, from the entire cluster to a single member galaxy, can now be used, allowing such small scale substructures to be resolved. Finally, the addition of weak lensing information to this method is investigated. While this clearly helps recover the environment around the strong lensing region, the mass sheet degeneracy may make a full strong and weak inversion difficult, depending on the quality of the ellipticity information at hand. We encounter ring-like structure at the boundary of the two regimes, argued to be the result of combining strong and weak lensing constraints, possibly affected by degeneracies.

The primary advantage of moderately superheated bubble chamber detectors is their simultaneous sensitivity to nuclear recoils from WIMP dark matter and insensitivity to electron recoil backgrounds. A comprehensive analysis of PICO gamma calibration data demonstrates for the first time that electron recoils in C$_3$F$_8$ scale in accordance with a new nucleation mechanism, rather than one driven by a hot-spike as previously supposed. Using this semi-empirical model, bubble chamber nucleation thresholds may be tuned to be sensitive to lower energy nuclear recoils while maintaining excellent electron recoil rejection. The PICO-40L detector will exploit this model to achieve thermodynamic thresholds as low as 2.8 keV while being dominated by single-scatter events from coherent elastic neutrino-nucleus scattering of solar neutrinos. In one year of operation, PICO-40L can improve existing leading limits from PICO on spin-dependent WIMP-proton coupling by nearly an order of magnitude for WIMP masses greater than 3 GeV c$^{-2}$ and will have the ability to surpass all existing non-xenon bounds on spin-independent WIMP-nucleon coupling for WIMP masses from 3 to 40 GeV c$^{-2}$.

We determine the cosmic expansion rate from supernovae of type Ia to set up a data-based distance measure that does not make assumptions about the constituents of the universe, i.e. about a specific parametrisation of a Friedmann cosmological model. The scale, determined by the Hubble constant $H_0$, is the only free cosmological parameter left in the gravitational lensing formalism. We investigate to which accuracy and precision the lensing distance ratio $D$ is determined from the Pantheon sample. Inserting $D$ and its uncertainty into the lensing equations for given $H_0$, esp. the time-delay equation between a pair of multiple images, allows to determine lens properties, esp. differences in the lensing potential ($\Delta \phi$), without specifying a cosmological model. We expand the luminosity distances into an analytic orthonormal basis, determine the maximum-likelihood weights for the basis functions by a globally optimal $\chi^2$-parameter estimation, and derive confidence bounds by Monte-Carlo simulations. For typical strong lensing configurations between $z=0.5$ and $z=1.0$, $\Delta \phi$ can be determined with a relative imprecision of 1.7%, assuming imprecisions of the time delay and the redshift of the lens on the order of 1%. With only a small, tolerable loss in precision, the model-independent lens characterisation developed in this paper series can be generalised by dropping the specific Friedmann model to determine $D$ in favour of a data-based distance ratio. Moreover, for any astrophysical application, the approach presented here, provides distance measures for $z\le2.3$ that are valid in any homogeneous, isotropic universe with general relativity as theory of gravity.

When light from a distant source object, like a galaxy or a supernova, travels towards us, it is deflected by massive objects that lie on its path. When the mass density of the deflecting object exceeds a certain threshold, multiple, highly distorted images of the source are observed. This strong gravitational lensing effect has so far been treated as a model-fitting problem. Using the observed multiple images as constraints yields a self-consistent model of the deflecting mass density and the source object. As several models meet the constraints equally well, we develop a lens characterisation that separates data-based information from model assumptions. The observed multiple images allow us to determine local properties of the deflecting mass distribution on any mass scale from one simple set of equations. Their solution is unique and free of model-dependent degeneracies. The reconstruction of source objects can be performed completely model-independently, enabling us to study galaxy evolution without a lens-model bias. Our approach reduces the lens and source description to its data-based evidence that all models agree upon, simplifies an automated treatment of large datasets, and allows for an extrapolation to a global description resembling model-based descriptions.

The galaxy-scale gravitational lens B0128+437 generates a quadrupole-image configuration of a background quasar that shows milli-arcsecond-scale subcomponents in the multiple images observed with VLBI. As this multiple-image configuration including the subcomponents has eluded a parametric lens-model characterisation so far, we determine local lens properties at the positions of the multiple images with our model-independent approach. Using PixeLens, we also succeed in setting up a global free-form mass density reconstruction including all subcomponents as constraints. We compare the model-independent local lens properties with those obtained by PixeLens and those obtained by the parametric modelling algorithm Lensmodel. A comparison of all three approaches and a model-free analysis based on the relative polar angles of the multiple images corroborate the hypothesis that elliptically symmetric models are too simplistic to characterise the asymmetric mass density distribution of this lenticular or late-type galaxy. In addition, the model-independent approach efficiently determines local lens properties on the scale of the quasar subcomponents, which are computationally intensive to obtain by free-form model-based approaches. As only 40% of the small-scale subcomponent local lens properties overlap within the 1-$\sigma$ confidence bounds, mass density gradients on milli-arcsecond scales cannot be excluded. Hence, aiming at a global reconstruction of the deflecting mass density distribution, increasingly detailed observations require flexible free-form models that allow for density fluctuations on milli-arcsecond scale to replace parametric ones, especially for asymmetric lenses or lenses with localised inhomogeneities like B0128.

We give a physical interpretation of the formalism intrinsic degeneracies of the gravitational lensing formalism that we derived on a mathematical basis in part IV of this series. We find that all degeneracies occur due to the partition of the mass density along the line of sight. Usually, it is partitioned into a background (cosmic) density and a foreground deflecting object. The latter can be further partitioned into a main deflecting object and perturbers. Weak deflecting objects along the line of sight are also added, either to the deflecting object or as a correction of the angular diameter distances, perturbing the cosmological background density. A priori, this is an arbitrary choice of reference frame and partition. They can be redefined without changing the lensing observables which are sensitive to the integrated deflecting mass density along the entire line of sight. Reformulating the time delay equation such that this interpretation of the degeneracies becomes easily visible, we note that the source can be eliminated from this formulation, which simplifies reconstructions of the deflecting mass distribution or the inference of the Hubble constant, $H_0$. Subsequently, we list necessary conditions to break the formalism intrinsic degeneracies and discuss ways to break them by model choices or including non-lensing observables like velocity dispersions along the line of sight with their advantages and disadvantages. We conclude with a systematic summary of all formalism intrinsic degeneracies and possibilities to break them.

We investigate to which precision local magnification ratios, $\mathcal{J}$, ratios of convergences, $f$, and reduced shears, $g = (g_{1}, g_{2})$, can be determined model-independently for the five resolved multiple images of the source at $z_\mathrm{s}=1.675$ in CL0024. We also determine if a comparison to the respective results obtained by the parametric modelling program Lenstool and by the non-parametric modelling program Grale can detect biases in the lens models. For these model-based approaches we additionally analyse the influence of the number and location of the constraints from multiple images on the local lens properties determined at the positions of the five multiple images of the source at $z_\mathrm{s}=1.675$. All approaches show high agreement on the local values of $\mathcal{J}$, $f$, and $g$. We find that Lenstool obtains the tightest confidence bounds even for convergences around one using constraints from six multiple image systems, while the best Grale model is generated only using constraints from all multiple images with resolved brightness features and adding limited small-scale mass corrections. Yet, confidence bounds as large as the values themselves can occur for convergences close to one in all approaches. Our results are in agreement with previous findings, supporting the light-traces-mass assumption and the merger hypothesis for CL0024. Comparing the three different approaches allows to detect modelling biases. Given that the lens properties remain approximately constant over the extension of the image areas covered by the resolvable brightness features, the model-independent approach determines the local lens properties to a comparable precision but within less than a second. (shortened)

We present a detailed description of the special procedures for calibration and quality assurance of Atacama Large Millimeter/submillimeter Array (ALMA) observations in Very Long Baseline Interferometry (VLBI) mode. These procedures are required to turn the phased ALMA array into a fully calibrated VLBI station. As an illustration of these methodologies, we present full-polarization observations carried out with ALMA as a phased array at 3mm (Band 3) and 1.3mm (Band 6) as part of Cycle-4. These are the first VLBI science observations conducted with ALMA and were obtained during a 2017 VLBI campaign in concert with other telescopes worldwide as part of the Global mm-VLBI Array (GMVA, April 1-3) and the Event Horizon Telescope (EHT, April 5-11) in ALMA Bands 3 and 6, respectively.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

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

With upcoming (continuum) surveys of high-resolution radio telescopes, detection rates of fast radio bursts (FRBs) might approach $10^5$ per sky per day by future extremely large observatories, such as the possible extension of the Square Kilometer Array (SKA) to a phase 2 array. Depending on the redshift distribution of FRBs and using the repeating FRB121102 as a model, we calculate a detection rate of multiply-imaged FRBs with their multiply-imaged hosts caused by the distribution of galaxy-cluster scale gravitational lenses of the order of $10^{-4}$ per square degree per year for a minimum total flux of the host of 10 $\mu$Jy at 1.4 GHz for SKA phase 2. Our comparison of estimated detection rates for quasars, supernovae, gamma ray bursts, and FRBs shows that multiple images of FRBs could be more numerous than those of gamma ray bursts and supernovae and as numerous as multiple images of quasars. Time delays between the multiple images of an FRB break degeneracies in model-based and model-independent lens reconstructions as other time-varying sources do, yet without a microlensing bias as FRBs are more point-like and have shorter duration times. We estimate the relative imprecision of FRB time-delay measurements to be $10^{-10}$ for time delays on the order of 100 days for galaxy-cluster scale lenses, yielding more precise (local) lens properties than time delays from the other time-varying sources. Using the lens modelling software Grale, we show the increase in accuracy and precision of the reconstructed scaled surface mass density map of a simulated cluster-scale lens when adding time delays for one set of multiple images to the set of observational constraints.

Based on the standard gravitational lensing formalism with its effective, projected lensing potential in a given background cosmology, we investigate under which transformations of the source position and of the deflection angle the observable properties of the multiple images, i.e. the time delay differences, the relative image positions, relative shapes, and magnification ratios, remain invariant. As these observables only constrain local lens properties, we derive general, local invariance transformations in the areas covered by the multiple images. We show that the known global invariance transformations, e.g. the mass sheet transformation or the source position transformation, are contained in our invariance transformations, when they are restricted to the areas covered by the multiple images and when lens-model-based degeneracies are ignored, like the freedom to add or subtract masses in unconstrained regions without multiple images. Hence, we have identified the general class of invariance transformations that can occur, in particular in our model-independent local characterisation of strong gravitational lenses.

In galaxy-galaxy strong gravitational lensing, Einstein rings are generated when the lensing galaxy has an axisymmetric lensing potential and the source galaxy is aligned with its symmetry centre along the line of sight. Using a Taylor expansion around the Einstein radius and eliminating the unknown source, I derive a set of analytic equations that determine differences of the deflection angle of the perturber weighted by the convergence of the axisymmetric lens and ratios of the convergences at the positions of the arcs from the measurable thickness of the arcs. In the same manner, asymmetries in the brightness distributions along an arc determine differences in the deflection angle of the perturber if the source has a symmetric brightness profile and is oriented parallel to or orthogonal to the caustic. These equations are the only model-independent information retrievable from observations to leading order in the Taylor expansion. General constraints on the derivatives of the perturbing lens are derived such that the perturbation does not change the number of critical curves. To infer physical properties such as the mass of the perturber or its position, models need to be inserted. The same conclusions about the scale of detectable masses and model-dependent degeneracies as in other approaches are then found and supported by analysing B1938 as an example. Yet, the model-independent equations show that there is a fundamental degeneracy between the main lens and the perturber that can only be broken if their relative position is known. This explains the degeneracies between lens models already found in simulations from a more general viewpoint. Depending on the properties of the pertuber, this degeneracy can be broken by characterising the surrounding of the lens or by measuring the time delay between quasar images embedded in the perturbed Einstein ring of the host galaxy.

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

Centaurus A (Cen A) is a bright radio source associated with the nearby galaxy NGC 5128 where high-resolution radio observations can probe the jet at scales of less than a light-day. The South Pole Telescope (SPT) and the Atacama Pathfinder Experiment (APEX) performed a single-baseline very-long-baseline interferometry (VLBI) observation of Cen A in January 2015 as part of VLBI receiver deployment for the SPT. We measure the correlated flux density of Cen A at a wavelength of 1.4 mm on a $\sim$7000 km (5 G$\lambda$) baseline. Ascribing this correlated flux density to the core, and with the use of a contemporaneous short-baseline flux density from a Submillimeter Array observation, we infer a core brightness temperature of $1.4 \times 10^{11}$ K. This is close to the equipartition brightness temperature, where the magnetic and relativistic particle energy densities are equal. Under the assumption of a circular Gaussian core component, we derive an upper limit to the core size $\phi = 34.0 \pm 1.8~\mu\textrm{as}$, corresponding to 120 Schwarzschild radii for a black hole mass of $5.5 \times 10^7 M_{\odot}$.

(shortened) We determine the transformation matrix T that maps multiple images with resolved features onto one another and that is based on a Taylor-expanded lensing potential close to a point on the critical curve within our model-independent lens characterisation approach. From T, the same information about the critical curve at fold and cusp points is derived as determined by the quadrupole moment of the individual images as observables. In addition, we read off the relative parities between the images, so that the parity of all images is determined, when one is known. We compare all retrievable ratios of potential derivatives to the actual ones and to those obtained by using the quadrupole moment as observable for two and three image configurations generated by a galaxy-cluster scale singular isothermal ellipse. We conclude that using the quadrupole moments as observables, the properties of the critical curve at the cusp points are retrieved to higher accuracy, at the fold points to lower accuracy, and the ratios of second order potential derivatives to comparable accuracy. We show that the approach using ratios of convergences and reduced shear is equivalent to ours close to the critical curve but yields more accurate results and is more robust because it does not require a special coordinate system like the approach using potential derivatives. T is determined by mapping manually assigned reference points in the images onto each other. If the assignment of reference points is subject to measurement uncertainties under noise, we find that the confidence intervals of the lens parameters can be as large as the values, when the uncertainties are larger than one pixel. Observed multiple images with resolved features are more extended than unresolved ones, so that higher order moments should be taken into account to improve the reconstruction.

We extend the model-independent approach to characterise strong gravitational lenses of Wagner & Bartelmann (2016) to its most general form to leading order by using the orientation angles of a set of multiple images with respect to their connection line(s) in addition to the relative distances between the images, their ellipticities and time-delays. For two symmetric images that straddle the critical curve, the orientation angle additionally allows to determine the slope of the critical curve and a second (reduced) flexion coefficient at the critical point on the connection line between the images. It also allows to drop the symmetry assumption that the axis of largest image extension is orthogonal to the critical curve. For three images almost forming a giant arc, the degree of assumed image symmetry is also reduced to the most general case, allowing to describe image configurations for which the source need not be placed on the symmetry axis of the two folds that unite at the cusp. For a given set of multiple images, we set limits on the applicability of our approach, show which information can be obtained in cases of merging images, and analyse the accuracy achievable due to the Taylor expansion of the lensing potential for the fold case on a galaxy cluster scale NFW-profile, a fold and cusp case on a galaxy cluster scale SIE-profile, and compare the generalised approach with the one of Wagner & Bartelmann (2016). The position of the critical points is reconstructed with less than 5" deviation for multiple images closer to the critical points than 30\% of the (effective) Einstein radius and the slope of the critical curve at a fold and its shape in the vicinity of a cusp deviate less than 20\% from the true values for distances of the images to the critical points less than 15\% of the (effective) Einstein radius.

As part of the Megamaser Cosmology Project (MCP), we present VLBI maps of nuclear water masers toward five galaxies. The masers originate in sub-parsec circumnuclear disks. For three of the galaxies, we fit Keplerian rotation curves to estimate their supermassive black hole (SMBH) masses, and determine (2.9 $\pm$ 0.3) $\times~10^{6}M_\odot$ for J0437+2456, (1.7 $\pm$ 0.1) $\times~10^{7}M_\odot$ for ESO 558$-$G009, and (1.1 $\pm$ 0.2) $\times~10^{7}M_\odot$ for NGC 5495. In the other two galaxies, Mrk 1029 and NGC 1320, the geometry and dynamics are more complicated and preclude robust black hole mass estimates. Including our new results, we compiled a list of 15 VLBI-confirmed disk maser galaxies with robust SMBH mass measurements. With this sample, we confirm the empirical relation of $R_{out} \propto 0.3 M_{SMBH}$ reported in Wardle & Yusef-Zadeh (2012). We also find a tentative correlation between maser disk outer radii and WISE luminosity. We find no correlations of maser disk size with X-ray 2-10 keV luminosity or [O III] luminosity.

The Galactic Center black hole Sagittarius A* (Sgr A*) is a prime observing target for the Event Horizon Telescope (EHT), which can resolve the 1.3 mm emission from this source on angular scales comparable to that of the general relativistic shadow. Previous EHT observations have used visibility amplitudes to infer the morphology of the millimeter-wavelength emission. Potentially much richer source information is contained in the phases. We report on 1.3 mm phase information on Sgr A* obtained with the EHT on a total of 13 observing nights over 4 years. Closure phases, the sum of visibility phases along a closed triangle of interferometer baselines, are used because they are robust against phase corruptions introduced by instrumentation and the rapidly variable atmosphere. The median closure phase on a triangle including telescopes in California, Hawaii, and Arizona is nonzero. This result conclusively demonstrates that the millimeter emission is asymmetric on scales of a few Schwarzschild radii and can be used to break 180-degree rotational ambiguities inherent from amplitude data alone. The stability of the sign of the closure phase over most observing nights indicates persistent asymmetry in the image of Sgr A* that is not obscured by refraction due to interstellar electrons along the line of sight.

We develop an approach to select families of lens models that can describe doubly and triply gravitationally lensed images near folds and cusps using the model-independent ratios of lensing-potential derivatives derived in Wagner & Bartelmann (2015). Models are selected by comparing these model-independent ratios of potential derivatives to (numerically determined) ratios of potential derivatives along critical curves for entire lens model families in a given range of parameter values. This comparison returns parameter ranges which lens model families can reproduce observation within, as well as sections of the critical curve where image sets of the observed type can appear. If the model-independent potential-derivative ratios inferred from the observation fall outside the range of these ratios derived for the lens model family, the entire family can be excluded as a feasible model in the given volume in parameter space. We employ this approach for the family of singular isothermal spheres with external shear to examples of lensing by a galaxy and two galaxy clusters (JVAS B1422+231, SDSS J2222+2745, and MACS J1149.5+2223) and show that the results obtained by our general method are in good agreement with results of previous model fits.

We develop a new approach to extracting model-independent information from observations of strong gravitational lenses. The approach is based on the generic properties of images near the fold and cusp catastrophes in caustics and critical curves. Observables used are the relative image positions, the magnification ratios and ellipticities of extended images, and time delays between images with temporally varying intensity. We show how these observables constrain derivatives and ratios of derivatives of the lensing potential near a critical curve. Based on these measured properties of the lensing potential, classes of parametric lens models can then easily be restricted to such parameter values compatible with the measurements, thus allowing fast scans of large varieties of models. Applying our approach to a representative galaxy (JVAS B1422+231) and a galaxy-cluster lens (MACS J1149.5+2223), we show which model-independent information can be extracted in those cases and demonstrate that the parameters obtained by our approach for known parametric lens models agree well with those found by detailed model fitting.

We report on 230 GHz (1.3 mm) VLBI observations of M87 with the Event Horizon Telescope using antennas on Mauna Kea in Hawaii, Mt. Graham in Arizona and Cedar Flat in California. For the first time, we have acquired 230 GHz VLBI interferometric phase information on M87 through measurement of closure phase on the triangle of long baselines. Most of the measured closure phases are consistent with 0$^{\circ}$ as expected by physically-motivated models for 230 GHz structure such as jet models and accretion disk models. The brightness temperature of the event-horizon-scale structure is $\sim 1 \times 10^{10}$ K derived from the compact flux density of $\sim 1$ Jy and the angular size of $\sim 40 $ $\rm \mu$as $\sim$ 5.5 $R_{{\rm s}}$, which is broadly consistent with the peak brightness of the radio cores at 1-86 GHz located within $\sim 10^2$ $R_{{\rm s}}$. Our observations occurred in the middle of an enhancement in very-high-energy (VHE) $\rm \gamma$-ray flux, presumably originating in the vicinity of the central black hole. Our measurements, combined with results of multi-wavelength observations, favor a scenario in which the VHE region has an extended size of $\sim$20-60 $R_{{\rm s}}$.

We report about a 230 GHz very long baseline interferometry (VLBI) fringe finder observation of blazar 3C 279 with the APEX telescope in Chile, the phased submillimeter array (SMA), and the SMT of the Arizona Radio Observatory (ARO). We installed VLBI equipment and measured the APEX station position to 1 cm accuracy (1 sigma). We then observed 3C 279 on 2012 May 7 in a 5 hour 230 GHz VLBI track with baseline lengths of 2800 M$\lambda$ to 7200 M$\lambda$ and a finest fringe spacing of 28.6 micro-arcseconds. Fringes were detected on all baselines with SNRs of 12 to 55 in 420 s. The correlated flux density on the longest baseline was ~0.3 Jy/beam, out of a total flux density of 19.8 Jy. Visibility data suggest an emission region <38 uas in size, and at least two components, possibly polarized. We find a lower limit of the brightness temperature of the inner jet region of about 10^10 K. Lastly, we find an upper limit of 20% on the linear polarization fraction at a fringe spacing of ~38 uas. With APEX the angular resolution of 230 GHz VLBI improves to 28.6 uas. This allows one to resolve the last-photon ring around the Galactic Center black hole event horizon, expected to be 40 uas in diameter, and probe radio jet launching at unprecedented resolution, down to a few gravitational radii in galaxies like M 87. To probe the structure in the inner parsecs of 3C 279 in detail, follow-up observations with APEX and five other mm-VLBI stations have been conducted (March 2013) and are being analyzed.

In recent years cosmology has undergone a revolution, with precise measurements of the microwave background radiation, large galaxy redshift surveys, and the discovery of the recent accelerated expansion of the Universe using observations of distant supernovae. In this light, the SKA enables us to do an ultimate test in cosmology by measuring the expansion rate of the Universe in real time. This can be done by a rather simple experiment of observing the neutral hydrogen (HI) signal of galaxies at two different epochs. The signal will encounter a change in frequency imprinted as the Universe expands over time and thus monitoring the drift in frequencies will provide a real time measure of the cosmic acceleration. Over a period of 12 years one would expected a frequency shift of the order of 0.1 Hz assuming a standard Lambda-CDM cosmology. Based on the sensitivity estimates of the SKA and the number counts of the expected HI galaxies, it is shown that the number counts are sufficiently high to compensate for the observational uncertainties of the measurements and hence allow a statistical detection of the frequency shift. [abstract abridged]

Aims. The phase scintillation of the European Space Agency's (ESA) Venus Express (VEX) spacecraft telemetry signal was observed at X-band (\lambda = 3.6 cm) with a number of radio telescopes of the European VLBI Network (EVN) in the period 2009-2013. Methods. We found a phase fluctuation spectrum along the Venus orbit with a nearly constant spectral index of -2.42 +/-0.25 over the full range of solar elongation angles from 0{\deg} to 45{\deg}, which is consistent with Kolmogorov turbulence. Radio astronomical observations of spacecraft signals within the solar system give a unique opportunity to study the temporal behaviour of the signal's phase fluctuations caused by its propagation through the interplanetary plasma and the Earth's ionosphere. This gives complementary data to the classical interplanetary scintillation (IPS) study based on observations of the flux variability of distant natural radio sources. Results. We present here our technique and the results on IPS. We compare these with the total electron content (TEC) for the line of sight through the solar wind. Finally, we evaluate the applicability of the presented technique to phase-referencing Very Long Baseline Interferometry (VLBI) and Doppler observations of currently operational and prospective space missions.

We report four new 22 GHz water masers found in a Green Bank Telescope search toward 37 nearby objects. Our goal was to find new maser galaxies, AGN disk-masers, and objects where both hydroxyl and water maser species coexist. We observed 37 sources within 250 Mpc that were selected by high X-ray luminosity (L_X > 10^40 W) and high absorbing column density (N_H ~ 10^22 cm^-2). Sources also included dual or triple AGN and interacting systems. We further included objects detected in hydroxyl (OH). The selection consisted of 16 new sources, 13 previous non-detections to follow up with a factor ten higher sensitivity, 10 OH masers and one deep OH absorber, of which 37 were observed. Water megamasers were detected towards the Sy 2 galaxy 2MFGC 13581 (a disk-maser), towards the 6 GHz OH absorber NGC 4261 (3C 270 with a twin-jet and dusty torus; broad water maser emission), and towards the two 1.6 GHz OH maser sources IRAS 17526+3253 and IRAS 20550+1656 (possible star formation water masers). We set upper limits on 33 non-detections. The search increased the number of known "dual-species" objects containing both OH and H2O masers to eight. The detection rate was 25 % in OH galaxies and 11 % overall. Combined with other searches, a total of 95 objects have now been searched for both OH and H2O masers. We found the overall dual-species detection rate (8 in 95) to be of the order of the joint probability of both species independently occurring in the same object (1 % lower bound). However, this needs to be verified by a more detailed analysis of selection criteria. Lastly, we see a lack of H2O kilomasers in OH megamaser objects already noted by Tarchi et al. (2011). This may be due to sensitivity bias rather than astrophysical reasons.

Global VLBI imaging at millimeter and sub-millimeter wavelength overcomes the opacity barrier of synchrotron self-absorption in AGN and opens the direct view into sub-pc scale regions not accessible before. Since AGN variability is more pronounced at short millimeter wavelength, mm-VLBI can reveal structural changes in very early stages after outbursts. When combined with observations at longer wavelength, global 3mm and 1mm VLBI adds very detailed information. This helps to determine fundamental physical properties at the jet base, and in the vicinity of super-massive black holes at the center of AGN. Here we present new results from multi-frequency mm-VLBI imaging of OJ287 during a major outburst. We also report on a successful 1.3mm VLBI experiment with the APEX telescope in Chile. This observation sets a new record in angular resolution. It also opens the path towards future mm-VLBI with ALMA, which aims at the mapping of the black hole event horizon in nearby galaxies, and the study of the roots of jets in AGN.