计算化学公社

标题: xtb中的GFN-xTB 1代在算原子团簇时依然有意义 [打印本页]

作者
Author: sobereva 时间: 2019-3-19 22:54
标题: xtb中的GFN-xTB 1代在算原子团簇时依然有意义
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文件是 (, 下载次数 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个）： (, 下载次数 Times of downloads: 18) 。而且得到了很理想的团簇结构，其中拿两个有代表性的：

(, 下载次数 Times of downloads: 68)

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

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

作者
Author: fhh2626 时间: 2019-3-19 23:32
请问具体的原因可能是什么呢？

作者
Author: zsu007 时间: 2019-3-20 08:39
谢谢社长的分享！

作者
Author: houhou 时间: 2019-3-20 19:48
我之前算Fe或者Fe的化合物团簇也基本不收敛，即使收敛结构也很怪，回去用GFN 1试试

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