Back to Academy

Optical Aberrations: The Five Seidel Terms, on a Real Lens

What each one does to an image, traced live on a real camera lens from an 1896 patent.

Key Takeaways

  • Even a perfect lens cannot focus a star to a point: diffraction alone makes a disk about 4.6 µm in radius here. This 1896 camera lens, on axis at full aperture, spreads it to a radius of 188 µm.
  • In one color, and to lowest order, rotational symmetry allows exactly five ways for a wavefront to go wrong: spherical aberration, coma, astigmatism, field curvature and distortion.
  • Each grows at its own rate with the aperture and with the angle off axis, and those rates are how a designer tells, from a ray fan or a spot diagram, which one is doing the damage.

The gap: four and a half micrometers, or a hundred and eighty-eight

Take a lens with a focal length of 98.27 mm and an aperture of 15.38 mm — the patent calls it f/6.5, and its printed radii give f/6.4 — and light it with a single color, the sodium D line at 589.3 nm. If the lens were perfect, a star on its axis would image as a bright disk with a first dark ring at a radius of 4.59 µm, the size diffraction alone would make it. That is the floor physics sets.

The floor, and the f-number

The first dark ring of the Airy pattern sits at 1.22 λ f/# from the center; at λ = 0.5893 µm and f/6.39 that is 4.59 µm, and every spot diagram on this page draws that circle so the blur has a scale. The patent rates this lens f/6.5, and that is the name it carries in the Ray-Optics Designer, the lens-design tool the widgets on this page run on. The entrance pupil — the width of the beam the lens actually accepts — that its printed numbers produce is 15.38 mm on a 98.27 mm focal length, which is f/6.4, and this page uses the number the lens actually has.

Now trace real rays through a real lens of exactly those numbers: the Cooke triplet from H. D. Taylor’s 1896 patent, US 568,052, read straight from the printed table. The 127 rays the Ray-Optics Designer sends through the lens land on the image plane in a blob whose outermost ray is 187.7 µm from the axis. Forty-one times the diffraction floor, on axis, in one color.

marginal raychief ray, 26° fieldlight in, from an object at infinitybeam in, 15.38 mm(the entrance pupil)aperture stop123456to the image plane,82 mm beyond surface 610 mmL1 crown, n = 1.6114L2 flint, n = 1.5679L3 crown, n = 1.6114
The lens every number on this page comes from: the Cooke triplet of US 568,052, drawn to scale from the same six surfaces the widgets trace, with the edge ray of an on-axis beam and the chief ray of the 26° field. The image plane is 82 mm beyond the last surface, off the figure.

What you are looking at. A Cooke triplet is three separate pieces of glass on one axis, a positive element, a negative one, and a second positive one, with the aperture stop, the ring that decides how wide a bundle of light the lens accepts, on the back surface of the middle element. It is the simplest lens with enough adjustable quantities, six curvatures and two air gaps, to hold the focal length and bring all seven primary aberrations under control at once: the five on this page and two color errors. That is why it is still the first lens a design textbook works through, and the ancestor of a large family of camera lenses.

Three pieces of glass mean six curved surfaces, numbered 1 to 6 from the object side in the diagram above, and every table on this page is a sum of one contribution per surface: 1 and 2 are the front crown, 3 and 4 the flint in the middle, 5 and 6 the rear crown. Crown and flint are the two broad families of optical glass, told apart by how strongly they spread colors, not by density or index: a flint fans white light into a wider spectrum than a crown does. That stronger color spread is what lets a negative flint cancel the color error of two positive crowns without canceling their focusing power. On this page, in one color, the middle element’s visible job is a different one: flattening the field, which the ledger below shows surface by surface.

The radii, the thicknesses and the two glasses’ refractive indices, the lens’s prescription, are the numbers printed in H. D. Taylor’s 1896 patent, US 568,052. Nothing on this page is fitted, tuned or rounded to look better.

The patent’s numbers, and which glass is the flint

The patent gives them with the focal length as the unit; here they are in millimeters, one hundred times larger, which is how a lens of 98 mm focal length comes out of a table rounded to three figures. And Taylor’s flint is in fact the lower-index glass of the two, which is the opposite of the way the word is usually read: the two families are told apart by how strongly they spread color, not by index.

Where did the factor of forty-one go? Into five named places. Not because a lens has five kinds of flaw, but because, to lowest order, rotational symmetry allows exactly five ways for a wavefront to depart from a sphere. They are spherical aberration, the blur that is there even on axis; coma, the comet tail off axis; astigmatism, two focus distances and no sharp plane; field curvature, a sharp surface that is not flat; and distortion, the image moves but does not blur. The rest of this page takes them one at a time: what each does to a point image, how each grows with aperture and with field angle, and how a designer reads them off a real prescription. The lens stays on the page throughout: every widget below re-traces those same six surfaces in your browser, with the same code the Ray-Optics Designer runs, and you move the sliders.


What an aberration is

