Skip to content

Material Library

Inspect a small set of sourced optical constants, compare their valid ranges and use a selected dataset in a film. Measured tables, self-consistent reconstructions and fitted formulas retain their source identities.

Where optical constants come from

The optics experiments here need two numbers for each material: how much it slows light down and how much it absorbs. For an ideal glass you can type 1.5 and move on. For silicon or gold the numbers change with wavelength, depend on how the sample was made, and were measured by somebody over a limited range. This page holds six such datasets from the refractiveindex.info database and keeps each one attached to its source, its conditions and the wavelengths it actually covers.

What n and k describe

The site uses the time dependence exp⁡(−iωt)\exp(-i\omega t) and a complex refractive index n~=n+ik\tilde n=n+ik. A plane wave travelling along +z+z is then

E(z,t)=E0 ei(n~k0z−ωt)=E0 e−kk0z ei(nk0z−ωt),k0=2πλ\begin{aligned}E(z,t)&=E_0\,e^{i(\tilde n k_0 z-\omega t)}\\ &=E_0\,e^{-k k_0 z}\,e^{i(n k_0 z-\omega t)},\\ k_0&=\frac{2\pi}{\lambda}\end{aligned}

with λ\lambda the vacuum wavelength. The real part nn sets the phase velocity c/nc/n and shortens the wavelength inside the material to λ/n\lambda/n. The imaginary part kk makes the wave decay; with this sign convention a passive material has k≥0k\ge 0. Intensity goes as ∣E∣2|E|^2, so it falls by a factor of ee over

ℓ=λ4πk\ell=\frac{\lambda}{4\pi k}

which is the inverse of the absorption coefficient α=4πk/λ\alpha=4\pi k/\lambda. The solvers work with permittivity instead. For a nonmagnetic material

ε=n~2,Re⁡ε=n2−k2,Im⁡ε=2nk\begin{aligned}\varepsilon&=\tilde n^{2},\\ \operatorname{Re}\varepsilon&=n^{2}-k^{2},\qquad \operatorname{Im}\varepsilon=2nk\end{aligned}

The Optical quantity menu below shows all four. The three metals here have kk larger than nn in the visible, so their real permittivity is negative. Silver at 550 nm has n=0.06n=0.06, which puts its phase velocity far above cc. That breaks nothing; phase velocity carries no signal.

Derivation from Maxwell's equations

Start from the macroscopic Maxwell equations for a linear, isotropic, nonmagnetic medium with no free charge, with any conduction current included in ε(ω)\varepsilon(\omega). Fields vary as exp⁡(−iωt)\exp(-i\omega t), units are SI, ε(ω)\varepsilon(\omega) is the relative permittivity, k0=ω/ck_0=\omega/c is the vacuum wave number and Y0=ε0/μ0Y_0=\sqrt{\varepsilon_0/\mu_0} is the vacuum admittance.

∇×E=iωμ0H,∇⋅(εE)=0,∇×H=−iωε0εE,∇⋅H=0.\begin{aligned}\nabla\times\mathbf E&=i\omega\mu_0\mathbf H, & \nabla\cdot(\varepsilon\mathbf E)&=0,\\ \nabla\times\mathbf H&=-i\omega\varepsilon_0\varepsilon\mathbf E, & \nabla\cdot\mathbf H&=0.\end{aligned}

In a uniform region ε\varepsilon is constant, so ∇⋅E=0\nabla\cdot\mathbf E=0. Taking the curl of Faraday's law and substituting the Ampère–Maxwell law gives the Helmholtz equation:

∇×∇×E=εk02 E,∇×∇×E=∇(∇⋅E)−∇2E=−∇2E,∇2E+εk02 E=0.\begin{aligned}\nabla\times\nabla\times\mathbf E&=\varepsilon k_0^2\,\mathbf E,\\ \nabla\times\nabla\times\mathbf E&=\nabla(\nabla\cdot\mathbf E)-\nabla^2\mathbf E\\ &=-\nabla^2\mathbf E,\\ \nabla^2\mathbf E+\varepsilon k_0^2\,\mathbf E&=0.\end{aligned}

