COF簇模型SCF不收敛

wangs · 发表于 Post on 2023-7-5 21:06:00

本人想要定性计算锂离子通过COF的迁移能垒，先是用VASP优化了COF的结构，其次从中取出一个簇模型，边缘用H原子饱和，在COF的孔道中心插入一个锂离子，想要高斯中进行单点能计算，输入文件如下，想问下各位老师，有没有什么好的解决方法吗？

sobereva · 发表于 Post on 2023-7-6 10:52:58

wangs 发表于 2023-7-6 10:11
多谢老师指点，我是想比较锂离子从COF中心传输和靠近磺酸根传输时的能垒不同，从而说明磺酸根促进锂离子 ...

这么算不靠谱，很容易遭质疑，相邻层的COF对计算结果产生的影响根本不可忽略。就算你用Gaussian来算至少也得把相邻层COF纳入进去，但这时候耗时已经挺高了

wangs · 发表于 Post on 2023-7-6 10:11:30

sobereva 发表于 2023-7-6 10:06
算这个体系的意义完全不明。你如果要算离子通过COF的迁移能垒，应当用CP2K等第一性原理程序先优化Li+在孔洞 ...

多谢老师指点，我是想比较锂离子从COF中心传输和靠近磺酸根传输时的能垒不同，从而说明磺酸根促进锂离子传输，以上方法是我们导师让我这么算的

sobereva · 发表于 Post on 2023-7-6 10:06:29

算这个体系的意义完全不明。你如果要算离子通过COF的迁移能垒，应当用CP2K等第一性原理程序先优化Li+在孔洞内的结构，然后跑诸如CI-NEB，用过渡态和极小点能量求差。拿Gaussian孤零零算这么个东西什么用也没有

wangs · 发表于 Post on 2023-7-6 09:33:06

Eudaimonia 发表于 2023-7-6 08:52
检查体系设置的电荷对不对

如果没问题，我猜测SCF不收敛的主要原因可能是初始构型离谱，Li离子不太可能 ...

好的，多谢老师指点，我再改改

Eudaimonia · 发表于 Post on 2023-7-6 08:52:48

检查体系设置的电荷对不对

如果没问题，我猜测SCF不收敛的主要原因可能是初始构型离谱，Li离子不太可能孤零零的在孔道正中心，要么放到旁边的SO3附近，要么加溶剂分子把它裹起来

wangs · 发表于 Post on 2023-7-5 21:52:50

niobium 发表于 2023-7-5 21:28
你中间那一圈SO3都没有用氢饱和，所以你的电荷设定得对不对？
另外，这么大的孔，感觉锂离子不会孤零零地 ...

我特意把中间的H去掉的，想看看负电性的SO3对于锂离子的影响，第二个是因为想比较在孔道正中心和吸附在磺酸根上时的锂离子扩散能垒不同，所以要分别设置其在中心和吸附在磺酸根，本人的COF是锂金属电池中的隔膜

niobium · 发表于 Post on 2023-7-5 21:28:17

本帖最后由 niobium 于 2023-7-5 21:30 编辑

你中间那一圈SO3都没有用氢饱和，所以你的电荷设定得对不对？
另外，这么大的孔，感觉锂离子不会孤零零地悬在中间吧，应该会吸附到某个SO3上吧

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

[Gaussian/gview] COF簇模型SCF不收敛