| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1601 | 回復: 4 | ||
| 【懸賞金幣】回答本帖問題,作者uspex將贈送您 15 個金幣 | ||
[求助]
wannier90計算得到_hr.dat文件,如何得到onsite和hopping能量 已有1人參與
|
||
|
最近看到一些文獻,分析了結構中每個原子每個軌道的onsite和hopping的能量,然后結合晶體場理論做出了一些分析,我想知道是如何得到這些能量的。求大佬們指點一下。 發(fā)自小木蟲手機客戶端 |
木蟲 (正式寫手)
|
“每個原子每個軌道的 onsite 能量和 hopping 能量”,通常是通過從第一性原理計算得到的能帶結構,投影到局域軌道基底,構建緊束縛模型(tight-binding model)或 Wannier 函數(shù)哈密頓量得到的。 1. 先從 DFT 得到哈密頓量 你先用 DFT(VASP、QE、WIEN2k、ELK…)算出自洽電荷密度和能帶。 DFT 本身得到的是平面波基底下的哈密頓量,不是直接的“原子軌道能量和躍遷矩陣”。 2. 局域化:投影到原子軌道基底 有兩條主流方法: 方法 A:直接投影到原子軌道 (Projected Tight Binding) 對平面波波函數(shù)做投影,得到每個原子軌道(如 d_xy, d_xz, p_x)的投影系數(shù)。 在這個基底上構造哈密頓量矩陣 𝐻𝑖𝑗。對角元 𝐻𝑖𝑖 就是所謂 onsite energy。 非對角元𝐻𝑖𝑗 (i≠j)就是原子 i 和 j 之間的 hopping 參數(shù)(躍遷能)。 這種方法計算比較直接,但對軌道正交性要求比較高。 方法 B:Wannierization (Maximally Localized Wannier Functions, MLWF) 先選擇一組目標能帶(比如價帶+導帶的 d/p 帶)。 用 Wannier90 等工具生成局域 Wannier 函數(shù)。 在 Wannier 基底下的哈密頓量是稀疏的,矩陣元就是: 對角元 = onsite energy 非對角元 = hopping (可能包括長程跳躍) 好處:Wannier 化的哈密頓量可以很好地再現(xiàn)原始 DFT 能帶,而且 Wannier 函數(shù)具有類似原子軌道的局域性,可以很好地和晶體場理論結合。 |
木蟲 (正式寫手)
|
“每個原子每個軌道的 onsite 能量和 hopping 能量”,通常是通過從第一性原理計算得到的能帶結構,投影到局域軌道基底,構建緊束縛模型(tight-binding model)或 Wannier 函數(shù)哈密頓量得到的。 1. 先從 DFT 得到哈密頓量 你先用 DFT(VASP、QE、WIEN2k、ELK…)算出自洽電荷密度和能帶。 DFT 本身得到的是平面波基底下的哈密頓量,不是直接的“原子軌道能量和躍遷矩陣”。 2. 局域化:投影到原子軌道基底 有兩條主流方法: 方法 A:直接投影到原子軌道 (Projected Tight Binding) 對平面波波函數(shù)做投影,得到每個原子軌道(如 d_xy, d_xz, p_x)的投影系數(shù)。 在這個基底上構造哈密頓量矩陣H_ij;。對角元 H_ii; 就是所謂 onsite energy。 非對角元H_ij; (i≠j)就是原子 i 和 j 之間的 hopping 參數(shù)(躍遷能)。 這種方法計算比較直接,但對軌道正交性要求比較高。 方法 B:Wannierization (Maximally Localized Wannier Functions, MLWF) 先選擇一組目標能帶(比如價帶+導帶的 d/p 帶)。 用 Wannier90 等工具生成局域 Wannier 函數(shù)。 在 Wannier 基底下的哈密頓量是稀疏的,矩陣元就是: 對角元 = onsite energy 非對角元 = hopping (可能包括長程跳躍) 好處:Wannier 化的哈密頓量可以很好地再現(xiàn)原始 DFT 能帶,而且 Wannier 函數(shù)具有類似原子軌道的局域性,可以很好地和晶體場理論結合。 有了 onsite 和 hopping,你就可以: 看 d 軌道 onsite energies 的相對大小 -> 分析晶體場劈裂(t₂g/eg 等) 看 hopping 強度和符號 -> 分析超交換路徑、帶寬、有效模型(Hubbard 模型、Heisenberg 模型) |
|
您好,我已經通過wannier90計算得到了hr.dat和wout文件,結構并非標準的正八面體,所以存在兩個問題 問題一 nsite energy 的eg和t2g軌道并不簡并(圖中1和4是eg軌道,235為t2g軌道)問題二:后續(xù)如何通過hr.dat得到hopping 向您請教一下計算流程,非常感謝您的答復 發(fā)自小木蟲手機客戶端 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學碩材料275調劑 +4 | 路三三 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境調劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[考研] 085600材料與化工調劑 280分 +10 | yyqqhh 2026-03-03 | 10/500 |
|
|
[考研] 江蘇省農科院招調劑1名 +5 | Qwertyuop 2026-03-01 | 5/250 |
|
|
[考研] 化工專碩348,一志愿985求調劑 +8 | 弗格個 2026-02-28 | 11/550 |
|
|
[考研] 284求調劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研] 268求調劑 +6 | 好運連綿不絕 2026-03-02 | 6/300 |
|
|
[考研] 中國科學技術大學材料與化工281求調劑,有科研和獲獎經歷 +6 | wsxw 2026-03-02 | 7/350 |
|
|
[考研] 求調劑院校 +6 | 云朵452 2026-03-02 | 8/400 |
|
|
[考研] 298求調劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 272求調劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
[考研] 材料復試調劑 +5 | 學材料的點 2026-03-01 | 6/300 |
|
|
[考研] 306分材料調劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 一志愿中石油(華東)本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-02 | 3/150 |
|
|
[考研] 一志愿東北大學材料專碩328,求調劑 +3 | shs1083 2026-03-02 | 3/150 |
|
|
[考研] 291 求調劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
[基金申請] 面上模板改不了頁邊距吧? +6 | ieewxg 2026-02-25 | 7/350 |
|
|
[基金申請] 此成果不能導入原因:元數(shù)據必填信息不完整,可 進行補充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
[考研] 313求調劑 +3 | 水流年lc 2026-02-28 | 3/150 |
|
|
[高分子] 求環(huán)氧樹脂研發(fā)1名 +3 | 孫xc 2026-02-25 | 11/550 |
|