pith. sign in

arxiv: 2601.00709 · v2 · submitted 2026-01-02 · 🧮 math-ph · math.MP

On the computation of the dyadic Green's functions of Maxwell's equations in layered media

Heng Yuan , Wenzhong Zhang , Bo Wang This is my paper

Pith reviewed 2026-05-16 18:02 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords dyadic Green's functionslayered mediaMaxwell's equationsvector potentialTE/TM decompositioninterface conditionsfar-field approximationselastic wave equation
0
0 comments X

The pith

A vector potential formulation with matrix basis yields dyadic Green's functions equivalent to TE/TM decomposition in layered media.

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

The paper compares two ways to compute dyadic Green's functions for Maxwell's equations in layered media. One relies on the familiar TE/TM decomposition; the other uses a vector potential and a specially chosen matrix basis. The authors simplify the second derivation and prove the two approaches give identical results because the vector potential directly separates the interface conditions. This equivalence also lets the matrix basis isolate non-symmetric terms, which streamlines far-field approximations and extends the method to elastic waves.

Core claim

The vector potential formulation using a matrix basis is equivalent to the TE/TM decomposition for dyadic Green's functions in layered media, but its derivation is more straightforward because the interface conditions are directly decoupled using the vector potential. The matrix basis splits out all non-symmetric factors in the density functions, which facilitates the derivation of far-field approximations for the dyadic Green's functions.

What carries the argument

Matrix basis that splits non-symmetric factors in the density functions while the vector potential decouples interface conditions directly.

If this is right

  • The two formulations can be used interchangeably in any layered-medium computation.
  • Far-field asymptotic expansions become easier to obtain from the split density functions.
  • The same matrix-basis construction applies without change to dyadic Green's functions of the elastic wave equation.
  • Implementation of fast multipole methods for layered media gains a simpler boundary-handling step.

Where Pith is reading between the lines

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

  • Numerical codes for electromagnetic scattering in stratified media could adopt the vector-potential route to shorten the boundary-condition coding.
  • The decoupling technique may transfer to other vector wave problems whose boundary conditions currently require coupled scalar potentials.
  • Explicit verification on a three-layer configuration with oblique incidence would test whether the equivalence holds under all polarizations.

Load-bearing premise

The chosen vector potential and matrix basis fully capture every boundary condition and polarization effect without adding or dropping terms that would break equivalence.

What would settle it

Numerical evaluation of both formulations at the same observation point inside a two-layer dielectric stack that produces visibly different field values would disprove equivalence.

read the original abstract

In this paper, two formulations for the computation of the dyadic Green's functions of Maxwell's equations in layered media are presented in details. The first formulation derived using TE/TM decomposition is well-known and intensively used in engineering community while the second formulation derived using vector potential and a matrix basis is recently used in establishing a fast multipole method. We significantly simplify the derivation of second formulation and show that it is equivalent to the first one while the derivation is more straightforward as the interface conditions are directly decoupled using the vector potential. The matrix basis is designed to split out all non-symmetric factors in the density functions which facilitates the derivation of far-field approximations for the dyadic Green's functions. Moreover, it can be applied to the computation of the dyadic Green's functions of elastic wave equation in layered media.

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

2 major / 3 minor

Summary. The manuscript presents two formulations for the dyadic Green's functions of Maxwell's equations in layered media. The first uses the standard TE/TM decomposition; the second employs a vector potential together with a specially constructed matrix basis. The authors simplify the derivation of the second formulation, assert its equivalence to the first, and note that the matrix basis isolates non-symmetric factors to aid far-field approximations. They also indicate that the same matrix-basis approach extends to the elastic-wave case.

Significance. If the equivalence is established rigorously, the work supplies a cleaner route to the same Green's functions that are already used in fast-multipole implementations. The explicit decoupling of interface conditions and the factorization of non-symmetric terms are concrete technical improvements that could reduce algebraic overhead in both analytic and numerical work. The noted applicability to elastic waves broadens the potential utility beyond electromagnetics.

major comments (2)
  1. [§3.2] §3.2, after Eq. (18): the claim that the chosen vector potential and matrix basis automatically enforce continuity of all tangential E and H components (and the associated normal D, B jumps) for arbitrary layer counts is asserted but not verified explicitly. The provided algebra is carried out for a single interface and then stated to generalize; an inductive step or explicit N=3 calculation is required to confirm that no cross terms appear when multiple interfaces are present.
  2. [§4] §4, Eq. (32)–(35): the equivalence mapping between the TE/TM amplitudes and the vector-potential coefficients is shown only for the spectral-domain representation. The passage to the spatial domain (Sommerfeld integrals) is not re-checked for preservation of the same boundary conditions; any truncation or contour deformation that differs between the two formulations could break the claimed equivalence.
