| 查看: 1182 | 回复: 2 | |||
[求助]
请教关于tight binding hopping参数的问题 已有1人参与
|
请问怎么从wannier90计算的hr.dat文件中得到tight binding的hopping和SOC强度。新手求助啊 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp+wannier90 |
» 猜你喜欢
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有213人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有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中可以对的上。 |
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