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.
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.
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.
Strong coupling of quantum emitters with confined electromagnetic modes of nanophotonic structures may be used to change optical, chemical, and transport properties of materials, with significant theoretical effort invested toward a better understanding of this phenomenon. However, a full theoretical description of both matter and light is an extremely challenging task. Typical theoretical approaches simplify the description of the photonic environment by describing it as a single mode or few modes. While this approximation is accurate in some cases, it breaks down strongly in complex environments, such as within plasmonic nanocavities, and the electromagnetic environment must be fully taken into account. This requires the quantum description of a continuum of bosonic modes, a problem that is computationally hard. We here investigate a compromise where the quantum character of light is taken into account at modest computational cost. To do so, we focus on a quantum emitter that interacts with an arbitrary photonic spectral density and employ the cumulant, or cluster, expansion method to the Heisenberg equations of motion up to first, second, and third order. We benchmark the method by comparing it with exact solutions for specific situations and show that it can accurately represent dynamics for many parameter ranges.
We develop an exact framework for describing the non-Markovian dynamics of an open quantum system interacting with an environment modeled by a generalized spectral density function.The approach relies on mapping the initial system onto an auxiliary configuration, comprising the original open system coupled to a small number of discrete modes, which in turn are each coupled to an independent Markovian reservoir.Based on the connection between the discrete modes and the poles of the spectral density function, we show how expanding the system using the discrete modes allows for the full inclusion non-Markovian effects within an enlarged open system whose dynamics is governed by an exact Lindblad master equation.Initially we apply this result to obtain a generalization of the pseudomode method [B.M. Garraway, Phys. Rev. A. 55, 2290 (1997)] in cases where the spectral density function has a Lorentzian structure.For many other types of spectral density function, we extend our proof to show that an open system dynamics may be modeled physically using discrete modes which admit a non-Hermitian coupling to the system and for such cases determine the equivalent master equation to no longer be of Lindblad form.For applications involving two discrete modes, we demonstrate how to convert between pathological and Lindblad forms of the master equation using the techniques of the pseudomode method.
In recent years, Born-Markov master equations based on tracing out the electromagnetic degrees of freedom have been extensively employed in the description of quantum optical phenomena originating from photon-mediated interactions in quantum emitter ensembles. The breakdown of these effective models, built upon assumptions such as ensemble spectral homogeneity, an unstructured photonic density of states, and weak light-matter coupling, has also recently attracted considerable attention. Here, we investigate the accuracy of this well-established framework beyond the most conventional, and extensively explored, spontaneous emission configuration. Specifically, we consider a system comprising two coherently driven and detuned quantum emitters, embedded within a hybrid photonic-plasmonic cavity, formed by a metallic nanorod integrated into a high-refractive-index dielectric microresonator. The local density of photonic states in this structure exhibits a complex frequency dependence, making it a compelling platform for exploring photon-mediated interactions beyond the assumptions above. We benchmark this modeling approach for the quantum dynamics of the emitter pair against exact calculations based on a macroscopic field quantization formalism, providing an illustrative assessment of its validity in significantly structured and dispersive photonic environments. Our analysis reveals four distinct regimes of laser driving and frequency splitting that lead to markedly different levels of accuracy in the effective model.
We provide a simple semi-classical formalism to describe the coupling between one or several quantum emitters and a structured environment. Describing the emitter by an electric polarizability, and the surrounding medium by a Green function, we show that an intuitive scattering picture allows one to derive a coupling equation from which the eigenfrequencies of the coupled system can be extracted. The model covers a variety of regimes observed in light-matter interaction, including weak and strong coupling, coherent collective interactions, and incoherent energy transfer. It provides a unified description of many processes, showing that different interaction regimes are actually rooted on the same ground. It can also serve as a basis for the development of more refined models in a full quantum electrodynamics framework.
While non-Markovian dynamics have been extensively studied in the single-excitation limit to predict non-trivial phenomena, this regime remains an idealization. Moving beyond it is essential, as optical nonlinearities and phase-error accumulation in multi-photon processes render the Markovian approximation fragile. In this work, we present a Green's function-based framework for modeling non-Markovian multi-emitter quantum electrodynamics within the two-excitation manifold. The modified Langevin noise (M-LN) formalism is employed for first-principles treatment of dissipative environments, while the emitter-centered mode (ECM) framework ensures computational tractability. Unlike conventional approaches that integrate out the reservoir, we construct a non-Markovian hierarchy of coupled differential equations by explicitly retaining photonic amplitudes. Within the two-excitation hierarchy, the formulation preserves total probability and retains phase information necessary to capture multi-photon interference. As numerical demonstrations, we investigate non-Markovian atom-field interactions in structured semi-infinite waveguide environments. We first consider a homogeneous waveguide as a baseline, observing enhanced Bell-state fidelity in selected configurations. Next, we examine collective decay of symmetric Dicke states in a waveguide with an embedded lossy dielectric slab, revealing selective stabilization and delayed excitation transfer induced by the structured reservoir. Finally, we analyze entanglement dynamics in the same setting, highlighting entanglement sudden birth and oscillatory revivals. In principle, the framework applies to arbitrary electromagnetic environments for which the dyadic Green's function can be obtained numerically, providing a versatile tool for investigating complex non-Markovian multi-photon phenomena beyond the single-excitation limit.