BeeTheory — Основы — Техническая заметка X

Анатомия остатков:
Линейная тенденция изменения размера диска

Слепой тест 94 галактик, приведенный в Примечании IX, показал систематическую остаточную тенденцию к увеличению размера диска. В этой заметке дается количественная характеристика этой тенденции, выделяются самые большие отклонения с каждой стороны и определяется структурное происхождение дисперсии.

1. Результат первый

Линейный остаток, две противоположные популяции

Ошибка предсказания линейно зависит от длины диска: $\text{error}\,(\%)\approx -31.7 + 12.8\,R_d$, с корреляцией Пирсона $r = +0.75$. Линия пересекает ноль при $R_d = 2,48$ кпк, что соответствует размеру диска Млечного Пути, на котором была основана калибровка. Две крайние точки этой регрессии соответствуют двум физически различным популяциям выбросов: крупные массивные спирали (завышенные прогнозы) на одном конце, компактные карлики (заниженные прогнозы) — на другом.

2. Остаток является линейным в $R_d$

При построении графика зависимости ошибки предсказания от $R_d$, где каждая точка окрашена в цвет типа Хаббла, сразу видна линейность тенденции. Красная линия — это линейная регрессия ошибки на $R_d$ по всем 94 слепым галактикам.

Ошибка предсказания в зависимости от размера диска — линейная в $R_d$, окрашенная по типу Хаббла 0.30.5123510 -50%-25%+0%+25%+50%+75% → Область завышенных прогнозов Диапазон ±30% ← Область с заниженным прогнозом Длина шкалы диска $R_d$ (кпк) — логарифмическая шкала Ошибка предсказания (%) F579-V1 (Sd-Im) Rd=3.20 Vf=105 err=-12%F583-1 (Sd-Im) Rd=1.80 Vf=83 err=-30%F583-4 (Sd-Im) Rd=1.40 Vf=67 err=-29%IC2574 (Sd-Im) Rd=2.80 Vf=69 err=+35%KK98-251 (Sd-Im) Rd=0.30 Vf=17 err=-23%M33 (Sc-Scd) Rd=1.40 Vf=100 err=-2%NGC0055 (Sd-Im) Rd=1.80 Vf=87 err=-1%NGC0100 (Sc-Scd) Rd=2.30 Vf=83 err=+13%NGC0247 (Sd-Im) Rd=2.40 Vf=90 err=+20%NGC0289 (Sb-Sbc) Rd=3.50 Vf=155 err=+29%NGC0300 (Sd-Im) Rd=1.50 Vf=76 err=+0%NGC0801 (Sc-Scd) Rd=5.80 Vf=208 err=+57%NGC0891 (Sb-Sbc) Rd=4.10 Vf=212 err=+7%NGC0925 (Sd-Im) Rd=3.10 Vf=105 err=+48%NGC1003 (Sc-Scd) Rd=2.80 Vf=115 err=+12%NGC1090 (Sb-Sbc) Rd=3.80 Vf=170 err=+17%NGC1705 (S0-Sa) Rd=0.60 Vf=54 err=-19%NGC2366 (Sd-Im) Rd=1.30 Vf=55 err=+14%NGC2403 (Sc-Scd) Rd=1.80 Vf=131 err=-4%NGC2683 (Sb-Sbc) Rd=2.90 Vf=175 err=+15%NGC2903 (Sb-Sbc) Rd=2.60 Vf=184 err=-0%NGC2915 (Sd-Im) Rd=0.50 Vf=85 err=-38%NGC2955 (Sb-Sbc) Rd=5.50 Vf=266 err=+53%NGC2976 (Sc-Scd) Rd=0.75 Vf=80 err=-37%NGC3109 (Sd-Im) Rd=1.40 Vf=68 err=-25%NGC3521 (Sb-Sbc) Rd=2.80 Vf=225 err=+5%NGC3621 (Sd-Im) Rd=2.10 Vf=149 err=+28%NGC3726 (Sc-Scd) Rd=3.00 Vf=152 err=+19%NGC3741 (Sd-Im) Rd=0.68 Vf=51 err=+26%NGC3769 (Sc-Scd) Rd=2.80 Vf=112 err=+25%NGC3877 (Sc-Scd) Rd=2.70 Vf=163 err=+12%NGC3893 (Sc-Scd) Rd=2.80 Vf=159 err=+14%NGC3949 (Sb-Sbc) Rd=1.40 Vf=125 err=-21%NGC3953 (Sb-Sbc) Rd=3.50 Vf=200 err=-6%NGC3972 (Sc-Scd) Rd=1.60 Vf=135 err=-27%NGC3992 (Sb-Sbc) Rd=3.80 Vf=242 err=-15%NGC4010 (Sc-Scd) Rd=1.80 Vf=128 err=-14%NGC4013 (Sc-Scd) Rd=2.20 Vf=185 err=+8%NGC4051 (Sb-Sbc) Rd=1.90 Vf=110 err=+3%NGC4085 (Sc-Scd) Rd=1.20 Vf=135 err=-41%NGC4088 (Sb-Sbc) Rd=1.90 Vf=175 err=-27%NGC4100 (Sb-Sbc) Rd=1.80 Vf=162 err=-28%NGC4138 (S0-Sa) Rd=1.30 Vf=150 err=-44%NGC4157 (Sb-Sbc) Rd=2.60 Vf=185 err=-1%NGC4183 (Sc-Scd) Rd=1.60 Vf=110 err=-36%NGC4214 (Sd-Im) Rd=0.50 Vf=68 err=-27%NGC4217 (Sb-Sbc) Rd=2.80 Vf=180 err=+5%NGC4389 (Sb-Sbc) Rd=1.20 Vf=110 err=-43%NGC4559 (Sc-Scd) Rd=3.20 Vf=123 err=+28%NGC5005 (Sb-Sbc) Rd=3.00 Vf=260 err=-8%NGC5033 (Sc-Scd) Rd=4.50 Vf=195 err=+44%NGC5055 (Sb-Sbc) Rd=3.50 Vf=180 err=+32%NGC5371 (Sb-Sbc) Rd=3.80 Vf=225 err=+14%NGC5585 (Sd-Im) Rd=1.50 Vf=87 err=-5%NGC5907 (Sc-Scd) Rd=4.20 Vf=210 err=+32%NGC5985 (Sb-Sbc) Rd=4.50 Vf=295 err=+2%NGC6015 (Sc-Scd) Rd=2.40 Vf=142 err=+6%NGC6195 (Sb-Sbc) Rd=5.20 Vf=260 err=+46%NGC6503 (Sc-Scd) Rd=2.40 Vf=121 err=+39%NGC6674 (Sb-Sbc) Rd=5.50 Vf=260 err=+46%NGC6789 (Sd-Im) Rd=0.30 Vf=60 err=-63%NGC6946 (Sc-Scd) Rd=2.60 Vf=180 err=+10%NGC7331 (Sb-Sbc) Rd=3.20 Vf=265 err=+4%NGC7793 (Sd-Im) Rd=1.80 Vf=118 err=-1%UGC00128 (Sd-Im) Rd=7.50 Vf=135 err=+80%UGC02259 (Sd-Im) Rd=1.60 Vf=90 err=+0%UGC02487 (S0-Sa) Rd=7.50 Vf=330 err=+41%UGC02885 (Sc-Scd) Rd=8.50 Vf=290 err=+52%UGC05716 (Sd-Im) Rd=2.00 Vf=75 err=-6%UGC05721 (Sd-Im) Rd=1.20 Vf=85 err=+0%UGC05750 (Sd-Im) Rd=4.50 Vf=80 err=+38%UGC05764 (Sd-Im) Rd=0.40 Vf=57 err=-46%UGC05829 (Sd-Im) Rd=1.60 Vf=69 err=-10%UGC06399 (Sd-Im) Rd=2.50 Vf=89 err=+11%UGC06446 (Sd-Im) Rd=1.80 Vf=87 err=+6%UGC06614 (S0-Sa) Rd=4.50 Vf=200 err=+19%UGC06628 (Sd-Im) Rd=2.50 Vf=75 err=+7%UGC06667 (Sd-Im) Rd=2.50 Vf=90 err=+16%UGC06917 (Sd-Im) Rd=2.50 Vf=110 err=-12%UGC06983 (Sc-Scd) Rd=2.50 Vf=113 err=+4%UGC07125 (Sd-Im) Rd=4.50 Vf=75 err=+37%UGC07151 (Sc-Scd) Rd=1.30 Vf=82 err=-32%UGC07261 (Sd-Im) Rd=1.10 Vf=72 err=-29%UGC07399 (Sd-Im) Rd=1.40 Vf=93 err=-19%UGC07690 (Sd-Im) Rd=0.70 Vf=62 err=-29%UGC08286 (Sc-Scd) Rd=1.30 Vf=84 err=-4%UGC08490 (Sd-Im) Rd=0.65 Vf=80 err=-29%UGC08550 (Sd-Im) Rd=1.50 Vf=67 err=-17%UGC09037 (Sc-Scd) Rd=3.50 Vf=160 err=-8%UGC11455 (Sc-Scd) Rd=5.50 Vf=275 err=-28%UGC11557 (Sd-Im) Rd=3.00 Vf=90 err=+1%UGC11820 (Sd-Im) Rd=4.50 Vf=90 err=+26%UGCA281 (Sd-Im) Rd=0.50 Vf=40 err=-36%UGCA442 (Sd-Im) Rd=1.00 Vf=57 err=-44% err ≈ -31.7 +12.8-Rd ноль при Rd=2,48 кпк Корреляция Пирсона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$ кпк — по сути, это размер диска, на котором основывалась оригинальная калибровка.

