Locally Purified Density Operators (LPDOs)

• LPDOs are tensor-network representations for mixed quantum states that ensure exact local positivity by using a purification ansatz.
• They bridge methods from matrix product states and matrix product density operators, enabling efficient simulation of open systems and noisy quantum circuits.
• Despite their physical validity through enforced Hermiticity and positivity, LPDOs face challenges like non-uniqueness and increasing bond dimensions under noise.

Searching arXiv for papers on locally purified density operators and related MPDO/MPρ literature. Locally Purified Density Operators (LPDOs) are tensor-network representations of mixed quantum states in which a density operator is expressed as a local purification and recovered by tracing out sitewise ancillary, purification, or Kraus degrees of freedom. In one dimension, this places LPDOs at the intersection of matrix product state (MPS) methods for pure states and matrix product density operator (MPDO) methods for mixed states. Their defining feature is positivity by construction: rather than representing a density matrix as a generic operator-valued tensor network, one represents a purification-like object $X$ or $|\Psi\rangle$ and forms $\rho = X X^\dagger$ or $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$. This makes LPDOs especially relevant for open-system simulation, noisy quantum circuits, tomography from local measurements, and the structural study of mixed-state phases. At the same time, LPDOs are not simply equivalent to generic MPDOs: exact local positivity can be intrinsically more expensive than operator-level representation, and the non-uniqueness of purification introduces substantial gauge structure and optimization challenges (Cuevas et al., 2013, Guo et al., 2023).

1. Definition and relation to MPDOs

An LPDO represents a mixed state through a local purification ansatz. In a standard one-dimensional form, one starts from an MPS-like purified state with one physical and one ancillary index per site,

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$

where $\tau_j$ is the physical index, $\kappa_j$ the local Kraus or ancilla index, and $\mu_j$ the virtual bond index. Tracing out all Kraus indices yields

$$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$

Equivalently, the representation can be viewed as a local factorization

$\rho = X X^\dagger,$

with $|\Psi\rangle$0 encoded as a tensor network. Because the density operator is quadratic in the local tensors, the LPDO is “guaranteed to be Hermitian and semidefinite positive by design” (Guo et al., 2023).

This construction is closely related to MPDOs. A generic MPDO is an operator written in matrix-product form,

$|\Psi\rangle$1

with doubled physical indices but without positivity built in. An LPDO induces an MPDO after contraction of the purification legs; typically the induced MPDO bond dimension satisfies $|\Psi\rangle$2 when $|\Psi\rangle$3 denotes the LPDO virtual bond dimension (Guo et al., 2023). The converse is subtler. The exact MPDO and exact local-purification descriptions are inequivalent: there is no universal bound $|\Psi\rangle$4 that controls purification rank $|\Psi\rangle$5 purely as a function of MPDO operator Schmidt rank $|\Psi\rangle$6 (Cuevas et al., 2013).

This distinction is central. MPDOs provide a compact operator-level language, but positivity is not locally manifest. LPDOs provide local positivity and Hermiticity by construction, but the purification bond dimensions and local ancilla structure can be substantially more demanding. In modern terminology, papers using the labels MP$|\Psi\rangle$7, CP-MPO, locally purified form, LPTN, or positive MPDO are often describing essentially the same underlying class, with differences mainly in gauge conventions, algorithmic emphasis, and physical interpretation of the extra local index (Jamadagni et al., 2024, Müller et al., 2024).

2. Structural complexity and the cost of local positivity

A foundational result for LPDOs is that exact local positivity is not complexity-neutral. For a mixed state $|\Psi\rangle$8, one may compare an MPDO representation with operator Schmidt rank

$|\Psi\rangle$9

to a purification-MPS representation with purification rank

$\rho = X X^\dagger$0

Although every purification yields an MPDO with $\rho = X X^\dagger$1, the reverse implication fails in general (Cuevas et al., 2013).

The separation is already visible on classical diagonal states. For a family

$\rho = X X^\dagger$2

the paper identifies

$\rho = X X^\dagger$3

Choosing $\rho = X X^\dagger$4 as the slack matrix of a regular $\rho = X X^\dagger$5-gon gives

$\rho = X X^\dagger$6

while

$\rho = X X^\dagger$7

Thus the MPDO bond dimension can remain constant while the purification rank diverges (Cuevas et al., 2013). The implication is precise: no exact mixed-state tensor-network description can be both universally efficient and locally positivity-certifying.

The same work also provides constructive remedies. The sum-of-squares polynomial method builds purifications through

$\rho = X X^\dagger$8

yielding approximate or exact LPDOs with purification-rank bounds

