Log-Periodic Spectral Hierarchies in a Boundary-Driven Electromagnetic Cavity: Evidence from FDTD Simulations

PaperLSH

byJill F. RankinPublished 5/13/2026AI Rating: 3.7/5

FDTD simulations show that explicit log-periodic time modulation of boundary permittivity (scaling ratio b = 13/8) in a one-dimensional electromagnetic cavity imprints a reproducible discrete-scale hierarchy in the probe-point power spectrum: after subtracting the smooth power-law envelope the residual oscillates periodically in ln ω with period ln(13/8). The effect is robust across parameter sweeps, absent in undriven/harmonic/random controls, statistically significant (r ≈ 0.81, p < 0.001), and persists with coherent propagation into the cavity interior.

Top 10% Internal Consistency

Top 10% Clarity

View Shareable Review Profile- permanent credential link for endorsements

Approved for Publication

Internal Consistency4/5

high confidence- spread 1- panel

The narrative is largely coherent: a boundary-layer ε(x,t) modulation periodic in ln t is applied, spectra are computed, a smooth envelope is removed, and residuals are fit to a fixed-period cosine in ln ω matching ln b. Controls (undriven/harmonic/random) are analyzed identically. Minor internal issues remain: (i) terminology tension between calling the domain a ‘cavity’ while using Mur absorbing boundaries (suggesting an open/absorbing domain rather than a reflective cavity), which affects interpretation of ‘modes’ and ‘hierarchy’; (ii) the residual is described as ‘after subtraction’ but then modeled as a relative residual R(ω) without a precise definition, leaving some ambiguity about what exactly is fit. These issues are patchable and do not on their own create a direct contradiction in the stated pipeline or reported statistics.

Mathematical Validity3/5

high confidence- spread 1- panel

Most explicit mathematics is syntactically and dimensionally plausible in the paper’s dimensionless units (Courant factor 0.5 is within the 1D stability limit; the ln((t+t0)/t0) argument is dimensionless). However, the main mathematical vulnerability is the load-bearing, unproven mapping from a drive periodic in ln t to a residual periodic in ln ω with the same ln-period (eqs. (1)–(2) are asserted as ‘predicted’/‘predicted form’ without derivation). Because the core conclusion is that the residual period ‘matches’ ln(13/8) consistent with an externally imposed discrete-scale symmetry, the absence of a transfer-function or even a reduced-model argument (e.g., parametric coupling/Floquet-in-log-time analysis) limits mathematical warrant. Additionally, the detrending step S_env(ω)∝ω^α and subsequent sinusoidal-in-lnω modeling are analysis assumptions that could generate apparent oscillations if the envelope model is slightly wrong. These gaps do not prove the result false, but they prevent the paper from reaching high mathematical validity as a claimed mechanism rather than an observed numerical correlation.

Falsifiability4/5

high confidence- spread 0- panel

The work makes a clear, differentiating prediction: imposing a boundary permittivity modulation periodic in ln(t) with scaling ratio b should produce, after envelope subtraction, a residual in the power spectrum that is periodic in ln(ω) with period ln(b), and this structure should be absent in undriven, harmonic, and random controls. That is a concrete observable with explicit fitting form, reported amplitudes, correlation metrics, and a stated null-hypothesis rejection criterion. The paper also gives several internal falsification handles: mismatch of extracted period with ln(13/8), comparable residuals in controls, or disappearance of the effect away from the boundary would undermine the proposed mechanism. This is strong by theory/simulation standards. The main reason this does not earn a 5 is that the falsification criteria are still somewhat analysis-dependent: the signal is defined after envelope subtraction and fixed-period fitting, and the manuscript does not fully explore alternative spectral estimators, different fitting windows, or possible numerical artifacts from finite duration, absorbing boundaries, and nonstationary driving. In addition, no experimental protocol with quantified uncertainty budget is developed, so near-term laboratory falsification is suggested rather than rigorously specified.

Clarity4/5

high confidence- spread 0- panel

