Моделирование излучения компактных рентгеновских источников

Прошло уже около 40 лет как было осознано, что падение вещества на компактный релятивистский объект (нейтронную звезду или черную дыру) приводит к выделению большого количества энергии. Запущенный в 1970 г. спутник Uhuru с рентгеновским детектором на борту открыл несколько сотен рентгеновских источников, часть из которых оказалась двойными, т. е. содержали кроме источника рентгеновского излучения еще и нормальную звезду. Так было окончательно установлено, что в процессе эволюции тесных двойных систем существуют фазы, во время которых более массивная звезда уже превратилась в компактный релятивистский объект, а вторая, менее массивная звезда, начала интенсивно терять массу либо в виде ветра, либо переполнив свою полость Роша. Часть теряемой массы перетекает (аккрецирует) на компактный объект.

Из-за орбитального движения аккрецируемое вещество обладает значительным угловым моментом и формирует вокруг компактного объекта геометрически тонкий аккреционный диск, вещество в котором вращается со скоростью, близкой к кеплеровской, и по туго закрученной спирали дрейфует к центру благодаря турбулентной вязкости. Теория таких аккреционных дисков была построена в работах Шакуры (1972) и Шакуры и Сюняева (1973, 1976). Поправки к теории за счет эффектов общей теории относительности впервые были рассмотрены в работе Новикова и Торна (1973).

Примерно в это же время «открылись глаза» на давно известный класс переменных звезд — катаклизмические, или взрывные переменные звезды. Было осознано, что это те же тесные двойные системы с перетеканием вещества, только компактным объектом является белый карлик. Поэтому максимум излучения катаклизмических переменных высокой светимости с оптически толстыми аккреционными дисками приходится на ультрафиолетовый диапазон спектра.

В настоящее время качественно и количественно (в смысле общей энергетики) процессы, происходящие при аккреции вещества в тесных двойных системах достаточно хорошо поняты. Наступил следующий этап исследований этих объектов, в течение которого научному сообществу предстоит научиться детально моделировать спектры излучения тесных двойных, систем и определять по наблюдаемым спектрам их фундаментальные параметры, подобно тому, как это делается для обычных звезд. Развитие техники рентгеновских наблюдений позволяет надеяться на получение в недалеком будущем рентгеновских спектров с хорошим разрешением и отношением сигнал/шум. Уже сейчас рентгеновский телескоп, установленный на орбитальной обсерватории Chandra способен получать спектры с разрешением АЕ/Е ~ 1000.

Вполне естественно, что для моделирования спектров излучения аккреционных потоков, по крайней мере, оптически толстых, необходимо использовать методы, применяемые при расчетах моделей атмосфер и спектров излучения обычных звезд (см. например Михалас 1982), модифицировав их применительно к рассматриваемым объектам.

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

1) Излучение при выделении гравитационной потенциальной энергии выпадающего вещества. В простейшем случае радиального падения гравитационная потенциальная энергия падающего вещества переходит в кинетическую, так что вещество перед падением на поверхность компактного объекта разгоняется вплоть до второй космической скорости. После удара о поверхность практически вся его кинетическая энергия переходит в тепло (внутреннюю энергию) и излучается.

В реальных двойных системах радиальное падение вещества осуществляется лишь при наличии сильного магнитного поля у компактного объекта, которое разрушает диск на радиусе Альвена (где сравнивается магнитное давление и динамическое давление вещества в диске). После этого радиуса вещество движется по магнитным силовым линиям и практически радиально выпадает на поверхность вблизи магнитных полюсов. Вблизи поверхности после фронта ударной волны формируется аккреционная колонка, в которой нагретое вещество, охлаждаясь, оседает на поверхность компактного объекта в дозвуковом режиме. Температура в этой колонке близка к вириальной лишь в случае магнитных белых карликов (поляры и промежуточные поляры), так как плазма в этом случае остается оптически тонкой.

В случае слабого магнитного поля компактного объекта аккреционный диск простирается вплоть до его поверхности, и энергия в диске выделяется и высвечивается постепенно, на всех радиусах аккреционного диска. Максимум энерговыделения находится вблизи внутренней границы диска. Эффективность энерговыделения в диске е довольно высока, е ~ 0.1 для аккреционных дисков вокруг нейтронных звезд и черных дыр, и е ~ Ю-4 для аккреционных дисков вокруг белых карликов.

