# Analyzing the Symmetry Properties of a Distribution in the Focal Plane for a Focusing Element with Periodic Angle Dependence of Phase

Authors: Svetlana N. Khonina; Andrey V. Ustinov
Journal: Advances in Optical Technologies
 (2012)

---

## Abstract

We analyze the symmetry properties of the focal plane distribution when light is focused with an element characterized by a periodic angular dependent phase, sin (mφ) or cos (mφ). The majority of wave aberrations can be described using the said phase function. The focal distribution is analytically shown to be a real function at odd values of m, which provides a simple technique for generating designed wave aberrations by means of binary diffractive optical elements. Such a possibility may prove useful in tight focusing, as the presence of definite wave aberrations allows the focal spot size to be decreased.
The analytical computations are illustrated by the numerical simulation, which shows that by varying the radial parameters the focal spot configuration can be varied, whereas the central part symmetry is mainly determined by the parity of m: for even the symmetry order is 2m and for odd is m.


---

## Body

## 1. Introduction

Various aberrations in the focusing system are known to result in a wider, distorted focal spot with disturbed axial symmetry [1]. Such an effect is normally considered to be a negative factor.However, it has been shown [2, 3] that some types of wave aberrations enable the central focal spot size to be decreased, providing tight focusing. Note that while only primary (axisymmetric) aberrations were dealt with in [2], aberrations associated with vortex phase components on the basis of Zernike polynomials were also discussed in [3].In Zernike polynomials, the radius dependence is polynomial and the angle dependence is trigonometric (periodic). Optical elements characterized by periodic angular changes were considered in [4, 5]. In [4], such an element was shown to form the zero central intensity, whereas [5] also looked into diffraction-free properties of the generated light beams.Based on the decomposition of a cosine angular dependent phase function in terms of angular harmonics, the transmission functionexp{iacos(mφ)} was shown [5] to produce a diffraction pattern composed of 2m light spots arranged on a circumference. The coaxial interference of two vortex beams with identical topological charges and opposite signs was shown to produce a similar result [6, 7].At the same time, the odd-order aberrations, such as distortion and coma, have been known to appear in distribution patterns with odd symmetry [1, 8, 9]. In particular, the presence of coma (m=3) results in distributions with the third-order symmetry [10], similar to the 2D Airy beams [11].It was also shown that the product of three 1D Airy functions, rotated by the angle of 120° relative to each other and characterized by the third-order symmetry, was transformed in the spectral plane into a function proportional toexp{iar3sin(3φ)} [12].Notice, the phase mask with sine but nonperiodic angular dependence [13] breaking symmetry can be used for selective edge enhancement.In this paper, we analyze symmetry properties of the distribution formed in the focal plane by an optical element characterized by a periodic angular dependent phase function in the form ofsin(mφ) or cos(mφ). Based on this phase relation, it is possible to describe the majority of wave aberrations, which then can be represented [14] as the decomposition in terms of Zernike functions [1].The focal distribution is analytically shown to be described by a real function at odd values ofm. Analytical estimates of the central part of the focal distribution are derived for an input circular aperture.The analytical computations are illustrated by the numerical simulation, which shows that by varying the radial parameters the focal spot configuration can be varied, whereas the central part symmetry is manily determined by the even value of the angular parameterm.


## 2. Fourier Transform of a Complex Distribution with Periodic Angular Dependence

