Nonperturbative Treatment of non-Markovian Dynamics of Open Quantum Systems

Quantum simulation of spin-boson models with structured bath

The spin-boson model, involving spins interacting with a bath of quantum harmonic oscillators, is a widely used representation of open quantum systems that describe many dissipative processes in physical, chemical and biological systems. Trapped ions present an ideal platform for simulating the quantum dynamics of such models, by accessing both the high-quality internal qubit states and the motional modes of the ions for spins and bosons, respectively. We demonstrate a fully programmable method to simulate dissipative dynamics of spin-boson models using a chain of trapped ions, where the initial temperature and the spectral densities of the boson bath are engineered by controlling the state of the motional modes and their coupling with qubit states. Our method provides a versatile and precise experimental tool for studying open quantum systems.

Variational Approach to Entangled Non-Hermitian Open Systems

The pseudomode model effectively captures the nonperturbative dynamics of open quantum systems with significantly reduced degrees of freedom. However, it is limited by the exponential growth of the Hilbert space dimension. To overcome the computational challenges, we propose a novel method that combines the multiple Davydov Ansatz with the Choi-Jamiolkowski isomorphism. Within this framework, the Lindblad equation is transformed into the non-Hermitian Schrödinger equation in a double Hilbert space, with its dynamics determined using the time-dependent variational principle. Three cases are calculated to demonstrate the effectiveness of the proposed method. We first discuss how the Davydov Ansatz works for the model with a single pseudomode. Extending the method to multiple pseudomodes, we show that the Ansatz effectively circumvents the exponential growth of the Hilbert space. Additionally, the method is also capable of addressing potential intersections that emerge in multibath scenarios. This approach offers potential applicability to various types of pseudomode models and other dissipative systems, providing a promising tool for the studies of open quantum dynamics.

Real-time propagation of adaptive sampling selected configuration interaction wave functions

We have developed a new time propagation method, time-dependent adaptive sampling configuration interaction (TD-ASCI), to describe the dynamics of a strongly correlated system. We employ the short iterative Lanczos method as the time-integrator, which provides a unitary, norm-conserving, and stable long-time propagation scheme. We used the TD-ASCI method to evaluate the time-domain correlation functions of molecular systems. The accuracy of the correlation function was assessed by Fourier transforming into the frequency domain to compute the dipole-allowed absorption spectra. The Fourier transform (FT) has been carried out with a short-time signal of the correlation function to reduce the computation time, using an efficient alternative FT scheme based on the ESPRIT signal processing algorithm. We have applied the TD-ASCI method to prototypical strongly correlated molecular systems and compared the absorption spectra to spectra evaluated using the equation of motion coupled cluster method with a truncation at the singles, doubles, and triples level.

Simulating many-body open quantum systems by harnessing the power of artificial intelligence and quantum computing

Simulating many-body open quantum systems (OQSs) is challenging due to the intricate interplay between the system and its environment, resulting in strong quantum correlations in both space and time. This Perspective presents an overview of recently developed theoretical methods using artificial intelligence (AI) and quantum computing (QC) to simulate the dynamics of these systems. We briefly introduce the dissipaton-embedded quantum master equation in second quantization, which provides a single master equation suitable for representation by neural quantum states or quantum circuits. The promising performance of AI- and QC-based approaches is demonstrated through preliminary research on simulating the quantum dissipative dynamics of many-body OQSs. We also discuss the limitations and future developments of these methods, which hold promise for overcoming the computational challenges associated with many-body OQS dynamics.

On-demand heralded MIR single-photon source using a cascaded quantum system

We propose a mechanism for generating single photons in the mid-infrared (MIR) using a solid-state or molecular quantum emitter. The scheme uses cavity quantum electrodynamics (QED) effects to selectively enhance a Frank-Condon transition, deterministically preparing a single Fock state of a polar phonon mode. By coupling the phonon mode to an antenna, the resulting excitation is then radiated to the far field as a single photon with a frequency matching the phonon mode. By combining macroscopic QED calculations with methods from open quantum system theory, we show that optimal parameters to generate these MIR photons occur for modest light-matter coupling strengths, which are achievable with state-of-the-art technologies. Combined, the cascaded system we propose provides a quasi-deterministic source of heralded single photons in a regime of the electromagnetic spectrum where this previously was not possible.

Efficient In Situ Generation of Photon-Memory Entanglement in a Nonlinear Cavity

