{"id":"20d6d640-3953-4c39-80e8-02e30fa822e6","arxiv_id":"2504.13029","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The electric field operator for a finite dispersive medium is expressed exactly through the classical Green tensor, and the Purcell factor reduces to the standard bulk formula through an exact cancellation of two continuum contributions.","lead":"This paper derives a complete quantum description of light interacting with a finite piece of metal or dielectric of arbitrary shape, expressing the electric field operator in terms of the classical Green tensor. It resolves a long-standing puzzle by showing that the standard bulk-medium formula for the Purcell factor remains valid for finite media because two omitted contributions cancel exactly.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The real Fourier basis in Sec. II.B is not orthonormal: Eq. (26) misses a factor 1/2, which propagates into the Green-tensor identities and the derivation of the central Purcell cancellation.","rationale":"The reader's weakest assumption was the unitary equivalence of Omega^2 and Omega_0^2, an external theorem from prior work. My reading finds a more direct, internally checkable problem: the real continuum basis used throughout the paper is not delta-normalized as claimed. Because the paper explicitly defines the integration measure dkappa without any compensating factor, the orthonormality and completeness relations (26)-(27) are wrong by a factor of 2. This error enters the derivation of the e-coefficient integral equation, the solution in terms of the medium Green tensor, and the LDOS identity that is essential to the Purcell calculation. The central claim may still be true, and the final cancellation is elegant, but the proof as written contains an internal inconsistency that must be resolved before the result can be considered reliable. The verdict should therefore be CONDITIONAL rather than ACCEPT: the paper is acceptable if the normalization is corrected and all subsequent factors are rederived consistently. I disagree with the reader's choice of the unitary-equivalence theorem as the weakest point; the basis normalization is more concrete and more directly tied to the paper's own equations.","tokens_in":24937,"tokens_out":16385,"duration_ms":170441,"concrete_test":"Evaluate the normalization integral: with Phi_kappa from Eq. (24) and kappa = kappa', compute integral_{R^3} d^3r Phi_kappa(r) . Phi_kappa(r) and compare with delta(0); the result is (1/2) delta(0), disproving Eq. (26). Then recompute Eq. (134) at coincident points, x = y = r, and compare with the known free-space value Im G_0+(r,r,omega) = (omega/(6 pi c)) I; the spectral sum as written gives (omega/(12 pi c)) I, i.e., half the correct value. Finally, repeat the derivation of Secs. III.F and V with the normalized basis Phi -> sqrt(2) Phi, and verify numerically or analytically that the cancellation leading to Eq. (153) still holds with the standard classical Green tensor.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The most load-bearing concern is an internal normalization inconsistency, not the external unitary-equivalence theorem. The basis functions are defined in Eq. (24) as Phi_kappa(r) = (2 pi)^(-3/2) epsilon cos(k.r) or sin(k.r) with k in a half-space. For kappa = kappa', the left side of Eq. (26) evaluates to (2 pi)^(-3) times the integral of cos^2(k.r) over R^3, which equals (1/2) delta(0), not delta(0). Thus Eq. (26) is false and the completeness relation (27) is off by a factor of 2, since the sum over zeta = c,s yields (2 pi)^(-3) cos(k.(r-r')) and the half-space integral gives half the full-space delta. This is not a harmless convention: the spectral identity (84), which rewrites a sum over Phi_kappa as (omega^2/c^2) G_perp_0+ + delta_perp, is used directly to derive the integral equation (95) for the e coefficients and hence the central solution (112). The same factor affects Eq. (134) for Im G_0 and therefore the LDOS identity (135)/(136) used in the Purcell calculation. If the basis is corrected by inserting a sqrt(2), many intermediate formulas change by factors of 2, and while the final cancellation in Sec. V may survive, the derivations as written are not internally consistent. The paper should be accepted only after this normalization error is fixed and all factors are rechecked.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a canonical quantization scheme for the electromagnetic field coupled to a finite, dispersive, dissipative, inhomogeneous dielectric in three dimensions. It writes the classical Hamiltonian as two oscillator continua, diagonalizes the squared-frequency operator via Lippmann-Schwinger equations, and quantizes the resulting modes in a bosonic Fock space. The electric field operator is split into an 'electromagnetic' term E_e and a 'medium' term E_m, whose coefficients are shown to satisfy Fredholm integral equations. The authors prove that these coefficients can be expressed through the classical Green tensor of the dielectric satisfying the Sommerfeld radiation condition, and they use this representation to derive a Green-tensor LDOS identity. Applying Fermi's golden rule to a two-level emitter, they obtain a decay rate in which the E_e contribution cancels exactly against a corresponding term in the E_m contribution, leaving the standard bulk Purcell formula Γ(r_a,ω_a) = (2ω_a^2/ħ ε_0 c^2) d·Im[G_m^+(r_a,r_a,ω_a)]·d. The stated goal is to justify the bulk formula for finite media and to provide exact Green-tensor expressions for quantum plasmonics.","tokens_in":25221,"tokens_out":10460,"duration_ms":112941,"significance":"The intended result is significant: if the derivation is correct, it extends the authors' one-dimensional exact solution [39] to three-dimensional finite media of arbitrary shape, addresses the long-standing difficulty that bulk formulas cannot be applied to finite media, and gives explicit Green-tensor expressions for all field operators and eigenfunctions. The paper is largely self-contained algebraically, with the Green-tensor solution, the LDOS identity, and the Purcell compensation derived in the text and appendices. The external input is the unitary equivalence theorem from [38], which is a reasonable reliance on prior work. The claim that the Purcell formula emerges as an exact consequence of the microscopic model, rather than as a postulate, is a genuine strength, and the cancellation mechanism between the e and m contributions is an instructive and publishable insight. However, the internal normalization consistency of the real Fourier basis must be repaired before the derivation can be considered correct as written.","major_comments":[{"comment":"The real basis functions are not orthonormal as stated. With Φ_κ(r) = (2π)^(-3/2) ε cos(k·r) or ε sin(k·r) and k restricted to half-space, evaluating Eq. (26) for κ=κ' gives (2π)^(-3) ∫ cos^2(k·r) d^3r = (1/2)δ(0), not δ(0), and the completeness relation (27) has the same factor-of-two error. This is not a harmless convention: the spectral identity (84) and the imaginary-part identity (134) both rely on this completeness, and these identities feed directly into the integral equation for e_κ (95) and the LDOS identity (135)-(136) used in the Purcell cancellation (151)-(153). As written, the derivation is internally inconsistent. The correction is straightforward, namely to insert a factor sqrt(2) in Φ_κ, but all subsequent factors in Sections III-V must then be rechecked.","section":"§II.B, Eqs. (24), (26), (27)"},{"comment":"The proof of the e-coefficient solution is not valid as printed. The text says the identity (113) is obtained by replacing r by x in Eq. (112), but Eq. (112) is the statement being proved; the starting point should be the integral equation (111), not its claimed solution. In addition, Eq. (113) has β(x) on the left-hand side although the integration variable is z; the factor should presumably be β(z). The intended argument can be repaired by multiplying Eq. (111) by β(x)G(r',x), integrating over x, using Eq. (B3), and then cancelling the common term, but the printed proof is circular and needs to be rewritten.","section":"§III.F.1, Eqs. (112)-(114)"},{"comment":"The double-continuum structure and the entire Fock-space quantization depend on the unitary equivalence of Ω^2 and Ω_0^2, cited from [38] and [44, Ch. 5] but not re-derived. Since the central claim would fail if this equivalence does not hold for some admissible finite geometry, the authors should state explicitly the hypotheses of that theorem (conditions on the coupling α and on the medium volume V) and confirm that the present model satisfies them. This is not an accusation of error, but the manuscript would be substantially easier to assess if this load-bearing input were stated.","section":"§II.D, paragraph after Eq. (43)"}],"minor_comments":[{"comment":"The phrase 'can by written' should read 'can be written'.","section":"Abstract"},{"comment":"All displayed subspaces are labelled F_B^0; the subscript should follow the Fock-space level n, e.g., F_n, otherwise the definition of the Fock space is garbled.","section":"§II.G, Eqs. (58a)-(58d)"},{"comment":"The frequency variable is used inconsistently: the prefactor and the Green-tensor arguments are sometimes written as ω and sometimes as ω_a. These equations should consistently use the atomic frequency ω_a.","section":"§V, Eqs. (149), (152), (153)"},{"comment":"The notation mixes ∫ d^3k with the multi-index integral ∫ dκ; using dκ consistently throughout would improve readability.","section":"Eq. (1) and Sec. II"},{"comment":"A short derivation of Eq. (B6) from Eq. (A7) would be helpful, since Eq. (134) is a central input to the LDOS identity and the factor conventions matter.","section":"Appendix B, Eq. (B6)"}],"recommendation":"major_revision","confidential_remarks":"The positive assessment of the reader is largely justified: the physical claim is plausible, nontrivial, and valuable if the algebra can be made internally consistent. The main obstacle is the normalization error in the real Fourier basis, which is local and fixable. I recommend requesting a revision that corrects the basis normalization, rechecks every factor in Sections III-V, and rewrites the proof of Eq. (112) so that it starts from Eq. (111). The paper need not be rejected if these fixes are carried out."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take on Semin, Jauslin, and Guerin. The paper is the natural 3D continuation of the authors' 1D exact solution, and the core conceptual move is sound: quantize the finite medium with a double continuum, express the field coefficients as solutions of Fredholm equations, and show they are given by the classical Green tensor. The exact cancellation in the Purcell factor is the payoff, and it does explain why the old bulk formulas with a small background work for finite media. The derivation is long but the structure is convincing, and the engagement with the literature is honest.\n\nThat said, the stress-test note is correct and it matters. The real basis in Eq. (24) is not orthonormal as claimed in Eq. (26): the cos/sin integrals over the half-space give (1/2)δ(k−k'), not δ(k−k'), so the completeness relation (27) is off by a factor of 2. This is not a harmless convention, because the spectral identity (84) uses that completeness to produce the δ⊥ term, and the same factor enters (134) and the LDOS identity. The fix is simple — insert a sqrt(2) into the basis or rescale the measure — but then every factor in the intermediate steps has to be rechecked. The final Purcell cancellation might survive, but the paper as written is internally inconsistent.\n\nThere are also smaller presentation problems. The proof of (112) says it substitutes Eq. (112) while the algebra actually uses Eq. (79), and Eq. (113) has a variable typo. The unitary equivalence theorem from [38] is a real dependency; it is plausible and cited, but the paper should either restate the theorem or make clear the result is conditional on it.\n\nBottom line: the paper deserves a serious referee. I would send it out, but I would not accept it in its current form. The normalization error must be fixed and the factors rechecked. Once that is done, this will be a useful reference for quantum plasmonics and macroscopic QED. It is the kind of paper I would bring to a reading group and cite.","headline":"A serious 3D extension with a correct-looking central cancellation, but a real factor-of-2 normalization error in the Fourier basis that must be fixed before acceptance.","tokens_in":25768,"tokens_out":3231,"would_cite":true,"duration_ms":32626,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that for a finite medium of arbitrary shape, the quantized electric-field operator is exactly expressible through the classical Green tensor of the medium, and that the atomic decay rate reduces to the standard…","keywords":["quantum plasmonics","canonical quantization","Green tensor","double continuum","Purcell effect","local density of states","finite-size medium","Lippmann-Schwinger equations"],"falsifier":"Compute the e-coefficients from Eq. (112) and the medium Green tensor from the Fredholm equation (108) for a small metallic sphere or slab using a standard numerical solver, and evaluate both sides of the LDOS identity (135); any discrepancy beyond numerical error would falsify the central claim. An equivalent check is to compute $\\Gamma_e$ and $\\Gamma_m$ separately at one emitter position and verify that their sum equals the standard formula while $\\Gamma_e$ alone does not.","tokens_in":24719,"feed_emoji":"⚛️","tokens_out":9415,"duration_ms":90867,"temperature":0.7,"pith_summary":"This paper aims to put quantum plasmonics for finite, arbitrarily shaped media on the same footing as bulk quantum optics. It constructs a canonical quantization of a dissipative and dispersive medium coupled to the electromagnetic field and proves that the electric-field observable is exactly expressible through the classical Green tensor of the medium. The central result is the spontaneous-emission decay rate of an atom near the medium, $\\Gamma(r_a,\\omega_a)=\\frac{2\\omega_a^2}{\\hbar\\varepsilon_0 c^2}\\,\\mathbf{d}\\cdot\\operatorname{Im}[\\bar{\\bar G}_m^+(r_a,r_a,\\omega_a)]\\cdot\\mathbf{d}$, which is identical to the standard bulk formula. The reason is an exact compensation: the decay contribution from the free-field-like continuum $\\Gamma_e$ is cancelled term by term by a matching piece inside $\\Gamma_m$. If correct, this justifies why bulk formulas with a small added dissipative background give the right Purcell factors for finite nanostructures.","feed_headline":"Decay near finite media follows the standard Green-tensor formula","feed_subtitle":"Two cancelling continuum terms explain why bulk formulas work for finite plasmonic nanostructures.","key_machinery":"The load-bearing object is the squared-frequency operator $\\Omega^2 = \\Omega_0^2 + V$ acting on the direct sum of field and medium oscillator degrees of freedom, together with its double-continuum eigenfunctions $\\psi^e_\\kappa,\\psi^m_\\mu$ obtained from Lippmann-Schwinger equations. The decisive identities are the integral equations (79) and (96) for the electric-field coefficients, whose unique solutions are written as $\\omega_\\kappa\\Phi_\\kappa$ plus a convolution with the medium Green tensor, and $-\\tilde\\alpha(x,\\nu)\\frac{\\nu^2}{c^2}\\bar{\\bar G}_m^+(r,x,\\nu)\\cdot n_j$; the Green-tensor LDOS identity then converts mode sums into $\\operatorname{Im}\\bar{\\bar G}_m^+$ and produces the exact cancellation in the decay rate.","core_discovery":"On the paper's own terms, the central discovery is that the canonical quantization of a finite, dissipative, dispersive medium coupled to the electromagnetic field has a double-continuous spectrum (one continuum reducing to free photons, one to medium oscillators in the uncoupling limit), and that the electric-field operator splits as $\\hat E=\\hat E_e+\\hat E_m$. The coefficients $e_\\kappa$ and $m_\\mu$ multiplying the two families of creation-annihilation operators satisfy Fredholm integral equations, and the paper proves that their unique solutions are exact functionals of the classical medium Green tensor $\\bar{\\bar G}_m^+$ satisfying the Sommerfeld radiation condition. Using the resulting Green-tensor LDOS identity, the spontaneous-emission rate of a two-level emitter is shown to be $\\Gamma(r_a,\\omega_a)=\\frac{2\\omega_a^2}{\\hbar\\varepsilon_0 c^2}\\,\\mathbf{d}\\cdot\\operatorname{Im}[\\bar{\\bar G}_m^+(r_a,r_a,\\omega_a)]\\cdot\\mathbf{d}$, identical to the standard bulk formula, because the $e_\\kappa\\otimes e_\\kappa^*$ contribution in $\\Gamma_m$ exactly cancels the whole of $\\Gamma_e$.","pith_inferences":["Because the cancellation is exact for arbitrary shapes, the standard bulk Purcell formula should also hold in geometries where the added-background prescription is least plausible, such as media with strong magnetic or anisotropic response; that is a testable extension of the paper's reasoning.","The free-field-like term $\\hat E_e$, although it cancels in the decay rate, could still leave observable signatures in other processes such as atom-surface dispersion forces or emitter-emitter interactions, where no analogous cancellation is guaranteed.","A practical recipe follows from the structure of the proof: solve the Fredholm equation for $\\bar{\\bar G}_m^+$ once numerically, then build every quantum observable from that tensor and the free modes, bypassing the direct diagonalization of the coupled Hamiltonian."],"forward_implications":["For a finite medium of arbitrary shape, the quantized electric-field operator is obtained exactly once the classical Green tensor of the medium is known, since both field coefficients are explicit functionals of that tensor.","The spontaneous-emission rate of an emitter at any position is $\\frac{2\\omega_a^2}{\\hbar\\varepsilon_0 c^2}\\,\\mathbf{d}\\cdot\\operatorname{Im}[\\bar{\\bar G}_m^+(r_a,r_a,\\omega_a)]\\cdot\\mathbf{d}$, identical to the formula used in bulk treatments.","In the uncoupling limit the model reproduces the free field: the medium part of the field operator vanishes and the medium Green tensor becomes the free Green tensor.","The exact term-by-term cancellation explains why the procedure of adding a small dissipative background extending to infinity, then letting the background tend to one, gives the right Purcell factor for finite media.","The same Green-tensor formulation carries over to anisotropic and to magnetic media without changing the structure of the derivation."],"supporting_citations":[{"why":"Supplies the canonical quantization of a finite medium and the proof that the coupled and uncoupled squared-frequency operators are unitarily equivalent, giving the double-continuum spectrum.","marker":"[38]"},{"why":"Provides the one-dimensional exact solution in terms of the classical Green function that this paper extends to three dimensions and arbitrary shapes.","marker":"[39]"},{"why":"Provides the standard bulk Purcell formula and the Fermi Golden Rule setup used as the comparison target for the decay rate.","marker":"[2]"},{"why":"States the conceptual difficulty of applying bulk formulas to finite media, which the exact compensation result resolves.","marker":"[34]"},{"why":"An example of the dissipative-background prescription whose predictions the finite-medium result is shown to reproduce.","marker":"[29]"},{"why":"Reports numerical simulations for finite lossy dielectrics whose agreement with the bulk decay formula is explained by the exact cancellation.","marker":"[35]"},{"why":"The spectral theorem cited for the unitary equivalence underlying the diagonalization and the Fock-space quantization.","marker":"[44]"},{"why":"Used to assert uniqueness of the Fredholm integral equation solution, so the Green tensor is genuinely the solution for the field coefficients.","marker":"[50]"}],"fun_headline_variants":["Finite media: double continuum cancels to standard LDOS","Exact Green-tensor Purcell formula from two cancelling continua","Finite media quantized: Purcell factor matches bulk Green tensor","Quantum plasmonics: two continua cancel to give standard decay","Exact solution: finite media Purcell factor from classical Green tensor"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole diagonalization rests on the theorem that $\\Omega^2=\\Omega_0^2+V$ is unitarily equivalent to the uncoupled $\\Omega_0^2$, a spectral result cited from earlier work and not re-derived here; if that equivalence fails for some finite geometry, the exact field-operator expressions and the decay-rate formula do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Finite media: double continuum cancels to standard LDOS","Exact Green-tensor Purcell formula from two cancelling continua","Finite media quantized: Purcell factor matches bulk Green tensor","Quantum plasmonics: two continua cancel to give standard decay","Exact solution: finite media Purcell factor from classical Green tensor"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000702,"raw_usage":{"total_tokens":3247,"prompt_tokens":1106,"completion_tokens":2141,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":722,"completion_tokens_details":{"reasoning_tokens":2053}},"tokens_in":722,"tokens_out":2141,"duration_ms":17104,"temperature":1.0,"reasoning_tokens":2053,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T12:18:13.069077+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the e-coefficients from Eq. (112) and the medium Green tensor from the Fredholm equation (108) for a small metallic sphere or slab using a standard numerical solver, and evaluate both sides of the LDOS identity (135); any discrepancy beyond numerical error would falsify the central claim. An equivalent check is to compute $\\Gamma_e$ and $\\Gamma_m$ separately at one emitter position and verify that their sum equals the standard formula while $\\Gamma_e$ alone does not.","supporting_citations":[{"cited_title":"Dorier, J","cited_arxiv_id":null,"evidence_quote":"Supplies the canonical quantization of a finite medium and the proof that the coupled and uncoupled squared-frequency operators are unitarily equivalent, giving the double-continuum spectrum."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the one-dimensional exact solution in terms of the classical Green function that this paper extends to three dimensions and arbitrary shapes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard bulk Purcell formula and the Fermi Golden Rule setup used as the comparison target for the decay rate."},{"cited_title":"Di Stefano, S","cited_arxiv_id":null,"evidence_quote":"Reports numerical simulations for finite lossy dielectrics whose agreement with the bulk decay formula is explained by the exact cancellation."},{"cited_title":"Semin, H.-R","cited_arxiv_id":null,"evidence_quote":"The spectral theorem cited for the unitary equivalence underlying the diagonalization and the Fock-space quantization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Used to assert uniqueness of the Fredholm integral equation solution, so the Green tensor is genuinely the solution for the field coefficients."}],"review_version":1}