计算化学公社

 找回密码 Forget password
 注册 Register
Views: 829|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] 优化结构的频率分析不收敛问题

[复制链接 Copy URL]

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

本帖最后由 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

361

帖子

0

威望

1989

eV
积分
2350

Level 5 (御坂)

所念皆星河

2#
发表于 Post on 2023-7-3 19:40:49 | 只看该作者 Only view this author
频率分析又不是结构优化,1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全原子频率分析。
心之所向,日复一日,必有精进

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-7-4 15:35:42 | 只看该作者 Only view this author
本帖最后由 Hillary 于 2023-7-4 15:41 编辑
含光君 发表于 2023-7-3 19:40
频率分析又不是结构优化,1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全 ...

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

487

帖子

1

威望

4039

eV
积分
4546

Level 6 (一方通行)

4#
发表于 Post on 2023-7-4 16:23:25 | 只看该作者 Only view this author
Hillary 发表于 2023-7-4 15:35
超级感谢老师的回复!!我只以为任务有问题就停了··· 然后把EDIFFG改成了-0.02 就又重新开始跑了,刚 ...

F才是固定,你就记 F=Frozen
现代化学以狄拉克的一句“一切化学问题业已解决”为嚆矢。滥觞于经验主义传统的期望正失去它们的借鉴意义。但面对看似不可达的通往天堂之阶梯,我想循伍德沃德“最好的模型是你底物的对映异构体”的信仰好过过早地振翮。
我们怀揣热忱的灵魂天然被赋予对第一性的追求,不屑于单一坐标的约束,钟情于势能面彼端的芬芳。但

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-7-4 16:53:53 | 只看该作者 Only view this author
北大-陶豫 发表于 2023-7-4 16:23
F才是固定,你就记 F=Frozen

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

361

帖子

0

威望

1989

eV
积分
2350

Level 5 (御坂)

所念皆星河

6#
发表于 Post on 2023-7-4 16:57:35 | 只看该作者 Only view this author
emmm感觉你还不是很清楚VASP中频率分析的流程和作用,建议先补习相应知识。另外,T是不固定F是固定。
心之所向,日复一日,必有精进

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-7-4 17:09:01 | 只看该作者 Only view this author
含光君 发表于 2023-7-3 19:40
频率分析又不是结构优化,1000多步离子步和收不收敛没有任何关系。大概率是忘了固定原子使得程序默认做了全 ...

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

361

帖子

0

威望

1989

eV
积分
2350

Level 5 (御坂)

所念皆星河

8#
发表于 Post on 2023-7-4 17:57:40 | 只看该作者 Only view this author
Hillary 发表于 2023-7-4 17:09
还想请问老师的一点是,固定的原子是固定我掺入过渡金属和附近的几个原子,还是?

不知道你对这个结构具体想研究的性质是什么,就不太好准确回答。
心之所向,日复一日,必有精进

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-7-5 13:11:54 | 只看该作者 Only view this author
含光君 发表于 2023-7-4 16:57
emmm感觉你还不是很清楚VASP中频率分析的流程和作用,建议先补习相应知识。另外,T是不固定F是固定。

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

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2023-7-5 13:24:39 | 只看该作者 Only view this author
含光君 发表于 2023-7-4 17:57
不知道你对这个结构具体想研究的性质是什么,就不太好准确回答。

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

361

帖子

0

威望

1989

eV
积分
2350

Level 5 (御坂)

所念皆星河

11#
发表于 Post on 2023-7-5 13:27:37 | 只看该作者 Only view this author
本帖最后由 含光君 于 2023-7-5 14:02 编辑

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

22

帖子

0

威望

213

eV
积分
235

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2023-7-5 13:44:30 | 只看该作者 Only view this author
含光君 发表于 2023-7-5 13:27
这种情况沸石骨架用不着算频率吧。一搬算吸附结构固定骨架,仅对吸附分子算频率。你们自己组的工作当然直接 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 10:17 , Processed in 0.363929 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list