Electromagnetic Waves in Materials

Introduction

Electromagnetic waves interact with matter through electric polarization, magnetic polarization and charge transport. These processes determine phase velocity, wavelength, attenuation, impedance, power deposition, reflection and refraction. Material-wave interaction therefore governs dielectric heating, microwave circuits, optical systems, antennas, radar absorbers, biomedical sensing, electromagnetic compatibility and wireless propagation.

This chapter adopts sinusoidal steady state with the phasor convention e^jωt. A wave travelling in the positive z direction consequently contains the factor e^-γz, where γ = α + jβ. Unless otherwise stated, materials are linear, homogeneous and isotropic.

Learning Objectives

After studying this chapter, the reader should be able to:

• derive the wave equation, propagation constant and intrinsic impedance of a material;
• calculate exact attenuation and phase constants in a lossy medium;
• distinguish low-loss dielectrics, good conductors and general lossy materials;
• evaluate electromagnetic power flow and power dissipation correctly;
• apply boundary conditions to reflection, transmission and impedance-matching problems;
• interpret Debye and Lorentz dispersion models;
• analyse skin depth, polarization, anisotropy and engineered electromagnetic media;
• design and critically assess a quarter-wave matching layer.

1. Constitutive Description of Materials

Macroscopic electromagnetic behaviour is described through constitutive relations. For a linear, homogeneous, isotropic medium,

Here E is electric field intensity (V·m⁻¹), D is electric flux density (C·m⁻²), H is magnetic field intensity (A·m⁻¹), B is magnetic flux density (T), J_c is conduction-current density (A·m⁻²), ε is permittivity (F·m⁻¹), μ is permeability (H·m⁻¹), and σ is conductivity (S·m⁻¹).

The relative parameters are ε_r = ε/ε₀ and μ_r = μ/μ₀, where ε₀ ≈ 8.854 × 10^-12 F·m⁻¹ and μ₀ = 4π × 10^-7 H·m⁻¹.

In anisotropic media, permittivity and permeability are tensors:

The field direction and flux direction need not then coincide. Different polarization eigenmodes can experience different refractive indices, producing birefringence.

Complex permittivity and loss tangent

For the adopted phasor convention, dielectric loss is represented by

where ε' stores electric energy and ε'' ≥ 0 represents passive dielectric loss. Conduction can be included through an effective complex permittivity:

A commonly used total loss tangent is

This definition is frequency dependent and is meaningful only when the loss mechanisms and phasor convention are clearly stated.

2. Wave Equation in a Conducting Medium

In a source-free homogeneous region, Maxwell’s curl equations are

Taking the curl of Faraday’s law and using ∇ × (∇ × E)=∇(∇·E)-∇²E, with ∇·E=0 in a homogeneous source-free medium, gives

where

The square-root branch is selected so that α ≥ 0 for a passive medium. The magnetic field satisfies the same wave equation.

For a uniform plane wave travelling in the +z direction,

Thus α is the attenuation constant in Np·m⁻¹ and β is the phase constant in rad·m⁻¹. The wavelength and phase velocity are

Figure: A plane wave in a passive lossy medium; its field envelope decreases as e^-αz.

3. Exact Attenuation and Phase Constants

For real, positive, frequency-local values of μ, ε and σ, define

Separating the real and imaginary parts of γ gives the exact expressions

These formulas describe the complete range from low-loss dielectrics to good conductors. For general complex, dispersive ε(ω) and μ(ω), one must instead evaluate γ(ω) directly and select its passive branch.

3.1 Lossless dielectric

If σ=0 and ε and μ are real,

For a nonmagnetic dielectric, v_p ≈ c/√ε_r and λ ≈ λ₀/√ε_r, provided material dispersion is negligible.

3.2 Low-loss dielectric

When σ/(ωε) \ll 1, expansion of the exact expressions yields

If loss is represented by a small total loss tangent, an equivalent approximation is α ≈ (β\tan δ)/2.

3.3 Good conductor

When σ/(ωε) \gg 1,

The skin depth, defined as the distance over which field amplitude falls by e^-1, is

The power density falls as e^-2αz, so its e-folding distance is 1/(2α), not the field skin depth.

Worked Example 1: General lossy dielectric

A material has ε_r=4, μ_r=1, σ=0.020 S·m⁻¹ at 1 GHz. Determine α, β and wavelength.

First,

Using the exact expressions,

Hence

