
книги из ГПНТБ / Каули, Ч. Теория звездных спектров
.pdf180 г л а в а е
Чтобы получить численные значения, нужно от операторной за писи перейти к матричной:
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) |
Мы опять использовали индексы «и» и «/» в гамильтониане воз мущения К, которые вводились отдельно для Vй и 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. Таким образом, дальние столкновения более эффективны, чем ближние, и если бы не плазменные взаимодействия, то имелась бы расходимость.