Parametrically driving an optical cavity that simultaneously couples to an atomic ensemble quantum memory enables in situ generation of multimode photon-memory entanglement. A high-rate bipartite photon-memory entanglement can be generated even after discarding one entangled optical mode. This protocol can be realized with existing technologies based on photonic resonators integrated with a rare-earth-ion doped quantum memory. The proposed scheme shows significant advantages in entanglement generation rates compared with prevailing quantum memory protocols and experiments, with theoretical ebit rates of tens of MHz without fine-tuned operating conditions. Such a photon-memory entanglement source offers a versatile resource for quantum networking and interconnect applications.

Non-Markovian quantum exceptional points

Exceptional points (EPs) are singularities in the spectra of non-Hermitian operators where eigenvalues and eigenvectors coalesce. Open quantum systems have recently been explored as EP testbeds due to their non-Hermitian nature. However, most studies focus on the Markovian limit, leaving a gap in understanding EPs in the non-Markovian regime. This work addresses this gap by proposing a general framework based on two numerically exact descriptions of non-Markovian dynamics: the pseudomode equation of motion (PMEOM) and the hierarchical equations of motion (HEOM). The PMEOM is particularly useful due to its Lindblad-type structure, aligning with previous studies in the Markovian regime while offering deeper insights into EP identification. This framework incorporates non-Markovian effects through auxiliary degrees of freedom, enabling the discovery of additional or higher-order EPs that are inaccessible in the Markovian regime. We demonstrate the utility of this approach using the spin-boson model and linear bosonic systems.

Trapped-ion quantum simulation of electron transfer models with tunable dissipation

Electron transfer is at the heart of many fundamental physical, chemical, and biochemical processes essential for life. The exact simulation of these reactions is often hindered by the large number of degrees of freedom and by the essential role of quantum effects. Here, we experimentally simulate a paradigmatic model of molecular electron transfer using a multispecies trapped-ion crystal, where the donor-acceptor gap, the electronic and vibronic couplings, and the bath relaxation dynamics can all be controlled independently. By manipulating both the ground-state and optical qubits, we observe the real-time dynamics of the spin excitation, measuring the transfer rate in several regimes of adiabaticity and relaxation dynamics. Our results provide a testing ground for increasingly rich models of molecular excitation transfer processes that are relevant for molecular electronics and light-harvesting systems.

Simulating Non-Markovian Dynamics in Multidimensional Electronic Spectroscopy via Quantum Algorithm

Including the effect of the molecular environment in the numerical modeling of time-resolved electronic spectroscopy remains an important challenge in computational spectroscopy. In this contribution, we present a general approach for the simulation of the optical response of multichromophore systems in a structured environment and its implementation as a quantum algorithm. A key step of the procedure is the pseudomode embedding of the system-environment problem resulting in a finite set of quantum states evolving according to a Markovian quantum master equation. This formulation is then solved by a collision model integrated into a quantum algorithm designed to simulate linear and nonlinear response functions. The workflow is validated by simulating spectra for the prototypical excitonic dimer interacting with fast (memoryless) and finite-memory environments. The results demonstrate, on the one hand, the potential of the pseudomode embedding for simulating the dynamical features of nonlinear spectroscopy, including lineshape, spectral diffusion, and relaxations along delay times. On the other hand, the explicit synthesis of quantum circuits provides a fully quantum simulation protocol of nonlinear spectroscopy harnessing the efficient quantum simulation of many-body dynamics promised by the future generation of fault-tolerant quantum computers.

Spectral density modulation and universal Markovian closure of fermionic environments

The combination of chain-mapping and tensor-network techniques provides a powerful tool for the numerically exact simulation of open quantum systems interacting with structured environments. However, these methods suffer from a quadratic scaling with the physical simulation time, and therefore, they become challenging in the presence of multiple environments. This is particularly true when fermionic environments, well-known to be highly correlated, are considered. In this work, we first illustrate how a thermo-chemical modulation of the spectral density allows replacing the original fermionic environments with equivalent, but simpler, ones. Moreover, we show how this procedure reduces the number of chains needed to model multiple environments. We then provide a derivation of the fermionic Markovian closure construction, consisting of a small collection of damped fermionic modes undergoing a Lindblad-type dynamics and mimicking a continuum of bath modes. We describe, in particular, how the use of the Markovian closure allows for a polynomial reduction of the time complexity of chain-mapping based algorithms when long-time dynamics are needed.

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.

Quantized embedding approaches for collective strong coupling-Connecting ab initio and macroscopic QED to simple models in polaritonics

Collective light-matter interactions have been used to control chemistry and energy transfer, yet accessible approaches that combine ab initio methodology with large many-body quantum optical systems are missing due to the fast increase in computational cost for explicit simulations. We introduce an accessible ab initio quantum embedding concept for many-body quantum optical systems that allows us to treat the collective coupling of molecular many-body systems effectively in the spirit of macroscopic quantum electrodynamics while keeping the rigor of ab initio quantum chemistry for the molecular structure. Our approach fully includes the quantum fluctuations of the polaritonic field and yet remains much simpler and more intuitive than complex embedding approaches such as dynamical mean-field theory. We illustrate the underlying assumptions by comparison to the Tavis-Cummings model. The intuitive application of the quantized embedding approach and its transparent limitations offer a practical framework for the field of ab initio polaritonic chemistry to describe collective effects in realistic molecular ensembles.

