arxiv: 2604.03798 · v1 · submitted 2026-04-04 · ⚛️ physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

From Wave Scattering to Bloch Bands: A Time-Domain Approach to Band Formation in Periodic Media

Amit Tanwar, Nishant Kashyap, Pragati Ashdhir, Vivek T. Ramamoorthy

Authors on Pith no claims yet

Pith reviewed 2026-05-13 17:34 UTC · model grok-4.3

classification ⚛️ physics.comp-ph

keywords Bloch bandstime-domain simulationperiodic mediawave scatteringtransmission spectrumfinite-difference time-domainband gapselastic waves

0 comments

The pith

Time-domain wave propagation through finite periodic layers reconstructs the Bloch dispersion relation of the infinite system.

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

The paper shows how to obtain band structures by sending a broadband pulse through a finite stack of layers and reading the transmission spectrum, instead of solving an eigenvalue problem in reciprocal space. A finite-difference time-domain calculation tracks the elastic wave as it travels, accumulates phase, and attenuates inside gaps, turning the abstract Bloch bands into visible consequences of scattering and interference. Students therefore see band gaps emerge as frequency windows of strong reflection rather than as static lines on a diagram, while the same simulation immediately displays finite-size shifts, disorder effects, and defect modes.

Core claim

Band formation is recovered directly from the transmission spectrum of a finite periodic structure, where the frequency-dependent phase delay and attenuation encode the dispersion relation of the corresponding infinite lattice through the accumulated wavevector and the spatial decay inside gaps.

What carries the argument

Staggered-grid finite-difference time-domain propagation of elastic waves through a finite layered medium, followed by extraction of the Bloch wavevector from the frequency-dependent transmission phase and identification of gaps from spatial attenuation.

If this is right

Band gaps appear as clear frequency intervals of strong spatial attenuation caused by destructive interference after many scattering events.
Finite-size corrections, weak disorder, and point defects can be added inside the same simulation to observe how they shift or fill the ideal bands.
Localized defect modes show up as sharp transmission peaks inside the gaps, directly linking scattering physics to mode confinement.
The identical workflow applies to any linear wave equation, supplying a single numerical route from scattering to band diagrams across acoustics, electromagnetics, and elasticity.

Where Pith is reading between the lines

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

The same finite-stack transmission data could be measured in a laboratory with ultrasonic pulses through stacked plates, providing an experimental test of the reconstruction.
Once the dispersion is obtained this way, small changes in layer thicknesses or material contrast can be varied to map design rules for desired gap locations without ever invoking reciprocal-space methods.
The approach supplies a concrete computational bridge that lets students move from time-domain intuition to reciprocal-space formalism rather than encountering the latter first.

Load-bearing premise

The transmission spectrum recorded in a finite-length stack faithfully reproduces the dispersion curve of the infinite periodic medium without large boundary or length-induced distortions.

What would settle it

Calculate the exact Bloch bands of the infinite layered medium with the standard transfer-matrix method and compare them with the dispersion extracted from the finite time-domain run; a mismatch larger than numerical error in band-edge locations or gap widths would falsify the reconstruction.

Figures

Figures reproduced from arXiv: 2604.03798 by Amit Tanwar, Nishant Kashyap, Pragati Ashdhir, Vivek T. Ramamoorthy.

