블랙홀

물리 계산물리 시뮬레이션 마지막 수정: 2026-08-30 04:14:22

1. 개요[편집]

블랙홀 Black Hole
정의하는 것	사건지평선 — 빠져나올 수 없는 인과 경계
슈바르츠실트 반지름	r_s = 2GM/c² (태양질량 ≈ 2.95 km)
광자구	r = 3GM/c² = 1.5 r_s
모발	질량 · 각운동량 · 전하, 끝
수치 형식	ADM 3+1 → BSSN · 일반화 조화 · Z4c
돌파 연도	2005 (이진 병합 안정 적분)
온도	T ∝ 1/M (태양질량 ≈ 6×10⁻⁸ K)

방정식이 무한대를 뱉으면 보통은 계산이 틀린 것이다. 여기서는 좌표가 틀린 것일 수도 있고, 시공간이 진짜로 끝난 것일 수도 있다. 그 둘을 구분하는 게 이 문서의 절반이다.

블랙홀은 중력이 너무 강해 어떤 신호도 바깥으로 탈출할 수 없는 시공간 영역이며, 그 영역의 경계인 사건지평선(event horizon)으로 정의된다. “빛보다 빠르지 않으면 못 나온다”는 흔한 설명은 결론만 맞다. 정확히는 지평선 안쪽에서 미래 광원뿔이 전부 안쪽을 향하므로, 바깥으로 나가는 것이 미래에서 과거로 가는 것과 같은 종류의 불가능이 된다.

심위키가 블랙홀을 다루는 이유는 천문학이 아니라 계산 때문이다. 아인슈타인 장방정식은 10개의 성분을 가진 비선형 결합 편미분방정식이고, 좌표 자유도(게이지)가 4개 남아 있으며, 해가 유한 시간에 특이점을 만든다. 수치해석이 싫어하는 성질을 전부 갖춘 방정식이라, 이걸 안정적으로 적분하는 데 인류가 40년을 썼다. 그리고 그 결과물이 지금 중력파 검출기의 템플릿 뱅크와 사건지평선망원경(EHT)의 영상 해석에 매일 쓰이고 있다.

2. 사건지평선과 슈바르츠실트 반지름[편집]

1916년 카를 슈바르츠실트가 구대칭 진공해를 찾았다. 이것이 아인슈타인 방정식의 첫 정확해다.

ds^2 = -\left(1-\frac{r_s}{r}\right)c^2dt^2 + \frac{dr^2}{1-r_s/r} + r^2 d\Omega^2, \qquad r_s = \frac{2GM}{c^2}

$r_s$ 가 슈바르츠실트 반지름이다. 태양질량이면 2.95 km, 지구질량이면 8.9 mm. 즉 태양을 반지름 3 km 짜리 공으로 압축하면 블랙홀이 된다.

주의할 점은 사건지평선의 정의 자체는 전역적이라는 것이다. “미래 널 무한대의 인과적 과거의 경계”이므로, 지평선이 어디인지 알려면 원리적으로 시공간 전체의 미래를 알아야 한다. 수치 시뮬레이션은 미래를 모른 채 전진하므로 이 정의를 그대로 쓸 수 없고, 대신 겉보기 지평선(apparent horizon, 바깥으로 향하는 널 측지선 다발의 팽창률이 0이 되는 국소적으로 정의되는 폐곡면)을 매 시간단계에 찾는다. 이것이 표준 진단 도구이고, 병합 순간의 판정 기준이며, 최종 블랙홀의 질량·스핀을 뽑아내는 통로다. 겉보기 지평선은 항상 사건지평선 안쪽이거나 일치한다.

측지선 쪽에서 나오는 또 하나의 반지름이 광자구 $r = 3GM/c^2 = 1.5\,r_s$ 다. 여기서 빛은 원궤도를 돈다(불안정 궤도라 밀면 떨어지거나 날아간다). 멀리서 본 임계 충돌계수는 $b_c = 3\sqrt3\,GM/c^2 \approx 5.20\,GM/c^2$ 이고, 이 값이 만드는 원이 “그림자”의 가장자리다. 지평선 크기가 아니라 광자구가 렌즈로 확대된 크기라는 점이 자주 오해되는 부분.

