THz Circular Waveguide Live Simulator
Live

Electromagnetic mode laboratory

THz Circular Waveguide Live Simulator

Explore Gaussian-beam coupling, cutoff, dispersion and copper-wall loss in an air-filled circular waveguide—directly in your browser.

Copper 2.0 mm ID 40 cm f = 100 mm
02

Transverse electric field

TE11 1.00 THz
x
y
Propagating
normalized Ex → transverse E direction copper boundary shown in amber
03

Spectral response

Drag or click the chart to set frequency
output power copper loss cutoff
05

Model used in the browser

fc = cχ / 2πa
β = √(k0² − kc²)
w0 = λf / πwL
Rs = √(πfμ0/σ)

Perfect-conductor circular eigenfields with first-order copper surface-resistance loss. Centered Gaussian overlap is evaluated numerically. Flanges, roughness, offsets and lens aberrations are not included.

Reference notes

Theory behind the simulation

The model combines ideal circular-waveguide eigenmodes, Gaussian beam optics, overlap-integral coupling, and a first-order copper surface-loss correction.

01

Circular-waveguide modes

In a uniform, air-filled circular guide, the longitudinal field satisfies the scalar Helmholtz equation. Separation in cylindrical coordinates gives radial solutions proportional to Bessel functions:

Ψmn(r,φ) = Jm(kcr) {cos(mφ), sin(mφ)}

TE modes have Ez = 0 and use roots χ′mn of J′m. TM modes have Hz = 0 and use roots χmn of Jm. The dominant circular-guide mode is TE11, with χ′11 = 1.84118.

02

Cutoff and propagation

For guide radius a, the transverse wavenumber is kc = χ/a. Each mode propagates only when the operating frequency exceeds its cutoff:

fc = cχ/(2πa),   β = √(k0² − kc²)

Above cutoff, β is real and power travels along the tube. Below cutoff, β is imaginary and the field is evanescent. Waveguide dispersion also makes the group velocity vg = c√(1 − (fc/f)²), producing delay relative to free space.

03

Gaussian focusing

A collimated Gaussian THz beam illuminating a thin lens is focused to an approximate 1/e electric-field waist:

w0 = λf/(πwL)

Here λ is the free-space wavelength, f is the lens focal length, and wL is the Gaussian radius on the lens. The simulator lets wL vary as a fraction of the lens radius. A shorter focal length or a more fully illuminated lens produces a smaller focused waist.

04

Mode coupling

Coupling is evaluated from the normalized transverse-field overlap between the centered, x-polarized Gaussian beam and the selected waveguide eigenmode:

η = |∬ EG·Emode* dA|² / (∬|EG|²dA ∬|Emode|²dA)

The result includes the fraction of Gaussian power entering the circular aperture. Symmetry can make the overlap of a centered beam with some modes zero even when those modes are above cutoff.

E(r,φ,L) = ATEETE11eiβTEL + ATMETM11ei(βTML+φ0)

In interference mode the two complex output fields are added coherently. Their output phase is φ0 + (βTM11 − βTE11)L, and the axial beat length is Lb = 2π/|βTE11 − βTM11|. Orthogonality makes total integrated modal power additive, while the local field pattern depends strongly on relative phase.

05

Copper-wall attenuation

Finite copper conductivity creates a surface resistance that increases with the square root of frequency:

Rs = √(πfμ0/σCu),   σCu = 5.8×10⁷ S m⁻¹

The wall-loss power is obtained from the tangential magnetic field at r = a. The resulting attenuation α gives conductor loss in dB m⁻¹ and output power Pout/Pin = η exp(−2αL).

06

Interpretation and limits

This is a rapid analytical design model, not a full 3D FEM or FDTD calculation. It assumes a straight, uniform circular tube, an ideal thin lens, an aligned centered beam, and smooth copper walls.

Not included: flange reflections, lens aberrations, beam tilt or offset, surface roughness, bends, discontinuities, dielectric windows, and mode conversion along the tube.

Use the results for parameter exploration, then validate a final design with measurement or a calibrated full-wave solver.

Symbols used

a waveguide radius L tube length χ, χ′ Bessel roots β propagation constant w0 focused waist wL beam radius at lens η modal coupling α field attenuation