При аккреции на нейтронную звезду или белый карлик со слабым магнитным полем аккреционный диск простирается вплоть до поверхности компактного объекта, его вещество движется вблизи поверхности с кеплеровской скоростью и должно потерять избыточную кинетическую энергию, прежде чем выпасть на компактный объект. Эта кинетическая энергия переходит в тепло и излучается благодаря трению вещества диска о компактный объект. При этом выделяется практически столько же энергии (эта величина зависит также от скорости вращения компактного объекта), как и в аккреционном диске. Такой переходный слой между аккреционным диском и компактной звездой носит название пограничного слоя. Такие пограничные слои, по всей видимости, дают значительный вклад в рентгеновское излучение маломассивных рентгеновских источников высокой светимости и в мягкое рентгеновское излучение некоторых катаклизмических переменных в высоком (ярком) состоянии (например, SS Лебедя).

Следовательно, гравитационная потенциальная энергия аккрецируемого вещества может выделяться в трех различных структурах — аккреционном диске, пограничном слое и аккреционной колонке.

2) Термоядерное горение вещества на поверхности компактного объекта.

В результате аккреции на поверхности компактного объекта скапливается вещество, образуя оболочку, богатую водородом. По мере накопления вещества в оболочке возникают условия (поверхностная плотность оболочки должна быть ~ 108 — 109 г/см2) для его термоядерного горения, что может привести к термоядерному взрыву скопившегося вещества. При достаточно высоком темпе аккреции такое горение может быть квазистационарным, как в слоевом источнике проэволюционировавших звезд.

Сравним эффективность энерговыделения при аккреции и при термоядерном горении, равную? ~ 0.007 при горении водорода, и е ~ 0.001 — 0.002 при горении гелия. Оказывается, что эффективность энерговыделения при дисковой аккреции на нейтронную звезду почти на 2 порядка выше, чем при термоядерном горении выпавшего вещества. В то же время, эффективность аккреции на белый карлик в 10−30 раз меньше, чем при термоядерном горении. Отсюда следует вывод, что при высоком темпе аккреции, достаточном для стационарного термоядерного горения (М > 10−7М (c)/год для белых карликов) энерговыделение за счет термоядерных реакций будет пренебрежимо мало для нейтронных звезд, и будет определять энергетику тесной двойной системы в случае белых карликов. Такие тесные двойные системы с белыми карликами, на поверхности которых происходит квазистационарное термоядерное горение водорода, возможно, наблюдаются как сверхмягкие рентгеновские источники.

При меньших темпах аккреции горение будет взрывным, квазипериодически повторяющимся по мере накопления критической массы. Такие термоядерные взрывы фиксируются в тесных двойных системах с нейтронными звездами и относительно низким темпом аккреции (М ~ 10~9М (c)/год) как рентгеновские вспышки длительностью около 10 секунд, при которых рентгеновский поток возрастает примерно в 10 раз. Такие источники носят название вспыхивающих рентгеновских источников или барстеров. Термоядерные взрывы на поверхности белых карликов хорошо известны как вспышки новых и повторных новых звезд. При этом их яркость увеличивается в тысячи раз.

Интересно отметить, что вспышки в системах с компактными объектами могут возникать и вследствие изменения темпа аккреции на него из-за ионизационной неустойчивости внешних областей аккреционного диска, где водород ионизован лишь частично. При этом наблюдается прямо противоположная картина (по сравнению с термоядерными вспышками) для систем с нейтронными звездами или черными дырами, и с белыми карликами. Вспышки в катаклизмических переменных не очень значительны (яркость увеличивается в несколько десятков раз) и такие системы носят название карликовых новых звезд, в то время как для рентгеновских двойных светимость в подобной вспышке увеличивается в тысячи раз, и такие явления носят название рентгеновских новых звезд (мягкие рентгеновские транзиенты).

Теоретическое моделирование спектров излучения и кривых блеска качественно описанных выше объектов должно помочь в решении ряда важных задач. Прежде всего, это подтверждение правильности наших представлений о происходящих процессах. Тогда из сравнения рассчитанных и наблюдаемых спектров оказывается возможным определить параметры компактных объектов.

