Tensor Network Theory: A formalization blueprint

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

Definition 16.1.1 MPO tensor
✓
#

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,

\begin{tenkz}[rows={wire}, cols=1, west=open, east=open]
            \tn[skin=mpo, ports={90:physical:$i$, 270:physical:$j$}]{M}
        \end{tenkz}
Definition 16.1.2 Doubled-index MPS view
✓
#

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.

Definition 16.1.3 Word evaluation
✓
#

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

\begin{align} M^{u,v} & :=M^{i_1j_1}M^{i_2j_2}\cdots M^{i_Nj_N}\in M_{D}(\mathbb {C}). \label{eq:mpdo_word_eval} \end{align}

The empty pair of words evaluates to the identity, \(M^{\varnothing ,\varnothing }=\mathbb {1}_D\), and word pairs of different lengths evaluate to \(0\).

Lemma 16.1.4 Word-evaluation factorization
✓

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}\).

Proof ▶

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

\begin{align} M^{iuu_2,\, jvv_2} & =M^{ij}M^{uu_2,\, vv_2} =M^{ij}M^{u,v}M^{u_2,v_2} =M^{iu,jv}M^{u_2,v_2}, \notag \end{align}

using associativity of matrix multiplication.

Lemma 16.1.5 Trace cyclicity of the closed MPO word
✓

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)\).

Proof ▶

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,

\begin{align} \operatorname{tr}\! \left(M^{ua,\, vb}\right) & =\operatorname{tr}\! \left(M^{u,v}M^{ab}\right) =\operatorname{tr}\! \left(M^{ab}M^{u,v}\right) =\operatorname{tr}\! \left(M^{au,\, bv}\right). \notag \end{align}
Definition 16.1.6 MPO operator family
✓
#

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

\begin{align} \rho ^{(N)}(M)_{\sigma ,\tau } & =\operatorname{tr}(M^{\sigma ,\tau }). \label{eq:mpdo_family} \end{align}

Diagrammatically,

\begin{tenkz}[rows={wire}, cols=4, boundary=periodic]
        \tn[skin=mpo, ports={90:physical:$\sigma_1$, 270:physical:$\tau_1$}]{M} &
        \tn[skin=mpo, ports={90:physical:$\sigma_2$, 270:physical:$\tau_2$}]{M} &
        \tn[skin=dots]{} &
        \tn[skin=mpo, ports={90:physical:$\sigma_N$, 270:physical:$\tau_N$}]{M}
    \end{tenkz}

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.

Proof ▶

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.

Definition 16.1.8 Direct sum of MPO tensors
✓
#

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,

\begin{align} (M\oplus N)^{ik} & =M^{ik}\oplus N^{ik}. \label{eq:mpdo_direct_sum} \end{align}

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}\).

Proof ▶

Induction on the common length of \(u,v\), expanding one letter at a time and using that block-diagonal sums multiply blockwise.

Theorem 16.1.10 Additivity of the operator family
✓

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)\).

Proof ▶

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.

Definition 16.1.11 Cyclic-shift reindexing
✓
#

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\} \).

Theorem 16.1.12 Translation invariance of the periodic MPDO
✓

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 }\).

Proof ▶

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:

\begin{align} \operatorname{tr}\bigl( M^{\sigma _2\tau _2}\cdots M^{\sigma _N\tau _N}M^{\sigma _1\tau _1} \bigr) & = \operatorname{tr}\bigl( M^{\sigma _1\tau _1}\cdots M^{\sigma _N\tau _N} \bigr), \notag \end{align}

so \(\rho ^{(N)}(M)_{r\sigma ,\, r\tau } =\rho ^{(N)}(M)_{\sigma ,\tau }\).

Definition 16.1.13 Hermitian MPO tensor
✓
#

An MPO tensor is Hermitian if \(M^{ij}=(M^{ji})^\dagger \) for all \(i,j\in \{ 0,\ldots ,d{-}1\} \).

Definition 16.1.14 Adjoint MPO tensor
✓

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,

\begin{align} \rho ^{(N)}(M^\dagger )_{\sigma ,\tau } & = \overline{ \rho ^{(N)}(M)_{\tau \circ \mathrm{rev}, \sigma \circ \mathrm{rev}} }, \label{eq:mpdo_coeff_adjoint_reflection} \end{align}

where \(\mathrm{rev}\) reverses the site order. For an MPDO and \(N\geq 1\), the conjugation disappears:

\begin{align} \rho ^{(N)}(M^\dagger )_{\sigma ,\tau } & = \rho ^{(N)}(M)_{\sigma \circ \mathrm{rev}, \tau \circ \mathrm{rev}}, \label{eq:mpdo_adjoint_mpdo} \end{align}

and the adjoint tensor is again an MPDO.

Proof ▶

At the level of word evaluations, the conjugate transpose of a product reverses the factors, so

\begin{align} (M^\dagger )^{\sigma ,\tau } & = \bigl( M^{\tau \circ \mathrm{rev}, \sigma \circ \mathrm{rev}} \bigr)^\dagger . \notag \end{align}

Taking traces gives the entrywise identity

\begin{align} \rho ^{(N)}(M^\dagger )_{\sigma ,\tau } & = \operatorname{tr}\! \bigl( (M^{\tau \circ \mathrm{rev}, \sigma \circ \mathrm{rev}})^\dagger \bigr) = \overline{ \rho ^{(N)}(M)_{\tau \circ \mathrm{rev}, \sigma \circ \mathrm{rev}} }. \notag \end{align}

For an MPDO, \(\rho ^{(N)}(M)\) is Hermitian, so

\begin{align} \overline{ \rho ^{(N)}(M)_{\tau \circ \mathrm{rev}, \sigma \circ \mathrm{rev}} } & = \rho ^{(N)}(M)_{\sigma \circ \mathrm{rev}, \tau \circ \mathrm{rev}}; \notag \end{align}

positive semidefiniteness is preserved under the simultaneous relabeling of rows and columns by the reflection.

16.2 Transfer maps

Definition 16.2.1 MPO transfer map
✓
#

The transfer map associated to \(M\) is the linear map \(\mathcal{E}_M:M_{D}(\mathbb {C})\to M_{D}(\mathbb {C})\) defined by

\begin{align} \mathcal{E}_M(X) & =\sum _{i,j=0}^{d-1}M^{ij}X(M^{ij})^\dagger . \label{eq:mpdo_transfer_map} \end{align}

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}}}\).

Proof ▶

After identifying the single physical index of the doubled tensor with a pair \((i,j)\), for every \(X\in M_{D}(\mathbb {C})\),

\begin{align} \mathcal{E}_{M^{\mathrm{MPS}}}(X) & = \sum _{(i,j)} M^{\mathrm{MPS}}_{(i,j)}X (M^{\mathrm{MPS}}_{(i,j)})^\dagger = \sum _{i,j=0}^{d-1}M^{ij}X(M^{ij})^\dagger = \mathcal{E}_M(X). \notag \end{align}
Theorem 16.2.3 Positivity of the MPO transfer map
✓

If \(X\ge 0\), then \(\mathcal{E}_M(X)\ge 0\).

Proof ▶

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

Definition 16.3.1 MPDO
✓
#

An MPO tensor is an MPDO if, for every \(N\geq 1\),

\begin{align} \rho ^{(N)}(M) & \ge 0. \label{eq:mpdo_positivity} \end{align}

This is the global positivity condition of [ 13 , Section 4 ] .

Definition 16.3.2 LPDO
✓
#

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\),

\begin{align} M^{ij} & = \left(\sum _{k=0}^{d_K-1} A^{(i,k)}\otimes \overline{A^{(j,k)}}\right)_{e,e}, \label{eq:mpdo_lpdo_factorization} \end{align}

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

\begin{tenkz}[rows={wire}, cols=1, west=open, east=open]
            \tn[skin=mpo, ports={90:physical:$i$, 270:physical:$j$}]{M}
        \end{tenkz} \( = \) \begin{tenkz}[rows={op,op}, cols=3]
            \tnfuse[at=(1,1), name=e, skin=pill]{e} \tnwire{open w}{e.180}
            \tn[at=(1,2), skin=box, ports={90:physical:$i$}]{A}
            \tn[at=(2,2), skin=box, ports={270:physical:$j$}]{\overline{A}}
            \tnfuse[at=(1,3), name=ei, skin=pill]{e^{-1}} \tnwire{ei.0}{open e}
        \end{tenkz}
Proof ▶

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\):

\begin{align} M^{\sigma _1\tau _1}\cdots M^{\sigma _N\tau _N} & = \left(\sum _{\kappa } P_\kappa \otimes \overline{Q_\kappa }\right)\big|_e, \notag \end{align}

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

\begin{align} \rho ^{(N)}(M)_{\sigma ,\tau } & = \sum _\kappa \operatorname{tr}(P_\kappa )\overline{\operatorname{tr}(Q_\kappa )} = \sum _\kappa \psi _\kappa (\sigma )\overline{\psi _\kappa (\tau )}, \notag \end{align}

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

