版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

wangjiang12

鐵蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 5
散金: 5
帖子: 104
在線: 20.8小時(shí)
蟲號(hào): 2317103
注冊(cè): 2013-03-04
性別: GG
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

[求助] VASP計(jì)算六方BN的缺陷時(shí)時(shí)進(jìn)行偶極修正，即dipole corrections的問題已有1人參與

各位VASP的高手
   大家好！小弟是VASP的新手，最近在計(jì)算六方BN的缺陷，總共196個(gè)原子，計(jì)算時(shí)覺得應(yīng)該要加上偶極修正，即dipole corrections。所以在INCAR里面嘗試分別設(shè)置 IDIPOL = 1,2,3,4進(jìn)行計(jì)算。（沒有設(shè)置DIPOL和LDIPOL，因?yàn)椴皇欠浅６?br />    計(jì)算結(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的計(jì)算結(jié)果比其他的大很多，不知道該采用哪個(gè)，請(qǐng)高手指導(dǎo)我如何正確進(jìn)行偶極修正的計(jì)算吧。

非常感謝

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士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博士申請(qǐng) 已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【求助完畢】硫酸根離子的結(jié)構(gòu)優(yōu)化，偶極修正及勢(shì)能修正。已經(jīng)有20人回復(fù)

1樓 2014-03-11 10:34:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

golddoushi

木蟲 (正式寫手)

應(yīng)助: 36 (小學(xué)生)
金幣: 4007.2
散金: 25
紅花: 15
帖子: 337
在線: 447.7小時(shí)
蟲號(hào): 2546370
注冊(cè): 2013-07-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
wangjiang12: 金幣+5, ★★★★★最佳答案 2014-03-11 18:31:05

原則上都不能用，非立方結(jié)構(gòu)dipole-dipole作用可以達(dá)到monopole-monopole 的30%以上，quardpole-dipole作用在半導(dǎo)體有缺陷的情況下也不可忽略
最重要的是，這里得出的數(shù)都是認(rèn)為介電常數(shù)=1的，手冊(cè)上強(qiáng)調(diào)了好幾次，要手動(dòng)除掉缺陷之間（超胞超胞之間）介質(zhì)的介電常數(shù)，然后手動(dòng)加到能量上
所以你可以
1.計(jì)算介電常數(shù)，或者查到
2.測(cè)試dipole的位置，因?yàn)檫吔鐥l件的關(guān)系任何方法都不能直接算出體系離子貢獻(xiàn)的電矩（你假定這個(gè)東西是無限大的）。dipole取正電中心，到負(fù)電中心之間的各個(gè)點(diǎn)，直到校正項(xiàng)最小。也可能直接落在缺陷位置。此外要確保dipole落在電荷FFT的節(jié)點(diǎn)上，這個(gè)通過NGFX/NGFY/NGFZ控制，或移動(dòng)缺陷到原點(diǎn)
3.測(cè)試超胞尺寸對(duì)偶極校正能量的影響，一直加大到能量變化不大（10%以內(nèi)），取可以接受的最大的超胞
4.把計(jì)算結(jié)果手動(dòng)除以介電常數(shù)，加到?jīng)]有修正的能量上。

================
假如你對(duì)能量不感興趣，或者你的電子結(jié)構(gòu)會(huì)因?yàn)檫@個(gè)內(nèi)部電矩有質(zhì)的變化（比如鐵磁和順磁的區(qū)別）
那這些都不用做了，因?yàn)関asp永遠(yuǎn)得不到正確的解

贊一下

回復(fù)此樓

2樓2014-03-11 13:05:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangjiang12

鐵蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 5
散金: 5
帖子: 104
在線: 20.8小時(shí)
蟲號(hào): 2317103
注冊(cè): 2013-03-04
性別: GG
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

引用回帖:

2樓: Originally posted by golddoushi at 2014-03-11 13:05:05
原則上都不能用，非立方結(jié)構(gòu)dipole-dipole作用可以達(dá)到monopole-monopole 的30%以上，quardpole-dipole作用在半導(dǎo)體有缺陷的情況下也不可忽略
最重要的是，這里得出的數(shù)都是認(rèn)為介電常數(shù)=1的，手冊(cè)上強(qiáng)調(diào)了好幾次， ...

非常感謝您的回復(fù)，看起來過程好復(fù)雜啊。那您認(rèn)為如果我不做dipole corrections, 直接對(duì)比不同缺陷排布下的能量大小，可以嗎

贊一下

回復(fù)此樓

3樓2014-03-11 15:17:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

golddoushi

木蟲 (正式寫手)

應(yīng)助: 36 (小學(xué)生)
金幣: 4007.2
散金: 25
紅花: 15
帖子: 337
在線: 447.7小時(shí)
蟲號(hào): 2546370
注冊(cè): 2013-07-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

3樓: Originally posted by wangjiang12 at 2014-03-11 15:17:00
非常感謝您的回復(fù)，看起來過程好復(fù)雜啊。那您認(rèn)為如果我不做dipole corrections, 直接對(duì)比不同缺陷排布下的能量大小，可以嗎...

通常情況約有 0.1eV ~ 1eV 數(shù)量級(jí)的誤差。如果你的材料介電常數(shù)較小，誤差還會(huì)更大。
dipole - dipole 作用的能量正比于 1/L3， monopole - monopole 正比于1/L，quardpole - dipole正比 1/L^5 或者更高次方。他們都反比于介電常數(shù)。
不過你不提這事也沒人會(huì)說你，這個(gè)比較偏

贊一下

回復(fù)此樓

4樓2014-03-11 15:57:47

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangjiang12

