Dynamics Of Open Quantum Systems Research Articles

General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation.

We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima-Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.

General framework for quantifying dissipation pathways in open quantum systems. II. Numerical validation and the role of non-Markovianity.

In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1-214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics.

Decoherence rate in random Lindblad dynamics

Open quantum systems undergo decoherence, which is responsible for the transition from quantum to classical behavior. The time scale in which decoherence takes place can be analyzed using upper limits to its rate. We examine the dynamics of open chaotic quantum systems governed by random Lindblad operators sourced from Gaussian and Ginibre ensembles with Wigner-Dyson symmetry classes. In these systems, the ensemble-averaged purity decays monotonically as a function of time. This decay is governed by the decoherence rate, which is upper-bounded by the dimension of their Hilbert space and is independent of the ensemble symmetry. These findings hold upon mixing different ensembles, indicating the universal character of the decoherence rate limit. Moreover, our findings reveal that open chaotic quantum systems governed by random Lindbladians tend to exhibit the most rapid decoherence, regardless of the initial state. This phenomenon is associated with the concentration of the decoherence rate near its upper bound. Our work identifies primary features of decoherence in dissipative quantum chaos, with applications ranging from quantum foundations to high-energy physics and quantum technologies. Published by the American Physical Society 2024

Entangled two-photon quantum heat engine: Dissipative nonequilibrium dynamics and correlated statistics

Recently, the thermodynamics of open quantum systems driven by light fields has been investigated in the framework of quantum heat engines, whose output work can be measured as a spectroscopic signal. In this work, we investigate a four-level quantum heat engine that generates cascaded entangled photon pairs, treating the hot bath as an incoherent thermal pump and focusing on the correlation statistics of the output work and photon indistinguishability. We show that the dissipative nonequilibrium dynamics of this thermodynamic open quantum system can be reconstructed by correlation measurements under carefully mediated optical control. Comparing our findings with the traditional coherently pumped model, we find that the thermal pumping has the potential to generate nonclassical correlations resulting in photon indistinguishability, displaying an advantage in probing the nonequilibrium effects induced by the baths at different temperatures. Our work also demonstrates that incoherent pumping can optimize such correlations in output work beyond the classical limit. Lastly, we show that nonclassical correlations and resonant power at steady state cannot be simultaneously maximized in the weak coupling regime. Published by the American Physical Society 2024

Bexcitonics: Quasiparticle approach to open quantum dynamics.

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.

Observation of Non-Hermitian Edge Burst Effect in One-Dimensional Photonic Quantum Walk.

Non-Hermitian systems can exhibit unique quantum phases without any Hermitian counterparts. For example, the latest theoretical studies predict a new surprising phenomenon that bulk bands can localize and dissipate prominently at the system boundary, which is dubbed the non-Hermitian edge burst effect. Here we realize a one-dimensional non-Hermitian Su-Schrieffer-Heeger lattice with bulk translation symmetry implemented with a photonic quantum walk. Employing time-resolved single-photon detection to characterize the chiral motion and boundary localization of bulk bands, we determine experimentally that the dynamics underlying the non-Hermitian edge burst effect is due to the interplay of non-Hermitian skin effect and imaginary band gap closing. This new non-Hermitian physical effect deepens our understanding of quantum dynamics in open quantum systems.

Open Quantum System Dynamics from Infinite Tensor Network Contraction.

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.

Distillation of quantum non-Markovianity

Non-Markovianty of open quantum systems dynamics is a physically relevant property which is usually associated with the backflow of (quantum) information. Using this paradigmatic marker, we develop an operational framework to investigate how non-Markovianity for qubit dynamics can be distilled when many copies of the channels are used, possibly allowing for a stronger effect on the backflow of information.

Entanglement Assisted Probe of the Non-Markovian to Markovian Transition in Open Quantum System Dynamics.

We utilize a superconducting qubit processor to experimentally probe non-Markovian dynamics of an entangled qubit pair. We prepare an entangled state between two qubits and monitor the evolution of entanglement over time as one of the qubits interacts with a small quantum environment consisting of an auxiliary transmon qubit coupled to its readout cavity. We observe the collapse and revival of the entanglement as a signature of quantum memory effects in the environment. We then engineer the non-Markovianity of the environment by populating its readout cavity with thermal photons to show a transition from non-Markovian to Markovian dynamics, ultimately reaching a regime where the quantum Zeno effect creates a decoherence-free subspace that effectively stabilizes the entanglement between the qubits.

Transient temperature dynamics of reservoirs connected through an open quantum system.

