Accepted: February 17, 2023 | Published Online: February 19, 2023

Citation: Xu Q, Yan S, Wu Z (2023) Focusing of Anoblique Incident Laguerre-Gaussian Vector Beam by a High Numerical-Aperture Lens. Int J Opt Photonic Eng 8:053.

Copyright: © 2023 Xu Q, et al. This This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

We describe diffractive and focusing properties of oblique incident vector vortex beams under tight focusing conditions. By dividing the effective pupil area of the focus lens into two parts for simple integration, using the electromagnetic wave polarization vector diffraction theory, the intensity distributions of oblique incident vector Laguerre-Gaussian beams focused by high numerical aperture lens are calculated and analyzed, the structural changes in the transverse and the axial direction near the focus plane are observed, numerical simulation show that, field intensity distributions are changed with the different oblique incident angle and topological charge of Laguerre-Gaussian beam, the distortion and asymmetry of intensity distribution are shown in different cases. This research can be applied for trapping and manipulating particles, super-resolution fluorescence microscopes.

Keywords

Laguerre-Gaussian vector beam, Vector diffraction theory, Oblique incident, High numerical-aperture Lens, Intensity distribution

Introduction

Recently, vortex beams have been studied extensively and have wide applications, such as optical communications [1-3], optical imaging [4]. Generating multiple vortex patterns beam by a diffraction grating [5], the conversion of polarized laser beams from a radial to an azimuthal polarization, or vice versa [6], the generation of beams with different polarization states in an anisotropic crystal [7]. Application of the vortex beam for trapping and manipulating particles, such as controlling the chirality of metal nanostructures [8], producing chiral nanoneedles on thin films [9]. The use of optical vortices to generate a controllable optical cage [10]. it is quite difficult to create 3D optical cage with identical in all directions light barriers at the usual focusing (from one side), since the axial resolution is several times worse than the transverse one, to decrease the axial extension of the focal spot, different techniques are used, including 4pi scheme with two opposing lenses [11-13]. Also, three-dimensional focus shaping may be realized by supplementing the focusing lens with annular, or choosing the suitable polarization states of vector beams [14], or more complex diffractive optical elements [15,16].

A tight-focusing vortex beam for achieving the smallest possible focal spot is an essential issue in a variety of applications. To reduce the size of the focal spot under sharp focusing, including high-aperture focusing system [17,18], stimulated emission depletion microscopy [19]. Laguerre-Gaussian beams being focused through a stratified medium [20].

The astigmatic transformations of focused vortex beam were studied by many scientists, such as an oblique incidence of an axially symmetric beam onto an optical element produce elliptic beams [21,22]. This effect is closely related to astigmatic transformations, used for the formation of Laguerre-Gauss modes from Hermite-Gauss modes [23,24] and for analyzing the structure of vortex beams [25-27]. Moreover, the sharp focusing of Laguerre-Gauss beams with astigmatic distortion, polarization state affects the longitudinal component of the electric vector of the light field [28].

But by high numerical aperture (NA) lenses, beam field of oblique incident vortex beam have some new character, can be used for trapping and manipulating particles, it need more study, Motivated by these facts, a basic goal of this paper is to examine the diffractive and focusing properties of oblique incident Laguerre-Gaussian (LG) beams focused by high numerical aperture (NA) lenses near the focal region.

Analysis and Methodology

An LG beam oblique is incident on high NA lens, which is schematically shown in Figure 1a. The angle between incident beam and z₀ axis is α, we rotate the coordinate system to make incident beam parallel to z axis, shown in Figure 1b.

We note the arrangement of the optical lens system in Figure 1, and use the x-axis polarized Laguerre-Gaussian beam [29].

$$\begin{array}{l}{\overrightarrow{E}}_0\left(x,y,z\right)\text{ = }\exp \left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m\text{ (1)}\\ \times L_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]\exp \left(im\varphi \right)\overrightarrow{i}\end{array}$$

$L_{\left(n-\left|m\right|\right)/2}^{\left|m\right|}$ are the generalized Laguerre polynomials, and the integers n, m obey the relation,

$$n\text{ = }\left|m\right|,\left|m\right|+2,\left|m\right|+4,\mathrm{...}\text{ (2)}$$

