---
title: Disordered MPS Decomposition
url: https://www.emergentmind.com/topics/disordered-matrix-product-state-decomposition
type: topic
---

# Disordered MPS Decomposition

Searching arXiv for the core paper and closely related work on disorder-adapted MPS/MPO representations.
Searching for the 2014 MBL–MPO paper.
Disordered matrix product state decomposition denotes a family of tensor-network representations in which disorder, inhomogeneity, or realization-dependent local structure is encoded directly into matrix product state (MPS) or matrix product operator (MPO) tensors. In the most specific and influential usage, the decomposition refers to the claim that in a one-dimensional many-body localized (MBL) phase the exact diagonalizing unitary \(U\) can be efficiently compressed as a low-bond-dimension MPO, so that the entire many-body spectrum is generated by binary local tensor choices \(G_i\) and \(E_i\) at each site [1410.2224]. In broader usage, the same phrase also covers site-dependent random MPS ensembles, fixed-bond-dimension MPS manifolds used for disordered dynamics, and disorder-sensitive tensor decompositions in phases without conventional local order [2103.02634]. Across these settings, the common idea is that disorder does not merely perturb a clean tensor network; it selects a nonuniform local basis, alters entanglement structure, and changes what constitutes an efficient decomposition.

## 1. Spectrum-wide decomposition in many-body localized phases

The most explicit formulation of a disordered MPS decomposition appears in the proposal to represent the full eigensystem of an MBL Hamiltonian by a single disorder-adapted MPO \(U\) [1410.2224]. The physical motivation is the analogy with Anderson insulators: in a noninteracting localized system, once the localized single-particle orbitals are known, the many-body eigenstates are product states in the occupation basis of those orbitals. The interacting analogue replaces localized orbital occupations by localized tensor choices.

The corresponding MBL ansatz is
\[
\Psi_\text{MBL}=\prod_i (\psi_i E_i^{\sigma_i}+\phi_i G_i^{\sigma_i})|\sigma_i\rangle,
\]
where \(G_i\) is the local tensor corresponding to the local ground-state choice on site \(i\), and \(E_i\) is the local tensor corresponding to the local excited-state choice [1410.2224]. For an \(n\)-site spin-\(\tfrac12\) chain, all \(2^n\) eigenstates are obtained by all binary choices of \(G_i\) or \(E_i\) across the chain. This is not merely a statement that each eigenstate is individually approximable by an MPS. Rather, the entire spectrum is organized by one common low-bond-dimension MPO.

The operator \(U\) is defined by
\[
U|p\rangle=|e_p\rangle,\qquad U=\sum_p |e_p\rangle\langle p|,
\]
where \(|p\rangle\) are product states in the microscopic \(\sigma^z\) basis and \(|e_p\rangle\) are the exact many-body eigenstates [1410.2224]. If \(U\) is represented as an MPO with local tensors \(O_i^{\chi_i\sigma_i}\), then acting on a product basis state selects one of two tensor slices, \(O_i^{\downarrow\sigma_i}\) or \(O_i^{\uparrow\sigma_i}\), which are relabeled as \(G_i\) and \(E_i\). In this way, a disordered tensor decomposition becomes a realization-dependent reorganization of the full spectrum into localized binary degrees of freedom.

This formulation is significant because MBL is not primarily a ground-state phenomenon. The relevant structural question concerns highly excited states at finite energy density. The proposal therefore extends tensor-network representability beyond the conventional gapped-ground-state setting and ties it to the existence of quasi-local conserved quantities, area-law entanglement, and failure of thermalization [1410.2224].

## 2. Tensor-network structure, l-bits, and bond-dimension diagnostics

The efficiency of the decomposition is controlled by the MPO bond dimension \(D\). Numerically, the MPO for \(U\) is obtained by gluing bra and ket indices into a doubled local index,
\[
\mu_i=2*\sigma_i+\chi_i,
\]
compressing the resulting object as an MPS, and then splitting \(\mu_i\) back into \(\sigma_i,\chi_i\) to recover the MPO [1410.2224]. The underlying expectation is that if \(U\) is quasi-local, then entanglement in the doubled-space representation remains short-ranged, so the required bond dimension stays small.

