Метод обобщенного потенциала нулевого радиуса в квантовой механике и его приложения в задаче рассеяния электрона на молекуле тема диссертации и автореферата по ВАК РФ 01.04.02, кандидат физико-математических наук Ялунин, Сергей Валерьевич

• Ялунин, Сергей Валерьевич
• кандидат физико-математических наук
• 2004

Оглавление диссертации кандидат физико-математических наук Ялунин, Сергей Валерьевич

I Введение

II Метод ПНР в задаче рассеяния с возбуждением молекулярных состояний

1 Многоканальное рассеяние на электрона на молекуле с учетом парциальных Е волн.

2 Рассеяние на двухатомных, гомоядерных молекулах с учетом двух каналов.

3 Приближение адиабатических ядер.

4 Дифференциальное сечение рассеяния для электронных переходов.

5 Дифференциальное сечение рассеяния для электронно-колебательных переходов.

6 Интегральное сечение для электронных переходов. it 7 Интегральное сечение для электронно-колебательных переходов.

Введение диссертации (часть автореферата) на тему «Метод обобщенного потенциала нулевого радиуса в квантовой механике и его приложения в задаче рассеяния электрона на молекуле»

Все важнейшие процессы в природе являются следствием взаимодействия атомов, молекул, элементарных частиц между собой, причем в большинстве случаев во взаимодействии участвует большое количество частиц. Важной частью описания такого взаимодействия является квантовая теория рассеяния. В этой теории понятие сечения вводится как наблюдаемая величина, связанная с вероятностью процесса, происходящего при парном столкновении. Сечение некоторого процесса столкновения с частицей определяется как число событий, происходящих в единицу времени, когда на одну частицу мишени приходится единичный поток налетающих частиц (см. книгу [37], стр. 309). Процесс, для которого находится сечение, можно описать полностью, если известны квантовые состояния, а также скорости продуктов реакции. В этом случае сечение называют дифференциальным. В большинстве экспериментальных методов дифференциальное сечение измеряется для хорошо определенных величин, таких как энергии налетающего Ei и рассеянного Ef электронов, угол рассеяния. Если энергетическое разрешение налетающего пучка и природа молекулярной мишени таково, что квантовые начальное (г ) и конечное {/) состояния являются определенными, то дифференциальное сечение рассеяния на угол в дается формулой (в атомных единицах me = h = qe = 1) f-i{9)=kf= где ki и kf - начальные и конечные импульсы электрона, F{9) - амплитуда рассеяния. Единица измерения сечения равна квадрату Боровского радиуса. Реально достижимое энергетическое разрешение современных спектрометров работающих с перекрестными пучками электронов и молекулярных мишеней составляет около 10 мэВ, что позволяет различать только вращательные состояния молекулы водорода

Н2), у которой энергия вращательного возбуждения порядка 45 мэВ [39]. Кроме того, количественные эксперименты не вовлекают поляризованные по спину пучки или мишени, которые, как часто говорится, были бы "одеты"в один тонкий или сверхтонкий уровень, не говоря уже об одиночном вращательном уровне. Поэтому измеренные сечения включают в себя усреднение по возможным начальным и суммирование по конечным состояниям.

Интегральное сечение для процесса г —► / получается интегрированием дифференциального сечения по всем углам рассеяния

Общее интегральное сечение равно сумме всех интегральных сечений, в которой суммирование выполняется по всем энергетически доступным процессам.

Мощным средством исследования межъчастичного взаимодействия и измерения дифференциальных сечений является метод перекрестных пучков [74]. В таких экспериментах сечение связано с измеряемой интенсивностью достаточно сложным выражением затрагивающим геометрические и приборные характеристики. Сечение усредненное по угловому и энергетическому разрешению, можно выразить как произведение измеряемой интенсивности на поправочный множитель. При идеальных инструментах и геометрических параметрах множитель представляет собой константу, то есть не зависит от энергии и угла рассеяния. Однако, на практике конструктивные особенности установки приводят к зависимости множителя от угла рассеяния. Кроме того, электронная оптика и, следовательно эффективность определения сигнала, во многих случаях зависит от энергии, что приводит к энергетической зависимости множителя от энергии. Определение этого множителя является нетривиальной задачей, поэтому часто используют технику относительного течения, позволяющую избежать проблем с определением этого множителя. Эта техника использует отношение интен-сивностей для мишени и стандартного газа и известное дифференциальное сечение для стандартного газа. В качестве стандартного газа обычно используется инертные газы Не, Ne, Аг, Кг, а также молекулярный водород Ffe. Критическое требование в этом методе заключается в том, что область взаимодействия, обычно определяемая как область перекрывания электронного и молекулярного пучков, должна быть одинакова для обоих газов. В типичном эксперименте по столкновению используются молекулярные пучки диаметра 1-5 мм при управляющем давлении 0,02-10 торр. Для получения пучков используют специальные массивы стеклянных капилляров.

Фундаментальную важность в физике лазеров, плазмы и химии имеют процессы взаимодействия электрона с возбужденными состояниями молекул. Несмотря на это, изучение таких процессов находится на ранней стадии развития. Положение дел в этой области отражает трудности, связанные с получением возбужденных молекул в нужном количестве и их определением. Один из первых количественных экспериментов по рассеянию на молекулах в возбужденном состоянии [75] был проведен на метастабильном а*Д состоянии молекулы кислорода (Ог). Основной причиной для такого выбора стало большое время жизни этого состояния с учетом радиационной и столкновительной дезактивации. Для получения молекулярного пучка содержащего около 7% молекулярного кислорода в состоянии а1 Дд использовался микроволновый разряд. В работе [76] для изучения процесса диссоциативного прилипания к колебательно-возбужденным состояниям СОг использовалось два электронных пучка. Вначале молекулярный пучок перекрещивали с первым электронным пучком, ^ который возбуждал молекулы, после чего его перекрещивали со вторым электронным пучком, в котором и наблюдалось диссоциативное прилипание. Другой оригинальный подход был использован в работе [77]. В ней изучалось сечение ионизации метастабильного состояния А32 молекулярного азота (N2) электронным ударом в энергетическом интервале от порога до 240 эВ. Для получения метастабильного состояния использовался процесс нейтрализации. Пучок ионов N с энергией 1 кэВ нейтрализовался N2 и окисью азота (N0). Обнаружено, что пучок, сформированный в процессе нейтрализации N+NO, содержит смесь основного и возбужденного (А3£) состояний N2 в соотношении один к одному. Сечение ионизации основного состояния азота Х*Е N2 было найдено по известными экспериментальным данным. Измеряя ионизацию пучка смешанных состояний и исключая вклад от ионизации основного ^ состояния удалось получить сечение ионизации возбужденного состояния в довольно широком интервале энергий.

