Сейфертовские галактики представляют собой редко встречающийся тип галактик, обладающих активным ядром. Светимость их яркого звездообразного ядра иногда превышает 10% от интегральной светимости галактики в видимой области спектра. Ядро характеризуется  необычайно широкими эмиссионными линиями, свидетельствующими о движении газа со скоростями в тысячи км/с.

В настоящей задаче исследуется спектрограмма ядра сейфертовской галактики, и по относительным интенсивностям спектральных линий определяются физические параметры излучающего газа.

Задача выполняется на ПК, работающем под Linux (операционная система, подобная unix), в графической оболочке X-windows с использованием системы обработки изображений MIDAS. Краткий обзор собственных команд Linux, которые могут быть полезны в работе, а также команд пакета MIDAS, приведены в Приложении.

Исходные спектры для задачи получены на фотопластинке, а затем оцифрованы. Поэтому обработку мы будем проводить с учетом особенностей фотоэмульсии как приемника излучения.

Для выполнения задачи необходимо иметь 4 файла: обработанные на автомикроденситометре спектрограммы галактики NGC 1068 (n5.bdf либо n16.bdf), звезды А0 (a0.bdf), фотографическую шкалку (klin.bdf), а также файл–таблицу (reference.tbl). Перед выполнением задачи обязательно ознакомьтесь с описанием простейших команд Linux и MIDAS, приведенных в Приложении.

1.            Войти в Linux с паролем, сообщенным преподавателем. Командой startx запустить X‑Windows [это может быть сделано преподавателем заранее]. Создать рабочую директорию (желательно, чтобы ее имя содержало фамилию и дату), в которой будут храниться рабочие файлы [например, mkdir ivan\_23022001]. Файлы, необходимые для работы, содержатся в директории initial\_data/11/. Cкопировать 4 файла в директорию, в которой будет выполняться работа, командой

