📚 Metamaterials Theory Handbook

Table of Contents

1. Fundamental Principles

1.1 Maxwell's Equations in Metamaterials

Metamaterials are artificial structures engineered to have properties not found in naturally occurring materials. The electromagnetic response is governed by Maxwell's equations:

Maxwell's Equations (Differential Form):
∇ · D = ρ
∇ · B = 0
∇ × E = -∂B/∂t
∇ × H = J + ∂D/∂t

Where the constitutive relations in metamaterials are:

Constitutive Relations:
D = ε₀εᵣ(ω)E
B = μ₀μᵣ(ω)H
J = σE
Key Insight: In metamaterials, both εᵣ(ω) and μᵣ(ω) can be simultaneously negative, leading to a negative refractive index n = -√(εᵣμᵣ).

1.2 Wave Propagation

The wave equation in a metamaterial medium:

Helmholtz Equation:
∇²E + k²E = 0
where k² = ω²εμ = (nω/c)²

For plane waves: E(r,t) = E₀ exp(i(k·r - ωt))

2. Effective Medium Theory

2.1 Homogenization

When the unit cell size a ≪ λ, the metamaterial can be treated as an effective medium with:

Effective Parameters:
εₑff = ⟨ε⟩ + corrections
μₑff = ⟨μ⟩ + corrections
Validity: a < λ/10

2.2 Lorentz Model for Magnetic Response

Lorentz Dispersion:
μ(ω) = 1 + Fω²/(ω₀² - ω² - iγω)

where:
F = filling factor
ω₀ = resonance frequency
γ = damping coefficient
Derivation:
Starting from the equation of motion for magnetic dipoles:
m̈ + γṁ + ω₀²m = αH₀e^(-iωt)
Solving in frequency domain yields the Lorentz form.

2.3 Drude Model for Electric Response

Drude Dispersion:
ε(ω) = 1 - ωₚ²/(ω² + iγω)

where:
ωₚ = plasma frequency
γ = collision frequency

3. Negative Index Materials

3.1 Conditions for Negative Index

Veselago's Conditions:
1. Re(ε) < 0 and Re(μ) < 0 simultaneously
2. Im(ε) > 0 and Im(μ) > 0 (passive medium)
3. Results in n = -√(εμ) (negative root chosen)

3.2 Reversed Physical Phenomena

Phenomenon Normal Material (n > 0) Metamaterial (n < 0)
Phase Velocity Same as group velocity Opposite to group velocity
Snell's Law n₁sinθ₁ = n₂sinθ₂ n₁sinθ₁ = -|n₂|sinθ₂
Doppler Shift Blue shift (approaching) Red shift (approaching)
Cherenkov Radiation Forward cone Backward cone

3.3 Perfect Lens

Pendry's Perfect Lens Conditions:
ε = μ = -1
Thickness: d
Resolution: Not limited by diffraction
Amplifies evanescent waves: exp(|kz|d)

4. Parameter Retrieval Methods

4.1 Nicolson-Ross-Weir (NRW) Method

S-Parameters to Material Parameters:
V₁ = S₂₁ + S₁₁
V₂ = S₂₁ - S₁₁

X = (1 - V₁V₂)/(V₁ - V₂)
Γ = X ± √(X² - 1)

T = (V₁ - Γ)/(1 - V₁Γ)

μ = (1 + Γ)/(Λ(1 - Γ))
ε = μ/Z²

where Λ = (1 - T²)/(2T), Z = √((1+S₁₁)²-S₂₁²)/√((1-S₁₁)²-S₂₁²)
Example Calculation:
Given: S₁₁ = 0.1 - 0.2i, S₂₁ = 0.8 - 0.3i, d = 200 μm, f = 1 THz
Result: n = -1.5 + 0.1i, Z = 0.8 - 0.05i
Extracted: ε = -2.3 + 0.2i, μ = -1.0 + 0.15i

4.2 Branch Selection and Phase Unwrapping

Critical Considerations:
1. Choose |Γ| < 1 for passive medium
2. Ensure Re(Z) > 0 for passive medium
3. Phase unwrapping: φₙ = φₙ₋₁ + Δφ (|Δφ| < π)
4. Kramers-Kronig consistency check

5. Transformation Optics

5.1 Coordinate Transformation

Transformation Relations:
x' = f(x, y, z)
y' = g(x, y, z)
z' = h(x, y, z)

Jacobian: Λ = ∂(x',y',z')/∂(x,y,z)

Transformed parameters:
ε' = ΛεΛᵀ/det(Λ)
μ' = ΛμΛᵀ/det(Λ)

5.2 Cylindrical Cloak Design

Cloak Parameters (2D):
r' = (b-a)r/b + a, for a < r < b

εᵣ = μᵣ = (r-a)/r
εθ = μθ = r/(r-a)
εz = μz = (b/(b-a))²(r-a)/r

where: a = inner radius, b = outer radius

6. Design Methodologies

6.1 Unit Cell Design Rules

  1. Size constraint: a < λ/10 for effective medium
  2. Symmetry: Use symmetry to reduce computational complexity
  3. LC resonance: f₀ = 1/(2π√LC)
  4. Quality factor: Q = ω₀/Δω determines bandwidth
  5. Bianisotropy: Minimize cross-coupling unless desired

6.2 Common Structures

Structure Response Frequency Range Key Parameters
Split-Ring Resonator Magnetic (μ < 0) MHz - THz Gap size, ring radius
Wire Array Electric (ε < 0) GHz - Optical Wire radius, spacing
Fishnet Both (n < 0) THz - Optical Hole size, periodicity
Cut-Wire Pairs Both (n < 0) GHz - THz Wire length, separation

6.3 Optimization Techniques

Figure of Merit (FOM):
FOM = |Re(n)|/Im(n)

Optimization objectives:
1. Maximize FOM > 3 for practical applications
2. Minimize losses: Im(ε), Im(μ) → 0
3. Maximize bandwidth: Δf/f₀ > 0.1
4. Impedance matching: Z → 1

7. Advanced Applications

7.1 Superlensing

Resolution Enhancement:
Conventional lens: Δx ≥ λ/(2NA)
Metamaterial superlens: Δx ~ λ/10

Transfer function:
T(kₓ) = exp(-2kz'd) for |kₓ| > k₀ (evanescent)
where kz' = √(k₀²ε'μ' - kₓ²)

7.2 Antenna Miniaturization

Size Reduction Factor:
L_meta/L_conv = 1/√|εᵣμᵣ|

Chu-Harrington Limit:
Q ≥ 1/(ka)³ + 1/(ka)
where k = 2π/λ, a = antenna radius

7.3 Perfect Absorption

Absorption Conditions:
A = 1 - |S₁₁|² - |S₂₁|²

For perfect absorption (A = 1):
1. Impedance matching: Z = Z₀
2. High losses: Im(n) >> 0
3. Zero transmission: |S₂₁| = 0
Design Example - THz Absorber:
Structure: Metal-Dielectric-Metal
Top layer: SRR array (200 nm Au)
Dielectric: 50 μm polyimide
Ground plane: 200 nm Au
Result: > 99% absorption at 1.5 THz

Key Takeaways