Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Горбацкий, В. Г. Новоподобные и новые звезды

.pdf
Скачиваний:
13
Добавлен:
19.10.2023
Размер:
6.44 Mб
Скачать

120

ГЛ. 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 в. Г. Горбацкий

Соседние файлы в папке книги из ГПНТБ