ХИМИЧЕСКАЯ ФИЗИКА
УДК 5139.194
ПРОБЛЕМЫ КОГЕРЕНТНОГО ОПИСАНИЯ РЕАКЦИЙ ПЕРЕНОСА ЭЛЕКТРОНА *
Д. В. Додин
Исследуются возможности модификаций стохастических уравнений теории Зусмана с целью учета квантовых интерференционных эффектов между базисными состояниями системы. Основываясь на теории Калдейры-Лэггетта (Caldeira-Leggett), рассмотрен ряд возможностей изменения уравнений Зусмана. А именно, учет неравновесного распределения по импульсам на малых временах и вариант приведения уравнения для матрицы плотности на терме к виду Линдблада (Lindblad). Получены модифицированные уравнения Зусмана в пределе сильного трения. В сравнении с оригинальными уравнениями Зусмана показано, что область применения предлагаемых моделей шире.
Часто реакции переноса электрона в конденсированных средах описываются в рамках стохастического подхода [1], [2]. Уравнения Зусмана для двухуровневой системы имеют вид:
Введение
(1)
(2)
(3)
где
диффузионный оператор для /-го состояния;
— продольное время диэлектрической релаксации; кв — константа Больцмана;
Т — температура; д — координата реакции;
Ег — энергия реорганизации;
оператор L — полусумма диффузионных операторов первого и вто-
рого состояний;
U¡ — энергетические поверхности базисных состояний.
* Работа была поддержана грантом РФФИ № 05-03-32680.
Данный подход применялся многочисленными авторами для описания задачи переноса электрона для различных моделей растворителей (см. [3]—[6] обзор [7] и приведенные там ссылки).
Обычно уравнения Зусмана решаются в приближении о локализации недиагональных элементов матрицы плотности в окрестности пересечения термов и исключением из исходной системы уравнения на недиагональные элементы [1]
(Рз1- Рп) = 2я*у(Ри ~ Р22)5^ - Чс), (4)
где — точка пересечения диабетических термов.
Хотя этот прием обладает мощными регуляризующими свойствами, оказывается потерянной всякая информация об эффектах квантовой интерференции состояний. При этом сколько-нибудь подробное описание когерентных квантовых эффектов становится невозможным. Отметим, что в связи с развитием точности и чувствительности современной экспериментальной техники описание этих эффектов становится особенно актуальным.
Кроме того, известно, что в некоторых частных случаях полная система уравнений Зусмана приводит к нефизичным результатам, таким как появление отрицательных населенностей и т. д.
В статье [8] Французов предложил использовать для исследования области применимости полной системы уравнений Зусмана соотношение Фон Неймана для элементов матрицы плотности:
Рп 0м)р22См) ^ |р,2См)|2 ■ (5)
Решая задачу на прямых термах
11{ = Рх, и2 = —Рх, К = О, (6)
где Р — наклон термов;
при начальных условиях
Р„(*, 0) = Р22(*, 0) = р12(х, 0) = 0,58(д - д0), (7)
легко получить
,(0= |р.*Мг =ехр|_д,^
р11(х,/)р„(д:,0 1 3 П2
2 Т2
(8)
Таким образом, уравнения Зусмана не пригодны для описания эволюции системы на [3 П
временах/<хуг= — [8].
2 квТ
2. Уравнения Вигнера — Фоккера — Планка
Заметим, что подход Зусмана основывается на уравнении Смолуховского, описывающем диффузию в потенциале. Известно, что это уравнение может использоваться только на временах существенно превышающих время свободного пробега частицы. Для описания движения частицы на меньших временах необходимо учитывать эволюцию распределения по координатам и импульсам и пользоваться уравнениями типа Фоккера — Планка — Крамерса или более общими [9]. Это соображение кажется разумным и, казалось бы, объясняет проблемы, возникающие при анализе решений на малых временах, которые выглядят как «плата» за потерю информации об импульсном распределении. Таким образом, мы можем принять в качестве рабочей гипотезы необходимость рассмотрения на временах / < \р эволюции импульсного распределения для не противоречивого описания эволюции системы. :
Такое рассмотрение можно провести на основе результатов работ [ 10], [ 11 ]. В случае гармонического или линейного потенциала в представлении Вигнера [12] система уравнений для двухуровневой системы имеет вид:
(9)
(10)
(П)
В работе [11] эта система уравнений получила название ВФП система (Wigner — Fokker — Planck equations).
Уравнения Зусмана (1)—(3) можно получить из этой системы в пределе сильного трения и больших температур (см. напр. [10], [11]).
Проводя для уравнений (9)—(11) анализ, аналогичный описанному в разделе 1, легко получить следующий результат. Отношения матричных элементов r(t), на временах t < хр, где т = 1/у возрастает экспоненциально с показателем 3F2/8ykBT. Таким образом, оказывается нарушенным условие фон Неймана на элементы матрицы плотности, и ее положительную определенность уже нельзя гарантировать.
Из этого мы делаем вывод о невозможности рассматривать систему уравнений (9)—(11) как существенное улучшение уравнений Зусмана.
3. Приведение уравнения Калдейры-Лэггетта к форме Линдблада
Исторически уравнения Зусмана были предложены из эвристических соображений. Однако в последующие годы был проведен ряд исследований [10], [13]—[15], в которых были предприняты попытки их вывода из «первых принципов». Таким образом, уместно обратится к этим работам для преодоления пробелов теории.
Так, в [10] авторы исходили из уравнений, полученных в работе [15]. В последней для системы, связанной с термостатом, было получено уравнение:
f=-j[*,.p]-£fe.[i.pj-£[ï.w)]. <i2)
Важно отметить, что данное уравнение не принадлежит к так называемому классу уравнений Линдблада [16], [17].
^=-^„р]-1^;4р-2^;+ред), (13)
где £. — операторы, описывающие взаимодействие системы со средой;
Я — гамильтониан системы, вообще говоря, перенормированный в результате взаимодействия со средой.
Данное уравнение обобщает уравнение Лиувилля и является наиболее общим видом уравнения, сохраняющим положительную определенность матрицы плотности [18].
Уравнение (12) приводится к желаемой форме добавлением дополнительного слагаемого (см. напр. [19])
0ууп[д,р,1)
dt
dW22(q,p,t) dt
= Ц і 'Щ, (?, Р’ 0 - 2 - ImЩ2 (q, Р, 0, п
L,, W„ (q, р, t) + 2—Im Wn (q, p,t), n
dWn^,P,t) = ^Wv(q,PJ)-2i^(W22(q,p,t)-Wu(q,p,t)), dt ' n
где
d d L. = -p— + — dq dp
d 17, +U2 dq 2
5р
я =-{[н1>р]-^У^,рЬ~ЬАр^р}}-^-{рАр’р1,
01 п п~ п 2
(14)
где 0о = у2/(2/)), В = 2квТ/-(.
В [19] предлагается рассматривать взаимодействие системы с резервуаром как процесс непрерывного измерения. В результате последовательных измерений положение частицы разрешается средой с все большей точностью. В силу принципа неопределенности Гейзенберга этот процесс должен вызывать изменение сопряженного координате импульса. Именно это и обеспечивает дополнительное слагаемое в (14). Таким образом, физический смысл добавленного слагаемого состоит в обеспечении согласования динамики частицы с принципом неопределенности Гейзенберга [19].
Переходя к представлению Вигнера с помощью соответствий Боппа (Ворр) [12], получим модифицированную систему уравнений (9)—(11), где операторы диффузии содержат дополнительное слагаемое:
дЦ ’
где х = 2Б20О.
Применяя к полученным уравнениям вышеизложенный анализ и вводя безразмерное время т = Ц, легко получить:
г(т) = ех{
2Р2
2Х
(15)
1+х +т
(16)
¡2(Х,т) =
(т (2 +т) +X2 (т2 + 2те~т - 2 + 2е х ))2 + X2(т -1 + е~т )2 2т (1 + А,2) - (1 - е~х )(3 -е~%)
(17)
Поведение при малых временах определяется значением управляющего параметра X = 2кТ/уЪ. При X < 1 получаем экспоненциальный рост отношения. Однако при
Х>\ (18)
получаем экспоненциальное затухание на малых временах. Отметим, что данный критерий отличается от полученного в [8] и носит условный характер (соотношение энергетики тепловых флуктуаций и диссипации), в то время как промежуток времени %/г не может быть устранен специальным выбором параметров модели.
4. Выводы
Предлагаемый подход можно использовать для описания когерентных явлений при соблюдении условия (18). Очевидно, что выполнение этого условия предполагает предел больших температур. Таким образом, время декогеренции будет малым и, как следствие, за характерное время декогеренции будет происходить разрушение недиагональных по координате элементов матрицы плотности на терме р^д, ц', 0- Это малосущественно вдали от области пересечения диабатических термов, так как именно в ее окрестности квантовые переходы эффективны и учет когерентности важен.
Очевидно, здесь скрыт существенный недостаток модели, так как фактически термическая декогеренция обеспечивает в этой модели своевременное разрушение недиагональных по диабатическому базису элементов матрицы плотности, в то время как их релаксация должна обеспечиваться независимым механизмом, учитывающим не только взаимодействие со средой, но и особенности взаимного расположения энергетических поверхностей базисных состояний. При таком описании, на наш взгляд, можно было бы ослабить требование высоких температур и изучать когерентные явления при менее жестких условиях.
Автор благодарит профессора А.И. Иванова за ценные замечания и советы, а также С.В. Фесь-кова за помощь в подготовке работы к публикации.
Summary
PROBLEMS OF COHERENT DESCRIPTION OF ELECTRON TRANSFER REACTIONS
D. V. Dodin
Possible modifications of the Zusman stochastic equations aimed to account for quantum interference effect between basic states of the system have been investigated. Starting with Caldeira-Leggett theory we obtain a set of equations in phase space for the Wigner function that includes account for the nonequilibrium distribution in a momentum space at short time scale and a variant of recasting the density matrix equation on a term to the Lindblad form. It has been shown that the applicability of the proposed models is broader compare to the original Zusman equations.
Список литературы
1. Zusman L.D. // Chem. Phys. 1980. № 49. P. 295.
2. Burshtein A.I., Cofman A.G. // Chem. Phys. 1979. № 40. P. 289.
3. Bagchi B., Gayathri N. // Adv. Chem. Phys. 1999. № 107. P. 1.
4. YoshiharaK. //Adv. Chem. Phys. 1999. № 107. P. 371.
5. Rips I., Jortner J. // J. Chem. Phys. 1988. № 86. P. 818.
6. Dodin D.V., Ivanov A.I. // Chem. Phys. 2004. № 303. P. 63.
7. Barzykin A.V., Frantsuzov P.A., Seki K., Tachya M. //Adv. Chem. Phys. 2002. № 123.
P. 511.
8. Frantsuzov P.A. // Chem. Phys. Lett. 1997. № 267. P. 427.
9. Kampen van N.G. Stochastic processes in physics and chemistry. Amsterdam, North Holland, 1990. 6-th ed.
10. Garg A., Onuchic J.N., Ambegaokar V. // J. Chem. Phys. 1985. № 83. P. 4491.
11. Zhang М., Zhang S., Poliak E. //J. Chem. Phys. 2003. № 119. P. 11864.
12. Hillery М., O’Connell R.F., Scully M.O., Wigner E.P. // Phys. Rep. 1984. № 106.
P. 121.
13. Alexandrov I.V. // Chem. Phys. 1980. № 51. P. 449.
14. Ivanov A.I. // Russ. Chem. Phys. 1991. № 10. P. 638.
15. Caldeira A.O., Leggett A.J. // Physica A. 1983. № 121. P. 587.
16. Lindblad G. // Commun. Math. Phys. 1976. № 48. P. 119.
17. Munro W.J., Gardiner C.W. // Phys. Rev. A. 1996. № 53. P. 2633.
18. Sandulescu A., Scutaru H. // Ann. Phys. 1987. № 173. P. 277.
19. Mensky M.B. // Russ. J. UFN. 2003. № 173. P. 1199.