Unit 6: Many-Electron Atoms & Quantum Chemical Bonding

In classical mechanics, identical particles can be distinguished by tracking their continuous trajectories through phase space. In quantum mechanics, the Heisenberg uncertainty principle prohibits continuous trajectory tracking. Identical particles are fundamentally and completely indistinguishable.

The Permutation Operator and Symmetrization Postulate

Let \(\hat{P}_{12}\) be the particle-exchange (permutation) operator that swaps all spatial and spin coordinates of particles 1 and 2:

Because the particles are identical, exchanging them cannot alter any physical observable. In particular, the probability density must remain unchanged:

Applying \(\hat{P}_{12}\) twice restores the original state:

This leads to the fundamental Symmetrization Postulate of quantum mechanics:

1. Bosons (integer spin, \(S = 0, 1, 2, \dots\)): Wavefunctions are strictly symmetric under particle exchange:

2. Fermions (half-integer spin, \(S = 1/2, 3/2, \dots\)): Wavefunctions are strictly antisymmetric under particle exchange:

Electrons are fermions (\(s = 1/2\)), so any valid electronic wavefunction must be antisymmetric with respect to the simultaneous exchange of spatial and spin coordinates of any two electrons.

The Pauli Exclusion Principle

A direct mathematical consequence of antisymmetry: if two electrons were to occupy the exact same spatial and spin quantum state \(\chi_a\), then exchanging them would yield:

The state cannot exist. Thus, no two electrons in an atom or molecule can possess the same set of four quantum numbers.

Slater Determinants

For an \(N\)-electron system with orthonormal spin-orbitals \(\{\chi_1, \chi_2, \dots, \chi_N\}\) (where \(\chi_i(\mathbf{x}) = \phi_i(\mathbf{r}) \sigma(s)\)):

Properties of the Slater determinant:

1. Automatic Antisymmetry: Swapping two electrons corresponds to interchanging two rows of the determinant, which naturally flips the sign: \(\det = -\det\).

2. Pauli Exclusion: If two spin-orbitals are identical (\(\chi_i = \chi_j\)), two columns are identical, causing the determinant to vanish identically (\(\det = 0\)).

3. Normalization: The factor \(\frac{1}{\sqrt{N!}}\) ensures \(\langle \Psi | \Psi \rangle = 1\) when the spin-orbitals are orthonormal.

The Hartree-Fock (HF) approximation seeks the single Slater determinant \(\Phi_0\) that minimizes the electronic energy expectation value via the variational principle.

The Electronic Hamiltonian

In atomic units:

where the one-electron core Hamiltonian is \(\hat{h}(i) = -\frac{1}{2}\nabla_i^2 - \sum_A \frac{Z_A}{r_{i A}}\).

Energy Expectation Value and Integrals

For a closed-shell system containing \(N\) electrons paired into \(N/2\) spatial orbitals \(\{\phi_1, \dots, \phi_{N/2}\}\), the Hartree-Fock energy is:

where the integrals are defined as:

1. One-Electron Core Integral:

2. Coulomb Integral \(J_{i j}\):

This represents the purely classical electrostatic repulsion between charge clouds \(|\phi_i|^2\) and \(|\phi_j|^2\). It is strictly positive: \(J_{i j} > 0\).

3. Exchange Integral \(K_{i j}\):

The exchange integral has no classical counterpart. It arises solely from the antisymmetry requirement for electrons of like spin. It is also strictly positive: \(K_{i j} > 0\), with \(K_{i i} = J_{i i}\).

The Hartree-Fock Equations

Applying the variational principle \(\delta E_{\text{HF}} = 0\) subject to orbital orthonormality constraints \(\langle \phi_i | \phi_j \rangle = \delta_{i j}\) yields the canonical Hartree-Fock eigenvalue equations:

where \(\hat{f}\) is the one-electron Fock operator:

