Terahertz and Infrared Probes

Terahertz and infrared probes measure how a material transmits, reflects, delays, rotates, and absorbs electromagnetic fields whose photon energies overlap many low-energy scales of quantum matter. Depending on frequency and sample, the response can contain mobile-carrier relaxation, superconducting gaps and stiffness, polar phonons, magnons, excitons, intersubband transitions, collective modes, and the onset of interband absorption.

The instrument does not directly record “the optical conductivity” or “the gap.” A Fourier-transform infrared spectrometer usually records detector power as a function of interferometer path difference. A terahertz time-domain spectrometer samples an electric-field waveform versus delay. Turning either record into a material response requires a reference measurement, optical boundary conditions, sample thickness and geometry, polarization calibration, and a causal model or inversion.

A useful evidence ladder is:

1. detector record: interferogram, reflected or transmitted power, or sampled electric field versus delay;
2. calibrated spectrum: frequency axis, instrument response, reference ratio, phase or polarization channel, and uncertainty;
3. optical coefficient: complex transmission, reflection, ellipsometric ratio, or surface impedance for a declared geometry;
4. material response: complex conductivity, dielectric function, refractive index, or sheet conductance from a Fresnel or multilayer forward model;
5. physical claim: carrier density and relaxation, energy gap, condensate stiffness, polar mode, polariton, or other excitation after alternatives and cross-probe checks.

Each step adds assumptions. A dip in transmission can be an intrinsic absorption, a substrate phonon, an etalon fringe, an aperture effect, or a reference mismatch. Trust begins by preserving those possibilities until the forward model separates them.

Canonical Scope

This page is the canonical home for terahertz and infrared measurement practice in quantum matter. It owns the relation among field or power records and complex electrodynamic functions; frequency-domain reflectance, transmission, and ellipsometry; field-resolved terahertz time-domain spectroscopy; thin-film and bulk inversion; and probe-specific evidence for Drude response, superconducting gaps and stiffness, polar phonons, and phonon polaritons.

Use Transport, Response, and Optics when the theoretical target has not yet been separated into conductivity, dielectric and loss, coherent, or nonlinear branches.

Drude Theory owns the relaxation-time derivation, Hall tensor, optical spectral weight, extended-Drude diagnostics, and failure modes of the model. Kubo Formula owns the causal current response, contact term, and order of limits. Sum Rules owns exact moment identities. This page uses those objects as outputs of an optical experiment rather than re-deriving them.

BCS Theory owns the pairing spectrum and coherence factors, while London Theory owns local superconducting stiffness and penetration. Phonons owns lattice eigenvectors and longitudinal–transverse mode splitting. Infrared Spectroscopy owns molecular dipole derivatives, rovibrational bands, ordinary FTIR, ATR, and fingerprint-region practice. The focus here is the electrodynamics of extended solids.

Time-domain THz spectroscopy of an equilibrium sample belongs here. An optical-pump–THz-probe experiment prepares a nonequilibrium state and therefore belongs conceptually with Pump–Probe Spectroscopy; this page supplies the equilibrium reference and the field-to-conductivity reduction needed for that analysis.

Spectral Coordinates and Energy Scales

Frequency labels overlap

Terahertz and infrared are laboratory regions, not distinct fundamental interactions. Boundaries depend on source, detector, community, and sample. One useful orientation is:

Region	Approximate span	Common quantum-matter content
microwave and millimetre wave	below roughly $0.1\ \mathrm{THz}$	narrow Drude response, resonances, stiffness, vortex dynamics
terahertz	roughly $0.1$ to $10\ \mathrm{THz}$	carrier relaxation, small gaps, magnons, soft modes, low polar phonons
far infrared	roughly $10$ to $400\ \mathrm{cm}^{-1}$	overlaps THz; polar phonons, gaps, collective modes
mid infrared	roughly $400$ to $4{,}000\ \mathrm{cm}^{-1}$	lattice and molecular vibrations, mid-infrared bands, low interband absorption

The overlap is deliberate. A $3\ \mathrm{THz}$ excitation may be described as terahertz, far infrared, $100\ \mathrm{cm}^{-1}$, or $12.4\ \mathrm{meV}$.

For ordinary frequency $\nu$, angular frequency $\omega$, vacuum wavenumber $\widetilde\nu$, and wavelength $\lambda_0$,

$$E = h\nu = \hbar\omega = hc\widetilde\nu, \qquad \lambda_0 = \frac{c}{\nu}.$$

Useful conversions are

$$\begin{aligned} 1\ \mathrm{THz} &= 4.13567\ \mathrm{meV}, \\ &= 33.3564\ \mathrm{cm}^{-1}, \\ &= 299.792\ \mathrm{\mu m} \quad \text{in vacuum}. \end{aligned}$$

The thermal frequency scale is

$$\frac{k_{\mathrm B}T}{h} \simeq 0.02084\ \mathrm{THz\,K^{-1}}\,T.$$

At $10\ \mathrm K$, this is about $0.208\ \mathrm{THz}$. Thermal occupation can therefore matter strongly for low-frequency modes even when it is negligible in visible spectroscopy.

Momentum and penetration

A freely propagating photon transfers little crystal momentum:

$$q_{\mathrm{ph}} \sim \frac{n\omega}{c}.$$

Terahertz and infrared far-field experiments usually probe the long-wavelength limit. They can identify zone-center dipole-active modes and uniform current response, but they do not map a full dispersion. Gratings, cavities, waveguides, attenuated-total-reflection geometries, and near-field tips can supply larger in-plane momentum.

Penetration is response dependent. An insulating crystal may transmit through a millimetre, a good metal may be sampled only over a skin depth, and a superconducting field may decay over a penetration depth. “Bulk” and “surface” are therefore frequency-, temperature-, and geometry-dependent statements.

From Optical Fields to Material Response

Convention ledger

Use time dependence $e^{-i\omega t}$ and write

$$\widetilde\sigma(\omega) = \sigma_1(\omega) + i\sigma_2(\omega).$$

$\sigma_1$ is the dissipative conductivity for a passive medium. $\sigma_2$ is the reactive part. With a background dielectric function $\varepsilon_\infty$,

$$\widetilde\varepsilon(\omega) = \varepsilon_1(\omega) + i\varepsilon_2(\omega) = \varepsilon_\infty + \frac{i\widetilde\sigma(\omega)} {\varepsilon_0\omega}.$$

Thus

$$\varepsilon_1 = \varepsilon_\infty - \frac{\sigma_2}{\varepsilon_0\omega}, \qquad \varepsilon_2 = \frac{\sigma_1}{\varepsilon_0\omega}.$$

For a nonmagnetic isotropic medium, define

$$\widetilde n = n+i\kappa, \qquad \widetilde n^2 = \widetilde\varepsilon.$$

The field varies as $e^{i\widetilde n\omega z/c-i\omega t}$, so $\kappa>0$ gives attenuation. The intensity absorption coefficient is

$$\alpha = \frac{2\omega\kappa}{c}.$$

These signs change if the Fourier convention changes. A reported $\sigma_2$ without the time convention is incomplete.

In an anisotropic crystal, $\widetilde\sigma$, $\widetilde\varepsilon$, and $\widetilde n$ are tensors. The measured projection depends on crystal axes, propagation direction, polarization, numerical aperture, and any domain average. A scalar inversion can fit an anisotropic sample while assigning the result to no actual tensor component.

Power measurements

At normal incidence from vacuum onto an optically thick isotropic half-space,

$$R(\omega) = \left| \frac{1-\widetilde n(\omega)} {1+\widetilde n(\omega)} \right|^2.$$

Reflectance is an intensity ratio and loses the complex reflection phase. Transmission through a finite slab additionally contains entrance and exit Fresnel factors, propagation phase, absorption, and usually internal reflections. For a smooth, nonscattering sample with all output channels collected,

$$R+T+A=1,$$

but this bookkeeping fails if diffuse scattering, finite numerical aperture, luminescence, diffraction, or an unmeasured polarization channel carries power away.

Fourier-transform infrared spectroscopy acquires an interferogram and reconstructs a power spectrum. A sample-to-reference ratio can give reflectance or transmittance, but not automatically an intrinsic absorption coefficient. The reference surface, purge, aperture, polarization, beam footprint, detector nonlinearity, and optical throughput must match the sample measurement.

Phase, causality, and model dependence

Real and imaginary response are linked by causality. For a causal susceptibility,

$$\chi_1(\omega) = \frac{2}{\pi} \mathcal P \int_0^\infty d\Omega\, \frac{\Omega\chi_2(\Omega)} {\Omega^2-\omega^2}.$$

A Kramers–Kronig transform of reflectance can recover phase only with a valid geometry and extrapolations outside the measured band. Conducting, superconducting, insulating, and interband regimes require different low- and high-frequency behavior. A fitted oscillator model can enforce causality, but its components are not unique merely because the total spectrum fits.

Ellipsometry measures the complex polarization ratio

$$\rho = \frac{r_p}{r_s} = \tan\Psi\,e^{i\Delta}.$$

It supplies relative phase information and is powerful for anisotropic or layered samples, but extracting dielectric tensors still requires orientation, roughness, thickness, and multilayer modeling.

Bulk, film, and substrate are different inverse problems

For a bulk slab, thickness uncertainty produces phase uncertainty, and internal reflections produce Fabry–Pérot fringes or delayed echoes. For a thin film on a substrate, the measured field traverses a coupled air–film–substrate stack. The natural low-frequency film observable is often the sheet conductance

$$\widetilde G(\omega) = d\,\widetilde\sigma(\omega),$$

where $d$ is film thickness.

In the thin-film limit, at normal incidence, for a film on a nonabsorbing substrate of index $n_s$, and after division by a matched bare-substrate reference,

$$\widetilde t_{\mathrm{rel}}(\omega) \simeq \frac{1+n_s} {1+n_s+Z_0\widetilde G(\omega)}.$$

Here $Z_0\simeq376.73\ \Omega$ is the vacuum impedance. This useful formula assumes that the film is electromagnetically thin, the film and reference substrates are optically matched, and unresolved multiple reflections can be neglected or removed consistently. Outside those conditions, use the full transfer matrix.