Наиболее актуальной в настоящее время представляется проблема определения масс и радиусов нейтронных звезд. Плотность вещества в их центральных областях превышает ядерную. Свойства вещества при таких условиях, выражаемые через связь между давлением и плотностью (уравнение состояния) невозможно установить экспериментально. В настоящее время существует целый спектр теоретических уравнений состояния вещества в недрах нейтронных звезд, от самых простых до экзотических, рассматривающих возможное образование в таких условиях тяжелых частиц — гиперонов. Не исключен переход вещества в кварковое состояние, когда кварки перестают быть связанными в адронах и становятся свободными (странные или кварковые звезды). Каждое из таких теоретических уравнений состояния предсказывает свою зависимость масса-радиус для нейтронных звезд. Поэтому одновременное определение массы и радиуса какой-либо нейтронной звезды позволит если и не определить точно, то хотя бы ограничить выбор между возможными уравнениями состояния. Все нейтронные звезды — радиопульсары в двойных системах, для которых возможно точное измерение масс, имеют массы, достаточно близкие к 1.4 Mq. Это дает основание надеяться, что большинство нейтронных звезд имеют массы, близкие к этой. Поэтому даже ограничение отношения M/R для нейтронных звезд дает значимую информацию об их радиусах. Одним из возможных способов получить отношение M/R является изучение спектров пограничных слоев в модели слоя растекания на поверхности нейтронной звезды. Результаты такого исследования представлены в диссертации.

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

Для сверхмягких рентгеновских источников определение масс и радиусов горячих белых карликов важно вдвойне. Во-первых, это позволит проверить правильность модели таких источников и определить их эволюционный статус. В частности, важен вопрос о том, являются ли белые карлики в сверхмягких источниках молодыми и горячими (с высокой внутренней температурой) объектами, или термоядерное горение происходит на поверхности остывших белых карликов. Во-вторых, сверхмягкие рентгеновские источники являются кандидатами в предшественники сверхновых типа SN 1а. Поэтому важно знать распределение по массам белых карликов в этих источниках для оценки того, какая их часть может накопить достаточное количество массы, чтобы превысить предел Чандрасекара.

Инструменты рентгеновских обсерваторий, работающих в мягком диапазоне спектра 40 — 100 А, в частности Chandra, калибруются в том числе и по излучению одиночных горячих водородных белых карликов. В последние годы появились указания на то, что спектр таких белых карликов может быть искажен под влиянием эффекта Комптона при рассеянии на свободных электронах. Поэтому важно путем точных расчетов проверить это предположение, чтобы подтвердить или подвергнуть сомнению точность калибровки инструментов рентгеновских обсерваторий в мягком рентгеновском диапазоне.

Хотя определение параметров черных дыр также важно, при моделировании излучения аккреционных дисков на первый план выходит проверка правильности и уточнение наших представлений о физике дисковой аккреции. Здесь можно выделить два аспекта. Стандартная теория дисковой аккреции является параметрической, тензор вязких напряжений в веществе диска считается пропорциональным давлению в диске (полному или газовому, в зависимости от вариантов), Wrcp — аР. Поэтому определение из наблюдений параметра, а крайне необходимо, поскольку позволит проверить теории, выводящие вязкость в диске «из первых принципов», в частности механизм магнито-ротационной неустойчивости (Балбус и Хоули 1991). Кроме того, наблюдениями установлено, что толщина внешних частей аккреционных дисков в ряде систем в несколько раз больше, чем это следует из стандартной теории. Поэтому необходимо найти дополнительные доказательства этому факту и попытаться исследовать возможные причины такого утолщения дисков.

Цель работы.

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

Научная новизна работы.

1. Определены эффективные температуры 10 сверхмягких рентгеновских источников по бланкетированиым моделям атмосфер горячих белых карликов.

2. Определены массы и радиусы 7 горячих белых карликов в сверхмягких рентгеновских источниках с использованием диаграммы «температурасветимость» .

3. Определены параметры атмосфер (Tefj и log д) трех сверхмягких источников, открытых обсерваторией Chandra в галактике М81.

4. Разработан метод расчета спектров слоев растекания на поверхности нейтронных звезд с учетом эффектов ОТО и рассчитана сетка спектров слоев растекания для различных масс и радиусов нейтронных звезд.

5. Определена полоса допустимых значений масс и радиусов нейтронных звезд из сравнения теоретических спектров слоев растекания и наблюдаемых рентгеновских спектров пограничных слоев в маломассивных двойных рентгеновских источниках высокой светимости.

