本文为由小强撰写的《VASP实用教程》第46篇，全系列约60篇，将在近期陆续更新。
1. 立方BaTiO3
- 概要：获取钙钛矿的MLWFs
- 目录：examples/example9/
- 输入文件
– batio3.scf用于基态计算的PWSCF输入文件。
– batio3.nscf 用于获取统一网格上的Bloch状态的PWSCF输入文件。
– batio3.pw2wan pw2wannier90的输入文件
– batio3.win wannier90的输入文件
首先，我们要获得氧2p态的MLWFs。从能带结构来看，这些形成了一组孤立的带。我们使用wannier90关键字exclude_bands从重叠和投影矩阵的计算中删除除2p带以外的所有带（非必须，但可以节省时间）。
【1】运行PWSCF获得BaTiO3的基态。
pw.x < BaTiO3.scf > scf.out
【2】运行PWSCF，在统一的k点网格上获得Bloch状态。
pw.x < BaTiO3.nscf > nscf.out
【3】运行wannier90来生成所需的重叠列表（写入BaTiO3.nnkp文件中）。
wannier90.x -pp BaTiO3
【4】运行pw2wannier90来计算Bloch状态和初猜的投影之间的重叠（写在BaTiO3.mmn和BaTiO3.amn文件中）。
pw2wannier90.x < BaTiO3.pw2wan > pw2wan.out
【5】运行wannier90来计算MLWFs。
wannier90.x BaTiO3
检查输出文件BaTiO3.wout。
通过向输入文件BaTiO3.win中添加以下关键字，绘制第二个MLWF。
wannier_plot = true
restart = plot
wannier_plot_list = 2
wannier_plot_supercell = 3
并重新运行wannier90。使用XCrySDen将其可视化。
xcrysden –xsf BaTiO3_00002.xsf
现在我们可以通过置换Ti原子来模拟铁电相。将其位置改为
Ti 0.505 0.5 0.5
并再生MLWF(即用PWSCF等计算地态电荷密度和布洛赫态)，再看第二个MLWF的变化。
2.石墨
- 概要：获得石墨的MLWFs（AB，Bernal）。
- 目录：examples/example10/
- 输入文件
– graphite.scf 用于基态计算的PWSCF输入文件。
– graphite.nscf 在统一网格上获得Bloch状态的PWSCF输入文件。
– graphite.pw2wan pw2wannier90的输入文件
– graphite.win wannier90的输入文件
【1】运行PWSCF，获得石墨的基态。
pw.x < graphite.scf > scf.out
【2】运行PWSCF，在统一的k点网格上获得Bloch状态。
pw.x < graphite.nscf > nscf.out
【3】运行wannier90来生成所需的重叠列表（写入graphite.nnkp文件中）。
wannier90.x -pp graphite
【4】运行pw2wannier90来计算Bloch状态和初猜的投影之间的重叠（写在graphite.mmn和graphite.amn文件中）。
pw2wannier90.x < graphite.pw2wan > pw2wan.out
【5】运行wannier90来计算MLWFs。
wannier90.x graphite
