Научная статья на тему 'Внутренняя кинематика сейфертовской галактики Mkn 938'

Внутренняя кинематика сейфертовской галактики Mkn 938 Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

CC BY
33
8
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МЕТОДЫ: НАБЛЮДАТЕЛЬНЫЕ / МЕТОДИКА: ФОТОМЕТРИЯ / ГАЛАКТИКИ: АКТИВНЫЕ / METHODS: OBSERVATIONAL / TECHNIQUES: PHOTOMETRIC / GALAXIES: ACTIVE

Аннотация научной статьи по наукам о Земле и смежным экологическим наукам, автор научной работы — Афанасьев В.Л., Моисеев А.В., Смирнова А.А.

В данной работе мы представляем результаты детального исследования центральной области сейфертовской галактики Mkn 938. Наблюдательный материал был получен на 6-м телескопе САО РАН с панорамным спектрографом MPFS и сканирующим интерферометром Фабри-Перо. Mkn 938 интересна тем, что является результатом слияния двух богатых газом галактик, причем мы наблюдаем уже финальную стадию слияния, сопровождаемую мощнейшей вспышкой звездообразования и активностью ядра. Анализ кинематики газа и звезд показал, что в околоядерной области Mkn 938 присутствует истечение газа со скоростями от -370 км с-1 до -480 км с-1, а на больших масштабах нами впервые картирован высокоскоростной галактический ветер в линии поглощения Na D.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INTERNAL KINEMATICS OF THE SEYFERT GALAXY MKN 938

We report the results of a detailed study of the central part of the Seyfert galaxy Mkn 938. Observational data were obtained with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences using integral-field spectrograph MPFS and a scanning Fabry-Perot interferometer. Mkn 938 is interesting for being a result of a merger of two gas-rich galaxies, and we observe the final stage of this merger accompanied with an extremely powerful burst of star formation and nuclear activity. Our analysis of the kinematics of gas and stars revealed the presence of gas outflow in the circumnuclear region Mkn 938 with velocities ranging from -370 to -480 kms-1, and allowed us for the first time to map the high-velocity galactic wind in Na D absorption line on large spatial scale.

Текст научной работы на тему «Внутренняя кинематика сейфертовской галактики Mkn 938»

АСТРОФИЗИЧЕСКИЙ БЮЛЛЕТЕНЬ, 2020, том 75, № 1, с. 13-22

УДК 524.74-32/75520.82

ВНУТРЕННЯЯ КИНЕМАТИКА СЕЙФЕРТОВСКОЙ ГАЛАКТИКИ

Mkn 938

© 2020 В.Л.Афанасьев1*, А.В.Моисеев1, А.А.Смирнова1

1Специальная астрофизическая обсерватория РАН, Нижний Архыз, 369167 Россия Поступила в редакцию 24 сентября 2019 года; после доработки 14 ноября 2019 года; принята к публикации 14 ноября 2019 года

В данной работе мы представляем результаты детального исследования центральной области сей-фертовской галактики Mkn938. Наблюдательный материал был получен на 6-м телескопе САО РАН с панорамным спектрографом MPFS и сканирующим интерферометром Фабри-Перо. Mkn938 интересна тем, что является результатом слияния двух богатых газом галактик, причем мы наблюдаем уже финальную стадию слияния, сопровождаемую мощнейшей вспышкой звездообразования и активностью ядра. Анализ кинематики газа и звезд показал, что в околоядерной области Mkn938 присутствует истечение газа со скоростями от —370 км с-1 до —480 км с-1, а на больших масштабах нами впервые картирован высокоскоростной галактический ветер в линии поглощения № D.

Ключевые слова: методы: наблюдательные — методика: фотометрия — галактики: активные

1. ВВЕДЕНИЕ

Исследования сливающихся галактик, богатых газом, позволяют понять роль диссипативных процессов наряду с гравитационными в эволюции галактик и их систем (Schweizer and Seitzer 2007). За последние десятилетия был достигнут прогресс в понимании общей роли слияний в образовании галактик раннего типа, и сейчас принято считать, что многие из них являются продуктом слияний богатых газом галактик. Однако не совсем ясно, как богатая газом система в конечном счете избавляется от газа и как на это влияет вспышка звездообразования в центре и наличие активного ядра. Яркие области звездообразовния и истечения формируют различные компоненты на луче зрения, не разрешаемые в картинной плоскости. Применение методов двухмерной (называемой так же 3D или интегрально-полевой) спектроскопии позволяет детально изучать подобные объекты. В настоящей работе представлены результаты двухмерной спектроскопии центральной части галактики Mkn 938. Основные параметры Mkn 938 приводятся в таблице 1.

