24小时热门版块排行榜    

查看: 4020  |  回复: 31

cumt蝈蝈

木虫 (正式写手)

博众彩之长

★ ★
identation: 金币+2, 奖励一下,多调试调试会有更多经验! 2013-09-21 10:48:05
引用回帖:
7楼: Originally posted by 漫天飘雪 at 2013-09-16 21:59:49
嗯,你可以对比一下...

我对比了下,把ISMEAR设为1和-5,他们的能量差还不到0.02。差别不大,应该可以忽略吧?
君子生非异也,善假于物也
11楼2013-09-17 08:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emilyoyang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
cumt蝈蝈: 金币+10, ★★★很有帮助 2013-09-17 13:45:10
identation: 金币+2, 应助指数+1, 谢谢应助 2013-09-21 10:48:34
引用回帖:
10楼: Originally posted by cumt蝈蝈 at 2013-09-16 23:10:20
我算的体系原子数是45-50不定。
我设置下ALGO 试试。
ismear=-5没有错误提示啊,你在manual的那部分看到的“z方向只有1个k点貌似不能用-5 manual”我没看到过?...

http://cms.mpi.univie.ac.at/vasp ... ethod_smearing.html
看一下这段内容,我是这么理解的。不知道对不对。你自己理解看一下
我很少用-5(DOS除外)一般都用0 or 1,然后再加上一个合适的sigma值
12楼2013-09-17 10:00:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbing9002

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
cumt蝈蝈: 金币+5, 有帮助 2013-09-20 10:52:50
identation: 金币+1, 谢谢应助 2013-09-21 10:48:52
把ISMEAR改为0
收敛算法用Fast
K点不用这么多,用5*5*1就足够了
穷则独善其身,达则兼济天下。
13楼2013-09-17 16:02:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cumt蝈蝈

木虫 (正式写手)

博众彩之长

引用回帖:
12楼: Originally posted by emilyoyang at 2013-09-17 10:00:04
http://cms.mpi.univie.ac.at/vasp/vasp/Number_k_points_method_smearing.html
看一下这段内容,我是这么理解的。不知道对不对。你自己理解看一下
我很少用-5(DOS除外)一般都用0 or 1,然后再加上一个合适的s ...

没看出来,感觉没有提z方向只有1个K点的情况。还是意在言外?
君子生非异也,善假于物也
14楼2013-09-20 10:48:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cumt蝈蝈

木虫 (正式写手)

博众彩之长

引用回帖:
13楼: Originally posted by lbing9002 at 2013-09-17 16:02:55
把ISMEAR改为0
收敛算法用Fast
K点不用这么多,用5*5*1就足够了

好像以前在哪看过,算过渡金属K点要求比较高,半导体、绝缘体低些,我这也是过渡金属啊?
君子生非异也,善假于物也
15楼2013-09-20 10:52:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbing9002

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2013-09-22 08:49:53
cumt蝈蝈: 金币+5, 有帮助 2013-09-22 17:02:18
引用回帖:
15楼: Originally posted by cumt蝈蝈 at 2013-09-20 10:52:27
好像以前在哪看过,算过渡金属K点要求比较高,半导体、绝缘体低些,我这也是过渡金属啊?...

我也算过渡金属,我之前做过k点测试,对于块体材料的小原胞,3*3*3的k点,能量就已经收敛了,取的再多,能量变化非常小。所以说表面取5*5*1已经算多的了。只要不算光学性质就够用的了。
穷则独善其身,达则兼济天下。
16楼2013-09-21 19:32:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李还乱

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
cumt蝈蝈: 金币+5, 有帮助 2013-09-22 17:02:11
franch: 金币+2, 谢谢回帖交流 2013-09-23 20:55:09
你的体系50左右个原子,也不算小了,4核也不算多,2天时间还觉得长啊,如果真想减少时间,精度调低,另外K点网格也太多了,你可以先用小体系测试下,3*3*3以后能量就变化很小了。ISMEAR=-5时 ,SIGMA=0.18这个参数应该没用吧
17楼2013-09-22 10:51:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cumt蝈蝈

木虫 (正式写手)

博众彩之长

引用回帖:
16楼: Originally posted by lbing9002 at 2013-09-21 19:32:18
我也算过渡金属,我之前做过k点测试,对于块体材料的小原胞,3*3*3的k点,能量就已经收敛了,取的再多,能量变化非常小。所以说表面取5*5*1已经算多的了。只要不算光学性质就够用的了。...

我也做过k带你测试‘k点为7*7*7、9*9*9、11*11*11、13*13*13、15*15*15时ΔE分别为:0.03828、0.01168、0.00373、0.000726、0.000328,我取在1E-3左右,所以选的13*13*13。
君子生非异也,善假于物也
18楼2013-09-22 16:59:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cumt蝈蝈

木虫 (正式写手)

博众彩之长

引用回帖:
17楼: Originally posted by 李还乱 at 2013-09-22 10:51:26
你的体系50左右个原子,也不算小了,4核也不算多,2天时间还觉得长啊,如果真想减少时间,精度调低,另外K点网格也太多了,你可以先用小体系测试下,3*3*3以后能量就变化很小了。ISMEAR=-5时 ,SIGMA=0.18这个参数应 ...

ISMEAR为-5时,SIGMA是应该不起作用的,我没去掉。我也不清楚精度为多少合适,k点测试时取得1E-3。
君子生非异也,善假于物也
19楼2013-09-22 17:01:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yayanoba

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
cumt蝈蝈: 金币+5, 有帮助 2013-09-23 18:43:15
franch: 金币+2, 谢谢回帖交流 2013-09-23 20:55:18
我觉得你用ISMEAR=-5的话,更多是算绝缘体和半导体比较准确,对于金属来说,你可以设ISMEAR=0,然后测试一下sigma,而且是不是可以把encut取大一些呢,或者你可以先用一个比较小的k点,然后得到它的CHGCAR,然后再后续用大k点计算的时候读入小k点时的CHGCAR,这样速度会不会快点!
20楼2013-09-23 13:27:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cumt蝈蝈 的主题更新
信息提示
请填处理意见