请选择 进入手机版 | 继续访问电脑版
第9届北京科音分子动力学与GROMACS培训班将于4月17~20日于北京举办,请点击此链接查看培训详情,欢迎参加和相互转告!

计算化学公社

 找回密码
 现在注册!
查看: 491|回复: 4

[Multiwfn功能更新] Multiwfn的轨道定域化和双正交化功能已支持从ORCA输出文件中读取Fock矩阵

[复制链接]

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

发表于 2020-10-14 07:54:52 | 显示全部楼层 |阅读模式
有很多Multiwfn用户做轨道定域化(http://sobereva.com/380)和双正交化(http://sobereva.com/448)分析的时候是基于ORCA产生的.molden格式的波函数文件来做的,但之前的Multiwfn没法直接从ORCA输出文件中读取Fock矩阵,导致没法算轨道能量。为了便民,从2020-Oct-14更新的版本开始,在Multiwfn的轨道定域化和双正交化功能中,在要求输入含有Fock矩阵的文件的那一步的时候用户可以直接写ORCA的输出文件路径,相应的ORCA输入文件里必须加上%output Print[P_Iter_F] 1 end关键词来要求ORCA每轮SCF迭代都输出Fock矩阵,Multiwfn会读取最后一轮的。计算后就可以得到定域化轨道/双正交化轨道能量了。

PS:ORCA输出的Fock矩阵很特殊,P壳层的xyz顺序和一般习俗不同,而且由于用的归一化条件特殊,F、G、H壳层的一些矩阵元需要翻转符号。Multiwfn在读取的时候都已经考虑了,得到的定域化和双正交化轨道能量和用Gaussian的情况对比过,完全一致。

评分

参与人数 2eV +8 收起 理由
cyx98 + 5 赞!
ghifi37 + 3 好功能!

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

104

帖子

0

威望

935

eV
积分
1039

Level 4 (黑子)

发表于 2020-10-14 09:01:26 | 显示全部楼层
前两天才遇到这个问题,自己都没还来得及调整轨道,社长就把问题给解决了!没见过这么贴近用户的软件和开发者,大赞!

1424

帖子

0

威望

854

eV
积分
2278

Level 5 (御坂)

发表于 2020-10-14 12:38:45 | 显示全部楼层
实际上从molden也可以读出Fock矩阵,因为molden给了轨道系数和轨道能,用F=(C^(-1))^T*epsilon*C^(-1)就出来了,这样的有效数字位数应该比从orca输出文件读要多。
类似地,从molden也可以读出overlap矩阵,方法是利用C^T*S*C=I,推出S=(C*C^T)^(-1)。
此外还有很多trick,比如让orca做Hcore初猜,NoIter,然后读molden,用(C^(-1))^T*epsilon*C^(-1)可以得到Hcore矩阵,Fock矩阵跟Hcore矩阵一减可以得到F2e+Fxc

2万

帖子

25

威望

3万

eV
积分
65392

管理员

公社社长+计算化学玩家

 楼主| 发表于 2020-10-15 05:38:25 | 显示全部楼层
wzkchem5 发表于 2020-10-14 12:38
实际上从molden也可以读出Fock矩阵,因为molden给了轨道系数和轨道能,用F=(C^(-1))^T*epsilon*C^(-1)就出 ...

牵扯到线性依赖基函数的时候会出现乱七八糟麻烦,所以我不考虑从轨道能量和组合系数反变换成Fock矩阵的做法。

重叠矩阵Multiwfn自身直接就会计算出来。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班。这些培训是计算化学快速入门以及全面系统性提升研究水平的最佳途径,培训各种相关信息见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436。3号:764390338,合计8000人,讨论范畴相同,可加入任意其一但不可都加入。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

1424

帖子

0

威望

854

eV
积分
2278

Level 5 (御坂)

发表于 2020-10-15 10:03:58 | 显示全部楼层
sobereva 发表于 2020-10-15 05:38
牵扯到线性依赖基函数的时候会出现乱七八糟麻烦,所以我不考虑从轨道能量和组合系数反变换成Fock矩阵的做 ...

有道理,线性相关的情况下一部分MO会输出成0
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )

GMT+8, 2021-4-11 18:38 , Processed in 0.211128 second(s), 25 queries .

快速回复 返回顶部 返回列表