Rates of tidal disruption events from constrained cosmological simulations of the local Universe: population properties and implications for transient surveys

기술 요약: 국부 우주의 제약된 우주론적 시뮬레이션으로부터 도출된 조석 파괴 사건(TDE) 발생률

문제 제로 (Problem Statement)
조석 파괴 사건(TDE)은 정적 상태의 초거대 블랙홀(SMBH) 집단을 연구하기 위한 핵심적인 탐사 도구이다. 그러나 이론적 예측과 관측적 측정치 사이에는 상당한 불일치가 존재한다. 이상적인 이체 이완(two-body relaxation) 모델은 $\gtrsim 10^{-4} \text{ yr}^{-1} \text{ gal}^{-1}$ 이상의 발생률을 예측하는 반면, 광학, X선 및 적외선 조사로부터 얻은 관측적 제약 조건은 일반적으로 이보다 한 차수 낮은 ($\sim 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$) 값을 나타낸다. 또한, 전통적인 해석적 접근 방식은 2D 외삽에 의존하거나 은하 밀도 프로파일(예: 특이 등온 구체)에 대한 단순화된 가정을 사용하는 경우가 많아, 우주 구조 형성, 병합 이력, 그리고 호스트 은하의 구체적인 형태학적 다양성(cuspy vs. cored) 사이의 복잡한 상호작용을 설명하지 못한다. 베라 루빈 천문대의 LSST(Legacy Survey of Space and Time)가 다가옴에 따라, 다가오는 과도 현상 데이터를 정확하게 해석하기 위해 공간적으로 해상된 정밀한 TDE 발생률 추정치가 요구된다.

방법론 (Methodology)
저자들은 국부 우주를 대상으로 한 최첨단 제약된 우주론적 시뮬레이션인 "Simulation of the LOcal Web (SLOW)"을 활용하였으며, 특히 코마(Coma), 헤라클레스(Hercules), 샤플리(Shapley), 처녀(Virgo), 페르세우스(Perseus) 초은하군과 포락스(Fornax) 은하단을 포함하는 6개의 고해외상(high-resolution) "줌인(zoom-in)" 영역에 집중하였다. 이 영역들은 55–92 Mpc의 방사형 범위를 가지며, 그룹 규모에서 거대 초은하단 규모($M_{500} \approx 3.5 \times 10^{13} - 1.7 \times 10^{15} M_\odot$)에 이르는 환경을 포괄한다.

방법론은 다음과 같은 엄격한 필터링 및 계산 단계를 거친다:

1. 블랙홀 추출: 시뮬레이션 볼륨에서 SMBH를 추출하고, 그 질량($M_\bullet$)과 무차원 스핀 매개변수($\tilde{a}$)를 기록한다.
2. 주계열 필터링 (Main-Sequence Filtering): 손실 원(loss-cone) 역학 계산의 타당성을 보장하기 위해, 엄격한 필터링 체계를 통해 역학적으로 안정적인 "주계열" 블랙홀만을 분리한다. 이는 다음 네 가지 기준을 통해 결합되지 않았거나, 회전이 지배적이거나, 혹은 조석적으로 박리된 시스템을 제외한다:
   • 오프셋 컷 (Offset Cut): BH는 서브할로 중심으로부터 1.5 kpc 이내에 있어야 한다.
   • 속도 컷 (Velocity Cut): 상대 속도는 국부 항성 속도 분산보다 낮아야 한다.
   • 분율 이방성 (Fractional Anisotropy, FA) 컷: 국부 항성 환경은 $FA < 0.4$여야 하며, 이는 표준 손실 원 형식에 필요한 근사적인 등방성 압력 지지 구체를 보장한다.
   • 암흑 물질 비율 컷: 박리된 핵을 제외하기 위해 암흑 물질 대 항성 질량 비율이 $> 1$이어야 한다.
3. 형태학적 분류: 1 kpc 이내에서 측정된 3D 항성 밀도 기울기를 기반으로, 호스트 은하를 "cuspy" ($g \ge 2$) 또는 "cored" ($g \le 1$)로 분류한다. 이 반지름은 수치적 연화(softening) 아티팩트를 피하면서 미해상 핵 구조의 대리 지표 역할을 하도록 선택되었다.
4. 상대론적 효율성 보정: Kesden (2012)으로부터 유도된 상대론적 효율 계수($\eta$)를 적용한다. 이 계수는 블랙홀의 스핀과 힐스 질량 한계(Hills mass limit)를 고려하여, 손실 원의 일부가 직접 포획(direct capture)이 아닌 관측 가능한 TDE로 이어지는 비율을 정량화한다.
5. 발생률 계산: Stone & Metzger (2016) 프레임워크를 사용하여 체적 TDE 발생률($\Gamma_{vol}$)을 계산하며, 이때 상대론적 효율 $\eta$에 의해 가중치가 부여된 cusp 및 core 환경에 대한 별도의 베이스라인 발생률을 적용한다.