Here the local Coulomb operator is \(\hat{J}_j \phi_i(\mathbf{r}_1) = \left[ \int \frac{|\phi_j(\mathbf{r}_2)|^2}{r_{12}} d\mathbf{r}_2 \right] \phi_i(\mathbf{r}_1)\), and the non-local exchange operator is \(\hat{K}_j \phi_i(\mathbf{r}_1) = \left[ \int \frac{\phi_j^*(\mathbf{r}_2)\phi_i(\mathbf{r}_2)}{r_{12}} d\mathbf{r}_2 \right] \phi_j(\mathbf{r}_1)\).

Because the Fock operator depends on its own eigenfunctions through \(\hat{J}_j\) and \(\hat{K}_j\), the equations must be solved iteratively until the orbital coefficients reach self-consistency: the Self-Consistent Field (SCF) procedure.

In multi-electron atoms, electrostatic electron-electron repulsions and spin-orbit couplings split electron configurations into discrete spectroscopic energy levels termed multiplets.

Russell-Saunders (\(L\)-\(S\)) Coupling Scheme

For light to medium atoms (\(Z \le 30\)), electrostatic repulsions dominate over spin-orbit coupling. Individual orbital angular momenta couple to form total orbital angular momentum \(\mathbf{L} = \sum_i \mathbf{l}_i\), and individual spins couple to form total spin \(\mathbf{S} = \sum_i \mathbf{s}_i\). The atomic state is designated by the Russell-Saunders term symbol:

where:

• \(2S + 1\) is the spin multiplicity (1: singlet, 2: doublet, 3: triplet, 4: quartet, etc.)
• \(L\) is the total orbital angular momentum designated by capital letters: \(S (L=0), P (L=1), D (L=2), F (L=3), G (L=4), \dots\)
• \(J\) is the total angular momentum, taking values \(J = |L - S|, |L - S| + 1, \dots, L + S\).

Hund's Rules for Ground-State Term Determination

For an equivalent electron subshell (e.g., \(p^2, p^3, d^4\)), Hund's empirical rules determine the lowest energy ground state:

1. Hund's First Rule (Maximum Multiplicity): The term with the maximum total spin \(S\) (maximum multiplicity \(2S+1\)) lies lowest in energy.

Physical Mechanism: Electrons with parallel spins must occupy different spatial orbitals by the Pauli exclusion principle, reducing Coulomb repulsion, while their exchange interaction lowers energy by \(-K_{i j}\).

2. Hund's Second Rule (Maximum \(L\)): For a given multiplicity, the term with the highest value of total orbital angular momentum \(L\) lies lowest in energy.

Physical Mechanism: Electrons revolving in the same orbital direction encounter each other less frequently than counter-rotating electrons.

3. Hund's Third Rule (Spin-Orbit Multiplet Ordering):

• For subshells that are less than half full (e.g., \(p^1, p^2, d^1 \dots d^4\)), the level with the lowest \(J\) lies lowest: \(J_{\text{ground}} = |L - S|\) (normal multiplet).
• For subshells that are more than half full (e.g., \(p^4, p^5, d^6 \dots d^9\)), the level with the highest \(J\) lies lowest: \(J_{\text{ground}} = L + S\) (inverted multiplet).
• For half-filled subshells (e.g., \(p^3, d^5\)), \(L = 0 \implies J = S\).

Molecules consist of multiple positively charged nuclei and multiple negatively charged electrons interacting through mutual Coulomb forces.

The Complete Molecular Hamiltonian

Physical Rationale of the Born-Oppenheimer Approximation

The mass of a proton or nucleus is at least 1,836 times greater than the electron mass (\(M_A / m_e \ge 1836\)). Consequently, electrons move roughly two orders of magnitude faster than nuclei (\(v_e \sim 100 v_N\)). On the timescale of electronic motion, the nuclei appear virtually stationary. Conversely, the heavy nuclei experience the time-averaged electronic probability distribution.

Mathematical Formulation

We approximate the total molecular wavefunction as a separable product:

1. Electronic Schrödinger Equation: For a fixed nuclear configuration \(\mathbf{R}\), solve for electronic eigenfunctions:

where the electronic energy \(U(\mathbf{R})\) includes nuclear-nuclear repulsion \(V_{N N}(\mathbf{R})\). As \(\mathbf{R}\) varies, \(U(\mathbf{R})\) traces out the Potential Energy Surface (PES).

