$ to each of them. Assume that every pair of neighboring virtual spins which correspond to different sites are initially in an (unnormalized) maximally entangled state $|I\rangle=\sum_{\alpha=1}^{D} |\alpha,\alpha\rangle$ often referred to as entangled *bond*. Then apply a map $$\label{curlyA}{\cal
A}=\sum_{i=1}^d\sum_{\alpha,\beta=1}^{D} A_{i,\alpha,\beta}
|i\rangle\langle\alpha,\beta|$$ to each of the $N$ sites. Here and in the following Greek indices correspond to the virtual systems. By writing $A_{i}$ for the $D\times D$ matrix with elements $A_{i,\alpha,\beta}$ we get that the coefficients of the final state when expressed in terms of a product basis are given by a matrix product $\operatorname{tr}\left[A_{i_1}A_{i_2}\cdots
A_{i_N}\right]$. In general the dimension of the entangled state $|I\rangle$ and the map ${\cal A}$ can both be site-dependent and we write $A_i^{[k]}$ for the $D_k\times D_{k+1}$ matrix corresponding to site $k\in\{1,\ldots,N\}$. States obtained in this way have then the form $$\label{MPS0}
|\psi\rangle=\sum_{i_1,\ldots,i_N=1}^d
\operatorname{tr}\left[A_{i_1}^{[1]}A_{i_2}^{[2]}\cdots
A_{i_N}^{[N]}\right]|i_1,i_2,\ldots,i_N\rangle\;,$$ and are called *matrix product states* [@MPSorigin]. As shown in [@Vi03] *every* state can be represented in this way if only the bond dimensions $D_k$ are sufficiently large. Hence, Eq.(\[MPS0\]) is a representation of states rather than the characterization of a specific class. However, typically states are referred to as MPS if they have a MPS-representation with small $D=\max_k D_k$ which (in the case of a sequence of states) does in particular not grow with $N$. Note that $\psi$ in Eq.(\[MPS0\]) is in general not normalized and that its MPS
$ to each of them. Assume that every pair of neighboring virtual spins which represent to unlike sites are initially in an (unnormalized) maximally entangled department of state $ |I\rangle=\sum_{\alpha=1}^{D } |\alpha,\alpha\rangle$ much referred to as entangle * shackle *. Then apply a function $ $ \label{curlyA}{\cal 
 A}=\sum_{i=1}^d\sum_{\alpha,\beta=1}^{D } A_{i,\alpha,\beta } 
 |i\rangle\langle\alpha,\beta|$$ to each of the $ N$ web site. Here and in the following Greek indices represent to the virtual systems. By writing $ A_{i}$ for the $ D\times D$ matrix with element $ A_{i,\alpha,\beta}$ we get that the coefficients of the final state when expressed in terminus of a product basis are given by a matrix intersection $ \operatorname{tr}\left[A_{i_1}A_{i_2}\cdots 
 A_{i_N}\right]$. In general the dimension of the entangled state $ |I\rangle$ and the function $ { \cal A}$ can both be site - subject and we write $ A_i^{[k]}$ for the $ D_k\times D_{k+1}$ matrix corresponding to site $ k\in\{1,\ldots, N\}$. States prevail in this way have then the form $ $ \label{MPS0 } 
 |\psi\rangle=\sum_{i_1,\ldots, i_N=1}^d 
 \operatorname{tr}\left[A_{i_1}^{[1]}A_{i_2}^{[2]}\cdots 
 A_{i_N}^{[N]}\right]|i_1,i_2,\ldots, i_N\rangle\;,$$ and are called * matrix product states * [ @MPSorigin ]. As shown in [ @Vi03 ] * every * state can be represented in this way if only the bond dimensions $ D_k$ are sufficiently large. Hence, Eq.(\[MPS0\ ]) is a representation of states preferably than the characterization of a specific class. However, typically department of state are refer to as MPS if they have a MPS - representation with small $ D=\max_k D_k$ which (in the shell of a sequence of states) does in particular not grow with $ N$. Note that $ \psi$ in Eq.(\[MPS0\ ]) is in cosmopolitan not normalized and that its MPS
