| 查看: 3172 | 回复: 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的计算结果比其他的大很多，不知道该采用哪个，请高手指导我如何正确进行偶极修正的计算吧。 非常感谢 |
回复此楼» 猜你喜欢
085600 英一数二272求调剂
已经有12人回复
-
材料284求调剂，一志愿郑州大学英一数二专硕
已经有12人回复
-
化学，材料，环境类求调剂
已经有3人回复
-
材料化工调剂
已经有14人回复
-
接收调剂
已经有3人回复
-
化工京区271求调剂
已经有4人回复
-
295求调剂。一志愿报考郑州大学化学工艺学硕，总分295分
已经有5人回复
-
求调剂
已经有10人回复
-
0856化工专硕求调剂
已经有15人回复
-
303求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助完毕】硫酸根离子的结构优化，偶极修正及势能修正。
已经有20人回复
golddoushi
木虫 (正式写手)
- 应助: 36 (小学生)
- 金币: 4007.2
- 散金: 25
- 红花: 15
- 帖子: 337
- 在线: 447.7小时
- 虫号: 2546370
- 注册: 2013-07-15
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
wangjiang12: 金币+5, ★★★★★最佳答案 2014-03-11 18:31:05
感谢参与，应助指数 +1
wangjiang12: 金币+5, ★★★★★最佳答案 2014-03-11 18:31:05
|
原则上都不能用，非立方结构dipole-dipole作用可以达到monopole-monopole 的30%以上，quardpole-dipole作用在半导体有缺陷的情况下也不可忽略 最重要的是，这里得出的数都是认为介电常数=1的，手册上强调了好几次，要手动除掉缺陷之间（超胞超胞之间）介质的介电常数，然后手动加到能量上 所以你可以 1.计算介电常数，或者查到 2.测试dipole的位置，因为边界条件的关系任何方法都不能直接算出体系离子贡献的电矩（你假定这个东西是无限大的）。dipole取正电中心，到负电中心之间的各个点，直到校正项最小。也可能直接落在缺陷位置。此外要确保dipole落在电荷FFT的节点上，这个通过NGFX/NGFY/NGFZ控制，或移动缺陷到原点 3.测试超胞尺寸对偶极校正能量的影响，一直加大到能量变化不大（10%以内），取可以接受的最大的超胞 4.把计算结果手动除以介电常数，加到没有修正的能量上。 ================ 假如你对能量不感兴趣，或者你的电子结构会因为这个内部电矩有质的变化（比如铁磁和顺磁的区别） 那这些都不用做了，因为vasp永远得不到正确的解 |
2楼2014-03-11 13:05:05
wangjiang12
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 5
- 散金: 5
- 帖子: 104
- 在线: 20.8小时
- 虫号: 2317103
- 注册: 2013-03-04
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
3楼2014-03-11 15:17:00
golddoushi
木虫 (正式写手)
- 应助: 36 (小学生)
- 金币: 4007.2
- 散金: 25
- 红花: 15
- 帖子: 337
- 在线: 447.7小时
- 虫号: 2546370
- 注册: 2013-07-15
- 专业: 凝聚态物性 II ：电子结构
4楼2014-03-11 15:57:47
wangjiang12
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 5
- 散金: 5
- 帖子: 104
- 在线: 20.8小时
- 虫号: 2317103
- 注册: 2013-03-04
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
5楼2014-03-11 16:28:57
golddoushi
木虫 (正式写手)
- 应助: 36 (小学生)
- 金币: 4007.2
- 散金: 25
- 红花: 15
- 帖子: 337
- 在线: 447.7小时
- 虫号: 2546370
- 注册: 2013-07-15
- 专业: 凝聚态物性 II ：电子结构
6楼2014-03-11 18:34:51
wangjiang12
铁虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 5
- 散金: 5
- 帖子: 104
- 在线: 20.8小时
- 虫号: 2317103
- 注册: 2013-03-04
- 性别: GG
- 专业: 金属非晶态、准晶和纳米晶
7楼2014-03-12 10:11:11
golddoushi
木虫 (正式写手)
- 应助: 36 (小学生)
- 金币: 4007.2
- 散金: 25
- 红花: 15
- 帖子: 337
- 在线: 447.7小时
- 虫号: 2546370
- 注册: 2013-07-15
- 专业: 凝聚态物性 II ：电子结构
8楼2014-03-12 10:59:37
9楼2014-03-19 13:27:21
10