Prism: TIR Critical-Angle Goniometry (Grazing Emergence)

Evaluates the closed-form analytical solution for extracting the refractive index of a prism using the critical angle (grazing emergence) method. In standard transmission goniometry, the extracted index depends symmetrically on the incident and exit angles. By targeting the grazing emergence condition where the exit beam grazes the prism face, the extraction becomes insensitive to first-order alignment errors in the exit angle. This simulation calculates the impact of systematic angular offsets on the global dispersion curve (Top Row) and plots the localized error sensitivity across the physical emergence domain (Bottom Row), evaluating the technique's error propagation characteristics.

Diagnostic Controls

Global Dispersion Errors
Simulates systematic calibration errors applied uniformly across all 5 discrete wavelength measurements to observe the resulting shift in the extracted dispersion curve.
Physical miscalibration of the prism's internal geometry.
Systematic mechanical misalignment of the incident normal reference.
Miscalibration or mechanical error in the detector arm reading.
Local Sensitivity Baseline
Defines the specific physical parameters used to evaluate the localized error derivatives (bottom panels). The angular perturbation slider (Δ) defines the magnitude of the simulated mechanical alignment error.
Selects the reference refractive index for the optical baseline.
True bulk refractive index. Overridden if a specific wavelength is selected.
Assumed physical internal geometry of the prism.
The targeted emergence angle (cursor for the continuous plot below).
The magnitude of the simulated mechanical alignment error (e.g., misreading the goniometer by Δ°) used to evaluate the resulting extraction error (Δn).
Governing Analytical Model
nglass = (nair / sin A) · [ sin2 θi + sin2 θexit + 2 sin θi sin θexit cos A ]1/2

Global Dispersion Shift

Sellmeier SF10 Curve vs. Offset Experimental Data

Numerical Evaluation (Fixed Dataset)

Evaluated at Base Geometry: A = 60.096°, θexit = 90°

Contrasts the theoretical Sellmeier index against the index extracted from the nominal geometry (Nominal nmeas) and the error-injected geometry (Perturbed nmeas).

λ (nm) Sellmeier nglass Nominal nmeas Perturbed nmeas Error (Δn)

Local Sensitivity Profile

Maps the absolute extraction error (Δn) induced by an isolated mechanical miscalibration. The magnitude of this simulated error is dictated by the Angular Perturbation (Δ) slider. The horizontal bar lengths compare the error sensitivity of each angular parameter at the currently selected target geometry.

Apex Angle (A)
0.00000
Incident Angle (θi)
0.00000
Exit Angle (θexit)
0.00000

Sensitivity vs. Target Emergence Geometry

Continuous Error Evaluation Domain: [60°, 90°]

Demonstrates that the extraction sensitivity to emergence angle miscalibration (∂n / ∂θexit) converges strictly to zero at the 90° grazing limit, while sensitivities to the incident angle and apex geometry remain non-zero.

Error Mitigation via the Critical Angle

1. Sellmeier Dispersion Baseline

The standard refractive index of SF10 glass across the visible spectrum is defined by the Sellmeier equation:

nglass = [ 1 + (B1 λ2) / (λ2 − C1) + (B2 λ2) / (λ2 − C2) + (B3 λ2) / (λ2 − C3) ]1/2

where λ is the wavelength in micrometers, and Bj and Cj are experimentally derived material coefficients.

2. The General Refractive Index Extraction Equation

By determining the exterior incident angle (θi) required to produce a specific exterior emergence angle (θexit), the refractive index of a prism with an apex angle A can be extracted using the general formulation:

nglass = (nair / sin A) · [ sin2 θi + sin2 θexit + 2 sin θi sin θexit cos A ]1/2

where nair is the refractive index of the ambient environment (modeled here as standard air, nair = 1.0003). The domain of valid physical solutions requires the sum within the brackets to be non-negative. A negative value yields an imaginary index, which corresponds to a non-physical propagation regime where the beam undergoes Total Internal Reflection (TIR) and cannot emerge from the prism.

3. Zero First-Order Sensitivity at Grazing Emergence

As the emergence angle approaches the 90° grazing limit, the system nears the TIR threshold. Physically, the Fresnel transmission coefficients approach zero, reducing the emergent beam intensity. Concurrently, the transverse cross-section of the beam compresses, causing diffraction and spatial broadening of the spectral line observed in the telescope. These factors introduce increased mechanical and visual uncertainty when defining the exact angular cut-off.

However, the critical angle method mitigates this physical uncertainty because the extraction function is stationary at the limit. Applying the chain rule to the extraction equation, taking the partial derivative of the terms inside the radical with respect to θexit yields:

∂ (sin2 θexit) / ∂θexit = 2 sin θexit cos θexit

∂ (2 sin θi sin θexit cos A) / ∂θexit = 2 sin θi cos θexit cos A

Because both derivative components share a factor of cos θexit, the partial derivative evaluates to zero at grazing emergence (cos 90° = 0). The overall sensitivity term ∂n / ∂θexit converges strictly to zero, meaning the extracted refractive index is insensitive to first-order alignment errors in the detector arm.

4. Error Propagation in Associated Parameters

The dynamic bar chart and continuous plot demonstrate that this zero first-order sensitivity does not extend to the incident angle (θi) or the prism's physical apex angle (A). The partial derivatives ∂n / ∂θi and ∂n / ∂A remain non-zero at the critical threshold. Therefore, while the experimenter can tolerate higher uncertainty in judging the emergence disappearance, the goniometer's incident arm reading and the initial calibration of the prism's fixed geometry require precise calibration to prevent systematic errors.

5. Computational Implementation

While the analytical derivative (∂n / ∂θexit = 0) provides the theoretical proof, the JavaScript kernel evaluates local sensitivity using a discrete finite-difference method. It computes the absolute index shift (Δn) resulting from a finite angular perturbation (Δ). Physically, this perturbation (Δ) represents the magnitude of the experimenter's mechanical uncertainty (e.g., misreading the goniometer's nonius scale by 0.5°). To render the continuous sensitivity profile, the algorithm iteratively cascades Snell's law backwards to find the required incident angle for each target emergence coordinate across the domain [60°, 90°] using a resolution of Δθ = 0.5°. If a calculation requires a sine argument exceeding ±1.0, it indicates TIR at the internal boundary; the algorithm halts the evaluation for that coordinate and flags the state as invalid to prevent non-physical NaN propagation.

Authored by E. U. Donev. HTML/JS generated in collaboration with Gemini 3.1 Pro (Extended) AI.