If the time-dependent HF orbitals, {ϕ i }, are expanded as linear combinations of a set of basis functions, ϕ k (r, t) = µ c µk (t)χ µ (r) ({c νi } being known as MO expansion coefficients), then we can define the corresponding density-matrix and Fock matrix elements to be: