Αριθμητική προσομοίωση του μοντέλου σκοτεινής μάζας της θεωρίας BeeTheory

Μια πλήρης, αναπαραγώγιμη περιγραφή της αριθμητικής ολοκλήρωσης, της αποσύνθεσης της βαρυονικής ταχύτητας, της ελαχιστοποίησης του χ² και των επιλογών εφαρμογής της προσομοίωσης σκοτεινής μάζας BeeTheory.

Αυτή η τεχνική σελίδα εξηγεί πώς να αναπαράγετε την αριθμητική προσομοίωση της κρυμμένης μάζας του Γαλαξία μας από την BeeTheory. Περιγράφει τα δεδομένα παρατήρησης, το βαρυονικό μοντέλο, την εξίσωση πυκνότητας με βάση τα κύματα, την αριθμητική ολοκλήρωση, τη μέθοδο προσαρμογής, τις δοκιμές σύγκλισης και τον κώδικα αναφοράς.

Ο στόχος είναι απλός: ξεκινήστε από έναν ορατό δίσκο του Γαλαξία μας, εφαρμόστε το μοντέλο κυματικής πυκνότητας BeeTheory, υπολογίστε την αποτελεσματική κρυμμένη μάζα και συγκρίνετε την προκύπτουσα καμπύλη κυκλικής ταχύτητας με τις παρατηρήσεις της εποχής Gaia.

0. Τι υπολογίζουμε και γιατί

Η καμπύλη περιστροφής του Γαλαξία μας είναι η κυκλική ταχύτητα Vc(R) των αστέρων σε συνάρτηση με την απόστασή τους R από το Γαλαξιακό Κέντρο. Μετριέται σήμερα με πολύ μεγαλύτερη ακρίβεια από ό,τι η συνολική κατανομή της μάζας που μπορούμε να δούμε άμεσα.

Το έλλειμμα μεταξύ της παρατηρούμενης ταχύτητας και της ορατής βαρυονικής ύλης είναι το πρόβλημα της κρυμμένης μάζας. Τα καθιερωμένα μοντέλα επικαλούνται μια αόρατη φωτοστέφανο σωματιδίων, συνήθως ψυχρή σκοτεινή ύλη. Η θεωρία BeeTheory προτείνει μια διαφορετική ερμηνεία: κάθε ορατό στοιχείο μάζας εκπέμπει ένα κυματικό πεδίο που διασπάται εκθετικά στις τρεις διαστάσεις, και το συσσωρευμένο πεδίο συμπεριφέρεται σαν κρυμμένη μάζα.

Η προσομοίωση κάνει τρία πράγματα:

1. Υπολογίζει το Vcbar(R) από τον ορατό δίσκο και τη διόγκωση χρησιμοποιώντας τον αναλυτικό τύπο δίσκου του Freeman.
2. Ολοκληρώνει αριθμητικά την πυκνότητα BeeTheory _ρdark(r; ℓ, λ) σε κάθε ακτίνα και στη συνέχεια τη μετατρέπει σε ^VcDM(R).
3. Ελαχιστοποιεί το χ²(ℓ, λ) έναντι της καμπύλης περιστροφής του Γαλαξία μας της εποχής Gaia για να βρει αντιπροσωπευτικές παραμέτρους με την καλύτερη δυνατή προσαρμογή.

Η προσομοίωση έχει σχεδιαστεί για να είναι αναπαραγώγιμη. Μπορεί να εκτελεστεί σε JavaScript ή Python χωρίς εξειδικευμένη βιβλιοθήκη αστροφυσικής.

Συμβολισμοί που χρησιμοποιούνται σε όλο το κείμενο

• R είναι η κυλινδρική γαλακτοκεντρική ακτίνα στο επίπεδο του δίσκου, σε kpc.
• r είναι η σφαιρική γαλακτοκεντρική ακτίνα.
• z είναι το ύψος πάνω από το δίσκο.
• ℓ είναι το μήκος συνοχής BeeTheory, σε kpc.
• λ είναι η χωρίς διάσταση σύζευξη κύματος-μάζας.
• Το G χρησιμοποιείται σε μονάδες kpc km² s-² M⊙-¹.

1. Δεδομένα παρατήρησης – Gaia 2024

Το σύνολο δεδομένων βασίζεται στην καμπύλη περιστροφής του Γαλαξία μας της εποχής Gaia. Χρησιμοποιεί ακτίνες R, κυκλικές ταχύτητες _Vc και αβεβαιότητες σ. Η αρχική τεχνική έκδοση χρησιμοποιούσε 16 σημεία δεδομένων που κάλυπταν R = 4-27,3 kpc.

Τα σημαντικά στοιχεία παρατήρησης είναι:

