16 Matrix Product Operators and Density Operators
Fix a physical dimension \(d\) and a bond dimension \(D\). An MPO tensor carries two physical indices, one ket index and one bra index, together with the virtual left and right indices. Matrix product density operators (MPDOs) and matrix product operators (MPOs) describe 1D mixed states and quantum operators with bounded virtual entanglement [ 13 , 14 , 17 ] . The renormalization fixed-point conditions and their algebraic consequences are developed in Chapters 20 and 21.
16.1 MPO tensors
An MPO tensor is a family of matrices \(\{ M^{ij}\} _{i,j=0}^{d-1}\) with \(M^{ij} \in M_{D}(\mathbb {C})\), indexed by a ket index \(i\) and a bra index \(j\). Diagrammatically,
By combining the pair \((i,j)\) into a single physical index in \(\{ 0,\ldots ,d^2-1\} \), one obtains an MPS tensor with physical dimension \(d^2\) and bond dimension \(D\), with \(M^{\mathrm{MPS}}_{(i,j)}=M^{ij}\), identifying \(\{ 0,\ldots ,d{-}1\} \times \{ 0,\ldots ,d{-}1\} \) with \(\{ 0,\ldots ,d^2{-}1\} \) via the product encoding.
For words \(u=(i_1,\ldots ,i_N)\) and \(v=(j_1,\ldots ,j_N)\) in \(\{ 0,\ldots ,d{-}1\} ^N\), the word evaluation is
The empty pair of words evaluates to the identity, \(M^{\varnothing ,\varnothing }=\mathbb {1}_D\), and word pairs of different lengths evaluate to \(0\).
For words \(u_1,v_1\) of equal length and arbitrary words \(u_2,v_2\), word evaluation factors as \(M^{u_1u_2,\, v_1v_2}=M^{u_1,v_1}M^{u_2,v_2}\).
We argue by induction on the common length of \(u_1\) and \(v_1\). For the empty prefix, \(M^{u_2,v_2}=\mathbb {1}_D M^{u_2,v_2}\). If \(u_1=iu\) and \(v_1=jv\), then
using associativity of matrix multiplication.
Let \(a,b\in \{ 0,\ldots ,d{-}1\} \), and let \(u,v\) be words of equal length. Writing \(au,bv\) for adding the letters at the beginning and \(ua,vb\) for adding them at the end, one has \(\operatorname{tr}\! \left(M^{au,\, bv}\right)=\operatorname{tr}\! \left(M^{ua,\, vb}\right)\).
Lemma 16.1.4 gives \(M^{ua,\, vb}=M^{u,v}M^{ab}\) and \(M^{au,\, bv}=M^{ab}M^{u,v}\). Hence, by cyclicity of the trace,
The operator generated by \(M\) on \(N\) sites is the matrix \(\rho ^{(N)}(M)\) indexed by configurations \(\sigma ,\tau \in \{ 0,\ldots ,d{-}1\} ^N\) and given by
Diagrammatically,
Let \(M\) be a matrix product operator tensor with periodic operators \(O_L\), and let \(c\in \mathbb {C}^d\). The tensor with letters \(c_jM^{ij}\) has periodic operators \(O_L\, \operatorname {diag}(c)^{\otimes L}\). If every \(c_j\) is nonzero, its doubled-index tensor is normal exactly when that of \(M\) is.
The trace of \(\prod _n c_{t_n}M^{s_nt_n}\) is \(\bigl(\prod _nc_{t_n}\bigr)\operatorname {tr}\prod _nM^{s_nt_n}\). Every word of the doubled-index tensor is multiplied by the product of the scalars of its letters, which is nonzero, so the span of the words of each length is unchanged.
The direct sum of two MPO tensors is block diagonal on the bond space \(\mathbb {C}^{D_1}\oplus \mathbb {C}^{D_2}\): every letter is the block-diagonal sum of the two letters,
Word evaluation of the direct sum is the block-diagonal sum of the word evaluations, \((M\oplus N)^{u,v}=M^{u,v}\oplus N^{u,v}\).
Induction on the common length of \(u,v\), expanding one letter at a time and using that block-diagonal sums multiply blockwise.
The operator generated by the direct sum is the sum of the two operators, at every system size: \(\rho ^{(N)}(M\oplus N)=\rho ^{(N)}(M)+\rho ^{(N)}(N)\).
The trace of a block-diagonal sum is the sum of the traces of the two blocks; apply Lemma 16.1.9 and take traces entrywise.
The cyclic shift of configurations on \(N\) sites is the bijection that relabels each site by one position, \(\sigma \mapsto \sigma \circ r\), where \(r\) is the cyclic rotation \(i\mapsto i+1\) of \(\{ 0,\ldots ,N-1\} \).
The operator \(\rho ^{(N)}(M)\) is invariant under the simultaneous cyclic shift of its bra and ket configurations: \(\rho ^{(N)}(M)_{r\sigma ,\, r\tau }=\rho ^{(N)}(M)_{\sigma ,\tau }\).
The entry is a closed bond trace, \(\rho ^{(N)}(M)_{\sigma ,\tau } =\operatorname{tr}(M^{\sigma _1\tau _1}\cdots M^{\sigma _N\tau _N})\). Shifting both configurations by \(r\) rotates the matrix product by one factor, and the trace is invariant under that cyclic rotation:
so \(\rho ^{(N)}(M)_{r\sigma ,\, r\tau } =\rho ^{(N)}(M)_{\sigma ,\tau }\).
An MPO tensor is Hermitian if \(M^{ij}=(M^{ji})^\dagger \) for all \(i,j\in \{ 0,\ldots ,d{-}1\} \).
The adjoint tensor \(M^\dagger \) of an MPO tensor \(M\) exchanges the ket and bra physical indices and conjugate-transposes each virtual matrix: \((M^\dagger )^{ij}=(M^{ji})^\dagger \). A tensor equals its adjoint tensor, \(M^\dagger =M\), precisely when it is Hermitian.
The adjoint tensor generates the adjoint operator family up to spatial reflection: entrywise,
where \(\mathrm{rev}\) reverses the site order. For an MPDO and \(N\geq 1\), the conjugation disappears:
and the adjoint tensor is again an MPDO.
At the level of word evaluations, the conjugate transpose of a product reverses the factors, so
Taking traces gives the entrywise identity
For an MPDO, \(\rho ^{(N)}(M)\) is Hermitian, so
positive semidefiniteness is preserved under the simultaneous relabeling of rows and columns by the reflection.
16.2 Transfer maps
The transfer map associated to \(M\) is the linear map \(\mathcal{E}_M:M_{D}(\mathbb {C})\to M_{D}(\mathbb {C})\) defined by
The transfer map of an MPO tensor agrees with the transfer map of its doubled-index MPS view: \(\mathcal{E}_M=\mathcal{E}_{M^{\mathrm{MPS}}}\).
After identifying the single physical index of the doubled tensor with a pair \((i,j)\), for every \(X\in M_{D}(\mathbb {C})\),
If \(X\ge 0\), then \(\mathcal{E}_M(X)\ge 0\).
By Lemma 16.2.2, \(\mathcal{E}_M\) is the transfer map of the doubled-index MPS tensor. The corresponding MPS transfer map is positive, so \(\mathcal{E}_M(X)\ge 0\) whenever \(X\ge 0\).
16.3 MPDOs and LPDOs
An MPO tensor is an MPDO if, for every \(N\geq 1\),
This is the global positivity condition of [ 13 , Section 4 ] .
An MPO tensor is an LPDO if there exist a Kraus dimension \(d_K\), an inner bond dimension \(D'\), a bond-space identification \(e:\{ 0,\ldots ,D-1\} \simeq \{ 0,\ldots ,D'-1\} \times \{ 0,\ldots ,D'-1\} \) (equivalently, \(D\) and \((D')^2\) have the same cardinality), and matrices \(A^{(i,k)}\in M_{D'}(\mathbb {C})\) such that, for all physical indices \(i,j\),
where \(\otimes \) denotes the Kronecker product, \(\overline{(\cdot )}\) denotes entrywise complex conjugation, and \((\cdot )_{e,e}\) denotes reindexing along \(e\) back to \(\{ 0,\ldots ,D-1\} \). See [ 13 , Section 4.3 ] . The local contraction is
Every LPDO tensor is an MPDO.
By the mixed-product property \((A\otimes B)(C\otimes D)=(AC)\otimes (BD)\), the \(N\)-site product decomposes as a sum of Kronecker products after transporting along the bond-space identification \(e\):
where \(P_\kappa =A^{(\sigma _1,\kappa _1)}\cdots A^{(\sigma _N,\kappa _N)}\) and \(Q_\kappa =A^{(\tau _1,\kappa _1)}\cdots A^{(\tau _N,\kappa _N)}\). Since the trace is invariant under the reindexing by \(e\), we apply \(\operatorname{tr}(X\otimes Y)=\operatorname{tr}(X)\operatorname{tr}(Y)\) to obtain
where \(\psi _\kappa (\sigma )=\operatorname{tr}(P_\kappa )\). Hence \(\rho ^{(N)}(M) =\sum _\kappa |\psi _\kappa \rangle \! \langle \psi _\kappa |\ge 0\).
16.4 Horizontal and Vertical Canonical Forms
An MPO tensor is in normalized BNT-refined horizontal form if its doubled-index MPS tensor is a BNT sector decomposition with an independent gauge on every repeated copy:
Here the \(A_j\) are the distinct minimal representatives, \(r_j{\gt}0\), and every weight \(\mu _{j,q}\) is nonzero, \(|\mu _{j,q}|\leq 1\), and at least one weight has modulus one. The BNT hypotheses include irreducibility, left-canonicality, normalized self-overlap, eventual linear independence of the representative MPV families, and the exclusion of similarity up to phase between distinct representatives. Repeated gauge-equivalent copies are retained over the same representative rather than separated as distinct blocks. The total bond dimension of the displayed direct sum is exactly \(D\). Equivalently, with the block-diagonal matrix \(X=\bigoplus _{j,q}X_{j,q}\),
This is the strengthened, normalized BNT-refined form of the canonical decomposition in [ 13 , Section 2.3 ] : unlike the literal predicate (Definition 7.6.7), the direct sum exhausts the full bond dimension \(D\), each copy weight is nonzero and normalized, and gauge-equivalent normal blocks are grouped over a minimal BNT representative with an independent gauge for each copy.
An MPO tensor has a horizontal-canonical-form MPV representation if its doubled-index MPS tensor generates the same positive-length MPV family as a BNT canonical form \(\bigoplus _{j=0}^{g-1}\bigoplus _{q=0}^{r_j-1}\mu _{j,q}A_j\). Its length-\(N\) matrix product vector is
This is the MPV-level decomposition in [ 13 , Section 2.3, lines 298–301 ] .
Every MPO tensor in normalized BNT-refined horizontal form has the associated positive-length BNT matrix-product-vector representation.
Each copy similarity cancels inside every positive-length closed matrix-product trace.
The vertical direction of an MPO tensor exchanges the notion of physical and virtual indices [ 13 , line 943 ] : the vertically viewed tensor \(\widetilde M\) has its letters indexed by the horizontal bond pair \((a,b)\), and each letter is the operator on the physical space with matrix elements \((\widetilde M_{ab})_{ij}=M^{ij}_{ab}\). Concatenating \(\widetilde M\) vertically generates the matrix product vectors of the vertical direction.
Adjunction exchanges the two oriented horizontal bond indices in the vertical view:
This follows entrywise from the exchange of the ket and bra indices and conjugate transposition of the virtual matrix.
An MPO tensor is in vertical canonical form if there are a basis of normal tensors \(\{ M_\alpha \} _\alpha \) for its vertically viewed tensor \(\widetilde M\), positive diagonal matrices \(\mu _\alpha =\operatorname{diag}(\omega _{\alpha ,0},\ldots ,\omega _{\alpha ,r_\alpha -1})\) with \(r_\alpha \geq 1\), and a matrix \(U\) from the physical space onto the retained nonzero sector space such that
This is the conclusion of [ 13 , Proposition 4.13, lines 1863–1870 ] . In this orientation \(UU^\dagger =I\) on the retained space, while \(U^\dagger U\) is its support projection in the original physical space. The reconstruction identity says that every letter of \(\widetilde M\) is supported on this projection, so the discarded orthogonal complement and the off-diagonal corners vanish. The orthogonal complement may be a zero sector: the source’s general canonical-form construction explicitly permits the sum of the nonzero block dimensions to be strictly smaller than the original dimension [ 13 , lines 214–225 ] . Since each \(\mu _\alpha \) is diagonal, the summand \(\mu _\alpha \otimes M_\alpha \) is the direct sum \(\bigoplus _{k=0}^{r_\alpha -1}\omega _{\alpha ,k}M_\alpha \) of weighted copies of \(M_\alpha \), so the right-hand side is a block-diagonal tensor over the pairs \((\alpha ,k)\) in which the multiplicities \(r_\alpha \) and the diagonal entries \(\omega _{\alpha ,k}\) stay explicit; the coefficients \(m_\alpha =\operatorname{tr}(\mu _\alpha )=\sum _{k=0}^{r_\alpha -1}\omega _{\alpha ,k}\) appearing in the renormalization fixed-point characterization are recovered from them.
Let \((A_k)_{k}\) be a finite family of block tensors of physical dimension \(d\), and fix a length \(L\); write \(A_k^w\) for the length-\(L\) word evaluation of \(A_k\) along a word \(w \in \{ 0,\ldots ,d{-}1\} ^L\). The finite-length spanning hypothesis is that the tuple-valued evaluations span the full product algebra of block matrices,
The corresponding finite-length separation property is that some length \(L\) makes the trace pairing against block words nondegenerate: there is a length \(L\) such that, for every tuple \((\Delta _k)_k\) of block matrices,
The spanning and separation conditions are equivalent after existentially quantifying the length. Spanning at length \(L\) gives separation at that length by nondegeneracy of the product trace pairing. Conversely, if the span were proper, a nonzero linear functional \(f\) would annihilate it. Nondegeneracy of the matrix trace pairing represents \(f\) as
for some tuple \((\Delta _k)_k\), contradicting separation. This finite-word property generalizes the one-letter simultaneous span in the source definition of block-injective canonical form to length-\(L\) words.
The spanning hypothesis has an equivalent entrywise form. For a block index \(k\) and matrix entry \((a,b)\), let \(w\longmapsto (A_k^w)_{a,b}\) be the corresponding scalar word family, indexed jointly by the triple \((k,a,b)\). Linear independence of this scalar family implies the product-algebra spanning condition, and hence trace separation.
The block-injectivity proposition of [ 13 , lines 340–345 ] gives a concrete sufficient criterion for the spanning condition: a common length \(L\) with every block \(L\)-block-injective, together with a length \(S\) at which, for every \(k\), there is a coefficient family \((c_w)_{w\in \{ 0,\ldots ,d{-}1\} ^S}\) with
Concatenating the injective prefix with this finite selector family produces the spanning condition at length \(L+S\), hence the product-algebra spanning and trace-separation properties.
The inverse tensor of a block-injective canonical form satisfies
Writing \(K=X(\bigoplus _j\mathcal K_j)X^{-1}\) in this contraction places the gauges on the virtual legs of \(K\) and identifies the block index of \(\mathcal K_j\) with the label \(j\):
Consequently, \(K^{-1}K=\bigoplus _j\mathbb {1}_j\). This is the inverse identity in [ 13 , Appendix C.2, lines 1666–1670 ] .
If a block family has a common injective word length and finite words that select each block from the others, then the family has finite-word trace separation.
If the common injective prefix has length \(L\) and the selector has length \(S\), concatenation gives simultaneous words of length \(L+S\). For every tuple \((X_j)_j\), the resulting coefficients satisfy \(\sum _{|w|=L+S}c_w A_j^w=X_j\). Hence the simultaneous word tuples span the product algebra, and nondegeneracy of the trace pairing gives finite-length trace separation.
Let \((A_j)_{j=0}^{g-1}\) be a family with \(g=1\). If its unique block is \(L\)-block injective, then the simultaneous length-\(L\) word evaluations span the full one-block product algebra.
For every matrix \(X\), block injectivity gives coefficients such that \(X=\sum _{|w|=L}c_wA_0^w\). Since the family has only the block \(A_0\), this is precisely the required simultaneous word-tuple spanning identity.
Let \((A_j)_{j=0}^{g-1}\) be the basis representatives of a sector decomposition. They have representative word-tuple separation at length \(L\) when
They have eventual representative word-tuple separation when this equality holds at every sufficiently large length. Notice that the family is indexed by the minimal representatives \(j\), and not by the repeated copies lying over each \(j\).
Let \(A_{\mathrm{tot}}\) be the tensor of a sector decomposition with basis representatives \((A_j)_{j=0}^{g-1}\), multiplicities \(r_j{\gt}0\), nonzero weights \(\mu _{j,q}\), and grouped coefficients \(c_j^{(N)}=\sum _{q=0}^{r_j-1}\mu _{j,q}^{N}\). Suppose that the first-site actions of \(Y,Z\in M_{d}(\mathbb {C})\) agree on the MPV family of \(A_{\mathrm{tot}}\). If representative word-tuple separation holds at length \(L\) and \(c_{j_0}^{(L+1)}\neq 0\), then \(YA_{j_0}=ZA_{j_0}\). Consequently, under eventual representative word-tuple separation, \(YA_j=ZA_j\) for every representative \(j\). The same conclusion holds if the first-site equality is initially given for any tensor with the same positive-length MPV family as \(A_{\mathrm{tot}}\).
The representative grouping is the argument of [ 13 , Appendix C.3, lines 1835–1858 ] .
In the first auxiliary assertion, eventual representative word-tuple separation is an explicit assumption. The later post-blocked theorem derives the needed finite representative separation from the full BNT canonical-form hypotheses. The sharp \(3D^5\) representative-separation bound is Theorem 8.6.1.24; its detailed source comparison is recorded in [ 108 ] .
At a fixed tail length \(L\), expand the first-site equality at system size \(L+1\) and group all copies lying over the same representative. For every length-\(L\) word \(w\) this gives
Representative word-tuple separation therefore forces each matrix \(c_j^{(L+1)}(YA_j-ZA_j)\) to vanish. A nonzero coefficient may then be cancelled.
For the eventual assertion, choose a threshold \(L_0\) and put \(M=L_0+1\). For a fixed representative \(j\), apply Corollary 8.5.2 to the nonzero numbers \(\mu _{j,q}^{M}\). It yields \(r\geq 1\) such that \(c_j^{(Mr)}=\sum _q(\mu _{j,q}^{M})^r\neq 0\). Since \(Mr-1\geq L_0\), the fixed-length argument applies with \(L=Mr-1\). Finally, equality of the positive-length MPV families gives the first-site hypothesis for the sector-decomposition tensor.
Let
be the tensor of a sector decomposition. Let \(s\) range over a second finite alphabet, and for each representative \(A_j\) let \(C_j^s\) be a matrix on its bond space. The corresponding direct sum with one marked letter is
The marked alphabet need not coincide with the alphabet of the unmarked tensor. This is the algebraic form of the marked matrices in the two diagrams of [ 13 , Proposition 4.13, lines 1909–1919 ] .
For every marked index \(s\) and every unmarked word \(w\),
The exponent \(|w|+1\) includes the weight at the marked site.
Expand both matrices as direct sums. Their product is the direct sum of \(\mu _{j,q}^{|w|+1}C_j^sA_j^w\). Taking the trace and summing first over \(q\) gives the stated coefficient \(c_j^{(|w|+1)}\).
Let \((C_j^s)\) and \((E_j^s)\) be two marked families over the same sector decomposition. Fix a tail length \(L\), and suppose that, for every marked index \(s\) and every unmarked word \(w\) of length \(L\),
If representative word-tuple separation holds at length \(L\) and \(c_{j_0}^{(L+1)}\neq 0\), then \(C_{j_0}=E_{j_0}\). Under eventual representative word-tuple separation, agreement for every marked index and every positive-length unmarked word implies \(C_j=E_j\) for every \(j\).
This extends the argument of [ 13 , Appendix C.3, lines 1848–1858 ] from first-site physical marks to general marked matrices in the marked slot.
Fix \(s\) and a tail length \(L\), and put \(\Delta _j=c_j^{(L+1)}(C_j^s-E_j^s)\). Lemma 16.4.13 gives \(\sum _{j=0}^{g-1}\operatorname{tr}(\Delta _jA_j^w)=0\) for every length-\(L\) word \(w\). Representative word-tuple separation implies \(\Delta _j=0\) for every \(j\). If \(c_{j_0}^{(L+1)}\neq 0\), cancellation gives \(C_{j_0}^s=E_{j_0}^s\). Since \(s\) was arbitrary, \(C_{j_0}=E_{j_0}\).
Under eventual separation, apply Corollary 8.5.2 to the nonzero weights above each fixed representative, exactly as in Theorem 16.4.11. This selects a sufficiently large \(L\) with \(c_j^{(L+1)}\neq 0\) and proves \(C_j=E_j\).
Let \(T_{\mathrm{grp}}\) be the full grouped tensor of Theorem 8.6.1.9, and let \(P_{\mathrm{rep}}\) be its representative sector decomposition. For a marked family \((C_j^s)_j\), define the marked tensor \(\widehat C_{\mathrm{grp}}\) on all grouped listed coordinates as follows. The copy \((j,q)\) carries the marked block \(\lambda _{j,q}C_j^s\), where \(\lambda _{j,q}=\nu _{k(j,q)}\zeta _{k(j,q)}\).
For every marked letter \(s\) and every unmarked word \(w\), one has
The exponent \(|w|+1\) includes the factor at the marked site.
Reindex the finite sum over listed coordinates by the class-copy pairs \((j,q)\) and use the bond-dimension identifications to replace each block by its unit-phase representative. The marked block supplies one factor \(\lambda _{j,q}\) and the word supplies \(|w|\) further factors. Summing first over \(q\) gives (20).
Let \((C_j^s)_j\) and \((E_j^s)_j\) be two marked families on the representatives. Form their full-coordinate marked tensors. If
for every \(s\) and every word \(w\) of positive length, then \(C_j=E_j\) for every representative \(j\).
In particular, let \(Y,Z\in M_{d}(\mathbb {C})\). If their first-site actions agree on the positive-length matrix product vectors of \(T_{\mathrm{grp}}\), then
for every representative \(j\) and every physical letter \(i\).
Theorem 16.4.15 rewrites (21) into marked-trace agreement for \(P_{\mathrm{rep}}\). Its representatives have eventual word-tuple separation, and every copy weight is nonzero. Theorem 16.4.14 therefore gives \(C_j=E_j\).
First-site agreement for \(T_{\mathrm{grp}}\) also holds for \(P_{\mathrm{rep}}\) because their matrix product vectors agree at every positive length. Applying Theorem 16.4.11 gives (22), namely
Let \(A\) be in literal CPSV canonical form. For canonical-form data and its BNT refinement, write
where \(B^i=GT_{\mathrm{grp}}^iG^{-1}\) and \(G\) is the listed block gauge. If marked matrices satisfy \(C=V^\dagger EV\), then for every word \(w\), including the empty word,
The inserted marked matrix ensures that the closed chain is nonempty.
For a linear physical mark \(f\), the linear marking of \(T_{\mathrm{grp}}\) is exactly its full grouped marked tensor. Each copy \((j,q)\) carries \(\lambda _{j,q}f(C_j)\), where \(\lambda _{j,q}=\nu _{k(j,q)}\zeta _{k(j,q)}\). For two linear marks \(f\) and \(g\), put \(F_A^u=\sum _i f_{ui}A^i\) and \(G_A^u=\sum _i g_{ui}A^i\). If
for every marked letter \(u\) and every word \(w\) of positive length, then \(F_A=G_A\) on the original bond space.
This positive-tail formulation corresponds to the argument in [ 13 , Appendix C.3, lines 1835–1858 ] : first-site agreement at positive system sizes suffices to separate representatives at a nonzero tail length without requiring empty-word data.
Induction on \(w\) in (24) gives \(CA^w=V^\dagger EB^wV\). Cyclicity of trace and \(VV^\dagger =\mathbb {1}\) give (25). Linearity of the physical mark gives the grouped marked tensor.
Apply (25) to both sides of (26), remove the block gauge by cyclicity of trace, and use Theorem 16.4.16 on the representatives. The resulting equality of the full grouped marked tensors, followed by (24), gives \(F_A=G_A\).
Let \(Y\in M_{d}(\mathbb {C})\), and write physical indices of an \((L+1)\)-site block as \(I=(i_0,\ldots ,i_L)\) and \(J=(j_0,\ldots ,j_L)\). The induced matrix on the blocked physical space is
This identity connects the block-injectivity argument [ 13 , lines 318–345 ] with Lemma L [ 13 , Appendix C.3, lines 1835–1858 ] .
Let \(A\) be a tensor, and let \(Y^{[L+1]}\) be the physical-block action induced by \(Y\in M_{d}(\mathbb {C})\). For \(I=(i_0,\ldots ,i_L)\), its insertion into the blocked tensor satisfies
Hence equality of the first-site actions of \(Y,Z\) on the MPV family of \(A\) implies equality of the first-block actions of \(Y^{[L+1]},Z^{[L+1]}\) on the MPV family of \(A^{[L+1]}\).
The first-site actions of \(Y,Z\) agree at every length if and only if
for every physical index \(i\) and every finite word \(w\).
Conversely, suppose that \(A\) is \(L\)-block injective. If the two induced insertion tensors on \(A^{[L+1]}\) are equal, then \(YA=ZA\).
For the forward implication, split the first block into its first letter and its remaining \(L\) letters. The insertion identity (28) reduces the blocked MPV equality to the original first-site equality, with the remaining letters of the first block prepended to the flattened tail.
For the converse, put \(\Delta ^i=(YA)^i-(ZA)^i\). Equality after blocking gives \(\Delta ^i A^w=0\) for every word \(w\) of length \(L\). Since the length-\(L\) words span \(M_{D}(\mathbb {C})\), the same equality holds with \(A^w\) replaced by the identity, and therefore \(\Delta ^i=0\) for every physical index \(i\).
Let \(P\) be a BNT canonical form with basis representatives \((A_j)_{j=0}^{g-1}\), and suppose that the representatives have already been blocked so that every \(A_j\) is injective. Then representative word-tuple separation holds at every length \(L\geq 6(g-1)+1\). Consequently, if the first-site actions of \(Y,Z\in M_{d}(\mathbb {C})\) agree on the positive-length MPV family of the tensor represented by \(P\), then \(YA_j=ZA_j\) for \(0\leq j{\lt}g\). The same conclusion holds when the first-site equality is given for any tensor with the same positive-length MPV family as the tensor represented by \(P\).
More generally, let \(P\) be a BNT canonical form before physical blocking, let \(p{\gt}0\), and suppose that every \(p\)-blocked representative \(A_j^{[p]}\) is injective. Then the blocked sector decomposition \(P^{[p]}\) has representative word-tuple separation at every blocked length \(L\geq 6(g-1)+1\).
Without changing the physical alphabet, a BNT canonical form also has eventual representative word-tuple separation. If \(p{\gt}0\) is a common block-injectivity length, then one may take pairwise separating words of length \(6p\), form selectors of length \(6p(g-1)\), and append homogeneous words of any length at least \(p\).
Finally, the BNT canonical-form hypotheses imply that there is a common \(L{\gt}0\) for which every unblocked representative \(A_j\) is \(L\)-block injective. Indeed, irreducibility and left-canonicality give a finite peripheral period \(m_j\). The self-overlap tends to \(m_j\) along multiples of \(m_j\), but by normalization it also tends to \(1\); hence \(m_j=1\). Thus every representative is normal, and finiteness gives a common positive block-injectivity length. If the first-site actions of \(Y,Z\) agree on the tensor represented by \(P\), then \(YA_j=ZA_j\) for \(0\leq j{\lt}g\). The same conclusion holds when the first-site equality is given for any tensor with the same positive-length MPV family. This composes the blocking conclusion of [ 13 , lines 318–345 ] with Appendix C.3, Lemma L [ 13 , lines 1835–1858 ] .
No additional common-length hypothesis is required.
Injectivity of \(A_j\) implies that its homogeneous words span \(M_{D_j}(\mathbb {C})\) at every positive length. Apply the three-block separation argument with source length \(1\). For each ordered pair of distinct representatives it gives a pairwise block-separating equation at length \(6\). Multiplying the equations for the \(g-1\) representatives different from \(A_j\) gives coefficients \(c_w^{(j)}\) satisfying
If \(L\geq 6(g-1)+1\), put \(m=L-6(g-1){\gt}0\). The length-\(m\) words of each \(A_j\) span \(M_{D_j}(\mathbb {C})\). Concatenating these words with the selectors in (32) shows that the simultaneous length-\(L\) word tuples span \(\prod _jM_{D_j}(\mathbb {C})\). This proves eventual representative word-tuple separation. The representative-grouped Lemma L gives the two first-site conclusions.
For the unblocked canonical form, propagate \(p\)-block injectivity to the lengths \(2p-1\), \(2p\), and \(6p\). Applying the direct-sum argument to the resulting injective tensors gives pair trace separation at original length \(6p\). The bijection between blocked words of length \(6\) and original words of length \(6p\) gives pair trace separation for the \(p\)-blocked representatives. The preceding selector argument then applies to \(P^{[p]}\).
The same pairwise equations at original length \(6p\) give selectors for the unblocked representatives directly. For every \(L\geq 6p(g-1)+p\), the remaining prefix has length at least \(p\) and hence spans each full representative matrix algebra. Concatenating the prefix with the selectors proves eventual representative word-tuple separation before blocking.
For the final assertion, first consider one representative \(A_j\). Irreducibility and left-canonicality give a positive peripheral period \(m_j\). The periodic self-overlap formula and the normalized self-overlap hypothesis give, along the same subsequence,
Hence \(m_j=1\). Primitivity then implies normality. Taking a common multiple of the finitely many normality lengths gives a positive length \(L\) at which every representative is block injective. Block \(L+1\) sites. Injectivity propagates from length \(L\) to length \(L+1\), so the blocked representatives are injective. Act on the first letter of each physical block and apply the blocked representative separation followed by the representative-grouped Lemma L. This gives equality of the two induced insertions on every blocked representative. The remaining \(L\) letters of the first block form a spanning family, so Lemma 16.4.19 recovers \(YA_j=ZA_j\) before blocking.
Let \((C_j^s)\) and \((E_j^s)\) be two marked families over a BNT canonical form. If their closed marked chains agree for every marked index \(s\) and every unmarked word \(w\) of positive length, then \(C_j^s=E_j^s\) for every representative \(j\) and every marked index \(s\).
The BNT canonical-form hypotheses give eventual representative word-tuple separation. Apply Theorem 16.4.14.
Let \(M\) be an MPO tensor in normalized BNT-refined horizontal form. For two coefficient families \(f_{u,z}\) and \(g_{u,z}\), put \(C^u=\sum _z f_{u,z}M^z\) and \(E^u=\sum _z g_{u,z}M^z\). If \(\operatorname{tr}(C^uM^w)=\operatorname{tr}(E^uM^w)\) for every marked index \(u\) and every word \(w\) of positive length, then \(C^u=E^u\) for every \(u\).
Restricting the marks to the span of the physical letters isolates the physical-insertion consequence of [ 13 , Appendix C.3, lines 1835–1858 ] needed for the contractions in [ 13 , Proposition 4.13, lines 1909–1919 ] .
Pass through the block-diagonal similarity in the normalized BNT-refined horizontal form. Linear combinations of physical letters are covariant under this similarity, and cyclicity of the trace removes it from each closed marked chain. On the representative-indexed sector decomposition one has
Corollary 16.4.21 identifies the two marks on every minimal representative. Taking their weighted direct sums and restoring the common block-diagonal similarity gives \(C^u=E^u\).
Let \(M\) be in normalized BNT-refined horizontal form. For doubled-index matrices \(Y,Z\), suppose that \(Y_1V^{(N)}(M)=Z_1V^{(N)}(M)\) for every \(N\geq 1\). Then the corresponding inserted tensors agree on the original bond space, \(YM=ZM\).
Now suppose in addition that \(M\) generates an MPDO, \(Q=Q^\dagger \), and that \(Q\widetilde M=Q\widetilde M Q\) on the vertically viewed tensor. Then, on every representative, the insertions corresponding to \(\widetilde M Q\) and \(Q\widetilde M Q\) agree: \(Y_R(Q)A_j=Y_C(Q)A_j\) for \(0\leq j{\lt}g\). This equality extends through every repeated weighted sector and the direct sum of the copy gauges, so that the original MPO tensor satisfies \(\widetilde M Q=Q\widetilde M Q\). The specialized assertion follows from the two entrywise identities \(Q_1H^{(N)}=Q_1H^{(N)}Q_1\) and \(Q_1H^{(N)}=H^{(N)}Q_1\) at every positive length.
The corresponding literal-CF assertion also holds. Let \(A\) be in literal CPSV canonical form. If the first-site actions of \(Y,Z\in M_{d}(\mathbb {C})\) agree on every positive-length matrix product vector of \(A\), then
on the original bond space. In particular, suppose that \(M\) generates matrix product density operators, \(P=P^\dagger \), and \(P\widetilde M=P\widetilde MP\). Positivity gives agreement of the first-site actions associated with \(\widetilde MP\) and \(P\widetilde MP\) at every positive length. Hence (33) gives the opposite one-sided identity on the original bond space.
This is the invariant-projection step of [ 13 , Proposition 4.13, lines 1873–1887 ] , grouping copies over each normal representative as in Lemma L [ 13 , Appendix C.3 ] ; Hermiticity of \(Q\) suffices for the implication.
For the normalized BNT-refined assertion, positive-length MPV equality gives the first-site identity for the BNT sector decomposition. The representative-grouped Lemma L gives \(YA_j=ZA_j\) for every representative. First-site insertion commutes with the weighted direct sum, and conjugating each copy by its gauge gives \(YM=ZM\) on the original bond space.
For the specialized assertion, positivity and Hermiticity give \(Q_1H^{(N)}=Q_1H^{(N)}Q_1=H^{(N)}Q_1\). These two identities give the corresponding first-site equalities for the BNT sector decomposition. Apply the representative-grouped Lemma L to identify first \(Y_L(Q)A_j\) with \(Y_C(Q)A_j\) and then with \(Y_R(Q)A_j\). Transitivity gives the stated equality for every representative, and the same direct-sum and copy-gauge argument gives the literal equality on \(\widetilde M\).
For literal CPSV canonical form, first-site agreement is a marked-trace identity on the original tensor. Theorem 16.4.17 needs only the identities with a nonempty unmarked tail; it identifies the two linear marks and gives (33). In the invariant-projection case, positivity first yields \(P_1H^{(N)}=P_1H^{(N)}P_1=H^{(N)}P_1\); its matrix coefficients are exactly the required first-site agreement.
Let \(\bigoplus _k \mu _k A_k\) be a block-diagonal tensor whose blocks are injective, left-canonical, weighted by nonzero scalars, and jointly block-injective in the sense that for some blocking length \(L\) the tuple of trace pairings against length-\(L\) block words is faithful across all blocks simultaneously. For \(Y\in M_{d}(\mathbb {C})\), write
for the physical-index insertion of \(Y\) into block \(A_k\), and for a configuration \(\sigma =(i_1,\ldots ,i_N)\) write
for \(Y\) acting on the first physical index of the MPV. If \(Y,Z\in M_{d}(\mathbb {C})\) satisfy \(Y_1V^{(N)}(A)=Z_1V^{(N)}(A)\) for every \(N\geq 1\), then \(YA_k=ZA_k\) for every block \(k\).
Scope restriction (block-injective specialization): The block-injective hypothesis excludes gauge-equivalent repeated blocks, specializing the source lemma [ 13 , Lemma L ] ; documented in [ 108 ] .
\(YH^{(N)}(\widetilde M)=ZH^{(N)}(\widetilde M)\)
\( = \)
\(YA_k=ZA_k\)
\( = \)
This is the per-block-separated form of Lemma L of [ 13 , Appendix C.3, lines 1835–1846 ] , using the finite criterion from Theorem 16.4.8 applied directly to the full block family rather than to a minimal basis of representatives. Expanding the hypothesis over a word of length \(L+1\) beginning with the fixed site and folding each block contribution into a trace pairing against the corresponding first-site insertion tensor yields, upon taking the \(Y\)-\(Z\) difference, a tuple \((\Delta _k)_k\) of block matrices. For every word \(w\) of length \(L\),
The finite block-separation hypothesis applied to the tuple \((\mu _k^{L+1}\Delta _k)_k\) forces each weighted block difference to be zero, and the nonzero weights give the blockwise equality of the two first-site insertions.
The trace-pairing step used above is the following local identity, \(\operatorname{tr}(WYA_k)=\operatorname{tr}(WZA_k)\):
For every \(W\) in the span of the length-\(N\) word evaluations of the canonical form, it states \(\operatorname{tr}(WA^Y)=\operatorname{tr}(WA^Z)\), where \(A^Y\) is obtained by applying \(Y\) to the physical index. Compare the trace-pairing identities in the proof of Lemma L, [ 13 , Appendix C.3, lines 1848–1858 ] .
Let \(H\) and \(P\) be Hermitian matrices. If \(PH=PHP\), then \((\mathbb {1}-P)HP=0\).
Taking adjoints gives \(HP=PHP\), and the conclusion follows by subtraction.
Let \(H\) be a nonempty-chain density operator generated by an MPDO, and let \(P=P^\dagger \). If \(PH=PHP\), then \((\mathbb {1}-P)HP=0\). This is the operator identity used in the proof of [ 13 , Proposition 4.13 ] .
Since \(H\) and \(P\) are Hermitian, taking adjoints gives \(HP=(PH)^\dagger =(PHP)^\dagger =PHP\). Therefore \((\mathbb {1}-P)HP=HP-PHP=0\).
Let \(\bigoplus _{j,q}\mu _{j,q}A_j\) be a BNT canonical form and let \(Y_L,Y_C,Y_R\in M_{d}(\mathbb {C})\), with \(Y_1V^{(N)}(A)\) denoting \(Y\) acting on the first physical index of the MPV. If, for every \(N\geq 1\),
then \(Y_RA_j=Y_CA_j\) for every minimal representative \(A_j\). This application of Lemma L specializes to the physical-index insertions of \(P\), \(PHP\), and \(HP\) in Theorem 16.4.31.
Apply the representative-grouped Lemma L first to the actions \(Y_L,Y_R\) and then to \(Y_L,Y_C\). The BNT canonical-form hypotheses supply the required common block-injectivity length. Transitivity gives the stated equality on every minimal representative; the repeated copies have already been combined through \(\sum _q\mu _{j,q}^{N}\).
Let \(P\in M_{d}(\mathbb {C})\) and write doubled indices as \(p=(a,b)\) and \(q=(c,e)\). The first-site matrices corresponding respectively to \(PH^{(N)}\), \(H^{(N)}P\), and \(PH^{(N)}P\) are
The horizontal operator \(H^{(N)}(\widetilde M)\) is the cyclic word \(H^{(N)}(\widetilde M)=\operatorname{tr}_{\mathrm v}(\widetilde M_1\cdots \widetilde M_N)\):
Its virtual indices are contracted around the row, while every ket and bra index remains external. Reading each ket–bra pair as one doubled index identifies its entries with the matrix product vector of the doubled-index tensor. This is the contraction in [ 13 , Proposition 4.13, lines 1876–1880 ] .
For ket and bra configurations \(\sigma ,\tau \), \(V^{(N)}(\widehat M)_{(\sigma ,\tau )} =H^{(N)}(M)_{\sigma ,\tau }\). If \(p_0=(a,b)\) is the first doubled index and \(\sigma ',\tau '\) are the ket and bra configurations on the remaining \(N\) sites, contraction with \(Y_L\), \(Y_R\), and \(Y_C\) gives, respectively,
Decode the doubled index at every site as a ket–bra pair. The defining trace of the resulting matrix word is exactly the corresponding entry of \(H^{(N)}(M)\). Summing the first doubled index against \(Y_L\), \(Y_R\), or \(Y_C\) leaves the single or double sums in (34).
Let \(\widehat M\) be the doubled-index tensor associated with \(M\). Then insertion of \(Y_L\), \(Y_R\), and \(Y_C\) gives, respectively, the doubled-index tensors associated with \(P\widetilde M\), \(\widetilde MP\), and \(P\widetilde MP\).
Write the output and summed doubled indices as \((a,b)\) and \((c,e)\). The Kronecker deltas in \(Y_L\) and \(Y_R\) remove, respectively, the sums over \(e\) and \(c\); \(Y_C\) leaves both multiplications. The resulting coefficients are \(\sum _c P_{a,c}M^{cb}\), \(\sum _e P_{e,b}M^{ae}\), and \(\sum _{c,e}P_{a,c}P_{e,b}M^{ce}\).
Let the doubled-index tensor \(\widehat M\) generate the same positive-length MPV family as a BNT canonical form \(\bigoplus _{j,q}\mu _{j,q}A_j\). If, entrywise at every positive length,
then for every minimal representative \(A_j\) the insertions associated with \(H^{(N)}P\) and \(PH^{(N)}P\) agree. Repeated copies over \(A_j\) are grouped through the coefficient \(\sum _q\mu _{j,q}^{N}\).
View the MPO tensor \(M\) vertically as the family of physical-space operators \(M_{ab}\) of Definition 16.4.4, indexed by the virtual indices, with matrix elements \((M_{ab})_{ij}=M^{ij}_{ab}\). A matrix \(P\in M_{d}(\mathbb {C})\) then acts on the tensor by
and on an \((N+1)\)-site chain by the first-spin operator \(P\otimes \mathbb {1}^{\otimes N}\). On the letters of the vertically viewed tensor these actions are the left and right matrix products, \((PM)_{ab}=PM_{ab}\) and \((MP)_{ab}=M_{ab}P\). The first-spin operators compose site by site, \(P_1Q_1=(PQ)_1\).
Let \(M\) generate an MPDO and let \(P=P^\dagger \) satisfy \(PM=PMP\) on the vertically viewed tensor. For every total length \(L\geq 1\), write \(P_1=P\otimes \mathbb {1}^{\otimes (L-1)}\). Then \(P_1H^{(L)}=P_1H^{(L)}P_1=H^{(L)}P_1\). This is the equality chain in the proof of [ 13 , Proposition 4.13 ] ; the proposition assumes \(P\) is an orthogonal projector, but only self-adjointness of \(P\) is used.
Put \(L=N+1\). Splitting off the first site of a density-operator entry gives
Summing the first ket index against \(P\) therefore pairs the operator \((PM)^{ij}\) against the remaining word; the hypothesis \(PM=PMP\) replaces it by \((PMP)^{ij}\), which unfolds to the entry of \(P_1H^{(N+1)}P_1\). This proves \(P_1H^{(N+1)}=P_1H^{(N+1)}P_1\). Since \(H^{(N+1)}\) is positive semidefinite and \(P\) is self-adjoint, Theorem 16.4.26 gives \((\mathbb {1}-P_1)H^{(N+1)}P_1=0\), hence \(H^{(N+1)}P_1=P_1H^{(N+1)}P_1=P_1H^{(N+1)}\).
Let \(M\) generate an MPDO whose doubled-index tensor \(\widehat M\) generates the same positive-length MPV family as a BNT canonical form \(\bigoplus _{j,q}\mu _{j,q}A_j\), and let \(P=P^\dagger \) satisfy \(PM=PMP\) on the vertically viewed tensor. Then on every minimal representative \(A_j\) the insertions of \(MP\) and \(PMP\) agree; equivalently, \((\mathbb {1}-P)MP=0\) on every representative.
For every tail length \(N\), Theorem 16.4.33 turns \(PM=PMP\) into the entrywise identities \(P_1H^{(N+1)}=P_1H^{(N+1)}P_1\) and \(P_1H^{(N+1)}=H^{(N+1)}P_1\). These are exactly the hypotheses of Theorem 16.4.31, whose conclusion is the equality of the insertions of \(MP\) and \(PMP\) on every minimal representative.
Let \(H\geq 0\), let \(p\geq 1\), and let \(QH^p=H^pQ\). Then \(QH=HQ\).
Diagonalize \(H\). The function \(x\mapsto x^p\) is injective on the non-negative reals, so a matrix commutes with the diagonal matrix of the \(p\)-th powers exactly when it commutes with the diagonal matrix itself.
A matrix commuting with a nonzero power of a nonempty-chain density operator generated by an MPDO commutes with the density operator itself. This is the final operator implication in the periodic-sector contradiction in the proof of [ 13 , Proposition 4.13 ] .
Diagonalize \(H=U\Lambda U^\dagger \) with \(U\) unitary and \(\Lambda \) the diagonal matrix of eigenvalues \(\lambda _i\geq 0\). Conjugating the hypothesis by \(U\) gives, entrywise for \(Q'=U^\dagger QU\),
Since \(x\mapsto x^p\) is injective on the non-negative reals, every entry \(Q'_{ij}\) with \(\lambda _i\neq \lambda _j\) vanishes, so \(Q'\Lambda =\Lambda Q'\). Conjugating back gives \(QH=HQ\).
Let \(M\) be an MPO tensor. A periodic-sector projector for \(M\) with exponent \(p\ge 1\) is an orthogonal projector \(Q\) on the physical space such that for every \(N\ge 0\), with \(Q_1=Q\otimes \mathbb {1}^{\otimes N}\) acting on the first site,
These are the commutation identities that the proof of [ 13 , Proposition 4.13 ] , lines 1889–1891, attaches to a nontrivial vertical period; the second family is the displayed identity at lines 1890–1891.
Let \(M\) generate matrix product density operators. Then \(M\) admits no periodic-sector projector.
Suppose \(Q\) were such a projector. Every \(H^{(N+1)}\) is positive semidefinite, so commutation of \(Q_1\) with \([H^{(N+1)}]^p\) implies commutation with \(H^{(N+1)}\) (Theorem 16.4.36), contradicting \(Q_1H^{(N+1)}\ne H^{(N+1)}Q_1\). This is the contradiction concluding the periodic-sector step in the proof of [ 13 , Proposition 4.13 ] , lines 1892–1893.
Say that periodic vectors of \(M\) supply projectors if for every nontrivial \(p\)-periodic vector of the vertically viewed tensor \(\widetilde M\) — an irreducible restriction \(B\) of \(\widetilde M\) to a minimal invariant subspace whose transfer map has a peripheral eigenvalue different from its spectral radius, as in Definition 7.7.4 — there is an exponent \(q\ge 1\) for which \(M\) admits a periodic-sector projector. The proof of [ 13 , Proposition 4.13 ] , lines 1889–1891, asserts this implication from the general theory of the peripheral spectrum of completely positive maps, with \(q\) equal to the period \(p\) of the vector; only the existence of some exponent enters the exclusion argument. The projector, its word invariance, and its single-letter displacement are constructed from the cyclic decomposition of the peripheral spectrum (Theorem 16.4.47); the single non-commutation instance suffices for the periodic-sector contradiction below.
Let \(M\) generate matrix product density operators and suppose that periodic vectors of \(M\) supply projectors. Then the vertically viewed tensor \(\widetilde M\) has no nontrivial \(p\)-periodic vectors. This is the periodic-sector step in the proof of [ 13 , Proposition 4.13 ] , lines 1888–1893, granted the projector that the source draws from the general theory of the peripheral spectrum.
A nontrivial \(p\)-periodic vector of \(\widetilde M\) supplies a periodic-sector projector for \(M\), and Theorem 16.4.38 shows that no such projector exists.
Placing one MPO layer on top of another and contracting the shared physical index gives the layer product of two MPO tensors \(M\) and \(M'\) with the same physical dimension, \((MM')^{ij}=\sum _kM^{ik}\otimes M'^{kj}\), an MPO tensor on the tensor product of the virtual spaces. The \(p\)-fold vertical stack \(M^{[p]}\) of \(M\) is the iterated layer product of \(p\) copies of \(M\); the empty stack \(M^{[0]}\) is the trivial-bond tensor with letters \(\delta _{ij}\), which generates the identity operator at every system size.
For every system size \(N\), the operator generated by a layer product is the product of the density operators generated by the layers, \(H^{(N)}(MM')=H^{(N)}(M)\, H^{(N)}(M')\). Consequently the \(p\)-fold vertical stack generates \(H^{(N)}(M^{[p]})=[H^{(N)}(M)]^p\) for all \(N\) and \(p\); this is the reading of the powers in the displayed identity of [ 13 , Proposition 4.13, lines 1890–1891 ] as density operators of stacked layers. In particular, the vertical stack of a tensor generating matrix product density operators again generates matrix product density operators on every nonempty chain.
The word product of layer-product letters along a ket–bra configuration expands, by induction on the length and the mixed-product property of the Kronecker product, as a sum over configurations \(\kappa \) of the contracted physical index of Kronecker products of the layers’ word products. Taking traces, each summand factors into the product of the two traces, so the entry of \(H^{(N)}(MM')\) at \((\sigma ,\tau )\) equals \(\sum _\kappa H^{(N)}(M)_{\sigma \kappa }H^{(N)}(M')_{\kappa \tau }\), the entry of the matrix product. The power identity follows by induction on \(p\), and positive semidefiniteness of \([H^{(N)}]^p\) follows from that of \(H^{(N)}\).
The vertically viewed letters of a layer product are the products of the layers’ vertically viewed letters, \((MM')_{(a,a'),(b,b')}=M_{ab}M'_{a'b'}\). Hence every vertically viewed letter of the \(p\)-fold stack \(M^{[p]}\) is the evaluation of a length-\(p\) word of \(\widetilde M\). In particular, if a matrix \(Q\) satisfies the one-sided invariance \(QA_w=QA_wQ\) for every product \(A_w\) of \(p\) letters of \(\widetilde M\), then \(QM^{[p]}=QM^{[p]}Q\) in the sense of Definition 16.4.32.
The letter identity is a direct computation from the layer product: contracting the shared physical index composes the physical-space operators. Induction on \(p\) writes each vertical letter of \(M^{[p]}\) as a length-\(p\) word. Thus, for every bond index \(v\), there is a word \(w\) of length \(p\) such that \(\widetilde{M^{[p]}}_v=A_w\), and \(QA_w=QA_wQ\) implies \(Q\widetilde{M^{[p]}}_v=Q\widetilde{M^{[p]}}_vQ\). Applying this identity to the two one-site actions gives \(QM^{[p]}=QM^{[p]}Q\).
Let \(M\) generate matrix product density operators, let \(p\ge 1\), and let \(Q\) be an orthogonal projector on the physical space satisfying \(QA_w=QA_wQ\) for every product \(A_w\) of \(p\) letters of \(\widetilde M\) and \(Q_1H^{(N)}\ne H^{(N)}Q_1\) for every length. Then \(Q\) is a periodic-sector projector for \(M\) with exponent \(p\).
By Theorem 16.4.43 the word invariance gives \(QM^{[p]}=QM^{[p]}Q\) on the stacked tensor, which generates matrix product density operators by Theorem 16.4.42. Theorem 16.4.33 applied to \(M^{[p]}\) yields \(Q_1H^{(N)}(M^{[p]})=H^{(N)}(M^{[p]})Q_1\) for every length, and \(H^{(N)}(M^{[p]})=[H^{(N)}]^p\) by the stacked-layers identity. Together with the assumed non-commutation family this is exactly the periodic-sector projector of Definition 16.4.37.
Say that periodic vectors of \(M\) supply cyclic projectors if for every nontrivial \(p\)-periodic vector of the vertically viewed tensor \(\widetilde M\) — presented as in Definition 7.7.4 — there are an exponent \(q\ge 1\) and an orthogonal projector \(Q\) on the physical space satisfying the one-sided invariance \(QA_w=QA_wQ\) for every product \(A_w\) of \(q\) letters of \(\widetilde M\), together with \(Q_1H^{(N)}\ne H^{(N)}Q_1\) for every length. This is the shape of the cyclic decomposition of the peripheral spectrum [ 244 , Theorem 6.6 ] applied to the vertical transfer map: each cyclic projection of a period-\(p\) sector returns to itself after \(p\) applications of the map, while its nontrivial cyclic displacement under one application prevents commutation with \(H^{(N)}\). Compared with Definition 16.4.39, the commutation of \(Q\) with \([H^{(N)}]^p\) is no longer assumed; it is derived through the stacked-layers identity.
Let \(M\) generate matrix product density operators. If periodic vectors of \(M\) supply cyclic projectors, then they supply periodic-sector projectors. This reduces the projector-supply hypothesis of the periodic-sector step in the proof of [ 13 , Proposition 4.13, lines 1888–1893 ] to the cyclic decomposition of the peripheral spectrum of the vertical transfer map.
A nontrivial \(p\)-periodic vector supplies an invariant non-commuting projector, which is a periodic-sector projector by Definition 16.4.44.
Let \(M\) be an MPO tensor and let \(\widetilde M\) carry a nontrivial \(p\)-periodic vector, presented as in Definition 7.7.4: an irreducible restriction \(B\) of \(\widetilde M\) to a minimal invariant subspace, reached through an isometry \(V\), with a positive definite transfer eigenvector \(\rho \) of eigenvalue \(r{\gt}0\) and a transfer eigenvalue \(\mu \ne r\) of modulus \(r\). Then there are an exponent \(q\ge 1\) and an orthogonal projector \(Q\) on the physical space such that \(QA_w=QA_wQ\) for every product \(A_w\) of \(q\) letters of \(\widetilde M\), while \(Q\widetilde M\ne Q\widetilde MQ\) in the sense of Definition 16.4.32. This realizes the projector that the proof of [ 13 , Proposition 4.13, lines 1889–1891 ] draws from the general theory of the peripheral spectrum.
The rescaled unital gauge \(K\) of the corner satisfies the hypotheses of the cyclic decomposition of the peripheral spectrum (Theorem D.3.5), and the eigenvalue \(\mu /r\ne 1\) of modulus one (Lemma D.3.6) forces the peripheral cyclic group to have order \(m\ge 2\). Theorem [ 216 , “Cyclic decomposition for irreducible finite Kraus maps” ] produces orthogonal projections \(P_0,\ldots ,P_{m-1}\) summing to the identity with \(\mathcal{E}_K(P_{k+1})=P_k\), and the letter-level shift \(K_vP_{k+1}=P_kK_v\) holds (Lemma 5.8.1.1). Set \(q=m\) and let \(Q\) be the complement of the support projection \(\pi \) of \(YY^\dagger \) for \(Y=V\rho ^{1/2}P_0\). A word \(A_w\) of \(m\) letters of \(\widetilde M\) satisfies \(A_wY=YG_w\), where \(G_w\) is a scalar multiple of the corresponding word in the letters of \(K\), which commutes with \(P_0\) by the full-period shift; so the range of \(Y\) is invariant and Lemma D.3.3 gives \(QA_w=QA_wQ\). If a single letter also satisfied the invariance, then \(\pi \) would fix \(A_vY\), whose columns lie in the sector \(V\rho ^{1/2}P_{m-1}\) by the cyclic shift; cancelling the injective factor \(V\rho ^{1/2}\) and using orthogonality of \(P_{m-1}\) and \(P_0\) would force \(P_{m-1}K_v=0\) for every letter, hence \(P_{m-1}=\mathcal{E}_K(P_0)=0\), contradicting the nonvanishing of the cyclic projections.
Let \(M\) be an MPO tensor whose doubled-index MPS tensor is in literal CPSV canonical form. If \(Q^2=Q\) and \(Q\widetilde M\ne Q\widetilde M Q\), then there is an integer \(N\geq 0\) such that, for \(Q_1=Q\otimes \mathbb {1}^{\otimes N}\),
This existential non-commutation statement supplies the single length needed for the contradiction in [ 13 , Proposition 4.13, lines 1888–1893 ] .
Suppose instead that \(Q_1H^{(N+1)}=H^{(N+1)}Q_1\) for every \(N\geq 0\). Idempotence gives
Taking matrix coefficients yields equality of the one-sided first-site contractions. The literal form of Lemma L in Theorem 16.4.23 identifies the two inserted tensors on the original bond space. By Lemma 16.4.30, this equality is \(Q\widetilde M=Q\widetilde M Q\), contradicting the displacement.
Let \(M\) generate matrix product density operators, and suppose that its doubled-index MPS tensor is in literal CPSV canonical form. Then the vertically viewed tensor \(\widetilde M\) has no nontrivial periodic vectors.
A nontrivial periodic vector supplies an integer \(p\geq 1\) and an orthogonal projector \(Q\) such that every word \(A_w\) of length \(p\) in the letters of \(\widetilde M\) satisfies \(QA_w=QA_wQ\), while \(Q\widetilde M\ne Q\widetilde M Q\) (Theorem 16.4.47). By Theorem 16.4.48, there is a length \(N+1\) for which
The word invariance and the vertical-stack identity give
Since \(H^{(N+1)}\) is positive semidefinite, Theorem 16.4.36 implies \(Q_1H^{(N+1)}=H^{(N+1)}Q_1\), a contradiction. This is the periodic exclusion in [ 13 , Proposition 4.13, lines 1888–1893 ] , with the single noncommuting length supplied by Lemma L.
Let \(M\) be in normalized BNT-refined horizontal form. If an idempotent \(Q\) satisfies \(Q\widetilde M\ne Q\widetilde MQ\), then for some \(N\ge 0\), \(Q_1H^{(N+1)}\ne H^{(N+1)}Q_1\). Consequently, if \(M\) generates matrix product density operators, then \(\widetilde M\) has no nontrivial periodic vectors.
Scope restriction (BNT-refined horizontal form): The normalized BNT-refined horizontal hypothesis is stronger than the literal CPSV canonical-form hypothesis; documented in [ 177 ] .
Suppose first that \(Q_1\) commutes with \(H^{(N+1)}\) for every \(N\). Since \(Q^2=Q\), this gives \(Q_1H^{(N+1)}=Q_1H^{(N+1)}Q_1\) at every positive length. Theorem 16.4.23 and Lemma 16.4.30 then imply \(Q\widetilde M=Q\widetilde MQ\), a contradiction. Thus a displaced idempotent has a noncommuting length \(N+1\).
A nontrivial periodic vector supplies the orthogonal projector \(Q\) of Theorem 16.4.47, invariant under every word of some full period \(p\ge 1\). At the noncommuting length just obtained, the stacked-layer identity and the invariant-projection calculation give \(Q_1[H^{(N+1)}]^p=[H^{(N+1)}]^pQ_1\). Positive semidefiniteness and Theorem 16.4.36 imply that \(Q_1\) commutes with \(H^{(N+1)}\), which yields the contradiction in [ 13 , Proposition 4.13, lines 1888–1893 ] using a single noncommuting length.
Let \(M\) be an MPO tensor whose doubled-index tensor is injective, and let \(Q\) be an idempotent on the physical space (Hermiticity is not needed). If \(Q_1H^{(2)}=H^{(2)}Q_1\), then \(Q\widetilde M=Q\widetilde MQ\).
Idempotence of \(Q\) upgrades the commutation to the one-sided reduction \(Q_1H^{(2)}=Q_1H^{(2)}Q_1\). Unfolding this identity against every trailing letter of \(M\) gives, for every pair of physical indices \(a,b\), a trace identity pairing \((Q\widetilde M)_{ab}\) against \((Q\widetilde MQ)_{ab}\) over every letter of the doubled-index tensor. Injectivity’s spanning property and nondegeneracy of the trace pairing force the two matrices to agree.
Let \(M\) generate matrix product density operators and let \(M\)’s doubled-index tensor be injective. Then the vertically viewed tensor \(\widetilde M\) has no nontrivial periodic vectors, unconditionally.
A nontrivial periodic vector supplies the invariant, displaced projector \(Q\) of Theorem 16.4.47. Its word invariance transfers to the stacked tensor, giving commutation with \([H^{(2)}]^p\) unconditionally through the invariant-projection step applied to the stack; positive semidefiniteness of the density operators removes the power, and Theorem 16.4.51 turns the resulting commutation with \(H^{(2)}\) into the letter-level invariance that the displacement forbids.
Theorem 16.4.50 gives the corresponding result for normalized BNT-refined horizontal form.
A nonzero positive semidefinite matrix has positive trace: the trace is a positive real.
The trace of a positive semidefinite matrix is a non-negative real, and it vanishes only if the matrix itself is zero.
Let \(H^{(N)}\), \(N\geq 1\), be a density operator generated by an MPDO and let \(P=P^\dagger \). Then \(PH^{(N)}P\ge 0\), and if \(PH^{(N)}P\ne 0\) its trace is a positive real. In particular this holds for the first-site compressions \(P_1H^{(N+1)}P_1\). Applied to the sector projectors \(P_{\alpha ,k}\) of the vertical decomposition, this is the positivity \(\operatorname{tr}(P_{\alpha ,k}H^{(N)}P_{\alpha ,k}){\gt}0\) from which the proof of [ 13 , Proposition 4.13 ] concludes \(\mu _{\alpha ,k}d_\alpha {\gt}0\); the identification of the trace with \(\mu _{\alpha ,k}d_\alpha \) is Theorem 16.4.62.
Compression by a Hermitian matrix preserves positive semidefiniteness, since \(PH^{(N)}P=PH^{(N)}P^\dagger \). By Theorem 16.4.53, a nonzero compression therefore has positive trace. The first-site case follows because \(P_1=P\otimes \mathbb {1}^{\otimes N}\) is Hermitian whenever \(P\) is.
Closing the vertical indices of a single site of an MPO tensor \(M\) gives the matrix \(E=\sum _i M^{ii}\in M_{D}(\mathbb {C})\), and inserting a matrix \(P\in M_{d}(\mathbb {C})\) into the loop gives \(E_P=\sum _{i,j}P_{ji}M^{ij}\), the closed loop of the vertically viewed product \(PM\), equivalently of \(MP\) by cyclicity of the single-site trace; inserting the identity recovers \(E\). On a chain of length \(N+1\), the compression of the density operator by \(P\) at the first site is \(P_1H^{(N+1)}P_1\) with \(P_1=P\otimes \mathbb {1}^{\otimes N}\), and for a matrix \(E_\alpha \in M_{D}(\mathbb {C})\) the chain trace with one marked loop is \(d_\alpha ^{(N+1)}=\operatorname{tr}(E_\alpha E^N)\). For a representative tensor \(M_\alpha \) whose letters are indexed by horizontal bond pairs, its closed vertical loop is defined entrywise by \((E_\alpha )_{ab}=\operatorname{tr}(M_\alpha ^{(a,b)})\). In this case the chain trace is the scalar \(d_\alpha \) displayed in the proof of [ 13 , Proposition 4.13 ] , resolved by chain length: one loop of \(M_\alpha \) followed by loops of \(M\), closed by the horizontal trace.
Summing the diagonal word evaluations of an MPO tensor over all length-\(N\) configurations closes the vertical loops site by site,
and hence \(\operatorname{tr}H^{(N)}=\operatorname{tr}(E^N)\).
By induction on \(N\): splitting off the first site factors the sum as \((\sum _a M^{aa})\) times the sum over the remaining \(N\) sites. Taking horizontal traces and summing the diagonal entries of \(H^{(N)}\) gives the trace identity.
For every matrix \(P\in M_{d}(\mathbb {C})\) and every chain length \(N+1\),
The marked site closes to the loop with \(P\) inserted, and the remaining \(N\) sites close to plain loops. If \(P\) is idempotent, then also \(\operatorname{tr}(P_1H^{(N+1)}P_1)=\operatorname{tr}(E_P E^N)\).
For a length-\(N\) configuration \(\rho \), write \(E_\rho ^{(N)}=M^{\rho _0\rho _0}\cdots M^{\rho _{N-1}\rho _{N-1}}\). Expanding at the first site gives
Summing over \(a\) and \(\rho \) and collecting closes the remaining loops by Lemma 16.4.56. For the two-sided compression, first-site actions compose site by site, \(P_1Q_1=(PQ)_1\), so cyclicity of the trace and idempotency reduce \(\operatorname{tr}(P_1H^{(N+1)}P_1)\) to \(\operatorname{tr}(P_1H^{(N+1)})\).
A matrix \(P\in M_{d}(\mathbb {C})\) is a sector projector for the weight \(\mu \in \mathbb {C}\) and the matrix \(E_\alpha \in M_{D}(\mathbb {C})\) if \(P\) is an orthogonal projector and inserting \(P\) into the single-site vertical loop selects the sector contribution, \(E_P=\mu E_\alpha \). In the vertical decomposition
of the vertically viewed tensor, carried by the physical space in the proof of [ 13 , Proposition 4.13 ] , the projector \(P_{\alpha ,k}\) selecting the sector \((\alpha ,k)\) has this property with \(\mu =\mu _{\alpha ,k}\) and \((E_\alpha )_{ab}=\operatorname{tr}(M_\alpha ^{(a,b)})\), independently of \(k\). Indeed, cyclicity of trace cancels the internal similarity \(X_{\alpha ,k}M_\alpha ^{(a,b)}X_{\alpha ,k}^{-1}\). Applied to the physical isometry of a grouped normal sector, the gauge-corner identity constructs the projector and proves its loop identity.
Let \(M\) be an MPO tensor whose doubled-index MPS tensor is in literal CPSV canonical form, and let \(P\in M_{d}(\mathbb {C})\). If one vertical letter has a nonzero two-sided corner, then some finite-chain first-site compression is nonzero:
where \(P_1=P\otimes \mathbb {1}^{\otimes N}\). This is the separation assertion at [ 13 , Proposition 4.13, line 1898 ] , obtained from Lemma L in [ 13 , Appendix C.3, lines 1835–1858 ] . It uses neither positivity of \(H^{(N+1)}\) nor a grouping of the vertical sectors.
Suppose that \(P_1H^{(N+1)}P_1=0\) for every \(N\geq 0\). By the coefficient formula for the doubled-index first-site contraction, the two-sided action \(Y_C(P)\) and the zero action agree on every positive-length matrix product vector of the doubled-index tensor. Literal Lemma L therefore identifies their inserted tensors. For every vertical letter \(v\), the insertion identity then gives
This contradicts the hypothesis in (39).
Let \(M\) be an MPO tensor of physical dimension \(d\) and bond dimension \(D\). Let \(A\) be a CPSV normal tensor of positive bond dimension \(n\) and physical dimension \(D^2\), let \(V\in \mathbb {C}^{d\times n}\) be an isometry, let \(X\in \mathbb {C}^{n\times n}\) be invertible, and let \(\omega \) be nonzero. Suppose that, for every vertical letter \(v\),
Then the range projector \(VV^\dagger \) has a nonzero vertical corner:
for some \(v\). This is the algebraic nonvanishing step used in [ 13 , Proposition 4.13, line 1898 ] .
Let \(M\) be in normalized BNT-refined horizontal form and let \(P\) be a matrix for which \(P\widetilde M_vP\) is nonzero for some vertical letter \(v\). Then the first-site compression \(P_1H^{(N+1)}P_1\) is nonzero for some \(N\). This is the separation step used in [ 13 , Proposition 4.13, line 1898 ] .
The normalized BNT-refined horizontal hypothesis provides the grouped vertical sectors used here; the corresponding literal separation statement is Theorem 16.4.59.
If every compression vanished, its matrix entries would say that the first-site insertion by \(P\) has the same matrix product vectors as zero. Lemma L on the normalized BNT-refined horizontal form would then force the inserted tensor itself to vanish, contrary to the assumed nonzero corner.
Let \(P\) be a sector projector for the weight \(\mu \) and the matrix \(E_\alpha \). Then for every chain length \(N+1\),
Once spectral normalization and grouping identify the reducing projectors with sector projectors \(P_{\alpha ,k}\) satisfying Definition 16.4.58, this specializes to the trace identification \(\operatorname{tr}(P_{\alpha ,k}H^{(N)}P_{\alpha ,k})=\mu _{\alpha ,k}d_\alpha \) used in the proof of [ 13 , Proposition 4.13 ] ; the source states it at the particular length where the compression does not vanish.
By Theorem 16.4.57 the trace of the compression is \(\operatorname{tr}(E_P E^N)\), and the sector-projector property replaces \(E_P\) by \(\mu E_\alpha \).
Let \(M\) generate an MPDO, let \(Q\) and \(P\) be sector projectors for the weights \(\mu _{\alpha ,0}\) and \(\mu _{\alpha ,k}\) and a common matrix \(E_\alpha \), and suppose \(\mu _{\alpha ,0}{\gt}0\). If the compression of \(H^{(N+1)}\) by \(P\) at the first site is nonzero, then
Granted such sector projectors, this yields the conclusion of the trace argument in the proof of [ 13 , Proposition 4.13 ] : since w.l.o.g. \(\mu _{\alpha ,0}{\gt}0\), the positive trace forces \(d_\alpha {\gt}0\) and hence \(\mu _{\alpha ,k}{\gt}0\).
By Theorem 16.4.54 the nonzero compression has positive trace, which Theorem 16.4.62 identifies with \(\mu _{\alpha ,k}d_\alpha ^{(N+1)}\). In particular \(d_\alpha ^{(N+1)}\ne 0\), so the trace of the compression by the sector projector of \((\alpha ,0)\) at the same length equals \(\mu _{\alpha ,0}d_\alpha ^{(N+1)}\ne 0\); that compression is therefore also nonzero, and its trace is a positive real. Dividing by \(\mu _{\alpha ,0}{\gt}0\) gives \(d_\alpha ^{(N+1)}{\gt}0\), and dividing the first inequality by \(d_\alpha ^{(N+1)}\) gives \(\mu _{\alpha ,k}{\gt}0\).
If \(\omega {\gt}0\) and \(X^\dagger X=\omega \mathbb {1}\), then \(\omega ^{-1/2}X\) is unitary.
Divide the Gram identity by \(\omega \).
Let \(M_\alpha \) be a normal tensor. Any matrix commuting with every matrix of \(M_\alpha \) is a scalar multiple of the identity. In particular, if \(X\ne 0\) and \(X^\dagger X\) commutes with every matrix of \(M_\alpha \), then \(X^\dagger X=\omega \mathbb {1}\) for a necessarily positive constant \(\omega \), and \(\omega ^{-1/2}X\) is unitary. This is the normalized form \(X_{\alpha ,k}^\dagger X_{\alpha ,k}=\omega _{\alpha ,k}\mathbb {1}\) of the Gram-matrix proportionality in the proof of [ 13 , Proposition 4.13 ] , with \(X_{\alpha ,0}=\mathbb {1}\) and \(U_{\alpha ,k}=\omega _{\alpha ,k}^{-1/2}X_{\alpha ,k}\). Deriving the commutation from self-adjointness of the sector compressions is a separate step.
Commutation with the matrices of \(M_\alpha \) extends to all word evaluations, and after blocking these span the full matrix algebra, so the commuting matrix is central and hence scalar. The Gram matrix \(X^\dagger X\) is positive semidefinite and nonzero, so its scalar value is a positive real \(\omega \); dividing \(X\) by \(\sqrt{\omega }\) makes the Gram matrix the identity.
Let \(M\) generate positive semidefinite operators on every nonempty chain, let \(A\) be a vertical tensor, and suppose that a matrix \(V\) selects the vertical corner \(V^\dagger \widetilde M_vV=cA^v\), where \(c{\gt}0\). For vertical indices \(r,s\), write \((\mathcal{R}(A)^{rs})_{ab}=(A^{(a,b)})_{rs}\), and set \((A^\sharp )^v=(A^{v^{\mathrm{op}}})^\dagger \). Then, for every pair of equal-length tail words \(\sigma ,\tau \),
Adjunction exchanges the oriented horizontal bond pair and reverses the unmarked tail, without identifying \(A\) with \(A^\sharp \) letterwise.
Compress the first site of the positive operator by \(V\). The resulting matrix is Hermitian. Expanding one matrix entry gives the marked chain on the left and the complex conjugate of the entry with \(r,s\) and the ket–bra words exchanged. Taking the conjugate transpose reverses the tail product, cyclicity of trace returns the marked factor to the first position, and positivity of \(c\) permits cancellation without a phase.
Let \(M\) generate positive semidefinite operators on every nonempty chain, and suppose that one vertical corner is a positive scalar multiple of a similarity of a tensor \(A\),
Then the marked-chain coefficient of the Gram-dressed tensor \((X^\dagger X)A^v(X^\dagger X)^{-1}\) equals the reflected-adjoint marked-chain coefficient of \(A\) in the adjoint tensor, with the two tail words reversed. This is the single-corner, open-index conjugation step from Figures 7–8 of [ 13 , Proposition 4.13, lines 1909–1919 ] .
Apply the reflected marked-chain identity to the similar tensor \(XA^vX^{-1}\). Conjugate its two open virtual indices by \(X^\dagger \) and \((X^{-1})^\dagger \). On the original side this gives the Gram dressing of \(A\); on the reflected side the gauge factors cancel, leaving the reflected adjoint of \(A\).
Let \(M\) be in normalized BNT-refined horizontal form and generate positive semidefinite operators on every nonempty chain. Suppose two vertical corners are positive scalar multiples of similarities of the same tensor \(A\):
Then the two Gram conjugations agree for every vertical letter,
Normality of \(A\) is not required for this equality; it enters only in the subsequent proportionality of the Gram matrices. This is the blockwise statement underlying the unnormalized pairwise conclusion of Figures 7–8 in [ 13 , Proposition 4.13, lines 1909–1919 ] .
Scope restriction (BNT-refined horizontal form): The normalized BNT-refined horizontal hypothesis is stronger than the literal CPSV canonical-form hypothesis; documented in [ 177 ] .
Conjugate the open indices in the reflected marked-chain identity for the first corner by \(X^\dagger \). Its left side becomes the marked chain of \((X^\dagger X)A^v(X^\dagger X)^{-1}\), while the gauge factors on the reflected side cancel and leave the marked chain of \(A^\sharp \). The same calculation with \(Y\) gives the identical reflected chain. Thus the two Gram-dressed marks have equal closed chains for every tail word. The reflected identity also holds for the empty tail, but marked separation uses only nonempty tails, as in the proof of Lemma L in [ 13 , Appendix C.3, lines 1848–1858 ] .
These marks lie in the physical-letter span. Indeed, with \(L_X=V_XX\) and \(R_X=V_X(X^{-1})^\dagger \), their coefficients are \(f^X_{(r,s),(i,j)} =c_X^{-1}\overline{(L_X)_{ir}}(R_X)_{js}\), and similarly for \(Y\). The corner identities show that these linear combinations are precisely the horizontal slices of the two Gram-dressed tensors. Physical-letter marked separation therefore identifies every slice, hence every vertical letter.
An MPO tensor \(M\) has the grouped-corner Gram-dressing property if the following holds. Let \(A\) be a tensor, let \(V_X,V_Y\) be two physical maps, let \(X,Y\) be invertible matrices, and let \(c_X,c_Y{\gt}0\). Suppose that, for every vertical letter \(v\), the following two corner identities hold:
Then, for every vertical letter \(v\),
Suppose a phase-class grouping has positive effective coefficients \(c_{\alpha ,q}\), corner identities
and distinguished gauges \(X_{\alpha ,0}=\mathbb {1}\). Fix a class \(\alpha \) and a copy \(q\). After identifying the representative bond space with that of copy \(q\), there are a physical map \(W_{\alpha ,q}\) and a scalar \(c_{\alpha ,0}{\gt}0\) such that, for every vertical letter \(v\),
Identify the distinguished bond space with that of copy \(q\) and set \(W_{\alpha ,q}=V_{\alpha ,0}\). Since the distinguished gauge is the identity, its corner equation reduces to
yielding (47) because \(c_{\alpha ,0}{\gt}0\).
Suppose a phase-class grouping has positive effective coefficients \(c_{\alpha ,q}\) and corner identities
Assume that the distinguished copy in each class has identity gauge and that \(M\) has the grouped-corner Gram-dressing property. Then for every class \(\alpha \), copy \(q\), and vertical letter \(v\),
Theorem 16.4.70 transports the distinguished identity-gauge corner to the bond space of copy \(q\). The grouped-corner Gram-dressing property gives
The left side is \(A_{\alpha ,q}^v\) because \(X_{\alpha ,0}=\mathbb {1}\). This gives (48).
Let \(X\) be invertible and suppose that \(B\) and a target \(C\) satisfy
Then conjugation by the Gram matrix carries \(B\) to \(C\),
Hence any two gauges satisfying this conditional identity with the same target have equal Gram conjugations, and the ratio of their Gram matrices commutes with \(B\). This gives the relative identity from an abstract common target in [ 13 , Proposition 4.13, lines 1909–1919 ] . In [ 13 , Proposition 4.13, lines 1909–1919 ] , the target appears via the reflected marked-chain argument without separate identification for either sector; the specializations \(C=B^\dagger \) and \(B^\dagger =B\) follow as corollaries.
Multiply the dressed-target identity by \(X^\dagger \) on the left and by \((X^{-1})^\dagger \) on the right. The left side becomes \((X^\dagger X)B(X^\dagger X)^{-1}\) and the right side collapses to \(C\). Eliminating the common target for two gauges and multiplying by the inverse Gram matrix gives the relative commutation relation.
Let \(M_\alpha \) be a normal tensor, and let \(X\) and \(Y\) be invertible gauges whose Gram matrices induce the same conjugation on each letter \(M_\alpha ^v\). Then there is a positive real number \(\omega \) such that
and \(\omega ^{-1/2}XY^{-1}\) is unitary. Theorem 16.4.68 gives the Figure 8 comparison in [ 13 , Proposition 4.13, lines 1909–1921 ] . This comparison supplies precisely the relative hypothesis for a sector and the distinguished sector. In the normalization \(Y=\mathbb {1}\), the resulting unitary is \(\omega ^{-1/2}X\). The common-target declarations attached to this entry are conditional algebraic corollaries, not the reflected marked-chain statement itself.
By Theorem 16.4.72, the ratio \((Y^\dagger Y)^{-1}(X^\dagger X)\) commutes with every letter of \(M_\alpha \). Theorem 16.4.65 makes this ratio scalar. When the matrix dimension is positive, positive definiteness of the two Gram matrices makes the scalar a positive real number \(\omega \). In dimension zero the equality is vacuous and one may choose \(\omega =1\). Conjugating the relative Gram identity by \(Y^{-1}\) shows that \(XY^{-1}\) has Gram matrix \(\omega \mathbb {1}\); division by \(\sqrt{\omega }\) makes it unitary.
Under the hypotheses of Theorem 16.4.71, suppose in addition that every representative \(A_{\alpha ,q}\) is normal and has positive bond dimension. Then for every class \(\alpha \) and copy \(q\) there is a real number \(\omega _{\alpha ,q}{\gt}0\) such that
Consequently \(U_{\alpha ,q}=\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}\) is unitary. This is the gauge normalization of [ 13 , Proposition 4.13, lines 1903–1921 ] . Literal CPSV canonical form and normalized BNT-refined horizontal form each supply the Gram-dressing hypothesis on their own.
Theorem 16.4.71 shows that conjugation by \(X_{\alpha ,q}^\dagger X_{\alpha ,q}\) fixes every letter of the normal representative. Apply the identity-gauge specialization of Theorem 16.4.73; it gives a real number \(\omega _{\alpha ,q}{\gt}0\) with \(X_{\alpha ,q}^\dagger X_{\alpha ,q}=\omega _{\alpha ,q}\mathbb {1}\). Dividing each gauge by the positive square root of its Gram scalar gives the unitary normalization.
After identifying the bond space of a copy with that of its representative, define its normalized physical map by
Let a phase-class grouping have gauges \(X_{\alpha ,q}\), physical maps \(V_{\alpha ,q}\), and effective coefficients \(c_{\alpha ,q}\). Suppose that, for every class \(\alpha \) and copy \(q\), \(\omega _{\alpha ,q}{\gt}0\) and
Suppose also that the maps \(V_{\alpha ,q}\) are isometries with pairwise orthogonal ranges, they intertwine \(\widetilde M\) with the corresponding gauged representatives, and their grouped corner sum reconstructs every letter of \(\widetilde M\). Then there are maps \(W_{\alpha ,q}\) such that, for every class \(\alpha \) and copy \(q\),
For distinct pairs \((\alpha ,q)\ne (\beta ,p)\),
For every vertical letter \(v\), these maps also satisfy
After identifying the copy bond space with the representative bond space, define
Since \(\omega _{\alpha ,q}{\gt}0\), the Gram identity implies
Hence \(W_{\alpha ,q}^\dagger W_{\alpha ,q}=\mathbb {1}\), and inserting the two normalized gauges on either side of the original orthogonality equation preserves pairwise orthogonality. The original intertwining identity gives
The same Gram identity and the positive square-root normalization give the corner equality
Multiplying (56) on the left by \(V_{\alpha ,q}\) and on the right by \(V_{\alpha ,q}^\dagger \) gives
Summing these physical-space identities over all copies and classes, then applying the original grouped reconstruction, proves (55).
Let a phase-class grouping have gauges \(X_{\alpha ,q}\), physical maps \(V_{\alpha ,q}\), and positive effective coefficients \(c_{\alpha ,q}\), with the distinguished gauge of every class equal to the identity. Assume that each representative is a normal tensor and that \(M\) has the grouped-corner Gram-dressing property. Put
after identifying the bond space of each copy with that of its chosen representative. Then
and every vertical letter is reconstructed exactly:
These are the normalized sector identities in [ 13 , Proposition 4.13, lines 1895–1921 ] . Literal CPSV canonical form and normalized BNT-refined horizontal form each supply the Gram-dressing hypothesis on their own.
For an MPO tensor \(M\) and a configuration \(\sigma \), the matrix product vector of the diagonal tensor \(i\mapsto M^{ii}\) equals the matrix product vector of the doubled-index tensor \(\widehat{M}\) evaluated on the diagonal-paired configuration \(k\mapsto (\sigma _k,\sigma _k)\): \(V^{(N)}(i\mapsto M^{ii})(\sigma ) =V^{(N)}(\widehat{M})(k\mapsto (\sigma _k,\sigma _k))\).
At each site value \(i\), the diagonal tensor agrees with the doubled-index tensor on the diagonal pair, \(M^{ii}=\widehat{M}^{\, (i,i)}\). The matrix product vector is multiplicative over sites, so replacing each factor \(M^{\sigma _k\sigma _k}\) by \(\widehat{M}^{\, (\sigma _k,\sigma _k)}\) reindexes the configuration from \(\sigma \) to \(k\mapsto (\sigma _k,\sigma _k)\), which is the stated equality.
Suppose the doubled-index tensor \(\widehat{M}\) shares its positive-length matrix product vector family with a block-diagonal direct sum \(\bigoplus _k\mu _kA_k\) of blocks \(A_k\) on the doubled index, as in a horizontal canonical form. Then the diagonal tensor of \(M\) generates the same positive-length matrix product vector family as the block-diagonal direct sum \(\bigoplus _k\mu _kA_k^{\mathrm{diag}}\) on the single-spin index, where \(A_k^{\mathrm{diag}}(i)=A_k(i,i)\) is the diagonal restriction of each block. For every \(N\geq 1\),
By Lemma 16.4.78, the diagonal matrix product vector equals that of \(\widehat{M}\) on the diagonal-paired configuration, which by hypothesis equals that of the block-diagonal direct sum there. Expanding the direct sum as the weighted sum \(\sum _k\mu _k^NV^{(N)}(A_k)\) and restricting each block to the diagonal pair gives \(\sum _k\mu _k^NV^{(N)}(A_k^{\mathrm{diag}})\), the matrix product vector of the single-spin block-diagonal direct sum.
There is a left-canonical, injective doubled-index tensor with a nonzero block weight satisfying the per-copy finite-length separation hypothesis whose restriction to the diagonal physical pairs is not normal. On bond space \(\mathbb {C}^2\), take
Then the four matrices \(A^{(a,b)}\) span \(M_{2}(\mathbb {C})\) and \(\sum _{a,b}(A^{(a,b)})^*A^{(a,b)}=\mathbb {1}\). For a single block, injectivity also gives finite-length block separation. The diagonal restriction, however, consists only of
Every product of these two matrices is diagonal, so no blocking length spans \(M_{2}(\mathbb {C})\).
The four nonzero scalar multiples of the matrix units span the full matrix algebra. The identity \((3/5)^2+(4/5)^2=1\) gives left-canonicality. After diagonal restriction, \(E_{ii}E_{jj}=\delta _{ij}E_{ii}\). Hence every word is either zero or a scalar multiple of some \(E_{ii}\), and therefore has vanishing \((0,1)\) entry. Since \(E_{01}\) does not, it lies outside the span at every blocking length.
For any word \(w\), the matrix product of the diagonal tensor along \(w\) equals the operator word evaluation with equal ket and bra words, since both multiply the matrices \(M^{w_kw_k}\) site by site.
By induction on the word.
For an MPO tensor \(M\) and a configuration \(\sigma \), the matrix product vector of the diagonal tensor equals the diagonal entry of the generated operator:
Both sides of (59) equal \(\operatorname{tr}(M^{\sigma _1\sigma _1}\cdots M^{\sigma _N\sigma _N})\).
Word evaluation of the diagonal tensor multiplies the matrices \(M^{\sigma _k\sigma _k}\), which is exactly the operator word evaluation with equal ket and bra indices; taking the trace gives the diagonal entry.
If \(\rho ^{(N)}(M)\) is positive semidefinite, then \(V^{(N)}(i\mapsto M^{ii})(\sigma )\geq 0\) for every configuration \(\sigma \).
By Lemma 16.4.82, the value is a diagonal entry of \(\rho ^{(N)}(M)\), and the diagonal entries of a positive semidefinite matrix are non-negative.
Let \(g\) be a block count, \(\dim ,\operatorname{mult}:\{ 0,\ldots ,g-1\} \to \mathbb {N}\), weights \(\omega _{\alpha j}\in \mathbb {C}\) for \(0\leq j{\lt}\operatorname{mult}_\alpha \), and blocks \(A_\alpha \) of bond dimension \(\dim _\alpha \). The vertical direct-sum tensor is the block-diagonal direct sum over the flattened index \((\alpha ,j)\) that repeats \(A_\alpha \) with weight \(\omega _{\alpha j}\), written \(\bigoplus _\alpha \bigoplus _j\omega _{\alpha j}A_\alpha \). Its matrix product vector decomposes as
The vertical direct-sum tensor is the block-diagonal direct sum over the single flattened index \(q\), with flattened weights \(\widetilde\omega _q\) and blocks \(\widetilde A_q\). By the block matrix-product-vector decomposition (Theorem 2.6.9),
The flattening identifies \(q\) with the pair \((\alpha ,j)\) via the bijection
under which \(\widetilde\omega _q=\omega _{\alpha j}\) and \(\widetilde A_q=A_\alpha \). Reindexing the sum along this bijection yields the double sum over \((\alpha ,j)\).
Let \((B_k)_{k\in I}\) be a finite family of blocks with weights \(\nu _k\), and let \((A_\alpha )_{\alpha \in J}\) be a finite family of representative blocks. Suppose there is a bijection
such that, writing \(e(k)=(\alpha ,j)\), one has \(\nu _k=\omega _{\alpha j}\) and, for every positive length \(N\), \(V^{(N)}(B_k)=V^{(N)}(A_\alpha )\). Then
for every positive length \(N\). This is the finite-sum reindexing after the vertical sectors have been grouped in [ 13 , Proposition 4.13 ] .
For \(N{\gt}0\) and every configuration \(\sigma \), the two sides expand as
The summand indexed by \(k\) in the first sum equals the summand indexed by \(e(k)\) in the second, so reindexing by \(e\) proves the equality.
Let \(M\) generate matrix product density operators, and suppose that its doubled-index MPS tensor is in literal CPSV canonical form. If an orthogonal projection \(P\) satisfies \(P\widetilde M_v=P\widetilde M_vP\) for every vertical letter \(v\), then
for every \(v\). Thus the vertically viewed tensor has invariant-projector closure.
The left-invariance hypothesis gives \(P_1H^{(N)}=P_1H^{(N)}P_1\) at every positive length. Since \(H^{(N)}\) is positive and \(P=P^\dagger \),
The matrix coefficients of (62) give agreement of the corresponding first-site actions. In particular,
The original-coordinate marked Lemma L gives equality of their inserted tensors on the original bond space. Reading this equality as a vertical one-site action yields
Let \(M\) be in normalized BNT-refined horizontal form and generate matrix product density operators. If an orthogonal projection \(P\) satisfies \(P\widetilde M_v=P\widetilde M_vP\) for every vertical letter \(v\), then it also satisfies \(\widetilde M_vP=P\widetilde M_vP\) for every vertical letter \(v\). Thus the vertically viewed tensor has invariant-projector closure.
Translate the left-invariance identity into the one-site ket and bra actions of the original tensor. The normalized BNT-refined version of Lemma L gives the opposite one-sided identity. Translating back to the vertically viewed tensor proves the claim.
Under the hypotheses of Theorem 16.4.87, there are positive integers \(d_k\), tensors \(B_k\) of bond dimension \(d_k\), and isometries \(V_k:\mathbb {C}^{d_k}\to \mathbb {C}^d\) such that
Every \(B_k\) is irreducible, and for every vertical letter \(v\),
The vertically viewed tensor \(\widetilde M\) and the unit-weight direct sum \(\bigoplus _kB_k\) have the same matrix product vectors at every positive length.
Apply the irreducible-corner decomposition under invariant-projector closure to \(\widetilde M\).
Under the same hypotheses, there are orthogonal projections \(P_k\) that are complete and pairwise orthogonal and, for every vertical letter \(v\), reduce the corresponding letter:
For the isometries of Theorem 16.4.88, set \(P_k=V_kV_k^\dagger \). The isometry, orthogonality, completeness, and intertwining identities give the four assertions directly.
Under the same hypotheses, there are positive integers \(d_k\), positive scalars \(\nu _k\), normal tensors \(B_k\) of bond dimension \(d_k\), and isometries \(V_k:\mathbb {C}^{d_k}\to \mathbb {C}^d\) with pairwise orthogonal ranges. For every sector \(k\) and every vertical letter \(v\),
Every vertical letter is reconstructed by
Retain exactly the irreducible corners containing a nonzero letter. The Perron–Frobenius eigenvalue of each retained corner is positive, and the absence of nontrivial periodic vectors makes the inverse-square-root rescaling normal. Put the removed square root into \(\nu _k\) and reindex the retained corners. The original complete reducing decomposition expresses each vertical letter as the sum of all compressed corners; every omitted summand vanishes, which gives (65). Thus the retained range projections need not sum to the identity, but their orthogonal complement is a zero sector.
Let \(M\) generate matrix product density operators, and suppose that its doubled-index MPS tensor is in literal CPSV canonical form. Then there are positive integers \(d_k\), positive scalars \(\nu _k\), normal tensors \(B_k\) of bond dimension \(d_k\), and isometries \(V_k:\mathbb {C}^{d_k}\to \mathbb {C}^d\) with pairwise orthogonal ranges. For every sector \(k\) and every vertical letter \(v\),
Every vertical letter is reconstructed by
The omitted orthogonal complement is an all-zero sector. This is the normal-corner intermediate between the periodic exclusion at lines 1889–1893 and the grouped vertical representation used at line 1898 of [ 13 , Proposition 4.13 ] .
Theorem 16.4.86 gives invariant-projector closure for \(\widetilde M\), and Theorem 16.4.49 excludes its nontrivial periodic vectors. Applying Theorem 7.7.5 then gives normal nonzero corners with positive spectral weights and orthogonal isometries, together with the stated intertwining and compression formulas. Their reconstruction is
An MPO tensor \(M\) has a grouped vertical BNT decomposition with physical isometries if there are positive bond dimensions \(d_k\), positive weights \(\mu _k\), normal tensors \(B_k\), and isometries \(V_k\colon \mathbb {C}^{d_k}\to \mathbb {C}^d\) with pairwise orthogonal ranges such that
Let \(k=(\alpha ,q)\) be the phase-class enumeration. There are equalities of the representative and copy bond dimensions, invertible gauges \(X_{\alpha ,q}\), and nonzero phases \(\zeta _{\alpha ,q}\) of unit modulus, with \(X_{\alpha ,0}=\mathbb {1}\) and \(\zeta _{\alpha ,0}=1\), such that
The representatives \(M_\alpha \) form a basis of normal tensors for the weighted block tensor and are pairwise inequivalent up to gauge and phase. Every effective coefficient \(c_{\alpha ,q}=\mu _{\alpha ,q}\zeta _{\alpha ,q}\) is nonzero and positive.
For \(P_{\alpha ,q}=V_{\alpha ,q}V_{\alpha ,q}^\dagger \), the triple consisting of \(P_{\alpha ,q}\), \(c_{\alpha ,q}\), and the representative loop of \(M_\alpha \) satisfies the sector-projector identity of Definition 16.4.58, and for every \((\alpha ,q)\) there is a finite chain length for which the first-site compression by \(P_{\alpha ,q}\) is nonzero. Finally, the grouped physical maps remain isometries with pairwise orthogonal ranges, satisfy both grouped intertwining identities and the grouped corner identity, and reconstruct each vertical letter exactly:
Under the same hypotheses, partition the normalized vertical sectors by equality of their matrix product vector families up to a fixed nonzero scalar per site. There are normal representatives \(M_\alpha \), an enumeration \(k=(\alpha ,q)\) of the original sectors, invertible matrices \(X_{\alpha ,q}\), and complex numbers \(\zeta _{\alpha ,q}\) of unit modulus, with \(X_{\alpha ,0}=\mathbb {1}\) and \(\zeta _{\alpha ,0}=1\), such that
The representatives are pairwise gauge-phase distinct and form a basis of normal tensors for the weighted block tensor \(\bigoplus _k\nu _kB_k\). Every grouped coefficient \(\nu _{\alpha ,q}\zeta _{\alpha ,q}\) is positive. With the physical isometries from Theorem 16.4.90 retained unchanged, one has both reducing identities, exact compression, and the literal reconstruction
For each pair \((\alpha ,q)\), the range projector of \(V_{\alpha ,q}\) obeys the representative-loop identity and gives a nonzero compression at some finite length. This is the decomposition and positivity argument of [ 13 , Proposition 4.13, lines 1898–1902 ] . The relative Gram identity and proportional-unitary normalization of lines 1903–1908 are not asserted here.
Take one representative from each phase class. Theorem 8.2.8 gives equality of bond dimensions and (69) for every copy. Distinct representatives are gauge-phase separated by construction, so Theorem 8.2.13 gives eventual linear independence. The phase-class regrouping of the finite sum gives the BNT expansion. The representative copy is assigned the identity gauge and phase one, while the self-overlap normalization of normal tensors shows that every phase has unit modulus. Substitution into the two intertwining identities, the compression identity, and the reconstruction (65) proves all remaining algebraic assertions without altering any physical isometry. Cyclicity of trace gives the sector-loop identities. The nonzero-corner separation theorem supplies a nonzero compression for every sector, and the sector trace identity and MPDO positivity then make every grouped coefficient positive.
Under the same hypotheses, partition the normalized vertical sectors by equality of their matrix product vector families up to a fixed nonzero scalar per site. There are normal representatives \(M_\alpha \), an enumeration \(k=(\alpha ,q)\) of the original sectors, invertible matrices \(X_{\alpha ,q}\), and complex numbers \(\zeta _{\alpha ,q}\) of unit modulus, with \(X_{\alpha ,0}=\mathbb {1}\) and \(\zeta _{\alpha ,0}=1\), such that
The representatives are pairwise gauge-phase distinct and form a basis of normal tensors for the weighted block tensor \(\bigoplus _k\nu _kB_k\). Every grouped coefficient \(\nu _{\alpha ,q}\zeta _{\alpha ,q}\) is positive. With the physical isometries from Theorem 16.4.91 retained unchanged, one has both reducing identities, exact compression, and the literal reconstruction
For each pair \((\alpha ,q)\), the range projector of \(V_{\alpha ,q}\) obeys the representative-loop identity and gives a nonzero compression at some finite length. This is the decomposition and positivity argument of [ 13 , Proposition 4.13, lines 1898–1902 ] . The relative Gram identity and proportional-unitary normalization of lines 1903–1908 are not asserted here.
Take one representative from each phase class. Theorem 8.2.8 gives equality of bond dimensions and (71) for every copy. Distinct representatives are gauge-phase separated by construction, so Theorem 8.2.13 gives eventual linear independence. The phase-class regrouping of the finite sum gives the BNT expansion. The representative copy is assigned the identity gauge and phase one, while the self-overlap normalization of normal tensors shows that every phase has unit modulus. Substitution into the two intertwining identities, the compression identity, and the reconstruction (66) proves all remaining algebraic assertions without altering any physical isometry. Cyclicity of trace gives the sector-loop identities. The nonzero-corner separation theorem supplies a nonzero compression for every sector, and the sector trace identity and MPDO positivity then make every grouped coefficient positive.
Let \(W_k:\mathbb {C}^{D_k}\to \mathbb {C}^d\) be a finite family of rectangular matrices. Their block-row map \(U:\mathbb {C}^d\to \bigoplus _k\mathbb {C}^{D_k}\) is defined by \(Ux=(W_k^\dagger x)_k\). This realizes the coisometry in the conclusion of [ 13 , Proposition 4.13, lines 1863–1870 ] once the sector isometries have been normalized.
Suppose that
and, for every letter \(v\) and every sector \(k\), \(T_vW_k=W_kB_k^v\). Then the block-row map \(U\) is a coisometry and conjugates \(T\) to the direct sum of the sector tensors:
If the letters also satisfy the literal sector decomposition
then
This last equality is a conditional consequence of (75); no additional assertion from the source is used.
The \((k,l)\) block of \(UU^\dagger \) is \(W_k^\dagger W_l\), so orthonormality yields (73). The corresponding block of \(UT_vU^\dagger \) is \(W_k^\dagger T_vW_l=W_k^\dagger W_lB_l^v\), which equals \(B_k^v\) for \(k=l\) and zero otherwise, yielding (74). Finally, multiplication of the block diagonal matrix by the block row gives the sum \(\sum _kW_kB_k^vW_k^\dagger \), which is \(T_v\) by hypothesis. The row indices of \(U\) form the dependent disjoint union \(\coprod _k\{ 0,\ldots ,D_k-1\} \). Its canonical order-preserving bijection with \(\{ 0,\ldots ,\sum _kD_k-1\} \) gives the coordinate reindexing required by the vertical canonical-form definition.
Let \(\{ A_\alpha \} _{\alpha =0}^{g-1}\) be a BNT for the vertically viewed tensor \(\widetilde M\). For each \(\alpha \), let \(r_\alpha {\gt}0\), let \(\omega _{\alpha ,q}{\gt}0\) for \(0\leq q{\lt}r_\alpha \), and let \(W_{\alpha ,q}:\mathbb {C}^{D_\alpha }\to \mathbb {C}^d\) satisfy
Suppose that, for every vertical letter \(v\) and every pair \((\alpha ,q)\), \(\widetilde M_vW_{\alpha ,q} =W_{\alpha ,q}(\omega _{\alpha ,q}A_\alpha ^v)\), and that the letters have the exact sector decomposition
Then \(M\) is in vertical canonical form. This is the final finite-dimensional implication in [ 13 , Proposition 4.13, lines 1863–1870 ] .
Form the block-row map \(U\) from the adjoints \(W_{\alpha ,q}^\dagger \). Theorem 16.4.96 gives
Identifying the dependent pairs \((\alpha ,q)\) and their bond indices with the standard coordinates of the two finite direct sums gives the required matrix dimensions. Positivity of the multiplicities and weights, together with the BNT hypothesis, supplies the remaining conditions of vertical canonical form.
Under the hypotheses of Theorem 16.4.93, the retained vertical decomposition has at least one sector and hence at least one phase class.
The unit-modulus weight witness in the normalized BNT-refined horizontal form, together with left canonicity of its representative, makes the doubled-index tensor nonzero. Hence its vertical reshaping \(\widetilde M\) is nonzero. If the retained decomposition had no sector, (70) would make every letter of \(\widetilde M\) zero, a contradiction. Passing from the nonempty finite set of retained sectors to its phase partition gives a retained representative.
If an MPO tensor \(M\) is in normalized BNT-refined horizontal form, then its vertical reshaping is nonzero: \(\widetilde M\ne 0\). This is the nonvanishing observation used in [ 13 , Proposition 4.13, lines 1895–1902 ] .
Write the doubled-index tensor in normalized BNT-refined horizontal form as \(\widehat M=GSG^{-1}\), where \(S\) is the weighted direct sum of its basis tensors. Choose a copy whose weight has unit modulus. Its basis tensor cannot have every letter equal to zero, since left canonicity gives \(\sum _i(A^i)^\dagger A^i=\mathbb {1}\). The corresponding nonzero weight and nonzero letter give \(S\ne 0\), hence \(\widehat M\ne 0\). If \(\widetilde M=0\), the bijection between its vertical index and the two horizontal bond indices would give \(\widehat M=0\), a contradiction.
Let \(M\) be an MPO tensor and put \(T=\widetilde M\). Let \(B_k\) have bond dimension \(d_k\), let \(\mu _k\in \mathbb {C}\), and let \(V_k:\mathbb {C}^{d_k}\to \mathbb {C}^d\) be an isometry for \(0\leq k{\lt}r\). Assume, for every vertical letter \(v\), that
Partition the tensors \(B_k\) into matrix-product-vector phase classes and choose one representative \(A_\alpha \) from each class. If \(\{ A_\alpha \} _\alpha \) is a basis of normal tensors in the spectral sense for the weighted direct sum \(A_{\mathrm{ret}}=\bigoplus _{k=0}^{r-1}\mu _kB_k\), then \(\{ A_\alpha \} _\alpha \) is an algebraic basis of normal tensors for \(T\).
This is the basis-of-normal-tensors conclusion in [ 13 , Proposition 4.13, lines 1863–1921 ] , stated for the grouped decomposition before the final coisometry is formed.
Put \(C_k^v=\mu _kB_k^v\). Induction on a word \(w\) gives \(V_k^\dagger T^w=C_k^wV_k^\dagger \). For a nonempty word \(w=(v,w')\), reconstruction of its first letter, cyclicity of trace, and \(V_k^\dagger V_k=\mathbb {1}\) therefore give
Thus \(T\) and \(A_{\mathrm{ret}}\) have the same matrix product vectors at every positive length. The spectral BNT hypothesis gives an algebraic BNT for \(A_{\mathrm{ret}}\) by Theorem 8.2.5; in particular, it supplies the required positive bond dimensions. Lemma E.1.5 then gives the same algebraic BNT for \(T\).
Let \(M\) be an MPO tensor whose doubled-index MPS tensor is in the literal CPSV canonical form, and suppose that \(M\) generates positive semidefinite operators on every nonempty chain. Let \(A\) be a tensor, let \(X\) and \(Y\) be invertible matrices, and suppose that two vertical corners satisfy
Then the Gram-dressed tensors agree:
for every vertical letter \(v\).
For \(Z=X,Y\), set
The corner identity identifies the mark with coefficients \(f^Z\) with the horizontal slice of \((Z^\dagger Z)A^v(Z^\dagger Z)^{-1}\). Self-adjointness of the two corner compressions identifies both closed marked chains with the same reflected chain. The original-coordinate marked Lemma L for literal canonical form therefore gives equality of the two marks, hence (76).
Let \(M\) be an MPO tensor whose doubled-index MPS tensor is in the literal CPSV canonical form, and suppose that \(M\) generates positive semidefinite operators on every nonempty chain. Let \(A\) be normal, let \(X\) and \(Y\) be invertible matrices, and suppose that
Then there is a real number \(\omega {\gt}0\) such that
Equation (76) makes the ratio of the two Gram matrices commute with every letter of \(A\). Normality makes this ratio scalar, and positive definiteness makes the scalar a positive real number.
Theorem 16.4.74 gives \(X_{\alpha ,q}^\dagger X_{\alpha ,q}=\omega _{\alpha ,q}\mathbb {1}\) with \(\omega _{\alpha ,q}{\gt}0\). Absorb \(\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}\) into the physical corner maps. The normalized maps are isometries with orthogonal ranges, intertwine \(\widetilde M\) with the undressed representatives, and preserve the exact reconstruction. The grouped representatives form a BNT, so Theorem 16.4.97 gives the vertical canonical form.
Let \(M\) be an MPO tensor whose doubled-index MPS tensor is in the literal CPSV canonical form. If \(M\) generates matrix product density operators, then \(M\) is in canonical form in the vertical direction. More precisely, there is a basis of normal tensors \(\{ M_\alpha \} _\alpha \) for the vertically viewed tensor \(\widetilde M\), together with positive diagonal matrices \(\mu _\alpha \) and a coisometry \(U\) from the physical space onto the retained nonzero sectors. Writing \(B=\bigoplus _\alpha \mu _\alpha \otimes M_\alpha \), one has
In [ 13 , Proposition 4.13 ] , \(U\) is called an isometry and (79) is displayed. In this rectangular orientation, \(U\) maps the physical space onto the retained sector space and satisfies (78), making \(U\) a coisometry and \(U^\dagger \) an isometric embedding. The reverse identity (80) supplies the exact reconstruction on the physical space.
A tensor in normalized BNT-refined horizontal form that generates matrix product density operators is also in canonical form in the vertical direction: there exist a basis of normal tensors \(\{ M_\alpha \} _\alpha \) for the vertically viewed tensor \(\widetilde M\), positive diagonal matrices \(\mu _\alpha \), and a coisometry \(U\) from the physical space onto the retained nonzero sectors. Writing \(B=\bigoplus _\alpha \mu _\alpha \otimes M_\alpha \), one has
This establishes vertical canonical form under the normalized BNT-refined horizontal hypothesis of Definition 16.4.1, complementing the literal CPSV formulation of Theorem 16.4.104.
16.5 Strong and weak symmetries
For a unitary on-site symmetry \(U\), a density operator \(\rho \) satisfies \([U,\rho ]=0\) (weak symmetry) or \(U\rho =\lambda \rho \) (strong symmetry) [ 190 , 198 ] . Sun formulates both notions for families of matrix product operators [ 241 ] . The density operators have a virtual boundary \(\mathcal X\) commuting with every tensor letter [ 241 ] .
For a virtual matrix \(\mathcal X\), define the operator family by
The admissible boundaries form the centralizer \(\{ \mathcal X:[\mathcal X,M^{ij}]=0\text{ for every }i,j\} \), a complex subalgebra. The family is an MPDO when \(\rho ^{(L)}(\mathcal X,M)\ge 0\) for every \(L{\gt}0\). Commutation does not imply this positivity condition. The vectorized coefficients are the same boundary-weighted traces on the doubled physical alphabet. Identity boundary recovers the ordinary periodic operator and vector families. This is the family of [ 241 , lines 175–182 ] .
An admissible boundary commutes with every product of virtual letters.
Commutation with each factor implies commutation with their product.
An operator \(O\) is a strong symmetry of a matrix \(\rho \) if \(O\rho =\lambda \rho \) for some \(\lambda \in \mathbb C\), and a weak symmetry if \([O,\rho ]=0\). A family of periodic matrix product operators \(O_a^{(L)}\) is a strong symmetry of \((\mathcal X,M)\), with eigenvalues \(\lambda _a^{(L)}\), if \(O_a^{(L)}\rho ^{(L)}=\lambda _a^{(L)}\rho ^{(L)}\) for every label \(a\) and every \(L\ge 1\). It is a weak symmetry if \([O_a^{(L)},\rho ^{(L)}]=0\) for every \(a\) and every \(L\ge 1\) [ 241 ] . The eigenvalues may depend on \(L\) and need not be phases. The definition requires neither positivity of \(\rho ^{(L)}\) nor fusion algebra relations on \(O_a\); these hypotheses are stated where used.
For a representation \(g\mapsto U_g\) on \(\mathbb C^d\), the on-site operators \(U_g^{\otimes L}\) are the periodic operators of the bond dimension one tensor with entries \((U_g)_{ij}\). Strong and weak on-site symmetry are the two notions of Definition 16.5.3 for this family. The eigenvalue \(1\) gives the strong symmetry \(U\rho =\rho \) [ 198 ] .
The periodic operator of length \(L\) of the bond dimension one tensor with entries \(u_{ij}\) is \(u^{\otimes L}\).
Each word evaluates to the scalar \(\prod _l u_{\sigma _l\tau _l}\) times the one-by-one identity, whose trace is that scalar.
With identity boundary, strong on-site symmetry with eigenvalue \(\lambda _g^{(L)}\) is the following identity of periodic contractions. By Theorem 16.5.5, the upper row is the bond dimension one tensor with entries \((U_g)_{ij}\).
Let \(\rho \) be Hermitian. If \(O\) and \(O^\dagger \) are both strong symmetries of \(\rho \), then \(O\) is a weak symmetry of \(\rho \). Hence a strong symmetry of an MPDO whose family contains the adjoint of each of its operators is a weak symmetry, and a strong on-site symmetry under a unitary representation is a weak on-site symmetry [ 241 ] .
From \(O^\dagger \rho =\mu \rho \) one gets \(\rho O=\bar\mu \rho \). Then \(\rho O\rho \) equals both \(\lambda \rho ^2\) and \(\bar\mu \rho ^2\). If \(\rho \neq 0\) then \(\rho ^2\neq 0\), so \(\lambda =\bar\mu \) and \(O\rho =\rho O\). For unitary \(O\), the relation \(\rho =O^\dagger O\rho =\lambda O^\dagger \rho \) shows that \(O^\dagger \) is a strong symmetry.
Let \(\rho ^{(L)}\neq 0\). A unitary strong symmetry of \(\rho ^{(L)}\) has an eigenvalue of modulus one. If the periodic operators obey fusion rules \(O_aO_b=\sum _c N_{ab}^cO_c\), the eigenvalues of a strong symmetry obey \(\lambda _a^{(L)}\lambda _b^{(L)}=\sum _cN_{ab}^c\lambda _c^{(L)}\), and they form a one-dimensional representation of the fusion ring when the eigenvalue of the unit is nonzero. For an on-site representation, \(g\mapsto \lambda _g^{(L)}\) is a character of the group.
For unitary \(O\), \(\rho ^\dagger \rho =(O\rho )^\dagger O\rho =|\lambda |^2\rho ^\dagger \rho \) with \(\rho ^\dagger \rho \neq 0\). The fusion relation follows by applying \(O_aO_b\) to \(\rho ^{(L)}\) in both forms.
Let \(\rho =\Phi \Phi ^\dagger \), where \(\Phi \) has system configurations as rows and ancilla configurations as columns, so that \(\Phi \) is a purification \(\Psi \) of \(\rho \). Then \(O\rho =\lambda \rho \) if and only if \((O\otimes \mathbb {1}_{\mathrm{anc}})\Psi =\lambda \Psi \). For a locally purified tensor with identity boundary the purification is the spin–ancilla matrix product state of the purifying tensor. For any virtual boundary and any supplied purification of \(\rho ^{(L)}(\mathcal X,M)\), the family is strongly symmetric under \(O_a\) if and only if every \(O_a^{(L)}\otimes \mathbb {1}_{\mathrm{anc}}\) acts on that state by \(\lambda _a^{(L)}\) [ 241 ] .
If \(O\rho =\lambda \rho \), the matrix \((O-\lambda )\Phi \) has Gram matrix \((O-\lambda )\rho (O-\lambda )^\dagger =0\), so it vanishes. The converse is immediate.
With identity boundary, let \(M\) be locally purified by \(A\). Suppose that \(\sum _j(U_g)_{ij}A^{(j,k)}=\zeta _gX_gA^{(i,k)}X_g^{-1}\) for invertible \(X_g\). Then \(M\) is strongly on-site symmetric under \(U\) with eigenvalues \(\zeta _g^L\). If instead \(\sum _{j,l}(U_g)_{ij}(V_g)_{kl}A^{(j,l)}=\zeta _gX_gA^{(i,k)}X_g^{-1}\) with \(U_g\) and \(V_g\) unitary, then \(M\) is weakly on-site symmetric under \(U\). These conditions are analogous to the criteria of de Groot, Turzillo and Schuch for channels with a symmetric unitary dilation, with trivial and with symmetric action on the ancilla [ 190 ] ; that work contains no statement about matrix product density operators.
The twisted purifying tensor is a gauge transform of \(A\) up to the scalar \(\zeta _g\), so the purifying state is multiplied by \(\zeta _g^L\). For the weak condition, \(U_g^{\otimes L}\rho U_g^{\otimes L\dagger } =|\zeta _g|^{2L}\rho \), and comparing traces gives \(|\zeta _g|^{2L}=1\) whenever \(\rho \neq 0\).
The two local conditions of Theorem 16.5.9 read diagrammatically as follows, the first for strong and the second for weak symmetry.
With identity boundary, let \(M\) be locally purified by \(A\), let the spin–ancilla tensor with entries \(A^{(i,k)}\) be normal, and let \(U\) be a unitary on-site representation. Then \(M\) is strongly on-site symmetric under \(U\) if and only if for every \(g\) there are an invertible \(X_g\) and a scalar \(\zeta _g\) with \(|\zeta _g|=1\) such that \(\sum _j(U_g)_{ij}A^{(j,k)}=\zeta _gX_gA^{(i,k)}X_g^{-1}\). In that case the eigenvalues are \(\lambda _g^{(L)}=\zeta _g^L\) at every length with \(\rho ^{(L)}\neq 0\). The symmetry thus acts on the purification by a virtual gauge transformation and trivially on the ancilla, in analogy with the strongly symmetric channel dilations of de Groot, Turzillo and Schuch [ 190 ] .
Let \(B_g\) be the twisted tensor with entries \(B_g^{(i,k)}=\sum _j(U_g)_{ij}A^{(j,k)}\). Acting with a unitary on the spin leg leaves the transfer map unchanged and preserves irreducibility, so \(B_g\) is again normal. By Theorem 16.5.8 there are scalars \(c_g(N)\) with
for every \(N{\gt}0\). Normality of \(B_g\) gives \(O_{B_gB_g}(N)\to 1\), so \(O_{B_gB_g}(N)\neq 0\) for all large \(N\), and then \(c_g(N)\neq 0\). Theorems 8.2.12 and 8.2.8 turn this proportionality into the gauge relation. The converse is Theorem 16.5.9.
Regard \(\rho ^{(L)}(\mathcal X,M)\) as the boundary-weighted vector \(|\rho \rangle \rangle \) of the tensor with doubled physical index \((i,j)\) and matrices \(M^{ij}\). Acting with \(u\otimes v\) on every doubled index sends \(\rho ^{(L)}\) to \(u^{\otimes L}\rho ^{(L)}(v^{\otimes L})^T\). For a unitary on-site representation \(U\), the operator \(M\) is weakly on-site symmetric if and only if this vector is invariant under \(g\mapsto U_g\otimes \bar U_g\), and strongly on-site symmetric with eigenvalues \(\lambda _g^{(L)}\) if and only if \((U_g\otimes \mathbb {1})^{\otimes L}|\rho \rangle \rangle =\lambda _g^{(L)}|\rho \rangle \rangle \) at every \(L\ge 1\) [ 198 ] . At identity boundary these vectors are the ordinary periodic matrix product vectors; weak symmetry is then on-site symmetry of the doubled tensor.
The boundary-weighted vector of the doubled tensor at the configuration \(((\sigma _1,\tau _1),\dots ,(\sigma _L,\tau _L))\) is the entry \(\rho ^{(L)}_{\sigma \tau }\). Expanding the trace of a product of linear combinations gives the twisted operator. For unitary \(u\), \([u,\rho ]=0\) is equivalent to \(u\rho u^\dagger =\rho \), and \((\bar u^{\otimes L})^T=(u^{\otimes L})^\dagger \).
For the operators \(\rho ^{(N)}=I^{\otimes N}+\sigma _z^{\otimes N}\) of Theorem 20.16.8.53, \(\sigma _z^{\otimes N}\) is a strong symmetry with eigenvalue \(1\) at every \(N\ge 1\). The operator \(\sigma _x^{\otimes N}\) is never a strong symmetry for \(N\ge 1\), and it is a weak symmetry exactly when \(N\) is even. In particular the family is not weakly \(\sigma _x\)-symmetric.
The operators are diagonal with entries \(1+s(\boldsymbol \sigma )\), where \(s(\boldsymbol \sigma )=\pm 1\) is the eigenvalue of \(\sigma _z^{\otimes N}\). Conjugation by \(\sigma _x^{\otimes N}\) flips every spin and multiplies \(s\) by \((-1)^N\).
Each site of the rescaling-stable example \(R\) carries two qubits. Its operators \(\rho ^{(N)}(R)\) are strongly symmetric under \((\sigma _z\otimes \sigma _z)^{\otimes N}\) with eigenvalue \(1\) and weakly symmetric under \((\sigma _x\otimes \sigma _x)^{\otimes N}\) at every \(N\ge 1\), and not strongly \(\sigma _x\otimes \sigma _x\)-symmetric.
On the support of \(\rho ^{(N)}(R)\) the second qubit of each site equals the first qubit of the next site, so the product of the \(\sigma _z\) signs is a square. Flipping every qubit preserves the support and the weight matrix. At \(N=1\) the entries of \(\rho ^{(1)}\) on the two supported configurations are in the ratio \(16:9\), and those of the flipped operator in the ratio \(9:16\).
Each site of the twisted quantum dimer \(T\) carries a left qubit, a right qubit and a flag qubit, with Pauli operators indexed by \(\ell \), \(r\) and \(f\). At every \(N\ge 1\) the operators \(\rho ^{(N)}(T)\) are strongly symmetric with eigenvalue \(1\) under \((Z_\ell Z_r)^{\otimes N}\) and under \((X_\ell X_rZ_f)^{\otimes N}\), and weakly symmetric under \((X_\ell X_r)^{\otimes N}\) and under \(Z_f^{\otimes N}\). Neither of the last two is a strong symmetry.
On the support of \(\rho ^{(N)}(T)\) the right qubit of each site equals the left qubit of the next, so the product of the \(Z_\ell Z_r\) signs is a square. Flipping the left and right qubits preserves the support and exchanges the two bond matrices, and the flag sign exchanges the two flag signs, which gives the second strong symmetry. The one-site coefficients are diagonal in the flag qubit, which gives the weak \(Z_f\) symmetry. At \(N=1\) the two flag sectors carry entries of opposite sign under both weak symmetries, so neither is strong.