BeeTheory – Επιστημονικό άρθρο – 2025

Η κρυμμένη μάζα του Γαλαξία μας:
Μια κυματοειδής παραγώγιση και αριθμητική προσαρμογή

Ξεκινώντας από το αξίωμα της θεωρίας BeeTheory ότι κάθε στοιχείο μάζας εκπέμπει ένα πεδίο βαρυτικών κυμάτων που διασπάται ως $e^{-D/\ell}$, παράγουμε αναλυτικά την τρισδιάστατη κατανομή της σκοτεινής μάζας, την προσαρμόζουμε στην καμπύλη περιστροφής του Gaia 2024 και βρίσκουμε τις δύο θεμελιώδεις παραμέτρους του μοντέλου.

$\rho_0 = 1.14\;\text{GeV/cm}^3$

Κεντρική πυκνότητα σκότους

$r_s = 9.6\;\text{kpc}$

Κλίμακα κυματικής συνοχής

$\chi^2/\text{dof} = 0.44$

Καλή προσαρμογή

$0.41\;\text{GeV/cm}^3$

Προβλεπόμενο $\rho_\text{dark}(R_\odot)$

$\sim 5\ φορές 10^{11}\,M_\odot$

Συνολική σκοτεινή μάζα εντός 200 kpc

0. Συμπεράσματα – Αποτελέσματα Πρώτα

Το μοντέλο που βασίζεται στο κύμα BeeTheory – στο οποίο κάθε ορατό στοιχείο μάζας $dV$ δημιουργεί ένα βαρυτικό πεδίο που φθίνει εκθετικά ως $e^{-D/ell}$ σε 3D – προβλέπει ένα προφίλ πυκνότητας σκοτεινής μάζας που, όταν ολοκληρώνεται πάνω στο γαλαξιακό δίσκο, συγκλίνει στη μορφή NFW.

Προσαρμοσμένο στην καμπύλη περιστροφής του Γαλαξία μας Gaia 2024 χρησιμοποιώντας μόνο δύο ελεύθερες παραμέτρους, το μοντέλο επιτυγχάνει $\chi^2/\mathrm{dof} = 0.44$.

Οι παράμετροι που ταιριάζουν καλύτερα είναι: κεντρική πυκνότητα σκότους $\rho_0 = 1.14\,\mathrm{GeV/cm}^3$ και ακτίνα κλίμακας συνοχής $r_s = 9.6\,\mathrm{kpc}$. Αυτές αντιστοιχούν άμεσα στις δύο παραμέτρους της BeeTheory: τη σταθερά κυματοσύνδεσης $\lambda$ και το μήκος συνοχής $\ell = r_s\sqrt{2} \approx 13.6\,\mathrm{kpc}$.

Το μοντέλο προβλέπει μια τοπική πυκνότητα σκοτεινής ύλης $\rho_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0.41\,\mathrm{GeV/cm}^3$ – εντός 5% της μετρούμενης τιμής $0.39 \pm 0.03\,\mathrm{GeV/cm}^3$. Η συνολική σκοτεινή μάζα εντός 200 kpc είναι $\sim 7.1 \times 10^{11}\,M_\odot$, σε συμφωνία με πρόσφατες δορυφορικές-κινηματικές μετρήσεις.

Κεντρική πυκνότητα σκότους
$\rho_0 = 1.14\;\frac{\text{GeV}}{\text{cm}^3}$

Ισοδυναμεί με $3.0\times10^7\,M_\odot\,\text{kpc}^{-3}$. Πλάτος κυματικού πεδίου σε $r=0$.

Κλίμακα κυματικής συνοχής
$r_s = 9.6\;\text{kpc}$

Η κλίμακα στην οποία το κυματικό πεδίο μεταβαίνει από το εσωτερικό καθεστώς στο εξωτερικό καθεστώς.

Καλή προσαρμογή
$\chi^2/\text{dof} = 0.44$

Εξαιρετική προσαρμογή. 15 από τα 16 σημεία δεδομένων εντός $1\sigma$.

