Abstract
We present a new formalism, based on Richards–Wolf theory, to rigorously model nonparaxial focusing of radially polarized electromagnetic beams with freeform wavefront. The beams can be expressed in terms of Zernike polynomials. Our approach is validated by comparing known results obtained by Richards–Wolf theory. Our integral representation is compliant with diffraction theory, is thoroughly discussed and solved for various freeform wavefront that, so far, have not been treated analytically. The extension of the method to other polarization states is straightforward.
© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
The modeling of nonparaxial beams is necessary in several area of science, including high-resolution microscopy [1–4], optical trapping [5], electron acceleration [6], and optical vortex knots [7,8]. Generally, papers have focused on beam shapes with prescribed polarization states that are critical for mentioned applications [9–16].
In focusing systems, the scalar approach can be considered if the angles of incidence are in the paraxial regime. However, if the numerical aperture increases the scalar approach is no longer valid and a vectorial treatment to the electromagnetic fields in the focal volume is needed. A well-know approach of the mentioned to this is the Richards–Wolf theory (RWT) [17–20]. In its initial formulation RWT assumes an input plane phase collimated beam and an output a spherical beam (stigmatic). The Richards–Wolf theory has been extended to radially and azimuthally polarized electromagnetic beams by axisymmetric aberrated systems [21,22].
The goal of this paper is to further extend RWT for non-paraxial focusing for a freeform wavefront in the optical system exit pupil. Notice that the original Richards-Wolf vector diffraction only considers perfect stigmatic optical systems. This generalization freeform wavefronts are also considered. To achieve this goal, we work directly with the propagation vector and the unit vector representing the polarization direction. We put emphasis that both must always be orthogonal. Also, we take account of the Jacobian metric during a change of variables that we will see in section 2. The manuscript is divided as follows, first an introduction to RWT is presented. Then, the new formalism is developed to compute integrals for an arbitrary wavefronts. Finally, examples are presented and conclusions are given.
2. Richard-Wolf theory
Richards–Wolf theory, consists in vector diffraction integral useful to analyze tightly focused beams. Richards–Wolf theory considers a plane wave incident on a focusing system that is transformed into a converging spherical wave. Under this paradigm, this spherical wave can be expanded into an angular spectrum of plane waves. Those plane waves contribute to the field in the focus the propagation directions of which correspond to the geometric optical rays. Therefore the field in the focus can be evaluated by superposing those plane waves. Based the aforementioned considerations and using cylindrical coordinates $(r,\phi,z)$ the electric field at the focus is,
As mention earlier RWT considers that the output waterfront is spherical and it is the superposition of several plane waves. Therefore the components of the wave vector ${\mathbf k}$ turn to,
where $k\equiv |{\mathbf k}|=2\pi /\lambda$ is the wave number, $\lambda$ is the wavelength, $\alpha$ and $\beta$ are the polar and the azimuthal angles of the incident plane wave component, respectively, the domain of the angles are $\alpha \in [\alpha _{\text {min}},\alpha _{\text {max}}]$ and $\beta \in [0,2\pi ]$. See Fig. 1. Tacking into account Eq. (2) and a radially polarized beam, the polarization vector becomes,Replacing Eqs. (2) (3) and (4) in Eq. (1), ${\mathbf E}(r,\phi,z)$ becomes,
Eq. (5) expresses the electric field near the focus when the system is stigmatic, which in other means, it has a single-point image. This means that the wavefront of the focusing beam is a perfect sphere. In the next section, we are going to consider the case when the wavefront is not necessary as a sphere.
3. Generalisation of the Richards-Wolf theory
The fundamental background of the RWT is that the focused spherical wave can be expanded into an angular spectrum of plane waves. Without loosing generalization we can also express the any other focused wavefront as sum of plane waves with the appropriate wave vector ${\mathbf k}$ and its Cartesian components $k_x$, $k_y$, and $k_z$. Therefore, given a wavefront $w(x,y)$ and according to the Malus-Dupin theorem, the wave vector ${\mathbf k}$ of the wavefront $w(x,y)$ has the direction of normal to $w(x,y)$ [23,24]. Here, $w(x, y)$ is an aberrated wavefront generated by an imperfect system. Then, the cosine directors of the wavefront wave vector ${\mathbf k}$ are,
Observe that Eq. (6) and Eq. (7) are perpendicular to each other, also notice that they are in terms of $(x,y)$, so we need to express $x$ and $y$ in terms of $\alpha$ and $\beta$. Let be $R$ the radius of the exit pupil then we can express, $(x,y)$ in terms or $(\beta,\alpha )$,
Since the output wavefront is not spherical Eq. (5) is not valid, thus to get the differentials $\text {d}\alpha \text {d}\beta$ we need to compute the Jacobian $J(\alpha,\beta )$,
where,To get Jacobian $J(\alpha,\beta )$ first we need to replace $x=R\sin \alpha \cos \beta$ and $y=R\sin \alpha \sin \beta$ in the cosine directors and then directly compute Eq. (10). Therefore, the Richard-Wolf generalized integral for arbitrary focusing wavefronts is given by,
Eq. (11) is the most important equation in this manuscript. It accurate describes the diffraction pattern generated by an arbitrary wavefront. Notice that all the parameters inside Eq. (11) should be expressed in terms of $(\alpha,\beta )$. Therefore, $\hat {{\mathbf a}}(\alpha,\beta )$ is the composite function between Eq. (7) and (8). The cosine directors $k_x$, $k_y$ and $k_z$ in Eq. (11) are given by the composite function between Eq. (6) and (8), respectively. Finally, $J(\alpha,\beta )$ is given by Eq. (10).
Observe that to get, azimuthally polarized components is straight forward. Minor modifications need to be added in Eq. (7). The components of Eq. (7) must be multiplied by minus one.
4. Examples and discussion
In this section we explore several diffraction patterns obtained directly from Eq. (11). As a first example, we are going to consider a spherical wavefront $w(x,y)=R-\sqrt {R^{2}-x^{2}-y^{2}}$, the purpose of this example is to show that the equations presented in section 3 naturally converge to the equations of section 2, which is the standard RWT. We start with the cosine directors, that for the spherical wavefront are,
Evaluating Eq. (12) with $x=R\sin \alpha \cos \beta$ and $y=R\sin \alpha \sin \beta$, we get Eq. (2). For the spherical wavefront $\mathbf {\hat {a}}$ is given by
The composite function between Eqs. (13) and (8) and becomes Eq. (3). Notice that during this procedure $R$ is cancelled in Eq. (3). The Jacobian $J(\alpha,\beta )$ can directly computed using Eq. (10) and Eq. (8) and divided by $k_z(\alpha,\beta )$, that turns to be Eq. (4). If we replacing replace Eqs. (12), (13), (10) with (8) in Eq. (11) we get the Richard-Wolf Integral, Eq. (1).
In the next example, the wavefront is give by the Zernike polynomial $Z_{2}^{0}=w(x,y)=c(-1+2x^{2}+2y^{2})$ where $c\in \mathbb {R}$. The cosine directors are computed with composite function of Eq. (6) and (8),
The polarization unit vector $\hat {{\mathbf a}}(\alpha,\beta )$ is the composite function between Eq. (7) and (8), where $w(x,y)=c(-1+2x^{2}+2y^{2})$,
Finally, we compute the Jacobian divided by $k_z(\alpha, \beta )$,
To compute the diffraction pattern generated by $Z_{2}^{0}=w(x,y)=c(-1+2x^{2}+2y^{2})$ we replace the Eqs. (14), (15), (16) in Eq. (11) and compute the integral. In Fig. 2 we present the intensity profile for a spherical wavefront in orange. In (a) is the axial profile and (b) is the lateral profile. In the figure but in black are the axial and lateral intensity profiles of $w(x,y)=c(-1+2x^{2}+2y^{2})$, respectively. The specifications of each example can be seen in the caption of Fig. 2.
The orange curve corresponds to an astigmatic diffraction pattern, while black has aberrations. The full width at half maximum (FWHM) of the orange profiles are smaller than the ones in black. FWHM for the orange curves are $0.52\lambda$ and $0.23\lambda$ in the $z$ and $r$ profiles, respectively. For the black curves the FWHM is $0.57\lambda$ in the $z$ profile and $0.24\lambda$ in for $r$ profile.
From the example of Fig. 2 (a) and (b) it is clear that the method presented in this manuscript converges to RWT when the wavefront is spherical (astigmatic).
Notice that in this computation we set the apodization factor for both wavefronts as one, $q(\alpha,\beta )\to 1$, since we are only considering the wavefront at the exit pupil and the optical system that produces the aforementioned wavefront. Thus, we do not consider the distribution of ray across the optical system. Also notice that the transverse electric field at the output is radially polarized; This means that the optical system does not produces any polarization aberration. But considering polarization aberrations is straight forward by directly modifying $\hat {{\mathbf a}}(\alpha,\beta )$. The illumination pattern $l_0(\alpha,\beta )$ has been taken as uniform as well.
The method presented here also can be used to model freeform wavefronts. As an example of these freeform wavefronts we choose $w(x,y)=a \left (b y^{2}+x^{2}\right )$, where $a,b\in \mathbb {R}$. The cosine directors of the normal of the wavefront can be expressed in terms of $\alpha,\beta$ are,
The Jacobian divided by $k_z(\alpha, \beta )$ for the $w(x,y)=a \left (b y^{2}+x^{2}\right )$, it is given by,
Inserting Eqs. (17), (18), (19) in (11) and directly computing the integral we get the diffraction patterns presented in Fig. 3. From the figure it is clear that the diffraction patterns is nor radially symmetric as is expected since $w(x,y)=a(b y^{2}+x^{2})$ is not radially symmetry. In Fig. 4, (a), (b) and (c) are the diffraction patterns in the $z$, $x$ and $y$ directions, respectively. Notice that the distribution of energy is wider than in Fig. 2, this is due the intrinsic aberration suffered by the example wavefront.
Finally, we computed the diffraction pattern generated by a spherical mirror. The diffraction pattern is displayed in Fig. 4. In the left side is the right side of the mentioned figure correspond to the diffraction pattern in the $z$ and $r$ axis, respectively. We set the apodization factor as one, $q(\alpha )\to 1$, we use a Gaussian illumination, $l_0(\alpha )\to \alpha\,\textrm{exp}(-\alpha ^{2}/\Delta \alpha ^{2})$ and we $w(x,y)$ the following wavefront,
It is important to remark the limitations of the method. The method begins with from the normal vector of the wavefront $w(x,y)$ which is in indeed, the propagation vector ${\mathbf k}$, if this vector is such that its respective Jacobian $J(\alpha,\beta )$ is not zero the computation can be done. An example of a wavefront with a zero Jacobian is a plane wave. This result, should not be surprising since the plane it is not a focusing beam. Therefore, the diffraction pattern of diverging beams can not be computed with the method presented in this paper.
Also notice that the method presented here considers the Jacobian $J(\alpha,\beta )$ as well the orthogonality of ${\mathbf k}$ and $\hat {{\mathbf a}}(\alpha,\beta )$. In the generalization presented in [11,21,22], both motioned aspects are not considered and only aberration is added in the axial axis. [11,21,22] ignore these aspects, thus as a consequence ${\mathbf k}$ and $\hat {{\mathbf a}}(\alpha,\beta )$ are no longer orthogonal and for the case $|{\mathbf k}|$ is no longer equal to $k$. Therefore, in the best knowledge of the authors, the generalization presented here is more robust than the previews ones [11,21,22]. We consider the method presented in this manuscript as a great bridge between the two paradigms, geometrical optics and diffraction theory.
5. Conclusions
In this manuscript, we presented a rigorous generalization of the RWT to arbitrary wavefronts. The standard RWT only considers spherical waterfronts (no aberration) but this new theory considers any arbitrary wavefronts. To do so, This generalization considers that ${\mathbf k}$ and $\hat {{\mathbf a}}(\alpha,\beta )$ should be always perpendicular as well considers that the Jacobian in the changes of variables for the first time. As a result, the method presented here successfully computes the diffraction patterns of freeform wavefronts, including spherical, Zernike polynomials and any freeform waterfronts. We also have discussed its limitations.
The natural step to follow in this formalism is the study of the diffraction pattern of freeforms wavefront. For example one potential application is to study the effect of optical tolerance like component decenters (coma), or off axis aberrations for non-paraxial focusing. Also the study diffraction patterns generated by on-axis and off-axis fields in microscope objectives or freeform telescopes, where phase and intensity spatial distributions of the beam should be considered.
Appendix
Here, we show how to get a wavefront generated by a spherical mirror when it reflects a plane wave aligned along with the mirror. From refraction’s law at the surface, the cosine directors of the output wavefront are,
Equaling Eq. (21) with Eq. (6) and solving for $\partial (w_x)$ and $\partial (w_y)$, we get,
From Eq. (22), If we integrate $\partial _xw$ respect to $x$, or we integrate $\partial _yw$ respect to $y$, we get the same result, the wavefront $w(x,y)$ expressed in Eq. (20).
Funding
Universite Laval (NSERQ RGPIN-2016-05962)).
Disclosures
The author declares no conflicts of interest.
Data availability
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the author upon reasonable request.
References
1. G. Thériault, Y. De Koninck, and N. McCarthy, “Extended depth of field microscopy for rapid volumetric two-photon imaging,” Opt. Express 21(8), 10095–10104 (2013). [CrossRef]
2. F. O. Fahrbach, V. Gurchenkov, K. Alessandri, P. Nassoy, and A. Rohrbach, “Light-sheet microscopy in thick media using scanned bessel beams and two-photon fluorescence excitation,” Opt. Express 21(11), 13824–13839 (2013). [CrossRef]
3. J. Liu, J. Tan, T. Wilson, and C. Zhong, “Rigorous theory on elliptical mirror focusing for point scanning microscopy,” Opt. Express 20(6), 6175–6184 (2012). [CrossRef]
4. J. Liu, ed., Elliptical Mirrors, 2053–2563 (IOP Publishing, 2018).
5. Q. Zhan, “Trapping metallic rayleigh particles with radial polarization,” Opt. Express 12(15), 3377–3382 (2004). [CrossRef]
6. C. Varin, M. Piché, and M. A. Porras, “Acceleration of electrons from rest to gev energies by ultrashort transverse magnetic laser pulses in free space,” Phys. Rev. E 71(2), 026603 (2005). [CrossRef]
7. D. Sugic, M. R. Dennis, F. Nori, and K. Y. Bliokh, “Knotted polarizations and spin in three-dimensional polychromatic waves,” Phys. Rev. Res. 2(4), 042045 (2020). [CrossRef]
8. M. F. Ferrer-Garcia, A. D’Errico, H. Larocque, A. Sit, and E. Karimi, “Polychromatic electric field knots,” Phys. Rev. Res. 3(3), 033226 (2021). [CrossRef]
9. K. Jahn and N. Bokor, “Solving the inverse problem of high numerical aperture focusing using vector slepian harmonics and vector slepian multipole fields,” Opt. Commun. 288, 13–16 (2013). [CrossRef]
10. R. Kant, “An analytical solution of vector diffraction for focusing optical systems,” J. Mod. Opt. 40(2), 337–347 (1993). [CrossRef]
11. A. April, P. Bilodeau, and M. Piché, “Focusing a tm 01 beam with a slightly tilted parabolic mirror,” Opt. Express 19(10), 9201–9212 (2011). [CrossRef]
12. S. Quabis, R. Dorn, M. Eberler, O. Glöckl, and G. Leuchs, “The focus of light–theoretical calculation and experimental tomographic reconstruction,” Appl. Phys. B 72(1), 109–113 (2001). [CrossRef]
13. H. Dehez, A. April, and M. Piché, “Needles of longitudinally polarized light: guidelines for minimum spot size and tunable axial extent,” Opt. Express 20(14), 14891–14905 (2012). [CrossRef]
14. A. April and M. Piché, “4π Focusing of TM01 beams under nonparaxial conditions,” Opt. Express 18(21), 22128 (2010). [CrossRef]
15. R. G. González-Acu na, “Stigmatic singlet with a user-defined apodization pupil function,” J. Opt. Soc. Am. A 39(2), 213–218 (2022). [CrossRef]
16. R. G. González-Acu na, J. Borne, and S. Thibault, “On the diffraction of a high-na aplanatic and stigmatic singlet,” J. Opt. Soc. Am. A 38(9), 1332–1338 (2021). [CrossRef]
17. H. Green and E. Wolf, “A scalar representation of electromagnetic fields,” Proc. Phys. Soc., London, Sect. A 66(12), 1129–1137 (1953). [CrossRef]
18. E. Wolf, “A scalar representation of electromagnetic fields: Ii,” Proc. Phys. Soc., London 74(3), 269–280 (1959). [CrossRef]
19. E. Wolf, “Electromagnetic diffraction in optical systems-i. an integral representation of the image field,” Proc. R. Soc. Lond. A 253(1274), 349–357 (1959). [CrossRef]
20. B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems, ii. structure of the image field in an aplanatic system,” Proc. R. Soc. Lond. A 253(1274), 358–379 (1959). [CrossRef]
21. D. Panneton, G. St-Onge, M. Piché, and S. Thibault, “Exact vectorial model for nonparaxial focusing by arbitrary axisymmetric surfaces,” J. Opt. Soc. Am. A 33(5), 801–810 (2016). [CrossRef]
22. J. Borne, D. Panneton, M. Piché, and S. Thibault, “Analytical inversion of the focusing of high-numerical-aperture aplanatic systems,” J. Opt. Soc. Am. A 36(10), 1642–1647 (2019). [CrossRef]
23. M. Born and E. Wolf, Principles of optics: electromagnetic theory of propagation, interference and diffraction of light (Elsevier, 2013).
24. R. G. González-Acu na and H. A. Chaparro-Romo, Stigmatic Optics (Institute of Physics Publishing, 2020).