| 查看: 1277 | 回复: 2 | |||
[求助]
请教关于tight binding hopping参数的问题 已有1人参与
|
请问怎么从wannier90计算的hr.dat文件中得到tight binding的hopping和SOC强度。新手求助啊 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp+wannier90 |
» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有81人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复
» 本主题相关价值贴推荐，对您同样有帮助:
【答案】应助回帖
|
这里有一个我写的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中可以对的上。 |
2楼2015-07-21 23:17:59
|
刚刚那个分号不对重发一遍 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}] |
3楼2015-07-21 23:21:09