Tensor Network Theory: A formalization blueprint

23 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. The renormalization fixed-point conditions and their algebraic consequences are developed in Chapters 26 and 27.

23.1 MPO tensors

Definition 23.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}[physical=updown]
            \tn[mpo, up=$i$, down=$j$]{M}
        \end{tenkz}
Definition 23.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 23.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 23.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 23.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 23.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 23.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}[physical=updown, periodic]
        \tn[mpo, up=$\sigma_1$, down=$\tau_1$]{M} &
        \tn[mpo, up=$\sigma_2$, down=$\tau_2$]{M} &
        \tndots &
        \tn[mpo, up=$\sigma_N$, down=$\tau_N$]{M}
    \end{tenkz}
Definition 23.1.7 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 23.1.8 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 23.1.9 Hermitian MPO tensor
#

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

Definition 23.1.10 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.

23.2 Transfer maps

Definition 23.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}
Lemma 23.2.2 Identification with the doubled-index MPS transfer map

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 23.2.3 Positivity of the MPO transfer map

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

Proof

By Lemma 23.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\).

23.3 MPDOs and LPDOs

Definition 23.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 [ CPGSV16 , Section 4 ] .

Definition 23.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 [ CPGSV16 , Section 4.3 ] . The local contraction is

\begin{tenkz}[physical=updown, tensor style=box]
            \tn[mpo, up=$i$, down=$j$]{M}
        \end{tenkz} \( = \) \begin{tenkz}[rows={op,op}, tensor style=box]
            \tnfuse{e} & \tn[up=$i$]{A} &
                \tnfuse[combined=east]{e^{-1}}\\
                       & \tn*[down=$j$]{A} &
        \end{tenkz}
Theorem 23.3.3 LPDOs are MPDOs

Every LPDO tensor is an MPDO.

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

23.4 Horizontal and Vertical Canonical Forms

Definition 23.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 a strengthened representative-grouped form of the canonical decomposition in [ CPGSV16 , Section 2.3 ] . In contrast with Definition 9.6.4, the retained direct sum has the full bond dimension, every copy weight is nonzero and normalized, and gauge-equivalent normal blocks are grouped over a BNT representative with a separate gauge for each copy. It is therefore not the literal CPSV canonical-form hypothesis.

Definition 23.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 [ CPGSV16 , Section 2.3, lines 298–301 ] .

Theorem 23.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}[physical=updown]
            \tn{K}
        \end{tenkz} \( = \) \begin{tenkz}[physical=updown]
            \tnX{X} & \tn{\bigoplus_j\mathcal K_j} & \tnX{X^{-1}}
        \end{tenkz}
Figure 23.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 [ CPGSV16 , Appendix C.2, lines 1660–1665 ] .
Definition 23.4.4 Vertically viewed tensor of an MPO

The vertical direction of an MPO tensor exchanges the notion of physical and virtual indices [ CPGSV16 , 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 23.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 [ CPGSV16 , 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 [ CPGSV16 , 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:none, op, op:none}, tensor style=box]
        \tn[pill, wide=3, legs at={1,2,3}]{U} & & \\
        \tn[wide=3]{\widetilde M} & & \\
        \tn[pill, wide=3, legs at={1,2,3}]{U^\dagger} & &
    \end{tenkz} \( = \bigoplus _\alpha \) \begin{tenkzfree}
        \tnput[box]{mu}{(10mm,7mm)}{\mu_\alpha}
        \tnput[box, ports={west:virtual, east:virtual, north:physical,
                    south:physical, 158:physical}]{M}{(34mm,0mm)}{M_\alpha}
        \tnput[dot]{d1}{(0mm,7mm)}{}
        \tnput[dot]{d2}{(20mm,7mm)}{}
        \tnput[dot, ports={east:physical}]{d3}{(d2 |- M.158)}{}
        \tnjoin{$(d1)+(0mm,5mm)$}{$(d1)+(0mm,-19mm)$}
        \tnjoin{$(d2)+(0mm,5mm)$}{$(d2)+(0mm,-19mm)$}
        \tnjoin{$(M.west)+(-10mm,0mm)$}{M.west}
        \tnjoin{M.east}{$(M.east)+(10mm,0mm)$}
        \tnjoin{M.north}{$(M.north)+(0mm,7mm)$}
        \tnjoin{M.south}{$(M.south)+(0mm,-7mm)$}
        \tnjoin{d1}{mu.west}
        \tnjoin{d2}{mu.east}
        \tnjoin{d3.east}{M.158}
    \end{tenkzfree}
Figure 23.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 [ CPGSV16 , Proposition 4.13, lines 1864–1870 ] .
Remark 23.4.7 Finite witnesses for block-injective canonical form
#

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}

Spanning at length \(L\) gives separation at the same length: a trace functional that vanishes on every tuple \((A_k^w)_k\) vanishes on their span, hence on the whole product algebra, so nondegeneracy of the product trace pairing forces every \(\Delta _k=0\). Combined with injectivity, left-canonicality, and nonzero block weights, this finite-length separation supplies the older per-copy trace-separation hypotheses for the block-diagonal tensor.

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 already implies the spanning condition, hence the separation property, and hence the older per-copy trace-separation hypotheses directly.

The duplicate-scalar counterexample shows that the separation property is strictly stronger than blockwise injectivity, left-canonicality, and nonzero weights alone. Take \(d=1\) and two blocks of bond dimension \(1\) with \(A_0^{0}=A_1^{0}=1\) and weights \(\mu _0=1\), \(\mu _1=2\): both blocks are injective and left-canonical, \((A_k^0)^\dagger A_k^0=1\), and the weights are nonzero, but at every length \(L\) there is a unique word \(w\) and \(A_0^w=A_1^w=1\), so \(\Delta _0=1\), \(\Delta _1=-1\) satisfies \(\sum _k \operatorname{tr}(\Delta _k A_k^w)=0\) while \(\Delta \neq 0\). Thus the separation property fails, and no length makes the entrywise scalar family linearly independent.