6. Показано, что локальные модели атмосфер оптически толстых аккреционных си-дисков мало отличаются от атмосфер звезд с такими же Теff и log д на поверхности.

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

8. Разработан метод расчетов. абсорбционных спектров самооблучаемых аккреционных дисков в рентгеновских новых звездах. Показано, что абсорбционные компоненты бальмеровских линий водорода должны наблюдаться при низкой относительной светимости диска.

9. Разработан метод расчета оптических и рентгеновских кривых блеска рентгеновских новых звезд на основе теории нестационарной дисковой аккреции Липуновой и Шакуры (2000, 2002) с учетом самооблучения и эффектов ОТО.

10. Из сравнения наблюдаемых кривых блеска рентгеновских новых Единорога 1975 г. и Мухи 1991 г. с теоретическими получены оценки значения параметра Керра, а для черных дыр 0.3 — 0.8 и < 0.2 соответственно) и величины параметра, а для аккреционных дисков (0.6 — 0.9 и 0.4 — 0.5) в этих системах.

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

12. Определены массы 14 белых карликов в промежуточных полярах с использованием спектров в жестком рентгеновском диапазоне, полученных обсерваториями ЯХТЕ/НЕХТЕ и INTEGRAL.

Достоверность полученных результатов.

Разработанные программные комплексы для моделирования атмосфер горячих звезд тестировались путем сравнения с результатами расчетов других авторов, в частности, с моделями Куруца и для белых карликовс моделями Т. Рауха (Университет г. Тюбингена, Германия).

Методика расчетов моделей атмосфер с учетом эффекта Комптона тестировалась путем сравнения с результатами расчетов атмосфер вспыхивающих нейтронных звезд — барстеров, выполненных Ю. Мадеем (Коперниковский центр, Варшава, Польша). Влияние эффекта Комптона на спектры водородных белых карликов рассчитывалось независимо с помощью программных комплексов автора и Мадея.

Результаты аппроксимации спектров сверхмягких источников сравнивались с результатами других авторов. Полученные массы белых карликов для промежуточных поляров сравнивались с массами, полученными другими авторами тем же методом, но только по относительно мягкой части рентгеновского спектра (3−20 кэВ).

Метод расчета кривых блеска рентгеновских новых тестировался на методе расчетов без учета эффектов ОТО и самооблучения, разработанной Липуновой и Шакурой (2002).

Во всех случаях, когда это было возможно, результаты расчетов сравнивались с доступными аналитическими моделями.

Научная, практическая и методическая ценность работы.

В диссертации рассчитаны спектры излучения слоев растекания на поверхности нейтронных звезд и получено ограничение на отношение M/R для нейтронных звезд из сравнения теоретических спектров с наблюдаемыми рентгеновскими спектрами пограничных слоев в маломассивных рентгеновских двойных высокой светимости. Полученное ограничение позволяет сделать вывод, что уравнение состояния вещества в недрах нейтронных звезд является жестким, а это свидетельствует против существования странных (кварковых) звезд. Этот вывод имеет большое научное значение для физики сверхплотного состояния материи.

Получены оценки масс белых карликов в ряде тесных двойных систем (промежуточные поляры и сверхмягкие рентгеновские источники). Эти данные позволят в дальнейшем уточнить теорию поздних стадий эволюции звезд и уже используются рядом исследователей. Научное значение имеет также оценка размеров белых карликов в сверхмягких рентгеновских источниках, что позволит уточнить их эволюционный статус. Методическое значение имеет впервые применённый для сверхмягких источников метод оценки масс белых карликов по их положению на диаграмме «температура — светимость». Предложенная формула для связи максимально возможной температуры белого карлика в сверхмягком рентгеновском источнике с его массой имеет практическое значение, т.к. позволяет получить ограничение снизу на массу белого карлика только из значения эффективной температуры. Последняя величина определяется относительно легко.

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

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

Созданы комплексы компьютерных программ для расчета моделей атмосфер горячих компактных звезд и аккреционных дисков, которые могут быть использованы при анализе спектров излучения таких объектов. Часть из них используется в настоящее время в КГУ, ГАИШ МГУ и САО РАН.

Полученные в работе результаты могут найти применение во всех астрономических учреждениях, где изучаются аккрецирующие компактные объекты, в частности, в КГУ, ГАИШ МГУ, ИКИ РАН, АКЦ ФИАН, ФТИ РАН (г. Санкт-Петербург), САО РАН, ГАО РАН, в Одесском университете, и других зарубежных институтах и обсерваториях.

