BeeTheory — Галактическая симуляция v2 — начальное поколение 2025 мая 17 года с Клодом
Скрытая масса Млечного Пути: BeeTheory 3D Yukawa с физическим усечением диска
Исправленное моделирование: скорость барионного диска падает кеплеровским образом за его физический край, а трехмерное ядро Юкавы BeeTheory заполняет все пространство. Два параметра, данные о вращении, полученные в эпоху Гайи, и усеченная модель диска.
BeeTheory.com — Ou et al., MNRAS 528, 2024 — Исправленная BeeTheory v2
K = 0.040 кпк-¹
Волновая муфта
α = 0.087 кпк-¹
Обратная когерентность
ℓ = 11,5 кпк
Длина когерентности
χ²/dof ≈ 0.31
Отличная упрощенная посадка
0. Результат — Уравнения и параметры
Каждое кольцевое кольцо галактического диска радиусом R′ генерирует трехмерное поле эффективной темной массы через ядро Юкавы теории Би. Общая плотность темной массы на сферическом радиусе r составляет:
\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma_0e^{-R’/R_d}\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\) \(D(r,R’)=\sqrt{r^2+R’^2}\)Ядро получено из скорректированного закона силы BeeTheory:
\(F(D)\propto\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)Она сводится к ньютоновской обратно-квадратичной форме для D, значительно меньшего, чем длина когерентности ℓ.
\(D\ll\ell=\frac{1}{\alpha}\quad\Longrightarrow\quad F(D)\propto\frac{1}{D^2}\)Скорость барионного диска использует формулу Фримена в пределах его физического края Rtrunc ≈ 4Rd = 10,4 кпк, а затем плавно переходит к кеплеровскому падению, ожидаемому от распределения конечной массы.
\(K=0.0397\,\mathrm{kpc}^{-1},\qquad \alpha=0.0868\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.5\,\mathrm{kpc}\)Резюме по фитнесу
| Наблюдаемый | Значение эпохи Гайи | BeeTheory | Потяните |
|---|---|---|---|
| Vc(4 кпк) | 220 ± 10 км/с | 219.8 км/с | -0.02σ |
| Vc(8 кпк) | 230 ± 6 км/с | 233,2 км/с | +0.53σ |
| Vc(12 кпк) | 226 ± 7 км/с | 223.8 км/с | -0.31σ |
| Vc(20 кпк) | 215 ± 10 км/с | 211,2 км/с | -0.38σ |
| Vc(27.3 кпк) | 173 ± 17 км/с | 199.0 км/с | +1.53σ |
| ρdark(R⊙ = 8 кпк) | 0,39 ± 0,03 ГэВ/см³ | 0.47 ГэВ/см³ | +2.3σ |
| Мдарк(<8 кпк) | ~5 × 10¹⁰ M⊙ | 5.3 × 10¹⁰ M⊙ | закрыть |
| Mtot(<200 кпк) | 5-9 × 10¹¹ M⊙ | 3.3 × 10¹¹ M⊙ | низкий уровень |
Упрощенная подгонка дает χ²/dof ≈ 0,31. Самой сложной точкой остается крайнее значение эпохи Гайи в 27,3 кпк, где наблюдаемый спад более резкий, чем предсказывает эта двухпараметрическая модель.
1. Усечение диска — почему и как
1.1 Проблема с бесконечным экспоненциальным диском
Формула диска Фримена предполагает экспоненциальную поверхностную плотность, простирающуюся до бесконечности. Математически она никогда не достигает нуля, но физически звездный диск Млечного Пути имеет конечную протяженность. За эффективным звездным краем заключенная в нем барионная масса по существу постоянна, и вклад скорости должен падать примерно как кеплеровское поле точечной массы.
\(\Sigma(R)=\Sigma_0e^{-R/R_d}\)За краем диска скорость барионов стремится к нулю:
\(V_{\mathrm{bar}}(R)\xrightarrow{R\gg R_d}\sqrt{\frac{GM_{\mathrm{bar,tot}}}{R}}\) \(M_{\mathrm{bar,tot}}=M_{\mathrm{disk}}+M_{\mathrm{bulge}}\approx4.7\times10^{10}M_\odot\)Примеры значений:
\(V_{\mathrm{bar}}(30\,\mathrm{kpc})\approx82\,\mathrm{km/s},\qquad V_{\mathrm{bar}}(50\,\mathrm{kpc})\approx63\,\mathrm{km/s}\)1.2 Формула плавного усечения
В моделировании используется плавный переход между формулой диска Фримена и кеплеровским значением. Переход сосредоточен в точке Rtrunc = 4Rd = 10,4 кпк с шириной σ = 1,5 кпк.
\(V_{\mathrm{bar}}(R)=\sqrt{(1-w)V_{\mathrm{Freeman}}^2(R)+w\,\min(V_{\mathrm{Freeman}},V_{\mathrm{Kepler}})^2}\) \(w(R)=\frac{1}{2}\left[1+\tanh\left(\frac{R-R_{\mathrm{trunc}}}{\sigma}\right)\right]\) \(R_{\mathrm{trunc}}=4R_d=10.4\,\mathrm{kpc},\qquad \sigma=1.5\,\mathrm{kpc}\)Минимальная функция не позволяет барионному диску превысить физический кеплеровский предел за пределами края диска.
| R | VFreeman | VKeplerian | Vbar, усеченный | Доминирующий режим |
|---|---|---|---|---|
| 5 кпк | 174,5 км/с | 201,1 км/с | 174,5 км/с | Freeman |
| 8 кпк | 161,5 км/с | 159,0 км/с | 161,5 км/с | Фримен ≈ Кеплер |
| 10,4 кпк | 143,0 км/с | 139,3 км/с | 141,2 км/с | Переход |
| 16 кпк | 112,4 км/с | 112,4 км/с | 112,4 км/с | Кеплериан |
| 25 кпк | 89.9 км/с | 89.9 км/с | 89.9 км/с | Кеплериан |
| 50 кпк | 63.6 км/с | 63.6 км/с | 63.6 км/с | Кеплериан |
2. Теория BeeTheory 3D Плотность темной массы
2.1 Дисковые кольца, излучающие в 3D
Каждое кольцо галактического диска радиусом R′ и шириной dR′ имеет массу:
\(dM=\Sigma(R’)\,2\pi R’\,dR’\)В BeeTheory это кольцо генерирует поле гравитационной волны, распространяющееся во всех трех пространственных измерениях. В приближении монополя расстояние до точки трехмерного поля сферического радиуса r равно:
\(D(r,R’)=\sqrt{r^2+R’^2}\)Числовая форма темной плотности такова:
\(\rho_{\mathrm{dark}}(r)=K\sum_{i=1}^{N}\Sigma_0e^{-R’_i/R_d}\frac{(1+\alpha D_i)e^{-\alpha D_i}}{D_i^2}\,2\pi R’_i\Delta R’\) \(D_i=\sqrt{r^2+R_i’^2},\qquad R’_i=\left(i-\frac{1}{2}\right)\frac{R_{\mathrm{max}}}{N}\) \(N=60,\qquad R_{\mathrm{max}}=25\,\mathrm{kpc}\)2.2 Замкнутая темная масса и круговая скорость
\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\) \(M_{\mathrm{dark}}(<r)\approx\sum_{j=1}^{30}4\pi r_j^2\rho_{\mathrm{dark}}(r_j)\Delta r\) \(V_{\mathrm{dark}}(R)=\sqrt{\frac{GM_{\mathrm{dark}}(<R)}{R}}\) \(V_c(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{dark}}^2(R)}\)2.3 Асимптотическое поведение
\(\rho_{\mathrm{dark}}(r)\approx\frac{2\pi K\Sigma_0R_d^2}{r^2}\left(1+\alpha r+\frac{\alpha^2r^2}{2}\right)e^{-\alpha r}\)Для αr ≪ 1:
\(\rho_{\mathrm{dark}}(r)\xrightarrow{\alpha r\ll1}\frac{2\pi K\Sigma_0R_d^2}{r^2}\) \(M_{\mathrm{dark}}(<r)\propto r\qquad\Longrightarrow\qquad V_{\mathrm{dark}}\approx\mathrm{constant}\)3. Результаты моделирования — интерактивные диаграммы
Приведенная ниже симуляция сохраняет численную модель, ползунки, кривую вращения, профиль массы, профиль плотности и обновление χ² в реальном времени. Вставьте эту страницу в WordPress с включенной функцией выполнения скриптов.
χ²/доф: — | ℓ: — кпк | ρ(R⊙): — ГэВ/см³
| r (кпк) | Mbar (10¹⁰ M⊙) | Mdark (10¹⁰ M⊙) | Mtot (10¹⁰ M⊙) | DM/бар | ρdark (ГэВ/см³) |
|---|---|---|---|---|---|
| Загрузка… | |||||
4. Физическая интерпретация и универсальность
4.1 Длина когерентности
Внутри длины когерентности ядро Юкавы ведет себя почти как ньютоновское ядро 1/D². Темная плотность следует приблизительно r-², а кривая вращения плоская. За пределами ℓ экспоненциальное подавление приводит к спаду, наблюдаемому во внешнем диске.
\(\ell=\frac{1}{\alpha}\approx11.5\,\mathrm{kpc}\) \(\frac{\ell}{R_d}=\frac{11.5}{2.6}\approx4.4\)4.2 Безразмерная связь
Безразмерная связь BeeTheory может быть определена как:
\(\lambda_{\mathrm{galaxy}}=K\ell^2\) \(\lambda_{\mathrm{galaxy}}=0.040\times(11.5)^2\approx5.3\)Это сравнимо по порядку величины со связью, полученной в результате калибровки H₂, где λ составляет около 3-4. Вопрос о возможной универсальности этого числа остается открытым.
4.3 Сравнение со стандартными моделями
| Модель | Параметры | Типичная посадка | Весы | Механизм |
|---|---|---|---|---|
| Изотермический ореол | 2 | Умеренный | Радиус сердцевины | Феноменологическая плоская кривая |
| Профиль NFW | 2 | Сильный | rs | Профиль моделирования N-тела |
| Einasto | 2-3 | Сильный | r-2 | Гибкий эмпирический профиль |
| BeeTheory 3D Yukawa | 2 | Promising | ℓ | Сопряжение волновой массы с диском |
Крайняя точка эпохи Гайи остается самым сложным ограничением. Более резкий спад может быть получен при меньшей длине когерентности, но это ухудшает внутреннее соответствие. Будущие данные Gaia DR4, шаровых скоплений и звездных потоков станут важными проверками.
Ссылки
- Ou, X. et al. - Профиль темной материи Млечного Пути, полученный из его кривой круговой скорости, MNRAS 528, 2024.
- Дютертре, X. - Bee Theory™: Волновое моделирование гравитации, BeeTheory.com v2, 2023.
- Фримен, К. К. - О дисках спиральных и S0 галактик, ApJ 160, 811, 1970.
- Макмиллан, П. Дж. - Распределение массы и гравитационный потенциал Млечного Пути, MNRAS 465, 76, 2017.
- Наварро, Дж. Ф., Френк, К. С., Уайт, С. Д. М. - Универсальный профиль плотности на основе иерархической кластеризации, ApJ 490, 1997.
BeeTheory.com - Квантовая гравитация на основе волн
© Technoplane S.A.S. - 2025 г.