Arı Teorisi Karanlık Kütle Modelinin Sayısal Simülasyonu

Sayısal entegrasyon, baryonik hız ayrıştırması, χ² minimizasyonu ve BeeTheory karanlık kütle simülasyonunun arkasındaki uygulama seçimlerinin eksiksiz, tekrarlanabilir bir açıklaması.

Bu teknik sayfa, Samanyolu’nun gizli kütlesinin BeeTheory sayısal simülasyonunun nasıl yeniden üretileceğini açıklamaktadır. Gözlemsel veriler, baryonik model, dalga tabanlı yoğunluk denklemi, sayısal entegrasyon, uydurma yöntemi, yakınsama testleri ve referans kodu açıklanmaktadır.

Amaç basit: görünür bir Samanyolu diskinden başlamak, BeeTheory dalga yoğunluğu modelini uygulamak, etkin gizli kütleyi hesaplamak ve sonuçta ortaya çıkan dairesel hız eğrisini Gaia dönemi gözlemleriyle karşılaştırmak.

0. Ne Hesaplıyoruz ve Neden Hesaplıyoruz

Samanyolu dönüş eğrisi, yıldızların Galaktik Merkez’e olan R uzaklıklarının bir fonksiyonu olarak dairesel hızları Vc(R)’dir. Günümüzde doğrudan görebildiğimiz toplam kütle dağılımından çok daha hassas bir şekilde ölçülmektedir.

Gözlenen hız ile görünür baryonik maddenin öngördüğü hız arasındaki açık, gizli kütle problemidir. Standart modeller görünmez bir parçacık halesine, genellikle soğuk karanlık maddeye başvurur. Arı Teorisi farklı bir yorum önermektedir: her görünür kütle elementi üç boyutta üstel olarak bozunan bir dalga alanı yayar ve biriken alan gizli kütle gibi davranır.

Simülasyon üç şey yapar:

1. Freeman’ın analitik disk formülünü kullanarak görünür disk ve şişkinlikten Vcbar(R) hesaplar.
2. ArıKuramı yoğunluğunu _ρdark(r; ℓ, λ) her yarıçapta sayısal olarak bütünleştirir, ardından ^VcDM(R)’ye dönüştürür.
3. Temsili en uygun parametreleri bulmak için Gaia dönemi Samanyolu dönüş eğrisine karşı χ²(ℓ, λ) değerini minimize eder.

Simülasyon tekrarlanabilir olacak şekilde tasarlanmıştır. Özel bir astrofizik kütüphanesi olmadan JavaScript veya Python’da çalışabilir.

Baştan sona kullanılan gösterim

• R, kpc cinsinden disk düzlemindeki silindirik galaktik merkez yarıçapıdır.
• r küresel galakosentrik yarıçaptır.
• z diskin üzerindeki yüksekliktir.
• ℓ, kpc cinsinden BeeTheory tutarlılık uzunluğudur.
• λ boyutsuz dalga-kütle bağıntısıdır.
• G, kpc km² s-² M⊙-¹ birimlerinde kullanılır.

1. Gözlemsel Veriler – Gaia 2024

Veri seti Gaia dönemi Samanyolu dönüş eğrisine dayanmaktadır. R yarıçapı, _Vc dairesel hızları ve σ belirsizliklerini kullanır. Orijinal teknik versiyonda R = 4-27.3 kpc’yi kapsayan 16 veri noktası kullanılmıştır.

Önemli gözlemsel gerçekler şunlardır:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, Güneş’in yörüngesine yakın bağlantı noktası.
• _Vc yaklaşık 5 ila 20 kpc arasında yaklaşık olarak düzdür.
• _Vc ölçülen dış diskte azalır ve referans alınan verilerde 27,3 kpc yakınında yaklaşık 173 ± 17 km/s’ye ulaşır.

Bu, orta yarıçaplarda güçlü bir şekilde yükselen, daha sonra daha büyük yarıçaplarda zayıflayan etkili bir gizli kütle dağılımı gerektirir.

En içteki Galaksiyi hariç tutuyoruz çünkü orada merkezi çubuk ve dairesel olmayan hareketler baskındır. Basitleştirilmiş bir eksenel simetrik model o bölge içinde güvenilir değildir.

3. Arı Teorisi Karanlık Yoğunluk Denklemi

3.1 Fiziksel Postülat

Galaktik diskteki r′ konumundaki her görünür kütle elemanı dV, r alan noktasında etkili bir kütle yoğunluğu yaratan bir yerçekimi dalga alanı oluşturur:

Herhangi bir noktadaki (R,z) toplam karanlık yoğunluk, tüm disk elemanlarının üst üste binmesidir. Kaynak bir disk olduğundan, hacim yoğunluğu terimi bir yüzey yoğunluğu terimi haline gelir:

Tam çift integral:

3.3 Asimptotik Limitin Analitik Olarak Doğrulanması

_Rd ≪ r ≪ ℓ için ifade basitleşir. Disk yarıçaptan çok daha küçüktür ve tutarlılık uzunluğu hala yarıçaptan çok daha büyüktür.

