计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: aristotleemma
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合交流] 轨道耦合常数计算

[复制链接 Copy URL]

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

16#
发表于 Post on 2016-11-12 03:55:28 | 只看该作者 Only view this author
janeli1995 发表于 2016-11-11 10:24
您好,谢谢!
Troisi, A., Orlandi, G., Troisi, A., & Orlandi, G. (2006). Dynamics of the Intermole ...

我觉得这个文章里的波动很大啊?基本上是在0 到2000 之间波动。难道你的波动还能更大?

107

帖子

0

威望

281

eV
积分
388

Level 3 能力者

17#
发表于 Post on 2016-11-12 05:47:53 | 只看该作者 Only view this author
wbn 发表于 2016-11-11 05:04
OK 我看了程序提到的JACS文章, 里面有说 the transfer integral (Jijeff ≡ t) can be estimated as one- ...

这个看2-level system,对角元的coupling J的数值(小可以看做微扰),在最后对角化后,导致的能量差2J

所以transfer integral的数值会取delta E的1/2
Marius 理论在能量转移方向的对应物就是 Foster theory,本质的数学结构是一样的。与此相关的还有物理里的spin-boson模型,在化学NMR和超快光谱研究分子体系动力学里模型本质都是一样的

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

18#
发表于 Post on 2016-11-14 06:11:36 | 只看该作者 Only view this author
wswrpd 发表于 2016-11-12 05:47
这个看2-level system,对角元的coupling J的数值(小可以看做微扰),在最后对角化后,导致的能量差2J
...

多谢。但我没明白哎。这个deltaE 是指哪两个态之间的能量差?看文章是HOMO 到HOMO-1,这两个轨道之间的能级跟这个电子转移过程有什么关系啊?

59

帖子

0

威望

206

eV
积分
265

Level 3 能力者

19#
发表于 Post on 2016-11-14 06:39:16 | 只看该作者 Only view this author
wbn 发表于 2016-11-12 03:55
我觉得这个文章里的波动很大啊?基本上是在0 到2000 之间波动。难道你的波动还能更大?

嗯……我的variance是他的五倍甚至更多……

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

20#
发表于 Post on 2016-11-15 04:07:52 | 只看该作者 Only view this author
wbn 发表于 2016-11-14 06:11
多谢。但我没明白哎。这个deltaE 是指哪两个态之间的能量差?看文章是HOMO 到HOMO-1,这两个轨道之间的能 ...

以中性分子的电子结构为准,如果是研究空穴转移,就是HOMO和HOMO-1。如果是电子转移,就是LUMO和LUMO+1。

107

帖子

0

威望

281

eV
积分
388

Level 3 能力者

21#
发表于 Post on 2016-11-15 04:39:40 | 只看该作者 Only view this author
wbn 发表于 2016-11-14 06:11
多谢。但我没明白哎。这个deltaE 是指哪两个态之间的能量差?看文章是HOMO 到HOMO-1,这两个轨道之间的能 ...

deltaE 可以认为是从电子(或者能量)转移的初始态的能量到结束态的能量差,因为你计算得到的结果都是Adiabatic states 的representation,而所谓的coupling J是在diabatic representation下的耦合系数。

文章是HOMO 到HOMO-1, HOMO应该是没有填满的吧,这样HOMO得到电子而HOMO-1失去电子形成空穴,如果得到的这样的一个激发态的电子分布与基态的电子分布有明显的空间上的区别,这就是一种电子transfer的现象了,而deltaE 就是激发态的能量和基态能量的差值(adiabatic representation)。当然由于没看这篇文章,所以上面只是自己根据以往的经验写的

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

22#
发表于 Post on 2016-11-16 00:01:11 | 只看该作者 Only view this author
wswrpd 发表于 2016-11-15 04:39
deltaE 可以认为是从电子(或者能量)转移的初始态的能量到结束态的能量差,因为你计算得到的结果都是Adi ...

嗯。如果是A-B-A',其中A和A'是对称的电子给/受体,B是桥(当然也可以是真空,这就变成Dimer了),那么HOMO和HOMO-1表示空穴的“转移积分”的绝对值的两倍。

所谓空穴就是中性分子敲掉一个电子得到的。如果假设Koopmans'定理成立,并不需要去真的计算那个正离子态,计算一个中性分子(或者dimer)的电子结构可以一次得到空穴的转移积分(HOMO/HOMO-1)和电子的转移积分(LUMO/LUMO+1)。

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

23#
发表于 Post on 2016-11-16 00:05:06 | 只看该作者 Only view this author
wswrpd 发表于 2016-11-15 04:39
deltaE 可以认为是从电子(或者能量)转移的初始态的能量到结束态的能量差,因为你计算得到的结果都是Adi ...

