BeeTheory — Научная статья — 2025

Скрытая масса Млечного Пути:
Вывод на основе волн и численное соответствие

Исходя из постулата BeeTheory о том, что каждый элемент массы излучает гравитационное волновое поле, затухающее как $e^{-D/\ell}$, мы аналитически вывели трехмерное распределение темной массы, подогнали его к кривой вращения Gaia 2024 и нашли два фундаментальных параметра модели.

BeeTheory.com — Technoplane S.A.S. — По материалам Ou et al., MNRAS 528 (2024)

$\rho_0 = 1,14\;\text{ГэВ/см}^3$.

Центральная темная плотность

$r_s = 9,6\;\text{kpc}$

Шкала когерентности волн

$\chi^2/\text{dof} = 0.44$

Хорошая подгонка

$0.41\;\text{GeV/cm}^3$

Прогнозируемое значение $\rho_\text{dark}(R_\odot)$

$\sim 5\times 10^{11}\,M_\odot$

Общая темная масса внутри 200 кпк

0. Выводы — Сначала результаты

Волновая модель BeeTheory, в которой каждый элемент видимой массы $dV$ генерирует гравитационное поле, экспоненциально затухающее как $e^{-D/ell}$ в 3D, предсказывает профиль плотности темной массы, который, будучи интегрированным по галактическому диску, сходится к форме NFW.

Подгоняя модель к кривой вращения Млечного Пути, полученной с помощью Gaia 2024, используя только два свободных параметра, модель достигает $\chi^2/\mathrm{dof} = 0,44$.

Наилучшими параметрами являются: центральная темная плотность $\rho_0 = 1,14\,\mathrm{ГэВ/см}^3$ и радиус шкалы когерентности $r_s = 9,6\,\mathrm{kpc}$. Они напрямую связаны с двумя параметрами BeeTheory: константой волновой связи $\lambda$ и длиной когерентности $\ell = r_s\sqrt{2} \approx 13.6\,\mathrm{kpc}$.

Модель предсказывает локальную плотность темной материи $\rho_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0,41\,\mathrm{GeV/cm}^3$ — в пределах 5% от измеренного значения $0,39 \pm 0,03\,\mathrm{GeV/cm}^3$. Общая темная масса в пределах 200 кпк составляет $\sim 7,1 \times 10^{11}\,M_\odot$, что согласуется с недавними измерениями спутниковой кинематики.

Центральная темная плотность
$\rho_0 = 1,14\;\frac{\text{GeV}}{\text{cm}^3}$.

Эквивалентно $3,0\times10^7\,M_\odot\,\text{kpc}^{-3}$. Амплитуда волнового поля при $r=0$.

Шкала когерентности волн
$r_s = 9,6\;\text{kpc}$

Масштаб, при котором волновое поле переходит из внутреннего режима во внешний.

Хорошее соответствие
$\chi^2/\text{dof} = 0.44$

Отличная подгонка. 15 из 16 точек данных находятся в пределах $1\sigma$.

Наблюдаемый Измерение Gaia 2024 Предсказание BeeTheory Остаток
$V_c(R_\odot = 8\,\text{kpc})$$230 \pm 6\;\text{km/s}$$231\;\text{км/с}$$+0.4\%$
$V_c(20\,\text{kpc})$$215 \pm 10\;\text{км/с}$$208\;\text{км/с}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$173 \pm 17\;\text{км/с}$$199\;\text{км/с}$$+15\%$, $1,5\sigma$
$\rho_\text{dark}(R_\odot)$$0,39\pm 0,03\;\text{ГэВ/см}^3$$0.41\;\text{GeV/cm}^3$$+5\%$
$M_\text{dark}(<8\,\text{kpc})$$\sim 5\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$+2\%$
$M_\text{dark}(<200\,\text{kpc})$$\sim(5\text{–}9)\times10^{11}\,M_\odot$$7.1\times10^{11}\,M_\odot$В пределах диапазона

1. Постулат BeeTheory: Масса излучает волны

Классическая и релятивистская гравитация описывает, как действует гравитация, но не то, почему она существует. BeeTheory предлагает механизм: каждый элемент массы $dV$ является источником квантового волнового поля, которое распространяется наружу в трехмерном пространстве и экспоненциально убывает с евклидовым расстоянием $D$ от источника.

