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

книги из ГПНТБ / Каули, Ч. Теория звездных спектров

.pdf
Скачиваний:
16
Добавлен:
21.10.2023
Размер:
9.13 Mб
Скачать

180 г л а в а е

Чтобы получить численные значения, нужно от операторной за­ писи перейти к матричной:

1

aß)).

(6.3.9)

a'ß' І (со — Wo) + Ф

ста'

ßß'

Выражение в двойных скобках в (6.3.9) является матричным элементом обратного оператора [соо как оператор равен

(#о — Но)/ h], так что сначала необходимо найти матричные эле­ менты /(со — соо)+Ф, а затем обратную матрицу. Эта операция может оказаться очень трудоемкой, так как в общем случае порядок рассматриваемых матриц равен произведению чисел элементарных состояний на верхнем и нижнем уровнях. При рас­ смотрении линий водорода пренебрегают спиновым вырожде­ нием, но для линии Н6 все же требуется обратить матрицу по­ рядка 22-62 = 144. Очевидно, что для более высоких членов бальмеровской серии нужно прибегнуть к аппроксимациям.

Приближенное рассмотрение обращения матриц высокого порядка было использовано Мозером [115] в его работе по во­ дороду. В одном случае обращение матрицы тривиально — это случай диагональной матрицы. Тогда элементы обратной ма­ трицы являются просто обратными матричными элементами. Можно воспользоваться этим фактом для обращения матриц, у которых недиагональные элементы малы. Пусть, следуя Мо­ зеру, В — такая матрица, а А — диагональная матрица. Тогда имеем

7 T T = T - i ßi + i ßi

e T -

<6-3'10>

В качестве А Мозер взял матрицу /Асо +

Ф', где Ф' — матрица,

составленная из диагональных элементов

Ф. Тогда для обраще­

ния

нужно взять В — Ф — Ф'.

Конечно,

гораздо проще иметь

дело

с правой частью (6.3.10),

чем с левой. Мозер нашел, что

в естественном представлении эффекта Штарка, т. е. в пред­ ставлении параболических координат, не только гамильтонианы

Но и Но диагональны, но и Ф «заметно» диагонально, так что можно использовать разложение (6.3.10).

Энергетические уровни Еи и Е1 могут расщепляться квази­ статическими полями, вызванными ионами. Следовательно, ус­ реднение по ансамблю должно проводиться по распределению электрического поля. Таким образом, окончательное выражение для контура линии будет

/ (Асо) сс

 

1

 

І (ü) —(Qq) Ф aß)M ^)-

о

cta'

Pß'

(6.3.11)

КВАНТОВОМЕХАНИЧЕСКОЕ РАССМОТРЕНИЕ УШИРЕНИЯ ДАВЛЕНИЕМ 18)

