计算化学公社

标题: 求助使用Multiwfn计算过渡金属氧化物的多中心键级卡住 [打印本页]

作者
Author: SS123 时间: 2020-9-10 11:25
标题: 求助使用Multiwfn计算过渡金属氧化物的多中心键级卡住
使用Multiwfn计算过度金属氧化物团簇的多中心键级，求问一直卡在这步是怎么回事？已经好几个小时了。。。

作者
Author: sobereva 时间: 2020-9-10 12:13
上传压缩后的输入文件

作者
Author: SS123 时间: 2020-9-10 14:31

sobereva 发表于 2020-9-10 12:13
上传压缩后的输入文件

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

作者
Author: sobereva 时间: 2020-9-10 16:06
纯粹是基组用得太大，而多中心键级耗时本身又高
本身基于NAO的分析对于基组就不敏感，Cu用SDD，O用6-31+G*就足够了

作者
Author: SS123 时间: 2020-9-10 17:07

sobereva 发表于 2020-9-10 16:06
纯粹是基组用得太大，而多中心键级耗时本身又高
本身基于NAO的分析对于基组就不敏感，Cu用SDD，O用6-31+G* ...

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

作者
Author: sobereva 时间: 2020-9-11 07:45
我不知道你说的小基组具体指什么。我之前自己用SDD结合6-31+G*算了下你的体系，也是基于NAO算的，归一化的多中心键级我记得是0.5多一点

另外，我刚刚发明了一个超快的计算多中心键级的方法，结果和原本的代码给出来的完全相同，而耗时低了N个数量级！我用新的代码基于你上传的log文件算了一下，结果是这样的：

The multicenter bond order from alpha density matrix: 0.0000064
The multicenter bond order from beta density matrix: 0.0074114
The sum of multicenter bond order from alpha and beta parts: 0.0074178
Above result in normalized form: 0.5417319
The multicenter bond order from mixed alpha&beta density matrix: 0.0000671
Above result in normalized form: 0.3008232

复制代码

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

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

作者
Author: SS123 时间: 2020-9-11 09:02

sobereva 发表于 2020-9-11 07:45
我不知道你说的小基组具体指什么。我之前自己用SDD结合6-31+G*算了下你的体系，也是基于NAO算的，归一化的 ...

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

作者
Author: sobereva 时间: 2020-9-11 09:56

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

对

作者
Author: SS123 时间: 2020-9-11 10:35

sobereva 发表于 2020-9-11 09:56
对

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

作者
Author: sobereva 时间: 2020-9-12 04:29

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

你都没读赝势啊

仔细看
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

作者
Author: SS123 时间: 2020-9-12 09:40

sobereva 发表于 2020-9-12 04:29
你都没读赝势啊

仔细看

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

作者
Author: SS123 时间: 2020-9-13 16:50

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基组不适用于计算多中心键级么？？附件为我上传的这两个不同基组的输入和输出文件，还请您指点。

作者
Author: sobereva 时间: 2020-9-14 06:50

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目前的波函数有一定问题。

作者
Author: SS123 时间: 2020-9-14 08:58

sobereva 发表于 2020-9-14 06:50
ma-TZVP的情况自洽场都没有收敛。要注意G16 C.01碰见自洽场没收敛时不会报错，所以需要自己留意。
一方 ...

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

作者
Author: sobereva 时间: 2020-9-16 09:37

SS123 发表于 2020-9-14 08:58
好的，感谢sob老师，如果我想对该体系进行GIMIC分析，那最好选用什么方法基组呢？我看博文里说开壳层体系 ...

用适合配合物的泛函，结合def2-TZVP就够了
算不动的话过渡金属用SDD，配体用6-311G*，耗时不太高结果又靠谱

作者
Author: SS123 时间: 2020-9-16 10:15

sobereva 发表于 2020-9-16 09:37
用适合配合物的泛函，结合def2-TZVP就够了
算不动的话过渡金属用SDD，配体用6-311G*，耗时不太高结果又 ...

噢噢，好的，感谢sob老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)