Документ взят из кэша поисковой машины. Адрес оригинального документа : http://ip.rsu.ru/~marsakov/paper4/article4.htm
Дата изменения: Fri Feb 2 12:03:51 2007
Дата индексирования: Mon Oct 1 21:26:17 2012
Кодировка: Windows-1251
V.Marsakov, Yu.Shevelev, F-stars in galactic disk

Свойства функции металличности F-звезд в галактическом диске

Ю.Шевелев & В.Марсаков

1998 Bull.Spec.Astrophys.Obs.V.46, P.101

[PostScript-формат, English, ]


Abstract

     На базе восстановленного статистическим методом (основанном на вычислении пространственных орбит) реального распределения по высоте F-звезд на солнечном галактоцентрическом расстоянии построены и исследованы свойства функций металличности (ФМ) звезд диска разного возраста. Показано, что на высоком уровне статистической значимости аппроксимация одной норалью ФМ звезд любого возраста отвергается против альтернативы описания ее суммой двух нормальных кривых с центрами в <[Fe/H]>&rovn;-0.2 и 0.0. При этом максимум распределения звезд моложе &rovn;3 млрд. лет находится вблизи солнечной металличности, а у более старых звезд он попадает в район второго выделенного значения металличности. Сделан вывод, что наблюдаемые в диске отрицательный тренд металличности с возрастом и вертикальный градиент металличности вызваны главным образом перераспределением относительной численности объектов между этими двумя выделенными значениями металличности.

     THE PROPERTIES OF THE METALLICITY FUNCTION OF F STARS IN THE GALACTIC DISC, by Shevelev Yu.G., Marsakov V.A. - We have constructed and investigated properties of metallicity function of disk stars with different ages, proceeding from the statistically restored real altitude distributions of stars at the solar galactocentric distance. It is shown that the metallicity distribution of stars in each age range is not of Gaussian form, but it may be considered as a superposition of two normal distributions centered at <[Fe/H]>&rovn;-0.2 and 0.0 dex. The maximum of the distribution function for stars younger than &rovn; 3 Gyr is located near the solar metallicity, and it falls within the region of the second specific metallicity for older stars. The conclusion is drawn that the observed negative metallicity trend with age and the verical metallicity gradient are mainly due to the redistribution of the relative number of objects between these two specific metallicity values in the Galactic disk.

.

