请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 141|回复: 15

[Multiwfn使用咨询] 求助使用Multiwfn计算过渡金属氧化物的多中心键级卡住

[复制链接]

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

发表于 2020-9-10 11:25:33 | 显示全部楼层 |阅读模式
使用Multiwfn计算过度金属氧化物团簇的多中心键级,求问一直卡在这步是怎么回事?已经好几个小时了。。。

计算八中心键级的求助

计算八中心键级的求助

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-10 12:13:52 | 显示全部楼层
上传压缩后的输入文件
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-10 14:31:19 | 显示全部楼层
sobereva 发表于 2020-9-10 12:13
上传压缩后的输入文件

sob老师,附件即为我的Multiwfn的输入文件,还请您帮忙看下

input.zip

613.34 KB, 下载次数: 4

输入文件

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-10 16:06:12 | 显示全部楼层
纯粹是基组用得太大,而多中心键级耗时本身又高
本身基于NAO的分析对于基组就不敏感,Cu用SDD,O用6-31+G*就足够了
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-10 17:07:27 | 显示全部楼层
sobereva 发表于 2020-9-10 16:06
纯粹是基组用得太大,而多中心键级耗时本身又高
本身基于NAO的分析对于基组就不敏感,Cu用SDD,O用6-31+G* ...

好的,谢谢sob老师,我已换用小基组重新计算了一下,但是这个结果我很难用其来说明我体系具有芳香性,请问“Above results in normlalized form”这行指代什么结果?

多中心键级计算结果

多中心键级计算结果

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-11 07:45:06 | 显示全部楼层
我不知道你说的小基组具体指什么。我之前自己用SDD结合6-31+G*算了下你的体系,也是基于NAO算的,归一化的多中心键级我记得是0.5多一点

另外,我刚刚发明了一个超快的计算多中心键级的方法,结果和原本的代码给出来的完全相同,而耗时低了N个数量级!我用新的代码基于你上传的log文件算了一下,结果是这样的:
  1. The multicenter bond order from alpha density matrix:    0.0000064
  2. The multicenter bond order from beta density matrix:     0.0074114
  3. The sum of multicenter bond order from alpha and beta parts:        0.0074178
  4. Above result in normalized form:    0.5417319
  5. The multicenter bond order from mixed alpha&beta density matrix:    0.0000671
  6. Above result in normalized form:    0.3008232
复制代码

单单从归一化的多中心键级上看,算是有8中心共轭作用的。

含有新的计算多中心键级的代码的Multiwfn 3.8(dev)我这两天就会传到Multiwfn官网上

北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-11 09:02:32 | 显示全部楼层
sobereva 发表于 2020-9-11 07:45
我不知道你说的小基组具体指什么。我之前自己用SDD结合6-31+G*算了下你的体系,也是基于NAO算的,归一化的 ...

噢噢,好的,感谢sob老师,“Above result in normalized form”这行显示的数值就是指归一化的多中心键级吧?

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-11 09:56:14 | 显示全部楼层
SS123 发表于 2020-9-11 09:02
噢噢,好的,感谢sob老师,“Above result in normalized form”这行显示的数值就是指归一化的多中心键级 ...

北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-11 10:35:41 | 显示全部楼层

sob老师,我用sdd结合6-31+G*算了多中心键级,归一化的多中心键级显示仍是0.3526854,我上传了输入和输出文件的附件,麻烦您看下,问题出在了哪儿?

output.zip

166.43 KB, 下载次数: 0

输出文件

input.gjf

665 Bytes, 下载次数: 1

输入文件

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-12 04:29:08 | 显示全部楼层
SS123 发表于 2020-9-11 10:35
sob老师,我用sdd结合6-31+G*算了多中心键级,归一化的多中心键级显示仍是0.3526854,我上传了输入和输出 ...

你都没读赝势啊

仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-12 09:40:34 | 显示全部楼层
sobereva 发表于 2020-9-12 04:29
你都没读赝势啊

仔细看

好的,sob老师,一时疏忽,犯了低级错误

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-13 16:50:08 | 显示全部楼层
sobereva 发表于 2020-9-10 16:06
纯粹是基组用得太大,而多中心键级耗时本身又高
本身基于NAO的分析对于基组就不敏感,Cu用SDD,O用6-31+G* ...

sob老师,我换用两个不同的基组分别基于NAO分析了该体系的多中心键级,结果Cu用SDD,O用6-31+G*的归一化键级为0.5322336;全用ma-TZVP基组时得出的归一化键级却只有0.1790230。为什么差别会这么大,ma-TZVP基组不适用于计算多中心键级么??附件为我上传的这两个不同基组的输入和输出文件,还请您指点。

ma-TZVP-output.zip

221.59 KB, 下载次数: 2

ma-TZVP-input.gjf

3.58 KB, 下载次数: 1

sdd-output.zip

151.83 KB, 下载次数: 0

sdd-input.gjf

694 Bytes, 下载次数: 2

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 2020-9-14 06:50:02 | 显示全部楼层
SS123 发表于 2020-9-13 16:50
sob老师,我换用两个不同的基组分别基于NAO分析了该体系的多中心键级,结果Cu用SDD,O用6-31+G*的归一化 ...

ma-TZVP的情况自洽场都没有收敛。要注意G16 C.01碰见自洽场没收敛时不会报错,所以需要自己留意。
一方面当前的ma-TZVP波函数质量可能不佳,另一方面也有收敛到了不稳定波函数的可能,做stable测一下才清楚。

我之前说的有新代码的Multiwfn 3.8(dev)我前两天已经上传了。
鉴于处理你最开始上传的文件,以及SDD结合6-31+G*的结果都合理,只能认为ma-TZVP目前的波函数有一定问题。
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

30

帖子

0

威望

121

eV
积分
151

Level 3 能力者

 楼主| 发表于 2020-9-14 08:58:53 | 显示全部楼层
sobereva 发表于 2020-9-14 06:50
ma-TZVP的情况自洽场都没有收敛。要注意G16 C.01碰见自洽场没收敛时不会报错,所以需要自己留意。
一方 ...

好的,感谢sob老师,如果我想对该体系进行GIMIC分析,那最好选用什么方法基组呢?我看博文里说开壳层体系选用CCSD最稳妥,基组应当选3-zeta加极化的。

2万

帖子

25

威望

3万

eV
积分
58600

管理员

公社社长

发表于 6 天前 | 显示全部楼层
SS123 发表于 2020-9-14 08:58
好的,感谢sob老师,如果我想对该体系进行GIMIC分析,那最好选用什么方法基组呢?我看博文里说开壳层体系 ...

用适合配合物的泛函,结合def2-TZVP就够了
算不动的话过渡金属用SDD,配体用6-311G*,耗时不太高结果又靠谱
北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2020-9-22 15:06 , Processed in 0.250323 second(s), 28 queries .

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