The block-injectivity proposition of [ CPGSV16 , 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 separation property and the older per-copy trace-separation hypotheses.

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{tenkzfree}[tensor style=box, compact]
        \tnput[box, ports={west:virtual,east:virtual,
                    north:physical,south:physical}]
        {biInvOne}{(0,8mm)}{K^{-1}}
        \tnput[box, ports={west:virtual,east:virtual,north:physical}]
        {biKOne}{(0,-5mm)}{K}
        \tnput[boundary, ports={east:virtual}]
        {biInvOneWest}{(-8mm,8mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biInvOneEast}{(8mm,8mm)}{}
        \tnput[boundary, ports={east:virtual}]
        {biKOneWest}{(-8mm,-5mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biKOneEast}{(8mm,-5mm)}{}
        \tnput[boundary, ports={south:physical}]
        {biOutOne}{(0,18mm)}{}
        \tnjoin{biInvOneWest.east}{biInvOne.west}
        \tnjoin{biInvOne.east}{biInvOneEast.west}
        \tnjoin{biKOneWest.east}{biKOne.west}
        \tnjoin{biKOne.east}{biKOneEast.west}
        \tnjoin{biOutOne.south}{biInvOne.north}
        \tnjoin{biInvOne.south}{biKOne.north}
    \end{tenkzfree} \( = \) \begin{tenkzfree}[tensor style=box, compact]
        \tnput[box, ports={west:virtual,east:virtual,
                    north:physical,south:physical}]
        {biInvTwo}{(18mm,15mm)}{K^{-1}}
        \tnput[box, ports={west:virtual,east:virtual,
                    south west:virtual,south east:virtual}]
        {biXTwo}{(0,0)}{X}
        \tnput[box, ports={west:virtual,east:virtual,
                    north:physical,south:virtual}]
        {biKTwo}{(18mm,0)}{\mathcal K_j}
        \tnput[box, ports={west:virtual,east:virtual,
                    south west:virtual,south east:virtual}]
        {biXinvTwo}{(36mm,0)}{X^{-1}}
        \tnput[boundary, ports={east:virtual}]
        {biInvTwoWest}{(10mm,15mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biInvTwoEast}{(26mm,15mm)}{}
        \tnput[boundary, ports={east:virtual}]
        {biGaugeTwoWest}{(-7mm,0)}{}
        \tnput[boundary, ports={west:virtual}]
        {biGaugeTwoEast}{(43mm,0)}{}
        \tnput[boundary, ports={south:physical}]
        {biOutTwo}{(18mm,27mm)}{}
        \tnput[boundary, label pos=west, ports={east:virtual}]
        {biJTwoWest}{(-7mm,-10mm)}{j}
        \tnput[dot, ports={west:virtual,east:virtual,
                    north:virtual,south:virtual}]
        {biXjTwo}{(biXTwo.south west |- 0,-10mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,north:virtual}]
        {biKjTwo}{(18mm,-10mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,
                    north:virtual,south:virtual}]
        {biXinvjTwo}{(biXinvTwo.south west |- 0,-10mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biJTwoEast}{(43mm,-10mm)}{}
        \tnput[boundary, label pos=west, ports={east:virtual}]
        {biQTwoWest}{(-7mm,-18mm)}{q}
        \tnput[dot, ports={west:virtual,east:virtual,north:virtual}]
        {biXqTwo}{(biXTwo.south east |- 0,-18mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,north:virtual}]
        {biXinvqTwo}{(biXinvTwo.south east |- 0,-18mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biQTwoEast}{(43mm,-18mm)}{}
        \tnjoin{biInvTwoWest.east}{biInvTwo.west}
        \tnjoin{biInvTwo.east}{biInvTwoEast.west}
        \tnjoin{biOutTwo.south}{biInvTwo.north}
        \tnjoin{biInvTwo.south}{biKTwo.north}
        \tnjoin{biGaugeTwoWest.east}{biXTwo.west}
        \tnjoin{biXTwo.east}{biKTwo.west}
        \tnjoin{biKTwo.east}{biXinvTwo.west}
        \tnjoin{biXinvTwo.east}{biGaugeTwoEast.west}
        \tnjoin{biXTwo.south west}{biXjTwo.north}
        \tnjoin{biXTwo.south east}{biXqTwo.north}
        \tnjoin{biKTwo.south}{biKjTwo.north}
        \tnjoin{biXinvTwo.south west}{biXinvjTwo.north}
        \tnjoin{biXinvTwo.south east}{biXinvqTwo.north}
        \tnjoin{biJTwoWest.east}{biXjTwo.west}
        \tnjoin{biXjTwo.east}{biKjTwo.west}
        \tnjoin{biKjTwo.east}{biXinvjTwo.west}
        \tnjoin{biXinvjTwo.east}{biJTwoEast.west}
        \tnjoin{biQTwoWest.east}{biXqTwo.west}
        \tnjoin{biXqTwo.east}{biXinvqTwo.west}
        \tnjoin{biXinvqTwo.east}{biQTwoEast.west}
    \end{tenkzfree} \( = \) \begin{tenkzfree}[tensor style=box, compact]
        \tnput[box, ports={west:virtual,east:virtual,
                    south west:virtual,south east:virtual}]
        {biXThree}{(0,0)}{X}
        \tnput[box, ports={west:virtual,east:virtual,
                    south west:virtual,south east:virtual}]
        {biXinvThree}{(36mm,0)}{X^{-1}}
        \tnput[boundary, ports={east:virtual}]
        {biGaugeThreeWest}{(-7mm,0)}{}
        \tnput[boundary, ports={west:virtual}]
        {biGaugeThreeEast}{(43mm,0)}{}
        \tnput[boundary, ports={east:virtual}]
        {biLeftUnit}{(9mm,10mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biRightUnit}{(27mm,10mm)}{}
        \tnput[boundary, ports={south:physical}]
        {biOutThree}{(18mm,14mm)}{}
        \tnput[boundary, label pos=west, ports={east:virtual}]
        {biJThreeWest}{(-7mm,-10mm)}{j}
        \tnput[dot, ports={west:virtual,east:virtual,
                    north:virtual,south:virtual}]
        {biXjThree}{(biXThree.south west |- 0,-10mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,north:physical}]
        {biDeltaThree}{(18mm,-10mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,
                    north:virtual,south:virtual}]
        {biXinvjThree}{(biXinvThree.south west |- 0,-10mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biJThreeEast}{(43mm,-10mm)}{}
        \tnput[boundary, label pos=west, ports={east:virtual}]
        {biQThreeWest}{(-7mm,-18mm)}{q}
        \tnput[dot, ports={west:virtual,east:virtual,north:virtual}]
        {biXqThree}{(biXThree.south east |- 0,-18mm)}{}
        \tnput[dot, ports={west:virtual,east:virtual,north:virtual}]
        {biXinvqThree}{(biXinvThree.south east |- 0,-18mm)}{}
        \tnput[boundary, ports={west:virtual}]
        {biQThreeEast}{(43mm,-18mm)}{}
        \tnjoin{biGaugeThreeWest.east}{biXThree.west}
        \tnjoin[route=arc,out=0,in=0]
        {biXThree.east}{biLeftUnit.east}
        \tnjoin[route=arc,out=180,in=180]
        {biXinvThree.west}{biRightUnit.west}
        \tnjoin{biXinvThree.east}{biGaugeThreeEast.west}
        \tnjoin{biOutThree.south}{biDeltaThree.north}
        \tnjoin{biXThree.south west}{biXjThree.north}
        \tnjoin{biXThree.south east}{biXqThree.north}
        \tnjoin{biXinvThree.south west}{biXinvjThree.north}
        \tnjoin{biXinvThree.south east}{biXinvqThree.north}
        \tnjoin{biJThreeWest.east}{biXjThree.west}
        \tnjoin{biXjThree.east}{biDeltaThree.west}
        \tnjoin{biDeltaThree.east}{biXinvjThree.west}
        \tnjoin{biXinvjThree.east}{biJThreeEast.west}
        \tnjoin{biQThreeWest.east}{biXqThree.west}
        \tnjoin{biXqThree.east}{biXinvqThree.west}
        \tnjoin{biXinvqThree.east}{biQThreeEast.west}
    \end{tenkzfree}

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

Theorem 23.4.8 Finite block-separation criterion

If a block family admits the finite-length data described in the block-injectivity proposition of [ CPGSV16 , lines 340–345 ] , then it is in block-injective canonical form.

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.

Lemma 23.4.9 Single-block word-tuple span

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 [ CPGSV16 , 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 10.6.1.17; its detailed source comparison is recorded in docs/paper-gaps/cpgsv17_bicf_block_separation.tex.

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 10.5.3 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 23.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 [ CPGSV16 , Proposition 4.13, lines 1909–1919 ] .

Lemma 23.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. Suppose that, for every marked index \(s\) and every unmarked word \(w\),

\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, one has \(C_j=E_j\) for every \(j\).

Lemma L is stated for marks induced by physical operators on the first site. After its closed-chain identity has been expanded, the separation calculation uses only the matrices in the marked slot. The theorem is the resulting algebraic extension of [ CPGSV16 , Appendix C.3, lines 1848–1858 ] ; it is not attributed to the paper as a separately stated arbitrary-mark theorem.

Proof

Fix \(s\) and a tail length \(L\), and put \(\Delta _j=c_j^{(L+1)}(C_j^s-E_j^s)\). Lemma 23.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 10.5.3 to the nonzero weights above each fixed representative, exactly as in Theorem 23.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 10.6.1.9, and let \(P_{\mathrm{act}}\) be its active 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 active 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)}\), and every inactive listed coordinate carries the zero block. Thus all original listed coordinates remain present.

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{act}}^sP_{\mathrm{act}}^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. Inactive coordinates contribute zero.

Proof

Split the finite sum over listed coordinates into active and inactive parts. The inactive marked blocks vanish. Reindex the active part by the class-copy pairs \((j,q)\) and use the bond-dimension identifications to replace each active 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 (??).

Let \((C_j^s)_j\) and \((E_j^s)_j\) be two marked families on the active representatives. Form their full-coordinate marked tensors by retaining the inactive listed coordinates as zero blocks. 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\), then \(C_j=E_j\) for every active 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 active representative \(j\) and every physical letter \(i\).

Proof

Theorem 23.4.15 rewrites (??) into marked-trace agreement for \(P_{\mathrm{act}}\). Its representatives have eventual word-tuple separation, and every copy weight is nonzero. Theorem 23.4.14 therefore gives \(C_j=E_j\).

First-site agreement for \(T_{\mathrm{grp}}\) also holds for \(P_{\mathrm{act}}\) because their matrix product vectors agree at every positive length. Applying Theorem 23.4.11 gives (??).

Let \(A\) be in literal CPSV canonical form. For canonical-form data and its active 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}\), \(G\) is the listed block gauge, and \(T_{\mathrm{grp}}\) retains every inactive listed coordinate with weight zero. 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 marked matrix makes the closed chain nonempty; no equality of unmarked empty-word matrix product vectors is asserted.

For a linear physical mark \(f\), the linear marking of \(T_{\mathrm{grp}}\) is exactly its full grouped marked tensor. Each active copy \((j,q)\) carries \(\lambda _{j,q}f(C_j)\), where \(\lambda _{j,q}=\nu _{k(j,q)}\zeta _{k(j,q)}\), while every inactive marked block is zero. 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\), then \(F_A=G_A\) on the original bond space.

Proof

Induction on \(w\) in (??) gives \(CA^w=V^\dagger EB^wV\). Cyclicity of trace and \(VV^\dagger =\mathbb {1}\) give (??). Linearity of the physical mark gives the grouped marked tensor, with inactive blocks killed by their zero weights.

Apply (??) to both sides of (??), remove the block gauge by cyclicity of trace, and use Theorem 23.4.16 on the active representatives. The resulting equality of the full grouped marked tensors, followed by (??), gives \(F_A=G_A\).

Definition 23.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 auxiliary statement supplies the identity needed to combine the block-injectivity argument [ CPGSV16 , lines 318–345 ] with Appendix C.3, Lemma L [ CPGSV16 , 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 (??) 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 [ CPGSV16 , lines 318–345 ] with Appendix C.3, Lemma L [ CPGSV16 , 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 (??) 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 verbatim 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 23.4.19 recovers \(YA_j=ZA_j\) before blocking.

Corollary 23.4.21 Marked-letter separation for a BNT canonical form

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\), 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 23.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\), then \(C^u=E^u\) for every \(u\).

The marks are required to lie in the span of the physical letters. No assertion is made for arbitrary bond-space matrices, since off-diagonal matrices between repeated sectors are invisible to closed chains. The corollary is the physical-insertion consequence of [ CPGSV16 , Appendix C.3, lines 1835–1858 ] used in the contractions of [ CPGSV16 , 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 23.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 (??) gives the opposite one-sided identity on the original bond space.

This is the representative-indexed conclusion of the invariant-projection step in [ CPGSV16 , Proposition 4.13, lines 1873–1887 ] . Copies with a common normal representative are grouped before separation, exactly as in Appendix C.3, Lemma L. The source assumes that \(Q\) is an orthogonal projector; this implication only uses its Hermiticity.

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 23.4.17 identifies the two linear marks and gives (??). 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.

Lemma 23.4.24 Blockwise equality from agreement on the MPV family

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 (per-block separation). The source’s Lemma L separates blockwise contributions at the level of a minimal basis of normal tensors, grouping gauge-equivalent copies of one representative before the trace-separation step; the block-injective hypothesis stated here instead demands trace separation simultaneously across every block, including gauge-equivalent copies, and so excludes decompositions with repeated normal tensors. Documented in docs/paper-gaps/cpgsv17_bicf_block_separation.tex.

\(YH^{(N)}(\widetilde M)=ZH^{(N)}(\widetilde M)\)

\begin{tenkz}[periodic, rows={wire, op:none, ket}, align=3]
            \tnskip & \tnskip & \tnskip & \tnskip \\
            \tn{Y} & \tnskip & \tnskip & \tnskip \\
            \tn{\widetilde M} & \tn{\widetilde M} & \tndots & \tn{\widetilde M}
        \end{tenkz}

\( = \)

\begin{tenkz}[periodic, rows={wire, op:none, ket}, align=3]
            \tnskip & \tnskip & \tnskip & \tnskip \\
            \tn{Z} & \tnskip & \tnskip & \tnskip \\
            \tn{\widetilde M} & \tn{\widetilde M} & \tndots & \tn{\widetilde M}
        \end{tenkz}

\(YA_k=ZA_k\)

\begin{tenkz}[rows={op:none, ket}]
            \tn{Y} \\
            \tn{A_k}
        \end{tenkz}

\( = \)

\begin{tenkz}[rows={op:none, ket}]
            \tn{Z} \\
            \tn{A_k}
        \end{tenkz}
Figure 23.3 Blockwise equality from agreement on the MPV family. Top: the hypothesis \(Y_1V^{(N)}(A)=Z_1V^{(N)}(A)\) for every length \(N\), 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, [ CPGSV16 , Appendix C.3, lines 1835–1846 ] .
Proof

This is the per-block-separated form of Lemma L of [ CPGSV16 , Appendix C.3, lines 1835–1846 ] , using the finite criterion from Theorem 23.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}[rows={wire, op, ket}, periodic]
         & \\
         & \tn{Y} \\
        \tn[no legs]{W} & \tn{A_k}
    \end{tenkz} \( = \) \begin{tenkz}[rows={wire, op, ket}, periodic]
         & \\
         & \tn{Z} \\
        \tn[no legs]{W} & \tn{A_k}
    \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, [ CPGSV16 , Appendix C.3, lines 1848–1858 ] .

Theorem 23.4.25 One-sided invariant corners of positive operators

Let \(H\geq 0\) and let \(P=P^\dagger \). If \(PH=PHP\), then \((\mathbb {1}-P)HP=0\). In particular, this conclusion holds for every nonempty-chain density operator generated by an MPDO. This is the operator identity used in the proof of [ CPGSV16 , 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\).

Theorem 23.4.26 First-site contraction reduction

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 is the representative-grouped application of Lemma L and anticipates the concrete instantiation with \(Y_L,Y_C,Y_R\) the physical-index insertions of \(P\), \(PHP\), \(HP\) in Theorem 23.4.30.

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 23.4.27 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}[physical=updown, periodic]
        \tn[up=$\sigma_1$, down=$\tau_1$]{\widetilde M_1} &
        \tn[up=$\sigma_2$, down=$\tau_2$]{\widetilde M_2} &
        \tndots &
        \tn[up=$\sigma_N$, down=$\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 [ CPGSV16 , Proposition 4.13, lines 1876–1880 ] .

Lemma 23.4.28 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 (??).

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}[rows={wire, op:none, op}, compact, periodic, align=3]
            \tnskip & \tnskip & \tnskip & \tnskip \\
            \tnskip & \tnskip & \tnskip & \tn{P} \\
            \tn[mpo]{\widetilde M} & \tn[mpo]{\widetilde M} & \tndots &
                \tn[mpo]{\widetilde M}
        \end{tenkz}
