Maxwell equations in matter

In macroscopic matter with free charge density $\rho$ and free current density $J$,The electric displacement $D$ and magnetic field strength $H$ encode the material response.

These are the Maxwell equations after separating free and bound material sources.

Constitutive relations specify how matter responds to fields. A homogeneous isotropic linear electromagnetic medium obeys $\mathbf{D}=\epsilon \mathbf{E}$ and $\mathbf{B}=\mu \mathbf{H}$.

Electromagnetic permeability relates magnetic flux density to magnetic field strength through $\mathbf{B}=\mu \mathbf{H}$ in an isotropic linear medium.

Vacuum permeability is the magnetic permeability appearing in the vacuum Maxwell equations. In SI, it relates vacuum magnetic field strength to magnetic flux density through $\mathbf{B}={\mu }_0\mathbf{H}$ and obeys ${\mu }_0{\epsilon }_0{c}^2=1$, where $c$ is the speed of light and ${\epsilon }_0$ the vacuum permittivity.

A linear anisotropic medium has constitutive relations $D_i={\epsilon }_{ij}{E}_j$ and $B_i={\mu }_{ij}{H}_j$. Symmetric, time-independent tensors make the field-energy density a quadratic form.

Here permittivity is a tensor, so the electric displacement need not be parallel to the electric field.

Birefringence is the existence of different optical propagation constants for different eigenpolarizations in an anisotropic material. It provides the propagation-mode separation needed for adiabatic optical following in a twisted nematic.

Two different refractive indices can produce two displaced images, as in the calcite double-image experiment, or a phase retardation between overlapping components, as in a waveplate. Spatial separation concerns where energy flux travels; phase retardation concerns the relative wave phase. A specimen need not show two visibly separated images to have birefringence.

A known retarder inserted in a polarization microscope to compare phase delay with an unknown specimen. Parallel slow directions add signed retardations; a slow direction paired with a fast one subtracts them. General relative angles require composing polarization transformations rather than scalar addition.

A quartz retarder used to identify specimen fast and slow directions through changes of interference colour or compensation. Its own retardation and orientation must be specified before unique final colours can be predicted.

For a lossless uniaxial dielectric, $n_{\mathrm{ex}}^{-2}={\cos }^2\vartheta /n_o^2+{\sin }^2\vartheta /n_e^2$, where $\vartheta$ is the angle between the wave vector and optic axis. Along the axis $n_{\mathrm{ex}}=n_o$, while perpendicular propagation gives $n_{\mathrm{ex}}=n_e$.

Colour produced when wavelength-dependent interference selectively transmits or reflects part of white light. Between crossed polarizers, a uniform birefringent plate gives a factor ${\sin }^2(\pi \Gamma /\lambda )$, where $\Gamma$ is its optical path difference. The observed colour also depends on illumination and the spectral response of the observer.

A birefringent plate with retardation $\Gamma =t\Delta n$ gives $I/I_P={\sin }^2(2\alpha ){\sin }^2(\pi \Gamma /\lambda )$. At non-extinction azimuths, wavelengths satisfying $\Gamma =m\lambda$ are removed by a crossed analyser. At a principal-axis azimuth every wavelength is extinguished.

Birefringence produced by subwavelength structure rather than solely by intrinsic crystalline anisotropy. For ideal dielectric lamellae, continuity of tangential electric field gives ${ϵ}_{\parallel }=f{ϵ}_1+(1-f){ϵ}_2$, while continuity of normal electric displacement field gives ${ϵ}_{\perp }^{-1}=f/{ϵ}_1+(1-f)/{ϵ}_2$. The unequal effective permittivities give different principal refractive indices.

Linear birefringence means that two orthogonal linear polarizations have different refractive indices. Their different propagation phases produce phase retardation. In a lossless waveplate, with fast-axis angle $\theta$ and slow-minus-fast retardation $\delta$, the Jones matrix, up to a common phase, isusing time dependence $e^{-i\omega t}$. Circular birefringence instead has circular eigenpolarizations.

Place a clear calcite specimen on a transparency bearing printed words on an overhead projector. The projected readout is two copies of the lettering. Put a rotatable linear polarizer after the specimen: suitable orientations suppress one copy or the other. These are the actual apparatus and observations in the University of Iowa demonstration.^[ref]

