| 查看: 1324 | 回复: 2 | ||||
[求助]
请教关于tight binding hopping参数的问题 已有1人参与
|
请问怎么从wannier90计算的hr.dat文件中得到tight binding的hopping和SOC强度。新手求助啊 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp+wannier90 |
» 猜你喜欢
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有265人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
-
3月 中国石油和化学工业联合会科学技术奖 马上开始，申获奖要点
已经有15人回复
-
中国石油和化学工业联合会科学技术奖的评审标准是什么？
已经有5人回复
-
经费如何使用，花不完怎么办？
已经有5人回复
-
兼职
已经有12人回复
-
海外Gustau Catalan课题组情况了解
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
【答案】应助回帖
|
这里有一个我写的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
10