Next Article in Journal
The Validation of FORMOSAT-3/COSMIC Measurements in the Middle Latitude Region of China with Ionosonde Observations during 2015–2018
Previous Article in Journal
A Covariant Polymerized Scalar Field in Semi-Classical Loop Quantum Gravity
Previous Article in Special Issue
An Updated View of the Milky Way from Maser Astrometry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Software and Techniques for VLBI Data Processing and Analysis

1
Max-Planck-Institut für Radioastronomie, 53121 Bonn, Germany
2
Department of Physics, University of Pretoria, Lynnwood Road, Hatfield, Pretoria 0083, South Africa
3
Jodrell Bank Centre for Astrophysics, University of Manchester, Oxford Road, Manchester M13 9PL, UK
4
National Institute for Theoretical and Computational Sciences (NITheCS), Matieland, Stellenbosch 7602, South Africa
*
Author to whom correspondence should be addressed.
Universe 2022, 8(10), 527; https://doi.org/10.3390/universe8100527
Submission received: 1 August 2022 / Revised: 8 September 2022 / Accepted: 13 September 2022 / Published: 10 October 2022
(This article belongs to the Special Issue VLBI Science Applications)

Abstract

:
Very-long-baseline interferometry (VLBI) is a challenging observational technique, which requires in-depth knowledge about radio telescope instrumentation, interferometry, and the handling of noisy data. The reduction in raw data is mostly left to the scientists and demands the use of complex algorithms implemented in comprehensive software packages. The correct application of these algorithms necessitates a good understanding of the underlying techniques and physics that are at play. The verification of the processed data produced by the algorithms demands a thorough understanding of the underlying interferometric VLBI measurements. This review describes the latest techniques and algorithms that scientists should know about when analyzing VLBI data.

1. Introduction

In recent years, several novel techniques and methods have been developed for the processing and analysis of very-long-baseline interferometry (VLBI) data, related primarily to the new needs and capabilities of novel VLBI instruments. For instance, the growing number of telescopes participating in large and global alliances of VLBI arrays (e.g., [1,2,3,4,5]), significant increases in bandwidths via recording rates with improved gigabit per second (Gbps) (e.g., [6,7]), multi-frequency receiver capabilities (e.g., [8]), advancements towards observing increasingly shorter wavelength emission to gain resolving power and counter self-absorption effects in the observed sources (e.g., [9,10]), and prospects of the future next-generation Very Large Array (ngVLA, [11,12]), VLBI capabilities of the Square Kilometre Array (SKA, [13,14]), next-generation Event Horizon Telescope (ngEHT, [15]), as well as millimeter space VLBI observatories (e.g., [16,17,18,19]).
Lately, the localization of a repeating fast radio burst [20] and black hole images [10,21] have been highlights of continuum total-intensity VLBI science. More complex VLBI methods, which require a few additional steps, enable the study of magnetic field signatures through polarized synchrotron radiation, spectral line emission from astronomical masers, and wide fields of view (FOVs) combined with the VLBI resolving power.
With polarimetric VLBI (e.g., [22,23,24,25,26]), the role that magnetic fields play in the acceleration, collimation, and morphology of active galactic nuclei (AGN) jets can be studied across a wide range of scales. At the smallest scales, parameters of the central object responsible for producing the jet can be inferred to test for the presence of an event horizon, measure the spin of the compact object, and assess the role of the magnetic field in the jet launching.
Through spectral line VLBI (e.g., [27,28,29,30,31,32,33,34,35,36]), maser emission from star-forming regions, H I absorption, ejected circumstellar envelopes from giant stars, and gas in the vicinity of AGN that may be part of inflows or outflows can be mapped. Through Doppler shifts, the kinematics and physical conditions of the sources can be studied. Furthermore, angular-diameter distances to megamaser disk galaxies can be used for cosmology.
Wide-field VLBI (e.g., [37,38,39,40,41,42]) enables powerful astrometric studies, the mapping of distances to pulsars in our galaxy, and the search for gravitationally lensed radio sources. Moreover, faint sources in the FOV can be detected, details in large extended structures such as AGN jet hotspots can be imaged when the instrument’s sensitivity is sufficient, and star formation radio emission can be distinguished from AGN activity in deep field studies.
Additionally, VLBI observations of spacecraft should be mentioned (e.g., [43], and references therein). With a sufficiently stable and strong onboard radio emitter, spacecraft can be tracked with high accuracy, allowing for a wide range of planetary science applications, gravimetric studies, and tests of general relativity, for example.
In this work, we review the typical sequence of VLBI processing steps and highlight the latest established open-source software packages and techniques. This information is targeted primarily at scientists who have obtained data from a VLBI observation and want an overview of the available tools for data processing and scientific analyses. We focus on the latest state-of-the-art algorithms and do not give a historical overview. We also do not discuss the scheduling of VLBI observations, which principal investigators (PIs) should do together with observatories and the help of the SCHED software 1.
For a broad theoretical VLBI background, the reader is referred to Thompson et al. [44]. Definitions of technical terms used in this work that new VLBI users might be unfamiliar with are given in Table 1.
VLBI measurements with the long baselines of the Low Frequency Array (LOFAR) are also possible but require unique calibration strategies that are beyond the scope of this work. For example, due to the prevalence of direction-dependent effects, very low observing frequencies, and large data volumes. For details, the reader is referred to Morabito et al. [45].
This review is structured as follows: Section 2 describes the correlation of VLBI data. Section 3 details the “signal stabilization”, which encompasses all data calibration steps needed before the data can be averaged in time and frequency. Section 4 outlines the amplitude gain corrections for telescopes that are needed for a flux density calibration in physical units. Section 5 explains how the sky brightness distribution of the observed source can be reconstructed from the calibrated data. Section 6 addresses how special-science cases of polarization, spectral line, and wide-field VLBI are to be handled. Use cases and algorithms for the synthetic generation of VLBI data are given in Section 7. A summary is presented in Section 8.

2. Correlation

PIs are usually provided with correlated data, but it is instructive for scientists to understand the fundamentals of interferometric observations and upstream processing applied to their data.
There are two types of astronomical interferometers: connected-element interferometers and very-long-baseline interferometric networks, which are used in astronomical and geodetic VLBI observations. Connected-element arrays can be “phased”, producing beamforming and phased-sum time-domain signal output (e.g., [3,47]). Such phased arrays can participate in VLBI and are equivalent to a single large dish with a narrow beam.
For interferometry, the recorded baseband data (Table 1) of all antennas are combined at a “correlator” computing system. In a connected-element interferometer, telescopes are connected over short (walkable) distances to a central correlator via RF-over-fiber or high-speed network links and share a common time and frequency standard. The short distances within an array permit very wide baseband signals (currently up to 64 GHz of bandwidth) from a large number of antennas (currently up to 10,000) to be transferred and processed at the array correlator in real time. The correlator usually employs hardware for the signal processing (FPGA, ASIC) due to power efficiency and operating cost considerations. With real-time correlation, the voluminous (often > 10 Tbit/s aggregate) baseband data do not need to be stored. Only the final time averaged data products are kept.
VLBI networks consist of relatively few (≲ 20) telescopes, and unlike in most connected-element arrays, they are often heterogeneous with different frequency tuning constraints, digital baseband bandwidths, baseband data formats, and no shared time and frequency standard. The large distances between telescopes located across the globe and in orbits around the Earth for space VLBI make real-time data transfer challenging.
Although there are recent VLBI hardware correlators, e.g., the 16-station × 8 Gbps Korean VLBI Network (KVN, [48,49,50]) and VLBI Exploration of Radio Astrometry (VERA, [51]) correlator at the Korea–Japan Correlation Center [52], software correlators are currently more common. The flexible signal processing in software can better address possible complications that arise from the heterogeneous set of telescopes. In addition, new observing modes are more rapidly implemented in software correlators.

2.1. Software Correlators

Two typical open-source software VLBI correlators are DiFX 2 [53,54] and SFXC 3 [55]. DiFX has a wide developer community and user base including the National Radio Astronomy Observatory (NRAO)/Very Long Baseline Array [VLBA, [56]], the Event Horizon Telescope (EHT, [57]), and geodetic services. SFXC is more native to the European VLBI Network (EVN; see [58], and references therein) environment.
Being “embarrassingly parallel”, the correlation process is highly suited for cluster computing. Recently, graphics processing units (GPUs) have started to be used. One production GPU correlator is the six-station × 16 Gbps near-realtime RasFX of the Quasar VLBI network [59,60]. Other GPU-based VLBI correlators are under development or in a demonstrator stage, including an Extended-KVN GPU correlator and a possible correlator for VLBI space science [61]. A GPU-accelerated version of DiFX is under development as well 4.

2.2. Station Clock Model

Typically, each VLBI telescope operates its own frequency standard, with an Allan deviation ([46], Table 1) of < 10 13  s/1 s and < 10 13  s/100 s to surpass atmospheric turbulence [62,63]. Receivers and backends are locked to the reference, and backends time-tag their digital baseband output data. An absolute time reference is provided by GPS, with an uncertainty of around ± 20  ns due to ionospheric diurnal variation. Possibilities of sharing clock signals over long distances are currently being investigated (e.g., [63,64]).
Without a shared clock signal, the time-frequency transfer is emulated by a clock model that is applied at the correlator. The model consists of station clock offsets and drifts over time of the station frequency standards. Clock offsets are the differences measured between the backend data time-tag and the absolute Universal Time (UT). Drifts are measured via a continuous comparison against the GPS absolute time, usually fit to first order (linear drift) or sometimes to a higher order (acceleration) in space VLBI [65].
Software correlators tend to be more flexible with their clock models than their hardware counterparts, handling offsets of several seconds at sub-microsecond granularity and drifts of 10 9  s/1 s or larger.

2.3. Processing at the Correlator

The inputs for the correlator data processing are the baseband signals E r ( t ) from all telescopes. Appendix A describes how the recorded data are shipped to the correlation facility. Processing involves trivial computations at a high data rate. It produces time averaged mutual coherences between pairs of telescopes, analytically given by
Γ r 1 , r 2 ( u , v , τ ) = lim T 1 2 T T T E r 1 ( t ) E r 2 * ( t τ ) d t ,
at τ = 0 , over an integration time T, for telescopes at locations r 1 and r 2 ; see also Equation (15.1) in [44]. Coherences between all telescope pairs and across their different spatial separations, or  ( u , v ) -plane spatial frequencies, are called interferometric “visibilities” (cf. § 4.1 of [44]).
According to the van Cittert-Zernike theorem, these visibilities are the Fourier counterpart of the angular source brightness distribution, I ( l , m ) , related through
Γ r 1 , r 2 ( u , v , 0 ) = I ( l , m ) exp 2 π i ( u l + v m ) d l d m ;
see Equations (3.7) and (15.7) in [44]. Here, the line-of-sight directions across the brightness distribution of the sky source are given by the unit vectors ( l , m , 1 l 2 m 2 ) . Visibility data produced by the correlator can thus be used for Fourier imaging after a series of calibration steps described in later sections.
In practice, the correlator determines Γ r 1 , r 2 ( u , v , n Δ τ ) over a discrete range of time “delay lags”, n Δ τ , or equivalently in the spectral domain, spectral channels N Δ ν across the baseband.

2.4. Correlator Delay Model

The correlator attempts to phase-align stations for maximum coherence to occur around τ = 0 . This requires a delay model that includes the station clock models (Section 2.2), and several predicted geometric and propagation-path-related delays. These stem from observing geometry, Earth orientation, tidal/loading effects on station coordinates, tropospheric and ionospheric propagation, near-field wavefront and relativistic aberration effects, telescope orbit and orientation parameters for space VLBI, etc. The delay model can be produced with the common NASA GSFC Calc/Solve package or software specific to space VLBI or near objects [65,66,67].
The a priori delay model cannot predict all encountered delay components. Visibility data correlated with a wide enough delay range n Δ τ and at sufficiently high temporal cadence allow unmodeled delay residuals, coming from random atmospheric turbulence, for example, to be calibrated out via fringe-fitting (Table 1) during post-processing (Section 3).

2.5. Correlator Parameters

Software correlators such as DiFX and SFXC are highly configurable and support various science and observing modes. Basic parameters are the PI-requested integration time and spectral resolution. Several additional data and settings flow into the correlator configuration, as illustrated in Figure 1.
In general, configuration files include a VEX file 5 derived from the observing schedule, edited as necessary to describe the actual rather than nominal frequency setups and baseband data layouts each station used during the observation.
Depending on correlator features, additional parameters may include, e.g., (1) pulsar ephemeris data and binning/gating setups, (2) additional interferometric phase centers for wide-field VLBI within the telescopes’ primary beam, and (3) settings for baseband spectral slicing or concatenation via DiFX “zoom” and “output” bands. These special band modes allow uniform-bandwidth visibility data records to be produced even for VLBI experiments where telescopes had incompatible frequency setups. In addition, (4) DiFX can also extract phase information of phase/pulse cal tone combs in the baseband data, if tones were injected during the observation. These combs are used for instrumental delay/phase calibration purposes, especially in geodetic VLBI. Furthermore, (5) at the VLBA, DiFX can extract on/off digital power levels of a synchronously switched calibration signal and produce continuous system temperature measurements (Section 4). Finally, (6) SFXC offers data windowing and weighted-overlap-add processing, which shape the frequency response of the spectral channels.
In networks, like the EVN and the VLBA, which provide frequent observations in standard modes and in a stable configuration, the preparatory steps and correlation can be largely automated and pipelined. In other instances and especially when there are technical changes at stations between sessions or there is no automated collection of metadata (station clock models, actual tunings, updated coordinates, etc), the step of preparing initial correlator VEX and configuration files can involve extensive manual editing. For VLBI arrays that do not provide regular observations with a consistent setup, multiple iterations of re-correlating plus data inspection on a small subset of scans are usually required until one arrives at a final correlator configuration.

2.6. Correlator Output