The electric-field attenuation is 8.686α ≈ 16.3 dB·m⁻¹. Because p is not extremely small, the exact calculation is preferable to an unqualified low-loss approximation.

4. Intrinsic Impedance and Field Relations

The intrinsic impedance of an isotropic medium is

For a forward uniform plane wave,

where \hat{k} is the propagation direction. In a lossless medium, η=√(μ/ε) is real and E and H are in phase. In a good conductor,

so the electric field leads the magnetic field by approximately 45° under the adopted convention.

5. Power Flow and Dissipation

The time-average Poynting vector is

For a forward plane wave with specified electric-field amplitude,

For a specified magnetic-field amplitude,

Thus it is not generally correct to say that \operatorname{Re}(η) alone determines power for an unspecified field amplitude: specified E requires \operatorname{Re}(1/η), whereas specified H requires \operatorname{Re}(η).

For a passive isotropic material with ε_c=ε'-jε'' and μ_c=μ'-jμ'', the average volumetric heating rate is

in W·m⁻³. This relation is central to induction heating, dielectric processing, specific absorption rate assessment and thermal multiphysics simulation.

6. Dispersion and Causal Material Models

A dispersive material has constitutive parameters that depend on frequency. Phase velocity and group velocity are then distinct:

In strongly absorbing media, group velocity may not represent signal or energy transport reliably; energy velocity or direct pulse analysis is preferable.

6.1 Debye relaxation

The single-pole Debye model describes orientational polarization without resonance:

Here ε_s is static permittivity, ε_∞ is high-frequency permittivity and τ is relaxation time. Separating components,

The loss component peaks at ωτ=1. Multiple-Debye models are frequently used for water-rich biological tissues and polar liquids.

6.2 Lorentz resonance

Bound-charge resonance is represented by the Lorentz model:

Here ω_0n is a resonance angular frequency, Γ_n is a damping coefficient, and Δε_n is the oscillator strength contribution. Lorentz models are important in infrared and optical materials, ferrites and engineered resonant structures.

Dispersion and absorption cannot be assigned independently in a causal passive material. Their real and imaginary response functions are connected by the Kramers–Kronig relations. Consequently, constant permittivity extrapolated over an unlimited bandwidth is generally nonphysical.

7. Reflection and Transmission at Interfaces

At an interface without imposed surface charge or surface current, tangential E and H are continuous. For normal incidence from medium 1 to medium 2, the electric-field reflection and transmission coefficients are

where T_E=E_t/E_i. For lossless media with real positive impedances,

These simple power formulas require lossless incident and transmitted media. In lossy media, power must be calculated from the real Poynting vector.

At oblique incidence, phase matching gives Snell’s law,

for isotropic lossless media with real refractive indices. TE and TM polarizations have different wave impedances and Fresnel coefficients. For nonmagnetic lossless dielectrics, TM reflection vanishes at the Brewster angle \tan θ_B=n₂/n₁. Total internal reflection occurs when a wave travels from higher to lower refractive index and θ_i>\sin^-1(n₂/n₁).

8. Quarter-Wave Impedance Matching

Consider lossless media 0 and 2 separated by a lossless matching layer 1 at normal incidence. The input impedance of the layer is

If d=λ₁/4, then β₁d=π/2 and

Zero reflection at the design frequency requires η_in=η₀, giving

This condition matches two generally different surrounding media; it does not require them to be equal. The result assumes normal incidence, lossless homogeneous media and a single design frequency.

Figure: A quarter-wave layer matching two generally different lossless media at normal incidence.

9. Polarization, Anisotropy and Engineered Media

Polarization specifies the locus traced by the electric-field vector. Linear, circular and elliptical polarization result from the relative amplitudes and phases of two orthogonal field components. An anisotropic material can split an incident wave into eigenmodes with different phase constants, producing polarization rotation or retardation.

Ferrites biased by a static magnetic field possess tensor permeability and support nonreciprocal propagation. Applications include isolators, circulators and phase shifters. Plasmas and metals may exhibit negative real permittivity over selected frequency ranges. Metamaterials use subwavelength structures to synthesize effective constitutive responses, including near-zero or simultaneously negative effective parameters. Such descriptions must respect homogenization limits, passivity, causality, spatial dispersion and unit-cell size.

Contemporary metasurfaces achieve thin-profile control of reflection, transmission and wavefront phase. Tunable implementations employ varactors, liquid crystals, graphene, ferroelectrics or semiconductor switches. AI-assisted inverse design can optimize geometry, but its solutions still require full-wave verification, tolerance analysis and physically causal material models.

