Резюме
Вероятность того, что вы стоите под звездой за пределами Млечного пути, составляет 1 на 500 миллиардов, вероятность того, что вы находитесь под звездой Млечного пути, составляет 1 на 3,3 миллиарда, а вероятность того, что вы находитесь под солнцем справа, равна 1 на 184 тысячи. в настоящее время.
Большой, толстый, вонючий, Внимание! Я сделал все возможное, чтобы сохранить свою математику, но это все, что я только что придумал. Я не даю никаких гарантий, что это абсолютно точно, но цифры, кажется, проходят проверку на исправность, поэтому я думаю, что мы хороши.
Предостережение Первое : числа звезд, отличных от Солнца, основаны на данных с большой степенью неопределенности, таких как количество звезд во вселенной и средний размер звезды. Числа, приведенные выше, могут легко отклоняться в 10 раз в любом направлении, и они просто предназначены для того, чтобы дать приблизительное представление о том, насколько пустое пространство.
Предостережение второе : числа Солнца и Млечного Пути основаны на предположении, что вы стоите (или летаете) в случайной точке на Земле. Любой, кто находится вне тропиков, никогда не будет иметь Солнца над головой. Люди в северном полушарии с большей вероятностью будут иметь над головой звезды Млечного Пути, при этом лучшими шансами будут люди около 36,8 ° с.ш., потому что на этой широте прямая линия вверх проходит через центр галактики один раз в день. 26
Примечание . В этом ответе вы можете в основном игнорировать все и просто посмотреть на телесный угол Солнца, чтобы получить тот же результат. Все остальные звезды действительно далеко и очень раскинулись. Разница в выраженном телесном угле на пять тысячных процента больше, когда мы добавляем остальную часть вселенной к Солнцу.
Фон
Давайте попробуем получить несколько реалистичный, жесткий номер. Для этого нам понадобятся некоторые предположения.
Как указано в ответе Майкла Уолсби 1 , если вселенная бесконечна (и однородна 2 ), существует только бесконечно малая вероятность того, что над головой не будет звезды, что в обычной математике рассматривается как абсолютно нулевой шанс. Итак, давайте предположим, что вселенная конечна.
Предположения
- В частности, давайте предположим, что вселенная состоит только из наблюдаемой вселенной. (Посмотрите расширение вселенной 3 для получения дополнительной информации.)
- Далее, давайте предположим, что содержимое наблюдаемой вселенной измеряется в их текущем (предполагаемом) положении, а не в том положении, в котором они находятся. (Если мы увидим свет от звезды через 400 миллионов лет после начала существования Вселенной, мы измерим его на расстоянии около 13,5 миллиардов световых лет, но мы рассчитываем, что он, вероятно, ближе к 45 миллиардам световых лет от расширения.)
- Мы примем число звезд в наблюдаемой вселенной равным . Оценка 4 в 2013 году составляла , оценка 5 в 2014 году - , а оценка 6 в 2017 году - , и каждая статья ожидала, что оценка будет увеличиваться по мере того, как со временем мы получим лучшие телескопы. Поэтому мы возьмем самое высокое значение и будем его использовать.1024 10 21 10 23 10 24102110231024
- Мы возьмем размер наблюдаемой Вселенной 7 будет , что дает площадь поверхности 8 из 9 и том 10 из 11 .8.8⋅1026m (diameter)2,433 × 10 54 м 2 3,568 × 10 80 м 32.433⋅1054m2 3.568⋅1080m3
- Мы возьмем средний размер звезды равным размеру Солнца, 12 . (Я не могу найти никаких источников для среднего размера звезды, просто Солнце - средняя звезда.)1.4⋅109m (diameter)
модель
Отсюда мы собираемся немного обмануть. Реально, мы должны моделировать каждую галактику отдельно. Но мы просто собираемся притвориться, что вся вселенная совершенно однородна (это достаточно верно, когда мы удаляемся от Земли в грандиозной схеме космоса). Далее, мы собираемся начать отсчет достаточно далеко, чтобы полностью игнорировать Млечный Путь и Солнце, а затем добавим их позже с другими вычислениями.
Учитывая вышеизложенные предположения, мы можем легко вычислить звездную плотность наблюдаемой вселенной как: 13 .δ=1024stars3.568⋅1080m3=2.803⋅10−57starsm3
Далее нам нужно вычислить телесный угол 14 , обозначенный звездой. Телесный угол сферы задается как 15 , где - телесный угол в стерадианах 16 (ср), - расстояние до сферы, а - радиус сферы. Используя в качестве диаметра, это преобразуется в . Учитывая предполагаемый выше средний диаметр ( ), это дает средний телесный уголΩ=2π(1−d2−r2√d) sr ΩdrDΩ=2π⎛⎝⎜1−d2−(D2)2√d⎞⎠⎟ sr1.4⋅109mΩ=2π(1−d2−4.9⋅1017m2√d) sr 17 .
На этом этапе мы могли бы установить правильный интеграл, но мое исчисление довольно ржавое и не очень резкое для начала. Поэтому я собираюсь приблизить ответ, используя серию концентрических оболочек, каждая из которых имеет толщину (около миллиона световых лет). Мы отложим нашу первую оболочку на , а затем оттуда доберемся.1022m1022m
Мы вычислим общий телесный угол каждой оболочки, затем сложим все оболочки вместе, чтобы получить телесный угол, охватываемый всей наблюдаемой вселенной.
Последняя проблема, которую нужно исправить, - это проблема перекрытия. Некоторые звезды в дальних оболочках перекрывают звезды в близлежащих оболочках, в результате чего мы переоцениваем общее покрытие. Таким образом, мы рассчитаем вероятность того, что любая звезда будет перекрываться, и изменим результат оттуда.
Мы будем игнорировать любое перекрытие внутри данной оболочки, моделируя, как будто каждая звезда в оболочке находится на фиксированном расстоянии, равномерно распределенном по всей оболочке.
Вероятность перекрытия
Чтобы данная звезда перекрывала более близкие звезды, она должна находиться в положении, уже покрытом более близкими звездами. Для наших целей мы будем рассматривать перекрытия как двоичные: либо звезда полностью перекрывается, либо не перекрывается совсем.
Вероятность будет определяться количеством телесного угла, уже существующего в предыдущих оболочках, деленным на общий телесный угол в небе ( ).4π sr
Давайте назовем вероятность того, что данная звезда перекрывается с , телесным углом, представленным этой звездой , и количеством звезд . Тогда количество неперекрывающихся телесных углов, представленных данной оболочкой, равно . Поскольку мы сказали, что звезды в оболочке не перекрывают друг друга, одинаково для всех в данной оболочке, что позволяет упростить приведенное выше уравнение до , гдеiPiΩinkΩkT=(1−P1)Ω1+(1−P2)Ω2+…+(1−Pn)Ωn srstarPiiΩkT=(1−Pk)(Ω1+Ω2+…+Ωn) srstarPkвероятность перекрытия для оболочки . Поскольку мы рассматриваем все звезды как имеющие одинаковый средний размер, это еще больше упрощается до , где - телесный угол звезды в оболочке .kΩkT=(1−Pk)Ωkn srstarΩkk
Расчет телесного угла
Количество звезд в оболочке определяется как объем оболочки, умноженный на звездную плотность этой оболочки. Для далеких оболочек мы можем рассматривать объем оболочки как площадь ее поверхности, умноженную на ее толщину. , где - расстояние до оболочки, а - ее толщина. Используя в качестве звездной плотности, количество звезд просто .Vshell=4πd2tdtδn=δVshell=δ4πd2t
Отсюда мы можем использовать расчет для телесного угла оболочки (из вероятности перекрытия , выше), чтобы получить .ΩkT=(1−Pk)Ωkδ4πd2t srstar
Обратите внимание, что задается частичной суммой телесного угла для всех предыдущих оболочек, деленной на общий телесный угол. И задается как (из модели выше).PkΩkΩk=2π(1−d2k−4.9⋅1017m2√dk) srstar
Это дает нам . Учитывая, что каждая оболочка находится на расстоянии , мы можем заменить на . Аналогично, можно заменить на . И мы уже вычислили (из модели выше).ΩkT=(1−Ω(k−1)T4π)2π(1−d2k−4.9⋅1017m2√dk)δ4πd2t sr1022mdkk1022mt1022mδ=2.803⋅10−57starsm3
Это дает нам
Ω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
Отсюда мы можем просто вставить числа в программу расчета.
ΩT=∑kmaxk=1ΩkT
Где - это радиус наблюдаемой вселенной, деленный на толщину данной оболочки. Таким образом,kmaxkmax=4.4⋅1026m1022m=4.4⋅104=44000
ΩT=∑44000k=1ΩkT
Результаты
Из-за большого числа вовлеченных в него программ трудно просто запустить это. Я прибег к написанию собственной программы на C ++ с использованием библиотеки 18 ttmath для больших чисел. В результате получилось или всего неба. И наоборот, вероятность того, что вы сейчас находитесь под звездой, составляет приблизительно 1 из 500 миллиардов.2.386⋅10−11 sr1.898⋅10−12
Обратите внимание, что мы игнорировали Млечный Путь и Солнце для этого.
Программу на C ++ можно найти на PasteBin 25 . Вы должны заставить ttmath работать должным образом. Я добавил несколько инструкций в начало кода C ++, чтобы вы могли начать, если хотите, чтобы он работал. Это не элегантно или что-то, просто достаточно, чтобы функционировать.
Солнце
Вольфрам Альфа любезно сообщил мне, что Солнце имеет телесный угол около , или примерно в 2,8 миллиона раз больше, чем у всех звезд во вселенной вместе взятых. Приведенная выше формула телесного угла дает тот же ответ 18, если мы предоставим расстояние в 150 Гб от Солнца и радиус в 0,7 Гб.6.8⋅10−5 sr
Млечный путь
Мы могли бы получить приближение для Млечного Пути, взяв его размер и плотность и выполнив те же вычисления, что и выше, за исключением меньшего масштаба. Однако галактика очень плоская, поэтому шансы сильно зависят от того, стоите ли вы в галактической плоскости или нет. Кроме того, мы отошли на одну сторону, так что звезд к галактическому центру гораздо больше, чем далеко.
Если мы приближаем галактику как цилиндр с радиусом (около 52000 световых лет) и высотой (примерно 2 световых года), мы получаем объем 20 .5⋅1020 m2⋅1016 m1.571⋅1058 m3
Текущие оценки радиуса галактики ближе к 100000 световых лет 21 22 , но я предполагаю, что подавляющее большинство звезд намного ближе к этому.
По оценкам, в Млечном Пути 21 от 100 до 400 миллиардов звезд . Давайте выберем 200 миллиардов для наших целей. Это определяет плотность Млечного Пути на уровне 22 , или примерно в 4,5 миллиарда раз плотнее, чем вселенная в целом.δ=200⋅109stars1.571⋅1058 m3=1.273⋅10−47starsm3
На этот раз мы возьмем снаряды толщиной (около 10 световых лет) и выйдем оттуда. Но нам нужно реорганизовать математику в сферическую форму, поэтому мы предположим, что галактика имеет такой же объем, но является сферой. Это дает радиус 24 или 155.4 оболочки. Мы округлим до 155 снарядов.1017 m1.554⋅1019 m
ΩT=∑155k=1ΩkT
Используя нашу формулу сверху ( Расчет телесного угла ), мы можем начать заменять числа.
Ω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
Включение этого в программу дает , что составляет от общего неба. Вероятность того, что вы окажетесь под звездой в Млечном Пути, составляет около 1 на 3,3 миллиарда.3.816⋅10−9 sr3.037⋅10−10
Твердые углы
Телесный угол:
- Солнце,6.8⋅10−5 sr
- Млечный Путь,3.816⋅10−9 sr
- Вселенная,2.386⋅10−11 sr
- Итого, (лишние цифры в основном бессмысленны, прибавляя около пяти тысячных процента к телесному углу Солнца) 6.800384⋅10−5 sr
- Млечный Путь плюс Вселенная, (примерно на 0,6% больше, чем просто Млечный Путь)3.840⋅10−9 sr
Рекомендации
1 Ответ Майкла Уолсби на этот вопрос , есть ли звезда над моей головой? , https://astronomy.stackexchange.com/a/33294/10678
2 Статья в Википедии , Космологический принцип . https://en.wikipedia.org/wiki/Cosmological_principle
3 Статья в Википедии , Расширение Вселенной . https://en.wikipedia.org/wiki/Expansion_of_the_universe
4 UCSB ScienceLine квест, о том , сколько звезд в космосе? , с 2013 года. https://scienceline.ucsb.edu/getkey.php?key=3775
5 AСтатья " Небо и телескоп" , сколько звезд во Вселенной? , с 2014 года. https://www.skyandtelescope.com/astronomy-resources/how-many-stars-are-there/
6 Статья на Space.com , сколько звезд во Вселенной? , с 2017 года. https://www.space.com/26078-how-many-stars-are-there.html
7 Статья в Википедии , « Наблюдаемая вселенная» . https://en.wikipedia.org/wiki/Observable_universe
8 Статья в Википедии , Сфера , раздел Вложенный том . https://en.wikipedia.org/wiki/Sphere#Enclosed_volume
9 WolframAlpha расчет, площадь поверхности сферы, диаметр 8,8 * 10 ^ 26 м . https://www.wolframalpha.com/input/?i=surface+area+of+a+sphere%2C+diameter+8.8*10%5E26+m
10 Статья из Википедии , Сфера , раздел Площадь поверхности . https://en.wikipedia.org/wiki/Sphere#Surface_area
11 WolframAlpha расчет, объем сферы, диаметр 8,8 * 10 ^ 26 м . https://www.wolframalpha.com/input/?i=volume+of+a+sphere%2C+diameter+8.8*10%5E26+m
12 Статья nineplanets.org , Солнце .https://nineplanets.org/sol.html
13 WolframAlpha Расчет, (10 ^ 24 звезды) / (3.568⋅10 ^ 80 м ^ 3) . https://www.wolframalpha.com/input/?i=%2810%5E24+stars%29+%2F+%283.568%E2%8B%8510%5E80+m%5E3%29
14 Статья в Википедии , Solid angle . https://en.wikipedia.org/wiki/Solid_angle
15 Ответ Хариша Чандры Раджпута на вопрос geometry.se , Расчет телесного угла для сферы в пространстве . https://math.stackexchange.com/a/1264753/265963
16 Статья в Википедии , Steradian .https://en.wikipedia.org/wiki/Steradian
17 WolframAlpha расчет, 2 * пи * (1-SQRT (д ^ 2- (1,4 * 10 ^ 9 м / 2) ^ 2) / г) . https://www.wolframalpha.com/input/?i=2*pi*%281-sqrt%28d%5E2-%281.4*10%5E9+m%2F2%29%5E2%29%2Fd%29
18 Веб-сайт для ttmath. https://www.ttmath.org/
19 WolframAlpha расчет, 2 * пи * (1 - SQRT (д ^ 2 - г ^ 2) / г), где d = 150 млрд, г = 0,7 млрд . https://www.wolframalpha.com/input/?i=2*pi*%281+-+sqrt%28d%5E2+-+r%5E2%29%2Fd%29%2C+where+d+%3D+150 +% 2C млрд + г +% 3D0.7 млрд
20 WolframAlpha расчет, пи * (5 * 10 ^ 20 м) ^ 2 * (2 * 10 ^ 16 м) .https://www.wolframalpha.com/input/?i=pi+*+%285*10%5E20+m%29%5E2+*+%282*10%5E16+m%29
21 Статья из Википедии , Млечный путь . https://en.wikipedia.org/wiki/Milky_Way
22 Статья Space.com, опубликованная в 2018 году, со скоростью света 200 000 лет, чтобы пересечь Млечный путь . https://www.space.com/41047-milky-way-galaxy-size-bigger-than-thought.html
23 WolframAlpha Расчет, (200 * 10 ^ 9 звезды) / (1,571 * 10 ^ 58 м ^ 3 ) . https://www.wolframalpha.com/input/?i=(200*10^9+stars)+%2F+(1.571*10^58+m^3)
24 WolframAlpha расчет,решить для r: (4/3) * pi * r ^ 3 = 1.571 * 10 ^ 58 м ^ 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 Моя программа на C ++ код на PasteBin . https://pastebin.com/XZTzeRpG
26 Physics Форум Начать новую тему ориентации Земли, Солнца и Солнечной системы в Млечном Пути . В частности, на рисунке 1 показаны углы 60,2 ° для Солнца и на 23,4 ° меньше, чем для Земли. https://www.physicsforums.com/threads/orientation-of-the-earth-sun-and-solar-system-in-the-milky-way.888643/