Definition 16.4.1 Normalized BNT-refined horizontal form for an MPO
✓

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:

\begin{align} M^{ab} & =\bigoplus _{j=0}^{g-1}\bigoplus _{q=0}^{r_j-1} \mu _{j,q}X_{j,q}A_j^{ab}X_{j,q}^{-1}. \label{eq:mpdo_horizontal_cf} \end{align}

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}\),

\begin{align} M^{ab} & =X\left[\bigoplus _{j,q}\mu _{j,q}A_j^{ab}\right]X^{-1}, \notag \\ X & =\bigoplus _{j,q}X_{j,q}. \notag \end{align}

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.

Definition 16.4.2 MPV representation associated with horizontal canonical form
✓

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

\begin{align} \sum _{j=0}^{g-1} \left(\sum _{q=0}^{r_j-1}\mu _{j,q}^{N}\right)V^{(N)}(A_j). \label{eq:mpdo_horizontal_mpv} \end{align}

This is the MPV-level decomposition in [ 13 , Section 2.3, lines 298–301 ] .

Theorem 16.4.3 BNT-refined horizontal form gives the BNT MPV representation
✓

Every MPO tensor in normalized BNT-refined horizontal form has the associated positive-length BNT matrix-product-vector representation.

Proof ▶

Each copy similarity cancels inside every positive-length closed matrix-product trace.

\(K=X\left(\bigoplus _j\mathcal K_j\right)X^{-1}\)
\begin{tenkz}[rows={wire}, cols=1, west=open, east=open]
            \tn[ports={90:physical, 270:physical}]{K}
        \end{tenkz} \( = \) \begin{tenkz}[rows={wire}, cols=3, west=open, east=open]
            \tn[skin=ring]{X} & \tn[ports={90:physical, 270:physical}]{\bigoplus_j\mathcal K_j} &
            \tn[skin=ring]{X^{-1}}
        \end{tenkz}
Figure 16.1 Horizontal canonical form of the doubled tensor: the gauge \(X\), the direct sum of basis tensors \(\mathcal K_j\) carrying the ket and bra indices, and \(X^{-1}\). In the source specialized simple-biCF/ZCL display the weights \(\mu _j\) are independent of the copy index and absorbed into the basis tensors, so the basis tensors depend only on the block index \(j\) while the gauge factors \(X_{j,q}\) depend on both block indices. This is [ 13 , Appendix C.2, lines 1660–1665 ] .
Definition 16.4.4 Vertically viewed tensor of an MPO
✓

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:

\begin{align} \widetilde{M^\dagger }_{a,b} & =(\widetilde M_{b,a})^\dagger . \label{eq:mpdo_adjoint_reflection} \end{align}
Proof ▶

This follows entrywise from the exchange of the ket and bra indices and conjugate transposition of the virtual matrix.

Definition 16.4.6 Vertical canonical form for an MPO
✓

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

\begin{align} U\widetilde MU^\dagger & =B, \notag \\ \widetilde M & =U^\dagger B U, \notag \\ B & =\bigoplus _\alpha \mu _\alpha \otimes M_\alpha . \label{eq:mpdo_vertical_cf} \end{align}

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.

\begin{tenkz}[rows={op,op,op}, cols=3]
        \tn[at=(1,1), skin=pill, wide=3, physical=up]{U}
        \tn[at=(2,1), skin=box, wide=3, name=M]{\widetilde M}
        \tn[at=(3,1), skin=pill, wide=3, physical=down]{U^\dagger}
        \tnwire{open w}{M.180} \tnwire{M.0}{open e}
    \end{tenkz} \( = \bigoplus _\alpha \) \begin{tenkz}[rows={wire,wire}, cols=4, bonds=none]
                \tn[at=(1,1), name=d1,
                    ports={90:physical, 270:physical, 0:physical}]{}
                \tn[at=(1,2), skin=box, name=mu,
                    ports={180:physical, 0:physical}]{\mu_\alpha}
                \tn[at=(1,3), name=d2,
                    ports={90:physical, 270:physical, 180:physical}]{}
                \tn[at={midway (1,3) and (2,3)}, name=d3,
                    ports={90:physical, 270:physical, 0:physical}]{}
                \tn[at=(2,4), skin=box, name=M,
                    ports={180:virtual, 0:virtual, 90:physical, 270:physical,
                            158:physical}]{M_\alpha}
                \tnwire{d1.90}{open n}
                \tnwire{d1.0}{mu.180}
                \tnwire{d1.270}{open s}
                \tnwire{d2.90}{open n}
                \tnwire{mu.0}{d2.180}
                \tnwire{d3.90}{d2.270}
                \tnwire{d3.270}{open s}
                \tnwire{d3.0}{M.158}
            \end{tenkz}
Figure 16.2 Vertical direct-sum form of the tensor: after the physical isometry \(U\), each block factorizes as the positive diagonal weight \(\mu _\alpha \) on the multiplicity factor times the basis tensor \(M_\alpha \), and only \(M_\alpha \) carries the virtual legs. This is the decomposition \(U\widetilde MU^\dagger =\bigoplus _\alpha \mu _\alpha \otimes M_\alpha \) of [ 13 , Proposition 4.13, lines 1864–1870 ] .

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,

\begin{align} \operatorname{span}_{\mathbb {C}}\bigl\{ \, (A_k^w)_k : w \in \{ 0,\ldots ,d{-}1\} ^L\, \bigr\} & =\textstyle \prod _k M_{D_k}(\mathbb {C}). \label{eq:mpdo_word_tuple_span} \end{align}

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,

\begin{align} \Bigl(\forall \, w \in \{ 0,\ldots ,d{-}1\} ^L, \textstyle \sum _k \operatorname{tr}(\Delta _k A_k^w) = 0\Bigr) & \implies \forall k,\ \Delta _k = 0. \label{eq:mpdo_finite_separation} \end{align}

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

\begin{align} f\bigl((M_k)_k\bigr) & =\sum _k\operatorname{tr}(\Delta _kM_k), \notag \end{align}

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

\begin{align} \sum _w c_w A_k^w & =\mathbb {1}, \notag \\ \sum _w c_w A_j^w & =0 \quad (j\neq k). \notag \end{align}

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

\begin{align} \sum _p(\mathcal K_j^p)_{\alpha ,\beta } (K^{-1,p}_{j'})_{\alpha ',\beta '} & =\delta _{\alpha ,\alpha '}\delta _{\beta ,\beta '}\delta _{j,j'}. \label{eq:mpdo_inverse_tensor} \end{align}

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\):

\begin{tenkz}[rows={op,op}, cols=1, west=open, east=open]
        \tn[skin=box, ports={90:physical}]{K^{-1}} \\
        \tn[skin=box]{K}
    \end{tenkz} \( = \) \begin{tenkz}[rows={wire,wire,wire}, cols=3, bonds=none]
        \tn[at=(1,1), skin=box, name=X,
            ports={180:virtual, 0:virtual, 270:virtual}]{X}
        \tn[at=(1,2), skin=box, name=K,
            ports={180:virtual, 0:virtual, 90:physical,
                    270:virtual}]{\mathcal K_j}
        \tn[at=(1,3), skin=box, name=Xinv,
            ports={180:virtual, 0:virtual, 270:virtual}]{X^{-1}}
        \tn[at={1 n of (1,2)}, skin=box, name=Kinv,
            ports={180:virtual, 0:virtual, 90:physical,
                    270:physical}]{K^{-1}}
        \tn[at=(2,1), name=jX,
            ports={180:virtual, 0:virtual, 90:virtual, 270:virtual}]{}
        \tn[at=(2,2), name=jK, ports={180:virtual, 0:virtual, 90:virtual}]{}
        \tn[at=(2,3), name=jXinv,
            ports={180:virtual, 0:virtual, 90:virtual, 270:virtual}]{}
        \tn[at=(3,1), name=qX, ports={180:virtual, 0:virtual, 90:virtual}]{}
        \tn[at=(3,3), name=qXinv,
            ports={180:virtual, 0:virtual, 90:virtual}]{}
        \tnwire{X.180}{open w}
        \tnwire{X.0}{K.180}
        \tnwire{K.0}{Xinv.180}
        \tnwire{Xinv.0}{open e}
        \tnwire{Kinv.90}{open n}
        \tnwire{Kinv.270}{K.90}
        \tnwire{X.270}{jX.90}
        \tnwire{K.270}{jK.90}
        \tnwire{Xinv.270}{jXinv.90}
        \tnwire{jX.270}{qX.90}
        \tnwire{jXinv.270}{qXinv.90}
        \tnwire[name=jw]{jX.180}{open w}
        \tnwire{jX.0}{jK.180}
        \tnwire{jK.0}{jXinv.180}
        \tnwire{jXinv.0}{open e}
        \tnwire[name=qw]{qX.180}{open w}
        \tnwire{qX.0}{qXinv.180}
        \tnwire{qXinv.0}{open e}
        \tnmark[form=label, label pos=w]{on jw 0.8}{$j$}
        \tnmark[form=label, label pos=w]{on qw 0.8}{$q$}
    \end{tenkz} \( = \) \begin{tenkz}[rows={wire,wire,wire}, cols=3, bonds=none]
        \tn[at=(1,1), skin=box, name=X,
            ports={180:virtual, 45:virtual, 270:virtual}]{X}
        \tn[at=(1,3), skin=box, name=Xinv,
            ports={135:virtual, 0:virtual, 270:virtual}]{X^{-1}}
        \tn[at=(2,1), name=jX,
            ports={180:virtual, 0:virtual, 90:virtual, 270:virtual}]{}
        \tn[at=(2,2), name=delta,
            ports={180:virtual, 0:virtual, 90:physical}]{}
        \tn[at=(2,3), name=jXinv,
            ports={180:virtual, 0:virtual, 90:virtual, 270:virtual}]{}
        \tn[at=(3,1), name=qX, ports={180:virtual, 0:virtual, 90:virtual}]{}
        \tn[at=(3,3), name=qXinv,
            ports={180:virtual, 0:virtual, 90:virtual}]{}
        \tnwire{X.180}{open w}
        \tnwire{Xinv.0}{open e}
        \tnwire{delta.90}{open n}
        \tnwire{X.270}{jX.90}
        \tnwire{Xinv.270}{jXinv.90}
        \tnwire{jX.270}{qX.90}
        \tnwire{jXinv.270}{qXinv.90}
        \tnwire[name=jw]{jX.180}{open w}
        \tnwire{jX.0}{delta.180}
        \tnwire{delta.0}{jXinv.180}
        \tnwire{jXinv.0}{open e}
        \tnwire[name=qw]{qX.180}{open w}
        \tnwire{qX.0}{qXinv.180}
        \tnwire{qXinv.0}{open e}
        \tnmark[form=label, label pos=w]{on jw 0.8}{$j$}
        \tnmark[form=label, label pos=w]{on qw 0.8}{$q$}
    \end{tenkz}

