BeeTheory — статистический анализ — 2025

Чтение фита

χ²/dof = 0.475

Число от 0 до ∞, которое говорит Вам, хорошо ли Ваша модель соответствует данным, слишком хорошо или совсем не соответствует. Вот что оно означает в галактической симуляции BeeTheory — и что нам нужно, чтобы получить действительно достоверный результат.

Что такое χ²/dof?

Когда модель предсказывает значения _Vmodel(_Ri), а данные дают наблюдаемые значения _Vobs(_Ri) с неопределенностями измерений _σi, сокращенный хи-квадрат равен:

\(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{model}}(R_i)-V_{\mathrm{obs}}(R_i)}{\sigma_i}\right)^2\)

где N — количество точек данных, а p — количество свободных параметров.

Здесь: N = 16 точек Gaia 2024, p = 2 установленных параметра, K и α, поэтому dof = 14.

Каждый член в сумме — это тяга: остаток, выраженный в единицах неопределенности измерения.

Отклонение в 0,5 означает, что модель расходится на полсигмы — отлично. Тяга в 2,0 означает расхождение в две сигмы — стоит исследовать.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{модель подходит к данным на уровне их неопределенностей}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{модель подходит слишком хорошо: неопределенности могут быть переоценены, или overfitting}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{модель неверна, или неопределенности недооценены}\)

Наше число: 0.475

0.475

χ² / dof — N = 16 — p = 2 — dof = 14

Статистически превосходно

Модель не ошибается

Погрешности могут быть несколько большими

Шкала χ²/доф

0,475 ← мы

1.0 идеальный

00.511.522.53+

< 0.3
чрезмерно подходящий

0.3-0.7
отлично

0.7-1.3
хорошо

1.3-2.0
приемлемо

> 2.0
плохо

Значение 0,475 находится в отличной зоне. Модель не ошибается — средний остаток составляет всего 0,69σ по всем 16 точкам данных.

Физически это означает, что BeeTheory предсказывает круговую скорость с точностью менее одной погрешности измерения практически на каждом наблюдаемом радиусе.

Но «превосходный» — это не то же самое, что «проверенный».

χ²/dof = 0,475 может также означать, что погрешности измерений _σi немного завышены. Если бы истинные погрешности были на 30% меньше, та же модель дала бы χ²/dof ≈ 1,0.

Мы не можем отличить «модель очень хороша» от «ошибки слегка щедрые» на основании одного лишь χ². Это стандартная статистическая неоднозначность.

The Residuals — Point by Point

Необработанное число 0,475 скрывает информацию. Взгляд на отдельные вытягивания раскрывает структуру подгонки: в каких точках модель приживается, а где испытывает трудности.

График: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) для каждой точки данных Gaia 2024

R (кпк)	Vobs	σ	Vmodel	Остаток	Потяните
Загрузка остатков…

15 / 16 очков в пределах 1σ

Каждая точка данных, кроме R = 27,3 кпк, имеет |pull| < 1,0. В хорошо откалиброванной модели с гауссовыми ошибками мы ожидаем около 68% точек в пределах 1σ - здесь мы имеем 94%.

Это говорит о том, что модель превосходна, или же планки погрешностей на внутренних точках немного великоваты.

Выброс в точке R = 27,3 кпк — притяжение = +1,53σ

Модель предсказывает _Vc = 195,3 км/с, в то время как Gaia измеряет 173 ± 17 км/с.

Расхождение составляет 22,3 км/с, или 1,53σ. Это не является статистически тревожным, но физически значимо: модель слишком медленно снижается при самом большом радиусе.

Наилучший подходящий параметр K = 0.01349

Константа связи K = 0,01349 кпк-¹ — это амплитуда поля темной массы BeeTheory, генерируемой на единицу видимой массы.

Он был определен путем минимизации χ² по 16 точкам данных Gaia с α = 0,0744 кпк-¹, зафиксированным в комбинированном подходе Newby + Gaia.

\(K=0.01349\,\mathrm{kpc}^{-1}\) \(\lambda=K\ell^2=\frac{K}{\alpha^2}=\frac{0.01349}{(0.0744)^2}=2.44\) \(\rho_{\mathrm{dark}}(R_\odot=8\,\mathrm{kpc})=0.37\,\mathrm{GeV/cm^3}\)

Наблюдаемая локальная плотность темной материи составляет:

\(\rho_{\mathrm{obs}}(R_\odot)=0.39\pm0.03\,\mathrm{GeV/cm^3}\)

Безразмерная связь λ = 2,44 находится в диапазоне λ ∈ [2, 7], наблюдаемом во всех вариантах BeeTheory, включая однокомпонентные, многокомпонентные и атомарные H₂ сравнения.

