Podsumowanie
Istnieje szansa 1 na 500 miliardów, że stoisz pod gwiazdą poza Drogą Mleczną, szansa 1 na 3,3 miliarda, że stoisz pod gwiazdą Drogi Mlecznej, i szansa 1 na 184 tysięcy, że stoisz pod Słońcem, w prawo teraz.
Duży, gruby, śmierdzący, Ostrzeżenie! Starałem się, aby matematyka była prosta, ale to wszystko, co właśnie wymyśliłem. Nie udzielam żadnych gwarancji, że jest całkowicie dokładny, ale liczby wydają się pozytywnie sprawdzone pod względem zdrowia psychicznego, więc myślę, że jesteśmy dobrzy.
Zastrzegamy pierwszy : Liczby gwiazd innych niż Słońce oparte są na danych z dużą niepewnością, takich jak liczba gwiazd we wszechświecie i średni rozmiar gwiazdy. Powyższe liczby mogą być z łatwością pomniejszone o współczynnik 10 w obu kierunkach i mają jedynie na celu przybliżenie, jak pusta jest przestrzeń.
Zastrzegaj drugi : liczby dla Słońca i Drogi Mlecznej oparte są na założeniu, że stoisz (lub unosisz się) w losowym punkcie na Ziemi. Nikt poza tropikami nigdy nie będzie miał Słońca nad głową. Ludzie na półkuli północnej częściej mają gwiazdy Drogi Mlecznej nad głową, przy czym największe szanse to ludzie w pobliżu 36,8 ° N, ponieważ na tej szerokości geograficznej prosto przez centrum galaktyki przechodzi raz dziennie. 26
Uwaga : Możesz w większości zignorować wszystko w tej odpowiedzi i po prostu spojrzeć na stały kąt Słońca, aby uzyskać ten sam wynik. Wszystkie pozostałe gwiazdy są naprawdę daleko i bardzo rozproszone. Różnica w podanym kącie bryłowym wynosi pięć tysięcznych procenta więcej, gdy dodamy resztę wszechświata do Słońca.
tło
Spróbujmy uzyskać nieco realistyczną, twardą liczbę. Aby to zrobić, potrzebujemy pewnych założeń.
Jak wskazano w odpowiedzi Michaela Walsby'ego 1 , jeśli wszechświat jest nieskończony (i jednorodny 2 ), istnieje tylko nieskończenie mała szansa, że nie będzie gwiazdy nad głową, co normalna matematyka traktuje jako dokładnie zerową szansę. Załóżmy więc, że wszechświat jest skończony.
Domniemania
- W szczególności załóżmy, że wszechświat składa się tylko z obserwowalnego wszechświata. (Więcej informacji na temat ekspansji wszechświata 3 ).
- Ponadto załóżmy, że zawartość obserwowalnego wszechświata jest mierzona w ich bieżących (zakładanych) pozycjach, a nie w pozycji, w jakiej się wydają. (Jeśli zobaczymy światło gwiazdy z 400 milionów lat po rozpoczęciu wszechświata, zmierzymy ją jako oddaloną o około 13,5 miliarda lat świetlnych, ale obliczamy, że z powodu ekspansji jest ona prawdopodobnie bliższa 45 miliardów lat świetlnych.)
- Przyjmiemy liczbę gwiazd we obserwowalnym wszechświecie jako . 2013 oszacowanie 4 było , A 2014 oszacowanie 5 było , a 2017 oszacowanie 6 było , z każdego artykułu spodziewa oszacowanie zwiększyć, ponieważ mamy lepsze teleskopy w czasie. Więc weźmiemy najwyższą wartość i wykorzystamy ją.1024 10 21 10 23 10 24102110231024
- Weźmiemy wielkości obserwowalnego świata 7 się , dając powierzchnię 8 o 9 i tom 10 od 11 .8.8⋅1026m (diameter)2,433 ⋅ 10 54 m 2 3,568 ⋅ 10 80 m 32.433⋅1054m2 3.568⋅1080m3
- Przyjmiemy, że średni rozmiar gwiazdy jest wielkością Słońca, 12 . (Nie mogę znaleźć żadnych źródeł średniej wielkości gwiazd, tylko to, że Słońce jest gwiazdą średnią).1.4⋅109m (diameter)
Model
Odtąd będziemy oszukiwać. Realistycznie powinniśmy modelować każdą galaktykę osobno. Ale zamierzamy udawać, że cały wszechświat jest idealnie jednolity (jest to wystarczające, gdy oddalamy się od Ziemi w wielkim schemacie kosmosu). Co więcej, zaczniemy liczyć na tyle daleko, aby całkowicie zignorować Drogę Mleczną i Słońce, a następnie dodać je później z różnymi obliczeniami.
Biorąc pod uwagę powyższe domniemania, możemy łatwo obliczyć gęstość gwiazdową obserwowalnego wszechświata, tak aby był 13 .δ=1024stars3.568⋅1080m3=2.803⋅10−57starsm3
Następnie musimy obliczyć kąt bryłowy 14 w oparciu o gwiazdę. Kąt bryłowy kuli podaje 15 , gdzie to kąt bryłowy w steradianach 16 (sr), to odległość od kuli, a to promień kuli. Użycie jako średnicy, która konwertuje na . Biorąc pod uwagę średnią średnicę zakładaną powyżej ( ), daje to średni kąt bryłowyΩ=2π(1−d2−r2√d) sr ΩdrDΩ=2π(1- √ΩdrDΩ=2π⎛⎝⎜1−d2−(D2)2√d⎞⎠⎟ sr1.4⋅109mΩ=2π(1−d2−4.9⋅1017m2√d) sr 17 .
W tym momencie moglibyśmy ustawić właściwą całkę, ale mój rachunek jest raczej zardzewiały i na początek niezbyt ostry. Przybliżę odpowiedź za pomocą serii koncentrycznych powłok, z których każda ma grubość (około miliona lat świetlnych). Odłożymy naszą pierwszą powłokę , a następnie wyruszymy stamtąd.1022m1022m
Obliczymy całkowity kąt bryły każdej skorupy, a następnie dodamy wszystkie skorupy razem, aby uzyskać kąt bryły objęty przez cały obserwowalny wszechświat.
Ostatnim problemem do rozwiązania jest nakładanie się. Niektóre gwiazdy w dalszych powłokach zachodzą na gwiazdy w pobliskich powłokach, co powoduje, że przeceniamy całkowity zasięg. Obliczymy zatem prawdopodobieństwo nakładania się danej gwiazdy i stamtąd zmodyfikujemy wynik.
Zignorujemy nakładanie się powłoki, modelując tak, jakby każda gwiazda w skorupie znajdowała się w stałej odległości, równomiernie rozmieszczonej w skorupie.
Prawdopodobieństwo nakładania się
Aby dana gwiazda zachodziła na gwiazdy bliższe, musi znajdować się w pozycji już pokrytej przez gwiazdy bliższe. Dla naszych celów będziemy traktować nakładanie się jako układ podwójny: albo gwiazda jest całkowicie zachodząca na siebie, albo wcale się nie nakłada.
Prawdopodobieństwo będzie podane przez ilość kąta bryłowego już uwzględnionego przez poprzednie pociski podzieloną przez całkowity kąt bryłowy na niebie ( ).4π sr
Nazwijmy prawdopodobieństwo, że dana gwiazda, , zachodzi na siebie , kąt bryłowy objęty przez tę gwiazdę , a liczba gwiazd . Ilość niezakładającego się kąta bryłowego objętego przez daną powłokę, , wynosi wtedy . Ponieważ powiedzieliśmy, że gwiazdy w powłoce nie nakładają się na siebie, jest takie samo dla wszystkich w danej powłoce, co pozwala nam uprościć powyższe równanie do , gdzieiPiΩinkΩkT=(1−P1)Ω1+(1−P2)Ω2+…+(1−Pn)Ωn srstarPiiΩkT=(1−Pk)(Ω1+Ω2+…+Ωn) srstarPkto prawdopodobieństwo nakładania się powłoki . Ponieważ traktujemy wszystkie gwiazdy jako mające ten sam, średni rozmiar, upraszcza to jeszcze bardziej do , gdzie to kąt bryły gwiazdy w powłoce .kΩkT=(1−Pk)Ωkn srstarΩkk
Obliczanie kąta bryłowego
Liczba gwiazd w skorupie jest podana przez objętość powłoki pomnożoną przez gęstość gwiazdową tej powłoki. W przypadku skorup odległych możemy traktować objętość skorupy jako jej pole powierzchni razy jej grubość. , gdzie to odległość od skorupy, a to jej grubość. Używając jako gęstości gwiazdowej, liczba gwiazd jest po prostu .Vshell=4πd2tdtδn=δVshell=δ4πd2t
Stąd możemy użyć obliczenia kąta bryłowego skorupy (z prawdopodobieństwa nakładania się powyżej), aby uzyskać .ΩkT=(1−Pk)Ωkδ4πd2t srstar
Zauważ, że jest podane jako częściowa suma kąta bryłowego dla wszystkich poprzednich powłok podzielona przez całkowity kąt bryłowy. A podaje (z modelu powyżej).PkΩkΩk=2π(1−d2k−4.9⋅1017m2√dk) srstar
To daje nam . Biorąc pod uwagę, że każda powłoka znajduje się w odległości , możemy zastąpić przez . Podobnie, można zastąpić . I już obliczyliśmy (z modelu powyżej).ΩkT=(1−Ω(k−1)T4π)2π(1−d2k−4.9⋅1017m2√dk)δ4πd2t sr1022mdkk1022mt1022mδ=2.803⋅10−57starsm3
To daje nam
ΩkT=(1−Ω(k−1)T4π)2π(1−(k1022m)2−4.9⋅1017m2√k1022m)2.803⋅10−57starsm34π(k1022m)21022m srstar
=(1−Ω(k−1)T4π)(1−k21044−4.9⋅1017√k1022)2.803⋅10−578π2k21066 sr
=(1−Ω(k−1)T4π)2.213⋅1011k2(1−k21044−4.9⋅1017√k1022) sr
Stąd możemy po prostu podłączyć liczby do programu obliczeniowego.
ΩT=∑kmaxk=1ΩkT
Gdzie to tylko promień obserwowalnego wszechświata podzielony przez grubość danej powłoki. Zatemkmaxkmax=4.4⋅1026m1022m=4.4⋅104=44000
ΩT=∑44000k=1ΩkT
Wyniki
Ze względu na dużą liczbę trudno jest po prostu uruchomić to w programie. Zacząłem pisać niestandardowy program C ++ przy użyciu biblioteki ttmath 18 dla dużych liczb. Wynik to lub całego nieba. Z drugiej strony istnieje około 1 na 500 miliardów szans, że staniesz teraz pod gwiazdą.2.386⋅10−11 sr1.898⋅10−12
Zauważ, że zignorowaliśmy w tym celu Drogę Mleczną i Słońce.
Program C ++ można znaleźć na PasteBin 25 . Będziesz musiał poprawnie działać. Dodałem kilka instrukcji na początku kodu C ++, aby zacząć, jeśli chcesz, aby działał. Nie jest elegancki ani nic, wystarczy, aby funkcjonować.
Słońce
WolframAlpha poinformował mnie, że Słońce ma stały kąt około , czyli około 2,8 miliona razy więcej niż wszystkie gwiazdy wszechświata razem wzięte. Powyższa formuła kąta bryłowego daje tę samą odpowiedź 18, jeśli podamy odległość Słońca wynoszącą 150 gigametrów i promień 0,7 gigametru.6.8⋅10−5 sr
Droga Mleczna
Możemy uzyskać przybliżenie Drogi Mlecznej, biorąc jej rozmiar i gęstość i wykonując te same obliczenia, co powyżej, z wyjątkiem mniejszej skali. Jednak galaktyka jest bardzo płaska, więc szanse w dużym stopniu zależą od tego, czy staniesz w płaszczyźnie galaktycznej, czy nie. Poza tym jesteśmy z boku, więc w kierunku centrum galaktyki jest znacznie więcej gwiazd niż w oddali.
Jeśli przybliżymy galaktykę jako walec o promieniu (około 52000 lat świetlnych) i wysokości (około 2 lata świetlne), otrzymujemy objętość 20 .5⋅1020 m2⋅1016 m1.571⋅1058 m3
Obecne oszacowania promienia galaktyki są bliższe 100 000 lat świetlnych 21 22 , ale przypuszczam, że ogromna większość gwiazd jest znacznie bliżej.
Szacuje się, że w Drodze Mlecznej znajduje się od 100 do 400 miliardów gwiazd 21 . Wybierzmy 200 miliardów na nasze cele. To stawia gęstość Drogi Mlecznej na 22 , czyli około 4,5 miliarda razy gęstsze niż cały wszechświat.δ=200⋅109stars1.571⋅1058 m3=1.273⋅10−47starsm3
Tym razem weźmiemy pociski o grubości (około 10 lat świetlnych) i stamtąd wyjdziemy. Ale musimy ponownie uporządkować matematykę w sferyczną formę, więc założymy, że galaktyka ma taką samą objętość, ale jest kulą. Daje to promień 24 lub 155,4 pocisków. Zaokrąglimy do 155 pocisków.1017 m1.554⋅1019 m
ΩT=∑155k=1ΩkT
Korzystając z powyższej formuły ( Obliczanie kąta bryłowego ), możemy zacząć zastępować liczby.
ΩkT=(1−Ω(k−1)T4π)2π(1−d2k−4.9⋅1017m2√dk)δ4πd2tsrstar
=(1−Ω(k−1)T4π)2π(1−(k⋅1017 m)2−4.9⋅1017 m2√k⋅1017 m)1.273⋅10−47starsm34π(k⋅1017 m)21017 msrstar
=(1−Ω(k−1)T4π)(1−k2⋅1034 m2−4.9⋅1017 m2√k1017 m)1.273⋅10−47starsm38π2k21051 m3srstar
=(1−Ω(k−1)T4π)⋅1.005⋅106k2(1−k2⋅1034−4.9⋅1017√k1017) sr
Podłączenie tego do programu daje , co stanowi całego nieba. Szanse, że stoisz pod gwiazdą w Drodze Mlecznej, wynoszą około 1 na 3,3 miliarda.3.816⋅10−9 sr3.037⋅10−10
Sumy z pełnym kątem
Kąt bryłowy wynosi:
- Sun,6.8⋅10−5 sr
- Droga Mleczna,3.816⋅10−9 sr
- Wszechświat,2.386⋅10−11 sr
- Łącznie, (dodatkowe cyfry są w zasadzie bez znaczenia, dodając około pięciu tysięcznych procent do stałego kąta Słońca) 6.800384⋅10−5 sr
- Droga Mleczna plus Wszechświat, (około 0,6% więcej niż tylko Droga Mleczna)3.840⋅10−9 sr
Bibliografia
1 Odpowiedź Michaela Walsby'ego na to pytanie , czy jest gwiazda nad moją głową? . https://astronomy.stackexchange.com/a/33294/10678
2 Artykuł w Wikipedii , zasada kosmologiczna . https://en.wikipedia.org/wiki/Cosmological_principle
3 Artykuł w Wikipedii , Ekspansja wszechświata . https://en.wikipedia.org/wiki/Expansion_of_the_universe
4 Zadanie UCSB ScienceLine , O ilu gwiazdach jest kosmos? , od 2013 roku. https://scienceline.ucsb.edu/getkey.php?key=3775
5 AArtykuł Sky and Telescope , Ile gwiazd jest we wszechświecie? , od 2014 roku. https://www.skyandtelescope.com/astronomy-resources/how-many-stars-are-there/
6 Artykuł Space.com , Ile gwiazd jest we wszechświecie? , od 2017 r. https://www.space.com/26078-how-many-stars-are-there.html
7 Artykuł w Wikipedii , Obserwowalny wszechświat . https://en.wikipedia.org/wiki/Observable_universe
8 Artykuł w Wikipedii , Kula , sekcja Tom zamknięty . https://en.wikipedia.org/wiki/Sphere#Enclosed_volume
9 WolframAlpha obliczeniowej powierzchni kuli o średnicy 8,8 * 10 ^ 26 m . https://www.wolframalpha.com/input/?i=surface+area+of+a+sphere%2C+diameter+8.8*10%5E26+m
10 Artykuł w Wikipedii , Kula , sekcja Powierzchnia . https://en.wikipedia.org/wiki/Sphere#Surface_area
11 WolframAlpha obliczeniowy objętości kuli o średnicy 8,8 * 10 ^ 26 m . https://www.wolframalpha.com/input/?i=volume+of+a+sphere%2C+diameter+8.8*10%5E26+m
12 nineplanets.org artykuł, Słońce .https://nineplanets.org/sol.html
13 A Obliczenia WolframAlpha , (10 ^ 24 gwiazdki) / (3,568⋅10 ^ 80 m ^ 3) . https://www.wolframalpha.com/input/?i=%2810%5E24+stars%29+%2F+%283.568%E2%8B%8510%5E80+m%5E3%29
14 Artykuł w Wikipedii , Kąt bryłowy . https://en.wikipedia.org/wiki/Solid_angle
15 Odpowiedź Harisha Chandry Rajpoot na pytanie geometry.se , Obliczanie kąta bryłowego dla kuli w przestrzeni . https://math.stackexchange.com/a/1264753/265963
16 Artykuł w Wikipedii , Steradian .https://en.wikipedia.org/wiki/Steradian
17 A Obliczenia WolframAlpha , 2 * pi * (1-sqrt (d ^ 2- (1,4 * 10 ^ 9 m / 2) ^ 2) / d) . https://www.wolframalpha.com/input/?i=2*pi*%281-sqrt%28d%5E2-%281.4*10%5E9+m%2F2%29%5E2%29%2Fd%29
18 Strona internetowa dla matematyki. https://www.ttmath.org/
19 A Obliczenie WolframAlpha , 2 * pi * (1 - sqrt (d ^ 2 - r ^ 2) / d), gdzie d = 150 miliardów, r = 0,7 miliarda . https://www.wolframalpha.com/input/?i=2*pi*%281+-+sqrt%28d%5E2+-+r%5E2%29%2Fd%29%2C+where+d+%3D+150 + miliard% 2C + r% 3D 0,7 + miliard
20 A Obliczenia WolframAlpha , pi * (5 * 10 ^ 20 m) ^ 2 * (2 * 10 ^ 16 m) .https://www.wolframalpha.com/input/?i=pi+*+%285*10%5E20+m%29%5E2+*+%282*10%5E16+m%29
21 Artykuł w Wikipedii , Droga Mleczna . https://en.wikipedia.org/wiki/Milky_Way
22 Artykuł Space.com z 2018 roku, zajęłoby 200 000 lat z prędkością światła, aby przekroczyć Drogę Mleczną . https://www.space.com/41047-milky-way-galaxy-size-bigger-than-thought.html
23 A Obliczenia WolframAlpha , (200 * 10 ^ 9 gwiazdek) / (1,571 * 10 ^ 58 m ^ 3 ) . https://www.wolframalpha.com/input/?i=(200*10^9+stars)+%2F+(1.571*10^58+m^3)
24 A Obliczenia WolframAlpha ,rozwiąż dla r: (4/3) * pi * r ^ 3 = 1,571 * 10 ^ 58 m ^ 3 . https://www.wolframalpha.com/input/?i=solve+for+r%3A++%284%2F3%29*pi*r%5E3+%3D+1.571*10%5E58+m%5E3
25 Mój program C ++ kod na PasteBin . https://pastebin.com/XZTzeRpG
26 Fizyka Forum zakładać, orientacji Ziemi, Słońca i Układu Słonecznego w Drodze Mlecznej . W szczególności, rysunek 1 , pokazujący kąty 60,2 ° dla Słońca i 23,4 ° mniejsze niż dla Ziemi. https://www.physicsforums.com/threads/orientation-of-the-earth-sun-and-solar-system-in-the-milky-way.888643/