Возмушенные колебательные состояния многоатомных молекул тема автореферата и диссертации по физике, 01.04.05 ВАК РФ

Гавва, Светлана Павловна АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Саратов МЕСТО ЗАЩИТЫ
2011 ГОД ЗАЩИТЫ
   
01.04.05 КОД ВАК РФ
Диссертация по физике на тему «Возмушенные колебательные состояния многоатомных молекул»
 
Автореферат диссертации на тему "Возмушенные колебательные состояния многоатомных молекул"

Гавва Светлана Павловна

ВОЗМУЩЕННЫЕ КОЛЕБАТЕЛЬНЫЕ СОСТОЯНИЯ МНОГОАТОМНЫХ МОЛЕКУЛ

01.04.05-Оптика

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

г 8 СЕН 2011

Саратов - 2011

4852773

Работа выполнена в Саратовском государственном техническом университете и в Саратовском государственном университете им. Н.Г.Чернышевского

Научный консультант: доктор физико-математических наук, профессор

Березин Кирилл Валентинович

Официальные оппоненты: доктор физико-математических наук, профессор

Бабков Лев Михайлович

доктор физико-математических наук, профессор Малыханов Юрий Борисович

доктор физико-математических наук, профессор Уманский Игорь Маркович

Ведущая организация: Учреждение Российской академии наук

Институт оптики атмосферы им. В.Е.Зуева Сибирского отделения РАН

Защита состоится 27 октября 2011 г. в 15 час. 30 мин. на заседании диссертационного совета Д 212.234.01 при Саратовском государственном университете имени Н.Г. Чернышевского по адресу: 410012, г. Саратов, ул. Астраханская, 83, СГУ, Физический факультет.

С диссертацией можно ознакомиться в зональной научной библиотеке имени В.А. Артисевич Саратовского государственного университета.

Автореферат разослан ■ 2011 г.

Ученый секретарь диссертационного совета

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

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

Большой вклад в разработку теории колебаний многоатомных молекул внесли М.А. Ельяшевич, Е. Вильсон, М.В. Волькенштейн, М.А.Ковнер, Л.М.Свердлов, Л.А.Грибов, Б.И.Степанов, В.А.Дементьев, В.И.Баранов. Дальнейшее развитие колебательной спектроскопии определяется разработкой теоретических методов, которые применялись для детального описания структуры низколежащих колебательных состояний, возбужденных колебательно-вращательным взаимодействием, и развитием теории интенсивностей многофотонных колебательных переходов. Результатам этих исследований посвящены работы многих авторов: С. Amata, H.H. Nielsena, G.Tarrago, A.Barbe, L.R. Browna, C. Camy-Peyret, J.-M.Flauda, T. Oka, M. Lóete, W. M. Smita, J. L. Teffo, J.K. G. Watson, M. Р. Алиева, B.M. Михайлова, Т.Д. Коломийцовой, Н.Д. Щепкина, А.Д. Быкова, Ю.С. Ма-кушкина, В.И. Перевалова, В.И. Старикова, В.Г. Тютерева и др.

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

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

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

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

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

В соответствии с поставленной целью решались следующие основные задачи:

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

- решение прямой спектроскопической задачи возмущенных колебаний молекул аналитическими методами дискретного спектрального анализа матричной теории возмущений;

- исследование влияния возмущений колебаний на изменения частот, элементов формы, внутримолекулярные, колебательно-вращательные и спектроскопические параметры многоатомных молекул;

- создание теоретической модели высоких обертонных колебаний многоатомных молекул, осуществление вычислительного эксперимента электрооптических параметров и интегральных интенсивно-стей в ИК спектрах поглощения высоких обертонов;

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

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

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

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

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

- впервые для решения прямой спектроскопической задачи возмущенных колебаний молекул применяются методы дискретного анализа матричной аналитической теории возмущений;

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

- применен процесс редукции матричной теории возмущений для описания более тонких эффектов и особенностей влияния возмущений высоких порядков на смещения частот, изменения элементов формы колебаний, внутримолекулярные и спектроскопические параметры;

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

- впервые для молекулы СН4 выполнены прямые расчеты ©?' и ©,,,

®к -параметров операторов центробежных дипольных моментов Мт и

М05 в основном и возбужденных колебательных состояниях, учет которых способствует активизации запрещенных вращательных переходов; для одноквантовых колебательных переходов вычислены Э - параметры операторов эффективного дипольного момента Мп и М13, в матричных элементах двухквантовых переходов рассчитаны ангармонические электрооптические параметры оператора М21 и 0 -

параметры оператора кориолисова типа М22; вычисленные значения параметров в операторах эффективного дипольного момента определяют вклады в интенсивности одно - и двухквантовых колебательных переходов, заимствованные из вращательных спектров за счет учета колебательно-вращательных взаимодействий;

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

- разработана методика решения фундаментальной задачи колебательной спектроскопии по определению ангармонических ЭОП многоатомных молекул, которая объединяет квантово-химическую модель электронного строения многоатомных молекул, реализованную с использованием полуэмпирического метода СИОО/2, метод сплайн -

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

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

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

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

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

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

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

Основные положения и результаты, выносимые на защиту

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

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

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

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

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

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

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

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

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

Апробация работы. Результаты исследований были доложены на: XX Всесоюзном Съезде по спектроскопии (Киев, 1988); XXI, XXII, XXIII Съездах по спектроскопии (Звенигород, 1995, 2001, 2005 ); III, IV, V, VI, VII, IX, X Всесоюзных симпозиумах по молекулярной спектроскопии высокого и сверхвысокого разрешения (Томск, 1976, 1978, 1980, 1982, 1985, 1989, 1991); Всесоюзных конференциях по квантовой химии (Новосибирск, 1990, Казань, 1991); Межвузовских конференциях «Спектроскопия и физика молекул» (Саратов, 1997, 1998, 1999); XXIV European Congress on Molecular Spectroscopy (Prague, 1998 ); IV Internationals Conferece on Advanctd Vibrational Spectroscopy (Turku, 2001); IV, V, VI, VII, VIII, IX, X, XI, XII, XIII Международной научной школе по оптике, лазерной физике и биофизике (Саратов, 2000, 2001, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010 ); XIII, XX Международной школе-семинаре «Спектроскопия молекул и кристаллов» (Киев, 2003, 2011); XI, XII, XIII, XIV, XV, XVI Международных симпозиумах по молекулярной спектроскопии высокого и сверхвысокого разрешения (Санкт-Петербург, 1994, Нижний Новгород, 1996, Томск, 1999, 2003, 2006, 2009 ); VII Международном семинаре «Физико-математическое моделирование систем» (Воронеж, 2010); XXIV Съезде по спектроскопии ( Москва, Троицк, 2010); Международном симпозиуме «Атмосферная радиация и динамика» (Санкт-Петербург -Петродворец, 2011).

Публикации. Список публикаций по теме диссертации содержит 82 наименования, из которых 15 статей опубликованы в рецензируемых журналах, 31 статья в других научных изданиях, 35 тезисов докладов на конгрессах, съездах, симпозиумах и конференциях, одно авторское свидетельство на разработанную расчетную программу.

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

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

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

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

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

