I've been trying to better understand matrix product states (in order to implement them in code in the near future), so I'm considering small examples. I was wondering if I could get some clarification on how to express states as an MPS for this specific case.
In the 2 site Ising model, we have a Hamiltonian in the form:
$$ H = -(\sigma^z_1 \sigma^z_2 + \sigma^z_2 \sigma^z_1) - h (\sigma^x_1 + \sigma^x_2) $$
where I have set the value of the interaction coefficient to 1. When $h=0$, I would expect my ground state to be some superposition of $| \uparrow \uparrow \rangle$ and $| \downarrow \downarrow \rangle$.
If I'm understanding this reference properly (page 3), if I want to express $| \psi \rangle = | \uparrow \uparrow \rangle + | \downarrow \downarrow \rangle$ (currently not normalized) as a matrix, my goal is to obtain two matrices that multiply to make:
$$ A = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} $$
The solution is to pick $\begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} $
So my questions:
Are matrix product states unique? It seems that $A$ is effectively $\begin{bmatrix} | \uparrow \uparrow \rangle & 0 \\ 0 & | \downarrow \downarrow \rangle \end{bmatrix} $, but I could easily change this to be $\begin{bmatrix} | \downarrow \downarrow \rangle & 0 \\ 0 & |\uparrow \uparrow \rangle \end{bmatrix} $
Physically, how do I interpret $\begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} $? What do the matrix elements mean?
In terms of actual implementation in a tensor network, can I express this MPS as a (2, 2, 2) tensor? Where the first two indices indicate that we have $2\times 2$ matrices, and the 3rd index shows that we have two of these $2 \times 2$ matrices?
I'm guessing I have a lot of misconceptions here.