Результаты, полученные при подготовке диссертации, частично использованы для создания и модернизации лекционных курсов и практических занятий, проводимых диссертантом на кафедре астрономии КГУ.

Основные положения, выносимые на защиту.

1. Метод расчета спектров слоев растекания на поверхности нейтронных звезд и ограничения на размеры нейтронных звезд (R = 15 ± 1.5 км при М = 1.4М©и солнечном химсоставе аккрецируемого вещества), полученные из сравнения теоретических спектров слоев растекания с наблюдаемыми спектрами пограничных слоев в маломассивных рентгеновских двойных высокой светимости.

2. Расчет моделей аккреционных а-дисков и спектров их излучения методами моделей звездных атмосфер с учетом реальной непрозрачности. Показано, что локальные (на данном радиусе) спектры оптически толстых аккреционных дисков практически не отличаются от спектров звезд с такими же Теg и log д на поверхности, если внешнее облучение на данном радиусе диска незначительно.

3. Метод расчета кривых блеска рентгеновских новых звезд с учетом эффектов ОТО и самооблучения диска на основе теории нестационарных аккреционных дисков. Применение метода к описанию рентгеновских и оптических кривых блеска Новой Единорога 1975 г. и Новой Мухи 1991 г. и оценка параметра, а в их аккреционных дисках (0.6 — 0.9 и 0.4 — 0.5 соответственно). Показано, что черная дыра в системе Новой Единорога обладает значительным угловым моментом (а ~ 0.3 — 0.8).

4. Предложен новый механизм переработки мягкого рентгеновского излучения в оптическое за счет многократного рассеяния между оптически толстыми облаками газа над аккреционным диском. Этот механизм способен объяснить высокую оптическую светимость сверхмягких рентгеновских источников.

5. Определены массы и размеры горячих белых карликов в семи сверхмягких рентгеновских источниках, полученные из анализа рентгеновских наблюдений (обсерватория ROSAT). Определены массы белых карликов в 14 промежуточных полярах из сравнения жестких рентгеновских спектров, наблюдаемых обсерваториями RXTE и ИНТЕГРАЛ, с рассчитанными теоретическими спектрами аккреционных колонок.

6. Показано, что учет некогерентности электронного рассеяния не влияет на выходящие спектры излучения водородных (DA) белых карликов в мягком рентгеновском диапазоне. Рассчитанные спектры моделей атмосфер DA белых карликов могут использоваться для калибровки инструментов рентгеновских обсерваторий.

Апробация работы.

Полученные в диссертации результаты представлялись на международных научных конференциях «Современные проблемы звездной эволюции» (Москва, 1998), «Астрофизика высоких энергий сегодня и завтра» (Москва, 2001, 2003, 2005), «Вселенная ИНТЕГРАЛа» (Мюнхен, Германия, 2004), «Жизнь нейтронных звезд» (Санкт-Петербург, 2005), «Изолированные нейтронные звезды» (Лондон, 2006), Собрании Европейского астрономическогообщества JENAM-2000 (Москва, 2000), Всероссийских астрономических конференциях (Санкт-Петербург, 2001; Москва 2004), международных рабочих семинарах «Нейтронные звезды и пульсары» (Бад-Хоннеф, Германия, 2006), «Лебедь Х-2+/-1» (Киттила, Финляндия, 2005), и «Проблемы современной Астрофизики» (Оулу, Финляндия, 2006) — на совещаниях рабочей группы «Звездные атмосферы» при Астросовете РАН (САО РАН, 1992; Киев, 1994), на совещании «Спектроскопические методы в современной астрофизике» (Москва, 2006), на ежегодных Итоговых конференциях КГУ (1990 — 2005), на семинарах Института астрофизики Общества Макса Планка (Гархинг, Германия, 1997), Института астрономии и астрофизики Мюнхенского университета (1997), центра космических полетов им. Маршалла НАСА (Хантсвилл, США, 1999, 2000), университета Оулу (Финляндия, 2003, 2004), кафедры астрономии КГУ (1990 — 2004).

Личный вклад автора.

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

Список публикаций по теме диссертации.

1. Сулейманов В. Ф., Липунова Г. В. и Шакура Н. И. Толщина аккреционных а-дисков: теория и наблюдения, Астрон. Журн. 84, 612−626 (2007).