$ to each of them. Assume thau every pair of ngithboriig virthal spinr which correspond to differxnt wites are initially in an (uvnormalizvd) maximaoly tntangled state $|I\czngle=\sum_{\alpha=1}^{S} |\alpka,\elpha\rangle$ oftgn referred do as entangleg *coud*. Then apply a map $$\label{curlyA}{\cal
A}=\stm_{i=1}^d\sum_{\sloha,\beta=1}^{D} A_{i,\alpra,\beua}
|i\wangmv\lcngle\alpha,\beta|$$ to each of the $N$ sites. Iere and in the following Greek indices clrredpond to the virtuwl systems. Vy wwuting $A_{i}$ for the $D\times D$ matrix wjth elements $A_{i,\alpha,\beta}$ we get that the coeffucuenhv of the fiial stwte when expvvssed it terms of a product nasis arw given by a matrix pcoduct $\operatorname{tt}\left[A_{i_1}A_{i_2}\cgocs
A_{i_N}\right]$. In general tye dikenshon ud tfe tntengmed stwte $|I\rangle$ ahd the map ${\xal A}$ can both be sotq-eependent and we wrytq $A_i^{[k]}$ for the $D_k\times D_{k+1}$ matrix correspmndjng to site $k\in\{1,\ldots,N\}$. Srates obtained in thid way havq then the form $$\label{MPS0}
|\psi\rangle=\sum_{i_1,\ldots,i_N=1}^d
\operdtornemd{tr}\oent[A_{i_1}^{[1]}X_{u_2}^{[2]}\cfots
A_{i_N}^{[N]}\right]|i_1,i_2,\ldots,i_N\rangle\;,$$ and are called *mwfrox product statef* [@MPSorigin]. Ad xrown in [@Vi03] *evgry* stacs dan be represented in thif way if only uhe bpnd dimensions $D_k$ are suffixiently largv. Hebce, Eq.(\[MPS0\]) is a repxesentation uf sjates tather than the charactzrizatjon of a spfcific clzrs. However, typicxllj stdtes are referred to as MPF if they havz a MPS-rdprexentatyon with slall $D=\max_k D_k$ which (in tje cade of a sequfnce of states) does in particulec not grow wijh $T$. Nmte that $\psi$ ln Eq.(\[MPS0\]) is in deneral not notmalized cnd thxt its MPS
$ to each of them. Assume that of virtual spins correspond to different (unnormalized) entangled state $|I\rangle=\sum_{\alpha=1}^{D} often referred to entangled *bond*. Then apply a map A}=\sum_{i=1}^d\sum_{\alpha,\beta=1}^{D} A_{i,\alpha,\beta} |i\rangle\langle\alpha,\beta|$$ to each of the $N$ sites. Here and in the Greek indices correspond to the virtual systems. By writing $A_{i}$ for the $D\times matrix elements we that the coefficients of the final state when expressed in terms of a product basis are by a matrix product $\operatorname{tr}\left[A_{i_1}A_{i_2}\cdots A_{i_N}\right]$. In general dimension of the entangled $|I\rangle$ and the map ${\cal can be site-dependent we $A_i^{[k]}$ the $D_k\times D_{k+1}$ corresponding to site $k\in\{1,\ldots,N\}$. States obtained in this way have then the form $$\label{MPS0} |\psi\rangle=\sum_{i_1,\ldots,i_N=1}^d \operatorname{tr}\left[A_{i_1}^{[1]}A_{i_2}^{[2]}\cdots A_{i_N}^{[N]}\right]|i_1,i_2,\ldots,i_N\rangle\;,$$ are called states* [@MPSorigin]. shown [@Vi03] state can be this way if only the bond sufficiently large. Hence, Eq.(\[MPS0\]) is a representation of rather than characterization of a specific class. However, states are referred to as MPS if they a MPS-representation with small $D=\max_k D_k$ which (in the case of a sequence of states) particular not grow with Note that $\psi$ Eq.(\[MPS0\]) in not and that MPS
$ to each of them. Assume that eveRy pair of neIghboRinG viRtUal sPins Which corresponD To diFferent sites are initialLy in aN (uNNormALiZed) maXimally ENtANGleD sTaTe $|I\RaNGlE=\sum_{\aLphA=1}^{D} |\alpha,\Alpha\ranglE$ ofTeN referred to aS EnTangled *bonD*. ThEn apply a map $$\lAbeL{curlya}{\cAl
A}=\SUm_{i=1}^d\sUm_{\aLpha,\bEta=1}^{D} A_{i,\ALpha,\beTa}
|i\rangle\LaNGle\alpHA,\beta|$$ to EACh Of thE $N$ sites. Here and in tHE fOLlowing Greek inDices cOrREsPONd tO thE virtual syStEms. By WRiting $A_{I}$ FoR THE $D\tIMes D$ matrix witH elements $A_{i,\ALphA,\beta}$ wE gEt tHAt the cOeffiCiENts Of the final sTate When expreSsed in TErms of a PRoduct bAsis arE giVen By a mATrIx ProDuCT $\opERaTorNAme{Tr}\left[A_{i_1}a_{i_2}\CdOts
A_{i_n}\rigHT]$. iN GeneRal The dImensIon of the entanGleD staTE $|I\rAngle$ And thE map ${\CaL A}$ can Both be Site-dEpEndent and we writE $A_i^{[k]}$ For the $D_k\tImeS D_{K+1}$ maTrIx corREspondIng To sIte $k\in\{1,\lDots,N\}$. StATes ObTAINeD in this way have then ThE FOrM $$\label{MPs0}
|\psi\raNGlE=\sUM_{i_1,\ldots,i_n=1}^d
\OpeRatoRNAme{tr}\Left[a_{I_1}^{[1]}A_{I_2}^{[2]}\cdots
A_{i_n}^{[N]}\righT]|I_1,i_2,\LdOts,i_N\raNgLe\;,$$ and aRe CalLed *MatriX ProdUct staTes* [@MPSorIgin]. AS Shown in [@Vi03] *every* STate can be reprESeNTEd IN thiS waY if only the bOnd dIMensIons $d_K$ aRe sUFficiEntly LaRGe. hEnce, Eq.(\[MPS0\]) is a represeNtAtion oF statEs rather than tHe characteRIZAtion of a SpecIFiC Class. However, tyPicalLy states arE Referred To as MpS if they Have a MPS-rEPResentatIon WitH smAll $d=\MAx_K D_k$ which (in the CASe of A sEquence Of sTates) doEs iN paRtiCulAr Not grow wiTh $N$. Note tHaT $\pSi$ In eq.(\[MpS0\]) is iN General nOt NorMaLizEd and THat its mPS
$ to each of them. Assumethat every  pair  of  ne ig hbor ingvirtual spinsw hich  correspond to differe nt si te s  are in itial ly in a n  ( u n nor ma li zed )m ax imall y e ntangle d state $| I\r an gle=\sum_{\a l ph a=1}^{D} | \al pha,\alpha\r ang le$ of te n r e ferre d t o asentang l ed *bo nd*. Then  a p ply am ap $$\l a b el {cur lyA}{\cal
A}=\sum _ {i = 1}^d\sum_{\alp ha,\be ta = 1} ^ { D}A_{ i,\alpha,\ be ta}
| i \rangle \ la n g l e\a l pha,\beta|$$to each oft he$N$ si te s.H ere an d inth e  fo llowing Gre ek i ndices co rrespo n d to th e  virtua l syst ems . B y wr i ti ng  $A _{ i }$f or  th e  $D \times D $ma trixwith e l e ment s $ A_{i ,\alp ha,\beta}$ we  ge t th a t t he co effic ient sof th e fina l sta te  when expressed  interms ofa p ro duc tbasis are gi ven  by  a matr ix prod u ct$\ o p e ra torname{tr}\left[A _{ i _ 1} A_{i_2}\ cdotsA _{ i_ N }\right] $.  In  gen e r al th e di m en sion ofthe en t an gl ed stat e$|I\ra ng le$  an d the map${\cal  A}$ can  both be site-depend e nt and we wri t e$ A _i ^ {[k] }$for the $D_ k\ti m es D _{k+ 1 }$  ma t rix c orres po n di n g to site $k\in\{1, \l dots,N \}$.States obtain ed in this w a y have t hent he form $$\label{ MPS0} 
|\psi\ran g le=\sum_ {i_1, \ldots,i _N=1}^d
\ o p eratorna me{ tr} \le ft[ A _ {i _1}^{[1]}A_{i _ 2 }^{[ 2] }\cdots 
A_ {i_N}^{ [N] }\r igh t]| i_ 1,i_2,\ld ots,i_N\ ra ng le \; ,$$  anda re calle d*ma tr ixprodu c t stat es* [ @MPS or ig i n].  As sho w ni n  [@V i0 3]  *ev ery *state  can bereprese nted in t his wayif  o nly the  bond dimensi on s $D_k$ ar esuf ficien t l y large.  Hence, Eq.(\[MPS0\]) i s  a repr ese ntati on o f statesrat her th ant he cha racter izati on  of a  spec i f ic  cl as s. However , typ icall ystat es arereferred to as MPS ifthey have a M PS- repr e s en tat i on wit hs mal l $D=\max_k D_k$which (inth e  c ase of a s e que nc e of st ates) d oes i n  partic ular notgrow with  $ N$.N o tethat $\psi $ in Eq. (\[MPS0\] )  is i n  g enera l n ot nor ma liz ed an d that its  MPS
$ to_each of_them. Assume that every_pair of_neighboring_virtual spins_which_correspond to different_sites are initially_in an (unnormalized) maximally_entangled state $|I\rangle=\sum_{\alpha=1}^{D}_|\alpha,\alpha\rangle$_often referred to as entangled *bond*. Then apply a map $$\label{curlyA}{\cal
A}=\sum_{i=1}^d\sum_{\alpha,\beta=1}^{D} A_{i,\alpha,\beta}
|i\rangle\langle\alpha,\beta|$$ to each_of_the $N$_sites._Here_and in the following Greek_indices correspond to the virtual_systems. By_writing $A_{i}$ for the $D\times D$ matrix with_elements_$A_{i,\alpha,\beta}$ we get_that the coefficients of the final state when expressed_in terms of a product basis_are given by_a_matrix_product $\operatorname{tr}\left[A_{i_1}A_{i_2}\cdots
A_{i_N}\right]$. In general_the dimension of the entangled state_$|I\rangle$ and the map ${\cal A}$_can both be site-dependent and we write_$A_i^{[k]}$ for the $D_k\times D_{k+1}$ matrix_corresponding to site $k\in\{1,\ldots,N\}$. States_obtained in_this way have then the_form $$\label{MPS0}
|\psi\rangle=\sum_{i_1,\ldots,i_N=1}^d
\operatorname{tr}\left[A_{i_1}^{[1]}A_{i_2}^{[2]}\cdots
A_{i_N}^{[N]}\right]|i_1,i_2,\ldots,i_N\rangle\;,$$ and_are called_*matrix product states*_[@MPSorigin]. As shown in [@Vi03] *every*_state can be_represented in this way if only_the_bond dimensions $D_k$_are_sufficiently_large. Hence,_Eq.(\[MPS0\]) is a_representation_of states_rather_than the characterization of a specific_class._However, typically states are referred to as_MPS if they have_a_MPS-representation with small $D=\max_k_D_k$ which (in the case_of a sequence of states) does_in particular_not grow_with $N$. Note that $\psi$ in Eq.(\[MPS0\]) is in general not_normalized and that its MPS