求助：在castep中使用相同参数对同一吸附模型(ZnO-NO2)进行几何优化时结果不同

巴里巴里巴里

各位老师好，我先对ZnO单胞进行几何优化，再进行超胞P（3x1），切（100）面，thickness为3，固定两层原子，表面四层原子弛豫，截断能为340eV，k点2x4x1，对这个模型进行几何优化和计算能量，之后将NO2分子放在相同格子里进行几何优化和计算能量，再将优化后的NO2分子放在ZnO表面，对吸附模型进行几何优化，图1是原始模型。
图2是使用的参数。
第一次对该模型进行几何优化，程序运行11个小时，未能结束，因为还要使用电脑，我主动停止，收敛趋势图震荡严重。 第二次对相同模型使用相同参数进行几何优化，运行约8.8小时得到结果，效率为100%。收敛趋势图为
两次几何优化时.castep文件中均有三个warning存在:
1.BFGS: WARNING - trial step suggests complex energy landscape
                 in which simple line minimization will fail.

               - This is usually an indication that the forces/streses
                 are not accurate enough. Consider making the following changes:

                 1) Increase cutoff_energy or basis_precision
                 2) Reduce elec_energy_tol and/or elec_force_tol
                 3) Increase grid_scale
                 4) Increase fine_grid_scale

               - Proceeding with a bisection search to find root instead.
2.BFGS: WARNING - looks like bisection search has failed?
BFGS:         - reverting to earlier configuration and reseting inverse Hessian
3. BFGS: WARNING - Repeated consecutive reset of inverse Hessian
BFGS:           without satisfying convergence criteria which
BFGS:           looks like BFGS has run out of search directions.
BFGS: WARNING - Lets try allowing some uphill steps and see if
BFGS:           we can get around this barrier.
BFGS:         - It is possible that the system may now converge to
BFGS:           a stationary point OTHER than the desired minimum.
BFGS: Hint    - this may be an indication that either:
BFGS:           a) you are using a poor guess at geom_frequency_est
BFGS:           and/or geom_modulus_est, or
BFGS:           b) you are using unrealistic convergence criteria.
BFGS:           Suggest therefore that you consider changing them!

之前也出现过这样的情况，白天开始优化，优化到第二天依旧优化不出来，晚上优化就可以优化出来，太奇怪了，每次用的模型和参数都一样但是结果不同，求助各位老师出现这种现象的原因。

WilliamH

可以增大截断能试试。我之前算过类似的体系，发现如果截断能太低的话，无论是收敛速度还是优化后的结果都很差。   建议你把色散校正打开。   此外，我看你设置了初始自选，自旋极化是忘记打开了吗？

巴里巴里巴里

WilliamH 发表于 2024-11-13 18:55
可以增大截断能试试。我之前算过类似的体系，发现如果截断能太低的话，无论是收敛速度还是优化后的结果都很 ...

谢谢老师的回答，我今天尝试一下把截断能设置到380eV，几何优化时打开自旋极化的话会非常耗时，然后发现不开也能算就没开了，算能量的时候开了色散校正跟自旋极化。

zhangpu@

您好，我再计算弹性常数的时候也出现了1和3报错，请问楼主您这个问题解决了么？

巴里巴里巴里

zhangpu@ 发表于 2025-8-28 19:56
您好，我再计算弹性常数的时候也出现了1和3报错，请问楼主您这个问题解决了么？

解决了，电脑配置不足，用了超算平台就没报错了