Bandstructure of SrVO3 in GW: Difference between revisions

From VASP Wiki
VisualWikitext
 
(98 intermediate revisions by 4 users not shown)
Line 1: Line 1:
Description: the GW bandstructure of SrVO<sub>3</sub> using VASP and [http://www.wannier.org WANNIER90].
{{Template:GW - Tutorial}}
== Task ==
 
Calculation of the GW bandstructure of SrVO<sub>3</sub> using VASP and [http://www.wannier.org WANNIER90].


----
----
Line 6: Line 9:
In this example we will also see how the results of the GW calculation may be postprocessed with [http://www.wannier.org WANNIER90] to obtain the dispersion of the bands along the usual high symmetry directions in reciprocal space.
In this example we will also see how the results of the GW calculation may be postprocessed with [http://www.wannier.org WANNIER90] to obtain the dispersion of the bands along the usual high symmetry directions in reciprocal space.


'''Practical issues''':
'''N.B.:''' This example involves quite a number of individual calculations.
#We have trimmed the calculations so that they run smoothly on 6 cores, to run VASP execute the following command: <pre>mpirun -np 6 /opt/tutorial/gw/vasp.5.3/vasp</pre>
The easiest way to run this example is to execute:
#To be able to use the scripts provided with this example to visualize (parts of) your results, you will have to enable ''X11-forwarding'' when you connect to your compute node with ''ssh'':<pre>ssh -X username@machinename</pre>
./doall.sh
 
== The DFT groundstate calculation ==
 
Everthing starts with a conventional DFT (in this case LDA) groundstate calculation:


*INCAR.DFT
And compare the output of the different steps (DFT, GW, HSE) by:
<pre>
  ./plotall.sh
System = SrVO3


NBANDS = 36
In any case, one can consider the <tt>doall.sh</tt> script to be an overview of the steps described below.


ISMEAR = -5
== DFT groundstate calculation ==
EMIN = -20 ; EMAX = 20 ; NEDOS = 1000  # usefull energy range for density of states
The first step is a conventional DFT (in this case PBE) groundstate calculation.


EDIFF = 1E-8                          # high precision for groundstate calculation
*{{TAG|INCAR}} (see INCAR.DFT)
 
KPAR = 3
{{TAGBL|SYSTEM}}  = SrVO3                        # system name
</pre>
{{TAGBL|NBANDS}} = 36                            # small number  of bands
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|EDIFF}} = 1E-8                          # high precision for groundstate calculation
{{TAGBL|KPAR}} = 2                              # parallelization of k-points in two groups


Copy the aforementioned file to {{TAG|INCAR}}:
Copy the aforementioned file to {{TAG|INCAR}}:
Line 32: Line 33:
  cp INCAR.DFT INCAR
  cp INCAR.DFT INCAR


*KPOINTS
The {{TAG|POSCAR}} file describes the structure of the system:
<pre>
Automatically generated mesh
      0
Gamma
4 4 4
</pre>
 
'''Mind''': this is definitely not dense enough for a high-quality description of SrVO<sub>3</sub>, but in the interest of speed we will live with it.


*POSCAR
*{{TAG|POSCAR}}
<pre>
<pre>
SrVO3
SrVO3
3.77706 #taken from 9x9x9 with sigma=0.2 ismear=2
3.84652 #cubic fit for 6x6x6 k-points
  +1.0000000000  +0.0000000000  +0.0000000000  
  +1.0000000000  +0.0000000000  +0.0000000000  
  +0.0000000000  +1.0000000000  +0.0000000000  
  +0.0000000000  +1.0000000000  +0.0000000000  
Line 58: Line 51:
  +0.0000000000  +0.5000000000  +0.5000000000
  +0.0000000000  +0.5000000000  +0.5000000000
</pre>
</pre>
The {{TAGBL|POSCAR}} can be visualized with [http://p4vasp.at p4v] or [http://www.jp-minerals.org/vesta/en/download.html VESTA] and remains unchanged in the following.


=== Analysis of the DOS ===
[[File:SrVO3_structure.png|340px]]
 
Add the following line to your {{FILE|INCAR}} file:
 
LORBIT = 11


and rerun VASP.
The {{TAG|KPOINTS}} file describes how the first Brillouin zone is sampled.
 
In the first step we use a uniform k-point sampling:
In addition to the total density-of-states (DOS), the {{FILE|DOSCAR}} file now contains blocks of information with the site-projected ''lm''-decomposed DOS as well, and the site-projected ''lm''-decomposed band character is written to the {{FILE|PROCAR}} file.
*{{TAG|KPOINTS}} (see KPOINTS.BULK)
 
To plot the total DOS and the Vanadium ''t<sub>2g</sub>'' and ''e<sub>g</sub>'' partial-DOS using ''gnuplot'', execute the following command:
 
./plotdos
 
'''Mind''': Check the {{FILE|OUTCAR}} file for the position of the Fermi level. These DOSs have not been shifted such that the Fermi level is at 0 eV.
 
=== Bandstructure using WANNIER90 ===
 
Add the following line to your {{FILE|INCAR}} to have VASP call [http://www.wannier.org WANNIER90]:
<pre>
<pre>
LWANNIER90_RUN = .TRUE.
Automatically generated mesh
      0
Gamma
4 4 4
0 0 0
</pre>
</pre>


[http://www.wannier.org WANNIER90] takes its input from the file {{FILE|wannier90.win}}.
'''Mind''': this is definitely not dense enough for a high-quality description of SrVO<sub>3</sub>, but in the interest of speed we will live with it.
To construct Wannier functions for the Vanadium ''t<sub>2g</sub>'' manifold in SrVO<sub>3</sub>, and plot the dispersion of the associated bands along R-G-X-M, one may use the following settings:
Copy the aforementioned file to {{TAG|KPOINTS}}:


*wannier90.win.dft
cp KPOINTS.BULK KPOINTS
<pre>
bands_plot = true


begin kpoint_path
and run VASP. If all went well, one should obtain a {{TAG|WAVECAR}} file containing the PBE wavefunction.
R  0.50000000  0.50000000  0.50000000  G  0.00000000  0.00000000  0.00000000
G  0.00000000  0.00000000  0.00000000  X  0.50000000  0.00000000  0.00000000
X  0.50000000  0.00000000  0.00000000  M  0.50000000  0.50000000  0.00000000
M  0.50000000  0.50000000  0.00000000  G  0.00000000  0.00000000  0.00000000
end kpoint_path


num_wann =   3
== Obtain DFT virtual orbitals and long-wave limit ==
Use following {{TAG|INCAR}} file to increase the number of virtual states and to determine the long-wave limit of the polarizability (stored in {{TAG|WAVEDER}}):


num_bands=   3
*{{TAG|INCAR}} (see INCAR.PBE)
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|KPAR}} = 2                              # parallelization of k-points in two groups
{{TAGBL|ALGO}} = Exact                          # exact diagonalization
{{TAGBL|NELM}} = 1                              # one electronic step suffices, since WAVECAR from previous step is present
{{TAGBL|NBANDS}} = 96                            # need for a lot of bands in GW
{{TAGBL|LOPTICS}} = .TRUE.                      # we need d phi/ d k  for GW calculations for long-wave limit


exclude_bands : 1-20, 24-36
Restart VASP.
At this stage it is a good idea to make a safety copy of the {{TAG|WAVECAR}} and {{TAG|WAVEDER}} files since we will repeatedly need them in the calculations that follow:  
cp {{TAGBL|WAVECAR}} WAVECAR.PBE
cp {{TAGBL|WAVEDER}} WAVEDER.PBE
Also make a backup of the charge density for later:
cp {{TAGBL|CHGCAR}} CHGCAR.PBE


begin projections
=== The dielectric function ===
V:dxy;dxz;dyz
end projections
</pre>


Copy the above to {{FILE|wannier90.win}}:
As a bonus, VASP determines the frequency dependent dielectric function in the independent-particle (IP) picture and writes the result to the {{FILE|OUTCAR}} and {{FILE|vasprun.xml}} files.
In the {{FILE|OUTCAR}} you should search for


cp wannier90.win.dft wannier90.win
  frequency dependent IMAGINARY DIELECTRIC FUNCTION (independent particle, no local field effects)


and restart VASP.
and


The Vanadium ''t<sub>2g</sub>'' band dispersion thus obtained, may conveniently be visualized with ''gnuplot'':
  frequency dependent      REAL DIELECTRIC FUNCTION (independent particle, no local field effects)


  gnuplot -persist plotme.dft
== GW Step ==
The actual GW calculation requires a set of one-electron energies and eigenstates. In this case we use the PBE solution obtained from previous step:
  cp WAVECAR.PBE {{TAGBL|WAVECAR}}
cp WAVEDER.PBE {{TAGBL|WAVEDER}}


'''Mind''': Here the eigenvalues have been shifted such that the Fermi level is a 0 eV.
The following {{TAG|INCAR}} file selects the 'single shot' GW calculation also known as G<sub>0</sub>W<sub>0</sub>:
*{{TAG|INCAR}} (see INCAR.GW0)
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|KPAR}} = 2                              # parallelization of k-points in two groups
{{TAGBL|ALGO}} = GW0                            # GW with iteration in G, W kept on DFT level
{{TAGBL|NELM}} = 1                              # one electronic step suffices, since WAVECAR from previous step is present
{{TAGBL|NBANDS}} = 96                            # need for a lot of bands in GW
{{TAGBL|PRECFOCK}} = Fast                        # fast mode for FFTs
{{TAGBL|ENCUTGW}} = 100                          # small energy cutoff for response function suffices for this tutorial
{{TAGBL|NOMEGA}} = 200                          # large number of real frequency points for Hilbert transforms of W and self-energy


== Obtain DFT virtual orbitals ==
Restarting VASP will overwrite the present {{TAGBL|WAVECAR}} and {{TAGBL|vasprun.xml}} file. Make a copy them for later.
cp {{TAGBL|WAVECAR}} WAVECAR.GW0
cp {{TAGBL|vasprun.xml}} vasprun.GW0.xml


*INCAR.DFT.all
== HSE hybrid functional ==
<pre>
System  = SrVO3


ISMEAR = -5
To illustrate the kind of results one would obtain for SrVO<sub>3</sub> using the [[:Category:Hybrid_functionals|DFT/Hartree-Fock hybrid functional]] HSE, without actually doing a full self-consistent calculation, we will recalculate the one-electron energies and DOS ({{TAG|ALGO}}=Eigenval) using the HSE functional with DFT orbitals as input
EMIN = -20 ; EMAX = 20 ; NEDOS = 1000 # usefull energy range for density of states
  cp WAVECAR.PBE {{TAGBL|WAVECAR}}


ALGO = Exact ; NELM = 1               # exact diagonalization one step suffices
Use the following {{TAG|INCAR}} file:
EDIFF = 1E-8                          # high precision for groundstate calculation
*{{TAG|INCAR}} (see INCAR.HSE)
NBANDS = 96                            # need for a lot of bands in GW
LOPTICS = .TRUE.                      # we need d phi/ d k for GW calculations
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|KPAR}} = 2                              # parallelization of k-points in two groups
{{TAGBL|ALGO}} = Eigenval                        # calulate eigenvalues
  {{TAGBL|NELM}} = 1                               # one electronic step suffices, since WAVECAR from previous step is present
{{TAGBL|NBANDS}} = 48                            # small number of bands suffice
{{TAGBL|PRECFOCK}} = Fast                        # fast mode for FFTs
{{TAGBL|LHFCALC}} = .TRUE.                      # switch on Hartree-Fock routines to calculate exact exchange
  {{TAGBL|HFSCREEN}} = 0.2                        # HSE06 screening parameter


