NVIDIA xBeam Moves Radio Interferometry Correlation to GPUs, But Does Not Perform Imaging
A NVIDIA technical blog post details how the cuPhoton.xBeam FX correlation pipeline maps to GPU parallelism: performing FFTs on voltage streams from each station to obtain frequency channels, then computing and accumulating cross-products of complex spectra for station pairs. The article reports that processing 10 TB of data on 64 GPUs in a Vera Rubin NVL72 system takes just 70 seconds, approximately 460 times faster than an x86 CPU baseline, but explicitly states that xBeam does not calculate delays based on station coordinates nor perform imaging calibration; downstream tasks still rely on DiFX for verification, rPICARD for calibration, and eht-imaging for reconstruction.
The Event Horizon Telescope (EHT) links 13 radio observatory stations located in Hawaii, Chile, Spain, Mexico, and Antarctica via very-long-baseline interferometry (VLBI). Each telescope records high-speed voltage streams in two polarizations at rates up to 64 Gbps, generating petabytes of data. Correlation processing—multiplying voltage signals from any two stations channel-by-channel and accumulating them—is one of the most time-consuming steps in the entire data pipeline. In a technical blog post released in November 2026, NVIDIA introduced how cuPhoton.xBeam moves this step onto GPUs.
xBeam adopts an FX correlator architecture, using the same mathematical framework as the open-source CPU correlator DiFX. Its pipeline consists of five steps: reading VDIF frames and building a .vdix index to support random access; performing integer alignment and fringe rotation on each station according to the delay model in the job file; executing an FFT on the voltage stream to split it into frequency channels; calculating the product of the complex spectrum of each station pair's channel with its conjugate and accumulating over the integration time; and finally writing visibilities in DiFX SWIN or xbeam-v1 format. Among these, the FFT and cross-multiplication are computationally intensive and independent across channels, making them naturally suitable for the GPU's parallel execution model.