Ошибка как функция размера диска

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

Линейная подгонка по 94 слепым галактикам, Pearson $r = +0,75$, RMSE остатков $= 18,4\%$.

Сравнение функциональных форм

Сравнивались несколько альтернативных параметров. Линейная форма статистически неотличима от логарифмических и квадратно-корневых альтернатив:

МодельПирсон $r$RMSEКомментарий
$\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 \approx -1.1$).

Поэтому линейная форма принимается как наиболее простое и верное описание данных.

Распределение типов Хаббла вдоль линии

Класс Хаббла $N$ Медиана $R_d$ (кпк) Медианная ошибка Позиция
S0-Sa (ранний тип)42.9$+0.0\%$Центр, около нулевого пересечения
Sb-Sbc (промежуточный)233.2$+3.9\%$Справа от центра; хвост в области завышенных прогнозов
Sc-Scd (поздняя спираль)272.5$+7.7\%$Распределите по диаграмме
Sd-Im (карликовый / неправильный)401.6$-3.2\%$Левая сторона; хвост в области заниженных прогнозов

Цветовая картина на рисунке не является независимым признаком линейного тренда — это тот же самый признак, который виден через ось морфологии. Последовательность Хаббла в дисковых галактиках коррелирует с размером диска: карлики позднего типа преимущественно компактны, а промежуточные спирали — преимущественно велики. Поэтому каждый цвет располагается на своем участке линии регрессии: Sd-Im — слева, Sc-Scd — в центре, а Sb-Sbc — справа.

