Experimental Simulation of Open Quantum System Dynamics via Trotterization.

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

The rapid development of quantum computers has enabled demonstrations of quantum advantages on various tasks. However, real quantum systems are always dissipative due to their inevitable interaction with the environment, and the resulting non-unitary dynamics make quantum simulation challenging with only unitary quantum gates. In this work, we present an innovative and scalable method to simulate open quantum systems using quantum computers. We define an adjoint density matrix as a counterpart of the true density matrix, which reduces to a mixed-unitary quantum channel and thus can be effectively sampled using quantum computers. This method has several benefits, including no need for auxiliary qubits and noteworthy scalability. Moreover, some long-time properties like steady states and the thermal equilibrium can also be investigated as the adjoint density matrix and the true dissipated one converge to the same state. Finally, we present deployments of this theory in the dissipative quantum XY model for the evolution of correlation and entropy with short-time dynamics and the disordered Heisenberg model for many-body localization with long-time dynamics. This work promotes the study of real-world many-body dynamics with quantum computers, highlighting the potential to demonstrate practical quantum advantages.

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.

Quantum error control and mitigation techniques help improve how quantum computers handle errors, making algorithms run more efficiently despite noisy hardware. These strategies work at the software level and are built into programs before they run on quantum machines. Unlike quantum error correction, which actively detects and fixes errors during computation, error mitigation does not use real-time corrections but instead reduces the impact of errors after computations. Error mitigation can also help reduce the extra resources needed for full error correction. This chapter explains how these techniques relate to each other, provides an overview of key error mitigation methods and their limitations, and highlights leading academic and technological players in this field. It also covers current trends, such as integrating error-aware programming, benchmarking performance, and ensuring access to quantum hardware for testing. To advance quantum computing, it is recommended that research ecosystems support specialized centers and contribute to key strategic areas in the field.

Simulating spin-boson dynamics with stochastic Liouville–von Neumann equations

We present a framework to simulate the dynamics of hard probes such as heavy quarks or jets in a hot, strongly coupled quark-gluon plasma (QGP) on a quantum computer. Hard probes in the QGP can be treated as open quantum systems governed in the Markovian limit by the Lindblad equation. However, due to large computational costs, most current phenomenological calculations of hard probes evolving in the QGP use semiclassical approximations of the quantum evolution. Quantum computation can mitigate these costs and offers the potential for a fully quantum treatment with exponential speed-up over classical techniques. We report a simplified demonstration of our framework on IBM Q quantum devices and apply the random identity insertion method to account for cnot depolarization noise, in addition to measurement error mitigation. Our work demonstrates the feasibility of simulating open quantum systems on current and near-term quantum devices, which is of broad relevance to applications in nuclear physics, quantum information, and other fields.

Quantum simulation of dynamics is an important goal in the noisy intermediate-scale quantum era, within which quantum error mitigation may be a viable path towards modifying or eliminating the effects of noise. Most studies on quantum error mitigation have focused on the resource cost due to its exponential scaling in the circuit depth. Methods such as probabilistic error cancellation rely on discretizing the evolution into finite time steps and applying the mitigation layer after each time step, modifying only the noise part without any Hamiltonian dependence. This may lead to Trotter-like errors in the simulation results even if the error mitigation is implemented ideally, which means that the number of samples is taken as infinite. Here we analyze the aforementioned errors which have been largely neglected before. We show that they are determined by the commutating relations between the superoperators of the unitary part, the device noise part, and the noise part of the open dynamics to be simulated. We include both digital quantum simulation and analog quantum simulation setups and consider defining the ideal error mitigation map both by exactly inverting the noise channel and by approximating it to first order in the time step. We use single-qubit toy models to numerically demonstrate our findings. Our results illustrate fundamental limitations of applying probabilistic error cancellation in a stepwise manner to continuous dynamics, thus motivating the investigations of truly time-continuous error cancellation methods. Published by the American Physical Society 2024

We propose a method that combines the quantum-classical mapping approach to surface hopping with the dissipative quantum dynamics of the Lindblad master equation. Like conventional surface-hopping methods, our approach is based on classical trajectories coupled to the dynamics of a quantum subsystem. However, instead of evolving the subsystem wavefunction according to the time-dependent Schrödinger equation, we use stochastic quantum trajectories derived from secular Redfield theory. This approach enables the simulation of open quantum systems coupled simultaneously to Markovian quantum baths and anharmonic, non-Markovian classical degrees of freedom. Applications to the spin-boson model and to the cavity-enhanced fluorescence of an electronically nonadiabatic molecule show excellent agreement with fully quantum-mechanical benchmarks.

Accurate simulation of the time evolution of a quantum system under the influence of an environment is critical to making accurate predictions in chemistry, condensed-matter physics, and materials sciences. Whereas there has been a recent surge in interest in quantum algorithms for the prediction of nonunitary time evolution in quantum systems, few studies offer a direct quantum analog to the Lindblad equation. Here, we present a quantum algorithm—utilizing a decomposition of nonunitary operators approach—that models dynamic processes via the unraveled Lindblad equation. This algorithm is employed to probe both a two-level system in an amplitude damping channel as well as the transverse field Ising model in a variety of parameter regimes; the resulting population dynamics demonstrate excellent agreement with classical simulation, showing the promise of predicting population dynamics utilizing quantum devices for a variety of important systems in molecular energy transport, quantum optics, and other open quantum systems.

