Yambo 光吸收计算后处理

在gengzi服务器计算，不load module，而是使用老师在个人环境下编译好的
首先`source ~/bin/intel/oneapi/setvars.sh`
后期所有命令直接  `~/bin/yambo-23-10-19/bin/yambo -x xxxxxx -x xxxx`
                 `~/bin/yambo-23-10-19/bin/ypp -x xxxxxx -x xxxx`
yambo下的所有脚本以此类推

![](/media/202403/2024-03-18_194940_0406190.7542342215214848.png)

**知识点：**

![](/media/202403/2024-03-18_210645_5341080.09096901440129934.png)
### 1.根据GW计算结果还原能带（Bands_interpolated）

输入文件的产生和修改：` ypp -s b -V all -F ypp_bands.in`

![](/media/202403/2024-03-18_193609_2982670.14095012392162987.png)
计算命令：`ypp -F ypp_bands.in`
输出文件：o.bands_interpolated（把这个文件放到origin里面作图，列名文件里给了，除了最后一列K点名称和倒数三列k点路劲，都拿去画图）
大致走向应和DFT能带差不多
![](/media/202403/2024-03-18_204558_6056720.797052946684198.png)

### 2.激子能量和强度
激子（准粒子，QP)可以有很多，因为在不同能带不同K点间的相互跃迁，会产生不同能量的激子，其实每一个激子都有激子结合能，这里我们一般说第一个激子。
QP带隙-激子能量=激子结合能
光学带隙，通过BSE计算得到，这里可以看到间接带隙更小，所以是间接带隙0.49 eV
但是这里只看直接跃迁，所以取0.57 eV 减去第一个激子的能量就是第一个激子的激子结合能
r-xxxxx_ppa文件中
![](/media/202403/2024-03-18_211719_2332210.8658536606465422.png)

### 3.激子能量和强度 
查看激子的振幅和位置，

![](/media/202403/2024-03-19_134508_6986980.8715818787161997.png)
这里o-*****E_sorted文件里可以看到第一个激子的能量是0.31 eV

![](/media/202403/2024-03-19_134810_2201120.8693409340372773.png)

也就是说，激子结合能是0.57-0.31=0.26 eV

![](/media/202403/2024-03-19_153536_9653470.1511115868037699.png)
查看weight文件可以看到第一个激子是21-22能带第一个k点的跃迁导致的，21-23的第一个k点跃迁能量简并，且gw_gpu/r-2D-********-corr文件中可以看见QP性质，虽然计算时算了18x18个k点，但是简并成37个k点，且第一个k点是 0 0 0，也就是gama点
可以把所有激子的大小和振幅放到图中，是对应的

![](/media/202403/2024-03-19_143659_7429410.5097870315011193.png)

### 3.激子波函数的空间分布

![](/media/202403/2024-03-19_152019_7797330.3025650887400445.png)

Yambo 光吸收计算后处理

服务器已完成数据迁移，数据丢失请及时反馈站长，旧服务器数据保留至11月30日。

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱
work@liuyaoze.com

Yambo 光吸收计算后处理

服务器已完成数据迁移，数据丢失请及时反馈站长，旧服务器数据保留至11月30日。

【温馨提示：本站文档可配置可见范围，如登录后可见、对特定群组可见等，看不到就是没权限】

注册码获取邮箱 work@liuyaoze.com

注册码获取邮箱
work@liuyaoze.com