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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】请指教VASP 总能，DOS，band-structure计算步骤
101/1返回列表
查看: 1444  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 08工学 309分求调剂 +3 Yin DY 2026-04-08 3/150 2026-04-09 16:24 by wp06
[考研] 288求调剂，一志愿华南理工大学071005 +15 ioodiiij 2026-04-08 15/750 2026-04-09 15:54 by gdupt2023
[考研] 求调剂 +3 猪肉墩粉条cc 2026-04-08 4/200 2026-04-09 10:05 by 猪肉墩粉条cc
[考研] 070300化学 求调剂 +10 73372112 2026-04-08 10/500 2026-04-09 07:59 by yuandd_2001
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 专硕310求调剂 +7 捞捞我…. 2026-04-04 8/400 2026-04-08 20:24 by yutian743
[考研] 0703调剂，一志愿天津大学319分 +23 haaaabcd 2026-04-05 26/1300 2026-04-08 16:19 by luoyongfeng
[考研] 305求调剂 +4 77Qi 2026-04-06 4/200 2026-04-07 20:06 by shanqishi
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 求调剂 +4 晟功? 2026-04-03 4/200 2026-04-04 21:58 by hemengdong
[考研] 11408 一志愿西电，277分求调剂 +4 zhouzhen654 2026-04-03 4/200 2026-04-04 18:10 by 猪会飞
[考研] 272求调剂 +4 松柏常青5 2026-04-03 4/200 2026-04-04 17:03 by babysonlkd
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
[考研] 286求调剂 +8 lim0922 2026-04-02 8/400 2026-04-03 20:19 by rzh123456
信息提示
关闭
请填处理意见
关闭