Pulse profile modelling of the accretion-powered millisecond pulsar SAX J1808.4-3658 using NICER data from its 2019 and 2022 outbursts

기술 요약: NICER 데이터를 이용한 SAX J1808.4−3658의 펄스 프로파일 모델링

문제 정의
펄스 프로파일 모델링(PPM)은 회전 동력 밀리초 펄서(RMP)에 성공적으로 적용되어 중성자별(NS)의 질량과 반경을 제약함으로써 고밀도 물질의 상태 방정식(EoS)을 조사하는 데 사용되는 상대론적 광선 추적 기법이다. 그러나 강착 동력 밀리초 펄서(AMP)에 PPM을 적용하는 것은 뚜렷한 과제를 안겨준다. RMP와 달리, AMP는 강착 원반과 강착 흐름으로부터 발생하는 상당한 비펄스 X선 방출을 보이며, 이는 배경 신호로부터 펄스 형태의 핫스팟 신호를 분리하는 것을 복잡하게 만든다. Salmi 등(2018)이 RXTE 데이터를 사용하여 수행한 이전의 AMP 적용 사례들은 컴프턴 산란(Comptonization)에 대한 경험적 모델과 아카이브 데이터에 의존했다. 본 연구는 새로운 NICER(Neutron star Interior Composition Explorer) 데이터, 구체적으로는 2019년과 2022년 폭발기 동안의 전형적인 AMP인 SAX J1808.4−3658(J1808)을 대상으로 물리적으로 동기화된 현대적 PPM 기법을 적용할 필요성을 다룬다. 주요 과제는 매개변수 추론을 편향시키는 퇴화(degeneracy)를 피하기 위해 강착 원반의 기여도와 핫스팟의 복잡한 기하학적 구조를 정확하게 모델링하는 것이다.

방법론
저자들은 J1808의 2019년 및 2022년 NICER 폭발 데이터의 펄스 프로파일을 피팅하기 위해 X-PSI 소프트웨어 패키지를 사용하여 베이지안 매개변수 추론을 채택하였다. 분석에는 편평한 슈바르츠칠트 근사(oblate Schwarzschild approximation), 중력 적색편이, 빛의 굴절 및 도플러 부스팅을 포함하는 상대론적 광선 추적 프레임워크가 활용되었다.

주요 방법론적 구성 요소는 다음과 같다:

• 핫스팟 모델: 단일 원형 핫스팟 모델(단일 온도, ST)과 두 개의 핫스팟 모델(ST-U, 매개변수 미공유)을 비교하였다. 핫스팟 대기는 뜨거운 전자의 슬래브 기하학 내에서의 컴프턴 산란을 가정하는 Bobrikova 등(2023)의 컴프턴 슬래브 대기를 사용하여 모델링되었다.
• 배경/원반 모델링: 비펄스 방출에 대해 세 가지 접근 방식을 테스트하였다:
  1. 원반 모델: 다색 원반 블랙바디(diskbb) 모델.
  2. 원반 + 라인 모델: 1 keV 부근에서 관찰되는 확장된 반사 특징을 설명하기 위해 diskbb 모델에 가우시안 라인을 추가한 모델.
  3. 한계적 배경(Marginalised Background): 각 에너지 빈(bin)에 대해 독립적으로 카운트를 한계화(marginalise)하여 배경 성분을 유연하게 처리하는 방식이며, 기준 플럭스와 곱셈 지지 계수($s$)에 의해 제약된다.
• 사전 분포(Priors): 질량($M$), 적도 반경($R_{eq}$), 거리, 경사각 및 핫스팟 기하학에 대한 사전 분포를 설정하였다. 넓고 무정보적인 사전 분포 외에도, 계산 효율성과 매개변수 일관성을 탐구하기 위해 Rutherford 등(2024)의 고밀도 물질 EoS 제약을 기반으로 한 EoS 정보 사전 분포를 테스트하였다.
• 데이터 선택: 충분한 광자 수를 확보하기 위해 폭발의 정점(2019년 8월 및 2022년 8월) 데이터를 선택하였으며, 업데이트된 타이밍 솔루션을 사용하여 펄스 프로파일을 폴딩(folding)하였다.

주요 결과

