本系列Quantum ESPRESSO教程將在計算設置講解的基礎上添加結合Pymatflow工具簡化計算流程的內容。話不多說，本文將為您帶來Quantum ESPRESSO系列教程中的電子局域函數計算部分。
簡介
Qauntum Espresso目前走的發展方向是提高模塊化程度。各個計算模塊進行獨立開發，但是相互之間也可能有以來關係，比如需要計算體系的電子局域函數就需要首先利用pwscf進行體系基態的計算，然後從輸出數據中利用pp.x處理得到體系的電子局域函數。
本文將介紹如何使用matflow完成此項工作。
使用Pymatflow高效完成任務
為了完成上述電子局域函數的計算任務，首先需要利用matflow完成靜態計算，然後再進行電子局域函數的計算。
靜態計算
<code>~$ matflow qe -r 0 --ecutwfc 40 --kpoints-mp-scf 5 5 5 0 0 0 --kpoints-mp-nscf 6 6 6 0 0 0 -d nacl-static --kpath-file kpath-cubic.txt/<code>
上述命令會在當前路徑創建nacl-static目錄，然後在其下自動生成輸入文件並進行作業任務的提交。默認使用PBS作業調度器。
你也可以通過--server指定其它服務器類型，比如呂梁天河二號超算系統請使用--server llhpc參數。如果你只是想在本地單節點運行，請使用一下參數組合: --runopt genrun、--auto 0。
電子局域函數的計算
<code>~$ mflow qe -r 9 --cif nacl.cif --plot-num 8 --runopt genrun --auto 0 -d nacl-static/<code>
結束後可以在nacl-static路徑下找到一個electron-local-function.3d.xsf文件，其可以通過VESTA等軟件進行可視化分析。下面是結果展示。
電子局域函數在[001]方向的投影:
電子局域函數在[111]方向的投影:
往期文章
Pymatflow是一個用於第一性原理模擬的工作流自動化管理軟件，目前支持的計算程序有VASP、Qutum ESPRESSO、Abinit、SIESTA、CP2K。
前往GitLab搜索pymatflow以查看。
此公眾號為Pymatflow開發者所運營，歡迎關注以獲取更多關於Pymatflow和第一性原理計算的精彩內容。
閱讀更多 開源科學計算 的文章
