<< 1. Введение | Оглавление | 3. История звездообразования карликовых ... >>
2. Звездные индикаторы
Разделы
- 2.1 Анализ диаграмм "Цвет-Звездная величина"
- 2.2 Точки поворота ветви Главной Последовательности
- 2.3 Звезды на стадии горения гелия в ядрах - Голубая Петля
- 2.4 Ветвь Красных Гигантов
- 2.5 Красное Сгущение / Горизонтальная Ветвь
- 2.6 Продленная Асимптотическая Ветвь Гигантов
- 2.7 Расстояние, Поглощение и Металличность
- 2.8 Наблюдаемые CMDs
- 2.9 Функции светимости звезд
- 2.10 Изохроны и функции распределения металличности
- 2.11 Модельные диаграммы "Цвет-Зв.Величина"
- 2.12 Сравнение наблюдаемых CMDs с модельными
- 2.13 Возможность проверки различных методов
- 2.14 Некоторые оставшиеся проблемы
Некоторые типы звезд могут служить хорошими индикаторами расстояний, возраста и/или металличности.
Так как некоторые типы звезд видны только в определенные моменты
своей эволюции, то их отсутствие или присутствие может быть использовано
для грубой оценки возрастов, как показано на Рис.1. Например, OB-звезды,
или звезды типа Вольфа-Райе2
(WR) являются несомненными признаками недавнего звездообразования
(например, Armandroff & Massey 1991, Maeder & Meynet 1994).
Присутствие же OB звезд при отсутствии WR звезд позволяет выделить районы
недавнего звездообразования - либо до образования звезд WR, либо старее
4 млн.лет, когда звезды WR уже не видны. Общее число
и доля Be-звезд зависит, в основном, от спектрального типа и т.о. от возраста
(Mermilliod 1982, Grebel 1997).
|
|
Рис. 1. Звезды и оценки возрастов. Условная схема, показывающая как в первом приближении можно оценить возраст звездного населения в галактике с помощью некоторых типов звезд. Интервалы возрастов взяты из работ: Maeder & Meynet (1994), Frogel et al. (1990), Caputo et al. (1995) и Ibata et al. (1997). Возраста звезд всегда зависят от их масс, светимостей и металличности. В конечном счете, присутствие определенного типа звезд - достаточный критерий для населения в соответствующем диапазоне возрастов, в то же время отсутствие этого типа звезд (например, звезд типа RR Лиры) не обязательно означает отсутствие в галактике звезд в данном диапазоне возрастов. Сокращения: AGB - асимптотическая ветвь гигантов, HB - горизонтальная ветвь, Myr - 106 лет, Gyr - 109 лет. |
Относительное число звезд определенного типа может быть переведено в долю
всего поколения звезд с таким же возрастом. Так как средние по данной галактике
светимости углеродных
звезд (C-звезды) определяется их возрастом, средняя доля вклада C-звезд в
видимую яркость галактики может использоваться для получения доли,
которую составляют звезды "среднего возраста" (Aaronson & Mould 1985).
Аномальные цефеиды (звезды типа 
Цефея) могут также быть указателями
на население "среднего возраста", если они не являются результатом переноса
массы в двойных системах или слипания звезд, будучи одного возраста с
основным звездным населением. В то же время, "голубые бродяги", по-видимому, не являются
признаком "молодого" населения, а на самом деле - результат эволюций
двойных звезд (например, Stryker 1993, Bailyn & Pinsonneault 1995).
Ниже будут описаны морфологические критерии возраста звезд, основывающиеся на специальных типах звезд, учитывающих разницы между звездами разных возрастов.
2.1 Анализ диаграмм "Цвет-Звездная величина"
Звездообразование в галактиках наиболее полно можно исследовать из изучения разрешаемого звездного населения3. Когда галактика находится достаточно близко от нас, то становится возможным изучать наблюдаемое распределение ее звезд по цветам и показателям цвета.
Эта информация наиболее полно представляется в виде диаграмм "цвет-зв.величина" (CMD). Так как температуру звезды можно оценить из ее показателя цвета, а светимость - из звездной величины, то CMD -- есть ни что иное как известная диаграмма Герцшпрунга-Рессела, позволяющая классифицировать звезды.
Современная теория эволюции звезд позволяет предсказывать те детали, которые мы должны наблюдать в разных частях диаграммы "Цвет-Зв.величина" (CMD) для звезд разной металличности и возраста (см. Рис. 2).
|
|
Рисунок 2: Изохроны для металличности (Z=0.001) и в диапазоне возрастов, указанных в миллиардах лет. [4], возле MSTOs. Эти изохроны были посчитаны для звезд шаровых скоплений, которые, как известно, одного возраста. Для интерпретации диаграмм со смешанным по возрасту звездным населением необходимо учитывать всю совокупность изохрон с разными возрастами, например, используя методы Монте-Карло [37]. |
Существуют чуткие индикаторы изменения скоростей звездообразования (sfr) для разных эпох эволюции галактик, комбинируя которые можно получить очень точную картину полной истории звездообразования (SFH) галактик.
2.2 Точки поворота ветви Главной Последовательности
Main Sequence Turnoffs (MSTOs)
Используя большие современные телескопы при достаточно глубоких экспозициях4 есть шанс профотометрировать звезды находящиеся в точках поворота ветви главной последовательности (Main Sequence Turnoffs - MSTOs) для достаточно близких галактик (МГ), светимость которых позволяет получить конкретную информацию о SFH. Как известно, звезды попадают на главную последовательность после того, как в их ядрах начинается сгорание водорода, и остаются на ней до сгорания около 12% его запасов. Более массивные звезды раньше покидают ГП по сравнению с медленно эволюционирующими маломассивными, что и образует MSTO. Вдоль Главной Последовательности (ГП) диаграммы Герцшпрунга-Рессела (Г-Р) звезды различных возрастов накладываются полностью друг на друга, делая сложной интерпретацию функции светимостей звезд ГП, особенно для старого звездного населения. Однако точки поворота не перекрываются и поэтому обеспечивают наиболее прямую и точную информацию о SFH в галактике. Используя MSTOs можно легко разделить по времени периоды вспышки звездообразования и спокойные стадии образования звезд, например [22]. Точность временного разрешения различна, она ухудшается для более древних времен. Наша способность распутать клубок изменений sfr от времени зависит от нескольких факторов: интенсивности прошлых колебаний; от того, как долго они происходили; и от того, какие фильтры использовались для наблюдений. Для звезд с возрастом менее 1.5 млрд.лет возможно построить SFH с разрешением 10-100 млн.лет. Как видно из Рисунка 2, возраста становятся неразличимыми для более старого звездного населения. Для возрастов больше 8 млрд.лет разрешение метода по времени около 1 млрд.лет. Это означает, что очень короткие бурные вспышки звездообразования в такие эпохи "размазываются" по периоду в млрд.лет и кажутся менее интенсивными и имеющими большую продолжительность. Однако, как будет описано в следующих пунктах, существуют другие индикаторы на CMD, которые позволяют сузить интервал возможных вспышек в истории звездообразования.
2.3 Звезды на стадии горения гелия в ядрах - Голубая Петля
Blue Loop (BL)
Звезды с определенной металличностью и массой (
4 масс Солнца)
на определенной стадии
своей эволюции (после того, как в их ядрах начинается термоядерное горение
гелия) "двигаются" на диаграмме Герцшпрунга-Рессела (Г-Р) по
эволюционным трекам, которые названы "путешествия по голубым петлям".
Звезды в BL-стадии располагаются вблизи главной последовательности (ГП),
но лежат на несколько величин выше ее точки поворота (для звезд одного
возраста). Использование этих более ярких звезд позволяет также достаточно
точно определять возраст и металличность молодого звездного населения
(в диапазоне 
1 млрд.лет) в близких галактиках с низкой
металличностью
[10,11].
Форма "петель" и масса звезд при которых они переходят в BL-стадию
сильно зависят от металличности и возраста, а их светимость постоянна
для данного возраста. Следующие поколения BL-звезд не перекрываются
друг с другом, в отличие от звезд ГП. Чем ниже металличность в галактике,
тем старее должны быть самые пожилые BL-звезды и тем самым более древнюю
историю звездообразования мы можем узнать. Часто BL-звезды ярче чем
звезды ГП такой же массы, на одну звездную величину и более.
2.4 Ветвь Красных Гигантов
Red Giant Branch (RGB)
Красные гиганты - это звезды высокой светимости
(
),
находящиеся на стадии эволюции после ГП,
в которых термоядерное горение водорода происходит в оболочке,
окружающей их гелиевые ядра. Несмотря на то что
RGB является хорошим индикатором расстояний (метод точки верхнего обрыва
RGB [Tip of the RGB - TRGB]), такое же нельзя сказать для оценки
возрастов. Дело в том, что для данной металличности голубая и
красная границы RGB определяется границами по возрасту самых молодых и
самых старых звезд, населяющих RGB (возраст 1 млрд.лет).
С возрастом RGs смещаются на диаграмме Г-Р в красную сторону,
при постоянной металличности. Голубой край ветви красных гигантов
определяется возрастом старейших звезд. Однако увеличение металличности
звездного населения приводит к такому же эффекту и также смещает RGB
в красную сторону. Это известная проблема
возраст-металличность(Рис.3a).
В результате, если в галактике происходило
изменение металличности со временем, то мы не можем однозначно разделить
эффекты, являющиеся следствием возраста и металличности на основе лишь
оптического цвета RGB.
|
|
Рисунок 3: Левая панель:
Эффекты металличности и возраста для ветви красных гигантов (RGB).
Все изохроны из Berbusch & VandenBerg (1992). Левое подмножество изохрон
сдвинуто на -0.5 зв.величин по V-I, для того, чтобы показать оба набора
на одной диаграмме. Так как место и наклон RGB чувствительны к
металличности,
ветвь является плохим индикатором возраста для "старого" населения.
Отметим ухудшение степени разрешения по металличности с уменьшением
металличности.
|
2.5 Красное Сгущение / Горизонтальная Ветвь
Red Clump/Horizontal Branch (RC/HB)
Звезды из Красного Сгущения [Red Clump (RC)] на CM-диаграмме и маломассивные родственные им звезды из Горизонтальной Ветви [Horizontal Branch (HB)] - звезды, в ядрах которых горит гелий (и водород в окружающей оболочке). Светимость этих звезд изменяется со временем, а также при изменении металличности и с потерей массы звездами [6]. Диапазон светимостей в наблюдаемом Красном Сгущении может использоваться для оценки возрастов звезд, образующих его [6], как показано на Рисунке 4. Этот индикатор возрастов звезд не зависит от абсолютных зв.величин и следовательно от расстояния до галактики, и в действительности эти свойства могут быть использованы для оценок расстояний с хорошей точностью по звездам из RC [7].
|
|
Рисунок 4: На верхнем графике показаны результаты Caputo, Castellani & Degl'Innocenti [6] для изменений размера Красного Сгущения (Red Clump - RC) в абсолютных зв.величинах MV с возрастом, для металличности Z=0.0004. Нарисовано изменение зв.величины верхнего и нижнего краев светимости RC в зависимости от возраста (в млрд.лет). Ясно видно, что этот размер сильно зависит от возраста звезд, населяющих его. Также нарисован график MV начальной горизонтальной ветви (Zero Age Horizontal Branch - ZAHB) в зависимости от возраста. На нижнем графике нарисованы результаты численного моделирования методом Монте-Карло [37], используя модели звездной эволюции при металличности Z=0.0004 [16], для относительного числа звезд RC/RGB в одной и той же части CM-диаграммы в зависимости от времени. |
Классические звезды из RC и RGB возникают в населении почти в одно и
время (
0.9-1.5 млрд.лет, в зависимости от различных начальных
условий в модели), где звезды из RGB являются предшественниками звезд
из RC. Время жизни звезды на RGB, tRGB, сильно уменьшается с ростом
массы звезды Mstar, но время жизни в RC, tRC, почти постоянно.
Следовательно отношение tRC / tRGB - есть убывающая функция
возраста доминирующего звездного населения в галактике, а отношение
числа звезд в RC и HB к числу звезд из RGB чувствительно к SFH галактики
[39,20]. Т.о. чем больше отношение N(RC)/N(RGB), тем моложе
доминирующее звездное население в галактике, как показано на
Рисунке 4.
Присутствие многочисленного населения Горизонтальной Ветви, с другой стороны, (большое N(HB)/N(RGB) или даже N(HB)/N(MS)) - есть свидетельство преобладания в галактике значительно более старого звездного населения (>10 млрд.лет). HB - самый яркий индикатор звезд самых малых масс в галактике (следовательно самых старых).
2.6 Продленная Асимптотическая Ветвь Гигантов
Extended Asymptotic Giant Branch (EAGB)
Температура и цвет EAGB-звезд в галактике определяется возрастом и металличностью звезд, которые ее образуют (см. Рисунок 2). Однако остается некоторое число неопределенностей при сравнении моделей и наблюдательных данных [18,26]. Очень важным фактом является то, что уже проделана большая работа по надежной калибровки этих очень ярких индикаторов событий в истории звездообразования. На Рисунке 5 теоретические изохроны звезд из EAGB [4] наложены на СМ-диаграмму галактики, находящейся на стадии после вспышки звездообразования в виде Голубого Компактного Карлика (BCD), на которой видно, что при большом числе EAGB-звезд последние являются ярким индикатором прошлой высокой скорости звездообразования (sfr), а их распределение по светимостям зависит от металличности и возраста прошлой вспышки звездообразования.
|
|
Рисунок 5: Изохроны EAGB-звезд [4] для металличностей, Z=0.001 и Z=0.004, наложены на наблюдаемую CM-диаграмму в галактике VII Zw403 [26]. Для каждой металличности приведены изохроны для звезд с возрастом 1.3, 2, 3, и 5 млрд.лет, так что самая молодая изохрона является ярчайшей. Видно возможное расхождение в будущих оценках металличности и возраста самых старых звезд, когда модели будут лучше прокалиброваны по галактикам с известной SFH, например, по галактике NGC 6822, в которой наблюдается большое число EAGB-звезд, а также видны старые MSTO. |
2.7 Расстояние, Поглощение и Металличность
Аккуратная интерпретация индикаторов, описанных выше зависит критически от наличия надежных оценок расстояния, поглощения и металличности галактики. В идеале, нам также хотелось бы знать, как меняется поглощение в разных частях галактики и как проходила эволюция металличности звездного населения со временем.
Основные параметры: расстояние, поглощение и металличность могут быть определены независимо по CM-диаграмме (CMD). Эти параметры, в сочетании с оценками ошибок наблюдений и неполноты данных вносят наиболее значительный вклад в получаемую CMD и следовательно в конечную модель SFH [37]. Получение надежных оценок этих величин - достаточно непростая задача, но она может быть разрешена путем высококачественных наблюдений и продуманных методов обработки и анализа.
2.7.1 Расстояние
- наиболее важный параметр для правильной оценки CMD, частично из-за того, что ошибки могут быть очень велики - несколько величин (например, если ошибочно принять молодые красные сверхгиганты (RSG) за RGB, при недостаточно глубокой для правильного отождествления CM-диаграмме, а именно: из-за отсутствия звезд из RC или HB). Если расстояние до галактики оценено неправильно, то это приведет к неправильной оценке масс индивидуальных звезд, и вследствии этого оценки возрастов звезд всех типов будут определены неверно. Теоретически, расстояния наиболее точно должны определяться с помощью первичных индикаторов расстояний (например, переменных типа RR Лиры или Цефеид), но существуют также и другие методы со сходной точностью: верхний конец ветви красных гигантов [Tip of the RGB (TRGB)] [25]; RC [7]; BLs [39]. Для тусклых карликовых галактик, которые составляют большую часть Местной Группы для оценки расстояний необходимо иметь достаточно глубокую CMD, захватывающую в своем слабом конце RC/HB.2.7.2 Поглощение,
как внутри исследуемой галактики, так и в диске нашей Галактики может повлиять на правильность анализа CMD. Если поглощение неправильно определено, то результат от этого будет такой же, как от ошибки в оценке расстояния и следовательно повлияет на надежность моделей SFH. Для учета поглощения света при прохождении сквозь нашу Галактику часто используют радио-карты распределения нейтрального водорода, а также инфракрасные карты (например, IRAS) на длине волны 100 микрон [17].2.7.3 Металличность.
Когда в галактиках рождаются звезды, то осколки от этого процесса (например, от вспышек сверхновых или вследствии звездного ветра) обогащают межзвездный газ тяжелыми элементами, т.е. изменяет металличность материала, который впоследствии идет на создание новых звезд5 [14]. Несмотря на отсутствие конкретных наблюдательных подтверждений этого процесса, отношения содержаний различных элементов позволяют нам сделать некоторые зависящие от модели предположения [29]. Например, недавно было показано, что в диске нашей Галактики есть общий тренд в эволюции металличности со временем, хотя и с большими вариациями [13]. Мы не понимаем всех деталей взаимодействия звезд с окружающей межзвездной средой, например, как текущий процесс звездообразования "подпитывает" металлами будущие поколения. Последние результаты изучения содержания цинка по линиям поглощения на расстояниях (z = 0.7 - 3.4) свидетельствуют о медленной эволюции металличности в галактиках на этом широком диапазоне красных смещений. Однако в любую из эпох (на одном z) наблюдается большое рассеяние в оценке величины металличности у разных галактик[30]. Изучения подобных линий поглощения позволяют получить наиболее надежные оценки металличности межзвездного газа в галактиках. Если подходящие источники фонового континуума будут найдены позади близких галактик, то это значительно улучшит наше понимание того, как межзвездный газ в различных галактиках эволюционирует и как на него воздействует соседство с текущим звездообразованием.
Трудно проводить расчеты эволюции металличности в близких галактиках по CMD. Невозможно определить одну единственную модель, основываясь только на RGB вследствии противоречия возраст-металличность. Однако, если пренебрегать эволюцией металличности в модели CMD, то наилучшие модели для галактик становятся моложе, чем если учитывать эволюцию металличности [39].
Понимание деталей эволюции металличности в галактиках - одна из наиболее важных областей дальнейшего изучения эволюций галактик.
2.8 Наблюдаемые CMDs
|
|
Рисунок 6: CM-диаграмма (V,B-V) для dSph галактики Sextans (левая панель взята из работы Mateo et al. 1995) и для dIr галактики NGC3109 (правая панель взята из работы Greggio et al. 1993). Показаны основные морфологические детали. Сокращения: MS = главная последовательность, BSG = голубые сверхгиганты, AGB = асимптотическая ветвь гигантов, RGB = ветвь красных гигантов, HB = горизонтальная ветвь, фон = фоновые звезды. Для CMD NGC3109 интерпретация сделана на основе работ Lee (1993) и Davidge (1993). |
Наблюдаемые CMDs близких галактик показывают смесь населений разных типов, что может значительно затруднить их анализ. Наиболее сложные CMDs у иррегулярных галактик. Для карликовых сфероидальных галактик наиболее выдающаяся деталь на CMD - широкая RGB, в паре с более или менее выраженной HB, или сгущением HB, или их обеих (Рис. 6). При достаточно глубоких экспозициях можно увидеть и MSTOs. Для CMDs иррегулярных галактик характерна концентрация звезд вблизи верхней части RGB/AGB, выше часто наблюдаются красные сверхгиганты (RSG), а также верхняя часть главной последовательности + голубые сверхгиганты (BSG). Вследствии зачастую сильной скученности звезд и увеличения ошибок для фотометрии для более слабых звезд происходит уширение структур на CMD: районов RGB/AGB/RSG, что делает интерпретацию CMDs еще более трудной (Рис. 6). Для анализа CMDs со смесью основного "среднего возраста"/"старого" населения с большой долей "молодых" звезд часто выбирают совершенно различные подходы.
2.9 Функции светимости звезд
Mighell (1990, Mighell& Butcher 1992) разработали метод изучения SFH населений "среднего возраста" путем вписывания в наблюдаемые функции светимостей звезд (SLF) модельные SLF, составленных для некоторого диапазона параметров и модулей расстояний до галактики. Основными параметрами для модельных SLF являются: возраст, металличность, содержание гелия и начальная функция масс (IMF). Правильная интерпретация запутанной истории звездообразования "старого" населения галактик требует большого числа звезд в SLF, для уменьшения противоречия "возраст-металличность". Если это условие выполнено, то этот метод позволяет последовательно получить возраст, металличность, и модуль расстояния до галактики.
Ярко выраженный перегиб кривой SLF для звезд из точки обрыва RGB (TRGB) может быть указателем на наличие "старого" населения (например, Minniti & Zijlstra 1996). Сравнение наблюдаемых SLFs шаровых скоплений с SLFs у dSph (Mateo et al. 1991,1995) показало избыток части звезд, находящихся вблизи RGB, возможно, вследствии попадания AGB-звезд от населения "среднего возраста". Наблюдаемый сверхизбыток звезд в SLF главных последовательностей у dSph галактик может быть объяснено скорее наличием в dSph "молодого" населения, чем необычайно большим содержанием "голубых бродяг". SLF "молодого" населения состоит из голубого "крыла" неразделимых массивных MS-звезд и BSG (Marconi et al. 1995), возможно, перемешанных со "среднего возраста" звездами "голубой петли" (BL).
2.10 Изохроны и функции распределения металличности
Grebel и др. (Grebel et al. 1994, Grebel 1995) разработали метод последовательного определения возраста, хим.состава, покраснения и модуля расстояния посредсвом последовательного вписывания изохрон в многоцветную фотометрию по крайней мере в четырех фильтрах. Четыре фильтра формально достаточно для исключения противоречия возраст-металличность. Вписывания производятся для диапазона величин по всем четырем параметрам и статистически оцениваются.
Использование особой системы фильтров, чувствительных к изменению металличности (например, широкополосной вашингтонской системы), позволяет определять содержание металлов у большого числа звезд (Harris & Canterna 1977, Geisler et al. 1991). Для определения металличности смешанного "старого" и "среднего возраста" населений в dSphs, можно использовать цвета красных гигантов M-T2 и C-M, независимо от возраста и изохрон. Результирующие функции рапределения металличности могут выделить поднаселения с разной металличностью, исключая проблемы с противоречием возраст-металличность.
|
|
Рисунок 7: Модельные CMDs для постоянной SFR, начиная с 15 млрд.лет назад до настоящего времени и линейно увеличивающейся Z(t) с Z=0.0001 до Z=0.0004. На каждой панели нанесены звезды с указанным (в млрд.лет) интервалом возрастов. Смоделированы наблюдательные эффекты, как описано в пункте 2.11. |
2.11 Модельные диаграммы "Цвет-Зв.Величина"
Процесс создания модельной CMD можно разбить на следующие два основных этапа:
- вычисление синтетической CMD и
- моделирование наблюдательных эффектов.
2.12 Сравнение наблюдаемых CMDs с модельными
Статистические оценки могут помочь уменьшить вырождение синтетических CMDs. Tolstoy (1995; Tolstoy & Saha 1996) разработали количественный метод сравнения наблюдаемых данных с синтетическими CMDs для смешанных населений, исследующих вероятность того, что модельная и наблюдаемая CMDs соответсвуют одному и тому же распределению.
2.13 Возможность проверки различных методов
Тривиальная проверка правильности какого-либо выбранного метода - сравнить его результаты с независимыми индикаторами возраста, такими как специальные типы звезд, либо с независимыми индикаторами металличности (содержание элементов), найденными из спектральных наблюдений или ветви сравнения шаровых скоплений. На Рис.3 (правая панель) показано сравнение теоретических моделей изохрон с наблюдаемыми RGB у шаровых скоплений с надежно установленной металличностью. Для всех методов, которые используют теоретические модели (эволюционные треки звезд и/или изохроны, модели Kurucz'а и т.п.) важно знать насколько эти модели соответсвуют реальным физическим свойствам звезд. Многие модели довольно часто перерабатываются и улучшаются.
Карликовая галактика в Киле (Carina dSph) - хорошо подходит для проверки чувствительности различных методов оценки возрастов звезд. Наблюдения показывают, что в Carina присутсвуют несколько MSTOs, в то время как разброс по металличности незначителен (например, Smecker-Hane et al 1996). Т.о. в ней хорошо выявлены возраста различных поколений звезд. Существующие данные по этой галактики могут быть ограничены по абсолютным звездным величинам с целью проверки различных методов анализа.
2.14 Некоторые оставшиеся проблемы
В целом, для надежной интерпетации разрешаемых звездных населений требуется как можно более качественные наблюдательные данные, полученные для как можно более слабых звездных величин. В этом деле требуется комбинация как высоких разрешений (наподобие получаемых на HST), так и наблюдения больших по угловым размерам площадок. Касательно теоретической базы, лежащей в основе анализа любым методом, требуется наличие высококачественных моделей эволюции звезд разнообразных типов.
<< 1. Введение | Оглавление | 3. История звездообразования карликовых ... >>
|
Публикации с ключевыми словами:
звездообразование - Местная группа галактик - ближайшие галактики
Публикации со словами: звездообразование - Местная группа галактик - ближайшие галактики | |
См. также:
Все публикации на ту же тему >> | |