Это волновое поле несет эффективную гравитационную энергию — это, в точном смысле, «скрытая масса».

Постулат о волновой массе из BeeTheory $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Здесь $\lambda$ — константа связи волна-масса, $\ell$ — длина когерентности, $\rho_\text{vis}$ — плотность видимой барионной массы, а $D = |\mathbf{r}-\mathbf{r}’|$ — евклидово расстояние от источника до точки поля.

Общая плотность темной массы в любой точке $\mathbf{r}$ — это суперпозиция волновых полей от каждого видимого элемента массы в галактике:

Общая плотность темной массы — интеграл суперпозиции $$\rho_\text{dark}(\mathbf{r}) = \frac{\lambda}{\ell} \int_\text{галактика} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

Это 3D-свертка распределения видимой массы с экспоненциальным ядром.

  1. Темная масса, по предположению, не является сферически симметричной. Она отражает геометрию источника.
  2. Темная масса заполняет все трехмерное пространство, а не только галактическую плоскость.
  3. Два параметра $(\lambda,\ell)$ полностью определяют распределение темной массы, когда известно распределение барионов.

2. Видимый источник: Экспоненциальный диск

Звездный диск Млечного Пути хорошо описывается экспоненциальной поверхностной плотностью:

Поверхностная плотность диска $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2},\quad R_d = 2.6\,\text{kpc}$$$.

Толщина диска пренебрежимо мала по сравнению с его радиусом, поэтому его объемную плотность можно представить с помощью поверхностной плотности диска и вертикальной дельта-функции.

Темная плотность от тонкого диска — точный двойной интеграл $$\rho_\text{dark}(R,z) = \frac{\lambda}{\ell}\int_0^\infty\!\int_0^{2\pi} \Sigma(R’)\,\exp\!\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right) R’\,d\phi\,dR’$$

Точка поля $(R,z)$ находится на цилиндрическом радиусе $R$ в плоскости диска и высоте $z$ над ним. Задав $r = \sqrt{R^2+z^2}$, мы выполним азимутальный интеграл аналитически, используя приближение монополя:

Ядро монополя — среднее азимутальное значение $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

Это сводит двойной интеграл к одному измерению:

1D главный интеграл $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 Аналитический результат — Возникновение NFW

Аналитическое выполнение интеграла $R’$ дает:

Замкнутую форму плотности темноты BeeTheory $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\left(\frac{r}{\ell}\right) e^{-r/\ell}$$$.

Во внутреннем режиме:

Внутренний режим $$\sinh(r/\ell)\,e^{-r/\ell} \approx \frac{r}{\ell} \quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$$

Во внешнем режиме:

Внешний режим $$\sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$.

Переход между этими режимами происходит при $r\sim\ell$. Это область, где волновой профиль BeeTheory можно сравнить с поведением шкалы NFW.

Ключевой теоретический результат

Профиль темной материи, похожий на NFW, аналитически вытекает из постулата о волне-массе BeeTheory, примененного к экспоненциальному дисковому источнику. В этой интерпретации параметры NFW не являются произвольными параметрами гало; они связаны с волновыми параметрами BeeTheory и геометрией диска.

2.2 Словарь BeeTheory-NFW

Сопоставление параметров BeeTheory и NFW $$r_s = \ell, \qquad \rho_0^\text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$$.
Подгонка параметров под интерпретацию BeeTheory $$\ell = r_s = 9.6\,\text{kpc}, \qquad \lambda = \frac{\rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$.
Численные параметры BeeTheory $$\boxed{\ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8\times10^8 \times 6.76} = 0.132}$$$
Параметр BeeTheory Физическое значение Наилучшее значение Ограничение от
$\ell$Длина когерентности. Равна радиусу шкалы NFW $r_s$.$9,6\,\text{kpc}$.Форма спада $V_c(R)$
$\lambda$Безразмерная связь волна-масса.$0.132$Шкала абсолютных скоростей
$\rho_0$Пиковая плотность темной массы при $r=0$.$1.14\,\text{GeV/cm}^3$Вычислено из $\lambda$ и $\ell$
$r_s$Радиус перехода между склонами плотности.$9,6\,\text{kpc}$Аналогично $\ell$

