E Basis of Normal Tensors: Supporting Results
The classification of matrix product states in a basis of normal tensors (Chapter 8) relies on algebraic properties of normal tensors, Gram-matrix criteria, change-of-basis relations, sector-blocking identities, and direct-sum block separation.
E.1 Elementary normal-tensor and Gram-matrix support
Throughout this section, \(V^{(N)}(A)\) denotes the length-\(N\) matrix product vector of a tensor \(A\), and \(O_{AB}(N)\) denotes the corresponding mixed overlap. Asymptotic Gram-matrix criteria govern the linear independence of matrix product vector families in the BNT characterization and matching arguments.
Every block in a coefficient-form basis of normal tensors has positive bond dimension.
Eventual linear independence makes the matrix product vector of each block nonzero at some positive length. A tensor with zero bond dimension has zero matrix product vector at every length, which is a contradiction.
Every tensor with a one-dimensional bond space is irreducible.
The only orthogonal projections on a one-dimensional complex vector space are zero and the identity.
A tensor with one-dimensional bond space and identity transfer map is a normal tensor.
Irreducibility follows from Lemma E.1.2. The identity transfer map has spectral radius one, and its only eigenvalue is one.
Let \(\sigma \) be positive definite. The transfer map of the Perron gauge determined by \(\sigma \) is primitive if and only if the original transfer map is primitive.
The two transfer maps are similar through the invertible congruence \(X\mapsto \sigma ^{-1/2}X\sigma ^{-1/2}\). Similar linear maps have the same spectrum, so the condition that \(1\) is the only eigenvalue on the unit circle is unchanged.
Let \(A\) and \(A'\) have the same matrix product vectors at every positive length. If \((A_j)_{j=1}^{g}\) is a basis of normal tensors for \(A\), then the same family is a basis of normal tensors for \(A'\).
Substitute the positive-length MPV equality into the spanning identity for \(A\). Normality and eventual linear independence depend only on the family \((A_j)_{j=1}^{g}\).
Let \(A\) and \(B\) be normal tensors of positive, possibly different, bond dimensions. Suppose that, for every positive \(N\), their matrix product vectors are proportional by a nonzero scalar. Then there is \(a\in \mathbb {C}\), with \(|a|=1\), such that \(|V^{(N)}(A)\rangle =a^N|V^{(N)}(B)\rangle \) for every positive \(N\). This is the positive-length one-normal-block consequence of [ 13 , Corollary A.3 and Theorem 2.10 ] .
Positive-length nonzero proportionality implies eventual nonzero proportionality. Apply Theorem 8.2.12.
Let \(I\) be a finite index set and let \((A_i)_{i \in I}\) be a finite family of tensors such that \(O_{A_iA_i}(N) \to 1\) for each \(i\), and \(O_{A_iA_j}(N) \to 0\) whenever \(i \neq j\). Then the MPV states \(\{ |V^{(N)}(A_i)\rangle \} _{i \in I}\) are linearly independent for all sufficiently large \(N\).
The Gram matrix of the family \(\{ |V^{(N)}(A_i)\rangle \} _{i \in I}\) converges entrywise to the identity matrix. Hence for all sufficiently large \(N\) it is invertible, so the MPV states are linearly independent.
Stating the threshold explicitly yields an existential formulation.
If the self-overlaps of a finite block family tend to \(1\) and the cross-overlaps of distinct blocks tend to \(0\), then the MPV states are linearly independent for all sufficiently large system sizes.
This is Lemma E.1.7, restated with an explicit threshold \(N_0\).
Let \((A_j)_{j=1}^{g_A}\) and \((B_k)_{k=1}^{g_B}\) be finite families of tensors. Suppose the overlaps inside each family are asymptotically orthonormal, and suppose every mixed overlap \(O_{A_jB_k}(N)\) tends to \(0\). Then the combined MPV family
is linearly independent for all sufficiently large \(N\).
Apply the finite-index Gram-matrix criterion to the disjoint union of the two index sets. The within-family hypotheses give the diagonal and off-diagonal limits inside each block, and the mixed overlap hypothesis gives the off-diagonal limits between the two blocks.
Eventual linear independence converts equality of linear combinations of MPV states into equality of their scalar coefficients, supplying the identification used in Chapter 9.
If a finite family of vectors is linearly independent for all sufficiently large \(N\), and two finite linear combinations of that family are equal for all sufficiently large \(N\), then the corresponding coefficients are equal for all sufficiently large \(N\).
Move the two linear combinations to one side. Linear independence forces each coefficient difference to vanish at every sufficiently large length.
When two BNT families must be compared against each other, the same Gram-matrix criterion applies to their disjoint union, provided the mixed overlaps also decay. This is the form used in the exact sector matching of Section 8.8.
E.2 Change of basis, power sums, and sector bookkeeping
Linear-algebraic change of basis, geometric power-sum extrapolation, and sector bookkeeping provide the tools for comparing alternative BNT decompositions. A sector has copy weights \(\mu _{j,q}\ne 0\), coefficient \(c_{N,j}=\sum _q\mu _{j,q}^N\), representative bond dimension \(D_j\), and total bond dimension \(D_{\mathrm{tot}}=\sum _j r_jD_j\).
Let \(v_1, \ldots , v_g\) and \(w_1, \ldots , w_g\) be two linearly independent families in a complex vector space with the same span. Then there exists an invertible matrix \(U \in \mathrm{GL}_g(\mathbb {C})\) such that \(w_j = \sum _{i=1}^g U_{ij} v_i\) for every \(j\).
Because the two families span the same subspace, each \(w_j\) has a unique expansion in the basis \(\{ v_i\} _{i=1}^g\); these coefficients form the columns of \(U\). If \(U a = 0\), then the corresponding linear relation among the \(w_j\) is trivial because the \(w_j\) are linearly independent. Thus \(U\) is injective, hence invertible.
If two finite families of blocks have asymptotically orthonormal overlaps within each family, and if their MPV spans agree for all sufficiently large system sizes \(N\), then for all sufficiently large \(N\) there exists an invertible matrix \(U_N\) expressing the MPV states of one family as linear combinations of the MPV states of the other.
Let \(w_1,\ldots ,w_n\) be nonzero complex numbers. If \(\sum _i c_i w_i^N=0\) for every sufficiently large \(N\), then the same equality holds for every non-negative exponent.
Induct on the number of terms. Subtracting \(w_1\) times the equality at exponent \(N\) from the equality at exponent \(N+1\) gives \(\sum _{i=2}^{n} c_i(w_i-w_1)w_i^N=0\). The induction hypothesis yields this equality at every exponent, so the original sum \(S_N=\sum _i c_iw_i^N\) satisfies \(S_{N+1}=w_1S_N\) for every \(N\). Since \(S_N\) eventually vanishes and \(w_1\ne 0\), it follows that \(S_0=0\) and hence \(S_N=0\) for every \(N\).
Let \(a_1,\ldots ,a_m\) and \(b_1,\ldots ,b_n\) be nonzero complex numbers. If their power sums agree for every sufficiently large exponent, then \(\sum _i a_i^k=\sum _j b_j^k\) for every non-negative exponent \(k\).
Set \(w_l=a_l\) for \(l\leq m\) and \(w_{m+l}=b_l\) for \(l\leq n\), with coefficients \(c_l=1\) for \(l\leq m\) and \(c_l=-1\) for \(l{\gt}m\). Then \(\sum _{l=1}^{m+n}c_lw_l^k=\sum _i a_i^k-\sum _j b_j^k\). The hypothesis makes this expression zero for every sufficiently large \(k\). Lemma E.2.3 makes it zero for every \(k\), which is the required power-sum equality.
Let \(\alpha , \beta : \{ 0, \ldots , n-1\} \to \mathbb {C}\). If \(\sum _i \alpha _i^k = \sum _i \beta _i^k\) for all \(k \ge 1\), then the multisets \(\{ \alpha _i\} \) and \(\{ \beta _i\} \) are equal.
Construct diagonal matrices \(\operatorname{diag}(\alpha )\) and \(\operatorname{diag}(\beta )\). The power-sum hypothesis gives \(\operatorname{tr}(\operatorname{diag}(\alpha )^k) = \operatorname{tr}(\operatorname{diag}(\beta )^k)\). By the corresponding theorem in the companion quantum-channel volume [ 216 , “Newton–Girard trace recursion” ] , the characteristic polynomials agree: \(\prod _i (X - \alpha _i) = \prod _i (X - \beta _i)\). Extracting roots gives multiset equality.
Let the representative \(A_j\) have bond dimension \(D_j\), and let \(D\) be the bond dimension of the flattened sector tensor. Then \(D_j\leq D\). If every \(D_j\) is positive and there are \(g\) representatives, then also \(g\leq D\).
Since \(r_j{\gt}0\), choose one copy of \(A_j\). Its contribution \(D_j\) is one non-negative summand in the total dimension of the flattened copy family, so \(D_j\leq D\).
For the second assertion, every representative has multiplicity \(r_j\geq 1\), and each copy has positive dimension \(D_j{\gt}0\). Therefore
The coefficient of the blocked decomposition satisfies
Moreover, blocking the tensor assembled from \(P\) gives the same MPV family as the tensor assembled from the blocked sector decomposition.
A nonempty finite family of blocks has at least one MPV phase class.
Choose one original block index. Its occurrence in the phase-class enumeration supplies a phase class containing it.
Let \(F_k\) be a finite family of square matrices indexed by the original blocks. If \(k=(j,q)\) enumerates the copies in the MPV phase classes, then \(\sum _j\sum _q F_{j,q}=\sum _k F_k\).
Evaluate both sides at a pair of matrix indices and apply the bijective enumeration of the finite phase-class quotient.
Given nonzero weights \((\mu _k)_{k=1}^r\) and blocks \((B_k)_{k=1}^r\), the single-copy sector decomposition has \(r\) sectors, one basis block \(B_k\) per sector, multiplicity \(1\), and sector weight \(\mu _k\). As a finite sum,
The single-copy sector decomposition of \((\mu _k,B_k)_k\) represents the same MPS family as \(\bigoplus _k\mu _kB_k\).
Expanding the sector decomposition gives the left-hand side of (2), which is the finite block-sum expression for \(\bigoplus _k\mu _kB_k\).
For any sector decomposition \(P\) and any sector \(j\), the coefficient \(c_{P,j}^{(N)}=\sum _q\mu _{j,q}^N\) is not eventually zero in \(N\). More generally, for every positive integer \(p\), the subsequence \(c_{P,j}^{(pN)}\) is not eventually zero.
If the coefficient vanished for every \(N\geq M\), the finite geometric-sum extrapolation would make it vanish at every non-negative exponent. However, at exponent zero, \(c_{P,j}^{(0)}=\sum _q\mu _{j,q}^0=r_j{\gt}0\), which is impossible. For the subsequence, apply the same argument to the nonzero bases \(\mu _{j,q}^p\).
E.2.1 Proportional matching and finite bijection extraction
The equal-MPV matching theorems in Chapter 8 specialize an eventual proportionality argument; matching in both directions produces a bijection between the sector sets.
Let \(P\) and \(Q\) be BNT canonical forms whose total MPV families are eventually nonzero and proportional. For every sector \(j\) of \(P\), there is a sector \(k\) of \(Q\) with equal bond dimension, gauge-phase-equivalent basis tensor after the dimension identification, and non-decaying cross-overlap.
Suppose the chosen sector decays against every sector of \(Q\), and separate all such \(P\)-sectors from the remaining ones. For each remaining sector, choose a \(Q\)-sector with non-decaying overlap. Both blocks are normal by the canonical-form hypotheses. Theorem 8.2.9 rules out its decay alternative and supplies equal bond dimensions and a gauge-phase relation with a unit phase. Gauge-phase covariance therefore expresses that \(P\)-sector MPV as a geometric scalar multiple of the chosen \(Q\)-sector MPV.
Substitute these relations into the eventual proportionality identity. The combined-family Gram criterion makes the decaying residue together with the \(Q\) family eventually linearly independent. Exact coefficient extraction then makes the chosen sector coefficient eventually zero, contradicting its power-sum nonvanishing. Thus a non-decaying partner exists. A second application of Theorem 8.2.9 rules out decay and gives the required dimension equality and gauge-phase equivalence; no modulus-one overlap limit is assumed.
Under the same hypotheses, every sector of \(Q\) has an equal-dimensional, gauge-phase-equivalent partner in \(P\) with non-decaying cross-overlap.
Reverse the eventual proportionality and apply Theorem E.2.1.1 to each sector of \(Q\); symmetry of the overlap converts the resulting witness to the stated orientation.
Suppose every sector of \(Q\) has an equal-dimensional, gauge-phase-equivalent partner in \(P\) with non-decaying overlap, and the converse holds from \(P\) to \(Q\). Then the forward partners may be chosen as a bijection between the two finite sector sets, preserving all three properties.
Basis separation makes each partner map injective: two sectors with the same partner would be gauge-phase equivalent to one another. The injections in both directions force equal finite cardinalities, so the forward map is bijective.
E.3 Direct-sum separation and prepared-family bookkeeping
Following Pérez-García et al. [ 228 ] , disjoint local MPS spaces yield pair-trace selectors whose products isolate individual BNT blocks. Grouping phase-equivalent blocks into sectors preserves total bond dimension and allows projective weight normalization.
Let \(A^1,\ldots ,A^r\) be irreducible left-canonical representatives with normalized self-overlaps, and suppose that distinct same-dimensional representatives are not gauge-phase equivalent. Fix \(L_0{\gt}0\) such that every block is injective at lengths \(L_0\), \(L_0+1\), and \(3(L_0+1)\). Then each ordered pair \(k\ne j\) is trace-separated at the common length \(S=3(L_0+1)\): the map
is onto \(M_{D_k}(\mathbb {C})\oplus M_{D_j}(\mathbb {C})\). In particular, one may prescribe an arbitrary matrix on either block while making the same word combination vanish on the other block.
Put \(L=L_0+1\) and \(S=3L\). For a block \(A\), write
Block injectivity at \(L\) and \(S\) makes both boundary maps injective; in particular,
Fix \(k\ne j\) and first order the pair so that \(D_j\leq D_k\). Suppose that a nonzero vector belongs to \(\mathcal G_S(A^k)\cap \mathcal G_S(A^j)\). There are boundary matrices \(X\) and \(Y\) such that
Injectivity at \(S\) gives \(X\ne 0\) and \(Y\ne 0\). Split every length-\(S\) word as \(uvw\) with \(|u|=|v|=|w|=L\). Relation (4) gives
Because the length-\(L\) words of \(A^k\) span \(M_{D_k}(\mathbb {C})\) and \(X\ne 0\), the two-sided ideal generated by \(X\) is all of \(M_{D_k}(\mathbb {C})\):
Indeed, expand arbitrary left and right matrix units in length-\(L\) words; a nonzero entry of \(X\) then produces every matrix unit. Thus, for each \(Z\in M_{D_k}(\mathbb {C})\), choose coefficients \(a_{uv}\) with \(Z=\sum _{u,v}a_{uv}A^k_uXA^k_v\) and set \(R(Z)=\sum _{u,v}a_{uv}A^j_uYA^j_v\). Cyclicity of the trace and (5), applied to the word \(vwu\), yield
Hence \(\Gamma _L^{A^k}(Z)=\Gamma _L^{A^j}(R(Z))\) for every \(Z\), and therefore
Together with \(D_j\leq D_k\), this forces \(D_j=D_k\) and equality in (7). This proves the unequal-dimension branch immediately: if \(D_j{\lt}D_k\), a nonzero intersection is impossible. Reversing the initial ordering treats \(D_k{\lt}D_j\) in the same way.
It remains to exclude the intersection when \(D_j=D_k\). Gauge-phase separation and Theorem 8.4.3 give, beyond any prescribed threshold, a chain length \(N{\gt}L\) for which the two periodic MPVs are not proportional. On the other hand, equality \(\mathcal G_L(A^k)=\mathcal G_L(A^j)\) gives identical translated length-\(L\) constraints on the \(N\)-site ring. The parent-Hamiltonian uniqueness mechanism identifies the corresponding chain ground spaces with the one-dimensional MPV spans:
Equality of the local spaces therefore forces these two MPV lines to be equal, contradicting the choice of \(N\). Consequently, for every distinct pair, including the unequal-dimensional pairs,
Finally, let \(\Phi _{kj}\) be the displayed length-\(S\) word-evaluation map from coefficient families to \(M_{D_k}(\mathbb {C})\oplus M_{D_j}(\mathbb {C})\). Pair the target by
A pair \((X,Y)\) annihilates \(\operatorname {im}(\Phi _{kj})\) precisely when
or equivalently \(\Gamma _S^{A^k}(X)=\Gamma _S^{A^j}(-Y)\). There is therefore a linear map
It is injective because both boundary maps are injective. Conversely, if \(\psi \) lies in the intersection, there are unique \(X\) and \(Z\) with \(\Gamma _S^{A^k}(X)=\psi =\Gamma _S^{A^j}(Z)\); the pair \((X,-Z)\) annihilates \(\operatorname {im}(\Phi _{kj})\) and maps to \(\psi \). This construction is the inverse of \(\Theta \), so the two spaces are isomorphic:
The right-hand side vanishes by (8), and hence so does the annihilator. The trace pairing is nondegenerate and all spaces are finite-dimensional, so \(\operatorname {im}(\Phi _{kj})^\perp =\{ 0\} \) implies that \(\Phi _{kj}\) is surjective. Surjectivity first gives an arbitrary value on one block and zero on an anchor block; the selector construction then chooses \((\mathbb {1},0)\) for the remaining pairs, and products of these homogeneous selectors give the simultaneous word-tuple span in Lemma 8.3.3.
Let \((B_k)_{k=1}^r\) be a finite family of positive-bond-dimension MPS tensors, each trace-preserving, primitive, and tensor-irreducible. Then there is a single \(L\ge 1\) such that each \(B_k^{[L]}\) is one-site injective, and this blocking preserves trace preservation, primitivity, and irreducibility for every block.
Suppose CF summands are grouped by BNT representative, every original copy weight is nonzero, and every summand assigned to a representative has that representative’s bond dimension. Then the total bond dimension of the resulting BNT sector decomposition is the sum of the original block dimensions.
The collapsed decomposition counts one representative dimension for each assigned copy. Let \(e(j,q)\) denote the original CF summand assigned to the \(q\)th copy of representative \(P_j\). Since \((j,q)\mapsto e(j,q)\) reindexes the original summands and each assigned summand has the representative’s bond dimension,
This is the equality between the collapsed total dimension and the original direct-sum dimension.
Let \(X\) and \(Y\) be trace-preserving, primitive, irreducible NTs of positive bond dimension. If their MPV families differ by a nonzero scalar power, then \(X\) and \(Y\) have the same bond dimension.
Choose \(\zeta \neq 0\) with \(|V^{(N)}(Y)\rangle =\zeta ^N|V^{(N)}(X)\rangle \) for every positive length \(N\). The normalized self-overlaps give \(|O_{XX}(N)|\to 1\), and the phase relation gives
Since the same normalization forces \(|\zeta |=1\), this gives \(|O_{YX}(N)|\to 1\). If the two bond dimensions differed, the rectangular overlap-decay theorem would force \(O_{YX}(N)\to 0\), contradicting the preceding limit.
In a prepared family of positive-bond-dimension trace-preserving primitive irreducible NTs, every CF summand assigned to a BNT representative has that representative’s bond dimension.
Apply Lemma E.3.4 to each representative and each assigned summand.
For prepared trace-preserving primitive irreducible NTs of positive bond dimension with all copy weights nonzero, grouping CF summands by BNT representative preserves the total bond dimension: the collapsed sector decomposition has total bond dimension equal to the sum of the original block dimensions.
For nonzero weights \(\mu _1,\ldots ,\mu _r\) with \(r\geq 1\) there is a positive real \(m\) such that every normalized weight satisfies
and at least one normalized weight has unit modulus, \(|\mu _{k_0}/m|=1\). This is the normalization that [ 13 , line 246 ] chooses for the canonical-form weights.
Take \(m=\max _k|\mu _k|\); the maximum over a nonempty finite family is attained, positive because every weight is nonzero, an upper bound for every modulus, and equal to the modulus of the attaining weight.
Multiplying every block weight in a block-diagonal MPS tensor by a scalar \(c\) multiplies every length-\(N\) MPV coefficient by \(c^N\).
Expand the block-diagonal MPV coefficient as a finite sum over blocks and factor \(c^N\) out of the sum.