Consequently, \(K^{-1}K=\bigoplus _j\mathbb {1}_j\). This is the inverse identity in [ 13 , Appendix C.2, lines 1666–1670 ] .

Theorem 16.4.8 Finite-word trace-separation criterion
✓

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.

Proof ▶

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.

Proof ▶

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

\begin{align} \operatorname{span}_{\mathbb {C}}\bigl\{ (A_j^w)_{j=0}^{g-1}: w\in \{ 0,\ldots ,d{-}1\} ^{L}\bigr\} & =\prod _{j=0}^{g-1}M_{D_j}(\mathbb {C}). \label{eq:mpdo_rep_tuple_sep} \end{align}

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 ] .

Proof ▶

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

\begin{align} \sum _{j=0}^{g-1} \operatorname{tr}\! \left(c_j^{(L+1)}(YA_j-ZA_j)A_j^w\right) & =0. \notag \end{align}

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.

Definition 16.4.12 Direct sum with one marked letter
✓

Let

\begin{align} A_{\mathrm{tot}}^i & =\bigoplus _{j=0}^{g-1}\bigoplus _{q=0}^{r_j-1} \mu _{j,q}A_j^i, \label{eq:mpdo_sector_direct_sum} \end{align}

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

\begin{align} C_{\mathrm{tot}}^s & =\bigoplus _{j=0}^{g-1}\bigoplus _{q=0}^{r_j-1} \mu _{j,q}C_j^s. \label{eq:mpdo_marked_direct_sum} \end{align}

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 ] .

Lemma 16.4.13 Trace expansion with one marked letter
✓

For every marked index \(s\) and every unmarked word \(w\),

\begin{align} \operatorname{tr}\! \left(C_{\mathrm{tot}}^s A_{\mathrm{tot}}^w\right) & =\sum _{j=0}^{g-1}c_j^{(|w|+1)} \operatorname{tr}\! \left(C_j^s A_j^w\right), \label{eq:mpdo_marked_trace_expansion} \\ c_j^{(N)} & =\sum _{q=0}^{r_j-1}\mu _{j,q}^N. \notag \end{align}

The exponent \(|w|+1\) includes the weight at the marked site.

Proof ▶

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\),

\begin{align} \operatorname{tr}\! \left(C_{\mathrm{tot}}^s A_{\mathrm{tot}}^w\right) & =\operatorname{tr}\! \left(E_{\mathrm{tot}}^s A_{\mathrm{tot}}^w\right). \notag \end{align}

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.

Proof ▶

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

\begin{align} \operatorname{tr}\! \left(\widehat C_{\mathrm{grp}}^sT_{\mathrm{grp}}^w\right) & =\operatorname{tr}\! \left(C_{\mathrm{rep}}^sP_{\mathrm{rep}}^w\right) \notag \\ & =\sum _j\left(\sum _q\lambda _{j,q}^{|w|+1}\right) \operatorname{tr}\! \left(C_j^sC_j^w\right). \label{eq:cpsv_full_coordinate_marked_trace_identity} \end{align}

The exponent \(|w|+1\) includes the factor at the marked site.

Proof ▶

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

\begin{align} \operatorname{tr}\! \left(\widehat C_{\mathrm{grp}}^sT_{\mathrm{grp}}^w\right) & =\operatorname{tr}\! \left(\widehat E_{\mathrm{grp}}^sT_{\mathrm{grp}}^w\right) \label{eq:cpsv_retained_marked_trace_agreement} \end{align}

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

\begin{align} (YC_j)^i & =(ZC_j)^i \label{eq:cpsv_retained_first_site_action} \end{align}

for every representative \(j\) and every physical letter \(i\).

Proof ▶

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

\begin{align} (YC_j)^i & =(ZC_j)^i. \label{eq:cpsv_prop_4_13_active_separation} \end{align}

Let \(A\) be in literal CPSV canonical form. For canonical-form data and its BNT refinement, write

\begin{align} A^i & =V^\dagger B^iV, \notag \\ VV^\dagger & =\mathbb {1}. \label{eq:cpsv_original_coisometry_reconstruction} \end{align}

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,

\begin{align} \operatorname{tr}(CA^w) & =\operatorname{tr}(EB^w). \label{eq:cpsv_original_marked_trace} \end{align}

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

\begin{align} \operatorname{tr}\! \left(F_A^uA^w\right) & =\operatorname{tr}\! \left(G_A^uA^w\right) \label{eq:cpsv_original_linear_marked_agreement} \end{align}

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.

Proof ▶

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\).

Definition 16.4.18 First-site action on a physical block
✓
#

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

\begin{align} Y^{[L+1]}_{I,J} & =Y_{i_0,j_0}\prod _{k=1}^{L}\delta _{i_k,j_k}. \label{eq:mpdo_first_site_block} \end{align}

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

\begin{align} (Y^{[L+1]}A^{[L+1]})^I & =(YA)^{i_0}A^{i_1}\cdots A^{i_L}. \label{eq:mpdo_blocked_insertion} \end{align}

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

\begin{align} \operatorname{tr}((YA)^i A^w) & =\operatorname{tr}((ZA)^i A^w), \label{eq:mpdo_marked_first_site_trace} \end{align}

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\).

Proof ▶

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.

Proof ▶

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

\begin{align} \sum _{|w|=6(g-1)}c_w^{(j)}A_k^w & = \begin{cases} \mathbb {1}, & k=j, \\ 0, & k\neq j. \end{cases}. \label{eq:mpdo_rep_selectors} \end{align}

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,

\begin{align} \lim _{k\to \infty }O_{A_jA_j}(m_jk) & =m_j, \notag \\ \lim _{k\to \infty }O_{A_jA_j}(m_jk) & =1. \notag \end{align}

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\).

Proof ▶

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 ] .

Proof ▶

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

\begin{align} C_{\mathrm{tot}}^u & =\bigoplus _{j,q}\mu _{j,q} \sum _z f_{u,z}A_j^z, \notag \\ E_{\mathrm{tot}}^u & =\bigoplus _{j,q}\mu _{j,q} \sum _z g_{u,z}A_j^z. \notag \end{align}

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

\begin{align} YA & =ZA \label{eq:cpsv_literal_first_site_inserted} \end{align}

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.

Proof ▶

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

\begin{align} (YA_k)^i & := \sum _{j=0}^{d-1}Y_{i,j}A_k^j \notag \end{align}

for the physical-index insertion of \(Y\) into block \(A_k\), and for a configuration \(\sigma =(i_1,\ldots ,i_N)\) write

\begin{align} (Y_1V^{(N)}(A))_{(i_1,\ldots ,i_N)} & := \sum _{j=0}^{d-1}Y_{i_1,j} V^{(N)}(A)_{(j,i_2,\ldots ,i_N)} \notag \end{align}

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)\)