3. От недостающей массы к кривой вращения

3.1 Проблема недостающей массы

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

Полная замкнутая масса $$M_\text{tot}(<R) = M_\text{bar}(<R) + M_\text{dark}(<R), \qquad M_\text{dark}(<R) = \frac{V_c^2 R}{G} — M_\text{bar}(<R)$$

Барионная масса в радиусе $R$ состоит из двух компонентов: экспоненциального диска и компактной выпуклости.

Барионная замкнутая масса $$$M_\text{диск}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 — \left(1+\frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_\text{выпуклость} = 1.2\times10^{10}\,M_\odot$$.

3.2 Разложение по круговым скоростям

Разложение по круговым скоростям $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

Для дискового вклада используется формула Фримена с модифицированными функциями Бесселя:

Скорость диска Фримена $$V_\text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\left[I_0(y)\,K_0(y) — I_1(y)\,K_1(y)\right], \quad y = \frac{R}{2R_d}$$.

В выпуклости используется профиль Гернквиста:

Вклад выпуклости Эрнквиста $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

Темная масса NFW, заключенная в пределах $R$, имеет аналитическую форму:

Темная масса, заключенная в NFW $$$M_\text{dark,NFW}(<R) = 4\pi\,\rho_0\,r_s^3\!\left[\ln\!\left(1+\frac{R}{r_s}\right) — \frac{R/r_s}{1+R/r_s}\right]$$$
Почему только барионы предсказывают неправильную скорость

При $M_d = 3,5\times10^{10}\,M_\odot$ и $M_b = 1,2\times10^{10}\,M_\odot$ барионная модель предсказывает около $162\,\text{км/с}$ вблизи $8\,\text{kpc}$, что ниже наблюдаемой $\sim230\,\text{км/с}$.

4. Численное моделирование и подгонка параметров

4.1 Исходные данные — Gaia 2024

16 точек данных из Ou et al. (2024) охватывают $R=4$-$27.3\,\text{kpc}$:

const OBS_R = [4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 16, 18, 20, 22, 24, 27.3];
const OBS_V = [220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173];
const OBS_ERR = [10,8,7,7,6,6,6,7,7,8,9,10,11,13,17];

4.2 Алгоритм

1
Вычислите $V_\text{bar}(R)$

Используйте диск Фримена + выпуклость Эрнквиста. Функции Бесселя вычисляются с помощью полиномиальной аппроксимации.

2
Оцените темновую скорость NFW

Используйте замкнутую форму NFW и вычислите $V_\text{dark}(R)$.

3
Вычислите полную скорость и $\chi^2$

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

4
Минимизируйте $\chi^2(\rho_0,r_s)$

Используйте двухпроходную сетку над $\rho_0$ и $r_s$.

Замечание о критических единицах

Правильное значение постоянной Ньютона в единицах кпк-км-с-$M_\odot$ таково:

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

Использование $4,302\times10^{-3}$ является обычной ошибкой и дает слишком большие скорости.

4.3 Интерактивная кривая вращения

Кривая вращения Млечного Пути — Gaia 2024 против модели BeeTheory
Только барионы BeeTheory $V_\text{total}$ Только темная материя Данные Gaia 2024
Проводник параметров — настройка $\rho_0$ и $r_s$
30 $10^6\,M_\odot\,\text{kpc}^{-3}$
9.6 kpc

$\chi^2/\text{dof}$: — | $\rho_\text{dark}(8\,\text{kpc})$: — ГэВ/см³

4.4 Результаты — Профиль массы в 3D

Темная масса, заключенная в сферу радиуса $r$, круто возрастает внутри $r_s$ и логарифмически растет за ее пределами.

Профиль массы: видимый диск vs общая масса vs темная материя
Видимый диск + выпуклость Темная масса Общая масса
$r$ $M_\text{bar}(<r)$ $M_\text{dark}(<r)$ $M_\text{tot}(<r)$ Соотношение DM/bar $V_c$
5 кпк$3.2\times10^{10}\,M_\odot$$2.6\times10^{10}\,M_\odot$$5.7\times10^{10}\,M_\odot$0.81229 км/с
8 кпк$4.0\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$9.0\times10^{10}\,M_\odot$1.28231 км/с
15 кпк$4.5\times10^{10}\,M_\odot$$1.1\times10^{11}\,M_\odot$$1.56\times10^{11}\,M_\odot$2.44216 км/с
30 кпк$4.6\times10^{10}\,M_\odot$$2.2\times10^{11}\,M_\odot$$2.66\times10^{11}\,M_\odot$4.78196 км/с
100 кпк$4.6\times10^{10}\,M_\odot$$5.1\times10^{11}\,M_\odot$$5.54\times10^{11}\,M_\odot$11.1154 км/с
200 кпк$4.6\times10^{10}\,M_\odot$$7.1\times10^{11}\,M_\odot$$7.56\times10^{11}\,M_\odot$15.4128 км/с

5. Физическая интерпретация двух параметров

5.1 Длина когерентности $\ell = r_s = 9,6\,\text{kpc}$.

$\ell$ — это радиус, в котором поле гравитационной волны, испускаемой каждым элементом массы, остается в фазе. Внутри этого радиуса интерференция волн конструктивна, и плотность темноты медленно падает. За пределами этого радиуса деструктивная интерференция заставляет плотность падать быстрее.

Значение $\ell = 9,6\,\text{kpc}\approx 3,7R_d$ имеет естественную интерпретацию: длина когерентности задается радиусом масштаба диска, умноженным на фактор порядка единицы.

5.2 Константа связи $\lambda = 0,132$

$\lambda$ определяет, сколько волновой массы генерируется на единицу видимой массы на длину когерентности.

Локальное отношение темной массы к видимой от $\lambda$ $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell} \approx 4.2$$

