Delayed radio emission in tidal disruption events from collisions of outflows driven by disk instabilities

기술 요약: 디스크 불안정성에 의해 구동되는 유출물 간의 충돌로부터 발생하는 지연된 라디오 방출

문제 제기
조석 파괴 사건(Tidal Disruption Events, TDEs)에서 광학 피크 이후 수백에서 수천 일 후에 나타나는 지연된 라디오 방출이 점점 더 많이 관측되고 있다. 즉각적인 라디오 방출은 일반적으로 파괴 직후 비상대론적 유출물과 주변 핵 주위 매질(CNM) 사이의 상호작용에 기인하지만, 이 지연된 약상대론적 유출물을 구동하는 메커니즘은 여전히 불분명하다. 기존의 가설에는 오프-액시스(off-axis) 제트, 스트림-스트림 충돌, 또는 불균일한 CNM 구름과의 상호작용 등이 포함된다. 그러나 약 100~1000일 후 파괴로부터 방출되는 유출물, 특히 다중 플레어나 뚜렷한 후기 피크를 보이는 유출물의 구체적인 기원은 추가적인 이론적 조사가 필요하다. 본 연구는 이러한 유출물이 TDE에 의해 형성된 강착 원반의 불안정성이나 상태 전이로부터 발생하는 질량 방출("플레어")과 그 유출물들이 서로 충돌하거나 주변 매질과 충돌함으로써 발생할 수 있다고 가정한다.

방법론
저자들은 질량 유출물 간의 충돌에 따른 수력학적 진화와 그에 따른 싱크로트론 라디오 방출을 모델링한다. 연구에서는 두 가지 주요 충돌 시나리오를 고려한다:

1. 플레어+플레어(Flare+Flare): 강착 원반에서 시간 간격 $\Delta t$를 두고 연속적으로 방출된 두 개의 플레어 간의 충돌.
2. 플레어+CNM(Flare+CNM): 단일 플레어와 주변 핵 주위 매질(CNM) 사이의 충돌.

유출물의 특성은 고강착 상태의 디스크에서 열적 불안정성이 질량 방출을 유도한다는 A. L. Piro & B. Mockler (2025)의 디스크 불안정성 시나리오를 바탕으로 한다. 플레어는 밀도 프로파일 $\rho_{fl} \propto r^{-2s-3}$ (여기서 $s=0.5$)와 속도 범위 $v_{min} < v < v_{max}$ (구체적으로 $0.04c$에서 $0.4c$)로 모델링된다. 저자들은 플레어 질량($M_{fl}$)을 $0.01$및$0.1 M_\odot$로 탐색한다.

CNM의 경우, 두 가지 밀도 프로파일이 테스트된다: 상수 밀도 ($\alpha=0$)와 거듭제곱 법칙 감소 ($\rho \propto r^{-2.5}$, $\alpha=2.5$)로, 이는 관측적으로 추론된 기울기들을 포괄한다. 충돌에 의해 형성된 충격파의 수력학적 진화는 얇은 쉘(thin shell)에 대한 질량 및 운동량 보존 방정식을 사용하여 계산된다. 결과적인 라디오 방출은 자유-자유 흡수(FFA)와 싱크로트론 자기 흡수(SSA)를 고려한 싱크로트론 복사 이론을 사용하여 계산된다. 모델은 전력 법칙 전자 에너지 분포(지수 $p=3$)와 등분배 파라미터 $\epsilon_B = \epsilon_e = 0.1$을 가정한다.

주요 결과
시뮬레이션 결과, 충격 충돌은 $\sim 10^{17}$ cm 반경에서 충격 속도가 $0.05c$에서 $0.3c$ 사이로 진화하며, 6 GHz에서 $10^{27}$–$10^{30}$ erg s$^{-1}$ Hz$^{-1}$ 범위의 라디오 광도를 생성할 수 있음을 보여준다. 광도 곡선의 형태는 충돌 유형과 환경 변수에 따라 크게 달라진다:

