 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 9 | 1/1 | 返回列表 |
| 查看: 3178 | 回復(fù): 8 | |||
wangjiang12鐵蟲 (小有名氣)
|
[求助]
VASP計算六方BN的缺陷時時進行偶極修正， 即dipole corrections的問題 已有1人參與
|
|
各位VASP的高手 大家好！小弟是VASP的新手，最近在計算六方BN的缺陷，總共196個原子，計算時覺得應(yīng)該要加上偶極修正，即dipole corrections。所以在INCAR里面嘗試分別設(shè)置 IDIPOL = 1,2,3,4進行計算。（沒有設(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)我如何正確進行偶極修正的計算吧。 非常感謝 |
回復(fù)此樓» 猜你喜歡
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
293求調(diào)劑
已經(jīng)有13人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 【求助完畢】硫酸根離子的結(jié)構(gòu)優(yōu)化，偶極修正及勢能修正。
已經(jīng)有20人回復(fù)
golddoushi
木蟲 (正式寫手)
- 應(yīng)助: 36 (小學(xué)生)
- 金幣: 4007.2
- 散金: 25
- 紅花: 15
- 帖子: 337
- 在線: 447.7小時
- 蟲號: 2546370
- 注冊: 2013-07-15
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wangjiang12: 金幣+5, ★★★★★最佳答案 2014-03-11 18:31:05
感謝參與，應(yīng)助指數(shù) +1
wangjiang12: 金幣+5, ★★★★★最佳答案 2014-03-11 18:31:05
|
原則上都不能用，非立方結(jié)構(gòu)dipole-dipole作用可以達到monopole-monopole 的30%以上，quardpole-dipole作用在半導(dǎo)體有缺陷的情況下也不可忽略 最重要的是，這里得出的數(shù)都是認為介電常數(shù)=1的，手冊上強調(diào)了好幾次，要手動除掉缺陷之間（超胞超胞之間）介質(zhì)的介電常數(shù)，然后手動加到能量上 所以你可以 1.計算介電常數(shù)，或者查到 2.測試dipole的位置，因為邊界條件的關(guān)系任何方法都不能直接算出體系離子貢獻的電矩（你假定這個東西是無限大的）。dipole取正電中心，到負電中心之間的各個點，直到校正項最小。也可能直接落在缺陷位置。此外要確保dipole落在電荷FFT的節(jié)點上，這個通過NGFX/NGFY/NGFZ控制，或移動缺陷到原點 3.測試超胞尺寸對偶極校正能量的影響，一直加大到能量變化不大（10%以內(nèi)），取可以接受的最大的超胞 4.把計算結(jié)果手動除以介電常數(shù)，加到?jīng)]有修正的能量上。 ================ 假如你對能量不感興趣，或者你的電子結(jié)構(gòu)會因為這個內(nèi)部電矩有質(zhì)的變化（比如鐵磁和順磁的區(qū)別） 那這些都不用做了，因為vasp永遠得不到正確的解 |
2樓2014-03-11 13:05:05
wangjiang12
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5
- 散金: 5
- 帖子: 104
- 在線: 20.8小時
- 蟲號: 2317103
- 注冊: 2013-03-04
- 性別: GG
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
3樓2014-03-11 15:17:00
golddoushi
木蟲 (正式寫手)
- 應(yīng)助: 36 (小學(xué)生)
- 金幣: 4007.2
- 散金: 25
- 紅花: 15
- 帖子: 337
- 在線: 447.7小時
- 蟲號: 2546370
- 注冊: 2013-07-15
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2014-03-11 15:57:47
wangjiang12
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5
- 散金: 5
- 帖子: 104
- 在線: 20.8小時
- 蟲號: 2317103
- 注冊: 2013-03-04
- 性別: GG
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
5樓2014-03-11 16:28:57
golddoushi
木蟲 (正式寫手)
- 應(yīng)助: 36 (小學(xué)生)
- 金幣: 4007.2
- 散金: 25
- 紅花: 15
- 帖子: 337
- 在線: 447.7小時
- 蟲號: 2546370
- 注冊: 2013-07-15
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
6樓2014-03-11 18:34:51
wangjiang12
鐵蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5
- 散金: 5
- 帖子: 104
- 在線: 20.8小時
- 蟲號: 2317103
- 注冊: 2013-03-04
- 性別: GG
- 專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
7樓2014-03-12 10:11:11
golddoushi
木蟲 (正式寫手)
- 應(yīng)助: 36 (小學(xué)生)
- 金幣: 4007.2
- 散金: 25
- 紅花: 15
- 帖子: 337
- 在線: 447.7小時
- 蟲號: 2546370
- 注冊: 2013-07-15
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
8樓2014-03-12 10:59:37
jiangge122
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 49
- 帖子: 32
- 在線: 1.9小時
- 蟲號: 3055054
- 注冊: 2014-03-15
- 性別: GG
- 專業(yè): 半導(dǎo)體材料
9樓2014-03-19 13:27:21
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 wangjiang12 的主題更新
| 9 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +8 | zjl的號 2026-03-16 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +9 | 負心者當(dāng)誅 2026-03-11 | 9/450 |
|
|
|
[考研] 化學(xué)工程321分求調(diào)劑 +11 | 大米飯！ 2026-03-15 | 14/700 |
|
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 271求調(diào)劑 +12 | 生如夏花… 2026-03-11 | 14/700 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家，不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 336求調(diào)劑 +6 | Iuruoh 2026-03-11 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +7 | 7712b 2026-03-13 | 7/350 |
|
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|