求助：windows下使用xtb6.5.1与6.6.0优化结构，均出现报错信息：forrtl: severe (172)

sophia527

windows下使用xtb6.5.1与6.6.0优化结构，均出现报错信息：forrtl: severe (172): Program Exception - exception code = 0x7e (126)，请问是为什么？如何调整？麻烦了！
计算单分子或者二聚体都是一样的报错。

# atoms in fragment 1/2:   216   216
fragment masses (1/2) :     1475.84     1475.84
CMA distance (Bohr)    :  27.252
constraining FC (au)   :  0.0500
forrtl: severe (172): Program Exception - exception code = 0x7e (126)
Image              PC                Routine            Line        Source
KERNELBASE.dll     00007FF9FB98CF19  Unknown               Unknown  Unknown
mkl_intel_thread.  00007FF95CB10A26  Unknown               Unknown  Unknown
mkl_intel_thread.  00007FF95CB103B5  Unknown               Unknown  Unknown
xtb.exe            00007FF677C3D756  Unknown               Unknown  Unknown
xtb.exe            00007FF67728B976  XTB_EEQ_mp_GOEDEC         482  eeq_model.f90
xtb.exe            00007FF676EED9CB  XTB_PROG_MAIN_mp_         557  main.F90
xtb.exe            00007FF676EFF51F  MAIN__                     57  primary.f90
xtb.exe            00007FF677CCCDCE  Unknown               Unknown  Unknown
xtb.exe            00007FF677CD6D00  Unknown               Unknown  Unknown
KERNEL32.DLL       00007FF9FC067614  Unknown               Unknown  Unknown
ntdll.dll          00007FF9FDCA26A1  Unknown               Unknown  Unknown

exity

可以把结构发出来看看。

sobereva

windows版毛病多，有问题时先尝试linux版

sophia527

sobereva 发表于 2023-5-15 19:50
windows版毛病多，有问题时先尝试linux版

老师，我用molclus_1.9.9.9win调用xtb计算，计算单个分子都没问题，很快就出结果了。但是计算由genmer生成的10个二聚体结构时，就会报错，如下：
  *** Configuration     1  ***
Current date: 2023-05-15   Time: 20:29:57
Loading geometry         1 from the inputted trajectory file
Generating xtb.xyz file...
Deleting xtbrestart
Running xtb: xtb xtb.xyz  --opt --gfn 0 --chrg 0 --uhf 0 > xtb.out
Error: Cannot find .xtboptok file!
The task is failed!
Wall clock time elapsed for calculating this configuration:      33 s
但是我对例子中dopamine.xyz的二聚体计算时就没有问题。
这是因为我的分子太大了吗？我一个分子中有220个原子，二聚体就是440个原子。

sophia527

exity 发表于 2023-5-15 17:00
可以把结构发出来看看。

结构不太好发，文章还没发呢

sobereva

sophia527 发表于 2023-5-15 20:32
老师，我用molclus_1.9.9.9win调用xtb计算，计算单个分子都没问题，很快就出结果了。但是计算由genmer生 ...

手动运行molclus提示的命令看xtb.out里显示什么报错相关信息没有

wzkchem5

有可能是bug，建议在xtb的github页面的issue区（https://github.com/grimme-lab/xtb/issues）提一下这个问题

sophia527

sobereva 发表于 2023-5-16 20:41
手动运行molclus提示的命令看xtb.out里显示什么报错相关信息没有

molecular fragmentation (1/2 indicates fragments):
111111111111111111111111111111111111111111111111111111111111111111111111
111111111111111111111111111111111111111111111111111111111111111111111111
111111111111111111111111111111111111111111111111111111111111111111111111
222222222222222222222222222222222222222222222222222222222222222222222222
222222222222222222222222222222222222222222222222222222222222222222222222
222222222222222222222222222222222222222222222222222222222222222222222222
# atoms in fragment 1/2:   216   216
fragment masses (1/2) :     1382.05     1382.05
CMA distance (Bohr)    :  31.521
constraining FC (au)   :  0.0500

           -------------------------------------------------
          |                 G F N 0 - x T B                 |
           -------------------------------------------------

        Reference                      10.26434/chemrxiv.8326202.v1
      * Hamiltonian:
        H0-scaling (s, p, d)           2.000000    2.486800    2.270000
        zeta-weighting                 1.000000
      * Dispersion:
        s8                             2.850000
        a1                             0.800000
        a2                             4.600000
        s9                             0.000000
      * Repulsion:
        kExp                           1.500000
        rExp                           1.000000
      * Coulomb:
        alpha                          0.000000
        third order                    atomic
        anisotropic                    false
      * Polar bond correction:
        rad-shift                      0.053700
        strength                      -0.012900
        en-exp                         3.484700
        en-scale                       0.509700


          ...................................................
          :                      SETUP                      :
          :.................................................:
          :  # basis functions                1288          :
          :  # atomic orbitals                1284          :
          :  # shells                          868          :
          :  # electrons                      1064          :
          :  electronic temp.          300.0000000     K    :
          :  accuracy                    1.0000000          :
          :  -> integral cutoff          0.2500000E+02      :
          :  -> integral neglect         0.1000000E-07      :
          ...................................................
到这就没了。

sophia527

wzkchem5 发表于 2023-5-16 20:57
有可能是bug，建议在xtb的github页面的issue区（https://github.com/grimme-lab/xtb/issues）提一下这个问 ...

好的，谢谢

lyj714

最开始的那个报错我记得是Windows上有的bug，那个6.6.0。最后这个报错，如果你还是用的Windows，那就是omp stack不足，运行时候需要先设置一下，比如cmd先输入set OMP_STACKSIZE=4G

1138711019

sophia527 发表于 2023-5-17 09:10
molecular fragmentation (1/2 indicates fragments):
1111111111111111111111111111111111111111111111 ...

你好，我也出现了这样的情况，面对大体系就直接没了，但是小体系还是可以计算，应该是哪里设置有问题吧？