BeeTheory – Στατιστική ανάλυση – 2025

Διαβάζοντας το Fit

χ²/dof = 0,475

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

Τι είναι το χ²/dof;

Όταν ένα μοντέλο προβλέπει τιμές _Vmodel(_Ri) και τα δεδομένα παρέχουν παρατηρούμενες τιμές _Vobs(_Ri) με αβεβαιότητες μέτρησης _σi, το μειωμένο τετράγωνο chi είναι:

\(\frac{\chi^2}{\mathrm{dof}}=\frac{1}{N-p}\sum_{i=1}^{N}\left(\frac{V_{\mathrm{model}}(R_i)-V_{\mathrm{obs}}(R_i)}{\sigma_i}\right)^2\)

όπου Ν είναι ο αριθμός των σημείων δεδομένων και p είναι ο αριθμός των ελεύθερων παραμέτρων.

Ορίστε: Ν = 16 σημεία Gaia 2024, p = 2 προσαρμοσμένες παράμετροι, K και α, οπότε dof = 14.

Κάθε όρος στο άθροισμα είναι ένα τράβηγμα: το υπόλειμμα εκφρασμένο σε μονάδες της αβεβαιότητας μέτρησης.

Ένα τράβηγμα 0,5 σημαίνει ότι το μοντέλο απέχει κατά μισό σίγμα – εξαιρετικό. Μια απόκλιση 2,0 σημαίνει απόκλιση δύο σίγμα – αξίζει να διερευνηθεί.

\(\frac{\chi^2}{\mathrm{dof}}\approx1\quad\Rightarrow\quad\text{μοντέλο ταιριάζει στα δεδομένα στο επίπεδο των αβεβαιοτήτων τους}\) \(\frac{\chi^2}{\mathrm{dof}}\ll1\quad\Rightarrow\quad\text{μοντέλο ταιριάζει πολύ καλά: οι αβεβαιότητες μπορεί να υπερεκτιμώνται, ή υπερπροσαρμογή}\) \(\frac{\chi^2}{\mathrm{dof}}\gg1\quad\Rightarrow\quad\text{το μοντέλο είναι λάθος, ή οι αβεβαιότητες είναι υποεκτιμημένες}\)

Ο αριθμός μας: 0.475

0.475

χ² / dof – N = 16 – p = 2 – dof = 14

Στατιστικά άριστα

Το μοντέλο δεν είναι λάθος

Οι αβεβαιότητες μπορεί να είναι ελαφρώς μεγάλες

χ²/dof κλίμακα

0,475 ← μας

1.0 ιδανικό

00.511.522.53+

< 0.3
over-fit

0.3-0.7
εξαιρετική

0.7-1.3
καλή

1.3-2.0
αποδεκτό

> 2.0
κακή

Η τιμή 0,475 βρίσκεται στη ζώνη της άριστης ποιότητας. Το μοντέλο δεν είναι λάθος – το μέσο υπόλοιπο είναι μόνο 0,69σ και για τα 16 σημεία δεδομένων.

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

Αλλά το “εξαιρετικό” δεν είναι το ίδιο με το “αποδεδειγμένο”.

χ²/dof = 0,475 θα μπορούσε επίσης να σημαίνει ότι οι αβεβαιότητες μέτρησης _σi είναι ελαφρώς υπερεκτιμημένες. Εάν τα πραγματικά σφάλματα ήταν 30% μικρότερα, το ίδιο μοντέλο θα έδινε χ²/dof ≈ 1,0.

Δεν μπορούμε να διακρίνουμε μεταξύ “το μοντέλο είναι πολύ καλό” και “τα σφάλματα είναι ελαφρώς γενναιόδωρα” μόνο από το χ². Αυτή είναι μια τυπική στατιστική ασάφεια.

The Residuals – Σημείο προς σημείο

Ο ακατέργαστος αριθμός 0,475 κρύβει πληροφορίες. Η εξέταση των επιμέρους έλξεων αποκαλύπτει τη δομή της προσαρμογής: ποια σημεία το μοντέλο καρφώνει και πού δυσκολεύεται.

Pull plot: \((V_{\mathrm{model}}-V_{\mathrm{obs}})/\sigma_i\) για κάθε σημείο δεδομένων Gaia 2024

