[论文解读] N-MODY: a code for collisionless N-body simulations in modified Newtonian dynamics

P. Londrillo, Carlo Nipoti|ArXiv.org|Mar 31, 2008

Scientific Research and Discoveries被引用 3

一句话总结

N-MODY 是一种用于模拟在修正牛顿动力学（MOND）下无碰撞系统的球坐标系并行粒子-网格 N 体代码，通过球谐函数展开的松弛法求解非线性 MOND 场方程。结果表明，与含暗物质的牛顿引力相比，剧烈弛豫、并合及动力摩擦等动力学过程在 MOND 中表现出显著差异，弛豫与并合的时间尺度更长，而摩擦力更强。

ABSTRACT

We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter.

研究动机与目标

• 开发一种用于修正牛顿动力学（MOND）下无碰撞 N 体模拟的数值代码，作为解释星系运动学中暗物质的替代方案。
• 使用球谐函数展开的松弛法，在球坐标系中高效求解非线性 MOND 场方程。
• 实现对孤立恒星系统及任意静态密度分布的 MOND 模拟，提供一种用于对比 MOND 与牛顿动力学的工具。
• 研究基本动力学过程（如剧烈弛豫、相空间混合、并合及动力摩擦）在 MOND 中与含暗物质的牛顿引力下的行为差异。
• 提供一个公开可用的、支持串行与 MPI 并行执行的并行 FORTRAN 90 代码，适用于高性能计算。

提出的方法

• N-MODY 采用球坐标系下的粒子-网格方案，通过网格上的线性或二次形状函数重构密度场。
• 通过使用球谐函数展开的松弛法，在网格上求解非线性场方程 ∇·[μ(|∇ϕ|/a₀)∇ϕ] = 4πGρ 计算 MOND 势。
• 代码采用松弛参数（ω⁻¹ ≈ 0.3–0.5）和收敛容差（ε = 10⁻³），实现数值解，使 ∥M̂(ϕ)∥∞ ≈ 0.1 且 ∥M̂(ϕ)∥rms ≈ 0.01。
• 粒子加速度通过相同的形状函数从网格插值得到，并转换为笛卡尔分量以进行时间积分。
• 时间积分采用二阶或四阶蛙跳格式，确保模拟中能量与动量守恒。
• 代码通过 MPI 实现并行化，粒子在处理器间均匀分布，并对 θ 坐标应用域分解，实现部分势场求解器并行化。

实验结果

研究问题

• RQ1 与含暗物质的牛顿引力相比，MOND 中剧烈弛豫与相空间混合的时间尺度有何不同？
• RQ2 MOND 系统中动力摩擦的效率相对于含暗物质的牛顿系统如何？
• RQ3 在 MOND 与牛顿引力中，星系并合的时间尺度有何差异？
• RQ4 非线性 MOND 场方程在多大程度上改变了孤立恒星系统结构与演化，相比牛顿势？
• RQ5 基于球坐标系的粒子-网格代码能否高效且准确求解天体物理相关系统中的 MOND 场方程？

主要发现

• N-MODY 中的 MOND 势求解器在 ε = 10⁻³ 容差下收敛于 5–10 次迭代，数值精度达到 ∥M̂(ϕ)∥∞ ≈ 0.1 且 ∥M̂(ϕ)∥rms ≈ 0.01。
• 在 MOND 模拟中，当 N ≈ 10Ng 时，串行模式下势求解器占总 CPU 时间的约 10%，但在并行运行中因并行化程度有限而成为主要瓶颈。
• MOND 求解器的并行化效率约为 70%，计算成本大致与网格点数 Ng = (Nr+1)NθNφ 成正比。
• 即使使用 16 个处理器，当 N = 10⁷ 且 Ng = 128×64×128 时，N-MODY 中每个粒子每时间步的总 CPU 时间仍低于 10⁻⁵ 秒。
• 在 MOND 中，剧烈弛豫、相空间混合及星系并合的发生显著慢于牛顿引力，表明动力学过程受到抑制。
• MOND 系统中的动力摩擦强于等效的牛顿系统（含暗物质），表明在 MOND 中大质量天体的轨道衰减速率更高。