Let us consider tetrahedron with vertices, labeled as 1, 2, 3, 4. Let us have $N$ electrons which can be placed in this vertices. Let the hamiltonian of this system satisfy such condition:
$$ \langle i|\hat{H}|j\rangle = \Delta, \: i \neq j, \: 1 \leq i,j \leq 4 .$$ Our task is to find the energy spectrum of this system in cases when 1) $N = 1$; 2) $N = 2$; 3) $N = 3$; 4) $N = 4$.
I have proceeded in case of $N=1$. My approach was to obtain eigenvalues of single-particle hamiltonian directly. My result is that we have two energy levels: ground state with energy $-\Delta$ which is triply degenerated and excited state with energy $3\Delta$.
Now I am stuck with the case of two electrons. The problem is that we should consider electrons as fermions with spin. I have two ideas how to proceed further, but misunderstand how to make it real:
- To construct a good basis in the space of states. It should contain antisymmetrized mixtures of single-particle states. It is something like: $$ |\text{state}_1\!\rangle = \frac{1}{\sqrt{2}} (|1 ;2\!\rangle- \:|2;1\!\rangle)$$ And so on. The first number in the brackets corresponds to the position of the first electron and the second number corresponds to the second electron. On this step I am not sure how to account for spin variables. If we add spin variable, then in formula above, it will look like: $$ |\text{state}_1\!\rangle = \frac{1}{\sqrt{2}} (|1,\sigma_1 ;2, \sigma_2\!\rangle-|2,\sigma_2 ; 1,\sigma_1 \!\rangle)$$ The hamiltonian should be rewritten in terms of this basis states. I tried to "extend" initial hamiltonian doing this: $$ \hat{H} = \hat{H_1}+\hat{H_2}$$ $$ \langle i,\sigma_1;k,\sigma_2|\hat{H_1}|j,\sigma_x; k,\sigma_2\rangle \:= \Delta, \: i \neq j, \: 1 \leq i,j \leq 4, \; \sigma_x = \pm\frac{1}{2} . $$ $$ \langle i,\sigma_1;k,\sigma_2|\hat{H_2}|i,\sigma_1; m,\sigma_x\rangle \:= \Delta, \: k \neq m, \: 1 \leq k,m \leq 4, \; \sigma_x = \pm\frac{1}{2} . $$ But this does not include Pauli exclusion principle. It can be inserted in this formula by imposing restrictions on indices like: $$ \langle i,\sigma_1;k,\sigma_2|\hat{H_1}|j,\sigma_x; k,\sigma_2\rangle \:= \begin{cases} \Delta, \: i \neq j, \: 1 \leq i,j \leq 4, \; \sigma_x = \pm\frac{1}{2}, \; \{j,\sigma_x\} \neq \{k, \sigma_2\} \\ 0, \; \text{otherwise} . \end{cases}$$ But this expression now looks too complicated to consider its matrix elements in new basis. Moreover, the Pauli exclusion principle makes this matrix elements dependent from spin. This confuses me about the dimensionality of this matrix. If we forget about spin, then we have $4 \times 3 = 12$ pairs of different indices and therefore 12 basic states. If we then consider the spin indices, then, as I see, we will get $12 \times 4 = 48$ basic states, because each electron in pair has two options of spin. It seems like I am to recollect this $48 \times 48$ matrix, iterating over all the indices and tracking all the restrictions. It is hard to be done for me. If I had somehow the resultant matrix, then I would analyze its eigenvalues (this problem also looks scary in perspective). Is such procedure correct in principle?
- Probably the problem may be simplified in secondary quantization representation. The hamiltonian can be considered as combination of annihilation operator in one vertice and creation operator in another: $$ \hat{H} = \Delta \sum_{i \neq j} a^\dagger_i a_j + \Delta \sum_{i \neq j} b^\dagger_i b_j$$ Here $a$ corresponds to the first electron and $b$ corresponds to the second. But I am not sure about this expression because it doesn't include spin and it doesn't forbid electrons to violate Pauli exclusion principle (electrons don't see each other). Moreover, even if I had the expression like this, it would be hard for me to calculate spectrum of this combination. Maybe it can be worked out with commutative relations for fermion operators. Is it correct way to approach the problem?
- The last idea that I have is to include "exchange interaction" adding the term like: $$ \widehat{H_{\text{ex}}} = J \sigma_1 \sigma_2 $$ I tried to consider position and spin variables separately, found eigenvalues of $16\times 16$ "positional hamiltonian" ($-\Delta$ and $15\Delta$) and then added the split caused by exchange interaction: $-\frac{3J}{4}, \frac{J}{4}$. Therefore 64 "stationary states" are obtained. It looks not very good for me, because it again doesn't account directly for antisymmetry, I don't know where I can get $J$ from and 64 is also as I think incorrect number, because I expect 48 from what was written above (Pauli principle sorts out redundant states, as I think). But may be something can be done in this manner?
I would be really grateful for any comments and remarks.