**Figure 2.** Figure 2: FIG. 2. Comparison of narrowband and broadband source [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. Normalized numerical phase velocity [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Workflow of the time-domain simulation. The lay [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Snapshots of the velocity field [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Two-unit-cell layered structure embedded in water. [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. Top: Input (dashed) and transmitted (solid) spectra [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8. Smoothed transmission (in dB) for a 15-unit-cell [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗

**Figure 9.** Figure 9: FIG. 9. Spatial amplitude inside a ten-unit-cell periodic s [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗

**Figure 10.** Figure 10: FIG. 10. Comparison between the analytical Rytov dispersio [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗

**Figure 11.** Figure 11: FIG. 11. Transmission spectrum of the quarter-wavelength [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗

**Figure 12.** Figure 12: FIG. 12. Transmission spectrum of a six-cell quarter [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗

read the original abstract

Band formation in periodic media is a central topic in undergraduate solid-state physics, typically introduced through Bloch's theorem as an eigenvalue problem in reciprocal space for infinitely periodic systems. While mathematically elegant, this formulation can appear abstract: it assumes an idealized infinite lattice, shifts attention away from real-space wave dynamics, and presents band structures as static results rather than emergent consequences of wave propagation. Consequently, students often struggle to relate band gaps to familiar physical phenomena such as reflection, transmission, and interference, leading to a disconnect between formal band theory and observable wave behavior. We present a computational framework that addresses this gap by reconstructing band formation directly from time-domain wave propagation in finite periodic systems. Using a staggered-grid finite-difference time-domain scheme for elastic waves, a broadband excitation is propagated through a layered medium to obtain its transmission spectrum. From this, students extract the Bloch dispersion relation and observe spatial attenuation in band-gap regions, revealing the roles of multiple scattering and phase coherence. This approach provides a physically transparent pathway to band theory and enables exploration of finite-size effects, disorder, and defect-localized modes within a unified computational framework. Implemented through compact code and guided exercises, the method offers an accessible and versatile pedagogical tool, while also equipping students with transferable skills in numerical modeling of wave phenomena across disciplines.

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.

Desk Editor's Note private letter to a colleague

This paper gives a practical FDTD-based teaching method to show band formation emerging from wave propagation in finite periodic stacks, but the extraction of clean infinite-system dispersion needs explicit checks against known results.

read the letter

The main point is that the authors describe a time-domain route to band theory: run staggered-grid FDTD on a finite layered medium, send a broadband pulse, measure the transmission spectrum, and pull out the dispersion relation plus the spatial attenuation in gaps. This is meant to make Bloch bands feel like a consequence of scattering and interference rather than an abstract eigenvalue problem. The framework also folds in finite-size effects, disorder, and defect modes in the same setup, and they supply compact code plus exercises. That combination is genuinely useful for classroom work. It lowers the barrier for students to connect real-space dynamics to the usual reciprocal-space picture. The code description and the unified treatment of several phenomena are the parts that stand out as well executed. The soft spot is the central mapping step. Transmission through a finite stack carries Fabry-Perot oscillations inside passbands and rounded gap edges, so recovering the exact infinite-lattice dispersion requires either large enough N or a correction procedure. The abstract gives no convergence data with period count and no side-by-side comparison to the transfer-matrix Bloch solution. If the full paper shows those checks and the extracted bands match the known analytic result to within a few percent, the method is reliable for teaching. If the match is only qualitative, the pedagogical claim still holds but the accuracy part is weaker than stated. This is written for instructors who want ready computational labs in wave physics or introductory solid-state courses. A reader looking for concrete exercises that link scattering to bands will get immediate value from the code and the guided examples. It deserves a serious referee because the implementation details determine whether the approach actually works in practice and because education-focused journals need this kind of material. I would send it out for review rather than desk-reject.

Referee Report

2 major / 2 minor

Summary. The paper claims to introduce a pedagogical computational framework that reconstructs Bloch band formation directly from time-domain elastic wave propagation in finite periodic layered media. Using a staggered-grid FDTD scheme, a broadband pulse is transmitted through a finite stack to obtain the transmission spectrum T(ω), from which the dispersion relation ω(k) is extracted; students are then shown spatial attenuation in gaps and finite-size effects, providing a real-space alternative to the conventional reciprocal-space Bloch eigenvalue problem.

Significance. If the extraction procedure is shown to recover the infinite-lattice dispersion with controlled error, the work could furnish a concrete, visually intuitive route from scattering and interference to band gaps that is currently missing from most undergraduate solid-state curricula. The emphasis on finite-N systems, disorder, and defect modes, together with compact code and guided exercises, would make the material immediately deployable across physics and engineering courses while also supplying transferable numerical skills.

major comments (2)

[Abstract and §3] Abstract and §3 (extraction procedure): the central claim that the transmission spectrum of a finite-N stack directly yields the Bloch dispersion of the corresponding infinite lattice is not supported by any convergence test with N, any comparison against the exact transfer-matrix Bloch solution, or an explicit formula linking T(ω) to k(ω) that is proven to eliminate Fabry–Perot oscillations and interface reflections.
[§4] §4 (numerical results): no error analysis or quantitative comparison to known analytical dispersion for a simple 1D bilayer is presented; without such benchmarks it is impossible to assess whether finite-size resonances distort the extracted bands at a level that would undermine the pedagogical goal.

minor comments (2)

[Abstract] The abstract would be clearer if it stated the precise algorithm (phase accumulation, fitting, or attenuation mapping) used to convert the transmission spectrum into the dispersion curve.
Figure captions should explicitly indicate the number of periods N used in each simulation so that readers can judge the scale of finite-size effects.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments correctly identify the need for stronger validation of the extraction procedure and quantitative benchmarks. We have revised the manuscript to incorporate explicit convergence tests, comparisons with the transfer-matrix method, and error analysis, as described in the point-by-point responses below.

read point-by-point responses

Referee: [Abstract and §3] Abstract and §3 (extraction procedure): the central claim that the transmission spectrum of a finite-N stack directly yields the Bloch dispersion of the corresponding infinite lattice is not supported by any convergence test with N, any comparison against the exact transfer-matrix Bloch solution, or an explicit formula linking T(ω) to k(ω) that is proven to eliminate Fabry–Perot oscillations and interface reflections.

Authors: We accept the criticism that the original manuscript lacked explicit validation. In the revised version, §3 now contains a derivation of the relation between the transmission phase and the Bloch wavevector k(ω) for large but finite N, together with a windowing procedure that suppresses Fabry–Perot oscillations. We also add convergence plots of the extracted dispersion versus N, directly compared to the exact transfer-matrix Bloch solution for the same layered stack. These additions furnish the missing support for the central claim while preserving the pedagogical emphasis on finite systems. revision: yes
Referee: [§4] §4 (numerical results): no error analysis or quantitative comparison to known analytical dispersion for a simple 1D bilayer is presented; without such benchmarks it is impossible to assess whether finite-size resonances distort the extracted bands at a level that would undermine the pedagogical goal.

Authors: We agree that quantitative error analysis is required. The revised §4 now includes a direct comparison of the FDTD-extracted dispersion for a 1D bilayer against the analytical Bloch relation obtained from the transfer-matrix method. We report the frequency deviation (maximum and RMS) as a function of N and include figures showing convergence to the infinite-lattice limit. For N ≳ 20 the deviation in the pass bands falls below 1 %, confirming that finite-size effects do not undermine the pedagogical objectives. revision: yes

Circularity Check

0 steps flagged

No significant circularity; band extraction is forward simulation plus post-processing

full rationale

The paper's central step is a forward FDTD simulation of wave propagation through a finite periodic stack to obtain a transmission spectrum T(ω), followed by an extraction step that maps features of T(ω) (phase accumulation or attenuation) onto the Bloch dispersion ω(k). This is a data-driven reconstruction from an independent numerical experiment, not a self-definition, a fitted parameter renamed as prediction, or a load-bearing self-citation. No equations in the abstract or description reduce the output dispersion relation to the input ansatz by construction; the method remains falsifiable against transfer-matrix benchmarks or larger-N limits. The derivation chain is therefore self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The framework rests on standard numerical wave propagation methods and the physical equivalence between finite and infinite periodic responses in the appropriate limit.

axioms (1)

domain assumption Finite periodic systems exhibit transmission spectra that directly correspond to the band structure of infinite periodic media.
This is the key link assumed to allow extraction of Bloch bands from finite simulations.

pith-pipeline@v0.9.0 · 5545 in / 1162 out tokens · 61093 ms · 2026-05-13T17:34:30.718538+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear
Using a staggered-grid finite-difference time-domain scheme for elastic waves, a broadband excitation is propagated through a layered medium to obtain its transmission spectrum. From this, students extract the Bloch dispersion relation
IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
cos(k_B a) = D(ω) ... k''_B = (1/a) arcosh(|D(ω)|)

Reference graph

Works this paper leans on

12 extracted references · 12 canonical work pages

[1]

A clearer physical picture emerges when the motion is expressed in terms of ﬁelds representing momentum transport and restor- ing forces

does not explicitly reveal the dynamical mech- anism underlying wave propagation. A clearer physical picture emerges when the motion is expressed in terms of ﬁelds representing momentum transport and restor- ing forces. To make this structure explicit, it is conve- nient to recast the second-order equation into an equiva- lent ﬁrst-order system. Introduci...

work page
[2]

pro- vides the natural starting point for the numerical formu- lation developed in this work. In the following section, we construct a staggered-grid ﬁnite-diﬀerence scheme that respects this coupled structure and enables eﬃcient sim- ulation of wave propagation in media with spatially vary- ing material properties. III. A STAGGERED-GRID FORMULATION FOR E...

work page
[3]

5 cos ( π (t − t0) t0 ) + 0

42 − 0. 5 cos ( π (t − t0) t0 ) + 0. 08 cos ( 2π (t − t0) t0 ) , |t − t0| ≤ t0, 0, otherwise. (14) The sinc factor produces an approximately ﬂat spec- trum up to fcut, while the Blackman window conﬁnes the pulse in time and suppresses spectral side lobes. In practice, the choice of excitation is guided by the desired spectral content: the Ricker wavelet p...

work page
[4]

Start- ing from the deﬁnition of the layered structure and the choice of numerical parameters, the algorithm advances the coupled velocity and stress ﬁelds in time while in- jecting the source and applying absorbing boundaries at the edges of the computational domain. This sequence of steps mirrors the logical progression followed by stu- dents when imple...

work page 1960
[5]

Examine how the larger impedance contrast modiﬁes the partition of energy between reﬂected and transmitted waves

and compare them with the nu- merical values obtained from the simulation. Examine how the larger impedance contrast modiﬁes the partition of energy between reﬂected and transmitted waves. VIII. FROM REPEATED SCATTERING TO FREQUENCY-SELECTIVE TRANSMISSION Having established the scattering properties of a single impedance interface, we now examine how repe...

work page
[6]

(36) When |D(ω )| ≤ 1, a real solution for kB exists and the frequency lies within a pass band

reduces to cos(kBa) = D(ω ). (36) When |D(ω )| ≤ 1, a real solution for kB exists and the frequency lies within a pass band. If instead |D(ω )|> 1, no real kB satisﬁes the equation. Taking the modulus of Eq. ( 36) and writing kB = ik′′ B, with cos( ix) = cosh( x), gives |cos(kBa)|= cosh(k′′ Ba) = |D(ω )|, (37) so that k′′ B = 1 a cosh− 1( |D(ω )| ) . (38)...

work page
[7]

In stop bands, where |D(ω )|> 1, the Bloch wavevector becomes complex, yielding the attenuation constant k′′ B given in Eq

admits real solutions for kB, corresponding to propagating Bloch modes and a deﬁ- nite phase advance kBa across each unit cell. In stop bands, where |D(ω )|> 1, the Bloch wavevector becomes complex, yielding the attenuation constant k′′ B given in Eq. ( 38). We now examine both manifestations directly within the time-domain simulation and compare them wit...

work page
[8]

The close agreement conﬁrms that the spatial attenuation observed in the ﬁnite stack is governed by the same Bloch attenuation constant predicted for an in- ﬁnite periodic medium

05%. The close agreement conﬁrms that the spatial attenuation observed in the ﬁnite stack is governed by the same Bloch attenuation constant predicted for an in- ﬁnite periodic medium. In other words, the exponential decay measured directly in the time-domain simulation is not merely a numerical artifact of multiple reﬂections, but a manifestation of the ...

work page
[9]

Compute the trans- mission spectra and observe how the defect-mode fre- quency shifts as the defect geometry changes

5 times the original layer width). Compute the trans- mission spectra and observe how the defect-mode fre- quency shifts as the defect geometry changes. XI. CONCLUSION We have presented a compact computational frame- work through which band formation in periodic media can be reconstructed directly from time-domain wave dy- namics in ﬁnite layered systems....

work page
[10]

The present solver readily supports a variety of additional one-dimensional studies. Once imple- mented, students can modify layer thicknesses, ma- terial properties, or introduce defects to investigate band-gap tuning, localized defect modes, and the inﬂuence of impedance contrast on transmission in periodic layered media. 6,35–38

work page
[11]

In the present context, it is ap- plied to the velocity and stress ﬁelds of elastody- namics

The staggered-grid update scheme follows the leapfrog ﬁnite-diﬀerence time-domain (FDTD) method originally introduced by Yee for electro- magnetic waves. In the present context, it is ap- plied to the velocity and stress ﬁelds of elastody- namics. This connection provides a natural path- way to extend the framework to electromagnetic wave simulations; a c...

work page
[12]

Quantum mechanics of electrons in crystal lattices,

The one-dimensional model can also be extended to two dimensions. Such simulations enable explo- ration of richer phononic or photonic crystal phe- nomena, including directional band gaps, waveg- uiding, and defect cavities. 20,40,41 Since phononic crystals are an active area of research with applica- tions in vibration control and wave manipulation, the ...

work page 1931