【资源之家】每日免费更新最热门的副业项目资源

本文为由小强撰写的《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.【资源之家】每日免费更新最热门的副业项目资源

版权声明 1 本网站名称：资源之家
2 本站永久网址：https://www.ziyuanzj.cn/
3 本网站的文章部分内容可能来源于网络，仅供大家学习与参考，如有侵权，请联系站长 QQ2332379进行删除处理。
4 本站一切资源不代表本站立场，并不代表本站赞同其观点和对其真实性负责。
5 互联网转载资源会有一些其他联系方式，请大家不要盲目相信，被骗本站概不负责！
6 本网站部分内容只做项目揭秘，无法一对一教学指导，每篇文章内都含项目全套的教程讲解，请仔细阅读。
7 本站分享的所有平台仅供展示，本站不对平台真实性负责，站长建议大家自己根据项目关键词自己选择平台。
8 因文章发布时间和您阅读文章时间存在时间差，有些项目红利期可能已经过了，能不能赚钱需要自己判断，本网站仅做资源分享，不做任何收益保障。
9 本站资源大多存储在云盘，如发现链接失效，请联系我们我们会第一时间更新。