A lens that images a point perfectly turns the light from that point into a spherical wavefront converging on one image point. Any real lens produces a wavefront that is not quite spherical. The wavefront aberration W is the gap between the two: pick a reference sphere centered on the ideal image point and passing through the center of the exit pupil (the lens aperture as seen from the image side), and measure, ray by ray, how far the actual wavefront sits from that sphere along the ray. W is a function of where the ray crosses the pupil.

Exit pupilImage planeaReference sphereWavefrontRchief rayyp = a·PyWεy
The wavefront aberration W is the gap between the actual wavefront and a reference sphere centered where the chief ray meets the image plane. A ray runs perpendicular to its wavefront, so the tilt W picks up across the pupil is the miss εy at the image. Schematic: the departure is exaggerated by several orders of magnitude.

A few words carry the rest of the page. A ray is labeled by where it crosses the pupil, from 0 at the center to 1 at the rim.

Field point
: a direction in object space, its angle θ from the axis, always set in the vertical plane here, so “up” in the image is the field direction; this lens covers fields out to 26°.
Chief ray
: the ray through the center of the pupil.
Marginal ray
: the ray at the rim, in the vertical plane.
Tangential
: the vertical section of the pupil. Sagittal is the horizontal one, and most aberrations treat the two differently.
Paraxial
: computed for rays so close to the axis that the lens is ideal.

The wavefront picture and the ray picture are the same information. A ray runs perpendicular to its wavefront, so a tilt in W is a bend in the ray, which is a miss at the image: the transverse ray aberration ε, where a ray meets the image plane minus where the chief ray of the same field meets it. With R the radius of the reference sphere, n′ the refractive index in image space, and (x_p, y_p) the ray’s position in the exit pupil,

εy=RnWyp,εx=RnWxp

Everything else on this page, the sums, the fans and the spots, is that one derivative applied to one polynomial.

Coordinates and signs, stated once

Because every sign on this page depends on them. Light travels toward +z; the lens is rotationally symmetric about z; the object is at infinity and a field point is specified by its angle θ from the axis in the y–z plane, so y is the tangential direction and x the sagittal one. A ray is labeled by its fractional pupil position P = (Pₓ, P_y), with ρ = |P| between 0 and 1 and azimuth φ measured from +y, so P_y = ρ cos φ and Pₓ = ρ sin φ; the chief ray is P = 0 and the marginal ray is P_y = 1. The field is normalized as h = tan θ / tan θ_max, 0 on axis and 1 at the edge of the field. Slopes are geometric, u = dy/dz; in image space the marginal ray converges, so its slope u′ is negative, and for this lens u′ = −(15.38/2)/98.268 = -0.07826 with n′ = 1 in air. ε = (εₓ, ε_y) is measured in the paraxial image plane relative to the chief ray of the same field, so ε(0) = 0. W is the optical path of the chief ray minus the optical path of the ray, both measured to the reference sphere; positive W means the ray’s path is shorter than the chief’s. In pupil coordinates of length, y_p = a·P_y with a the exit-pupil semi-diameter; on axis R/a is exactly −1/u′, so in fractional coordinates the derivative reads ε_y = (1/n′u′) · ∂W/∂P_y, and the engine reproduces that relation on the real lens on axis to 1.3 %, the residual being the finite-difference step used in the check. Off axis the exact relation carries obliquity factors from the tilted chief ray, which is one of the reasons the third-order formulas on this page are a small-field theory.


Why exactly five

A rotationally symmetric system cannot tell x from y, or a field point on one side of the axis from its mirror image on the other. So W can depend on pupil and field only through three combinations that survive a rotation: ρ², how far out in the pupil; h², how far out in the field; and h·ρ·cos φ, their product, with φ the angle round the pupil from the field direction. Expand W in those three. The first-order terms are a constant, a focus shift and a change of magnification, absorbed into the choice of image plane and scale. The next order, fourth-order in the coordinates, is where the aberrations begin, and there are exactly five:

W(ρ,φ;h)=18SIρ4+12SIIhρ3cosφ+12SIIIh2ρ2cos2φ+14(SIII+SIV)h2ρ2+12SVh3ρcosφ

This is the third-order, or Seidel, wavefront aberration, and SI to SV are the Seidel sums, in millimeters on this page: spherical aberration, coma, astigmatism, field curvature (the Petzval term) and distortion, in that order. Third-order, because the ray misses that follow from this W are cubic in the coordinates. Sums, because each surface contributes to each of the five separately and the contributions simply add, which is what makes the ledger below readable.

How the sums are computed

Trace a paraxial marginal ray (y, u) and a paraxial chief ray (ȳ, ū) through the system and define, at each surface of curvature c between index n and n′: A = n (u + y c), the refraction invariant of the marginal ray, and Ā = n (ū + ȳ c) for the chief ray; δ(u/n) = u′/n′ − u/n and Δ(1/n) = 1/n′ − 1/n; and H = n (ū y − u ȳ), the Lagrange invariant, the same number at every surface. Then, summing over surfaces, SI = −Σ A² y δ(u/n); SII = −Σ A Ā y δ(u/n); SIII = −Σ Ā² y δ(u/n); SIV = −Σ H² c Δ(1/n); SV = −Σ (Ā/A) [² y δ(u/n) + H² c Δ(1/n)]. These are Welford’s forms in his normalization and what the Ray-Optics Designer evaluates for its Seidel view; they reproduce the rayoptics reference implementation term by term with a scale factor of exactly one. Move the stop and only the chief ray changes: at every surface ȳ becomes ȳ + ε·y for one number ε, so Ā becomes Ā + εA, and substituting that into the five sums gives the stop-shift equations of Moving the stop.

