Electronic correlation

Last updated

Electronic correlation is the interaction between electrons in the electronic structure of a quantum system. The correlation energy is a measure of how much the movement of one electron is influenced by the presence of all other electrons.

Contents

Atomic and molecular systems

Electron correlation energy in terms of various levels of theory of solutions for the Schrodinger equation. Electron correlation.svg
Electron correlation energy in terms of various levels of theory of solutions for the Schrödinger equation.

Within the Hartree–Fock method of quantum chemistry, the antisymmetric wave function is approximated by a single Slater determinant. Exact wave functions, however, cannot generally be expressed as single determinants. The single-determinant approximation does not take into account Coulomb correlation, leading to a total electronic energy different from the exact solution of the non-relativistic Schrödinger equation within the Born–Oppenheimer approximation. Therefore, the Hartree–Fock limit is always above this exact energy. The difference is called the correlation energy, a term coined by Löwdin. [1] The concept of the correlation energy was studied earlier by Wigner. [2]

A certain amount of electron correlation is already considered within the HF approximation, found in the electron exchange term describing the correlation between electrons with parallel spin. This basic correlation prevents two parallel-spin electrons from being found at the same point in space and is often called Fermi correlation. Coulomb correlation, on the other hand, describes the correlation between the spatial position of electrons due to their Coulomb repulsion, and is responsible for chemically important effects such as London dispersion. There is also a correlation related to the overall symmetry or total spin of the considered system.

The word correlation energy has to be used with caution. First it is usually defined as the energy difference of a correlated method relative to the Hartree–Fock energy. But this is not the full correlation energy because some correlation is already included in HF. Secondly the correlation energy is highly dependent on the basis set used. The "exact" energy is the energy with full correlation and complete basis set.

Electron correlation is sometimes divided into dynamical and non-dynamical (static) correlation. Dynamical correlation is the correlation of the movement of electrons and is described under electron correlation dynamics [3] and also with the configuration interaction (CI) method. Static correlation is important for molecules where the ground state is well described only with more than one (nearly-)degenerate determinant. In this case the Hartree–Fock wavefunction (only one determinant) is qualitatively wrong. The multi-configurational self-consistent field (MCSCF) method takes account of this static correlation, but not dynamical correlation.

If one wants to calculate excitation energies (energy differences between the ground and excited states) one has to be careful that both states are equally balanced (e.g., Multireference configuration interaction).

Methods

In simple terms, the molecular orbitals of the Hartree–Fock method are optimized by evaluating the energy of an electron in each molecular orbital moving in the mean field of all other electrons, rather than including the instantaneous repulsion between electrons.

To account for electron correlation, there are many post-Hartree–Fock methods, including:

One of the most important methods for correcting for the missing correlation is the configuration interaction (CI) method. Starting with the Hartree–Fock wavefunction as the ground determinant, one takes a linear combination of the ground and excited determinants as the correlated wavefunction and optimizes the weighting factors according to the Variational Principle. When taking all possible excited determinants, one speaks of Full-CI. In a Full-CI wavefunction all electrons are fully correlated. For non-small molecules, Full-CI is much too computationally expensive. One truncates the CI expansion and gets well-correlated wavefunctions and well-correlated energies according to the level of truncation.

Perturbation theory gives correlated energies, but no new wavefunctions. PT is not variational. This means the calculated energy is not an upper bound for the exact energy. It is possible to partition Møller–Plesset perturbation theory energies via Interacting Quantum Atoms (IQA) energy partitioning (although most commonly the correlation energy is not partitioned). [4] This is an extension to the theory of Atoms in Molecules. IQA energy partitioning enables one to look in detail at the correlation energy contributions from individual atoms and atomic interactions. IQA correlation energy partitioning has also been shown to be possible with coupled cluster methods. [5] [6]

There are also combinations possible. E.g. one can have some nearly degenerate determinants for the multi-configurational self-consistent field method to account for static correlation and/or some truncated CI method for the biggest part of dynamical correlation and/or on top some perturbational ansatz for small perturbing (unimportant) determinants. Examples for those combinations are CASPT2 and SORCI.

