应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp 频率计算虚频
北京石油化工学院2026年研究生招生接收调剂公告
141/2返回列表12下一页
查看: 3601  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

dummy1

木虫 (正式写手)

[交流] 【求助】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 ]
回复此楼

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

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

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2009-06-08 16:07:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangyc

木虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+1,VIP+0):感谢交流 (*^__^*) …… 6-8 19:31
甚至上千个虚频？你体系多少原子，不得500个以上啊。这么大的体系优化起来很困难，你仅优化晶胞参数再做频率肯定有会问题的，离子也得优化的
赞一下(21人) 回复此楼
2楼2009-06-08 19:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):^_^ 3Q！鼓励答疑！ 6-9 22:08
首先可能是你的自恰有问题。
其次，你的虚频的值大概是多少，如果小的话，没什么问题。你计算的是G点，虚点正常。
赞一下(11人) 回复此楼
真理是一点点接近的！
3楼2009-06-08 21:34:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):^_^ 3Q！鼓励答疑！ 6-9 22:08
楼主，请问VASP里面虚频是 f/i 表示吗？你说虚频上千，至少你的原子上了300个。你full relaxation试试
赞一下(11人) 回复此楼
4楼2009-06-08 23:06:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
谢谢各位的回复，

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

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

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

谢谢各位！

[ Last edited by dummy1 on 2009-6-9 at 09:32 ]
赞一下 回复此楼
5楼2009-06-09 09:30:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xh512

金虫 (正式写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
spur(金币+1,VIP+0): 3Q！！谢谢参与！欢迎常来！ 6-9 12:51
对于出现几千的虚频为什么要300多原子数以上呢？
`````````````
是因为你的文字描述很容易让人误解
可能你想说的是你的虚频值为 1000以上
可是我们看起来，像是 有1000个虚频

1000/3=333 , 当然至少有这么多原子了
赞一下(11人) 回复此楼
6楼2009-06-09 12:44:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
仔细看了下自己发的帖子，的确说的有问题，不好意思了。

我是想表达的意思是：单胞计算频率产生了好几个虚频，采用超胞结果虚频数目更多，有几十个之多，最大虚频都达到了有上千个波数。

表述有误，给各位造成误会了，实在不好意思。见谅！
赞一下 回复此楼
7楼2009-06-09 15:12:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zxzj05(金币+1,VIP+0):^_^ 3Q！鼓励答疑！ 6-9 22:09
引用回帖:
Originally posted by dummy1 at 2009-6-9 09:30:
谢谢各位的回复，

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

对于超胞，我是 ...

频率的数量是弛豫原子的3倍。如果你N个原子弛豫，那么最后VASP给出的gamma点的频率就是3N个。你既然快到1000个虚频了，所以你的原子个数至少有300

[ Last edited by stractor on 2009-6-9 at 15:50 ]
赞一下(11人) 回复此楼
8楼2009-06-09 15:48:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
看来我的语文都成小学水平了，我是说虚频的数目是几十个，最大虚频值为1000波数

大家在计算超胞的时候是怎么做的呢？是要先优化超胞吗？还是也是用的单胞优化结果呢？

单胞我改了potim，prec，ISMEAR,都没有效果呢，大家怎么做的呢？

谢谢
引用回帖:
Originally posted by stractor at 2009-6-9 15:48:

频率的数量是弛豫原子的3倍。如果你N个原子弛豫，那么最后VASP给出的gamma点的频率就是3N个。你既然快到1000个虚频了，所以你的原子个数至少有300

[ Last edited by stractor on 2009-6-9 at 15:50 ]

Sample Text
赞一下 回复此楼
9楼2009-06-09 16:12:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dummy1

木虫 (正式写手)
单胞时候，原子数为40， 虚频几到130个波数都有，有七八个虚频之多。
超胞：采用单胞优化结果，原子数80，虚频几到1000波数，有几十个虚频。

另外，对于超胞，能不能先对超胞优化，再做频率呢？（由于在网上看到说结构优化似乎采用超胞还不如单胞，而且本身体系已经很大了，担心超胞优化困难，所以仍然采用了单胞优化的构型）

谢谢！
引用回帖:
Originally posted by gavinliu7390 at 2009-6-8 21:34:
首先可能是你的自恰有问题。
其次，你的虚频的值大概是多少，如果小的话，没什么问题。你计算的是G点，虚点正常。

赞一下 回复此楼
10楼2009-06-09 16:19:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dummy1 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿211，化学310分，本科重点双非，求调剂 +10 努力奋斗112 2026-04-08 10/500 2026-04-08 09:45 by 上岸上岸，，
[考研] 一志愿深大085601材料工程专业（专硕）300分可以调剂去哪 +9 10160315 2026-04-02 9/450 2026-04-08 09:41 by 小葡萄阿星
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 调剂 +17 不逢春 2026-04-07 18/900 2026-04-08 06:18 by lijunpoly
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 326求调剂 +5 9ahye 2026-04-02 6/300 2026-04-07 21:37 by lijunpoly
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 315求调剂 +3 TUZEIQAQ 2026-04-02 3/150 2026-04-07 17:32 by chenp123
[考研] 277求调剂 +13 倪建设 2026-04-06 13/650 2026-04-07 15:46 by 啊俊！
[考研] 334分机械专硕求调剂 +3 蛋花紫菜汤 2026-04-03 3/150 2026-04-07 14:49 by 逍遥cocoa
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 277求调剂 +5 考研调剂lxh 2026-04-05 5/250 2026-04-05 19:03 by chy09050039
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 一志愿上海大学生物学346 +3 上海大学346调剂 2026-04-03 3/150 2026-04-04 20:20 by dongzh2009
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 085601一志愿北理325分求调剂 +6 找调剂，， 2026-04-02 6/300 2026-04-03 22:20 by 柟�?
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
信息提示
关闭
请填处理意见
关闭