Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2302.06145 · Physical Review A

Numerical framework for modeling quantum electromagnetic systems involving finite-sized lossy dielectric objects in free space

Dong-Yeop Na, Thomas E. Roth, Jie Zhu, Weng Cho Chew, Christopher Jayun Ryu

Published Jun 8, 2023• 16 citations
arXivHTMLPDF

Abstract

The modified Langevin noise formalism [A. Drezet, Phys. Rev. A 95, 023831 (2017); O. D. Stefano, S. Savasta, and R. Girlanda, J. Mod. Opt. 48, 67 (2001)] has been proposed for the correct characterization of quantum electromagnetic fields in the presence of finite-size lossy dielectric objects in free space. The main modification to the original one [T. Gruner and D.-G. Welsch, Phys. Rev. A 53, 1818 (1996); H. T. Dung, L. Kn\"oll, and D.-G. Welsch, Phys. Rev. A 57, 3931 (1998)] (also known as the Green's function approach for a bulk inhomogeneous lossy dielectric medium) was to introduce another fluctuating source in reaction to the radiation loss. Consequently, the resulting electric field operator is now determined by (i) boundary-assisted and (ii) medium-assisted fields on an equal footing, which originate from radiation and medium losses, respectively. However, due to the lengthy mathematical manipulation and complicated concepts, the validity of the modified Langevin noise formalism has not been clearly confirmed yet. In this work, we propose and develop a numerical framework for the modified Langevin noise formalism by exploiting computational electromagnetic methods. Specifically, we present utilization of the finite-element method to numerically solve plane-wave-scattering and point-source-radiation problems whose solutions are boundary-assisted and medium-assisted fields, respectively. We numerically validated the modified Langevin noise model calculating the Purcell factor of a two-level atom inside or outside a lossy dielectric slab. The proposed numerical framework is particularly useful for analyzing the dynamics of multilevel atoms near plasmonic structures or metasurfaces in the open space.

DOI
10.1103/physreva.107.063702

Related papers

Numerical Framework for Multimode Jaynes- and Tavis-Cummings Models Incorporating the Modified Langevin Noise Formalism: Non-Markovian Analysis of Atom-Field Interactions in Dissipative Electromagnetic Environments

Choi, Hyunwoo, Weng Cho Chew, Dong-Yeop Na

We present a novel numerical framework that integrates the modified Langevin noise formalism into the multimode Jaynes- and Tavis-Cummings models, enabling a first-principles, non-Markovian analysis of atom-field interactions in dissipative electromagnetic (EM) environments that account for both radiative losses and absorptive dissipation in lossy dielectric media (or satisfying general inhomogeneous causal media exhibiting both dispersion and absorption effects). In the modified Langevin noise formalism, the boundary- and medium-assisted (BA and MA) fields, which constitute a continuum set of EM modes in dissipative EM environments, are numerically obtained using the finite-element method (FEM). Specifically, BA field modes are extracted by solving plane-wave scattering problems, while MA field modes are determined through point-source radiation problems. These numerically obtained BA and MA field modes are then incorporated into the multimode Jaynes- and Tavis-Cummings models such that the coupling strength between atoms and BA-MA field modes can be calculated for the study of atom-field interactions in dissipative EM environments. The proposed methodology captures non-Markovian atomic dynamics that cannot be described by traditional quantum master equations under the Markovian approximation. To validate the accuracy of the proposed numerical framework, we present four numerical examples: (i) a two-level system (TLS) in a perfect electric conductor (PEC) half-space; (ii) dissipative cavity electrodynamics with two limiting cases approaching spontaneous emission in free space and ideal Rabi oscillations; (iii) super-radiance in TLS arrays; and (iv) entanglement sudden death of two TLSs inside dissipative cavities. The proposed methodology can serve as a ground-truth numerical simulator for studying atom-field interactions in general dissipative EM environments.

arXiv:2504.02183•Apr 3, 2025•0 citations

Non-Markovian analysis of atom-field interactions in dissipative electromagnetic environments

Hyunwoo Choi, Thomas E. Roth, Weng Cho Chew +1 more

We present a numerical framework that integrates a modified Langevin noise formalism into multimode Jaynes- and Tavis-Cummings models for first-principles, non-Markovian analysis of atom-field interactions in dissipative electromagnetic environments under weak coupling. This formalism enables the quantization of electric field operators in terms of boundary- and medium-assisted modes with real-valued eigenfrequencies, accounting for both radiative and absorptive losses. To efficiently handle the continuum nature of the boundary- and medium-assisted mode spectrum, we apply adaptive frequency refinement and mode-space coarse-graining based on finite-element simulations to obtain a tractable discrete set of field modes. These modes are then used to construct matrix representations of the Schr\"odinger equation, enabling direct computation of atomic observables and field quantities---such as single-photon amplitudes and higher-order correlation functions---beyond the Markovian approximation. Importantly, this framework does not rely on any phenomenologically designed memory kernel. Instead, the spectrum of boundary- and medium-assisted modes interacting with atoms inherently and rigorously captures memory effects arising from the open and dissipative electromagnetic environment. We present four numerical examples to validate the proposed framework and demonstrate its ability to model non-Markovian atomic dynamics in open and dissipative electromagnetic environments: (i) a two-level system near the surface of a perfect electric conductor, (ii) dissipative cavity quantum electrodynamics for a single two-level system, (iii) super- and subradiance in two-level system arrays, and (iv) entanglement sudden death between two two-level systems. This framework provides a versatile simulation tool for the design of quantum photonic devices, including nitrogen-vacancy center single-photon sources coupled to plasmonic structures, entangled photon-pair generation near lossy media, and decoherence modeling of solid-state quantum memories based on rare-earth ions and color centers. It is particularly well suited to dissipative quantum electrodynamics in photonic integrated circuits, where accurate modeling of scattering and diffraction is essential.

Physical Review Applied•Sep 22, 2025•1 citation

Quantum emitter interacting with a dispersive dielectric object: a model based on the modified Langevin noise formalism

Giovanni Miano, Loris Maria Cangemi, Carlo Forestiere

In this paper, we model the interaction of a quantum emitter with a finite-size dispersive dielectric object in an unbounded space within the framework of macroscopic quantum electrodynamics, using the modified Langevin noise formalism, without any restrictions on the emitter level structure or dipole operator. The quantized electromagnetic field consists of two contributions: the medium-assisted field, which accounts for the electromagnetic field generated by the noise polarization currents of the dielectric, and the scattering-assisted field, which takes into account the electromagnetic field incoming from infinity and scattered by the dielectric. We show that the emitter couples to two distinct bosonic baths: a medium-assisted bath and a scattering-assisted bath, each characterized by its own spectral density. We identify the conditions under which the electromagnetic environment composed of these two baths can be effectively replaced by a single bosonic bath, so that the reduced dynamics of the quantum emitter remain unchanged. In particular, when the initial states of the medium- and scattering-assisted baths are thermal states with the same temperature, we find that a single bosonic bath with a spectral density equal to the sum of the medium-assisted and scattering-assisted spectral densities is equivalent to the original electromagnetic environment.

arXiv:2410.10730•Oct 14, 2024•0 citations

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

Hyunwoo Choi, Jisang Seo, Weng C. Chew +1 more

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.

quant-pharXiv:2511.03561•Nov 6, 2025•0 citations

Quantum emitter interacting with a dispersive dielectric object: a model based on the modified Langevin noise formalism

Giovanni Miano, Loris Maria Cangemi, Carlo Forestiere

In this paper, we model the interaction of a quantum emitter with a finite-size dispersive dielectric object in unbounded space within the framework of macroscopic quantum electrodynamics, using the modified Langevin noise formalism. The quantized electromagnetic field consists of two contributions: the medium-assisted field, which accounts for the electromagnetic field generated by the noise polarization currents of the dielectric, and the scattering-assisted field, which takes into account the electromagnetic field incoming from infinity and scattered by the dielectric. We show that the emitter couples to two distinct bosonic reservoirs: a medium-assisted reservoir and a scattering-assisted reservoir, each characterized by its own spectral density. We then use emitter-centered modes to reduce the degrees of freedom of both reservoirs. Eventually, we identify the conditions under which the electromagnetic environment composed of these two reservoirs can be effectively replaced by a single bosonic reservoir so that the reduced time evolution of the quantum emitter remains unchanged. In particular, when the initial states of the medium- and scattering-assisted reservoirs are thermal quantum states at the same temperature, we find that a single bosonic bath with a spectral density equal to the sum of the medium- and scattering-assisted spectral densities is equivalent to the original electromagnetic environment.

Nanophotonics•Mar 28, 2025•7 citations

Universal and rigorous asymptotic convergence of standard macroscopic quantum electrodynamics to the modified Langevin noise formalism

A. Ciattoni

The standard Langevin noise formalism (LNF) is a well-established approach for macroscopic quantum electrodynamics in unbounded lossy media, and accordingly, its application to finite objects requires endowing the surrounding vacuum with a fictitious absorption that is taken to asymptotically vanish. Existing justifications for such a singular limiting procedure have been mainly restricted to phenomena governed by vacuum fluctuations and medium-induced interactions (Casimir-Polder forces, thermal radiation, quantum emitter dynamics, etc.), thereby strictly excluding scenarios that involve the scattering of external quantum radiation. Conversely, the modified Langevin noise formalism (MLNF) has recently been established as the complete canonical quantum theory for arbitrary macroscopic systems and electromagnetic processes, entirely bypassing the need for limiting procedures, treating incident quantum radiation fields as independent degrees of freedom, and exactly reducing to the LNF for unbounded media. In this work, we demonstrate that the LNF limiting procedure is universally valid by showing that it recovers any physical prediction of the MLNF across all regimes, including arbitrary quantum scattering processes. Specifically, we introduce a suitable mapping associating any MLNF state and observable with a corresponding entity within the LNF augmented by fictitious vacuum absorption. We then rigorously prove that, as the fictitious absorption tends to zero, every LNF matrix element constructed via this mapping exactly converges to its MLNF original, so that the limiting procedure establishes the asymptotic convergence of the LNF to the MLNF.

quant-pharXiv:2610.05202•Oct 6, 2026
The modified Langevin noise formalism [A. Drezet, Phys. Rev. A 95, 023831 (2017); O. D. Stefano, S. Savasta, and R. Girlanda, J. Mod. Opt. 48, 67 (2001)] has been proposed for the correct characterization of quantum electromagnetic fields in the presence of finite-size lossy dielectric objects in free space. The main modification to the original one [T. Gruner and D.-G. Welsch, Phys. Rev. A 53, 1818 (1996); H. T. Dung, L. Kn"oll, and D.-G. Welsch, Phys. Rev. A 57, 3931 (1998)] (also known as the Green's function approach for a bulk inhomogeneous lossy dielectric medium) was to introduce another fluctuating source in reaction to the radiation loss. Consequently, the resulting electric field operator is now determined by (i) boundary-assisted and (ii) medium-assisted fields on an equal footing, which originate from radiation and medium losses, respectively. However, due to the lengthy mathematical manipulation and complicated concepts, the validity of the modified Langevin noise formalism has not been clearly confirmed yet. In this work, we propose and develop a numerical framework for the modified Langevin noise formalism by exploiting computational electromagnetic methods. Specifically, we present utilization of the finite-element method to numerically solve plane-wave-scattering and point-source-radiation problems whose solutions are boundary-assisted and medium-assisted fields, respectively. We numerically validated the modified Langevin noise model calculating the Purcell factor of a two-level atom inside or outside a lossy dielectric slab. The proposed numerical framework is particularly useful for analyzing the dynamics of multilevel atoms near plasmonic structures or metasurfaces in the open space.
We present a novel numerical framework that integrates the modified Langevin noise formalism into the multimode Jaynes- and Tavis-Cummings models, enabling a first-principles, non-Markovian analysis of atom-field interactions in dissipative electromagnetic (EM) environments that account for both radiative losses and absorptive dissipation in lossy dielectric media (or satisfying general inhomogeneous causal media exhibiting both dispersion and absorption effects). In the modified Langevin noise formalism, the boundary- and medium-assisted (BA and MA) fields, which constitute a continuum set of EM modes in dissipative EM environments, are numerically obtained using the finite-element method (FEM). Specifically, BA field modes are extracted by solving plane-wave scattering problems, while MA field modes are determined through point-source radiation problems. These numerically obtained BA and MA field modes are then incorporated into the multimode Jaynes- and Tavis-Cummings models such that the coupling strength between atoms and BA-MA field modes can be calculated for the study of atom-field interactions in dissipative EM environments. The proposed methodology captures non-Markovian atomic dynamics that cannot be described by traditional quantum master equations under the Markovian approximation. To validate the accuracy of the proposed numerical framework, we present four numerical examples: (i) a two-level system (TLS) in a perfect electric conductor (PEC) half-space; (ii) dissipative cavity electrodynamics with two limiting cases approaching spontaneous emission in free space and ideal Rabi oscillations; (iii) super-radiance in TLS arrays; and (iv) entanglement sudden death of two TLSs inside dissipative cavities. The proposed methodology can serve as a ground-truth numerical simulator for studying atom-field interactions in general dissipative EM environments.
We present a numerical framework that integrates a modified Langevin noise formalism into multimode Jaynes- and Tavis-Cummings models for first-principles, non-Markovian analysis of atom-field interactions in dissipative electromagnetic environments under weak coupling. This formalism enables the quantization of electric field operators in terms of boundary- and medium-assisted modes with real-valued eigenfrequencies, accounting for both radiative and absorptive losses. To efficiently handle the continuum nature of the boundary- and medium-assisted mode spectrum, we apply adaptive frequency refinement and mode-space coarse-graining based on finite-element simulations to obtain a tractable discrete set of field modes. These modes are then used to construct matrix representations of the Schr\"odinger equation, enabling direct computation of atomic observables and field quantities---such as single-photon amplitudes and higher-order correlation functions---beyond the Markovian approximation. Importantly, this framework does not rely on any phenomenologically designed memory kernel. Instead, the spectrum of boundary- and medium-assisted modes interacting with atoms inherently and rigorously captures memory effects arising from the open and dissipative electromagnetic environment. We present four numerical examples to validate the proposed framework and demonstrate its ability to model non-Markovian atomic dynamics in open and dissipative electromagnetic environments: (i) a two-level system near the surface of a perfect electric conductor, (ii) dissipative cavity quantum electrodynamics for a single two-level system, (iii) super- and subradiance in two-level system arrays, and (iv) entanglement sudden death between two two-level systems. This framework provides a versatile simulation tool for the design of quantum photonic devices, including nitrogen-vacancy center single-photon sources coupled to plasmonic structures, entangled photon-pair generation near lossy media, and decoherence modeling of solid-state quantum memories based on rare-earth ions and color centers. It is particularly well suited to dissipative quantum electrodynamics in photonic integrated circuits, where accurate modeling of scattering and diffraction is essential.
In this paper, we model the interaction of a quantum emitter with a finite-size dispersive dielectric object in an unbounded space within the framework of macroscopic quantum electrodynamics, using the modified Langevin noise formalism, without any restrictions on the emitter level structure or dipole operator. The quantized electromagnetic field consists of two contributions: the medium-assisted field, which accounts for the electromagnetic field generated by the noise polarization currents of the dielectric, and the scattering-assisted field, which takes into account the electromagnetic field incoming from infinity and scattered by the dielectric. We show that the emitter couples to two distinct bosonic baths: a medium-assisted bath and a scattering-assisted bath, each characterized by its own spectral density. We identify the conditions under which the electromagnetic environment composed of these two baths can be effectively replaced by a single bosonic bath, so that the reduced dynamics of the quantum emitter remain unchanged. In particular, when the initial states of the medium- and scattering-assisted baths are thermal states with the same temperature, we find that a single bosonic bath with a spectral density equal to the sum of the medium-assisted and scattering-assisted spectral densities is equivalent to the original electromagnetic environment.
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.
In this paper, we model the interaction of a quantum emitter with a finite-size dispersive dielectric object in unbounded space within the framework of macroscopic quantum electrodynamics, using the modified Langevin noise formalism. The quantized electromagnetic field consists of two contributions: the medium-assisted field, which accounts for the electromagnetic field generated by the noise polarization currents of the dielectric, and the scattering-assisted field, which takes into account the electromagnetic field incoming from infinity and scattered by the dielectric. We show that the emitter couples to two distinct bosonic reservoirs: a medium-assisted reservoir and a scattering-assisted reservoir, each characterized by its own spectral density. We then use emitter-centered modes to reduce the degrees of freedom of both reservoirs. Eventually, we identify the conditions under which the electromagnetic environment composed of these two reservoirs can be effectively replaced by a single bosonic reservoir so that the reduced time evolution of the quantum emitter remains unchanged. In particular, when the initial states of the medium- and scattering-assisted reservoirs are thermal quantum states at the same temperature, we find that a single bosonic bath with a spectral density equal to the sum of the medium- and scattering-assisted spectral densities is equivalent to the original electromagnetic environment.
The standard Langevin noise formalism (LNF) is a well-established approach for macroscopic quantum electrodynamics in unbounded lossy media, and accordingly, its application to finite objects requires endowing the surrounding vacuum with a fictitious absorption that is taken to asymptotically vanish. Existing justifications for such a singular limiting procedure have been mainly restricted to phenomena governed by vacuum fluctuations and medium-induced interactions (Casimir-Polder forces, thermal radiation, quantum emitter dynamics, etc.), thereby strictly excluding scenarios that involve the scattering of external quantum radiation. Conversely, the modified Langevin noise formalism (MLNF) has recently been established as the complete canonical quantum theory for arbitrary macroscopic systems and electromagnetic processes, entirely bypassing the need for limiting procedures, treating incident quantum radiation fields as independent degrees of freedom, and exactly reducing to the LNF for unbounded media. In this work, we demonstrate that the LNF limiting procedure is universally valid by showing that it recovers any physical prediction of the MLNF across all regimes, including arbitrary quantum scattering processes. Specifically, we introduce a suitable mapping associating any MLNF state and observable with a corresponding entity within the LNF augmented by fictitious vacuum absorption. We then rigorously prove that, as the fictitious absorption tends to zero, every LNF matrix element constructed via this mapping exactly converges to its MLNF original, so that the limiting procedure establishes the asymptotic convergence of the LNF to the MLNF.