В связи с возникшим интересом к рассеянию на колебательно-возбужденных молекулах представляется разумным использовать тепловые процессы для заселения колебательных уровней. В этом случае, при фиксированной температуре окружающей среды, степень заселения возбужденных колебательных уровней будет определяться распределением Больцмана. Для простых двухатомных молекул энергия возбуждения колебательных уровней велика (>150 мэВ), поэтому при комнатной температуре процент молекул находящихся в возбужденном состоянии очень мал. Однако, уже у простых многоатомных молекул энергия колебательного возбуждения много меньше и, поэтому, нижние колебательные уровни будут заселены даже при 4 комнатной температуре. Первое количественное экспериментальное изучение электронного рассеяния на колебательно-возбужденных молекулах было осуществлено в работе Бакмана и др. [78]. Заселение колебательных уровней осуществлялось прямым нагреванием области взаимодействия. При температуре 310 К только 3% молекул СОг находилось в нижней колебательной моде с энергией возбуждения Е(0щ=83 мэВ, соответствующей изгибу молекулы, но уже при 573 К заселенность моды увеличивалась в три раза. Возбужденные колебательные и электронные состояния могут проявлять новые свойства, что приводит к существенному изменению сечения с увеличением температуры. Например, нижняя колебательная мода ССЬ обладает ненулевым усредненным электрическим дипольным моментом, в то время как основное состояние не обладает дипольным моментом. Поэтому, в результате нагревания в сечении рассеяния электронов (Е<2 эВ) начинает доминировать дипольный вклад, приводящий к резкому его увеличению. В работе Джонстоуна и др. [79] впервые сообщается об измерении дифференциального сечения электронного рассеяния. В работе использовался 280 Вт газовый СОг-лазер, который легко позволял получать температуру до 773 К. Похожие эксперименты на изоэлектронном партнере CQ - окиси азота (N2O) - показали менее драматическое, по сравнению с CCfe, изменение сечения с изменением температуры. Это связано с тем, что молекула Г^О нелинейная и имеет постоянный дипольный момент в основном состоянии. Поэтому индуцированный дипольный момент будет играть менее значительную роль, чем в ССЬ и, соответственно приведет к менее драматическому изменению сечения при увеличении температуры.

Измерение сечений возбуждения представляет отдельную задачу. Впервые, интегральное сечение возбуждения молекул электронным ударом было получено оптическим методом. В этом методе число возбужденных молекул, производимых электронным пучком в единицу времени, оценивается по интенсивности спектральных линий фотоэмиссионного спектра. Возбужденные молекулы, как правило, нестабильны по отношению к фотоэмиссии, поэтому через короткое время они, испустив квант света, переходят либо в основное, либо в другое, разрешенное правилом отбора, состояние. Например, в работе Зубка [86] изучалось сечение возбуждения С?Пи состояния N2. Переход в основное состояние запрещен, поэтому 7-распад осуществляется в несколько стадий: СЯП,,—>В3Пд и В3П9—»A3EU. В результате распада происходит заселение метастабильного состояния. Число зафиксированных в эксперименте фотонов (интенсивность линий) можно вычислить по формуле v'v" "v' в которой NV>V»(E) - интенсивность линии перехода из колебательного состояния г/ в состояние v" , NV>(E) = Ienlcrvi(E) - число молекул возбуждаемых в единицу времени пучком электронов с энергией Е и интенсивностью 1е,п- концентрация молекул, I -длина области возбуждения, Av>v>> - вероятность перехода из колебательного состояния v' в состояние v", tv> - время жизни электронно-колебательного состояния к', а г](\) представляет квантовую эффективность определяющей оптической системы. Определение множителя 1еп1г}{\) представляет весьма сложную задачу, поэтому, для нормировки сечений, в работе [86] использовалось максимальное значение сечения фотоэмиссии NV'V»(E) = Ienlrj(X)a^n(E) для линии v' = v" = 0 при энергии электронов Е = 14.25 эВ. Значительную проблему при таком подходе представляет собой выбор параметров и tv>, который может приводить к значительному увеличению погрешности. В этой работе использовались не сами параметры, а их произведение, которое было известно с хорошей точностью. Если эти параметры неизвестны то, ^ обычно, ограничиваются изучением возбуждения одного колебательного состояния.

Например, в работе Зубка и др. [87] изучалось сечение возбуждения нижнего колебательного состояния Ь3£+ угарного газа (СО). Результат был нормирован по сечению возбуждения при 14 эВ полученному в работе [88].

Нужно заметить, что неучтенный каскадный вклад от распада энергетически более высоких электронных состояний в изучаемое состояние может привести к переоценке сечения электронного возбуждения.

В отличие от интегрального, дифференциальное сечение нельзя определить с помощью оптического метода. В настоящее время, единственный известный путь определения дифференциальных сечений (упругих и неупругих) заключается в измерении интенсивности электрического тока при разных углах рассеяния. Как указыва-f1 ется в обзорах Брангера и др. [69, 70], наиболее логично и эффективно определять абсолютные дифференциальные сечения неупругого рассеяния на газовых мишенях в ходе измерения относительных (по отношению к упругим) неупругих интенсивно-стей с последующей их нормировкой на дифференциальное сечение упругого рассеяния. Самый сложный этап в этом подходе состоит в измерении интенсивностей электрического тока в прямом и обратном направлениях. В прямом направлении поток, практически невозможно корректно отличить рассеянные электроны от падающих, поэтому втт = 5°-10°. В обратном направлении этому мешают конструктивные особенности установки, поэтому, обычно, максимальный угол не превышает0гаах = 140°.

Чтобы определить суммарное сечение по дифференциальному нужно знать его при всех углах. В связи с этим, особое внимание уделялось развитию различных 4 способов экстраполяции дифференциального сечения на весь диапазон углов рассеяния. При таком подходе, точность определения суммарного сечения в существенной степени зависит от способа экстраполяции. Первый способ экстраполяции основывается на сшивании с численными результатами. При этом подразумевается, что эти результаты доступны. Второй способ, применяется для аппроксимации дифференциального сечения упругого рассеяния [80]. Этот путь привлекает формализм фазовых сдвигов, развитый в теории упругого рассеяния электронов на атомах, а также квантовый [81] или квазиклассический [82] методы обратной задачи рассеяния. Процедура экстраполяции предложенная в этих работах, предполагает, что потенциал взаимодействия электрона с молекулой является сферически симметричным, что заведомо неправильно. Однако, существуют случаи, когда такое приближение является разумным. Например, для высших парциальных волн, когда рассеивающий потенциал при больших г является почти сферически симметричным. Это вызвано тем, что высшие парциальные волны определяются, главным образом, взаимодействием на больших расстояниях. Вклад высших парциальных волн доминирует, когда дифференциальное сечение сосредоточено в области малых углов. В другой работе (Сан и др. [83]) приближение сферического потенциала не используется. Для описания угловой зависимости парциальных волн предлагается использовать первое Борнов-кое приближение для фиксированных в пространстве ядер молекулы. Вся динамическая информация о рассеянии содержится в вещественной /С-матрице, которая легко параметризуется через фазовые сдвиги и параметры поворота, а также обеспечивает унитарность ^-матрицы. Однако, в отличие от приближения сферического поля, число параметров (элементов /^-матрицы) сильно увеличивается, что приводит к переопределенной задаче. Чтобы избежать противоречий, в работе [83] делается предположение, что вклад всех матричных элементов, соответствующих высшим парциальным волнам определяется первым Борновским приближением. В частности, для e-N-i рассеяния было найдено, что все недиагональные элементы можно положить равными нулю. Исключение составляют элементы, для которых \1—1'\ = 2. В итоге, остается несколько свободных параметров, соответствующих низшим парциальным волнам.

Значительное продвижение было сделано в последние годы группой экспериментаторов из Манчестера. Эта группа разработала методику, позволяющую измерять сечение рассеяния в обратном направлении. В отличии от большинства спектроскопистов, которые стараются избегать магнитного поля как напасти, они приложили локальное магнитное поле к области столкновений электронного спектрометра, чтобы отклонить траектории налетающих и рассеянных электронов. В этом подходе магнитное поле индуцируется двумя коаксиальными соленоидами. Во внешней и ^ внутренней обмотках токи текут в противоположных направлениях, поэтому снаружи соленоида магнитное поле отсутствует. Внутри соленоидов магнитное поле имеет довольно сложную радиальную зависимость, в частности, при пересечении обмотки внутреннего соленоида магнитное поле меняет направление. В результате, траектории электронов подкручиваются вначале в одну сторону, а после прохождения внутреннего соленоида в другую сторону. Подбирая радиусы соленоидов и силу тока, удается вывести поток электронов рассеянных в прямом и обратном направлениях под нужными углами. Ряд точных экспериментов, проведенных на атомах инертных газов (Зубек и др. [89]), в которых было измерено сечение вплоть до0 = 180°, показали систематическое расхождение с аналогичными теоретическими результатами. Причем, если дифференциальное сечение имеет минимум при 180°, то результаты могут ^ различаться в два раза. В качестве примера, можно привести результаты измерений

90], стр. 93) для атомов неона при энергии электронов в 5 эВ. Дифференциальное сечение достигает максимума при 50°, после чего монотонно убывает, вплоть до 140°, что великолепно согласуется с известными теоретическими результатами. После этого экспериментальное сечение начинает увеличиваться, а теоретическое продолжает' уменьшаться. В итоге, при 180° экспериментальное значение в два раза превосходит теоретические результаты. Нужно заметь, что теоретические результаты, полученные разными группами находится в прекрасном согласии между собой при всех углах рассеяния.

Для измерения дифференциального сечения возбуждения используется спектральный анализ энергетических потерь электронов. Одной из главных проблем в примене-Ф нии этой техники является перекрытие нескольких колебательных линий с близкими колебательными энергиями. Характерным примером является сильное перекрытие особенностей спектров энергетических потерь для электронных Е?П и W3A состояний азота N2. Разность энергий возбуждения этих состояний меньше 10 мэВ. Для извлечения сечения из спектра энергетических потерь электрона был разработан метод спектральной деконволюции [84]. Использование этого метода базируется на приближении Франка-Кондона. В этом приближении сечение возбуждения на каждый колебательный подуровень v некоторого электронно-возбужденного состояния представляется в виде произведения сечения электронного возбуждения в приближении неподвижных ядер и некоторого множителяqv (фактор Франка-Кондона). Множитель qv представляет собой вероятность перехода из начального колебательного 4 состояния в соответствующее колебательное состояние v электронно-возбужденного состояния. Чтобы извлечь из спектра энергетических потерь интенсивность возбуждения Хп, соответствующую электронному состоянию п, спектр аппроксимируется . по методу наименьших квадратов функцией

