Title: Gravitational radiation from a photon on the light ring: multipole fluxes and the logarithmic divergence

Motivation

light ring 근처 원궤도의 입자가 내는 GW 플럭스는 궤도가 light ring에 다가갈수록 발산하지만 입자의 에너지는 유한하다. 1970년대 'geodesic synchrotron radiation' 연구는 질량 있는 입자의 스펙트럼이 $\omega^{-1}$로 가다가 $\omega_{\rm crit}\propto\gamma^2$에서 지수적으로 잘린다는 것을 보였지만, 고정 에너지에서 정지 질량을 0으로 보내는 극한은 취하지 않았다.

저자가 참여한 Barausse et al.(2021)은 $\ell\sim400$까지 수치로 풀어 다중극당 플럭스가 $\kappa E^2/(M^2\ell)$, $\kappa\simeq0.064$로 간다는 것을 보였고 원인을 상대론적 beaming으로 추정했다. 형식적 증명, $\kappa$의 닫힌 형태, 정확한 원인은 남아 있었다.

Contribution

• 광자를 직접 다룸: 질량 있는 입자의 $\delta\to0$ 극한 대신 null 응력-에너지 텐서를 Zerilli–Regge–Wheeler와 Teukolsky 방정식의 원천으로 넣음.
• $\kappa$의 닫힌 형태: $\kappa=2\sum_j\kappa_j=0.063653$($m=\ell-j$; $\kappa_0=0.0277668$, $\kappa_1=0.0037670$, …, 대략 $e^{-\pi j}$로 감소).
• 발산의 두 원인: (i) 원천 진동수 $m\Omega$가 eikonal QNM에서 $\ell$과 무관하게 $(j+\tfrac12)\lambda_{\rm L}$만큼 떨어져 있어 어떤 다중극도 지수적으로 억제되지 않는다. (ii) $T_{tt}$와 $T_{\phi\phi}$에서 오는 원천의 선행 두 항이 $p^\mu p_\mu=0$ 때문에 정확히 상쇄돼, 다중극당 플럭스가 $\ell$로 커지는 대신 $1/\ell$로 준다.
• 무한원 = 지평: 파동이 장벽 정점에서 생기고 정점에서 보면 장벽이 양쪽으로 대칭이라 두 플럭스가 같다. 장벽은 $m=\ell$ 파동의 $96\%$를 반사하므로 "절반이 투과한다"는 뜻은 아니다. $O(\ell^{-1/2})$ 비대칭 계수 $\bar\sigma=0.61$도 계산.
• Kerr: 모든 적도 light ring에서 $\Omega_\theta=\lambda_{\rm L}$이라 $\eta_j=j+\tfrac12$가 스핀에 무관하고, 스핀 의존성 전체가 하나의 인자 $g(a)$로 줄어든다. 지평/무한원 비가 모든 스핀에서 정확히 1임을 증명.
• 1973년 공식의 오류: Breuer et al.의 홀수 패리티 항이 4배 작다는 것을 찾음.

Novelty

발산의 원인을 beaming 하나로 보던 이전 설명을 고쳤다. 근공명만으로는 다중극당 플럭스가 오히려 $\ell$에 비례해 커지며, $1/\ell$로 만드는 것은 질량 없는 원천에 특유한 상쇄다. 질량 있는 입자에서 정지 질량이 하는 일은 장벽 매개변수 $\eta$를 $\tfrac12+\tfrac32\ell\delta$로 옮기는 것뿐이고, 이 이동이 $\ell\sim\gamma^2$에서 지수 절단을 켠다.

Kerr 스핀 인자를 닫힌 형태로 얻은 것도 새롭다. Chrzanowski–Misner(1974)의 WKB 결과는 스핀 의존성은 정확히 같지만 상수가 0.94배인데, 저자는 이것이 두 오류(텐서 공식의 $(s!)^2$ 인자와 Stirling 근사)가 서로 상쇄된 우연이라고 분석한다.

Method

• 원천: null 입자의 $T^{\mu\nu}$는 질량 있는 입자의 형태에서 $\mu u^\alpha\to p^\alpha$로 바꾼 것과 같아, 표준 점입자 원천을 $E_0\to E$, $L_0\to bE$로 재사용($b=3\sqrt3M$). Teukolsky 원천으로도 같은 플럭스를 $\ell\le5$에서 9–17자리까지 재현.
• 장벽 정점 해: $m=\ell-j$로 두고 $1/\ell$ 전개. 정점 근처 포물선 장벽에서 Weber 방정식의 포물선 원통 함수 해를 쓰고, 정점 값 $u(0)$와 $\partial_{r_*}u(0)$만으로 플럭스를 구함. 삼차 항은 $O(\ell^{-1/2})$ 비대칭으로 따로 계산.
• 수치 검증: 9차 Runge–Kutta로 동차 방정식을 근사 없이 적분해 $\ell=12800$까지(Schwarzschild), Kerr는 8차 Dormand–Prince Teukolsky 적분으로 $\ell=800$까지. 코드·노트북·데이터 공개.

Results