Overall the paper is readable and organized in a conventional way: introduction, methods, post-processing, robustness checks, results, discussion, and conclusion. The central observable is defined clearly, the fitting model is explicit, and the controls are easy to understand. A graduate-level reader in computational electromagnetism or wave physics should be able to follow the main claim without difficulty. The main clarity weaknesses are localized but meaningful. First, there is an apparent text corruption in the introduction ('under [REDACTED].he boundary-conditioned modulation'), which interrupts the narrative and should be fixed. Second, some methodological details essential for evaluating robustness are only asserted briefly: for example, the exact source pulse parameters, how the boundary update is implemented in the Yee scheme, and how the randomized-phase surrogate test is applied to a nonstationary driven signal. Third, figures and supplemental material are referenced for key support, but several decisive details are not present in the text itself. These issues reduce transparency somewhat, but the manuscript remains substantially followable.

Novelty3/5

high confidence- spread 1- panel

Log-periodic phenomena and discrete scale invariance are well-established (Sornette refs cited). Time-varying boundary conditions and dynamical Casimir effects are also known. The novel contribution is the specific numerical demonstration that explicit log-periodic boundary modulation imprints a measurable log-frequency hierarchy that survives propagation through a finite cavity, distinguishable from harmonic/random drives. This is a useful but incremental synthesis rather than a fundamentally new mechanism. The appendix showing generic behavior across b values further reduces novelty of the specific 13/8 choice — it suggests the 13/8 is arbitrary rather than physically motivated within this paper, raising questions about why this ratio is highlighted.

Completeness4/5

high confidence- spread 1- panel

The work is well-structured and addresses its stated goals comprehensively. The FDTD implementation is clearly specified with grid parameters, boundary conditions, and time stepping. The post-processing pipeline is detailed, including windowing, envelope fitting, and residual extraction. Statistical analysis includes correlation coefficients, significance testing via surrogate data, and systematic parameter sweeps. Control simulations (undriven, harmonic, random) provide proper baselines. Energy conservation verification adds robustness. The methodology section provides sufficient detail for reproduction, supported by promised code deposit. Minor gaps include: some technical details about the Mur absorbing boundaries are abbreviated, and the physical interpretation of why this specific scaling ratio was chosen could be stronger. However, these don't affect the core argument's validity.

Publication criteria: All dimensions must score at least 2/5 with an overall average of 3/5 or higher. The AI recommendation badge above is advisory - publication is determined by the numerical scores.

This paper presents a methodologically rigorous numerical study demonstrating that explicit log-periodic boundary permittivity modulation in a 1D FDTD electromagnetic cavity produces a measurable spectral hierarchy that is periodic in log-frequency with the imposed scaling period ln(13/8). The work is internally consistent, employs proper controls (undriven, harmonic, random), and provides strong statistical validation through surrogate testing (r ≈ 0.81, p < 0.001). The mathematical setup using standard 1D Yee FDTD with stable Courant factor is sound, though the central theoretical gap is the absence of an analytic derivation linking the time-domain log-periodic modulation to the predicted frequency-domain residual form. The specialists flagged specific mathematical risks where derivations are compressed or assumed rather than proven, particularly Equations 1-2 which present the geometric peak spacing and residual form as 'predicted' without showing the transfer function from boundary modulation to spectral response. The work succeeds as a numerical phenomenology demonstrating a controlled electromagnetic effect, with comprehensive robustness checks including energy conservation verification, probe-position dependence, and windowing sensitivity tests. While the effect appears generic across scaling ratios (Appendix A), this somewhat undermines the specific focus on b = 13/8, whose physical motivation is never explained.

This review was generated by AI for research and educational purposes. It is not a substitute for formal peer review. All analyses are advisory; publication decisions are based on numerical score thresholds.

Key Equations (3)

\displaystyle \omega_n = \omega_0\, b^{n}, \quad b=\frac{13}{8}

Predicted geometric (geometric progression) spacing of spectral peaks under discrete scale invariance; equivalent to constant spacing \Delta\ln\omega=\ln b in log-frequency.

