计算化学公社

标题: LiBH4（010）表面计算CI-NEB收敛标准达不到的问题 [打印本页]

作者
Author: gehan 时间: 2023-6-16 11:06
标题: LiBH4（010）表面计算CI-NEB收敛标准达不到的问题
各位老师好，我想尝试计算一下体系LiBH4（010）表面的过渡态结构，计算一下该结构的脱氢能垒。首先我拿到初始结构并进行结构优化，同时也构造最终想要的结构（让LiBH4(010)表面脱出一个氢气分子，选择的z方向最高的两个氢原子），并且也进行了结构优化，并且这两个结构的结构优化固定了最底下三层的B,Li原子（因为主要是想看看氢气脱出情况，但并不清楚这样合不合理，因为Li离子迁移率也比较高，固定后会不会就不合理了？）。最后结构优化完按照讲义的步骤将restart文件用Multiwfn转为xyz文件，接着构造CI-NEB任务的输入文件，我选择的是晶胞比较大的一个restart文件作为创建neb.inp文件的输入文件。由于服务器我使用的是36核，所以就创建了9个副本同时运行，也就是总共是九个点。但是在NEB-BAND9.out文件中，到63步时，四个收敛标准的值突然增大（出现异常），我本以为可能只是这一步有问题，继续跑。但是后面四个标准就会出现星号（*）,这是不是因为结构有问题，还请各位老师指点一下。图一对应是第62步的收敛标准结果，图二是第63步的结果，图三为78步的结构（有*号出现）。以下附上了对应的初始结构A-opted.xyz，以及预测的最终结构B-opted.xyz，和输入文件。NEB-BAND9.out文件由于文件类型限制无法上传，我将其压缩后上传。

作者
Author: sobereva 时间: 2023-6-16 17:41
尝试提供自己猜的过渡态结构作为插点。并且考虑利用&BAND / &OPTIMIZE_BAND设置步长上限

作者
Author: gehan 时间: 2023-6-17 11:07

sobereva 发表于 2023-6-16 17:41
尝试提供自己猜的过渡态结构作为插点。并且考虑利用&BAND / &OPTIMIZE_BAND设置步长上限

好的，谢谢老师，我更改一下试试

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