Quantum ESPRESSO教程之电子局域函数计算(Pymatflow篇)

本系列Quantum ESPRESSO教程将在计算设置讲解的基础上添加结合Pymatflow工具简化计算流程的内容。话不多说，本文将为您带来Quantum ESPRESSO系列教程中的电子局域函数计算部分。

简介

Qauntum Espresso目前走的发展方向是提高模块化程度。各个计算模块进行独立开发，但是相互之间也可能有以来关系，比如需要计算体系的电子局域函数就需要首先利用pwscf进行体系基态的计算，然后从输出数据中利用pp.x处理得到体系的电子局域函数。

本文将介绍如何使用matflow完成此项工作。

使用Pymatflow高效完成任务

为了完成上述电子局域函数的计算任务，首先需要利用matflow完成静态计算，然后再进行电子局域函数的计算。

静态计算

~$ 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

上述命令会在当前路径创建nacl-static目录，然后在其下自动生成输入文件并进行作业任务的提交。默认使用PBS作业调度器。

你也可以通过--server指定其它服务器类型，比如吕梁天河二号超算系统请使用--server llhpc参数。如果你只是想在本地单节点运行，请使用一下参数组合: --runopt genrun、--auto 0。

电子局域函数的计算

~$ mflow qe -r 9 --cif nacl.cif --plot-num 8 --runopt genrun --auto 0 -d nacl-static

结束后可以在nacl-static路径下找到一个electron-local-function.3d.xsf文件，其可以通过VESTA等软件进行可视化分析。下面是结果展示。

电子局域函数在[001]方向的投影:

电子局域函数在[111]方向的投影:

往期文章

Quantum ESPRESSO教程之态密度能带结构计算(Pymatflow篇)

Cylon-Arch Linux新手的维护帮手

Linux使用图形界面连接WiFi

Linux下挂载使用Windows的NTFS分区

Pymatflow是一个用于第一性原理模拟的工作流自动化管理软件，目前支持的计算程序有VASP、Qutum ESPRESSO、Abinit、SIESTA、CP2K。

前往GitLab搜索pymatflow以查看。

此公众号为Pymatflow开发者所运营，欢迎关注以获取更多关于Pymatflow和第一性原理计算的精彩内容。