Astronomy:cosmology:simulation

Sub-resolution modelling of astrophysical processes 在宇宙学流体动力学模拟中，涉及星系形成和演化的天体物理过程跨越了巨大的动力学范围。它们从约$\sim$Mpc尺度——引力不稳定性驱动暗物质组分演化与暗晕等级式并合之处——一直延伸到约$\sim$pc和$\sim$亚pc尺度——黑洞吸积和恒星形成等过程发生之处，途经约$\sim$kpc尺度——例如星系风将恒星反馈能量分配到周围介质之处。 因此，对次网格物理的需求至关重要：次分辨率模型描述的是发生在宇宙学流体动力学模拟分辨率极限以下的过程，但这些过程会影响在显式解析尺度上演化的模拟结构。次分辨率方案通常借助相当简单的解析或理论模型，辅以恰当选择并经校准以重现观测的参数，对相当复杂的过程进行唯象描述（另见§3.16小节）。The astrophysical processes involved in the formation and evolution of galaxies in cosmological hydrodynamical simulations span a huge dynamical range of scales. They indeed vary from the $\sim$Mpc scales where gravitational instabilities drive the evolution of the DM component and the hierarchical assembly of haloes, down to the $\sim$parsec and $\sim$sub-parsec scales, where processes like BH accretion and star formation take place, going through the $\sim$kpc scales, where e.g. galactic winds distribute the stellar feedback energy to the ambient medium. The call for sub-grid physics is thus essential: sub-resolution models account for processes that occur below the resolution limit of cosmological hydrodynamical simulations, but that affect the evolution of the simulated structure on scales that are explicitly resolved. Sub-resolution prescriptions usually resort to rather simple analytical or theoretical models, and/or to the phenomenological description of rather complex processes, through a suitable choice of parameters that are calibrated to reproduce observations (see also sub-section §3.16). ...

Gas cooling 热气体冷却出现在方程equation:firstlaw的右端，同时也出现在描述亚网格模型中可能存在的不同气体相之间质量与能量流动的方程中。Hot gas cooling enters the right-hand side of equation equation:firstlaw, as well as the equations describing mass and energy flows among different gas phases which may be present in sub-resolution models. 在宇宙学模拟中，关注的重点通常是维里温度超过$\sim 10^4$ K的结构。冷却函数$\Lambda(u,\rho)$标准实现中的常见假设是气体为光学薄且处于电离平衡。 通常还假设三体冷却过程不重要，从而将处理限定在二体过程。对于具有H和He原初组分的等离子体，这些过程包括：H$^0$和He$^+$的碰撞激发，H$^0$、He$^0$和He$^{+}$的碰撞电离，H$^+$、He$^+$和He$^{++}$的标准复合，He$^+$的双电子复合，以及自由-自由发射（轫致辐射）。 碰撞电离和复合速率仅依赖于温度。因此，若不存在电离背景辐射，所得速率方程可解析求解，由此得到的冷却函数$\Lambda(u)/\rho^2$如Fig. 9左图所示。 反之，若存在电离背景辐射，速率方程则可通过迭代求解。注意，对于典型的宇宙学辐射背景（例如，来自恒星形成星系和类星体的紫外背景，见[118] 、[119] 、[120] 、[121] ），冷却函数的形状可能发生显著改变，特别是在低密度下。更详细的讨论见星系形成章节和[17] 。In cosmological simulations, the focus is usually on structures whose virial temperature exceeds $\sim 10^4$ K. Common assumptions in standard implementations of the cooling function $\Lambda(u,\rho)$ is that the gas is optically thin and in ionization equilibrium. It is also usually assumed that three-body cooling processes are unimportant so as to restrict the treatment to two-body processes. For a plasma with primordial composition of H and He, these processes are: collisional excitation of H$^0$ and He$^+$, collisional ionization of H$^0$, He$^0$, and He$^{+}$, standard recombination of H$^+$, He$^+$, and He$^{++}$, dielectric recombination of He$^+$, and free--free emission (bremsstrahlung). The collisional ionization and recombination rates depend only on temperature. Therefore, should an ionizing background radiation be absent, the resulting rate equations can be solved analytically. This leads to a cooling function $\Lambda(u)/\rho^2$ as illustrated in the left panel of Fig. 9. On the other hand, in the presence of ionizing background radiation, the rate equations can be solved iteratively. Note that for a typical cosmological radiation background (e.g., UV background from star-forming galaxies and quasars, see [118] ,[119] ,[120] ,[121] ), the shape of the cooling function can be significantly altered, especially at low densities. For a more detailed discussion, see, for example, 星系形成章节 and [17] . ...

