Tutorial 3b - 1st & 3rd Order Aberrations
This tutorial illustrates how various aberration coefficients are computed.
[1]:
from optiland.samples.objectives import TripletTelescopeObjective
[2]:
lens = TripletTelescopeObjective()
lens.draw()
[2]:
(<Figure size 1000x400 with 1 Axes>, <Axes: xlabel='Z [mm]', ylabel='Y [mm]'>)
[3]:
print("Seidel Aberrations:")
for k, seidel in enumerate(lens.aberrations.seidels()):
print(f"\tS{k + 1}: {seidel:.3e}")
Seidel Aberrations:
S1: -2.707e-03
S2: -1.412e-03
S3: -1.479e-03
S4: -5.568e-04
S5: 3.606e-05
[4]:
print("Third-order transverse spherical aberration:")
for k, value in enumerate(lens.aberrations.TSC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order transverse spherical aberration:
Surface 1: -5.087e-01
Surface 2: -2.442e-01
Surface 3: 2.466e-02
Surface 4: -2.759e+00
Surface 5: 3.481e+00
Surface 6: -5.623e-04
[5]:
print("Third-order longitudinal spherical aberration:")
for k, value in enumerate(lens.aberrations.SC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order longitudinal spherical aberration:
Surface 1: -2.849e+00
Surface 2: -1.367e+00
Surface 3: 1.381e-01
Surface 4: -1.545e+01
Surface 5: 1.949e+01
Surface 6: -3.149e-03
[6]:
print("Third-order sagittal coma:")
for k, value in enumerate(lens.aberrations.CC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order sagittal coma:
Surface 1: -2.491e-02
Surface 2: 2.011e-02
Surface 3: 4.053e-03
Surface 4: 8.930e-02
Surface 5: -9.345e-02
Surface 6: 9.330e-04
[7]:
print("Third-order tangential coma:")
for k, value in enumerate(lens.aberrations.TCC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order tangential coma:
Surface 1: -7.473e-02
Surface 2: 6.034e-02
Surface 3: 1.216e-02
Surface 4: 2.679e-01
Surface 5: -2.803e-01
Surface 6: 2.799e-03
[8]:
print("Third-order transverse astigmatism:")
for k, value in enumerate(lens.aberrations.TAC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order transverse astigmatism:
Surface 1: -1.220e-03
Surface 2: -1.657e-03
Surface 3: 6.659e-04
Surface 4: -2.890e-03
Surface 5: 2.509e-03
Surface 6: -1.548e-03
[9]:
print("Third-order longitudinal astigmatism:")
for k, value in enumerate(lens.aberrations.AC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order longitudinal astigmatism:
Surface 1: -6.831e-03
Surface 2: -9.280e-03
Surface 3: 3.729e-03
Surface 4: -1.618e-02
Surface 5: 1.405e-02
Surface 6: -8.670e-03
[10]:
print("Third-order transverse Petzval sum:")
for k, value in enumerate(lens.aberrations.TPC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order transverse Petzval sum:
Surface 1: -1.850e-03
Surface 2: -9.425e-05
Surface 3: -1.636e-03
Surface 4: -5.416e-04
Surface 5: 1.167e-03
Surface 6: 1.395e-03
[11]:
print("Third-order longitudinal Petzval sum:")
for k, value in enumerate(lens.aberrations.PC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order longitudinal Petzval sum:
Surface 1: -1.036e-02
Surface 2: -5.278e-04
Surface 3: -9.159e-03
Surface 4: -3.033e-03
Surface 5: 6.537e-03
Surface 6: 7.812e-03
[12]:
print("Third-order distortion:")
for k, value in enumerate(lens.aberrations.DC()):
print(f"\tSurface {k + 1}: {value:.3e}")
Third-order distortion:
Surface 1: -1.503e-04
Surface 2: 1.443e-04
Surface 3: -1.593e-04
Surface 4: 1.111e-04
Surface 5: -9.870e-05
Surface 6: 2.540e-04
[13]:
print("First-order transverse axial color:")
for k, value in enumerate(lens.aberrations.TAchC()):
print(f"\tSurface {k + 1}: {value:.3e}")
First-order transverse axial color:
Surface 1: -1.893e-01
Surface 2: -1.120e-01
Surface 3: -5.758e-02
Surface 4: -2.546e-01
Surface 5: 7.011e-01
Surface 6: -1.541e-02
[14]:
print("First-order longitudinal axial color:")
for k, value in enumerate(lens.aberrations.LchC()):
print(f"\tSurface {k + 1}: {value:.3e}")
First-order longitudinal axial color:
Surface 1: -1.060e+00
Surface 2: -6.274e-01
Surface 3: -3.224e-01
Surface 4: -1.426e+00
Surface 5: 3.926e+00
Surface 6: -8.627e-02
[15]:
print("First-order lateral color:")
for k, value in enumerate(lens.aberrations.TchC()):
print(f"\tSurface {k + 1}: {value:.3e}")
First-order lateral color:
Surface 1: -9.272e-03
Surface 2: 9.230e-03
Surface 3: -9.461e-03
Surface 4: 8.240e-03
Surface 5: -1.882e-02
Surface 6: 2.556e-02