N М(п')

N(E, в, W) = ]Г Nn,(E, в) J2 qv>F(Wn,v, -W) + В(Е, в, W), п'=0 v'—O в которой N - сигнал, полученный спектрометром при определенном угле рассеяния 1д и энергии электронного пучка Е, зависящий от энергии W, теряемой электроном при рассеянии. Величины Nn' пропорциональны дифференциальным сечениям возбуждения состояний п', N - общее число электронных переходов, М(п') - число колебательных подуровней, Wniv> - энергия каждого колебательного подуровня. Нормированная на единицу функция F описывает экспериментальный профиль линии и характеризует разрешающую способность спектрометра. Часто она выбирается в виде гауссовой функции с полушириной Д. Функция В представляет фоновый вклад в измеренном сигнале, который аппроксимируется полиномом по W. В процессе вычислений величины Wn>v' и qv> считаются известными, а значения Nn>, В и Д рассматриваются как вариационные параметры. После их определения можно вычислить дифференциальное сечение. Величина сигнала N сложным образом зависит от плотности молекулярного и тока электронного пучков, поэтому для определения дифференциального сечения возбуждения используют нормировку на сечение упругого рассеяния da (й\ Мп'(в)d<T (й\

При этом предполагается, что сечение упругого рассеяния известно. Нужно заметить, что в случае, когда приближение Франка-Кондона нарушается, использование техники спектральной деконволюции может приводить к значительным ошибкам. Кроме этого, описанная техника выделения сигнала, сама по себе, вносит ошибки деконволюции, которые для сильно перекрывающихся спектров могут оказаться значительными. Характерным примером является попытка выделения сигнала из спектра энергетических потерь для состояний В?П и \¥3Д, проведенная двумя независимыми группами экспериментаторов (см. обзор [70]). Несмотря на то, что результаты измерения суммарного дифференциального сечения Е?П+\¥3Д согласуются дифференциальные сечения отдельных состояний существенно различаются. Сильное перекрытие спектральных особенностей характерно и для вращательных переходов.

Прилипание электронов с низкой энергией к молекулам является важным процессом в газовых диэлектриках и других средах включая эксимерные лазеры и земную атмосферу. Необходимым условием для осуществления диссоциативного прилипания является наличие устойчивых отрицательных ионов для атомов образующих молекулу. Характерным примером, таких атомов, являются галогены. Благодаря большому сродству галогенов к электрону сечение диссоциативного прилипания для галоген-содержащих молекул при малых энергиях может становиться очень большим, как и эффективность производства анионов галогенов, а так же атомов галогенов и галоге-носодержащих радикалов, которые являются очень важными компонентами многих химических реакций. В частности, эти продукты используются при травлении с помощью разряда. Также, процесс прилипания электрона к молекуле с последующим ее развалом на отдельные атомы или фрагменты возможен, если существуют устойчивые отрицательные ионы фрагментов. Примером таких процессов является диссоциативное прилипание к органическим молекулам. В связи с большой промышленной важностью, в последнее время идет активное экспериментальное изучение прилипания к простейшим органическим нитратам (Масон и Илленберг и др. [71], стр. 103), таким как нитрометан (CH3NO2) и нитроэтан (C2H5NO2). Характерной особенностью таких процессов является наличие нескольких каналов распада. Например, в результате диссоциативного прилипания к нитрометану образуется ряд устойчивых отрицательных ионов: NO2, О-, ОН~, CN~ и CNO". В связи с большим числом источников загрязнения окружающей среды (промышленность, автомобильные выхлопы, сигаретный дым) большую важность имеет изучение ацетонитрила (Cf^CN). Захват электрона ацетонитрилом приводит к анионам CtbCN"", CN~, CHCN~, C2N~ и CHJ. В 2001 году было обнаружено присутствие муравьиной (НСООН) и уксусной (СН3СООН) кислот в межзвездном пространстве и коме кометы Хейли-Боппа. Предполагается, что эти кислоты могут играть ключевую роль в образовании биомолекул в межзвездном пространстве (Масон, Илленберг и др. [71], стр. 102). Уксусная кислота и аммиак необходимы для производства простейшей биологически значимой аминокислоты - глицина. Понимание механизма образования этих аминокислот является важным в связи с тем, что они являются строительными блоками протеинов и ДНК. Присутствие этих молекул в пространстве указывает области возможного возникновения жизни.

Диссоциативное прилипание электрона к молекуле протекает через образование квазистационарного составного состояния [91, 92], которое, посредством диссоциации, приводит развалу молекулы на несколько фрагментов. Сечение диссоциативного прилипания (Бардсли [93]) с развалом молекулы на два фрагмента можно выразить через асимптотическую форму волновой функции фрагментов Ф( Я), в которой расстояние между фрагментами обозначено через R. Расходящийся поток пар фрагментов в единицу телесного угла равен г>|Ф(Я)|2, где относительная скорость фрагментов v выражается через относительный импульс и приведенную массу фрагментов, выраженную в массах электрона. Поток налетающих электронов в единицу площади равен импульсу к, поэтому дифференциальное сечение прилипания находится по формуле

Суммарное сечение диссоциативного прилипания получается интегрированием по всему телесному углу.

Как и в экспериментах по рассеянию, очень часто для изучения диссоциативного прилипания используется техника скрещеных пучков. В последнее время (Клар [94]) была предложена новая техника, позволяющая получать монохроматический пучок электронов с низкой энергией и экстремально низким разбросом по энергии (до 50 мкэВ). Техника основывается на использовании резонансной лазерной фотоионизации метастабильных атомов аргона (Аг). Для получения фотоэлектронов использовались метастабильные атомы аргона в состоянии (4s) 3Рг и переход с возбуждением ионизационного континуума Аг+ 2Рз/2, осуществляемый через промежуточное состояние (4р) 30з. Вначале, пучок метастабильных атомов возбуждался диодным лазером, с одиночной модой (А=811,75 нм), на промежуточный уровень (4р) 30з. После чего, небольшая фракция возбужденных атомов в состоянии (4р)3Бз ионизировалась . перенастраиваемым многорежимным лазером. С помощью двулучепреломляющих фильтров разброс фотонов и, следовательно, фотоэлектронов по энергии уменьшается до 50-150 мкэВ. К сожалению, эта технология не позволяет варьировать энергию пучка в широком пределе. Например, в работе [94] максимальная энергия пучка составляла около 170 мэВ. В этой работе пучок электронов направлялся в область взаимодействия, где происходило прилипание к молекулам CCfeBr. Выход анионов выражается через интенсивность /е электронного пучка, концентрацию п молекул СС1зВг, длину области взаимодействия/, эффективность масс-спектрометра^ и сечение диссоциативного прилипания ое{Е) по формуле Ne(E) = Ienlr]<j(E). Множитель Ienlr] определялся через известное значение для коэффициента кинетики прилипания при тепловых энергиях ке(Т) = (vae). Скобки означают усреднение по энергиям тепловых электронов. Предполагается, что распределение тепловых электронов подчиняется закону Максвелла f(E) ос \/Ёе~Е/кТ. Исследования с экстремально монохроматическими пучками электронов позволили проверить некоторые теоретические предсказания. В качестве примера, можно привести ясное экспериментальное доказательство (Барсотти, Раф, Хотоп [90], стр. 43) р-волнового характера (сге ос у/Ё) порогового поведения электронного прилипания к молекуле хлора (СЬ). Ранее было предсказано, что при энергиях порядка 50 мэВ диссоциативное прилипание к СЬ и F2 протекает, главным образом, через анионное состояние с симметрией2Е„, в то время, как основное состояние обладает симметрией X1Eff. Из симметрии следует, что, при малых энергиях, электронный захват с образованием аниона возможен только для электронов в р-состоянии.

За прошедшие более чем полвека развития квантовой теории, были сделаны значительные успехи в развитии ab initio подходов в изучении электронной структуры молекул. В настоящее время вычислительные возможности позволяют проводить Л вычисления методом полного конфигурационного взаимодействия (FCI) с учетом до Ю10 определителей Слейтера [45] и тем самым судить о степени эффективности теории возмущения Меллера-Плессета четвертого порядка (МР-4) или, например, теории конфигурационного взаимодействия (CI) без сравнения с экспериментом. Несмотря на это, успехи в развитии методов для изучения электрон-молекулярных столкновений никак нельзя назвать большими. На это указывает, хотя бы то, что все они существенно различаются как по "философии", так и по численной реализации. Следует отметить, что, несмотря на все различия, все они имеют одну общую особенность. Все они используют современные технологии вычисления связанных состояний молекул в качестве вспомогательного инструмента. Многие идеи теории молекул, такие, как разложение iV-частичной волновой функции по детерминантам ^ Слейтера, перекочевали в теорию рассеяния. Например, во всех ab initio методах

