24小时热门版块排行榜    

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

juliantiger

铜虫 (初入文坛)

[求助] vasp计算态密度结果出问题

我计算的结构是钛酸锶,在计算态密度后,发现锶的分波态密度很小,几乎为零,与我看的文献上的完全不一样,文献上的都是在-15ev处,有很高的峰值,请问这是怎么回事?是不是我哪里设置有问题,下面是我的INCAR:
System = STO
ISTART = 1
ICHARG = 11
PREC = Accurate
LREAL = .FALSE.
EDIFF = 1E-5
EDIFFG = 2E-5
IBRION = -1
ENCUT = 400
EMAX = 10
EMIN = -20
POTIM = 0.05
ISMEAR = -5
LORBIT = 1
NEDOS = 501
RWIGS = 2.138 1.476 0.741
我也使用了LORBIT=11,不设置RWIGS,出来的结果一样。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

没什么说的
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bjmaowei25

铜虫 (小有名气)

【答案】应助回帖

★ ★ ★
御剑江湖: 金币+3, 谢谢回帖 2012-06-15 09:00:46
LZ 对系统可能没有事先进行晶格优化,在你的INCAR 中
设置ICHARG= 11, 就是利用优化后的电荷密度。
所以你有两种改法:
1. 直接把 ICHARG = 11 ,改为ICHARG= 2. ISTART = 1 改为ISTART = 0
其实如果你仔细看侯老师的入门手册,当设置 ICHARG = 11 是,
利用了前面自洽的电荷密度的结果。再运算看看。
2. 先进行 自洽计算,此时ICHARG= 2, 得到优化的晶格参数,用得到的平衡晶格常数代入POSCAR;
修改 ISTART = 1 ,ICHARG = 11, 也就是你上面INCAR内容不变。

另外,有一点提醒LZ, 你的O radius 好像不对哦,
追求真理
6楼2012-06-14 22:17:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

likai0106

金虫 (小有名气)

【答案】应助回帖

★
感谢参与,应助指数 +1
WDD880227: 金币+1, 感谢您的提示 2012-05-31 15:20:17
楼主注意观察一下-15ev处是不是d电子提供的?如果是那么可能是你分波的脚本有问题,检查一下总态密度,如果没有问题就看看脚本吧,我遇到过很多脚本在处理d电子时写的有问题。
2楼2012-05-31 07:47:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juliantiger

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by likai0106 at 2012-05-31 07:47:50
楼主注意观察一下-15ev处是不是d电子提供的?如果是那么可能是你分波的脚本有问题,检查一下总态密度,如果没有问题就看看脚本吧,我遇到过很多脚本在处理d电子时写的有问题。

态密度的-15的地方也没有,今天我换了赝势,-15的地方是有了,但是很异常,非常的大。
没什么说的
3楼2012-05-31 12:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

likai0106

金虫 (小有名气)

引用回帖:
3楼: Originally posted by juliantiger at 2012-05-31 12:29:23
态密度的-15的地方也没有,今天我换了赝势,-15的地方是有了,但是很异常,非常的大。...

楼主的计算中没有加自旋啊,不知道文献中有没有给出自旋部分的计算。还有楼主对体系进行了优化么?原子坐标和晶格常数。
4楼2012-05-31 15:47:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:45 by Z9YWQ5EAO3qp
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:26 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +4 tianxiaoxian 2026-09-28 7/350 2026-09-30 08:34 by 88817753
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
信息提示
请填处理意见