2. Ревнивцев М. Г., Лутовинов А. А., Сулейманов В. Ф., Мольков С. В., и Сюняев Р. А. Широкополосный рентгеновский спектр промежуточного поляра V2400 Oph, Письма в Астрон. Журн. 30, 848−854 (2004).

3. Ибрагимов А. А., Сулейманов В. Ф., Вихлинин А. А., и Сахибуллин Н. А. Сверхмягкие рентгеновские источники. Параметры атмосфер, Астрон. Журн. 80, 212−222 (2003).

4. Сулейманов В. Ф. и Ибрагимов А. А. Сверхмягкие рентгеновские источники. Фундаментальные параметры, Астрон. Жури. 80, 223−231 (2003).

5. Сулейманов В. Ф., Гош К., Остин Р., и Рамсай Б. Рентгеновское излучение аккреционных дисков в активных ядрах галактик, Письма в Астрон. Журн. 28, 829−839 (2002).

6. Сахибуллин Н. А., Сулейманов В. Ф., Шиманский В. В., и Сулейманова СЛ. Формирование бальмеровских линий в спектрах рентгеновских новых звезд Письма в Астрон. Журн. 24, 26−38 (1998).

7. Сулейманов В. Ф. Могут ли наблюдаться оптические линии поглощения в спектрах рентгеновских новых звезд?, Письма в Астрон. Журн. 22, 107−123 (1996).

8. Сулейманов В. Ф. и Шакура Н. И. Эффекты отражения в катаклизмических переменных звездах. V603 Aql, Письма в Астрон. Журн. 20, 28−31 (1994).

9. Сулейманов В. Ф. Возможности программного комплекса STARDISK, Изв. АОЭ, N58, 11−26 (1994).

10. Сулейманов В. Ф. Моделирование аккреционных дисков и спектров катаклизмических переменных звезд. I. V603 Aql, Письма в Астрон. Журн. 18, 255−265 (1992).

11. Сулейманов В. Ф. О возможной природе источника рентгеновского излучения в квазарах, Письма в Астрон. Журн. 17, 575−582 (1991).

12. Suleimanov V., Werner К. Importance of Compton scattering for radiation spectra of isolated neutron stars with weak magnetic field, Astron. Astrophys. 466, 661−666 (2007).

13. Suleimanov V., Madej J., Dracke J.J., Rauch Т., Werner K. On the relevance of Compton scattering for the soft X-ray spectra of hot DA white dwarfs, Astron. Astrophys. 455, 679−684 (2006).

14. Suleimanov V., Poutanen J. Spectra of the spreading layers on the neutron star surface and constraints on the neutron star equation of state, MNRAS 369, 2036;2048 (2006).

15. Suleimanov V., Revnivtsev M, Ritter H. RXTE broadband X-ray spectra of intermediate polars and white dwarf mass estimates, Astron. Astrophys. 435, 191−199 (2005).

16. Falanga M., Bonnet-Bidaud J.M., Suleimanov V. INTEGRAL broadband, X-ray spectrum of the intermediate polar V709 Cassiopeiae, Astron. Astrophys. 444, 561−564 (2005).

17. Revnivtsev M., Lutovinov A., Suleimanov V., Sunyaev R., Zheleznyakov V. Broadband X-ray spectrum of intermediate polar V1223 Sgr, Astron. Astrophys. 426, 253−257 (2004).

18. Suleimanov V., Shakura N., Lipunova G. Light curve modeling for time-dependent accretion disks in X-ray novae with General Relativity effects taken into account in: 5-th INTEGRAL Workshop Proceedings on «The INTEGRAL Universe», pp.403−407 (2004).

19. Lipunova G., Suleimanov V. Radial Dependences of Physical Parameters in alpha-disks as a Consequence of Two Vertical Structure Solutions, Baltic Astronomy 13, 167−170 (2004).

20. Suleimanov V., Meyer F., Meyer-Hofmeister E. High efficiency of soft X-ray radiation reprocessing in supersoft X-ray sources due to multiple scattering, Astron. Astrophys. 401, 1009−1015 (2003).

21. Swartz D., Ghosh K., Suleimanov V., Tennant A., Kinwan Wu CHANDRA discovery of luminous supersoft X-ray sources in M81, Astrophys. J. 574, 382−397 (2002).