The measurement workflow. Power spectra require phase recovery or a causal forward model, while THz time-domain sampling retains field amplitude and phase. Free carriers, superconducting redistribution, and polar modes can overlap in frequency. A thin film is naturally reduced to sheet conductance only under declared boundary and thickness assumptions.

Low-Energy Electrodynamics

Drude response as a hypothesis

For one isotropic relaxation channel,

$$\widetilde\sigma_{\mathrm D}(\omega) = \frac{\varepsilon_0\omega_p^2} {\gamma-i\omega}, \qquad \gamma = \frac{1}{\tau}.$$

The full derivation and parameter meanings belong to Drude Theory. Experimentally, the key signatures are:

• $\sigma_1$ approaches a finite dc value and rolls over near $\omega\tau\sim1$;
• $\sigma_2/\sigma_1=\omega\tau$ for the one-component model;
• the integrated intraband weight is tied to $\omega_p^2$;
• the same parameters should be compatible with dc transport over the appropriate limit.

A narrow Drude peak can lie partly below the measured band. Fitting only its high-frequency tail makes $\omega_p$ and $\gamma$ strongly correlated. Conversely, a broad background can be an additional band, incoherent response, localization, a phonon tail, or an interband process. One Drude term is a model choice, not a definition of metallicity.

For several carrier channels,

$$\widetilde\sigma = \sum_a \frac{\varepsilon_0\omega_{p,a}^2} {\gamma_a-i\omega} + \widetilde\sigma_{\mathrm{reg}}.$$

The decomposition may be nonunique over a limited band. Temperature, field, doping, polarization, Hall data, and a wider spectral range can break parameter degeneracies.

Spectral weight and correlation

Define a partial optical weight

$$W(\Omega_c,T) = \int_0^{\Omega_c} d\omega\, \sigma_1(\omega,T).$$

Changing temperature or phase can move weight between a narrow Drude peak, a mid-infrared band, interband transitions, and a zero-frequency condensate. A partial weight depends on cutoff. Calling its variation a violation of a sum rule is unjustified unless omitted tails, tensor direction, units, and the corresponding exact operator identity are controlled.

An extended-Drude inversion can define frequency-dependent optical scattering and mass functions. These are compact representations of a complex response after choosing a plasma frequency. They are not uniquely measured microscopic self-energies, especially in multiband systems or where interband absorption overlaps the intraband response.

Gaps and thresholds outside superconductors

A suppression of $\sigma_1$ may indicate a band gap, density-wave gap, hybridization gap, localization scale, Coulomb gap, or superconducting gap. The threshold need not equal twice a single-particle gap:

• matrix elements can suppress the onset;
• indirect transitions need momentum assistance;
• excitons can appear below a continuum;
• disorder produces tails;
• coherence factors reshape thresholds;
• collective modes can lie inside a gap;
• a finite experimental floor can mimic zero absorption.

Optical conductivity is a two-particle current response. Compare it with tunneling or photoemission rather than treating all three as identical densities of states.

Superconducting Gap and Stiffness

What an ideal threshold means

In a local, isotropic, weak-coupling $s$-wave superconductor, a photon with negligible momentum can create two quasiparticles once

$$\hbar\omega \gtrsim 2\Delta$$

at low temperature. The Mattis–Bardeen dirty-limit response gives a sharp model for how $\sigma_1$ and $\sigma_2$ evolve across that scale. More general calculations include arbitrary purity, strong coupling, anisotropy, several bands, pair breaking, and nonlocality.

The threshold is kinematic; its optical weight depends on the current vertex, band structure, and momentum-relaxing processes. In a perfectly clean, translationally invariant single-band limit, regular absorption at negligible photon momentum can be strongly suppressed. The phrase “optical gap equals $2\Delta$” is therefore conditional. A nodal superconductor has low-energy quasiparticles rather than a hard threshold. Thermal quasiparticles create subgap absorption. Disorder can broaden or fill the gap, and a multiband material can display several weighted onsets. Phonons, substrate resonances, and Fabry–Pérot structure can overlap the same frequency range.

Reactive response and condensate weight

At low frequency, an ideal local superconducting contribution has

$$\sigma_{2,s}(\omega) \simeq \frac{D_s}{\omega}, \qquad \lambda^{-2} = \mu_0D_s.$$

$D_s$ is the electromagnetic stiffness in this convention and $\lambda$ is the local penetration depth. A $1/\omega$ rise in $\sigma_2$ is stronger evidence for coherent inductive response than a transmission increase alone, but regular quasiparticle and dielectric terms must be subtracted.

The finite-frequency weight lost from $\sigma_1$ enters the zero-frequency condensate. With a positive-frequency convention,

$$D_s = \frac{2}{\pi} \int_{0^+}^{\Omega_c} d\omega\, \left[ \sigma_{1,n}(\omega) - \sigma_{1,s}^{\mathrm{reg}}(\omega) \right],$$

provided $\Omega_c$ is high enough for the normal and superconducting weights to converge and the normal reference is valid. This is the operational Ferrell–Glover–Tinkham test. If the integral has not converged, the inferred stiffness remains cutoff dependent.

A defensible superconducting inference

Evidence should include:

1. complex $\widetilde\sigma$ or surface impedance, not only one transmission curve;
2. a documented normal-state reference and temperature calibration;
3. a low-frequency $\sigma_2$ consistent with stiffness and penetration measurements;
4. missing-area recovery within uncertainty;
5. a gap model appropriate to cleanliness, symmetry, band count, and scattering;
6. substrate, thickness, etalon, and phonon controls;
7. bulk superconductivity from magnetic or thermodynamic evidence.

An optical threshold constrains an electrodynamic excitation scale. It does not alone determine pairing symmetry or establish a bulk superconducting phase.

Polar Phonons and Polaritons

Infrared-active crystal modes

An infrared-active zone-center phonon carries a nonzero mode effective charge. Its coherent ionic displacement produces macroscopic polarization and couples to the electric field. A convenient dielectric model is

$$\widetilde\varepsilon(\omega) = \varepsilon_\infty + \sum_j \frac{ \Delta\varepsilon_j\omega_{\mathrm{TO},j}^2 }{ \omega_{\mathrm{TO},j}^2 - \omega^2 - i\gamma_j\omega }.$$

$\omega_{\mathrm{TO},j}$ is a transverse optical frequency in the model, $\gamma_j$ is damping, and $\Delta\varepsilon_j$ is dielectric strength. Coupled or low-symmetry modes may require factorized, tensor, or nonlocal forms; a sum of independent Lorentzians is not universal.

Transverse optical modes appear near poles of the transverse dielectric response. Longitudinal modes are associated with zeros of the relevant longitudinal dielectric function. For one ideal isotropic polar mode,

$$\frac{\varepsilon(0)} {\varepsilon_\infty} = \frac{\omega_{\mathrm{LO}}^2} {\omega_{\mathrm{TO}}^2}.$$

This Lyddane–Sachs–Teller relation generalizes for several modes and anisotropic crystals, where tensor determinants and mode mixing matter.

Between an ideal TO pole and LO zero, $\varepsilon_1$ can be negative and bulk propagation is strongly suppressed. The resulting high-reflectance interval is a Reststrahlen band. A Raman-active mode need not be infrared active: Raman coupling requires a polarizability derivative, while infrared coupling requires a dipole or mode-effective-charge channel.

Bulk and surface phonon polaritons

A transverse electromagnetic wave coupled to a polar optical phonon obeys, in an isotropic local medium,

$$q^2 = \widetilde\varepsilon(\omega) \frac{\omega^2}{c^2}.$$

The mixed eigenmodes are phonon polaritons. Their photon-like and phonon-like character varies along the dispersion. A far-field spectrum at small $q$ samples only a narrow part of that dispersion.

At a planar interface with a dielectric of permittivity $\varepsilon_d$, an electrostatic surface phonon-polariton condition is approximately

$$\varepsilon(\omega) + \varepsilon_d = 0.$$

The exact retarded dispersion includes both dielectric functions and the in-plane momentum. Near-field tips, prisms, gratings, resonators, or patterned structures provide the momentum needed to launch confined modes.

In anisotropic crystals, principal dielectric components can have opposite signs within a Reststrahlen band. The resulting hyperbolic phonon polaritons have strongly directional, high-momentum propagation. Assigning such a mode requires the dielectric tensor and boundary geometry, not just a reflectance peak.

Soft modes and phase transitions

A polar mode that softens can strongly enhance the static dielectric response. Frequency, damping, oscillator strength, central relaxation, domains, and conductivity should be fitted together. A low-frequency dielectric anomaly can come from electrodes, defects, domain walls, Maxwell–Wagner polarization, or mobile carriers rather than an intrinsic soft phonon.

Cross-check infrared mode symmetry with diffraction, Raman spectroscopy, isotope or pressure evolution, and calculated Born effective charges. A peak position alone rarely determines the atomic displacement pattern.

Time-Domain THz Spectroscopy

Sampling the electric field

A common THz time-domain spectrometer uses an ultrashort optical pulse to generate a broadband THz transient by a photoconductive antenna, optical rectification, or another coherent source. A synchronized gate samples the transmitted or reflected THz field with a photoconductive receiver or electro-optic crystal. Scanning delay produces

$$E_{\mathrm{ref}}(t), \qquad E_{\mathrm{sam}}(t).$$

Fourier transformation gives complex spectral fields:

$$\widetilde E(\omega) = \int_{-\infty}^{\infty} dt\, E(t)e^{i\omega t}.$$

The measured transfer function is

$$\widetilde T(\omega) = \frac{ \widetilde E_{\mathrm{sam}}(\omega) }{ \widetilde E_{\mathrm{ref}}(\omega) }.$$

Because the waveform is coherently sampled, both amplitude and phase survive. This often permits direct recovery of complex refractive index or conductivity over the usable bandwidth without a reflectance-only phase transform. It does not remove the need for a boundary model.

Thin-film inversion

For the thin-film expression above,

$$\widetilde G(\omega) \simeq \frac{1+n_s}{Z_0} \left[ \frac{1}{\widetilde t_{\mathrm{rel}}(\omega)} - 1 \right].$$

Dividing by an independently measured thickness gives bulk conductivity. If the active layer is atomically thin or thickness is ill-defined, sheet conductance is the less assumption-heavy result.