Dynamics of a small quantum system open to a bath with thermostat

We investigate dynamics of a small quantum system open to a bath with thermostat. We introduce another bath, called a superbath, weakly coupled with the bath to provide it with a thermostat, which has either the Lindblad or Redfield type. We treat the interaction between the system and bath via a rigorous perturbation theory. Due to the thermostat, the bath behaves dissipative and stochastic, for which the usual Born-Markov assumption is not needed. We consider a specific example of a harmonic oscillator system of interest and a photonic bath in a large container, and a superbath of the Caldeira-Legget oscillators distributed on the inner surface of the container. After taking the trace over the superbath states, we use the P representation for the total harmonic system of the system and bath. We derive the reduced time-evolution equation for the system by explicitly finding the correlation between the system and bath beyond the product state that was not obtainable in the previous theory for the system and bath isolated from environment, and marginalizing bath degrees of freedom. Remarkably, the associated dynamic equation for the system density matrix is of the same form as the Redfield master equation with different coefficients depending on thermostat used. We find the steady state does not depend on the thermostat but the time-dependent state does, which agrees with common expectation. We expect to apply our theory to general systems. Unlike the usual quantum master equations, our reduced dynamics allows investigation for time-dependent protocols, and nonequilibrium quantum stochastic dynamics will be investigated in the future.

Unifying methods for optimal control in non-Markovian quantum systems via process tensors

The large dimensionality of environments is the limiting factor in applying optimal control to open quantum systems beyond the Markovian approximation. Various methods exist to simulate non-Markovian systems, which effectively reduce the environment to a number of active degrees of freedom. Here, we show that several of these methods can be expressed in terms of a process tensor in the form of a matrix-product-operator, which serves as a unifying framework to show how they can be used in optimal control and to compare their performance. The matrix-product-operator form provides a general scheme for computing gradients using back propagation and allows the efficiency of the different methods to be compared via the bond dimensions of their respective process tensors.

Quantum neural network approach to Markovian dissipative dynamics of many-body open quantum systems

Numerous variational methods have been proposed for solving quantum many-body systems, but they often face exponentially increasing computational complexity as the Hilbert space dimension grows. To address this, we introduce a novel approach using quantum neural networks to simulate the dissipative dynamics of many-body open quantum systems. This method combines neural-network quantum state representation with the time-dependent variational principle, both implemented via quantum algorithms. This results in accurate open quantum dynamics described by the Lindblad quantum master equation, exemplified by the spin-boson and transverse field Ising models. Our approach avoids the computational expense of classical algorithms and demonstrates the potential advantages of quantum computing for many-body simulations. To reduce measurement errors, we introduce a projection reset procedure, which could benefit other quantum simulations. In addition, our approach can be extended to simulate non-Markovian quantum dynamics.

Photo-induced dynamics with continuous and discrete quantum baths

The ultrafast quantum dynamics of photophysical processes in complex molecules is an extremely challenging computational problem with a broad variety of fascinating applications in quantum chemistry and biology. Inspired by recent developments in open quantum systems, we introduce a pure-state unraveled hybrid-bath method that describes a continuous environment via a set of discrete, effective bosonic degrees of freedom using a Markovian embedding. Our method is capable of describing both, a continuous spectral density and sharp peaks embedded into it. Thereby, we overcome the limitations of previous methods, which either capture long-time memory effects using the unitary dynamics of a set of discrete vibrational modes or use memoryless Markovian environments employing a Lindblad or Redfield master equation. We benchmark our method against two paradigmatic problems from quantum chemistry and biology. We demonstrate that compared to unitary descriptions, a significantly smaller number of bosonic modes suffices to describe the excitonic dynamics accurately, yielding a computational speed-up of nearly an order of magnitude. Furthermore, we take into account explicitly the effect of a δ-peak in the spectral density of a light-harvesting complex, demonstrating the strong impact of the long-time memory of the environment on the dynamics.

Entropy production in the mesoscopic-leads formulation of quantum thermodynamics