Plane wave. Try E=x^E0 ei(n~k0z−ωt)\mathbf E=\hat{\mathbf x}E_0\,e^{i(\tilde n k_0 z-\omega t)}. The Helmholtz equation gives −n~2k02+εk02=0-\tilde n^2k_0^2+\varepsilon k_0^2=0, so n~2=ε\tilde n^2=\varepsilon. Of the two roots, take the one with k≥0k\ge0 (and n>0n>0 when ε\varepsilon is real and positive), so the wave decays, or at least does not grow, as it travels along +z+z. For a passive medium with Im⁡ε≥0\operatorname{Im}\varepsilon\ge0 that is the principal square root the solvers use; for real negative ε\varepsilon it is n~=i∣ε∣\tilde n=i\sqrt{|\varepsilon|}. Writing ε=ε1+iε2\varepsilon=\varepsilon_1+i\varepsilon_2 and expanding (n+ik)2(n+ik)^2 gives ε1=n2−k2\varepsilon_1=n^2-k^2 and ε2=2nk\varepsilon_2=2nk.

Power flow. With E\mathbf E along xx, the magnetic field points along yy and the flux along zz. Faraday's law gives iωμ0Hy=∂zEx=in~k0Exi\omega\mu_0H_y=\partial_zE_x=i\tilde nk_0E_x, and since k0/(ωμ0)=Y0k_0/(\omega\mu_0)=Y_0, Hy=n~ Y0ExH_y=\tilde n\,Y_0E_x. The time-averaged Poynting flux ⟨Sz⟩=12Re⁡(ExHy∗)\langle S_z\rangle=\tfrac12\operatorname{Re}(E_xH_y^*) is then

⟨Sz⟩=12 n Y0 ∣E0∣2 e−2kk0z,ℓ=12kk0=λ4πk,α=1ℓ\begin{aligned}\langle S_z\rangle&=\tfrac12\,n\,Y_0\,|E_0|^2\,e^{-2kk_0z},\\ \ell&=\frac{1}{2kk_0}=\frac{\lambda}{4\pi k},\qquad \alpha=\frac1\ell\end{aligned}

The power falls by ee over ℓ\ell. The field amplitude falls by ee over twice that, λ/(2πk)\lambda/(2\pi k), so check which of the two a quoted "penetration depth" means.

Where the power goes. Poynting's theorem for time-harmonic fields says the flux lost per unit volume is the power dissipated there, −∇⋅⟨S⟩=12ωε0ε2∣E∣2-\nabla\cdot\langle\mathbf S\rangle=\tfrac12\omega\varepsilon_0\varepsilon_2|\mathbf E|^2. For the plane wave both sides agree:

−d⟨Sz⟩dz=2kk0⟨Sz⟩=nk k0Y0∣E∣2=12 ωε0 (2nk) ∣E∣2=12 ωε0 ε2∣E∣2\begin{aligned}-\frac{d\langle S_z\rangle}{dz}&=2kk_0\langle S_z\rangle=nk\,k_0Y_0|E|^2\\ &=\tfrac12\,\omega\varepsilon_0\,(2nk)\,|E|^2\\ &=\tfrac12\,\omega\varepsilon_0\,\varepsilon_2|E|^2\end{aligned}

using k0Y0=ωε0k_0Y_0=\omega\varepsilon_0. So ε2>0\varepsilon_2>0 means the material takes energy from the wave, and ε2=0\varepsilon_2=0 means it takes none.

Why the k=0k=0 choice matters. With ε\varepsilon real and positive, n~\tilde n is real, the flux is constant and nothing is absorbed at any thickness. The transparent-band model sets ε2=0\varepsilon_2=0 for silica. That is a statement about the model, not about the glass.

With ε\varepsilon real and negative, n~\tilde n is purely imaginary: the field still decays, but n=0n=0 so it carries no power and nothing is absorbed. A thick enough film of such a material is opaque because it reflects, not because it absorbs.