The simple inversion fails when:

• the film is not electromagnetically thin;
• substrate and reference thicknesses differ enough to create a phase ramp;
• substrate absorption or dispersion is neglected;
• internal reflections overlap the main pulse;
• the beam samples different apertures or positions;
• the film is laterally inhomogeneous or anisotropic;
• magnetic permeability or spatial dispersion matters.

A transfer-matrix fit to the complex field is preferable when several layers or echoes contribute.

Time window, bandwidth, and dynamic range

For a sampled window of duration $T_{\mathrm{win}}$, the natural Fourier grid is of order

$$\Delta\nu \sim \frac{1}{T_{\mathrm{win}}}.$$

Zero padding makes a smoother plotted curve but does not create finer physical resolution. Truncating the trace before a substrate echo shortens the usable window and broadens sharp features. Keeping the echo requires modeling its multiple-reflection physics.

The high-frequency limit is set by generation, detection, optics, sampling interval, and signal-to-noise. The low-frequency limit is affected by finite window, beam diffraction, apertures, and drift. Near frequencies where the reference field is small, division amplifies noise and can produce convincing but meaningless structure.

For a substrate of thickness $d_s$ and refractive index $n_s$, the round-trip echo delay is approximately

$$\Delta t_{\mathrm{echo}} \simeq \frac{2n_sd_s}{c}.$$

This time-domain separation is a strength of THz spectroscopy, but only if pulse duration, absorption, and window length permit distinct echoes.

Phase errors are material errors

A sample–reference displacement $\delta z$ produces an approximate phase error

$$\delta\phi(\omega) \simeq \frac{\omega\,\delta z}{c}$$

for a simple free-space offset. The resulting linear phase ramp biases refractive index and $\sigma_2$. Thickness uncertainty has a similar effect. Repeat insertion, fixed mounts, independent profilometry, and joint uncertainty propagation are therefore central, not administrative.

Atmospheric water has sharp rotational absorption throughout much of the THz band. Purge or vacuum, matched acquisition times, and inspection of known gas lines are standard controls. Numerical removal cannot recover information at a frequency where the sample signal has fallen below the noise floor.

Equilibrium versus driven response

Ordinary THz-TDS characterizes a stationary sample with a weak probe. In optical-pump–THz-probe spectroscopy, a pump changes the state and the measured quantity depends on pump–probe delay. Converting the differential field to a transient conductivity requires the excitation-depth profile, pump and probe penetration mismatch, nonstationarity during the THz waveform, and a multilayer nonequilibrium model.

Strong-field THz pulses can also drive nonlinear current, phonon, magnon, Josephson, or order-parameter dynamics. Once the response depends on field amplitude or waveform, a linear $\widetilde\sigma(\omega)$ is no longer a complete description.

Experimental Workflow

Before acquisition

1. Name the target response. Decide whether the claim requires sheet conductance, bulk conductivity, dielectric tensor, surface impedance, or only a calibrated optical coefficient.
2. Choose the frequency and temperature window. Include both sides of the expected relaxation, gap, or mode when possible.
3. Characterize the sample. Record thickness, roughness, area, crystallographic orientation, substrate, interfaces, contacts, and environment.
4. Write the forward stack. List every layer, index, thickness, polarization, incidence angle, and internal reflection retained in the model.
5. Plan references and controls. Use a matched substrate, mirror, open aperture, normal state, field-suppressed state, or empty cryostat as appropriate.

During acquisition

• monitor source power or field, detector linearity, purge, temperature, and drift;
• repeat sample and reference scans in an interleaved order;
• measure both polarization axes where anisotropy is possible;
• vary aperture, spot position, and sample orientation to expose diffraction or inhomogeneity;
• retain raw interferograms or time traces and acquisition metadata;
• measure thickness and substrate dispersion independently where feasible.

During reduction

1. transform sample and reference with identical windows and phase conventions;
2. mask frequencies below a declared dynamic-range threshold;
3. fit complex amplitude and phase together;
4. include substrate dispersion, echoes, roughness, and finite thickness at the level required by the claim;
5. propagate reference, thickness, alignment, baseline, and model uncertainty;
6. enforce causality and test sum-rule convergence over a wider band;
7. compare with dc, microwave, Raman, tunneling, thermodynamic, or scattering data.

Observable and Claim Ledger

Observed feature	Immediate inference	Additional requirement
low-frequency $\sigma_1$ rolloff and dispersive $\sigma_2$	finite current-relaxation scale	dc consistency and multichannel tests
increased subgap transmission below $T_c$	reduced dissipation in that band	complex conductivity, stiffness, and bulk superconductivity
$1/\omega$ contribution to $\sigma_2$	inductive spectral weight	regular-background subtraction and missing-area closure
absorption onset	allowed optical excitation scale	matrix element, disorder, and indirect-process model
Lorentz-like polar mode	dipole-active resonance	symmetry, tensor direction, and substrate separation
high reflectance between TO and LO scales	negative dielectric interval	causal dielectric fit and mode assignment
delayed pulse replica	internal optical path	thickness and index consistency
narrow atmospheric lines	gas-path mismatch	purge or matched reference, not a material fit

Common Mistakes