Взаимодействующая галактика Mkn 938 (=NGC34 = VV850), имеющая пекулярную морфологию и активное ядро, была классифицирована нами как Sy2 (Afanasev et al. 1980), что было в дальнейшем подтверждено в работах Dahari

E-mail: vafan@sao.ru

(1985), Veron-Cetty and Veron (1986). Напротив, в Osterbrock and Dahari (1983) она классифицирована как эмиссионная галактика, в которой доминируют области звездообразования, а в Mulchaey et al. (1996) на изображениях в узких фильтрах обнаружено, что в Mkn 938 эмиссия в линии [OIII] А5007 слаба по сравнению с большинством известных сейфертовских галактик. Кроме того, там же обнаружено, что яркая эмиссия Ha видна по всей галактике. Это предполагает, что ионизация большей части газа не связана с какой-либо сейфертовской активностью. Морфологически галактика, по-видимому, переживает слияние, о чем свидетельствует наличие приливных хвостов. Это подтверждается инфракрасными изображениями Mkn 938 на 8.8 и 12.5 мкм (Miles et al. 1996), на которых в ядерной области галактики был обнаружен двойной источник с разделением примерно 1'!2. Предварительные результаты интегрально-полевой спектроскопии центральной части этой галактики, полученные на 6-м телескопе САО РАН (Rafanelli et al. 2000), показали, что в центре есть как признаки звездообразования, так и истечения газа. В частности, в линии Na D были обнаружены два динамических центра на расстоянии 1", которые совпадают с двойным источником, обнаруженным в ИК. Обнаружение в этой галактике на VLA протяженной структуры в H I (Fernandez et al. 2010) и околоядерного диска в CO на ALMA (Xu et al. 2014) свидетельствует о том, что она является богатой газом взаимо-

АФАНАСЬЕВ и др. Таблица 1. Основные параметры Mkn938

Тип активности Sy2

Морфологический тип (NED) Pec

Системная скорость (Rothberg and Joseph 2006) 5881 kmc"1

Принятое расстояние 85.2 Мпк

Соответствующая шкала 395 пк/"

Абсолютная звездная величина Му (Schweizer and Seitzer2007) -21.42

Масса нейтрального водорода (Kandalyan 2003) 5.3 х 109 М0

Максимальная скорость вращения

(наши измерения для г = 38°) 225 ±6

действующей системой. Глубокие изображения в оптическом диапазоне (Schweizer and Seitzer 2007), также указывают на то, что Mkn 938 является вероятным продуктом слияния галактик с соотношением масс 1: 3. Высокая ИК-светимость объекта предполагает, что звездообразование является доминирующим в центре с незначительным вкладом от AGN (Gonealves et al. 1999), что подтверждается детальным исследованием инфракрасного спектра Mkn 938 (Esquej et al. 2012), где болометрический вклад AGN в общую ИК-светимость оценивается в 2%. Кроме того, ИК-изображения показывают, что звездообразование происходит внутри 0.5—2 кпк, где имеются следы слияния. Однако эти авторы считают, что наличие AGN необходимо для объяснения светимости в жестком рентгене.

2. НАБЛЮДЕНИЯ

Mkn 938 наблюдалась нами в оптическом диа-

пазоне на 6-м телескопе Специальной астрофи-

зической обсерватории РАН в рамках программы

спектрофотометрии галактик Сейферта. Гелио-

центрическая скорость центра галактики, определенная по абсорбционным спектрам (Rothberg

and Joseph 2006), составляет 5881 ± 2 км с-1, что дает расстояние до галактики 85.2 Мпк для H0 = 73 км с-1 Мпк-1. Это соответствует масштабу в картинной плоскости 395 пк/".

Журнал наблюдений Mkn 938 на 6-м телескопе приведен в таблице 2, где указаны: дата наблюдений, суммарная экспозиция, размер построенного куба данных, качество изображений, используемый прибор, дисперсер или фильтр, спектральный диапазон и спектральное разрешение. Описание используемых приборов (интегрально-полевого волоконного спектрографа MPFS и редуктора светосилы с интерферометром Фабри-Перо — ИФП)

можно найти в Afanasiev et al. (2001). Обработка полученных обоими приборами данных выполнялась по стандартной методике, описанной в наших предыдущих работах Smirnova and Moiseev(2010), Smirnova et al. (2006).

2.1. Интерферометр Фабри-Перо

При наблюдениях со сканирующим ИФП было последовательно получено 12 интерферограмм, равномерно заполняющих свободный спектральный диапазон (АЛ w 29 A для данного порядка интерференции). Спектральная область вокруг крас-носмещённой линии Ha выделялась узкополосным фильтром. Для экономии времени считывания и повышения отношения сигнал/шум ПЗС детектор TK 1025 1K.x 1K считывался в режиме аппаратного биннинга 2x2. Поле зрения составляло 5i8 при масштабе 0'!68 на пиксель.