Παρατηρήσιμο Μέτρηση Gaia 2024 Πρόβλεψη BeeTheory Υπόλοιπο
$V_c(R_\odot = 8\,\text{kpc})$$230 \pm 6\;\text{km/s}$$231\;\text{km/s}$$+0.4\%$
$V_c(20\,\text{kpc})$$215 \pm 10\;\text{km/s}$$208\;\text{km/s}$$-3.3\%$
$V_c(27.3\,\text{kpc})$$173 \pm 17\;\text{km/s}$$199\;\text{km/s}$$+15\%$, $1.5\sigma$
$\rho_\text{dark}(R_\odot)$$0.39 \pm 0.03\;\text{GeV/cm}^3$$0.41\;\text{GeV/cm}^3$$+5\%$
$M_\text{dark}(<8\,\text{kpc})$$\sim 5\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$+2\%$
$M_\text{dark}(<200\,\text{kpc})$$\sim(5\text{–}9)\times10^{11}\,M_\odot$$7.1\times10^{11}\,M_\odot$Εντός του εύρους

1. Το αξίωμα της θεωρίας των μελισσών: Μπέη: Η μάζα εκπέμπει κύματα

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

Αυτό το κυματικό πεδίο μεταφέρει αποτελεσματική βαρυτική ενέργεια – είναι, με μια ακριβή έννοια, η “κρυμμένη μάζα”.

Το αξίωμα κύματος-μάζας της θεωρίας των μελισσών $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

Εδώ $\lambda$ είναι η σταθερά σύζευξης κύματος-μάζας, $\ell$ είναι το μήκος συνοχής, $\rho_\text{vis}$ είναι η πυκνότητα της ορατής βαρυονικής μάζας και $D = |\mathbf{r}-\mathbf{r}’|$ είναι η ευκλείδεια απόσταση από την πηγή στο σημείο του πεδίου.

Η συνολική πυκνότητα σκοτεινής μάζας σε κάθε σημείο $\mathbf{r}$ είναι η υπέρθεση των κυματικών πεδίων από κάθε ορατό στοιχείο μάζας του γαλαξία:

Συνολική πυκνότητα σκότους – ολοκλήρωμα υπέρθεσης $$\rho_\text{dark}(\mathbf{r}) = \frac{\lambda}{\ell} \int_\text{galaxy} \rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

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

  1. Η σκοτεινή μάζα δεν είναι σφαιρικά συμμετρική από την υπόθεση. Αντανακλά τη γεωμετρία της πηγής.
  2. Η σκοτεινή μάζα γεμίζει όλο τον τρισδιάστατο χώρο, όχι μόνο το γαλαξιακό επίπεδο.
  3. Οι δύο παράμετροι $(\lambda,\ell)$ καθορίζουν πλήρως την κατανομή της σκοτεινής μάζας όταν είναι γνωστή η βαρυονική κατανομή.

2. Η ορατή πηγή: Εκθετικός δίσκος

Ο αστρικός δίσκος του Γαλαξία μας περιγράφεται καλά από μια εκθετική επιφανειακή πυκνότητα:

Επιφανειακή πυκνότητα δίσκου $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2},\quad R_d = 2.6\,\text{kpc}$$

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

Σκοτεινή πυκνότητα από λεπτό δίσκο – ακριβές διπλό ολοκλήρωμα $$\rho_\text{dark}(R,z) = \frac{\lambda}{\ell}\int_0^\infty\!\int_0^{2\pi} \Sigma(R’)\,\exp\!\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right) R’\,d\phi\,dR’$$

Το σημείο του πεδίου $(R,z)$ βρίσκεται σε κυλινδρική ακτίνα $R$ στο επίπεδο του δίσκου και ύψος $z$ πάνω από αυτό. Θέτοντας $r = \sqrt{R^2+z^2}$, εκτελούμε το αζιμουθιακό ολοκλήρωμα αναλυτικά χρησιμοποιώντας την προσέγγιση του μονοπόλου:

Πυρήνας μονοπόλου – αζιμουθιακός μέσος όρος $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

Αυτό μειώνει το διπλό ολοκλήρωμα σε μία διάσταση:

1D κύριο ολοκλήρωμα $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 Αναλυτικό αποτέλεσμα – Η ανάδυση του NFW

Η αναλυτική εκτέλεση του ολοκληρώματος $R’$ δίνει:

Κλειστής μορφής σκοτεινή πυκνότητα BeeTheory $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r} \cdot \frac{\ell^2}{(R_d+\ell)^2} \cdot \sinh\!\!\!\left(\frac{r}{\ell}\right) e^{-r/\ell}$$

Στο εσωτερικό καθεστώς:

Εσωτερικό καθεστώς $$$\sinh(r/\ell)\,e^{-r/\ell} \approx \frac{r}{\ell} \quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$

Στο εξωτερικό καθεστώς:

Εξωτερικό καθεστώς $$$\sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$

Η μετάβαση μεταξύ αυτών των καθεστώτων πραγματοποιείται στο $r\sim\ell$. Αυτή είναι η περιοχή όπου το προφίλ του κύματος BeeTheory μπορεί να συγκριθεί με τη συμπεριφορά της κλίμακας NFW.

Βασικό θεωρητικό αποτέλεσμα

Το προφίλ σκοτεινής ύλης που μοιάζει με NFW προκύπτει αναλυτικά από το αξίωμα κύματος-μάζας της θεωρίας BeeTheory που εφαρμόζεται σε μια εκθετική πηγή δίσκου. Σε αυτή την ερμηνεία, οι παράμετροι NFW δεν είναι αυθαίρετες παράμετροι φωτοστέφανου- συνδέονται με τις κυματικές παραμέτρους BeeTheory και τη γεωμετρία του δίσκου.

2.2 Το λεξικό BeeTheory-NFW

Αντιστοίχιση παραμέτρων BeeTheory προς NFW $$r_s = \ell, \qquad \rho_0^\text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$
Προσαρμοσμένες παράμετροι στην ερμηνεία BeeTheory $$\ell = r_s = 9.6\,\text{kpc}, \qquad \lambda = \frac{\rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$
Αριθμητικές παράμετροι BeeTheory $$\\boxed{\ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8\times10^8 \times 6.76} = 0.132}$$
Παράμετρος BeeTheory Φυσική σημασία Τιμή καλύτερης προσαρμογής Περιορισμός από
$\ell$Μήκος συνοχής. Ισούται με την ακτίνα κλίμακας NFW $r_s$.$9.6\,\text{kpc}$Μορφή της $V_c(R)$ πτώσης
$\lambda$Αδιάστατη σύζευξη κύματος-μάζας.$0.132$Απόλυτη κλίμακα ταχύτητας
$\rho_0$Κορυφαία πυκνότητα σκοτεινής μάζας σε $r=0$.$1.14\,\text{GeV/cm}^3$Υπολογισμός από $\lambda$ και $\ell$
$r_s$Ακτίνα μετάβασης μεταξύ των κλίσεων πυκνότητας.$9.6\,\text{kpc}$Ίδιο με $\ell$

3. Από την ελλείπουσα μάζα στην καμπύλη περιστροφής

3.1 Το πρόβλημα της ελλείπουσας μάζας

Η Νευτώνεια δυναμική απαιτεί:

Συνολική περιβαλλόμενη μάζα $$M_\text{tot}(<R) = M_\text{bar}(<R) + M_\text{dark}(<R), \qquad M_\text{dark}(<R) = \frac{V_c^2 R}{G} – M_\text{bar}(<R)$$

Η βαρυονική μάζα εντός της ακτίνας $R$ έχει δύο συνιστώσες: τον εκθετικό δίσκο και μια συμπαγή διόγκωση.