Figure 23.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 [ CPGSV16 , Proposition 4.13, lines 1881–1887 ] .
Theorem 23.4.30 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 23.4.28 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 23.4.26 with \((Y_L,Y_C,Y_R)\).

View the MPO tensor \(M\) vertically as the family of physical-space operators \(M_{ab}\) of Definition 23.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 [ CPGSV16 , 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 23.4.25 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 23.4.32 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 23.4.30, whose conclusion is the equality of the insertions of \(MP\) and \(PMP\) on every minimal representative.

Theorem 23.4.34 Commutation with a power of a positive semidefinite matrix

Let \(H\geq 0\), let \(p\geq 1\), and let \(Q\) satisfy \(QH^p=H^pQ\). Then \(QH=HQ\). In particular, 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 [ CPGSV16 , 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 23.4.35 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 [ CPGSV16 , Proposition 4.13 ] , lines 1889–1891, attaches to a nontrivial vertical period; the second family is the displayed identity at lines 1890–1891.

Theorem 23.4.36 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 23.4.34), contradicting \(Q_1H^{(N+1)}\ne H^{(N+1)}Q_1\). This is the contradiction concluding the periodic-sector step in the proof of [ CPGSV16 , Proposition 4.13 ] , lines 1892–1893.

Definition 23.4.37 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 9.7.3 — there is an exponent \(q\ge 1\) for which \(M\) admits a periodic-sector projector. The proof of [ CPGSV16 , 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 23.4.45); the passage from the displacement to the all-length non-commutation family is not needed 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 [ CPGSV16 , 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 23.4.36 shows that no such projector exists.

Definition 23.4.39 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.

Theorem 23.4.40 Stacked layers generate powers of the density operator

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 [ CPGSV16 , 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 23.4.31.

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 23.4.41 the word invariance gives \(QM^{[p]}=QM^{[p]}Q\) on the stacked tensor, which generates matrix product density operators by Theorem 23.4.40. Theorem 23.4.32 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 23.4.35.

Definition 23.4.43 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 9.7.3 — 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 [ Wol12 , 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 23.4.37, the commutation of \(Q\) with \([H^{(N)}]^p\) is no longer assumed; it is derived through the stacked-layers identity.

Let \(M\) generate matrix product density operators. If periodic vectors of \(M\) supply cyclic projectors, then they supply periodic-sector projectors, and the vertically viewed tensor \(\widetilde M\) has no nontrivial \(p\)-periodic vectors. This reduces the hypothesis of the periodic-sector step in the proof of [ CPGSV16 , 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 23.4.42. The exclusion of periodic vectors is then Theorem 23.4.38.

Let \(M\) be an MPO tensor and let \(\widetilde M\) carry a nontrivial \(p\)-periodic vector, presented as in Definition 9.7.3: 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 23.4.31. This realizes the projector that the proof of [ CPGSV16 , 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 E.3.5), and the eigenvalue \(\mu /r\ne 1\) of modulus one (Lemma E.3.6) forces the peripheral cyclic group to have order \(m\ge 2\). Theorem 7.9.1.1 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 7.9.1.2). 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 E.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 23.4.46 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}

Thus a displaced idempotent has a noncommuting chain length. The proof of [ CPGSV16 , Proposition 4.13, lines 1888–1893 ] only requires this existential conclusion, not noncommutation at every length.

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 23.4.23 identifies the two inserted tensors on the original bond space. By Lemma 23.4.29, 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 23.4.45). By Theorem 23.4.46, 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 23.4.34 implies \(Q_1H^{(N+1)}=H^{(N+1)}Q_1\), a contradiction. This is the periodic exclusion in [ CPGSV16 , 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.

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 23.4.23 and Lemma 23.4.29 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 23.4.45, 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 23.4.34 imply that \(Q_1\) commutes with \(H^{(N+1)}\), which is the required contradiction. This uses the pointwise contradiction in [ CPGSV16 , Proposition 4.13, lines 1888–1893 ] ; the stronger all-length inequality printed there is unnecessary.

Definition 23.4.49 Displaced projectors never commute

Say that displaced projectors of \(M\) never commute if every orthogonal projector \(Q\) on the physical space with \(Q\widetilde M\ne Q\widetilde MQ\) satisfies \(Q_1H^{(N)}\ne H^{(N)}Q_1\) at every length. This is the canonical-form step of the periodic-sector argument: the proof of [ CPGSV16 , Proposition 4.13 ] assumes the tensor is in canonical form in the horizontal direction, and the argument of lines 1874–1887 turns all-length commutation of \(Q_1\) with \(H^{(N)}\) into the letter-level invariance \(Q\widetilde M=Q\widetilde MQ\). This argument uses the normalized BNT-refined horizontal hypothesis. The all-length condition is stronger than is needed for either Theorem 23.4.47 or Theorem 23.4.48; each uses a single noncommuting length supplied by the contrapositive of Lemma L.

Let \(M\) generate matrix product density operators. If displaced projectors of \(M\) never commute, then periodic vectors of \(M\) supply cyclic projectors, and the vertically viewed tensor \(\widetilde M\) has no nontrivial \(p\)-periodic vectors. This narrows the hypothesis of the periodic-sector step in the proof of [ CPGSV16 , Proposition 4.13, lines 1888–1893 ] to the canonical-form implication of lines 1874–1887: every other ingredient of the step is proved.

Proof

A nontrivial \(p\)-periodic vector supplies the invariant, displaced projector of Theorem 23.4.45; the non-commutation hypothesis upgrades the displacement to the all-length non-commutation family, so periodic vectors supply cyclic projectors and Theorem 23.4.44 applies.

Theorem 23.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 23.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 23.4.45. 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 23.4.51 turns the resulting commutation with \(H^{(2)}\) into the letter-level invariance that the displacement forbids.

Theorem 23.4.48 gives the corresponding result for normalized BNT-refined horizontal form.

Theorem 23.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 [ CPGSV16 , Proposition 4.13 ] concludes \(\mu _{\alpha ,k}d_\alpha {\gt}0\); the identification of the trace with \(\mu _{\alpha ,k}d_\alpha \) is Theorem 23.4.62.

Proof

Compression by a Hermitian matrix preserves positive semidefiniteness, since \(PH^{(N)}P=PH^{(N)}P^\dagger \). By Theorem 23.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 [ CPGSV16 , Proposition 4.13 ] , resolved by chain length: one loop of \(M_\alpha \) followed by loops of \(M\), closed by the horizontal trace.

Lemma 23.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 23.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 23.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 [ CPGSV16 , 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 [ CPGSV16 , Proposition 4.13, line 1898 ] , obtained from Lemma L in [ CPGSV16 , 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 (??).

Theorem 23.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 [ CPGSV16 , 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 (??) is nonzero. If (??) vanished, multiplying it by \(V^\dagger \) and \(V\) and using \(V^\dagger V=1\) would contradict (??).

Theorem 23.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 [ CPGSV16 , Proposition 4.13, line 1898 ] .

The normalized BNT-refined horizontal hypothesis is stronger than the literal CPSV canonical-form hypothesis. The literal separation statement is Theorem 23.4.59; it does not by itself construct the grouped vertical sectors used here.

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 23.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 23.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 [ CPGSV16 , Proposition 4.13 ] ; the source states it at the particular length where the compression does not vanish.

Proof

By Theorem 23.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 ,1}\) and \(\mu _{\alpha ,k}\) and a common matrix \(E_\alpha \), and suppose \(\mu _{\alpha ,1}{\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 [ CPGSV16 , Proposition 4.13 ] : since w.l.o.g. \(\mu _{\alpha ,1}{\gt}0\), the positive trace forces \(d_\alpha {\gt}0\) and hence \(\mu _{\alpha ,k}{\gt}0\).

Proof

By Theorem 23.4.54 the nonzero compression has positive trace, which Theorem 23.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 ,1)\) at the same length equals \(\mu _{\alpha ,1}d_\alpha ^{(N+1)}\ne 0\); that compression is therefore also nonzero, and its trace is a positive real. Dividing by \(\mu _{\alpha ,1}{\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\).

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 [ CPGSV16 , Proposition 4.13 ] , with \(X_{\alpha ,1}=\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}

Thus the adjoint exchanges the oriented horizontal bond pair and reverses the unmarked tail. In particular, this is not a same-tail identity and does not identify \(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\) 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 pairwise conclusion of Figures 7–8 in [ CPGSV16 , 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.

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.

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.

Theorem 23.4.67 Gram conjugation for an actual grouped sector

Consider the grouped vertical decomposition furnished by Theorem 23.4.89. For a phase class \(\alpha \) and a copy \(q\), transport its representative tensor to the bond space of that copy. If \(X_{\alpha ,q}\) is the corresponding gauge, then

\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:mpdo_grouped_gram_conj} \end{align}

for every vertical letter \(v\). This is Figure 8 for the actual witnesses of the grouped decomposition.

Proof

Transport the distinguished corner of the phase class along the equality of its bond dimension with that of the copy \(q\). Its distinguished gauge is the identity. Apply Theorem 23.4.66 to this corner and the corner indexed by \(q\). The Gram conjugation belonging to the distinguished corner is the identity, which gives the stated equality.

Theorem 23.4.68 Literal Gram conjugation for an actual grouped sector

Consider the actual phase-class grouping furnished under literal CPSV canonical form. For a class \(\alpha \) and copy \(q\), identify its normal representative with the bond space of that copy using the equality of dimensions. Then

\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:cpsv_literal_grouped_gram_conj} \end{align}

for every vertical letter \(v\).

Proof

The distinguished copy has identity gauge. Use the equality of bond dimensions to regard its positive corner equation on the bond space of the copy \(q\). The literal pairwise Figure 8 theorem compares the two corners. Its distinguished Gram conjugation is the identity, which gives (??).

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 [ CPGSV16 , Proposition 4.13, lines 1909–1919 ] . The reflected marked-chain argument in the source does not identify the target separately for either sector. The same-letter choice \(C=B^\dagger \) and the further specialization \(B^\dagger =B\) are useful corollaries, but the latter is not assumed by the source.

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}[physical=updown, compact]
            \tnX{X_\alpha^{\dagger}} & \tnX{X_\alpha} & \tn{M_\alpha} &
            \tnX{X_\alpha^{-1}} & \tnX{X_\alpha^{-1\dagger}}
        \end{tenkz} \( = \) \begin{tenkz}[physical=updown, compact]
            \tnX{Y_\alpha^{\dagger}} & \tnX{Y_\alpha} & \tn{M_\alpha} &
            \tnX{Y_\alpha^{-1}} & \tnX{Y_\alpha^{-1\dagger}}
        \end{tenkz}
Figure 23.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 ,1}^{\dagger }X_{\alpha ,1})M_\alpha (X_{\alpha ,1}^{\dagger }X_{\alpha ,1})^{-1}\). Theorem 23.4.66 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 [ CPGSV16 , 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 23.4.66 gives the Figure 8 comparison in [ CPGSV16 , 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 23.4.69, the ratio \((Y^\dagger Y)^{-1}(X^\dagger X)\) commutes with every letter of \(M_\alpha \). Theorem 23.4.64 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.

In the grouped vertical decomposition of a matrix-product density operator in normalized BNT-refined horizontal form, let \(X_{\alpha ,k}\) be the gauge of a copy of the normal representative \(M_\alpha \). Fix the distinguished gauge to \(X_{\alpha ,1}=\mathbb {1}\). Then there is a positive real number \(\omega _{\alpha ,k}\) such that

\begin{align} X_{\alpha ,k}^\dagger X_{\alpha ,k} & =\omega _{\alpha ,k}\mathbb {1}, \notag \\ U_{\alpha ,k} & =\omega _{\alpha ,k}^{-1/2}X_{\alpha ,k}, \notag \end{align}

and \(U_{\alpha ,k}\) is unitary. This is the normalization in [ CPGSV16 , Proposition 4.13, lines 1903–1921 ] .

Scope restriction (BNT-refined horizontal form): the normalized BNT-refined horizontal hypothesis is stronger than the literal CPSV canonical-form hypothesis.

Proof

Transport the normal representative to the bond space of the chosen copy. Theorem 23.4.67 states that conjugation by \(X_{\alpha ,k}^\dagger X_{\alpha ,k}\) fixes every letter of the transported tensor. Apply Theorem 23.4.70 with the second gauge equal to the identity. The resulting Gram identity has a positive real scalar, and division by its positive square root gives the stated unitary.

In the actual phase-class grouping under literal CPSV canonical form, every gauge \(X_{\alpha ,q}\) has a positive real Gram scalar \(\omega _{\alpha ,q}\) such that

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

Consequently \(U_{\alpha ,q}=\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}\) is unitary. This proves the gauge-normalization step of [ CPGSV16 , Proposition 4.13, lines 1903–1921 ] and supplies the positive Gram scalars used to normalize the physical sector maps.

Proof

The representative on this bond space is normal. Apply normal Gram rigidity to (??) with the distinguished gauge equal to the identity. This gives (??) with \(\omega _{\alpha ,q}{\gt}0\). Dividing by its positive square root gives a unitary matrix.

In the grouped vertical decomposition, 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}

