计算化学公社

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

[综合交流] 求助双能级公式中关键态的确认

[复制链接 Copy URL]

15

帖子

0

威望

39

eV
积分
54

Level 2 能力者

最近在研究双能级公式,有几个问题,迷惑几天了,烦请各位老师解惑。

Ps:上传的4个图片是按照sob老师博文《使用Multiwfn对第一超极化率做双能级和三能级模型分析》,将SOS.txt文件导入Multiwfn中后,输入指令24-2-20以后,得到的。每个图片右侧区域的红色字体数据,是通过Gaussian,输入指令#p wB97XD/aug-cc-pvdz polar计算得到的。分析数据的目的,是想要得到一个双能级公式得到的β=(9*f*∆μ)/(E)3,4个β排列顺序,就像红色字体中Gaussian计算得到的Beta一样,得到47000>9000>3200>2900这样的顺序。

1、现在4个图片中,不知道该确认哪个是关键态,上述博文看过了,但还是思路不清晰,烦请各位老师解惑。

2、确认关键态以后,图片中左侧的Trans.dipole moment一栏中的Tot,是否就是β=(9*f*∆μ)/(E)3中的∆μ,Trans.dipole moment一栏中的E,是否就是β=(9*f*∆μ)/(E)3中的∆E?如果不是,该如何去确认∆μ∆E(查阅了相关方向硕士、博士们论文的双能级相关章节,发现有的∆μ直接解释的就是跃迁偶极矩,所以有这个疑问)


3、有一篇论文,4个化合物,对Beta分量以后,2个是Beta_Y主导,2个是Beta_X主导(主导的意思就是,比如算出来是4000,Beta_X=3980),接下来,作者使用了β=(9*f*∆μ)/(E)3。是不是说不一定非要4个化合物都是Y方向主导,也可以使用双能级公式进行排序?比如,4个化合物不管它们自己计算得到的Beta是哪个方向主导的,不限制于只关注一个方向,按照它自己实际主导的方向分析(比如1是X方向,2是Y方向,3是Z方向,4是X方向),各自使用双能级公式分析出各自的对应的β,然后再对这4个β排序即可?

4、假设只关注X方向的(输入指令#p wB97XD/aug-cc-pvdz polar,计算得到Beta是6000,Beta_X=5980),在图片中左侧的Trans.dipole moment一栏得到了30个激发态的X分量,比较30个激发态各自的X绝对值大小即可,绝对值大的,就意味着可以确认它是关键态?

5、还是说双能级公式,不止用SOS方法进行分析?还有其它方法的双能级公式?通过Gaussian,输入指令#p wB97XD/aug-cc-pvdz polar计算得到的Beta排序,和用TDDFT方法的SOS双能级公式分析得出的β排序,结论是否有必然的一致性?

6b728946234d2347cc5cffe5f6edc7fb.png (175.24 KB, 下载次数 Times of downloads: 1)

6b728946234d2347cc5cffe5f6edc7fb.png

360eca354005e308cef1de34fe408288.png (164.67 KB, 下载次数 Times of downloads: 3)

360eca354005e308cef1de34fe408288.png

833f72b54285e200d4802df3fc6b053e.png (175.29 KB, 下载次数 Times of downloads: 2)

833f72b54285e200d4802df3fc6b053e.png

a8f2a3f9547d9c78cc5ab1926df95d3a.png (174.15 KB, 下载次数 Times of downloads: 2)

a8f2a3f9547d9c78cc5ab1926df95d3a.png

6万

帖子

99

威望

5万

eV
积分
124668

管理员

公社社长

2#
发表于 Post on 2026-1-15 02:31:00 | 只看该作者 Only view this author
你引用的我的博文里的例子已经明确说了,beta evaluated by two-level model下面直接就是Multiwfn按双能级模型给你的各个激发态按双能级模型算的beta,直接看这个多好

“∆μ直接解释的就是跃迁偶极矩”这种硕博论文直接删了。没事甭看国内的硕博论文,净是错误,作者自己都没搞清楚基本情况就乱写

双能级公式是怎么从SOS简化来的,下文已经说得、推导得清楚得没法更清楚。不存在其它的双能级公式。但凡看到其它的,要么你没理解对等价性关系,要么就是错的
谈谈计算第一超极化率的双能级公式
http://sobereva.com/361http://bbs.keinsci.com/thread-5102-1-1.html

双能级模型适用的情况,其给出的beta顺序和polar按照CPHF方式算的结果能对应(除非本身不同体系间beta相差就很小,一点定量误差都可能影响顺序)。本身当考虑所有态的情况下,原理上SOS和CPHF算的结果就是等价的,双能级模型适用的情况只不过在严格的SOS结果基础上做了一定程度影响定量准确度的近似。
北京科音自然科学研究中心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

15

帖子

0

威望

39

eV
积分
54

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2026-1-15 09:14:58 | 只看该作者 Only view this author
sobereva 发表于 2026-1-15 02:31
你引用的我的博文里的例子已经明确说了,beta evaluated by two-level model下面直接就是Multiwfn按双能级 ...

老师,再麻烦您解答,那按照我发的图片中polar计算出来的Beta值:

第1个,Beta=47000;Beta_X=-11000,Beta_Y= 45000,Beta_Z= -2100;可以认为是X主导吗?因为它在三个分量数值中绝对值最大。

第2个,Beta=9000;Beta_X=6000,Beta_Y= -5600,Beta_Z= -3600;可以认为是X主导吗?因为它在三个分量数值中绝对值最大。

第3个,Beta=3200;Beta_X=1200,Beta_Y= -500,Beta_Z= -3000;可以认为是Z主导吗?因为它在三个分量数值中绝对值最大。

第4个,Beta=2900;Beta_X=0,Beta_Y= -2400,Beta_Z= -1700;可以认为是Y主导吗?因为它在三个分量数值中绝对值最大。

1、如果按照我理解的,目前4个图片中的物质,1和2是X分量主导,3是Z分量主导,4是Y分量主导,直接对它们用双能级公式,依次算出来一个数β,对它们排序,是否合理?(疑惑的点在于,它们并没有同时由一个方向主导)

2、之前有篇帖子,我咨询过您,如果旋转坐标,让3和4的Beta都平行于Z轴,得出来的结果尽可能是Z轴方向的绝对值大(体系还在计算中,暂时没出结果)。是否有必要按照调整过的坐标,再次使用TDDFT计算,得出来的beta evaluated by two-level model是否会有数值大小的变化,还是说只是XYZ方向矢量的变化(比如原来是X方向变得多,现在是Z方向变得多,但不影响数值大小)

按照老师您的建议,数据都在计算中,自身刚接触双能级公式,麻烦老师了。

本版积分规则 Credits rule

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

GMT+8, 2026-1-24 00:59 , Processed in 0.163895 second(s), 23 queries , Gzip On.

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