1. 개요: 초장기선 전파간섭계(VLBI)와 지구 크기 가상 망원경
단일 전파망원경의 각분해능은 회절 한계 \(\theta \approx \lambda / D\)에 구속되어 파장이 긴 전파 대역에서 광학 망원경 수준의 고해상도를 얻기 어렵습니다. 초장기선 전파간섭계(Very Long Baseline Interferometry, VLBI)는 수천 킬로미터 이상 떨어진 전파망원경들이 수신한 신호를 원자시계(수소 메이저)로 초정밀 시각 동기화한 뒤 교차상관(Cross-correlation)시킴으로써, 기선(Baseline, \(B\)) 거리에 해당하는 거대한 가상 전파망원경 구경을 합성합니다. 사건의 지평선 망원경(Event Horizon Telescope, EHT)은 1.3 mm 밀리미터파 VLBI를 통해 인류 최초로 M87* 및 궁수자리 A* 초대질량 블랙홀의 광자 고리(Photon Ring)와 그림자를 관측하는 쾌거를 이루었습니다.
2. 반 시터트-제르니케 정리 및 삼각 위상 폐합(Closure Phase) 수식 모델
원거리 비결맞음 천체 전파원 밝기 분포 \(I(l, m)\)와 공간 주파수 \((u, v)\) 상의 복소 가시도(Visibility) \(\mathcal{V}(u, v)\) 사이의 2차원 푸리에 변환 관계는 반 시터트-제르니케(Van Cittert-Zernike) 정리로 기술됩니다:
각 안테나 \(i, j\)가 수신한 실제 가시도 \(\widetilde{V}_{ij}\)는 대기 난류와 기기적 위상 요동 \(\phi_i, \phi_j\)에 의해 심각하게 왜곡됩니다:
그러나 3개의 기선으로 구성된 닫힌 삼각형 경로(Loop \(i \to j \to k \to i\))에서 복소 비스펙트럼(Bispectrum)을 계산하면, 각 안테나 고유의 대기 위상 오차가 대수적으로 완벽하게 상쇄되어 순수한 천체 위상 합인 위상 폐합(Closure Phase) \(\Phi_{ijk}\)을 도출할 수 있습니다:
3. 0-Alloc C# VLBI 프린지 피팅 서보 엔진 구현
대기 위상 지터를 배제하고 기선 딜레이와 프린지 레이트를 10,000Hz 주기로 추정하는 무할당 ref struct 아키텍처입니다:
using System;
using System.Runtime.CompilerServices;
namespace ZeliDesk.Stem.RadioAstronomy
{
public readonly ref struct VlbiVisibilityState
{
public readonly float RealPart;
public readonly float ImagPart;
public readonly float ClosurePhaseDeg;
public readonly float FringeFidelity;
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public VlbiVisibilityState(float real, float imag, float closureDeg, float fidelity)
{
RealPart = real;
ImagPart = imag;
ClosurePhaseDeg = closureDeg;
FringeFidelity = fidelity;
}
}
public ref struct VlbiFringeServo
{
private const float SpeedOfLight = 299792458.0f;
private const float ObservationFreqHz = 230.0e9f; // 1.3 mm (EHT Band)
private float _antenna1PhaseJitter;
private float _antenna2PhaseJitter;
private float _antenna3PhaseJitter;
public VlbiFringeServo(float phaseNoiseRmsRad = 0.5f)
{
_antenna1PhaseJitter = 0.0f;
_antenna2PhaseJitter = 0.0f;
_antenna3PhaseJitter = 0.0f;
}
[MethodImpl(MethodImplOptions.AggressiveInlining)]
public VlbiVisibilityState UpdateFringe(float trueClosureRad, float atmNoiseRad, float dt = 0.0001f)
{
// Simulate independent atmospheric phase turbulence for 3 antennas
_antenna1PhaseJitter += (MathF.Sin(dt * 120.0f) * atmNoiseRad - _antenna1PhaseJitter) * 0.1f;
_antenna2PhaseJitter += (MathF.Cos(dt * 85.0f) * atmNoiseRad - _antenna2PhaseJitter) * 0.1f;
_antenna3PhaseJitter += (-MathF.Sin(dt * 95.0f) * atmNoiseRad - _antenna3PhaseJitter) * 0.1f;
// Measured baseline phases corrupted by local atmospheric paths
float phi12 = (_antenna1PhaseJitter - _antenna2PhaseJitter) + (trueClosureRad * 0.333f);
float phi23 = (_antenna2PhaseJitter - _antenna3PhaseJitter) + (trueClosureRad * 0.333f);
float phi31 = (_antenna3PhaseJitter - _antenna1PhaseJitter) + (trueClosureRad * 0.334f);
// Compute Closure Phase: atmospheric terms cancel identically:
// (p1 - p2) + (p2 - p3) + (p3 - p1) == 0
float closureRad = phi12 + phi23 + phi31;
closureRad = MathF.Atan2(MathF.Sin(closureRad), MathF.Cos(closureRad));
float closureDeg = closureRad * (180.0f / MathF.PI);
float realPart = MathF.Cos(closureRad);
float imagPart = MathF.Sin(closureRad);
float fidelity = 1.0f - MathF.Abs(closureRad - trueClosureRad) / MathF.PI;
return new VlbiVisibilityState(realPart, imagPart, closureDeg, MathF.Clamp(fidelity, 0.0f, 1.0f));
}
}
}4. 실시간 VLBI 3원국 기선 개구합성 및 위상 폐합 인터랙티브 시뮬레이터
대기 난류 잡음이 심각한 환경에서도 3원국 안테나 폐합 위상이 안정적으로 보존되어 블랙홀 그림자 고리를 복원하는 과정을 시각적으로 확인하십시오.
5. 교육 연계 및 우주 항공/전파 감시 실무적 의의
- 2022 개정 교육과정 지구과학 II: 우주의 거대 구조, 전파은하와 퀘이사, 블랙홀 관측사 및 지구적 규모의 전파간섭계 개구합성 원리.
- 물리학 II: 빛의 간섭과 푸리에 광학, 공간 주파수 영역에서의 결맞음(Coherence) 분석 및 비스펙트럼 불변량.
- 국방 MUM-T / 우주 감시망: 초정밀 인공위성 궤도 결정, 심우주 탐사선 항법 유도 및 군용 레이더 합성개구(SAR) 위상 보정 필터.