---
title: "怎样读反馈、暗晕与MOND的条件比较"
description: "逐图阅读固定2026年Smith等人论文，分清核—尖峰、重子反馈、几何先验和弱透镜平均曲线，并完成一个势场变化的轨道算例。"
canonical_url: "https://sandboxphysics.com/zh/stories/galaxy-rotation-curves/reading-feedback-and-gravity-comparisons/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/galaxy-rotation-curves/reading-feedback-and-gravity-comparisons/"
last_modified: "2026-11-05"
---

[星系旋转曲线与质量分解](https://sandboxphysics.com/zh/stories/galaxy-rotation-curves/)

星系旋转曲线 · L4

# 怎样读反馈、暗晕与MOND的条件比较

Sandbox Physics · 2026-11-05

比较暗晕与修改动力学，先要看用了哪一组观测、哪些参数被固定，以及误差怎样处理。Smith等人的一份2026年论文把这些选择摆在同一研究里，适合沿着图表逐项读。

## 开始前需要什么

读完前三章及SPARC方法教程，了解有符号引力分解、共同几何误差和工作似然。旁边打开arXiv:2610.33157v1的固定PDF。

## 先记版本与147个对象怎样选出来

主源是Amy Smith等人的“Dark Matter and MOND successes and degeneracies in reproducing SPARC rotation curves and weak lensing data”，固定arXiv:2610.33157v1。水印日期为2026年10月23日，摘要页注MNRAS已接收；以下按实际核读的26页预印版本展开。PDF的2025排版年不作为发表年份。

第2节从SPARC选出 $147$ 个晚型星系，排除质量标记为3、低倾角、数据过少等情况，并说明了总速度平方无效时的点选择。读统计图前先看这些条件；它们定义了作者比较的是哪批星系。

本课程的固定NGC2403分支保留自己的73个样点，不复制这份多星系研究的删选。一个对象的NFW与平核拟合，和一整批星系的模型比较，回答的尺度也不同。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Lelli, McGaugh & Schombert (2016), SPARC, AJ152,157

## 图1画的是密度，反馈来自另一套模拟经验

图1横轴是千秒差距半径，纵轴是暗晕密度，都是对数坐标。用图例中的模型名辨认曲线，再看图注给定的示例晕质量与光度。密度有平核、尖峰或过渡，不能直接读成同样形状的速度曲线；速度还要先积分质量。

第3.1节的NFW曲线来自纯无碰撞暗物质引力模拟的经验形状；DC14则用包含重子反馈的流体模拟结果来规定密度轮廓。反馈是恒星形成、爆炸和气体运动改变引力势的过程。暗物质可以在势场变化中改变轨道，并不需要被气体风直接撞开。

DC14把轮廓形状与恒星—晕质量比联系起来。图1展示模型预测；为它提供依据的模拟计算了密度轮廓怎样变化。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

## 图5先看内侧残差，再回到表1与表3

图5给出UGC01281：上排采用NFW，下排采用DC14，每排又有总曲线与残差。比较时先找同一半径处的偏离，再看恒星、气体与晕怎样分担引力。作者用它展示DC14怎样缓解这个对象的内侧失配。

表3把比较扩到样本中位数。固定距离与倾角时，NFW工作约化卡方的中位数是 $1.387$，DC14为 $0.934$；允许几何变化的对应列是 $1.485$ 与 $1.157$。列里的参数数目也从3变为5。不能把中位数的变化只解释成一条曲线“变好多少”：拟合条件和自由度都在变。

再翻表1：质量光度比、尺度与几何都有指定的先验或范围。先验表示拟合前愿意给哪些参数值多少权重。对数正态宽度的dex单位与线性标准差不同。重做数值比较时，要把分布形式、宽度单位和参数范围完整写进拟合程序。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

## 图6的四格没有共同的零点模型

第4.1.6节采用 $\mathrm{BIC}=-2\ln L+p\ln N$，其中 $p$ 是自由参数数目，$N$ 是一个星系的测量点数。BIC用拟合改善与参数代价作近似权衡；它有大样本与模型规则等适用条件，也不是暗物质粒子概率。

图6横轴是 $\Delta\mathrm{BIC}=\mathrm{BIC}_X-\mathrm{BIC}_{\rm ref}$，纵轴是累积分布：在每个横轴值处，显示比较值不超过它的星系比例。正值表示参考模型的BIC更小。上排参考是DC14；下排参考改成相应的固定加速度尺度MOND插值模型。把上下两排的曲线位置直接相减，便混用了参考模型。这里按第4.1.6节正文，把两种插值函数各自对照其固定加速度尺度版本。图6右下图注只提Simple MOND，与正文的表述不完全一致；尚未取得实现，便保留这处文字差异，不另作裁决。

MOND是修改牛顿动力学的缩写，尝试改变低加速度下的运动关系。这里采用作者声明的修改惯性实现，比较两种连接高、低加速度行为的插值函数，以及外部环境引力影响内部运动的一维近似。它与重新求解所有修改引力场方程并不是同一件事。作者还把每星系自由加速度尺度作为退化演示，不能把它默认为一个普适常数理论本来就有的自由度。

小算例：对同一组 $73$ 点和同一误差坐标，多一个参数若只让 $\chi^2$ 降低 $2$，BIC变化为 $-2+\ln73\simeq2.29$，旧模型反而更低。这个代数例子说明分数怎样工作；真实跨理论比较还需共同的数据定义、合理先验和模型适用性检查。参数的物理名称不同，本身并不是禁止统计比较的数学定理。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Independent saved-observation rotation-curve calculation

## 图10的百万秒差距曲线来自叠加透镜信号

到图10，横轴单位变成 $\mathrm{Mpc}$。近处蓝线是按重子质量分组后平均的HI旋转曲线，远处点则由弱引力透镜的平均表面质量信号去投影成加速度，再表示成等效圆周速度。远处没有一圈被直接追踪到的气体测试粒子。

第4.3节与图注说明了分组、至少三个对象才取平均的HI径向箱、剔除的内侧两个透镜点，以及一百万秒差距之外的孤立性限制。比较曲线还使用不同的晕质量—集中度条件。每个选择都影响这次特定比较。

读图时把总拟合与外侧拟合分开。表6中，有些暗晕模型的全范围约化卡方不大，外侧分数却明显较高；较多的内侧数据可以掩盖外侧失配。图14再检查这些拟合参数是否符合所选的晕标度关系。

作者发现其所试模型中，MOND能较自然地维持大尺度平坦曲线，并讨论束缚晕之外物质等进一步解释。这里完整读过固定版本，但未取得其私下提供的完整弱透镜协方差，也没有重跑147星系MCMC或透镜去投影。这个结果应连同上述条件一起阅读。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

## 例题：气体离开后，无碰撞粒子的轨道怎样改变

用一个小模型理解反馈为什么可能影响暗物质：让测试粒子起初绕中心质量作圆周运动，然后中心质量瞬间减少十分之一。瞬间速度来不及变，角动量也保持原值；引力势却变浅，粒子开始沿更伸长的轨道运动。

初始质量为 $M$、圆轨半径为 $r_0$，有 $v_0^2=GM/r_0$。移走比例 $f$ 后，单位质量能量成为 $E'=(f-\tfrac12)GM/r_0$。单位质量角动量为 $h=r_0v_0$。到达远心点时，径向速度为零，因此 $E'=h^2/(2r_a^2)-G(1-f)M/r_a$。代入上式可得两个转向点：原来的 $r_0$，以及 $0<f<1/2$ 时的远心半径 $r_a=r_0/(1-2f)$。取 $f=0.1$，得到 $r_a=1.25r_0$。

独立材料用转向点能量方程的数值求根，核对三个移走比例下的解析结果。它展示时变引力势怎样改变无碰撞粒子的轨道。真实反馈还需要气体何时流走、怎样返回、重复多少次以及自引力等条件；这个点质量例题没有生成DC14轮廓。

对应阅读：Independent saved-observation rotation-curve calculation · Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1

## 把比较条件写在结论旁边

读完后，先用自己的话说清：比较的是个体旋转还是叠加透镜；核与尖峰怎样由模型进入；距离与倾角是否固定；每格分数相对于谁；外侧物质与协方差怎样处理。这样读出的结论，才保留了原研究真正回答的问题。

下一次观测能做什么，也由这些细节决定。更好的空间分辨率能检查内侧曲线，独立几何与恒星质量资料能缩小退化，透镜则把问题带到更远尺度。几种证据一起推进，才有机会区分相近曲线背后的物理。

对应阅读：Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1 · Independent saved-observation rotation-curve calculation

## 练习与解释性答案

1. 图6某星系对NFW的上排比较值为 $\Delta\mathrm{BIC}=+8$。这个正号支持谁？能否由此读出“MOND概率”？

   展开解释

   上排参考是DC14，正号表示NFW的BIC更大，在这次BIC比较里DC14较受支持。下排另有MOND内部的参考模型，不能把这一个数换成跨理论概率。

2. 为什么图10远处的“速度”需要与个体星系气体速度分开说明？

   展开解释

   远处点来自叠加弱透镜信号、几何去投影和动力学换算，涉及不同星系群体、选择与协方差。它不是在同一个星系百万秒差距处直接量出的气体轨道速度。

3. 点质量例题中若瞬间移走 $f=0.2$，远心半径是多少？为何不能由此宣称已经算出一张真实星系的平核？

   展开解释

   $r_a/r_0=1/(1-0.4)=5/3$。这只计算了一个测试粒子在规定中心势变化后的轨道。真实密度轮廓还需大量轨道、自引力、气体历程与重复反馈共同演化。

## 动手检验预测

### [比较恒星盘与条件暗晕](https://sandboxphysics.com/zh/labs/galaxy-rotation-curves/)

1. 分别查看尖峰与平核的手动正演，用内侧分量曲线描述差别。
2. 展开几何设置，观察同一份数据的残差如何一起变化。
3. 把Lab作为牛顿质量模型的有限示例，再回到原论文图6与图10，列出那里新增的模型、观测与统计条件。

## 独立科学复算材料

[下载科学复算ZIP](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/galaxy-rotation-science-v1.zip) · [阅读安装与计算范围](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/README.md)

离线计算需要Python、NumPy与SciPy，默认合成数据会完整生成并重新拟合。可选的具名SPARC文件只在自己的电脑取得；ZIP不含原始观测或浏览器交互引擎。

单件脚本、字节与SHA-256

- [README.md](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/README.md) (6785 B)
  `5fbf0d29e006fc06f35d2290fd8fee7d5abc2301b601814165495331f9c37a84`
- [LICENSE.txt](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/LICENSE.txt) (1325 B)
  `4f15facc82fd0218a74907ae83f1716569277d9c40f55b249451dd8d89191b2f`
- [requirements.txt](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/requirements.txt) (185 B)
  `359f4c650ca6dd5a2705535e32377aa86db27d8d6ebd9b47333070850b9c5980`
- [source-identities.json](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/source-identities.json) (724 B)
  `409cf4931c19dc97049f943eb005165a59f88d1842a3dd2596eaa1dabcfbd458`
- [own-observation.json](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/own-observation.json) (9009 B)
  `48cd32150a8eff02b4c3e870cb82f8df008d90c2b5aed3f9bca743506852a5e1`
- [analyze.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/analyze.py) (2750 B)
  `5baed281117d3980af0f0e19e428c4163f79a97a84d1db71955ecbee7993ac3a`
- [conditional\_fit.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/conditional_fit.py) (3618 B)
  `37700d28edc852a10d09f949a7270a652b61f1c9d30cc2d8d0e99d4347e92d70`
- [fetch\_named\_source.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/fetch_named_source.py) (5014 B)
  `240b29a1c6e929895b5ead31269b495b1223be0615da62c703e2962028e6de5a`
- [generate\_observation.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/generate_observation.py) (849 B)
  `c5442424f9cecdcd81022ec31cae0b81bdc21400f06eea6dde4108c7062de743`
- [observation\_record.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/observation_record.py) (2693 B)
  `2d791013b5c8b289e6bfdc5ba9ed408e18b9fc4334d44ba32091546ed297e752`
- [physical\_units.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/physical_units.py) (255 B)
  `2db6448e5f4e657ef48ffd9d2a58c7487c6f7655a8033dc001327fe43192507c`
- [plot.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/plot.py) (5734 B)
  `c30e8e160f1f73f6451ebb903402a316ecba858814e6b1d5059ad31926abb4e1`
- [rotation\_model.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/rotation_model.py) (2003 B)
  `97c2c5ca56373f20996b26ef179267549ad3a341817b29666ce368029cfa96a1`
- [save\_native.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/save_native.py) (655 B)
  `bbf0881e28cf1dd4d13105af656cc4072e3a42a62a5f7013b9f48075b764c082`
- [spherical\_halos.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/spherical_halos.py) (1459 B)
  `4678e06bf172167f3a8a3a5d199447658eb51860507405b80110807165c5ef83`
- [synthetic\_record.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/synthetic_record.py) (1538 B)
  `6285cda07f2b9324806d25ecc89712eb6b94479e3caaa1f24a88c004d449a063`
- [thin\_disk\_reference.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/thin_disk_reference.py) (1298 B)
  `67ba0512e47c10d10c5b8d17f435d808812eec25f0edce296b6c80eaf72e0bc0`
- [validate.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/validate.py) (5709 B)
  `7cfbf116c5e51e9a3e1aa0ecd113a24f7a069706336ae01221de67dbf331cb6d`
- [verify\_analytic\_limits.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/verify_analytic_limits.py) (3073 B)
  `68bca45f2d3d8a1f8cda804bb45c3e165025effa0221c2f45d7c1416009bae6a`
- [verify\_nuisance\_mapping.py](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/verify_nuisance_mapping.py) (4184 B)
  `04b8587eb21a257b809fe930c33d43ddb2e008e4ebf2026134ee0a1d35b4015e`
- [SHA256SUMS.json](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/SHA256SUMS.json) (3306 B)
  `2fb52c937f2026d872f053e17343e55595fb3568862c20573d59123f3d897b49`
- [galaxy-rotation-science-v1.zip](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/galaxy-rotation-science-v1.zip) (29081 B)
  `f389268ed651be3a66ad7d87ee6760ebe39d87f6bd2498b2fb8904c9f10893a3`

[上一篇：从具名SPARC文件重新算质量模型](https://sandboxphysics.com/zh/stories/galaxy-rotation-curves/recomputing-a-named-sparc-galaxy/)

## 来源与计算范围

1. [Lelli, McGaugh & Schombert (2016), SPARC, AJ152,157](https://astroweb.cwru.edu/SPARC/)

   官方固定Table1与Rotmod\_LTG.zip，只使用具名NGC2403成员。主论文第3.1至3.3及5.1至5.3节、原Table2说明用于核读气体符号、单位恒星质量光度比、有限厚度与误差语义。新取源与缓存复算分别记录；原观测与实测派生图留在本地。

2. [Smith et al., Dark Matter and MOND successes and degeneracies, arXiv:2610.33157v1](https://arxiv.org/abs/2610.33157v1)

   完整核读固定26页预印本，并视觉检查图1/5/6/10/14与表3。水印2026年10月23日；摘要页注MNRAS已接收，PDF仍留接收占位符和2025排版年。指南沿第2至6节，区分固定与可变几何、BIC各自参考与叠加透镜条件。未重现147对象MCMC或私下提供的协方差，不再发布原PDF或原图。

3. [Independent saved-observation rotation-curve calculation](https://sandboxphysics.com/stories/galaxy-rotation-curves/materials/galaxy-rotation-science-v1.zip)

   自有MIT Python实际生成合成观测、六起点拟合、轮廓、共同几何场景、正密度/圆环势/Jacobian参照与SVG。可选具名源获取、73点保存转换和拟合均实际执行。默认离线，ZIP不含Lab引擎、原SPARC行、论文或实测派生图。