주요 기여 및 결과 (Key Contributions and Results)

• 체적 및 블랙홀당 발생률: 본 연구는 6개 환경 전체에 대해 평균 약 $\approx 600 \text{ Gpc}^{-3} \text{ yr}^{-1}$의 평균 체적 TDE 발생률을 도출하였다. 블랙홀당 평균 TDE 발생률은 $\approx 4.5 \times 10^{-5} \text{ yr}^{-1}$로 안정적인 것으로 나타났으며, 이는 현재의 관측 제약 조건과 일치하지만 이상적인 이론적 예측보다는 유의미하게 낮은 수치이다.
• 공간 분포 vs. 총 질량: 본 논문은 TDE 산출물이 블랙홀 인구 통계 및 공간적 집중도에 의해 결정되며, 총 클러스터 질량에 의해 결정되는 것이 아님을 입증한다. 절대적인 TDE 예산은 확장된 클러스터 헤일로 내의 cuspy 위성 은하들에 의해 압도적으로 지배되지만, 체적 발생률($\Gamma_{vol}$)은 블랙홀의 수 밀도($n_\bullet$)에 의해 크게 영향을 받는다.
• 코어 고갈 및 효율성: 중심부 코어 발생률은 역학적 고갈(병합이 저질량 서브할로를 파괴함) 및 직접 포획 제약(Hills limit)으로 인해 크게 감소한다. 결과적으로, 전체 TDE 예산은 거대한 중심 은하가 아니라 클러스터 외곽에 존재하는 저질량 cuspy 위성 은하 집단에 의해 주도된다.
• 환경 의존성:
  • 헤라클레스 (Hercules): 높은 블랙홀 수 밀도를 가진 거대하고 활발히 조립 중인 초은하군임에도 불구하고, 병합에 의한 블랙홀 질량 성장으로 인해 많은 BH가 $\eta \to 0$인 힐스 한계로 밀려나면서 가장 낮은 블랙홀당 효율($\tilde{\Gamma}$)을 보인다.
  • 포락스 (Fornax): 이 저질량 환경은 힐스 컷오프를 피하는 병합되지 않은 저질량 블랙홀의 비율이 불균형적으로 높기 때문에 블랙홀당 가장 높은 효율($\tilde{\Gamma} \approx 1.07$)을 보인다. 결정적으로, 포락스는 6개 환경 중 가장 낮은 블랙홀 수 밀도($n_\bullet$)와 가장 큰 반수 밀도 반지름($R_{1/2,\bullet} = 31.26$ Mpc)을 나타내며, 이는 블랙홀의 공간적 분포가 가장 분산되어 있음을 의미한다.
• 2D 외삽과의 불일치: 본 연구는 전통적인 2D 분석적 외삽이 중심 표면 밀도를 과대평가한다는 점을 강조한다. 시뮬레이션된 2D 표면 발생 밀도는 문헌에 보고된 2D 투영 값보다 약 1 차수 낮다. 이는 실제 TDE 예산이 중심핵이 아닌 확장된 헤일로 전체에 걸쳐 분포되어 있기 때문이다.

의의 (Significance)
본 논문은 국부적인 항성 환경을 해상하는 완전한 3D 우주론적 프레임워크를 활용함으로써 이론적 예측과 관측적 측정 사이의 간극을 메운다고 주장한다. 단순화된 2D 외삽과 이상적인 밀도 프로파일을 넘어섬으로써, 저자들은 복잡한 우주 구조의 조립 이력을 고려한 더 현실적인 "TDE 예산"을 제공한다.

연구 결과는 역학적 코어 고갈과 상대론적 포획 한계를 고려할 때, 관측된 TDE 발생률 $\sim 3 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$가 우주론적 시뮬레이션에 의해 견고하게 뒷받침됨을 보여준다. 도출된 $4.5 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$의 발생률은 시뮬레이션의 하위 격자 물리(subgrid physics)가 절대적 산출량을 지배하는 저질량 왜소 은하를 약간 과잉 생성할 수 있음을 인정하면서도, 향후 과도 현상 조사에 대한 견고한 이론적 상한선으로 제시된다. 이 연구는 TDE 발생률이 보편적인 상수가 아니라, 블랙홀의 공간적 집중도와 호스트 환경의 병합 이력과 같은 거시적 특성에 의해 근본적으로 결정된다는 점을 강조한다.