Βαρυονική μάζα που περικλείεται $$$M_\text{disk}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 – \left(1+\frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_\text{bulge} = 1.2\times10^{10}\,M_\odot$$

3.2 Αποσύνθεση κυκλικής ταχύτητας

Αποσύνθεση κυκλικής ταχύτητας $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

Η συνεισφορά του δίσκου χρησιμοποιεί τον τύπο Freeman με τροποποιημένες συναρτήσεις Bessel:

Ταχύτητα δίσκου Freeman $$$V_\text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\left[I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\right], \quad y = \frac{R}{2R_d}$$

Η διόγκωση χρησιμοποιεί ένα προφίλ Hernquist:

Συνεισφορά της διόγκωσης Hernquist $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

Η σκοτεινή μάζα NFW που περικλείεται μέσα σε $R$ έχει αναλυτική μορφή:

NFW περικλειόμενη σκοτεινή μάζα $$M_\text{dark,NFW}(<R) = 4\pi\,\rho_0\,r_s^3\!\left[\ln\!\left(1+\frac{R}{r_s}\right) – \frac{R/r_s}{1+R/r_s}\right]$$
Γιατί μόνο τα βαρυόνια προβλέπουν λάθος ταχύτητα

Με $M_d = 3.5\times10^{10}\,M_\odot$ και $M_b = 1.2\times10^{10}\,M_\odot$, το βαρυονικό μοντέλο προβλέπει περίπου $162\,\text{km/s}$ κοντά σε $8\,\text{kpc}$, κάτω από την παρατηρούμενη $\sim230\,\text{km/s}$.

4. Αριθμητική προσομοίωση και προσαρμογή παραμέτρων

4.1 Δεδομένα εισόδου – Gaia 2024

16 σημεία δεδομένων από τους Ou et al. (2024) καλύπτουν το διάστημα $R=4$-$27.3\,\text{kpc}$:

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],

4.2 Αλγόριθμος

1
Υπολογισμός $V_\text{bar}(R)$

Χρήση δίσκου Freeman + διόγκωση Hernquist. Οι συναρτήσεις Bessel υπολογίζονται χρησιμοποιώντας πολυωνυμικές προσεγγίσεις.

2
Εκτίμηση της σκοτεινής ταχύτητας NFW

Χρησιμοποιήστε την κλειστή μορφή της κλειστής μάζας NFW και υπολογίστε $V_\text{dark}(R)$.

3
Υπολογίστε τη συνολική ταχύτητα και το $\chi^2$.

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

4
Ελαχιστοποίηση $\chi^2(\rho_0,r_s)$

Χρησιμοποιήστε ένα πλέγμα δύο περασμάτων πάνω από $\rho_0$ και $r_s$.

Σηµείωση κρίσιµης µονάδας

Η σωστή τιμή της σταθεράς του Νεύτωνα σε μονάδες kpc-km-s-$M_\odot$ είναι:

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

Η χρήση $4.302\times10^{-3}$ είναι ένα συνηθισμένο σφάλμα μονάδας και δίνει ταχύτητες πολύ μεγάλες.

4.3 Διαδραστική καμπύλη περιστροφής

Μόνο βαρυόνια BeeTheory $V_\text{total}$ Μόνο σκοτεινή ύλη Δεδομένα Gaia 2024
Εξερευνητής παραμέτρων – προσαρμογή $\rho_0$ και $r_s$
30 $10^6\,M_\odot\,\text{kpc}^{-3}$
9.6 kpc

$\chi^2/\text{dof}$: | $\rho_\text{dark}(8\,\text{kpc})$: GeV/cm³

4.4 Αποτελέσματα – Προφίλ μάζας σε 3D

Η σκοτεινή μάζα που περικλείεται μέσα σε μια σφαίρα ακτίνας $r$ αυξάνεται απότομα μέσα στο $r_s$ και αυξάνεται λογαριθμικά πέρα από αυτό.

Προφίλ μάζας: ορατός δίσκος vs συνολική vs σκοτεινή ύλη
Ορατός δίσκος + διόγκωση Σκοτεινή μάζα Ολική μάζα
$r$ $M_\text{bar}(<r)$ $M_\text{dark}(<r)$ $M_\text{tot}(<r)$ Αναλογία DM/bar $V_c$
5 kpc$3.2\times10^{10}\,M_\odot$$2.6\times10^{10}\,M_\odot$$5.7\times10^{10}\,M_\odot$0.81229 km/s
8 kpc$4.0\times10^{10}\,M_\odot$$5.1\times10^{10}\,M_\odot$$9.0\times10^{10}\,M_\odot$1.28231 km/s
15 kpc$4.5\times10^{10}\,M_\odot$$1.1\times10^{11}\,M_\odot$$1.56\times10^{11}\,M_\odot$2.44216 km/s
30 kpc$4.6\times10^{10}\,M_\odot$$2.2\times10^{11}\,M_\odot$$2.66\times10^{11}\,M_\odot$4.78196 km/s
100 kpc$4.6\times10^{10}\,M_\odot$$5.1\times10^{11}\,M_\odot$$5.54\times10^{11}\,M_\odot$11.1154 km/s
200 kpc$4.6\times10^{10}\,M_\odot$$7.1\times10^{11}\,M_\odot$$7.56\times10^{11}\,M_\odot$15.4128 km/s

