Materials Project Documentation
Return to materialsproject.org
  • Introduction
  • Frequently Asked Questions (FAQ)
    • Glossary of Terms
  • Changes and Updates
    • Database Versions
    • Website Changelog
  • Documentation Credit
  • Community
    • Getting Help
    • Getting Involved
      • Contributor Guide
      • Potential Collaborators
      • MP Community Software Ecosystem
    • Community Resources
    • Code of Conduct
  • Services
    • MPContribs
  • Methodology
    • Materials Methodology
      • Overview
      • Calculation Details
        • GGA/GGA+U Calculations
          • Parameters and Convergence
          • Hubbard U Values
          • Pseudo-potentials
        • r2SCAN Calculations
          • Parameters and Convergence
          • Pseudopotentials
      • Thermodynamic Stability
        • Energy Corrections
          • Anion and GGA/GGA+U Mixing
          • GGA/GGA+U/r2SCAN Mixing
        • Phase Diagrams (PDs)
        • Chemical Potential Diagrams (CPDs)
        • Finite Temperature Estimation
      • Electronic Structure
      • Phonon Dispersion
      • Diffraction Patterns
      • Aqueous Stability (Pourbaix)
      • Magnetic Properties
      • Elastic Constants
      • Piezoelectric Constants
      • Dielectric Constants
      • Equations of State (EOS)
      • X-ray Absorption Spectra (XAS)
      • Surface Energies
      • Grain Boundaries
      • Charge Density
      • Suggested Substrates
      • Related Materials
      • Optical absorption spectra
      • Alloys
    • Molecules Methodology
      • Overview
      • Calculation Details
      • Atomic Partial Charges
      • Atomic Partial Spins
      • Bonding
      • Metal Coordination and Binding
      • Natural Atomic and Molecular Orbitals
      • Redox and Electrochemical Properties
      • Molecular Thermodynamics
      • Vibrational Properties
      • Legacy Data
    • MOF Methodology
      • Calculation Parameters
        • DFT Parameters
        • Density Functionals
        • Pseudopotentials
        • DFT Workflow
  • Apps
    • Explore and Search Apps
      • Materials Explorer
        • Tutorial
      • Molecules Explorer
        • Tutorial
        • Legacy Data
      • Battery Explorer
        • Background
        • Tutorial
      • Synthesis Explorer
        • Background
        • Tutorial
      • Catalysis Explorer
        • Tutorial
      • MOF Explorer
        • Downloading the Data
        • Structure Details
          • QMOF IDs
          • Structure Sources
          • Finding MOFs by Common Name
          • Structural Fidelity
        • Property Definitions
          • SMILES, MOFid, and MOFkey
          • Pore Geometry
          • Topology
          • Electronic Structure
          • Population Analyses and Bond Orders
          • Symmetry
        • Version History
        • How to Cite
    • Analysis Apps
      • Phase Diagram
        • Background
        • Tutorials
        • FAQ
      • Pourbaix Diagram
        • Background
        • Tutorial
        • FAQ
      • Crystal Toolkit
        • Background
        • Tutorial
        • FAQ
      • Reaction Calculator
      • Interface Reactions
    • Characterization Apps
      • X-ray Absorption Spectra (XAS)
    • Explore Contributed Data
  • Downloading Data
    • How do I download the Materials Project database?
    • Using the API
      • Getting Started
      • Querying Data
      • Tips for Large Downloads
      • Examples
      • Advanced Usage
    • Differences between new and legacy API
    • Query and Download Contributed Data
    • AWS OpenData
  • Uploading Data
    • Contribute Data
  • Data Production
    • Data Workflows
    • Data Builders
Powered by GitBook
On this page
  • Methodology
  • 1) Anion corrections
  • 2) GGA / GGA+U Mixing Corrections
  • Accuracy of Total Energies
  • Estimating errors in calculated reaction energies
  • Sources of error
  • GGA errors on reaction energies between chemically similar compounds
  • Citations
  • References
  • Acknowledgments

Was this helpful?

Edit on GitHub
Export as PDF
  1. Methodology
  2. Materials Methodology
  3. Thermodynamic Stability
  4. Energy Corrections

Anion and GGA/GGA+U Mixing

Details on anion and GGA/GGA+U mixing scheme corrections

PreviousEnergy CorrectionsNextGGA/GGA+U/r2SCAN Mixing

Last updated 2 years ago

Was this helpful?

