Computational framework for non-Markovian multi-emitter dynamics beyond the single-excitation limit
Pith reviewed 2026-05-13 20:15 UTC · model grok-4.3
The pith
A Green's function-based framework models non-Markovian multi-emitter dynamics in the two-excitation manifold by retaining photonic amplitudes explicitly.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
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.
What carries the argument
A non-Markovian hierarchy of coupled differential equations obtained from the dyadic Green's function, with explicit retention of photonic amplitudes inside the two-excitation manifold, implemented via the modified Langevin noise formalism and emitter-centered mode reduction.
If this is right
- Collective decay of symmetric Dicke states and selective stabilization become simulable in waveguides with embedded lossy structures.
- Enhanced Bell-state fidelity appears in selected homogeneous waveguide configurations under non-Markovian conditions.
- Entanglement sudden birth and oscillatory revivals are captured in the same structured-reservoir setting.
- The formulation applies in principle to any electromagnetic environment whose dyadic Green's function is available numerically.
Where Pith is reading between the lines
- The same hierarchy construction could be extended to three or more excitations provided sufficient computational resources are available.
- The retained phase information may prove useful for designing multi-photon quantum gates or sensors that operate in dispersive or lossy media.
- Hybridization with other open-system techniques could enlarge the accessible system sizes while preserving non-Markovian features.
Load-bearing premise
The dyadic Green's function of an arbitrary electromagnetic environment can be computed numerically, and the two-excitation manifold plus explicit photonic amplitudes are sufficient to describe the relevant phenomena without higher-order corrections.
What would settle it
A direct numerical or experimental comparison of predicted entanglement sudden birth and oscillatory revivals for symmetric Dicke states in a waveguide containing a lossy dielectric slab against an exact solution or measurement in the identical geometry.
Figures
read the original abstract
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.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents a Green's function-based computational framework for non-Markovian multi-emitter quantum electrodynamics restricted to the two-excitation manifold. It combines the modified Langevin noise (M-LN) formalism for first-principles treatment of dissipative environments with the emitter-centered mode (ECM) approach to maintain computational tractability. A hierarchy of coupled differential equations is constructed by explicitly retaining photonic amplitudes rather than integrating out the reservoir. The framework is claimed to preserve total probability and phase information. Numerical demonstrations are given for atom-field interactions in semi-infinite waveguides, including enhanced Bell-state fidelity in homogeneous waveguides, selective stabilization and delayed transfer in collective Dicke decay with an embedded lossy dielectric slab, and entanglement sudden birth with oscillatory revivals.
Significance. If the two-excitation truncation proves valid, the framework would provide a practical extension of non-Markovian studies to multi-photon regimes where Markovian approximations break down due to nonlinearities and phase accumulation. It enables first-principles modeling in arbitrary electromagnetic environments once the dyadic Green's function is available numerically, with explicit retention of photonic amplitudes supporting interference phenomena. This could be useful for designing quantum optical systems involving collective effects and structured reservoirs.
major comments (2)
- [Abstract / Numerical demonstrations] Abstract and numerical demonstrations section: The central claim that the two-excitation hierarchy captures all relevant non-Markovian multi-photon effects without higher-order corrections lacks supporting evidence. No convergence tests with manifold size, comparisons to full Hilbert-space numerics, or perturbative estimates of truncation error are provided for the reported Bell fidelity, Dicke decay rates, or entanglement revivals in the lossy dielectric slab geometry.
- [Abstract] Abstract: The assertion that the formulation preserves total probability and retains phase information is stated but not quantitatively validated through error analysis or benchmarks against exact methods for the structured-reservoir cases (e.g., waveguide with lossy slab).
minor comments (2)
- [Framework description] Clarify the precise definition and numerical implementation of the dyadic Green's function retrieval step, including any discretization or approximation details used in the waveguide examples.
- [Introduction / Methods] Add explicit references to prior M-LN and ECM literature when stating how the hierarchy is constructed independently.
Simulated Author's Rebuttal
Thank you for the opportunity to respond to the referee's report. We address the major comments point by point below, clarifying the scope of our framework and providing additional validation where feasible in the revision.
read point-by-point responses
-
Referee: [Abstract / Numerical demonstrations] Abstract and numerical demonstrations section: The central claim that the two-excitation hierarchy captures all relevant non-Markovian multi-photon effects without higher-order corrections lacks supporting evidence. No convergence tests with manifold size, comparisons to full Hilbert-space numerics, or perturbative estimates of truncation error are provided for the reported Bell fidelity, Dicke decay rates, or entanglement revivals in the lossy dielectric slab geometry.
Authors: We acknowledge the absence of explicit convergence tests in the original manuscript. Full numerical comparisons to the complete Hilbert space are intractable for the considered systems due to the continuous nature of the reservoir. In the revised version, we include perturbative estimates of higher-order corrections based on the low excitation density in our parameter regimes and add a dedicated paragraph discussing the validity of the two-excitation truncation for the reported observables. revision: partial
-
Referee: [Abstract] Abstract: The assertion that the formulation preserves total probability and retains phase information is stated but not quantitatively validated through error analysis or benchmarks against exact methods for the structured-reservoir cases (e.g., waveguide with lossy slab).
Authors: We have revised the manuscript to include quantitative validation: numerical results now show that the total probability is conserved up to an error of order 10^{-5} or better across all simulations. Phase retention is evidenced by the accurate capture of oscillatory revivals and interference patterns. While exact benchmarks for the full non-Markovian case are not feasible, we compare to limiting cases (e.g., Markovian approximation) to support the claims. revision: yes
Circularity Check
Derivation chain is self-contained with no circular reductions
full rationale
The paper constructs its non-Markovian hierarchy by starting from the dyadic Green's function (assumed obtainable numerically) and the established M-LN and ECM formalisms cited from prior literature, then explicitly retaining photonic amplitudes to form coupled differential equations within the two-excitation manifold. No equation is defined in terms of its own output, no fitted parameter is relabeled as a prediction, and no load-bearing uniqueness theorem or ansatz is imported solely via self-citation. The probability preservation and phase retention follow directly from the explicit retention of amplitudes in the hierarchy, independent of the numerical demonstrations. The two-excitation truncation is stated as an assumption rather than derived, so it does not create circularity.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The modified Langevin noise (M-LN) formalism provides first-principles treatment of dissipative environments.
- domain assumption The emitter-centered mode (ECM) framework ensures computational tractability for the multi-emitter hierarchy.
Reference graph
Works this paper leans on
-
[1]
0.2474 Effective Schmidt numberK eff 4.0418 von Neumann entropyS vN (nats) 1.6891 von Neumann entropyS vN (bits) 2.4368 whereαandβdenote the combined channel- frequency indices. The eigenvalue distribution ofρ 1 determines the Schmidt-mode content of the bipho- ton state: a nearly separable state would yield one dominant eigenvalue, whereas a broader dist...
- [2]
-
[3]
M. S. Tame, K. R. McEnery, K. ¨Ozdemir, J. Lee, S. A. Maier, and M. S. Kim, Quantum plasmonics, Nature Physics9, 329–340 (2013)
work page 2013
- [4]
-
[5]
M. Ye, Y. Tian, J. Lin, Y. Luo, J. You, J. Hu, W. Zhang, W. Chen, and X. Li, Universal quantum optimization with cold atoms in an optical cavity, Phys. Rev. Lett.131, 103601 (2023)
work page 2023
-
[6]
A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakinskiy, and A. N. Poddubny, Waveguide quantum electrodynamics: Collective radiance and photon-photon correlations, Rev. Mod. Phys.95, 015002 (2023)
work page 2023
-
[7]
S. Mahmoodian, M. ˇCepulkovskis, S. Das, P. Lo- dahl, K. Hammerer, and A. S. Sørensen, Strongly correlated photon transport in waveguide quan- tum electrodynamics with weakly coupled emitters, Phys. Rev. Lett.121, 143601 (2018)
work page 2018
-
[8]
S. Slussarenko and G. J. Pryde, Photonic quantum information processing: A concise review, Applied Physics Reviews6, 041303 (2019)
work page 2019
- [9]
-
[10]
A manufacturable platform for photonic quantum computing, Nature641, 876–883
-
[11]
A. Gonz´ alez-Tudela and J. I. Cirac, Markovian and non-markovian dynamics of quantum emitters cou- pled to two-dimensional structured reservoirs, Phys. Rev. A96, 043811 (2017)
work page 2017
-
[12]
I. de Vega and D. Alonso, Dynamics of non- markovian open quantum systems, Rev. Mod. Phys. 89, 015001 (2017)
work page 2017
-
[13]
M. O. Scully and M. S. Zubairy,Quantum Optics (Cambridge University Press, 1997)
work page 1997
-
[14]
L. Mandel and E. Wolf,Optical Coherence and Quantum Optics(Cambridge University Press, 1995)
work page 1995
- [15]
- [16]
- [17]
-
[18]
Jin,The Finite Element Method in Electro- magnetics, 3rd ed
J.-M. Jin,The Finite Element Method in Electro- magnetics, 3rd ed. (Wiley-IEEE Press, Hoboken, NJ, 2014)
work page 2014
-
[19]
A. Taflove and S. Hagness,Computational electro- dynamics: the finite-difference time-domain method. 2nd ed, Vol. 67–106 (2000)
work page 2000
-
[20]
W. C. Chew, B. Zhang, and J. Zhu, Some selected unsolved problems in classical and quantum electro- magnetics, Progress In Electromagnetics Research 180, 79–87 (2024)
work page 2024
-
[21]
S. Moon, G. Khan, S. T. Elkin, C. J. Ryu, D.-Y. Na, and T. E. Roth, Computational quantum elec- tromagnetics: Basic concepts and emerging trends, IEEE Antennas and Propagation Magazine , 2 (2025)
work page 2025
-
[22]
S. T. Elkin, G. Khan, E. Forati, B. W. Langley, D. Timucin, R. Molavi, S. Sussman, and T. E. Roth, Opportunities and challenges of computa- tional electromagnetics methods for superconduct- ing circuit quantum device modeling: A practical 20 review (2025), arXiv:2511.20774 [quant-ph]
-
[23]
T. E. Roth and W. C. Chew, Macroscopic circuit quantum electrodynamics: A new look toward devel- oping full-wave numerical models, IEEE Journal on Multiscale and Multiphysics Computational Tech- niques6, 109 (2021)
work page 2021
-
[24]
T. E. Roth and W. C. Chew, Full-wave methodology to compute the spontaneous emission rate of a trans- mon qubit, IEEE Journal on Multiscale and Multi- physics Computational Techniques7, 92 (2022)
work page 2022
-
[25]
D.-Y. Na, J. Zhu, and W. C. Chew, Diagonaliza- tion of the hamiltonian for finite-sized dispersive media: Canonical quantization with numerical mode decomposition, Phys. Rev. A103, 063707 (2021)
work page 2021
-
[26]
D.-Y. Na, T. E. Roth, J. Zhu, W. C. Chew, and C. J. Ryu, Numerical framework for modeling quantum electromagnetic systems involving finite-sized lossy dielectric objects in free space, Phys. Rev. A107, 063702 (2023)
work page 2023
-
[27]
S. Moon, D.-Y. Na, and T. E. Roth, Analytical quantum full-wave solutions for a 3-d circuit quan- tum electrodynamics system, IEEE Transactions on Antennas and Propagation72, 6702 (2024)
work page 2024
-
[28]
C. Forestiere and G. Miano, Operative approach to quantum electrodynamics in dispersive dielec- tric objects based on a polarization-mode expansion, Phys. Rev. A106, 033701 (2022)
work page 2022
-
[29]
C. Forestiere and G. Miano, Integral formulation of the macroscopic quantum electrodynamics in dis- persive dielectric objects, Phys. Rev. A107, 063705 (2023)
work page 2023
-
[30]
C. J. Ryu, D.-Y. Na, W. C. Chew, and E. Kudeki, Efficient tensor-network simulation for the few-atom multimode dicke model via coupling-matrix trans- formation, Phys. Rev. A108, 043707 (2023)
work page 2023
-
[31]
C. J. Ryu, D.-Y. Na, and W. C. Chew, Matrix prod- uct states and numerical mode decomposition for the analysis of gauge-invariant cavity quantum elec- trodynamics, Phys. Rev. A107, 063707 (2023)
work page 2023
-
[32]
C. Huang, H. Ge, Y. Cheng, Z. He, F. Liu, and W. E. I. Sha, Modeling of far-field quantum co- herence by dielectric bodies based on the volume integral equation method (2025), arXiv:2508.16471 [quant-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2025
- [33]
-
[34]
H. Choi, T. E. Roth, W. C. Chew, and D.-Y. Na, Non-markovian analysis of atom-field interactions in dissipative electromagnetic environments, Phys. Rev. Appl.24, 044056 (2025)
work page 2025
-
[35]
H. Choi, J. Seo, W. C. Chew, and D.-Y. Na, Atom- field non-markovian dynamics in open and dissipa- tive systems: An efficient memory-kernel approach linked to dyadic greens function and cem treatments (2025)
work page 2025
- [36]
-
[37]
V. I. Yudson and P. Reineker, Multiphoton scatter- ing in a one-dimensional waveguide with resonant atoms, Phys. Rev. A78, 052713 (2008)
work page 2008
-
[38]
A. Reiserer, N. Kalb, G. Rempe, and S. Ritter, A quantum gate between a flying optical photon and a single trapped atom, Nature508, 237–240 (2014)
work page 2014
-
[39]
J.-T. Shen and S. Fan, Strongly correlated two- photon transport in a one-dimensional waveguide coupled to a two-level system, Phys. Rev. Lett.98, 153003 (2007)
work page 2007
-
[40]
D. Roy, Two-photon scattering by a driven three- level emitter in a one-dimensional waveguide and electromagnetically induced transparency, Phys. Rev. Lett.106, 053601 (2011)
work page 2011
-
[41]
W. Gu, T. Li, Y. Tian, Z. Yi, and G.-x. Li, Two- photon dynamics in non-markovian waveguide qed with a giant atom, Phys. Rev. A110, 033707 (2024)
work page 2024
-
[42]
W. Xie, I. M. Mirza, and J. C. Schotland, Two- photon superradiance and subradiance, Phys. Rev. A112, 043712 (2025)
work page 2025
-
[43]
O. D. Stefano, S. Savasta, and R. Girlanda, Mode expansion and photon operators in dispersive and absorbing dielectrics, Journal of Modern Optics48, 67 (2001)
work page 2001
-
[44]
Drezet, Quantizing polaritons in inhomogeneous dissipative systems, Phys
A. Drezet, Quantizing polaritons in inhomogeneous dissipative systems, Phys. Rev. A95, 023831 (2017)
work page 2017
-
[45]
A. Ciattoni, Quantum electrodynamics of lossy magnetodielectric samples in vacuum: Modified langevin noise formalism, Phys. Rev. A110, 013707 (2024)
work page 2024
-
[46]
A. Ciattoni, Direct derivation of the modified langevin noise formalism from the canonical quan- tization of macroscopic electromagnetism (2026), arXiv:2603.04336 [quant-ph]
-
[47]
T. Gruner and D.-G. Welsch, Green-function ap- proach to the radiation-field quantization for ho- mogeneous and inhomogeneous kramers-kronig di- electrics, Phys. Rev. A53, 1818 (1996)
work page 1996
-
[48]
H. T. Dung, L. Kn¨ oll, and D.-G. Welsch, Sponta- neous decay in the presence of dispersing and ab- sorbing bodies: General theory and application to a spherical cavity, Phys. Rev. A62, 053804 (2000)
work page 2000
-
[49]
S. Scheel and S. Y. Buhmann, Macroscopic quantum electrodynamics – concepts and applications, Acta Physica Slovaca58, 675 (2008)
work page 2008
- [50]
-
[51]
M. S´ anchez-Barquilla, F. J. Garc´ ıa-Vidal, A. I. Fern´ andez-Dom´ ınguez, and J. Feist, Few-mode field quantization for multiple emitters, Nanophotonics 11, 4363 (2022)
work page 2022
- [52]
- [53]
- [54]
-
[55]
J.-C. Zheng, X.-W. Zheng, X.-L. Hei, Y.-F. Qiao, X.-Y. Yao, X.-F. Pan, Y.-M. Ren, X.-W. Huo, and P.-B. Li, Non-markovian dynamics withλ-type atomic systems in a single-end photonic waveg- uide, Advanced Quantum Technologies8, 2500133 (2025)
work page 2025
-
[56]
L. Mazzola, S. Maniscalco, J. Piilo, K.-A. Suomi- nen, and B. M. Garraway, Sudden death and sudden birth of entanglement in common structured reser- voirs, Phys. Rev. A79, 042302 (2009)
work page 2009
-
[57]
H. T. Dung, L. Kn¨ oll, and D.-G. Welsch, Three- dimensional quantization of the electromagnetic field in dispersive and absorbing inhomogeneous di- electrics, Phys. Rev. A57, 3931 (1998)
work page 1998
-
[58]
A. Miguel-Torcal, A. Gonz´ alez-Tudela, F. J. Garc´ ıa- Vidal, and A. I. Fern´ andez-Dom´ ınguez, Photon- mediated interactions and dynamics of coherently driven quantum emitters in complex photonic envi- ronments, Phys. Rev. B112, 235419 (2025)
work page 2025
-
[59]
R. H. Dicke, Coherence in spontaneous radiation processes, Phys. Rev.93, 99 (1954)
work page 1954
- [60]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.