计算化学公社

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

[量化理论] 跃迁偶极矩的大小是否与原点选取有关?

[复制链接 Copy URL]

230

帖子

0

威望

935

eV
积分
1165

Level 4 (黑子)

跃迁偶极矩的定义是<i|r|j>。我觉得理论上它应该是一个定义良好的物理量而不应该与坐标原点的选取有关。但是假如用高斯分布来模拟一维的电子波函数,结果却似乎与原点位置有关,比如我用如下的matlab代码来模拟psi_i(x)*x*psi_j(x):
  1. sigma=1;
  2. distance=6;
  3. center=0;
  4. x=-10:0.1:10;
  5. y1=-normpdf(x,center-distance/2,sigma);title(['center at x=',num2str(center)]);
  6. y2=normpdf(x,center+distance/2,sigma);
  7. y=-normpdf(x,center-distance/2,sigma).*x.*normpdf(x,center+distance/2,sigma);
  8. subplot(311);plot(x,y1);ylabel('\psi_i(x)');
  9. subplot(312);plot(x,y2);ylabel('\psi_j(x)');
  10. subplot(313);plot(x,y);ylabel('\psi_i(x)x\psi_j(x)');
  11. xlabel('x');
复制代码
normpdf是归一化的高斯分布。其中psi_i<0,psi_j>0。这样可以保证<i|j>=0,满足正交归一条件。
然后计算psi_i(x)*x*psi_j(x)。可以看到当center=0的时候,被积函数是奇函数,跃迁偶极矩为0。但是当center不为零(相当于坐标平移了一下或者另取零点的时候)的时候,被积函数则不是奇函数,所模拟的跃迁偶极矩积分也不为零。


这个结果我一直没想明白。不知道真正的波函数里有什么机制可以保证跃迁偶极矩不随原点变化而变化呢?






评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
sobereva + 3 有意义的讨论

查看全部评分 View all ratings

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

2#
发表于 Post on 2015-6-7 15:13:06 | 只看该作者 Only view this author
有个关键问题是这两个函数并不正交,只不过是相位相反而已。你做一个这两个函数乘积的图,就会发现只是在负值区域有个峰,积分显然不为0。

对于真正正交的情况,不依赖于原点是肯定满足的。比如原先是<1|x|2>,函数原点位移R(等价于x位移R),变成<1|x-R|2>=<1|x|2>-R*<1|2>,由于<1|2>=0,因此还是原先的<1|x|2>。
北京科音自然科学研究中心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

296

帖子

1

威望

2588

eV
积分
2904

科音成员

3#
发表于 Post on 2015-6-7 23:25:50 | 只看该作者 Only view this author
给初学者的提示:楼上推导的关键是,因为R是一个常数,所以<1|R|2>=R<1|2>。
华北电力大学数理学院,理论与计算化学,团簇、表面的结构与反应机理。(招第一性原理计算,量子化学计算方向的教师、硕士/博士研究生)

230

帖子

0

威望

935

eV
积分
1165

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2015-6-7 23:58:35 | 只看该作者 Only view this author
本帖最后由 stecue 于 2015-6-8 00:03 编辑
sobereva 发表于 2015-6-7 15:13
有个关键问题是这两个函数并不正交,只不过是相位相反而已。你做一个这两个函数乘积的图,就会发现只是在负 ...

了解!Thanks!我是把那两个高斯函数想象成了狄拉克delta函数(两个不同位置的狄拉克函数根据应该是正交的吧?)……

6万

帖子

99

威望

5万

eV
积分
120062

管理员

公社社长

5#
发表于 Post on 2015-6-8 00:24:09 | 只看该作者 Only view this author
stecue 发表于 2015-6-7 23:58
了解!Thanks!我是把那两个高斯函数想象成了狄拉克delta函数(两个不同位置的狄拉克函数根据应该是正 ...

是的
北京科音自然科学研究中心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

230

帖子

0

威望

935

eV
积分
1165

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2015-6-8 03:33:21 | 只看该作者 Only view this author
本帖最后由 stecue 于 2015-6-9 04:04 编辑

重复编辑掉(不能自己删掉回复的帖子吗?)

230

帖子

0

威望

935

eV
积分
1165

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2015-6-9 04:02:52 | 只看该作者 Only view this author

对了那对于任意一个厄米算符O,|1>和|2>的正交归一性是<1|O|2>具有平移不变性,也就是<1|O|2>=<1(x-R)|O|2(x-R)>的充要条件吗?简单想了一下好像没想出来……

本版积分规则 Credits rule

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

GMT+8, 2025-8-12 21:54 , Processed in 0.625768 second(s), 24 queries , Gzip On.

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