Bolometric correction factor and radiative efficiency for the super-Eddington accretion flow in tidal disruption events

기술 요약: 조석 파괴 이벤트(TDE)에서의 초에딩(Super-Eddington) 강착 흐름에 대한 볼로미터 보정 계수 및 복사 효율

문제 제기
볼로미터 광도(bolometric luminosity)와 복사 효율(radiative efficiency)을 추정하는 것은 거대 질량 블랙홀(BH) 주변의 강착 흐름을 이해하는 데 매우 중요하다. 조석 파괴 이벤트(TDE)에서, 관측으로부터 유도된 강착 질량($\Delta M$)과 태양형 별의 파괴에 대한 이론적 기대치인 $\sim 0.5 M_\odot$ 사이에는 상당한 불일치가 존재한다. 이 "잃어버린 에너지(missing energy)" 문제는 기존의 계산들이 흔히 일정한 X선 볼로미터 보정 계수($k_{bol} \approx 1$)와 높은 복사 효율($\eta = 0.1$)을 가정한다는 점 때문에 발생한다. 이전 연구들은 이러한 가정이 특히 X선보다는 관측 불가능한 극자외선(EUV) 대역이 방출을 주도하고, 이론적으로 강착 흐름의 내점(advection)으로 인해 복사 효율이 낮아질 것으로 예상되는 초기 초에딩 단계에서 결함이 있을 수 있음을 시사한다.

방법론
이러한 불확실성을 해결하기 위해, 저자들은 ATHELA++ 코드를 사용하여 2차원 축대칭 방사-유체 역학 시뮬레이션을 수행하였다. 이 연구는 $10^6, 10^{6.5}, 10^7 M_\odot$의 블랙홀에 의해 파괴된 포물선 궤도의 태양형 별($M_* = M_\odot, R_* = R_\odot$)에 의한 강착 흐름을 모델링했다.

• 시뮬레이션 설정: 저자들은 슈바르츠칠트 블랙홀 효과를 모사하기 위해 의사-뉴턴 퍼텐셜(pseudo-Newtonian potential)을 사용하였다. 질량은 원형화 반경($R_c = 2R_t$)에서 $\dot{M}_{inject} \propto t^{-5/3}$ 형태의 낙하율 프로파일을 따라 지속적으로 주입되었다. 계산 영역은 $2R_S \le r \le 10^5 R_S$이며 해상도는 $768 \times 256$이다.
• 후처리: $t = 16$일(안정적인 초에딩 강착이 확립된 시점으로 선택됨)에 저자들은 PYTHON을 이용한 몬테카를로 복사 전달 계산을 수행하여 방출 스펙트럼을 생성하였다.
• 매개변수: 스펙트럼은 시선 각도($\theta$)가 $0^\circ$(face-on)에서 $90^\circ$(edge-on)까지인 경우에 대해 계산되었다. X선 대역은 0.3–2 keV로 정의되었다.
• 계산: 볼로미터 보정 계수($k_{bol} = L_{bol}/L_X$)와 겉보기 복사 효율($\eta = L_{bol}/\dot{M}_{in}c^2$)을 도출하기 위해 스펙트럼으로부터 등방성 등가 X선 광도($L_X$)와 볼로미터 광도($L_{bol}$)를 적분하였다.

주요 결과

1. 질량 및 시선 각도에 따른 의존성: $k_{bol}$과 $\eta$는 모두 블랙홀 질량($M_{BH}$)과 시선 각도($\theta$)에 강하게 의존한다.
   • $k_{bol}$: 값은 수십에서 수천까지 범위가 넓다. $M_{BH} = 10^6 M_\odot$의 경우, $k_{bol}$은 $0^\circ$에서 $\sim 25$에서 $90^\circ$에서 $\sim 1000$까지 증가한다. 질량이 더 높은 경우($10^7 M_\odot$), $k_{bol}$은 모든 각도에 대해 $\sim 1000$ 근처에서 비교적 일정하게 유지된다. 이 값들은 흔히 채택되는 값인 1보다 실질적으로 훨씬 크다.
   • $\eta$: 도출된 복사 효율은 $\sim 10^{-3}$에서 $10^{-1}$ 사이이다. 구체적으로 $10^6 M_\odot$의 경우 $\eta$는 $\sim 0.001-0.02$이며, $10^7 M_\odot$의 경우 $\sim 0.13$에 도달하지만 대부분의 각도에서 일반적으로 0.1 미만이다. $\eta$는 시선 각도가 증가함에 따라 감소하며, 에딩도(Eddington-scaled) 강착률이 더 높은 낮은 블랙홀 질량에서 일반적으로 더 낮다.