\displaystyle R(\omega)=A\cos\!\left(\frac{2\pi}{\ln(13/8)}\ln\!\left(\frac{\omega}{\omega_0}\right)+\phi\right)

Fitting model for the log-frequency residual after subtraction of the smooth power-law envelope; predicts oscillations periodic in ln(\omega) with period \ln(13/8).

\displaystyle \varepsilon(x,t)=1+\beta_{*}\cos\!\left(\frac{2\pi}{\ln(13/8)}\ln\!\left(\frac{t+t_0}{t_0}\right)+\phi\right)

Time-dependent log-periodic permittivity modulation applied in the boundary layers (cells 0--39 and 560--599) that drives the system; t_0 is the reference time and \beta_{*} the modulation depth.

Other Equations (2)

\displaystyle S(\omega)=|\tilde{E}(\omega)|^{2}

Definition of the probe-point power spectrum used as the primary observable.

\displaystyle S_{\rm env}(\omega)\propto\omega^{\alpha}

Fitted smooth spectral envelope (power-law) subtracted from S(\omega) to obtain the residual R(\omega).

Testable Predictions (3)

A boundary log-periodic permittivity modulation with scaling ratio b = 13/8 imprints a discrete-scale hierarchy on the probe-point power spectrum such that the residual after power-law subtraction oscillates periodically in ln(\omega) with period ln(13/8).

otherpending

Falsifiable if: If repeated FDTD runs with the same boundary modulation and analysis fail to produce a residual whose fitted period matches ln(13/8) within the reported fit uncertainty, or if identical control protocols (undriven/harmonic/random) produce residuals with comparable amplitude/correlation, the claim is falsified.

The log-periodic spectral imprint is robust across modulation depths and phase (within the reported sweep), and statistical significance (r \approx 0.81, p < 0.001) is attained relative to randomized-phase surrogates and control runs.

otherpending

Falsifiable if: If parameter sweeps over modulation depth \beta_{*} and phase \phi systematically fail to produce significantly larger residual amplitudes or correlations than controls, or if surrogate tests yield comparable |r| values (p \gtrsim 0.05), the robustness and statistical-significance claim is falsified.

The discrete-scale spectral hierarchy propagates coherently into the cavity interior: residual amplitude decays with distance from the driven boundary and plateaus near the cavity center at roughly 37% of the boundary value.

otherpending

Falsifiable if: If probe measurements at increasing interior distances show no coherent decay/plateau trend (e.g., residual amplitude falls to noise or does not plateau), or if interior probes display identical behavior in undriven controls, the claim of coherent propagation is falsified.

Tags & Keywords

classical electromagnetism(domain)discrete scale invariance(physics)dynamic boundary conditions(physics)finite-difference time-domain (FDTD)(methodology)log-periodic modulation(physics)power-law envelope fitting(math)spectral analysis(methodology)

Keywords: finite-difference time-domain (FDTD), log-periodic modulation, discrete scale invariance, boundary-driven electromagnetism, power spectrum, log-frequency residual, dynamic boundary conditions, Yee grid

Log-Periodic Spectral Hierarchies in a Boundary-Driven Electromagnetic Cavity: Evidence from FDTD Simulations Jill F. Rankin Independent Researcher (Dated: May 12, 2026) 1

