 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2270 | 回復(fù): 4 | ||||||
cenwanglai榮譽版主 (知名作家)
老和山猥瑣派九段
|
[交流]
【交流】從一個實際例子解讀IBZKPT 已有4人參與
|
|
IBZKPT可以用kpoints程序生成.如果KPOINTS中k-mesh是自動生成的,IBZKPT也會在vasp運行過程中生成.有幾個概念不懂,求助: (1)k-point weights (2)tetrahedron method (3)tetrahedron connection tables tetrahedron與ISMEAR=-4/-5有關(guān)嗎? tetrahedron connection tables與DOS計算和Bandstructrue計算有關(guān)嗎?如果有是怎樣相關(guān)的? 相關(guān)問題vasp manual上的說法見 http://www.gaoyang168.com/bbs/viewthread.php?tid=2114071&fpage=0&view=&highlight=&page=4 下面是一個實際自動生成的IBZKPT: Automatically generated mesh 7 Reciprocal lattice 0.00000000000000 0.00000000000000 0.25000000000000 不可約k點坐標 2 簡并度(重復(fù)個數(shù)) 0.20000000000000 0.00000000000000 0.25000000000000 6 0.40000000000000 0.00000000000000 0.25000000000000 6 -0.40000000000000 0.00000000000000 0.25000000000000 6 -0.20000000000000 0.00000000000000 0.25000000000000 6 0.40000000000000 0.20000000000000 0.25000000000000 12 -0.40000000000000 0.20000000000000 0.25000000000000 12 !綠色數(shù)字表示該行kpoint的weight,設(shè)置相對大小后會自動歸一化.它是什么物理意義?根據(jù)什么來定義比例大小?有什么用? Tetrahedra 31 不可約四面體個數(shù) 0.00333333333333 不可約四面體體積 8 后面四個k點連成的不可約四面體簡并度 1 1 2 5 四個不可約k點做為頂點組成四面體 12 1 2 5 5 12 2 3 5 6 12 2 5 5 6 12 3 4 5 6 8 4 5 6 6 8 2 5 6 6 12 3 4 4 7 24 3 4 6 7 12 4 6 6 7 8 3 4 7 7 12 4 6 7 7 12 3 6 7 7 8 2 3 6 6 12 3 6 6 7 12 1 2 2 5 12 2 2 5 6 2 2 2 6 6 2 1 1 2 2 12 2 4 5 6 4 2 4 6 6 12 2 3 4 6 12 6 6 7 7 24 6 7 7 7 12 7 7 7 7 12 3 3 4 7 2 3 3 7 7 2 5 5 6 6 2 1 1 5 5 4 3 5 6 6 2 4 4 7 7 第一列個數(shù)為31,和為300.意思是7個不可約k點,組成31個不可約四面體,這些四面體簡并度之和為300.每個四面體體積0.00333333,總體積為300*0.00333333=1 結(jié)合樓下的解釋,看了看IBZKPT的數(shù)字特征,猜想應(yīng)該這樣解釋(藍色). 不可約k點的產(chǎn)生和權(quán)重(重復(fù)數(shù)/簡并度)由POSCAR和KPOINTS兩個文件來決定.點群學(xué)的好的話,可以自己來手寫.也可以借助程序來生成.MS,XcrysDen都可以產(chǎn)生這樣的k點,而且XcrysDen似乎可以圖形畫顯示倒空間的四面體劃分.不知道顯示的是不是不可約布里淵區(qū)。 布里淵區(qū)似乎是可變的，超胞的布里淵區(qū)似乎可以分解為不可約布里淵區(qū)。不理解，呵呵。 [ Last edited by cenwanglai on 2010-7-1 at 13:14 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 VASP | VASP相關(guān) | 第一性原理和電化學(xué) |
» 猜你喜歡
0856材料專碩353求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有6人回復(fù)
-
293求調(diào)劑
已經(jīng)有4人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有10人回復(fù)
-
【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】
已經(jīng)有7人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
材料與化工328分調(diào)劑
已經(jīng)有4人回復(fù)
-
279分求調(diào)劑 一志愿211
已經(jīng)有20人回復(fù)
-
一志愿華東理工大學(xué)081700，初試分數(shù)271
已經(jīng)有6人回復(fù)
-
085404電子信息284分求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- vasp對布里淵區(qū)積分問題
已經(jīng)有10人回復(fù)
- 【交流】有關(guān)布里淵區(qū)、K空間的問題
已經(jīng)有11人回復(fù)
- 【求助完畢】IBZKPT,KPOINTS與DOS計算
已經(jīng)有4人回復(fù)
★ ★
zzy870720z(金幣+2):謝謝指導(dǎo) 2010-07-01 08:13:00
cenwanglai(金幣+3):呵呵,不可約布里淵區(qū)，理解。謝謝。這個程序可以根據(jù)POSCAR自動生成。 2010-07-01 08:13:33
zzy870720z(金幣+2):謝謝指導(dǎo) 2010-07-01 08:13:00
cenwanglai(金幣+3):呵呵,不可約布里淵區(qū)，理解。謝謝。這個程序可以根據(jù)POSCAR自動生成。 2010-07-01 08:13:33
|
IBZKPT應(yīng)該是不可約布里淵區(qū)的的K點分布。 為了減少計算量，量化軟件會根據(jù)晶體的對稱性，把對第一布里淵區(qū)的積分簡化成對不可約布里淵區(qū)中不等價K點的積分。不可約布里淵區(qū)可通過對稱性操作還原第一布里淵區(qū)的所有信息。類似于超胞和元胞的關(guān)系。在第一布里淵區(qū)，如果兩個K點在同一類對稱操作中等價，那我們就認為權(quán)重為2，在不可約布里淵區(qū)計算中只計算1次，然后乘2就可以了。以此類推。四面體積分類似于有限元方法，具體原理我也不清楚。 語言表述不夠準確，個人看法，歡迎指正。 |
2樓2010-07-01 08:10:10
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強帖: 3
- 應(yīng)助: 51 (初中生)
- 貴賓: 0.481
- 金幣: 1564.7
- 散金: 2975
- 紅花: 34
- 帖子: 1278
- 在線: 241.4小時
- 蟲號: 675997
- 注冊: 2008-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
★ ★ ★
cenwanglai(金幣+3):呵呵,謝謝回復(fù)! 2010-07-01 15:00:59
zzy870720z(金幣+3):謝謝指點 2010-07-01 17:08:57
cenwanglai(金幣+3):呵呵,謝謝回復(fù)! 2010-07-01 15:00:59
zzy870720z(金幣+3):謝謝指點 2010-07-01 17:08:57
|
根據(jù)能帶理論， E是k的連續(xù)函數(shù)。 在解ks方程的時候，需要求band energy。 才能進行自洽計算。 band energy就是布呂淵區(qū)對所有占據(jù)態(tài)k點的積分。 所以就需要對布呂淵區(qū)進行撒點。例如，給10*10*10的k點就1000個k點。在解ks方程的時候，是應(yīng)該對這1000個k點都求出它的能量本征值。 但是由于晶體對稱性的關(guān)系，有些k點是對稱的，能量本征值是一樣的， 這就大大簡化了計算量。 之所以要歸一，是因為能量對k點的積分可以化成對離散k點的求和。（找找相關(guān)文獻，應(yīng)該有公式。） ismear=-4或者-5 肯定是有區(qū)別的。 區(qū)別在于-5是布洛赫對k點的取法進行了修正，防止在計算中有錯誤產(chǎn)生。具體的還是看文獻吧。 |
3樓2010-07-01 14:16:08
cenwanglai
榮譽版主 (知名作家)
老和山猥瑣派九段
- 1ST強帖: 5
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
- 管轄: 仿真模擬
4樓2010-07-01 20:40:57
對馬拉二胡
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 27 (小學(xué)生)
- 金幣: 25945.6
- 散金: 244
- 紅花: 13
- 帖子: 4071
- 在線: 938.7小時
- 蟲號: 2471734
- 注冊: 2013-05-19
- 性別: GG
- 專業(yè): 無機非金屬材料
5樓2015-05-27 14:58:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cenwanglai 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +6 | holy J 2026-03-21 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 081700 調(diào)劑 267分 +8 | 迷人的哈哈 2026-03-23 | 8/400 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+4
|
孅華 2026-03-22 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4 | 晨昏線與星海 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研]
320求調(diào)劑0856
+3
|
不想起名字112 2026-03-19 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|