搜MAX_DR可以改几何优化收敛的力与位移判据。一共4项,含义可参照Gaussian的判据理解。这里可以适当放松,例如可以将MAX_DR、RMS_DR、MAX_FORCE、RMS_FORCE分别改为1E-2、5E-3、1E-3、7E-4。 然后就可以运行了。在我的电脑上用我自己编译的cp2k.ssmp用了1941秒算完。作为对比,用EDIFF=5E-5和EDIFFG=-1E-...
&FORCE_EVAL小节中包括了体系的相关参数,如原子的晶胞、坐标、元素、基组(BASIS_SET)和赝势(POTENTIAL)、磁性、K点、截断能、泛函、收敛精度及收敛方法等信息; 在SCF部分,我们可以对自洽循环圈数以及收敛标准进行修改。 MAX_SCF 128 EPS_SCF 1.0E-06 #Convergence threshold of density matrix of inner SCF 需要...
MAX_SCF 128 &END SCF &POISSON &EWALD EWALD_TYPE SPME ALPHA 0.35 GMAX 25 O_SPLINE 6 &END...
结构优化的inp文件和scf的文件如下。结构优化精度是MAX_FORCE 2*10^-9,EPS_SCF是1E-7。k点是单...
MAX_SCF 300&DIAGONALIZATION ON ALGORITHM STANDARD&END DIAGONALIZATION&MIXING T METHOD BROYDEN_MIXING ALPHA 0.4 NBROYDEN 8&END MIXING&END SCF&END DFT&PRINT&FORCES ON&END FORCES&END PRINT&END FORCE_EVAL 输入文件中的主模块包括: ”GLOBAL":包含CP2K运行的一般选项,例如作业的名称、运行的类型等。
计算能量,需要设置RUN_TYPE为ENERGY;如果还要进行梯度(受力)计算,则设置为ENERGY_FORCE。 使用DFT方法进行能量计算,使用OT算法,并开启OUTER_SCF,示例如下: 1&SCF23EPS_SCF1.0E-645SCF_GUESS RESTART67MAX_SCF10089&OT T1011PRECONDITIONER FULL_ALL1213MINIMIZER DIIS1415LINESEARCH 3PNT1617&END OT1819&OUTER_SCF ...
MAX_ITER200MAX_FORCE0.0009725OPTIMIZER BFGS&END GEO_OPT 很明显,博主进行的是几何优化计算,是静态的。当时回看手册,EFIELD要进行的是含时的,非静态的,非周期性的计算。 CP2K官网对于EFIELD模块的注释 所以施加匀强电场时,尤其对于静态计算,是不能使用EFIELD的,这就是第一个误区。
计算能量,需要设置RUN_TYPE为ENERGY;如果还要进行梯度(受力)计算,则设置为ENERGY_FORCE。 使用DFT方法进行能量计算,使用OT算法,并开启OUTER_SCF,示例如下: 1&SCF23EPS_SCF1.0E-645SCF_GUESS RESTART67MAX_SCF10089&OT T1011PRECONDITIONER FULL_ALL1213MINIMIZER DIIS1415LINESEARCH 3PNT1617&END OT1819&OUTER_SCF ...
几何优化可采用最小化局部能量(MINIMIZATION)或寻找过渡态的鞍点(TRANSITION_STATE),通过设置MAX_DR,MAX_FORCE,RMS_DR和RMS_FORCE阈值来确定终止条件。在CG参数中,禁止在共轭梯度下降之前执行最速梯度下降,并在角度的余弦值低于0.9时重置共轭梯度。通过MOTION部分的CONSTRAINT设定,可以对特定原子进行...
Max.stepsize=0.0099827301Conv.limitforstepsize=0.0100000000Convergenceinstepsize=YESRMSstepsize=0.0007891138Conv.limitforRMSstep=0.0050000000ConvergenceinRMSstep=YESMax.gradient=0.0061057290Conv.limitforgradients=0.0100000000Conv.ingradients=YESRMSgradient=0.0002605894Conv.limitforRMSgrad.=0.0050000000Conv.inRMS...