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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】请指教VASP 总能，DOS，band-structure计算步骤
101/1返回列表
查看: 1542  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

leoaqu

金虫 (正式写手)

[交流] 【求助】请指教VASP 总能，DOS，band-structure计算步骤

各位虫友，各位大虾：

    我对VASP计算体系总能，计算体系DOS以及band-structure的步骤有些疑问。以ZnO为例，如果我要计算2个Co掺杂的ZnO体系（332的supercell）的磁学性子，那么我就要比较FM和AFM结构的体系总能，所以，我的步骤是这样：
一、对Co-ZnO体系进行结构优化（FM和AFM结构分别优化），K点用小点，比如332。得到CONTCAR 和CHGCAR。
二、增加K点到443，用第一步得到的CONTCAR做位这步的POSCAR，NSW＝0， ICHARG=1。这样得到DOSCAR，CHGCAR和OUTCAR用来做结果分析，即DOSCAR－－DOS，CHGCAR－－spin density， OUTCAR－－系统总能E。
三、还是用第一步得到的CONTCAR 和CHGCAR来做位其实文件，并设置NSW＝0， ICHARG=11， KPIONTS沿某些特殊方向取。这样得到的EIGENVAL用来做band分析。

我一直是以为总能的计算与DOS的计算是一样的，我的方法对吗？有什么需要改进的请不吝赐教！

我的问题还有，就是我用我的方法，有的时候会遇到：
1.在结构优化的时候，本来是FM很稳定（能量低），而我做完第二步的时候，发现变成很AFM稳定了。
2.同一结构（e.g. FM 排列的Co－ZnO）在优化的时候，即KPT＝332的能量反而比KPT＝443的DOS计算时候得到的能量稳定很多。

谢谢!
回复此楼

» 猜你喜欢
1楼 2009-10-07 22:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)
★ ★ ★ ★
leoaqu(金币+2,VIP+0): 10-8 22:34
leoaqu(金币+1,VIP+0): 1-4 18:35
leoaqu(金币+1,VIP+0): 1-4 18:36
严重建议看一下vasp的manu和handbook，呵呵
其实不麻烦，dos和能量的INCAR不一样，需要设置几个参数，而且设置的时候根据情况还要有点变化
另外kpoints里面，不能用能量衡量ktp的取值，应该反过来，无论哪个能量低，这里涉及准确性问题
赞一下(6人) 回复此楼
有问题可以联系我QQ676869820
2楼2009-10-08 10:04:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linggang87

银虫 (小有名气)
★ ★ ★ ★ ★
leoaqu(金币+2,VIP+0): 10-8 22:34
leoaqu(金币+2,VIP+0): 1-4 18:35
leoaqu(金币+1,VIP+0): 1-4 18:36
楼主所说的
“2.同一结构（e.g. FM 排列的Co－ZnO）在优化的时候，即KPT＝332的能量反而比KPT＝443的DOS计算时候得到的能量稳定很多”。
看来是把能量的高低作为稳定的标准了，实际上对K点的测试，是以K点位横坐标，能量为纵坐标做曲线，曲线趋于平缓时的K即可选取，不是单纯以能量高低判别。
赞一下(6人) 回复此楼
3楼2009-10-08 11:08:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoaqu

金虫 (正式写手)
我有很仔细的看vasp的manu，但是还是不太清楚。我的INCAR是这样的：
一，结构优化

INCAR file：
ENCUT=520
ISTART=0
ICHARG=2     
INIWAV=1
GGA=PE
LWAVE= .FALSE.
LCHARG=.TRUE.         
NSW=100               
NELM=100
IBRION=2               
ISIF=3                  
EDIFF=1.0E-04        
EDIFFG = -0.01
ISMEAR= 0               
SIGMA = 0.1
LORBIT=11   
LREAL=Auto

KPONITS file：
Monkhorst Pack
0
Monkhorst Pack
3 3 2
0 0 0

二 DOS和系统能量
INCAR file：
ENCUT=520
ISTART=0
ICHARG=1     
INIWAV=1
GGA=PE
LWAVE= .FALSE.
LCHARG=.TRUE.         
NSW=0               
NELM=100
IBRION=－1               
ISIF=0                  
EDIFF=1.0E-06        
＃EDIFFG = -0.01
ISMEAR= 0               
SIGMA = 0.1
LORBIT=11   
LREAL=Auto
EMAX=20                 
EMIN=-20              
NEDOS=2000

KPONITS file：
Monkhorst Pack
0
Monkhorst Pack
4 4 3
0 0 0

请问我的INCAR需要做什么样的改动？再就是能量计算跟这个DOS计算是一样的吗？

关于KPT的问题，我是觉得对332的supercell而言，332的应该已经可以收敛（文献里面的也是用这个kpt），所以，在DOS计算的时候就增加到443（也是别人跟我说的）。
请帮我看看，谢谢。
引用回帖:
Originally posted by 化学小工 at 2009-10-8 10:04:
严重建议看一下vasp的manu和handbook，呵呵
其实不麻烦，dos和能量的INCAR不一样，需要设置几个参数，而且设置的时候根据情况还要有点变化
另外kpoints里面，不能用能量衡量ktp的取值，应该反过来，无论哪个能量 ...

赞一下 回复此楼
4楼2009-10-08 12:04:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoaqu

