Траекторное моделирование динамики столкновительно-индуцированной диссоциации и ионной рекомбинации тема автореферата и диссертации по физике, 01.04.17 ВАК РФ

РОССИЙСКАЯ АКАДЕМИЯ НАУК ИНСТИТУТ ЭНЕРГЕТИЧЕСКИХ ПРОБЛЕМ ХИМИЧЕСКОЙ ФИЗИКИ

ТРАЕКТОРНОЕ МОДЕЛИРОВАНИЕ I, СТОЖНОВИТЕЛЬНО-ИНДУЦИРОВАННОЙ ДИССОЦИАЦИИ И ИОННОЙ РЕКОМБИНАЦИИ

01.04.17- химическая физика, в том числе физика горения и взрыва

АВТОРЕФЕРАТ диссертации на соискание ученой степени доктора физико-математических наук

На правах рукописи

АЗРИЕЛЬ Владимир Михайло!

Москва 2008 г.

003454062

Работа выполнена в лаборатории динамики элементарных процессов Института энергетических проблем химической физики Российской Академии наук

Научный консультант:

доктор физико-математических наук Л.Ю. Русин

Официальные оппоненты:

доктор физико-математических наук, профессор A.M. Борисов

доктор физико-математических наук, профессор А.И. Осипов

доктор физико-математических наук, профессор С.Я Уманский

Ведущая организация:

Институт проблем механики РАН

Защита диссертации состоится

"■/v^" ^-е^^Х 2008 г. в ^

es,

заседании диссертационного совета Д 002.012.02 при Институте химической физики им. H.H. Семенова Российской Академии наук по адресу: 119991 ГСП-1, Москва, ул. Косыгина, 4.

С диссертацией можно ознакомиться в библиотеке Института химической физики им. H.H. Семенова Российской Академии наук.

Автореферат разослан

1

2008 г.

Ученый секретарь

диссертационного совета Д 002.012.02, доктор физико-математических наук

С.М. Фролов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы

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

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

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

В помощь исследователю приходит компьютерное моделирование, адекватность которого реальному процессу проверяется сравнением результатов расчета с теми характеристиками процесса, которые непосредственно могут быть измерены в эксперименте. Конечной целью такого моделирования является восстановление поверхности потенциальной энергии (ППЭ), необходимой для решения квантово-механической задачи рассеяния в приближении Борна-Оппенгеймера о разделении квантовой системы на быструю и медленную подсистемы. Однако во многих случаях физически оправданным и точным является классический или квазиклассический подходы к математическому моделированию динамического поведения системы без решения чрезвычайно трудоемких квантово-мехалических задач. Так как длина волны де Бройля (*.= И/ту) мала по сравнению с расстояниями, на которые перемещаются атомы в самом акте реакции, то классические и квазиклассические методы работают достаточно точно, особенно для больших масс и энергий. И хотя поведение атомов и молекул при химических превращениях подчиняются строгим законам квантовой теории, подавляющее большинство таких процессов качественно правильно описываются уравнениями классической механики с квантовыми поправками, касающимися, как правило, внутренних состояний реагентов и продуктов взаимодействия.

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

Г

описывающей элементарный акт химического взаимодействия. Траекторное моделирование динамики элементарного процесса занимает промежуточное место между неэмпирическими расчетами и экспериментом и опирается как на теоретические представления о природе взаимодействия частиц, так и на эмпирические значения характеристик взаимодействия. Этот вид моделирования в настоящее время представляет собой наиболее распространенную, обладающую широкими пределами применимости процедуру выяснения связи непосредственно измеряемых характеристик рассеяния с динамическими особенностями каналов взаимодействия. Соответствующая связь устанавливается тем, что классические гамильтоновы уравнения движения ядер в заданном канале решаются на поверхностях потенциальной энергии взаимодействия в электронных состояниях, через которые система проходит между исходными (реагенты) и конечными (продукты) состояниями. Хорошее совпадение результатов динамических расчетов и экспериментальных данных свидетельствует о правильности используемой ППЭ и может служить гарантией достоверности расчетных характеристик процесса, не поддающихся прямому экспериментальному измерению.

Как отмечено выше, одной из важнейших задач, определяющих корректность результатов траекторного моделирования элементарных процессов, является построение ППЭ, адекватной физической картине исследуемого явления. Однако, определение поверхности потенциальной энергии, адекватной рассматриваемому процессу, представляет одну из наиболее трудных задач квантовой химии, если речь идет об "ab initio" расчетах ППЭ. В настоящее время, несмотря на огромные достижения в развитии вычислительной техники, известны лишь несколько таких ППЭ, построенных для наиболее простых систем, состоящих не более чем из трех атомов. Кроме того, полученные неэмпирическими методами потенциалы обычно представляют собой таблицы, что создает определенные неудобства при их использовании в вычислительных процедурах. В большинстве случаев эти табличные значения аппроксимируются гладкими математическими функциями, удобными в применении, но не позволяющими впоследствии проанализировать, какие именно характеристики поверхности потенциальной энергии ответственны за те или иные динамические особенности поведения исследуемой системы.

Поэтому большое распространение получили более доступные полуэмпирические и эмпирические методы построения ППЭ. Эти методы имеют два существенных преимущества в применении к динамическим исследованиям. Первое, они используют доступные экспериментальные данные для калибровки ППЭ, что позволяет воспроизвести топографические особенности поверхности, оказывающие наиболее существенное влияние на динамику процесса. И второе, они, как правило, допускают аналитическое представление ППЭ в виде функций, значительно более удобных с вычислительной точки зрения, что позволяет реализовывать более эффективные и быстрые алгоритмы и программы динамических расчетов. Поверхности потенциальной энергии, построенные с помощью таких методов, включают ряд подбираемых параметров, изменение которых позволяет существенным образом менять топологию ППЭ и, таким образом, исследовать влияние тех или иных топологических особенностей ППЭ на динамику элементарного взаимодействия.

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

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

Главные задачи, поставленные автором, могут быть сформулированы следующим образом:

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

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

3. Определить происхождение основных каналов, возникающих при столкновительной диссоциации в четырехатомной системе, и их взаимосвязь.

4. Разработать методологию и провести, основываясь на принципе микроскопической обратимости, исследование динамики прямой трехгельной рекомбинации ионов на поверхностях потенциальной энергии, количественно описывающих обратные процессы СИД в соответствующих системах.

5. Построить траекторные модели рассматриваемых процессов СИД и рекомбинации, позволяющие понять механизмы химических превращений и их зависимость от динамических и кинематических особенностей взаимодействия.

Предмет исследования

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

Достоверность

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

Научная новизна работы состоит в следующем:

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

2. Исследовано влияние параметров, определяющих отталкивание в потенциалах взаимодействия одноименных ионов, на топологию ППЭ и на динамику многоканального взаимодействия, условия реализации и сечения отдельных элементарных процессов.

3. Построены траекторные модели процессов СИД в четырехатомной системе, учитывающие как динамические (структура и топология ППЭ), так и кинематические (энергетические параметры, ориентационные и массовые эффекты) характеристики взаимодействия.

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

5. Создан высокопроизводительный и удобный в использовании пакет вычислительных средств для компьютерного моделирования процессов столкновительно-индуцированной диссоциации и рекомбинации.

Автор выносит на защиту:

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

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

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

Практическая значимость работы:

1. Результаты и выводы, полученные в работе, позволяют сделать заключение о высокой эффективности траекторных методов для исследования динамики элементарных процессов. Разработаны пакеты вычислительных средств для траекторного моделирования процессов СИД и трехтельной рекомбинации, обеспечивающие высокую производительность и точность расчетов при реализации на персональных компьютерах.

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

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

4. Использованная методология построения ППЭ для трех- и четырехатомных систем с успехом может быть применена для конструирования поверхностей потенциальной энергии для систем с большим числом атомов.

5. Построенные модели СИД и рекомбинации могут быть использованы для понимания основных закономерностей процессов, протекающих в верхних слоях атмосферы, низкотемпературной плазме и других практически важных областях естествознания.

Апробация результатов работы

Результаты, вошедшие в диссертацию, были представлены на следующих конференциях и симпозиумах: XII Международный симпозиум по молекулярным пучкам (Италия, 1989), III Международный симпозиум по элементарным процессам и химическим реакциям (Чехословакия, 1989), VIII Симпозиум по атомной физике и физике поверхности (Италия, 1992), IX Европейская конференция по динамике молекулярных столкновений (Чехословакия, 1992), Конференция НАТО по ориентационным и поляризационным эффектам в динамике химических реакций (Италия, 1992), XV Международный симпозиум по молекулярным пучкам (Германия, 1993), XVI Международный симпозиум по молекулярным пучкам (Израиль, 1995), XI Европейская конференция по динамике молекулярных столкновений (Дания, 1996), Международный симпозиум по применению компьютеров в химических исследованиях (Москва, 1996), XVII Международный симпозиум по молекулярным пучкам (Франция, 1997), Общероссийская конференция по теоретической химии (Казань, 1997), XII Европейская конференция по динамике молекулярных столкновений (Великобритания, 1998), X Конференция по современным тенденциям в вычислительной химии (США, 2001), VI Всемирный конгресс теоретически ориентированных химиков (Швейцария, 2002), XX Международный симпозиум по молекулярным пучкам (Португалия, 2003), XVIII Международная конференция по передаче молекулярной энергии (Испания, 2003), XIV Симпозиум по атомной физике, физике кластеров и поверхности (Италия, 2004), VIII Европейская конференция по атомной и молекулярной физике (Франция, 2004), XXTV Международный симпозиум по динамике разреженного газа (Италия, 2004), Международная конференция по динамике молекулярных систем (Нидерланды, 2004), VII Конгресс всемирной ассоциации теоретически ориентированных химиков (ЮАР, 2005), XXI Международный симпозиум по молекулярным пучкам (Греция, 2005), VI Международная конференция по химической кинетике (США, 2005), XV Симпозиум по атомной физике, физике кластеров и поверхности (Австрия, 2006), XIX Международный симпозиум по газовой кинетике (Франция, 2006), Европейская конференция по динамике молекулярных систем (Италия, 2006), IX Европейская конференция по атомам, молекулам и фотонам (Греция, 2007), XX Международная конференция по передаче молекулярной энергии (Франция, 2007), XXII Международный симпозиум по молекулярным пучкам (Германия, 2007), XXV Международная конференция по фотонным, электронным и атомным столкновениям (Германия, 2007), Симпозиум по диссоциативной рекомбинации (Нидерланды, 2007), XVI конференция по современным тенденциям в вычислительной химии (США, 2007), XVI Симпозиум по атомной физике, физике кластеров и поверхности (Швейцария, 2008).

Личное участие автора

Автор принимал личное определяющее участие в постановке задач и разработке подходов для их решения, создании пакетов вычислительных программ и проведении расчетов, а также анализе полученных результатов.

Публикации

Основные результаты диссертации опубликованы в 70 работах, список которых приведен в конце автореферата.

Объем и структура диссертации

Диссертация состоит из введения, пяти глав, заключения и списка литературы. Общий объем - 299 страниц, в том числе 211 рисунков и 5 таблиц. Список литературы содержит 417 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

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

В главе 1 диссертации приведено сравнение результатов траекторного моделирования СИД в системе СэС1 + ЯЫ с данными эксперимента в скрещенных молекулярных пучках и проанализировано влияние некоторых особенностей топологии ППЭ на основные характеристики рассеяния. Следует отметить, что рассматриваемая система является модельной и отражает основные особенности взаимодействия в широком круге систем, включающих молекулы с ионной связью.

Как показывает траекгорный счет, взаимодействие двух молекул с ионной связью, таких как СэС1 и ЛЫ, может протекать по 12 каналам, включающим образование как ионных, так и нейтральных продуктов [23,24,27,32,67]:

1. СэС1 + Ш -> СэС1 + Ш

2. ~>С5+ + СГ + НЬ+ + Г

3. -> Сэ+ + СГ + Ш

4. -» ЯЬ+ + Г + СвС1

5. -> Сэ:++сг

6. ->КЬС1 + Сз+ + Г

7. -»С^ + ЯЮ

8. -> СвСШЬ*+ Г

9. ->СвСи" + КЬ+

10. -»СзЖЬ^ + СГ

11. -»ШСГ + Сз+

12. СэСШЫ

Очевидно, что преобладающим по сечению каналом в этой схеме будет канал 1, включающий упругие и неупругие взаимодействия, поскольку каналы, в которых происходят химические превращения, реализуются лишь в определенных и достаточно узких диапазонах прицельных параметров, взаимных ориентаций реагентов и т.д. Два других канала, приводящих к образованию нейтральных продуктов (каналы 7 и 12), представляют большой интерес с точки зрения их динамики. Канал 12 соответствует образованию долгоживущего промежуточного комплекса, который затем распадается на исходные молекулы или продукты СвД и ЯЬС1. Такие же продукты образуются в канале 7 при непосредственном взаимодействии исходных молекул, причем без участия долгоживущего комплекса. По этим каналам отсутствуют экспериментальные данные, и динамика их реализации может быть получена только из адекватных траекторных моделей. В каналах 2-6 и 8-11 происходит образование ионных продуктов. Экспериментально полные сечения и угловые распределения были измерены для положительно заряженных ионных продуктов в области энергий столкновения (Ео™) до 9 эВ. Динамическая информация о рассеянии этих продуктов послужила основой для траекторного моделирования взаимодействия в системе СэС1 + ЯЫ Следует отметить, что канал 2, соответствующий диссоциации обеих молекул реагентов, имеет высокий энергетический порог и при энергиях столкновения ниже 9 эВ закрыт. Таким образом, в области энергий до 9 эВ регистрируемые в эксперименте со скрещенными молекулярными пучками положительно заряженные ионы Се* образуются в каналах 3,6 и 11, ионы Ш>+ - в каналах 4, 5 и 9, а положительно заряженные ионные комплексы СБСтЬ* и ЯЫСз+ образуются соответственно в каналах 8 и 10.

Используемая в расчетах поверхность потенциальной энергии построена как аддитивная функция шести попарных взаимодействий, каждое из которых описывается в рамках усеченного потенциала Риттнера. Эта модель является основной и наиболее часто используемой формой представления потенциалов взаимодействия ионов щелочных металлов и галогенов. С использованием атомной системы единиц (е=Й=т=1) усеченный потенциал Риттнера (Шйпег) имеет вид:

и(Ю = Е(Я) ± 1/Я-(а|+а2)/2114 - С/Я6 (1)

где Е(Я) - энергия отталкивания, принимаемая обычно в форме модели Борна-Майера (Вогп-Мауег)

Е(Я) = А-ехр(-Я/р) (2)

Здесь А - калибровочный множитель и р - параметр жесткости отталкивательной стенки, а| и аг - поляризуемости взаимодействующих ионов, С - дисперсионная постоянная и Я - расстояния между ионами. Знак "+" перед кулоновским членом в выражении (1) соответствует взаимодействию одноименно заряженных ионов, а знак "-" -взаимодействию ионов разного знака. Значения параметров парных потенциалов взаимодействия, взятые из литературных источников [см. ссылки в работах 23, 24, 32], приведены в таблице 1, а на рисунке 1 показана произвольная конфигурация столкновения и обозначены расстояния между взаимодействующими ионами.

Таблица 1. Параметры потенциалов парных взаимодействий для системы СбС1 + ЯЫ

Ионная пара 1 А, Р, С,

Сэ+ - СГ 1 136.266 0.6881 138.0

ЯЬ+-Г 2 172.9586 0.7127 200.0

СГ - Г 3 41.5247 0.995 527.18

Св+-ЯЬ+ 4 877.69 0.462 115.678

Сэ+ - У 5 179.6108 0.7448 305.0

ЯЬ+-СГ 6 135.3389 0.6541 92.0

Поляризуемости ионов: (Хс5=16.48, аа=24.93, ащ,=9.595, а]=52.8.

Использование значений параметров, указанных в таблице 1, приводит к качественно правильному описанию функций возбуждения всех измеряемых в эксперименте положительных ионов, как это можно видеть из рисунка 2(а-г). Хорошее качественное согласие рассчитанных и экспериментальных данных показывает, что использованные значения параметров парных потенциалов, равно как и концепция построения полной поверхности потенциальной энергии из парных потенциалов взаимодействия адекватно описывают рассеяние в четырехатомной системе, включая каналы перегруппировок ионов.

Из этого рисунка следует, однако, что сходство динамических расчетов с экспериментальными данными действительно только качественное. Рис. 2(а) показывает экспериментальную функцию возбуждения образования ионов Се* вместе с рассчитанной для суммы каналов 3+6+11, в которых образуются ионы Сэ+. Как можно видеть, обе функции заметно различаются по поведению во всем диапазоне энергий взаимодействия. Совпадение существенно улучшается, если исключить вклад канала 11 (см. рис. 3), однако практически нет существенных резонов считать, что этот канал не возбуждается в реакции, в то время как симметричный ему канал 9 образования ионов ЯЬ+ нельзя исключить из сравнения с экспериментом. Очевидно, что такое разногласие является следствием недостаточности динамической модели из-за неточности потенциала взаимодействия. Улучшение топологии поверхности потенциальной энергии может быть достигнуто уточнением параметров парных потенциалов, составляющих ППЭ.

0123456789 10 Энергия столкновения, эВ

0123456789 10 Энергия столкновения, эВ

0123456789 10 Энергия столкновения, эВ

СвЖЬ . •

Эксперимент

01 23456789 10 Энергия столкновения, эВ

Рис. 2. Сравнение экспериментальной и полученной из траекторных расчетов функций возбуждения: а) ионов Се*, б) ионов ЯЬ+, в) ионов СгСШЬ*, г) ионов СэЖЬ*.

Рис. 3. Сравнение экспериментальной и полученной из траекторных расчетов функции возбуждения ионов Сэ+ при учете вкладов каналов 3, 6 и 11 (сплошная линия) и только каналов 3 и 6 (пунктир).

Параметры четырех из шести парных потенциалов, соответствующих устойчивым молекулам, известны с высокой степенью точности из спектроскопических данных молекул и не могут варьироваться в значительных пределах. К наименее достоверным величинам, входящим в описание поверхности потенциальной энергии, относятся параметры отгалкивательной стенки, задаваемой уравнением (2), в потенциалах взаимодействия одноименных ионов и, в частности, значения параметров жесткости рз и Р4 в парных потенциалах взаимодействия СГ-Г и Сз+-ЯЬ+ соответственно.

Расчеты показывают, что увеличение рэ приводит к резкому уменьшению сечения канала 11. При этом сечение канала 3 заметно увеличивается, а сечение канала б меняется в небольших пределах. Функции возбуждения каналов 3, 6 и 11 при различных значениях параметра рз приведены на рис. 4.

О ' ................................о1Л,—.....Г .—................о1 ч I-.' .Тг"?-^ ,1,1—_

2 4 б В 10 12 14 16 18 20 22 24 2 4 6 1 10 12 14 16 18 20 2! 24 0 2 4 6 8 10 12 14 16 18 20 22 24

Энергия столкновения, эВ Энергия столкновения, эВ Энергия столкновения, эВ

Рис. 4. Функции возбуждения каналов 3 (а), б (б) и 11 (в) при значениях рз=0,995; 1,3; 1,7; 1,99 и фиксированном значении р4=0,462.

Наилучшее согласие с экспериментально измеренной функцией возбуждения ионов достигается при удвоенном значении параметра жесткости в потенциале взаимодействия СГ-Г, т.е. рз=1,99. При этом значении рз достигается практически количественное согласие экспериментально измеренной функции возбуждения с рассчитанной, представляющей сумму всех трех каналов, в которых образуются ионы Св+. Увеличение в 2 раза параметра р4 не оказывает заметного влияния на форму этой функции

возбуждения (см. рис. 5), что связано, по-видимому, с небольшим числом ведущих к образованию ионов Св+ конфигураций столкновения, вероятность реализации которых зависит от величины этого параметра.

Улучшение количественного согласования экспериментальной и рассчитанной функций возбуждения ионов Ш>+ потребовало удвоения параметров жесткости сразу в обоих парных потенциалах взаимодействия одноименных ионов; увеличение этого параметра только в одном из потенциалов оказалось неэффективным (см. рис. б.)

Следует отметить, что формы и положения максимумов функций возбуждения комплексных ионов СвСШЬ* и СзЖЬ+ не чувствительны к изменению значений параметров рз и р4, хотя абсолютные величины сечений меняются существенно.

0123456789 10 Энергия столкновения, эВ

О 1 2 3 4 5 б 7 8 9 10 Энергия столкновения,эВ

0123456789 10 Энергия столкновения, эВ

I 4

С5Т г/

р,=199,р4=0.924 #

Траекг. расчет ^»Эксперимент

0 1 2 3 4 5 6 7 8 9 10 Энергия столкновения, эВ

Рис. 5. Сравнение экспериментальной и полученных из траекторных расчетов функций возбуждения ионов Св+ для различных наборов параметров жесткости в потенциалах взаимодействия одноименно заряженных ионов.

На рисунке 7 показаны расчетные и измеренные в эксперименте в лабораторной системе координат угловые распределения атомных ионов Св+ и Шз+. Несмотря на некоторые количественные расхождения для ионов Ш>+ можно отметить хорошее совпадение угловых диаграмм. Это служит еще одним аргументом адекватности используемой ППЭ реальному взаимодействию молекул в условиях эксперимента, поскольку угловые распределения продуктов обладают высокой чувствительностью даже к небольшим изменениям топологии поверхности потенциальной энергии.

р1=0,995, р4=0 462

0123456789 10 Энергия столкновения, эВ

1 23456789 10 Энергия столкновения, эВ

123456789 10 Энергия столкновения, эВ

123456789 10 Энергия столкновения, эВ

Рис. 6. Сравнение экспериментальной и полученных из траекторных расчетов функций возбуждения ионов ЛЬ+ для различных наборов параметров жесткости в потенциалах взаимодействия одноименно заряженных ионов.

0 10 20 30 40 50 60 70 80 90 Угол рассеяния, град.

14

12

3 10

г 8

о

о 6

X

и 3- 4

и 2

0

_,Траект. расчет

^Эксперимент

"0 10 20 30 40 50 60 70 80 90 Угол рассеяния, град.

Рис. 7. Сравнение экспериментальных и полученных в траекторном счете угловых распределений ионов Сэ+ (а) и Шэ+ (б) в лабораторной системе координат при энергии столкновения 6,0 эВ.

Приведенные результаты показывают, что траекторная модель с удвоенными значениями рз и р4 практически количественно описывает реальное взаимодействие в основных каналах в системе СвСЛ + ИЫ и может служить основой для расчета динамических характеристик взаимодействия, не поддающихся прямому экспериментальному измерению.

Увеличение значений параметров рз и р4 соответствует более раннему отталкиванию взаимодействующих ионов одного знака и реализации отталкивания одноименных ионов на большем межъядерном расстоянии. Практически это означает, что при удвоенных значениях параметров жесткости взаимодействующие ионы не могут сблизиться до столь малых расстояний, как при начальных значениях этих параметров. Сравнение парных потенциалов взаимодействия СГ-Г и Сз+-ЯЬ+ при двух различных наборах параметров жесткости показано на рис. 8.

15

12

»

г?

е в- 6

в

о

3

0

р,=0,995

"О 2 4 6 8 10 12 14 16 18 20 Межатомное расстояние, а.е.

0 2 4 6 8 10 12 14 16 18 20 Межатомное расстояние, а е.

Рис. 8. Парные потенциалы взаимодействия СГ - Г (а) и Се* - Юэ+ (б) для начальных и удвоенных значений параметров жесткости рз и Р4.

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

Канал 2, как было отмечено выше, открывается при энергиях столкновения, превышающих сумму энергий диссоциации обеих молекул реагентов, равную примерно 9 эВ. Далее сечение увеличивается с ростом энергии столкновения, но для каждого значения энергии сечение тем выше, чем больше величина любого из параметров рз и р4. Объясняется это тем, что необходимое для диссоциации сильное внутреннее возбуждение обеих молекул за счет сил отталкивания одноименных ионов возникает при больших величинах межмолекулярных расстояний, и, в частности, при больших значениях прицельного параметра.

