计算化学公社

标题: 使用CP2K进行GFN1-xTB计算不收敛问题 [打印本页]

作者
Author: northern_mar 时间: 2026-1-12 19:43
标题: 使用CP2K进行GFN1-xTB计算不收敛问题
各位老师同学好，近期我在使用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的无法收敛还是计算参数设置问题导致的？问题出在哪呢？
谢谢！

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

作者
Author: northern_mar 时间: 2026-1-12 20:57

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

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

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

作者
Author: northern_mar 时间: 2026-1-13 09:17

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

好的，谢谢卢老师！

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