Результатом обработки являлся куб данных, где каждый пиксель в поле зрения содержит 12-канальный спектр. Астрометрическая привязка куба данных выполнялась средствами программы Aladin (Boch and Fernique 2014). С помощью гаусс-аппроксимации профиля эмиссионной линии были построены распределения яркости в линии Ha и поле лучевых скоростей ионизованного газа, показанные на рис. 1. По техническим причинам нам пришлось ограничиться небольшим числом интеферограмм, дающих дискретизацию спектров, близкую к ширине аппаратного контура прибора, поэтому карт распределения дисперсии скоростей в эмиссионной линии построить не удалось.

Эмиссионная линия детектируется внутри яркого звездного диска галактики (r < 15"—20"), а также в эмиссионном «островке», находящемся на проекционном расстоянии 30"—55" к северо-востоку от ядра. Этот «островок» хорошо заметен и

Таблица 2. Log of observations

Date oí Total exp, Size of Seeing, Device Disperser Spectral Spectral

observation sec. data cube arcsec or Filter coverage, A resolution, A

1999/12/04 2400 512x512x12 2 focal reducer FP260 6675-7005 2.5

2009/10/18 9600 16x16x3000 1.5 MPFS 1200/17.7 4790-7740 3

30

20

и 0)

S 10

о

«о <1

0

-10

-20

На velocity

40 30 20 10 0 -10 -20 Да, arcsec

Рис. 1. (a) Изображение Mkn938, полученное на HST с камерой ACS в фильтре F435W Квадратом отмечено поле, наблюдавшееся с MPFS. (b) и (c) — результаты наблюдений со сканирующим ИФП: соответственно распределение яркости в линии Ha и поле лучевых скоростей в линии. Шкала скоростей в кмс-1. Крестиком показан центр изофот континуума вблизи Ha.

на широкополосных изображениях с космического телескопа Хаббла (HST), полученных с камерой ACS и опубликованных в Kim et al. (2013), а также на более ранних снимках из Schweizer and Seitzer (2007). Снимок в синем фильтре F438W показан на рис. 1 внизу. Эта эмиссионная деталь, скорее всего, является областью интенсивного звездообразования, связанного с протяженным приливным рукавом. Возможно, это приливная галактика размером около 6 кпк.

2.2. Панорамный спектрограф МРРБ

В случае МРРБ одновременно регистрируются спектры от 256 пространственных элементов, формируемых при помощи прямоугольного линзового растра 16" х 16", который помещается в фокальную плоскость. Угловой размер каждой линзы в наших наблюдениях был 1"х1". Этот линзовый растр формирует массив микрозрачков, которые реформируются световодами в псевдощель на входе спектрографа. Пробы спектра фона неба берутся другими световодами, расположенным на расстоянии 4' от линзового растра. В качестве детектора использовалась ПЗС матрица Е2У42-90 формата 2кх4.5к. Расположение растра МРРБ на изображении галактики показано на рис. 1.

После первичной обработки результаты были представлены в виде кубов данных с полем зрения 16х 16 элементов, каждый из которых содержал спектр, включающий яркие эмиссионные линии, типичные для сейфертовских галактик — Нв, [О III] А4959,5007, На, [N11] А6548,6584 и [Б II] А6717,6731, и линии поглощения Нв, MgI и N8 0.

Для измерения потоков и лучевых скоростей эмиссионные линии аппроксимировались гаусси-анами. Вдали от центра галактики их профили хорошо описываются однокомпонентной моделью, но в центральной области (г < 4'', см. далее раздел 3.2) наблюдается второй широкий компонент с амплитудой до 1/4 от основного. Далее поля скоростей и карты потоков приводятся для основного (узкого) компонента. Для учета вклада контрастных абсорбций от бальмеровских линий звездного населения абсорбция в Нв также описывалась гауссианой. В области На+^ II] скорость и ширина абсорбционной детали принимались такими же, как в области Нв, а амплитуда выбиралась такой, чтобы сохранить ту же эквивалентную ширину линий.

Поле скоростей звездного населения строилось методом кросс-корреляции (Мо^ееу 2001) с использованием шаблонных спектров звезд в спектральном диапазоне 5015—5400 А, содержащем

