Phonon主循环
https://www.quantum-espresso.org/Doc/developer_man/node5.html
主程序：phonon.f90
控制聲子計算的參數：phcom.f90=control_ph
控制電聲耦合的參數：elph.f90
trans = .t. 計算聲子的指標
elph = .t. 計算電聲耦合的指標
輸入文件中控制這兩個的量：
trans=.true.控制聲子計算
electron_phonon=’ ’
讀取輸入并設置要計算的量的標志
1.1) 讀取 pw.x 寫入的所有量
1.2) 讀取贗勢數據
確定必須計算的內容。
2.1) 如果還沒有在disk上，計算 q 個點的網格和
所有 q 點的所有modes并保存在disk上 (SD means save on disk)
2.2) 如果要求image 并行化，則在images之間劃分工作
進行恢復運行檢查，檢查 .xml 文件中已有的內容，并將以完成的部分設置為TRUE
開始對于q點的主循環：do_phonon.f90
4.1) 計算所有不依賴于系統響應的量 initialize_ph()
4.2）檢查是否需要進行能帶計算并進行。
注意：以下幾點僅在 q 為 Gamma 時執行。
4.3) 在頻率上開始循環
4.3.1) 計算極化率作為 iu SD 的函數
4.4) 計算對電場的響應
4.5) 計算 epsilon 和 SD
4.6）計算zeu和SD
4.7) 計算電光系數和SD
4.8) 計算 E 的二階響應
4.9) 計算拉曼張量和 SD
結束注
在不可約表示上開始循環
5.1) 計算對不可約表示的響應
代碼解析
8)繼續步驟4的循環，直到所有q點都計算完全。
總結
- 上一篇： 西游记中孙悟空所尊敬的神佛
- 下一篇： php版临时邮箱,Forsaken Ma