QE计算光吸收

组内没有vasp版权，所以所有计算只能用免费的第一性原理软件，我比较常用的是QE。不得不承认，QE有的方面确实没有vasp方便，或者说影响力没有vasp强，所以vasp有了很便捷的vaspkit后处理软件，而qe很多东西只能自己手动处理。 前段时间用QE计算光吸收的过程中遇到了很多问题，上网上搜也没有一个全面的解答，查了很多资料并与老师讨论很久后得到了一个比较系统得过程，分享出来给大家参考。（所有内容都是自己收动验证过的，但人无完人，学生一个，不保证没有错误）。 QE计算光吸收时要先计算出介电函数。计算流程为：relax.in、scf.in、nscf.in、epsilon.in。详细内容在

http://yyyu200.gitee.io/dftbook/blogs/2019/04/01/HandsOn/

网址的3.1章节由具体描述。得到介电函数虚实部后可根据下列公式计算得到光吸收。

得到光吸收后基本计算已完成，但还有一个很重要的部分：单位。vaspkit处理得到的光吸收纵坐标单位为cm-1，横坐标单位为nm或eV。那么QE按上述操作完成后的横纵坐标单位分别是什么呢？ 我常用的单位为：横坐标eV, 纵坐标cm-1。所以整理了由上述公式向该单位变化的一个过程。 首先是横坐标转换：根据以下公式将w与E进项变换

经过转换后的光吸收公式为

括号内为相对介电函数，没有单位，公式推导过程中E单位为eV，h为ev·s，c为m/s，最终化简得纵轴单位为m-1. 此时横坐标单位为eV。纵坐标单位为m-1.从m-1转换为cm-1应该不用详细介绍了。 再次声明：所有内容是网上查找并与老师讨论得到的，不排除有错误的存在，若你发现错误，请指出，我一定改正。感谢。（自己写了个python脚本来处理上述所有过程，之后有时间分享给大家） 吐槽：硕士生活很开心，与硕导关系也很亲近。只在网上看见过学生与导师间有矛盾的新闻。因此考虑过读博后可能会遇到比较难相处的导师。但没想到所有事物都是相关的：学生与导师关系不和，学生逐渐增多，导师为了便于管理开始挑拨学生间关系，随着时间推移，学生间果然产生了不可调和的矛盾，面子工程都做不下去了直接开吵。可真是开了眼了。 共同进步才能整体进步，目前组里各干各的，互帮互助精神很薄弱，所有进展都很慢，很不喜欢这个环境。 第一次投稿不能只说这些不开心的，发个自拍开心一下