Swift‑BAT가 밝힌 초거대 질량 X‑ray 이진성 IGR J16493‑4348의 궤도 주기와 전식 현상

초록

Swift‑BAT 54개월 데이터에서 IGR J16493‑4348는 21σ 이상의 검출을 보였으며, 주기 분석을 통해 약 6.78 일의 궤도 주기와 완전 식(전식) 현상을 확인했다. 평균 광도 변동은 약 4배, 단시간(≈300 s) 변동은 최대 20배에 달한다. 케플러 법칙을 이용한 계산으로 반장축은 약 55 R⊙(≈2 R★)이며, 궤도 이심률은 0.15 이하로 추정된다. 이러한 특성은 풍에 의해 물질을 공급받는 풍‑급(supergiant wind‑fed) 시스템임을 시사한다.

상세 분석

본 논문은 Swift‑BAT가 제공하는 15–50 keV 하드 X‑ray 전천구 데이터를 54개월 동안 누적한 뒤, IGR J16493‑4348의 장기 변광 특성을 정밀 분석하였다. 먼저 BAT 이미지에서 20.8σ의 검출 신호를 확인하고, 300 s 시간 해상도로 추출한 광도 시계열을 태양계 질량 중심(SSB) 기준으로 보정하였다. 변동성은 평균 카운트율 1.03 × 10⁻⁴ cnt s⁻¹에 대해 7σ 수준까지 상승했으며, 이는 평균 대비 약 20배의 급격한 플럭스 상승을 의미한다.

주기 탐색은 0.5–100 일 구간을 16위상 구간으로 나누어 folding 기법을 적용했으며, 주기 해상도는 ΔP ≈ P²/(NΔT) 로 정의하였다(ΔT ≈ 1.4 × 10⁸ s). χ² 검정에서 가장 뚜렷한 피크는 P₀ = 6.782 ± 0.005 일이며, χ²≈240이라는 높은 값을 보였다. 피크 주변의 χ² 분포를 가우시안으로 피팅해 오류를 추정하고, 다중 피크(P₁ = P₀/2, P₂, P₃)는 P₀와의 비선형 결합(에클립스와 비트 주파수)에서 유도된 것으로 확인했다.

통계적 유의성을 검증하기 위해 두 가지 방법을 사용하였다. 첫째, χ²–Fχ(P) 형태의 트렌드 함수를 차감한 뒤 z값(=χ²–Fχ) 히스토그램을 구축하고, z>25 구간을 지수함수로 피팅해 z₀ = 196.6 이상의 확률을 4.6 × 10⁻¹⁸로 추정, 이는 8.5σ에 해당한다. 둘째, 원본 시계열의 카운트율을 무작위 재배열한 1000개의 모의 데이터에 대해 동일한 주기 탐색을 수행했으며, 최대 χ²값은 97.9에 불과했다. 따라서 실제 관측값(≈240)의 우연 발생 확률은 <2 × 10⁻⁸, 즉 5.5σ 수준이다.

위상으로 접어본 평균 광도 프로파일은 거의 평탄하지만, 위상 0.319 ± 0.015 근처에서 완전 식을 보이며, 식의 지속시간은 전체 주기의 약 12%에 해당한다. 식 중심 시각은 Tₑcl = 54175.92 ± 0.10 MJD이며, 이는 고궤도 기울기(i≈90°)와 낮은 이심률(e<0.15) 하에서 초거대 별의 반지름(R★≈32 R⊙)에 의해 발생하는 차폐와 일치한다.

Kepler 제3법칙을 적용해 M★≈47 M⊙(B0.5 Ib), MX≈1.4 M⊙(NS) 가정 시 반장축 a≈55 R⊙, 즉 R★의 약 2배가 된다. 이는 전형적인 sgHMXB(초거대 고질량 X‑ray 이진성)의 궤도 규모와 일치한다. 풍‑급 시스템이라는 가정 하에 라그랑주점 L₁이 별의 표면보다 안쪽에 머무르도록 e≤0.15가 제한된다.

Swift‑XRT 0.2–10 keV 관측에서는 평균 카운트율 0.376 cnt s⁻¹, 변동 폭은 약 3배였으며, 5–1000 s 구간에서 추가적인 펄스 신호는 검출되지 않았다. 이는 현재 데이터로는 중성자 별의 스핀 주기를 확인할 수 없음을 의미한다.

결론적으로, BAT 장기 모니터링은 IGR J16493‑4348의 궤도 주기와 전식 현상을 확정했으며, 변동성 특성과 궤도 파라미터는 풍‑급 sgHMXB 모델에 부합한다.