Theory

Optical Concepts

Refractive index, polarization, angle, birefringence, and ellipsometry

This chapter defines the optical quantities that appear most often in this app and maps them to the relevant fields and result pages.

Scope of this chapter

ConceptFocus
Complex refractive indexRoles of n and k, and app-side representation
Polarization and angleEffect of s / p and pRatio on reflection and transmission
Cone angle averagingInput modes, sampling, distribution weights, and scope
BirefringenceMeaning of nExt / kExt and validity range
Ellipsometric quantitiesPhysical meaning of Psi / Delta and validity range

Complex refractive index and input models

The most basic material quantity is the complex refractive index:

N = n + i k

where:

  • n controls phase velocity and refraction;
  • k controls absorption and is typically required to satisfy k >= 0.

The current app supports three input models.

Input modelMain fieldsTypical useKey limit
Constantn, kMaterials that can be approximated as wavelength-independent over the relevant rangeIsotropic only
Constant Birefringencen, k, nExt, kExtSimplified birefringent materials with separate ordinary and extraordinary constantsNot available for incoherent layers
FileSampled wavelength-dependent dataMeasured dispersion, tabulated materials, or imported database dataIf the file is birefringent, it is not available for incoherent layers

For engineering interpretation, prioritize these three rules:

  1. Increasing n changes both optical path length and interface contrast, so peak position and peak height may both move.
  2. Increasing k increases absorption, which usually raises A and lowers transmission.
  3. If dispersion is important, prefer File; otherwise fringe locations and color results can drift from the physical material.

Polarization, pRatio, and incident angle

At oblique incidence, the electric field must be defined relative to the plane of incidence.

QuantityDefinitionApp representationEngineering meaning
s polarizationElectric field perpendicular to the plane of incidencepRatio = 0Often shows a different reflectance response from p
p polarizationElectric field parallel to the plane of incidencepRatio = 1Can show markedly reduced reflection near specific angles
Mixed polarizationLinear combination of s and p0 < pRatio < 1Useful for partially polarized light
Unpolarized approximationEqual weighting of s and ppRatio = 0.5Good default for first-pass evaluation

Incident Angle is the external incidence angle. In the solver, it changes all of the following:

  1. the propagation angle inside each layer;
  2. the Fresnel coefficients at each interface;
  3. the effective optical path length through each layer.

As a result, changing the angle usually affects both spectral position and contrast, and it strongly affects Psi / Delta. For most introductory cases, start from or a moderate oblique angle such as 60° ~ 70°.

Using Snell's law to interpret refraction through the stack

Use the following relation to interpret refraction through the stack:

n_i sin(theta_i) = n_j sin(theta_j)

In the context of this app, that means:

ObservationPhysical reasonPractical use when reading results
Spectral features shift as angle increasesOblique incidence changes the effective optical pathExplains angle-sweep peak movement
Propagation angles are smaller in higher-index layersRays bend toward the normal in higher-n mediaHelps interpret phase-thickness changes across the stack
s and p curves separateFresnel coefficients respond differently for the two polarizationsExplains polarization-sensitive spectra and ellipsometry

The app does not ask you to enter internal layer angles manually. The important point is that changing Incident Angle changes the propagation condition of the full stack, not just the first interface.

Cone angle averaging

Standard TMM assumes ideal plane-wave incidence, but real optical systems use convergent or divergent beams with a finite angular spread. Cone angle averaging models this by computing results for multiple rays distributed over a cone centered on the nominal incident angle and returning the weighted average.

Input modes and conversion formulas

The software accepts three ways to specify the cone size. All are converted internally to the half-cone angle θ_cone.

Input modeUser-entered quantityConversion formula
Half-angleθ_cone (degrees)Used directly
F-number (F/#)F/#θ_cone = arctan(1 / (2 × F/#))
Numerical Aperture (NA)NAθ_cone = arcsin(NA / n_medium)

Here n_medium is the refractive index of the incidence medium. When the incidence medium is air (n ≈ 1), NA is equivalent to sin(θ_cone).

Ring sampling

Discrete sampling within the cone uses a concentric-ring layout:

  • 1 central ray at the nominal incident angle
  • 12 equally spaced rays per ring

Effective ray count = 1 + 12 × ringCount

ringCount ranges from 2 to 20, giving 25 to 241 effective rays. Increasing ringCount improves angular integration accuracy at linearly increasing computational cost.

Angular distribution weights

DistributionWeighting ruleTypical use
UniformEqual weight per solid-angle elementUniformly divergent beams
LambertianWeight proportional to cos(θ)Diffuse illumination

Scope of applicability

Cone angle averaging applies only to certain result types:

Supported result typesUnsupported result types
R, T, A, Layer AbsorptionPsi / Delta (ellipsometry)
Reflection / Transmission / Absorption spectraDepth distributions (Poynting Vector, Electric Field, Absorption Density)
Reflection / Transmission / Absorption colorDispersion detectors

Unsupported result types are still computed at the single nominal incident angle.

Validity constraint

The cone must not extend to grazing incidence:

incidentAngle + θ_cone < 89.9°

If this condition is violated, the software reports a validation error and blocks computation.

Birefringence and extended parameters

Birefringent materials have different complex refractive indices along different principal directions.

ParameterMeaningTypical use
n, kReal and imaginary parts for the ordinary axisBase optical constants
nExt, kExtReal and imaginary parts for the extraordinary axisDirection-dependent response

In the current app, birefringence can be represented in two ways:

  1. Constant Birefringence: enter one set of ordinary and extraordinary constants directly;
  2. File: import wavelength-dependent birefringent data.

The implementation limits are explicit:

  • surrounding media do not expose birefringent parameters;
  • incoherent layers cannot use the constant birefringent model;
  • incoherent layers also cannot use birefringent refractive-index files.

So in practice, birefringence is intended for coherent thin-film modeling, not for thick-substrate or incoherent approximations.

Physical definition of Psi and Delta

Ellipsometric outputs are not independent quantities. They are derived from the ratio of reflection coefficients:

  • tan(Psi) = |r_p / r_s|
  • Delta = arg(r_p) - arg(r_s)

This gives the two curves distinct roles:

QuantityPhysical meaningSensitivity
PsiAmplitude ratio between p and s reflectionSensitive to index contrast, thickness, and angle
DeltaPhase difference between p and s reflectionEspecially sensitive to interference and phase changes

In engineering work, Psi / Delta are commonly used for thickness and refractive-index inversion. The app follows the same physical requirement: the stack must remain coherent. If any enabled layer is incoherent, Psi / Delta are blocked.

Mapping the concepts to UI fields

Physical quantityMain field or result pageWhat to watch first
Complex refractive indexStructure page: indexType, n, k, file dataPeak position, peak height, absorption strength
Polarization mixOptics page: pRatioWhether s / p behavior transitions smoothly
Incident angleOptics page: Incident AngleSpectral shifts and ellipsometric curve changes
Cone angle averagingOptics page: Cone Angle settingsWhether R / T / A curves are smoothed by angular spread
BirefringenceStructure page: nExt / kExt or birefringent fileWhether polarization-dependent separation increases
Psi / DeltaEllipsometry result pageSensitivity of the curves to small thickness or index changes
Copyright © 2026 Dreapex