This correction scheme assumes independent, linear corrections associated with each corrected element. For example, VO2\text{VO}_2VO2​ would receive both a 'V\text{V}V' and an 'oxide' correction (as explained below), while elemental V\text{V}V would receive no corrections. For complete details of our correction scheme, refer to Wang et al.

Methodology

1) Anion corrections

For many elements that take on negative oxidation states in solids, differences in electron localization between the elements and the solid can result in substantial errors in formation energies computed from DFT calculations. This is especially true for elements that are gaseous in their standard state - O2\text{O}_2O2​, N2\text{N}_2N2​, Cl2\text{Cl}_2Cl2​, F2\text{F}_2F2​, and H2\text{H}_2H2​.

To address this, we adjust the energies of materials containing certain elements by applying a correction to anionic species, as explained in ref . Specifically, we apply energy corrections to 14 anion species -- 'oxide', 'peroxide', 'superoxide', S\text{S}S, F\text{F}F, Cl\text{Cl}Cl, Br\text{Br}Br, I\text{I}I, N\text{N}N, H\text{H}H, Se\text{Se}Se, Si\text{Si}Si, Sb\text{Sb}Sb, and Te\text{Te}Te. In the case of oxygen-containing compounds, separate corrections are applied to oxides, superoxides, and peroxides based on the specific bonding environment of oxygen in the material, as determined from nearest-neighbor bond lengths (e.g., <1.35 Å for superoxide, <1.49 Å for 'peroxide', and 'oxide' otherwise). Thus, Na2O\text{Na}_2\text{O}Na2​O receives an 'oxide' correction while NaO2\text{NaO}_2NaO2​ receives a `superoxide' correction.

Anion corrections are applied to a material only when it contains a corrected element as an anion. For example, the 'H\text{H}H' correction is applied to LiH\text{LiH}LiH but not to H2O\text{H}_2\text{O}H2​O. A specie is classified as an anion if its estimated oxidation state (when available) is negative, or if it is the most electronegative element in the formula.

2) GGA / GGA+U Mixing Corrections

Some compounds are better modeled with a U correction term to the density functional theory Hamiltonian while others are better modeled without (i.e., regular GGA). Energies from calculations with the +U correction are not directly comparable to those without. To obtain better accuracy across chemical systems, we use GGA+U when appropriate, GGA otherwise, and mix energies from the two calculation methodologies by adding an energy correction term to the GGA+U calculations to make them comparable to the GGA calculations.

Specifically, we use GGA+U for oxide and fluoride compounds containing any of the transition metals V\text{V}V, Cr\text{Cr}Cr, Mn\text{Mn}Mn, Fe\text{Fe}Fe, Co\text{Co}Co, Ni\text{Ni}Ni, W\text{W}W, and Mo\text{Mo}Mo, and GGA for everything else. More details on this method can be found in refs.

Accuracy of Total Energies

To estimate the accuracy of our total energy calculations, we compute reaction data and compare against experimental data. Note that this data set was compiled using a lower k-point mesh and pseudopotentials with fewer electrons than the current Materials Project parameter set.

Estimating errors in calculated reaction energies

The accuracy of calculated reaction energies depends on the chemical system investigated. In general, GGA calculations have similar errors among chemically similar systems. Hence, reaction energies between chemically similar systems (e.g., a reaction where the reactants and products are all oxides, such as MgO + Al2O3→MgAl2O4\text{MgO + Al2O3} \rightarrow\text{MgAl2O4}MgO + Al2O3→MgAl2O4 tend to have smaller errors than reactions between chemically dissimilar systems (e.g., between metals and insulators).

Figure 1: Errors in Calculated Formation Energies for 413 binaries in the Kubaschewski Tables. Energies are normalized to per mol atom.

To provide a quantitative indicator of the error we may expect from the reaction calculator, we have computed the reaction energies of 413 binaries in the Kubaschewski Tables formed with Group V, VI and VII anions. Figure 1 shows the errors in the calculated formation energies (compared to the experimental values) for these compounds. The mean absolute error (MAE) is around 14 kJ mol−1^{-1}−1. 75% of the calculated formation energies are within 20 kJ mol−1^{-1}−1. We also found that compounds of certain elements tend to have larger errors. For example, Bi\text{Bi}Bi, Co\text{Co}Co, Pb\text{Pb}Pb, Eu\text{Eu}Eu, U\text{U}U, Tl\text{Tl}Tl and W\text{W}W compounds often have errors larger than 20 kJ mol−1^{-1}−1.

Sources of error

GGA errors on reaction energies between chemically similar compounds

The main conclusions are:

  • The error in reaction energies for the binary oxide to ternary oxides reaction energies are an order of magnitude lower than for the more often reported formation energies from the element. An error intrinsic to GGA (+U) is estimated to follow a normal distribution centered in zero (no systematic underestimation or overestimation) and with a standard deviation around 24 meV/at.

  • When looking at phase stability (and for instance assessing if a phase is stable or not), the relevant reaction energies are most of the time not the formation energies from the elements but reaction energies from chemically similar compounds (e.g., two oxides forming a third oxide). Large cancellation of errors explain this observation.

  • The +U is necessary for accurate description of the energetics even when reactions do not involve change in formal oxidation states

Citations

To cite the calculation methodology, please reference the following works:

References

[2]: A. Jain, G. Hautier, S.P. Ong, C. Moore, C.C. Fischer, K.A. Persson, G. Ceder, Formation Enthalpies by Mixing GGA and GGA+U calculations, Physical Review B, vol. 84 (2011), 045115.

[3]: J.B. Foresman, A.E. Frisch, Exploring Chemistry With Electronic Structure Methods: A Guide to Using Gaussian, Gaussian. (1996).

[4]: A. Jain, S.-a Seyed-Reihani, C.C. Fischer, D.J. Couling, G. Ceder, W.H. Green, Ab initio screening of metal sorbents for elemental mercury capture in syngas streams, Chemical Engineering Science. 65 (2010) 3025-3033.

[5]: S. Lany, Semiconductor thermochemistry in density functional calculations, Physical Review B. 78 (2008) 1-8.

[6]: G. Hautier, S.P. Ong, A. Jain, C. J. Moore, G. Ceder, Accuracy of density functional theory in predicting formation energies of ternary oxides from binary oxides and its implication on phase stability, Physical Review B, 85 (2012), 155208

Acknowledgments

Thank you to the original authors of this page:

  1. Anubhav Jain

  2. Shyue Ping Ong

  3. Geoffroy Hautier

  4. Charles Moore

It should be noted that while an MAE of 14 kJ mol−1^{-1}−1 is significantly higher than the desired chemical accuracy of 4 kJ mol−1^{-1}−1, it compares fairly well with the performance of most quantum chemistry calculations . Other than the most computationally expensive model chemistries such as G1-G3 and CBS, the reaction energy errors of most computational chemistry model chemistries are well above 10 kJ mol−1^{-1}−1.

For oxidation of the elements into binary compounds, an average error of ~4% or 33 kJ/mol-O2\text{O}_2O2​ is typical.[^9] For conventional ternary oxide formation from the elements, we have found a mean relative absolute error of about 2%.

The largest contribution to the error comes from the inability of the GGA to fully describe electronic exchange and correlation effects. In addition, there is some error associated with neglecting zero-point effects and with comparing 0K, 0atm computations with room-temperature enthalpy experiments. The latter effect was estimated to contribute less than 0.03 eV/atom by Lany. The stability of antiferromagnetic compounds may be underestimated, as the majority of our calculations are performed ferromagnetically only. The effect of magnetism may be small (under 10 meV/atom) or large (100 meV/atom or greater), depending on the compound. For compounds with heavy elements, relativistic effects may lead to greater-than-expected errors.

We recently conducted a more in-depth study comparing GGA (+U) reaction energies of ternary oxides from binary oxides on 135 compounds.

A. Jain, G. Hautier, C. Moore, S.P. Ong, C.C. Fischer, T. Mueller, K.A. Persson, G. Ceder., A High-Throughput Infrastructure for Density Functional Theory Calculations, Computational Materials Science, vol. 50, 2011, pp. 2295-2310.

A. Jain, G. Hautier, S.P. Ong, C. Moore, C.C. Fischer, K.A. Persson, G. Ceder, Accurate Formation Enthalpies by Mixing GGA and GGA+U calculations, Physical Review B, vol. 84, 2011, p. 045115.

[1]: Wang, A., Kingsbury, R.S., Horton, M., Jain, A., Ong, S.P., Dwaraknath, S., Persson, K. A framework for quantifying uncertainty in DFT energy corrections. Scientific Reports 11 (2021), 15496.

DOI:10.1016/j.commatsci.2011.02.023
DOI:10.1103/PhysRevB.84.045115
DOI: 10.1038/s41598-021-94550-5
[1]
[1]
[1,2]
[3]
[4]
[5]
[6]