| 5 | 1/1 | 返回列表 |
| 查看: 3213 | 回復: 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的計算結果比其他的大很多,不知道該采用哪個,請高手指導我如何正確進行偶極修正的計算吧。 非常感謝 |
鐵蟲 (小有名氣)
木蟲 (正式寫手)
|
原則上都不能用,非立方結構dipole-dipole作用可以達到monopole-monopole 的30%以上,quardpole-dipole作用在半導體有缺陷的情況下也不可忽略 最重要的是,這里得出的數都是認為介電常數=1的,手冊上強調了好幾次,要手動除掉缺陷之間(超胞超胞之間)介質的介電常數,然后手動加到能量上 所以你可以 1.計算介電常數,或者查到 2.測試dipole的位置,因為邊界條件的關系任何方法都不能直接算出體系離子貢獻的電矩(你假定這個東西是無限大的)。dipole取正電中心,到負電中心之間的各個點,直到校正項最小。也可能直接落在缺陷位置。此外要確保dipole落在電荷FFT的節(jié)點上,這個通過NGFX/NGFY/NGFZ控制,或移動缺陷到原點 3.測試超胞尺寸對偶極校正能量的影響,一直加大到能量變化不大(10%以內),取可以接受的最大的超胞 4.把計算結果手動除以介電常數,加到沒有修正的能量上。 ================ 假如你對能量不感興趣,或者你的電子結構會因為這個內部電矩有質的變化(比如鐵磁和順磁的區(qū)別) 那這些都不用做了,因為vasp永遠得不到正確的解 |
木蟲 (正式寫手)
鐵蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 313求調劑 +6 | 賣個關子吧 2026-03-31 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 調劑申請 +5 | 張張張張zy 2026-03-31 | 5/250 |
|
|
[考研] 085600材料與化工329分求調劑 +5 | 這是有太陽哇 2026-03-25 | 5/250 |
|
|
[考研] 看看我!301求調劑 081700 +3 | 一個小明同學 2026-03-24 | 3/150 |
|
|
[考研] 一志愿中海洋320化學工程與技術學碩求調劑 +8 | 披星河 2026-03-30 | 8/400 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調劑 +6 | 邱gl 2026-03-30 | 10/500 |
|
|
[考研] 083000環(huán)境科學與工程調劑,總分281 +4 | 橙子(勝意) 2026-03-30 | 4/200 |
|
|
[考研] 一志愿華東師范大學有機化學專業(yè),初試351分,復試被刷求調劑! +7 | 真名有冰 2026-03-29 | 8/400 |
|
|
[考研] 085601材料工程找調劑 +17 | oatmealR 2026-03-29 | 18/900 |
|
|
[考研] 322求調劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 310求調劑 +10 | 爭取九點睡 2026-03-30 | 10/500 |
|
|
[考研] 考研調劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
[考研] 求調劑 +10 | 張zz111 2026-03-27 | 11/550 |
|
|
[考研] 086000生物與醫(yī)藥調劑 +5 | Feisty。 2026-03-28 | 9/450 |
|
|
[考研] 343求調劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 0856求調劑 +13 | zhn03 2026-03-25 | 14/700 |
|
|
[考研] 調劑考研 +3 | 王杰一 2026-03-29 | 3/150 |
|
|
[考研] 299求調劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 286求調劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|