A central implication is the entanglement bound
\[
S\lesssim \log D,
\]
since the entanglement entropy across a cut is bounded by the logarithm of the Schmidt rank and hence by \(\log D\) [1410.2224]. If \(D\) remains bounded in the thermodynamic limit, then all eigenstates generated by the MPO obey an area law and cannot satisfy ETH volume-law scaling. This makes bond-dimension growth itself a localization diagnostic.

The l-bit structure is represented directly in MPO language. The operator
\[
U\sigma_i^+U^\dagger
\]
is identified as the l-bit raising operator, and a product-basis operator such as
\[
\rho_\textrm{product}=I_1\otimes\cdots\otimes I_{k-1}\otimes \sigma_k^z\otimes I_{k+1}\otimes\cdots\otimes I_n
\]
is mapped to a conserved quasi-local operator
\[
U\rho_\textrm{product}U^\dagger=\sum_i \alpha_i |e_i\rangle\langle e_i|
\]
in the physical basis [1410.2224]. This provides a tensor-network realization of emergent localized integrals of motion.

The same work argues that if the localization length is \(l\), then the MPO bond dimension scales roughly as
\[
D\lesssim 4^l.
\]
For strong disorder, \(l\) is finite and \(D\) should saturate with system size; in a delocalized phase \(D\) should grow strongly with \(L\) [1410.2224]. The numerical fit
\[
\log_4[D(L)] = a\tanh(L/\xi)
\]
introduces a length scale \(\xi\) interpreted as a saturation length or localization length, and this \(\xi\) was reported to agree with the localization length extracted from effective l-bit couplings up to an overall rescaling factor of about \(8.7\) [1410.2224].

## 3. Disorder, rare regions, and the limits of compactness

The canonical model used to test the decomposition is the random-field Heisenberg chain
\[
H=\sum_{\langle i,j\rangle}\sigma_i\cdot \sigma_j+\sum_i h_i \sigma_i^z,\qquad h_i\in[-\Delta,\Delta],
\]
which is noted to have an MBL transition near \(\Delta\sim 3.5\) [1410.2224]. The construction proceeds by exact diagonalization in each fixed total \(S_z\) sector, matching product states to eigenstates using the Hungarian algorithm to maximize \(\sum_i |\langle M(i)|i\rangle|^2\), and compressing the resulting unitary into an MPO. Because the matching is heuristic rather than globally optimal, the resulting bond dimensions are upper bounds on the optimal compression [1410.2224].

The reported finite-size behavior separates three regimes. For weak disorder, roughly \(\Delta\lesssim 3\), the bond dimension grows linearly with \(L\). For strong disorder, roughly \(\Delta\gtrsim 15\), the bond dimension has essentially saturated already by \(L=12\). For intermediate disorder \(3\lesssim \Delta\lesssim 15\), accessible sizes do not support a definitive asymptotic conclusion [1410.2224]. This establishes the decomposition as efficient only in the localized regime.

An important qualification is that the localized phase is not represented by a uniformly bounded bond dimension over all disorder samples. The paper emphasizes strong Griffiths effects: rare weak-disorder regions can resemble thermal inclusions and require much larger local bond dimension [1410.2224]. Accordingly, the bond-dimension distribution in the localized phase is peaked at small \(D\), around \(D\sim 4^\xi\), but has a power-law tail; near the transition the distribution becomes extremely broad. This means that the decomposition is typical rather than uniform. A plausible implication is that disorder-adapted MPS descriptions of MBL should be understood probabilistically across realizations rather than as deterministic worst-case bounds.

The same theme reappears in later hybrid and dynamical contexts. In disordered electron–phonon systems treated by a hybrid MPS–multi-trajectory Ehrenfest scheme, disorder enters as quenched onsite potentials \(\epsilon_\ell\in[-W,W]\) plus a trajectory-dependent field \(-\sqrt{2}\gamma x_\ell(t)\), so that along one trajectory the effective onsite potential is
\[
\epsilon_{\mathrm{total},\ell}(t)=\epsilon_\ell-\sqrt{2}\gamma x_\ell(t).
\]
The MPS decomposition itself remains standard, but the numerics become dominated by averaging over both disorder realizations and classical trajectories, with \(N_\mathrm{traj}=200\) and \(N_\mathrm{dis}=200\) to \(2000\) in the main disordered calculations [2512.10899]. This suggests that in disordered tensor decompositions the primary complication may shift from ansatz design to sampling and realization dependence.