LOG(Flux, 10"15, erg/s/cm2) 0.5 1.0 1.5 2.0 2.5 3.0

8 6 4

о

° п от 2 о

о

о о

ш о < -2

-4

-6

-8 -6 -4 -2 0 2 4 6 ДРА.агсэес

Рис. 2. Поле скоростей Мкп938 в линии На, наложенное на изображение центральной части галактики в континууме в полосе V.

контрастные абсорбции Mg I и Ре I. Этим же методом анализировались профили лучевых скоростей в абсорбционной линии дублета N8 О, поскольку здесь лучевые скорости сильно отличаются от наблюдаемых как в звездном населении, так и ионизованном газе (см. далее раздел 3.2). В качестве шаблонного спектра N8 О использовались эмиссионные линии ночного неба, регистрируемые в поле зрения МРРБ.

3. РЕЗУЛЬТАТЫ

3.1. Кинематика газа и звезд

Согласно данным МРРБ, распределение лучевых скоростей ионизованного газа в области г < 5''—7'' в целом согласуется с картиной, ожидаемой для вращающегося плоского галактического диска (рис. 2). Это подтверждается и наблюдениями с ИФП (рис. 1) в линии На. В этой области поле скоростей хорошо совпадает с радиоданными по Ш (Регпаndez е1 а1. 2010). Но за пределами области, наблюдавшейся с МРРБ, поле скоростей ионизованного газа демонстрирует значительные отличия от кругового вращения, особенно к западу от ядра. В то же время, лучевые скорости внешней эмиссионной области к северо-востоку от ядра (вероятной приливной галактики) соответствуют системной скорости центра Мкп 938.

Анализ полей лучевых скоростей методом «наклонных колец» (ИН^-г^б) позволяет определить в узком кольце на заданном расстоянии от ядра г следующие параметры: скорость вращения У^, позиционный угол кинематической большой

250

200

'от

Е 150

100

>

50

> оо

20

сл

си 0

X)

<

-20

-40

А

О

О ф

л

♦ +

о

А

На IFP О На MPFS Л stars •

I | | | I | | | I | | | I

I | | | I

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Н-

2 4 6 8 10 12 14 г, arcsec

Рис. 3. Результаты анализа полей скоростей газа и звезд: кривые вращения (верхний рисунок) и позиционный угол кинематической большой оси (внизу). Пунктиром показано изменение ориентации большой оси изофот в полосе K (Rothberg and Joseph 2004), штрихпунктиром — ориентация большой оси внутреннего диска согласно Schweizer and Seitzer (2007).

структуре, хорошо заметной на снимках HST (см. рис. 1). На больших расстояниях от ядра РАк\п в ионизованном газе начинает уже сильно отличаться от фотометрической большой оси, что вызвано либо ростом величины некруговых (радиальных) движений газа, либо изгибом газового диска. Но последнее менее вероятно, так как кривая вращения на г = 3"—10" остается практически плоской. Параметры вращения звезд на г > 4" заметно отличаются от газа: кинематическая ось отклоняется в другую сторону от фотометрической, а скорость вращения для принятой величины г резко падает. Отметим, что и во внутренних областях Угоъ для звезд почти в два раза меньше, чем для ионизованного газа. Такое различие сложно объяснить только эффектом асимметричного дрейфа (большей дисперсией скоростей) в звездном диске. Скорее всего, в поле скоростей звезд мы видим минимум две динамические подсистемы (основную галактику и ядро сливающегося с ней спутника). Сложение их скоростей по лучу зрения усложняет наблюдаемую картину и не дает ее точно интерпретировать в рамках модели кругового вращения плоского диска. В то же время в диссипативной газовой системе мы наблюдаем уже установившееся квази-круговое движение газовых облаков.

оси PAkin, наклон к лучу зрения i и системную скорость Fsys. Подробнее о приложении этого метода к данным ИФП и MPFS рассказано в наших предыдущих работах (Smirnova and Moiseev 2010, Smirnova et al. 2018). Наклон принимался фиксированным i = 38°, согласно фотометрическим оценкам Rothberg and Joseph (2004) для изофот внутреннего диска в полосе K, свободной от пылевого поглощения. Fsys фиксировалась в соответствии со средней величиной по полю.

Сперва из соображений симметрии поля лучевых скоростей определялось положение кинематического центра вращения. В ионизованном газе (измерения в Ha с MPFS и ИФП) центр вращения совпадает с центром изофот континуума. Однако в поле скоростей звезд кинематический центр смещен от фотометрического примерно на 1" к северу.

На рис. 3 показано распределение кинематических параметров вдоль радиуса. Видно, что во внутренней области r < 4" (1.6 кпк) PAkin ~ —20°, что близко к имеющимся оценкам ориентации внутреннего диска галактики в оптическом (Schweizer and Seitzer 2007) и ближнем ИК (Rothberg and Joseph 2004) диапазонах: PAphot ~ —10°. Небольшое различие позиционных углов естественно объясняется радиальными движениями вещества в околоядерной спиральной