\begin{tenkz}[cols=4, boundary=periodic, physical=up]
            \tn[name=M, ports={90:physical}]{\widetilde M} & \tn{\widetilde M} &
            \tn[skin=dots]{} & \tn{\widetilde M}
            \tn[at=1 n of M, name=Y, label pos=w, ports={270:physical, 90:physical}]{Y}
            \tnwire{Y.270}{M.90}
        \end{tenkz} \( = \) \begin{tenkz}[cols=4, boundary=periodic, physical=up]
            \tn[name=M, ports={90:physical}]{\widetilde M} & \tn{\widetilde M} &
            \tn[skin=dots]{} & \tn{\widetilde M}
            \tn[at=1 n of M, name=Z, label pos=w, ports={270:physical, 90:physical}]{Z}
            \tnwire{Z.270}{M.90}
        \end{tenkz}

\(YA_k=ZA_k\)

\begin{tenkz}[rows={op,ket}, cols=1]
            \tn[label pos=w, ports={90:physical}]{Y} \\
            \tn[ports={180:virtual, 0:virtual}]{A_k}
        \end{tenkz} \( = \) \begin{tenkz}[rows={op,ket}, cols=1]
            \tn[label pos=w, ports={90:physical}]{Z} \\
            \tn[ports={180:virtual, 0:virtual}]{A_k}
        \end{tenkz}

Figure 16.3 Blockwise equality from agreement on the MPV family. Top: the hypothesis \(Y_1V^{(N)}(A)=Z_1V^{(N)}(A)\) for every positive chain length \(N\geq 1\), with the operator inserted on the physical leg of the first site. Bottom: the conclusion that the two insertions agree on every block \(A_k\) of the canonical form. These are the two displays of Lemma L, [ 13 , Appendix C.3, lines 1835–1846 ] .
Proof ▶

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\),

\begin{align} \sum _k \operatorname{tr}(\mu _k^{L+1}\Delta _k A_k^w) & =0. \notag \end{align}

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)\):

\begin{tenkz}[cols=2, boundary=periodic]
        \tn{W} & \tn[name=A, ports={90:physical}]{A_k}
        \tn[at=1 n of A, name=Y, label pos=w, ports={270:physical, 90:physical}]{Y}
        \tnwire{Y.270}{A.90}
    \end{tenkz} \( = \) \begin{tenkz}[cols=2, boundary=periodic]
        \tn{W} & \tn[name=A, ports={90:physical}]{A_k}
        \tn[at=1 n of A, name=Z, label pos=w, ports={270:physical, 90:physical}]{Z}
        \tnwire{Z.270}{A.90}
    \end{tenkz}

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 ] .

Theorem 16.4.25 One-sided invariant corners of Hermitian operators
✓

Let \(H\) and \(P\) be Hermitian matrices. If \(PH=PHP\), then \((\mathbb {1}-P)HP=0\).

Proof ▶

Taking adjoints gives \(HP=PHP\), and the conclusion follows by subtraction.

Theorem 16.4.26 One-sided invariant corners of positive operators
✓

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 ] .

Proof ▶

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\),

\begin{align} (Y_L)_1V^{(N)}(A) & =(Y_C)_1V^{(N)}(A), \notag \\ (Y_L)_1V^{(N)}(A) & =(Y_R)_1V^{(N)}(A), \notag \end{align}

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.

Proof ▶

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}\).

Definition 16.4.28 Doubled-index forms of the three first-site contractions
✓

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

\begin{align} (Y_L)_{p,q} & =P_{a,c}\delta _{b,e}, \notag \\ (Y_R)_{p,q} & =\delta _{a,c}P_{e,b}, \notag \\ (Y_C)_{p,q} & =P_{a,c}P_{e,b}. \notag \end{align}

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)\):

\(H^{(N)}(\widetilde M) = \) \begin{tenkz}[cols=4, boundary=periodic]
        \tn[label pos=225, ports={90:physical:$\sigma_1$, 270:physical:$\tau_1$}]{\widetilde M_1} &
        \tn[label pos=225, ports={90:physical:$\sigma_2$, 270:physical:$\tau_2$}]{\widetilde M_2} &
        \tn[skin=dots]{} &
        \tn[label pos=225, ports={90:physical:$\sigma_N$, 270:physical:$\tau_N$}]{\widetilde M_N}
    \end{tenkz}

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 ] .

Lemma 16.4.29 Coefficient form of the three first-site contractions
✓

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,

\begin{align} \sum _c P_{a,c}H^{(N+1)}_{(c,\sigma '),(b,\tau ')} & , \label{eq:mpdo_contraction_coefficients} \\ \sum _e H^{(N+1)}_{(a,\sigma '),(e,\tau ')}P_{e,b} & , \notag \\ \sum _{c,e}P_{a,c}H^{(N+1)}_{(c,\sigma '),(e,\tau ')}P_{e,b} & . \notag \end{align}
Proof ▶

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\).

Proof ▶

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}\).

\(PH^{(N+1)} = \) \begin{tenkz}[cols=4, boundary=periodic, physical=updown, size=s]
        \tn[skin=mpo]{\widetilde M} & \tn[skin=mpo]{\widetilde M} & \tn[skin=dots]{} &
        \tn[skin=mpo, name=M, physical=down, ports={90:physical}]{\widetilde M}
        \tn[at=1 n of M, name=P, label pos=w, ports={270:physical, 90:physical}]{P}
        \tnwire{P.270}{M.90}
    \end{tenkz}
Figure 16.4 The periodic contraction representing \(PH^{(N+1)}\). The analogous contractions with \(Y_R\) and \(Y_C\) place \(P\) on the bra leg or on both legs and represent \(H^{(N+1)}P\) and \(PH^{(N+1)}P\). These are the three first-site contractions in [ 13 , Proposition 4.13, lines 1881–1887 ] .
Theorem 16.4.31 Representative reduction for the MPO contractions
✓

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,

\begin{align} PH^{(N)} & =PH^{(N)}P, \notag \\ PH^{(N)} & =H^{(N)}P, \notag \end{align}

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}\).

Proof ▶

Lemma 16.4.29 turns the two operator identities into equality of the corresponding first-site contractions. Equality of matrix product vectors gives these equalities for the chosen canonical-form decomposition. Apply Theorem 16.4.27 with \((Y_L,Y_C,Y_R)\).

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

\begin{align} (PM)^{ij} & =\sum _k P_{ik}\, M^{kj}, \notag \\ (MP)^{ij} & =\sum _k P_{kj}\, M^{ik}, \notag \end{align}

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.

Proof ▶

Put \(L=N+1\). Splitting off the first site of a density-operator entry gives

\begin{align} H^{(N+1)}_{(i,\sigma ),(j,\tau )} & =\operatorname{tr}\bigl(M^{ij}M^{\sigma _0\tau _0}\cdots M^{\sigma _{N-1}\tau _{N-1}}\bigr). \notag \end{align}

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.

Proof ▶

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.

Theorem 16.4.35 Commutation with a power of a positive semidefinite matrix
✓
#

Let \(H\geq 0\), let \(p\geq 1\), and let \(QH^p=H^pQ\). Then \(QH=HQ\).

Proof ▶

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.

Theorem 16.4.36 Commutation with a power of a positive semidefinite matrix
✓

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 ] .

Proof ▶

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\),

\begin{align} Q’_{ij} \lambda _j^{\, p} & =\lambda _i^{\, p}\, Q’_{ij}. \notag \end{align}

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\).

Definition 16.4.37 Periodic-sector projector
✓

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,

\begin{align} Q_1H^{(N+1)} & \ne H^{(N+1)}Q_1, \label{eq:mpdo_periodic_noncommute} \\ Q_1[H^{(N+1)}]^p & =[H^{(N+1)}]^pQ_1. \label{eq:mpdo_periodic_commute_power} \end{align}

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.

Theorem 16.4.38 No periodic-sector projector for a matrix product density operator
✓

Let \(M\) generate matrix product density operators. Then \(M\) admits no periodic-sector projector.

Proof ▶

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.

Definition 16.4.39 Periodic vectors supply projectors
✓

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.

Proof ▶

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.

Definition 16.4.41 Vertically stacked layers
✓

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.

Proof ▶

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.

Proof ▶

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\).

Proof ▶

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.

Definition 16.4.45 Periodic vectors supply cyclic projectors
✓

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.

Theorem 16.4.46 Reduction of the periodic-sector step to the cyclic decomposition
✓

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.

Proof ▶

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.

Proof ▶

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.

Theorem 16.4.48 A displaced idempotent fails to commute at some length
✓

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}\),

\begin{align} Q_1H^{(N+1)} & \ne H^{(N+1)}Q_1. \label{eq:cpsv_literal_displaced_noncommutation} \end{align}

This existential non-commutation statement supplies the single length needed for the contradiction in [ 13 , Proposition 4.13, lines 1888–1893 ] .

Proof ▶

Suppose instead that \(Q_1H^{(N+1)}=H^{(N+1)}Q_1\) for every \(N\geq 0\). Idempotence gives

\begin{align} Q_1H^{(N+1)} & =Q_1^2H^{(N+1)} =Q_1H^{(N+1)}Q_1. \notag \end{align}

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.

Proof ▶

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