• Schwarzschild: $400\le\ell\le12800$의 평균 플럭스 fit이 단일 $m$ 합 $0.0318263$을 주어 해석값 $0.0318265$와 맞는다. $\ell=12800$에서 $\kappa_{j\le4}$를 $10^{-3}$ 이내로 재현. 기존 수치 fit $\kappa=0.064\pm0.001$과 일치(Fig. 1).
• 수렴 속도: $\ell\simeq100$에서도 두 플럭스는 약 $12\%$ 다르다. $\ell_{\max}$까지 합했을 때 지평으로 간 비율은 로그적으로만 $1/2$에 다가가며, $\ell_{\max}=100$–$140$에서 $40\%$다. 이것이 $\ell\simeq140$에서 자른 Gundlach et al.(2012)이 $42\%$를 얻은 이유라고 설명한다.
• Kerr: $g(a)$는 $a/M\simeq0.834$($r_0=\sqrt3M$)에서 최대 $6\sqrt3-9\simeq1.392$, 극한 순행에서 0, 극한 역행($r_0=4M$)에서 $81/112\simeq0.723$. 수치 Richardson 외삽이 다섯 스핀 모두에서 해석값과 $1.5\times10^{-5}$ 이내로 맞는다(Fig. 2). 순행 고스핀에서는 장벽이 평평해져 점근 영역이 $\ell\gg800$으로 밀려난다.
• 1973년 공식: 짝수 공식은 맞지만 홀수 공식으로 정확한 플럭스를 나누면 비가 $m=40$–$400$에서 $3.15\to3.77$로 가고 fit 극한이 $4.06$이다. 1973년 홀수 항을 쓰면 $\kappa=0.0580$이 되어 수치값과 맞지 않는다.
• 질량 있는 입자: 식 (24)가 $\gamma=20$에서 $1\%$ 이내로 수치와 맞는다. 절단 $\ell\sim\gamma^2$까지 합하면 각 플럭스가 $\simeq0.127\ln\gamma$로, Barausse et al.의 fit 기울기 $0.12\pm0.01$과 일치.

Significance

• 독립적 이중 확인: 2610.20708(Trestini, Gonzo, Leather)은 질량 있는 입자로 같은 문제를 풀었고, 그 논문 스스로 두 결과가 $\alpha=\kappa/3$으로 일치한다고 밝힌다. 이 논문에는 부선행 계수 $\beta$가 없고, 그 논문에는 광자의 직접 처리와 Kerr 유도가 없다. 두 편을 함께 읽어야 그림이 완성된다.
• Kerr 인자 비교: 2610.20708은 Kerr를 Chrzanowski–Misner 인자 $\Phi(a)$로 추측했고 모드별 수치 비로만 시험했다. 두 논문의 식을 대조하면 $\Phi=g\cdot2\sqrt{3Mr_0}/(r_0+3M)$로, Schwarzschild($r_0=3M$)에서 1이 되는 인자만큼 다르다(극한 역행에서 $81/112\times0.9897\simeq0.7158$로 그 논문의 $\Phi(-1)$과 맞는다). 고정 광자 에너지(이 논문)와 고정 $m\delta(a)$(그 논문)라는 정규화 차이로 보이며, 이 논문의 유도가 그 추측의 스핀 의존성을 뒷받침한다. 두 논문 모두 이 변환을 명시하지는 않는다.
• 응용: 임계 충격 매개변수 근처 초상대론적 산란의 궤도당 에너지 손실 $\Delta E\simeq4.2(E^2/M)\ln\ell_{\max}$(절반은 흡수), light ring 근처 느린 모드 합의 외삽, eikonal 2차 QNM 계산의 입력.

Limitations

• light ring 위 광자 하나: 주된 결과는 정확히 원형 null 궤도에 대한 것이다. 질량 있는 입자로의 확장(Sec. VI)은 선행 차수 공식과 수치 비교이며, 부선행 $\beta$ 계수는 다루지 않는다.
• Kerr의 부선행 항: 지평/무한원 플럭스의 $O(\ell^{-1/2})$ 차이는 Kerr에서 해석적으로 구하지 않았고 $\ell=800$ 수치값($a/M=0.5$, $0.9$, $-0.5$, $-0.9$에서 $0.78$, $1.10$, $0.51$, $0.45$)만 있다. $|a|/M=0.99$는 수치로 돌리지 않았다.
• 1973년 오류의 위치는 미확정: 홀수 패리티의 4배 차이를 수치로 확정했지만, 1973년 유도의 어느 단계에서 생겼는지는 추적하지 않았다고 밝힌다.
• 물리적 절단은 정성적: 발산은 원천의 유한 크기가 자른다는 설명(실제 광자는 파장 $\ell_{\max}\sim Er_0/\hbar$)은 휴리스틱이다. M87* 예시의 에너지 손실 $\Delta E/E\sim10^{-77}$이 보여 주듯 천체물리적 광자에는 영향이 없다.
• 관측과의 거리: 결과는 섭동론 내부의 정밀한 점근식이다. 직접 관측 가능한 양이 아니라 산란·self-force 계산의 입력으로 쓰인다.

Takeaway

light ring 위 광자의 GW 플럭스는 다중극마다 $1/\ell$로 줄어 로그 발산하며, 이는 QNM과의 근공명과 null 조건에 의한 원천 상쇄가 함께 만든 결과다. 계수 $\kappa=0.063653$과 Kerr 스핀 인자 $g(a)$는 닫힌 형태로 구해졌다.