| 查看: 1952 | 回復: 10 | ||||
[求助]
Silicane靜態(tài)計算之后利用Bader分析發(fā)現(xiàn)兩個Si的電荷數(shù)居然不等,還相差很大…… 已有1人參與
|
||||
|
前幾天計算Silicane(單層Silicene兩邊用H飽和的結構,Silicene是一種二維材料),發(fā)現(xiàn)一個問題…… 結構優(yōu)化,靜態(tài)計算之後,進行Bader分析,結果相當出人意料,兩個Si的電荷數(shù)居然不等! 要知道,這兩個Si是完全等價的…… 分析結果如下: # X Y Z CHARGE MIN DIST ATOMIC VOL -------------------------------------------------------------------------------- 1 2.2286 0.0000 9.5609 3.5469 0.7023 20.9265 2 4.4571 0.0000 10.2849 3.2780 0.7023 19.4710 3 2.2286 0.0000 8.0619 1.5876 0.7280 107.8755 4 4.4571 0.0000 11.7838 1.5876 0.7280 107.8053 -------------------------------------------------------------------------------- VACUUM CHARGE: 0.0000 VACUUM VOLUME: 0.0000 NUMBER OF ELECTRONS: 10.0000 一個Si是3.5469,另一個是3.2780,這兩個數(shù)值相差的還挺大 開始以為是精度不夠,但是加大K點到45X45X1,還是不等…… 嘗試加大NGX, NGY, NGZ 還是不行 將原胞擴展成2X2超胞之后,也還是不行,Si原子之間的電子數(shù)就是不等…… 不知道有沒有前輩遇到過這個問題,還請各位多多賜教…… |

木蟲 (著名寫手)
|
這也正常,因為Bader以電荷密度極小處做為原子邊界分割點,這種方法其實是不適合共價體系。Si-Si鍵當中有較大電荷密度,如以下文章圖1(b): http://journals.aps.org/prb/abstract/10.1103/PhysRevB.91.205433 所以Bader方法在分割原子邊界時,可定會避開Si-Si中間電荷極大點,而選擇靠近其中一個原子的某個地方做為邊界,所以會導致兩個原子不均等了。 |

木蟲 (著名寫手)
|
應該很多人已經在VASP上實現(xiàn)了瑪麗肯,但我不知道是否有人公開他們的代碼,樓主可以搜索一下看看。 巴德也是UT Ausdin那幫人做的,只是他們公開了他們的代碼而已。個人淺見,對于量子化學鄰域,最常用的還是原子軌道投影的一些方法,只是不同方法所用的原子軌道和處理細節(jié)不一樣罷了,比如:Mulliken和Loewdin population analyses是最有名的兩種。用原子軌道感覺比較物理,無論得到的結果是否符合預期,都有跡可循。用Mulliken和Loewdin方法得到的軌道投影,還可以進一步用來計算bond order,這個可以用來指示材料的共價性。比如如下文章里的公式(2)和其中參考文獻[68-73]: http://pubs.acs.org/doi/abs/10.1021/acs.jpcc.5b04669 PS: 謝謝樓主的慷慨金幣!可以多下半本書了。 |


送紅花一朵
|

金蟲 (小有名氣)
|
呵呵,我計算了硅稀與硅烷(silicene,silicane)異質結,即將硅稀層放在硅烷層的上面。后來我也用了那個bader charge 分析,結果在最后一步bader CHGRG -ref CHGCAR_sum卻出現(xiàn)了錯誤,should be no new maxima in edge refinement。不知道樓上大神遇見過這種情況沒有,拜求解惑 發(fā)自小木蟲Android客戶端 |
專家顧問 (著名寫手)
Ab Initio Amateur
![]() |
專家經驗: +224 |
|
http://schmidt.chem.wisc.edu/nbosoftware可以將平面波投影到AO上,該組提供了程序,然后另一個程序可以進行NBO。Mulliken可以在投影之后做,目測得自己寫程序,理論上Mulliken少了正交化,還是Hirshfeld靠得住一些。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 318求調劑 +4 | plum李子 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 306求調劑 +4 | 來好運來來來 2026-03-22 | 4/200 |
|
|
[考研] 生物學071000 329分求調劑 +4 | 我愛生物生物愛?/a> 2026-03-17 | 4/200 |
|
|
[考研] 298求調劑 +4 | 上岸6666@ 2026-03-20 | 4/200 |
|
|
[考研] 0805材料320求調劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
[考研] 求調劑 +3 | 白QF 2026-03-21 | 3/150 |
|
|
[考研] 316求調劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 299求調劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 307求調劑 +10 | 冷笙123 2026-03-17 | 10/500 |
|
|
[考研] 一志愿武理材料305分求調劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿蘇州大學材料求調劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 329求調劑 +9 | 想上學吖吖 2026-03-19 | 9/450 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 廣西大學家禽遺傳育種課題組2026年碩士招生(接收計算機專業(yè)調劑) +3 | 123阿標 2026-03-17 | 3/150 |
|
|
[考研] 0856調劑,是學校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 081700化工學碩調劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 085601材料工程專碩求調劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 085600材料與化工求調劑 +6 | 緒幸與子 2026-03-17 | 6/300 |
|
|
[考研] 材料考研調劑 +3 | xwt。 2026-03-19 | 3/150 |
|