24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
散金: 111
红花: 2
帖子: 322
在线: 102.6小时
虫号: 860161
注册: 2009-09-29
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】VASP，算的是CaFe2As2，结果DOSCAR中全部是零已有6人参与

老师让我先算CaFe2As2单胞的DOS，我算的结果根本就不对，DOSCAR中数据基本上都是0。大致过程是参看http://old.blog.edu.cn/user3/xiaomagh/archives/2007/1855608.shtml的《VASP计算DOS的步骤》来的。
这是INCAR：

CODE:

System = CaFe2As2

PREC = accurate

ISPIN = 1

ISTART = 0

ICHARG = 2

EDIFF = 10E-4

ISIF = 2

ISMEAR = -5

SGIMA = 0.2

这是POSCAR：

CODE:

CaFe2As2

1

912680000000000 0.000000000000000 0.000000000000000

000000000000000 3.912680000000000 0.000000000000000

000000000000000 0.000000000000000 11.667450000000001

4 4

Direct

0000000000000000   0.0000000000000000   0.0000000000000000

5000000000000001   0.5000000000000001   0.5000000000000000

5000000000000001   0.0000000000000000   0.2500000000000000

0000000000000000   0.5000000000000001   0.2500000000000000

5000000000000001   0.0000000000000000   0.7500000000000000

0000000000000000   0.5000000000000001   0.7500000000000000

0000000000000000   0.0000000000000000   0.3581399917602537

5000000000000001   0.5000000000000001   0.8581399917602538

0000000000000000   0.0000000000000000   0.6418600082397463

5000000000000001   0.5000000000000001   0.1418600082397463

这是KPOINTS：

CODE:

Automatic  generation

0

Monhkorst Pack

11 11 11

0.0 0.0 0.0

POTCAR使用的是pot中的势库文件。
DOSCAR中，给出了K点坐标之后，下面的全部为0。
请问，这个问题出在哪里了？是INCAR文件哪里设置不对么？
谢谢

[ Last edited by gleerat on 2010-4-30 at 18:20 ]

回复此楼

明察、慎思、笃行

1楼2010-04-30 17:53:57

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

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
散金: 111
红花: 2
帖子: 322
在线: 102.6小时
虫号: 860161
注册: 2009-09-29
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

Originally posted by 锐利的碎片 at 2010-04-30 19:12:30:
INCAR很有问题啊
RWIGS没给，计算要读自洽的电荷密度。ISMEAR=-5不用设SIGMA。加了自旋要给原子磁矩。
[ Last edited by 锐利的碎片 on 2010-4-30 at 19:13 ]

我没有加自旋，ISPIN＝1。我在静态自洽计算中给了RWIGS的值，这是我做完计算时的INCAR：

CODE:

System = CaFe2As2

PREC = accurate

ISPIN = 1

ISTART = 0

#ICHARG = 11

#ISIF = 2

ISMEAR = 1

SGIMA = 0.2

#NSW = 0

#IBRION = 2

#NELMDL = -6

#LORBIT = 2

#RWIGS = 2.5 2.17 1.92