Введение

     Важнейшим источником информации об истории обогащения межзвездной среды являются распределения звезд разного возраста по содержанию в них тяжелых элементов. Благодаря довольно высокой стабильности относительного содержания химических элементов в космических объектах различной природы о содержании тяжелых элементов (т.е. элементов тяжелее водорода и гелия) зачастую судят по содержанию в них железа - элементу, наиболее богатому линиями поглощения в звездах спектральных классов A-G. Величину [Fe/H]=lg(Fe/H)star-lg(Fe/H)&sun;, где (Fe/H)star - отношение числа атомов железа к числу атомов водорода в звезде, а (Fe/H)&sun; - такое же отношение на Солнце, называют металличностью. Распределение звезд по [Fe/H] называется функцией металличности (ФМ). Уже первая построенная ФМ для полутора сотен ближайших к Солнцу звезд обнаружила "парадокс G-карликов" (см. van den Berg, 1962), заключающийся в дефиците малометалличных звезд по сравнению с "простой" моделью химической эволюции Галактики (подразумевающей непрерывные процессы обогащения межзвездной среды и звездообразования и мгновенный круговорот вещества, см., например, Audouze & Tinsley, 1976). Более подробные исследования ФМ примерно по пятистам F5-K5 карликам и субкарликам показали, что распределение тяжелых элементов среди звездных населений Галактики имеет дискретный характер (Марсаков & Сучков, 1977). Оказалось, что объекты гало отделяются в целом от населений диска областью значений [Fe/H]=-(0.5-0.7) dex, где обнаружен существенный дефицит звезд. Такая же структурная деталь (провал на ФМ), но значительно менее выраженная наблюдается у G-карликов диска и красных гигантов, лежащих в направлениях галактических полюсов, в окрестности [Fe/H]=-(0.1-0.2) dex (Марсаков & Сучков, 1980). В итоге звездное население диска оказывается разделенным по признаку металличности на две группы. Примечательно, что группы населений диска с разной металличностью демонстрируют различие и в кинематике и в пространственном положении (Марсаков & Сучков, 1977; 1978; 1980).

     В настоящей работе на основе созданного нами каталога фундаментальных астрофизических параметров близких звезд, содержащего около 5500 объектов (Marsakov & Shevelev, 1995a), исследованы свойства ФМ F-звезд диска: построены распределения по металличности звезд разного возраста и рассмотрена структура ФМ на разных расстояниях от плоскости Галактики. Ни в одной из предыдущих работ совершенно не учитывался тот факт, что в выборке звезд, находящихся в окрестности Солнца и вместе с последним - вблизи галактической плоскости, существует возрастная селекция. Селекция обусловлена более высокой вертикальной компонентой скорости более старых звезд, из-за чего они при движении по орбите лишь незначительную часть времени находятся вблизи плоскости Галактики, а затем удаляются от нее на большие расстояния. Отсюда и доля старых (в среднем менее металличных) звезд в выборке из околосолнечной окрестности будет меньше. В данной работе для учета описанного селекционного эффекта мы использовали разработанный нами в (Марсаков & Шевелев, 1995b) метод статистического восстановления реального распределения по высоте объектов на солнечном галактоцентрическом расстоянии. Не учитывался ранее и другой тип селекции, который связан с металличностью. Действительно, звезды, находящиеся в настоящий момент в окрестности Солнца, родились, вероятнее всего, вблизи апогалактических радиусов своих орбит (Rmax), т.е. на расстояниях зачастую значительно более далеких, чем их нынешнее местоположение. Поскольку же в каждый момент времени в галактическом диске существует отличный от нуля отрицательный градиент металличности (см., например, Шевелев & Марсаков, 1993), мы учли этот градиент (d[Fe/H]/dRmax) и привели значения [Fe/H] используемых здесь звезд к солнечному галактоцентрическому расстоянию.

