--- title: Random Matrix Product States (RMPS) url: https://www.emergentmind.com/topics/random-matrix-product-states-rmps type: topic --- # Random Matrix Product States (RMPS) Random Matrix Product States (RMPS) are ensembles of quantum many-body wavefunctions constructed by sequentially applying Haar-random unitaries to a set of local sites and ancillary degrees of freedom, leading to tensor network states that efficiently capture core statistical properties of fully random states. RMPS provide a physically motivated, computationally efficient framework whose statistical and entanglement features closely mimic Haar-random vectors in Hilbert space while requiring only polynomially many parameters. RMPS play a foundational role in the study of quantum statistical mechanics, quantum information, and the simulation of chaotic or thermalizing quantum systems, and serve as a controlled interpolation between product states, low-entanglement ground states, and genuinely random pure states. ## 1. Formal Definition and Ensemble Construction An RMPS is defined for a 1D quantum chain of $N$ sites, each of local dimension $d$, with maximal bond dimension $D$ (or $\chi$). The canonical construction proceeds as follows: - Attach a $D$-dimensional ancillary space to the system. - For each site $k=1,\dots,N$, apply a Haar-random unitary $U^{[k]}$ on the site and the ancilla. The output is mapped to a collection of MPS site tensors $A^{i}[k]_{\alpha\beta} = \langle i, \alpha | U^{[k]} | 0, \beta \rangle$, with physical index $i$ and ancilla indices $\alpha,\beta$. - The RMPS wavefunction is then \[ |\psi\rangle = \sum_{i_1,\dots,i_N} \langle \phi_F | A^{i_N}[N] \cdots A^{i_1}[1] | \phi_I \rangle\, |i_1 \dots i_N\rangle \] where $|\phi_I\rangle$, $|\phi_F\rangle$ are fixed boundary ancilla vectors. Sampling each $U^{[k]}$ independently from the Haar measure induces the "non-homogeneous" RMPS ensemble. The homogeneous ensemble fixes $U^{[k]} = U$ for all $k$. Alternative constructions utilize right-normalized or left-canonical MPS forms or block-unitary encodings, embedding each local tensor into a larger Haar-random unitary and projecting along fixed ancilla directions [1003.5253, 2404.18751, 2505.00073]. ## 2. Statistical and Typicality Properties ### Ensemble Average and Concentration The key feature of RMPS is that, for any such ensemble, \[ \overline{\psi} = \int_{\mathrm{Haar}} |\psi(U)\rangle\langle\psi(U)|\, dU = \frac{I}{d^N} \] independent of the ancilla dimension $D$ or $χ$ [1003.5253, 1302.3980]. This establishes that all first-moment (mean) expectation values of observables in RMPS exactly reproduce the Haar (microcanonical) ensemble. Concentration-of-measure techniques show that for any observable $O$ with support on $\ell$ contiguous sites, the variance of $\langle\psi|O|\psi\rangle$ over the RMPS ensemble is exponentially suppressed in $D$ and the local dimension. Explicitly, \[ \Pr[|\langle\psi|O|\psi\rangle - \overline{f}| \geq \epsilon] \leq c_1 \exp(-c_2\,\epsilon^2\, dD/N^2) \] for suitable constants $c_1, c_2$ and subsystem size $\ell \ll N$ [1003.5253, 1302.3980]. For subsystem reduced density matrices $\rho_S$ of size $L$, similar concentration yields \[ \Pr[\|\rho_S - \overline{\rho_S}\|_1 \geq \delta ] \leq 4^L c_1 \exp(-c_2\,\epsilon^2\, dD/N^2) \] with $\delta$ scaling as $4^{3L/2} \epsilon$ [1003.5253]. ### Entanglement Structure and Maximal Entropy Small subsystems (blocks) of a large RMPS are nearly maximally mixed provided the bond dimension $D$ satisfies $D \gg d^\ell$, where $\ell$ is the block size. The normalized reduced state $\sigma_\ell$ satisfies \[ \|\sigma_\ell - \tfrac{1}{d^\ell}\mathbb{I}\|_\infty \leq (d^\ell-1)\sqrt{d^\ell}\,O(D^{-1/10}) \] and the von Neumann entropy of the block approaches $\ell \log d$ up to corrections vanishing in $D$ [1201.6324, 2103.02634]. There is a volume-law for 2-Rényi entropy in sufficiently disconnected subsystems and almost-maximal entropy for small, connected blocks [2103.02634]. ### Higher Moments and Magic For large $D$, RMPS reproduce not just the first moment (uniform average) but can approximate higher-moment properties of Haar states. Specifically, RMPS ensembles are approximate 2-designs with corrections of $O(1/D)$ in the operator norm for $\mathrm{E}[|\psi\rangle\langle\psi|^{\otimes2}]$ [1309.0851, 2404.18751]. The so-called "magic" (nonstabilizerness) of RMPS, measured via the $L_1$-norm over the Pauli coefficients of $|\psi\rangle$, is extensive with system size $n$: $\log_d M(\psi)$ scales linearly in $n$ with overwhelming probability. Numerical and analytic results confirm that for qubits and moderate $D$, RMPS are far more magical than product or Clifford states, achieving magic approaching the maximum entropic bound [2211.10350]. ## 3. Sampling, Computational Methods, and Bias ### Sequential and Unbiased Sampling The canonical procedure for generating RMPS involves product sampling of local Haar-random unitaries (the "sequential RMPS"). However, this measure is not uniform with respect to the restriction of the Fubini–Study measure on the full Hilbert space to the MPS manifold. The sequential method exhibits an entanglement asymmetry under inversion, with marginal distributions of Schmidt eigenvalues distinct for different chain cuts [2505.00073]. An unbiased measure can be constructed using the left-canonical form and associated gauge-fixed coordinates. The correct Fubini–Study measure is then \[ d\mu_\mathrm{FS}(\psi) = d\mu_\mathrm{RMPS}(\psi)\, \prod_{\ell=1}^N |\Gamma_\ell|^{D(d-1)}, \] where $\Gamma_\ell$ is the right-environment matrix at bond $\ell$. Numerical and analytic studies confirm that the unbiased ensemble eliminates entanglement asymmetry and modifies the entanglement spectrum compared to the Haar-RMPS case [2505.00073]. ### Metropolis–Hastings Algorithm Sampling from the unbiased ensemble is achieved via a Metropolis–Hastings scheme, which proposes updates that preserve left-canonical form, with acceptance probability determined by the changes in the measure's weighting factor $w(\psi) = \prod_\ell |\Gamma_\ell|^{D(d-1)}$ [2505.00073]. ### Monte Carlo and Observables Monte Carlo sampling of RMPS can be used to compute thermal or microcanonical expectation values. Typical observables (e.g., magnetization, correlations) are estimated by evaluating expectation values in sampled RMPS and averaging. For each trial, with $N$ sites and bond dimension $\chi$, the computational cost per sample is $O(N\chi^3)$ and $M = O(1/(\chi\epsilon^2\delta))$ samples suffice for $(\epsilon, \delta)$ accuracy, where the error scales as $1/\chi$ [1309.0851]. ## 4. Applications and Physical Implications ### Quantum Statistical Mechanics RMPS efficiently model foundational aspects of statistical physics, including the emergence of microcanonical, canonical, and generalized canonical ensembles from pure states with low circuit complexity. Subsystem thermalization and typicality, normally attributed to highly entangled Haar-random states, are reproduced within the RMPS framework by entanglement and large bond dimension alone [1003.5253, 2103.02634, 1201.6324]. ### Simulation and Computation RMPS enable efficient stochastic simulation of thermal and out-of-equilibrium quantum dynamics, as only polynomial resources in $N, D$ are required to reproduce subsystem properties and expectation values that would conventionally require exponential resources for Haar sampling [1302.3980, 1309.0851]. RMPS have proven utility in tensor network algorithms and variational ansätze in ground-state and thermal-state approximations. ### Quantum Information and Magic The extensive magic of RMPS implies that they provide a rich resource for quantum computation beyond the stabilizer (Clifford) regime; RMPS can evade expressivity obstacles seen in variational states limited to Clifford or low-magic subspaces [2211.10350, 2404.18751]. The action of global Clifford unitaries on RMPS realizes an ensemble of Clifford-enhanced MPS (C MPS) that exactly forms a 3-design and approximates 4-designs to error $O(N/\chi^2)$ [2404.18751]. This demonstrates that moderate $\chi$ and Clifford post-processing suffice to approximate high-order Haar randomness and magic, relevant for randomized benchmarking and classical simulation tasks. ### Gravity and Topological Expansions RMPS have been leveraged as exact finite-$N$ models for gravitationally prepared states in 2D quantum gravity path integrals. The transfer-matrix formalism and spectral properties of RMPS directly reproduce key gravitational phenomena such as the bra–ket wormhole phase transition and allow explicit computation of higher topology and off-shell wormhole contributions to all orders [2512.11966]. Continuous RMPS models (cRMPS) generalize these results to the continuum limit. ## 5. Limitations, Extensions, and Open Directions ### Entanglement Scaling and Design Properties While RMPS are approximate 2-designs for sufficiently large $D$, they generally cannot reproduce all properties of Haar-random states at finite $D$. The 2-stabilizer Rényi entropy of RMPS converges to the Haar value as $O(N/\chi^2)$, and Clifford enhancement is required to reach approximate 4-designs [2404.18751]. For finite $D$, edge effects and lack of uniform measure in the sequentially generated ensemble may manifest, although unbiased sampling corrects these issues [2505.00073]. ### Connection to Tensor Networks Beyond 1D Current theory and numerics focus primarily on 1D chains. Extensions to higher-dimensional tensor network states, such as random PEPS or MERA, remain an active research area with preliminary evidence suggesting analogous typicality phenomena but different scaling of fluctuations and correlations [2211.10350]. ### Algorithmic and Experimental Implementation Practical sampling of RMPS is efficient for moderate $N$ and $D$, but the construction of unbiased ensembles and scaling to higher dimensions involve algorithmic overhead. Nevertheless, sequential application of Haar-random gates is implementable in quantum circuits, motivating experimental exploration, especially for many-body thermalization and benchmarking tasks [2103.02634]. ## 6. Representative Results and Quantitative Benchmarks Table: Scaling Properties and Key Results for RMPS | Property | Scaling / Value | Reference | |---------------------------------|---------------------------------------------------------------------|-------------------| | Mean state $\overline{\psi}$ | $I/d^N$ | [1003.5253] | | Typical variance of observables | $\exp(-c D)$ (for $D$ large) | [1003.5253] | | 2-design approximation error | $O(1/D)$ | [1309.0851] | | Purity for block size $\ell$ | $1/d^\ell + O(D^{-1/5})$ | [1201.6324] | | Magic growth | $\exp(\alpha n)$ with $\alpha \approx 0.5$--$1$ for small $D$ | [2211.10350] | | Convergence of subsystem entropy| $S \geq \ell \log d - O(D^{-1/10})$ | [1201.6324] | | 4-design error with Clifford | $O(N/\chi^2)$ | [2404.18751] | | MC sample cost per state | $O(N\chi^3)$ | [1309.0851] | | Metropolis unbiased sampling | Accepts with weighting $\prod_\ell |\Gamma_\ell|^{D(d-1)}$ | [2505.00073] | These results collectively establish that RMPS provide a tractable, flexible, and mathematically rigorous testbed for studying quantum typicality, entanglement, quantum resource theory, and emergent thermodynamics in low- and intermediate-entanglement quantum systems. Source: https://www.emergentmind.com/topics/random-matrix-product-states-rmps