10. Material Characterization

Material parameters are usually inferred rather than measured directly. Common methods include:

Reliable characterization requires temperature control, uncertainty reporting, surface preparation and a model compatible with the sample’s anisotropy and dispersion. Fitting independent permittivity values at each frequency may reproduce data but violate causality; Debye or Lorentz parameterizations are often more defensible.

11. Applied Case: Microwave Radome Matching Layer

Background

An illustrative engineering team must reduce normal-incidence reflection from air into a nonmagnetic radome substrate having ε_r2=4 at 10 GHz. Loss is initially neglected to establish a first-order design.

Problem

Air has η₀≈377 Ω, while the substrate impedance is

Without a coating,

so the reflected power fraction is R=|Γ|²≈0.111, or 11.1%.

Approach and solution

The required matching-layer impedance is

For a nonmagnetic coating, η₁=377/√ε_r1, giving

At 10 GHz, the free-space wavelength is approximately 30.0 mm. The wavelength in the coating is 30.0/√2≈21.2 mm, so

Outcome and critical evaluation

The ideal model predicts zero reflection at 10 GHz and normal incidence. A production design must additionally account for coating loss, oblique incidence, TE/TM polarization, curved geometry, adhesive layers, moisture, temperature dependence and manufacturing tolerances. A multilayer or graded-index coating may be chosen when bandwidth is more important than minimum reflection at one frequency.

Learning points

• Matching is between two different media through the geometric-mean impedance.
• Electrical thickness, not merely physical thickness, determines quarter-wave behaviour.
• A mathematically perfect single-frequency design may be operationally inadequate.
• Material uncertainty should be propagated into predicted reflection and bandwidth.

12. Engineering Practice, Safety and Governance

Electromagnetic material design is inseparable from thermal and safety considerations. Absorbed RF energy can cause tissue heating, dielectric breakdown, arcing or thermal runaway. Applicable exposure limits, specific absorption rate requirements, electromagnetic compatibility rules and product safety standards must be identified for the jurisdiction and frequency range concerned.

Material databases and AI-generated constitutive models should not be accepted without provenance. Best practice includes documenting measurement methods, calibration, uncertainty, frequency and temperature ranges, tensor orientation, fitted-model assumptions and extrapolation limits. Full-wave simulation should be validated through analytical limiting cases and representative measurements.

Chapter Summary

Electromagnetic propagation in materials follows from Maxwell’s equations and constitutive relations. The complex propagation constant γ=α+jβ determines attenuation and phase, while intrinsic impedance relates electric and magnetic fields. Exact formulas for α and β cover general conducting dielectrics; low-loss and good-conductor equations are controlled approximations.

Average power depends on the specified field: \operatorname{Re}(1/η) applies for specified electric field, whereas \operatorname{Re}(η) applies for specified magnetic field. Interfaces produce reflection and transmission governed by impedance and boundary conditions. A quarter-wave layer matches two generally different lossless media when η₁=√(η₀η₂). Debye and Lorentz models connect microscopic polarization mechanisms with dispersion and absorption. Advanced material systems—including anisotropic media, ferrites and metamaterials—require attention to causality, passivity, homogenization and measurement uncertainty.

Key Terms

Attenuation constant: Real part of the propagation constant, measured in Np·m⁻¹.

Birefringence: Different phase responses experienced by polarization eigenmodes in an anisotropic medium.

Complex permittivity: Frequency-domain permittivity containing energy-storage and loss components.

Debye model: Relaxation model for frequency-dependent orientational polarization.

Dispersion: Dependence of electromagnetic material response or propagation on frequency.

Group velocity: Reciprocal of dβ/dω in a sufficiently weakly attenuating propagation regime.

Intrinsic impedance: Complex ratio relating electric and magnetic fields in a uniform plane wave.

Lorentz model: Damped-resonance model for bound-charge polarization.

Loss tangent: Ratio of effective loss response to electric energy-storage response.

Metamaterial: Structured medium engineered to exhibit an effective electromagnetic response.

Phase constant: Imaginary part of the propagation constant, measured in rad·m⁻¹.

Phase velocity: Rate of advance of a constant-phase surface, ω/β.

Poynting vector: Vector representing electromagnetic power-flow density.

Propagation constant: Complex quantity γ=α+jβ governing amplitude and phase variation.

