参考下文解决。用outerSCF时MAX_SCF甭设得那么大 CP2K中遇到SCF难收敛时的解决方法 http://sobereva....
topic/cp2k/IIp-D6AA7ME 3.2 使用CP2K进行能量计算 3.2.1 OUTER_SCF 计算能量,需要设置RUN_TYPE为ENERGY;如果还要进行梯度(受力)计算,则设置为ENERGY_FORCE。 使用DFT方法进行能量计算,使用OT算法,并开启OUTER_SCF,示例如下: 1&SCF23EPS_SCF1.0E-645SCF_GUESS RESTART67MAX_SCF10089&OT T1011PRECONDITIONER F...
较大的N_DIIS值可以提高SCF收敛速度,但同时也会增加计算成本和内存需求。选择合适的N_DIIS值需要在计算效率和资源消耗之间找到平衡。 &OUTER_SCF:开启outer SCF可令难收敛的SCF收敛 MAX_SCF :意义和前面一样 EPS_SCF :意义和前面一样 &XC:对交换—相关泛函中的纯泛函进行部分设置杂化,泛函和双杂化泛函另需...
使用OT设置中的OUTER_SCF加速收敛。 一、32个水分子体系的构建 直接套用这个网页里的32个水分子结构,用VESTA可视化如下: 二、SCF计算 1.直接套用第(1)节的代码,把KIND部分的元素替换成H和O 完整的代码见文末,因为CP2K的输入较为冗长,所以常见的构造.inp文件的方式是套用之前的模板。(备注:后来发现这里的基组最...
• 如果之前没用outer SCF则开启之,这对解决SCF不收敛非常有效,也即加入&SCF / &OUTER_SCF字段。由于其重要性,Multiwfn产生的输入文件里开OT后默认就是启用outer SCF的。此时每次到达&SCF / MAX_SCF所设的inner SCF的步数上限时,就会进入outer SCF过程更新OT的preconditioner,然后进入下一次的inner SCF迭代过程。
另外,CP2K 的 SCF 循环可以用 OUTER_SCF 来帮助收敛,所以 MAX_SCF 不用设置为 128 这么大,一般 ...
qs_outer_scf USE cp_control_types, ONLY ddapc_restraint_type,& dft_controltype,& s2_restraint_type USE cp_log_handling, ONLY: cp_to_string USE input_constants, ONLY: & broyden_type_1, broyden_type_1_explicit, broyden_type_1_explicit_ls, broyden_type_1_ls, & broyden_type_2, ...
每次OUTER_SCF中优化的次数依然是100次,最多可以进行5次OUTER_SCF。所以,最多可以进行500次SCF计算。更多的细节清参考:https:〃/?fromgrouDs=#!topic/cp2k/要输出每个原子上的受力,GLOBAL部分的RUN_TYPE必须设置为ENERGY_FORCE或者GEOJDPT。要输出受力,需要在FORCE_EVAL部分开启选项:&PRINT&FORCESON...
一步完成所有插件编译,没有任何报错。 plumed=installMPIis detected and it appears to be OpenMPI Compilingwith56processes.===FindingGCCfromsystem paths===path to gcc is/usr/bin/gcc path to g++is/usr/bin/g++path to gfortran is/usr/bin/gfortran Found include directory/usr/include Step gcc ...
在CP2K中使用B3LYP泛函时,关键输入文件如下: DFT BASIS_SET_FILE_NAME ./BAS1S_MOLOPT POTENTIAL_FILE_NAME ./POTENTIAL CHARGE MULTIPLICITY 0 1 SCF SCF_GUESS EPS_SCF MAX_SCF ATOMIC 1.0E-6 50 OUTER_SCF MAX_SCF 10 END OUTER_SCF 0T # My scheme PRECONDITIONER FULL.SINGLEJNVERSE MINIMIZER N_D...