N + 1)-частичная волновая функция представляется, согласно принципу Паули, в виде антисимметричной суммы произведений iV-частичных волновых функций молекулы Ф„(х), в которой координаты электронов обозначены через ж = (xi,. и состояний рассеяния переходящих на бесконечности в плоскую волну.

Обычно в качестве Ф„(х) выбираются состояния, возбуждение которых возможно с энергетической точки зрения. Для их представления используются различные приближения, такие, как метод Хартри-Фока (HF) или метод конфигурационного взаимодействия (CI).

Широкой популярностью пользуется вариационные методы, такие, как многоканальный метод Швингера. Многочисленные работы показали, его применимость для Ч вычисления упругого рассеяния на молекулах в газовой фазе, таких как 14, N2, метан

СН4) [46], силан (SiHj) [47], вода (Н20), а также неупругого рассеяния на небольших молекулах, таких как Н2 [48] и СО [49, 50]. В отличие от ранних методов, метод Швингера избегает прямого использования волновых функций сплошного спектра (фп) и включает граничные условия состояний рассеяния через интегральное уравнение Липпмана-Швиягера. При этом уравнение проектируется, с помощью проектора Р = ]Гт |Фт)(Фт|, на подпространство нескольких молекулярных состояний. Таким образом, метод Швингера учитывает короткодействующие обменное взаимодействие и эффект электрон-молекулярной корреляции в рамках ab initio подхода. Недавно, на примере 2П3 резонанса формы в упругом рассеянии электрона на молекуле N2 была показана способность метода Швингера учитывать эффекты дальнодействия (поляризации) для малых молекул. Приближенная волновая функция и, следовательно, сечения рассеяния находятся из условия стационарности некоторого функционала. В методе Швингера пробная волновая функция (фп) записывается в виде суммы гауссовых функций и плоских волн. Это позволяет вычислить все необходимые матричные элементы (за исключением интеграла от VGV, который содержит свободную функцию Грина G) в аналитическом виде. Для приближенного вычисления этого интеграла используется технический прием, который в литературе называется "а-вложением", в котором оператор VGV приближенно представляется в виде

В этом разложении через |о;) обозначены гауссовы функции. Это представление является приближенным, поскольку суммирование всегда ведется по конечному набору (базису) гауссовых функций. Эффективность а-вложения зависит от того, насколько удачно выбран базис гауссовых функций. Для вычисления матричных элементов от функции Грина G, она представляется в виде интеграла по плоским волнам. Этот прием называется "fc-вложением". В ряде работ отмечается, что для получения удовлетворительного результата достаточно вычислить главное значение интеграла[49]. Вариационный метод Швингера напрямую позволяет вычислить Т-матрицу перехода. Другие величины, такие как /('-матрица реакции и 5*-матрица могут быть получены из матрицы перехода. Решением задачи на собственные значения для /^-матрицы являются тангенсы фазовых сдвигов. Изучение зависимости фазовых сдвигов от

VGV « ^I/HHGVKck'IV. а,а' энергии позволяет определить резонансные структуры. Нужно заметить, что применение многоканального метода Швингера для изучения рассеяния и возбуждения на полиатомных молекулах весьма затруднительно. Поэтому ряд работ используют модельные потенциалы, которые в литературе носят название "псевдопотенциалов". Псевдопотенциалы использовались при изучения упруго рассеяния на полиатомных молекулах, таких как формальдегид СН20 [51], фторид азота NF3 [52], метан СН4, а также на галогенсодержащих соединениях СНДт, SiH„Xm в которых п + т = 4 и X=F, С1, Вг [53]. В этом методе псевдопотенциал представляет собой гладкую, часто конечную в нуле, функцию, которая совпадает с реальным потенциалом за пределами сферы некоторого радиуса R. Внутри сферы псевдопотенциал подбирается так, чтобы валентные псевдоволновые функции совпадали с действительными за пределами молекулы. Их интегралы от 0 до Я должны совпадать (сохранение нормы). Как правило, в качестве псевдопотенциалов используют функции, для которых вычисление матричных элементов в методе Швингера, выполняется аналитически.

Для многоканального рассеяния предложен алгебраический вариационный метод (Шнайдер, Ресчигно [54, 55]), основанный на применении принципа Конна к Г-матрице. В отличие от многоканального метода Швингера, многоканальный вариационный метод Конна имеет дело не с интегральным уравнением, а с дифференциальным уравнением Шредингера. Поэтому при реализации этого метода требуется выполнить расчет матричных элементов гамильтониана, тем самым отпадает необходимость вычислять интегралы от VGV. Пробная функция (фп) записывается через гауссов базис и "свободные"функции ( регулярные функции, которые на бесконечности переходят в решения уравнения Гельмгольца). Основной особенностью метода является технический прием позволяющий исключить вычисление обменных интегралов содержащих свободные функции. Это достигается в результате ортогонализа-ции базиса. При этом происходит "выедание дыр"в свободных функциях при малых г, при которых обменное взаимодействие наиболее ощутимо. Указанная особенность приводит к тому, что метод Конна для малых молекул, таких как Нг, оказывается немного быстрее метода Швингера. По-видимому, более очевидные преимущества метода Конна должны проявиться при вычислении рассеяния на более крупных молекулах. Кроме этого метод Конна допускает возможность учета эффектов корреляции и поляризации через введение оптического потенциала.

В 1983 году был предложен итеративный метод решения уравнения Липпмана-Швингера - метод цепных дробей для задачи рассеяния электрона на атоме. Вычисления рассеяния позитрона на атоме водорода (Н), электрона на атоме и ионе гелия

Не) с использованием модельных и реалистичных оптических потенциалов показали быструю сходимость метода даже для сильного взаимодействия. В других работах аналогичный подход применялся для изучения динамики ядерного рассеяния, упругого рассеяния электрона на молекуле Н2 [56] и диссоциативного прилипания к молекуле бромоводорода (НВг). Метод цепных дробей был применен и для неупругого рассеяния на молекуле Н2 (Ли и др. [57, 58]). В этом методе, в ходе каждой итерации производится коррекция волновой функции, матрицы реакции и потенциала. Подразумевается, что на каждом шаге происходит ослабление потенциала, а последовательность матриц реакции оказывается сходящейся. При изучении упругого рассеяния на Н2 трех итераций оказалось достаточно для получения результатов, которые очень хорошо согласовывались с экспериментом и другими теоретическими вычислениями. После трех итераций сходимость оказалась равной одному проценту. При низких энергиях (Е< 8 эВ) электронов существенный вклад в рассеяние оказывает поляризация молекул. Существует два различных способа учесть влияние поляризации молекулы при электронном рассеянии. В первом, поляризация объясняется виртуальными возбуждениями энергетически недоступных возбужденных состояний молекулы. Поэтому эффект поляризации трактуется как эффект влияния закрытых каналов в приближении сильной связи каналов. Другой путь - использовать метод поляризованных орбиталей - был реализован в рамках метода цепных дробей в работе [56]. В этом подходе поляризационное взаимодействие возникает из-за искажения электронных оболочек и, следовательно, наведения дипольного момента в присутствии возмущающего электрического поля электрона. На примере упругого рассеяния на Н2 было показано, что для получения хороших результатов достаточно учесть только дипольную часть кулоновского взаимодействия приг»<гт/+1, где через Ti и rjv+i обозначены расстояния до молекулярного и рассеиваемого электронов. При rN+i<Ti взаимодействие можно положить равным нулю. Обычно время релаксации электронного облака мало, поэтому искажение оболочек можно описывать в адиабатическом приближении (в котором электрон покоится и сосредоточен в некоторой точке rjv+i), причем форма самих оболочек определяется положением электрона. Коррекция искаженных оболочек с учетом движения электрона возможна в рамках теории возмущений. В качестве возмущения выступает оператор кинетической энергии электрона. В работе [58] рассматривается возбуждение состояния t?E„ молекулы Н2 в области энергий 10.5-40 эВ. Этот дипольно-запрещенный переход особенно чувствителен к применяемому методу, поскольку лишь несколько симметрий вносят вклад в сечение, причем осуществляется это исключительно через обменное взаимодействие. Отмечается хорошее согласие результатов метода цепных дробей как с экспериментом так и с аналогичными вычислениями в рамках многоканального метода Швингера. Из-за недостаточной точности эксперимента, небольшое расхождение с результатами метода Швингера не позволяет делать выбор в пользу какого-либо метода.

