# Derivation of Time Dependent Hartree-Fock (TDHF) Equations

This derivation is largely based off the TD-DFT derivation (TDDFT and TDHF are very similar), found in the wonderful review, “Single-Reference ab Initio Methods for the Calculation of Excited States of Large Molecules” by Dreuw and Head-Gordon (Chem Reviews, 2005).

Now,

Substituting the first two expressions into the right hand side, these two expressions yields the Dirac form of the time-dependent Hartree-Fock equations,

where ${\mathbf{F}}$ is the Fock matrix and ${\mathbf{P} = \mathbf{CC}^{\dagger}}$ is the density matrix, same as in the time independent Hartree-Fock equations. Before a time dependent perturbation is applied, we assume that the system is in its electronic ground state, such that

and

with ${\mathbf{P}^{(0)}}$ and ${\mathbf{F}^{(0)}}$ as the unperturbed density and Fock matrix, respectively. The first condition comes from the fact that commuting matrices share a common set of eigenvectors, and the second comes from the fact that at convergence the eigenvectors of the Fock matrix are orthonormal. (Recall that the density matrix may be written as ${\mathbf{P} = \mathbf{CC}^{\dagger}}$, where the columns of ${\mathbf{C}}$ form an orthonormal set, e.g. ${\mathbf{C}^{\dagger}\mathbf{C} = \mathbf{1}}$). The Fock matrix elements are given by

Note the Fock matrix dependence on the density matrix. At convergence, the Fock matrix is diagonal, with the elements corresponding to the orbital energies ${\epsilon}$.

Furthermore, the density matrix enjoys the following relations:

Following convention, indices ${\{i,j,...\}}$ correspond to occupied orbitals, ${\{a,b,...\}}$ correspond to unoccupied orbitals, and ${\{p,q,...\}}$ correspond to general orbitals. In general perturbation theory, a perturbed wavefunction (or density matrix, for our purposes) can be decomposed, to first order, as

and likewise for the density-dependent Fock matrix

where the superscript (1) indicates the first-order time dependent change. Inserting the expressions 13 and 14 into the time-dependent equation 6 and collecting only the first order terms yields

The time-dependent perturbation can be described by a single Fourier component,

The matrix ${f}$ is a one-electron operator, which describes the applied perturbation. This perturbation acts on the electron density, resulting in a first order Fock matrix response

so that the overall first order change in the Fock matrix is

Similarly, the first order density response is given as

with ${d_{pq}}$ as the perturbation densities. Plugging the above expressions for ${\mathbf{F}^{(1)}}$ and ${\mathbf{P}^{(1)}}$ into the first order TDHF equation 15, and collecting the terms containing ${e^{-i\omega t}}$ gives (after some algebra):

The terms multiplied by ${e^{i \omega t}}$ give the complex conjugate of the above expression. Because the density matrix must be idempotent, we can show that

If ${\mathbf{PP} = \mathbf{P}}$, and we expand ${\mathbf{P}}$ as a perturbation with arbitrary scalar ${\lambda}$, we have:

Equating like powers of ${\lambda}$ yields the expression 21, as well as our original expression 8. Using 21 allows us to restrict the elements of ${d_{pq}}$, utilizing the nature of the zeroth order density matrix given in 11. As an example, consider the element ${d_{ii}}$:

which is only true if ${P_{ii}^{(1)} = 0}$. Similar arguments show that the only contributing blocks of ${d_{pq}}$ are ${d_{ia}}$ and ${d_{ai}}$. In other words, the only contributions to the TDHF equations are the occupied-virtual and virtual-occupied blocks. All virtual-virtual and occupied-occupied blocks are necessarily zero. With this in mind, along with the diagonal nature of the unperturbed Fock matrix, yields

and

with ${x_{ai} = d_{ai}}$ and ${y_{ai} = d_{ia}}$, as is convention. The derivatives inside the expression evaluate as follows:

Finally, we make the assumption that the perturbation is infinitely small, which sends ${f_{ai} = f_{ia} \rightarrow 0}$, and we recover (recognizing that ${F_{pp}^{(0)} = \epsilon_p}$ and ${P_{ii}^{(0)} = 1}$):

Which is the non-Hermitian eigenvalue TDHF equation. The elements of ${\mathbf{A}}$ and ${\mathbf{B}}$ are

and