Tutorial:DFT+U

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.

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