R (kpc)	Vobs	σ	Vmodel	Υπόλοιπο	Τραβήξτε
Φόρτωση υπολειμμάτων…

15 / 16 σημεία εντός 1σ

Κάθε σημείο δεδομένων εκτός από το R = 27.3 kpc έχει |pull| < 1.0. Σε ένα καλά βαθμονομημένο μοντέλο με γκαουσιανά σφάλματα, περιμένουμε περίπου 68% των σημείων εντός 1σ - εδώ έχουμε 94%.

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

Το ακραίο σημείο στο R = 27,3 kpc – pull = +1,53σ

Το μοντέλο προβλέπει _Vc = 195,3 km/s, ενώ το Gaia μετρά 173 ± 17 km/s.

Η απόκλιση είναι 22,3 km/s ή 1,53σ. Αυτό δεν είναι στατιστικά ανησυχητικό, αλλά είναι φυσικά σημαντικό: το μοντέλο μειώνεται πολύ αργά στη μεγαλύτερη ακτίνα.

Η παράμετρος καλύτερης προσαρμογής K = 0,01349

Η σταθερά σύζευξης K = 0,01349 kpc-¹ είναι το πλάτος του πεδίου σκοτεινής μάζας BeeTheory που δημιουργείται ανά μονάδα ορατής μάζας.

Προσδιορίστηκε ελαχιστοποιώντας το χ² στα 16 σημεία δεδομένων Gaia με α = 0,0744 kpc-¹ που καθορίστηκε από τη συνδυασμένη προσαρμογή Newby + Gaia.

\(K=0.01349\,\mathrm{kpc}^{-1}\) \(\lambda=K\ell^2=\frac{K}{\alpha^2}=\frac{0.01349}{(0.0744)^2}=2.44\) \(\rho_{\mathrm{dark}}(R_\odot=8\,\mathrm{kpc})=0.37\,\mathrm{GeV/cm^3}\)

Η παρατηρούμενη τοπική πυκνότητα σκοτεινής ύλης είναι:

\(\rho_{\mathrm{obs}}(R_\odot)=0.39\pm0.03\,\mathrm{GeV/cm^3}\)

Η άνευ διαστάσεων σύζευξη λ = 2,44 βρίσκεται στο εύρος λ ∈ [2, 7] που παρατηρείται σε όλες τις προσαρμογές της BeeTheory, συμπεριλαμβανομένων συγκρίσεων ενός συστατικού, πολλαπλών συστατικών και ατομικού H₂.

Αυτή η καθολικότητα – η ίδια τάξη μεγέθους από την κλίμακα του femtometre έως την κλίμακα του kiloparsec – είναι ο ισχυρότερος έλεγχος εσωτερικής συνέπειας του πλαισίου BeeTheory.

Γιατί το Κ είναι μικρότερο εδώ από ό,τι στις προσαρμογές ενός δίσκου

Προηγούμενες προσαρμογές που χρησιμοποιούσαν μόνο τον λεπτό δίσκο ως πηγή έδωσαν_Ksingle = 0,038 kpc-¹. Με τις έξι γαλαξιακές συνιστώσες να συνεισφέρουν στο σκοτεινό πεδίο, η συνολική βαρυονική πηγή είναι 2,8 φορές μεγαλύτερη.

Επομένως, το Κ πρέπει να είναι αναλογικά μικρότερο για να παράγει την ίδια σκοτεινή μάζα.

\(K_{\mathrm{multi}}\approx\frac{K_{\mathrm{single}}}{2.8}=\frac{0.038}{2.8}=0.0136\,\mathrm{kpc}^{-1}\)

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

Γιατί το 27.3 kpc είναι δύσκολο για κάθε μοντέλο

Το απώτατο σημείο δεδομένων του Gaia 2024, _Vc(27,3 kpc) = 173 ± 17 km/s, δεν είναι μόνο δύσκολο για την BeeTheory. Είναι το δυσκολότερο σημείο για κάθε μοντέλο σκοτεινής ύλης που εφαρμόζεται στην καμπύλη περιστροφής του Γαλαξία μας.

