Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Chapter 1: Worked Examples and Learning Path

Question: what can you predict before opening an electromagnetic solver?

Read Chapter 1, make predictions in the browser lab, and then inspect or run the notebook. The notebook and browser calculate an ideal sheet response; neither assigns a physical geometry.

Learning path

StageActivityEvidence to retain
1Draw the structure and illuminationElement, pitch, period, aperture, and reference planes
2Interpret complex coefficientsAmplitude, phase, and appropriate power normalization
3Predict four sheet limitsA written prediction before using each preset
4Change detuning and strengthComplex response and conservation checks
5Examine frequency dependencePhase plus amplitude, and a stated spectral metric
6Compare evidence levelsAn analytical claim, a full-wave claim, and their limits

For a handmade demonstration, draw two complex arrows on movable cards. Arrange them head-to-tail to make zero, then align them to make a nonzero sum. Relate the construction to the model’s reflected and transmitted channel combinations. The cards illustrate complex addition; they do not measure the radiated fields of a physical meta-atom.

Example 1: wavelength in the relevant medium

At 200 THz200\,\mathrm{THz},

λ0=299.792458200 μm≈1.499 μm.\lambda_0=\frac{299.792458}{200}\,\mu\mathrm m\approx1.499\,\mu\mathrm m.

In a real-index medium with n=1.5n=1.5, the wavelength is approximately 0.999 μm0.999\,\mu\mathrm m. A dimension of 0.8 μm0.8\,\mu\mathrm m is therefore about 0.53λ00.53\lambda_0 but 0.80λ0.80\lambda in that medium. Always name the wavelength used in a “subwavelength” statement.

This conversion does not provide a homogenization theorem. Period, illumination, resonance, coupling, and the observable still determine whether an effective description is useful.

Example 2: the same power can have a different phase

For equal exterior impedances at normal incidence, t=1t=1 and t=−1t=-1 both give T=1T=1. Their phases differ by 180∘180^\circ at the same reference plane. Interfering either field with a fixed reference can distinguish them.

The coefficient t=0.8it=0.8i gives T=0.64T=0.64 and phase 90∘90^\circ in the same normalization. At t=0t=0, transmitted power is zero and phase is undefined. Reporting a phase there cannot establish a useful transmitted-wave design.

Prediction: if a solver reports only transmitted intensity, can it distinguish the first two coefficients? The intensity alone does not retain that relative phase.

Example 3: an electric-like resonance

At f=fef=f_e and nonzero electric-like strength, Ce=−1C_e=-1. With the magnetic-like channel disabled, Cm=1C_m=1, so

r=−1−12=−1,t=−1+12=0.r=\frac{-1-1}{2}=-1,\qquad t=\frac{-1+1}{2}=0.

The ideal sheet reflects all incident power. Its reflected field has a sign change, while its transmitted phase is undefined. The notebook stores the missing phase as null/blank rather than treating it as a physical zero-degree phase.

Example 4: a balanced pair at resonance

Now match the two channels, so Ce=Cm=−1C_e=C_m=-1:

r=0,t=−1,R=0,T=1.r=0,\qquad t=-1,\qquad R=0,\qquad T=1.

The two reflected contributions cancel, and the transmitted coefficient retains a 180∘180^\circ phase relative to the transparent reference. The model can therefore illustrate unit transmission with nonzero phase. It does not identify a ridge width, material, or fabricated sample that achieves the response.

Change the magnetic-like resonance to 230 THz230\,\mathrm{THz} while keeping illumination at 200 THz200\,\mathrm{THz}. Predict whether cancellation remains, then inspect the detuned preset. The model remains lossless, so the response redistributes power between its two ports.

Example 5: two-sided illumination is a stronger check

For S=[rttr]S=\begin{bmatrix}r&t\\t&r\end{bmatrix}, take the normalized even input (1,1)T/2(1,1)^T/\sqrt2. Its output is multiplied by r+t=Cer+t=C_e. The odd input (1,−1)T/2(1,-1)^T/\sqrt2 is multiplied by r−t=−Cmr-t=-C_m.

Both eigenvalues have unit magnitude. Thus this model preserves power for arbitrary combinations of the two incoming waves, not only the one-sided reference. The audit verifies the full unitary matrix and a general complex two-sided excitation.

The equality R+T=1R+T=1 for one incident direction would be insufficient by itself to establish every matrix constraint for a general proposed model.

Example 6: phase slope and delay

For the balanced model at resonance, the exact channel derivative gives

τg(f0)=1000×2πaf0[THz]  fs.\tau_g(f_0)=\frac{1000\times2}{\pi a f_0[\mathrm{THz}]}\;\mathrm{fs}.

At f0=200 THzf_0=200\,\mathrm{THz} and a=0.25a=0.25, the group delay is about 12.73 fs12.73\,\mathrm{fs}. Reducing aa makes the spectral phase transition sharper and increases its peak delay. The notebook compares a numerical derivative with the closed form and keeps amplitude beside phase.

The derivative concerns this imposed balanced analytical channel. It does not establish material dispersion, an achromatic geometry, or a geometry-specific bandwidth.

Unwrapped phase and exact group delay for two strengths of the matched analytical channel.

The executed notebook’s two balanced cases retain unit transmitted power by construction. Their different phase widths and delays concern this model, rather than a measured or simulated geometry bandwidth.

A bridge to saved full-wave data

For a worked comparison, inspect the first accepted row in the existing Chapter 3 TE ridge library. Its exterior indices are ni=1n_i=1 and nt=1.5n_t=1.5, with real nonmagnetic media, normal incidence, and only the specular order open. For the electric-field transmission coefficient,

T=ntni∣t∣2.T=\frac{n_t}{n_i}|t|^2.

This follows from a normally incident propagating wave’s average flux ∣E∣2/(2η)|E|^2/(2\eta) and η=η0/n\eta=\eta_0/n. The 0.08 μm ridge has ∣t∣2≈0.63996|t|^2\approx0.63996 but T≈0.95994T\approx0.95994. Its reflected power is about 0.04006. Adding R+∣t∣2R+|t|^2 would incorrectly report missing power. These values can be checked directly from the saved complex coefficient in that Chapter 3 dataset: t≈−0.797799−0.058944it\approx-0.797799-0.058944i.

The ridge has height 0.60 μm, pitch 0.65 μm, wavelength 1.55 μm, and ridge/gap indices 2.4/1. Its transmission reference is at the layer exit, whereas the introductory sheet uses its single sheet plane. This worked arithmetic uses a previously executed full-wave result; it does not rerun Maxwell’s equations or retrieve a geometry-to-sheet mapping. Oblique incidence and more general media need the appropriate admittance factors developed in Chapter 2.

Reading a research result

Choose a metasurface figure or abstract and record the illumination, structure/material model, target quantity, normalization, and evidence level. Distinguish analytical, simulated, and experimental results. Check whether “efficiency” refers to incident power, transmitted power, one polarization, or a finite collection region.

For physical context, read Decker et al. on dielectric Huygens surfaces or Arbabi et al. on phase and polarization control. Their structures and measured/simulated results are external references, not outputs of this introductory model.

Report checklist

Submit predictions, parameter changes, complex coefficients, power and phase plots, assumptions, checks, and a paragraph describing a claim the model cannot support. A successful report can conclude that more geometry/material information is needed. The lab provides the runnable source and a bridge to existing full-wave results.

Use the self-check and mini-project page for answer guidance, a suggested grading rubric, a glossary, and an editable report worksheet.