\begin{align} Q_1H^{(N+1)} & \ne H^{(N+1)}Q_1. \notag \end{align}

The word invariance and the vertical-stack identity give

\begin{align} Q_1[H^{(N+1)}]^p & =[H^{(N+1)}]^pQ_1. \notag \end{align}

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 ] .

Proof ▶

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.

Theorem 16.4.51 Commutation at chain length two forces letter invariance, for an injective tensor
✓

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\).

Proof ▶

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.

Theorem 16.4.52 No periodic vectors for an injective matrix product density operator
✓

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.

Proof ▶

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.

Theorem 16.4.53 Positive trace of nonzero positive semidefinite matrices
✓
#

A nonzero positive semidefinite matrix has positive trace: the trace is a positive real.

Proof ▶

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.

Proof ▶

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.

Lemma 16.4.56 The trace of the density operator is a chain of loops
✓

Summing the diagonal word evaluations of an MPO tensor over all length-\(N\) configurations closes the vertical loops site by site,

\begin{align} \sum _{\sigma } M^{\sigma _0\sigma _0}\cdots M^{\sigma _{N-1}\sigma _{N-1}} & =E^N, \notag \end{align}

and hence \(\operatorname{tr}H^{(N)}=\operatorname{tr}(E^N)\).

Proof ▶

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\),

\begin{align} \operatorname{tr}(P_1H^{(N+1)}) & =\operatorname{tr}(E_P E^N). \label{eq:mpdo_first_site_trace} \end{align}

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)\).

Proof ▶

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

\begin{align} (P_1H^{(N+1)})_{(a,\rho )(a,\rho )} & =\sum _i P_{ai} \operatorname{tr}(M^{ia}E_\rho ^{(N)}). \notag \end{align}

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)})\).

Definition 16.4.58 Sector projectors of the vertical decomposition
✓

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

\begin{align} \widetilde M & =\bigoplus _{\alpha ,k} \mu _{\alpha ,k}X_{\alpha ,k}M_\alpha X_{\alpha ,k}^{-1}, \notag \end{align}

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:

\begin{align} \left(\exists v,\; P\widetilde M_vP\ne 0\right) & \Longrightarrow \left(\exists N\geq 0,\; P_1H^{(N+1)}P_1\ne 0\right), \label{eq:cpsv_literal_sector_compression_separation} \end{align}

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.

Proof ▶

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

\begin{align} P\widetilde M_vP & =0. \label{eq:cpsv_literal_sector_compression_zero_corner} \end{align}

This contradicts the hypothesis in (39).

Theorem 16.4.60 A gauge corner gives a nonzero range-projector corner
✓

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\),

\begin{align} \omega X A_v X^{-1} & =V^\dagger \widetilde M_vV. \label{eq:mpdo_gauge_corner} \end{align}

Then the range projector \(VV^\dagger \) has a nonzero vertical corner:

\begin{align} (VV^\dagger )\widetilde M_v(VV^\dagger ) & \ne 0 \label{eq:mpdo_range_projector_corner} \end{align}

for some \(v\). This is the algebraic nonvanishing step used in [ 13 , Proposition 4.13, line 1898 ] .

Proof ▶

Normality supplies a letter \(v\) for which \(A_v\ne 0\). Since \(\omega \) is nonzero and \(X\) is invertible, the left-hand side of (41) is nonzero. If (42) vanished, multiplying it by \(V^\dagger \) and \(V\) and using \(V^\dagger V=1\) would contradict (41).

Theorem 16.4.61 A nonzero vertical corner has a nonzero sector compression
✓

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.

Proof ▶

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.

Theorem 16.4.62 The sector trace identity
✓

Let \(P\) be a sector projector for the weight \(\mu \) and the matrix \(E_\alpha \). Then for every chain length \(N+1\),

\begin{align} \operatorname{tr}(P_1H^{(N+1)}P_1) & =\mu \, d_\alpha ^{(N+1)}. \label{eq:mpdo_sector_trace} \end{align}

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.

Proof ▶

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

\begin{align} \mu _{\alpha ,k}\, d_\alpha ^{(N+1)} & {\gt}0, \notag \\ d_\alpha ^{(N+1)} & {\gt}0, \notag \\ \mu _{\alpha ,k} & {\gt}0. \notag \end{align}

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\).

Proof ▶

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\).

Theorem 16.4.64 Unitary normalization of a scalar Gram matrix
✓

If \(\omega {\gt}0\) and \(X^\dagger X=\omega \mathbb {1}\), then \(\omega ^{-1/2}X\) is unitary.

Proof ▶

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.

Proof ▶

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 \),

\begin{align} \operatorname{tr}\! \left(\mathcal{R}(A)^{rs} M^{\sigma _1\tau _1}\cdots M^{\sigma _L\tau _L}\right) & = \operatorname{tr}\! \left(\mathcal{R}(A^\sharp )^{rs} (M^\dagger )^{\sigma _L\tau _L}\cdots (M^\dagger )^{\sigma _1\tau _1}\right). \label{eq:mpdo_reflected_chain} \end{align}

Adjunction exchanges the oriented horizontal bond pair and reverses the unmarked tail, without identifying \(A\) with \(A^\sharp \) letterwise.

Proof ▶

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\),

\begin{align} V^\dagger \widetilde M^vV & =c\, X A^vX^{-1}, & c & {\gt}0. \notag \end{align}

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 ] .

Proof ▶

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\):

\begin{align} V_X^\dagger \widetilde M^vV_X & =c_X X A^vX^{-1}, & c_X & {\gt}0, \notag \\ V_Y^\dagger \widetilde M^vV_Y & =c_Y Y A^vY^{-1}, & c_Y & {\gt}0. \notag \end{align}

Then the two Gram conjugations agree for every vertical letter,

\begin{align} (X^\dagger X)A^v(X^\dagger X)^{-1} & =(Y^\dagger Y)A^v(Y^\dagger Y)^{-1}. \label{eq:mpdo_pairwise_gram_conj} \end{align}

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 ] .

Proof ▶

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.

Definition 16.4.69 Pairwise Gram dressing for grouped corners
✓

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:

\begin{align} V_X^\dagger \widetilde M^vV_X & =c_X X A^vX^{-1}, \notag \\ V_Y^\dagger \widetilde M^vV_Y & =c_Y Y A^vY^{-1} \notag \end{align}

Then, for every vertical letter \(v\),

\begin{align} (X^\dagger X)A^v(X^\dagger X)^{-1} & =(Y^\dagger Y)A^v(Y^\dagger Y)^{-1}. \label{eq:grouped_corner_gram_dressing} \end{align}
Theorem 16.4.70 Distinguished reference corner of a grouped phase class
✓

Suppose a phase-class grouping has positive effective coefficients \(c_{\alpha ,q}\), corner identities

\begin{align} V_{\alpha ,q}^\dagger \widetilde M^vV_{\alpha ,q} & =c_{\alpha ,q}X_{\alpha ,q}A_{\alpha ,q}^vX_{\alpha ,q}^{-1}, \notag \end{align}

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\),

\begin{align} W_{\alpha ,q}^\dagger \widetilde M^vW_{\alpha ,q} & =c_{\alpha ,0}A_{\alpha ,q}^v. \label{eq:distinguished_grouped_reference_corner} \end{align}
Proof ▶

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

\begin{align} W_{\alpha ,q}^\dagger \widetilde M^vW_{\alpha ,q} & =c_{\alpha ,0}A_{\alpha ,q}^v, \notag \end{align}

yielding (47) because \(c_{\alpha ,0}{\gt}0\).

Theorem 16.4.71 Gram conjugation from grouped-corner dressing
✓

Suppose a phase-class grouping has positive effective coefficients \(c_{\alpha ,q}\) and corner identities

\begin{align} V_{\alpha ,q}^\dagger \widetilde M^vV_{\alpha ,q} & =c_{\alpha ,q}X_{\alpha ,q}A_{\alpha ,q}^v X_{\alpha ,q}^{-1}. \notag \end{align}

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\),

\begin{align} (X_{\alpha ,q}^\dagger X_{\alpha ,q})A_{\alpha ,q}^v (X_{\alpha ,q}^\dagger X_{\alpha ,q})^{-1} & =A_{\alpha ,q}^v. \label{eq:grouped_sector_gram_conj_of_dressing} \end{align}
Proof ▶

Theorem 16.4.70 transports the distinguished identity-gauge corner to the bond space of copy \(q\). The grouped-corner Gram-dressing property gives

\begin{align} (X_{\alpha ,0}^\dagger X_{\alpha ,0})A_{\alpha ,q}^v (X_{\alpha ,0}^\dagger X_{\alpha ,0})^{-1} & =(X_{\alpha ,q}^\dagger X_{\alpha ,q})A_{\alpha ,q}^v (X_{\alpha ,q}^\dagger X_{\alpha ,q})^{-1}. \end{align}

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

\begin{align} XBX^{-1} & =(X^{-1})^\dagger C X^\dagger . \notag \end{align}

