Title: Quasinormal Mode Excitation in the Horizon Shear of Near-Extremal Kerr Black Holes

Motivation

EMRI가 Kerr 블랙홀 지평을 조석 변형시키는 연구(O'Sullivan & Hughes 2014, 2016; Penna, Hughes, O'Sullivan 2017)에서, O'Sullivan–Hughes(2016)는 근극한 스핀에서 근점 통과 후 지평 shear에 **작고 빠르게 감쇠하는 진동(wiggle)**이 남는 것을 발견했다. 그들은 진동수를 궤도와 연결하지 못했고, 조석 구동장 $\psi_0^{\rm HH}$에서는 wiggle을 찾지 못해 Kerr QNM이 아니라고 결론지었다.

반면 Thornburg, Wardell, van de Meent(2020)는 무한원 EMRI 파형의 비슷한 wiggle이 낮은·중간 스핀에서 QNM 들뜸임을 보였다. 지평 쪽 wiggle의 정체가 남은 질문이다. 근극한에서는 QNM이 실수축 근처로 몰려 수치적으로 구하기 어려웠는데, 최근 방법으로 정확한 값을 얻을 수 있게 됐다.

Contribution

• 5개 근극한 스핀에서 GSN 코드로 지평 shear를 주파수 영역 모드 합($\ell\le9$, $|n|\le150$)으로 계산.
• 전체 시간 영역 신호에 "QNM 2개 + 2차 다항식 배경" 피팅, Mahalanobis 거리로 Kerr QNM과 비교.
• 다중극 $(\ell,m)$별로 신호를 분리해 다시 피팅. 다른 다중극의 오염 없이 $(2,1,0)$, $(3,2,0)$를 확인하고 $(3,1,0)$, $(4,1,0)$, $(4,2,0)$의 흔적을 찾음.
• 전달 함수 $\sigma_{\ell mn} = i\psi_{0,\ell mn}/(p_{mn}-i\kappa)$에 QNM 구조가 없으므로, shear의 wiggle은 구동장에도 반드시 있어야 함을 논증하고 실제로 보임.

Novelty

• 근극한 Kerr QNM은 감쇠가 0으로 가는 zero-damping mode(ZDM)와 유한하게 남는 damped mode(DM) 두 가지로 갈라진다. wiggle이 DM 가지라는 것을 정량 판정한 것이 핵심이다. 스핀이 바뀌어도 감쇠 시간이 거의 같다는 관찰이 그 단서다.
• 진동수 격자 위 진폭에서 극점을 직접 피팅하는 대신, 다중극별 역 Fourier 변환으로 시간 영역 감쇠 사인 피팅 문제로 바꾼 점이 실용적이다. 격자 간격($\Omega_r\approx0.0205$)이 극점 폭과 비슷해 직접 피팅이 불안정했기 때문이다.

Method

• 배경과 궤도: $a/M\in\{0.999, 0.9995, 0.9999, 0.99995, 0.99999\}$, 적도 순행 궤도 $p=3M$, $e=0.85$. 첫 궤도 burst만 분석. 관측 각 $(\theta,\psi)=(\pi/2,0)$에서 cross 편극이 사라져 plus만 사용.
• 기준 QNM: Cook–Zalutskiy 방법(qnm 패키지)으로 DM 진동수 계산. 패키지가 두 가지를 구분하지 못해 $(2,2,0)$, $(3,3,0)$, $(4,3,0)$, $(4,4,0)$, $(2,2,1)$, $(3,2,1)$ 등은 DM 값을 얻지 못했다($m=\ell$은 DM 가지가 원래 없음).
• 피팅: 최소제곱, 기준 창 $v\in[377.5,487.5]M$. 창 끝점($[365,390]$, $[475,500]$)과 다항식 차수를 바꿔 체계 오차를 공분산에 포함.
• 판정: 2차원 Mahalanobis 거리 $D_2^2$를 $\chi^2_2$ 분위수로 $n\varepsilon$에 대응($D^2=2.296, 6.180, 11.829$ ↔ $1,2,3\varepsilon$). 잡음 없는 자료라 명목적 척도임을 저자도 밝힌다.

Results

