蜜蜂理论 – 基础 – 技术说明 X

剖析残差:
磁盘大小的线性趋势

注释 IX 中的 94 个星系盲测表明,随着星盘大小的增大,会出现系统性的残差趋势。本注释定量地描述了这一趋势,分离出每一侧的最大偏差,并确定了分散的结构来源。

1.第一项结果

一个线性残差,两个相反的种群

预测误差与星盘尺度长度成线性比例:$text{error}\, (\%) \approx -31.7 + 12.8\,R_d$, 皮尔逊相关性为 $r = +0.75$。回归线在 $R_d = 2.48$ kpc 时与零点相交,这基本上是校准所依据的银河系圆盘大小。回归线的两个极端对应于两个物理上截然不同的异常群:一端是大质量螺旋体(预测过高),另一端是紧凑矮星(预测过低)。

2.残差与 $R_d$ 呈线性关系

将预测误差与 $R_d$ 的关系绘制成图,每个点都按哈勃类型着色,这样就能立刻看出趋势的线性关系。红线是所有 94 个盲星系的误差与 $R_d$ 的线性回归。

预测误差与磁盘大小的关系 – 与 $R_d$ 成线性关系,按哈勃类型着色 0.30.5123510 -50%-25%+0%+25%+50%+75% → 预测过高区域 ±30% 带 ← 预测不足区域 盘尺度长度 $R_d$ (kpc) – 对数尺度 预测误差 (%) F579-V1(Sd-Im) Rd=3.20 Vf=105 err=-12F583-1(Sd-Im) Rd=1.80 Vf=83 err=-30F583-4(Sd-Im) Rd=1.40 Vf=67 错误=-29IC2574(Sd-Im) Rd=2.80 Vf=69 误差=+35%。KK98-251(Sd-Im) Rd=0.30 Vf=17 err=-23M33(Sc-Scd) Rd=1.40 Vf=100 err=-2NGC0055(Sd-Im) Rd=1.80 Vf=87 err=-1NGC0100(Sc-Scd) Rd=2.30 Vf=83 err=+13NGC0247(Sd-Im) Rd=2.40 Vf=90 err=+20NGC0289(Sb-Sbc) Rd=3.50 Vf=155 err=+29NGC0300(Sd-Im) Rd=1.50 Vf=76 err=+0NGC0801(Sc-Scd) Rd=5.80 Vf=208 err=+57NGC0891(Sb-Sbc) Rd=4.10 Vf=212 err=+7NGC0925(Sd-Im) Rd=3.10 Vf=105 err=+48NGC1003(Sc-Scd) Rd=2.80 Vf=115 err=+12NGC1090(Sb-Sbc) Rd=3.80 Vf=170 err=+17NGC1705(S0-Sa) Rd=0.60 Vf=54 err=-19NGC2366(Sd-Im) Rd=1.30 Vf=55 err=+14NGC2403(Sc-Scd) Rd=1.80 Vf=131 err=-4NGC2683(Sb-Sbc) Rd=2.90 Vf=175 err=+15NGC2903(Sb-Sbc) Rd=2.60 Vf=184 err=-0NGC2915(Sd-Im) Rd=0.50 Vf=85 err=-38NGC2955(Sb-Sbc) Rd=5.50 Vf=266 err=+53NGC2976(Sc-Scd) Rd=0.75 Vf=80 err=-37NGC3109(Sd-Im) Rd=1.40 Vf=68 err=-25NGC3521(Sb-Sbc) Rd=2.80 Vf=225 err=+5NGC3621(Sd-Im) Rd=2.10 Vf=149 err=+28NGC3726(Sc-Scd) Rd=3.00 Vf=152 err=+19NGC3741(Sd-Im) Rd=0.68 Vf=51 err=+26NGC3769(Sc-Scd) Rd=2.80 Vf=112 err=+25NGC3877(Sc-Scd) Rd=2.70 Vf=163 err=+12NGC3893(Sc-Scd) Rd=2.80 Vf=159 err=+14NGC3949(Sb-Sbc) Rd=1.40 Vf=125 err=-21NGC3953(Sb-Sbc) Rd=3.50 Vf=200 err=-6NGC3972(Sc-Scd) Rd=1.60 Vf=135 err=-27NGC3992(Sb-Sbc) Rd=3.80 Vf=242 err=-15NGC4010(Sc-Scd) Rd=1.80 Vf=128 err=-14NGC4013(Sc-Scd) Rd=2.20 Vf=185 err=+8NGC4051(Sb-Sbc) Rd=1.90 Vf=110 err=+3NGC4085(Sc-Scd) Rd=1.20 Vf=135 err=-41NGC4088(Sb-Sbc) Rd=1.90 Vf=175 err=-27NGC4100(Sb-Sbc) Rd=1.80 Vf=162 err=-28NGC4138(S0-Sa) Rd=1.30 Vf=150 err=-44NGC4157(Sb-Sbc) Rd=2.60 Vf=185 err=-1NGC4183(Sc-Scd) Rd=1.60 Vf=110 err=-36NGC4214(Sd-Im) Rd=0.50 Vf=68 err=-27NGC4217(Sb-Sbc) Rd=2.80 Vf=180 err=+5NGC4389(Sb-Sbc) Rd=1.20 Vf=110 err=-43NGC4559(Sc-Scd) Rd=3.20 Vf=123 err=+28NGC5005(Sb-Sbc) Rd=3.00 Vf=260 err=-8NGC5033(Sc-Scd) Rd=4.50 Vf=195 err=+44NGC5055(Sb-Sbc) Rd=3.50 Vf=180 err=+32NGC5371(Sb-Sbc) Rd=3.80 Vf=225 err=+14NGC5585(Sd-Im) Rd=1.50 Vf=87 err=-5NGC5907(Sc-Scd) Rd=4.20 Vf=210 err=+32NGC5985(Sb-Sbc) Rd=4.50 Vf=295 err=+2NGC6015(Sc-Scd) Rd=2.40 Vf=142 err=+6NGC6195(Sb-Sbc) Rd=5.20 Vf=260 err=+46NGC6503(Sc-Scd) Rd=2.40 Vf=121 err=+39NGC6674(Sb-Sbc) Rd=5.50 Vf=260 err=+46NGC6789(Sd-Im) Rd=0.30 Vf=60 err=-63NGC6946(Sc-Scd) Rd=2.60 Vf=180 err=+10NGC7331(Sb-Sbc) Rd=3.20 Vf=265 err=+4NGC7793(Sd-Im) Rd=1.80 Vf=118 err=-1UGC00128(Sd-Im) Rd=7.50 Vf=135 err=+80UGC02259 (Sd-Im) Rd=1.60 Vf=90 err=+0% (Sd-Im) Rd=1.60 Vf=90 err=+0% (Sd-Im)UGC02487 (S0-Sa) Rd=7.50 Vf=330 err=+41UGC02885 (Sc-Scd) Rd=8.50 Vf=290 err=+52UGC05716(Sd-Im) Rd=2.00 Vf=75 err=-6UGC05721 (Sd-Im) Rd=1.20 Vf=85 err=+0% UGC05750 (Sd-Im) Rd=8.50 Vf=290 err=+52UGC05750 (Sd-Im) Rd=4.50 Vf=80 err=+38% UGC05764 (Sd-Im) Rd=2.00 Vf=75 err=-6UGC05764(Sd-Im) Rd=0.40 Vf=57 err=-46UGC05829(Sd-Im) Rd=1.60 Vf=69 err=-10UGC06399(Sd-Im) Rd=2.50 Vf=89 err=+11UGC06446(Sd-Im) Rd=1.80 Vf=87 err=+6UGC06614 (S0-Sa) Rd=4.50 Vf=200 err=+19UGC06628 (Sd-Im) Rd=2.50 Vf=75 err=+7UGC06667(Sd-Im) Rd=2.50 Vf=90 err=+16UGC06917(Sd-Im) Rd=2.50 Vf=110 err=-12UGC06983(Sc-Scd) Rd=2.50 Vf=113 err=+4UGC07125(Sd-Im) Rd=4.50 Vf=75 err=+37UGC07151 (Sc-Scd) Rd=1.30 Vf=82 err=-32UGC07261(Sd-Im) Rd=1.10 Vf=72 err=-29UGC07399(Sd-Im) Rd=1.40 Vf=93 err=-19UGC07690(Sd-Im) Rd=0.70 Vf=62 err=-29UGC08286(Sc-Scd) Rd=1.30 Vf=84 err=-4UGC08490(Sd-Im) Rd=0.65 Vf=80 err=-29UGC08550(Sd-Im) Rd=1.50 Vf=67 err=-17UGC09037(Sc-Scd) Rd=3.50 Vf=160 err=-8UGC11455(Sc-Scd) Rd=5.50 Vf=275 err=-28UGC11557(Sd-Im) Rd=3.00 Vf=90 err=+1UGC11820(Sd-Im) Rd=4.50 Vf=90 err=+26UGCA281(Sd-Im) Rd=0.50 Vf=40 err=-36UGCA442 (Sd-Im) Rd=1.00 Vf=57 err=-44% 误差≈-31.7+31.7%。 err ≈ -31.7 +12.8-Rd 在 Rd=2.48 kpc 时为零 皮尔逊相关性r = +0.749 S0-Sa (T=0-2)Sb-Sbc (T=3-4)Sc-Scd (T=5-6)Sd-Im(T=7-10)
94 个盲星系与星盘大小的关系图,按哈勃类型着色。红线是误差与 $R_d$ 的线性回归。它在 $R_d = 2.48$ kpc 时与零点相交–这基本上就是最初校准时的圆盘大小。

