Dynamics Of Open Quantum Systems Research Articles

Symmetry-resolved dynamical purification in synthetic quantum matter

When a quantum system initialized in a product state is subjected to either coherent or incoherent dynamics, the entropy of any of its connected partitions generically increases as a function of time, signalling the inevitable spreading of (quantum) information throughout the system. Here, we show that, in the presence of continuous symmetries and under ubiquitous experimental conditions, symmetry-resolved information spreading is inhibited due to the competition of coherent and incoherent dynamics: in given quantum number sectors, entropy decreases as a function of time, signalling dynamical purification. Such dynamical purification bridges between two distinct short and intermediate time regimes, characterized by a log-volume and log-area entropy law, respectively. It is generic to symmetric quantum evolution, and as such occurs for different partition geometry and topology, and classes of (local) Liouville dynamics. We then develop a protocol to measure symmetry-resolved entropies and negativities in synthetic quantum systems based on the random unitary toolbox, and demonstrate the generality of dynamical purification using experimental data from trapped ion experiments [Brydges et al., Science 364, 260 (2019)]. Our work shows that symmetry plays a key role as a magnifying glass to characterize many-body dynamics in open quantum systems, and, in particular, in noisy-intermediate scale quantum devices.

Simulation of open quantum systems by automated compression of arbitrary environments

The central challenge for describing the dynamics in open quantum systems is that the Hilbert space of typical environments is too large to be treated exactly. In some cases, such as when the environment has a short memory time or only interacts weakly with the system, approximate descriptions of the system are possible. Beyond these, numerically exact methods exist, but these are typically restricted to baths with Gaussian correlations, such as non-interacting bosons. Here we present a numerically exact method for simulating open quantum systems with arbitrary environments which consist of a set of independent degrees of freedom. Our approach automatically reduces the large number of environmental degrees of freedom to those which are most relevant. Specifically, we show how the process tensor -- which describes the effect of the environment -- can be iteratively constructed and compressed using matrix product state techniques. We demonstrate the power of this method by applying it to problems with bosonic, fermionic, and spin environments: electron transport, phonon effects and radiative decay in quantum dots, central spin dynamics, anharmonic environments, dispersive coupling to time-dependent lossy cavity modes, and superradiance. The versatility and efficiency of our automated compression of environments (ACE) method provides a practical general-purpose tool for open quantum systems.

Decoherence-Induced Exceptional Points in a Dissipative Superconducting Qubit.

Open quantum systems interacting with an environment exhibit dynamics described by the combination of dissipation and coherent Hamiltonian evolution. Taken together, these effects are captured by a Liouvillian superoperator. The degeneracies of the (generically non-Hermitian) Liouvillian are exceptional points, which are associated with critical dynamics as the system approaches steady state. We use a superconducting transmon circuit coupled to an engineered environment to observe two different types of Liouvillian exceptional points that arise either from the interplay of energy loss and decoherence or purely due to decoherence. By dynamically tuning the Liouvillian superoperators in real time we observe a non-Hermiticity-induced chiral state transfer. Our study motivates a new look at open quantum system dynamics from the vantage of Liouvillian exceptional points, enabling applications of non-Hermitian dynamics in the understanding and control of open quantum systems.

Improving the efficiency of open-quantum-system simulations using matrix product states in the interaction picture

Modeling open quantum systems---quantum systems coupled to a bath---is of value in condensed-matter theory, cavity quantum electrodynamics, nanosciences, and biophysics. The real-time simulation of open quantum systems was advanced significantly by the recent development of chain mapping techniques and the use of matrix product states that exploit the intrinsic entanglement structure in open quantum systems. The computational cost of simulating open quantum systems, however, remains high when the bath is excited to high-lying quantum states. We develop an approach to reduce the computational costs in such cases. The interaction representation for the open quantum system is used to distribute excitations among the bath degrees of freedom so that the occupation of each bath oscillator is ensured to be low. The interaction picture also causes the matrix dimensions to be much smaller in a matrix product state of a chain-mapped open quantum system than in the Schr\"odinger picture. Using the interaction representation accelerates the calculations by one to two orders of magnitude over the existing matrix-product-state method. In the regime of strong system-bath coupling and high temperatures, the speedup can be as large as three orders of magnitude. The approach developed here is especially promising to simulate the dynamics of open quantum systems in high-temperature and strong-coupling regimes.

Open quantum system dynamics and the mean force Gibbs state

