24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1451  |  回复: 9
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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计算时候得到的能量稳定很多。

谢谢!
回复此楼

» 猜你喜欢

已阅   回复此楼   关注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的回帖
查看全部 10 个回答

化学小工

木虫 (正式写手)

★ ★ ★ ★
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的取值,应该反过来,无论哪个能量低,这里涉及准确性问题
有问题可以联系我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即可选取,不是单纯以能量高低判别。
3楼2009-10-08 11:08:19
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085404 22408 309分求调剂 +8 lzmk 2026-04-14 9/450 2026-04-15 13:23 by 黑科技矿业
[考研] 290调剂生物0860 +36 哇哈哈,。 2026-04-11 42/2100 2026-04-15 13:13 by 西北望—风沙
[考研] 求调剂学校 +13 不会吃肉 2026-04-13 15/750 2026-04-15 12:49 by 西北望—风沙
[考研] 289 分105500药学专硕求调剂(找B区学校) +3 白云123456789 2026-04-13 3/150 2026-04-15 12:42 by 西北望—风沙
[考研] 求调剂 +9 小聂爱学习 2026-04-11 13/650 2026-04-14 19:18 by Art1977
[考研] 材料085601调剂 +32 何润采123 2026-04-10 34/1700 2026-04-14 08:47 by 木木mumu~
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 085410 273求调剂 +10 X1999 2026-04-09 10/500 2026-04-12 09:24 by 逆水乘风
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 270求调剂 +14 杨乐369 2026-04-11 14/700 2026-04-11 20:16 by 蓝云思雨
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
[考研] 296求调剂 +6 汪!?! 2026-04-09 6/300 2026-04-11 11:25 by zhq0425
[考研] 求调剂 +6 archer.. 2026-04-09 8/400 2026-04-11 10:55 by zhq0425
[考研] 297求调剂 +9 Kwgyz 2026-04-09 9/450 2026-04-11 10:09 by zhq0425
[考研] 一志愿中科大070300化学,314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 材料专硕(0856) 339分求调剂 +9 哈哈哈鹅哈哈哈 2026-04-09 10/500 2026-04-09 20:01 by Orcid
信息提示
请填处理意见