Как было показано ранее, сечение канала 3 заметно увеличивается с ростом параметра рз и, вместе с тем, уменьшается, хотя и в меньшей степени, при увеличении параметра р4. Сечение симметричного канала 4 оказывается практически не чувствительным к изменению рз, но увеличивается с ростом р4. Совместное действие этих факторов приводит к тому, что при удвоенных значениях обоих параметров рз и р4 рост сечений этих каналов с увеличением энергии столкновения происходит быстрее, и при Есп,.=25 эВ сечение канала 3 оказывается больше в 3,5 раза, а канала 4 - в 2 раза по сравнению с величинами сечений при начальных значениях параметров жесткости.

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

Интересные изменения с ростом параметров жесткости происходят в поведении функций возбуждения каналов 8-11, в которых образуются трехатомные ионные комплексы. Все функции возбуждения имеют пороговый характер, проходят через максимум и далее сечения уменьшаются с ростом энергии столкновения. Для всех исследованных наборов параметров жесткости функции возбуждения каналов 8 и 10 образования положительно заряженных ионных комплексов имеют максимум при одних и тех же значениях энергии столкновения около 4,0 эВ. При этом сечение канала 8 систематически уменьшается с ростом любого из параметров рз или р4, в то время как сечение симметричного канала 10 увеличивается с ростом рз и уменьшается при увеличении р4. В результате сечение в максимуме функции возбуждения канала 8 при удвоенных значениях р3 и р4 примерно в 7 раз меньше сечения для начального набора параметров жесткости, а для канала 10 - в 7 раз больше.

Иначе обстоит дело для каналов 9 и 11, в которых образуются положительные атомные ионы и отрицательно заряженные ионные комплексы. Сечение канала 9 заметно уменьшается с ростом рз и относительно столь же быстро растет при увеличении р4. При этом энергия, соответствующая положению максимума функции возбуждения, в первом случае увеличивается, а во втором случае максимум несколько смещается в сторону меньших значений энергии столкновения. В результате, при увеличении в 2 раза значений Рз и Р4 максимальное сечение канала 9 практически не изменилось, но положение максимума функции возбуждения сместилось с 4,0 эВ на 5,5 эВ. Как было отмечено выше, сечение канала 11 обладает очень высокой чувствительностью к значениям параметра рз и резко уменьшается при его увеличении. При этом положение максимума функции возбуждения смещается в сторону более высоких значений энергии столкновения. С увеличением значения параметра р4 сечение канала также уменьшается, хотя и в значительно меньшей степени, сохраняя при этом положение максимума функции возбуждения. В итоге, при удвоении значений рз и Р4 максимум функции возбуждения канала 11 сместился с 4,0 эВ на 9,0 эВ, а максимальное значение сечения уменьшилось почти в 30 раз.

Следует отметить, что при удвоении значений параметров жесткости значительно увеличилось сечение канала 10, в котором образуется наиболее легкий ион СГ и тяжелый комплекс СбЖЬ+. Сечение канала 9 образования среднего по массе иона Ш>+ практически не изменилось, а сечения каналов 8 и 11, в которых образуются тяжелые ионы Г и Се*, сильно уменьшились. Аналогичный эффект наблюдается и для группы каналов 3-6. В наибольшей степени здесь увеличилось сечение в максимуме функции возбуждения для канала 5, где образуются два самых лёгких иона и СГ. В меньшей степени увеличилось сечение для каналов 3 и 4, а сечение канала 6 образования тяжелых ионов Сз+ и Г и легкой молекулы ЮэС1 уменьшилось.

Наиболее драматичные различия в поведении функции возбуждения при изменении значений р3 и р4 наблюдаются для канала 7, причем значительно более сильное влияние оказывает параметр рз. Как видно из сравнения рисунков 9(а) и 9(в), для начального значения р3 = 0,995 даже удвоение р4 не меняет качественного поведения функции возбуждения, которая характеризуется монотонным уменьшением сечения с ростом энергии столкновения во всем просчитанном диапазоне. В то же время удвоение р3 при сохранении начального значения р4 = 0,462 качественно меняет вид функции возбуждения, которая приобретает колоколообразную форму с явно выраженным порогом (рисунок 9,6). Для начальных значений рз = 0,995 и р4 = 0,462 сечение этого канала было максимальным при близких к нулю значениях энергии и далее плавно уменьшалось с

ростом энергии столкновения вплоть до 25,0 эВ (при изменении Ео™ от 0,1 до 3,0 эВ сечение уменьшалось в 10 раз). При удвоенных значениях параметров жесткости канал открывается при Бол, =0,3 5 эВ, сечение достигает максимума при 3,3 эВ и далее уменьшается с ростом энергии столкновения.

Энергия столкновения, эВ Энергия столкновения, эВ

Рис. 9. Функции возбуждения канала 7 для начальных и удвоенных значений каждого из параметров жесткости рз и р4.

Следует отметить, что при удвоенных значениях рз и р4 даже в максимуме при Его, =3,3 эВ сечение этого канала в 3 раза меньше, чем при начальных значениях параметров жесткости для той же энергии столкновения. Изменение функции возбуждения от экспоненциально убывающей зависимости без порога к качественно иной форме с явно выраженным пороговым поведением свидетельствует об изменении основного механизма взаимодействия в этом канале. Качественные изменения формы функции возбуждения этого канала для ряда промежуточных значений параметров жесткости в потенциалах взаимодействия одноименных ионов приведены на рисунке 10.

Наблюдаемые зависимости функций возбуждения различных каналов взаимодействия в системе С&С1 + ЯЫ от значений параметров жесткости в потенциалах взаимодействия одноименных ионов определяются совместным действием динамических и кинематических факторов, определяющих основные механизмы взаимодействия в каждом из реализующихся каналов. Подробному рассмотрению этих механизмов на основе анализа отдельных траекторий посвящена глава 3 диссертации.

7 канал (р^О.924) СбС! + RbJ -> СвД + ЯЬС!

О 2 4 б 8 10 Энергия столкновения, эВ

100

80

В

я* 60

5

Й 40

Я

б 20

0-

7 канал (р,=1.7) СвС1 + ЯЫ -> Сб: + ЯЬС1

0 2 4 6 8 10 Энергия столкновения, эВ

Рис. 10, Функции возбуждения канала 7 при различных значениях рз (а) и р4 (б).

В главе 2 рассматривается динамика конверсии энергии столкновения во внутреннюю энергию молекул реагентов.

Процессы конверсии поступательной энергии столкновения молекул во внутреннюю энергию партнеров играют большую роль в механизмах активации реагентов при бимолекулярных превращениях и, в частности, в процессах индуцированной столкновениями диссоциации молекул на атомы и ионы. Эффективность передачи энергии, как показано в ряде исследований [11-13, 17, 18], сильно зависит от ориентации партнеров по столкновению, соотношения масс атомов, участвующих в реакции, параметров поверхности потенциальной энергии и ряда других динамических и кинематических особенностей процесса.

Хорошее совпадение экспериментально измеренных функций возбуждения образования атомных и молекулярных ионов, а также угловых распределений продуктов с результатами траекторного моделирования, как следует из предыдущей главы, дает основание использовать рассматриваемую ППЭ для определения характеристик взаимодействия, не поддающихся прямым измерениям. К таким характеристикам относятся, в частности, распределения продуктов при энергиях вблизи порога и под порогом, а также процессы передачи энергии Т —» V, Я и перераспределение энергии в двухатомной молекуле. Эти характеристики дают представление о динамике столкновений, сопровождающихся передачей энергии.

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

В следующей части главы приводятся результаты трехмерного траекгорного моделирования процессов внутримолекулярного перераспределения и межмолекулярного обмена энергией для различных значений энергии столкновения и начальных внутренних состояний молекул реагентов. Обнаружено, что эффективность конверсии поступательной энергии столкновения во внутреннюю энергию (Евн) молекул уменьшается с ростом Ео™, что можно видеть на рисунке 11, показывающем гистограммы распределений кинетической энергии разлета молекул для трех различных значений энергии столкновения. Так, для Ео™=0,2 эВ в 75% столкновений конечная энергия разлета молекул отличается не более чем на 10% от начальной энергии столкновения. Для Е(ГО1=1,0 эВ этот показатель составляет 86%, а при энергии столкновения 10,0 эВ разность конечной энергии разлета молекул и энергии столкновения находится в пределах 10% уже для 95,5% просчитанных траекторий.

CíCI♦ ИЫ -» CiC]

Т^-ЮООК

T_,-IOOOK

Праду]

Ж.

0.12 0.14 0.16 0.19 0J0 0,22 0.24 0,26

Энергия, эВ

в;

& 1001

90

§ S0'

о 70

5 60

я т S0'

8 40

5 30

20'

1 10'

| 0.4

Csd + Rbj->C|CI T^-IOOOK т -юоок

Продукта

J k

0.6 0,8 1,0 1.2 Энергия, эВ

£-100.

s 50 § »:

в 70

5 60

§ JO'

u 40'

I 30-

á 20-

cici»Rbj->ceCi * ru T^-IOOOK ■1000К

Реагенты

Энергия, эВ

Рис. 11. Распределения кинетической энергии разлета молекул в конце траектории при температурах источников молекулярных пучков ЮООК для трех фиксированных значений начальной энергии столкновения - ОД эВ (а), 1,0 эВ (б) и 10,0 эВ (в).

Для более подробного исследования особенностей Т V, R обмена энергии расчеты проведены для заданного начального колебательного (V,nl) и вращательного (J'"') состояния молекул в широком диапазоне энергий столкновения. Следует отметить, что конечные внутренние состояния партнеров столкновения определялись по формулам классической механики через радиальную и тангенциальную составляющие импульса относительного движения ионов в молекуле в предположении независимости колебательного и вращательного состояний. На рисунке 12 показаны распределения конечных колебательных квантовых чисел Vcscifin молекул CsCl после столкновения для Ест, =0,2 эВ и начальных состояний молекул с VCsci"" = VRbJinl = 15, Jcsci"" = W' = 0 (а) и Vcsci'"1 = VRbj"" = 15, Jcsci'"' = Wm = 150 (б). В случае отсутствия начального вращательного возбуждения основная часть молекул, а точнее 72,5%, сохраняют значение начального колебательного квантового числа и еще 7,4% молекул оказываются на соседних уровнях, отличающихся на единицу от Vcsci"11. Совершенно иная ситуация наблюдается в случае Jcsci"" =150. Здесь лишь для 17% молекул Vcsci"" = VCsa,n' =15, еще 13,5% имеют Veso™ =14 и 52% столкновений приводят к увеличению колебательной энергии молекул CsCl до уровней Vcscif'n =16 и Vcscifin =17. Следует отметить, что в обоих случаях возбуждаются достаточно высокие колебательные состояния вплоть до Vcsc/™ =30, однако с ростом вращательной энергии заселенность высоких колебательных уровней с Vcscifin >20 значительно снижается. Возбуждение столь высоких колебательных уровней при низкой энергии столкновения и Jcsci"" =0 возможно лишь при эффективном V—>V обмене энергии между партнерами столкновения, т.е. когда одна из молекул теряет большую часть своей колебательной энергии.

СвС1 +

т"

12

-> сэс! + ЯЬ; Ет=0,2 эВ

>с,с""=15 ;с,с'",= 0

IV "'=15 Г,""=0

яы

лы

16 20 24 Колебательное квантовое число

100 8060 40 20 0

СяС!

Ш -> СвС! + ЯЫ Е =0,2 эВ

оти '

уе*> 15 '«> 150

V, "=15 ■!„,.,""= 150

12 16 20 24

Колебательное квантовое число

Рис. 12. Распределения колебательных квантовых чисел молекул СвС1 после столкновения для Бот =0,2 эВ и начальных состояний молекул с Х^а"11 = УаЫт' =15, ^а'"' = ■Гш™ = 0 (а) и УС5сГ = = 15, Зааш = :КьУп1 =150 (б).

С ростом энергии столкновения доля молекул, оставшихся на прежнем колебательном уровне, увеличивается, что свидетельствует о снижении эффективности Т—»V обмена. Соответствующие распределения для энергий столкновения 1,0 эВ и 5,0 эВ приведены на рисунке 13. Доли молекул, сохранивших начальное колебательное состояние, составляют соответственно 79% и 84%.

Колебательное квантовое число Колебательное квантовое число

Рис. 13. Распределения колебательных квантовых чисел молекул СэС1 после столкновения для начальных состояний молекул с УЪа"" = Уи,Уш = 15, ^Ба'"' = •'яы'"' = 0 при энергиях столкновения 1,0 эВ (а) и 5,0 эВ (б).

Обнаружено, что эффективность внутримолекулярного V—>Я обмена энергии весьма низка. На рисунке 14 показаны распределения вращательных квантовых чисел молекул СвС1 после столкновения для Ео™ =0,2 эВ и различных начальных колебательных состояний реагентов. Видно, что на начальном нулевом вращательном уровне после столкновения молекул нет вообще. Конечные вращательные квантовые числа имеют равновесное распределение, причем из сравнения рисунков 14 (а) и 14 (б) можно сделать вывод о том, что начальное колебательное возбуждение не оказывает сильного влияния на это распределение, лишь в незначительной степени смещая максимум в сторону более высоких значений 1с8С1й"-

Вращательное квантовое число Вращательное квантовое число

Рис. 14. Распределения вращательных квантовых чисел молекул CsCl после столкновения для Еот„,=0,2 эВ и начальных состояний молекул с Vcsci"" = Jcsci""= "V^j1"' = JRbj"" = 0 (а) и Vcsaini - VRbJini = 15, JCsCiini = = 0 (б).

При значительном начальном вращательном возбуждении молекул с Jcsci'"' =150 максимумы конечных распределений Jcsci"" совпадают с начальным состоянием молекул. Гистограммы соответствующих распределений для трех различных значений энергии столкновения приведены на рисунке 15. Видно, что эффективность Т—>R обмена увеличивается с ростом энергии столкновения. Для Еотн=0,2 эВ около 16% молекул сохраняют начальное вращательное состояние с Jcsci"1' =150, для энергий 1,0 эВ и 5,0 эВ этот показатель составляет соответственно 4,7% и 3,8%.

1 80 g

5 60'

ж г

о 40-

CsCl + RbJ -> CsCl + RM Em=0,2 зВ

150

-o

■о J..

■ ISO

0 40 80 [20 160 200

Вращательное квантовое число

Вращательное квантовое число

Рис. 15. Распределения вращательных квантовых чисел молекул CsCl после столкновения для начальных состояний молекул с Vcsci'"'= Vru'"1 = 0, Jcsci"" = Jrw'"' = 150 при энергиях столкновения 0,2 эВ (а), 1,0 эВ (б) и 5,0 эВ (в).

Глава 3 посвящена исследованию динамики появления и реализации каналов химического рассеяния в системе CsCl + RbJ.

По характеру протекающих процессов и набору образующихся продуктов все каналы химического рассеяния можно разделить на три группы. Первая группа включает каналы 2-4, в которых реализуется типичная столкновительно-индуцированная диссоциация с распадом на ионы одной или обеих сталкивающихся молекул. Ко второй группе следует отнести каналы 5-7, в которых наряду с диссоциацией обоих партнеров столкновения происходит образование новых молекул. И, наконец, третья группа каналов с 8 по 11 ведет к образованию трехатомных ионных комплексов.

Функции возбуждения каналов первой группы имеют форму, типичную для процессов СИД, и характеризуются размытым порогом и быстрым ростом сечения при

увеличении энергии столкновения выше порога с последующим насыщением. В качестве примера на рисунке 16 показана функция возбуждения канала 2 в диапазоне энергий столкновения от 0 до 100,0 эВ и более подробно вблизи порога.

Для каналов 3 и 4 характерно достижение максимумов функциями возбуждения при энергиях столкновения соответственно 33,0 эВ и 25,0 эВ и дальнейшее плавное уменьшение сечения с ростом энергии столкновения.

Обнаружено, что для всех трех каналов наиболее предпочтительными являются отличные от нуля значения прицельных параметров (Ь), а максимальные их величины несколько увеличиваются с ростом энергии столкновения. Канал 2 в широком диапазоне энергий столкновения с максимальной вероятностью реализуется в диапазоне прицельных параметров от 1 до 6 а.е., функции прозрачности канала 3 имеют максимум при Ь=8 а.е., а для канала 4 наиболее характерные значения прицельных параметров лежат в интервале 2-10 а.е.

СзС1 + Ш -> Сз* + С1' + Ю>* + Г 2 канал

"0 20 40 60 80 100 Энергия столкновения, эВ

3

5 б

0 и

1 *

и т

а 2

о

СэС1 + Ш -> С5* + С1 + ш>* + г , 2 канал

8 10 12 14

Энергия столкновения, эВ

Рис. 16. Функция возбуждения канала 2 в диапазоне энергий столкновения от 0 до 100,0 эВ (а) и вблизи порога (б).

Анализ распределений внутренних состояний молекул реагентов в каждом из рассматриваемых каналов показывает, что при близких к порогу энергиях столкновения эти распределения характеризуются более высокими температурами по сравнению со значением Т=1000К, выбранным для розыгрыша начальных состояний всей совокупности молекул, охватывающей траектории полного набора реализующихся при этой энергии каналов. С ростом Е<т, температуры распределений уменьшаются. При этом для канала 2 распределения внутренней энергии СвО и ЯЫ соответствуют примерно одинаковым значениям температур (Тасг1273К и Тш =1255К при Ео™=10,0 эВ и уменьшаются до значений Тс5сгЮ75К и Тад =105бК при Еотн =13,0 эВ), а для каналов 3 и 4 более высокая температура распределения внутренней энергии характерна для той молекулы, которая в результате столкновения испытывает распад на ионы. Так, в канале 3 при Еотн =4,5 эВ распределения колебательной энергии реагентов имеют температуры ТЪсг! 575К и Т1Ш=1302К, уменьшаясь до значений Тйа=1158К и ТцЫ=1047К уже при энергии столкновения 5,0 эВ. Для Бот, =7,0 эВ температуры составляют Тс5сгЮ55К и Тш=1019К. Для канала 4 более высокая температура характерна для распределения внутренней энергии молекул ЯЫ.

Показано, что для всех трех каналов подавляющее большинство траекторий соответствует механизму СИД через преимущественное взаимодействие в паре Сб* - Г. Именно это межионное расстояние оказывается минимальным среди попарных расстояний принадлежащих разным молекулам реагентов ионов. Канал 2 характеризуется прямым механизмом при всех энергиях столкновения, т.е. после сближения реагентов и их диссоциацией образующиеся продукты реакции сразу покидают зону взаимодействия.

В то же время в каналах 3 и 4 при близких в порогу энергиях столкновения значительное число траекторий реализуется через образование промежуточного комплекса молекулы-снаряда с более тяжелым ионом диссоциирующей молекулы. Время жизни такого промежуточного состояния составляет, как правило, 1,2-10'1 -3,6-10"2 с, по прошествии которого происходит распад комплекса на конечные продукты.

Динамика реализующихся в каналах 5-7 взаимодействий до настоящего времени в литературе не рассматривалась. Набор продуктов каналов 5 и 6, аналогично рассмотренным выше каналам 3 и 4, составляют три частицы - два атомных иона и молекула. Вследствие этого функции возбуждения всех 4 каналов качественно идентичны и характеризуются размытым пороговым поведением, а энергия появления каждого канала определяется его эндоэргичностью. Максимальные значения сечений каналов 5 и 6 достигаются при меньших по сравнению с каналами 3 и 4 энергиях столкновения, равных соответственно 11,0 эВ и 17,0 эВ, и быстрее уменьшаются при дальнейшем росте Em,. Интересно отметить, что в отличие от первой группы каналов в широком диапазоне энергий столкновения каналы 5 и 7 реализуются с наибольшей вероятностью при близких к нулю значениях прицельного параметра. Форма функции прозрачности канала 6 заметно меняется с ростом энергии столкновения, и положение ее максимума с Ь=7 а.е. при Еотн =5,0 эВ смещается на Ь=4 а.е. при Е<л„=13,0 эВ и соответствует близким к нулю значениям прицельного параметра при Ецт =25,0 эВ (см. рисунок 17).

1,2

1.0

d 0,8

ё 0,6 о*

1 0,4 ° 0,2 0,0

CsCI + RbJ ■> RbCl + Cs* + J" 6 канал

\

уЕ„-25,0зВ

V- E -13,0 эВ

. \ /

\ .

=7,0 3B-V. E =5,0'B

0 4 8 12 16 20 Прицельный параметр, а.е.

Рис. 17. Функции прозрачности канала 6 при энергиях столкновения 5,0 эВ, 7,0 эВ, 13,0 эВ и 25,0 эВ.

Распределения кинетической энергии разлета продуктов канала 5 относительно центра масс системы при энергии столкновения 7,0 эВ приведены на рисунке 18. Как следует из рисунка, минимальной скоростью разлета относительно центра масс обладает молекула Се!

Представленные на рисунке 18 гистограммы распределений энергий разлета и угловые распределения продуктов, характеризующиеся преимущественным рассеянием ионов Ю>+ вперед (за положительное направление принимается направление первоначального вектора скорости молекул СвС!), молекул в заднюю полусферу и ионов СГ в широкий диапазон углов, находятся в соответствии с полученными из анализа отдельных траекторий возможными механизмами взаимодействия реагентов. Почти во всех рассмотренных траекториях в момент наибольшего сближения молекул минимальными оказываются расстояния в парах Сэ+ - Г и Сэ+ - ЯЬ+. За счет полученной энергии молекула ЯЫ диссоциирует и, благодаря отталкиванию от более тяжелого иона Се*, освободившийся ион ЯЬ+ рассеивается вперед, т.е. в сторону, противоположную направлению своего первоначального движения. При этом молекула с невысокой скоростью рассеивается назад.

£>00, О

о .. § 60-М

СаС) + ЯЫ -> С*) + С!" + ЯЬ'

Кинетическая энергия иона СГ относительно центре масс

I4 100,

С|С1 + ЯЫ -> си + СГ + ЯЬ'

3 канал Кинетическая энергия ЯЬ относительно центра масс

к"т

-7,0 эВ

\

Энергия, эВ

Энергия, эВ

СяС1 + ЯЫ •> Са + СГ * 5 какал

Кинетическая энергия молекулы С^ относительно

Энергия, эВ

Рис. 18. Распределение кинетической энергии разлета продуктов канала 5 относительно центра масс системы при энергии столкновения 7,0 эВ: СГ (а), ЯЬ+ (б) и Св.! (в).

Даже при небольшой энергии столкновения в 5,0 эВ более половины всех траекторий носит прямой характер, т.е. после сближения реагентов и их взаимодействия образовавшиеся продукты сразу удаляются друг от друга. Существует также небольшое количество траекторий, в которых при низкой кинетической энергии разлета за счет кулоновского притяжения один из образовавшихся ионов остается некоторое время вблизи молекулы, причем это может быть как ион СГ, так и ион ЯЬ+.

С ростом энергии столкновения преимущественный механизм взаимодействия молекулы ЯЫ с ионом Сб+ не меняется, но образовавшаяся молекула обладает намного более высокой внутренней энергией, причем в большей степени увеличивается ее колебательная составляющая. На рисунке 19 показаны распределения колебательной энергии молекул при энергиях столкновения 4,5 эВ, 7,0 эВ и 13,0 эВ. При Ест, =13,0 эВ имеется значительное число траекторий, в которых молекула Сз.1 обладает колебательной энергией вблизи порога диссоциации молекулы.

£ 10Оп

1 £

СаС| + !1Ы->С!Л-С1Ч 5 канал Е =4,5 эВ

Колебательная энергия молекулы С&Л

О.,=4,0320 эВ

ч

I

Энергия, эВ

Энергия, эВ

Рис. 19. Распределения колебательной энергии молекул Сэ1 в конце траектории канала 5 при энергиях столкновения 4,5 эВ (а), 7,0 эВ (б) и 13,0 эВ (в).

Канал б, в котором образующаяся молекула ЯЬС1 немного уступает по массе обоим ионам Сэ+ и Г, имеет больше различных механизмов реализации, меняющихся, кроме того, с ростом энергии столкновения.