The dynamical convergence of a system to the thermal distribution, or Gibbs state, is a standard assumption across all of the physical sciences. The Gibbs state is determined just by temperature and energies of the system. However, at decreasing system sizes, i.e., for nanoscale and quantum systems, the interaction with their environments is not negligible. The question then arises: Is the system's steady state still the Gibbs state? If not, how may the steady state depend on the interaction details? Here, we provide an overview of recent progress on answering these questions. We expand on the state of the art along two general avenues: First, we take the static point-of-view, which postulates the so-called mean force Gibbs state. This view is commonly adopted in the field of strong coupling thermodynamics, where modified laws of thermodynamics and nonequilibrium fluctuation relations are established on the basis of this modified state. Second, we take the dynamical point of view, originating from the field of open quantum systems, which examines the time-asymptotic steady state within two paradigms. We describe the mathematical paradigm, which proves return to equilibrium, i.e., convergence to the mean force Gibbs state, and then discuss a number of microscopic physical methods, particularly master equations. We conclude with a summary of established links between statics and equilibration dynamics and provide an extensive list of open problems. This comprehensive overview will be of interest to researchers in the wider fields of quantum thermodynamics, open quantum systems, mesoscopic physics, statistical physics, and quantum optics and will find applications whenever energy is exchanged on the nanoscale, from quantum chemistry and biology to magnetism and nanoscale heat management.

Autoregressive Neural Network for Simulating Open Quantum Systems via a Probabilistic Formulation.

The theory of open quantum systems lays the foundation for a substantial part of modern research in quantum science and engineering. Rooted in the dimensionality of their extended Hilbert spaces, the high computational complexity of simulating open quantum systems calls for the development of strategies to approximate their dynamics. In this Letter, we present an approach for tackling open quantum system dynamics. Using an exact probabilistic formulation of quantum physics based on positive operator-valued measure, we compactly represent quantum states with autoregressive neural networks; such networks bring significant algorithmic flexibility due to efficient exact sampling and tractable density. We further introduce the concept of string states to partially restore the symmetry of the autoregressive neural network and improve the description of local correlations. Efficient algorithms have been developed to simulate the dynamics of the Liouvillian superoperator using a forward-backward trapezoid method and find the steady state via a variational formulation. Our approach is benchmarked on prototypical one-dimensional and two-dimensional systems, finding results which closely track the exact solution and achieve higher accuracy than alternative approaches based on using Markov chain MonteCarlo method to sample restricted Boltzmann machines. Our Letter provides general methods for understanding quantum dynamics in various contexts, as well as techniques for solving high-dimensional probabilistic differential equations in classical setups.

Revealing higher-order light and matter energy exchanges using quantum trajectories in ultrastrong coupling

The dynamics of open quantum systems is often modelled using master equations, which describe the expected outcome of an experiment (i.e., the average over many realizations of the same dynamics). Quantum trajectories, instead, model the outcome of ideal single experiments -- the ``clicks'' of a perfect detector due to, e.g., spontaneous emission. The correct description of quantum jumps, which are related to random events characterizing a sudden change in the wave function of an open quantum system, is pivotal to the definition of quantum trajectories. In this article, we extend the formalism of quantum trajectories to open quantum systems with ultrastrong coupling (USC) between light and matter by properly defining jump operators in this regime. In such systems, exotic higher-order quantum-state- and energy-transfer can take place without conserving the total number of excitations in the system. The emitted field of such USC systems bears signatures of these higher-order processes, and significantly differs from similar processes at lower coupling strengths. Notably, the emission statistics must be taken at a single quantum trajectory level, since the signatures of these processes are washed out by the ``averaging'' of a master equation. We analyze the impact of the chosen unravelling (i.e., how one collects the output field of the system) for the quantum trajectories and show that these effects of the higher-order USC processes can be revealed in experiments by constructing histograms of detected quantum jumps. We illustrate these ideas by analyzing the excitation of two atoms by a single photon~[Garziano et al., Phys. Rev. Lett.117, 043601 (2016)]. For example, quantum trajectories reveal that keeping track of the quantum jumps from the atoms allow to reconstruct both the oscillations between one photon and two atoms, as well as emerging Rabi oscillations between the two atoms.

Open quantum dynamics for plant motions

Stochastic Schrödinger equations that govern the dynamics of open quantum systems are given by the equations for signal processing. In particular, the Brownian motion that drives the wave function of the system does not represent noise, but provides purely the arrival of new information. Thus the wave function is guided by the optimal signal detection about the conditions of the environments under noisy observations. This behaviour is similar to biological systems that detect environmental cues, process this information, and adapt to them optimally by minimising uncertainties about the conditions of their environments. It is postulated that information-processing capability is a fundamental law of nature, and hence that models describing open quantum systems can equally be applied to biological systems to model their dynamics. For illustration, simple stochastic models are considered to capture heliotropic and gravitropic motions of plants. The advantage of such dynamical models is that they allow for the quantification of information processed by the plants. By considering the consequence of information erasure, it is argued that biological systems can process environmental signals relatively close to the Landauer limit of computation, and that loss of information must lie at the heart of ageing in biological systems.

Quantum correlations of localized atomic excitations in a disordered atomic chain

