MOLOPT是CP2K的特制基组，一般与GTH赝势联用。DZVP基组的精度相当于6-31G*（据https://www.cp2k.org...
定义计算里出现的元素 ELEMENT Si BASIS_SET DZVP-MOLOPT-SR-GTH !与BASIS_MOLOPT文件对应 POTENTIAL GTH-PBE-q4 !与GTH_POTENTIALS文件对应 &END KIND BASIS_SET后写什么可以在BASIS_MOLOPT文件里找到: 划线部分为基组名字,末尾q4表示要用价带电子数为4的赝势 同理,POTENTIAL后写什么可以在GTH_POTENTIALS文件里...
(备注:后来发现这里的基组最好不用-SR-的版本,而是直接用DZVP-MOLOPT-GTH,根据这个网址,SR(Short Range)更加适合于固体材料,对于分子的计算则可以不用SR。) &KIND H ELEMENT H BASIS_SET DZVP-MOLOPT-SR-GTH POTENTIAL GTH-PBE-q1 &END KIND &KIND O ELEMENT O BASIS_SET DZVP-MOLOPT-SR-GTH POTENTIAL GT...
DZVP-MOLOPT-SR-GTH太小，改成TZVP-MOLOPT-GTH 你既然不用杂化泛函，算这个还不如考虑k点，超胞耗时...
errorThe requested basis set <DZVP-MOLOPT-GTH> for element <MG> was not found in the basis set files <BASIS_MOLOPT> For cp2k input file, one may need to change some parameters. &FORCE_EVAL METHOD QMMM &PRINT &FORCES ON &END FORCES ...
有些人对CP2K使用的基组可谓近乎一无所知,就随便网上看几个零碎的(往往还是不适合自己情况或者有严重误导性的)例子,拿比如TZVP-GTH、DZVP-MOLOPT-SR-GTH这种BSSE挺大的基组算弱相互作用能,也不知道counterpoise是何物,发现算出来结果差(一般会严重高估)就直接赖CP2K算不准,这些人太令人无语了。我不得不说,参加...
尝试利用CP2K的OT自洽计算方法对八个原子的体心立方Si晶胞进行能量和每个原子受力的静态计算。启动CP2K计算需准备基组和赝势文件，通常选用DZVP-MOLOPT-SR-GTH基组和GTH-PBE赝势，适用于常规固态计算并具有较好的精度。文件可在指定网站获取。从Materials Project网站获取Si的结构信息。选择cif格式下载较大的...
对于TDDFT目的用TZVP-GTH赝势基组就不错,精度足够,而且也不很贵。再好一点还可以用TZV2P-GTH。用CP2K里常用的DZVP-MOLOPT-SR-GTH赝势基组算TDDFT时还能比TZVP-GTH更便宜一些,但是由于此基组收缩度太高,导致之后用Multiwfn做很多分析的时候在计算格点数据时会非常昂贵,所以本文用TZVP-GTH。
有这种错误印象的人,要么是被以讹传讹,要么是缺乏CP2K最基本常识导致了严重错觉(诸如用区区DZVP-MOLOPT-SR-GTH这种低档次基组算结合能却还不知道要考虑counterpoise校正,或者只用很小晶胞却不知道考虑k点就去算晶体属性之类,这能算准就怪了。这和Gaussian用户只用6-31G*算相互作用能、VASP只用很低的截断能算...
test.inp内容(来自计算化学公社论坛,2*2*2金刚石超胞做PBE/DZVP-MOLOPT-SR-GTH单点计算的输入文件) #Generated by Multiwfn &GLOBAL PROJECT test PRINT_LEVEL LOW RUN_TYPE ENERGY &END GLOBAL &FORCE_EVAL METHOD Quickstep &SUBSYS &CELL A 7.13358000 0.00000000 0.00000000 ...