Ви є тут

Излучательная динамика атомных систем

Автор: 
Безуглов Николай Николаевич
Тип роботи: 
докторская
Рік: 
1999
Кількість сторінок: 
232
Артикул:
1000259850
129 грн
Додати в кошик

Вміст

Оглавление
Оглавление 4
Введение 5
1 Квантовое описание сил самодействия излучением в атоме 25
1.1 Принцип наименьшего действия при исследовании спонтанного излучения 26
1.1.1 Фсйнмаиовский формализм ....................................... 27
1.1.2 Разложение радиационного самодействия по степеням с-1.......... 30
1.1.3 Квантовые члены................1............................... 35
1.1.4 Краткое резюме................................................. 40
1.2 Квазиклассические правила квантования радиационной ширины............. 41
1.2.1 Принцип соответствия для спонтанного излучения................. 44
1.2.2 Вариационный принцип........................................... 46
1.2.3 Перемещение переменных действия в комплексную плоскость ... 47
1.3 Особенности рассмотрения 5-состояний ................................. 50
1.3.1 Процедура регуляризации расходящихся параметров................ 52
2 "Излучательная" модель квантования осциллятора с диссипацией 59
2.1 Эффективный комплексный гамильтониан осциллятора ..................... 59
2.2 Кваптовьге члены, стабилизирующие основное состояние.................. 61
2.2.1 Обобщение гамильтониана (2.2) для нелинейного осциллятора .. 61
2.2.2 Динамика диссипативного нелинейного осциллятора в представлении интегралов по траекториям......................................... 64
2.2.3 Сдвиг уровней, обусловленный диссипацией энергии............... 65
2.2.4 Квантовые поправки ............................................ 68
2.3 Временная динамика и комплексный спектр энергии осциллятора .... 70
3 Радиационные времена жизни атомов и ионов 75
3.1 Эмпирические закономерности распределения радиационных времен жизни по возбужденным состояниям ............................................ 76
3.2 Обзор различных методов вычисления радиационных времен жизни тя . 77
3.3 Одноканальная теория естественной ширины Ап .......................... 83
3.3.1 Обозначения и допущения........................................ 83
3.3.2 Главные члены разложения Ая по энергии е....................... 87
3.3.3 Иизколежащие возбужденные состояния............................ 90
3.3.4 Критерии выполнимости правила Бете и их связь с подавлением
радиационных процессов......................................... 91
2
3.3.5 Квантовые поправки к кваэиклассическим формулам........... 98
3.4 Радиационные ширины 5-состояний.........................................100
3.4.1 Правила подобия.......................................... 100
3.4.2 Вариационная процедура для нахождения расходящихся параметров 102
3.5 Экстраполяционная схема для продолжения радиационных времен жизни
с нижних возбужденных состояний до верхних, включая континуум ... 103
3.5.1 Одноканальное приближение.................................104
3.5.2 Пример: времена жизни всех возбужденных состояний атомов щелочных металлов........................................................108
3.5.3 Пример парциальных сечений фоторекомбинации для медленных
электронов.................................................108
4 Квазиклассическое представление дипольных матричных элементов 113
4.1 Квазиклас.сичсские матричные элементы..................................114
4.1.1 Обобщенные правила соответствия...........................115
4.1.2 Дипольный матричный элемент для длинных переходов.........117
4.2 Квазиклассические сечения фотоионизации в области куперовского минимума ....................................................................121
4.3 Спектр фоторекомбинации на примере модельного потенциала Зоммер-
фельда ............................................................... 123
5 Метод обобщенных волновых уравнений в задачах радиационной кинетики газовых сред 127
5.1 Интегральное уравнение пленения излучения..............................128
5.2 Радиационная кинетика как разновидность задачи о квантовом бильярде 130
5.3 Обоснование концепции ассоциированной квазичастицы в рамках формализма контипуального интегрирования........................................135
5.4 Эталонные уравнения для нахождения решений в окрестности особых
поверхностей...........................................................139
5.5 Резонансные условия и разрешенные конфигурации волновых фронтов . 144
5.6 Рассеяние квазичастицы на границе полубесконечной среды............... 147
5.6.1 Свойства скачка фазы при нормальном падении ....................148
5.6.2 Особенности скачка фазы при наклонпом падении . ...............150
5.7 Практические правила квантования в многомерных задачах переноса излучения .............................................................. 151
5.7.1 Параллелепипед с длинами сторон Нх, Ну, Н2................153
5.7.2 Конечный цилиндр радиуса Я и высоты Н ..........................153
5.7.3 Шар радиуса К.............................................155
5.7.4 Призма-иодобные ^геометрии...............................,157
5.8 Случай эллиптических геометрий.........................................161
5.8.1 Эллиптический цилиндр.....................................162
5.8.2 Вытянутая эллипсоидазьиая кювета : = #< ; Яг = 72>..........167
5.8.3 Сплющенная эллипсоидальная кювета: Я± = /£> ; Кг — Л< .... 169
5.9 Результаты расчетов эффективных констант радиационного раепада . . 169
5.9.1 "Одномерные" геометрии....................................170
5.9.2 Многомерные геометрии.................................... 173
3
5.10 Критерий применимости приближения изолированных особых поверхностей ....................................................................178
6 Стохастическая динамика Ридберговского электрона в переменном электрическом поле квазимолекулы 185
6.1 Внутреннее дипольное поле квазимолекулы в модели Думана-Шматова-Михайлова-Янева (ДШМЯ) .................................................187
6.1.1 Резонансы в сталкивающейся паре ридберговский атом-атом в основном состоянии................................................... 187
6.1.2 Параметр нелинейности модели ДШМЯ и эффекты насыщения . 189
6.1.3 Обсуждение результатов расчета констант скоростей, полученных
в рамках модели ДШМЯ..........................................192
6.2 Диффузия Арнольда ридберговского электрона в однократном соударении 194
6.2.1 Уравнение диффузии РЭ по кулоновскому спектру энергий .... 195
6.2.2 Приближение по Вейскопфу для сечения ударной ионизации
в диффузионной модели........................................ 197
Заключение 201
А Метод геометрического квантования 206
А.1 Метод стационарной фазы в континуальном интеграле....................206
А.2 Многомерные правила квантования для Лаграижевых поверхностей . . . 207 А.З Свойства условно-периодического движения в переменных у гол-действие 209
A.4 Вариационный принцип Персиваля.......................................210
В Математические аспекты проблемы силы трения излучением 212
B.1 Пространственно-временное представление хронологической свертки свободного электромагнитного поля..........................................212
В.2 Разложение сверток в ряд по степеням параметра запаздывания..........214
B.З Процедура регуляризации интеграла (1.27) ............................215
С Определение фазовых факторов для обобщенных волновых уравнений217
C.1 Фаза отражения квазичастицы от плоской потенциальной стенки .... 217 С.2 Среднее значение скачка фазы для эллиптических поверхностей.......220
Библиография 224
4
Введение
Исследования радиационной кинетики газовых и плазменных сред составляют обширный раздел физики, имеющий свою собственную историю развития и накопивший большое число ставших теперь классических методов решения возникающих задач. Быстрое развитие лазерной техники открывает в последние десятилетия новые возможности для экспериментов с оптически возбужденными средами. Типичный пример в этом плане представляет фотоплазма, в которой, как это было показано сравнительно недавно, можно получать рекордно большую (приближающуюся к 100%) степень ионизацию вещества [1, 2, 3]. При этом реализуется целый ряд радиационных [4, о) и столкновительных (6) процессов, ведущих к перераспределению начальной энергии оптического возбуждения по различным каналам, конкуренция между которыми в конечном итоге и определяет специфические свойства фотоплазмы [3, 7]. Широкое внедрение лазерных источников возбуждения в практику современного эксперимента создало также техническую базу для начала комплексного изучения свойств ридбергов-ских состояний атомов, ионов и молекул. Интерес к последним во многом стимулирован такими практическими приложениями, как современные энергосберегающие методы разделения изотопов, создание эффективных методов непосредственного преобразования световой энергии в электрическую, детектирование ифракрасного и микроволнового излучения. Однако в силу целого ряда причин [8, 9] существовавшие до последнего времени теоретические и экспериментальные методы исследования атомных систем были не в состоянии дать надежных сведений о такой фундаментальной характеристике высоковозбужденных состояний как радиационное время
5
жизни.
Таким образом достижения современного эксперимента значительно расширили круг изучаемых явлений, связанных с процессами испускания, переноса и поглощения лучистой энергии ансамблем атомов. Соответственно возросли требования к теории как в плане углубления фундаментальных знаний о новых проявлениях физики взаимодействия излучения с веществом, так и в плане создания базовых методик для эффективного решения задач радиационной кинетики газовых сред. Последнее во многом и определило главную цель исследований, составивших основу диссертации: развитие нетрадиционных представлений о внутренних механизмах, определяющих динамику излучающих атомных систем, с последующим построением на их основе новых математических приемов для расчета констант скоростей различных радиационных процессов.
Материал, представленный в диссертации, посвящен преимущественно решению трех классов задач. Во первых, это описание временной эволюции различных физических процессов (главы 1, 2), сопутствующих спонтанному излучению, и определение вероятностей порождающих его оптических переходов для изолированных атомных систем (главы 3, 4). По мере роста оптической плотности газовых сред выход излучения начинает осуществляться в условиях обмена световой энергией между отдельными атомами. Это определило второй класс проблем, связанный с исследованиями излу-чательной релаксации резонансно возбужденных сред в условиях пленения излучения глава 5). Наконец третий класс задач обусловлен практической важностью изучения элементарных столкновительных реакций типа хемоионизации (глава 6), эффективность протекания которых определяется внутренними переменными электрическими полями в квазимолекуле, составленной из партнеров по соударению. Остановимся более подробно на месте рассматриваемых в диссертации тематик в контексте основных этапов развития современной физики излучательных процессов.
Основная трудность теоретического осмысления экспериментальных
6
данных по радиационным временам жизни тд связана со сложной структурой естественной ширины Ая возбужденных состояний:
Тд = Ад = £д < £д (1)
<7
Особенно это проявляется для ридберговских уровней, порождающих большое число индивидуальных оптических переходов с парциальными вероятностями АЯ)д. Между тем в классической электродинамике эффекты радиационного затухания сравнительно легко учитываются с помощью концепции сил трения излучением Лоренца [10]. Реакция системы на потерю энергии описывается силами радиационного самодействия ~Ра [11]-(16], которые действуют на а-частицу заряда га, и определяются полным диполь-ным моментом ~б(г) атомной системы
1 N
(г) = Ь' га7^{1) ’ ^ = (2)
а=1
Скорость потери энергии (Ш/(И пропорциональна модулю Т^(£) второй производной по времени от /?(£), а уносимая светом энергия \¥(Т) за характерный период времени Т дается усреднением <Ш/(И по классической траектории движения фс/(£) = {®ь системы N частиц:
в? 2
й 2
= —
<И Зс3
, 1Г(Т) = I (3)
Заметим, что радиационное самодействие традиционно определяется через силы Абрагама—Лоренца (2). Последние получаются при подстановке в уравнение Ньютона для заряженных частиц собственных магнитных и электрических напряженностей нолей, выраженных через динамические переменные вещества. Обобщая это наблюдение, естественно радиационное самодействие связать с любой процедурой исключения полевых переменных из системы «вещество плюс электромагнитное поле». В дальнейшем мы будем придерживаться подобного широкого толкования для сил реакции на излучение, имея в вид}' любую классическую или квантовомеханическую ситуацию, когда движение системы излучающих зарядов удается
описать с помощью собственных переменных. Подобным же образом иод термином «излучающие» для частиц будет пониматься факт их взаимодействия (обмена энергией) с электромагнитным полем (в том числе и с физическим вакуумом). Дело в том, что на вопрос, какую часть поля следует трактовать как «излучение» [16], можно дать несколько нетривиальных ответов [13, 16]. Это обусловлено существованием широкого класса решений однородного уравнения Максвелла, необходимостью их доопределения с помощью тех или иных граничных условий [17, 13, 14].
Говоря о силах торможения излучением, надо отдавать себе отчет о степени огрубления реально решаемых задач. Строгое описание движения заряженной частицы в физическом вакууме подразумевает введение в теорию бесконечномерной системы, поскольку число полевых переменных нео-граничено. За те упрощения, которые предоставляет исследователю работа с локальными силами (2), следует "платить". В теории самодействия это проявляется в возникновении целого ряда парадоксов и противоречий типа самоускорения зарядов [10, 11,12,15], преодоление которых иногда требует даже отказа от традиционных физических представлений [13|. Сразу отметим, что сказанное в полной мере относится и к концепции самодействия применительно к квантовым системам. Наличие большого резервуара оттока лучистой энергии, который представляет из себя физический вакуум, приводит к расходимости некоторых параметров в конечных формулах теории. В результате приходиться привлекать тонкие физические [17, 20] и математические приемы [19, 21] для устранения возникающих сингулярностей.
Силы Лоренца (2) локальны во времени, а уносимая фотонами энергия IV в соотношении (3) определяется усреднением работы этой силы вдоль классической траектории движения. Это разительно отличается от ситуации с квантовомеханической вероятностью радиационного распада Ад в (1), представляющую сложное образование из индивидуальных оптических переходов на все нижележащие уровни. Очевидно, что радиационную
8
ширину Ад в таком виде невозможно, в отличие от положения энергии исходного уровня бд, непосредственно связать с характеристиками излучающего состояния д. В связи с этим дискуссия по отысканию аналогов реакции на излучение в квантовой теории неоднократно поднималась в отечественной и зарубежной научной литературе. Однако в большинстве работ (см., например, [11, 12]) проблематика взаимодействия зарядов с физическим вакуумом в конечном итоге ограничивалась исследованием узкого класса динамических характеристик движения вещества (импульсов и координат), усредненных по квантовым амплитудам вероятности.
Ричард Фейнман указал оригинальный подход для решения проблем квантового радиационного затухания, который должен был обрести преимущества концепции самодействия при расчетах радиационных ширин
[17]. В рамках формализма интегралов по траекториям ему удалось исключить из теории переменные квантового электромагнитного поля и сформулировать временную эволюцию взаимодействующих через вакуум зарядов в терминах переменных вещества. Таким образом, совокупность заряженных частиц стало возможным трактовать как динамическую систему с комплексным гамильтонианом Я. Трение за счет излучения включается в мнимую часть Н и количественно описывается соответствующими энергетическими уширениями уровней. Определение радиационных ширин Ад возбужденных состояний в данном подходе сводится к нахождению комплексных уровней энергии гамильтониана Ы
6д -> Ед = Ед ~ ШАд / 2 . (4)
Важное свойство гамильтоновых систем связано с существованием альтернативой формулировки их временной эволюции в терминах теории вариационного исчисления. Обращение к вариационному принципу как к возможному математическому инструменту для расчета комплексных значений Ед встречается, однако, в рамках подхода Фейнмана с двумя принципиальными затруднениями. Во-первых, мнимая часть гамильтониана Я,
9
обусловленная радиационным самодсйствием заряженных частиц, оказывается не Ê-локальной по времени [17] и не вписывается в стандартную схему канонического квантования [18, 19]. Во-вторых, расходимости, связанные с бесконечностью электромагнитной массы точечного заряда, не позволили Фейнману создать релятивистски-инвариантную самосогласованную теорию при стремлении радиуса заряда к нулю ([17], стр. 283).
Такая последовательная теория развивается в первой главе диссертации, где, исходя из представлений Фейнмана, формулируется гамильто-новский принцип наименьшего действия для подверженных самодействию излучающих частиц. В отличие от [17], где релятивистски-инвариантная фейнмановская перенормировка не приводит к однозначно определенному гамильтониану, здесь (раздел 1.1.2) осуществляется по аналогии с классическим случаем 110] разложение полного лагранжиана Ltot по параметру запаздывания взаимодействия между зарядами [22, 21]. В полученном разложении мнимой части Ltot указываются члены, несущие в себе черты классической картины самодействия [22], а также выделяются факторы, обусловленные спецификой некоммутативности физических величин в волновой механике [21]. Везде последовательно используется нерелятивистская кулоповская калибровка, в рамках которой стационарно движущиеся атомные частицы создают в основном собственное магнитное поле ~Ё. В отличие от классических сил Лоренца квантовая теория приводит к нелокальным по времени физическим величинам. Прежде всего это сказывается на векторном потенциале самодействия /?, в формировании которого принимают участие все точки, лежащие на траектории движения зарядов (см. формулу (1.18) на стр. 34). Определение энергетического спектра ^-нелокальных гамильтонианов выходит за пределы применимости традиционных канонических схем квантования, что однако не мешает воспользоваться методами квазиклассического приближения. Для этого формулируются многомерные правила квантования (раздел 1.2 главы 1) и осуществляется построение квазиклассической вариационной процедуры [23] применительно
к расчетам радиационных времен жизни. Наличие комплексного поля са-модействия ~Й приводит к появлению мнимых добавок у атомных переменных действия, что в конечном итоге и вызывает сдвиг энергетических термов в комплексную область. В полных вероятностях радиационного распада выделяются два члена Ад = + Адди\ первый из которых обусловлен
собственно самодействием, а второй возникает в силу нулевых флуктуаций элетромагнитного поля в физическом вакууме и имеет квантовую природу своего происхождения. Для вероятности Адс1^ формулируется принцип соответствия (соотношение (1.44) на стр. 45), позволяющий пересчитывать мощность классических диссипативных сил в вероятности квантовомехаии-ческих переходов. Вариационный принцип оказывается эффективным техническим приемом для нахождения регуляризованных значений радиационной ширины 5-состояний водородоподобных атомов и ионов в ситуации когда применимость нерелятивистского приближения нарушается (раздел 1.3). В частности, показано, как с помощью законов подобия Кеплера для кулоновских траекторий можно перейти от условного принципа минимума Персиваля к безусловной вариационной процедуре по определению радиационных времен жизни 5-состояний.
Концепция сил самодействия, развитая в первой главе, имеет много общего с проблематикой квантования классических систем с диссипацией. Силы трения, действующие на движущиеся частицы, не являются потенциальными, и уравнения движения механической системы с диссипацией оказываются за пределами применимости Гамильтон о веко го формализма. Несмотря на более чем полувековую историю своего исследования, проблема описания сил трения на квантовом уровне, остается во многом еще открытой. Даже простейшее уравнение движения
= ’ и(г) = '^Аг2 (5)
И
для трехмерного гармонического осциллятора с линейной по скорости силой трения порождает до сих пор множество идеологически различающихся работ. Среди подходов, претендующих на построение последовательной квантовой теории диссипативных систем, можно выделить два основных направления. Для первого из них характерно моделирование процессов диссипации за счет введения внешней среды ("термостат") В, с которой квантовая частица С обменивается энергией. Взаимодействие внутри замкнутой системы В+С задается при помощи феноменологических констант связи. Кроме того, задается некоторая статистика для резервуара В. Возникающая на этом пути квантовая динамика частицы С формулируется в терминах Гейзенберговских операторов развития, причем результирующий Гамильтониан II существенным образом зависит от времени (см.,например, одну из недавних работ [24] на эту тему и многочисленные ссылки в ней). Однако, в литературе существует мнение, что получаемые таким образом эффективные гамильтонианы типа Калдирола-Канаи [24, 25] описывают скорее не эффекты диссипации, а имитируют их появление с помощью временной зависимости у массы частицы. Так, полученный в |24] результат о полной остановке двигающейся частицы находится в полном количественном согласии с экспоненциальным временным ростом ее массы.
Авторы второго направления предпочитают давать более абстрактную формулировку квантовой динамики, стараясь избегать каких-либо модельных представлений о природе сил трения. В качесгве типичных примеров построения диссипативных несамосопряжснных Гамильтонианов следует указать работы [27, 28]. Слабая физическая мотивация подобных работ приводит к целому ряду недостатков в полученных результатах: это или появление нелинейных членов в уравнении Шредингера [27] (что исключает фундаментальный принцип суперпозиции в квантовой механике), или полное исчезновение стабильных состояний в квантовой системе [28].
Во второй главе диссертации обсуждается один из путей решения проблемы квантования энергии частиц с диссипацией на примере гармоничес-
12
кого осциллятора (5). Предлагаемая здесь "излунательная" модель квантования не требует какой-либо конкретизации характера взаимодействия с термостатом и основывается на использовании аналогии с радиационными силами трения, введенными Лоренцем. Соотношения для пересчета мощности диссипативных сил в вероятностную скорость квантовых переходов "заимствованы" из принципа соответствия, который сформулирован на стр. 45 диссертации (раздел 1.2.1) для классической части спонтанного излучения. Нетривиальный момент нашего подхода связан с возможностью естественной реконструкции комплексного тамильтониана для системы (5) (с произвольным потенциалом и (г)), мнимая часть спектра энергий которого дает количественную информацию о скорости диссипативных процессов. Исследование случая гармонического осциллятора, приведенного в разделе 2.3, позволяет продемонстрировать стабильность основного состояния соответствующей квантовой системы и проследить предельный переход к классической механике.
Содержание третий главы диссертации в отличие от предыдущих двух глав имеет более прикладной характер и посвящен вопросам комплексного исследования радиационных времен жизни тя возбужденных состояний атомов и ионов. В первых двух разделах главы проведена систематизация имеющихся в литературе данных по тя атомов, сопоставлены разные теоретические и формульные методы их расчета тя. Далее общая теория радиационного самодействия "адаптируется" к конкретному случаю расчета радиационных ширин для возбужденных состояний в условиях выполнимости одноканального приближения (раздел 3.3.1). В рамках квазиклассичес-кого приближения устанавливается однозначная связь между резонансной конгруэнцией классических траекторий движения, соответствующих волновой функции | д ), и энергетическим спектром Ея. В результате естественная ширина Лд, будучи мнимой частью полной энергии ЕЯ) становится совместно с положением энергии ея характеристиками только излучающего
13
состояния 4. Подчеркнем, что в нашем подходе отсутствуют индивидуальные переходы между возбужденными уровнями, и при определении Тд = Л“1 отпадает необходимость знать волновые функции всех нижележащих состояний. С формальной точки зрения энергия Ед является диагональным матричным элементом комплексного гамильтониана, что обеспечивает высокую точность квазиклассической схемы расчета. Этот факт демонстрируется на примере естественных ширин возбужденных состояний атома водорода. Полученные для них аналитические формулы (раздел 3.3.4) описывают численные квантово-механические расчеты с точностью до доли процента. Нестандартная ситуация возникает для 5-состояний серий. Соответствующие им орбиты движения оптического электрона проходят вблизи кулоновского центра притяжения где нарушаются условия нерелятивистского приближения. Формально это отражается в появлении двух расходящихся параметров в формулах теории. Последний, однако, удается восстановить исходя из вариационного принципа и используя спектроскопические данные по метастабильным состояниям [33]. Интересный круг явлений связан с соотношением классических (самодействие) и квантовых (нулевые флуктуации вакуума) факторов в излучении. В разделе 3.3.5 формулируется критерии выполнимости правил Бете для оптических переходов и прослеживается их связь с возникновением особенностей в спектре спонтанного излучения. Иллюстрация выводов дается на примере движения электрона в модельном потенциале Зоммерфельда К
На основе полученных аналитических выражений предлагается экстраполяционная схема продолжения радиационных времен жизни с нижних уровней серии на высшие вплоть до ионизационного предела (раздел 3.5). Здесь полезным методическим приемом оказалась возможность сопоставления естественной ширины уровня с мнимой частью квантового дефекта серии. В частности подобная аналогия позволила восстановить парциальные сечения фоторекомбинационных процессов для медленных тепловых
’кулоновский потенциал с добавкой типа центробежного отталкивания: ир(г) = —г/т + р2/(2г2).
14
электронов на счет информации по радиационным временам жизни (раздел 3.5.3).
Материал первых трех глав диссертации касался в основном интегральных характеристик интенсивности спонтанного излучения. Между тем большой интерес представляют вероятности оптических переходов ЛЯ)д между отдельными возбужденными состояниями атомов и ионов. Информация о дипольных матричных элементах Ядя> = ( | г | ), появляющихся
в квантово-механических формулах для Ад>я, необходима также для решения широкого круга задач. Наряду с тем, что величины Я(]({> определяют эффективность всевозможных радиационных процессов для связанно-связанных, связанно -свободных и свободно-свободных электронных переходов, они входят также в большое число соотношений для расчетов сечений разного рода столкновительных реакций. Точные аналитические представления известны только для единственного случая атома водорода и водородоподобных ионов |34). Соответствующая формула Гордона содержит гипергеометрические функции с набором индексов, сложным образом зависящих от главных и орбитальных квантовых чисел исследуемых состояний. Поэтому даже она трудно обозрима и позволяет разумные упрощения для ограниченного круга условий [35, 36, 39). В связи с этим ряд авторов, начиная с работ [37, 38|, предлагают для расчета дипольных матричных элементов использовать квазиклассическое приближение. Ведущий член квазиклассических формул для Я^ определяется Фурье-разложением радиальной координаты электрона при его классическом движении но некоторой средней кеплеровской орбите дс с энергией £с и орбитальным моментом Ьс. Для переходов между состояниями с близкими квантовыми числами д = {п*,/} , д = {п*,^} параметры орбиты дс = { ес Ьс} легко находятся но д (атомная система единиц):
_ г2 е + £ Т_Ь'\-Ь 1 + 1+1
£с ~ ~2(+У - 2 ’ с ~ ~ 2 ’ ( )
где ^-эффективный заряд ионного атомного остова.
Квазиклассический орбитальный момент Ь электрона отличается от его соответствующего квантового значения на поправку Лангера Ь = I + 0.5
[18]. Поэтому стандартное определение Ьс для средней орбиты qc берется как среднее арифметическое значение для Ь и Ь'. Главным препятствием для применимости квазиклассического подхода в случае далеко отстоящих уровней при 5 = \п * — п*| >> 1 является большая неопределенность в выборе энергии ес (или соответствующего эффективного квантового числа п*) для средне-взвешенной орбиты которая должна вбирать в себя характеристики каждого из исходных состояний с[. Указанная неопределенность для п* породила большое количество различных рецептов (правил соответствия) как найти п*. Отметим, что все они были сформулированы без какой-либо серьезной физической мотивации и основаны во многом на эмпирических данных. Ситуация с орбитальными квантовыми числами проще: поскольку I и I' отличаются только на единицу, указанный выше выбор Ьс (6) является оптимальным (40, 41].
В четвертой главе диссертации представлены исследования работы
[42], отвечающие на вопрос существует ли фундаментальный принцип для определения квантового числа дс средней кеплеровской орбиты. В качестве такого принципа берется требование ортогональности квазиклассических волновых функций, принадлежащих одной серии. Однако условие ортогональности фиксирует целый набор значений п*, что обусловлено возможностью восстановления фазовых факторов только с точностью до 27г. Возникшая неопределенность снимается за счет привлечения первого принципа соответствия Бора (18). Полученное квантовое число п* (соотношение
(4.9) на стр. 116) позволяет, после подстановки его в квазиклассические формулы Буреевой-Давыдкина-Зона [37, 38], решить вопрос о расчетах сил осцилляторов для длинных связано-связанных переходов. Более того, структура соответствующих формул оказалась приспособленной также для их продолжения в континуум энергий серий, т.е. для количественного
16
исследования вероятностей в случае связано-свободных переходов. Далее (раздел 4.1.2) указана возможность существенного улучшения полученных соотношений за счет введения в них дальнейших квантовомехаиических поправок, полученных в работах [35, 39), что демонстрируется на примере расчета сечений фотоионизации возбужденных состояний 5-серий для атомарного натрия и однократного иона стронция в условиях возникновения куперовского минимума (раздел 4.2). Наконец в последнем разделе четвертой главы представлены результаты расчета фоторекомбинационно-го спектра медленных тепловых электронов при моделировании атомного потенциала для связанных состояний потенциалом Зоммерфельда. Особое внимание уделяется анализу случая появления траекторий типа "глории"
[43) когда угол рассеяния электронов обращается в нуль. Здесь в соответствии с критерием подавления классических факторов в радиационных процессах наблюдается существенная деформация спектра фоторекомби-нации. Последнее явление выражается в возникновении тенденции к созданию инверсной заселенности нескольких нижних возбужденных состояний с главными квантовыми числами п ~ 3 + 4.
Приближение изолированных атомов перестает выполняться по мере роста плотности газовых сред. Резонансный квант света может с заметной вероятностью быть поглощен (пленен) атомами в нормальном состоянии и покидает объем с газом лишь после нескольких актов своего захвата с последующим испусканием [44). Указанное явление играет существенную роль во многих разделах физики, начиная от астрофизических исследований [45], где перенос лучистой энергии является одним из важных факторов, формирующих спектры небесных тел, и кончая современным лабораторным экспериментом [46, 47, 48), количественный анализ результатов которого встречается с необходимостью корректного учета явления пленения излучения |49, 50, 51). Для теории переноса излучения характерно развитие двух направлений. Прежде всего, это изучение явления пленения
17
как самостоятельного физического феномена, вывод основных уравнений, описывающих ансамбль возбужденных состояний и взаимодействующего с ним ноля излучения (см., например, недавние работы [52, 53] и библиографию в них). Необходимо отметить, что получение решений возникающих интегро-дифференциальных уравнений сталкивается с большими математическими трудностями даже в простейшей двухуровневой модели атома и предположении о полном перераспределении по частотам испускаемых фотонов [54, 49]. Главным образом это обусловлено возможностью обмена световой энергией между удаленными атомами за счет высвечивания квантов света в крыльях линии, прозрачных для газовых сред. Формальным следствием существования длиннопролетной части спектра в уравнениях типа Бибермаиа-Холстейна является расходимость средней длины свободного пробега фотона, что делает невозможным приближение интегральных уравнений пленения локальными диффузионными уравнениями типа Фоккера-Планка.
В связи с этим в научной литературе выделилось самостоятельное, второе направление работ, посвященных разработкам специальных методов, как аналитических, так и численных, решения уравнений переноса [55, 51]. Однако, несмотря на интенсивную деятельность многих исследователей на протяжении последних пятидесяти лет, аналитические достижения оказались весьма скромными [56]. Основные неприятности содержатся в пространственной части уравнений, что обуславливается дальнодействующим характером передачи энергии возбуждения за счет обмена фотонами. Даже рассмотрение простейших одномерных геометрий поглощающего объема (бесконечный слой, бесконечный цилиндр, шар) породило огромное число весьма специфических методик и приемов решений возникающих задач [45, 55, 51]. Учет квантово-механических аспектов проблемы в некотором условном смысле оказывается проще. Так например уравнения Дьяконова-Псреля [57], которые описывают широкую совокупность физических параметров, сопутствующих переносу излучения (наличие зеемановской структуры уровней, поляризацию излучения, перераспределение возбужденных
18
атомов по скоростям, эффекты корреляции между частотами испущенного и поглощенного фотонов), имеют явное аналитическое решение [58] в бесконечном пространстве (см. также [47]). Между тем именно факт конечности размеров газовых ячеек определяет важные особенности спектра эффективных радиационных констант и их влияние на измеряемые в эксперименте физические величины (см. , например, [59]).
В связи со сказанным представляет практический и научный интерес разработка принципиально новых аналитических подходов, которые сочетают в себе, с одной стороны, простоту выполнения вычислений, а, с другой стороны, являются достаточно универсальными. Изложению такого метода, получившего название метода геометрического квантования [60, 61, 62], посвящена пятая глава диссертации. Метод геометрического квантования (МГК) представляв!' собой различные модификации многомерного квазик-лассического приближения, которые приспособлены к построению решений интегральных уравнений в довольно сложных конфигурациях объема. Здесь уместно отметить, что хотя возможности МГК продемонстрированы в диссертации для класса задач, формулируемых в рамках уравнения Бибермана-Холстейна, он имеет хорошие перспективы как стать базой для разработки новых численных методов типа волновых пакетов [63], так и для проведения дальнейших аналитических исследований более реальных физических задач, включающих в себя описание эффектов частичного перераспределения по частотам.
Основная идея МГК [60] заключается в проведении аналогии между интегральными уравнениями пленения и волновыми уравнениями квантовой механики. Уравнение Бибермана-Холстейна предлагается рассмотреть как разновидность уравнения Шредингера для некоторой трехмерной классической гамильтоиой системы, которую условно можно назвать ассоциированной квазичастицей. Кинетическая энергия квазичастицы однозначным образом определяется спектральными характеристиками газовой среды и имеет сложную зависимость от импульса, в то время как потенциальная
19
энергия совпадает с вероятностью тушения уровней (рассмотрение ведется в системе единиц с единичным значением постоянной Планка). В рамках функционального интегрирования (раздел 5.3) доказывается, что дискретный спектр квантовых значений энергии квазичастицы совпадает с полным наборов эффективных радиационных констант распада возбужденных состояний атомов [64]. При этом волновые функции квазичастицы задают пространственную структуру соответствующих мод. Здесь же формулируется принцип "дополнительности”: факт прозрачности стенок дГ2 кюветы
для выходящего излучения вызывает бесконечный перепад потенциала (тушения) на поверхности кюветы. Последний вынуждает квазичастицу упруго отражаться от границы сЮ. Таким образом, проблема определения факторов пленения в П сводится к нахождению уровней энергии квазичастицы, запертой в потенциальном ящике О - типичная задача о квантовом бильярде (см. рис. 5.1 на стр. 140).
В последующих разделах пятой главы диссертации дается собственно обоснование метода геометрического квантования применительно к специфическим условиям задач радиационной кинетики. Формулируются многомерные правила квантования (раздел 5.5), которые фиксируют разрешенные (резонансные) конфигураций волновых фронтов и нолей импульсов квазичастицы для целого ряда конфигураций кювет (раздел 5.7): конечные цилиндры, трехмерные параллелепипеды, призма-подобные геометрии, сферы, эллиптические цилиндры, эллипсоиды вращения. Радиационные константы рассчитываются как кинетическая энергия квазичастицы для "резонансных" значений ее импульсов. Таким образом наш подход позволяет аналитически определять факторы пленения для всех фактически встречающихся контуров линии, включая Фойхтовский при наличии СТС структуры, при произвольном значении оптической плотности поглощающей среды.
К важной особенности развиваемой теории относится идентичность между траекториями движения к в аз и частицы и бильярдными траекториями
20
внутри газовой ячейки. Это позволяет без изменения переносить все результаты традиционных квазиклассических теорий для определения фазовых факторов при прохождении казичастицей каустических поверхностей. Особенности нашего подхода возникают по существу от граничных условий, которые включены в скачки фазы после отражения волны де Бройля квазичастицы от стенок кюветы Ш. Явная формула для Д5 получена в Приложении С. В разделе 5.6 демонстрируется сильная зависимость Д*9 от угла отражения в квазичастицы. В эллиптических бильярдных задачах угол отражения, как оказывается, имеет различные значения на поверхности бильярда. Это обстоятельство учитывается в разделе 5.8 за счет внесения модификаций в правила геометрического квантования. Последние находятся в рамках квазиклассической вариационной процедуры, разобранной в Приложении С (раздел С.2). Приведенные в разделе 5.9 сравнения с численными расчетами факторов пленения указывают на высокую точность МГК (в пределах 5 % для фундаментальных мод и в долях процента для всех остальных высших мод). Обоснование последнего обстоятельства дано в разделе 5.9 на основе анализа поведения волновых функций (мод) квазичастицы в окрестности квазиклассически особых поверхностей (каустик и поверхностей отражения).
В вопросах кинетики радиационного возбуждения и ионизации газовых сред в последнее время большое значение приобретают исследования ударной ионизации с участием ридберговских состояний атомов. Как показывают оценки, относительный вклад реакций типа ассоциативной ионизации в процесс образование первичной ионизации импульсной фотоплазмы может достигать 20 % (4]. В связи с этим в работе |65| были получены аналитические формулы для расчета скоростей ионизации при симметричном столкновении пары: ридберговский атом (А**)-атом в основном состоянии (А). Однако предсказанные теорией эффективности ионизации оказались сильно заниженными по сравнению с экспериментально полученными, и авторы
21