Star formation 恒星形成是一个关键的天体物理过程，涉及一系列层次递进的物理过程：来自大尺度（$\sim$Mpc）的气体吸积、气体冷却产生中性氢（HI）和分子氢（H${2}$）（$\sim$kpc尺度）、巨分子云（GMCs；$\sim$10$\div$100 pc）的形成、H${2}$ 碎裂并吸积为团块（$\sim$1 pc）和核心（$\sim$0.1 pc），以及随后核心收缩形成恒星（$\sim$10$^{-8}$ pc）。大量综述对这些复杂过程以及多相恒星形成 ISM 的性质和不同恒星形成模式的描述进行了详尽的讨论，例如[127] 、[128] 、[129] 、[130] 。 恒星形成是致密冷气体云引力坍缩的结果。由于恒星形成实际发生的亚pc尺度无法为最先进的宇宙学流体动力学模拟所分辨，因此这是一个亚分辨率过程。从宇宙学模拟所涉及的质量和尺度的数量级来估计，气体元素的质量通常在 $10^4 \div 10^7$ M$_{\odot}$ 范围内（它们无法分辨GMCs），其相关的软化长度也与数百pc相当或更大（取决于分辨率）。恒星粒子在质量和力分辨率上与气体元素相当，并被处理为SSPs。总而言之，宇宙学模拟中的恒星形成旨在有效捕捉经历了辐射冷却的冷气体向无碰撞恒星粒子的转化。 若满足若干判据，气体粒子或网格便可形成恒星：流动必须局部汇聚（$\nabla \cdot \mathbf{v} < 0$）且Jeans不稳定（即声波穿越时标超过动力学时标）；气体密度需超过临界密度阈值（以防止在极高红移的低密度气体中形成恒星）；以及氢原子的局部数密度需超过阈值数密度（见下文）。上述条件的详细信息见星系形成章节（另见[17] ）。 [131] 的奠基性工作推动了宇宙学模拟中对恒星形成 ISM 的描述和恒星形成数值建模的发展。他们引入了一个多相 ISM 模型，该模型考虑了恒星形成和恒星反馈。该模型认为，热气体和冷气体成分在致密气体粒子内以压强平衡的方式共存，即所谓的多相粒子。由于分子云无法被分辨，冷相和热相的性质被假定为代表 ISM 小体积上的平均值。多相粒子内的冷相为恒星形成提供了储库，同时通过热气体冷却得到补充。假设驱动恒星形成的冷气体中有一部分被瞬时用于加热热相并蒸发一部分冷相，以此模拟短寿命大质量恒星以超新星形式爆发的过程。通过积分描述热气体和冷气体相能量演化的方程，便可计算出恒星形成率。模型的参数经过校准以重现 Schmidt-Kennicutt 关系[132] 、[133] 。 该模型假设冷气体成分具有恒定温度（$\sim 10^3$ K）。此外，它以宁静、自调节的恒星形成为特征，因为冷却作用对冷相的补充被超新星加热导致的蒸发所平衡。在此模型中，ISM 可用有效状态方程来描述：两相气体粒子的有效压强具有由热相和冷相性质所约束的函数形式，并且由于自调节恒星形成，预期在时间上保持恒定。 在最先进的宇宙学模拟中广泛采用的若干亚分辨率模型都受到上述模型的启发，或大致基于该模型。 该模型的一个有趣变体认为，气体粒子能够形成恒星的密度阈值依赖于金属丰度，其假设是气体金属丰度越高，暖中性气体转变为冷分子气体所需的密度便越低[134] 、[41] 。 Springel 和 Hernquist 模型的一个关键特征是采用恒星形成气体的有效状态方程。该思想源于这样一种认识：调节 ISM 的小尺度效应（例如热不稳定性、湍流、热传导）在短时标内建立了 ISM 的自调节平衡态。在类似的机制下，ISM 的平均温度或能量可近似为仅依赖于暖/冷气体密度的函数（不依赖状态方程或不依赖多相分辨率元素的模拟示例分别见[135] 、[136] 和[137] ）。 如上所述，宇宙学背景下大多数结构形成的数值模拟仍远不能分辨 GMCs 的结构，也远不能从第一性原理出发建模恒星形成。它们反而依赖于各种将 SFR 与可用恒星形成气体联系起来的亚分辨率处方。这通常体现为考虑冷而致密的气体（$T \lesssim \mbox{a few} \times 10^4$ K；$n\gtrsim 0.1$-$10$ cm$^{-3}$），并假设局部 SFR 正比于其密度除以一个时标（该时标接近冷气体的动力学时标）。上述用于选择恒星形成气体的密度和温度阈值实际上反映了 ISM 中 HI 的物理性质。而 SFR 与密度的正比关系则大致受 Schmidt-Kennicutt 定律[132] 、[133] 的启发。 采用这一方法的最先进宇宙学模拟包括 Illustris[40] 、Magneticum[138] 、NIHAO[137] 、Illustris-TNG[139] 、Massive-Black[39] 和Fable[140] 模拟、Horizon-AGN系列[141] 、[142] 、DIANOGA模拟[143] 、[144] ，以及SLOW模拟[145] 。 在通过[146] 的理论模型估算分子气体丰度、进而推算 SFR 的宇宙学模拟中，值得提及的有 MUFASA 和 SIMBA 模拟[147] 、[44] 、FIRE 系列[148] 、[149] 以及 FIREBOX[150] 。 其他对 H$_{2}$ 进行亚分辨率处理的宇宙学模拟示例包括[151] 、[122] 、[152] 、[153] 、[154] 、[155] 、[156] 、[157] 、[158] 、[159] 、[160] 。由于上述模拟针对单个暗晕或较小体积，它们能够对恒星形成分子气体实现更精细的建模。 接下来，我们简要概述最先进宇宙学模拟为选择恒星形成气体和建模恒星形成而采用的不同处方。下列条件通常作为本小节开头所述条件的补充。 Illustris、Magneticum、Illustris-TNG 和 SLOW 模拟均基于一个密度阈值（$n_{\ast} \simeq 0.1$ cm$^{-3}$）来选择恒星形成气体，并假设一个恒定的耗散时标将冷气体转化为恒星并估算 SFR（基于[131] ）。 NIHAO 模拟以及 ERIS 和 GIGA-ERIS 模拟[161] 、[162] 则同时采用密度和温度阈值来选择恒星形成气体。 依赖 MUPPI 亚分辨率模型（例如[163] 、[155] 、[135] 、[164] 、[165] 、[166] 、[160] ）选择多相气体的模拟也采用了密度和温度阈值：恒星形成随后从冷气体中的分子部分进行，分子比例通过[167] 或[168] 的处方估算；最后，耗散时标取为冷气体的动力学时标。MUPPI 中更先进的建模利用了对尘埃形成和演化的处理[165] 、[166] 来预测 H$_{2}$ 的演化[169] ，进而用于估算局部 SFR。 类似地，MUFASA 和 SIMBA 模拟采用基于 H$_{2}$ 的恒星形成处方，遵循[146] 公式（认为密度超过阈值的气体是恒星形成的），并假设耗散时标为局部动力学时标。FIRE 模拟中也采用了类似的假设来估算冷相中的分子气体。 EAGLE 和 FLAMINGO 模拟系列[41] 、[6] 假设了一个金属丰度依赖的恒星形成阈值[134] 以及一个温度阈值：为了将恒星形成与分子气体联系起来，它们利用 Cloudy 辐射转移代码[124] 来确定从暖原子气体相到冷中性气体相的转变。在这些模拟中，SFR 依赖于压强而非密度。 有趣的是，HORIZON 模拟系列（例如[142] ）——它从高于给定密度阈值的气体网格中选择恒星形成气体——具有局部变化的恒星形成效率，这与几乎所有上述模拟不同，后者假设冷/分子气体转化为恒星粒子的比例是恒定的。 [170] 、[171] 讨论了选择恒星形成气体的不同判据如何对模拟星系产生强烈影响。 恒星形成通常根据随机模型来实现：该过程不是对每颗恒星分别描述，而是以与实际 SFR 一致的期望值随机进行。例如，若概率 $p = m_{gas}/m_{\star} (1- {exp}(m_{\star}/m_{gas}))$ 超过一个随机数，则生成一个新的恒星粒子[131] 。 作为改进，可以假设存在若干恒星代[172] ，这样新生成的恒星粒子以其来源气体元素质量的一定份额作为它们的初始质量。 系列导航 ← 上一篇：§3.11 → 下一篇：§3.13 ...