cp/initial\_data/11/* /ваша\_рабочая\_директория.

В X–Windows открыть окно и запустить MIDAS командой inmidas -p 0x , где x — цифра от 0 до 9. Данная команда позволяет запускать одновременно несколько сессий MIDAS на одной и той же машине.

2.            Плотности почернения, получаемые на фотопластинке, нелинейно зависят от упавшего на них потока. Поэтому для калибровки на аналогичной фотопластинке, обработанной вместе с исходной, впечатывается т.н. фотометрическая шкалка (или клин), содержащая несколько площадок (ступеней), равномерно засвеченных изменяющимся от ступени к ступени потоками. Определив, как меняется плотность почернения шкалки в зависимости от падающего потока , мы сможем перевести почернения на фотопластинке в относительные интенсивности. Изображение калибровочной шкалки находится в файле klin.bdf. Потоки, падавшие на ступени шкалки в относительных единицах (в логарифмах интенсивности) известны:

     

      2.0     1.83      1.68      1.54      1.36      1.18      0.99      0.80      2.0

 

Цифры расположены в порядке убывания интенсивности (почернения). Вычтем плотность вуали по всему изображению. Для этого воспользуйтесь командами load/image и fit/flatsky с ключом cursor, и отметьте площадки, не попадающие на ступени шкалки. Можно отмечать произвольное число площадок. По отмеченным областям будет построена поверхность и вычтена из начального изображения. В результате получается исправленная за фон плотность пропускания .

3.            С помощью команд load/image и statiscics/image (c ключом cursor) можно определить среднюю плотность почернения в каждой ступени клина, выделяя нужную область курсором. Эти данные следует занести в таблицу, для чего ее нужно сначала создать (командой create/table), затем проименовать колонки таблицы (с помощью create/column) и внести в нее данные (команда edit/table). Для дальнейшей обработки переведем данные из плотностей пропускания  в беккеровские плотности (W), связанные с разностью плотностей простым соотношением:

 

 

Эта промежуточная операция выполняется командой compute/table. Она нужна для того, чтобы дальнейшая полиномиальная аппроксимация сохраняла физический смысл величины: необходимо, чтобы  при .

4.            Следующий шаг — построение зависимости  от W и аппроксимация ее полиномом 4‑ой или 5‑ой степени (полученный полином должен быть монотонно возрастающим !). Для поиска коэффициентов регрессии воспользуйтесь командой regresssion/polynomial, затем сразу save/regression и compute/regression, чтобы записать в таблицу значение полинома с найденными коэффициентами в заданных точках (см. пример в Приложении). Вывести на экран график полученной зависимости (plot/table). Сюда же можно вывести исходные точки командой overplot/table. С помощью команды copy/graphisc скопируйте файл‑таблицу reference.tbl и отредактируйте ее так, чтобы можно было посмотреть на графике полученную полиномиальную функцию с достаточно малым шагом (используйте команды create/column, name/column, compute/column, overplot/table). Получив и проверив аппроксимирующий полином, мы можем перевести изображения спектров галактики и звезды А0 в интенсивности. Для этого необходимо вывести их на экран, (команды clear/display и load/image), вычесть вуаль (fit/flatsky, см. выше), перевести отсчеты в беккеровские плотности по формуле (1), а затем, используя характеристичкескую кривую, в интенсивности (переводить именно в интенсивности, а не в логарифмы интенсивностей!).

Примечание: Величина  должна быть положительной в каждом пикселе, т.к. мы имеем дело с логарифмами, поэтому рекомендуется воспользоваться в команде comp/image операцией max(10**(D–Df),10e-6).

Примечание: Зависимость характеристической кривой от длины волны не учитывается. Чтобы окончательно подготовить спектры звезды и галактики к дальнейшей обработке, вычтем из них фон (load/image, fit/flatsky. Не перепутайте фон с непрерывным спектром!).

5.            Построить кривую спектральной реакции аппаратуры. Для этого вывести изображение спектра звезды А0 (файл a0.bdf) на экран, определить номер строки, проходящей через середину спектра (get/cursor) и, используя  команды create/graphics, plot/row, set/graphics и get/gcursor, определить положения максимумов бальмеровских длин волн в пикселях и мировых координатах (см. Табл.1 и Приложение).

 

Таблица 1. Длины волн линий бальмеровской серии.

Линия	Hβ	Hγ	Hδ	Hε	Hζ	Hη	Hθ
λ, Å	4861	4340	4102	3970	3889	3835	3798

 

Построить зависимость между мировыми координатами на изображении (вдоль дисперсии) и длинами волн (plot/table и regress/polynomial). Оценить дисперсию (она будет несколько различаться для синей и красной частей спектра). Заметим, что в задаче мы исследуем относительно близкую галактику, поэтому полученную зависимость можно использовать для определения длин волн в спектре ее ядра напрямую. В общем же случае обычно учитывается красное смещение спектральных линий объекта.

6.            Используя полученную зависимость, а также изображение разреза спектра звезды A0 в графическом окне, определить значения интенсивности непрерывного  спектра (игнорируя линии поглощения) в диапазоне 3800 – 5000 Å Измерения проводить через каждые 100 Å. Теоретические величины  интенсивностей  непрерывного спектра звезды А0 приводятся в табл.2:

 

Таблица 2. Теоретические величины интенсивностей непрерывного спектра звезды класса А0.

λ, Å	I	λ, Å	I	λ, Å	I	λ, Å	I
3800	5.6	4200	7.25	4500	5.7	4800	5.11
3900	7.97	4300	6.8	4600	5.7	4900	4.85
4000	8.0	4400	5.70	4700	5.37	5000	4.6
4100	7.80

 

Интенсивности даны в произвольных единицах, поскольку распределение энергии в спектре здесь используется только для учета изменения чувствительности приемника (фотоэмульсии) в зависимости от длины волны. Сравнение результатов измерений с данными, полученными при абсолютных измерениях стандартных звезд А0, позволяет построить функцию спектральной реакции аппаратуры. Она получается путем деления истинных интенсивностей непрерывного спектра звезды А0 (см. табл.2) на регистрируемые. Полученную функцию следует вывести на экран (командами create/table, create/column, edit/table, plot/table), аппроксимировать полиномом (команда regression/polynomial), вывести полиномиальную функцию на экран (так же, как и в п.4), а затем сохранить в postscript–файл (set/graphics и copy/graphics). Это нужно для дальнейшей обработки.

7.            Вывести изображение спектра галактики на экран, определить номер строки, соответствующей разрезу через центр галактики (get/cursor). Провести командой filter/median медианную фильтрацию спектра галактики с шагом  пикселя. Выделить центральную строку спектра (команда extract/line) и скопировать ее в таблицу (команда copy/it). Добавить в полученную таблицу колонку из таблицы reference.tbl (команда copy/tt). Отождествить любую из линий спектра (напр. [OIII] 5007Å или H_β) и найти нуль–пункт координатной сетки. Изменить координатную сетку таким образом, чтобы одна мировая координата соответствовала одному ангстрему (команда compute/table). Исправить данные в таблице за дисперсию и функцию спектральной реакции (команда compute/table). Создать новую таблицу с единственной колонкой — интенсивностями спектра и преобразовать таблицу в новое одномерное изображение (команды copy/table, delete/column, copy/ti). Используя команды plot/table, set/graphics и get/gcursor, вывести в графическом окне спектр ядра галактики из таблицы, а командой load/image — одномерный спектр. Произвести отождествление ярких эмиссионных линий в спектре ядра: [OIII] 5007Å, [OIII] 4959Å, H_β, HeII 4686Å, [OIII] 4363Å, H_γ, [SII] 4069Å, [SII] 4076Å, H_δ, [NeIII]+H_ε, [NeIII]+H_ζ.

Примечание: Самая яркая линия [OIII] 5007Å передержана, поэтому ее вид (в интенсивностях) может иметь М–образную форму — в центральной части линии будет провал, что связано с переполнением отсчетов в тех пикселях ПЗС–матрицы, куда подал очень интенсивный поток.

Примечание: Дальнейшая работа ведется одновременно с преобразованным из таблицы спектром (в окне дисплея) и с таблицей (в графическом окне). Последняя нужна для отождествления (связи длин волн, выраженными в ангстремах, с пикселями и мировыми координатами). Изображение необходимо для нахождения светимостей спектральных линий. Спектр ядра галактики необходимо распечатать (используя команду copy/graphics), нанести отождествленные линии и представить в качестве одного из результатов. Отождествив эмиссионные линии, определить их относительные интенсивности. Для этого с помощью команды get/gcursor определите границы линий. Используя полученные данные, измерьте интенсивности линий I_ls с помощью команды statistics/image (можно определять границы линий непосредственно по изображению и сразу получать их интенсивности при помощи команды statistics/image с ключом cursor; это — более быстрый, но менее точный метод). Измерьте интенсивность непрерывного спектра слева и справа от линии. Среднее значение последних двух величин дает нам значение непрерывного спектра I_s. Для получения интенсивности отдельно в линияхI_l вычитаем I_s из I_ls.

Примечание: Линия [OIII] 4363Å блендируется с H_γ, поэтому измерения могут привести к завышению ее интенсивности. Для потока в передержанной линии [OIII] 5007Å следует принять его теоретическое значение: I([OIII] 5007) = 3·I([OIII] 4959), затем найти отношение

 

 

Аналогичное отношение получить для линий серы:

 

Примечание: Интенсивность красных линий серы известна из независимых наблюдений и составляет 1.5·I(H_β ).

 

 

8.            Полученные интенсивности линий и, в частности, их комбинации f₁([OIII]) и f₂([SII]), можно использовать для получения плотности и температуры газа в ядре галактики. Воспользуйтесь для этого номограммами, расчитанными в Крымской обсерватории [2] (см. рис.1). Выбирая из семейства кривых f([OIII]) и f([SII]) те, которые соответствуют наблюденным значением, найти по пересечению кривых величины n_e и T_e. В принципе, n_e и T_e можно найти и по интенсивностям отдельных линий или дуплетов по отношению к H_β:

 

где  — отношение концентрации соответствующих ионов к концентрации протонов (оно обычно считается равным «нормальному» содержанию данного элемента по отношению к Н), а  — функции, вид которых дан в [2].

9.            Измерить эквивалентную ширину линии H_β. Для этого с помощью команды statistics/image с ключом cursor найти истинную интенсивность этой линии I_l (как и в п.7), и разделить ее на высоту непрерывного спектра под линией. Последняя находится делением интенсивности непрерывного спектра I_s (полученной с тем же окном, что и линия) на размер окна вдоль дисперсии (размер в пикселях выводится на дисплей при выполнении statistics/image). Переведите, использую таблицу, эквивалентную ширину в ангстремы.

10.        Определить доплеровскую скорость движения газа. По линиям [OIII] 4959Å и H_β найти полуширину линии на половине интенсивности (сначала в пикселях, затем перевести в ангстремы). По формуле Доплера найти скорость движения газа (в километрах в секунду) вдоль луча зрения. Взять среднее значение по двум линиям.

11.        Определить параметры газа в ядре сейфертовской галактики. Для этого найти сначала расстояние D до галактики по скорости ее удаления (см. любой каталог галактик). Значение постоянной Хаббла принимаем равным 75 (км/с)/Мпк. Известно, что звезда величины  дает поток 3.4·10^–9 erg cm^–2Å ^–1 c^–1 на длине волны 5000Å. Звёздная величина ядра галактики NGC 1068 равна примерно  (для другой галактики она указывается преподавателем). Найти по формуле Погсона, какой поток дает объект  в спектральном интервале, равном измеренной выше эквивалентной ширине H_β. Умножая это значение на 4π D², где D — расстояние до галактики, вычислить светимость ядра L в линии H_β. Далее определить коэффициент излучения оптически тонкого газа по известным значениям его плотности и температуры:

 

Параметр A зависит от температуры:

 

T/10000K	1	2	3	4
A	1.24	0.56	0.42	0.20

 

Тогда объем газа, излучающего в линии H_β, можно вычислить по формуле:

 

 

Умножая объём на известное значение электронной плотности  и массу атома водорода , получаем оценку массы светящегося газа:

 

Кинетическая энергия газа может быть оценена по доплеровской скорости u, определенной по полуширине эмиссионных линий

 

Коэффициент 3 учитывает, что мы измеряем лишь лучевую, а не пространственную скорость (здесь распределение скоростей принимается изотропным). Порядок величины размера r излучающей области равен . Можно оценить также характерное кинематическое время движения газа как t=r/u.

Полученные данные необходимо представить в таблице, куда войдут следующие параметры:

·         электронная плотность газа n_e;

·         электронная температура Т_e;

·         светимость в линии H_β;

·         объем газа, излучающий в линии H_β;

·         масса газа, излучающего в линии H_β;

·         кинетическая энергия газа, излучающего в линии H_β;

·         характерное кинематическое время t.

Скачать эту задачу в формате Ms.Word