The atom-waveguide interface mediates significant and long-range light-matter interactions through guided modes. In this one-dimensional system, we theoretically investigate the excitation localization of multiple atomic excitations under strong position disorder. Deep in the localization side, we obtain the time evolutions of quantum correlations via Kubo cumulant expansions, which arise initially and become finite and leveled afterward, overtaking those without disorder. This indicates two distinct regimes in time: before the onset of excitation localization, the disorder engage the disturbance of quantum correlations, which is followed by disorder-assisted buildup of quantum correlations that maintain at a later stage owing to the absence of excitations diffusion. The crossing of distinct regimes is pushed further in time for longer-range correlations, which indicates a characteristic timescale needed for disorder to sustain them. We also explore the effect of directionality of couplings and resonant dipole-dipole interactions, which can drive the system toward the delocalized side when it is under chiral couplings or large dipole-dipole interaction strengths. The time-evolved quantum correlations can give insights into the studies of few-body localization phenomenon and nonequilibrium dynamics in open quantum systems.

Time-periodic interaction between a spin-pair: A quantum master equation approach

Quantum master equations describe the dynamics of open quantum systems. Such systems might be subjected to external drives and internal interactions such as direct dipolar couplings. Both the external drives and the internal interactions could be oscillatory functions of time. In particular, a spin pair undergoing magic angle sample spinning experiences such oscillatory interactions in time. In the present work, we analyze a spin pair with a time-periodic interaction by using our recently-proposed fluctuation-regulated quantum master equation [PRA, 97, 063837 (2018)] (FRQME). FRQME helps calculate the first as well as the second-order effects of all time-dependent interactions Hamiltonians. Using FRQME, we find the usual sideband pattern and also have shown the linewidths of the sidebands as a function of the coupling strength, the period of oscillation, and an environmental parameter. We show that at high spinning speed, the sidebands would have a Lorentzian shape. We predict that the dipolar contributions to the sidebands’ linewidths are smaller for a shorter correlation time of the environmental fluctuations.

Many-Body Quantum State Diffusion for Non-Markovian Dynamics in Strongly Interacting Systems.

Capturing non-Markovian dynamics of open quantum systems is generally a challenging problem, especially for strongly interacting many-body systems. In this Letter, we combine recently developed non-Markovian quantum state diffusion techniques with tensor network methods to address this challenge. As a first example, we explore a Hubbard-Holstein model with dissipative phonon modes, where this new approach allows us to quantitatively assess how correlations spread in the presence of non-Markovian dissipation in a 1D many-body system. We find regimes where correlation growth can be enhanced by these effects, offering new routes for dissipatively enhancing transport and correlation spreading, relevant for both solid state and cold atom experiments.

Digital Quantum Simulation of Open Quantum Systems Using Quantum Imaginary–Time Evolution

Quantum simulation on emerging quantum hardware is a topic of intense interest. While many studies focus on computing ground state properties or simulating unitary dynamics of closed systems, open quantum systems are an interesting target of study owing to their ubiquity and rich physical behavior. However, their non-unitary dynamics are also not natural to simulate on digital quantum devices. Here, we report algorithms for the digital quantum simulation of the dynamics of open quantum systems governed by a Lindblad equation using adaptations of the quantum imaginary time evolution (QITE) algorithm. We demonstrate the algorithms on IBM Quantum's hardware with simulations of the spontaneous emission of a two level system and the dissipative transverse field Ising model. Our work advances efforts to simulate the dynamics of open quantum systems on quantum hardware.

Capacity of non-Markovianity to boost the efficiency of molecular switches

Quantum resource theory formulations of thermodynamics offer a versatile tool for the study of fundamental limitations to the efficiency of physical processes, independently of the microscopic details governing their dynamics. Despite the ubiquitous presence of non-Markovian dynamics in open quantum systems at the nanoscale, rigorous proofs of their beneficial effects on the efficiency of quantum dynamical processes at the bio-molecular level have not been reported yet. Here we combine the resource theory of athermality with concepts from the theory of divisibility classes for quantum channels, to prove that memory effects can increase the efficiency of photoisomerization to levels that are not achievable under a purely thermal Markovian (i.e. memoryless) evolution. This provides rigorous evidence that memory effects can provide a resource in biological quantum dynamics, and, more generally, quantum thermodynamics at the nanoscale.

Dynamics of Open Quantum Systems I, Oscillation and Decay

