VASP实用教程：立方BaTiO3和石墨的MLWFs

本文为由小强撰写的《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