2. Nuclear Schrödinger Equation: The nuclei move on the effective potential energy surface generated by the electrons:

The nuclear motion separates into center-of-mass translation, overall molecular rotation, and internal nuclear vibrations. The Born-Oppenheimer approximation is exceptionally accurate for ground electronic states, breaking down only near conical intersections or avoided crossings where non-adiabatic vibronic coupling becomes substantial.

The hydrogen molecule-ion \(\text{H}_2^+\) consists of two protons (A and B) separated by internuclear distance \(R\), and a single electron. It is the simplest chemical bond in nature.

Electronic Hamiltonian of \(\text{H}_2^+\)

In atomic units:

Linear Combination of Atomic Orbitals (LCAO)

We approximate the molecular orbital \(\psi_{\text{MO}}\) as a linear combination of hydrogen \(1s\) atomic orbitals centered on protons A and B:

Because the two protons are identical, the probability density must possess inversion symmetry through the midpoint (\(R/2\)): \(|c_A|^2 = |c_B|^2 \implies c_B = \pm c_A\).

1. Bonding Molecular Orbital (\(\sigma_g 1s\)):

2. Antibonding Molecular Orbital (\(\sigma_u^* 1s\)):

where \(S = \langle 1s_A | 1s_B \rangle = e^{-R} (1 + R + R^2/3)\) is the atomic overlap integral.

Electronic Energy Expectation Values

The resulting energies are:

where \(J' = \langle 1s_A | -1/r_B | 1s_A \rangle\) is the Coulomb attraction of electron cloud \(A\) to nucleus \(B\), and \(K' = \langle 1s_A | -1/r_B | 1s_B \rangle\) is the resonance (exchange) integral.

• For \(\psi_+\) (bonding): Constructive quantum interference increases electron probability density in the internuclear region between the two protons, screening their mutual repulsion and lowering total energy, forming a stable potential well with equilibrium bond length \(R_e \approx 1.32\text{ \AA}\) and dissociation energy \(D_e \approx 1.76\text{ eV}\).
• For \(\psi_-\) (antibonding): Destructive interference creates a nodal plane midway between the protons, depleting electron density between the nuclei and resulting in a purely repulsive potential curve for all \(R\).

The hydrogen molecule \(\text{H}_2\) represents the classic two-electron covalent chemical bond. Two competing theoretical frameworks describe its electronic structure.

Molecular Orbital (MO) Theory

In simple MO theory, both electrons are placed into the bonding spatial orbital \(\sigma_g 1s\):

Expanding the spatial product in terms of atomic orbitals:

Catastrophic MO Dissociation Failure: As \(R \rightarrow \infty\), the molecule should dissociate into two neutral hydrogen atoms (\(\text{H} + \text{H}\)). However, the simple MO wavefunction predicts a 50% probability of dissociating into ions (\(\text{H}^+ + \text{H}^-\)), severely overestimating electron correlation energy at large distances.

Heitler-London Valence Bond (VB) Theory

In 1927, Walter Heitler and Fritz London formulated the Valence Bond wavefunction by assigning one electron to each atomic orbital and symmetrizing the spatial part:

1. Singlet Ground State (\(^1\Sigma_g^+\), Covalent Bonding):

2. Triplet Excited State (\(^3\Sigma_u^+\), Purely Repulsive):

The VB wavefunction is 100% covalent and dissociates with exact correct physics to \(\text{H}(1s) + \text{H}(1s)\) as \(R \rightarrow \infty\).

Comparison of Energies

• Experimental \(\text{H}_2\): \(R_e = 0.741\text{ \AA}\), \(D_e = 4.75\text{ eV}\).
• Heitler-London VB: \(R_e = 0.869\text{ \AA}\), \(D_e = 3.14\text{ eV}\).
• Simple MO: \(R_e = 0.850\text{ \AA}\), \(D_e = 2.68\text{ eV}\).

Both theories capture the essence of covalent bonding; their equivalence is restored when full configuration interaction is included.

The Hartree-Fock method provides an exceptional single-determinant mean-field approximation, but by definition it neglects instantaneous electron-electron interactions.

The Correlation Energy

