Advanced Quantum Chemistry: Problem Workshop

Introduction

Advanced quantum-chemistry problems often fail because one convention changes silently midway. A point group is assigned for one geometry, orbital labels are copied from another axis system, a Hückel β sign is left unstated, or a Hartree–Fock eigenvalue is treated as a total energy. This workshop practices a disciplined ledger: specify geometry, symmetry, basis, electron count and approximation before calculating. Several short problems then connect symmetry-adapted orbitals, cyclic levels and SCF interpretation.

Core explanation

Problem A begins with ideal bent water in the yz plane, z along the angle bisector. Its point group is C₂v, with E, C₂(z), σ(xz) and σ(yz). The two H 1s orbitals yield characters (2,0,0,2) and reduce to A₁ + B₂. Oxygen 2s and 2p z are A₁; 2p y is B₂; 2p x is B₁. The six-function valence matrix therefore separates into 3A₁ + 2B₂ + B₁. The identity dimension check gives 3+2+1=6. A proposed 2p x–H sigma interaction is zero by symmetry in this basis, even though the atoms are spatially near.

Problem B uses ideal planar BF₃. Its three equivalent in-plane F sigma functions form A₁′ + E′ under D₃h. B 2s is A₁′ and B 2p x/2p y is E′, so matching sigma interactions are possible. B 2p z is A₂″ and has no partner among those in-plane sigma SALCs. This conclusion is basis-qualified: adding out-of-plane F p functions creates other possible pi-type combinations, and actual interaction strengths still depend on energy and overlap.

Problem C concerns a uniform four-site open Hückel chain with β<0 and four pi electrons. Energies are α+1.618β, α+0.618β, α−0.618β and α−1.618β, from lowest to highest. The bottom two spatial orbitals are doubly occupied, so the HOMO–LUMO gap is 1.236 β . If a student reports 2 β , they used the ethene two-site result. If they report a negative gap, they likely sorted levels without respecting β's sign. The gap remains a one-electron model quantity, not an exact absorption energy.

Problem D concerns benzene. Six p sites generate six MOs, though only four different energies occur because two pairs are degenerate. Six pi electrons fill one lowest orbital and the next degenerate pair, a closed shell. Comparing its occupied Hückel sum 6α+8β with three isolated ethene units at 6α+6β gives a model energy difference 2β. Since β<0, benzene's model value is lower. Calling that number a measured thermochemical resonance energy would exceed the model's scope.

Problem E concerns a Hartree–Fock calculation that reports a converged density and ε HOMO=−0.35 hartree. Koopmans' frozen-orbital estimate is a vertical ionisation of about +0.35 hartree, or 9.52 eV using 27.2114 eV per hartree. The total Hartree–Fock energy is not the sum of occupied eigenvalues. Before using the result, check the target spin, ⟨S²⟩ for an unrestricted state, numerical residuals and whether a different geometry or correlated calculation is required.

Across these examples, three checks recur. First, representation dimensions must equal basis size. Second, total electron populations must equal electron count. Third, an energy comparison must use one parameter and geometry convention on both sides. These checks do not replace physical reasoning but expose many mistakes before a polished answer is written.

Step-by-step reasoning

Read the question twice and write a one-line ledger of nuclear geometry, point group, axes, basis functions, electron count and β convention. Solve symmetry before energetic mixing, and solve eigenvalues before electron filling. For a numerical calculation, name the output quantity—orbital energy, total energy, vertical IP or excitation—then check units and model limitations.

Visual explanation

Make a four-panel worksheet. Panel one shows a point-group decision path, panel two a ligand SALC phase drawing, panel three a Hückel energy ladder with electron arrows, and panel four an SCF loop with a separate total-energy box. Put a red check mark next to dimension, electron and unit checks. The panels keep different types of answer from being conflated.

Real-world analogy

A laboratory notebook records sample identity, units and instrument settings before reporting a number. Quantum-chemistry work benefits from the same discipline: an orbital gap without β sign or an irrep without axis convention is like a concentration without units. The analogy highlights reproducibility, while the calculations still require chemistry-specific rules.

Real-world example

A chemist comparing substituted conjugated dyes might use a Hückel-like model for a first gap trend and a more complete calculation for quantitative wavelengths. If a dye's steric substituent twists its pi network, the simple equal-β prediction can fail. A transparent ledger makes it clear whether disagreement arises from geometry, parameterisation, selection-rule intensity or the many-electron excitation method.

Why?

Why combine symmetry and numerical methods in one workshop? Symmetry can remove matrix elements before diagonalisation, while eigenvalues and occupancies provide energy and electron-density information the point group cannot. A chemically meaningful answer needs both constraint and scale, plus a clear statement of which effects were approximated away.

Common misconception

An exact symmetry zero is different from a small energy or a weak allowed transition. Likewise, a Hückel orbital-gap estimate is different from a correlated excitation and a Fock orbital energy is different from total energy. Mixing these categories can produce numerically plausible but conceptually incorrect answers.

Worked example

Solve a two-site heteroatomic pi problem with H=[[−8,−2],[−2,−10]] eV and S=I. The average diagonal is −9 eV; half their difference has magnitude 1 eV. Eigenvalues are −9±√(1²+2²)=−9±√5, or about −11.236 and −6.764 eV. Two pi electrons occupy the lower level, giving a model gap about 4.472 eV. Unequal α values mean the lower MO coefficients are not equal-magnitude ethene combinations; one must solve the eigenvectors before assigning site populations.

Quick check

1. What does the identity character of a reducible orbital representation check? Answer: It equals the number of independent basis functions, so the dimensions of its irreducible components must sum to that count. 2. Is a converged Hartree–Fock HOMO value alone sufficient to report an adiabatic ionisation energy? Answer: No. Koopmans gives a frozen-orbital vertical estimate, while adiabatic ionisation includes cation geometry relaxation and other effects.

Exam focus

For every result, identify its quantity and scope in one sentence. Use explicit operation order, class sizes and beta sign. Show at least one independent consistency check and do not report extra precision beyond model assumptions.

Advanced insight

One can often predict qualitative block sizes, node counts or degeneracies before any numerical integral is evaluated. This precomputation becomes especially useful in larger calculations: if numerical output violates a symmetry-enforced degeneracy or basis-dimension count, the input geometry, program settings or state assignment deserves inspection.

Summary

The workshop's reliable pattern is to define geometry and basis, apply symmetry, solve the appropriate matrix, fill electrons and label the output quantity precisely. Water and BF₃ illustrate symmetry blocks; butadiene and benzene illustrate Hückel energies; Hartree–Fock illustrates why orbital and total energies differ. Dimension, electron and unit checks keep the answers coherent.

Practice questions

1. An ideal C₂v water sigma-basis calculation gives two B₁ ligand SALCs from two H 1s functions. What check identifies the error? Answer: The two-H basis reduces to A₁+B₂ in the stated axes, with dimension two. A B₁ H SALC is absent, so the phase or operation classification is wrong. 2. A four-site Hückel chain with β=−1.5 eV is assigned a HOMO–LUMO gap of 3.0 eV. What is the correct simple gap? Answer: Four-site uniform-chain gap is 1.236 β ≈1.854 eV. The 3.0 eV answer uses the two-site ethene formula 2 β . 3. A UHF calculation converges with Nα−Nβ=1 but ⟨S²⟩/ħ²=1.30. What should be reported? Answer: It targets S z=1/2, but the value exceeds the pure-doublet 0.75 benchmark and suggests spin contamination. Convergence and energy should be reported with that limitation.