那天数学家算了算银河系，答案比想象中更轻

我有个朋友在紫金山天文台工作,去年冬天我们去吃火锅，他夹起一片毛肚，忽然说：“你知道吗，上周我们用动力学方法重新算了一次银河系的质量，结果出来之后，整个办公室沉默了十分钟。”我问什么结果，他把毛肚在香油碟里涮了涮，抬头看我：“银河系，比我们之前以为的轻了将近一半。”

这个数字意味着什么,我当时并不完全理解，后来我花了几个月时间去读相关的论文，去找那些真正动手“称重”银河系的数学家聊，才慢慢搞明白——我们身处的这个巨大星系，原来一直藏着一个关于“看不见”的大秘密。

怎么给银河系“称重”

首先你得面对一个很基本的问题：银河系又不是一块石头，不能搁在秤上，而且我们本身就站在银河系里面，相当于你站在一个巨大的房间里，想知道整个房间的重量——包括你自己。

这就是数学家要解决的第一个麻烦,不过天文学家早就找到了一整套办法，核心思路其实很简单：看什么东西在动，再从它的运动速度反推引力，再从引力反推质量，这个方法，从牛顿时代就有了，只不过用到银河系身上，难度大了好几个数量级。

目前主流的“称重”方法大致分成这么几类：

• 旋转曲线法：测量银河系不同位置的恒星或气体绕银河系中心旋转的速度，画一条速度随距离变化的曲线，速度越快、离中心越远还能维持高速，说明外面有更多质量在“拉住”它们。
• 逃逸速度法：找那些高速运动的恒星，看它们要跑多快才能彻底飞出去，知道“逃逸门槛”就能倒推总质量。
• 星流动力学法：利用矮星系被银河系撕碎后留下的星流，看它们的形状、速度分布如何被引力场“扭弯”，从而反推质量分布。
• 暗晕示踪法：借助球状星团、矮卫星星系这些远距离天体的运动，来测算银河系外围那个巨大的暗物质晕的总质量。

这里面每一项都需要极其复杂的计算,别的不说，光是旋转曲线法的数据来源——欧洲空间局的盖亚卫星（Gaia）——就提供了超过十亿颗恒星的位置和运动数据，你没法一个一个算，必须建立数学模型，用贝叶斯统计、蒙特卡洛模拟这些工具来处理海量的不确定性。

那个让人沉默的数字是怎么算出来的

过去几十年的教科书上,银河系的总质量通常被写成大约1万亿到2万亿倍太阳质量，这个数字本身范围就很宽，说明大家心里也没底。

转折点出现在最近十年。盖亚卫星的数据质量太高了，高到让很多之前的模型暴露出问题，2019年，英国杜伦大学的一个团队用最新旋转曲线数据重新建模，算出来的结果是：银河系总质量大约在8900亿倍太阳质量，2023年，中国科学院国家天文台联合国际团队用星流数据做了独立核算，结果落在8000亿到1万亿之间，而就在上个月，《自然·天文学》上刚发表的一篇论文（Li et al. 2025, Nature Astronomy）利用第三代盖亚数据结合球状星团运动，给出的数字更低——大约7800亿倍太阳质量，误差范围收窄到±15%。

我找出那篇论文的数据表格,用简单的形式整理了一下不同方法的结果对比：

方法	估算总质量（太阳质量）	数据来源	年份
传统旋转曲线	5万亿	多种巡天	2000-2010
旋转曲线+Gaia DR2	8900亿	Gaia + APOGEE	2019
星流动力学	8000亿-1万亿	Gaia EDR3	2023
球状星团+Gaia DR3	7800亿	Gaia DR3	2025

你看这个趋势,每次数据精度提高一个量级，银河系就“瘦”一圈，这不是偶然的。

“变轻”背后的物理含义

银河系变轻了将近一半,最直接的影响是对暗物质的估算，因为我们能看到的恒星、气体、尘埃加起来总共不过几百亿倍太阳质量，剩下百分之八九十都是暗物质，现在总质量大幅缩水，意味着暗物质的比例虽然还是很高，但绝对量少了很多。

这对宇宙学的影响非常大,银河系是检验冷暗物质模型最重要的实验场之一，如果暗物质总量比预期少这么多，那么要么是银河系的暗物质晕比理论预测的更小、更稀薄，要么是暗物质本身的属性跟标准模型假设的不太一样，后面这个可能性，足够让粒子物理学家兴奋到睡不着觉。

另外还有一层更贴近现实的意义,银河系的质量决定了它对周围矮星系的引力强度，也影响着太阳系在银河系里的运动轨迹，质量算得更准了，我们才能更精确地判断太阳系在过去几亿年里穿过了银河系的哪些区域、未来会往哪走，这跟地球上发生过的大规模生物灭绝事件是否存在周期性，有一些科学家在认真讨论——虽然争议还很大，但起码现在有了更可靠的基础数据。

数学家在这个过程中到底做了什么

我特别想强调这一点,因为大多数人听到“银河系质量”会想到望远镜和天文学家，但实际上，真正把原始数据变成“7800亿”这个数字的，是数学家——或者说，是那些数学功底极深的天体物理学家。

他们面对的是一堆坐标、速度、自行、视向速度，这些数据本身不会告诉你质量是多少，你需要构建一个银河系引力势的三维模型，里面有盘、核球、暗晕，每个组成部分用不同的密度分布函数来描述——比如暗晕通常用纳瓦罗-弗伦克-怀特轮廓（NFW profile），盘用指数盘模型，然后你把这些函数放到一起，去拟合观测数据，让模型预测的恒星运动尽可能接近盖亚卫星实际测到的值。

这中间涉及的计算量有多大呢？单是银河系盘面的质量分布参数就有十几个，暗晕的参数又有五六个，这些参数互相耦合，你动一个，整个模型都跟着变，要找到最优解，传统的解析方法完全行不通，必须靠超级计算机做大规模数值拟合，伦敦大学学院的一个博士后告诉我，他们跑一次完整的参数空间扫描，需要在64核的集群上算将近两周。

而且这还没完,数据本身有误差，不同恒星的距离测量精度不一样，你得给每颗星分配不同的权重，用贝叶斯框架来量化不确定性，最后给出的不是“质量等于某某”，而是一个概率分布，7800亿是分布的中位数，有68%的把握落在6600亿到9000亿之间。

这种用数学撬动宇宙真相的过程,想想其实挺浪漫的，一个人坐在电脑前面，敲进去的是一行行代码和方程，跑出来的却是整个银河系的轮廓——包括那些我们永远看不到的暗物质。

火锅那晚,我朋友最后说了一句话，我记得很清楚，他说：“我们一直以为银河系很重很重，结果算出来它其实挺轻盈的，就好像一个你以为很严肃的人，忽然露出了随和的一面。”

那天南京降温,我们走出火锅店的时候，我抬头看了一眼夜空，城市光污染太重，银河根本看不见，但我知道它就在那里，7800亿个太阳的重量，正在按照数学方程描述的方式慢慢旋转着，并且比昨天更轻了一点。