Abstract We demonstrate that explicit log-periodic permittivity modulation at the boundaries of a one- dimensional electromagnetic cavity, implemented with fixed scaling ratio b = 13/8, imprints a reproducible hierarchy in the electromagnetic energy spectrum. The primary observable is the probe-point power spectrum S(ω) =| ̃ E(ω)| 2 . After subtraction of the smooth power-law envelope, the log-frequency residual exhibits oscillations whose period matches ∆ lnω = ln(13/8), consis- tent with the externally imposed drive period. This structure is absent in undriven, harmonically driven, and randomly modulated control simulations. The framework is strictly mode-accessibility engineering within linear electromagnetism. It does not modify fundamental constants or the un- derlying quantum field theory and reduces exactly to standard Maxwell theory for static boundaries and thermal equilibrium. I. INTRODUCTION Dynamic boundary conditions in classical wave systems can redistribute spectral energy without altering the underlying field equations. Time-dependent boundaries are central to cavity quantum electrodynamics, photonic structures, plasma sheaths, metamaterials, and dynamical Casimir systems [1–6]. The central question addressed here is whether explicitly scale-structured temporal boundary conditions can leave a robust and distinguishable spectral imprint after finite-wave propagation, modal interference, absorbing-boundary evolution, and Fourier reconstruction. Here we investigate a log-periodic boundary modulation with fixed scaling ratio b = 13/8 ≈ 1.625, framed strictly as mode-accessibility engineering within linear electromag- netism. The frequency-domain hierarchy emerges dynamically from time-domain propagation un- der the boundary-conditioned modulation. The log-periodic drive is predicted to generate peaks spaced geometrically as ω n = ω 0 b n , b = 13 8 ,(1) which appears as constant spacing ∆ lnω = lnb = ln(13/8) in logarithmic frequency. The strongest diagnostic is the residual after subtraction of the smooth spectral envelope. Discrete scale invariance and log-periodic corrections are well known in systems with hier- 2

FIG. 1. Example power spectrum from a log-periodically driven run, illustrating the emergence of a discrete-scale hierarchy (geometric spacing in ω) that becomes approximately periodic when plotted versus lnω. archical or self-similar structure [7, 8]. The present work studies a controlled numerical ques- tion: whether an explicitly imposed log-periodic boundary drive can imprint a measurable, reproducible, and control-distinguishable log-frequency hierarchy in a finite electromagnetic cavity. 3

II. METHODS A. FDTD Implementation A one-dimensional Yee-grid finite-difference time-domain (FDTD) solver is implemented on a uniform grid of N = 600 cells with spatial step ∆x = 1 and time step ∆t = 0.5 in dimensionless code units. This gives Courant factor c∆t/∆x = 0.5, satisfying the one- dimensional stability bound. First-order Mur absorbing boundaries are applied at both ends [9? ]. A spatially inhomogeneous, time-varying permittivity is applied only in the left and right boundary layers, cells 0–39 and 560–599 (Fig. 2): ε(x,t) = 1 + β ∗ cos 2π ln(13/8) ln t + t 0 t 0

φ , with t 0 = 1. The system is excited by a broadband Gaussian pulse injected near the left boundary, and the electric field E z is recorded at multiple interior probes for 16000 timesteps. FIG. 2. Schematic of the one-dimensional FDTD cavity and boundary-driven modulation protocol. The time-varying permittivity is applied only within the left and right boundary layers, while the field is recorded at interior probe locations. B. Post-Processing Pipeline A Hann window is applied to the recorded time series prior to Fourier transformation. Fits are performed over the interval ω ∈ [0.1, 20] in code units. The smooth envelope S env (ω) ∝ ω α is obtained by least-squares fitting in log-log space. The fitted exponent α is 4

stable to ±0.05 across all runs and is insensitive to modest changes in the fitting-interval endpoints (tested over [0.05, 25] versus [0.1, 20]); the resulting residual amplitude A changes by less than 5%. Spectral leakage robustness was confirmed with a Kaiser window (β = 8); the modulated-run residual amplitude changed by less than 4%. The relative residual R(ω) is fit to the predicted form R(ω) = A cos 2π ln(13/8) ln ω ω 0