Understanding the entropy production of systems strongly coupled to thermal baths is a core problem of both quantum thermodynamics and mesoscopic physics. While many techniques exist to accurately study entropy production in such systems, they typically require a microscopic description of the baths, which can become numerically intractable to study for large systems. Alternatively an open-systems approach can be employed with all the nuances associated with various levels of approximation. Recently, the mesoscopic leads approach has emerged as a powerful method for studying such quantum systems strongly coupled to multiple thermal baths. In this method, a set of discretized lead modes, each locally damped, provide a Markovian embedding. Here we show that this method proves extremely useful to describe entropy production of a strongly coupled open quantum system. We show numerically, for both noninteracting and interacting setups, that a system coupled to a single bath exhibits a thermal fixed point at the level of the embedding. This allows us to use various results from the thermodynamics of quantum dynamical semigroups to infer the nonequilibrium thermodynamics of the strongly coupled, non-Markovian central systems. In particular, we show that the entropy production in the transient regime recovers the well-established microscopic definitions of entropy production with a correction that can be computed explicitly for both the single- and multiple-lead cases.

A mixed perturbative-nonperturbative treatment for strong light-matter interactions

The full information about the interaction between a quantum emitter and an arbitrary electromagnetic environment is encoded in the so-called spectral density. We present an approach for describing such interaction in any coupling regime, providing a Lindblad-like master equation for the emitter dynamics when coupled to a general nanophotonic structure. Our framework is based on the splitting of the spectral density into two terms. On the one hand, a spectral density responsible for the non-Markovian and strong-coupling-based dynamics of the quantum emitter. On the other hand, a residual spectral density including the remaining weak-coupling terms. The former is treated nonperturbatively with a collection of lossy interacting discrete modes whose parameters are determined by a fit to the original spectral density in a frequency region encompassing the quantum emitter transition frequencies. The latter is treated perturbatively under a Markovian approximation. We illustrate the power and validity of our approach through numerical simulations in three different setups, thus offering a variety of scenarios for a full test, including the ultra-strong coupling regime.

Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

Dynamics of a Quantum System Interacting with White Non-Gaussian Baths: Poisson Noise Master Equation

Quantum systems are unavoidably open to their surrounding degrees of freedom. The theory of open quantum systems is thus crucial to understanding the fluctuations, dissipation, and decoherence of a quantum system of interest. Typically, the bath is modeled as an ensemble of harmonic oscillators, which yields Gaussian statistics of the bath influence on the quantum systems. However, there are also phenomena in which the bath consists of two-state systems, spins, or anharmonic oscillators; therefore, the non-Gaussian properties of the bath become important. Nevertheless, a theoretical framework to describe quantum systems under the influence of such non-Gaussian baths is not well established. Here, we develop a theory to describe quantum dissipative systems affected by Poisson noise properties of the bath, because the Lévi-Itô decomposition theorem asserts that Poisson noise is fundamental in describing arbitrary white noise beyond Gaussian properties. We introduce a quantum bath model that allows for the consistent description of dissipative quantum systems. The obtained master equation reveals non-Gaussian bath effects in the white noise regime, and provides an essential step toward describing open quantum dynamics under the influence of generic baths.

Lindblad Master Equation Capable of Describing Hybrid Quantum Systems in the Ultrastrong Coupling Regime

Despite significant theoretical efforts devoted to studying the interaction between quantized light modes and matter, the so-called ultrastrong coupling regime still presents significant challenges for theoretical treatments and prevents the use of many common approximations. Here we demonstrate an approach that can describe the dynamics of hybrid quantum systems in any regime of interaction for an arbitrary electromagnetic (EM) environment. We extend a previous method developed for few-mode quantization of arbitrary systems to the case of ultrastrong light-matter coupling, and show that even such systems can be treated using a Lindblad master equation where decay operators act only on the photonic modes by ensuring that the effective spectral density of the EM environment is sufficiently suppressed at negative frequencies. We demonstrate the validity of our framework and show that it outperforms current state-of-the-art master equations for a simple model system, and then study a realistic nanoplasmonic setup where existing approaches cannot be applied.

Systematic Coarse Graining of Environments for the Nonperturbative Simulation of Open Quantum Systems

Conducting precise electronic-vibrational dynamics simulations of molecular systems poses significant challenges when dealing with realistic environments composed of numerous vibrational modes. Here, we introduce a technique for the construction of effective phonon spectral densities that capture accurately open-system dynamics over a finite time interval of interest. When combined with existing nonperturbative simulation tools, our approach can reduce significantly the computational costs associated with many-body open-system dynamics.

Certifying Multimode Light-Matter Interaction in Lossy Resonators

