求助，在2000，4500，8000的原子体系中，初步的结构优化采用GFN1-xTB还是分子力学

Ying_Zhang

想请教大家一下，遇到大体系，例如我建立的结构，原子数分别是2000，4500和8000。请问在CP2K中，初步弛豫使用半经验的GFN1-xTB还是分子力学，后续计算会使用QMMM计算。

1）目前，采用半经验的GFN1-xTB方法，开了对角，K点1 1 2，但是SCF不收敛，文件如下。

2）想采用分子力学，但是我材料是Mg3Bi2，势函数为BornMayerHuggins (BMH)，从文献提取，不知道是否可靠，文献如下。

3）目前该3个模型为完美晶体，后续还需要插入位错结果，所以综上所述，初步弛豫使用半经验的GFN1-xTB还是分子力学？

麻烦大家提一下意见，感谢。