# Trouble getting the matrix representation of a 4-state Hamiltonian

$\newcommand{\bra}[1]{\left\langle #1 \right|} \newcommand{\ket}[1]{\left| #1 \right\rangle}$I have a simple 4-state Hamiltonian and am trying to find the matrix representation (in order to determine the eigenvalues of energy for the system); however, I keep ending up with the wrong matrix, so I'd appreciate some help with spotting what concept I'm misunderstanding.

We have that the Hamiltonian is given by:

$$\hat{\mathbf{H}}=\sum_{n=1}^{4}E_{0}\left|n\right\rangle\left\langle n \right|+\sum_{n=1}^{4}W\{\left| n \right\rangle\left\langle n+1 \right| + \left|n+1\right\rangle\left\langle n \right|-\left|2\right\rangle\left\langle 3\right|-\left|3\right\rangle\left\langle 2 \right| \}$$

Using the distributivity of the sum over the space formed by the tensor product $\mathbb{C}^{4} \otimes \mathbb{C}^{4}$, we have that this can be written by:

$$\hat{\mathbf{H}}=E_{0}\sum_{n=1}^{4}\ket{n}\bra{n}+W\left\{\sum_{n=1}^{4}\ket{n}\bra{n+1}+\sum_{n=1}^{4}\ket{n+1}\bra{n}-\sum_{n=1}^{4}\ket{2}\bra{3}-\sum_{n=1}^{4}\ket{3}\bra{2}\right\}$$

We further can reduce this by noting that: $\sum_{n=1}^{4}k = 4k$, $\forall k \in \mathbb{C}^{4}\otimes \mathbb{C}^{4}$, and thus:

$$\hat{\mathbf{H}}=E_{0}\sum_{n=1}^{4}\ket{n}\bra{n}+W\sum_{n=1}^{4}\ket{n}\bra{n+1}+W\sum_{n=1}^{4}\ket{n+1}\bra{n}-4W\ket{2}\bra{3}-4W\ket{3}\bra{2}$$ If we work in the space where: $$\ket{n}=\begin{cases}\hat{\mathbf{e}}_{n} & n \in \{1,2,3,4\} \\ \mathbf{0} & n \not\in \{1,2,3,4\}\end{cases}$$

Then we have:

$$\sum_{n=1}^{4}\ket{n}\bra{n}=\mathbf{I}_{4\times 4} \quad\land\quad \ket{2}\bra{3} = \begin{pmatrix}0 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0\end{pmatrix} \quad\land\quad \ket{3}\bra{2}=\begin{pmatrix}0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0\end{pmatrix}$$

And:

$$\sum_{n=1}^{4}\ket{n}\bra{n+1}=\begin{pmatrix}0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ 0 & 0 & 0 & 0\end{pmatrix} \quad\land\quad \sum_{n=1}^{4}\ket{n+1}\bra{n}=\begin{pmatrix}0 & 0 & 0 & 0 \\ 1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \end{pmatrix}$$

Thus we have:

$$\hat{\mathbf{H}}=E_{0}\mathbf{I}_{4\times 4}+\begin{pmatrix}0 & W & 0 & 0 \\ 0 & 0 & W & 0 \\ 0 & 0 & 0 & W \\ 0 & 0 & 0 & 0\end{pmatrix}+\begin{pmatrix}0 & 0 & 0 & 0 \\ W & 0 & 0 & 0 \\ 0 & W & 0 & 0 \\ 0 & 0 & W & 0\end{pmatrix}-\begin{pmatrix}0 & 0 & 0 & 0 \\ 0 & 0 & 4W & 0 \\ 0 & 4W & 0 & 0 \\ 0 & 0 & 0 & 0\end{pmatrix}$$

Expanding, we have:

$$\hat{\mathbf{H}} = \begin{pmatrix}E_{0} & W & 0 & 0 \\ W & E_{0} & -3W & 0 \\ 0 & -3W & E_{0} & W \\ 0 & 0 & W & E_{0}\end{pmatrix}$$

Which is incorrect, as I know that the eigenvalues are $\{E_{0}\pm W\}$, and the Hamiltonian can be written in matrix form as:

$$\hat{\mathbf{H}} = \begin{pmatrix}E_{0} & W & 0 & 0 \\ W & E_{0} & 0 & 0 \\ 0 & 0 & E_{0} & W \\ 0 & 0 & W & E_{0}\end{pmatrix}$$

However, I'm not sure where I've made my error, each step seems valid to me?

If you want $$\hat{\mathbf{H}} = \begin{pmatrix}E_{0} & W & 0 & 0 \\ W & E_{0} & 0 & 0 \\ 0 & 0 & E_{0} & W \\ 0 & 0 & W & E_{0}\end{pmatrix},$$ the corresponding Hamiltonian should be $$\hat{\mathbf{H}}=\sum_{n=1}^{4}E_{0}\left|n\right\rangle\left\langle n \right|+\sum_{n=1}^{4}W\{\left| n \right\rangle\left\langle n+1 \right| + \left|n+1\right\rangle\left\langle n \right| \}-W\left|2\right\rangle\left\langle 3\right|-W\left|3\right\rangle\left\langle 2 \right|.$$ Here I suppose $\left| 5 \right\rangle$ does not exist.
If the Hamiltonian is as same as the one in your question, i.e., $$\hat{\mathbf{H}}=\sum_{n=1}^{4}E_{0}\left|n\right\rangle\left\langle n \right|+\sum_{n=1}^{4}W\{\left| n \right\rangle\left\langle n+1 \right| + \left|n+1\right\rangle\left\langle n \right|-\left|2\right\rangle\left\langle 3\right|-\left|3\right\rangle\left\langle 2 \right| \},$$ you will indeed get $$\hat{\mathbf{H}} = \begin{pmatrix}E_{0} & W & 0 & 0 \\ W & E_{0} & -3W & 0 \\ 0 & -3W & E_{0} & W \\ 0 & 0 & W & E_{0}\end{pmatrix}$$ Because you sum $[-\left|2\right\rangle\left\langle 3\right|-\left|3\right\rangle\left\langle 2 \right|]$ for 4 times.