Quantum models based on few-mode master equations have been a central tool in the study of resonator quantum electrodynamics, extending the seminal single-mode Jaynes-Cummings model to include loss and multiple modes. Despite their broad application range, previous approaches within this framework have either relied on a Markov approximation or a fitting procedure. By combining ideas from pseudomode and quasinormal mode theory, we develop a certification criterion for multi-mode effects in lossy resonators. It is based on a witness observable, and neither requires a fitting procedure nor a Markov approximation. Using the resulting criterion, we demonstrate that such multi-mode effects are important for understanding previous experiments in x-ray cavity QED with Mössbauer nuclei and that they allow one to tune the nuclear ensemble properties.

Spin-Dependent Momentum Conservation of Electron-Phonon Scattering in Chirality-Induced Spin Selectivity

The elucidation of the mechanisms underpinning chirality-induced spin selectivity remains an outstanding scientific challenge. Here we consider the role of delocalized phonon modes in electron transport in chiral structures and demonstrate that spin selectivity can originate from spin-dependent energy and momentum conservation in electron-phonon scattering events. While this mechanism is robust to the specific nature of the vibrational modes, the degree of spin polarization depends on environmental factors, such as the specific temperature and phonon relaxation rates, as well as the presence of external driving fields. This parametric dependence is used to present experimentally testable predictions of our model.

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

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.

Anomalous spontaneous emission dynamics at chiral exceptional points

An open quantum system operated at the spectral singularities where dimensionality reduces, known as exceptional points (EPs), demonstrates distinguishing behavior from the Hermitian counterpart. Here, we present an analytical description of local density of states (LDOS) for microcavity featuring chiral EPs, and unveil the anomalous spontaneous emission dynamics from a quantum emitter (QE) due to the non-Lorentzian response of EPs. Specifically, we reveal that a squared Lorentzian term of LDOS contributed by chiral EPs can destructively interfere with the linear Lorentzian profile, resulting in the null Purcell enhancement to a QE with special transition frequency, which we call EP induced transparency. While for the case of constructive interference, the squared Lorentzian term can narrow the linewidth of Rabi splitting even below that of bare components, and thus significantly suppresses the decay of Rabi oscillation. Interestingly, we further find that an open microcavity with chiral EPs supports atom-photon bound states for population trapping and decay suppression in long-time dynamics. As applications, we demonstrate the advantages of microcavity operated at chiral EPs in achieving high-fidelity entanglement generation and high-efficiency single-photon generation. Our work unveils the exotic cavity quantum electrodynamics unique to chiral EPs, which opens the door for controlling light-matter interaction at the quantum level through non-Hermiticity, and holds great potential in building high-performance quantum-optics devices.

Fingerprint and Universal Markovian Closure of Structured Bosonic Environments

We exploit the properties of chain mapping transformations of bosonic environments to identify a finite collection of modes able to capture the characteristic features, or fingerprint, of the environment. Moreover we show that the countable infinity of residual bath modes can be replaced by a universal Markovian closure, namely, a small collection of damped modes undergoing a Lindblad-type dynamics whose parametrization is independent of the spectral density under consideration. We show that the Markovian closure provides a quadratic speedup with respect to standard chain mapping techniques and makes the memory requirement independent of the simulation time, while preserving all the information on the fingerprint modes. We illustrate the application of the Markovian closure to the computation of linear spectra but also to nonlinear spectral response, a relevant experimentally accessible many body coherence witness for which efficient numerically exact calculations in realistic environments are currently lacking.

Few-mode field quantization for multiple emitters

The control of the interaction between quantum emitters using nanophotonic structures holds great promise for quantum technology applications, while its theoretical description for complex nanostructures is a highly demanding task as the electromagnetic (EM) modes form a high-dimensional continuum. We here introduce an approach that permits a quantized description of the full EM field through a small number of discrete modes. This extends the previous work in ref. (I. Medina, F. J. García-Vidal, A. I. Fernández-Domínguez, and J. Feist, "Few-mode field quantization of arbitrary electromagnetic spectral densities," Phys. Rev. Lett., vol. 126, p. 093601, 2021) to the case of an arbitrary number of emitters, without any restrictions on the emitter level structure or dipole operators. The low computational demand of this method makes it suitable for studying dynamics for a wide range of parameters. We illustrate the power of our approach for a system of three emitters placed within a hybrid metallodielectric photonic structure and show that excitation transfer is highly sensitive to the properties of the hybrid photonic-plasmonic modes.

Unveiling atom-photon quasi-bound states in hybrid plasmonic-photonic cavity

