Revisiting GW170817 at milliarcsecond scale: high-precision constraints on jet geometry and $H_0$

기술적 요약: GW170817 의 밀리초각 규모 재검토

문제 제기
이중 중성자별 병합 GW170817 의 탐지는 허블 상수 ($H_0$) 에 대한 '표준 사이렌' 측정의 첫 기회를 제공했습니다. 그러나 이 측정의 정밀도는 기본적으로 중성자별 쌍성계의 광도 거리 ($D_L$) 와 경사각 ($\iota$) 간의 강한 공변성에 의해 근본적으로 제한됩니다. 중력파 (GW) 데이터만으로는 GW170817 과 같은 정면 (face-on) 사건의 경사각을 정밀하게 제약할 수 없지만, 병합으로 방출된 상대론적 제트는 이 기하학적 구조에 대한 전자기적 단서를 제공합니다. 잔광 광곡선을 사용하여 관측각 ($\theta_v$) 을 제약하려는 이전 시도들은 관측각과 제트 개구각 ($\theta_{cp}$) 간의 기하학적 퇴화성으로 인해 제한되었습니다. 반면, 제트의 초광속 고유 운동에 대한 초장기선 간섭계 (VLBI) 관측은 직교하는 제약을 제공하지만, 역사적으로 이산적이며 반분석적인 제트 모델 피팅이나 매개변수 공간의 제한된 부분집합에 의존해 왔으며, 종종 완전한 공간 정보 (예: 소스 크기) 나 체계적 불확실성을 엄밀하게 고려하지 않았습니다.

방법론
본 연구는 GW170817 의 관측 기하학 및 이를 통해 얻어지는 $H_0$에 대한 견고한 제약을 도출하도록 설계된 포괄적인 베이지안 가시성면 (visibility-plane) 모델 피팅 프레임워크를 제시합니다. 방법론은 다음 구성 요소들을 통합합니다:

1. 데이터 통합: 분석은 사용 가능한 네 가지 밀리초각 규모 측량 데이터 세트를 모두 포함합니다:

   • Mooley 등 (2018) 의 병합 후 75 일 및 230 일 두 시점의 고감도 어레이 (HSA) VLBI 데이터.
   • Ghirlanda 등 (2019) 의 병합 후 207 일 시점의 글로벌-VLBI 데이터 (저자들은 HSA 보정과의 일관성을 보장하기 위해 이를 독립적으로 재감산함).
   • 병합의 기준점 (Day 0) 으로 작용하는 병합 후 8 일의 킬로노바 위치에 대한 허블 우주 망원경 (HST) 측정치.

2. 유체역학 모델링: 이산적이고 사전 계산된 모델에 의존하는 대신, 저자들은 2 차원 상대론적 시뮬레이션에 기반한 새로운 유체역학적 제트 잔광 모델의 연속체에 적용되는 스케일링 형식을 사용합니다 (Govreen-Segal & Nakar 2023). 이러한 모델들은 네 가지 제트 에너지 분포 (상자형 및 지수 $b=3, 4, 6$인 멱함수) 를 다룹니다. 스케일링은 모델이 특정 매개변수 비율 하에서 불변성을 갖는 특성을 활용하여 임의의 관측각 ($\theta_v$), 개구각 ($\theta_{cp}$), 그리고 관측자 시간에 맞게 모델을 조정할 수 있게 합니다.

3. 베이지안 추론: 프레임워크는 GPU 가속 소프트웨어 GALARIO와 중첩 샘플러 Dynesty( Bilby 를 통해) 를 사용하여 관측된 전파 간섭계 가시성 (visibilities) 에 직접 작동합니다. 이 접근법은 이미지면 디컨볼루션과 관련된 정보 손실을 피합니다. 모델은 각 시점에 대한 측량 오프셋 ($\delta_{ra}, \delta_{dec}$) 및 플럭스 스케일 불확실성 ($\delta_F$) 을 포함한 체계적 불확실성을 엄밀하게 처리하기 위해 방해 변수 (nuisance parameters) 를 포함합니다.

4. 우주론적 제약: $H_0$를 측정하기 위해 프레임워크는 GW 가능도 (경사각과 거리에 대한 커널 밀도 추정기를 통해) 및 고유 속도 보정을 통합합니다. 저자들은 28 가지 은하군 목록과 고유 속도 재구성 지도의 조합을 고려하며, 이러한 체계적 불확실성에 대해 주변화하기 위해 베이지안 모델 평균화 (BMA) 를 적용합니다.

주요 기여