При низких энергиях столкновения большая часть траекторий реализуется через промежуточный комплекс ЯЫСГ со временем жизни 1,2'10" -5-10"'2 с. Комплекс образуется через механизм, близкий к механизму "срыва", т.е. при пролете молекул с большим прицельным параметром молекула ЯЫ за счет притяжения в паре ЯЬ+ - СГ отрывает ион СГ. При этом освободившийся ион Сэ+ продолжает движение в первоначальном направлении, а обладающий избыточной энергией комплекс распадается

на молекулу 11ЬС1 и ион Г, причем оба эти близкие по массе продукта обладают небольшой кинетической энергией относительно центра масс.

Значительно меньшее число траекторий реализуется симметричным описанному выше механизмом, когда вначале происходит диссоциация молекулы 1Ш с освобождением иона Г и образованием промежуточного комплекса СвСШЬ* с более длинным, как правило, временем жизни, который впоследствии распадается на ион Сб+ и молекулу ЯЬС1. Существует также небольшое число траекторий, взаимодействие в которых имеет прямой механизм, когда образование новой молекулы ЯЬС1 происходит одновременно с распадом обеих молекул реагентов. При этом механизме оба освободившихся тяжелых иона Сэ+ и Г продолжают движение в первоначальных направлениях, а молекула ЮэС1 рассеивается в широкий диапазон углов.

С ростом энергии столкновения все траектории становятся прямыми, а преимущественный механизм взаимодействия меняется. В значительной части траекторий обе молекулы диссоциируют сразу после сближения, причем молекула СзС1 распадается в большинстве случаев несколько раньше, чем 11Ы. При этом расстояние в паре Се* -1' на некоторое время становится близким к равновесному и оба иона движутся в направлении первоначального вектора скорости молекулы СэО, т.е. вперед. Образовавшаяся молекула Ш)С1 рассеивается в направлении первоначального движения более тяжелого из составляющих ее ионов рубидия, т.е. в заднюю полусферу. Однако, полная энергия пары Сэ+ - Г положительна, вследствие чего ионы начинают удаляться друг от друга, соответствующее попарное расстояние меняется медленнее, чем расстояния от этих ионов до молекулы ЛЬС1.

Следует отметить, что при всех энергиях столкновения молекулы ЯЬС1 обладают низкой вращательной энергией, не превышающей 1,0 эВ даже при £„„ =13,0 эВ. В то же время распределения колебательной энергии меняются драматическим образом, и при энергиях выше 10,0 эВ заселенность колебательных квантовых чисел становится инверсной, а значительная часть молекул ЯЬС1 обладают колебательной энергией, близкой к энергии диссоциации молекулы (см. рисунок 20).

С^ + Ш-ЖЬСКСв' + Г 6 канал Е =4.5 ЗВ

Колебательная энергия молекулы КЬС1

Энергия, эВ

Энергия, эВ

Рис. 20. Распределение колебательной энергии молекул ЯЬС1 в конце траектории канала 6 при энергиях столкновения 4,5 эВ (а), 7,0 эВ (б) и 13,0 эВ (в).

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

Канал является почти термонейтральным, и, как отмечено выше, открывается при энергии 0,35 эВ, далее сечение увеличивается с ростом Е0тн. до 3,3 эВ, и уменьшается при дальнейшем увеличении энергии столкновения. Продукты этого канала рассеиваются в широкий диапазон углов при всех значениях Ео™., однако направления

Энергия, эВ

СзС! + 1?Ы -> ЯЬС) + Са* + Г

6 канал п "4.7357,1 E_.-13.05B \

преимущественного рассеяния изменяются неожиданным образом с ростом энергии столкновения. Расчеты показывают, что при £^,=1,0 эВ образующиеся молекулы С^ рассеиваются преимущественно в заднюю полусферу, т.е. в сторону, противоположную первоначальному направлению движения молекул СэС1 (см. рисунок 21,а). В диапазоне энергий столкновения 2,0-4,0 эВ - симметрично вперед и назад и в меньшей степени в стороны (рисунок 21,6). При Ео™ =5,0 эВ большая часть молекул Се,! начинает рассеиваться вперед, тенденция эта усиливается с ростом энергии до 7,0 эВ, а потом опять ослабевает при дальнейшем увеличении энергии столкновения до 10,0 эВ (см. рисунок 21,в-д). При энергиях столкновения выше 10,0 эВ образующиеся молекулы Сз1 рассеиваются преимущественно назад (рисунок 21,е). Изменение угловых гистограмм рассеяния продуктов свидетельствует об изменении преимущественного механизма взаимодействия в этом канале с ростом энергии столкновения.

Расчеты показывают, что при энергии столкновения 1,0 эВ все траектории реализуются через промежуточный комплекс со временем жизни от 1,2-10"12 до 1,2-10 с. В большинстве траекторий взаимодействие начинается с разрыва связи в молекуле ЯЫ и формирования новой связи между ионами Шз+ и СГ. Образующийся комплекс Сз+--СГ~ Шз+ имеет линейную структуру с расстояниями в парах С$+ - СГ и СГ - ШГ, близкими к равновесным. Ион Г вращается вокруг комплекса на расстоянии 20-30 а.е., после чего происходит перераспределение энергии, ион Сэ+ покидает комплекс и образуется молекула С&1.

С ростом энергии столкновения комплексный механизм взаимодействия меняется на прямой и при энергии 3,0 - 4,0 эВ около половины траекторий имеют прямой характер, а при Ео™=5,0 эВ и выше практически все траектории получаются прямыми. Механизм взаимодействия при этом также может оказаться различным. Примерно в половине траекторий при Ео™ =5,0 эВ механизм аналогичен рассмотренному выше и характеризуется первоначальным образованием прочной связи между сблизившимися на расстояние, близкое к равновесному, ионами ЯЬ+ и СГ. Однако теперь одновременно происходит и диссоциация обеих молекул реагентов с последующим формированием молекулы Сб-Г, которая при таком механизме оказывается, как правило, весьма сильно возбужденной. Значительно меньшая часть траекторий при энергии столкновения 5,0 эВ характеризуется минимальным расстоянием между сталкивающимися молекулами в паре - Г. Как и ранее, одновременно с образованием этой новой связи происходит диссоциация обеих сталкивающихся молекул, но теперь уже образующаяся последней молекула ШэС1 обладает более высокой внутренней энергией.

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

Следует отметить, что обе образующиеся в канале 7 молекулы продуктов CsJ и Ш)С1 обладают исключительно высокой внутренней энергией, реализующейся преимущественно на колебательных степенях свободы. При этом распределения колебательной энергии молекул СбЗ могут бьпъ аппроксимированы равновесными лишь до энергии столкновения 7,0 эВ, а при более высоких значениях Е»™ распределения становятся бимодальными, и значительное число молекул находится в максимально возможных колебательных состояниях с энергией, близкой к энергии связи в молекуле. Распределение колебательной энергии молекул ЯЬС1 становится инверсным уже при Еотн.=5,0 эВ, причем наиболее вероятное значение колебательной энергии находится в диапазоне 3,8-4,5 эВ. При дальнейшем росте энергии столкновения гистограммы распределения колебательной энергии качественно не меняются и имеют максимум при Екол.=4.5 эВ.

СбС1 + Ш -> С51 + ЯЬС1 Д 7 канал

30 60 90 120 150 180 Угол рассеяния, град.

5 4

"< 3

£ 8 2

1

СзС1 + КЫ->С51 + ЯЬС1 7 канал Ет=17,0эВ

, ЯЬС1

30 60 90 120 150 Угол рассеяния, град.

180

Рис. 21. Угловые распределения образующихся в канале 7 молекул С&1 и ЯЬС1 при энергиях столкновения 1,0 эВ (а), 3,0 эВ (б), 5,0 эВ (в), 7,0 эВ (г), 10,0 эВ (д) и 17,0 эВ (е).

Третью группу каналов с 8 по 11 объединяет то, что в каждом из них образуется всего лишь 2 частицы - атомный ион и трехатомный ионный комплекс. Все каналы этой группы открываются при энергиях столкновения 1,5-3,0 эВ, что существенно ниже энергии диссоциации каждого из реагентов. Это возможно только в том случае, когда диссоциация одной из молекул реагентов происходит одновременно с образованием новой связи освобождающегося иона с входящим в состав второй молекулы ионом

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

Расчеты показывают, что в подавляющем большинстве траекторий канала 8 минимальным расстоянием по ходу траектории является расстояние в паре СГ - ЛЬ*. В совокупности с угловыми распределениями продуктов взаимодействия, показывающими при всех энергиях столкновения в системе центра масс рассеяние ионного комплекса вперед, т.е. по направлению первоначального движения молекул СбС1, а образовавшегося иона Г - назад, и достаточно большими прицельными параметрами это позволяет сделать вывод о том, что в этом канапе реализуется так называемый "срывной" механизм. Молекулы СбС1 и ЯЫ при сближении ориентированы друг к другу ионами СГ и ЯЬ+, за счет сильного кулоновского притяжения первая молекула отрывает ион ЯЬ+ и образующийся комплекс продолжает двигаться в направлении первоначального вектора скорости молекулы СэСГ

В симметричном канале 10, в котором также образуется положительно заряженный ионный комплекс и отрицательный атомный ион, происходит изменение основного механизма взаимодействия с ростом энергии столкновения. При Еотн<5,0 эВ механизм аналогичен реализующемуся в канале 8, взаимодействие происходит при больших прицельных параметрах и образующийся атомный ион рассеивается в направлении своего первоначального движения. При дальнейшем увеличении энергии столкновения преимущественный механизм взаимодействия постепенно меняется со "срывного" на "рикошетный", наиболее предпочтительные значения прицельных параметров заметно уменьшаются, а направления рассеяния продуктов меняются на противоположные.

При всех энергиях столкновения "рикошетный" механизм взаимодействия является основным для канала 9. В подавляющем большинстве траекторий этого канала при сближении молекул реагентов наиболее сильное взаимодействие реализуется в парах Се* -Ш>+ и Св+ - Г, вследствие чего молекула ЯЫ диссоциирует, и за счет отталкивания в паре Сэ+ - ЯЬ+ образующийся ион ЯЬ+ рассеивается в направлении, противоположном вектору начальной скорости молекул ЯЫ, т.е. в переднюю полусферу. Одновременно формируется новая связь Сб+ - Г и образующийся ионный комплекс рассеивается назад по отношению к первоначальному направлению движения молекулы СэС1.

Характер взаимодействия реагентов в канале 11 для всех энергий столкновения включает два равновероятных механизма взаимодействия. Примерно половина траекторий соответствует "срывному" механизму, когда минимальным расстоянием между сближающимися молекулами оказывается расстояние в паре ЯЬ+ - СГ, и формирование этой новой связи происходит одновременно с разрывом связи Св+ - СГ в исходной молекуле реагента. При этом освободившийся ион Сэ+ сохраняет первоначальное направление своего движения, т.е. рассеивается вперед, а ионный комплекс рассеивается в заднюю полусферу. В другой части траекторий реализуется иной механизм взаимодействия, при котором распад молекулы СвС1 происходит вследствие притяжения в паре Сэ+ - Г. Однако, связь в молекуле ЯЫ остается не возбужденной и обладающий высокой кинетической энергией тяжелый ион Сэ+ пролетает на небольшом расстоянии от йода и продолжает движение в первоначальном направлении. Ион СГ присоединяется к молекуле ЯЫ, при этом новая связь ЯЬ+ - СГ в начальный период жизни комплекса оказывается весьма сильно возбужденной.

На рисунке 22 приведены гистограммы распределений кинетической энергии разлета продуктов в каналах 8-11 при энергии столкновения 5,0 эВ. Видно, что наиболее вероятное значение энергии разлета максимально для канала 8, откуда следует, что образующийся в нем ионный комплекс обладает минимальной внутренней энергией. Это соответствует отмеченным выше представлениям о "срывном" механизме взаимодействия в этом канале, при котором существующая до столкновения молекул связь в ионном комплексе остается не возбужденной и после взаимодействия. В трех других каналах сочетаются различные механизмы, при которых как "старая", так и вновь образующаяся связь в ионном комплексе оказываются весьма сильно возбужденными.

100

5

! 80-

г

I 60-

8 40-

в

го-

С5С1 + Ш -> С$С1ЮГ + У 8 канал Энергия разлета продуктов

\

(С«С1ИЬ + Г)

Е -5,0 эВ

ОТО

2 3 4 5 Энергия, эВ

СвС1 + Ш -> СзСи' + ЛЬ* 9 канал Энергия разлета продукте] / (СбС1Г + ЯЬ*)

1)11",-

Е -5,0 эВ \

б

2 3 4 5 Энергия, эВ

С5С1 + Ш->СзаьГ+ СГ

10 канал Энергия разлета продуктов / (CsR.br+ СГ)

Ет-5,0эВ

1 2 3 4 5

Энергия, эВ

100

80

60

§ 40

Л

I 20

СзС1 + Ш -> ШС1' + Сэ* 11 канал

Энергия разлета продуктов ' (ЯЫС1 + Сз*)

Е =5,0 эВ

1 2 3 4 5 Энергия, эВ

Рис. 22. Гистограммы распределений кинетической энергии разлета продуктов взаимодействия в каналах: 8 (а), 9 (б), 10 (в) и 11 (г) при энергии столкновения 5,0 эВ.

В главе 4 исследуется роль внутренней энергии реагентов и влияние соотношения масс сталкивающихся молекул на примере динамики взаимодействия в симметричных системах СвС1 + Ш и Св.! + Ш>С1.

Одной из наиболее важных проблем динамики элементарных процессов является определение роли внутренней энергии реагентов. Исследование этого вопроса в процессах СИД для трехатомных систем [3, 6,10] показало, что роль внутренней энергии реагентов значительна только при небольших, близких к порогу энергиях столкновения. Также было обнаружено, что для этих систем не существует сколько-нибудь заметного разделения внутренних степеней свободы по эффективности воздействия и сечение процесса определяется суммарной внутренней энергией молекул реагентов.

Аналогичные расчеты выполнены для системы СвС1 + ЯЫ при относительных энергиях столкновения 4,0 эВ и 17,0 эВ. Розыгрыш начальных внутренних состояний молекул реагентов при температурах сопел источников пучков 1000К, что соответствует условиям эксперимента в скрещенных молекулярных пучках, показывает, что реально максимально достижимые значения колебательных (V) и вращательных (I) квантовых чисел составляют:

УтмС5С| = 20, 1т„С5С, = 300, Утахкы = 30, :тикы = 400

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

Расчеты показывают, что при энергии столкновения 4,0 эВ и невысоком внутреннем возбуждении реагентов (внутренняя энергия ~ 0,13 эВ) колебательное возбуждение оказывается заметно более существенным по сравнению с вращением молекул для всех реализующихся при этой энергии каналов, в то время как для энергии 17,0 эВ столь не высокая внутренняя энергия не оказывает существенного влияния на сечения. При высоких начальных значениях V и 1, соответствующих внутренней энергии 1,5 эВ, эта тенденция еще заметнее при энергии столкновения 4,0 эВ для каналов 4-11, в то время как для канала 3 предпочтительнее оказывается вращение любой из молекул реагентов. Для энергии столкновения 17,0 эВ явно более чувствительными к высокому колебательному возбуждению оказываются каналы 3 и 6, остальные каналы не обнаруживают преимущественного роста сечения в соответствии с возбуждением колебательных или вращательных степеней свободы партнеров столкновения.

Таким образом, четырехатомная система показывает более сложную зависимость сечения отдельных каналов от внутренней энергии молекул реагентов по сравнению с трехатомными системами, причем роль колебательных и вращательных степеней свободы для разных каналов неодинакова. Для более подробного исследования влияния внутренней энергии реагентов на динамику взаимодействия в системе СзС1 + ЯЫ траекторные расчеты были выполнены для энергий столкновения 4,0 эВ, 6,0 эВ, 13,0 эВ и 25,0 эВ при различных температурах выходной щели источников молекулярных пучков. Значения температур составляли 1000К, а также 300К и 3000К каждого источника пучка в отдельности, и, кроме того, расчеты были выполнены для случая, когда температуры выходных щелей обоих источников пучков равны 1000К, но с добавкой 0,172341 эВ к относительной энергии столкновения, что соответствует увеличению температуры пучка на 2000К.

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

Сечения каналов 2-6 уменьшаются при понижении любой из температур источников пучков молекул реагентов и растут при их увеличении. При этом сечение канала 3 увеличивается в большей степени с ростом внутренней энергии молекул СвО, а сечения остальных каналов этой группы - при увеличении внутренней энергии молекул КЫ.

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

внутренней энергии любой из молекул реагентов вызывало рост сечения этого канала, причем более эффективным оказалось возбуждение ИЫ. Так, при Еотн =4,0 эВ повышение Тс»с1 с 1000К до 3000К привело к росту сечения на 10,6%, в то время как аналогичное увеличение Тш вызвало рост сечения на 26,1%. При увеличении энергии столкновения до 6,0 эВ различие в эффехтивностях начального возбуждения молекул реагентов становится еще более заметным и соответствующие цифры роста сечения составляют 7,5% и 28,7%.

Среди каналов образования трехатомных ионных комплексов наибольшую чувствительность к внутренней энергии реагентов обнаруживает канал 8. При относительной энергии столкновения 4 эВ сечение этого канала увеличивается на треть при снижении температуры источника пучка молекул СзС1 с 1000К до 300К и более чем в два раза при увеличении температуры источника пучка 1Ш с 1000К до ЗОООК. Интересно, что сечение симметричного канала 10 при аналогичных изменениях температур источников пучков исходных молекул меняется менее чем на 10%. Сечения каналов 9 и 11 образования отрицательно заряженных ионных комплексов обнаруживают зависимость от температур источников пучков реагентов лишь при низких значениях энергии столкновения, причем снижение любой из температур приводит к уменьшению сечений, а рост температур - к их увеличению.

Во второй части главы 4 сравниваются результаты траекторного моделирования динамики взаимодействия в симметричных четырехатомных системах СгС1 + ЮхГ (система 1) и Сб.1 + ЯЬС1 (система 2). Обе системы включают одинаковый набор ионов и управляются, следовательно, одной и той же ППЭ, но система 2 отличается от исходной системы эффективностью передачи энергии вследствие иного соотношения масс партнеров столкновения.

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

Массы молекул реагентов в исходной системе отличались в 1,26 раза, в то время как молекула С^ тяжелее ЯЬС1 в 2,15 раза. Можно ожидать, что сечения каналов, в которых происходит диссоциация тяжелой молекулы Се.!, уменьшатся, а связанные с диссоциацией легкой молекулы ЮзС1 каналы будут реализовываться с большей вероятностью.

Как было показано ранее, взаимодействие двух молекул с ионной связью может протекать по 12 каналам, которые, сохраняя для обеих систем симметрию по характеру взаимодействия в каждом канале, могут быть записаны следующим образом:

1. СзС1 + Ш СвСИ-КЫ

СбЛ-ЮЗС! -> СвЛ-КЬС!

2.

3.

4.

5.

6.

7.

8.

9.

10. И. 12.

-> Сз+ + С1' + Ш)+ + Г

Сз+ + СГ + Ш -» ЯЬ+ + Г + СзС1 Сз1 + ЯЬ+ + СГ -> ЯЬС1 + Сэ+ + г -> Сэ1 + ЯЬО

СбС1ШЭ+ + У -> СэСи' + ЯЬ+ -* СбЖЬ* + СГ ЯЫСГ + Сэ+

-> сваш

С5+ + СГ + ЯЬ+ + Г + .Г + ЛЬС1 -» + СГ + СБ] -> СэС1 + Юэ+ +Г -> Ш + Св+ + СГ СэС1 + Ш

->. сэжь*+сг

->• СвГСГ + ИЬ*

свст^+г яьси" + Сэ+ СзЖЬС!

Из этой схемы видно, что с точки зрения образующихся продуктов каналы 2, 9,11 и 12 для систем 1 и 2 совпадают, в то время как остальные каналы изменили своё положение. Далее анализируются различия функций возбуждения в симметричных каналах для систем 1 и 2.

Канал 2, в котором происходит образование четырех ионов, имеет для обеих систем примерно одинаковый порог около 9 эВ, равный сумме энергий диссоциации молекул реагентов. Функции возбуждения для обеих систем демонстрируют сходное поведение, но во всем просчитанном диапазоне энергий столкновения сечения этого канала для взаимодействия + ЯЬС1 несколько выше (см. рисунок 23).

Энергия столкновения, эВ

Рис. 23. Функции возбуждения канала 2 для взаимодействия в системах СвС1 + 1Ш (1) и Сз.Г + ЯЬС1 (2).

Каналы 3 и 4 соответствуют диссоциации на ионы первой или второй молекулы реагентов соответственно. Для системы 1 сечение канала 3 (диссоциация более легкой молекулы СбС1) в широком диапазоне энергий столкновения превышает сечения канала 4 примерно в 3 раза. Для системы 2 каналы 3 и 4 имеют примерно то же значение порога, но сечения канала 3 (в данном случае это диссоциация очень тяжелой молекулы С^) меньше сечений канала 4 (диссоциация ЯЬС1) в значительно большей степени, что соответствует представлениям об эффективности передачи энергии в зависимости от масс сталкивающихся партнеров. Положения максимумов функций возбуждения обоих каналов остались прежними (рис. 24).

Рис. 24. Функции возбуждения каналов 3 (а) и 4 (б) для взаимодействия в системах СэС1 + Ш(1)иС8Л-Ю>С1(2).

Каналы 5 и 6, в которых одновременно с процессами СИД образуется новая молекула, для системы 1 имеют колоколообразные функции возбуждения с максимумами соответственно при 11,0 эВ и 17,0 эВ. Сечение образования более тяжелой новой молекулы Сз1 в канале 5 в максимуме примерно в 1,5 раза превышает наибольшее значение сечения канала 6. Для взаимодействия в системе С^+ЯЬС1 максимумы функций возбуждения этих каналов поменялись местами и соответствуют 13,0 эВ и 10,0 эВ, а значения сечений в максимумах имеют близкие величины (см. рис. 25).

Рис. 25. Функции возбуждения каналов 5 (а) и б (б) для взаимодействия в системах СвС1 + ЛЫ (1) и СэЛ + ШзС1 (2).

Функции возбуждения обменной реакции в канале 7 качественно совпадают для обеих систем, но для взаимодействия в системе Сэ-Г + ЯЬС1 сечения процесса несколько выше:

Энергия столкновения, эВ

Рис. 26. Функции возбуждения канала 7 для взаимодействия в системах СяС1 + ЯЫ (1) и Сэ1 + ЯЬС1(2).

Функции возбуждения каналов 9 и 11, которые совпадают по продуктам для обеих систем, приведены на рис. 27. Если для системы 1 сечения этих каналов, соответствующих диссоциации одной из молекул реагентов с присоединением отрицательного иона к

оставшейся молекуле, были одного порядка величины, то для системы 2 они отличаются более чем в 10 раз. Причем сечение канала 9 для системы 2 превышает аналогичную величину для системы 1, а максимальное значение сечения канала II для системы 2 в 12 раз меньше сечения для системы 1.

Энергия столкновения, эВ Энергия столкновения, эВ

Рис. 27. Функции возбуждения каналов 9 (а) и 11 (б) для взаимодействия в системах СбС1 + ЯЫ (1) и С^ + ЯЬС1 (2).

С точки зрения образующихся продуктов каналы 8 и 10 для систем 1 и 2 поменялись местами, их функции возбуждения остались колоколообразными с максимумами при энергии столкновения около 4 эВ и для обеих систем по крайней мере на порядок выше сечение того канала, в котором связующим для двух положительных ионов является большой по размерам ион Г (соответственно канал 10 для системы 1 и канал 8 для системы 2). Функции возбуждения для каналов 8 и 10 показаны на рисунке 28.

Энергия столкновения, эВ Энергия столкновения, эВ

Рис. 28. Функции возбуждения каналов 8 (а) и 10 (б) для взаимодействия в системах СвС1 + ЯЫ (1) и С&1 + ЯЬС1 (2).

Сравнение динамики взаимодействия в симметричных системах СбС1 + ЯМ и С&1 + ЯЬС1 показывает, что поверхность потенциальной энергии является гладкой и обладает

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

Глава 5 содержит результаты траекторного моделирования динамики прямой трехтельной рекомбинации ионов Св+ и Вг" с участием атомов Кг, Хе и Щ в качестве стабилизирующей частицы [49,58,64,66].

