spin polarized case

This page describes the basic workflow of using the VASP2WANNIER90 interface to generate Wannier functions.

Steps

• SCF
• Band_analysis
• Wannierzation

SCF

Keywords needed:

ISPIN = 2
LCHARG = .TRUE.
LWAVE = .TRUE.

Keywords optional:

• smearing:

ISMEAR = 0
SIGMA = xxx

• orbital projections:

LORBIT = 11

Run a spin-polarized SCF calculation with ISPIN = 1. This can be done with the standard version vasp_std or the non-collinear version vasp_ncl

The goal of this step is to generate self-consistent Bloch wavefunctions (in WAVECAR) and charge density (in CHGCAR).

Band_analysis

keywords needed:

ISPIN = 2
ICHARG = 11
ISMEAR = 0
SIGMA = xxx
LORBIT = .TRUE.

Then, using p4vasp or pyband you can plot the fat-bands and analysis which bands you want to wannierize and, additionally, the energy windows to use.

Wannierzation

• keywords needed:

ISPIN = 2
ISTART = 1
ICHARG = 2
LWANNIER90 = .TRUE.

• files needed:
  • wannier90.up.win
  • wannier90.dn.win
  • WAVECAR
  • CHGCAR
  • INCAR
  • KPOINTS
  • POSCAR
  • POTCAR

The wannier90.up.win and wannier90.dn.win files should at least contain:

num_wann        =  ??

begin projections
??
end projections

ISTART=1 and ICHARG=1 tell VASP to read in the Bloch wavefunctions and charge density files. LWANNIER90=.TRUE. tells VASP to run the interface with default specs.

After the job is done, there should be a block in the stdout look like:

 *---------------------------------------*
 *             VASP2WANNIER90            *
 *---------------------------------------*
 *          Collinear spin -> 'up'       *
 *            MMN calculated.            *
 *          Projection   1 done.         *
 *          Projection   2 done.         *
 *          Projection   3 done.         *
 *          Projection   4 done.         *
 *          Projection   5 done.         *
 *          Projection   6 done.         *
 *          Projection   7 done.         *
 *          Projection   8 done.         *
 *          Projection   9 done.         *
 *          Projection  10 done.         *
 *          Projection  11 done.         *
 *          Projection  12 done.         *
 *            AMN calculated.            *
 *            EIG written.               *
 *---------------------------------------*
 *          Collinear spin -> 'dn'       *
 *            MMN calculated.            *
 *          Projection   1 done.         *
 *          Projection   2 done.         *
 *          Projection   3 done.         *
 *          Projection   4 done.         *
 *          Projection   5 done.         *
 *          Projection   6 done.         *
 *          Projection   7 done.         *
 *          Projection   8 done.         *
 *          Projection   9 done.         *
 *          Projection  10 done.         *
 *          Projection  11 done.         *
 *          Projection  12 done.         *
 *            AMN calculated.            *
 *            EIG written.               *
 *---------------------------------------*

And you should have:

• wannier90.up.mmn and wannier90.dn.mmn
• wannier90.up.amn and wannier90.dn.amn
• wannier90.up.eig and wannier90.dn.eig

in your directory.

VASP will automatically fill in other necessary blocks in the wannier90.up.win and wannier90.dn.win files.