For a useful control, remove the specimen while keeping the transparency, projector and screen fixed: the doubling should disappear. Replace it, hold it stationary and rotate only the polarizer; complementary brightness changes distinguish two linear polarizations from an accidental duplicate transparency. These are suggested controls, rather than additional reported measurements. Ideally the two copies obey complementary Malus's law curves, $I_o(\alpha )=I_{o,\max }{\cos }^2\alpha$ and $I_e(\alpha )=I_{e,\max }{\sin }^2\alpha$, with angle $\alpha$ measured from the exiting ordinary wave's linear polarization. Unequal maxima do not invalidate the test. Finite polarizer extinction, poor focus and specimen defects limit visibility; image separation alone does not determine either refractive index without the specimen geometry.

A transparent, nonmagnetic linear anisotropic dielectric is uniaxial when two principal permittivities agree. With the optic axis along $z$, its permittivity tensor has the displayed diagonal matrix form. Here $n_o$ and $n_e$ are principal refractive indices at a specified frequency; $n_e$ is not the extraordinary wave's index in every direction.^[ref]

The symmetry direction of a uniaxial optical medium. For a wavevector along it, the two refractive indices coincide: there is no birefringence for that propagation direction, even though $n_e\ne n_o$.

The ordinary wave's electric field is perpendicular to the plane containing the wavevector and optic axis. Its refractive index $n_o$ is direction independent at fixed frequency.

The extraordinary wave's electric displacement field lies in that plane and is perpendicular to the wavevector; its electric field need not be. For angle $\theta$ between the wavevector and optic axis,Thus $n_{\mathrm{ex}}(0)=n_o$ and $n_{\mathrm{ex}}(\pi /2)=n_e$.^[ref]

A short derivation starts from the existing plane wave in an anisotropic dielectric equation. Put $\mathbf{k}=(k_x,0,k_z)$ and $k_0=\omega /c$, with angular frequency $\omega$ and vacuum speed of light $c$. The $y$ component gives the ordinary wave relation $k_x^2+k_z^2=n_o^2{k}_0^2$. The determinant of the $x,z$ system givesSubstituting $k_x=n_{\mathrm{ex}}{k}_0\sin \theta$ and $k_z=n_{\mathrm{ex}}{k}_0\cos \theta$ yields the displayed directional refractive index.

The extraordinary wave's Poynting vector can differ in direction from its wavevector. The latter specifies wave phase propagation, while the former specifies energy flux. This angular separation is walk-off.^[ref]

In the lossless, nondispersive model above, take the gradient of $\omega =c\sqrt{k_x^2/n_e^2+k_z^2/n_o^2}$ to obtain the group velocity, parallel to the Poynting vector. If its angle to the optic axis is $\varphi$, this gives $\tan \varphi =(n_o^2/n_e^2)\tan \theta$. Consequently the extraordinary wave generally carries energy away from its wavevector direction. An internal wavevector normal to a slab face therefore does not guarantee a normal energy path. Birefringent walk-off vanishes at $\theta =0$ and $\theta =\pi /2$; the latter direction still has unequal refractive indices. This is a useful waveplate geometry: overlapping components can acquire phase retardation without spatial splitting. A thin specimen, insufficient imaging resolution, or excitation of only one linear polarization can also hide a double image, for different reasons.

For a lossless nonmagnetic uniaxial optical medium, the electric displacement field of a plane electromagnetic wave is perpendicular to its wavevector. The magnetic field is perpendicular to the plane containing the electric field and wavevector, while the Poynting vector is perpendicular to the electric field in that same plane. Thus the angle between $\mathbf{E}$ and $\mathbf{D}$ equals the walk-off angle between energy flow and phase propagation. At normal incidence on a parallel slab of thickness $L$, a lateral separation $\Delta$ of the ordinary and extraordinary exiting rays gives ${\alpha }_e=\arctan (|\Delta |/L)$, while ${\alpha }_o=0$. The output rays can be parallel even though their internal energy paths differ.

A waveplate uses birefringence to change polarization of an electromagnetic wave through phase retardation. Its fast axis has the smaller refractive index, hence the larger phase velocity; its slow axis has the larger index. Consider normal incidence with the optic axis in the plate face, so both components follow the same straight path. A linear polarization at $45^{\circ }$ to these axes excites equal wave amplitudes.

A complementary experiment uses two sheet polarizers with a cellophane sheet between them. First cross the polarizers and remove the sheet: the field is dark. Insert and rotate the sheet about the beam: light returns at intermediate orientations, but extinction recurs at two perpendicular orientations. The sheet has changed the polarization of an electromagnetic wave, without requiring a resolved double image. With white light, different wavelengths can produce different transmitted colors.^[ref]. Brightness alone does not establish circular polarization: an ideal rotating linear polarizer also gives orientation-independent light intensity for unpolarized light.

