MS里Dmol3优化不同自旋态的相同结构能量差距很大是什么问题？

LZLLLL · 发表于 Post on 2020-7-29 12:41:28

对金属铁配合物，用了两种不同的自旋多重度三重态和单重态，但是发现输出文件的能量相差大概4 Hartree, 这是什么原因？
我打开输出文件发现单重态和三重态里spin polarization 里相差很大，这是什么原因？

zjxitcc · 发表于 Post on 2020-7-29 13:19:02

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

LZLLLL · 发表于 Post on 2020-7-29 13:43:01

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

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

zjxitcc · 发表于 Post on 2020-7-29 14:03:03

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

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

LZLLLL · 发表于 Post on 2020-7-29 17:11:08

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

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

sobereva · 发表于 Post on 2020-7-30 00:34:11

用Dmol3莫名其妙。图速度，用ORCA不仅更快，功能还强大得多得多，还免费，还能结合Multiwfn做各种后处理分析

卡开发发 · 发表于 Post on 2020-7-30 08:39:57

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

LZLLLL · 发表于 Post on 2020-7-31 09:06:28

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

好的sob老师，下次不用了

LZLLLL · 发表于 Post on 2020-7-31 09:07:13

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[其它量化程序] MS里Dmol3优化不同自旋态的相同结构能量差距很大是什么问题？

浏览过的版块