Title: The structured pseudospectrum of Schwarzschild black hole

Motivation

QNM 문제는 지평에서 들어가는 파, 무한원에서 나가는 파라는 경계조건 때문에 비에르미트(non-Hermitian)이고, 그래서 작은 섭동에도 스펙트럼이 크게 움직일 수 있다. Jaramillo–Panosso Macedo–Al Sheikh(2021) 이후 pseudospectrum(주어진 노름에서 크기 $\varepsilon$ 이하의 섭동으로 도달 가능한 진동수 영역)이 이 불안정성의 표준 진단 도구가 됐고, 특히 overtone이 극도로 불안정하다는 결론이 널리 인용된다.

문제는 기존(비구조) pseudospectrum이 임의의 복소 연산자 섭동을 허용한다는 점이다. 그중 상당수는 실제 물리적 변형(물질 분포, 수정 중력 등)에서 나올 수 없는 방향이다. 실제로 퍼텐셜 bump가 기본 모드까지 불안정하게 만든다는 결과(Cheung et al. 2022)에 대해, 그런 bump가 물리적으로 허용되지 않는 물질 배치에 해당한다는 반론(Cardoso, Kastha, Panosso Macedo 2024; Della Rocca, Rosato, Pani 2026)이 나왔다. 그렇다면 물리적으로 허용된 섭동만 남겼을 때 불안정성 중 얼마가 남는가.

Contribution

• 블랙홀 섭동 이론에 구조화된 pseudospectrum(structured pseudospectrum)을 처음 도입(저자 주장).
• 실수 퍼텐셜 변형 $\delta V \in \mathcal{S} = \mathrm{span}\{\phi_1,\dots,\phi_m\}$으로 섭동을 제한하는 두 공간을 구성:
  • $\mathcal{S}_{\rm EFT}$: Riemann 3차·4차 곡률 EFT 보정에서 나오는 $\phi_j = (1-\sigma)\sigma^j$, $j\ge10$ (축 중력 섭동, $s=2$).
  • $\mathcal{S}_{\rm RZ}$: 지평과 거북이 좌표를 유지하는 RZ 계량 부분공간에서 유도 (무질량 스칼라, $s=0$).
• 에너지 노름과 Hilbert–Schmidt 기반 노름 두 가지로 구조화 거리 $d_{\rm str,E}$, $d_{\rm str,\mathcal{S}}$를 정의하고 부등식 $d_{\rm unstr}\le d_{\rm str,E}\le d_{\rm str,\mathcal{S}}\le\sqrt{N_c}\,d_{\rm str,E}$를 증명.
• 작은 $\varepsilon$에서 $d_{\rm str,\mathcal{S}}$ 등고선이 점근적으로 타원이 된다는 정리(부록 D)와 수치 확인.

Novelty

비구조 pseudospectrum은 고립된 단순 QNM 주변에서 작은 $\varepsilon$일 때 원형 등고선을 갖는다. 즉 스펙트럼이 어느 방향으로든 같은 비용으로 움직인다. 이 논문은 섭동을 물리적 공간으로 제한하면

1. 도달 가능한 이동 범위 자체가 수 자릿수 줄고,
2. 선호 방향이 생긴다(이동이 특정 방향으로 쉬움)

는 것을 보였다. 두 번째가 질적으로 새로운 점이다. 수정 중력이나 환경 효과가 QNM을 어느 방향으로 옮기는지에 대한 모델 의존적 정보를 pseudospectrum 틀 안에서 다룰 수 있게 된다.

Method

• 배경: 슈바르츠실트, 콤팩트 hyperboloidal 좌표($\sigma\in[0,1]$, $\sigma=0$이 $\mathscr{I}^+$, $\sigma=1$이 지평)에서 1계 고유값 문제 $(\omega I - L)U = 0$.
• 노름: 섭동 방정식에서 나오는 에너지 노름. Chebyshev–Lobatto 이산화 후 Cholesky 분해로 2-노름 문제로 변환.
• 비구조 거리: $d_{\rm unstr}(\omega) = \sigma_{\min}[\mathbf{W}(\omega\mathbf{I}-\mathbf{L})\mathbf{W}^{-1}]$.
• 구조 거리: $\delta V = \sum_k c_k\phi_k$ (실수 $c_k$)로 제한한 뒤 최소 노름 문제를 최적화로 풂(부록 B).
• 수치 설정: $N_c = 120$, 주파수 격자 $150\times150$, 섭동 공간 차원 $m=4$.