взаимодействия, ^ - матрица силовых постоянных, Ь - матрица элементов формы колебаний, Л - диагональная матрица собственных значений А, =4л-2сг(»12, пропорциональных квадратам частот нормальных колебаний (0,(5=1, 2,... ,ЗМ-6), N - число атомов в молекуле. Фундаментальное уравнение теории колебаний молекул, переход от него к вековому уравнению и способы его решения определяют целесообразность применения в молекулярной спектроскопии при исследовании возмущенных колебаний аналитических методов матричной теории.

Первый применяемый метод основан на теории Дж.Уилкинсона, в которой проблема возмущения частот и элементов формы колебаний решается для матрицы вида Б (е) = (Э+е здесь невозмущен-

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

Л^е) = Я,- + + £+... при / = 1,2, ...ЗАГ - 6 , А,-(г) А,- для е -» О и соответствующих элементов в разложении тензора формы колебаний ЬДг) = Ь, + (й^,, + ¿г/я'|, + +... + + + ...)Ь„, причем каждая скобка также является сходящимся степенным рядом.

В работе решена прямая спектроскопическая задача возмущенных колебаний молекул при аналитическом представлении матрицы динамических коэффициентов 0(е)= Х£и£>(п). Начиная со второго порядка тео-

л = 0

рии возмущений, получены новые формулы для вкладов в изменения коэффициентов к„х' и которые имеют вид

=/?«(л)+32"6(А/и-%л' + - + Р а (1)' л»-1) > ^ ]=\ '

. (2)

ЗN-6

Ь',"' = X ; , (1 = 1,2,..¿N-6). (3)

] = 1

В формулах (1) - (3) учитываются вклады всех членов до п-го порядка аналитической теории возмущений. Предложенный аналитический метод прямых расчетов изменений частот и элементов формы возмущенных

колебаний молекул является эффективным и более точным, так как погрешность вычислений имеет порядок о (е"+'). Этот метод хорошо алгоритмизируется.

Второй метод решения прямой спектроскопической задачи возмущенных колебаний молекул основывается на теории проекторов применительно к матрице динамических коэффициентов Б(е) . На основании теоретических утверждений этого метода существует спектральное разложение для невозмущенной матрицы Д и оно имеет вид: Б = А]Р]+Л2Р2+... + Аъ„_6Р}„_6, здесь Рк (к=1,2,...ЗЫ-б) - проекторы матрицы £>. В рамках этого метода с использованием резольвенты 2 5 определенной через проекторы, получены выражения в разложении собственных векторов матрицы Б(е), координаты которых равны возмущенным элементам формы к-то колебания и 5-го порядка возмущений

5 ~ 1 Я

ь ь = г (£ як„ь ь_я - £ д„ь ь_„), (4)

п = 1 п = 1

вместе с соответствующими возмущениями частот

®ь=(Чо1мь-») ,(5 = 1,2,3...),(^ = 1,2,...,3^-6). (5)

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

Рк{е) = 1^£пРк с членами возмущений первого, второго и п-го порядив "

ков Рк1, Рк2,...,Рь,, (к=1,2,...,ЗМ-6), которые применяются для исследований особенностей и тонких эффектов в изменениях частот и внутримолекулярных параметров при возмущениях более высоких порядков.

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

0 <" > = (( Ь -1 )' ( И <" ) ь , +

Ь ,<я> = Е (( Д (|> - Л/1)1)Ь + ...

... + (Д("-" - Я1(в",)/)Ь,(1) + Д('0ЬД (7)

(5 = 1,2,...,3^-6), (и = 1,2,...).

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

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

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

Di(£) = (D(£)-AE)P(£)/£ = £ S"Dn + l. (8)

и = 0

Редукционный процесс сначала применяется для получения аналитического представления подпроектора Рл(£) = Pj\ + sPjU) + s2Pf2) +..., возмущенного проектора P¡(s). Доказательство существования разложения в ряд PjX{s) сводится к нахождению явных выражений коэффициентов в этом разложении

р;п) = -PnA2Sf - PjAS + PjASAJ™ -

-SyA2PjX-SA,Pn + S^AMiPn. (9)

Определены выражения для последующих возмущенных подпроекторов.

На основании процесса редукции получены дополнительные члены в разложение частот

юу\е)=а>™ + е<о}х1)+е2а>™+..........(10)

Они получились соответственно равными

íy<I2»=(ír[Z)2JP.1-Z)i5Z)1P.1])1/2,

«,<»> = (tr[D,Pn -D[SD1PJ] - D.SD.Pj, + + DXSDXS'D,^., - ÁlJ)D1S2DlPJ1 -- D2sfD2PjX + DlSDiSf)D2PJl + (11)

+ D^fD.SD.Pj, - D.SD.SfD.SD^})1'2 ,(j = 1,2,...,«).

Формулы (10) и (11) определяют дополнительные члены при возмущениях второго и третьего порядков, которые вместе с <у'2> и <у'3' позволяют учесть изменения частот для возмущенных колебаний молекул и записать окончательное выражение

a>j(e) = coj + so)f + £2{О){Р + <у'12))+ (12)

+ 8г{(о]Ъ) + <3)) +..., (j=l,2,...,3iV-6).

Получены формулы для дополнительных членов в изменения элементов матрицы формы L(s) возмущенных колебаний молекул. Элементы столбцов этой матрицы являются координатами собственных векторов

L^f) {к = 1,2,...,37V - 6 ) матрицы D(e). В рамках редукционного процесса возмущенной частоте co(k\s)t представленной разложением (10), соответствует собственный вектор координатами которого являются элементы формы ¿-го возмущенного колебания. Для него существует аналитическое представление

= + = 1.....ЗЛГ - 6). (13)

Выражения для вкладов второго L ,1 и третьего L 1 порядков получаются равными

= -(^о +DlLkl)), (14)

Ц-f = Z|1)(41)L*2 + 42)Lu - (ALi3 + D2Ltl +D3Lk0)). (15)

Представление к -го собственного вектора матрицы D{s) с учетом возмущенных членов, полученных в рамках процесса редукции, определяется сходящимся рядом

Lt(e) = Lk0+eLtl + e2(Lk2 + L™) + £3(hk3 +L(^)+..... (16)

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

имеют вид . Аналитическое представление координат

/1 = 0

симметрии определяется рядом

ЗД = 1/(2 e'LM)Q = 2 s"(uLw)q . (17)

/1=0 n=0

Влияние возмущений на матрицу коэффициентов кинематического

взаимодействия G ( г ) = ^ £ " G " исследуется на основании мат/1 = 0

ричного соотношения G {s) = L {s ) L {s ) , из которого вклад в возмущение »-го порядка матрицы кинематического взаимодействия получается равным

G (п) = LL(n) + L0)Lin~]) + ... + + L(n) L . (18)

Формулы, определяющие аналитические изменения декартовых координат колеблющихся атомов, записываются следующим образом

Аа(а) = М~т1а Ш = ЛГ'^^У(Г1)0'^, (19)

и=0

здесь Аа(£) = Дог2(^),...,АоГд,^)) - строка изменений а-й

(а = х, у, г) декартовой координаты каждого атома при возмущенных колебаниях в молекуле. Приведены выражения, определяющие влияние возмущений на параметры центробежного искажения а"" (е), а"р (/;), А™. {е), {е), {а,р,у=кул = ~ 6)-

Для и-го возмущенного члена в аналитическом представлении коэффициентов кориолисова взаимодействия имеет место равенство