minor comments (3)
  1. [§5] The abstract states that the matrix basis 'facilitates the derivation of far-field approximations,' yet §5 only sketches the leading-order term; a short explicit far-field expression (or reference to where it appears) would strengthen the claim.
  2. [§2–3] Notation for the layer indices (superscript (l) versus subscript l) is used interchangeably in §2 and §3; a single consistent convention would improve readability.
  3. A brief comparison table listing the number of independent scalar functions retained in each formulation would help readers see the claimed simplification at a glance.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive suggestions. We address each major comment below and will revise the manuscript accordingly to strengthen the rigor of the derivations.

read point-by-point responses

  1. Referee: [§3.2] §3.2, after Eq. (18): the claim that the chosen vector potential and matrix basis automatically enforce continuity of all tangential E and H components (and the associated normal D, B jumps) for arbitrary layer counts is asserted but not verified explicitly. The provided algebra is carried out for a single interface and then stated to generalize; an inductive step or explicit N=3 calculation is required to confirm that no cross terms appear when multiple interfaces are present.

    Authors: We agree that an explicit verification for multiple interfaces is needed to confirm the absence of cross terms. In the revised manuscript we will insert an inductive argument immediately after Eq. (18) in §3.2. The induction proceeds by noting that the matrix basis decouples each interface independently; the continuity conditions at the k-th interface depend only on the local layer parameters and the already-satisfied conditions from the previous interface, with no additional coupling introduced by the vector-potential representation. revision: yes

  2. Referee: [§4] §4, Eq. (32)–(35): the equivalence mapping between the TE/TM amplitudes and the vector-potential coefficients is shown only for the spectral-domain representation. The passage to the spatial domain (Sommerfeld integrals) is not re-checked for preservation of the same boundary conditions; any truncation or contour deformation that differs between the two formulations could break the claimed equivalence.

    Authors: The equivalence is established pointwise in the spectral domain, and both formulations are converted to the spatial domain by exactly the same Sommerfeld integral. Consequently the boundary conditions, once satisfied spectrally, remain satisfied after the identical transformation. We will add a clarifying sentence in §4 stating that the mapping commutes with the integral operator and that no distinct contour or truncation is applied to either formulation. revision: yes

Circularity Check

0 steps flagged

Derivation chain is self-contained from Maxwell's equations with no reduction to inputs by construction

full rationale

The paper derives the dyadic Green's functions using two independent approaches: the established TE/TM decomposition and a vector-potential formulation with a designed matrix basis. Equivalence is established by showing that the vector-potential approach directly decouples interface conditions without additional correction terms, which follows from the standard continuity requirements on tangential E and H fields in Maxwell's equations. No parameter is fitted to data and then relabeled as a prediction, no self-citation supplies a uniqueness theorem that forces the result, and the matrix basis is explicitly constructed to factor non-symmetric terms rather than being smuggled in from prior work. The derivation remains independent of the target result and does not collapse to a renaming or self-definition.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on standard Maxwell's equations and vector calculus in piecewise-homogeneous media; no new entities or fitted parameters are introduced in the abstract.

axioms (1)
  • domain assumption Maxwell's equations govern electromagnetic fields in each homogeneous layer with appropriate continuity conditions at interfaces.
    Invoked implicitly as the starting point for both formulations.

pith-pipeline@v0.9.0 · 5433 in / 1162 out tokens · 25742 ms · 2026-05-16T18:02:37.044493+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