Two conventions are pinned: the chief ray is launched through the center of the entrance pupil (ȳ = −ū·z_EP at the first vertex), which is what makes the per-surface SII, SIII and SV rows meaningful; and a positive SI produces ε_y = +SI ρ³/(2 n′u′), which with u′ negative puts the marginal ray inside the paraxial focus, the classic undercorrected spherical aberration of a positive lens.

Differentiating W along the two sections of the pupil gives the third-order shapes of the two ray fans at any field, with the sums evaluated at that field: tangential ε_y(P) = [SI P³ + 3 SII P² + (3 SIII + SIV) P] / (2 n′u′) with P = P_y, and sagittal εₓ(P) = [SI P³ + (SIII + SIV) P] / (2 n′u′) with P = Pₓ. Distortion is absent from both because the fans are measured relative to the chief ray, which is the ray distortion moves; coma is absent from the sagittal fan because the sagittal rays’ coma displacement is in y, which that fan does not plot. The per-surface algebra and the suppression rule for tilted systems are in the Ray-Optics Designer manual, §6.

Look at how the pupil and field enter each term, because that is the whole art of reading a Seidel table:

TermWavefront WTransverse εPupil powerField power
SI sphericalρ⁴ρ³30
SII comah ρ³h ρ²21
SIII astigmatismh² ρ²h² ρ12
SIV field curvatureh² ρ²h² ρ12
SV distortionh³ ρ03

Every row sums to three powers in ε. Spherical aberration is all aperture and no field; distortion is all field and no aperture; the other three sit between. The exponents are exact in the engine, and the widget prints them.

The Seidel ledger of a real lens

Slide the field to 0° and watch four of the five columns vanish.

SurfaceSISphericalSIIComaSIIIAstigmatismSIVField curvatureSVDistortion
10.16760.13480.10840.31400.3398
20.2365-0.34370.49960.05648-0.8082
3-0.28880.3738-0.4838-0.090990.7439
4-0.09211-0.1300-0.1835-0.3205-0.7114
5-5.853e-8-2.773e-5-0.013140.014730.7563
60.01128-0.028090.069940.06932-0.3468
Sum0.034510.006704-0.0024900.04306-0.02636
with apertureEPDEPD³EPD²EPD²EPD¹
with fieldconstanttan θtan² θtan² θtan³ θ

Focal length

98.268 mm

Here is the ledger at the full 26° field and the patent’s aperture. One row per surface, numbered as in the diagram at the top, each with the radius it is ground to and the glass it sends light into; all values in mm:

SurfaceSISIISIIISIVSV
1 (R +17.0, into crown)+0.1676+0.1348+0.1084+0.3140+0.3398
2 (R −94.5, into air)+0.2365−0.3437+0.4996+0.0565−0.8082
3 (R −56.0, into flint)−0.2888+0.3738−0.4838−0.0910+0.7439
4 (R +15.9, into air, stop)−0.0921−0.1300−0.1835−0.3205−0.7114
5 (R +362.3, into crown)−0.0000−0.0000−0.0131+0.0147+0.7563
6 (R −77.0, into air)+0.0113−0.0281+0.0699+0.0693−0.3468
Sum+0.0345+0.0067−0.0025+0.0431−0.0264

Read down any column and you see what a triplet is for. The astigmatism column has individual surfaces contributing ±0.5 mm and a sum of −0.0025 mm: a 200-to-1 cancellation, arranged by hand in 1896. The field-curvature column shows the two crowns (surfaces 1–2 and 5–6) putting in +0.37 and +0.08 mm and the flint (surfaces 3–4) taking out −0.41 mm. That is the job of the negative element in the middle, and it is why a triplet has three pieces of glass and not two. Surface 5 contributes almost nothing to SI and SII because the marginal ray meets that nearly flat surface (R = 362 mm) close to its normal, so the refraction invariant A — n times the angle at which the ray meets the surface — is near zero there.

Want the whole ledger, editable, on the lens this page is about?

Open this lens in the Ray-Optics Designer
About the numbers themselves

The preset is the patent’s printed table at three significant figures, with the patent’s own indices (1.6114 for the crowns, 1.5679 for the flint, at the D line). Our camera-triplet guide found that a small change to the first radius alone halves the spot, so the absolute aberration values on this page belong to the printed prescription, not to whatever Taylor’s workshop actually ground. The physics of how each term scales does not depend on that. The guide begins from exactly that change, and lets you find it yourself.


Spherical aberration

One sentence: rays through the edge of the pupil focus at a different distance from rays through its center, so a point on the axis becomes a disk, and it is the only one of the five that does not vanish on axis.

