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

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

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

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

經(jīng)過(guò)轉(zhuǎn)換后的光吸收公式為

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