## 4. Random and site-dependent MPS as disordered tensor ensembles

A different meaning of disordered matrix product state decomposition arises in the theory of random matrix product states with site-dependent local tensors [2103.02634]. Here the state is written as a periodic MPS
\[
|\psi\rangle=\sum_{i_1,\dots,i_n}\operatorname{tr}\!\left[A_{i_1}^{(1)}A_{i_2}^{(2)}\cdots A_{i_n}^{(n)}\right]|i_1,\dots,i_n\rangle,
\]
with \(A_i^{(k)}\in\mathbb C^{D\times D}\) explicitly dependent on site \(k\) [2103.02634]. Disorder means that the tensors are not translation invariant and are generated from independent Haar-random local unitaries \(U^{(k)}\in U(dD)\). The resulting ensemble is denoted \(\mu_{d,n,D}\).

This framework is not an algorithm for decomposing a given state. It is instead a formal ensemble of disordered tensor decompositions. The local tensors are obtained by feeding a fixed state \(|0\rangle\in\mathbb C^d\) into a Haar-random unitary acting on \(\mathbb C^d\otimes\mathbb C^D\), equivalently by drawing tensor cores uniformly from the Stiefel manifold of isometries [2103.02634]. Because the disorder is local and i.i.d., global moments factorize into products of local Haar integrals.

The analytical machinery is built from Weingarten calculus and a mapping of moment tensor networks to one-dimensional classical partition functions. For example, moments such as
\[
\mathbb E|\langle\psi|\phi\rangle|^4
\]
are reduced to sums over local permutation variables in \(S_2=\{\mathbb 1,\mathbb F\}\), while subsystem purities become effective one-dimensional statistical mechanics models with local plaquette weights [2103.02634]. This converts a disordered tensor-network problem into a tractable transfer-weight problem.

Several structural consequences follow. The norm concentrates near one,
\[
\mathrm{Pr}\left(|\langle\psi|\psi\rangle-1|\ge \varepsilon\right)\le \varepsilon^{-2} d^{-n},
\]
and for Hamiltonians with non-degenerate spectrum and non-degenerate spectral gap the normalized state equilibrates exponentially well with overwhelming probability [2103.02634]. For disconnected subsystems, the Rényi-2 entropy is generically extensive, while for connected subsystems of fixed size and sufficiently large bond dimension the reduced state is almost maximally mixed [2103.02634]. These results identify a disordered MPS ensemble as a model of generic states in the trivial phase of matter, with disorder expressed entirely through site-dependent random local tensors.

## 5. Fixed-bond MPS manifolds and disordered dynamics

In dynamical problems, disordered MPS decomposition often means not a static factorization theorem but a restriction of quantum evolution to a fixed-bond-dimension MPS manifold. The time-dependent variational principle realizes this as
\[
\frac{d}{dt}\ket{\Psi[A]}=-iP_{\mathcal M}\hat H\ket{\Psi[A]},
\]
where \(P_{\mathcal M}\) projects the Schrödinger vector field onto the tangent space of the fixed-bond-dimension MPS manifold \(\mathcal M_{\mathbf r}\) [1710.09378]. The relevant decomposition is therefore the split of \(\hat H|\Psi[A]\rangle\) into tangent and discarded components.

The principal result is methodological rather than structural: apparent convergence of long-time observables with bond dimension can be misleading [1710.09378]. In benchmark cases, TDVP produces incorrect asymptotic behavior even when observables appear stable. For the clean XX chain, exact ballistic transport crosses over under TDVP to apparently diffusive behavior with diffusion constant \(\sim 2\); for the disordered XX chain, exact Anderson localization is missed, and the mean-square displacement fails to plateau; for the clean XX ladder, the correct \(D\sim0.95\) crosses over to an incorrect \(D\sim0.2\) [1710.09378].

