Chemistry:Davidson correction

From HandWiki

The Davidson correction is an energy correction often applied in calculations using the method of truncated configuration interaction, which is one of several post-Hartree–Fock ab initio quantum chemistry methods in the field of computational chemistry. It was introduced by Ernest R. Davidson.[1] It allows one to estimate the value of the full configuration interaction energy from a limited configuration interaction expansion result, although more precisely it estimates the energy of configuration interaction up to quadruple excitations (CISDTQ) from the energy of configuration interaction up to double excitations (CISD). It uses the formula[2]

[math]\displaystyle{ \Delta E_Q = (1 - a_0^2)(E_{\rm CISD} - E_{\rm HF}), \ }[/math]
[math]\displaystyle{ E_{\rm CISDTQ} \approx E_{\rm CISD} + \Delta E_Q, \ }[/math]

where a0 is the coefficient of the Hartree–Fock wavefunction in the CISD expansion, ECISD and EHF are the energies of the CISD and Hartree–Fock wavefunctions respectively, and ΔEQ is the correction to estimate ECISDTQ, the energy of the CISDTQ wavefunction. Such estimation is based on perturbation theory analysis.[3] Therefore, CISD calculations with the Davidson correction included are frequently referred to as CISD(Q).

Application

The Davidson correction is very popular due to its low computational cost. The correction improves the contribution of electron correlation to the energy. The size-consistency and size-extensivity problems of truncated CI are alleviated but still exist. In small molecules, accuracy of the corrected energies can be similar to results from coupled cluster theory calculations.

The Davidson correction does not give information about the wave function. Therefore, it cannot be used to correct wave-function-dependent quantities such as dipole moment, charge density and vibronic couplings. Analytical gradients for Davidson corrections are in general not available in quantum chemistry programs.

As with other perturbative approaches, the Davidson correction is not reliable when the electronic structure of CISD and the reference Hartree–Fock wave functions are significantly different (i.e. when [math]\displaystyle{ a_0^2 }[/math] is not close to 1). This happens when multi-reference character is significant or when CISD is used to calculate a state that is not the reference state, for example, an excited state or a state with different spin multiplicity.

Size-consistency and size-extensivity problem

Davidson correction improves both size consistency and size extensivity of CISD energies.[2][4] Therefore, Davidson correction is frequently referred to in literature as size-consistency correction or size-extensivity correction.

However, neither Davidson correction itself nor the corrected energies are size-consistent or size-extensive. This is especially the case in larger molecules, where contribution from higher than quadruple excitations becomes more significant.

Corrections for Multi-reference CISD

Similar corrections exist for MR-CISD energies, including multi-reference Davidson correction, Pople correction, and others. These methods can be used to correct excited state energies.

See also

References

  1. Langhoff, Stephen R.; Davidson, Ernest R. (1 January 1974). "Configuration interaction calculations on the nitrogen molecule". International Journal of Quantum Chemistry 8 (1): 61–72. doi:10.1002/qua.560080106. 
  2. 2.0 2.1 Meissner, L. (1 May 1988). "Size-consistency corrections for configuration interaction calculations". Chemical Physics Letters 146 (3–4): 204–210. doi:10.1016/0009-2614(88)87431-1. Bibcode1988CPL...146..204M. 
  3. Sherrill, C. David. "Some Comments on the Davidson Correction". Georgia Institute of Technology. http://vergil.chemistry.gatech.edu/notes/davidcorr/davidcorr.html. Retrieved 11 November 2012. 
  4. Duch, Wl̸odzisl̸aw; Diercksen, Geerd H. F. (1 January 1994). "Size-extensivity corrections in configuration interaction methods". The Journal of Chemical Physics 101 (4): 3018. doi:10.1063/1.467615. Bibcode1994JChPh.101.3018D.