What it looks like: rays from the rim of the pupil cross the axis short of where the central rays do, and the further out in the pupil a ray enters, the sooner it crosses; the figure below is that crossing, drawn from real rays. On any plane you choose the star is a disk with a bright core and a halo. In the ray fan — a plot of how far each ray in the pupil misses, which the last section of this page draws — it is a cubic through the origin, ε ∝ ρ³; in the spot diagram, where every ray is a dot, it is round, dense at the center and thin at the rim.

Where the rays cross the axis

Halve the aperture: every crossing slides toward the plane and the outermost ray comes in sevenfold.

20 %40 %60 %80 %rimlight in, on the axis(fractions of the pupil radius)focal region,enlarged
-3-2-10-200-1000100200Distance from the paraxial image plane (mm)Height (µm)the disk, edge-on: the rays fill it out to the outermost ray20 %ρ = 0.240 %ρ = 0.460 %ρ = 0.680 %ρ = 0.8rimρ = 1.0

The focal region enlarged, with the vertical scale ×2 so the crossings can be told apart; the horizontal scale is millimeters of real lens.

What is drawn on it

The heavy rule at zero is the paraxial image plane, and the bracket on it is the disk. The rim rays cross first and the innermost pair last, and nowhere do all ten meet.

Outermost ray

187.7 µm

Airy radius

4.59 µm

Rim ray crosses the axis

2.33 mm before the plane

How it scales: pupil power 3, field power 0: the two exponents, one for the aperture and one for the field, that tell the five apart. The blur grows as the cube of the aperture and does not care where in the field you look. On this lens, SI converts to a 220 µm miss for the marginal ray, third order; the real marginal ray misses by 188 µm, the number the page opened with, forty-one times the 4.6 µm floor.

Stop the lens down and the blur collapses far faster than the diffraction floor grows. Near 6 mm the geometric blur and the Airy radius are about equal, 13.0 and 11.8 µm: below that the lens is diffraction-limited on axis, above it geometry wins. Halve the patent’s aperture and the outermost ray comes in from 188 to 27 µm: the cube law says eightfold, the rays say sevenfold. The slider above stops at the patent’s aperture; our camera-triplet guide opens it further.

Opening up past the patent’s aperture

Opening from f/6.4 to f/4.4 (entrance pupil 22.22 mm, ×1.445 in diameter) multiplies SI by 1.445⁴ = 4.36 and the third-order blur by 1.445³ = 3.02, while the Airy radius shrinks only by 1.445. The slider above stops at the patent’s aperture because beyond about 16 mm the tilted bundles begin running off the edges of the glass (36 of 127 rays on axis at 22.22 mm) and the spot statistics stop being a clean aberration story; our camera-triplet guide opens it further with the apertures re-set. That guide’s remark that one stop of speed can turn a sharp lens into a soft one is this same physics — one stop being one doubling of the light the lens lets in, which is a different “stop” from the diaphragm in the middle of this lens.


Coma

One sentence: off axis, each ring of the pupil images to a circle of a different size and center, so a point becomes a comet, its tail stretched along the line between the image point and the axis.

What it looks like: take a ring of rays around the pupil and follow them to the image. Off axis they do not land in a ring around the chief ray’s image point; they land on a circle pushed to one side of it, toward the axis on this lens.

A ring of the pupil becomes a circle

Drag the azimuth once round the pupil and watch each dot go twice round its circle; drag the field and the whole comet scales with it.

The pupil

1/32/3rimtopsideone point on each named ring, at the azimuth you set

The image plane · micrometers from the chief ray

toward the edge of the field1/32/3rimthe chief rayleft and right of the pupil land here, 1×top and bottom of the pupil land here, 3×toward the axis

Rim ring, third order

-46.5 and -15.5 µm from the chief ray, toward the axis

Top and bottom of the pupil against left and right: 3 to 1.

The rays, rim ring

-170.9 and -37.4 µm from the chief ray, toward the axis

Top and bottom of the pupil against left and right: 4.57 to 1.

Where the whole pattern sits, on average

Average position, third order

-17.9 µm toward the axis

Averaged over the same 127 rays.

Average position, the rays

-46.3 µm toward the axis

Root-mean-square spot radius 142 µm.

These circles are drawn from this lens’s own third-order coma; the numbers the rays actually give are printed beside them, and the last section of the page says why they differ.

Smaller rings give smaller circles, pushed less far, all tangent to a pair of lines 60° apart: that wedge is the comet. The rays from the top and bottom of the pupil land three times as far from the chief ray as the rays from its left and right, so tangential coma is three times sagittal coma, and because the whole pattern lies to one side of the chief ray, its average position is displaced from it.

The circles in symbols

