24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1849  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学070300 求调剂 +25 哈哈哈^_^ 2026-04-12 25/1250 2026-04-16 10:47 by Espannnnnol
[考研] 一志愿沪9,生物学326求调剂 +8 刘墨墨 2026-04-15 8/400 2026-04-16 10:29 by 圆心602
[考研] 294求调剂 +11 淡然654321 2026-04-15 11/550 2026-04-16 10:23 by beilsong20
[考研] 初试324 中药学 一志愿天中医 求调剂 +4 李的Lucky 2026-04-10 4/200 2026-04-16 09:22 by fqwang
[考研] 恳请有学校收留 +6 柯淮然 2026-04-12 6/300 2026-04-16 07:23 by 学员JpLReM
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 085500求调剂材料 +11 易11122 2026-04-09 11/550 2026-04-14 17:59 by lhj2009
[考研] 335求调剂 +19 想上岸呀!! 2026-04-12 21/1050 2026-04-14 16:23 by Art1977
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 求调剂 +12 璃茉一定上岸 2026-04-10 13/650 2026-04-14 00:08 by Equinoxhua
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 求调剂,一志愿材料科学与工程985,365分, +8 材化李可 2026-04-11 10/500 2026-04-12 08:42 by 852137818
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 化学工程与技术324调剂 +23 孙常华 2026-04-09 25/1250 2026-04-11 00:07 by 骑牛渡寒江
[考研] 考研调剂 +26 硕星赴 2026-04-09 27/1350 2026-04-10 22:24 by 猪会飞
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 085601初试330分找调剂 +10 流心奶黄包l 2026-04-09 10/500 2026-04-10 08:14 by Sammy2
信息提示
请填处理意见