Hubbard U Values

Details on Hubbard U corrections used by the Materials Project

Hubbard U Values

It is well-known that first principles calculations within the local density approximation (LDA) or generalized gradient approximation (GGA) lead to considerable error in calculated redox reaction energies of many transition metal compounds. This error arises from the self-interaction error in LDA and GGA, which is not canceled out in redox reactions where an electron is transferred between significantly different environments, such as between a metal and a transition metal or between a transition metal and oxygen or fluorine. Extensive discussion of this issue can be found in the following works. [1-4]

In the Materials Project, we have calibrated UU values for many transition metals of interest using the approach outlined in Wang et al.'s work [5]. At the present moment, UU values have only been calibrated for transition metal oxide systems. UU values were calibrated for the following elements: Co\text{Co}, Cr\text{Cr}, Fe\text{Fe}, Mn\text{Mn}, Mo\text{Mo}, Ni\text{Ni}, V\text{V} and W\text{W}. The choice of systems to which we apply UU was largely determined by our experience and by systematic benchmarking. It is very likely that we will expand calibration of UU values to more chemical systems in the future.

In the Materials Project, for an oxide or fluoride material with a transition element listed previously, with the VASP input settings constructed according to the logic defined in pymatgen.

Note that for fluorides, the UU value gets set to the one calibrated from the oxide system, although in principle our architecture allows different UU values to be set for oxides and fluorides respectively.

Calibration and Values

The UU values were obtained by fitting to experimental binary formation enthalpies as described in Wang et al.'s work. This method is simple yet accurately reproduces phase stabilities. A least squares method of obtaining the correct UU value was used, as follows:

  1. For each non-overlapping formation energy reaction considered, we find the region where the formation energy error passes zero. For the V-O\text{V-O} system, this includes the following:

    • 2V2O3+O24VO22\text{V}_2\text{O}_3 + \text{O}_2 \rightarrow 4 \text{VO}_2

    • 4VO2+O22V2O54 \text{VO}_2 + \text{O}_2 \rightarrow 2\text{V}_2\text{O}_5

  2. For each formation energy region identified, we fit the linear equation \begin{align} \mbox{Error/redox} & = m U + c \end{align} to the final UU range. In the case of V\text{V}, we will have two sets of (m,c)(m,c).

  3. We find the U value that minimizes the sum of square Error / Redox.

  4. In the case of V\text{V}, we get a UU value of 3.25.

The full list of U values used is described in the table below. For oxides and fluorides containing any of the elements, only GGA+U calculations are performed.

ElementSystemFitting ReactionRedox CoupleCalibrated U (eV)Comments


























The U values are calibrated for phase stability analyses, and should be used with care if applied to obtain other properties such as band structures. Also, the U values depend on the pseudopotential used. Further, typically, U values should be site specific, however in our approach, U values were applied to all sites with an element listed above, and only to the d-orbitals. A discussion of the pseudopotentials used in the Materials Project can be found here.


[1]: F. Zhou, M. Cococcioni, C. A. Marianetti, D. Morgan and G. Ceder. First-principles prediction of redox potentials in transition-metal compounds with LDA+U. Physical Review B, 2004, 70, 235121. doi:10.1103/PhysRevB.70.235121

[2]: M. Cococcioni, S. de Gironcoli, Linear response approach to the calculation of the effective interaction parameters in the LDA+U method. Physical Review B, 2005, 71, 035105. doi:10.1103/PhysRevB.71.035105

[3]: L. Wang, T. Maxisch, & G. Ceder. Oxidation energies of transition metal oxides within the GGA+U framework. Physical Review B. 2006, 73, 195107, doi:10.1103/PhysRevB.73.195107

[4]: A. Jain, G. Hautier, S. P. Ong, C. Moore, C. Fischer, K. A. Persson, & G. Ceder. Formation enthalpies by mixing GGA and GGA + U calculations. Physical Review B, 2011, 84(4), 045115. doi:10.1103/PhysRevB.84.045115

[5]: M. Wang, A. Navrotsky Enthalpy of formation of LiNiO2, LiCoO2 and their solid solution, LiNi1-xCoxO2, Solid State Ionics, vol. 166, no. 1-2, pp. 167-173, Jan. 2004.

Last updated