| 5 | 1/1 | 返回列表 |
| 查看: 3207 | 回復(fù): 8 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
wangjiang12鐵蟲 (小有名氣)
|
[求助]
VASP計算六方BN的缺陷時時進(jìn)行偶極修正, 即dipole corrections的問題 已有1人參與
|
|
|
各位VASP的高手 大家好!小弟是VASP的新手,最近在計算六方BN的缺陷,總共196個原子,計算時覺得應(yīng)該要加上偶極修正,即dipole corrections。所以在INCAR里面嘗試分別設(shè)置 IDIPOL = 1,2,3,4進(jìn)行計算。(沒有設(shè)置DIPOL和LDIPOL,因為不是非常懂) 計算結(jié)果如下: 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的計算結(jié)果比其他的大很多,不知道該采用哪個,請高手指導(dǎo)我如何正確進(jìn)行偶極修正的計算吧。 非常感謝 |
鐵蟲 (初入文壇)
木蟲 (正式寫手)
|
原則上都不能用,非立方結(jié)構(gòu)dipole-dipole作用可以達(dá)到monopole-monopole 的30%以上,quardpole-dipole作用在半導(dǎo)體有缺陷的情況下也不可忽略 最重要的是,這里得出的數(shù)都是認(rèn)為介電常數(shù)=1的,手冊上強(qiáng)調(diào)了好幾次,要手動除掉缺陷之間(超胞超胞之間)介質(zhì)的介電常數(shù),然后手動加到能量上 所以你可以 1.計算介電常數(shù),或者查到 2.測試dipole的位置,因為邊界條件的關(guān)系任何方法都不能直接算出體系離子貢獻(xiàn)的電矩(你假定這個東西是無限大的)。dipole取正電中心,到負(fù)電中心之間的各個點,直到校正項最小。也可能直接落在缺陷位置。此外要確保dipole落在電荷FFT的節(jié)點上,這個通過NGFX/NGFY/NGFZ控制,或移動缺陷到原點 3.測試超胞尺寸對偶極校正能量的影響,一直加大到能量變化不大(10%以內(nèi)),取可以接受的最大的超胞 4.把計算結(jié)果手動除以介電常數(shù),加到?jīng)]有修正的能量上。 ================ 假如你對能量不感興趣,或者你的電子結(jié)構(gòu)會因為這個內(nèi)部電矩有質(zhì)的變化(比如鐵磁和順磁的區(qū)別) 那這些都不用做了,因為vasp永遠(yuǎn)得不到正確的解 |
鐵蟲 (小有名氣)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +6 | 爭取九點睡 2026-03-28 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 340求調(diào)劑 +5 | jhx777 2026-03-27 | 5/250 |
|
|
[考研] 0703化學(xué)求調(diào)劑,各位老師看看我。! +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研]
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
[考研] 281求調(diào)劑 +6 | Koxui 2026-03-24 | 7/350 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 0854人工智能方向招收調(diào)劑 +4 | 章小魚567 2026-03-24 | 4/200 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|