误差是磁盘大小的函数

$$text{error}\,(\%) \;\approx\; -31.7 \;+\; 12.8 \times R_d \,[\text{kpc}]$$

对 94 个盲星系进行线性拟合,Pearson $r = +0.75$,残差 RMSE $= 18.4\%$.

函数形式的比较

对几种备选参数进行了比较。线性形式与对数和平方根形式在统计上没有区别:

模型皮尔逊 $r$均方根误差评论
$text{err} = a + b\,R_d$ (线性)$+0.749$$18.4\%$最简洁的分析形式
$text{err} = a + b\,\log_{10}R_d$$+0.748$$18.4\%$统计当量
$text{err} = a + b\,\sqrt{R_d}$$+0.768$$17.7\%$略有改善,没有实际收益
$text{err} = a + b\,R_d + c\,R_d^2$$17.8\%$二次项非常小 ($c大约 -1.1$)

因此,采用线性形式是对数据最简单忠实的描述。

哈勃类型沿线分布

哈勃级 $N$ 中位数 $R_d$(kpc) 中位误差 职位
S0-Sa(早期型)42.9$+0.0\%$中心,靠近零点交叉点
锑-钴(中间体)233.2$+3.9\%$居中偏右;尾部位于预测过高区域
Sc-Scd(后期螺旋)272.5$+7.7\%$遍布图表
Sd-Im(矮/不规则)401.6$-3.2\%$左侧;尾部位于预测不足的区域

