Milgromian 引力下恒星形成主序星系的模拟

我们利用自适应网格加密代码 \textsc{phantom of ramses} (\textsc{por})（一种带有 Milgromian 引力求解器的 \textsc{ramses} 代码改版）对恒星质量范围 $M_{\star}/M_\odot = 10^7 - 10^{11}$ 内的孤立盘星系进行了流体动力学 MOND 模拟。尺度长度和气体分数基于观测星系确定，模拟运行时间为 5~Gyr。主要目的是检验现有的恒星形成与恒星反馈子网格物理方案能否重现观测到的主序，并合理匹配刻画局部与全局恒星形成率如何与其他性质相关联的 Kennicutt-Schmidt 关系。模型中的恒星形成在初始化后不久即开始，并随模型演化持续进行。初始化后的星系确实演化到位于观测主序上的状态，并合理匹配了 Kennicutt-Schmidt 关系。因此，现有的子网格物理公式足够恰当，可得到大体上表现如同观测星系、半径增长且具有平坦转动曲线的星系——前提是我们使用 Milgromian 引力。此外，棒的强度往往与星系恒星质量呈反相关，而棒长与恒星质量强相关。无论质量如何，棒图样速度随时间保持恒定，表明动力学摩擦不影响棒动力学。这些模型展示了 Renzo 法则，并在大半径处形成结构，与真实星系非常相似。因此，在该框架下，重子物理已得到充分理解，不会对我们全局星系性质建模带来主要不确定性。