CODEV Ch8
CODEV Ch8
CODEV Ch8
Search
Chapter 7
Non-Spherical Surfaces
Many optical systems use non-spherical surfaces, including aspheres, cylinders, and
others. CODE V has a number of types available, including diffractive surfaces.
Contents
Non-Spherical Surfaces
Non-spherical surfaces are widely used in modern optical systems. Conic sections
appear in many reflecting systems (including the common Newtonian telescope
with its parabolic primary mirror). Polynomial aspheres are becoming quite
widespread due to modern manufacturing methods such as glass or plastic molding
and numerically controlled methods such as diamond turning (modern diamond-
turning methods can directly produce precision aspherics on metal, glass, plastic,
and special materials such as germanium). Many infrared systems use aspherics,
and most optical disc systems use a doubly aspheric singlet for the disc-reading
objective lens. Diffractive optics have also become more common. Ruled
diffraction gratings have been used for many years in spectrometers and other
instruments. Holographic elements are used in some applications. Non-holographic
diffractives such binary optics are becoming more common.
CODE V has two classes of surfaces, rotationally symmetric and non-rotationally
symmetric. These are listed in the section below. In general, when using rotationally
symmetric surfaces, you can use Y-fields and circular apertures (maintaining
bilateral symmetry), whereas for non-rotationally symmetric surfaces you may need
to define multiple field or object points with X- and Y components, and use
rectangular or elliptical apertures.
Non-spherical surfaces are all identified by descriptive names in the user interface
and by a three-letter command code (SPH, CON, etc.), used in command inputs and
in many output listings. In this section, we will introduce several of the more
common surface types. Other types are defined and modified in exactly the same
way, although they will typically have more and different coefficients and
parameters defining them. Do a search in the Help or in the CODE V Reference
Manual for complete information on any special surface type (the three-letter code
is a good search string, although you can use other key words as well).
A Non-Spherical Example
Search
Flat-Field Schmidt
As an example of a system with a non-spherical surface, take a look at the CODE V
sample lens Flat-field Schmidt. This is a well-known catadioptric (reflecting and
refracting) design form. It has an aspheric glass element of very low optical power
in order to correct the spherical aberration of the spherical primary mirror, allowing
a fairly fast system (f/1.3 in this case). Field curvature is corrected by a lens element
very close to the image plane. We won’t do much with this example except look at
the user interface elements that you will use in working with all non-spherical
surface types.
1. Choose the File > New menu to launch the New Lens Wizard. Click the Next
button to go to the second page in the New Lens Wizard.
2. Select CODE V Sample Lens (default) and click Next.
3. Scroll down to find the file cv_lens:schmidt.len, and click on it.
4. Click the Next button on the remaining New Lens Wizard pages to view and
accept the system data values for the lens, or click the Finish button to accept
all values and view the lens.
5. Choose the Display > View Lens menu. In the View Lens dialog box, click
Hatch back of mirrors on the Lens Drawing tab, then click OK.
The lens picture is shown below (the label was added with text and line annotation
tools that are available in all graphics windows).
Here is the LDM spreadsheet for this lens. Note the word Asphere in the Sur-
face Type column for surface 1 (stop).
Search
If you double-click the word Asphere on surface 1, you will see the following
list.
This is one way to change a surface to a non-spherical type. It is also handy for
changing to a sphere, which requires no additional data beyond the always-pres-
ent radius or curvature. But other surface types require additional data found on
the surface properties dialog.
Tip: Note that changing any special surface to a sphere throws away all the
special coefficient data (but you can always use Undo to get it back if you make
a mistake).
This will open the surface properties dialog at the appropriate page, as shown
below.
Search
On the Surface Type page, you can change the type of surface from the drop-
down list, and you can also change the values or the status (variable, frozen,
pickup, zoom) of any of the coefficients associated with the surface type (right-
click on the value of interest). These are shown in the spreadsheet to the right.
If you choose the Display > List Lens Data > Surface Data menu, a window is
displayed containing a text-style listing of all the surface information. The infor-
mation is listed in compact form labeled with the command names of the vari-
ous parameters. These command names are usually listed in the Surface
Properties window as well, such as Conic Constant (K). The command names
are usually derived from the variable name used in equations for the various sur-
face types, some of which you will see later in this chapter.
Search
How important is the aspheric surface corrector for this system? Take a look at a
spot diagram for the lens with and without the asphere.
7. Choose the Analysis > Geometrical > Spot Diagram menu. In the Spot
Diagram dialog box, click OK to run this option.
8. Zoom in on the resulting output with the magnifier icon and note the scale size
(it should be 0.05 mm).
9. Restore the scale (click the 1:1 button) and tear off a copy of the plot for later
comparison.
10. In the LDM spreadsheet window, double-click Asphere on s1 and select
Sphere from the dropdown list.
11. Click the Modify Settings button on the Spot Diagram output window you just
created. Click Aberration Scaling, enter 0.05 for the scale size, and click OK.
The Spot Diagram output window is regenerated.
Search
The before and after Spot Diagram plots are shown above. With the asphere, the on-
axis (zero field) RMS spot diameter is 0.0049 mm, while without the aspheric
corrector, the RMS spot diameter is 0.57 mm (shown on the Text tab of the Spot
Diagram output window).
It’s clear that a spherical mirror is not a good solution. Even with the corrector, this
lens is far from diffraction limited. Undo the sphere change and run the supplied
macro SPOTDET.SEQ (Tools > Macro Manager, in the Geometrical Analysis
category) and enter .05 for the scale and AIRY for the detector X width—you will
have to zoom in very close to see the circle for the f/1.3 Airy disk size.
Chromatic aberration is a problem for this lens due to the two glass elements
(note the separation of the red and blue light, especially evident in the classical
coma-shaped off-axis spot). Lack of such chromatic aberration problems is one
advantage of all-reflecting designs.
A Note on Obscurations
Search You may have noticed that there is no central obscuration anywhere in this system.
Strictly speaking, this is a mistake, because the detector and field lens would block
the central 33% radial zone of the entrance pupil. This is not important for
explaining the user interface for non-spherical surfaces, but it would be important
for detailed analysis, especially for MTF and other diffraction-based analysis.
For the sake of completeness, here are the commands (used for brevity) to enter a
dummy surface to hold the required obscuration. Database items from CODE V are
used to show the logic of the calculation, but you could also calculate or estimate
the values by hand and insert the numbers in the spreadsheets. THI is thickness,
OAL is overall length, and SD is semi-diameter (default aperture radius). You can
type these commands at the CODE V prompt in the command window.
Conic Surfaces
Search
The conic surface is used to specify pure conics, such as parabolas, ellipses, and
hyperbolas. True conics are characterized by two foci; any ray passing through one
of the foci will pass perfectly through the other focus with no aberration. Conics are
defined by their conic constant k. The equation for a conic is given by the equation
shown below.
c( r2 )
z = ------------------------------------------------------ where r2 = x2 + y2 and c = vertex curvature
1 + 1 – ( k + 1 )c 2 r 2
k=0 sphere
-1 <k <0 ellipse
k = -1 parabola
k < -1 hyperbola
k>1 oblate sphere (not a true conic)
Ellipses
An ellipse, shown in the following figure, can also be defined by its major and
minor axes 2a and 2b. It is defined as the locus of points whose sum of the distances
from two points (foci) is constant. With the local coordinate system at the surface
vertex, the equation is given by
z – a ) 2- + ----
(----------------- r 2- = 1
a2 b2
The vertex radius of curvature is given by R = b2/a and the conic constant is given
by k = (b2 - a2)/a2. The distance between the two foci is 2F, where F2 = a2 - b2, and
the distance from the surface vertex to the foci is given by a ± F.
Foci
Surface vertex a
(radius R = 1/C)
Hyperbolas
Search A hyperbola, shown in the following figure, also has major and minor axes a and b,
and is defined as the locus of points whose difference of the distances from two
points (foci) is constant. The equation can be given as
( z – a )2 r2
------------------ + ----- = 1
a2 b2
The radius of curvature is R = b2/a and the conic constant k = -(a2 + b2)/a2. The
distance from the center of the hyperbola to the foci is c = (a2 + b2)1/2 with the
distance from the surface vertex to the foci given by F = -a ± c. The hyperbola
curves extend to be asymptotic to lines at ±θ from the Z axis. The conic constant k
can also be given by k = -(1+tan2θ). As such, a cone can be approximated in
CODE V by a hyperbola with an appropriate conic constant and a very small radius
(e.g., 0.000001).
b
Z
Focus
Surface vertex
(radius R = 1/C)
Parabolas
A parabola is a degenerate ellipse or hyperbola with one of the foci being at
infinity. It is defined as the locus of points whose distances from a point (the focus)
and a line are the same. Since the conic constant is k = -1, it is described by an exact
second order equation, given by the following, where R is the vertex radius of
curvature and F is the focal length (F = R/2).
r 2- = -----
r 2-
z = ------
2R 4F
Entering an Ellipse
Search To show how these equations are used, you can enter an ellipsoidal reflector with a
major axis of 100 (a = 50) and a minor axis of 60 (b = 30). Using the equations from
the previous section, you get the vertex radius R = 18 and the conic constant k = -0.64.
The distance from the surface vertex to the foci is 50 ± 40 = 90, 10.
Place the object at the distant focus (object distance = 90) and the image at the closer
focus (image distance = 10). Since you are using the “right hand side” of the ellipse,
the vertex radius will be negative (its center of curvature is to the left of the surface).
Since the light reflects, the image distance is negative also.
1. Choose the File > New menu to launch the New Lens Wizard, and click Next to
bypass the welcome screen.
2. Click the Blank Lens button and click Next.
3. On the pupil screen, accept the default for now (EPD 2 – we will later change
this to object numerical aperture, but this is not available now because the
default object distance is infinite). Click Next.
4. Click Next on the next four screens, accepting the default wavelength and
object angle (zero), and click Done on the final screen.
5. In the LDM spreadsheet window, change the object thickness to 90.0, the
Y radius of surface 1 (stop) to –18.0, the thickness of surface 1 to –10, and the
refract mode for s1 to Reflect (double-click the word Refract).
6. Choose the Lens > System Data menu to display the System Data window. On
the Pupil page, change the pupil specification to Object Numerical Aperture,
and the value to 0.4, then click the Commit Changes button.
7. Go to the System Settings page in the System Data window and change the title
to something like Ellipse.
8. Choose the Display > View Lens menu and on the Plot Parameters tab,
change the start surface in the Surface Span spreadsheet to Object, then click
the Ray Definition tab.
9. Click in the first cell of the Fan of Rays spreadsheet to get the default ray fan,
then click OK.
Search
The following View Lens output window is displayed.
What’s wrong with the picture? Only some of the inner zone rays are drawn.
The default spherical reflector can’t trace the marginal rays for NAO 0.4—you
need to make the surface a conic. Leave the View window open so you can
rerun it in a moment.
10. In the LDM spreadsheet, right-click on the surface type for surface 1 (stop) and
choose Surface Properties from the shortcut menu.
11. On the Surface Type page in the Surface Properties window, choose Conic
from the dropdown list labeled Type, then click another cell to commit this
change.
12. Change the value for the parameter Conic Constant (K) to –0.64 and tab or
click elsewhere to commit this change.
Search
Now the picture is what you might expect. All of the rays from the left focus
(object point) are perfectly focused to the right (close to vertex) focus. Note that
this perfect image will degrade pretty quickly if you move the object point off
the axis. Try changing the object height to 1 mm and see what happens to a
quick spot diagram or ray aberration plot.
Polynomial Aspheres
Search
The polynomial asphere is a rotationally symmetric surface whose surface profile
(sag) is given by the following equation.
cr 2 - + Ar 4 + Br 6 + Cr 8 + Dr 10 + Er 12 + Fr 14 + Gr 16 + Hr 18 + Jr 20
z = -------------------------------------------------
1 + 1 – ( 1 + k )c r 2 2
where:
z(r) 2 2
r is the radial distance = x + y
c is the curvature at the pole of the
r surface.
k is the conic constant.
z A through J are the 4th through 20th
order deformation terms.
It is a 20th-order deformation to a base conic (starting at the 4th order; the 2nd
order term is essentially equivalent to the curvature). If the aspheric coefficients
A-J are all zero, the surface reduces to a pure conic. If the conic constant k=0, the
surface is a 20th order deformation to a sphere.
Since the equation of a parabola is a pure 2nd order equation, if the conic constant
for a polynomial asphere is -1, the equation is an exact 20th order polynomial
(technically, spheres and general conics are of infinite order).
Fresnel Surface
The ASP type surface can also be used to represent a Fresnel surface (with
infinitesimally narrow zones). Such a Fresnel surface is specified by making the
Fresnel base curvature (CUF) non-zero (for a flat substrate, use a very small
curvature, such as 0.000001). The surface sag then follows the base substrate, but
the surface slope is defined by the derivative of the equation above.
In the list of surface types, you will find this listed as the type Thin Fresnel. This
entry really creates an ASP type surface with access to the Fresnel base curvature.
You will also find other types of Fresnel surfaces in list of surface types. The type
Fresnel Planar Substrate (SPS FRS) is the finite-facet-size version of the basic
(ASP-based) Fresnel. Search in the CODE V Reference Manual for SPS FRS or
Search
Fresnel planar to find more information on these surfaces.
Entering an Asphere
This procedure applies to every type of non-spherical surface.
1. In the LDM spreadsheet window, right-click on the Surface Type cell of the
surface you wish to change and choose Surface Properties from the shortcut
menu. This brings up the correct page of the Surface Properties window.
2. Choose Asphere in the Type field, and click in another field to commit this
change (the correct coefficients for this type will now appear).
3. Enter the appropriate special surface data (if known) in the labeled data fields.
Most often you will not know a value but will use optimization to vary it.
4. To make a parameter variable (such as the conic constant or an aspheric
coefficient), right-click on the parameter and choose Vary from the shortcut
menu.
Tip: Aspheric “sag” is the deviation from a sphere caused by the aspheric
deformation. To display tables of sag-related data for any non-spherical surfaces
in the lens, choose the Analysis > Fabrication Support > Fabrication Data
Tables menu and click the Sag Data tab. Select the Display sag table checkbox
and any special options.
It’s a good idea to check with the fabricators before designing a system with
aspheres or other non-spherical surface profiles. Molding, diamond turning, and
other methods have limitations in achievable sag, element size, aperture shape,
part thickness, materials, etc. You can constrain these things in CODE V
optimization if you are aware of the requirements.
Now you can analyze this lens as a starting point. Use the Quick Lens Drawing
and Quick Ray Aberration Plot buttons on the toolbar, shown in the following
Search
illustrations.
The spherical aberration is quite large. Note the 2 mm scale on the ray aberra-
tion plot.
Now you can make the front surface aspheric, and allow the fourth and sixth
order coefficients to be variables for optimization.
8. In the LDM spreadsheet window, right-click the Surface Type for surface 1
(stop) and choose Surface Properties from the shortcut menu.
9. In the dropdown Type list, choose Asphere, then click in another cell.
Search 10. To make the coefficients variable, right-click in the 4th Order Coefficient field
and choose Vary Parameter (leave the value at the default, 0.0). Do the same
with the 6th Order Coefficient field.
11. Choose the Optimization > Automatic Design menu. Since you are only
varying aspherics and defocus, the gross properties of the lens such as focal
length cannot change, and you won't need any specific constraints. Click OK to
start this default AUTO run. Look at the Surface Properties window again.
12. Use the Quick toolbar buttons again to make a new lens picture and ray
aberration curve and compare them to the starting versions. Note that the
automatic scale factor on the Quick Ray Aberration Plot is now 0.0045 mm
(was 2 mm), and the curve has multiple crossing points and residual higher
order aberration.
Search
If you return to the Surface Properties window and make the 8th-order coefficient
variable, then re-run the optimization window (Optimization > Automatic Design
> OK) and aberration curve, the scale factor will be 0.0005 mm and the lens will be
essentially perfect (on-axis only, one wavelength only, SA corrected to 8th order).
13. Choose the File > Save Lens As menu to save your optimized aspheric singlet
in a file for possible later use.
Tip: When optimizing an asphere, you should vary either the 4th-order
coefficient (A) or the conic constant (K), but not both. They are often roughly
equivalent in optical effect, but if the vertex curvature of the surface is very weak
(as in the Schmidt corrector), the conic will not be an effective variable, and you
should use the 4th-order term.
Diffractive Surfaces
Search
Color-Corrected IR Lens
Diffractive surfaces have a number of special features and therefore they have their
own page in the Surface Properties window. If you right-click on the surface type
for s2 (Sphere) and choose Surface Properties, you will find only the Y Radius
listed in the Surface Properties window. Click the Diffractive Properties item in
the navigation tree on the left side of the window and you will see the following
page:
Search
In this case, construction points have been specified as if the surface were actually
holographic, although the two points are located at the same coordinates and thus
contribute no power to the diffractive. You could change the Diffractive Surface
Type to Phase Polynomial (Kinoform/Binary) and it will have no effect in this case
(this actually zeroes out the construction point coordinates).
Note that the polynomial type is rotationally symmetric, and that there are phase
terms from R**2 to R**20. These phase terms can be varied in optimization, much
like aspheric coefficients (it is best to introduce lower order terms first, adding
higher orders in later optimization runs).
The output is a plot showing the ray distribution. The plot will give you an idea of
how many rays the program will use, as well as their positions. Smaller values
Search
(denser grids) mean more rays. If you use high order aspheres or diffractives, you
should consider using a smaller DEL (e.g., 0.15 which gives 70 rays), and also run
the AUTOGRID macro to visualize the grid that AUTO will use.
Search