Введение к работе
актуальность темы. Мазерные источники НгО могут быть двух типов. Одни вязаны с областями, где идет процесс звездообразования, а другие, напротив, вязаны со звездами поздних спектральных классов, которые находятся на авершающей стадии своей эволюции.
I настоящей работе приведены результаты исследования мазерных источников, оторые ассоциируются с областями активного звездообразования. Интерес ызван тем, что наблюдения мазеров НгО в этих областях являются в настоящее ремя единственной возможностью изучения поведения звезд на ранней стадии їх эволюции, поскольку они в этот период не видны в оптическом диапазоне з—за сильного поглощения в газо-пылевой оболочке, окружающей звезду, 'егулярные наблюдения в течение длительного времени позволяют выявить юлгопериодические процессы в нестационарном характере аккреции вещества а формирующуюся звезду. Согласно теоретическим исследованиям не-тационарная аккреция вещества может иметь периодический характер с ериодом 5-10 лет [1, 2]. Несмотря на то, что наблюдения мазеров НгО, вязанных с областями звездообразования, проводились неоднократно, они іьіли выполнены для небольшого числа источников и, как правило, в течение 'граниченного периода времени (обычно несколько лет) [3, 4]. 'формулируем наиболее актуальные задачи исследования мазерных істочников НгО в областях звездообразования.
. Изучение переменности полного излучения спектра НгО в течение периода юлее 10 лет. Выделение компонентов с различным характером и разными іериодами переменности.
!. Исследование эволюции отдельных эмиссионных деталей в спектрах НгО. Ізмерение вариаций потока, лучевой скорости и ширины линии, і. Исследование мазеров НгО как физических объектов. Здесь можно выделить ри аспекта проблемы:
-
насыщенное и ненасыщенное состояния мазеров;
-
явление антикорреляции потоков различных деталей в спектрах НгО в лучае частично насыщенного мазера при переменном темпе накачки [5, 6],
-
место локализации области мазерного излучения в оболочке [7].
L Построение моделей мазерных источников - определение характера труктуры, в которой располагаются мазерные конденсации (кеплеровский диск їли оболочка), ее размера, скорости расширения и вращения, массы іентральной звезды.
1ель работы. Первой задачей работы являлись разработка и создание іногоканального анализатора спектра с высоким спектральным разрешением: ),1 км/с для получения полных спектров мазерных источников НгО и 34 м/с іля исследования отдельных, наиболее узких спектральных деталей. Другая задача состоит в проведении регулярных наблюдений большого числа істочников и получении каталога спектров с высоким спектральным разре-
шением, а также в разработке метода анализа этих данных
Третья задача состоит в интерпретации полученных наблюдательных даннь
(построение моделей мазерных источников НгО).
Новизна. Новыми являются следующие результаты, полученные в диссертаци
-
С помощью 96-канального анализатора спектра на радиотелескопе РТ-2 РАС ФИАН в Пущино получен каталог спектров НгО с 1979 года (дл большинства источников с 1980 г.) с высоким спектральным разрешением.
-
Разработан и применен новый метод обработки данных, позволяют^ определить основные параметры оболочек, на которых находятся мазерны конденсации (скорости расширения и вращения, размер).
-
Обнаружены два компонента переменности светимости звезд в начально стадии их формирования: долгопериодический с характерным периодом 4-1 лет (это согласуется с результатами теоретических исследований) и вспышеч ный, повторяющийся каждые 2-4,5 года. Существование вспышечного компо нента не было теоретически предсказано.
-
Интерпретация мазерного излучения НгО источника S140 с триплетно структурой спектра отождествлением его с кеплеровским диском (протопла нетным диском). Наблюдаемый характер переменности потока и лучево скорости компонентов спектра объяснен движением мазерных конденсаци (возможных протопланет) по кеплеровским орбитам.
-
Впервые наблюдался процесс перехода мазера из ненасыщенного состояии в насыщенное и обратное. Этот процесс сопровождался изменением все характеристик излучения - потока, ширины линии, ее формы и лучевой скорості
6. Показано, что' в оболочечных структурах мазерных источников явлени
антикорреляции также может проявляться, но при более жестких условиях, чеі
в кеплеровских дисках.
Научная и практическая значимость. Наиболее ценными результатам диссертации нам представляются: созданный анализатор спектра, которы функционирует в составе штатной аппаратуры на радиотелескопе РТ-22 РА< ФИАН в Пущино; полученный каталог спектров НгО мазерных источникої методы обработки экспериментальных данных; результаты интерпретаци полученных данных Научные результаты, полученные в диссертации, могу/ быть использованы для развития теории космических мазеров. Методы анализа, изложенные в диссертации, будут применены для анализ данных по другим источникам из нашей программы долговременных наблюдений мазеров НгО. Представляет научный интерес построение моделе; большого числа источников с последующей их классификацией на основі эволюционного статуса.
Апробация. Результаты работы докладывались на Всесоюзных конференция по "Галактической и внегалактической радиоастрономии", на научных семинара в ГАИШ, ФИАН, РАС ФИАН в Пущино и на совещаниях рабочей rpynnt "межзвездная среда".
се основные результаты опубликованы в периодических астрономических зданиях и сборниках
а защиту выносятся:
. Разработка и создание многоканального анализатора спектра с разрешением
,1 км/с и 34 м/с.
, Разработка нового метода анализа спектров путем выделения отдельных
эмпонентов и последующего построения различных гистограмм, связывающих
ва из трех полученных параметров: поток в максимуме излучения, лучевую
іазе излучения.
. Результаты исследования мазерных источников НгО ОІМ1, W31A, S252A,
/75N, W75S и S128, отождествляемых с газо-пылевыми оболочками.
. Результаты исследования мазера с триплетной структурой спектра в S140,
гождествляемого с кеплеровским диском.
. Результаты исследования мазерных источников как физических объектов.
труктура и объем диссертации. Диссертация состоит из введения, шести іав, заключения и списка цитируемой литературы из 87 наименований. Общий бъем диссертации составляет 189 "страниц, включая 51 рисунок и 2 таблицы.
о введении показана актуальность проблемы звездообразования и роль сследования мазерного излучения НгО в решении этой проблемы. Дана 5основанность постановки данной работы и сформулированы некоторые, аиболее актуальные задачи.
первой главе дано описание аппаратурного комплекса диапазона 1,35 см, помощью которого были проведены регулярные наблюдения мазерных ггочников НгО [8\ Наблюдения были выполнены нами на радиотелескопе Т-22 Радиоастрономической станции ФИАН в Пущино в период с 1979 по 995 г. Антенна на волне 1,35 см имеет угловое разрешение 2,6'. Полная іумовая температура системы составляла 200-250 К. ходным усилителем служил малошумящий парамагнитный усилитель бегущей элны. Приемник супергетеродинного типа с тройным преобразованием частоты, еобходимая стабильность частоты достигается применением в качестве эдающего перестраиваемого генератора синтезатора частоты 46-31. В /істеме умножения частоты используются три кольца фазовой автоподстройки зстоты. С помощью 46-31 осуществляется настройка приемника на нужный пектральный участок
нализ сигнала проводился 9б-канальным анализатором спектра фильтрового та, разработанного и изготовленного автором. Спектральное разрешение эставляет 7,5 кГц, или 0,101 км/с по лучевой скорости. Полная полоса анализа
составляет 9,7 км/с. Для анализа отдельных узкополосных пиков излучени применялся 32-канальный анализатор с разрешением 2,5 кГц (34 м/с). Для наблюдений источников используется режим диаграммной модуляциі осуществляемой с помощью ферритового переключателя, при этом ко вход приемника происходит поочередное подключение сигнала сначала одноп облучателя антенны, а затем другого. Частота переключений около 400 П. Управление работой телескопа и радиометра, а также регистрация сигнал, осуществлялись с помощью ЭВМ.
Для реализации необходимой чувствительности время записи каждой спектрограммы составляло 15-30 минут. Антенная температура в 1 К для точечной источника с неполяризованным излучением соответствует плотности потока 2! Ян. Наблюдения каждого источника в нашей программе почти всегд, проводились на одинаковых позиционных углах Это необходимо для избежани: дополнительного фактора переменности сигнала при существовании линейної поляризации излучения НгО.
В разделе 1.4 этой главы описана методика, используемая нами при проведенні программы многолетних наблюдений мазерных источников НгО. Мазерные источники имеют большой динамический диапазон значений потокг Для наблюдений не сильных источников используется метод ГШ - ON - ОГ* Во время циклов ГШ - ON телескоп наводится так, что сигнал от источник, поступает в приемник через правый облучатель, а от области сравнения чере левый. Калибровочный сигнал вводится в тракт, связанный с левым облучателег* Во время второго цикла ON сигнал от источника поступает в приемник чере левый облучатель. Для сильных источников на РТ-22 используется метод наведения ГШ - OFI
- ON, поскольку нельзя проводить калибровку на фоне достаточно сильной
сигнала от источника. Во время первых двух циклов телескоп наводится ні
области сравнения левым и правым облучателями. Затем телеском наводите:
так, чтсг'сигнал от источника поступает на вход приемника через левьп
облучатель. В этом методе Обычно берется равное время накопления для все
трех циклов наблюдений.
В этом разделе также рассмотрены факторы, влияющие на точность калибровкі сигнала. Основными из них являются:
точность наведения антенны на источник,
влияние погодных условий,
наличие линейной поляризации принимаемого излучения. Проведен тщательный анализ этих факторов.
Во второй главе описана достаточно важная часть работы - методы анализ; и обработки данных Для восьми источников, включенных в диссертацию, бьілі получено свыше 800 спектров. Каждый спектр обрабатывался как целое і как суперпозиция отдельных "пиков излучения". В первом случае бьілі вычислены значения полных потоков путем интегрирования сигнала по всем; спектру, а в некоторых случаях и значение средневзвешенной лучевой скорое определяемой соотношением
Vo^FjViAV/ZFikV,
і і
\e Fj есть значение потока* на лучевой скорости // , a AV - шаг итерирования.
ряде случаев для выявления спектральных интервалов мазерной активности эти получены средние спектры источников как для всего периода наших эблюдений, так и для некоторых периодов времени, в течение которых эоисходили какие-либо важные процессы в эволюции мазерного источника, утем наложения всех спектров (для OIN1 использован 21-ин спектр, которые эти равномерно распределены по временной шкале) на один график для N1, S128 и S252A получена суперпозиция спектров, где наибольшее сгущение іний дает наивероятнейший спектр излучения данного источника, ля отдельных "пиков излучения" были определены следующие параметры: I поток и его вариации со временем; I лучевая скорость;
I дрейф по лучевой скорости, характер дрейфа; I ширина линии (по уровню 0,5 по мощности), ее временные изменения; I форма линии;
I функциональная связь между какими-либо параметрами линии; например, ежду интенсивностью и шириной линии, между интенсивностью и лучевой соростью; ) время нахождения в активной фазе излучения.
разделе 2.4 описан разработанный нами метод (для источников с достаточно южными спектрами), основанный на статистическом анализе данных По всем іектрам путем вписывания гауссовских профилей были выделены "пики тучения" и для каждого из них измерены три параметра - поток в максимуме шучения, лучевая скорость (в случае дрейфа компонента по лучевой скорости іачение скорости бралось в максимуме излучения) и время нахождения энного компонента в активной фазе, когда постоянно без перерывов эблюдалось его излучение. По этим трем параметрам для каждого источника, Зработанного по данной методике, были построены по три гистограммы и три іафика. Для построения гистограммы распределения числа компонентов по /чевой скорости весь интервал лучевых скоростей, в котором наблюдалось азерное излучение НгО данного источника, был разбит на отдельные эследовательные участки протяженностью 1 км/с. Для каждого такого участка это подсчитано число компонентов, которые были выделены за весь период эших наблюдений.
істограммьі распределения числа компонентов по потоку и по времени жизни ыли получены аналогично с шагом ДІ nF = 0,4 и Дг = 1 мес. соответственно, писанный метод статистического анализа данных позволяет учитывать все те азерные конденсации, которые за период наших регулярных наблюдений хоть жое-то время находились в активной фазе. Полученные гистограммы ззволяют подобрать для источника вполне определенную модель. Вид стограммы зависит от спектра источника, который в свою очередь опреде-
ляется структурой, где располагаются мазерные конденсации. Самой важно на наш взгляд, гистограммой является распределение числа компонентов г лучевой скорости. Одновременный анализ такой гистограммы с другил графиками, со средними спектрами и т.д. позволяет определить конкретнь параметры оболочки (тороида): скорости расширения и вращения. Третья глава посвящена исследованию трех источников (ONI, W31A и S252/ для которых обнаружен более менее периодический характер переменное мазерного излучения НгО.
Источник ON1 связан с достаточно молодой областью, где только начинает! процесс звездообра.зования и исключены "спусковые" механизмы звездообрг зования (столкновения облаков, ударные волны и т.д.) [9]. Особенностью спеют. Н?0 ON1 является разделение пиков излучения на две основные группы: і 9 до 12 и от 13 до 18 км/с. Эмиссионные пики этих двух групп имеют разнь характер переменности. Наибольший дрейф по лучевой скорости имеют пиі второй группы спектра. Согласно модели Стрельницкого [10] эти групп отождествляются с левой и правой частями оболочки (кольца, тороид; Расстояние между центроидами этих групп равно удвоенной скорости вращені кольца. Отсюда имеем Увращ = 2,5 км/с. Более слабое излучение на скорост? от 3 до 8 и от 20 до 22 км/с отождествляется с центральными частяг* оболочки. Следовательно, скорость расширения равна ~8 км/с. Таким образом, для каждой мазерной конденсации можно написать уравнен движения. Используя найденную из наблюдений величину изменения лучевс скорости мазерной конденсации вычисляются угловая .скорость и соответств> ющее ей расстояние (-10 см). При условии кеплеровского движенр конденсации для массы центральной звезды получим значение ~6Ме. Кривая изменения полной энергии в спектре НгО ON1 содержит два компонент долгопериодический с почти синусоидальным характером переменности (перис -10 лет) и вспышечный с очень небольшой амплитудой и отсутствие какой-либо периодичности. Такой характер переменности полного поток по-видимому, связан с незначительными изменениями светимости центральне звезды в ON1.
Источники W31A и S252A также имеют два компонента переменности полної потока. Долгопериодический компонент для обоих источников имеет одинаковь период, равный около б-ти лет. Вспышки в W31A повторялись каждые 4-4 года, а в S252A в два раза чаще. Для S252A нами найдена составляющг с предположительным периодом порядка 30-40 лет. Также обнаружь достаточно регулярный и существенный дрейф средневзвешенной лучевс скорости (3 км/с) за период с 1977 по 1994 г. Наиболее подходящей модель для обоих источников является расширяющийся и вращающийся тороид. Ь исключена также модель прецессирующего околозвездного диска для S252 В четвертой главе рассмотрены мазерные источники, связанные с взаимс действующими молекулярными облаками. Внутри таких облаков возникающг ударная волна может привести к неустойчивости с возникновением отдельнь фрагментов повышенной плотности. Взаимодействие облаков может являться "спусковым" механизмом процесс
звездообразования [10]. Наличие градиента скорости в облаке может оказать существенное влияние на характер переменности мазерного излучения. Например, в W75N излучение в диапазоне скоростей 20-28 км/с всегда наблюдалось в режиме "включено-выключено". В S128 мазерное излучение пиков, близких по скорости к облаку с большим градиентом скорости, имеет более сильные вариации потока, чем остальных пиков. Наиболее интенсивное излучение всех трех источников происходит на скоростях, близких к скоростям своих молекулярных облаков, т.е. мазерные источники имеют небольшие пекулярные скорости относительно молекулярных облаков. Проведенные нами исследования показали, что в целом мазерное излучение в W75S имеет три вида переменности:
-
допгопериодические изменения полного потока,
-
вспышки групп деталей,
-
вспышки одиночных деталей.
Первые два вида переменности отражают состояние мазера как системы и могут быть следствием изменения светимости центральной звезды на ранней стадии ее эволюции. Долгопериодическая переменность проявлялась как чередование максимумов и минимумов активности с характерным временем от 2-х до 4-х лет. Не исключено, что она является суперпозицией двух составляющих с разными периодами переменности (4 и 11 лет). Вспышечная активность имела более менее периодический характер с интервалом между вспышками -2,5 года и проявлялась она независимо от фазы долгопериодической составляющей, но интенсивность вспышек мазерного излучения все же коррелировала с ней.
В периоды активной стадии мазера происходили одновременно изменения потока и лучевой скорости некоторых деталей спектра. В рамках оболочечной .модели мазерного источника это объясняется локальными изменениями скорости расширения оболочки в периоды наибольшей активности центральной звезды.
В некоторые периоды времени наблюдалась антикорреляция потоков между группами деталей, либо между отдельными деталями спектра НгО. В этих случаях более слабый компонент спектра всегда подавлялся более сильным. Аналогичное явление наблюдалось и в W31A, поэтому дальнейшее изложение материала будет касаться обоих источников. Антикоррепяция постулируется неоднородностью условий в среде и может быть связана с анизотропией накачки мазера формирующейся звездой, либо с неоднородностью плотности, температуры, содержанием активных молекул. В кеплеровском диске (S255, S140) антикорреляция потоков между компонентами осуществляется в непрерывной среде.
В случае мазера, в котором отдельные конденсации образуют кольцо или тороид, среда не является непрерывной, но она содержит достаточное число сгустков вещества, в которых при благоприятных условиях может возникнуть мазерное излучение. В такой "квазинепрерывной" среде антикорреляция потоков происходит при более жестких условиях, чем в кеплеровском диске. Сначала во время роста активности звезды происходит активизация среды, а затем в
период вспышечной активности звезды могут возникнуть благоприятные условия как для мазерного излучения отдельных конденсаций НгО, так и для процесса конкуренции хордовых мод за накачку в случае частично насыщенного мазера. Для источников W75N и S128, как и в случае всех других источников, рассмотренных в диссертации, также было обнаружено существование двух компонентов переменности мазерного излучения - долгопериодического и вспышечного. Долгопериодическая составляющая в W75N имеет довольно сложный вид. Скорее всего в этом источнике происходит чередование максимумов и минимумов активности с временным интервалом между максимумами около б лет. В S128 эта составляющая напоминает синусоиду с периодом 7 лет. Вспышечная активность в обоих источниках проявлялась нерегулярно.
На основании разработанного нами метода статистического анализа данных построена модель мазеров НгО в W75S и W75N в виде вращающегося и расширяющегося кольца. Сделана оценка массы центральной звезды. В направлении источника S128 наблюдается очень сложная структура области, где компактная зона НИ и мазерный источник НгО расположены вдоль границы раздела двух взаимодействующих облаков. В результате этого возможно сильное отклонение в расположении конденсаций от оболочечной структуры. Не исключено также, что мазерные конденсации расположены в двух, близких гнездах с единой возбуждающей звездой. Каждое гнездо может быть связано в большей степени с одним из облаков и в меньшей степени с другим. В пятой главе приведены результаты исследования мазерного источника G43.8-0.1 как астрономического, так и физического объекта. Мазер НгО находится на расстоянии -7" от сверхкомпактной области НИ [121 Многочисленные мазерные конденсации образуют оболочку радиусом -10 см [13} Мазер НгО G43.8-0.1 уникален тем, что в спектре имеется эмиссионная деталь на лучевой скорости 42.2 км/с, излучение которой наблюдается постоянно с момента открытия мазера G43.8-0.1 в 1976 г. В среднем один раз в три года происходили сильные вспышки излучения отдельных деталей спектра. Таких вспышек наблюдалось пять и все они не были похожи друг на друга. Самая узкая деталь имела минимальную ширину 0,36 км/с при амплитуде 3500 Ян (конец 1981 г. - начало 1982 г.). Найдена корреляция между потоком и шириной линии. Самая мощная вспышка (4500 Ян в максимуме) произошла в начале 1987 г. на лучевой. скорости 40 км/с. В течение всего периода эволюции этой детали ширина линии росла с 0.85 до 0,95 км/с. Излучение компонента 42,2 км/с имеет ряд особенностей по отношению к другим: долгое время нахождения в активной фазе, медленные изменения потока со временем, незначительные изменения лучевой скорости. Для этого компонента найдена функциональная связь между шириной линии и потоком излучения на протяжении достаточно длительного интервала времени. По характеру связи можно выделить три разных участка. На первом происходит сужение линии с 0,70 до 0,55 км/с при непрерывном росте потока: на втором участке сохраняется ширина линии на уровне -0.55 км/с при дальнейшем росте потока: наконец, на третьем участке происходит уширение линии при
уменьшении потока. Последующие вариации ширины линии и плотности потока имели очень сложный характер.
Более тщательный анализ позволяет рассмотреть шесть стадий эволюции мазерной конденсации 42,2 км/с и определить степень насыщения мазера. Это можно сделать по исследованию формы линии (ее крыльев). При достаточном коэффициенте усиления среды (г = 35 в центре, линии) отличие потока на 8% от гаусианы в крыле линии может быть аргументом в пользу насыщенного состояния мазера. За весь период наших наблюдений с 1981 по 1994 г. дважды происходил переход мазера из ненасыщенного состояния в насыщенное и обратно.
Как и для всех других источников, рассмотренных в диссертации, нами было получено, что кривая переменности полного потока состоит из двух компонентов, которые отражают вспышечный и долгопериодический характер процесса формирования центральной звезды.
Вспышечный компонент полного потока проявлял себя необычно. Резкий скачок его вверх в конце 1986 г. был вызван серией последовательных вспышек нескольких деталей в течение чуть более трех лет. Одиночные сильные вспышки мало повлияли на изменения полного потока. С учетом этого кривая, отражающая долгопериодическую переменность, имеет достаточно простой вид На ней имеются два минимума с интервалом между ними 8-9 лет и один максимум. Период наблюдений этого источника, 13 лет, оказался недостаточным, чтобы утверждать о существовании строго периодических процессов мазерной активности в G43.8-0.1.
Шестая глава. Мазерное излучение НгО в S140 сосредоточено в трех узких интервалах лучевых скоростей, разделенных примерно одинаковыми промежутками. Такие триплетные структуры могут возникать в околозвездных кеплеровских дисках, ориентированных к наблюдателю ребром. Нами наблюдалась сильная переменность интенсивностей, лучевых скоростей и ширин линий отдельных компонентов. Также наблюдалась слабая антикорреляция интенсивностей боковых компонентов и более четко выраженная антикорреляция интенсивностей центрального и боковых компонентов. Одновременно все три компонента наблюдались относительно редко - лишь в 1977, 1984-1985 и 1987 гг. В остальное время либо два компонента, либо только один. Во всех случаях, когда какой-либо из трех компонентов достигал своего "абсолютного" максимума, оба других отсутствовали. Если в качестве модели источника принять однородный кеплеровский диск, то боковые компоненты спектра НгО должны иметь одинаковые интенсивности. Однако, наблюдения S255 [5] и S140 показали, что в действительности это не так и, более того, наблюдается антикорреляция интенсивностей. Для объяснения антикорреляции необходима неоднородность условий в диске (анизотропия накачки мазера центральной звездой). Анизотропия накачки может быть заменена какой-либо другой неоднородностью физических условий в диске. Уточнение характера неоднородностей в диске связана с необходимостью объяснить наблюдавшиеся вариации лучевой скорости и ширины компонентов спектра. Любая неоднородность плотности с размерами, превы-
шающими несколько процентов от длины хорды и плотностью, заметно отличающейся от средней (т.е. способная исказить общую оптическую длину на несколько процентов) может заметно изменить усиление вдоль этой хорды. Мы наблюдали три типа поведения лучевой скорости левого компонента, которые могут быть объяснены тремя траекториями движения плотной конденсации. Входя в "коридор" усиления такие конденсации увеличивают оптическую длину пути усиления. По наблюдаемому темпу изменения лучевой скорости были вычислены периоды обращения конденсаций, расстояния до звезды и масса звезды (<5Мо).
Аналогичные вычисления были выполнены и для центрального компонента в период его активности.
Эволюция правого компонента была более сложной. В интервале скоростей 4.4-9.7 км/с наблюдалось последовательное появление субисточников на непродолжительное время. Все они имели дрейф по лучевой скорости. У одних скорость возрастала, у других падала. Оказалось, что в группе "восходящих" субисточников лучевая скорость в целом также возрастает, а в группе "нисходящих" - убывает. Описанная картина объясняется наложением на луче зрения двух скоплений газовых сгустков, которые могли бы быть протоспутнико-выми системами двух гигантских протопланет.
В заключении перечислены основные результаты работы, сделаны некоторые обобщения и намечены основные направления будущих исследований.