| 查看: 3976 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 有关零点能计算的一些帖子 | 零点能 |
» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关商家推荐: (我也要在这里推广)
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
4楼2009-06-08 23:06:37
huangyc
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2242.1
- 散金: 12
- 红花: 2
- 帖子: 411
- 在线: 312.5小时
- 虫号: 624241
- 注册: 2008-10-12
- 专业: 理论和计算化学
2楼2009-06-08 19:15:26
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
3楼2009-06-08 21:34:47
dummy1
木虫 (正式写手)
- 应助: 13 (小学生)
- 金币: 5841.7
- 散金: 32
- 红花: 1
- 帖子: 312
- 在线: 295.8小时
- 虫号: 270440
- 注册: 2006-08-06
- 专业: 理论和计算化学
|
谢谢各位的回复， 我是固定的晶格参数，仅仅做的粒子驰豫。 对于单胞，出现这么多的小的虚频也正常吗？因为我是希望想到准确点的频率，所以担心出现虚频不准确。各位有什么好的办法解决吗？ 对于超胞，我是优化的单胞然后用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
