计算化学公社

标题: ORCA 建议部分使用相对论计算代替赝势是否合理? [打印本页]

作者
Author:
liyuanhe211    时间: 2015-10-27 15:18
标题: ORCA 建议部分使用相对论计算代替赝势是否合理?
本帖最后由 liyuanhe211 于 2015-10-27 15:31 编辑

在ORCA Input Library 中,对ECP的使用(https://sites.google.com/site/or ... ive-core-potentials)有如下解释:

General recommendations
- If your molecule does not contain elements heavier than Kr, then you typically should not use an ECP. Stick with an all-electron approach so as not to not sacrifice accuracy for marginal speed-up.
- If your molecule contains one heavy element (heavier than Kr), e.g. in a typical transition metal complex, there is probably not a big advantage to using an ECP. You may get more accurate results and almost as fast, using an all-electron relativistic approximation instead.
- If the molecule contains many heavy elements (heaver than Kr), then ECPs have their greatest advantage and efficient geometry optimizations can be performed. Single-point calculations using an all-electron scalar relativistic approach might then be a good idea for even more accurate energies or property calculations.
- There are also cases where one is dealing with so heavy elements that one runs into numerical problems when using an all-electron scalar relativistic method (ZORA, DKH). Using an ECP might then be preferred.

即:
第一过渡系用全电子基组(用ECP节省不了多少时间)
有一个2/3过渡系原子,推荐使用全电子相对论近似(用ECP不会有很大优势)
有很多2/3过渡系原子,用ECP计算,最后用全电子标量相对论计算可以提高精度。
                                如果标量相对论效应有数值问题,则换用ECP

这样的说法是否合理?
考虑ORCA建议对相对论计算使用最高的积分格点Grid7,并自行增加积分精度至IntAcc 10~20,对有一个2/3过渡系原子的情况(如Rh配合物)是否会因积分格点增加导致计算时间大量增加?(而不是“not a big advantage to using an ECP”)

-------

另外相对论计算方法ZORA和DKH2有什么区别,如何选择?



作者
Author:
sobereva    时间: 2015-10-27 15:47
说法基本合理。不过对于第一周期过渡金属,用赝势也还是能省不少时间的。
用全电子标量相对论计算本身比赝势耗时高,再加上又得用更好积分格点,所以耗时更高。
ZORA和DKH2都是常用的把四分量相对论计算方程近似变换成二分量的做法。DKH2远比原始的ZORA更好,不过ORCA用的是修正的ZORA。实际计算就用DKH2就行了。
作者
Author:
liyuanhe211    时间: 2015-10-27 16:28
sobereva 发表于 2015-10-27 15:47
说法基本合理。不过对于第一周期过渡金属,用赝势也还是能省不少时间的。
用全电子标量相对论计算本身比赝 ...

明白了。十分感谢~~
作者
Author:
清微    时间: 2018-1-5 22:56
老师,在高斯中DKH2和DKH是同义词,但DKH和DK含义完全不同,据我的理解,DK是把Dirac方程转化为二分量,属于A2C方法,但DKH属于纯标量相对论方法,一字之差,含义天壤之别。
作者
Author:
sobereva    时间: 2018-1-6 07:42
清微 发表于 2018-1-5 22:56
老师,在高斯中DKH2和DKH是同义词,但DKH和DK含义完全不同,据我的理解,DK是把Dirac方程转化为二分量,属 ...


一般也不那么较真,诸如Dirac程序,写.DKH2即做二分量相对论计算,做标量的时候只不过要再写.SPINFREE

作者
Author:
rayinyin    时间: 2018-5-4 10:46
楼上的讨论是否主要针对DFT,如果是波函数的计算,还有相应grid的设置的必要吗?
作者
Author:
sobereva    时间: 2018-5-4 20:41
rayinyin 发表于 2018-5-4 10:46
楼上的讨论是否主要针对DFT,如果是波函数的计算,还有相应grid的设置的必要吗?

对于波函数方法的计算,grid设置完全不影响结果,除非是那种利用RI加速的情况
作者
Author:
rayinyin    时间: 2018-5-7 09:50
rijcosx方法需要设置一个gridx参数,如果是RI-JK加速于波函数方法的话,还需要设置grid相关的参数否?
作者
Author:
sobereva    时间: 2018-5-7 10:45
rayinyin 发表于 2018-5-7 09:50
rijcosx方法需要设置一个gridx参数,如果是RI-JK加速于波函数方法的话,还需要设置grid相关的参数否?

不需要




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3