We present a matrix product state (MPS) method for including decoherence processes in calculations involving waveguide quantum electrodynamics (waveguide QED) using density matrices. The approach is based on collision quantum optics, where the many-body state of the waveguide is represented using discrete time bins, which are then efficiently encoded using an MPS chain. Our method is a generalization of previous MPS methods, and we demonstrate how one can efficiently extend the method to density matrices, allowing for the inclusion of various loss processes in the form of Lindblad terms in the Liouvillian superoperator responsible for the relevant dissipation dynamics. As an application of the theory, we study various waveguide QED systems and the influence of emitter pure dephasing (which is one of the most important processes in real systems) on the light-matter interactions, including a two-level system (TLS) in a semi-infinite waveguide with time-delayed feedback, two spatially separated TLSs with finite delays, and finally the scattering of few-photon Fock pulses on a TLS. In addition to emitter pure dephasing, we also show how to include off-chip radiative decay, and show how it differs qualitatively from pure dephasing.
We present a matrix product state (MPS) method for including decoherence processes in calculations involving waveguide quantum electrodynamics (waveguide QED) using density matrices. The approach is based on collision quantum optics, where the many-body state of the waveguide is represented using discrete time bins, which are then efficiently encoded using an MPS chain. Our method is a generalization of previous MPS methods, and we demonstrate how one can efficiently extend the method to density matrices, allowing for the inclusion of various loss processes in the form of Lindblad terms in the Liouvillian superoperator responsible for the relevant dissipation dynamics. As an application of the theory, we study various waveguide QED systems and the influence of emitter pure dephasing (which is one of the most important processes in real systems) on the light-matter interactions, including a two-level system (TLS) in a semi-infinite waveguide with time-delayed feedback, two spatially separated TLSs with finite delays, and finally the scattering of few-photon Fock pulses on a TLS. In addition to emitter pure dephasing, we also show how to include off-chip radiative decay, and show how it differs qualitatively from pure dephasing.
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.
Squeezed light states have a special place in quantum optics with potentially profound applications in emerging quantum technologies. We present a numerically-exact matrix product states (MPS) approach to model quantum pulses of squeezed light, at the few-photon level, interacting with a two-level system in a waveguide environment. We represent the squeezed state as a coherent superposition of Fock states and explore the nonlinear population dynamics as well as multi-photon correlation functions that emerge. We show how the squeezed pulse can create quantum correlations that are unique to squeezed pulses, including ⟨b(t)b(t+t′)⟩ for transmitted fields as well as ⟨b†(t)b†(t+t′)b(t+t′)b(t)⟩. We also demonstrate how ⟨b(t)b(t+t′)⟩, a first-order correlation function, shows nonlinear photon correlations that are similar to those known and measured for two-photon scattering states. Finally, we also study the squeezed spectra of the pulse before and after interacting with the two-level system, and highlight the role of the spectral bandwidth of the incident pulse. The MPS theory allows the modeling of arbitrary bandwidth squeezing without making any Markov and Born approximations for the light-matter interaction processes, and can easily be extended to waveguide systems with multiple emitters and time-delayed feedback.
We present a quantum dynamical study of pulsed few-photon scattering from a single artificial atom, consisting of a two-level system (TLS) or qubit, in a waveguide QED system, directly comparing and contrasting two different quantum theoretical simulation methods: (i) an input-output scattering approach that uses frequency-dependent scattering matrices, and (ii) a matrix product states (MPS) approach, which uses quantum noise operators in time bins and a tensor network technique to solve the time-dependent waveguide function for the entire system. Beginning with pulsed excitation using one-photon and two-photon Fock state pulses, we first show how to compute time-dependent observables with the scattering matrix approach, in terms of frequency integrals that encode the pulse spectrum, including how to extract the population dynamics of the excited quantum emitter, as well as the linear and nonlinear contributions. We present solutions for both symmetric and chiral TLS coupling. We then show how to compute the qubit and field observables in a more direct way using MPS, and obtain the characteristic bird-like shape for the two-photon correlation function at two times, which has been observed in recent experiments. We compare and contrast both of these methods, for one and two-photon excitation pulses, and show excellent agreement. We also present a study of the linear and nonlinear contributions, which can easily be calculated using scattering theory, and show the important role of pulse duration. Finally, we demonstrate the clear advantages of MPS by easily going to higher N-photon excitations, and show selected example population dynamics of up to eight-photon Fock-state pulses, manifesting in clear nonlinear population oscillations during the pulse interaction, similar to classical Rabi oscillations, but with quantum input fields that have a vanishing electric field expectation value.
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 detailed quantum theory and simulations of few-photon Fock state pulses interacting with a two-level system (TLS) in an open waveguide-QED system. For a rectangular pulse shape, we present an exact temporal scattering theory to derive analytical expressions for the TLS population, using one-photon and two-photon pulses, for both chiral and symmetric emitters. We also derive the stationary (long-time) and time-dependent spectra for one-photon excitation, and show how these differ at a fundamental level when considering TLS population effects. Numerically, we also present matrix product state (MPS) simulations, which allow us to compute more general photon correlation functions for arbitrary quantum pulses, and we use this approach to also show results for Gaussian quantum pulses and to confirm the accuracy of our analytical solutions. In the case of a chiral TLS, we show how a one-photon pulse, of any temporal shape, yields a transmitted long-time spectrum that is identical to the input pulse, despite significant TLS population effects. However, the dynamical spectra and spectral intensity show rich population effects, allowing for their detection in spectroscopic experiments. We also show the striking differences between one-photon and two-photon excitation, where the latter shows clear signatures of nonlinear saturation effects. In addition, we demonstrate how significant population TLS dynamics also occur for pulses that are relatively long compared to the radiative decay time (showing that a weak excitation approximation, which neglects finite TLS population effects, cannot be made) and investigate the population signatures, nonlinear features and dynamical behavior as a function of quantum pulse length.
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.