计算化学公社

标题: 优化结构的频率分析不收敛问题 [打印本页]

作者
Author: Hillary 时间: 2023-7-3 13:56
标题: 优化结构的频率分析不收敛问题
本帖最后由 Hillary 于 2023-7-3 16:11 编辑

计化新新新人一枚，在对已有的沸石结构进行过渡金属掺杂并完成结构优化后，再进行频率分析的时候，发现离子步已经跑了1200多步，还是没有完成计算，不知道是哪个参数出了问题。希望各位老师和有过同样问题的同学可以能够帮助解答，非常感谢，以下是我的INCAR参数。
GGA = PE
PREC = Accurate
EDIFF  = 1E-7
ENCUT  = 400.0
NELMIN = 5
NSW =  500
NELM = 200
EDIFFG =  -0.015
ISIF = 0

ISMEAR = 0

IBRION = 5
NFREE  = 2
ISYM = 0
LVDW = .TRUE.

LCHARG = .TRUE.
LWAVE  = .TRUE.
ISTART = 1

INSPIN = 2
MAGMOM = 144*0 71*0 1*8 1*0

LDAU = .T.
LDAUL  = -1 -1 2 -1
LDAUU  = 0.0 0.0 2.5 0.0
LDAUJ  = 0.0 0.0 0.0 0.0

作者
Author: 含光君 时间: 2023-7-3 19:40
频率分析又不是结构优化，1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全原子频率分析。

作者
Author: Hillary 时间: 2023-7-4 15:35
本帖最后由 Hillary 于 2023-7-4 15:41 编辑

含光君发表于 2023-7-3 19:40
频率分析又不是结构优化，1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全 ...

超级感谢老师的回复！！我只以为任务有问题就停了··· 然后把EDIFFG改成了-0.02 就又重新开始跑了，刚刚去看了我的POSCAR，发现原子都是T固定的，所以老师所说的固定原子是什么意思嘞？

作者
Author: 北大-陶豫 时间: 2023-7-4 16:23

Hillary 发表于 2023-7-4 15:35
超级感谢老师的回复！！我只以为任务有问题就停了··· 然后把EDIFFG改成了-0.02 就又重新开始跑了，刚 ...

F才是固定，你就记 F=Frozen

作者
Author: Hillary 时间: 2023-7-4 16:53

北大-陶豫发表于 2023-7-4 16:23
F才是固定，你就记 F=Frozen

感谢老师的指正！我去试一试！

作者
Author: 含光君 时间: 2023-7-4 16:57
emmm感觉你还不是很清楚VASP中频率分析的流程和作用，建议先补习相应知识。另外，T是不固定F是固定。

作者
Author: Hillary 时间: 2023-7-4 17:09

含光君发表于 2023-7-3 19:40
频率分析又不是结构优化，1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全 ...

还想请问老师的一点是，固定的原子是固定我掺入过渡金属和附近的几个原子，还是？

作者
Author: 含光君 时间: 2023-7-4 17:57

Hillary 发表于 2023-7-4 17:09
还想请问老师的一点是，固定的原子是固定我掺入过渡金属和附近的几个原子，还是？

不知道你对这个结构具体想研究的性质是什么，就不太好准确回答。

作者
Author: Hillary 时间: 2023-7-5 13:11

含光君发表于 2023-7-4 16:57
emmm感觉你还不是很清楚VASP中频率分析的流程和作用，建议先补习相应知识。另外，T是不固定F是固定。

是的，我是属于那种自己摸进门的，有好多东西确实一知半解

作者
Author: Hillary 时间: 2023-7-5 13:24

含光君发表于 2023-7-4 17:57
不知道你对这个结构具体想研究的性质是什么，就不太好准确回答。

其实我是接手师兄师姐的工作，他们原来研究的这个沸石结构大体是在研究烯烃在该沸石上的过渡态，然后老师想把Al换成过渡金属，让我来做这个尝试，但是并没有具体说想研究什么性质，请问老师有没有什么建议，然后我看到POSCAR里的原子都是T的状态，那么频率分析算不完还是原子固定的问题吗？

作者
Author: 含光君 时间: 2023-7-5 13:27
本帖最后由含光君于 2023-7-5 14:02 编辑

这种情况沸石骨架用不着算频率吧。一般算吸附结构固定骨架，仅对吸附分子算频率。你们自己组的工作当然直接请教师兄师姐最为方便可靠。

作者
Author: Hillary 时间: 2023-7-5 13:44

含光君发表于 2023-7-5 13:27
这种情况沸石骨架用不着算频率吧。一搬算吸附结构固定骨架，仅对吸附分子算频率。你们自己组的工作当然直接 ...

是这样没错，确实也是师兄师姐不是很明白为什么，似乎他们一贯优化完最初的沸石结构会频率分析一下，看是否有虚频？我想问一下老师我换完原子的结构优化完以后，除了看是否有没有虚频，是否还有别的的办法看一下我这个结构是否合理（优化完的结构已从VESTA里看过了，没有啥大问题），还是直接既可以往上吸附分子了？

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