24小时热门版块排行榜    

查看: 5114  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

园中园

新虫 (正式写手)

[求助] vasp2wannier90计算hr.dat里面数据怎么理解

假设 lattice vectors:  a    0     0
                                   0    b    0
                                    0    0    c

据图片可得知 0 0 -2 是 R的分量  ,  怎么计算R 在原胞的值

vasp2wannier90计算hr.dat里面数据怎么理解
1.jpg
回复此楼
追求,奋斗,再追求,再奋斗
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

园中园

新虫 (正式写手)

引用回帖:
2楼: Originally posted by liqizuiyang at 2017-09-16 10:46:53
tight-binding模型只需要<0,i|H|R,j>的值。前三个整数是R,后两个是i和j,最后两个是哈密顿量矩阵元的实部和虚部。

谢谢您了 ,我上面说错了, 图片中第三行数字 17 意思是  一共有17个 Wigner-Seitz grid-points ;这个怎么与tight-binding联系起来啊? 建tight-binding模型时,我需要取哪一个 Wigner-Seitz grid-poin
追求,奋斗,再追求,再奋斗
4楼2017-09-16 13:24:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

liqizuiyang

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
蛋蛋小童鞋: 应助指数+1 2017-09-16 11:21:47
园中园: 金币+5, ★★★很有帮助 2017-09-16 12:48:19
tight-binding模型只需要<0,i|H|R,j>的值。前三个整数是R,后两个是i和j,最后两个是哈密顿量矩阵元的实部和虚部。

» 本帖已获得的红花(最新10朵)

2楼2017-09-16 10:46:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

园中园

新虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by liqizuiyang at 2017-09-16 10:46:53
tight-binding模型只需要<0,i|H|R,j>的值。前三个整数是R,后两个是i和j,最后两个是哈密顿量矩阵元的实部和虚部。

谢谢了,您太厉害,就是要建tight-binding模型,前三个整数是R,那么R 与晶格基矢的关系是什么啊?

wannier 输出的 hr.dat是 实空间的 Ri 点出的哈密顿量矩阵元,我该怎么让它变成倒格矢空间的哈密顿量矩阵元
追求,奋斗,再追求,再奋斗
3楼2017-09-16 13:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

引用回帖:
3楼: Originally posted by 园中园 at 2017-09-16 13:01:17
谢谢了,您太厉害,就是要建tight-binding模型,前三个整数是R,那么R 与晶格基矢的关系是什么啊?

wannier 输出的 hr.dat是 实空间的 Ri 点出的哈密顿量矩阵元,我该怎么让它变成倒格矢空间的哈密顿量矩阵元...

倒空间矩阵元等于实空间矩阵元乘以相位因子exp(i*K*R)知乎求和,公式是Hij(K) = Sum[ exp(i*K*R) * <0,i|H|R,j> ],i,j,k取固定值,R是求和变量,取所有近邻原子所在的原胞。

相位因子算例:

假设R = [1,1,-1],代表 R = a1 + a2 - a3。再假设K = [-1/3, 1/3, 0],代表K = -1/3 * b1 + 1/3 * b2 + 0 * b3。
相位因子等于 exp(i * R * K) = exp( i * 2 * pi * (1 * -1/3 + 1 * 1/3) ) = exp(i * 2 * pi * 0) = exp(0) = 1。
5楼2017-09-16 13:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见