xtb中的GFN-xTB 1代在算原子团簇时依然有意义

sobereva · 发表于 Post on 2019-3-19 22:54:48

Grimme的xtb程序如今默认的已经是GFN2-xTB了，虽然还可以用--gfn 1来使用之前的GFN-xTB 1代，但看似已经没有价值了。今天发现对xtb用--gfn 1选项改为1代，对于计算电子结构古怪的情况很有用。

今天用molclus里的genmer组件（molclus和genmer介绍见http://www.keinsci.com/research/molclus.html）通过如下设置构建了一批中空的Li团簇（rmin选项从molclus 1.8里的genmer才开始有）
30 rmin 4 rmax 5.5
examples/Li.xyz

产生的.xyz文件是

traj.xyz (24.14 KB, 下载次数 Times of downloads: 24)

用molclus调用最新版xtb，居然发现一个都收敛不了，最后都提示#WARNING! SCF not converged, aborting...

遂尝试用GFN-xTB 1代，即在molclus的settings.ini的xtb_arg=里面加上--gfn 1选项，molclus运行过程中所有体系优化都顺利收敛了，结果见此文件（只算了traj.xyz前10个）：

isomers.xyz (16.94 KB, 下载次数 Times of downloads: 18) 。而且得到了很理想的团簇结构，其中拿两个有代表性的：

上图左边是genmer产生的初始结构，非常无序，而通过molclus调用xtb做GFN-xTB 1代计算，一个体系几秒钟就优化完，得到了上图右边的结构，非常规整有序。

此例体现xtb算原子团簇效果惊人地好，和molclus相结合获得笼状团簇结构很容易。如前面强调的，原子团簇本来就不容易SCF收敛，初始结构比较诡异时更是难收敛，GFN2-xTB虽然好，但是对这种情况容易碰到不收敛问题，而改用GFN-xTB 1代就解决了。经常用molclus结合xtb做原子团簇结构搜索的人，或者直接用xtb算各种原子团簇的人应注意这一点。

fhh2626 · 发表于 Post on 2019-3-19 23:32:42

请问具体的原因可能是什么呢？

zsu007 · 发表于 Post on 2019-3-20 08:39:00

谢谢社长的分享！

houhou · 发表于 Post on 2019-3-20 19:48:15

我之前算Fe或者Fe的化合物团簇也基本不收敛，即使收敛结构也很怪，回去用GFN 1试试

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

[xtb] xtb中的GFN-xTB 1代在算原子团簇时依然有意义

评分 Rate

浏览过的版块