• Treating THz and infrared as sharply separated physical regimes.
• Calling transmittance, absorbance, loss function, $\sigma_1$, and $\varepsilon_2$ interchangeable.
• Reporting $\sigma_2$ without the Fourier-sign convention.
• Applying a half-space Fresnel formula to a thin film on a substrate.
• Dividing by film thickness when sheet conductance is the actually constrained quantity.
• Ignoring a substrate thickness mismatch because the power spectra look similar.
• Fitting only amplitude when phase is measured.
• Assuming a reflectance Kramers–Kronig inversion is independent of extrapolation.
• Calling every low-frequency peak a phonon or every suppression a gap.
• Equating an optical gap threshold with a unique single-particle gap.
• Using a one-component Drude fit without testing dc and wider-band consistency.
• Interpreting a fitted extended-Drude scattering rate as a unique self-energy.
• Identifying superconductivity from reduced absorption without stiffness or bulk evidence.
• Reading a surface polariton directly from a far-field peak without a momentum-coupling model.
• Zero-padding a short waveform and claiming improved resolution.
• Removing water lines numerically where the underlying sample field is below noise.
• Time-gating an echo without reporting the resulting spectral broadening.
• Treating an optical-pump–THz-probe differential spectrum as an equilibrium conductivity.

Exercises

1. Convert a terahertz scale

Convert $2.50\ \mathrm{THz}$ to meV, reciprocal centimetres, and vacuum wavelength.

Solution

Using the conversion ledger,

$$\begin{aligned} E &= (2.50)(4.13567)\ \mathrm{meV} \\ &= 10.34\ \mathrm{meV}, \\[4pt] \widetilde\nu &= (2.50)(33.3564)\ \mathrm{cm}^{-1} \\ &= 83.39\ \mathrm{cm}^{-1}, \\[4pt] \lambda_0 &= \frac{299.792\ \mathrm{\mu m}}{2.50} \\ &= 119.9\ \mathrm{\mu m}. \end{aligned}$$

This lies in the overlap commonly called terahertz and far infrared.

2. Invert a thin-film transmission

A $100\ \mathrm{nm}$ film lies on a substrate with $n_s=3.0$. At one frequency the complex film-plus-substrate transmission relative to a matched bare substrate is

$$\widetilde t_{\mathrm{rel}} = 0.80-0.10i.$$

Use the thin-film expression to estimate sheet conductance and bulk conductivity.

Solution

First,

$$\frac{1}{0.80-0.10i} = \frac{0.80+0.10i}{0.65} = 1.23077+0.15385i.$$

Therefore

$$\begin{aligned} \widetilde G &= \frac{4}{376.73\ \Omega} \left( 0.23077+0.15385i \right) \\ &\simeq \left( 2.45+1.63i \right) \times10^{-3}\ \mathrm S. \end{aligned}$$

With $d=1.00\times10^{-7}\ \mathrm m$,

$$\widetilde\sigma = \frac{\widetilde G}{d} \simeq \left( 2.45+1.63i \right) \times10^4\ \mathrm{S\,m^{-1}}.$$

The result inherits the thin-film, matched-substrate, normal-incidence, and time-convention assumptions.

3. Read a Drude point

A one-component Drude response has $\sigma_0=4.00\times10^4\ \mathrm{S\,m^{-1}}$ and $\tau=0.250\ \mathrm{ps}$. Find $\sigma_1$, $\sigma_2$, and their ratio at $\nu=0.800\ \mathrm{THz}$.

Solution

The dimensionless frequency is

$$\omega\tau = 2\pi (0.800\times10^{12}\ \mathrm{s^{-1}}) (0.250\times10^{-12}\ \mathrm s) = 1.257.$$

For

$$\widetilde\sigma = \frac{\sigma_0}{1-i\omega\tau},$$

one has

$$\begin{aligned} \sigma_1 &= \frac{\sigma_0}{1+(\omega\tau)^2} \simeq 1.55\times10^4\ \mathrm{S\,m^{-1}}, \\ \sigma_2 &= \frac{\sigma_0\omega\tau}{1+(\omega\tau)^2} \simeq 1.95\times10^4\ \mathrm{S\,m^{-1}}, \\ \frac{\sigma_2}{\sigma_1} &= \omega\tau \simeq 1.257. \end{aligned}$$

A measured ratio inconsistent with this value would reject the one-component parameter set even if one component alone looked plausible.

4. Infer a penetration depth

After subtracting regular dielectric and quasiparticle terms, a superconductor has $\sigma_2=3.00\times10^5\ \mathrm{S\,m^{-1}}$ at $\nu=0.500\ \mathrm{THz}$. Assuming $\sigma_2=D_s/\omega$, estimate $\lambda$.

Solution

The stiffness is

$$\begin{aligned} D_s &= \omega\sigma_2 \\ &= 2\pi (0.500\times10^{12}\ \mathrm{s^{-1}}) (3.00\times10^5\ \mathrm{S\,m^{-1}}) \\ &= 9.42\times10^{17} \ \mathrm{S\,m^{-1}\,s^{-1}}. \end{aligned}$$

Then

