Open quantum systems (OQS's) are ubiquitous in non-equilibrium quantum dynamics and in quantum science and technology. Solving the dynamics of an OQS in a quantum many-body bath has been considered a computationally hard problem because of the dimensionality curse. Here, considering that full knowledge of the bath dynamics is unnecessary for describing the reduced dynamics of an OQS, we prove a polynomial complexity theorem, that is, the number of independent equations required to fully describe the dynamics of an OQS increases at most linearly with the evolution time and polynomially with the bath size. Therefore, efficient computational algorithms exist for solving the dynamics of a small-sized OQS (such as a qubit or an atom). We further prove that, when the dynamics of an OQS and the bath is represented by a tensor network, a tensor contraction procedure can be specified such that the bond dimension (i.e., the range of tensor indices contracted in each step) increases only linearly (rather than exponentially) with the evolution time, providing explicitly efficient algorithms for a wide range of OQS's. We demonstrate the theorems and the tensor-network algorithm by solving two widely encountered OQS problems, namely, a spin in a Gaussian bath (the spin-boson model) and a central spin coupled to many environmental spins (the Gaudin model). This work provides approaches to understanding dynamics of OQS's, learning the environments via quantum sensors, and optimizing quantum information processing in noisy environments.
Open quantum systems (OQS's) are ubiquitous in non-equilibrium quantum dynamics and in quantum science and technology. Solving the dynamics of an OQS in a quantum many-body bath has been considered a computationally hard problem because of the dimensionality curse. Here, considering that full knowledge of the bath dynamics is unnecessary for describing the reduced dynamics of an OQS, we prove a polynomial complexity theorem, that is, the number of independent equations required to fully describe the dynamics of an OQS increases at most linearly with the evolution time and polynomially with the bath size. Therefore, efficient computational algorithms exist for solving the dynamics of a small-sized OQS (such as a qubit or an atom). We further prove that, when the dynamics of an OQS and the bath is represented by a tensor network, a tensor contraction procedure can be specified such that the bond dimension (i.e., the range of tensor indices contracted in each step) increases only linearly (rather than exponentially) with the evolution time, providing explicitly efficient algorithms for a wide range of OQS's. We demonstrate the theorems and the tensor-network algorithm by solving two widely encountered OQS problems, namely, a spin in a Gaussian bath (the spin-boson model) and a central spin coupled to many environmental spins (the Gaudin model). This work provides approaches to understanding dynamics of OQS's, learning the environments via quantum sensors, and optimizing quantum information processing in noisy environments.
Approaching the long-time dynamics of non-Markovian open quantum systems presents a challenging task if the bath is strongly coupled. Recent proposals address this problem through a representation of the so-called process tensor in terms of a tensor network. We show that for Gaussian environments highly efficient contraction to a matrix product operator (MPO) form can be achieved with infinite MPO evolution methods, leading to significant computational speed-up over existing proposals. The result structurally resembles open system evolution with carefully designed auxiliary degrees of freedom, as in hierarchical or pseudomode methods. Here, however, these degrees of freedom are generated automatically by the MPO evolution algorithm. Moreover, the semigroup form of the resulting propagator enables us to explore steady-state physics, such as phase transitions.
The difficulty to simulate the dynamics of open quantum systems resides in their coupling to many-body reservoirs with exponentially large Hilbert space. Applying a tensor network approach in the time domain, we demonstrate that effective small reservoirs can be defined and used for modeling open quantum dynamics. The key element of our technique is the timeline reservoir network (TRN), which contains all the information on the reservoir's characteristics, in particular, the memory effects timescale. The TRN has a one-dimensional tensor network structure, which can be effectively approximated in full analogy with the matrix product approximation of spin-chain states. We derive the sufficient bond dimension in the approximated TRN with a reduced set of physical parameters: coupling strength, reservoir correlation time, minimal timescale, and the system's number of degrees of freedom interacting with the environment. The bond dimension can be viewed as a measure of the open dynamics complexity. Simulation is based on the semigroup dynamics of the system and effective reservoir of finite dimension. We provide an illustrative example showing the scope for new numerical and machine learning-based methods for open quantum systems.
We present an efficient, nearly optimal quantum algorithm for solving linear matrix differential equations, with applications to the simulation of open quantum systems and beyond. For unitary or dissipative dynamics, the algorithm computes an entry of the solution matrix with query complexity O(νLt/ε), where the constant ν depends on the problem parameters, L involves a time integral of upper bounds on the norms of evolution operators, and ε is the error. In particular, νL is linear in t for unitary dynamics and can be a constant for dissipative dynamics. Our result contrasts prior quantum approaches for differential equations that typically require exponential time for this problem due to the encoding in a quantum state, which can lead to exponentially small amplitudes. We demonstrate the utility of the algorithm through an end-to-end application, namely the simulation of dissipative dynamics for non-interacting fermions, which can be extended to other quantum and classical systems. We compare with classical algorithms and give evidence of polynomial quantum speedups for systems in a lattice, which become more pronounced for systems with long-range interactions and can be shown to be exponential in general. We also provide a lower bound of Ω(νLt/ε) for unitary or dissipative dynamics that proves our algorithm is optimal up to logarithmic factors.
Open quantum systems are ubiquitous in the physical sciences, with widespread applications in the areas of chemistry, condensed-matter physics, material science, optics, and many more.Not surprisingly, there is significant interest in their efficient simulation.However, direct classical simulation quickly becomes intractable with coupling to an environment with an effective dimension that grows exponentially.This raises the question: Can quantum computers help model these complex dynamics?A first step in answering this question requires an understanding of the computational complexity of this task.Here, we map the temporal complexity of a process to the spatial complexity of a many-body state using a computational model known as the process-tensor framework.With this, we are able to explore the simulation complexity of an open quantum system as a dynamic sampling problem: a system coupled to an environment can be probed at successive points in time-accessing multitime correlations.The complexity of multitime sampling, which is an important and interesting problem in its own right, contains the complexity of master equations and stochastic maps as a special case.Our results show how the complexity of the underlying quantum stochastic process corresponds to the complexity of the associated family of master equations for the dynamics.We present both analytical and numerical examples where the multitime sampling of an open quantum system is as complex as sampling from a many-body state that is classically hard.This also implies that the corresponding family of master equations is classically hard.Our results pave the way for studying open quantum systems from a complexity-theoretic perspective, highlighting the role quantum computers will play in our understanding of quantum dynamics.
The study of interacting quantum many-body systems is one of the central challenges in modern physics. Interestingly, by understanding the physical properties of such systems from a microscopic description, one can explain many complex quantum phenomena such as the high-temperature superconductivity, fractional quantum Hall effect, etc. However, due to high complexity of these systems, theoretical investigations such as analytical and perturbation theories are limited to studying only a small subset of problems. In the past few decades, numerical simulation has emerged as an indispensable method for studying many phenomena, and has led to immense success and progress in recent years. The main obstacle in addressing quantum many-body systems is that the number of parameters required to represent a state of the system scales exponentially with the system size. However, one hardly needs the full Hilbert space but only a small subspace to describe physical properties of the system. From the point of view in quantum information theory, quantum entanglement is one of the main factors that is responsible for the complexity of studying quantum many-body systems. In particular, the amount of quantum entanglement in such system quantifies the correlation between constituents and therefore it determines how big the subspace is required to represent the system. In this thesis, we are interested in using tensor network states, a new emerging numerical technique to study dynamical and ground state properties of the system. The key ingredient is that a tensor network state which consists of a fairly small number of parameters can represent very well the complete wave function of a system. In one spatial dimension, Matrix Product State (MPS) has become a standard tool, combined with related algorithms such as Time Evolving Block Decimation algorithm (TEBD) to study dynamical properties of a locally-perturbed quantum lattice system. In this thesis, we introduce a technique called ``infinite boundary conditions'' that allows us to use a finite-size MPS to represent the state of an infinite-size system. As a result, one can simulate dynamics of a system for a longer time with a significant computational cost reduction compared to traditional methods. In two spatial dimensions, the tensor network known as Projected Entangled Pair States (PEPS) is emerging as a useful approach. However algorithms for optimizing the matrix elements of these tensors are complex and not very efficient. In particular, the efficiency of the optimization depends critically on the quality of the so-called ``environment'', which is expensive to compute. We have introduced a new method to accelerate the convergence of calculating the PEPS ground state by recycling the environment tensors for many optimization steps, vastly reducing the number of times that the full environment needs to be calculated.
In this paper we present a quantum algorithm that uses noise as a resource. The goal of our quantum algorithm is the calculation of operator averages of an open quantum system evolving in time. Selected low-noise system qubits and noisy bath qubits represent the system and the bath of the open quantum system. All incoherent qubit noise can be mapped to bath spectral functions. The form of the spectral functions can be tuned digitally, allowing for the time evolution of a wide range of open-system models at finite temperature. We study the feasibility of this approach with a focus on the solution of the spin-boson model and assume intrinsic qubit noise that is dominated by damping and dephasing. We find that classes of open quantum systems exist where our algorithm performs very well, even with gate errors as high as 1%. In general the presented algorithm performs best if the system-bath interactions can be decomposed into native gates.