22. Suleimanov V., Meyer-Hofmeister E., Meyer F. Reprocessing of soft X-ray radiation in supersoft X-ray sources, Astron. Astrophys. 350, 63−72 (1999).

23. Suleimanov V., R. Rebolo R. The Li I line in spectra of X-ray Novae accretion discs, Astron. Astrophys. 330, 999−1004 (1998).

Заключение

.

Представленная диссертация посвящена, в основном, моделированию спектров излучения аккреционных течений на компактные объекты: белые карлики, нейтронные звезды и черные дыры. Рассматривались почти все возможные виды аккреционных течений: аккреционные диски, пограничные слои в модели слоя растекания, и аккреционные колонки. Моделирование производилось на основе построения моделей атмосфер излучающих аккрецонных течений с учетом того, что энергия генерируется, как правило, во всем объеме течения. Для построения моделей атмосфер использовались методы и модифицированные программные комплексы, предназначенные для моделирования атмосфер звезд. В том случае, когда аккреционное течение нельзя охарактеризовать единой эффективной температурой (например, аккреционный диск или слой растекания), оно разбивалось на некоторое количество колец, и строились локальные модели атмосфер каждого кольца. Интегральный спектр суммировался из локальных спектров колец.

В результате моделирования был сделан вывод о том, что для оптически толстых аккреционных течений высокой светимости локальные атмосферы и спектры мало отличаются от моделей и спектров излучения звездных атмосфер с такими же значениями эффективной температуры и ускорения силы тяжести на поверхности. Связано это с тем, что выходящий спектр излучения формируется в относительно тонком поверхностном слое течения, где генерируется достаточно малая часть общей энергии. Поэтому появление дополнительного потока излучения (по сравнению с классической звездной атмосферой, в которой вообще нет генерации энергии и поток излучения постоянен) не оказывает заметного влияния на локальную атмосферу и спектр ее излучения.

Сравнение рассчитанных спектров излучения с наблюдаемыми рентгеновскими спектрами позволило получить оценки масс и радиусов 7 горячих белых карликов в сверхмягких рентгеновских источниках. Полученные величины согласуются с моделью (квази-) устойчивого термоядерного горения на поверхности. Кроме того, получены оценки масс 14 белых карликов в промежуточных полярах. Однако, наиболее интересным следует признать получение зависимости масса-радиус для нейтронных звезд из сравнения теоретических спектров слоев растекания с наблюдаемыми рентгеновскими спектрами пограничных слоев в маломассивных рентгеновских двойных высокой светимости. Из нее, в частности, следует, что массе нейтронной звезды 1.4 М©соответствует радиус ~ 15 км. Этот результат свидетельствует в пользу жесткого уравнения состояния нейтронных звезд и против существования кварковых звезд.

Проведено сравнение рассчитанных спектров излучения аккреционных дисков с наблюдаемыми оптическими и ультрафиолетовыми спектрами избранных катаклизмических переменных и рентгеновских новых звезд. Было подтверждено наличие избыточного оптического излучения в спектре повоподобной V603 Орла. Такое же избыточное излучение, но намного более значительное, наблюдается в спектрах сверхмягких рентгеновских источников. Оно может быть объяснено только переработкой жесткого излучения центрального источника и/или внутренних частей диска внешними частями аккреционного диска. Для этого требуется увеличить видимую толщину внешних частей диска. Моделирование рентгеновских и оптических кривых блеска двух рентгеновских новых звезд (Новой Единорога 1975 г. и Новой Мухи 1991 г.) подтверждает увеличенную толщину аккреционных дисков в этих системах. Предложена гипотеза, что такое утолщение может быть вызвано наличием вещества над диском, вероятнее всего — короны. Для объяснения высокой эффективности переработки мягкого рентгеновского излучения в оптическое в сверхмягких рентгеновских источниках предложена гипотеза неоднородной короны, с погруженными в нее относительно холодными оптически толстыми облаками. Многократное рассеяние излучения между облаками позволяет увеличить коэффициент переработки до 0.3−0.5.

Подводя итог, можно заключить, что в диссертации разработаны и программно реализованы методы расчета моделей атмосфер и спектров излучения аккреционных течений на компактные объекты. Сравнение рассчитанных и наблюдаемых спектров позволило оценить параметры ряда компактных объектов и исследовать физику дисковой аккреции вблизи внешнего радиуса аккреционных дисков.