Рекомбинация атомов, радикалов и ионов играет большую роль в кинетике разнообразных газофазных процессов и особенно в цепных и разветвленно-цепных реакциях, химии верхней атмосферы, горении и в низкотемпературной плазме. Редактор фундаментальной "Энциклопедии низкотемпературной плазмы" В.Е. Фортов в предисловии к первому тому отметил, что "физика низкотемпературной плазмы составляет научную основу создания масштабных исследовательских и технологических стендов, производственных линий и заводов". Свойства природных и искусственно создаваемых неравновесных сред в значительной степени определяются концентрацией заряженных частиц и их взаимодействием между собой, т.е. элементарными процессами появления заряженных пар и их последующего взаимодействия с нейтральными частицами и между собой, в том числе их нейтрализацией и рекомбинацией.

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

- прямая трехтельная рекомбинация, протекающая в одну стадию

А + В + С —АВ + С (3)

- механизм Линдемана с передачей энергии от возбужденного ассоциата третьему

телу

А + В 4-1- АВ* (4а)

АВ* + С —>АВ + С (4Ь)

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

В + С-*ВС* (5а)

ВС* + С - ВС + С (5Ь)

А + ВС —► АВ + С (5с)

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

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

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

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

Альтернативная перспектива исследования динамики рекомбинации вытекает из одного из первых принципов - принципа микроскопической обратимости. Законы классической и волновой механики допускают обращение времени, что собственно отражает принцип микроскопической обратимости. Инвариантность уравнений движения относительно обращения времени означает, что если при данном наборе начальных условий рассчитывать траекторию до того момента времени, когда сталкивающиеся молекулы разлетятся после реакции или без нее, то, обратив направление всех моментов и снова начав интегрирование, система вернется к первоначальным условиям. Это означает, что каждой траектории, ведущей из долины реагентов в долину продуктов, всегда можно сопоставить точную траекторию, ведущую в обратном направлении.

Поскольку прямая трехтельная рекомбинация и СИД двухатомных молекул являются точно обратными друг другу процессами, то динамика элементарного процесса рекомбинации может быть исследована на той же ППЭ, которая управляет процессом прямой диссоциации. На основе принципа микроскопической обратимости можно утверждать, что если СИД протекает как прямой процесс, то трехтельная рекомбинация также является прямой реакцией и ее сечение должно быть соответствующим сечению СИД. Последнее утверждение вытекает из статистического аналога принципа микроскопической обратимости - принципа детального равновесия. Рассматриваемый подход к исследованию динамики рекомбинации в настоящее время является по существу единственно возможным, поскольку полные квантовые расчеты таких процессов пока что недоступны. Он очень перспективен для исследований динамики рекомбинации ионов.

Динамика бимолекулярных процессов СИД нейтральных молекул с ионной связью с образованием ионов типа

М*Х" + 11-»М++Х" + 11, (6)

где М*Х" - молекула гапогенида металла с ионной связью и Я - снаряд (атом или молекула), подробно исследована в молекулярных пучках и практически количественно описывается квазиклассическими траекторными моделями на основе достаточно хорошо установленных поверхностей потенциальной энергии. На рисунке 29 приведено сравнение экспериментальных и полученных из траекторного счета функций возбуждения процессов СИД молекул СэВг в столкновениях с атомами Кг, Хе и Н(». Результаты расчетов совпадают количественно с экспериментальными данными.

Рис. 29. Сравнение экспериментальных и полученных из траекторного счета функций возбуждения процессов СИД: СбВг + Я Сб+ + Вг" + Я, где Я = Кг (а), Хе (б) и (в).

Использованные в этих расчетах ППЭ явились основой траекторного исследования прямой рекомбинации ионов Св+ и Вг' в присутствии атомов Кг, Хе и Щ в качестве третьего тела. На рисунке 30 показана использованная в расчетах система координат.

Рис. 30. Система координат для представления прямой трехтельной рекомбинации. Атом Я движется параллельно оси X, Ьтах - максимальный прицельный параметр, 0 и Ф -соответственно полярный и азимутальный углы, определяющие взаимное пространственное расположение линии сближения ионов и вектора скорости третьего атома.

Одной из основных характеристик динамического поведения системы взаимодействующих атомов является зависимость вероятности процесса от энергии столкновения, называемая функцией возбуждения. В отличие от бимолекулярных реакций на процесс трехтельной рекомбинации влияют две энергии - энергия сближения рекомбинирующих ионов Е1Ш и кинетическая энергия атома относительно центра масс ионной пары Еге!.. На рис. 31 показана гистограмма зависимости вероятности рекомбинации ионов с Хе в качестве третьего тела от этих двух энергий. Из рисунка следует, что процесс рекомбинации в диапазоне энергий сближения и относительных энергий от 0 до 10,0 эВ протекает без заметных особенностей. Вероятность рекомбинации

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

Сравнение соответствующих вероятностей рекомбинации как функций Еге| для различающихся по массе нейтральных атомов Кг, Хе и Н§ при Е™=1,0 эВ показано на рис. 32. Как видно из этого рисунка, вероятность рекомбинации для Кг всегда ниже, чем для двух других более тяжелых атомов. В области энергий Еге|. ниже 2,0 эВ атом Хе несколько более эффективен как третья частица по сравнению с атомом Для более высоких энергий столкновения эффективность тяжелого атома Ь^ выше.

Рис. 31. Гистограмма зависимости вероятности рекомбинации Се*" и Вг' в присутствии Хе от энергии сближения ионов Е|п|_ и энергии столкновения Еге]. ионной пары с атомом Хе.

Сб* + Вг' + Я -> СзВг" + Я

Относительная энергия атома И, эВ

Рис. 32. Зависимость вероятности рекомбинации Сб+ + Вг" + Я —» СэВг* + Я от кинетической энергии атома Я при Еы =1,0 эВ, где Я = Кг, Хе, Щ.

Функция прозрачности, показывающая зависимость вероятности рекомбинации от прицельного параметра нейтрального атома относительно центра масс ионной пары, имеет сложную бимодальную структуру, вид которой существенно зависит от обеих энергий столкновения, а также от свойств атома-акцептора энергии. На рис. 33 показаны функции прозрачности процесса рекомбинации с атомом Кг в качестве третьего тела для четырех комбинаций Е|Ш И Пге].

§_ 0,010

Се* + Вг" + Кг -> СаВг + Кг а

Е =1,0 эВ 101 '

Е_,=1,0эВ

0 4 8 12 16 20 Прицельный параметр, а е.

0,020 I 0,015

I

| 0,010 ^ 0,005

В.

ш 0,000

Сз* + Вг' + Кг -> СэВг + Кг б Е =1,0 эВ

им '

Е =5,0 эВ

0 4 8 12 16 20 Прицельный параметр, а е.

0,020

а

| 0,015 •8

10,010

I 0,005 |

И 0,000

Сэ+ + Вг" + Кг -> СзВг + Кг в ЕШ1=1,0эВ Е =7,0 эВ

4ШЦ

0 4 8 12 16 20 Прицельный параметр, а.е.

0,020

I

§ 0,015 к ю х

| 0,010

| 0,005 8.

Ш 0,000

Се* + Вг' + Кг -> СэВг + Кг Г Е^-1,0 эВ Е .=10,0 эВ

0 4 8 12 16 Прицельный параметр, а.е

20

Рис.-33. Функции прозрачности при рекомбинации в присутствии атомов Кг для 4 различных сочетаний Е,ш. и Еге1.

Обнаружено, что группы траекторий, соответствующие двум областям прицельных параметров, характеризуются существенно различными распределениями углов 0 и Ф, определяющих взаимное пространственное расположение линии сближения ионов и вектора скорости третьего атома. Для комбинации энергий Е,ш =1,0 эВ и ЕГС| =7,0 эВ эти распределения показаны на рисунке 34.

Рисунки 33 и 34 свидетельствуют о том, что процесс рекомбинации имеет, по крайней мере, два динамических механизма стабилизации образующейся молекулы, реализуемых в различных областях изменения прицельных параметров и при различных пространственных ориентациях векторов скорости партнеров.

0,0030 I 0,0025

I 0,0020

I 0,0015 § 0,0010

1 0,0005' &

® 0,0000

Cs* + Br + Kr CsBr + Kr Еы = 1,0эВ E ,=7,0 эВ

b>5 a.e.

60 90 120 150 180 Угол ©, град.

0,0020 I 0,0015

I

£ 0,0010 I 0,0005

I

10 0,0000

Cs* + Br + Kr -> CsBr + Kr 5 E =1,0 эВ

w

E ,-7,0 эВ

b>5 n.c. b<=5 a.e. ./

0 60 120 180 240 300 360 Угол Ф, град.

Рис. 34. Распределения ориентационных углов 0 (а) и Ф (б) при рекомбинации в присутствии атомов Кг при Е1Щ =1,0 эВ и Е^ -7,0 эВ для двух областей прицельных параметров.

Один из этих механизмов реализуется в первой области и состоит в стабилизации молекулы CsBr столкновением атома ксенона с обоими ионами при их сближении в центре масс пары. В первой области прицельных параметров большую роль играет не только взаимная ориентация векторов скорости сближения ионов и скорости стабилизирующего атома, но и фаза сближения ионов в момент столкновения с третьим телом. Если столкновение атома с ионной парой происходит в конфигурации, близкой к коллинеарной, на фазе сближения ионов, то вместо стабилизации ионной пары происходит передача ей дополнительной энергии от третьего тела. Вследствие этого, как видно из рисунка 34(а), вероятность рекомбинация крайне низкая при близких к нулю значениях угла 0. Для первой области прицельных параметров оптимальным, как следует из рисунка, является значение 0=60°. Очевидно, что при малых прицельных параметрах влияние угла Ф не столь значительно, поскольку конфигурация трехтельного столкновения мало меняется в зависимости от значения этого угла. Однако следует отметить все же несколько более высокую вероятность рекомбинации при значениях угла Ф вблизи 90° и 270°, т.е. при реализации прицельного параметра перпендикулярно линии сближения ионной пары, что увеличивает взаимодействие атома Хе с обоими ионами.

Другой механизм рекомбинации осуществляется в области больших прицельных параметров, где столкновение в области центра тяжести невозможно. В этом случае передача энергии от одного из ионов нейтральному атому может происходить на больших расстояниях от центра тяжести ионной пары при разлете последней после сближения и отталкивания. Этот механизм предполагает, что претерпевший отталкивание с цезием ион Вг" (или наоборот) сталкивается с атомом Я, который успевает оказаться на расстоянии и в конфигурации, благоприятной для передачи ему части поступательной энергии движения иона. В этой области рекомбинация наиболее вероятна при значениях угла Ф, соответствующих плоской ориентации векторов скорости частиц. И в этом случае взаимодействие Хе с Вг" значительно эффективней, чем с Cs+, что находит отражение в более высокой вероятности рекомбинации для углов Ф=0° и Ф=360°, чем для Ф=180°. Зависимость вероятности рекомбинации от угла 0 для области больших прицельных параметров похожа на аналогичную зависимость для первой области, однако ее максимум смещен в сторону больших значений угла 0 в силу того, что взаимодействие атома с ионом происходит на разлете ионной пары.

Одним из очень важных вопросов динамики прямой рекомбинации является распределение внутренних энергий продуктов рекомбинации и их зависимости от динамических параметров рекомбинации. Разделение траекторий на две группы по прицельным параметрам позволяет различить распределения колебательных и вращательных состояний молекул, образующихся в обеих областях прицельных параметров. На рис. 35 показаны эти распределения для обеих групп траекторий для Л=Кг с комбинацией энергий Б,,¡=1,0 эВ и Е^ =7,0 эВ. Из этого рисунка видно, что оба механизма рекомбинации приводят примерно к одинаковым распределениям по колебательной энергии. Существенное различие наблюдается только в области колебательных энергий СбВг в диапазоне 3,5-4,3 эВ, т.е. рекомбинация с большими прицельными параметрами несколько менее эффективна в плане отвода колебательной энергии ассоциата и перевода его в стабильное состояние. Распределения по вращательным энергиям практически одинаковы для обеих возможных областей прицельных параметров.

0 1 2 3 4 5 Колебательная энергия, эВ

7000

