Matrix product density operators (MPDOs) are tensor network representations of locally purified density matrices where each physical degree of freedom is associated to an environment degree of freedom. MPDOs have interesting properties for mixed state representations: guaranteed positivity by construction, efficient conservation of the trace and computation of local observables. However, they have been challenging to use for noisy quantum circuit simulation, as the application of noise increases the dimension of the environment Hilbert space, leading to an exponential growth of bond dimensions. MPDOs also lack a unique canonical form, due to the freedom in the choice of basis for the environment Hilbert space, which leads to a vast variation of bond dimensions. In this work, we present a systematic way to reduce the bond dimensions of MPDOs by disentangling the purified state. We optimize the basis for the environment Hilbert space by performing density matrix renormalization group (DMRG)-like sweeps of local 2-qubit basis optimization. Interestingly, we find that targeting only the disentanglement of the purified state leads to a reduction of the environment dimension. In other words, a compact MPDO representation requires a low-entanglement purified state. We apply our compression method to the emulation of noisy random quantum circuits. Our technique allows us to keep bounded bond dimensions, and thus bounded memory, contrary to previous works on MPDOs, while keeping reasonable truncation fidelities.
Matrix product density operators (MPDOs) are tensor network representations of locally purified density matrices where each physical degree of freedom is associated to an environment degree of freedom. MPDOs have interesting properties for mixed state representations: guaranteed positivity by construction, efficient conservation of the trace and computation of local observables. However, they have been challenging to use for noisy quantum circuit simulation, as the application of noise increases the dimension of the environment Hilbert space, leading to an exponential growth of bond dimensions. MPDOs also lack a unique canonical form, due to the freedom in the choice of basis for the environment Hilbert space, which leads to a vast variation of bond dimensions. In this work, we present a systematic way to reduce the bond dimensions of MPDOs by disentangling the purified state. We optimize the basis for the environment Hilbert space by performing density matrix renormalization group (DMRG)-like sweeps of local 2-qubit basis optimization. Interestingly, we find that targeting only the disentanglement of the purified state leads to a reduction of the environment dimension. In other words, a compact MPDO representation requires a low-entanglement purified state. We apply our compression method to the emulation of noisy random quantum circuits. Our technique allows us to keep bounded bond dimensions, and thus bounded memory, contrary to previous works on MPDOs, while keeping reasonable truncation fidelities.
Abstract Open quantum system simulations are essential for exploring novel quantum phenomena and evaluating noisy quantum circuits. In this Letter, we investigate whether mixed states generated from noisy quantum circuits can be efficiently represented by locally purified density operators (LPDOs). We map an LPDO of N qubits to a pure state of size 2 × N defined on a ladder and introduce a unified method for managing virtual and Kraus bonds. We numerically simulate noisy random quantum circuits with depths of up to d = 40 using fidelity and entanglement entropy as accuracy measures. The LPDO representation is effective in describing mixed states in both the quantum and classical regions; however, it encounters significant challenges at the quantum-classical critical point, restricting its applicability to the quantum region. In contrast, matrix product operators (MPO) successfully characterize the entanglement trend throughout the simulation, while the truncation in MPOs breaks the positivity condition required for a physical density matrix. This work advances our understanding of efficient mixed-state representations in open quantum systems and provides insights into the entanglement structure of noisy quantum circuits.
Simulating quantum circuits with classical computers requires resources growing exponentially in terms of system size.Real quantum computer with noise, however, may be simulated polynomially with various methods considering different noise models.In this work, we simulate random quantum circuits in one dimension with matrix product density operators (MPDOs), for different noise models such as dephasing, depolarizing, and amplitude damping.We show that the method based on matrix product states (MPSs) fails to approximate the noisy output quantum states for any of the noise models considered, while the MPDO method approximates them well.Compared with the method of matrix product operators (MPOs), the MPDO method reflects a clear physical picture of noise (with inner indices taking care of the noise simulation) and quantum entanglement (with bond indices taking care of two-qubit gate simulation).Consequently, in case of weak system noise, the resource cost of the MPDO will be significantly less than that of the MPO due to a relatively small inner dimension needed for the simulation.In case of strong system noise, a relatively small bond dimension may be sufficient to simulate the noisy circuits, indicating a regime that the noise is large enough for an "easy" classical simulation, which is further supported by a comparison with the experimental results on an IBM cloud device.Moreover, we propose a more effective tensor updates scheme with optimal truncations for both the inner and the bond dimensions, performed after each layer of the circuit, which enjoys a canonical form of the MPDO for improving simulation accuracy.With truncated inner dimension to a maximum value κ and bond dimension to a maximum value χ , the cost of our simulation scales as ∼NDκ 3 χ 3 , for an N-qubit circuit with depth D.
We aim to rigorously address how local noise affects classical simulability of quantum dynamics benchmarked by tensor-network methods. Using operator entanglement entropy (OEE), we prove the following: (1) For single-qubit depolarizing noise on arbitrary circuits, tensor networks with poly(n) bond dimension suffice for fixed absolute Hilbert-Schmidt error after $\order{1}$ depth, while relative error demands $\order{\log n}$ depth; and this bound is optimal. (2) For single-qubit depolarizing noise on 1D local circuits, the existence of whole-trajectory error-bounded matrix product operator (MPO) of poly(n) bond dimension at all depths. (3) For general single-qubit noise in 1D brickwall circuits, random two-design gates with contraction coefficient c<1/3 yield an $\order{1}$ OEE plateau with probability 1−Te−Ω(n), while arbitrary gates with c<1/48 give $\order{\log n}$ OEE in the worst case. (4) In higher dimensions, these bounds yield uniform-in-depth poly(n) average boundary-bond dimensions for projected entangled pair operators~(PEPO) across every cut -- under depolarizing noise at either absolute or relative accuracy, and under general noise with strong contraction at absolute accuracy. Our results establish a rigorous connection between certain noise models, circuit types, and their classical simulability.
In recent years, efficient quantum circuit simulations incorporating ideal noise assumptions have relied on tensor network simulators, particularly leveraging the matrix product density operator (MPDO) framework. However, experiments on real noisy intermediate-scale quantum (NISQ) devices often involve complex noise profiles, encompassing uncontrollable elements and instrument-specific effects such as crosstalk. To address these challenges, we employ quantum process tomography (QPT) techniques to directly capture the operational characteristics of the experimental setup and integrate them into numerical simulations using MPDOs. Our QPT-assisted MPDO simulator is then applied to explore a variational approach for generating noisy entangled states, comparing the results with standard noise numerical simulations and demonstrations conducted on the Quafu cloud quantum computation platform. Additionally, we investigate noisy MaxCut problems, as well as the effects of crosstalk and noise truncation. Our results provide valuable insights into the impact of noise on NISQ devices and lay the foundation for enhanced design and assessment of quantum algorithms in complex noise environments.
Locally Purified Density Operators (LPDOs) are state-of-the-art tensor network ansatze candidates that efficiently represent mixed quantum states at scale. However, given their non-uniqueness, their representational complexity is generally sub-optimal in practical computations. In this work we perform a comprehensive numerical and analytical analysis and resolve this issue in the experimentally relevant limit where noise depolarizes the density operator into a maximally mixed state. To resolve the sub-optimality issue, we analyze two numerical tools, one analytic method, and detail the relations between them. The numerical tools used are fidelity-preserving truncations and isometric gauge transformations leveraging Riemannian optimizations over entropic objective functions. In addition, by invoking the injectivity and symmetry constraints of the maximally mixed LPDO, we also present analytical closed-form expressions for the disentangler and discuss their relation to numerical optimizers. Further, away from the maximally mixed state, our simulations highlight how the truncation threshold smoothly interpolate, as a function of depolarization, between established matrix product results and our new results. Our work shows how, by minimizing the resources required to represent key states of practical interest in experiment, the efficiency of tensor network algorithms can be substantially increased. This paves the path for uncovering tensor network’s fundamental scalability limits and latent potential in representing the wide locus of mixed quantum states that are accessible on near-term quantum devices.
Entangled possibly mixed states are an essential resource for quantum computation, communication, metrology, and the simulation of many-body systems. It is important to develop and improve preparation protocols for such states. One possible way to prepare states of interest is to design an open system that evolves only towards the desired states. A Markovian evolution of a quantum system can be generally described by a Lindbladian. Tensor networks provide a framework to construct physically relevant entangled states. In particular, matrix product density operators (MPDOs) form an important variational class of states. MPDOs generalize matrix product states to mixed states, can represent thermal states of local one-dimensional Hamiltonians at sufficiently large temperatures, describe systems that satisfy the area law of entanglement, and form the basis of powerful numerical methods. In this work we develop an algorithm that determines for a given linear subspace of MPDOs whether this subspace can be the stable space of some frustration free k-local Lindbladian and, if so, outputs an appropriate Lindbladian. We proceed by using machine learning with networks of quantum channels, also known as quantum neural networks (QNNs), to train denoising post-processing devices for quantum sources. First, we show that QNNs can be trained on imperfect devices even when part of the training data is corrupted. Second, we show that QNNs can be trained to extrapolate quantum states to, e.g., lower temperatures. Third, we show how to denoise quantum states in an unsupervised manner. We develop a novel quantum autoencoder that successfully denoises Greenberger-Horne-Zeilinger, W, Dicke, and cluster states subject to spin-flip, dephasing errors, and random unitary noise. Finally, we develop recurrent QNNs (RQNNs) for denoising that requires memory, such as combating drifts. RQNNs can be thought of as matrix product quantum channels with a quantum algorithm for training and are closely related to MPDOs. The proposed preparation and denoising protocols can be beneficial for various emergent quantum technologies and are within reach of present-day experiments.