图中的颜色模式并不是线性趋势的独立特征,而是通过形态轴看到的相同特征。盘状星系的哈勃序列与盘状星系的大小相关:晚期矮星系主要是紧凑型,中期螺旋星系主要是大型。因此,每种颜色都位于回归线的不同位置,Sd-Im 在左边,Sc-Scd 在中间,Sb-Sbc 在右边。

结构性残差,而非随机噪声

与单个物理参数线性相关的散度,在校准点处越过零,是模型关系中一个缺失的加常数,而不是随机观测散度。这种偏差是可以纠正的:它可以被相干长度定律中的一个额外自由度所吸收。

3.十个被过度预测的星系

这些星系BeeTheory 预测其平面旋转速度高于观测值的星系。按残差大小排序:

银河系哈勃类型R_d$ (kpc)$M_\star/10^{10}$$f_\text{gas}$$\Sigma_d$V_f$$V_\text{tot}$错误
UGC00128Sd-Im7.501.060.3960135243+80.0%
NGC0801Sb-Sbc5.802.010.32190208326+56.6%
NGC2955Sb-Sbc5.503.990.23420266406+52.7%
UGC02885Sc-Scd8.503.400.41150290441+52.0%
NGC0925Sc-Scd3.100.220.7572105155+48.0%
NGC6195Sb-Sbc5.203.400.26400260380+46.3%
NGC6674Sb-Sbc5.503.330.29350260380+46.2%
NGC5033Sb-Sbc4.501.270.46200195280+43.7%
UGC02487S0-Sa7.505.300.23300330465+40.8%
NGC6503Sc-Scd2.400.380.55210121168+38.9%
财产中位值范围比较
$R_d$4.5 千兆位点2.4 – 8.5比中位数大 2 美元/次
$M_\star$1.3美元乘以10^{10}\,M_\odot$2.2 \times 10^{9}$ – 5.3 \times 10^{10}$8 美元/次,更大规模
$f_\text{gas}$$0.41$$0.23$ – $0.87$低于中位数 (0.64)
哈勃 $T$5 美元(Sbc)$1$ – $8$集中在中间螺旋
V_f$195 美元 km/s$69$ – $330$样品中最快的旋转器