3. 좌표특이점과 진짜 특이점[편집]

슈바르츠실트 계량은 $r=r_s$ 에서 $g_{tt}=0$ , $g_{rr}=\infty$ 로 폭발한다. 여기가 무슨 일이 벌어지는 곳인가? 아무 일도 벌어지지 않는다. 좌표가 나쁠 뿐이다.

이걸 판별하는 표준 도구가 곡률 스칼라다. 좌표에 의존하지 않는 양이므로, 이게 유한하면 기하는 멀쩡한 것이다. 슈바르츠실트의 크레치만 스칼라는

R_{abcd}R^{abcd} = \frac{48G^2M^2}{c^4 r^6}

이고 $r=r_s$ 에서 완벽하게 유한하다. 반면 $r=0$ 에서는 발산한다. $r=r_s$ 는 좌표특이점, $r=0$ 은 곡률특이점이라는 결론이 나온다.

좌표를 갈아 끼우면 지평선이 사라진다는 것을 명시적으로 보인 것이 에딩턴-핑켈스타인 좌표(널 좌표를 쓰면 계량이 지평선을 매끄럽게 통과한다)와 크루스칼-셰케레시 좌표(1960, 크루스칼과 셰케레시가 독립적으로)다. 크루스칼 좌표는 최대 확장을 준다 — 두 개의 외부 영역, 블랙홀 내부, 화이트홀 내부로 이루어진 유명한 마름모 그림이 그것이다.¹

계산 실무에서 이건 취향 문제가 아니다. 슈바르츠실트 좌표나 등방 좌표를 쓰면 지평선 위에서 계량 성분이 발산해 코드가 즉시 죽는다. 그래서 강착 시뮬레이션과 영상 계산은 거의 전부 지평선 관통 좌표(커-슈바르츠실트 계열)를 쓴다. 격자를 지평선 안쪽까지 깔아 두면, 안쪽 경계에서 정보가 바깥으로 나올 수 없다는 사실 자체가 물리적으로 정당한 유출 경계조건이 된다. 이게 이 분야에서 경계 조건을 다루는 가장 우아한 트릭이다.

4. 커 해 — 각운동량, 에르고 영역, 무모정리[편집]

현실의 블랙홀은 거의 다 회전한다. 각운동량을 가진 진공해는 1963년 로이 커가 찾았고, 스핀 매개변수 $a = J/(Mc)$ 로 쓰면(기하단위 $G=c=1$ )

r_\pm = M \pm \sqrt{M^2-a^2}, \qquad r_E(\theta) = M + \sqrt{M^2 - a^2\cos^2\theta}

가 각각 지평선과 정적 한계면이다. $a>M$ 이면 지평선이 사라지고 벌거벗은 특이점이 되는데, 이런 일은 안 일어난다는 것이 우주 검열 가설(펜로즈, 1969)이다 — 여전히 미해결.

두 면 사이의 영역이 에르고 영역이다. 여기서는 시간 방향 킬링 벡터가 공간꼴이 되어, 어떤 관측자도 정지해 있을 수 없다. 틀 끌림(frame dragging)이 광속을 넘어서기 때문이 아니라, 정지 상태의 세계선이 더 이상 시간꼴이 아니기 때문이다. 에르고 영역에서는 음의 에너지 궤도가 존재하므로 회전 에너지를 뽑아낼 수 있고(펜로즈 과정), 최대 추출량은 극단 커 블랙홀에서 총 질량에너지의 약 29%다. 실제 천체물리에서 중요한 것은 이 역학적 버전이 아니라 자기장을 통한 전자기적 추출(블랜포드-즈나예크 기구, 1977)이며, 뒤에 나올 GRMHD 시뮬레이션이 재현하려는 제트의 동력원이 바로 이것이다.

