&FORCE_EVAL小节中包括了体系的相关参数,如原子的晶胞、坐标、元素、基组(BASIS_SET)和赝势(POTENTIAL)、磁性、K点、截断能、泛函、收敛精度及收敛方法等信息; 在SCF部分,我们可以对自洽循环圈数以及收敛标准进行修改。 MAX_SCF 128 EPS_SCF 1.0E-06 #Convergence threshold of density matrix of inner SCF 需要...
MAX_STEPS 200 &END DIIS &END OPTIMIZE_BAND &CONVERGENCE_CONTROL !收敛判据,过渡态可以比结构优化稍微宽松一点 MAX_FORCE 0.003 RMS_FORCE 0.005 MAX_DR 0.002 RMS_DR 0.005 &END &REPLICA !分别输入初末态的结构,CP2K会自动在这些结构和SUBSYS里的结构之间插点 COORD_FILE_NAME ini_opted.xyz &END &REPL...
结构优化的inp文件和scf的文件如下。结构优化精度是MAX_FORCE 2*10^-9,EPS_SCF是1E-7。k点是单...
计算能量,需要设置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_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运行的一般选项,例如作业的名称、运行的类型等。
MAX_SCF 128 EPS_SCF 1.0E-06 #Convergence threshold of density matrix of inner SCF # SCF...
MAX_ITER200MAX_FORCE0.0009725OPTIMIZER BFGS&END GEO_OPT 很明显,博主进行的是几何优化计算,是静态的。当时回看手册,EFIELD要进行的是含时的,非静态的,非周期性的计算。 CP2K官网对于EFIELD模块的注释 所以施加匀强电场时,尤其对于静态计算,是不能使用EFIELD的,这就是第一个误区。
几何优化可采用最小化局部能量(MINIMIZATION)或寻找过渡态的鞍点(TRANSITION_STATE),通过设置MAX_DR,MAX_FORCE,RMS_DR和RMS_FORCE阈值来确定终止条件。在CG参数中,禁止在共轭梯度下降之前执行最速梯度下降,并在角度的余弦值低于0.9时重置共轭梯度。通过MOTION部分的CONSTRAINT设定,可以对特定原子进行...
计算能量,需要设置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_scf=None, cutoff=None, xc=None, force_eval_method=None, basis_set=None, pseudo_potential=None, basis_set_file=None, potential_file=None, stress_tensor=False, poisson_solver=None, print_level=None, command="./cp2k_shell.ssmp --shell", ) # %% # We now load a new struture, add...