Take a ring of rays at fractional pupil radius ρ and go once round it in azimuth φ, measured from the field direction. Third-order coma sends them to ε_y = Kρ²(2 + cos 2φ), ε_x = Kρ² sin 2φ, with K = SII/(2n′u′) and SII the coma sum at that field (which is what carries the h). That is a circle of radius |K|ρ², its center displaced 2Kρ² from the chief ray, traversed twice as φ goes round once. The tangential rays (φ = 0 and 180°) land at 3Kρ² and the sagittal ones (φ = ±90°) at Kρ², which is the three-to-one rule; the tangent from the chief-ray point to any of the circles makes an angle whose sine is |K|ρ²/(2|K|ρ²) = ½, so all the circles nest inside one 60° wedge. On this lens SII is positive and n′u′ negative at every field, so K is negative and the wedge points toward the axis: K = -15.5 µm at 10° and -42.8 µm at 26°. In the tangential fan these same rays are the even part, ½[ε_y(+P) + ε_y(−P)] = 3KP², which at the rim is the 3K the ring test below compares with the real fan.

How it scales: pupil power 2, field power 1. Double the field and coma doubles, which the field slider shows; double the aperture and it quadruples. On this lens at 10° the rays’ average position sits 46.3 µm inside the chief ray, toward the axis, and the RMS spot radius, the root-mean-square distance of the same 127 rays from their average position, has grown from 119 µm on axis to 142 µm. Taylor got the third-order coma small (it moves the rim of the pupil 46 µm at 10°, against 220 µm of spherical), so the tail you see is mostly the next order up, whose own top-to-side ratio is five rather than three, and which the ledger cannot see.

The ring test

At 10° the sum is SII = 0.00242 mm against SI = 0.0345 mm, so K = 0.00242/(2 × -0.07826) = -15.5 µm and the third-order tangential coma at the pupil edge is 3K = -46 µm. Take the even part of the real tangential fan at the edge, ½[ε_y(+1) + ε_y(−1)], and the engine returns -171, not -46. Run the same ring test on the spot diagram, comparing the y-displacement of the top-and-bottom pair of edge rays against the left-and-right pair, and the ratio is 4.57, not 3. Both say the same thing: the coma-like blur that remains is mostly fifth-order linear coma, the term ∝ h ρ⁴ cos φ, whose tangential-to-sagittal ratio is 5 rather than 3. The third-order SII is small because it was corrected; what is left is the next order. The Ray-Optics Designer’s own fan label at 10° reads “Spherical-dominant, with fifth-order”, the same diagnosis from a polynomial fit.


Astigmatism

One sentence: off axis, the tangential fan of rays and the sagittal fan focus at two different distances, each to a short line, so there is no plane where the point is a point.

What it looks like: in the wavefront it is a cylinder added to the sphere, so the wavefront curves more in the tangential section than in the sagittal one, and the two sections come to a focus at two different distances from the lens.

Line, blur, line

Slide the focus from one line to the other; halfway between, the blur is round. Slide the field in and the whole thing shrinks as the square of the field.

sagittal focus

midway

tangential focus

At each line focus one fan of rays has come together and the other is still spread; midway between them the blur is round — the circle of least confusion.

Airy radius

where the slider is: −3.31 mm from the sensor

sagittal focusmidwaytangential focus← toward the lenssensor

Third order puts the two foci at

−2.91 and −3.31 mm

407 µm apart.

The rays put them at

+1.87 and −0.45 mm

2.32 mm apart, and +1.87 mm is past the right-hand end of the rule.

How long the line is, and what diffraction allows

Line half-length

31.8 µm

Blur radius halfway between the two foci, 15.9 µm.

Airy radius

4.59 µm

These three shapes are drawn from this lens’s own third-order numbers; the foci the rays actually have are printed beside them, and the last section of the page says why they differ.

At the tangential focus the point is a short line perpendicular to the field direction (the tangential rays have come together in y but are still spread in x); at the sagittal focus it is a short line the other way round; halfway between sits a round blur, the circle of least confusion, and nowhere is it a point.

The two foci in symbols

In the wavefront, astigmatism is the ½ SIII h² ρ² cos² φ term: a cylinder added to the sphere, curving the tangential section (φ = 0) more than the sagittal (φ = 90°). With the field-curvature term ¼(SIII + SIV) h² ρ² alongside it, the two sections focus at δz_T = −(3SIII + SIV)h²/(2n′u′²) and δz_S = −(SIII + SIV)h²/(2n′u′²) from the paraxial image plane, positive δz meaning downstream, further from the lens; the sums here are taken at the field, so h is already inside them. The line at either focus has half-length |SIII|/|n′u′| and the circle of least confusion, halfway between, has radius half that. The picture above is exactly these three formulas at the field you set; the foci the rays actually have, which it prints beside them, come from tracing pairs of rays a hair apart in each section and finding where each pair crosses.

How it scales: pupil power 1, field power 2, so doubling the field separates the two foci four times as far and makes each line twice as long; the field slider shows both. At 26° the rays put the foci somewhere else entirely, one of them on the far side of the sensor, and that is the story of Where third order stops.


Field curvature and the Petzval surface

One sentence: even with astigmatism removed, the surface on which a lens forms sharp images is curved, and how curved depends only on the curvatures and refractive indices of its surfaces, never on where the stop is.