3.2. Истечение ионизованного газа

По данным MPFS профили эмиссионных линий как в самом ядре, так и в пределах г < 4" к югу от фотометрического центра демонстрируют заметную синюю асимметрию. Чтобы изучить ее с лучшим отношением сигнал/шум, мы просуммировали спектры в этой области, приведя их в покоящуюся систему отсчета в соответствии с полем скоростей в линии Ha. На рис. 2 показано, что интегральный спектр в каждой из эмиссионных линий (На и дублеты [ЫП] А6548,6583, [Э II]А6717,6731, [ОШ]А4959,5007) можно описать минимум двумя гауссовыми компонентами — узким, с шириной близкой к аппаратному контуру спектрографа, и в несколько раз более широким (ЕШИИ к 700-1000 кмс-1). При этом во всех запрещенных линиях широкий компонент заметно смещен в синюю сторону, относительные лучевые скорости составляют около —370 кмс-1, — —380 км с-1 и —480 км с-1 соответственно в линиях [Ы II], [Э II] и [О III].

В тоже время, широкий компонент в бальме-ровской линии На имеет почти нулевую скорость относительно узкой. К сожалению, сделать такой анализ в эмиссии Н@ не удалось из-за ее низкой яркости и блендирования звездной абсорбцией Нв. Однако наши данные позволяют построить карты распределения интегральной яркости линий [О III]

4940 4960 4980 5000 5020 Wavelength, Ä

6550 6600 6650

Wavelength, Á

6700

Рис. 4. Интегральный спектр области, в которой по данным МРРБ наблюдается асимметрия профилей спектральных линий. Слева — эмиссионные линии дублета [О III], справа — линии На, ^ II] и [Б II]. Показано разложение профилей на узкий и широкий гауссовские компоненты, показаны остатки от вычитания компонент. Для удобства отображения спектры смещены вдоль вертикальной оси.

-8 -6 -4 -2 0 2 4 6 ARA, arcsec

Рис. 5. Изофоты изображения в эмиссии Нв, наложенные на изображение в линии [O III] 5007.

и деблендированной эмиссии Hß, что показано на рис. 5. На рисунке видно заметное различие распределения яркости в линиях [O III] и Hß. Если максимум яркости в [O III] находится в с центре галактики, то в Hß он смещен на запад на 2".

Изофоты изображения в имеют форму дуги, центр которой смещен на юго-запад на 1.5—2". Аналогичное смещение наблюдается также в линии Ha по данным ИФП и MPFS. Средняя лучевая скорость, определенная по линии поглощения H^ равна 5890 км с-1, что совпадает с системной скоростью 5880 км с-1. Напротив, узкая эмиссия H^ показывает смещение на —200 км с-1.

Наблюдаемую картину можно интерпретировать как ветровое истечение из ядра со скоростями не менее 400—500 км с-1, что часто наблюдается в галактиках с мощным звездообразованием. При этом в галактическом ветре доминирует ударная ионизация (Heckman et al. 1990, Westmoquette et al. 2012), на что также указывает смещение узкой линии H@. Наличие в центре несмещенного широкого компонента в линии Ha указывает на слабую активность ядра и связано с областью образования широких линий (BLR).

Альтернативным объяснением наблюдаемой асимметрии линий может быть взаимодействие джета с окружающей его средой (jet-cloud interaction, см. работу Smirnova et al. (2007)). Однако детальные наблюдения с высоким пространственным разрешением в радиодиапазоне не обнаружили присутствия протяженных структур, которые можно было бы интерпретировать как джет (Fernandez et al. 2014).

3.3. Галактический ветер в протяженной линии поглощения Na D

