當(dāng)前位置：首頁 > 编程资源 > 编程问答 >内容正文

编程问答

Phonon主循环

發(fā)布時間：2023/12/15 编程问答 39 豆豆

生活随笔收集整理的這篇文章主要介紹了 Phonon主循环小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.

https://www.quantum-espresso.org/Doc/developer_man/node5.html
主程序：phonon.f90
控制聲子計算的參數(shù)：phcom.f90=control_ph
控制電聲耦合的參數(shù)：elph.f90
trans = .t. 計算聲子的指標(biāo)
elph = .t. 計算電聲耦合的指標(biāo)
輸入文件中控制這兩個的量：
trans=.true.控制聲子計算
electron_phonon=’ ’

讀取輸入并設(shè)置要計算的量的標(biāo)志
1.1) 讀取 pw.x 寫入的所有量
1.2) 讀取贗勢數(shù)據(jù)

確定必須計算的內(nèi)容。
2.1) 如果還沒有在disk上，計算 q 個點的網(wǎng)格和
所有 q 點的所有modes并保存在disk上 (SD means save on disk)
2.2) 如果要求image 并行化，則在images之間劃分工作

進行恢復(fù)運行檢查，檢查 .xml 文件中已有的內(nèi)容，并將以完成的部分設(shè)置為TRUE

開始對于q點的主循環(huán)：do_phonon.f90
4.1) 計算所有不依賴于系統(tǒng)響應(yīng)的量 initialize_ph()
4.2）檢查是否需要進行能帶計算并進行。
注意：以下幾點僅在 q 為 Gamma 時執(zhí)行。
4.3) 在頻率上開始循環(huán)
4.3.1) 計算極化率作為 iu SD 的函數(shù)
4.4) 計算對電場的響應(yīng)
4.5) 計算 epsilon 和 SD
4.6）計算zeu和SD
4.7) 計算電光系數(shù)和SD
4.8) 計算 E 的二階響應(yīng)
4.9) 計算拉曼張量和 SD
結(jié)束注

在不可約表示上開始循環(huán)
5.1) 計算對不可約表示的響應(yīng)
代碼解析

phqscf.f90:SUBROUTINE phqscfnpert(:)每個不可約表示的擾動數(shù)對于每個不可約表示，我們計算波函數(shù)的變化 5.1.1）積累對電子-聲子的貢獻(xiàn)5.1.2) 累積對動力學(xué)矩陣的貢獻(xiàn)5.1.3）積累對zue的貢獻(xiàn)5.1.4) SD 這個對動態(tài)矩陣的貢獻(xiàn)和 zue 繼續(xù)循環(huán)

直到當(dāng)前 q 點的所有表示都被計算出來

對動力學(xué)陣進行對角化（僅當(dāng)此 q 的所有表示都已計算時）

對 k 求和，對能帶上的電子-聲子耦合求和以計算 gamma_mat（僅當(dāng)此 q 的所有表示都已計算時）
8)繼續(xù)步驟4的循環(huán)，直到所有q點都計算完全。

總結(jié)

以上是生活随笔為你收集整理的Phonon主循环的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。

phonon