Математические основы галактической недостающей массы: Диск, сфера, плотность, потенциал и радиальное масштабирование
TL;DR: Проблема недостающей массы возникает, когда масса, вычисленная по кривым галактического вращения, превышает массу, непосредственно наблюдаемую в звездах, газе и пыли. С математической точки зрения, для этого необходимо связать поверхностную плотность на диске, объемную плотность в трех измерениях, гравитационный потенциал, радиальное ускорение и вложенную массу.
1. Радиальные координаты и геометрия
Мы различаем две геометрии:
- Геометрия диска: видимая галактическая материя в основном распределена в тонком вращающемся диске.
- Сферическая геометрия: темная или отсутствующая масса часто моделируется как примерно сферический ореол.
Элемент дисковой области:
\[ dA = R\,dR\,d\phi \]
Элемент сферического объема:
\[ dV = 4\pi r^2\,dr \]
Тот же символ \(r\) часто используется для галактоцентрического радиуса, но его значение зависит от геометрии. В диске \(R\) — это цилиндрический радиус. В гало \(r\) — это, как правило, сферический радиус.
2. Видимая масса на галактическом диске
Видимый диск часто аппроксимируется экспоненциальной поверхностной плотностью:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Масса кольца между \(R\) и \(R+dR\) равна:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Таким образом, суммарная видимая масса диска составляет:
\[ M_{\rm disk}(R)=2\pi\int_0^R R’\Sigma_0e^{-R’/R_d}\,dR’ \]
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
При большом радиусе:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Масса видимого диска приближается к конечному значению.
3. Сферическая масса и объемная плотность
Для сферического распределения массы объемная плотность \(\rho(r)\) определяет объемную массу:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Обратное отношение:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Это соотношение является центральным в проблеме недостающей массы. Если предполагаемая масса линейно растет с радиусом, то соответствующая сферическая плотность уменьшается как \(1/r^2\).
4. Динамическая масса от кругового движения
Для кругового движения гравитационное ускорение удовлетворяет:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Поэтому:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Для плоской кривой вращения:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Это дает стандартное масштабирование:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Определение недостающей массы
Недостающая масса — это разница между динамической и видимой массой:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Для экспоненциального видимого диска:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} — 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
Для \(v(r)\approx v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r — 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\frac{r}{R_d} \right) \right] \]
При большом радиусе дисковый член насыщается, а динамический член продолжает расти примерно как \(r\).
6. Гравитационный потенциал в 3D
Ньютоновский гравитационный потенциал, создаваемый точечной массой, равен:
\[ \Phi(r)=-\frac{GM}{r} \]
Соответствующее гравитационное поле является радиальной производной от потенциала:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Это объясняет взаимосвязь между \(1/r\) и \(1/r^2\): потенциал локализованной массы падает как \(1/r\), а сила или ускорение падают как \(1/r^2\).
7. Уравнение Пуассона
Плотность массы и гравитационный потенциал связаны через уравнение Пуассона:
\[ \nabla^2\Phi=4\pi G\rho \]
В сферической симметрии это выглядит так:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Это уравнение связывает три величины:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Потенциал расширенной 3D-плотности
Для общего трехмерного распределения плотности потенциал имеет вид:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Ядро \(1/|\mathbf{x}-\mathbf{x}’|\) — это математическое происхождение \(1/r\) потенциала в трех измерениях.
9. Потенциал тонкого диска
Для тонкого диска с поверхностной плотностью \(\Sigma(R’)\) гравитационный потенциал в плоскости диска может быть записан как:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’\,dR’\,d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Расстояние между точкой поля с радиусом \(R\) и точкой источника с радиусом \(R’\) равно:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Радиальное ускорение в диске можно получить, продифференцировав потенциал:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Скорость вращения зависит от:
\[ v^2(R)=R\,|g_R(R)| \]
10. Проекция 3D взаимодействия на диск
Если взаимодействие распространяется в трех измерениях, но оценивается в плоскости диска, радиальная проекция вносит геометрический фактор. Для двух точек диска, разделенных расстоянием \(d\), коэффициент радиальной проекции равен:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Таким образом, общее трехмерное радиальное ядро \(K(d)\), спроецированное на диск, выглядит так:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Например, ядро, подобное ньютоновской силе, имеет:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm диск}\propto \frac{R-R’\cos\phi}{d^3} \]
Экспоненциальное 3D-ядро можно записать как:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Экспоненциальные ядра в трех измерениях
Чисто экспоненциальный радиальный коэффициент имеет вид:
\[ e^{-r/\lambda} \]
В трехмерной теории поля экспоненциально экранированный потенциал часто появляется в форме, подобной Юкаве:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Соответствующая радиальная сила содержит \(1/r^2\) и экспоненциальные члены:
\[ g_Y(r) = -\frac{d}{dr} \left( \frac{e^{-r/\lambda}}{r} \right) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \right) \]
Это показывает, почему экспоненциальное радиальное поведение в 3D не зависит от \(1/r\) геометрии. Экспоненциальный показатель контролирует затухание, в то время как \(1/r\) и \(1/r^2\) возникают из-за трехмерного распространения.
12. Законы радиального масштабирования
Проблема недостающей массы сильно связана с радиальным масштабированием. Несколько важных радиальных законов проявляются неоднократно:
| Количество | Типичное масштабирование | Значение |
|---|---|---|
| Потенциал точечной массы | \(\Phi(r)\sim 1/r\) | 3D функция Грина гравитации |
| Сила, действующая на точечную массу | \(g(r)\sim 1/r^2\) | Производная от \(1/r\) |
| Скорость плоского вращения | \(v(r)\sim постоянная\) | Наблюдается во внешних галактических дисках |
| Динамическая масса | \(M(r)\sim r\) | Требуется для плоского вращения |
| Плотность ореола | \(\rho(r)\sim 1/r^2\) | Дает \(M(r)\sim r\) |
| Экспоненциальный диск | \(\Sigma(R)\sim e^{-R/R_d}\) | Видимый диск быстро исчезает |
| Экранированный 3D потенциал | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Экспоненциальное затухание плюс 3D-распространение |
13. От плотности к кривой вращения
Для сферического ореола с плотностью:
\[ \rho(r)\propto \frac{1}{r^2} \]
Приложенная масса составляет:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Затем:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approx постоянная \]
\[ v(r)\approx константа \]
Это математический мост между \(1/r^2\) плотностью гало и плоской кривой галактического вращения.
14. От массы диска к отсутствующей массе
Масса видимого диска сначала быстро растет, а затем насыщается:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Динамическая масса, выведенная на основании плоской кривой вращения, продолжает расти:
\[ M_{\rm dyn}(r)\propto r \]
Поэтому недостающая масса ведет себя примерно так:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
При достаточно большом радиусе:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Математическое предупреждение: диск и сфера не являются взаимозаменяемыми
Уравнение
\[ M(r)=\frac{v(r)^2r}{G} \]
является точным для сферической симметрии. Для сплюснутого диска необходимо вычислить потенциал путем интегрирования по диску, а затем вывести радиальное ускорение. Сферическое выражение часто используется в качестве эффективного приближения, особенно при обсуждении массы, необходимой для поддержания заданной кривой вращения.
16. Сводка ключевых уравнений
Плотность поверхности диска:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Элемент массы диска:
\[ dM_{\rm диск}=2\pi R\Sigma(R)dR \]
Видимая масса диска:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\frac{R}{R_d} \right) \right] \]
Сферическая объемная масса:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Плотность от заключенной массы:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Динамическая масса:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Недостающая масса:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Ньютоновский потенциал:
\[ \Phi(r)=-\frac{GM}{r} \]
Уравнение Пуассона:
\[ \nabla^2\Phi=4\pi G\rho \]
Трехмерный потенциальный интеграл:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Потенциал тонкого диска:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’dR’d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Проецируемое радиальное ядро:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Экспоненциальный трехмерный экранированный потенциал:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Заключение
Проблема недостающей массы — это математическое несоответствие между двумя радиальными моделями поведения. Видимый диск имеет экспоненциальную поверхностную плотность и достигает конечной суммарной массы. Динамическая масса, выведенная из приблизительно плоских кривых вращения, растет примерно линейно с радиусом. Если интерпретировать это как сферическое гало, то это соответствует плотности, уменьшающейся примерно как \(1/r^2\). Уравнения интегрирования диска, сферических оболочек, гравитационного потенциала и радиальной проекции предоставляют математический язык, необходимый для анализа этого несоответствия.