Simulates the extraction of the refractive index of a glass prism using variable angle reflectometry. In an ideal textbook scenario, the bulk index is directly extracted from the Brewster angle, the precise coordinate where p-polarized reflection reaches an absolute intensity null. However, experimental metrology is constrained by physical realities. This interactive simulation explores how nanometer-scale surface contamination films, finite hardware extinction ratios, azimuthal polarization leakage, and mechanical goniometer offsets modify the ideal Fresnel scenario. By modulating these parameters, users can observe the physical mechanisms that mask the intensity minimum, distort ellipsometric phase transitions, and systematically skew the analytically extracted index.
The phase thickness (β) of the surface contamination film incorporates the wavelength (λ) and the optical path difference:
The complex 3-layer reflection coefficient (r) for each polarization state (p, s) is calculated by superimposing the Fresnel boundaries (0=air, 1=film, 2=bulk) with the film's phase term:
Hardware polarization leakage mixes the absolute intensity states. The measured reflectances incorporating the extinction ratio (ε) and analyzer misalignment (γ) are:
The ellipsometric parameters—the amplitude ratio (Ψ) and the phase difference (Δ)—are traditionally defined by the complex reflectance ratio (ρ):
When hardware leakage corrupts the intensity, the extracted Ψ is derived directly from the measured macroscopic reflectances. Conversely, Δ remains governed by the fundamental phase arguments (δ) of the complex interfacial coefficients:
Standard educational derivations model the Brewster angle as a two-layer system (air-to-bulk), predicting a mathematically absolute zero in the p-polarized reflectance (Rp → 0) accompanied by an abrupt 180° to 0° phase transition in Δ. However, real optical glass possesses a nanometer-scale hydration or oxidation overlayer. To accurately model this physical reality, we evaluate the complex 3-layer Airy formulation. As the film thickness (d) increases, the single boundary breaks down. The absolute zero lifts off the axis, forming a shallow pseudo-Brewster minimum, and the perfectly vertical Δ jump visibly relaxes into a smooth sigmoid. Crucially, the angular position of the minimum shifts, resulting in a systemic miscalculation of the bulk refractive index if analyzed using the ideal tangent equation.
The metrological search for the Brewster angle fundamentally relies on detecting an intensity minimum. However, two distinct hardware limitations mask this null:
Both errors physically truncate the bottom of the logarithmic intensity plot. The sharp mathematical turning point is replaced by a broad, flat valley. At the Brewster angle, the ideal p-polarized reflectance approaches zero (|rp|2 → 0), and the measured intensity is dominated by the orthogonal s-polarized state leaking through the optics:
Applying the small-angle approximation (sin γ ≈ γ and cos γ ≈ 1) reduces this response to:
Because γ2 and ε act as mathematically identical additive constants to the signal floor, they create nearly perfect parameter covariance (degeneracy). An inverse regression algorithm cannot independently distinguish a poor extinction ratio from a mechanical misalignment based solely on the intensity minimum. Note that while these errors degrade the intensity null, they do not alter the underlying interfacial phase shift (the Δ curve remains unaffected).
If the rotation stage is improperly zeroed against the normal reference plane (θoff), the entire optical response function is rigidly translated across the angular domain. Unlike multi-point curve-fitting methods (which often flag systematic offsets by analyzing residual bowing), the standard Brewster extraction is a single-point operation utilizing the tangent equation:
Therefore, a zero-offset error bypasses statistical suppression and propagates linearly into a skewed refractive index measurement.
Algorithmic Note: The observed null angle (θi,obs) is extracted via a discrete computational grid search executing 2000 evaluations over the restricted domain [50°, 70°] with a resolution of Δθ = 0.01°. The corresponding uncorrected bulk index is then analytically extracted using the explicitly defined tangent formulation.