24小时热门版块排行榜    

查看: 3559  |  回复: 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的回帖

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 ]
8楼2009-06-09 15:48:49
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 【0856】化学工程(085602)313 分,本科学科评估A类院校化学工程与工艺,诚求调剂 +7 小刘快快上岸 2026-03-11 8/400 2026-03-17 16:57 by ruiyingmiao
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 281求调剂 +9 Koxui 2026-03-12 11/550 2026-03-13 20:50 by Koxui
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
信息提示
请填处理意见