Fresnel approximation of the classical scalar diffraction theory can only be applied to the paraxial optical system or the low numerical aperture. It has lost its effectiveness in the high NA system. In order to accurately predict the performance of high NA, we use the electromagnetic wave polarization vector diffraction theory. The space dependent electric vectors $\overrightarrow{E}\left(\overrightarrow{r}\right)$ in the image space of the optical system have the integral representation [30,31]

$$\overrightarrow{E}(\overrightarrow{r})\text{ = }\frac{e^{ikr\text{'}}}{i\lambda r\text{'}}{\displaystyle \underset{\Sigma }{\iint }{\overrightarrow{E}}_S}{e}^{ik\widehat{N}\cdot \overrightarrow{r}}ds\text{ (3)}$$

Where $\lambda$ is the wavelength, Σ is the wavefront surface over the exit pupil of a lens, $\widehat{N}$ is a unit ray vector (normal to the wavefront), ${\overrightarrow{E}}_s$ is electric field at the exit pupil, k is the wave number in medium. O is geometrical focus, X is an observation point; Y is a point on the exit pupil surface, the distance OX is $\left|\overrightarrow{r}\right|$, the distance OY is $\left|\overrightarrow{r\text{'}}\right|$, the distance YX is $R=\left|\overrightarrow{r}-\overrightarrow{r\text{'}}\right|$ ( shown in Figure 1).

The evaluation of the vector Debye integral requires a mathematical explanation of the effective pupil area in Figure 1. As illustrated in Figure 2a, the original circular pupil region and the oblique region overlap the portion to form an effective pupil. The normalized pupil function $P\left(\theta ,\varphi \right)$ is derived from the exit pupil of the image space of the spherical coordinates ($\theta$ : polar angle, $\varphi$: azimuthal angle). In order to facilitate the calculation of the integral, the circular asymmetric pupil in the oblique plane imaging is divided into ${\Sigma }_1$ (rotationally symmetric part) and ${\Sigma }_2$ (the rest area). The pupil functions for ${\Sigma }_1$ and ${\Sigma }_2$ are [30]

$$P_{{\Sigma }_1}\left(\theta ,\varphi \right)\text{ = }\left\{\begin{array}{cc}1,& \theta \in \left[0,{\theta }_c\right],\varphi \in \left[0,2\pi \right],y_c<0\\ 0,& otherwise\end{array}\right.\text{ (4)}$$

$$P_{{\Sigma }_2}\left(\theta ,\varphi \right)\text{ = }\left\{\begin{array}{cc}1,& \theta \in \left[{\theta }_c,{\theta }_{\max }\right],\varphi \in \left[{\varphi }_1(\theta ),{\varphi }_2(\theta )\right]\\ 0,& otherwise\end{array}\right.\text{ (5)}$$

Where ${\theta }_C$, ${\theta }_{\max }$ , ${\varphi }_1\left(\theta \right)$ , and ${\varphi }_2\left(\theta \right)$ are defined in follow. The values of $\theta$, $\varphi$ are expressed as a function of $\alpha$, NA and n.

In Figure 2, the equation of the plane ABC is $y-z\mathrm{c}tg\alpha +\sqrt{1-{\left(NA/n\right)}^2}\sin \alpha \text{ = }0$, plugging this equation into the unit sphere equation give the coordinate values of the point C as

$${\theta }_c\text{ = }{\cos }^{-1}\left|z_c\right|\text{ = }{\cos }^{-1}\left|\frac{NA}{n}\sin \alpha +\sqrt{1-{\left(\frac{NA}{n}\right)}^2}\cos \alpha \right|\text{ (6)}$$

$${\theta }_{\max }\text{ = }{\sin }^{-1}\left(\frac{NA}{n}\right)\text{ (7)}$$

$$\begin{array}{l}\left(x_c,y_c,z_c\right)\text{ = }\left(0,\frac{NA}{n}\cos \alpha +\sqrt{1-{\left(\frac{NA}{n}\right)}^2}\sin \alpha ,\right.\text{ (8)}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}\left.\frac{NA}{n}\sin \alpha -\sqrt{1-{\left(\frac{NA}{n}\right)}^2}\cos \alpha \right)\end{array}$$

$$\begin{array}{l}{\varphi }_1\left(\theta \right)\text{ = }arc\sin \left(\frac{y}{\sin \theta }\right)\left|{}_{y=z\cot \alpha +\frac{\sqrt{1-{\left(NA/n\right)}^2}}{\sin \alpha }}\right.\text{ (9)}\\ \begin{array}{cc}& \end{array}\begin{array}{c}\\ \end{array}=arc\sin \left|-\cot \alpha \cot \theta +\sqrt{1-{\left(\frac{NA}{n}\right)}^2}\csc \alpha csc\theta \right|\end{array}$$

$${\varphi }_1\left(\theta \right)+{\varphi }_2\left(\theta \right)\text{ = }\pi \text{ (10)}$$

The electric field near the focal area of the oblique lens can be obtained from equation (1) to be

$$\begin{array}{l}\left[\begin{array}{c}{E}_x\\ E_y\\ E_z\end{array}\right]\text{ = }{\displaystyle \underset{\Sigma }{\iint }\exp \left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}\cdot \exp \left(im\varphi \right)\left[\begin{array}{c}{\sin }^{\text{2}}\varphi +\cos \theta {\cos }^{\text{2}}\varphi \\ -\sin \varphi \cos \varphi +\cos \theta \sin \varphi \cos \varphi \\ \cos \varphi \sin \theta \end{array}\right]\text{ (11)}\\ \\ \times \exp \left[-ik\left(x\sin \theta \cos \varphi +y\sin \theta \sin \varphi -z\cos \theta \right)\right]\cdot \left[P_{{\Sigma }_1}\left(\theta ,\varphi \right)+P_{{\Sigma }_2}\left(\theta ,\varphi \right)\right]{\cos }^{\frac{1}{2}}\theta \sin \theta d\varphi d\theta \\ \end{array}$$

Here $x\text{ = }r\sin \theta \cos \varphi$ , $\text{y = }r\sin \theta \sin \varphi$ , $z\text{ = }-r\cos \theta$ , (a Cartesian-to-spherical coordinate relation).

The total intensity can be obtained

$$I\text{ = }{I}_x+I_y+I_z\text{ = }{\left|E_x\right|}^2+{\left|E_y\right|}^2+{\left|E_z\right|}^2\text{ (12)}$$

Then the two-dimensional electric field over an oblique plane with a oblique angle α is

$$\begin{array}{l}E\left(x_{\alpha },y_{\alpha }\right)\text{ = }{\displaystyle \underset{{\Sigma }_1+{\Sigma }_2}{\iint }\exp {\left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]}^{}{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}\cdot \exp \left(im\varphi \right)\left[\begin{array}{c}{\sin }^{\text{2}}\varphi +\cos \theta {\cos }^{\text{2}}\varphi \\ -\sin \varphi \cos \varphi +\cos \theta \sin \varphi \cos \varphi \\ \cos \varphi \sin \theta \end{array}\right]\text{ (13)}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}\times \exp \left\{-ik\left[x_{\alpha }\sin \theta \cos \varphi \right.\right.\left.\left.+y_{\alpha }\left(cos\alpha \sin \theta \sin \varphi -\sin \alpha \cos \theta \right)\right]\right\}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}\times \left(P_{{\Sigma }_1}\left(\theta ,\varphi \right)+P_{{\Sigma }_2}\left(\theta ,\varphi \right)\right){\cos }^{\frac{1}{2}}\theta \sin \theta d\varphi d\theta \end{array}$$

