Atom-Field Non-Markovian Dynamics in Open and Dissipative Systems: An Efficient Memory-Kernel Approach Linked to Dyadic Greens Function and CEM Treatments

TL;DR

This paper introduces a computational framework that models non-Markovian atom-field interactions in open systems using Green's functions, compatible with standard electromagnetic solvers, enabling accurate quantum emitter simulations.

Contribution

It develops a Green's function-based method for non-Markovian dynamics, integrated with FDTD and FEM, advancing quantum light-matter interaction modeling in classical electromagnetic simulations.

Findings

Reconstructed Greens function via modal expansion in 3D.

Memory kernel determined by Greens function's imaginary part.

Validated approach with lossy mirror and cavity scenarios.

Abstract

In this work, we present a numerical framework for modeling single photon emission from a two level system in open and dissipative systems beyond the Markovian approximation. The method can be readily integrated into standard computational electromagnetic (CEM) solvers such as finite difference time domain (FDTD) and finite element method (FEM). We numerically verify the completeness of boundary and medium assisted modes in the modified Langevin noise formalism by reconstructing the imaginary part of the dyadic Greens function through modal expansion in three dimensions. This reconstruction enables a first principles description of atom field interaction via the multi mode Jaynes Cummings model in open and dissipative environments. Within the single excitation manifold, we show that the memory kernel of a two level system is determined by the imaginary part of the Greens function, implying that radiative modes alone govern the relevant dynamics. The proposed framework thus provides a Greens function based approach for describing atomic population and single photon dynamics, directly compatible with Maxwell solvers. We then present concrete strategies for implementing our method in both FDTD and FEM frameworks, demonstrating its practical applicability. We further verify numerical results for a lossy Lorentz Drude type mirror, including both the case of a TLS near a finite sized metallic mirror and that of a TLS centered in a Fabry Perot cavity. This work establishes a rigorous foundation for incorporating quantum emitter dynamics into computational electromagnetics, thereby extending classical solvers toward quantum light matter interactions.