| 5 | 1/1 | 返回列表 |
| 查看: 1201 | 回復: 2 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
[求助]
請教關于tight binding hopping參數(shù)的問題 已有1人參與
|
|||
請問怎么從wannier90計算的hr.dat文件中得到tight binding的hopping和SOC強度。新手求助啊 |
vasp+wannier90 |

|
剛剛那個分號不對重發(fā)一遍 de = Flatten[ Take[Import["D:\\wannier\\silicon_hr.dat"], {4, 10}]]; hr = {Take[#, 3], Take[#, {4, 5}], #[[6]] + I #[[7]]} & /@ Drop[Import["D:\\wannier\\silicon_hr.dat"], 10]; hk = Table[{Sum[ E^(I hr[[l ;; -1 ;; 64]][[i, 1]].a.k) (hr[[l ;; -1 ;; 64]][[i, -1]]/de), {i, 93}], hr[[l ;; -1 ;; 64]][[1, 2]]}, {l, 64}]; Re@Eigenvalues[SparseArray[Table[hk[[i, 2]] -> hk[[i, 1]], {i, 64}]]; Normal[hkk] /. k -> {0, 0, 0}] |
|
這里有一個我寫的Mathematica程序 是example03下的 你看看吧 de = Flatten[Take[Import["D:\\wannier\\silicon_hr.dat"], {4, 10}]]; hr = {Take[#, 3], Take[#, {4, 5}], #[[6]] + I #[[7]]} & /@ Drop[Import["D:\\wannier\\silicon_hr.dat"], 10]; hk = Table[{Sum[ E^(I hr[[l ;; -1 ;; 64]][[i, 1]].a.k) (hr[[l ;; -1 ;; 64]][[i, -1]]/de[]), {i, 93}], hr[[l ;; -1 ;; 64]][[1, 2]]}, {l, 64}]; Re@Eigenvalues[SparseArray[Table[hk[[i, 2]] -> hk[[i, 1]], {i, 64}]]; Normal[hkk] /. k -> {0, 0, 0}] 輸出Gamma點的能帶 {9.70555, 8.79934, 8.79933, 8.79932, 6.22852, 6.22851, 6.2285, \ -5.82185} 和band.dat中可以對的上。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招08考數(shù)學 +5 | laoshidan 2026-03-20 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +6 | Zhao anqi 2026-03-22 | 6/300 |
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|
|
[考研] 材料與化工085600,總分304,本科有兩篇sci參與,求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
[考研] 一志愿 西北大學 ,070300化學學碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 中南大學化學學碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
[考研] 0703化學調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|