24小时热门版块排行榜    

查看: 2121  |  回复: 5

傻孩阿布

铁虫 (初入文坛)

[求助] vasp小白一枚,求大神指教,单原子能量计算问题,内附计算文件

vasp小白一枚,计算形成能,需要各单原子的化学势,于是我就对各单原子的能量进行了计算,计算文件如下。还请大神指教,那些地方不妥当。
incar
ENCUT =500

ISPIN=2

MAGMOM=1*3

ISTART=0

ISPIN=2

IALGO=38

GGA=PE

LREAL=.F.

ISMEAR=0

SIGMA=0.1

ICHARG=2

ISIF=3

IBRION=2

NSW=30

EDIFF=1E-6

EDIFFG=-1E-4
LWAVE  =F
poscar
atom

1.00

10.000000000000000       0.000000000000000       0.000000000000000
0.000000000000000      10.000000000000000       0.000000000000000
0.000000000000000       0.000000000000000      10.000000000000000

Cr

1

Direct

0.5000000000000000   0.5000000000000000   0.5000000000000000
K点
Automatic

0

Gamma

5 5 5

0 0 0
。一个是请教大神,上述文件的写入是否正确。第二个是计算时是否要考虑自旋极化的问题,因为发现加极化与否,对结果的影响还是很大的。先谢过各位大侠了。在线等。谢谢。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

Germanman

金虫 (小有名气)

【答案】应助回帖


sunyang1988: 金币+1, 谢谢交流 2013-10-11 15:55:31
对于单个原子的计算,INCAR不需要这么多参数,设置几个必要的,其他的用默认即可。自旋极化采取默认值。
2楼2013-10-11 14:25:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egg19900805

新虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2013-10-13 09:17:24
形成能是计算各原子行程所需的能量吗?
像是中:1/5E(Cr2O3)-1/2E(O2)= E(Cr)(1/5是因为有5个原子)
那你就只需要把Cr2O3的O2平衡取final energy
再将其相减就能够得到
3楼2013-10-11 17:15:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

米色空白

金虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2013-10-13 09:17:32
可以参考一下VASP官方举办的研讨会的讲义,还有里面的一些案例,对初学者挺有帮助。
教程:http://www.vasp.at/index.php/documentation
案例:http://cms.mpi.univie.ac.at/wiki/index.php/The_VASP_Manual
4楼2013-10-11 23:41:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Brilliantly

捐助贵宾 (正式写手)

厉害!

[ 发自小木虫客户端 ]
加油!
5楼2013-10-13 13:26:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鹰羽龙

木虫 (著名写手)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流 2013-10-15 09:30:05
据说SPIN必须要开,K点只用gamma点,晶胞3边不等长,具体效果如何还不知道,我正在测试
6楼2013-10-13 14:44:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 傻孩阿布 的主题更新
信息提示
请填处理意见