A detected lower-path photon with vertical polarization evolves under an angle-dependent Jones matrix. For fast-axis angle $\eta$, its horizontal amplitude is $(1-i)\sin \eta \cos \eta$ and its vertical amplitude is ${\sin }^2\eta +i{\cos }^2\eta$. A later vertical-path probability therefore depends on the plate orientation and on whether the earlier lower-path event has been conditioned on.

The relative phase accumulated by orthogonal eigenpolarizations is $\Delta =2\pi L(n_2-n_1)/{\lambda }_0$. A quarter-wave plate has relative phase $\pi /2$ modulo $2\pi$, while a half-wave plate has $\pi$ modulo $2\pi$.

For the overlapping monochromatic wave components, let $d$ be the thickness and ${\lambda }_0$ the vacuum wavelength. Define $\Delta n=n_s-n_f>0$ using the slow and fast refractive indices. With complex amplitude convention $e^{i(kz-\omega t)}$, the accumulated slow-minus-fast phase difference isThis is the optical path difference $\Delta nd$ expressed as a wave phase. It is a signed accumulation convention; with the opposite time convention the sign reverses. The outgoing polarization of an electromagnetic wave depends on $\delta$ modulo $2\pi$. Here the positive slow-minus-fast $\Delta n$ differs from the signed principal-index $n_e-n_o$ used in the parent birefringence entry when $n_e<n_o$.

For an original worked example, stipulate $\Delta n=0.0100$ at ${\lambda }_0=600\mathrm{nm}$ and a uniform thickness $d=15.0\mu \mathrm{m}$. These are design assumptions, not measured material data. The optical path difference is $150\mathrm{nm}$, soIt is a quarter-wave plate. At the same wavelength and refractive indices, doubling $d$ to $30.0\mu \mathrm{m}$ gives $\delta =\pi$, a half-wave plate. A $1.0\mu \mathrm{m}$ thickness increase changes $\delta$ by $\pi /30$, or $6^{\circ }$; retaining the same birefringence at $500\mathrm{nm}$ would instead give $\delta =3\pi /5$ for the $15.0\mu \mathrm{m}$ plate. Real refractive indices also vary with wavelength, so this last calculation explicitly neglects optical dispersion.

For these equal transmitted wave amplitudes, orthogonal projection onto a polarizer crossed with the initial $45^{\circ }$ linear polarization gives a simple independent readout. After removing a common wave phase, the component complex amplitudes are proportional to $(1,e^{i\delta })/\sqrt{2}$; orthogonal projection onto $(1,-1)/\sqrt{2}$ gives $(1-e^{i\delta })/2$. Squaring its modulus yields $I_{\perp }/I_{\mathrm{in}}={\sin }^2(\delta /2)$: zero at $d=0$, one half at $15\mu \mathrm{m}$, and one at $30\mu \mathrm{m}$. This assumes ideal polarizers, negligible specular reflection and wave absorption, overlapping components and stable relative wave phase. It measures a periodic function of phase retardation, so one brightness measurement cannot determine the thickness uniquely.

A zero-order quarter-wave plate satisfies ${\omega }_{0}d\Delta n({\omega }_0)/c=\pi /2$. Its retardance error at nearby angular frequency is $ϵ=\omega d\Delta n(\omega )/c-\pi /2$. Taylor expansion gives the title expression, including material dispersion. Neglecting that dispersion yields $ϵ=(\pi /2)\Delta \omega /{\omega }_0$. A half-wave plate made from twice the same material thickness has exactly twice this phase error. Common propagation phase can be removed from its Jones matrix without changing polarization.

A quarter-wave plate gives a quarter-cycle phase retardation. Equal transmitted wave amplitudes initially in linear polarization at $45^{\circ }$ become circular polarization; unequal wave amplitudes generally give elliptical polarization. Reversing the sign reverses the traversal of the same circle.

Take a half-wave plate with fast axis at $15^{\circ }$ to a linearly polarized input, followed by a quarter-wave plate at a further $60^{\circ }$, so the second axis is at $75^{\circ }$ to the input. With shared quarter-wave retardance error $ϵ$ and half-wave error $2ϵ$, use $Q=\mathrm{diag}(1,i{e}^{iϵ})$ and $H=\mathrm{diag}(1,-e^{2iϵ})$. Rotating each matrix into the input axes and multiplying gives $\mathbf{J}(0)=e^{i75^{\circ }}(1,-i{)}^{\mathsf{T}}/\sqrt{2}$ and the title expansion. Its first-order change is common phase only, so the output remains circular to first order. This input-specific local compensation does not prove a frequency-independent quarter-wave operator for arbitrary input across a finite spectral band.

