Актуальность темы
диссертации
Нейтронные звезды (НЗ) это компактные астрофизические объекты с массой порядка солнечной, и радиусом около 10 км. Плотность вещества в центре звезды может доходить до недостижимых в лабораторных условиях значений в несколько ядерных плотностей ро (ро ~ 2.8×1014 г см-3). Свойства такого вещества (состав, уравнение состояния, критические температуры перехода барионов в сверхтекучее состояние) можно исследовать, в частности, сравнивая теоретические модели пульсирующих НЗ с данными наблюдений. Колебания НЗ могут приводить к модуляции электро-магнитного излучения с поверхности, а также к излучению гравитационных волн. В настоящее время уже наблюдаются высокочастотные осцилляции электромагнитного излучения во время гигантских рентгеновских вспышек [1−4], которые интерпретируются как результат колебаний НЗ. Кроме того, проектируются детекторы, которые, согласно оценкам, смогут зафиксировать гравитационное излучение от колеблющихся одиночных НЗ [5−8]. Для интерпретации имеющихся (и будущих) наблюдательных данных необходимо иметь хорошо развитую теорию пульсаций НЗ. Разработку такой теории осложняет тот факт, что нуклоны и гипероны во внутренних слоях нейтронных звезд могут находиться в сверхтекучем состоянии. Тем самым, для описания пульсирующей нейтронной звезды необходимо понимать процессы, происходящие в сверхтекучих, сильно неидеальных смесях ферми-жидкостей при отличных от нуля температурах. Этой задаче и посвящена настоящая диссертация. В частности, в диссертации исследованы динамические явления и диссипативные процессы в сверхтекучих НЗ, внесен вклад в разработку теории пульсаций НЗ с учетом сверхтекучести барионов в их ядрах.
Цели работы
1. Исследование динамических свойств сверхтекучего нуклон-гиперонного вещества ядер массивных НЗ и диссипативных процессов в этом веществе.
2. Расчет для этих целей коэффициентов эффективной объемной вязкости релятивистской сверхтекучей нуклон-гиперонной смеси, а также расчет для того же состава вещества релятивистской матрицы сверхтекучих плотностей основного параметра релятивистской гидродинамики сверхтекучих смесей.
3. Исследование влияния температурных эффектов на спектры пульсаций сверхтекучих НЗ.
4. Разработка приближенного метода расчета спектров пульсаций сверхтекучих НЗ в рамках общей теории относительности (ОТО).
Новизна работы
Большинство результатов, полученных в работе, носит приоритетный характер. В частности:
1. Впервые сформулирована диссипативная релятивистская гидродинамика сверхтекучего нуклон-гиперонного вещества, состоящего из нейтронов, протонов, электронов, мюонов, а также, А и ?~ гиперонов. Показано, что она содержит шестнадцать коэффициентов объемной вязкости.
2. Для нуклон-гиперонного вещества впервые рассчитаны все шестнадцать коэффициентов эффективной объемной вязкости, генерируемые неравновесными реакциями взаимного превращения частиц.
3. Для того же состава вещества впервые рассчитана релятивистская матрица сверхтекучих плотностей, являющаяся необходимым элементом при рассмотрении динамических процессов в сверхтекучих НЗ. Опеределена ее температурная зависимость.
4. Впервые в рамках релятивистской гидродинамики сверхтекучих смесей исследованы динамические и диссипативные свойства сверхтекучего нуклон-гиперонного вещества. В частности, рассчитаны скорости всех трех существующих в таком веществе звуковых мод. Вычислены времена затухания звуковых мод за счет сдвиговой и эффективной объемной вязкости.
5. Разработан новый приближенный метод, позволяющий кардинально упростить расчеты спектров пульсирующих сверхтекучих НЗ в рамках ОТО. С использованием этого метода сделан ряд важных выводов о свойствах сверхтекучих пульсационных мод. В частности, впервые показано, что гравитационное излучение от сверхтекучих мод должно быть подавлено по сравнению с нормальными модами. Кроме того, показано, что сверхтекучие моды практически не выходят на поверхность НЗ. Тем самым, детектирование этих мод может быть затруднено.
6. Впервые изучено влияние двух основных температурных эффектов на спектр пульсаций сверхтекучих НЗ. Рассмотрение проведено на примере радиальных пульсаций. Учитывалась как зависимость от температуры релятивистской матрицы сверхтекучих плотностей, так и температурная зависимость размеров сверхтекучей области.
Показано, что модели, учитывающие влияние температуры только на размер сверхтекучей области (приближение, использованное в некоторых современных работах, например в [9−11]), дают нереалистичные спектры.
Достоверность научных результатов
Результаты диссертации получены аналитически или численно при помощи компьютера. Их достоверность подтверждается использованием адекватных математических и численных методов в рамках физически разумных приближений, применимость которых ограничена четко сформулированными критериями. Там, где это было возможно, результаты численных расчетов сверялись с теоретическими формулами, полученными в предельных случаях. Результаты работы опубликованы в реферируемых международных журналах, представлялись на всероссийских и международных конференциях, а также на семинарах сектора теоретической астрофизики ФТИ им. А. Ф. Иоффе (Санкт-Петербург).
Практическая значимость работы
Результаты диссертации имеют большое значение для исследования фундаментальных свойств вещества сверхъядерной плотности в ядрах НЗ, для численного моделирования колебаний сверхтекучих НЗ, адекватной интерпретации наблюдений пульсирующих НЗ и определения целей для будущих поколений детекторов гравитационного излучения.
Основные положения, выносимые на защиту
1. Формулировка диссипативной релятивистской гидродинамики сверхтекучего нуклон-гиперонного вещества, состоящего из нейтронов, протонов, электронов, мюонов, а также Л и Египеронов.
2. Расчет коэффициентов эффективной объемной вязкости, генерируемой неравновесными реакциями взаимного превращения частиц в сверхтекучих нуклон-гиперонных ядрах массивных нейтронных звезд.
3. Самосогласованный расчет релятивистской матрицы сверхтекучих плотностей нуклон-гиперонной смеси в рамках релятивистской теории ферми-жидкости Ландау, а также с использованием сгс^р-модели среднего поля. Определение температурной зависимости этой матрицы.
4. Исследование звуковых мод в сверхтекучем нуклон-гиперонном веществе ядер массивных нейтронных звезд. о. Разработка нового приближенного метода для расчета пульсаций сверхтекучих нейтронных звезд в рамках общей теории относительности.
6. Анализ влияния температурных эффектов на спектр пульсаций сверхтекучих нейтронных звезд (на примере радиальных пульсаций).
Апробация работы и публикации
Результаты, вошедшие в диссертацию, были получены в период с 2006 по 2010 гг. и изложены в 16 печатных работах [12−27] (включая 5 статей в реферируемых журналах [12−16]). Результаты диссертационной работы представлялись и обсуждались на всероссийских и международных конференциях: Всероссийская конференция «Астрофизика высоких энергий сегодня и завтра» (Москва, 2007, 2008, 2010), Всероссийская конференция «Физика нейтронных звезд» (Санкт-Петербург, 2008), International Workshop «Compstar 2009: the crust of compact stars and beyond» (Coimbra, 2009), Всероссийская астрономическая конференция (CAO, 2010), Всероссийская конференция по физике и астрономии для молодых ученых Санкт-Петербурга и Северо-запада (Санкт-Петербург, 2010), International Workshop «MODE-SNR-PWN» (Bordeaux, 2010) — семинарах сектора теоретической астрофизики ФТИ им. А. Ф. Иоффе (Санкт-Петербург).
7.8 Выводы
В этой главе в рамках ОТО было изучено влияние температурных эффектов на пульсации сверхтекучих нейтронных звезд. Был рассмотрен наиболее простой случай — радиальные пульсации невращающейся нейтронной звезды, состоящей из нейтронов, протонов и электронов. Было рассмотрено два температурных эффекта. Первый эффект связан с тем, что релятивистская матрица сверхтекучих плотностей Угк, входящая в уравнения сверхтекучей гидродинамики, зависит от температуры (этот эффект обсуждался в работе СА06). Второй эффект заключается в том, что размер сверхтекучей области падает с ростом температуры (см. работы [9−11]).
Было рассмотрено четыре модели, в которых учитывался один, или оба температурных эффекта и рассматривались различные профили критической температуры нейтронов. Для этих моделей были рассчитаны спектры радиальных пульсаций. Было показано, что во всех моделях существует два типа колебаний — нормальные и сверхтекучие. Частоты нормальных колебаний практически совпадают с частотами несверхтекучей звезды. Частоты сверхтекучих колебаний сильно зависят от температуры, и имеют различное поведение для каждой модели. Было продемонстрировано, что в модели, учитывающей влияния температуры только на размер сверхтекучей области, но не на матрицу сверхтекучих плотностей, получаются качественно неверные спектры.
Используя приближенные уравнения описывающие пульсации сверхтекучей нейтронной звезды, полученные в работе [16], мы нашли асимптотики для собственных частот при температурах близких к критической температуре нейтронов [Т°° —"¦ Т^(0)], которые хорошо согласуются с численными результатами.
Кроме того, были проанализированы собственные функции для первых шести мод колебаний. Было показано, что в сверхтекучих модах практически не возмущается ток барионов, поэтому они локализованы в ядре и не практически выходят на поверхность. Как следствие, они вряд ли могут быть зарегистрированы по наблюдениям модуляции электромагнитного излучения с поверхности звезды.
Также были рассчитаны характерные времена затухания г радиальных пульсаций за счет сдвиговой вязкости и неравновесных урка-процессов и проанализированы асимптотики т при Г°° —> 0). Было показано, что времена затухания рассчитанные в одно-жидкостной гидродинамике несверхтекучей жидкости могут в разы отличаться от времен затухания нормальных мод (хотя в тоже время частоты нормальных мод практически совпадают с частотами несверхтекучей звезды). Было показано, что при переходе звезды (например, при остывании) с нормального типа колебаний на сверхтекучий время затухания уменьшается на 1−3 порядка. Таким образом, сверхтекучие колебания затухают быстрее нормальных.
ЗАКЛЮЧЕНИЕ
1. Сформулирована диссипативная релятивистская гидродинамика сверхтекучего нуклон-гиперонного вещества ядер массивных нейтронных звезд, состоящего из нейтронов, протонов, электронов, мюонов, а также, А и Египеронов. Показано, что она содержит шестнадцать коэффициентов объемной вязкости.
2. В рамках этой гидродинамики для нуклон-гиперонного вещества рассчитаны все шестнадцать коэффициентов эффективной объемной вязкости, генерируемой неравновесными слабыми безлептонными реакциями взаимного превращения частиц. Показано, что в этом случае только три коэффициента являются независимыми.
3. Для того же состава вещества для произвольных температур рассчитана релятивистская матрица сверхтекучих плотностей являющаяся необходимым элементом при рассмотрении динамических процессов в массивных сверхтекучих НЗ. Расчет проведен аналитически и самосогласованно, в рамках релятивистской теории ферми жидкости Ландау, обобщенной на случай смесей. Показано, что матрица Ylk выражается через параметры Ландау f{k и универсальные функции температуры Фг. Численно релятивистская матрица сверхтекучих плотностей рассчитана для одной из релятивистских ашр-моделей среднего поля со скалярным самодействием.
4. На примере звуковых волн впервые исследованы динамические и диссипативные свойства сверхтекучего нуклон-гиперонного вещества. В частности, в рамках релятивистской гидродинамики сверхтекучих смесей показано, что в сверхтекучем нуклон-гиперонном веществе может распространяться 3 звуковые моды. Рассчитаны скорости этих мод, вычислены их времена затухания за счет сдвиговой вязкости и эффективной объемной вязкости, генерируемой неравновесными реакциями. На этом примере показана важность использования релятивистской гидродинамики сверхтекучих смесей при исследовании динамических и диссипативных процессов — времена затухания, рассчитанные в такой гидродинамике и в одножидкостной несверхтекучей гидродинамике (часто используемое приближение) могут отличаться на несколько порядков.
5. Разработай новый приближенный метод расчета спектров пульсаций сверхтекучих релятивистских НЗ, позволяющий существенно упростить вычисления. Метод основан на том, что уравнения, описывающие пульсации таких звезд, можно представить в виде двух слабо связанных систем. Одна из них описывает нормальные моды, другая — сверхтекучие. Показано, что для реалистичных уравнений состояния параметр связи этих систем мал (в ~ 0.01 — 0.05), что позволяет уже в нулевом приближении по й (когда системы считаются полностью независимыми) рассчитывать спектры с точностью в несколько процентов. Это утверждение было проверено на примере радиальных пульсаций. Также было показано, что благодаря такому расщеплению уравнений (г) гравитационное излучение от сверхтекучих мод должно быть существенно подавлено по сравнению с излучением от нормальных мод- (гг) сверхтекучие моды практически не выходят на поверхность НЗ и вряд ли могут являться причиной наблюдаемой модуляции электромагнитного излучения с поверхности звезды [1].
6. Исследовано влияние температурных эффектов на спектр осцилляций сверхтекучих НЗ на примере радиальных пульсаций. При этом учтены два температурных эффекта: зависимость от температуры матрицы сверхтекучих плотностей, и температурная зависимость размеров сверхтекучей области. Сопоставлены различные модели, учитывающие как один, так и оба температурных эффекта. Показано, что модели, учитывающие влияние температуры только на размер сверхтекучей области (приближение использованное в некоторых современных работах, например, [9]), дают нереалистичные спектры. Кроме того, на примере радиальных пульсаций подтверждено, что сверхтекучие моды практически не выходят на поверхность звезды и не возбуждают ток барионов. Рассчитаны времена затухания радиальных пульсаций. Показано, что сверхтекучие моды затухают быстрее нормальных, и время затухания при переходе (например, в ходе тепловой эволюции) с нормального на сверхтекучий тип пульсаций падает на 1−3 порядка.