Построение состояний рассеяния (iftn) при прямом интегрировании уравнения Шредингера реализовано в обобщенном "одноцентровом" подходе [59]. В этом методе молекулярные орбитали, входящие в волновые функции молекулы Фп(х), а также состояний рассеяния (фп) разлагаются по сферическим гармоникам относительно центра масс молекулы

1,т

После подстановки разложений в уравнение Шредингера и интегрирования по угловым переменным получается система сцепленных интегро-дифференциальных уравнений для радиальных функций (ynim). Одноцентровый подход плохо передает поведение волновой функции вблизи ядер, что приводит к его неэффективности при оценке общей энергии и других наблюдаемых, чувствительных к поведению волновой функции вблизи ядер. Однако, для задач рассеяния ядерные сингулярности менее важны, поэтому одноцентроный подход может приводить к хорошим результатам. Построение радиального матричного потенциала представляет самый напряженный этап метода, особенно трудоемким является вычисление обменной его части. Интегрирование системы осложняется интегральной частью, поэтому решение осуществляется с помощью итераций [60, 61]. После каждого шага вычисляется /('-матрица и оценивается сходимость процедуры. Обычно для удовлетворительной сходимости требуется 20-30 итераций. Одноцентровый подход с удовлетворительными результатами был применен к упругому рассеянию электрона на метане (СИ) [60], силане (SiH4) [59], а также на CF4 и SF6.

Одноцентровый подход слегка напоминает метод Д-матрицы, предложенный в 1947 году Вигнером. В этом методе все пространство делится сферой на две области. Внутри Д-матричной сферы рассматриваются кулоновское и обменное взаимодействия между электронами. Во внешней области взаимодействие записывается в виде мультипольного разложения. Предполагается, что в этой области волновая функция молекулы обращается в ноль. Во внутренней области решается задача на собственные значения и функции с дополнительным граничным условием на сфере (обычно полагают радиальную производную равной нулю). Для решения задачи используются стандартные алгоритмы диагонализации гамильтониана в некотором базисе. Собственные значения и функции нужны для построения Я-матрицы на сфере как функции энергии, связывающей радиальную производную многоканальной функции У = (Уыт) с самой многоканальной функцией

Ddy y = RTr•

Во внешней области система дифференциальных уравнений решается с учетом R-матричного граничного условия. Во многих случаях достаточно учесть только один мультипольный член - поляризационный. Метод Д-матрицы применялся для исследования упругого и неупругого рассеяния на молекулах Цз [62, 63], N2, 02 [63] и СО [64]. В работе [62] рассматривалось рассеяние на молекуле водорода с одновременным учетом шести нижних возбужденных состояний (от t?E„ до С1!!,,). Несмотря на значительное расхождение с экспериментом присутствует качественное согласие. Особенно плохое согласие с экспериментом было отмечено для дифференциального сечения возбуждения состояния Ь3Е*. Идея Д-матрицы оказалась плодотворной для модельного описания различных процессов. Например, в 1991 году Фабрикант предложил полуэмпирическую версию резонансной теории Д-матрицы для описания диссоциативного прилипания электрона к молекуле. Этот подход позволяет изучать температурную зависимость сечения диссоциативного прилипания, а также судить о влиянии поляризации и дипольного взаимодействия на прилипание электрона к полярным молекулам. В рамках этой теории было рассмотрено диссоциативное прилипание к метил галогенам СЕзВг, СН31 [65, 66], а также к CF3CI [67].

Другой полуэмпирический подход, вовлекающий понятие граничного условия -метод потенциала нулевого радиуса - был фактически использован Ферми в 1934 году [1, 2] при описании движения нейтронов в водородосодержащих средах. Годом позже подход был использован при изучении фоторасщепления дейтрона в работе Бете и Пайерлса [3]. В этой работе потенциал взаимодействия нейтрона и протона впервые был заменен граничным условием на волновую функцию. Другая важная работа, в которой явно выписывалось граничное условие, выполнена Брейтом[4]. В ней описывается математически корректная постановка задачи и решение. Применение метода при исследовании дифракции нейтронов на кристаллах описывается в работе Гольдбергером и Зейтцем [5]. В этой же работе впервые получено выражение для энергии зоны в импульсном представлении. В ранних работах для вычисления амплитуды рассеяния использовалось приближение Борна и поправки к нему. Для нуклонов приближение Борна, обычно, дает хорошие результаты, поскольку длина волны нуклонов очень мала по сравнению с характерными расстояниями между ними. Важным достоинством метода является то, что он позволяет находить точные решения и, тем самым, учитывать перерассеяние, что делает его особенно актуальным в теории электронных и атомных столкновений. Исследование многократного рассеяния на двух потенциалах провели в 1952 году Чу и Вик в работе [6]. Большое количество интересных приложений метода в физике атомных процессов было выполнено в Санкт-Петербургском университете и изложено в книге Демкова и Островского [7]. Приложение метода нужно начинать с рассмотрения асимптотики волновой функции. При рассеянии электрона на атоме при низких энергиях(Е = \к2) можно ограничиться рассмотрением ^-состояния. Вдали от атома асимптотика волновой функции хорошо описывается формулойгф = sin(fcr+5), где фазовый сдвиг обозначен как rj. Если речь идет об анионе, то асимптотика волновой функции гф = В ехр(—/Зг). В ряде случаев, для решения задачи достаточно знать не саму волновую функцию, а только ее асимптотику. Конечно, это не относится к такой наблюдаемой, как энергия, поскольку асимптотика относительно слабо влияет на величину функционала энергии. Тем не менее, существуют физические величины, такие, как квадруполь-ный момент, поляризуемость, которые определяются, в основном, областью больших расстояний электрона от атома. В качестве характерного примера, можно привести работу Армстронга [8], в которой получено сечение фотоотрыва для иона Н~ (атомном аналоге дейтрона). Сечение фотоотрыва пропорционально квадрату матричного элемента дипольного момента для перехода из связанного (анионного) состояния в состояние континуума. В области максимума метод дает результат, который не уступает по точности ab initio расчету. В то же время, указанные асимптотики волновых функций, по существу, являются решениями Гельмгольца, которые удовлетворяют граничному условию в нуле г-,/lV fcctg(<5). гф)' гф г=0

В случае связанного состояния fcctg(<S) нужно заменить на -/? и добавить еще одно граничное условие гф = В. Основная проблема метода - подбор параметров таких, как kctg(5), а, В. При малых энергиях можно воспользоваться аппроксимацией fcctg(<5) ~ —а, где через а обозначена обратная длина рассеяния. Параметр /? можно выразить через энергию сродства электрона к атому, а В взять из ab initio расчета. Аналогично, чтобы описать рассеяние на нескольких атомах нужно каждый атом заменить своим граничным условием. Хотя рассеяние на нескольких атомах не является сферически симметричным, метод позволяет легко находить парциальные волны, фазовые сдвиги, матрицу рассеяния. Причем, в работе [9] показано, что если сущеш ствует симметрия, то парциальные волны и фазы классифицируются по представлениям группы симметрий, а вырождение определяется размерностью представления. Показывается, что при малых энергиях рассеяние на любом числе атомов аналогично сферически симметричному случаю, поэтому парциальные волны и фазы можно классифицировать по азимутальному числу I. В качестве примера рассматривается рассеяние на кубе. В работе [10] парциальные волны используются для анализа дифракции электрона на тетраэдре. При рассмотрении рассеяния на нескольких атомах нужно иметь в виду, что атомы могут перекрываться, тогда параметры в граничных условиях меняются.

В классической формулировке метода потенциалы действуют только на волны с нулевым орбитальным моментом. Это обстоятельство делает возможным примене-^ ние модели только при низких энергиях, при которых основной вклад дают5-волны.

Более того, существуют случаи, когда уже при очень малых энергиях основной вклад в сечение осуществляется Р-волнами. Например, в рассеянии электрона на атомах щелочных металлов, таких как К, Rb, Cs при очень низких энергиях (до 100 мэВ) существует Р-резонанс формы. Первая попытка обобщения метода на высшие орбитальные моменты была предпринята в 1957 году в работе Хуанга и Янга[11, 12], в которой, однако, не выписывается граничное условие в явном виде, а рассматривается псевдопотенциал пропорциональный J-функции и тангенсу парциальной фазы Si. В рамках метода и теории возмущений, в этой же работе изучается Ферми и Бозе статистики для многочастичных систем в конечном объеме. Выражения для амплитуды рассеяния для псевдопотенциалов приводится в работе [13]. В этих работах не выписывается граничное условие в явной форме, хотя его можно записать так г\ (2i+1) O'/l г=0

