---
title: Matrix Product State Approximation
url: https://www.emergentmind.com/topics/matrix-product-state-mps-approximation
type: topic
---

# Matrix Product State Approximation

A matrix product state (MPS) approximation is a cornerstone technique in quantum many-body physics, tensor networks, and quantum information theory for efficiently representing and simulating states of extended 1D systems. By expressing highly correlated wavefunctions as products of low-rank tensors with controlled entanglement, MPS provides both an explicit variational ansatz and an algorithmic framework with theoretically quantifiable error guarantees. The approximation properties of MPS are pivotal for describing ground states of gapped, local Hamiltonians, Gaussian/Bogoliubov vacuum states, critical models, infinite systems, and even quantum fields.

## 1. Structure and Formal Statement of the Approximation Problem

Consider an infinite or finite one-dimensional quantum lattice with a local Hilbert space of dimension $d$ at each site. An MPS of bond dimension $D$ is defined by site tensors $A^{[i]\,s_i}$ (with virtual indices of size up to $D$), and the many-body Hilbert space vector is assembled as
$$
|\Psi[A]\rangle = \sum_{s_1,\dots,s_L} \mathrm{tr}(A^{[1]s_1}A^{[2]s_2}\cdots A^{[L]s_L})\,|s_1s_2\ldots s_L\rangle
$$
for open or periodic boundary conditions. In the thermodynamic (infinite) limit, translationally invariant MPS (iMPS) are defined by a single tensor $A^s_{\alpha\beta}$ repeated at each site.

The MPS approximation problem is to find, given a target state $|\Psi\rangle$ (e.g., the ground state of a local Hamiltonian or a correlated Gaussian state), an MPS $|\Phi_D\rangle$ such that a specified cost function is minimized. Standard cost functions include
- the $L^2$-norm $\| |\Psi\rangle - |\Phi_D\rangle\|_2$ for pure states,
- the trace-norm $\| \sigma_L - \rho_L \|_\mathrm{tr}$ between reduced density matrices on $L$ contiguous sites,
- or the maximum deviation of local observable expectation values $|\langle\Psi|O|\Psi\rangle - \langle \Phi_D |O| \Phi_D \rangle|$ for all $O$ supported on $L$ sites.

The central theoretical result for infinite, gapped, translation-invariant 1D spin chains with unique ground states is that any local observable on $L$ sites can be approximated to accuracy $\epsilon$ by an iMPS of bond dimension $D = O( (L-1)/\epsilon )$, up to subpolynomial corrections, with rigorous constants depending on the spectral gap and local dimension [1711.06559].

## 2. Main Theorems and Error Bounds

The precise quantitative theorem for infinite 1D gapped systems [1711.06559] states:

*Given a nearest-neighbor, translationally invariant, gapped Hamiltonian $H = \sum_n h_{n,n+1}$ on an infinite chain with local dimension $d$ and spectral gap $\Delta > 0$, and with a unique translation-invariant ground state $|\Psi\rangle$, for any $L$ and $\epsilon > 0$, there exists a translationally invariant iMPS $|\Phi_D\rangle$ with*
$$
D(L,\epsilon) \leq \frac{6(L-1)}{\epsilon}\;\exp\left(
\frac{2}{c_2^{3/4} \Delta^{1/4} [\log\frac{(L-1)\,\tilde c_1}{\epsilon^3}]^{3/4}}
\right)
$$
*such that for any local region of $L$ consecutive spins, the reduced density matrices $\sigma_L$ (of $|\Psi\rangle$) and $\rho_L$ (of $|\Phi_D\rangle$) satisfy*
$$
\| \sigma_L - \rho_L \|_\mathrm{tr} \leq \epsilon.
$$
This guarantees local indistinguishability up to error $\epsilon$ for observables supported on $L$ sites. The scaling in $D$ is dominated by $D=O((L-1)/\epsilon)$; the exponential correction is subpolynomial in $1/\epsilon$ for fixed $L,\Delta$.

The proof structure leverages area-law bounds on the decay of the Schmidt coefficients $\{\lambda_i\}$ of $|\Psi\rangle$ across any cut (specifically, $\log\sum_{i>\chi}\lambda_i^2 \leq \log c_1 - c_2 \Delta^{1/3} (\log \chi)^{4/3}$), truncates the chain over a finite block of size $k$, assembles a local MPO approximation, and restores translation invariance by averaging. The two error sources—tail truncation and block-tiling—are balanced to optimize $D$.

For ground states of general 1D gapped Hamiltonians (not necessarily translationally invariant), the "locally accurate" MPS approximation result [1903.10241] gives that for any $k$ and $\epsilon$ there exists an MPS of bond dimension $\chi=(k/\epsilon)^{1+o(1)}$ that is $\epsilon$-close in trace-norm for every contiguous $k$-site block, independently of the total system size.

## 3. Algorithmic Constructions and Variational Methods