• 플레어+플레어 충돌: 이들은 $\sim 10$–$100$일의 상승 및 쇠퇴 시간 규모를 가진 빠르게 진화하는 광도 곡선을 생성한다. 피크는 일반적으로 SSA 광학적 두께가 1로 떨어질 때 결정된다. 이 모델들은 다중의 날카로운 라디오 피크를 재현할 수 있다.
• 플레어+CNM 충돌 (상수 밀도, $\alpha=0$): 이 모델들은 $\sim 10^2$–$10^3$일 동안 지속될 수 있는 광도 곡선을 생성한다. 피크 시간 규모는 일반적으로 세돕 시간(Sedov time, 쓸어 담은 질량이 플레어 질량과 같아지는 시점)에 의해 결정된다. 높은 CNM 밀도는 더 빠른 감속과 더 짧은 피크 시간 규모를 초래한다.
• 플레어+CNM 충돌 (거듭제곱 법칙, $\alpha=2.5$): 이 모델들은 중간적인 거동을 보인다. 낮은 밀도에서는 더 일찍 피크에 도달하고 더 희미하며, 높은 밀도에서는 더 오래 지속되고 더 밝다. 피크 시간 규모는 Flare+Flare 케이스와 유사하게 SSA 광학적 두께에 의해 결정되는 경우가 많다.

저자들은 모델 예측치를 지연된 라디오 방출을 보이는 선택된 관측된 TDE들과 비교하였다:

• PS16dtm: $M_{fl}=0.1 M_\odot$ 및 $n_0=30$ cm$^{-3}$인 Flare+$\alpha=0$ CNM 모델로 잘 재현된다.
• ASASSN-15oi: 두 번째 라디오 피크는 Flare+$\alpha=0$ CNM 모델($M_{fl}=0.1 M_\odot$, $n_0=10^3$ cm$^{-3}$)과 일치하지만, 첫 번째 피크는 다른 메커니즘(예: 스트림-스트림 상호작용)을 필요로 할 가능성이 높다.
• AT2020vwl: 두 번째 피크는 Flare+$\alpha=0$ CNM 모델($M_{fl}=0.01 M_\odot$, $n_0=10^2$ cm$^{-3}$)과 일치한다.
• AT2024tvd: 두 개의 날카롭고 빠르게 진화하는 피크는 Flare+Flare 모델에 매우 적합하며, 이는 다중 방출 에피소드를 시사한다.
• AT2019azh 및 AT2018zr: 이 사건들은 혼합된 특성을 보이며, 단일 모델이 상승 및 쇠퇴 시간 규모를 동시에 맞추는 데 어려움을 겪는다. 이는 복잡한 기저 물리 또는 다중 구성 요소를 시사한다.

의의 및 주장
본 논문은 디스크 불안정성에 의해 구동되는 유출물의 충돌이 상당수의 TDE에서 관측되는 지연된 라디오 방출을 일으키는 유효한 메커니즘이라고 주장한다. 구체적으로:

1. 메커니즘 검증: 모델은 여러 관측된 TDE의 광도 곡선과 분광 에너지 분포(SED)를 성공적으로 재현하며, 이는 지연된 약상대론적 유출물이 강착 디스크의 불안정성 또는 상태 전이에서 기원한다는 가설을 뒷받침한다.
2. 시나리오 차별화: 본 연구는 광도 곡선 진화를 통해 서로 다른 물리적 기원을 구분할 수 있는 프레임워크를 제공한다. Flare+Flare 충돌은 독보적으로 밝고 빠르게 진화하는(시간 규모 $\sim 100$일) 다중 피크를 설명할 수 있는 반면, Flare+CNM 충돌(특히 상수 밀도)은 더 느리고 장기간 지속되는 쇠퇴($\sim 10^3$일)를 더 잘 설명한다.
3. 관측적 제약: 저자들은 모델이 많은 사건을 잘 맞추지만, 일부 사건(예: ASASSN-14ae, AT2018hyz)의 SED(예: 전자 스펙트럼 지수 $p$) 및 광도 곡선 상승 시간에서의 불일치는 모든 지연된 방출이 이 단일 메커니즘으로 설명될 수는 없음을 시사하며, 추가적인 물리(예: 역 컴프턴 냉각 또는 다른 밀도 프로파일)가 필요할 수 있음을 언급한다.

본 연구는 디스크 불안정성에 의한 유출물 충돌이 특히 다중 피크를 가진 특정 클래스의 지연된 라디오 방출을 설명하는 데 유망한 설명을 제공하지만, 유출 기원의 다양성을 완전히 규명하기 위해서는 추가적인 후기 라디오 및 다파장 모니터링이 필요하다고 결론짓는다.