Dissipation, often associated with plasmons, leads to decoherence and is generally considered fatal for quantum nonlinearities and entanglement. Counterintuitively, by introducing a dissipative plasmonic nanoantenna into a typical cavity quantum electrodynamics (QED) system, we unveil the wide existence of the atom-photon quasi-bound state (qBS), a kind of exotic eigenstate with anomalously small decay, in the hybrid plasmonic-photonic cavity. To derive the analytical condition of atom-photon qBS, we formulate a quantized two-mode model of the local density of states by connecting the interacting uncoupled cavity modes to the macroscopic QED. With resonant plasmon-photon coupling, we showcase the single-atom qBS that improves the efficiency of single-photon generation over one order of magnitude; and the two-atom qBS that significantly enhances spontaneous entanglement generation compared with a bare photonic cavity. Notably, such single-atom and multi-atom qBS can be simultaneously accessed in realistic plasmonic-photonic cavities, providing a versatile platform for advanced quantum technologies, such as quantum light sources, quantum computation, and quantum information.

A statistical quasi-particles thermofield theory with Gaussian environments: System-bath entanglement theorem for nonequilibrium correlation functions

For open quantum systems, environmental dissipative effect can be represented by statistical quasi-particles, namely dissipatons. We exploit this fact to establish the dissipaton thermofield theory. The resulting generalized Langevin dynamics of absorptive and emissive thermofield operators are effectively noise-resolved. The system-bath entanglement theorem is then readily followed between a important class of nonequilibrium steady-state correlation functions. All these relations are validated numerically. A simple corollary is the transport current expression, which exactly recovers the result obtained from the nonequilibrium Green's function formalism.

Exact simulation of pigment-protein complexes unveils vibronic renormalization of electronic parameters in ultrafast spectroscopy

The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra.

Multimode vibronic mixing in model photosynthetic systems revealed by numerically exact simulations is shown to strongly modify linear and non-linear optical responses and facilitate the persistence of coherent dynamics.

Universal time-domain Prony fitting decomposition for optimized hierarchical quantum master equations

In this work, we propose the time-domain Prony fitting decomposition (t-PFD) as an accurate and effcient exponential series method, applicable to arbitrary bath correlation functions. The resulting numerical effciency of hierarchical equations of motion (HEOM) formalism is greatly optimized, especially in low temperature regimes that would be inaccessible with other methods. For demonstration, we calibrate the present t-PFD against the celebrated Padé spectrum decomposition (PSD) method, followed by converged HEOM evaluations on the single-impurity Anderson model system.

Correlations, Information Backflow, and Objectivity in a Class of Pure Dephasing Models

We critically examine the role that correlations established between a system and fragments of its environment play in characterising the ensuing dynamics. We employ a dephasing model with different initial conditions, where the state of the initial environment represents a tunable degree of freedom that qualitatively and quantitatively affects the correlation profiles, but nevertheless results in the same reduced dynamics for the system. We apply recently developed tools for the characterisation of non-Markovianity to carefully assess the role that correlations, as quantified by the (quantum) Jensen-Shannon divergence and relative entropy, as well as changes in the environmental state, play in whether the conditions for classical objectivity within the quantum Darwinism paradigm are met. We demonstrate that for precisely the same non-Markovian reduced dynamics of the system arising from different microscopic models, some exhibit quantum Darwinistic features, while others show that no meaningful notion of classical objectivity is present. Furthermore, our results highlight that the non-Markovian nature of an environment does not a priori prevent a system from redundantly proliferating relevant information, but rather it is the system's ability to establish the requisite correlations that is the crucial factor in the manifestation of classical objectivity.

Quantum-parametric-oscillator heat engines in squeezed thermal baths: Foundational theoretical issues

In this paper we examine some foundational issues of a class of quantum engines where the system consists of a single quantum parametric oscillator, operating in an Otto cycle consisting of four stages of two alternating phases: the isentropic phase is detached from any bath (thus a closed system) where the natural frequency of the oscillator is changed from one value to another, and the isothermal phase where the system (now rendered open) is put in contact with one or two squeezed baths of different temperatures, whose nonequilibrium dynamics follows the Hu-Paz-Zhang (HPZ) master equation for quantum Brownian motion. The HPZ equation is an exact non-Markovian equation which preserves the positivity of the density operator and is valid for (1) all temperatures, (2) arbitrary spectral density of the bath, and (3) arbitrary coupling strength between the system and the bath. Taking advantage of these properties we examine some key foundational issues of theories of quantum open and squeezed systems for these two phases of the quantum Otto engines. This includes (1) the non-Markovian regimes for non-Ohmic, low-temperature baths, (2) what to expect in nonadiabatic frequency modulations, (3) strong system-bath coupling, as well as (4) the proper junction conditions between these two phases. Our aim here is not to present ways for attaining higher efficiency but to build a more solid theoretical foundation for quantum engines of continuous variables covering a broader range of parameter spaces that we hope are of use for exploring such possibilities.

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.

Efficient Exploration of Hamiltonian Parameter Space for Optimal Control of Non-Markovian Open Quantum Systems

