Simulation Complexity of Open Quantum Dynamics: Connection with Tensor Networks

Exact Hidden Markovian Dynamics in Quantum Circuits

Characterizing nonequilibrium dynamics in quantum many-body systems is a challenging frontier of physics. In this Letter, we systematically construct solvable nonintegrable quantum circuits that exhibit exact hidden Markovian subsystem dynamics. This feature thus enables accurately calculating local observables for arbitrary evolution time. Utilizing the influence matrix method, we show that the influence of the time-evolved global system on a finite subsystem can be analytically described by sequential, time-local quantum channels acting on the subsystem with an ancilla of finite Hilbert space dimension. The realization of exact hidden Markovian property is facilitated by a solvable condition on the underlying two-site gates in the quantum circuit. We further present several concrete examples with varying local Hilbert space dimensions to demonstrate our approach.

Transitions in Computational Complexity of Continuous-Time Local Open Quantum Dynamics

We analyze the complexity of classically simulating continuous-time dynamics of locally interacting quantum spin systems with a constant rate of entanglement breaking noise. We prove that a polynomial time classical algorithm can be used to sample from the state of the spins when the rate of noise is higher than a threshold determined by the strength of the local interactions. Furthermore, by encoding a 1D fault tolerant quantum computation into the dynamics of spin systems arranged on two or higher dimensional grids, we show that for several noise channels, the problem of weakly simulating the output state of both purely Hamiltonian and purely dissipative dynamics is expected to be hard in the low-noise regime.

Multipartite Correlations in Quantum Collision Models

Quantum collision models have proved to be useful for a clear and concise description of many physical phenomena in the field of open quantum systems: thermalization, decoherence, homogenization, nonequilibrium steady state, entanglement generation, simulation of many-body dynamics, and quantum thermometry. A challenge in the standard collision model, where the system and many ancillas are all initially uncorrelated, is how to describe quantum correlations among ancillas induced by successive system-ancilla interactions. Another challenge is how to deal with initially correlated ancillas. Here we develop a tensor network formalism to address both challenges. We show that the induced correlations in the standard collision model are well captured by a matrix product state (a matrix product density operator) if the colliding particles are in pure (mixed) states. In the case of the initially correlated ancillas, we construct a general tensor diagram for the system dynamics and derive a memory-kernel master equation. Analyzing the perturbation series for the memory kernel, we go beyond the recent results concerning the leading role of two-point correlations and consider multipoint correlations (Waldenfelds cumulants) that become relevant in the higher-order stroboscopic limits. These results open an avenue for the further analysis of memory effects in collisional quantum dynamics.

Time-Dependent Variational Principle for Open Quantum Systems with Artificial Neural Networks

We develop a variational approach to simulating the dynamics of open quantum many-body systems using deep autoregressive neural networks. The parameters of a compressed representation of a mixed quantum state are adapted dynamically according to the Lindblad master equation by employing a time-dependent variational principle. We illustrate our approach by solving the dissipative quantum Heisenberg model in one dimension for up to 40 spins and in two dimensions for a 4×4 system and by applying it to the simulation of confinement dynamics in the presence of dissipation.

Constructing tensor network influence functionals for general quantum dynamics

We describe an iterative formalism to compute influence functionals that describe the general quantum dynamics of a subsystem beyond the assumption of linear coupling to a quadratic bath. We use a space-time tensor network representation of the influence functional and investigate its approximability in terms of its bond dimension and time-like entanglement in the tensor network description. We study two numerical models, the spin-boson model and a model of interacting hard-core bosons in a 1D harmonic trap. We find that the influence functional and the intermediates involved in its construction can be efficiently approximated by low bond dimension tensor networks in certain dynamical regimes, which allows the quantum dynamics to be accurately computed for longer times than with direct time evolution methods. However, as one iteratively integrates out the bath, the correlations in the influence functional can first increase before decreasing, indicating that the final compressibility of the influence functional is achieved via non-trivial cancellation.

Model-Independent Simulation Complexity of Complex Quantum Dynamics

We present a model-independent measure of dynamical complexity based on simulation of complex quantum dynamics using stroboscopic Markovian dynamics. Tools from classical signal processing enable us to infer the Hilbert space dimension of the complex quantum system evolving under a time-independent Hamiltonian via pulsed interrogation. We illustrate this using simulated third-order pump-probe spectroscopy data for exciton transport in a toy model of a coupled dimer with vibrational levels, revealing the dimension of the singly excited manifold of the dimer. Finally, we probe the complexity of excitonic transport in light harvesting 2 (LH2) and Fenna-Matthews-Olson (FMO) complexes using data from two recent nonlinear ultrafast optical spectroscopy experiments. For the latter we make model-independent inferences that are commensurate with model-specific ones, including the estimation of the fewest number of parameters needed to fit the experimental data and identifying the spatial extent, i.e., delocalization size, of quantum states participating in this complex quantum dynamics.