Stellar evolution and chemical enrichment 化学演化是宇宙学流体动力学模拟中星系演化的自然结果。化学演化模型已纳入宇宙结构形成的宇宙学模拟中，以恰当地研究化学增丰过程：星际介质（ISM）和星系周介质（CGM）中金属的分布记录了关于过去恒星形成和反馈历史的宝贵信息，是流体动力学模拟的一项关键特征（例如 [173] ,[174] ）。与明确分辨单颗恒星的观测不同，宇宙学模拟通过星粒子对星族进行粗略采样。大多数最先进的宇宙学流体动力学模拟（涵盖大宇宙体积和单个星系）中的星粒子质量，取决于分辨率，通常介于 $10^8$ 到 $10^3$ M$_{\odot}$ 之间（例如 [36] ,[40] ,[39] ,[175] ,[41] ,[176] ,[139] ,[6] ,[177] ,[178] ,[179] ,[180] ,[155] ,[135] ,[181] ）。在这些模拟中，星粒子是分辨率元素：每个星粒子代表一个简单星族（SSP），即一组具有相同初始金属丰度的同时代恒星的集合。每个星粒子的初始化学组分与其母气体元素相同，并由初始质量函数（IMF）表征。Chemical evolution is a natural outcome of galaxy evolution in cosmological hydrodynamical simulations. Models of chemical evolution have been included in cosmological simulations of cosmic structure formation to properly address the study of the chemical enrichment process: the distribution of metals in the ISM and in the circum-galactic medium (CGM) encodes valuable information of the past history of star formation and feedback, and is a crucial feature of hydrodynamical simulations (e.g. [173] ,[174] ). At variance with observations that explicitly resolve individual stars, cosmological simulations provide a coarse sampling of stellar populations by means of star particles. The majority of state-of-the-art cosmological hydrodynamical simulations of both large cosmological volumes and individual galaxies have star particles whose mass typically ranges between $10^8$ and $10^3$ M$_{\odot}$ according to resolution (e.g. [36] ,[40] ,[39] ,[175] ,[41] ,[176] ,[139] ,[6] ,[177] ,[178] ,[179] ,[180] ,[155] ,[135] ,[181] ). In these simulations, star particles are resolution elements: each of them represents a simple stellar population (SSP), i.e. an ensemble of coeval stars that share the same initial metallicity. Every stellar particle initially shares the chemical composition of the gas element from which it has been originated, and is characterized by an initial mass function (IMF). ...