무모정리(no-hair theorem)는 “정상 상태의 점근적 평탄 블랙홀은 질량·각운동량·전하 세 개로 완전히 결정된다”는 주장이다. 정확히는 하나의 정리가 아니라 이스라엘·카터·로빈슨·호킹의 유일성 정리 묶음이고, 각각 해석성·정상성·전기진공 같은 가정을 깔고 있다. 그래서 “털이 없음이 증명됐다”보다는 “넓은 조건에서 유일성이 증명됐다”가 정직한 서술이다. 계산 쪽 함의는 뚜렷하다 — 병합 후 남는 잔해가 커 해 하나로 수렴하므로, 감쇠 진동(ringdown)의 주파수와 감쇠율이 오직 최종 질량과 스핀으로만 결정된다. 여러 개의 준정상 모드를 독립적으로 측정해 같은 $(M,a)$ 를 가리키는지 보는 것이 이른바 블랙홀 분광학이고, 무모정리의 관측적 시험이다.

5. 수치상대론 — 3+1 분해[편집]

아인슈타인 방정식 $G_{\mu\nu}=8\pi T_{\mu\nu}$ 는 시공간 전체에 대한 관계식이라 그대로는 컴퓨터에 넣을 수 없다. 초기값 문제로 바꾸는 절차가 ADM 3+1 분해(아노윗-데저-미즈너, 1959~62)다. 시공간을 공간꼴 초곡면들로 얇게 썰고

ds^2 = -\alpha^2 dt^2 + \gamma_{ij}\left(dx^i+\beta^i dt\right)\left(dx^j+\beta^j dt\right)

로 적는다. $\alpha$ 는 랩스(고유시간이 얼마나 흐르는가), $\beta^i$ 는 시프트(격자점이 옆으로 얼마나 미끄러지는가)이며, 둘 다 물리가 아니라 좌표 선택이다. 이것이 이 분야의 첫 번째 문화 충격 — 게이지를 잘못 고르면 물리적으로 멀쩡한 시공간에서 코드가 폭발한다.

방정식은 두 갈래로 갈린다.

구속조건 — 시간 미분이 없는 4개. 해밀토니안 구속 $R + K^2 - K_{ij}K^{ij} = 16\pi\rho$ 와 운동량 구속 $D_j(K^{ij}-\gamma^{ij}K)=8\pi S^i$ . 초기 데이터가 반드시 만족해야 하고, 이를 푸는 것 자체가 포아송 방정식 꼴의 타원형 문제다.
발전방정식 — $\gamma_{ij}$ 와 외재곡률 $K_{ij}$ 의 시간 발전 12개.

이론적으로 구속조건은 발전에 의해 보존되지만, 수치적으로는 보존되지 않는다. 반올림 오차가 구속조건 위반 모드를 깨우고, 그 모드가 지수적으로 자란다. 1990년대 내내 이진 블랙홀 코드가 궤도 한 바퀴도 못 돌고 죽은 이유가 이것이다.

돌파구는 방정식을 다시 쓰는 것이었다. 표준 ADM 형식은 약쌍곡형이라 잘 정립된 초기값 문제가 아니다.

BSSN(시바타-나카무라 1995, 바움가르테-샤피로 1999): 계량을 등각 분해하고 자취를 분리한 뒤, 등각 접속함수 $\tilde\Gamma^i$ 를 독립 변수로 승격한다. 이 승격 덕에 계량의 공간 2계 미분이 라플라시안 꼴로 정리되어 시스템이 강쌍곡형이 된다. 유한차분 + 적응 격자 세분화와 궁합이 좋다.
일반화 조화 형식: 좌표 자체에 파동방정식 $\Box x^\mu = H^\mu$ 를 부과하면 장방정식이 각 성분마다 파동방정식 꼴이 된다. 여기에 구속조건 감쇠항을 넣어 위반 모드를 지수적으로 끌어내리는 것이 결정타였다.
Z4c / CCZ4: 구속조건 위반을 아예 전파하는 장으로 승격해 감쇠시키는 계열. 요즘 중성자별 코드에서 많이 쓴다.

