24小时热门版块排行榜    

查看: 1689  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 288求调剂 +11 于海海海海 2026-03-19 11/550 2026-03-19 16:30 by 30660438
[考研] 一志愿北京化工大学0703化学318分,有科研经历,求调剂 +3 一瓶苯甲酸 2026-03-14 3/150 2026-03-19 15:17 by 尽舜尧1
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 281求调剂(0805) +9 烟汐忆海 2026-03-16 19/950 2026-03-19 11:42 by laoshidan
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
信息提示
请填处理意见