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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】vasp 频率计算虚频
51/1返回列表
查看: 4209  |  回复: 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的回帖

dummy1

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

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

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

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的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +9 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-09 10:50 by Equinoxhua
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-09 10:18 by Equinoxhua
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 13/650 2026-11-09 06:50 by l48hFd8Nwae2
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 06:39 by l48hFd8Nwae2
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5 5sJFDKp2kwuI 2026-11-06 9/450 2026-11-08 22:46 by l48hFd8Nwae2
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 09:01 by 2lRNluKrJpwL
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 kvewwGtKktqB 2026-11-04 11/550 2026-11-06 05:10 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 04:47 by qL1cU46YIdK0
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-06 02:15 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 01:42 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
关闭
请填处理意见
关闭