I want to make sure whether I do understand the transverse Ising model correctly or not. The classical Ising model describes the interaction between spins in a grid and the state of spins can be either +1 or -1. The transverse-field Ising model which is the quantum version of the classical Ising model. The possible outcomes of a single spin $\sigma_i$ can either be +1 and -1 therefore it can be represented as a Hamiltonian of Z-axis. And now if we want to present the coupling between two spins (based on the definition of the Ising model, see below)
$$H=-\sum_{ij}J_{ij}\sigma_{i}^{z}\sigma_{j}^{z}-\sum_{i}h_{i}\sigma_{i}^{z},$$
then we can present it, namely $\sigma_{i}^{z} \otimes \sigma_{j}^{z}$ in a matrix form as the following:
$$ \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 &0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{bmatrix} \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & -1 & 0 &0 \\ 0 & 0 & -1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{bmatrix} \begin{bmatrix} 1 & 0 & 0 & 0 \\ 0 & 1 & 0 &0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ \end{bmatrix}$$
And this is a product of two Pauli Z-matrix. Do I understand the transverse Ising model correct?
And why do we need unit matrixes here (based my university lecture, $\sigma_{i}^{z}:=I\otimes\ldots\otimes I \otimes\sigma_{z} \otimes I\otimes\ldots\otimes I$, where $\sigma_{z}$ is in the i-th position of the tensor product. I think the $\sigma_{z}$ is a single spin, isn't it?)? Are unit matrixes necessary? Thank you for your answer!