Problemas para obtener la representación matricial de un hamiltoniano de 4 estados

$\newcommand{\bra}[1]{\left\langle #1 \right|} \newcommand{\ket}[1]{\left| #1 \right\rangle}$Tengo un hamiltoniano simple de 4 estados y estoy tratando de encontrar la representación de la matriz (para determinar los valores propios de energía para el sistema); sin embargo, sigo terminando con la matriz incorrecta, por lo que agradecería alguna ayuda para detectar qué concepto estoy malinterpretando.

Tenemos que el hamiltoniano viene dado por:

$$\begin{equation} \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| \} \end{equation}$$

Usando la distributividad de la suma sobre el espacio formado por el producto tensorial$\mathbb{C}^{4} \otimes \mathbb{C}^{4}$, tenemos que esto puede ser escrito por:

$$\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\}$$

Podemos reducir aún más esto observando que:$\sum_{n=1}^{4}k = 4k$,$\forall k \in \mathbb{C}^{4}\otimes \mathbb{C}^{4}$, y por lo tanto:

$$\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}$$ Si trabajamos en el espacio donde: $$\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}$$

Entonces nosotros tenemos:

$$\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}$$

Y:

$$\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}$$

Así tenemos:

$$\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}$$

Ampliando tenemos:

$$\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}$$

Lo cual es incorrecto, ya que sé que los valores propios son$\{E_{0}\pm W\}$, y el hamiltoniano se puede escribir en forma matricial como:

$$\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}$$

Sin embargo, no estoy seguro de dónde he cometido mi error, ¿cada paso me parece válido?