计算化学公社

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

[VASP] 求助:VASP的CI-NEB计算均存在小虚频或较大虚频,如何消除?

[复制链接 Copy URL]

44

帖子

0

威望

207

eV
积分
251

Level 3 能力者

本帖最后由 红烧肉盖饭 于 2023-12-19 10:10 编辑

VASP消虚频的方法,我在论坛中搜索过了,像加上“打破对称性,ISYM=-1”“使用vaspkit的507程序”等等,均未成功。下面是我的一些结果数据。
1.我的CI-NEB计算INCAR文件以及频率计算文件。频率计算的POSCAR均是对应结构的CONTCAR。
INCAR-freq (152 Bytes, 下载次数 Times of downloads: 11)
INCAR-NEB (537 Bytes, 下载次数 Times of downloads: 4)

2.这里贴的是我一个初态结构的频率计算结果。

这是我在这个文件夹下,使用vaspkit的507程序得到的输出。可以发现,有“Frequencies have excluded the rotational and translational ones! No imaginary-freq to be corrected!”的输出结果。

我的问题如下:
1.我列出的这个结构计算的虚频,作为初态有点过于大了。有的结构算出的初态,只有不到10i的虚频。这个虚频是可以忽略的还是必须消除掉?
2.如果是必须消除的,有什么具体的解决方法呢?我有试过jmol查看OUTCR文件,但是不知道如何保存它在虚频上的扭转结构。请问能给出具体的解决方法吗?
谢谢各位老师。





96

帖子

0

威望

894

eV
积分
990

Level 4 (黑子)

2#
发表于 Post on 2023-12-18 09:54:29 | 只看该作者 Only view this author
1,没看到过具体文献说可以忽略小的虚频,所以我个人认为还是不可忽视的,即便确实很小。
这部分请见
http://bbs.keinsci.com/thread-13344-3-1.html
2,可以试试用主要虚频作为MODECAR来跑个dimer,从经验上来讲这样跑个dimer小的虚频都可以消除掉。
dimer的使用方法请看
http://bbs.keinsci.com/forum.php ... amp;highlight=dimer

44

帖子

0

威望

207

eV
积分
251

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-18 10:16:17 | 只看该作者 Only view this author
白皛白 发表于 2023-12-18 09:54
1,没看到过具体文献说可以忽略小的虚频,所以我个人认为还是不可忽视的,即便确实很小。
这部分请见
htt ...

非常感谢,我仔细看看博文

2

帖子

0

威望

97

eV
积分
99

Level 2 能力者

4#
发表于 Post on 2023-12-28 17:10:46 | 只看该作者 Only view this author
白皛白 发表于 2023-12-18 09:54
1,没看到过具体文献说可以忽略小的虚频,所以我个人认为还是不可忽视的,即便确实很小。
这部分请见
htt ...

您好,请问如何将主要虚频提取出来作为MODECAR来跑dimer

96

帖子

0

威望

894

eV
积分
990

Level 4 (黑子)

5#
发表于 Post on 2023-12-29 11:22:49 | 只看该作者 Only view this author
bingo.0 发表于 2023-12-28 17:10
您好,请问如何将主要虚频提取出来作为MODECAR来跑dimer

1,在计算频率的文件夹下,grep cm-1 OUTCAR就可以看到有哪些频率,像图1这样。
2,以图1为例子,grep “12 f/i” OUTCAR -A n+2  > freq    实际选择哪个虚频看你的实际情况,n是你的体系的原子数,这一步做完请检查下是否包含了全部原子。
3,把freq文件中最右侧三列的全部数据如图2红色框起来的那些,复制到一个新的文本文件里,将该文本文件命名为MODECAR即可。

202312291119312026..png (48.78 KB, 下载次数 Times of downloads: 25)

图2

图2

202312291114068918..png (28.76 KB, 下载次数 Times of downloads: 17)

图1

图1

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
123wo + 3 赞!

查看全部评分 View all ratings

86

帖子

0

威望

1777

eV
积分
1863

Level 5 (御坂)

碳基生物,地表动物

6#
发表于 Post on 2023-12-29 15:27:35 | 只看该作者 Only view this author
qvasp有消虚频的功能,另外,你可以拿这个结构做个IDM或者Dimer试试

2

帖子

0

威望

97

eV
积分
99

Level 2 能力者

7#
发表于 Post on 2024-1-2 20:01:12 | 只看该作者 Only view this author
白皛白 发表于 2023-12-29 11:22
1,在计算频率的文件夹下,grep cm-1 OUTCAR就可以看到有哪些频率,像图1这样。
2,以图1为例子,grep  ...

谢谢老师,我去试一下

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 11:52 , Processed in 0.535872 second(s), 24 queries , Gzip On.

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