Title: Constraining the Quadratic-mode Amplitude Coupling in GW250114

Motivation

ringdown 분석은 선형 모드 식별을 넘어 2차(quadratic) QNM 탐색으로 넘어가고 있다. GW250114에서는 LVK가 $221$ overtone의 결정적 증거와 $|m|=\ell=4$ 다중극 기여를 보고했고, 이후 연구들은 그 $\ell=4$ 기여가 2차 모드일 가능성을 보았다(ringdown만 쓴 분석에서 조건부 Bayes factor $\simeq10$, IMR 사전분포와 2차 모드 10개를 묶은 분석에서 $62$).

지금까지의 탐색은 "2차 모드의 진동수는 부모 진동수의 합"이라는 관계로 후보를 찾았다. 그런데 이 관계는 단순한 비선형 파동 모형에도 있고 섭동적 수정 중력에서도 유지될 것으로 보인다. 반면 2차 모드 진폭이 부모 진폭의 곱에 비례하는 결합 계수는 GR의 2차 섭동론이 remnant 스핀의 함수로 예측하는 양이라 더 구별력 있는 시험이 된다. 이 계수는 아직 관측 자료에서 직접 측정된 적이 없었다.

Contribution

• GW250114 ringdown에서 $220\times220$ 모드의 진폭 비 $\mathcal R=|\tilde A_{220\times220}|/|\tilde A_{220}|^2$를 처음으로 측정.
• 여러 시작 시각에서 측정해, 스핀 무관 상수 예측(Cheung et al. 2023), NR 기반 선형·비선형 스핀 fit(Cheung et al. 2024; Redondo-Yuste et al. 2024), 2차 섭동론의 스핀 의존 예측(Khera, Ma, Yang 2025)과 비교.
• 고정한 외재 매개변수, 사전분포 폭, $440$ 모드의 존재가 측정에 주는 영향을 점검.

Novelty

이전 GW250114 분석은 2차 모드가 있는지를 Bayes factor로 물었다. 이 논문은 그 진폭이 GR이 예측하는 크기인지를 묻는다. 결합 비는 부모와 자식 모드가 시작 시각 이동에 같은 위상·감쇠 인자를 얻으므로 기준 시각에 무관하다. 그래서 서로 다른 시작 시각의 분석이 같은 양을 측정한다는 점을 측정의 기반으로 삼는다.

Method

• 모형: $220$, $221$, $220\times220$의 세 모드. 2차 모드의 진동수와 감쇠율은 부모의 합이고, 각도 의존성은 $(4,4,0)$ Kerr 스핀 가중 구면파 함수로 근사. 매개변수는 $M_z$, $\chi_f$, 세 모드의 진폭·위상 등 8개이며, 2차 모드 위상은 병합 시각 정의에 민감해 측정 대상에서 뺐다.
• 자료와 추론: 시간 영역 ringdown 분석(bilby + dynesty, live point 1000개). 시작 시각을 피크에서 $t_{M_z}=0.337$ ms 간격으로 바꿔 반복.
• 고정 매개변수: 하늘 위치·편광각은 IMR 최대우도 값으로, 진폭과 강하게 겹치는 광도거리 $D_L=411$ Mpc와 경사각 $\iota=0.76$ rad도 고정.
• 점검: $D_L$과 $\iota$를 IMR $1\sigma$ 범위로 바꾼 재분석, 균일 사전분포 상한을 바꾼 KL 정보량, 모드 구성을 바꾼 GW250114 유사 주입 시험(실제 잡음 10개 실현), $440$을 더한 4모드 재분석.

Results

