Μαθηματικά θεμέλια της γαλαξιακής ελλείπουσας μάζας: Δίσκος, σφαίρα, πυκνότητα, δυναμικό και ακτινική κλιμάκωση
TL;DR: Το πρόβλημα της χαμένης μάζας εμφανίζεται όταν η μάζα που προκύπτει από τις γαλαξιακές καμπύλες περιστροφής υπερβαίνει τη μάζα που παρατηρείται άμεσα στα άστρα, το αέριο και τη σκόνη. Μαθηματικά, αυτό απαιτεί τη σύνδεση της επιφανειακής πυκνότητας σε ένα δίσκο, της πυκνότητας όγκου σε τρεις διαστάσεις, του βαρυτικού δυναμικού, της ακτινικής επιτάχυνσης και της εγκλωβισμένης μάζας.
1. Ακτινικές συντεταγμένες και γεωμετρία
Διακρίνουμε δύο γεωμετρίες:
- Γεωμετρία δίσκου: η ορατή γαλαξιακή ύλη κατανέμεται κυρίως σε έναν λεπτό περιστρεφόμενο δίσκο.
- Σφαιρική γεωμετρία: η σκοτεινή ή ελλείπουσα μάζα συχνά μοντελοποιείται ως ένα περίπου σφαιρικό φωτοστέφανο.
Στοιχείο περιοχής δίσκου:
\[ dA = R\,dR\,d\phi \]
Στοιχείο σφαιρικού όγκου:
\[ dV = 4\pi r^2\,dr \]
Το ίδιο σύμβολο \(r\) χρησιμοποιείται συχνά για τη γαλακτοκεντρική ακτίνα, αλλά η σημασία του εξαρτάται από τη γεωμετρία. Σε ένα δίσκο, \(R\) είναι μια κυλινδρική ακτίνα. Σε ένα φωτοστέφανο, \(r\) είναι συνήθως μια σφαιρική ακτίνα.
2. Ορατή μάζα σε γαλαξιακό δίσκο
Ο ορατός δίσκος συχνά προσεγγίζεται από μια εκθετική επιφανειακή πυκνότητα:
\[ \Sigma(R)=\Sigma_0 e^{-R/R_d} \]
Η μάζα ενός δακτυλίου μεταξύ \(R\) και \(R+dR\) είναι:
\[ dM_{\\rm disk}=2\pi R\Sigma(R)\,dR \]
\[ dM_{\rm disk}=2\pi R\Sigma_0e^{-R/R_d}\,dR \]
Η αθροιστική μάζα του ορατού δίσκου είναι επομένως:
\[ M_{\\rm disk}(R)=2\pi\int_0^R R’\Sigma_0e^{-R’/R_d}\,dR’ \]
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\\frac{R}{R_d} \right) \right] \]
Σε μεγάλη ακτίνα:
\[ M_{\\rm disk}(R)\rightarrow 2\pi\Sigma_0R_d^2 \]
Η μάζα του ορατού δίσκου προσεγγίζει μια πεπερασμένη τιμή.
3. Σφαιρική μάζα και πυκνότητα όγκου
Για μια σφαιρική κατανομή μάζας, η πυκνότητα όγκου \(\rho(r)\) καθορίζει την περιβαλλόμενη μάζα:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2\,dr’ \]
Η αντίστροφη σχέση είναι:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Η σχέση αυτή είναι κεντρική στο πρόβλημα της χαμένης μάζας. Αν η μάζα που προκύπτει αυξάνεται γραμμικά με την ακτίνα, τότε η αντίστοιχη σφαιρική πυκνότητα μειώνεται ως \(1/r^2\).
4. Δυναμική μάζα από κυκλική κίνηση
Για κυκλική κίνηση, η επιτάχυνση της βαρύτητας ικανοποιεί:
\[ \frac{v(r)^2}{r}=\frac{GM(r)}{r^2} \]
Επομένως:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Για μια επίπεδη καμπύλη περιστροφής:
\[ v(r)\approx v_0 \]
\[ M_{\\rm dyn}(r)\approx \frac{v_0^2}{G}r \]
Αυτό δίνει την τυπική κλιμάκωση:
\[ M_{\\rm dyn}(r)\propto r \]
\[ \rho_{\rm dyn}(r)\propto \frac{1}{r^2} \]
5. Ορισμός ελλείπουσας μάζας
Η μάζα που λείπει είναι η διαφορά μεταξύ της δυναμικής μάζας και της ορατής μάζας:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Για έναν εκθετικό ορατό δίσκο:
\[ M_{\\rm miss}(r)= \frac{v(r)^2r}{G} – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\\frac{r}{R_d} \right) \right] \]
Για \(v(r)\approx v_0\):
\[ M_{\\rm miss}(r)\approx \frac{v_0^2}{G}r – 2\pi\Sigma_0R_d^2 \left[ 1-e^{-r/R_d} \left( 1+\\frac{r}{R_d} \right) \right] \]
Σε μεγάλες ακτίνες, ο όρος του δίσκου κορεστεί, ενώ ο δυναμικός όρος συνεχίζει να αυξάνεται περίπου ως \(r\).
6. Βαρυτικό δυναμικό σε 3D
Το Νευτώνειο βαρυτικό δυναμικό που δημιουργείται από μια σημειακή μάζα είναι:
\[ \Phi(r)=-\frac{GM}{r} \]
Το αντίστοιχο βαρυτικό πεδίο είναι η ακτινική παράγωγος του δυναμικού:
\[ g(r)=-\frac{d\Phi}{dr} \]
\[ g(r)=-\frac{GM}{r^2} \]
Αυτό εξηγεί τη σχέση μεταξύ \(1/r\) και \(1/r^2\): το δυναμικό μιας εντοπισμένης μάζας πέφτει ως \(1/r\), ενώ η δύναμη ή η επιτάχυνση πέφτει ως \(1/r^2\).
7. Εξίσωση Poisson
Η πυκνότητα μάζας και το βαρυτικό δυναμικό συνδέονται μέσω της εξίσωσης του Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Σε σφαιρική συμμετρία, αυτό γίνεται:
\[ \frac{1}{r^2}\frac{d}{dr} \left( r^2\frac{d\Phi}{dr} \right) = 4\pi G\rho(r) \]
Η εξίσωση αυτή συνδέει τρεις ποσότητες:
\[ \rho(r) \longrightarrow M(r) \longrightarrow \Phi(r) \longrightarrow v(r) \]
8. Δυνατότητα μιας εκτεταμένης 3D πυκνότητας
Για μια γενική τρισδιάστατη κατανομή πυκνότητας, το δυναμικό είναι:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} \,d^3x’ \]
Ο πυρήνας \(1/|\mathbf{x}-\mathbf{x}’|\) είναι η μαθηματική προέλευση του δυναμικού \(1/r\) στις τρεις διαστάσεις.
9. Δυναμικό ενός λεπτού δίσκου
Για ένα λεπτό δίσκο με επιφανειακή πυκνότητα \(\Sigma(R’)\), το βαρυτικό δυναμικό στο επίπεδο του δίσκου μπορεί να γραφεί ως εξής:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’\,dR’\,d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Η απόσταση μεταξύ ενός σημείου πεδίου σε ακτίνα \(R\) και ενός σημείου πηγής σε ακτίνα \(R’\) είναι:
\[ d(R,R’,\phi) = \sqrt{R^2+R’^2-2RR’\cos\phi} \]
Η ακτινική επιτάχυνση στο δίσκο προκύπτει από τη διαφοροποίηση του δυναμικού:
\[ g_R(R)=-\frac{\partial \Phi}{\partial R} \]
Η ταχύτητα περιστροφής προκύπτει από:
\[ v^2(R)=R\,|g_R(R)| \]
10. Προβολή μιας τρισδιάστατης αλληλεπίδρασης στο δίσκο
Εάν μια αλληλεπίδραση διαδίδεται σε τρεις διαστάσεις αλλά αξιολογείται στο επίπεδο του δίσκου, η ακτινική προβολή εισάγει έναν γεωμετρικό παράγοντα. Για δύο σημεία στο δίσκο που απέχουν απόσταση \(d\), ο παράγοντας ακτινικής προβολής είναι:
\[ \cos\theta = \frac{R-R’\cos\phi}{d(R,R’,\phi)} \]
Έτσι, ένας γενικός τρισδιάστατος ακτινικός πυρήνας \(K(d)\), που προβάλλεται στο δίσκο, εμφανίζεται ως εξής:
\[ K_{\\rm disk}(R,R’,\phi) = K(d) \frac{R-R’\cos\phi}{d} \]
Για παράδειγμα, ένας πυρήνας που μοιάζει με Νευτώνεια δύναμη έχει:
\[ K(d)\propto \frac{1}{d^2} \]
\[ K_{\rm disk}\propto \frac{R-R’\cos\phi}{d^3} \]
Ένας εκθετικός τρισδιάστατος πυρήνας μπορεί να γραφτεί ως εξής:
\[ K(d)\propto e^{-d/\lambda} \]
\[ K_{\rm disk}\propto e^{-d/\lambda} \frac{R-R’\cos\phi}{d} \]
11. Εκθετικοί πυρήνες σε τρεις διαστάσεις
Ένας καθαρά εκθετικός ακτινικός παράγοντας έχει τη μορφή:
\[ e^{-r/\lambda} \]
Στην τρισδιάστατη θεωρία πεδίου, ένα εκθετικά θωρακισμένο δυναμικό εμφανίζεται συχνά σε μορφή Yukawa-like:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Η σχετική ακτινική δύναμη περιέχει τόσο \(1/r^2\) όσο και εκθετικούς όρους:
\[ g_Y(r) = -\frac{d}{dr} \left( \frac{e^{-r/\lambda}}{r} \right) \]
\[ g_Y(r)\propto e^{-r/\lambda} \left( \frac{1}{r^2} + \frac{1}{\lambda r} \right) \]
Αυτό δείχνει γιατί η εκθετική ακτινική συμπεριφορά σε 3D δεν είναι ανεξάρτητη από τη γεωμετρία \(1/r\). Η εκθετική ελέγχει την εξασθένηση, ενώ τα \(1/r\) και \(1/r^2\) προκύπτουν από την τρισδιάστατη εξάπλωση.
12. Νόμοι ακτινικής κλιμάκωσης
Το πρόβλημα της ελλείπουσας μάζας συνδέεται στενά με την ακτινική κλιμάκωση. Αρκετοί σημαντικοί ακτινικοί νόμοι εμφανίζονται επανειλημμένα:
| Ποσότητα | Τυπική κλιμάκωση | Σημασία |
|---|---|---|
| Δυναμικό σημειακής μάζας | \(\Phi(r)\sim 1/r\) | Τρισδιάστατη πράσινη συνάρτηση της βαρύτητας |
| Δύναμη σημειακής μάζας | \(g(r)\sim 1/r^2\) | Παράγωγος της \(1/r\) |
| Επίπεδη ταχύτητα περιστροφής | \(v(r)\sim constant\) | Παρατηρείται στους εξωτερικούς γαλαξιακούς δίσκους |
| Δυναμική μάζα | \(M(r)\sim r\) | Απαιτείται από την επίπεδη περιστροφή |
| Πυκνότητα φωτοστέφανου | \(\rho(r)\sim 1/r^2\) | Δίνει \(M(r)\sim r\) |
| Εκθετικός δίσκος | \(\Sigma(R)\sim e^{-R/R_d}\) | Ο ορατός δίσκος εξασθενεί γρήγορα |
| Διαλεγμένο τρισδιάστατο δυναμικό | \(\Phi(r)\sim e^{-r/\lambda}/r\) | Εκθετική εξασθένηση συν 3D εξάπλωση |
13. Από την πυκνότητα στην καμπύλη περιστροφής
Για ένα σφαιρικό φωτοστέφανο με πυκνότητα:
\[ \rho(r)\propto \frac{1}{r^2} \]
η περικλειόμενη μάζα είναι:
\[ M(r)=4\pi\int_0^r \rho(r’)r’^2dr’ \]
\[ M(r)\propto r \]
Τότε:
\[ v^2(r)=\frac{GM(r)}{r} \]
\[ v^2(r)\approx constant \]
\[ v(r)\approx constant \]
Αυτή είναι η μαθηματική γέφυρα μεταξύ μιας πυκνότητας \(1/r^2\) της άλω και μιας επίπεδης γαλαξιακής καμπύλης περιστροφής.
14. Από τη μάζα του δίσκου στη μάζα που λείπει
Η μάζα του ορατού δίσκου αυξάνεται γρήγορα στην αρχή και στη συνέχεια κορεστεί:
\[ M_{\\rm disk}(R)\rightarrow M_d \]
Η δυναμική μάζα που προκύπτει από μια επίπεδη καμπύλη περιστροφής συνεχίζει να αυξάνεται:
\[ M_{\\rm dyn}(r)\propto r \]
Επομένως, η ελλείπουσα μάζα συμπεριφέρεται περίπου ως εξής:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm disk}(r) \]
\[ M_{\rm miss}(r)\approx \frac{v_0^2}{G}r-M_d \]
Σε αρκετά μεγάλη ακτίνα:
\[ M_{\\rm miss}(r)\sim r \]
\[ \rho_{\rm miss}(r)\sim \frac{1}{r^2} \]
15. Μαθηματική προειδοποίηση: ο δίσκος και η σφαίρα δεν είναι εναλλάξιμες.
Η εξίσωση
\[ M(r)=\frac{v(r)^2r}{G} \]
είναι ακριβές για σφαιρική συμμετρία. Για έναν πεπλατυσμένο δίσκο, θα πρέπει να υπολογιστεί το δυναμικό με ολοκλήρωση πάνω στο δίσκο και στη συνέχεια να εξαχθεί η ακτινική επιτάχυνση. Η σφαιρική έκφραση χρησιμοποιείται συχνά ως αποτελεσματική προσέγγιση, ειδικά όταν συζητείται η μάζα που απαιτείται για την υποστήριξη μιας δεδομένης καμπύλης περιστροφής.
16. Σύνοψη βασικών εξισώσεων
Πυκνότητα επιφάνειας δίσκου:
\[ \Sigma(R)=\Sigma_0e^{-R/R_d} \]
Στοιχείο μάζας δίσκου:
\[ dM_{\rm disk}=2\pi R\Sigma(R)dR \]
Μάζα ορατού δίσκου:
\[ M_{\rm disk}(R)=2\pi\Sigma_0R_d^2 \left[ 1-e^{-R/R_d} \left( 1+\\frac{R}{R_d} \right) \right] \]
Σφαιρική μάζα όγκου:
\[ M(r)=4\pi\int_0^r\rho(r’)r’^2dr’ \]
Πυκνότητα από την περικλειόμενη μάζα:
\[ \rho(r)=\frac{1}{4\pi r^2}\frac{dM}{dr} \]
Δυναμική μάζα:
\[ M_{\rm dyn}(r)=\frac{v(r)^2r}{G} \]
Λείπουσα μάζα:
\[ M_{\rm miss}(r)=M_{\rm dyn}(r)-M_{\rm vis}(r) \]
Νευτώνειο δυναμικό:
\[ \Phi(r)=-\frac{GM}{r} \]
Εξίσωση Poisson:
\[ \nabla^2\Phi=4\pi G\rho \]
Ολοκληρωμένο τρισδιάστατο δυναμικό:
\[ \Phi(\mathbf{x}) = -G\int \frac{\rho(\mathbf{x}’)}{|\mathbf{x}-\mathbf{x}’|} d^3x’ \]
Δυναμικό λεπτού δίσκου:
\[ \Phi(R) = -G \int_0^\infty \int_0^{2\pi} \frac{\Sigma(R’)R’dR’d\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Προβαλλόμενος ακτινικός πυρήνας:
\[ \cos\theta= \frac{R-R’\cos\phi} {\sqrt{R^2+R’^2-2RR’\cos\phi}} \]
Εκθετικό τρισδιάστατο διαλεγμένο δυναμικό:
\[ \Phi_Y(r)\propto -\frac{e^{-r/\lambda}}{r} \]
Συμπέρασμα
Το πρόβλημα της ελλείπουσας μάζας είναι μια μαθηματική αναντιστοιχία μεταξύ δύο ακτινικών συμπεριφορών. Ο ορατός δίσκος ακολουθεί μια εκθετική επιφανειακή πυκνότητα και φτάνει σε μια πεπερασμένη αθροιστική μάζα. Η δυναμική μάζα που συνάγεται από τις περίπου επίπεδες καμπύλες περιστροφής αυξάνεται περίπου γραμμικά με την ακτίνα. Αν ερμηνευθεί ως σφαιρική άλω, αυτό αντιστοιχεί σε μια πυκνότητα που μειώνεται περίπου ως \(1/r^2\). Οι εξισώσεις της ολοκλήρωσης του δίσκου, των σφαιρικών κελυφών, του βαρυτικού δυναμικού και της ακτινικής προβολής παρέχουν τη μαθηματική γλώσσα που απαιτείται για την ανάλυση αυτής της αναντιστοιχίας.