Taming Quantum Noise for Efficient Low Temperature Simulations of Open Quantum Systems.

We develop a quasiparticle approach to capture the dynamics of open quantum systems coupled to bosonic thermal baths of arbitrary complexity based on the Hierarchical Equations of Motion (HEOM). This is done by generalizing the HEOM dynamics and mapping it into that of the system in interaction with a few bosonic fictitious quasiparticles that we call bexcitons. Bexcitons arise from a decomposition of the bath correlation function into discrete features. Specifically, bexciton creation and annihilation couple the auxiliary density matrices in the HEOM. The approach provides a systematic strategy to construct exact quantum master equations that include the system-bath coupling to all orders even for non-Markovian environments. Specifically, by introducing different metrics and representations for the bexcitons it is possible to straightforwardly generate different variants of the HEOM, demonstrating that all these variants share a common underlying quasiparticle picture. Bexcitonic properties, while unphysical, offer a coarse-grained view of the correlated system-bath dynamics and its numerical convergence. For instance, we use it to analyze the instability of the HEOM when the bath is composed of underdamped oscillators and show that it leads to the creation of highly excited bexcitons. The bexcitonic picture can also be used to develop more efficient approaches to propagate the HEOM. As an example, we use the particle-like nature of the bexcitons to introduce mode-combination of bexcitons in both number and coordinate representation that uses the multi-configuration time-dependent Hartree to efficiently propagate the HEOM dynamics.

Stochastic evolution underpins several approaches to the dynamics of open quantum systems, such as random modulation of Hamiltonian parameters, the stochastic Schrödinger equation (SSE), and the stochastic Liouville equation (SLE). These approaches replace the explicit system–environment coupling with an effective system‐only dynamics, where dissipative behavior emerges from ensemble averaging. Stochastic Hamiltonians, in particular, have long served as phenomenological tools in physical chemistry to include environmental effects without recourse to an explicit microscopic derivation. In this work, we aim at a self‐contained and accessible presentation of these approaches to further elaborate on their common roots in essential concepts of stochastic calculus and to delineate the conditions under which they are equivalent. We also discuss how different formulations naturally lead to different numerical time‐integration schemes, better suited for either classical simulation platforms, based on finite‐difference approximations, or quantum algorithms, that employ random unitary maps. Our analysis aims at providing a unified perspective and actionable recipes for classical and quantum implementations of stochastic evolution in the simulation of open quantum systems.

Non-Markovian effects in open quantum systems are central to understanding spectral lineshape. Here, we quantify the non-Markovianity associated with both overdamped and underdamped vibrations in terms of information flow between the bath and the system and compare this with the broadening and ellipticity of two-dimensional spectra. Using the Breuer Laine Piilo (BLP) measure, we link the well-known stochastic models for spectral lineshape with modern quantum information theory. Specifically, we study the effect of non-Markovianity in a system in contact with underdamped vibrations and examine the differences observed on increasing the damping to the overdamped limit. The open quantum system dynamics are evolved using the hierarchical equations of motion, efficiently terminated with a Markovian cutoff, where separate hierarchies are derived for the underdamped and overdamped environments. It is shown that the BLP measure is quantitatively correlated with the ellipticity of two-dimensional spectra and memory effects are more pronounced in underdamped environments, due to the long-lived feedback of information between the system and its bath, compared to overdamped environments. Environmental signatures in spectral lineshapes emerge as a result of information flow from the bath back into the system.

Information scrambling refers to the phenomenon in which local quantum information in a many-body system becomes dispersed throughout the entire system under unitary evolution. It has been extensively studied in closed quantum systems, where it is quantified by operator size growth, revealing deep connections between condensed matter physics, high-energy physics, and quantum information. However, when extending the study of operator size dynamics to open quantum systems, two different definitions of operator size distributions emerge. These definitions are based on different treatments of the bath. In this work, we aim to establish a unified picture for operator size dynamics in open quantum systems, using the solvable Brownian SYK models at generic system size. In particular, we provide the conditions under which the signature of scrambling transition — discovered using one particular definition — appears in operator size dynamics under the other definition. Additionally, we extend previous studies by exploring finite-size effects that are not captured by the scramblon theory. Our results provide a refined understanding of operator size dynamics in open quantum systems.