Entanglement phase transitions in hybrid quantum circuits describe individual quantum trajectories rather than the measurement-averaged ensemble, despite the fact that results of measurements are not conventionally used for feedback. Here, we numerically demonstrate that a class of generalized measurements with identical measurement-averaged dynamics give different phases and phase transitions. We show that measurement-averaged destruction of Bell state entanglement is a useful proxy for determining which hybrid circuit yields the lowest-entanglement dynamics. We use this to argue that no unfolding of our model can avoid a volume law phase, which has implications for simulation of open quantum systems.

Methodological aspects of using the driven Liouville-von Neumann (DLvN) approach for simulating dynamical properties of molecular junctions are discussed. As a model system we consider a non-interacting resonant level uniformly coupled to a single Fermionic bath. We demonstrate how a finite system can mimic the depopulation dynamics of the dot into an infinite band bath of continuous and uniform density of states. We further show how the effects of spurious energy resolved currents, appearing due to the approximate nature of the equilibrium state obtained in DLvN calculations, can be avoided. Several ways to approach the wide band limit, which is often adopted in analytical treatments, using a finite numerical model system are discussed including brute-force increase of the lead model bandwidth as well as efficient cancellation or direct subtraction of finite-bandwidth effect. These methodological considerations may be relevant also for other numerical schemes that aim to study non-equilibrium thermodynamics via simulations of open quantum systems.

We consider the quantum-control task of obtaining a target unitary operation via control fields that couple to the quantum system and are chosen to best mitigate errors resulting from time-dependent noise. We allow for two sources of noise: fluctuations in the control fields and those arising from the environment. We address the issue of error mitigation by means of a formulation rooted in the Martin-Siggia-Rose (MSR) approach to noisy, classical statistical-mechanical systems. We express the noisy control problem in terms of a path integral, and integrate out the noise to arrive at an effective, noise-free description. We characterize the degree of success in error mitigation via a fidelity, which characterizes the proximity of the sought-after evolution to ones achievable in the presence of noise. Error mitigation is then accomplished by applying the optimal control fields, i.e., those that maximize the fidelity subject to any constraints obeyed by the control fields. To make connection with MSR, we reformulate the fidelity in terms of a Schwinger-Keldysh (SK) path integral, with the added twist that the `forward' and `backward' branches of the time-contour are inequivalent with respect to the noise. The present approach naturally allows the incorporation of constraints on the control fields; a useful feature in practice, given that they feature in real experiments. We illustrate this MSR-SK approach by considering a system consisting of a single spin $s$ freedom (with $s$ arbitrary), focusing on the case of $1/f$ noise. We discover that optimal error-mitigation is accomplished via a universal control field protocol that is valid for all $s$, from the qubit (i.e., $s=1/2$) case to the classical (i.e., $s \to \infty$) limit. In principle, this MSR-SK approach provides a framework for addressing quantum control in the presence of noise for systems of arbitrary complexity.

The hierarchical equations of motion (HEOM), derived from the exact Feynman-Vernon path integral, is one of the most powerful numerical methods to simulate the dynamics of open quantum systems. Its applicability has so far been limited to specific forms of spectral reservoir distributions and relatively elevated temperatures. Here we solve this problem and introduce an effective treatment of quantum noise in frequency space by systematically clustering higher order Matsubara poles, equivalent to an optimized rational decomposition. This leads to an elegant extension of the HEOM to arbitrary temperatures and very general reservoirs in combination with efficiency, high accuracy, and long-time stability. Moreover, the technique can directly be implemented in other approaches such as Green's function, stochastic, and pseudomode formulations. As one highly nontrivial application, for the subohmic spin-boson model at vanishing temperature the Shiba relation is quantitatively verified which predicts the long-time decay of correlation functions.

Quantum simulation is a powerful tool to study the properties of quantum systems. The dynamics of open quantum systems are often described by completely positive (CP) maps, for which several quantum simulation schemes exist. Such maps, however, represent only a subset of a larger class of maps: the general dynamical maps which are linear, Hermitian preserving, and trace preserving but not necessarily positivity preserving. Here we present a simulation scheme for these general dynamical maps, which occur when the underlying system-reservoir model undergoes entangling (and thus non-Markovian) dynamics. Such maps also arise as the inverse of CP maps, which are commonly used in error mitigation. We illustrate our simulation scheme on an IBM quantum processor, demonstrating its ability to recover the initial state of a Lindblad evolution. This paves the way for a novel form of quantum error mitigation. Our scheme only requires one ancilla qubit as an overhead and a small number of one and two qubit gates. Consequently, we expect it to be of practical use in near-term quantum devices.

We inquire into the time evolution of quantum systems associated with pseudo-or quasi-Hermitian Hamiltonians. We obtain, in the pseudo-Hermitian case, a generalized Liouville-von Neumann equation for closed systems. We show that quantum systems with quasi-Hermitian Hamiltonians admit the proper interpretation in terms of open quantum system and derive a generalized Lindblad-Kossakowski equation. Finally, we extend such formalism to the study of decaying systems.