특이점 처리도 두 갈래다. 엑시전은 지평선 안쪽 영역을 격자에서 통째로 도려낸다(내부는 인과적으로 바깥에 영향을 못 주므로 정당하다). 펑처는 계량을 해석적으로 알려진 특이 부분과 매끄러운 부분으로 쪼개고 특이 부분을 고정한 채 나머지만 발전시킨다. 초기에는 펑처를 격자에 고정해 놨는데, 2005년의 발견은 1+log 슬라이싱과 감마-드라이버 시프트를 쓰면 펑처가 격자 위를 그냥 돌아다녀도 된다는 것이었다. 게이지가 슬라이스를 특이점에서 알아서 밀어내기 때문이다.

6. 2005년, 그리고 파형 카탈로그[편집]

2005년 프란스 프레토리우스가 일반화 조화 형식 + 구속조건 감쇠 + 엑시전으로 이진 블랙홀의 궤도·병합·감쇠진동을 처음으로 끝까지 적분했다. 몇 달 뒤 브라운스빌 그룹과 NASA 고다드 그룹이 독립적으로 움직이는 펑처 기법을 발표하며 같은 결과를 훨씬 단순한 코드로 재현했고, 그 순간 분야가 열렸다. 40년 묵은 문제가 1년 사이에 두 가지 방법으로 뚫린, 과학사에서도 드문 장면이다.

지금의 표준 파이프라인은 이렇다.

단계	하는 일	함정
초기 데이터	구속조건 타원 방정식을 풀어 궤도 상태 구성	등각평탄 가정이 남기는 정크 복사
발전	BSSN/GHG + 고차 유한차분 또는 다영역 유사스펙트럼	게이지 선택, 구속조건 표류
지평선 진단	겉보기 지평선 탐색, 스핀·질량 추출	병합 직후 탐색 실패
파형 추출	뉴먼-펜로즈 $\Psi_4$ 를 스핀가중 구면조화 함수로 분해	유한 반경 추출 → 널 무한대 외삽
후처리	$\Psi_4$ 를 두 번 시간적분해 $h_+,h_\times$ 로	적분 상수 표류(고정주파수 적분으로 억제)

계산량은 잔인하다. 질량비가 크거나 스핀이 크면 궤도 수가 늘어 한 파형에 수십만~수백만 코어시간이 든다. 그래서 파형을 매번 시뮬레이션할 수 없다는 것이 LIGO 데이터 분석의 출발점이다. 실제 탐색은 정합 필터에 템플릿 뱅크 — 질량·스핀 공간을 촘촘히 덮는 수십만 개의 파형 — 를 물려 돌린다. 이 뱅크를 채우는 것은 수치상대론 파형 자체가 아니라, 그것에 교정된 준해석 모형이다.

EOB(유효 일체 접근): 이체 문제를 유효 시험입자 운동으로 사상하고, 미지 계수를 수치상대론 파형에 맞춘다. 병합 전 궤도 진화를 싸고 정확하게 준다.
IMRPhenom 계열: 주파수 영역에서 진폭·위상을 조각별 함수로 피팅한 현상론 모형. 필터 하나 돌리는 비용이 극히 싸서 대량 탐색용.
NR 대리모형(surrogate): 수치 파형 집합에 직접 축소기저·보간을 적용. 교정 영역 안에서는 가장 정확하다.

뱅크 설계는 순수한 수치해석 문제다. 파형 공간에 계량을 정의하고(두 파형의 내적이 곧 거리), 최소 정합(보통 0.97) 이상을 보장하도록 격자 또는 확률적 배치로 덮는다. 정합이 0.97이면 신호대잡음비 3% 손실, 즉 탐지 가능 부피 약 10% 손실. “수치상대론 → 준해석 모형 → 템플릿 뱅크 → 정합 필터”가 2015년 첫 검출로 이어진 사슬 전체이고, 이 사슬 어디를 끊어도 사건은 잡히지 않았다.

7. 강착과 영상 — GRMHD와 널 측지선 적분[편집]

블랙홀은 그 자체로는 안 보인다. 보이는 것은 주변 플라스마이므로, 영상 예측은 두 단계 계산이다.