Здесь функция Н описывает плотность распределения вероятно­ стей электрического поля (разд. 5.7). Полезно переписать это уравнение в форме, которая показывает усредненное со0 = = (Еи — Е1)Ш в функции электрического поля &. Если S Ф О, то соо не является истинным центром линии. Обозначим через со0о истинный центр линии, а через io0 — со0о = А а/— смещение, вы­ званное (%. Тогда можно усреднять по функции распределения Я (Асо') Для Асі/, и (6.3.11) примет вид

 

і (Ди — Ди') + Ф

где

(6.3.12)

 

Асо = (о — со00.

(6.3.13)

6.4. ОЦЕНКА Ф ДЛЯ ВОДОРОДОПОДОБНЫХ СПЕКТРОВ

Прежде чем оценивать Ф, сделаем несколько замечаний от­ носительно физической картины, принимаемой для усреднения

вуравнении (6.3.4). Для примера рассмотрим усреднение по ан­ самблю совсем простой «системы», которую нужно модифициро­ вать, если рассматривается действительная физическая ситуа­ ция. Однако принимаемая картина исключительно близка к той, на которой часто основаны практические вычисления.

Пусть система состоит из одного атома, подверженного воз­ мущениям одной частицей, и только одно соударение происходит

винтервале времени t = —оо, t — + о о . Такое соударение опи­ сывается на рис. 6.4.1 прямой линией — классической траекто­ рией возмущающей частицы. Ансамбль состоит из всех возмож­ ных систем такого вида.

Будем считать, что нет предпочтительного направления в про­ странстве или предпочтительной ориентации траектории возму­ щающей частицы. Стоит только рассмотреть все возможные столкновения, т. е. выполнить усреднение по ансамблю, как ста­ нет очевидно, что результат должен обладать сферической сим­ метрией. Можно наглядно показать это, детально проводя некоторые этапы усреднения по ансамблю (см., например, [31]). Здесь же мы приведем аргументацию, основанную на физиче­ ских предпосылках, согласно которой все нечетные функции для

возмущающей частицы и атомные координаты обратятся в нуль при усреднении по ансамблю, а все четные функции должны быть сферически симметричными. Таким образом, если х н у — координаты атома, то при усреднении по ансамблю матричные элементы х, у, ху и т. п. обратятся в нуль, а выражения вида х2

182 ГЛАВА 6

или у2, которые получились при данной ориентации, показанной на рис. 6.4.1, станут равными г2/3. Аналогично нечетные функции ( R) v = Р и (Ѵ)х = V должны обратиться в нуль при усреднении

по ансамблю, так как знаки плюс и минус у р и ѵ будут встре­ чаться в ансамбле с одинаковой частотой. Хотя р и ѵ будут иногда появляться в наших формулах, но только в смысле зна­ чений I р| и IVI, которые не исчезают при усреднении.

Иногда мы будем обращаться к среднему по сфере, чтобы описать усреднение по ансамблю, при котором х2 становится г2/3 и т. д. Несколько результатов усреднения по сфере пред­ ставлено в табл. 6.4.1.

 

 

В е л и ч и н а

X, у . Z, х у , х у 2 и т. д.

X 2,

у 2,

Z 2

X4,

у 4 ,

Z 4

х 2у г , x 2z 2 и т. д.

к р и р ' и 2

Таблица 6.4.1

З а м е н я е т с я п ри с ф е р и ч е с к о м у с р е д н е н и и на

0

г213

г4! 5

г4/ 1 5

( 2 / + 1 ) _ І 2 \{aJM \r\a'J'M')\2

ММ '

1(Рк1 1Р')|*

( 2 / + I ) - 1 2 l ( a / A f 1 г2 1 аЧ'М')\2

 

ММ '

Из табл. 6.4.1 следует, что матричные элементы г нужно ин­ терпретировать как усреднение по 2/ + 1-кратно вырожденному состоянию ß — aJM. Этот результат включен в табл. 6.4.1 для справки. Он наиболее полезен при рассмотрении отдельных ли­ ний (ср. разд. 6.9).

Усреднение также проводится по возможным значениям ско­ рости возмущающей частицы и по всем значениям прицельных параметров. При этом используются функция распределения скоростей (распределение Максвелла) и прицельные пара­ метры, которые будут введены в дальнейшем.

Оператор Ф определяется соотношением (6.3.4). Если пере­ нести начало отсчета времени с t — 0 на t — т, то экспонента исчезнет*). Поскольку усреднение будет проводиться по множе­ ству столкновений, в котором времена фактических столкновений не коррелированы, то можно сосредоточить внимание на интер­ вале Ат, в течение которого совершается столкновение. Далее, поскольку возмущения быстро спадают к нулю, можно распро­ странить интегрирование в выражении для U от — оо до -J- оо.

*) См. разд. III. 5. При эволюции U(оо, —оо) столкновение считается «сосредоточенным» в t = 0. Это означает, что малыми фазами порядка ДЯр/ufi пренебрегают.

КВЛНТОВОМЕХАНИЧЕСКОЕ РАССМОТРЕНИЕ УШИРЕНИЯ ДАВЛЕНИЕМ 183

Пусть Q — полное поперечное сечение столкновений, которые совершаются через интервалы времени Ат и удовлетворяют ус­ ловиям ударного приближения. Тогда для частоты таких стол­ кновений с возмущающими частицами, имеющими скорости ме­ жду V и V + dv, можно записать

 

 

 

1/Ат = Nf(v)Qvdv,

(6.4.1)

где

N — полное

число возмущающих

частиц,

f ( v ) — максвел­

ловское распределение скоростей.

 

 

 

Таким образом, получаем

 

 

 

оо

 

 

 

 

Ф

= |

f { v ) d v \

^ ^ - N Q ü { U u{oo,

— оо) [/'(оо, — оо) — 1}ср.

 

о

 

 

 

(6.4.2)

 

 

 

 

 

Отметим, что поперечное сечение Q, которое определено соотно­

шением

(6.4.1),

сокращается. Мы ввели его,

чтобы наглядно

показать усреднение по прицельному параметру р. Пределы ин­ тегрирования по р умышленно опущены (см. ниже). Индекс «ср4» сохраняется в (6.4.2), поскольку усреднение по ансамблю должно также распространяться на ориентацию атома и воз­ мущающей частицы.

Оператор эволюции во времени в картине взаимодействия

разлагается по возмущениям (см.

приложение III

и приведен­

ные в нем ссылки) следующим образом:

 

и (Ат, 0) = 1 — (і/А) JДт V' (/,) dtx

 

 

о

 

 

Дт

tt

 

~ { \ т J V (f,) dtl j V' (t2) dt2+

(6.4.3)

о

0

 

причем оператор V'(t) связан с

гамильтонианом

возмущения

V(/) соотношением

 

 

V' (t) — exp (iH0t/h) V (t) exp (— ІН4ІЩ.

(6.4.4)

Одним из упрощений при рассмотрений водородоподобных спек­ тров является возможность пренебречь экспонентами в (6.4.4) вследствие почти полного вырождения верхнего и нижнего со­ стояний. Это верно только в нашем конкретном приближении, когда верхние и нижние состояния можно считать независимыми. Отметим, что все операторы V в (6.2.9) обозначены индексами и п І я что они действуют только на верхние или нижние состояния

184 ГЛАВА 8

соответственно. Таким образом, матричные элементы ІІи бе­ рутся только по отношению к верхним состояниям и т. д. Если мы выберем обычное диагональное представление для энергии (ср. разд. III. 3), то в матричных элементах все экспоненты исчезнут *).

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

делы интегрирования в

(6.4.3)

до — оо

и +

оо.

 

 

*

Мы допускаем раздельное

разложение

для Uu и

U1 вида

(6.4.3). Поскольку гамильтониан возмущений

V действителен,

U1 берется с заменой знаков у нечетных

членов

(которые со­

держат і) в (6.4.3).

записать (ср.

[86,

§

67;

92,

р. 124])

Возмущения

можно

в форме

 

 

 

 

 

 

 

 

F =

f V F +

w f 3 (r ‘ R)2“

r2jR2]+

••••

(6Л5)

Первое слагаемое представляет собой взаимодействие атомного дипольного момента с возмущением, а второе — квадрупольное взаимодействие, причем второе слагаемое приводится лишь для дальнейших ссылок, а здесь мы ис-

 

______ пользуем

только

дипольный член

 

 

 

разложения.

 

 

 

 

 

 

/>*/>¥

Переменные и координатная си­

frvtzyy

 

стема,

использованные

для описа­

 

 

ния

столкновения,

показаны

на

 

 

 

 

 

 

рис. 6.4.1.

R — вектор,

направлен­

 

 

 

ный из центра излучателя к возму­

Рис. 6.4.1.

Переменные,

ис­

щающей

частице,

г — координата

оптического электрона,

р— вектор,

пользуемые

при рассмотрении

направленный в точку

максималь­

столкновения, х и у — единич­

ные векторы.

 

ного

приближения

возмущающей

 

 

 

частицы, а V — ее скорость.

за­

Чтобы оценить выражение в фигурных скобках

(6.4.2),

писываются отдельно

разложения

для

Uu и

Ü1 вида (6.4.3)

и

результаты перемножаются. Все дипольные члены первого по­ рядка обращаются в нуль, когда берется усреднение матричных элементов по всем углам **). Поэтому первыми не обращающи-

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

**) Т. е. по всем ориентациям R относительно і\

КВАНТОВОМЕХАНЙЧЕСКОЕ РАССМОТРЕНИЕ УШИр ЕМИЯ ДАВЛЕНИЕМ 185

мися в нуль в результате усреднения членами станут члены вто­ рого порядка и один «перекрестный» член:

£/“(/' — 1 = — Ъ2

I'00

Vu(t1)dtl

Jt x V“(t2)dt2 +

 

 

+

^T

ooJ

 

 

oo

J

W(*,)<«,-

 

 

Й2

Joo

V l

( f \ ) d h

t xI

(t2) dt2+

. . . .

(6.4.6)

Мы опять использовали индексы «и» и «/» в гамильтониане воз­ мущения К, которые вводились отдельно для и U1, чтобы по­ казать, что матричные элементы этих операторов могут браться соответственно только по отношению к верхнему и нижнему со­ стояниям.

Рассмотрим первый член в правой части уравнения (6.4.6).

Мы пишем просто «интеграл» для одного из двойных интегра­ лов, поскольку его коэффициент обратится в нуль при усредне­ нии по ансамблю.

Двойные интегралы можно упростить следующим путем. Ин- . тегралы берутся по полуплоскости, заштрихованной на рис. 6.4.2. Однако подынтегральные выражения инвариантны по отноше­

нию к подстановке

t\ -> 1\,

Следовательно,

можно

расширить пределы

интегрирования

на всю плоскость

(ti,t2),

умножив на У2. Двойные интегралы сводятся к квадратам одно­ кратных интегралов. Третье слагаемое в (6.4.7) содержит инте­ грал нечетной функции от — оо до + оо и потому обращается в нуль. Первый интеграл вычисляется при помощи формул

186 ГЛАВА б

разд. 5.5, и в результате имеем

2е4

3 Ь2ѵ2р2 Г2,

(6.4.8)

Ъ2ѵ2р2 У2

поскольку у2 = г2/3 при усреднении по сфере.

аналогично.

Другие интегралы в (6.4.6)

можно вычислить

Перекрестный член первого порядка после усреднения по сфере равен

4е4

4е4

 

Ь2ѵ2р2 У У

* ЗЬ2ѵ2р2

( 6 . 4 . 9 )

 

Таким образом, выражение для Ф можно записать в следующем виде:

 

00

 

ф =

f(v) ~ $

[(г“)2 — 2r“ • r' + (r')2]. (6.4.10)

 

о

 

Интегралы по p в (6.4.10) логарифмически расходятся по обоим пределам. Эта трудность имеет место всегда, когда необ­ ходимо проводить усреднение по прицельному параметру, и с ней встречаются исследователи звездной динамики. В приведенном виде теория ошибочна и для больших, и для малых р, так что имеет смысл ограничить пределы интегрирования значениями р, начиная с которых теория больше не верна.

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

Ртах « Pd = {kTISnNe2)'!*.

(6 .4 . 1 1 )

Второе соображение заключается в том, что длительность столк­ новения должна быть такой, чтобы оставалась применимой кон­ цепция уширения вследствие соударений в противоположность квазистатическому приближению (ср. разд. 6.7). Таким образом, смещение частоты Дсо линии должно быть меньше частоты со­ ударения, откуда следует

P m a x ^ W A ö .

( 6 . 4 . 1 2 )

Согласно [63], нужно брать меньшее из значений ртах, даваемых

(6.4.11) и (6.4.12).

б. Для все меньших И Меньших прицельных параметров на­ рушается предположение о том,'что среднее соударение является слабым или что теория возмущений применима. Пусть в урав­

нении (6.3.3) t2 oo, ti

— ос, и мы требуем,

чтобы А f a 1.

Вместо Т можно использовать U, и если мы ограничимся первым

порядком в разложении

(6.4.3)

для U,

то будем

иметь

I (а I £/ — 11а) J

1

е2 І (« I

г I а) I Ршіп

(6.4.13)

Ѣ

2

 

 

 

 

Pmin

V

 

 

 

 

 

 

КВАНТОВОМЕХАНИЧЕСКОЕ РАССМОТРЕНИЕ УШИРЕНИЯ ДАВЛЕНИЕМ [87

При

(а [ г I а) ~ п2аа= п2Н2/те2

имеем

Pmiti ^ n2h/mv.

(6.4.14)

Сильное столкновение (р < pmin) можно приближенно учесть, вводя простое допущение (см. заключительные замечания в разд. 5.5) о том, что

у (сильное) = — 2Nnp2minv.

(6.4.15)

Знак минус приходится брать для совместимости с уравнением

(6.4.10). С учетом (6.4.14) находим

у (сильное) = 2 М -^ -п 4 Jсо

f(v)-^~.

(6.4.16)

о

 

 

Эту величину можно ввести как поправку, учитывающую силь­ ные столкновения, к результату интегрирования (6.4.10) от рт!П

ДО Р та х -

Тогда имеем

<D = - f g ^ j V { f(v) ^ { s + [(r“)2- 2 r “ - r 4 ( r ;)2] l n - ^ } ,

 

(6.4.17)

где S — слагаемое, учитывающее сильные столкновения:

 

S = ЗН4п4/4е4т2.

(6.4.18)

При этом нужно подставлять главное квантовое число п верх­ него состояния, поскольку расщепление из-за эффекта Штарка растет пропорционально п2.

Чтобы оценить оператор Ф для ионизованного гелия и дру­ гих водородоподобных атомов, нужно анализировать орбиты заряженных частиц и учитывать отклонение траектории от пря­ молинейной. Мы не будем подробно останавливаться на этом анализе, так как принципы его те же, что и при расчете водо­ рода, а алгебраические выкладки хорошо объясняются в [63].

Выводом формулы (6.4.17) мы закончили теоретический рас­ чет штарковского контура для водорода. Чтобы получить чис­ ленный результат, фактически нужно произвести вычисления матричных элементов Ф, которые включают отыскание матрич­ ных элементов г и г2. Для системы с одним электроном опубли­ кованы аналитические выражения (см., например, [12]). Окон­ чательные матричные элементы объединяются затем с «триви­ альными» матричными элементами і (со — соо), и полная матрица

188

ГЛАВА 6

должна обращаться методом, рассмотренным в разд. 6.3. Нако­ нец, контур линии вычисляется по уравнению (6.3.12). Ряд таких контуров опубликован (см. разд. 6.6). Детальные вычисления приводят к необходимости получения приближенной формулы для крыльев линии, которая рассматривается в следующем раз­ деле.

6.5. МОДИФИЦИРОВАННОЕ УДАРНОЕ ПРИБЛИЖЕНИЕ ФОРМУЛЫ ДЛЯ КРЫЛЬЕВ ЛИНИЙ

Пусть ос (ДА)— зависящий от длины волны коэффициент по­ глощения, обусловленный действием электронов и ионов, а ссДДА)— аналогичный коэффициент, обусловленный действием одних только ионов (ср. разд. 6.7). В ударном приближении отношение а а)/а; (ДА,) является простой функцией длины волны, температуры Т и электронной концентрации N. Тогда можно записать [68]

а(ДА) = ссг(ДА) {і + R{N, Г) ] / äÄ}.

(6.5.1)

Функция R(N, Т) затабулирована в [62, 68] для нескольких во­ дородных линий и двух линий Hell.

Формула (6.5.1) содержит два члена. Первый пропорциона­ лен (ДА)~5А, т. е. имеет зависимость квазистатического уширения, вызванного ионами, а второй имеет зависимость (ДА)'2, характерную для уширения вследствие столкновений. Рассмо­ трим более подробно второй член.

Из общего рассмотрения гл. 5 следует, что в далекой обла­ сти крыльев электроны будут уширять квазистатически. Дей­ ствительно, в разд. 5.6 мы оценили расстояние Дсоь в крыльях линий, за которым целесообразно применение квазистатической теории.

В разд. 5.6 Дсог, было выражено через постоянную взаимо­ действия ^ „ и параметр цп{Ці = я). ^2 можно найти приравни­ ванием изменения энергии АЕ атома в электрическом поле & величине йДсо, где

Асо = 2я'ё’г/г2.

(6.5.2)

Выражение для ДЕ выведено во многих книгах по спектроско­ пии и атомной физике, довольно полное рассмотрение его приво­ дится в [12]. Если электрическое поле записать как е/г2, то

Aa= w n‘ - k '

<6-5-3>

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

КВАН ТО ВО М ЕХА Н И Ч ЕС К О Е РАССМ ОТРЕН ИЕ У Ш И РЕН И Я Д А ВЛ Е Н И Е М 189

Следуя Унзольду [158, § 82], определим среднее значение nh по­ средством соотношения

«* = 2 % /* /£ /* ,

(6.5.4)

где /ft — интенсивность k-й штарковской компоненты. Таким об­ разом, можно записать среднее

2л^2 = ^кпк/2/п.

(6.5.5)

Подстановка в (5.6.4) дает критическую границу для круговой частоты:

Лш0 = 2v2mr\lj3qlhnk.

(6.5.6)

Перепишем (6.5.6) таким образом, чтобы привести его в

соот­

ветствие с аналогичным выражением, данным в [63]. Установле­ но, что для оценки (6.5.4) с точностью до множителя 2, а обычно

даже более

высокой,

можно брать ііь — п2/3, где

п — главное

квантовое число верхнего

состояния. Учитывая это и полагая

V2 = 3kT/m,

получим

для

критической границы,

выраженной

в длинах волн,

 

 

(6.5.7)

 

АХь = (6цЦп2)ХаІгТ/2лсПп2.

Это выражение можно согласовать с аналогичным выражением,

полученным Гримом, если

выбрать тіо таким, чтобы величина

в круглых скобках (6.5.7)

стала равной единице. К сожалению,

нельзя оценить (6.5.7) точнее чем до порядка величины. Мы не можем точно сказать, нужно ли выбирать г)о таким, чтобы вы­ ражение в круглых скобках равнялось единице, или я, или, мо­ жет быть, 1/я! Экспериментальные работы заставляют считать, что в некоторых интересных для астрофизики случаях АХь лежит достаточно близко к ядру линии, так что при низкой электрон­ ной концентрации можно считать электронное и ионное уширения линии полностью квазистатическими. Читателя, желающего подробнее познакомиться с этим вопросом, отсылаем к статье [46] и приведенной в ней библиографии. См. также [64].

Кроме границы АХъ или Аюь, нужно также рассматривать границу АХр, которая появляется из-за нарушения столкновительного приближения вследствие плазменных взаимодействий при значениях р, превышающих дебаевский радиус. Важно под­ черкнуть, что многократные соударения (рассматриваемые так, как если бы они происходили по одному во времени) при боль­ ших р вызывают большее уширение, чем столкновение с бли­ жайшей частицей. Это происходит потому, что для линейного эффекта Штарка сдвиг фазы г) уменьшается как р_І, в то время как число возмущающих частиц растет как р3. Таким образом, дальние столкновения более эффективны, чем ближние, и если бы не плазменные взаимодействия, то имелась бы расходимость.

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