KPAR = 3
Restart VASP and make a copy of the wavefunction for post-processing
</pre>
cp {{TAGBL|WAVECAR}} WAVECAR.HSE


Copy the aforementioned file to {{TAG|INCAR}}:
== Post-processing: Density of states and Bandstructure for PBE, GW and HSE ==


  cp INCAR.DFT.all INCAR
=== Density of States ===
The DOS of the PBE, GW and HSE solution can be calculated in a post-processing step with
*{{TAG|INCAR}} (see INCAR.DOS)
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = -5                            # Bloechl's tetrahedron method (requires at least 3x3x3 k-points)
{{TAGBL|ALGO}} = NONE                            # no electronic changes required
{{TAGBL|NELM}} = 1                              # one electronic step suffices, since WAVECAR from previous step is present
{{TAGBL|NBANDS}} = 48                            # number of bands used
{{TAGBL|EMIN}} = -20 ; {{TAGBL|EMAX}} = 20                # smallest/largest energy included in calculation
{{TAGBL|NEDOS}} = 1000                          # sampling points for DOS
{{TAGBL|LORBIT}} = 11                            # calculate l-m decomposed DOS
{{TAGBL|LWAVE}} = .FALSE.                        # do not overwrite WAVECAR
{{TAGBL|LCHARG}} = .FALSE.                       # do not overwrite CHGCAR