Against this cautionary backdrop, the disordered nonintegrable XXZ chain
\[
\hat H=J_{xy}\sum_{i=1}^{N-1}(\hat S_i^x\hat S_{i+1}^x+\hat S_i^y\hat S_{i+1}^y)+\Delta\sum_i \hat S_i^z\hat S_{i+1}^z+\sum_{i=1}^N h_i\hat S_i^z,
\qquad h_i\in[-W,W],\quad J_{xy}=1,
\]
with \(\Delta=1\) and \(W=1.5\), behaves differently [1710.09378]. The mean-square displacement
\[
\sigma^2(t)=\sum_{i=1}^{L}\left(\frac{L}{2}-i\right)^2\left\langle \hat S_i^z(t)\hat S_{L/2}^z(0)\right\rangle
\]
obeys \(\sigma^2(t)\sim t^\alpha\) with fitted exponent \(\alpha=0.54\), in excellent agreement with short-time exact and previously numerically exact results [1710.09378]. This does not establish a new disordered MPS decomposition, but it does delimit when a fixed-complexity disordered MPS representation may be asymptotically faithful for transport.

A common misconception is that weak bond-dimension dependence of long-time observables suffices to validate the method. The explicit counterexamples show that the limits \(t\to\infty\) and \(\chi\to\infty\) need not commute [1710.09378]. In this sense, disordered MPS manifolds can be reliable dynamical coarse-grainings in some interacting disordered regimes, but not as a universal principle.

## 6. Other uses of the term and broader methodological context

Outside quenched-random spin chains, the phrase also appears in settings where “disordered” refers to absence of local order or to general site dependence rather than to random couplings. In the one-dimensional extended quantum compass model, the disordered phases are efficiently represented by low-bond-dimension MPS and are characterized by a doubly degenerate entanglement spectrum and nonzero string order parameters rather than by conventional local order parameters [1201.1435]. In the limiting disordered phases with \(L_1=0\), the exact ground states have bond dimension \(\chi=2\) and doubly degenerate diagonal bond matrix \(\Lambda=I_{2\times 2}\), providing a simple example in which the MPS decomposition exposes hidden order [1201.1435].

In multimode cavity QED, an MPS-compatible chain decomposition is constructed not for a random many-body state but for an inhomogeneous bosonic environment. Numerical mode decomposition yields arbitrary cavity mode frequencies and couplings, which are then transformed from star geometry to a chain Hamiltonian by an exact orthogonal map
\[
\hat b_n=\sum_{k=1}^{M} U_{n,k}\hat a_k,
\qquad
\sum_k \omega_k \hat a_k^\dagger \hat a_k
=
\sum_{n=1}^{M}\left[\xi_n\hat b_n^\dagger \hat b_n+t_n(\hat b_n^\dagger \hat b_{n+1}+\hat b_{n+1}^\dagger \hat b_n)\right].
\]
The resulting state is then represented as an MPS along the chain [2212.01935]. This is not a disorder paper, but it shows how nonuniform mode structure is converted into an MPS-compatible decomposition.

More abstractly, site-dependent MPS can also be viewed through algebraic and compilation lenses. The algebraic-geometric theory of translation-invariant MPS provides an instructive contrast: many simplifications there depend on cyclic symmetry and fail in genuinely disordered settings [1210.2812]. By contrast, shallow-circuit decompositions of generic site-dependent MPS do not assume translation invariance and therefore extend naturally to inhomogeneous inputs; the best-performing protocol was a sequential analytic circuit growth plus global constrained unitary optimization strategy [2209.00595]. A plausible implication is that disorder primarily obstructs symmetry-based reductions, not tensor-network compilation in principle.

Taken together, these works support a broad but technically precise interpretation of disordered matrix product state decomposition. In one line of work, it is a disorder-adapted full-spectrum factorization of localized eigenstates by an MPO of the diagonalizing unitary [1410.2224]. In another, it is a formal ensemble of site-dependent random tensors with analyzable entropic and equilibration properties [2103.02634]. In a third, it is a constrained manifold for long-time disordered dynamics whose reliability is regime dependent [1710.09378]. The unifying content is that disorder changes the natural tensor basis, the relevant bond-dimension diagnostics, and the distinction between typical and worst-case representability.

Source: https://www.emergentmind.com/topics/disordered-matrix-product-state-decomposition