pith. sign in

arxiv: 2606.27445 · v1 · pith:7EAGVEWOnew · submitted 2026-06-25 · ⚛️ physics.optics · cs.NA· math.NA

Analysis of Nonlinear Random Polarization in Dispersive Dielectrics

Nathan L. Gibson , Emmanuel E. Oguadimma This is my paper

Pith reviewed 2026-06-29 01:19 UTC · model grok-4.3

classification ⚛️ physics.optics cs.NAmath.NA
keywords nonlinear Debye polarizationpolynomial chaos expansionYee discretizationdispersive dielectricsuncertainty quantificationfinite-difference time-domainelectromagnetic wave propagationsensitivity analysis
0
0 comments X

The pith

Polynomial chaos expansions turn random nonlinear Debye polarization into a deterministic coupled system solved by a second-order Yee scheme.

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

The paper models electromagnetic wave propagation in dispersive dielectrics whose nonlinear Debye polarization carries random perturbations. It applies polynomial chaos expansions to convert the stochastic equations into a larger but fully deterministic system of coupled equations. The standard Yee finite-difference scheme is extended to this system and shown to retain second-order spatial accuracy, with numerical tests confirming the expected convergence rates. Sensitivity studies then quantify how uncertainty in the nonlinear response grows with input signal amplitude, addressing the practical problem that manufacturing variations can shift optimal material parameters and spoil intended nonlinear behavior.

Core claim

The paper claims that random perturbations in a nonlinear Debye polarization model can be represented by a polynomial chaos expansion whose resulting deterministic system admits a second-order-accurate Yee discretization; numerical verification establishes convergence, and the same framework reveals that nonlinear properties become markedly more sensitive to parameter uncertainty at large input amplitudes, thereby supporting simulation-based identification of realizable materials whose desired effects survive fabrication variations.

What carries the argument

The polynomial chaos expansion of the random nonlinear Debye polarization, which produces an enlarged deterministic system to which the Yee finite-difference time-domain discretization extends while preserving second-order spatial accuracy.

If this is right

  • The extended Yee scheme solves the transformed deterministic system to second-order accuracy in space.
  • Numerical experiments confirm the expected convergence rates for the coupled system.
  • Nonlinear response properties exhibit greater sensitivity to uncertainty at large input signal amplitudes.
  • The framework supports model-based selection of material parameters that remain effective despite random manufacturing variations.

Where Pith is reading between the lines

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

  • The same expansion-plus-Yee approach could be tested on other dispersive nonlinear models such as Kerr or Raman media to check whether second-order accuracy persists.
  • Large-amplitude regimes may require adaptive choice of chaos expansion order to keep truncation error from contaminating sensitivity estimates.
  • The method supplies a concrete route to quantify how tight manufacturing tolerances must be to preserve a target nonlinear optical effect.
  • Extension to three-dimensional domains or to inverse design problems would follow directly from the deterministic coupled system already derived.

Load-bearing premise

The truncation error of the polynomial chaos expansion remains small enough that it does not degrade the second-order spatial accuracy of the extended Yee scheme.

What would settle it

A convergence study in which the observed order of accuracy falls below two when the polynomial chaos truncation order is deliberately lowered or when the input amplitude is increased to the regime where nonlinear sensitivity is reported to be highest.

Figures

Figures reproduced from arXiv: 2606.27445 by Emmanuel E. Oguadimma, Nathan L. Gibson.

