求助一系列含不同金属的配合物电极电势计算与对比相关问题

树独Comche

老师们好，最近在做一系列分别包含主族金属、过渡金属和镧系金属的150原子以上的配合物的电极电势的对比。根据sob老师的方法建立了一个热循环的方法来计算其氧化势,其中的溶解自由能的计算按照http://sobereva.com/327的方法计算。计算是在b3lyp/def2svp条件用高斯做结构优化，用双杂化泛函用orca算气相下的单点（之前也尝试在g09下用def2tzvp水平算单点，但担心精度不够），再利用m062x/6-31g*用高斯算溶解自由能。其中由于含镧系金属的结构无法用b3lyp/def2svp得到，所以我们采用大核赝势计算后用b3lyp/def2svp计算不同重态确定其能量最低的重态。后算溶解自由能时对Gd使用SDD。但这样就存在不同金属使用的计算水平不同，无法对一系列分别含主族金属元素、过渡金属和镧系金属的结构相似的配合物的电势进行对比，请问这个问题可以这么解决吗？

wzkchem5

不同元素的基组只需要相似（zeta数、弥散、极化、赝势、相对论）即可，不需要严格一致。原因很简单，即便都是叫一个名字的基组，也可能是不同人开发的，比如甲给H-Ar的元素开发了某基组，乙给K-Kr开发了另一基组，然后丙发现两个基组的大小类似，就给它们起了同一个名字。你去看基组需要引什么文献也会发现，很多时候不同周期元素需要引的文献不一样。所以你对不同元素用同样的基组，引入的理论不一致性可能不比你对不同元素用不同的基组那样更小。
如果实在担心基组不统一导致问题，就把所有金属元素的基组都选大一些，因为当基组趋于完备时，所有基组都是等价的。

树独Comche

wzkchem5 发表于 2022-7-14 18:02
不同元素的基组只需要相似（zeta数、弥散、极化、赝势、相对论）即可，不需要严格一致。原因很简单，即便都 ...

大概懂了，谢谢老师！