
книги из ГПНТБ / Горбацкий, В. Г. Новоподобные и новые звезды
.pdf120 |
ГЛ. V. ОКОЛОЗВЕЗДНЫЕ ОБОЛОЧКИ НОВЫХ |
Угол fl „ |
определяет границу околозвездной оболочки |
(рис. 34).
Рассматривая захват вещества главной оболочкой, будем, как и выше, считать ее вначале, при г = г0, сфери ческой со значением массы т (г0) и скорости v (г0), а так же примем, что к ней присоединяется вещество, состоящее
Кнпфиодтиелш
Рис. 34. Схема околозвездной оболочки.
из невзаимодействующих между собой частиц. Таким об разом, пренебрегаем газодинамическими эффектами, воз никающими при захвате вещества главной оболочкой.
Масса dr-ih (г;й) da, присоединяющаяся в пределах телесного угла dсо к элементу главной оболочки, движу
щемуся под углом ----й к |
орбитальной |
плоскости на |
пути от г до г -|- dr, равна |
|
|
dr т (г;й) da |
= r2p da dr. |
(20.5) |
Для величины массы т (г;й) da, присоединившейся в уг ле da к моменту, когда главная оболочка удалится на рас стояние г, имеем
Г
т (г; й) da — da т (г; й) dr.
Л
Принимая во внимание (18.5) и (20.5), находим
т (г;й) da — {г — га) Ф (й; й 0) da,
где обозначено:
Ф (й; й„) = |
Ai Lg flo (1 — ctg2 ft • tgftp) |
|
sin2 ft |
(21.5)
(22.5)
(23.5)
Введем величину т' (г0) = — |
. Из закона сохрапеиия |
количества движения вытекает равенство: |
|
Ы' (го) -f т (г;-0 )1 г? (г; ft) |
= т (r0) v (г0), (24.5) |
§ 2. ВЛИЯНИЕ ОКОЛОЗВЕЗДНОЙ ОБОЛОЧКИ |
121 |
откуда получается выражение радиальной скорости в на правлении ft:
v (г; ft) = |
________ т ' (го) v (го)_________ |
(25.5) |
|
Ф ("O'; {hi) (г — го) + т' (г0) |
’ |
соответствующее формуле для скорости (12.5) в случае сферической околозвездной оболочки. Согласно (25.5) скорость зависит от направления движения в главной оболочке и формы околозвездной оболочки. Поэтому пер воначально сферическая главная оболочка в результате взаимодействия с околозвездной деформируется. Ее фор му в зависимости от времени, т. е. функцию г (ft; t) нахо дим, интегрируя (24.5). С учетом (22.5) имеем
-у- Ф № ®а) + (г — г0) |
— г0 (г — 4 “ г°) Ф № ^ = |
= ^ P - v ( r 0) ( t ~ t 0). (26.5)
Обозначив количество вещества, поступающего в еди ницу времени в околозвездную оболочку с уровня г0, че рез Qs, нетрудно найти [82] путем интегрирования (18.5),
что в случае, когда угол-^----ft„ мал (т. е. и0 |
иразл) |
(?s = - g - • |
(27.5) |
Тогда из (26.5) при помощи (23.5) и (27.5) получается такое уравнение для функции г (й1, £):
1 |
— ctg3 ft tg- ■ |
r — го — |
|
|
|
|
||
-VPi |
sin'- ft |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
|
1 — ctg2 ft tgftp |
= |
v (r0) (t — t0), |
(28.5) |
||
— гn г — -52 -гэ P |
sin3 ft |
|
||||||
в котором содержится |
важный параметр (Зх: |
|
||||||
|
|
|
|
3QStg ftp |
|
|
(29.5) |
|
|
|
|
|
2 и о т (го) |
|
|
||
|
|
|
|
|
|
|
||
Когда |
t^> t, и, |
|
соответственно, |
г^ > г 0, уравнение |
(28.5) |
|||
записывается в более простом виде: |
|
|
|
|||||
|
г3 |
„ |
1 — Ctg3ftlg 2ftn |
г = |
v (го) t. |
(30.5) |
||
|
2 |
^ |
|
sin2 ft |
||||
|
|
|
|
|
|
122 ГЛ. V. ОКОЛОЗВЕЗДНЫЕ ОБОЛОЧКИ НОВЫХ
Форма сечения оболочки г (в) в данный момент t зави сит только от значения (Зх и угла f>0. Величину рх уместно назвать поэтому «параметром формы», в отличие от [3, яв ляющегося «кинематическим параметром». Значения р вы ше были определены по наблюдаемому значению скорости главной оболочки новой, а рх находится но наблюдениям ее формы. Зная рх, можно оценить массу околозвездной оболочки и скорость потери вещества двойной системой.
У некоторых из новых звезд — DQ Hei, Т Aur — оболочки, наблюдаемые через десятки лет после вспышки, представляются овальными. Кроме того, по снимкам, по лученным с фильтрами, обнаружена неоднородность свече ния оболочки. В ней имеется так называемый экватори
альный |
пояс — область, излучение |
которой в частотах |
||||
линий [N II] и Н повышено по отношению к излучению |
||||||
в линиях [О III]. |
У центральных областей оболочки из |
|||||
лучение |
в |
линиях [О III] |
интенсивнее, чем в линиях |
|||
N [II]. |
Наиболее |
отчетливо |
эти особенности |
выявились |
||
в оболочке |
новой |
DQ Her |
(1934). |
Приняв во |
внимание |
двойственность системы, Э. Р. Мустель и А. А. Боярчук установили [83], что экваториальный пояс расположен в орбитальной плоскости системы. Большая ось овала перпендикулярна к лучу зрения и благодаря этому обо лочка DQ Her особенно пригодна для изучения ее формы. В другом, много обсуждавшемся, случае оболочки новой V 603 Aql, где по спектральным данным можно предполо жить структуру, аналогичную имеющейся у оболочки DQ Her, ось оболочки составляет малый угол с лучом зрения и поэтому вытянутость оболочки на снимках неза метна .
Для объяснения наблюдаемой вытянутости оболочек у некоторых новых была выдвинута гипотеза о решающем влиянии на форму оболочки очень сильного магнитного поля дипольного характера [84]. Предполагается, что это поле заставляет газ оболочки двигаться преимущественно вдоль магнитных силовых линий, т. е. вдоль оси диполя, препятствуя движению в перпендикулярном направле нии. Как нетрудно видеть, для заметного воздействия на движение главной оболочки, кинетическая энергия кото рой порядка 1048 эрг, поле должно обладать сравнимой энергией, а это приводит при размерах системы около 1011 см к средней по занимаемому ею объему напряжен-
|
$ 2. ВЛИЯНИЕ ОКОЛОЗВЕЗДНОЙ ОБОЛОЧКИ |
123 |
носты |
10Gгс. На поверхности же новой звезды напря |
|
женность |
поля должна достигать 108 гс. Присутствие |
столь сильных магнитных полей в тесных двойных систе мах звезд карликов пока ничем не доказано. Сама возмож ность такого воздействия поля на движение оболочки, ко торое приведет к наблюдаемой ее форме, количественно пе продемонстрирована.
С нашей точки зрения можно объяснить, причем более естественно, наблюдаемую вытянутость оболочек новых и структуры оболочек на основе сделанных выводов о ре зультатах взаимодействия главной и околозвездной обо лочек. Очевидно, что в направлении, перпендикулярном к орбитальной плоскости (и близких к нему), газ, состав ляющий главную оболочку, движется беспрепятственно. Вблизи орбитальной плоскости движение газа затормажи вается. Поэтому по прошествии достаточно долгого вре мени после вспышки, когда главной оболочкой будет за хвачена достаточно большая масса газа, оболочка должна быть вытянутой в направлении, перпендикулярном к орби тальной плоскости системы.
Рассмотрим случай, когда большая ось овальной обо лочки приблизительно перпендикулярна к лучу зрения. При условии, что вначале главная оболочка была сфериче ской и ее деформация обусловлена только взаимодействием с главной оболочкой, по наблюдаемому отношению малой оси овала 26 к большой оси 2а определяется величина
параметра |
Учитывая, что при й = |
значение г = 6, |
а величина а = v (г0) t (r^ > r0), из (30.5) |
находим соотно |
|
шение |
|
|
|
Pi ~2— \-Ь — а, |
(31.5) |
которое дает следующую величину рх:
(32.5)
Значение v (г0) определяется по смещению в первое время после вспышки линий поглощения, принадлежащих главной оболочке. Для Новой DQ Her (1934) по наблю дениям 1935 г. она составляла 3,9-107 см/сек. В 1943 г. (при t = fP> = 2,7- 10s сек) по фотографиям оболочки

124 |
ГЛ. V. |
ОКОЛОЗВЕЗДГТЫЕ |
ОБОЛОЧКИ |
НОВЫХ |
|
|
DQ Her |
[85] |
получается |
5 = 0,9, а |
в 1956 |
г. (при |
|
t = <(2) == |
6,8-108 сек) величина (yV J =0,8. Для отноше |
|||||
на) |
|
№ |
|
2,5. Оно совпадает |
со зна |
|
ния —- т- = |
-т— имеем значение |
чением этого же отношения, находимым по формуле (32.5). Величина ^ = 2,3-10-17 см, если принять указанное зна чение v (г0).
Будем считать, что скорость потери массы системой DQ Нет оставалась постоянной до вспышки и оценим при помощи найденного значения массу Ат вещества, при соединившегося к моменту t к главной оболочке. По-види мому, экваториальный пояс оболочки возник вследствие захвата ею вещества. Тогда ширина пояса Н 2Ъ ctg й 0. По фотографиям оболочки [85] получаем значение tg ^^ S !,
а величина т (г0) из (29.5) равна: |
|
||
т (г0) |
3QS |
(33.5) |
|
bt0(3i |
|||
|
|
Количество Ат вещества, присоединившегося к главной
оболочке за время t (t |
£„), составляет |
|
|
|||
Ат |
JL — п |
v М 1 |
Ъ |
(34.5) |
||
а |
||||||
|
ио |
^ s |
ио |
|
Поэтому отношение массы Ат к первоначальной массе главной оболочки равно
Дт |
- L p lV(r0) . t± - . |
(35.5) |
|
т (го) |
|||
|
|
При указанных значениях наблюдаемых величин v (г0), Ыа и Pj получаем, что через 25 лет после вспышки Ат ^ ^ 0,3 т (г0). Присоединившееся вещество распределено не по всей оболочке, а сосредоточено в ее экваториальной области. Поэтому количество газа вдоль радиуса в оболоч ке DQ Her вблизи орбитальной плоскости больше, чем в других направлениях. Различия такого рода неизбежно должны сказываться на наблюдаемом излучении оболочки. По-видимому, существование экваториального пояса, т. е. области, излучающей в линиях [N II], сильнее, чем дру гие части оболочки, связано со специфическими физиче скими условиями в этой области. Плотность, температура и
§ 3. КОРОНАЛЬНЫЕ ЛИНИИ |
125 |
условия возбуждения и ионизации в ней должны быть иными, чем там, где газ не испытал столкновения с диско образной оболочкой [86].
В работе [86] показано, что наиболее вероятной причи ной наблюдаемых особенностей свечения экваториального
пояса является |
большая |
оптическая толщина оболочки |
в направлениях, |
близких к орбитальной плоскости, для |
|
излучения, ионизующего |
атомы 0 + и N+. Концентрация |
атомов 0 ++в области экваториального пояса низка, а ато мов N+ высока, и поэтому там излучение в линиях [N II] преобладает.
В заключение этого раздела заметим, что и при учете газодинамических эффектов закон движения главной оболочки не отличается сколько-нибудь существенно от (13.5) или (25.5) [86]. Однако разогрев газа ударной вол ной, образующейся при столкновении, сильно влияет на свечение оболочки, приводя, в частности, к появлению корональных линий в спектре. Подробно эти вопросы рас смотрены в следующем разделе.
§ 3. Корональные линии в спектрах новых звезд
Появление корональных линий в спектрах многих но вых и повторных новых после вспышки служит одним из главнейших свидетельств присутствия вокруг них еще до вспышки больших масс разреженного газа, т. е. околозвездных оболочек. В этом газе, при взаимодействии его с быстро расширяющейся главной оболочкой, сорванной в момент вспышки, проявляются специфические газоди намические эффекты, которые и приводят к интенсивному свечению в частотах корональных линий. Ниже излагают ся результаты теоретического исследования процессов, приводящих к образованию корональных линий [80], и расчетов их интенсивностей [87], а также определения на этой основе, по наблюдаемым интенсивностям линий, масс околозвездных оболочек.
Плотность в околозвездной оболочке р0 (г) мала по сравнению с плотностью главной оболочки ргл (г) при всех значениях г. Этот вывод следует, в частности, из зна чений параметра р, определяемого формулой (14.5) для конкретных новых. Даже если принять, что масса главной оболочки распределена равномерно с плотностью ргл
12fi |
ГЛ. V. ОКОЛОЗВЕЗДНЫЕ ОБОЛОЧКИ |
НОВЫХ |
|
внутри сферы радиуса г„, то при значениях |
г0 |
5• 1013 см, |
соответствующих радиусу главной оболочки в момент максимального блеска и Р = 10_0 -н10-5, р0/ргл < 1- Так как на самом деле толщина главной оболочки гораздо меньше ее радиуса, то р0<^; рглПоэтому главную оболочку можно рассматривать как сферический поршень, движущийся по газу околозвездной оболочки со скоростью, во много раз превышающей скорость звука. Перед поршнем распро страняется сильная ударная волна, нагревающая газ до температуры в несколько миллионов градусов. Концентра ция атомов в газе мала — около 10е—107 атомов/слл3 и, значит, условия в нем сходны с существующими в солнеч ной короне. Поэтому спектр излучения газа и содержит запрещенные линии очень высокого возбуждения — [Fe X], [Fe XIV] и другие, свойственные спектру солнеч ной короны.
Для вычисления интенсивностей корональных линий, последующего сравнения теории с наблюдениями и опре деления этил! путем физических условий в околозвездной оболочке, потребовалось предварительно решить газоди
намическую задачу о движении |
сферического поршня |
в среде с убывающей плотностью |
[80]. Это решение для |
случая сферически-симметричного распределения газа при водится ниже.
При распространении по газу ударной волны перед поршнем образуется слой сжатого газа, движущегося в ту же сторону, что и поршень. Если главная оболочка (поршень) приобрела скорость только в момент ее срыва под действием взрыва, то в дальнейшем суммарное коли чество движения оболочки и приведенного ею в движение газа не меняется. Как легко показать, газодинамическая задача о движении сферического поршня с постоянным количеством движения по среде, плотность в которой р ~ г-2,— а именно так меняется, согласно (9.5) плотность в околозвездных оболочках,— люжет быть сведена к авто
модельной, если |
пренебречь противодавлением. Пусть |
Р — количество |
движения поршня, приходящееся на |
единицу его площади и р0 (г) — начальное распределение плотности в газе, определяемое (9.5).
Для Р и р0 (г) имеем |
|
Ргг = Р 0Го0 si а; р0 (г)г2= р00г20 = Ь, |
(36.5) |
|
§ |
3. К О РО Н А Л ЬН Ы Е |
Л И Н И И |
127 |
где р00 = |
р0 (г00), |
Р 0 — Р (г00). Величины а и b являются |
||
в данном |
случае единственными параметрами с независи |
|||
мыми размерностями, а тогда, как |
известно |
[88], задача |
о движении газа под действием поршня автомодельна. Решение задачи о распространении ударной волны по
околозвездной оболочке позволит найти, как меняется температура газа со временем. В дальнейшем будем счи тать, что энергия газа за фронтом волны быстро высвечи вается, т. е. уменьшается вследствие излучения, причем закон высвечивания определим условием
- f £ = °- |
(37.5) |
Это условие приближенно выполняется во многих случаях распространения ударной волны [50; 89].
Система уравнений газодинамики в сферически-симмет- ричном случае при условии (37.5) имеет такой вид:
|
w + ” ! |
— |
0’ |
|
<38-5> |
|
д In р |
д In р . |
dv |
, |
2v |
(39.5) |
|
- a r + v ~ d i - + |
5Г + — = |
|||||
|
где v — скорость газа. Она сводится к системе обыкно венных дифференциальных уравнений для функций V (£) и G (£), связанных с v и р следующим образом:
v = J-V(l)-, |
р = р0(г)б(£). |
(40.5) |
Безразмерная переменная £ записывается в такой форме:
ь_ _гМ‘/» _ /_Роо_ |
т- \7« |
(41.5) |
|
a t — I Ро |
t |
||
|
Скорость ударной волны D получается из (41.5) и равна
О = |
- 5 - , |
(42.5) |
где г — координата фронта |
|
|
у Г |
Ра "I Vi |
(43.5) |
|
|
и — значение величины £, соответствующее фронту. Величины v и р зависят, помимо указанных параметров,
128 |
ГЛ. V. ОКОЛОЗВЕЗДНЫЕ ОБОЛОЧКИ НОВЫХ |
также |
и от величины скачка плотности на фронте |
волны G0. Функции V (£) и G (£) находятся сравнительно просто путем решения системы обыкновенных дифферен циальных уравнений [80]. Анализ выражений v и р, ■определяемых согласно (40.5), которые здесь из-за их громоздкости не приводятся, показал, что в результате высвечивания энергии за фронтом волны плотность газа неограниченно возрастает при приближении к поверхно сти поршня [80]. Часть захваченного ударной волной ве щества присоединяется к поршню — главной оболочке, а слой сжатого газа, движущегося перед поршнем, состав-
ляет лишь долю, равную 1 — у |
от всего возмущен |
|
ного газа. |
При значениях G., ^ |
4, характерных для силь |
ной волны, |
эта доля менее 1/4. |
|
Давление газа д(2> непосредственно за фронтом волны связано со скоростью волны в случае, когда газ одноатом
ный, известной зависимостью: |
|
|
||
|
Р™ = -тРоО'1. |
|
(44.5) |
|
Используя (42.5), |
а также учитывая, что величина |
|||
pW связана с температурой Г(2>следующим образом: |
||||
р(2) = |
ЗР»_М"(*), |
|
(45.5) |
|
находим из (44.5) температуру за фронтом волны: |
|
|||
r to |
_3__Р|L |
1 |
(46.5) |
|
|
16 |
рои liGu |
t |
|
|
|
При помощи величины [3, определяемой формулой (14.5),
подставив |
значения |
$ ^ 2 , |
G0 = 4, записываем зависи |
|
мость Т^2) |
от £ в виде |
|
|
|
|
r (s)»2 .10^V -5 -"T - |
(47.5) |
||
В случае |
оболочки |
системы |
RS Орй по (47.5) |
получаем |
7Ч2>^ 3,5-107 ц °К при t — 10° сек и Т<2>^ 4• 10е ц °К
нри t — 8 -106 сек.
Соотношение (47.5) получено при условии (37.5); если закон, определяющий сток энергии за фронтом, другой,
§ 3. КОРОНАЛЬНЫЕ ЛИНИИ |
129 |
то коэффициент в (47.5) может быть несколько иным. За висимость же ТW от времени останется той же самой, по скольку скорость волны D ~ t~x.
Если газ за фронтом волны состоит преимущественно из ионизованного водорода и величины электронной тем пературы Те и температуры тяжелых частиц T t совпадают,
то надо принять ц ^ V2 и значения Г*2) оказываются 1,7-107°К и 2 -10°°К соответственно.
Для того чтобы газ излучал в линиях [Fe XIV], элект ронная температура в нем Тедолжна быть не менее 2-106°К, а для свечения в линиях [Fe X] нужно, чтобы выпол
нялось условие ГГ ^ 1 ,2 -1 0 6 К° [90]. Поэтому следует ожидать, что на первом этапе вспышки новой, когда ско
рость волны и, соответственно, значения T f5велики, наи более интенсивными должны быть линии [Fe XIV]. По
мере уменьшения со временем величины Tf* относитель ная интенсивность линий [Fe X] должна увеличиваться. Вычисления интенсивностей корональных линий [87] под твердили эти ожидания.
Энергия Elt излучаемая в частоте корональной линии [Fe X], определяется выражением
(48.5)
где пре — концентрация атомов железа в газе и Фх (Те) — функция, значения которой приводятся в работе [90]. Интегрирование производится по всему объему излучаю щего газа. Аналогичной формулой определяется и энергия Е е, излучаемая в линии [Fe XIV]. В случае сферической симметрии слоя излучающего газа и в предположении о постоянстве химического состава во всем его объеме имеем
£т,2 (0 |
(Г’ 1) ф 1.2 [Те {п *)] г2 dr, |
(49.5) |
где /у— координата того уровня в оболочке, на котором значение Те становится недостаточным для того, чтобы обеспечить излучение в корональных линиях. При вычис лениях [87] принималось,. что процесс температурной
5 в. Г. Горбацкий