and restart VASP.
and requires the apropriate {{TAG|WAVECAR}} file from one of the previous steps. Copy
cp WAVECAR.PBE WAVECAR
or
cp WAVECAR.GW0 WAVECAR
or
cp WAVECAR.HSE WAVECAR
and restart VASP. The density of states is written to {{TAG|DOSCAR}}, make a copy of this file
cp {{TAGBL|DOSCAR}} DOSCAR.XXX
where XXX is either PBE, GW0 or HSE. Visualize the projected DOS for the V-t2g, V-eg and O-p states with the scriptfile
./plotdos_comp.sh DOSCAR.???
This requires gnuplot to be installed.


At this stage it is a good idea to make a safety copy of the {{FILE|WAVECAR}} and {{FILE|WAVEDER}} files since we will repeatedly need them in the calculations that follow:
[[File:DOS_SrVO3_comparison.png|480px]]


  cp WAVECAR WAVECAR.DFT.96bands
=== Bandstructure with <tt>wannier90</tt>===
  cp WAVEDER WAVEDER.DFT.96bands
The bandstructure can be calculated via Wannier interpolation using <tt>wannier90</tt> in the library mode
*{{TAG|INCAR}} (see INCAR.WAN.XXX)
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|ALGO}} = NONE                            # no electronic changes required
{{TAGBL|NELM}} = 1                              # one electronic step suffices, since WAVECAR from previous step is present
{{TAGBL|NBANDS}} = 48                            # number of bands used
{{TAGBL|LWAVE}} = .FALSE.                        # do not overwrite WAVECAR
{{TAGBL|LCHARG}} = .FALSE.                      # do not overwrite CHGCAR
{{TAGBL|LWANNIER90_RUN}} = .TRUE.                # run wannier90 in library mode
# As of vasp-6.4.3 define wannier90 interpolation from INCAR as follows:
# For older versions copy wannier90.win.XXX to wannier90.win and use INCAR.WAN as input for vasp
{{TAGBL|NUM_WANN}} = 3
{{TAGBL|WANNIER90_WIN}} = "
# use this energy window for t2g states for GW
dis_win_min = 7.4
dis_win_max = 9.95
begin projections
V:dxy;dxz;dyz
end projections
# plot bands
bands_plot = true
begin kpoint_path
R  0.50000000  0.50000000  0.50000000  G  0.00000000  0.00000000  0.00000000
G  0.00000000  0.00000000  0.00000000  X  0.50000000  0.00000000  0.00000000
X  0.50000000  0.00000000  0.00000000  M  0.50000000  0.50000000  0.00000000
  M  0.50000000  0.50000000  0.00000000  G  0.00000000  0.00000000  0.00000000