Μοντέλο	Vc(27,3 kpc) προβλεπόμενη	Υπόλοιπο	Τραβήξτε	χ²/dof
BeeTheory	195,3 km/s	+22.3	+1.53σ	0.475
Προφίλ NFW	~198 km/s	+25	+1.5σ	0.44
Einasto α = 0,91	~191 km/s	+18	+1.1σ	0.38
Ισόθερμο φωτοστέφανο	~218 km/s	+45	+2.6σ	2.6

Κανένα τυπικό μοντέλο δύο παραμέτρων δεν αναπαράγει το _Vc(27.3 kpc) = 173 km/s εντός 1σ.

Υπάρχουν τρεις εύλογες εξηγήσεις.

Η ίδια η μέτρηση: R = 27,3 kpc είναι το πιο μακρινό σημείο, με τους λιγότερους κινηματικούς ιχνηλάτες και τις μεγαλύτερες συστηματικές αβεβαιότητες. Το Gaia DR4 μπορεί να μετατοπίσει την τιμή.
Συμβολή του δίσκου αερίων: και συμβάλλει στη βαρυονική μάζα σε μεγάλες ακτίνες. Η συμπερίληψή του ως ξεχωριστή συνιστώσα θα μπορούσε να αυξήσει την απότομη πτώση.
Ένα μικρότερο μήκος συνοχής: α = 0,12 kpc-¹, ή ℓ = 8,3 kpc, ταιριάζει καλύτερα στο εξωτερικό σημείο, αλλά επιδεινώνει το εσωτερικό οροπέδιο.

Πώς θα έμοιαζε ένα πραγματικά σταθερό αποτέλεσμα

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

0.475

Τρέχον χ²/dof, 16 σημεία, 2 παράμετροι

< 0.8

Στόχος: καλή προσαρμογή στο εκτεταμένο σύνολο δεδομένων

~1.0

Ιδανικό: βαθμονομημένο σε δείγματα γαλαξιών

Απαιτήσεις για μια σταθερή επικύρωση

✓ Καλό χ²/dof στο σύνολο δεδομένων εκπαίδευσης: Επιτεύχθηκε. χ²/dof = 0,475 στο Gaia 2024.
✓ Σωστή τοπική πυκνότητα σκοτεινής ύλης: Επιτεύχθηκε. 0,37 GeV/cm³ προβλεπόμενο έναντι 0,39 ± 0,03 παρατηρούμενο.
✓ Συνεπής σύζευξη χωρίς διαστάσεις: λ = 2,44 σε σύγκριση με λ ≈ 3,5 από τη βαθμονόμηση H₂.
! Αναπαράγετε το εξώτατο σημείο Gaia εντός 1σ: Δεν έχει επιτευχθεί πλήρως. Το υπόλοιπο στο R = 27,3 kpc είναι +1,53σ.
✗ Τυφλή πρόβλεψη σε άλλους γαλαξίες: Δεν έχει γίνει ακόμα. Ο κατάλογος SPARC παρέχει τη φυσική δοκιμή.
✗ Παραγωγή του ℓ από τις πρώτες αρχές: Δεν έχει γίνει ακόμη. Το ℓ είναι προς το παρόν προσαρμοσμένο, όχι εξαγόμενο.
✗ Επικύρωση σε κλίμακα συστάδας: Δεν έχει γίνει ακόμα. Ο φακός του σμήνους σφαίρας είναι μια βασική δοκιμή.
→ Εκτεταμένη καμπύλη περιστροφής πέραν των 30 kpc: Gaia DR4 θα πρέπει να παρέχει μια κρίσιμη βραχυπρόθεσμη δοκιμή.

Η ειλικρινής επιστημονική κατάσταση

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

Δεν είναι ακόμη θεωρία με την αυστηρή έννοια του όρου, επειδή τα K και ℓ είναι προσαρμοσμένα και όχι παραγόμενα.

Ο δρόμος προς μια πλήρη θεωρία είναι σαφής: παράγουμε το ℓ = f(_Lsource) από το αξίωμα της κυματικής μάζας, δοκιμάζουμε καθολικά στο SPARC και κάνουμε μια τυφλή πρόβλεψη για το Gaia DR4.

Δεδομένα: κ.ά., MNRAS 528, 2024. Μοντέλο: Dutertre 2023, προσαρμογή πολλαπλών συνιστωσών. Αναφορά SPARC: Lelli, McGaugh, Schombert, AJ 152, 2016. Bullet Cluster: Clowe et al., ApJL 648, 2006.