| 查看: 1206 | 回复: 2 | |||
[求助]
请教关于tight binding hopping参数的问题 已有1人参与
|
请问怎么从wannier90计算的hr.dat文件中得到tight binding的hopping和SOC强度。新手求助啊 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp+wannier90 |
» 猜你喜欢
计算与实验科学工程国际会议(8月7-11日,香港) 二维材料分会S5与声学换能器分会S8
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有4人回复
-
物理学I论文润色/翻译怎么收费?
已经有197人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有28人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有4人回复
-
欢迎调剂到赣南师范大学智能制造与未来能源学院
已经有26人回复
-
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有4人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有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
