이벤트 호라이즌을 통과하는 널 좌표 시뮬레이션과 전하 스칼라장 적용

초록

저자들은 널 좌표계에서 $-2G,du,(dx+B,du)+R^{2}\dots$ 형태의 계량을 이용해, $R$ 혹은 $G$ 중 하나를 진화시키고 다른 하나는 레이차우디 방정식으로 제약하는 세 가지 ‘포뮬레이션’을 제시한다. 이를 임의의 방사형 게이지와 결합해 사건 지평선을 통과하는 시뮬레이션을 구현하고, 구형 전하 스칼라장을 대상으로 수렴성 검증, 임계 붕괴의 질량·전하 스케일링 지수와 미세구조, 그리고 극한 전하 Reissner‑Nordström 해의 섭동을 분석한다.

상세 분석

이 논문은 널 좌표계에서 중력 붕괴를 다루는 새로운 수치적 프레임워크를 제시한다. 기본 메트릭을 $ds^{2}=-2G,du,(dx+B,du)+R^{2}\gamma_{ij}(d\theta^{i}+\beta^{i}du)(d\theta^{j}+\beta^{j}du)$ 로 두고, $u$를 시간, $x$를 방사형 좌표로 삼는다. 여기서 $G>0$는 광선 생성자의 정규화 인자이며, $B$는 방사형 게이지를 조정하는 자유도이다. 중요한 점은 $R$(면적 반경)과 $G$ 사이의 관계를 레이차우디 방정식(광선 팽창률 $\rho_{+}=d\ln R/d\lambda$)을 통해 제약할 수 있다는 점이다.

세 가지 포뮬레이션은 다음과 같이 정의된다.

1. eR (evolve‑R): $R$을 진화시키고 레이차우디 방정식으로 $G$를 구한다. 이 경우 $R_{,x}>0$이 필요해 사건 지평선 근처에서 $R$이 감소하면 시뮬레이션이 중단된다.
2. eG (evolve‑G): $G$를 진화시키고 $R$을 레이차우디 방정식으로 구한다. $R_{,x}$의 부호 제한이 없으므로 사건 지평선을 통과할 수 있다.
3. fe (free evolution): $R$과 $G$를 모두 진화시키고 레이차우디 방정식은 초기 데이터 구성에만 사용한다.

이 포뮬레이션들은 방사형 게이지 선택($x\to x’(u,x,\theta)$)과 독립적으로 결합될 수 있다. 저자들은 특히 local shifted affine (lsa) 게이지($G=G(u,x)$)와 local shifted Bondi (lsB) 게이지($R=R(u,x)$)를 이용해, 기존 연구에서 발생했던 “$B$가 발산” 문제를 회피한다.

구형 전하 스칼라장 시스템을 테스트베드로 삼아, 두 종류의 계산 영역을 설정한다. 첫 번째는 $R=0$인 정규 중심에서 시작하는 외향 널 원뿔들, 두 번째는 입사 널 경계와 교차하는 “널 직사각형” 영역이다. 각각에서 저자들은 다음을 수행한다.

• 수렴성 검증: 동일한 초기 데이터와 동일한 해상도에서 세 포뮬레이션(eR, eG, fe)을 비교하고, 전반적으로 2차 수렴을 확인한다. 특히 널 직사각형에서는 점별 수렴이 깨끗하게 나타난다.
• 사건 지평선 통과: eG와 fe 포뮬레이션을 이용해, 사건 지평선 내부까지 진화한 뒤 간단한 특이점 배제(singularity excision) 기법을 적용한다. 배제 영역을 $R$이 일정 임계값 이하가 되는 지점으로 정의하고, 해당 영역을 계산에서 제외함으로써 안정적인 장기 진화를 얻는다.
• 임계 붕괴: 한 파라미터($p$)를 조정해 블랙홀 형성 임계점 $p_{}$를 탐색한다. 질량 $M$과 전하 $Q$는 $M\sim (p-p_{})^{\gamma_{M}}$, $Q\sim (p-p_{*})^{\gamma_{Q}}$ 형태의 스케일링을 보이며, $\gamma_{M}\approx0.37$, $\gamma_{Q}\approx0.68$ 정도의 지수를 얻는다. 또한 로그 주기적 미세구조가 존재함을 확인하고, 전하에 대한 미세구조는 질량에 비해 보편성이 낮다는 새로운 물리적 통찰을 제공한다.
• 극한 전하 Reissner‑Nordström 해: 초기 데이터를 약간 변형해 극한 전하(Extremal) RN 해를 만들고, 작은 섭동을 가해 진화시켰다. eG와 fe 포뮬레이션 모두 섭동이 사라지면서 해가 안정적으로 극한 RN 형태로 수렴함을 보여준다.

이 논문은 널 좌표계에서 사건 지평선을 넘는 시뮬레이션이 가능하도록 하는 이론적·수치적 토대를 마련했으며, 특히 전하를 포함한 복소 스칼라장 시스템에서 임계 현상의 새로운 정량적 결과를 제공한다. 포뮬레이션과 게이지의 독립적 조합 가능성은 향후 비구형(예: 축대칭) 임계 붕괴 연구에 중요한 전진을 의미한다.