我对当前这个近收敛结构又单独做了一次 ENERGY_FORCE 检查。输出里 SCF 正常收敛，warnings 为 0，也正常结束：
SCF run converged in 5 steps
ENERGY| Total FORCE_EVAL ( QS ) energy [hartree] -20401.861195213728934
The number of warnings for this run is : 0
PROGRAM ENDED
grep 没有发现 SCF run NOT converged、ERROR、ABORT、FATAL、OOM、OUT OF MEMORY、TIME LIMIT 等异常。
这个结构的力大致是：
E = -20401.861195213728934 Ha
Fmax_norm = 2.1425056445E-03 Ha / Bohr
Frms_norm = 2.4250835498E-04 Ha / Bohr
Top force atom = O635
可以看到 RMS step size 和 RMS gradient 已经满足，但是 Maximum step size 和 Maximum gradient 还没有满足：
Maximum step size = 0.0094486307, NO
RMS step size = 0.0007242304, YES
Maximum gradient = 0.0019521460, NO
RMS gradient = 0.0001388723, YES
后面我又从这个结构继续做过小步 BFGS 和 trust radius 测试，其中有一个中间 block 是三项满足，只有 Maximum gradient 不满足，
Maximum step size = 0.0028345892, YES
RMS step size = 0.0001863851, YES
Maximum gradient = 0.0018455107, NO
RMS gradient = 0.0001306425, YES
不过这个继续优化测试没有很稳定地把最大力降下去，所以我没有直接把对应的测试结构当成最终基体。
我还做了一个考虑周期性边界的短接触检查，没有发现明显短接触：
O-H too short cutoff 0.600 count 0
H-H too short cutoff 0.800 count 0
O-O too short cutoff 1.500 count 0
Si-O too short cutoff 1.200 count 0
Al-O too short cutoff 1.200 count 0
Ca-O too short cutoff 1.800 count 0
最大力对应的 O635 附近距离如下：
O635-H1282 = 0.977884 A
O635-H1277 = 1.608170 A
O635-Al364 = 1.846455 A
O635-O651 = 2.618448 A
Nearest Ca: Ca135 = 3.018074 A, Ca69 = 3.153868 A从距离上看，我感觉这个最大力可能和局部 Al-OH 以及氢键环境有关，没有明显的短接触，但最大力确实还没有消掉。
想请教一下：
对这种约 1500 原子的周期性 C-A-S-H 体系，如果 RMS force 已经比较低，但是 Fmax 还在 2.1E-3 Ha / Bohr 左右，而且最大力集中在一个局部 Al-OH 和氢键区域，这种结构是否仍然不建议作为后续缺陷和吸附模型的基体？
如果后续 Cu 模型尽量避开这个最大力所在区域，是否可以把这个结构作为一个受限的基体结构使用，并在文章里明确说明它不是 fully relaxed structure？还是说只要基体没有严格收敛，就不建议用来做后续相对能量比较？
如果后面要比较 Cu 替代、吸附、层间配位或者局部络合等构型，是否必须先把这个基体继续优化到更低 Fmax？或者可以通过做几个 parent sensitivity check 来判断这个残余力对能量排序的影响？