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

Источник финансирования НИР

грант РФФИ

Этапы НИР

# Сроки Название
1 1 января 2015 г.-31 декабря 2015 г. Прецизионное моделирование оптических свойств двухатомных газов при высоких температурах
Результаты этапа: (1) Разработаны прецизионные методы решения радиального урвнения Шредингера с учетверенной машинной точностью, позволившие впервые достоверно оценить вероятности очень слабых обертонных переходов молекулы CO вплоть до энергии ее диссоциации. (2) Разработан универсальный алгоритм и осуществлена его програмная реализация для моделирования электронно-колебательно-вращательных энергий, колебательных волновых функций и вероятностей ровибронных переходов двухатомных молекул как в адиабатическом, так и неадиабатическом приближении (3) Предложен траекторный (классический) анализ вращательной динамики многоатомных молекул позволяющий предсказать энергии состояний в условиях сильного вращательного возбуждения (4) Выполнен критический обзор литературы по исследованиям поведения атомов натрия в экзосферах Луны и Меркурия. Проведены термодинамические расчеты химического состава ударно-образованного облака. Получено равновесное содержание NaO в широком диапазоне температуры и давления закалки для нескольких вариантов химического состава ударников (углистые хондриты, кометы) и скоростей столкновения ударников с телами Солнечной системы. Показано, что одним из основных источников атомов натрия в планетных экзосферах во время активности сильных ежегодных метеорных потоков Персеид, Квадрантид и Геминид является фотолиз NaO солнечными фотонами. Оценено время фотолиза оксида натрия солнечными фотонами на основании корреляции между временем фотолиза и энергией диссоциации атмосферных двухтомных молекул (5) Получено новое доказательство того, что в спектрах звезд типа Миры Кита эмиссионные линии водорода и металлов образуются ниже области, в котором формируются абсорбционные молекулярные полосы
2 1 января 2016 г.-31 декабря 2016 г. Прецизионное моделирование оптических свойств двухатомных газов при высоких температурах
Результаты этапа: 1. Путем сравнительного анализа различных вариантов аналитической замены радиальной координаты на ее приведенные аналоги установлено, что максимальную эффективность демонстрируют аналитические формы, основанные на использовании так называемых приведенных Морзовских координат: y(r)=1-exp(-beta*(r-re)); y(r)=1/[1+exp(-beta*r)] и y(r)=exp[-exp(-beta*r)]. Показано, что, и в случае решения системы связанных радиальных уравнений конечно-разностным методом, абсолютная величина коррекции к полной энергии (заданной в интегральной форме) может быть эффективно использована для повышения точности решения без практического увеличения вычислительных затрат. 2. Методом математического моделирования исследовано влияние интерполяционной формы представления собственного дипольного момента на интенсивности высоких обертонов колебательного спектра. Установлено, что если функция дипольного момента задана в табличном виде (как, например, результат неэмпирического расчёта) и если она интерполирована кубическим сплайном с нормальными граничными условиями – то теоретическте интенсивности высших обертонных переходов демонстрируют качественно неправильное поведение так, как выходят на «насыщение» (плато) при больших v. В случае использования для интерполяции обычного (гладкого) полинома соответствующей степени, артефакта «насыщения» не наблюдается, и рассчитанные интенсивности точно следуют нормальному закону распределения вплоть до предела диссоциации за исключением нескольких аномалий. 3. Предложены приближенные модели термической и фотоиндуцированной диссоциации оксидов щелочных и щелочноземельных металлов, образовавшихся в результате химических реакций в неравновесных, возбужденных колебательно-вращательных состояниях. Необходимые для моделирования энергетических и радиационных свойств колебательных уровней кривые потенциальной энергии рассматриваемых молекул (LiO, NaO, KO, MgO, CaO) получены в рамках неэмпирического квантово-химического метода теории возмущений второго порядка (XMC-QDPT2). Наиболее близким по энергии к основному электронному состоянию молекул LiO и NaO оказалось состояние A2Σ+, поэтому возможны два конкурирующих механизма развала молекулы на атомы: 1) столкновительно-индуцированное возбуждение в рамках основного электронного состояния X2Π, и 2) оптически разрешенный электронный переход X2Π → A2Σ+ в континуальную часть колебательного спектра. Пересечение X2Π и A2Σ+ термов локализуется вблизи первого колебательного обертона возбужденного терма A2Σ+, качественно характеризуется большим значением фактора Франка – Кондона и поэтому соответствующий безызлучательный переход же из основного Х2Π терма в состояние A2Σ+ почти не требует термической активации. Электронные конфигурации этих состояний отличаются только заселенностью π- и σ-орбиталей. Обе МО локализованы в значительной степени на атоме кислорода, поэтому матричный элемент СОВ близок к величине 0.5ξ, где ξ = 153 см–1 – константа СОВ для атома кислорода. Структура электронных конфигураций нижних состояний молекул LiO, NaO и KO близка. Однако, в отличие от LiO и NaO первое возбужденное состояние молекулы KO двукратно вырождено (A2Π). Электронный терм A2Π молекулы KO также очень близок к основному терму X2Σ+ в области потенциальной ямы. Диссоциация молекулы KO может проходить по двум каналам: в рамках многоступенчатого столкновительного возбуждения основного электронного состояния X2Σ+ и посредством неадиабатического спин-орбитального перехода из основного состояния X2Σ+ в первое возбужденное состояние A2Π. Второй механизм более вероятен, так как требует меньшей энергии активации из-за большой величины СОВ. Для оксидов щелочно-земельных металлов CaO и MgO возможен распад молекул в рамках основного электронного состояния X1Σ+, а также посредством запрещенного по спину (интеркомбинационного) перехода X1Σ+→a3Π в наиболее близкое по энергии триплетное состояние a3Π за счет значительного эффекта спин-орбитального взаимодействия. 4. Выполнено новое климатическое моделирование палеоатмосферы Марса на основании данных квантовохимических и полинейных спектроскопических расчетов слабосвязанных молекулярных комплексов. Показано, что столкновительно-индуцированное поглощение CO2-H2 и CO2-CH4 молекулярных пар до сих пор было недооценено. Установлено, что метан может выступать в качестве активного парникового газа в палеоатмосфере Марса за счет индуцированного поглощения CO2-CH4 в окне 250-500~cм-1. В атмосфере, содержащей 0.5 бар CO2, достаточно содержания H2 или CH4 на уровне процентов, чтобы средняя температура поверхности повысилась на десятки градусов. Температура поверхности в 273 K достигается при давлении СО2 1.25-2 бар в смеси с 2-10% H2 или CH4. Предполагается, что комбинация метана и водорода, попадающих в атмосферу при геохимической эволюции поверхностных отложений, с вулканическим СО2 могла составить переходную атмосферу Марса около 3.5-4.5 миллиардов лет назад. 5. Методом квазистационарного кинетического моделирования исследованы особенности протекания процессов возбуждения, ионизации, рекомбинации и образования атомарного водорода за фронтом ударной волны в условиях атмосферы долгопериодической переменной звезды типа Миры Кита. Показано, что, при высвечивании за фронтом ударной волны в условиях атмосфер звезд поздних спектральных классов, заселенности дискретных уровней водорода определяются не только текущими значениями температуры и электронной плотности, но зависят от всей предыстории процесса, начиная с локального нагрева на фронте волны.
3 1 января 2017 г.-31 декабря 2017 г. Прецизионное моделирование оптических свойств двухатомных газов при высоких температурах
Результаты этапа: Создан универсальный пакет прикладных программ для прецизионного моделирования электронно-колебательно-вращательных энергий и вероятностей ровибронных переходов двухатомных молекул в широком интервале их термического возбуждения. В основе численного алгоритма лежат различные варианты аналитической замены радиальной координаты на ее приведённые аналоги, что значительно повышает точность решения спектроскопической задачи как в адиабатическом, так и в нетрадиционном для спектроскопии неадиабатическом приближении. Адиабатический вариант программы эффективно работает с учетверенной машинной точностью, что позволяет оценить вероятности очень слабых обертонных переходов молекулы вплоть до энергии ее диссоциации. Методом математического моделирования изучено критическое влияние интерполяционной формы представления собственного дипольного момента молекулы на достоверность предсказания интенсивности высоких обертонов ее колебательно-вращательного спектра. Разработаны приближенные квантовохимические модели термической и фотоиндуцированной диссоциации оксидов щелочных и щелочноземельных металлов, представляющих астрофизический интерес. Впервые выполнены неэмпирические расчеты сечений фотодиссоциации оксида натрия из основного и низколежащих электронных состояний с учетом их взаимного спин-орбитального и электронно-вращательного взаимодействия. В видимой области выполнено моделирование наблюдаемых эмиссионных спектров метеора Бенешова, в результате которого впервые обнаружена оранжевая система оксида кальция. В экзосферах Луны и Меркурия оценено время фотолиза оксида натрия солнечными фотонами. Сравнение с экспериментом указывает на неравновесный характер заселения начальных уровней. Для ряда слабосвязанных молекулярных пар методами квантовой химии, классической механики и статистической физики найдены температурные зависимости их спектральных и термодинамических характеристик. Проанализировано значение рассчитанных физико-химических параметров для моделирования свойств атмосферы Титана. Выполнено моделирование палеоатмосферы Марса на основании данных квантовохимических и полинейных спектроскопических расчетов слабосвязанных молекулярных комплексов. При этом впервые в рамках климатического расчета удалось воспроизвести условия, позволявшие в прошлом существовать на поверхности Марса жидкой воде. Методом квазистационарного кинетического моделирования исследованы процессы возбуждения, ионизации, рекомбинации и образования атомарного водорода за фронтом ударной волны в условиях атмосферы долгопериодической переменной звезды типа Миры Кита. Показано, что, при высвечивании за фронтом ударной волны в условиях атмосфер звезд поздних спектральных классов, заселенности дискретных уровней водорода определяются не только текущими значениями температуры и электронной плотности, но зависят от всей предыстории процесса, начиная с локального нагрева на фронте ударной волны.

Прикрепленные к НИР результаты

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