本文为由小强撰写的《VASP实用教程》第56篇，全系列约60篇，将在近期陆续更新。
- 概要：使用SCDM方法计算钨（W）的Wannier内插带结构，以计算初猜（详见例27）。SCDM方法中的自由参数，即µ和σ，是通过拟合一个可投影性的互补误差函数来获得的。MLWFs的数量由赝原子轨道（PAO）的数量决定，本例中为21个。本例中的所有步骤都是在AiiDA[21]工作流程中自动完成的，可以从MaterialsCloud网站下载[1]。
- 目录：examples/example31/
- 输入文件
– W.scf 用于基态计算的PWSCF输入文件– W.nscf 在统一网格上获得Bloch状态的PWSCF输入文件– W.pw2wan pw2wannier90的输入文件– W.proj projwfc的输入文件– generate_weights.sh 从projwfc输出中提取projectabilities的bash脚本– W.win wannier90的输入文件【1】运行PWSCF以获得钨的基态pw.x -in W.scf > scf.out【2】在10×10×10的均匀k点网格上运行PWSCF来获得Bloch状态pw.x -in W.nscf > nscf.out【3】运行wannier90生成所需的重叠列表（写入W.nnkp文件中）wannier90.x -pp W【4】运行projwfc来计算Bloch状态对Bloch总和的投射能力，从赝势中的PAO得到的Bloch的和projwfc.x -in W.proj > proj.out【5】运行generate_weights，将proj.out中的projectabilitites提取为适合Xmgrace或gnuplot读取的格式./generate_weights.sh【6】绘制投影图，用互补的误差函数拟合数据
我们要用Xmgrace来绘制投影并进行拟合。打开Xmgracexmgrace要导入p_vs_e.dat文件，请从顶栏点击数据，然后Import -> ASCII….。这时，一个新的窗口Grace: 读取集应弹出。在文件部分选择p_vs_e.dat，点击底部的OK，然后关闭窗口。现在你应该可以看到一个相当粗糙的函数，它的边界在1和0之间。你可以通过点击顶部栏中的Plot，然后点击Set appearance….，来修改该图的外观。在弹出窗口的主要部分，将符号类型从无改为圆。将线条类型从直线改为无，因为Xmgrace默认添加的线条是无意义的。关于拟合，请进入Data -> Transformations -> Non-linear curve fitting。在这个窗口中，在设置框中选择源，在公式框中插入以下内容y = 0.5 * erfc( ( x – A0 ) / A1 )选择2作为参数数，给A0的初始条件是40，给A1的初始条件是7。点击应用。一个新的窗口会弹出，显示拟合的统计数据。特别是应该发现一个相关系数为0.96，A0的值为39：9756，A1的值为6：6529。这些是我们将用于SCDM方法的拟合值和适合值。具体而言，µSCDM = µfit-3σfit = 20.0169 eV和σSCDM = σfit = 6:6529 eV。这个特定的µfit和σfit选择的动机可以在文献[23]中找到。其中作者还展示了在200种材料的数据集上验证了这种方法。现在应该看到拟合函数了，以及投影能力，见图10-（a））。【7】打开W.pw2wan并添加以下几行scdm_entanglement = 'erfc'scdm_mu = 20.0169scdm_proj = .true.scdm_sigma = 6.6529/【8】运行pw2wannier90来计算Bloch状态和初猜的投影之间的重叠部分（写在W.mmn和W.amn文件中）pw2wannier90.x -in W.pw2wan > pw2wan.out【9】运行wannier90以获得插值的带状结构（见图10-(b)）。wannier90.x W在任何出版物中采用本例中概述的程序获得σ和µ，请引用参考文献[2]。
图10:a）每个蓝点表示可投射性，如参考文献[23]中等式（22）所定义的状态|nk〉，作为钨的相应能量nk的函数。黄线表示拟合的互补误差函数。垂直红线表示σfit的值，垂直绿线表示µSCDM的最佳值，即µSCDM=µfit-3σfit。b）钨的能带结构Γ-H-N-Γ，从DFT计算得到的路径（实心黑色）和使用SCDM方法构造初始猜测（红点）的Wannier插值。参考文献：V. Vitale, G. Pizzi, A. Marrazzo, J. R. Yates, N. Marzari, and A. A. Mostofi, Materials Cloud Archive (2019), doi:10.24435/materialscloud:2019.0044/v2.V. Vitale, G. Pizzi, A. Marrazzo, J. Yates, N. Marzari, and A. Mostofi, “Automated high throughput wannierisation,” (2019), arXiv:1909.00433.【资源之家】每日免费更新最热门的副业项目资源