Results

• 전역적 억제: 복소 진동수 평면 대부분에서 구조 거리가 비구조 거리보다 훨씬 크다. 즉 같은 $\varepsilon$에서 구조 pseudospectrum 영역이 훨씬 작다.
• 기본 모드 $n=0$ ($|\omega-\omega_0|\sim10^{-2}$ 이동에 필요한 섭동 크기):
  • EFT: 비구조 $\sim10^{-6}$ → 구조 $\sim10^{-2}$
  • RZ: 비구조 $\sim10^{-5}$ → 구조 $\sim10^{-2}$
• 첫 overtone $n=1$: 비구조 등고선이 $10^{-12}$(EFT), $10^{-12}$–$10^{-11}$(RZ) 수준인데 구조 등고선은 $\sim10^{-2}$에서도 닫혀 있다. 즉 기본 모드와 overtone 사이의 불안정성 위계가 크게 줄어든다.
• 두 구조 노름의 차이: $d_{\rm str,\mathcal{S}}/d_{\rm str,E}$는 약 $1$–$2.2$로 이론 상한($\sqrt{120}\approx10.95$)보다 훨씬 작다. 노름 선택은 등고선 모양을 조금 바꿀 뿐 전체적인 억제는 같다.
• 비등방성: EFT에서 $n=0$ 등고선은 비스듬한 방향으로 늘어나고 $n=1$은 다른 방향·종횡비를 갖는다. RZ에서 $n=1$은 주로 허수축(감쇠율) 방향으로 늘어난다.

Significance

• overtone 불안정성 논쟁의 재해석: "작은 섭동으로 overtone이 크게 움직인다"는 결론이 섭동 공간 선택에 크게 의존한다는 것을 정량적으로 보였다. 링다운 분광학에서 overtone을 쓰는 것의 정당성 논의에 직접 연결된다.
• 모델 의존적 안정성 평가의 틀: "어떤 섭동 공간과 어떤 노름에서의 안정성인가"를 명시하도록 만든다. 수정 중력이나 환경 모델마다 자기 섭동 공간을 넣어 쓸 수 있는 일반적인 도구다.
• 방향 정보: 타원 등고선의 방향은 특정 이론군이 QNM을 주로 진동수 쪽으로 옮기는지, 감쇠율 쪽으로 옮기는지를 알려준다. 분광학 검정에서 어느 측정량이 민감한지 판단하는 데 쓸 수 있다.

Limitations

• 섭동 공간이 매우 작다: 결과는 $m=4$차원 공간에 대한 것이다. 제약이 강할수록 불안정성이 줄어드는 것은 정의상 당연한 방향이고(부등식 (36)), $m$을 키우면 억제가 얼마나 유지되는지는 보이지 않았다. 이 점이 결론의 강도를 판단하는 핵심인데 논문에는 $m$ 의존성 분석이 없다.
• 공간 선택의 자의성: EFT 공간에서 $j=10$만 특정 곡률 연산자에 대응하고, $j>10$ 기저는 저자도 "특정 연산자와 일대일 대응하지 않는 확장"이라고 명시한다. RZ 공간도 $\epsilon=0$, $B=A$라는 특수 부분공간이다.
• 섭동 유형이 분리됨: EFT는 축 중력 섭동($s=2$), RZ는 스칼라($s=0$)로만 계산해 두 공간을 같은 섭동 유형에서 비교하지 않았다.
• 수치 수렴: 해상도 $N_c=120$ 하나만 보고했다. pseudospectrum은 이산화 해상도에 민감한 것으로 알려져 있어 수렴 확인이 아쉽다.
• 모드 범위: 자세한 분석은 $n=0, 1$에 집중됐다. 더 높은 overtone에서 같은 억제가 유지되는지는 큰 그림(Fig. 1, 2)으로만 보인다.
• 슈바르츠실트만: Kerr 확장과 total transmission mode는 향후 과제로 남겼다.
• pseudospectrum은 고유값 이동 범위의 진단이지 파형(시간영역 링다운)에 미치는 영향은 아니다. 파형 수준의 함의는 따로 따져야 한다.

Takeaway

물리적으로 허용된 퍼텐셜 변형만 남기면 슈바르츠실트 QNM, 특히 overtone의 pseudospectrum 불안정성은 수 자릿수 줄고 방향성을 갖게 되지만, 이 결론은 4차원 섭동 공간이라는 강한 제약 위에 서 있다.