This approach includes a term depending on interelectron distance into wavefunction. This leads to faster convergence in terms of basis set size than pure gaussian-type basis set, but requires calculation of more complex integrals. To simplify them, interelectron distances are expanded into a series making for simpler integrals. The idea of R12 methods is quite old, but practical implementations begun to appear only recently.

Crystalline systems

In condensed matter physics, electrons are typically described with reference to a periodic lattice of atomic nuclei. Non-interacting electrons are therefore typically described by Bloch waves, which correspond to the delocalized, symmetry adapted molecular orbitals used in molecules (while Wannier functions correspond to localized molecular orbitals). A number of important theoretical approximations have been proposed to explain electron correlations in these crystalline systems.

The Fermi liquid model of correlated electrons in metals is able to explain the temperature dependence of resistivity by electron-electron interactions. It also forms the basis for the BCS theory of superconductivity, which is the result of phonon-mediated electron-electron interactions.

Systems that escape a Fermi liquid description are said to be strongly-correlated. In them, interactions plays such an important role that qualitatively new phenomena emerge. [7] This is the case, for example, when the electrons are close to a metal-insulator transition. The Hubbard model is based on the tight-binding approximation, and can explain conductor-insulator transitions in Mott insulators such as transition metal oxides by the presence of repulsive Coulombic interactions between electrons. Its one-dimensional version is considered an archetype of the strong-correlations problem and displays many dramatic manifestations such as quasi-particle fractionalization. However, there is no exact solution of the Hubbard model in more than one dimension.

The RKKY interaction can explain electron spin correlations between unpaired inner shell electrons in different atoms in a conducting crystal by a second-order interaction that is mediated by conduction electrons.

The Tomonaga–Luttinger liquid model approximates second order electron-electron interactions as bosonic interactions.

Mathematical viewpoint

For two independent electrons a and b,

where ρ(ra,rb) represents the joint electronic density, or the probability density of finding electron a at ra and electron b at rb. Within this notation, ρ(ra,rb) dra drb represents the probability of finding the two electrons in their respective volume elements dra and drb.

If these two electrons are correlated, then the probability of finding electron a at a certain position in space depends on the position of electron b, and vice versa. In other words, the product of their independent density functions does not adequately describe the real situation. At small distances, the uncorrelated pair density is too large; at large distances, the uncorrelated pair density is too small (i.e. the electrons tend to "avoid each other").

See also

Related Research Articles

Density-functional theory (DFT) is a computational quantum mechanical modelling method used in physics, chemistry and materials science to investigate the electronic structure of many-body systems, in particular atoms, molecules, and the condensed phases. Using this theory, the properties of a many-electron system can be determined by using functionals, i.e. functions of another function. In the case of DFT, these are functionals of the spatially dependent electron density. DFT is among the most popular and versatile methods available in condensed-matter physics, computational physics, and computational chemistry.

Electron density or electronic density is the measure of the probability of an electron being present at an infinitesimal element of space surrounding any given point. It is a scalar quantity depending upon three spatial variables and is typically denoted as either or . The density is determined, through definition, by the normalised -electron wavefunction which itself depends upon variables. Conversely, the density determines the wave function modulo up to a phase factor, providing the formal foundation of density functional theory.

Coupled cluster (CC) is a numerical technique used for describing many-body systems. Its most common use is as one of several post-Hartree–Fock ab initio quantum chemistry methods in the field of computational chemistry, but it is also used in nuclear physics. Coupled cluster essentially takes the basic Hartree–Fock molecular orbital method and constructs multi-electron wavefunctions using the exponential cluster operator to account for electron correlation. Some of the most accurate calculations for small to medium-sized molecules use this method.

In computational physics and chemistry, the Hartree–Fock (HF) method is a method of approximation for the determination of the wave function and the energy of a quantum many-body system in a stationary state.

Psi is an ab initio computational chemistry package originally written by the research group of Henry F. Schaefer, III. Utilizing Psi, one can perform a calculation on a molecular system with various kinds of methods such as Hartree-Fock, Post-Hartree–Fock electron correlation methods, and density functional theory. The program can compute energies, optimize molecular geometries, and compute vibrational frequencies. The major part of the program is written in C++, while Python API is also available, which allows users to perform complex computations or automate tasks easily.

Møller–Plesset perturbation theory (MP) is one of several quantum chemistry post-Hartree–Fock ab initio methods in the field of computational chemistry. It improves on the Hartree–Fock method by adding electron correlation effects by means of Rayleigh–Schrödinger perturbation theory (RS-PT), usually to second (MP2), third (MP3) or fourth (MP4) order. Its main idea was published as early as 1934 by Christian Møller and Milton S. Plesset.

Configuration interaction (CI) is a post-Hartree–Fock linear variational method for solving the nonrelativistic Schrödinger equation within the Born–Oppenheimer approximation for a quantum chemical multi-electron system. Mathematically, configuration simply describes the linear combination of Slater determinants used for the wave function. In terms of a specification of orbital occupation (for instance, (1s)2(2s)2(2p)1...), interaction means the mixing (interaction) of different electronic configurations (states). Due to the long CPU time and large memory required for CI calculations, the method is limited to relatively small systems.

Multi-configurational self-consistent field (MCSCF) is a method in quantum chemistry used to generate qualitatively correct reference states of molecules in cases where Hartree–Fock and density functional theory are not adequate. It uses a linear combination of configuration state functions (CSF), or configuration determinants, to approximate the exact electronic wavefunction of an atom or molecule. In an MCSCF calculation, the set of coefficients of both the CSFs or determinants and the basis functions in the molecular orbitals are varied to obtain the total electronic wavefunction with the lowest possible energy. This method can be considered a combination between configuration interaction and Hartree–Fock.

In computational chemistry, post–Hartree–Fock (post-HF) methods are the set of methods developed to improve on the Hartree–Fock (HF), or self-consistent field (SCF) method. They add electron correlation which is a more accurate way of including the repulsions between electrons than in the Hartree–Fock method where repulsions are only averaged.

Koopmans' theorem states that in closed-shell Hartree–Fock theory (HF), the first ionization energy of a molecular system is equal to the negative of the orbital energy of the highest occupied molecular orbital (HOMO). This theorem is named after Tjalling Koopmans, who published this result in 1934.

Full configuration interaction is a linear variational approach which provides numerically exact solutions to the electronic time-independent, non-relativistic Schrödinger equation.

Local-density approximations (LDA) are a class of approximations to the exchange–correlation (XC) energy functional in density functional theory (DFT) that depend solely upon the value of the electronic density at each point in space. Many approaches can yield local approximations to the XC energy. However, overwhelmingly successful local approximations are those that have been derived from the homogeneous electron gas (HEG) model. In this regard, LDA is generally synonymous with functionals based on the HEG approximation, which are then applied to realistic systems.

Restricted open-shell Hartree–Fock (ROHF) is a variant of Hartree–Fock method for open shell molecules. It uses doubly occupied molecular orbitals as far as possible and then singly occupied orbitals for the unpaired electrons. This is the simple picture for open shell molecules but it is difficult to implement. The foundations of the ROHF method were first formulated by Clemens C. J. Roothaan in a celebrated paper and then extended by various authors, see e.g. for in-depth discussions.

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.

<span class="mw-page-title-main">Spartan (chemistry software)</span>

Spartan is a molecular modelling and computational chemistry application from Wavefunction. It contains code for molecular mechanics, semi-empirical methods, ab initio models, density functional models, post-Hartree–Fock models, and thermochemical recipes including G3(MP2) and T1. Quantum chemistry calculations in Spartan are powered by Q-Chem.

Ab initio quantum chemistry methods are computational chemistry methods based on quantum chemistry. The term ab initio was first used in quantum chemistry by Robert Parr and coworkers, including David Craig in a semiempirical study on the excited states of benzene. The background is described by Parr. Ab initio means "from first principles" or "from the beginning", implying that the only inputs into an ab initio calculation are physical constants. Ab initio quantum chemistry methods attempt to solve the electronic Schrödinger equation given the positions of the nuclei and the number of electrons in order to yield useful information such as electron densities, energies and other properties of the system. The ability to run these calculations has enabled theoretical chemists to solve a range of problems and their importance is highlighted by the awarding of the Nobel prize to John Pople and Walter Kohn.