We present a general method to efficiently design optimal control sequences for non-Markovian open quantum systems, and illustrate it by optimizing the shape of a laser pulse to prepare a quantum dot in a specific state. The optimization of control procedures for quantum systems with strong coupling to structured environments-where time-local descriptions fail-is a computationally challenging task. We modify the numerically exact time evolving matrix product operator (TEMPO) method, such that it allows the repeated computation of the time evolution of the reduced system density matrix for various sets of control parameters at very low computational cost. This method is potentially useful for studying numerous optimal control problems, in particular in solid state quantum devices where the coupling to vibrational modes is typically strong.

Few-Mode Field Quantization of Arbitrary Electromagnetic Spectral Densities

We develop a framework that provides a few-mode master equation description of the interaction between a single quantum emitter and an arbitrary electromagnetic environment. The field quantization requires only the fitting of the spectral density, obtained through classical electromagnetic simulations, to a model system involving a small number of lossy and interacting modes. We illustrate the power and validity of our approach by describing the population and electric field spatial dynamics in the spontaneous decay of an emitter placed in a complex hybrid plasmonic-photonic structure.

Simulating Quantum Vibronic Dynamics at Finite Temperatures With Many Body Wave Functions at 0 K

For complex molecules, nuclear degrees of freedom can act as an environment for the electronic “system” variables, allowing the theory and concepts of open quantum systems to be applied. However, when molecular system-environment interactions are non-perturbative and non-Markovian, numerical simulations of the complete system-environment wave function become necessary. These many body dynamics can be very expensive to simulate, and extracting finite-temperature results—which require running and averaging over many such simulations—becomes especially challenging. Here, we present numerical simulations that exploit a recent theoretical result that allows dissipative environmental effects at finite temperature to be extracted efficiently from a single, zero-temperature wave function simulation. Using numerically exact time-dependent variational matrix product states, we verify that this approach can be applied to vibronic tunneling systems and provide insight into the practical problems lurking behind the elegance of the theory, such as the rapidly growing numerical demands that can appear for high temperatures over the length of computations.

Matrix Product State Simulations of Non-Equilibrium Steady States and Transient Heat Flows in the Two-Bath Spin-Boson Model at Finite Temperatures

Simulating the non-perturbative and non-Markovian dynamics of open quantum systems is a very challenging many body problem, due to the need to evolve both the system and its environments on an equal footing. Tensor network and matrix product states (MPS) have emerged as powerful tools for open system models, but the numerical resources required to treat finite-temperature environments grow extremely rapidly and limit their applications. In this study we use time-dependent variational evolution of MPS to explore the striking theory of Tamascelli et al. (Phys. Rev. Lett. 2019, 123, 090402.) that shows how finite-temperature open dynamics can be obtained from zero temperature, i.e., pure wave function, simulations. Using this approach, we produce a benchmark dataset for the dynamics of the Ohmic spin-boson model across a wide range of coupling strengths and temperatures, and also present a detailed analysis of the numerical costs of simulating non-equilibrium steady states, such as those emerging from the non-perturbative coupling of a qubit to baths at different temperatures. Despite ever-growing resource requirements, we find that converged non-perturbative results can be obtained, and we discuss a number of recent ideas and numerical techniques that should allow wide application of MPS to complex open quantum systems.

Analytic expressions for the steady-state current with finite extended reservoirs

Open-system simulations of quantum transport provide a platform for the study of true steady states, Floquet states, and the role of temperature, time dynamics, and fluctuations, among other physical processes. They are rapidly gaining traction, especially techniques that revolve around "extended reservoirs," a collection of a finite number of degrees of freedom with relaxation that maintains a bias or temperature gradient, and have appeared under various guises (e.g., the extended or mesoscopic reservoir, auxiliary master equation, and driven Liouville-von Neumann approaches). Yet, there are still a number of open questions regarding the behavior and convergence of these techniques. Here, we derive general analytical solutions, and associated asymptotic analyses, for the steady-state current driven by finite reservoirs with proportional coupling to the system/junction. In doing so, we present a simplified and unified derivation of the non-interacting and many-body steady-state currents through arbitrary junctions, including outside of proportional coupling. We conjecture that the analytic solution for proportional coupling is the most general of its form for isomodal relaxation (i.e., relaxing proportional coupling will remove the ability to find compact, general analytical expressions for finite reservoirs). These results should be of broad utility in diagnosing the behavior and implementation of extended reservoir and related approaches, including the convergence to the Landauer limit (for non-interacting systems) and the Meir-Wingreen formula (for many-body systems).

Universal non-Markovianity detection in hybrid open quantum systems