超预测组的轮廓

大质量中型螺旋星系。这些星系位于回归线的右侧,远远高于零交叉点。模型的相干长度定律 $\ell = c_text\{disk}\,R_d$产生的$\ell$值在这一范围内超过了20 kpc,产生的波场质量超过了观测到的旋转所需的质量。

4.预测最不充分的十个星系

这些是BeeTheory预测的平旋转速度低于观测值的星系。按残差大小排序:

银河系哈勃类型R_d$ (kpc)$M_\star/10^{10}$$f_\text{gas}$$\Sigma_d$V_f$$V_\text{tot}$错误
NGC6789Sd-Im0.300.010.532506022-63.0%
UGC05764Sd-Im0.400.000.86805731-45.6%
UGCA442Sd-Im1.000.000.85155732-44.2%
NGC4138S0-Sa1.300.130.3325015085-43.6%
NGC4389Sb-Sbc1.200.070.3715011062-43.4%
NGC4085Sb-Sbc1.200.090.4220013579-41.1%
NGC2915Sd-Im0.500.010.841608553-38.2%
NGC2976Sb-Sbc0.750.040.292208050-37.4%
NGC4183Sc-Scd1.600.030.814011070-36.3%
UGCA281Sd-Im0.500.010.63804026-36.1%
财产中位值范围比较
$R_d$1.1 千兆位点0.30 – 1.802美元/次$ 小于中位数
$M_\star$2.7 *times 10^{8}\,M_\odot$$4 times 10^{7}$ – $1.3 times 10^{9}$6 美元/次,规模较小
$f_\text{gas}$$0.58$$0.29$ – $0.86$低于中位数 (0.64)
哈勃 $T$8$ (Sd)$1$ – $10$集中于晚型矮星
V_f$82$ 千米/秒$40$ – $150$慢速旋转器

预测不足的星团的轮廓

紧凑的低质量矮星系和小螺旋星系。这些星系位于回归线的左侧,远远低于零交叉点。相干长度定律 $\ell = c_\text{disk}\,R_d$ 在这种情况下产生的 $\ell$ 数量级为 1$$-$$3 kpc,可能太短了,不足以收集波场的全部范围。

5.三个组别的并排比较

