Simulasi Numerik Model Massa Gelap BeeTheory

Penjelasan lengkap dan dapat direproduksi tentang integrasi numerik, dekomposisi kecepatan baryonik, minimalisasi χ², dan pilihan implementasi di balik simulasi massa gelap BeeTheory.

Halaman teknis ini menjelaskan cara mereproduksi simulasi numerik BeeTheory untuk massa tersembunyi Bima Sakti. Di dalamnya dijelaskan mengenai data observasi, model baryonik, persamaan kerapatan berbasis gelombang, integrasi numerik, metode fitting, uji konvergensi, dan kode referensi.

Tujuannya sederhana: mulai dari piringan Bima Sakti yang terlihat, terapkan model kerapatan gelombang BeeTheory, hitung massa tersembunyi yang efektif, dan bandingkan kurva kecepatan melingkar yang dihasilkan dengan pengamatan era Gaia.

0. Apa yang Kami Komputasi dan Mengapa

Kurva rotasi B imasakti adalah kecepatan edar Vc(R) bintang-bintang sebagai fungsi dari jarak R dari Pusat Galaksi. Saat ini, kurva ini diukur jauh lebih akurat daripada distribusi massa total yang bisa kita lihat secara langsung.

Perbedaan antara kecepatan yang teramati dan yang diprediksi oleh materi baryonik yang tampak adalah masalah massa yang tersembunyi. Model standar menyebutkan adanya lingkaran cahaya partikel yang tidak terlihat, biasanya berupa materi gelap yang dingin. BeeTheory mengusulkan interpretasi yang berbeda: setiap elemen massa yang terlihat memancarkan medan gelombang yang meluruh secara eksponensial dalam tiga dimensi, dan medan yang terakumulasi berperilaku seperti massa yang tersembunyi.

Simulasi ini melakukan tiga hal:

1. Menghitung Vcbar (R) dari cakram dan tonjolan yang terlihat dengan menggunakan rumus cakram analitik Freeman.
2. Secara numerik mengintegrasikan kerapatan BeeTheory _ρdark(r; ℓ, λ) pada setiap radius, kemudian mengubahnya menjadi ^VcDM(R).
3. Meminimalkan χ²(ℓ, λ) terhadap kurva rotasi Bima Sakti era Gaia untuk menemukan parameter-parameter yang paling sesuai.

Simulasi ini dirancang untuk dapat direproduksi. Simulasi ini dapat dijalankan dalam JavaScript atau Python tanpa pustaka astrofisika khusus.

Notasi yang digunakan di seluruh bagian

• R adalah jari-jari galaksi silindris pada bidang piringan, dalam kpc.
• r adalah radius galaksi sferis.
• z adalah ketinggian di atas disk.
• ℓ adalah panjang koherensi BeeTheory, dalam kpc.
• λ adalah kopling gelombang-massa yang tidak berdimensi.
• G digunakan dalam satuan kpc km² s-² M⊙-¹.

1. Data Pengamatan – Gaia 2024

Dataset ini didasarkan pada kurva rotasi Bima Sakti era Gaia. Kurva ini menggunakan jari-jari R, kecepatan edar _Vc, dan ketidakpastian σ. Versi teknis aslinya menggunakan 16 titik data yang mencakup R = 4-27,3 kpc.

Fakta-fakta pengamatan yang penting adalah:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/detik, titik jangkar di dekat orbit Matahari.
• _Vc kira-kira datar dari sekitar 5 hingga 20 kpc.
• _Vc menurun pada piringan terluar yang terukur, mencapai sekitar 173 ± 17 km/detik di dekat 27,3 kpc pada data yang dirujuk.

Hal ini membutuhkan distribusi massa tersembunyi yang efektif yang naik dengan kuat pada jari-jari menengah, kemudian melemah pada jari-jari yang lebih besar.

Kami tidak menyertakan Galaksi terdalam karena batang pusat dan gerakan non-lingkaran mendominasi di sana. Model axisymmetric yang disederhanakan tidak dapat diandalkan di dalam wilayah tersebut.

3. Persamaan Kepadatan Gelap Teori Lebah

3.1 Postulat Fisik

Setiap elemen massa yang terlihat dV pada posisi r′ di piringan galaksi menghasilkan medan gelombang gravitasi yang menciptakan kerapatan massa efektif pada titik medan r:

Kepadatan gelap total pada titik mana pun (R, z) adalah superposisi dari semua elemen cakram. Karena sumbernya adalah cakram, istilah densitas volume menjadi istilah densitas permukaan:

Integral ganda penuh adalah:

3.3 Verifikasi Analitis dari Batas Asimtotik

