Сейфертовские галактики представляют собой редко встречающийся тип галактик, обладающих активным ядром. Светимость их яркого звездообразного ядра иногда превышает 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_l