What it looks like: SIV is the one sum with no ray heights and no ray slopes in it, only each surface’s curvature and the index change across it. Take the astigmatism away and both line foci coincide on one curved surface, the Petzval surface, which on this lens bows toward the glass. A flat sensor behind a curved image surface is sharp at the center and soft at the corners, or the reverse if you refocus; no focus setting fixes both.

The bowl against the sensor

Drag the field in toward the center and back out: the three surfaces are drawn out to the field you set, and the bracket is how far the sharp surface has left the sensor there.

-4-3-2-1012-40-2002040Focus shift (mm)Image height (mm)flat sensor← toward the lens
  • the Petzval surface, from the sums (not where the rays go)
  • traced rays, tangential
  • traced rays, sagittal
  • the bracket: the sag at the field you set

Petzval sag

-3.52 mm from the sensor

The rays

tangential 1.87 mm, sagittal -0.45 mm

Petzval radius

-327 mm

Curving toward the lens.

Image space from the side: the sensor is edge-on, the lens is off to the left, and the horizontal axis is stretched exactly 10 times. The heavy rule is the flat sensor, the paraxial image plane.

How deep the bowl really is

At true scale the bowl is 3.5 mm deep across a 96 mm image, and you would barely see it. The bowl itself is the third-order Petzval surface; where the traced surfaces go instead is the subject of the last section of the page.

For thin lenses in air this is the Petzval theorem: every positive element curves the field inward, toward the lens, in proportion to its power divided by its index, and the only way to flatten it with spherical surfaces is a negative element.

How it scales: pupil power 1, field power 2, exactly like astigmatism. On this lens the Petzval surface has a radius of -327 mm, curving toward the lens, a little over three focal lengths; the slider’s readout is its sag at the field you set. The traced-ray surfaces the picture also draws tell a different story at 26°, told in Where third order stops, because it is not a story about SIV.

The Petzval radius, two ways

Set the astigmatism to zero (SIII = 0) and both line foci coincide on the Petzval surface, at δz_P = −SIV h²/(2n′u′²) from the paraxial image plane, the sag the picture above brackets. Since H = n′u′ η′, with η′ the full-field image height, the Petzval sag can be rewritten as a curvature: R_P = −H²/(n′ SIV) = 1/[n′ Σ c Δ(1/n)], and for thin lenses in air 1/R_P = −Σ φ_j / n_j. On the Cooke, H at 26° is 3.7507 mm and SIV = 0.04306 mm, so R_P = −(3.7507)²/0.04306 = -326.7 mm; the surface sum Σ c Δ(1/n) = -0.003061 mm⁻¹ gives the same number. The thin-lens check agrees: element powers +0.04243, −0.04586 and +0.00963 mm⁻¹ (focal lengths 23.6, −21.8 and 103.9 mm), and Σ φ/n = 0.003061 mm⁻¹ = −1/R_P. The middle element’s power is almost as large as the front element’s and negative; that is the Petzval correction, bought at the cost of the aberrations the flint’s strongly curved surfaces add to every other column of the ledger.


Distortion

One sentence: the magnification changes with field angle, so straight lines bow, but each point is still a point; distortion moves the image, it does not blur it.

What it looks like: the distortion term of the wavefront is a pure tilt across the pupil, and a tilt does not spread a bundle, it steers it. So every point is still a point; it is just not where a perfect lens would have put it, which for a lens focused at infinity is the focal length times the tangent of the field angle, f·tan θ (the photographic convention, and the Ray-Optics Designer’s).

A square grid, magnified

The grid is drawn at ×25 to start; drag the magnifier down to ×1 and it goes square, because at true scale the corners move about one pixel.

dashed: a perfect lens’s grid, f·tan θ

pincushion: this lens, third order, at ×25

barrel: the opposite sign, not this lens

The corner, third order

168 µm outside

Its ideal height is 47.93 mm, so +0.35 %. The grid above is drawn from this, at ×25.

The corner, the traced chief ray

59 µm outside

+0.12 %, pincushion, at the 26° corner.

The grid is third order at every point; the traced chief ray’s own corner is printed beside it.

When the outer points land too far out, the corners of a square grid are pulled outward and its sides curve inward: pincushion. When they land too far in, the corners are pulled in and the sides bulge outward: barrel.

Distortion in symbols

The wavefront term ½ SV h³ ρ cos φ is a tilt across the pupil, and a tilt steers the whole bundle without spreading it. The chief ray lands at Δy = SV h³/(2n′u′) away from the ideal image height f·tan θ, the photographic convention for a lens focused at infinity (the sums here are taken at the field, so h is inside them). As a fraction of the image height that is SV/(2n′u′ f tan θ), which grows as h². Positive Δy at the edge is pincushion, negative is barrel. On this lens at 26° third order gives 168 µm and the traced chief ray 59 µm; the grid above is drawn from the third-order term at every point of the grid, and the magnifier multiplies that shift by the number you set, nothing else. The traced ray gives less because the next order pulls the other way past about 15°: its distortion peaks near 23° and eases toward the corner, which is why a grid drawn from it would inflate rather than bow.

How it scales: pupil power 0, field power 3. Stopping down does nothing to it, and as a fraction of the image height it grows as the square of the field. One part in eight hundred at the corner of the plate is invisible to the eye, which is why the picture above has a magnifier, and it was one of the triplet’s selling points.

