计算化学公社

标题: MS里Dmol3优化不同自旋态的相同结构能量差距很大是什么问题？ [打印本页]

作者
Author: LZLLLL 时间: 2020-7-29 12:41
标题: MS里Dmol3优化不同自旋态的相同结构能量差距很大是什么问题？
对金属铁配合物，用了两种不同的自旋多重度三重态和单重态，但是发现输出文件的能量相差大概4 Hartree, 这是什么原因？
我打开输出文件发现单重态和三重态里spin polarization 里相差很大，这是什么原因？

作者
Author: zjxitcc 时间: 2020-7-29 13:19
这种分子体系最好的办法是用高斯算，还能检验波函数稳定性，特殊情况下可以构造复杂的波函数初猜。用Dmol3直接这样算，结果波函数稳定与否都不知道，比较能量高低意义不大。

作者
Author: LZLLLL 时间: 2020-7-29 13:43

zjxitcc 发表于 2020-7-29 13:19
这种分子体系最好的办法是用高斯算，还能检验波函数稳定性，特殊情况下可以构造复杂的波函数初猜。用Dmol ...

谢谢，已经用高斯提交上任务了，那请问如果还要用Dmol3算这个结构的话，把收敛限设的严一点会不会有效果？

作者
Author: zjxitcc 时间: 2020-7-29 14:03

LZLLLL 发表于 2020-7-29 13:43
谢谢，已经用高斯提交上任务了，那请问如果还要用Dmol3算这个结构的话，把收敛限设的严一点会不会有效果 ...

没有用的。3d过渡金属化合物的波函数稳定性是首要的。波函数稳定性与收敛限没有直接关系。

作者
Author: LZLLLL 时间: 2020-7-29 17:11

zjxitcc 发表于 2020-7-29 14:03
没有用的。3d过渡金属化合物的波函数稳定性是首要的。波函数稳定性与收敛限没有直接关系。

好的老师，谢谢您，那我还是老老实实用高斯吧

作者
Author: sobereva 时间: 2020-7-30 00:34
用Dmol3莫名其妙。图速度，用ORCA不仅更快，功能还强大得多得多，还免费，还能结合Multiwfn做各种后处理分析

作者
Author: 卡开发发 时间: 2020-7-30 08:39
问题大概率出在截断半径上，排查的理由：初猜给出的能量就已经差别很大了，说明原子的初始密度差别很大，因此看原子自洽部分的信息，可以看到不少原子的e_ref值很大，表示当前截断半径原子计算的能量和参考值差异比较大。通过调整截断半径到Fine后，两者差异在在1e-2Ha（我自己是用Fine精度直接算的单点，你可以自行再排查）。速度方面由计算方法决定，数值原子轨道的积分方式指前因子很大，小体系DMol3不会看得出什么优势。

作者
Author: LZLLLL 时间: 2020-7-31 09:06

sobereva 发表于 2020-7-30 00:34
用Dmol3莫名其妙。图速度，用ORCA不仅更快，功能还强大得多得多，还免费，还能结合Multiwfn做各种后处理分 ...

好的sob老师，下次不用了

作者
Author: LZLLLL 时间: 2020-7-31 09:07

卡开发发发表于 2020-7-30 08:39
问题大概率出在截断半径上，排查的理由：初猜给出的能量就已经差别很大了，说明原子的初始密度差别很大，因 ...

好的老师，我调整一下截断半径看看有没有区别

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