Figure 2.1
Figure 2.1. Figure 2.1: Complex permittivity of the nonlinearly forced Debye model under relaxation-time uncertainty. [PITH_FULL_IMAGE:figures/full_fig_p006_2_1.png] view at source ↗
Figure 2.2
Figure 2.2. Figure 2.2: Complex permittivity of the linear Debye model and the nonlinearly forced Debye model under [PITH_FULL_IMAGE:figures/full_fig_p007_2_2.png] view at source ↗
Figure 2.3
Figure 2.3. Figure 2.3: Cole–Cole plots of the dielectric response for the nonlinearly forced Debye model under relaxation [PITH_FULL_IMAGE:figures/full_fig_p007_2_3.png] view at source ↗
Figure 5.1
Figure 5.1. Figure 5.1: Error over spatial and temporal domain with Uniform random variable by Legendre Chaos. [PITH_FULL_IMAGE:figures/full_fig_p013_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Snapshots of the Electric Field in the Material at [PITH_FULL_IMAGE:figures/full_fig_p015_5_2.png] view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: The electric field recorded at a distance of 0.016m and 0.02m from the antenna [PITH_FULL_IMAGE:figures/full_fig_p015_5_3.png] view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: The electric field recorded at a distance of 0.0016m and 0.0025m from the antenna [PITH_FULL_IMAGE:figures/full_fig_p016_5_4.png] view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: Power Spectrum of Signal at the antenna and at a distances of 0.015m and 0.016m from the [PITH_FULL_IMAGE:figures/full_fig_p016_5_5.png] view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: Power Spectrum of Signal at a distance of 0.0015m and 0.0025m from the antenna [PITH_FULL_IMAGE:figures/full_fig_p017_5_6.png] view at source ↗
Figure 6.1
Figure 6.1. Figure 6.1: Optimization results for deterministic (left) and random (right) models [PITH_FULL_IMAGE:figures/full_fig_p017_6_1.png] view at source ↗
Figure 6.2
Figure 6.2. Figure 6.2: Optimization results for linear (left) and nonlinear (right) models [PITH_FULL_IMAGE:figures/full_fig_p018_6_2.png] view at source ↗
Figure 6.3
Figure 6.3. Figure 6.3: Optimization results for linear deterministic (left) and nonlinear random (right) models [PITH_FULL_IMAGE:figures/full_fig_p018_6_3.png] view at source ↗
read the original abstract

We present a study on the time-domain propagation of electromagnetic waves in dielectric materials modeled by a nonlinear Debye medium with random perturbations. Polynomial Chaos Expansions are employed to transform the random nonlinear Debye polarization model into a deterministic framework. We extend the Yee discretization to the resulting coupled system, establish second order accuracy, and verify convergence numerically. We investigate the sensitivity of nonlinear properties to uncertainty, particularly when the amplitude of the input signal is large. Given the challenges in manufacturing where uncertainties can cause optimal parameters to vary and potentially disrupt nonlinear effects, our approach incorporates these uncertainties within the simulation. This can enable the model-based design identification of realizable materials that maintain their desired effects despite variations. The findings from this study contribute to a deeper understanding of wave propagation in complex media, with potential implications for applications in optical communications, material science, and electromagnetic wave control.

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 / 2 minor

Summary. The paper claims to study time-domain EM wave propagation in nonlinear Debye dielectrics subject to random perturbations by applying Polynomial Chaos Expansions to obtain a deterministic coupled system, extending the Yee discretization to this system while establishing and numerically verifying second-order accuracy, and examining the sensitivity of nonlinear properties to uncertainty (especially at large input amplitudes) for applications in robust material design.

Significance. If the claimed second-order accuracy of the extended Yee scheme is rigorously supported and not degraded by PCE truncation in the nonlinear setting, the work would supply a practical framework for uncertainty quantification in dispersive nonlinear optics, enabling model-based identification of realizable materials whose desired effects persist under manufacturing variations.

major comments (2)
  1. [Abstract] Abstract: the claim that second-order accuracy is established and verified numerically is not supported by any a-priori error analysis or explicit bounds relating the PCE truncation remainder to the spatial discretization error; in the nonlinear polarization equations the remainder enters as a forcing term whose magnitude is not shown to remain O(h²).
  2. [Numerical results / convergence verification] Convergence verification (implied in the numerical results section): without details on how the nonlinear Debye polarization couples to the random expansion, it is unclear whether the observed convergence rate reaches 2 or saturates at the level set by the chosen PCE truncation order.
minor comments (2)
  1. Specify the exact PCE truncation order and the distributions chosen for the random material parameters in all reported experiments.
  2. Clarify the precise form of the extended Yee update equations for the coupled deterministic system obtained after PCE transformation.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive comments. We address each major comment below and will revise the manuscript accordingly to improve clarity on the error analysis and numerical coupling details.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that second-order accuracy is established and verified numerically is not supported by any a-priori error analysis or explicit bounds relating the PCE truncation remainder to the spatial discretization error; in the nonlinear polarization equations the remainder enters as a forcing term whose magnitude is not shown to remain O(h²).

    Authors: We acknowledge that the manuscript does not contain a complete a-priori error analysis with explicit bounds on the PCE truncation remainder. The second-order accuracy follows from extending the standard Yee-scheme analysis to the deterministic augmented system; the additional polarization-coefficient equations are discretized with the same centered differences. For a fixed PCE truncation order the remainder appears as a smooth, h-independent forcing term whose contribution does not alter the O(h²) spatial rate. In the revision we will add a short paragraph in the methods section making this argument explicit and will rephrase the abstract to state that second-order accuracy is verified numerically for the spatial discretization once the PCE order is chosen sufficiently high. revision: yes

  2. Referee: [Numerical results / convergence verification] Convergence verification (implied in the numerical results section): without details on how the nonlinear Debye polarization couples to the random expansion, it is unclear whether the observed convergence rate reaches 2 or saturates at the level set by the chosen PCE truncation order.

    Authors: We agree that the coupling mechanism and its effect on observed rates require more explicit description. The nonlinear term (product of field and polarization) is projected onto the PCE basis, producing a closed deterministic system whose nonlinear interactions are handled by the triple-product coefficients of the orthogonal polynomials. In the reported experiments a PCE order of 4 was used after verifying that order 5 produced changes below the discretization error. The revision will add a dedicated subsection describing this projection and will include supplementary convergence tables for two different PCE orders, confirming that the measured rate reaches 2 when the truncation error lies below the spatial error. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation applies standard Polynomial Chaos Expansion to convert the random nonlinear Debye polarization into a deterministic coupled system, then extends the classical Yee finite-difference scheme to that system. Second-order accuracy is asserted via discretization analysis of the extended scheme and confirmed by numerical convergence tests; neither step reduces by construction to a fitted parameter, self-defined quantity, or self-citation chain internal to the paper. The central claims rest on external, independently verifiable numerical methods rather than on any renaming, ansatz smuggling, or uniqueness theorem imported from the authors' prior work.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the applicability of standard PCE truncation to a nonlinear polarization equation and on the preservation of Yee-scheme order under that transformation. No new physical entities are introduced.

free parameters (2)
  • PCE truncation order
    Chosen to balance accuracy and computational cost; not derived from first principles.
  • Random variable distributions for material parameters
    Assumed forms (e.g., Gaussian or uniform) that are fitted or chosen to represent manufacturing uncertainty.
axioms (2)
  • domain assumption Polynomial Chaos Expansion converges in L2 for the nonlinear Debye polarization operator
    Invoked when transforming the random model into a deterministic coupled system (abstract).
  • domain assumption Extended Yee scheme retains second-order accuracy on the PCE-augmented system
    Stated as established without proof details in the abstract.

pith-pipeline@v0.9.1-grok · 5681 in / 1368 out tokens · 28550 ms · 2026-06-29T01:19:33.374323+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

29 extracted references

  1. [1]

    Albanese, J

    R. Albanese, J. Penn, and R. Medina. Ultrashort pulse response in nonlinear dispersive media. In H. Bertoni, L. Carin, and L. Felsen, editors,Ultra-Wideband, Short-Pulse Electromagnetics, pages 259–264. Plenum Press, 1993

    1993

  2. [2]

    Electromagnetic inverse problems involving distributions of dielectric mechanisms and parameters.Quarterly of Applied Mathematics, 64(4):749–795, 2006

    H Banks and N Gibson. Electromagnetic inverse problems involving distributions of dielectric mechanisms and parameters.Quarterly of Applied Mathematics, 64(4):749–795, 2006

    2006

  3. [3]

    Banks and Gabriella A

    H.T. Banks and Gabriella A. Pinter. High frequency pulse propagation in nonlinear dielectric materials. Nonlinear analysis: real world applications, Elsevier, 2003

    2003

  4. [4]

    Ross Macdonald, editors.Impedance Spectroscopy: Theory, Experiment, and Applications

    Evgenij Barsoukov and J. Ross Macdonald, editors.Impedance Spectroscopy: Theory, Experiment, and Applications. Wiley-Interscience, Hoboken, NJ, 2 edition, 2005

    2005

  5. [5]

    Generalized polynomial chaos and dispersive dielectric media.REU Program at Oregon State University Proceedings, 2010

    ERIN Bela, ERIK Hortsch, and NL Gibson. Generalized polynomial chaos and dispersive dielectric media.REU Program at Oregon State University Proceedings, 2010

    2010

  6. [6]

    C. J. F. Böttcher, O. C. Van Belle, P. Bordewijk, A. Rip, and D. D. Yue. Theory of electric polarization. Journal of The Electrochemical Society, 121:211C, 1974

    1974

  7. [7]

    Boyd.Nonlinear Optics

    Robert W. Boyd.Nonlinear Optics. Academic Press, 3 edition, 2008

    2008

  8. [8]

    Cole and Robert H

    Kenneth S. Cole and Robert H. Cole. Dispersion and absorption in dielectrics i. alternating current characteristics.The Journal of Chemical Physics, 9(4):341–351, 1941

    1941

  9. [9]

    Cole and Robert H

    Kenneth S. Cole and Robert H. Cole. Dispersion and absorption in dielectrics ii. direct current characteristics.The Journal of Chemical Physics, 10(2):98–105, 1942. 19

    1942

  10. [10]

    The dielectric properties of biological tissues: III

    Christoph Gabriel, Stuart Gabriel, and Elaine Corthout. The dielectric properties of biological tissues: III. parametric models for the dielectric spectrum of tissues.Physics in Medicine & Biology, 41(11):2271– 2293, 1996

    1996

  11. [11]

    Ghanem and S

    R. Ghanem and S. Dham. Stochastic finite element analysis for multiphase flow in heterogeneous porous media.Porous Media, 32:239–262, 1998

    1998

  12. [12]

    Nathan L. Gibson. A polynomial chaos method for dispersive electromagnetics.Journal of Computational Physics, 270:105–125, 2014

    2014

  13. [13]

    Calculation of gauss quadrature rules.Mathematics of computation, 23(106):221–230, 1969

    Gene H Golub and John H Welsch. Calculation of gauss quadrature rules.Mathematics of computation, 23(106):221–230, 1969

    1969

  14. [14]

    Multiterm debye dispersion relations for permittivity of muscle.IEEE Transactions on Biomedical Engineering, (1):60–64, 2007

    William D Hurt. Multiterm debye dispersion relations for permittivity of muscle.IEEE Transactions on Biomedical Engineering, (1):60–64, 2007

    2007

  15. [15]

    John Wiley & Sons, New York, 3rd edition, 1998

    John David Jackson.Classical Electrodynamics. John Wiley & Sons, New York, 3rd edition, 1998

    1998

  16. [16]

    Springer, Berlin, Heidelberg, 2002

    Friedrich Kremer and Andreas Schönhals, editors.Broadband Dielectric Spectroscopy. Springer, Berlin, Heidelberg, 2002

    2002

  17. [17]

    L. D. Landau and E. M. Lifshitz.Electrodynamics of Continuous Media. Pergamon Press, Oxford, 2nd edition, 1984

    1984

  18. [18]

    A nonlinear wave equation associated with a nonlinear integral equation

    Nguyen Thanh Long. A nonlinear wave equation associated with a nonlinear integral equation. 2011

    2011

  19. [19]

    Generalized polynomial chaos and random oscillators

    Didier Lucor, C-H Su, and George Em Karniadakis. Generalized polynomial chaos and random oscillators. International Journal for Numerical Methods in Engineering, 60(3):571–596, 2004

    2004

  20. [20]

    Oguadimma

    Emmanuel E. Oguadimma. Analysis of nonlinear random polarization in dispersive dielectrics. Master’s thesis, Oregon State University, Corvallis, OR, 2025

    2025

  21. [21]

    Oughstun and G

    Kurt E. Oughstun and G. C. Sherman.Electromagnetic Pulse Propagation in Causal Dielectrics. Springer-Verlag, 1994

    1994

  22. [22]

    Pleshko and I

    P. Pleshko and I. Palocz. Experimental observation of sommerfeld and brillouin precursors in the microwave domain.Physical Review Letters, 22:1201–1204, 1969

    1969

  23. [23]

    SIAM, 2024

    Ralph C Smith.Uncertainty quantification: theory, implementation, and applications. SIAM, 2024

    2024

  24. [24]

    Hagness.Computational Electrodynamics: The Finite-Difference Time- Domain Method

    Allen Taflove and Susan C. Hagness.Computational Electrodynamics: The Finite-Difference Time- Domain Method. Artech House, 3 edition, 2005

    2005

  25. [25]

    The homogeneous chaos.American Journal of Mathematics, 60(4):897–936, 1938

    Norbert Wiener. The homogeneous chaos.American Journal of Mathematics, 60(4):897–936, 1938

    1938

  26. [26]

    W. H. Hunter Woodward. Broadband dielectric spectroscopy—a practical guide. InBroadband Dielectric Spectroscopy: A Modern Analytical Technique, volume 1375 ofACS Symposium Series, pages 3–59. American Chemical Society, Washington, DC, 2021

    2021

  27. [27]

    Princeton university press, 2010

    Dongbin Xiu.Numerical methods for stochastic computations: a spectral method approach. Princeton university press, 2010

    2010

  28. [28]

    The Wiener–Askey polynomial chaos for stochastic differential equations.SIAM journal on scientific computing, 24(2):619–644, 2002

    Dongbin Xiu and George Em Karniadakis. The Wiener–Askey polynomial chaos for stochastic differential equations.SIAM journal on scientific computing, 24(2):619–644, 2002

    2002

  29. [29]

    Mathematical analysis of high frequency pulse propagation in nonlinear dielectric materials.Nonlinear Analysis: Real World Applications, 11(5):3453–3462, 2010

    Younis Zaidan. Mathematical analysis of high frequency pulse propagation in nonlinear dielectric materials.Nonlinear Analysis: Real World Applications, 11(5):3453–3462, 2010. 20

    2010