计算化学公社

标题: xtb这样的问题怎么解决？ [打印本页]

作者
Author: aaq2800 时间: 2023-12-4 14:51
标题: xtb这样的问题怎么解决？
本帖最后由 aaq2800 于 2023-12-6 08:02 编辑

Molclus+xTB进行结构搜索，都是Win10系统，一台电脑上跑，一点问题都没有，但学生的电脑上就不行，gfn1部分成功，部分失败，一半对一半错，（gfn0全部成功，gfn2全部失败），不明白为什么会这样。
错误如下：
*** convergence criteria cannot be satisfied within 250 iterations ***

      # Occupation          Energy/Eh          Energy/eV
   -------------------------------------------------------------
      1       2.0000          -0.6283511          -17.0983
   ...          ...                ...                ...
      12       2.0000          -0.4024432          -10.9510
      13       1.9987          -0.3780246          -10.2866
      14       1.9842          -0.3756165          -10.2210
      15       1.9739          -0.3751298          -10.2078
      16       1.9464          -0.3744237          -10.1886
      17       1.6907          -0.3725407          -10.1373
      18       1.0727          -0.3707807          -10.0895 (HOMO)
      19       0.2575          -0.3684453          -10.0259 (LUMO)
      20       0.0742          -0.3671091             -9.9895
      21       0.0016          -0.3634027             -9.8887
      22                      -0.3419257             -9.3043
      23                      -0.3419248             -9.3042
   ...                               ...                ...
      41                      -0.1148795             -3.1260
   -------------------------------------------------------------
               HL-Gap          0.0023354 Eh          0.0635 eV
         Fermi-level          -0.3706038 Eh       -10.0846 eV

SCC (total)                0 d,  0 h,  0 min,  0.279 sec
SCC setup                   ...       0 min,  0.000 sec (  0.000%)
Dispersion                   ...       0 min,  0.000 sec (  0.000%)
classical contributions       ...       0 min,  0.000 sec (  0.000%)
integral evaluation          ...       0 min,  0.001 sec (  0.358%)
iterations                   ...       0 min,  0.277 sec ( 99.283%)
molecular gradient          ...       0 min,  0.001 sec (  0.358%)
printout                      ...       0 min,  0.000 sec (  0.000%)

last geometry written to: xtblast.xyz

########################################################################
[ERROR] Program stopped due to fatal error
-9- Geometry optimization failed
-8- xtb_calculator_singlepoint: Electronic structure method terminated
-7- scf: Self consistent charge iterator did not converge
-6- xtb_geoopt: Trying to recover from failed geometry optimization
-5- xtb_geoopt: Geometry optimization did not converge
-4- optimizer_ancopt: Could not relax structure
-3- optimizer_relax: SCF not converged, aborting...
-2- xtb_calculator_singlepoint: Electronic structure method terminated
-1- scf: Self consistent charge iterator did not converge
########################################################################

补充一下，搜索xtb的帖子，发现有类似问题还不少。
如果genmer产生随机结构时候，采用共价键半径，没有这问题了，全部成功。采用范德瓦尔斯键，可能是产生的结构不大合适，原子间距比较大，结构有点散，导致很难收敛

作者
Author: GoldenBaby 时间: 2023-12-4 17:50
个人推荐，检查一下你的xtb有没有正确安装。另外，有些体系不适合用GFN2-xtb，也会出现报错，总体上来说，GFN1-xtb通用性更好一点。

作者
Author: sobereva 时间: 2023-12-5 02:44
GFN0-xTB不牵扯SCF迭代，所以肯定不会出现SCF不收敛的报错。GFN2-xTB形式比GFN1-xTB更为复杂，SCF不收敛的可能性往往更高。实际中GFN2-xTB的精度比GFN1-xTB没太显著的优势，GFN1-xTB能跑就用这个就完了。

作者
Author: aaq2800 时间: 2023-12-5 08:39
本帖最后由 aaq2800 于 2023-12-5 08:43 编辑

GoldenBaby 发表于 2023-12-4 17:50
个人推荐，检查一下你的xtb有没有正确安装。另外，有些体系不适合用GFN2-xtb，也会出现报错，总体上来说，G ...

两台电脑都按照同样的方式安装，所以就搞不明白问题出在哪了，有一台linux系统的电脑，使用gfn1也会出现少量的错误，这个还可以忍受

作者
Author: aaq2800 时间: 2023-12-5 10:23
发现一个问题，如果genmer产生随机结构时候，采用共价键半径，没有这问题了，全部成功。采用范德瓦尔斯键，可能是产生的结构不大合适，原子间距比较大，结构有点散，导致很难收敛

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