Untuk_Rd ≪ r ≪ ℓ, ekspresinya menjadi lebih sederhana. Piringan jauh lebih kecil daripada jari-jari, dan panjang koherensi masih jauh lebih besar daripada jari-jari.

Karena:

Densitas asimtotik ini berperilaku sebagai ρ ∝ r-², yang memberikan M( [lateks]\rho(r)\propto r^{-2}\Panah kanan atas M(<r)\propto r\Panah kanan bawah V_c=\sqrt{\frac{GM(<r)}{r}}\pendek\mathrm{konstanta}[/lateks]

Memasukkan nilai yang representatif memberikan kerapatan lokal yang mendekati nilai Bimasakti yang diamati:

5. Minimisasi χ² dan Penyesuaian Parameter

5.1 Fungsi Tujuan

Kecocokan ini meminimalkan chi-kuadrat yang berkurang:

dengan N = 16 titik data dan p = 2 parameter bebas, ℓ dan λ. Hal ini memberikan 14 derajat kebebasan.

5.2 Pencarian Kisi Dua Lintasan

Pencarian grid digunakan daripada gradient descent karena lanskap memiliki lembah degenerasi yang panjang dan melengkung antara λ dan ℓ.

• Lulus 1: kisi-kisi kasar pada ℓ dan λ.
• Lulus2: penyempurnaan lokal di sekitar minimum kasar.

Perwakilan wilayah yang paling sesuai adalah:

5.3 Lembah Kemunduran

Lanskap χ² bukanlah mangkuk simetris. Ini membentuk lembah yang memanjang karena normalisasi densitas utama sangat bergantung pada kekuatan kopling dan hanya lemah pada panjang koherensi di dalam rezim rotasi datar.

Penurunan luar kurva rotasi mematahkan kemunduran ini karena nilai ℓ yang lebih kecil menekan densitas efektif lebih awal.

Data di atas 30 kpc, termasuk gugus bola, aliran bintang, bintang halo, dan galaksi satelit, sangat penting untuk membatasi ℓ dengan lebih ketat.

7. Kode Referensi Lengkap

Implementasi referensi JavaScript berikut ini mereproduksi logika simulasi utama. Ini dimaksudkan untuk validasi teknis dan harus ditempatkan di lingkungan skrip yang tepat, tidak langsung di dalam blok konten WordPress standar kecuali jika skrip kustom diizinkan.

// Konstanta fisik
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;

// Data rotasi era Gaia
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,6,6,7,7,8,9,10,11,13,17];

// Penampung kecepatan baryonik
function vBaryonic(R) {
  // Dalam implementasi penuh, ini menggunakan rumus disk Freeman
  // dengan fungsi Bessel yang dimodifikasi I0, I1, K0, K1.
  const vb2 = G * Mbulge / Math.max(R, 0.2);
  return Math.sqrt(Math.max(0, vb2));
}

// Kepadatan gelap BeeTheory
function rhoDark(r, ell, lam) {
  const N = 40;
  const dRp = 30.0 / N;
  biarkan jumlah = 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));

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

  return (lam / ell) * sum;
}

// Massa gelap tertutup
function tertutupMassaGelap(R, ell, lam) {
  const Nr = 30;
  const dr = R / Nr;
  biarkan 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;
  }

  mengembalikan M;
}

// Kecepatan gelap dan total
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);
}

// Chi-kuadrat
function chiKuadrat(ell, lam) {
  let 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);
}

Versi yang lengkap harus menyertakan implementasi yang akurat dari fungsi Bessel yang dimodifikasi _I0, _I1, _K0, dan _K1 untuk kecepatan disk Freeman.

Referensi

• Abramowitz, M., Stegun, I. A. - Buku Pegangan Fungsi Matematika, Dover, 1972.
• Bovy, J., Rix, H.-W. - Pengukuran Dinamis Langsung Profil Kerapatan Permukaan Piringan Bima Sakti, ApJ 779, 115, 2013.
• Freeman, K. C. - Pada piringan galaksi spiral dan galaksi S0, ApJ 160, 811, 1970.
• McMillan, P. J. - Distribusi massa dan potensi gravitasi Bimasakti, MNRAS 465, 76, 2017.
• Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - Profil materi gelap Bima Sakti yang disimpulkan dari kurva kecepatan melingkarnya, MNRAS 528, 693, 2024.
• Pato, M., Iocco, F., Bertone, G. - Kendala dinamik pada distribusi materi gelap di Bima Sakti, JCAP 12, 001, 2015.
• Portail, M. dkk. - Pemodelan dinamis tonjolan dan batang Galaksi, MNRAS 465, 1621, 2017.