Hydrodynamic Simulations of Tidal Disruption Encores

기술 요약: 조석 파괴 앙코르(Tidal Disruption Encores)의 유체역학 시뮬레이션

문제 정의
핵성단(Nuclear star clusters, NSCs)은 중심의 거대 블랙홀(MBH)을 공전하는 별들과 항성 질량 블랙홀(sBH)의 밀집된 집단을 수용한다. Ryu 등[13]의 이전 분석 연구는 sBH가 별을 조석 파괴할 때(마이크로-TDE), 결합되지 않은 파편의 일부가 중심 MBH에 중력적으로 묶여 있을 수 있다고 제안했다. 이 파편는 나중에 강착되어 "조석 파괴 앙코르(Tidal Disruption Encore, TDEE)"라 불리는 2차 플레어를 생성할 수 있다. Ryu 등은 탄도 궤도를 가정하여 이 사건들을 분석적으로 모델링했으나, 스트림의 자기 교차(self-intersection), 충격파 형성, 원형화(circularization)와 같은 복잡한 유체역학적 상호작용은 탐구되지 않은 채 남아 있었다. TDEE의 형태학, 열역학 및 관측 신호를 이해하는 것은 비정상적인 TDE 유사 플레어를 해석하고 NSC 내의 sBH 집단을 조사하는 데 매우 중요하다.

방법론
저자들은 이동 격자 코드인 AREPO를 사용하여 TDEE에 대한 최초의 유체역학 시뮬레이션을 제시한다. 본 연구는 시스템을 서브-이진계(sBH와 별)가 중심 MBH를 공전하는 구조로 취급하는 "2+1체" 근사를 사용한다.

• 항성 모델: 1 $M_\odot$의 주계열 항성(태양과 유사한 금속성 $Z=0.02$)을 중심 수소 함량이 0.3이 될 때까지 MESA를 사용하여 진화시켰다.
• 시뮬레이션 설정: 별은 비회전 sBH($10 M_\odot$)에 의해 파괴되며, 동시에 sBH-별 이진계는 MBH를 공전한다. 시뮬레이션은 MBH 질량($10^3 M_\odot$인 직접 앙코르; $10^6 M_\odot$인 링 앙코르)과 이진계-MBH 분리 거리($R_{bin}$)를 변화시킨다. MBH에 대한 이진계 근점(periapsis)의 방향은 각도 $\theta$로 매개변수화된다.
• 해상도 및 물리: 시뮬레이션은 $N = 5 \times 10^5$개의 셀을 사용한다. 상태 방정식에는 국소 열역학적 평형 하의 복사압이 포함된다. 명시적인 복사 전달은 포함되지 않지만, 광자 확산 시간($t_{diff}$)과 광학적 깊이를 고려하여 가상의 원통형 격자를 통해 복사 플럭스를 적분함으로써 볼로미터 휘도를 사후 처리(post-processing)를 통해 추정한다.
• 분석: 저자들은 서로 다른 TDEE 결과들을 구분하기 위해 파편의 형태학, 점성 시간($t_{visc}$), 그리고 광도 곡선을 분석한다.

주요 기여 및 결과
시뮬레이션 결과, TDEE의 결과는 균일하지 않으며 파괴 기하학, MBH 질량 및 분리 거리에 매우 민감하게 의존함을 보여준다. 저자들은 두 가지 뚜렷한 형태학적 클래스를 식별하였다:

1. 직접 앙코르 (Direct Encores):

   • 조건: MBH 질량이 낮거나($\sim 10^3 M_\odot$, 중간 질량 블랙홀), 분리 거리가 작을 때($R_{bin} \sim 10^{-5}$ pc), 또는 MBH에 대한 파편의 에너지가 음수이지만 0에 가까울 때 발생한다.
   • 형태학: 파편 스트림은 MBH를 향해 매우 편심된 탄도 궤적을 따른다. 이들은 자유 낙하 시간($t_{ff}$) 동안 자기 교차 및 충격 가열을 겪는다.
   • 관측 특성: 이 사건들은 $t \sim t_{ff}$(약 며칠에서 몇 주)에서 정점에 도달하는 2차 플레어를 생성한다. 휘도는 $10^{40}$에서 $10^{41}$ erg/s 범위이다. 광도 곡선은 충격으로 인한 급격한 상승을 보인 후 감소한다. $M_{MBH} \sim 10^3 M_\odot$인 경우, 강착률은 초-에딩턴(hyper-Eddington) 수준에 달할 수 있으며, 이는 유출물(outflow)이나 제트를 구동할 수 있다.

2. 링 앙코르 (Ring Encores):

   • 조건: MBH가 더 무겁고($\sim 10^6 M_\odot$), 분리 거리가 더 클 때($R_{bin} \sim 10^{-4}$ pc) 발생한다. 이 경우 마이크로-TDE로부터의 에너지 편차는 sBH의 원래 경로에 비해 파편의 궤도를 크게 변화시키기에 불충분하다.
   • 형태학: 파편은 직접 낙하하는 대신 약 20번의 궤도 주기 동안 토러스(torus) 또는 링 형태로 원형화된다. 링은 약 한 번의 궤도 주기($P$)의 지연 후에 형성된다.
     기초적인 마이크로-TDE 이후, 링 형성에 해당하는 휘도의 고원(plateau) 현상이 나타난다. 초기 휘도는 점성 시간($t_{visc}$)보다 훨씬 짧은 광자 확산 시간($t_{diff}$) 동안 복사되는 열에너지에 의해 구동된다. 링이 MBH로 점성 강착됨에 따라 수년에서 수십 년 후에 2차 재밝기(rebrightening)가 발생할 수 있다. 정점 휘도는 $\sim 10^{42}$ erg/s에 달한다.

의의 및 주장
본 논문은 이러한 유체역학 시뮬레이션이 기존의 분석 모델에 비해 TDEE 광도 곡선의 예측력을 크게 향상시킨다고 주장한다.

• 형태학적 구분: 본 연구는 TDEE가 단일한 현상이 아니라 서로 다른 시간 척도(자유 낙하 vs 궤도 주기)와 휘도 진화를 가진 뚜렷한 "직접" 및 "링" 징후를 나타낸다는 것을 입증하였다.
• NSC 역학 조사: 저자들은 TDEE가 NSC 내의 sBH 집단을 조사하는 프로브 역할을 한다고 상정한다. 일차 마이크로-TDE와 앙코르 사이의 시간 지연은 MBH 질량과 sBH-MBH 분리 거리에 대한 정보를 인코딩한다.
• 이상 현상 설명: TDEE가 예측하는 다양한 감쇠 기울기와 지연된 재밝기는 표준 $t^{-5/3}$ 거듭제곱 법칙에서 벗어나는 관측된 광학/UV TDE 후보들을 설명할 수 있다.
• IMBH 탐지: 특히 중간 질량 블랙홀($10^3 M_\odot$)을 포함하는 직접 앙코르는, 다른 방법으로는 제약하기 어려운 NSC 내의 IMBH를 탐지하기 위한 독립적인 진단 도구를 제공한다.
• 관측 후보: 저자들은 $\sim 10^6 M_\odot$ MBH를 포함하는 TDEE의 예측된 특성과 일치하는, 주요 정점 이전에 광학/UV 범프를 보였던 사건 AT 2023lli를 제시한다.

결론적으로, TDEE는 새로운 종류의 핵 트랜지언트(nuclear transient)를 나타낸다. 현재의 시뮬레이션은 뉴턴 역학 및 특정 항성 질량으로 제한되어 있으나, 다파장 트랜지언트를 해석하고 Rubin, ULTRASAT, 그리고 X선 미션(eROSITA, Athena)과 같은 미래의 관측 시설을 안내하는 강력한 프레임워크를 제공한다.