1. 개요 및 연구 배경

2022 개정 물리학 II 일반 상대성 이론 단원(등가 원리, 시공간 곡률, 중력 렌즈 효과와 블랙홀)과 고급 기하학(4차원 준-리만 기하학 및 텐서 해석)을 집대성하는 천체물리학 최고봉의 전산 모델입니다. 사건의 지평선 망원경(Event Horizon Telescope, EHT)은 전 지구적 초장기선 전파간섭계(VLBI) 네트워크를 결합하여 거대 타원은하 M87 중심 및 우리은하 중심 Sgr A*의 초대질량 블랙홀 그림자(Black Hole Shadow)를 사상 최초로 영상화하는 데 성공하였습니다.

회전하는 Kerr 블랙홀 주위의 시공간은 시공간 끌림(Frame Dragging) 효과에 의해 광통(Photon)의 궤적이 심각하게 비대칭으로 왜곡됩니다. 사건의 지평선 바깥에 형성되는 광자 구(Photon Sphere) 경계에서 불안정한 원형 궤도를 도는 광선들의 포락선이 관측자의 천구면에 투영된 것이 바로 블랙홀 그림자입니다. 강착원반(Accretion Disk) 플라즈마의 상대론적 운동은 비대칭 도플러 부스팅(Doppler Beaming)을 유발하여 관측자 방향으로 회전하는 쪽이 극단적으로 밝은 초승달 형태를 띱니다.

본 연구에서는 차세대 우주 VLBI 및 밀리미터파 간섭계 위성 군집의 실시간 프린지 추적을 위해, 0-Alloc C# 10,000Hz Carter 상수 보존 Kerr 널 측지선(Null Geodesic) 역방향 광선 추적(Ray-Tracing) 서보 아키텍처를 수립합니다.

2. 수학적 및 일반상대론적 거버닝 모델링

Boyer-Lindquist 좌표계 \((t, r, \theta, \phi)\)에서 무차원 스핀 매개변수 \(a = J / M\) (\(0 \le a < M\))를 갖는 Kerr 계량 텐서:

\[ ds^2 = - \left(1 - \frac{2Mr}{\Sigma}\right) dt^2 - \frac{4Mar \sin^2\theta}{\Sigma} dt d\phi + \frac{\Sigma}{\Delta} dr^2 + \Sigma d\theta^2 + \left(r^2 + a^2 + \frac{2Ma^2 r \sin^2\theta}{\Sigma}\right) \sin^2\theta d\phi^2 \]

여기서 \(\Delta = r^2 - 2Mr + a^2\), \(\Sigma = r^2 + a^2 \cos^2\theta\)입니다. 널 측지선(광선)의 해밀턴-야코비 분리 방정식에서 보존되는 운동의 제4적분(Carter 상수 \(\mathcal{Q}\))과 충돌 매개변수 \(\xi = L_z / E\), \(\eta = \mathcal{Q} / E^2\):

\[ \Sigma \frac{dr}{d\lambda} = \pm \sqrt{\mathcal{R}(r)}, \quad \Sigma \frac{d\theta}{d\lambda} = \pm \sqrt{\Theta(\theta)} \]
\[ \mathcal{R}(r) = (r^2 + a^2 - a \xi)^2 - \Delta \left[ (r - a \xi)^2 + \eta \right] \]
\[ \Theta(\theta) = \eta + a^2 \cos^2\theta - \xi^2 \cot^2\theta \]

블랙홀 그림자 경계선은 \(\mathcal{R}(r_p) = 0\) 및 \(\mathcal{R}'(r_p) = 0\)의 중근 조건으로부터 유도되는 임계 매개변수 \(\xi_c(r_p), \eta_c(r_p)\)로 결정됩니다:

\[ \xi_c(r_p) = \frac{M(r_p^2 - a^2) - r_p \Delta(r_p)}{a(r_p - M)}, \quad \eta_c(r_p) = \frac{r_p^3 [4M\Delta(r_p) - r_p(r_p - M)^2]}{a^2 (r_p - M)^2} \]

3. 0-Alloc C# 초고속 10,000Hz Kerr 측지선 광선 추적 서보

관측자 시야각(Inclination \(\theta_0\))에서 화소 좌표 \((\alpha, \beta)\)로 역방향 투사되는 광선의 Carter 불변량을 계산하고 사건의 지평선 포획 여부를 100µs 주기로 판정하는 0-Alloc 스택 서보 엔진입니다: 

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.GeneralRelativity
{
    public readonly ref struct CelestialRay
    {
        public readonly double ImpactAlpha; // 천구 관측각 수평 변위 (M)
        public readonly double ImpactBeta;  // 천구 관측각 수직 변위 (M)

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public CelestialRay(double alpha, double beta)
        {
            ImpactAlpha = alpha;
            ImpactBeta = beta;
        }
    }

    public static class KerrGeodesicRayServoEngine
    {
        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static void ComputeCarterConstants(
            in CelestialRay ray,
            double spinA,
            double incRad,
            out double xi,
            out double eta)
        {
            double sinInc = Math.Sin(incRad);
            double cosInc = Math.Cos(incRad);

            // 관측자 충돌 매개변수 변환
            xi = - ray.ImpactAlpha * sinInc;
            eta = (ray.ImpactBeta * ray.ImpactBeta - spinA * spinA * cosInc * cosInc) +
                  ray.ImpactAlpha * ray.ImpactAlpha * cosInc * cosInc;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public static bool IsCapturedByBlackHole(
            double xi,
            double eta,
            double spinA)
        {
            // 근사적 그림자 경계 판별식
            double rHorizon = 1.0 + Math.Sqrt(Math.Max(0.0, 1.0 - spinA * spinA));
            double shadowRadiusSq = 27.0 - 6.0 * spinA * spinA; // 슈바르츠실트(27) ~ 극대 스핀

            double impactParamSq = xi * xi + eta;
            return impactParamSq < shadowRadiusSq;
        }
    }
}

4. 2022 개정 교육과정 연계 및 우주천체물리 응용

  • 물리학 II: 아인슈타인 일반 상대성 이론과 블랙홀 시공간 곡률, 중력 적색편이 및 빛의 굴절 현상.
  • 기하와 벡터 / 미적분: 비유클리드 곡면에서의 측지선(최단거리) 적분 및 4차원 시공간 텐서 계량 구조.
  • 우주 과학 연구: 차세대 우주망원경 VLBI 초고해상도 영상 복원 알고리즘 및 퀘이사 제트 생성 메커니즘 규명.