core40 min studyLesson 19 in Physics

Dispersion and causal response

Connect a material's frequency-dependent polarization to delayed response, absorption, and the limits of pulse-velocity language.

On this page

What you’ll understand

  • Relate a causal time-response kernel to complex susceptibility.
  • Calculate a Lorentz-oscillator response with a stated Fourier convention.
  • Distinguish phase and group velocities from causal signal propagation.

A dielectric does not follow an applied electric field instantaneously at every frequency. Its charges have inertia, restoring forces, and damping, so polarization depends on earlier driving. The same delayed response produces both frequency-dependent phase and absorption. Modeling those effects consistently requires a complex susceptibility and an explicit time convention; otherwise an innocent sign error can turn a passive material into an apparent amplifier.

The material remembers earlier fields

In a linear, time-translation-invariant, spatially local medium, polarization is a convolution of the field with a response kernel. Causality requires the kernel to vanish for negative delay: an output cannot depend on a future input. Linearity is an approximation valid for sufficiently small fields, and spatial locality can fail when microscopic structure matters. These assumptions are separate from causality and should be stated rather than hidden inside one symbol.

Choose fields proportional to e−iωte^{-i\omega t}. Fourier transformation turns the convolution into multiplication by complex χ(ω)\chi(\omega). Its real part contributes to in-phase polarization, while its imaginary part describes the quadrature response and dissipation. For an ordinary passive medium at positive frequency, this convention gives nonnegative imaginary susceptibility. Reversing the time-exponential convention reverses the corresponding imaginary sign without changing the physical absorption.

P(t)=ϵ0∫0∞g(τ)E(t−τ) dτ,χ(ω)=∫0∞g(τ)eiωτ dτ\mathbf P(t)=\epsilon_0\int_0^\infty g(\tau)\mathbf E(t-\tau)\,d\tau,\qquad \chi(\omega)=\int_0^\infty g(\tau)e^{i\omega\tau}\,d\tau
The scalar kernel g has units s⁻¹ in this isotropic local model, making χ dimensionless. Tensor and spatially nonlocal responses generalize the formula. The lower integration limit enforces a retarded response.↔ Scroll wide equations horizontally.

A damped bound charge gives a concrete model

Model a bound charge as a driven oscillator with positive damping and a restoring frequency. Its displacement produces a dipole moment, and a density of such oscillators produces polarization. The steady-state response has a resonance rather than an instantaneous constant permittivity. Damping keeps the real-frequency response finite and transfers energy irreversibly to degrees of freedom omitted from the oscillator model.

Near resonance the imaginary part can be substantial even when the real susceptibility crosses through zero. A vanishing real part therefore does not mean the material is optically indistinguishable from vacuum. A material can contain multiple resonances, free-carrier terms, anisotropy, or nonlinearities; the one-resonance model is a useful interpretable fit only over a range where those omissions are acceptable.

χ(ω)=Fω02ω02−ω2−iγω,⟨pabs⟩=ωϵ02Im⁡χ(ω)∣E0∣2\chi(\omega)=\frac{F\omega_0^2}{\omega_0^2-\omega^2-i\gamma\omega},\qquad \langle p_{\mathrm{abs}}\rangle=\frac{\omega\epsilon_0}{2}\operatorname{Im}\chi(\omega)|E_0|^2
F is a dimensionless oscillator strength, γ > 0 is a damping rate in s⁻¹, and ω₀ is an angular resonance frequency. The second expression gives local average electric absorption in W/m³ for a peak field phasor E₀ and the e⁻ⁱωᵗ convention.↔ Scroll wide equations horizontally.

Causality links phase to absorption across frequency

A causal kernel makes susceptibility analytic in the appropriate complex-frequency half-plane. With sufficient high-frequency decay, its real and imaginary parts satisfy Kramers–Kronig relations. Changing absorption across a frequency band therefore constrains dispersion elsewhere; the two curves cannot be specified independently while retaining this causal linear model. Applying the relation to finite measurement data needs extrapolations, and those assumptions can dominate an inferred refractive index.

For a weakly absorbing narrow-band pulse, group velocity follows the derivative of its dispersion relation. Strong absorption or rapid dispersion can reshape a pulse so severely that the location of its maximum is a poor marker for transmitted information. A superluminal or negative formal group velocity in such a region does not establish faster-than-light signaling. The earliest causal front depends on the full broadband response, not one derivative at a carrier frequency.