• 전체 신호 피팅: 첫 모드는 모든 스핀에서 $(2,1,0)$과 $3\varepsilon$ 안($D^2=5.4$–$6.6$), 둘째 모드는 $(3,2,0)$에 가장 가깝다. 다음 후보와의 거리는 $D^2\sim10^2$–$10^3$으로 배정은 분명하다. 다만 $a/M\ge0.9999$에서 둘째 모드가 Kerr 값에서 벗어난다($D^2=15.7$, $28.1$, $38.6$).
• 다중극 분리 피팅: $(2,1)$, $(3,2)$는 모든 스핀에서 기본 모드와 $2\varepsilon$ 안이며 대부분 $1\varepsilon$ 안이다. $(3,1)$, $(4,1)$, $(4,2)$도 $2\varepsilon$ 안이지만, 추정 진폭이 2-모드 잔차 이하라 전체 신호에서는 분해되지 않는다.
• 편향의 원인: 전체 신호에서 $(2,1,0)$ 감쇠율이 더 크게 나오는 것은 초기에 덜 감쇠한 다른 다중극이 섞인 탓으로 해석.
• $(3,2,0)$ 표류: 허수부가 ZDM 쪽으로 움직이는 듯 보이지만, 실수부가 $2\Omega_H$로 가지 않고 ZDM을 고정 모드로 넣으면 진폭이 불안정해 ZDM은 아니라고 판단. 5-모드 피팅에서는 Kerr 값 쪽으로 돌아온다.
• 구동장의 wiggle: $\psi_{0,+}^{\rm HH}$에도 wiggle이 보이며, 전달 함수 때문에 shear보다 진폭이 작아 이전에 놓쳤다고 설명.
• $m=0$ 미들뜸: $\ell\le4$의 $m=0$ 다중극에서는 DM이 존재할 수 있는데도 들뜬 흔적이 없다. 저자도 설명하지 못한다.

Significance

• 지평 쪽 QNM 들뜸의 확인: 무한원 파형의 QNM 들뜸(Thornburg et al. 2020)과 같은 현상이 지평 shear에도 나타난다. 지평 동역학을 무한원 신호와 같은 언어로 읽을 수 있다는 그림을 지지하며, 지평 shear와 news의 유사성을 다룬 학위논문 2611.01275와 같은 방향이다.
• 근극한 스펙트럼 검증: 두 가지로 갈라지는 근극한 QNM 중 어느 쪽이 실제 원천에 의해 들뜨는지를 보여 주는 드문 시간 영역 사례다.
• 방법론: 다중극 분리와 Mahalanobis 판정은 잡음 없는 섭동 자료에서 QNM을 식별하는 재사용 가능한 절차다.

Limitations

• 단일 궤도 구성: $p=3M$, $e=0.85$, 적도 순행 하나만 분석했다. 저자는 결과가 $p$, $e$에 매끄럽게 의존한다고 쓰지만 다른 구성의 결과는 보여 주지 않는다. 경사 궤도나 역행 궤도는 없다.
• 배경 모형이 임의적: EMRI 원천 배경의 함수 형태를 몰라 2차 다항식으로 근사했고, 그림 3의 재구성은 데이터와 눈에 띄게 어긋난다. 저자는 이 차이를 배경 탓으로 돌리지만 독립 검증은 없다.
• 기준값의 한계: qnm 패키지가 ZDM과 DM을 구분하지 못해 일부 모드는 기준값이 없다(특히 $(3,2,1)$). $(3,2,0)$ 표류의 원인으로 저자는 e-fold 수 부족(~2)과 qnm·GSN 정확도를 꼽지만 확정하지 않았다.
• 통계 해석: 오차 막대는 창·배경 변화로 만든 체계 오차가 대부분이므로, "$2\varepsilon$ 안"은 확률적 의미보다 명목적 일치도로 읽어야 한다.
• 약한 모드는 '흔적' 수준: $(3,1,0)$, $(4,1,0)$, $(4,2,0)$은 다중극 분리 신호에서만 보이고, 전체 신호에서는 잔차에 묻힌다. 초록의 "evidence"는 그 정도 의미다. $\ell\ge5$는 분석하지 않았다.
• 관측과의 거리: 지평 shear는 관측량이 아니며, LISA EMRI 파형에 대한 함의는 논의하지 않는다.

Takeaway

근극한 Kerr 지평 shear의 근점 후 wiggle은 궤도와 무관한 damped-mode Kerr QNM(주로 $(2,1,0)$, $(3,2,0)$)이며, 조석 구동장에도 같은 성분이 더 작은 진폭으로 들어 있다.