Эволюция дисковых галактик: исследование иерархии структур
МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТИМЕНИ М.В. ЛОМОНОСОВА
Государственный астрономический институт
имени П.К. Штернберга
На правах рукописи
УДК 524.5-7;
52-17
Хоперсков Сергей Александрович
ЭВОЛЮЦИЯ ДИСКОВЫХ ГАЛАКТИК:
ИССЛЕДОВАНИЕ ИЕРАРХИИ
СТРУКТУР
01.03.02 – Астрофизика и звездная астрономия
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
Москва – 2013
Работа выполнена в отделе физики звездных систем Федерального государственного бюджетного учреждении науки Института астро номии Российской академии наук.
Научный руководитель: член-корреспондент РАН, доктор физико-математических наук, профессор, Шустов Борис Михайлович директор Института астрономии РАН
Официальные оппоненты: доктор физико-математических наук, профессор, Расторгуев Алексей Сергеевич зав. отдела Изучения Галактики и переменных звезд ГАИШ МГУ доктор физико-математических наук, Моисеенко Сергей Григорьевич ведущий научный сотрудник отдела физики космической плазмы ИКИ РАН
Ведущая организация: Южный федеральный университет
Защита состоится «4» апреля 2013 года в 14:00 на заседании дис сертационного совета Д501.001.86 при Государственном астрономи ческом институте им. П.К.Штернберга (ГАИШ МГУ), расположен ном по адресу: Университетский проспект, дом 13, 119991, Москва
С диссертацией можно ознакомиться в Научной библиотеке Мос ковского государственного университета имени М.В. Ломоносова (119991, г. Москва, Ломоносовский проспект, д.27, Фундаменталь ная библиотека).
Автореферат разослан «4» марта 2013 года
Ученый секретарь диссертационного совета доктор физико-математических наук С.О. Алексеев
Общая характеристика работы
Актуальность работы Галактики представляют собой сложные, многокомпонентные системы, характеризующиеся разнообразием физических процес сов, протекающих в них. На ранних этапах формирование структу ры галактик происходит в основном под действием внешних факто ров — в вириализованных гало темной материи (ТМ) остывая и на капливая угловой момент скапливается первичный барионный газ, с началом звездообразования происходит формирование протога лактик и первых галактик. Наблюдения таких объектов на больших красных смещениях 5 10 стали блестящим подтверждением современных теоретических представлений, основанных в основном на CDM модели (см., например, [1]). В более близкой Вселенной 0.5 взаимодействие с окружением проявляется в виде столкно вений галактик сопоставимых размеров и выпадении карликовых галактик-спутников на более массивные галактики [2]. Во всех этих процессах важнейшая роль в управлении галактической эволюцией отводится темному веществу. Влияние эволюции темной материи на изменение структуры галактик проявляется, например, в регулиро вании толщины галактических дисков [3] и образовании различных морфологических особенностей в них — баров и/или спиралей [4, 5].
Существенными факторами секулярной эволюции галактических дисков являются формирование и разрушение бара и спиральной структуры [6], связанная с этим радиальная миграция звезд [7] и пе рераспределение химического состава межзвездной среды [8]. Раз личные коллективные процессы в галактических дисках в свою оче редь влияют на строение галактик. Перераспределение газа и его концентрация в отдельных областях под действием как глобаль ной галактической динамики, так и неоднородностей магнитного поля, потока излучения и других факторов межзвездной среды, определяют эффективность звездообразования. При определенных условиях динамические процессы в дисковых компонентах способ ны влиять на структуру темной материи в пределах оптического радиуса галактики [9].
Принципиальная невозможность проведения натурного экс перимента в физике галактик и наблюдений большинства астрофи зических процессов в галактиках в реальном времени неизбежно приводит к повышению роли результатов численных эксперимен тов. Данные наблюдений, полученные в рамках многочисленных об зоров, таких как SDSS (Sloan Digital Sky Survey), GALEX (GALaxy Evolution EXplorer), 2MASS (2 Micron All Sky Survey), настолько обширны и детальны, что имеется определенное отставание в раз витии численных моделей с точки зрения их разрешения и учета значительного числа сложных физических факторов, необходимых для адекватного описания свойств реальных объектов. Прогресс в теоретическом изучении наблюдаемых структур и эволюции га лактик и других астрофизических систем немалой степени связан с развитием компьютерной техники и широким распространением суперкомпьютеров. Показательным является трехкратное относи тельное увеличение числа астрономических публикаций, основан ных на численных газодинамических и -body расчетах, в период с 1995 по 2012 год (данные SAO/NASA Astrophysics Data System).
Благодаря ресурсам вычислительных кластеров не просто сокраща ется время решения конкретной динамической задачи, а становится возможным учет новых факторов и физических процессов.
Наличие темного вещества в галактиках и характер его про странственного распределения восстанавливаются по косвенным при знакам. Самые различные кинематические и морфологические осо бенности видимого вещества позволяют получать определенные огра ничения на те или иные параметры темного вещества в галактиках.
Постоянство кривой вращения на больших расстояниях от центра галактики и устойчивость галактических дисков к гравитационным возмущениям позволяют оценивать массу гало темной материи [10].
Кинематика и морфология приливных галактических структур [11] и полярных колец также дают ограничения на форму гало [12].
Несмотря на это однозначной картины, описывающей глобальные свойства ТМ внутри оптического радиуса спиральных галактик, еще нет. Ограниченность наших знаний о свойствах ТМ делает важ ной задачей является поиск новых проявлений темного вещества в галактиках через особенности в эволюции барионного вещества.
Спиральная структура является важным фактором, опреде ляющим эволюцию галактик. Понимание природы этого феноме на позволяет с одной стороны оценить вклад этих структур в на блюдаемую дисперсию скоростей звезд в современную эпоху [13].
С другой стороны дает возможность понять на какой эволюцион ной стадии находится сама галактика в целом [14]. Несмотря на более чем полувековую историю исследования феномена галакти ческого спирального узора, до сих пор не предложено универсаль ной теории, объясняющей все основные наблюдаемые проявления этих образований. Проблема определения времени жизни спираль ного узора в галактиках является одной из ключевых [15]. Время существования спиралей очевидным образом оценивается в числен ных экспериментах, однако, прямое сравнение этого параметра с данными наблюдений не представляется возможным. Анализ кос венных наблюдаемых характеристик спирального узора, таких как угол закрутки, скорость вращения и амплитуда волны, тем не ме нее позволяет делать ограничения на теоретические модели. Наибо лее перспективным подходом является воспроизведение структуры конкретных галактик, в том числе и нашей Галактики, в численных экспериментах, основанных на согласовании результатов моделиро вания с фотометрическими и кинематическими данными наблюде ний.
В основе представлений о структуре нашей Галактики и ее спиральном узоре лежит анализ кинематики различных галактиче ских компонент и определение расстояний между ними. Данные по кинематике ассоциаций молодых ОВ-звезд позволяют определять пространственные масштабы в галактическом диске [16]. Одна из оценок расстояния от Солнца до центра Галактики, по собствен ным движениям цефеид равна 7.1 кпк [17], при этом разброс в определении этой величины другими методами лежит в диапазоне от 7 до 10 кпк. Измерения тригонометрических паралаксов галак тических мазерных источников являются хорошими источниками данных при построении кривой вращения [18]. Морфология нашей Галактики неоднократно пересматривалась и по настоящее время является предметом дискуссий. Различные предположения о коли честве спиральных рукавов и их геометрии [19], проблема наличия кольца в центре [20] и другие свойства Галактики нуждаются в на блюдательных подтверждениях. Стоит ожидать, что наши знания о пространственном распределении вещества и кинематике спираль ного узора нашей Галактики в ближайшее время будут расширены благодаря работе космического телескопа Gaia (ESO), запуск кото рого намечен на 2013 год.
Цель работы — детальное исследование взаимосвязанных процессов в дисках спиральных галактик на основе многомерных численных экспериментов с применением технологий параллель ных вычислений. Достижение поставленной цели предусматривает решение следующих основных задач:
1. Разработка и реализация универсальных численных алгорит мов моделирования динамики газовых, звездных и звездно газовых галактических структур с учетом физико-химических особенностей межзвездной среды.
2. Исследование газодинамических, гравитационных и тепловых неустойчивостей галактической ударной волны в спиральном галактическом потенциале.
3. Исследование и обоснование различных механизмов формиро вания галактических спиральных структур.
4. Изучение эволюции дисковых галактик с учетом взаимодей ствия барионного вещества с эволюционирующей темной ма терией.
Научная новизна работы Все результаты, выносимые на защиту являются новыми.
Созданы новые универсальные комплексы программ для по строения динамических моделей газовых, звездных и звездно-газо вых галактических дисков с учетом самогравитации газа, процессов радиационного нагрева и охлаждения, химической кинетики моле кулярного водорода.
С использованием трехмерных газодинамических расчетах впервые был определен механизм, ответственный за образование по лигональных структур (ПС) в галактиках. Впервые проведено ста тистическое сравнение параметров полигональных структур (ПС), полученных по результатам численного моделирования, с данны ми наблюдений для выборки из 200 галактик. Показано, что необ ходимым требованием для формирования ПС является медленное вращение галактического спирального узора.
Построена новая модель образования и динамики гигантских молекулярных облаков (ГМО) в нашей Галактике, позволившая провести сравнительный количественный анализ между результа тами теоретических расчетов и наблюдаемыми свойствами облач ной структуры. Получены доводы, указывающие на отсутствие так называемого «Молекулярного кольца» в диске нашей Галактики на расстоянии 3 5 кпк от центра.
Проведено обобщение теории формирования спиральной струк туры, основанной на нелинейной стадии развития гравитационно неустойчивых глобальных мод в диске. Вместо политропных мо делей газовых дисков рассмотрен неадиабатический газ, а также исследованы модели звездных дисков в бесстолкновительном при ближении ( -body) и самосогласованные звездно-газовые модели.
Впервые детально исследована динамика трехмерных газово го и звездного дисков в гравитационном потенциале триаксиально го темного гало. Спиральная структура в таких моделях существу ет на временах порядка возраста дисковых галактик и ее форми рование не приводит к динамическому разогреву звездного диска.
Впервые показан нестационарный характер вращения спиральной структуры в гравитационном поле неосесимметричного гало.
В рамках новой самосогласованной модели образования дис ковой галактики в течении 10 млрд лет было показано влияние эво люции барионного вещества на форму гало темной материи, при водящее к сглаживанию профиля плотности гало в центральной области галактики.
Практическая и научная значимость работы Научная ценность диссертации состоит в создании пакетов программ для моделирования многокомпонентной динамики спи ральных галактик. Результаты, изложенные в диссертации, могут иметь применение для постановки наблюдательных задач по иссле дованию морфологии и кинематики галактик. В частности, сформу лирована задача определения параметров неосесимметричного тем ного гало по вариациям угловой скорости вращения спирального узора на основе спектральных наблюдений.
Критерии формирования полигональных структур и их гео метрические параметры позволят получить оценки кинематики спи рального узора, а также лучше понять физику галактических удар ных волн (ГУВ). Появляется дополнительная возможность получе ния новых ограничений на сценарии развития как газодинамиче ских, так и гравитационных неустойчивостей в окрестности ГУВ.
Исследование пространственного распределения и кинемати ки ГМО в моделях и наблюдениях позволит выявить ограничения на форму кривой вращения и крупномасштабную морфологию Га лактики.
Разработанная методика для моделирования отдельных ре альных галактик с заданными физическими параметрами, в рамках звездно-газовых моделей, позволяет дополнить список галактик с модальной природой галактической спиральной структуры.
Проблема центрального каспа в космологических CDM мо делях по всей видимости связана с недостаточным пространствен ным разрешением и ограниченностью набора физических процес сов в моделях. Проведенные расчеты позволяют объяснить сглажи вание профиля плотности гало темной материи без привлечения экзотических сценариев (например, специфического спектра возму щений в ранней Вселенной, наличия смешанных ароматов у частиц темной материи и др.).
В разные годы работа поддерживалась грантами Российского фонда фундаментальных исследований (07-02-01204-а, 12-02-00685-а, 11-02-12247-офи-м-2011, 09-02-97021-р_поволжье_а, 10-02-00231), про граммами Отделения общей физики и астрономии РАН, фондом некоммерческих программ «Династия» и федеральной целевой про граммой «Научные и научно-педагогические кадры инновационной России».
Основные положения, выносимые на защиту:
1. Универсальный комплекс программ для моделирования само согласованной динамики бесстолкновительной системы -тел, описывающей звездные подсистемы и темное вещество, и га зодинамических процессов в галактиках и межзвездной сре де с учетом самогравитации, внешних гравитационных полей, химической трансформации вещества, радиационных процес сов нагрева и охлаждения газа, звездообразования и взрывов сверхновых звезд.
2. Многомерные газодинамические расчеты формирования по лигональных галактических структур, позволившие впервые доказать газодинамический механизм образования спрямлен ных участков спиральных галактических рукавов. Было пока зано, что длинноволновые возмущения фронта галактической ударной волны при развитии гофрировочной неустойчивости образуют полигональные структуры. Ключевым фактором, влияющим на спрямление фронта галактической ударной вол ны, является медленное вращения спирального узора.
3. Новая модель формирования облачной структуры Млечного Пути с учетом термодинамических свойств газа, химической кинетики молекулярного водорода, самогравитации и реали стичного галактического потенциала. В рамках динамической модели впервые дано объяснение наблюдаемым особенностям диаграммы «галактическая долгота — лучевая скорость» без учета кольцевой структуры на расстоянии 3–5 кпк от центра нашей Галактики.
4. Механизм формирования крупномасштабной спиральной струк туры в галактике NGC 5247, основанный на развитии гравита ционно неустойчивых глобальных мод в диске. Показана вос производимость в численных моделях наблюдаемого спираль ного узора в пределах неопределенности данных наблюдений о кривой вращения, радиальном профиле звездной плотности и дисперсии скоростей звезд. Расчеты в рамках самосогласо ванных звездно-газовых моделей позволяют сделать вывод о значительном увеличении времени жизни спирального узора с учетом холодной газовой компоненты.
5. Механизм генерации долгоживущей глобальной спиральной структуры в звездно-газовых галактических дисках за счет взаимодействия с неосесимметричным темным гало. Вывод о нестационарном вращении спиральной волны плотности в по ле неосесимметричного массивного гало.
6. Результаты расчетов разрушения центрального каспа плотно сти темной материи за счет динамического взаимодействия с барионной компонентой дисковой галактики в процессе ее формирования и эволюции.
Основные публикации по теме диссертации Статьи в рецензируемых изданиях:
1. Khoperskov S.A., Vasiliev E.O., Sobolev A.M., Khoperskov A.V.
The simulation of molecular clouds formation in the Milky Way // Monthly Notices of the Royal Astronomical Society.
2013, v. 428, № 3, p. 2311-2320.
2. Khoperskov S.A., Khoperskov A.V., Khrykin I.S., Korchagin V.I., Casetti-Dinescu D.I., Girard T., van Altena W., Maitra D. Global gravitationally-organized spiral waves and the structure of NGC // Monthly Notices of the Royal Astronomical Society.
2012, v. 427, № 3, p. 1983-1993.
3. Khoperskov A.V., Khoperskov S.A., Zasov A.V., Bizyaev D.V., Khrapov S.S. Interaction between collisionless galactic discs and nonaxissymmetric dark matter haloes // Monthly Notices of the Royal Astronomical Society. 2013, arXiv:1302. 4. Хоперсков С.А., Хоперсков А.В., Еремин М.А., Бутенко М.А.
Полигональные структуры в газовом диске: численные экс перименты // Письма 2011, т.37, в Астроном. Журнал.
№ 8, с. 614-627.
5. Khoperskov S.A., Eremin M.A., Khoperskov A.V. Polygonal Structures in the Gaseous Disks // Astronomical and Astro physical Transactions. 2012, v. 27, № 2, p. 245-250.
6. Хоперсков А.В., Еремин М.А., Хоперсков С.А., Бутенко М.А., Морозов А.Г. Динамика газового диска в неосесимметричном темном гало // Астрономический журнал. 2012, т. 89, № 1, с. 19-31.
7. Хоперсков C.А., Шустов Б.М., Хоперсков А.В. Взаимодей ствие каспа темного вещества с барионной составляющей в // 2012, Астрономический журнал.
дисковых галактиках т. 89, № 9, с. 736-744.
8. Бутенко М.А., Еремин М.А., Корчагин В.И., Морозов А.Г., Хоперсков С.А. Определение собственных мод для гравита ционной неустойчивости в газовом диске // Вестник Вол гоградского государственного университета. Серия 1:
2009, № 12, с. 70-72.
Математика. Физика.
9. Еремин М.А., Хоперсков А.В., Хоперсков С.А. Конечно-объем ная схема интегрирования уравнений гидродинамики // Из вестия Волгоградского государственного техническо го университета. Сер.: Актуальные проблемы управ ления, вычислительной техники и информатики в тех 2010, т. 6, № 8, c. 24-27.
нических системах.
10. Засов А.В., Хоперсков А.В., Тюрина Н.В., Еремин М.А., Хопер сков С.А. Использование суперкомпьютеров для моделирова ния внегалактических объектов // Суперкомпьютерные технологии в науке, образовании и промышленности.
Выпуск 2 / Под редакцией: академика В.А. Садовничего, ака демика Г.И. Савина, чл.-корр. РАН Вл.В. Воеводина. — М.:
Издательство Московского университета. 2010, c. 93- 11. Хоперсков А.В., Васильев Е.О., Хоперсков С.А., Соболев А.М., Еремин М.А. Модель образования молекулярных облаков в на шей Галактике. Роль темного гало. // Вестник Волгоград ского государственного университета. Серия 1: Мате 2011, т. 14, № 1, c. 93-98.
матика. Физика.
12. Хоперсков С.А., Хоперсков А.В., Засов А.В., Бутенко М.А.
Параллельный алгоритм для моделирования динамики газа // Вестник в сильно неоднородных гравитационных полях Уфимского государственного авиационного техниче ского университета. Серия: Управление, вычислитель 2012, т. 16, № 3, c. 108-114.
ная техника и информатика.
13. Бутенко М.А., Хоперсков С.А., Хоперсков А.В. Численное мо делирование внешних газовых спиралей в галактиках // Вест ник Волгоградского государственного университета.
2012, т. 15, № 1, с.49-56.
Серия 1: Математика. Физика.
В сборниках трудов конференций:
14. Хоперсков А.В., Бутенко М.А., Хоперсков С.А. Сверхбыстрые звезды и свойства темного галактического гало // Сборник трудов IV Всероссийского научного семинара «Физи 2008, Элиста, Изд-во КалмГУ. 2008, ка Солнца и звезд».
с. 106-113.
15. Бутенко М.А., Еремин М.А., Хоперсков C.A., Корчагин В.И., Хоперсков А.В. Глобальные спиральные моды, порождаемые гравитационной неустойчивостью в галактическом диске // Сборник трудов V Всероссийского научного семинара 2009, Элиста, Изд-во КалмГУ.
«Физика Солнца и Звезд».
2010, с. 37-39.
16. Хоперсков С.А., Еремин М.А., Хоперсков А.В. Шпуры в дис ках спиральных галактик // Труды конференции «Хими ческая и динамическая эволюция галактик» (28-30 сен тября, 2009, Ростов-на-Дону, ЮФУ) под ред. Бочкарева, Ще кинова, Ростов-на-Дону, Изд-во РСЭИ, 2010, с. 60-68.
17. Еремин М.А., Хоперсков С.А., Хоперсков А.В. Моделирование полигональных структур // Труды конференции «Хими ческая и динамическая эволюция галактик» (28-30 сен тября, 2009, Ростов-на-Дону, ЮФУ) под ред. Бочкарева, Ще кинова. Ростов-на-Дону, Изд-во РСЭИ, 2010, с. 51-59.
18. Хоперсков C.А., Еремин М.А., Хоперсков А.В. Параллельный код для моделирования динамики сверхзвукового газа с само // Сборник статей участников Всероссий гравитацией ского конкурса научных работ студентов и аспиран тов «Телематика’2010: телекоммуникации, веб-техно — СПб: СПбГУ ИТМО, 2010.
логии, суперкомпьютинг».
— c. 159-163.
19. Хоперсков С.А., Еремин М.А., Хоперсков А.В., Засов А.В., Тюрина Н.В., Применение высокопроизводительных вычисле ний для моделирования гидродинамических течений в сильно // Труды междуна неоднородных гравитационных полях Челябинск, родной научной конференции ПАВТ2011.
Издательский центр ЮУрГУ. 2011, c. 334-342.
Апробация работы. Основные результаты диссертации до кладывались на международных конференциях:
1. Nonstationary Phenomena and Instabilities in Astrophysics (NPIA), Volgograd, Russia, 8-12 September (2009).
2. Dynamics and Evolution of Disc Galaxies, Pushchino, Russia, May-4 June (2010).
3. Параллельные вычислительные технологии (ПАВТ2011), Москва, 28 марта – 1 апреля (2011).
4. Galaxies origin, dynamics, structure & astrophysical disks, Sochi, Russia, 14-18 May (2012).
5. European week of Astronomy and Space Science, Rome, Italy, 1- July (2012).
6. Galactic scale star formation, Heidelberg, Germany, 30 July – August (2012).
Всероссийских конференциях:
1. Химическая и динамическая эволюция галактик, Ростов-на Дону, 28-30 сентября (2009).
2. Галактические и аккреционные диски, Нижний Архыз, 21- сентября (2009).
3. Телематика-2010: телекоммуникации, веб-технологии, супер компьютинг, Санкт-Петербург, 21-24 июня (2010).
4. Всероссийская астрономическая конференция (ВАК-2010) «От эпохи Галилея до наших дней», Нижний Архыз, 13-18 сентяб ря (2010).
5. Физика космоса, Коуровская АО (2008, 2010, 2011).
6. Астрофизика высоких энергий сегодня и завтра, Москва (2010, 2011, 2012).
7. Актуальные проблемы внегалактической астрономии, Пущи но, 17-19 апреля (2012).
8. Наблюдательные проявления эволюции звезд, Нижний Архыз, 15-19 октября (2012).
Выступления на конференциях молодых ученых Института астро номии РАН (2011, 2012) и астрофизических семинарах:
1. Семинар ГАИШ МГУ по звездной астрономии (13.02.2013).
2. Семинар отдела Внегалактической астрономии ГАИШ МГУ (29.04.2010, 30.03.2011).
3. Объединенный семинар Института астрономии РАН (22.04.2010, 29.11.2012).
Работа по реализации параллельных технологий и исполь зованию суперкомпьютеров, описанная в Главе 1, получила Спе циальный приз от суперкомпьютерной программы «СКИФ-ГРИД»
на конкурсе Intel-РОСНАНО (Хоперсков С.А., Засов А.В., Тюрина Н.В. «Применение высокопроизводительных вычислений при моде лировании мелкомасштабных структур газовых подсистем спираль ных галактик», Москва, 2010 г.).
Личный вклад автора в совместных работах По теме диссертации опубликовано 19 работ. Из них 8 опуб ликовано в рецензируемых научных журналах из перечня ВАК.
Содержание диссертации и основные положения, выносимые на за щиту, отражают персональный вклад автора в опубликованные ра боты. Подготовка к публикации полученных результатов проводи лась совместно с соавторами, причем вклад диссертанта был опре деляющим. Все представленные в диссертации результаты получе ны лично автором. Соискатель в равной степени с другими соавто рами участвовал в постановке задач, имея определяющую роль на этапах разработки методов, их тестирования, проведения расчетов, получения и представления результатов и выводов.
Разработаны численные газодинамические алгоритмы, опи санные в диссертации. Адаптирован алгоритм TreeCode для рас четов газодинамических и гибридных -body / газодинамических моделей. Реализованы параллельные OpenMP и MPI алгоритмы.
Подавляющая часть численных экспериментов была проведена дис сертантом с использованием суперкомпьютеров НИВЦ МГУ («Че бышев» и «Ломоносов»).
Структура и объем диссертации Диссертация состоит из введения, четырех глав, заключения и библиографии. Общий объем диссертации 223 страницы, включая 69 рисунков, 4 таблицы и 4 приложения. Библиография включает 336 наименований на 38 страницах.
Основное содержание работы
Во Введении обсуждается актуальность работы, цели и за дачи исследования, научная новизна, научная и практическая цен ность полученных результатов. Также формулируются основные результаты и положения, выносимые на защиту, и приводится спи сок работ, в которых опубликованы основные результаты диссерта ции.
Первая глава посвящена описанию базовых идей числен ных алгоритмов, лежащих в основе численных экспериментов. В параграфе 1.1 обсуждается методика построения численной газо динамической схемы типа Годунова, описывается приближенный метод решения задачи Римана для уравнений газовой динамики.
В параграфе 1.2 приведены уравнения Гамильтона, описывающие динамику системы гравитационно взаимодействующих -тел. Кро ме того рассмотрены два наиболее распространенных подхода к решению уравнения Пуассона: с помощью быстрого дискретного преобразования Фурье и метод, основанный на иерархическом ал горитме TreeCode. Обсуждаются особенности метода объединения газодинамических и -body моделей. В параграфе 1.3 дано описа ние использованных в работе методов распараллеливания газодина мического алгоритма, рассмотренного выше в параграфе 1.1. При ведены результаты тестирования эффективности численного кода для различных постановок задач на суперкомпьютере «Чебышев»
(НИВЦ МГУ). Показано, что наиболее эффективное использова ние вычислительных мощностей достигается в случае трехмерной декомпозиции расчетной области и максимальном учете факторов в МЗС. При этом эффективность достигает 70% при расчетах на 200 потоках.
Во второй главе исследованы динамические процессы в га зовой компоненте дисковых галактик, связанные с проявлением спиральной волны плотности и галактической ударной волны (ГУВ).
В параграфе 2.1 приведен краткий обзор современных представле ний о структуре и наблюдаемых свойствах ГУВ. В п. 2.2.1 собраны многочисленные наблюдаемые свидетельства спрямленных сегмен тов спиральных рукавов внешних галактик — так называемых по лигональных структур (ПС) [21]. К настоящему моменту известно более 200 спиральных галактик, обладающих такой особенностью.
Причем полигональные структуры получены как в галактиках ти па «grand-design» (например, NGC 5194), так и во взаимодействую щих галактиках (NGC 1512), а также в галактиках с флоккулент ной спиральной структурой (M 101). В п. 2.2.2 построена газоди намическая модель формирования полигональных структур. Пока зано, что этот феномен связан с неустойчивостью ударной волны, которая на мелких масштабах проявляется в виде гофрировочной неустойчивости. Определены параметры способствующие развитию ПС. В п. 2.2.3 проведено сравнение наблюдаемых параметров ПС, с полученными в численных экспериментах. В параграфе 2.3 изу чены наиболее актуальные, с точки зрения динамической эволю ции Галактики, физические характеристики молекулярных обла ков. В п. 2.4.1 построена многокомпонентная химико-динамическая модель образования ГМО с учетом внешнего поля Галактики, само гравитации, тепловой структуры газа и химической кинетики мо лекулярного водорода. В п. 2.4.2 обсуждаются закономерности рас пределения физических параметров ГМО по результатам числен ных экспериментов: распределение облаков по массам, соотношение «масса-размер» и дисперсия скоростей газа. В п. 2.4.3 приводится обоснование идеи, об отсутствии материальной кольцевой структу ры в диске Галактики на расстоянии 35 кпк от центра. Анализ син тетических карт «лучевая скорость — галактическая долгота» для молекулярного газа обнаруживает особенность, часто интерпрети руемую как «Молекулярное кольцо», однако кольцевой структуры в диске не наблюдается. В тоже время нейтральный в построенных в работе моделях газ более равномерно заполняет галактических диск.
Третья глава посвящена изучению феномена и механиз мов образования галактической спиральной структуры. В парагра фе 3.1 описаны современные представления о процессе формирова ния и устойчивости спиральной структуры в дисковых галактиках за счет внутренних факторов без учета приливных сил. Обсужда ются доводы в пользу модальной природы и длительного, в тече ние миллиардов лет, времени существования спирального узора со слабо меняющейся морфологией. В п. 3.1.1 рассмотрена спираль ная галактика NGC 5247 и на ее примере изучено формирование спиральных волн за счет гравитационной неустойчивости глобаль ных спиральных мод. В п. 3.1.2 сравниваются результаты линей ного анализа устойчивости гравитирующего диска в гидродинами ческом приближении с эволюцией нелинейных газодинамических расчетов. В п. 3.1.3 обсуждаются результаты численных -body и -body/газодинамических экспериментов, подтверждающих воз можность продолжительного времени жизни спирального узора, формирующегося на нелинейной стадии развития неустойчивой мо ды. Пункт 3.1.4 посвящен изучению влияния газовой компоненты на время существования крупномасштабного регулярного спираль ного узора типа «grand design». Численные эксперименты указы вают на увеличение времени жизни спиралей в несколько раз в самосогласованном звездно-газовом диске. Дальнейшая часть гла вы (параграф 3.2) посвящена исследованию динамики дисковых подсистем в гравитационном потенциале массивного неосесиммет ричного темного гало. Основной упор сделан на изучение механиз мов генерации спиральных узоров и физических свойств формиру ющихся спиральных волн. В п. 3.2.1 проанализированы многочис ленные факты, подтверждающие наличие массивного темного гало в галактиках. В п. 3.2.2 подробно рассмотрены различные наблю дательные и теоретические указания на триаксиальный характер распределения темной материи внутри сферы 1 3 оптических ра диуса галактики. Собраны оценки параметров трехосности, полу ченные различными методами на основе космологических моделей и наблюдаемых свойств галактических структур (полярные коль ца, изгибы газовых и пылевых дисков, их утолщения, приливные потоки, карликовые спутники, гипербыстрые звезды и др.).
Рассмотрена гипотеза о том, что неосесимметричность (три аксиальность) темного гало может быть ответственна за формиро вание спиральных структур в галактических дисках по крайней ме ре у части галактик. Пункт 3.2.3 содержит в себя описание моделей звездных и газовых дисков в поле трехосного гало и базовых пара метров, определяющих ход эволюции систем. В п. 3.2.4 и п. 3.2. приведены результаты расчетов динамики газовых и звездных си стем соответственно. Главным результатом является вывод о неиз бежном формировании долгоживущего спирального узора в галак тических звездных и газовых дисках. Отличительной особенностью возникающего спирального узора является его нестационарный ха рактер, что проявляется, в частности, в нестационарном и нетвер дотельном вращении волны плотности. Параметры вращения спи ралей испытывают квазипериодические изменения. В п. 3.2.6 об суждается механизм нестационарности спирального узора в неосе симметричном гравитационном поле темного гало. Описываются перспективы использования этой особенности для измерения пара метров формы гало реальных галактик. В п. 3.2.7 изучены эффек ты, возникающие при взаимодействии спирального узора с динами ческим темным гало. В частности, рассмотрен процесс изменения параметров формы гало в результате динамической эволюции. По казано, что перераспределение энергии между гало и спиральными волнами приводит к сферизации гало в центральной области с со хранением основных особенностей спиралей структуры в диске.
В четвертой главе в рамках численных экспериментов изу чены процессы образования дисковой галактики и эффекты взаи модействия видимого вещества (газа и звезд) со скрытой компо нентой галактики в результате динамической эволюции галакти ческого диска. В параграфе 4.1 обсуждаются ключевые проблемы теории формирования галактик в рамках современной CDM кон цепции. Особое внимание здесь уделено проблеме центрального кас па плотности темной материи. В параграфе 4.2 построена модель с учетом основных физических факторов, необходимых для изу чения долговременной эволюции дисковой галактики. В парагра фе 4.3 рассмотрены различные сценарии эволюции темного веще ства под действием барионной компоненты. Детально рассмотрен процесс трансформации профиля темной материи при учете про цессов звездообразования и звездной эволюции в центре галактиче ского диска. В качестве главного результата исследования выделим полученный эффект сглаживания профиля плотности темного ве щества в модели с учетом звездообразования и звездной эволюции, по сравнению с расчетами, основанными только на моделях -тел.
Таким образом для дисковых галактик именно звездообразование и эволюция звезд приводит к эффективному оттоку вещества (в том числе темного) из центральной области галактического диска, и, как следствие, к сглаживанию каспа плотности ТМ.
В Заключении перечислены основные результаты диссерта ции.
В Приложении А приведены основные тесты для числен ной газодинамической схемы типа TVD MUSCL, реализованной в рамках диссертации и использованной для решения поставленных задач. Помимо стандартной задачи о распаде произвольного скачка давления, большое внимание уделено тестированию адекватности описания процесса распространения газа в вакуум. Продемонстри рована универсальность реализованной схемы и сохранение пози тивности при решении задачи о взрыве в вакууме. Кроме того рас смотрена эволюция начального состояния, приводящего к развитию неустойчивости Кельвина-Гельмгольца. В Приложении Б приве дена модель гравитационного потенциала Галактики, использован ная при химико-динамических расчетах формирования облачной структуры. Описан вклад темного гало, балджа и звездного диска в кривую вращения, а также определены параметры спирального четырех-рукавного потенциала. В Приложении В перечислены процессы нагрева и охлаждения в модели газовой компоненты Га лактики, необходимые для адекватного моделирования двухфазной среды. В Приложении Г описана методика построения равновес ного начального состояния для звездного и звездно-газового галак тических дисков на основе совместного решения уравнений Джинса и Пуассона.
Цитированная литература 1. Bouwens R. J., Illingworth G. D., Franx M., Ford H. z 7-10 Galax ies in the HUDF and GOODS Fields: UV Luminosity Functions // Astrophys. J. 2008. Vol. 686. Pp. 230–250.
2. Klypin A., Kravtsov A. V., Valenzuela O., Prada F. Where Are the Missing Galactic Satellites? // Astrophys. J. 1999. Vol. 522.
Pp. 82–92.
3. Khoperskov A., Bizyaev D., Tiurina N., Butenko M. Numerical modelling of the vertical structure and dark halo parameters in disc galaxies // Astronomische Nachrichten. 2010. Vol. 331.
Pp. 731–745.
4. Athanassoula E. Bar-Halo Interaction and Bar Growth // Astro phys. J. Let. 2002. Vol. 569. Pp. L83–L86.
5. Bekki K., Freeman K. C. Extended H I Spiral Structure and the Figure Rotation of Triaxial Dark Halos // Astrophys. J. Let. 2002.
Vol. 574. Pp. L21–L24.
6. Minchev I., Famaey B., Quillen A. C. et al. Evolution of galac tic discs: multiple patterns, radial migration, and disc outskirts // Astron. and Astrophys. 2012. Vol. 548. Pp. 126–149.
7. Rokar R., Debattista V. P., Quinn T. R. et al. Riding the Spi s ral Waves: Implications of Stellar Migration for the Properties of Galactic Disks // Astrophys. J. Let. 2008. Vol. 684. Pp. L79–L82.
8. de Avillez M. A., Mac Low M.-M. Mixing Timescales in a Super nova-driven Interstellar Medium // Astrophys. J. 2002. Vol. 581.
Pp. 1047–1060.
9. Abadi M. G., Navarro J. F., Fardal M. et al. Galaxy-induced trans formation of dark matter haloes // Monthly Notices Roy. Astron.
Soc. 2010. Vol. 407. Pp. 435–446.
10. Athanassoula E., Bosma A., Papaioannou S. Halo parameters of spiral galaxies // Astron. and Astrophys. 1987. Vol. 179. Pp. 23–40.
11. Bournaud F., Duc P.-A., Masset F. The large extent of dark matter haloes probed by the formation of tidal dwarf galaxies // Astron.
and Astrophys. 2003. Vol. 411. Pp. L469–L472.
12. Snaith O. N., Gibson B. K., Brook C. B. et al. The halo shape and evolution of polar disc galaxies // Monthly Notices Roy. Astron.
Soc. 2012. Vol. 425. Pp. 1967–1979.
13. Minchev I., Quillen A. C. Radial heating of a galactic disc by mul tiple spiral density waves // Monthly Notices Roy. Astron. Soc.
2006. Vol. 368. Pp. 623–636.
14. Sil’chenko O. K., Proshina I. S., Shulga A. P., Koposov S. E. Ages and abundances in large-scale stellar discs of nearby S0 galaxies // Monthly Notices Roy. Astron. Soc. 2012. Vol. 427. Pp. 790–805.
15. Sellwood J. A. The lifetimes of spiral patterns in disc galaxies // Monthly Notices Roy. Astron. Soc. 2011. Vol. 410. Pp. 1637–1646.
16. Dambis A. K., Mel’Nik A. M., Rastorguev A. S. Trigonometric Parallaxes and a Kinematically Adjusted Distance Scale for OB Associations // Astronomy Letters. 2001. Vol. 27. Pp. 58–64.
17. Dambis A. K., Mel’Nik A. M., Rastorguev A. S. Rotation curve of the system of classical Cepheids and the distance to the galactic center // Astronomy Letters. 1995. Vol. 21. Pp. 291–307.
18. Bobylev V. V., Bajkova A. T. Galactic parameters from masers with trigonometric parallaxes // Monthly Notices Roy. Astron. Soc.
2010. Vol. 408. Pp. 1788–1795.
19. Valle J. P. The Spiral Arms and Interarm Separation of the Milky e Way: An Updated Statistical Study // Astron. J. 2005. Vol. 130.
Pp. 569–575.
20. Dobbs C. L., Burkert A. The myth of the molecular ring // Monthly Notices Roy. Astron. Soc. 2012. Vol. 421. Pp. 2940–2946.
21. Chernin A. D., Zasov A. V., Arkhipova V. P., Kravtsova A. S.
Vorontsov-Velyaminov Rows: Straight Segments in the Spiral Arms of Galaxies // Astronomy Letters. 2000. Vol. 26. Pp. 285–296.