Assume a complex distribution with periodic angular dependence of the general form:(1)f(r,φ)=ψ(r)exp[-iβ(αr)q{sinmφcosmφ}],
where ψ(r) is an axisymmetric function, m is integer, q is a positive real number, and α, β are real parameters (with α having the dimension of mm−1).Note that various combinations of the functionssin(mφ) and cos(mφ) in the phase function of (1) can be represented as a product of functions given by (1).The complex distribution in (1) can describe the majority of wave aberrations, which, then, can be represented [14] by the decomposition in terms of Zernike functions [1].For the sake of specificity, choosesin(mφ) in (1). Note that for the cos(mφ) function the results will be analogous.The spatial Fourier spectrum for the function (1) is given by
(2)F(ρ,θ)=-ikexp(ikf)2πf×∫0∞ψ(r){∫02πexp[-iβ(αr)qsinmφ]×exp[-ikfrρcos(θ-φ)]dφ∫02πexp[-iβ(αr)qsinmφ]}rdr,
where k=2π/λ is the wave number, λ is the wavelength of laser light, f is the focal length, (r,φ) and (ρ,θ) are the polar coordinates  in the input and output planes, respectively.Let us analyze the integral over the angle in braces in (2):
(3)Φ(ρ,θ,r)=∫02πexp[-iβ(αr)qsinmφ]exp[-ikfrρcos(θ-φ)]dφ.The axial distribution of light (ρ=0) is described by the relation independent of m:
(4)Φ(0,0;r)=∫02πexp[-iβ(αr)qsinmφ]dφ=2πJ0[β(αr)q].
Designating
(5)a=β(αr)q,b=krρf
and following the transformation in (3) at m=1(6)asinφ+bcos(φ-θ)=(a+bsinθ)sinφ+(bcosθ)cosφ=a2+b2+2absinθcos(φ-θ′),
where θ′=arctg[(a+bsinθ)/bcosθ], we obtain
(7)Φ(ρ,θ;r)=2πJ0(β2(αr)2q+(krρf)2+2β(αr)q(krρf)sinθ).
Replacing sinφ in (1) by cosφ (at m=1) yields an analogous relation to (7), following the relevant substitutions.Let us consider a more general case of integerm. The first term in (3) can be decomposed as
(8)exp(-iasinmφ)=J0(a)+2∑p=1∞J2p(a)cos(2pmφ)-2i∑p=0∞J2p+1(a)sin((2p+1)mφ).
Equation (8) employs the following relations [15]:(9a)cos(zsint)=J0(z)+2∑p=1∞J2p(z)cos(2pt),(9b)sin(zsint)=2∑p=1∞J2p+1(z)sin[(2p+1)t].To calculate the integral in (3), we can employ the well-known relations:
(10)∫02π{cosmφsinmφ}exp[-ibcos(θ-φ)]dφ=2π(-i)m{cosmθsinmθ}Jm(b).
Thus, we obtain
(11)Φ(ρ,θ;r)=2πJ0(a)J0(b)+4π∑p=1∞(-1)pmJ2p(a)J2pm(b)cos(2pmθ)-4πi∑p=0∞(-1)pm(-i)mJ2p+1(a)J(2p+1)m(b)×sin((2p+1)mθ).
It should be noted that for the odd values of m, (11) is a real function, because the coefficients in the third sum i(-i)m=(-i)m-1 are real at odd m.Takingcosmφ in (1), we obtain an analogous relation:
(12)exp(-iacosmφ)=J0(a)+2∑p=1∞(-1)pJ2p(a)cos(2pmφ)-2i∑p=0∞(-1)pJ2p+1(a)cos((2p+1)mφ),
which utilizes the relations [15]: (13a)cos(zcost)=J0(z)+2∑p=1∞(-1)pJ2p(z)cos(2pt),(13b)sin(zcost)=2∑p=1∞(-1)pJ2p+1(z)cos[(2p+1)t].Thus,
(14)∫02πexp(-iacosmφ)exp[-ibcos(θ-φ)]dφ=2πJ0(a)J0(b)+4π∑p=1∞(-1)p(m+1)J2p(a)J2pm(b)cos(2pmθ)-4πi∑p=0∞(-1)p(m+1)(-i)mJ2p+1(a)J(2p+1)m(b)×cos((2p+1)mθ).
For the odd values of m, the function in (14) is real.Substitute (11) in (2):
(15)F(ρ,θ)=-ikfexp(ikf)×∫0∞ψ(r){J0[β(αr)q]J0(kfrρ)-2i(-i)mJ1[β(αr)q]Jm×(kfrρ)sin(mθ)+2(-1)mJ2[β(αr)q]J2m×(kfrρ)cos(2mθ)-⋯}rdr.
Equation (15) can be rearranged as follows:
(16)F(ρ,θ)=-ikfexp(ikf){Ψ0(ρ)-2(-i)m-1Ψ1(ρ)sin(mθ)+2(-1)mΨ2(ρ)cos(2mθ)-⋯(-i)m-1},
where
(17)Ψp(ρ)=∫0∞ψ(r)Jp[β(αr)q]Jpm(kfrρ)rdr.Let us obtain approximate estimates of the integrals in (17), putting
(18)ψ(r)={1,r≤R,0,r>R,
and assuming J0(z)≈1-z2/4,  Jn(z)≈(z/2)n/n!.Then,(19a)Ψ0(ρ)=∫0RJ0[β(αr)q]J0(kfrρ)rdr≈∫0R[1-β2(αr)2q4-(krρ)24f2]rdr=R22-(βR)2(αR)2q8(q+1)-(kρR24f)2,(19b)Ψp(ρ)=∫0RJp[β(αr)q]Jpm(kfrρ)rdr≈1p!(pm)!∫0R[β(αr)q2]p(k2frρ)pmrdr=βpαpq(kρ)pmp!(pm)!2p(2f)pm·Rp(q+m)+2p(q+m)+2,p>0.It can be seen from (19a) and (19b) that the major contribution to the central part of the focal plane (at ρ→0) is provided by the terms with lower-order p. Truncating (16) at the first two terms, the intensity in the central part of the focal plane can be estimated for the even values of  m as(20a)|F(ρ,θ)|2=(kf)2[Ψ02(ρ)+4Ψ12(ρ)sin2(mθ)],
and for odd m as
(20b)|F(ρ,θ)|2=(kf)2[4(-1)((m-1))/2Ψ02(ρ)+4Ψ12(ρ)sin2(mθ)-4(-1)(m-1)/2Ψ0(ρ)Ψ1(ρ)sin(mθ)],with the first term being independent of the angle and the second term resulting in 2m uniformly distributed intensity peaks of equal height (at sin(mθ)=±1). The third term is only present for odd m, introducing corrections into the peaks: depending on the sign of the term containing sin(mθ), half of the peaks are enhanced, while the other half are suppressed, which results in m clearly expressed maximums.Thus, (20a) and (20b) suggest that the symmetry of the central part of the focal distribution is determined by the parity of m: for even m, the symmetry order is 2m, for odd m, the symmetry order is m.By varying the parametersβ, α, q, and m, the values of the coefficients in (20a) and (20b) can be varied, thus changing the symmetry of the intensity pattern in the focus.


## 3. Numerical Simulation

The numerical simulation based on (2) was conducted for the following parameters: wavelength of incident light, λ=633 nm, focal length, f=100 mm.Table1 gives the simulation results based on (1), for sin(mφ), for different values of mand different parameter’s values of α, β, and q, with the ψ(r) function being given by (18), describing a circular diaphragm of radius R=1 mm, and a Gaussian beam exp(-r2/σ2) of waist radius σ=0.5 mm.Table 1
Simulated focusing results for simple functionsψ(r) and different parameters.
Parameters
Input (phase) distribution
Light distribution in the focal plane  (intensity and phase),0.5×0.5 mm
m
=
1, ψ(r) from (18),q=0,  α=1 mm−1, β=27
m
=
2, ψ(r) from (18),q=0,  α=1 mm−1, β=5
m
=
3, ψ(r) is the Gaussian beam,q=0, α=1 mm−1, β=3
m
=
3, ψ(r) is the Gaussian beam,q=3, α=3 mm−1, β=1
m
=
4, ψ(r) from (18),q=0,  α=1 mm−1, β=3Notice that the input phase distributions (shown in the second column of the table) have the same values in opposite points on diameter line for evenm and negative values for odd m. It is easy to confirm by sin(2(x+π))=sin(2x), sin(x+π)=-sin(x). Both of cases differ from a situation with a spiral phase exp(imφ) where a phase difference at a symmetric position is mπ [13], so the spiral phase is used for the radial Hilbert transform of mth order. In our case a phase difference in symmetric points can explain focal distribution. Since at odd m in focus there are identical phases with opposite sign the field will be real valued (exp(iψ)+exp(-iψ)~cos(ψ)). For even m, identical values are summarized, therefore the phase structure is regenerated in focal domain.Table2 presents the focusing results for q=0 and β=5, with the values of m and ψ(r) being varied.Table 2
Simulation results for different parameters and the form of theψ(r,φ) function, q=0.
Parameters
Input distribution (intensity and phase)
Focal distribution (intensity and phase),0.5×0.5 mm
m
=
1,   ψ(r) describes the Laguerre-Gaussian mode (3,0), σ=0.2 mm
m
=
3,  ψ(r) describes the Laguerre-Gaussian mode (3,0), σ=0.2 mm
m
=
1, ψ(r,φ) describes the Laguerre-Gaussian mode (3,1), σ=0.2 mm
m
=
3, ψ(r,φ) describes the Laguerre-Gaussian mode (3,3), σ=0.2 mm
m
=
3, ψ(r)=J0(γr), γ=10 MM−1
m
=
1, ψ(r,φ)=J1(γr)exp(iφ), γ=10 MM−1The effect of replacing the axisymmetric functionψ(r) with a more general function ψ(r,φ) with a harmonic angle dependence, exp(ilφ), was also studied.The four top lines are for the functionψ(r,φ) described by the Fourier-transform-invariant Laguerre-Gaussian modes [15] of different orders (n,l), including those with a vortex function exp(ilφ). It is noteworthy that the vortex angle component does not exert such an essential influence on the focal pattern structure as the periodic function in question sin(mφ) does. The comparison of lines 1 and 3 shows that a phase vortex introduced in the input plane results in a similar pattern in the focal plane, without changing the major distribution structure (except for the zero intensity in the singular point).Whenm=1, the resulting distribution is similar to that associated with the coma-type aberration introduced in the focusing system; when m=3, the distribution is described by the product of three 1D Airy functions, rotated relative to each other [12]. Note, however, that in the latter case the product is actually composed of somewhat different functions.In two bottom lines of Table2, the function ψ(r,φ) is given by the Bessel modes Jl(γr)exp(ilφ), which produce a Fourier spectrum in the form of a narrow ring  of radius proportional to γ. The simulation results show that the annular focal structure produced by the Bessel modes undergoes essential variations caused by the periodic angle change.


## 4. Conclusions

We have analyzed symmetry properties of the light distribution generated in the focal plane by focusing the light beam with a periodic angular dependent phase given bysin(mφ) or cos(mφ).The focal distribution of light is analytically shown to be a real function when the value ofm is odd. Thus, in a way similar to [12], phase distributions associated with the corresponding specified aberration types can be formed in the focal plane of a conventional spherical lens with the aid of binary diffractive optical elements. Such a possibility may prove useful in tight focusing, where the presence of certain wave aberrations enables the focal spot size to be reduced [2, 3]. In particular, for a circularly polarized beam the spot size can be reduced by the presence of coma (m=3), corresponding to the 2D superpositions of Airy beams [11, 12].The analytical computations have been illustrated by the numerical simulation, which has shown that by varying the radial parameters, the focal distribution configuration can be varied; meanwhile, the central part symmetry is mainly determined by the parity of the angle parameterm: when m is even, the central part has the 2m-order symmetry, and when m is odd—the symmetry is m order. Note that although the introduced vortex component exp(ilφ) is preserved in the focal plane, it does not have an essential effect on the focal distribution structure when compared with the periodic function in question sin(mφ).



---

*Source: 918298-2012-10-10.xml*