Quarter-wave transformer: Layer of one-quarter wavelength used for impedance matching at a design frequency.

Skin depth: Distance over which field amplitude decreases by a factor of e^-1 in a conductor.

Spatial dispersion: Dependence of constitutive response on wave vector as well as frequency.

Wave impedance: Ratio of transverse electric and magnetic fields for a specified propagation mode.

Review Questions

Short Answer Questions

1. Distinguish attenuation constant from phase constant.
2. State the assumptions underlying the exact formulas for α and β presented in this chapter.
3. Why are electric and magnetic fields not generally in phase in a lossy medium?
4. What physical process is represented by the Debye relaxation time?
5. Define skin depth and distinguish field attenuation from power attenuation.

Long Answer Questions

1. Derive the electromagnetic wave equation and propagation constant for a homogeneous conducting medium.
2. Starting from the general propagation constant, derive the good-conductor approximations for α, β, η and skin depth.
3. Explain reflection and transmission at normal and oblique incidence, including the roles of polarization and impedance.
4. Compare Debye and Lorentz material models with respect to physical mechanism, mathematical form and applications.
5. Discuss the opportunities and limitations of anisotropic materials, ferrites, metamaterials and metasurfaces.

Analytical / Application Questions

1. A dielectric has ε_r=6, μ_r=1, σ=0.04 S·m⁻¹ at 2 GHz. Calculate its exact α, β, wavelength and electric-field attenuation in dB·m⁻¹.
2. Design a nonmagnetic quarter-wave coating that matches air to a lossless substrate with ε_r=9 at 12 GHz. Determine coating permittivity and thickness.
3. Develop an experimental plan to estimate a material’s complex permittivity over a broad microwave band, including calibration, causal model fitting and uncertainty assessment.

Multiple Choice Questions with Answers

Question 1

For the phasor convention e^jωt, which factor represents a wave travelling in the positive z-direction in a passive medium?

A. e^+αze^+jβz

B. e^-αze^-jβz

C. e^-αze^+jβz

D. e^+αze^-jβz

Correct Answer: B. e^-αze^-jβz

Question 2

For real positive μ, ε and σ, the exact attenuation constant is:

A. ω√{(με/2)(√{1+(σ/ωε)²}-1)}

B. ω√{(με/2)(√{1+(σ/ωε)²}+1)}

C. σ√{με}

D. ωμε/σ

Correct Answer: A. ω√{(με/2)(√{1+(σ/ωε)²}-1)}

Question 3

For a forward plane wave with specified electric-field amplitude in a lossy medium, average power density depends on:

A. Re(η) only

B. Im(η) only

C. Re(1/η)

D. |η|²

Correct Answer: C. Re(1/η)

Question 4

A lossless nonmagnetic dielectric has ε_r=4. At 3 GHz, its wavelength is approximately:

A. 0.025 m

B. 0.050 m

C. 0.100 m

D. 0.200 m

Correct Answer: B. 0.050 m

Question 5

A quarter-wave layer matches two lossless media of impedances η₀ and η₂ at normal incidence when its impedance is:

A. (η₀+η₂)/2

B. η₀η₂

C. √(η₀η₂)

D. |η₀-η₂|

Correct Answer: C. √(η₀η₂)

Question 6

In the single-pole Debye model, dielectric loss reaches its maximum when:

A. ωτ=0

B. ωτ=1

C. ωτ=π

D. ωτ→∞

Correct Answer: B. ωτ=1

Question 7

For a good conductor, which approximation is correct?

A. α≈0

B. α≈β≈√(πfμσ)

C. β≈σ/(ωε)

D. α≈ω√(με)

Correct Answer: B. α≈β≈√(πfμσ)

Question 8

A normally incident wave travels from a 377 Ω lossless medium into a 188.5 Ω lossless medium. Its electric-field reflection coefficient is approximately:

A. +0.333

B. −0.333

C. +0.500

D. −0.500

Correct Answer: B. −0.333

Question 9

Which model explicitly represents damped bound-charge resonance?

A. Debye model

B. Lorentz model

C. Ohm’s law

D. Snell’s law

Correct Answer: B. Lorentz model

Question 10

If the attenuation constant is α=2 Np·m⁻¹, the electric-field amplitude after travelling 0.5 m is what fraction of its initial value?

A. e^-0.5

B. e^-1

C. e^-2

D. 1/2

Correct Answer: B. e^-1