В спектре центральной части галактики наблюдается широкая линия поглощения Na D, смещенная на 620 км с-1 в голубую область относительно системной скорости 5880 км с-1 (Schweizer

EW, А

5000 6000 Velocity, km/s

Isophotes of continuum + map EW(NaD)

5000 6000 Velocity, km/s

Рис.6. Карта распределения эквивалентных ширин линии № D с наложенными на нее изофотами распределения яркости континуума в диапазоне 5500—6300 А и примеры спектров и их анализа в различных частях изображения.

and Seitzer 2007). Такие линии, демонстрирующие наличие высокоскоростного истечения из центра, наблюдаются у многих ультраярких в ИК-диапазоне активных галактик с интенсивным звездообразованием (Rupke et al. 2005). По данным MPFS мы проанализировали изменение наблюдаемых профилей дублета натрия. рис. 6 иллюстрирует идею нашего анализа. На центральной панели показана карта распределения эквивалентных ширин линии поглощения Na D, которые вычислялись в интервале 5970—6020 A. На левой и правой панелях рисунка приведены профили линий Na D и результат кросс-корреляции с линиями Na D ночного неба. На оси абсцисс указана шкала скоростей. Рядом с профилями линий объекта показаны негативные профили линий ночного неба, смещенные по Z. Пунктирной вертикальной линией отмечена системная скорость галактики. Кривые кросс-корреляции, по-существу, являются профилями распределения лучевых скоростей в линии Na D в соответствующей точке изображения. Декомпозиция этих профилей показала, что их форма хорошо описывается двумя гауссианами приблизительно одинаковой ширины в пределах FWHM = 600-900 км с-1. На рис. 7 показаны карты лучевых скоростей каждой из компонент. Сильно смещенный в голубую область компонент мы будем далее называть высокоскоростным, а другой — низкоскоростным. На рисунке изофоты показывают распределение эквивалентных ширин NaD. Высокоскоростной компонент в центре показывает скорость —800 км с-1 относительно си-

стемной скорости и ширину порядка 700 км с 1, а низкоскоростной —200 км с-1 и 900 км с-1 соответственно.

4. ОБСУЖДЕНИЕ

В Mkn 938 мы наблюдаем финальную стадию слияния двух богатых газом галактик с доминирующим звездообразованием в центре. Наши наблюдения обнаружили многочисленные особенности кинематики газа и звезд в этой галактике. Во-первых, профили эмиссионных линий как в самом ядре, так и в пределах r < 4" к югу от фотометрического центра демонстрируют заметную синюю асимметрию. И если во всех запрещенных линиях широкий компонент заметно смещен в синюю сторону (относительные лучевые скорости составляют от —370 км с-1 до —480 км с-1 в разных линиях), то широкий компонент в бальмеровской линии Ha имеет почти нулевую скорость относительно системной. Наличие в центре несмещенного широкого компонента в линии Ha связанно с областью образования широких линий (BLR) вокруг активного ядра Mkn938.

Наблюдаемую картину можно интерпретировать как ветровое истечения из ядра со скоростями не менее 400-500 км с-1, что часто наблюдается в галактиках с мощным звездообразованием, при этом в галактическом ветре доминирует ударная ионизация (Heckman et al. 1990, Westmoquette et al. 2012), на что также указывает смещение узкой линии H@. Наличие в центре несмещенного широкого компонента в линии Ha указывает на слабую

Velocity, km/s 4600 4800 5000 5200 5400 5600 5800 6000

-6 -4 -2 0 2 4 6 ARA, arcsec

о <

-2 -

-4 -

-6 -4 -2 0 2 4 6 ARA, arcsec

Рис. 7. Карты лучевых скоростей высокоскоростного (вверху) и низкоскоростного (внизу) компонент линии NaD. Распределение эквивалентных ширин NaD показано изофотами.

активность ядра и связано с областью образования широких линий (BLR).

Асимметрия линий скорее всего вызвана истечением из активного ядра. Эта гипотеза подтверждается и другими авторами, оценивавшими скорость движения газа из центра Mkn938 по NaD от —620 км с-1 до —1050 км с-1 (Schweizer and Seitzer 2007). Подобного рода истечения могут наблюдаться и в галактиках с сильным

темпом звездообразования, но даже для темпа 70 ± 20 Mq год-1 такая скорость истечения слишком велика (Schweizer and Seitzer 2007), и объяснить ее можно только воздействием активного ядра. Например, наличием джета, пока не обнаруженного в радионаблюдениях.

Высокоскоростные истечения нейтрального газа в абсорбции Na D наблюдаются в галактических ветрах достаточно часто. Как правило, эти наблюдения выполняются методом длиннощелевой спектроскопии, так что пространственная структура и протяженность истечений остается не ясной. Поля лучевых скоростей истечений в Na D пока построены лишь для небольшого числа галактических ветров (смотри, к примеру, нашу работу про Mkn334 (Smirnova and Moiseev 20l0) или недавнее исследование ветра в NGC 5394 (Martín-Fernandez et al. 2016)). При этом пространственный размер истечения нейтрального газа в Mkn 938 значительно больше, высокоскоростной газ наблюдается на проекционном расстояниях до 3 кпк от центра. Построение подробной геометрической модели этого истечения мы планируем выполнить в отдельной работе с привлечением новых наблюдательных данных.

