二十四旋涡星系旋转曲线的严格计算与暗物质问题.DOC

半径有限薄圆盘引力场的精确数值计算和银河系的旋转曲线与暗物质 梅晓春 沈世银 徐宽俞平（1）福州原创物理研究所（2）福州大学物理系，（3）上海天文台（4）清华大学核能与新能源研究院 （5）Cognitech Computer Inc, C A, U. S. A 内容摘要 有限薄圆盘均匀分布的引力场是引力理论的一个基本问题，然而由于数学上的困难，物理学家至今对这个问题仍缺乏的认识。本文严格按照牛顿引力理论，采用计算机数值计算方法，计算了半径有限物质薄盘均匀分布的引力场，证明它与物质球对称均匀分布引力场有很大的区别。在此基础上计算了有限盘密度指数衰减引力场的旋转曲线将它与无限密度指数衰减截断计算的旋转曲线进行比较结果显示存在大的差别有限盘的数值方法，可以漩涡星系的引力和旋转曲线，正确估计漩涡星系暗物质的量文用这种方法计算银河系的，结果暗物质。1．半径无限薄圆盘的引力场和旋转曲线 虽然牛顿引力理论已经建立三百多年，由于数学上的困难，我们至今对物质薄圆盘状分布引力场的特性实际上仍然缺乏很好的了解。为了能进行简单计算，目前天体力学中常用质量球对称分布代替盘状分布，近似地估计漩涡星系的旋转曲线。按牛顿引力理论，当质量分布具有球对称性，且物体在球外绕球心稳定旋转时，物体的运动速度与其到球心的距离之间满足以下关系： 或 （1） 对于质量的球状星系，假设物质均匀分布在半径的球内，按（1）式计算得到图1中所示的旋转曲线。在球半径以内曲线上升，在球半径以外曲线下降。然而对漩涡星系中星体光谱多普勒频移的大量观察结果发现，在相当大的范围内星体运动速度与其离星系中心的距离无关，大致呈现常数的关系。然而一般漩涡星系的旋转曲线大体上都是平坦的或向上弯曲的，这种现象是用（1）式无法解释的。 为了解释漩涡星系中星体运动速度与距离关系，物理学家提出两种方法。一是对牛顿引力理论进行修正，认为在宇宙尺度上牛顿引力理论应当做修改，这就是所谓的MOND理论。二是假设漩涡星系中存在大量暗物质。暗物质包括两种，一种是重子暗物质，如星际稀薄气体、不发光的行星和已经死亡冷却了的恒星等等。另一种是非重子暗物质，除了不参与电磁相互作用和不发光外，它们的具体特性不详。根据宇宙学理论的估计，宇宙中非重子暗物质大约是正常重子物质的6倍以上。重子暗物质的存在是一个基本事实，根据天文观察，星系中不发光重子暗物质的量发光重子物质的倍。而非重子暗物质至今仅是一种假设，天体物理学家努力了30多年，至今一直没有找到非重子暗物质。宇宙中是否存在如此多的非重子暗物质已经成为现代天体物理学和宇宙学的最大谜团之一。 现有天体力学上更为精确的方法是，常将漩涡星系看成一个边界无限的薄圆盘。先设法求出盘引力势的表达式，然后再通过微分计算求引力。令代表薄圆盘面上某点的坐标，代表空间任意观察点的坐标，一个半径无穷大的薄圆盘在空间任意点产生的引力势为： （2） 计算上式时涉及到椭圆函数的积分，目前按这种方法只得到少数几种质量分布的引力势。一种简单的形式是，假设时薄圆盘的质量面密度为： （3） 其中是常数面密度，是具有长度量纲的常数。通过椭圆积分的计算，得到旋转曲线公式为： （4） 按上式旋转速度是一个常数，然而问题在于实际漩涡星系的质量密度分布不满足（3）式。同时由于时，（3）式也是不合适的。按（3）式薄圆盘的总质量为： （5） 这也是不合适的。如果假设薄盘质量面密度分布具有指数衰减的形式，即令： （6） 式中参数，是盘的长度标度。在时，按椭圆函数计算方法可得： （7） 式中，和为第一类修正的贝赛尔函数，和为第二类修正的贝赛尔函数。 注意到按（6）式的物质分布形式，圆盘的半径是无限大的。然而漩涡星系的半径都是有限的，对于边界有限的薄圆盘，目前没有找到解析解。为了能利用（7）式计算漩涡星系的旋转速度，目前一般采用半径无限大的薄圆盘加上人为截断的方法，认为圆盘半径以外的区域物质密度很小，可以忽略不计。因此我们只需考虑半径范围内薄圆盘的质量，就有： （8） 时为有限值。选择适当的，得到图1中曲线，其中。曲线表示如果是质量球对称分布，且在半径为的球内质量也是时的旋转曲线。则是物质全部集中在球心（坐标系原点）时的开普勒旋转曲线。可见旋转曲线与质量球对称分布的旋转曲