In the last article in the series we focused on Quantum Many Body Systems and how we try to work on the low energy physics using a few powerful intuitive algorithms like DMRG (Density Matrix Renormalization Group). We also further looked into the dessimation procedure of DMRG which helps in approximation a Dd basis block to a D basis block thus reducing the exponential hurdle. Before going into the link between DMRG and Tensor Network states we will be looking at one of the most simple Tensor Network states which is the Matrix Product State (MPS) and how we decompose an arbitrary quantum state into one.
Let's say we have a system with n particles each with d degrees of freedom, the Hilbert space seems to be exponentially large where the number of states ranges to dn. But when we consider a local Hamiltonian where there is a gap between the ground state and the excited state, the pool of states is not very big compared to the full Hilbert Space. Here comes the Matrix product States which help in parametrization and DMRG earlier being difficult to work on the system now accepts it beautifully.
You will be seeing a lot of SVD (Singular Value Decomposition) and Schmidt Decomposition in the articles so you can refer to the previous articles for it or can google them for better understanding.
We will be further looking into the details of decomposing an arbitrary quantum state into MPS. Let's examine a lattice of N sites with every site having d dimensionality. The sites being represented by α In general we will be looking at a one dimensional lattice though it can be extended to higher dimensions too but they are not easily manageable. The state being :
∣ψ⟩=∑α1,...αNc∣α1,....αN⟩
The vector space being exponential here with the coefficients ranging to dN. We would like to bring this state to a more local self without affecting the non-locality of the whole state. We will be looking into the left canonical notation of the Matrix Product State. We will see a heavy usage of SVD for the given decomposition.
The state vector above has dN components. So what we are trying to do here is to represent the coefficient set into a form of matrix. The above state being normalized. We start by reshaping our state vector in a matrix. As the total number of components are dN, we can represent it with a matrix of dimensionality d×dN−1.
The matrix being ψα1,(α2...αN) and the coefficients cα1...αN being related to the matrix in the same way.
On applying SVD to the above matrix we get the given decomposition :
ψα1,(α2...αN)=∑a1r1Uσ1,a1Sa1,a1(V†)(a1,α2...αN)
Thus we decompose the amplitude or the coefficients using the above SVD criteria decomposing it into a product of matrices.
Next we combine the last two matrices (S and V†) and reshape it back into the vector form as can be seen in the diagram below.
Copy
<IPython.core.display.Image object>
Image Source : imsc.res.in
Copy
<IPython.core.display.Image object>
Copy
Image Source : imsc.res.in
The matrix U can be decomposed in a set of d row vectors and can be represented as Aaα where Uα1,a1 = Aa1α1 , α1 ranging from 1 to d. We reshape the above vector ψ2 into matrix form which has a dimensionality of r1d×d(N−2) where r1 is the rank of the singular matrix in the SVD equation and r1≤d. The equation thus comes out to be:
Again we apply SVD to the ψ matrix and get:
ψα1,(α2...αN)=∑a1r1∑a2r2AaαUa1σ2,a2Sa2,a2(V†)(a2,α3...αN)
Copy
<IPython.core.display.Image object>
In the first equation we replaced the unitary U by a set of row vectors of the form Aa1α1 while here we replace the matrix by a set of d matrices of dimensions r1×r2 and subsequently S and V are multiplied to forma matrix of the dimensionality r2d×dN−3
Thus the coefficient:
Copy
<IPython.core.display.Image object>
can be conveniently written as:
Copy
<IPython.core.display.Image object>
The quantum state can be written as:
cα1....αN=∑a1....aN−1Aα1Aα2....AαN−1AαN∣α1....αN⟩