| 查看: 1590 | 回复: 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 部分没看懂。 谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | soft study |
» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有5人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
★ ★ ★ ★ ★ ★ ★ ★ ★
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
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
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★ ★ ★ ★
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
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
chuzhaonan
铁杆木虫 (著名写手)
第一性原理计算大师
- 应助: 3 (幼儿园)
- 金币: 7049.8
- 散金: 120
- 红花: 3
- 帖子: 2282
- 在线: 129小时
- 虫号: 515926
- 注册: 2008-03-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2009-08-18 16:37:35
