24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2986  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wangjiang12

铁虫 (小有名气)

[求助] VASP计算六方BN的缺陷时时进行偶极修正, 即dipole corrections的问题已有1人参与

各位VASP的高手
       大家好!小弟是VASP的新手,最近在计算六方BN的缺陷,总共196个原子,计算时觉得应该要加上偶极修正,即dipole corrections。所以在INCAR里面尝试分别设置 IDIPOL = 1,2,3,4进行计算。(没有设置DIPOL和LDIPOL,因为不是非常懂)
       计算结果如下:
  For IDIPOL = 1:
  DIPCOR: dipole corrections for dipol
  direction  1 min pos    76,
dipolmoment        -164.963149      1.868515      0.000000 electrons x Angstroem
Tr[quadrupol]    -21343.950543
energy correction for charged system         0.000000 eV
dipol+quadrupol energy correction          309.867790 eV
added-field ion interaction          0.000000 eV  (added to PSCEN)

For IDIPOL = 2:
DIPCOR: dipole corrections for dipol
direction  2 min pos   137,
dipolmoment           0.384939    -34.809948      0.000000 electrons x Angstroem
Tr[quadrupol]    -60708.131283
energy correction for charged system         0.000000 eV
dipol+quadrupol energy correction           19.814583 eV
added-field ion interaction          0.000000 eV  (added to PSCEN)

For IDIPOL = 3:
DIPCOR: dipole corrections for dipol
direction  3 min pos    81,
dipolmoment           0.000000      0.000000      0.000000 electrons x Angstroem
Tr[quadrupol]    -79324.968578
energy correction for charged system         0.000000 eV
dipol+quadrupol energy correction            0.000000 eV
added-field ion interaction          0.000000 eV  (added to PSCEN)

For IDIPOL=4:
DIPCOR: dipole corrections for dipol
direction  1 min pos    87, direction  2 min pos  112, direction  3 min pos    97,
dipolmoment           3.159018     94.771489      0.002912 electrons x Angstroem
Tr[quadrupol]     -4892.150186
energy correction for charged system         0.000000 eV
dipol+quadrupol energy correction           47.526385 eV
added-field ion interaction          0.000000 eV  (added to PSCEN)

IDIPOL = 1的计算结果比其他的大很多,不知道该采用哪个,请高手指导我如何正确进行偶极修正的计算吧。

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

wangjiang12

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by golddoushi at 2014-03-11 13:05:05
原则上都不能用,非立方结构dipole-dipole作用可以达到monopole-monopole 的30%以上,quardpole-dipole作用在半导体有缺陷的情况下也不可忽略
最重要的是,这里得出的数都是认为介电常数=1的,手册上强调了好几次, ...

非常感谢您的回复,看起来过程好复杂啊。那您认为如果我不做dipole corrections, 直接对比不同缺陷排布下的能量大小,可以吗
3楼2014-03-11 15:17:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangjiang12

铁虫 (小有名气)

引用回帖:
4楼: Originally posted by golddoushi at 2014-03-11 15:57:47
通常情况约有 0.1eV ~ 1eV 数量级的误差。如果你的材料介电常数较小,误差还会更大。
dipole - dipole 作用的能量正比于 1/L3, monopole - monopole 正比于1/L,quardpole - dipole正比 1/L^5 或者更高次方。他们 ...

哦,好的,我明白了,谢谢。对了,请问您一下,除了dipole corrections外,您觉得用VASP计算六方NB还有什么其他需要修正或者注意的地方吗?
5楼2014-03-11 16:28:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangjiang12

铁虫 (小有名气)

引用回帖:
6楼: Originally posted by golddoushi at 2014-03-11 18:34:51
具体不太清楚,我见过有人用半芯态的赝势算的比较准...

半芯态的赝势是指什么啊?

谢谢。
7楼2014-03-12 10:11:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangjiang12 的主题更新
信息提示
请填处理意见