若无特殊指定,这个系列此后所有文章出现的基组均为DZVP-MOLOPT-SR-GTH,赝势均用GTH-PBE,并使用PBE泛函。使用这套计算等级可以应对绝大多数常规的固态计算,并且有着可观的精度。 可以在这个网站里下载到MOLOPT基组和GTH赝势。 2.fcc-Si的结构 可以在Materials Project里获得Si的结构,打开网站并登录后,直接搜索Si,...
ELEMENT Li BASIS_SET DZVP-MOLOPT-SR-GTH POTENTIAL GTH-PBE-q3 &END KIND &KIND F ELEMENT F BASIS_SET DZVP-MOLOPT-SR-GTH POTENTIAL GTH-PBE-q7 &END KIND &CELL ABC 7.82371 7.82371 7.82371 &END CELL &TOPOLOGY COORD_FILE_FORMAT XYZ COORD_FILE_NAME ini_opted.xyz &END TOPOLOGY &END SUBSYS &...
DZVP-MOLOPT-SR-GTH太小，改成TZVP-MOLOPT-GTH 你既然不用杂化泛函，算这个还不如考虑k点，超胞耗时...
(3)DZVP-MOLOPT-GTH是带有指数很小的高度弥散的高斯函数的，带SR（short-range）的DZVP-MOLOPT-SR-GTH...
启动CP2K计算需准备基组和赝势文件，通常选用DZVP-MOLOPT-SR-GTH基组和GTH-PBE赝势，适用于常规固态计算并具有较好的精度。文件可在指定网站获取。从Materials Project网站获取Si的结构信息。选择cif格式下载较大的晶体结构，适合CP2K处理。创建Si_SC.inp文件，作为CP2K执行任务的输入。在文件中进行全局设置...
BaTiO3晶体的cif文件是本文文件包里的BaTiO3.cif。用Multiwfn创建其4*4*4超胞在PBE/DZVP-MOLOPT-SR-GTH下做单点的CP2K计算的输入文件。当前的超胞在每个方向有16埃,足够大了。CP2K计算后得到molden文件,在里面加入恰当的[Cell]和[Nval]字段后,启动Multiwfn并将之载入,然后输入 ...
在当前用的赝势下,各种元素的有效核电荷数是多少可以去看基组定义,比如在CP2K的data目录下的BASIS_MOLOPT文件中可以看到有一行Cl DZVP-MOLOPT-SR-GTH DZVP-MOLOPT-SR-GTH-q7,其中q后面的数字就是有效核电荷数(价电子数)。在目前版本Multiwfn产生的使用GTH赝势的输入文件中,直接看BASIS_SET后面的基组名的q后面的...
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 ...
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 ...
and changeDZVP-MOLOPT-GTHtoDZVP-MOLOPT-SR-GTH, and qmmm-cp2k-active = true qmmm-cp2k-qmgroup = QMatoms qmmm-cp2k-qmmethod = PBE qmmm-cp2k-qmcharge = -3 qmmm-cp2k-qmmultiplicity = 1 to qmmm-cp2k-active = true qmmm-cp2k-qmgroup = QMatoms ...