24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dummy1

木虫 (正式写手)

应助: 13 (小学生)
金币: 5841.7
散金: 32
红花: 1
帖子: 312
在线: 295.8小时
虫号: 270440
注册: 2006-08-06
专业: 理论和计算化学

[交流] 【求助】vasp 频率计算虚频

最近在计算一个单胞是三十多个分子的有机体系，
先采用PBE优化（固定晶格参数），而后做频率分析，结果出现了很多虚频，
单胞时候虚频有六七个，最大是不到100，
超胞选取的是2*1*1(因为体系太大，很慢没有多取，参数由单胞优化得到)，结果虚频更多，甚至上千。
各位有什么解决的办法吗？尤其是超胞是怎么回事阿？
贴出单胞INCAR：
PREC = high
NWRITE =    2; LPETIM=F write-flag & timer
ISTART =    0 job : 0-new  1-cont  2-samecut
ICHARG =    2 charge: 1-file 2-atom 10-const
ISPIN  =    1 spin polarized calculation?
ALGO = normal
RWIGS= 0.37 0.863 1.164
NELM =    200; NELMIN=  8; NELMDL= 10    # of ELM steps
EDIFF  = 0.1E-06 stopping-criterion for ELM
LREAL  = auto

# Ionic relaxation
NFREE  = 2
NSW =    300 number of steps for IOM
IBRION =    5 ionic relax: 0-MD 1-quasi-New 2-CG
POTIM  = 0.01 time-step for ionic-motion
# DOS related values:
ISMEAR =  0 ;  broadening in eV -4-tet -1-fermi 0-gaus
SIGMA=  0.05

单胞Kpoints:

Monkhorst Pack
0
Monkhorst Pack
1 1 1
0  0 0

由于刚刚用vasp，希望各位帮帮忙了，谢谢了！

[ Last edited by dummy1 on 2009-6-8 at 16:09 ]

回复此楼

» 收录本帖的淘帖专辑推荐

有关零点能计算的一些帖子

零点能

» 猜你喜欢

寻求一种能扛住强氧化性腐蚀性的容器密封件已经有7人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼2009-06-08 16:07:11

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dummy1

木虫 (正式写手)

应助: 13 (小学生)
金币: 5841.7
散金: 32
红花: 1
帖子: 312
在线: 295.8小时
虫号: 270440
注册: 2006-08-06
专业: 理论和计算化学

谢谢各位的回复，

我是固定的晶格参数，仅仅做的粒子驰豫。
对于单胞，出现这么多的小的虚频也正常吗？因为我是希望想到准确点的频率，所以担心出现虚频不准确。各位有什么好的办法解决吗？

对于超胞，我是优化的单胞然后用MS扩为2*1*1的超胞的，因为我看到论坛上，还有手册上介绍说VASP
计算频率需要计算超胞，但是不知道是不是在计算超胞频率前也应该先对超胞进行粒子驰豫优化呢？
不驰豫就出现那么大的虚频，实在很困惑！

另外，对于出现几千的虚频为什么要300多原子数以上呢？虚频的大小和原子数有关系吗？我总觉得应该是与元素类型有关系阿，请指教！（另，虚频是用f/i表示，应该没错的）

谢谢各位！

[ Last edited by dummy1 on 2009-6-9 at 09:32 ]