GFN1-xTB不支持含卤素原子的PBC计算么

乘风万里

xTB 6.3.0版本就声称支持了GFN1的PBC计算
然而我在6.4.0版本(conda 安装)下对用这个分子构建的周期性原胞用 --gfn 1 计算就报如下错误

1. ########################################################################
   
2. [ERROR] Program stopped due to fatal error
   
3. -4- Single point calculation terminated
   
4. -3- xtb_calculator_singlepoint: Electronic structure method terminated
   
5. -2- scf: Halogen bond correction not available with PBC
   
6. -1- restart_readRestart: Version number missmatch in restart file.
   
7. ########################################################################

复制代码

但是 --gfn 0 就不会报错，用其它简单分子构建的原胞 --gfn 1 也不报错。
是因为Halogen bond correction的原因么，我看了xtb-docs跟github的releases似乎都没提这个事。
这个分子中都还算是比较常见的主族元素，我不知道是因为这种体系GFN1-xTB还不支持，还是我软件没用对的缘故...

sobereva

建议拿CP2K来跑GFN1-xTB
如果已经有了cif文件，用Multiwfn敲几下键盘就能得到做GFN1-xTB各种任务的输入文件，十分方便
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）
相比CP2K来说，xtb自己的PBC功能太弱了

乘风万里

sobereva 发表于 2021-5-11 05:53
建议拿CP2K来跑GFN1-xTB
如果已经有了cif文件，用Multiwfn敲几下键盘就能得到做GFN1-xTB各种任务的输入文 ...

谢谢社长，今天装了CP2K，的确可以跑完GFN1-xTB下的PBC结构优化，且优化完的结构还算合理。输出文件显示整个过程中的 Correction for halogen bonding 的能量为0 。