^ = £ [и . (20)

г = о

Колебательный момент импульса я(£) с использованием определяющей формулы преобразуется к аналитическому представлению

= (21)

5 $' /7 = 0

Аналитическое возмущение обратного для модифицированного тензора инерции /'(<£■) будет проявляться через возмущение соответствующих величин /л"\е), (£■)• При возмущении /^'(¿О имеет место равенство

Л) = -(/0ГЧ«Ш~\ (22)

в котором о^е)выражается через элементы а"" (¿г), (е). Возмущение

второго порядка (г.*) определяется через аналитические представления

элементов второго порядка.

В третьей главе разработана математическая модель обертонных колебаний молекул. Моделирование колебаний атомов в молекуле с кратными частотами 2vs, Зу„ 4у„ 5у„... (.?=/, 2.....ЗИ-6) основывается на теории

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

них <р{А\В)= 2 С у А' , где С,у-численные коэффициенты. Эта связь

'>.7=0

существует на основании теоретического утверждения о том, что если

собственные значения матрицы А и Хх,...,Л„- собственные значения матрицы В, то собственными значениями матрицы (р(Л\В) будут тп следующих чисел q>(jur\Xs\(r = \,...,m\s=\,...,ri).

Осуществлен вывод уравнения обертонных колебаний на основании модификации основного уравнения нормальных колебаний молекул. Для преобразований выбирается вторая матрица Р размерности (т ут), все элементы которой равны единице. Результатом прямого произведения (POD) будет блочная матрица, каждый из т" блоков равен D. Характеристическое уравнение матрицы Р записывается следующим образом: цт-тцт'х =0. Только один корень его имеет ненулевое значение /ui=m. На основании теоретических утверждений и свойств прямого произведения получено уравнение колебаний с кратными частотами

(Р ® £>)(Р, ® Z,) = ягЛ(Р, ® L), (23)

в котором т = 1с , к- порядок обертона. Собственными значениями матрицы (P®D) уравнения (23) являются числа тЛ1,тЛ2,...,тЛ3!/_6, которые пропорциональны квадратам частот kcos (s=],2,...,3N-6), а собственные векторы равны L, / Vw, L2 /4т, ...,L3jV_61yfm . Их координаты являются элементами столбцов матрицы формы обертонных колебаний молекул. Результаты решения уравнения (23) позволяют утверждать, что колебания с кратными частотами kws относятся к тем же неприводимым группам симметрии, к которым относятся нормальные колебания с частотами cos (s=l,2,...,3N-6).

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

d^a/dQks JNf\d^JdS)\dS)ldQks)=\lkN±('Lid^aldSkj, (24) м м

в которой L{ - элементы матрицы формы колебаний, к -порядок обертона; jua(a = x,y,z)- координаты молекулярного дипольного момента. Производные df ia / dQk = Мь определяются по формуле (24) для обертонов любого высокого порядка.

В рамках построенной математической модели выполнен прямой расчет электрооптических параметров обертонных колебаний до шестого порядка молекулы Н21бО и двух изотопов Н2П0 и Н7180. Выбор молекулы воды объясняется тем, что в теоретических исследованиях она является объектом, для которого существует большое количество экспериментальных данных по интегральным интенсивностям, но значительно меньше для изотопов. Результаты численного эксперимента для обертонов 2vs,3vs,4vs ,5vs,6vs молекулы H20 и изотопов Н2'70 и H2lsO в сравнении с существующими экспериментальными данными для первых обертонов

2у,,2у2, 2у, и высоких обертонов 3v2,4v2 Н20 и обертона 2v2 изотопа Н2180 подтверждаются вполне хорошим совпадением с экспериментом.

Вычислены электрооптические параметры СН4 для колебаний с частотами 2^,3^,4^,5^,(^ = 3,4). Сравнение с имеющимися экспериментальными данными для первых обертонов 2vs (.? = 3,4) показывает достаточно хорошее совпадение, а для более высоких обертонов к>3 такой эксперимент отсутствует, и полученные результаты являются предсказательными. Осуществлен численный эксперимент для электрооптических параметров и интегральных интенсивностей в ИК спектрах обертонов до шестого порядка четырехфтористого углерода 12CF4 и изотопа 13CF4 . Несмотря на высокую симметрию этих соединений, обертонный спектр молекулы 12CF4 исследован недостаточно, а для изотопа l3CF4 эксперимент отсутствует. Численный эксперимент осуществлен с использованием комплекса программ «CNDO - SPLINE». Полученные результаты ЭОП обертонных колебаний этих молекул, особенно для высоких обертонов, обосновывают целесообразность проведенных теоретических исследований, надежность алгоритма и эффективность расчетной программы.

Вторая разработанная математическая модель этой главы описывает возбужденные колебания молекул с составными и разностными частотами. Основу ее составляет вывод уравнения колебаний молекул с частотами (соа±соь), (a,b = l,...ßN-6;a*b), которое получается в результате преобразований двух уравнений нормальных колебаний с использованием утверждений и свойств прямого произведения матриц

(Dk ®/„ -/„, ®Dsf\Lk ®Lsf] =(Ак ®/„ -/„, ®Л,)[2](4 ®1,)[2]. (25) Диагональная матрица собственных значений среди своих элементов имеет произведения, пропорциональные квадратам составных и разност-

'У "У 'У "У 'У 'У

ных частот 4л с (cokl+cosi) Атг с (aki-cosj) . Из уравнения (25) элементы тензора формы колебаний определяются в результате прямого произведения Lks = (Lk ® Ls). Для составной частоты (o)kj + a>SJ) соответствующий столбец элементов формы колебаний запишется в следующем виде zj = LtML, Lik,2ls ....Lkk,XII •

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

Gh = (Lk ®Ls)(Lk ®LJ = (Lk ®Ls)(L'k ®V) = (LkL'k)®{LsLs') = Gk ®GS, (26) в которой Gkn Gs — блоки, стоящие на диагонали матрицы G нормальных колебаний молекул.

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

Lh(е) = Lk(е) ® Ls(е) = {Lk®Ls) + s{Lk ® 4> + 4> ® L,) +... . (27)

Экспериментальными методами двумерной лазерной спектроскопии К.Окшшга, О.М.Ропаэ. У.Таттига установлено смещение мод гармонически связанных осцилляторов и рассмотрено его влияние на механическую ангармоничность в ИК спектрах. В связи с этим актуальными являются исследования влияний взаимодействий нормальных колебаний с комбинационными частотами на внутримолекулярные и электрооптические параметры. На основании аналитических представлений матриц Ок(£),От(е),1т(е),1к(£) и

элементов Лк (е), Лт (а), (к,т = 1,2,3,...,ЗЛГ — 6) получено уравнение возмущенных колебаний с комбинационными частотами, при выводе которого учитываются взаимодействия, определяемые взаимным возмущением одного колебания другим по схеме аналитического метода теории возмущений с использованием проекторов.

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

Ч-* =(Ь,®Ьт/)(1+ I ^4). (28)

ГФ!

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

а т к ц / ( ( %

1=\г=\Ойкг Ойт1

через возмущенные элементы тензора формы колебаний

= (^к(г,0 + 4'о-Л ) И = + )'

Осуществлены расчеты производной второго порядка 2-й координаты дипольного момента молекулы метана ,2СН4 и изотопа 13СН4. Для |2СН4 в гармоническом приближении получили величину ¡л:з4 = - 0,040070, а экспериментальное значение М. Лоете оказалось равным (¡/¡дж,?.- - 0,040460. Производная для изотопа 13СН4 по расчетам - 0,038350, но в литературе для нее отсутствует сравнительный эксперимент. Применяя формулы (28) и (29) с учетом взаимного влияния колебательных мод, ангармонические электрооптические параметры получились равными следующим значениям: для молекулы метана /О, = - 0,040870 и для изотопа /Г3-1 = - 0,038740. Анализ результатов показывает, что вклады в значения вычисленных ангармонических параметров с учетом взаимного влияния колебательных мод сравнимы по величине с изменением производной Д//», определяемым эффектом изотопического замещения. Поэтому их необходимо учитывать при прямом расчете ангармонических электрооптических параметров.

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

SU = I к* ■ W - Ч* (W)+ fojt/U - V,uYILJ I (3°)

где lfk¡ - элементы матрицы lk; (l'smj)y- элементы возмущенной матрицы изменения координат атомов l'm,(a,ß,y = x,y,z,).

Обобщением выполненных исследований этой главы является математическое моделирование высоких возбужденных колебаний молекул с комбинационными частотами (kсоа ± 1 соь) при условии (к+1)> 3, где к, 1 -порядки обертонов, а,Ь - номера колебаний. На основании разработанных математических моделей высоких обертонов и колебаний с комбинационными частотами второго порядка, а также операции прямого произведения и ее особых свойств, получено уравнение высоких возбужденных колебаний молекул с частотами (kсоа ± \а>ь). Записаны определяющие выражения для внутримолекулярных параметров, элементов матрицы коэффициентов кинематического взаимодействия и расчетные формулы для колебательно-вращательных и ангармонических электрооптических параметров высоких возбужденных колебаний с комбинационными частотами.

В четвертой главе выполнены исследования возбужденных колебательных состояний многоатомных молекул методами операторной теории возмущений. В квантовой теории колебаний молекул теоретическое описание возбужденных колебательных состояний сводится к построению математической модели, которая основывается на определении корректного колебательного гамильтониана. Таким гамильтонианом в теории колебаний молекул с учетом колебательно-вращательных взаимодействий является гамильтониан, который предложили E.B.Wilson и J.B.Hovard, а записал в преобразованной форме J.K.G.Watson. При решении задач колебательной спектроскопии с таким гамильтонианом используются методы операторной теории возмущений. Самым распространенным среди них в последнее время является метод контактных преобразований (КП), впервые предложенный Ван Флеком. Теорию метода контактных преобразований в приложении к задачам колебательно-вращательной спектроскопии молекул разработали С. Amat, H.H. Nielsen, G. Tarrago. В дальнейшем метод КП успешно развивали и использовали в теории колебательных и вращательных спектров молекул М.Р.Алиев, J.K.G.Watson, В.М.Михайлов, Ю.С.Макушкин, В.И.Перевалов, В. И.Стариков, В. Г.Тютерев.

Метод контактных преобразований является основным при построении теоретической модели эффективного дипольного момента молекул сферической симметрии Td на примере молекулы метана СН4 для возбужденных колебательных состояний и переходов между ними с использованием квантово-химической модели электронного строения молекул и ее электрического дипольного момента в рамках полуэмпирического метода CNDO/2 и метода аппроксимации и численного дифференцирования сплайн-функциями. В колебательной спектроскопии молекула метана СН4 занимает особое место, что подтверждается большим числом публика-

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

Осуществлены прямые расчеты параметров операторов однокванто-вых колебательно-вращательных переходов молекулы метана: соответственно 0"^- параметров оператора эффективного дипольного М,2 и ©"^ - коэффициентов оператора М]3. Для них были вычислены электрооптические параметры (и = 3,4; а,/3,у - х, у, г). Они применялись при проверке равенства:

©з>/ = -@4>4 , (31)

0.2317 = 0.2352 .

Численные значения правых и левых частей в (31) позволяют сделать вывод о том, что для прямых расчетов электрооптических параметров и

® и (и — 3,4) коэффициентов оператора М12 результаты показывают достаточно хорошее совпадение. В работах 1.0г1ег, К.Рох были введены также параметры <?3 и С4, которые соответственно равны с„ = 0;/т /(и ~ 3,4). Вычисленные значения для них и сравнение с экспериментом показывают (теор.)=0.00364, с3 (эксп.)=0.00413,

<^4 (теор.)=0.00474, (эксп.)=0.00492, ^

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

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

Вычислены коэффициенты 0? и ©J , &к в операторах центробежного типа (М03 +М05), ангармонического М21 и кориолисова М22типов молекулы метана. Однако экспериментальные значения существуют пока только для параметров 0„„ (и =3,4), 0? . Они приведены в работах 1.0г!ег и А.С.Ьигйг. Сравнительный анализ теоретических значений 0 -параметров с имеющимся экспериментом показывает хорошее согласие:

©зз (теор.)= 0.018211), ©33(эксп.)= 0.02041); (33)

©: (теор.)= 1.97679-10"5/) , ©У(эксп.)= 2.410-10"51>;

©44(теор.)= 0.01641 В , ©44(эксп.)= 0.013821).

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

Измерения штарковских сдвигов, выполненные Н.8азас1а, К.Зигитига, С. ^ЫЬазЫ в полосе 2\>3 молекулы СН4, позволили из эксперимента определить -параметр оператора эффективного дипольного момента ко-риолисового типа М22, который получился равный (-2,2-10 I)). Вычисленное значение этого параметра ©£ = (-2,8096-10~4£>) и вполне хорошо согласуется с экспериментальным. Следует отметить, что значение параметра 0^ было предсказано В.М. Михайловым и автором работы на девять лет раньше, чем осуществлен спектроскопический эксперимент. Результаты численного эксперимента для 0 -параметров операторов эффективных дипольных моментов одно- и двухквантовых колебательных переходов в основном и возбужденных колебательных состояниях определяют вклады в интенсивности ИК полос, заимствованные из вращательных спектров, за счет учета колебательно-вращательных взаимодействий.

На основании метода контактных преобразований во втором порядке теории возмущений получены выражения матричных элементов эффективного дипольного момента (оаиь |м|и0 + 2,иь ± 1) и вероятностей трехбайтовых колебательных переходов, расчетные формулы для абсолютных интенсивностей спектральных линий отдельных переходов и интегральных интенсивностей ИК полос поглощения многоатомных молекул с комбинационными частотами (2уа±уь),в которых учитывается влияние вкладов ангармонических механических и электрооптических параметров до четвертого порядка, а также температурной зависимости на их формирование.

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

функции ^ (аз) и обратной для нее (ге), с помощью которых осуществляется преобразование подобия волновых функций невозмущенного колебательного гамильтониана к=1,2,...,ЗИ-6} в искомые волновые функции колебательно-вращательного гамильтониана {хУк(х),к=1,2,...,ЗМ-6} и вывод определяющих выражений для возмущений волновых функций и уровней энергии возбужденных колебательных состояний многоатомных молекул.

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

хрО) = = ^'Ч',, = = (34)

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

Ч/О) = = > (35)

-5,2Н(1)РН(1)Р-^РН("52Н("Р)4/А И Х.Д. (36)

В формулах (35) и (36) используются обозначения: Р - тотальный проектор, 5 - приведенная резольвента невозмущенного гамильтониана.

Энергия возмущенных колебательных состояний определяется формулой

Е (ге)= и-2у ( ас)Н(эг) К (х), (37)

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

Е1 = 1Г(РН™Р + РН20^(1) + 20Р). (38)

Возмущение энергии второго порядка определяется выражением Ег =^(РН<2)Р + 2<2)Н20Р + РН20Г„(2) +

