使用CP2K进行GFN1-xTB计算不收敛问题

northern_mar · 发表于 Post on 2026-1-12 19:43:45

各位老师同学好，近期我在使用CP2K计算时想像卢老师课上说的在结构优化前先用GFN1-xTB进行优化。但是使用的时候发现128轮基本无法收敛，也没有收敛趋势。文件使用Multiwfn创建，条件如下：
1. 使用卢老师课上说的，人为设定闭壳层，并且用对角化+smearing来替代开壳层计算。载入输入文件后，-9，将自旋多重度设为1；开启对角化，开启smearing。
2. 使用3,3,1的k点。
3. 体系为H, C, O, Mo, S, Au
但是使用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH使用相同的方式就可以收敛（如图）。请问这是体系问题导致的GFN1-xTB的无法收敛还是计算参数设置问题导致的？问题出在哪呢？
谢谢！

UW_0728. · 发表于 Post on 2026-1-12 20:15:54

若用GFN1-xTB，几乎总是应当使用OT
我建议能用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH正常优化好就不要多此一举

northern_mar · 发表于 Post on 2026-1-12 20:57:32

UW_0728. 发表于 2026-1-12 20:15
若用GFN1-xTB，几乎总是应当使用OT
我建议能用PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH正常优化好就不要多此一举

好的，谢谢您！我是想着节约时间才想试试，因为PBE-D3(BJ)/DZVP-MOLOPT-SR-GTH可能得十来个小时才优化好。

sobereva · 发表于 Post on 2026-1-13 00:02:50

这么小体系用GFN1-xTB预优化毫无必要（如果你就是拿来做正式优化，文章也发不出去，本身就是DFT很容易算的体系），也别高估GFN1-xTB的普适性。就算GFN1-xTB真的收敛了，得到的结构说不定还不如你初始建的。
我课上也明确说了GFN1-xTB等半经验级别的方法应当用OT（而且本身需要靠便宜的GFN1-xTB来算的往往也都是较大晶胞的情况，此时也不需要考虑k点）。

northern_mar · 发表于 Post on 2026-1-13 09:17:32

sobereva 发表于 2026-1-13 00:02
这么小体系用GFN1-xTB预优化毫无必要（如果你就是拿来做正式优化，文章也发不出去，本身就是DFT很容易算的 ...

好的，谢谢卢老师！

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 使用CP2K进行GFN1-xTB计算不收敛问题