and identify 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 \end{align}

Writing \(c_{\alpha ,q}=\nu _{\alpha ,q}\zeta _{\alpha ,q}{\gt}0\), the normalized maps satisfy

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

and retain the literal reconstruction

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

These are the normalized sector identities in [ CPGSV16 , Proposition 4.13, lines 1895–1921 ] .

Proof

Each \(Q_{\alpha ,q}\) is unitary, so right multiplication by it preserves the isometry equation and the pairwise orthogonality of the physical sector maps. In the intertwining relation, the scalar normalization commutes with every letter and the adjacent factors \(X_{\alpha ,q}^{-1}X_{\alpha ,q}\) cancel. The Gram identity gives \(X_{\alpha ,q}^\dagger =\omega _{\alpha ,q}X_{\alpha ,q}^{-1}\); hence

\begin{align} Q_{\alpha ,q}M_\alpha ^vQ_{\alpha ,q}^\dagger & =X_{\alpha ,q}M_\alpha ^vX_{\alpha ,q}^{-1}. \notag \end{align}

Substitution into every summand of the grouped reconstruction proves the final identity. Transport along the equality of bond dimensions changes only the coordinate labels.

Theorem 23.4.74 Literal normalized physical maps of the grouped vertical sectors

In the phase-class grouping obtained from literal CPSV canonical form, let \(Q_{\alpha ,q}=\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}\) and \(W_{\alpha ,q}=V_{\alpha ,q}Q_{\alpha ,q}\), after identifying the bond space of each copy with that of its 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}

