Functional variation in derivation of the Hartree-Fock equations
Is it not just because of the fact that $[aa|bb]$ is in chemists notation, meaning that $[aa|bb] = [a(1)a(1)|b(2)b(2)]$ and so you can just switch them around?