end kpoint_path
"


== The GW calculation ==
{{NB|mind|Prior VASP.6.4.3 a proper wannier90.win.XXX file should be used instead.}}
Use the corresponding INCAR.WAN.XXX file as input for <tt>wannier90</tt>
cp INCAR.WAN.XXX INCAR
where XXX=PBE, GW0 or HSE.


*INCAR.GW0
Use the corresponding WAVECAR.XXX file as input
<pre>
cp WAVECAR.XXX {{TAGBL|WAVECAR}}
System = SrVO3
and restart VASP. If all went well, the Vanadium t2g band dispersion thus obtained, may conveniently be visualized with gnuplot:
  gnuplot -persist ./wannier90_band.gnu


ISMEAR = -5
:'''N.B.:''' Most modern versions of <tt>gnuplot</tt> will respond with an error message unless you remove the first line of <tt>wannier90_band.gnu</tt> (some deprecated syntax issue).
EMIN = -20 ; EMAX = 20 ; NEDOS = 1000  # usefull energy range for density of states


NBANDS = 96                            # need for a lot of bands in GW
=== The preferred way to calculate the PBE bandstructure ===


ALGO = GW0                            #
Provided one has a self-consistent charge density ({{FILE|CHGCAR}}) file of sufficient quality (generated using a regular grid of k-points of sufficient density) one may read this charge density and keep it fixed ({{TAG|ICHARG}}=11).
NELM = 1                              # one step so this is really G0W0
For density functional calculations this charge density completely defines the Hamiltonian and using this Hamiltonian one may non-selfconsistently determine the orbitals and corresponding eigenenergies at arbitrary k-points.
PRECFOCK = Fast                        # select fast mode for FFT's
This is a very convenient way to calculate the bandstructure.
ENCUTGW = 100                          # energy cutoff for response function
NOMEGA = 200                          # metal, we need a lot of frequency points
MAXMEM = 2500                          # memory per core
NKRED = 2                              # sample down the GW to a coarse 2x2x2 grid


KPAR = 3
First we copy the self-consistent charge density of one of our previous calculations:
</pre>
 
  cp CHGCAR.PBE CHGCAR
  cp WAVECAR.PBE WAVECAR


Copy the aforementioned file to {{TAG|INCAR}}:
The bandstructure is conventionally plotted along lines of high symmetry in the 1st Brillouin zone.
The easiest way to specify these is by means of the so-called ''linemode'':
*{{TAG|KPOINTS}} (see KPOINTS.BSTR)
Auto
15
Linemode
reciprocal
0.50000000  0.50000000  0.50000000  !R
0.00000000  0.00000000  0.00000000  !G
0.00000000  0.00000000  0.00000000  !G
0.50000000  0.00000000  0.00000000  !X
0.50000000  0.00000000  0.00000000  !X
0.50000000  0.50000000  0.00000000  !M
0.50000000  0.50000000  0.00000000  !M
0.00000000  0.00000000  0.00000000  !G


cp INCAR.GW0 INCAR
'''N.B.:''' using these ''k''-points for a ''self-consistent'' calculation ({{TAG|ICHARG}}<10) would be a very bad idea since such and irregular sampling of the 1st Brillouin zone will not yield sensible charge densities.