For the circular symmetry region ${\sum }_1$, the double integration is changed into a single integral. For , integration can be reduced by half. These considerations lead to

$$I\left(x_{\alpha },y_{\alpha }\right)\text{ = }{\left|I_a+I_b\right|}^2\text{ (14)}$$

$$\begin{array}{l}{I}_a\text{ = }{\displaystyle {\int }_{{\theta }_c}^{{\theta }_{\max }}{\displaystyle {\int }_{{\varphi }_1\left(\theta \right)}^{\frac{\pi }{2}}\exp {\left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]}^{}{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}}\text{ (15)}\\ \begin{array}{cc}& \end{array}\times \exp \left(im\varphi \right)\left[\begin{array}{c}2\left({\sin }^2\varphi +\cos \theta {\cos }^2\varphi \right)\cos \left(k{x}_{\alpha }\sin \theta \cos \varphi \right)\\ -2i\left(\cos \theta -1\right)\sin \varphi \cos \varphi \sin \left(k{x}_{\alpha }\sin \theta \cos \varphi \right)\\ -2i\sin \theta \cos \varphi \sin \left(k{x}_{\alpha }\sin \theta \cos \varphi \right)\end{array}\right]\\ \begin{array}{cc}& \end{array}\times \exp \left[ik{y}_{\alpha }\left(\sin \alpha \cos \theta -\cos \alpha \sin \theta \sin \varphi \right)\right]{\cos }^{\frac{1}{2}}\theta \sin \theta d\varphi d\theta \end{array}$$