Исходные данные

     Мы воспользовались нашей выборкой звезд F2-G2 (Marsakov & Shevelev, 1995a), где на основе однородных компилятивных uvby-данных, а также опубликованных положений и собственных движений определены металличности, изохронные возрасты, фотометрические расстояния, тангенциальные скорости и другие параметры примерно для 5500 звезд главной последовательности, лежащих в окрестности 80 пс от Солнца. Примерно для трети звезд, для которых удалось найти в литературе лучевые скорости, вычислены элементы галактических орбит. В работе (Марсаков & Шевелев, 1995b) мы подробно обосновали репрезентативность данной выборки по отношению к F-звездам дисковой подсистемы. Более того, в пределах 50 пс от Солнца выборка в интересующем нас спектральном диапазоне F2-G2 оказалась практически полной (см. гистограмму по расстояниям в Марсаков & Шевелев, 1995b). В той же работе показано, что репрезентативность сохраняется и для выборки звезд с измеренными лучевыми скоростями, которые вместе с тангенциальными компонентами скоростей дают возможность вычислить полные скорости и элементы орбит. Для удобства расчетов мы взяли в данной работе в качестве исходного объема куб со стороной 80 пс и центром в Солнце, внутри которого полнота выборки F-звезд практически сохроняется.

     Для всех используемых здесь звезд мы определили возрасты по Новым Йельским изохронам (Green et al., 1987) с помощью разработанного нами в (Шевелев & Марсаков, 1993) экспресс-метода. Поскольку в данной работе не рассматриваются конкретные значения индивидуальных возрастов звезд (а только средние возрасты звездных групп), мы посчитали возможным использовать формальные определения изохронных возрастов и для звезд, лежащих в непосредственной близости от главной последовательности нулевого возраста (ГПНВ). Все эти звезды попадают только в самую молодую группу и поэтому на физические характеристики, вычисленные по более старым выборкам звезд, они не влияют, а только увеличивают статистическую достоверность результатов у молодых звезд. Анализ показал, что в исследуемом спектральном диапазоне представлены как самые молодые (лежат на ГПНВ), так и самые старые звезды диска - именно в окрестности G2 находятся точки поворота металличных и малометалличных звезд (см. рис. 1 из Marsakov et al., 1990 и обсуждение в (Шевелев & Марсаков, 1993)). В итоге мы получили возможность построить неискаженные ФМ для F-звезд диска любого возраста.

     Любая выборка, составленная из объектов ближайшей околосолнечной окрестности, отягощена кинематическим селекционным эффектом,вызванным преимущественным нахождением здесь более молодых звезд, вертикальные составляющие скоростей которых малы, и орбиты лежат вблизи галактической плоскости. Поэтому для целей настоящей работы корректнее использовать ФМ звезд, находящихся на солнечном галактоцентрическом расстоянии в достаточно высоком вертикальном цилиндре (чтобы внутри его объема попали все объекты используемой подсистемы на данном расстоянии). При этом радиус основания цилиндра должен быть довольно малым, чтобы результат не искажался существующим радиальным градиентом металличности. В статье (Марсаков & Шевелев, 1995b) мы описали разработанный нами метод статистического восстановления по репрезентативной выборке не только полной численности звезд в вертикальном цилиндре на солнечном галактоцентрическом расстоянии, но и их реального распределения по высоте. Понятно, что объем исходной выборки должен быть достаточен для формулирования статистически обоснованных выводов. Идея метода заключается в следующем. Среди звезд, находящихся в настоящий момент вблизи галактической плоскости, имеются представители любой подсистемы Галактики. Мы полагаем, что полная выборка звезд из ограниченной околосолнечной окрестности является репрезентативной по отношению к подсистемам разного возраста в том смысле, что селекционные эффекты, кроме упомянутого выше кинематического, в ней отсутствуют. Поэтому с течением времени их ФМ меняться не будет, хотя в выделенном объеме одни звезды будут заменяться другими. Через некоторое время звезды каждой подсистемы распределятся по своим орбитам и займут тот же объем (по координате Z), который занимают в настоящее время на данном галактоцентрическом расстоянии все остальные звезды, принадлежащие той же подсистеме. Благодаря тому, что меридиональная орбита почти каждой звезды диска за достаточно большое число оборотов представляет собой полностью заполненный "ящик" (см. рис. 4 из (Marsakov & Shevelev, 1995b)), любая звезда при своем движении по орбите рано или поздно попадет в каждую точку перпендикуляра к галактической плоскости, проходящего через нынешнее положение звезды. Длина перпендикуляра ограничивается величиной максимального удаления орбиты звезды от плоскости Галактики, т.е. |Zmax| (см. рис. 4 из (Marsakov & Shevelev, 1995b)). Мы положили, что число звезд на перпендикуляре с такими же, как у данной звезды астрофизическими параметрами (металличностью, возрастом, скоростью и т.д.) равно отношению периода ее обращения вокруг центра Галактики к времени пребывания ее в течение этого периода в пределах ±40 пс от галактической плоскости. Восстановленное распределение звезд подсистемы по высоте получается в результате обработки распределения звезд полной выборки из кубического объема по Zmax вычисленными весами и функцией плотности вероятности. (Последняя операция необходима для получения "случайного" распределения звезды на орбите от-Zmax до +Zmax.) Подробнее о восстановлении по Z см. в (Марсаков & Шевелев, 1995b).