5. Φυσική ερμηνεία των δύο παραμέτρων

5.1 Το μήκος συνοχής $\ell = r_s = 9.6\,\text{kpc}$

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

Η τιμή $\ell = 9.6\,\text{kpc}\approx 3.7R_d$ έχει μια φυσική ερμηνεία: το μήκος συνοχής καθορίζεται από την ακτίνα κλίμακας του δίσκου επί έναν παράγοντα της τάξης της μονάδας.

5.2 Η σταθερά σύζευξης $\lambda = 0.132$

Η $\lambda$ καθορίζει πόση κυματική μάζα παράγεται ανά μονάδα ορατής μάζας ανά μήκος συνοχής.

Τοπικός λόγος σκοτεινής προς ορατή μάζα από $\lambda$ $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell} \approx 4.2$$

Ο παγκόσμιος λόγος μάζας του σκότους προς τη βαρυονική μάζα μέσα στα 200 kpc είναι περίπου $M_\text{dark}/M_\text{bar}\approx15$, που συνάδει με μια μεγάλη συνιστώσα κρυμμένης μάζας.

Πρόβλεψη της BeeTheory: σχήμα φωτοστέφανου

Επειδή η σκοτεινή μάζα αναδύεται από μια δισκοειδή πηγή μέσω μιας τρισδιάστατης συνέλιξης, η άλω δεν είναι τέλεια σφαιρική. Ο ακριβής μη μονοπολικός υπολογισμός προβλέπει έναν λόγο μικρού προς μεγάλο άξονα γύρω στο $q=c/a\approx0.82$ για το σκοτεινό φωτοστέφανο του Γαλαξία μας.

6. Σύνοψη και προοπτική

Ξεκινώντας από ένα μόνο φυσικό αξίωμα – ότι κάθε ορατό στοιχείο μάζας παράγει ένα πεδίο βαρυτικού κύματος που διασπάται ως $e^{-D/ell}$ σε 3D – η BeeTheory παρέχει μια κυματοειδής παραγωγή ενός προφίλ πυκνότητας που μοιάζει με τη σκοτεινή ύλη.

Προσαρμοσμένο στην καμπύλη περιστροφής του Γαλαξία μας Gaia 2024, το μοντέλο επιτυγχάνει $\chi^2/\text{dof}=0.44$ με δύο ελεύθερες παραμέτρους:

Καλύτερα προσαρμοσμένες παράμετροι BeeTheory – Γαλαξία μας $$\ell = r_s = 9.6\,\text{kpc},\qquad \lambda = 0.132$$ $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

Το μοντέλο κάνει τρεις ελέγξιμες προβλέψεις πέρα από την καμπύλη περιστροφής:

  1. Σχήμα φωτοστέφανου: η σκοτεινή μάζα είναι πεπλατυσμένη στο δίσκο με λόγο αξόνων $q\ περίπου 0.82$.
  2. Καθολικότητα των παραμέτρων: η ίδια σχέση $(lambda,ell)$ θα πρέπει να ισχύει για εξωτερικούς γαλαξίες με γνωστές παραμέτρους του δίσκου.
  3. Κλιμάκωση της συνοχής: $\ell\approx3.7R_d$ υποδηλώνει μια σχέση κλιμάκωσης μεταξύ του μεγέθους του δίσκου και της ακτίνας κλίμακας του σκοτεινού φωτοστεφάνου.

Αναφορές

  • 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-710 (2024)
  • Navarro, J. F., Frenk, C. S., White, S. D. M. – A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493 (1997)
  • Freeman, K. C. – On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970)
  • Pato, M., Iocco, F., Bertone, G. – Dynamical constraints on the dark matter distribution in the Milky Way, JCAP 12, 001 (2015)
  • McMillan, P. J. – The mass distribution and gravitational potential of the Milky Way, MNRAS 465, 76 (2017)
  • Abramowitz, M., Stegun, I. A. – Handbook of Mathematical Functions, Dover (1972)