$\rho = X X^\dagger$9

Its approximate version is formulated as a semidefinite program. A second route, the eigenbasis method, starts from the standard purification

$$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$0

and gives exact and approximate purification-rank bounds scaling with the rank or truncated rank of $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$1 (Cuevas et al., 2013). The overarching conclusion is not that LPDOs fail, but that they should be regarded primarily as approximate tools unless additional spectral or entanglement structure is present.

A complementary theoretical result concerns mixed-state area laws. One-dimensional mixed states with sufficiently small $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$2-Rényi entanglement of purification can be approximated efficiently by MPDOs in trace norm with polynomial bond dimension and super-polynomially decaying error (Jarkovsky et al., 2020). The theorem is stated in MPDO language, but its controlling quantity,

$$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$3

is defined by minimizing over purifications. This strongly suggests that LPDO-type descriptions are natural in the same low-$$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$4 regime, even though the paper does not prove existence of a single globally compatible local purification with the same scaling (Jarkovsky et al., 2020).

3. Gauge freedom, canonicalization, and representation non-uniqueness

LPDOs are highly non-unique. As in ordinary MPS, there is gauge freedom on virtual bonds, but LPDOs also possess additional freedom on the local purification legs. If $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$5 is an isometry acting on the $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$6-subspace, then one can transform one factorization $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$7 into another $$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$8 without changing the physical density operator (Jamadagni et al., 19 Sep 2025). Similarly, in the MP$$\rho=\operatorname{Tr}_{\kappa}|\Psi\rangle\langle\Psi|$$9 formulation,

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$0

so internal transformations on the $$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$1-index leave $$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$2 invariant (Jamadagni et al., 2024).

This non-uniqueness has direct computational consequences. An LPDO can retain large virtual bond dimensions even when the physical state becomes simple. A particularly sharp example is the maximally mixed state. The optimal LPDO for qubits has

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$3

since

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$4

Yet if one reaches the maximally mixed state by applying noise to an initially entangled pure-state LPDO, the resulting representation can still carry large inherited $$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$5-bonds despite unit fidelity with the known maximally mixed operator (Jamadagni et al., 19 Sep 2025). The discrepancy is representation-theoretic rather than physical.

To manage this, several works develop mixed-state canonical forms. In the MP$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$6 framework, one reshapes local tensors $$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$7 and performs SVDs,

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$8

leading to a mixed-state isometric condition

$$\ket{\psi} = \sum_{\{\bm{\tau},\bm{\kappa}\}} \sum_{\{\bm{\mu}\}} \prod_{j=1}^N [A_j]^{\tau_j,\kappa_j}_{\mu_{j-1},\mu_j} \ket{\tau_1,\kappa_1,\cdots,\tau_N,\kappa_N},$$9

The orthogonality-center concept from MPS is generalized so that coherent quantum correlations are carried by $\tau_j$0-bonds and positive classical-mixture correlations by $\tau_j$1-structure (Jamadagni et al., 2024).

A related proposal for LPDO compression in noisy-circuit settings uses a unified treatment of virtual and Kraus bonds. Rather than truncating Kraus legs first and then virtual bonds, one first performs QR and LQ decompositions to construct left and right gauges $\tau_j$2, then builds projectors for both Kraus and virtual indices, and finally applies all projectors simultaneously. For Kraus compression, one absorbs the gauges into a local tensor $\tau_j$3, performs

$\tau_j$4

and sets

$\tau_j$5

For virtual compression, one SVDs the product $\tau_j$6 and defines projectors

$\tau_j$7

The point is not merely numerical convenience; it reflects the fact that virtual and Kraus resources jointly determine LPDO expressivity (Guo et al., 2023).

4. Dynamics, noisy circuits, and open-system simulation

LPDOs are particularly natural for Markovian open-system dynamics because CPTP channels act directly on the purification operator. In one formulation, the state is written as

$\tau_j$8

with $\tau_j$9 encoded by local rank-4 tensors $\kappa_j$0,

$\kappa_j$1

Applying a local channel enlarges the local Kraus dimensions but preserves positivity by construction (Godinez-Ramirez et al., 2024).

A major application is noisy quantum-circuit simulation. One influential study maps an LPDO on $\kappa_j$2 qubits to a pure-state supervector on a $\kappa_j$3 ladder and argues that efficient LPDO representation is controlled by the entanglement structure of this ladder state (Guo et al., 2023). For brick-wall circuits with Haar-random two-qubit gates and local noise, the paper derives entanglement-growth bounds such as

$\kappa_j$4

so that if

$\kappa_j$5