<span class="mw-page-title-main">Nuclear structure</span> Structure of the atomic nucleus

Understanding the structure of the atomic nucleus is one of the central challenges in nuclear physics.

In computational chemistry, spin contamination is the artificial mixing of different electronic spin-states. This can occur when an approximate orbital-based wave function is represented in an unrestricted form – that is, when the spatial parts of α and β spin-orbitals are permitted to differ. Approximate wave functions with a high degree of spin contamination are undesirable. In particular, they are not eigenfunctions of the total spin-squared operator, Ŝ2, but can formally be expanded in terms of pure spin states of higher multiplicities.

In quantum chemistry, Brillouin's theorem, proposed by the French physicist Léon Brillouin in 1934, relates to Hartree–Fock wavefunctions. Hartree–Fock, or the self-consistent field method, is a non-relativistic method of generating approximate wavefunctions for a many-bodied quantum system, based on the assumption that each electron is exposed to an average of the positions of all other electrons, and that the solution is a linear combination of pre-specified basis functions.

<span class="mw-page-title-main">Complete active space perturbation theory</span>

Complete active space perturbation theory (CASPTn) is a multireference electron correlation method for computational investigation of molecular systems, especially for those with heavy atoms such as transition metals, lanthanides, and actinides. It can be used, for instance, to describe electronic states of a system, when single reference methods and density functional theory cannot be used, and for heavy atom systems for which quasi-relativistic approaches are not appropriate.

References

  1. ↑ Löwdin, Per-Olov (March 1955). "Quantum Theory of Many-Particle Systems. III. Extension of the Hartree–Fock Scheme to Include Degenerate Systems and Correlation Effects". Physical Review . 97 (6): 1509–1520. Bibcode:1955PhRv...97.1509L. doi:10.1103/PhysRev.97.1509.
  2. ↑ Wigner, E. (December 1, 1934). "On the Interaction of Electrons in Metals". Physical Review. 46 (11): 1002–1011. Bibcode:1934PhRv...46.1002W. doi:10.1103/PhysRev.46.1002.
  3. ↑ J.H. McGuire, "Electron Correlation Dynamics in Atomic Collisions", Cambridge University Press, 1997
  4. ↑ McDonagh, James L.; Vincent, Mark A.; Popelier, Paul L.A. (October 2016). "Partitioning dynamic electron correlation energy: Viewing Møller-Plesset correlation energies through Interacting Quantum Atom (IQA) energy partitioning". Chemical Physics Letters. 662: 228–234. Bibcode:2016CPL...662..228M. doi: 10.1016/j.cplett.2016.09.019 .
  5. ↑ Holguín-Gallego, Fernando José; Chávez-Calvillo, Rodrigo; García-Revilla, Marco; Francisco, Evelio; Pendás, Ángel Martín; Rocha-Rinza, Tomás (July 15, 2016). "Electron correlation in the interacting quantum atoms partition via coupled-cluster lagrangian densities". Journal of Computational Chemistry. 37 (19): 1753–1765. doi:10.1002/jcc.24372. ISSN   1096-987X. PMID   27237084. S2CID   21224355.
  6. ↑ McDonagh, James L.; Silva, Arnaldo F.; Vincent, Mark A.; Popelier, Paul L. A. (April 12, 2017). "Quantifying Electron Correlation of the Chemical Bond" (PDF). The Journal of Physical Chemistry Letters. 8 (9): 1937–1942. doi:10.1021/acs.jpclett.7b00535. ISSN   1948-7185. PMID   28402120. S2CID   24705205.
  7. ↑ Quintanilla, Jorge; Hooley, Chris (2009). "The strong-correlations puzzle" (PDF). Physics World . 22 (6): 32–37. Bibcode:2009PhyW...22f..32Q. doi:10.1088/2058-7058/22/06/38. ISSN   0953-8585.
This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.