φ ,(2) with fixed period and ω 0 = 1 (arbitrary reference scale). The fitted amplitude A, uncertainty, and Pearson correlation coefficient r are reported. Control runs use identical procedures. A systematic 5× 5 parameter sweep over β ∗ ∈ [0.01, 0.20] and φ∈ [0, 2π] is performed for both modulated and control cases. C. Numerical Robustness The log-periodic structure persists throughout later stages of the simulation. To quantify propagation effects, the residual amplitude was extracted at multiple probe locations (cells 50, 150, 300, 450). The hierarchy amplitude decays continuously with distance from the driven boundary layer and plateaus near the cavity center. Statistical significance of the observed correlations was assessed via 1000 randomized- phase surrogate time series generated from the modulated runs (preserving the power spec- trum but randomizing phases to destroy coherent log-periodic structure). None of the sur- rogates produced |r| > 0.25 under the identical fixed-period fit, establishing p < 0.001 for the modulated-run correlations (r = 0.81± 0.04). Energy conservation is verified by tracking the total electromagnetic energy and cumula- tive boundary work (see Results). All simulations were performed using dimensionless code units and fully deterministic update rules. D. Null Hypothesis The null hypothesis—no effective discrete-scale spectral imprint—is rejected if the resid- ual R(ω) exhibits a reproducible oscillation exceeding control-run amplitudes, the extracted 5

period matches ln(13/8) within fit uncertainty, and the structure is absent in all control simulations. III. RESULTS Figure 5 shows the power spectrum and log-frequency residual for a representative run with β ∗ = 0.12 and φ = π/4. The fitting interval spans ≈ 11 full periods of the predicted residual. Quantitative residual statistics from the full 5× 5 sweeps are summarized in Table I. Modulated runs (β ∗ ≥ 0.08) yield mean Pearson r = 0.81± 0.04; all control runs across identical sweeps yield mean |r| = 0.13± 0.06 with maximum |r| = 0.21. TABLE I. Residual fit statistics (full 5× 5 sweeps). Drive TypeMean Amplitude A Mean Pearson r (max |r|) Log-periodic (β ∗ ≥ 0.08)0.079± 0.0120.81± 0.04 Undriven (β ∗ = 0)0.008± 0.0110.12± 0.05 (0.19) Harmonic0.017± 0.0140.15± 0.06 (0.21) Random0.010± 0.0130.11± 0.05 (0.18) Figure 3 visualizes the 5× 5 parameter sweep, showing that strong, coherent residuals occur only in the log-periodically driven case and grow with modulation depth β ∗ , while control protocols remain near the noise floor. Energy conservation is explicitly verified (U (t) + W boundary (t) constant to better than 0.1%; see Supplemental Material for the time series). IV. DISCUSSION The observed spectral hierarchy demonstrates that explicit log-periodic boundary modu- lation can impose a reproducible geometric structure on the electromagnetic mode spectrum after full cavity propagation. The hierarchy amplitude decays continuously with probe dis- tance from the boundary layer (Fig. 6) and plateaus at ∼ 37% of the direct-drive value, confirming coherent transmission rather than localized near-field effects. 6

FIG. 3. Summary of the full 5× 5 sweep over modulation depth β ∗ and phase φ. The log-periodic drive produces consistently larger fitted residual amplitudes and correlations than undriven, har- monic, or random controls. While the presence of the hierarchy is generic across scaling ratios (Supplemental Ma- terial), the fitted amplitude is only weakly sensitive to the precise value of b in the tested range. These results are consistent with the proposed boundary-driven spectral imprint- ing mechanism. No alternative mechanism producing a fixed log-periodic residual in the absence of a log-periodic drive has been identified. The present study is purely numerical and phenomenological; no closed-form analytic derivation of the residual transfer function is claimed. Future work with longer simulations or narrower-band analysis may improve specificity. Physical realizations could include plasma density modulation, varactor-loaded metama- terial boundaries, or optomechanical systems [10–13][14]. For example, the representative 7

