Численное моделирование модели темной массы BeeTheory

Полный, воспроизводимый отчет о численном интегрировании, разложении барионных скоростей, минимизации χ² и вариантах реализации, лежащих в основе моделирования темной массы в BeeTheory.

Эта техническая страница объясняет, как воспроизвести численное моделирование скрытой массы Млечного Пути в BeeTheory. Здесь описаны данные наблюдений, барионная модель, уравнение плотности на основе волн, численное интегрирование, метод подгонки, тесты на сходимость и справочный код.

Цель проста: начните с видимого диска Млечного Пути, примените модель волновой плотности BeeTheory, вычислите эффективную скрытую массу и сравните полученную кривую круговой скорости с наблюдениями эпохи Gaia.

0. Что мы вычисляем и почему

Кривая вращения Млечного Пути — это круговая скорость Vc(R) звезд в зависимости от их расстояния R от Галактического центра. Сегодня она измеряется гораздо точнее, чем общее распределение массы, которое мы можем непосредственно наблюдать.

Дефицит между наблюдаемой скоростью и тем, что предсказывает видимая барионная материя, — это проблема скрытой массы. Стандартные модели предполагают наличие невидимого гало частиц, обычно холодной темной материи. BeeTheory предлагает другую интерпретацию: каждый элемент видимой массы излучает волновое поле, которое экспоненциально затухает в трех измерениях, и накопленное поле ведет себя как скрытая масса.

Симуляция делает три вещи:

1. Вычислите Vcbar(R) из видимого диска и выпуклости, используя аналитическую формулу диска Фримена.
2. Численно проинтегрируйте плотность BeeTheory _ρdark(r; ℓ, λ) на каждом радиусе, а затем преобразуйте ее в ^VcDM(R).
3. Минимизирует χ²(ℓ, λ) по отношению к кривой вращения Млечного Пути эпохи Гайи, чтобы найти репрезентативные параметры для наилучшей подгонки.

Моделирование разработано таким образом, чтобы его можно было воспроизвести. Она может выполняться на JavaScript или Python без специализированной библиотеки по астрофизике.

Условные обозначения, используемые во всем тексте

• R — цилиндрический галактоцентрический радиус в плоскости диска, в кпк.
• r — сферический галактоцентрический радиус.
• z — это высота над диском.
• ℓ — длина когерентности BeeTheory, в кпк.
• λ — безразмерная связь между волной и массой.
• G используется в единицах кпк км² с-² M⊙-¹.

1. Данные наблюдений — Gaia 2024

Набор данных основан на кривой вращения Млечного Пути времен Гайи. В нем используются радиусы R, круговые скорости _Vc и неопределенности σ. В первоначальной технической версии использовалось 16 точек данных, охватывающих R = 4-27,3 кпк.

Важными фактами наблюдения являются:

• _Vc(R⊙ = 8 кпк) ≈ 230 км/с, точка привязки вблизи орбиты Солнца.
• _Vc приблизительно плоский в диапазоне от 5 до 20 кпк.
• _Vc уменьшается во внешнем измеряемом диске, достигая около 173 ± 17 км/с вблизи 27,3 кпк в данных, полученных по ссылке.

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

Мы исключили самую внутреннюю часть Галактики, потому что там преобладает центральный бар и некруговые движения. Упрощенная осесимметричная модель не является надежной внутри этой области.

3. Уравнение темной плотности BeeTheory

3.1 Физический постулат

Каждый видимый элемент массы dV в положении r′ в галактическом диске генерирует поле гравитационных волн, которое создает эффективную плотность массы в точке поля r:

Полная плотность темноты в любой точке (R,z) — это суперпозиция по всем элементам диска. Поскольку источник представляет собой диск, член объемной плотности становится членом поверхностной плотности:

Полный двойной интеграл:

3.3 Аналитическая проверка асимптотического предела

Для _Rd ≪ r ≪ ℓ выражение упрощается. Диск намного меньше радиуса, а длина когерентности все еще намного больше радиуса.

Потому что:

Эта асимптотическая плотность ведет себя как ρ ∝ r-², что дает M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}}\approx\mathrm{constant}\)

Подставив репрезентативные значения, Вы получите локальную плотность, близкую к наблюдаемому значению Млечного Пути:

5. χ² Минимизация и подгонка параметров

5.1 Целевая функция

Подгонка минимизирует уменьшенный хи-квадрат:

с N = 16 точками данных и p = 2 свободными параметрами, ℓ и λ. Это дает 14 степеней свободы.

5.2 Двухпроходной поиск по сетке

Вместо градиентного спуска используется сеточный поиск, поскольку ландшафт имеет длинную, изогнутую долину вырождения между λ и ℓ.

• Проход 1: грубая сетка на ℓ и λ.
• Проход 2: локальное уточнение вокруг грубого минимума.

Репрезентативная область наилучшего соответствия — это:

5.3 Долина вырождения

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

Внешний спад кривой вращения нарушает это вырождение, поскольку меньшие значения ℓ раньше подавляют эффективную плотность.

Данные за пределами 30 кпк, включая шаровые скопления, звездные потоки, звезды гало и галактики-спутники, необходимы для более точного определения ℓ.

7. Полный код ссылки

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

// Физические константы
const G = 4.302e-3; // кпк км² с-² M☉-¹
const Sig0 = 800.0; // M☉ pc-²
const Rd = 2.6; // кпк
const Mdisk = 3.5e10; // M☉
const Mbulge = 1.2e10; // M☉
const CONV_RHO = (1.989e30 * 5.61e26) / (3.086e21)**3;

// Данные о вращении в эпоху Гайи
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];

// Заполнитель барионной скорости
function vBaryonic(R) {
  // В полной реализации здесь используется формула диска Фримена
  // с модифицированными функциями Бесселя I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  return Math.sqrt(Math.max(0, vb2));
}

// Плотность темноты BeeTheory
function rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30.0 / N;
  let sum = 0;

  for (let i = 0; i < N; i++) {
    const Rp = (i + 0.5) * dRp;
    const Sig = Sig0 * Math.exp(-Rp / Rd);
    const K = (2 * Math.PI * ell / r)
      * Math.sinh(Math.min(r / ell, 30))
      * Math.exp(-Math.min((r + Rp) / ell, 30));

    sum += Sig * Rp * K * dRp;
  }

  return (lam / ell) * sum;
}

// Замкнутая темная масса
function enclosedDarkMass(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  пусть M = 0;

  for (let j = 0; j < Nr; j++) {
    const rj = (j + 0.5) * dr;
    M += 4 * Math.PI * rj * rj * rhoDark(rj, ell, lam) * dr;
  }

  return M;
}

// Темновая и полная скорость
function vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

function vTotal(R, ell, lam) {
  const vb = vBaryonic(R);
  const vd = vDM(R, ell, lam);
  return Math.sqrt(vb * vb + vd * vd);
}

// Хи-квадрат
function chiSquared(ell, lam) {
  пусть s = 0;

  for (let i = 0; i < OBS_R.length; i++) {
    const dv = (vTotal(OBS_R[i], ell, lam) - OBS_V[i]) / OBS_ERR[i];
    s += dv * dv;
  }

  return s / (OBS_R.length - 2);
}

Полная версия должна включать точную реализацию модифицированных функций Бесселя _I0, _I1, _K0 и _K1 для скорости диска Фримена.

Ссылки

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