Структурный остаток, а не случайный шум

Разброс, линейно зависящий от одного физического параметра и пересекающий ноль в точке калибровки, является признаком отсутствия аддитивной константы в одном из соотношений модели, а не случайного разброса наблюдений. Отклонение можно исправить: оно может быть поглощено одной дополнительной степенью свободы в законе длины когерентности.

3. Десять галактик с самыми завышенными прогнозами

Это те галактики, для которых BeeTheory предсказывает скорость плоского вращения выше, чем наблюдается. Отсортированы по величине остатка:

GalaxyТип Хаббла$R_d$ (кпк)$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 \times 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$ км/с$69$ — $330$Самые быстрые ротаторы в образце

Профиль группы с завышенными прогнозами

Крупные, массивные спирали промежуточного типа. Эти галактики находятся в правой части линии регрессии, намного выше нулевого пересечения. Закон длины когерентности модели $\ell = c_\text{disk}\,R_d$ приводит к значениям $\ell$ выше 20 кпк в этом режиме, генерируя больше массы волнового поля, чем требуется для наблюдаемого вращения.

4. Десять самых недооцененных галактик

Это те галактики, для которых BeeTheory предсказывает плоскую скорость вращения ниже наблюдаемой. Отсортированы по величине остатка:

GalaxyТип Хаббла$R_d$ (кпк)$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.80В $2\ раза$ меньше медианы
$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$ кпк в этом режиме, возможно, слишком короткий, чтобы собрать весь объем волнового поля.

