The Project window shown below contains active interactive GUI and simulation tools. Interested in doing fun with terahertz light, alright! Go ahead!

Interactive scientific tools in terahertz light and a live platform to see the change in terahertz transient developed by-KARARWal Deepak! Have fun with the simulators!

PhD Projects in Terahertz optics

A collection of interactive simulations, theoretical tools, and educational resources developed around terahertz time-domain spectroscopy, polarization, structured vector beams, spatiotemporal electric-field analysis, and circular-waveguide propagation. Select a project from the vertical menu to read the physics, purpose, and main capabilities.

THz Optics Polarization Vector Beams Spatiotemporal Fields Waveguides Interactive Simulation
01

THz Time-Domain Pulse Explorer

From E(t) to spectrum, phase, and group delay

This project is an interactive introduction to the information contained in a terahertz time-domain waveform. A measured THz pulse is not only a pulse-amplitude trace: because coherent THz detection preserves the sign and temporal phase of the electric field, the waveform also contains information about arrival time, polarity, temporal width, ringing, delayed echoes, chirp, reflections, and noise.

The project connects the time-domain field E(t) to its Fourier-transform representation. The frequency-domain view reveals spectral amplitude, phase, bandwidth, and group-delay behavior. It also explains how the delay-stage step and total scan window determine the Nyquist frequency and the frequency-bin spacing.

Time domainElectric-field polarity, arrival time, peak field, pulse width, oscillations, ringing, echoes, and noise.
Frequency domainFourier amplitude, spectral phase, bandwidth, group delay, chirp, and spectral fringes.
SamplingTime step, total observation window, Nyquist limit, and spectral resolution.
Educational goalConnect features seen in experimental THz-TDS waveforms to their physical meaning.
E(t)  → Fourier transform →  Ē(f) = |Ē(f)| exp[iφ(f)]
02

Light Polarization Simulator

Build a polarization state from two orthogonal electric-field components

Polarization emerges from the relative amplitudes and phase of two orthogonal electric-field components. This simulator allows the user to vary the x- and y-component amplitudes and their relative phase and immediately observe how the resulting polarization trajectory changes between linear, circular, and elliptical states.

The same electromagnetic field is displayed through face, side, and top views so that the polarization ellipse can be related directly to the propagating wave. The simulator also evaluates normalized Stokes parameters and displays the polarization orientation, ellipticity, and handedness.

Amplitude controlChange Ax and Ay to control the principal axes of the polarization ellipse.
Phase controlChange the relative phase δ to move continuously from linear to elliptical or circular polarization.
Stokes analysisVisualize s1, s2, and s3 for the selected polarization state.
Propagation viewsCompare the xy polarization trajectory with xz and yz views of the propagating field.
Ex = Ax cos(kz − ωt),    Ey = Ay cos(kz − ωt + δ)
03

Broadband THz Radial and Azimuthal Vector Beams

Structured polarization across a broadband THz beam

This interactive project explores radially and azimuthally polarized terahertz vector beams. Unlike a uniformly polarized beam, a vector beam has a polarization state that changes across the transverse plane. The local electric-field orientation therefore carries spatial structure in addition to the usual temporal and spectral structure of a broadband THz pulse.

The simulator presents spatial electric-field vectors together with Stokes-parameter maps, allowing the user to connect the local vector field to S0, s1, s2, and s3. It is designed to illustrate how radial and azimuthal polarization appear in spatially resolved full-Stokes measurements.

Radial beamLocal linear polarization points along the radial direction from the beam center.
Azimuthal beamLocal linear polarization follows the tangential direction around the beam center.
Stokes mapsInspect spatial variations of intensity and normalized polarization parameters.
Broadband contextConnect spatial polarization structure with frequency-resolved THz characterization.
The project is intended as an interactive visualization of structured THz polarization and full-Stokes mapping rather than a replacement for a complete experimental model.
04

Spatiotemporal THz Electric-Field Explorer

From delay-stage scanning to E(x,y,t) and E(x,y,f)

This project explains how a sequence of two-dimensional THz electric-field images becomes a three-dimensional spatiotemporal dataset. At each delay position the experiment records a field map; arranging those frames in temporal order constructs E(x,y,t). A Fourier transform along the time axis at every pixel then produces the complex frequency-domain field E(x,y,f).

The explorer also calculates practical sampling quantities such as the mechanical delay-stage scan length, time step, total time window, frequency resolution, Nyquist frequency, number of temporal frames, number of positive-frequency frames, pulse width, and spectral bandwidth.