Stellar feedback 恒星反馈是结构形成宇宙学模拟的关键组成部分。在其众多关键作用中，它可防止高红移处的过度冷却，在星际介质中分配能量，将金属驱出恒星形成区，并触发星系外流，从而确保形成中的星系与其周围星系周介质之间的持续相互作用。 鉴于驱动和发射星系风的pc尺度物理过程尚不完全清楚，且远未达到能在宇宙学模拟中直接实现的程度，这些模拟不得不借助唯象预设来捕捉星系外流的效应。从驱动外流并支配其运动学的能量来源来看，有两种常见途径：星系风可以是能量驱动的，也可以是动量驱动的。关于这两种情景的细节，建议读者参考星系形成章节，此文为本小节的具体主题提供了补充信息。 本概述远非完备：更深入的讨论可见于[213] 、[155] 、[214] 、[215] 等文献。Stellar feedback is a key component of cosmological simulations of structure formation. Among its many crucial roles, it prevents overcooling at high redshift, distributes energy in the ISM, drives metals out of the star formation sites, and triggers galactic outflows, which guarantee a continuous interaction between the forming galaxy and its surrounding CGM. Since the understanding of the pc-scale physics responsible for launching and driving winds is still partial and however far from being implemented directly in cosmological simulations, these simulations have to resort to phenomenological prescriptions to capture the effects of galactic outflows. As for the sources of energy that power outflows and govern their kinematic, there are two commonly pursued approaches: winds can be either energy-driven or momentum-driven. For details about these two possible scenarios, we refer the reader to 星系形成章节, which provides complementary information on the specific topic of this sub-section. This outline is far from being complete: deeper insight can be gained e.g. from [213] ,[155] ,[214] ,[215] . ...