Çünkü:

Bu asimptotik yoğunluk ρ ∝ r-² olarak davranır, bu da M( \(\rho(r)\propto r^{-2}\Longrightarrow M(<r)\propto r\Longrightarrow V_c=\sqrt{\frac{GM(<r)}{r}}\approx\mathrm{constant}\)

Temsili değerlerin yerleştirilmesi, gözlemlenen Samanyolu değerine yakın bir yerel yoğunluk verir:

5. χ² Minimizasyonu ve Parametre Uydurma

5.1 Amaç Fonksiyonu

Uyum, indirgenmiş ki-kare değerini en aza indirir:

N = 16 veri noktası ve p = 2 serbest parametre, ℓ ve λ ile. Bu 14 serbestlik derecesi verir.

5.2 İki Geçişli Izgara Arama

Gradyan inişi yerine ızgara araması kullanılır çünkü manzara λ ve ℓ arasında uzun, kavisli bir dejenerasyon vadisine sahiptir.

• Geçiş 1: ℓ ve λ üzerinde kaba ızgara.
• Geçiş 2: kaba minimum etrafında yerel iyileştirme.

Temsili en uygun bölge şudur:

5.3 Yozlaşma Vadisi

χ² manzarası simetrik bir kase değildir. Uzatılmış bir vadi oluşturur, çünkü öncü yoğunluk normalizasyonu bağlanma gücüne güçlü bir şekilde bağlıdır ve düz dönme rejimi içindeki tutarlılık uzunluğuna sadece zayıf bir şekilde bağlıdır.

Dönme eğrisinin dış düşüşü bu dejenerasyonu kırar çünkü daha küçük ℓ değerleri etkin yoğunluğu daha erken bastırır.

Küresel kümeler, yıldız akıntıları, halo yıldızları ve uydu galaksiler de dahil olmak üzere 30 kpc’nin ötesindeki veriler ℓ’yi daha sıkı bir şekilde kısıtlamak için gereklidir.

7. Tam Referans Kodu

Aşağıdaki JavaScript referans uygulaması ana simülasyon mantığını yeniden üretir. Teknik doğrulama için tasarlanmıştır ve özel komut dosyalarına izin verilmediği sürece doğrudan standart bir WordPress içerik bloğunun içine değil, uygun bir komut dosyası ortamına yerleştirilmelidir.

// Fiziksel sabitler
const G = 4.302e-3; // kpc km² s-² M☉-¹
const Sig0 = 800.0; // M☉ pc-²
const Rd = 2.6; // kpc
const Mdisk = 3.5e10; // M☉
const Mbulge = 1.2e10; // M☉
const CONV_RHO = (1.989e30 * 5.61e26) / (3.086e21)**3;

// Gaia dönemi rotasyon verileri
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];

// Baryonik hız yer tutucusu
function vBaryonic(R) {
  // Tam uygulamada, bu Freeman'ın disk formülünü kullanır
  // değiştirilmiş Bessel fonksiyonları I0, I1, K0, K1 ile.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  Math.sqrt(Math.max(0, vb2)) döndür;
}

// BeeTheory karanlık yoğunluk
function rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30.0 / N;
  toplam = 0 olsun;

  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;
  }

  (lam / ell) * toplam döndürür;
}

// Kapalı karanlık kütle
function enclosedDarkMass(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  M = 0 olsun;

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

  M'yi döndür;
}

// Karanlık ve toplam hız
function vDM(R, ell, lam) {
  return Math.sqrt(Math.max(0, G * enclosedDarkMass(R, ell, lam) / R));
}

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

// Chi-squared
function chiSquared(ell, lam) {
  s = 0 olsun;

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

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

Tam bir sürüm, Freeman disk hızı için değiştirilmiş Bessel fonksiyonları _I0, _I1, _K0 ve _K1 'in doğru uygulamalarını içermelidir.

Referanslar

• Abramowitz, M., Stegun, I. A. - Handbook of Mathematical Functions, Dover, 1972.
• Bovy, J., Rix, H.-W. - Samanyolu'nun Disk Yüzey Yoğunluk Profilinin Doğrudan Dinamik Ölçümü, ApJ 779, 115, 2013.
• Freeman, K. C. - Spiral ve S0 galaksilerinin diskleri üzerine, ApJ 160, 811, 1970.
• McMillan, P. J. - The mass distribution and gravitational potential of the Samanyolu, MNRAS 465, 76, 2017.
• Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - Samanyolu'nun dairesel hız eğrisinden çıkarılan karanlık madde profili, MNRAS 528, 693, 2024.
• Pato, M., Iocco, F., Bertone, G. - Samanyolu'ndaki karanlık madde dağılımı üzerine dinamik kısıtlamalar, JCAP 12, 001, 2015.
• Portail, M. ve diğerleri - Galaktik şişkinlik ve çubuğun dinamik modellemesi, MNRAS 465, 1621, 2017.