+ +2(1>Н20Ж„(1) + РН(,)^(1)) и т.д. (39)

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

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

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

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

Затем для каждого , принадлежащего сетке А ^ ( вычисляется значение координаты дипольного момента по формуле метода €N00/2 и получается следующий набор значений Р"(а = х,у,г) относительно положения равновесия

Построение интерполяционного кубического сплайна Р" (А, удовлетворяющего условию

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

Для вычисления смешанных производных второго порядка диполь-ный момент молекулы Р необходимо рассматривать как функцию двух переменных Р = Р(5Я,, на прямоугольной сетке значений координат симметрии и ^ . Дальнейшие исследования основываются на определении и свойствах двумерных кубических сплайн - функций. Способ построения двумерного сплайна для функции дипольного момента Р = РО?^,состоит в последовательном определении кубических сплайнов от одной переменной Яр- при фиксированном значении другой переменной . Для этого

задается изменение первой координаты симметрии с шагом А,, а значения второй координаты симметрии изменяются с шагом к2.

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

ра ра ра ра ра ра р

■ ■ ■ г-3 > -2 > г-\ 1 г0 1 г\ ■> 2 > ' 3

3 >••

(41)

а

(42)

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

Для учета вкладов в изменения ангармонических электрооптических параметров применяются определяющие выражения для элементов тензора формы колебаний, полученные на основании аналитических методов матричной теории возмущений. Разработанная методика аппроксимации функции дипольного момента в возбужденных колебательных состояниях многоатомных молекул реализована с использованием обобщенного вычислительного алгоритма, составленного на основании результатов, полученных методами матричной теории возмущений и комплекса расчетных программ «СНОО-БРЬШЕ».

Основные результаты и выводы

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

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

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

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

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

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

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

6. На основании построенной модели обертонных колебаний осуществлен вычислительный эксперимент электрооптических параметров молекулы Н2О и изотопов Н? 170 и Н2 180 для обертонов до шестого порядка. Выполнен прямой расчет электрооптических параметров обертонных колебаний молекулы СН4 до пятого порядка. Реализован численный эксперимент для электрооптических параметров и интегральных интенсивностей в ИК спектрах поглощения обертонов до шестого порядка четырехфто-ристого углерода 12CF4 и изотопа 13CF4, обертонный спектр молекулы 12CF4 исследован недостаточно, а для изотопа 13CF4 он отсутствует. Сравнение вычисленных электрооптических параметров молекулы 1ЬО и ее изотопов Н2 170, Н2 180, молекулы СН4, четырехфтористого углерода ~CF4 и изотопа 13CF4 с имеющимися экспериментальными данными подтверждает достоверность разработанной математической модели, надежность вычислительного алгоритма и эффективность комплекса расчетных программ «CNDO-SPLINE». Однако большая часть результатов численного эксперимента является предсказательной, особенно для высоких обертонных колебаний.

7. Разработана математическая модель возбужденных колебаний молекул с составными частотами (соа ± ojb) основанная на выводе уравнения, полученного в результате преобразования фундаментального уравнения нормальных колебаний с использованием операции прямого произведения матриц и его важных свойств. Исследовано влияние межмодового взаимодействия при колебаниях с составными частотами на колебательно-вращательные и ангармонические электрооптические параметры. Получено уравнение высоких возбужденных колебаний молекул с комбинационными частотами (k соа ± 1«/,).

8. На основании метода контактных преобразований операторной теории возмущений осуществлены расчеты вкладов в матричные элементы операторов эффективного дипольного момента для одно- и двухбайтовых колебательных переходов с учетом колебательно-вращательных взаимодействий молекул сферической симметрии типа метана СН4. Выполнены прямые расчеты параметров для операторов эффективного дипольного момента М12 и М13 одноквантовых колебательно-вращательных переходов. Вычислены 0? и 0А- -параметры соответствующих центробежных дипольных моментов М03 и М05 в основном и возбужденных колебательных состояниях, активизирующих запрещенные вращательные переходы. Определены коэффициенты ангармонического оператора дипольного момента М21. Осуществлен прямой расчет параметров для оператора дипольного момента кориолисова типа М22, для которого значение параметра было предсказано на девять лет раньше, чем осуществлен спектроскопический эксперимент. Достаточно хорошее согласие вычисленных значений 0 - параметров с имеющимися экспериментальными данными подтверждает достоверность разработанной теоретической модели функции дипольного момента для молекул сферической симметрии и определения дополнительных вкладов в интенсивности колебательных переходов, заимствованных из вращательных спектров.

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

и обратная для нее ^(эг), с помощью которых осуществляются преобразования подобия волновых функций невозмущенного колебательного гамильтониана {, к = 1,2,...,37^ - 6} в волновые функции возмущенного гамильтониана {Ч'Даг}, к = 1,2,...,ЗЛГ -6}. Получены явные выражения для вкладов в возмущенные волновые функции колебательно-вращательного гамильтониана, преобразованного по схеме Алиева-Уотсона, до четвертого порядка теории возмущений. Определены соответствующие возмущения энергии возбужденных колебательных состояний многоатомных молекул.

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

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

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

12. Составлен алгоритм и создана расчетная программа «ОВЕЯТОЫ», предназначенная для реализации численного эксперимента интенсивностей многофотонных колебательных переходов из основного и возмущенных колебательных состояний с учетом влияния механической и электрооптической ангармоничности, колебательно-вращательного взаимодействия и анализа влияния температуры на их формирование.

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

Основные публикации по теме диссертации Статьи в рецензируемых журналах

1. Гавва С.П. Аналитические возмущения волновых функций и уровней энергии колебательных состояний молекул // Оптика и спектроскопия. 2009. Т. 104. №4. С. 533 - 536.

2. Гавва С.П. Влияние взаимодействий возбужденных колебаний на молекулярные параметры // Известия вузов. Физика. 2008. № 3. С. 27 - 32.

3. Гавва С.П. Влияние взаимодействия мод на спектроскопические параметры молекулярных колебаний с составными частотами // Журнал прикладной спектроскопии. 2006. Т. 73. № 4. С. 421 - 425.

4. Гавва С. П. Моделирование высоких обертонных колебаний молекул // Известия вузов. Физика. 2005. №3. С. 48 - 52.

5. Гавва С. П. Высокие возбужденные колебания молекул с составными и разностными частотами // Оптика и спектроскопия. 2005. Т. 98. №6. С. 906-909.

6. Гавва С.П. Возбужденные колебания молекул с составными и разностными частотами // Оптика и спектроскопия. 2003. Т. 95. №2. С. 225-229.

7. Гавва С.П. Спектроскопические параметры высоких обертонных колебаний молекул // Журнал прикладной спектроскопии. 2003. Т.70. №5. С. 598 - 602.

8. Гавва С.П. Обертонные колебания молекул // Оптика и спектроскопия. 2002. Т. 93. № 5. С.735 - 738.

9. Гавва С.П. Коэффициенты кинематического взаимодействия возбужденных колебательных состояний // Журнал прикладной спектроскопии. 2002. Т. 69. № 1. С.120- 123.

10. Гавва С.П. Ангармонические молекулярные параметры возбужденных колебательных состояний // Известия вузов. Физика. 2002. №5. С.45 - 47.

11. Гавва С.П. Алгебраический метод исследований возбужденных колебательных переходов// Известия вузов. Физика. 2001. №6. С. 58 - 61.

12. Гавва С. П. Алгебраический подход к исследованию возбужденных колебательных состояний // Оптика и спектроскопия. 2001. Т. 91. №4. С. 611-613.

13. Гавва С.П., Элькин М.Д., Свердлов Л.М. О расчетных формулах для интенсивностей составных частот и обертонов в спектрах KP и ИК // Оптика и спектроскопия. 1981. Т.51. № 4. С. 728-731.

14. Гавва С. П., Болотина Э.Н, Свердлов Л. М., Элькин М.Д. Электрооптические параметры и абсолютные интенсивности ИК полос молекулы воды и ее изотопозамещенных // Оптика и спектроскопия . 1982. Т. 53. № 1. С. 46-50.

15. Гавва С.П., Свердлов Л.М. Формулы для интенсивностей полос 2va ±ve в ИК спектрах многоатомных молекул // Оптика и спектроскопия. 1984. Т.57. № 4. С. 623-627.

Публикации в других изданиях

16. Гавва С.П. Аналитические возмущения колебательных состояний многоатомных молекул // Проблемы оптической физики и биофотоники. Саратов: изд-во «Новый ветер», 2010. С. 261-267.

17. Гавва С.П., Сальников А.Н. Аналитические представления частот и элементов тензора формы возбужденных колебаний молекул // Проблемы оптической физики и биофотоники. Саратов: изд-во «Новый ветер», 2009. С.132- 138.

18. Гавва С.П. Математическое моделирование взаимодействий колебательных мод возбужденных состояний молекул // Проблемы оптической физики и биофотоники. Саратов: изд-во «Новый ветер», 2008. С. 151-156.

19. Гавва С.П. Особенности кинематического взаимодействия при колебаниях молекул с составными частотами // Проблемы оптической физики. Саратов: изд-во «Новый ветер», 2006. С. 205 - 209.

20. Гавва С.П. Математическое моделирование высоких обертон-ных колебаний молекул // Оптические проблемы физики. Саратов: изд-во СГУ, 2005. 4.1. С. 22-28.

21. Гавва С.П. Аналитический подход в исследовании молекулярных параметров возбужденных колебательных состояний // Оптические проблемы физики: кн.1. Саратов: изд-во СГУ, 2004. С.15-21.

22. Gawa S.P. Analitical theory of the exited vibrational states of molecules // Proceedings of the Society of Photo-optical instrumentation engineers (SPIE), 2004. V. 5476. P. 160 - 164.

23. Gawa S.P. Projectors theoiy usu for research into molecular parameters of exited vibrational states // Proceedings of the Society of Photo-optical instrumentation engineers (SPIE), 2003. V. 5311. P. 121 - 125.

24. Gawa S.P. Spectroscopic parameters of high exited vibrational states of molecules // Proceedings of the Society of Photo-optical instrumentation engineers (SPIE), 2002. V. 5067. P. 265 - 270.

25. Гавва С.П. Математическое моделирование возмущенных колебаний молекул // Проблемы оптической физики. Саратов: изд-во СГУ, 2002. С. 53 - 58.

26. Гавва С.П. Тензор коэффициентов формы колебаний молекул сферической симметрии // Актуальные вопросы научных исследований. Спектроскопия молекул: межвуз. сб. Саратов: изд-во СГУ, 1997. 4.1. С. 16 - 20.

27. Гавва С.П. Параметры кинематической ангармоничности // Актуальные вопросы научных исследований: Спектроскопия молекул: межвуз. сб. науч. тр. Саратов: изд-во СГУ, 1998. 4.2. С. 40 - 47.

28. Гавва С.П. Элементы высшего порядка тензоров формы валентных колебаний молекул // Актуальные вопросы научных исследований. Теоретическая и прикладная спектроскопия: межвуз. сб. науч. тр. Саратов: изд-во СГУ, 1998. 4.2. С. 56 - 59.

29. Гавва С.П., Землянухина И.Я., Болотина Э. Н. Сплайн аппроксимация функции дипольного момента. Вращательные спектры молекул // Сборник ИС АН СССР. М.: Наука, 1986. С. 188 - 211.

30. Гавва С.П., Михайлов В.М. Теоретическое определение функции дипольного момента и дипольных моментов колебательно-вращательных переходов молекулы метана. Спектроскопия высокого разрешения малых молекул // Сборник ИС АН СССР. М: Наука, 1988. С. 148 -195.

31. Гавва С.П. Элементы формы высших порядков деформационных колебаний молекул // Спектроскопия и физика молекул: межвуз. сб. науч. тр. Саратов: изд-во СГУ, 1998. С. 9 - 13.

32. Гавва С.П. Тензоры высшего порядка формы колебаний молекул / Сарат. гос. техн. ун-т. Саратов, 1998. 9 с. Деп. в ВИНИТИ.

33. Гавва С.П. Тензоры высокого порядка формы деформационных колебаний молекул // Актуальные вопросы научных исследований. Теоретическая и прикладная спектроскопия: межвуз. сб. науч. тр. Саратов: изд-во СГУ, 1999. Ч. 2. С. 23 - 29.

34. Gawa S.P. Elements of tensors of molecules exited vibration form // SPIE Proceedings. 1999. V. 4063. P.251- 254.

35. Гавва С.П., Болотина Э.Н., Землянухина И.Я.Теоретические исследования влияний температуры на формирование полос ИК спектров молекулы S02 // Известия вузов. Физика. 1989. 7 с. Деп. в ВИНИТИ, №7398-В89.

36. Gawa S.P., Mikhailov V.M. Unharmonic electrooptical parameters and dipole moment for methane // SPIE Proceedings. Omsk, 1991. V.1811. P. 222 - 225.

37. Гавва С.П. Параметры дипольных моментов возбужденных колебательных переходов // Актуальные вопросы научных исследований. Теоретическая и прикладная спектроскопия: межвуз. сб. науч. тр. Саратов: изд-во СГУ, 1999.4.2. С.35 - 40.

38. Гавва С.П. Моделирование возбужденных колебательных переходов // Спектроскопия и физика молекул: тр. III межвуз. конф. Саратов: изд-во СГУ,1998. С.15-19.

39. Гавва С.П. Центробежный дипольный момент молекулы метана//Проблемы оптической физики. Саратов: изд-во СГУ, 2000. С. 106-109.

40. Гавва С.П. Ангармонические коэффициенты оператора центробежного дипольного момента // Проблемы оптической физики. Саратов: изд-во СГУ, 2001. С.73 - 77.

41. Gawa S.P. Higher order tensors of molecules vibration form//Book of Abstracts XXIV European congress on molecular spectroscopy. Prague. 1998. P. 132.

42. Gawa S.P. Quantum chemical study of molecular spectroscopic parameters //Book of Abstracts XXIV European congress on molecular spectroscopy. Prague. 1998. P. 133.

43. Гавва С.П., Болотина Э.Н., Землянухина И.Я. Формирование ин-тенсивностей ИК полос 2v, 3v частот молекулы СО в интервале температур 300-3000 К // Известия вузов. Физика. 1987. 8 с. Деп. в ВИНИТИ, №7299-В87.

44. Гавва С.П. Алгебраическая теория возбужденных колебательных состояний молекул // Тезисы докладов XXII Съезда по спектроскопии. Звенигород. 2001. С.45.

45. Гавва С.П.Ангармонические параметры операторов колебательно-вращательного взаимодействия // Тезисы докладов XXII Съезда по спектроскопии. Звенигород. 2001. С.46.

46. Гавва С.П. Влияние возмущений на спектроскопические параметры молекул // Тезисы докладов XXIII Съезда по спектроскопии. Звенигород. 2005. С.159.

47. Гавва С.П. Высокие возбужденные колебательные состояния молекул // Тезисы докладов XXIII Съезда по спектроскопии. Звенигород. 2005. С.160.

48. Гавва С.П. Аналитические возмущения волновых функций и уровней энергии колебательных состояний молекул // Тезисы докладов XXIV Съезда по спектроскопии. Москва, Троицк. 2010. С.165 - 166.

49. Гавва С.П., Сальников А.Н. Аналитические модели в колебательной спектроскопии многоатомных молекул // Вестник СГТУ. Саратов: СГТУ, 2010. №1 (44). С. 7 - 17.

50. Гавва С.П., Сальников А.Н. Квантово-механическая аналитическая модель возбужденных колебательных состояний многоатомных молекул // Вестник СГТУ. Саратов: СГТУ, 2010. №2(45). С. 7-19.

51. Гавва С.П. Математическое моделирование возмущенных колебательных состояний многоатомных молекул // Физико-математическое моделирование систем. Воронеж: изд-во ВГТУ, 2011. 4.2. С. 68 - 74.

Авторское свидетельство

Гавва С.П. Программа «ОВЕЯТОЫ» для расчета абсолютных интен-сивностей спектральных линий и интегральных интенсивностей обертонов и составных частот многоатомных молекул в зависимости от температуры. Номер авторского свидетельства № 2011613491. Зарегистрировано в Реестре программ для ЭВМ 04. 05. 2011.

Полппсаио в печать 18.08.11 г. Формат 60x84/16 Усл.печ. л. 1.2 Тираж 120 -жз. Заказ № 167-Т. Отпечатано с оригинал-макета в типографии СГУ 410012, г. Саратов, ул Б. Казачья 112 Л, 8-й корпус, тел. 27-33-85

 
Текст научной работы диссертации и автореферата по физике, доктора физико-математических наук, Гавва, Светлана Павловна, Саратов

ГОУ ВПО «Саратовский государственный технический университет»

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

052011 52225

Гавва Светлана Павловна

ВОЗМУЩЕННЫЕ КОЛЕБАТЕЛЬНЫЕ СОСТОЯНИЯ МНОГОАТОМНЫХ МОЛЕКУЛ

01.04.05 - Оптика

Диссертация на соискание учёной степени доктора физико-математических наук

Научный консультант: доктор физико-математических наук, профессор Березин Кирилл Валентинович

Саратов—2011

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ..................................................................................9

Глава 1

АНАЛИТИЧЕСКИЕ МЕТОДЫ В ТЕОРИИ ВОЗМУЩЕННЫХ КОЛЕБАНИЙ МНОГОАТОМНЫХ МОЛЕКУЛ

Введение.....................................................................................54

1.1. Основные положения алгебраической теории возмущений.................57

1.2. Алгебраический метод определения изменений частот

и элементов формы колебаний молекул.........................................61

1.3. Представление обратной матрицы элементов формы колебаний..........66

1.4. Аналитические возмущения частот и элементов тензора

формы колебаний молекул.......................................................69

1.5. Методы построения проекторов в матричной теории........................72

1.6. Метод матричной теории возмущений с использованием проекторов..............................................................................74

1.7. Аналитический метод исследования возмущенных колебаний

молекул с использованием проекторов........................................76

1.8. Аналитические возмущения проекторов матрицы колебательного

уравнения..............................................................................80

1.9. Компонентные матрицы и их связь с резольвентой

в матричной теории возмущений.................................................83

1.10. Аналитические возмущения компонентных матриц

и резольвенты.........................................................................86

1.11. Аналитический метод исследования возмущенных колебаний

молекул через компонентные матрицы..............................................88

1.12. Сравнительный анализ результатов применения разных аналитических методов матричной теории возмущений....................92

1.13. Алгоритм определения аналитических возмущений

собственных значений и собственных векторов..............................98

Глава 2

АНАЛИТИЧЕСКИЕ ПРЕДСТАВЛЕНИЯ ВНУТРИМОЛЕКУЛЯРНЫХ

И СПЕКТРОСКОПИЧЕСКИХ ПАРАМЕТРОВ ВОЗМУЩЕННЫХ КОЛЕБАНИЙ

Введение....................................................................................102

2.1. Разделение нормальных колебаний и вращательного движения молекул................................................................................103

2.2. Уравнение колебаний молекул с учетом колебательно-вращательных взаимодействий......................................................................107

2.3. Редукционный процесс в матричной аналитической теории возмущений...........................................................................114

2.4. Возмущения высших порядков частот и элементов

формы колебаний молекул, полученные с применением процесса редукции..............................................................................120

2.5. Аналитические возмущения структурных

параметров многоатомных молекул............................................124

2.6. Аналитические возмущения матрицы кинематического взаимодействия......................................................................132

2.7. Аналитические возмущения параметров колебательно-вращательных взаимодействий.....................................................................................137

Глава 3

КОЛЕБАНИЯ МОЛЕКУЛ С КРАТНЫМИ И КОМБИНАЦИОННЫМИ ЧАСТОТАМИ

Введение....................................................................................143

3.1. Колебания многоатомных молекул с кратными частотами................147

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

и собственных вектров..............................................................148

3.1.2. Уравнение колебаний молекул с кратными частотами............152

3.2. Структурно-динамические и кинематические параметры

колебаний молекул с кратными частотами....................................154

3.3. Параметры колебательно-вращательного взаимодействия обертонных колебаний молекул.................................................155

3.4. Электрооптические параметры высоких обертонных

колебаний молекул..................................................................159

3.5. Алгоритм численного эксперимента для спектроскопических параметров обертонных колебаний молекул...................................163

3.6. Возмущенные колебания молекул с составными и

разностными частотами............................................................166

3.6.1. Уравнение колебаний многоатомных молекул с

составными и разностными частотами........................................167

3.7. Ангармонические электрооптические параметры колебаний

молекул с комбинационными частотами.......................................171

3.8. Высокие возмущенные колебания молекул с составными

и разностными частотами.....................................................178

3.9. Молекулярные параметры высоких возмущенных колебаний

с комбинационными частотами........................................................182

3.10. Влияние взаимодействий колебательных мод на спектроскопические параметры молекул.....................................186

3.10.1. Преобразованное уравнение возмущенных колебаний молекул с составными и разностными частотами....................187

3.10.2. Влияние межмодового взаимодействия на спектроскопические параметры...................................................................189

Глава 4

АНАЛИТИЧЕСКИЕ ОПЕРАТОРНЫЕ МЕТОДЫ ИССЛЕДОВАНИЙ ВОЗМУЩЕННЫХ КОЛЕБАТЕЛЬНЫХ СОСТОЯНИЙ

МНОГОАТОМНЫХ МОЛЕКУЛ

Введение....................................................................................197

4.1. Квантово-механическая теория колебательных

состояний молекул..................................................................199

4.2. Квантово-механический колебательно-вращательный гамильтониан...204

4.3. Редуцированный колебательно-вращательный гамильтониан............206

4.4. Метод контактных преобразований операторной

теории возмущений..................................................................212

4.5. Эффективный дипольный момент многоатомных молекул................217

4.6. Интенсивности полос колебательно-вращательных переходов

в ИК спектрах поглощения многоатомных молекул..........................220

4.7. Одноквантовые и двухквантовые колебательно-вращательные переходы молекул..................................................................224

4.8. Определение ©-параметров операторов дипольных моментов колебательно-вращательных переходов молекулы метана.................227 '

4.8.1. Квантово-химический расчет параметров функции

дипольного момента молекулы метана...................................227

4.8.2. Операторы эффективного дипольного момента одноквантовых колебательных переходов молекулы метана...........................230

4.8.3. Операторы эффективного дипольного момента колебательно-вращательных переходов в основном и возбужденных состояниях молекулы метана..............................................234

4.9. Формулы для интенсивностей полос (2 + уь ) и (2 уа — V/,)

в РЖ спектрах многоатомных молекул.........................................241

4.10. Аналитические возмущения волновых функций и

собственных значений линейных операторов.................................249

4.10.1. Общая характеристика аналитической теории возмущений

линейных операторов.....................................................249

4.10.2. Преобразование подобия собственных векторов возмущенного оператора...........................................................................................252

4.10.3. Алгоритм моделирования аналитических возмущений уровней энергии и волновых функций возбужденных состояний многоатомных молекулярных систем.................260

4.11. Аналитические возмущения волновых функций колебательно-вращательного гамильтониана..................................................262

4.12. Аналитические возмущения уровней энергии колебательных состояний молекул................................................................264

Глава 5

АППРОКСИМАЦИЯ ФУНКЦИИ ДИПОЛЬНОГО МОМЕНТА МОЛЕКУЛ В ВОЗБУЖДЕННЫХ КОЛЕБАТЕЛЬНЫХ СОСТОЯНИЯХ

Введение....................................................................................269

5.1. Полуэмпирические квантово-химические методы расчета молекул. Приближение нулевого дифференциального перекрывания...............271

5.2. Метод полного пренебрежения дифференциальным перекрыванием......................................................................273

5.3. Характеристика квантово-химических моделей молекул..................278

5.4. Методика прямого расчета электрооптических параметров многоатомных молекул...........................................................282

5.5. Формулы связи электрооптических параметров с производными дипольного момента по геометрическим координатам.......................286

5.6. Дипольный момент многоатомных молекул в приближении нулевого дифференциального перекрывания.............................................296

5.7. Аппроксимация функции дипольного момента молекул полиномиальными кубическим одномерными и двумерными сплайнами.............................................................................298

ЗАКЛЮЧЕНИЕ........................................................................310

СПИСОК ЛИТЕРАТУРЫ..........................................................315

Приложение 1

КОМПЛЕКС ПРОГРАММ «CNDO - SPLINE» ДЛЯ ПРЯМОГО РАСЧЕТА ЭЛЕКТРООПТИЧЕСКИХ ПАРАМЕТРОВ МОЛЕКУЛ РАЗЛИЧНЫХ КЛАССОВ.................................................................347

Приложение 2

ПРОГРАММА «OBERTON» ДЛЯ РАСЧЕТА ИНТЕГРАЛЬНЫХ ИНТЕНСИВНОСТЕЙ ОБЕРТОНОВ, СОСТАВНЫХ ЧАСТОТ И АБСОЛЮТНЫХ ИНТЕНСИВНОСТЕЙ ВОЗМУЩЕННЫХ КОЛЕБАТЕЛЬНЫХ ПЕРЕХОДОВ В ЗАВИСИМОСТИ ОТ ТЕМПЕРАТУРЫ..................................................................373

ВВЕДЕНИЕ

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

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

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

Теория колебательных спектров молекул, изложенная в [1-10], основывается на динамической модели молекулы как системы, состоящей из N

числа материальных точек, на которые наложены голономные связи, а действующие между ними силы являются потенциальными. Число N совпадает с количеством атомов в молекуле, совершающих малые гармонические колебания относительно положения равновесия. Процесс колебаний атомов в молекулах моделируется уравнениями Лагранжа второго рода. Для этой системы уравнений М. А. Ельяшевичем [1] и независимо от него Е. Вильсоном [2] впервые было предложено в качестве обобщенных координат использовать естественные колебательные координаты, которые определяются« изменением межатомных расстояний по химическим связям, валентных углов и углов внутреннего вращения. В системе естественных колебательных координат было получено фундаментальное уравнение теории колебаний молекул [1,2]. В рамках этой физической модели М.А. Ельяшевичем [1] и М. В. Волькенштейном [5,6] была разработана валентно-оптическая теория ин-тенсивностей инфракрасных спектров поглощения многоатомных молекул, в развитие которой большой вклад внесли Л. М. Свердлов, М. А. Ков-нер [9] и Л. А. Грибов [7,10].

Все теоретические исследования в колебательной спектроскопии молекул основываются на решении фундаментальных задач, которые разделяются на прямую механическую задачу по вычислению частот и элементов формы колебаний, а также обратные спектральные задачи по определению силовых постоянных на основании экспериментальных частот и расчету электрооптических параметров при использовании измеренных интегральных интен-сивностей инфракрасных (ИК) спектров поглощения молекул [8-20].

Прямая механическая задача определяется решением фундаментального колебательного уравнения, которое относится к типу уравнений по определению собственных значений X 5 и соответствующих собственных векторов (з=1,2,...,ЗЫ-6), а решается оно методами численного спектрального анализа матричной теории [72, ВО, 81, 86, 95]. Это обеспечивает разработку надежных алгоритмов по решению прямой спектроскопической задачи для вычисления частот, элементов формы колебаний и смещений ато-

мов относительно положения равновесия. Такие алгоритмы описаны в [8, 1418, 20, 23] и реализованы в виде программ, с использованием которым выполнены расчеты при решении прямой механической задачи для большого числа органических и неорганических многоатомных молекул Л. М. Свердловым и М. А. Ковнером [9], Л. А. Грибовым и В. А. Дементьевым с сотрудниками [8, 14-18, 20, 23]. Таким образом, теория нормальных колебаний молекул, методы решения прямых и обратных спектральных задач вместе с библиотеками алгоритмов и расчетных программ уже в 70-е годы прошлого века определяли успешное применение математического моделирования как эффективного метода фундаментальных исследований в колебательной спектроскопии многоатомных молекул. Однако все теоретические исследования и расчеты* в колебательной спектроскопии молекул, содержащиеся в [7-20], были выполнены для фундаментальных нормальных колебаний многоатомных молекул.

Современные теоретические и высокочувствительные экспериментальные методы молекулярной спектроскопии являются мощными, эффективными и исключительно разнообразными по возможностям и широкому диапазону их применения. Они обеспечивают получение качественно новой спектроскопической информации для многоатомных молекул. Состояние экспериментальных исследований характеризуется тенденциями быстрого увеличения разрешающей способности и перехода к наблюдению более тонких эффектов и особенностей в инфракрасных (ПК) спектрах поглощения молекул. Применение методов внутрирезонаторной, лазерной и оптико-акустической спектроскопии [39-41] позволяет использовать спектрометры с пороговой чувствительностью коэффициента поглощения 10"8 см"1, а погрешность определения интенсивностей линий составляет 2-3%. Поэтому в последнее время благодаря интенсивному развитию инфракрасной и микроволновой спектроскопии, разработке новых экспериментальных методов стали возможными исследования слабых спектров молекул, определяемых возбужденными и особенно высоковозбужден-

ными переходами, интенсивности которых на несколько порядков меньше интенсивностей фундаментальных нормальных колебаний, а для частот наблюдаются изменения в тысячные доли обратных сантиметров. Результаты таких исследований выполнены Ю.Н. Пономаревыми и J1.H. Синицей с сотрудниками и приведены в [39-40].

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

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

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

щих и позволяющих моделировать процессы, происходящие в молекулах, которые проявляются в изменении частот, структурных, динамических, колебательно-вращательных и электрооптических параметров и интенсивностей в ИК спектрах поглощения при возмущении колебаний атомов в молекулярных системах. Отметим, что в динамической теории колебаний многоатомных молекул [1-3, 9,.10]; такие исслед