In a film. At incidence angle θ\theta from a medium of index n0n_0, the same Helmholtz equation gives the normal wave number q=ε−n02sin⁡2θq=\sqrt{\varepsilon-n_0^2\sin^2\theta} in units of k0k_0, which is how the thin-film page uses the permittivity from this library.

Tables and formulas

A table lists nn and kk at discrete wavelengths. Gold from Johnson and Christy has 49 points between 187.9 and 1937 nm; silicon from Green has 121, every 10 nm from 250 to 1450 nm. Between points this page interpolates nn and kk linearly. Beyond the first and last point it refuses. A table also records choices: the TiO₂ film table lists k=0k=0 exactly for its last 198 points, from 1014.7 to 1689.8 nm. The file does not say why; the step from 10−610^{-6} to zero looks like rounding or a model cut-off. Read them as no reported loss, not as a measurement of perfect transparency.

Two tables are reconstructions rather than direct measurements. Green's silicon combines datasets while enforcing Kramers–Kronig consistency, the causality relation that ties nn to kk (see Kramers–Kronig). Rakić's aluminium comes from an iterative, self-consistent algorithm that combines Kramers–Kronig analysis with a semiquantum model, applied to optical and electron-energy-loss data for evaporated films. The source records keep those labels.

A formula stores a few coefficients. The fused silica record has the Sellmeier form

n2=1+∑i=13Bi λ2λ2−Ci2n^{2}=1+\sum_{i=1}^{3}\frac{B_i\,\lambda^{2}}{\lambda^{2}-C_i^{2}}

with Bi=0.6961663, 0.4079426, 0.8974794B_i=0.6961663,\ 0.4079426,\ 0.8974794 and Ci=0.0684043, 0.1162414, 9.896161 μmC_i=0.0684043,\ 0.1162414,\ 9.896161\ \mu\text{m}. Each term behaves like a fitted resonance at λ=Ci\lambda=C_i: two in the vacuum ultraviolet (68 and 116 nm), one in the infrared (9.9 µm). A formula will return a number almost anywhere you ask, which is the problem. Malitson fitted it between 0.21 and 3.71 µm, and the source record says Tan verified it up to 6.7 µm. Evaluated at 9 µm, next to the infrared term, the same coefficients give n2=−2.2n^2=-2.2. The page only evaluates it from 210 to 6700 nm.

The formula also has no kk. That does not make silica lossless; loss was simply not part of the fit. Malitson's index is relative to air and his wavelengths are in air. Polyanskiy's 2024 descriptor uses this record as its example of the fields n_is_absolute: false\texttt{n\_is\_absolute: false} and wavelength_is_vacuum: false\texttt{wavelength\_is\_vacuum: false}; the pinned file bundled here has neither. The page does not convert. Multiplying by the index of dry standard air at 15 °C from Ciddor's equation, nair≈1.00028n_\text{air}\approx1.00028, gives an absolute index about 0.03% higher: 1.4589 instead of 1.4585 at 587.6 nm. Malitson measured at 20 °C; warmer air or the small wavelength correction changes that by less than 10−510^{-5}, too little to show in the film results here.

Three things to try

Silica without a k. The page opens on fused silica. It gives n=1.4585n=1.4585 at 587.6 nm and 1.4440 at 1550 nm. The film buttons stay disabled until you tick the transparent-band model, because k=0k=0 is your assumption, not the source's. With it, the 100 nm film sent to the thin-film page (air above, a substrate of index 1.5, normal incidence) reflects 3.0% at 550 nm.