and restart VASP.
Use the following {{TAGBL|INCAR}} file:
*{{TAG|INCAR}} (see INCAR.BSTR)
{{TAGBL|SYSTEM}} = SrVO3                        # system name
{{TAGBL|ISMEAR}} = 0                            # Gaussian smearing
{{TAGBL|EDIFF}} = 1E-7                          # tight convergence criterion
{{TAGBL|NBANDS}} = 36                            # 36 bands are sufficient
{{TAGBL|LWAVE}} = .FALSE.                        # do not overwrite WAVECAR
{{TAGBL|LCHARG}} = .FALSE.                      # do not overwrite CHGCAR
{{TAGBL|ICHARG}} = 11                            # read the charge density from the CHGCAR file and keep it fixed
{{TAGBL|LORBIT}} = 11                            # compute lm-decomposed states
{{TAGBL|EMIN}} = -20 ; {{TAGBL|EMAX}} = 20                # smallest/largest energy included in calculation
{{TAGBL|NEDOS}} = 1000                          # sampling points for DOS


=== Analysis of the DOS ===
'''N.B.:''' Mind that this approach works only for density functional calculations (''e.g.'' PBE or LDA) and is not applicable to orbital dependent functionals (like hybrid functionals) or in case of GW calculations.


Again, add the following line to your {{FILE|INCAR}} file:
This PBE bandstructure and the Wannier-interpolated structures of the PBE, HSE and GW calculation can be compared via
./plotbands.sh


LORBIT = 11
[[File:SrVO3 bandstructure.png|600px]]


and rerun VASP.
== Download ==
 
[[Media:SrVO3 GW band.zip| SrVO3_GW_band.zip]]
To plot the total DOS and the Vanadium ''t<sub>2g</sub>'' and ''e<sub>g</sub>'' partial-DOS using ''gnuplot'', execute the following command:
 
./plotdos
 
'''Mind''': Check the {{FILE|OUTCAR}} file for the position of the Fermi level. These DOSs have not been shifted such that the Fermi level is at 0 eV.
 