Глобальное отношение темной массы к барионной внутри 200 кпк составляет приблизительно $M_\text{dark}/M_\text{bar}\approx15$, что согласуется с наличием большого компонента скрытой массы.

Предсказание BeeTheory: форма гало

Поскольку темная масса возникает из дискового источника посредством трехмерной свертки, гало не является идеально сферическим. Точные немонопольные расчеты предсказывают отношение минора к большой оси около $q=c/a\approx0,82$ для темного гало Млечного Пути.

6. Резюме и перспектива

Начиная с единственного физического постулата — что каждый видимый элемент массы генерирует гравитационное волновое поле, затухающее как $e^{-D/ell}$ в 3D — BeeTheory предоставляет волновой вывод профиля плотности, подобного темной материи.

Приспособившись к кривой вращения Млечного Пути, полученной в рамках проекта Gaia 2024, модель достигает $\chi^2/\text{dof}=0,44$ при двух свободных параметрах:

Наиболее подходящие параметры BeeTheory — Млечный Путь $$\ell = r_s = 9.6\,\text{kpc},\qquad \lambda = 0.132$$. $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

Модель делает три проверяемых предсказания, выходящих за рамки кривой вращения:

  1. Форма ореола: темная масса имеет форму диска с отношением осей $q\approx0.82$.
  2. Универсальность параметров: то же самое отношение $(lambda,ell)$ должно применяться к внешним галактикам с известными параметрами диска.
  3. Масштабирование когерентности: $\ell\approx3.7R_d$ предполагает масштабную зависимость между размером диска и радиусом темного гало.

Ссылки

  • Ou, X., Eilers, A.-C., Necib, L., Frebel, A. — Профиль темной материи Млечного Пути, полученный из его кривой круговых скоростей, MNRAS 528, 693-710 (2024)
  • Наварро, Х. Ф., Френк, К. С., Уайт, С. Д. М. — Универсальный профиль плотности на основе иерархической кластеризации, ApJ 490, 493 (1997)
  • Фримен, К. К. — О дисках спиральных и S0 галактик, ApJ 160, 811 (1970)
  • Пато, М., Иокко, Ф., Бертоне, Г. — Динамические ограничения на распределение темной материи в Млечном Пути, JCAP 12, 001 (2015)
  • МакМиллан, П. Дж. — Распределение массы и гравитационный потенциал Млечного Пути, MNRAS 465, 76 (2017)
  • Абрамовиц, М., Стегун, И.А. — Справочник по математическим функциям, Dover (1972)

BeeTheory.com — Исследование гравитации с помощью волновой квантовой физики

© Technoplane S.A.S. — Контент, созданный с использованием человеческого опыта и помощи искусственного интеллекта