1단계, GRMHD. 커 시공간(고정 배경, 원반 질량은 무시)에서 이상 자기유체 방정식을 보존형으로 푼다. 자기유체역학의 표준 도구에 상대론적 보정이 붙는 셈이지만, 실무적 고통 지점은 따로 있다.

보존량 → 원시변수 역변환. 상대론에서는 $(\rho, u^i, p)$ 를 보존량에서 얻는 것이 닫힌 형태가 아니라 비선형 방정식 풀이다. 강자기장·저밀도 영역에서 이 뉴턴 반복이 실패하는 것이 코드가 죽는 1순위 원인.
바닥값(floor). 제트 깔때기 영역은 밀도가 사실상 0으로 가는데, 수치적으로 음수 밀도가 나오면 끝이다. 최소 밀도·내부에너지를 강제로 깔아 주는데, 이 인공 질량 주입이 제트 동역학에 실제로 영향을 준다는 것이 이 분야가 인정하는 체계 오차다.
$\nabla\cdot\mathbf B=0$ 유지를 위한 구속 수송 또는 발산 세척.

강착 상태는 크게 SANE(자기장이 난류에 눌린 표준 원반)과 MAD(자기 플럭스가 포화해 강착을 저지하는 상태)로 갈리고, EHT의 제트 동력 논의는 대체로 MAD 쪽 손을 들어 준다.

2단계, 일반상대론적 복사전달(GRRT). 카메라 화소마다 광선을 시간 역방향으로 쏘아 널 측지선을 적분한다.

\frac{d^2x^\mu}{d\lambda^2} + \Gamma^\mu_{\ \alpha\beta}\frac{dx^\alpha}{d\lambda}\frac{dx^\beta}{d\lambda}=0

를 룽게-쿠타법으로 적분하면서, 그 위에서 편광 복사전달 방정식을 함께 적분한다. 핵심 수치 트릭은 불변량으로 적분하는 것 — $I_\nu/\nu^3$ 과 $\nu\,\alpha_\nu$ 는 로런츠 불변이라, 적색편이를 매 단계 손으로 보정하는 대신 불변량을 나르고 마지막에 관측자 진동수로 환산하면 된다. 커 시공간에서는 카터 상수 덕분에 측지선이 완전 적분가능해서 구적법으로도 풀 수 있지만, 일반 배경·수치 시공간에는 그냥 상미분방정식 적분기를 쓴다.²

여기서 레이 트레이싱과의 관계가 분명해진다. 알고리즘의 뼈대는 컴퓨터그래픽스의 역방향 광선추적과 같다. 다른 점은 광선이 직선이 아니라 곡률에 따라 휘는 측지선이고, 매질이 상대론적으로 움직이며, 적분해야 할 방정식이 4개가 아니라 8개(위치 4 + 운동량 4)라는 것뿐이다.

마지막 단계는 비교다. GRMHD 파라미터(스핀, 자기 상태, 전자 온도 처방, 관측 각도)를 훑어 수만 장의 합성 영상을 만들고, 각각을 전파간섭계의 실제 $(u,v)$ 표본에 투영해 관측 데이터와 점수를 매긴다. 2019년의 M87* 영상과 2022년의 궁수자리 A* 영상은 그렇게 걸러진 결과다. 즉 EHT의 “사진”은 관측이자 동시에 거대한 시뮬레이션 카탈로그와의 적합이며, 이 사실을 빼고 읽으면 그림의 의미를 오해하게 된다.

8. 호킹 복사[편집]

1974년 호킹은 휘어진 시공간의 양자장론에서 블랙홀이 정확히 열복사를 방출한다는 것을 보였다. 온도는

T_H = \frac{\hbar c^3}{8\pi G M k_B}

이고 엔트로피는 지평선 넓이에 비례한다( $S = k_B c^3 A / 4G\hbar$ ). 핵심은 $T \propto 1/M$ , 즉 비열이 음수라는 것이다. 방출할수록 가벼워지고, 가벼워질수록 뜨거워지며, 마지막에는 폭주한다.