FIG. 4. Representative comparison between the log-periodically modulated run and control simu- lations (undriven, harmonic, and random modulation) under identical analysis. The log-frequency residual oscillation is prominent only for the log-periodic drive. modulation depth β ∗ ≈ 0.12 corresponds to density fluctuations of order 10–20% in a GHz- range plasma sheath or capacitance tuning of ∼ 10% in a varactor-loaded metamaterial wall—both within reach of current RF and microwave technology. V. CONCLUSION Boundary-conditioned FDTD simulations with explicit log-periodic modulation at ratio 13/8 produce a reproducible spectral hierarchy periodic in lnω. This structure is absent under conventional drives, providing a falsifiable numerical demonstration of spectral redis- tribution via dynamic scaling symmetry within standard electromagnetism. ACKNOWLEDGMENTS The complete Python implementation, raw data, analysis scripts, and full reproduction package (including all parameter sweeps, probe scans, surrogate tests, and auxiliary runs) will be deposited in a public Zenodo repository prior to publication. This ensures immediate 8

FIG. 5. Power spectrum S(ω) and log-frequency residual R(ω) from the log-periodically modulated FDTD simulation with β ∗ = 0.12 and φ = π/4. The residual is fit to the cosine form with fixed period ln(13/8). independent reproducibility. [1] E. M. Purcell, Phys. Rev. 69, 681 (1946). [2] L. Novotny and B. Hecht, Principles of Nano-Optics, 2nd ed. (Cambridge University Press, 2012). [3] E. N. Economou, Green’s Functions in Quantum Physics, 3rd ed. (Springer, 2006). [4] I. Carusotto et al., Nat. Phys. 16, 268 (2020). [5] V. V. Dodonov, Physics 2, 67 (2020). [6] C. M. Wilson et al., Nature 479, 376 (2011). 9

FIG. 6. Residual amplitude A as a function of probe position (cells 50, 150, 300, 450). The hierarchy decays continuously with distance from the driven boundary layer (cells 0–39) and plateaus near the cavity center, demonstrating coherent propagation through the finite dispersive cavity rather than near-field imprinting. [7] D. Sornette, Phys. Rep. 297, 239 (1998). [8] D. Sornette, Critical Phenomena in Natural Sciences: Chaos, Fractals, Selforganization and Disorder: Concepts and Tools, 2nd ed. (Springer, 2006). [9] D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, 2nd ed. (IEEE Press, 2013). [10] J. B. Pendry, D. Schurig, and D. R. Smith, Science 312, 1780 (2006). [11] C. Caloz and T. Itoh, Electromagnetic Metamaterials: Transmission Line Theory and Mi- crowave Applications (Wiley-IEEE Press, 2006). [12] A. M. Shaltout, V. M. Shalaev, and M. L. Brongersma, Science 364, eaat3100 (2019). [13] E. Lustig, Y. Sharabi, and M. Segev, Optica 5, 1390 (2018). [14] D. N. Basov, M. M. Fogler, and F. J. Garc ́ıa de Abajo, Science 354, aag1992 (2016). Appendix A: Genericity of the scaling ratio b (Supplemental Material) To demonstrate that the observed hierarchy is generic, auxiliary simulations were performed with irrational b = √ 2 ≈ 1.414 and b = 1.9 (identical grid, source, duration, β ∗ = 0.12, and post-processing). Interior-probe residuals yield r = 0.79 and r = 0.82, respectively, statistically indistinguishable from the main-text result. 10

While the presence of the hierarchy is generic, the fitted amplitude is only weakly sensitive to the precise value of b in the range [1.4, 1.9] (∆r < 0.05). All code required to reproduce these runs is included in the Zenodo deposit. Appendix B: Energy Conservation (Supplemental Material) Time evolution of total electromagnetic energy U (t) and cumulative boundary work W boundary (t) confirms that their sum remains constant to better than 0.1% for both mod- ulated and control runs (plot available in Zenodo repository). ... 11

You Might Also Find Interesting

Semantically similar papers and frameworks on TOE-Share

Finding recommendations...

← All Papers

Log-Periodic Spectral Hierarchies in a Boundary-Driven Electromagnetic Cavity: Evidence from FDTD Simulations

AI Review Scores

Derivation Flags

Overall Assessment

Strengths & Improvements

Specialist Agent Reports

Science Highlights

Key Equations (3)

Testable Predictions (3)

Tags & Keywords

Full Content

Related Papers

You Might Also Find Interesting

Community Discussion