If I want to minimize the energy of a Slater determinant subject to the constraint that the spin orbitals are orthonormal (as in the Hartree-Fock approximation), I can use Lagrange's method of undetermined multiplier, i.e.
$$L[\{\chi_{a}\}] = E_{0}[\{\chi_{a}\}]-\sum_{a=1}^{N}\sum_{b=1}^{N}\varepsilon_{ba}([a|b]-\delta_{ab})$$
where $\{\chi_{a}\}$ are the spin orbitals, $E_{0}$ is the ground state energy, $[a|b]$ is the overlap integral between spin-orbitals $\chi_{a}$ and $\chi_{b}$ and $\varepsilon_{ba}$ is a Langrange multiplier. By setting the first variation of $L$ to zero I can proceed to minimize the energy of a single determinant $|\Psi_0\rangle=|\chi_{1}\chi_{2}\cdots\chi_{a}\chi_{b}\cdots\chi_{N}\rangle$.
A problem in Szabo and Ostlund's "Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory" asks me to prove that the Lagrange multipliers are elements of a Hermitian matrix, i.e. $\varepsilon_{ba} = \varepsilon_{ab}^{*}$, given that $L$ is real and $[a|b] = [b|a]^{*}$.
How would one go about proving this?
Where I am starting at the moment is that since $L$ is real, then $L = L^{*}$. If I find $L^{*}$ then making this equal to $L$ will show that it is necessary for the Lagrange multipliers to be elements of a Hermitian matrix. However, following through with this doesn't seem to get me to a point where I can say that $\varepsilon_{ba} = \varepsilon_{ab}^{*}$.