
Electronic structure calculations
233
where H
e
is the electronic Hamiltonian operator given in full in equation (2.4) and E
e
is
the corresponding energy (it includes the nuclear–nuclear repulsion). Equation (B.1)isthe
Schr¨odinger equation for a fixed set of nuclear positions. Equation (B.2) describes the effect
of nuclear motion, with T
n
being the nuclear kinetic energy operator and E being the total
energy of the molecule. Notice that the potential energy ‘operator’ in (B.2)isthe energy
from solution of (B.1), so equation (B.1) must be solved before tackling (B.2).
In the case of a diatomic molecule, solution of (B.1)atvarious internuclear separations
gives the potential energy curve for that molecule. In polyatomic molecules consisting of
N atoms the energy E
e
is a function of 3N − 6or3N − 5 internal nuclear coordinates,
depending on whether the molecule is non-linear or linear, and it constitutes the potential
energy surface. The potential energy curve or surface defines the vibrational motion of a
molecule and therefore in order to predict vibrational frequencies equation (B.1) can be
solved at a variety of nuclear configurations to generate the potential energy surface, and
then (B.2)issubsequently solved. In fact in the majority of calculations equation (B.2)
is rarely solved explicitly to extract vibrational frequencies: a quicker route, based on the
evaluation of first and second derivatives of the total electronic energy with respect to the
internal nuclear coordinates, is usually employed [2].
An important point is that the wavefunction must satisfy the Pauliprinciple. In its simplest
form, this says that each electron in an atom or molecule has a unique set of quantum
numbers. In formal quantum mechanics, this corresponds to the insistence that the total
electronic wavefunction,
e
, must be antisymmetric with respect to the exchange of any
two electrons. A simple product wavefunction, one for each electron, of the type shown in
equation (2.5), will not satisfy the Pauli principle.
Take,asanexample, the case of H
2
in its ground electronic state, where the two electrons
are paired up in the 1σ
+
g
orbital. The wavefunctions for each electron are different, the
difference being not the spatial distributions of the two electrons, which are the same, but
the spins, which are opposite. We could therefore factor the wavefunction for each electron
into a common spatial part, which will be written as σ
+
g
, and a spin part, which is designated
as either α or β depending on whether the spin is ‘up’ or ‘down’. Notice that the spatial
wavefunction represents what is commonly referred to as an orbital,inthis case a molecular
orbital. The total electronic wavefunction can therefore be written as
e
= σ
+
g
(1)α(1)σ
+
g
(2)β(2) (B.3)
Unfortunately, this doesn’t satisfy the Pauli principle since an exchange of electrons 1 and 2
(equivalent to just switching the ‘1’ and ‘2’ labels in (B.3)) does not change the sign of the
wavefunction. However, the following function is antisymmetric with respect to electron
exchange:
e
= σ
+
g
(1)σ
+
g
(2)[α(1)β(2) − α(2)β(1)] (B.4)
This is an acceptable form of the wavefunction for a spin singlet since it satisfies the
Pauli principle and it retains, albeit in a slightly more complicated manner, the concept of
molecular orbitals.
Can similar antisymmetrized electronic wavefunctions be constructed for more compli-
cated molecules? The answer is yes, but written out in full algebraic form the expressions are