Suppose the unknown ideal elements are a linear polarizer and a quarter-wave plate, and the two unknown sources are unpolarized and circularly polarized. Test both element orders while rotating their relative axes and monitoring output intensity. Only circular input followed by the quarter-wave plate and then the polarizer gives complete modulation: the plate makes linear light and the analyzer obeys Malus's law. Polarizer-first output is constant half-intensity because the following lossless plate changes polarization but not intensity. Unpolarized input is constant in either order. This identifies both elements and both sources without comparing source brightnesses.

Two rotatable ideal quarter-wave plates can transform a fixed circular polarization into any pure polarization represented by a Jones vector. Use time dependence $e^{-i\omega t}$, $Q(\theta )=R(-\theta )\mathrm{diag}(1,i)R(\theta )$ and input $J_L=(1,i{)}^T/\sqrt{2}$. Every target polarization ellipse can be written, up to a common phase, as $R(-\beta )(\cos \eta ,i\sin \eta {)}^T$, with signed ellipticity angle $-\pi /4\le \eta \le \pi /4$. The first quarter-wave plate at ${\theta }_1=\beta +\eta +\pi /4$ produces linear polarization at $\beta +\eta$. The second, at ${\theta }_2=\beta$, therefore produces the required target. The two angles supply the two independent parameters of pure polarization.

For a reciprocal ideal quarter-wave plate followed by a mirror at normal incidence, use the same fixed transverse laboratory axes for the forward and returning Jones vectors. If the axes coincide with the plate's fast and slow axes, each pass has Jones matrix $J=\mathrm{diag}(1,i)$ up to a common phase factor. A polarization-independent mirror has scalar reflection coefficient $r_m$, so the returned Jones vector is $r_m{J}^2\mathbf{e}=r_m\mathrm{diag}(1,-1)\mathbf{e}$. The plate acts twice as a half-wave plate, rather than undoing its phase retardation.

For input linear polarization at angle $\alpha$ to the fast axis, the returned line is at $-\alpha$ modulo $\pi$. At $\alpha =45^{\circ }$ it is perpendicular to the original line. An axis-aligned input instead returns along the same line. General elliptical polarization requires its full input Jones vector; no unique returned state follows from the word “polarized” alone. The propagation direction reverses, so comparing descriptions of circular polarization must also specify the viewing convention.

Take a linear polarizer $P=\mathrm{diag}(1,0)$ and a reciprocal quarter-wave plate $Q=2^{-1/2}\left(\begin{array}{cc}1& i\\ i& 1\end{array}\right)$ in fixed physical transverse axes. A near-normal isotropic mirror has $M=-I_2$ up to global phase. If the source/detector side encounters P before Q, the round-trip matrix is $PQMQP=-P{Q}^{2}P=0$. Reversing the filter gives $QPMPQ=-QPQ$, which transmits half an initially unpolarized intensity. The result is double-pass rejection, not optical nonreciprocity; changing direction without changing the input polarization does not preserve the same physical input state.

A half-wave plate reverses the relative sign of the two components. For initial linear polarization at angle $\alpha$ to the fast axis, the outgoing line is at $-\alpha$ modulo $\pi$. The $45^{\circ }$ input therefore rotates by $90^{\circ }$.

Circular birefringence means different refractive indices for opposite circular polarizations. In a transparent isotropic medium, their relative wave phase changes rotate the direction of transmitted linear polarization; optical activity gives the solution measurement and convention.

The rotatory power is the rotation of transmitted linear polarization per unit path length, in $\mathrm{rad}{\mathrm{m}}^{-1}$. With equally transmitted circular eigenmodes $({\mathbf{e}}_1\mp i{\mathbf{e}}_2)/\sqrt{2}$ and time convention $e^{-i\omega t}$, their accumulated phase difference gives $\theta =({\varphi }_{+}-{\varphi }_{-})/2$. For uniform circular birefringence, $\rho =\pi (n_{+}-n_{-})/{\lambda }_0$. Generic linear birefringence need not admit one rotation angle common to all incident polarization directions. In a weak longitudinal Faraday effect, $\rho =V{B}_{\parallel }$, with signed Verdet constant $V$.