$$I_b\text{ = }\frac{1}{2}\cdot \left[\begin{array}{c}2\pi i^m\exp \left(im{\varphi }_p\right)\cdot I_m^n-\pi i^{m+2}\exp \left[i\left(m+2\right){\varphi }_p\right]\cdot I_{m+2}^n-\pi i^{m-2}\exp \left[i\left(m-2\right){\varphi }_p\right]\cdot I_{m-2}^n\\ -i\pi i^{m-2}\exp \left[i\left(m-2\right){\varphi }_p\right]\cdot I_{m-2}^n+i\pi i^{m+2}\exp \left[i\left(m+2\right){\varphi }_p\right]\cdot I_{m+2}^n\\ -2\pi i^{m+1}\exp \left[i\left(m+1\right){\varphi }_p\right]\cdot I_{m+1}^n-2\pi i^{m-1}\exp \left[i\left(m-1\right){\varphi }_p\right]\cdot I_{m-1}^n\end{array}\right]\text{ (16)}$$

$$\begin{array}{l}{I}_m^n\left(x_{\alpha },y_{\alpha }\right)\text{ = }{\displaystyle {\int }_0^{{\theta }_c}\exp {\left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]}^{}{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}\text{ (17)}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}\times {\cos }^{\frac{1}{2}}\theta \sin \theta \left(1+\cos \theta \right)J_m\left(k\sqrt{x_{\alpha }^2+y_{\alpha }^2{\cos }^2\alpha }\sin \theta \right)\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}\times \exp \left(ik{y}_{\alpha }\sin \alpha \cos \theta \text{+}im\varphi \right)d\theta \end{array}$$

$$\begin{array}{l}{I}_{m\pm 1}^n\left(x_{\alpha },y_{\alpha }\right)\text{ = }{\displaystyle {\int }_0^{{\theta }_c}\exp {\left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]}^{}{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}\text{ (18)}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}& \end{array}\times {\cos }^{\frac{1}{2}}\theta {\sin }^2\theta J_{m\pm 1}\left(k\sqrt{x_{\alpha }^2+y_{\alpha }^2{\cos }^2\alpha }\sin \theta \right)\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}& \end{array}\times \exp \left(ik{y}_{\alpha }\sin \alpha \cos \theta \text{+}im\varphi \right)d\theta \end{array}$$

$$\begin{array}{l}{I}_{m\pm 2}^n\left(x_{\alpha },y_{\alpha }\right)\text{ = }{\displaystyle {\int }_0^{{\theta }_c}\exp {\left[\text{-}{\left(\frac{\rho }{w_0}\right)}^2\right]}^{}{\left(\sqrt{\text{2}}\frac{\rho }{w_0}\right)}^m{L}_{\frac{n-\left|m\right|}{2}}^{\left|m\right|}\left[2{\left(\frac{\rho }{w_0}\right)}^2\right]}\text{ (19)}\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}& \end{array}\times {\cos }^{\frac{1}{2}}\theta \sin \theta \left(1-\cos \theta \right)J_{m\pm 2}\left(k\sqrt{x_{\alpha }^2+y_{\alpha }^2{\cos }^2\alpha }\sin \theta \right)\\ \begin{array}{cc}\begin{array}{cc}\begin{array}{cc}\begin{array}{cc}& \end{array}& \end{array}& \end{array}& \end{array}\times \exp \left(ik{y}_{\alpha }\sin \alpha \cos \theta \text{+}im\varphi \right)d\theta \end{array}$$

