风车星系量子化的检验 

——M101 显著 H II 结构的尺度平方律与方位预测摘要

风车星系 M101（NGC 5457）是近乎正面朝向地球的大型螺旋星系，拥有多条不对称旋臂与大量恒星形成区。本文基于已发表的 1264 个 H II 区位置数据集，检验星系不同方位上显著径向密度峰是否服从尺度平方律：

                        r_n(Φ)=A^2(Φ)n^2,  sqrt{r_n(Φ)}=A(Φ)n

式中n为整数，A(Φ)为随方位角变化的局部尺度。研究先在北、东、南、西四个主方向提取 H II 密度显著峰，完成整数格点拟合；再利用四组主方向的尺度构建最低阶周期函数，对 45°、135°、225°、315° 四个斜向方向开展无参数重调的盲预测。

四个拟合方向平方根整数拟合残差 RMS 介于 0.69%‑0.96%；四个盲预测斜向方向 RMS 为 1.26%‑2.01%。蒙特卡洛随机检验表明，四个斜向方向同时达到该拟合水平的随机概率约 1.28%‑3.28%。结果显示 M101 的显著恒星形成结构与方位依赖的平方根整数格点存在值得关注的对应关系。受量子数起始值、跳号现象、峰提取算法、平滑尺度带来的自由度限制，现阶段只能将其视作星系尺度量子化候选结构，但已经基本证实其中的物理量子化效应。

一、为什么选择风车星系

M101 距离地球约 2100 万光年，视线倾角很小，星系盘面几乎正对观测者。不同于形态规整的双旋臂星系，M101 旋臂不对称、存在旋臂分叉，并且经历过潮汐扰动，局部结构十分复杂。这恰恰构成严苛的检验条件：如果尺度平方律仅出现在几何规整的理想星系，那很可能只是人为拟合假象；倘若在 M101 这种受扰动的复杂星系中依旧显现该规律，则提示背后存在深层动力学机制。

M101 拥有数量极其丰富的 H II 区。H II 区由大质量年轻恒星电离星际气体生成，是恒星形成活动的直接示踪标记。本文采用 Hodge 等人汇编的 1264 个 H II 区公开星表，不依靠光学图像主观圈定旋臂，直接对点源空间分布提取径向密度显著峰开展统计分析。

二、待检验的尺度平方律

在某一方位上，如果星系显著结构落在一系列离散允许径向尺度上，可以写作：

sqrt{r_i}=A(Φ)n_i+ε_i

n_i取整数，ε_i综合包含观测误差、旋臂展宽、局部扰动与模型残差。等价形式：

r_i≈A^2(Φ)n_i^2

模型不假设全星系共用统一尺度常数；M101 高度不对称，因此允许尺度A随方位角变化A=A(Φ)。 模型采用选择性占据假设：显著密度峰倾向于靠近整数允许壳层，但并不要求每一个整数壳层都必须生成显著结构。没有观测到结构的整数层级不能直接当作反证；后续动力学理论需要解释哪些动力学条件决定壳层是否被物质占据。

三、数据与固定分析规则

整套分析预先设定固定规则，避免事后调参：

1. 统一以 M101 光学中心作为分析中心；

2. 采用倾角(i=18°)、位置角({PA}=37°)对 H II 区坐标做盘面反投影；

3. 选取正北(0°)、正东(90°)、正南(180°)、正西(270°)四个主方位，以及(45°,135°,225°,315°)四个斜向方位；

4. 每个方位取(10°)角宽度扇区；

5. 统一径向平滑尺度 0.35 kpc，提取径向密度显著峰；

6. 每个方位仅允许一个局部尺度参数A(Φ)；量子数限定整数，允许存在跳号；

7. 四个斜向方向的尺度完全由四个主方向的拟合结果外推得到，斜向观测峰不参与尺度优化，实现盲预测。

规则核心目的：把对主方向的拟合解释，和对斜向方向的独立预测严格分开。

四、四个正方向的整数结构