The dynamics of open quantum systems connected with several reservoirs attract great attention due to their importance in quantum optics, biology, quantum thermodynamics, transport phenomena, etc. In many problems, the Born approximation is applicable, which implies that the influence of the open quantum system on the reservoirs can be neglected. However, in the case of long-time dynamics or mesoscopic reservoirs, the reverse influence can be crucial. In this paper, we investigate the transient dynamics of several bosonic reservoirs connected through an open quantum system. We use an adiabatic approach to study the temporal dynamics of temperatures of the reservoirs during relaxation to thermodynamic equilibrium. We show that there are various types of temperature dynamics that strongly depend on the values of dissipative rates and initial temperatures. We demonstrate that temperatures of the reservoirs, including the hottest and coldest ones, can exhibit nonmonotonic behavior. Moreover, there are moments of time during which the reservoir with an initially intermediate temperature becomes the hottest or coldest reservoir. The obtained results pave the way for managing energy flows in mesoscale and nanoscale systems.

Non-Markovian dynamics of time-fractional open quantum systems

The Time-Fractional Schrödinger Equation (TFSE) is suitable to describe the non-Markovian dynamics of a quantum system exposed to an environment, which is instructive for understanding and characterizing the time behavior of actual physical systems. Three popular TFSEs, namely Naber’s TFSE I, Naber’s TFSE II, and XGF’s TFSE, have been introduced in the form of ∂β∂tβ with β∈0,1. However, they suffer from the drawbacks. In their respective descriptions, the total probability for finding a particle in the entire space equals one only when β=1, implying that time-fractional quantum mechanics violates quantum mechanical probability conservation. By applying the three TFSEs to a basic single-qubit open system model, we discover that in describing the non-Markovian evolution of the system, the three TFSEs are limited by their own applicable ranges of β. This indicates that the three TFSEs cannot effectively describe the non-Markovian quantum dynamics. To address these issues, we introduce a well-performed TFSE by constructing a new analytic continuation of time without iβ combined with the conformable fractional derivative. In its description, for one thing, the total probability for finding the system in any single-qubit state equals one when β∈0,1. For another, the system evolves correctly in the non-Markovian manner at all values of β. Furthermore, we study the performances of the four TFSEs applying to a two-qubit open system model and show that our TFSE still possesses the above two advantages compared with the other three TFSEs.

Simple and General Unitarity Conserving Numerical Real-Time Propagators of the Time-Dependent Schrödinger Equation Based on Magnus Expansion.

Magnus expansion (ME) provides a general way to expand the real-time propagator of a time-dependent Hamiltonian within the exponential such that the unitarity is satisfied at any order. We use this property and explicit integration of Lagrange interpolation formulas for the time-dependent Hamiltonian within each time interval and derive approximations that preserve unitarity for the differential time evolution operators of general time-dependent Hamiltonians. The resulting second-order approximation is the same as using the average of Hamiltonians for two end points of time. We identify three fourth-order approximations involving commutators of Hamiltonians at different times and also derive a sixth-order expression. A test of these approximations along with other available expressions for a two-state time-dependent Hamiltonian with sinusoidal time dependences provides information on the relative performance of these approximations and suggests that the derived expressions can serve as useful numerical tools for time evolution in time-resolved spectroscopy, quantum control, quantum sensing, real-time ab initio quantum dynamics, and open system quantum dynamics.

Variational Quantum Simulation of Lindblad Dynamics via Quantum State Diffusion.

Quantum simulation of dynamics in open quantum systems is crucial but poses a significant challenge due to the non-Hermitian nature leading to nonunitary evolution and the limited quantum resources on current quantum computers. Here we introduce a variational hybrid quantum-classical algorithm designed for simulating the time evolution governed by the Lindblad master equation. Our approach involves on a stochastic unveiling of the density matrix, transforming the Lindblad equation into a wave function-based quantum state diffusion (QSD) method with the aim of reducing qubit requirements. We then apply variational quantum simulation (VQS) to efficiently capture the nonunitary evolution in QSD. We demonstrate our QSD-VQS algorithm by investigating the quantum dynamics in a two-level system subjected to an amplitude damping channel and a four-level transverse field Ising model within a dissipative environment including time-independent and periodic Hamiltonian cases. The results reveal its promising utility with upcoming hardware in the near future.

Dynamics of a strongly coupled quantum heat engine-Computing bath observables from the hierarchy of pure states.

We present a fully quantum dynamical treatment of a quantum heat engine and its baths based on the Hierarchy of Pure States (HOPS), an exact and general method for open quantum system dynamics. We show how the change of the bath energy and the interaction energy can be determined within HOPS for arbitrary coupling strength and smooth time dependence of the modulation protocol. The dynamics of all energetic contributions during the operation can be carefully examined both in its initial transient phase and, also later, in its periodic steady state. A quantum Otto engine with a qubit as an inherently nonlinear work medium is studied in a regime where the energy associated with the interaction Hamiltonian plays an important role for the global energy balance and, thus, must not be neglected when calculating its power and efficiency. We confirm that the work required to drive the coupling with the baths sensitively depends on the speed of the modulation protocol. Remarkably, departing from the conventional scheme of well-separated phases by allowing for temporal overlap, we discover that one can even gain energy from the modulation of bath interactions. We visualize these various work contributions using the analog of state change diagrams of thermodynamic cycles. We offer a concise, full presentation of HOPS with its extension to bath observables, as it serves as a universal tool for the numerically exact description of general quantum dynamical (thermodynamic) scenarios far from the weak-coupling limit.

