24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1312  |  回复: 3
当前主题已经存档。

leoaqu

金虫 (正式写手)

[交流] 【求助】VASP怎么模拟带电体系(charged system)

各位大侠请赐教:
      我看文献里面总经常提给一个体系(比如说2*2*2的ZnO supercell)带上n
个(0.5e, 1e,。。。)电荷,然后计算这个带电体系的性质(例如磁性),文章中一般会提到同时引入一个带正电的compensating background来达到保持体系电中性。
    我的问题是:
         一,我觉得这个compensating background是VASP自动引入的吧,就是说你只要在INCAR中设置好NELECT=q+n 就可以了吧?
         二,体系带电了,我只想计算磁性,不计算别的,还需要把IDIPOL tag打开吗?要是打开的话,我需要设置IDIPOL=多少?
         三,体系带电之后,体系的能量怎么计算?

PS:看VASP手册上面的Monopole, Dipole and Quadrupole corrections 部分没看懂。

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

gavinliu7390

木虫 (著名写手)

叶落鹰飞

★ ★ ★ ★ ★ ★ ★
yjr(金币+1,VIP+0):谢谢! 8-18 11:49
leoaqu(金币+2,VIP+0):详细 8-18 13:04
leoaqu(金币+1,VIP+0): 8-20 20:11
leoaqu(金币+3,VIP+0): 9-12 15:15
first, you think it is right.
second, IDIPOL= direction of electric field. 1 is x direction , 2 is y,3 is z.
you also should set EFIELD= xxx (units eV/A)
LDIPOL=.TRUE.
third, do it over, you can find free energy of OUTCAR .
真理是一点点接近的!
2楼2009-08-18 10:36:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hlshi

金虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★
leoaqu(金币+2,VIP+0): 8-18 13:04
yjr(金币+1,VIP+0):谢谢!! 8-18 15:32
leoaqu(金币+2,VIP+0): 8-20 20:11
leoaqu(金币+4,VIP+0): 9-12 15:15
vasp 计算带电体系时:

先grep NELECT OUTCAR (中性的OUTCAR)得到中性体系的价电子数n
然后
0.5e 时NELECT=n+0.5
0.5h 时NELECT=n-0.5
然后再驰誉,使用力收敛
静态计算得到能量,计算时一般不用打开你提到的那些tag
可以参考 JOURNAL OF APPLIED PHYSICS 106, 023910 (2009)

[ Last edited by hlshi on 2009-8-18 at 13:31 ]
3楼2009-08-18 12:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chuzhaonan

铁杆木虫 (著名写手)

第一性原理计算大师

★ ★ ★ ★ ★ ★
leoaqu(金币+2,VIP+0): 8-18 21:44
leoaqu(金币+1,VIP+0): 8-20 20:11
leoaqu(金币+3,VIP+0): 9-12 15:15
支持三楼,如果计算在外电场作用下的性质需要考虑IDIPOL , LDIPOL and EFIELD
玩第一性原理就如看大片
4楼2009-08-18 16:37:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 leoaqu 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见