Charge Density Matrix

jameson2

Firstly, I have been able to find almost nothing on this kind of question in text books or online anywhere. Most places (including my lecture notes) give at most the definition of the operator and that's all. One page if you're lucky out of a whole book. I'd kill for some examples, if you could direct me anywhere it'd be great. Here's my attempt anyway...

1. The problem statement, all variables and given/known data
(i)Consider a diatomic molecule, of atoms A and B. It's electronic structure can be given by the 1s orbitals of A and B, |A> and |B> (assumed orthonormal). Calculate the associated density matrix expanded over the basis set |A>, |B> as a function of number of electrons N=1,2,3,4. Use the nearest neighbour approximation with on site energies [tex] \epsilon_A=-4eV , \epsilon_B=4eV [/tex] and hopping parameter [tex] \gamma=-3eV [/tex]
(ii) By using ρ from part 1 evaluate the total energy for N=1,2,3,4.


2. Relevant equations
(i) [tex] \rho=\sum_{j}^{occupied} |\psi_j><\psi_j| [/tex]
(ii)[tex] \begin{bmatrix} \epsilon_A & \gamma \\ \gamma & \epsilon_B \end{bmatrix}\left[ \begin{array}{c} \psi_A \\ \psi_B \end{array} \right]= E\left[ \begin{array}{c} \psi_A \\ \psi_B \end{array} \right] [/tex]
(iii)[tex] L=Tr(L\rho) [/tex] for an operator L.

3. The attempt at a solution
(i) I don't know how to go about this, mainly the fact that the question brings in the number of electrons is confusing me. Here's what I've reasoned so far: Since it's a molecule, all the electrons are in the bonding state. Solving the matrix equation to find the bonding state, I get bonding energy of -5eV which leads to a bonding state of [tex] \frac{1}{\sqrt{10}} (3,1) [/tex]. Applying the formula for the density matrix for N electrons in the bonding state, I get [tex] \rho= N\frac{1}{10}\begin{bmatrix} 9 & 3 \\ 3 & 1 \end{bmatrix} [/tex]

(ii) Using the formula above, I guess the answer comes from [tex] E=Tr(H\rho) [/tex], but I need the density matrix first. Assuming I have the right density matrix, plugging in the Hamiltonian [tex] \begin{bmatrix} -4 & -3 \\ -3 & 4 \end{bmatrix} [/tex] I get the energy to be -5N for N=1,2,3,4.

Thanks alot for any help, or pointers towards somewhere that can help me. As I said, I'm having trouble finding info on this topic, you're my last chance!