ФМ звезд диска разного возраса

     На рисунке 1 приведены для сравнения распределения по металличности F-звезд диска в нескольких диапазонах возрастов, построенные по всей выборке звезд из сферы радиуса 80 пс вокруг Солнца из (Marsakov & Shevelev, 1995a) (левый столбец) и по восстановленной с помощью описанной выше методике выборке звезд из вертикального цилиндра с квадратным основанием 80 X 80 пс (правый столбец). Радиальный градиент металличности учтен только в последнем случае. Для выборки из 80 пс на оси ординат отложены наблюдаемые численности звезд (n), тогда как для восстановленной выборки приведена доля в процентах от вычисленного полного числа звезд в вертикальном столбе. Из рисунка видно, что между одноименными распределениями имеются систематические структурные различия. Введение весовых коэффициентов увеличило численности звезд с большими вертикальными составляющими скорости. Однако высокоскоростные звезды, как правило, являются и менее металличными и более старыми (Marsakov et. al., 1990). Поэтому ФМ звезд промежуточного (t=(3-5) млрд. лет) и большого возраста (t=(5-12) млрд. лет) несколько расширились в малометалличную сторону (см. соответствующие пары на рисунке 1). Учет же радиального градиента металличности (согласно формуле (1) из (Шевелев & Марсаков, 1995)) увеличивает металличность звезд с большими апогалактическими радиусами орбит, что приводит к небольшому смещению всех ФМ в положительную сторону по [Fe/H] и некоторому перераспределению объектов внутри каждой. Кривые в левом столбце на рисунке представляют собой аппроксимации гистограмм полиномами десятой степени. Высокая степень полинома способствует выявлению структурных деталей на распределениях. Из рисунка видно, что уже исходные ФМ во всех возрастных диапазонах демонстрируют, хотя и одновершинные, но довольно вытянутые в малометалличную сторону распределения, которые одной гауссианой не описать. После же учета показанных выше селекционных эффектов оказалось, что начиная с возраста &rovn;3 млрд. лет ФМ приобретает сложную структуру в центре распределения.

Рис. 1.—Распределения по металличности F-звезд диска в разных диапазонах возрастов: левый столбец - для выборки звезд в окрестности 80 пс от Солнца; правый столбец - для восстановленной полной выборки звезд на солнечном галактоцентрическом расстоянии с учетом радиального градиента металличности. Цифрами указаны диапазоны возрастов в млрд. лет. В правом столбце по оси ординат отложен процент от восстановленного полного числа звезд. Кривые на рисунках представляют собой аппроксимации гистограм полиномами десятой степени (слева) и суммами двух гауссиан (справа).

     Остановимся немного на методе, позволяющем оценить статистическую значимость аппроксимации гистаграмм суммой нормальных распределений. Представим ФМ на рис. 1 суммой гауссовых кривых:

&form1;

где m=[Fe/H], ak и &sigma_k; - соответственно среднее значение и дисперсия k-ой нормали, &alfa_k; - нормировочный параметр k-ой нормали, удовлетворяющий условию &form3;, N - полное число звезд. Параметры распределения находятся методом максимального правдоподобия

&form2;

где l - число классовых интервалов в исследуемой гистограмме, ni -число звезд в i-ом классовом интервале, mi-значение металличности в середине i-го интервала. Величину s можно определить, используя статистику

&form4;

где Ls и Ls+1 берутся в точке максимума. Известно, что величина &lamda; распределена по закону &ksi;2 c одной степенью свободы (Martin, 1971). Нас в первую очередь интересует вопрос, можно ли описать ФМ нормальным законом (нулевая гипотеза H0: s=1), или же распределение металличности имеет более сложную структуру (альтернативная гипотезаH1: s>1). Расчеты проведем для области -0.5<[Fe/H]<+0.2. Определенные методом максимального правдоподобия параметры гауссиан и соответствующие значение &lamda; показали, что вероятность ошибочно отвергнуть гипотезу об описании распределения гауссовой кривой против альтернативы представления ее суммой двух гауссиан для звезд промежуточного возроста оказались <12%, а для самых старых звезд <13% < 12%. (При пользовании методом максимального правдоподобия в данной работе для каждой гистограммы мы принимали в формуле (1) в качестве нормирующего множителя N полное наблюдаемое число звезд, использованное для построения данного распределения, а не полное восстановленное, полученное из первого путем обработки его весами. Тем самым мы ушли от неоправданно завышенных величин вероятностей.) Очень высокими в указанном диапазоне [Fe/H] оказались значимости описания суммами двух гауссиан двух верхних гистограмм правого столбца. Это связанно с существованием у них протяженных малометалличных "хвостов". При сужении диапазона металличности с малометалличной стороны гипотизу об описании распределений одной гауссианой отбросить уже нельзя. Нас, однако, интерисует именно весь диапазон металличности, характерный для диска. И мы видим, что в его пределах ФМ звезд любого возроста биноминальна. В то же время статистическая достоверность гипотизы s>2 оказывается во всех случиях крайне низкой. Аппроксимация распределений суммами двух гауссиан приведены в правом столбце рисунка 1 сплошными кривыми линиями.

Связь между возрастом и металличностью

     Рассмотрим как выглядит связь между возрастом и металличностью в галактическом диске, полученная по полной выборке F-звезд в ограниченном объеме околосолнечной окрестности после учета кинематического селекционного эффекта и радиального градиентов металличности. На рис. 2 приведена зависимость t-[Fe/H], для построения которой наша выборка F-звезд из кубического объема ±40 пс с корректно определенными изохронными возрастами (т.е. использовались только звезды уже слегка проэволюционировавшие от ГПНВ начиная с δMV>0.3m) была разделена по возрасту на 8 диапазонов, равных 1 млрд. лет (кроме самого старого диапазона, который из-за малой численности пришлось увеличить до 4 млрд. лет). Барами нанесены дисперсии металличности в каждом возрастном диапазоне. Ошибки определения средних значений [Fe/H], найденные по реальным численностям (т.е. без учета весовых коэффициентов, которые увеличивают первоначальную численность) лежат в диапазоне (0.01-0.03) dex и почти везде меньше различия между соседними группами звезд. Как видно из рисунка, тренд металличности с возрастом звездной группы хорошо ложится на прямую линию: прямая регрессия имеет вид [Fe/H]=-(0.05±0.03)t+(0.11±0.02) при коэффициенте корреляции r=0.98±0.01. Штриховой линией для сравнения нанесена аналогичная зависимость, построенная по F-звездам выборки с корректно определенными возрастами, лежащим вокруг Солнца в сфере с радиусом 80 пс (см. Шевелев & Марсаков, 1993). В последнем случае никакой коррекции за кинематическую селекцию и градиент металличности не делалось. Видно, что зависимости лишь нескоько сместились параллельно самой себе, оставив наклон неизменными. Качественно сохранились и описанное в работе (Шевелев & Марсаков, 1993) повидение дисперсии металличности от возроста, т.е. с увеличением возроста звездной группы до &rovn;7 млрд. лет разброс [Fe/H] у звезд данного возроста увеличивается, а затем снова уменьшается. Это можно проследить по уменьшающимся в обе стороны от этого возраста величинам баров дисперсий на рис. 2 и сравнить с разбросом звезд по [Fe/H] в узких возростных диапазонах на диаграмме возраст-металличность на рис. 4 из (Шевелев & Марсаков, 1993).

Рис. 2.—Зависимость среднего значения металличности от возраста звезд, построенная по восстановленной полной выборке с учетом радиального градиента металличности. Бары равны дисперсиям металличности для соответствующего возрастного бина. Величина центрального кружка равна средней ошибке вычисления значения [Fe/H]. Сплошная линия - прямая регрессия для полной восстановленной выборки; пунктирная - для исходной выборки F-звезд из окрестности 80 пс от Солнца (без учета радиального градиента металличности и весовых множетелей).