We present a renormalization group (RG) method which allows for an analytical study of the transient dynamics of open quantum systems on all time scales. Whereas oscillation frequencies and decay rates of exponential time evolution follow from the fixed point positions, the long-time behavior of pre-exponential functions is related to the scaling behavior around the fixed points. We show that certain terms of the RG flow are only cut off by inverse time, which leads to a difference between infrared and ultraviolet scaling. An evaluation for the ohmic spin boson model at weak damping reveals significant deviations from previous predictions in the long-time regime. We propose that weak coupling problems for stationary quantities can in principle turn into strong coupling ones for the determination of the long-time behavior.

A new approach to theory and simulation of the non-Markovian dynamics of open quantum systems is presented. It is based on identification of a parameter which is uniformly bounded on wide time intervals: the occupation of the virtual cloud of quanta. By ‘virtual’ we denote those bath excitations which were emitted by the open system, but eventually will be reabsorbed before any measurement of the bath state. A useful property of the virtual cloud is that the number of its quanta is expected to saturate on long times, since physically this cloud is a (retarded) polarization of the bath around the system. Therefore, the joint state of open system and virtual cloud (we call it dressed state) can be accurately represented in a truncated basis of Fock states, on a wide time scale. At the same time, there can be an arbitrarily large number of the observable quanta (which survive up to measurement), especially if the open system is under driving. However, it turns out that the statistics of the bath-measurement outcomes is classical (in a suitable measurement basis): one can employ a Monte Carlo sampling of these outcomes. Therefore, it is possible to efficiently simulate the dynamics of the observable quantum field. In this work we consider the bath measurement with respect to the coherent states, which yields the Husimi function as the positive (quasi)probability distribution of the outcomes. The joint evolution of the dressed state and the corresponding outcome is called the dressed quantum trajectory. The Monte Carlo sampling of these trajectories yields a stochastic simulation method with promising convergence properties on wide time scales.

This paper discusses several methods for describing the dynamics of open quantum systems, where the environment of the open system is infinite-dimensional. These are purifications, phase space forms, master equation and liouville equation forms. The main contribution is in using Feynman-Kac formalisms to describe the infinite-demsional components.

Digital quantum simulators provide a diversified tool for solving the evolution of quantum systems with complicated Hamiltonians and hold great potential for a wide range of applications. Although much attention is paid to the unitary evolution of closed quantum systems, dissipation and noise are vital in understanding the dynamics of practical quantum systems. In this work, we experimentally demonstrate a digital simulation of an open quantum system in a controllable Markovian environment with the assistance of a single ancillary qubit. By Trotterizing the quantum Liouvillians, the continuous evolution of an open quantum system is effectively realized, and its application in error mitigation is demonstrated by adjusting the simulated noise intensities. High-order Trotter for open quantum dynamics is also experimentally investigated and shows higher accuracy. Our results represent a significant step toward hardware-efficient simulation of open quantum systems and error mitigation in quantum algorithms in noisy intermediate-scale quantum systems.

The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent effects of the bath, the reservoir correlation functions are usually decomposed into summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. This method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultralow temperature and general band structures. We demonstrate the efficiency, accuracy, and long-time stability of this decomposition scheme by applying it to the Anderson impurity model in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.

