假如把太阳系里每颗行星、卫星和足够大的小行星都当成一枚棋子，按下“快进一亿年”，谁会撞上谁？现实中你没法重来，但在超级计算机里，天文学家可以让这场引力棋局从头播放。这里的关键工具，就是N体模拟。
把引力变成慢放电影
一句话定义：N体模拟，是把N个天体（这里的“体”可以是行星、恒星、暗物质粒子）看作有质量的点，让它们只按引力相互拉扯，再由计算机一小步一小步推演未来或回溯过去。
你可以把它想成一场没有裁判的台球赛：每颗球都同时被其他所有球轻轻拉一把。计算机每隔极短时间算一次所有拉力，再让每颗球挪动一点点。步子越小，轨迹越准，但算得也越慢。
两三个天体时，方程还算“听话”；到了成千上万颗恒星，或上亿个暗物质粒子，笔算就彻底失效。N体模拟不是望远镜观测，而是理论预言和数值实验：它告诉我们“如果引力规则成立，系统会变成什么样”。
为什么算起来会“爆炸”
麻烦来自平方增长。100个天体两两配对约有4950对；100万个天体则有约5000亿对相互作用。若每对只花一纳秒，算一步也要约500秒，而真实模拟往往要算几十亿步。到宇宙学尺度，直接硬算几乎不可能。
于是天文学家换算法：远处的一团天体可以被近似成“一个总质量”，近处才逐个细算；空间也被划成网格，用快速傅里叶方法处理大尺度引力。这些技巧不改变物理目标，只是把算力花在刀刃上。
2005年的“千年模拟”用约100亿个粒子追踪宇宙结构。把100亿个粒子分给当时地球上的每个人，每人能分到一个还多。它让人类第一次看清暗物质如何织成纤维状网络——这是理论计算给出的图景，不是直接拍到的照片。
屏幕里的星系像不像真的
在宇宙学里，很多模拟以N体模拟为骨架，再加入气体、恒星形成、超新星反馈等“配料”。如果模拟出的星系分布、旋转曲线和观测对得上，模型就获得支持；对不上，就得回头修改物理假设。
但要分清三层：观测证实的是星系在旋转、宇宙在膨胀；理论预言的是暗物质晕如何合并；假说推测的则是暗物质粒子到底是什么。这类模拟能复现大尺度结构，却不能单独证明暗物质是什么。小尺度上的一些差异，例如卫星星系数量、星系中心密度，目前仍有争论，尚无定论。
一亿年后的那颗小行星
回到开头：一亿年后谁撞上谁？对太阳系长期演化，这类计算能给出概率图景，但初始条件的一点点误差会被放大，所以答案不是宿命，而是一组可能。你看到的每一条轨道线，都是“在已知规则下的一种未来”。
更远的问题是：如果引力之外还有未知力量，或者暗物质并不像我们想的那样“冷”，屏幕里的宇宙会不会换一副面孔？下一次你看到星系模拟动画，不妨想一想：那不仅是图像，更是一场关于我们能否读懂引力的长期考试。
延伸阅读
- 书籍：James Binney & Scott Tremaine，《Galactic Dynamics》，Princeton University Press。
- 纪录片：《宇宙时空之旅》（Cosmos: A Spacetime Odyssey，2014）。
- 机构网站：Illustris TNG 项目官网（tng-project.org），可查看宇宙学模拟图像。