Математические особенности спектральной теории для псевдопотенциалов рассматривается в работе [14]. Обобщение метода на случай Р-волны упоминается в работе [16]. Интересное приложения метода возникает в связи с изучением адиабатических кривых для Ридберговских состояний молекул Rt&, Cs2. В 2000 году было предсказано, что состояния Ридберга сформированные в результате столкновения Ридбер-говского атома с атомом в основном состоянии могут быть стабильными на очень больших расстояниях, сравнимых с размером Ридберговского атома. Хотя энергия связи этих состояний очень мала, они характеризуются огромным дипольным моментом, что позволяет их производить и управлять ими во внешнем поле в конденсате Бозе-Эйнштейна. Изучение влияния внешних электрических и магнитных полей на стабильность слабосвязанных состояний рассмотрено в работе Друкарева[17]. В частности, показано как электрическое и магнитное поля позволяют управлять связанными состояниями.

Другое ограничение метода связано с возбуждением атомов и молекул, а также их ионизацией. Первоначальная формулировка метода ограничивается рассмотрением упругого рассеяния. Необходимость в многоканальной постановке задачи возникает уже при попытке учесть спин взаимодействующих частиц. Например, хорошо известно, что рассеяние электрона на атомах первого ряда таблицы Менделеева (Н, Li, К, Rb, Cs) существенно зависит от общего спина системы, который равен либо нулю (синглетное рассеяние), либо единице (триплетное рассеяние). Во многих случаях учет вращательных и колебательных степеней свободы молекулы можно осуществить в рамках приближения Борна-Опенгеймера. Вычисления по методу потенциалов нулевого радиуса с учетом этого приближения проведены Друкаревым и Юровой в работе [18]. Получены сечения возбуждения вращательно-колебательных переходов для двухатомных молекул. В рамках метода установлено, что при небольших энергиях сечение вращательного перехода из j в j' пропорционально квадрату коэффициенту Клебша-Гордана C(j 21 j', ОО)2, в котором, обычно, полагают 1 = 1 (см. [19], стр. 64). При энергиях меньше первого порога возбуждения на взаимодействие электрона с молекулой могут влиять высшие виртуальные молекулярные состояния, приводящие к появлению предпороговых резонансов, которые нельзя обнаружить в рамках приближения Борна-Опенгеймера, Такое влияние можно рассматривать как эффект закрытых каналов в многоканальной постановке задачи. В 1983 году Островский и Устимов [20] показали, что в рамках метода учет вращения молекулы можно выполнить точно. Для этого они рассмотрели рассеяние на жестком ротаторе и получили аналитическое выражение для амплитуды вращательных переходов в виде суммы бесконечного ряда. Изучая полученный результат они показали, что при очень низких энергиях эффекты неадиабатичности могут приводить к захвату электрона в долгоживущие вращательные состояния молекулярного иона. Эти состояния приводят к появлению в сечении рассеяния очень узких резонансов (в диапазоне 20200 мэВ). В работе отмечается, что в действительности, благодаря взаимодействию с колебательными степенями свободы, время жизни резонансов может быть значительно меньшим. В 1970 году Демков и Островский [21], с целью описать возбуждение высших электронных состояний, предложили использовать матричное гранич-ffH ное условие на векторную волновую функцию, компоненты которой соответствует некоторому возбужденному состоянию. При этом они ограничились рассмотрением ы й'-волн. В качестве примера рассмотрено возбуждение t?Eu состояния молекулы Н2 в приближении неподвижных ядер. К настоящему моменту, выполнено несколько достаточно надежных экспериментов и теоретических исследований этого дипольно-запрещенного перехода. В частности, установлено, что в области максимума сечение возбуждения примерно равно о max ~ 2 а. е., в то же время расчет по методу потенциалов нулевого радиуса дает сильно заниженное значение (« 0.7 а.е.) Это расхождение указывает на необходимость дальнейшего развития метода.

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

Д) спектра для конечных систем до зонального спектра для бесконечных систем). Метод позволяет изучать влияние размера и формы атомной системы на положение энергетических уровней, зон, расположение энергии Ферми (из расчета, что каждый уровень, согласно принципу Паули, заполняется двумя электронами), ширины запрещенных зон [22]. Последние работы в этом направлении были выполнены де Прунелем [23, 24]. В них используется сепарабельные потенциалы, которые, в некотором смысле, можно рассматривать как частный случай метода потенциалов нулевого радиуса (см. Главу 2). Интересные результаты [25] получены для непериодических атомных систем в форме цепочек Фибоначчи, Коха, мозаики Пенроуза. Для цепочек Фибоначчи и Коха (последняя представляет самую простую фрактальную фигуру) при увеличении числа атомов происходит образование трех энергетических зон. Для цепочки Коха положение энергии Ферми как раз ниже главной зоны. В отношении волновой функции основного состояния установлено, что она локализуется в окрестности нескольких атомов.

Ч.

Заключение диссертации по теме «Теоретическая физика», Ялунин, Сергей Валерьевич

Заключение

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

Выводятся основные формулы для дифференциального и интегрального сечений возбуждения молекулы из колебательного состояния vq в электронно-колебательные состояния n,v с учетом усреднения по начальным вращательным состояниям и суммирования по конечным вращательным состояниям. Выведена формула для сечения возбуждения чистых электронных состояний п - сумма сечений возбуждения состояний п, v для всевозможных v. При выводе использовалось условие полноты колебательных гармоник. В рамках метода, выводятся формулы для дифференциального и интегрального сечений возбуждения чистых электронных состояний с простой аналитической зависимостью от угла рассеяния и энергии столкновения. Учет колебательного движения ядер, в рамках приближения Борна-Опенгеймера, сводится к простому усреднению сечений по начальному колебательному состоянию 0, v. Далее выводится выражение для дифференциального сечения для электронно-колебательных переходов. Ответ представляется в виде бесконечной суммы. При выводе формулы для интегрального сечения использовался подход позволяющий избежать вычисления 3jm символов или коэффициентов Клебша-Гордана. Конечный результат для интегрального сечения электронно-колебательных переходов представлен в виде двойной бесконечной суммы.

В третьем параграфе рассматривается приложение метода к реакции электронно-колебательного возбуждения триплетного состояния п =а3Е+ молекулы водорода.

Возбуждение этого триплетного состояния носит нерезонансный характер, поэтому вычисления проводились в приближении двух каналов. Часть параметров (параметры адиабатических кривых основного и возбужденного состояний) взяты из справочной литературы или эксперимента. При подборе параметров из граничного условия использовалось усредненное по спину значение длины рассеяния, которое ранее уже использовалось в работе Друкарева и Юровой. Приведено обоснование его использования. Два параметра подбирались так, чтобы получить максимальное согласие с дифференциальными сечениями в диапазоне 13-17 эВ, полученными в методе MCF. В этом диапазоне интегральное сечение достигает максимума. При 15 эВ наблюдается хорошее количественное согласие с другими расчетами. Используемые параметры позволяют найти интегральные и дифференциальные сечения в диапазоне энергий 12-24 эВ. Результаты метода можно сравнить с другими расчетами, а также с экспериментом при 20 эВ.

Во второй главе формулируется и решается задача многоканального рассеяния на потенциале конечного радиуса действия. Рассеяние описывается матричной волновой функцией, которая удовлетворяет многоканальному уравнению Липпмана-Швингера. Чтобы описать рассеянную волну в области, в которой потенциал равен нулю, вводится вспомогательная матричная функция

1 iKr °° ^ / ' \ п / 7 I М

1=0 п=0 4 ' v ' представляющая бесконечную сумму по полиномам Лежандра Pi (х). Функция Н играет важную роль и при нахождении матрицы амплитуды рассеяния на нескольких потенциалах. Выведена система уравнений, обобщающая уравнения метода потенциалов нулевого радиуса на случай системы неперекрывающихся рассеивателей. Система позволяет обобщить метод на случай матричных потенциалов, позволяющих рассматривать переходы с изменением симметрии. Выводится обобщенное матричное граничное условие, позволяющее рассматривать переходы с изменением орбитального момента.

Решение системы позволяет найти матричные амплитуды рассеяния на фиксированных ядрах. Результат выражается через коэффициенты Клебша-Гордана. Конкретный вид амплитуды возбуждения молекулы в лабораторной системе зависит от четности молекулярного состояния и орбитального момента. Выводятся формулы для дифференциальных и интегральных сечений для электронно-колебательных переходов. При статистическом усреднении учитывается вырождение состояний с ненулевым значением проекции орбитального момента |т|. Полученные формулы для электронно-колебательных возбуждений, полученные с учетом усреднения по колебательному движению, позволяют изучать не-Франк-Кондоновское поведение возбуждения связанное с многократным рассеянием. Рассматривается возбуждение с3Пи состояния молекулы водорода медленными электронами.

