은하 규모 시뮬레이션을 위한 초에디턴 제한 블랙홀 서브그리드 모델

초록

본 논문은 초에디턴 제한 상태에서 회전하는 블랙홀(BH)의 질량·스핀 진화를 추적할 수 있는 새로운 서브그리드 모델을 제시한다. GIZMO 코드에 구현된 이 모델은 내부 광자 포획 영역과 외부 얇은 α‑디스크(세 구역)로 구성된 가상 디스크를 사용하고, 질량 유입률에 따라 Bardeen‑Petterson 정렬 효과와 두꺼운 디스크의 프리세션을 전환한다. 이상적인 원형 가스 원반과 구형 별 성분을 포함한 실험에서, 가스 유입이 디스크와 BH의 각운동량을 급격히 감소시킬 때 초에디턴 비율이 크게 상승함을 확인하였다. 각운동량 불일치가 초에디턴 유입을 촉발하는 핵심 메커니즘임을 강조하며, 고‑레드쉬프트 은하의 동적 환경에서 이러한 현상이 자연스럽게 발생할 가능성을 제시한다.

상세 분석

이 연구는 초에디턴 제한(𝑓_Edd≫1) 상황을 물리적으로 구현하기 위해 두 단계의 서브그리드 디스크 구조를 도입한다. 내부 영역(R_ISCO ≤ R < R_trap)은 방사선‑수송 시뮬레이션에서 도출된 광자 포획(Photon‑trapping) 관계식(Eq.10‑11)을 사용해 표면밀도 Σ와 각운동량 L을 계산하고, 외부 영역은 전통적인 얇은 α‑디스크 이론(K08)을 네 개의 방사압·가스압 지배 구역(a, b, c)으로 나눈다(Eq.12‑15). 각 구역마다 온도·압력·두께 H를 α‑파라미터(α=0.1)와 연관시켜 연속적인 구조를 보장한다.

질량 유입률 𝑓_Edd,16은 디스크 전체 질량 M_disc와 총 각운동량 J_disc을 입력으로 뉴턴‑라프슨 방법을 통해 역으로 구한다(Eq.20). 이는 기존의 Bondi‑HL 혹은 단순한 베타‑디스크 모델과 달리, 디스크 자체의 저장·전달 능력을 고려한 자기일관적 접근이다. 특히, f_Edd,16의 상한값을 명시함으로써 물리적으로 허용 가능한 초에디턴 한계를 제시한다.

스핀 진화는 두 가지 모드로 전환된다. (1) 𝑓_Edd,16 ≤ ˆf_Edd,16(≈0.3)에서는 Bardeen‑Petterson 효과가 지배해 Lense‑Thirring 토크가 디스크를 BH 스핀 축에 맞추며, 토크식 Eq.36을 통해 스핀 방향·크기가 변화한다. (2) 𝑓_Edd,16 > ˆf_Edd,16에서는 내부 두꺼운 디스크가 프리세션을 일으키며, Eq.41에 기반한 프리세션 토크가 적용된다. 이 전환은 초에디턴 흐름이 디스크 두께를 급격히 증가시켜 수직 점성 ν_2가 얇은 디스크 가정에서 붕괴하기 때문이며, 모델이 이러한 비선형 전이를 자연스럽게 포착한다는 점이 혁신적이다.

시뮬레이션 설정은 중앙 BH(10^6 M_⊙)와 가스 원반(Σ∝R^−1)·구형 별 성분을 포함한다. 가스 유입을 의도적으로 비정렬(θ≈45°)시키면, 외부 가스가 디스크에 충돌하면서 디스크의 전체 각운동량이 급감하고, 이에 따라 R_warp이 축소되어 디스크가 BH와 거의 동조화된다. 결과적으로 𝑓_Edd,16이 10–100 수준까지 상승하고, 질량 성장률이 Eddington 제한을 크게 초과한다. 반면, 정렬된 경우에는 토크가 약해 초에디턴 비율이 𝑓_Edd≈1 이하에 머문다.

모델의 강점은 (i) 광자 포획 영역을 직접 포함해 실제 초에디턴 흐름의 효율을 재현, (ii) 스핀·각운동량 상호작용을 물리적으로 연결, (iii) GIZMO와 같은 현대 수치코드에 모듈형으로 삽입 가능하다는 점이다. 한계로는 (a) 1D 디스크 구조 가정이 강한 비축퇴성(예: 강한 자기장·극심한 비대칭) 상황을 완전히 포착하지 못함, (b) 방사선 피드백이 별도 서브그리드 모듈에 의존해 구현돼 실제 광자‑가스 상호작용을 단순화, (c) 파라미터 α·ξ 등에 대한 민감도 분석이 부족해 고‑레드쉬프트 은하 환경에 대한 일반화가 아직 제한적이다. 향후 3D 방사선‑MHD 시뮬레이션과의 교차 검증이 필요하다.