Criteria for Davies irreducibility of Markovian quantum dynamics

The dynamics of Markovian open quantum systems are described by Lindblad master equations, generating a quantum dynamical semigroup. An important concept for such systems is (Davies) irreducibility, i.e. the question whether there exist non-trivial invariant subspaces. Steady states of irreducible systems are unique and faithful, i.e. they have full rank. In the 1970s, Frigerio showed that a system is irreducible if the Lindblad operators span a self-adjoint set with trivial commutant. We discuss a more general and powerful algebraic criterion, showing that a system is irreducible if and only if the multiplicative algebra generated by the Lindblad operators La and the operator iH+∑aLa†La , involving the Hamiltonian H, is the entire operator space. Examples for two-level systems, show that a change of Hamiltonian terms as well as the addition or removal of dissipators can render a reducible system irreducible and vice versa. Examples for many-body systems show that a large class of spin chains can be rendered irreducible by dissipators on just one or two sites. Additionally, we discuss the decisive differences between (Davies) reducibility and Evans reducibility for quantum channels and dynamical semigroups which has lead to some confusion in the recent physics literature, especially, in the context of boundary-driven systems. We give a criterion for quantum reducibility in terms of associated classical Markov processes and, lastly, discuss the relation of the main result to the stabilization of pure states and argue that systems with local Lindblad operators cannot stabilize pure Fermi-sea states.

Non-Markovian dynamics of open quantum systems via auxiliary particles with exact operator constraint

Non-Markovian dynamics of open quantum systems via auxiliary particles with exact operator constraint

Adaptive variational simulation for open quantum systems

Emerging quantum hardware provides new possibilities for quantum simulation. While much of the research has focused on simulating closed quantum systems, the real-world quantum systems are mostly open. Therefore, it is essential to develop quantum algorithms that can effectively simulate open quantum systems. Here we present an adaptive variational quantum algorithm for simulating open quantum system dynamics described by the Lindblad equation. The algorithm is designed to build resource-efficient ansatze through the dynamical addition of operators by maintaining the simulation accuracy. We validate the effectiveness of our algorithm on both noiseless simulators and IBM quantum processors and observe good quantitative and qualitative agreement with the exact solution. We also investigate the scaling of the required resources with system size and accuracy and find polynomial behavior. Our results demonstrate that near-future quantum processors are capable of simulating open quantum systems.

The influence of Ohmic noise on the dynamics of three-spin open quantum system

The dynamics of open quantum systems under decoherence effects remain a hotly debated topic in the case of the practical deployment of quantum information processing. In this regard, we examine how the three-qubit mixed state is affected by changes in the cut-off frequency of the Ohmic spectral density in a harmonic reservoir. Specifically, the degree of quantum memory-assisted entropic uncertainty (QMA-EU), along with the entanglement, coherence, and purity of the system, is to be demonstrated. A thorough exercise is done to investigate any underlying relationship between the three-qubit quantum characteristics. We show that a harmonic reservoir controlled by Ohmic noise prevails a monotonic-like decay in the current case where, sooner or later, the state becomes completely disentangled, decoherent, and mixed. The QMA-EU has always been found to have an increasing function causing the quantum resourcefulness to be reduced. Although, there is no pathway to avoid the Ohmic noise consequences and complete decay, however, we provided parameterization which would lead to prolonged preservation of quantum correlations with time. Finally, we provide various settings for the tuning of cut-off frequency regarding the Ohmic type bath and state parameters on the initial as well as final levels of quantum features.

Digital quantum simulation of non-perturbative dynamics of open systems with orthogonal polynomials

Classical non-perturbative simulations of open quantum systems' dynamics face several scalability problems, namely, exponential scaling of the computational effort as a function of either the time length of the simulation or the size of the open system. In this work, we propose the use of the Time Evolving Density operator with Orthogonal Polynomials Algorithm (TEDOPA) on a quantum computer, which we term as Quantum TEDOPA (Q-TEDOPA), to simulate non-perturbative dynamics of open quantum systems linearly coupled to a bosonic environment (continuous phonon bath). By performing a change of basis of the Hamiltonian, the TEDOPA yields a chain of harmonic oscillators with only local nearest-neighbour interactions, making this algorithm suitable for implementation on quantum devices with limited qubit connectivity such as superconducting quantum processors. We analyse in detail the implementation of the TEDOPA on a quantum device and show that exponential scalings of computational resources can potentially be avoided for time-evolution simulations of the systems considered in this work. We applied the proposed method to the simulation of the exciton transport between two light-harvesting molecules in the regime of moderate coupling strength to a non-Markovian harmonic oscillator environment on an IBMQ device. Applications of the Q-TEDOPA span problems which can not be solved by perturbation techniques belonging to different areas, such as the dynamics of quantum biological systems and strongly correlated condensed matter systems.

Path integral quantum algorithm for simulating non-Markovian quantum dynamics in open quantum systems

Path integral quantum algorithm for simulating non-Markovian quantum dynamics in open quantum systems