an area law for the ladder supervector is expected and efficient LPDO simulation should be possible (Guo et al., 2023).

The same work identifies three dynamical regimes: a quantum region, a classical region, and a difficult quantum-classical critical point. LPDOs represent states well in both the quantum and classical regions as sets of states, but actual time evolution tends to fail at the crossover. A projection experiment from MPO to LPDO shows that classical-region states may still be representable, yet an LPDO trajectory propagated from the beginning often cannot traverse the critical region accurately enough to reach them (Guo et al., 2023). This is a structural limitation of the positivity-preserving manifold, not merely a bond-dimension issue.

A related line of work studies positive tensor networks for deep noisy-circuit simulation and uses the term MPDO for a locally purified density matrix. The central bottleneck is that each application of noise increases local purification dimensions, leading otherwise to exponential growth. The proposed remedy is to exploit purification gauge freedom through DMRG-like sweeps of local two-qubit basis optimization. The local ancilla unitary $\kappa_j$6 is optimized to minimize purification entanglement,

$\kappa_j$7

where $\kappa_j$8 is a Rényi entropy of the singular values across the virtual bond (Müller et al., 2024). The striking empirical result is that disentangling the purified state also reduces local purification dimensions. In noisy random circuits, this keeps both $\kappa_j$9 and $\mu_j$0 bounded with depth while maintaining reasonable truncation fidelity (Müller et al., 2024).

For continuous-time open-system evolution, another work treats purely dissipative Lindbladian dynamics with nearest-neighbor jump operators,

$\mu_j$1

and improves standard second-order Trotter-Suzuki splitting by optimizing Kraus gauges on a Stiefel manifold (Godinez-Ramirez et al., 2024). Each local two-site channel is converted to a Kraus-isometry representation

$\mu_j$2

and the global product of local channels is optimized to approximate the exact propagator,

$\mu_j$3

The method preserves positivity, improves splitting error by orders of magnitude in favorable cases, and also serves as a compression mechanism by controlling the effective local Choi rank $\mu_j$4 (Godinez-Ramirez et al., 2024).

5. Tomography, local data, and reconstruction from measurements

LPDOs provide a natural variational family for mixed-state tomography because positivity is enforced by the ansatz rather than by post hoc constraints. A tomography framework based on local measurements reconstructs an unknown state by fitting local reduced density matrices of an LPDO ansatz (Guo et al., 2023). For contiguous windows $\mu_j$5, the local loss is

$\mu_j$6

with each term

$\mu_j$7

The gradient is

$\mu_j$8

and tensors are updated locally via

$\mu_j$9

with Adam (Guo et al., 2023).

The reconstruction uses informationally complete local Pauli measurements. For a block of length $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$0,

$$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$1

Although there are $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$2 local Pauli observables in principle, the number of distinct global measurement settings can be reduced to

$$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$3

by using translation-invariant basis configurations (Guo et al., 2023).

This LPDO tomography scheme has been demonstrated numerically for one-dimensional pure and mixed states and for two-dimensional pure states up to $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$4, and experimentally on IBM and Quafu hardware (Guo et al., 2023). One-dimensional mixed-state tests under depolarizing, bit-flip, amplitude-damping, and phase-damping noise show that the method works well, but increasing Kraus dimension alone is insufficient when only short-range local measurements are used. Increasing the local measurement range $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$5 is substantially more effective; with $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$6, reconstruction fidelity above $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$7 was reported for depolarizing noise up to $$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$8 (Guo et al., 2023).

A limitation emerges for locally indistinguishable states such as GHZ-like families. With only local loss terms, the relative phase remains undetermined, so additional global observables must be added to the objective. This is not specific to LPDOs, but it highlights the difference between local positivity and global identifiability (Guo et al., 2023).

6. Phases, locality, and dissipative parent structures

Beyond numerical representation, LPDOs have become a language for mixed-state structure and phases. One direction studies whether mixed-state fixed points can be realized as steady states of local open-system dynamics. For MPDO renormalization fixed points (RFPs), local, frustration-free parent Lindbladians have been constructed analytically, with minimal steady-state degeneracy and, in some classes, necessarily noncommuting local terms (Liu et al., 17 Jan 2025). The work is phrased in MPDO language rather than LPDO language, but the implications carry over whenever the MPDO admits a local purification. The parent Lindbladian acts on the physical mixed state itself; the paper explicitly notes that preparing a purification and tracing ancillas is not an adequate substitute, both because not every MPDO has a local purification and because partial trace is not itself generated by a Lindbladian (Liu et al., 17 Jan 2025).

