 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2268 | 回復(fù): 4 | ||||||
cenwanglai榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
|
[交流]
【交流】從一個實(shí)際例子解讀IBZKPT 已有4人參與
|
|
IBZKPT可以用kpoints程序生成.如果KPOINTS中k-mesh是自動生成的,IBZKPT也會在vasp運(yùn)行過程中生成.有幾個概念不懂,求助: (1)k-point weights (2)tetrahedron method (3)tetrahedron connection tables tetrahedron與ISMEAR=-4/-5有關(guān)嗎? tetrahedron connection tables與DOS計(jì)算和Bandstructrue計(jì)算有關(guān)嗎?如果有是怎樣相關(guān)的? 相關(guān)問題vasp manual上的說法見 http://www.gaoyang168.com/bbs/viewthread.php?tid=2114071&fpage=0&view=&highlight=&page=4 下面是一個實(shí)際自動生成的IBZKPT: Automatically generated mesh 7 Reciprocal lattice 0.00000000000000 0.00000000000000 0.25000000000000 不可約k點(diǎn)坐標(biāo) 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點(diǎn)連成的不可約四面體簡并度 1 1 2 5 四個不可約k點(diǎn)做為頂點(diǎn)組成四面體 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點(diǎn),組成31個不可約四面體,這些四面體簡并度之和為300.每個四面體體積0.00333333,總體積為300*0.00333333=1 結(jié)合樓下的解釋,看了看IBZKPT的數(shù)字特征,猜想應(yīng)該這樣解釋(藍(lán)色). 不可約k點(diǎn)的產(chǎn)生和權(quán)重(重復(fù)數(shù)/簡并度)由POSCAR和KPOINTS兩個文件來決定.點(diǎn)群學(xué)的好的話,可以自己來手寫.也可以借助程序來生成.MS,XcrysDen都可以產(chǎn)生這樣的k點(diǎn),而且XcrysDen似乎可以圖形畫顯示倒空間的四面體劃分.不知道顯示的是不是不可約布里淵區(qū)。 布里淵區(qū)似乎是可變的，超胞的布里淵區(qū)似乎可以分解為不可約布里淵區(qū)。不理解，呵呵。 [ Last edited by cenwanglai on 2010-7-1 at 13:14 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 VASP | VASP相關(guān) | 第一性原理和電化學(xué) |
» 猜你喜歡
306求調(diào)劑
已經(jīng)有8人回復(fù)
-
招08考數(shù)學(xué)
已經(jīng)有12人回復(fù)
-
接收2026碩士調(diào)劑(學(xué)碩+專碩)
已經(jīng)有4人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有9人回復(fù)
-
306求0703調(diào)劑一志愿華中師范
已經(jīng)有6人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有8人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有3人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有3人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有3人回復(fù)
-
336化工調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- vasp對布里淵區(qū)積分問題
已經(jīng)有10人回復(fù)
- 【交流】有關(guān)布里淵區(qū)、K空間的問題
已經(jīng)有11人回復(fù)
- 【求助完畢】IBZKPT,KPOINTS與DOS計(jì)算
已經(jīng)有4人回復(fù)
燈塔守望者
鐵桿木蟲 (正式寫手)
- 1ST強(qiáng)帖: 2
- 應(yīng)助: 28 (小學(xué)生)
- 貴賓: 0.015
- 金幣: 5815
- 散金: 51
- 紅花: 40
- 帖子: 511
- 在線: 943.5小時
- 蟲號: 577747
- 注冊: 2008-06-27
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★
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點(diǎn)分布。 為了減少計(jì)算量，量化軟件會根據(jù)晶體的對稱性，把對第一布里淵區(qū)的積分簡化成對不可約布里淵區(qū)中不等價K點(diǎn)的積分。不可約布里淵區(qū)可通過對稱性操作還原第一布里淵區(qū)的所有信息。類似于超胞和元胞的關(guān)系。在第一布里淵區(qū)，如果兩個K點(diǎn)在同一類對稱操作中等價，那我們就認(rèn)為權(quán)重為2，在不可約布里淵區(qū)計(jì)算中只計(jì)算1次，然后乘2就可以了。以此類推。四面體積分類似于有限元方法，具體原理我也不清楚。 語言表述不夠準(zhǔn)確，個人看法，歡迎指正。 |
2樓2010-07-01 08:10:10
gavinliu7390
木蟲 (著名寫手)
葉落鷹飛
- 1ST強(qiáng)帖: 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):謝謝指點(diǎn) 2010-07-01 17:08:57
cenwanglai(金幣+3):呵呵,謝謝回復(fù)! 2010-07-01 15:00:59
zzy870720z(金幣+3):謝謝指點(diǎn) 2010-07-01 17:08:57
|
根據(jù)能帶理論， E是k的連續(xù)函數(shù)。 在解ks方程的時候，需要求band energy。 才能進(jìn)行自洽計(jì)算。 band energy就是布呂淵區(qū)對所有占據(jù)態(tài)k點(diǎn)的積分。 所以就需要對布呂淵區(qū)進(jìn)行撒點(diǎn)。例如，給10*10*10的k點(diǎn)就1000個k點(diǎn)。在解ks方程的時候，是應(yīng)該對這1000個k點(diǎn)都求出它的能量本征值。 但是由于晶體對稱性的關(guān)系，有些k點(diǎn)是對稱的，能量本征值是一樣的， 這就大大簡化了計(jì)算量。 之所以要?dú)w一，是因?yàn)槟芰繉點(diǎn)的積分可以化成對離散k點(diǎn)的求和。（找找相關(guān)文獻(xiàn)，應(yīng)該有公式。） ismear=-4或者-5 肯定是有區(qū)別的。 區(qū)別在于-5是布洛赫對k點(diǎn)的取法進(jìn)行了修正，防止在計(jì)算中有錯誤產(chǎn)生。具體的還是看文獻(xiàn)吧。 |
3樓2010-07-01 14:16:08
cenwanglai
榮譽(yù)版主 (知名作家)
老和山猥瑣派九段
- 1ST強(qiáng)帖: 5
- 應(yīng)助: 46 (小學(xué)生)
- 貴賓: 8.842
- 金幣: 7440.4
- 散金: 9654
- 紅花: 118
- 沙發(fā): 23
- 帖子: 5306
- 在線: 1961.4小時
- 蟲號: 537452
- 注冊: 2008-04-01
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
- 管轄: 仿真模擬
4樓2010-07-01 20:40:57
對馬拉二胡
至尊木蟲 (職業(yè)作家)
- 應(yīng)助: 27 (小學(xué)生)
- 金幣: 25940.6
- 散金: 244
- 紅花: 13
- 帖子: 4070
- 在線: 938.7小時
- 蟲號: 2471734
- 注冊: 2013-05-19
- 性別: GG
- 專業(yè): 無機(jī)非金屬材料
5樓2015-05-27 14:58:02
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 cenwanglai 的主題更新
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +5 | vv迷 2026-03-22 | 5/250 |
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | 十三加油 2026-03-21 | 6/300 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[基金申請] 山東省面上項(xiàng)目限額評審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 14/700 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +9 | 何潤采123 2026-03-18 | 11/550 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 329求調(diào)劑 +9 | 想上學(xué)吖吖 2026-03-19 | 9/450 |
|
|
|
[考研] 材料學(xué)求調(diào)劑 +4 | Stella_Yao 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|