Notice that the per-surface SV entries in the ledger were the largest numbers in the table (±0.8 mm) and the sum is 30 times smaller. Distortion is not corrected surface by surface; it is corrected by balancing large contributions of opposite sign, which is what the stop position controls, and that is the next section.


Moving the stop

Here is the single most useful theorem in the whole subject. Move the aperture stop along the axis and the marginal ray does not change at all, so SI is unchanged. SIV knows only curvatures and indices, so it is unchanged too. Only the chief ray moves, and it moves by a fixed multiple of the marginal ray; the multiple is ε, the stop-shift parameter, the change in chief-ray height at the first surface divided by the marginal height there (the substitution is in How the sums are computed, above). Write S* for a sum after the move; the stop-shift equations are one line per sum:

  • SI* = SI
  • SII* = SII + ε SI
  • SIII* = SIII + 2ε SII + ε² SI
  • SIV* = SIV
  • SV* = SV + ε (3 SIII + SIV) + 3ε² SII + ε³ SI

On the widget, look at the bars first, SI and SIV pinned while the other three slide along their curves; then the field surfaces; then the distortion.

Move the stop

Slide the stop back to the tick at 2.27 mm, where third-order coma passes through zero.

The five sums as the stop moves (mm, 26° field)

0246810-0.0500.05Stop position behind the flint (mm)Seidel sum (mm)
SISpherical0.03451pinned
SIIComa0.006704
SIIIAstigmatism-0.002490
SIVField curvature0.04306pinned
SVDistortion-0.02636

Curves: the stop-shift equations from the patent row. Dots: the engine’s own Seidel table for the shifted system. Spherical aberration and field curvature do not move at all, and neither does the focal length (98.268 mm).

Real-ray field surfaces (mm from the image plane)

0510152025-202Field angle (deg)Focus (mm)

Solid: tangential. Dashed accent: sagittal. Gray dashed: the third-order Petzval surface, which the stop cannot move.

Distortion (% of the ideal image height)

051015202500.050.10.15Field angle (deg)Distortion (%)

Spot at 13°

147.2 µm RMS

Stop-shift parameter

0.00000

Third-order coma zero

2.2657 mm

ε = -0.19426: the patent row's coma sum over its spherical aberration sum, negated.

Why 13° and not 26°

Where third order is still a fair guide. At 26° the tilted bundle starts running off the edges of the glass from about 5 mm, so that spot is not one quantity across this slider and is not shown.

Coma moves linearly with the stop, astigmatism quadratically, distortion cubically, and the two aperture-only terms not at all. So if a lens has spherical aberration, there is always a stop position that removes its third-order coma: ε = −SII/SI. This is exactly what the position of the diaphragm is for in a triplet, and why Taylor’s patent specifies it.

On the Cooke, the patent puts the stop “as closely as possible behind” the middle element, which the saved lens this page uses puts at the rear surface of the flint. Slide it backward through the 11.2 mm airspace toward the third element, and the engine reports (26° field, patent aperture, sums in mm):

Stop behind flintεSISIISIIISIVSV
0 mm (patent)0.0000.0345+0.00670−0.002490.0431−0.0264
2.25 mm-0.1930.0345+0.00005−0.003790.0431−0.0327
11.2 mm-0.9820.0345−0.02719+0.017640.0431−0.0746

SI and SIV do not move, and neither does the focal length (98.268 mm at every position). Third-order coma passes through zero at 2.27 mm behind the flint. Astigmatism dips and then climbs; distortion grows monotonically, because the ε(3SIII + SIV) term is dominated by SIV and never changes sign.

Checking the table

Every entry is the engine’s Seidel table for the shifted stop, and every entry matches the stop-shift equations evaluated from the patent row to the last printed digit: over the whole airspace the largest disagreement between the two, in any of the five sums, is 7.41e-15 mm. At the coma zero, ε = −0.00670/0.0345 = −0.194, which is the ε = −SII/SI of the paragraph below evaluated at the patent row.

The real rays follow, roughly. At 13°, the spot RMS goes from 147 µm at the patent’s stop to 121 µm at the third-order coma zero, keeps improving to 106 µm at 5 mm, then degrades to 277 µm at the far end of the airspace. The ledger tells you the direction and roughly where; the rays tell you where.


Where third order stops

The Seidel sums are the first term of a series, and this lens is a good place to watch the series stop converging. Here are the real-ray field surfaces beside the third-order ones, in mm from the paraxial image plane, negative meaning inside, toward the lens:

FieldTangential, realTangential, third orderSagittal, realSagittal, third order
5°−0.088−0.093−0.103−0.107
10°−0.283−0.380−0.382−0.433
26°+1.870−2.905−0.451−3.312

