Quantum plasmonics with N emitters: bright hybrid continuum selection
Pith reviewed 2026-05-09 22:09 UTC · model grok-4.3
The pith
N quantum emitters reduce the plasmon-polariton continuum to N independent one-dimensional modes, one per emitter.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The quantum plasmon-polariton field supported by a finite dielectric medium consists of two families of continuum modes, each with infinite degeneracy, obtained as deformations of the free electromagnetic field and the free medium via the Lippmann-Schwinger equations. The interacting part of this field with N emitters reduces to N non-degenerate one-dimensional continua, one per emitter. The representation in terms of a single hybrid continuum spectrum coincides exactly with the macroscopic Langevin model with bulk medium due to an exact compensation between one term in the coupling Hamiltonian and another term in a Green tensor identity.
What carries the argument
The bright-dark decomposition of the plasmon-polariton field, which isolates N non-degenerate one-dimensional continua (one per emitter) that carry all the interaction.
If this is right
- The multi-emitter system can be described by N separate one-dimensional bosonic continua rather than the original doubly infinite continuum.
- The effective Hamiltonian for any number of emitters coincides with the one obtained from the macroscopic Langevin approach for bulk media.
- The equivalence holds for finite media because the bright-dark decomposition and the Green tensor cancellation are independent of specific geometry.
- Emission spectra, coherence, and entanglement dynamics can be computed from the reduced set of N continua.
Where Pith is reading between the lines
- The same reduction may apply to other hybrid light-matter systems that support continuum modes, such as photonic crystals or waveguides with multiple atoms.
- Numerical simulations of larger emitter arrays become feasible because only N continua need to be discretized instead of the full field.
- One could measure the predicted cancellation by comparing the spontaneous emission rate of a single emitter versus two closely spaced emitters in a controlled dielectric environment.
Load-bearing premise
The separation into bright modes that couple to emitters and dark modes that do not remains valid for any finite dielectric medium, and the Green tensor identity produces exact cancellation independent of emitter positions or medium details.
What would settle it
A full numerical diagonalization of the Hamiltonian for two emitters placed near a small dielectric structure that reveals the effective interaction requires more than two independent one-dimensional continua would falsify the claimed reduction.
Figures
read the original abstract
We construct mode-selective effective models describing the interaction of the quantum plasmon-polariton field supported by a finite dielectric medium and one or several quantum emitters. The construction of the effective model is based on the decomposition of the field into bright modes relevant to the interaction with the emitters and dark modes, which do not interact with the emitters. We show that the quantum plasmon-polariton field can be represented equivalently by a double-continuum spectrum or by a single hybrid continuum spectrum for each emitter. The system of the electromagnetic field coupled to a finite medium is composed of two families of continuum modes, each of them with an infinite degeneracy. The two families are deformations of the free electromagnetic field and the free medium, induced by the interaction between them, as described by the Lippmann-Schwinger equations. We show that if there are $N$ emitters interacting with this plasmon-polariton field, the effective interaction involves a much smaller set of bosonic continuum modes: the interacting part of the continuum can be described by $N$ non-degenerate one-dimensional continua, one for each emitter. The representation of the interaction in terms of a single hybrid continuum spectrum coincides with the one within the macroscopic Langevin model with bulk medium. This coincidence is explained by an exact compensation of two terms, one in the coupling term of the Hamiltonian and the other one in a Green tensor identity.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript constructs mode-selective effective models for the interaction of a quantum plasmon-polariton field supported by a finite dielectric medium with one or several quantum emitters. It decomposes the field into bright modes (relevant to emitters) and dark modes (non-interacting) via Lippmann-Schwinger equations, shows that for N emitters the interacting part of the continuum reduces to N non-degenerate one-dimensional continua, and claims that the resulting hybrid continuum spectrum coincides exactly with the macroscopic Langevin model with bulk medium due to cancellation of one term from the emitter-field coupling Hamiltonian against another from a Green tensor identity.
Significance. If the central claims hold, the work would be significant for quantum plasmonics by reducing the effective description of N-emitter interactions from two families of infinitely degenerate continua to a minimal set of N independent 1D continua. It also provides an explicit microscopic-to-macroscopic bridge via the asserted exact compensation, which could streamline theoretical modeling and device design in hybrid quantum systems involving emitters and finite plasmonic structures.
major comments (1)
- [Derivation of the hybrid continuum spectrum and Green tensor identity (following the abstract statement of the N-emitter] The central claim that the hybrid continuum representation coincides with the macroscopic Langevin model via exact compensation of a Hamiltonian coupling term against a Green tensor identity, and that this holds independently of emitter positions and for finite media, is load-bearing but insufficiently demonstrated. The construction begins from the same continuum-mode framework, and the Lippmann-Schwinger deformation for finite dielectrics incorporates geometry-dependent boundary scattering in the Green tensor; no explicit verification shows the opposing terms remain equal for arbitrary positions without additional restrictions on the medium response.
minor comments (1)
- The distinction between the double-continuum spectrum and the single hybrid continuum spectrum would be clearer with an explicit schematic or table comparing the mode degeneracies before and after the bright/dark decomposition.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable feedback on our manuscript arXiv:2604.21560. We appreciate the recognition of the work's potential significance and address the major comment below.
read point-by-point responses
-
Referee: [Derivation of the hybrid continuum spectrum and Green tensor identity (following the abstract statement of the N-emitter] The central claim that the hybrid continuum representation coincides with the macroscopic Langevin model via exact compensation of a Hamiltonian coupling term against a Green tensor identity, and that this holds independently of emitter positions and for finite media, is load-bearing but insufficiently demonstrated. The construction begins from the same continuum-mode framework, and the Lippmann-Schwinger deformation for finite dielectrics incorporates geometry-dependent boundary scattering in the Green tensor; no explicit verification shows the opposing terms remain equal for arbitrary positions without additional restrictions on the medium response.
Authors: We thank the referee for identifying this as a load-bearing claim requiring clear demonstration. The Green tensor identity follows directly from the general resolvent properties of the electromagnetic operator in an arbitrary finite dielectric, as obtained from the Lippmann-Schwinger equation; this identity is independent of emitter locations. The position-dependent term in the emitter-field coupling Hamiltonian is constructed precisely to cancel the corresponding contribution in the Green tensor, yielding the hybrid continuum that matches the macroscopic Langevin model. This algebraic cancellation holds for any emitter positions and any finite geometry without further restrictions on the medium response. While the manuscript outlines the compensation, we agree that an explicit step-by-step verification would strengthen the presentation. In the revised manuscript we will add a dedicated appendix that performs the cancellation algebraically using the general form of the Green tensor, confirming the result for arbitrary positions. revision: yes
Circularity Check
Derivation via Lippmann-Schwinger deformation yields independent reduction to N continua with explicit cancellation to Langevin model
full rationale
The paper derives the bright/dark decomposition from the Lippmann-Schwinger equations applied to the coupled field-medium system, treating both families as deformations of the free electromagnetic and free medium continua. The reduction of the interacting continuum to exactly N non-degenerate one-dimensional continua follows directly from associating one bright continuum per emitter via the position-dependent coupling. The claimed coincidence with the macroscopic Langevin model is obtained by exhibiting an exact algebraic cancellation between one term in the emitter-field Hamiltonian and a term in the Green-tensor identity; this cancellation is presented as a mathematical identity that holds after the deformation step and does not rely on fitting parameters, self-referential definitions, or prior results by the same authors. No load-bearing step equates a derived quantity to an input by construction, and the central claims remain self-contained against the stated assumptions.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The system of the electromagnetic field coupled to a finite medium is composed of two families of continuum modes, each with infinite degeneracy, as deformations of the free electromagnetic field and the free medium induced by their interaction via the Lippmann-Schwinger equations.
Reference graph
Works this paper leans on
-
[1]
The phase of Ω m ν (r0) can be chosen arbi- trarily
DBM decomposition for the m-component The term−d eg · ˆEm+(r0) in the interaction part−ˆσx ⊗ deg · ˆEm+(r0) of ˆH m (23b) can be written in the form: −deg · ˆEm+(r0) =ℏ Z +∞ 0 dν X dµ gm ν,dµ(r0) ˆC m ν,dµ =ℏ Z +∞ 0 dνΩ m ν (r0)ˆbm ν (r0) (24) with thebright annihilation operatorsassociated with the m-components defined as ˆbm ν (r0) = X dµ hm ν,dµ(r0) ˆC...
-
[2]
DBM decomposition for the e-component The term−d eg · ˆEe+(r0) from the interaction part −ˆσx ⊗d eg · ˆEe+(r0) of ˆH e (23a) can be written in the form −deg · ˆEe+(r0) =ℏ Z +∞ 0 dω c X dκ ge ω,dκ(r0) ˆC e ω,dκ =ℏ Z +∞ 0 dωΩ e ω(r0)ˆbe ω(r0).(33) One can construct similarly as before a normalizedbright annihilation operatorfor the e-component, that sums ov...
-
[3]
Effective Hamiltonian with the double continuum After finding the components of the effective Hamilto- nian related to the m-components and e-components, we can construct the effective component of the full Hamil- tonian. By replacing the e- and m-terms in (22) with ˆH e → ˆH e eff and ˆH m → ˆH m eff and substituting (32) and (39), the effective Hamilton...
-
[4]
i =ℏ NX k=1 ˆσ(k) x ⊗ Z +∞ 0 dω × X dκ ge(k) ω,dκ ˆC e ω,dκ + X dµ gm(k) ω,dµ ˆC m ω,dµ + h.c
Double-continuum bright operators Keeping the double-continuum structure of the Hamil- tonian (59), we rewrite its interaction part by substitut- ing the electric-field operator (21) ˆHint =− NX k=1 ˆσ(k) x ⊗ h d(k) eg · ˆE+(rk) + h.c. i =ℏ NX k=1 ˆσ(k) x ⊗ Z +∞ 0 dω × X dκ ge(k) ω,dκ ˆC e ω,dκ + X dµ gm(k) ω,dµ ˆC m ω,dµ + h.c. =ℏ NX k=1 ˆσ(k) x ⊗ Z +∞ 0...
-
[5]
Double-continuum orthogonalization We next construct two sets of bright operators that are mutually orthonormal by taking suitable linear combina- tions of the ˆae(i) ω : ˆbe(j) ω = NX i=1 βe(ji) ω ˆae(i) ω forj∈[1, N e ind],(65a) ˆbe(j) ω = 0,forj∈[N e ind + 1, N],(65b) and for the ˆam(i) ω : ˆbm(j) ω = NX i=1 βm(ji) ω ˆam(i) ω ,forj∈[1, N m ind] (65c) ˆ...
-
[6]
The explicit expression is given in Appendix A
or by other orthonormalization procedures, such as the L¨ owdin orthonormalization [53, 66–68]. The explicit expression is given in Appendix A. By inversion of (65), one can write ˆae(i) ω = N e indX j=1 γe(ij) ω ˆbe(j) ω , i∈[1, N],(68a) ˆam(i) ω = N m indX j=1 γm(ij) ω ˆbm(j) ω , i∈[1, N].(68b) Theγ e(ij) ω andγ m(ij) ω coefficients are related to theβ ...
-
[7]
Direct construction of a single hybrid continuum Instead of the orthogonalization of the bright modes operators of the e- and m-continua separately, one can introduce a bright operator of ahybrid continuum, which includes both bright operators of e- and m-continua. We rewrite the interaction part of the Hamiltonian as follows ˆHint =− NX k=1 ˆσ(k) x ⊗ h d...
-
[8]
Orthogonalization of the hybrid bright operators We next construct a set of bright operators that are mutually orthonormal by taking suitable linear combina- tions of the ˆaω: ˆB(j) ω = NX i=1 β(ji) ω ˆa(i) ω , j∈[1, N ind], ˆB(j) ω = 0, j∈[N ind + 1, N],(90) whereN ind ≤Nis the number of linearly independent operators ˆa(i) ω (ri) and the coefficientsβ (...
-
[9]
or by other orthonormalization procedures, such as the L¨ owdin orthonormalization [53, 66–68], as described in Appendix A. By inversion of (90), one can write ˆa(i) ω = NindX j=1 γ(ij) ω ˆB(j) ω , i∈[1, N],(93) whereγ (ij) ω is defined as γ(ij) ω = NX l=1 β(jl)∗ ω M(il) ω (94) that, by (93), satisfies NX k=1 β(ik) ω γ(kj) ω =δ ij.(95) The interaction ter...
-
[10]
L. Novotny and B. Hecht,Principles of Nano-Optics (Cambridge University Press, 2012)
work page 2012
-
[11]
J. Pitarke, V. Silkin, E. Chulkov, and P. Echenique, The- ory of surface plasmons and surface-plasmon polaritons, Reports on progress in physics70
-
[12]
M. S. Tame, K. R. McEnery, S ¸. K.¨Ozdemir, J. Lee, S. A. Maier, and M. S. Kim, Quantum plasmonics, Nat. Phys. 9, 329 (2013). 14
work page 2013
- [13]
-
[14]
B. Rousseaux, D. Dzsotjan, G. Colas Des Francs, H. R. Jauslin, C. Couteau, and S. Gu´ erin, Adiabatic passage mediated by plasmons: A route towards a decoherence- free quantum plasmonic platform, Phys. Rev. B93, 045422 (2016); Erratum Phys. Rev. B94, 199902 (2016)
work page 2016
-
[15]
S. I. Bozhevolnyi and J. B. Khurgin, The case for quan- tum plasmonics, Nature Photonics11, 398 (2017)
work page 2017
- [16]
-
[17]
N. Westerberg and R. Bennett, Perturbative light–matter interactions; from first principles to inverse design, Physics Reports1026, 1
-
[18]
A. Gonz´ alez-Tudela, A. Reiserer, J. J. Garc´ ıa-Ripoll, and F. J. Garc´ ıa-Vidal, Light–matter interactions in quan- tum nanophotonic devices, Nature Reviews Physics6, 166 (2024)
work page 2024
-
[19]
A. F. Koenderink, Plasmon nanoparticle array waveg- uides for single photon and single plasmon sources, Nano Letters9, 4228 (2009)
work page 2009
-
[20]
H. Siampour, S. Kumar, and S. I. Bozhevolnyi, Nanofab- rication of plasmonic circuits containing single photon sources, ACS Photonics4, 1879 (2017)
work page 2017
-
[21]
M. Sugawara, Y. Xuan, Y. Mitsumori, K. Edamatsu, and M. Sadgrove, Plasmon-enhanced single photon source di- rectly coupled to an optical fiber, Physical Review Re- search4(2022)
work page 2022
-
[22]
D. E. Chang, A. S. Sørensen, P. R. Hemmer, and M. D. Lukin, Quantum optics with surface plasmons, Physical Review Letters97, 053002 (2006)
work page 2006
-
[23]
C. Lee, B. Lawrie, R. Pooser, K. G. Lee, C. Rockstuhl, and M. Tame, Quantum plasmonic sensors, Chemical Re- views121, 4743 (2021)
work page 2021
- [24]
-
[25]
J. Fregoni, F. J. Garcia-Vidal, and J. Feist, Theoreti- cal challenges in polaritonic chemistry, ACS Photonics 9, 1096 (2022)
work page 2022
-
[26]
B. Huttner, J. J. Baumberg, and S. M. Barnett, Canoni- cal quantization of light in a linear dielectric, Europhysics Letters16, 177 (1991)
work page 1991
-
[27]
B. Huttner and S. M. Barnett, Quantization of the elec- tromagnetic field in dielectrics, Physical Review A46, 4306 (1992)
work page 1992
-
[28]
B. Huttner and S. M. Barnett, Dispersion and loss in a Hopfield dielectric, Europhysics Letters18, 487 (1992)
work page 1992
-
[29]
S. M. Barnett, B. Huttner, and R. Loudon, Spontaneous emission in absorbing dielectric media, Physical Review Letters68, 3698 (1992)
work page 1992
-
[30]
J. J. Hopfield, Theory of the contribution of excitons to the complex dielectric constant of crystals, Physical Re- view112, 1555 (1958)
work page 1958
-
[31]
M. Wubs and L. G. Suttorp, Transient QED effects in ab- sorbing dielectrics, Physical Review A63, 043809 (2001)
work page 2001
-
[32]
L. G. Suttorp and M. Wubs, Field quantization in inho- mogeneous absorptive dielectrics, Physical Review A70, 013816 (2004)
work page 2004
-
[33]
L. G. Suttorp and A. J. Van Wonderen, Fano diagonal- ization of a polariton model for an inhomogeneous ab- sorptive dielectric, Europhysics Letters67, 766 (2004)
work page 2004
-
[34]
M. Wubs, L. G. Suttorp, and A. Lagendijk, Multiple- scattering approach to interatomic interactions and su- perradiance in inhomogeneous dielectrics, Physical Re- view A70, 053823 (2004)
work page 2004
-
[35]
N. A. Bhat and J. Sipe, Hamiltonian treatment of the electromagnetic field in dispersive and absorptive struc- tured media, Physical Review A73, 063808 (2006)
work page 2006
-
[36]
Friedrichs, On the perturbation of continuous spectra, Commun
K. Friedrichs, On the perturbation of continuous spectra, Commun. Pure Appl. Math.1, 361 (1948)
work page 1948
-
[37]
K. O. Friedrichs,Perturbation of spectra in Hilbert space, Vol. 3 (American Mathematical Soc., 1965)
work page 1965
-
[38]
U. Fano, Effects of configuration interaction on intensities and phase shifts, Physical Review124, 1866 (1961)
work page 1961
-
[39]
T. Gruner and D. G. Welsch, Correlation of radiation- field ground-state fluctuations in a dispersive and lossy dielectric, Physical Review A51, 3246 (1995)
work page 1995
-
[40]
T. Gruner and D. G. Welsch, Green-function approach to the radiation-field quantization for homogeneous and inhomogeneous Kramers-Kronig dielectrics, Physical Re- view A53, 1818 (1996)
work page 1996
-
[41]
H. T. Dung, L. Kn¨ oll, and D.-G. Welsch, Three- dimensional quantization of the electromagnetic field in dispersive and absorbing inhomogeneous dielectrics, Physical Review A57, 3931 (1998)
work page 1998
- [42]
-
[43]
S. Scheel and S. Buhmann, Macroscopic quantum elec- trodynamics — Concepts and applications, Acta Physica Slovaca58, 675 (2008)
work page 2008
- [44]
-
[45]
R. Matloob, R. Loudon, S. M. Barnett, and J. Jef- fers, Electromagnetic field quantization in absorbing di- electrics, Physical Review A52, 4823 (1995)
work page 1995
-
[46]
R. Matloob, R. Loudon, M. Artoni, S. M. Barnett, and J. Jeffers, Electromagnetic field quantization in amplify- ing dielectrics, Physical Review A55, 1623 (1997)
work page 1997
-
[47]
A. Drezet, Dual-Lagrangian description adapted to quan- tum optics in dispersive and dissipative dielectric media, Physical Review A94, 053826 (2016)
work page 2016
-
[48]
A. Drezet, Quantizing polaritons in inhomogeneous dis- sipative systems, Physical Review A95, 023831 (2017)
work page 2017
-
[49]
A. Drezet, Equivalence between the Hamiltonian and Langevin noise descriptions of plasmon polaritons in a dispersive and lossy inhomogeneous medium, Physical Review A96, 033849 (2017)
work page 2017
- [50]
- [51]
-
[52]
D. Y. Na, J. Zhu, and W. C. Chew, Diagonalization of the hamiltonian for finite-sized dispersive media: Canonical quantization with numerical mode decomposition, Phys- ical Review A103, 063707 (2021)
work page 2021
-
[53]
D. Y. Na, T. E. Roth, J. Zhu, W. C. Chew, and C. J. Ryu, Numerical framework for modeling quantum electromag- netic systems involving finite-sized lossy dielectric objects 15 in free space, Physical Review A107, 063702 (2023)
work page 2023
-
[54]
G. Semin, H.-R. Jauslin, and S. Gu´ erin, Canonical quan- tum plasmonics for finite nanostructures: Exact solution in one dimension, Physical Review A110, 043514 (2024)
work page 2024
-
[55]
G. Semin, H.-R. Jauslin, and S. Gu´ erin, Three- dimensional canonical quantum plasmonics for finite me- dia: Exact solution in terms of the classical green tensor, Physical Review A112(2025)
work page 2025
- [56]
-
[57]
A. Ciattoni, Quantum electrodynamics of lossy magne- todielectric samples in vacuum: Modified Langevin noise formalism, Physical Review A110, 013707 (2024)
work page 2024
- [58]
- [59]
- [60]
-
[61]
T. H¨ ummer, F. Garc´ ıa-Vidal, L. Mart´ ın-Moreno, and D. Zueco, Weak and strong coupling regimes in plasmonic QED, Physical Review B87, 115419 (2013)
work page 2013
-
[62]
A. Castellini, H. R. Jauslin, B. Rousseaux, D. Dzsotjan, G. Colas des Francs, A. Messina, and S. Gu´ erin, Quan- tum plasmonics with multi-emitters: application to stim- ulated Raman adiabatic passage, The European Physical Journal D72, 223 (2018)
work page 2018
- [63]
-
[64]
L. Kn¨ oll, S. Scheel, and D. Welsch,Coherence and statis- tics of photons and atoms(Wiley New York, 2001)
work page 2001
-
[65]
S. Y. Buhmann and D.-G. Welsch, Casimir-polder forces on excited atoms in the strong atom-field coupling regime, Physical Review A77, 012110 (2008)
work page 2008
-
[66]
D. Dzsotjan, A. S. Sørensen, and M. Fleischhauer, Quan- tum emitters coupled to surface plasmons of a nanowire: A Green’s function approach, Physical Review B82, 075427 (2010)
work page 2010
-
[67]
C. V. Vlack, P. T. Kristensen, and S. Hughes, Spon- taneous emission spectra and quantum light-matter in- teractions from a strongly coupled quantum dot metal- nanoparticle system, Physical Review B85, 075303 (2012)
work page 2012
- [68]
-
[69]
D. Dzsotjan, B. Rousseaux, H. R. Jauslin, G. C. des Francs, C. Couteau, and S. Guerin, Mode-selective quan- tization and multimodal effective models for spherically layered systems, Physical Review A94, 023818 (2016)
work page 2016
-
[70]
H. Varguet, B. Rousseaux, D. Dzsotjan, H.-R. Jauslin, S. Gu´ erin, and G. C. des Francs, Non-hermitian hamil- tonian description for quantum plasmonics: from dissi- pative dressed atom picture to fano states, Journal of Physics B: Atomic, Molecular and Optical Physics52, 055404 (2019)
work page 2019
-
[71]
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, Nanophotonics11, 4363 (2022)
work page 2022
- [72]
- [73]
-
[74]
S. Y. Buhmann,Dispersion Forces I, II(Springer, 2013)
work page 2013
-
[75]
P.-O. L¨ owdin, On the non-orthogonality problem con- nected with the use of atomic wave functions in the the- ory of molecules and crystals, The Journal of Chemical Physics18, 365 (1950)
work page 1950
-
[76]
L¨ owdin, Quantum theory of cohesive properties of solids, Advances in Physics5, 597 (1956)
P.-O. L¨ owdin, Quantum theory of cohesive properties of solids, Advances in Physics5, 597 (1956)
work page 1956
-
[77]
L¨ owdin,Linear Algebra for Quantum Theory(John Wiley & Sons, 1998)
P.-O. L¨ owdin,Linear Algebra for Quantum Theory(John Wiley & Sons, 1998)
work page 1998
-
[78]
C. Sauvan, J.-P. Hugonin, I. S. Maksymov, and P. Lalanne, Theory of the spontaneous optical emission of nanosize photonic and plasmon resonators, Physical Review Letters110, 237401 (2013)
work page 2013
- [79]
-
[80]
R. C. Ge and S. Hughes, Quantum dynamics of two quan- tum dots coupled through localized plasmons: An intu- itive and accurate quantum optics approach using quasi- normal modes, Physical Review B92, 205420 (2015)
work page 2015
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.