Supermassive black holes in cosmological simulations 观测表明，几乎每个星系在其最内层区域都寄宿着一个超大质量（$10^8 \div 10^{10}$ M$_{\odot}$）黑洞（SMBH）（例如，[231] ,[232] ,[233] ,[234] ,[235] ,[236] ,[237] ,[238] ,[239] ）。其中一部分黑洞表现出持续的活动性，称为活动星系核（AGN）。过去和当前活动的证据可见于椭圆星系、星系群和星系团的X射线图像，其中AGN的印记通常表现为凹陷和涟漪。在众多例子中，一个众所周知的代表是MS 0735+7421星系团的复合（X射线、射电和光学）图像（例如[240] ）：巨大的X射线空洞充满射电辐射，并被Chandra图像中清晰可见的椭圆形茧状激波包围。Observations suggest that almost every galaxy hosts a supermassive ($10^8 \div 10^{10}$ M$_{\odot}$) black hole (SMBH) in its innermost regions (e.g., [231] ,[232] ,[233] ,[234] ,[235] ,[236] ,[237] ,[238] ,[239] ). A fraction of these BHs exhibits ongoing activity and is called AGN (active galactic nuclei). Evidence of past and ongoing activity is observed for instance in the X-ray images of elliptical galaxies and galaxy groups and clusters, where AGN imprints often appear as depressions and ripples. Among many, a well-known example is represented by the composite (X-ray, radio, and visual) image of the MS 0735+7421 galaxy cluster (e.g. [240] ): giant X-ray cavities are filled with radio emission, and surrounded by a cocoon shock clearly visible in the Chandra image as an elliptical edge. ...

AGN feedback in state-of-the-art cosmological simulations 观测使我们得以理解AGN反馈如何以不同的机制（如膨胀气泡或发射外流）和不同的表现形式，在各种系统中发展。此外，这是一个反复发生的过程，每一次AGN爆发都在系统中留下清晰的印记。星系内部及周围多相气体的存在，进一步增加了AGN反馈理解和建模的复杂性。多波长观测揭示星系中存在跨越广泛密度、温度和电离状态范围的气体。多相气体不仅存在于富含冷气体的旋涡星系中，也存在于椭圆星系以及星系群和星系团的最内部区域——这些环境以X射线辐射的热气体为主导。Observations allow us to appreciate how AGN feedback develops with different mechanisms (for instance inflating bubbles or launching outflows), with different appearances, in a variety of systems. Besides, it is a recurrent process, with each AGN burst leaving a clear signature in the system. Additional evidence which adds complexity to the comprehension and modelling of AGN feedback is the presence of multiphase gas in and around galaxies. Multiwavelength observations reveal the presence of gas spanning a wide range of densities, temperatures, and ionisation states in galaxies. Multiphase gas is present not only in spiral galaxies, which are systems rich in cold gas, but also in ellipticals and in the innermost regions of galaxy groups and clusters, which are environments known to be dominated by X-ray emitting, hot gas. ...

Current state and perspectives 即使在今天，用作超大型巡天（如EUCLID、DESI、LSST）理论对照的宇宙学模拟，由于需要覆盖极大的体积，仍然基于纯引力物理（例如，[10] ,[14] ）。这些模拟通常辅以星系形成的半解析模型（SAMs）（例如，[314] ）。虽然SAMs提供了星系族群性质的逼真描述，但它们充其量只能给出ISM/IGM性质的间接信息。事实上，它们并未自洽地处理气体动力学，也只能粗略捕捉重子对结构形成的影响，而这对环境效应的研究至关重要。Even nowadays, the cosmological simulations used as theoretical counterparts for very large surveys (like EUCLID, DESI, LSST), for which extremely large volumes need to be sampled, are based on pure gravitational physics (e.g., [10] ,[14] ). These simulations are usually complemented by running semianalytic models (SAMs) of galaxy formation (e.g., [314] ). While SAMs provide a realistic description of the properties of galaxy populations, they bring at best indirect information on the properties of the ISM/IGM. In fact, they do not include a self-consistent treatment of gas dynamics and loosely capture the effects of baryons on structure formation, which is highly relevant for the study of environmental effects. ...