Correlation (VLBI)
From Recorded Voltages to Complex Visibilities
In a connected-element array the antenna signals travel down cables to a central building where they are combined in real time. VLBI cannot do this: the baselines are too long, so each station digitizes its intermediate-frequency band, time-stamps every sample against a local hydrogen-maser clock, and records to disk or transfers over fiber. The correlator is the offline machine that finally brings these streams together. Its job is to estimate, for each pair of stations, the time-averaged product of their electric-field voltages, which by the Van Cittert-Zernike theorem equals a single Fourier component of the sky brightness distribution sampled at the projected baseline coordinate (u, v).
The challenge is that the two recordings are not aligned in time. The wavefront from the source reaches the more distant antenna later by the geometric delay τg = −(b · s) / c, where b is the baseline vector and s the unit vector to the source. This delay can reach tens of milliseconds on intercontinental baselines and changes continuously as the Earth rotates. The correlator computes τg from a model that also folds in station clock offsets, polar motion, solid-Earth tides, and tropospheric and ionospheric path lengths, then applies an integer-sample shift plus a fractional-sample delay interpolation to one stream.
Once delay-aligned, the residual interferometer phase still spins at the fringe rate, the time derivative of the delay scaled by the observing frequency. The correlator multiplies one station by a complex sinusoid at the negative fringe rate so that, over the accumulation period, the cross-product phase is nearly stationary and survives integration. The accumulated complex number, after a quantization (Van Vleck) correction, is the calibrated visibility that downstream imaging and fringe-fitting software consume.
The Cross-Correlation and Fringe Relations
τg = −(b · s) / c (seconds)
Fringe rate (phase spin to be removed):
ffringe = ν × (dτg / dt) ≈ ν × (ωE b⊥ / c) (Hz)
Measured complex visibility on baseline ij:
Vij(u,v) ≈ 〈 vi(t) × vj*(t − τg) e−j2πffringet 〉
Coherence integration limit:
Δt < 1 / (2π ffringe) so the phase turns < 1 rad per accumulation
Where ν = observing frequency, ωE = Earth rotation rate (7.29 × 10−5 rad/s), b⊥ = projected baseline, c = speed of light, and vj* denotes the delayed, conjugated station-j voltage. Example: at ν = 8.4 GHz with a 6,000 km baseline, ffringe can exceed 10 kHz, demanding millisecond-class fringe rotation before accumulation.
Correlator Architectures and Quantization
| Parameter | XF (lag) correlator | FX correlator | Notes |
|---|---|---|---|
| Order of operations | Cross-multiply lags, then FFT | FFT per station, then cross-multiply | FX dominates modern systems |
| Compute scaling | ∝ N2 × Nlag | ∝ N log N + N2 per channel | FX wins for large N stations |
| Typical platform | Custom ASIC / FPGA | CPU & GPU cluster (DiFX) | Software correlation now standard |
| 1-bit sampling loss | ≈ 1.96 dB (recovers 0.637 of correlation) | Van Vleck correction applied | |
| 2-bit sampling loss | ≈ 0.55 dB (recovers 0.881) | Common modern compromise | |
| Recorded data rate | 0.5 to 32 Gbps per station | VLBA, EVN, EHT class systems | |
Frequently Asked Questions
What is the difference between an FX and an XF VLBI correlator?
An XF design cross-multiplies lagged time samples first and accumulates a lag spectrum, then Fourier transforms it into channels; its cost grows as N2 × Nlag. An FX design Fourier transforms each station first, then cross-multiplies the per-channel spectra, scaling roughly as N log N plus N2 per channel. FX wins decisively for large station counts, which is why software correlators like DiFX and the EHT and VLBA systems are FX machines on CPU or GPU clusters.
Why must a delay model and fringe rotation be applied before integrating?
The wavefront reaches the two antennas at different times, and Earth rotation makes that geometric delay τg change continuously. The correlator shifts one station by the integer-plus-fractional delay and counter-rotates its phase at the fringe rate ν(dτg/dt), which can exceed several kHz at centimeter wavelengths. Without this, the spinning phase averages to zero over even a 1-second integration and the fringe is lost. Models include clock offsets, Earth orientation, and atmospheric path terms.
How does coarse one-bit or two-bit quantization affect correlation sensitivity?
To hold down data volume, samples are quantized coarsely. One-bit (two-level) Nyquist sampling recovers only 0.637 of the correlation amplitude, a 1.96 dB hit, while two-bit (four-level) sampling with optimum thresholds recovers about 0.881, only 0.55 dB. The correlator applies the Van Vleck correction to convert the measured digital coefficient back to true correlation. Two-bit is the usual modern choice because the wider recordable bandwidth it enables outweighs its small extra loss.