A universal characterization of non-Markovianity for any open hybrid quantum systems is presented. This formulation is based on the negativity volume of the generalized Wigner function, which serves as an indicator of the quantum correlations in any composite quantum systems. It is shown, that the proposed measure can be utilized for any single or multi-partite quantum system, containing any discrete or continuous variables. To demonstrate its power in revealing non-Markovianity in such quantum systems, we additionally consider a few illustrative examples.

Universal Anti-Kibble-Zurek Scaling in Fully Connected Systems

We investigate the quench dynamics of an open quantum system involving a quantum phase transition. In the isolated case, the quench dynamics involving the phase transition exhibits a number of scaling relations with the quench rate as predicted by the celebrated Kibble-Zurek mechanism. In contact with an environment however, these scaling laws break down and one may observe an anti-Kibble-Zurek behavior: slower ramps lead to less adiabatic dynamics, increasing thus nonadiabatic effects with the quench time. In contrast to previous works, we show here that such anti-Kibble-Zurek scaling can acquire a universal form in the sense that it is determined by the equilibrium critical exponents of the phase transition, provided the excited states of the system exhibit singular behavior, as observed in fully connected models. This demonstrates novel universal scaling laws granted by a system-environment interaction in a critical system. We illustrate these findings in two fully connected models, namely, the quantum Rabi and the Lipkin-Meshkov-Glick models. In addition, we discuss the impact of nonlinear ramps and finite-size systems.

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.

Pushing the limits of the reaction-coordinate mapping

The reaction-coordinate mapping is a useful technique to study complex quantum dissipative dynamics into structured environments. In essence, it aims to mimic the original problem by means of an "augmented system," which includes a suitably chosen collective environmental coordinate-the "reaction coordinate." This composite then couples to a simpler "residual reservoir" with short-lived correlations. If, in addition, the residual coupling is weak, a simple quantum master equation can be rigorously applied to the augmented system, and the solution of the original problem just follows from tracing out the reaction coordinate. But, what if the residual dissipation is strong? Here, we consider an exactly solvable model for heat transport-a two-node linear "quantum wire" connecting two baths at different temperatures. We allow for a structured spectral density at the interface with one of the reservoirs and perform the reaction-coordinate mapping, writing a perturbative master equation for the augmented system. We find that (a) strikingly, the stationary state of the original problem can be reproduced accurately by a weak-coupling treatment even when the residual dissipation on the augmented system is very strong, (b) the agreement holds throughout the entire dynamics under large residual dissipation in the overdamped regime; and (c) such a master equation can grossly overestimate the stationary heat current across the wire, even when its nonequilibrium steady state is captured faithfully. These observations can be crucial when using the reaction-coordinate mapping to study the largely unexplored strong-coupling regime in quantum thermodynamics.

Dissipation-Assisted Matrix Product Factorization

Charge and energy transfer in biological and synthetic organic materials are strongly influenced by the coupling of electronic states to a highly structured dissipative environment. Nonperturbative simulations of these systems require a substantial computational effort, and current methods can only be applied to large systems if environmental structures are severely coarse grained. Time evolution methods based on tensor networks are fundamentally limited by the times that can be reached due to the buildup of entanglement in time, which quickly increases the size of the tensor representation, i.e., the bond dimension. In this Letter, we introduce a dissipation-assisted matrix product factorization (DAMPF) method that combines a tensor network representation of the vibronic state within a pseudomode description of the environment where a continuous bosonic environment is mapped into a few harmonic oscillators under Lindblad damping. This framework is particularly suitable for a tensor network representation, since damping suppresses the entanglement growth among oscillators and significantly reduces the bond dimension required to achieve a desired accuracy. We show that dissipation removes the "time-wall" limitation of existing methods, enabling the long-time simulation of large vibronic systems consisting of 10-50 sites coupled to 100-1000 underdamped modes in total and for a wide range of parameter regimes. For these reasons, we believe that our formalism will facilitate the investigation of spatially extended systems with applications to quantum biology, organic photovoltaics, and quantum thermodynamics.

Efficient Simulation of Finite-Temperature Open Quantum Systems

Chain-mapping techniques in combination with the time-dependent density matrix renormalization group are a powerful tool for the simulation of open-system quantum dynamics. For finite-temperature environments, however, this approach suffers from an unfavorable algorithmic scaling with increasing temperature. We prove that the system dynamics under thermal environments can be nonperturbatively described by temperature-dependent system-environmental couplings with the initial environment state being in its pure vacuum state, instead of a mixed thermal state. As a consequence, as long as the initial system state is pure, the global system-environment state remains pure at all times. The resulting speed-up and relaxed memory requirements of this approach enable the efficient simulation of open quantum systems interacting with highly structured environments in any temperature range, with applications extending from quantum thermodynamics to quantum effects in mesoscopic systems.