• 尝试将&QS / EPS_DEFAULT减小几个数量级,默认为1E-10。EPS_DEFAULT这个参数控制一堆涉及数值计算精度的参数,如果设得太大,有可能造成SCF难收敛。若这个数值当前已经很小了,比如都1E-12了,尝试继续减小则没必要。 • 改用其它泛函、基组。如果原先的计算级别就是极难收敛,而有其它级别性价比、精度相仿佛...
输入文件的精度如下：EPS_SCF 1.0E-4 ，EPS_SCHWARZ 1.0E-4，EPS_DEFAULT 1.0E-5，不知是否...
我不爱用cp2k-1.restart是因为里面格式变得和我设置的不一样了，比如EPS_SCF 1E-05到cp2k-1....
123&DFT45BASIS_SET_FILE_NAME ./BASIS_MOLOPT67POTENTIAL_FILE_NAME ./POTENTIAL89CHARGE01011MULTIPLICITY11213&SCF1415SCF_GUESS ATOMIC1617EPS_SCF1.0E-61819MAX_SCF502021&OUTER_SCF2223MAX_SCF102425&END OUTER_SCF2627&OT2829# My scheme3031PRECONDITIONER FULL_SINGLE_INVERSE3233MINIMIZER DIIS3435N_DIIS73637&...
EPS_SCF 1.0E-06 #Convergence threshold of density matrix of inner SCF 需要注意的是,在&FORCE_EVAL小节中的SCF收敛中,我们对于Cu体系采用的算法是基于对角化的算法(DIAG)。另外一种算法是基于轨道变换的算法(OT)。通常,对于金属体系,我们使用对角化算法,而对于具有较大带隙的体系,如半导体或绝缘体,则推荐使用...
阈值通过&DFT/&XC/&HF/&SCREENING中EPS_SCHWARZ设置,此值越大,被忽略的ERI数目越多、计算精度越差。默认的1E-10虽然足够保证精度,但此时太贵,没必要。通常设为1E-6就可以了,是精度和耗时的较好权衡,想要精度更好点用1E-8。 (2)结合密度矩阵做积分屏蔽:ERI对能量的贡献取决于它与特定密度矩阵元的乘积。
123&DFT45BASIS_SET_FILE_NAME ./BASIS_MOLOPT67POTENTIAL_FILE_NAME ./POTENTIAL89CHARGE01011MULTIPLICITY11213&SCF1415SCF_GUESS ATOMIC1617EPS_SCF1.0E-61819MAX_SCF502021&OUTER_SCF2223MAX_SCF102425&END OUTER_SCF2627&OT2829# My scheme3031PRECONDITIONER FULL_SINGLE_INVERSE3233MINIMIZER DIIS3435N_DIIS73637...
EPS_DEFAULT 1.0E-10&END QS&MGRID NGRIDS 4 CUTOFF 300 REL_CUTOFF 60&END MGRID&XC&XC_FUNCTIONAL PADE&END XC_FUNCTIONAL&END XC&SCF SCF_GUESS ATOMIC EPS_SCF 1.0E-7 MAX_SCF 300&DIAGONALIZATION ON ALGORITHM STANDARD&END DIAGONALIZATION&MIXING T ...
CP2K计算频率使用的是数值算法,即对每个原子向+x, -x, +y, -y, +z, -z 6个方向分别进行移动,用数值的方法得到能量的二阶导(即力常数),然后计算频率。所以,如果有N个原子要进行移动,总共要进行6N+1次SCF收敛计算。 要计算部分原子的振动频率,有两种办法。一种是直接在MOTION中使用CONSTRAINT对不需要进行...
EPS_SCF 1.0E-05 &END OUTER_SCF &END SCF &XC &XC_FUNCTIONAL PBE &END XC_FUNCTIONAL &END XC &END &SUBSYS &CELL ABC 19.133 19.949 19.98 PERIODIC NONE &END CELL &COORD 坐标 &END COORD &KIND C #研究体系的原子种类,有多少个就写多少个,每个都要对应地写基组和赝势 ELEMENT C BASIS_SET DZVP...