鐵蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 5
散金: 5
帖子: 104
在線: 20.8小時(shí)
蟲號(hào): 2317103
注冊(cè): 2013-03-04
性別: GG
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

引用回帖:

4樓: Originally posted by golddoushi at 2014-03-11 15:57:47
通常情況約有 0.1eV ~ 1eV 數(shù)量級(jí)的誤差。如果你的材料介電常數(shù)較小，誤差還會(huì)更大。
dipole - dipole 作用的能量正比于 1/L3， monopole - monopole 正比于1/L，quardpole - dipole正比 1/L^5 或者更高次方。他們 ...

哦，好的，我明白了，謝謝。對(duì)了，請(qǐng)問您一下，除了dipole corrections外，您覺得用VASP計(jì)算六方NB還有什么其他需要修正或者注意的地方嗎？

贊一下

回復(fù)此樓

5樓2014-03-11 16:28:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

golddoushi

木蟲 (正式寫手)

應(yīng)助: 36 (小學(xué)生)
金幣: 4007.2
散金: 25
紅花: 15
帖子: 337
在線: 447.7小時(shí)
蟲號(hào): 2546370
注冊(cè): 2013-07-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

5樓: Originally posted by wangjiang12 at 2014-03-11 16:28:57
哦，好的，我明白了，謝謝。對(duì)了，請(qǐng)問您一下，除了dipole corrections外，您覺得用VASP計(jì)算六方NB還有什么其他需要修正或者注意的地方嗎？...

具體不太清楚，我見過有人用半芯態(tài)的贗勢(shì)算的比較準(zhǔn)

贊一下

回復(fù)此樓

6樓2014-03-11 18:34:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wangjiang12

鐵蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 5
散金: 5
帖子: 104
在線: 20.8小時(shí)
蟲號(hào): 2317103
注冊(cè): 2013-03-04
性別: GG
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

引用回帖:

6樓: Originally posted by golddoushi at 2014-03-11 18:34:51
具體不太清楚，我見過有人用半芯態(tài)的贗勢(shì)算的比較準(zhǔn)...

半芯態(tài)的贗勢(shì)是指什么�。�

謝謝。

贊一下

回復(fù)此樓

7樓2014-03-12 10:11:11

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

golddoushi

木蟲 (正式寫手)

應(yīng)助: 36 (小學(xué)生)
金幣: 4007.2
散金: 25
紅花: 15
帖子: 337
在線: 447.7小時(shí)
蟲號(hào): 2546370
注冊(cè): 2013-07-15
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

7樓: Originally posted by wangjiang12 at 2014-03-12 10:11:11
半芯態(tài)的贗勢(shì)是指什么��？

謝謝。...

就是把N，和B內(nèi)層的s電子當(dāng)不那么活潑的價(jià)電子處理。

贊一下

回復(fù)此樓

8樓2014-03-12 10:59:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

jiangge122

鐵蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 49
帖子: 32
在線: 1.9小時(shí)
蟲號(hào): 3055054
注冊(cè): 2014-03-15
性別: GG
專業(yè): 半導(dǎo)體材料

引用回帖:

您好，我對(duì)這個(gè)dipole correction也非常感興趣。那個(gè)計(jì)算過程有點(diǎn)復(fù)雜，請(qǐng)問有例子可以參考一下嗎？還有，如果我用DFTB+計(jì)算，是否還需要考慮這個(gè)因素呢？謝謝

贊一下

回復(fù)此樓

9樓2014-03-19 13:27:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wangjiang12 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考博] 環(huán)境領(lǐng)域全國重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3	QGZDSYS 2026-03-13	3/150	2026-03-18 00:13 by PY-Hzb
[考研] 268求調(diào)劑 +3	簡(jiǎn)單點(diǎn)0 2026-03-17	3/150	2026-03-17 23:15 by wangjihu
[考研] 341求調(diào)劑 +5	搗蛋豬豬 2026-03-11	7/350	2026-03-17 19:09 by 搗蛋豬豬
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6	大火山小火山 2026-03-16	8/400	2026-03-17 15:05 by 無懈可擊111
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 211本，11408一志愿中科院277分，曾在中科院自動(dòng)化所實(shí)習(xí) +6	Losir 2026-03-12	7/350	2026-03-17 12:09 by danranxie
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3	增銳漏人 2026-03-17	4/200	2026-03-17 09:26 by xs74101122
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[教師之家] 焦慮 +7	水冰月月野兔 2026-03-13	9/450	2026-03-16 10:00 by Quakerbird
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6	zlingli 2026-03-13	6/300	2026-03-15 20:00 by ryzcf
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 工科，求調(diào)劑 +3	我887 2026-03-11	3/150	2026-03-13 21:39 by JourneyLucky
[考研] 315求調(diào)劑 +9	小羊小羊_ 2026-03-11	10/500	2026-03-13 21:13 by SXNU李老師
[考研] 281求調(diào)劑 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么
[考研] 310求調(diào)劑 +3	【上上簽】 2026-03-11	3/150	2026-03-13 16:16 by JourneyLucky
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8	努力奮斗112 2026-03-12	9/450	2026-03-13 15:41 by JourneyLucky
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

24小時(shí)熱門版塊排行榜

[求助] VASP計(jì)算六方BN的缺陷時(shí)時(shí)進(jìn)行偶極修正， 即dipole corrections的問題 已有1人參與

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

【答案】應(yīng)助回帖

[求助] VASP計(jì)算六方BN的缺陷時(shí)時(shí)進(jìn)行偶極修正，即dipole corrections的問題已有1人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: