Astronet Астронет: П. В. Бакланов/ГАИШ Теоретический каталог сверхновых второго типа
http://variable-stars.ru/db/msg/1174722/node3.html
Теоретический каталог сверхновых второго типа
<< 1. Моделирование взрыва. | Оглавление | 3. Заключение. >>

Разделы


2. Результаты моделирования.

2.1 Полученные модели для наблюдателя.

Начиная эту работу, для теоретического каталога планировалось расчитать 27 моделей, т.е. сетку , где масса принимает значения 8.41, 16.41, 20.41 , радиус -- 50, 500, 5000 , а энергия -- , , эрг. Такое большое число моделей требовало приличного количества времени для счета, но зато обеспечивало любую комбинацию из значений массы, радиуса и энергии. Выполнение работы шло активно, поэтому появилась возможность расчитать еще несколько моделей. Добавлены были модели с массой и энергией эрг. Эти модели давали возможность перекрыть диапазон значений параметров из работ ([3],[4]) и проводить сравнение полученных теоретических соотношений более корректно.

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

2.2 Пример работы с каталогом.

Покажем на примере как работать с каталогом. Для примера выбрана хорошо известная и яркая кривая блеска сверхновой SN1969L ([15]). Построенную кривую блеска SN1969L в трех! фильтрах сравнил с кривыми блеска из каталога (рис.2.2) ( учтя межзвездное поглошение, для SN1969L , , ) и нашел, что наиболее близкая теоретическая кривая у модели r05m3ni2, т.е. предсверхновая была сверхгигант с , , .

На примере другой сверхновой 1999em можно показать влияние никеля на кривую блеска. Неплохое согласие кривых блеска для модели r500m3ni2E10 (рис.2.2) говорит нам, что предсверхновая SN1999em была сверхгигант с , , . Начиная от 80 дней и примерно до 140 дней на теоретической кривой блеска виден горб, вызванный переработкой жесткого излучения, образующегося при распаде . Полная масса во всех моделях была , типичная для "средней" сверхновой II типа. Для SN1999em подобная деталь отсутствует, что позволяет предположить обедненность в SN1999em. Таким образом можно судить о хим. составе предсверхновой.

Рис.: Кривая блеска SN 1969L в фильтрах U,B,V треугольники, квадраты и кружки соответственно. Сплошные линии - численный расчет. Модель r05m3ni2.

Рис.: Кривая блеска SN 1999em в фильтрах U,B,V крестики, квадраты и кружки соответственно. Сплошные линии - численный расчет. Модель r500m3ni2E10.

2.3 Рекомендации по сравнению кривых блеска.

Наилучшее согласие с реальными кривыми блеска достигается в первые 4-6 недель. Поэтому для поиска подходящей кривой надо ориентироваться на этот интервал времени. Дальнейший ход кривой блеска определяется количеством никеля , который проявляется ввиде горба на кривой блеска.

При сравнении надо учитывать межзвездное поглошение.

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

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

2.4 Полученные модели для теоретика.

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

Далее используются обозначения:
масса предсверхновой, в массах Солнца
энергия предсверхновой, в эрг
радиус фотосферы предсверхновой, в радиусах Солнца
длительность стадии плато на кривой блеска, в сут.
абсолютная зв. величина в фильтре V на середине плато кривой блеска
скорость вещества фотосферы в момент, соответствующий середине плато кривой блеска, в см.
(Для упрощения сравнения все величины приводятся по анологии со статьей [4]

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

В таблице (2.1) собраны характеристики полученных моделей. Не все модели вошли в таблицу, так как на некоторых модельных кривых блеска невозможно было отождествить плато.


Таблица. Параметры полученных моделей
No Имя модели M () R () E,
0 R5000_M20_Ni008_E3 15.63 5000 3 87 -16.60 1.98
1 R5000_M16_Ni008_E3 11.63 5000 3 108 -16.50 1.84
2 r5000m8ni2E10 3.63 5000 10 62 -17.50 3.72
3 R5000_M16_Ni008_E17 11.63 5000 17 68 -18.00 4.32
4 R5000_M8_Ni008_E3 3.63 5000 3 84 -16.70 2.54
5 R5000_M8_Ni008_E1 3.63 5000 2 57 -17.70 2.16
6 r05m3ni2 15.63 5000 17 90 -18.00 3.45
7 R500_M16_Ni008_E17 11.63 500 17 78 -16.50 3.76
8 R5000_M8_Ni008_E17 3.63 5000 17 57 -17.70 4.37
9 R50_M16_Ni008_E17 11.63 50 17 33 -17.00 10.17
10 r500m3ni2E10 15.63 500 10 110 -16.20 2.68
11 R500_M8_Ni008_E3 3.63 500 3 70 -15.90 2.24
12 r500m3ni2 15.63 500 17 92 -16.00 3.44
13 R5000_M20_Ni008_E10 15.63 5000 10 95 -17.40 2.91
14 R500_M16_Ni008_E3 11.63 500 3 121 -15.00 1.78
15 R5000_M16_Ni008_E10 11.63 5000 10 82 -17.40 3.11
16 R500_M16_Ni008_E10 11.63 500 10 89 -15.80 2.93

С помощью метода МНКа были обработаны данные из таблицы (2.1). В результате получены теоретические соотношения, связывающие V, , с M, R, E. Однако для М и Е ошибка параметров превысила сами значения, поэтому формулы для М и Е не приводятся. Отличие logM от истиного значения даже на 50% приводит к ошибке в 4 раза, а знание массы с точностью до порядка бессмыслено(тоже для Е).

Только для радиуса получилась надежная аппроксимация:


Максимальное отклонение от апроксимационной кривой составило 35 процентов. Таким образом можно ожидать, что оценка R по полученному соотношению не превысит 35 процентов.

Кривые блеска, расчитанные с помощью STELLA, ближе к реальным, чем аналогичные из работ ([4],[3]), причем во всех трех фильтрах. В программе STELLA учтенно много нового, чего не было в упомянутых выше работах. Так, например, нестационарность процесса выхода ударной волны на поверхность, гидродинамика, корректный учет поглощения при расширении и многое др. Не удивительно, что в соотношениях (2.3) и аналогичном из работы ([4]) коэффициенты при , и отличаются друг от друга.


<< 1. Моделирование взрыва. | Оглавление | 3. Заключение. >>

Rambler's Top100 Яндекс цитирования