We develop a framework to analyze the dynamics of a finite-dimensional quantum systemSin contact with a reservoirR. The full, interactingSRdynamics is unitary. The reservoir has a stationary state but otherwise dissipative dynamics. We identify a main part of the full dynamics, which approximates it for small values of theSRcoupling constant, uniformly for all timest≥0. The main part consists of explicit oscillating and decaying parts. We show that the reduced system evolution is Markovian for all times. The technical novelty is a detailed analysis of the link between the dynamics and the spectral properties of the generator of theSRdynamics, based on Mourre theory. We allow forSRinteractions with little regularity, meaning that the decay of the reservoir correlation function only needs to be polynomial in time, improving on the previously required exponential decay.In this work we distill the structural and technical ingredients causing the characteristic features of oscillation and decay of theSRdynamics. In the companion paper [27] we apply the formalism to the concrete case of anN-level system linearly coupled to a spatially infinitely extended thermal bath of non-interacting Bosons.

Controlling the dynamics of dissipationless localized bound states in open quantum systems with periodic driving fields

In this paper, we study the exact dynamics of open quantum systems in the presence of dissipationless localized bound states and periodic driving fields. We show that different from the static adjustment of the system parameters to control the existence of localized bound states in open quantum systems, the periodic driving can either modulate the dynamics of the existing localized bound states or make some of them disappear due to the sideband generation from the driving that makes electrons transit between the localized bound states and continuous states. This analysis is also different from the widely used Floquet theory in the study of the driving effect on open system dynamics. We also find the conditions for the protection of the dissipationless localized bound states from the driving and for the manipulation of quantum coherence between localized bound states using the driving. The results for the dynamics of localized bound states have a potential application in controlling the quantum state against decoherence for the sake of its sensitivity to the fundamental frequency of the driving field and its strength.

Convergence Guarantees for Discrete Mode Approximations to Non-Markovian Quantum Baths.

Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. The non-Markovian environment is often approximated by discrete bosonic modes, thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been previously proposed, the resulting approximation lacks rigorous and general convergence guarantees. In this Letter, we show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that this approximation error typically falls off polynomially with the number of modes. Our results lend rigor to classical and quantum algorithms for approximating non-Markovian dynamics.

Diverging Quantum Speed Limits: A Herald of Classicality

When is the quantum speed limit (QSL) really quantum? While vanishing QSL times often indicate emergent classical behavior, it is still not entirely understood what precise aspects of classicality are at the origin of this dynamical feature. Here, we show that vanishing QSL times (or, equivalently, diverging quantum speeds) can be traced back to reduced uncertainty in quantum observables and thus can be understood as a consequence of emerging classicality for these particular observables. We illustrate this mechanism by developing a QSL formalism for continuous variable quantum systems undergoing general Gaussian dynamics. For these systems, we show that three typical scenarios leading to vanishing QSL times, namely large squeezing, small effective Planck's constant, and large particle number, can be fundamentally connected to each other. In contrast, by studying the dynamics of open quantum systems and mixed states, we show that the classicality that emerges due to incoherent mixing of states from the addition of classical noise typically increases the QSL time.

Adiabatic terminator for fermionic hierarchical equations of motion

The hierarchical equation of motion method has become one of the most popular numerical methods for describing the dissipative dynamics of open quantum systems linearly coupled to environment. However, its applications to systems with strong electron correlation are largely restrained by the computational cost, which is mainly caused by the high truncation tier L required to accurately characterize the strong correlation effect. In this work, we develop an adiabatic terminator by decoupling the principal dissipation mode with the fastest dissipation rate from the slower ones. The adiabatic terminator leads to substantially enhanced convergence with respect to L as demonstrated by the numerical tests carried out on a single impurity Anderson model. Moreover, the adiabatic terminator alleviates the numerical instability problems in the long-time dissipative dynamics.

Quantum optimal control for pure-state preparation using one initial state

This paper presents a framework for solving the pure-state preparation problem using numerical optimal control. As an example, we consider the case where a number of qubits are dispersively coupled to a readout cavity. We model open system quantum dynamics using the Markovian Lindblad master equation driven by external control pulses. The main result of this paper develops a basis of density matrices (a parameterization) where each basis element is a density matrix by itself. Utilizing a specific objective function, we show how an ensemble of the basis elements can be used as a single initial state throughout the optimization process—independent of the system dimension. We apply the general framework to the specific application of ground-state reset of one and two qubits coupled to a readout cavity.

Capturing non-Markovian dynamics with the reaction coordinate method

The reaction coordinate (RC) technique is emerging as a significant tool in the study of quantum dissipative dynamics and quantum thermodynamics. With the objective to further establish this tool, here we explore to what extent the method can capture non-Markovian dynamics of open quantum systems. As a case study, we focus on the pure decoherence model of a spin coupled to a harmonic reservoir. We compare the spin dynamics and measures for non-Markovianity from the exact analytical solution to simulations based on the RC method at the level of a second-order quantum master equation. We find that the RC method can quantitatively capture non-Markovian effects at strong system-bath coupling and for structured baths. This is rationalized by the fact that the collective RC bath mode, which is made part of the system, maintains system-bath correlations. Lastly, we apply our RC method and study the spin-boson model in the non-Markovian regime.