1. 단일 핫스팟 모델의 불충분성: 단일 핫스팟 모델(ST)은 모든 배경 처리 방식에서 데이터를 적절히 피팅하는 데 실패하였다. 특히 0.8–1.1 keV 부근의 스펙트럼 범프와 위상 의존적 잔차가 남았는데, 이는 단일 원형 핫스팟이 관측된 펄스 프로파일을 설명할 수 없음을 나타낸다.
2. 원반 배경을 포함한 두 핫스팟 모델: 두 번째 핫스팟(ST-U)을 도입하면 위상 잔차는 개선되지만, 물리적 원반 모델(diskbb 또는 diskbb + Gaussian)을 사용할 때 스펙트럼 잔차를 해결하는 데는 실패한다. 이 모델은 여전히 특정 밴드에서 카운트를 과대평가하고 1 keV 근처의 좁은 스펙트럼 특징을 과소평가한다. 추론된 강착율 또한 독립적인 추정치보다 현저히 낮았다.
3. 한계적 배경의 성공: 가장 좋은 통계적 피팅은 두 핫스팟 모델(ST-U)과 한계적 배경 접근법을 결 함께 사용했을 때 달성되었다. 이 구성은 적절한 피팅을 보여주는 $\chi^2$/DOF 값(2019년과 2022년에 대해 각각 0.17 및 0.30의 p-value)을 산출하였다.
4. 매개변수 퇴화 및 견고성: 개선된 피팅에도 불구하고, 추론된 매개변수(질량, 반경, 경사각)는 비펄스 배경과 핫스팟의 비펄스 성분 사이의 강한 퇴화로 인해 견고하지 않은 것으로 나타났다.
   • 한계적 배경의 지지 계수($s$)를 100에서 1000으로 증가시켰을 때, 추론된 질량은 약 $2.0 M_\odot$에서 $1.8 M_\odot$으로 감소했고, 반경은 약 8.7 km에서 8.1 km로 이동했다.
   • 지지 경계의 선택이 사후 분포를 크게 변화시켰으며, 이는 배경에 대한 물리적 상한선 없이 매개변수가 유일하게 제약될 수 없음을 시사한다.
5. EoS 정보 사전 분포: EoS 정보 사전 분포를 적용하면 계산 비용이 약 2.5배 감소하고 추론된 반경이 사전 분포와 일치하는 더 큰 값($\sim 12-13$ km)으로 이동하지만, 베이지안 증거를 크게 개선하거나 퇴화 문제를 해결하지는 못했다.

의의 및 주장
저자들은 본 연구를 고품질 NICER 데이터에 대한 PPM 기법의 적용을 입증하며, AMP의 PPM을 위한 중요한 다음 단계로 자리매김한다. 논문의 주장은 다음과 같다:

• 모델 복잡성의 필요성: J1808에 단일 핫스팟 모델은 불충분하며, 펄스 프로파일을 설명하기 위해서는 두 개의 핫스팟 기하학이 필요하다.
• 원반 모델링의 결함: 현재의 단순한 강착 원반 모델(blackbody + Gaussian line)은 AMP의 스펙트럼 특징, 특히 1 keV 주변의 확장된 반사선을 포착하기에 부적절하다.
• 유연한 배경 모델의 한계: 한계적 배경은 우수한 통계적 피팅을 제공하지만, 비펄스 복사의 기원(원반 vs 핫스팟)을 분리할 수 없는 퇴화를 유발한다. 결과적으로, 배경에 대한 더 엄격한 물리적 제약 없이는 추론된 질량과 반경이 견고하지 않다.
• 미래 요구 사항: 본 논문은 AMP에 대한 견고한 매개변수 제약을 위해서는 더 정확한 강착 원반 및 핫스팟 모델링이 필요하다고 결론짓는다. 특히, 고에너지 데이터와 편광 데이터(예: IXIXE)의 포함이 경사각, 핫스팟 위도, 크기에 대한 보완적인 제약을 제공하여 현재의 퇴화를 해소할 것으로 기대된다.

저자들은 보고된 매개변수가 최종적인 것이 아니며, 강착 원반 및 핫스팟 기하학의 물리적 모델링 개선에 따라 재조정될 수 있음을 명시적으로 밝히고 있다.