Then conjugation by the Gram matrix carries \(B\) to \(C\),

\begin{align} (X^\dagger X)B(X^\dagger X)^{-1} & =C. \label{eq:mpdo_gram_target} \end{align}

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.

Proof ▶

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.

\(X_\alpha ^\dagger X_\alpha =Y_\alpha ^\dagger Y_\alpha \)
\begin{tenkz}[rows={wire}, cols=5, west=open, east=open, size=s]
            \tn[skin=ring]{X_\alpha^{\dagger}} & \tn[skin=ring]{X_\alpha} &
            \tn[label pos=s, ports={90:physical, 270:physical}]{M_\alpha} &
            \tn[skin=ring]{X_\alpha^{-1}} & \tn[skin=ring]{X_\alpha^{-1\dagger}}
        \end{tenkz} \( = \) \begin{tenkz}[rows={wire}, cols=5, west=open, east=open, size=s]
            \tn[skin=ring]{Y_\alpha^{\dagger}} & \tn[skin=ring]{Y_\alpha} &
            \tn[label pos=s, ports={90:physical, 270:physical}]{M_\alpha} &
            \tn[skin=ring]{Y_\alpha^{-1}} & \tn[skin=ring]{Y_\alpha^{-1\dagger}}
        \end{tenkz}
Figure 16.5 Equality of the Gram conjugations for two copies of the same vertical normal tensor: \((X_{\alpha ,k}^{\dagger }X_{\alpha ,k})M_\alpha (X_{\alpha ,k}^{\dagger }X_{\alpha ,k})^{-1} =(X_{\alpha ,0}^{\dagger }X_{\alpha ,0})M_\alpha (X_{\alpha ,0}^{\dagger }X_{\alpha ,0})^{-1}\). Theorem 16.4.68 gives this equality from self-adjointness of the corresponding sector compressions and horizontal marked separation. Normality is used only in the subsequent proportionality argument. This is [ 13 , Proposition 4.13, lines 1915–1921 ] .

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

\begin{align} X^\dagger X & =\omega \, Y^\dagger Y, \label{eq:mpdo_gram_proportionality} \end{align}

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.

Proof ▶

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.

Theorem 16.4.74 Positive scalar Gram matrix and unitary normalization from grouped-corner dressing
✓

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

\begin{align} X_{\alpha ,q}^\dagger X_{\alpha ,q} & =\omega _{\alpha ,q}\mathbb {1}. \label{eq:grouped_sector_gram_scalar_of_dressing} \end{align}

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.

Proof ▶

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.

Definition 16.4.75 Normalized grouped-sector map
✓
#

After identifying the bond space of a copy with that of its representative, define its normalized physical map by

\begin{align} W_{\alpha ,q} & =V_{\alpha ,q}\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}. \notag \end{align}

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

\begin{align} X_{\alpha ,q}^\dagger X_{\alpha ,q} & =\omega _{\alpha ,q}\mathbb {1}. \notag \end{align}

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\),

\begin{align} W_{\alpha ,q}^\dagger W_{\alpha ,q} & =\mathbb {1}. \end{align}

For distinct pairs \((\alpha ,q)\ne (\beta ,p)\),

\begin{align} W_{\alpha ,q}^\dagger W_{\beta ,p} & =0. \end{align}

For every vertical letter \(v\), these maps also satisfy

\begin{align} \widetilde M_vW_{\alpha ,q} & =W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v), \notag \\ \widetilde M_v & =\sum _\alpha \sum _q W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v)W_{\alpha ,q}^\dagger . \label{eq:normalized_grouped_sector_maps_of_gram} \end{align}
Proof ▶

After identifying the copy bond space with the representative bond space, define

\begin{align} Q_{\alpha ,q} & =\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}, \notag \\ W_{\alpha ,q} & =V_{\alpha ,q}Q_{\alpha ,q}. \notag \end{align}

Since \(\omega _{\alpha ,q}{\gt}0\), the Gram identity implies

\begin{align} X_{\alpha ,q}^\dagger & =\omega _{\alpha ,q}X_{\alpha ,q}^{-1}, \notag \\ Q_{\alpha ,q}^\dagger Q_{\alpha ,q} & =\mathbb {1}. \notag \end{align}

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

\begin{align} \widetilde M_vW_{\alpha ,q} & =W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v). \notag \end{align}

The same Gram identity and the positive square-root normalization give the corner equality

\begin{align} Q_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v)Q_{\alpha ,q}^\dagger & =c_{\alpha ,q}X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}. \label{eq:normalized_grouped_sector_corner} \end{align}

Multiplying (56) on the left by \(V_{\alpha ,q}\) and on the right by \(V_{\alpha ,q}^\dagger \) gives

\begin{align} W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v)W_{\alpha ,q}^\dagger & =V_{\alpha ,q} (c_{\alpha ,q}X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}) V_{\alpha ,q}^\dagger . \end{align}

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

\begin{align} Q_{\alpha ,q} & =\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}, \notag \\ W_{\alpha ,q} & =V_{\alpha ,q}Q_{\alpha ,q}, \notag \end{align}

after identifying the bond space of each copy with that of its chosen representative. Then

\begin{align} W_{\alpha ,q}^\dagger W_{\alpha ,q} & =\mathbb {1}, \notag \\ W_{\alpha ,q}^\dagger W_{\beta ,p} & =0 \quad ((\alpha ,q)\ne (\beta ,p)), \notag \\ \widetilde M_vW_{\alpha ,q} & =W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v), \notag \end{align}

and every vertical letter is reconstructed exactly:

\begin{align} \widetilde M_v & =\sum _\alpha \sum _q W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v) W_{\alpha ,q}^\dagger . \label{eq:normalized_grouped_sector_maps_of_dressing} \end{align}

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.

Proof ▶

Theorem 16.4.74 supplies the positive Gram scalars. Apply Theorem 16.4.76 to the phase-class grouping. Its corner identity preserves the grouped reconstruction and gives (58).

Lemma 16.4.78 Matrix product vector of the diagonal tensor
✓
#

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))\).

Proof ▶

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.

Lemma 16.4.79 Diagonal tensor as a block decomposition under horizontal form
✓

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\),

\begin{align} V^{(N)}(i\mapsto M^{ii})(\sigma ) & =V^{(N)}\Bigl(\textstyle \bigoplus _k \mu _kA_k^{\mathrm{diag}}\Bigr)(\sigma ). \notag \end{align}
Proof ▶

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.

Theorem 16.4.80 Per-copy horizontal hypotheses do not control diagonal normality
✓

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

\begin{align} A^{(a,b)} & =c_aE_{ab}, & c_0 & =\frac35, & c_1 & =\frac45. \notag \end{align}

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

\begin{align} A^{\mathrm{diag}}(0) & =\frac35E_{00}, & A^{\mathrm{diag}}(1) & =\frac45E_{11}. \notag \end{align}

Every product of these two matrices is diagonal, so no blocking length spans \(M_{2}(\mathbb {C})\).

Proof ▶

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.

Lemma 16.4.81 Diagonal tensor word evaluation
✓
#

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.

Proof ▶

By induction on the word.

Lemma 16.4.82 Diagonal tensor evaluates the density-operator diagonal
✓

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:

\begin{align} V^{(N)}(i\mapsto M^{ii})(\sigma ) & =\langle \sigma |\rho ^{(N)}(M)|\sigma \rangle . \label{eq:mpdo_diagonal_density} \end{align}

Both sides of (59) equal \(\operatorname{tr}(M^{\sigma _1\sigma _1}\cdots M^{\sigma _N\sigma _N})\).

Proof ▶

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.

Lemma 16.4.83 Nonnegativity of the diagonal matrix product vector
✓

If \(\rho ^{(N)}(M)\) is positive semidefinite, then \(V^{(N)}(i\mapsto M^{ii})(\sigma )\geq 0\) for every configuration \(\sigma \).

Proof ▶

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.

Lemma 16.4.84 Matrix product vector of the vertical direct-sum tensor
✓

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

\begin{align} V^{(N)}\Bigl(\textstyle \bigoplus _\alpha \bigoplus _j \omega _{\alpha j}A_\alpha \Bigr)(\sigma ) & =\sum _\alpha \sum _j(\omega _{\alpha j})^N V^{(N)}(A_\alpha )(\sigma ). \label{eq:mpdo_vertical_mpv_sum} \end{align}
Proof ▶

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),

\begin{align} V^{(N)}\Bigl(\textstyle \bigoplus _q \widetilde\omega _q\widetilde A_q\Bigr)(\sigma ) & =\sum _q(\widetilde\omega _q)^N V^{(N)}(\widetilde A_q)(\sigma ). \notag \end{align}

The flattening identifies \(q\) with the pair \((\alpha ,j)\) via the bijection

\begin{align} \{ (\alpha ,j):0\leq \alpha {\lt}g,\ 0\leq j{\lt}\operatorname{mult}_\alpha \} & \simeq \{ 0,\ldots ,\textstyle \sum _\alpha \operatorname{mult}_\alpha -1\} , \notag \end{align}

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)\).

