A soap film is a sheet of water between two layers of soap molecules, often thinner than a micrometer. A vertical film drains under gravity and is thinner at the top, so in white light it shows bands of color ordered by thickness. Newton watched them on bubbles kept under a glass cover (Opticks, Book II, Part I, Observation 17): colored rings formed at the top, widened and moved down as the water drained, and a black spot grew at the top before the bubble broke.
The color at thickness follows from two ingredients: the fraction of light the film reflects at each wavelength , and the three color-matching functions that the CIE (the International Commission on Illumination) uses to reduce a spectrum to a color. The first is a single cosine in , so the color is a fixed white minus the cosine transform of the color-matching functions, evaluated at , the extra optical path traveled by light reflected from the film’s back surface, where is the refractive index of water. The computed colors agree with the table of thicknesses that Newton compiled for them.
The reflected spectrum is a cosine
Consider light of wavelength that travels normal to the film, perpendicular to its two surfaces. Two beams return toward the eye. One is reflected by the front surface. The other crosses the film, is reflected by the back surface, and crosses the film again.
Light is an electromagnetic wave: an electric field and a magnetic field that oscillate together as the wave travels. The amplitude of a beam is the size of its electric field’s oscillation, and the brightness of the beam is proportional to the square of its amplitude. In the figure, each beam is labeled with its amplitude as a multiple of the incoming beam’s amplitude.
All the beams in the figure travel along the same line, normal to the film. They are drawn side by side so that each one can be labeled. The right half shows the two returning waves at one instant. The back-surface wave is the front-surface wave turned upside down, and it trails the front-surface wave by a distance , derived below. The light that reaches the eye is the sum of the two waves.
How much each surface reflects. When light travels normal to a surface, its electric and magnetic fields both lie parallel to the surface. In a medium of refractive index , the amplitude of the magnetic field is proportional to times the amplitude of the electric field. Suppose a beam of amplitude crosses a boundary from a medium of index into a medium of index . It splits into a reflected beam of amplitude and a transmitted beam of amplitude . Maxwell’s equations require each field to take the same value on both sides of the boundary. For the electric field, this condition reads . The reflected beam travels in the opposite direction, so its magnetic field points the opposite way, and the condition for the magnetic field reads . Solving these two equations gives the Fresnel coefficients for normal incidence,
In the soap film, the index is for air and for water. At the front surface, light passes from air into water, so and . For , . At the back surface, light passes from water into air. The two indices trade places, and the reflection coefficient becomes . The two reflections therefore have opposite signs. On its way out, the back-surface beam passes through the front surface from water into air, with transmission coefficient . This beam is transmitted, reflected and transmitted again, so it returns with amplitude
The last step uses . For water, . From here on I approximate , which changes the back-surface beam’s amplitude by , so the two returning beams have amplitudes and .
How far the second beam trails the first. The back-surface beam travels an extra distance inside the water, where the wavelength shortens to . It therefore falls behind the front-surface beam by wavelengths, where
is the optical path difference between the two beams. In air, where the wavelength is , this is a distance , as in the figure. Each wavelength of delay is one full cycle of the oscillation, so the back-surface wave is shifted in its cycle by the angle . This angle is called the phase difference.
Two waves with different phases are added by writing each one as a complex number whose magnitude is the amplitude and whose angle is the phase, as in the Euler’s identity post. A shift by the angle multiplies the complex number by . The light reflected by the film therefore has amplitude
The reflectance is the fraction of the incoming light that the film reflects. Brightness is proportional to the square of the amplitude, so is the squared magnitude of this complex number. Since ,
Taking constant across the visible band, , so lies between and . In the wavenumber , is a single cosine with period . Light can also bounce back and forth inside the film several times before it leaves. Summing all of these reflections gives the exact Airy formula (Born and Wolf, §7.6). Its higher harmonics , for integers , shrink by a factor from each to the next, and the two-beam formula stays within about of it.
Three numbers from a spectrum
For the CIE 1931 standard observer, light with spectral power is described by three tristimulus values
where are the tabulated color-matching functions. Two lights with the same look identical to the standard observer, and since is the eye’s sensitivity to brightness, is the luminance. Light reflected by the film has for an illuminant with spectrum . I use D65, the CIE’s standard representation of average daylight (Newton viewed his bubbles “by Reflexion of the Skies when of a white Colour”), normalized so that a perfect reflector has the white point . The sRGB standard specifies the matrix and the encoding curve that take to the color a screen displays.
The color is a cosine transform
Substitute into the three integrals and change variables to . Define three weight curves on the wavenumber axis by
and , likewise, so that their integrals are the three components of . With , the tristimulus values of the film are
The vector is the cosine transform (the real part of the Fourier transform) of the three weight curves, evaluated at the path difference.
The black film. At the cosine equals at every wavelength, so and the film reflects nothing: the two reflections cancel exactly. For small , is proportional to , so the first light reflected by a thin film is weighted toward short wavelengths. Its chromaticity, in the CIE coordinates , is , bluer than the white point at . This bluish gray is the faint blue that Newton lists right after black. At nm the film reflects of the light, less than a tenth of the reflected by a single water surface but not zero, and Newton could see “the Images of some things (as of a Candle or the Sun) very faintly reflected” from the black spot.
The orders. The luminance is smallest when the wavelengths at the center of are extinguished, that is, when is a multiple of nm, the reciprocal of the mean wavenumber of . The computed minima are at , , and nm, spaced by nm. With the middle of the spectrum removed, red and blue light remain, and these purple bands separate the orders: the second and third orders run through blue, green, yellow and red, and from the fourth order on, greens and pinks alternate and grow paler.
The widget evaluates the formula from the CIE tables at nm steps: the reflectance with the products , , dashed, the chromaticity diagram with the film’s path up to the current , and the color for every up to m.
In the first order the peak of , at , moves from the ultraviolet into the visible band. At nm it sits at nm, stays above half its maximum across the band, and the film is white. The zero of at enters the band at nm, which turns the film yellow and then red, and reaches the green at nm. The chromaticity path goes once around the white point in each of the first three orders. After that, it swings between the green and pink sides of the white point, at a distance that shrinks in every order after the second.
The colors fade to white
As grows, by the Riemann–Lebesgue lemma (the Fourier transform of an integrable function tends to zero at infinity), so a thick film reflects the illuminant’s own color at , the reflectance of two water surfaces whose reflections add without interfering. The cosine in has period in , and the visible band holds of its periods: at m and at m. A weight curve that spans several periods averages the cosine to nearly zero.
The widths of the weight curves set the rate of decay. Replacing by a Gaussian with the same mean and standard deviation gives
the oscillation that produces the orders, under a Gaussian envelope. For the tabulated D65 and , , and the envelope falls to at m, a thickness of m. The modulus of the exact transform of the tabulated curve stays within of this envelope up to m. The transform of a narrower weight curve would decay more slowly. This is the width trade-off of the Fourier transform from the Heisenberg post. Computed from the formula, the chromaticity stays within of the white point for every m and within for every m.
Newton’s table
Newton measured the rings of color in the thin film of air between two lenses pressed together, and compiled a table of the thickness “at which each Colour is most intense and specifick”, in millionths of an inch, for films of air, water and glass (Book II, Part II). Selected entries of his water column, converted at nm per millionth of an inch, next to the computed color at that thickness:
| order | Newton’s color | Newton (10⁻⁶ in) | (nm) | computed |
|---|---|---|---|---|
| 1st | white | 3⅞ | 98 | |
| 1st | yellow | 5⅓ | 135 | |
| 1st | red | 6¾ | 171 | |
| 2nd | violet | 8⅜ | 213 | |
| 2nd | blue | 10½ | 267 | |
| 2nd | green | 11⅔ | 296 | |
| 3rd | purple | 15¾ | 400 | |
| 3rd | green | 18⁹⁄₁₀ | 480 | |
| 4th | green | 26½ | 673 | |
| 7th | ruddy white | 57¾ | 1467 |
The brightest computed film, at nm, is Newton’s white of the first order at nm. His violet of the second order, at nm, sits next to the first luminance minimum at nm, and his last entry, the ruddy white of the seventh order at nm, comes out as a pinkish white near the white point.
Newton also viewed the bubble obliquely (Observation 19). At angle to the normal inside the film the path difference is , so a color seen at thickness at normal incidence appears where the film is thick. Snell’s law gives , so the largest factor, at grazing incidence, is for water. Newton’s ratio for his most oblique view: “By the best of my Observations it was between 15 and 15-1/2 to 10.”
The Michel-Lévy chart
A plate of a birefringent crystal between crossed polarizers, at to them, transmits the fraction of the polarized light, where the retardation is the optical path difference between the plate’s two polarizations. That is with , so the strip in the widget is also the interference color chart that mineralogists use to identify minerals in thin section, with as the retardation (for the same illuminant, and with the birefringence taken constant across the band). Sørensen (2013) recomputed that chart from first principles.
References
- Isaac Newton. Opticks: or, a Treatise of the Reflections, Refractions, Inflections and Colours of Light. London, 1704. Fourth edition, 1730. Book II, Part I, Observations 17–19 on soap bubbles, and Part II, the table of thicknesses of colored plates.
- Max Born and Emil Wolf. Principles of Optics, 7th edition. Cambridge University Press, 1999. Chapter 7: two-beam interference in a plane-parallel plate (§7.5) and multiple-beam interference (§7.6).
- CIE. Colour-matching functions of CIE 1931 standard colorimetric observer. Dataset, 2019, from ISO/CIE 11664-1:2019.
- CIE. CIE standard illuminant D65. Dataset, 2019, from ISO/CIE 11664-2:2022.
- IEC 61966-2-1:1999. Multimedia systems and equipment – Colour measurement and management – Part 2-1: Colour management – Default RGB colour space – sRGB.
- Bjørn Eske Sørensen. A revised Michel-Lévy interference colour chart based on first-principles calculations. European Journal of Mineralogy, 25(1):5–10, 2013.