[? 6000

| 5000

5 4000 в-

зооо 1 2000 | 1000

Се* + Вг" + Кг -> СвВг + Кг Е1ш-1,0эВ

Ь>5 а.е Е.=7,0эВ

0,0 0,1 0.2 0,3 0,4 Вращательная энергия, эВ

0,5

Рис. 35. Распределения колебательных (а) и вращательных (б) состояний молекул СбВг при рекомбинации с участием атомов Кг при Е,„( =1,0 эВ и Еге| =7,0 эВ для двух диапазонов прицельных параметров.

ВЫВОДЫ

1. Разработаны комплексы программных средств для траекторного моделирования столкновительно-индуцированной диссоциации в трех- и четырехатомных системах, а также обратного процесса трехтельной рекомбинации ионов. Комплексы реализованы на языке Фортран и при использовании на персональных ЭВМ обеспечивают высокую производительность и точность расчетов.

2. Методом траекторного моделирования построена математическая модель для описания СИД с участием двух молекул с ионной связью в скрещенных молекулярных пучках. Сравнение полученных на основе предложенной модели данных с экспериментом показывает, что модель адекватно описывает экспериментальные результаты и может быть использована для детального исследования динамики столкновений.

3. Представленные результаты показывают полную применимость использованной диабетической ППЭ и ее параметров для описания различных каналов' СИД при взаимодействии двух молекул с ионной связью. Предложенная методология построения

ППЭ может быть использована для описания взаимодействия в системах с большим числом атомов.

4. Подробно исследовано влияние положения и формы отталкивательной стенки в потенциалах взаимодействия одноименно заряженных ионных пар на соотношения сечений различных каналов химического взаимодействия. Это позволяет рассмотреть основные динамические особенности появления и реализации каналов в многоканальном процессе с точки зрения топологических особенностей поверхности потенциальной энергии.

5. Полученные данные показывают, что динамика передачи энергии в системе СбС1 + ЯЫ является сложным процессом и не может быть адекватно представлена в рамках импульсной модели. Для коллинеарных конфигураций столкновения эффективность конверсии поступательной энергии во внутреннюю энергию молекул достаточно низка, и поэтому основной вклад в реализацию каналов химического рассеяния вносят взаимодействия с неколлинеарными конфигурациями столкновения. При низких энергиях столкновения межмолекулярный V—«V обмен играет заметную роль, в то время как V— обмен мало эффективен.

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

7. Сравнение результатов динамических расчетов в симметричных четырехатомных системах СзС1 + ЯЫ и СвД + ШзС1, управляемых одной и той же ППЭ, показывает, что поверхность потенциальной энергии является гладкой и обладает сходной топологией для движения изображающей точки из различных начальных положений во входной долине к области активного взаимодействия двух молекул и далее в долину продуктов каждого канала взаимодействия. Количественные отличия в динамическом поведении двух систем обусловлены разницей в соотношении масс молекул реагентов и, как следствие, различной эффективностью передачи энергии при заданных ориентациях столкновения.

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

9. Обнаружено существенное влияние взаимной ориентации вектора скорости стабилизирующего атома и линии сближения ионной пары на исход процесса рекомбинации. Для плоских конфигураций столкновения трех частиц при рекомбинации с участием атома Хе наиболее предпочтительным является коллинеарное столкновение атома с ионом брома, а рекомбинация с участием атома ртути с максимальной вероятностью происходит в случае, когда вектор его скорости перпендикулярен линии сближения ионной пары. В случае атома Кг рекомбинация происходит в том же диапазоне энергий сближения как для перпендикулярной конфигурации столкновения, так и для коллинеарной конфигурации со стороны иона Вг', в то время как коллинеарная конфигурация столкновения со стороны тяжелого иона цезия вообще не приводит к рекомбинации. Значительное влияние массы третьего тела на процесс рекомбинации обнаружено также по величинам сечений рекомбинации, которые оказываются тем выше, чем больше масса стабилизирующего атома.

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

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

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ В СЛЕДУЮЩИХ РАБОТАХ:

1. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Траекторное моделирование столкновительной диссоциации солей щелочных металлов в молекулярных пучках. // Деп. в ВИНИТИ 30.05.1989, № 3574-В89,56 с.

2. V.M.Akimov, V.M.Azriel, L.Yu.Rusin. Dynamics of the ionic complexes formation in collisional dissociation processes. // In: 12-th International Symposium on Molecular Beams, Italy, Perugia, 1989, p.41.

3. Азриель B.M., Акимов B.M., Грико Я., Русин Л.Ю. Динамика образования ионов Cs+ при столкновительной диссоциации CsBr с Хе в скрещенных молекулярных пучках. // Деп. в ВИНИТИ 12.09.1989, № 5811-В89,73 с.

4. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Динамика передачи энергии в столкновениях Хе с CsBr в диапазоне энергий от 0,35 до 11 эВ. И Деп. в ВИНИТИ 31.10.1989, №6583-В89, 83 с.

5. V.M.Akimov, V.M.Azriel, L.Yu.Rusin. Dynamics of CsBr dissociation by the ionic channel in collision with Xe. // In: 3-th International Symposium on Elementary Processes and Chemical Reactivity, Czechoslovakia, Prague, 1989, p.47.

6. Азриель B.M., Акимов B.M., Русин Л.Ю. Динамика образования комплексных ионов в столхновительно-индуцированной диссоциации. // Деп. в ВИНИТИ 3.08.1990, № 4478-В90,76 с.

7. Азриель В.М., Акимов В.М., Русин Л.Ю. Рассеяние ионов Cs+ при диссоциации CsBr в скрещенных молекулярных пучках. Эксперимент. // Хим. Физика, 1990, т.9, Ка 9, с.1224-1230.

8. Азриель В.М., Акимов В.М., Ленин Л.В., Русин Л.Ю. Рассеяние ионов Cs+ в процессах диссоциации молекул CsBr. II Отчет ИНЭП ХФ АН СССР, Москва, 1990, № 02900037121, 315 с.

9. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Рассеяние ионов Cs+ при диссоциации CsBr в скрещенных молекулярных пучках. Сравнение экспериментальных данных с траекторной моделью. // Хим. Физика, 1990, т.9, № 10, с. 1306-1310.

10. Азриель В М., Акимов В.М., Грико Я., Русин Л.Ю. Динамика реакции образования атомных ионов при диссоциации молекул CsBr в столкновениях с ксеноном. // Хим, Физика, 1990, т.9, № 11, с.1463-1470.

11. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Траеюорное моделирование элементарных процессов передачи поступательной энергии в столкновениях молекул CsBr с атомами Хе. // Хим. Физика, 1990, т.9, № 11, с.1471-1479.

12. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Orientation effects on energy transfer in atomdiatom collisions for the case of multiple impacts. // In: Proc. 8-th Symposium on Atomic and Surface Physics, Italy, Trentino, 1992, p.2.14.

13. Азриель B.M., Русин Л.Ю., Севрюк М.Б. Ориентационные эффекты в процессах передачи энергии при многократных соударениях в системе трех частиц: сравнение траекторного счета и вычислений в импульсном приближении. // Деп. в ВИНИТИ 3.04.1992, №1135.В92, 69 с.

14. V.M.Azriel, L.Yu.Rusin. Dynamics of four-atom interactions. // In: 9-th European Conference on Dynamics of Molecular Collisions (MOLEC - IX), Czechoslovakia, Prague,

1992, p.141.

15. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Impulsive approximation in energy transfer calculations. // In: 9-th European Conference on Dynamics of Molecular Collisions (MOLEC -IX), Czechoslovakia, Prague, 1992, p.50.

16. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. A description of energy transfer in the hard sphere mechanics: the case of multiple impacts. // In: Orientation and Polarization Effects in Chemical Reaction Dynamics (a NATO Advanced Workshop), Italy, Assisi, 1992, p.5.

17. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Multiple impacts and energy transfer in a three-body system for noncollinear collisions. // Theor. Chim. Acta, 1993, v.87, p.195-213.

18. Азриель B.M., Русин Л.Ю., Севрюк М.Б. Ориентационные эффекты в процессах передачи энергии при многократных соударениях в системе трех частиц. // Хим. Физика,

1993, т.12, № 11, с.1427-1443.

19. V.M.Azriel, L.Yu.Rusin. Dynamics of the interaction of two ionic diatomic molecules. // In: 15-th International Symposium on Molecular Beams, Germany, Berlin, 1993, D.21.

20. V.M.Akimov, V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Dissociation dynamics in the collisions of two diatomic molecules with ionic bonds. // In: 16-th International Symposium on Molecular Beams, Israel, Ma'ale Hachamisha, 1995, p.lA.

21. Азриель B.M., Русин Л.Ю., Севрюк М.Б. Импульсная модель взаимодействия пары двухатомных молекул с ионной связью. // Хим. Физика, 1995, т. 14, № 7, с.28-37.

22. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. A hard sphere model for the reaction of two diatomic molecules with ionic bond. // Chem. Phys., 1995, v. 199, p. 195-206.

23. Азриель B.M., Акимов B.M., Русин Л.Ю. Функции возбуждения основных каналов ионной диссоциации при столкновении двухатомных молекул с ионной связью. CsCl + RbJ.//Хим.Физика, 1996,т.15,№3, с.3-19.

24. V.M.Akimov, V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Excitation functions of cation formation in collisions of diatomic molecules with an ionic bond: CsCl + RbJ. // J. Chem. Soc. Faraday Trans., 1996, v.92, p. 1683-1688.

25. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk and J.P.Toennies. Construction of PESs with prescribed characteristics for triatomic systems by an extended LEPS method. II In: 17-th International Symposium on Molecular Beams, France, Orsay, 1997, P3A.

26. V.M.Akimov, V.M.Azriel and L.Yu.Rusin. Product angular distributions in the four-atom system CsCl + RbJ for the collision energy ranging from 4 to 6 eV. // In: 17-th International Symposium on Molecular Beams, France, Orsay, 1997, P15B.

27. V.M.Azriel and L.Yu.Rusin. Channel branching in the collisions between CsCl and RbJ. // In: All Russian Conference on Theoretical Chemistry (ARCTC97), Russia, Kazan, 1997, P8.

28. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk and J.P.Toennies. Effect of the angular dependence of the barrier height on the features of the F + H2 reaction. // Chem. Phys., 1998, v.232, p.307-320.

29. Азриель B.M., Русин Л.Ю., Севрюк М.Б., Тоннес Я.П. Связь особенностей топографии поверхности потенциальной энергии с динамикой реакции атомов фтора с молекулярным водородом.//Журн. Физ. Хим., 1998,т.72,№ 8,с.1351-1359.

30. V.M.Akimov, V.M.Azriel and L.Yu.Rusin. Dynamics of collision-induced dissociation in four-atom systems. // In: 12-th European Conference on Dynamics of Molecular Collisions (MOLEC ХП), Bristol, UK, 1998, p.28.

31. V.M.Azriel and L.Yu.Rusin. Influence of parameters in repulsive potentials of the ions with the same sign on the dynamics of the interaction of two molecules with ionic bond. // In: 10-th Conference on Current Trends in Computational Chemistry (CCTCC), Jackson, Mississippi, USA, 2001, p 33-35 (PI).

32. Азриель B.M., Акимов B.M., Русин Л.Ю. Траекторное моделирование динамики столкновительной диссоциации двухатомных молекул в 4-х атомных системах. // Хим. Физика, 2002, т.21, № 4, с.18-27.

33. V.M.Azriel and L.Yu.Rusin. Quasiclassical trajectory simulation of the dynamics of intermolecular interactions in the system CsCl + RbJ for the collision energy ranging from 0,1 to 25 eV. // In: 6-th World Congress of Theoretically Oriented Chemists (WATOC'02), Lugano, Switzerland, 2002, PB277.

34. V.M.Azriel and L.Yu.Rusin. Influence of the topography of potential energy surface on the channels formation in the collisions of two molecules with ionic bond. // In: 6-th World Congress of Theoretically Oriented Chemists (WATOC'02), Lugano, Switzerland, 2002, PB280.

35. V.M.Azriel and L.Yu.Rusin. Branching in the ion pairs formation in the collisions of CsCl and RbJ for the energy ranging from 0,1 to 25 eV. // In: 20-th International Symposium on Molecular Beams, Portugal, Lisbon, 2003, p.139-141 (P9A).

36. V.M.Azriel and L.Yu.Rusin. Exchange reaction in the interaction of CsCl and RbJ for the energy ranging from 0,1 to 17 eV. II In: 20-th International Symposium on Molecular Beams, Portugal, Lisbon, 2003, p.142-144 (P10A),

37. V.MAzriel and L.Yu.Rusin. Dynamics of energy transfer for the interaction of CsCl and RbJ in cross molecular beams. // In: 18-th International Conference on Molecular Energy Transfer (COMET_XVIII), Spain, San Lorenzo de El Escorial, 2003, p.40 (10-TH).

38. V.M.Azriel and L.Yu.Rusin. Trajectory simulation of the recombination reaction Cs+ + Br" + Xe -> CsBr + Xe. // In: 14-th Symposium on Atomic, Cluster and Surface Physics (SASP 2004), Italy, La Thuile, 2004, PC-4.

39. V.MAzriel and L.Yu.Rusin. Dynamics of four-atom complex formation in the collisions of CsCl and RbJ in cross molecular beams. // In: 8-th European Conference on Atomic and Molecular Physics (ECAMP VIII), France, Rennes, 2004,3-7 (A-3-7).

40. V.M.Azriel and L.Yu.Rusin. Quasiclassical trajectory simulation of recombination process Cs+ + Br' + Hg —► CsBr + Hg. II In: 8-th European Conference on Atomic and Molecular Physics (ECAMP VIII), France, Rennes, 2004,11-3 (B-7-3).

41. V.M.Azriel, D.B.Kabanov, L.Yu.Rusin. Dynamics of collision-induced dissociation of two molecules with ionic bond. // In: 24-th International Symposium on Rarefied Gas Dynamics (RGD24), Ban, Italy, 2004, AIP Conference Proceedings 762, p.926.

42. V.M.Azriel, L.Yu.Rusin. Three-atom complexes formation in the collisions of two molecules with ionic bond. // In: International Conference on Dynamics of Molecular Systems (MOLEC 2004), Nunspeet, Netherlands, 2004, p.65 (РОЗ).

43. V.M.Azriel, L.Yu.Rusin, M.B.Seviyuk. The complete T—>V,R energy conversion in hard sphere A+BC collisions. // In: International Conference on Dynamics of Molecular Systems (MOLEC 2004), Nunspeet, Netherlands, 2004, p.66 (P04).

44. V.M.Azriel, L.Yu.Rusin. Trajectory simulation of the influence of internal energy of reagents on the interaction in the system CsCl + RbJ. // In: International Conference on Dynamics of Molecular Systems (MOLEC 2004), Nunspeet, Netherlands, 2004, p.64 (P02).

45. V.M.Azriel, L.Yu.Rusin, M.B.Seviyuk. The complete T—>V,R energy conversion in three-body collisions within the hard sphere model. II J. Chem. Phys. 122,074322 (2005).

46. V.M.Azriel, L.Yu.Rusin. Influence of internal energy of reagents on the dynamics of exchange reaction CsCl + RbJ -» CsJ + RbCl. // In: 7-th Congress of the World Association of Theoretically Oriented Chemists (WATOC'05), Cape Town, South Africa, 2005, p. 183 (RM-P13).

47. V.M.Azriel, L.Yu.Rusin. Conversion of translation energy in the collisions of two diatomic molecules. // In: 7-th Congress of the World Association of Theoretically Oriented Chemists (WATOC'05), Cape Town, South Africa, 2005, p. 184 (RM-P14).

48. V.M.Azriel, L.Yu.Rusin. Influence of energy of approach of ionic pair on internal excitation of a formed molecule in reaction of recombination Cs+ + Br' + Xe -> CsBr* + Xe. II In: 21-th International Symposium on Molecular Beams, Hersonissos, Crete, Greece, 2005, p.52 (P09).

49. V.M.Azriel, L.Yu.Rusin. Dynamic features of recombination of ionic pair Cs+ - Br at the presence of atoms Xe and Hg. // In: 21-th International Symposium on Molecular Beams, Hersonissos, Crete, Greece, 2005, p.53 (PI 1).

50. V.M.Azriel, L.Yu.Rusin. Influence of internal energy of molecules CsCl and RbJ on the formation of three-atomic ionic complexes. II In: 21-th International Symposium on Molecular Beams, Hersonissos, Crete, Greece, 2005, p.54 (P13).

51. V.M.Azriel, L.Yu.Rusin. Influence of weight of third body on probability of recombination in the system Cs+ + Br" + R —> CsBr* + R. II In: 6-th International Conference on Chemical Kinetics, Gaithersburg, Maryland, USA, 2005, F20.

52. V.M.Azriel, L.Yu.Rusin. The effect of internal energy of molecules CsCl and RbJ on dynamics of their collision-induced dissociation. // In: б-th International Conference on Chemical Kinetics, Gaithersburg, Maryland, USA, 2005, C20

53. V.M.Azriel, L.Yu.Rusin. Trajectory simulation of the effect of internal excitation of reagents on dynamics of long-living four-atom complex formation in the system CsCl + RbJ. // In: 6-th International Conference on Chemical Kinetics, Gaithersburg, Maryland, USA, 2005, L21.

54. V.M.Azriel, L.Yu.Rusin. Dynamics of ionic recombination in the system Cs+ + Br' + > CsBr* + R, where R = Hg, Xe and Kr. // In: 15-th Symposium on Atomic, Cluster and Surface Physics (SASP 2006), Austria, Obergurgl, 2006, p.156-159.

55. V.M.Azriel, L.Yu.Rusin. Trajectory simulation of the interaction of CsCl and RbJ in cross molecular beams at high collision energies. // In: 15-th Symposium on Atomic, Cluster and Surface Physics (SASP 2006), Austria, Obergurgl, 2006, p.153-155.

56. V.M.Azriel, L.Yu.Rusin. Collision-induced dissociation in the systems CsCl + RbJ and CsJ + RbCl. II In: 19-th International Symposium on Gas Kinetics (GK 2006), Orleans, France, 2006, p.383.

57. V.M.Akimov, V.M.Azriel, D.B.Kabanov, L.Yu.Rusin. Dynamics of collision-induced dissociation in the symmetric systems CsCl + RbJ and CsJ + RbCl. II In: European Conference

on Dynamics of Molecular Systems (MOLEC XVI 2006), Levico Termo (TN), Italy, 2006, p. 127-129.

58. Азриель B.M., Русин Л.Ю. Квазиклассическое траекторное моделирование трехтельной рекомбинации ионов Cs+ и Вг" в присутствии атомов Кг, Хе и Hg. // Физико-химическая кинетика в газовой динамике, 2006, т.4., http://www.chemphys.edu.ru/pdf/2006-10-23-002.pdf.

59. V.M.Azriel, L.Yu.Rusin. Dynamics of direct three-body recombination of the heavy ions. // In: 22-th International Symposium on Molecular Beams, Freiburg, Germany, 2007, p.122-125

60. V.M.Azriel, L.Yu.Rusin. The effect of hardness parameters in the pairwise potentials СГ -J' and Cs+ - Rb+ on the excitation function of the exchange reaction CsCl + RbJ —»CsJ + RbCl. //In: 22-th International Symposium on Molecular Beams, Freiburg, Germany, 2007, p.126-127.

61. V.M.Azriel, L.Yu.Rusin. The effect of vibrational and rotational excitation of the reagents on dynamics of the interaction in the system CsCl + RbJ. // In: 9-th European Conference on Atoms, Molecules & Photons, (ECAMPIX), Crete, Greece, 2007, Mo5-2.

62. V.M.Azriel, L.Yu.Rusin. Internal energy of three-atomic ionic complexes for the interaction of CsCl and RbJ in crossed molecular beams. // In: 20-th International Conference on Molecular Energy Transfer, (COMET XX), Arcachon, France, 2007, P-33.

63. V.M.Azriel, L.Yu.Rusin. Trajectory simulation of T—>V energy transfer for the collinear collisions in the system CsCl + RbJ. // In: 20-th International Conference on Molecular Energy Transfer, (COMET XX), Arcachon, France, 2007, P-l.

64. Азриель B.M., Кабанов Д.Б., Колесникова Л.И., Русин Л.Ю. Динамика рекомбинации ионов в низкотемпературной плазме. // Известия Академии Наук, серия Энергетика, 2007, №5, с. 50-69.

65. V.M. Azriel, L.Yu. Rusin. Dynamics of intra- and intermolecular energy exchange in the collisions of two molecules with ionic bond. II In: 16-th Symposium on Atomic, Cluster and Surface Physics (SASP 2008), Les Diablerets, Switzerland, 2008, p.95.

66. Азриель B.M., Русин Л.Ю. Динамика прямой трехтельной рекомбинации ионов. // Хим. Физика, 2008, т.27, № 7, с.5-17.

67. Азриель В.М., Акимов В.М., Русин Л.Ю. Траекторное моделирование динамики образования ионных пар при столкновительной диссоциации двухатомных молекул. II Отчет ИНЭПХФ РАН, Москва, 2008, № 02200803601,40 с.

68. Азриель В.М., Акимов В.М., Кабанов Д.Б., Русин Л.Ю. Детальная динамика элементарных процессов образования и рекомбинации ионов в области энергий столкновения от 0 до 10 эВ. // Отчет ИНЭП ХФ РАН, Москва, 2008, № 02200803603,78 с.

69. V.M. Azriel, L.Yu. Rusin. T-»V, R energy exchange for inelastic collisions of cold molecules in the system CsBr + CsBr. // In: 4-th Conference on Elemental Processes in Atomic Systems (CEPAS 2008), Cluj-Napoca, Romania, 2008, p.59 (We-1).

70. V.M. Azriel, L.Yu. Rusin. The effect of internal state of diatomic molecules on the dynamics of energy exchange. // In: 4-th Conference on Elementary Processes in Atomic Systems (CEPAS 2008), Cluj-Napoca, Romania, 2008, p.\32 (Th-14).

Отпечатано в ООО «Компания Спугник+» ПД № 1-00007 от 25.09.2000 г. Подписано в печать 29.09.08. Тираж 100 экз. Усл. п. л 2,75 Печать авторефератов: 730-47-74,778-45-60

ВВЕДЕНИЕ

Глава 1. СТОЛКНОВИТЕЛЬНО-ИНДУЦИРОВАННАЯ ДИССОЦИАЦИЯ В 12 СИСТЕМЕ СэО + ЫМ. СРАВНЕНИЕ С ЭКСПЕРИМЕНТОМ И УТОЧНЕНИЕ ПАРАМЕТРОВ ППЭ

1.1. Метод траекторного моделирования и используемая ППЭ

1.2. Сравнение результатов траекторного моделирования с экспериментом 19 и уточнение параметров ППЭ

1.3. Исследование влияния параметров жесткости в потенциалах 31 взаимодействия одноименно заряженных ионов на функции возбуждения различных каналов взаимодействия

Глава 2. ДИНАМИКА КОНВЕРСИИ ЭНЕРГИИ СТОЛКНОВЕНИЯ ВО 41 ВНУТРЕННЮЮ ЭНЕРГИЮ МОЛЕКУЛ С ИОННОЙ СВЯЗЬЮ

2.1. Коллинеарные столкновения. Сравнение с импульсной моделью

2.2. Распределения поступательной, вращательной и колебательной 54 энергий в неупругих столкновениях молекул

2.3. Динамический анализ траекторий

Глава 3. ДИНАМИКА ПОЯВЛЕНИЯ И РЕАЛИЗАЦИИ КАНАЛОВ 76 ХИМИЧЕСКОГО РАССЕЯНИЯ В СИСТЕМЕ СэС! + ИМ

3.1. Каналы столкновительно-индуцированной диссоциации (каналы 2-4)

3.1.1. Канал 2: СбС1 + ЯЫ Сэ+ + СГ + ЯЬ+ + Г

3.1.2. Канал 3: СэС1 + ЯЫ -» Сб+ + СГ + ЯЫ

3.1.3. Канал 4: СэС1 + ЯЫ ЯЬ+ + У + СзС

3.2. Каналы 5-7, СИД с одновременным образованием новых продуктов

3.2.1. Канал 5: СбО + ЯЫ -» СэД + ЯЬ+ + СГ

3.2.2. Канал 6: СэО + ЯЫ -> ЯЬС1 + Сз+ + Г

3.2.3. Канал 7: СзС1 + ЯЫ -> Сэ1 + ЯЬС

3.3. Каналы образования трехатомных ионных комплексов (каналы 8-11)

3.3.1. Канал 8: СэС1 + ЯЫ -> СбС1ЯЬ+ +1"

3.3.2. Канал 9: СбС! + ЯЫ -> СвСИ" + ЯЬ+

3.3.3. Канал 10: СэС1 + ЯЫ -> СзЖЬ+ + СГ

3.3.4. Канал 11: СбС1 + ЯЫ -> ЯЫСГ + Сэ+

Глава 4. РОЛЬ ВНУТРЕННЕЙ ЭНЕРГИИ РЕАГЕНТОВ. ВЛИЯНИЕ 168 МАССОВОГО ФАКТОРА: ДИНАМИКА ВЗАИМОДЕЙСТВИЯ В СИММЕТРИЧНЫХ СИСТЕМАХ СэС1 + И + Ш>С

4.1. Влияние внутренней энергии реагентов на динамику взаимодействия 168 в системе СэСИНЫ

4.1.1. Канал

4.1.2. Каналы 3,

4.1.3. Каналы 5,

4.1.4. Канал

4.1.5. Каналы 8

4.1.6. Канал

4.2. Динамика столкновительной диссоциации в симметричных 183 четырехатомных системах СэС1 + ЯЫ и Сн.Г + ЯЬС

Глава 5. ДИНАМИКА ПРЯМОЙ ТРЕХТЕЛЬНОЙ РЕКОМБИНАЦИИ ИОНОВ 192 Сэ+ И Вг С УЧАСТИЕМ АТОМОВ Кг, Хе И Щ

5.1. Математическая постановка задачи и технология расчетов

5.2. Результаты и обсуждение

5.2.1. Рекомбинация в системе Сэ+ + Вг" + Хе

5.2.2. Рекомбинация в системе Сэ+ + Вг" + Щ

5.2.3. Рекомбинация в системе Сб+ + Вг" + Кг 249 ЗАКЛЮЧЕНИЕ. ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ 266 СПИСОК ЛИТЕРАТУРЫ

Актуальность проблемы

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

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

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

В помощь исследователю приходит компьютерное моделирование, адекватность которого реальному процессу проверяется сравнением результатов расчета с теми характеристиками процесса, которые непосредственно могут быть измерены в эксперименте. Конечной целью такого моделирования является восстановление поверхности потенциальной энергии (ППЭ), необходимой для решения квантово-механической задачи рассеяния в приближении Борна-Оппенгеймера о разделении квантовой системы на быструю и медленную подсистемы. Однако во многих случаях физически оправданным и точным является классический или квазиклассический подходы к математическому моделированию динамического поведения системы без решения чрезвычайно трудоемких квантово-механических задач. Так как длина волны де Бройля (А= к/ту) мала по сравнению с расстояниями, на которые перемещаются атомы в самом акте реакции, то классические и квазиклассические методы работают достаточно точно, особенно для больших масс и энергий. И хотя поведение атомов и молекул при химических превращениях подчиняются строгим законам квантовой теории, подавляющее большинство таких процессов качественно правильно описывается уравнениями классической механики с квантовыми поправками, касающимися, как правило, внутренних состояний реагентов и продуктов взаимодействия.

Наиболее подробные сведения об элементарном процессе можно извлечь из траекторного моделирования этого процесса, разумеется, при наличии ППЭ, адекватно описывающей элементарный акт химического взаимодействия. Траекторное моделирование динамики элементарного процесса занимает промежуточное место между неэмпирическими расчетами и экспериментом и опирается как на теоретические представления о природе взаимодействия частиц, так и на эмпирические значения характеристик взаимодействия. Этот вид моделирования в настоящее время представляет собой наиболее распространенную, обладающую широкими пределами применимости процедуру выяснения связи непосредственно измеряемых характеристик рассеяния с динамическими особенностями каналов взаимодействия. Соответствующая связь устанавливается тем, что классические гамильтоновы уравнения движения ядер в заданном канале решаются на поверхностях потенциальной энергии взаимодействия в электронных состояниях, через которые система проходит между исходными (реагенты) и конечными (продукты) состояниями. Хорошее совпадение результатов динамических расчетов и экспериментальных данных свидетельствует о правильности используемой ППЭ и может служить гарантией достоверности расчетных характеристик процесса, не поддающихся прямому экспериментальному измерению.

Как отмечено выше, одной из важнейших задач, определяющих корректность результатов траекторного моделирования элементарных процессов, является построение ППЭ, адекватной физической картине исследуемого явления. Однако, определение поверхности потенциальной энергии, адекватной рассматриваемому процессу, представляет одну из наиболее трудных задач квантовой химии, если речь идет об "ab initio" расчетах ППЭ. В настоящее время, несмотря на огромные достижения в развитии вычислительной техники, известны лишь несколько таких ППЭ, построенных для наиболее простых систем, состоящих не более чем из трех атомов. Кроме того, полученные неэмпирическими методами потенциалы обычно представляют собой таблицы, что создает определенные неудобства при их использовании в вычислительных процедурах. В большинстве случаев эти табличные значения аппроксимируются гладкими математическими функциями, удобными в применении, но не позволяющими впоследствии проанализировать, какие именно характеристики поверхности потенциальной энергии ответственны за те или иные динамические особенности поведения исследуемой системы.

Поэтому большое распространение получили более доступные полуэмпирические и эмпирические методы построения ППЭ. Эти методы имеют два существенных преимущества в применении к динамическим исследованиям. Первое, они используют доступные экспериментальные данные для калибровки ППЭ, что позволяет воспроизвести топографические особенности поверхности, оказывающие наиболее существенное влияние на динамику процесса. И второе, они, как правило, допускают аналитическое представление ППЭ в виде функций, значительно более удобных с вычислительной точки зрения, что позволяет реализовывать более эффективные и быстрые алгоритмы и программы динамических расчетов. Поверхности потенциальной энергии, построенные с помощью таких методов, включают ряд подбираемых параметров, изменение которых позволяет существенным образом менять топологию Г1ПЭ и, таким образом, исследовать влияние тех или иных топологических особенностей ППЭ на динамику элементарного взаимодействия.

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

Цель работы

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

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

3. Определить происхождение основных каналов, возникающих при столкновительной диссоциации в четырехатомной системе, и их взаимосвязь.

4. Разработать методологию и провести, основываясь на принципе микроскопической обратимости, исследование динамики прямой трехтельной рекомбинации ионов на поверхностях потенциальной энергии, количественно описывающих обратные процессы СИД в соответствующих системах.

5. Построить траекторные модели рассматриваемых процессов СИД и рекомбинации, позволяющие понять механизмы химических превращений и их зависимость от динамических и кинематических особенностей взаимодействия.

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

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

Научная новизна работы

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

2. Исследовано влияние параметров, определяющих отталкивание в потенциалах взаимодействия одноименных ионов, на топологию ППЭ и на динамику многоканального взаимодействия, условия реализации и сечения отдельных элементарных процессов.

3. Построены траекторные модели процессов СИД в четырехатомной системе, учитывающие как динамические (структура и топология ППЭ), так и кинематические (энергетические параметры, ориентационные и массовые эффекты) характеристики взаимодействия.

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

5. Создан высокопроизводительный и удобный в использовании пакет вычислительных средств для компьютерного моделирования процессов столкновительно-индуцированной диссоциации и рекомбинации.

На защиту выносятся следующие положения

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

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

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

Практическая значимость работы

1. Результаты и выводы, полученные в работе, позволяют сделать заключение о высокой эффективности траекторных методов для исследования динамики элементарных процессов. Разработаны пакеты вычислительных средств для траекторного моделирования процессов СИД и трехтельной рекомбинации, обеспечивающие высокую производительность и точность расчетов при реализации на персональных компьютерах.

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

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

4. Использованная методология построения ППЭ для трех- и четырехатомных систем с успехом может быть применена для конструирования поверхностей потенциальной энергии для систем с большим числом атомов.

5. Построенные модели СИД и рекомбинации могут быть использованы для понимания основных закономерностей процессов, протекающих в верхних слоях атмосферы, низкотемпературной плазме и других практически важных областях естествознания.

Апробация результатов работы

Результаты, вошедшие в диссертацию, были представлены на следующих конференциях и симпозиумах: XII Международный симпозиум по молекулярным пучкам (Италия, 1989), III Международный симпозиум по элементарным процессам и химическим реакциям (Чехословакия, 1989), VIII Симпозиум по атомной физике и физике поверхности (Италия, 1992), IX Европейская конференция по динамике молекулярных столкновений (Чехословакия, 1992), Конференция НАТО по ориентационным и поляризационным эффектам в динамике химических реакций (Италия, 1992), XV Международный симпозиум по молекулярным пучкам (Германия, 1993), XVI Международный симпозиум по молекулярным пучкам (Израиль, 1995), XI Европейская конференция по динамике молекулярных столкновений (Дания, 1996), Международный симпозиум по применению компьютеров в химических исследованиях (Москва, 1996), XVII Международный симпозиум по молекулярным пучкам (Франция, 1997), Общероссийская конференция по теоретической химии (Казань, 1997), XII Европейская конференция по динамике молекулярных столкновений (Великобритания, 1998), X Конференция по современным тенденциям в вычислительной химии (США, 2001), VI Всемирный конгресс теоретически ориентированных химиков (Швейцария, 2002), XX Международный симпозиум по молекулярным пучкам (Португалия, 2003), XVIII Международная конференция по передаче молекулярной энергии (Испания, 2003), XIV Симпозиум по атомной физике, физике кластеров и поверхности (Италия, 2004), VIII Европейская конференция по атомной и молекулярной физике (Франция, 2004), XXIV Международный симпозиум по динамике разреженного газа (Италия, 2004), Международная конференция по динамике молекулярных систем (Нидерланды, 2004), VII Конгресс всемирной ассоциации теоретически ориентированных химиков (ЮАР, 2005), XXI Международный симпозиум по молекулярным пучкам (Греция, 2005), VI Международная конференция по химической кинетике (США, 2005), XV Симпозиум по атомной физике, физике кластеров и поверхности (Австрия, 2006), XIX Международный симпозиум по газовой кинетике (Франция, 2006), Европейская конференция по динамике молекулярных систем (Италия, 2006), IX Европейская конференция по атомам, молекулам и фотонам (Греция, 2007), XX Международная конференция по передаче молекулярной энергии (Франция, 2007), XXII Международный симпозиум по молекулярным пучкам (Германия, 2007), XXV Международная конференция по фотонным, электронным и атомным столкновениям (Германия, 2007), Симпозиум по диссоциативной рекомбинации (Нидерланды, 2007), XVI Конференция по современным тенденциям в вычислительной химии (США, 2007), XVI Симпозиум по атомной физике, физике кластеров и поверхности (Швейцария, 2008).

Личное участие автора

Автор принимал личное определяющее участие в постановке задач и разработке подходов для их решения, создании пакетов вычислительных программ и проведении расчетов, а также анализе полученных результатов.

Публикации

Основное содержание диссертации опубликовано в 70 работах, список которых приведен в конце автореферата.

Объем и структура диссертации

Диссертация состоит из введения, пяти глав, заключения и списка литературы. Общий объем - 299 страниц, в том числе 211 рисунков и 5 таблиц. Список литературы содержит 417 наименований.

ЗАКЛЮЧЕНИЕ. ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ

1. Разработаны комплексы программных средств для траекторного моделирования столкновительно-индуцированной диссоциации в трех- и четырехатомных системах, а также обратного процесса трехтельной рекомбинации. Комплексы реализованы на языке Фортран и при использовании на персональных ЭВМ обеспечивают высокую производительность и точность расчетов.

2. Методом траекторного моделирования построена математическая модель для описания СИД с участием двух молекул с ионной связью в скрещенных молекулярных пучках. Сравнение полученных на основе предложенной модели данных с экспериментом показывает, что модель адекватно описывает экспериментальные результаты и может быть использована для детального исследования динамики столкновений.

3. Представленные результаты показывают полную применимость использованной диабатической ППЭ и ее параметров для описания различных каналов СИД при взаимодействии двух молекул с ионной связью. Предложенная методология построения ППЭ может быть использована для описания взаимодействия в системах с большим числом атомов.

4. Подробно исследовано влияние положения и формы отталкивательной стенки в потенциалах взаимодействия одноименно заряженных ионных пар на соотношения сечений различных каналов химического взаимодействия. Это позволяет рассмотреть основные динамические особенности появления и реализации каналов в многоканальном процессе с точки зрения топологических особенностей поверхности потенциальной энергии.

5. Полученные данные показывают, что динамика передачи энергии в системе СэС1 + ЯЫ является сложным процессом и не может быть адекватно представлена в рамках импульсной модели. Для коллинеарных конфигураций столкновения эффективность конверсии поступательной энергии во внутреннюю энергию молекул достаточно низка, и поэтому основной вклад в реализацию каналов химического рассеяния вносят взаимодействия с неколлинеарными конфигурациями столкновения. При низких энергиях столкновения межмолекулярный V—>У обмен играет заметную роль, в то время как V—>-11 обмен мало эффективен.

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

7. Сравнение результатов динамических расчетов в симметричных четырехатомных системах СэС1 + ИЫ и Сз.Г + КЬСЛ. управляемых одной и той же ППЭ, показывает, что поверхность потенциальной энергии является гладкой и обладает сходной топологией для движения изображающей точки из различных начальных положений во входной долине к области активного взаимодействия двух молекул и далее в долину продуктов каждого канала взаимодействия. Количественные отличия в динамическом поведении двух систем обусловлены разницей в соотношении масс молекул реагентов и, как следствие, различной эффективностью передачи энергии при заданных ориентациях столкновения.

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

9. Обнаружено существенное влияние взаимной ориентации вектора скорости стабилизирующего атома и линии сближения ионной пары на исход процесса рекомбинации. Для плоских конфигураций столкновения трех частиц при рекомбинации с участием атома Хе наиболее предпочтительным является коллинеарное столкновение атома с ионом брома, а рекомбинация с участием атома ртути с максимальной вероятностью происходит в случае, когда вектор его скорости перпендикулярен линии сближения ионной пары. В случае атома Кг рекомбинация происходит в том же диапазоне энергий сближения как для перпендикулярной конфигурации столкновения, так и для коллинеарной конфигурации со стороны иона Вг", в то время как коллинеарная конфигурация столкновения со стороны тяжелого иона цезия вообще не приводит к рекомбинации. Значительное влияние массы третьего тела на процесс рекомбинации обнаружено также по величинам сечений рекомбинации, которые оказываются тем выше, чем больше масса стабилизирующего атома.

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

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

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

1. Kumar M., Shanker J. Electronic polarizabilities, potential fonctions, and spectroscopic constants for diatomic molecules of alkali halides and alkali hydrides. J. Chem. Phys., 1992, v.96, № 7, p.5289.

2. Каплан И.Г. Введение в теорию межмолекулярных взаимодействий. М., "Наука", 1982.

3. Chauhan R.S., Sharma S.С., Sharma S.В., Sharma В.S. Analysis of polarizabilities, potentials, and geometries of alkali halide dimers. J. Chem. Phys., 1991, v.95, № 6, p.4397.

4. Varandas A.J.C. Intermolecular and intramolecular potentials. Advances in Chemical Physics. Eds. I.Prigogine and S.A.Rice, 1988, v.74, p.255.

5. Hirst D.M. Potential energy surfaces: Molecular structure and reaction dynamics. London, 1985.

6. Maltz C. Palarizable-ion model for vibrationally excited alkali halides. Chem. Phys. Lett., 1969, v.3, p.707.

7. Степанов Н.Ф. Поверхности потенциальной энергии. Современное естествознание: Энциклопедия. М.: Магистр-пресс, 2000, т.1, с.21.

8. Chang-Hwan Rhee. Dissociation energies and dipole moments of alkali halides. Journal of the Korean Chemical Society, 1997, v.41, № 9, p.449.

9. Ohwada K. Application of potential constants: molecular chemical potentials of heteronuclear diatomic molecules. J. Phys. Chem., 1992, v.96, p.5825.

10. Никитин E.E. Теория элементарных атомно-молекулярных процессов в газах. М.:Химия, 1970.

11. Alexander A.J., Zare R.N. Anatomy of elementary chemical reactions. J. Chem. Edu., 1998, v.75, № 9, p.l 105.

12. Miller W.H. The semiclassical nature of atomic and molecular collisions. Acc. Chem. Res., 1971, v.4, p.161.

13. Genry R.W. Atom-molecule collision theory. A guide for experimentalist. Plenum Press, NY, 1979, p.391.

14. Hirschfelder J.O., Curtiss C.F., Bird R.B. Molecular theory of gases and liquids. Wiley, NY, 1964.

15. Cashion J.K. Testing of diatomic potential-energy functions by numerical methods. J. Chem. Phys., 1963, v.39, p. 1872.

16. Present R.D. Kinetic theory of gases. McGraw-Hill, NY, 1958.

17. Diestler D.J. Collision-induced dissociation. I. Quantal treatment. In: Atom-molecule collision theory. Ed. R.B.Bernstein. Plenum Press, 1979, p.655.

18. Pauling L., Wilson E.B. Introduction to quantum mechanics. McGraw-Hill, NY, 1935.

19. Empedocles P.B. Force constants of diatomic molecules. I. A general semitheoretical approach. J. Chem. Phys., 1967, v.46, p.4474.

20. Duff J.W., Truhlar D.G. Effect of curvature of the reaction path on dynamic effects in endothermic chemical reactions and product energies in exothermic reactions. J. Chem. Phys., 1975, v.62, p.2477.

21. Bernstein R.B., Levine R.D. Role of energy in relative molecular scattering: an information theoretic approach. Adv. At. Mol. Phys., 1975, v.l 1, p.215.

22. Herschbach D.R. Molecular dynamics of chemical reactions. Pure Appl. Chem., 1976, v.47, p.61.

23. Mishra M., Guzman R., Rabitz H. Role of potential structure in nonadiabatic collisions. Phys. Rev. A, 1987, v.36, p.l 124.

24. Goodisman J. Diatomic interaction potential theory. Academic Press, NY, 1973.

25. Waldman M., Gordon R.G. Scaled electron gas approximation for intennolecular forces. J. Chem. Phys., 1979, v.71, p.1325.

26. Agrawal G.G., Shanker J. Theory of short range interactions and static polarizabilities of alkali halides. J. Chem. Phys., 1978, v.68, № 11, p.5244.

27. Полуэмпирические методы расчета электронной структуры. Ред. Дж. Сигал, т.1, 2. М.: Мир, 1980.

28. Никитин Е.Е., Уманский С.Я. Неадиабатические переходы при медленных атомных столкновениях. М.:Атомиздат, 1979, 272 с.

29. Born M., Oppenheimer J.R. Zur Quantentheorie der Molekeln. Ann. der Phys., 1927, v.84, p.457.

30. Judson R., Rabitz H. On understanding the relationship between structure in the potential surface and observables in classical dynamics: a functional sensitivity analysis approach. J. Chem. Phys., 1987, v.86, p.3886.

31. Kaur A.J., Gupta R.K., Bakhshi P.S., Shanker J. Molecular properties of diatomic alkali halides. Indian Journal of Chemistry. 1983, v.22A, p.969.

32. Gupta R.K., Kaur A.J., Bakhshi P.S., Shanker J. Evaluation of molecular properties of diatomic alkali halides. Indian Journal of Physics, 1982, V.56B, p.344.

33. Van Veen N.J.A., De Vries M.S., De Vries A.E. Accurate determination of the dissociation energies of alkali halides. Chem. Phys. Lett., 1979, v.64, № 2, p.213.

34. Pathak R.K., Thakkar A.J. Very short-range interatomic potentials. J. Chem. Phys., 1987, v.87, № 4, p.2186.

35. Lacmann K. Collisional ionization. In: Advances in Chemical Physics, 1980, v.XLII, p.513.

36. Berry R.S. Interaction of vibrational and electronic motion in alkali halide molecules. J. Chem. Phys., 1957, v.27, № 6, p. 1288.

37. Слэтер Дж. Электронная структура молекул. М.:Мир, 1965.

38. Schaefer H.F. III. Methods of electronics structure theory. Plenum, NY, 1977.

39. Schaefer H.F. III. The electronic structure of atoms and molecules. Addison-Wesley Reading, Massachusetts, 1972.

40. Wall F.T., Porter R.N. Sensitivity of exchange-reaction probabilities to the potential energy surface. J. Chem. Phys., 1963, v.39, p.3112.

41. Kleyn A.W., Los J., Gislason E.A. Vibronic coupling and intersection of covalent and ionic states. Phys. Rev., 1982, v.90, p.l.

42. Мак-Вини P., Сатклиф Б. Квантовая механика молекул. М.:Мир, 1972.

43. Степанов Н.Ф. Потенциальные поверхности и химические реакции. Соросовский образовательный журнал, 1996, № 10, с.33.

44. Ghosh D.C., Biswas R. Theoretical calculation of absolute radii of atoms and ions. Part 1. The atomic radii. Int. J. Mol. Sci., 2002, v.3, p.87.

45. Wilson J.N., Curtis R.M. Dipole polarizabilities of ions in alkali halide crystals. J. Phys. Chem., 1970, v.74, № 1, p.187.

46. Mulliken R.S. Low electronic states of simple heteropolar diatomic molecules. III. Hydrogen and univalent metal halides. Physical Review, 1937, v.51, p.310.

47. Herschbach D.R. Reactive scattering. J. Faraday Discuss. Chem. Soc., 1973, v.55, p.233.

48. Ahlrichs R., Boehm H.J., Brode S., Tang K.T., Toennies J.P. Interaction potentials for alkali ion-rare gas and halogen ion-rare gas systems. J. Chem. Phys., 1988, v.88, № 10, p.6290.

49. Story T.L., Hebert A.J. Dipole moments of KJ, RbBr, RbJ, CsBr, and CsJ by the electric deflection method. J. Chem. Phys., 1976, v.64, № 2, p.855.

50. Aquilanti V., Cappelletti D., Pirani F. Range and strength of interatomic forces: dispersion and induction contributions to the bond of dications and of ionic molecules. Chem. Phys., 1996, v.209, p.299.

51. Solomonik V.G. Dissociation of alkali halides to ion pairs. J. Struct. Chem., 1978, v. 19, p.860.

52. Sato S. On a new method of drawing the potential energy surface. J. Chem. Phys., 1955, v.23, p.592.

53. Sato S. Potential energy surface of the system of three atoms. J. Chem. Phys., 1955, v.23, p.2465.

54. Mason E.A., Schamp H.W. Potential energy surface for alkali halide dissociation. Ann. Phys., NY, 1958, v.4.p.233.

55. Grice R., Herschbach D.R. Long-range configuration interaction of ionic and covalent states. Mol. Phys., 1974, v.27, p. 159.

56. Shukla A.K., Futrell J.H. Impulsive excitation in collision-induced dissociation reactions of polyatomic cations. European Mass Spectrometry, 1997, v.3, p.259.

57. Parks E.K., and Wexler S. Bond energies of the cesium halides determined by collision-induced dissociation. J. Phys. Chem., 1984, v.88, p.4492.

58. Inouye H., Noda K., Kita S. Repulsive potentials derived from beam scattering of Rb+ and Cs+ ions by rare gas atoms. J. Chem. Phys., 1979, v.71, p.2135.

59. Varshni Y.P., Bloore F.J. Rydberg function as inter-atomic potential for metals. Phys. Rev., 1963, v. 129, p. 115.

60. Kita S., Noda K., Inouye H. Repulsive potentials for Cl"-R and Br"-R (R = He, Ne and Ar) derived from beam experiments. J. Chem. Phys., 1976, v.64, p.3446.

61. Pauly H., Toennies J.P. The study of intermolecular potentials with molecular beams at thermal energies. In: Advances in Atom and Molecular Physics, 1965, v.l, p. 195.

62. Adelman S.A., Szabo A. Coulomb approximation for multiple polarizabilities and dispersion forces: analytic static polarizabilities for ground and excited state atoms. J. Chem. Phys., 1973, v.58, p.687.

63. Stwalley W.C. Semiempirical correlation of potential parameters. J. Chem. Phys., 1971, v.55, p.170.

64. Kramer H.L., Herschbach D.R. Combination rules for van der Waals force constants. J. Chem. Phys., 1970, v.53, p.2792.

65. Schatz G.C. Fitting potential energy surfaces. In: Reaction and molecular dynamics. Lecture Notes in Chemistry, Vol. 14, Eds. Lagana A. and Riganelli A., Springer, Berlin, 2000, p. 15.

66. Schatz G.C. The analytical representation of electronic potential energy surfaces. Rev. Mod. Phys., 1989, v.61, p.669.

67. Rodgers M.T., Ervin K.M., Armentrout P.B. Statistical modeling of collision-induced dissociation thresholds. J. Chem. Phys., 1997, v. 106, p.4499.

68. Russell D.H., ed. Fundamental studies of collision-induced dissociation of ions. In: Experimental Mass Spectrometry, Plenum Publishing Corp., New York, 1994, p.71.

69. Dalgarno A. Determination of polarizabilities parameters of the inert gas atoms. Adv. Phys., 1962, v.l 1, p.281.

70. Bender C.F., Davidson E.R. A natural orbital based energy calculation for helium hydride and lithium hydride. J. Phys. Chem., 1966, v.70, p.2675.

71. Levine R.D., Jortner J. Molecular energy transfer. Wiley, NY, 1976.

72. Wilson A.D., Levine, R.D. Simple models of vibrational excitation in energy transfer molecular collisions. Mol. Phys., 1974, v.27, № 5, p.l 197.

73. Schechter I., Levine R.D., Gordon R.G. Kinematic constraints in reactive collisions. J. Phys. Chem., 1991, v.95, p.8201.

74. Goldberger M.L., Watson K.M. Collision theory. John Wiley&Sons, Inc., NY, 1964.

75. Gianturco F.A., Delgado-Barrio G., Roncero O., Villarreal P. Dynamical coupling and energy transfer in weakly bound molecular complexes. Int. Rev. Phys. Chem., 1988, v.7, № l,p.l.

76. Harris R.M., Herschbach D.R. Translational-to-vibrational energy transfer: impulsive approximation for half-collisions. J. Chem. Phys., 1971, v.54, p.3652.

77. Bowman. J.M. Molecular collision dynamics. Springer-Verlag, Berlin, 1983.

78. Grice R. Effect of translational energy on reaction dynamics. J. Faraday Discuss. Chem. Soc., 1979, v.67, p.16.

79. Urena A.G. Influence of translational energy upon reactive scattering cross section: neutral-neutral collisions. Adv. Chem. Phys., 1987, v.66, p.214.

80. Agrawal P.M., Agrawal N.C., Garg V. Dynamical constraints on the angular momentum transfer in rotationally inelastic molecular collisions. J. Chem. Phys., 1985, v.83, № 9, p.4444.

81. Schatz G.C. Overview of reactive scattering. In.: Potential energy surfaces and dynamics calculations for chemical reactions and molecular energy transfer. Ed. D.G.Truhlar. Plenum Press, 1981, p.287.

82. Clary D.C. The theory of chemical reaction dynamics. Reidel, Boston, 1986.

83. Азриель B.M., Акимов B.M., Ерико Я., Русин Л.Ю. Динамика передачи энергии в столкновениях Хе с CsBr в диапазоне энергий от 0,35 до 11 эВ. М.:1989, 83 с. Деп. в ВИНИТИ 31.10.1989, №6583-В89.

84. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. Multiple impacts and energy transfer in a three-body system for noncollinear collisions. Theor. Chem. Acta, 1993, v.87, p. 195.

85. Levine R.D., Bernstein R.B. Post-threshold energy dependence of the cross section for endoergic processes: vibrational excitation and reactive scattering. J. Chem. Phys., 1972, v.56, p.2281.

86. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Траекторное моделирование элементарных процессов передачи поступательной энергии в столкновениях молекул CsBr с атомами Хе. Хим. Физика, 1990, т.9, № 11, с. 1471.

87. Азриель В.М., Акимов В.М., Русин Л.Ю. Динамика образования ионов Cs+ при столкновительной диссоциации CsBr с Хе в скрещенных молекулярных пучках. М.:1989, 73с. Деп. В ВИНИТИ 12.09.89, № 5811-В89.

88. Азриель В.М., Русин Л.Ю., Севрюк М.Б. Ориентационные эффекты в процессах передачи энергии при многократных соударениях в системе трех частиц. Хим. Физика, 1993, т. 12, № 11, с.1427.

89. Flynn G.W., Parmenter C.S., Wodtke A.M. Vibrational energy transfer. J. Phys. Chem., 1996, v.100, p.12817.

90. Никитин E.E., Осипов А.И., Уманский С.Я. Колебательно-поступательный обмен энергией при столкновениях гомоядерных двухатомных молекул. Химия плазмы. М.: Энергоатомиздат, 1989, вып.15.

91. Dickinson A.S. Non-reactive heavy particle collision calculations. Сотр. Phys. Commun., 1979, v,17,p.51.

92. Fisk G.A., Crim F.F. Single collision studies of vibrational energy transfer mechanics. Acc. Chem. Res., 1977, v. 10, p.73.

93. Gentry W.R. Vibrational excitation: classical and semiclassical methods. In: Atom-molecule collision theory. A guide for the experimentalist. Ed. R.B.Bernstein. Plenum Press, NY, 1979.

94. Rapp D., Kassal T. The theory of vibrational energy transfer between simple molecules in nonreactive collisions. Chem. Rev., 1969, v.69, p.61.

95. Armentrout P.В., Baer T. Gas-phase ion dynamics and chemistry. J. Phys. Chem., 1996, v.100, p.12866.

96. Levine R.D., Wulfman C.E. Energy transfer to a Morse oscillator. Chem. Phys. Lett., 1979, v.60, p.372.

97. Liu K. Crossed-beam studies of neutral reactions: state-specific differential cross sections. Annu. Rev. Phys. Chem., 2001, v.52, p. 139.

98. Faubel M., Toennies J.P. Scattering studies of rotational and vibrational excitation of molecules. Advances in Atomic and Molecular Physics. Eds. D.R.Bates, B.Bederson. N.Y.: Academic Press, 1977, v. 13, p.229.

99. Everitt H.O., De Lucia F.C. Rotational energy transfer in small polyatomic molecules. Advances in Atomic and Molecular Physics. Eds. Benderson B. and Walther H. San Diego, Academic, 1995, v.35,p.331.

100. Mott N.F. On the theory of excitation by collisions with heavy particles. Proc. Cambridge Philos. Soc., 1931, v.27, p.553.

101. Bates D.R. Collision processes not involving chemical reactions. Discuss. Faraday Soc., 1962, v.33, p.7.

102. Зембеков А.А., Маергойз А.И. Исследование процессов диссоциативной ионизации при столкновениях CsJ + Хе вблизи порога методом классических траекторий. Хим. Физика, 1984, т.З, с.495.

103. Азриель В.М., Акимов В.М., Грико Я., Русин Л.Ю. Траекторное моделирование столкновительной диссоциации солей щелочных металлов в молекулярных пучках. М.:1989, 56 с. Деп. в ВИНИТИ 30.05.1989, № 3574-В89.

104. Азриель В.М. Квазиклассическое траекторное моделирование рассеяния молекулярных пучков, сопровождающегося химической реакцией. Дисс. . канд. хим. наук. М.:ИНЭП ХФ РАН, 1990, 167 с.

105. Акимов В.М., Зембеков А.А., Ломакин Л.А., Маергойз А.И., Никитин Е.Е., Русин Л.Ю., Стельмах Л.С. Диссоциативная ионизация молекул CsBr при столкновениях с атомами Хе в скрещенных молекулярных пучках. Доклады АН СССР, 1980, т.253, с.633.

106. Акимов В.М., Маергойз А.И., Русин Л.Ю. Динамика элементарного процесса образования ионных пар при столкновении атома с двухатомной молекулой. VIII Всесоюзная конференция по динамике разреженных газов, 1985, Москва.

107. Parks Е.К., Hansen N.J., Wexler S. Collision-induced ion pair formation of thalium halides: absolute cross sections. J. Chem. Phys., 1973, v.58, p.5489.

108. Sheen S.H., Diplomon G., Parks E.K., Wexler S. Collision-induced ion pair formation on CsCl and Cs2Cl2. J. Chem. Phys., 1978, v.68, p.4950.

109. Ленин Л.В., Русин Л.Ю. Функция возбуждения столкновительно индуцированной атомами ртути диссоциации CsBr. Хим. Физика, 1990, т.9, с.895.

110. Lenin L.V., Rusin L.Yu. Ionic dissociation of CsBr induced by collisions with Hg: molecular beam investigation. Chem. Phys. Lett., 1990, v. 170, p.502.

111. Akimov V.M., Lenin L.V., Rusin L.Yu. Ionic dissociation of CsBr induced by collisions with Hg: molecular ions formation. Chem. Phys. Lett., 1991, v.180, p.541.

112. Parks E.K., Wagner A., Wexler S. Collision-induced ion pair formation of thalium halides: threshold behavior. J. Chem. Phys.,1973, v.58, p.5502.

113. Parks E.K., Kuhry J., Wexler S. Collision-induced ion pair formation of the thalium halides TIF and T12F2. J. Chem. Phys.,1977, v.67, p.3014.

114. Parks E.K., Inoue M., Wexler S. Collision-induced dissociation of CsJ and Cs2J2 to ion pairs by Kr, Xe and SF6. J. Chem. Phys., 1982, v.76, p. 1357.

115. Tully F.P., Lee Y.T., Berry R.S. Crossed molecular beam study of collision-induced dissociation of alkali halides. Chem. Phys. Lett., 1971, v.9, p.80.

116. Parks E.K., Sheen S.H., Wexler S. Threshold behavior for collision-induced dissociation of CsCl and Cs2Cl2 by SF6 molecules. J. Chem. Phys., 1978, v.69, p.l 190.

117. Piper L.B., Hellemans L., Sloan J., Ross J. Total cross section for formation of ions from CsBr by collision with Ar, Xe and NaBr (Ar). J. Chem. Phys., 1972, v.57, p.4742.

118. Los J., Kleyn A.W. Ion-pair formation. In: Alkali halide vapors. Eds. P.Davidovits and McFadden. Academic Press, NY, 1979, p.275.

119. Berry R.S., Cernoch T., Coplan M., Ewing J.J. Inverted population in dissociation of CsBr molecules. J. Chem. Phys., 1968, v.49, p. 127.

120. Fan H. Calculation of collisional dissociation of alkali halide molecules by classical one-dimentional models. J. Chem. Phys., 1971, v.55, p.4628.

121. Brumer P., Karplus M. Collision complex dynamics in alkali halide exchange reactions. J. Faraday Discuss. Chem. Soc., 1973, v.55, p.80.

122. Aoiz F.J., Banares L., Herrero V.J. Recent results from quasiclassical trajectory computations of elementary chemical reactions. J. Chem. Soc. Faraday Trans., 1998, v.94, № 17, p.2483.

123. Rusin L.Yu. Kinematic mechanism for the formation of two-body ionic products in collision-induced dissociation of cesium halides. J. Phys. Chem., 1995, v.99, p. 15502.

124. V.M.Azriel and L.Yu.Rusin. Channel branching in the collisions between CsCl and RbJ. In: All Russian Conference on Theoretical Chemistry (ARCTC97), Russia, Kazan, 1997, P8.

125. V.M.Azriel and L.Yu.Rusin. Branching in the ion pairs formation in the collisions of CsCl and RbJ for the energy ranging from 0,1 to 25 eV. In: 20-th International Symposium on Molecular Beams, Portugal, Lisbon, 2003, p. 139 (P9A).

126. Levine R.D., Bernstein R.B. Branching ratios in reactive molecular collisions. Chem. Phys. Lett., 1974, v.29, p.l.

127. Bernstein R.B. Branching ratios in reactive collisions: Theory vs. experiment. Int. J. Quantum Chem. Symp., 1975, v.9. p.385.

128. Akimov V.M., Azriel V.M., Rusin L.Yu. Product angular distributions in the four-atom system CsCl + RbJ for the collision energy ranging from 4 to 6 eV. In: 17-th International Symposium on Molecular Beams, France, Orsay, 1997, P.15B.

129. Азриель В.М., Акимов В.М., Русин Л.Ю. Функции возбуждения основных каналов ионной диссоциации при столкновении двухатомных молекул с ионной связью. CsCl + RbJ. Хим. Физика, 1996, т. 15, № 3, с.З.

130. Akimov V.M., Azriel V.M. and Rusin L.Yu. Dynamics of collision-induced dissociation in four-atom systems. In: 12-th European Conferenceon Dynamics of Molecular Collisions (MOLEC XII), Bristol, UK, 1998, p.28.

131. Akimov V.M., Azriel V.M., Rusin L.Yu., Sevryuk M.B. Excitation functions of cation formation in collisions of diatomic molecules with an ionic bond: CsCl + RbJ. J. Chem. Soc. Faraday Trans., 1996, v.92, № 10, p. 1683.

132. Азриель B.M., Акимов B.M., Русин Л.Ю. Траекторное моделирование динамики столкновительной диссоциации двухатомных молекул в 4-х атомных системах. Хим. Физика, 2002, т.21, № 4, с. 18.

133. Kuntz P.J. Collision-induced dissociation. II. Trajectories and models. In: Atom-molecule collision theory. Ed. R.B.Bernstein. Plenum Press, NY, 1979, p.669.

134. Эйринг Г., Лин C.E., Лин С.М. Основы химической кинетики. М.: Мир, 1983.

135. Банкер Д. Метод классических траекторий. В кн.: Вычислительные методы в физике атомных и молекулярных столкновений. М.: Мир, 1974, с.277.

136. Поллак Л.С., Еольденберг М.Я., Левицкий А.А. Вычислительные методы в химической кинетике. М.: Наука, 1984, 280 с.

137. Haselgrove С.В. The choice of initial parameters in classical trajectory technique. Math. Computation, 1961, v.15, p.323.

138. Miller W.H. Classical-limit quantum mechanics and the theory of molecular collisions. Adv. Chem. Phys., 1974, v.25, p.69.

139. Polanyi J.C. Some Concepts in Reaction Dynamics. Acc. Chem. Res., 1972, v.5, p.161.

140. Polanyi J.C., Schreiber J.L. The dynamics of bimolecular reactions. In: Physical chemistry an advanced treatise, Vol. 6A, "Kinetics of gas reactions", Chap. 6, p.383-487, Academic Press, Inc., NY, 1974.

141. Porter R.N. Molecular trajectory calculations. Annu. Rev. Phys. Chem., 1974, v.25, p.317.

142. Porter R.N., Raff L.M. Classical trajectory methods in molecular collisions. In: Dynamics of molecular collisions, Part B, Ed. W.H.Miller, Plenum Press, NY, 1976, v.l, p.l.

143. Raff L.M., Thompson D.L. Classical trajectory approach to reactive scattering. In: The theory of chemical reaction dynamics. Ed. M.Baer. CRC Press, Boca Raton, 1985.

144. Thompson D.L. Quasiclassical trajectoiy studies of reactive energy transfer. Acc. Chem. Res., 1976, v.9, p.338.

145. Thompson D.L. On a classical trajectory study of energy transfer in some atom-diatomic molecule systems. J. Chem. Phys., 1972, v.56, p.3570.

146. Kwei G.H. Classical trajectory studies of alkali halide collision dynamics. In: Alkali halide vapors: structure, spectra, and reaction dynamics. NY, 1979, p.441.

147. Скарборо Дж. Численные методы математического анализа. М.:Мир, 1934.

148. Gill S. Mathematical methods in modern physics. Proc. Cambridge Phil. Cos.,1951, v.47, p.96.

149. Хемминг P.B. Численные методы. M.: Наука, 1968, 400 с.

150. Корн В.Г., Корн Т.Г. Справочник по математике. М.: Наука, 1984.

151. Бахвалов Н.С. Численные методы. М.: Наука, 1987.

152. Rittner E.S. Binding energy and dipole moment of alkali halide molecules. J. Chem. Phys., 1951, v.l9, № 8, p. 1030.

153. Brumer P., Karplus M. Perturbation theory and ionic models for alkali halide systems. I. Diatomics. J. Chem. Phys., 1973, v.58, № 9, p.3903.

154. Kumar M., Kaur A.J., Shanker J. Spectroscopic constants of alkali halides and hydrides using modified T-Rittner model. J. Chem. Phys., 1986, v.84, p.5735

155. Patil S.H. Interionic potentials in alkali halides. J. Chem. Phys., 1987, v.86, № 1, p.313.

156. Patil S.H. Interaction of halogen anions with some cations and inert gas atoms. J. Chem. Phys., 1988, v.89, № 10, p.6357.

157. Gilbert T.L., Simpson O.C., Williamson M.A. Relation between charge and force parameters of closed-shell atoms and ions. J. Chem. Phys., 1975, v.63, № 9, p.4061.

158. Brumer P. Combination rules and correlations in repulsive potential parameters for alkali halide diatomics. Phys. Rev. A, 1974, v. 10, № 1, p.l.

159. Miller W.B., Safron S.A., Hershbach D.R. Exchange reactions of alkali atoms with alkali halides: a collision complex mechanism. Disc. Faraday Soc., 1967, № 44, p.292.

160. Safron S.A., Weinstein N.D., Herschbach D.R., Tully J.C. Transition state theory for collision complexes: Product translational energy distributions. Chem. Phys. Lett., 1972, v.12, p.564.

161. Miller W.B., Safron S.A., Herschbach D.R. Molecular beam kinetics: Four-atom collision complexes in exchange reactions of CsCl with KC1 and KJ. J. Chem. Phys., 1972, v.56, p.3581.

162. Kwei G.H., Boffardi B.P., Sun S.F. Classical trajectory studies of long-lived collision complexes. I. Reaction of К atoms with NaCl molecules. J. Chem. Phys., 1973, v.58, p. 1722.

163. Marcus R.A. On the theory of energy distributions of products of molecular beam reactions involving transient complexes. J. Chem. Phys., 1975, v.62, p. 1372.

164. Patil S.H. Alkali ion-inert gas potentials. J. Chem. Phys., 1987, v.86, № 12, p.7000

165. Nyeland C., Tang K.T., Toennies J.P. Repulsive ion-atom and ion-ion potentials from charge density overlap integrals. Chem. Phys., 1990, v. 147, p.229.

166. Gowda B.T., Benson S.W. Empirical potential parameters for alkali halide molecules and crystals, hydrogen halide molecules, alkali metal dimers, and hydrogen and halogen molecules. J. Phys. Chem., 1982, v.86, № 5, p.847.

167. Douglas T.B. Modification of the Born-Mayer potential function as applied to the crystalline alkali halides. J. Chem. Phys., 1963, v.38, p.2461.

168. Spears F.G. Repulsive potentials of atomic ions, atoms and molecules. J. Chem. Phys., 1972, v.57, p. 1842.

169. Калинин А.П. О влиянии крутизны отталкивательного потенциала взаимодействия атомных частиц на ход упругого дифференциального сечения рассеяния на малые углы. Журнал технической физики, 1982, т.52, № 3, с.549.

170. Bohm H.J., Ahlrichs R. A study of short-range repulsions. J. Chem. Phys., 1982, v.77, № 4, p.2028.

171. Калинин А.П., Леонас В.Б., Сермягин А.В. Короткодействующие силы отталкивания между атомами и молекулами атмосферных газов. В сб. Межпланетная среда. М.:Наука, 1972. с. 196.

172. Леонас В.Б. Исследования короткодействующих межмолекулярных сил. Успехи Физ. Наук, 1972, т. 107, № 1, с.29.

173. Tossi М.Р. In: Advances in Solid State Physics, 1964, v. 16, p.l.

174. Jordan K.D. Structure of alkali halides: theoretical methods. In: Alkali halide vapours. Structure, spectra and reaction dynamics. New York, 1979, p.479.

175. Varshni Y.P., Shukla R.C. New potential function for alkali halide molecules. J. Chem. Phys., 1961, v.35,p.528.

176. Redington R.L. Internuclear potential energy functions for alkali halide molecules. J. Phys. Chem., 1970, v.74, № l,p,181.

177. Hafmeister D.W., Zahrt J.D. Exchange-charge-model calculation of Born-Mayer repulsive potential in ionic gases and crystals. J. Chem. Phys., 1967, v.47, p.1428.

178. Ding A.M.G., Kirsch L.J., Perry D.S., Polanyi J.C., Schreiber J.L. Effect of changing reagent energy on reaction probability and product energy distribution. J. Faraday Discuss. Chem. Soc., 1973, v.55, p.252.

179. Семенов H.H. Цепные реакции. М.:Наука, 1986, 534 с.

180. Мотт Н., Месси Г. Теория атомных столкновений. М.:Мир, 1969.

181. Levine R.D., Bernstein R.B. Molecular reaction dynamics and chemical reactivity. NY, Oxford: Qxford Univ. Press, 1987, 535 p.

182. Gianturco F.A. The transfer of molecular energies by collisions: recent quantum treatments. Lect. Notes Chem., v.ll, (Springer, NY, 1979).

183. Никитин Е.Е. Элементарные процессы передачи энергии в системе трех атомов. Сер. Радиационная химия. Фотохимия, т.4, М.: ВИНИТИ, 1985.

184. Levine R.D. Energy transfer in molecular collisions. Israel J. Chem., 1969, v.7, p.237.

185. Secrest D. Theory of rotational and vibrational energy transfer in molecules. Annu. Rev. Phys. Chem., 1973, v.24, p.379.

186. Faubel. M. Vibrational and rotational excitation in molecular collisions. Adv. At. Mol. Phys., 1983, v.19, p.345.

187. Stevens B. Collisional activation in gases. Pergamon, Oxford, UK, 1967.

188. Yardley J.T. Introduction to molecular energy transfer. Academic Press, NY, 1975, v.28, p.171.

189. Callear A.B. An overview of molecular energy transfer in gases. Spec. Per. Rep., 1978, v.3, p.82.

190. Steinfeld J.I. Energy-transfer processes. Int. Rev. Sci. Phys. Chem., 1972, v.9, p.247.

191. Fisk G.A., Crim F.F. Single collision studies of vibrational energy transfer mechanisms. Acc. Chem. Res., 1977, v. 10, p.73.

192. Amme R.C. Vibrational and rotational excitation in gaseous collisions. Adv. Chem. Phys., 1972.

193. Azriel V.M., Rusin L.Yu. Dynamics of energy transfer for the interaction of CsCl and RbJ in cross molecular beams. In: 18-th International Conference on Molecular Energy Transfer (COMETXVIII), Spain, San Lorenzo de El Escorial, 2003, p.40 (10-TH).

194. Azriel V.M., Rusin L.Yu. Conversion of translation energy in the collisions of two diatomic molecules. In: 7-th Congress of the World Association of Theoretically Oriented Chemists (WATOC'05), Cape Town, South Africa, 2005, p. 184 (RM-P14).

195. Levine R.D. Non-reactive molecular encounters. J. Chem. Phys., 1967, v.46, p.331.

196. Verster N.F., Beiierinck H.C.W., Henrichs J.M., Van der Kam P.M.A. Multiple collisions in molecular beam scattering experiments. Chem. Phys., 1981, v.55, № 2, p.169.

197. Kafri A., Pollak E., Kosloff R., Levine R.D. Translational energy disposal in molecular collisions. Chem. Phys. Lett., 1975, v.33, p.201.

198. Moore C.B. Vibration-vibration energy transfer. Adv. Chem. Phys., 1973, v.23, p.41.

199. Pekeris C.L. The rotation-vibration coupling in diatomic molecules. Phys. Rev., 1934, v.45, p.98.

200. Кондратьев B.H., Никитин E.E., Резников А.И., Уманский С.Я. Термические бимолекулярные реакции в газах. М.: Наука, 1976.

201. Осипов А.И. Термическая диссоциация двухатомных молекул при высоких температурах. Теоретическая и экспериментальная химия, 1966, т.2, № 5, с.649.

202. Смехов Т.Д., Лосев С.А. О роли колебательного и вращательного возбуждения в процессе диссоциации двухатомных молекул. Теоретическая и экспериментальная химия, 1979, т. 15, № 5, с.492.

203. Svedung Н., Krems R., Markovic N., Nordholm S. Pair-potential model for simulation of collisional energy transfer: quantum effects and hardness dependence. Phys. Chem. Chem. Phys., 2001, v.3, p.2216.

204. Levine R.D., Manz J. The effect of reagent vibrational excitation on chemical reaction rates: an information theoretic analysis. J. Chem. Phys., 1975, v.63, p.4280.

205. De Presto А.Е., Rabitz Н. Vibrational and rotational collision processes. Adv. Chem. Phys., 1980, v.42, p.271.

206. Ioup J.W., Russek A. Vibrational-rotational excitation in atom-diatomic molecule collisions. Phys. Rev. A, 1973, v. 8, № 6, pp. 2898-2914.

207. Nyman G., Clary D.C., Levine R.D. Potential energy surface effects on differential cross sections for polyatomic reactions. Chem. Phys., 1995, v.191, p.223.

208. Rusin L.Yu. Orientation effects on collision dissociation with ion pair formation. J. Chem. Biochem. Kinetics, 1991, v. 1, № 3, p.205.

209. Schechter I., Prisant M.G., Levine R.D. A classical kinematic model for direct reactions of oriented reagents. J. Phys. Chem., 1987, v.91, p. 5472.

210. Clark A.P., Dickinson A.S., Richards D. The correspondence principle in heavy-particle collisions. Adv. Chem. Phys., 1977, v.36, p.63.

211. Ben-Shaul A. Product state distributions in chemical reactions: vibrational temperature and rotational distributions. Chem. Phys., 1973, v.l, p.244.

212. Billing. G.D. The semiclassical treatment of molecular roto/vibrational energy transfer Сотр. Phys. Rep., 1984, v.l, № 5, p.239.

213. Schatz G.C., Kupperman A. Comment on "Theory of Collisions Between an Atom and a Diatomic Molecule in the Body-Fixed Coordinate System,". J. Chem. Phys., 1979, v.70, p.3151.

214. Lendvay G., Schatz G.C. Observation of highly energetic collisions in classical trajectory studies of collisional energy transfer. J. Phys. Chem., 1990, v.94, p.8864.

215. Иванов Г.К., Саясов Ю.С. Теория рассеяния в импульсном приближении. Журнал экспериментальной и теоретической физики, 1963, т.45, с.456.

216. Иванов Г.К., Саясов Ю.С. Теория колебательного возбуждения в импульсном приближении. Доклады АН СССР, 1964, т. 154, с. 1314.

217. Иванов Т.К., Саясов Ю.С. Теория прямых атомно-молекулярных реакций. Теоретическая и экспериментальная химия, 1967, т.З, с.331.

218. Саясов Ю.С., Иванов Т.К. Теория диссоциации молекул нейтронами. Журнал экспериментальной и теоретической физики, 1961, т.40, с.513.

219. Герасименко В.И., Оксюк Ю.Д. Диссоциация двухатомных молекул при быстрых столкновениях. Журнал экспериментальной и теоретической физики, 1965, т.48, с.499.

220. Bates R., Cook С.J., Smith F.J. Classical theory of ion-molecular rearrangement collisions at high impact energies. Proc. Phys. Soc., 1964, v.83, p.49.

221. Ландау Л.Д., Лифшиц E.M. Механика. M.: Наука, 1973.

222. Kelley J.D., Wolfsberg M. Comparison of approximate translational-vibrational energy-transfer formulas with exact classical calculation. J. Chem. Phys., 1966, v.44, p.324.

223. Dickinson A.S. Non-reactive heavy particle collision calculations. Сотр. Phys. Commun., 1979, v. 17, p.51.

224. Fisanick-Englot G., Rabitz H. Studies in inelastic molecular collisions using impact parameter methods. I. Model calculations. J. Chem. Phys., 1975, v.62, p.1409.

225. Rapp D., Kassal T. The theory of vibrational energy transfer between simple molecules in nonreactive collisions. Chem. Rev., 1969, v.69, p.61.

226. Schatz G.C., Mulloney T. Collisional energy transfer in polyatomic molecules: a study of anharmonicity effects in Kr + C02. J. Chem. Phys., 1979, v.71, № 12, p.5257.

227. Никитин E.E., Овчинникова М.Я. Аналитическая модель для колебательно-вращательной передачи энергии при сильно анизотропном взаимодействии атома и двухатомной молекулы. Хим. Физика, 1989, т.5, № 3, с.413.

228. Loeser J.G., Rabitz Н., DePristo А.Е., Rohlfmg Е.А. A simple model for inelastic scattering. Am. J. Phys., 1981, v.49, p. 1046.

229. Liu M., Nordholm S., Nyman G., Davidsson J. Monte-Carlo simulation study of statistical limits for energy transfer in binary molecular collisions. Chem. Phys. Lett., 1993, v.211, № 2, p. 189.

230. Квливидзе B.A., Красильников С.С. Введение в физику атомных столкновений. Изд-во Московского университета, 1985.

231. Micha D.A., Yuan J.-M. Collision dynamics of three interacting atoms. The multiple collision expansion. J. Chem. Phys., 1975, v.63, p.5462.

232. Sohlberg K.W. Energy transfer in diatom/diatom molecular collisions. Thesis (Ph.D.) Newark, DE (United States), Univ. of Deleware, US, 1992, 121 p.

233. Bowman J.M., Gazdy В., Sun Q. A method to constrain vibrational energy in quasiclassical trajectory calculations. J. Chem. Phys.,1989, v.91, p.2859.

234. Loesch H.J., Herschbach D.R. Ballistic mechanism for vibrational and rotational energy transfer in Ar + CsJ collisions. J. Chem. Phys., 1972, v.57, p.2038.

235. Liu J., Song K., Hase W.L., Anderson S.L. Direct dynamics study of energy transfer and collision-induced dissociation: Effects of impact energy, geometry, and reactant vibrational mode in H2CO+-Ne collisions. J. Chem. Phys., 2003, v.l 19, № 6, p.3040.

236. Bradley K.S., Schatz G.C., Balint-Kurti G.G. Wavepacket methods for the direct calculation of energy transfer moments in molecular collisions. J. Phys. Chem., 1999, v. 103, p.947.

237. Schatz G.C. Collisional energy transfer as a probe of ergodicity in molecular vibrational motions. Chem. Phys. Lett., 1979, v.67, p.248.

238. Mulloney T., Schatz G.C. Classical rotational and centrifugal sudden approximations for atom-molecule collisional energy transfer. Chem. Phys., 1980, v.45, p.213.

239. King F.W., Schatz G.C. Stochastic theory of vibrational energy transfer in collinear atom-diatom collisions: The role of Non-Markovian effects. Mol. Phys., 1979, v.38, p.257.

240. Rusinek I., Roberts R.E. Semiclassical calculation for collision induced dissociation. J. Chem. Phys., 1976, v.65, № 3, p.872.

241. Schatz G.C. Quasiclassical trajectory studies of state to state collisional energy transfer in polyatomic molecules. In: Molecular collision dynamics, Topics in Current Physics, Vol. 33, Ed. Bowman J.M., Springer, Heidelberg, 1983, p.25.

242. Зембеков А.А., Маергойз А.И., Никитин Е.Е., Русин Л.Ю. Импульсная модель для диссоциативной ионизации при столкновениях атомов инертных газов с молекулами галоидов щелочных металлов. Теоретическая и экспериментальная химия, 1981, т. 17, № 5, с.579.

243. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. A hard sphere model for the reaction of two diatomic molecules with ionic bond. Chem.Phys., 1995, v. 199, p. 195.

244. Акимов B.M., Маергойз А.И., Русин Л.Ю. Динамика индуцированной столкновениями диссоциации в скрещенных молекулярных пучках. CsBr + Хе. Эксперимент. Импульсная модель. Хим. Физика, 1986, т.5, № 11, с.1514.

245. Ленин Л.В., Русин Л.Ю., Севрюк М.Б. Процессы диссоциации и комплексообразования при столкновениях молекул галогенидов цезия с атомами инертных газов и ртути: импульсная модель. М.: 1991. 189 с. Деп. В ВИНИТИ 10.12.1991, №4561-В91.

246. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. The complete T —> V, R energy conversion in hard sphere A+BC collisions. In: International Conference on Dynamics of Molecular Systems (MOLEC 2004), Nunspeet, Netherlands, 2004, p.66 (P04).

247. Benson S.W., Berend G.C., Wu J.С. Classical model for vibrational and rotational excitation of diatomic molecules by collision. I. Hard-sphere collision. J. Chem. Phys., 1963, v.38, p.25.

248. Азриель B.M., Русин Л.Ю., Севрюк М.Б. Импульсная модель взаимодействия пары двухатомных молекул с ионной связью. Хим. Физика, 1995, т. 14, № 7, с.28.

249. V.M.Azriel, L.Yu.Rusin, M.B.Sevryuk. The complete T —> V, R energy conversion in three-body collisions within the hard sphere model. J. Chem. Phys., 122, 074322 (2005).

250. Севрюк М.Б. Динамический анализ атомно-молекулярных столкновений. Дисс. . д-ра физ.-мат. наук. М.: ИНЭП ХФ РАН, 2003, 312 с.

251. Mahan В.Н. Refined impulse approximation for the collisional excitation of the classical anharmonic oscillator. J. Chem. Phys., 1970, v.52, № 10, p.5221.

252. Mahan B.H. Collinear collision chemistry. I. A simple model for inelastic and reactive collision dynamics. J. Chem. Educ., 1974, v.51, № 5, p.308.

253. Mahan B.H. Collinear collision chemistry. II. Energy disposition in reactive collisions. J. Chem. Educ., 1974, v.51, № 6, p.377.

254. Cheng M.H., Chiang M.H., Gislason E.A., Mahan B.H., Tsao C.W., and Werner A.S. Collision induced dissociation of molecular ions. J. Chem. Phys., 1970, v.52, № 12, p.6150.

255. Mahan B.H., Ruska W.E., Winn J.S. Sequential impuls model of direct reactions. J. Chem. Phys., 1976, v.65, p.3888.

256. Shin H.K. Effects of the multiplicity of impacts on vibration translation energy transfer in collinear collisions. J. Chem. Phys., 1975, v.62, № 10, p.4130.

257. Shin H.K. Vibrational energy transfer. In: Dynamics of molecular collisions, Part A, Modem Theoretical Chemistry, edited by W.H. Miller (Plenum, NY, 1976), v.l, p. 131.

258. Shin H.K. Vibrational energy transfer for an impulsive ВС + A collisions. Chem. Phys. Lett., 1969, v.4, № 5, p.297.

259. Jepsen D.W., Hirschfelder J.O. Idealized theory of the recombinations of atoms by three-body collision. J. Chem. Phys., 1959, v.30, № 4, p. 1032.

260. Eyring H., Polanyi M. Über einfache Gasreaktionen. Z. Phys. Chem. B., 1931, B.12, s.279.

261. Rebick C., Levine R.D., and Bernstein R.B. Energy requirements and energy disposal: Reaction probability matrices and a computational study of a model system. J. Chem. Phys., 1974, v.60, p.4977.

262. Hijazi N.N., Laidler K.J. Dynamics of collinear A + BC systems. J. Chem. Phys., 1973, v.58, p.349.

263. Kuntz P.J. Analytical properties of a direct interaction model for gas-phase chemical reactions A + BC —> AB + C. Trans. Faraday Soc., 1970, v.66, p.2980.

264. Marcus R.A. On the analytical mechanics of chemical reactions: Classical mechanics of linear collisions. J. Chem. Phys., 1966, v.45, p.4500.

265. Fan H. Calculation of collisional dissociation of alkali halide molecules by classical one-dimensional models. J. Chem. Phys., 1971, v.55, p.4628.

266. Truhlar D.G. Quasiclassical predictions of final vibrational state distributions in reactive and nonreactive collisions. Int. J. Quantum Chem. Symp., 1976, v. 10, p.239.

267. Levine R.D., Bernstein R.B. Energy disposal and energy requirements for elementary chemical reactions. J. Faraday Discuss. Chem. Soc., 1973, v.55, p. 100.

268. Levine R.D., Bernstein R.B. Energy disposal and energy consumption in elementary chemical reactions: The information theoretic approach. Acc. Chem. Res., 1974, v.7, p.393.

269. Anlauf K.G., Kuntz P.J., Maylotte D.H., Pacey P.D., Polanyi J.C. Energy distribution among reaction products. Discuss. Faraday Soc., 1967, v.44, p. 183.

270. Bernstein R.B., Rulis A.M. Translational energy dependence of product energy and angular distribution for the K + CH3I —» KI + CH3 reaction. J. Faraday Discuss. Chem. Soc., 1973, v.55, p.293.

271. Anlauf K.G., Charters P.E., Home D.S., McDonald R.G., Maylotte D.H., Polanyi J.C., Skrlac W.J., Tardy D.C., Woodall K.B. Translational energy distribution in the products of some exothermic reactions. J. Chem. Phys., 1970, v.53, p.4091.

272. Kuntz P.J., Nemeth E.M., Polanyi J.C., Rosner S.D., Young C.E. Energy distribution among products of exothermic reactions. II. Repulsive, mixed, and attractive energy release. J. Chem. Phys., 1966, v.44, p.l 168.

273. Azriel V.M., Rusin L.Yu. The effect of internal energy of molecules CsCl and RbJ on dynamics of their collision-induced dissociation. In: 6-th International Conference on Chemical Kinetics, Gaithersburg, Maryland, USA, 2005, C20.

274. Azriel V.M., Rusin L.Yu. The effect of vibrational and rotational excitation of the reagents on dynamics of the interaction in the system CsCl + RbJ. In: 9-th European Conference on Atoms, Molecules & Photons, (ECAMP IX), Crete, Greece, 2007, Mo5-2.

275. Pack R.T., Walker R.B., Kendrick B.K. Three-body collision to recombination and collision-induced dissociation. I. Cross sections. J. Chem. Phys., 1998, v. 109, p.6701.

276. Бенсон С. Основы химической кинетики. Пер. с англ., М.:Мир, 1964.

277. Мелвин-Хьюз Э.А. Физическая химия. Пер. с англ., М., Изд-во Иностр. литературы, 1962.

278. Кондратьев В.Н., Никитин Е.Е. Кинетика и механизм газофазных реакций. М.:Наука, 1974.

279. Смит Ф. Тройные столкновения и скорости тримолекулярных реакций. В сб. "Кинетические процессы в газах и плазме", М.:Атомиздат, 1972, с.277.

280. Pack R.T., Walker R.B., Kendrick B.K. Three-body collision contributions to recombination and collision-induced dissociation. II. Kinetics. J. Chem. Phys., 1998, v. 109, p.6714.

281. Smith D., Adams N.G. Ionic reactions in thermal plasmas. J. Chem. Soc. Faraday Trans. II., 1987, v.83.p,149.

282. Bates D.R. Dissociative recombination. Phys. Rev., 1950, v.78, p.492.

283. Гейдон А.Г., Вольфгард Х.Г. Пламя, его структура, излучение и температура. М.: Металлургиздат, 1959, 333 с.294. "Энциклопедия низкотемпературной плазмы", под ред. В.Е.Фортова, т.1., М.:Наука, 2000.

284. Pack R.T., Walker R.B., Kendrick B.K. Mechanisms of recombination and collision-induced dissociation. Theoretical chemistry and molecular physics, 1997, v. 12, p.l.

285. Wei G.W., Saman Alavi, Snider R.F. Equal reaction rates for all recombination pathways. J. Chem. Phys., 1997, v. 106, № 4, p. 1463.

286. Pack R.T., Butcher E.A., Parker G.A. Accurate three-dimensional quantum probabilities and collision lifetimes of the H + 02 combustion reaction. J. Chem. Phys., 1995, v,102,№ 15, p.5998.

287. Azriel V.M., Rusin L.Yu. Trajectory simulation of the recombination reaction Cs++ Br" + Xe —> CsBr* + Xe. In: 14-th Symposium on Atomic, Cluster and Surface Physics (SASP 2004), Italy, La Thuile, 2004, PC-4.

288. Azriel V.M., Rusin L.Yu. Quasiclassical trajectory simulation of recombination process Cs + Br" + Hg —> CsBr + Hg. In: 8-th European Conference on Atomic and Molecular Physics (ECAMP VIII), France, Rennes, 2004, 11-3 (B-7-3).

289. Azriel V.M., Rusin L.Yu. Dynamic features of recombination of ionic pair Cs+ Br" at the presence of atoms Xe and Hg. In: 21-th International Symposium on Molecular Beams, Hersonissos, Crete, Greece, 2005, p.53 (PI 1).

290. D by collisions with slow electrons. J. Chem. Phys., 1987, v.86, p.726.

291. Graham W.G., Fritsch W., Hahn Y., Tanis J.A. Recombination of atomic ions. NATO Science series: B, Springer US, 1992, 356 p.

292. Nahar S.N. Electron-ion recombination rate coefficients, photoionization cross sections, and ionization fractions for astrophysically abundant elements. II. Oxygen ions. Astrophys. J. Suppl., 1999, v.120, p.131.

293. Larsson M. Experimental studies of the dissociative recombination of Нз+. Phil. Trans. R. Soc. Lond. A, 2000, v.358, p.2433.

294. Иванов B.A. Диссоциативная рекомбинация молекулярных ионов в плазме инертных газов. Успехи физических наук, 1992, т.162, № 1, с.35.

295. Strasser D., Lammich L., Krohn S., Lange M., Kreckel H., Levin J., Schwalm D., Vager Z., Wester R., Wolf A., Zajfman D. Two- and three-body kinematical correlation in the dissociative recombination of H3+. Phys. Rev. Lett., 2001, v.86, № 5, p.779.

296. Guberman S.L. Ab initio studies of dissociative recombination. In: Dissociative recombination: theory, experiment and applications. Eds. Mitchell J.B.A., Guberman S.L. Word Scientific, Singapore, 1989, p.45.

297. Gousset G., Sayer В., Berlande J. Electron Cs+ ion recombination in the presence of neutral helium atoms. Phys. Rev., A, 1977, v.16, № 3, p. 1070.

298. Nahar S.N., Pradhan A.K. Unified treatment of electron-ion recombination in the close coupling approximation. Phys. Rev. A, 1994, v.49, p. 1816.

299. McCall B.J. Recombination cool and fast. Nature, 2006, V.440, № 9, p.l 57.

300. Whitten B.L., Morgan W.L., Bardsley J.N. Monte-Carlo calculations of two- and three-body ionic recombination. J. Phys. В., 1982, v. 15, p.319.

301. Mitchell J.B.A. The dissociative recombination of molecular ions. Phys. Rep., 1990, v.186, p.215.

302. Howard R.E., Planck T.L., Trussell S.R., Saadevandi B. Quasiclassical trajectory calculation of argon-argon recombination rates. Chem. Phys. Lett., 1987, v. 142, p.33.

303. Feibelman P.J. Numerical study of the three-body ion-ion recombination rate. J. Chem. Phys., 1965, v.42, p.2462.

304. Shepler B.C., Balabanov N.B., Peterson K.A. Hg + Br —> HgBr recombination and collision-induced dissociation dynamics. J. Chem. Phys. 127, 164304 (2007).

305. Mitchell J.B.A., Guberman S.L. Dissociative recombination: theory, experiment and applications. Singapore, World Scientific, 1989.

306. Clarke A.G., Burns G. Trajectory studies of atomic recombination reactions. J. Chem. Phys., 1971, v.55,№ 10,p.4717.

307. Глестон С, Лейдер К., Эйринг Г. Теория абсолютных скоростей реакций. М., Изд. Ин. Лит., 1948,583 с.

308. Семиохин И.А., Страхов Б.В., Осипов А.И. Кинетика химических реакций. М. Изд. МГУ, 1995,347 с.

309. Azriel V.M., Rusin L.Yu. Influence of weight of third body on probability of recombination in the system Cs+ + Br" + R —> CsBr* + R. In: 6-th International Conference on Chemical Kinetics, Gaithersburg, Maryland, USA, 2005, F20.

310. Muckerman J.Т., Faist М.В. Rate constants from Monte Carlo quasi-classical trajectory calculations. A procedure for importance sampling. J. Phys. Chem., 1979, v.83, № 1, p.79.

311. Aoiz F.J., Banares L., Herrero V.J. Classical versus quantum mechanical calculations of the dynamics of elementary reactions: From state-resolved cross sections to rate constants. Advances in Classical Trajectory Methods, 1998, v.3, p.121.

312. Lammich L., Kreckel H., Krohn S., Lange M., Schwalm D., Strasser D., Wolf A., Zajfman D.

313. Breakup dynamics in the dissociative recombination of Нз+ and its isotopomers. Rad. Phys. Chem., 2003, v.68,p,175.

314. Guberman S.L. Dissociative recombination mechanisms. In: Dissociative recombination of molecular ions with electrons. Ed. Guberman S.L., Kluver, Plenum Academic Press, NY, 2003, p.l.

315. Nahar S.N., Pradhan A.K. Electron-ion recombination in the close coupling approximation. Phys. Rev. Lett., 1992, v.68, p.1488.

316. Scrinzi A., Elander N., Wolf A. Laser induced recombination to excited stated of hydrogen-like ions. Z. Phys. D, 1995, v.34, p.l85.

317. Krohn S., Lange M., Grieser M., Knoll L., Kreckel H., Levin J., Repnow R., Schwalm D., Wester R., Witte P., Wolf A., Zajfman D. Rate coefficients and Final States for the Dissociative Recombination of LiH+. Phys. Rev. Lett., 2001, v.86, p.4005.

318. Kreckel H., Mikosch J., Wester R., Glosik J., Plasil R„ Motsch M., Gerlich D., Schwalm D., Zajfman D., Wolf A. Towards state selective measurements of the Нз+ dissociative recombination rate coefficient. J. Phys.: Conf. Ser.,2005, v.4, p. 126.

319. Valcu В., Schneider I.F., Rault F., Stromholm S., Larsson M., Suzor-Weiner A. Rotational effects in low energy dissociative recombination of diatomic ions. Eur. Phys. J. D„ 1998, v.l, p.78.

320. Larsson M., Thomas R. Three-body reaction dynamics in electron-ion dissociative recombination. Phys. Chem. Chem. Phys., 2001, v.3, p.4471.

321. Colavecchia F.D., Mrugala F., Parker G.A., Pack R.T. Accurate quantum calculations on three-body collisions in recombination and collision-induced dissociation. II. An improved propagator. J. Chem. Phys., 2003, v.l 18, p. 10387.

322. Geppert W.D., Thomas R.D., Ehlerding A., Semaniak J., Ósterdahl F., Af Ugglas M., Djuric N., Paal A., Larsson M. Extraordinary branching ratios in astrophysically important dissociative recombination reactions. Faraday Discuss., 2004, v. 127, p.425.

323. Azriel V.M., Rusin L.Yu. Dynamics of ionic recombination in the system Cs+ + Br" + R —> CsBr* + R, where R = Hg, Xe and Kr. In: 15-th Symposium on Atomic, Cluster and Surface Physics (SASP 2006), Austria, Obergurgl, 2006, p. 156.

324. Азриель B.M., Кабанов Д.Б., Колесникова Л.И., Русин Л.Ю. Динамика рекомбинации ионов в низкотемпературной плазме. Известия Академии Наук, серия Энергетика, 2007, № 5, с.50.

325. Bunker D.L. Theory of elementary gas reaction rates. Pergamon Press, Inc., Elmsford, NY, 1966.

326. Laidler K.J. Theories of chemical reaction rates. McGraw-Hill Book Company, NY, 1969.

327. Casavecchia P. Chemical reaction dynamics with molecular beams. Rep. Prog. Phys., 2000, v.63, p.355.

328. Larsson M. Molecular ion recombination in merged beams: experimental results on small systems and future perspectives. J. Phys. Conf. Ser., 2005, v.4, p.50.

329. Casavecchia P., Balucani N., Volpi G.G. Crossed beam studies of reaction dynamics. Annu. Rev. Phys. Chem., 1999, v.50, p.347.

330. Scoles G. Atomic and molecular beam methods. Oxford Univ. Press, NY, 1988.

331. Rusin L.Yu., Sevryuk M.B., Toennies J.P. Comparison of experimental time-of-flight spectra of the HF products from the F+H2 reaction with quantum mechanical calculations. J. Chem. Phys., 2005, 122 (134314).

332. Fernandez-Ramos A., Miller J.A., Klippenstein S.J., Truhlar D.G. Modeling the kinetics of bimolecular reactions. Chem. Rev., 2006, v.106, p.4518.

333. Herschbach D.R. Analysis of reactive scattering in crossed molecular beams. J. Chem. Phys., 1960, v.33, p.1870.

334. Herschbach D.R. Chemical reactions in crossed molecular beams. The Vortex, 1961, v.22, p.348.

335. Herschbach D.R. Reactive collisions in crossed molecular beams. Disc. Faraday Soc., 1962, v.33, p. 149.

336. Ramsey N.F. Molecular beams. Oxford University Press, NY, 1986.

337. Grice R. Reactions studied by molecular beam technique. Sp. Per. Rep., 1981, v.4, p.l.

338. Bernstein R.B. Molecular reaction dynamics by the crossed molecular beam technique: a few recent examples. Pure & Appl. Chem., 1978, v.50, p.781.

339. Buss R.J., Lee Y.T. Crossed beam studies related to gas kinetics. J. Phys. Chem., 1979, v.83, № 1, p.34.

340. Lee Y.T., Shen Y.R. Studies with crossed laser and molecular beams. Phys. Today, 1980, v.33, № 11, p.52.

341. Toennies J.P. Molecular beam scattering experiments on elastic, inelastic and reactive collisions. In: Phys. Chem. Adv. Treat. Vol.4. Kinetics of gas reactions. Academic Press, NY, 1972, p.227.

342. Акимов B.M. Динамика индуцированного столкновениями образования ионных пар в скрещенных молекулярных пучках. Дисс. . канд. хим. наук. М.: ИНЭП ХФ, 1990, 156 с.

343. Herschbach D.R. Chemical dynamics. Science Year, 1970, v. 1971, p.281.

344. Mariella R.P., Jr., Neoh S.K., Herschbach D.R., Klemperer W. Vibrational relaxation in seeded supersonic alkali halide beams. J. Chem. Phys., 1977, v.67, p.2981.

345. Rajan M.S., Gislason E.A. Determination of cesium ion-rare gas potentials from total cross section measurements. J. Chem. Phys., 1983, v.78, № 5, p.2428.

346. Lobb W.S., McCartan D.G. The effective range of interatomic potentials in collision broadening. J. Phys. В., 1988, v.21, № 24, p.723.

347. Futrell J.H. Crossed molecular beam studies of state-to-state reaction dynamics. Adv. Chem. Phys., 1992, v.82, p.501.

348. Азриель B.M., Акимов B.M., Грико Я., Русин Л.Ю. Рассеяние ионов Cs+ при диссоциации CsBr в скрещенных молекулярных пучках. Сравнение экспериментальных данных с траекторной моделью. Хим. Физика, 1990, т.9, № 10, с.1306.

349. Азриель В.М., Русин Л.Ю., Севрюк М.Б., Тоннес Я.П. Некоторые особенности топологии ППЭ и их связь с динамикой химического канала взаимодействия атомов фтора с молекулярным водородом. Журн. Физ. Хим., 1998, т.72, № 8, с. 1351.

350. Леонас В.Б. Изучение элементарных химических реакций с помощью метода молекулярного пучка. Успехи Химии, 1966, Вып.12, т.35, с.2105.

351. Леонас В.Б. Современное состояние и некоторые новые результаты метода молекулярного пучка. Успехи Физ. Наук, 1972, т. 107, № 2, с.287.

352. Леонас В.Б. Межмолекулярные взаимодействия и столкновения атомов и молекул. Итоги науки и техники. Физика атома и молекулы, оптика, магнитный резонанс, т.1, 1980.

353. Grice R. Molecular beam studies of chemical reaction dynamics. Endeavour. New Series, 1993, v. 17, № 4, p. 173.

354. Liu K. Recent advances in crossed-beam studies of bimolecular reactions. J. Chem. Phys., 2006, 125 (132307).

355. Herschbach D.R. Chemical reaction dynamics and electronic structure. In: The chemical bond: structure and dynamics. Ed. Zewail A. Academic Press, 1992, p. 175.

356. Herschbach D.R. The shape of molecular collisions. In: Science and Society. Ed. Moskovits M. House of Anansi Press, Ltd., Concord, Ontario, 1995, p. 13.

357. Liu K. Excitation functions of elementary chemical reactions: a direct link from crossed-beam dynamics to thermal kinetics. Int. Rev. Phys. Chem., 2001, v.20, № 2, p. 189.

358. Datz S., Herschbach D.R., Taylor E.H. Collision mechanics in crossed maxwellian molecular beams. J. Chem. Phys., 1961, v.35, p.1549.

359. Русин Л.Ю. Молекулярные пучки новый путь исследования элементарных химических процессов в неравновесных условиях. В сб.: Химия плазмы. Вып.1. Под ред. Б.М.Смирнова. М.: Атомиздат, 1974, с.260.

360. Bernstein R.B. Chemical dynamics via molecular beams and laser techniques. New York, 1982,562 p.

361. Polanyi J.C. Molecular beam scattering. J. Faraday Discuss. Chem. Soc., 1973, v.55, p.389.

362. Lee Y.T. Crossed molecular beam studies and dynamics of decomposition of chemically activated radicals. Berichte der Bunsen-Gesellschaft, 1974, v.78, № 2, p. 135.

363. Balucani N., Capozza G., Leonori F., Segoloni E., Cassavechia P. Crossed molecular beam reactive scattering: from simple triatomic to multichannel polyatomic reactions. Int. Rev. Phys. Chem., 2006, v.25, p. 109.

364. Fluendy M.A.D., Lawley K.P. Chemical applications of molecular beam scattering. Chapman&Hall, London, 1973.

365. Clough P.N., Geddes J. Some aspects of chemical interaction of molecular beams. J. Phys. E: Sci. Instrum., 1981, v,14,p.519.

366. Karplus M. Molecular beams and Reaction kinetics. Acad. Press, NY, 1970.

367. Johnston H.S. Gas phase reaction rate theory. The Ronald Press Company, NY, 1966.

368. Herschbach D.R. Reactive scattering in molecular beams. Adv. Chem. Phys., 1966, v.10, p.319.

369. Herschbach D.R. Atomic and molecular beams in chemical physics: a continuing odyssey. In: Atomic and molecular beams. Ed. Campargue R. Springer-Verlag, Berlin, 2000, p.3.

370. Birely J.H., Herm R.R., Wilson K.R., Herschbach D.R. Molecular beam kinetics: Reactions of K, Rb and Cs with Br2 and J2. J. Chem. Phys., 1967, v.47, p.993.

371. Grice R., Empedocles P.B. Molecular beam kinetics: Reactions of K, Rb and Cs with Cl2. J. Chem. Phys., 1968, v.48, p.5352.

372. Wilson K.R., Kwei G.H., Norris J.A., Herm R.R., Birely J.H., Herschbach D.R. Reactive scattering in molecular beams: evidence for a stripping mechanism in reactions of alkali atoms with halogens. J. Chem. Phys., 1964, v.41, p.l 154.

373. Kauzmann W. Kinetic theory of gases. Benjamin, NY, 1966.

374. Karplus M., Porter R., Sharma R.D. Exchange reaction with activation energy. I. Simple barrier potential for (Н,Нг). J. Chem. Phys., 1965, v.43, p.3259.

375. Базь А.И., Зельдович Я.Б., Переломов A.M. Рассеяние, реакции и распады в нерелятивистской квантовой механике. М.:Наука, 1966, 340 с.

376. Лайт Дж., Росс Дж., Шулер К. Сечения реакций, константы скорости и микроскопическая обратимость. Кинетические процессы в газах и плазме. М.: Атомиздат, 1972, с.241.

377. Фейнман, Лейнон, Сэндс. Фейнмановские лекции по физике. Том 4. Кинетика. Теплота. Звук. М.: Мир, 1967, 261 с.

378. Базаров И.П., Ееворкян Э.В., Николаев П.Н. Неравновесная термодинамика и физическая кинетика. М.: Изд. МЕУ, 1989. 240 с.

379. Романовский Б.В. Основы химической кинетики. М.: Изд. Экзамен, 2006. 416 с.

380. Gallavotti G., Cohen E.G.D. Dynamical ensembles in nonequilibrium statistical mechanics. Phys. Rev. Lett., 1995, v.74, p.2694.

381. Gallavotti G., Cohen E.G.D. Dynamical ensembles in stationary states. J. Stat. Phys., 1995, v.80, p.931.

382. Kurchan J. Fluctuation theorem for stochastic dynamics. J. Phys. A, 1998, v.31, p.3719.

383. Gibbs J.W. Elementary principles in statistical mechanics, developed with especial reference to the rational foundation of thermodynamics. Yale Univ. Press, 1902.

384. Jaynes E.T. Information theory and statistical mechanics. Phys. Rev., 1957, v. 106, p.620.

385. Manousiouthakis V.I., Deem M.W. Strict detailed balance is unnecessaiy in Monte Carlo simulation. J. Chem. Phys., 1999, v.l 10, p.2753.

386. Dirac P.A.M. The conditions for statistical equilibrium between atoms, electrons and radiation. Proc. Rol. Soc. A, 1924, v.106, p.581.

387. Eyink G.L. Action principle in statistical dynamics. Prog. Theor. Phys. Suppl., 1998, v.130, p.77.

388. Newman M.E.J., Barkema G.T. Monte Carlo methods in statistical physics. Clarendon Press, 1999.

389. Hoover W.G. Time reversibility, computer simulation, and chaos. Advanced Series in Nonlinear Dynamics, 1999, v. 13, 280 p.

390. Alberty R.A. Principle of detailed balance in kinetics. J. Chem. Edu., 2004, v.81. № 8, p.1206.

391. Маергойз А.И., Никитин Е.Е., Русин Л.Ю. Динамика образования ионов при столкновительной диссоциации двухатомных молекул. В сб. "Химия плазмы", М.: Энергоиздат, 1985, т. 12, с.З.

392. Русин Л.Ю. Динамика образования ионных пар, индуцированного столкновениями тяжелых частиц. Изв. Академии наук. Энергетика, 1997, № 1, с.41.

393. Wexler S., Parks Е.К. Molecular beam studies of collisional ionisation and ion-pair formation. Ann. Rev. Phys. Chem., 1979, v.30, p. 179.

394. Азриель В.M., Акимов В.М., Ленин Л.В., Русин Л.Ю. Рассеяние ионов Cs+ в процессах диссоциации молекул CsBr. Отчет ИНЭП ХФ АН СССР, Москва, 1990, № 0290.0037121.

395. Русин Л.Ю. Рассеяние атомов и молекул при диабатических и неадиабатических взаимодействиях: Дисс. д-ра физ.-мат. наук. М.: ИНЭП ХФ РАН, 1991, 425 с.

396. Parks Е.К., Pobo L.G., Wexler S. Collision-induced dissociation of the cesium halides. J. Chem. Phys., 1984, v.80, p.5003.

397. Азриель В.M., Акимов В.M., Грико Я., Русин Л.Ю. Динамика реакции образования атомных ионов при диссоциации молекул CsBr в столкновениях с ксеноном. Хим. Физика, 1990, т.9, № 11, с. 1463.

398. Азриель В.М., Акимов В.М., Русин Л.Ю. Рассеяние ионов Cs+ при диссоциации CsBr в скрещенных молекулярных пучках. Эксперимент. Хим. Физика, 1990, т.9, № 9, с.1224.

399. Tully F.P., Cheung N.H., Haberland H., Lee Y.T. Mass and orientation effects in dissociative collisions between rare gas atoms and alkali halide molecules. J. Chem. Phys., 1980, v.73, № 9, p.4460.

400. Lenin L.V., Rusin L.Yu. Ionic dissociation of CsBr induced by collisions with Hg: trajectory simulation. Chem. Phys. Lett., 1990, v. 175, p.608.