Abstract
We obtain simple exact formulas for the refracted wavefronts through plano-convex aspheric lenses with arbitrary aspheric terms by considering an incident plane wavefront propagating along the optical axis. We provide formulas for the zero-distance phase front using the Huygens’ Principle and the Malus-Dupin theorem. Using the fact that they are equivalent, we have in the second method found a way to use an improper integral, instead of the usual evaluated integral, to arrive at these formulas. As expected, when the condition of total internal reflection is satisfied, there is no contribution to the formation of the refracted wavefront.
© 2013 Optical Society of America
1. Introduction
The aspheric lens can help simplify optical system design by minimizing the number of elements required and yields sharper images than conventional lenses. Aspherical elements are particularly useful for correcting distortion in wide-angle lenses. In summary aspheric optical surfaces deliver higher performing, more compact, and lighter systems in a wide range of applications. There are three ways to make an aspherical lens. The basic way is to grind a piece of glass down to the right shape. This is quite difficult due to the extreme precision required to achieve the complex geometry. Another way is to mould a glass lens element. Finally one can cement a plastic resin aspheric surface to the top of a glass spherical element [1].
The caustic can be defined as the locus of the principal centers of curvature of a wavefront [2]. In other words, if equal optical paths are measured along each ray from the source, the surface constructed by the end points will be normal to all the rays. These surfaces are the phase fronts of the wave system, for which the rays are defined within the geometrical optics approximation. Although the caustics and wavefronts either by reflection or refraction are part of a well known subject [3–5], the contribution in this work it is to provide simple formulas of the wavefront surfaces caused by refraction on plano-convex aspheric lenses exclusively in a meridional plane by using the Huygens’ principle and the theorem of Malus-Dupin. In this manuscript we consider the aspheric equation given according to [5], although recently there have been defined new formulas to represent this kind of surface [6, 7].
The treatment in this work is physically equivalent to that presented by Shealy and Hoffnagle [3, 4] in their solution to the wavefront of an aspheric lens, and the results presented here illustrate the equivalence of the Huygens’ principle and the eikonal treatment. On a practical level, this Huygens’ approach has advantages: since the treatment here involves performing a sum of wavelets whose envelope will form the wavefront and does not imply solving a differential equation, we do not need to construct initial boundary conditions. The final result is a relatively simple solution that can be used to directly compute the wavefront produced by an arbitrary aspheric lens, and vice versa, starting with the aspheric coefficients used in practical applications. From an optical design point of view, this is applicable to aspheric lens design and more generally to the design of arbitrary lenses for waveform shaping.
2. Huygens’ Principle and its phase fronts
Throughout this manuscript lower cases for (z, y) will be used to designate either formulas for the caustic surfaces or refracted rays through lenses and upper cases for (Z, Y) will be used to designate the geometrical wavefront. So, we define that the Z axis is parallel to the optical axis, we assume that the Y − Z-plane is the plane of incidence, which is a cross section of a plano-convex aspheric lens with an arbitrary number of aspheric coefficients, whose paraxial radius is R, and the origin of the system 𝒪 is placed at the vertex of the plano-convex lens. We assume that a plane wave is incident on the lens parallel to the optical axis, crossing the plane face of the lens without being deflected, and this is propagated to the aspheric surface. In this way the Huygens’ principle assumes that a wavefront is the envelope of an aggregate of wavelets centered on a previous wave front in the wave-front train. The idea here is to find the family of wavelets that are centered on the aspheric refracting surface and whose envelope yields the refracted wavefront, which also is called the phase front, zero-distance phase front or alternatively the archetype wavefront. We follow a procedure similar to [8], in such a way that it is easy to see that the wavelets are formed by a disturbance of circles centered along the aspheric surface with radii varying as are shown in Fig. 1(a), and they are represented by
where H is the entrance aperture, t is the thickness of the lens, ni the index of refraction of the lens for a predefined wavelength which is immersed in a medium with index of refraction na, (ni > na), and where we have assumed that Sh represents the aspheric equation given by where c = 1/R is the paraxial curvature, k is the conic constant, A4, A6, A8, ... A2N, are the aspheric order terms with N arbitrary and h represents either the height for an arbitrary incident ray or the height for the centers of the wavelets along Y-axis lying on the meridional plane. It is worth commenting that we are using the sign convention suggested in [9]. In order to obtain the wavefront, which is the envelope of the wavelets represented by Eq. (1), we differentiate Eq. (1) with respect to h and reducing further we have where S′h is the first derivative with respect to h, thus from Eq. (2) we obtain Finally, solving for Z and Y by using Eqs. (1) and (3) we get where the subscript (+) sign provides the equation of a wavefront progressing in a retrograde direction and is ignored, and for (−) sign it yields a wavefront progressing in the forward direction as shown in Fig. 1(b), both wavefronts are tangent to each of the wavelets centered along the aspheric surface passing through the vertex V of the lens. Note that the continuous circles in the figures indicate the envelope of the forward propagating wave and the dotted circles form the envelope of a backward propagating wave that is a nonphysical solution and is not taken into account in the subsequent analysis. We consider exclusively from Eq. (5) a wavefront progressing in the forward direction and reducing further it becomes where the subscript pc means plano-convex. It is important to say that Eq. (6) gives the coordinates of the locus of points that parametrically represent the zero-distance phase front or the archetype wavefront, which is produced by a plano-convex aspheric lens in a meridional plane as a function of h for a plane wavefront incident on the lens. From Eq. (6) assuming that the radical , and solving for h whose solution is defined by hc, the rays undergo Total Internal Reflection (TIR), see for example [5]. Physically, if hc > H, then all the wavelets contribute to the formation of the refracted wavefront, being tangent to the wavefront as is shown in Fig. 2(a). On the other hand, if hc < H there are wavelets which do not contribute to the formation of the wavefront, and therefore they are not tangent to the wavefront as is shown in Fig. 2(b).3. An alternative method
A congruence is formed by rays propagating in an homogeneous isotropic medium, which originate from a point source being the congruence perpendicular to the phase front. The theorem of Malus and Dupin states that the propagating phase front in a uniform medium, is everywhere normal to the congruence of rays [10]. Since the caustic is the locus of the principal centers of curvature of a wavefront, also the caustic is the evolute of the zero-distance phase front, in others words the phase front is the involute of the caustic. Mathematically, if a curve 𝒜 is the involute of curve ℬ, then ℬ is the evolute of 𝒜, and vice versa. Let 𝒞 be a curve in parametric form 𝒞 = [f(h), g(h)], in such a way that the equation of the involute is given by
where for simplicity the constant of integration has been considered null. According to Eq. (7) in reference [5], the caustic surface (zpc, ypc) for plano-convex aspheric lenses was written as where S″h is the second derivative from Eq. (2). It is important to say that Eq. (8) gives the coordinates of the locus of points that parametrically represents the envelope of the family of refracted rays produced by an aspheric lens in a meridional plane as a function of h when the point source is placed at infinity. Additionally, Eq. (8) is similar to Eq. (29a) for C+ which provides the radii of curvature of the wavefront along the tangential plane according to reference [3], by using the follows considerations: r = h, z(r) = Sh, n1 = ni, n2 = na and t = 0. By substituting Eq. (8) into Eq.(7) for zpc → f and ypc → g and reducing further we obtain Eq. (6). It is worth stating that the constant of integration omitted in Eq. (7) simply translates the vertex of the zero-distance phase front along the optical axis.4. Exact wavefront for plano-convex aspheric lenses
Since either the zero-distance phase front or the archetype wavefront can also be considered as a source of such a congruence, all ensuing phase fronts are parallel in the geometrical sense. That is, points on each phase surface are equidistant where the distances are measured along the common normal between points on each of the phase front surfaces, according to Huygens’ principle as are shown in Fig. 3(a). Let 𝒲0 be a point on a phase front whose profile is given parametrically by 𝒲0 = [F(h), G(h)], then 𝒲 will be on a parallel phase front at a distance ℒ according to the following equation
In this way, Eq. (9) for 𝒲(−,+) provides a retrogade wavefront and for 𝒲(+,−) provides the forward wavefront. Substituting Eq. (6) into Eq. (9) for Zpc → F and Ypc → G and reducing further we get where the subscript ‖pc means wavefronts propagating in parallel form produced by plano-convex lenses, where (Zpc, Ypc) are defined in Eq. (6). We can see that for ℒ = 0, Eq. (10) is reduced to Eq. (6) providing the zero-distance phase front or archetype wavefront as expected, and is shown in Fig. 3(b).5. Wavefront for aspheric convex-plano lenses
Let H be the height for the marginal wavelet, so that the family of wavelets centered along the first aspheric surface of the lens can be represented by the following formula
where SH is evaluated at H from Eq. (2), considering in this configuration that c > 0. For example from Eq. (11), for h = Y = 0, then Sh = 0 and therefore it reduces to: niZ = naSH as is shown in Fig. 4(a). Following all the steps explained above to obtain the envelopes to the wavelets of Eq. (11), these can be represented by as is shown in Fig. 4(b). Finally, from Eq. (12) we choose the wavefront progressing in the forward direction (Z−, Y−) inside the lens and reducing further it becomes It is simple to see that Eq. (13) represents the wavefront refracted inside the aspheric lens, in such a way that we can consider a new family of wavelets centered along the wavefront (Zin, Yin) and propagate it inside the lens up to the marginal wavelet reaching the plane face of the lens, propagating along the marginal ray instead of marginal wavelet touching the border of the lens at (t, H) as is shown in Fig. 5(a). In other words, we obtain the length ℒH between a point on aspheric surface (SH, H) and the plane face of the lens (t, yi[H, t]), where yi[H, t] represents the height of the marginal ray impinging on the plane face of the lens. By using the Snell’s law we provide the equation for the rays propagating inside the lens which can be written as where h and z are variables and we have considered that all the parameters of the lens are constants. Substituting h → H and z → t into Eq. (14), we calculate the distance along the refracted marginal ray according to ℒH = {(yi [H, t] − H)2 + (t − SH)2}1/2, and reducing further we get S′H is given from Eq. (4) evaluated at H, we can see that ℒH is constant and depends on all parameters involved in the process of refraction. By substituting Eq. (13) into Eq. (9) and reducing further we obtain a wavefront translated a distance ℒH as is shown in Fig. 5(b), given by Finally, we can see that Eq. (16) provides the wavefront which will be refracted outside of the aspheric lens. Substituting h → 0 into Eq. (16) provides the distance covered between the vertex 𝒪 of the lens and the translated wavefront along the optical axis, it becomes To obtain the remaining distance to refract the translated wavefront outside of the lens we subtract the length Z‖in [0] from the thickness of the lens t, and we get or alternatively, we can provide the distance between the plane face of the lens and the wave-front Z‖in, Y‖in for arbitrary heights h obtaining In order to obtain the refracted wavefront outside of the aspheric lens, we consider that the Optical Path Length, OPL = naℒa + niℒi, is null, thus we have that ℒa = −(ni/na)ℒi, therefore the wavelets for the refracted wavefront outside of the lens are represented by where yi[h, t] is given from Eq. (14) evaluated at z = t and ℒi is given from Eq. (17) as is shown in Fig. 6(a). Following the steps explained above, we can obtain the zero-distance phase for the refracted wavefront outside of the aspheric lens. An alternative method is given as follows. We need to provide a length ℒa outside of the lens between the plane face of the lens at (t, yi[h, t]) and an arbitrary length z propagating along the refracted rays at (z, yo[h, z]), where yo[h, z] represents the equation for the rays propagating outside of the lens with an arbitrary height h, in such a way that by using Snell’s law, this can be written as so that the length ℒa is calculated as ℒa = {(z −t)2 + (yo[h, z] − yi[h, t])2}1/2, by using Eq. (19) and reducing further we get Assuming the condition ℒa = −(ni/na)ℒi and substituting Eqs. (17) and (20) it is reduced to where we have defined . By solving from Eq. (21) for z and substituting this value into Eq. (19) we obtain the zero-distance phase where the subscript cp means convex-plano. It is important to say that Eq. (22) provides the coordinates of the locus of points that parametrically represents the zero-distance phase front, which is produced by an aspheric convex-plano lens in a meridional plane as a function of h for a plane wavefront incident on the lens as is shown in Fig. 6(b). Substitution of Eq. (22) into Eq. (9) yields the exact wavefront which is propagated an arbitrary distance ℒ, according to where for simplicity and brevity we have omitted writing this equation in its complete form. In Fig. 7(a) we have represented several wavefronts at different distances, and a zoom near the back focal length area is shown in Fig. 7(b).For the ray tracing we have considered a lens with F/# = 0.8, by using the following parameters in both plano-convex and convex-plano lens configurations: na = 1, ni = 1.5111, R = 30.67mm, diameter D = 75mm, t = 35.5mm, k = −0.905 and entrance aperture H = ±D/2, whose as-pheric coefficients are given in Table 1.
It is worth noting that since Eq. (23) defines the exact relation between aspheric terms and the resultant wavefront (and caustic), it can contribute to both the optimization of the focal properties of aspheric lenses and to the solution of the converse problem: the design of refractive elements for the synthesis of arbitrary wavefronts. In other words, knowing the wavefront, we know the aspheric terms. This provides an alternative route to wavefront engineering where a plane wave needs to be transformed into a complex waveform and vice versa. For example, in microscopy, one can imagine using arbitrary aspheric terms in point spread function engineering [11] or aberration correction [12]. In applications of current interest, this formalism may be of more general interest where the form of a lens can be controlled in autofocus applications and even in refractive elements for terahertz applications [13].
Conclusions
We have obtained simple formulas for the wavefronts produced by positive convex-plano and plano-convex aspheric lenses having an arbitrary number of aspheric terms, considering a plane wavefront propagating along the optical axis. The shape of the wavefronts can be modified by changing the parameters of the lens and the distance where they are observed. When the condition of total internal reflexion is satisfied the wavelets do not contribute to the formation of the refracted wavefront, since these wavelets are unlimited therefore they are not perpendicular to the refracted rays. Upon inspection one can see that our wavefront coincides with the tangential solution according to [3]. We do not treat the sagittal plane that they does, but we could if needed extend our results by using differential geometry to obtain the wavefront for the sagittal plane. The net result is a simplification in the formalism that gives us the simple set of equations 10 and 23 in this manuscript that allow a direct calculation. We believe that this method for obtaining the wavefronts reported here are straightforward, giving a relationship between caustics and wavefronts.
Acknowledgments
This work has been partially supported by a Programa de Apoyo a Proyectos de Investigación e Innovación Tecnológica – Universidad Nacional Autónoma de México (PAPIIT-UNAM) under project # IN112612, and Consejo Nacional de Ciencia y Tecnología (CONACyT) under project # 168570. The corresponding author is grateful to I. Goméz-García for their valuable assistance and comments.
References and links
1. A. Epple and H. Wang, “Design to manufacture from the perspective of optical design and fabrication,” Optical Fabrication and Testing, OSA Technical Digest, OFB1. (2008).
2. O. N. Stavroudis, The Mathematics of Geometrical and Physical Optics, The K-function and its Ramifications (Wiley-VCH Verlag GmbH & Co.KGaA, 2006), Chap. 12, 179–186. [CrossRef]
3. D. L. Shealy and J. A. Hoffnagle, “Wavefront and caustics of a plane wave refracted by an arbitrary surface,” J. Opt. Soc. Am. A , 25, 2370–2382 (2008). [CrossRef]
4. J. A. Hoffnagle and D. L. Shealy, “Refracting the k-function: Stavroudis’s solution to the eikonal equation for multielement optical systems,” J. Opt. Soc. Am. A , 28, 1312–1321 (2011). [CrossRef]
5. M. Avendaño-Alejo, “Caustics in a meridional plane produced by plano-convex aspheric lenses,” J. Opt. Soc. Am. A , 30, 501–508 (2013). [CrossRef]
6. G. W. Forbes, “Robust, efficient computational methods for axially symmetric optical aspheres,” Opt. Exp. , 18, 19700–19712 (2010). [CrossRef]
7. G. W. Forbes, “Shape specification for axially symmetric optical surfaces,” Opt. Exp. , 15, 5218–5226 (2007). [CrossRef]
8. M. Avendaño-Alejo and O. N. Stavroudis, “Huygens’s Principle and Rays in Uniaxial Anisotropic Media II. Crystal Axis with Arbitrary Orientation,” J. Opt. Soc. Am. A. 19, 1674–1679 (2002). [CrossRef]
9. F. A. Jenkins and H. E. White, Fundamentals of Optics (Mc. Graw-Hill, 1976), Chap. 3, 44–57.
10. S. Cornbleet, Microwave and Optical Ray Geometry (Wiley, London, 1984), 11–35.
11. M. A. A. Neil, T. Wilson, and R. Juskaitis, “A wavefront generator for complex pupil function synthesis and point spread function engineering,” Journal of Microscopy , 197, 219–223 (2000). [CrossRef] [PubMed]
12. M. J. Booth, M. A. A. Neil, and T. Wilson, “Aberration correction for confocal imaging in refractive index mismatched media,” Journal of Microscopy , 192, 90–98 (1998). [CrossRef]
13. B. Scherger, C. Jordens, and M. Koch, “Variable-focus terahertz lens,” Opt. Exp. , 19, 4528–4535 (2011). [CrossRef]