where \(c_{\alpha ,q}=\nu _{\alpha ,q}\zeta _{\alpha ,q}{\gt}0\). Moreover, 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:cpsv_literal_normalized_reconstruction} \end{align}
Proof

For each copy, put \(Q_{\alpha ,q}=\omega _{\alpha ,q}^{-1/2}X_{\alpha ,q}\) and \(W_{\alpha ,q}=V_{\alpha ,q}Q_{\alpha ,q}\). Unitarity of \(Q_{\alpha ,q}\) preserves the isometry equation and pairwise orthogonality. Substituting \(W_{\alpha ,q}=V_{\alpha ,q}Q_{\alpha ,q}\) into the intertwining identity gives

\begin{align} \widetilde M_vW_{\alpha ,q} & =V_{\alpha ,q}\bigl( c_{\alpha ,q}\omega _{\alpha ,q}^{-1/2} X_{\alpha ,q}M_\alpha ^v X_{\alpha ,q}^{-1}X_{\alpha ,q}\bigr) \notag \\ & =W_{\alpha ,q}(c_{\alpha ,q}M_\alpha ^v). \notag \end{align}

Thus the adjacent factors \(X_{\alpha ,q}^{-1}X_{\alpha ,q}\) cancel and leave the undressed representative \(M_\alpha ^v\). The Gram identity \(X_{\alpha ,q}^\dagger X_{\alpha ,q} =\omega _{\alpha ,q}\mathbb {1}\) also gives

