The Very Late Time Afterglow of GW170817 Favors a Wobbling Jet

기술 요약: GW170817 의 매우 늦은 시점 잔광은 요동하는 제트를 지지함

문제 제기
이중 중성자별 (BNS) 병합 사건 GW170817 은 중력파, 감마선 폭발 (GRB 170817A), 그리고 킬로노바가 동시에 관측된 유일한 사건으로 남아 있습니다. 정렬된 구조화된 제트 모델은 잔광의 상승 단계와 겉보기 초광속 운동을 성공적으로 설명해 왔으나, 매우 늦은 시점의 잔광 (약 10 년간 관측됨) 에서 결정적인 불일치에 직면해 있습니다. 표준 모델은 제트 붕괴 후 정렬된 제트가 측면으로 퍼지면서 $F_\nu \propto t^{-p}$ (여기서 $p > 2$는 전자 파워 지수임) 의 시간적 감쇠를 예측합니다. 그러나 GW170817 의 관측 결과는 훨씬 더 완만한 감쇠율 ($\alpha \lesssim 1.8$) 을 보여줍니다. 관측되지 않은 추가 구성 요소 (킬로노바 잔광이나 리프레시드 충격파 등) 나 극단적인 물리적 가정을 도입하지 않고 이 불일치를 해결하는 것은 여전히 과제입니다.

방법론
콜랩사 (collapsars) 와 중성자별 - 블랙홀 병합에 대한 일반상대론적 자기유체역학 (GRMHD) 시뮬레이션은 낙하하는 가스의 상호작용이나 비대칭적인 각운동량 방출로 인해 강착 원반이 기울어질 수 있음을 보여주는데, 이에 영감을 받아 저자들은 GW170817 의 제트가 '요동 (wobbling)'하고 있다고 제안합니다.

• 이론적 모델: 저자들은 요동하는 제트를 하늘 위에서 연속적인 띠를 끌고 가는 고리 모양의 구조로 근사화합니다. 그리고 초기 감속, 가시 호 (arc) 길이가 늘어나는 중간 단계, 그리고 전체 고리가 넓은 제트로 합쳐지는 최종 단계 등 세 단계를 고려하여 잔광의 시간적 진화를 해석적으로 유도합니다.
• 해석적 유도: 중간 단계에서 가시 방출 영역은 한 차원 (호 길이) 에서 증가하는 반면, 빔 콘은 두 차원에서 확장됩니다. 이로 인해 고유한 감쇠 지수 $F_\nu \propto t^{-3(2p-1)/7}$이 도출됩니다. 전형적인 $p$ 값 (2 에서 2.5) 에 대해 이는 1.3 에서 1.7 사이의 감쇠 지수를 산출하여 관측된 완만한 기울기와 일치합니다.
• 수치 시뮬레이션 및 베이지안 분석: 저자들은 반해석적 모델보다 고리 구조 내의 복잡한 퍼짐을 더 정확하게 반영하는 유체역학 코드 jetsimpy 를 사용하여 충격파 역학을 모델링합니다. 그들은 전파 및 X 선 광곡선, 겉보기 초광속 운동 측정치, 그리고 경사와 광도 거리에 대한 중력파 사후 표본을 포함한 전체 다중신호 데이터셋에 대해 베이지안 분석을 수행합니다. 그리고 표준 파워-법칙 구조화된 제트와 제안된 고리 모양 (요동) 제트 두 모델을 비교합니다.

주요 기여 및 결과

• 모델 검증: 고리 모양 제트 모델은 겉보기 초광속 운동 제약 조건 (정렬된 '날개' 구조 필요) 과 일관성을 유지하면서 완만한 늦은 시점 감쇠 ($\alpha \approx 1.3-1.7$) 를 성공적으로 재현합니다.
• 통계적 유의성: 베이지안 모델 비교는 $\log_{10} \text{BF} = 6.16$의 로그 베이지스 인자를 산출하여 $4.8\sigma$의 유의수준에 해당합니다. 이는 표준 파워-법칙 구조화된 제트보다 고리 모양 제트 모델을 강력하게 지지합니다.
• 모수 추정: 분석은 $\theta_{wo} \approx 27^\circ$의 요동 각도를 추정합니다. 고리 모양 제트에 대한 최적 적합 모수는 반개방각 $\theta_c \approx 0.03$ rad 와 관측자 각도 $\theta_{obs} \approx 0.73$ rad 를 포함합니다.
• 물리적 함의: 추정된 큰 요동 각도는 강착 원반의 상당한 기울기를 시사하며, 이는 BNS 병합 중 원반 - 낙하 가스 상호작용이나 비대칭적인 각운동량 방출에서 비롯될 수 있습니다.

의의 및 주장
본 논문은 요동하는 제트 모델이 추가 구성 요소나 비표준 물리 과정을 도입하지 않고 GW170817 다중신호 데이터에 대한 자기일관적인 설명을 제공한다고 주장하며, 구조화된 제트 모델에 대한 최소한의 확장으로서 오컴의 면도날 원칙을 준수합니다.

• "누락된 제트 붕괴"의 해결: 저자들은 이 모델이 감쇠 기울기가 1.4–1.8 인 다른 감마선 폭발들에서 나타나는 '누락된 제트 붕괴' 문제를 해결할 수 있음을 시사하며, 요동하는 제트가 흔할 수 있음을 내포합니다.
• 우주론적 영향: 이 모델은 정렬된 제트 모델에 비해 더 큰 관측 각도를 산출하여 광도 거리 제약 ($d_L \approx 36.4$ Mpc) 을 줄입니다. 이는 $H_0 \approx 80.8$ km/s/Mpc 의 허블 상수 추정을 초래하여 Ia 형 초신성 추론을 약간 지지하지만, 저자들은 요동 축이 궤도 축과 정렬되지 않을 경우 이 값이 체계적 오차를 겪을 수 있음을 지적합니다.
• 병합율 추론: 고리 모양 구조는 정렬된 제트 ($\pi\theta_c^2$) 에 비해 더 큰 입체각 ($4\pi\theta_c\theta_{wo}$) 을 의미합니다. 저자들은 이것이 중력파에서 추론된 BNS 병합율과 짧은 감마선 폭발에서 추론된 병합율 사이의 불일치를 해결할 수 있으며, 후자가 과대평가되었을 수 있음을 시사한다고 주장합니다.
• 한계: 저자들은 고리 모양 제트가 '토이 모델 (toy model)' 근사임을 명시적으로 밝힙니다. 그들은 시뮬레이션이 밝혀낸 실제 유출 구조가 매우 비대칭적이고 복잡함을 인정하며, 이러한 신호를 더 깊이 연구하기 위해 효율적인 3 차원 잔광 모델링 도구의 개발을 요구합니다.