В третьей главе рассматривается возможность изменения свойств потенциала с помощью процедуры "одевания". Решение задачи рассеяния для потенциалов нулевого радиуса приводит к элементам ^-матрицы при к Е С л>- ч в»-»***1 ехр{2щ) = ац + ik2l+1'

Одевание решений приводит к модификации граничного условия и, следовательно, обобщенного потенциала нулевого радиуса. Также, это выражается в появлении конечного радиуса эффективного взаимодействия. В третьем параграфе подымается вопрос о связи рассеяния на отдельном атоме с рассеянием на молекуле типа YXn. Получены выражения для обратных длин рассеяния. Построены графики иллюстрирующие изложенную процедуру. Рассматривается приложение метода к молекуле SiH4. Вклад высших парциальных волн при увеличения энергии столкновения оценивается при помощи разложения ctg^ в ряд Лорана при к = 0.

Список литературы диссертационного исследования кандидат физико-математических наук Ялунин, Сергей Валерьевич, 2004 год

1. Fermi Е., Sopra lo spostamento per pressione delle righe elevate delle serie spettrali (1. Italian), Nuovo Cim, 11 (1934) 157-166

2. Fermi E., Sul moto dei neutroni nelle sostanze idrogenate (In Italian), Ricerca Scientifica, 7 (1936) 13-52

3. Bethe H., Peierls R., Quantum theory of the diplon, Proc. Roy. Soc. (London), 148A (1935) 146-156

4. Breit G., The scattering of slow neutrons by bound protons I. Methods of calculation, Phys. Rev., 71 (1947) 215-231

5. Goldberger M. L., Seitz F., Theory of the refraction and the diffraction of neutrons by crystals, Phys. Rev., 71 1947 294-310

6. Chew G. F., Wick G. C., The impulse approximation, Phys. Rev., 85 1952 636-642

7. Демков Ю. H., Островский В. H., "Метод потенциалов нулевого радиуса в атомной физике", 1975 (Ленинград: ЛГУ) 240 е.; Demkov Yu. N. and Ostrovsky V. N., "Zero-Range Potentials and their Applications in Atomic Physics", 1988 (New York: Plenum)

8. Armstrong В. H., Empirical analysis of the H- photodetachment cross section, Phys. Rev. 131 (1963) 1132-1137

9. Демков Ю. H., Рудаков В. С., Метод парциальных волн для несферического рассеивателя, ЖЭТФ 59 (1970) 2035-2047

10. Szmytkowski R., Szmytkowski С., Scalar wave diffraction from zero-range scatterers, J. Math. Chem., 26 (1999) 243-254

11. Huang K., Yang C. N., Quantum-mechanical many-body problem with hard-sphere interaction, Phys. Rev. 105 (1957) 767

12. Хуанг К., "Статистическая механика", 1966 (Москва: Мир) 520 с.

13. Baltenkov A. S., The application of general zero-range potentials to multi-center problems, J. Phys. Lett. A 286 (2000) 92-99

14. Андреева Т. Г., Рудаков В. С., Рассеяние потенциалами нулевого радиуса парциальных волн с ненулевыми моментами, Вестник ЛГУ N22 (1977) 12-18

15. Khuskivadze A. A., Chibisov М. I., Fabrikant 1.1., Adiabatic energy levels and elecric dipole moments of Rydberg states of Rt>2 and Cs2 dimers, Phys. Rev. A 66 (2002) 042709, смотри ссылку Omont A., J. Phys. (Paris) 38 (1977) 1343

16. Drukarev G. F., The zero-range potential model and its application in atomic and molecular physics, Adv. Quantum Chem. 11 (1978) 251

17. Drukarev G. F. and Yurova I. Yu., Multiple-scattering approach to the vibrational-rotational excitation of molecules by slow electrons, J. Phys. B: At. Mol. Phys., 10 (1977) 3551-8

18. Lane N. F., The theory of electron-molecule collisions, Rev. Mod. Phys. 52 (1980) 29-119

19. Ostrovsky V. N., Ustimov V. I., Exact solution of the particle-rigid-rotor scattering problem, J. Phys B: At. Mol. Phys., 16 (1983) 991-1008

20. Albeverio S., Gesztesy F., H0egh-Krohn R and Holden H, Solvable Models in Quantum Mechanics, 1988 (New York: Springer-Verlag)

21. Bouju X., Prunele E., Size and shape effects on electronic energy levels: from infinite to nanoscopic systems in three-dimensional space, Phys. stat. sol. (b) 217 (2000) 819-832

22. De Prunele E., Bouju X., Fibonacci, Koch, and Penrose Structures: Spectrum of finite subsystems in three-dimensional space, Phys. stat. sol. (b)225 (2001) 95-114

23. Leble S. В., Yalunin S., Zero-range potentials in multi-channel diatomic molecule scattering, Phys. Lett. A 306 35-44; preprint quant-ph/0205110

24. Blatt J. M., Wiesskopf V. F., "Theoretical Nuclear Physics", 1952 (New York: Wiley), Section 2.3.С

25. Островский В. H., Модель сепарабельных потенциалов и термы Щ, Вестник ЛГУ, 10 (1971) 10-19

26. Базь А. И., Зельдович Я. Б., Переломов А. М., "Рассеяния, реакции и распады в нерелятивистской квантовой механике", 1971 (Москва: Наука) 544 с.

27. Drozdov S. I., Scattering of fast neutrons by nonspherical nuclei, Sov. Phys. JETP 1 (1955) 591-2

28. Chase D. M., Adiabatic approximation for scattering processes, Phys. Rev. 104 (1956) 838-42

29. Oksyuk Yu. D., An adiabatic approximation treatment of exitation of rotational levels in diatomic molecules during electron impact, Sov. Phys. JETP 49 (1965) 1261-73

30. Oksyuk Yu. D., Sov. Phys. JETP 22 (1966) 873-85

31. Chang E. S., Temkin A., Rotational excitation of diatomic molecules by electron impact, Phys. Rev. Lett. 23 (1969) 399

32. Варшалович Д. А., Москалев A. H. и Херсонский В. К., "Квантовая теория углового момента", 1975 (Ленинград: Наука)

33. Ландау Л. Д., Лифшиц Е. М., "Квантовая механика. Нерелятивистская теория", 1989 (Москва: Наука) 752 с.

34. Бруйар Ф., Мак-Гоуэн Дж., "Физика ион-ионных и электрон-ионных столкновений", 1986 (Москва: Мир) 432

35. Юрова И. Ю., Иванов В. Е., "Сечения рассеяния электронов на атмосферных газах", 1989 (Ленинград: Наука) 144 с.

36. Golden D. Е., Lane N. F., Temkin A., Gerjuoy Е., Low Energy Electron-Molecule Scattering Experiments and the Theory of Rotational Excitation, Rev. Mod. Phys. 43(4) (1971) 642-679

37. Zubek M., Mielewska В., King G., J. Phys. B: At. Mol. Opt. Phys. 33 (2000) L527-L532.

38. Szmytkowski R. and Szmytkowski C., Scalar wave diffraction from zero-range scatterers, J. Math. Chem. 26 (1999) 243-254

39. Jain A. K., Tripathi A. N. and Jain A., Low-energy electron-SiHt collisions using a parameter-free spherical optical potential, J. Phys. В 20 (1987) 389

40. Wan H. X., Moore J. H. and Tossel J. A., Electron scattering cross sections and negative ion states of silane and halide derivatives of silane, J. Chem. Phys. 91 (1989) 1340

41. Chung S., Lin С. C., Phys. Rev. A 17 (1978) 1874

42. Rossi E., Bendazzoli G. L., Evangelisti S., Maynau D., A full-configuration benchmark for the N2 molecule, Chem. Phys. Lett. 310 (1999) 530-536

43. Brescansin L. M., Lima M. A. P., McKoy V., Cross sections for rotational excitation of CH4 by 3-20 eV electrons, Phys. Rev. A 40 (1989) 5577-5582

44. Winstead C., McKoy V., Low-energy electron scattering by silane (SiHi), Phys. Rev. A 42 (1990) 5357

45. Lima M. A. P., Gibson T. L., McKoy V., Huo W. M., Cross sections for excitation of the b3£+, a3E+, and с3П„ states of H2 by low-energy electrons, Phys. Rev. A 38 (1988) 4527

46. Weatherford C. A., Huo W. M., Application of the Schwinger multichannel formulation to electron-impact excitation of the b?L+ state of CO, Phys. Rev. A 41 (1990) 186-196

