Galaktik Kayıp Kütlenin Matematiksel Temelleri: Disk, Küre, Yoğunluk, Potansiyel ve Radyal Ölçeklendirme
TL;DR: Galaktik dönüş eğrilerinden çıkarılan kütle, yıldızlarda, gazda ve tozda doğrudan gözlemlenen kütleyi aştığında kayıp kütle sorunu ortaya çıkar. Matematiksel olarak bu, bir disk üzerindeki yüzey yoğunluğunu, üç boyuttaki hacim yoğunluğunu, yerçekimi potansiyelini, radyal ivmeyi ve kapalı kütleyi birbirine bağlamayı gerektirir.
1. Radyal koordinatlar ve geometri
İki geometri ayırt ediyoruz:
- Disk geometrisi: görünür galaktik madde esas olarak ince bir dönen disk içinde dağılmıştır.
- Küresel geometri: karanlık veya kayıp kütle genellikle kabaca küresel bir hale olarak modellenir.
Disk alanı elemanı:
\[ dA = R\,dR\,d\phi \]
Küresel hacim elemanı:
\[ dV = 4\pi r^2\,dr \]
Aynı sembol \(r\) genellikle galakosentrik yarıçap için kullanılır, ancak anlamı geometriye bağlıdır. Bir diskte \(R\) silindirik bir yarıçaptır. Bir halede, \(r\) genellikle küresel bir yarıçaptır.
2. Galaktik bir disk üzerinde görünür kütle
Görünür disk genellikle üstel bir yüzey yoğunluğu ile yaklaştırılır:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
(R\) ile \(R+dR\) arasındaki bir halkanın kütlesi:
\[ dM_{\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Bu nedenle kümülatif görünür disk kütlesi:
\[ 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 \sol[ 1-e^{-R/R_d} \sol( 1+\frac{R}{R_d} \sağ) \right] \]
Büyük yarıçapta:
\[ M_{\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Görünür disk kütlesi sonlu bir değere yaklaşır.
3. Küresel kütle ve hacim yoğunluğu
Küresel bir kütle dağılımı için, hacim yoğunluğu \(\rho(r)\) kapalı kütleyi belirler:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Ters ilişki şöyledir:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Bu ilişki kayıp kütle probleminin merkezinde yer alır. Eğer çıkarılan kütle yarıçapla doğrusal olarak büyüyorsa, buna karşılık gelen küresel yoğunluk \(1/r^2\) kadar azalır.
4. Dairesel hareketten dinamik kütle
Dairesel hareket için, yerçekimi ivmesi aşağıdakileri sağlar:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Bu yüzden:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Düz bir dönüş eğrisi için:
\[ v(r)\approx v_0 \]
\[ M_{\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Bu standart ölçeklendirmeyi verir:
\[ M_{\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Eksik kütle tanımı
Kayıp kütle, dinamik kütle ile görünür kütle arasındaki farktır:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Üstel görünür bir disk için:
\[ M_{\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \sol[ 1-e^{-r/R_d} \sol( 1+\frac{r}{R_d} \sağ) \right] \]
Çünkü \(v(r)\yaklaşık v_0\):
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \sol[ 1-e^{-r/R_d} \sol( 1+\frac{r}{R_d} \sağ) \right] \]
Büyük yarıçapta, disk terimi doygunluğa ulaşırken, dinamik terim yaklaşık olarak \(r\) kadar büyümeye devam eder.
6. 3B’de yerçekimi potansiyeli
Noktasal bir kütle tarafından üretilen Newton yerçekimi potansiyeli:
\[ \Phi(r)=-\frac{GM}{r} \]
Karşılık gelen yerçekimi alanı, potansiyelin radyal türevidir:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Bu \(1/r\) ve \(1/r^2\) arasındaki ilişkiyi açıklar: lokalize bir kütlenin potansiyeli \(1/r\) olarak düşerken, kuvvet veya ivme \(1/r^2\) olarak düşer.
7. Poisson denklemi
Kütle yoğunluğu ve yerçekimi potansiyeli Poisson denklemi aracılığıyla birbirine bağlıdır:
\[ \nabla^2\Phi=4\pi G\rho \]
Küresel simetride, bu şu hale gelir:
\[ \frac{1}{r^2}\frac{d}{dr} \sol( r^2\frac{d\Phi}{dr} \sağ) = 4\pi G\rho(r) \]
Bu denklem üç büyüklüğü birbirine bağlar:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Genişletilmiş 3D yoğunluk potansiyeli
Genel bir 3D yoğunluk dağılımı için potansiyel şudur:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
\(1/|\mathbf{x}-\mathbf{x}’|\) çekirdeği, üç boyuttaki \(1/r\) potansiyelinin matematiksel kökenidir.
9. İnce bir diskin potansiyeli
Yüzey yoğunluğu \(\Sigma(R’)\) olan ince bir disk için, disk düzlemindeki yerçekimi potansiyeli şu şekilde yazılabilir:
\[ \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}} \]
Yarıçap \(R\)’deki bir alan noktası ile yarıçap \(R’\)’deki bir kaynak noktası arasındaki mesafe:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Diskteki radyal ivme, potansiyelin farklılaştırılmasıyla elde edilir:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Dönüş hızı aşağıdakilerden kaynaklanır:
\[ v^2(R)=R\,|g_R(R)| \]
10. 3D etkileşimin disk üzerine yansıtılması
Bir etkileşim üç boyutta yayılıyor ancak disk düzleminde değerlendiriliyorsa, radyal izdüşüm geometrik bir faktör getirir. Diskte \(d\) mesafesiyle ayrılmış iki nokta için radyal projeksiyon faktörü şudur:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Böylece, disk üzerine yansıtılan genel bir 3B radyal çekirdek \(K(d)\) olarak görünür:
\[ K_{\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Örneğin, Newton kuvveti benzeri bir çekirdek vardır:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Üstel bir 3D çekirdek şu şekilde yazılabilir:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Üç boyutta üstel çekirdekler
Tamamen üstel bir radyal faktör formuna sahiptir:
\[ e^{-r/\lambda} \]
Üç boyutlu alan teorisinde, üstel olarak perdelenmiş bir potansiyel genellikle Yukawa benzeri bir biçimde ortaya çıkar:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
İlgili radyal kuvvet hem \(1/r^2\) hem de üstel terimler içerir:
\[ g_Y(r) = -\frac{d}{dr} \sol( \frac{e^{-r/\lambda}}{r} \sağ) \]
\[ g_Y(r)\propto e^{-r/\lambda} \sol( \frac{1}{r^2} + \frac{1}{\lambda r} \sağ) \]
Bu, 3B’deki üstel radyal davranışın neden \(1/r\) geometrisinden bağımsız olmadığını göstermektedir. Üstel zayıflamayı kontrol ederken \(1/r\) ve \(1/r^2\) üç boyutlu yayılmadan kaynaklanır.
12. Radyal ölçeklendirme yasaları
Kayıp kütle problemi radyal ölçeklendirmeye güçlü bir şekilde bağlıdır. Birkaç önemli radyal yasa tekrar tekrar ortaya çıkar:
| Miktar | Tipik ölçeklendirme | Anlamı |
|---|---|---|
| Noktasal kütlenin potansiyeli | \(\Phi(r)\sim 1/r\) | Yerçekiminin 3D Green fonksiyonu |
| Noktasal kütlenin kuvveti | \(g(r)\sim 1/r^2\) | ‘nin türevi (1/r\) |
| Düz dönüş hızı | \(v(r)\sim sabiti\) | Dış galaktik disklerde gözlemlendi |
| Dinamik kütle | \(M(r)\sim r\) | Düz rotasyon için gerekli |
| Halo yoğunluğu | \(\rho(r)\sim 1/r^2\) | \(M(r)\sim r\) verir |
| Üstel disk | \(\Sigma(R)\sim e^{-R/R_d}\) | Görünür disk hızla kaybolur |
| Elenmiş 3D potansiyel | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Üstel zayıflama artı 3D yayılma |
13. Yoğunluktan dönme eğrisine
Yoğunluğa sahip küresel bir halo için:
\[ \rho(r)\propto \frac{1}{r^2} \]
kapalı kütle:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Sonra:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\yaklaşık sabit \]
\[ v(r)\yaklaşık sabit \]
Bu, \(1/r^2\) halo yoğunluğu ile düz bir galaktik dönüş eğrisi arasındaki matematiksel köprüdür.
14. Disk kütlesinden kayıp kütleye
Görünür disk kütlesi ilk başta hızla büyür ve sonra doygunluğa ulaşır:
\[ M_{\rm disk}(R)\rightarrow M_d \]
Düz bir dönüş eğrisinden çıkarılan dinamik kütle büyümeye devam eder:
\[ M_{\rm dyn}(r)\propto r \]
Bu nedenle kayıp kütle yaklaşık olarak şu şekilde davranı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 \]
Yeterince büyük yarıçapta:
\[ M_{\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Matematiksel uyarı: disk ve küre birbirinin yerine kullanılamaz
Denklem
\[ M(r)=\frac{v(r)^2r}{G} \]
küresel simetri için kesindir. Düzleştirilmiş bir disk için, disk üzerinde integral alarak potansiyel hesaplanmalı ve ardından radyal ivme türetilmelidir. Küresel ifade, özellikle belirli bir dönme eğrisini desteklemek için gereken kütleyi tartışırken, genellikle etkili bir yaklaşım olarak kullanılır.
16. Anahtar denklemler özeti
Disk yüzey yoğunluğu:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Disk kütle elemanı:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Görünür disk kütlesi:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \sol[ 1-e^{-R/R_d} \sol( 1+\frac{R}{R_d} \sağ) \right] \]
Küresel hacim kütlesi:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Kapalı kütleden elde edilen yoğunluk:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Dinamik kütle:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Kayıp kütle:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Newton potansiyeli:
\[ \Phi(r)=-\frac{GM}{r} \]
Poisson denklemi:
\[ \nabla^2\Phi=4\pi G\rho \]
3D potansiyel integral:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
İnce disk potansiyeli:
\[ \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}} \]
Öngörülen radyal çekirdek:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Üstel 3D ekranlı potansiyel:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Sonuç
Kayıp kütle problemi iki radyal davranış arasındaki matematiksel bir uyumsuzluktur. Görünür disk üstel bir yüzey yoğunluğu izler ve sonlu bir kümülatif kütleye ulaşır. Yaklaşık olarak düz dönüş eğrilerinden çıkarılan dinamik kütle yarıçapla kabaca doğrusal olarak büyür. Küresel bir halo olarak yorumlanırsa, bu yaklaşık olarak \(1/r^2\) kadar azalan bir yoğunluğa karşılık gelir. Disk entegrasyonu, küresel kabuklar, yerçekimi potansiyeli ve radyal projeksiyon denklemleri bu uyumsuzluğu analiz etmek için gereken matematiksel dili sağlar.