Holstein–Herring method

Although the exchange energy becomes elusive at large internuclear systems, it is of prominent importance in theories of molecular binding and magnetism.

This splitting results from the symmetry under exchange of identical nuclei (Pauli exclusion principle).

The basic idea pioneered by Theodore Holstein, Conyers Herring and Boris M. Smirnov in the 1950-1960.

We consider states that are represented by even or odd functions with respect to behavior under space inversion.

This is denoted with the suffixes g and u from the German gerade and ungerade and are standard practice for the designation of electronic states of diatomic molecules, whereas for atomic states the terms even and odd are used.

For the hydrogen molecular ion, this is: For any gerade (or even) state, the electronic Schrödinger wave equation can be written in atomic units (

) as: For any ungerade (or odd) state, the corresponding wave equation can be written as: For simplicity, we assume real functions (although the result can be generalized to the complex case).

, located at the nuclei and write: This is similar to the LCAO (linear combination of atomic orbitals) method used in quantum chemistry, but we emphasize that the functions

are in general polarized i.e. they are not pure eigenfunctions of angular momentum with respect to their nuclear center, see also below).

collapse into the well-known atomic (hydrogenic) psi functions

as the mid-plane located exactly between the two nuclei (see diagram for hydrogen molecular ion for more details), with

representing the unit normal vector of this plane (which is parallel to the Cartesian

By considerations of symmetry: This implies that: Also, these localized functions are normalized, which leads to: and conversely.

Integration of the above in the whole space left to the mid-plane yields: and From a variation of the divergence theorem on the above, we finally obtain: where

From the latter, Conyers Herring was the first to show[3] that the lead term for the asymptotic expansion of the energy difference between the two lowest states of the hydrogen molecular ion, namely the first excited state

(as expressed in molecular notation—see graph for energy curves), was found to be: Previous calculations based on the LCAO of atomic orbitals had erroneously given a lead coefficient of

While it is true that for the Hydrogen molecular ion, the eigenenergies can be mathematically expressed in terms of a generalization of the Lambert W function, these asymptotic formulae are more useful in the long range and the Holstein–Herring method has a much wider range of applications than this particular molecule.

The Holstein–Herring formula had limited applications until around 1990 when Kwong-Tin Tang, Jan Peter Toennies, and C. L. Yiu[8] demonstrated that

can be a polarized wave function, i.e. an atomic wave function localized at a particular nucleus but perturbed by the other nuclear center, and consequently without apparent gerade or ungerade symmetry, and nonetheless the Holstein–Herring formula above can be used to generate the correct asymptotic series expansions for the exchange energies.

Subsequently, it has been applied with success to one-active electron systems.

Later, Scott et al. explained and clarified their results while sorting out subtle but important issues concerning the true convergence of the polarized wave function.

[9][10][11] The outcome meant that it was possible to solve for the asymptotic exchange energy splittings to any order.

The Holstein–Herring method has been extended to the two-active electron case i.e. the hydrogen molecule for the two lowest discrete states of

[13] The Holstein–Herring formula can be physically interpreted as the electron undergoing "quantum tunnelling" between both nuclei, thus creating a current whose flux through the mid-plane allows us to isolate the exchange energy.

Related to the tunnelling effect, a complementary interpretation from Sidney Coleman's Aspects of Symmetry (1985) has an "instanton" travelling near and about the classical paths within path integral formulation.

Energies (E) of the two lowest discrete states of the hydrogen molecular ion , as a function of inter-nuclear distance (R) in atomic units.