수치를 넣어 보면 이 효과가 왜 관측되지 않는지 즉시 보인다. 태양질량 블랙홀의 온도는 약 $6\times10^{-8}$ K로, 2.7 K인 우주배경복사보다 압도적으로 차갑다. 먹는 게 뱉는 것보다 많으므로 실제로는 계속 자란다. 증발까지 걸리는 시간은 $\propto M^3$ 이라 태양질량이면 $10^{67}$ 년 규모. 우주 나이의 $10^{57}$ 배쯤 된다.

계산과학 쪽 파장은 정보 역설이다. 순수상태에서 출발한 물질이 완전한 열복사가 되어 나온다면 유니터리성이 깨지고, 이는 양자역학과 일반상대론 중 하나가 틀렸다는 뜻이다. AdS-CFT 대응이 이 논쟁에서 특별한 위치를 갖는 이유는, 경계 이론이 명백히 유니터리한데 벌크에 블랙홀이 있는 상황을 제공하기 때문이다. 최근의 얽힘 엔트로피 계산(페이지 곡선의 재현)은 이 방향에서 나온 결과다. 다만 아직 아무도 “블랙홀에서 정보가 어떻게 나오는가”를 구체적 기구로 답하지는 못했다.³

9. 여담[편집]

슈바르츠실트는 제1차 세계대전 동부전선 복무 중에 이 해를 찾아 아인슈타인에게 보냈고, 몇 달 뒤 병으로 사망했다. 정작 아인슈타인 본인은 지평선이 물리적으로 실현될 수 없다고 오래 믿었다.
“블랙홀”이라는 이름은 존 휠러가 1967년 강연에서 쓰면서 굳었다. 그 전에는 “얼어붙은 별(frozen star)“이라고도 불렀는데, 이건 멀리 있는 관측자 기준으로 붕괴가 지평선 근처에서 무한히 느려 보이기 때문이다. 좌표 선택이 작명까지 바꾼 사례.
영화 인터스텔라의 블랙홀 영상은 실제 커 시공간의 널 측지선을 적분해 만든 것이고, 화소당 광선 하나가 아니라 광선 다발을 나르는 방식으로 노이즈를 줄였다. 제작 과정이 학술지 논문으로 정리되어 나왔다.
이 분야의 오픈소스 생태계는 의외로 튼튼하다. 아인슈타인 툴킷, SpEC/SpECTRE, GRChombo, 그리고 GRMHD 쪽 HARM 계열이 공개돼 있어서, GPU 컴퓨팅 자원만 있으면 학생도 이진 병합을 돌려 볼 수 있다. 물론 “돌려 볼 수 있다”와 “수렴한 파형을 얻는다”는 다른 이야기다.

10. 관련 문서[편집]

마름모 그림에 화이트홀과 “다른 우주”가 등장하니 흥분하기 쉬운데, 그건 영원히 존재해 온 진공 블랙홀의 최대 확장일 뿐이다. 별이 붕괴해 실제로 만들어진 블랙홀은 붕괴 이전 시공간이 물질로 채워져 있어 그 영역이 존재하지 않는다. 즉 웜홀로 옆 우주에 가는 계획은 그림을 잘못 읽은 데서 출발한다. ↩
카터 상수(1968)는 커 측지선을 완전 적분가능하게 만드는 네 번째 보존량인데, 대칭성에서 직접 나오는 것이 아니라 킬링 텐서에서 나온다. “숨은 대칭”이라 불리는 이유다. 덕분에 커 배경의 광선추적은 원리적으로 4차원 ODE 적분이 아니라 1차원 구적법으로 환원되고, 정밀도 검증용 기준해로 애용된다. ↩
그래서 이 분야 세미나의 질의응답은 대체로 “그래서 정보가 어떻게 나온다는 겁니까”에서 끝난다. 유니터리성이 회복된다는 계산은 있고, 그 계산이 무엇을 의미하는지에 대한 합의는 없다. 물리학이 40년째 같은 질문을 붙들고 있는 흔치 않은 사례이고, 그동안 배운 것은 적어도 “블랙홀은 열역학계다”라는 사실이 우연이 아니라는 점이다. ↩

블랙홀