我觉得并不能说“deltaE 可以认为是从电子(或者能量)转移的初始态的能量到结束态的能量差”,因为Marcus型的电子传递在真正的“电子传递/隧穿”这一步是等能量的,没有原子核构象的变化。这一步之前之后的原子核构象变化导致的能量变化属于“重组能”的范畴。

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

24#
发表于 Post on 2016-11-16 00:08:49 | 只看该作者 Only view this author
wbn 发表于 2016-11-11 05:04
OK 我看了程序提到的JACS文章, 里面有说 the transfer integral (Jijeff ≡ t) can be estimated as one- ...

用GMH要先算出一个电荷转移激发态吧?以前试过但感觉不是很容易找到……

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

25#
发表于 Post on 2016-11-16 00:57:37 | 只看该作者 Only view this author
stecue 发表于 2016-11-15 04:07
以中性分子的电子结构为准,如果是研究空穴转移,就是HOMO和HOMO-1。如果是电子转移,就是LUMO和LUMO+1。

OK I see,谢谢。因为我以前只接触过闭壳层体系,所以只知道HOMO 到 LUMO 的转移。

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

26#
发表于 Post on 2016-11-16 01:17:31 | 只看该作者 Only view this author
wbn 发表于 2016-11-16 00:57
OK I see,谢谢。因为我以前只接触过闭壳层体系,所以只知道HOMO 到 LUMO 的转移。

呃……HOMO到LUMO一般应该是 charge transfer excitation 吧?

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

27#
发表于 Post on 2016-11-16 01:27:49 | 只看该作者 Only view this author
wswrpd 发表于 2016-11-15 04:39
deltaE 可以认为是从电子(或者能量)转移的初始态的能量到结束态的能量差,因为你计算得到的结果都是Adi ...

多谢指点。那么这个delta E 大致就相当于Marcus 理论中的活化自由能,就是那个二项式吗?其实我对这些概念一直有一点混乱不清, 正好二位在这里我们一起讨论下。首先这个adiabatic 和 non-adiabatic过程是否就如这篇文章中描述的http://www.annualreviews.org/doi ... schem.012809.103324 adiabatic 过程有一个隧穿的过程,而non-adiabatic 是直接从交点过去。而Marcus 方程(上面文章方程13)就是处理non-adiabatic 过程,也就是electron coupling 比较弱的过程,这么理解对不对?

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

28#
发表于 Post on 2016-11-16 01:33:44 | 只看该作者 Only view this author
stecue 发表于 2016-11-16 01:17
呃……HOMO到LUMO一般应该是 charge transfer excitation 吧?

嗯。因为我做的是photo-induced electron transfer, 首先激发到一个 local excited state, 然后在到CT 态。我的体系CT态相当容易找,只要把激发态各片段电荷加一加就找到了,或者用电荷转移距离delta R 判断,或者用激发态的dipole moment判断,都非常清晰。(而且我所有的case第一激发态就是CT态)。但是不激发的能量转移过程我就不会了。您能帮我看一下我上一楼关于non-adiabatic 过程的理解对不对吗?

282

帖子

0

威望

3020

eV
积分
3302

Level 5 (御坂)

29#
发表于 Post on 2016-11-16 01:35:06 | 只看该作者 Only view this author
wbn 发表于 2016-11-16 01:27
多谢指点。那么这个delta E 大致就相当于Marcus 理论中的活化自由能,就是那个二项式吗?其实我对这些概 ...

那为什么marcus 方程却有electron coupling 这一项呢?

230

帖子

0

威望

919

eV
积分
1149

Level 4 (黑子)

30#
发表于 Post on 2016-11-16 04:52:50 | 只看该作者 Only view this author
本帖最后由 stecue 于 2016-11-22 03:55 编辑
wbn 发表于 2016-11-16 01:33
嗯。因为我做的是photo-induced electron transfer, 首先激发到一个 local excited state, 然后在到CT 态 ...

你的体系很理想哈,不错不错。我们实验上用的体系就比较复杂,理论算一下往往越算越confusing了,呼呼。

Adiabatic 和 non-adiabatic 有多种相关但区别很微妙的含义。从根本上说,量子力学里 adibatic 的意思是缓变的,所谓“缓”的意思是,外场变化足够慢,以至于体系每一时刻的“不含时薛定谔”的本征态,都可以良好定义。这也经常翻译成“绝热的”。对于量子化学体系,所谓的“缓变”,就是原子核(作为外场),它们的运动速度比电子运动速度慢很多,所以可以在任意时刻应用“波恩-奥本海默”近似。这意味着什么呢?就是在热致电子传递的整个过程中“势能面”总是良好定义的,而且永不交叉。——一般来说,电子传递一开始总是在一个非简并的基态吧?那么只要符合“原子核总是运动的足够慢”的情况,我们就有“量子绝热定理”保证,这个体系将一直处于基态,而不会跑到激发态去!

