应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp计算态密度结果出问题
121/2返回列表12下一页
查看: 2722  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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，出来的结果一样。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
没什么说的
1楼 2012-05-30 22:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

juliantiger

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by likai0106 at 2012-05-31 15:47:53
楼主的计算中没有加自旋啊，不知道文献中有没有给出自旋部分的计算。还有楼主对体系进行了优化么？原子坐标和晶格常数。...

我的没有加自旋，文献中也没加自旋，我在不停测试我的INCAR，发现NEBAND的大小对结果影响很大，这个不知道什么原因？
赞一下 回复此楼
没什么说的
5楼2012-05-31 18:30:50
已阅   回复此楼   关注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的回帖

huangguiyang

新虫 (小有名气)

【答案】应助回帖

★ ★
御剑江湖: 金币+2, 谢谢回帖 2012-06-15 09:00:52
可能是ISMEAR=-5的原因。

如果k点不是很多的话，结果是有问题的，加起来都不等于总的。

改用
ISMEAR=0, SIGMA = 0.05

就不会再有问题了。
赞一下 回复此楼
7楼2012-06-14 23:59:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liufeila

新虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢回帖 2012-06-17 18:45:19
ISMEAR=-5的话，K点要取得足够大，EDOS为1000，如果不行试试2000.
赞一下 回复此楼
8楼2012-06-15 09:27:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juliantiger

铜虫 (初入文坛)
引用回帖:
6楼: Originally posted by bjmaowei25 at 2012-06-14 22:17:11
LZ 对系统可能没有事先进行晶格优化，在你的INCAR 中
设置ICHARG= 11, 就是利用优化后的电荷密度。
所以你有两种改法：
1. 直接把 ICHARG = 11 ，改为ICHARG= 2. ISTART = 1 改为ISTART = 0
其实如果你仔细看侯 ...

前面的优化和静态都做了，还是有问题，O的radius是从potcar中提出来
赞一下 回复此楼
没什么说的
9楼2012-06-17 16:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

umklapp

银虫 (小有名气)

【答案】应助回帖

★ ★
fzx2008: 金币+2, 鼓励交流 2012-06-23 09:57:41
引用回帖:
9楼: Originally posted by juliantiger at 2012-06-17 16:45:03
前面的优化和静态都做了，还是有问题，O的radius是从potcar中提出来...

如果你的体系大，ISMEAR=-5 只会让你的DOS 看上去更peaky, 但是 weight 应该是没问题的。换用其他的smearing method 或者增加K 点 也只是让你的dos看上去更好看一些。
最好能贴出你的DOS；
既然你不relax 就不要设置EDIFFG = 2E-5； 而且如果你要优化结果，EDIFFG 也应该设置为负值；
赞一下 回复此楼
10楼2012-06-23 02:55:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 juliantiger 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:53 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 68gx6bh3ei 2026-07-09 3/150 2026-07-10 06:33 by pbuc7it0i8
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ph7bxbho8f 2026-07-09 3/150 2026-07-10 05:53 by j4frgf6zcw
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +4 al9t2f50x3 2026-07-09 4/200 2026-07-10 05:33 by j4frgf6zcw
[基金申请] b口会评 +6 当空照 2026-07-08 6/300 2026-07-09 20:49 by tonyliu0401
[基金申请] 信息学部会评时间 +3 xuel2011 2026-07-06 3/150 2026-07-09 19:05 by mollyzhang_2003
[基金申请] 祈祷青基必中 50+3 hadengry 2026-07-09 9/450 2026-07-09 17:37 by 文杰猪猪
[考研] 在职考研 +6 wu5d 2026-07-08 7/350 2026-07-09 16:40 by 化工小霸王呀
[论文投稿] MSER送审了还被拒稿 +3 S778950906 2026-07-08 5/250 2026-07-09 15:46 by S778950906
[基金申请] 关于如何从代码看上不上会 +11 且听虎啸 2026-07-09 14/700 2026-07-09 15:27 by libeny
[有机交流] 硝基还原 20+3 暮年飞雪 2026-07-07 5/250 2026-07-09 13:30 by tfang
[基金申请] 化学口会评？应用基础研究和基础研究，希望能分开 +3 Tide man 2026-07-08 3/150 2026-07-09 11:36 by dede5556
[基金申请] 今年E04面上 +15 布布和一二 2026-07-05 21/1050 2026-07-09 10:25 by 学员fOzRO9
[硕博家园] 大龄残疾硕士的一点执念 +14 3776737374 2026-07-06 22/1100 2026-07-09 06:07 by lovechizi
[基金申请] 生命口会评 +4 wuke100666 2026-07-07 4/200 2026-07-08 16:38 by 逍遥不够自在
[教师之家] 网络举报学术不端成为了男女之间情杀的手段了，真荒唐。。。 +3 zju2000 2026-07-04 3/150 2026-07-07 21:00 by hyx552200
[论文投稿] 刑法学论文投稿求助 +3 高源678 2026-07-04 5/250 2026-07-06 12:20 by 平天尺
[论文投稿] Journal of Environmental Chemical Engineering +3 雪花天涯 2026-07-05 3/150 2026-07-06 09:00 by bobvan
[基金申请] 最后一年，祈求好运 +3 pyw1992 2026-07-06 3/150 2026-07-06 08:53 by 苏_frp
信息提示
关闭
请填处理意见
关闭