Per Löwdin's exact definition, the correlation energy \(E_{\text{corr}}\) is the difference between the exact non-relativistic energy \(E_{\text{exact}}\) and the Hartree-Fock limit energy \(E_{\text{HF}}\):

Because \(E_{\text{HF}}\) is variational, \(E_{\text{corr}}\) is always negative. Although \(E_{\text{corr}}\) typically constitutes only 1% to 2% of the total electronic energy, it is comparable in magnitude to chemical bond energies (typically \(1\text{ to }5\text{ eV}\) per bond) and is therefore crucial for thermochemical accuracy.

Configuration Interaction (CI)

To capture electron correlation, the exact electronic wavefunction is expanded as a linear combination of the Hartree-Fock reference determinant \(\Phi_0\) and excited determinants formed by promoting electrons from occupied orbitals \(i, j\) to virtual (unoccupied) orbitals \(a, b\):

• \(\Phi_i^a\): Singly excited determinants (Singles, S)
• \(\Phi_{i j}^{a b}\): Doubly excited determinants (Doubles, D)
• \(\Phi_{i j k}^{a b c}\): Triply excited determinants (Triples, T)

Resolution of the \(\text{H}_2\) Dissociation Problem via CI

In a minimal basis set, the two molecular orbitals for \(\text{H}_2\) are \(\sigma_g\) and \(\sigma_u^*\). There are two closed-shell singlet configurations:

1. Ground configuration: \(\Phi_0 = |\sigma_g \bar{\sigma}_g|\)
2. Doubly excited configuration: \(\Phi_1 = |\sigma_u^ \bar{\sigma}_u^|\)

By Brillouin's theorem, singly excited determinants do not interact directly with \(\Phi_0\) (\(\langle \Phi_0 | \hat{H} | \Phi_i^a \rangle = 0\)). Constructing the two-configuration CI wavefunction:

At equilibrium bond length (\(R = R_e\)), \(c_1 \approx 0.99\) and \(c_2 \approx -0.11\) (dominantly ground state). At infinite separation (\(R \rightarrow \infty\)), \(\sigma_g\) and \(\sigma_u^*\) become degenerate. The CI secular equation yields \(c_1 = 1/\sqrt{2}\) and \(c_2 = -1/\sqrt{2}\):

The ionic terms cancel identically! Two-configuration CI eliminates the unphysical ionic dissociation artifact and restores correct physical dissociation to two neutral hydrogen atoms.

Hückel Molecular Orbital (HMO) theory provides a simple yet predictive quantum framework for planar conjugated \(\pi\)-electron systems (such as ethylene, butadiene, benzene, and annulenes).

Fundamental Approximations of HMO Theory

1. \(\sigma\)-\(\pi\) Separability: The planar \(\sigma\)-bond framework is treated as a rigid, frozen electrostatic core. Only delocalized \(\pi\)-electrons formed from out-of-plane \(2p_z\) atomic orbitals are explicitly treated.

2. Matrix Element Simplifications:

In the linear combination \(\psi_\pi = \sum_i c_i 2p_{z, i}\):

• Coulomb Integral \(\alpha\): Energy of an electron in an isolated carbon \(2p_z\) orbital:

• Resonance Integral \(\beta\): Coupling between directly bonded adjacent carbon atoms:

• Zero Differential Overlap (ZDO): Overlap matrix is approximated as the identity:

HMO Solution for Benzene (\(\text{C}_6\text{H}_6\))

For cyclic benzene, the secular determinant is:

Letting \(x = \frac{\alpha - E}{\beta}\), the Frost-Musulin circle mnemonic yields the 6 orbital energy roots:

• \(k = 0\): \(E_1 = \alpha + 2\beta\) (non-degenerate bonding)
• \(k = 1, 5\): \(E_2 = E_3 = \alpha + \beta\) (doubly degenerate bonding)
• \(k = 2, 4\): \(E_4 = E_5 = \alpha - \beta\) (doubly degenerate antibonding)
• \(k = 3\): \(E_6 = \alpha - 2\beta\) (non-degenerate antibonding)

Total \(\pi\)-Electronic Energy and Delocalization Energy

