組內沒有vasp版權，所以所有計算只能用免費的第一性原理軟件，我比較常用的是QE。不得不承認，QE有的方面確實沒有vasp方便，或者說影響力沒有vasp強，所以vasp有了很便捷的vaspkit后處理軟件，而qe很多東西只能自己手動處理。 前段時間用QE計算光吸收的過程中遇到了很多問題，上網(wǎng)上搜也沒有一個全面的解答，查了很多資料并與老師討論很久后得到了一個比較系統(tǒng)得過程，分享出來給大家參考。（所有內容都是自己收動驗證過的，但人無完人，學生一個，不保證沒有錯誤）。 QE計算光吸收時要先計算出介電函數(shù)。計算流程為：relax.in、scf.in、nscf.in、epsilon.in。詳細內容在

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

網(wǎng)址的3.1章節(jié)由具體描述。得到介電函數(shù)虛實部后可根據(jù)下列公式計算得到光吸收。

得到光吸收后基本計算已完成，但還有一個很重要的部分：單位。vaspkit處理得到的光吸收縱坐標單位為cm-1，橫坐標單位為nm或eV。那么QE按上述操作完成后的橫縱坐標單位分別是什么呢？ 我常用的單位為：橫坐標eV, 縱坐標cm-1。所以整理了由上述公式向該單位變化的一個過程。 首先是橫坐標轉換：根據(jù)以下公式將w與E進項變換

經過轉換后的光吸收公式為

括號內為相對介電函數(shù)，沒有單位，公式推導過程中E單位為eV，h為ev·s，c為m/s，最終化簡得縱軸單位為m-1. 此時橫坐標單位為eV。縱坐標單位為m-1.從m-1轉換為cm-1應該不用詳細介紹了。 再次聲明：所有內容是網(wǎng)上查找并與老師討論得到的，不排除有錯誤的存在，若你發(fā)現(xiàn)錯誤，請指出，我一定改正。感謝。（自己寫了個python腳本來處理上述所有過程，之后有時間分享給大家） 吐槽：碩士生活很開心，與碩導關系也很親近。只在網(wǎng)上看見過學生與導師間有矛盾的新聞。因此考慮過讀博后可能會遇到比較難相處的導師。但沒想到所有事物都是相關的：學生與導師關系不和，學生逐漸增多，導師為了便于管理開始挑撥學生間關系，隨著時間推移，學生間果然產生了不可調和的矛盾，面子工程都做不下去了直接開吵?？烧媸情_了眼了。 共同進步才能整體進步，目前組里各干各的，互幫互助精神很薄弱，所有進展都很慢，很不喜歡這個環(huán)境。 第一次投稿不能只說這些不開心的，發(fā)個自拍開心一下