| 5 | 1/1 | 返回列表 |
| 查看: 3433 | 回復(fù): 8 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
雜化泛函計(jì)算能帶的k點(diǎn)選擇問(wèn)題 已有4人參與
|
|||
|
使用雜化泛函計(jì)算能帶的時(shí)候,權(quán)重為零的k點(diǎn)是怎么生成的呢?我之前用dft計(jì)算的時(shí)候都是syml和gk.f生成k點(diǎn)但是計(jì)算產(chǎn)生的能帶感覺(jué)都不對(duì)啊,一點(diǎn)都不漂亮而且圖形都不對(duì)。用line-mode可以算出正確的圖但是雜化泛函計(jì)算的時(shí)候要把所有的點(diǎn)都寫出來(lái),總不可能一個(gè)一個(gè)點(diǎn)自己算出來(lái)寫上吧。有沒(méi)有什么方法可以查看line-mode模式產(chǎn)生的詳細(xì)k點(diǎn)呢? 發(fā)自小木蟲(chóng)Android客戶端 |
vasp 第一性原理 |
金蟲(chóng) (正式寫手)

金蟲(chóng) (正式寫手)
|
本帖內(nèi)容被屏蔽 |
金蟲(chóng) (正式寫手)
|
不好意思。樓上有一個(gè)回復(fù)是回復(fù)錯(cuò)帖子了。 你的這個(gè)情況,詳細(xì)可以看帖子: http://www.gaoyang168.com/bbs/viewthread.php?tid=6105659&fpage=1&target=blank 利用HSE方法計(jì)算能帶 復(fù)制過(guò)來(lái): “ 主要是k點(diǎn)的設(shè)置,由于HSE方法無(wú)法從電荷密度進(jìn)行非自恰,所以只能用自恰的方法計(jì)算能帶。 通常是先把IBZKPT 拷貝成KPOINTS, Automatically generated mesh 75 Reciprocal lattice 0.08333333333333 0.08333333333333 0.08333333333333 96 0.25000000000000 0.08333333333333 0.08333333333333 288 0.41666666666667 0.08333333333333 0.08333333333333 288 -0.41666666666667 0.08333333333333 0.08333333333333 288 -0.25000000000000 0.08333333333333 0.08333333333333 288 -0.08333333333333 0.08333333333333 0.08333333333333 288 0.25000000000000 0.25000000000000 0.08333333333333 288 0.41666666666667 0.25000000000000 0.08333333333333 576 -0.41666666666667 0.25000000000000 0.08333333333333 576 -0.25000000000000 0.25000000000000 0.08333333333333 576 -0.08333333333333 0.25000000000000 0.08333333333333 576 0.41666666666667 0.41666666666667 0.08333333333333 288 -0.41666666666667 0.41666666666667 0.08333333333333 576 -0.25000000000000 0.41666666666667 0.08333333333333 576 -0.08333333333333 0.41666666666667 0.08333333333333 576 -0.41666666666667 -0.41666666666667 0.08333333333333 288 -0.25000000000000 -0.41666666666667 0.08333333333333 576 -0.25000000000000 -0.25000000000000 0.08333333333333 288 0.25000000000000 0.25000000000000 0.25000000000000 96 0.41666666666667 0.25000000000000 0.25000000000000 288 -0.41666666666667 0.25000000000000 0.25000000000000 288 -0.25000000000000 0.25000000000000 0.25000000000000 288 0.41666666666667 0.41666666666667 0.25000000000000 288 -0.41666666666667 0.41666666666667 0.25000000000000 576 -0.25000000000000 0.41666666666667 0.25000000000000 576 -0.41666666666667 -0.41666666666667 0.25000000000000 288 0.41666666666667 0.41666666666667 0.41666666666667 96 -0.41666666666667 0.41666666666667 0.41666666666667 288 0.50000000 0.25000000 0.75000000 0.0 0.50000000 0.27777778 0.72222222 0.0 0.50000000 0.30555556 0.69444444 0.0 0.50000000 0.33333333 0.66666667 0.0 0.50000000 0.36111111 0.63888889 0.0 0.50000000 0.38888889 0.61111111 0.0 0.50000000 0.41666667 0.58333333 0.0 0.50000000 0.44444444 0.55555556 0.0 0.50000000 0.47222222 0.52777778 0.0 0.50000000 0.50000000 0.50000000 0.0 0.44444444 0.44444444 0.44444444 0.0 0.38888889 0.38888889 0.38888889 0.0 0.33333333 0.33333333 0.33333333 0.0 0.27777778 0.27777778 0.27777778 0.0 0.22222222 0.22222222 0.22222222 0.0 0.16666667 0.16666667 0.16666667 0.0 0.11111111 0.11111111 0.11111111 0.0 0.05555556 0.05555556 0.05555556 0.0 0.00000000 0.00000000 0.00000000 0.0 0.00000000 0.00000000 0.00000000 0.0 0.05555556 0.00000000 0.05555556 0.0 0.11111111 0.00000000 0.11111111 0.0 0.16666667 0.00000000 0.16666667 0.0 0.22222222 0.00000000 0.22222222 0.0 0.27777778 0.00000000 0.27777778 0.0 0.33333333 0.00000000 0.33333333 0.0 0.38888889 0.00000000 0.38888889 0.0 0.44444444 0.00000000 0.44444444 0.0 0.50000000 0.00000000 0.50000000 0.0 0.50000000 0.02777778 0.52777778 0.0 0.50000000 0.05555556 0.55555556 0.0 0.50000000 0.08333333 0.58333333 0.0 0.50000000 0.11111111 0.61111111 0.0 0.50000000 0.13888889 0.63888889 0.0 0.50000000 0.16666667 0.66666667 0.0 0.50000000 0.19444444 0.69444444 0.0 0.50000000 0.22222222 0.72222222 0.0 0.50000000 0.25000000 0.75000000 0.0 0.48611111 0.26388889 0.75000000 0.0 0.47222222 0.27777778 0.75000000 0.0 0.45833333 0.29166667 0.75000000 0.0 0.44444444 0.30555556 0.75000000 0.0 0.43055556 0.31944444 0.75000000 0.0 0.41666667 0.33333333 0.75000000 0.0 0.40277778 0.34722222 0.75000000 0.0 0.38888889 0.36111111 0.75000000 0.0 0.37500000 0.37500000 0.75000000 0.0 這里第二行75,意味著一共有75個(gè)k點(diǎn),接下來(lái)會(huì)有這75個(gè)k點(diǎn),后面跟著的是權(quán)重(如果大家不理解什么是權(quán)重,可以聯(lián)系我)。這些權(quán)重為0的,即為想計(jì)算的能帶k點(diǎn)。大家仔細(xì)觀察,就能體會(huì)出它是能帶的k點(diǎn)了,因?yàn)閗點(diǎn)是連續(xù)的。這里權(quán)重是0,意味著在自恰計(jì)算的時(shí)候,不會(huì)計(jì)算這些k點(diǎn),但是計(jì)算完后,會(huì)輸出這些本征值。達(dá)到計(jì)算能帶的目的。之后,我的方法是在EIGENVAL里刪除自恰的k點(diǎn)信息,之后就可以用正常的能帶處理程序處理了! 而你的問(wèn)題: 權(quán)重為零的k點(diǎn)是怎么生成的呢?我之前用dft計(jì)算的時(shí)候都是syml和gk.f生成k點(diǎn)但是計(jì)算產(chǎn)生的能帶感覺(jué)都不對(duì)啊,一點(diǎn)都不漂亮而且圖形都不對(duì)。用line-mode可以算出正確的圖但是雜化泛函計(jì)算的時(shí)候要把所有的點(diǎn)都寫出來(lái),總不可能一個(gè)一個(gè)點(diǎn)自己算出來(lái)寫上吧。有沒(méi)有什么方法可以查看line-mode模式產(chǎn)生的詳細(xì)k點(diǎn)呢? >>>>>我都是用syml gk.x來(lái)做能帶計(jì)算的產(chǎn)生K點(diǎn)工作,我覺(jué)得還行,不清楚你說(shuō)的不漂亮和圖形不對(duì)的具體意思。沒(méi)有用過(guò)line-mode產(chǎn)生k點(diǎn)過(guò) |

