The calculation of the non-Markovian dynamics of waveguide-QED systems with coherent time-delayed feedback within the Heisenberg picture gives access to the microscopic dynamics and allows the inclusion of dissipation channels. However, typically a hierarchy of multi-time correlations arises. We propose to unravel these multi-time correlations via a projection onto a complete set of states in the Hilbert space. In particular, we consider a two-level system inside a semi-infinite waveguide subjected to quantum light pulses. Our approach makes it possible to calculate the dynamics in a numerically exact and efficient manner for multiphoton pulses of arbitrary shape with memory requirements known a priori.
The calculation of the non-Markovian dynamics of waveguide-QED systems with coherent time-delayed feedback within the Heisenberg picture gives access to the microscopic dynamics and allows the inclusion of dissipation channels. However, typically a hierarchy of multi-time correlations arises. We propose to unravel these multi-time correlations via a projection onto a complete set of states in the Hilbert space. In particular, we consider a two-level system inside a semi-infinite waveguide subjected to quantum light pulses. Our approach makes it possible to calculate the dynamics in a numerically exact and efficient manner for multiphoton pulses of arbitrary shape with memory requirements known a priori.
The dynamics of waveguide-QED systems involving coherent time-delayed feedback give rise to a hierarchy of multitime correlations within the Heisenberg picture due to the induced non-Markovianity. We propose to perform a projection onto a complete set of states in the Hilbert space to decompose the multitime correlations into single-time matrix elements. To illustrate the procedure, we consider the paradigmatic example of a two-level system that couples to a semi-infinite waveguide and interacts with quantum light pulses. Our approach complements the range of available methods as it allows calculating the dynamics under the inclusion of additional dissipation channels in a numerically exact and efficient manner for multiphoton pulses of arbitrary shape where memory requirements are known in advance.
We present a Heisenberg-operator-based formulation of coherent quantum feedback and Pyragas control. This model is easy to implement and allows for an efficient and fast calculation of the dynamics of feedback-driven observables as the number of contributing correlations grows in systems with a fixed number of excitations only linearly in time. Furthermore, our model unravels the quantum kinetics of entanglement growth in the system by explicitly calculating non-Markovian multi-time correlations, e.g., how the emission of a photon is correlated with an absorption process in the past. Therefore, the time-delayed differential equations are expressed in terms of insightful physical quantities. Another considerable advantage of this method is its compatibility with typical approximation schemes, such as factorization techniques and the semi-classical treatment of coherent fields. This allows the application on a variety of setups, ranging from closed quantum systems in the few excitation regimes to open systems and Pyragas control in general.
We study the nonlinear spectra and multiphoton correlation functions of a waveguide system coupled to a two-level system (including realistic dissipation channels) with a time-delayed coherent feedback. We compute these observables by extending a recent quantum trajectory discretized-waveguide (QTDW) approach which exploits quantum trajectory simulations and a collisional model for the waveguide to efficiently simulate the quantum dynamics. Following a description of the general technique, we show how to calculate the first- and second-order quantum correlation functions, in the presence of a coherent pumping field. With a short delay time, we show how feedback can be used to filter out the central peak of the Mollow triplet or switch the output between bunched and antibunched photons by proper choice of round trip phase. We further show how the loop length and round trip phase affect the zero-time second-order quantum correlation function, an indicator of bunching or antibunching. New resonances introduced through the feedback loop are also shown through their appearance in the incoherent output spectrum from the waveguide. We explain these results in the context of the waiting time distributions of the system output and individual trajectories, uniquely stochastic observables that are easily accessible with the QTDW model.
Within the matrix product state framework, we study the non-Markovian feedback dynamics of a two-level system interacting with the electromagnetic field inside a semi-infinite waveguide where the formation of an atom-photon bound state is possible. Taking the steady-state excitation of the emitter as a figure of merit, we compare the trapped excitation for an initially excited quantum emitter and an emitter prepared via quantized pulses. In the latter case, we find that for certain feedback delay times, multi-photon pulses yield a significantly higher steady-state excitation than possible with an initially excited emitter. In the system we consider, we find an increase of around 90 percent. Thus, the quantum optical initialization is found to be more efficient due to the self-consistent inclusion of emerging quantum noise. Furthermore, we propose a recursive Heisenberg ansatz which shows good qualitative agreement with the numerically expensive matrix product state method and requires only numerical resources growing quadratically with the number of excitations.
We present two different methods for modeling non-Markovian quantum light-matter interactions in waveguide QED systems, using matrix product states (MPSs) and a space-discretized waveguide (SDW) model.After describing the general theory and implementation of both approaches, we compare and contrast these methods directly on three topical problems of interest in waveguide-QED, including (i) a two-level system (TLS) coupled to an infinite (one-dimensional) waveguide, (ii) a TLS coupled to a terminated waveguide with a time-delayed coherent feedback, and (iii) two spatially separated TLSs coupled within an infinite waveguide.Both approaches are shown to efficiently describe multiphoton nonlinear dynamics in highly non-Markovian regimes, and we highlight the advantages and disadvantages of these methods for modeling waveguide QED interactions, including their implementation in PYTHON, computational run times, and ease of conceptual understanding.We explore both vacuum dynamics as well as regimes of strong optical pumping, where a weak excitation approximation cannot be applied.The MPS approach scales better when modeling multiphoton dynamics and long delay times and explicitly includes non-Markovian memory effects.In contrast, the SDW model accounts for non-Markovian effects through space discretization and solves Markovian equations of motion, yet rigorously includes the effects of retardation.The SDW model, based on an extension of recent collisional pictures in quantum optics, is solved through quantum trajectory techniques and can more easily add in additional dissipation processes, including off-chip decay and TLS pure dephasing.The impact of these processes is shown directly on feedback-induced population trapping and TLS entanglement between spatially separated TLSs.
Phase-space methods are powerful tools for the treatment of Markovian open quantum systems: they map the reduced dynamics of a system S, in interaction with an environment E, exactly onto Langevin equations for c-number stochastic variables, as opposed to Heisenberg-Langevin equations for operators. Langevin equations provide analytical insight in key regimes and excel at handling strong nonlinearities and couplings, where other methods often falter. Extending phase-space methods to non-Markovian dynamics, however, has remained a long-standing challenge. Here we address this gap by applying phase-space representations to the full S+E system; integrating out the environmental degrees of freedom then yields a general Langevin framework for S that incorporates both deterministic and stochastic contributions from E. Normally ordered representations, such as the Glauber-Sudarshan P representation and its positive variant due to Drummond and Gardiner, lead to Langevin equations in which (i) non-Markovian effects emerge exclusively in the deterministic terms, via a memory kernel, and (ii) noise contributions vanish when E is initially in the vacuum state. To demonstrate the power of this framework, we address the paradigmatic problem of delayed coherent feedback, in which the system is driven by its own past state, and study its impact on the laser linewidth: we recover the narrowing observed well above threshold and predict an enhanced narrowing just above it. Crucially, the number of stochastic variables scales linearly with the system size, making the framework suitable for problems ranging from a few degrees of freedom to genuinely many-body systems. This opens the way to the systematic study of non-Markovian driven-dissipative quantum systems using the same analytical and numerical tools that have long made phase-space methods so successful in the Markovian regime.