Full-scatter vector field analysis of an overmoded and periodically-loaded cylindrical structure for the transportation of THz radiation

Highly overmoded and periodically loaded structures, such as the iris-line waveguide, offer an attractive solution for the efficient transportation of diffraction-prone THz pulses over long distances (hundreds of meters). This paper presents the full-scatter field theory that allows us to analytically derive all the spectral (modal) coefficients on the discontuities of the iris line. The spectral analysis uses vector fields, superseding scalar field descriptions, to account for diffraction loss as well as polarization effects and ohmic loss on practical conductive surfaces. An advanced application of Lorentz’s reciprocity theory, using a generalized guided-field configuration, is developed to reduce complexity of the mode-matching problem over nonuniform sections. The used technique is quite general and applies to a wide class of structures, as it only assumes a paraxial incidence (i.e. a parabolic wave equation) along the axis of the structure. It removes the traditional assumption of very thin screens, allowing for the study of thicker screens in the high-frequency limit, while formulating the problem efficiently by scattering matrices whose coefficients are found analytically. The theory agrees with and expands previously established techniques, including Vainstein’s asymptotic limit and the forward-scatter approximation. The used formulation also facilitates accurate visualization of the transient regime at the entrance of the structure and how it evolves to reach steady state.

💡 Research Summary

This paper presents a comprehensive analytical framework for the transport of terahertz (THz) radiation through highly over‑moded, periodically‑loaded cylindrical structures, commonly known as iris‑line waveguides. The authors address a critical challenge in long‑distance THz delivery: the simultaneous presence of strong diffraction, finite‑thickness metallic screens, and ohmic losses, which render conventional single‑mode or scalar‑field models inadequate.

Starting from Maxwell’s equations, the total electric and magnetic fields in both the narrow “waveguide” sections (radius a) and the wider “cavity” sections (radius r₀) are expanded in complete sets of transverse‑electric (TE) and transverse‑magnetic (TM) vector modes. Each mode is expressed with Bessel functions J₁ and its derivative, and the forward (+z) and backward (‑z) traveling components are retained explicitly, leading to equations (1)–(12). By normalizing each modal field (denoted with a hat) the authors obtain compact forms (13)–(18) that make the symmetry between forward and backward waves evident.

The central difficulty lies in matching these modal expansions at the discontinuities where the radius changes (step‑in and step‑out). Traditional approaches either assume identical cross‑sections on both sides of the discontinuity or invoke a forward‑scatter approximation that neglects backward reflections. To overcome this, the paper combines three powerful theorems:

1. Lorentz Reciprocity – generalized to a guided‑field configuration, allowing the inner product of two field solutions to be expressed as a surface integral over the discontinuity.
2. Schelkunoff’s Equivalence Theorem – replaces the physical discontinuity by an equivalent electric (and possibly magnetic) surface current J = n̂ × H₀ δ(z). This converts the original boundary‑value problem into two independent scattering problems on the left and right of the plane z = 0.
3. Uniqueness Theorem – guarantees that, given the equivalent surface currents, the fields in each region are uniquely determined, so the modal coefficients can be solved analytically.

Applying these theorems, the authors derive a scattering matrix S that relates the forward and backward modal amplitudes on either side of each discontinuity. The matrix elements are expressed in closed form as functions of the screen thickness δ, the radii a and r₀, the operating frequency (through k = ω/c), and the material conductivity (through the surface impedance). Importantly, the formulation does not require the thin‑screen assumption; it remains valid when δ ≈ λ or larger, thereby capturing both diffraction‑dominated loss (thin screens) and ohmic loss (thick, finite‑conductivity screens).

Numerical examples illustrate the theory. By varying the normalized thickness δ/λ from 0.1 to 1.0, the authors show that diffraction loss dominates for thin screens, reproducing Vainstein’s asymptotic results, while ohmic loss becomes significant as the screen thickens, leading to a rapid drop in transmission efficiency. The transient regime at the entrance of the iris line is visualized: the first few cells (tens of centimeters) exhibit strong mode conversion and field reshaping, after which the field settles into a steady‑state modal distribution that propagates with minimal additional loss.

The paper validates the new full‑scatter model against previously published forward‑scatter and eigen‑mode approaches, confirming that the latter are special cases of the present theory. It also demonstrates that the analytical scattering matrices dramatically reduce computational effort compared with full‑wave finite‑element simulations, while retaining high accuracy.

In conclusion, the authors deliver a versatile, vector‑field‑based scattering formalism that bridges the gap between scalar diffraction theory and full‑wave numerical methods for over‑moded THz waveguides. By incorporating both forward and backward waves, finite screen thickness, and material losses, the model provides a reliable tool for the design and optimization of long‑distance THz transport systems, such as those required for pump‑probe experiments at large accelerator facilities. Future extensions could address non‑paraxial incidence, non‑cylindrical geometries, and nonlinear material effects, further broadening the applicability of the framework.