Delay-stage physicsConvert mechanical displacement into temporal delay, including double-pass retroreflector geometry.
Sampling limitsCalculate Δt, total time window, Δf, Nyquist frequency, and number of FFT frames.
Field cubeUnderstand how camera frames form Ex(x,y,t) and Ey(x,y,t).
Spectral imagingApply a pixel-wise FFT to obtain frequency-resolved amplitude and phase maps.
Δt = mΔL/c,    Δf = 1/(NΔt),    fNyquist = 1/(2Δt)
05

Temporal Evolution of a Terahertz Electric Field

Mathematical theory of coherent spatiotemporal THz detection

This theory-focused resource develops the mathematical link between a mechanical optical delay line and the measured evolution of a broadband THz electric field. It explains why each translation-stage position corresponds to one temporal sample and how a sequence of sampled spatial frames reconstructs Ex(x,y,t) and Ey(x,y,t).

The treatment continues into the frequency domain, where pixel-wise Fourier transformation provides complex amplitude and phase. It also introduces polarization analysis and shows how frequency-resolved electric fields can be interpreted using spatial Stokes parameters and polarization quantities.

Optical delayStage displacement, optical-path difference, and conversion to relative time.
Temporal samplingRelationship among sample spacing, observation window, frequency spacing, and Nyquist bandwidth.
Complex fieldFourier amplitude and phase evaluated independently at each spatial coordinate.
PolarizationConnection from Ex and Ey to Stokes parameters and spatially resolved polarization.
τ(x) = 2x/c,    Ēx,y(x,y,f) = Ft{Ex,y(x,y,t)}
06

THz Circular Waveguide Live Simulator

Gaussian-beam coupling and modal propagation in a circular copper waveguide

This browser-based simulator models the coupling of a focused Gaussian terahertz beam into an air-filled circular copper waveguide. The default configuration uses a 2 mm inner diameter, a 40 cm copper tube, a 100 mm focal-length lens, a 50.8 mm clear aperture, and a 0.10–2.00 THz frequency range.

The simulator evaluates circular-waveguide TE and TM modes, their cutoff frequencies, propagation constants, focused Gaussian waist, Gaussian-to-mode overlap, first-order copper conductor loss, transmitted modal power, relative group delay, and spectral response. Controls update the results directly in the browser.

ModesExplore TE11, TM01, TE21, TE01, TM11, TE12 and their frequency-dependent behavior.
CouplingRelate lens geometry and Gaussian beam size to focal waist and modal overlap.
LossEstimate conductor loss using the copper surface-resistance model.
DispersionCompare propagation constants and relative group delays of supported modes.
fc = cχ/(2πa),    β = √(k0² − kc²),    w0 = λf/(πwL)
The model is a fast analytical/semi-analytical tool for physical interpretation and parameter exploration. It is not a full 3D FDTD/FEM solver and does not include all experimental imperfections.
07

THz Circular Copper Waveguide Modes

Semi-analytical broadband modal model for an air-filled copper guide

This repository provides a reproducible semi-analytical treatment of a linearly polarized Gaussian THz beam coupled through a lens into a circular copper waveguide. The default system uses a 2.0 mm inner diameter, 40 cm guide length, copper conductivity of 5.8 × 10⁷ S m⁻¹, a 100 mm focal-length lens, and a 50.8 mm aperture.

The model calculates TE and TM cutoff frequencies from Bessel-function roots, frequency-dependent propagation constants and relative group delays, full vector eigenfields, diffraction-limited Gaussian focusing, aperture clipping, Gaussian-to-mode power overlap, copper-wall conductor loss, modal power at the input and output, and coherent output field maps at selected THz frequencies.

EigenmodesFull vector Ex, Ey, Ez, Hx, Hy, and Hz fields for circular-waveguide TE/TM solutions.
Cutoff physicsBessel-function roots determine which modes can propagate at each frequency.
Power evolutionTrack coupled modal power and its reduction after propagation through copper.
Output fieldsConstruct coherent multimode field maps at representative THz frequencies.
Rs = √(πfμ0/σ)    — copper surface resistance used in the first-order conductor-loss model
The repository is intended for fast interpretation and parameter sweeps. Real insertion loss can be larger because flanges, roughness, oxidation, misalignment, lens aberration, and other experimental effects are not included in the idealized model.
These projects are designed for scientific visualization, education, and physical interpretation. For research use, the assumptions and limitations of each model should be considered together with the corresponding experimental conditions.