Silicon across the visible. Uncheck silica and select silicon, then switch to kk (the record's kk covers band-to-band absorption only). It drops from 0.296 at 400 nm to 0.0280 at 550 nm and 0.0054 at 800 nm. The decay lengths are 108 nm, 1.6 µm and 11.8 µm. Sent to the thin-film page, the 100 nm film absorbs 37.7% at 400 nm (set the wavelength on the thin-film page), 3.0% at 550 nm and 1.0% at 800 nm. That is not simply 1−e−d/ℓ1-e^{-d/\ell} (60%, 6.2% and 0.85%), with dd the film thickness: the front surface reflects part of the light first, and interference between passes redistributes the rest. At 800 nm the film absorbs slightly more than one pass would.

Three metals. Untick every selected dataset, then select gold, silver and aluminium, in that order. The chart uses their common range, 187.9–1937 nm, even though the aluminium table runs from 0.124 nm to 200 µm. Switch to real permittivity: at 550 nm (interpolated between table points) it is −5.9 for gold, −12.9 for silver and −42.9 for aluminium. The intensity decay lengths there are 17.7, 12.2 and 6.6 nm, so a 100 nm film is effectively opaque. Send gold, the first selection, to the thin-film page: at 550 nm the 100 nm film reflects 78.6%, absorbs 20.7% and transmits 0.7%, and at 400 nm it absorbs 59%. Come back to this page, untick everything, select silver on its own and send it: it reflects 98.3%.

Choose source datasets

Select up to three datasets. The first selected entry is used for handoff.

Optical constants

Comparison uses the common source range 210.00–6700.0 nm. Curves interpolate n and k linearly between original knots; formulas are evaluated directly. No extrapolation.

Selected material data
1.541.150.7690.3850210.06700Vacuum wavelength (nm)n

Fused silica · Malitson

Sample values
Vacuum wavelength (nm)Fused silica · Malitson (n)
210.001.5384
426.331.4675
642.671.4567
859.001.4524
1075.31.4495
1291.71.4470
1508.01.4445
1724.31.4419
1940.71.4389
2157.01.4357
2373.31.4321
2589.71.4281
2806.01.4236
3022.31.4187
3238.71.4133
3455.01.4072
3671.31.4005
3887.71.3932
4104.01.3850
4320.31.3760
4536.71.3659
4753.01.3548
4969.31.3423
5185.71.3284
5402.01.3128
5618.31.2952
5834.71.2752
6051.01.2524
6267.31.2261
6483.71.1956
6700.01.1596

Silica supplies n only. Missing k and derived complex permittivity are unavailable in this comparison; they are not filled with zero.

Source records

Fused silica · Malitson

I. H. Malitson, C. Z. Tan. Interspecimen comparison of the refractive index of fused silica (1965); infrared-range verification by Tan (1998). Fused silica at 20 °C. Malitson fit 0.21–3.71 µm; Tan extension to 6.7 µm. No loss data. Sellmeier fit to refractive-index data. Source files use µm; wavelength converted to nm. No extrapolation.

Valid range: 210–6700 nm. Original dataset file. Database revision: c5c2f188e848453def5970e347399d653df2ffc2. License: CC0-1.0 (database).

Legacy research tables

The original thesis/paper dispersion tables remain available in the coherence experiment, with their existing source attribution and unsupported negative-loss regions. They are not silently repaired or relabelled as the curated datasets above. Tensor/anisotropic data imports are not supported by these scalar builders.

Acquisition and references

These files are pinned from the refractiveindex.info database. Original files, SHA-256 hashes and acquisition metadata are in the import manifest; the database’s CC0 license is retained.

Mikhail N. Polyanskiy, Refractiveindex.info database of optical constants, Scientific Data 11, 94 (2024) describes the compilation. Each selected source record above credits the original work. Inspection of the imported data and metadata supports the numerical inputs; these entries are not detailed summaries of the original papers. Source conditions matter: a deposited film does not automatically represent a bulk material.

W. Sellmeier, Zur Erklärung der abnormen Farbenfolge im Spectrum einiger Substanzen, Annalen der Physik 219(6), 272–282 (1871). Read in full, in German, from the OCR text (its tables are partly garbled) of the Internet Archive scan.

W. Sellmeier, Ueber die durch die Aetherschwingungen erregten Mitschwingungen der Körpertheilchen und deren Rückwirkung auf die ersteren, besonders zur Erklärung der Dispersion und ihrer Anomalien, Annalen der Physik 223(11), 386–403 (1872), one part of a series that also appeared in volume 221. §8 was read from the OCR text of the Internet Archive scan; its printed equations are not legible there.

I. H. Malitson, Interspecimen comparison of the refractive index of fused silica, Journal of the Optical Society of America 55(10), 1205 (1965). Only the abstract was read; the full text was not consulted.

C. Z. Tan, Determination of refractive index of silica glass for infrared wavelengths by IR spectroscopy, Journal of Non-Crystalline Solids 223(1–2), 158–163 (1998). Database record and metadata only; cited for the source record's statement that the silica formula holds to 6.7 µm.

P. B. Johnson and R. W. Christy, Optical constants of the noble metals, Physical Review B 6(12), 4370–4379 (1972). Only the abstract was read; the full text was not consulted.

M. A. Green, Self-consistent optical parameters of intrinsic silicon at 300K including temperature coefficients, Solar Energy Materials and Solar Cells 92(11), 1305–1310 (2008). Database record and metadata only.

A. D. Rakić, Algorithm for the determination of intrinsic optical constants of metal films: application to aluminum, Applied Optics 34(22), 4755–4767 (1995). Database record and metadata only; the method description comes from the comment in that record, which appears to reproduce the paper's abstract.

S. V. Zhukovsky, A. Andryieuski, O. Takayama, E. Shkondin, R. Malureanu, F. Jensen and A. V. Lavrinenko, Experimental demonstration of effective medium approximation breakdown in deeply subwavelength all-dielectric multilayers, Physical Review Letters 115(17), 177402 (2015). Database record and metadata only; numerical data provided to the database by Osamu Takayama.

P. E. Ciddor, Refractive index of air: new equations for the visible and near infrared, Applied Optics 35(9), 1566–1573 (1996). Not read; used only for the index of standard air, evaluated from the coefficients in the refractiveindex.info air record (dry air, 15 °C, 101.325 kPa, 450 ppm CO₂).

History of the method

Sellmeier's 1871 article starts from the observations by Christiansen and Kundt that a few substances, among them a fuchsine solution, spread the colours of the spectrum in an abnormal order. Published in the Annalen der Physik, it says that he had expected this. His explanation was that light travels at a different speed in matter because it sets the material's particles vibrating along with it. That links refraction to absorption: particles that absorb light of one period raise the index for light of longer period and lower it for shorter periods, more and more strongly near the absorbed period. He had looked for the effect in a fuchsine solution in 1866 and not found it.

In a series of papers in 1872 he used conservation of energy, in the elastic-ether picture of the time, to write the refracting power n2−1n^2-1 as a sum over the particles' vibrations, each with its own natural period. Whether a particle's contribution adds or subtracts depends on whether the light's period is longer or shorter than that natural period. The resonance form used for silica on this page carries his name. The equations in the scan available here were not legible, so this page does not claim its notation matches his term for term.

In 1965 I. H. Malitson measured the index of optical-quality fused silica at 60 wavelengths from 0.21 to 3.71 µm at 20 °C, compared 12 specimens, and published a three-term dispersion equation. Its coefficients are the ones bundled here. In 1972 P. B. Johnson and R. W. Christy obtained nn and kk for copper, silver and gold from reflection and transmission of vacuum-evaporated films 185–500 Å thick, between 0.5 and 6.5 eV. Three measurements were inverted to give the film thickness as well as nn and kk, with an estimated error in nn and kk below 0.02 over most of that range.

Polyanskiy's descriptor places refractiveindex.info after printed compilations such as Palik's Handbook of Optical Constants of Solids, dates the database to 2008, and says its first list of sources came from the references in the Handbook of Optics. It is also direct about the limits: the database produces no measurements and cannot verify every record, so it relies on peer review and checks its own transcription instead.

Calculation source code · Implementation and validation notes