=== Bandstructure using WANNIER90 ===
Again, add the following line to your {{FILE|INCAR}} to have VASP call [http://www.wannier.org WANNIER90]:
<pre>
LWANNIER90_RUN = .TRUE.
</pre>
 
and use the following [http://www.wannier.org WANNIER90] input:
 
*wannier90.win.gw
<pre>
bands_plot = true
 
begin kpoint_path
R  0.50000000  0.50000000  0.50000000  G  0.00000000  0.00000000  0.00000000
G  0.00000000  0.00000000  0.00000000  X  0.50000000  0.00000000  0.00000000
X  0.50000000  0.00000000  0.00000000  M  0.50000000  0.50000000  0.00000000
M  0.50000000  0.50000000  0.00000000  G  0.00000000  0.00000000  0.00000000
end kpoint_path
 
num_wann =    3
 
num_bands=    3
 
exclude_bands : 1-20, 24-96
 
begin projections
V:dxy;dxz;dyz
end projections
</pre>
 
Copy the above to {{FILE|wannier90.win}}:
 
cp wannier90.win.gw wannier90.win
 
and redo the GW step.
 
To compare the Vanadium ''t<sub>2g</sub>'' band dispersion in the GW approximation with the LDA bandstructure, run the following command:
 
gnuplot -persist plotme.gw
 
'''Mind''': Here the eigenvalues have been shifted such that the Fermi level is a 0 eV.
 
=== A more accurate GW calculation ===
 
As you might have noticed in the previous example the Vanadium ''t<sub>2g</sub>'' bands look a bit ''wobbly'' along G-X and X-M.
In the present example, this turns out to be an artifact of the downsampling of the GW.
Try removing (or comment out) the line
 
NKRED = 2
 
from the {{FILE|INCAR}} file, and redo the GW step.
 
== A comparison to the HSE hybrid functional ==
 
To illustrate the kind of results one would obtain for SrVO<sub>3</sub> using the DFT/Hartree-Fock hybrid functional HSE, without actually doing a full selfconsistent calculation, we will recalculate the one-electron energies and DOS ({{TAG|ALGO}}=Eigenval) using the HSE functional with DFT orbitals as input:
 
*INCAR.HSE
<pre>
System  = SrVO3
 
ISMEAR = -5
EMIN = -20 ; EMAX = 20 ; NEDOS = 1000  # usefull energy range for density of states
 
EDIFF = 1E-8                          # high precision for groundstate calculation
 
KPAR = 3
LHFCALC = .TRUE.  ; HFSCREEN = 0.2  ; NBANDS = 48
PRECFOCK = Fast  ; NELM = 1
ALGO = Eigenval
 
LWAVE = .FALSE.                        # do not write the wave functions
</pre>
 
Copy the aforementioned file to {{TAG|INCAR}}:
 
cp INCAR.HSE INCAR
 
and restart VASP.
 
'''Mind''': This calculation (and the ones following below) needs to restart from a set of converged DFT wave functions, therefore:
 
cp WAVECAR.DFT.96bands WAVECAR
 
=== Analysis of the DOS ===
 
Again, add the following line to your {{FILE|INCAR}} file:
 
LORBIT = 11
 
and rerun VASP.
 
To plot the total DOS and the Vanadium ''t<sub>2g</sub>'' and ''e<sub>g</sub>'' partial-DOS using ''gnuplot'', execute the following command:
 
./plotdos
 
'''Mind''': Check the {{FILE|OUTCAR}} file for the position of the Fermi level. These DOSs have not been shifted such that the Fermi level is at 0 eV.
 
=== Bandstructure using WANNIER90 ===
 
As before, add the following line to your {{FILE|INCAR}} to have VASP call [http://www.wannier.org WANNIER90]:
<pre>
LWANNIER90_RUN = .TRUE.
</pre>
 
and use the following [http://www.wannier.org WANNIER90] input:
 
*wannier90.win.hse
<pre>
bands_plot = true
 
begin kpoint_path
R  0.50000000  0.50000000  0.50000000  G  0.00000000  0.00000000  0.00000000
G  0.00000000  0.00000000  0.00000000  X  0.50000000  0.00000000  0.00000000
X  0.50000000  0.00000000  0.00000000  M  0.50000000  0.50000000  0.00000000
M  0.50000000  0.50000000  0.00000000  G  0.00000000  0.00000000  0.00000000
end kpoint_path
 
num_wann =    3
 
num_bands=    3
 
exclude_bands : 1-20, 24-48


begin projections
{{Template:GW - Tutorial}}
V:dxy;dxz;dyz
end projections
</pre>


Copy the above to {{FILE|wannier90.win}}:
Back to the [[The_VASP_Manual|main page]].
 
cp wannier90.win.hse wannier90.win
 
and redo the HSE calculation.
 
To compare the HSE Vanadium ''t<sub>2g</sub>'' band dispersion with the LDA bandstructure, run the following command:
 
gnuplot -persist plotme.hse
 
'''Mind''': Here the eigenvalues have been shifted such that the Fermi level is a 0 eV.
 
== Download ==
[http://www.vasp.at/vasp-workshop/examples/SrVO3.tgz SrVO3.tgz]
 
----
[[VASP_example_calculations|To the list of examples]] or to the [[The_VASP_Manual|main page]]


[[Category:Examples]]
[[Category:Examples]]

Latest revision as of 10:44, 12 June 2024

Task

Calculation of the GW bandstructure of SrVO3 using VASP and WANNIER90.

Performing a GW calculation with VASP is a 3-step procedure: a DFT groundstate calculation, a calculation to obtain a number of virtual orbitals, and the actual GW calculation itself. In this example we will also see how the results of the GW calculation may be postprocessed with WANNIER90 to obtain the dispersion of the bands along the usual high symmetry directions in reciprocal space.

N.B.: This example involves quite a number of individual calculations. The easiest way to run this example is to execute: 

./doall.sh

And compare the output of the different steps (DFT, GW, HSE) by: 

./plotall.sh

In any case, one can consider the doall.sh script to be an overview of the steps described below.

DFT groundstate calculation

The first step is a conventional DFT (in this case PBE) groundstate calculation. 

SYSTEM  = SrVO3                        # system name
NBANDS = 36                            # small number  of bands
ISMEAR = 0                             # Gaussian smearing
EDIFF = 1E-8                           # high precision for groundstate calculation
KPAR = 2                               # parallelization of k-points in two groups

Copy the aforementioned file to INCAR: 

cp INCAR.DFT INCAR

The POSCAR file describes the structure of the system: 

SrVO3
3.84652  #cubic fit for 6x6x6 k-points
 +1.0000000000  +0.0000000000  +0.0000000000 
 +0.0000000000  +1.0000000000  +0.0000000000 
 +0.0000000000  +0.0000000000  +1.0000000000 
Sr V O
 1 1 3
Direct
 +0.0000000000  +0.0000000000  +0.0000000000 
 +0.5000000000  +0.5000000000  +0.5000000000 
 +0.5000000000  +0.5000000000  +0.0000000000 
 +0.5000000000  +0.0000000000  +0.5000000000 
 +0.0000000000  +0.5000000000  +0.5000000000

The POSCAR can be visualized with p4v or VESTA and remains unchanged in the following.

The KPOINTS file describes how the first Brillouin zone is sampled. In the first step we use a uniform k-point sampling: 

Automatically generated mesh
       0
Gamma
 4 4 4
 0 0 0

Mind: this is definitely not dense enough for a high-quality description of SrVO3, but in the interest of speed we will live with it. Copy the aforementioned file to KPOINTS: 

cp KPOINTS.BULK KPOINTS

and run VASP. If all went well, one should obtain a WAVECAR file containing the PBE wavefunction.

Obtain DFT virtual orbitals and long-wave limit

Use following INCAR file to increase the number of virtual states and to determine the long-wave limit of the polarizability (stored in WAVEDER): 

SYSTEM = SrVO3                         # system name
ISMEAR = 0                             # Gaussian smearing
KPAR = 2                               # parallelization of k-points in two groups
ALGO = Exact                           # exact diagonalization
NELM = 1                               # one electronic step suffices, since WAVECAR from previous step is present
NBANDS = 96                            # need for a lot of bands in GW
LOPTICS = .TRUE.                       # we need d phi/ d k  for GW calculations for long-wave limit

Restart VASP. At this stage it is a good idea to make a safety copy of the WAVECAR and WAVEDER files since we will repeatedly need them in the calculations that follow: 

cp WAVECAR WAVECAR.PBE
cp WAVEDER WAVEDER.PBE

Also make a backup of the charge density for later: 

cp CHGCAR CHGCAR.PBE

The dielectric function

As a bonus, VASP determines the frequency dependent dielectric function in the independent-particle (IP) picture and writes the result to the OUTCAR and vasprun.xml files. In the OUTCAR you should search for 

 frequency dependent IMAGINARY DIELECTRIC FUNCTION (independent particle, no local field effects)

and 

 frequency dependent      REAL DIELECTRIC FUNCTION (independent particle, no local field effects)

GW Step

The actual GW calculation requires a set of one-electron energies and eigenstates. In this case we use the PBE solution obtained from previous step: 

cp WAVECAR.PBE WAVECAR
cp WAVEDER.PBE WAVEDER

The following INCAR file selects the 'single shot' GW calculation also known as G0W0: 

SYSTEM = SrVO3                         # system name
ISMEAR = 0                             # Gaussian smearing
KPAR = 2                               # parallelization of k-points in two groups
ALGO = GW0                             # GW with iteration in G, W kept on DFT level
NELM = 1                               # one electronic step suffices, since WAVECAR from previous step is present
NBANDS = 96                            # need for a lot of bands in GW
PRECFOCK = Fast                        # fast mode for FFTs
ENCUTGW = 100                          # small energy cutoff for response function suffices for this tutorial
NOMEGA = 200                           # large number of real frequency points for Hilbert transforms of W and self-energy

Restarting VASP will overwrite the present WAVECAR and vasprun.xml file. Make a copy them for later. 

cp WAVECAR WAVECAR.GW0
cp vasprun.xml vasprun.GW0.xml

HSE hybrid functional

To illustrate the kind of results one would obtain for SrVO3 using the DFT/Hartree-Fock hybrid functional HSE, without actually doing a full self-consistent calculation, we will recalculate the one-electron energies and DOS (ALGO=Eigenval) using the HSE functional with DFT orbitals as input 

cp WAVECAR.PBE WAVECAR

Use the following INCAR file: 

SYSTEM = SrVO3                         # system name
ISMEAR = 0                             # Gaussian smearing
KPAR = 2                               # parallelization of k-points in two groups
ALGO = Eigenval                        # calulate eigenvalues
NELM = 1                               # one electronic step suffices, since WAVECAR from previous step is present
NBANDS = 48                            # small number of bands suffice
PRECFOCK = Fast                        # fast mode for FFTs
LHFCALC = .TRUE.                       # switch on Hartree-Fock routines to calculate exact exchange
HFSCREEN = 0.2                         # HSE06 screening parameter

Restart VASP and make a copy of the wavefunction for post-processing 

cp WAVECAR WAVECAR.HSE

Post-processing: Density of states and Bandstructure for PBE, GW and HSE

Density of States

The DOS of the PBE, GW and HSE solution can be calculated in a post-processing step with 

SYSTEM = SrVO3                         # system name
ISMEAR = -5                            # Bloechl's tetrahedron method (requires at least 3x3x3 k-points)
ALGO = NONE                            # no electronic changes required
NELM = 1                               # one electronic step suffices, since WAVECAR from previous step is present
NBANDS = 48                            # number of bands used
EMIN = -20 ; EMAX = 20                 # smallest/largest energy included in calculation
NEDOS = 1000                           # sampling points for DOS
LORBIT = 11                            # calculate l-m decomposed DOS
LWAVE = .FALSE.                        # do not overwrite WAVECAR
LCHARG = .FALSE.                       # do not overwrite CHGCAR

and requires the apropriate WAVECAR file from one of the previous steps. Copy 

cp WAVECAR.PBE WAVECAR

or 

cp WAVECAR.GW0 WAVECAR

or 

cp WAVECAR.HSE WAVECAR

and restart VASP. The density of states is written to DOSCAR, make a copy of this file 

cp DOSCAR DOSCAR.XXX

where XXX is either PBE, GW0 or HSE. Visualize the projected DOS for the V-t2g, V-eg and O-p states with the scriptfile 

./plotdos_comp.sh DOSCAR.???

This requires gnuplot to be installed.

Bandstructure with wannier90

The bandstructure can be calculated via Wannier interpolation using wannier90 in the library mode

  • INCAR (see INCAR.WAN.XXX)
SYSTEM = SrVO3                         # system name
ISMEAR = 0                             # Gaussian smearing
ALGO = NONE                            # no electronic changes required
NELM = 1                               # one electronic step suffices, since WAVECAR from previous step is present
NBANDS = 48                            # number of bands used
LWAVE = .FALSE.                        # do not overwrite WAVECAR
LCHARG = .FALSE.                       # do not overwrite CHGCAR
LWANNIER90_RUN = .TRUE.                # run wannier90 in library mode

# As of vasp-6.4.3 define wannier90 interpolation from INCAR as follows:
# For older versions copy wannier90.win.XXX to wannier90.win and use INCAR.WAN as input for vasp
NUM_WANN = 3
WANNIER90_WIN = "

# use this energy window for t2g states for GW
dis_win_min = 7.4
dis_win_max = 9.95

begin projections
V:dxy;dxz;dyz
end projections

# plot bands
bands_plot = true

begin kpoint_path
R  0.50000000  0.50000000  0.50000000  G  0.00000000  0.00000000  0.00000000
G  0.00000000  0.00000000  0.00000000  X  0.50000000  0.00000000  0.00000000
X  0.50000000  0.00000000  0.00000000  M  0.50000000  0.50000000  0.00000000
M  0.50000000  0.50000000  0.00000000  G  0.00000000  0.00000000  0.00000000
end kpoint_path
"


Mind: Prior VASP.6.4.3 a proper wannier90.win.XXX file should be used instead.

Use the corresponding INCAR.WAN.XXX file as input for wannier90 

cp INCAR.WAN.XXX INCAR

where XXX=PBE, GW0 or HSE.

Use the corresponding WAVECAR.XXX file as input 

cp WAVECAR.XXX WAVECAR

and restart VASP. If all went well, the Vanadium t2g band dispersion thus obtained, may conveniently be visualized with gnuplot: 

gnuplot -persist ./wannier90_band.gnu
N.B.: Most modern versions of gnuplot will respond with an error message unless you remove the first line of wannier90_band.gnu (some deprecated syntax issue).

The preferred way to calculate the PBE bandstructure

Provided one has a self-consistent charge density (CHGCAR) file of sufficient quality (generated using a regular grid of k-points of sufficient density) one may read this charge density and keep it fixed (ICHARG=11). For density functional calculations this charge density completely defines the Hamiltonian and using this Hamiltonian one may non-selfconsistently determine the orbitals and corresponding eigenenergies at arbitrary k-points. This is a very convenient way to calculate the bandstructure.

First we copy the self-consistent charge density of one of our previous calculations: 

 cp CHGCAR.PBE CHGCAR
 cp WAVECAR.PBE WAVECAR

The bandstructure is conventionally plotted along lines of high symmetry in the 1st Brillouin zone. The easiest way to specify these is by means of the so-called linemode: 

Auto
15
Linemode
reciprocal
0.50000000  0.50000000  0.50000000   !R
0.00000000  0.00000000  0.00000000   !G

0.00000000  0.00000000  0.00000000   !G
0.50000000  0.00000000  0.00000000   !X

0.50000000  0.00000000  0.00000000   !X
0.50000000  0.50000000  0.00000000   !M 

0.50000000  0.50000000  0.00000000   !M
0.00000000  0.00000000  0.00000000   !G

N.B.: using these k-points for a self-consistent calculation (ICHARG<10) would be a very bad idea since such and irregular sampling of the 1st Brillouin zone will not yield sensible charge densities.

Use the following INCAR file: 

SYSTEM = SrVO3                         # system name
ISMEAR = 0                             # Gaussian smearing
EDIFF = 1E-7                           # tight convergence criterion
NBANDS = 36                            # 36 bands are sufficient
LWAVE = .FALSE.                        # do not overwrite WAVECAR
LCHARG = .FALSE.                       # do not overwrite CHGCAR
ICHARG = 11                            # read the charge density from the CHGCAR file and keep it fixed
LORBIT = 11                            # compute lm-decomposed states
EMIN = -20 ; EMAX = 20                 # smallest/largest energy included in calculation
NEDOS = 1000                           # sampling points for DOS

N.B.: Mind that this approach works only for density functional calculations (e.g. PBE or LDA) and is not applicable to orbital dependent functionals (like hybrid functionals) or in case of GW calculations.

This PBE bandstructure and the Wannier-interpolated structures of the PBE, HSE and GW calculation can be compared via 

./plotbands.sh

Download

SrVO3_GW_band.zip

Back to the main page.

Retrieved from "https://www.vasp.at/wiki/index.php?title=Bandstructure_of_SrVO3_in_GW&oldid=25105"