\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}. \notag \end{align}

Hence conjugation by \(W_{\alpha ,q}\) reproduces the original grouped corner. Summing these equalities gives (??).

Lemma 23.4.75 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 23.4.76 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 23.4.75, 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 23.4.77 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 older 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 23.4.78 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 23.4.79 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 (??) 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 23.4.80 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 23.4.79, the value is a diagonal entry of \(\rho ^{(N)}(M)\), and the diagonal entries of a positive semidefinite matrix are non-negative.

Lemma 23.4.81 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.5.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 23.4.82 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 [ CPGSV16 , 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 (??) give agreement of the corresponding first-site actions. The original-coordinate marked Lemma L gives equality of their inserted tensors on the original bond space. Reading this equality as a vertical one-site action yields (??).

Let \(M\) be in normalized BNT-refined horizontal form and generate matrix product density operators. If an orthogonal projection \(P\) satisfies \(P\widetilde M_v=P\widetilde M_vP\) for every vertical letter \(v\), then it also satisfies \(\widetilde M_vP=P\widetilde M_vP\) for every vertical letter \(v\). Thus the vertically viewed tensor has invariant-projector closure.

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 23.4.85 Irreducible reducing sectors of the vertically viewed tensor

Under the hypotheses of Theorem 23.4.84, 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}