ФМ звезд диска на разной высоте над галактической плоскостью

     Рассмотрим изменение формы ФМ F-звезд с высотой Z над плоскостью Галактики для всей выборки в целом и для трех узких возрастных диапазонов по отдельности: (t<12) млрд. лет, (t<3) млрд. лет, (3<t<5) млрд. лет и (5<t<12) млрд. лет. С этой целью мы построили распределения по металличности звезд в полосах шириной 10 пс из восстановленной в вертикальном столбе полной выборки через каждые 50 пс вплоть до Z=250 пс. (На более далеких Z даже у полной выборки не хватаетчисленности объектов для получения статистически достоверных результатаов.) В качестве иллюстрации на Шевелев & Марсаков (1995). Сплошными кривыми на рисунке нанесены аппроксимации гистограмм суммой двух гауссиан. По оси ординат отложена доля звезд в процентах от полной численности восстановленной выборки, ФМ которых приведены в правом столбце Марсаков & Шевелев (1995).) Анализ показывает, что самые молодые звезды вплоть до 200 пс уверено демонстрируют одновершинные распределения с максимумом вблизи [Fe/H]&rovn;0.0, но с не большим "хвостом" в малометалличную сторону, увеличивающимся с ростом Z. (На более далеких Z численность самых молодых звезд резко падает и на структуру ФМ начинают оказывать заметное влияние выбросы на гистограммах, связанные с высокими весовыми коэффициентами отдельных звезд, как это видно, например, на Z=250 пс.) Именно этот "хвост", как показывают оценки максимального правдоподобия, и обуславливает невозможность (уровень значимости <5%) описания распределений при Z=0 и 100 пс одной гауссианой против альтернативы аппроксимации их суммой двух гауссовых кривых, приведенных на рисунке. Если же ограничить рассматриваемый диапазон металличности интервалом [-0.25 - +0.20], (т.е. отбросить классовые интервалы в хвостах с малой численностью звезд), то окажется, что центральная часть этих распределений лучше описывается одной гауссианой. У звезд промежуточного возраста структура ФМ более сложная и, начиная с Z=50 пс, распределение становится явно бинормальным. В частности, аппроксимация приведенного на рисунке 3 распределения для 100 пс одной гауссианой отвергается по критерию отношения правдоподобий в пользу альтернативы описания ее суммой двух гауссиан на уровне значимости <17% (см. пунктирную кривую на соответствующей гистограмме рис. 3). Но уже на Z = 250 пс ФМ вновь принимает одновершинный характер (вероятность ошибочно отвергнуть гипотезу об описании ее одной гауссианой ~35%). Примечательно, что максимум распределения смещается здесь в малометалличную сторону примерно на -0.2 dex (по отношению к максимуму на ФМ самых молодых звезд), оставляя некоторый избыток объектов на металличном крыле. В верхнем ряду рис. 3). На более далекихрасстояниях именно эта группа начинает доминировать (на рисунке не приведено). В результате с ростом Z средняя металличность ФМ уменьшается, а дисперсия металличности несколько увеличивается.

Fig. 3.—Распределения по металличности F-звезд восстановленной выборки в разных диапазонах возрастов на трех расстояниях от галактической плоскости. Кривые - аппроксимации гистограмм суммами двух гауссиан.

     Отметим, что точно также меняют форму с высотой поднятия над галактической плоскостью и ФМ, построенные по "звездам" восстановленной выборки без учета радиального градиента металличности, только все они оказываются смещенными в малометалличную сторону на d[Fe/H]dR&rovn; -(0.02-0.04) dex. Таким образом, свойства ФМ F-звезд свидетельствуют, что существует по крайней мере два выделенных значения металличности, которыми преимущественно обладают звезды дисковой подсистемы в окрестности Солнца.

Вертикальный градиент металличности

     В работе (Марсаков & Шевелев, 1995a на основе полной выборки из Шевелев & Марсаков, 1995 мы определили величины вертикальных градиентов металличности звезд диска в нескольких узких возрастных диапазонах по наблюдаемым значениям металличности и вычисленным максимальным удалениям их орбит от галактической плоскости. Было обнаружено, что величины градиентов систематически уменьшаются с переходом к более старым звездам. Представляется важным проверить, что станет с вычисляемой величиной вертикального градиента и как изменится характер ее эволюции со временем после учета обсуждаемого в п.2 кинематического селекционного эффекта. С этой целью мы разделили все звезды восстановленной выборки на 3 возрастных диапазона и каждую на 6 групп по Z (как и в п.5) и построили для них зависимости средних значений металличности на данной высоте от самой величины Z. Эти зависимости для трех узких возрастных диапазонов, а также для звезд смешанных возрастов (т.е. для всех звезд восстановленной выборки) приведены на рис. 4. Барами обозначены ошибки определения <[Fe/H]>, вычисленные по реальным численностям звезд (т.е. без учета весовых коэффициентов). Видно, что вс