Simulating physical systems on quantum devices is one of the most promising applications of quantum technology. Current quantum approaches to simulating open quantum systems are still practically challenging on NISQ-era devices, because they typically require ancilla qubits and extensive controlled sequences. In this work, we propose a hybrid quantum-classical approach for simulating a class of open system dynamics called random-unitary channels. These channels naturally decompose into a series of convex unitary evolutions, which can then be efficiently sampled and run as independent circuits. The method does not require deep ancilla frameworks and thus can be implemented with lower noise costs. We implement simulations of open quantum systems up to dozens of qubits and with large channel ranks. Published by the American Physical Society 2024

Non-Markovian open quantum systems represent the most general dynamics when the quantum system is coupled with a bath environment. The quantum dynamics arising from many important applications are non-Markovian. Although for special cases, such as Hamiltonian evolution and Lindblad evolution, quantum simulation algorithms have been extensively studied, efficient quantum simulations for the dynamics of non-Markovian open quantum systems remain underexplored. The most immediate obstacle for studying such systems is the lack of a universal succinct description of their dynamics. In this work, we fulfill the gap of studying such dynamics by (1) providing a succinct representation of the dynamics of non-Markovian open quantum systems with quantifiable error, and (2) developing an efficient quantum algorithm for simulating such dynamics with cost $${\mathcal {O}}(t\, \textrm{polylog}(t/\epsilon ))$$ for evolution time t and precision $$\epsilon $$ . Our derivation of the succinct representation is based on stochastic Schrödinger equations, which could lead to new alternatives to deal with open quantum systems as well.

We present a general analysis of the role of initial correlations between the open system and an environment on quantum dynamics of the open system.

This paper proposes the concept of subsystem controllability for both closed and open quantum dynamical systems, and the existence conditions of controllable subsystems are derived accordingly. It is demonstrated that a controllable subsystem may exist for a closed quantum dynamical system even if the overall system is not controllable. Furthermore, we also demonstrate that, for a class of decoherence-free subsystems of Markovian open quantum dynamical systems, the subsystems are controllable under open-loop coherent control if the initial states of the open systems are factorized. To a certain extent, this paper gives an alternative control theoretical interpretation on why decoherence-free subsystems are useful for quantum computation.

The repeated interaction model provides a framework for emulating and analyzing the dynamics of open quantum systems. We explore here the dynamics generated by this protocol in a system that is simultaneously coupled to two baths through noncommuting system operators. One bath is made to couple to nondiagonal elements of the system, thus it induces dissipative dynamics, while the other couples to diagonal elements, and by itself it generates pure dephasing. By solving the problem analytically exactly, we show that when both baths act concurrently, a strong systembath coupling gives rise to nonadditive effects in the dynamics. A prominent signature of this nonadditivity is the characteristic slowing down of population relaxation, driven by the influence of the dephasing bath. Beyond dynamics, we investigate the thermodynamic behavior of the model. Previous studies, using quantum master equations, showed that strong system-bath coupling created bath-cooperativity in this model, allowing heat exchange to the dephasing (diagonally coupled) bath. We find instead that, under the repeated interaction scheme, heat flows exclusively to the dissipative bath (coupled through nondiagonal elements). Our results highlight the need for a deeper understanding of the types of open quantum system dynamics and steady-state phenomena that emerge within the repeated interaction framework and the relation of this protocol to other common open quantum system techniques.

A photonic circuit is generally described as a structure in which light propagates by unitary exchange and transfers reversibly between channels. In contrast, the term ‘diffusive’ is more akin to a chaotic propagation in scattering media, where light is driven out of coherence towards a thermal mixture. Based on the dynamics of open quantum systems, the combination of these two opposites can result in novel techniques for coherent light control. The crucial feature of these photonic structures is dissipative coupling between modes, via an interaction with a common reservoir. Here, we demonstrate experimentally that such systems can perform optical equalisation to smooth multimode light, or act as a distributor, guiding it into selected channels. Quantum thermodynamically, these systems can act as catalytic coherent reservoirs by performing perfect non-Landauer erasure. For lattice structures, localised stationary states can be supported in the continuum, similar to compacton-like states in conventional flat-band lattices.