HLX‑1 변광의 비밀: 디스크 불안정성보다 별의 궤도 질량이송이 주인공

초록

ESO 243‑49 은하 외곽에 위치한 초고광도 ULX HLX‑1은 수백 일 주기의 강렬한 X‑ray 변광을 보인다. 저자들은 디스크 열‑점성 불안정성(DIM) 모델을 적용했지만, 필요한 디스크 크기와 점성 시간척도가 관측된 10 일‑수주 스케일과 맞지 않음을 확인한다. 대신, 고궤도 이심률을 가진 별이 근일점에서 질량을 떼어내는 ‘조석 질량이송’이 변광을 일으킬 가능성을 제시하고, 이러한 시스템이 형성될 확률을 별밀도와 동역학적 충돌을 고려해 1 %–10 % 수준으로 추정한다.

상세 분석

본 논문은 HLX‑1의 장기 변광 메커니즘을 두 가지 물리적 모델로 검증한다. 첫 번째는 전통적인 디스크 열‑점성 불안정성(DIM)이다. DIM은 디스크가 충분히 넓어 온도가 수소 재결합 온도(≈ 10⁴ K) 이하가 되는 영역을 포함해야 하며, 그 경우 임계 질량이송률 (\dot M_{\rm crit}) ≈ 3 × 10⁻⁶ M⊙ yr⁻¹(α ≈ 0.1, R ≈ 10¹³ cm) 를 초과하면 급격한 FRED 형태의 폭발이 일어난다. HLX‑1의 관측된 피크 광도(≈ 10⁴² erg s⁻¹)와 변광 주기(≈ 380 d)로부터 추정되는 디스크 반지름은 최소 5 × 10¹³ cm이며, 이는 DIM이 요구하는 최소 반지름보다 크다. 그러나 DIM이 예측하는 점성 시간척도 (t_{\rm vis}\sim) 수년에서 수십년 수준이며, 이는 실제 관측된 급격한 상승(≈ 1 주)과 90 일 내의 감쇠와 크게 불일치한다. 또한, 디스크가 방사선 압력에 의해 지배되는 경우( (L/L_{\rm Edd}>0.01) ) 내부 영역에서만 열‑점성 불안정성이 발생하고, 재발 주기가 수백 년에 달한다는 기존 시뮬레이션 결과와도 모순된다. 따라서 저자들은 DIM이 HLX‑1의 변광을 설명하기에 물리적으로 부적절하다고 결론짓는다.

두 번째 모델은 고이심률 궤도를 가진 별이 근일점에서 조석적으로 물질을 잃어버리는 ‘조석 질량이송’이다. 질량비 (q=M_{\bullet}/M_{\star}\sim10^{-4})–10⁻³인 경우 로체 한계 반지름 (R_L\approx a(q/3)^{1/3}) 로 표현되며, 별이 적색거성 혹은 초거성 단계에 이르면 평균 밀도 (\bar\rho\approx10^{-4}) g cm⁻³에 해당한다. 이 경우 궤도 주기 380 d는 반지름 (a\approx2.3\times10^{13}) cm에 해당하고, 근일점에서 별의 외피가 로체 한계보다 크게 팽창하면 급격한 질량이송이 발생한다. 전이된 물질은 즉시 블랙홀 주변에 원반을 형성하고, 짧은 자유낙하 시간(≈ 수일) 내에 고밀도 흐름을 만들어 관측된 급격한 상승을 재현한다. 저자들은 이 과정이 디스크의 열‑점성 불안정성보다 훨씬 빠른 시간척도로 변광을 일으킬 수 있음을 강조한다.

시스템 형성 확률을 추정하기 위해, 저자들은 젊은 고밀도 별단지에서 런어웨이 별 합병으로 10⁴ M⊙ 이하의 중간질량 블랙홀(IMBH)이 형성되는 과정을 모델링했다. 이후 별단지 내에서 1 M⊙ 수준의 주계열 혹은 AGB 별이 고이심률 궤도로 끌려가 로체 한계에서 질량을 잃는 경우는 전체 IMBH 중 1 %–10 % 정도에 불과하다고 계산한다. 충돌 파괴나 2체 상호작용에 의한 궤도 재배열을 고려하면 실제 관측 가능한 HLX‑1과 같은 시스템은 매우 드물지만, 존재 자체는 통계적으로 배제되지 않는다.

마지막으로, 저자들은 궤도 이심률이 거의 포물선에 가까운 경우( (e\rightarrow1) ) 질량이송이 순간적으로 급증하고, 디스크가 즉시 고온 상태로 전이되면서 관측된 FRED 형태와 짧은 재발 주기를 완벽히 재현할 수 있다고 제안한다. 다만, 이러한 극단적인 이심률 궤도는 형성 확률을 더욱 낮추므로, 추가적인 중력 포텐셜 변동에 의한 열원(예: 강제 진동에 의한 마찰 가열) 등을 고려해야 할 필요가 있다.