• 스핀: 피크에 가까운 시작에서는 추론된 $\chi_f$가 IMR 값에서 벗어나고, $3t_{M_z}$ 이후 IMR과 맞는다. 저자들은 이른 구간이 세 모드로 충분히 기술되지 않는다고 해석한다(Fig. 1).
• 결합 비: 사후분포가 넓지만 모든 시작 시각에서 이론 예측을 포함한다. 각 예측의 중앙값이 관측 사후분포에서 차지하는 백분위를 가우스 유의도로 바꾸면 최대 $\approx1.3\sigma$(양측 꼬리 확률 약 $19\%$)다. $(\chi_f,\mathcal R)$ 평면에서는 이론 곡선보다 위쪽으로 퍼지는 경향이 있지만 결정적이지 않다(Fig. 2).
• 정보량: 사전분포 $U(0,5)$에서 이른 시작의 KL 정보량은 약 4 bit이고, $\Delta t=8t_{M_z}$ 근처에서 최소가 된 뒤 일부 회복한다(Fig. 3).
• 외재 매개변수: $\mathcal R$는 $D_L$에 반비례하고 $\iota$가 커질수록 작아진다. 시작 시각에 따른 경향은 유지되지만 전체 크기가 체계적으로 이동한다(Fig. 4).
• $440$ 모드: 주입 시험에서 2차 모드가 없으면 회수된 $\mathcal R$이 0 쪽으로 크게 이동하고, $440$을 넣으면 큰 쪽으로 이동한다(Fig. 5). 실제 자료에 $440$을 더하면 $\Delta t=5t_{M_z}$에서 $\mathcal R=0.207^{+0.086}_{-0.092}\to0.209^{+0.142}_{-0.129}$로 중앙값은 그대로이고 구간만 넓어진다.

Significance

• 비선형 GR 시험의 출발점: 진동수 합 관계보다 구별력 있는 진폭 결합을 관측으로 측정하는 틀을 만들었다. 지금은 예측 넷을 가르지 못하지만, 고SNR 사건이 쌓이면 여러 사건을 결합해 스핀 의존성을 제약할 수 있다.
• 이론과의 연결: 비교 대상인 2차 섭동론 예측(Khera, Ma, Yang)의 저자 중 Khera와 Yang이 3차까지 확장한 2610.20547은, $(2,2)$ 다중극에 cubic 응답이 크게 섞이면 선형 overtone 진폭이 편향된다고 보였다. 이 논문의 모형은 $220$, $221$에 2차 모드 하나만 넣으므로, 그런 고차 응답이 $221$이나 $\mathcal R$ 측정에 어떻게 스며드는지가 다음 질문이 된다.
• 시작 시각 선택의 문제: 스핀이 IMR 값과 맞는 시작 시각과 정보량이 큰 시작 시각이 다르다. 초기 ringdown의 비-QNM 성분을 다룬 QNM 시리즈의 2610.39058과 같은 맥락의 문제다.

Limitations

• 외재 매개변수 고정에 조건부: $D_L$과 $\iota$를 고정했고, 저자들도 측정이 이 값, 특히 경사각에 민감하다고 인정한다. 결과는 "그 값을 가정할 때의 $\mathcal R$"이다. inspiral 정보를 함께 쓰는 반결맞음 분석이 해법으로 제시된다.
• 모형 불일치 가능성: 주입 시험은 $440$ 모드가 $\mathcal R$을 위로 끌어올릴 수 있음을 보여 주며, 관측값이 이론보다 위쪽으로 퍼지는 경향도 이것으로 설명될 수 있다. 4모드 분석은 중앙값을 바꾸지 않았지만 $440$과 $220\times220$을 실제로 분리하지는 못한다. 2차 모드 위상도 측정 대상에서 뺐다.
• 이른 구간의 모형 부족: 스핀 편향이 보이는 피크 근처 시작 시각의 $\mathcal R$도 함께 비교에 쓰였다. 그 구간이 세 모드로 기술되지 않는다면, 그 시작 시각의 결합 비는 다른 성분을 흡수했을 수 있다.
• 유의도 산정 방식: "최대 $1.3\sigma$"는 이론 분포의 중앙값을 관측 사후분포의 백분위로 바꾼 값이다. 이론 예측의 불확실성과 시작 시각 사이의 상관은 이 수치에 들어가지 않는다.
• 단일 사건: 스핀 의존성에 대한 제약은 GW250114의 좁은 스핀 범위($\chi_f\approx0.68$ 근처)에 한정된다.

Takeaway

GW250114에서 처음 측정한 $220\times220$ 2차 모드의 진폭 결합 비는 GR 예측과 $1.3\sigma$ 안에서 맞지만, 예측들을 구별할 정밀도는 없다. 측정은 고정한 거리·경사각과 $440$ 모드 처리에 민감하다.