📐 Theory & Design of Metamaterials

Fundamental Theory

Electromagnetic Metamaterials Principle

Metamaterials derive their properties from their engineered structure rather than their chemical composition. By designing sub-wavelength resonant elements, we can achieve effective material parameters not found in nature.

Effective Medium Condition:

a ≪ λ → Homogenization valid
where a = unit cell size, λ = wavelength

Typically: a < λ/10

Negative Index Materials

Veselago's Prediction (1968):

If ε < 0 and μ < 0 simultaneously:
n = -√(εμ) < 0

Phase velocity: v_p = c/n (opposite to group velocity)
Poynting vector: S = E × H (opposite to k)

Key Consequences of Negative Index:

Reversed Snell's Law: n₁sinθ₁ = -n₂sinθ₂
Reversed Doppler Effect
Reversed Cherenkov Radiation
Amplification of evanescent waves

Design Methodology

Requirements

Define frequency, bandwidth, losses

Topology

Select SRR, fishnet, or hybrid

Simulation

FDTD/FEM optimization

Fabrication

Lithography & validation

Design Parameters Optimization

Parameter	Symbol	Typical Range	Effect on Response
Unit Cell Size	a	λ/10 - λ/4	Determines homogenization limit
Ring Radius	r	0.3a - 0.45a	Controls resonance frequency
Gap Width	g	0.05a - 0.2a	Affects capacitance & bandwidth
Metal Width	w	0.05a - 0.15a	Influences inductance & losses
Substrate Height	h	0.1a - 0.5a	Modifies coupling strength
Metal Thickness	t	30nm - 200nm	Determines ohmic losses

Design Rules for Optimal Performance

1 Resonance Frequency Tuning

f₀ = 1/(2π√(LC))
L ∝ μ₀πr²/w
C ∝ ε₀wt/g

2 Quality Factor Optimization

Q = f₀/Δf = ω₀L/R
Target: Q > 10 for sharp resonance
Trade-off: High Q → Narrow bandwidth

3 Loss Minimization

• Use high-conductivity metals (Au, Ag, Cu)
• Minimize surface roughness (< λ/100)
• Optimize metal thickness (> 3× skin depth)
• Select low-loss dielectrics (tan δ < 0.001)

4 Coupling Control

• Inter-unit coupling: Adjust periodicity
• Layer coupling: Control spacer thickness
• Near-field coupling: Optimize gap dimensions
• Far-field interaction: Consider array effects

Resonator Design Theory

Split-Ring Resonator (SRR) Analysis

LC Circuit Model:

Inductance: L = μ₀πr²(ln(8r/w) - 2)
Capacitance: C = ε₀εᵣ(πr/g)
Resistance: R = ρl/(wt)

Resonance: ω₀ = 1/√(LC)
Bandwidth: Δω = R/L

Fishnet Structure Analysis

Transmission Line Model:

Z_mesh = jωL_mesh/(1 - ω²LC_mesh)
Y_gap = jωC_gap + G_gap

Effective index: n_eff = √(Z_mesh × Y_gap)

Gap Size (nm)

FOM

Material Selection Guide

Gold (Au)

Conductivity 4.1×10⁷ S/m

Plasma freq 2.18×10¹⁵ Hz

Best for NIR-Visible

Silver (Ag)

Conductivity 6.3×10⁷ S/m

Plasma freq 2.18×10¹⁵ Hz

Best for Visible

Copper (Cu)

Conductivity 5.8×10⁷ S/m

Cost Low

Best for RF-THz

Silicon (Si)

Permittivity 11.7

Bandgap 1.12 eV

Best for Substrate

SiO₂

Permittivity 3.9

Loss tan 0.0001

Best for Spacer

Rogers RO4003

Permittivity 3.38

Loss tan 0.0027

Best for RF PCB

Advanced Design Concepts

Transformation Optics

Transformation optics allows the design of metamaterials with spatially varying properties to control electromagnetic wave propagation in unprecedented ways.

Coordinate Transformation:

ε'ᵢⱼ = (det J)⁻¹ Jᵢₖ εₖₗ Jⱼₗ
μ'ᵢⱼ = (det J)⁻¹ Jᵢₖ μₖₗ Jⱼₗ

where J = ∂x'/∂x (Jacobian matrix)

Metasurfaces

2D Metamaterials (Metasurfaces):

Phase Control: Gradient metasurfaces for beam steering
Polarization: Quarter/half-wave plates
Amplitude: Perfect absorbers
Holography: Computer-generated holograms

Topology Optimization

Objective Function:

min F(ρ) = ∫Ω f(E(ρ), H(ρ)) dΩ
subject to: ∇×E = iωμH
∇×H = -iωεE
∫Ω ρ dΩ ≤ V_max

Design Example: Invisible Cloak at 1.5 THz

Design Specifications:

Operating frequency: 1.5 THz (λ = 200 μm)
Cloak radius: 1 mm
Number of layers: 10
Material: Silver rings on silicon

Layer	Radius (μm)	εᵣ	μᵣ	Ring Size (nm)
1 (inner)	100	0.1	0.1	150
2	200	0.25	0.25	180
3	300	0.4	0.4	210
4	400	0.55	0.55	240
5	500	0.7	0.7	270

Design Performance Metrics

Critical Design Metrics:

Figure of Merit (FOM): |Re(n)|/Im(n) > 3
Bandwidth: Δf/f₀ > 0.1 for broadband
Insertion Loss: < 3 dB acceptable
Impedance Matching: |Z - Z₀|/Z₀ < 0.2
Fabrication Tolerance: ±10% dimension variation

Performance Equations:

FOM = |Re(n)|/Im(n)
Bandwidth = f₂ - f₁ where Re(n) < 0
Loss = -20 log₁₀|S₂₁| (dB)
Efficiency = |S₂₁|²/(1 - |S₁₁|²)