|
使用雜化泛函計(jì)算能帶的時(shí)候,權(quán)重為零的k點(diǎn)是怎么生成的呢? 靜態(tài)計(jì)算IBZKBT里復(fù)制過(guò)去 我之前用dft計(jì)算的時(shí)候都是syml和gk.f生成k點(diǎn)但是計(jì)算產(chǎn)生的能帶感覺(jué)都不對(duì)啊,一點(diǎn)都不漂亮而且圖形都不對(duì)。用line-mode可以算出正確的圖但是雜化泛函計(jì)算的時(shí)候要把所有的點(diǎn)都寫出來(lái),總不可能一個(gè)一個(gè)點(diǎn)自己算出來(lái)寫上吧。有沒(méi)有什么方法可以查看line-mode模式產(chǎn)生的詳細(xì)k點(diǎn)呢? HSE能帶計(jì)算的KPOINTS是拼接的,前面是靜態(tài)計(jì)算IBZKBT里復(fù)制過(guò)去,后面是PBE能帶計(jì)算的outcar里的K點(diǎn)復(fù)制過(guò)去,希望我沒(méi)記錯(cuò) |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿華理,數(shù)一英一285求A區(qū)調(diào)劑 +8 | AZMK 2026-03-25 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動(dòng)285求調(diào)劑 +3 | AZMK 2026-03-28 | 6/300 |
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +6 | 1428151015 2026-03-27 | 6/300 |
|
|
[考研] 一志愿北京化工大學(xué) 070300 學(xué)碩 336分 求調(diào)劑 +9 | vv迷 2026-03-22 | 9/450 |
|
|
[考研] 復(fù)試調(diào)劑,一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研]
|
鐘llll 2026-03-26 | 4/200 |
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 340求調(diào)劑 +3 | Amber00 2026-03-26 | 3/150 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研]
|
酥酥魚(yú).. 2026-03-21 | 4/200 |
|