Наши данные позволяют построить карты распределения интегральной яркости линий [OIII] и деблендированной эмиссии Hß (см. рис. 5). И если максимум яркости в [OIII] совпадает с изображением ядра галактики в континууме, то в Hß он смещен на запад на 2". При этом эмиссия в линии Hß сдвинута от системной скорости на —200 км с-1. Наиболее вероятным объяснением наблюдаемого эффекта является истечение газа из ядра со скоростями, достигающими величины 400—500 км с-1. Истечение газа со скоростью приблизительно 400 км с-1 было также обнаружено на ALMA (см. работу Xu et al. 2014).

Наблюдения Mkn 938 в эмиссионной линии CO (6—5) на ALMA обнаружили динамическую структуру, которую авторы проинтерпретировали как околоядерный вращающийся диск (Xu et al. 2014). Однако, этот диск оказывается сильно смещенным по скорости (скорость центра диска приблизительно 5700 км с-1), что почти на 200 км с-1 меньше системной скорости. Данное смещение динамического центра можно объяснить тем, что на самом деле в линии CO (6—5) мы наблюдаем совсем не диск, а истечение со скоростями порядка 200 км с-1, и ядро галактики тогда соответствует скоростям около 5900 км с-1. Эта величина системной скорости хорошо согласуется с многочисленными наблюдениями других авторов (VSyS = 5931 ± 11 км с-1 — RC3/NED, 5881 ± 2 км с-1 — Rothberg and Joseph (2006) и т.д.).

5. ЗАКЛЮЧЕНИЕ

Выполнено детальное исследование активной галактики Mkn938 с применением методов панорамной спектроскопии. В ходе работы удалось изучить особенности структуры Mkn 938 на разных пространственных масштабах, а именно:

1. На больших расстояния от центра (r > 4") параметры вращения звезд заметно отличаются от газа. В поле скоростей звезд мы видим минимум две динамические подсистемы: основную галактику и ядро сливающегося с ней спутника, поэтому сложение их скоростей по лучу зрения усложняет наблюдаемую картину и не дает её точно интерпретировать в рамках модели кругового вращения плоского диска. В то же время в газовой системе мы наблюдаем уже установившееся квази-круговое движение газовых облаков, позволяющее построить реалистичную кривую вращения.

2. Нами впервые картирован высокоскоростной галактический ветер в линии поглощения NaD. Декомпозиция профилей этой линии показала, что их форма хорошо описывается двумя гаус-сианами приблизительно одинаковой ширины в пределах FWHM = 600-900 км с-1. Нами были построены поля скоростей для обоих компонент: высокоскоростной и низкоскоростной. Высокоскоростной компонент в центре показывает скорость —800 км с-1 относительно системной скорости и ширину порядка 700 км с-1, а низкоскоростной —200 км с-1 и 900 км с-1 соответственно.

3. В самой внутренней области галактики, в пределах r < 4" к югу от фотометрического центра, также наблюдается второй компонент в эмиссионных линиях ионизованного газа, смещенный в синюю сторону до —500 км с-1. Кроме того, центр яркости в бальмеровских линиях заметно (почти на 2", т.е. около 0.8 кпк) смещен относительно центра изофот [O III] и континуума. Все вместе это может интерпретироваться либо как ионизованная составляющая галактического ветра, картируемого в Na D, либо как результат воздействия джета из активного ядра на окружающую межзвездную среду. К сожалению, имеющиеся радионаблюдения Mkn 938 пока джета не обнаруживают.

БЛАГОДАРНОСТИ

Авторы благодарят анонимного рецензента за замечания, позволившие улучшить статью. В процессе работы над статьей использовалась база внегалактических данных ЫАБА/1РАС

(NED), управляемая Лабораторией реактивного движения Калифорнийского технологического института по контракту с NASA(США). В статье приводится снимок из базы данных HLA (Hubble Legacy Archive) оперирующей наблюдения с телескопа NASA/ESA Hubble Space Telescope по коллаборации между Space Telescope Science Institute (STScl/NASA), Space Telescope European Coordinating Facility (ST-ECF/ESA), Canadian Astronomy Data Centre (CADC/NRC/CSA).

ФИНАНСИРОВАНИЕ

Работа выполнена при финансовой поддержке Российского научного фонда (проект 17-12-01335 «Ионизованный газ в галактических дисках и за пределами оптического радиуса»). Наблюдения на телескопах САО РАН выполняются при поддержке Министерства науки и высшего образования Российской Федерации (включая соглашение № 05.619.21.0016, уникальный идентификатор проекта RFMEFI61919X0016).

КОНФЛИКТ ИНТЕРЕСОВ

Авторы заявляют об отсутствии конфликта интересов.

СПИСОК ЛИТЕРАТУРЫ

1.V. L. Afanasev, V. A. Lipovetskii, B. E. Markarian, and D. A. Stepanian, Astrofizika 16, 193 (1980).

