[论文解读] FORGE'd in FIRE: Resolving the End of Star Formation and Structure of AGN Accretion Disks from Cosmological Initial Conditions

Philip F. Hopkins, Michael Y Grudić+9 more|arXiv (Cornell University)|Sep 22, 2023

Astrophysical Phenomena and Observations被引用 4

一句话总结

本研究首次实现了将FIRE尺度的星系物理与STARFORGE尺度的恒星形成和反馈自洽耦合的宇宙学辐射磁流体动力学模拟，解析度达到超大质量黑洞周围<100 au。结果表明，引力扭矩和强磁场抑制了亚秒差距尺度内的恒星形成，从而通过由磁压稳定、磁冻结的引力湍流结构，实现了10–100 M⊙ yr⁻¹的持续吸积率。

ABSTRACT

It has recently become possible to zoom-in from cosmological to sub-pc scales in galaxy simulations to follow accretion onto supermassive black holes (SMBHs). However, at some point the approximations used on ISM scales (e.g. optically-thin cooling and stellar-population-integrated star formation [SF] and feedback [FB]) break down. We therefore present the first cosmological radiation-magnetohydrodynamic (RMHD) simulation which self-consistently combines the FIRE physics (relevant on galactic/ISM scales where SF/FB are ensemble-averaged) and STARFORGE physics (relevant on small scales where we track individual (proto)stellar formation and evolution), together with explicit RMHD (including non-ideal MHD and multi-band M1-RHD) which self-consistently treats both optically-thick and thin regimes. This allows us to span scales from ~100 Mpc down to &lt;100 au (~300 Schwarzschild radii) around a SMBH at a time where it accretes as a bright quasar, in a single simulation. We show that accretion rates up to $\sim 10-100\,{ m M_{\odot}\,yr^{-1}}$ can be sustained into the accretion disk at $\ll 10^{3}\,R_{ m schw}$, with gravitational torques between stars and gas dominating on sub-kpc scales until star formation is shut down on sub-pc scales by a combination of optical depth to cooling and strong magnetic fields. There is an intermediate-scale, flux-frozen disk which is gravitoturbulent and stabilized by magnetic pressure sustaining strong turbulence and inflow with persistent spiral modes. In this paper we focus on how gas gets into the small-scale disk, and how star formation is efficiently suppressed.

研究动机与目标

• 解析星系模拟中从宇宙学尺度到亚秒差距尺度的过渡，其中标准ISM尺度近似失效。
• 研究在类星体阶段实现超大质量黑洞持续高吸积率的物理机制。
• 理解活跃星系核附近极端条件下恒星形成被抑制的原因。
• 为未来高分辨率吸积盘模拟提供自洽的初始和边界条件。

提出的方法

• 采用自适应网格加密技术，从~100 Mpc的宇宙学尺度逐步放大至超大质量黑洞周围<100 au的尺度，进行zoom-in宇宙学模拟。
• 自洽耦合FIRE物理模型（在星系/ISM尺度）与STARFORGE物理模型（在恒星/亚秒差距尺度），以实现恒星形成与反馈的统一描述。
• 引入非理想磁流体动力学（NIMHD），包括双向扩散、欧姆电阻率和各向异性热传导，以模拟光学厚与光学薄区域。
• 采用多波段M1辐射流体动力学求解器，包含详细的热化学过程、尘埃消光截面和宇宙射线输运。
• 实现一个可动态演化的超大质量黑洞吸积粒子，其吸积过程可追踪至300个史瓦西半径以内。
• 采用自适应加密技术，将气体和恒星结构的解析度提升至<1 M⊙，并保持重力与流体动力学求解器的一致性。

实验结果

研究问题

• RQ1 气体如何从星系尺度向类星体宿主星系的最内区运输并吸积？
• RQ2 在超大质量黑洞周围亚秒差距环境内，何种物理机制抑制了恒星形成？
• RQ3 磁场与引力扭矩如何塑造亚秒差距尺度吸积盘的结构与稳定性？
• RQ4 当同时解析大尺度星系物理与小尺度恒星反馈时，自洽的吸积率与盘面特性为何？
• RQ5 辐射、磁场与冷却消光截面如何协同调控湍流、恒星形成气体向稳定、吸积盘的转变？

主要发现

• 尽管处于极端条件，10–100 M⊙ yr⁻¹的吸积率仍能持续维持在距离<10³个史瓦西半径的吸积盘内。
• 由于冷却光学厚度极高，且环向磁场强度使磁临界质量超过盘面质量，亚秒差距尺度内的恒星形成几乎完全被抑制。
• 在中等尺度形成引力湍流、磁冻结的盘结构，由磁压稳定，持续维持螺旋模态与内流。
• 在亚千秒差距尺度，恒星与气体之间的引力扭矩占主导地位，是恒星形成停止前角动量传输的主要驱动力。
• 盘结构由磁压稳定，非理想MHD效应在小尺度维持湍流与内流中起关键作用。
• 本模拟为未来吸积盘模拟提供了可预测的内边界条件，揭示了此前未探索的强磁场、湍流吸积盘区域。