Список литературы
1. Грибов Л. А., Баранов В. И., Зеленцов Д. Ю. Электронно-колебательные спектры многоатомных молекул. М.: Наука, 1997. 476 с.
2. Baker J., Jarzecki A. A., Pulay P. Direct Scaling of Primitive Valence Force Constants: An Alternative Approach to Scaled Quantum Mechanical Force Fields. // J. Phys. Chem. A.1998. V. 102, № 8. P. 1412 -1424.
3. Березин К. В., Нечаев В. В. Расчет частот колебаний и интенсивностей полос в ИК-спектре тетраазапорфи-на методом функционала плотности // Журн. прикл. спектр. 2003. Т. 70, № 2. С. 182-188.
4. Березин К. В., Нечаев В. В., Кривохижина Е. В. При-
УДК 539.182/.184, 519.677
М. К. Березин, М. Л. Чернавина, А. В. Новоселова, В. И. Березин
Саратовский государственный университет E-mail: nusena1975@mail.ru; berezin@sgu.ru
Описаны алгоритмы и программа для расчета электронно-колебательных (вибронных) спектров поглощения и флуоресценции молекул с учетом зависимости интенсивностей от температуры. Алгоритмы основаны на построении квантовых моделей молекул в основном и возбужденном электронных состояниях с применением методов ab initio и теории функционала плотности (DFT-методов). Теоретическое изучение электронно-возбужденных состояний для многоатомных молекул. Предсказание электронных спектров поглощения и флуоресценции многоатомных молекул. Интенсивности полос в спектрах зависят от температуры. Все исходные данные для расчета формируются на основе решения задач молекулярной динамики в рамках квантовых моделей. Предложен параметрический метод учета эффекта Герцберга-Теллера для расчета интенсивностей вибронных полос. Ключевые слова: эффект Герцберга-Теллера, интегралы Франка-Кондона, вибронное взаимодействие, SVL-флуоресцен-ция, матрица Душинского, спектры поглощения.
Temperature Effects in Vibronic Intensity of Absorption and Resonance Fluorescence Spectra
M. K. Berezine, M. L. Chernavina, A. V. Novoselova, V. I. Berezine
Algorithms and program for calculation of the vibronic absorption and fluorescence molecular spectra are described taking into account the dependence of intensities on temperature. The algorithms are based on the construction of molecule's quantum models in ground and excited electronic states with application of ab initio and density functional methods. The theoretic study of excited electronic states for polyatomic
менение метода линейного масштабирования частот в расчетах нормальных колебаний многоатомных молекул // Оптика и спектр. 2003. Т. 94, № 3. С. 398-401.
5. Березин В. И., Богачев Н. В., Абахаева З. М. Применение ЭВМ в молекулярной спектроскопии : учеб. пособие. Саратов : Изд-во Сарат. ун-та, 1989. 92 с.
6. Березин М. К., Тен Г. Н., Березин К. В., Березин В. И. Моделирование спектров 8УЬ-флуоресценции // Проблемы оптической физики и биофотоники : материалы 14-й междунар. молодежной науч. школы по оптике, лазерной физике и биофизике, Саратов, 5-8 октября 2010 г. / под ред. В. В. Тучина, Г. В. Симоненко. Саратов : Изд-во «Новый ветер», 2010. С.212-217.
molecules, prediction of electronic absorption and fluorescence spectra of polyatomic molecules, the temperature dependence of intensities are described. All the initial data for calculation are formed on the basis of molecular dynamics within the framework of quantum model. The parametric technique of accounting Herzberg-Teller effect for calculation vibronic intensities is offered.
Key words: Herzberg-Teller effect, Franck-Condon integrals, vibronic interaction, SVL-fluorescence, Duschinsky matrix, absorption spectra.
Введение
В статье описаны алгоритмы и программа для расчета электронно-колебательных (вибронных) спектров поглощения и флуоресценции с учетом зависимости интенсивностей от температуры. Алгоритмы основаны на построении квантовых моделей молекул в основном и возбужденном электронных состояниях с применением методов ab initio и теории функционала плотности (DFT-методов) [1-3] и современных методов решения задач молекулярной динамики, таких как масштабирование квантовомеханических силовых полей [4, 5], использование шкал для масштабирования частот [6, 7], вычисление эффективных силовых полей [8], оценка ангармонических эффектов [9] и других. Предложен параметрический метод учета эффекта Герцберга-Теллера для расчета интенсивностей вибронных полос.
ВЛИЯНИЕ ТЕМПЕРАТУРЫ НА ИНТЕНСИВНОСТИ ВИБРОННЫХ ПОЛОС В СПЕКТРАХ ПОГЛОЩЕНИЯ И РЕЗОНАНСНОЙ ФЛУОРЕСЦЕНЦИИ
1. Учет температурной зависимости
При расчете вибронных спектров использовались интегралы Франка-Кондона (ИФК) вида
(вдко^сладвд).
Г.,ПГ.(;) Г(/)Г.тГ,.,:И>, . (1)
■о ¡■¡;}Г.(МГЛЛ1!Г:(,! ,
которые с помощью метода перестановочных матриц [10], не выводя дополнительных формул для ИФК, можно привести к виду
{вд^отвдадж«),
и Л/У.-ЛПСЛ^'-М-) ГЦ)ГЛ]).. (2)
■:Г:(/,Г.(./К ..[.,')Г:(п (ОГ.
В такой форме с их помощью можно описывать «горячие» переходы в спектре поглощения. Указанный факт был использован авторами статьи при создании программ для расчета относительных интенсивностей «холодных» и «горячих» переходов и абсолютных интенсивностей линий
в спектрах флуоресценции и поглощения с учетом зависимости интенсивностей от температуры.
В состоянии термодинамического равновесия интенсивность «горячей» полосы поглощения может быть оценена по формуле из [11]
I
(V , ) =
3ке
<
/ (О >
в которой N следует заменить на населенность уровня N , с которого начинается вибронный переход [12]:
Ь п А* ^Д**
,
(-1 п=\
(3)
где п - число осцилляторов (нормальных колебаний) в основном электронном состоянии, участвующих в вибронном переходе; у", со — квантовые числа и частоты колебаний этих осцилляторов; М - число нормальных колебаний молекулы. Подставляя (4) в (3), получаем для абсолютных интенсивностей вибронных полос поглощения на частотах у00 и у в приближении Франка -Кондона:
^ , _ 8тг пужN
3/гс
8
,_, м V
{(о^коо) м#-(о)Па-« *г
<И
,
Ц'
Ъкс
а для относительной интенсивности
кт
7=1
Для компьютерных расчетов удобно выразить абсолютную интенсивность поглощения (5) через относительную (6):
*погл\у} 'Г.
где
:и-),((),| ¡(П ['|с г ■ (7)
(8)
(4)
(5)
(6)
деление по колебательным уровням возбужденной молекулы [12]. Поэтому число вибронных переходов с излучением определяется населенностью этих колебательных состояний N:
м
,
а сам расчет абсолютных интенсивностей выполнять в единицах с1.
Что касается спектров обычной флуоресценции, то большое число фактов и, в частности, независимость контура полосы люминесценции от частоты возбуждающего света [13] показывают, что между актами поглощения и испускания света происходит очень быстрое перераспределение энергии по колебательным степеням свободы. В результате этого при стационарном режиме возбуждения устанавливается температурное распре-
\ А^г ^Г'г/ ] [ :1 , >, (9)
С лГ 1ш\ У
где N — число возбужденных молекул в единице объема; т — число осцилляторов возбужденной молекулы, участвующих в вибронном переходе с излучением; и\, со1 - квантовые числа и частоты колебаний этих осцилляторов.
Подставляя (9) в формулу для интенсивности флуоресценции из [14], получаем для абсолютных интенсивностей вибронных полос в спектре флуоресценции на частотах у00 и у в приближении Франка - Кондона
64тг*п\4 N , ->—2 м
,
64 jrVv4fl Зс3
_,, и п у_
¡r;i ¡г;, ■ .ц,<о),-хр( ;■ У(■■<■>■){ ](i < ■■
ц=\
(10)
Агм
а для относительной интенсивности (по отношению к линии излучения на частоте v00)
-----
см= —
'оо j
{(О') |(0")}
ехр
(11)
Компьютерный расчет абсолютной интенсивности вибронной полосы в спектре флуоресценции выполняется через относительную интенсивность
hi/
r iw кг, '[],i , ,, (12)
I
а сама интенсивность вычисляется в единицах с2:
(13)
Из формул (6) и (11) видно, что относительные интенсивности «холодных» переходов в спектрах поглощения и флуоресценции не зависят от температуры, так как такие переходы начинаются из нулевых колебательных состояний (V = U), и экспоненциальный множитель в этих формулах равен 1.
Соотношения (4)-(5) и (10)—(11) легко обобщаются на случай, когда расчет выполняется с учетом эффекта Герцберга-Теллера (ЭГТ). Для этого в них нужно сделать замену:
(u's)\(VD) м2(0)^MP
(14)
гдеМР - момент перехода, состоящий в общем случае из двух вкладов - франк-кондоновского РК, учет которого описан выше, и герцберг-тел-леровского ГТ.
Для моделирования вибронных спектров поглощения и флуоресценции с учетом зависимости интенсивностей от температуры создан программный комплекс, состоящий из трех программ - РЬШЛББ, РЬШШЛББ, АВБОЯР.
Основные параметры, типы матриц и их структура сохраняются такими же, как и при моделировании спектров БУЬ-флуоресценции [14]. Дополнительно введен параметр ТАВ — абсолютная температура и два управляющих параметра — БУАВ и 1АВБ, которые позволяют по одной и той же программе проводить раздельно расчет интенсивностей в спектрах флуоресценции и поглощения. Программы РЬШ, РЬиШ2 для спектров поглощения рассчитывают интенсивности только «горячих» полос, возникающих
с возбужденных уровней всех осцилляторов, когда их квантовые числа меняются от 1 до 6, т.е. с уровня тона и пяти обертонов всех колебаний (FLU1), а также с уровней их комбинационных парных тонов (FLU1U2).
Для расчета интенсивностей переходов в поглощении с нулевого уровня выделена специальная программа ABSORP.
Программный комплекс позволяет рассчитывать как относительные (по отношению к 0-0 полосе или полосе сравнения), так и абсолютные интенсивности (в единицах с1 (8) или с2 (13)) вибронных полос в спектрах поглощения и флуоресценции многоатомных молекул.
Исходные данные программ
d - смещения нормальных координат, А(а.е.м)^1, (002 — частоты нормальных колебаний в основном электронном состоянии, 103см 001 — частоты нормальных колебаний в возбужденном электронном состоянии, 103см 1, Z - матрица Душинского,
i00 = 1(10, 100...) - шкала для интенсивности 0-0-перехода,
ТАВ - абсолютная температура, °К,
V00 - частота 0-0-перехода, см 1,
I0 - ограничитель минимальной интенсивности,
выводимой на печать.
Управляющие параметры
ND = 1 или 0 - расчет с учетом или без учета эффекта Душинского,
NU = 1 или 0 - вводится или не вводится матрица частотных интервалов,
SVAB = 1 или 0 - расчет «горячих» полос в спектрах флуоресценции или поглощения, IABS = 1 или 0 - расчет абсолютных или относительных интенсивностей в спектрах флуоресценции и поглощения.
Все необходимые исходные данные формируются на основе решения задач молекулярной динамики в рамках квантовых моделей по базовой программе [3]. Избранными являются квантовые модели, построенные с применением методов ab initio и теории функционала плотности (DFT-методов).
Рассмотрим блок-схему программного комплекса для компьютерного моделирования ви-бронных спектров поглощения и флуоресценции многоатомных молекул с учетом зависимости интенсивностей от температуры.
Программы
РЬШЛББ РЬШШЛББ АВБОИР
Задающие матрицы для отдельных типов
®
переходов
NS11 NS22 NS33
1 1 1 1 1 1
2 1 2 1 2 1
3 1 3 0 3 1
4 0 4 0 4 1
5 0 5 0 5 1
6 0 6 0 6 1
7 0 7 1 7 1
8 1 8 0 8 0
9 0 9 0 9 0
10 0 10 1 10 0
11 0 11 0 11 1
12 0 12 0 12 1
13 0 13 1 13 0
Матрицы для задания линий сравнения
Матрицы частотных интервалов для выделения отдельных линий
|nui[
NU2
NU3
Матрицы частот и вычисленных интенсивностей выделенных линий
Матрицы параметров и вычисленных интенсивностей вибронных переходов
Построение вибронного спектра с помощью программы моделирования и визуализации контуров спектральных линий 1ХСОКТОиК.
Задающие матрицы определяют 13 типов вибронных переходов (первый столбец), каждый из которых или несколько из них одновременно могут быть актуализированы подстановкой 1 в правый столбец.
2. Параметрический учет эффекта
Герцберга-Теллера
Программный комплекс можно использовать и для моделирования вибронных спектров в учетом ЭГТ. С этой целью (14) можно привести к виду [15]
МР(у) = ((0,)|(0"))Щ(0)(Б, +|Б,) (15)
и ввести для относительного момента ЭГТ параметр
5:
MJ0)
с. ■
(16)
в который входит St— производная от момента перехода по нормальной координате.
Рассмотрим случай, когда за счет ЭГТ в ви-бронном спектре проявляются антисимметричные колебания представления Г.
Если электронный переход разрешен на ком -поненте Т, а «заимствование» интенсивности происходит по компоненте с, то относительные интенсивности полносимметричного и антисимметричного фундаментальных колебаний будут равны
OJ = — ^
00
С:"'.,)- ! К' !,
ПиЛ > гГ V(;'-(i 2.1!..,
где
G, =
(17)
(18)
(19)
(20)
0)
а суммирование проводится по колебаниям представления Г, число которых равно к.
Из (17) и (18) видно, что относительную интенсивность антисимметричного колебания можно вычислять так же, как и полносимметричного колебания, если ввести матрицу Б"' с компонентами
V
Bas =
В-
(21)
Приведенные смещения нормальных координат А для полносимметричных колебаний, когда матрица В считается известной, находятся по соотношению
Д' = Л(\ + рах)~1 В = &{\-2А)~1 В. (22)
Введем матрицу эффективных приведенных смещений для антисимметричных колебаний представления Г, заменив в (22) В на Б" :
ДГ (23)
Введение матрицы А1, позволяет при компьютерном моделировании учесть вклад антисимметричных колебаний в вибронный спектр поглощения. Для этого следует выполнить два расчета.
Первый расчет выполняется в приближении Франка - Кондона, когда приведенные смещения для всех антисимметричных колебаний равны 0. В результате будут вычислены интенсивности всех полносимметричных колебаний, включая четные обертоны антисимметричных колебаний.
Второй расчет выполняется, когда для антисимметричных колебаний представления Г вводятся приведенные смещения Д1 . Результаты такого расчета будут содержать не только интенсивности переходов для антисимметричных колебаний, активных по ЭГТ, но и интенсивность комбинационных переходов с участием полносимметричных колебаний и колебаний представления Г.
Выделяя эти колебания из расчетных данных и суммируя с результатами первого расчета, получаем полный вибронный спектр поглощения, который теперь содержит наиболее интенсивные вибронные полосы.
Если электронный переход запрещен М^ (0) = = 0, то в (20) следует положить М^ (0) = 1 и выполнить расчет по тому же алгоритму, что и для разрешенного перехода, но с одним исключением, состоящим в том, что относительная интенсивность вибронного перехода теперь должна вычисляться по отношению к линии сравнения, которая может быть задана с помощью матрицы
1БТЗ.
Для спектров флуоресценции расчет полного вибронного спектра может быть выполнен аналогичным образом. Для этого определим относительную интенсивность антисимметричного фундаментального колебания по соотношению
Г V V
/"""(к )= —
погя \ ш
где
Введем матрицу Das с компонентами
КГ
г>; -Vt,; /-..
(24)
(25)
(26)
и с ее помощью вычислим матрицу эффективных приведенных смещений:
V' V V/J- ■ /) ■. (27)
Для получения полного вибронного спектра излучения расчет выполняется, как и в случае спектра поглощения, в два шага - в приближении Франка - Кондона и с учетом эффективных смещений Д2 для антисимметричных колебаний представления Г.
Если активными по ЭГТ оказываются полносимметричные колебания, то в этом случае нельзя проводить два альтернативных расчета, как в случае антисимметричных колебаний, так как в (15) оба члена Bs и Bs при расчете интенсивностей могут быть одновременно отличными от нуля. При наличии экспериментального вибронного спектра можно ограничиться расчетом интенсивностей в приближении Франка - Кондона, если путем вариации смещений нормальных координат удается воспроизвести опытное распределение относительных интенсивностей в спектре. Это касается прежде всего прогрессий. Если такого согласия не наблюдается, то следует учесть вклад ЭГТ. Интенсивность полносимметричных колебаний в поглощении в этом случае может быть представлена в частном случае, когда электронный переход разрешен по компоненте Т , а переход в третье электронное состояние z разрешен по компоненте с, в виде
(28)
где
а Of определено в (20).
Для относительных интенсивностей в спектре флуоресценции аналогичным образом получаем:
где
С" _ 00 -
1 +
» V Уь-Л .
j=I J
(31)
Чтобы исключить влияние параметров S00 и S00 в (28) и (30), удобнее при моделировании вычислять относительные интенсивности не по отношению к 0-0 полосе, а по отношению к линии сравнения.
Если электронный переход запрещен (Mif (0) = 0), то в (28) и (30) следует приравнять
0 относительный ИФК (^(0') | (0")} | =0), в G'sa, G"sa (20) положить Mf (0) = 1, в (29) и (31) опустить
1 в скобках. Такая опер ация исключает франк-кондоновский вклад FK в интенсивность.
Выводы
Представлена блок-схема программного комплекса для компьютерного моделирования вибронных спектров поглощения и флуоресценции многоатомных молекул с учетом зависимости интенсивностей от температуры. Все необходимые исходные данные для комплекса формируются на основе решения задач молекулярной динамики в рамках квантовых моделей. Предложен параметрический метод учета эффекта Герцберга - Теллера для расчета интенсивностей вибронных полос.
Список литературы
1. Кон В. Электронная структура вещества - волновые функции и функционалы плотности // Успехи физических наук (Нобелевские лекции по химии - 1998). 2002. Т. 172, № 3 (март). С. 336-348.
2. Попл Дж. Квантово-химические модели PDF // Успехи физических наук. 2002. Т. 172, № 3 (март). С. 349-356.
3. Frisch M. J., Trucks G. W., Schlegel H. B, Colthup N. B, Daly L. H., Wiberley S. E. Introduction to Infrared and Raman spectroscopy. Gaussian 03, Revision B. -Gaussian Inc., Pittsburgh PA, 2003.
4. Baker J., Jarzecki A.A., Pulay P. Direct Scaling of Primitive Valence Force Constants: An Alternative Approach to Scaled Quantum Mechanical Force Fields // J. Phys. Chem. A. 1998. Vol. 102, № 8. P. 1412-1424.
(30)
5. Березин К. В. Нечаев В. В., Кривохижина Е. В. Применение метода линейного масштабирования частот в расчетах нормальных колебаний многоатомных молекул // Оптика и спектр. 2003. Т. 94, № 3. С. 398-401.
6. Yoshida H., Takeda K., Okamura J., Ehara A., Matsuu-ra H. New Approach to Vibrational Analysis of Large Molecules by Density Functional Theory : Wavenum-ber-Linear Scaling Method // Chem. Phys. Lett. 2000. Vol. 325, № 4. Р. 477-483.
7. Березин К. В., Нечаев В. В., Кривохижина Е. В. Применение метода линейного масштабирования частот в расчетах нормальных колебаний многоатомных молекул // Оптика и спектр. 2003. Т. 94, № 3. С. 398- 401.
8. Березин К. В., Черняев С. Н., Кирносов Н. А., Березин В. И. Учет базисной суперпозиционной ошибки при расчетах межмолекулярных расстояний и термодинамических характеристик водородосвя-занных комплексов // Проблемы оптической физики и биофотоники (Optical physics and biophotonics) : материалы 12-й междунар. молодежной научной школы по оптике, лазерной физике и биофизике. 23-26 сентября 2008 г., Саратов. SFM-2008 / под ред. В. В. Тучина, Г. В. Симоненко. Саратов : Изд-во «Новый ветер», 2008. С. 176-180.
9. Березин К. В., Нечаев В. В., Элькин П. М. Ангармонический анализ колебательных состояний пиримидина методом функционала плотности // Оптика и спектр. 2004. Т. 97, № 2. С. 221-231.
10. Березин В. И., Богачев Н. В., Абахаева З. М. Применение ЭВМ в молекулярной спектроскопии : учеб. пособие. Саратов : Изд-во Сарат. ун-та, 1989. 92 с.
11. Березин М. К., Тен Г. Н., Березин К. В., Березин В. И. Моделирование спектров SVL-флуо-ресценции // Проблемы оптической физики и биофотоники (Optical physics and biophotonics) : материалы 14-й междунар. молодежной научной школы по оптике, лазерной физике и биофизике, 5-8 октября 2010 г., Саратов. SFM-2010 /под ред. В. В. Тучина, Г. В. Симоненко. Саратов : Изд-во «Новый ветер», 2010. С. 212-217.
12. Ельяшевич М. А. Атомная и молекулярная спектроскопия. М. : ГИФМЛ, 1962. 891 с.
13. Степанов Б. И., Грибковский В. П. Введение в теорию люминесценции. Минск : Изд-во АН БССР, 1963. 352 с.
14. Березин М. К., Тен Г. Н., Березин К. В., Бере-зин В. И. Решение обратной электронно-колебательной задачи по спектрам SVL-флуоресцен-ции // Проблемы оптической физики и биофото-ники (Optical physics and biophotonics) : материалы 14-й междунар. молодежной научной школы по оптике, лазерной физике и биофизике, 5-8 октября
2010 г., Саратов. БЕМ-2010 /под ред. В. В. Тучина, Г. В. Симоненко. Саратов : Изд-во «Новый ветер», 2010. С. 218-224.
15. Грибов Л. А., Баранов В. И., Зеленцов Д. Ю. Электронно-колебательные спектры многоатомных молекул. М. : Наука, 1997. 476 с.