计算化学公社

标题: 截断能测试请教 [打印本页]

作者
Author: 冰河洗剑 时间: 3 day ago
标题: 截断能测试请教
Cutoff |  Energy (Hartree)  | delte E | NG on grid 1 | NG on grid 2 | NG on grid 3 | NG on grid 4
   50.0    -472.2819672557                         0             0             0             0
   100.0    -470.6751934136 1.6067738421          0             0             0             0
   150.0    -470.4968160865 0.1783773271          0             0             0             0
   200.0    -470.4782294308 0.0185866557          0             0             0             0
   250.0    -470.4747374882 0.0034919426          0             0             0             0
   300.0    -470.4733675802 0.0013699079          0             0             0             0
   350.0    -470.4698608261 0.0035067541          0             0             0             0
   400.0    -470.4696340860 0.0002267401          0             0             0             0
   450.0    -470.4689384287 0.0006956573          0             0             0             0
   500.0    -470.4693119620  -0.0003735333          0             0             0             0
   550.0    -470.4689087165 0.0004032455          0             0             0             0
   600.0    -470.4689612296  -0.0000525131          0             0             0             0
   650.0    -470.4690423747  -0.0000811451          0             0             0             0
   700.0    -470.4691843913  -0.0001420166          0             0             0             0
   750.0    -470.4691015262 0.0000828651          0             0             0             0
   800.0    -470.4690758475 0.0000256788          0             0             0             0

在尝试使用卢老师CP2K课程中的脚本做截断能测试，存在如下问题：
① NG on grid 无论是在Medium还是High下都读不到。
② ΔE始终无法很好收敛，达不到卢老师建议的1E-6收敛限。
烦请指正！感恩！

作者
Author: UW_0728. 时间: 3 day ago
本帖最后由 UW_0728. 于 2026-4-10 23:24 编辑

这么多原子数，建议按照1E-6 Hartree/atom的标准来判断

# IGNORE_CONVERGENCE_FAILURE #Continue calculation even if SCF not converged, works for version >= 2024.1

复制代码

这行前面的注释应当去掉

PRINT_LEVEL设置为HIGH没啥意义，用MEDIUM

REL_CUTOFF也太大了吧，60管够了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)