pith. machine review for the scientific record. sign in

arxiv: 2605.12346 · v1 · submitted 2026-05-12 · ⚛️ physics.optics · math-ph· math.MP

Recognition: 2 theorem links

· Lean Theorem

General and concise operator approach to the dyadic Green's function of layered media

Aliaksandr Arlouski, Andrey Novitsky, Dongliang Gao, Lei Gao

Pith reviewed 2026-05-13 03:10 UTC · model grok-4.3

classification ⚛️ physics.optics math-phmath.MP
keywords dyadic Green's functionlayered mediaanisotropic mediaevolution operatorssurface impedance tensorscomputational photonicsnanophotonics
0
0 comments X

The pith

An operator method derives the dyadic Green's function for anisotropic layered media using evolution operators and surface impedance tensors.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an operator-based derivation for the dyadic Green's function in generic anisotropic planarly-layered media, covering both electric and magnetic fields. The function is built directly from the evolution operators of the individual layers together with surface impedance tensors. This construction automatically isolates the singular term that describes local field response at the source. A reader would care because the method replaces lengthy component-wise algebra with a compact operator formalism that simplifies both analysis and use in photonics calculations.

Core claim

The central claim is that the dyadic Green's function of a generic anisotropic planarly-layered medium for electric and magnetic fields is expressed through the evolution operators of the comprising layers and the surface impedance tensors, with the singular term naturally separated from the remaining terms.

What carries the argument

Evolution operators of the layers combined with surface impedance tensors, which together construct the full Green's function and isolate the singular term without separate regularization steps.

If this is right

  • The Green's function applies equally to electric and magnetic fields in anisotropic planar layers.
  • The singular term separates automatically, avoiding separate handling in applications.
  • The same operator structure extends to spherical and cylindrical layered geometries.
  • Conceptual understanding and practical computations in nanophotonics become simpler.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The operator form may integrate more readily into numerical solvers for multilayer devices than traditional component expansions.
  • Generalization to curved layers suggests the method could apply to modeling fields around spherical nanoparticles or cylindrical waveguides.

Load-bearing premise

Evolution operators of the layers can be combined with surface impedance tensors to yield the complete Green's function while automatically separating the singular term without extra ad-hoc regularization or component-wise algebra.

What would settle it

A direct comparison for a single isotropic layer in which the operator-derived Green's function fails to reproduce the known closed-form expression or violates the required boundary conditions at the interfaces.

Figures

Figures reproduced from arXiv: 2605.12346 by Aliaksandr Arlouski, Andrey Novitsky, Dongliang Gao, Lei Gao.

Figure 1
Figure 1. Figure 1: FIG. 1. The layered medium under study. There are [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. A typical form of the integrand in the expression [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. LDOS enhancement compared to its vacuum value [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: εg εg * g = 1 -2 -1 0 1 2 3 4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 z / λ Tr Im  GE / Tr Im  GE0 (a) εg εg * g = 4 -2 -1 0 1 2 3 4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 z / λ Tr Im  GE / Tr Im  GE0 (b) FIG. 5. Tr Im GˆE as a function of distance for the bi-layered gyrotropic medium, the permittivity tensors of the two lay￾ers being complex-conjugate of one another, with parameters ε0 = 4 and (a) g = 1; (b) g = 4. IV. SUMM… view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The surface impedance tensor (C7) is not defined at [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
read the original abstract

Dyadic Green's function is an important tool of computational photonics, giving deeper insights into light-matter interaction. We present an operator approach to the derivation of the dyadic Green's function of a generic anisotropic planarly-layered medium for both electric and magnetic fields. The resulting Green's function is expressed through the evolution operators (a kind of transfer matrices) of the comprising layers and the surface impedance tensors, the singular term being naturally separated from other terms. The operator approach to the Green's function simplifies both the conceptual understanding of the problem and the subsequent practical applications, some of which are demonstrated here. The proposed approach can be easily generalized to the case of spherical and cylindrical layers. The obtained results can be applied in nanophotonics engineering problems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper presents an operator approach to derive the dyadic Green's function (for both electric and magnetic fields) of a generic anisotropic planarly-layered medium. The Green's function is expressed via the evolution operators of the individual layers combined with surface impedance tensors, with the singular (delta-function) term emerging naturally from the algebra without additional regularization. The method is claimed to simplify both conceptual understanding and practical computations, with suggested extensions to spherical and cylindrical geometries.

Significance. If the derivation is valid, the work supplies a compact, unified framework for Green's functions in layered media that avoids component-wise algebra and ad-hoc steps. The explicit operator construction, impedance matching at interfaces, and automatic separation of the singular term are clear strengths that could streamline applications in nanophotonics and computational electromagnetics. The internal consistency for the planar case and the direct path to curved-geometry generalizations add value for engineering problems.

minor comments (2)
  1. The notation for evolution operators, impedance tensors, and field components should be introduced with a dedicated definitions subsection or table to improve accessibility for readers new to transfer-matrix methods.
  2. While the abstract and introduction mention demonstrated applications, the manuscript would benefit from a brief, self-contained example (with explicit numerical comparison to a known case) to illustrate the practical simplification.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our work on the operator approach to the dyadic Green's function in anisotropic layered media. The recognition of the method's conceptual clarity, natural separation of the singular term, and potential for generalization to curved geometries is appreciated. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

Derivation is self-contained; no circularity detected

full rationale

The manuscript presents a direct operator-based derivation of the dyadic Green's function for anisotropic layered media, starting from Maxwell's equations and composing layer evolution operators with interface impedance tensors. The singular delta-function term is stated to emerge automatically from the operator algebra rather than being inserted by hand or fitted. No step reduces the target Green's function to a parameter fit, a self-citation chain, or a redefinition of the input quantities. The approach re-uses standard transfer-matrix concepts but does not presuppose the final Green's function expression; the construction is therefore independent of the claimed result and remains falsifiable against known planar-media benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard Maxwell equations and the existence of well-defined evolution operators for each layer; no free parameters, invented entities, or ad-hoc axioms are mentioned in the abstract.

axioms (2)
  • standard math Maxwell's equations hold inside each anisotropic layer and at interfaces.
    Implicit background for any electromagnetic Green's function derivation.
  • domain assumption Evolution operators (transfer matrices) exist and can be multiplied across layers.
    Core building block invoked to express the Green's function.

pith-pipeline@v0.9.0 · 5430 in / 1413 out tokens · 105579 ms · 2026-05-13T03:10:04.110407+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    restore the EM field vectors in the physical space us- ing an analog of (A11) and compare the result with the definition of the DGF (3). A suitable orthogonal transformation for cylindrically symmetric layers would be E(r) = ∞X m=−∞ Z ˜Em(ξ, bz)ei(k0bzz+mφ)d(k0bz),(D1) and for spherically symmetric layers it would be E(r) = ∞X l=0 lX m=−l Flm(θ, φ)˜Elm(ξ)...

    work page

  2. [2]

    as functions of polarθand azimuthalφangles (same for ˜Hand ˜j). After substituting these into Maxwell’s equations, one obtains the evolution equations of the same form as (8) (ξnow playing the role of thezco- ordinate used for planar media), with ˜Ucyl = 4π c   1 εq   q×ε(q⊗q) −Pt b× q + m k0ξ e× φ (q⊗q)   − q× 0   ˜j,(D3) ˜Usph = 4π c   1 εq  ...

    work page

  3. [3]

    A. L. Kuzemsky,Statistical Mechanics and The Physics of Many-Particle Model Systems(World Scientific Pub- lishing Company, 2017)

    work page 2017

  4. [4]

    Rother,Green’s Functions in Classical Physics (Springer International Publishing, 2017)

    T. Rother,Green’s Functions in Classical Physics (Springer International Publishing, 2017)

    work page 2017

  5. [5]

    Tai,Dyadic Green Functions in Electromagnetic Theory(IEEE, 1994)

    C.-T. Tai,Dyadic Green Functions in Electromagnetic Theory(IEEE, 1994)

    work page 1994

  6. [6]

    Faryad and A

    M. Faryad and A. Lakhtakia,Infinite-Space Dyadic Green Functions in Electromagnetism(Morgan & Clay- pool Publishers, 2018)

    work page 2018

  7. [7]

    W. L. Barnes et al., J. Opt.22, 1 (2020)

    work page 2020

  8. [8]

    Luomei, T

    Y. Luomei, T. Yin, K. Tan, , and X. Chen, Inverse Prob- lems40, 045002 (2024)

    work page 2024

  9. [9]

    Bordo, J

    V. Bordo, J. Opt. Soc. Am. B29, 1799 (2012)

    work page 2012

  10. [10]

    K. Jain, L. Ruks, F. le Kien, and T. Busch, Phys. Rev. Research6, 033311 (2024)

    work page 2024

  11. [11]

    J. Ren, S. Franke, B. Van Drunen, and S. Hughes, Phys. Rev. A109, 013513 (2024)

    work page 2024

  12. [12]

    C. Qi, S. Zeng, N. J. A. Egarguin, D. Onofrei, and J. Chen, Journal of Electromagnetic Waves and Applica- tions38(1), 66 (2023)

    work page 2023

  13. [13]

    Novotny and B

    L. Novotny and B. Hecht,Principles of Nano-Optics, 2nd ed. (Cambridge University Press, 2012)

    work page 2012

  14. [14]

    W. C. Chew,Waves and Fields in Inhomogeneous Media (Wiley-IEEE Press, 1999)

    work page 1999

  15. [15]

    M. H. Cho and W. Cai, J. Sci. Comp.71, 1319–1350 (2016)

    work page 2016

  16. [16]

    Paulus, P

    M. Paulus, P. Gay-Balmaz, and O. J. F. Martin, Phys. Rev. E62, 5797 (2000)

    work page 2000

  17. [17]

    J. E. Sipe, J. Opt. Soc. Am. B4, 481 (1987)

    work page 1987

  18. [18]

    Hohenester,Nano and Quantum Optics(Springer Cham, 2021)

    U. Hohenester,Nano and Quantum Optics(Springer Cham, 2021)

    work page 2021

  19. [19]

    Koh, IEEE Transactions on Antennas and Propaga- tion72, 4410 (2024)

    I.-S. Koh, IEEE Transactions on Antennas and Propaga- tion72, 4410 (2024)

    work page 2024

  20. [20]

    Ground Penetrating Radar: From Theoretical Endeavors to Com- putational Electromagnetics, Signal Processing, Antenna Design and Field Applications

    S. Lambot, K. Wu, A. Sluyters, and J. Vanderdonckt, The Full-Wave Radar Equation for Wave Propagation 12 in Multilayered Media and Its Applications. In “Ground Penetrating Radar: From Theoretical Endeavors to Com- putational Electromagnetics, Signal Processing, Antenna Design and Field Applications” (eds. M. Serhir and D. Lesselier)(Wiley, 2024)

    work page 2024

  21. [21]

    S. M. Ali, T. M. Habashy, and J. A. Kong, J. Opt. Soc. Am. A9, 413 (1992)

    work page 1992

  22. [22]

    Okhmatovski and S

    V. Okhmatovski and S. Zheng,Theory and Computation of Electromagnetic Fields in Layered Media(Wiley-IEEE Press, 2024)

    work page 2024

  23. [23]

    Koh, IET Microwaves, Antennas & Propagation16, 405 (2022)

    I.-S. Koh, IET Microwaves, Antennas & Propagation16, 405 (2022)

    work page 2022

  24. [24]

    A. S. Potemkin, A. N. Poddubny, P. A. Belov, and Y. S. Kivshar, Phys. Rev. A86, 023848 (2012)

    work page 2012

  25. [25]

    F. I. Fedorov,Optics of Anisotropic Media(Akademia Nauk BSSR, 1958)

    work page 1958

  26. [26]

    F. I. Fedorov,Theory of Gyrotropy(Nauka i Tehnika, 1976)

    work page 1976

  27. [27]

    L. M. Barkovsky and A. N. Furs,Operator methods for describing optical fields in complex environments(Be- laruskaya Navuka, 2003)

    work page 2003

  28. [28]

    tangential

    Here and below, “tangential” means “orthogonal toz”

    work page

  29. [29]

    Hassani,Mathematical Physics: A Modern Introduc- tion to Its Foundations, 2nd ed

    S. Hassani,Mathematical Physics: A Modern Introduc- tion to Its Foundations, 2nd ed. (Springer, 2013)

    work page 2013

  30. [30]

    Stone and P

    M. Stone and P. Goldbart,Mathematics for Physics: A Guided Tour for Graduate Students, 1st ed. (Cambridge University Press, 2009)

    work page 2009

  31. [31]

    (8), which has the formg(z, z ′) =θ(z−z ′) Ω z z′, θbeing the Heaviside step function: R z z0 Ωz z′ ˜U(z′)dz′ =R g(z, z ′) ˜U(z′)dz′

    The inhomogeneous term can be obtained either by the variation of parameter, or by using the Green’s function of Eq. (8), which has the formg(z, z ′) =θ(z−z ′) Ω z z′, θbeing the Heaviside step function: R z z0 Ωz z′ ˜U(z′)dz′ =R g(z, z ′) ˜U(z′)dz′. The usual definition of the Green’s functiong(z, z ′) isLg(z, , z ′) =δ(z−z ′), whereLis the differential ...

    work page

  32. [32]

    We could have expressed ˜W(N+1) zf instead of ˜W(0) zin and use it in the expressions to follow, which would result in a slightly different form of the DGF, but it would be just as good for practical purposes

    work page

  33. [33]

    Shaposhnikov, M

    L. Shaposhnikov, M. Mazanov, D. A. Bobylev, F. Wilczek, and M. A. Gorlach, Phys. Rev. B108, 115101 (2023)

    work page 2023

  34. [34]

    G. N. Borzdov, J. Math. Phys.38, 6328 (1997)

    work page 1997

  35. [35]

    Born and E

    M. Born and E. Wolf,Principles of Optics, 4th ed. (Perg- amon Press, 1970)

    work page 1970

  36. [36]

    F. J. Dyson, Phys. Rev.75, 486 (1949)

    work page 1949

  37. [37]

    N. P. Armitage, E. J. Mele, and A. Vishwanath, Rev. Mod. Phys.90, 015001 (2018)

    work page 2018

  38. [38]

    C. Guo, V. S. Asadchy, B. Zhao, and S. Fan, eLight3, 2 (2023)

    work page 2023

  39. [39]

    Gantmacher,The Theory of Matrices, 2nd ed

    F. Gantmacher,The Theory of Matrices, 2nd ed. (Nauka, Moscow, 1966)

    work page 1966

  40. [40]

    Novitsky and L

    A. Novitsky and L. Barkovsky, Phys. Rev. A77, 033849 (2008)

    work page 2008

  41. [41]

    A. V. Novitsky and L. M. Barkovsky, J. Phys. A38, 391 (2005)

    work page 2005