Эта универсальность — один и тот же порядок величин от фемтометровых масштабов до масштабов килопарсека — является самой сильной проверкой внутренней согласованности системы BeeTheory.

Почему K здесь меньше, чем в однодисковых посадках

Предыдущие расчеты, использовавшие в качестве источника только тонкий диск, дали _Ksingle = 0,038 кпк-¹. С учетом того, что все шесть галактических компонентов вносят свой вклад в темное поле, общий барионный источник в 2,8 раза больше.

Поэтому K должен быть пропорционально меньше, чтобы создать ту же самую темную массу.

\(K_{\mathrm{multi}}\approx\frac{K_{\mathrm{single}}}{2.8}=\frac{0.038}{2.8}=0.0136\,\mathrm{kpc}^{-1}\)

Это именно то, что дает подгонка. Это проверка соответствия, а не совпадение.

Почему 27,3 кпк — это сложно для любой модели

Самая крайняя точка данных Gaia 2024, _Vc(27,3 кпк) = 173 ± 17 км/с, сложна не только для BeeTheory. Это самая сложная точка для любой модели темной материи, применяемой к кривой вращения Млечного Пути.

Модель	Vc(27.3 кпк) предсказано	Остаток	Потяните	χ²/доф
BeeTheory	195,3 км/с	+22.3	+1.53σ	0.475
Профиль NFW	~198 км/с	+25	+1.5σ	0.44
Эйнасто α = 0.91	~191 км/с	+18	+1.1σ	0.38
Изотермический ореол	~218 км/с	+45	+2.6σ	2.6

Ни одна стандартная двухпараметрическая модель не воспроизводит _Vc(27,3 кпк) = 173 км/с в пределах 1σ.

Есть три правдоподобных объяснения.

Само измерение: R = 27,3 кпк — самая удаленная точка, с наименьшим количеством кинематических трассеров и наибольшей систематической неопределенностью. Gaia DR4 может сместить это значение.
Вклад газового диска: Диск HI простирается примерно до 25 кпк и вносит свой вклад в барионную массу на больших радиусах. Включение его в качестве отдельного компонента может сделать спад более крутым.
Меньшая длина когерентности: α = 0,12 кпк-¹, или ℓ = 8,3 кпк, лучше соответствует самой внешней точке, но ухудшает состояние внутреннего плато.

Как должен выглядеть по-настоящему надежный результат

Текущее соответствие статистически хорошо. Но хорошее соответствие в модели с 2 свободными параметрами и 16 точками данных — это не то же самое, что доказанное.

0.475

Текущий χ²/dof, 16 точек, 2 параметра

< 0.8

Цель: хорошее соответствие расширенному набору данных

~1.0

Идеальный вариант: калибровка по образцам галактик

Требования к надежной проверке

✓ Хорошее χ²/dof на обучающем наборе данных: Достигнуто. χ²/dof = 0,475 на Gaia 2024.
✓ Исправьте локальную плотность темной материи: Достигнута. 0,37 ГэВ/см³ предсказано против 0,39 ± 0,03 наблюдается.
✓ Согласованное безразмерное соединение: Достигнуто. λ = 2,44 по сравнению с λ ≈ 3,5 из калибровки H₂.
! Воспроизведите крайнюю точку Gaia в пределах 1σ: Не полностью достигнуто. Остаток при R = 27,3 кпк составляет +1,53σ.
✗ Слепое предсказание по другим галактикам: Еще не сделано. Каталог SPARC обеспечивает естественную проверку.
✗ Вывод ℓ из первых принципов: Еще не сделано. В настоящее время ℓ устанавливается, а не выводится.
✗ Проверка в масштабах кластера: Еще не все сделано. Линзирование пулевых скоплений — ключевой тест.
→ Расширенная кривая вращения за пределами 30 кпк: Gaia DR4 должна обеспечить критический тест в ближайшем будущем.

Честный научный статус

BeeTheory в ее нынешнем виде — это физически мотивированная феноменологическая схема, которая соответствует кривой вращения Млечного Пути, а также лучшим эмпирическим профилям темной материи, с тем преимуществом, что у нее есть физический механизм.

Это еще не теория в строгом смысле слова, потому что K и ℓ скорее подогнаны, чем выведены.

Путь к полной теории ясен: выведите ℓ = f(_Lsource) из постулата о волновой массе, проведите универсальное тестирование на SPARC и сделайте слепое предсказание для Gaia DR4.

Данные: Ou, X. et al., MNRAS 528, 2024. Модель: BeeTheory v2, Dutertre 2023, многокомпонентная подгонка. Ссылка на SPARC: Lelli, McGaugh, Schombert, AJ 152, 2016. Кластер пули: Clowe et al., ApJL 648, 2006.