Tutorial:DFT+U

From OctopusWiki
Jump to navigation Jump to search

The objective of this tutorial is to give a basic idea of how DFT+U in Octopus works.

Input

As a prototypical example for DFT+U, we will consider bulk NiO in its anti-ferromagnetic configuration. We will neglect the small lattice distortion and consider only its cubic cell.

 CalculationMode = gs
 PeriodicDimensions = 3 
 BoxShape = parallelepiped
 ExperimentalFeatures = yes
 PseudopotentialSet=hscv_pbe
 a = 7.8809
 %LatticeParameters
   a | a | a
 %
 %LatticeVectors
  0.0 | 1/2 | 1/2
  1/2 | 0.0 | 1/2
  1.0 | 1.0 | 0.0
 %
 
 %Species
 "Ni" | species_pseudo | hubbard_l | 2 | hubbard_u | 5.0*eV
 %
 
 DFTULevel = dft_u_empirical
 
 %ReducedCoordinates
  "Ni" | 0.0 | 0.0 | 0.0
  "Ni" | 0.0 | 0.0 | 0.5
  "O"  | 0.5 | 0.5 | 0.25
  "O"  | 0.5 | 0.5 | 0.75
 %
 Spacing = 0.5
 %KPointsGrid
 2 | 2 | 2
 %
 ParDomains = no
 ParKPoints = auto
 
 SpinComponents = polarized
 GuessMagnetDensity = user_defined
 %AtomsMagnetDirection
  8.0
 -8.0
  0.0
  0.0
 %
 
 OutputLDA_U = occ_matrices

As we are interested by the antiferromagnetic order, the primitive cell is doubled along the last lattice vector. To help the convergence, an initial guess should be added, by adding to the input file the variables GuessMagnetDensity and AtomsMagnetDirection.

In order to perform a calculation with a U of 5eV on the 3d orbitals (corresponding to the quantum number l=2), we define a block Species, where hubbard_l specifies the orbitals (l=0 for s orbitals, l=1 for p orbitals, ...) and hubbard_u is used to set the value of the effective Hubbard U.

In order to activate the DFT+U part, one finally needs to specify the level of DFT+U used . This is done using the variable DFTULevel. At the moment there is three possible options for this variable, which correspond to no +U correction (dft_u_none), an empirical correction (dft_u_empirical) or the ab initio U correction based on the ACBN0 functional[1] (dft_u_acbn0).

Some specific outputs can then be added, such as the density matrix of the selected localized subspaces.

Output

References

  1. Agapito, Luis A. and Curtarolo, Stefano and Buongiorno Nardelli, Marco, Reformulation of $\mathrm{DFT}+U$ as a Pseudohybrid Hubbard Density Functional for Accelerated Materials Discovery, Phys. Rev. X 5 011006 (2015)

Back to Tutorials