24小时热门版块排行榜    

查看: 3532  |  回复: 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

木虫 (正式写手)

单胞时候,原子数为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的回帖
查看全部 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考博] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:54 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 07:55 by lotyj5cz79
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:08 by lotyj5cz79
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 3/150 2026-02-18 06:53 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +4 瞬息宇宙 2026-02-15 5/250 2026-02-15 22:28 by baiboxie
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见