№ 11.2Взрыв нечерного тела
Вспышки новых звезд объясняются взрывным синтезом гелия из водорода, происходящим на поверхности белого карлика, являющегося одним из компонент двойной звезды. Пусть вспышка произошла на белом карлике с массой, равной массе Солнца, и радиусом, равным радиусу Земли, в результате чего видимая звездная величина звезды изменилась на 10 и в максимуме блеска составила +10. Известно, что в результате вспышки 1% массы водородной оболочки превратился в гелий, общая энергия, высвеченная при вспышке, составила 1038 Дж, а вспышка продолжалась 11 суток, в течение которых блеск новой практически не менялся. Оцените:
- A. ускорение свободного падения на поверхности белого карлика;
- B. массу взорвавшейся водородной оболочки белого карлика (в массах Солнца);
- C. расстояние до новой звезды;
- D. характерную температуру водородной оболочки непосредственно перед вспышкой;
- E. характерную толщину водородной оболочки непосредственно перед вспышкой.
В реакциях синтеза гелия из водорода 0.7% массы прореагировавшего вещества переходит в энергию. Все болометрические поправки считать равными нулю, межзвездным поглощением пренебречь.
Решение
A. Вращением белого карлика (БК) вокруг оси можно пренебречь: он входит в состав тесной двойной системы (в противном случае аккреция вещества со второго компонента не происходила бы) и в системе должна была произойти приливная синхронизация. Как следствие, период вращения БК вокруг оси примерно совпадает с орбитальным периодом системы, который не может быть малым.
Поэтому ускорение свободного падения можно вычислить просто как
$$g = \frac{G\mathfrak{M}}{R^2},\tag{1}$$
где — гравитационная постоянная, M — масса карлика, — его радиус. Зная массу Солнца и радиус Земли, получаем (здесь и далее в промежуточных вычислениях все численные данные выражены в СИ)
$$g = \frac{6.7 \cdot 10^{-11} \cdot 2 \cdot 10^{30}}{(6.4 \cdot 10^6)^2} \approx 3 \cdot 10^6 \text{ m/c}^2.$$ (2)
B. При превращении 1% массы оболочки Mоб из водорода в гелий должна была выделиться энергия = 0.01Mоб 2 , где = 7 · 10−3 , — скорость света. Тогда
$$\mathfrak{M}_{06} = \frac{E}{\alpha \, c^2 \cdot 0.01} = \frac{10^{38}}{7 \cdot 10^{-3} \cdot (3 \cdot 10^8)^2 \cdot 10^{-2}} \approx 2 \cdot 10^{25} \text{ Kp} = 10^{-5} \, \mathfrak{M}_{\odot}. \tag{3}$$
C. Оценим светимость новой в максимуме вспышки. За время = 11 суток была высвечена энергия , поэтому светимость
$$L = \frac{E}{t} = \frac{10^{38}}{11 \cdot 9 \cdot 10^4} = 10^{32} \text{ B}_{\text{T}}$$ (4)
(с достаточной для нас точностью в сутках 9 · 104 секунд).
Тогда абсолютная звездная величина новой в максимуме вычисляется из соотношения
$$M - M_{\odot} = -2.5 \lg \frac{L}{L_{\odot}} \tag{5}$$
и
$$M = 4^{m}.7 - 2.5 \lg \frac{10^{32}}{3.8 \cdot 10^{26}} \approx -8^{m}.8$$ (6)
(отметим, что она будет болометрической — у нас нет информации о том, в каких диапазонах высвечивалась энергия вспышки; впрочем, по условию задачи разница между оптическими и болометрическими звездными величинами отсутствует).
Для оценки расстояния воспользуемся соотношением = − 5 lg + 5, откуда
$$r = 10^{\frac{m-M+5}{5}} = 10^{(10+8.8+5)/5} \approx 58 \cdot 10^3 \text{ пк} = 58 \text{ кпк}.$$ (7)
Это за пределами Галактики, но странным такой результат не является: вспышки новых успешно наблюдаются и в других галактиках. Эта вспышка, по-видимому, произошла в каком-то из Магеллановых облаков.
D. Напрашивается идея, что непосредственно перед вспышкой БК представлял собой шар, излучающий только за счет того, что он нагрет до некоторой температуры . Поскольку его радиус нам известен, то если мы знаем его светимость 0, то температуру можно будет определить из соотношения 0 = 424 (тут — постоянная Стефана-Больцмана).
Верхнюю возможную оценку светимости 0 можно найти, считая, что видимая звездная величина белого карлика (а значит и абсолютная тоже) изменилась на 10 при вспышке. Фактически это эквивалентно предположению, что излучением второго компонента системы можно пренебречь, что более-менее правдоподобно только для максимума вспышки. Все же воспользуемся им, понимая, что учет наличия второго компонента понизит оценку светимости 0 и, как следствие, найденную оценку температуры.
Такое изменение абсолютной звездной величины соответствует изменению светимости на 4 порядка, откуда 0 = 1028 Вт. Тогда можно вычислить температуру (точнее, также ее верхнюю оценку):
$$T = \sqrt[4]{\frac{L_0}{4\pi\sigma R^2}} = \sqrt[4]{\frac{10^{28}}{4\pi \cdot 5.67 \cdot 10^{-8} \cdot (6.4 \cdot 10^6)^2}} \approx 1.4 \cdot 10^5 \text{ K}.$$ (8)
Однако полученный результат явно нелеп: при этой температуре должны были «запуститься» реакции синтеза гелия из водорода, а она для этого слишком мала (несмотря на то, что мы получили ее оценку сверху). Причина в том, что водородная оболочка
БК, как мы увидим далее, является достаточно тонкой и, как следствие, она прозрачна для собственного излучения. Поэтому модель абсолютно черного тела для описания ее излучения непригодна.
Это означает, что оценку температуры надо получить другим путем. Поскольку температура должна быть достаточной для запуска термоядерных реакций, можно использовать в качестве оценки минимальную температуру в центре звезд главной последовательности — это несколько миллионов кельвинов. Более точной оценка получится, если учесть, что реакции синтеза гелия из водорода в новых — это реакции CN-цикла, для запуска которых нужна большая температура (этот цикл становится основным в звездах главной последовательности с массой, несколько большей солнечной, поэтому и температура должна быть несколько больше, чем температура в центре Солнца). Лучшей оценкой будет $T \approx 2 \cdot 10^7 \ \mathrm{K}$ .
Е. Для оценки толщины оболочки воспользуемся известным выражением для толщины однородной атмосферы (она же экспоненциальный масштаб атмосферы)
$$h_0 = \frac{\Re T}{\mu g},\tag{9}$$
где $\mathfrak{R}$ — универсальная газовая постоянная, $\mu$ — молярная масса газа. Поскольку водород при таких температурах полностью ионизован, молярная масса окажется равной $0.5 \cdot 10^{-3}$ кг/моль, а все остальные данные у нас уже имеются. Тогда
$$h_0 = \frac{\Re T}{\mu g} = \frac{8.3 \cdot 2 \cdot 10^7}{5 \cdot 10^{-4} \cdot 3 \cdot 10^6} \approx 10^5 \text{ M}. \tag{10}$$