$$\begin{aligned} \lambda &= \frac{1}{\sqrt{\mu_0D_s}} \\ &\simeq 9.19\times10^{-7}\ \mathrm m \\ &= 0.919\ \mathrm{\mu m}. \end{aligned}$$

If regular terms were not removed, this procedure would misassign their reactive weight to the condensate.

5. Locate a polariton frequency

An ideal one-mode polar dielectric has the factorized response

$$\varepsilon(\nu) = 4\, \frac{3^2-\nu^2}{2^2-\nu^2},$$

where all frequencies are in THz, the TO pole is at $2\ \mathrm{THz}$, and the LO zero is at $3\ \mathrm{THz}$. Determine $\varepsilon(0)$, identify the negative-permittivity interval, and solve the vacuum-interface surface condition $\varepsilon=-1$.

Solution

At zero frequency,

$$\varepsilon(0) = 4\,\frac{9}{4} = 9.$$

For $2<\nu<3\ \mathrm{THz}$, the denominator is negative while the numerator is positive, so $\varepsilon<0$. This is the ideal Reststrahlen interval.

At a vacuum interface,

$$4\frac{9-\nu^2}{4-\nu^2} = -1,$$

so

$$36-4\nu^2 = -4+\nu^2, \qquad \nu = \sqrt{8}\ \mathrm{THz} \simeq 2.83\ \mathrm{THz}.$$

It lies between the $2\ \mathrm{THz}$ TO pole and $3\ \mathrm{THz}$ LO zero.

6. Echoes and resolution

A substrate is $0.500\ \mathrm{mm}$ thick with $n_s=3.40$. Estimate its round-trip echo delay. Compare the Fourier-grid scales for a trace gated just before the echo and for an $80\ \mathrm{ps}$ trace in which echoes are modeled.

Solution

The echo delay is

$$\begin{aligned} \Delta t_{\mathrm{echo}} &= \frac{2n_sd_s}{c} \\ &= \frac{ 2(3.40)(0.500\times10^{-3}\ \mathrm m) }{ 2.998\times10^8\ \mathrm{m\,s^{-1}} } \\ &\simeq 11.3\ \mathrm{ps}. \end{aligned}$$

A pre-echo gate has a frequency scale no finer than roughly

$$\Delta\nu_{\mathrm{gate}} \sim \frac{1}{11.3\ \mathrm{ps}} \simeq 0.088\ \mathrm{THz}.$$

An $80\ \mathrm{ps}$ modeled trace has

$$\Delta\nu_{80} \sim \frac{1}{80\ \mathrm{ps}} = 0.0125\ \mathrm{THz}.$$

The longer trace can resolve narrower structure only if the echo model, timing stability, and dynamic range are trustworthy.

7. Audit a superconducting-gap claim

Below a resistive transition, a thin film transmits more strongly below $0.8\ \mathrm{THz}$. A report identifies $2\Delta/h=0.8\ \mathrm{THz}$ from this observation alone. What is missing?

Solution

The transmission increase establishes reduced net attenuation in that sample–substrate stack. It does not yet isolate a superconducting pair-breaking threshold. A defensible analysis should add:

1. complex field or phase-sensitive data and a multilayer inversion to $\sigma_1$ and $\sigma_2$;
2. thickness, substrate, temperature, and reference uncertainties;
3. tests for phonons, etalons, apertures, and normal-state carrier narrowing;
4. a $1/\omega$ inductive response and missing-area consistency;
5. magnetic or thermodynamic confirmation of bulk superconductivity;
6. a Mattis–Bardeen or more appropriate gap model, including temperature, disorder, anisotropy, and multiple bands;
7. enough bandwidth to show both sides of the proposed onset.

The original observation is compatible with a gap, but it does not determine one uniquely.

Research Status

• Established: causal optical response, Fresnel and multilayer electrodynamics, Drude and Lorentz models in their stated domains, infrared-active phonons, Lyddane–Sachs–Teller relations, superconducting missing-area transfer, coherent THz field sampling, and transfer-matrix reduction.
• Model dependent: multicomponent Drude decompositions, extended-Drude self-energy language, gap extraction in anisotropic or multiband materials, oscillator decomposition of overlapping modes, effective-medium fits, and substrate-corrected thin-film conductivity.
• Active: nonlinear and multidimensional THz spectroscopy, extreme-condition THz access, nanoscale and on-chip THz probes, cavity-modified quantum matter, hyperbolic and nonlocal phonon polaritons, and uncertainty-aware inversion of heterogeneous multilayers.
• Unsupported without controls: assigning a new gap, condensate, exotic collective mode, or polariton from one intensity feature without complex response, boundary modeling, and complementary evidence.

Cross-Links