After a correlation run (bottom of Figure 1), the output data are typically converted from a correlator-native format into the widely adopted FITS-IDI 6, Common Astronomy Software Applications (CASA) Measurement Set (MS) 7, Haystack Observatory Postprocessing System (HOPS)/MkIV 8, or vgosDB NetCDF 9 Metadata or single-dish measurements (top of Figure 1) may further be attached to the visibility dataset as calibration tables. These may include system temperatures, antenna gain-elevation characteristics, weather data, flagging tables, and phase cal tones.
VLBI arrays typically provide observations in a circular polarization basis (through the use of quarter-waveplates), as the calibration of data from linear feeds is coupled to the different feed-angle rotations at different antennas in a complicated way. If observations had some telescopes observing in linear rather than circular polarization, the visibilities would have a mixed polarization basis, which complicates the analysis of the data and is not supported by the FITS-IDI format. The VLBI Global Observing System (VGOS), some EVN, and also (sub)mm-VLBI observations that include the phased Atacama Large Millimeter/submillimeter Array [ALMA, [3,68]] have to deal with this problem. Before being released to the PI, the mixed polarization visibilities are converted into a purely circular basis with the PolConvert [69] 10 program. The conversion requires telescope D-term and cross-polarization phase and gain characterization (Section 6.1. The ALMA Phasing Project quality-assurance (QA2) deliverables provide these parameters a priori for ALMA [70]. Alternatively, parameters can be derived from the visibilities themselves, via a PolConvert calibration pass on unpolarized sources.

3. Signal Stabilization

As described above, the correlator delay models are not perfect, which leads to residual delays and rates being present in the visibility data. Additionally, corrections for atmospheric phase fluctuations and instrumental bandpasses, where telescopes exhibit phase and amplitude variations over their frequency response, are usually not applied at the correlator. Post-correlation processing steps are therefore used to stabilize the source signal, ideally to remove all errors in the data down to the thermal noise. This allows us to average the data significantly in time and frequency, with limits set only by the structure and variability of the observed source in the presence of thermal noise only. In practice, the signal-to-noise ratio ( S / N ) of the data limits the degree to which data errors can be solved for. Without the signal stabilization steps, baseline-based errors will be frozen into the data due to phase decoherence, delay errors, and imperfect bandpass responses, for example.
Instrumental errors such as amplitude and phase bandpasses, phase and delay offsets between analog frequency bands, and polarization imperfections of the receivers (Section 6.1) are stable or drift only slowly over the course of observations. VLBI scans on bright calibrators can thus be used to solve for these effects. Some telescopes inject and log phase/pulse cal signals into their signal chain (Section 2.5). These signals appear as weak periodic spectral lines (and should therefore not be used for spectral line observations) and can be used to estimate instrumental delay errors.
A common fringe-fitting method is to pick a reference station, transform the visibilities on baselines to that station into the delay-rate space with a fast Fourier transform (FFT), and to pick delay and rate estimates at the peak location in the two-dimensional FFT space [71]. The delay, rate, and phase is set to zero for the reference station, and the FFT peaks are used as starting guesses for the station-based delays and rates of all other antennas. Additionally, the FFT S/Ns (Table 1) can serve as a threshold to distinguish between detections and non-detections. Finally, the data from all baselines on which the source is detected are used to refine the starting guesses and obtain final delay, rate, and phase solutions with a least-squares solver based on an assumed source model [72]. The default assumption of a point-source model does not usually lead to noticeable errors compared to the thermal noise [73]. Exceptions are highly sensitive arrays with a low thermal noise or geodesy experiments, where source structures might prevent reaching the required accuracy on the delays. If necessary, an iterative approach can be employed, where an image is first reconstructed (Section 5) using data calibrated with the assumption of a point source; that image is then used as a source model to fringe-fit the data again, and the data with this improved calibration are then used to create a final image. Alternatively, it is possible, although computationally very expensive, to solve for the combined fringes and source structure using a Bayesian approach [74].
Unlike instrumental errors, the residual post-correlation delays, rates, and atmospheric phase fluctuations are source/direction-dependent. In cm-VLBI experiments, atmospheric errors are moderate and calibration solutions can be transferred from nearby calibrators to the science target. The closer the calibrator is to the science target, the better the results are; ideally, both are within the same primary beams of the telescope. If the in-beam calibrator is weak, a “bi-gradient phase referencing” [75] can be employed, where a very bright source further awayis used to detect a closer weak calibrator via phase-referencing, from which phase-referencing is then used to detect the very weak science target. Two collinear aligned calibrators on opposite sides of a target source can be used for an improved correction of atmospheric errors [76]. With three or more calibrators, a two-dimensional phase screen of the atmosphere can be corrected for via two-dimensional interpolation of antenna phases using a “MultiView” technique [77,78]. At high observing frequencies, dominant path delay uncertainties are caused by the troposphere. Akin to geodesy VLBI, “geodetic blocks” of ∼ 10 sources can be observed in short succession to accurately determine the zenith tropospheric delays for each telescope [79]. Phase-referencing allows for astrometric measurements of the science target with respect to calibrator sources. At mm wavelengths, the science target itself has to be fringe-fitted and must therefore be bright enough to be detected within atmospheric coherence times, or incoherent averaging has to be used [80].
With multi-frequency receivers, it is possible to transfer phase solutions across bands to extend the coherence time at the higher frequencies and enable the detection of weaker sources. A simple frequency-phase-transfer (FPT) removes short-timescale tropospheric phase errors but does not allow for astrometric source registration. This technique can be combined with phase referencing between multiple targets as “source frequency phase referencing” [SFPR), [81,82,83,84]. SFPR works at mm wavelengths, removes slowly varying residual ionospheric and instrumental effects not captured by the FPT, enables astrometry, and allows for longer on-source VLBI scans as well as longer telescope slewing times to calibrator sources that are further away, since the fast tropospheric phase errors have been removed by the FPT. Rioja and Dodson [85] provide an outlook of high-precision astrometry measurements that will be enabled with future instruments and surveys.
The complex post-correlation data in multi-dimensional polarization (correlation products), frequency, time, and baseline space have a considerable size and are challenging to handle. Particular difficulties can arise in heterogeneous VLBI arrays, where different stations may require different treatments, results may be dominated by the measurements (and potential systematics) of the most sensitive telescope(s), and varying station sensitivities may make a uniform data reduction difficult. For example, calibration solutions with a high time resolution and/or frequency resolution can be obtained for the sensitive telescopes, while large solution windows might be required to accumulate enough S / N for the less sensitive telescopes, even when all baselines are used. For these considerations, array geometry also plays an important role; high S / N data from the large scale emission of resolved sources can be obtained for smaller dishes if they have short baselines to more sensitive telescopes.
Errors in the data, which can typically be traced back to single antennas, are often not trivial to find; a variety of diagnostic tools and visualizations that make different cuts through the data are needed. Some—but typically not all—errors can be identified at the signal stabilization step. Appendix B provides an overview of how bad data can be identified and dealt with at the different data-reduction stages. Calibration errors are easily baked into the data and can have substantial effects on the final scientific results. A typical order of VLBI calibration operations are shown in Figure 2. Comprehensive and well-tested software suites are generally used for the primary calibration steps. The most common software packages are described below. Note that tools for the flux density calibration (Section 4), imaging (Section 5), and advanced calibration methods for polarimetry (Section 6.1), spectral lines (Section 6.2), and wide fields of view (Section 6.3) are implemented in several of these packages.

3.1. AIPS

The Astronomical Image Processing System (AIPS) [86] 11 is currently the most widely used tool for the post-correlation calibration of VLBI data. The common visibility input data format that is produced by modern VLBI correlators is FITS-IDI. The calibrated output data are saved as a UVFITS 12 file to disk.
With AIPS, ionospheric corrections based on a priori maps of the global total electron content (TEC) can be applied to the data for an improved astrometry. Furthermore, the applied fringe corrections from the correlator can be improved with a posteriori determined IERS Earth Orientation Parameters (EOPs). Finally, pulse cal signals can be used to calibrate instrumental effects.
Several hundred tasks are available within AIPS, covering all standard data calibration and analysis steps as well as advanced methods, such as computing the statistics on Allan Variances (ALVPR task), modeling of gravitational lenses (GLENS), and optimizing configurations of telescope arrays (CONFI), to name a few examples.
The AIPS cookbook 13 can serve as a detailed guide to typical VLBI data-reduction pathways, where the needed AIPS tasks and their parameters are described for each step. Additionally, several generically useful VLBI tips and tricks are given.

ParselTongue

ParselTongue [87] 14 is a Python interface for AIPS. It allows the user to run AIPS tasks, access calibration and visibility data, and  script their data processing. A ParselTongue pipeline is currently used for EVN pipelining operations.

3.2. CASA

The CASA software [88,89] 15 and its underlying suite of C++ “casacore” libraries are widely used for the processing of connected-element-interferometer radio data (e.g., from ALMA and the VLA) and also single-dish radio observations. (e.g., [13,68,90,91]). Initiated by the “BlackHoleCam” project [92], CASA has been recently upgraded with VLBI capabilities, which are now supported and further developed by the Joint Institute for VLBI ERIC (JIVE) and NRAO [93]. The input visibility can be in the native, locally stored MS format or FITS-IDI files, which will be loaded and converted to an MS. The final calibrated visibilities are stored together with the raw data in the MS and can also be exported as UVFITS files.
The CASA calibration framework is based on the Hamaker–Bregman–Sault measurement equation, a powerful mathematical framework for interferometric calibration problems, including direction-dependent effects [94,95,96,97,98]. Unlike AIPS, CASA stores the visibility data (in MS format) and calibration tables (in an MS-like format) separately on disk in the working directories.
The fringe-fitter and all standard CASA tasks can make use of a source model to improve calibration accuracy. The source model information is stored in the MS when the imaging is performed within CASA (Section 5) or when an external FITS image is loaded. Ionospheric TEC corrections can also be performed in CASA. Additionally, the fringe-fitter has the capability of solving for dispersive delays from the ionosphere.
Development plans for the short-term future include the handling of pulse cal data, an improved handling of irregularly spaced frequency channels when performing multi-band fringe-fitting, EOPs correction, and a general accountability of the correlator model [93].
The next-generation CASA (ngCASA) project 16 will improve the hardware scalability significantly, particularly the speed with which large data volumes from future instruments can be accessed and processed, by using off-the-shelf Python libraries such as Dask [99] and xarray [100].

3.2.1. rPICARD

CASA has a convenient toolkit for extensive visbility and calibration data access, which allows users to easily automate their data processing. This led to the development of a generic VLBI calibration pipeline, the Radboud PIpeline for the Calibration of high Angular Resolution Data (rPICARD), which is based on CASA [73,101,102]. rPICARD employs flexible calibration strategies that work for any array and also synthetic VLBI data [103,104].
rPICARD uses MPI for scalability and a custom robust scalar bandpass calibration method. The ionospheric dispersive delay and Faraday rotation correction is performed with standard CASA tasks using TEC files from NASA’s Crustal Dynamics Data Information System. Fringe-fitting is performed with the CASA Schwab-Cotton algorithm [72], and a simple single-source phase referencing method is implemented in rPICARD. Instrumental phase, delay, and bandpass effects are corrected by utilizing all scans on calibrator sources. For corrections of atmospheric phase turbulence at high observing frequencies, the fringe-fitting is performed on a segmentation time, which is optimized based on the expected atmospheric coherence time at the observing frequency and the S / N of the fringe detections. Figure 3 shows this two-step fringe-fitting approach, where fringe detections are first obtained over VLBI scan durations and residual atmospherically-induced errors are subsequently calibrated out on short segmentation times.
Next to the aforementioned methods, which are implemented through standard CASA methods, rPICARD also has a few custom capabilities: (1) Exhaustive baseline searches to an arbitrary order are performed for non-detections to a primary reference antenna [Appendix D in [102]]. (2) The same rate solutions are applied to all polarization channels, and a simple nearest-neighbor algorithm is used to re-reference all fringe solutions to a single antenna for the stability of polarized source signals. (3) Robust solutions for atmospheric opacities are derived with estimations of the atmospheric temperature from the ATM software [105].
rPICARD can be obtained from a Bitbucket repository 17 and for every version of the pipeline, Docker and Singularity containers are available (details are given in the README file of the pipeline repository).

3.2.2. VPIPE

The VLBI PIPEline (VPIPE) pipeline is an end-to-end CASA-based VLBI pipeline designed for mid-frequency ( GHz ) observations 18. It currently supports the EVN, VLBA and LBA networks. The pipeline employs similar steps to the rPICARD pipeline (Section 3.2.1) including TEC corrections, phase referencing via fringe fitting, and self-calibration.
Originally designed for the high data volumes associated with wide-field VLBI data (typically 10–50 TB), the pipeline utilises MPI-CASA and HPC parallelism provided by cluster scheduling software, such as SLURM and PBS Pro, for enhanced scalability. This allows the pipeline to process data quickly and efficiently and permits the processing of “multiple phase centre observing” correlated data within short timescales (see Section 6.3) [54]. VPIPE is designed to be highly modular and supports containerization technologies, such as Singularity. This enables the user to easily add extra processing/analysis steps that use custom scripts or software packages.
The pipeline incorporates advanced direction-dependent calibration steps needed for wide-field VLBI data. These include multi-source self-calibration (MSSC; see Section 6.3), and primary beam correction schemes. This can be conducted in the ( u , v ) or image plane. Correcting for these allows the pipeline to reduce noise levels, improve the dynamic range, increase the effective field of view, and ensure consistent flux density scales across the primary beam.

3.2.3. The e-MERLIN CASA Pipeline

The e-MERLIN CASA Pipeline (eMCP) is a Python pipeline that is designed for use with e-MERLIN [106] data 19. It currently supports continuum and spectral line observations. The pipeline uses CASA-based VLBI tools, such as baseline and global fringe-fitting routines, to perform phase-referencing and self-calibration with parameters that are tailored for the e-MERLIN array and the observing frequency. The pipeline adopts other software tools such as AOFlagger [107] 20 and WSClean [108] 21 to address the challenging RFI environment and large image sizes often required.
Since 2019, eMCP has been incorporated as the official pipeline for the processing and delivery of e-MERLIN data to the users. As a result, it contains a suite of advanced quality assessment tools, displayed through a web interface, that allow the user to easily assess the quality of the calibration and images of the science targets.

3.2.4. JIVE EVN Continuum Jupyter Notebooks

EVN users are provided with a ParselTongue-based pipeline-calibrated version of their data by JIVE. Additionally, a CASA-based pipeline for EVN continuum is under development at JIVE and will run in a Jupyter notebook [93,109] 22. With this notebook, users will be able to start with a default set of parameters and steps used to calibrate EVN data, which can then be adjusted and refined to produce science-ready data.

3.3. EHT-HOPS

HOPS 23 is a collection of basic VLBI functions such as fringe-fitting that are sufficient for the current needs of geodetic analyses [110]. However, HOPS is not as flexible as AIPS and CASA and lacks a few capabilities, which prevents a full scientific data analysis. The missing features will be implemented in an updated next-generation HOPS (ngHOPS) software package.
In the meantime, an “EHT-HOPS” [73,111] pipeline has been designed, which augments HOPS with custom Python scripts 24 to enable the full reduction of data from the EHT and Global mm-VLBI Array (GMVA) when ALMA is participating in the observations [73,112]. A least-squares solver is used to find station-based delay and rate solutions from the baseline-based estimates that HOPS can solve for [71]. High S / N detections on bright calibrators are used to solve for bandpass effects. A piece-wise polynomial phase model is fitted directly to the visibilities to correct for atmospheric phase turbulence. To avoid over-fitting, phase corrections for each spectral window are obtained from a fit to all other spectral windows in the data. Relative complex gains between polarization channels are fitted to facilitate polarimetric analyses.
EHT-HOPS and rPICARD are used to produce science-ready data from the EHT [73,113,114,115].

3.4. PIMA

PIMA [116] 25 is a specialized data reduction program that can be used to accurately distinguish fringe detections from non-detections. PIMA is primarily employed for space VLBI experiments, because its fringe-fitter can solve for an acceleration (second-order rate) term. As space VLBI is beyond the scope of this work, we will not consider PIMA further.

4. Flux Density Calibration

After the signal stabilization, the correlation coefficients in units of thermal noise are normalized to unity autocorrelations and scaled corresponding to idealized analog correlation amplitudes. Ignoring baseline-dependent amplitude losses (such as de-correlation from data averaging), visibilities V i j can be obtained in flux density units of Jansky (Jy) by calibrating the complex correlation coefficients r i j of the baseline connecting stations i and j with the a priori estimated system equivalent flux density (SEFD) sensitivities of our antennas,
V i j = 1 η Q SEFD i SEFD j e τ i + τ j r i j .
Here, η Q is a correction factor for the digital quantization efficiency [44]; for 2-bit sampling, we have η Q 0.88 . To correct for atmospheric attenuation of ground-based millimeter and submillimeter observations and thereby measure V i j “above the atmosphere”, e τ correction factors are applied. Here, τ i is the mean atmospheric opacity along the line of sight signal path of antenna i.
The SEFDs are based on measured system temperatures T sys , the telescope gain G as a function of elevation, and the phasing efficiency η ph for phased arrays: SEFD η ph 1 G 1 T sys . Hot-load calibration scans measure the effective system temperature T sys e τ directly and thereby immediately satisfy Equation (3) [117]. For noise-diode-based T sys measurements, τ must be estimated after the fact. A common approach is to first approximate the system temperature as T sys T rx + 1 e τ T atm , where T rx and T atm are the two dominant noise terms, corresponding to receiver and atmospheric temperatures, respectively. Then, T rx is estimated by extrapolating T sys to zero airmass and an atmospheric model is used to estimate T atm from weather parameters measured by the telescope. The  T sys measurements themselves can subsequently be used to estimate τ (Section 4.3 in [102]). Examples of opacity-corrected system temperatures are shown in Figure 4.
The calibration metadata, gathered as a priori gain and system temperature information from all telescopes, are commonly stored in simple ASCII text files following the ANTAB format 26. Note that invalid entries are often marked with a value of 999 or 999 . Auxiliary weather information that might be needed to estimate T atm for an opacity correction is usually attached to the visibility data format. Otherwise, it can be loaded from an ASCII text file (by AIPS and rPICARD as described in [102]).

4.1. Software Implementations

AIPS, CASA/rPICARD, and EHT-HOPS post-processing scripts can be used to extract SEFDs from ANTAB tables and subsequently perform the a priori flux density calibration. AIPS reads external calibration information with the ANTAB task and creates SEFD calibration tables with APCAL, optionally with a simple fit for additional atmospheric opacity corrections 27. Here, the measured ambient temperature T amb is used to estimate the sky temperature via T atm = 1.12 T amb 50 K . CASA and rPICARD use custom Python modules 28 to load the ANTAB data if they are not already attached to the FITS-IDI files. The CASA task gencal is used to create SEFD calibration tables (separately for telescope gains and system temperatures). With a custom module, rPICARD can perform a robust atmospheric opacity correction for a user-specified selection of antennas. Here, the ATM code [105] is used to find T atm . The EHT Analysis Toolkit 29 provides post-processing scripts for the flux density calibration of EHT-HOPS data.

4.2. Advanced Methods Based on Array Redundancy and a Priori Source Assumptions

  • A network calibration [73,111] can be employed if the unresolved flux density on large scales seen by short baselines in the VLBI array is known. The method only works for telescope sites that are close enough (almost co-located or in “walking distance”) to effectively form a zero-baseline interferometer. The flux density measurement used to calibrate the gains can be obtained from (quasi-)simultaneous single-dish or connected-element-interferometry observations and allows for an absolute amplitude calibration. A least-squares approach with all baselines to pairs of redundant sites are used to robustly constrain the gains of the two neighboring antennas. A convenient implementation of the network calibration method can be found in the eht-imaging software (Section 5.1.5).
  • Cross-track calibration [118] is based on the fact that redundant baselines anywhere in the ( u , v ) space should measure the same source properties’ modulo intrinsic source variability. A threshold can be set for how close two baselines should be to be considered identical, given how “quickly” the source structure varies in the Fourier space. Baselines might cross within a small region or stay closely parallel for long ( u , v ) -tracks. Network calibration is a special case of the cross-track calibration, where a known total flux density can be used for an absolute gain calibration. Generally, we have no a priori knowledge about the resolved source structure at larger ( u , v ) -spacings. Here, the least-squares solver can be used for all redundant baselines to tighten the gains of the involved stations by constraining the amplitude ratio of crossing-track baselines to unity. In the absence of accurate and independent flux density information, the product of the solved gains are enforced to be unity to only solve for relative gains without adjusting the total flux density. Within AIPS, the UVCRS task can be used as a diagnostic tool to identify anomalous gains of stations that are involved in crossing ( u , v ) -tracks. The only real cross-track calibration method implementation that the authors are aware of is the UVCROSS Caltech VLBI Analysis Program [119] 30.
  • Second-moment source size calibration [120] can be employed under the assumption that the short baselines in the array sample a simple large-scale source structure such as a Gaussian. From the assumed large-scale image, model amplitudes can be computed and gains from the stations connected by short baselines can be obtained by self-calibration (Section 5). For a single baseline, a common application is to keep the gains fixed for the station that has the more accurate SEFD-based a priori flux density calibration. This method can be employed by all imaging software packages’ self-calibration routines. A convenient implementation can be found in eht-imaging.

5. Imaging and Geometric Model-Fitting

When imaging VLBI data, one tries to reconstruct the sky brightness distribution that fits the data best and satisfies additional assumptions imposed to deal with the ill-posed imaging problem of incompletely sampled Fourier data from a sparse array. Strictly speaking, “imaging” refers to obtaining a pixel-based model (convolved with a resolving beam) through inverse modeling or forward modeling [121]. Inverse modeling techniques such as the CLEAN algorithm [122,123] gradually reconstruct a model in the image domain of the inverse Fourier transform of the visibilities. There are a couple of different CLEAN implementations in different software packages that are worth mentioning. The original Högbom CLEAN method [122] uses iterations of delta image model components convolved with the point spread function (PSF) that are subtracted from the “dirty” image of the inverse Fourier transform. The process converges when the residual map consists only of noise. The collection of delta components form the final image model. When the residual image is added and the delta components are convolved with the VLBI resolving beam, the final image is obtained. Knowledge about interferometric images is needed to identify “false” sidelobe emission regions and to place boxes (“CLEAN windows”), which limit the locations in which CLEAN can put delta components: the emission peaks within the CLEAN windows. Clark [123] invented a faster variant of the Högbom CLEAN for large images. Clark CLEAN operates in two cycles; the minor cycle works like the Högbom CLEAN, except that only the inner “beam patch” of the PSF is used. In major cycles, the collection of delta components are then subtracted with the full PSF. Multi-scale CLEAN [MSC, [124]] allows the use of extended Gaussian components for the model next to the delta components. This proved to be useful in accurately deconvolving extended/diffuse emission.
These techniques are used to average the frequency bandwidth into a single channel to improve S / N . However, the development of multi-frequency synthesis [MFS, [125,126,127]] allowed the gridding of visibilities measured across multiple frequencies onto the same ( u , v ) plane, thus improving the fidelity of the resultant image. With the widening bandwidths of modern instruments, the in-band source spectral indices limited the deconvolution accuracy when using MFS. To solve this, the multi-term multi-frequency synthesis (MT-MFS, [128]) algorithm models the frequency dependence of the source using a Taylor expansion. These various algorithms can be used in combination with each other depending on the observational configuration and source structures. For example, the Multi-scale multi-frequency synthesis [MS-MFS, [128]] combines MSC with MT-MFS.
Forward modeling techniques such as the maximum entropy method [129] fit the Fourier transform of an image to the visibility data. Additional assumptions can be incorporated as regularizers. The fitting of geometric forms such as two-dimensional Gaussian components is also performed with the forward modeling approach in the visibility domain. Figure 5 shows an example VLBI image reconstruction together with the underlying visibility amplitude measurements. Projected along the jet direction, the visibility amplitudes are indicative of a smooth Gaussian-like structure. Perpendicular to the jet direction, bounding amplitudes are indicative of strong intensity gradients across the transverse jet profile.
While the VLBI scientific analysis is usually performed in the image domain, it is also possible to compare physics-based theoretical source models directly to the visibility measurements. The direct data-model comparison method will likely become more prominent in the future with the availability of increasingly advanced AGN jet models (e.g., see [130] for a recent review) and machine learning methods [131,132]. For the analysis of M87* and Sgr A* EHT observations, horizon-scale simulations [121] are directly compared to the visibility data with (proprietary) software packages such as GENA [133] and THEMIS [134].
Specific combinations of the baseline-based visibilities can be formed as “closure quantities”, where station-based data errors drop out [135,136,137,138]. For example, the sum of baseline phases in a closed triangle form a robust closure phase. While uncorrected residual gain errors negatively impact reconstructed images when the full set of visibilities is utilized, with only closure quantities, one does not make use of the full information content of the data. All imaging software packages listed below (Section 5.1) are able to perform self-calibration, where, for sufficiently bright sources, within given solution intervals, per-station phase and/or amplitude gains can be found with a least-squares solver based on a comparison between the observed data and model visibilities. Typically, an initial self-calibration is performed based on a simple starting model or an initial image reconstruction where only closure quantities have been used. Subsequently, one can iteratively image and self-calibrate until the gain solutions are converged. This method of solving for the observed source structure together with calibration gains is referred to as “hybrid-mapping” [139,140].
All imagers listed below can make use of the UVFITS visibility data format as input and produce FITS images as output. Full Stokes polarimetric imaging can be performed with each imager. For details, the reader is referred to ([141], and Martí-Vidal et al., in prep.).

5.1. Software Implementations

5.1.1. AIPS

The AIPS IMAGR task implements the Clark CLEAN method.

5.1.2. CASA

The CASA tclean task implements the Högbom, Clark, MSC, and MS/MT-MFS deconvolver options for CLEAN. Next to UVFITS input, the CASA imager works with its native MS visibility and metadata file format.
CASA is commonly used for spectral-line and wide-field observations, for which the corresponding special tclean features are described in Section 6.2 and Section 6.3, respectively.

5.1.3. Difmap

Difmap [142] 31 is a widely used VLBI imaging and geometric model fitting software owing to its computational efficiency, intuitive user interface, and convenient data plotting plus flagging procedures. Difmap implements Högbom CLEAN.

5.1.4. WSClean

WSClean [108] 32 has been developed with the primary goal of computationally efficient wide-field imaging, for which the corresponding special software features are descried in Section 6.3. The latest version of the software can use the Clark CLEAN implementation and reconstruct MS-MFS images computationally more efficient than tclean [143].
WSClean uses the MS file format. UVFITS input can easily be converted with the CASA importuvfits task.

5.1.5. eht-Imaging

The eht-imaging software package [144,145] 33 is a collection of Python modules for VLBI data analysis, with a focus on regularized maximum likelihood (RML) imaging. It includes routines for detailed mock data simulation, gain and polarimetric calibration of visibilities, image and data plotting and statistical analysis, and RML image reconstruction, as well as geometric model fitting from visibilities and closure products (including full-Stokes and multi-frequency data).

5.1.6. SMILI

SMILI [146,147] 34 is a collection of sparse sampling RML forward modeling libraries, similar to eht-imaging. SMILI can be used through a Python interface.

5.1.7. UVMULTIFIT

UVMULTIFIT [148] 35 is a flexible geometric model fitting library that is integrated in the CASA ecosystem.

5.1.8. Comrade and DPI

The Julia-based [149] Comrade [150] and Python-based Deep Probabilistic Imaging/Inference (DPI/ α -DPI, [151,152]) frameworks are versatile geometric model fitting tools, which are used by the EHT [153].

6. Advanced Scientific Applications

6.1. Polarization Calibration

For polarization considerations, we consider the sky signal to be split into independent right-handed and left-handed circular polarizations (RCP and LCP, respectively). Such a signal split is commonly obtained in VLBI receiving systems by placing quarter-waveplates in front of linear polarization feeds. Advanced calibration methods, which have traditionally been developed for the inclusion of ALMA in VLBI experiments [69], enable routine conversions of recorded linear polarization signals to circular ones during correlation (Section 2).

6.1.1. Additional Polarization Signal Stabilization Steps

The first step is to properly align the RCP and LCP signals. The geometric feed rotation angle evolution between the RCP and LCP data (commonly referred to as parallactic angle rotation) depends on a telescope’s focus and mount configuration (e.g., see Appendix C of [102]) and is usually taken out when the data are fringe-fitted (Section 3). 36 When all RCP and LCP fringe solutions are (re-)referenced to a common station in a stable system, a global cross-hand R-L delay will be present in the data. This is the instrumental delay between the RCP and LCP signal chains of the chosen reference station, which can be determined by fringe-fitting the RL and LR visibilities. The R-L phase is similarly affected by instrumental effects of the reference station. The true absolute R-L phase sets the orientation of the electric vector polarization angle (EVPA) of the polarized radio emission, which has to be bootstrapped from other (e.g., single-dish) observations.
If the fringe-fitting is conducted separately for RCP and LCP, the R-L alignment can be performed as final calibration step. If polarizations are combined for the fringe-fitting, the alignment should be performed as part of the instrumental calibration steps (Figure 2).

6.1.2. Solving for Polarization Leakage Effects

Subsequently, the mutual leakage of signals between the RCP and LCP signal paths at every station should be calibrated. If left uncorrected, these leakages or “D-terms” will impose small errors in total intensity data and add polarization signals that are purely instrumental (i.e., Stokes Q and U signals can occur for unpolarized sources). D-terms can have a frequency dependence, are usually stable in time, and are small enough that the signal stabilization calibration does not have to be revised with leakage-corrected visibilities. Before leakages can be determined, a good image (Section 5) of the observed source is required. This is, firstly, to have gain errors removed by self-calibration and, secondly, because solving for D-terms for resolved sources with complex polarization structures is not trivial [141]. For a thorough and recent overview of D-term calibration methods and software implementations, the reader is referred to [[141], and Martí-Vidal et al., in prep.].

6.1.3. Circular Polarization

For Stokes V studies of VLBI data in a circular polarization basis, self-calibration solutions obtained from the RR and LL visibilities are needed to estimate the relative R/L amplitude gains of every antenna in the array. The problem is that an intrinsic Stokes V source polarization will confound the true residual R/L gains. A common remedy is to plan an observation where R/L gains from many observed sources based on Stokes I images can be obtained. Under the assumption that V = 0 on average for the sample, as there is no intrinsic preference for positive or negative intrinsic circular polarization, the mean R/L gains should be applied to the data [154].

6.2. Spectral Line Observations

For telescopes equipped with spectrometers and a sufficient frequency resolution set at the correlator, emission and absorption features from spectral lines can be studied.

6.2.1. Spectral Line Signal Stabilization

For spectral line observations, a few special considerations and adjustments of certain steps have to be taken into account for the signal stabilization (e.g., [155]). Fringe-fit delay solutions and bandpass corrections should be obtained only on continuum sources. Depending on the reference frame used during correlation, Doppler shifts of spectral lines have to be corrected to keep the position of the line constant with respect to the correlated frequency channels after delays and the bandpass have been corrected. The center of the Earth is usually chosen as reference point during correlation, and small corrections that are common for all antennas due to the Earth’s orbit around the Sun relative to the observed source are applied when shifting to the Local Standard of Rest. Such effects become increasingly relevant for longer observations. Fringe rates can be determined on the science target, selecting one or a few channels with strong line emission. Strong line emission can also be used for phase-referencing.
The AIPS and CASA fringe-fitting tasks are able to handle spectral line observations. The Doppler shifts can be corrected with the combined SETJY and CVEL tasks in AIPS and the mstransform task in CASA. The CASA imager is able to reconstruct spectral cubes (Section 5).

6.2.2. Template Spectrum Flux Density Calibration

Bright lines will show up in telescope’s auto-correlation spectra, exceeding the atmospheric noise. A template spectrum can be obtained from the measured bandpass-corrected total power response of the line from a sensitive antenna in the VLBI array. Amplitude gains in the form of time-dependent SEFDs can then be obtained without using system temperatures (Section 4) by fitting the total power spectral of all antennas to the template ([156], Lecture 12).
This calibration method is implemented in the AIPS task ACFIT.

6.2.3. Spectral Line Imaging

For spectral lines, the self-calibration is best performed with the image of the peak line emission channel(s). Subsequently, spectral cube images can be made from the channels/velocities where line emission is present. If necessary, CASA tclean can perform Doppler tracking to keep the positions of spectral lines stable in the Kinematic Local Standard of Rest. For an example of a particularly extensive spectral line VLBI imaging work, the reader is referred to Matthews et al. [157].

6.3. Wide-Field VLBI

In wide-field VLBI experiments, multiple sources are observed within single pointing fields that cover the full or a significant fraction of the participating telescopes’ primary beams 37.

6.3.1. Wide-Field Correlation

For us to observe sources across the entire primary beam using VLBI, we require correlation to be conducted at very high time (typically millisecond) and frequency (∼kHz) resolutions in order to restrain smearing at the edge of the primary beam. However, the result is that a single large and often unwieldy dataset is produced, which makes calibration, imaging, and subsequent analysis difficult without significant computing resources. This method has been used for many of the early wide-field VLBI experiments (e.g., [159]), but the computational complexity and correlation resources to map the entire primary beam made wide-field VLBI observations an unattractive concept.
An alternative correlation method for wide-field VLBI was developed with the DiFX software correlator [160] and has been implemented in the JIVE SFXC correlator [55]. Named “multiple simultaneous phase centre observing”, this uses a two-step correlation approach. Firstly, an internal wide-field correlation produces data at the required time and frequency resolution to restrain smearing to acceptable levels at the edge of the primary beam. These data are copied and phase-rotated to multiple positions across the primary beam FOV, which retains high-precision astrometry as an internal geometric delay model within the correlator applies the required phase rotation. Each of these datasets, which are centered at different positions, are then averaged to coarse time and frequency resolution, producing multiple pencil-beam datasets across the primary beam. These pencil beams can be arranged to cover multiple sources discovered in accompanying low-resolution surveys or be mosaicked to cover a region of interest or even the whole primary beam FOV. The result from this mode of correlation is multiple small and manageable datasets (often ∼GB-sized), one for each position.

6.3.2. Wide-Field Signal Stabilization

Signal stabilization steps (Section 3) can be applied to wide-field observations as usual, by applying the calibration solutions to the data from all phase centers. Calibrator sources may be found within the phase centers as in-beam calibrators or from another area of the sky covered in a different pointing direction.

6.3.3. Wide-Field Imaging: Primary Beam Correction, the W Term, and Multi-Source Phase Self Calibration

When the beams of the participating telescopes are well determined, primary beam corrections can be employed across the FOV. For homogeneous arrays, such as the VLBA, the primary beam effectively adjusts the flux densities with distance from the pointing centre as each baseline sees the same apparent sky brightness distribution that is modified by the primary beam power envelope 38.
However, most VLBI arrays contain heterogeneous elements, meaning that the primary beam has a significant effect even within the main lobe of the primary beam. This manifests in the data as a direction-dependent and antenna-independent gain error (see Radcliffe et. al. in prep). Corrections for this effect must be conducted in the ( u , v ) -plane and can be achieved in two ways. The first is to produce a direction-independent gain table with the corrections corresponding to the center phase position (implemented through the task CLVLB in AIPS and internally within the VPIPE pipeline). This is computationally inexpensive, can account for the frequency and time dependence of the primary beam, and can be applied together with the phase-referencing corrections.
The second method is to correct during the gridding of the image reconstruction using the convolution of the complex conjugates of the primary beam voltage responses. This is more computationally complex but any residual errors are only dependent upon the accuracy of the primary beam models. Prototypes have been implemented in WSClean through the image-domain gridder (IDG) [161] through the application of a direction-dependent diagonal gain correction (seeRadcliffe et al. in prep).
This correction is applied together with corrections for a non-negligible w term [44] that causes increasingly larger errors in images with distance from the phase center due to sky curvature and non-coplanar baselines. Due to the w term, the Fourier-transform between the visibilities and the sky brightness distribution is no longer two-dimensional. A full three-dimensional FFT would result in strong aliasing along the w direction, and a direct transform along w is computationally too expensive for a cube that is mostly devoid of emission.
The w term effects are typically corrected with one of the following methods:
  • Faceting [162] performs a 2D Fourier-transform for a number of different phase centers across the FOV. The resultant pieces of image facets are small enough that the 2D FFT approximation works. The facets are stitched together and deconvolved.
  • w-projection [163] uses the fact that the visibility V as a function of w can be calculated from V at w = 0 through a convolution operation: V ( u , v , w ) = G ˜ ( u , v , w ) V ( u , v , w = 0 ) . Here, G ˜ ( u , v , w ) is the Fourier transform of G ( u , v , w ) = exp 2 π i w 1 l 2 m 2 1 . This method is faster than faceting because the visibilities need to be gridded only once.
  • w-stacking [108] uses a w-dependent grid that is FFT’ed for each gridded w value and phase-shifted by G 1 ( u , v , w ) . For imaging, all grids are summed with appropriate scaling factors (see [108], for details). Compared to w-projection, w-stacking is faster when the visibility gridding is computationally more expensive than the FFT computations.
WSClean implements the w-stacking method and tclean implements the w-projection algorithm and can perform faceting.
The phase stability of VLBI arrays is worse than connected short-baseline arrays due to the differing atmospheric paths between antenna and source. This means that self-calibration (Section 5) is essential in achieving low-noise and high-dynamic ranges. Self-calibration requires a source with sufficient S/N to estimate the gain corrections along the line of sight to the target field. However, for VLBI observations, its lack of sensitivity to large-scale, diffuse emission means that the total flux available for self-calibration of a single source can be insufficient, restricting VLBI observations to those fields with in-beam, or nearby phase calibrators.
However, for wide-field VLBI, the targeting of multiple sources allows self-calibration to proceed, even when each individual target is faint, by using the combined flux density of the detected sources [164,165]. This technique, called “multi-source self-calibration”, images each source detected after phase referencing and divides the visibilities by the model to produce a point source at the phase center. These datasets are stacked in the ( u , v ) -plane and self-calibration solutions obtained using the combined S/N of all target sources. This technique is currently being developed into a fully directional-dependent technique through the derivation of self-calibration solutions along different line-of-sights across the primary beam (Harth et al. in prep).

7. Synthetic Data

Forward-modeling methods, where the ( u , v ) -coverage, thermal noise, and simple telescope gain errors are taken into account to create simulated/fake visibility data based on a model sky brightness distribution, have been available for a long time (e.g., the DTSIM AIPS task, simobserve in CASA, and MeqTrees [166] 39). Recently, very accurate visibility simulations have been developed, namely the data-simulation toolkit within eht-imaging and the SYnthetic Measurement creator for long Baseline Arrays (SYMBA) [167] 40. Both methods make use of sophisticated data-corruption models, are able to include the effects of interstellar scattering, can use arbitrary-source models defined in text files or images from HDF5 and FITS data, and use the exact ( u , v ) -coverage from (observational) UVFITS input. SYMBA and eht-imaging are routinely used by the EHT to generate synthetic data for dedicated tests to answer specific scientific questions (e.g., about the robustness of specific features in reconstructed images), optimize calibration procedures, produce imaging parameter surveys, create test suites for software packages, and identify how to optimally upgrade VLBI arrays (e.g., [10,114,121,168,169,170]).
The versatile eht-imaging package bootstraps data-corruption effects from a calibrated UVFITS file of observational data and produces a UVFITS output. The CASA-based SYMBA software produces both UVFITS and MS output. Data-corruption effects are simulated from first principles with MeqSilhouette to emulate the correlator output in a real observation [103,104] 41. For example, physical parameters for atmospheric or antenna-pointing models can be adjusted and their effects studied. The suite of possible MeqSilhouette data-corruption effects is continuously being expanded; Faraday rotation effects from the interstellar medium and atmospheric effects from the ionosphere will soon be added, for example. SYMBA passes the MeqSilhouette output to rPICARD, which results in science-ready calibrated synthetic data akin to those from a real observation. Work is ongoing to convert SYMBA to simulate low-frequency VLBI and enable feasability studies of new instruments such as SKA-VLBI.
The ability to simulate VLBI data based on known ground-truth models allows non-experts to better understand the fundamentals of VLBI data, for example, by looking at the amplitude and phase signatures of simple geometric models or by studying how specific calibration errors manifest themselves in the data. Examples of synthetic EHT observations based on state-of-the-art theoretical models are shown in Figure 6.
In addition to the full forward-modeling software packages, there are also dedicated VLBI observing simulators, such as the Very long baseline interferometry Network SIMulator (VNSIM) [172] 42. These tools assist users and engineers with the design and scheduling of VLBI experiments, as well as the evaluation of the array performance under various conditions.

8. Summary

In this work, we have presented current open-source VLBI methods and their implementations in various software. Table 2 provides and overview of the major software packages’ capabilities. Several of these packages can be easily obtained through the KERN repository [173] 43 alongside many additional radio astronomical software libraries.

Author Contributions

Conceptualization, M.J.; methodology, M.J., J.F.R. and J.W.; software, M.J., J.F.R. and J.W.; writing—original draft preparation, M.J., J.F.R. and J.W.; writing—review and editing, M.J., J.F.R. and J.W.; visualization, M.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Research Council (ERC) Advanced Grant “M2FINDERS—Mapping Magnetic Fields with INterferometry Down to Event hoRizon Scales” (Grant Agreement No. 101018682).

Acknowledgments

The authors thank Ciriaco Goddi, Freek Roelofs, and Ilse van Bemmel for their valuable comments to the manuscript. We also thank the anonymous referees for helpful comments that improved the paper.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
AGNActive galactic nuclei
ALMAAtacama Large Millimeter/submillimeter Array
ASCIIAmerican Standard Code for Information Interchange
ASICApplication-specific integrated circuit
CLEANImaging algorithm for an incomplete Fourier coverage
EHTEvent Horizon Telescope
EVNEuropean VLBI Network
EVPAElectric vector polarization angle
EOPEarth orientation parameter
FFTFast Fourier transform
FPGAField-programmable gate array
FOVField of view
FPTFrequency-phase-transfer
GbpsGigabit per second
GMVAGlobal mm-VLBI array
GPSGlobal Positioning System
GPUGraphics processing unit
GRMHDGeneral relativistic magnetohydrodynamics
HPCHigh-performance computing
IDGImage-domain gridder
IERSInternational Earth Rotation and Reference Systems Service
JIVEJoint Institute for Very Long Baseline Interferometry European Research Infrastructure Consortium
KVNKorean VLBI Network
LCPLeft circularly polarized
LOFARLow Frequency Array
MADMagnetically arrested accretion disc
MERLINMulti-Element Radio Linked Interferometer Network
MFSMulti-frequency synthesis
MPIMessage Passing Interface
MSCMulti-scale CLEAN
MSSCMulti-source self-calibration
MS-MFSMulti-scale multi-frequency synthesis (combines MSC with MT-MFS)
MT-MFSMulti-term multi-frequency synthesis
ngEHTNext-generation Event Horizon Telescope
ngVLANext-generation Very Large Array
NRAONational Radio Astronomy Observatory
RCPRight circularly polarized
RMLRegularized maximum likelihood
PIPrincipal investigator
PSFPoint spread function
RFRadio frequency
RFIRadio frequency interference
SANEStandard and normal evolution accretion state
SEFDSystem equivalent flux density
SFPRSource frequency phase referencing
SKASquare Kilometre Array
S / N Signal-to-noise ratio
TCPTransmission Control Protocol
TECTotal electron content
UDPUser Datagram Protocol
UTUniversal Time
VERAVLBI Exploration of Radio Astrometry
VGOSVLBI Global Observing System
VLAVery Large Array
VLBAVery-Long Baseline-Array
VLBIVery-long-baseline interferometry

Appendix A. Telescope Baseband Data Transport to the VLBI Correlator

Baseband data from all stations need to be available simultaneously for correlation. One option is real-time streaming to the VLBI correlator. This is an observing mode in e-MERLIN, e-EVN, and future VLBA after ongoing upgrades.
Real-time correlation is possible if the network paths from the VLBI telescopes (connected at typically 1 G/10 G/100 G) to the correlator (typically 10 G/100 G) offer sufficient capacity and can maintain real-time rate reliably with an acceptably low data loss fraction. On paths that cross the Internet, this can be challenging. A drawback is also that re-correlations with different settings are not possible, unless the baseband data are stored.
In most VLBI experiments, the baseband data are captured by the telescopes onto a recorder such as FlexBuff 44, Mark 6 45, or Octadisk2 46. The recording media can be physically shipped by courier mail. The proprietary form factor of the media requires a corresponding playback unit at the correlator.
Alternatively, baseband data files can be transferred at a slower pace over the Internet to the correlator. Popular high-speed data transfer software for this purpose includes JIVE jive5ab 47 and etc/etd 48, Tsunami UDP transfer 49, and Globus GridFTP 50. Transfer with slower scp or rsync is also used as a fallback. In an “e-shipping” mode, the transferred data are stored on a cluster or cloud file system until later correlation. In contrast, real-time correlation of the transferred data is commonly referred to as “e-VLBI”.

Appendix B. Data Flagging

The flagging (removal) of “bad data” is common practice in interferometry. With enough antennas in the array, one can be rather liberal with the data flagging. VLBI arrays, however, are often sparse. Care must be taken to distinguish bad data from correctable calibration errors and measurements with low S / N . Common uncorrectable problems are issues throughout the telescope data recoding path that may render data at specific frequency channels or polarizations unusable, telescopes that are late at their source or experience severe weather, bandpass fall-offs (“frequency edge-channels”), and radio frequency interference (RFI) 51. In the unaveraged data, baseline-based problems are rare, and the flagging should be performed on a per-antenna basis. There are some algorithms that can automatically identify some bad data, e.g., for RFI [107], but to the authors’ knowledge, there is no software that is fully sufficient for VLBI data. The VLBI data flagging therefore remains a subjective task.
Flags are usually applied at different data-reduction stages. A priori flags based on telescope operator logs are made available to PIs. Low S / N -fringe non-detections render the underlying data flagged, but care must be taken here as the fringe-fitting segmentation time does not necessarily match the time ranges over which data issues are present. Partially affected data over VLBI scan durations should be removed before fringe-fitting. Finally, manifestations of bad data and calibration errors can easily be identified in image reconstructions as spurious features. Figure A1 shows two typical examples of bad data that should be flagged.
Figure A1. Examples of flagged bad data indicated with crosses. These plots are produced automatically with jiveplot (Available online: https://github.com/haavee/jiveplot, accessed on 30 July 2022) by rPICARD and show visibility amplitudes and phases per baseline and per scan, color-coded by polarization. Top: Scan-averaged EVN data on the Effelsberg–Westerbork baseline. The Effelsberg RCP channel and first four spectral windows contain no signal and have been flagged. Bottom: Frequency-averaged VLBA data on the Brewster–Hancock baseline. Brewster has been a few seconds late on source, resulting in bad phases and amplitudes at the start of the scan.
Figure A1. Examples of flagged bad data indicated with crosses. These plots are produced automatically with jiveplot (Available online: https://github.com/haavee/jiveplot, accessed on 30 July 2022) by rPICARD and show visibility amplitudes and phases per baseline and per scan, color-coded by polarization. Top: Scan-averaged EVN data on the Effelsberg–Westerbork baseline. The Effelsberg RCP channel and first four spectral windows contain no signal and have been flagged. Bottom: Frequency-averaged VLBA data on the Brewster–Hancock baseline. Brewster has been a few seconds late on source, resulting in bad phases and amplitudes at the start of the scan.
Universe 08 00527 g0a1

Notes

1
The sched user manual (Available online: http://www.aoc.nrao.edu/software/sched (accessed on 30 July 2022)). See also pySCHED ( Available online: https://github.com/jive-vlbi/sched (accessed on 30 July 2022)) and the EVN Observation Planner (Available online: https://planobs.jive.eu (accessed on 30 July 2022)).
2
DiFX. Available online: https://www.atnf.csiro.au/vlbi/dokuwiki/doku.php/difx/documentation, accessed on 30 July 2022.
3
SFXC. Available online: https://www.jive.eu/jivewiki/doku.php?id=sfxc, accessed on 30 July 2022).
4
GPU Acceleration of the DiFX Software Correlator. Available online: https://adacs.org.au/project/gpu-acceleration-of-the-difx-software-correlator, accessed on 30 July 2022.
5
VEX definition. Available online: https://vlbi.org/wp-content/uploads/2019/03/vex-definition-15b1.pdf, accessed on 30 July 2022.
6
FITS-IDI definition. Available online: https://fits.gsfc.nasa.gov/registry/fitsidi/AIPSMEM114.PDF, accessed on 30 July 2022.
7
MS definition version 2.0. Available online: https://casa.nrao.edu/Memos/229.html, accessed on 30 July 2022.
8
Mk4 File Format Definitions. Available online: https://www.haystack.mit.edu/wp-content/uploads/2020/07/docs_hops_002_mk4_files.txt, accessed on 30 July 2022.
9
vgosDB NetCDF. Available online: https://vievswiki.geo.tuwien.ac.at/VLBI-Analysis/Input-data, accessed on 30 July 2022.) visibility data formats.
10
PolConvert. Available online: https://github.com/marti-vidal-i/PolConvert, accessed on 30 July 2022.
11
Astronomical Image Processing System. Available online: http://www.aips.nrao.edu, accessed on 30 July 2022.
12
The UVFITS format is following AIPS conventions. AIPS Memo 117 revised. Available online: ftp://ftp.aoc.nrao.edu/pub/software/aips/TEXT/PUBL/AIPSMEM117.PS, accessed on 30 July 2022. The AIPS FITS format can sometimes change. The format is therefore not well defined and different software packages use slightly different conventions
13
AIPS cookbook. Available online: http://www.aips.nrao.edu/cook.html, accessed on 30 July 2022.
14
The Parseltongue Wiki: https://www.jive.eu/jivewiki/doku.php?id=parseltongue:parseltongue, accessed on 30 July 2022.
15
CASA. Available online: https://casa.nrao.edu, accessed on 30 July 2022.
16
ngCASA. Available online: https://pypi.org/project/ngcasa, accessed on 30 July 2022.
17
rPICARD pipeline. Available online: https://bitbucket.org/M_Janssen/picard, accessed on 30 July 2022.
18
Generic VLBI pipeline. Available online: https://github.com/jradcliffe5/VLBI_pipeline, accessed on 30 July 2022.
19
eMERLIN CASA pipeline. Available online: https://github.com/e-merlin/eMERLIN_CASA_pipeline, accessed on 30 July 2022.
20
AOFlagger. Available online: https://gitlab.com/aroffringa/aoflagger, accessed on 30 July 2022.
21
WSCLEAN. Available online: https://gitlab.com/aroffringa/wsclean, accessed on 30 July 2022.
22
EVN_CASA_pipeline. Available online: https://code.jive.eu/bemmel/EVN_CASA_pipeline, accessed on 30 July 2022. See also the EVN data reduction guide. Available online: https://www.evlbi.org/evn-data-reduction-guide, accessed on 30 July 2022.
23
Haystack Observatory Postprocessing System. Available online: https://www.haystack.mit.edu/haystack-observatory-postprocessing-system-hops, accessed on 30 July 2022.
24
The Event Horizon Telescope Analysis Toolkit. Available online: https://github.com/sao-eht/eat, accessed on 30 July 2022.
25
VLBI processing software PIMA. Available online: http://astrogeo.org/pima, accessed on 30 July 2022.
26
ANTAB format. Available online: http://www.aips.nrao.edu/cgi-bin/ZXHLP2.PL?ANTAB, accessed on 30 July 2022.
27
Opacity Correction for High Frequency VLBI Observations. Available online: https://library.nrao.edu/public/memos/vlba/sci/VLBAS_01.pdf, accessed on 30 July 2022.
28
CASA-VLBI Scripts. Available online: https://github.com/jive-vlbi/casa-vlbi, accessed on 30 July 2022.
29
https://github.com/sao-eht/eat, accessed on 30 July 2022.
30
The Caltech VLBI Programs. Available online:https://sites.astro.caltech.edu/~tjp/citvlb, accessed on 30 July 2022, have been discontinued and are therefore not considered further in this work
31
Difmap. Available online: ftp://ftp.astro.caltech.edu/pub/difmap/difmap.html, accessed on 30 July 2022.
32
WSCLEAN. Available online: https://gitlab.com/aroffringa/wsclean, accessed on 30 July 2022.
33
ehtim. Available online: https://github.com/achael/eht-imaging, accessed on 30 July 2022.
34
SMILI. Available online: https://github.com/astrosmili/smili, accessed on 30 July 2022.
35
UVMultiFit’s Documentation. Available online: http://mural.uv.es/imarvi/docums/uvmultifit, accessed on 30 July 2022.
36
Note that the telescope focus and mount configuration information are not always correctly specified in the visibility data files produced by the correlators. For all software packages described here, there are methods to overwrite this information manually.
37
See Strom [158] for the primary beam of a heterogeneous interferometer array.
38
We note that this is not strictly the case outside the main primary beam lobe, especially for non-equatorial mounts, where the different parallactic angles at each antenna result in different side-lobe structure
39
MeqTrees. Available online: http://meqtrees.net, accessed on 30 July 2022.
40
SYMBA. Available online: https://bitbucket.org/M_Janssen/symba
41
MeqSilhouete. Available online: https://github.com/rdeane/MeqSilhouette, accessed on 30 July 2022.
42
VNSIM. Available online: https://github.com/ZhenZHAO/VNSIM, accessed on 30 July 2022.
43
KERN. Available online: https://kernsuite.info, accessed on 30 July 2022.
44
JIVE FlexBuff RAID/JBOD recorders, 4  Gbps, scalable, stationary. Available online: https://www.jive.eu/technical-operations-rd-group, accessed on 30 July 2022.
45
Conduant/MIT Haystack Mark 6, max. 32 Gbps, shippable. Available online: https://www.haystack.mit.edu/mark-6-vlbi-data-system/, accessed on 30 July 2022.
46
Elecs/NAOJ Octadisk2, max. 32 Gbps, shippable. Available online: https://www.elecs.co.jp/product/removable_storage.html, accessed on 30 July 2022.
47
Available online: https://github.com/jive-vlbi/jive5ab accessed on 30 July 2022.
48
Available online: https://github.com/jive-vlbi/etransfer, accessed on 30 July 2022.
49
Available online: https://tsunami-udp.sourceforge.net, accessed on 30 July 2022.
50
Available online: https://www.globus.org/, accessed on 30 July 2022.
51
RFI affects autocorrelation-based calibration steps and strong RFI can also affect the cross-correlations.

References

  1. Gaylard, M.J.; Bietenholz, M.F.; Combrinck, L.; Booth, R.S.; Buchner, S.J.; Fanaroff, B.L.; MacLeod, G.C.; Nicolson, G.D.; Quick, J.F.H.; Stronkhorst, P.; et al. An African VLBI Network of radio telescopes. arXiv 2014, arXiv:1405.7214. [Google Scholar]
  2. Deane, R. Extragalactic VLBI surveys in the MeerKAT era. In Proceedings of the MeerKAT Science: On the Pathway to the SKA, Stellenbosch, South Africa, 25–27 May 2016; p. 17. [Google Scholar]
  3. Matthews, L.D.; Crew, G.B.; Doeleman, S.S.; Lacasse, R.; Saez, A.F.; Alef, W.; Akiyama, K.; Amestica, R.; Anderson, J.M.; Barkats, D.A.; et al. The ALMA Phasing System: A Beamforming Capability for Ultra-high-resolution Science at (Sub)Millimeter Wavelengths. Publ. Astron. Soc. Pac. 2018, 130, 015002. [Google Scholar] [CrossRef] [Green Version]
  4. Venturi, T.; Paragi, Z.; Lindqvist, M.; Bartkiewicz, A.; Beswick, R.; Bogdanović, T.; Brisken, W.; Charlot, P.; Colomer, F.; Conway, J.; et al. VLBI20-30: A scientific roadmap for the next decade—The future of the European VLBI Network. arXiv 2020, arXiv:2007.02347. [Google Scholar]
  5. Barbosa, D.; Coelho, B.; Antón, S.; Bergano, M.; Boekholt, T.; Correia, A.C.M.; Maia, D.; Pandeirada, J.; Ribeiro, V.; Adams, J.; et al. Radio astronomy and Space science in Azores: Enhancing the Atlantic VLBI infrastructure cluster. Adv. Space Res. 2021, 68, 3064–3078. [Google Scholar] [CrossRef]
  6. Whitney, A.R.; Beaudoin, C.J.; Cappallo, R.J.; Corey, B.E.; Crew, G.B.; Doeleman, S.S.; Lapsley, D.E.; Hinton, A.A.; McWhirter, S.R.; Niell, A.E.; et al. Demonstration of a 16 Gbps Station-1 Broadband-RF VLBI System. Publ. Astron. Soc. Pac. 2013, 125, 196. [Google Scholar] [CrossRef]
  7. Vertatschitsch, L.; Primiani, R.; Young, A.; Weintroub, J.; Crew, G.B.; McWhirter, S.R.; Beaudoin, C.; Doeleman, S.; Blackburn, L. R2DBE: A Wideband Digital Backend for the Event Horizon Telescope. Publ. Astron. Soc. Pac. 2015, 127, 1226. [Google Scholar] [CrossRef]
  8. Han, S.T.; Lee, J.W.; Kang, J.; Je, D.H.; Chung, M.H.; Wi, S.O.; Sasao, T.; Wylde, R. Millimeter-wave Receiver Optics for Korean VLBI Network. Int. J. Infrared Millim. Waves 2008, 29, 69–78. [Google Scholar] [CrossRef]
  9. Boccardi, B.; Krichbaum, T.P.; Ros, E.; Zensus, J.A. Radio observations of active galactic nuclei with mm-VLBI. A&A Rev. 2017, 25, 4. [Google Scholar] [CrossRef] [Green Version]
  10. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole. Astrophys. J. Lett. 2019, 875, L1. [Google Scholar] [CrossRef]
  11. Selina, R.J.; Murphy, E.J.; McKinnon, M.; Beasley, A.; Butler, B.; Carilli, C.; Clark, B.; Durand, S.; Erickson, A.; Grammer, W.; et al. Science with an ngVLA: The ngVLA Reference Design. arXiv 2018, arXiv:astro-ph.IM/1810.08197. [Google Scholar]
  12. Murphy, E.J.; Bolatto, A.; Chatterjee, S.; Casey, C.M.; Chomiuk, L.; Dale, D.; de Pater, I.; Dickinson, M.; Francesco, J.D.; Hallinan, G.; et al. Science with an ngVLA: The ngVLA Science Case and Associated Science Requirements. arXiv 2018, arXiv:astro-ph.IM/1810.07524. [Google Scholar]
  13. Dewdney, P.E.; Hall, P.J.; Schilizzi, R.T.; Lazio, T.J.L.W. The Square Kilometre Array. IEEE Proc. 2009, 97, 1482–1496. [Google Scholar] [CrossRef]
  14. Paragi, Z.; Godfrey, L.; Reynolds, C.; Rioja, M.; Deller, A.; Zhang, B.; Gurvits, L.; Bietenholz, M.; Szomoru, A.; Bignall, H.; et al. Very Long Baseline Interferometry with the SKA. Advancing Astrophysics with the Square Kilometre Array (AASKA14). arXiv 2015, arXiv:astro-ph.IM/1412.5971. [Google Scholar]
  15. Doeleman, S.; Blackburn, L.; Dexter, J.; Gomez, J.L.; Johnson, M.D.; Palumbo, D.C.; Weintroub, J.; Farah, J.R.; Fish, V.; Loinard, L.; et al. Studying Black Holes on Horizon Scales with VLBI Ground Arrays. Bull. Am. Astron. Soc. 2019, 51, 256. [Google Scholar]
  16. Kardashev, N.S.; Novikov, I.D.; Lukash, V.N.; Pilipenko, S.V.; Mikheeva, E.V.; Bisikalo, D.V.; Wiebe, D.S.; Doroshkevich, A.G.; Zasov, A.V.; Zinchenko, I.I.; et al. Review of scientific topics for the Millimetron space observatory. Phys. Uspekhi 2014, 57, 1199–1228. [Google Scholar] [CrossRef] [Green Version]
  17. Johnson, M.; Haworth, K.; Pesce, D.W.; Palumbo, D.C.M.; Blackburn, L.; Akiyama, K.; Boroson, D.; Bouman, K.L.; Farah, J.R.; Fish, V.L.; et al. Studying black holes on horizon scales with space-VLBI. Bull. Am. Astron. Soc. 2019, 51, 235. [Google Scholar]
  18. Fish, V.L.; Shea, M.; Akiyama, K. Imaging black holes and jets with a VLBI array including multiple space-based telescopes. Adv. Space Res. 2020, 65, 821–830. [Google Scholar] [CrossRef]
  19. Gurvits, L.I.; Paragi, Z.; Casasola, V.; Conway, J.; Davelaar, J.; Falcke, H.; Fender, R.; Frey, S.; Fromm, C.M.; Miró, C.G.; et al. THEZA: TeraHertz Exploration and Zooming-in for Astrophysics. Exp. Astron. 2021, 51, 559–594. [Google Scholar] [CrossRef]
  20. Marcote, B.; Paragi, Z.; Hessels, J.W.T.; Keimpema, A.; van Langevelde, H.J.; Huang, Y.; Bassa, C.G.; Bogdanov, S.; Bower, G.C.; Burke-Spolaor, S.; et al. The Repeating Fast Radio Burst FRB 121102 as Seen on Milliarcsecond Angular Scales. Astrophys. J. Lett. 2017, 834, L8. [Google Scholar] [CrossRef]
  21. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. Astrophys. J. Lett. 2022, 930, L12. [Google Scholar] [CrossRef]
  22. Gabuzda, D. Evidence for Helical Magnetic Fields Associated with AGN Jets and the Action of a Cosmic Battery. Galaxies 2018, 7, 5. [Google Scholar] [CrossRef] [Green Version]
  23. Gabuzda, D.C. Polarization VLBI observations of AGN jets now and into the future. Adv. Space Res. 2020, 65, 731–738. [Google Scholar] [CrossRef]
  24. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. VII. Polarization of the Ring. Astrophys. J. Lett. 2021, 910, L12. [Google Scholar] [CrossRef]
  25. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. VIII. Magnetic Field Structure near The Event Horizon. Astrophys. J. Lett. 2021, 910, L13. [Google Scholar] [CrossRef]
  26. Gabuzda, D.C. Inherent and Local Magnetic Field Structures in Jets from Active Galactic Nuclei. Galaxies 2021, 9, 58. [Google Scholar] [CrossRef]
  27. Habing, H.J. Circumstellar envelopes and Asymptotic Giant Branch stars. Astron. Astrophys. Rev. 1996, 7, 97–207. [Google Scholar] [CrossRef]
  28. Momjian, E.; Romney, J.D.; Carilli, C.L.; Troland, T.H. Sensitive VLBI Continuum and H I Absorption Observations of NGC 7674: First Scientific Observations with the Combined Array VLBA, VLA, and Arecibo. Astrophys. J. 2003, 597, 809–822. [Google Scholar] [CrossRef]
  29. Goddi, C.; Moscadelli, L.; Alef, W.; Tarchi, A.; Brand, J.; Pani, M. Kinematics of H2O masers in high-mass star forming regions. Astron. Astrophys. 2005, 432, 161–173. [Google Scholar] [CrossRef]
  30. Hachisuka, K.; Brunthaler, A.; Menten, K.M.; Reid, M.J.; Imai, H.; Hagiwara, Y.; Miyoshi, M.; Horiuchi, S.; Sasao, T. Water Maser Motions in W3(OH) and a Determination of Its Distance. Astrophys. J. 2006, 645, 337–344. [Google Scholar] [CrossRef] [Green Version]
  31. Reid, M.J.; Braatz, J.A.; Condon, J.J.; Greenhill, L.J.; Henkel, C.; Lo, K.Y. The Megamaser Cosmology Project. I. Very Long Baseline Interferometric Observations of UGC 3789. Astrophys. J. 2009, 695, 287–291. [Google Scholar] [CrossRef] [Green Version]
  32. Brunthaler, A.; Reid, M.J.; Menten, K.M.; Zheng, X.W.; Bartkiewicz, A.; Choi, Y.K.; Dame, T.; Hachisuka, K.; Immer, K.; Moellenbrock, G.; et al. The Bar and Spiral Structure Legacy (BeSSeL) survey: Mapping the Milky Way with VLBI astrometry. Astron. Nachrichten 2011, 332, 461. [Google Scholar] [CrossRef] [Green Version]
  33. Goddi, C.; Moscadelli, L.; Sanna, A. Infall and outflow within 400 AU from a high-mass protostar. 3D velocity fields from methanol and water masers in AFLG 5142. Astron. Astrophys. 2011, 535, L8. [Google Scholar] [CrossRef]
  34. Reid, M.J.; Braatz, J.A.; Condon, J.J.; Lo, K.Y.; Kuo, C.Y.; Impellizzeri, C.M.V.; Henkel, C. The Megamaser Cosmology Project. IV. A Direct Measurement of the Hubble Constant from UGC 3789. Astrophys. J. 2013, 767, 154. [Google Scholar] [CrossRef]
  35. Goddi, C.; Surcis, G.; Moscadelli, L.; Imai, H.; Vlemmings, W.H.T.; van Langevelde, H.J.; Sanna, A. Measuring magnetic fields from water masers in the synchrotron protostellar jet in W3(H2O). Astron. Astrophys. 2017, 597, A43. [Google Scholar] [CrossRef] [Green Version]
  36. Moscadelli, L.; Sanna, A.; Goddi, C.; Krishnan, V.; Massi, F.; Bacciotti, F. Protostellar Outflows at the EarliesT Stages (POETS). III. H2O masers tracing disk-winds and jets near luminous YSOs. Astron. Astrophys. 2019, 631, A74. [Google Scholar] [CrossRef] [Green Version]
  37. Garrett, M.A.; Muxlow, T.W.B.; Garrington, S.T.; Alef, W.; Alberdi, A.; van Langevelde, H.J.; Venturi, T.; Polatidis, A.G.; Kellermann, K.I.; Baan, W.A.; et al. AGN and starbursts at high redshift: High resolution EVN radio observations of the Hubble Deep Field. Astron. Astrophys. 2001, 366, L5–L8. [Google Scholar] [CrossRef]
  38. Herrera Ruiz, N.; Middelberg, E.; Deller, A.; Norris, R.P.; Best, P.N.; Brisken, W.; Schinnerer, E.; Smolčić, V.; Delvecchio, I.; Momjian, E.; et al. The faint radio sky: VLBA observations of the COSMOS field. Astron. Astrophys. 2017, 607, A132. [Google Scholar] [CrossRef]
  39. Radcliffe, J.F.; Garrett, M.A.; Muxlow, T.W.B.; Beswick, R.J.; Barthel, P.D.; Deller, A.T.; Keimpema, A.; Campbell, R.M.; Wrigley, N. Nowhere to Hide: Radio-faint AGN in GOODS-N field. I. Initial catalogue and radio properties. Astron. Astrophys. 2018, 619, A48. [Google Scholar] [CrossRef] [Green Version]
  40. Spingola, C.; Mckean, J.P.; Deller, A.; Moldon, J. Gravitational lensing at milliarcsecond angular resolution with VLBI observations. In Proceedings of the 14th European VLBI Network Symposium & Users Meeting (EVN 2018), Granada, Spain, 8–11 October 2018; p. 33. [Google Scholar]
  41. Spingola, C.; McKean, J.P.; Lee, M.; Deller, A.; Moldon, J. A novel search for gravitationally lensed radio sources in wide-field VLBI imaging from the mJIVE-20 survey. Mon. Not. R. Astron. Soc. 2019, 483, 2125–2153. [Google Scholar] [CrossRef]
  42. Deller, A.T.; Goss, W.M.; Brisken, W.F.; Chatterjee, S.; Cordes, J.M.; Janssen, G.H.; Kovalev, Y.Y.; Lazio, T.J.W.; Petrov, L.; Stappers, B.W.; et al. Microarcsecond VLBI Pulsar Astrometry with PSRπ II. Parallax Distances for 57 Pulsars. Astrophys. J. 2019, 875, 100. [Google Scholar] [CrossRef] [Green Version]
  43. Duev, D.A.; Molera Calvés, G.; Pogrebenko, S.V.; Gurvits, L.I.; Cimó, G.; Bocanegra Bahamon, T. Spacecraft VLBI and Doppler tracking: Algorithms and implementation. Astron. Astrophys. 2012, 541, A43. [Google Scholar] [CrossRef] [Green Version]
  44. Thompson, A.R.; Moran, J.M.; Swenson, G.W., Jr. Interferometry and Synthesis in Radio Astronomy, 3rd ed.; Springer International Publishing: Berlin/Heidelberg, Germany, 2017. [Google Scholar] [CrossRef] [Green Version]
  45. Morabito, L.K.; Jackson, N.J.; Mooney, S.; Sweijen, F.; Badole, S.; Kukreti, P.; Venkattu, D.; Groeneveld, C.; Kappes, A.; Bonnassieux, E.; et al. Sub-arcsecond imaging with the International LOFAR Telescope. I. Foundational calibration strategy and pipeline. Astron. Astrophys. 2022, 658, A1. [Google Scholar] [CrossRef]
  46. Allan, D.W. Statistics of atomic frequency standards. IEEE Proc. 1966, 54, 221–230. [Google Scholar] [CrossRef] [Green Version]
  47. Primiani, R.A.; Young, K.H.; Young, A.; Patel, N.; Wilson, R.W.; Vertatschitsch, L.; Chitwood, B.B.; Srinivasan, R.; MacMahon, D.; Weintroub, J. SWARM: A 32 GHz Correlator and VLBI Beamformer for the Submillimeter Array. J. Astron. Instrum. 2016, 5, 1641006–1641810. [Google Scholar] [CrossRef] [Green Version]
  48. Kim, H.G.; Han, S.T.; Sohn, B.W.; Oh, S.J.; Je, D.H.; Wi, S.O.; Song, M.G. Construction of the Korean VLBI Network (KVN). In Proceedings of the European VLBI Network on New Developments in VLBI Science and Technology, Toledo, Spain, 12–15 October 2004; pp. 281–284. [Google Scholar]
  49. Lee, S.S.; Byun, D.Y.; Oh, C.S.; Han, S.T.; Je, D.H.; Kim, K.T.; Wi, S.O.; Cho, S.H.; Sohn, B.W.; Kim, J.; et al. Single-Dish Performance of KVN 21 m Radio Telescopes: Simultaneous Observations at 22 and 43 GHz. Publ. Astron. Soc. Pac. 2011, 123, 1398. [Google Scholar] [CrossRef] [Green Version]
  50. Lee, S.S.; Petrov, L.; Byun, D.Y.; Kim, J.; Jung, T.; Song, M.G.; Oh, C.S.; Roh, D.G.; Je, D.H.; Wi, S.O.; et al. Early Science with the Korean VLBI Network: Evaluation of System Performance. Astron. J. 2014, 147, 77. [Google Scholar] [CrossRef]
  51. Kobayashi, H.; Sasao, T.; Kawaguchi, N.; Manabe, S.; Omodaka, T.; Kameya, O.; Shibata, K.M.; Miyaji, T.; Honma, M.; Tamura, Y.; et al. VERA: A New VLBI Instrument Free from the Atmosphere. Astron. Soc. Pac. Conf. Ser. 2003, 306, 367. [Google Scholar]
  52. Lee, S.S.; Oh, C.S.; Roh, D.G.; Oh, S.J.; Kim, J.; Yeom, J.H.; Kim, H.R.; Jung, D.G.; Byun, D.Y.; Jung, T.; et al. A New Hardware Correlator in Korea: Performance Evaluation Using KVN Observations. J. Korean Astron. Soc. 2015, 48, 125–137. [Google Scholar] [CrossRef] [Green Version]
  53. Deller, A.T.; Tingay, S.J.; Bailes, M.; West, C. DiFX: A Software Correlator for Very Long Baseline Interferometry Using Multiprocessor Computing Environments. Publ. Astron. Soc. Pac. 2007, 119, 318–336. [Google Scholar] [CrossRef] [Green Version]
  54. Deller, A.T.; Brisken, W.F.; Phillips, C.J.; Morgan, J.; Alef, W.; Cappallo, R.; Middelberg, E.; Romney, J.; Rottmann, H.; Tingay, S.J.; et al. DiFX-2: A More Flexible, Efficient, Robust, and Powerful Software Correlator. Publ. Astron. Soc. Pac. 2011, 123, 275. [Google Scholar] [CrossRef] [Green Version]
  55. Keimpema, A.; Kettenis, M.M.; Pogrebenko, S.V.; Campbell, R.M.; Cimó, G.; Duev, D.A.; Eldering, B.; Kruithof, N.; van Langevelde, H.J.; Marchal, D.; et al. The SFXC software correlator for very long baseline interferometry: Algorithms and implementation. Exp. Astron. 2015, 39, 259–279. [Google Scholar] [CrossRef] [Green Version]
  56. Napier, P.J.; Bagri, D.S.; Clark, B.G.; Rogers, A.E.E.; Romney, J.D.; Thompson, A.R.; Walker, R.C. The Very Long Baseline Array. IEEE Proc. 1994, 82, 658–672. [Google Scholar] [CrossRef]
  57. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. II. Array and Instrumentation. Astrophys. J. Lett. 2019, 875, L2. [Google Scholar] [CrossRef]
  58. Porcas, R. A history of the EVN. In Proceedings of the 10th European VLBI Network Symposium and EVN Users Meeting: VLBI and the New Generation of Radio Arrays, Manchester, UK, 20–24 September 2010; Volume 10, p. 11. [Google Scholar]
  59. Surkis, I.; Ken, V.; Kurdubova, Y.; Mishina, N.; Mishin, V.; Shantyr, V.; Zhuravov, D.; Zimovsky, V. The RASFX VGOS GPU Based Software Correlator. Trans. Inst. Appl. Astron. Ras 2017, 41, 123–126. [Google Scholar]
  60. Shuygina, N.; Ivanov, D.; Ipatov, A.; Gayazov, I.; Marshalov, D.; Melnikov, A.; Kurdubov, S.; Vasilyev, M.; Ilin, G.; Skurikhina, E.; et al. Russian VLBI network “Quasar”: Current status and outlook. Geod. Geodyn. 2019, 10, 150–156. [Google Scholar] [CrossRef]
  61. Zhang, F.; Zhao, C.; Han, S.; Ma, F.; Xiang, D. GPU-Based Parallel Implementation of VLBI Correlator for Deep Space Exploration System. Remote Sens. 2021, 13, 1226. [Google Scholar] [CrossRef]
  62. Rioja, M.; Dodson, R.; Asaki, Y.; Hartnett, J.; Tingay, S. The Impact of Frequency Standards on Coherence in VLBI at the Highest Frequencies. Astrophys. J. 2012, 144, 121. [Google Scholar] [CrossRef] [Green Version]
  63. Clivati, C.; Aiello, R.; Bianco, G.; Bortolotti, C.; De Natale, P.; Di Sarno, V.; Maddaloni, P.; Maccaferri, G.; Mura, A.; Negusini, M.; et al. Common-clock very long baseline interferometry using a coherent optical fiber link. Optica 2020, 7, 1031–1037. [Google Scholar] [CrossRef]
  64. Krehlik, P.; Buczek, Ł.; Kołodziej, J.; Lipiński, M.; Śliwczyński, .; Nawrocki, J.; Nogaś, P.; Marecki, A.; Pazderski, E.; Ablewski, P.; et al. Fibre-optic delivery of time and frequency to VLBI station. Astron. Astrophys. 2017, 603, A48. [Google Scholar] [CrossRef] [Green Version]
  65. Likhachev, S.F.; Kostenko, V.I.; Girin, I.A.; Andrianov, A.S.; Rudnitskiy, A.G.; Zharov, V.E. Software Correlator for Radioastron Mission. J. Astron. Instrum. 2017, 6, 1750004–1750131. [Google Scholar] [CrossRef] [Green Version]
  66. Gordon, D. CALC: The Next Upgrade. In Proceedings of the IVS 2004 General Meeting Proceedings, Ottawa, ON, Canada, 9–11 February 2004; pp. 265–266. [Google Scholar]
  67. Sekido, M.; Fukushima, T. A VLBI Delay Model for Radio Sources at a Finite Distance. J. Geod. 2006, 80, 137–149. [Google Scholar] [CrossRef]
  68. Wootten, A.; Thompson, A.R. The Atacama Large Millimeter/Submillimeter Array. IEEE Proc. 2009, 97, 1463–1471. [Google Scholar] [CrossRef] [Green Version]
  69. Martí-Vidal, I.; Roy, A.; Conway, J.; Zensus, A.J. Calibration of mixed-polarization interferometric observations. Tools for the reduction of interferometric data from elements with linear and circular polarization receivers. Astron. Astrophys. 2016, 587, A143. [Google Scholar] [CrossRef] [Green Version]
  70. Goddi, C.; Martí-Vidal, I.; Messias, H.; Crew, G.B.; Herrero-Illana, R.; Impellizzeri, V.; Rottmann, H.; Wagner, J.; Fomalont, E.; Matthews, L.D.; et al. Calibration of ALMA as a Phased Array. ALMA Observations During the 2017 VLBI Campaign. Publ. Astron. Soc. Pac. 2019, 131, 075003. [Google Scholar] [CrossRef] [Green Version]
  71. Alef, W.; Porcas, R.W. VLBI fringe-fitting with antenna-based residuals. Astron. Astrophys. 1986, 168, 365–368. [Google Scholar]
  72. Schwab, F.R.; Cotton, W.D. Global fringe search techniques for VLBI. AJ 1983, 88, 688–694. [Google Scholar] [CrossRef]
  73. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. III. Data Processing and Calibration. Astrophys. J. Lett. 2019, 875, L3. [Google Scholar] [CrossRef]
  74. Natarajan, I.; Deane, R.; van Bemmel, I.; van Langevelde, H.J.; Small, D.; Kettenis, M.; Paragi, Z.; Smirnov, O.; Szomoru, A. A probabilistic approach to phase calibration - I. Effects of source structure on fringe-fitting. Mon. Not. R. Astron. Soc. 2020, 496, 801–813. [Google Scholar] [CrossRef]
  75. Doi, A.; Fujisawa, K.; Habe, A.; Honma, M.; Kawaguchi, N.; Kobayashi, H.; Murata, Y.; Omodaka, T.; Sudou, H.; Takaba, H. Bigradient Phase Referencing. Publ. Astron. Soc. Jpn. 2006, 58, 777–785. [Google Scholar] [CrossRef]
  76. Fomalont, E.B.; Kopeikin, S. Phase Referencing Using Several Calibrator Sources. In Proceedings of the 6th EVN Symposium, Bonn, Germany, 25–28 June 2002; p. 53. [Google Scholar]
  77. Rioja, M.J.; Dodson, R.; Orosz, G.; Imai, H.; Frey, S. MultiView High Precision VLBI Astrometry at Low Frequencies. AJ 2017, 153, 105. [Google Scholar] [CrossRef] [Green Version]
  78. Hyland, L.J.; Reid, M.J.; Ellingsen, S.P.; Rioja, M.J.; Dodson, R.; Orosz, G.; Masson, C.R.; McCallum, J.M. Inverse Multiview. I. Multicalibrator Inverse Phase Referencing for Microarcsecond Very Long Baseline Interferometry Astrometry. Astrophys. J. 2022, 932, 52. [Google Scholar] [CrossRef]
  79. Reid, M.; Honma, M. Microarcsecond Radio Astrometry. Annu. Rev. Astron. Astrophys. 2014, 52, 339–372. [Google Scholar] [CrossRef] [Green Version]
  80. Rogers, A.E.E.; Doeleman, S.S.; Moran, J.M. Fringe Detection Methods for Very Long Baseline Arrays. AJ 1995, 109, 1391. [Google Scholar] [CrossRef]
  81. Rioja, M.J.; Dodson, R.; Kamohara, R.; Colomer, F.; Bujarrabal, V.; Kobayashi, H. Relative Astrometry of the J = 1–˃0, v = 1 and v = 2 SiO Masers toward R Leonis Minoris Using VERA. Publ. Astron. Soc. Jpn. 2008, 60, 1031. [Google Scholar] [CrossRef] [Green Version]
  82. Rioja, M.; Dodson, R. High-precision Astrometric Millimeter Very Long Baseline Interferometry Using a New Method for Atmospheric Calibration. AJ 2011, 141, 114. [Google Scholar] [CrossRef] [Green Version]
  83. Rioja, M.J.; Dodson, R.; Jung, T.; Sohn, B.W.; Byun, D.Y.; Agudo, I.; Cho, S.H.; Lee, S.S.; Kim, J.; Kim, K.T.; et al. Verification of the Astrometric Performance of the Korean VLBI Network, Using Comparative SFPR Studies with the VLBA at 14/7 mm. AJ 2014, 148, 84. [Google Scholar] [CrossRef]
  84. Rioja, M.J.; Dodson, R.; Jung, T.; Sohn, B.W. The Power of Simultaneous Multifrequency Observations for mm-VLBI: Astrometry up to 130 GHz with the KVN. Astron. J. 2015, 150, 202. [Google Scholar] [CrossRef] [Green Version]
  85. Rioja, M.J.; Dodson, R. Precise radio astrometry and new developments for the next-generation of instruments. Astron. Astrophys. Rev. 2020, 28, 6. [Google Scholar] [CrossRef]
  86. Greisen, E.W. AIPS; VLA; VLBA. In Information Handling in Astronomy—Historical Vistas; Heck, A., Ed.; Astrophysics and Space Science Library; Springer: Berlin/Heidelberg, Germany, 2003; Volume 285, p. 109. [Google Scholar] [CrossRef]
  87. Kettenis, M.; van Langevelde, H.J.; Reynolds, C.; Cotton, B. ParselTongue: AIPS Talking Python. In Astronomical Data Analysis Software and Systems XV; Gabriel, C., Arviset, C., Ponz, D., Enrique, S., Eds.; Astronomical Society of the Pacific Conference Series; Astronomical Society of the Pacific: San Francisco, CA, USA, 2006; Volume 351, p. 497. [Google Scholar]
  88. McMullin, J.P.; Waters, B.; Schiebel, D.; Young, W.; Golap, K. CASA Architecture and Applications. In Astronomical Data Analysis Software and Systems XVI; Shaw, R.A., Hill, F., Bell, D.J., Eds.; Astronomical Society of the Pacific Conference Series; Astronomical Society of the Pacific: San Francisco, CA, USA, 2007; Volume 376, p. 127. [Google Scholar]
  89. Bean, B.; Bhatnagar, S.; Castro, S.; Donovan Meyer, J.; Emonts, B.; Garcia, E.; Garwood, R.; Golap, K.; Gonzalez Villalba, J.; Harris, P.; et al. CASA, the Common Astronomy Software Applications for Radio Astronomy. Publ. Astron. Soc. Pac. 2022, in press. Available online: https://arxiv.org/abs/2210.02276 (accessed on 30 July 2022).
  90. Thompson, A.R.; Clark, B.G.; Wade, C.M.; Napier, P.J. The Very Large Array. ApJS 1980, 44, 151–167. [Google Scholar] [CrossRef]
  91. van Haarlem, M.P.; Wise, M.W.; Gunst, A.W.; Heald, G.; McKean, J.P.; Hessels, J.W.T.; de Bruyn, A.G.; Nijboer, R.; Swinbank, J.; Fallows, R.; et al. LOFAR: The LOw-Frequency ARray. Astron. Astrophys. 2013, 556, A2. [Google Scholar] [CrossRef] [Green Version]
  92. Goddi, C.; Falcke, H.; Kramer, M.; Rezzolla, L.; Brinkerink, C.; Bronzwaer, T.; Davelaar, J.R.J.; Deane, R.; De Laurentis, M.; Desvignes, G.; et al. BlackHoleCam: Fundamental physics of the galactic center. Int. J. Mod. Phys. D 2017, 26, 17300014. [Google Scholar] [CrossRef]
  93. Van Bemmel, I.M.; Kettenis, M.; Small, D.; Janssen, M.; Moellenbrock, G.A.; Petry, D.; Goddi, C.; Linford, J.D.; Rygl, K.L.J.; Liuzzo, E.; et al. CASA on the fringe – Development of VLBI processing capabilities for CASA. Publ. Astron. Soc. Pac. 2022, in press. Available online: https://arxiv.org/abs/2210.02275 (accessed on 30 July 2022).
  94. Hamaker, J.P.; Bregman, J.D.; Sault, R.J. Understanding radio polarimetry. I. Mathematical foundations. A&AS 1996, 117, 137–147. [Google Scholar]
  95. Smirnov, O.M. Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism. Astron. Astrophys. 2011, 527, A106. [Google Scholar] [CrossRef] [Green Version]
  96. Smirnov, O.M. Revisiting the radio interferometer measurement equation. II. Calibration and direction-dependent effects. Astron. Astrophys. 2011, 527, A107. [Google Scholar] [CrossRef]
  97. Smirnov, O.M. Revisiting the radio interferometer measurement equation. III. Addressing direction-dependent effects in 21 cm WSRT observations of 3C 147. Astron. Astrophys. 2011, 527, A108. [Google Scholar] [CrossRef] [Green Version]
  98. Smirnov, O.M. Revisiting the radio interferometer measurement equation. IV. A generalized tensor formalism. Astron. Astrophys. 2011, 531, A159. [Google Scholar] [CrossRef] [Green Version]
  99. Rocklin, M. Dask: Parallel Computation with Blocked algorithms and Task Scheduling. In Proceedings of the 14th Python in Science Conference, Austin, TX, USA, 6–12 July 2015; pp. 126–132. [Google Scholar] [CrossRef] [Green Version]
  100. Hoyer, S.; Hamman, J. xarray: N-D labeled Arrays and Datasets in Python. J. Open Res. Softw. 2017, 5, 10. [Google Scholar] [CrossRef] [Green Version]
  101. Janssen, M.; Goddi, C.; Falcke, H.; van Rossum, D.; van Bemmel, I.; Kettenis, M.; Small, D.; Marti-Vidal, I. RPICARD: A CASA-based Calibration Pipeline for VLBI Data. In Proceedings of the 14th European VLBI Network Symposium & Users Meeting (EVN 2018), Granada, Spain, 8–11 October 2018; p. 80. [Google Scholar]
  102. Janssen, M.; Goddi, C.; van Bemmel, I.M.; Kettenis, M.; Small, D.; Liuzzo, E.; Rygl, K.; Martí-Vidal, I.; Blackburn, L.; Wielgus, M.; et al. rPICARD: A CASA-based calibration pipeline for VLBI data. Calibration and imaging of 7 mm VLBA observations of the AGN jet in M 87. Astron. Astrophys. 2019, 626, A75. [Google Scholar] [CrossRef]
  103. Blecher, T.; Deane, R.; Bernardi, G.; Smirnov, O. MEQSILHOUETTE: A mm-VLBI observation and signal corruption simulator. Mon. Not. R. Astron. Soc. 2017, 464, 143–151. [Google Scholar] [CrossRef] [Green Version]
  104. Natarajan, I.; Deane, R.; Martí-Vidal, I.; Roelofs, F.; Janssen, M.; Wielgus, M.; Blackburn, L.; Blecher, T.; Perkins, S.; Smirnov, O.; et al. MeqSilhouette v2: Spectrally resolved polarimetric synthetic data generation for the event horizon telescope. Mon. Not. R. Astron. Soc. 2022, 512, 490–504. [Google Scholar] [CrossRef]
  105. Pardo, J.R.; Cernicharo, J.; Serabyn, E. Atmospheric transmission at microwaves (ATM): An improved model for millimeter/submillimeter applications. IEEE Trans. Antennas Propag. 2001, 49, 1683–1694. [Google Scholar] [CrossRef] [Green Version]
  106. Davies, J.G.; Anderson, B.; Morison, I. The Jodrell Bank radio-linked interferometer network. Nature 1980, 288, 64–66. [Google Scholar] [CrossRef]
  107. Offringa, A.R.; van de Gronde, J.J.; Roerdink, J.B.T.M. A morphological algorithm for improving radio-frequency interference detection. Astron. Astrophys. 2012, 539, A95. [Google Scholar] [CrossRef] [Green Version]
  108. Offringa, A.R.; McKinley, B.; Hurley-Walker, N.; Briggs, F.H.; Wayth, R.B.; Kaplan, D.L.; Bell, M.E.; Feng, L.; Neben, A.R.; Hughes, J.D.; et al. WSCLEAN: An implementation of a fast, generic wide-field imager for radio astronomy. Mon. Not. R. Astron. Soc. 2014, 444, 606–619. [Google Scholar] [CrossRef]
  109. Keimpema, A.; Kettenis, M.; Small, D.; Dijkema, T.J.; Szomoru, A. Efficient Remote Interactive Pipelines Using CASA and Jupyter. In Astronomical Data Analysis Software and Systems XXIX; Pizzo, R., Deul, E.R., Mol, J.D., de Plaa, J., Verkouter, H., Eds.; Astronomical Society of the Pacific Conference Series; Astronomical Society of the Pacific: San Francisco, CA, USA, 2020; Volume 527, p. 579. [Google Scholar]
  110. Whitney, A.R.; Cappallo, R.; Aldrich, W.; Anderson, B.; Bos, A.; Casse, J.; Goodman, J.; Parsley, S.; Pogrebenko, S.; Schilizzi, R.; et al. Mark 4 VLBI correlator: Architecture and algorithms. Radio Sci. 2004, 39, RS1007. [Google Scholar] [CrossRef]
  111. Blackburn, L.; Chan, C.k.; Crew, G.B.; Fish, V.L.; Issaoun, S.; Johnson, M.D.; Wielgus, M.; Akiyama, K.; Barrett, J.; Bouman, K.L.; et al. EHT-HOPS Pipeline for Millimeter VLBI Data Reduction. Astrophys. J. 2019, 882, 23. [Google Scholar] [CrossRef]
  112. Issaoun, S.; Johnson, M.D.; Blackburn, L.; Brinkerink, C.D.; Mościbrodzka, M.; Chael, A.; Goddi, C. Martí-Vidal, I.; Wagner, J.; Doeleman, S.S.; et al. The Size, Shape, and Scattering of Sagittarius A* at 86 GHz: First VLBI with ALMA. Astrophys. J. 2019, 871, 30. [Google Scholar] [CrossRef] [Green Version]
  113. Kim, J.Y.; Krichbaum, T.P.; Broderick, A.E.; Wielgus, M.; Blackburn, L.; Gómez, J.L.; Johnson, M.D.; Bouman, K.L.; Chael, A.; Akiyama, K.; et al. Event Horizon Telescope imaging of the archetypal blazar 3C 279 at an extreme 20 microarcsecond resolution. Astron. Astrophys. 2020, 640, A69. [Google Scholar] [CrossRef]
  114. Janssen, M.; Falcke, H.; Kadler, M.; Ros, E.; Wielgus, M.; Akiyama, K.; Baloković, M.; Blackburn, L.; Bouman, K.L.; Chael, A.; et al. Event Horizon Telescope observations of the jet launching and collimation in Centaurus A. Nat. Astron. 2021, 5, 1017–1028. [Google Scholar] [CrossRef]
  115. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. II. EHT and Multiwavelength Observations, Data Processing, and Calibration. Astrophys. J. Lett. 2022, 930, L13. [Google Scholar] [CrossRef]
  116. Petrov, L.; Kovalev, Y.Y.; Fomalont, E.B.; Gordon, D. The Very Long Baseline Array Galactic Plane Survey—VGaPS. AJ 2011, 142, 35. [Google Scholar] [CrossRef] [Green Version]
  117. Ulich, B.L.; Haas, R.W. Absolute calibration of millimeter-wavelength spectral lines. ApJS 1976, 30, 247–258. [Google Scholar] [CrossRef]
  118. Pearson, T.J.; Readhead, A.C.S. Image Formation by Self-Calibration in Radio Astronomy. ARA&A 1984, 22, 97–130. [Google Scholar]
  119. Pearson, T.J. Caltech VLBI Analysis Programs, California Institute of Technology. Bull. Am. Astron. Soc. 1991, 23, 991–992. [Google Scholar]
  120. Issaoun, S.; Johnson, M.D.; Blackburn, L.; Mościbrodzka, M.; Chael, A.; Falcke, H. VLBI imaging of black holes via second moment regularization. Astron. Astrophys. 2019, 629, A32. [Google Scholar] [CrossRef]
  121. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. IV. Imaging the Central Supermassive Black Hole. Astrophys. J. Lett. 2019, 875, L4. [Google Scholar] [CrossRef]
  122. Högbom, J.A. Aperture Synthesis with a Non-Regular Distribution of Interferometer Baselines. A&AS 1974, 15, 417. [Google Scholar]
  123. Clark, B.G. An efficient implementation of the algorithm ‘CLEAN’. Astron. Astrophys. 1980, 89, 377. [Google Scholar]
  124. Cornwell, T.J. Multiscale CLEAN Deconvolution of Radio Synthesis Images. IEEE J. Sel. Top. Signal Process. 2008, 2, 793–801. [Google Scholar] [CrossRef] [Green Version]
  125. Conway, J.E.; Cornwell, T.J.; Wilkinson, P.N. Multi-frequency synthesis: A new technique in radio interferometrie imaging. Mon. Not. R. Astron. Soc. 1990, 246, 490. [Google Scholar]
  126. Sault, R.J.; Wieringa, M.H. Multi-frequency synthesis techniques in radio interferometric imaging. A&AS 1994, 108, 585–594. [Google Scholar]
  127. Likhachev, S. Multi-Frequency Imaging for VLBI. In Future Directions in High Resolution Astronomy; Romney, J., Reid, M., Eds.; Astronomical Society of the Pacific Conference Series; Astronomical Society of the Pacific: San Francisco, CA, USA, 2005; Volume 340, p. 608. [Google Scholar]
  128. Rau, U.; Cornwell, T.J. A multi-scale multi-frequency deconvolution algorithm for synthesis imaging in radio interferometry. Astron. Astrophys. 2011, 532, A71. [Google Scholar] [CrossRef]
  129. Cornwell, T.J.; Evans, K.F. A simple maximum entropy deconvolution algorithm. Astron. Astrophys. 1985, 143, 77–83. [Google Scholar]
  130. Mizuno, Y. GRMHD Simulations and Modeling for Jet Formation and Acceleration Region in AGNs. arXiv 2022, arXiv:2201.12608. [Google Scholar] [CrossRef]
  131. van der Gucht, J.; Davelaar, J.; Hendriks, L.; Porth, O.; Olivares, H.; Mizuno, Y.; Fromm, C.M.; Falcke, H. Deep Horizon: A machine learning network that recovers accreting black hole parameters. Astron. Astrophys. 2020, 636, A94. [Google Scholar] [CrossRef]
  132. Yao-Yu Lin, J.; Pesce, D.W.; Wong, G.N.; Uppili Arasanipalai, A.; Prather, B.S.; Gammie, C.F. VLBInet: Radio Interferometry Data Classification for EHT with Neural Networks. arXiv 2021, arXiv:2110.07185. [Google Scholar]
  133. Fromm, C.M.; Younsi, Z.; Baczko, A.; Mizuno, Y.; Porth, O.; Perucho, M.; Olivares, H.; Nathanail, A.; Angelakis, E.; Ros, E.; et al. Using evolutionary algorithms to model relativistic jets. Application to NGC 1052. Astron. Astrophys. 2019, 629, A4. [Google Scholar] [CrossRef] [Green Version]
  134. Broderick, A.E.; Gold, R.; Karami, M.; Preciado-López, J.A.; Tiede, P.; Pu, H.Y.; Akiyama, K.; Alberdi, A.; Alef, W.; Asada, K.; et al. THEMIS: A Parameter Estimation Framework for the Event Horizon Telescope. Astrophys. J. 2020, 897, 139. [Google Scholar] [CrossRef]
  135. Jennison, R.C. A phase sensitive interferometer technique for the measurement of the Fourier transforms of spatial brightness distributions of small angular extent. Mon. Not. R. Astron. Soc. 1958, 118, 276. [Google Scholar] [CrossRef] [Green Version]
  136. Blackburn, L.; Pesce, D.W.; Johnson, M.D.; Wielgus, M.; Chael, A.A.; Christian, P.; Doeleman, S.S. Closure Statistics in Interferometric Data. Astrophys. J. 2020, 894, 31. [Google Scholar] [CrossRef]
  137. Broderick, A.E.; Pesce, D.W. Closure Traces: Novel Calibration-insensitive Quantities for Radio Astronomy. Astrophys. J. 2020, 904, 126. [Google Scholar] [CrossRef]
  138. Thyagarajan, N.; Nityananda, R.; Samuel, J. Invariants in copolar interferometry: An Abelian gauge theory. Phys. Rev. D 2022, 105, 043019. [Google Scholar] [CrossRef]
  139. Readhead, A.C.S.; Wilkinson, P.N. The mapping of compact radio sources from VLBI data. Astrophys. J. 1978, 223, 25–36. [Google Scholar] [CrossRef]
  140. Cotton, W.D. A method of mapping compact structure in radio sources using VLBI observations. AJ 1979, 84, 1122–1128. [Google Scholar] [CrossRef]
  141. Martí-Vidal, I.; Mus, A.; Janssen, M.; de Vicente, P.; González, J. Polarization calibration techniques for the new-generation VLBI. Astron. Astrophys. 2021, 646, A52. [Google Scholar] [CrossRef]
  142. Shepherd, M.C. Difmap: An Interactive Program for Synthesis Imaging. In Astronomical Data Analysis Software and Systems VI; Hunt, G., Payne, H., Eds.; Astronomical Society of the Pacific Conference Series; Astronomical Society of the Pacific: San Francisco, CA, USA, 1997; Volume 125, p. 77. [Google Scholar]
  143. Offringa, A.R.; Smirnov, O. An optimized algorithm for multiscale wideband deconvolution of radio astronomical images. Mon. Not. R. Astron. Soc. 2017, 471, 301–316. [Google Scholar] [CrossRef]
  144. Chael, A.A.; Johnson, M.D.; Bouman, K.L.; Blackburn, L.L.; Akiyama, K.; Narayan, R. Interferometric Imaging Directly with Closure Phases and Closure Amplitudes. Astrophys. J. 2018, 857, 23. [Google Scholar] [CrossRef]
  145. Chael, A.A.; Johnson, M.D.; Narayan, R.; Doeleman, S.S.; Wardle, J.F.C.; Bouman, K.L. High-resolution Linear Polarimetric Imaging for the Event Horizon Telescope. Astrophys. J. 2016, 829, 11. [Google Scholar] [CrossRef]
  146. Akiyama, K.; Ikeda, S.; Pleau, M.; Fish, V.L.; Tazaki, F.; Kuramochi, K.; Broderick, A.E.; Dexter, J.; Mościbrodzka, M.; Gowanlock, M.; et al. Superresolution Full-polarimetric Imaging for Radio Interferometry with Sparse Modeling. AJ 2017, 153, 159. [Google Scholar] [CrossRef] [Green Version]
  147. Akiyama, K.; Kuramochi, K.; Ikeda, S.; Fish, V.L.; Tazaki, F.; Honma, M.; Doeleman, S.S.; Broderick, A.E.; Dexter, J.; Mościbrodzka, M.; et al. Imaging the Schwarzschild-radius-scale Structure of M87 with the Event Horizon Telescope Using Sparse Modeling. Astrophys. J. 2017, 838, 1. [Google Scholar] [CrossRef]
  148. Martí-Vidal, I.; Vlemmings, W.H.T.; Muller, S.; Casey, S. UVMULTIFIT: A versatile tool for fitting astronomical radio interferometric data. Astron. Astrophys. 2014, 563, A136. [Google Scholar] [CrossRef] [Green Version]
  149. Bezanson, J.; Karpinski, S.; Shah, V.B.; Edelman, A. Julia: A Fast Dynamic Language for Technical Computing. arXiv 2012, arXiv:1209.5145. [Google Scholar]
  150. Tiede, P.; Broderick, A.E.; Palumbo, D.C.M. Variational Image Feature Extraction for the Event Horizon Telescope. Astrophys. J. 2022, 925, 122. [Google Scholar] [CrossRef]
  151. Sun, H.; Bouman, K.L.; Tiede, P.; Wang, J.J.; Blunt, S.; Mawet, D. α-deep Probabilistic Inference (α-DPI): Efficient Uncertainty Quantification from Exoplanet Astrometry to Black Hole Feature Extraction. Astrophys. J. 2022, 932, 99. Available online: http://xxx.lanl.gov/abs/2201.08506 (accessed on 30 July 2022). [CrossRef]
  152. Sun, H.; Bouman, K.L.; Tiede, P.; Wang, J.J.; Blunt, S.; Mawet, D. alpha-Deep Probabilistic Inference (alpha-DPI): Efficient uncertainty quantification from exoplanet astrometry to black hole feature extraction. arXiv 2022, arXiv:2201.08506. [Google Scholar]
  153. [Event Horizon Telescope Collaboration]. First Sagittarius A* Event Horizon Telescope Results. III. Imaging of the Galactic Center Supermassive Black Hole. Astrophys. J. Lett. 2022, 930, L14. [Google Scholar] [CrossRef]
  154. Homan, D.C.; Wardle, J.F.C. Detection and Measurement of Parsec-Scale Circular Polarization in Four AGNS. AJ 1999, 118, 1942–1962. [Google Scholar] [CrossRef]
  155. Goddi, C.; Moscadelli, L.; Alef, W.; Brand, J. EVN observations of H2O masers towards the high-mass young stellar object in AFGL 5142. Astron. Astrophys. 2004, 420, 929–936. [Google Scholar] [CrossRef]
  156. Felli, M.; Spencer, R.E. Very Long Baseline Interferometry: Techniques and Applications: Proceedings of the NATO Advanced Study Institute on VLBI; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1988. [Google Scholar]
  157. Matthews, L.D.; Greenhill, L.J.; Goddi, C.; Chandler, C.J.; Humphreys, E.M.L.; Kunz, M.W. A Feature Movie of SiO Emission 20-100 AU from the Massive Young Stellar Object Orion Source I. Astrophys. J. 2010, 708, 80–92. [Google Scholar] [CrossRef]
  158. Strom, R. What is the primary beam response of an interferometer with unequal elements? In Proceedings of the European VLBI Network on New Developments in VLBI Science and Technology, Toledo, Spain, 12–15 October 2004; pp. 273–274. [Google Scholar]
  159. Chi, S.; Barthel, P.D.; Garrett, M.A. Deep, wide-field, global VLBI observations of the Hubble deep field north (HDF-N) and flanking fields (HFF). Astron. Astrophys. 2013, 550, A68. [Google Scholar] [CrossRef] [Green Version]
  160. Morgan, J.S.; Mantovani, F.; Deller, A.T.; Brisken, W.; Alef, W.; Middelberg, E.; Nanni, M.; Tingay, S.J. VLBI imaging throughout the primary beam using accurate UV shifting. Astron. Astrophys. 2011, 526, A140. [Google Scholar] [CrossRef] [Green Version]
  161. van der Tol, S.; Veenboer, B.; Offringa, A.R. Image Domain Gridding: A fast method for convolutional resampling of visibilities. Astron. Astrophys. 2018, 616, A27. [Google Scholar] [CrossRef]
  162. Cornwell, T.J.; Perley, R.A. Radio-interferometric imaging of very large fields. The problem of non-coplanar arrays. Astron. Astrophys. 1992, 261, 353–364. [Google Scholar]
  163. Cornwell, T.J.; Golap, K.; Bhatnagar, S. The Noncoplanar Baselines Effect in Radio Interferometry: The W-Projection Algorithm. IEEE J. Sel. Top. Signal Process. 2008, 2, 647–657. [Google Scholar] [CrossRef] [Green Version]
  164. Middelberg, E.; Deller, A.T.; Norris, R.P.; Fotopoulou, S.; Salvato, M.; Morgan, J.S.; Brisken, W.; Lutz, D.; Rovilos, E. Mosaiced wide-field VLBI observations of the Lockman Hole/XMM. Astron. Astrophys. 2013, 551, A97. [Google Scholar] [CrossRef] [Green Version]
  165. Radcliffe, J.F.; Garrett, M.A.; Beswick, R.J.; Muxlow, T.W.B.; Barthel, P.D.; Deller, A.T.; Middelberg, E. Multi-source self-calibration: Unveiling the microJy population of compact radio sources. Astron. Astrophys. 2016, 587, A85. [Google Scholar] [CrossRef]
  166. Noordam, J.E.; Smirnov, O.M. The MeqTrees software system and its use for third-generation calibration of radio interferometers. Astron. Astrophys. 2010, 524, A61. [Google Scholar] [CrossRef] [Green Version]
  167. Roelofs, F.; Janssen, M.; Natarajan, I.; Deane, R.; Davelaar, J.; Olivares, H.; Porth, O.; Paine, S.N.; Bouman, K.L.; [The Event Horizon Telescope Collaboration]; et al. SYMBA: An end-to-end VLBI synthetic data generation pipeline. Simulating Event Horizon Telescope observations of M 87. Astron. Astrophys. 2020, 636, A5. [Google Scholar] [CrossRef]
  168. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. VI. The Shadow and Mass of the Central Black Hole. Astrophys. J. Lett. 2019, 875, L6. [Google Scholar] [CrossRef]
  169. Wielgus, M.; Akiyama, K.; Blackburn, L.; Chan, C.-k.; Dexter, J.; Doeleman, S.S.; Fish, V.L.; Issaoun, S.; Johnson, M.D.; Krichbaum, T.P.; et al. Monitoring the Morphology of M87* in 2009-2017 with the Event Horizon Telescope. Astrophys. J. 2020, 901, 67. [Google Scholar] [CrossRef]
  170. Roelofs, F.; Fromm, C.M.; Mizuno, Y.; Davelaar, J.; Janssen, M.; Younsi, Z.; Rezzolla, L.; Falcke, H. Black hole parameter estimation with synthetic very long baseline interferometry data from the ground and from space. Astron. Astrophys. 2021, 650, A56. [Google Scholar] [CrossRef]
  171. [Event Horizon Telescope Collaboration]. First M87 Event Horizon Telescope Results. V. Physical Origin of the Asymmetric Ring. Astrophys. J. Lett. 2019, 875, L5. [Google Scholar] [CrossRef]
  172. Zhao, Z.; An, T.; Lao, B. VLBI Network SIMulator: An Integrated Simulation Tool for Radio Astronomers. J. Korean Astron. Soc. 2019, 52, 207–216. [Google Scholar]
  173. Molenaar, G.; Smirnov, O. Kern. Astron. Comput. 2018, 24, 45. [Google Scholar] [CrossRef]
Figure 1. VLBI observation and data processing. Per-telescope baseband and calibration data are collected during VLBI observations. Calibration data are bundled with the final visibility data released to the PI. Delay models for far-field VLBI are commonly derived via NASA GSFC CALC 9/10, for near-field VLBI via c5++ or others. Dashed boxes indicate optional processing, such as pulsar processing, space VLBI antenna support, multi phase center for wide-field VLBI, spectral windowing (in SFXC), or custom channelization (zoom bands in DiFX). Correlation itself can run on hardware ranging from Raspberry PI and HPC clusters, to Cloud platforms.
Figure 1. VLBI observation and data processing. Per-telescope baseband and calibration data are collected during VLBI observations. Calibration data are bundled with the final visibility data released to the PI. Delay models for far-field VLBI are commonly derived via NASA GSFC CALC 9/10, for near-field VLBI via c5++ or others. Dashed boxes indicate optional processing, such as pulsar processing, space VLBI antenna support, multi phase center for wide-field VLBI, spectral windowing (in SFXC), or custom channelization (zoom bands in DiFX). Correlation itself can run on hardware ranging from Raspberry PI and HPC clusters, to Cloud platforms.
Universe 08 00527 g001
Figure 2. Flowchart of typical continuum VLBI data reduction steps shown in the order as implemented in the rPICARD pipeline (Section 3.2.1). Feed rotation angle effects are calibrated on-the-fly here.
Figure 2. Flowchart of typical continuum VLBI data reduction steps shown in the order as implemented in the rPICARD pipeline (Section 3.2.1). Feed rotation angle effects are calibrated on-the-fly here.
Universe 08 00527 g002
Figure 3. Example fringe solution diagnostics from rPICARD of M87 observed at 43 GHz with the VLBA (project code BW0106). Dashed lines show scan boundaries. Crosses and circles correspond to RCP and LCP, respectively. Top: Gross per-scan solutions. Bottom: Residual atmospheric corrections on adaptive segmentation timescales after the solutions from the top panel have been applied to the data.
Figure 3. Example fringe solution diagnostics from rPICARD of M87 observed at 43 GHz with the VLBA (project code BW0106). Dashed lines show scan boundaries. Crosses and circles correspond to RCP and LCP, respectively. Top: Gross per-scan solutions. Bottom: Residual atmospheric corrections on adaptive segmentation timescales after the solutions from the top panel have been applied to the data.
Universe 08 00527 g003
Figure 4. Example of effective (opacity-corrected) system temperatures measured by three VLBA antennas from observations of M87 at 43 GHz (project code BW0106). Pie Town system temperatures show slight increases at the beginning and end of the observing track, where the source is at a low elevation. The source us set at the end of the track for the North Liberty station, which causes an increases in system temperature. Bad weather is responsible for the highly variable system temperature measurements at Saint Croix.
Figure 4. Example of effective (opacity-corrected) system temperatures measured by three VLBA antennas from observations of M87 at 43 GHz (project code BW0106). Pie Town system temperatures show slight increases at the beginning and end of the observing track, where the source is at a low elevation. The source us set at the end of the track for the North Liberty station, which causes an increases in system temperature. Bad weather is responsible for the highly variable system temperature measurements at Saint Croix.
Universe 08 00527 g004
Figure 5. VLBI image reconstruction and underlying visibility amplitude data reproduced from Janssen et al. [114]. The results are from EHT observations of Centaurus A in 2017 at 228 GHz. The image shown in the left panel was reconstructed with the eht-imaging regularized maximum-likelihood method. The right panel shows the visibility amplitudes as a function of ( u , v ) -distance projected along the jet direction (top) and perpendicular (bottom) to the jet direction.
Figure 5. VLBI image reconstruction and underlying visibility amplitude data reproduced from Janssen et al. [114]. The results are from EHT observations of Centaurus A in 2017 at 228 GHz. The image shown in the left panel was reconstructed with the eht-imaging regularized maximum-likelihood method. The right panel shows the visibility amplitudes as a function of ( u , v ) -distance projected along the jet direction (top) and perpendicular (bottom) to the jet direction.
Universe 08 00527 g005
Figure 6. Synthetic EHT observations of M87*. Top: Ray-traced general relativistic magnetohydrodynamics (GRMHD) simulations [171] from the EHT library with different black hole spins, electron temperature parameters, and magnetic field states (“MAD” versus “SANE”). Bottom: SYMBA synthetic observations corresponding to these models. Reproduced from Event Horizon Telescope Collaboration and et al. [10].
Figure 6. Synthetic EHT observations of M87*. Top: Ray-traced general relativistic magnetohydrodynamics (GRMHD) simulations [171] from the EHT library with different black hole spins, electron temperature parameters, and magnetic field states (“MAD” versus “SANE”). Bottom: SYMBA synthetic observations corresponding to these models. Reproduced from Event Horizon Telescope Collaboration and et al. [10].
Universe 08 00527 g006
Table 1. Basic terms, definitions, and concepts are listed here, which new VLBI users might be unfamiliar with or which might have slightly different meanings in other related works.
Table 1. Basic terms, definitions, and concepts are listed here, which new VLBI users might be unfamiliar with or which might have slightly different meanings in other related works.
NameMeaning
Telescope frontendEquipment “directly attached to the front of a telescope”, used to detect and encode the sky signal. Receivers (e.g., bolometers) are frontend equipment.
Telescope backendEquipment in the telescope that is used to process and store the signal from the frontend. Block Downconverters, which mix the sky signal down to a lower frequency range (heterodyning); analog-to-digital converters, which digitize the data; and data recorders, which store the digitized measurements on hard drives, are backend equipment.
Baseband dataThe recorded data at a telescope that will be used for the correlation. More precisely, the filtered, down-converted, sampled, and quantized electric field measurements stored in the backends.
Signal stabilizationDescribed in Section 3: The collection of all post-correlation calibration measures, excluding the a priori flux density calibration (Section 4). The signal stabilization is often referred to as fringe-fitting, but it also involves additional steps, e.g., corrections for bandpass responses and corrections for atmospheric phase turbulence.
DelayResidual post-correlation phase-slopes as a function of frequency (e.g., due to atmospheric path length differences). To be corrected in the signal stabilization step.
RateResidual post-correlation phase-slopes as a function of time (e.g., due to Doppler shifts). To be corrected in the signal stabilization step.
Fringe-fit FFTThe fast Fourier transform step of common fringe-fitting algorithms. Transforms the visibilities from time, frequency space into a rate, delay space, where the peaks are to be found. The height of the peaks corresponds to the strength of the signal.
Low ν An observing frequency below 20 GHz.
High ν An observing frequency above 20 GHz.
Allan deviationA measure of frequency stability [46]. An easy-to-follow derivation and description of the Allan deviation equation is given in Section 9.5.1 of Thompson et al. [44].
VEX file“VLBI EXperiment” file, which describes the VLBI setups and observing schedules in a standardized text format (VEX File Definition/Example. Available online: https://vlbi.org/wp-content/uploads/2019/03/vex-definition-15b1.pdf, accessed on 30 July 2022).
FITS-IDI file“FITS Interferometry Data Interchange Convention” standardized file format for visibility and VLBI metadata built upon the upon the standard FITS format (FITS-IDI format. Available online: https://fits.gsfc.nasa.gov/registry/fitsidi/AIPSMEM114.PDF, accessed on 30 July 2022).
ANTAB file“Antenna table”, which contains station gain and system temperature information in a simple text file format (ANTAB format. Available online: http://www.aips.nrao.edu/cgi-bin/ZXHLP2.PL?ANTAB, accessed on 30 July 2022).
Table 2. Major VLBI software packages and their capabilities. AIPS and Difmap are still supported by their originators; requested features and bug fixes can be implemented. All other packages are actively being developed and supported by larger communities and institutions.
Table 2. Major VLBI software packages and their capabilities. AIPS and Difmap are still supported by their originators; requested features and bug fixes can be implemented. All other packages are actively being developed and supported by larger communities and institutions.
NameCapabilitiesLinks and References
DiFX Correlationhttps://ascl.net/1102.024,
https://svn.atnf.csiro.au/difx, [53,54]
SFXC Correlationhttps://svn.astron.nl/sfxc, [55]
AIPS,
ParselTongue
Signal stabilization, flux density calibration, imaging, synthetic data generationhttps://ascl.net/9911.003, https://ascl.net/1208.020,
http://www.aips.nrao.edu,
https://www.jive.eu/jivewiki/doku.php?id=parseltongue:parseltongue, [86,87]
CASASignal stabilization, flux density calibration, imaginghttps://ascl.net/1107.013, https://casa.nrao.edu
   - VPIPE https://github.com/jradcliffe5/VLBI_pipeline
   - rPICARD https://ascl.net/1905.015,
https://bitbucket.org/M_Janssen/picard,
https://hub.docker.com/r/mjanssen2308/casavlbi, [88,89,93,101,102]
(EHT-)HOPS Signal stabilization,
flux denisity calibration
https://www.haystack.mit.edu/tech/vlbi/hops.html,
https://github.com/sao-eht/eat, [110,111]
Difmap Imaging, model fittinghttps://ascl.net/1103.001,
ftp://ftp.astro.caltech.edu/pub/difmap/difmap.html, [142]
eht-imaging Imaging, model fitting,
synthetic data generation
https://ascl.net/1904.004,
https://github.com/achael/eht-imaging, [144,145]
SMILI Imaginghttps://ascl.net/1904.005,
https://github.com/astrosmili/smili, [146,147]
WSClean Imaginghttps://ascl.net/1408.023,
https://gitlab.com/aroffringa/wsclean, [108]
UV-MULTIFIT Geometric model fittinghttps://ascl.net/1402.017,
https://launchpad.net/uvmultifit, [148]
SYMBA Synthetic data generationhttps://bitbucket.org/M_Janssen/symba, [167]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Janssen, M.; Radcliffe, J.F.; Wagner, J. Software and Techniques for VLBI Data Processing and Analysis. Universe 2022, 8, 527. https://doi.org/10.3390/universe8100527

AMA Style

Janssen M, Radcliffe JF, Wagner J. Software and Techniques for VLBI Data Processing and Analysis. Universe. 2022; 8(10):527. https://doi.org/10.3390/universe8100527

Chicago/Turabian Style

Janssen, Michael, Jack F. Radcliffe, and Jan Wagner. 2022. "Software and Techniques for VLBI Data Processing and Analysis" Universe 8, no. 10: 527. https://doi.org/10.3390/universe8100527

APA Style

Janssen, M., Radcliffe, J. F., & Wagner, J. (2022). Software and Techniques for VLBI Data Processing and Analysis. Universe, 8(10), 527. https://doi.org/10.3390/universe8100527

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop