24小时热门版块排行榜    

查看: 3980  |  回复: 13
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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 ]
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 14 个回答

huangyc

木虫 (正式写手)

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

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zxzj05(金币+1,VIP+0):^_^ 3Q!鼓励答疑! 6-9 22:08
首先可能是你的自恰有问题。
其次,你的虚频的值大概是多少,如果小的话,没什么问题。你计算的是G点,虚点正常。
真理是一点点接近的!
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试试
4楼2009-06-08 23:06:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于Filecode分析方法 +8 majunge000 2026-08-10 9/450 2026-08-12 22:17 by 虫友是什么虫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见