BeeTheory – Gravitaatio- ja aaltofysiikka
Linnunradan piilomassan säteittäiset yhtälöt
Tiheysprofiileista rengasintegraaleihin ja kiertokäyriin – piilomassan matemaattinen käsittely galaktisen säteen R funktiona.
Tällä sivulla esitellään säteittäiset yhtälöt, joita käytetään Linnunradan piilomassan kuvaamiseen. Siinä verrataan klassisia pimeän aineen tiheysprofiileja, rengas- ja kuoriintegraaleja, suljetun massan yhtälöitä, kiertokäyriä ja BeeTeorian tulkintaa puuttuvasta massasta mahdollisena aalto-interferenssi-ilmiönä.
~10¹² M⊙
Klassinen arvio Linnunradan halon kokonaismassasta.
0,39 GeV/cm³
Tyypillinen paikallinen pimeän aineen tiheys lähellä Aurinkoa.
R⊙ ≈ 8 kpc R⊙ ≈ 8 kpc
Auringon likimääräinen etäisyys galaktisesta keskuksesta.
~200 kpc
Linnunradan haloarvioiden arvioinnissa käytetty likimääräinen ulompi asteikko.
Sisältö
- Miksi massa puuttuu?
- Tiheysprofiilit ρ(R)
- Rengas- ja kehämassa dM/dR
- Suljetun pimeän aineen massa M(<R)
- Kiertokäyrä V(R)
- Nykyiset havainnolliset arviot
- Kilpailevat hypoteesit
- BeeTheory-näkökulma
1. Miksi massa puuttuu?
Vuonna 1933 Fritz Zwicky huomasi, että Coma-joukon galaksit liikkuivat aivan liian nopeasti pysyäkseen koossa pelkästään näkyvän massansa avulla. 1970-luvulla Vera Rubin ja Kent Ford mittasivat spiraaligalaksien kiertokäyriä ja havaitsivat jotain yhtä hämmästyttävää: suurilla säteillä olevat tähdet kiertävät lähes yhtä nopeasti kuin lähellä keskustaa olevat tähdet, vaikka näkyvän aineen Newtonin painovoiman mukaan niiden pitäisi hidastua.
Yksinkertaisen kepleriläisen radan, joka kiertää keskusmassaa, odotamme:
\(V(R)\propto \frac{1}{\sqrt{R}}\)Sen sijaan havaitaan suunnilleen tasainen tai vain hitaasti laskeva kiertokäyrä:
\(V(R)\approx V_{\infty}=\mathrm{const}\qquad \mathrm{for}\ R\gtrsim 5\,\mathrm{kpc}\)Näiden tosiseikkojen sovittaminen yhteen Newtonin painovoiman kanssa edellyttää ylimääräistä näkymätöntä massakomponenttia, jonka tiheys laskee suunnilleen seuraavasti:
\(\rho(r)\propto r^{-2}\)Näin saadaan suljetun kappaleen kokonaismassa suhteessa säteeseen:
\(M(<R)\propto R\)ja siksi:
\(V\propto \sqrt{\frac{M}{R}}=\mathrm{const}\)Keskeinen määrällinen pulma
Linnunradan valoisa baryoninen massa on noin 5 × 10¹⁰ M⊙. Kinematiikan perusteella päätelty dynaaminen kokonaismassa noin 200 kpc:n etäisyydelle on noin 10¹² M⊙. Tämä merkitsee, että pimeän ja valon massasuhde on noin 10:1.
2. Tiheysprofiilit ρ(R)
Tiheysprofiili on matemaattinen funktio, joka kuvaa, miten pimeän aineen tiheys ρ vaihtelee galaktosentrisen säteen r tai galaktisessa tasossa olevan sylinterin säteen R mukaan.
2.1 NFW-profiili
Navarron, Frenkin ja Whiten esittämä NFW-profiili on johdettu N-kappaleen kosmologisista simulaatioista. Se tuottaa tyypillisen kaksinkertaisen potenssilain, jossa on keskeinen kärki.
\(\rho_{\mathrm{NFW}}(r)=\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}\)| Parametri | Symboli | Linnunradan arvio | Rooli |
|---|---|---|---|
| Mittakaavan säde | rs | 15-25 kpc | Sisä- ja ulkorinteiden välinen siirtymä |
| Ominaistiheys | ρ0 | Kalibroitu paikalliseen pimeän aineen tiheyteen | Yleinen normalisointi |
| Sisäinen rinne | γ | -1 | Nöyrä käytös |
| Ulompi rinne | – | -3 | Nopea lasku suurella säteellä |
2.2 Einaston profiili
Einasto-profiilissa vältetään tiukka keskipisteen divergenssi ja käytetään muotoparametria α, jonka avulla tiheyden kaltevuus muuttuu tasaisesti säteen myötä.
\(\rho_{\mathrm{Ein}}(r)=\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)| Parametri | Symboli | Linnunradan arvio | Rooli |
|---|---|---|---|
| Muotoindeksi | α | Mallista riippuvainen | Säätää kuinka nopeasti kaltevuus muuttuu |
| Mittakaavan säde | r-2 | ~18-22 kpc | Säde, jonka logaritminen kaltevuus on -2. |
| Tiheys r-2:ssa | ρ-2 | Kalibroitu paikallisen tiheyden mukaan | Normalisointi |
Viimeaikaiset havainnolliset jännitteet
Viimeaikaiset Gaia-pohjaiset tutkimukset viittaavat siihen, että Linnunradan kiertokäyrä saattaa laskea nopeammin Auringon säteen ulkopuolella kuin tavanomainen NFW-halo ennustaa. Tämä tekee Einaston kaltaisista ydin- tai tasaisesti vaihtelevista profiileista erityisen tärkeitä nykyisissä keskusteluissa.
2.3 Pseudo-isoterminen profiili
Pseudo-isotermistä profiilia käytetään usein yksinkertaisena analyyttisenä approksimaationa ydinhalolle.
\(\rho_{\mathrm{iso}}(r)=\frac{\rho_0}{1+\left(\frac{r}{r_s}\right)^2}\)Pienellä säteellä tiheys lähestyy vakioarvoa. Suurella säteellä se laskee r-²:n suuruisena ja tuottaa litteän kiertokäyrän.
\(V_{\infty}=\sqrt{4\pi G\rho_0 r_s^2}\)Huipun ja ytimen välinen ongelma
N-kappale-simulaatiot ennustavat usein kuoppaisia NFW-profiileja, kun taas monissa havaituissa kääpiögalakseissa tiheysprofiilit näyttävät olevan mieluummin sydänmaisia. Tämä kärki-ydin -ongelma on edelleen yksi tärkeimmistä pimeän aineen fysiikan ratkaisemattomista kysymyksistä.
3. Rengas- ja rengaskehän massa – dM/dR
Jotta voimme laskea, kuinka paljon pimeää ainetta on galaksin kussakin säteittäisessä viipaleessa, integroimme tiheyden ohuen kuoren tai kehän yli. Geometria riippuu siitä, käsitelläänkö haloa pallomaisena vai litteänä.
3.1 Pallomainen ohut kuori
Pallosymmetrisen halon massa kuoren paksuudessa dr säteen r kohdalla on:
\(\frac{dM_{\mathrm{shell}}}{dr}=4\pi r^2\rho(r)\)3.2 Levytason rengasmainen rengas
Galaksin tasossa sijaitsevalle renkaalle, jonka sylinterin säde on R ja efektiivinen puolipituus H(R), rengasmassa on:
\(dM_{\mathrm{ann}}=2\pi R\,\rho(R,0)\,2H(R)\,dR\)Pallomaisen halon tapauksessa tämä voidaan kirjoittaa integraalina korkeuden z yli:
\(\frac{dM}{dR}=2\pi R\int_{-\infty}^{+\infty}\rho\left(\sqrt{R^2+z^2}\right)dz\)Pallomaisessa approksimaatiossa tämä liittyy takaisin:
\(\frac{dM}{dR}\approx4\pi R^2\rho(R)\)3.3 NFW:n massa kuorta kohti
\(\frac{dM_{\mathrm{NFW}}}{dr}=4\pi r^2\frac{\rho_0}{\left(\frac{r}{r_s}\right)\left(1+\frac{r}{r_s}\right)^2}=\frac{4\pi\rho_0 r_s r}{\left(1+\frac{r}{r_s}\right)^2}\)Tämä funktio saavuttaa huippunsa asteikon säteenrs ympärillä, mikä tarkoittaa, että suuri osa pimeän aineen massasta kuorta kohti on kerrostunut pikemminkin välihaloon kuin vain sen keskelle tai laitamille.
3.4 Einaston massa kuorta kohti
\(\frac{dM_{\mathrm{Ein}}}{dr}=4\pi r^2\rho_{-2}\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}\)Einaston koteloitu massa arvioidaan yleensä numeerisesti.
Fyysinen merkitys
Funktio dM/dr kertoo, mikä galaktinen säde vaikuttaa eniten piilomassan budjettiin. Jyrkempi ulompi profiili vähentää halon kokonaismassaa, kun taas matalampi profiili lisää sitä.
4. Suljetun pimeän aineen massa M(
Integroimalla kuorielementti 0:sta R:ään saadaan pimeän aineen kokonaismassa säteen R sisällä:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)
4.1 NFW Suljettu massa
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)
4.2 Einasto Suljettu massa
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)
4.3 Kokonaismassan hajoaminen
Suljettu dynaaminen kokonaismassa voidaan jakaa näkyviin ja piilotettuihin komponentteihin:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)
~3 × 10¹⁰ M⊙
Arvioitu kätketty massa aurinkokehän sisällä.
~1-2 × 10¹¹ M⊙
Arvioitu kätketty massa 20 kpc:n sisällä.
5-12 × 10¹¹ M⊙
Likimääräinen piilotettu massa viriaalialueen sisällä, vahvasti mallista riippuvainen.
Massaprofiili riippuu edelleen mallista.
Linnunradan halon massa-arviot riippuvat suuresti siitä, miten ulompi halo ekstrapoloidaan sen alueen ulkopuolelle, jolla on vahvat havaintorajoitukset.
Integroimalla kuorielementti 0:sta R:ään saadaan pimeän aineen kokonaismassa säteen R sisällä:
\(M_{\mathrm{DM}}(<R)=\int_0^R4\pi r^2\rho(r)\,dr\)4.1 NFW Suljettu massa
\(M_{\mathrm{NFW}}(<R)=4\pi\rho_0r_s^3\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]\)4.2 Einasto Suljettu massa
\(M_{\mathrm{Ein}}(<R)=4\pi\rho_{-2}\int_0^R r^2\exp\left\{-\frac{2}{\alpha}\left[\left(\frac{r}{r_{-2}}\right)^\alpha-1\right]\right\}dr\)4.3 Kokonaismassan hajoaminen
Suljettu dynaaminen kokonaismassa voidaan jakaa näkyviin ja piilotettuihin komponentteihin:
\(M_{\mathrm{tot}}(<R)=M_{\mathrm{bulge}}(<R)+M_{\mathrm{disk}}(<R)+M_{\mathrm{DM}}(<R)\)~3 × 10¹⁰ M⊙
Arvioitu kätketty massa aurinkokehän sisällä.
~1-2 × 10¹¹ M⊙
Arvioitu kätketty massa 20 kpc:n sisällä.
5-12 × 10¹¹ M⊙
Likimääräinen piilotettu massa viriaalialueen sisällä, vahvasti mallista riippuvainen.
Massaprofiili riippuu edelleen mallista.
Linnunradan halon massa-arviot riippuvat suuresti siitä, miten ulompi halo ekstrapoloidaan sen alueen ulkopuolelle, jolla on vahvat havaintorajoitukset.
5. Kiertokäyrä V(R)
Kokonaismassa määrää ympyrän nopeuden säteellä R, joka määräytyy painovoiman ja keskipakokiihtyvyyden tasapainon avulla:
\(V_c(R)=\sqrt{\frac{G\,M_{\mathrm{tot}}(<R)}{R}}\)Koska painovoimapotentiaaliin vaikuttavat toisistaan riippumattomat massakomponentit, niiden nopeusosuudet lasketaan usein yhteen kvadratuurisesti:
\(V_c^2(R)=V_{\mathrm{bulge}}^2(R)+V_{\mathrm{thin\,disk}}^2(R)+V_{\mathrm{thick\,disk}}^2(R)+V_{\mathrm{gas}}^2(R)+V_{\mathrm{DM}}^2(R)\)5.1 Baryonisen levyn osuus
Tähtien ohut kiekko noudattaa eksponentiaalista pintatiheysprofiilia:
\(\Sigma(R)=\Sigma_0\exp\left(-\frac{R}{R_d}\right)\)Vastaava ympyränopeus eksponentiaaliselle levylle on:
\(V_{\mathrm{disk}}^2(R)=\frac{2GM_d}{R_d}y^2\left[I_0(y)K_0(y)-I_1(y)K_1(y)\right],\qquad y=\frac{R}{2R_d}\)TässäIn jaKn ovat muunnettuja Besselin funktioita. Tyypilliset Linnunradan ohuen kiekon parametrit ovatRd ≈ 2,6 kpc ja Md ≈ 3,5 × 10¹⁰ M⊙.
5.2 Pimeän aineen osuus
\(V_{\mathrm{DM,NFW}}(R)=\sqrt{\frac{4\pi G\rho_0r_s^3}{R}\left[\ln\left(1+\frac{R}{r_s}\right)-\frac{R/r_s}{1+R/r_s}\right]}\)5.3 Baryoninen Tully-Fisher-suhde
Baryoninen Tully-Fisher-suhde yhdistää galaksin litteän pyörimisnopeuden sen baryoniseen kokonaismassaan:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0,\qquad a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)~230 km/s
Kiertonopeus lähellä auringon sädettä.
~170-180 km/s
Mahdollinen laskeva arvo ulommassa levyssä, riippuen merkkiainetiedoista.
~150 km/s
Lähtökohtainen ulomman halon nopeusasteikko halon merkkiaineista.
6. Nykyiset havainnolliset arviot
Alla olevassa taulukossa on yhteenveto pimeän aineen tiheyden ja massan edustavista arvoista keskeisillä galaktisilla säteillä. Tarkat arvot vaihtelevat tietokokonaisuuden, merkkiainepopulaation ja halomallin mukaan.
| Säde R | Pimeän aineen tiheys | Suljettu tumma massa | Menetelmä |
|---|---|---|---|
| Keskusta | Eroava NFW:ssä, äärellinen ydinmalleissa. | Mallista riippuvainen | N-kappale-simulaatiot ja galaksin sisäisen galaksin mallintaminen |
| R⊙ ≈ 8 kpc | ~0.39 GeV/cm³ | ~3 × 10¹⁰ M⊙ | Kiertokäyrä ja pystysuora kinematiikka |
| 20 kpc | ~0.05 GeV/cm³ | ~1-2 × 10¹¹ M⊙ | Gaia ja spektroskooppiset merkkiaineet |
| 50 kpc | ~5 × 10-³ GeV/cm³ | ~3-5 × 10¹¹ M⊙ | Palloparvet ja halotähdet |
| 100-200 kpc | ≤10-³ GeV/cm³ | ~5-12 × 10¹¹ M⊙ | Satelliittigalaksit ja pakenemisnopeusmenetelmät |
Palloparvien kinematiikan, halotähtien, satelliittigalaksien ja Gaia-ajan astrometrian yhdistäminen osoittaa, että Linnunradan ulomman halon profiili on edelleen epävarma. Tämä epävarmuus on keskeistä piilomassaongelman kannalta.
7. Kilpailevat hypoteesit puuttuvasta massasta
Useat suuret selitysperheet ovat edelleen aktiivisia. Mitään niistä ei ole lopullisesti vahvistettu tai suljettu pois kaikilla havaintoskaaloilla.
7.1 Kylmän pimeän aineen hiukkaset
Kylmä pimeä aine on edelleen johtava paradigma. Ehdokkaita hiukkasia ovat WIMPit, steriilit neutriinot ja muut standardimallin ulkopuoliset mahdollisuudet. Nämä ehdokkaat muodostavat laajoja haloja, joita usein mallinnetaan NFW- tai Einasto-profiileilla.
\(m_{\chi}\sim10\text{–}1000\,\mathrm{GeV}\)Suurin jännite on kokeellinen: suoralla havaitsemisella ei ole vielä löydetty vahvistettua pimeän aineen hiukkasta.
7.2 Ultrakevyt tai sumea pimeä aine (Fuzzy Dark Matter)
Sumeassa pimeässä aineessa käytetään ultrakevyitä hiukkasia, joiden de Broglie -aallonpituus voi kasvaa astrofyysisesti suureksi, jolloin pienen mittakaavan rakenne tukahdutetaan.
\(m_a\sim10^{-22}\,\mathrm{eV}\) \(\lambda_{\mathrm{dB}}\sim\mathrm{kpc}\)Tämä kehys tuottaa luonnollisesti tasaisempia sisäisiä tiheysytimiä, mutta sitä rajoittavat Lyman-alfa-metsätiedot ja kääpiögalaksien rakenne.
7.3 Modifioitu newtonilainen dynamiikka
MOND muuttaa efektiivistä gravitaatiokiihtyvyyttä tietyn asteikon alapuolella:
\(a_0\approx1.2\times10^{-10}\,\mathrm{m/s^2}\)Syvässä MOND-järjestelmässä tehollinen kiihtyvyys on:
\(g_{\mathrm{eff}}=\sqrt{g_Na_0}\)MOND ennustaa baryonisen Tully-Fisher-suhteen:
\(V_{\infty}^4=G\,M_{\mathrm{bar}}\,a_0\)Se toimii hyvin monille galaksien pyörimisliikkeille, mutta galaksijoukkojen ja kosmologian tutkiminen on edelleen vaikeaa.
7.4 Itse vuorovaikutuksessa oleva pimeä aine
Itsevuorovaikutteisen pimeän aineen mukaan pimeän aineen hiukkaset vuorovaikuttavat toistensa kanssa niin voimakkaasti, että ne muokkaavat halon sisäisiä tiheysprofiileja.
\(\frac{\sigma}{m}\sim1\text{–}100\,\mathrm{cm^2/g}\)Tämä saattaa auttaa selittämään halosydänten monimuotoisuuden, mutta mitään tiettyä hiukkasehdokasta ei ole vielä vahvistettu.
7.5 Alkuaikojen mustat aukot
Varhaisessa maailmankaikkeudessa syntyneet alkumustat aukot voisivat muodostaa osan piilossa olevasta massasta. Monet massaikkunat ovat vahvasti rajattuja mikrolisenssi-, kosmisen mikroaaltotaustan ja gravitaatioaaltohavaintojen perusteella.
\(10^{-16}\text{–}10^{-11}\,M_\odot\)Ne ovat edelleen spekulatiivisia, koska ne eivät ole täydellinen selitys Linnunradan piilomassalle.
8. Mehiläisteorian näkökulma
Mehiläisteoriassa ehdotetaan, että painovoima voidaan ymmärtää aaltokäyttäytymisestä johtuvana emergenttinä vaikutuksena eikä niinkään perustavanlaatuisena voimana, jota vain hiukkanen kantaa tai jonka tuottaa vain avaruusajan kaarevuus.
Tässä kehyksessä jokaiseen massiiviseen systeemiin liittyy aaltofunktio ψ(r,t). Kvanttitutkimuksen peruslähtökohta on kolmiulotteinen Schrödingerin yhtälö:
\(i\hbar\frac{\partial\psi}{\partial t}=-\frac{\hbar^2}{2m}\nabla^2\psi+V(\mathbf{r})\psi\)Kun kaksi massajakaumaa lähestyy toisiaan, niiden aaltofunktiot menevät päällekkäin. Näiden aaltofunktioiden konvoluutio voidaan kirjoittaa seuraavasti:
\(\Psi_{12}(\mathbf{r})=(\psi_1*\psi_2)(\mathbf{r})=\int\psi_1(\mathbf{r}’)\psi_2(\mathbf{r}-\mathbf{r}’)\,d^3r’\)BeeTeoriassa gravitaatiovetovoima tulkitaan strukturoitujen aaltojen päällekkäisyyden, resonanssin ja kenttäkoherenssin laajamittaiseksi ilmentymäksi.
8.1 Mehiläisteorian uudelleentulkinta piilomassasta
BeeTeoriassa se, mitä yleensä kutsutaan pimeäksi aineeksi, voidaan tulkita monien galaktisen halon alueelle jakautuneiden värähtelevien järjestelmien aiheuttaman aaltojen interferenssin kumulatiiviseksi gravitaatiovaikutukseksi.
\(\rho_{\mathrm{eff}}(R)=\rho_{\mathrm{bar}}(R)+\Delta\rho_{\mathrm{wave}}(R)\)Tässä Δρwave(R) edustaa ylimääräistä tehokasta gravitaatiotiheyttä, joka johtuu pikemminkin koherentista aaltokentän rakenteesta kuin suoraan näkyvästä baryonisesta aineesta.
Tämän termin olisi toistettava säteittäinen käyttäytyminen, joka tavallisesti liitetään pimeään aineeseen. Erityisesti sen olisi tuotettava suunnilleen tasaiset pyörimisnopeuskäyrät galaktisella alueella.
\(\rho_{\mathrm{wave}}(R)\propto R^{-2}\)Avoin määrällinen haaste
BeeTeorian on osoitettava, voiko aaltopohjainen interferenssimalli toistaa havaittujen pyörimiskäyrien edellyttämän tarkan säteittäisen tiheysprofiilin. Sen on myös selitettävä, miksi efektiivinen piilomassa on usein paljon suurempi kuin näkyvä baryoninen massa.
Viitteet
- Ou, X., Eilers, A.-C., Necib, L., Frebel, A. – The dark matter profile of the Milky Way inferred from its circular velocity curve (Linnunradan pimeän aineen profiili sen ympyränopeuskäyrästä pääteltynä), MNRAS 528, 693-710, 2024.
- Navarro, J. F., Frenk, C. S. S., White, S. D. M. – A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493, 1997.
- Einasto, J. – Galaksin komposiittimallin rakentamisesta, Trudy 5, 87, 1965.
- Watkins, L. L., van der Marel, R. P. et al. – Evidence for an Anticorrelation between the Masses of the Milky Way and Andromeda Galaxies, ApJ 873, 111, 2019.
- Milgrom, M. – A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis, ApJ 270, 365, 1983.
- McGaugh, S. S. et al. – Radial Acceleration Relation in Rotationally Supported Galaxies, PRL 117, 201101, 2016.
Huomautus: tuoreet tai tulevaisuudessa päivitettävät viittaukset on tarkistettava ennen julkaisemista, jos sivua käytetään tieteellisenä lähteenä.
BeeTheory näkökulma
Piilomassaongelma ei ole vain kysymys siitä, kuinka paljon ainetta puuttuu. Kyse on siitä, millainen fyysinen rakenne tuottaa painovoiman galaktisessa mittakaavassa.
Klassiset pimeän aineen mallit tulkitsevat puuttuvan massan näkymättömäksi aineeksi. BeeTheory tutkii täydentävää mahdollisuutta: osa piilevästä gravitaatiovaikutuksesta voi syntyä rakenteellisesta aaltokoherenssista.
Seuraava vaihe on matemaattinen: määritetään radiaalinen aaltotiheystermi, johdetaan sen pyörimisliike ja verrataan sitä suoraan Gaia-ajan Linnunradan tietoihin.