优化用DZVP-MOLOPT-SR-GTH基组,单点用TZVP-MOLOPT-GTH 计算能正常进行,对于计算结果,我的理解是初猜...
53 一个体系内DZVP-MOLOPT-SR-GTH和TZVP-MOLOPT-SR-GTH可以混用吗?可以,这就是混合基组 ...
虽然算大体系之间弱相互作用会比小体系间的BSSE更大,但只要用到TZV2P-MOLOPT-GTH结合counterpoise档次依然误差微乎其微。 有些人对CP2K使用的基组可谓近乎一无所知,就随便网上看几个零碎的(往往还是不适合自己情况或者有严重误导性的)例子,拿比如TZVP-GTH、DZVP-MOLOPT-SR-GTH这种BSSE挺大的基组算弱相互作用能,也...
CP2K默认的DZVP-MOLOPT-SR-GTH基组适合大部分分子体系,处理金属体系时建议换用更密集的MOLOPT基组。比如研究铂催化剂表面吸附,将基组从DZVP升级到TZVP级别,电子态描述更精细,但计算耗时可能增加2-3倍。周期性体系中平面波截断能需特别关注,500Ry适用于多数半导体材料,处理过渡金属氧化物等强关联体系时建议提升至800Ry...
对于TDDFT目的用TZVP-GTH赝势基组就不错,精度足够,而且也不很贵。再好一点还可以用TZV2P-GTH。用CP2K里常用的DZVP-MOLOPT-SR-GTH赝势基组算TDDFT时还能比TZVP-GTH更便宜一些,但是由于此基组收缩度太高,导致之后用Multiwfn做很多分析的时候在计算格点数据时会非常昂贵,所以本文用TZVP-GTH。
之后给出具体的考虑k点的杂化泛函计算例子,包括PBE0/pob-TZVP-rev2计算Li晶体,以及HSE06/TZVP-MOLOPT-GTH结合admm-dzp计算硅原胞,同时会讨论基组的弥散特征是如何对耗时产生关键影响的。特别值得一提的是,本培训里会介绍我们自行探索出的通过修改CP2K源代码实现降低HSE06几倍耗时而几乎完全不牺牲精度的重要技巧!
前文中我提到了CP2K的performance是非常好的。这篇文章我简单介绍一下为什么这么好。也就是不同于其他code的一个implementation,Quickstep重点Orbital Transformation。 CP2K作为一个周期性(或非周期性)的量子化学软件,不同于VASP, QE这些,在表示波函数和电子密度的时候使用了另外一种形式 --- Gaussian and ...
用CP2K里常用的DZVP-MOLOPT-SR-GTH赝势基组算TDDFT时还能比TZVP-GTH更便宜一些,但是由于此基组收缩度...
我也很想知道 TDDFT下用杂化泛函能算多少数目的原子 (均在PBE0 and TZVP-MOLOPT-SR-GTH 精度下)
在当前用的赝势下,各种元素的有效核电荷数是多少可以去看基组定义,比如在CP2K的data目录下的BASIS_MOLOPT文件中可以看到有一行Cl DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q7,其中q后面的数字就是有效核电荷数(价电子数)。在目前版本Multiwfn产生的使用GTH赝势的输入文件中,直接看BASIS_SET后面的基组名的q后面的...