47. Sun Q., Winstead С., McKoy V., Electronic excitation of carbon monoxide by low-energy electron impact, Phys. Rev. A 46 (1992) 6987-6994

48. Sobrinho A. A., Machado L. E., Michelin S. E., Mu-Tao L., Brescansin L. M., Elastic and excitation cross sections for electron-formaldehyde collisions, J. Mol. Str. 539 (2001) 65-74

49. Joucoski E., Bettega M. H. F., Elastic scattering of low-energy electrons by NF3, J. Phys. B: At. Mol. Opt. Phys. 35 (2002) 783-793;

50. Natalense A. P. P., Bettega M. H. F., Varella M. T. do N. Oliveira A. J. S., Azevedo D. L., Ferreira L. G., Lima M. A. P., Norm-conserving pseudopotential in electronic excitation of molecules by electron impact

51. Schneider В. I., Rescigno T. N., Complex Kohn variational method: Application to low-energy electron-molecule collisions, Phys. Rev. A 37 (1988) 3749-3754

52. Schneider В. I., Rescigno T. N., Electron-impact excitation of the t?E* state of H2 using the complex Kohn method: R dependence of the cross section, J. Phys. В 21 (1988) L691-L695

53. Lee M.-T., Fujimoto M. M., Iga I., Application if the method of continued fractions to low-energy electron scattering by the hydrogen molecule, J. Mol. Str. 394 (1997) 117-125

54. Lee M. Т., Fujimoto M. M., Iga I., Electronic excitation of thea3E* and с3!!,, states of Hi by electron impact using the method of continued fractions, J. Mol. Str. 432 (1998) 197-209

55. Lee M.-T., Fujimoto M. M., Kroin Т., Iga I., Electronic excitation of the tfE* state of H2 by electron impact using the method of continued fractions, J. Phys. B: At. Mol. Opt. Phys. 29 (1996) L425-431

56. Branchett S. E., Tennyson J., Morgan L. A., Differential cross sections of electronic excitation of molecular hydrogen using the /2-matrix method, J. Phys. B: At. Mol. Opt. Phys. 24 (1991) 3479-3490

57. Collins L. A., Schneider В. I., "Recent advances in the theory of electron impact excitation of molecules", Electronic and atomic collisions, H. B. Gilbody, W. R. Mewell, F. H. Read, А. С. H. Smith (eds.) (1988) 57-72

58. Morgan L. A., Sarpal В. K., Tennyson J., "Я-matrix calculations of inelastic electron scattering by diatomic molecules", Electron Collisions with Molecules, Clasters, and Surfaces, ed. H. Ehrhardt and L. A. Morgan 1994 (New York: Plenum Press) 15-29

59. Wilde R. S., Gallup G. A., Fabrikant 1.1., Comparative studies of dissociative electron attachment to methyl halides, J. Phys. B: At. Mol. Opt. Phys. 33 (2000) 5479-5492

60. Schramm A., Fabrikant I. I., Weber J. M., Leber E., Ruf M.-W., Hotop H., Vibrational resonance and threshold effects in inelastic electron collisions with methyl iodide molecules, J. Phys. B: At. Mol. Opt. Phys. 32 (1999) 2153-2171

61. Wilde R. S., Gallup G. A., Fabrikant I. I., Semiempirical i?-matrix theory of low energy electron-CF3Cl inelastic scattering, J. Phys. B: At. Mol. Opt. Phys. 32 (1999) 663-673

62. Khakoo M., Trajmar S., Electron-impact excitation of the с3П„, and СЧГи states of H2, Phys. Rev. A 34 (1986) 146

63. Brunger M. J., Teubner P. J. O., Differential cross sections for electron-impact excitation of the electronic states of N2, Phys. Rev. A 41 (1990) 1413-1426

64. Mason N. J., Newell W. R., The total excitation cross section of the <?Пи state of H2, J. Phys. В 17 (1986) L587-L591; Electron resonance structure observed in the с3Пи and BXS+ cross sections of H2, J. Phys. В 17 (1986) L203-L207

65. Lee Mu-Tao, Lucchese R. R., McKoy V., Electron-impact excitation and dissociation processes in H2, Phys. Rev. A 26, (1982) 3240-3248

66. Hall R.I., Trajmar S., Scattering of 4.5 eV electrons by ground X?Ej state and metastable a}&g oxygen molecules, J. Phys В 8 (1975) L293-L296

67. Srivastava S. K., Orient O. J., Double e-beam technique for collision studies from excited states: Application to vibrationally excited GQj, Phys. Rev. 27 (1983) 12091212

68. Armentrout P. В., Tarr S.M., Dori A., Freund R. S., Electron impact ionization cross section of metastable N2(A), J. Chem. Phys. 75 (1981) 2786-2794

69. Buckman S. J., Elford M. Т., Newman D. S., Electron scattering from vibrationally excited CO, J . Phys. В 20 (1987) 5175-5182

70. Johnstone W. M., Mason N. J., Newell W. R., Electron scattering from vibrationally excited carbon dioxide, J. Phys. В 26 (1993) L147-L152

71. Lun A., Chen X. J., Allen L. J., Amos K., Inversion of electron-water elastic-scattering data, Phys. Rev. A 49 (1994) 3788-3798

72. Lipperheide R, Fiedeldey H. Inverse problem for potential scattering at fixed energy. Journal Paper] Zeitschrift fur Physik a (Atoms and Nuclei), vol.286, no.l, 1978, pp.45-56

73. Fiedeldey H., Lipperheide R., Naidoo K., Sofianos S.A., Semiclassical and quantal inversion of nuclear scattering at fixed energy, Phys. Rev DC 30 (1984) P434

74. Sun W., Morrison M. A., Isaacs W. A., Trail W. K., Alle D. Т., Gulley R. J., Brennan M. J., Buckman S. J., Detailed theoretical and experimental analysis of low-energy electron-N2 scattering, Phys. Rev. A 52 (1995) 1229-1256

75. Trajmar S., Cartwright D. C., Williams W., Differential and Integral Cross Sections for the Electron-Impact Excitation of the al Ag and bl States of O2, Phys. Rev. A 4 (1971) 1482-1492

76. Nishimura H. (1988) (как указано в 48])

77. Zubek M., Excitation of the СЯПц state of N2 by electron impact in the near-threshold region, J. Phys. B: At. Mol. Opt. Phys. 27 (1994) 573-581

78. Zubek M., Olszewski R and Woliriski P, Electron impact excitation of the tf state of carbon monoxide in the energy region from thresold to 17 eV, J. Phys. B: At. Mol. Opt. Phys. 30 (1997) L791-L795

79. Zobel J., Mayer U., Jung K., Ehrhardt H., Pritchard H., Winstead C. and McKoy V., J. Phys. B: At. Mol. Opt. Phys. 29 839 (как указано в работе 87])

80. Christophorou L. G., Electron-Molecule Interaction and Their Applications, Vol. 1 (1984) Academic, Orlando FL

81. Illenberger E. and Momigny J., Gaseous Molecular Ions. An introduction to Elementary Processes Induced by Ionization (1992) Springer-Verlag, New York

82. Bardsley J. N., Configuration interaction in the continuum states of molecules, J. Phys. В (Proc. Phys. Soc.) ser.2 vol.1 (1968) 349-364

83. Fraga S., Ransil B. J., Studies in molecular structure, J. Chem. Phys. 35 (1961) 1967

84. Демков Ю. H., "Вариационные принципы в теории столкновений", 1958 (Москва: Физматгиз) 168 с.

85. Yamaguchi Y., Two-nucleon problem when the potential is nonlocal but separable. I, Phys. Rev., 95 (1954) 1628-1634

86. Тейлор Дж., "Квантовая теория нерелятивистских столкновений", 1975 (Москва: Мир) 566 с.

87. Martin P. A., Multiple scattering and the Rehr-Albers-Fritzsche formula for the propagator matrix, J. Phys. A: Math. Gen. 31 (1998) 8923-8932.

88. Zeller R., "Band-structure methods Unoccupied Electronic States"(Topics in Applied Physics 69) ed. J С Fuggle and J E Inglesfield 1992 (Berlin; Springer) p 24-49

89. Gonis A., "Green Functions for Ordered and Disordered Systems"1992 (Amsterdam: North-Holland)

90. Миллер У., "Симметрия и разделение переменных" 1981 (Москва: Мир) 342 стр.

91. Blatt J. М., Biedenharn L. С., The angular distribution of scattering and reaction cross sections, Rev. Mod. Phys. 24 (1952) 258-272

92. Giantureo F. A., Thompson D. G., The Ramsauer-Townsend effect in methane, J. Phys. В 2 N12 (1976) L383-L385