At 5° they agree to within 6.4 %. At 10° third order overshoots the tangential focus by a quarter and the sagittal by an eighth, and every “third order says, the rays say” clause in the sections above is this row. By 26° the real tangential surface has swung to the far side of the paraxial plane while third order says it should be 2.9 mm inside: the two have opposite signs. That is fifth-order astigmatism and oblique spherical aberration, and it is not a defect of the lens: it is how the lens works. Taylor did not have the Seidel sums at zero and the field flat; he had them small and the higher orders bending the field surfaces back so that the whole 26° plate was tolerably sharp. The real tangential surface recrosses the image plane near 20°, and the spot RMS across the field falls with it, from 142 µm at 10° to 114 at 20°.

Where the blur is smallest

The four samples the field ladder gives: 119 µm on axis, 142 at 10°, 114 at 20° and 183 at 26°. The fall does not bottom out at the 20° sample: swept finely, the minimum is 94 µm near 22.6°, below both neighbors, which is the fifth-order story this section is telling — the higher orders bend the field surfaces back and the sharpest zone of the plate is out near the corner, not at the center.

Here are the fans and the spot of the whole lens at once, at 26° with the third-order prediction switched on.

Ray fans and spot diagram

Solid is what the rays do, dashed is what the ledger predicts; the gap is everything third order cannot see. Slide the field toward 0° and watch it shrink.

Tangential fan · εy vs Py

-1-0.500.51-5000500Fractional pupil coordinateε (µm)

Coma-dominant, with fifth-order

Sagittal fan · εx vs Px

-1-0.500.51-5000500Fractional pupil coordinateε (µm)

Spherical-dominant, with fifth-order

Spot diagram · 127 rays

One box for the whole slider range, 1081 µm across, plotted about the rays’ average position. Dashed circle: the Airy radius. Cross: the chief ray.

RMS radius

182.7 µm

Outermost ray

514.7 µm

Airy radius

4.59 µm

Dashed: what the Seidel sums predict. Where the dashed curve leaves the real fan, that gap is the higher-order aberration; see Where third order stops.

Two things to read off the fans: the two leave the origin at visibly different slopes, the tangential climbing and the sagittal dipping, and that difference is astigmatism, 5.7 times what third order drew, though with the same sign; that the tangential fan climbs at all, where its dashed twin dips, is the field-curvature sign flip of the table above. And the tangential fan is no longer odd about the origin, the top and bottom rays missing the chief ray on the same side, which is coma.

The Seidel view in the Ray-Optics Designer is the third-order ledger and says so; the spot diagrams and ray fans are exact real rays through the exact surfaces, and where they disagree with the ledger, the rays are right. The fan labels under each plot are fitted to the real fan and name what dominates: at 26° on this lens “Coma-dominant, with fifth-order” on the tangential side, which is why the comet in the coma figure, dragged to the full field, points the way it does and the drawn circles do not. Learn to read the sums for what to change and the rays for whether it worked.


What this page left out: color

Everything above was at one wavelength, because the patent prints only one index per glass. Real glass has an index that falls with wavelength, so a single lens has a different focal length in blue than in red (axial color) and a different magnification (lateral color), and both are as large as anything in the Seidel table. Correcting them is the reason the Cooke’s two outer elements are crown glass and the middle one is flint, and the reason the Abbe number is the first thing a designer reads off a glass map. That is the next tutorial: chromatic aberration and the achromatic doublet.


Try it on the tool

Open the Ray-Optics Designer, load the Cooke Triplet f/6.5 preset from the preset list, and look at the Seidel aberrations view, then Ray fans and Spot diagrams at 0°, 10° and 26°. Then do the three things this page did: halve the entrance pupil diameter and watch the on-axis spot fall about sevenfold; insert a plane surface after surface 4 (the back of the middle element — see the diagram at the top) and move the stop to it; and nudge the first radius, which is where the camera-triplet guide begins. Loading this lens, editing it and seeing it drawn need no account. The Seidel table, spot diagrams and ray fans need a free one, and so does the seventh surface the stop shift adds, because without an account the editor holds six; nothing here needs more than that.


References

  • W. T. Welford, Aberrations of Optical Systems (Adam Hilger, 1986), chapters 7 and 8: the Seidel sums in the A, Ā, δ(u/n) form used here, the stop-shift equations, and the wavefront expansion in this normalization.
  • R. Kingslake and R. B. Johnson, Lens Design Fundamentals, 2nd ed. (Academic Press, 2010), chapters 6 to 11 and 14: the aberration-by-aberration treatment, the Petzval theorem and the triplet.
  • W. J. Smith, Modern Optical Engineering, 4th ed. (McGraw-Hill, 2008), chapter 3: the transverse and longitudinal pictures of each aberration and the ray-fan reading.
  • M. Born and E. Wolf, Principles of Optics, 7th ed. (Cambridge, 1999), chapter 5: the symmetry argument for the five primary aberrations and the wavefront-to-ray relation.
  • H. H. Hopkins, Wave Theory of Aberrations (Oxford, 1950): the wavefront formulation and the insensitivity of W to the reference-sphere radius.
  • H. D. Taylor, “Lens”, US Patent 568,052 (1896), Series III prescription, the lens on this page. The preset’s provenance, including the convexity-to-directional radius-sign conversion, is recorded in the Ray-Optics Designer manual.