24小时热门版块排行榜    

查看: 2130  |  回复: 28
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lison4

金虫 (小有名气)

[求助] 初学VASP 测试2*2*1的ZnO超胞的ENCUT但是得不到结果 已有2人参与

下面是计算结果和所用的INCAR   
SYSTEM=ZnO
IBRION=2
NSW=100
ENCUT=360.0 eV
NELM=60
ISPIN=2
SIGMA  =   0.05
ISMEAR=0              
LORBIT=11            
ALGO  = 48
NSIM  =4
LPLANE = .TRUE.
NPAR = 12
LSCALU=.FALSE.

POSCAR的内容是直接由materials studio建好模型后导出的,K点为4*4*4自动生成的MP特殊点,赝势均从GGA里的PW91内选取的;现在截断值已经取到700了总能还是没有很稳定,想请问一下这可能是什么原因导致的。谢谢!

初学VASP 测试2*2*1的ZnO超胞的ENCUT但是得不到结果
计算结果
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
lison4: 金币+20, ★★★很有帮助 2016-11-06 10:09:43
引用回帖:
7楼: Originally posted by lison4 at 2016-11-04 21:24:45
谢谢您啦,我又更改了INCAR如下,然后测试了ENCUT和K如下  您帮我看看有问题吗 ,我发现总能还是没有很稳定啊

SYSTEM=ZnO
ISTART=0
ICHARG=2
ENCUT=360.0 eV
ISMEAR=-5           
PRE ...

做截断测试时,自己设置K点网格,建议用肯定能够收敛的密集网格算,比如对于氧化锌原胞,取个19*19*17这样的,然后逐步提高截断能。得到截断能后,用这个截断能放到超胞里手动一点点加k点,直到k点也收敛。

» 本帖已获得的红花(最新10朵)

思想重于技巧,内涵重于表象
8楼2016-11-04 23:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 29 个回答

wyding

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
蛋蛋小童鞋: 金币+2, 谢谢参与,小木虫鼓励积极回答问题者 2016-11-04 17:56:11
lison4: 金币+10, ★★★很有帮助 2016-11-06 10:11:23
NPAR = 12
设置不合理,不知道是否按VASP推荐的设置,不过不影响计算结果。可以看这个来学习   
http://blog.sciencenet.cn/blog-588243-741859.html

对于你遇到的问题
2*2*1的ZnO的单胞,4*4*4的k点够不够?
ISMEAR等于0合适?半导体绝缘体不是要-5吗?
一般建议找下O和Zn各自的Enmax,POTCAR里面有。取最大的,在这个附近测试。

» 本帖已获得的红花(最新10朵)

向强者看齐,即使高度不够也是昂首挺胸!
2楼2016-11-04 17:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

没必要这么复杂,截断能是跟着赝势走的,拿个HCO的氧化锌单胞做做足够了

» 本帖已获得的红花(最新10朵)

思想重于技巧,内涵重于表象
3楼2016-11-04 17:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lison4

金虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by wyding at 2016-11-04 17:23:03
NPAR = 12
设置不合理,不知道是否按VASP推荐的设置,不过不影响计算结果。可以看这个来学习   
http://blog.sciencenet.cn/blog-588243-741859.html

对于你遇到的问题
2*2*1的ZnO的单胞,4*4*4的k点够不够? ...

谢谢,我再试试其他的

发自小木虫IOS客户端
4楼2016-11-04 18:06:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +8 要乐观耀哥 2026-08-10 23/1150 2026-08-13 00:35 by 要乐观耀哥
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 好奇怪的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
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 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 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
信息提示
请填处理意见