Quantum Information and Algorithms for Correlated Quantum Matter

Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.

Rate Operator Unraveling for Open Quantum System Dynamics

Stochastic methods with quantum jumps are often used to solve open quantum system dynamics. Moreover, they provide insight into fundamental topics, such as the role of measurements in quantum mechanics and the description of non-Markovian memory effects. However, there is no unified framework to use quantum jumps to describe open-system dynamics in any regime. We solve this issue by developing the rate operator quantum jump (ROQJ) approach. The method not only applies to both Markovian and non-Markovian evolutions, but also allows us to unravel master equations for which previous methods do not work. In addition, ROQJ yields a rigorous measurement-scheme interpretation for a wide class of dynamics, including a set of master equations with negative decay rates, and sheds light on different types of memory effects which arise when using stochastic quantum jump methods.

Machine Learning Non-Markovian Quantum Dynamics

Machine learning methods have proved to be useful for the recognition of patterns in statistical data. The measurement outcomes are intrinsically random in quantum physics, however, they do have a pattern when the measurements are performed successively on an open quantum system. This pattern is due to the system-environment interaction and contains information about the relaxation rates as well as non-Markovian memory effects. Here we develop a method to extract the information about the unknown environment from a series of projective single-shot measurements on the system (without resorting to the process tomography). The method is based on embedding the non-Markovian system dynamics into a Markovian dynamics of the system and the effective reservoir of finite dimension. The generator of Markovian embedding is learned by the maximum likelihood estimation. We verify the method by comparing its prediction with an exactly solvable non-Markovian dynamics. The developed algorithm to learn unknown quantum environments enables one to efficiently control and manipulate quantum systems.

Exploiting the Causal Tensor Network Structure of Quantum Processes to Efficiently Simulate Non-Markovian Path Integrals

In the path integral formulation of the evolution of an open quantum system coupled to a Gaussian, noninteracting environment, the dynamical contribution of the latter is encoded in an object called the influence functional. Here, we relate the influence functional to the process tensor-a more general representation of a quantum stochastic process-describing the evolution. Then, we use this connection to motivate a tensor network algorithm for the simulation of multitime correlations in open systems, building on recent work where the influence functional is represented in terms of time evolving matrix product operators. By exploiting the symmetries of the influence functional, we are able to use our algorithm to achieve orders-of-magnitude improvement in the efficiency of the resulting numerical simulation. Our improved algorithm is then applied to compute exact phonon emission spectra for the spin-boson model with strong coupling, demonstrating a significant divergence from spectra derived under commonly used assumptions of memorylessness.

Variational Autoencoder Reconstruction of Complex Many-Body Physics

Thermodynamics is a theory of principles that permits a basic description of the macroscopic properties of a rich variety of complex systems from traditional ones, such as crystalline solids, gases, liquids, and thermal machines, to more intricate systems such as living organisms and black holes to name a few. Physical quantities of interest, or equilibrium state variables, are linked together in equations of state to give information on the studied system, including phase transitions, as energy in the forms of work and heat, and/or matter are exchanged with its environment, thus generating entropy. A more accurate description requires different frameworks, namely, statistical mechanics and quantum physics to explore in depth the microscopic properties of physical systems and relate them to their macroscopic properties. These frameworks also allow to go beyond equilibrium situations. Given the notably increasing complexity of mathematical models to study realistic systems, and their coupling to their environment that constrains their dynamics, both analytical approaches and numerical methods that build on these models show limitations in scope or applicability. On the other hand, machine learning, i.e., data-driven, methods prove to be increasingly efficient for the study of complex quantum systems. Deep neural networks, in particular, have been successfully applied to many-body quantum dynamics simulations and to quantum matter phase characterization. In the present work, we show how to use a variational autoencoder (VAE)—a state-of-the-art tool in the field of deep learning for the simulation of probability distributions of complex systems. More precisely, we transform a quantum mechanical problem of many-body state reconstruction into a statistical problem, suitable for VAE, by using informationally complete positive operator-valued measure. We show, with the paradigmatic quantum Ising model in a transverse magnetic field, that the ground-state physics, such as, e.g., magnetization and other mean values of observables, of a whole class of quantum many-body systems can be reconstructed by using VAE learning of tomographic data for different parameters of the Hamiltonian, and even if the system undergoes a quantum phase transition. We also discuss challenges related to our approach as entropy calculations pose particular difficulties.