24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1997  |  回复: 4
【悬赏金币】回答本帖问题,作者uspex将赠送您 15 个金币

uspex

新虫 (初入文坛)

[求助] wannier90计算得到_hr.dat文件,如何得到onsite和hopping能量 已有1人参与

最近看到一些文献,分析了结构中每个原子每个轨道的onsite和hopping的能量,然后结合晶体场理论做出了一些分析,我想知道是如何得到这些能量的。求大佬们指点一下。
发自小木虫手机客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

河马杀手

新虫 (初入文坛)

您好 我最近也在做这个东西 想向您请教一下 感谢
发自小木虫手机客户端
2楼2025-09-23 13:08:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

至尊木虫 (正式写手)

【答案】应助回帖

“每个原子每个轨道的 onsite 能量和 hopping 能量”,通常是通过从第一性原理计算得到的能带结构,投影到局域轨道基底,构建紧束缚模型(tight-binding model)或 Wannier 函数哈密顿量得到的。
1. 先从 DFT 得到哈密顿量

你先用 DFT(VASP、QE、WIEN2k、ELK…)算出自洽电荷密度和能带。

DFT 本身得到的是平面波基底下的哈密顿量,不是直接的“原子轨道能量和跃迁矩阵”。

2. 局域化:投影到原子轨道基底

有两条主流方法:

方法 A:直接投影到原子轨道 (Projected Tight Binding)
对平面波波函数做投影,得到每个原子轨道(如 d_xy, d_xz, p_x)的投影系数。
在这个基底上构造哈密顿量矩阵 𝐻𝑖𝑗。对角元 𝐻𝑖𝑖 就是所谓 onsite energy。
非对角元𝐻𝑖𝑗 (i≠j)就是原子 i 和 j 之间的 hopping 参数(跃迁能)。
这种方法计算比较直接,但对轨道正交性要求比较高。

方法 B:Wannierization (Maximally Localized Wannier Functions, MLWF)
先选择一组目标能带(比如价带+导带的 d/p 带)。
用 Wannier90 等工具生成局域 Wannier 函数。
在 Wannier 基底下的哈密顿量是稀疏的,矩阵元就是:
对角元 = onsite energy
非对角元 = hopping (可能包括长程跳跃)
好处:Wannier 化的哈密顿量可以很好地再现原始 DFT 能带,而且 Wannier 函数具有类似原子轨道的局域性,可以很好地和晶体场理论结合。
3楼2025-09-24 00:09:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

至尊木虫 (正式写手)

【答案】应助回帖

“每个原子每个轨道的 onsite 能量和 hopping 能量”,通常是通过从第一性原理计算得到的能带结构,投影到局域轨道基底,构建紧束缚模型(tight-binding model)或 Wannier 函数哈密顿量得到的。
1. 先从 DFT 得到哈密顿量
你先用 DFT(VASP、QE、WIEN2k、ELK…)算出自洽电荷密度和能带。
DFT 本身得到的是平面波基底下的哈密顿量,不是直接的“原子轨道能量和跃迁矩阵”。

2. 局域化:投影到原子轨道基底
有两条主流方法:

方法 A:直接投影到原子轨道 (Projected Tight Binding)
对平面波波函数做投影,得到每个原子轨道(如 d_xy, d_xz, p_x)的投影系数。
在这个基底上构造哈密顿量矩阵H_ij;。对角元 H_ii; 就是所谓 onsite energy。
非对角元H_ij; (i≠j)就是原子 i 和 j 之间的 hopping 参数(跃迁能)。
这种方法计算比较直接,但对轨道正交性要求比较高。

方法 B:Wannierization (Maximally Localized Wannier Functions, MLWF)
先选择一组目标能带(比如价带+导带的 d/p 带)。
用 Wannier90 等工具生成局域 Wannier 函数。
在 Wannier 基底下的哈密顿量是稀疏的,矩阵元就是:
对角元 = onsite energy
非对角元 = hopping (可能包括长程跳跃)
好处:Wannier 化的哈密顿量可以很好地再现原始 DFT 能带,而且 Wannier 函数具有类似原子轨道的局域性,可以很好地和晶体场理论结合。

有了 onsite 和 hopping,你就可以:
看 d 轨道 onsite energies 的相对大小 -> 分析晶体场劈裂(t₂g/eg 等)
看 hopping 强度和符号 -> 分析超交换路径、带宽、有效模型(Hubbard 模型、Heisenberg 模型)
4楼2025-09-24 00:19:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

河马杀手

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by rlafite at 2025-09-24 00:19:01
“每个原子每个轨道的 onsite 能量和 hopping 能量”,通常是通过从第一性原理计算得到的能带结构,投影到局域轨道基底,构建紧束缚模型(tight-binding model)或 Wannier 函数哈密顿量得到的。
1. 先从 DFT 得到 ...

您好,我已经通过wannier90计算得到了hr.dat和wout文件,结构并非标准的正八面体,所以存在两个问题
问题一nsite energy 的eg和t2g轨道并不简并(图中1和4是eg轨道,235为t2g轨道)
问题二:后续如何通过hr.dat得到hopping
向您请教一下计算流程,非常感谢您的答复
wannier90计算得到_hr.dat文件,如何得到onsite和hopping能量


发自小木虫手机客户端
5楼2025-09-24 10:33:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 uspex 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 3/150 2026-05-17 00:45 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 v9tggjlwd0 2026-05-15 5/250 2026-05-17 00:32 by xiangfeng
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[基金申请] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 00:10 by ue3ir18jc3
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 v9tggjlwd0 2026-05-15 3/150 2026-05-16 23:57 by ue3ir18jc3
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 l7k6xnh0yc 2026-05-14 5/250 2026-05-16 23:05 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-16 22:25 by lindoudou
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-16 18:45 by j6b2pdz07o
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 6/300 2026-05-16 16:17 by 0i5p09z61n
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-16 13:57 by vcdazktkjx
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +11 dd921ww 2026-05-12 13/650 2026-05-16 09:16 by Howard28
[教师之家] 上海大学实验技术岗位非升即走 +5 嘻嘻哈哈乐呵呵 2026-05-15 5/250 2026-05-16 00:17 by caiyun
[文学芳草园] 风把牡丹吹跑了 +5 myrtle 2026-05-12 9/450 2026-05-15 15:27 by myrtle
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见