2. V. L. Afanasiev, S. N. Dodonov, and A. V. Moiseev,

in L. P. Ossipkov and I. I. Nikiforov (eds.), Stellar Dynamics: from Classic to Modern, p. 103 (2001).

3. T. Boch and P. Fernique, ASP Conf. Ser. 485, 277

(2014).

4. O. Dahari, Astrophys. J. Suppl. 57, 643(1985).

5. P. Esquej, A. Alonso-Herrero, A. M. Perez-Garcia,

et al., Monthly Notices Royal Astron. Soc. 423 (1), 185 (2012).

6. X. Fernandez, A. O. Petric, F. Schweizer, and J. H. van

Gorkom, Astron. J. 147 (4), 74 (2014).

7. X. Ferna ndez, J. H. van Gorkom, F. Schweizer, and J. E.

Barnes, Astron. J. 140 (6), 1965(2010).

8. A. C. Goncalves, M. P. Veron-Cetty, and P. Veron,

Astron. and Astrophys. Suppl. 135,437(1999).

9. T. M. Heckman, L. Armus, and G. K. Miley,

Astrophys. J. Suppl. 74, 833 (1990).

10. R. A. Kandalyan, Astron. and Astrophys. 398, 493 (2003).

11. D. C. Kim, A. S. Evans, T. Vavilkin, et al., Astrophys. J. 768 (2), 102(2013).

12. P. Martin-Fernandez, J. Jimenez-Vicente, A. Zurita, et al., Monthly Notices Royal Astron. Soc. 461 (1), 6 (2016).

13. J. W. Miles, J. R. Houck, T. L. Hayward, and M. L. N. Ashby, Astrophys. J. 465, 191 (1996).

14. A. V. Moiseev, Bulletin of the Special Astrophysical Observatory 51, 11 (2001).

15. J. S. Mulchaey, A. S. Wilson, and Z. Tsvetanov, Astrophys. J. Suppl. 102, 309(1996).

16. D. E. Osterbrock and O. Dahari, Astrophys. J. 273, 478(1983).

17. P. Rafanelli, A. Rifatto, V. Afanasiev, et al., in W. van Breugel and J. Bland-Hawthorn (eds.), Imaging the Universe in Three Dimensions, Astronomical Society of the Pacific Conference Series, vol. 195, p. 232 (2000).

18. B. Rothbergand R. D. Joseph, Astron. J. 128 (5), 2098 (2004).

19. B. Rothberg and R. D. Joseph, Astron. J. 131 (1), 185 (2006).

20. D. S. Rupke, S. Veilleux, and D. B. Sanders, Astrophys. J. 632 (2), 751 (2005).

21. F. Schweizer and P. Seitzer, Astron. J. 133 (5), 2132 (2007).

22. A. Smirnova and A. Moiseev, Monthly Notices Royal Astron. Soc. 401 (1), 307(2010).

23. A. A. Smirnova, N. Gavrilovic, A. V. Moiseev, et al., Monthly Notices Royal Astron. Soc. 377 (2), 480 (2007).

24. A. A. Smirnova, A. V. Moiseev, and V. L. Afanasiev, Astronomy Letters 32 (8), 520 (2006).

25. A. A. Smirnova, A. V. Moiseev, and S. N. Dodonov, Monthly Notices Royal Astron. Soc. 481 (4), 4542 (2018).

26. M. P. Veron-Cetty and P. Veron, Astron. and Astrophys. Suppl. 65, 241 (1986).

27. M. S. Westmoquette, D. L. Clements, G. J. Bendo, and S. A. Khan, Monthly Notices Royal Astron. Soc. 424 (1), 416(2012).

28. C. K. Xu, C. Cao, N. Lu, et al., Astrophys. J. 787 (1), 48(2014).

Internal Kinematics of the Seyfert Galaxy Mkn 938

V. L. Afanasiev, A. V. Moiseev, and A. A. Smirnova

We report the results of a detailed study of the central part of the Seyfert galaxy Mkn 938. Observational data were obtained with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences using integral-field spectrograph MPFS and a scanning Fabry—Perot interferometer. Mkn 938 is interesting for being a result of a merger of two gas-rich galaxies, and we observe the final stage of this merger accompanied with an extremely powerful burst of star formation and nuclear activity. Our analysis of the kinematics of gas and stars revealed the presence of gas outflow in the circumnuclear region Mkn 938 with velocities ranging from —370 to —480 km s~1, and allowed us for the first time to map the high-velocity galactic wind in Na D absorption line on large spatial scale.

Keywords: methods: observational—techniques: photometric—galaxies: active

i Надоели баннеры? Вы всегда можете отключить рекламу.