Создание средств компьютерного моделирования фотохимических изомерных превращений сложных молекулярных систем тема автореферата и диссертации по химии, 02.00.04 ВАК РФ
Завалий, Максим Владимирович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Москва
МЕСТО ЗАЩИТЫ
|
||||
2005
ГОД ЗАЩИТЫ
|
|
02.00.04
КОД ВАК РФ
|
||
|
На правах рукописи
Завалий Максим Владимирович
СОЗДАНИЕ СРЕДСТВ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ ФОТОХИМИЧЕСКИХ ИЗОМЕРНЫХ ПРЕВРАЩЕНИЙ СЛОЖНЫХ МОЛЕКУЛЯРНЫХ СИСТЕМ
02.00.04 - физическая химия, 02.00.02 - аналитическая химия
АВТОРЕФЕРАТ
диссертации на соискание ученой степени кандидата физико-математических наук
Москва - 2005
Работа выполнена в лаборатории молекулярного моделирования и спектроскопии Института геохимии и аналитической химии им. В.И. Вернадского Российской академии наук
Научные руководители:
член-корреспондент РАН
доктор физико-математических наук, профессор
Грибов Л.А. Баранов В.И.
Официальные оппоненты:
член-корреспондент РАН
доктор физико-математических наук, профессор
Чибисов А. К. Юхневич Г.В.
Ведущая организация: ГНЦ РФ "Научно-
исследовательский физико-химический институт им. Л.Я. Карпова"
Защита состоится < ОЬ. 05, > в </£_> час. на заседании диссертационного совета Д.002.109.01 в Институте геохимии и аналитической химии им. В.И. Вернадского РАН по адресу: 119991, ГСП-1, Москва, В-334, ул. Косыгина, 19.
С диссертацией можно ознакомиться в библиотеке Института геохимии и аналитической химии им. В И. Вернадского РАН.
Автореферат разослан «¿У2005 г.
Учёный секретарь совета доктор химических наук
Кубракова И.В.
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы.
Сегодня оптическая спектроскопия, благодаря открытиям в области техники коротких и сверхкоротких лазерных импульсов, дала в руки химикам качественно новый инструмент - импульсы света фемтосекундной длительности. Они принесли в химию новые возможности и стали толчком для развития нового направления науки - фемтосекундной спектроскопии, основанной на резонансном возбуждении отдельных высоколежащих уровней энергии молекул с последующим наблюдением их высвечивания. Это позволило получать очень развитые спектры с временным разрешением, в которых можно наблюдать как изменения интенсивностей излучения или поглощения во времени на определенных длинах волн, так и сразу весь спектр в определенные интервалы времени.
Фемтосекундная технология позволяет изучать химические реакции с высоким временным и пространственным разрешением. Стало возможным наблюдение динамики реакций и целенаправленное управление их ходом с высоким временным разрешением, что открывает значительные перспективы в области качественного и количественного анализа, синтеза новых химических соединений и предсказания их реакционной способности. Импульсы фемтосекундной длительности дают мощный инструмент для исследования элементарных превращений на уровне процессов разрыва и образования связей, переноса заряженных частиц в мономолекулярных (например, фото диссоциация, изомеризация) и бимолекулярных реакциях. Фемтосекундная спектроскопия позволяет регистрировать временные зависимости интенсивностей излучения или поглощения света молекулой в ходе протекания внутримолекулярных и межмолекулярных процессов (в том числе и безызлучательных), что дает возможность изучать изомер-изомерные преобразования, а они служат основой для передачи энергии и информации, сигнала во внутримолекулярном пространстве.
Применение техники фемтосекундной спектроскопии в аналитической и физической химии является нетривиальной и мало изученной задачей. Высока сложность интерпретации динамических спектров, которые имеют не только тонкую структуру вдоль шкалы частот, но и сложным образом зависят от времени. При этом вероятны случаи изменения интенсивности излучения или поглощения света молекулой в ходе протекания внутримолекулярных и межмолекулярных процессов. Для выбора подходящих условий принципиально использовать возможности компьютерных экспериментов, так как только с их помощью можно заранее, варьируя заданную ситуацию, выбрать наиболее оптимальный путь постановки натурных экспериментов.
Попыток проведения моделирования динамики фотохимических изомер-изомерных молекулярных преобразований, возникающих в результате облучения молекулы короткими фемтосекундными лазерными импульсами, для реальных моделей сложных молекул никогда раньше не предпринималось. Между тем такая задача весьма актуальна и есть все основания для ее постановки. Успех работ по расчету динамических спектров для одиночных молекул без возможных их фотохимических преобразований и по применению этих спектров к вопросам
качественного и количественного анализа в аналитической химии (в частности, без применения стандартных образцов состава), а также теоретических работ по описанию межизомерных превращений, подтверждает это [1,2].
Цель настоящей работы состоит в создании средств моделирования динамики фотохимических процессов, связанных с изомерными превращениями молекул, применимых в физической и аналитической химии.
Для достижения этой цели были поставлены следующие задачи:
1. Развитие теории и разработка методов компьютерного моделирования фотохимических изомерных преобразований и квантовых биений состояний сложных молекул на базе детального описания их моделей, обеспечивающих проведение предсказательных расчетов в масштабе реального времени
2. Проведение на основе разработанных методов модельных расчетов и анализ получаемой спектральной картины и кинетики образования продуктов фотохимической реакции с целью выявления общих закономерностей спектральных проявлений процесса передачи сигнала и энергии внутри сложных молекул с учетом изомеризации.
Научная новизна работы состоит в создании теоретической базы качественного и количественного анализа в фотохимии:
1. Развита теория и разработано эффективное программное обеспечение для моделирования динамики фотохимических изомер-изомерных преобразований.
2. Выявлены общие закономерности межизомерных преобразований, в том числе:
- зависимость кинетики межизомерных преобразований и их спектральных проявлений от начальных условий;
- условия возможности наблюдения спектров флуоресценции образующихся изомеров и специфика проявления квантовых биений в колебательной структуре спектров;
- временные и количественные характеристики образования изомеров.
3. Впервые исследованы процессы цепочечных изомер-изомерных перестроек. связанные с передачей сигналов и энергии во внутримолекулярном пространстве.
4. Проанализирована возможность применения развитых методов и программ для качественного и количественного спектрального анализа
Практическая значимость. Разработан сервисный, пригодный для широкого использования программный комплекс, обеспечивающий проведение моделирования динамики фотохимических процессов, включающих межизомерные переходы, и их спектральных проявлений в масштабе реального времени для сложных молекулярных систем. Создан теоретический базис для проведения компьютерных экспериментов в физической и аналитической химии, опережающих натурный эксперимент и позволяющих осуществлять его целенаправленное планирование.
На защиту выносятся: 1. Формулировка и решение задачи о моделировании фотохимических межизомерных превращений молекул, возникающих в результате произвольного воз-
мущения, действующего на молекулу и приводящего к начальному заселению ее возбужденных электронно-колебательных уровней энергии.
2. Вычислительные методы, алгоритмы и программы для персонального компьютера и суперкомпьютера, обеспечивающие проведение моделирования рассматриваемых процессов в масштабе реального времени.
3. Результаты конкретных расчетов для изомер-изомерных преобразований сложных молекул (пентадиен, гептадиенилбензол) с демонстрацией и анализом динамики населенностей уровней и динамических спектров.
4. Прогноз возможности расширения методов спектрального анализа и распространения их в область фотохимических превращений вещества.
Апробация работы.
Материалы диссертации доложены на 3-й Всероссийской конференции «Молекулярное моделирование» (г. Москва, 2003 г.) и на Конференции «Информационно-вычислительные технологии в фундаментальных и прикладных физико-математических исследованиях» (г. Москва. 2004 г.).
Публикации.
Основное содержание диссертации опубликовано в трех научных статьях и двух тезисах докладов на научных конференциях.
Структура и объём работы.
Диссертация состоит из Введения, четырёх глав, выводов, списка литературы и двух приложений. Общий объём работы составляет 239 страниц текста, включая 29 рисунков, 3 таблицы, список литературы, содержащий 136 наименований, приложения А и Б, объём которых составляет 92 страницы.
СОДЕРЖАНИЕ РАБОТЫ
Во Введении обосновываются актуальность и новизна темы исследования, формулируются основные задачи диссертационной работы и положения, выносимые на защиту, кратко излагается содержание диссертации.
Глава 1. Обзор литературы
На основе обзора литературы и проведенного анализа современных методов экспериментальной и теоретической электронно-колебательной спектроскопии и состояния теории качественного и количественного спектрального анализа выяснено следующее.
1. В экспериментальной спектроскопии разработан новый метод регистрации вибронных спектров сложных молекул с временным разрешением, использующий в качестве возбуждающего излучения фемтосекундные лазерные импульсы. Применение метода имеет огромное значение в фотохимии, поскольку становится возможным наблюдение динамики реакций с высоким временным разрешением и целенаправленное управление их ходом. Фемтосекундная технология позволяет управлять выходом продуктов элементарной химической реакции, исследовать элементарные превращения на уровне процессов разрыва и образования связей, переноса заряженных частиц в мономолекулярных и бимолекулярных реакциях.
2. Существуют все предпосылки для создания методов моделирования динамики фотохимических превращений и их спектральных проявлений. Базовыми являются теоретические разработки по определению вероятностей безызлучатель-ных и межизомерных переходов, путей химических реакций, а также развитые методы моделирования динамических вибронных спектров без учета фотохимических превращений. Для решения этих задач вполне приемлем параметрический метод расчета тонкоструктурных электронно-колебательных спектров многоатомных молекул и определения моделей молекул в возбужденных состояниях
[3].
3. Сложность постановки динамического фотохимического эксперимента с пи-ко- и фемтосекундным временным разрешением требует его разумного планирования, выбора оптимальных условий и т. д. Это может быть осуществлено только на основе компьютерных экспериментов предсказательного характера (более дешевых и быстрых). Необходимо развитие соответствующих методов моделирования, компьютерных программ, призванных обеспечить проведение такого вычислительного эксперимента.
4. Применение теоретических методов динамической вибронной спектроскопии с учетом межмолекулярных преобразований сложных молекул для решения аналитических задач является весьма перспективным, может существенно расширить возможности анализа вещества и распространить эти методы в область фотохимических превращений веществ.
Показано, что моделирование динамики фотохимических процессов является весьма актуальной научной проблемой, для решения которой необходимо, во-первых, развитие теории, разработка компьютерных методов предсказательного, в масштабе реального времени моделирования изомерных преобразований сложных молекул с учетом квантовых биений, во-вторых, выявление закономерностей спектральных проявлений процесса передачи сигнала и энергии внутри сложных молекул с учетом изомеризации и анализа получаемой спектральной картины и кинетики образования продуктов фотохимической реакции.
Собственно оригинальный материал диссертации изложен в главах 2-4.
Глава 2. Развитие теории моделирования фотохимических межизомерных преобразований с учетом квантовых биений для сложных молекул
Развит подход к описанию изомер-изомерных преобразований многоатомных молекул, основывающийся на том, что оно становится возможным и наиболее вероятным при резонансе уровней изомеров. Под изомерами могут пониматься и структуры, образующиеся на первой стадии реакции разложения. В результате резонанса появляются дублеты с разностью энергетических состояний Д£ = 2£„5„, где £м - энергия резонирующих электронно-колебательных состояний подсистем и - интеграл перекрывания соответствующих волновых функций. При наложении на систему возмущения возникают квантовые биения и формируется волновой пакет, который при значительно меньшей энергии любого перехода передается функцией
Здесь функции у/х и относятся к стационарным состояниям резонирующих подсистем, сомножитель ехр(-г'й>„ ?) обеспечивает возможность резонансного поглощения световой энергии или спонтанного излучения, величины Аа> = АЕ/ 2Й, а>й = (й>, + о)2) / 2, (0,= Е11Ь и Ех=Ег (резонанс уровней). Этот волновой пакет, несмотря на осциллирующие во времени сомножители, в хорошем приближении (АЕ мало) отвечает стационарному по энергии состоянию Величина характеризует частоту квантовых биений между состояниями комбинирующих подсистем (в дальнейшем для упрощения записи будем обозначать ее просто , Значения частот квантовых биений для конкретных структур сложных молекул могут быть вычислены прямыми методами.
л Сказанное позволяет ввести
модельное представление, согласно которому обобщенная задача
рассматривается для системы уровней невзаимодействующих состояний пар изомеров, определяемых обычным образом, дополненной общими уровнями, соответствующими резонансам и объединяющими пары изомеров (см рис 1а). Если для резонирующих уровней и биения, т.е. переходы из одной резонансной формы в другую наблюдаться не будут. Если же то
и возникнет состояние, которое будет характеризовать безызлучательный (и без изменения энергии) переход от начальной изомерной формы (состояние
к другой (состояние Квадрат
модуля волновой функции 4х этого со-
Изомер 1
Изомер 1
Изомер 3
Рис 1 Пример системы уровней и переходов для трех изомерных форм молекулы для моделей изомеризации общего характера (модель 1) (а) и в форме бе-зызлучательного перехода (модель 2)
(б).
СТОЯНиет имеет имя*
М2 = (у/,2 соя2 со { +у/2г ьт2 со0 (1)
Опираясь на данный теоретический подход, предложен способ решения задачи о динамике населенностей состояний и интенсивностей спектров сложных молекул с учетом межизомерных преобразований.
Ансамбль невзаимодействующих совершенно одинаковых молекул можно рассматривать как квантовую систему, обладающую тем же. что и одиночная молекула, набором уровней энергии, но характеризующихся еще и их заселенностью. При возбуждении (поглощение света) молекул ансамбля или их высвечивании (люминесценция) будут совершаться переходы между уровнями энергии, совпадающими с уровнями одиночной молекулы, приводящие к изменению их заселенностей. Кинетические уравнения для населенностей уровней п,(Х), пред-
ставляющие собой дифференциальные уравнения первого порядка и образующие систему дифференциальных уравнений, будут иметь вид:
(2)
где w, = X + wi0 — полные суммарные вероятности переходов из i-го состояния во все нижележащие возбужденные и основное (wj0); некоторые коэффициенты м>д могут зависеть от времени осцилляторно. Начальные условия (л,(0), i =1,2,...,N) могут быть различными в зависимости от характера возбуждения молекулы (резонансное, широкополосное и др ), вероятностей поглощения и спектрального состава падающего излучения.
Предположим, что импульс возбуждения обладает длительностью, меньшей характерных времен релаксации возбужденных состояний (обратных значениям полных суммарных вероятностей и достаточно малой интенсивностью, позволяющей не принимать во внимание нелинейные эффекты. При этом после окончания акта импульсного фотовозбуждения создается неравновесное распределение населенностей, характеризуемое начальными условиями и,(0), временная динамика которого при отсутствии внешнего поля при t>0 описывается системой (2).
Для решения системы дифференциальных уравнений (2), прежде всего, необходимо определить величины вероятностей переходов (поглощения, излучения и безызлучательных). Ограничимся моделью изолированной молекулы (разреженный газ, сверхзвуковая струя), учитывая переходы в другие изомерные формы и не принимая во внимание межмолекулярные взаимодействия (и обусловленные ими безызлучательные переходы). Следует отметить, что принятые допущения носят чисто модельный характер и не лишают сам метод расчёта спектров с временным разрешением общности, а лишь отвечают современному состоянию теории и позволяют решать задачу о вычислении вероятностей в рамках хорошо отработанных методов.
Данная модель (модель 1) содержит лишь стационарный набор уровней; возможны оптические переходы на общий уровень с других уровней с вероятностями wn, определяемыми квадратами обычных матричных элементов диполь-ных переходов между внутренними для каждого из изомеров состояниями,
умноженными на
На основе этой модели общего характера (смешанные состояния изомеров) построена вторая модель, интерпретирующая эффект квантовых биений состояний как безызлучательный переход (модель 2). Получено выражение для вероятности безызлучательного перехода wlu . Например, для w'J (см. рис. 16) оно принимает вид:
dn _ л _
—L = -Wini+ Y i = 1,2,...,N,
dt
= в>а sin 2a>01\ 1 + 0.2sf
/ "
> \
'31
+ /*0)И< 12)4 яп2й>„Д. (3) ®„(«i,,+wiJ1)J "j
Вероятность безызлучательного перехода, будучи прямо пропорциональной частоте биений, что отражает физический смысл этой частоты как характеристики времени перехода системы между связанными состояниями, зависит от времени осцилляторно, моделируя таким образом эффект биений, и включает в себя вероятности оптических внутримолекулярных переходов, что обеспечивает согласованность заселенностей связанных состояний в процессе квантовых биений. Качественно новые особенности вероятности безызлучательного перехода по сравнению с обычными подходами обусловлены учетом в данной
модели взаимосвязи резонирующих уровней, внутренне им присущей. Собственно величина вероятности безызлучательного перехода не является параметром модели молекулярной системы. В качестве такового выступает лишь частота биений, непосредственно связанная с параметром, характеризующим внутримолекулярные межизомерные взаимодействия (соответствующий недиагональный элемент энергетической, матрицы).
Данная модель 2 адекватна более общей модели 1. Вопрос о степени эквивалентности этих моделей и выборе наиболее оптимальной для моделирования сложных многоуровневых и многоизомерных молекулярных структур может быть решен на основе модельных расчетов реальных систем.
Зная вероятности электронно-колебательных переходов и величины частот квантовых биений, можно решить систему кинетических уравнений для на-селенностей уровней молекул и определить динамику системы, а также построить электронно-колебательный спектр флуоресценции с временным разрешением, который может быть представлен в виде трехмерной поверхности -зависимости интенсивности / от частоты излучения и времени наблюдения t:
Показано, что в рамках данного подхода для решения задачи о динамике населенностей состояний и интенсивностей спектров сложных молекул с учетом межизомерных преобразований можно построить соответствующие вычислительные методы на базе уже существующих в теории традиционных (стационарных) и динамических электронно-колебательных спектров молекул. Тем самым обеспечивается необходимая преемственность в теории вибронных (и колебательных) спектров, основанная на единых общепринятых и наиболее полных молекулярных моделях. Соответствующие средства моделирования динамических фотохимических межизомерных преобразований включаются в единый программный комплекс.
Результаты выполненного анализа составляют необходимую базу для программной реализации средств моделирования динамики фотохимических межизомерных преобразований и их спектральных проявлений и проведения соответ-
/,(в-,0 = л®г -W, и,(О-
(4)
ствующих компьютерных экспериментов для моделей реальных молекул, чему посвящены главы 3 и 4
Глава 3. Создание средств моделирования фотохимических превращений сложных молекул
Моделирование динамики фотохимических превращений основывается на решении самостоятельных фундаментальных задач молекулярной спектроскопии, а именно нахождение электронных уровней молекулярных моделей изомеров, определение колебательных подуровней молекулярных моделей в основном и электронно возбужденных состояниях, вычисление вероятностей электронно-колебательных переходов Только решив все эти задачи, можно переходить к моделированию динамики фотохимических превращений (формированию и решению систем кинетических уравнений, расчету и построению многомерных динамических спектров и временных зависимостей населенностей уровней изомеров) Общая схема расчета показана на рис 2
На основе разработанных методов и алгоритмов создан комплекс программ сервисного типа для моделирования межизомерных фотохимических преобразований сложных молекул на ПК и суперкомпьютере В программном комплексе реализовано два модельных подхода к описанию изомер-
изомерных преобразований Комплекс обеспечивает построение временных зависимостей заселенно-стей всех электронно-колебательных V ровней многоизомерной молекулярной системы и полных трехмерных (частота, интенсивность, время) спектров поглощения и флуоресценции, а также любых их двумерных сечений при любом (резонансном, широкополосном и др) характере возбуждения системы
Рис 2 Схема моделирования динамики фотохимических изомер-изомерныч превращений
Возможности комплекса характеризуются следующими параметрами: размер молекул - до 52 атомов, число нормальных колебаний каждого типа симметрии не более 80, число блоков симметрии не более 15: число изомерных форм не ограничено; число электронно-колебательных состояний не ограничено: количество переходов не ограничено. Ограничения на размер молекул не принципиальны и легко снимаются при использовании программного комплекса LevWindows для подготовки исходных данных об уровнях энергии и вероятностях переходов между ними. Такая возможность предусмотрена.
Временные затраты на моделирование динамики процессов изомеризации и их спектральных отображений во многом определяются сложностью конкретной решаемой задачи и характеризуются следующими величинами - количество изомеров в задаче, количество электронно-колебательных уровней молекулярной модели, частоты квантовых биений, различия времен «экспоненциальной» релаксации состояний и периодов квантовых биений. Типичное время расчета для крупных молекул (40-50 атомов) колеблется в пределах от нескольких секунд (для одного электронно-возбужденного состояния и без учета колебательной релаксации) до нескольких минут. Заложенные в основу комплекса специальные алгоритмы позволяют решать задачу в масштабе реального времени на персональных компьютерах, а в более сложных случаях - на суперкомпьютерах типа МВС-1000.
Отметим некоторые особенности разработанных алгоритмов и программ.
Использовались приближенные численные методы для интегрирования систем кинетических уравнений очень высокой размерности (104 и более) - метод Эйлера для модели 1 и метод Рунге-Кутта 4-го порядка для модели 2. Это обеспечило возможность решения задачи в общем случае и дало существенный выигрыш по времени.
Заметим, что число электронно-колебательных состояний велико (заведомо N >104 для молекул, содержащих порядка 30 и более атомов), а при моделировании динамики межизомерных преобразований, и в особенности при решении обратных задач, необходимо многократно осуществлять расчет при различных параметрах молекулярных моделей, вероятностях электронно-колебательных переходов и безызлучательных, частотах квантовых биений, начальных условиях, временных интервалах и т.д. Поэтому метод решения кинетических уравнений должен обладать высокой скоростью для проведения быстрых расчетов в случае достаточно сложных систем. Видно, что задача сложна и требует оптимизаций и разработки специальных программных алгоритмов.
Проведенный ранее анализ показал, что наиболее перспективным является приближенный численный метод интегрирования кинетических уравнений. Для решения может быть использован простейший метод Эйлера (при условии выбора достаточно малого шага интегрирования). Он хорошо зарекомендовал себя при моделировании динамических спектров одиночных молекул. Однако учет межизомерных переходов приводит к качественному изменению системы кинетических уравнений - появляются дополнительные осциллирующие во времени (из-за квантовых биений) коэффициенты, причем временные характеристики ос-
цилляций и релаксации возбужденных состояний могут различаться на порядки величин. В ряде случаев (модель 2) это приводит к расходимости алгоритма даже при относительно малом шаге итераций и требует поиска иного, более точного метода, который снял бы возникающую проблему и не приводил бы к существенному уменьшению эффективности метода моделирования (прежде всего по затратам времени). Показано, что этим требованиям удовлетворяет метод Рунге-Кутта высокого порядка. Именно он был реализован в программном комплексе при использовании модели 2 описания изомеризации. Для первой модели изомеризации, как следует из выполненного анализа, вполне приемлем (при соответствующем выборе шага итерационной процедуры) более простой метод Эйлера
Особенность моделирования многомерных спектров и решения обратных спектральных задач заключается в необходимости многократно осуществлять расчет при различных параметрах молекулярной модели. В ходе решения задачи о моделировании межизомерных превращений молекул на персональном компьютере было выяснено, что такое моделирование возможно, но для многоатомных молекул и большого числа изомеров время расчета значительно возрастает, что не позволяет производить анализ и оценку результатов с высокой скоростью Поэтому была поставлена и решена задача создания версии программ комплекса для суперкомпьютера (типа МВС-1000). Это потребовало принципиального изменения программ и перехода к параллельным вычислениям. При распараллеливании вычислений принято, что каждый из зарезервированных процессоров проводит вычисления только для одной изомерной формы, закрепленной за ним. Как было показано, такой подход является оптимальным, поскольку физика моделируемого процесса такова, что все участвующие в задаче изомеры имеют почти одинаковое количество активных уровней и, следовательно, все процессоры будут приходить почти одновременно к точке синхронизации (обмен информацией о заселенностях резонансных уровней) независимо от особенностей конкретной физической модели.
Реализован специальный разработанный алгоритм автоматического исключения уравнений с нулевыми решениями, основанный на имитации физического процесса распада возбужденного состояния. Такая процедура крайне актуальна в случае цепочечной изомеризации. Это позволило существенно понизить размерность задач и увеличить быстродействие метода.
Использована специальная форма представления данных в виде структурированных списков, которые по ходу решения задачи заполняются и хранятся в оперативной памяти многопроцессорного вычислителя, минуя длительную процедуру записи во внешнюю память. В результате достигнута существенная экономия вычислительных ресурсов и сняты ограничения на число изомеров в цепи межизомерных преобразований и размеры молекул.
Разработана специальная процедура автоматического определения оптимального шага при численном интегрировании системы кинетических уравнений. Он основан на анализе вероятностей всех переходов, величин частот квантовых биений и выборе шага и времени наблюдений, оптимальных с точки зрения точности и производительности расчета и информативности получаемых
спектральных и кинетических характеристик. Тем самым сняты ограничения. возникающие из-за возможных существенных (на порядки величин) различий во временных характеристиках экспоненциальной релаксации возбужденных состояний изомеров и квантовых биений их состояний.
Алгоритмическим путем преодолена проблема возможного возникновения расходящихся решений для модели 2 (в критических точках, отвечающих нулевым заселенностям уровней-акцепторов) из-за ошибок округления при вычислениях.
Разработанный комплекс программ обеспечивает построение временных зависимостей заселенностей всех вибронных уровней многоизомерных молекулярных систем, полных трехмерных спектров поглощения и флуоресценции каждого изомера и суммарных, а также любых (по времени и частоте) их двумерных сечений. Показано, что в силу сложности решаемых задач, с одной стороны, и оптимально решенной в программном комплексе проблемы распараллетива-ния вычислений с учетом специфики задач, с другой, суперкомпьютер становится инструментом очень высокой эффективности при моделировании фотохимических превращений. Многие из большой совокупности модельных задач, поставленных в процессе отладки методики, алгоритмов и программ и при анализе закономерностей изомеризации, удалось решить именно благодаря использованию суперкомпьютера. Эта версия комплекса программ не имеет каких-либо ограничений по размерности задачи; ограничения накладывает только быстродействие ЭВМ, поэтому выигрыш за счет использования суперкомпьютера весьма велик (см. табл. 1).
Таблица 1. Время полного расчета (мин.) фотохимических превраще-
ний ряда моделей молекул
Тип компьютера Модельные примеры *
1 (2; 100; 100) 2 (4; 800; 1000) 3 (4; 800: 5000)
Стандартный ПК базовой конфигурации 1 15 100
Суперкомпьютер типа МВС-1000 1 1 5
* В скобках приведены характеристики примеров модельных расчетов -
число изомеров; число уравнений; число итераций. Например, на решение системы, состоящей из 5000 уравнений, с использованием суперкомпьютера потребовалось менее 1 часа, в то время как на ПК временные затраты для этой же задачи почти на два порядка больше и составляют около 3-х недель. Весьма существенно, что разработанный программный комплекс позволяет осуществлять моделирование очень сложных процессов изомеризации (для длинных цепочек изомерных превращений, многоуровневых молекулярных систем) в масштабе реального времени. На основе испытаний можно утверждать, что разработаны эффективные средства моделирования межизомерных фотохимических преобразований достаточно сложных соединений, обеспечи-
вающие проведение массовых компьютерных экспериментов в масштабе реального времени. Тем самым можно ставить и решать задачи выбора оптимальных условий натурного эксперимента методами компьютерного моделирования. Разработанный комплекс в совокупности с современными возможностями спектрального эксперимента с фемтосекундным временным разрешением (включая и случай межизомерных перестроек) делает возможным и весьма перспективным применение методов высокоскоростной динамической спектроскопии при качественном и количественном химическом анализе.
Глава 4. Моделирование фотохимических превращений сложных молекул
Разработка комплекса программ позволила перейти к моделированию фотохимических межизомерных превращений достаточно сложных систем. Представляло интерес, в частности, изучение зависимости кинетики межизомерных превращений и ее спектральных проявлений от начальных условий задачи (характер возбуждения, начальная изомерная форма и ее возбужденное состояние и др.); определение условий, при которых возможно наблюдение спектров флуоресценции для образовавшихся изомерных форм, по интенсивности сравнимых со спектром исходного изомера; выявление условий, обеспечивающих проявление эффекта квантовых биений в структуре спектров и его особенностей; исследование временных и количественных характеристик образования изомеров. Особенно интересовали процессы цепочечных изомер-изомерных перестроек, связанные с передачей сигналов и энергии во внутримолекулярном пространстве. Анализировались способы применения развитых методов и реализованных программ к качественному и количественному анализу.
Отдельной задачей является определение критериев выбора наиболее оптимальной с точки зрения эффективности моделирования модели описания межизомерных переходов, а также аппаратных средств решения задач (ПК или суперкомпьютера) .
В качестве примеров, иллюстрирующих эти закономерности, остановимся на модельных расчетах изомер-изомерных преобразований пентадиена (пента-диен-1,3 пентадиен-1,4, соответственно, изомеры 1, 2) и гептадиенилбензола (гепта-1,6-диенилбензол О гепта-1,5-диенилбензол гепта-1,4-диенилбензол гепта-1,3-диенилбензол, соответственно, изомеры 1, 2, 3, 4) при различных начальных условиях и параметрах моделей.
Результаты компьютерных экспериментов показали следующее. Вид графиков временного поведения интенсивностей колебательной структуры виброн-ных спектров (см. рис. 3) во многом зависит от соотношения вероятностей электронно-колебательного перехода и межизомерного (частоты квантовых биений); его изменения весьма значительны даже при относительно небольших изменениях 6). При а)»\*> (ус - вероятность оптического электронно-колебательного перехода) и спектры изомеров подобны и отличаются лишь частотой бие-
ний, которая явно выражена во временной зависимости интенсивностей обоих изомеров. При этом интенсивности спектров изомеров близки по величине.
отличаются менее чем на 1%, в то время как при интенсивность для изомера 2 на 4 порядка меньше, чем для изомера 1. Условия возможности наблюдения спектра изомера 2 на фоне 1-го -Осцилляции ин-тенсивностей линий
неполные (коэффициент модуляции <15 %; рис. 3. 4). что связано с их сложным характером (суперпозиция нескольких). Такая особенность, появляющаяся даже при полном смешивании взаимодействующих уровней при резонансе, очень важна, поскольку этот эффект наблюдается в эксперименте и его интерпретация неоднозначна. В модельных расчетах проявляется смещение во времени максимума спектра люминесценции, что связано с постепенным образованием изомерных форм (см., например, рис. 4).
Графики временных зависимостей населенностей (рис. 5) демонстрирую! общий характер уменьшения их величин для первоначально возбужденных электронно-колебательных уровней энергий изомера 1 и увеличения для основных состояний всех изомеров. Суммарная населенность колебательных уровней основного электронного состояния 1-го изомера не достигает значения первоначального заселения этого изомера из-за пе-
рехода молекул в изомерные формы 2, 3 и 4. Количество молекул в разных изомерных формах, образовавшихся в результате межизомерных переходов, при цепочечной изомеризации существенно зависит от параметров моделей.
Временная зависимость населенностей уровней изомеров носит сложный характер. Отличие суммарных вероятностей переходов с резонирующего уровня двух изомеров на уровни каждого из них приводит к модуляции монотонной
Рис 4 Рассчитанные спектры флуоресценции изомеров 1 (а), 2 (б), 3 (в), 4 (г) гептадие-нилбензола с учетом межизомерных переходов при широкополосном возбуждении изо-меоа 1
временной зависимости заселенности этого уровня в виде хорошо выраженных «ступенек» (см., например, кривые 3 на рис. 5). Если происходит также и заселение этих резонансных уровней за счет переходов молекулярной системы из более высоких по энергии состояний, то модуляция временной зависимости заселенности приобретает осцилляторный характер, причем тем более ярко выраженный (коэффициент модуляции), чем сильнее различия в суммарных вероятностях переходов (см., например, кривые 4, 5 на рис. 5). В целом обнаруживается существенная не только количественная, но и качественная зависимость трехмерных спектров и заселенностей уровней изомеров от условий постановки эксперимента и параметров молекулярных моделей. Сложность получаемой картины натурных спектров однозначно указывает, что интерпретация крайне затруднена (даже невозможна) без проведения компьютерного эксперимента. Его задача состоит, с одной стороны, в определении оптимальных условий постановки натурного эксперимента для выявления и качественной регист-
рации искомых эффектов, а с другой. - в интерпретации натурного эксперимента (в частности, идентификация продуктов фотохимических реакций, оценка времени их накопления, определение квантовых выходов и т д) и опредетении величин и параметров молекулярных моделей, адекватных реальным внутримолекулярным взаимодействиям
Важнейшим является вопрос о цепочечной изомеризации в сложных молекулах. Специфика процесса цепочечной (последовательной) изомеризации изучалась на примере последовательной изомеризации гептадиенилбензола
На рис 6 приведена зависимость конечных засе-ленностей основных состояний изомеров от величины вероятности оптических внутриизомерных переходов Хорошо видно быстрое возрастание доли четвертого изомера, те конечного в цепи и наиболее сильно структурно отличающегося от первого Принципиально важно, что качественные изменения ситуации - доля 4-го изомера увеличивается от пренебрежимо малой до превалирующей - наблюдаются при сравнительно незначительных (в пределах одного порядка) изменениях вероятностей переходов Возникает вполне выраженный эффект, который в химическом эксперименте по внешним признакам (без детального анализа во времени заселенностей возбужденных состояний изомеров цепочки) может быть
Рис 5 Рассчитанные изменения населенностей уровней изомеров 1 (а), 2 (б), 3 (в), 4 (г) гептадиенилбензола с учетом межизомерных переходов при широкополосном возбуждении изомера 1 / - населенность первого возбужденного чисто электронного состояния изомера 1, 2 — основных состояний изомеров, 3-5 — резонансных уровней, соответственно, 1-го и 2-го, 2-го и 3-го, 3-го и 4-го изомеров Населенности нормированы на максимальное значение для соответствующих резонансных уровней
ошибочно интерпретирован как прямой переход через возбужденное состояние изомера 1 в форму 4 с весьма существенной структурной перестройкой Временная характеристика данного процесса изомеризации, как показывает вычислительный эксперимент, мало зависит от длины цепочки изомерных преобразований Процессы заселения возбужденных состояний изомеров 2-4 носят практически синхронный характер Эти состояния промежуточных изомерных форм 2 и 3 выполняют роль "передаточного звена" изомерного преобразования 1 - > 4 поч-
ти без задержки во времени и "потери сигнала", т.е. без образования устойчивых состояний изомеров 2, 3. Такая ситуация характерна для изомерных превращений через триплетные состояния.
Рисунок 7 иллюстрирует «инверсию» количеств 3-го и 4-го изомеров при малых изменениях вероятности оптической релаксации четвертого. Обнаружение этой своеобразной неустойчивости направления хода фотохимической реакции представляется очень важным, так как типичным для фотохимических процессов, особенно сложных, является как раз их мобильность и довольно сильная зависимость от межмолекулярных взаимодействий. Известно, что именно последние часто приводят к существенным перераспределениям интенсивностей в спектрах и, следовательно, изменениям вероятностей оптических переходов. Oi-сюда - большая неопределенность в интерпретации натурных экспериментов. Видно, что компьютерное моделирование может очень помочь в этом. Кроме планирования спектральных экспериментов с временным разрешением и интерпретации их результатов можно также изучать процессы передачи сигналов и энергии во внутримолекулярном пространстве за счет последовательных структурных перестроек и процессы, связанные с логической переработкой информации, лежащие в основе нового научного направления, иногда называемого молекулярной логикой.
Выше отмечалось, что для моделирования фотохимических межизомерных преобразований были предложены две модели изомер-изомерных преобразований. Проведенные модельные расчеты для сложных многоуровневых, многоизомерных систем показали их полную эквивалентность. Различия в получаемых кинетических и спектральных кривых пренебрежительно малы, менее 1%. Они обусловлены ошибками округления при компьютерных вычислениях. Выбор наиболее- приемлемого подхода при молекулярном моделировании определяется личными пристрастиями исследователя, наглядностью и удобством описания молекулярной модели и интерпретации результатов. Следует, однако, заметить, что расчеты с использованием модели 2 требуют несколько больших временных затрат, чем для модели 1. В тоже время, модель 2. возможно, более предпочтительна с точки зрения наглядности и привычности описания молекулярных моделей.
п, ото ед
5 10 15 20 23 30 38 40 45 50
И-, X Ю® с'1
Рис. 6. Конечные заселенности основных состояний изомеров I (1), 2 (2), 3 (3) и 4 (4) гептадиенилбензола при компьютерном эксперименте, отвечающем варьированию вероятностей оптических переходов.
Обратим внимание на общие закономерности, обнаруженные в эксперименте:
1. при очень коротких импульсах возбуждения наблюдается периодические изменения на фоне общего временного поведения спектральных линий
2. Временные зависимости отдельных линий могут сильно различаться даже качественно.
3. При развитии спектральной картины наблюдаются запаздывания, связанные с фотоизомерией.
4. Накопление изомерных форм при фотоизомеризации очень чувствительно к особенностям структуры объекта и условиям проведения эксперимента
Все эти закономерности и другие эффекты чётко воспроизводятся в компьютерных экспериментах. Это позволяет утверждать, что развитые методы моделирования адекватно описывают реальные природные процессы.
Остановимся, на возможных аналитических приложениях динамической спектроскопии изомерных превращений молекул и, соответственно, разработанных средств их компьютерного моделирования.
Прежде всего, временная спектроскопия сложных молекул высокого разрешения может быть использована как эффективный инструмент при количественном анализе молекул и их смесей. Это было ранее показано для случая отсутствия межизомерных переходов при фотовозбуждении молекул [1] Фотохимические превращения существенно усложняют динамику спектров, появляются дополнительные параметры, характеризующие вероятности (и динамику) межизомерных переходов. Разработанный ранее [1] метод анализа применим в полной мере, поскольку эти дополнительные параметры молекулярных моделей (частоты квантовых биений) легко определяются по осцилляторной модуляции интенсивностей линий в экспериментальных спектрах. Единственное требование к эксперименту -это выбор таких условий его проведения, чтобы биения проявлялись в спектрах. Таким образом, в данной работе сделан следующий, очень важный шаг в теории молекулярного качественного и количественного спектрального анализа - по-
казаны возможности распространения метода на случай структурных фотохимических превращений.
Динамическая спектроскопия фотохимических преобразований молекул и разработанные средства их компьютерного моделирования открывают новые, существенные и очень важные возможности. Появляется новый спектральный сигнал для химического анализа. Таковым может служить, во-первых, сам факт появления другой изомерной формы молекулы и соответствующего спектра и, во-вторых, частота квантовых биений и соответствующая осцилляторная модуляция интенсивностей спектральных линий. Эта вторая характеристика уже используется в литературе для анализа. Разработанные средства моделирования позволяют провести этот анализ на более корректный и количественный уровень. Надо заметить, что экспериментальные методы фемтосекундной спектроскопии и были развиты для наблюдения процессов фотохимических превращений, химических реакций при фотовозбуждении. Разработанные средства компьютерного моделирования восполняют пробел в теоретической части этой области и создают теоретическую базу количественных методов анализа в фотохимии.
Другое важное аналитическое приложение разработанных теоретических методов - это изучение процессов передачи сигнала во внутримолекулярном пространстве, например, при цепочечной изомеризации, рассмотренной выше. Как следствие, - возможность создания аналитических реагентов нового типа с разделением области реакции и области образования аналитического сигнала с существенным повышением чувствительности и избирательности [4]. Для поиска молекулярных систем с такими свойствами необходим целенаправленный компьютерный эксперимент. Разработанные средства моделирования как раз и являются базовыми для постановки задач таких задач.
Проведённые исследования показали перспективы и реальные возможности применения новейших методов фемтосекундной спектроскопии в физической и аналитической химии.
ВЫВОДЫ
1. Разработан метод моделирования фотохимических изомерных превращений сложных молекул при произвольном по спектральной форме возбуждающем импульсе света и для двух модельных подходов для описания межизомерных переходов. Показано, что эти модели изомеризации полностью эквивалентны. На этой основе предложена новая система кинетических уравнений, полностью описывающая развитие процессов во времени.
2. На базе развитых методов решения систем кинетических уравнений очень высокого порядка с зависящими от времени коэффициентами, разработанных алгоритмов оптимизации и распараллеливания вычислений, создан программный комплекс сервисного типа, реализующий разработанный метод в двух версиях - на ПК и суперкомпьютере. Комплекс программ не имеет принципиальных ограничений на характеристики молекулярных моделей и позволяет осуществ-
лять компьютерное моделирование фотохимических превращений сложных многоизомерных систем в масштабе реального времени
3 Изучены общие закономерности изомерных преобразований и их спектральных проявлений
- зависимость кинетики межмолекулярных превращений и их спектральных проявлений от особенностей возбуждения и величин параметров молекулярной системы,
- условия, обеспечивающие возможность регистрации спектров изомеров, образовавшихся в результате фотохимического превращения,
- временные и количественные характеристики образования изомеров
- специфические особенности цепочечных, последовательных изомерных перестроек и, как следствие, передачи сигналов и энергии во внутримолекулярном пространстве В частности, показана своеобразная «неустойчивость» направления хода фотохимических процессов, высокая чувствительность к условиям эксперимента, роль промежуточных изомерных форм как «передаточного звена» в цепочке преобразований без задержки во времени и «потери сигнала»
4 Разработанные средства молекулярного моделирования впервые обеспечивают возможность целенаправленного поиска оптимальных условий постановки натурных экспериментов в области фотохимических изомерных превращений и позволяют проводить детальные теоретические исследования таких процессов В целом, создана теоретическая база количественных методов анализа в фотохимии
Материалы диссертации изложены в следующих работах:
1. Баранов В.И., Завалий М.В., Грибов Л.А. Моделирование и расчет динамических спектров с учетом изомеризации сложных молекул // Журн. прикл. спектр. 2003. Т. 70. С. 626-634.
2. Баранов В.И., Завалий М.В., Грибов Л.А. Моделирование процессов изомеризации и спектров с временным разрешением сложных молекул // Журн. прикл. спектр. 2004. Т. 71, № 3. С. 295-301.
3. Грибов Л.А., Баранов В.И., Завалий М.В. Метод моделирования динамических вибронных спектров и состояний с учетом межизомерных переходов сложных молекул // ДАН. 2004. Т. 397. № 6. С. 447-450.
4. Завалий М.В. Моделирование динамических вибронных спектров сложных молекул с учетом межизомерных переходов // 3-я Всероссийская конференция «Молекулярное моделирование», Тезисы докладов, 2003, М., С. 67.
5. Баранов В.И., Завалий М.В. Программный комплекс для моделирования процессов изомеризации сложных молекул на суперкомпьютере типа МВС-1000 // Конференция "Информационно-вычислительные технологии в фундаментальных и прикладных физико-математических исследованиях", Сборник .материалов, 2004, М., С. 50.
Цитированная литература:
1. Астахов С.А., Баранов В.И., Грибов Л.А. //Журн. аналитич. химии. 2001. Т. 56. №7. С. 703-713.
2. Грибов Л.А. От теории спектров к теории химических превращений. - М.: Эдиториал УРСС, 2001. - 368 с.
3. Baranov V.I., Gribov LA, DjenjerV.O.. Zelent'sov D.Yu. // J. Mol. Struct. 1997. V. 407. No 2-3. P. 177-216.
4. Грибов Л.А., Сударушкин С.К. // Журн. аналитич. химии. 2004. Т. 59. № 4. С.343-348.
Отпечатано на ризографе в ОНТИ ГЕОХИ РАН Заказ Лг5> 2-й Тираж 120 экз
02.00
1 1 - Г uJ
Введение.
Глава 1. Обзор литературы.
1.1 Развитие экспериментальных и теоретических методов электронно-колебательной спектроскопии.
1.2 Состояние теории качественного и количественного анализа без 19 использования стандартных образцов состава.
1.3 Выводы к главе 1.
Глава 2. Развитие теории моделирования фотохимических межизомерных преобразований с учетом квантовых биений для сложных молекул.
2.1 Теоретическая модель для описания изомер-изомерных превращений.
2.2 Система кинетических уравнений для населенностей уровней.
2.3 Модель описания изомерных преобразований как безызлучательных переходов между состояниями изомеров.
2.4 Выводы к главе 2.
Глава 3. Создание средств моделирования фотохимических превращений сложных молекул.
3.1 Общая структура метода моделирования и программного комплекса.
3.2 Способ решения системы кинетических уравнений.
3.3 Программная реализация методов расчета изомеризации молекул на персональном компьютере.
3.4 Программная реализация методов расчета изомеризации молекул на суперкомпьютере типа МВС-1000.
3.5 Особенности программного решения системы кинетических уравнений для населенностей уровней для модели 2 изомер-изомерных преобразований.
3.6 Выводы к главе 3.
Глава 4. Моделирование фотохимических превращений сложных молекул.
4.1 Моделирование кинетики внутримолекулярных превращений и спектров сложных систем с учетом межизомерных переходов.
4.2 Моделирование кинетики и спектральных проявлений процесса цепочечной изомеризации.
4.3 Сравнение моделей описания межизомерных преобразований.
4.4 О возможных приложениях развитых методов.
4.5 Выводы к главе 4.
Аналитическая химия в полном объеме уже взяла на вооружение спектральные методы исследования вещества. Нет необходимости в перечислении всех методов спектрального анализа, которые нашли широкое применение в аналитической химии. Их число велико.
В настоящие время в экспериментальных методах молекулярной электронно-колебательной спектроскопии достигнут существенный прогресс. Стало возможным получение спектров с высоким разрешением колебательной структуры при помощи методики сверхзвуковой струи, наблюдения флуоресценции с одиночных возбужденных электронно-колебательных уровней и динамических (трехмерных) спектров с временным разрешением при возбуждении короткими пико-, фемтосекундными импульсами света. Все это предоставило в руки химиков-аналитиков новые возможности.
Методы, основанные на получении и исследовании трехмерных спектров, давно заняли свое важное место в аналитической химии. Эти методы включают в себя синхронные спектры, получаемые при одновременном изменении частот возбуждающей линии и наблюдении спектра люминесценции, а также спектры, где в качестве дополнительного измерения учитывается время (техника импульсного фотолиза (флеш-фотолиза) и др.). Еще в первых работах с использованием техники флеш-фотолиза обнаружилось, что использование спектров с временным разрешением предоставляет новые возможности для анализа, т. к. люминесценция различного рода компонентов, входящих в смесь веществ, происходит с некоторой «разверткой» по времени. Поэтому можно указать такие временные диапазоны, в пределах которых выгоднее изучать спектр либо одного, либо другого компонента.
В последние 20 лет лазерная спектроскопия, благодаря открытиям в области техники коротких (фемтосекундных) лазерных импульсов, дала в руки химикам качественно новый инструмент - импульсы света фемтосекундной длительности. Они принесли в химию новые возможности и воплотились в фемтохимии, стали толчком для развития нового направления науки - фем-тосекундной спектроскопии, основанной на резонансном возбуждении отдельных высоколежащих уровней энергии молекул с последующим наблюдением их высвечивания. Это позволило получать очень развитые спектры с временным разрешением, в которых можно наблюдать как изменения интен-сивностей излучения или поглощения во времени на определенных длинах волн, так и сразу весь спектр в определенные интервалы времени.
Фемтосекундная технология позволяет изучать химические реакции с высоким временным и пространственным разрешением. Стало возможным наблюдение динамики реакций и целенаправленное управление их ходом с высоким временным разрешением, что открывает значительные перспективы в области качественного и количественного анализа, синтеза новых химических соединений и предсказания их реакционной способности. Импульсы фемтосекундной длительности дают мощный инструмент для исследования элементарных превращений на уровне процессов разрыва и образования связей, переноса заряженных частиц в мономолекулярных (например, фотодиссоциация, изомеризация) и бимолекулярных реакциях. Фемтосекундная спектроскопия позволяет регистрировать временные зависимости интенсив-ностей излучения или поглощения света молекулой в ходе протекания внутримолекулярных и межмолекулярных процессов (в том числе и безызлуча-тельных), что дает возможность изучать изомер-изомерные преобразования.
Все перечисленные выше возможности использования импульсов фемтосекундной длительности на сегодняшний день еще до конца не реализованы. Прежде всего, из-за сложностей экспериментальной техники и того, что сама постановка эксперимента требует предварительно выбора оптимальных условий, так как чисто эмпирический поиск наиболее выгодной с точки зрения окончательного результата постановки эксперимента, является очень громоздким.
Спектральные методы как качественного, так и количественного анализа, применяемые в аналитической химии, основаны на использовании эталонов, или стандартных образцов состава. Стандартные образцы состава позволяют решать качественные и количественные задачи, используя разнообразные методы оптического наблюдения - от эмиссии отдельных атомов до спектров поглощения сложных соединений. Необходимость использования эталонов для получения окончательного ответа составляет очень большое ограничение в применении спектроскопии как таковой. Вопрос о том, как можно обойти использования эталонов, возник не так давно, воплотившись в теорию качественного и количественного анализа без использования натурных стандартных образцов состава. Такой анализ может быть проведен, если воспользоваться теоретическими расчетными спектрами разного рода, в том числе и динамическими спектрами, получаемыми с использованием фемто-секундных импульсов.
Применение техники фемтосекундной спектроскопии в аналитической и физической химии является в некоторой степени новым направлением в науке. Но как уже упоминалось, использование ее является нетривиальной и мало изученной при решении задач аналитической и физической химии. Высока сложность интерпретации динамических спектров, которые имеют не только тонкую структуру вдоль шкалы частот, но и сложным образом зависят от времени. При этом вероятны случаи изменения интенсивности излучения или поглощения света молекулой в ходе протекания внутримолекулярных и межмолекулярных процессов. Для выбора подходящих условий крайне желательно было бы воспользоваться возможностями компьютерного эксперимента, так как совершенно ясно, что с помощью компьютерных экспериментов можно заранее, варьируя заданную ситуацию, выбрать наиболее оптимальный путь постановки экспериментов, после чего уже реализовывать их в натурном виде.
Попыток проведения моделирования динамики фотохимических изомер-изомерных молекулярных преобразований, возникающих в результате облучения молекулы короткими фемтосекундными лазерными импульсами, для реальных моделей сложных молекул никогда раньше не предпринималось. Между тем есть все основания для постановки такой задачи. Успех работ по расчету динамических спектров для одиночных молекул без возможных их фотохимических преобразований и применению этих спектров к вопросам качественного и количественного анализа в аналитической химии, а также теоретических работ по описанию межизомерных превращений, подтверждает это.
Следует отметить также, что изомерные преобразования молекул служат основой для передачи энергии и информации, сигнала во внутримолекулярном пространстве. Как было показано, на этом эффекте могут быть построены реагенты нового типа с пространственным разделением в молекуле реакционной области и области появления аналитического сигнала. Как следствие, существенно повышается их аналитическая чувствительность, избирательность и другие характеристики. Теории таких внутримолекулярных процессов не было, а исследовать их можно лишь методами динамической спектроскопии. Развитие такой теории весьма актуально, так как позволяет выйти в совершенно новую область науки - молекулярную информатику.
Из сказанного выше видно, что моделирование динамики фотохимических процессов является весьма перспективным для аналитической и физической химии. Для обеспечения такого моделирования необходимо решить следующие задачи.
1. Развитие теории и разработка методов компьютерного моделирования фотохимических изомерных преобразований и квантовых биений состояний сложных молекул на базе детального описания их моделей, обеспечивающих проведение предсказательных расчетов в масштабе реального времени.
2. Проведение на основе разработанных методов модельных расчетов и анализ получаемой спектральной картины и кинетики образования продуктов фотохимической реакции с целью выявления общих закономерностей спекхимии, опережающих натурный эксперимент и позволяющих осуществлять его целенаправленное планирование. На защиту выносятся:
1. Формулировка и решение задачи о моделировании фотохимических межизомерных превращений молекул, возникающих в результате произвольного возмущения, действующего на молекулу и приводящего к начальному заселению ее возбужденных электронно-колебательных уровней энергии.
2. Вычислительные методы, алгоритмы и программы для персонального компьютера и суперкомпьютера, обеспечивающие проведение моделирования рассматриваемых процессов в масштабе реального времени.
3. Результаты конкретных расчетов для изомер-изомерных преобразований сложных молекул (пентадиен, гептадиенилбензол) с демонстрацией и анализом динамики населенностей уровней и динамических спектров.
4. Прогноз возможности расширения методов спектрального анализа и распространения их в область фотохимических превращений вещества.
Все содержание диссертации отличается достаточной степенью новизны, так как все рассмотренные проблемы, в ходе диссертационного исследования, ранее в литературе не ставились и не решались. Решена важная и новая задача для аналитической и физической химии, показывающая перспективность такого моделирования для решения качественных и количественных задач анализа поведения и физической природы передачи сигналов и информации в молекулярном мире сложных молекулярных структур при возбуждении.
Созданные методы расчета и комплексы программ позволяют целенаправленно выбирать наиболее оптимальные условия для постановки натурных экспериментов еще до начала их проведения.
Материалы диссертации доложены на 3-й Всероссийской конференции «Молекулярное моделирование» (г. Москва, 2003 г.) и на Конференции «Информационно-вычислительные технологии в химии, биологии, фармацевтике и медицине» (г. Москва, 2004 г.).
Результаты исследований опубликованы в трех научных статьях и двух тезисах докладов.
Диссертация состоит из введения, четырех глав, выводов, списка литературы и приложений. Общий объем работы составляет 239 страниц текста, включая 29 рисунков, 3 таблицы, список литературы, содержащий 136 наименований, приложений А и Б, объем которых составляет 92 страницы.
ОСНОВНЫЕ ВЫВОДЫ
1. Разработан метод моделирования фотохимических изомерных превращений сложных молекул при произвольном по спектральной форме возбуждающем импульсе света и для двух модельных подходов для описания межизомерных переходов. Показано, что эти модели изомеризации полностью эквивалентны. На этой основе предложена новая система кинетических уравнений, полностью описывающая развитие процессов во времени.
2. На базе развитых методов решения систем кинетических уравнений очень высокого порядка с зависящими от времени коэффициентами, разработанных алгоритмов оптимизации и распараллеливания вычислений, создан программный комплекс сервисного типа, реализующий разработанный метод в двух версиях - на ПК и суперкомпьютере. Комплекс программ не имеет принципиальных ограничений на характеристики молекулярных моделей и позволяет осуществлять компьютерное моделирование фотохимических превращений сложных многоизомерных систем в масштабе реального времени.
3. Изучены общие закономерности изомерных преобразований и их спектральных проявлений:
- зависимость кинетики межмолекулярных превращений и их спектральных проявлений от особенностей возбуждения и величин параметров молекулярной системы;
- условия, обеспечивающие возможность регистрации спектров изомеров, образовавшихся в результате фотохимического превращения;
- временные и количественные характеристики образования изомеров;
- специфические особенности цепочечных, последовательных изомерных перестроек и, как следствие, передачи сигналов и энергии во внутримолекулярном пространстве. В частности, показана своеобразная «неустойчивость» направления хода фотохимических процессов, высокая чувствительность к условиям эксперимента, роль промежуточных изомерных форм как «передаточного звена» в цепочке преобразований без задержки во времени и «потери сигнала».
4. Разработанные средства молекулярного моделирования впервые обеспечивают возможность целенаправленного поиска оптимальных условий постановки натурных экспериментов в области фотохимических изомерных превращений и позволяют проводить детальные теоретические исследования таких процессов. В целом, создана теоретическая база количественных методов анализа в фотохимии.
1. Шпольский Э.В. Линейные спектры флуоресценции органических соединений и их применение// Успехи физ.наук. - 1960. - Т.71. - С.215-242.
2. Шпольский Э.В. Новые данные о природе квазилинейных спектров органических соединений // Успехи физ.наук. 1963. - Т.80. - С.255-279.
3. Никитина А.Н., Пономарева Н.А., Тер-Саркисян Г.С. и др. Спектры Шпольского соединений класса линейных полиенов // Опт. и спектроск. -1976. Т.40. - №2. - С.251-257.
4. Durand Y., Bloeb A., A.M. van Oijen, Kohler U., E.J.J.Groenen, Schmidt J. An optical study of a single pentacene molecules in n-tetradecane // Chem. Phys. Lett. -2000.- V.317. -P.232-237.
5. Walla P.J., Jelezko P., Tamarat Ph., Lounis В., Orrit M. Perylene in biphenyl and anthracene crystals: an example of the influence of the host on single-molecule signals // Chem. Phys. 1998. - V.233. - P.l 17-125.
6. Amirav A., Even U., Jortner J. Analytical applications of supersonic jet spectroscopy // Anal. Chem. 1982. - V. 54. №11. - P. 1666-1673.
7. Beck S.M., Powers D.E., Hopkins J.B., Smally R.E. Jet-cooled naphthalene. II. Single vibronic level fluorescence spectra // Chem. Phys. 1981. - V. 74. - № 1. -P.43-52.
8. Современные аспекты тонкоструктурной и селективной спектроскопии: Межвузовский сб. научн. тр // Моск. педагогич. ин-т. М.:МГПИ, 1984. - 204 с.
9. Syage J.A., Felker P.M., Semmes D.H., Adel F.A1., Zewail A.H. Picosecond excitation and selective intramolecular rates in supersonic molecular beams. IV. Alkylanthracenes // J. Chem. Phys. 1985. - V.82. -P.2896-2908.
10. Kaimoto O., Hayami S., Shizuka H. The LIF spectra and the structure of jet-cooled ADMA // Chem. Phys. Lett. 1991. - V. 177. - P.219-222.
11. Cheatham C.M., Fluang M., Meinander N. Kelly M.B., HalerK., Chiang W.Y., Laane J. A computer-controlled apparatus for laser induced fluorescence spectroscopy in a supersonic jet // J. Mol. Struct. 1996. - V.377. - P.81-89.
12. Mukamel S. Direct vs redistributed fluorescence of supercooled molecules: the onset of intramolecular vibrational redistribution // J. Chem. Phys. 1985. - V.82. P.2867- 2876.
13. Patwari G.N., Doraiswami S., Wategaonkar S. IVR in the Si state of jet-cooled cis- and trans-p-dimethoxybenzene // Chem. Phys. Lett. 2000. - V.316. - P.433-441.
14. Patwari G.N., Doraiswamy S., Wategaonkar S. Electronic relaxation of the Si state of jet-cooled p-alkoxyanilines // Chem. Phys. Lett. 1999. - V.305. - P.381-388.
15. Syage J.A., Felker P.M., Zewail A.H. Picosecond dynamics and photoisomerization of stilbene in supersonic beams. I-II. // J. Chem. Phys. 1984. -V.81. - P.4685-4723.
16. Prochorow J., Deperasinska I., Morawski O. Fluorescence excitation and fluorescence spectra of jet-cooled acridine molecules: acridine dimmer formation and structure // Chem. Phys. Lett. 2000. V.316. - P.24-30.
17. Harms G.S., Irngartinger Т., Reiss D., Renn A., Wild U.P. Fluorescence lifetimes of terrylene in solid matrices // Chem. Phys. Lett. 1999. V.313. - P.533-538.
18. Грибов JI.А., Баранов В.И., Зеленцов Д.Ю. Электронно-колебательные спектры многоатомных молекул: Теория и методы расчета. М.: Наука, 1997. - 475 с.
19. Грибов Л.А., Баранов В.И., Новосадов Б.К. Методы расчета электронно-колебательных спектров многоатомных молекул. М.: Наука, 1994. - 325 с.
20. Kinsey J.L. Laser-induced fluorescence // Ann. Rev. Phys. Chem. 1977. -V.28. - P. 349-372.
21. Лазерная и когерентная спектроскопия. Под ред. Дж. Стейнфелда. -М.:Мир, 1982.-629 с.
22. Демтредер В. Лазерная спектроскопия. Основные принципы и техника эксперимента: Пер. с англ. / Под ред. И.И. Собельмана. М.:Наука. Гл. ред. физ.-мат. лит, 1985. - 608 с.
23. Heikal А.А., Baskin J.S., Banares L., Zewail A.H. Structural effects of the isomerization dynamics of trans-stilbenes: IVR, microcanonical reaction rates, and the nature of the transition state // J. Phys. Chem. A. 1997. - V. 101. - P.572-590.
24. Tsuji K., Terauchi C., Shibuya K., Tsuchiya S. Trans-cis isomerization of thea}~ A state as studied by dispersed fluorescence spectroscopy // Chem. Phys. Lett. 1999.-V. 306. -P.41-47.
25. Patwari G.N., Wategaonkar S. Does inversion symmetry influence IVR. // Chem. Phys. Lett. 2000. - V. 323. - P.460-466.
26. Weiner A.M. Femtosecond pulse processing // Opt. Quant. Electr. 2000. -V.32. - P.473-487.
27. Zewail A.H. Laser femtochemistry // Science. 1988. - V.242. - P. 1645.
28. Sundstrom Ed.V. Femtochemistry and Femtobiology. World Scientific. Singapore. 1998.
29. Zewail A.H. Femtochemistry // Adv. in Chem. Phys. 1997. - V. 101. - P.892.
30. Manz J, Woste L. (Eds.) Femtosecond chemistry. VCH. New York. 1995.
31. Zewail A.H. Femtochemistry ultrafast dynamics of the chemical bond. World Scientific. Singapore. 1994.
32. Бучаченко A.JI. Химия на рубеже веков: свершения и прогнозы. // Успехи химии 1999 - Т. 68 - Вып. 10 - С. 99 - 118.
33. Mukamel S. Femtosecond optical spectroscopy: a direct look at elementary chemical events //Ann. Rev. Phys. Chem. 1990. - V.41. - P.647-681.
34. Takayoshi Kabayashi, Akira Shirakawa. Sub-10-fs tunable pulses in visible and NIR and visible sub-5-fs pulses generated by noncollinear OPA // J. Lumin. 2000. -V.87-89.-P.119-120.
35. Baskin J.S., Bañares L., Pedersen S., Zewail A.H. Femtosecond real time probing of reactions. 20. Dynamics of twisting, alignment, and IVR in trans-stilbene isomerization reaction // J. Phys. Chem. 1996. - V.100. - P. 11920-11933.
36. Assion A., Baumert Т., Bergt M., Brixner Т., Kiefer В., Seyfried V., Strehle M. and Gerber G. Control of chemical reactions by feedback-optimized phase-shaped femtosecond laser pulses // Science. 1998. - V.282.- P.919-922.
37. Kasaai M.R., Arul J., Chim S.L., Charlet G. The use of intense femtosecond laser pulses for fragmentation of chitosan // J. Photochem. Photobiol. A: Chem. -1999.-V.120.-P.201-205.
38. Чибисов A.K. Применение импульсного фотовозбуждения для исследования трипетных состояний органических веществ // Успехи химии. -1970. Т. 39. - Вып. 10. - С.1886-1911.
39. Чибисов А.К. Кинетическая спектроскопия и ее применение в аналитической химии // Ж. аналитич. химии. 1983. - Т.38. - Вып. 6. - С.1104-1130.
40. Atabekyan L.S., Chibisov А.К. Acid effect on photochromism of spiropyrans: a study by microsecond and nanosecond flash photolysis // J.Photochem. 1985. -V.31, - P.123-130.
41. Атабекян Jl.С., Чибисов А.К. Лазерная кинетическая спектроскопия фотохромных спиропиранов и их комплексов. Итоги науки и техники. Современные проблемы лазерной физики. М.: ВИНИТИ, 1990. - Т.1. - С. 152-196.
42. Чибисов А.К., Захарова Г.В., Шаповалов В.Л., Толмачев А.И., Брике Ю.Л., Сломинский Ю.Л. Импульсный фотолиз полиметиновых красителей с различными противоионами в неполярных средах // Химия высоких энергий.- 1995. Т.29. - № 3. - С.211-218.
43. Chibisov А.К., Zakharova G.V., Shapiro B.I. Flash-photolysis study of dye positive holes in spectrally sensitized emulsions // J.Information Record. Materials.- 1996. -V.22. P. 313-326.
44. Chibisov A.K., Gorner H. Photoprocesses in spiropyran-derived merocyanines // J.Phys. Chem.A. 1997. - V. 101. - № 24. - P.4305-4312.
45. Khimenko V., Chibisov A.K., Gorner H. Effects of alkyl substituents in the polymethine chain on the photoprocesses in thiacarbocyanine dyes // J.Phys. Chem.A. 1997. - V. 101. - № 39. - P.7304-7310.
46. Атабекян JI.C., Чибисов A.K., Алфимов M.B. Импульсный фотолиз краунсодержащих стириловых красителей и их комплексов с ионами металлов // Химия высоких энергий. 1997. - Т.31. - № 5. - С.381-385.
47. Chibisov A.K., Gorner H. Complexes of spiropyran-derived merocyanines with metal ions: relaxation kinetics, photochemistry and solvent effects II Chem. Phys. -1998.-V. 237. P.425-442.
48. Gorner H., Chibisov A.K. Complexes of spiropyran-derived merocyanines with metal ions. Thermally activated and light-induced processes // J.Chem.Soc. Faraday Trans. 1998. - V.94. - № 17. - P.2557-2564.
49. Chibisov A.K., Gorner H. Photoprocesses in spirooxazines and their merocyanines // J.Phys. Chem. A. 1999. - V. 103. - № 26. - P.5211-5216.
50. Chibisov A.K., Gorner H. Complexes of spiropyran-derived merocyanines with metal ions: relaxation kinetics, photochemistry and solvent effects // Chem. Phys. -1998.-V. 237.-P.425-442.
51. Felker P.M., Zewail A.H.//Chem. Phys. Lett. 1983.-V. 109. - P. 113-119.
52. Грибов Л.А., Баранов В.И. О сопоставлении экспериментальных и вычисленных оптических молекулярных спектров и о постановке обобщенной обратной спектральной задачи // Опт. и спектроск. 1998. - Т.85. - №1. - С.45-52.
53. Грибов JI.A. Новая постановка обратной спектральной задачи в теории колебательных спектров многоатомных молекул. // Опт. и спектроск. 1997. -Т.82. - №1. - С.23-25.
54. Баранов В.И., Грибов JI.A., В.О Дженжер. Параметрический полуэмпирический метод в теории электронно-колебательных спектров многоатомных молекул. // Журн. структ. химии 1996. - Т.37. - №3. - С.419-431.
55. Баранов В.И., Грибов JI.A., Дженжер В.О., Зеленцов Д.Ю. Параметрический метод расчета электронно-колебательных спектров сложных молекул. Дефинил пол иене. // Журн. структ. химии 1996. - Т.37. -№6. - С. 1041-1049.
56. Baranov V.I., Zelent'sov D.Yu. Methods for calculation of electronic-vibrational spectra of polyatomic molecules // J.Mol.Struct. 1994. - V.328. -P. 179-199.
57. Baranov V. I., Gribov L. A., Zelent'sov D. Yu. Methods of calculation of electronic-vibrational spectra of polyatomic molecules. Part 2. Vibrational solution of the electronic-vibrational problem. J. Mol. Struct., 1994. - V.328. - P. 189-197.
58. Baranov V. I., Zelent'sov D. Yu. J.Mol.Struct., 1994. - V.328. - P. 199-210.
59. Baranov V.I., Gribov L.A., Zelent'sov D.Yu. Refined approach for matrix element calculation in the theory of electronic-vibrational spectra of polyatomic molecules // J.Mol.Struct. 1996. - V.376. - P.475-493.
60. Baranov V.I., Gribov L.A., DjenjerV.O., Zelent'sov D.Yu. Adiabatic semi-empirical parametric method for computing electronic-vibrational spectra of complex molecules // J. Mol. Struct. 1997. - V.407. - №2-3. - P.177-216.
61. Wurzer A.J., Wilhelm Т., Piel J., Riedle E. Comprehensive measurement of the Si azulene relaxation dynamics and observation of vibrational wavepacket motion // Chem. Phys. Lett. 1999. - V.299. - P.296-302.
62. Varnavsky O., Goodson T. Femtosecond fluorescence dynamics and molecular interactions in water-soluble nonlinear optical polymeric dye // Chem. Phys. Lett. -2000. V.320. - P.688-696.
63. Paa W., Yang J.-P., Helbig M., Hein J., Rentsch S. Femtosecond time-resolved measurements terthiophene: fast singlet-triplet intersystem crossing // Chem. Phys. Lett. 1998. - V.292. - P.607-614.
64. Клищенко А.П., Сенюк М.А., Тихомиров С.А., Толсторожев Г.Б. Пикосекундная динамика внутримолекулярного переноса заряда и аномальная люминесценция аминозамещенных оксадиазолов и кетонов // Опт. и спектроск. 1997. - Т.83. - №4. - С.664-668.
65. Бушук Б.А., Рубинов А.Н. Влияние специфических межмолекулярных взаимодействий на динамику спектра флуоресценции растворов красителей // Опт. и спектроск. 1997. - Т.82. - №4. - С.573-576.
66. M.J. van der Meer, Zhang H., Rettig W., Glasbeek M. Femto- and picosecond fluorescence studies of solvation and non-radiative deactivation of ionic styryl dyes in liquid solution // Chem.Phys.Lett. 2000. - V.320. - P.673-680.
67. Jurczok M., Plaza P., Martin M.M., Meyer Y.H., Rettig W. Excited state relaxation paths in 9,9'-bianthryl and 9-carbazolyl-anthracene: a sub ps transient absorption study // Chem. Phys. 2000. - V.253. - P.339-349.
68. Azuma J., Tamai N., Shishido A., Ikeda T. Femtosecond dynamics and stimulated emission from S2 state of a liquid crystalline trans-azobenzene // Chem. Phys. Lett. 1998. - V.288. - P.77-82.
69. Yang J.-P., Paa W., Rentsch S. Femtosecond investigations of photophysics of ultrafast intersystem crossing in terthiophen by wavelength dependent excitation // Chem. Phys. Lett. 2000. - V.320. - P.665-672.
70. Braatz H., Hecht S., Seifert H., Helm S., Bendig J., Rettig W. Photochemistry and photophysics of donor-acceptor-polyenes. I: all-trans-4-dimethylamino-4'-cyano-1,4-diphenylbutadiene (DCB) // J. Photochem. Photobiol. A: Chem. 1999. -V.123.-P. 99-108.
71. Hirata Y., Okada Т., Nomoto T. Photophysical properties of diphenylacetylene derivatives in the solution phase: comparison with the dynamical properties of diphenylbutadiyne // Chem. Phys. Lett. 1998. - V.293. - P.371-377.
72. Morlet-Savary F., Ley C., Jacques P., Wielder F., Fouassier J.P. Time-dependent solvent effects on the Ti-Tn absorption spectra of thioxanthone: a picosecond investigation // J. Photochem. Photobiol. A: Chem. 1999. - V.126. -P.7-14.
73. Ohta K., Kang T.J., Tominaga K., Yoshihara K. Ultrafast relaxation processes from a higher excited electronic state of a dye molecule in solution: a femtosecond time-resolved fluorescence study // Chem. Phys. 1999. - V.242. - P. 103-114.
74. Nishimura Y., Yasuda A., Speiser S., Yamazaki I. Time-resolved analysis of intramolecular electronic energy transfer in methylene-linked naphthalene-anthracene compounds // Chem. Phys. Lett. 2000. - V. 323. - P. 117-124.
75. Астахов C.A., Баранов В.И. Моделирование тонкоструктурных электронно-колебательных спектров многоатомных с временным разрешением. Постановка задачи и анализ кинетических уравнений // Опт. и спектроск.- 2000. Т.90. - №2. - С.237-245.
76. Astakhov S.A., Baranov V.I., Gribov L.A. Advances in the theory methods of computional vibronic spectroscopy.// J. Mol. Struct. 2003. - V. 655. - № 1. - P. 97-123.
77. Баранов В.И., Соловьев А.Н. Параметрический метод в теории вибронных спектров. Структура стильбена в возбужденном состоянии и спектр флуоресценции // Журн. структ. химии 2000. - Т.41. - №2. - С.370-378.
78. Баранов В.И. Расчет параметрическим методом электронно-колебательных спектров и структуры метилстирола в возбужденном состоянии // Журн. структ. химии 1999. - Т.40. - №2. - С.242-250.
79. Баранов В.И. Параметрический метод в теории вибронных спектров сложных молекул. Спектры поглощения и флуоресценции и структура стирола в возбужденном состоянии // Опт. и спектроск. 2000. - Т.88. - №2. -С.216-223.
80. Баранов В.И., Соловьев А.Н. Параметрический метод расчета возбужденных состояний и вибронных спектров сложных молекул. Спектры поглощения и флуоресценции перилена // Журн. прикл. спектр. 2000. - Т.67. - №2. - С.148-153.
81. Грибов Л.А., Баранов В.И. По поводу одного метода определения потенциальных поверхностей и соотношения Душинского в теории электронно-колебательных спектров молекул // Журн. прикл. спектр. 1986. -Т.44. - №2. - С.341-343.
82. Баранов В.И., Грибов JI.A. О возможности построения полуэмпирической адиабатической теории электронно-колебательных спектров сложных молекул // Журн. прикл. спектр. 1988. - Т.48. - №6. - С.963-967.
83. Баранов В.И., Грибов JI.A. К вопросу о построении полуэмпирической неадиабатической теории электронно-колебательных спектров сложных молекул // Опт. и спектроск. 1989. - Т.67. - №1. - С.32-38.
84. Gribov L.A., Baranov V.I. A non-adiabatic semiempirical method for calculating the electron-vibrational spectra of large molecules. A possible approach for future studies //J. Mol. Struct. 1990. - V.224. - №17. - P.45-60.
85. Баранов В.И. Исследование переносимости параметров в ряду родственных соединений в теории электронно-колебательных спектров многоатомных молекул // Журн. прикл. спектр. 1989. - Т.51. - №4. - С.625-628.
86. Эляшберг М.Е., Грибов JI.A., Серов В.В. Молекулярный спектральный анализ и ЭВМ. М.: Наука, 1980. - 307 с.
87. Математические методы и ЭВМ в аналитической химии. Проблемы аналитической химии. T. IX М.: Наука, 1989. - 302 с.
88. Грибов JI.A., Дементьев В.А. Методы и алгоритмы вычислений в теории колебательных спектров молекул. М.: Наука, 1981. - 356 с.
89. Баранов В.И., Савин Ф.А., Грибов JI.A. Программы расчета электронно-колебательных спектров многоатомных молекул. М.: Наука, 1983. - 192 с.
90. Грибов JI.A., Дементьев В.А. Моделирование колебательных спектров сложных соединений на ЭВМ. М.: Наука. Гл. ред. физ.-мат. лит., 1989. - 158 с.
91. Баранов В.И., Грибов JI.A. К теории безызлучательных переходов при оптическом возбуждении газовых сред // Журн. прикл. спектроск. 2000. -Т.67. - №3. - С.289-295.
92. Баранов В.И., Грибов JI.A. Моделирование кинетики внутримолекулярных процессов и нестационарных спектров с учетом безызлучательных переходов. // Журн. прикл. спектроск. 2004. - Т.71. - №4. -С.421-428.
93. Медведев Э.С., Ошеров В.И. Теория безызлучательных переходов в многоатомных молекулах. М.: Наука, 1983. - 280 с.
94. Грибов Jl.А., Маслов И.В. Об одном возможном подходе к моделированию бимолекулярных химических реакций // Журн. физич. химии. 2000. - Т.74. - №3. - С.441-448.
95. Грибов J1.A. О некоторых направлениях изомеризации и образования разветвленных структурных форм ненасыщенных соединений // Журн. физич. химии 2001. - Т.75. - №3. С.547-550.
96. Gribov L.A., Maslov I.V. Methods of estimation of probabilities of photoinduced isomer-isomer structural transformations of polyatomic molecules // J.Mol.Struct. 2000. - V.521. - №1-3. - P.107-120.
97. Грибов JI.А. Высшие приближения в теории изомер-изомерных фотопреобразований молекул // Журн. физич химии. 2001. - Т.75. - №3. С.484-489.
98. Грибов Л.А. К вопросу о поиске путей химических реакций изомеризации и энергий активации // Журн. физ. химии. 2001. - Т.75. - №1. -С.83-88.
99. Грибов Л.А. Колебания молекул и определение путей химических реакций. Изомер-изомерное преобразование молекулы бензола в линейные структуры того же атомного состава // Журн. физич. химии. 2001. - Т.75. -№2 - С.290-293.
100. Грибов Л.А. Возможные направления изомеризации и образование разветвленных структурных форм ненасыщенных соединений. // Журн. физич. химии. 2001. - Т.75. - №3 - С.573-576.
101. Грибов Л.А. Колебания молекул и определение путей реакций присоединения и разложения. // Журн. физич. химии. 2001. Т.75. - №8 -С.1403-1408.
102. Грибов Л.А. Колебания молекул и изомеризация. // Журн. физич. химии. 2001. Т.75. - №9. - С.1638-1641.
103. Грибов Л.А. Колебания большой амплитуды и химическое превращение молекул. // Журн. физич. химии. 2001. - Т.75. - №10 - С.1775-1781.
104. Грибов JI.А. Теоретические модели для расчета кинетики спектров люминесценции с учетом эффекта квантовых биений при изомеризации молекул. // Журн. прикл. спектроск. 2002. - Т.69. - №3. - С.312-317.
105. Грибов Л. А. От теории спектров к теории химических превращений. М.: Эдиториал УРСС 2001 - 368 с.
106. Основы аналитической химии. В 2 кн. Кн.2. Методы химического анализа. Под ред. Золотова Ю.А. 2-е изд., перераб. и доп. - М.: Высш. шк., 1999.-494 с.
107. Грибов Л.А. Безэталонный количественный молекулярный анализ объектов окружающей среды по инфракрасным спектрам // Журн. аналитич. химии 1995. - Т.50. - №6. - С.589-594.
108. Gribov L.A., Elyashberg М.Е., Karasev Yu.Z. // Anal. Chim. Acta. 1995. V. 316. P.217-224.
109. Elyashberg M.E., Gribov L.A., Karasev Yu.Z., Martirosian E.R. IR quantitative analysis of organic mixtures without any calibration // Anal. Chim. Acta. 1997. V. 353. P.105-114.
110. Баранов В.И., Грибов Л.А. О возможности анализа вещества методами спектроскопии в ультрафиолетовой и видимой областях без использования образцов стандартного состава // Журн. аналитич. химии 1999. - Т.54. - №4. - С.350-358.
111. Грибов Л.А., Баранов В.И., Эляшберг М.Е. Безэталонный молекулярный спектральный анализ. Теоретические основы. М.: Эдиториал УРСС-2002-317 с.
112. Грибов Л.А., Баранов В.И., Астахов С.А. О возможности безэталонного анализа вещества методами спектроскопии с временным разрешением в УФ ивидимой областях // Доклады Академии наук. 2000. - Т.374. - №4. - С.493-498.
113. Астахов С.А., Баранов В.И., Грибов J1.A. О безэталонном анализе многокомпонентных смесей веществ методами вибронной спектроскопии с временным разрешением // Журн. аналитич. химии. 2001. - Т.56. - №7. -С.703-713.
114. Арнольд В.И. Обыкновенные дифференциальные уравнения. М.: Наука, 1975. - 240 с.
115. Гантмахер Ф.Р. Теория матриц. М.: Наука, 1967. - 575 с.
116. Желтовский А.В., Данилов В.И. Квантово-механическое изучение электронной структуры оснований нуклеиновых кислот методом CNDO/2. Препринт ИТФ-73-80Р. Киев. 1973. 38 с.
117. Snir М., Otto S., Huss-Lederman S., Walker D., Dongarra J. MPI: The Complete Reference. — Massachusetts.: Massachusetts Institute of Technology, 1996.
118. Gropp B. Tutorial on MPI: The Message Passing Interfase. —* Massachusetts.: Massachusetts Institute of Technology, 1996.
119. Баранов A.B., Лацис A.O., Сажин C.B., Храмцов М.Ю. Руководство пользователя системы МВС-1000/16. М: ИПМ им. М.В. Келдыша РАН.,2001.- 18 с.
120. Warshel A., Karplus М. Calculation of ground and excited state potential surfaces of conjugated molecules. I. Formulation and parameterization // J. Amer. Chem. Soc. 1972. - V.94. - №.12. - P.5612-5626.
121. S. Ami, M. Hliwa, C. Joachim. // Chemical Physics Letters. 2003. - V.367 -P. 662-668.
122. С.S. Lent, В. Isaksen, M. Lieberman. // J. Chem. Soc. 2003 -№ 125 - P. 1056-1063.
123. F.M. Raymo, R.J. Alvarado, S. Giordani, M.A. Cejas. // J. Chem. Soc. 2003 -V.125-P. 2361-2364.
124. A. Saghatelian, N.H. Volcker, K. M. Guckian, V. S-Y. Lin, M.R. Ghadiri. // J. Chem. Soc. 2003 - V.125 - P. 346-347.
125. Грибов JI.А., Сударушкин C.K. Особенности структур и действия логических молекулярных элементов.// Журн. прикладной спектроскопии. -2004. Т.71. - №2. - С. 264-266.
126. Грибов Л.А., Сударушкин С.К. О возможности создания органических реагентов принципиально нового типа.//ЖАХ. 2003. - Т.58. - №5. - С. 480485.
127. Грибов Л.А., Сударушкин С.К. Прогноз возможности создания специфических реагентов, одновременно взаимодействующих с несколькими элементами.//ЖАХ. 2004. - Т.59. - №4. - С. 343-348.
128. Лен Ж.-М. Супрамолекулярная химия: Концепции и перспективы — Пер. с англ. Новосибирск: Наука. Сиб. предприятие РАН, 1998.
129. Материалы диссертации отражены в следующих публикациях:
130. Баранов В. И., Завалий М. В., Грибов Л. А. Моделирование и расчет динамических спектров с учетом изомеризации сложных молекул И Журн. прикл. спектр., 2003. Т. 70. С. 626-634.
131. Баранов В. И., Завалий М. В., Грибов Л. А. Моделирование процессов изомеризации и спектров с временным разрешением сложных молекул//Журн. прикл. спектр., 2004. Т. 71, № 3. С. 295-30.
132. Грибов Л. А., Баранов В. И., Завалий М. В. Метод моделирования динамических вибронных спектров и состояний с учетом межизомерных переходов сложных молекул // Д.А.Н., 2004. Т. 397. № 6. С. 447-450.
133. Завалий М. В. Моделирование динамических вибронных спектров сложных молекул с учетом межизомерных переходов // 3-я Всероссийская конференция «Молекулярное моделирование» апреля, 2003, Москва С. 67.
134. Баранов В. И., Завалий М. В. Моделирование изомеризации сложных молекул на суперкомпьютере типа МВС-1000 // ИВТН, 2004 С. 50.