第3O卷第5期
2016年1O月
黑龙江工程学院学报
Journal of Heilongjiang Institute of Technology
Vo1.3O,NO.5
0ct.,2016
考虑Duschinsky转动效应的电荷转移速率理论
陈九菊
(黑龙-;x- ̄程学院电气与信息工程学院,黑龙江哈尔滨150050)
摘要:根据费米黄金规则,在位移简谐振子近似下,采用路径积分方法推导了一个高斯类型的关联函数形式,并最
终发展了一个考虑Duschinsky转动效应的电荷转移速率理论。文中方法的优点在于不需要计算复杂的弗兰克一康
登因子,从而将Dusehinsky转动效应自然地包括到理论形式之中。在位移简谐振子近似下,所发展的电荷转移速率
形式可以回归到已有的电荷转移速率理论。
关键词:电荷转移;路径积分;Duschinsky转动;费米黄金规则;关联函数
中图分类号:0649 文献标识码:A 文章编号:1671—4679(2016)05—0034—04
Charge transfer rate formalism with Duschinsky rotation effect
CHEN Jiuju
(College of Electrical and Information Engineering,Heilongjiang Institute of Technology,Harbin 150050,China)
Abstract:Based on the Fermi—golden rule,the charge transfer rate in the displaced oscillator case has ever
been derived,but the Franck-Condon factor has been explicitly evaluated.In the present paper,a new
formalism for charge transfer rate with Duschinsky rotation effect is developed through the path integral of
Gaussian type correlation function.The advantage of the present path—integral formalism is that it is no
Jonger to compute the Franck—Condon factor and the Duschinsky rotation effect can be easily included.The
new charge transfer rate formalism can easily go back to the previous formalism under the displaced
oscillator approximation.
Key words:charge transfer;path integral;Duschinsky rotation;Fermi—golden rule;correlation function
理论描述电荷转移过程一直是凝聚态领域一
个研究方向[1]。从给体到受体的标准电荷转移图像
需要考虑这两个态之间的电子耦合,并受到核运动
的影响。通常情况下,核自由度被描述为简谐振
子,从而形成热库。基于这种图像,一些数值方法
被发展出来,其中包括直接路径积分方法l_2。]、质心
近似[4]以及其他非微扰方法[5瑁]。这些方法通常被
用于处理一维反应事件,这种事件会与简谐振子构
成的热库发生相互作用。然而,实际化学反应过程
中的反应隧穿需要包括几个自由度,因此需要全量
子力学来处理。在这种情况下,费米黄金规则对于
速率理论是一个有效的工具Lg ,其常常被用于分
析实验上的电荷转移速率[1 。
收稿日期:2016—02—29
基金项目:黑龙江工程学院青年科学基金资助项目(2014o3 10)
作者简介:陈九菊(1981一),女,讲师,硕士,研究方向:半导体器件
费米黄金规则方法在计算体系的电荷转移速
率时具有一定的局域性。首先,因为费米黄金规则
是基于一级含时微扰理论,所以其要求初态与末态
之间的电子耦合必须足够小。其次,初始势能面的
布局在非绝热弛豫过程中要始终保持在热平衡。
在这两个条件下,费米黄金规则方法将电荷转移速
率描述为对某一个函数的时间积分,但很难精确求
解[1 。随后,Lin等人进一步发展了该电荷转移速
率的表达形式,使之可以精确用于计算电荷转移速
率常数 引。
电荷转移是化学反应过程中的一个重要的动
力学过程,因此,发展电荷转移速率理论对于理解
电荷转移过程至关重要,也是理论化学一个研究方
向。理论工作者发展了一些适用于任意电子耦合
强度的电荷转移理论,这些理论可以将非绝热极限
和绝热极限下的电荷转移速率联系起来。然而,在
电荷转移过程中,势能面会在电荷转移时发生转
第5期 陈九菊:考虑Duschinsky转动效应的电荷转移速率理论
动,这被称之Duschinsky转动效应。
Duschinsky转动效应在无辐射过程中至关重
要l_1 ,而电荷转移过程就是其中一个无辐射过程。
如何在电荷转移速率理论中考虑Duschinsky转动
效应是一个难题。本论文的目的就是根据路径积
分方法将基于费米黄金规则的电荷转移速率进行
重新推导,并把Duschinsky转动效应包含到该电荷
转移速率形式中。
1 考虑Duschinsky转动效应的电荷
转移速率
根据费米黄金规则,一般的量子力学电荷转移
速率具有以下形式[¨ :
尼口==:等l yl ∑∑P I( ,l >『。 ( , ). ,‘
u tJ
(1)
其中: 为初态和末态之间的电子耦合,
『<@而,I >l。为弗兰克一康登因子。初态和末态的
振动波函数可以分别写为
『 一 })( ,l )一 I >・(2)
其中,I >和『 知 >一维简谐振子哈密顿量的本
证态,它们相应的本证能量分别为
E ===( +吉)充 , , 一(u +丢) cu .
(3)
同时,初态的分布函数可以表达为
P 一[ 唧( kBT 1]].J-1expI- E.)一
[ 唧f生kBT/ -] ̄xpi k-7E ̄ ]一
I ]7 2sinh exp( I I.c4
当把初态的分布函数重新定义为P 一
ⅡP ,并把狄拉克函数表达为对时间的傅里叶积
形式,方程(])就被萤新表达为
志盯一 1 舻 I IP f …p it +丢) , 一( + 1) ))=:=
1
一 ̄V
j
[Gi㈩・ (5)
其中Gj具有形式
一 xp( + 一( + 1) ㈤
壳∞ 是初态和末态之间的能量差。在Lin等人的工
作中[坞],方程(5)是在位移简谐振子近似下,通过
明确计算弗兰克一康登因子得到的,即△Q—Q 一
QJ和叫 一∞ 。因此,Duschinsky转动效应很难被
包含进他们的理论形式中。本文将采用路径积分
形式的高斯类型关联函数来推导公式(5),从而可
以很自然地把Duschinsky转动效应包含进来。
将方程(4)代入到方程(6), (z)可以被写为
一 唧(斋) -1 I- E.) 唧[砉c 一 ]一
[ 唧(斋) 】expE— (~ + t,E 1 expI— (一砉) 一
[ exp k=T)l~ J exp[~ (一 +舌)自 ]exp[一 (一砉)方 ]l ,一
1I G 一 (2sinh ) 旧,.
这里,关联函数X(r ,rf)具有以下形式:
(7)
(8)
(9)
黑龙江工程学院学报 第30卷
其中,r 一一 +砉, 一一舌。
当考虑Duschinsky转动效应时,初态的正则坐
标Q与末态的正则坐标Q 可以通过Duschinsky转
动矩阵S和振动模位7f'CD联系起来,即
Q 一SQ+D. (10)
这里,具有下划线的变量代表列矢量,而黑体变量
代表方矩阵。下面,初态的正则坐标用兰,
一
Y来表示,
末态的正则坐标用z ,yt来表示。当把关联函数对
于初态坐标求迹时,可以得到
x(r , )一I d <兰I e-iriH ̄e-irf'Ft 1 >一
I d—zd—yd—x d—y <兰1 ly一)(一y l—yl>< 1 f ̄tfu Ix一 >< I >. (11)
当a(r)和b(r)表示对角矩阵,并且其对角矩阵元
分别为
bs(r)一面c壳oj , (13)
就可以得到以下两个数学形式
<
_
x l e-胡u l >一 f
(
d
2
e
7c
ta
方
(r
)
)
・
exp{吾[丢( 6(r)三+ (r) )~兰 n(r) ]).
E===b(r r)一a(r,). (18)
并把方程(14)和(15)带人方程(11)后,关联函数的
形式就可以表达为
(Ti r,)::= /d
(2
et a(
疗
r
)
l )√ j-电 ・
eXp{砉l专( Bx一-F—y By一)~ A +
D ( +y)+DTED]). (19)
这里,每一个振动模式的初态和末态频率被认为是
一致的。
当把矩阵K以及列矢F和z分别定义为
(14) K::=l _A l, (20)
< I z>一 r ,一fS z+D、] f1 5 L—A B _J
当定义矩 。一 一 一~
fT一[里ES,DES3, (21
A—n( )q-STa(rr)s, (16) 三 E—
xT
,
一
y -7・ (22)
B一6(ri) ̄-S b(r/ )S, (17) 关联函数就变为
( , )一 ̄/d
、eta(z'i) √ 蓑 .『d兰exp{ i L 1兰 + 兰+ E ])一
deta(r
d ̄)
e
d eta(rf)
・exp{吾[一号 一F+一D ]}. (23)
在把方程(8)N1(23)代人到方程(5)后,最终的考虑Duschinsky转动效应的电荷转移速率形式为
愚 一 II(2si )×
J。d  ̄/deta(rdi)ed eta(ri)exp{拓∞ +吾[一号 一F+~D ED]). (24)
2 没有考虑Duschinsky转动效应的 关联函数就变为
电荷转移速率形式 ( , )===√ § ・
:
略DuS出 sk ’ exp{吾[里TE(B—A)1G旦]). (27)D 子近似时
,uschinsky转动矩阵S变为单位矩阵, ¨。一 一
初态正则坐标Q和末态正则坐标Q,满足 为了获得Lin等人的电荷转移速率形式,这里
铁,方稗
Q' = Qff -D.
(16 钳阵当
需要弓l入新的变量SJ一 一麦挪么, 显然,方程 )~(18)变为对角矩阵。当定义一个 一
新的对角矩阵,即 exp{ ̄-ED E(B—A)一 G ])一
南
第5期 陈九菊:考虑Duschinsky转动效应的电荷转移速率理论 ・ 37・
Ⅱexp(吾 EJ(B 一Aj) }=
exp{一∑SJ[(2 +1)-n一,e一 一( +1) ]),
(28)
/deta(Ti)deta( )
detBdet(B—All--1A)
8一 一( +1)P ]). (31)
3 结 论
本文采用路径积分技术,推导了一个高斯类型
的关联函数,从而最终发展了一个考虑Duschinsky
转动效应的电荷转移速率理论。本文发展的方法
不再需要计算复杂的弗兰克一康登因子在位移简
谐振子近似下,可以回归到文献上已有的没有包含
Duschinsky转动效应的电荷转移速率理论。因此,
该方法更具有普适性。该方法为今后数值研究
Duschinsky转动效应对电荷转移的影响提供了理
论基础。
[1]
E2]
[3]
参考文献
ULSTRUP J.Charge Transfer Processes in Condensed
Media[M].Berlin:Springer-Verlag,1 979.
MAKRI N.Quantum dissipative dynamics:A numeri—
cally exact methodology FJ].J Phys Chem A,1998,
102:4414—4427.
G0L0S0V A A,FRIESNER R A,Peehukas.Reduced
dynamics in spin-boson models:A method for both
slow and fast bath[J].J Chem Phys,2000,112:
2095一Zl05.
[43 GEVA E,sHI Q,VOTH G A.Quantum-mechanical
reaction rate constants from centroid molecular dynam—
ics simulations[J].J Chem Phys,2001,115:
9209—9222.
[5]HANGGI P,TALKNER P,BORKOVEC M.Reaction-
rate theory:fifty years after Kramers EJ].Rev Mod
Phys,1990,62:251—341.
[6]P0MYAI OV A,TANNOR D J.The non-Markovian
quantum master equation in the collective-mode repre—
sentation:Application to barrier crossing in the inter—
mediate friction regime[J].J Chem Phys,2005,
123:2041l1.
[7]MAKR1 N.Time-dependent quantum methods for
large systems[J].Annu Rev Phys Chem,1999,50:167—
191.
[8] MILLER W H.Spiers memorial lecture quantum and
semiclassical theory of chemical reaction rates[J].Far—
aday Discuss,1998,l10:1—21.
[9]FAIN 13.Theory of Rate Processes in Condensed Media
[M].Berlin:Springer-Verlag,1980.
Elo]WEISS U..Quantum Dissipative Systems.2nd edEM].
Singapore:W0rld Scientific,1 999.
[11]BIXON M,J()RTNER J.Electron transfer-from isola—
ted molecules tO biomolecules[J].Adv Chem Phys,
1999,106:35—202.
r12]NERIA董:,NITZAN A.Semiclassical evaluation of
nonadiabatic rates in condensed phases[J].J Chem
Phys,1993,99:1109—1123.
[13]LIN S H,CHANG C H,LIANG K K,et a1.Ultrafast
dynamics and spectroscopy of bacterial photosynthetic
reaction centers[J].Adv Chem Phys,2001,121:1-88.
r14]LIN S H.Rate of interconversion of electronic and vi—
brational energy[J].J Chem Phys,1966,44:
3759—3767.
[15]JORTNER J.Temperature dependent activation energy
for electron transfer between biological molecules[J].J
Chem Phys,1976,64:4860—4867.
[责任编辑:刘文霞]
本文发布于:2023-01-03 15:03:54,感谢您对本站的认可!
本文链接:http://www.wtabcd.cn/fanwen/fan/90/85167.html
版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。
留言与评论(共有 0 条评论) |