Particle-Hole Excitations
Quantum Chemistry
What you need to know first 10 concepts, 6 layers
The requisite-knowledge inventory for this page, bottom-up: the primitives at the base, combined upward until you reach what this page assumes. Skim the layers you already own; start wherever the ground gets unfamiliar.
- base
- Linear algebraconcept
- Quantum mechanics (states & operators)
- L1
- L2
- L3
- L4
- L5
- ↳you are here
1 of these are concepts without a dedicated page yet — the grey chips. Following the linked ones first makes the rest land.
Take a converged mean field, lift one electron from an occupied orbital to a virtual orbital , and you have made a particle-hole excitation — the elementary move of molecular spectroscopy. The set of all such moves, one per occupied-virtual pair, is the coordinate system that Casida's equation is written in. This page is about that coordinate system: what a pair is, how many there are, what energy it costs at zeroth order, and why the hole deserves equal billing with the particle.
The ground state of a closed-shell molecule, at the mean-field level, is one Slater determinant : the lowest orbitals, each doubly filled. The simplest thing you can do to it is a single substitution — remove the electron in occupied orbital , place it in virtual orbital — producing the singly-excited determinant
Read the operators right to left: digs a hole in the occupied sea, drops the particle into an empty level. That ordered pair — hole index, particle index — is the whole identity of the excitation. Every absorption spectrum on this site starts life as a list of these substitutions.
How many single excitations does a molecule have? Exactly — a rectangular grid, one axis for the hole, one for the particle. Water in 6-31G: 13 basis functions, 5 occupied, 8 virtual, so the grid is pairs. Not "about forty" — exactly forty, fixed the moment you choose the basis, before any physics happens.
This grid is the vector space low-energy spectroscopy lives in. A general excited state, in linear response, is a vector over the grid — so much of , a little of , and so on — and the Casida matrix for water/6-31G is a matrix indexed by pairs on both sides. When quantum chemists say "the singles space," this grid is what they mean; grow the basis and the grid grows with the virtual axis.
Why insist on "particle-hole" instead of just "excited electron"? Because the vacancy left behind is not passive bookkeeping — it is a positively-charged quasiparticle with its own dynamics. The promoted electron and its hole attract each other Coulombically, and that attraction is enormous: in the run below, the bare cost of the promotion is 19.2 eV, but the actual lowest excitation lands at 9.4 eV. The electron-hole attraction claws back almost 10 eV — half the zeroth-order energy.
A bound electron-hole pair is an exciton, and the language is deliberately shared with solid-state physics: same object, wildly different binding scale (eV in a small molecule, meV in a semiconductor, because the pair delocalizes). The interaction terms that the Casida coupling matrix adds on top of the bare grid are precisely this electron-hole physics — restoring what the frozen mean field left out.
If the electrons truly didn't interact, the excitation spectrum would cost nothing to compute: promoting would cost exactly the orbital-energy difference,
and the full spectrum would be the sorted list of all forty differences. That list is the bare particle-hole ladder, and it is genuinely the zeroth order of the real theory: the Casida matrix has exactly these differences on its diagonal, with interaction supplying everything off it (and part of the diagonal too). For water's lowest pairs the bare ladder reads 19.18, 20.80, 21.80 eV — wrong by a factor of two, but wrong systematically, and in the right order for the first few states.
This is the same trust-the-differences logic as on the orbital ladder page: individual eigenvalues are unreliable, but differences within one consistent ladder are a well-conditioned starting point that the response machinery then corrects.
One spatial pair is secretly two excitations. The electron carries spin, and after the promotion the open-shell pair can couple to total spin 0 or 1:
The singlet is what light makes (photons don't flip spins); the triplet is dark but always lies lower, because parallel spins keep the electrons apart and reduce their repulsion — Hund's rule, relocated to excited states. At zeroth order the two are degenerate (both cost ); the interaction splits them. That is why the Casida problem for a closed-shell molecule falls into separate singlet and triplet blocks of size each — same grid, two spin channels.
Three ladders, one molecule
Water in 6-31G at the Hartree-Fock level. Ladder one is the bare grid — orbital-energy differences, no interaction. Ladder two diagonalizes the coupling matrix over the grid (the Tamm-Dancoff approximation, identical to CIS here). Ladder three is the full response problem with de-excitations included.
import numpy as np
from pyscf import gto, scf, tdscf
HA_EV = 27.211386
mol = gto.M(
atom="""O 0.0000 0.0000 0.1173
H 0.0000 0.7572 -0.4692
H 0.0000 -0.7572 -0.4692""",
basis="6-31g", verbose=0)
mf = scf.RHF(mol).run()
nocc = mol.nelectron // 2
nvir = mf.mo_coeff.shape[1] - nocc
eps = mf.mo_energy
# Ladder 1: bare particle-hole gaps, no interaction at all
pairs = [(i, a, eps[nocc+a] - eps[i])
for i in range(nocc) for a in range(nvir)]
pairs.sort(key=lambda t: t[2])
# Ladder 2: TDA/CIS (diagonalize A). Ladder 3: full TDHF (A and B).
tda = tdscf.TDA(mf); tda.nstates = 5; tda.kernel()
tdhf = tdscf.TDHF(mf); tdhf.nstates = 5; tdhf.kernel()
# Composition of the lowest excitation over the p-h grid
x = tda.xy[0][0] # X_ia, shape (nocc, nvir)
w = 2.0 * x**2 # weights, sum to 1 for a singlet H2O / 6-31G: nocc = 5, nvir = 8 -> 40 particle-hole pairs
(the Casida singlet block is 40x40)
# bare (i->a) Delta eps TDA/CIS TDHF (eV, singlets)
1 5->6 19.184 9.421 9.364
2 4->6 20.796 11.358 11.285
3 5->7 21.799 11.865 11.783
4 4->7 23.411 13.949 13.859
5 3->6 24.857 15.535 15.482
Lowest singlet (TDA, 9.421 eV) -- p-h pair weights |X_ia|^2:
5 -> 6 weight 0.984 (bare gap 19.184 eV)
5 -> 10 weight 0.013 (bare gap 45.939 eV)
5 -> 13 weight 0.003 (bare gap 59.798 eV) Two things are true at once, and they pull in opposite directions. The interaction is huge: every excitation drops by 9-10 eV from its bare gap, the electron-hole attraction from framing 3. And yet the interaction barely mixes: the lowest excitation is 98.4% the single pair , with the leftover 1.6% scattered over pairs that share the same hole. The big correction sits on the diagonal of the coupling matrix — each pair's energy shifts a lot, but the pairs keep their identities. That is what makes the particle-hole grid a good coordinate system rather than just a basis: the answer is not merely expressible in it, it is nearly diagonal in it. (The last step, from TDA to TDHF, is worth only ~0.06 eV here — the de-excitation block is a garnish on this particular molecule.)
Try First
Each prompt asks a checkable question about the working code or math above — predict an output, derive a sign, state an invariant, find a bug. Commit to an answer before clicking "reveal." That commitment is the whole point: if your answer matched, you understand the piece you were looking at; if it didn't, that's the part worth re-reading.
Where this goes
You now have the basis (the pair grid), the zeroth order (the bare ladder), and the physics the mean field forgot (the electron-hole interaction). Assembling those three into an eigenvalue problem — with the interaction as a kernel and the de-excitations riding along in the block — is exactly Casida's equation, and the time-dependent story of where that kernel comes from is time-dependent DFT (with the matrix-free large-scale version at Liouville-Lanczos). Upstream, the rungs themselves are explained on Kohn-Sham orbitals and built by Hartree-Fock.