• 연속 매개변수 공간: 이 연구는 이산 모델을 피팅하는 것을 넘어 제트 기하학 ($\theta_v$ 및 $\theta_{cp}$) 의 연속적인 범위를 탐색함으로써 특정 제트 구조를 선택함으로써 발생하는 편향을 피합니다.
• 가시성면 피팅: 이미지면 중심점에 직접 피팅하는 대신 가시성에 직접 피팅함으로써 소스 크기 제약 등을 포함한 데이터의 완전한 공간 정보를 활용하여 퇴화성을 깨는 데 결정적인 역할을 합니다.
• 체계적 불확실성 처리: 측량 및 플럭스 보정 불확실성을 위한 방해 변수의 포함과 고유 속도 보정을 위한 BMA 접근법은 고정된 보정을 가정했던 이전 연구들보다 더 견고한 오차 예산을 제공합니다.
• 독립적인 글로벌-VLBI 재감산: 저자들은 HSA 데이터와의 프레임 일관성을 보장하기 위해 글로벌-VLBI 데이터를 독립적으로 재감산하여 이전 감산 결과에 대한 측량 오프셋을 식별하고 보정했습니다.

결과

• 관측각 제약:
  • 표면 밝기 요동 측정을 기반으로 한 고정된 광도 거리 $D_L = 40.7$ Mpc 를 가정할 때, 광곡선 기반 제트 기하학 사전 정보를 통합하면 관측각은 $\theta_v = 18.3^\circ - 20.3^\circ$(68% 신뢰 구간) 로 도출됩니다.
  • $D_L$과 고유 속도 보정에 대해 주변화할 때, 선호되는 관측각은 $\theta_v = 16.8^\circ - 19.2^\circ$로 약간 이동합니다.
• 광도 거리: 거리에 대한 직접 피팅은 $D_L = 44.0 \pm 1.6$ Mpc 를 산출합니다. 이는 표면 밝기 요동 방법에서 유도된 $40.7 \pm 2.4$ Mpc 값보다 약 $1\sigma$ 크게, 관측각 제약이 가정된 거리에 얼마나 민감한지를 보여줍니다.
• 허블 상수 ($H_0$): 28 가지 고유 속도 보정 시나리오에 대한 베이지안 모델 평균화 측정은 $H_0 = 65.5 \pm 4.4$ km s$^{-1}$ Mpc$^{-1}$입니다.
  • 이 값은 초기 우주 플랑크 (Planck) 협력의 값 ($67.4 \pm 0.5$ km s$^{-1}$ Mpc$^{-1}$) 과 $0.5\sigma$ 이내입니다.
  • 후기 우주의 SH0ES 측정치 ($73.04 \pm 1.04$ km s$^{-1}$ Mpc$^{-1}$) 보다 $1.7\sigma$ 낮습니다.

의의 및 주장
본 논문은 이 독립적인 후기 우주 $H_0$ 측정이 SH0ES 결과보다 초기 우주 플랑크 결과와 더 일관된다고 주장합니다. 저자들은 이를 있는 그대로 받아들일 경우, 허블 긴장 (Hubble tension) 을 해결하기 위한 비표준 우주론 모델의 필요성에 도전할 수 있음을 시사합니다. 그러나 저자들은 이 주장에 대해 겸손하게 접근하며 다음과 같은 몇 가지 주의점을 명시합니다:

• 모델 오지정: 잔광 광곡선의 후기 시간 진화 ($\gtrsim 350$일) 는 표준 싱크로트론 모델로 잘 피팅되지 않으며, 과잉 방출의 원인은 아직 알려지지 않았습니다. 저자들은 자신의 분석 (초기 시점에 초점) 이 모델의 부정확성보다는 기하학에 의해 지배된다고 주장하지만, 모델링되지 않은 물리가 결과를 편향시킬 수 있음을 인정합니다.
• 제트 정렬: 분석은 제트가 궤도 각운동량 축을 따라 방출된다고 가정합니다. 만약 상당한 정렬 오차가 존재한다면, $H_0$ 추정에 비선형적으로 편향을 일으킬 수 있습니다.
• 고유 속도 우세: 저자들은 관측각 제약 자체보다는 고유 속도 보정을 $H_0$ 측정의 주요 불확실성 원인 (약 5% 의 분수 불확실성 기여) 으로 식별합니다.

본 연구는 제트 잔광에 대한 고정밀 VLBI 관측이 당분간 경사각을 제약하는 데 있어 독보적으로 강력한 도구임을 결론지으며, 허블 긴장에 대한 결정적인 해결책을 달성하려면 더 많은 밝은 사이렌 표본과 제트 물리 및 고유 속도에 대한 개선된 모델링이 필요하다고 강조합니다.