Another foundational direction investigates when locally purified MPDO-type states are approximately Gibbs states of quasi-local Hamiltonians. A class of MPDOs built from chains of 1-input/2-output “Y-shaped” completely positive maps is studied as a local-purification/Stinespring construction (Chen et al., 2020). The locality diagnostic is exponential decay of the quantum conditional mutual information

$$\rho = \mathrm{Tr}_{\kappa}\!\left[\ket{\psi}\!\bra{\psi}\right].$$9

For bistochastic channels with trivial correctable algebra, the paper proves

$\rho = X X^\dagger,$0

with $\rho = X X^\dagger,$1 iff the correctable algebra is trivial (Chen et al., 2020). It also proves exponential decay for strictly positive channels and channels with a forgetful component. This strongly supports the view that generic locally purified mixed states are approximately thermal states of quasi-local Hamiltonians, although the most general formulation remains conjectural (Chen et al., 2020).

LPDOs have also entered the classification of mixed-state symmetry-protected topological phases. In $\rho = X X^\dagger,$2 and $\rho = X X^\dagger,$3, LPDOs are used to distinguish weak injectivity and strong injectivity, corresponding respectively to decay of ordinary linear correlators and decay of Rényi-2 correlators in the doubled-state formalism (Guo et al., 2024). For $\rho = X X^\dagger,$4, weak injectivity is defined through injectivity of the purified MPS tensor

$\rho = X X^\dagger,$5

while strong injectivity is defined through injectivity of the doubled tensor

$\rho = X X^\dagger,$6

This produces a mixed-state analogue of the MPS/PEPS classification framework, with virtual symmetry actions on both virtual and Kraus spaces and explicit fixed-point LPDO constructions for decohered average SPT phases (Guo et al., 2024).

A different but related structural development concerns local-to-global compatibility of density operators under Markov constraints. A recent study of the marginal problem for density operators introduces a canonical logarithmic reconstruction

$\rho = X X^\dagger,$7

for clique marginals on a chordal graph and proves that the trace-one condition

$\rho = X X^\dagger,$8

is equivalent to existence of a unique quantum Markov completion, which is also the maximum-entropy completion (Lauritzen et al., 19 May 2026). Although this is not an LPDO paper, it is directly relevant to LPDO-style reconstruction from local data: it characterizes when local marginals support a low-complexity global mixed state with a prescribed Markov structure and when quantum noncommutativity obstructs such gluing.

7. Advantages, limitations, and current outlook

LPDOs offer three persistent advantages. First, they preserve positivity and Hermiticity exactly through the ansatz. Second, they are naturally adapted to local channels and open-system dynamics, because Kraus indices can be absorbed into local purification structure. Third, they permit a physically meaningful separation between coherent virtual correlations and mixture correlations stored in purification degrees of freedom (Guo et al., 2023, Jamadagni et al., 2024).

These same features create the main difficulties. LPDOs are less expressive than generic MPOs at fixed effective compression, and this gap is particularly pronounced near dynamical crossovers in noisy circuits (Guo et al., 2023). Exact LPDO complexity cannot be controlled solely by MPDO bond dimension (Cuevas et al., 2013). The representation is highly non-unique, so practical bond dimensions may reflect gauge artifacts rather than intrinsic state complexity (Jamadagni et al., 19 Sep 2025). Under repeated noisy evolution, local Kraus dimensions tend to grow rapidly, which historically limited scalability until gauge-optimization and purification-disentangling methods were introduced (Müller et al., 2024, Godinez-Ramirez et al., 2024).

The comparison with MPOs remains nuanced. MPOs can often track entanglement trends across longer noisy-circuit evolutions, but truncation may violate positivity and produce unphysical density operators. LPDOs maintain physicality exactly, which is essential for tasks where the resulting state must remain interpretable as a bona fide density matrix (Guo et al., 2023). The tradeoff is therefore not between accuracy and inaccuracy simpliciter, but between unconstrained operator compression and positivity-preserving compression.

Taken together, the literature suggests a stable division of labor. LPDOs are especially well suited to regimes where positivity, CPTP compatibility, and weak-to-moderate mixed-state complexity matter more than maximal operator-space expressivity. This includes low-noise open-system simulation, physically valid tensor-network tomography, fixed-point and phase-structure analysis, and settings where local purification is itself the relevant organizing principle. MPOs remain useful as more flexible operator ansätze, but LPDOs define the positivity-preserving branch of mixed-state tensor-network theory, with an increasingly rich body of results on their approximation properties, gauge structure, dissipative dynamics, and role in mixed-state many-body physics (Jarkovsky et al., 2020, Guo et al., 2023, Müller et al., 2024, Liu et al., 17 Jan 2025).