Lemma 16.4.85 Regrouping a block family by repeated representatives
✓

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

\begin{align} e:I & \simeq \{ (\alpha ,j):\alpha \in J,\ 0\leq j{\lt}m_\alpha \} \notag \end{align}

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

\begin{align} V^{(N)}\left(\textstyle \bigoplus _{k\in I}\nu _kB_k\right) & = V^{(N)}\Bigl( \textstyle \bigoplus _{\alpha \in J} \bigoplus _{j=0}^{m_\alpha -1}\omega _{\alpha j}A_\alpha \Bigr) \notag \end{align}

for every positive length \(N\). This is the finite-sum reindexing after the vertical sectors have been grouped in [ 13 , Proposition 4.13 ] .

Proof ▶

For \(N{\gt}0\) and every configuration \(\sigma \), the two sides expand as

\begin{align} \sum _{k\in I}\nu _k^NV^{(N)}(B_k)(\sigma ), \notag \\ \sum _{\alpha \in J}\sum _{j=0}^{m_\alpha -1} \omega _{\alpha j}^NV^{(N)}(A_\alpha )(\sigma ). \notag \end{align}

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

\begin{align} \widetilde M_vP & =P\widetilde M_vP \label{eq:cpsv_literal_vertical_projector_closure} \end{align}

for every \(v\). Thus the vertically viewed tensor has invariant-projector closure.

Proof ▶

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 \),

\begin{align} P_1H^{(N)} & =P_1H^{(N)}P_1=H^{(N)}P_1. \label{eq:cpsv_literal_first_site_commutation} \end{align}

The matrix coefficients of (62) give agreement of the corresponding first-site actions. In particular,

\begin{align} P_1H^{(N)} & =H^{(N)}P_1. \label{eq:cpsv_prop_4_13_first_site} \end{align}

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

\begin{align} \widetilde M_vP & =P\widetilde M_vP. \label{eq:cpsv_prop_4_13_projector_closure} \end{align}

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.

Proof ▶

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.

Theorem 16.4.88 Irreducible reducing sectors of the vertically viewed tensor
✓

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

\begin{align} V_k^\dagger V_k & =\mathbb {1}, \notag \\ \sum _kV_kV_k^\dagger & =\mathbb {1}, \notag \\ V_k^\dagger V_l & =0\quad (k\ne l). \notag \end{align}

Every \(B_k\) is irreducible, and for every vertical letter \(v\),

\begin{align} \widetilde M_vV_k & =V_kB_k^v, \notag \\ V_k^\dagger \widetilde M_v & =B_k^vV_k^\dagger , \notag \\ B_k^v & =V_k^\dagger \widetilde M_vV_k. \notag \end{align}

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.

Proof ▶

Apply the irreducible-corner decomposition under invariant-projector closure to \(\widetilde M\).

Theorem 16.4.89 Complete vertical reducing projectors
✓

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:

\begin{align} \sum _kP_k & =\mathbb {1}, \notag \\ P_kP_l & =0\quad (k\ne l), \notag \\ \widetilde M_vP_k & =P_k\widetilde M_v. \notag \end{align}
Proof ▶

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.

\(U\mathcal K U^\dagger =\bigoplus _\alpha M_\alpha \).
Figure 16.6 The vertically viewed tensor decomposes through complete orthogonal physical sectors. Each \(V_k\) embeds the irreducible tensor \(B_k\), while \(P_k=V_kV_k^\dagger \) is its reducing range projection. The decomposition \(U\mathcal K U^\dagger =\bigoplus _\alpha M_\alpha \) is equivalent, for every vertical letter \(v\), to \(\widetilde M_v=\sum _kV_kB_{k,v}V_k^\dagger \). This is the iterated invariant-subspace construction in [ 13 , lines 200–219 ] , applied to the reduction in [ 13 , Proposition 4.13, lines 1873–1887 ] .

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\),

\begin{align} \widetilde M_vV_k & =V_k(\nu _kB_k^v), \notag \\ V_k^\dagger \widetilde M_v & =(\nu _kB_k^v)V_k^\dagger , \notag \\ \nu _kB_k^v & =V_k^\dagger \widetilde M_vV_k. \notag \end{align}

Every vertical letter is reconstructed by

\begin{align} \widetilde M_v & =\sum _kV_k(\nu _kB_k^v)V_k^\dagger . \label{eq:mpdo_vertical_normal_reconstruction} \end{align}
Proof ▶

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\),

\begin{align} \widetilde M_vV_k & =V_k(\nu _kB_k^v), \notag \\ V_k^\dagger \widetilde M_v & =(\nu _kB_k^v)V_k^\dagger , \notag \\ \nu _kB_k^v & =V_k^\dagger \widetilde M_vV_k. \notag \end{align}

Every vertical letter is reconstructed by

\begin{align} \widetilde M_v & =\sum _kV_k(\nu _kB_k^v)V_k^\dagger . \label{eq:mpdo_literal_normal_vertical_reconstruction} \end{align}

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 ] .

Proof ▶

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

\begin{align} \widetilde M_v & =\sum _kV_k(\nu _kB_k^v)V_k^\dagger . \label{eq:cpsv_prop_4_13_normal_corners} \end{align}
Definition 16.4.92 Grouped vertical BNT decomposition with physical isometries
✓

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

\begin{align} \widetilde M_vV_k & =V_k(\mu _kB_k^v), \notag \\ V_k^\dagger \widetilde M_v & =(\mu _kB_k^v)V_k^\dagger , \notag \\ \mu _kB_k^v & =V_k^\dagger \widetilde M_vV_k, \notag \\ \widetilde M_v & =\sum _kV_k(\mu _kB_k^v)V_k^\dagger . \notag \end{align}

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

\begin{align} B_{\alpha ,q}^v & =\zeta _{\alpha ,q}X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}. \notag \end{align}

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:

\begin{align} \widetilde M_v & =\sum _\alpha \sum _q V_{\alpha ,q} \bigl(c_{\alpha ,q}X_{\alpha ,q}M_\alpha ^v X_{\alpha ,q}^{-1}\bigr)V_{\alpha ,q}^\dagger . \label{eq:vertical_bnt_grouping_with_isometry} \end{align}

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

\begin{align} B_{\alpha ,q}^v & =\zeta _{\alpha ,q}X_{\alpha ,q} M_\alpha ^vX_{\alpha ,q}^{-1}. \label{eq:mpdo_vertical_gauge_phase} \end{align}

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

\begin{align} \widetilde M_v & =\sum _\alpha \sum _q V_{\alpha ,q} \bigl(\nu _{\alpha ,q}\zeta _{\alpha ,q} X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}\bigr) V_{\alpha ,q}^\dagger . \label{eq:mpdo_vertical_grouped_reconstruction} \end{align}

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.

Proof ▶

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

\begin{align} B_{\alpha ,q}^v & =\zeta _{\alpha ,q}X_{\alpha ,q} M_\alpha ^vX_{\alpha ,q}^{-1}. \label{eq:cpsv_literal_vertical_gauge_phase} \end{align}

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

\begin{align} \widetilde M_v & =\sum _\alpha \sum _q V_{\alpha ,q} \bigl(\nu _{\alpha ,q}\zeta _{\alpha ,q} X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}\bigr) V_{\alpha ,q}^\dagger . \label{eq:cpsv_literal_vertical_grouped_reconstruction} \end{align}

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.

Proof ▶

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.

Definition 16.4.95 Block-row map of sector isometries
✓
#

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.

Theorem 16.4.96 Conjugation and reconstruction from orthogonal sectors
✓

Suppose that

\begin{align} W_k^\dagger W_k & =\mathbb {1}, \notag \\ W_k^\dagger W_l & =0\quad (k\ne l), \notag \end{align}

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:

\begin{align} UU^\dagger & =\mathbb {1}, \label{eq:mpdo_u_coisometry} \\ UT_vU^\dagger & =\bigoplus _kB_k^v. \label{eq:mpdo_u_conj_direct_sum} \end{align}

If the letters also satisfy the literal sector decomposition

\begin{align} T_v & =\sum _kW_kB_k^vW_k^\dagger , \label{eq:mpdo_sector_literal_decomposition} \end{align}

then

\begin{align} T_v & =U^\dagger \! \left(\bigoplus _kB_k^v\right)U. \notag \end{align}

This last equality is a conditional consequence of (75); no additional assertion from the source is used.

Proof ▶

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.

Theorem 16.4.97 Vertical canonical form from orthogonal grouped sectors
✓

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

\begin{align} W_{\alpha ,q}^\dagger W_{\alpha ,q} & =\mathbb {1}, \notag \\ W_{\alpha ,q}^\dagger W_{\beta ,p} & =0 \quad ((\alpha ,q)\ne (\beta ,p)). \notag \end{align}

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

\begin{align} \widetilde M_v & =\sum _{\alpha =0}^{g-1}\sum _{q=0}^{r_\alpha -1} W_{\alpha ,q}(\omega _{\alpha ,q}A_\alpha ^v) W_{\alpha ,q}^\dagger . \notag \end{align}

