Лабораторное моделирование высокоэнергетических газо-плазменных процессов представляющих аэрономический и астрохимический интересНИР

Laboratory modeling of high-energy gas-plasma processes of aeronomic and astrochemical interest

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

грант РНФ

Этапы НИР

# Сроки Название
1 31 марта 2023 г.-31 декабря 2023 г. Лабораторное моделирование высокоэнергетических газо-плазменных процессов представляющих аэрономический и астрохимический интерес
Результаты этапа: Установлено, что слабо связанные, квази-связанные и континуальные состояния молекул, локализованные вблизи порога диссоциации, вносят заметный вклад в суммы по состояниям при повышенных температурах, особенно в отсутствии локального термодинамического равновесия. Слабосвязанные молекулярные состояния активно участвуют в образовании устойчивых ансамблей молекул при радиационной ассоциации сталкивающихся атомов, что приводит к эффективному охлаждению начальной среды. В сочетании с прецизионными квантово-химическими расчетами электронной структуры и глобальным неадиабатическим анализом квази-вырожденных состояний, сходящихся к одному и тому же диссоциационному пределу, лазерно-индуцированные спектры позволяют детально изучать структурно-динамические свойства газо-плазменных сред в широком интервале их электронно-колебательно- Проект № 23-13-00207/2023 Страница 11 из 40 вращательного возбуждения. На основании численного эксперимента оценена эффективность различных алгоритмов реконструкции профилей излучательной способности в лазерной плазме и сформулированы требования к разрешению системы сбора излучения для примененеия к получаемым данных обратного преобразования Абеля. Собрана новая собирающая система, сконструирован ввод вращения для устранения прецессии мишени в вакуумной камере. В диапазоне экспериментальных условий (давление в вакуумной камере, задержка наблюдения, энергия испаряющего лазерного импульса), близких к таковым при сгорании конструкционных материалов в верхних слоях земной атмосферы были зарегистрированы спектры лазрено-индуцированной плазмы при испарении мишени из металлического титана. В этих условиях наблюдается интенсивная эмиссия как самих атомов титана, так и компонент электронной системы B3Π-X3Δ молекул TiO. Апробировано 8 схем флуоресценции для пространственно-разрешенной диагностики лазерной плазмы. Проведена оценка температуры лазерной плазмы с помощью флуоресценции при давлении 100 Торр и задержке 5 мкс. Получены карты распределения температуры в лазерной плазме. Было обнаружено существенное различие с данными, получаемыми по методу графика Больцмана, причина которго требует дополнительных исследований на следующих этапах выполнения проекта. Наиболее горячая часть плазмы располагается на перифирии, что может быть связано с прохождением ударной волны. С использованием оригинальных алгоритмов участников проекта по моделированию спектров плазмы создан обучающий набор из 100 000 синтетических спектров плазмы, расссчитанных для различного химического состава в широком диапазоне параметров лазерной плазмы и регистрирующей системы. Отработана методология аппроксимации спектров красной системы CaO (A1Σ+ –X1Σ+) с возможностью использования нескольких температур и учетом эффекта самопоглощения. Для этого имеющийся алгоритм моделирования спектров плазмы в предположении ЛТР был расширен на случай трех различных температур: электронной, колебательной и вращательной (при моделировании первая приравнивалась к вращательной температуре). В рамках одноканальной теории квантового дефекта и аналитической модели остовно-поляризационного потенциала исследовано асимптотическое поведение потенциальных кривых возбужденных электронных состояний молекулярного водорода на малых и средних межъядерных расстояниях R. Показано, что электронные термы, обладающие gerade и ungerade симметрией, гладко сходятся в пределе объединенного атома к синглетным и триплетным состояниям атома гелия 2HeI, находящегося в S,D,… и P,F,… состояниях, соответственно. Проведено однозначное отнесение рассматриваемых состояний по главному квантовому числу n и угловому моменту ридберговского электрона l, как на диссоциационном пределе, так при R->0. Из анализа степени гладкости построенных функций квантового дефекта и их 1-ой производной по R удалось установить места некорректного отнесения энергий, полученных в рамках прецизионных ab initio расчетов для состояний с высоким значением l. Получены аналитические выражения для неадиабатических матричных элементов электронно-колебательного (радиального) и электронно-вращательного (углового) взаимодействия между электронными состояниями двухатомных молекул, принадлежащими двум ридберговским сериям с угловым моментом электрона l отличающимся на 2. На примере молекулы водорода изучен конкурирующий вклад связанных и континуальных ридберговских электронных состояний в величину адиабатической коррекции и неадиабатического сдвига ровибронных уровней энергии низколежащих ридберговских состояний двухатомных молекул. Получены приближенные (модельные) оценки констант скоростей нейтраль-нейтральных и ион-нейтральных химических реакций с участием атомов N, C+, S+ и Si+ и заселенных неравновесно молекул водорода, рассчитанные в рамках теории активированного комплекса и активных столкновений, соответственно. Выполнен критический анализ данных о строении и динамике астрохимически важных молекул полициклических ароматических углеводородов (ПАУ) в S1 и Т1 электронных состояниях. Проведен анализ методик получения спектров возбуждения флуоресценции и фосфоресценции молекул, охлажденных в сверхзвуковой струе инертного газа, с целью определения оптимальной конструкции приставки для наблюдения сенсибилизированной фосфоресценции. Созданы чертежи новой вакуумной камеры и азотной ловушки к ней с учетом временных характеристик фосфоресценции исследуемых объектов, опубликованных в литературе. Новая вакуумная камера с азотной ловушкой и второй ФЭУ для наблюдения фосфоресценции были изготовлены и смонтированы в установке. Проведена модернизация вакуумного клапана, являющегося одним из ключевых узлов установки для получения спектров возбуждения флуоресценции молекул, охлажденных в сверзвуковой струе. Выполнена модернизация электронной схемы сбора данных для подключения второго ФЭУ. Зарегистрированы тестовые спектры возбуждения флуоресценции молекулы нафталина, охлажденного в сверхзвуковой струе аргона.
2 1 января 2024 г.-31 января 2024 г. Лабораторное моделирование высокоэнергетических газо-плазменных процессов представляющих аэрономический и астрохимический интерес
Результаты этапа: 1. Состав мишени из смеси оксида железа (II, III) и карбоната кальция , при испарении которой спектры лазерной плазмы и болида Бенешев схожи, определенный с помощью калибровки по сигналам Ca и Fe, соответствует соотношению Fe/Ca 84:1 ,что существенно отличается от LL-хондритов и оценок для болида Бенешев (14:1-17:1). Это может быть связано с селективным испарением компонентов метеороида по мере нагрева его поверхности, что приводит к значительному обогащению плазменного облака при его движении в атмосфере Земли легкоатомизуемыми элементами, в частности кальцием. 2. На основании молекулярных постоянных рассчитаны вероятности переходов оранжевой системы CaO (≈24000 переходов). Качественное сопоставление спектра CaO высокого разрешения с рассчитанным модельным спектром в приближении ЛТР позволяет говорить только об общем совпадении положений соответствующих полос, что говорит о необходимости расширения количества рассматриваемых переходов данной системы как минимум на порядок для более реалистичного её описания. 3. При поисковых работах по обнаружению полос SiO2 зарегистрированы диффузные полосы вблизи 310 нм, в диапазонах 427-455 нм, 610-620 нм, об отнесении которых нет литературных данных, а энергетические характеристики могут говорить об их принадлежности к SiO2. На следующем этапе выполнения проекта будет изучаться влияние давления кислорода на изменение интенсивности этих полос. 4. Показана возможность предсказания с помощью нейронной сети температуры и электронной плотности по калиброванному спектру лазерной плазмы. 5. Неадиабатическая модель связанных колебательных каналов была использована для генерации полного списков линий радиационных переходов между X2Σ+, B2Σ+ и A2Π состояниями наиболее распространенных изотопологов CN радикала. Рассчитанные для A−X, B−X и B−A электронных переходов коэффициенты спонтанного излучения Эйнштейна были использованы для оценки радиационных времен жизни возбужденных ровибронных уровней и прогнозирования коэффициентов ветвления вибронных B→A/X переходов. Необходимые для расчета радиационных характеристик собственные дипольные моменты электронных состояний и дипольные моменты электронных переходов были получены многореференcными методами конфигурационного взаимодействия. 6. На основании функций квантового дефекта электронных состояний молекулы H2 с угловым моментом ридберговского электрона l<5 и его проекцией Λ=0-3 получены оценки неадиабатических матричных элементов радиального и углового связывания между ридберговскими состояниями одинаковой симметрии, а также радиальной функции адиабатической коррекции с использованием аналитических выражений молекулярной теории квантового дефекта. Рассчитаны регулярные неадиабатические сдвиги и энергии экспериментально наблюдаемых колебательных уровней низколежащих a3Σg+ и e3Σu+ состояний H2, D2 и T2 изотопологов в не вращающейся молекуле. Многоссылочным методом конфигурационного взаимодействия с использованием остовно-поляризационного потенциала и насыщенных диффузными функциями полноэлектронных базисных наборов выполнены расчеты кривых потенциальной энергии и собственных дипольных моментов ридберговских состояний эксиплексной молекулы NaHe. Полученные ab initio результаты количественно подтвердили волнообразное поведение межатомных потенциалов, предсказанное ранее в рамках теории неупругого рассеяния. 7. Выполнены квантово-химические расчеты энергетических и радиационных свойств синглетных и триплетных термов молекулы СаО, спин-разрешенные и интеркомбинационные электронные переходы, между которыми образуют «оранжевую» систему полос. Полученные результаты использованы для генерации списка линий, включающие частоты ровибронных переходов и коэффициенты Эйнштейна спонтанной эмиссии, который использовался для моделирования высокотемпературных спектров ЛИЭС. 8. Получены «поуровневые» константы сечения неупругих столкновений для переходов из всех связанных уровней основного электронного состояния молекулы водорода в ряд низколежащих колебательно-вращательных уровней молекулы HF. 9. В рамках по модификации установки регистрации спектров возбуждения флуоресценции молекул ПАУ выбрана оптимальная геометрия и материалы штока (полиацетат, фторопласт), обеспечивающие полную герметичность клапана в закрытом состоянии и хорошую стабильность работы в импульсном режиме. Изготовлен и встроен в газовую магистраль перед импульсным вакуумным клапаном подогреваемый резервуар для твердых образцов исследуемых соединений. Зарегистрирован и проанализирован тестовый спектр возбуждения флуоресценции молекулы нафталина, охлажденной в сверхзвуковой струе аргона. 10. В рамках теоретического исследования реакций образования молекул ПАУ предложен многоступенчатый механизм последовательного роста молекул ПАУ в низкотемпературных условиях МЗС через электронно-возбужденный мономер бензола, его эксимер, а также их катионов, рассмотренных в качестве промежуточных продуктов фотоиндуцированной десорбции бензольной ледяной мантии, образовавшейся на частицах космической пыли. Выполнены предварительные расчеты электронной структуры возбужденных электронных состояний мономера и димера бензола, выступающих в качестве прекурсоров реакций роста ПАУ.
3 1 января 2025 г.-31 января 2025 г. Лабораторное моделирование высокоэнергетических газо-плазменных процессов представляющих аэрономический и астрохимический интерес
Результаты этапа: 1. Наиболее отчётливо рост сигналов в спектрах при испарении кремния с ростом давления кислорода проявляется в области 445 нм, что с учетом немногочисленных литературных данных, указывает на возможные полосы молекул SiO2. Вблизи длины волны 610 нм можно наблюдать слабый сигнал, который значительно усиливается с увеличением давления кислорода. Как и в случае с диоксидом кремния, данная спектральная особенность может быть отнесена к диоксиду железа FeO2. 2. Имеющиеся данные о состояниях и переходах TiO из базы ExoMol позволяет оценивать равновесную температуру при атмосферном давлении на поздних временах наблюдения плазмы (18 до 30 мкс) с точностью на уровне 25 %. Основными факторами, объясняющими расхождение экспериментальных и модельных спектров являются наличие интенсивных атомных линий. Систематическое занижение интенсивности синтетических спектров в области 550-560 нм относительно экспериментальных данных может указывать на необходимость уточнения данных по вероятностям переходов для системы B(3PI) - X(3DELTA). 3. При работе с металлической мишенью при пониженном давлении эмиссионные спектры, содержащие атомных и/или ионных линий как правило оказываются интенсивнее сигналов рассеяния. Однако, для давлений до 10 Торр и ранних временах наблюдения плазмы результаты определения температуры и электронной плотности плазмы по спектрам некогерентного Томсоновского рассеяния коррелируют с результатами диагностики по собственному излучению плазмы. 4. Анализ спектров флуоресценции молекулы CaO в лазерной плазме показал, что в области «синей» полосы (360-440 нм) сигналы можно отнести к колебательно-вращательным уровням B1Π состояния; переходам между, предположительно, состояниями C1Σ и X1Σ (без идентификации). Также для части сигналов недостаточно литературных данных для отнесения даже к электрнным состояниям. 5. Улучшенная архитектура нейросети и расширенный обучающий набор позволили получить модель, точности предсказания которой близка по порядку величины к погрешности определения экспериментальных сигналов, используемых для соответствующих расчетов. 6. Из-за значительного уширения линии Li надежная оценка поступательной температуры возможна на относительно поздних временах существования лазерной плазмы (1-5 мкс). В этот момент поступательная температура может достигать 40 000 К, далее быстро уменьшаясь по мере расширения плазмы до 10 000–15 000 К. Можно отметить, что с ростом давления наблюдается более выраженное снижение температуры при временной эволюции плазмы, что отражает более быстрый теплообмен и торможение частиц при ее расширении. 7. Построена единая неадиабатическая модель описания энергетических, радиационных и магнитных характеристик X~A~B комплекса для всех изотопологов молекулы CN, которая учитывает в явном виде внутримолекулярные спин-орбитальные и электронно-вращательные взаимодействия между дублетными состояниями комплекса и с близлежащими квартетными Σ+ и Π термами, а также диагональную Борн-Оппенгеймеровскую коррекцию. Предложенная масс-инвариантная модель позволяет существенно расширить энергетическую область адекватного описания оптических свойств X~A~B комплекса и тем самым увеличить точность моделирования эмиссионных спектров радикала СN при высоких температурах и в неравновесных условиях. 8. На примере ЛИЭС спектров радикала CN, полученных при испарении углерод-азот-содержащей твердой мишени на воздухе, обнаружено разделение колебательной и вращательной температуры (разница до 2000 К) на средних и поздних временах наблюдения после лазерного импульса. Учет эффекта самопоглощения позволил существенно улучшить качество описания экспериментальных данных, в особенности в области наиболее “сильных” B2Σ+-X2Σ+ переходов полосы 0-0. 9.В рамках аналитической одноканальной теории квантового дефекта исследованы вклады дискретных и континуальных ридберговских состояний в адиабатическую коррекцию и регулярные неадиабатические сдвиги вибронных уровней низколежащих термов двухатомной молекулы. Установлено, что парциальный вклад электронного континуума может достигать нескольких десятков процентов даже для первых членов ридберговской серии. Показано, что модельные оценки суммарных вкладов в адиабатическую коррекцию e3Σ_u^+(3pσ_u) и a3Σ_g^+ (2sσ_g) состояний Н2 хорошо согласуются с высокоточными ab initio данными на малых и средних межъядерных расстояниях, а рассчитанные с учетом адиабатической коррекции и неадиабатического сдвига энергии вибронных уровней с экспериментальными термами его гомоядерных изотопологов. 10.Выполнены неэмпирические расчеты адиабатических функций потенциальной энергии, собственных дипольных моментов и дипольных моментов разрешенных по спину электронных переходов для основного и ряда низколежащих возбужденных синглетных состояний изолированной молекулы монооксида кальция. На основании систематического анализа полученных в настоящей работе ab initio результатов и имеющихся в литературе альтернативных теоретических и экспериментальных данных, сделан вывод, что наиболее адекватные результаты расчета электронной структуры СаО можно получить только в рамках метода MR-AQCC причем для низколежащих состояний и в относительно узком интервале межъядерных расстояний. 11.Метод произведения неортогональных волновых функций (NOPA) объединен с теорией возмущений в полном активном пространстве (CDAS-PT2) для исследования возбужденных электронных состояний молекулы димера бензола. Предложен общий механизм фотоиндуцированного формирования эксимерных состояний, с их последующей стабилизацией за счет колебательной релаксации до минимума оптически неактивных состояний, а также параллельного преддиссоциативного распада, вызванного внутримолекулярной связью с множеством отталкивательных состояний. 12.Систематические кватовохимические расчеты электронной структуры уровня TD-B3LYP/TZVP показали заметное влияние пространственной топологии и включения атомов азота в полициклические ароматические углеводороды (ПАУ) на особенности строения их УФ– и ИК-спектров. 13. Методами астрохимического моделирования исследовано влияние неравновесного колебательно-вращательного возбуждения молекулярного водорода на итоговую концентрацию молекул HF, синтезируемых в реакции F + H2 → HF + H в фотодиссоциационной области межзвездной среды.

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

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

Прикрепленные файлы


Имя Описание Имя файла Размер Добавлен