2. 스펙트럼 특성: 방출 스펙트럼은 EUV 대역에 의해 지배된다. $\theta$가 증가함에 따라 X선 방출은 광학적으로 두꺼운 유출물(outflow)에 의해 점점 더 많이 차폐되고 재처리되어 광학 및 자외선 대역으로 변환된다. 결과적으로, $\theta$가 $0^\circ$에서 $90^\circ$로 변함에 따라 $L_X$는 수 차수(최대 $10^6 M_\odot$에서 $\sim 1000$배) 감소하는 반면, $L_{bol}$은 완만하게($\sim 15-20$배) 감소한다.
3. 보정 계수 비율: 이전의 추정치와 비교하여 추론된 강착 질량을 스케일링하는 비율인 $(k_{bol}/1.1) / (\eta/0.1)$은 시뮬레이션된 매개변수 전반에 걸쳐 $\sim 90$에서 $8 \times 10^4$까지의 범위를 갖는다. 이는 이전의 $\Delta M$ 추정치들이 2~5 자릿수만큼 과소평가되었을 가능성이 높음을 나타낸다.

관측된 TDE에 대한 적용
저자들은 도출된 $k_{bol}$ 및 $\eta$ 값을 몇몇 특정 X선 TDE(예: NGC 5905, NGC 3599, SDSS J132341.97+482701.3, ASASSN-14li)에 적용하였다. 적절한 $M_{BH}$와 $\theta$(일반적으로 TDE의 통합 모델과 일치하는 낮은 경사각)를 선택함으로써, 이 사건들에 대해 재계산된 강착 질량 $\Delta M$은 $0.25 - 0.5 M_\odot$ 범위 내에 들어온다. 이 결과는 태양형 주계열성의 파괴에 대한 이론적 기대치와 일치하며, 따라서 이러한 특정 사례들에 대한 "잃어버린 에너지" 문제를 해결하는 방안을 제시한다.

의의 및 한계
본 논문은 도출된 $k_{bol}$ 및 $\eta$ 값이 관측자의 시선 각도가 작을 때 TDE의 잃어버린 에너지 문제를 상당히 완화할 수 있는 물리적으로 동기화된 보정을 제공한다고 주장한다. 저자들은 X선으로 선택된 이벤트들이 낮은 경사각에서 관측된다는 TDE의 통합 모델과 이들의 해결책이 일치함을 강조한다.

그러나 저자들은 자신들의 연구 결과가 갖는 보편성에 대해 신중한 입장을 유지한다:

• 시선 각도 불확실성: 해결책은 $\theta$에 강력하게 의존하지만, $\theta$는 관측만으로는 정확하게 결정될 수 없다.
• 시간적 진화: 계산은 $t=16$일을 기준으로 한다. 저자들은 초에딩 단계가 수백 일 동안 안정적이라고 주장하지만, 강착률이 에딩 한계 아래로 떨어짐에 따라 $k_{bol}$과 $\eta$가 어떻게 장기적으로 진화하는지에 대해서는 추가 연구가 필요하다.
• 단순화: 본 연구는 태양형 별을 가정하고, 슈바르츠칠트 블랙홀(의사-뉴턴 퍼텐셜 사용)을 모델링하며, 블랙홀의 스핀과 자기장을 무시한다. 침투 계수(penetration factor), 항성 구조, 그리고 상대론적 효과와 같은 요소들은 테스트되지 않았으며 결과에 영향을 미칠 수 있다.

결론적으로, 이 연구는 X선 TDE와 같이 작은 시선 각도를 가진 경우에 적용될 수 있는, 관측된 복사 에너지와 이론적 질량 파괴 한계 사이를 일치시키는 구체적인 시뮬레이션 기반 매개변수 세트를 제공한다. 다만, TDE 물리학의 전체 매개변수 공간에 대한 완전한 해결을 위해서는 추가적인 조사가 필요함을 인정하고 있다.