(3)DZVP-MOLOPT-GTH是带有指数很小的高度弥散的高斯函数的,带SR(short-range)的DZVP-MOLOPT-SR-GTH...
培训里的各种好几百原子体系的AIMD的例子在96核双路7R32上用PBE/DZVP-MOLOPT-SR-GTH跑,平均每步的...
CP2K默认的DZVP-MOLOPT-SR-GTH基组适合大部分分子体系,处理金属体系时建议换用更密集的MOLOPT基组。比如研究铂催化剂表面吸附,将基组从DZVP升级到TZVP级别,电子态描述更精细,但计算耗时可能增加2-3倍。周期性体系中平面波截断能需特别关注,500Ry适用于多数半导体材料,处理过渡金属氧化物等强关联体系时建议提升至800Ry...
BASIS_SET选择DZVP-MOLOPT-GTH(固体用DZVP-MOLOPT-SR-GTH,这里是非周期性分子,不需要Short Range) 以上为非周期性边界条件的计算相对常规的计算需要修改的全部内容。 三、乙醇分子的振动频率计算 通过上一步的结构优化,我们获得了一个可以用于频率计算的乙醇结构,接下来我们计算乙醇分子的振动频率和强度。 我们要用到...
若无特殊指定,这个系列此后所有文章出现的基组均为DZVP-MOLOPT-SR-GTH,赝势均用GTH-PBE,并使用PBE泛函。使用这套计算等级可以应对绝大多数常规的固态计算,并且有着可观的精度。 可以在这个网站里下载到MOLOPT基组和GTH赝势。 2.fcc-Si的结构 可以在Materials Project里获得Si的结构,打开网站并登录后,直接搜索Si,...
启动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并将之载入,然后输入 ...
对于TDDFT目的用TZVP-GTH赝势基组就不错,精度足够,而且也不很贵。再好一点还可以用TZV2P-GTH。用CP2K里常用的DZVP-MOLOPT-SR-GTH赝势基组算TDDFT时还能比TZVP-GTH更便宜一些,但是由于此基组收缩度太高,导致之后用Multiwfn做很多分析的时候在计算格点数据时会非常昂贵,所以本文用TZVP-GTH。
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 ...
之后给出杂化泛函计算实际体系的重要例子,包括PBE0/DZVP-MOLOPT-SR-GTH结合ADMM辅助基组admm-dzp计算MOF-5晶体、B3LYP和HSE06结合DZVP-MOLOPT-SR-GTH(利用ADMM)和6-31G*(不利用ADMM)以超胞方式计算Si晶体的gap。最后介绍CP2K中Hartree-Fock的计算、适合搭配的赝势和基组,并示例用HF/DZVP-GTH计算64个水的盒子。