The time-averaged electrical energy density or intensity can be obtained by the modulus squared of the electric fields. The total intensity distribution at the focal plane is the sum of the longitudinal and transverse components.

Numerical Simulation Results

In order to probe further the field structure at the focal plane, we have performed numerical calculations for the intensity distributions. Assuming that the light is an x-axis polarized Laguerre-Gaussian beam (the topological charge m = 1, n = 1), the numerical aperture of the lens NA = 1.4, the wavelength of the beam ${\lambda }_0\text{ = 519 nm}$, the diameter of optical aperture is 2.0 cm, the FWHM (full width at half maxima) of the laser beam is 0.5 cm.

In Figure 3a, Figure 3b, Figure 3c and Figure 3d, α = 0°, the intensity is equal to the conventional intensity of the circular aperture system, intensity profile of x component (I_x) has two peaks, intensity profile of y component (I_y) has eight peaks (four secondary peaks around four main peaks), I_y is about 1% of the total intensity, intensity of z components (I_z) is about one third of the total intensity, I_z has three peaks (two secondary peaks around a main peak), the x-direction intensity is the main part of the total intensity. As oblique angle α = 10°, the intensity distribution is distorted, two peaks of I_x become asymmetrical, most peaks of I_y degenerate, Intensity I_z decreases slightly, and the distortion is bigger when the oblique angle α = 20°.

Figure 4a, Figure 4b, Figure 4c, Figure 4d, Figure 4e, Figure 4f, Figure 4g, Figure 4h, Figure 4i, Figure 4j, Figure 4k and Figure 4l shows the oblique incident beam with different the topological charge (m = 1, n = 3; m = 3, n = 3; m = 3, n = 5), the beam incident angle α = 10°, NA = 1.4, other parameters are same as which in Figure 3. The intensity distribution of first mode (the topological charge m = 1, n = 3) is shown in Figure 4a, Figure 4b, Figure 4c and Figure 4d, The intensity I_x is about 75% of total intensity I. The intensity I_y is about 0.1% of total intensity I, there are four main peaks at the center and sub-peaks around. The intensity I_z is about 25% of total intensity I, there is one main intensity peak surrounded by two sub-peaks. The second modes (the topological charge m = 3, n = 3) are shown in Figure 4e, Figure 4f and Figure 4g, I_x is nearly 80% of total intensity I, the intensity I_x distribute in a circle, have eight main lobes, upper (according to y axis direction) two lobes are overlapped. I_y is about 2% of total intensity I, The intensity distribution have ten main lobes. For I_z, is about 20% of total intensity I, some of lobes are overlapped. The intensity distributions of third mode (the topological charge m = 3, n = 5) are shown in Figure 4i, Figure 4j, Figure 4k and Figure 4l, which is similar to that of the second mode. For all these modes, the structure of the field intensity near the focal area significantly changes to an asymmetrical shape, and distribution of I is mainly depended on I_x.

Conclusions

We analyzed focusing of an oblique incident Laguerre-Gaussian vector beam by a high numerical-aperture lens. Based on the Debye-Wolf vector diffraction theory, we formulate the electrical field distribution. Numerical simulation has also been performed, three-dimension intensity distributions of an oblique incident Laguerre-Gaussian beam on the focal plane of the lens were gotten, and the total intensity distribution on the focal plane, which is the sum of the longitudinal and transverse components, were analyzed in detail, the intensity profile of the focused Laguerre-Gaussian beam has several peaks related with topological charge. For normal incident beams, intensity profiles are symmetric, when the beam is obliquedly incident on the lens, value of some intensity peaks increase, as value of others decrease, the distribution profile is asymmetric. The distortion of light beams intensity distribution get greater with the increasing oblique incident angle. We also analyze the distortions of light beams intensity distribution with different topological charge. In practical application, it is an important reference for trapping and manipulating particles, super-resolution fluorescence microscopes.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 61975158), State Key Laboratory fund (The State Key Laboratory of Transient Optics and Photonics，Xi'an Institute of Optics and Precision Mechanics).

Data Availability Statement

Data underlying the results presented in this paper can be obtained from the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

---