Benzene has \(6\) \(\pi\)-electrons. Filling the lowest 3 bonding orbitals:

For three localized ethylene double bonds:

The Delocalization (Resonance) Energy is:

Hückel's \(4n + 2\) Rule

Planar, monocyclic, fully conjugated systems with \(4n + 2\) \(\pi\)-electrons possess closed-shell electronic configurations with special thermodynamic stability (Aromaticity). Systems with \(4n\) \(\pi\)-electrons possess partially filled non-bonding or antibonding degenerate orbitals, leading to destabilization (Antiaromaticity).

## Advanced Mathematical Supplement: Second Quantization Algebra

Second Quantization Algebra: Creation, Annihilation & Field Operators

In many-body quantum chemistry, tracking explicit coordinate labels \((\mathbf{r}_1, \dots, \mathbf{r}_N)\) and Slater determinants becomes unwieldy. The formalism of Second Quantization replaces coordinate wavefunctions with an algebraic representation on Fock Space:

Fermionic Creation and Annihilation Operators

Let \(\{ |\phi_p\rangle \}\) be a complete orthonormal basis of single-particle spin-orbitals. We define the creation operator \(a_p^\dagger\) and annihilation operator \(a_p\):

• \(a_p^\dagger |0\rangle = |\phi_p\rangle\) (creates an electron in spin-orbital \(p\))
• \(a_p |\phi_p\rangle = |0\rangle\) (annihilates the electron in spin-orbital \(p\))
• \(a_p |0\rangle = 0\) (annihilation on the vacuum state yields zero)

The Pauli exclusion principle and antisymmetry are enforced identically by the Canonical Anticommutation Relations (CAR):

Notice that \(\{ a_p^\dagger, a_p^\dagger \} = 2 (a_p^\dagger)^2 = 0 \implies (a_p^\dagger)^2 = 0\): attempting to create two electrons in the exact same spin-orbital gives zero, encoding the Pauli exclusion principle into the operator algebra!

Second Quantized Hamiltonian

The full non-relativistic molecular electronic Hamiltonian transforms into:

where the one-electron core integrals are:

and the two-electron repulsion integrals (in physicist's notation) are:

All many-electron Slater determinants are written as operator strings acting on vacuum:

Wick's theorem then evaluates all matrix elements purely through algebraic contractions, forming the mathematical backbone of modern coupled cluster and quantum computing algorithms.

## Research Monograph: Coupled Cluster Theory & The Gold Standard of Quantum Chemistry

The Exponential Ansatz

While Configuration Interaction (CI) truncated at doubles (CISD) suffers from size-extensivity failure (error scales linearly with system size \(N\), rendering it unsuited for thermochemistry), Coupled Cluster (CC) theory uses an exponential ansatz:

where \(|\Phi_0\rangle\) is the reference Hartree-Fock determinant, and the cluster operator \(\hat{T}\) is:

with excitation operators \(\hat{T}_1 = \sum_{i, a} t_i^a a_a^\dagger a_i\) (singles) and \(\hat{T}_2 = \frac{1}{4} \sum_{i, j, a, b} t_{i j}^{a b} a_a^\dagger a_b^\dagger a_j a_i\) (doubles).

Size Extensivity and Disconnected Clusters

Expanding the exponential:

Even when \(\hat{T}\) is truncated at doubles (\(\hat{T} \approx \hat{T}_1 + \hat{T}_2\), CCSD), the quadratic term \(\frac{1}{2} \hat{T}_2^2\) automatically generates quadruple excitations corresponding to two simultaneous, non-interacting pair excitations. This guarantees strict size extensivity and size consistency for arbitrary molecular sizes.

CCSD(T): The Gold Standard

In CCSD(T), single and double excitation cluster amplitudes are solved self-consistently to infinite order, while triple excitations (\(\hat{T}_3\)) are evaluated via a non-iterative fourth-order perturbation correction. CCSD(T) reliably achieves sub-chemical accuracy (\(< 1\text{ kcal/mol} \approx 0.043\text{ eV}\)) for equilibrium bond distances, vibrational frequencies, and reaction barrier heights, earning its designation as the universal benchmark standard in molecular quantum chemistry.