ORCA 建议部分使用相对论计算代替赝势是否合理？

liyuanhe211

在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有什么区别，如何选择？

sobereva

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

liyuanhe211

sobereva 发表于 2015-10-27 15:47
说法基本合理。不过对于第一周期过渡金属，用赝势也还是能省不少时间的。
用全电子标量相对论计算本身比赝 ...

明白了。十分感谢~~

清微

老师，在高斯中DKH2和DKH是同义词，但DKH和DK含义完全不同，据我的理解，DK是把Dirac方程转化为二分量，属于A2C方法，但DKH属于纯标量相对论方法，一字之差，含义天壤之别。

sobereva

清微 发表于 2018-1-5 22:56
老师，在高斯中DKH2和DKH是同义词，但DKH和DK含义完全不同，据我的理解，DK是把Dirac方程转化为二分量，属 ...

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

rayinyin

楼上的讨论是否主要针对DFT，如果是波函数的计算，还有相应grid的设置的必要吗？

sobereva

rayinyin 发表于 2018-5-4 10:46
楼上的讨论是否主要针对DFT，如果是波函数的计算，还有相应grid的设置的必要吗？

对于波函数方法的计算，grid设置完全不影响结果，除非是那种利用RI加速的情况

rayinyin

rijcosx方法需要设置一个gridx参数，如果是RI-JK加速于波函数方法的话，还需要设置grid相关的参数否？

sobereva

rayinyin 发表于 2018-5-7 09:50
rijcosx方法需要设置一个gridx参数，如果是RI-JK加速于波函数方法的话，还需要设置grid相关的参数否？

不需要