Then \(M\) is in vertical canonical form. This is the final finite-dimensional implication in [ 13 , Proposition 4.13, lines 1863–1870 ] .

Proof ▶

Form the block-row map \(U\) from the adjoints \(W_{\alpha ,q}^\dagger \). Theorem 16.4.96 gives

\begin{align} UU^\dagger & =\mathbb {1}, \notag \\ U\widetilde M_vU^\dagger & =\bigoplus _{\alpha ,q}\omega _{\alpha ,q}A_\alpha ^v, \notag \\ \widetilde M_v & =U^\dagger \! \left(\bigoplus _{\alpha ,q}\omega _{\alpha ,q}A_\alpha ^v\right)U. \notag \end{align}

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.

Lemma 16.4.98 A retained vertical phase class exists
✓

Under the hypotheses of Theorem 16.4.93, the retained vertical decomposition has at least one sector and hence at least one phase class.

Proof ▶

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 ] .

Proof ▶

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

\begin{align} V_k^\dagger T_v & =(\mu _kB_k^v)V_k^\dagger , \notag \\ T_v & =\sum _{k=0}^{r-1}V_k(\mu _kB_k^v)V_k^\dagger . \notag \end{align}

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.

Proof ▶

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

\begin{align} \operatorname{tr}(T^w) & =\sum _k\operatorname{tr}\! \left(V_kC_k^vV_k^\dagger T^{w'}\right) \notag \\ & =\sum _k\operatorname{tr}(C_k^vC_k^{w'}) =\sum _k\mu _k^{|w|}\operatorname{tr}(B_k^w) =V^{(|w|)}(A_{\mathrm{ret}})_w. \notag \end{align}

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

\begin{align} V_X^\dagger \widetilde M^vV_X & =c_X X A^vX^{-1}, & c_X & {\gt}0, \notag \\ V_Y^\dagger \widetilde M^vV_Y & =c_Y Y A^vY^{-1}, & c_Y & {\gt}0. \notag \end{align}

Then the Gram-dressed tensors agree:

\begin{align} (X^\dagger X)A^v(X^\dagger X)^{-1} & =(Y^\dagger Y)A^v(Y^\dagger Y)^{-1} \label{eq:cpsv_literal_pairwise_gram_conj} \end{align}

for every vertical letter \(v\).

Proof ▶

For \(Z=X,Y\), set

\begin{align} f^Z_{(r,s),(i,j)} & =c_Z^{-1}\overline{(V_Z Z)_{ir}} \bigl(V_Z(Z^{-1})^\dagger \bigr)_{js}. \notag \end{align}

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

\begin{align} V_X^\dagger \widetilde M^vV_X & =c_X X A^vX^{-1}, & c_X & {\gt}0, \notag \\ V_Y^\dagger \widetilde M^vV_Y & =c_Y Y A^vY^{-1}, & c_Y & {\gt}0. \notag \end{align}

Then there is a real number \(\omega {\gt}0\) such that

\begin{align} X^\dagger X & =\omega \, Y^\dagger Y. \label{eq:cpsv_literal_positive_relative_gram} \end{align}
Proof ▶

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.

If \(M\) has a grouped vertical BNT decomposition with physical isometries in the sense of Definition 16.4.92, and \(M\) has the grouped-corner Gram-dressing property of Definition 16.4.69, then \(M\) is in vertical canonical form.

Proof ▶

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

\begin{align} UU^\dagger & =\mathbb {1}. \label{eq:cpsv_prop_4_13_coisometry} \end{align}
\begin{align} U\widetilde M U^\dagger & =B. \label{eq:cpsv_prop_4_13_vertical_form} \end{align}
\begin{align} \widetilde M & =U^\dagger BU. \label{eq:cpsv_prop_4_13_reconstruction} \end{align}

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.

Proof ▶

Theorem 16.4.94 supplies the grouped vertical BNT decomposition with physical isometries. The literal pairwise Figure 8 theorem supplies the grouped-corner Gram-dressing property. Apply Theorem 16.4.103.

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

\begin{align} UU^\dagger & =\mathbb {1}, \notag \\ U\widetilde MU^\dagger & =B, \notag \\ \widetilde M & =U^\dagger BU. \notag \end{align}

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.

Proof ▶

Theorem 16.4.93 supplies the grouped vertical decomposition. The pairwise Figure 8 theorem supplies its grouped-corner Gram-dressing property. Apply Theorem 16.4.103.

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

\begin{align} \rho ^{(L)}(\mathcal X,M)_{\sigma ,\tau } & =\operatorname{tr}\bigl(\mathcal X M^{\sigma _1\tau _1}\cdots M^{\sigma _L\tau _L}\bigr). \label{eq:mpdo_boundary_operator} \end{align}

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 ] .

Theorem 16.5.2 Commutation of boundaries with virtual words
✓

An admissible boundary commutes with every product of virtual letters.

Proof ▶

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.

Definition 16.5.4 On-site symmetry
✓

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 ] .

Theorem 16.5.5 Periodic operators of an on-site tensor
✓

The periodic operator of length \(L\) of the bond dimension one tensor with entries \(u_{ij}\) is \(u^{\otimes L}\).

Proof ▶

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}\).

\begin{tenkz}[rows={op,op}, cols=4, west=trace, east=trace, bonds=none]
        \tn[skin=box, ports={90:physical:$i_1$, 270:physical}]{U_g} &
        \tn[skin=box, ports={90:physical:$i_2$, 270:physical}]{U_g} &
        \tn[skin=dots]{} &
        \tn[skin=box, ports={90:physical:$i_L$, 270:physical}]{U_g} \\
        \tn[skin=mpo, ports={90:physical, 270:physical:$j_1$}]{M} &
        \tn[skin=mpo, ports={90:physical, 270:physical:$j_2$}]{M} &
        \tn[skin=dots]{} &
        \tn[skin=mpo, ports={90:physical, 270:physical:$j_L$}]{M}
        \tnwire{(1,1)}{(1,2)} \tnwire{(1,2)}{(1,3)} \tnwire{(1,3)}{(1,4)}
        \tnwire{(2,1)}{(2,2)} \tnwire{(2,2)}{(2,3)} \tnwire{(2,3)}{(2,4)}
        \tnwire{(1,1)}{(2,1)} \tnwire{(1,2)}{(2,2)} \tnwire{(1,4)}{(2,4)}
    \end{tenkz}  \(=\lambda _g^{(L)}\) \begin{tenkz}[rows={op}, cols=4, boundary=periodic]
        \tn[skin=mpo, ports={90:physical:$i_1$, 270:physical:$j_1$}]{M} &
        \tn[skin=mpo, ports={90:physical:$i_2$, 270:physical:$j_2$}]{M} &
        \tn[skin=dots, ports={180:virtual, 0:virtual}]{} &
        \tn[skin=mpo, ports={90:physical:$i_L$, 270:physical:$j_L$}]{M}
    \end{tenkz}

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 ] .

Proof ▶

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.

Proof ▶

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 ] .

Proof ▶

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.

Theorem 16.5.9 Local conditions for strong and weak symmetry
✓

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.

Proof ▶

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.

\begin{tenkz}[rows={op,op}, cols=1]
        \tn[skin=box, ports={90:physical:$i$, 270:physical}]{U_g} \\
        \tn[skin=mpo,
            ports={180:virtual, 0:virtual, 90:physical, 270:physical:$k$}]{A}
    \end{tenkz}  \(=\zeta _g\) \begin{tenkz}[cols=3]
        \tn[skin=ring, ports={180:virtual, 0:virtual}]{X_g} &
        \tn[skin=mpo,
            ports={180:virtual, 0:virtual, 90:physical:$i$, 270:physical:$k$}]{A} &
        \tn[skin=ring, ports={180:virtual, 0:virtual}]{X_g^{-1}}
    \end{tenkz}     \begin{tenkz}[rows={op,op,op}, cols=1]
        \tn[skin=box, ports={90:physical:$i$, 270:physical}]{U_g} \\
        \tn[skin=mpo,
            ports={180:virtual, 0:virtual, 90:physical, 270:physical}]{A} \\
        \tn[skin=box, ports={90:physical, 270:physical:$k$}]{V_g}
    \end{tenkz}  \(=\zeta _g\) \begin{tenkz}[cols=3]
        \tn[skin=ring, ports={180:virtual, 0:virtual}]{X_g} &
        \tn[skin=mpo,
            ports={180:virtual, 0:virtual, 90:physical:$i$, 270:physical:$k$}]{A} &
        \tn[skin=ring, ports={180:virtual, 0:virtual}]{X_g^{-1}}
    \end{tenkz}

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 ] .

Proof ▶

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

\begin{align} |V^{(N)}(B_g)\rangle & =c_g(N)|V^{(N)}(A)\rangle \notag \end{align}

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.

Proof ▶

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.

Proof ▶

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.

Proof ▶

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.

Proof ▶

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.