18 extracted references · 18 canonical work pages

  1. [1]

    Electromagnetic fields of buried sources in stratified anisotropic media.IEEE Transactions on Antennas and Propagation, 27(5):671–678, 1979

    S Ali and S Mahmoud. Electromagnetic fields of buried sources in stratified anisotropic media.IEEE Transactions on Antennas and Propagation, 27(5):671–678, 1979

    work page 1979

  2. [2]

    Arbel and L

    E. Arbel and L. B. Felsen. Theory of radiation from sources in anisotropic media, part i: General sources in stratified media.Electromagnetic Theory and Antennas, pages 391–420, 1963

    work page 1963

  3. [3]

    S Barkeshli. On the electromagnetic dyadic green’s functions for planar multi-layered anistropic uniaxial material media.International journal of infrared and millimeter waves, 13(4):507–527, 1992

    work page 1992

  4. [4]

    Cambridge University Press, 2013

    Wei Cai.Computational Methods for Electromagnetic Phenomena: electrostatics in solvation, scattering, and electron transport. Cambridge University Press, 2013

    work page 2013

  5. [5]

    A matrix-friendly formulation of layered medium green’s function.IEEE Antennas and Wireless Propagation Letters, 5:490–494, 2006

    WC Chew, JL Xiong, and MA Saville. A matrix-friendly formulation of layered medium green’s function.IEEE Antennas and Wireless Propagation Letters, 5:490–494, 2006

    work page 2006

  6. [6]

    John Wiley & Sons, 1999

    Weng Cho Chew.Waves and fields in inhomogenous media. John Wiley & Sons, 1999

    work page 1999

  7. [7]

    SIAM, 2013

    David Colton and Rainer Kress.Integral equation methods in scattering theory. SIAM, 2013

    work page 2013

  8. [8]

    Nonuniqueness of resolution of hertz vector in presence of a boundary, and the horizontal dipole problem.IEEE Transactions on Antennas and Propagation, 17(3):376–378, 1969

    A Erteza and B Park. Nonuniqueness of resolution of hertz vector in presence of a boundary, and the horizontal dipole problem.IEEE Transactions on Antennas and Propagation, 17(3):376–378, 1969

    work page 1969

  9. [9]

    Chapman and Hall/CRC, 2021

    Walton C Gibson.The method of moments in electromagnetics. Chapman and Hall/CRC, 2021

    work page 2021

  10. [10]

    John Wiley & Sons, 2015

    Jian-Ming Jin.Theory and computation of electromagnetic fields. John Wiley & Sons, 2015

    work page 2015

  11. [11]

    J. A. Kong. Electromagnetic fields due to dipole antennas over stratified anisotropic media.Geophysics, 37(6):985–996, 1972

    work page 1972

  12. [12]

    On the scalar potential of a point charge associated with a time- harmonic dipole in a layered medium.IEEE transactions on antennas and propagation, 35(11):1299–1301, 2003

    K Michalski. On the scalar potential of a point charge associated with a time- harmonic dipole in a layered medium.IEEE transactions on antennas and propagation, 35(11):1299–1301, 2003

    work page 2003

  13. [13]

    Multilayered media green’s functions in integral equation formulations.IEEE Transactions on Antennas and Propagation, 45(3):508–519, 1997

    Krzysztof A Michalski and Juan R Mosig. Multilayered media green’s functions in integral equation formulations.IEEE Transactions on Antennas and Propagation, 45(3):508–519, 1997

    work page 1997

  14. [14]

    Electromagnetic scattering and radiation by surfaces of arbitrary shape in layered media

    Krzysztof A Michalski and Dalian Zheng. Electromagnetic scattering and radiation by surfaces of arbitrary shape in layered media. i. theory.IEEE Transactions on Antennas and propagation, 38(3):335–344, 1990

    work page 1990

  15. [15]

    ¨Uber die ausbreitung der wellen in der drahtlosen telegraphie.Ann

    Arnold Sommerfeld. ¨Uber die ausbreitung der wellen in der drahtlosen telegraphie.Ann. Physik, 28:665–737, 1909

    work page 1909

  16. [16]

    Dyadic green function for the elec- tromagnetic field in multilayered isotropic media: An operator approach

    T Sphicopoulos, V Teodoridis, and FE Gardiol. Dyadic green function for the elec- tromagnetic field in multilayered isotropic media: An operator approach. InIEE Pro- ceedings H (Microwaves, Antennas and Propagation), volume 132, pages 329–334. IET, 1985

    work page 1985

  17. [17]

    A newly developed formulation suitable for matrix manipulation of layered medium green’s functions.IEEE transactions on antennas and propagation, 58(3):868–875, 2009

    Jie L Xiong and Weng Cho Chew. A newly developed formulation suitable for matrix manipulation of layered medium green’s functions.IEEE transactions on antennas and propagation, 58(3):868–875, 2009

    work page 2009

  18. [18]

    A matrix basis formulation for the dyadic green’s functions of maxwell’s equations in layered media.SIAM Journal on Applied Mathematics, 82(5):1710–1732, 2022

    Wenzhong Zhang, Bo Wang, and Wei Cai. A matrix basis formulation for the dyadic green’s functions of maxwell’s equations in layered media.SIAM Journal on Applied Mathematics, 82(5):1710–1732, 2022

    work page 2022