计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9271|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] MECP两个态能量相同的条件

[复制链接 Copy URL]

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

看到关于MECP的blog,“设两个态波函数分别为|1>和|2>,体系哈密顿算符为H,则这样的两态模型构成2*2哈密顿矩阵,其四个矩阵元包括两个对角元H11=<1|H|1>、H22=<2|H|2>,以及两个等价的耦合矩阵元H12=<1|H|2>=H21=<2|H|1>。让这两个态的能量E1、E2相同需要同时满足两个条件:(1) H11=H22 (2) H12=0。” 为什么H12要等于0,麻烦老师指导一下。

176

帖子

0

威望

1452

eV
积分
1628

Level 5 (御坂)

7#
发表于 Post on 2017-2-18 21:34:29 | 只看该作者 Only view this author
卡开发发 发表于 2016-8-30 23:33
唉。。。我重新找了下资料,发觉原本的推导涉及到的条件和我想的不一样,这个你可以在Landau & Lifshitz ...

等于0了不就是不考虑旋轨耦合了?

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
卡开发发 + 1

查看全部评分 View all ratings

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

6#
发表于 Post on 2016-8-30 23:33:38 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2016-8-31 10:24 编辑

唉。。。我重新找了下资料,发觉原本的推导涉及到的条件和我想的不一样,这个你可以在Landau & Lifshitz的QM书上的Chapter11的$79小节找到(中英文版本一致)。

书上写的很详细,在这只是简要做些说明:

实际上|1>和|2>并非H的本征态(我之前考虑的差异就在这了),而是H0的本征态,此时假想两者本征值比较接近但并不相等,但只要对系统的原子核稍作移动δr,系统的势能项的增量H'就会变成dH/dr*δr,这两个本征值就会变为E1和E2。通过简并微扰的方法,令|ψ>=c1|1>+c2|2>H|ψ>=(H0+H')|ψ>=E|ψ> --> (H-E)|ψ>=0

依次左乘<1|和<2|,令H11=<1|H|1>、H22=<2|H|2>,H12=<1|H|2>=H21=<2|H|1>,分别得到
c1(H11-E)+c2H12=0
c1H12+c2(H22-E)=0

求解后有
E1=1/2(H11+H22)-1/2√[(H11-H22)^2+4H12^2]
E2=1/2(H11+H22)+1/2√[(H11-H22)^2+4H12^2]
仅当H12=0的时候E1=E2.

不过只考察交叉点处的话,之前的推导也没什么问题。



日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2016-8-30 15:17:02 | 只看该作者 Only view this author
卡开发发 发表于 2016-8-30 12:42
|1>和|2>应当是正交归一的,
H12==E2
H21==E1=E1(Hermite性)

谢谢大神

588

帖子

0

威望

3310

eV
积分
3898

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2016-8-30 15:14:11 | 只看该作者 Only view this author
谢谢sob老师

6万

帖子

99

威望

5万

eV
积分
120201

管理员

公社社长

3#
发表于 Post on 2016-8-30 15:08:19 | 只看该作者 Only view this author
你去求那个2*2哈密顿矩阵的本征值,也就是两个绝热态的能量E1和E2,会有这样的公式

显然必须让H12为0才有可能让E1=E2。否则光是H11=H22,即透热态能量相等,再考虑它们的耦合H12的话,就又不相等了。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

3754

帖子

3

威望

1万

eV
积分
19681

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2016-8-30 12:42:01 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2016-8-30 12:50 编辑

|1>和|2>应当是正交归一的,
H12=<1|H|2>=E2<1|2>
H21=<2|H|1>=E1<2|1>=E1<1|2>(Hermite性)
上下两式相减,H12=H21 => (E2-E1)<1|2>=0,考虑条件E1=E2,有两种情形:
(1)<1|2>=1,此时|1>=|2>,H11=H22=H12=H21;
(2)<1|2>=0,E1=E2但|1>和|2>不同,所以此时|1>和|2>是简并态,H12=H21=E1(or E2)*<1|2>=0。
写法不是完全严谨,仅供参考吧。不过说起来直觉上,本征值相同的两个态,不是相同大概也就是简并了。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

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

GMT+8, 2025-8-18 03:28 , Processed in 0.222699 second(s), 25 queries , Gzip On.

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