物业(中位数) 预测过高
(err > +30%, $N = 15$)
预测良好
(|err| ≤ 30%, $N = 67$)
预测不足
(err < -30%, $N = 12$)
R_d$ (kpc)4.52.41.1
$M_\star / 10^{10}$1.270.150.027
$M_\text{gas}/ 10^{10}$0.930.270.04
$f_\text{gas}$0.410.640.58
$\Sigma_d$200140115
哈勃 $T$5 (Sbc)6 (Sc)8 (Sd)
V_f$ (千米/秒)19511382

从左到右,每个属性都是单调变化的。预测过高的星团体积更大、质量更大、恒星更多,旋转速度更快;预测过低的星团体积更小、质量更轻、气体更多,旋转速度更慢;预测良好的大多数星团则介于两者之间。银河系($R_d = 2.6$ kpc,$V_f 约 230$ km/s)自然属于校准锚定的良好预测范围。

6.口译

该模型有一个耦合参数$\lambda$和三个通用几何常数$(c_\text{disk}, c_\text{sph}, c_\text{arm})$。这些常数是在一个中等大小的星系(银河系,$R_d = 2.6$ kpc)上确定的,并在22个大小范围相似的星系上得到了验证。注释 IX 的盲测结果表明,它们的通用性相当好,但残差随圆盘大小呈线性漂移。

仿射修正就足够了

在 $R_d$ 中的残差线性–在 $R_d = 2.48$ kpc 处有一条直线过零–是相干长度关系中缺失的加法偏移的特征。当前的定律$\ell = c_\text{disk}\,R_d$严格地将波相干长度与磁盘尺度成比例地联系在一起。用仿射关系$ell = c_\text{disk}(R_d – R_0)$来代替它,其中$R_0$是一个约为2.5$ kpc的小偏移量,将产生一个在校准点消失并在两侧线性增长的残差–这正是观测到的模式。

预测良好的大多数具有广泛代表性

三分之二的样本属于预测良好的范围。这 67 个星系涵盖了所有的哈勃类型,恒星质量的系数为 100 美元。模型的有效范围并不狭窄:它涵盖了 SPARC 星群的大部分,偏差集中在星盘大小的两个极端,这与线性残差 $R_d$ 的结果完全一致。

7.摘要

1.94个星系盲测的预测误差在磁盘尺度长度上呈明显的线性趋势:$\text{error}(\%) \approx -31.7 + 12.8\,R_d$, Pearson $r = +0.75$,残差的RMSE $= 18.4\%$.

2.线性回归线在 $R_d = 2.48$ kpc 处交叉为零,这基本上是校准所依据的银河系圆盘大小。线的两端对应两个物理上截然不同的离群点。

3.15个高估超过+30%$的星系都是大质量的中间型螺旋星系:中位数$R_d = 4.5$ kpc,$M_star (约10^{10}\,M_\odot$),$V_f (约200$ km/s)。

4.这12个低估了30%以上的星系都是紧凑的低质量矮星系:中位数$R_d = 1.1$ kpc,$M_star (约为10^{8}的3倍),$M_odot (约为80$ km/s)。

5.这种偏差可以通过对相干长度定律进行仿射修正来吸收,即 $\ell = c_text\{disk}(R_d – R_0)$ ,其中 $R_0 \approx 2.5$ kpc – 引入一个新常数。


参考文献Lelli, F., McGaugh, S. S., Schombert, J. M. –SPARC:Mass Models for 175 DiskGalaxies with Spitzer Photometry and Accurate Rotation Curves, AJ 152, 157 (2016).- de Vaucouleurs, G. et al. –Third Reference Catalogue of Bright Galaxies, Springer (1991).- McGaugh, S. S. –The third law of galactic rotation, Galaxies 2, 601 (2014).- Dutertre, X. –Bee Theory™:基于波的引力建模,v2,BeeTheory.com(2023 年)。

BeeTheory.com – 波基量子引力 – SPARC残差 – © Technoplane S.A.S. 2026