Коллективные ридберговские состояния в верхней атмосфере Земли тема автореферата и диссертации по физике, 01.04.17 ВАК РФ
Голубков, Максим Геннадиевич
АВТОР
|
||||
доктора физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Москва
МЕСТО ЗАЩИТЫ
|
||||
2014
ГОД ЗАЩИТЫ
|
|
01.04.17
КОД ВАК РФ
|
||
|
На правах рукописи
ГОЛУБКОВ Максим Геннадиевич
КОЛЛЕКТИВНЫЕ РИДБЕРГОВСКИЕ СОСТОЯНИЯ В ВЕРХНЕЙ АТМОСФЕРЕ ЗЕМЛИ
01.04.17 — химическая физика, горение и взрыв, физика экстремальных
состояний вещества
автореферат диссертации на соискание ученой степени доктора физико-математических наук
005557467
15 ;;:;з 2015
Москва 2014
005557467
Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте химической физики им. Н.Н.Семенова Российской академии наук
Научный консультант: Авакян Сергей Вазгенович
доктор физико-математических наук Санкт-Петербургский государственный политехнический университет, профессор
Официальные оппоненты: Бычков Владимир Львович
доктор физико-математических наук Московский государственный университет им. М.В.Ломоносова, Физический факультет ведущий научный сотрудник
Лушников Алексей Алексеевич
доктор физико-математических наук Геофизический центр Российской академии наук главный научный сотрудник
Юрова Инна Юрьевна
доктор физико-математических наук Санкт-Петербургский государственный университет, Физический факультет профессор
Ведущая организация: Институт энергетических проблем
химической физики им. В.Л.Тальрозе Российской академии наук
Защита состоится 18 марта 2015 г. в 11 часов 00 минут на заседании Диссертационного совета Д.002.012.02 при Федеральном государственном бюджетном учреждении науки Институте химической физики им. Н.Н.Семенова Российской академии наук по адресу: Москва, Ленинский проспект, 38-1.
С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Института химической физики им. Н.Н.Семенова Российской академии наук.
Автореферат разослан 18 декабря 2014 года.
Автореферат размещен на сайте Высшей аттестационной комиссии Министерства образования и науки Российской Федерации 17 декабря 2014 года.
Ученый секретарь
Диссертационного совета Д.002.012.02
доктор химических наук, профессор Корчак В.Н.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы исследования. В последние два десятилетия всплеск широкомасштабных исследований физико-химических процессов в атмосфере Земли объясняется возросшим интересом к процессам распространения электромагнитных сигналов со спутников и радиолокационных станций. Фундаментальную роль в формировании спектров радиосигналов играют столкновительные и радиационные процессы с участием коллективных ридберговских состояний в нейтральной среде, что определяет актуальность проведенного в диссертации исследования.
Цели и задачи исследования. Целью данной диссертационной работы было установить причину сверхфонового сверхвысокочастотного (СВЧ) и инфракрасного (ИК) излучения, формирующегося в £>-слое ионосферы Земли, в моменты повышения солнечной активности, сопровождающиеся сильными геомагнитными возмущениями. В соответствии с целью работы были поставлены следующие задачи:
1. создать метод расчета поверхностей потенциальной энергии (ППЭ) и волновых функций системы "ридберговский атом А** + нейтральный атом В "в широком диапазоне изменения межатомных расстояний
2. рассчитать полосы эффективного излучения для переходов между уровнями высоковозбужденных квазимолекул А** N2 и А**0-2 в диапазоне частот 0.8—100 ГГц
3. рассчитать заселенности ридберговских квазимолекул А**N2 и А**02 в неравновесной двухтемпературной ионосферной плазме
4. рассчитать частотные профили мощности потока СВЧ излучения
5. проанализировать возможности стимулирования и подавления сверхфонового СВЧ излучения
Научная новизна. В диссертации впервые было показано, что в периоды повышения солнечной активности и сильных геомагнитных возмущений в Е- и £)-слоях ионосферы, расположенных на высоте 50—110 км от поверхности Земли, формируется сверхфоновое СВЧ излучение, интенсивность которого многократно превышает типичные уровни всплесков Солнца. Было показано, что это излучение образуется в результате формирования неравновесной двухтемпературной плазмы за счет переходов между орбитально вырожденными состояниями квазимолекул А**N2 и А**С>2- Их
заселенности определяются концентрацией и температурой среды, а также концентрацией, температурой и потоком свободных электронов.
Впервые установлена определяющая роль нейтральных молекул азота и кислорода верхней атмосферы Земли в формировании спектра спонтанного СВЧ излучения в дециметровом диапазоне. Обнаружено, что спектр является существенно неоднородным и содержит три диапазона частот, в которых происходит заметное уменьшение интенсивности излучения.
Кроме того, было показано, что радиационные переходы с изменением главного квантового числа должны приводить к сверхфоновому ИК излучению в диапазоне длин волн от 15 до 60 мкм, что было подтверждено результатами измерения спектрометра FIRST. Стоит отметить, что наблюдаемое ИК излучение может быть использовано для детектирования ридберговских состояний.
Теоретическая и практическая значимость работы. В диссертации подробно описаны физические причины появления сверхфонового СВЧ излучения в D-слое ионосферы Земли в периоды повышения солнечной активности, сопровождающиеся вспышками на Солнце и сильными геомагнитными возмущениями ионосферы. Такие возмущения, как правило, приводят к сбоям в работе измерительной аппаратуры, искажению или исчезновению на приемниках спутниковых сигналов.
Традиционные модели ионосферных процессов, опирающиеся на полное содержание электронов и волновую оптику, в этом случае оказываются малоэффективными. Правильный учет оптических резонансных квантовых свойств нейтральной среды нижней ионосферы позволяет выделить и использовать такие частотные диапазоны, в которых нарушения работы приборов оказываются минимальными. Кроме того, в диссертации рассмотрены возможности локального усиления и подавления сверхфонового СВЧ излучения.
Методы исследования. Процедура определения одноцентровых операторов рассеяния слабосвязанного электрона ридберговского атома на изолированных центрах развита в рамках обощенного метода псевдопотенциала конечного радиуса. Анализ особенностей поведения поверхностей потенциальной энергии системы проводился с использованием интегрального варианта теории многоканального квантового дефекта (МКД). Сравнение с результатами квантовохимических расчетов методом MRD CI показало, что для небольших значений главных квантовых чисел эти методы взаимно дополняют друг друга, так как первый оказывается более надежным на больших межатомных расстояниях, а второй — на малых. Для корректного описания ППЭ подобных систем предлагается использовать метод сшивки.
Для описания взаимодействия медленно сталкивающихся ридбергов-ского атома А** с нейтральным атомом В, находящимся в основном электронном состоянии, развит метод комплексного оптического потенциала. В основе метода лежит асимптотический подход с привлечением теории МКД, использующий формализм перенормированных уравнений Липпмана-Швингера.
Расчет низкотемпературной реакции диссоциативной рекомбинации, протекающей во внешнем поле монохроматического лазерного излучения, проведен в рамках теории МКД с использованием стационарного формализма матрицы радиационных столкновений.
Положения, выносимые на защиту. На защиту выносятся:
1. Метод расчета поверхностей потенциальной энергии и волновой функции системы, состоящей из высоковозбужденного атома и нейтрального атома с заполненной электронной оболочкой.
2. Механизмы формирования двухтемпературной квазистационарной неравновесной плазмы, которая образуется в верхней атмосфере Земли под действием потока электронов, высыпающих из ионосферы в периоды повышенной солнечной активности.
3. Расчет спектра некогерентного излучения в диапазоне 0.8—100 ГГц, возникающего за счет радиационных переходов между возмущенными средой орбитально вырожденными состояниями ридберговских комплексов, происходящих без изменения главного квантового числа.
4. Решение задачи о генерации сверхфонового СВЧ излучения в О-слос ионосферы Земли, наведенного точечным импульсным 7-источником.
5. Теория лазерного управления реакцией диссоциативной рекомбинации медленных электронов и молекулярных ионов, отвечающей за кинетику исчезновения зарядов в ионосферной плазме. Математический аппарат теории основан на формализме матрицы многоканального рассеяния в представлении квазиэнергетических состояний.
Степень достоверности полученных результатов. Достоверность полученных результатов обеспечена использованием известных апробированных теоретических подходов. Правильность сделанных предположений относительно параметров ионосферы подтверждается достаточным количеством экспериментальных данных, измеренных с помощью баллистиче-
ских ракет. Результаты расчета мощности потока СВЧ излучения подтверждаются данными спутниковых измерений.
Апробация результатов. Результаты исследований, представленных в диссертации, докладывались и обсуждались на: IV International Conference "Dissociative Recombination: Theory, Experiment and Applications" (г. Стокгольм, Швеция 1999); International Conference "Highly Excited Electronic States" (г. Барселона, Испания 1999); International Seminar on Physics of Electronic and Atomic Collisions (г. Клязьма, Россия 2001); International Symposium "Dissociative Recombination of Molecules with Electrons "(г. Чикаго, США 2001); International seminar "Molecular Interactions and differential scattering"(г. Андреасберг, Германия 2002); International Conference "Electron Excited Complexes "(г. Афины, Греция 2002); XIV—XXVI Всероссийских симпозиумах "Современная химическая физика"(г. Туапсе, Россия 2002— 2014); VII Voevodsky Conference "Physics and chemistry of elementary chemical processes "(г. Черноголовка, Россия 2007); XIX и XXI International Conference on Spectral Line Shapes (г. Вальядолид, Испания 2008 и г. Санкт-Петербург, Россия 2012); I—IV International Conference "Atmosphere, ionosphere, safety"(r. Калининград, Россия 2008, 2010, 2012, 2014); European Conference on Dynamics of Molecular Systems (г. Санкт-Петебург, Россия 2008); 40th Scientific Assembly COSPAR (г. Москва, Россия 2014).
Публикации. По материалам диссертации опубликовано 30 печатных работ, включающих 3 коллективные монографии. Работ, опубликованных в рецензируемых научных изданиях, рекомендованных ВАК — 26.
Личный вклад автора. Все исследования, описанные в диссертации, выполнены лично автором. Во многих случаях эти работы проводились в рамках сотрудничества с другими российскими и зарубежными учеными. Постановка задач и целей исследований, обсуждение и обобщение полученных результатов, формулировка выводов и подготовка публикаций принадлежит лично автору данной работы. Многие работы автора получили широкое международное признание. Лично автором была создана вычислительная программа "Ридберг", с помощью которой были получены основные результаты диссертации.
Объем и структура диссертации. Диссертация состоит из введения, пяти глав, формулировки основных результатов и выводов, списка сокращений и списка цитируемой литературы. Работа изложена на 270 страницах и содержит 44 рисунка, 17 таблиц и библиографию из 172 наименований.
ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ
В первой главе рассматриваются системы "ридберговскии атом ^4** -атом В". При описании подобных систем обычно применяется асимптотическое приближение с использованием моделей потенциала нулевого радиуса (ПНР) или конечного радиуса (ПКР), котрые не учитывают дально-действующие потенциалы: поляризационный (в случае электрон-атомного взаимодействия) и электродипольный или электроквадрупольный (в случае взаимодействия с молекулами) потенциалы. Они имеют простую форму только на больших расстояниях электрона от рассеивающего центра. На промежуточных расстояниях возникает нетривиальная проблема учета взаимодействия иона с возмущающими центрами и устранения неопределенности в разложении действующего на слабосвязанный электрон эффективного потенциала на короткодействующие и дальнодействующие части.
Предложен метод построения таких одноцентровых операторов рассеяния, которые позволяют рассчитывать поверхности потенциальной энергии (ППЭ) системы "ридберговский атом А** - атом В" на "промежуточных" расстояниях между ними с учетом их реальной структуры. С помощью уравнений теории многократного рассеяния полученные результаты естественно распространяются и на многоцентровые системы А** — Bs, где индекс s предполагает суммирование по возмущающим нейтральным центрам. Под "промежуточными"следует понимать такие межъядерные расстояния, на которых внутренние электронные орбиты взаимодействующих центров не перекрываются.
Высоковозбужденный атом A**(nlm) является электронейтральной частицей и представляет собой ион А+, окруженный протяженным электронным облаком, размеры которого существенно превышают размеры невозбужденных атомов. Здесь п - главное квантовое число, 1\\т- орбитальный момент ридберговского электрона и его проекция на ось, направленную вдоль вектора R, соединяющего центры тяжести частиц А+ и В. Энергия связи этого электрона Еп неисонзмеримо мала по сравнению с потенциалом ионизации Iq атома А.
При достаточно высоких возбуждениях (п2 1) существует значение углового момента электрона /*, которое за счет центробежного барьера делит ридберговские состояния на сильнопроникающие и слабопроникающие в ионный остов. Для большинства атомов I* = 3. Состояния с малыми угловыми моментами I < I* сильно взаимодействуют с ионным остовом и характеризуются наличием квантовых дефектов ¡ц, которые приводят к сдвигу Дп; ридберговских уровней относительно энергии Eni кулоновского
центра, т.е. в атомной системе единиц (Л = е = те — 1)
Еп1 = -2(п-щУ=Еп + Ап1> (1)
Состояния с большими угловыми моментами 1^1* с хорошей точностью можно считать кулоновскими (орбитально вырожденными), поскольку квантовые дефекты /х; с ростом I убывают как щ ~ I"5.
Решение задачи о нахождении поверхностей потенциальной энергии (ППЭ) системы А**(п1т)+В основано на формализме перестроенных уравнений Липпмана-Швингера для оператора сдвига уровней
т = Ке-_В(СА» - Со)т, (2)
где за счет разделения функции Грина ридберговского атома А** на сильно и слабо зависящие от энергии части
Сл..(г, г', Е) = тги\Е) ]Г |Ф^(г))(ф№)1 ^тг(и(Е) + щ) + С0(г, г'),
1т
(3)
в теорию вводится матрица Ке--в рассеяния, которая строится на стоячих волнах и удовлетворяет интегральному уравнению
К е—в = Уе—В + Уе—вС0Ке-_в, (4)
где — оператор е~ — В взаимодействия. Элементы Ке-_д -матрицы
выражаются через характеристики рассеяния свободного электрона на атоме В и могут быть рассчитаны независимо.
Вследствие вырождения уровней при I ^ I* полюсную часть функции Грина (3) можно разбить на два слагаемых:
СА..-Со = тш3 £ |Ф^(г)>(Ф^(г)|^7г(1/ + №)+
1<1',т
+ £ |ф^(г)><ф^(г)|, (5)
1^1*,т
где Ф^(г) — волновая функция ридберговского атома. Полюса первого слагаемого в (5) описывают сильно проникающие в ионный остов ридбер-говские состояния (с щ 0), а второго - слабопроникающие кулоновские состояния (с /х; = 0). Сильно проникающие ридберговские состояния воз-
мущаются следующим образом:
AEvl = 2ira\<S>BJ0d(R)\\ (6)
где а — длина рассеяния электрона.
От группы орбитально вырожденных состояний отщепляется ковалент-ный терм, который, согласно (5), определяется уравнением
tg7Ti/ = 27ri/3a5>№)|2- (7)
i^i'
Заметим, что ковалентные состояния, которые формируются из орбитально вырожденных компонент состояний ридберговского электрона, могут классифицироваться по значениям углового момента L электрона относительно возмущающего атома В.
В общем виде Ке-_в-опреатор записывается в виде:
К е-_в(р, R) = Y, Kf](p, R) + U(p, R), (8)
L
где первое слагаемое в правой части выражения (8) соответствует короткодействующей, а второе — дальнодействующей части псевдопотенциала рассеяния.
В случае L = 0, отвечающего модели ПНР, считается, что сильное взаимодействие Ve--B электрона с атомом В сосредоточено в окрестности точки R радиуса р = |г — R| < ро и волновая функция Ф^(г) слабо меняется на масштабе ро. Взаимодействие е~ — В можно представить в виде псевдопотенциала Ферми:
Ve—B = 2тга<5(г - R), (9)
который позволяет выйти за рамки приближения ограниченного базиса, что выгодно отличает излагаемый ниже подход от традиционно используемых в квантовой химии.
Из условия обращения в нуль определителя системы уравнений для элементов оператора сдвига уровней (2) получается секулярное уравнение
'tg7r(l/ + №)-40)'
которое описывает поведение термов системы A**(nlm) + В на промежуточных расстояниях. Матричные элементы оператора К^_дзаписываются
(tgvv - кт = (кЩ , (10)
как
= («Ж £ ^ = 2 (и)
Функция определяется как
1риЬ = ч/^зф^ = Е а„Ро1р*$, (12)
О/'
где а^/о — коэффициенты разложения. Если ось квантования направлена вдоль вектора И, то для взаимодействия (9) эффективно представлены только состояния с т = 0. При этом вырожденные состояния с отличными от нуля проекциями т являются кулоновскими.
Основные достижения теории взаимодействия ридберговских атомов с нейтральными частицами связаны с применением модели ПНР, соответствующей учету только первого слагаемого ( Ь = 0) в (8). Выходя за рамки этого приближения, следует включить более высокие угловые моменты Ь рассеивающегося электрона. Обобщение на гармоники более высокого порядка возможно в расширенном варианте теории с учетом конечных размеров радиуса взаимодействия ро в области действия силового поля атома В. Учет состояний с I ^ 1 должен приводить к снятию вырождения для состояний с отличными от нуля проекциями углового момента электрона т на направление вектора И.. При этом на промежуточных расстояниях И, ~ п сдвиги уровней с \т\ = 1 будут соизмеримыми со сдвигами уровней, создаваемыми Р-гармониками рассеяния в состояниях с ш = 0.
При рассмотрении эффектов дальнодействия в системе А** — В следует уделить внимание роли второго слагаемого в правой части выражения (8). Представим эффективное взаимодействие е~ — В в виде суммы двух слагаемых
Уе-В(р, К) = У^_в(р) + иА+е-в(р, И), (13)
Т г(0)
где сильное Уе__в-взаимодеиствие сосредоточено внутри ограниченной области р ^ ро- Далыюдействующая часть силового поля (/л+е-в будет содержать поляризационные слагаемые, включая поляризацию атома В электроном, и может быть учтена по теории возмущений. Для случая ЬБ-связи в основном ¿'-состоянии атома В (с нулевым значением орбитального мо-
мента электрона) второе слагаемое в (13) выглядит следующим образом:
&иА+е-в(р,К) =
ррК
(рД)3'
(14)
Применение теории возмущения по взаимодействию и^+е-в и устранение неопределенности в разделении Уе--в{р, К) на короткодействующую и дальнодействующую части предполагают выполнение следующих условий:
(/3|Я|,/3/Д2)«1, /3/^«1. (15)
В результате матричные элементы И) оператора принимают вид
К,
тг1т,п'1'т
ф
,Яу<1
и1т
Е*
где первое слагаемое в (16), отвечающее Ь = О , записывается как
2па0Ф^(К)Ф^т(К) (17)
Здесь ао отвечает длине рассеяния на короткодействующей части потенциала (13). Так как электронные волновые функции слабо меняются на масштабе ро, то интеграл от поляризационного члена —/?/2р4 в потенциале (14) равен -2п{Р / р0)Ф^{В,)Ф^т{П) и матричные элементы К„1т,п1'т{*1) принимают вид
где а = ао — /3/ро — наблюдаемая длина рассеяния и АКп1тХ1,т(11) = 2тг(Д//0о)Ф^(Н)ФХ(а)^о+
,ПТТП5
(18)
+ ( Ф
ру<1 и 1т
Ь> 1
+ (19)
С точностью до членов ~ /3/ро <С а этот результат не зависит от величины Ро и неопределенность в ее выборе исчезает.
В окрестности возможных сближений или квазипересечений термов рид-
берговских и ковалентного состояний секулярное уравнение принимает вид:
Е - U™(R) [Е - UffiJR)] = {Кп1т,п'Пш? ,
(20)
т.е. взаимодействовать могут только состояния с одинаковыми значениями т. В результате для ридберговского и^ и ковалентного термов
получаем:
тад) 1
27Г а
№
6т0 + А К,
nlm,nl'm
(.R)
J_ 2 /г4'
___
2 [п-(1пиЩ2 2Д4'
HnLm(R) =--arctg [7W3KnLm,nLm(R)] ,
(21)
где элементы ДKnimt„i>m задаются выражением (19), а величина Knim^nim равна:
KnLm,nLn{R) = J2 + AKnlm,nlrn(R)} . (22)
W
Взаимодействие между ридберговским и ковалентным термами определяется как
AVnlm,n'Lm = 2 |KnimyLm\ cos2 ni/(E), (23)
где величина
^nim.n'Lm = J2 + Д^т,„Тт(Д)] • (24)
Орбитально вырожденные термы находятся из системы уравнений:
= Ke-_B™3Ctg7n(Ф.
Ryd\_
Wm ' '
(25)
условие однозначного разрешения которой дает
Det 18iv tgTTv(E) - nnzKnim>nVm\ = 0.
(26)
В случае L = 0 недиагональные элементы обладают свойством сепарабельности
-(о) тАо)
(к^ = Kv
V nlm,nl'ml nlm,nlm nl'm,nl'm'
и
при наличии которого от кулоновских уровней отщепляется только один вырожденный уровень. Учет следующих гармоник (Ь ^ 1) приводит к появлению ветвей из N отщепленных термов, отвечающих смешанным по I группам состояний для заданного значения т.
В качестве конкретного примера в диссертации представлены результаты расчетов ППЭ систем Ыа**—Не и К**—Не с использованием интегрального варианта теории МКД и метода псевдопотенциала конечного радиуса. Проведено подробное сравнение с результатами расчетов других авторов. Отмечено, что в области расстояний Я ^ п наиболее предпочтительными являются результаты, полученные в рамках обобщенного метода ПКР, поскольку точность расчетов здесь контролируется точностью определения амплитуд рассеяния, которые обычно рассчитываются независимо или могут быть определены из эксперимента. С ростом п классификация адиабатических термов без предварительного построения диабатической картины становится затруднительной вследствие увеличения числа возможных пересечений. Существующие же квантовохимические программы, которые используются другими авторами ориентированы только на расчеты адиабатических термов. Достоинством обобщенного метода ПКР заключается в том, что он позволяет одновременно находить как адиабатические, так и диабатические потенциальные кривые и строго устанавливать структуру адиабатической картины термов со всем многообразием квазипересечений.
Для построения правильной картины потенциальных кривых квазимолекулы А**(п1) — В следует осуществлять коррекцию квантовохимических расчетов прямым совмещением (сшивкой) с результатами обобщенного метода ПКР на тех промежуточных расстояниях, где отсутствуют квазипересечения ридберговских и орбитально вырожденных ковалентных термов. В итоге происходит нужное смещение рассматриваемой потенциальной кривой на малые расстояния по энергетической шкале. Найденные таким образом термы на малых межъядерных расстояниях должны строго переходить в потенциальные кривые соответствующих ридберговских молекул. В результате получается согласованная структура адиабатических кривых, на основании которых полностью восстанавливается диабатическая картина термов во всей области изменения межцентровых координат (со всем многообразием квазипересечений). Использование диабатических термов требуется для проведения корректных расчетов зависимостей квантовых дефектов от межатомных расстояний, которые входят в определение матричных элементов.
В заключительной части первой главы строятся полные волновые функции рассматриваемой двухцентровой системы А** — В в зависимости от межатомного расстояния Я. Исходя из соотношения формальной теории
рассеяния, имеем
ФГ) =
= V,
е~—В
|Фыт>
Точную волновую функцию можно представить в виде
№п1т) = V-1_в Т
=
(28)
(29)
При этом функция Грина ридберговского атома для Е < 0 записывается в виде суммы двух слагаемых
(г, г', Е) = С£.(г, г', Е) + С£'.(г, г', Е),
(2)
(30)
где
по+Лп
Е — Е„1
с« (г, г >,е)= Е
по—Дп 1т
С%(г,г',Е) = 60%(т,г',Е)+
У1т[~7),
(31)
— Еи1
ь ■
Здесь п0 соответствует энергии .Е,,; ридберговского атома, — часть
функции Грина, неучтенная в при суммировании по дискретным состояниям, г) — радиальная волновая функция дискретного спектра. Функция в(г, г', Е) задается выражением
С« (г, г\Е) =
1т
ctg7г (^(£) + -
, (32)
которое получается при суммировании по интервалу [п0 — Ап,щ + Дп]. При этом величина Дп зависит от значения щ. Волновые функции непрерывного спектра нормированы как
(ДгНДл) = 7Г^(£ - е')
(33)
и представляют собой "стоячие" волны
Rd(r) = Ms) I F£l(r)
(-1)' 2Л(е)
Tu(e)e-^QEl(r)+
+ TfcyMQMr)
v=ijk
Здесь Ai(e) = cos ¿/(e) — коэффициент, ответственный за нормировку волновой функции при наличии ионного остова, FEi(r) — кулоновская волновая функция сплошного спектра, Тц — амплитуда упругого рассеяния электрона на ионном остове:
Ти(е) = - =
tg7r[i/(e) -> i/k]+tg6t(e)
где
Ti =
О
"*/*) Sin<5;(£)
<pi = arctg
/к + е-2тт/к _ 2 cos 2Si{e)' g^A -j- e_7rA
э7Г/А: _ р—7г/к
tg5i(e)
(35)
(36)
¿¡(б) — фаза рассеяния, связанная при малых энергиях с квантовым дефектом уровня соотношением Ситона = пщ. В результате Д£;(г) равно
R£i(r) = cos6l(e){F£l(r)-
2kltgSl(e)
гл/е2ж/к + е~2ж!к - 2 cos 5г(е)
=Re
ilWi/Kl+1/2(-2ikr)
V(i/k — I)
, (37)
где \¥1/<1+1!2{2г/и) — функция Уиттекера.
Выражение для аналогично С^1,'. и включает суммирование по
нижележащим состояниям с п = 1 Ч- (щ — Ап — 1) и вышележащим с п = (щ + Ап + 1) Ч- оо. Перенося второе слагаемое в правой части (32) из в гладкую по энергии часть г') для функции во окончательно
имеем
С0(г,г\Е) = -^>2(Я) |Ф^(г)) (Ф^(г')| + (38)
1тп
Если атом В находится вдали от области псевдопересеченнй термов с различными I < Г, то для определения волновой функции возмущенных
ридберговских состояний можно воспользоваться однополюсным приближением
Ф:
/фПт|ФГ\
Яу/1\ _ \ у1т I и1т/ п1т/ ~ Е- Ей
ФГ) + СО{Е)т ФГ). (39)
" и1т !
ЯуЛ\
и1т.! '
Парциальные ненормированные волновые функции вырожденных состояний определяются по аналогии с (39) как
|ФХ> = (Ф&1ФХ) \фУт) + с°т И™) - V > п
1'>1*
(40)
где число слагаемых при суммировании по V зависит от главного квантового числа п. На промежуточных расстояниях Я ~ п в каждой из п групп взаимодействующих по I состояний при заданном значении т наиболее эффективно выделяется только одно вырожденное ковалентное состояние с энергией Е = ^„¿^(.й), а остальные близки к кулоновским.
Материал второй главы посвящен диффузии ридберговских атомов в атмосфере инородного буферного газа. Если один из атомов сильно возбужден, то наиболее эффективным является использование "трехчастич-ной" модели, когда оптический электрон находится на орбите, размеры которой существенно превышают характерные размеры ионного остова и атомов среды. Тогда при описании упругого рассеяния ридберговских атомов на атомах мишени помимо определения потенциальной энергии их взаимодействия необходимо учитывать и неупругие виртуальные переходы во все возможные дискретные и континуальные состояния слабосвязанного электрона. Естественным подходом к решению этой задачи является введение оптического потенциала, мнимая часть которого характеризует полную вероятность перехода в указанные состояния.
Переход от многочастичного к двухчастичному описанию упругого столкновения атомных частиц в представлении оптического потенциала является одним из традиционных методов в квантовой теории рассеяния. Обычно в теорию вводится нелокальный оператор Л^ и строится многочастичное уравнение для его определения. Такая постановка задачи носит формальный характер и не может быть разрешена без дополнительных предположений о взаимодействующей системе. Для решения задачи необходимо перейти к представлению невзаимодействующих сталкивающихся частиц, определить базис его собственных состояний \д) и вычислить соответствующий матричный элемент Этот матричный элемент и будет представлять собой искомый оптический потенциал Уор(, мнимая часть которого
должна включать полный набор квантовых чисел я, описывающих внутренние состояния взаимодействующей системы, а также зависеть от полной энергии Е и орбитального момента сталкивающихся частиц.
Рассмотрим медленное упругое столкновение ридберговского атома А**(п 1) с атомом В в основном электронном состоянии (здесь п — главное квантовое число ридберговского уровня) и определим оптический потенциал для данной системы. Задача решается в рамках асимптотической теории с использованием интегрального варианта метода многоканального квантового дефекта .
Полная энергия системы Е в атомной системе единиц (Н = е = те = 1)
Е = + (41)
где VI = п — щ — эффективное главное квантовое число, щ — квантовый дефект уровня, I — угловой момент электрона, Е^ — начальная относительная энергия сталкивающихся частиц.
Будем считать потенциал взаимодействия иа+в иона А+ с атомом В известным. Для решения задачи на собственные значения воспользуемся интегральным уравнением для т-оператора сдвига уровней, которое в этом случае принимает вид:
г = ил+.вС(£)т. (42)
Трехчастичный оператор взаимодействия и а+в, который описывает взаимодействующую пару А+ — В плюс свободный электрон е-, является локальным и представляется как
иЛ+в(Н, В/; р, р') = (2тг)6иА+в(Я) £ №. р)> р')\, (43)
ч
где Ы и р — координаты иона А+ и электрона, отсчитанные от центра атома В. Базисная волновая функция \д) выглядит следующим образом:
|<г(К'р)) = (¿)з 6ХР ['(кК + рер)]' (44)
Здесь к = \J1McEk — относительный импульс сталкивающихся частиц А+ и В, Мс — их приведенная масса, ре — импульс электрона.
Тогда для оператора взаимодействия Цд+ в с учетом (44) получим следующее выражение:
ил+в(К, р, р') = (2тг)6{/л+в(ВДК - Н-'Ж? - Р% (45)
где 6(х) — дельта-функция Дирака. При этом матричный элемент
(д\иА+в(В.,В!;р,рГ)\д) = иА+в(Я). (46)
Гриновский оператор О(Е) в (42) описывает систему А** + В с выключенным иА+в взаимодействием, т.е. атом В в уравнении (42) считается взаимодействующим только со слабосвязанным электроном. Для построения этого оператора воспользуемся уравнением Дайсона:
С(Е) = С Л..В(Е) + СА..в(£^е-вС(£), (47)
в котором Уе-В представляет собой оператор е~ — В взаимодействия, а гриновский оператор Са^в описывает невзаимодействующую систему А**+В с заданной кинетической энергией Е);.
Зависящий от полной энергии системы Е нелокальный оператор оптического взаимодействия У^ можно ввести в теорию следующим образом:
т = Чор1(Е)Оа..в(Е)т. (48)
Он удовлетворяет интегральному уравнению
Уор1 = иА+в + иЛ+вСЛ..В [(1 + гКе-в)-1Ке-в] Сл-вУ^, (49)
где Ке-в — вещественная матрица, строящаяся на стоячих волнах. По определению оператор Сл*'в представляет собой свертку:
Сл..в(В.,В!,Е) = ^ У ехр(гкК)Сл» (е - ехр(-гк11')<Лс,
(50)
где Сл" — гриновский оператор изолированного ридберговского атома. Интеграл (50) сводится к интегралу вида:
Сл..в(В, К', Е) = ¿з ^ Сл» (Е - Ек) (51)
где Ек = к2/2Мс. Заметим, что в случае Е < 0 это выражение является строго вещественным, так как в соответствии с (41) функция Грина ридберговского атома Сл» здесь определена при отрицательной энергии и отвечает связанному состоянию электрона во всей области изменения координат электрона и атома В. Для положительной полной энергии системы (Е > 0) функция Сл*. в (51) является, вообще говоря, комплексной, поскольку при Ек < Е она описывает движение электрона в сплошном спек-
тре. Эта область изменения импульса к сталкивающихся частиц включает виртуальные переходы в ионизационный континуум и переходы между его состояниями. Поскольку основной вклад в (51) вносит классически разрешенная область движения, максимальное значение импульса выбирается из условия:
и равно
йи«=у2Мв(£? + 1)> (53)
Ящгп — ' Втт ^ Я ^ 2^ ,
где 1а —потенциал ионизации атома А.
Для функции Грина высоковозбужденного атома Сд-. имеет место следующее общее представление
С л..(р,р',Я,£) = С(0с}(р,р',£) + 2 22Н(з'19*ЛЯ£), (54) 555'(К,е)= ре(е)^{ж1>(£))533>+а1[и{Т1,£)5м05м'а ■ (55)
Здесь
— квазиклассический импульс электрона в поле кулоновского центра,
(57)
— эффективное главное квантовое число, — электронная волновая функция, имеющая вид:
\з) = \lmLM) = МРеР)Г1тп(К/Я)Уьм(в, ф), (58)
где ]ь{х) — сферическая функция Бесселя первого рода порядка Ь, У/л/(0, ф)
— шаровая функция, в — угол между направлениями векторов ре(е) и р, Ь и М — орбитальный момент электрона относительно атома В и его проекция на направление вектора II, — гладкая часть кулоновской функции
Грина:
Матрица «¿¿/(К, г:) в выражении (55) определяется как
= 2*3у/(2Ь + 1){21/ + 1) х (е)х
эш [7гг/(е]]
х ^¿(е, Д)|Цт(К/Д)|2^ (60)
Соответствующие радиальные волновые функции в выражении (60) имеют вид:
Ь = 2к> к- 0 12 Г6П
Они различаются сдвигом фазы на ж/2 и выражаются через функции Уит-теккера:
и{£)) = (б2)
Г(х) — гамма-функция.
Для нахождения в дальнейшем явного вида оптического потенциала требуется разложить (51) и (59) по сферическим гармоникам, т.е.
ьм
^ркгпах
' ^(кЯ)^(кЛ')СА^(Е - Ек)кЧк, (63) о
где Ь и М — орбитальный момент и его проекция на направление К — К'.
Воспользуемся первым порядком теории возмущений и заменим оператор Vорг в правой части уравнения (49) на локальный оператор (45). Тогда оптический потенциал определяется как матричный элемент
Ус* = (д^д), (64)
рассчитанный с базисными волновыми функциями (44). Интегрируя по межатомным координатам и пользуясь свойством симметрии оператора
<5|Ке-в|5') = (Ке-В^А*, (65)
получим
адя, п, I, Ь, Я) = иА+в(Я) + 2итг1и^в(11)\Уь(е = 0)| х
ктах
II
к2йкк'2йк']1{кЯ)^ь{к'Я)
(Ке-в)
ьь
1 + г(Ке-в)ьь.
х
тпМ
Уьп ( д
Уьм ( ~ Р
(66)
где £к = Е—к2/2Мс, Ь — начальный орбитальный момент сталкивающихся частиц А** и В. Шаровую функцию У1т(И/Я) в (66), которая входит также в матрицу д33'(Я,£), следует заменить на величину У5т(0 = 0). Это связано с тем, что при I <С и^2 движение электрона вблизи атома В описывается плоской волной с импульсом, направленным вдоль вектора К.
Поскольку вблизи атома В ридберговский электрон ведет себя как свободная частица, диагональные элементы Ке-д-матрицы в (66), зависящие от полной энергии системы Е и кинетической энергии падающего ее и рассеянного е'е электронов, можно заменить на соответствующие элементы К®в-матрицы рассеяния свободного электрона, взятые при
„ 1 к2 £е = Е+Я~Ш.
(67)
В этом случае выражение (66), определенное для заданной начальной энергии Ек, принимает следующий вид:
УорЬ = иА+в + Д - г-,
(68)
где сдвиг и уширение ионного терма равны
'(0)
А(Ек, п, I, Ц Я) = ^¿¿(".Д) ^(Е,, Я)
1
(69)
Г (Ек,п,1,Ь,Я) =
2тг
1 +
К?1(п,Я)
2 К)
(70)
и зависят
от начального возбуждения атома А**(п, I), орбитального момен-
Е, эВ
Рис. 1. Зависимость ширины оптического потенциала Г от полной энергии Е для системы ]Уа**(10«) + Не.
та Ь и расстояния Я между частицами А+ и В. Выражение для фактора 5,2 в (69) и (70) определяется знаком полной энергии системы Е и записывается в виде
8.г(Ек,Я) =
для Е < 0 и
'"тах
64тгС/Л+в(Д)У|о(0 = О)Уга(0 = 0) У з1{кЯ)д83(Я,ек)кЧк
(71)
.й) = 647г[Ъ+в(Д)уЗ(0 = 0)
х
I %(кЯ)Кед33(Я,ек)к2<1к + г ^ ]\(кЯ)1т двв{Я, ек)кЧк о о
(72)
для Е > 0.
В качестве иллюстрации рассмотрим медленные столкновения ридбер-говского атома натрия с атомом гелия. На Рис.1 изображена зависимость ширины ионного терма Г от полной энергии системы Е, построенная для следующих значений начальных данных:
Ек = 0.19 эВ, п0 = 10, 1 = Ь = 0,
(73)
(7е(, СМ2
ек, эВ
Рис. 2. Зависимость полного сечения упругого рассеяния Na**(10s) + Не от кинетической энергии сталкивающихся частиц. Тонкая линия результаты расчета с использованием оптического потенциала (08); гладкая линия сечение упругого Ма+ — Не рассеяния, рассчитанная на поляризационной части потенциала.
Как и следовало ожидать, осцилляции сходятся к пределу При иеболь-
(г)
ших значениях п в окрестности ридберговских Еп -рсзонапсов, положение которых определяется полюсами матрицы де8(Е)7 наблюдается расщепление соответствующих пиков. С ростом п картина заметно усложняется.
На Рис.2 приведена зависимость сечения упругого рассеяния Ма**(пв)+ Не от кинетической энергии Ек (для фиксированного главного квантового числа По = 10). Полное сечение рассеяния сге;(^.), изображенное на Рис.2, является сильно осциллирующей функцией от кинетической энергии.
Мы ограничились анализом простейшей ситуации, когда атом В является бесструктурной частицей. Для атомов инертных газов это приближение вполне оправдано. Однако, при рассмотрении более сложных атомов необходимо учитывать многоканальный характер движения электрона, связанный с возбуждением ионного остова, а также возможность его захвата с образованием ионной пары.
В третьей главе диссертации рассматриваются высоковозбужденные состояния атомов и молекул, которые образуются в нижних Е- и Д-слоях ионосферы Земли в периоды повышения солнечной активности и возникновения сильных геомагнитных возмущений. Это связано с тем, что в Е- и Д-слоях атмосферы, расположенных на высоте 50-110 км от поверхности Земли, формируется свсрхвысокочастотное (СВЧ) излучение, интенсивность которого многократно превышает типичные уровни всплесков Солнца. Анализ различных возможностей генерации обнаруженного СВЧ из-
Е
Е > О
Я
Е = 0
2
Е.
Рис. 3. Энергетическая схема положения горлышка стока.
лучения установил, что оно обусловлено радиационными переходами между ридбсрговскими состояниями атомов и молекул, которые возбуждаются потоками электронов, выбрасываемых из ионосферы.
Ридбсрговскими называются такие высоковозбужденные состояния атомов и молекул, которые расположены вблизи границы ионизации и характеризуются наличием бесконечной последовательности уровней энергии, сходящихся к порогу ионизации. Ридбсрговские атомы и молекулы обладают одним возбужденным слабосвязанным электроном, состояние которого характеризуется энергией уровня с угловым моментом I относительно ионного остова. Энергии уровней с большими угловыми моментами не зависят от I, т.е. являются орбитально вырожденными. Статистически наиболее устойчивыми оказываются именно эти состояния, в которых основное время электрон проводит на больших расстояниях от ионного остова. Процесс, приводящий к формированию таких состояний, называется /-перемешиванием. В верхней атмосфере он протекает быстро и является необратимым, т.е. практически все ридбсрговские частицы оказываются вырожденными. В результате различия между высоковозбужденными атомами пропадают и спектр излучения должен зависеть только от нейтральных компонент верхней атмосферы.
Как показано на Рис.3 ридбсрговские состояния заселяются выше некоторой энергии связи Е„, формирующейся в условиях двухтемпературпой квазистациопарпой неравновесной плазмы, которая создается в верхней атмосфере под действием потоков электронов, высыпающих из ионосферы
при сильных геомагнитных возмущениях. В этих условиях распределение заселенностей по уровням Еп, расположенным вблизи границы ионизации, является практически равновесным и характеризуется температурой, близкой к температуре свободных электронов Те. Для больших энергий связи Еп ~ Е* на энергетическом интервале АЕ Е* равновесное распределение сильно нарушается. Этот интервал называется "горлышком стока"или "узким местом рекомбинационного потока". Выше горлышка стока (т.е. при Еп ^ Е*) преобладают столкновительные переходы между связанными состояниями и континуумом. Ниже горлышка доминируют излуча-тельные переходы, приводящие к равновесному заселению низколежащих состояний с температурой среды Та.
Прохождение электронов через горлышко стока шириной АЕ «С Е, по энергетической шкале является наиболее медленной стадией процесса, которая определяет кинетические особенности неравновесной плазмы. Положение горлышка стока находится из условия минимума константы скорости тушения за счет переходов в нижележащие состояния и определяется следующим образом:
Я, = Щ я 27.21
п{
ч 21 !/4
пе /11600Ч2
4.5 • 1013 \ Те )
(74)
Здесь концентрация электронов пе выражена в см-3, температура электронов Те — в Кельвинах. Вблизи ионизационного континуума распределение ридберговских частиц по состояниям с большими главными квантовыми числами п можно определять по формуле Больцмана:
¡гр \ № ч ехр [-(Ек - Еп)/Те}
тпп(Те) = тпк{Те)-, (75)
9к
где дп— стерень вырождения или статистический вес данного уровня, т.е. число состояний с энергией Еп. С этой целью достаточно выделить уровни, расположенные выше горлышка стока и определить эффективную заселенность шп(Те), которая для рекомбинационной плазмы будет связана с концентрацией свободных электронов пе(Те) соотношением:
тп(Те) = Рп(Те)п2е(Те)Кп(Те),
(76>
= \ 2,л;
где дп = 2п2 — статистический вес ридберговского состояния, множитель
Рп{Те) характеризует степень отклонения от равновесной плазмы и определяется потоком электронов. В случае Р„(Те) = 1 выражение (76) переходит в формулу Саха. Полная статистическая сумма положительных ионов Е,-в (76) определяется как
где б?* — статистический вес молекулярного иона А+ для заданного полного электронного момента Ь и спина Б, V — колебательное квантовое число, — частота колебаний, — вращательная постоянная. Суммирование в (77) выполняется по всем возможным значениям полного момента
Процесс ¿-перемешивания является важнейшим актом в ходе формирования резонансных квантовых свойств возмущенного потоком электронов .О-слоя атмосферы. Одновременно с /-перемешиванием происходит расщепление уровней энергии орбитально вырожденных состояний за счет взаимодействия с нейтральными молекулами азота и кислорода, что приводит к образованию дискретных уровней промежуточных ридберговских квазимолекул А**АТ2 и А**02 в классически разрешенной области движения слабосвязанного электрона. Радиационные переходы между ними, происходящие без изменения главного квантового числа (Дп = 0), соответствуют сантиметровому и дециметровому диапазону СВЧ излучения.
Фотоионизационная плазма создается под действием широкополосного излучения, приходящего от солнечной вспышки, в течение 20-30 минут. Этот процесс обусловлен многоквантовым возбуждением электронных состояний атомов и молекул, в которых спиновые запреты для соответствующих радиационных переходов снимаются вследствие взаимодействия с молекулами среды М. Отличие засоленностей высоколежащих ридберговских состояний от случая неравновесной рекомбинационной плазмы заключается в том, что они дополнительно обедняются за счет фотоионизации. Для небольших п опустошение ридберговских и низколежащих возбужденных состояний происходит за счет процессов предиссоциации, включая неадиабатические переходы через промежуточные валентные конфигурации, и резонансной (нерезонансной) передачи внутренней энергии в результате столкновительных процессов с последующей термолизацией среды. На это указывает рост ее температуры с увеличением высоты в интервале 40-60 км. Следует отметить, что для этих состояний процесс /-перемешивания является сильно подавленным и влияние среды заметно уменьшается. Последнее особенно важно для формирования частотного профиля спектра ИК-излучения (Дп = 1), изображенного на Рис.4. Прямым указанием на
Ел + Та
ЗгТа
(77)
0.20
О. 0.15 U
0.00
0 500 1000 1500 2000
волновое число, см-1
Измерения спутника FIRST 7 июня 2005 г., 14:25 UTC
Рис. 4. Спектр ИК-излучения, записанный спутником FIRST 7 июня 2005 года, это обстоятельство является наличие характерного спада интенсивности с увеличением частоты излучения от положения первого максимума, отвечающего длине волны около 20 мкм, до минимума вблизи 15 мкм, что связано с процессом /-перемешивания. Заметим также, что ИК-излучение, образующееся на высотах выше 120 км, принципиально не может обладать такой рельефной структурой. Таким образом, наблюдаемое ИК-излучение может быть использовано для детектирования ридберговских состояний.
В случае возбужденных состояний, возмущенных нейтральными атомами или молекулами газовой среды, образуется ридберговская квазимолекула, которая имеет два (а в достаточно плотных средах и более) разделенных центра рассеяния слабосвязанного электрона: ионный остов и возмущающая частица М. Ионным остовом возмущаются состояния с малыми моментами / ^ /*. Поле возмущающей частицы М в свою очередь оказывает влияние только на такие суперпозиции состояний кулоновского центра, которые обладают малыми значениями углового момента L относительно М. Эти две разноцентровые группы термов будем сокращенно называть ридбергоескими /- и вырожденными ¿-термами.
В орбитально вырожденных состояниях (при / > /*) изолированная
ридберговская молекула XY** отличается от атома А** тем, что обладает ядерной подсистемой, которая совершает колебательное и вращательное движения. Так как температура среды Та в рассматриваемых нами условиях порядка тепловой, переходами между различными сериями вырожденных состояний молекулы XY** вследствие неадиабатической связи электронного и ядерного движений можно пренебречь и ограничиться рассмотрением только одной серии основного состояния ионного остова. В этом приближении электронные спектры высоковозбужденных состояний атома А** и молекулы XY** полностью совпадают.
Таким образом, задача сводится к анализу особенностей ППЭ орбиталь-но вырожденных состояний квазимолекулы А**М, составленной из атома А и частицы М. Отметим, что взаимодействие электрона с ионным остовом и возмущающей частицей М является сильным и не может быть учтено по теории возмущений, что является принципиальной особенностью этой задачи.
На Рис.5 изображены зависящие от главного квантового числа п частотные линии излучения ридберговских квазимолекул А**N2 и А**02. Линии излучения начинаются от положения горлышка стока п* = 19, которое отвечает сильному геомагнитному возмущению с плотностью электронов пе = 107см-3 и температурой Те = 3500 К. В диапазоне частот до 10 ГГц распределение зависящих от п линий излучения, отвечающих значениям L = 0 — 3, для каждой из молекул содержит по четыре серии линий переходов L —> L\ которые сходятся с ростом L' к пределу перехода L —>■ п. На Рис. 5(а) для молекул N2 и О2 эти линии обозначены соответственно N^jj и Olu- Сдвиг пределов серий Nm и Ощ по частотам относительно друг друга приводит к существенной неоднородности результирующей картины спектра излучения. Здесь выделяется три спектральные области изменения частот, в которых при небольших значениях п линии переходов отсутствуют. К первой относится диапазон, расположенный между близколежащими линиями N23, N12 и линией Оз„, частоты переходов для которых равны = 1-17 ГГц и oj(N23) = 1.706 ГГц, т.е. ширина безызлучательного диапазона составляет Aoj = 0.59 ГГц.
Второй диапазон соответствует разности частот переходов üj(Nin) = 4.31 ГГц и w(Oin) = 6.094 ГГц с шириной Aoj = 2.63 ГГц. Третий диапазон, изображенный на Рис.5, определяется разностью частот переходов u{Oln) = 7.27 ГГц и üj{N01) = 57.10 ГГц, равной До; = 49.83 ГГц. Ширины и положения этих диапазонов по шкале п связаны с величиной горлышка стока, обозначенной штриховыми линиями и зависящей от уровня геомагнитного возмущения, который отвечает определенным значениям
п. = 10' см1 п. = 10г' см * п, = 10* см1
Т,= 3500 К Т.= 2000 К Т,= 2000 К
Главное квантовое число, п
а)
п„= 10й см'3 п„ = 105 см'3
Главное квантовое число, п
б)
Рис. 5. Зависимость частоты излучения высоковозбужденных квазимолекул А**N2 и Л**С?2 от главного квантового числа п. Жирными линиями обозначены границы диапазонов ослабления интенсивности излучения. Вертикальными штриховыми линиями обозначены положения горлышка стока для различных значений температуры Те и кон-центраци пе свободных электронов: (а) полоса частот от 0.5 до 7.5 ГГц. (б) полоса частот от 7 до 100 ГГц.
плотности и температуры электронов.
Интенсивность некогерентного СВЧ излучения с частотой V в единице объема записывается следующим образом:
= Е а(г)юп(() Е Ре//(п)тпДг)Щ^Н5), (78)
¿=1,2 пв
где а(г) — весовые множители молекул N2 и О2, ~~ соответствующие парциальные вероятности столкновительиых процессов. Индекс 5 в формуле (78) принимает значения I/ > Ь{р) или п в соответствии с порядком пересечений с линиями излучения Ь(у) для заданной частоты излучения и (см. Рис.5).
Полная мощность СВЧ излучения ридберговских частиц с заданной частотой V на единичную площадь поверхности Земли дается выражением:
н/¡М
1ш{и,пе,Те) = К I Шм{р,пе{К),Те{К))ак, (79)
о
где зависимость пе(К) берется из эксперимента, К — соответствующий коэффициент ослабления. Фактор исравновеспости плазмы Рец учитывает поток высыпающих из ионосферы электронов и тушение ридберговских частиц. Оптическая толщина излучающего слоя отсчитанная от его верхней границы Но = 110 км, определяется предельным значением плотности среды р*а{Ня{и)), начиная с которой величина \¥м(и,пе,Те) практически обращается в пуль.
Рассмотрим диапазон частот 0.8 - 1.8 ГГц. На Рис.6 приведена зависимость мощности потока некогерентного СВЧ излучения при концентрации электронов па верхней границе слоя пе(ра) = 104 см~,! для различных значений температуры электронов Те. Температура Те в данном случае изменяется от 1000 К до 2000 К и входит как средняя величина при интегрировании в (79) по интервалу 90-110 км. Видно, что при температуре Те = 1000 К мощность потока излучения 1ш(и,пе,Те) имеет локальный минимум на частоте и = 1.6 ГГц, который соответствует столкновительпому провалу в вероятности процесса ¿-перемешивания вследствие неадиабатической связи с вращением в ридберговской молекуле азота. На частоте V = 1.44 ГГц формируется точка перетяжки. Причем ее положение по частоте практически не зависит от температуры Те, а мощность СВЧ излучения на частоте и = 1.6 ГГц с ростом температуры Те заметно возрастает.
На Рис.7 представлена зависимость мощности потока излучения
Рис. 6. Зависимость мощности потока излучения от частоты V для различных температур электронов Те при заданной концентрации электронов пе = 104см_3. Черные треугольники — Те = 1000 К; белые кружки — Те = 1200 К; черные квадраты — Те = 1500 К; черные кружки — Те = 2000 К.
I,o,, 10'28 Вт/см 2
Рис. 7. Зависимость мощности потока излучения Itot от частоты и для различных концентраций электронов пе при заданной температуре электронов Те = 2000 К. Белые кружки — пе = 0.6-104 см"3 ; черные треугольники — пе = 0.8-104 см-3; черные кружки — пе = 1.0 • 104 см-3; черные квадраты — пе = 1.2 • 104 см-3.
от частоты для различных концентраций электронов пе с фиксированной температурой электронов, равной 2000 К. На частоте v = 1.6 ГГц наблюдаются изломы, и мощность СВЧ излучения с увеличением пе растет. На частоте v = 1.44 ГГц мощность излучения также возрастает. Причем для каждого значения пе с увеличением температуры Те точка перетяжки для семейства кривых на данной частоте смещается вертикально вверх. Это означает, что положение точки перетяжки на вертикальной оси зависит только от концентрации пе, т.е. здесь имеет место однопараметрическая зависимость.
Обратимся теперь к диапазону частот 4.0-6.0 ГГц, где происходит необычное поведение мощности потока Itotiv). В отличие от диапазона 0.8-1.8 ГГц с минимумом vi = 1.57 ГГц поведение мощности потока здесь обусловлено радиационными переходами и является более устойчивым к изменениям параметров плазмы. При небольших концентрациях электронов пе в спокойных геомагнитных условиях мощность потока Itot(v) для частот v > 4.0 ГГц мала, а с увеличением пе заметно возрастает. Эволюция этого роста для различных температур электронов Те продемонстрирована на Рис.8. Видно, что в зависимости от Те графики для мощности потока излучения претерпевают инверсию, и соответствующие кривые изменяются с вогнутых на выпуклые (см. Рис.8 (а-в)), сходящиеся в точке, расположенной вблизи частоты г/2 = 5.0 ГГц. Причем прямое наложение этих графиков показывает, что они практически симбатны (т.е. однотипны). Другими словами, наблюдается эффект стабилизации в общей картине зависимости мощности потока излучения от частоты.
Такую форму графиков в радиофизике называют "бутылочным горлышком" ("bottle neck"). Свойства этого горлышка напрямую связаны со вторым диапазоном уменьшения интенсивности некогерентного СВЧ излучения (4.31 - 6.09 ГГц). Причем в широкой области изменения Те положение узкого места вплоть до концентраций пе = 2 • 104 см оказывается фиксированным. Это означает, что частотный профиль ПК спектра, формирующийся за счет ридберговских состояний, излучающих в D-слое на частотах, близких к щ = 5 ГГц, будет зависеть только от концентрации электронов пе.
Заметим, что в интервале 5 • 103 < пе < 2 • 104 см-3, где все кривые симбатны, при определенном выборе изменения температуры Те формируется практически линейная зависимость Ziof(г/2) от пе. Это связано с тем, что квадратичный рост мощности СВЧ-излучения с увеличением концентрации электронов компенсируется ее падением с ростом температуры Те, а также зависимостью положения горлышка стока (74) от температуры электронов. Положение бутылочного горлышка на частоте ъ>-2 зависит от
Ito,, 10 Вт/см
1,5*10*
/X.
4,0 4,S 5,0 5,5
a)
v, ГГ1
Itot, 10 Вт/см 5x10*-
4x10*-
3x10*-
2x10* •
1x10*-
4,0 4,5 5,0 5,5
6)
6,0
V, ГГц
4,0 4,5 5.0 5,5
в)
Itoi, 10' Вт/см' 2,0x1 o'-i
1,5x10'
1,0x107 5,0*10*
4.0 4,5
5,0 5.5
r)
v, ГГц
I ta, 10*" Вт/см 2
1,5x10'-i
Д)
4,0 4,5 5,0 5.5
V, ГГц
Рис. 8. Зависимость мощности потока излучения Itot от частоты v в диапазоне 4.0-6.0 ГГц для различных концентраций пе и температур Те электронов: а) пе = 5 • 103 см-3, квадраты — Те = 1000 К, кружки — Те = 1500 К; 6) пе = 104 см"3, квадраты — Те = 1500 К, кружки — Те= 2000 К; в) пе = 2 • 104 см"3, квадраты — Те = 2000 К, кружки — Те = 2500 К, треугольники -Те = 3000 К; г) тге = 5 • 104 см~3, квадраты — Те = 2500 К, кружки — Те = 3000 К; д) пе = 105 см"3, квадраты — Те = 2500 К, кружки — Те = 3000 К.
величины пе и не зависит от температуры электронов Те (см. Рис.8(а-в)).
При большей концентрации 2-104 < пе < 5-104 см~3 наблюдается переходный процесс (см. Рис.8(г)). Далее, начиная со значения пе = 105 см-3, восстанавливается симбатность кривых зависимости Itot{v) на фоне возрастающей зависимости графиков от частоты (см. Рис.8(д)).
Так как ширина бутылочного горлышка зависит от концентрации электронов пе, можно полагать, что мощность потока излучения в диапазоне частот 4.0 — 6.0 ГГц, за исключением переходной области значений пе, зависит только от вариации АТе. При больших концентрациях электронов появляется нелинейная зависимость мощности потока излучения от температуры электронов. Это связано с проявлением совокупности физических факторов: эффективности процесса ¿-перемешивания, потока электронов и зависимости Itotiv) от их температуры.
В четвертой главе рассмотрена генерация СВЧ излучения в атмосфере Земли, наведенного точечным импульсным гамма-источником с длительностью импульса Дт 100 не, испускающим за один импульс 1022 фотонов с энергией = 1 — 3 МэВ. Вследствие взаимодействия испускаемых 7-квантов с атмосферой образуются высоковозбужденные атомы и молекулы, которые излучают в длиноволновом диапазоне. Проанализированы кинетика и механизмы их заселения. Кроме того, определены амплитуды напряженности электрической компоненты поля и потоки излучения в заданном СВЧ диапазоне в зависимости от расположения источника и приемника.
Известно, что в области энергии Еу ^ 10 МэВ основными процессами при взаимодействии 7-кваитов с веществом, приводящими к образованию заряженных и высоковозбужденных частиц являются: фотоэффект, эффект Комптона и рождение электрон-позитронных пар. В фотоэффекте, который наиболее существен при энергиях Е1 < 1 МэВ, фотон поглощается атомом. В результате рождается электрон-ионная пара, причем дополнительно могут также освобождаться несколько свободных электронов. Эффект Комптона становится доминирующим в области энергий E~t > 1МэВ. Рождение электрон-позитронных пар эффективно в области энергий Е-у ^ 3 МэВ.
Для решения поставленной задачи необходимо оценить число электронов, образованных 7-излучением при прохождении фотонов от точечного источника через слой воздуха толщиной L в заданный элемент на плоской поверхности. Вблизи поверхности Земли при нормальном давлении длина пробега 7-кванта составляет ¿7 ~ 300 м. В этих условиях поток
7-излучения точечного источника через сферу радиуса г убывает как:
где N-[(0) — полное число испущенных источником 7-квантов за один импульс, а7(р) — полный линейный коэффициент ослабления, зависящий от плотности среды р, величина /3 есть безразмерный фактор накопления.
Будем считать движение 7-квантов прямолинейным, полагая, что распределение образующихся электронов изотропно. Тогда для концентрации ионизованных электронов получим
Пе(р) =
3х7(р)дг7(0)/3
тах
2ехр[-а7(р)г]
¿г,
(81)
где х7 — ионизирующая способность одного 7-кванта. Величина Ьтах есть минимальный верхний предел интегрирования, начиная с которого интеграл в (81) не меняется. По физическому смыслу эта величина представляет собой зависящую от плотности р эффективную длину пробега 7-квантов от источника, т.е. Ьтах = Ь-у(р). В Таблице 1 представлены результаты расчета длин пробега гамма-квантов и концентраций ионизованных электронов в зависимости от высоты источника 7-квантов.
Таблица 1. Зависимости плотности ра и температуры Та атмосферы, длин пробега Ь1
Ь, км Ра см 3 Та, К М пе, см 3
0 2.7 ■1019 300 300 4.0 •1011
10 6.1 •1018 220 690 7.5 •Ю10
20 1.6 • ю18 230 1260 2.3 • ю10
30 5.4 • ю17 250 1360 1.9 • ю10
40 2.6 • ю17 280 1420 1.8 •ю10
50 1.4 • ю17 310 1530 1.5 • ю10
60 5.1 ■ ю16 310 1730 1.2 • ю10
70 1.5 • ю15 230 1550 1.5 • ю10
80 2.8 • ю14 160 1540 1.5 • ю10
90 1.0 • ю14 180 1540 1.5 • ю10
100 3.1 • ю13 230 1540 1.5 • ю10
Видно, что от поверхности Земли до высоты 20 км наблюдается сильное увеличение длины пробега 7-квантов и рост концентрации ионизованных электронов, что объясняется одновременным уменьшением плотности и температуры атмосферы с ростом высоты. На интервале 20-50 км увеличение Ьп и падение пс демонстрируют плавную зависимость от высоты вследствие роста атмосферной температуры Та. Пик в и коррелиро-
ванный с ним провал в пе{Ъ) вблизи высоты /г = 60 км отвечает максимуму температуры атмосферы. Последующее значительное уменьшение плотности выводит эти зависимости на плато.
Согласно результатам Таблицы 1 кинетические энергии электронов, ионизованных 7-квантом на различных высотах к, в среднем порядка 10 эВ. Поэтому каждый такой электрон может возбудить не более одной ридбер-говской частицы, так как его энергия ненамного превышает ее потенциал ионизации. Наиболее эффективным процессом образования ридбергов-ских атомов и молекул в этом случае является околопороговое возбуждение электронами с энергиями е, близкими к потенциалу ионизации I, т.е. когда энергия рассеянного электрона е' I.
Ударный механизм возбуждения состояний с большими значениями I подробно обсуждался нами выше. С учетом дополнительного вклада процесса /-перемешивания можно утверждать, что при ударном заселении возбужденных ридберговских состояний атомов и молекул должны полностью доминировать орбитально вырожденные состояния.
Радиационный механизм возбуждения ридберговских состояний в рассматриваемом случае может происходить только под действием 7-квантов. Для заданных энергий Е1 к, 1 — 3 МэВ механизм фотоионизации атомов и молекул вносит малый вклад в процесс образования первичных электронов. По этой причине радиационным механизмом можно пренебречь.
Рождение Ие я» 4.5 • 1025 электронов, ионизованных 7-квантами за время импульса порядка Ю-7 с, приводит к формированию слабоионизованной низкотемпературной плазмы с локальным возмущением равновесия среды внутри Ь1 сферы. За времена порядка Ю-10 с здесь установится неравновесное стационарное распределение по энергиям возбуждения, которое зависит от температуры среды Та и температуры электронов Те. Поскольку концентрации возникших электронов малы по сравнению с плотностью атмосферы, то изменения ее температуры фактически не наблюдается. При этом спонтанное излучение произойдет значительно раньше, чем установится термодинамическое равновесие.
Концентрацию свободных электронов п*(Ге) и заселенности возбужденных состояний ридберговских атомов и молекул тпп(Те) можно определить из уравнений баланса с учетом процессов рекомбинации, ионизации и излу-
чения. В этих условиях формируется два локальных больцмановских распределения по энергиям дискретных состояний атомов и молекул. Первое образуется за счет быстрого обмена энергией между связанными и свободными электронами для высоковозбужденных электронных состояний, расположенных выше горлышка стока. Второе обусловлено столкновениями между тяжелыми частицами среды и относится к низколежащим электронным состояниям. Заселенность уровней, расположенных между этими областями спектра мала, так как они эффективно тушатся за счет радиационных переходов в видимом и ближнем УФ-диапазонах.
Вблизи ионизационного континуума с ростом главного квантового числа тг интенсивность радиационных процессов заметно уменьшается, в то время как интенсивность столкновительных процессов существенно возрастает. Соответствующие сечения переходов а ~ п4. Это позволяет разделить интервал энергии на две области: выше горлышка стока Е > Е„, где доминируют процессы столкновения, и ниже горлышка стока Е < Е* где возбуждения происходят за счет радиационных переходов. Состояния, относящиеся к первой области спектра, являются практически равновесными между собой и ионизационным континуумом.
В общем случае концентрации электронов п*е{Те) и ионов п*(Та) не совпадают по величине. Уравнение для стационарной концентрации электронов имеет следующий вид:
Суммирование проводится по всем дискретным ¿-состояниям с заселенностью тпк; ги^е и ъиек — вероятности переходов из дискретного ¿-состояния в ионизационный континуум и обратно. Величина = пе/Ат представляет собой стационарный поток ионизованных 7-квантами электронов, Се(Т,,) описывает уход электронов из объема 1/7-сферы за счет пространственной диффузии и составляет порядка 1 /то, где характерное время тд ~ Ю-9 с.
В области спектра, расположенной выше горлышка стока стационарный ионизационный потоку « 0, поэтому для концентрации электронов п*е(Те) внутри сферы имеем
<И
= Зе + Ре~ П*еве И 0.
(82)
Ионизационный поток в (82) равен
(83)
к
Пе Т0 Пе
- ~ 71»- ~ -.
ДтСе Дт 10
Регистрируемый на приемнике электромагнитного излучения частотный диапазон составляет Аи = 0.8 — 1.0 ГГц. Для переходов с изменением главного квантового числа Ап = 1 этот диапазон отвечает главным квантовым числам п = 150 — 200. Для процессов, протекающих без изменения главного квантового числа (Ап = 0), существенны два типа дипольно-разрешенных переходов:
1) между ридберговскими состояниями с небольшими угловыми моментами электронов I <1* = 3;
2) между орбитально вырожденными (ковалентными) состояниями с большими моментами 3 ^ I ^ п — 1.
Первый тип отвечает случаю, когда влияние среды мало. Для основных атмосферных атомов и молекул в среднем эти переходы отвечают интервалу главных квантовых чисел п = 100 — 200.
Переходы второго типа обусловлены взаимодействием с окружающей средой и наиболее эффективны, когда нейтральные частицы попадают в область классически разрешенного движения слабосвязанного электрона радиуса 2п2. Ионным остовом возмущаются состояния с малыми угловыми моментами I относительно кулоновского центра, т.е. I = 0,1,2. Полное число вырожденных состояний равно п2 — I2, где I2 — число ридберговских состояний. Разность энергий начального г- и конечного /-состояний в этом случае оценивается по формуле:
{Щ) (85)
■кпл 7ггг4
где а — длина рассеяния электрона на возмущающей частице среды, (ре) л; 1 /п — характерный импульс электрона, движущегося в кулоновском потенциале. Интенсивность перехода здесь не зависит от электронной структуры возбужденных частиц и может быть представлена в виде
и шШ\2 « °4к"г|/)'а (86)
<Ш бтгс3 К 1 |/Л бтгЗДб ' ^
где с — скорость света, г — координата выделенного электрона.
Интенсивность излучения за счет переходов между вырожденными ридберговскими состояниями определяется как
^ 47ГI? х ^
1 = 2^ ^i^fIif = —з (87)
где — разность энергий начального г и конечного / состояний (см.
(85)); д{ и — кратности вырождения этих состояний; тг — заселенность начального состояния; A¡f — коэффициенты Эйнштейна, связанные с г-компонентой вектора дипольного момента перехода (г\г\/) соотношением:
Мтг^Р2
Д/ = (88)
Тогда для амплитуды электрической компоненты поля излучения ридбер-говских состояний запишем:
-^■ ^ £ т
где £о и Цо — электрическая и магнитная постоянные, е — заряд электрона.
Полная интенсивность регистрируемого излучения зависит от взаимного расположения источника и приемника. Если они пространственно разнесены (обозначим это как случай (а)), то поток излучения на единицу поверхности приемника ]г для условия Н > Н + Ь7 (Н — высота приемника над поверхностью Земли) равен
Здесь I определяется выражением (87), ДГ2Г — телесный угол от источника на приемник, 5 — площадь приемника, Я = Н — к — расстояние от источника до приемника. Умножением на 1/2 мы учитываем только те фотоны, которые движутся в направлении приемника. Если источник испускает Щ(0) = 1022 гамма-квантов за один импульс, то поток излучения
7\
(а)
на высоте приемника Н = 20000 км составляет $ ' = 6.4 • 10 18 Вт/м2.
В случае (б), когда приемник расположен на поверхности Земли непосредственно под источником, возможны две ситуации. Первая соответствует условию К < Тогда для потока на приемнике имеем
,(б>
31 Г
= 21? У ю(гУ2(1г> (91)
7
о
где "ш{г) = 1/2г2 — вероятность прихода испущенных фотонов на единицу поверхности приемника.
Вторая ситуация отвечает условию К 1>7, когда интенсивность опи-
Таблица 2. Зависимости интенсивности I, амплитуды электрической компоненты по-
к, км I, Вт Етах, В/м > Вт/м2 Зг ) Вт/м2
0 1.7 10"2 6.2 10"1 6.4 ю-18 4.3 •10"8
10 1.0 10"2 4.8 10"1 3.8 10"18 1.5 Ю-11
20 5.3 10"3 3.5 ю-1 2.0 Ю-18 2.0 10"12
30 4.6 Ю-3 3.2 ю-1 1.7 10"18 7.6 10"13
40 4.0 10~3 3.0 10"1 1.5 10-18 3.6 Ю-13
50 3.2 10~3 2.7 ю-1 1.2 10"18 1.9 Ю-13
60 2.8 10"3 2.5 10"1 1.1 Ю-18 1.2 Ю-13
70 4.3 Ю-4 3.1 10"1 1.6 10"18 1.3 Ю-13
80 6.3 Ю-4 3.7 ю-1 2.3 10"18 1.5 10"13
90 5.6 Ю-4 3.5 ю-1 2.1 10-18 1.0 Ю-13
100 4.3 Ю-4 3.1 ю-1 1.6 10"18 6.5 10"14
сывается формулой (90), в которой нужно положить Н = 0 и К = /¿, т.е.
# - (92)
В Таблице 2 представлены интенсивность излучения светящегося шара (87) и амплитуда электрической компоненты поля (89), а также даны рассчитанные по формулам (90)-(92) зависимости потоков ридберговского излучения на приемнике от высоты источника /г. Отметим для сравнения, что фоновое излучение верхней атмосферы, приходящее сверху, порядка
« 10~24 — Ю-23 Вт/м2, т.е. существенно меньше по величине регистрируемого на приемнике излучения.
В пятой главе диссертации обсуждается возможность лазерного управления реакцией диссоциативной рекомбинации (ДР) медленных электронов и положительно заряженных молекулярных ионов. Реакция
е" + ХУ+ —>Х + У* (93)
является важнейшим процессом в атмосфере, отвечающим за кинетику исчезновения зарядов и играющим важную роль на высотах, превышающих 100 км от поверхности Земли.
Среди многообразия положительных ионов, принимающих участие в
реакции (93), особое место занимают ионы кислорода О^ ■ Максимум их концентрации соответствует высоте порядка 200 км, лежащей в ^-области ионосферы. Возбужденные продукты этой реакции, 0(15) и 0{1Б), вносят заметный вклад в свечение атмосферы.
Лазерное управление элементарными процессами, включая реакцию (93), является одной из фундаментальных проблем химической физики. Хорошо известно, что внешнее поле не приводит к дипольно-разрешенным переходам в изолированном ионе ХУ+ и не влияет на продукты этой реакции. Оно воздействует на состояния ХУ**, образующиеся на промежуточном этапе процесса, так как полевое воздействие на систему "электрон + мишень" возможно только при наличии ее внутренней структуры. Наиболее эффективно это реализуется, когда падающие электроны не взаимодействуют с полем и обладают определенной энергией Ее, т.е. при выполнении условия (Н = те = е = 1)
ре/а;/2 <1, Ре = л/2Ее, (94)
где — импульс падающего электрона, / и ш/ — амплитуда напряженности и частота внешнего электромагнитного поля. Дополнительным ограничением на величину напряженности поля является требование
/Я « 1, (95)
что позволяет при описании реакции (93) ограничиться учетом трех типов радиационных переходов. К ним относятся безызлучательный (к = 0) и два вынужденных (к = ±1) перехода, отвечающих излучению или поглощению фотона. Индекс к обозначает изменение числа фотонов в системе, Б — дипольный момент перехода.
Расчеты, представленные в диссертации, показывают, что сечение реакции
' 0(3Р) + 0(3Р) + 6.96 эВ, 0{3Р) + 0(1£>) + 4.99 эВ, О^Б) + 0(1£>) + 3.03 эВ, (96) 0(3Р) + 0(х5) + 2.77 эВ, ^ 0(1£') + 0(15) + 0.80 эВ
с индуцированными радиационными переходами в видимой области спектра при определенном выборе параметров может возрастать на несколько порядков. Такое поведение сечения указывает на прямую возможность лазерного стимулирования этой реакции.
е~ + 0^(Х2Пд,у = 0) —<
Е, эВ
Рис. 9. Схематическое изображение диабатических потенциальных кривых диссоциативных термов, которые ответственны за реакцию (96) в отсутствие внешнего поля.
В отсутствие внешнего поля основной вклад в реакцию (96) вносят только такие диссоциативные состояния, потенциальные кривые которых лежат вблизи минимума ионного потенциала. Исходя из расчитанных величин матричных элементов конфигурационной связи, было выделено шесть основных диссоциативных термов (см. Рис. 9). При наличии поля в реакцию е~ 4- 0% + kuJf подключаются дополнительные индуцированные полем состояния, которые заселяются за счет поглощения или излучения фотонов.
Считая внешнее электромагнитное поле классическим (полагая, что среднее число фотонов велико, N0 1), взаимодействие системы е~+ХУ+ с этим полем будем описывать периодической по времени зависимостью:
и/ = совш^, У/ = —. (97)
Тогда при условии (94) исследование может быть выполнено в рамках стационарной теории многоканального квантового дефекта (МКД), в которой наличие поля учитывается введением "квазиэнергетических состояний" и связанных с ними новых каналов движения.
Нулевой гамильтониан Но определяется на кулоновском базисе таким образом, чтобы все взаимодействия в изолированной диабатической конфигурации X + У* учитывались точно. Тогда в полном гамильтониане системы
Н = Н0 + V (98)
оператор V является суммой
V = Vе + V', (99)
где оператор электростатического взаимодействия
Vе = УпС + Ус/ (100)
включает в себя некулоновскую часть взаимодействия с ионным остовом УпС и конфигурационное взаимодействие Vе7, связывающее конфигурации е~ +ХУ+ и ХУ**, которым мы присвоим индексы д и /3, соответственно.
Теория МКД дает следуещее фундаментальное уравнение для Т-матрицы радиационного многоканального рассеяния:
т = * +* Т. - а 1рк)(рк\т, (101)
г,к 0,к
где = 1/%/—2£Ук — эффективное главное квантовое число, соответствующее энергии электрона
Еук = Ее — УШу - кш/. (102)
При этом энергия падающего электрона считается малой и лежит в диапазоне 0 < Ее < ш, где и — частота колебаний иона.
Матрица реакций 1 описывает взаимодействие электрона с ионным остовом в присутствии электромагнитного поля и удовлетворяет интегральному уравнению
t = V + УС1;. (103)
Гриновский оператор в является гладкой функцией от энергии Е и включает ридберговские состояния е~ + ХУ+ за вычетом полюсов от дискретных уровней кулоновского центра, диссоциативные состояния X + У* и состояния валентных неридберговских конфигураций, взаимодействующих с ридберговскими и диссоциативными состояниями. Матрица t в уравнении (101) представляет собой сумму двух операторов, характеризующих электростатическое и полевое взаимодействия
t = te + t/. (104)
Второе слагаемое в (104), ответственное за взаимодействие с внешним полем с точностью до квадратичных по взаимодействию У^ членов имеет
вид
Ь* = [1 + Vе С (В)] [V1 + [1 + С(Е)Vе]. (105)
Полное сечение реакции диссоциативной рекомбинации для перехода из начального состояния г/о есть
а^Ее, ве) = £ °ЯоЛ(Ее, 0е), (106)
к,0
где начальное состояние системы до описывает падающий электрон, движущийся вдоль оси г в лабораторной системе координат с начальным импульсом ре = \/2Ее, аЯофк(Ее,ве) — парциальное сечение, которое определяется как
8?Г2 ^ (5) аяоРк = х
i! expmD^iO, в,ir- ч>Шв„ Ve)T^Pk{l\ -> b\p) \ ,
lm,\ /
St = arg[r(l + 1 - i/pe)}.
(107)
Здесь — спиновое вырождение для /3-конфигурации, I и А — угловой момент падающего электрона и его проекция в системе координат, связанной с молекулой, D^x(a, /3,7) — обобщенная сферическая функция, в и (р, ве и <ре — углы, задающие соответственно направление оси молекулы и пучка падающих электронов относительно направления вектора f, символ (...) означает усреднение по ориентациям молекул. При этом все промежуточные ридберговские состояния с учетом правил отбора автоматически включены в определение амплитуды перехода Т^ (IX —»■ IpXp).
Поскольку поле эффективно воздействует только на состояния промежуточного комплекса, то с учетом (95) здесь можно выделить следующие четыре механизма реакции:
(М1) — безызлучательный переход (к = 0);
(М2) — свободносвязанный переход с излучением (поглощением) кванта
поля (к = ±1) на стадии образования комплекса; (МЗ) — связанно-связанный переход с излучением (поглощением) света на
стадии его распада; (М4) — свободно-связанный переход, отвечающий поглощению фотона (к = — 1) с заселением дважды возбужденных промежуточных ридбергов-
Рис. 10. Зависимость сечения реакции е~ + О% —> 0(3Р) + 0[}И) от энергии падающего электрона Ее для различных значений напряженности / внешнего поля: топкая сплошная кривая соответствует / = 0, квадраты — / = 5-Ю5 В/см, кружки — / = 5-Ю6 В/см, толстая сплошная кривая — / = 5 • 107 В/см, треугольники — / = 108 В/см.
ских серий Пк(1~ки{Ук ^ 0), 3Е~(а4П„) на стадии образования комплекса.
В качестве конкретного примера приведем результаты расчетов динамики распада по диссоциативному терму в диссоциативный континуум Шумана-Рунге 3Р + 1П. Расчет сечения был выполнен в 7-канальном приближении (пять открытых каналов (и = 0), один закрытый (у = 1) и один открытый диссоциативный /3-канал). В интересующей нас области энергии Ее ^ 0.6 эВ резонансные вибронные состояния отсутствуют, и использование такого приближения вполне оправдано.
Зависимость сечения реакции распада в нанал 3Р + 1Б от энергии падающего электрона Ее для различных значений напряженности поля /, рассчитанная для величин ш/ = 24040 см-1 и 9е = тг/2, представлена на Рис.10. В отсутствие поля (/ = 0) наблюдается характерная резонансная фано-фешбаховская структура, образованная в результате безыз-лучательного перехода с заселением промежуточной ридберговской серии 7гоР7Ги(3Е~, г>о = 1). Скобками (по,г>о) и стрелками обозначены соответствующие резонансные линии. Первый резонанс с провалом в окрестности Ее = 0.3 эВ соответствует главному квантовому числу щ = 9. Частота и>у здесь выбирается таким образом, чтобы индуцированные полем резонан-сы П1ртти (3Е~, У\ = 0) от входных каналов движения электрона ва, с?сг, ¿-к и ¿8 были расположены вблизи этого провала. Как и следовало ожидать,
и>/,см 1
Рис. 11. Зависимость сечения реакции е~ +0% —► 0(3Р) + 0(10) от частоты и/ внешнего поля. Показаны парциальные вклады трех механизмов реакции. Толстая сплошная кривая соответствует полному сечению реакции. Расчеты выполнены для величины напряженности внешнего поля / = 5 -107 В/см, угла между направлениями вектора напряженности и электронного пучка ве = 7г/2 и энергии падающего электрона Ее = 0.03 эВ.
полевые резонанеы достаточно четко проявляются, начиная с / ^ 5 ■ Ю-6 В/см. Они соответствуют П\ =ь 3 и расщеплены за счет интенсивного полевого взаимодействия. Наибольший эффект полевого воздействия в этой области спектра достигается при / ~ 5 • 108 В/см.
На Рис.11 изображена зависимость сечения реакции —>• 0(3Р)+
0(ХВ) от частоты u>f внешнего поля. Согласно этим данным механизм (М1) является наиболее эффективным в диапазоне частот 13000 см-1 ^ а;/ ^ 22000 см-1, роль механизма (М2) является значительной в диапазоне частот 22000 см-1 ^ ^ 25000 см-1. Вклад механизма (МЗ) в диапазоне видимого света оказывается незначительным. Полевое воздействие оказывается наиболее эффективным для энергии электрона Ее = 0.03 эВ при частоте и??ах = 24040 см-1.
Зависимость сечения реакции е- + —> 0(3Р) + 0(10) от напряженности внешнего поля / приведена на Рис.12 для следующих значений параметров:
£е = 0.03эВ, = 24040 см-1, ве = тг/2, (108)
при которых эффект является максимальным и сечение реакции возрастает на пять порядков (см. Рис.10). Расчитанная кривая правильно передает
<70/3,10 13 см2
Рис. 12. Зависимость сечения реакции е 4- О$ —> 0(3Р) Н- 0(1!)) от напряженности / внешнего поля.
основные особенности полевого воздействия на реакцию диссоциативной рекомбинации: квадратичный рост сечения при малых значениях /, достижение максимума в точке /тах = 1.15-108 В/см и спад, пропорциональный I//2) ПРИ больших значениях /. Такая ситуация является типичной для радиационных и столкновительных процессов с участием промежуточных ридберговских состояний.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
1. Разработан метод расчета вибронных поверхностей потенциальной энергии и полной волновой функции квантовой системы "ридбергов-ский атом + молекула"в широкой области изменения координат.
2. Разработан метод расчета оптического потенциала для описания динамики взаимодействия ридберговского атома с нейтральной частицей при медленных столкновениях, включая процесс тушения.
3. Разработан метод расчета заселенностей ридберговских состояний квазимолекул А**Х2 (X = N,0) в неравновесной двухтемпературной рекомбинационной плазме в зависимости от потока, плотности и температуры свободных электронов.
4. Разработан метод расчета принимаемой на приемнике полной мощности потока некогерентного дециметрового излучения ридберговских квазимолекул в неравновесной плазме верхней атмосферы Земли для различных классов солнечной активности на основе созданной программы "Ридберг"и впервые дано объяснение экспериментально наблюдаемых закономерностей.
5. Впервые рассчитан спектр некогерентного сверхфонового излучения верхней атмосферы Земли в диапазоне 0.8-100 ГГц и обнаружены особенности типа "осиная талия "вблизи частоты 1.4 ГГц и "бутылочное горлышко "в окрестности частоты 5 ГГц в широкой области изменения параметров двухтемпературной плазмы, которые активно используются в дистанционной пассивной локации поверхности Земли и спутниковой связи.
6. Впервые рассмотрена низкотемпературная реакция диссоциативной рекомбинации медленных электронов и молекулярных ионов, протекающая в поле монохроматического лазерного излучения в видимом диапазоне. В качестве конкретного примера рассчитаны зависимости сечения реакции е~+02 от энергии падающего электрона, напряженности и частоты внешнего электромагнитного поля. Показано, что при определенном выборе этих параметров сечение реакции возрастает на несколько порядков, что позволяет говорить о возможности ее лазерного стимулирования.
Список работ, опубликованных по теме диссертации
1. Golubkov G.V., Golubkov M.G., Ivanov G.K. Difficulties inherent in the theory of low-temperature dissociative recombination and possible ways of their obviations / Dissociative Recombination: Theory, Experiment and Application. Ed. by M. Larsson, J.B.A. Mitchell, I.F. Schneider // Singapore, New Jersey, London, Hong Kong: World Scientific, 2000. — P. 106-110. - 298 P.
2. Golubkov G.V., Devdariani A.Z., Golubkov M.G. Collision of Rydberg atom A** with ground-state atom B. Optical potential // Письма в журнал экспериментальной и теоретической физики. — 2002. — Т. 75. — N. 7. - С. 385-387.
3. Голубков Г.В., Голубков М.Г., Романов А.Н. Диссоциативная рекомбинация электронов и молекулярных ионов 02 в поле интенсивного лазерного излучения в видимой области спектра // Журнал экспе-
Риментальной и теоретической физики. — 2002. — Т. 121. — N. 3. — . 573-582.
4. Голубков Г.В., Девдариани А.З., Голубков М.Г. Столкновение ридбер-говского атома А** с атомом В, находящимся в основном электронном состоянии. Оптический потенциал // Журнал экспериментальной и теоретической физики. - 2002. - Т. 122. — N. 6. - С. 1146-1157.
5. Golubkov G.V., Golubkov M.G., Romanov A.N. Dissociative Recombination of Slow Electrons and 02 Molecular Ions under the Action of Monochromatic Laser Radiation in the Frequency Range 13000-22000 cm-1 // Russian Journal of Physical Chemistry. - 2002. — V. 76. - Suppl. 1. - P. S115-S125.
6. Golubkov M.G., Golubkov G.V., Romanov A.N. Dissociative Recombination of Slow Electrons and Molecular Oxygen Ions in a Strong Laser Field
/ Dissociative Recombination of Molecular Ions with Electrons. Ed. by
S.L.Guberman // New York: Kluwer Academic/Plenum Publishers, 2003.
- P. 197-202. - 473 P.
7. Golubkov G.V., Golubkov M.G., Romanov A.N., Buenker R.J. Dissociative Recombination e~ -f 02 —>■ 0(lD) + 0(3P) in a strong laser field // Physical Chemistry and Chemical Physics. — 2003. - V. 5. - N. 1. - P. 3174-3182.
8. Golubkov G.V., Golubkov M.G., Buenker R.J. Dissociative recombination of slow electronswith molecular oxygen ions in the intensive laser field // Химическая физика. - 2004. - Т. 23. - N. 2. - С. 53-60.
9. Голубков Г.В., Иванов Г.К., Голубков М.Г. Возмущение состояний ридберговского атома А*" нейтральными частицами на промежуточных расстояниях // Химическая физика. — 2005. — Т. 24. — N. 6. — С. 3-13.
10. Голубков Г.В., Иванов Г.К., Голубков М.Г. Высоковозбужденные электронные состояния наномасштабных систем // Химическая физика.
- 2005. - Т. 24. - N. 9. - С. 3-18.
11. Golubkov M.G., Adamson S.O., Apukhtina N.V., Golubkov G.V., Dementiev A.I., Ryabinkin I.G. The Reaction of Dissociative Recombination in a Strong Light Field / Spectral Line Shapes. V. 15. N. 1058.19th International Conference on Spectral Line Shapes. Ed. by M.A.Gigosos, M.A.Gonzalez
// Mellville, New York: American Institute of Physics Conference Proceedings, 2008. - P. 137-139. - 390 P.
12. Golubkov G.V., Golubkov M.G., Ivanov G.K. Rydberg Atom A** in a Field of Neutral Atom В / Spectral Line Shapes. V. 15. N. 1058. 19th International Conference on Spectral Line Shapes. Ed. by M.A.Gigosos, M.A.Gonzalez // Mellville, New York: American Institute of Physics Conference Proceedings, 2008. - P. 140-142. - 390 P.
13. Адамсон С.О., Бюнкер Р.Дж., Голубков Г.В., Голубков М.Г., Дементьев А.И. Лазерное стимулирование низкотемпературной реакции диссоциативной рекомбинации электронов и молекулярных ионов кислорода // Химическая физика. - 2009. - Т. 28. - N. 4. — С. 26-42.
14. Golubkov G.V., Golubkov M.G., Ivanov G.K. Rydberg States of Atoms and Molecules in a Field of Neutral Particles / The Atmosphere and Ionosphere: Dynamics, Processes and Monitoring. Ed. by V.L.Bychkov, G.V.Golubkov, A.I.Nikitin // Dordrecht, Heidelberg, London, New York: Springer, 2010. - P. 1-67. - 375 P.
15. Голубков Г.В., Голубков М.Г., Чудновский Л.С., Карпов И.В. Обнаружение радиоизлучения, индуцированного импульсным источником гамма-квантов // Радиотехника. — 2010. — N. 2. — С. 81-85.
16. Голубков Г.В., Голубков М.Г. Потенциальные энергии квазимолекулы Na** + He // Химическая физика. - 2010. - Т. 29. - N. 3. - С. 14-26.
17. Голубков Г.В., Голубков М.Г., Бюнкер Р.Дж. Лазерное управление низкотемпературной реакцией е~ + О? // Журнал экспериментальной и теоретической физики. - 2011. - Т. 139. - N. 2. - Р. 221-227.
18. Бессараб Ф.С., Кореньков Ю.Н., Голубков М.Г. Роль возбужденных газовых компонент в атмосфере Земли. Горячий кислород // Хими-
ческая физика. - 2011. - Т. 30. - N. 5. - С. 16-23.
19. Голубков Г.В., Голубков М.Г., Карпов И.В. Микроволновое излучение атмосферы, индуцированное импульсным гамма-источником // Химическая физика. - 2011. - Т. 30. - N. 5. — С. 61-74.
20. Голубков Г.В., Голубков М.Г., Карпов И.В. Девдариани А.З. Излучение высоковозбужденных атомов и молекул в верхних слоях атмосферы // Химическая физика. - 2011. - Т. 30. - N. 5. - С. 75-79.
21. Golubkov G.V., Golubkov M.G., Buenker R.J. Laser control of the low-temperature e~ + Oq reaction // Химическая физика. - 2011. — Т. 30.
- N. И. - С. 46-49.
22. Малышев Н.С., Голубков Г.В., Голубков М.Г., Бюнкер Р.Дж., Либер-ман Х.П. Потенциальные энергии квазимолекулы К** + Не // Химическая физика. - 2011. - Т. 30. - N. 11. - С. 57-66.
23. Голубков Г.В., Голубков М.Г., Манжелий М.И. Микроволновое излучение верхней атмосферы Земли в периоды сильных геомагнитных возмущений // Химическая физика. — 2012. — Т. 31. — N. 2. — С. 31-47.
24. Голубков Г.В., Голубков М.Г., Манжелий М.И. Микроволновое и инфракрасное излучения верхней атмосферы в периоды повышения солнечной активности // Доклады академии наук. — 2012. — Т. 447. — N. 5. - С. 503-506.
25. Buenker R.J., Golubkov G.V., Golubkov M.G., Karpov I.V., Manzheliy M.I. Relativity Laws for the Variation of Rates of Clocks Moving in Free Space and GPS Positioning Errors Caused by Space-Weather Events / Global Navigation Satellite Systems — From Stellar to Satellite Navigation. Ed. by A.H.Mohamed // Berlin: InTech, 2013. - P. 3-48.
26. Golubkov G.V., Golubkov M.G., Devdariani A.Z. Quenching of Rydberg states in slow collisions with neutral atoms and molecules of medium // Вестник Санкт-Петербургского государственного университета. Серия 4,- 2013. - N. 4. - С. 149-152.
27. Голубков Г.В., Голубков М.Г., Манжелий М.И. Сверхфоновое излучение D-слоя атмосферы в диапазоне 0.8 - 6.0 ГГц // Доклады академии наук. - 2013. - Т. 452. - N. 5. - С. 510-513.
28. Голубков Г.В., Голубков М.Г. Возмущение высоковозбужденных состояний атома полем нейтральной частицы // Химическая физика.
- 2014. - Т. 33. - N. 2. - С. 42-51.
29. Голубков Г.В., Голубков М.Г., Манжелий М.И. Ридберговские состояния в D-слое атмосферы и ошибки позиционирования системы GPS // Химическая физика. - 2014. — Т. 33. - N. 2. — С. 64-77.
30. Golubkov G.V., Golubkov M.G., Manzhelii M.I., Karpov I.V. Optical Quantum Properties of GPS Signal Propagation Medium — D Layer / The Atmosphere and Ionosphere: Elementary Processes, Monitoring, and Ball Lightning. Ed. by V.L.Bychkov, G.V.Golubkov, A.I.Nikitin // Dordrecht, Heidelberg, London, New York: Springer, 2014. — P. 1-68. —
Заказ № 61-Р/12/2014 Подписано в печать 17.12.14 Тираж 120 экз. Усл. пл. 2,4
ООО "Цифровичок", г. Москва, Большой Чудов пер., д.5 тел. (495)649-83-30 www.cfr.ru ; e-mail: zakpark@cfr.ru