• _Vc(R⊙ = 8 kpc) ≈ 230 km/s, το σημείο αγκύρωσης κοντά στην τροχιά του Ήλιου.
• Η _Vc είναι περίπου επίπεδη από περίπου 5 έως 20 kpc.
• Η _Vc μειώνεται στον εξωτερικό μετρούμενο δίσκο, φτάνοντας περίπου 173 ± 17 km/s κοντά στα 27,3 kpc στα δεδομένα που αναφέρονται.

Αυτό απαιτεί μια αποτελεσματική κατανομή κρυμμένης μάζας που αυξάνεται έντονα σε ενδιάμεσες ακτίνες και στη συνέχεια εξασθενεί σε μεγαλύτερες ακτίνες.

Εξαιρούμε τον εσωτερικότερο Γαλαξία, επειδή εκεί κυριαρχούν η κεντρική μπάρα και οι μη κυκλικές κινήσεις. Ένα απλουστευμένο αξονοσυμμετρικό μοντέλο δεν είναι αξιόπιστο στο εσωτερικό αυτής της περιοχής.

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 Αντικειμενική συνάρτηση

Η προσαρμογή ελαχιστοποιεί το μειωμένο τετράγωνο chi:

με Ν = 16 σημεία δεδομένων και p = 2 ελεύθερες παραμέτρους, ℓ και λ. Αυτό δίνει 14 βαθμούς ελευθερίας.

5.2 Αναζήτηση πλέγματος δύο περασμάτων

Χρησιμοποιείται αναζήτηση με πλέγμα και όχι κάθοδος κλίσης, επειδή το τοπίο έχει μια μεγάλη, καμπυλωτή κοιλάδα εκφυλισμού μεταξύ λ και ℓ.

• Πέρασμα 1: χονδροειδές πλέγμα πάνω από το ℓ και το λ.
• Πέρασμα 2: τοπική βελτίωση γύρω από το χονδροειδές ελάχιστο.

Η αντιπροσωπευτική περιοχή βέλτιστης προσαρμογής είναι:

5.3 Η κοιλάδα εκφυλισμού

Το τοπίο χ² δεν είναι ένα συμμετρικό μπολ. Σχηματίζει μια επιμήκη κοιλάδα επειδή η ομαλοποίηση της κορυφαίας πυκνότητας εξαρτάται έντονα από την ισχύ της σύζευξης και μόνο ασθενώς από το μήκος συνοχής μέσα στο καθεστώς επίπεδης περιστροφής.

Η εξωτερική πτώση της καμπύλης περιστροφής σπάει αυτόν τον εκφυλισμό, επειδή οι μικρότερες τιμές ℓ καταστέλλουν την αποτελεσματική πυκνότητα νωρίτερα.

Δεδομένα πέραν των 30 kpc, συμπεριλαμβανομένων σφαιρικών σμηνών, αστρικών ρευμάτων, αστέρων φωτοστέφανου και δορυφορικών γαλαξιών, είναι απαραίτητα για τη στενότερη δέσμευση του ℓ.

7. Πλήρης κωδικός αναφοράς

Η ακόλουθη εφαρμογή αναφοράς JavaScript αναπαράγει την κύρια λογική της προσομοίωσης. Προορίζεται για τεχνική επικύρωση και θα πρέπει να τοποθετείται σε κατάλληλο περιβάλλον δέσμης ενεργειών, όχι απευθείας μέσα σε ένα τυπικό μπλοκ περιεχομένου του WordPress, εκτός αν επιτρέπονται προσαρμοσμένες δέσμες ενεργειών.

// Φυσικές σταθερές
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
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],

// Baryonic velocity placeholder
function vBaryonic(R) {
  // Στην πλήρη εφαρμογή, χρησιμοποιείται ο τύπος του δίσκου του Freeman
  // με τροποποιημένες συναρτήσεις Bessel 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,
  let 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),
}

// Chi-squared
function chiSquared(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),
}

Μια πλήρης έκδοση θα πρέπει να περιλαμβάνει ακριβείς εφαρμογές των τροποποιημένων συναρτήσεων Bessel _I0, _I1, _K0 και _K1 για την ταχύτητα του δίσκου Freeman.

Αναφορές

• Abramowitz, M., Stegun, I. A. - Handbook of Mathematical Functions, Dover, 1972.
• Bovy, J., Rix, H.-W. - A Direct Dynamical Measurement of the Milky Way's Disk Surface Density Profile, ApJ 779, 115, 2013.
• Freeman, K. C. - On the disks of spiral and S0 galaxies, ApJ 160, 811, 1970.
• McMillan, P. J. - The mass distribution and gravitational potential of the Milky Way, MNRAS 465, 76, 2017.
• Ou, X., Eilers, A.-C., Necib, L., Frebel, A. - The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693, 2024.
• Pato, M., Iocco, F., Bertone, G. - Dynamical constraints on the dark matter distribution in the Milky Way, JCAP 12, 001, 2015.
• Portail, M. et al. - Dynamical modelling of the Galactic bulge and bar, MNRAS 465, 1621, 2017.