5. Боковое сравнение трех групп

Недвижимость (медиана) Завышенные прогнозы
(err > +30%, $N = 15$)
Хорошо предсказанный
(|err| ≤ 30%, $N = 67$)
Недостаточно предсказанный
(err < -30%, $N = 12$)
$R_d$ (кпк)4.52.41.1
$M_\star / 10^{10}$1.270.150.027
$M_\text{газ} / 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$ кпк, $V_f около 230$ км/с) естественным образом попадает в хорошо предсказанный режим, на котором была закреплена калибровка.

6. Интерпретация

Модель имеет единственный параметр связи $\lambda$ и три универсальные геометрические константы $(c_\text{disk}, c_\text{sph}, c_\text{arm})$. Они были определены на галактике промежуточного размера ( Млечный Путь, $R_d = 2,6$ кпк) и подтверждены на двадцати двух галактиках аналогичного диапазона размеров. Слепой тест в Примечании IX показывает, что они достаточно хорошо обобщаются, но с остатком, который линейно дрейфует с размером диска.

Достаточно аффинной коррекции

Линейность остатка в $R_d$ — хорошо подогнанная единственная прямая линия, пересекающая ноль при $R_d = 2,48$ кпк — является признаком отсутствия аддитивного смещения в соотношении длины когерентности. Текущий закон $\ell = c_\text{диск}\,R_d$ связывает длину когерентности волн строго пропорционально масштабу диска. Замена его аффинным соотношением $\ell = c_\text{disk}(R_d — R_0)$, где $R_0$ — небольшое смещение около 2,5$ кпк, приведет к остатку, который исчезает в точке калибровки и линейно растет по обе стороны — именно такая картина наблюдается.

Хорошо предсказанное большинство является в целом репрезентативным

Две трети выборки попадают в хорошо предсказанный диапазон. Эти 67 галактик охватывают весь диапазон хаббловских типов и фактор $\sim 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$ кпк, что соответствует размеру диска Млечного Пути, на котором основывалась калибровка. Два конца линии соответствуют двум физически различным популяциям выбросов.

3. 15 галактик, предсказания которых превышены более чем на $+30\%$, являются крупными, массивными спиралями промежуточного типа: медиана $R_d = 4,5$ кпк, $M_\star \approx 10^{10}\,M_\odot$, $V_f \approx 200$ км/с.

4. 12 галактик, недооцененных более чем на $-30\%$, являются компактными карликами с низкой массой: медиана $R_d = 1,1$ кпк, $M_\star \approx 3 \times 10^{8}\,M_\odot$, $V_f \approx 80$ км/с.

5. Отклонение устраняется аффинной поправкой к закону длины когерентности, $\ell = c_\text{disk}(R_d — R_0)$, при $R_0 \approx 2,5$ kpc — введение одной новой константы.


Ссылки. Lelli, F., McGaugh, S. S., Schombert, J. M. — SPARC: модели масс для 175 дисковых галактик с помощью спитцеровской фотометрии и точных кривых вращения, AJ 152, 157 (2016). — Де Вокулерс, Г. и др. — Третий справочный каталог ярких галактик, Шпрингер (1991 г.). — McGaugh, S. S. — Третий закон галактического вращения, Galaxies 2, 601 (2014). — Дютертр, X. — Bee Theory™: Волновое моделирование гравитации, v2, BeeTheory.com (2023).

BeeTheory.com — Квантовая гравитация на основе волн — остатки SPARC — © Technoplane S.A.S. 2026