Run GPU-Powered Black Hole Interferometry with NVIDIA cuPhoton - NVIDIA Developer

Imaging black holes presents unique challenges due to their distance and tiny apparent size in the sky. For instance, the black hole M87* spans just 42 microarcseconds. Since black holes emit no light, telescopes capture radio emissions from superheated gas swirling outside the event horizon. The intense gravity of black holes bends these light paths, creating the characteristic bright ring that frames the central shadow.

The Event Horizon Telescope (EHT) addresses these challenges by linking 13 radio observatory stations across Hawaii, Chile, Spain, Mexico, and the South Pole through very-long-baseline interferometry (VLBI). These stations record millimeter and submillimeter radio signals as synchronized voltage streams that are combined to reconstruct detailed images. The next-generation EHT (ngEHT) aims to add more observing sites to improve image quality and track changes around black holes over time.

GPU Acceleration Breakthrough

NVIDIA's cuPhoton.xBeam performs GPU-based interferometry, organizing recorded data for parallel reads and correlating signals from telescope pairs. This process produces complex visibility measurements for each pair of stations, frequency channel, and time interval. These measurements serve as inputs for downstream calibration and image reconstruction processes like Gaussian splatting, which models emitting gas as glowing blobs near black holes.

The Imaging Challenge

The EHT records radio signals at up to 64 gigabits per second per telescope, generating petabytes of data. Correlating signals between telescope pairs produces visibilities that sample Fourier components of the sky's brightness. Earth's rotation creates sparse measurement tracks, and reconstruction algorithms combine these samples into an image. For targets like Sagittarius A*, the source changes during observation, complicating reconstruction methods that assume a static source.

Light passing near a black hole can orbit it before escaping, forming a thin photon ring near the shadow's edge. The proposed Black Hole Explorer (BHEX) aims to measure this ring, with future ngEHT observations at 345 GHz potentially detecting its signature through polarization measurements.

How cuPhoton.xBeam Works

cuPhoton.xBeam performs five key steps: ingesting VDIF frames and building an index for efficient access; applying delay and fringe corrections to align wavefronts; channelizing data through FFT; cross-multiplying and accumulating station pairs; and writing visibilities in DiFX's SWIN format or xBeam's own format. The system uses delays specified in the job file rather than calculating them from station positions.

Data moves from disk to GPU memory through either standard POSIX reads or the GPUDirect Storage API. The phase correction, FFT, and cross-multiplication all run on the GPU. For jobs using multiple GPUs, the system splits scans in time and merges partial outputs. This architecture enables unprecedented processing speeds while maintaining compatibility with existing telescope data formats and processing pipelines.