The explicit construction of an MPS approximation typically proceeds by:
- variational energy minimization over the MPS ansatz (DMRG/variational MPS), usually using two-site or one-site sweeps, with local tensor updates and SVD-based truncation controlling the bond dimension and truncation error [1310.4118];
- for Gaussian states (Hartree–Fock–Bogoliubov/Bogoliubov vacua), constructing the optimal MPS by sequential site-wise Schmidt decompositions, with overlaps given by Pfaffians for fermionic Gaussian states, resulting in closed-form MPS tensors and guaranteed optimality for a given bond dimension [2111.09101, 1504.07701].

For infinite, translation-invariant systems, tangent-space-based variational algorithms project gradient steps onto the MPS manifold, yielding optimal truncations of infinite MPS under the manifold geometry [2001.11882]. These methods control the overlap per site (i.e., the fidelity density) and maintain optimality relative to the fixed bond dimension.

Continuous quantum fields are addressed with continuous MPS (cMPS) ansätze, using path-ordered exponentials of parameter matrices and direct energy-minimization via gradient or conjugate-gradient updates in the central canonical gauge [1611.03779].

## 4. Applications, Extensions, and Sample Complexity

MPS approximations are fundamental for:
- efficient simulation of low-energy physics in 1D lattice field theories (e.g., mass spectra and condensates in the Schwinger model), where variational MPS/DMRG methods yield systematic error control via bond dimension, system size, and lattice spacing extrapolations [1310.4118];
- tomographic and quantum information protocols, e.g., Sketch Tomography, which reconstructs the MPS/tensor-train representation of a quantum state from classical shadow data with $O(n^2\,\mathrm{poly}(1/\epsilon))$ sample complexity, outperforming generic shadow methods on observables involving moderate or large subsystems [2512.03333].

In quantum state preparation for quantum devices, MPS with moderate bond dimension give a route to load structured classical data (piecewise polynomials, wavelet-compressed images) using low-depth quantum circuits constructed via Matrix Product Disentangler and subsequent tensor-network optimization [2502.16464].

In models with exact MPS ground states ("MPS skeletons") such as Onsager-integrable quantum chains, dense networks of explicit analytical MPS populate entire gapped phases, enabling exponential convergence of the MPS approximation in bond dimension to the true ground-state energy and correlators [2511.07212].

## 5. Theoretical Significance: RG Interpretation and Limitations

The MPS approximation naturally emerges from the area law for entanglement entropy in 1D systems. Truncating bond dimension corresponds to discarding Schmidt values below a cutoff, paralleling Wilson's numerical renormalization group in imaginary time: the virtual (bond) direction as an RG axis, with physical sites as impurities in the transfer matrix [1509.01522]. MPS compresses short-scale entanglement, with the transfer-matrix fixed points characterizing effective long-range behavior.

This RG picture generalizes to MERA and hybrid MPO schemes, and illuminates the scaling of errors and the nature of excitations—perturbing only the upper MPO layers of a layered MPS restricts variational access to low-energy excitations [1509.01522].

Key limitations of the MPS approximation arise at criticality: polynomial decay of Schmidt values implies that the bond dimension must grow polynomially (or faster) in inverse error to keep local errors small [1711.06559]. In two or higher dimensions, the area law generalizes to more complex entanglement structures and the advantages of 1D MPS are lost.

## 6. Advanced Extensions: MPS for Field Theories and CFTs

MPS-based approximations extend to 2D conformal field theories (CFTs) [1601.00470, 1509.07414]. For rational chiral CFTs, the vacuum $n$-point functions can be expressed exactly as MPS correlation functions of transfer operators derived from vertex operator algebras, with truncation in mode number $N$ giving error $\mathcal{O}(q^{N/4})$ and total bond dimension polynomial in $1/\epsilon$. Similar constructions yield finitely correlated state approximations for full (chiral and anti-chiral) theories, with rigorous error bounds and explicit scaling of bond dimension versus ultraviolet cutoff and target accuracy.

## 7. Excited-State and Dynamic MPS Approximation

Beyond ground states, the MPS approximation framework generalizes to excited states via tangent-space parametrizations, random-phase approximations (RPA-MPS), and post-DMRG methods (TDA-MPS, CISD-MPS, CCSD-MPS) [1305.1761, 1103.2155]. The tangent-space formalism allows systematic exploration of low-lying excitations, construction of non-redundant bases, and direct analogues of the Thouless theorem for MPS. Time evolution and linear response within the MPS manifold can be realized by projecting the equations of motion, with dynamic error control at each time step.

---

In summary, the MPS approximation is a quantitatively controlled, theoretically rigorous, and algorithmically efficient framework for representing and simulating quantum many-body states in one dimension, with robust generalizations to infinite and continuous systems, noninteracting Gaussian states, critical models, and quantum information protocols. Its approximation properties are determined by the entanglement structure of the target state and the scaling of bond dimension with error and subsystem size, enabling high-precision studies across condensed matter, field theory, and quantum computation [1711.06559, 1310.4118, 2111.09101, 1509.01522, 2502.16464, 2511.07212, 1903.10241, 1601.00470, 1504.07701, 2512.03333, 1509.07414, 2001.11882, 1611.03779, 1103.2155, 1305.1761, 1008.4667].

Source: https://www.emergentmind.com/topics/matrix-product-state-mps-approximation