理解了 adibatic,也就理解了non-adiabatic。简单的说,既然adibatic是无交叉点,那non-adibatic显然就是有交叉点的了。当然你细想一下这其实有些问题:所谓的“交叉点”,本来就是定义在“势能面”上的。而按照上面的说法, adiabtic 情况才有良好定义的势能面,那non-adiabatic情况下,BO近似就失效了,电子和原子核强烈耦合在一起(vibronic coupling很大),根本就没法定义“势能面”,又如何定义“交叉点”呢?——没错你可以强行把原子核固定住计算电子态,但是按照Marcus的理论,你将得到两个完全简并的电子态:电子传递前和电子传递后。但是,如果真的是完全简并,那这两个态之间的相互作用就必须完全是零。如果用微扰论算一下,这就完全不可能发生电子传递。实际上,由于 adiabatic states 是能量本征态,也就说是哈密度算符(SCF平均场时就是Fock算符)的本征态,你如果硬去算它们之间的 electronic coupling (<S0|F|S1>),结果总是零。

看到这里你就可以去看 Van Voorhis 那篇综述了。那篇文章告诉大家,人们早就发现,虽然BO近似会失效,但是因为量子力学的基可以随便换,我们总可以找到一组基来描述系统状态,使得这组基中每一对基的vibronic coupling尽可能的小。那么,在这组基中,我们就总能定义势能面,“交叉点”也就终于可以明确定义了。这组基就是所谓的“diabatic states”。所以你看,真正交叉的其实是diabatic 的势能面。diabatic势能面和adiabatic势能面在基态和激发态差别很大的时候区别非常小。只有在“交叉点”附近才有较大区别。在交叉点时,两个diabatic states由于他们之间的coupling,发生了等能量的跃迁,也就是电子传递。而 adiabatic 势能面并没有交叉,相反基态和激发态的能量差异相当于两个diabatic states之间electronic coupling的两倍。

现在,你也许会问,按照以上的推导,一切电子传递似乎都可以在diabatic states的视角中描述,一切电子传递都只在 diabatic 势能面的交叉点发生,为什么还要分为non-adiabatic和adiabatic呢?其实我们必须指出,所谓diabatic states并不是能量本征态。这意味着,假如electronic coupling很大,也就是在基态和激发态能量分得很开, adiabatic 势能面可以良好定义的情况下,所谓的“diabatic states”将完全不能描述真实的物理状态!——因为这时体系必须满足量子绝热定理,一个良好的近似状态必定是不含时薛定谔方程的本征态(能量本征态),或与之足够接近。也就是说,在经过交叉点之前,已经发生了部分的电子传递(比如说可能达到了“桥”的某一位置了),但这并未在diabatic picture中体现,所以按照diabatic picture进行的计算,很可能就出现了偏差。也许有人会说,能量本征态也并不比叠加态更“真实”,而diabatic state又可以看做基态和激发态的叠加态,怎么就“不能描述真实的物理状态”了呢?我的理解是,单纯的看某一个叠加态,它确实可能在假象的“圆锥交叉点(conical intersection)”处出现。但是,假如系统的初态就是基态和激发他差别很大的热平衡态开始的,那么几乎可以肯定,初态为50%基态,50%激发态叠加态的概率几乎为零,几乎所有的初态都是基态。那么,如果体系满足量子绝热定理(也就是adiabatic electron transfer),初态如果是基态,在“圆锥交叉点”的态,毫无疑问也应该是基态,而不是50%的基态与50%激发态叠加形成的一个diabatic state。

所以,non-adibatic过程确实是diabatic states的electronic coupling比较弱的过程。但是要注意,"adiabatic"这个词本身并不是“强耦合”的意思。实际上假如只考虑“vibronic coupling”,adiabatic 描述反而要求这种耦合很弱。

还有,adibatic 虽然是“缓变”,但是adiabatic electron transfer的速率一般比non-adiabatic electron transfer快很多。因为实验测定的电子传递速率是一个宏观量。对于单次electron transfer,虽然adibatic ET比non-adiabatic ET的耗时要多,但是adiabatic ET的“成功率”比non-adiabatic ET要高非常多。也就是说,虽然adiabatic ET微观很慢,但是很容易发生,non-adiabatic虽然“快”,但是极其罕见。总的效果就是,adiabatic ET在宏观上的速率要快很多。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
wbn + 5 太谢谢了!

查看全部评分 View all ratings

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 02:23 , Processed in 0.172726 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list