Re⁡χ(ω)=2π P ⁣∫0∞ω′Im⁡χ(ω′)ω′2−ω2 dω′,vg=cn+ω dn/dω\operatorname{Re}\chi(\omega)=\frac{2}{\pi}\,\mathcal P\!\int_0^\infty\frac{\omega'\operatorname{Im}\chi(\omega')}{\omega'^2-\omega^2}\,d\omega',\qquad v_g=\frac{c}{n+\omega\,dn/d\omega}
The first formula assumes a real causal response and sufficiently decaying χ; a nonzero high-frequency limit requires subtraction. The second uses a real-index, weak-loss narrow-band approximation and must not be applied uncritically at an absorbing resonance.↔ Scroll wide equations horizontally.

Resonant absorption with a clear sign convention

A local Lorentz medium has F=0.40F=0.40, γ=0.10ω0\gamma=0.10\omega_0, and ω0=109 s−1\omega_0=10^9\,\mathrm{s^{-1}}. A local peak electric field of 10 V/m drives it at resonance. Calculate susceptibility and average absorbed power per volume.

  1. At ω=ω0\omega=\omega_0, the real denominator terms cancel. Thus χ=Fω02/(−iγω0)=iFω0/γ=4i\chi=F\omega_0^2/(-i\gamma\omega_0)=iF\omega_0/\gamma=4i. Its positive imaginary part is consistent with passive absorption for e−iωte^{-i\omega t}.

  2. The relative permittivity in this nonmagnetic electric model is ϵr=1+χ=1+4i\epsilon_r=1+\chi=1+4i. It is not equal to vacuum's permittivity, even though the real susceptibility is zero at resonance.

  3. Insert the local field into the loss formula: ⟨pabs⟩=(109)(8.85×10−12)(4)(102)/2=1.77 W/m3\langle p_{\mathrm{abs}}\rangle=(10^9)(8.85\times10^{-12})(4)(10^2)/2=1.77\,\mathrm{W/m^3}. The specified field is inside the material, so no separate interface-transmission factor belongs in this step.

The result

The susceptibility is 4i and the average electric absorption is 1.77 W/m³. Finding the internal field from an externally incident beam would require an additional boundary and propagation calculation.

Explore the model

Interactive lab

Interactive lab · oscillator

Damping in a related free oscillator

Adjust damping ratio ζ and natural frequency ω₀ for the analytic oscillator released from x(0)=1 and velocity(0)=0. It shows underdamped, critical, and overdamped relaxation. This is free motion rather than driven steady-state response: the lab does not sweep a driving frequency, calculate susceptibility, or evaluate Kramers–Kronig integrals.

Damping in a related free oscillatorUnit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Current readings: Displacement x(t): 1; Velocity v(t): 0; Regime: underdamped. All parameters are adjustable with the labeled controls below.02.557.510−1.15−0.5700.571.15Time · sDisplacementunderdamped
Displacement x(t)Current time
Displacement x(t)
1
Velocity v(t)
0
Regime
underdamped
t = 0 s

Unit mass, x(0) = 1 and v(0) = 0. The analytic solution includes underdamped, critically damped and overdamped cases; time is measured in seconds. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forA causal relaxation model supplies intuition, but a frequency response requires an explicit drive and Fourier analysis.

Check your understanding

With the e⁻ⁱωᵗ convention at positive frequency, what sign of Im χ describes passive electric absorption?

Try it yourself

For F = 0.40 and γ = 0.10ω₀, evaluate χ at ω = ω₀/2. Explain whether this frequency is more strongly absorbed than exact resonance for the same local electric amplitude.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Divide numerator and denominator by ω₀², then multiply by the complex-conjugate denominator.

Compare with the solution

The dimensionless denominator is 0.75−0.05i0.75-0.05i, so χ=0.4(0.75+0.05i)/(0.752+0.052)=0.531+0.0354i\chi=0.4(0.75+0.05i)/(0.75^2+0.05^2)=0.531+0.0354i. The absorption factor ωIm⁡χ\omega\operatorname{Im}\chi is 0.0177ω00.0177\omega_0, compared with 4ω04\omega_0 at resonance. It is therefore much smaller, about 0.00442 of the resonant absorption for the same internal peak amplitude. This comparison includes the changed frequency as well as the changed imaginary susceptibility.

Take it with you

  • Causality restricts which past fields can influence present polarization.
  • A Fourier convention fixes the sign that represents material absorption.
  • Dispersion, loss, and pulse reshaping must be interpreted together.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.