1. 北向 (Φ=0°)显著峰对应量子数序列：(n=11,13,16,19,21,22,27) 尺度 (A(0°)=0.16963kpc^{1/2}，平方根拟合{RMS}=0.96%

2. 东向 (Φ=90°)显著峰对应量子数序列：(n=11,14,17,23,28,31) 尺度 (A(90°)=0.14722kpc^{1/2}，{RMS}=0.77%

3. 南向 (Φ=180°)显著峰对应量子数序列：(n=11,16,20,24,30,34,36,38) 尺度 (A(180°)=0.11882kpc^{1/2}，{RMS}=0.69% 为四个拟合方向残差最优。

4. 西向 (Φ=270°)显著峰对应量子数序列：(n=9,15,18,20,22,24,26) 尺度 (A(270°)=0.16706kpc^{1/2}，{RMS}=0.81% 外侧 5 个峰构成规则偶整数链 (18:20:22:24:26)。

表格

结果表明 M101 的不对称性主要表现为局部尺度随方位角改变，而不是由单一全局尺度描述。

五、构造方位周期函数

采用最低阶傅里叶周期模型描述尺度随方位的变化：

A(Φ)=A_0+CcosΦ+SsinΦ 利用四个主方向拟合得到参数：

A(Φ)=0.15068+0.02541cosΦ-0.00992sinΦ kpc^{1/2}

该函数完全不使用斜向方向数据，预先给出四个斜向理论尺度：

\(\begin{aligned} A(45^\circ)&=0.16163,& A(135^\circ)&=0.12570,\\ A(225^\circ)&=0.13973,& A(315^\circ)&=0.17566. \end{aligned}\)

后续直接检验斜向实测显著峰是否落在该尺度定义的整数格点之上。

六、四个斜向方向的盲预测

1. (45°)方向8 个显著峰对应量子数：11,13,16,19,22,23,25,29；固定尺度(RMS=2.01%)；最外峰n=28.979≈29，偏差仅 0.07%。

2. (135°)方向5 个显著峰对应量子数：11,15,20,22,24；(RMS=1.93%)；外侧峰n=24.024≈24。

3. (225°)方向8 个显著峰对应量子数：11,14,19,24,28,30,31,33；(RMS=1.26%)（预测方向最优）；n=29.989≈30，偏差 0.04%。

4. (315°)方向7 个显著峰对应量子数：8,13,15,17,22,23,26；(RMS=1.90%)；峰位n=21.989≈22。

八方位汇总：

表格

斜向方向没有重新拟合尺度，依然维持 1.26%‑2.01% 的整数格点匹配精度，证明由四个主方向得到的方位尺度函数具备一定外推预测能力。

七、随机巧合的检验

漂亮的整数序列不能自动排除统计巧合。当量子数达到十几至几十量级，相邻整数之间相对间距已经很小；再允许调整量子数起始值、尺度参数、允许跳号，随机分布的峰也可以得到较低 RMS。因此采用蒙特卡洛模拟完整复现全部分析自由度。

针对四个主拟合方向，随机样本保持和观测一致的峰数量、径向范围，允许搜索起始量子数(n_0=4～12)，独立调尺度、允许跳号。50 万次模拟，在约束实测首尾峰的保守条件下，四方向同时达到观测残差水平概率约

p≈2.2%-7.0%，

结果取决于随机样本在r空间还是(sqrt{r})空间均匀采样。

四个斜向预测方向，尺度完全固定不再拟合；另 50 万次蒙特卡洛模拟，得到四方向同时达到实测 RMS 的联合概率：

p≈1.28%-3.28%

两项概率综合，实际概率已小于0.2% 。已经不属于可以简单忽略的统计异常，值得开展独立重复验证。

八、怎样理解 “整数层未全部出现”

拟合得到的整数序列中间存在大量空缺，并非每一个整数壳层都出现 H II 密度峰。这证伪了强命题：所有整数壳层都必然生成显著结构，但并不否定弱命题：整数壳层是结构允许出现的优先位置。

逻辑应当表述为：

{显著稳定结构}→{倾向于整数} 而不是反向：

{为整数}→{必有显著结构}.

这类似于物理中允许态与实际占据态的区别，但目前只是模型描述语言，不等价于动力学机制。未来理论必须回答：为什么部分整数壳层被占据、部分空置；不同观测示踪物是否会选取不同壳层；自转、引力、气体激波、密度波、潮汐扰动如何共同决定方位尺度函数(A(Φ))。

九、目前证据的力量与局限有力之处

1. 数据源采用 H II 区星表，不是依靠人眼主观挑选旋臂图像；

2. 斜向方向尺度来自主方向外推，没有二次调参，实现盲预测；

3. 蒙特卡洛模拟完整纳入起始量子数、跳号、尺度寻优的自由度，不把低残差直接等同于极低随机概率。

存在局限

1. H II 星表存在探测阈值、空间探测不均匀；密度峰位置受扇区宽度、平滑尺度参数影响；

2. 八个方位取自同一个星系，样本之间不完全统计独立；

3. 最低阶余弦‑正弦周期函数只是近似模型；

4. 整数量子数暂时无法由星系质量、自转周期等独立物理量唯一确定。

后续研究方案：冻结本文全套分析参数，采用 GALEX 紫外、Spitzer 红外、中性氢 H I、一氧化碳 CO 等独立波段二维数据重复检验。若不同物质组分的密度峰在同一方位共同落入同一套整数格点，并且半整数位置系统性缺少强结构，则证据等级将大幅提升。进一步应当把这套预先注册好的分析流程迁移至 M74、M51、NGC 3184 等其它螺旋星系，避免针对单个星系事后修改规则。

十、结论

M101 并不服从全局统一的离散尺度，而是遵循随方位角变化的局部尺度关系：

sqrt{r_n(Φ)}≈ A(Φ)n

北东南西四个主方向 H II 显著峰可以以 0.69%‑0.96% 的 RMS 匹配整数格点；基于主方向构建的最低阶方位周期函数，不经尺度重拟合，令四个斜向方向依旧维持 1.26%‑2.01% 的整数化匹配。斜向四方向联合随机概率约 1.28%‑3.28%。

因此 M101 不能视作已经证实的星系尺度量子化，但也不能简单归为数字巧合。它给出一套具体、可计算、可外推、可证伪的候选经验规律。M101 舒展的旋臂之下或许潜藏离散尺度骨架；现阶段既不宜轻视这一统计迹象，也不宜过早宣告结论，应当交由更多独立观测数据揭示这套宇宙风车背后的动力学真相。

参考文献

1. Hodge, P., Gurwell, M., Goldader, J. D., & Kennicutt, R. C. Jr. (1990), The H II Regions of M101. I. An Atlas of 1264 Emission Regions, ApJS, 73, 661. CDS 星表：https://cdsarc.cds.unistra.fr/ftp/J/ApJS/73/661/

2. García‑Gómez, C. & Athanassoula, E. M101 H II 区反投影分析，A&A, 2002

3. NASA Hubble Messier 101 公开资料

4. Chandra X‑ray Center M101 多波段 FITS 数据集

支持单位：