• How Quantum Matter Is Measured gives the general detector-to-response-to-claim contract.
• Drude Theory derives the free-carrier response and its responsible fitting limits.
• Kubo Formula derives causal electrical conductivity and the required contact term.
• Sum Rules owns exact spectral moments and cutoff audits.
• BCS Theory develops the quasiparticle gap and pairing evidence.
• London Theory owns electromagnetic stiffness and penetration depth.
• Superconducting Proximity Effect predicts layer-dependent gaps, stiffness transfer, and inverse suppression in superconducting multilayers; this page retains the Fresnel or transfer-matrix inversion, substrate, thickness, bandwidth, and background ledger.
• Vortex Matter, Pinning, and Flux Flow owns the pinning, creep, and driven-vortex interpretation of microwave and low-frequency mixed-state response; this page retains field-to-conductivity reduction, multilayer optics, backgrounds, and bandwidth limits.
• Phonons develops lattice eigenvectors, mode charges, and LO–TO splitting.
• Raman and Optical Spectroscopy gives complementary polarizability and symmetry channels.
• Infrared Spectroscopy owns molecular vibrational, FTIR, and ATR practice.
• Transport Measurements provides the dc and tensor limits against which optical transport is checked.

References

1. M. Dressel and G. Grüner, Electrodynamics of Solids: Optical Properties of Electrons in Matter (Cambridge University Press, 2002), doi:10.1017/CBO9780511606168.
2. F. Wooten, Optical Properties of Solids (Academic Press, 1972).
3. D. N. Basov, R. D. Averitt, D. van der Marel, M. Dressel, and K. Haule, “Electrodynamics of Correlated Electron Materials,” Reviews of Modern Physics 83, 471–541 (2011), doi:10.1103/RevModPhys.83.471.
4. A. B. Kuzmenko, “Kramers–Kronig Constrained Variational Analysis of Optical Spectra,” Review of Scientific Instruments 76, 083108 (2005), doi:10.1063/1.1979470.
5. R. E. Glover III and M. Tinkham, “Conductivity of Superconducting Films for Photon Energies between $0.3$ and $40k_{\mathrm B}T_c$,” Physical Review 108, 243–256 (1957), doi:10.1103/PhysRev.108.243.
6. R. A. Ferrell and R. E. Glover III, “Conductivity of Superconducting Films: A Sum Rule,” Physical Review 109, 1398–1399 (1958), doi:10.1103/PhysRev.109.1398.
7. M. Tinkham and R. A. Ferrell, “Determination of the Superconducting Skin Depth from the Energy Gap and Sum Rule,” Physical Review Letters 2, 331–333 (1959), doi:10.1103/PhysRevLett.2.331.
8. D. C. Mattis and J. Bardeen, “Theory of the Anomalous Skin Effect in Normal and Superconducting Metals,” Physical Review 111, 412–417 (1958), doi:10.1103/PhysRev.111.412.
9. W. Zimmermann, E. H. Brandt, M. Bauer, E. Seider, and L. Genzel, “Optical Conductivity of BCS Superconductors with Arbitrary Purity,” Physica C 183, 99–104 (1991), doi:10.1016/0921-4534(91)90771-P.
10. R. H. Lyddane, R. G. Sachs, and E. Teller, “On the Polar Vibrations of Alkali Halides,” Physical Review 59, 673–676 (1941), doi:10.1103/PhysRev.59.673.
11. K. Huang, “On the Interaction between the Radiation Field and Ionic Crystals,” Proceedings of the Royal Society A 208, 352–365 (1951), doi:10.1098/rspa.1951.0166.
12. J. D. Caldwell et al., “Low-Loss, Infrared and Terahertz Nanophotonics Using Surface Phonon Polaritons,” Nanophotonics 4, 44–68 (2015), doi:10.1515/nanoph-2014-0003.
13. D. N. Basov, M. M. Fogler, and F. J. García de Abajo, “Polaritons in van der Waals Materials,” Science 354, aag1992 (2016), doi:10.1126/science.aag1992.
14. D. H. Auston, K. P. Cheung, J. A. Valdmanis, and D. A. Kleinman, “Cherenkov Radiation from Femtosecond Optical Pulses in Electro-Optic Media,” Physical Review Letters 53, 1555–1558 (1984), doi:10.1103/PhysRevLett.53.1555.
15. D. Grischkowsky, S. Keiding, M. van Exter, and C. Fattinger, “Far-Infrared Time-Domain Spectroscopy with Terahertz Beams of Dielectrics and Semiconductors,” Journal of the Optical Society of America B 7, 2006–2015 (1990), doi:10.1364/JOSAB.7.002006.
16. Q. Wu and X.-C. Zhang, “Free-Space Electro-Optic Sampling of Terahertz Beams,” Applied Physics Letters 67, 3523–3525 (1995), doi:10.1063/1.114909.
17. P. U. Jepsen, D. G. Cooke, and M. Koch, “Terahertz Spectroscopy and Imaging: Modern Techniques and Applications,” Laser & Photonics Reviews 5, 124–166 (2011), doi:10.1002/lpor.201000011.
18. W. Withayachumnankul, B. M. Fischer, H. Lin, and D. Abbott, “Uncertainty in Terahertz Time-Domain Spectroscopy Measurement,” Journal of the Optical Society of America B 25, 1059–1072 (2008), doi:10.1364/JOSAB.25.001059.
19. M. van Exter, C. Fattinger, and D. Grischkowsky, “Terahertz Time-Domain Spectroscopy of Water Vapor,” Optics Letters 14, 1128–1130 (1989), doi:10.1364/OL.14.001128.
20. Joint Committee for Guides in Metrology, Evaluation of Measurement Data: Guide to the Expression of Uncertainty in Measurement, JCGM 100:2008.