金虫 (正式写手)
谢谢回帖。
是的，我知道那个K点测试是怎么怎么回事，但是当K点达到收敛之后，你继续增高K点之后，能量是应该降低，或者基本一致的才对的。而我之所以对这个443的K点得到的能量较高有疑惑，是因为按照文献上讲332应该已经收敛，所以，443得到的能量应该在那个平滑的线段上，不应该有这么大的跳跃吧。

PS：我没有做K点的详细测试，我只是按照文献里面的值来取的，我觉得好多人也差不多是这样取K点的吧。
谢谢
引用回帖:
Originally posted by linggang87 at 2009-10-8 11:08:
楼主所说的
“2.同一结构（e.g. FM 排列的Co－ZnO）在优化的时候，即KPT＝332的能量反而比KPT＝443的DOS计算时候得到的能量稳定很多”。
看来是把能量的高低作为稳定的标准了，实际上对K点的测试，是以K点位横坐 ...

赞一下 回复此楼
5楼2009-10-08 12:12:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)
★ ★
leoaqu(金币+1,VIP+0): 10-8 22:35
leoaqu(金币+1,VIP+0): 1-4 18:36
引用回帖:
Originally posted by leoaqu at 2009-10-8 12:12:
谢谢回帖。
是的，我知道那个K点测试是怎么怎么回事，但是当K点达到收敛之后，你继续增高K点之后，能量是应该降低，或者基本一致的才对的。而我之所以对这个443的K点得到的能量较高有疑惑，是因为按照文献上讲33 ...

问一下，为什么你认为K点高了以后能量应该降低的？
赞一下(3人) 回复此楼
有问题可以联系我QQ676869820
6楼2009-10-08 22:12:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)
★ ★ ★ ★
leoaqu(金币+2,VIP+0): 10-8 22:34
leoaqu(金币+1,VIP+0): 1-4 18:34
leoaqu(金币+1,VIP+0): 1-4 18:36
目前看你的设置没什么问题好像，你自己感觉到有什么问题了吗？只有band计算的KPOINTs不知道你怎么设置的，呵呵，其他的好像没什么大问题，大家来看看呢？
赞一下(5人) 回复此楼
有问题可以联系我QQ676869820
7楼2009-10-08 22:22:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoaqu

金虫 (正式写手)
我基本不做band计算，主要看能量和DOS，Band的KPT我是可以从网上查到的。如果体系的能量可以用DOS计算得到的那个能量的话，我就没有问题了。谢谢。
引用回帖:
Originally posted by 化学小工 at 2009-10-8 22:22:
目前茨愕纳柚妹皇裁次侍夂孟瘢阕约焊芯醯接惺裁次侍饬寺穑恐挥衎and计算的KPOINTs不知道你怎么设置的，呵呵，其他的好像没什么大问题，大家来看看呢？

赞一下 回复此楼
8楼2009-10-08 22:29:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leoaqu

金虫 (正式写手)
我没有说一定降低。我是说如果当332的kpt已经收敛了的话，那么用443的kpt的时候，能量即使不继续降低，也不应该有太大的波动吧，否则，就不能说332已经收敛了吧。那么你觉得那？

我也不是每个case都这样，只是遇到一两次这样的。
引用回帖:
Originally posted by 化学小工 at 2009-10-8 22:12:

问一下，为什么你认为K点高了以后能量应该降低的？

赞一下 回复此楼
9楼2009-10-08 22:33:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

化学小工

木虫 (正式写手)
★ ★ ★ ★ ★ ★ ★
zdhlover(金币+2,VIP+0):欢迎常来交流 10-9 00:33
leoaqu(金币+3,VIP+0): 10-11 16:48
leoaqu(金币+2,VIP+0): 1-4 18:34
是不一样的，看你的lx，晶格常数不一样以后K也要相应调整
另外低K应该说明周期性的近似降低了，所以它收敛不代表高k的时候也收敛，所以我们经常采用低K优化结构，高k算dos。
赞一下(15人) 回复此楼
有问题可以联系我QQ676869820
10楼2009-10-08 22:46:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leoaqu 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:50 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 9zw1gfyi7x 2026-08-01 8/400 2026-08-01 18:50 by Kddn1pB6pAzX
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 9zw1gfyi7x 2026-08-01 7/350 2026-08-01 18:43 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 6/300 2026-08-01 18:37 by Kddn1pB6pAzX
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 9zw1gfyi7x 2026-08-01 9/450 2026-08-01 18:30 by Kddn1pB6pAzX
[论文投稿] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 11/550 2026-08-01 18:17 by Kddn1pB6pAzX
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 9zw1gfyi7x 2026-07-31 9/450 2026-08-01 18:10 by Kddn1pB6pAzX
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 9zw1gfyi7x 2026-07-31 10/500 2026-08-01 18:03 by Kddn1pB6pAzX
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 10/500 2026-08-01 17:51 by Kddn1pB6pAzX
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 7/350 2026-08-01 17:51 by Kddn1pB6pAzX
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:50 by Kddn1pB6pAzX
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 jzg3xk1po8 2026-07-31 9/450 2026-08-01 17:43 by Kddn1pB6pAzX
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 jzg3xk1po8 2026-07-31 12/600 2026-08-01 17:43 by Kddn1pB6pAzX
[基金申请] 面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！ +16 低垂的野花 2026-07-31 22/1100 2026-08-01 16:39 by Cathyau
[基金申请] 微信指数没变化，科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[基金申请] 今年的WR进展到哪一步了？ +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了，我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
信息提示
关闭
请填处理意见
关闭