Moreover, \(\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 23.4.86 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 23.4.85, 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 23.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 displayed identity 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 [ CPGSV16 , lines 200–219 ] , applied to the reduction in [ CPGSV16 , Proposition 4.13, lines 1873–1887 ] .
Theorem 23.4.87 Spectrally normalized vertical sectors with physical isometries

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}

Moreover, every vertical letter is reconstructed literally 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 (??). 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 [ CPGSV16 , Proposition 4.13 ] .

Proof

Theorem 23.4.83 gives invariant-projector closure for \(\widetilde M\), and Theorem 23.4.47 excludes its nontrivial periodic vectors. Applying Theorem 9.7.4 then gives normal nonzero corners with positive spectral weights and orthogonal isometries, together with the stated intertwining, compression, and reconstruction formulas.

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 ,1}=\mathbb {1}\) and \(\zeta _{\alpha ,1}=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 23.4.87 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 [ CPGSV16 , 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 10.2.6 gives equality of bond dimensions and (??) for every copy. Distinct representatives are gauge-phase separated by construction, so Theorem 10.2.11 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 (??) 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 ,1}=\mathbb {1}\) and \(\zeta _{\alpha ,1}=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 23.4.88 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 [ CPGSV16 , 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 10.2.6 gives equality of bond dimensions and (??) for every copy. Distinct representatives are gauge-phase separated by construction, so Theorem 10.2.11 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 (??) 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 23.4.91 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 [ CPGSV16 , Proposition 4.13, lines 1863–1870 ] once the sector isometries have been normalized.

Theorem 23.4.92 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}, \notag \\ UT_vU^\dagger & =\bigoplus _kB_k^v. \notag \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 (??); no additional assertion from the source is used.

Proof

The \((k,l)\) block of \(UU^\dagger \) is \(W_k^\dagger W_l\), hence the first identity follows from orthonormality. The corresponding block of \(UT_vU^\dagger \) is \(W_k^\dagger T_vW_l=W_k^\dagger W_lB_l^v\), which is \(B_k^v\) for \(k=l\) and zero otherwise. 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 23.4.93 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 [ CPGSV16 , Proposition 4.13, lines 1863–1870 ] .

Proof

Form the block-row map \(U\) from the adjoints \(W_{\alpha ,q}^\dagger \). Theorem 23.4.92 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 23.4.94 A retained vertical phase class exists

Under the hypotheses of Theorem 23.4.89, 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, (??) 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.

Theorem 23.4.95 Nonvanishing of the vertical reshaping

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 [ CPGSV16 , 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.

Theorem 23.4.96 Algebraic BNT from a grouped vertical decomposition

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 [ CPGSV16 , 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 10.2.3; in particular, it supplies the required positive bond dimensions. Lemma F.1.5 then gives the same algebraic BNT for \(T\).

Theorem 23.4.97 Pairwise Figure 8 identity under literal canonical form

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 (??).

Theorem 23.4.98 Positive relative Gram scalar under literal canonical form

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 (??) 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.

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}

The source explicitly displays (??) and calls \(U\) an isometry [ CPGSV16 , Proposition 4.13 ] . In this rectangular orientation, however, \(U\) maps the physical space onto the retained sector space and satisfies (??); thus \(U\) is a coisometry and \(U^\dagger \) is the isometric inclusion. The reverse identity (??) is not displayed separately in the proposition. It records the exact reconstruction implicit in the preceding assertion that \(\widetilde M\) is itself in canonical form, including the permitted zero orthogonal complement.

Proof outline

Write \(A\) for the doubled-index MPS tensor of \(M\). If \(A=0\), take no retained sectors; the unique map onto the zero-dimensional sector space gives the conclusion. Suppose henceforth that \(A\ne 0\). Apply Theorem 10.6.1.9 to its literal canonical-form data. The active listed indices are identified with pairs \((\alpha ,k)\) consisting of a phase class and a copy, and every active copy retains its original nonzero coefficient \(\nu _{\alpha ,k}\), its bond dimension, a unit phase, and an invertible gauge \(X_{\alpha ,k}\). The chosen tensors \(C_\alpha \) form a basis of normal tensors. Adjoin the inactive zero-weight indices to these pairs and call the resulting grouped index set \(\mathcal G\). There is an explicit equivalence from the dependent coordinates of \(\mathcal G\) to the original flattened listed coordinates. If \(T_{\mathrm{grp}}\) is the grouped direct sum, reindexed along this equivalence, then there are a block-diagonal gauge \(G\) and the original coisometry \(V\) such that

\begin{align} \bigoplus _\ell \mu _\ell A_\ell ^i & =GT_{\mathrm{grp}}^iG^{-1}, \label{eq:cpsv_prop_4_13_active_bnt} \end{align}

and

\begin{align} A^i& =V^\dagger GT_{\mathrm{grp}}^iG^{-1}V. \label{eq:cpsv_prop_4_13_active_reconstruction} \end{align}
\begin{align} VV^\dagger & =\mathbb {1}. \label{eq:cpsv_prop_4_13_original_coisometry} \end{align}

The active summands of \(T_{\mathrm{grp}}\) are \(\nu _{\alpha ,k}\zeta _{\alpha ,k}C_\alpha ^i\), with \(|\zeta _{\alpha ,k}|=1\). Every inactive summand is the original listed block with coefficient zero. Thus the grouped coordinates are exact without deleting the inactive complement or introducing a selecting coisometry.

Let \(H^{(N)}\) be the positive operator generated horizontally by \(\widetilde M\). If an orthogonal projector \(P\) on the physical space satisfies \(P\widetilde M=P\widetilde MP\), and \(P_1\) acts on the first site, then self-adjointness gives

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

Write \(Y_R(P)\) and \(Y_C(P)\) for the right and central first-site contractions determined by \(P\). The active representative sector decomposition has weights \(\lambda _{\alpha ,k}=\nu _{\alpha ,k}\zeta _{\alpha ,k}\), all nonzero, and its representatives have simultaneous word-tuple separation at every sufficiently large length. It has the same matrix product vectors as \(T_{\mathrm{grp}}\) at every positive length. Hence Theorem 23.4.16 gives

\begin{align} Y_R(P)C_\alpha ^i& =Y_C(P)C_\alpha ^i \label{eq:cpsv_prop_4_13_active_separation} \end{align}

on every active representative. Theorem 23.4.17 then uses the listed block gauge and the ambient coisometry, with \(VV^\dagger =\mathbb {1}\), to give the original-space inserted-tensor equality. Theorem 23.4.83 therefore yields

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

Thus every one-sided invariant subspace is reducing. The inactive listed coordinates remain present with weight zero throughout this argument; no insertion identity is asserted for an individual inactive block.

Theorem 23.4.47 excludes nontrivial periodic vertical sectors. Its proof uses the displaced cyclic projector \(Q\) and one length \(N+1\) for which \(Q_1H^{(N+1)}\ne H^{(N+1)}Q_1\); full-period invariance gives commutation with \([H^{(N+1)}]^p\), and positivity gives the contradictory commutation with \(H^{(N+1)}\) itself.

Combining this exclusion with invariant-projector closure gives the normal corners of Theorem 23.4.88. Thus there are positive weights \(\nu _k\), normal tensors \(B_k\), and isometries \(V_k\) with pairwise orthogonal ranges such that

\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 \\ \widetilde M_v& =\sum _kV_k(\nu _kB_k^v)V_k^\dagger . \label{eq:cpsv_prop_4_13_normal_corners} \end{align}

This completes periodic exclusion and the normal vertical block decomposition under the literal canonical-form hypothesis.

Theorem 23.4.90 groups these normal corners into actual vertical representative families over tensors \(M_\alpha \). It retains the physical isometries, both intertwining identities, exact corner compression, and the literal reconstruction. Every effective coefficient \(\mu _{\alpha ,k}=\nu _{\alpha ,k}\zeta _{\alpha ,k}\) is positive. Thus the phase-class grouping and the positivity argument of source lines 1898–1902 are complete under the literal canonical-form hypothesis.

On these actual grouped corners, Theorem 23.4.68 gives

\begin{align} (X_{\alpha ,k}^\dagger X_{\alpha ,k})M_\alpha ^v (X_{\alpha ,k}^\dagger X_{\alpha ,k})^{-1} & =M_\alpha ^v. \label{eq:cpsv_prop_4_13_gram_conjugation} \end{align}

Theorem 23.4.72 then gives a real number \(\omega _{\alpha ,k}{\gt}0\) such that

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

and \(\omega _{\alpha ,k}^{-1/2}X_{\alpha ,k}\) is unitary. These are the actual-grouped Figure 8 and gauge-normalization conclusions of Proposition 4.13. They are distinct from the tensor-attached two-site marked problem in Appendix C.4.

Theorem 23.4.74 absorbs these unitary bond gauges into the mutually orthogonal physical sector maps. The resulting maps intertwine \(\widetilde M\) with the undressed representatives and reconstruct every vertical letter exactly. Theorem 23.4.96 identifies the representatives as a basis of normal tensors for \(\widetilde M\). Taking the adjoints of the normalized sector maps as the block-row map \(U\), Theorem 23.4.93 gives (??), and the sector intertwinings give (??). Invariant-projector closure makes the retained support reducing, while the complement is the all-zero sector. Summing the sector reconstructions therefore gives the exact identity (??).

The preceding paragraphs recast the source proof on the active retained coordinates required by the literal canonical-form definition.

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}

The conclusion includes the basis-of-normal-tensors condition in [ CPGSV16 , Proposition 4.13 ] , together with the positive weights and the coisometric realization. Its hypothesis is the normalized BNT-refined form of Definition 23.4.1, rather than the literal CPSV canonical form assumed in Theorem 23.4.99.

Proof

Apply Theorem 23.4.89 once and keep its representatives \(M_\alpha \), multiplicities \(r_\alpha \), exact positive weights \(\mu _{\alpha ,q}\), and physical sector maps. For these same witnesses, Theorem 23.4.96 proves that \(\{ M_\alpha \} _\alpha \) is a BNT for \(\widetilde M\).

Theorem 23.4.73 absorbs each normalized bond gauge into its physical map and gives maps \(W_{\alpha ,q}\) such that

\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}(\mu _{\alpha ,q}M_\alpha ^v), \notag \\ \widetilde M_v & =\sum _\alpha \sum _q W_{\alpha ,q}(\mu _{\alpha ,q}M_\alpha ^v) W_{\alpha ,q}^\dagger . \notag \end{align}

The canonical equivalence between the dependent disjoint union \(\coprod _\alpha \{ 0,\ldots ,r_\alpha -1\} \) and its standard finite index set reindexes the nested reconstruction sum without changing any weight or sector. The hypotheses of Theorem 23.4.93 now hold, and that theorem gives the vertical canonical form of \(M\).