“Know how to solve every problem that has been solved.” “What I cannot create, I do not understand.” — Richard Feynman

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.

  1. base
  2. L1
  3. L2
  4. L3
  5. L4
  6. L5
  7. 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.

read it as

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.

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.

predict
Close the calculator. Water again, two other bases: STO-3G (a minimal basis) and cc-pVDZ (24 basis functions). How many particle-hole pairs in each — and what does the STO-3G number tell you about doing spectroscopy in a minimal basis?
why does this work
The output shows the lowest excitation dropping from 19.18 eV (bare) to 9.42 eV (TDA) — a 10 eV correction — while remaining 98.4% pure . If the interaction is strong enough to halve the energy, why is it too weak to mix the pairs?
invariant
The forty singlet TDA eigenvalues of water/6-31G sum to some number . The forty bare gaps sum to another number . What is the exact relation between and , and what everyday linear-algebra fact is it?

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.