Digital Sorting of Optical Vortices in Perturbed Singular Beams

*Alexander Volyar, Mikhail Bretsko, Yana Akimova and Yuriy Egorov*

## **Abstract**

The chapter provides a brief overview of shaping and measuring techniques of the vortex spectra (squared amplitudes and initial phases of vortex modes) including radial indices. The main physical mechanisms causing the formation of laser beams with a complex vortex composition, in particular, in biological media, are indicated, and the need for a digital analysis of vortex spectra is substantiated. It is the analysis of vortex spectra that allows us to find the orbital angular momentum and informational entropy (Shannon's entropy) of perturbed laser beams in real time. In the main part of the chapter, we consider in detail a new approach for measuring vortex spectra without cuts and gluing of the wavefront, based on digital analyzing high-order intensity moments of complex beams and sorting the vortex beam in computer memory sells. It is shown that certain types of weak local inhomogeneities cause a vortex avalanche causing a sharp dips and bursts of the orbital angular momentum spectra and quick ups and downs of the informational entropy. An important object of analysis is also the vortex spectra of beams scattered by simple opaque obstacles such as a hole, a disk, and a sector aperture.

**Keywords:** optical vortex, moment's intensity, orbital angular momentum, medical optics

## **1. Introduction**

As is well known, optical vortices [1] accompany light scattering processes due to both simple and complex medium inhomogeneities (see, e.g., [2] and references therein). Scattered light carries a huge array of information, both on the composition and structure of scattering (diffraction) centers, and on the structure of the initial light beam. At the same time, optical information can be read off both by analyzing the frequency spectrum [3–5] and the spectrum of optical vortices [6, 7]. An important aspect of this problem is the study of biomedical objects [8, 9], for example, the composition of the blood or the reflection of a vortex beam from the skin surface [10] for express diagnostics of skin diseases. The fact is that a vortex beam scattered by the skin is transformed into a speckle-like structure resembling one that occurs when light passes through a medium with weak turbulence [11, 12]. At the same time, the speckle structure is formed by the skeleton of optical vortex array [13–15]. In turn, the analysis of such a complex vortex structure is conveniently carried out on the basis of vortex fractal techniques [16] (see also [17] and references therein). There are a variety of approaches for the fractal vortex models

of laser beams scattered by biological tissues [18, 19] based on the light scattering by nonspherical particles [20] that involves the representation of the wave field in terms of Legendre polynomials. However, a real optical experiment for measuring the vortex spectrum requires the use of particular approaches to the orthogonal basis for the scattered beam expansion in terms of special functions, which can differ significantly from the corresponding theoretical models.

The problem of using the properties of optical vortices in various areas of science and technology requires the development of reliable but simple techniques for measuring the spectrum of optical vortices in complex beams scattered by various objects. Therefore, the attention of many researchers is drawn to measuring the orbital angular momentum (OAM) of vortex beams [7, 21–30] that is directly related to the spectrum of optical vortices. Really, the complex amplitude Ψð Þ *x*, *y*, *z* of a composite paraxial beam can be written in the form of a superposition of the orthogonal vortex modes *ψ m*ð Þ *x*, *y*, *z*

$$\Psi(\mathbf{x}, \mathbf{y}, \mathbf{z}) = \sum\_{m=-N}^{N} \mathbf{C}\_{m} \Psi\_{m}(\mathbf{x}, \mathbf{y}, \mathbf{z}),\tag{1}$$

of the combined beam as a whole at the focal plane of the cylindrical lens without cutting and gluing the wavefront while the first two techniques distorts radically the initial beam structure. In the technique [26], the measurement of second-order intensity matrix elements was used that enables the authors to devoid breaking down the beam structure. However, it makes sense to use such an approach only if the modes of a complex beam have axial symmetry. The authors of [34] improved the technique by using two cylindrical lenses whose axes are rotated by the angle *π=*2 relative to each other. This made it possible to measure the fractional OAM for any

At the same time, the measurement of the fractional OAM does not give complete information about the combined beam content. As can be seen from expression (2), the same value of the OAM can correspond to a different composition of

medium [37] or optical fibers [38], there is an intense energy exchange between the vortex modes of the singular beam. It makes to change both the magnitude and the spectrum of the vortex modes including their initial phases. To measure the vortex spectrum, a special technique was developed based on the analysis of higher-order intensity moments [25] at the focal plane of a spherical lens. A characteristic feature of this technique is the ability to measure not only amplitudes but also the initial phases of the vortex modes. The point is that the intensity moments at the focal plane of the spherical lens are degenerate with respect to the sign of the vortex topological charge. Therefore, the authors of [25] note that the technique is applicable only for nondegenerate spectra of combined singular beams. In [39], the authors expanded the possibilities of measuring the vortex spectrum, also covering

optical vortices with different signs of topological charges due to additional

transformations of the Laguerre-Gauss beams (LG) into Hermite-Gauss ones (HG) via astigmatic conversions at a cylindrical lens [40]. However, the method of measuring optical vortices in LG beams of higher orders remained uncovered. In this chapter, we consider in detail the technique of measuring the vortex spectrum based on the analysis of high-order intensity moments that excludes cuts and gluing of the beam without losses of information on initial mode phases. Unlike the method of holographic gratings that transforms a combined beam into vortex modes with different propagation directions, just as a prism converts "white" light into a spatial frequency spectrum, we will try to demonstrate how a perturbed singular beam can be sorted into vortex modes sited in computer memory cells and then reproduce the beam main characteristics: OAM, information entropy, and initial topological charge in real time. Moreover, knowing the digital spectrum, the beam can be recovered again, and by adjusting the spectral vortex amplitudes we

As far as we know, the first theoretical and experimental studies of optical vortex arrays refer to 1991, when the authors of [41] succeeded in reproducing holographically individual letters and words due to ordering optical vortex array in typical phase skeleton on the base of the technique that had been developed back in the early 1980s [29]. However, only after the article [42] by Berry did the studies of the vortex array properties become widespread. As a result, it was shown that the diffracted beam turns into a combined beam containing a large number of optical vortices with integer topological charges. These are sometimes called the beams with fractional topological charges. Using Eq. (2), we can verify that beams with a fractional topological charge also have fractional OAM ℓ*z*. However, in the process

*m*. For example, in a weakly turbulent

type of the mode beam symmetry.

vortex modes with squared amplitudes *C*2

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

can improve the structure of the transmitted field.

**2. Preliminary remarks**

**165**

where *Cm* ¼ h j Ψ *ψ m*i stands for the mode amplitudes specified by the normalizing condition P*N m*¼�*N*j j *Cm* 2 ¼ 1, and 2*N* is a total number of vortex modes, Then the average on photon OAM of the composite beam is found by the formula [31].

$$\mathcal{E}\_z = \sum\_{m=-N}^{N} m \left| \mathbf{C}\_m \right|^2. \tag{2}$$

The basis of this expression is an implicit connection between the topological charge *m* and the OAM ℓ*z m* of a single mode in the beam. However, expression (2) includes the squared modulus j j *Cm* 2 of a complex number while we are talking about measuring the complex number *Cm* itself, i.e., its amplitude and the initial phase. The known for us approaches to measuring OAM sorting modes can be divided into four groups according to the methods used. All known for us approaches to measuring OAM spectrum and mode sorting can be conditionally divided into four groups according to the methods used and the effect on the internal beam structure. The most widely used method of holographic spatial filtering of the Laguerre-Gaussian, Hermite-Gaussian, Bessel-Gaussian, Airy-Gaussian and other mode types [22–24, 30]. At that, this method allows modes sorting, both by topological charge (azimuthal numbers *m*) and by radial numbers *n*. The holographic mode sorting resembles the effect of a "white light" decomposition into a color spectrum by means of a prism or diffraction grating. The second type of measurements requires the use of a so-called log-polar transparency and diaphragms [6, 7, 24] placed along the path of a complex beam. Such a structured transparency converts the beam into a set of horizontal or vertical fringes, the pattern analysis of which allows one to obtain the OAM spectrum. The direct application of digital processing of interference patterns of a composite vortex beam [32] and a collinear phase-shifting holography [33] can be considered as the third interference technique. However, measurements in the listed above approaches lead to complete or partial damage of the combined beam and losses of useful information. At the same time, an original approach presented in the papers [26, 27, 34] to measuring the OAM in the vortex beams with a fractional topological charge enables one to avoid the beam damage in the result of employing cylindrical lenses and analyzing the second order intensity moments [35]. Unfortunately, such a simple approach is not applicable for both sorting and measuring the spectrum of vortex modes. In [26, 27], the technique analyzes a fractional OAM [36]

### *Digital Sorting of Optical Vortices in Perturbed Singular Beams DOI: http://dx.doi.org/10.5772/intechopen.91419*

of laser beams scattered by biological tissues [18, 19] based on the light scattering by nonspherical particles [20] that involves the representation of the wave field in terms of Legendre polynomials. However, a real optical experiment for measuring the vortex spectrum requires the use of particular approaches to the orthogonal basis for the scattered beam expansion in terms of special functions, which can

The problem of using the properties of optical vortices in various areas of science

and technology requires the development of reliable but simple techniques for measuring the spectrum of optical vortices in complex beams scattered by various objects. Therefore, the attention of many researchers is drawn to measuring the orbital angular momentum (OAM) of vortex beams [7, 21–30] that is directly related to the spectrum of optical vortices. Really, the complex amplitude Ψð Þ *x*, *y*, *z* of a composite paraxial beam can be written in the form of a superposition of the

*N*

*m*¼�*N*

average on photon OAM of the composite beam is found by the formula [31].

*m*¼�*N*

The basis of this expression is an implicit connection between the topological charge *m* and the OAM ℓ*z m* of a single mode in the beam. However, expression (2) includes the squared modulus j j *Cm* 2 of a complex number while we are talking about measuring the complex number *Cm* itself, i.e., its amplitude and the initial phase. The known for us approaches to measuring OAM sorting modes can be divided into four groups according to the methods used. All known for us approaches to measuring OAM spectrum and mode sorting can be conditionally divided into four groups according to the methods used and the effect on the internal beam structure. The most widely used method of holographic spatial filtering of the Laguerre-Gaussian, Hermite-Gaussian, Bessel-Gaussian, Airy-Gaussian and other mode types [22–24, 30]. At that, this method allows modes sorting, both by topological charge (azimuthal numbers *m*) and by radial numbers *n*. The holographic mode sorting resembles the effect of a "white light" decomposition into a color spectrum by means of a prism or diffraction grating. The second type of measurements requires the use of a so-called log-polar transparency and diaphragms [6, 7, 24] placed along the path of a complex beam. Such a structured transparency converts the beam into a set of horizontal or vertical fringes, the pattern analysis of which allows one to obtain the OAM spectrum. The direct application of digital processing of interference patterns of a composite vortex beam [32] and a collinear phase-shifting holography [33] can be considered as the third interference technique. However, measurements in the listed above approaches lead to complete or partial damage of the combined beam and losses of useful information. At the same time, an original approach presented in the papers [26, 27, 34] to measuring the OAM in the vortex beams with a fractional topological charge enables one to avoid the beam damage in the result of employing cylindrical lenses and analyzing the second order intensity moments [35]. Unfortunately, such a simple approach is not applicable for both sorting and measuring the spectrum of vortex modes. In [26, 27], the technique analyzes a fractional OAM [36]

ℓ*z* ¼ X *N*

where *Cm* ¼ h j Ψ *ψ m*i stands for the mode amplitudes specified by the normaliz-

*m*¼�*N*j j *Cm* 2 ¼ 1, and 2*N* is a total number of vortex modes, Then the

*m C*j j *m* 2

*Cm*Ψ*m*ð Þ *x:y*, *z* , (1)

*:* (2)

differ significantly from the corresponding theoretical models.

Ψð Þ¼ *x*, *y*, *z* X

orthogonal vortex modes *ψ m*ð Þ *x*, *y*, *z*

*Nonlinear Optics - From Solitons to Similaritons*

ing condition P*N*

**164**

of the combined beam as a whole at the focal plane of the cylindrical lens without cutting and gluing the wavefront while the first two techniques distorts radically the initial beam structure. In the technique [26], the measurement of second-order intensity matrix elements was used that enables the authors to devoid breaking down the beam structure. However, it makes sense to use such an approach only if the modes of a complex beam have axial symmetry. The authors of [34] improved the technique by using two cylindrical lenses whose axes are rotated by the angle *π=*2 relative to each other. This made it possible to measure the fractional OAM for any type of the mode beam symmetry.

At the same time, the measurement of the fractional OAM does not give complete information about the combined beam content. As can be seen from expression (2), the same value of the OAM can correspond to a different composition of vortex modes with squared amplitudes *C*2 *m*. For example, in a weakly turbulent medium [37] or optical fibers [38], there is an intense energy exchange between the vortex modes of the singular beam. It makes to change both the magnitude and the spectrum of the vortex modes including their initial phases. To measure the vortex spectrum, a special technique was developed based on the analysis of higher-order intensity moments [25] at the focal plane of a spherical lens. A characteristic feature of this technique is the ability to measure not only amplitudes but also the initial phases of the vortex modes. The point is that the intensity moments at the focal plane of the spherical lens are degenerate with respect to the sign of the vortex topological charge. Therefore, the authors of [25] note that the technique is applicable only for nondegenerate spectra of combined singular beams. In [39], the authors expanded the possibilities of measuring the vortex spectrum, also covering optical vortices with different signs of topological charges due to additional transformations of the Laguerre-Gauss beams (LG) into Hermite-Gauss ones (HG) via astigmatic conversions at a cylindrical lens [40]. However, the method of measuring optical vortices in LG beams of higher orders remained uncovered.

In this chapter, we consider in detail the technique of measuring the vortex spectrum based on the analysis of high-order intensity moments that excludes cuts and gluing of the beam without losses of information on initial mode phases. Unlike the method of holographic gratings that transforms a combined beam into vortex modes with different propagation directions, just as a prism converts "white" light into a spatial frequency spectrum, we will try to demonstrate how a perturbed singular beam can be sorted into vortex modes sited in computer memory cells and then reproduce the beam main characteristics: OAM, information entropy, and initial topological charge in real time. Moreover, knowing the digital spectrum, the beam can be recovered again, and by adjusting the spectral vortex amplitudes we can improve the structure of the transmitted field.

## **2. Preliminary remarks**

As far as we know, the first theoretical and experimental studies of optical vortex arrays refer to 1991, when the authors of [41] succeeded in reproducing holographically individual letters and words due to ordering optical vortex array in typical phase skeleton on the base of the technique that had been developed back in the early 1980s [29]. However, only after the article [42] by Berry did the studies of the vortex array properties become widespread. As a result, it was shown that the diffracted beam turns into a combined beam containing a large number of optical vortices with integer topological charges. These are sometimes called the beams with fractional topological charges. Using Eq. (2), we can verify that beams with a fractional topological charge also have fractional OAM ℓ*z*. However, in the process

of shaping such beams, Eq. (2) requires careful use. Indeed, it is assumed that the functions Ψ*m* of mode beams in Eq. (2) are normalized; otherwise, misunderstandings arise. For example, in [43] considers a Bessel beam diffracted by a holographic grating with a fractional topological charge (see the callout in **Figure 1**).

If the Bessel beam is represented in the form of a conical beam of plane waves with fractional phase bypath *p*, then after the inverse Fourier transform we get the non-normalized superposition of Bessel-Gauss modes, including vortices with both positive and negative topological charges:

$$\Psi = 2\text{NG}e^{ip\pi}\sin\left(p\pi\right)\sum\_{m=-\infty}^{\infty}i^m f\_m(\text{KR})e^{im\phi}/(p-m),\tag{3}$$

choice of the normalization of vortex modes endows the combined beam with

For a simplest model for our studies, we choose a scalar field of wave in the form

*m* > 0 and *n* ¼ 0 are the azimuthal and radial indices, respectively, so that all the vortex modes in superposition have the same waist radius *w*0 at *z* ¼ 0. Let us

> *n*¼0 ¼ X *N*�1

For this, we make use of the intensity moments approach [35] in the form

*Mp*,*q*ð Þ *r*, *φ ℑ*ð Þ *r*, *φ r dr dφ=*

*N* X�1

*m*, *m*0 ¼0, *m* > *m*'

sin *m* � *m*0 ½ � ð Þ*φ* sin *βm*,*m*0*G r*ð Þ2

*Mp*,*q*ð Þ *r*, *φ* stands for the intensity moments function, *p*, *q* ¼ 0, 1, 2, … . Now the problem is to choose a combination of intensity moments *Jp*,*q* in such a way as to exclude all terms of the last two sums and leave only the first one in Eq. (7). The equations for the squared amplitudes and the initial phases are separated if the moment function is written as *Mp*ð Þ¼ *r*, *φ rp*. Then the equations for the squared

*m*¼0

amplitude, *βm* is the initial mode phase. We consider the nondegenerate case of the field representationð Þ *m* ≥0 or*m* ≤ 0 due to the axial symmetry; combined beams with different topological signs �*m* are indistinguishable in intensity moments. Our goal is to analyze the distribution of the wave field intensity *ℑ*ð Þ *r*, *φ*, *z* ¼ 0 in such a way as to express the squared amplitudes and the initial phases of the vortex modes in terms of physically measured quantities in the region of the beam waist.

*Cmr me*

*m*!*π* p stands for the normalization factor, *Cm* is the mode

ðð

*S*

*r m*þ*m*0

,

*Cm*0*Cm NmNm*0

*n* of the lowest order with the *N* monochromatic beams, where

*i m*ð Þ *φ*þ*βm G r*ð Þ , *z =Nm*, (6)

*ℑ*ð Þ *r*, *φ r dr dφ*, (7)

cos *m* � *m*0 ½ � ð Þ*φ* cos *βm*,*m*0*G r*ð Þ2

(8)

**3. Theoretical background of the digital vortex sorting and**

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

consider the wave field in the waist plane *z* ¼ 0 in the form [25].

*CmLGm*

*N* X�1 *m*¼0

special properties.

**experimental results**

Ψð Þ¼ *r*, *φ*, *z* ¼ 0

where *Nm* ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi

2�*m*�1

*Jp*,*q* ¼

*ℑ* ¼ Ψ∗ Ψ ¼ X

�2 X *N*�1

**167**

*m*, *m*0 ¼0, *m* > *m*'

*N*�1

*C*2 *m N*2 *m r*

> *r m*þ*m*0

*m*¼0

*Cm*0*Cm NmNm*0

amplitudes take the form

ðð

*S*

where the intensity distribution is written as

2*mG r*ð Þ2 þ 2

**3.1 Nondegenerate case**

of superposition *LGm*

where *Jm*ð Þ *x* stands for a Bessel function of the first kind and *m* order, *R* ¼ *r=w*0, *r* and *φ* are polar coordinates, *N* ¼ exp *iK*2 *=*2*kZ* � �, *G* ¼ exp *ikr*2 *=*2*Z* � �, *Z* ¼ *z* � *i z*0, *w*0 is a beam waist at *z* ¼ 0, *K* is a scale parameter, *k* is a wavenumber. After the corresponding transformations, we find the OAM of the perturbed beam

$$\mathcal{E}\_z = \sum\_{m = -\infty}^{\infty} m I\_m \left( |K|^2 w\_0^2 / 4 \right) / (p - m)^2 / \sum\_{m = -\infty}^{\infty} I\_m \left( |K|^2 w\_0^2 / 4 \right) / (p - m)^2. \tag{4}$$

The result of the plotting is represented by the curve 1 in **Figure 1(a)**. The OAM oscillates at large values of the topological charge *p*. The integer values of the topological charge *p* correspond to sharp bursts of the OAM ℓ*z*. However, a small deviation of the parameter *p* from the integer value causes sharp OAM dips. A completely different situation occurs if the elementary beams in Eq. (3) are normalized, as the authors of [44] do, the OAM oscillations disappear (see **Figure 1(b)** and curve 2 in **Figure 1(c)**). In this case, OAM ℓ*z*ð Þ *p* obeys a simple relation (see also [45])

$$\ell\_x = p - \sin\left(2\pi p\right)/\pi. \tag{5}$$

The OAM oscillations disappear (see **Figure 1(c)**, curve 2). As the authors of [45] revealed, a gradual increase in OAM is observed only at small values of the fractional topological charge. In fact, we are dealing with different beams, although the basis for their shaping is the same physical process. As we will see later, the

#### **Figure 1.**

*(a) Dependence of the OAM* ℓz *on the fractional topological charge* p*: curve 1 is calculated according to Eq. (3) with K* ¼ *2* � *104m*�*1, curve 2 presents Eq. (5), leader—diffraction grating for p* ¼ *5=2. (b and c) Computer modeling of the OAM curves ℓz*ð Þ *m related with combined beams in Eq. (4): (b) total value of the OAM curves at Ω* ¼ *10; (c) the characteristic spectral bursts of OAM; curve 1—the OAM bursts* ð Þ *Ω* ¼ *0:4 , curve 2—dips of the OAM* ð Þ *Ω* ¼ *14 .*

choice of the normalization of vortex modes endows the combined beam with special properties.

## **3. Theoretical background of the digital vortex sorting and experimental results**

#### **3.1 Nondegenerate case**

of shaping such beams, Eq. (2) requires careful use. Indeed, it is assumed that the functions Ψ*m* of mode beams in Eq. (2) are normalized; otherwise, misunderstandings arise. For example, in [43] considers a Bessel beam diffracted by a holographic

If the Bessel beam is represented in the form of a conical beam of plane waves with fractional phase bypath *p*, then after the inverse Fourier transform we get the non-normalized superposition of Bessel-Gauss modes, including vortices with both

> *m*¼�∞ *i*

*w*0 is a beam waist at *z* ¼ 0, *K* is a scale parameter, *k* is a wavenumber. After the corresponding transformations, we find the OAM of the perturbed beam

> 2 *=* X∞ *m*¼�∞

The result of the plotting is represented by the curve 1 in **Figure 1(a)**. The OAM

The OAM oscillations disappear (see **Figure 1(c)**, curve 2). As the authors of [45] revealed, a gradual increase in OAM is observed only at small values of the fractional topological charge. In fact, we are dealing with different beams, although the basis for their shaping is the same physical process. As we will see later, the

*(a) Dependence of the OAM* ℓz *on the fractional topological charge* p*: curve 1 is calculated according to Eq. (3) with K* ¼ *2* � *104m*�*1, curve 2 presents Eq. (5), leader—diffraction grating for p* ¼ *5=2. (b and c) Computer modeling of the OAM curves ℓz*ð Þ *m related with combined beams in Eq. (4): (b) total value of the OAM curves at Ω* ¼ *10; (c) the characteristic spectral bursts of OAM; curve 1—the OAM bursts* ð Þ *Ω* ¼ *0:4 , curve 2—dips*

*=*ð Þ *p* � *m*

oscillates at large values of the topological charge *p*. The integer values of the topological charge *p* correspond to sharp bursts of the OAM ℓ*z*. However, a small deviation of the parameter *p* from the integer value causes sharp OAM dips. A completely different situation occurs if the elementary beams in Eq. (3) are normalized, as the authors of [44] do, the OAM oscillations disappear (see **Figure 1(b)** and curve 2 in **Figure 1(c)**). In this case, OAM ℓ*z*ð Þ *p* obeys a simple relation (see

where *Jm*ð Þ *x* stands for a Bessel function of the first kind and *m* order, *R* ¼ *r=w*0,

*mJm*ð Þ *KR e*

*=*2*kZ* � �, *G* ¼ exp *ikr*2

*Im* j j *K* 2

*w*2 0*=*4 � �

ℓ*z* ¼ *p* � sin 2ð Þ *πp =π:* (5)

*imφ=*ð Þ *p* � *m* , (3)

*=*ð Þ *p* � *m*

*=*2*Z* � �, *Z* ¼ *z* � *i z*0,

2

*:* (4)

grating with a fractional topological charge (see the callout in **Figure 1**).

Ψ ¼ 2*NG eipπ* sin ð Þ *pπ* X∞

positive and negative topological charges:

*Nonlinear Optics - From Solitons to Similaritons*

*r* and *φ* are polar coordinates, *N* ¼ exp *iK*2

*mIm* j j *K* 2

*w*2 0*=*4 � �

ℓ*z* ¼ X∞

also [45])

**Figure 1.**

**166**

*of the OAM* ð Þ *Ω* ¼ *14 .*

*m*¼�∞

For a simplest model for our studies, we choose a scalar field of wave in the form of superposition *LGm n* of the lowest order with the *N* monochromatic beams, where *m* > 0 and *n* ¼ 0 are the azimuthal and radial indices, respectively, so that all the vortex modes in superposition have the same waist radius *w*0 at *z* ¼ 0. Let us consider the wave field in the waist plane *z* ¼ 0 in the form [25].

$$\Psi(r,\varphi,z=0) = \sum\_{m=0}^{N-1} \mathcal{C}\_{m} L G\_{n=0}^{m} = \sum\_{m=0}^{N-1} \mathcal{C}\_{m} r^{m} e^{j(m\varphi+\beta\_{n})} G(r,z)/N\_{m},\tag{6}$$

where *Nm* ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi 2�*m*�1 *m*!*π* p stands for the normalization factor, *Cm* is the mode amplitude, *βm* is the initial mode phase. We consider the nondegenerate case of the field representationð Þ *m* ≥0 or*m* ≤ 0 due to the axial symmetry; combined beams with different topological signs �*m* are indistinguishable in intensity moments.

Our goal is to analyze the distribution of the wave field intensity *ℑ*ð Þ *r*, *φ*, *z* ¼ 0 in such a way as to express the squared amplitudes and the initial phases of the vortex modes in terms of physically measured quantities in the region of the beam waist. For this, we make use of the intensity moments approach [35] in the form

$$J\_{p,q} = \iint\_{\mathcal{S}} \mathcal{M}\_{p,q}(r,\rho) \mathcal{G}(r,\rho) r \, dr \, d\rho / \iint\_{\mathcal{S}} \mathcal{G}(r,\rho) r \, dr \, d\rho,\tag{7}$$

where the intensity distribution is written as

$$\mathcal{B} = \Psi^\* \Psi = \sum\_{m=0}^{N-1} \frac{C\_m^2}{N\_m^2} r^{2m} G(r)^2 + 2 \sum\_{m, m'=0, \atop m, m'=0, \atop m \sim m'}^{N-1} \frac{C\_{m'} C\_m}{N\_m N\_{m'}} r^{m+m'} \cos \left[ (m - m') \wp \right] \cos \beta\_{m, m'} G(r)^2$$

$$-2 \sum\_{m, m'=0, \atop m, m'=0, \atop m, m'=0, \atop m, m'=0}^{N-1} \frac{C\_{m'} C\_m}{N\_m N\_{m'}} r^{m+m'} \sin \left[ (m - m') \wp \right] \sin \beta\_{m, m'} G(r)^2,$$

$$m > m' \tag{8}$$

*Mp*,*q*ð Þ *r*, *φ* stands for the intensity moments function, *p*, *q* ¼ 0, 1, 2, … . Now the problem is to choose a combination of intensity moments *Jp*,*q* in such a way as to exclude all terms of the last two sums and leave only the first one in Eq. (7). The equations for the squared amplitudes and the initial phases are separated if the moment function is written as *Mp*ð Þ¼ *r*, *φ rp*. Then the equations for the squared amplitudes take the form

$$\begin{aligned} \sum\_{m=0}^{N-1} \frac{(N-1+m)!}{m!} \mathcal{C}\_m^2 &= \sum\_{j=0}^p \binom{p}{j} J\_{2j, 2(p-j)}, p \le N-2, \\\\ \sum\_{m=0}^{N-1} \frac{m+1/2}{4\sqrt{2}m!} \Gamma(m+1/2) \mathcal{C}\_m^2 &= \sqrt{J\_{2,0} + J\_{0,2}}. \end{aligned} \tag{9}$$

second CCD2 camera (Thorlabs DCC1645) then the second-order intensity moment *Jxy* was measured and the average OAM per photon was calculated in accordance

*(a) Model of experimental setup for real-time vortex measurement and value of the OAM, P—polarizer, FF light filter of space, SLM—space light modulator, L1, L2—spherical lenses with focal length fsh, BS—beam splitter, CL—cylindrical lens with a focal length* f*cyl, CCD1,2—CCD-color cameras; (b and c) the vortex*

mð Þ m , *(c) β*mð Þ m *the OAM: theory* ℓz ¼ *0:39, experiment*

Before proceeding to the current measurements, it was necessary to adjust the experimental setup. For this, a number of calibration measurements were carried out. The digitization of the intensity distribution *ℑ*ð Þ *x*, *y* required normalization of the transverse coordinates *x*, *y* in units of the Gaussian beam waist at the plane of the CCD camera. To measure the initial beam waist radius *w*0, the second-order intensity moments *J*12 were used in accordance with the beam quality measurement method [46]. In order to ensure calibration measurements, a special computer program was developed that assigned random values of the amplitudes and initial phases of the vortex modes in Eq. (1) due to the random-number generator. Then, an appropriate diffraction grating was formed on the liquid crystal element of the SLM modulator and the calibration beam was restored. A typical intensity distribu-

tion of the calibration beam is shown in the callout of **Figure 2(c)**. The

corresponding spectra of the squared amplitudes j j *Сm* 2 and initial phases *βm* are presented in **Figure 2(b)** and **(c)**. Red and blue colors in **Figure 2(a)** and **(b)** present experimental and theoretical data, respectively. Note that the OAM measured in the second arm was ℓ*z* ¼ 4*:*2 and the first arm calculated according to Eq. (2) was ℓ*z* ¼ 4*:*5, while the theory gives ℓ*z* ¼ 4*:*9. The measurement error for *N* ¼ 10 beams did not exceed 3–4% for amplitudes and 5–6% for phases. The most interesting effect is observed in the region of OAM oscillations in **Figure 1(a)** with a characteristic perturbation of the holographic grating in **Figure 2(c)**. This effect is accompanied by a sharp restructuring of the vortex spectrum. For small perturbations of the singular beam, additional vortex modes arise only near the initial value of the topological charge *m* ¼ *M* (see [47]). However, in general case of the hologram perturbation, vortex modes with opposite signs of topological charges appear.

Such a vortex spectrum reconstruction cannot be detected using the above

We consider the case when the initial LG beam with a zero radial index *n* ¼ 0, *m* 6¼ 0 is subjected to local perturbation at the central region of the holographic grating. As a result of the perturbation, vortex modes with various types of topological charges appear in the beam. The method of the beam expansion in a series over the vortex modes depends on the type of the grating perturbation. We suppose

approach. It is required to expand the measurement technique.

**3.2 Degenerate case and vortex avalanche**

**169**

with the technique described in detail in the paper [26].

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

ℓ*z* ¼ *0:42; callout: the intensity distribution and corresponding diffraction grating.*

*spectrum of the beam array with N* ¼ *10: (b) С*2

**Figure 2.**

The right-hand sides of the linear Eq. (9) are easy to calculate, and the left-hand sides are directly measured experimentally using the Eq. (7) that gives the desired values of the squared amplitudes *С*2 *m*.

As intensity of array of the beams depends only on the difference of phase between two pairs of modes, in computation of the initial phases we will assume that one of beams phase is given (let us say *β*0 ¼ 0). Now it is important to select the combination of moments of intensity *Jp*,*q* that so either the second or the third sum in Eq. (2) is considered in the calculation. It is natural to calculate not all differences of phase *βm*,*m*0, but only *β*0,*m*. However, we did not find such combinations *Jp*,*q* that could filter out only terms with *β*0,*m*. Generally speaking, to calculate the phase difference *βm*,*m*0 it is required to find the *M* equations for variables *Xm*,*m*0 ¼ *CmCm*0 cos *βm m*0 and *Ym*,*m*0 ¼ *CmCm*0 sin *βm*,*m*0 in Eq. (2) the second and third sum, number of which is equal to the number of 2-combinations of *N* elements. It turns out that the number of equations can be significantly decreased if using the moments *J*2ð Þ *p*þ1 ,1 for *Xm*,*m*0 variables and *J*1,2ð Þ *p*þ1 for *Ym:m*0 variables. The system of linear equations for the phase difference are written in the form

$$J\_{2(p+1),1} = \frac{1}{2^{2(p+1)}} \sum\_{m=0}^{N-1} \sum\_{k=0}^{p} \binom{2p}{k} \left[ \frac{M\_{m,n\_1}}{N\_n N\_{n\_1}} X\_{m,n\_1} + \frac{M\_{m,n\_2}}{N\_m N\_{n\_2}} X\_{m,n\_2} \right],\tag{10}$$

$$J\_{1,2(p+1)} = \frac{1}{2^{2(p+1)}} \sum\_{m=0}^{N-1} \sum\_{k=0}^{p} (-1)^{p-k} \binom{2p}{k} \left[ \frac{M\_{m,n\_1}}{N\_m N\_{n\_1}} Y\_{m,n\_1} + \frac{M\_{m,n\_2}}{N\_m N\_{n\_2}} Y\_{m,n\_2} \right]. \tag{11}$$

where *n*1 ¼ j j *m* � 2ð Þ� *p* � *k* 3 , *n*2 ¼ j j *m* � 2ð Þ� *p* � *k* 1 , *Mm*,*n*1,2 ¼ Ð ∞ 0 *rm*þ*n*1,2þ1*G*2 *dr*. The number of linear equations in each system Eqs. (10) and (11) is *K* ¼ 3ð Þ *N* � 3 , *N* ≥6. It is noteworthy that Eqs. (10) and (11) contain only the terms with an odd difference of indices, including *β*0,*m*, so that a finite solution enables us to obtain all phases of the partial beams in the form tan *βm*,*m*0 ¼ *Ym*,*m*0*=Xm*,*m*0 .

A key element of the experimental setup in **Figure 2(a)** was a spatial light modulator SLM (Thorlabs EXULUS-4K1), which converted the fundamental *TEM*00 mode of the He-Ne laser (wavelength *λ* ¼ 0*:*6328 *mcm* and power 1 mW) into a composite beam with mode amplitudes *Cm* and initial phases *βm*. To make this possible, the laser beam was additionally filtered by the system FF. The beam splitter BS formed two working arms of the experimental setup. The beam in the first arm projected by a spherical lens onto the input pupil of the CCD1 camera (Thorlabs DCC1645) is subjected to the image computer processing. The result of the computer processing is a digitization of the beam intensity distribution and calculations of high-order intensity moments *Jpq*. Additional computer software made it possible to compose a system of linear Eqs. (9)–(11) and calculate both the mode amplitudes *Cm* and their initial phases *βm* in real time. The second arm was in use for the control measurement of the OAM of the composite beam. For this purpose, the beam was focused by a cylindrical lens CL onto the input pupil of the

*Digital Sorting of Optical Vortices in Perturbed Singular Beams DOI: http://dx.doi.org/10.5772/intechopen.91419*

#### **Figure 2.**

X *N*�1

*Nonlinear Optics - From Solitons to Similaritons*

ð Þ *N* � 1 þ *m* ! *m*! *C*2

*m* þ 1*=*2 4 ffiffi 2 p *m*!

*m* ¼ X *p*

Γð Þ *m* þ 1*=*2 *C*2

*m*.

linear equations for the phase difference are written in the form

2*p k*

ð Þ �1 *p*�*k* 2*p*

where *n*1 ¼ j j *m* � 2ð Þ� *p* � *k* 3 , *n*2 ¼ j j *m* � 2ð Þ� *p* � *k* 1 , *Mm*,*n*1,2 ¼ Ð ∞

obtain all phases of the partial beams in the form tan *βm*,*m*0 ¼ *Ym*,*m*0*=Xm*,*m*0 .

� � *Mm*,*n*1

*k*

*NnNn*1

� � *Mm*,*n*1

*K* ¼ 3ð Þ *N* � 3 , *N* ≥6. It is noteworthy that Eqs. (10) and (11) contain only the terms with an odd difference of indices, including *β*0,*m*, so that a finite solution enables us to

A key element of the experimental setup in **Figure 2(a)** was a spatial light modulator SLM (Thorlabs EXULUS-4K1), which converted the fundamental *TEM*00 mode of the He-Ne laser (wavelength *λ* ¼ 0*:*6328 *mcm* and power 1 mW) into a composite beam with mode amplitudes *Cm* and initial phases *βm*. To make this possible, the laser beam was additionally filtered by the system FF. The beam splitter BS formed two working arms of the experimental setup. The beam in the first arm projected by a spherical lens onto the input pupil of the CCD1 camera (Thorlabs DCC1645) is subjected to the image computer processing. The result of the computer processing is a digitization of the beam intensity distribution and calculations of high-order intensity moments *Jpq*. Additional computer software made it possible to compose a system of linear Eqs. (9)–(11) and calculate both the mode amplitudes *Cm* and their initial phases *βm* in real time. The second arm was in use for the control measurement of the OAM of the composite beam. For this purpose, the beam was focused by a cylindrical lens CL onto the input pupil of the

*NmNn*1

*dr*. The number of linear equations in each system Eqs. (10) and (11) is

*Xm*,*n*1 þ

� �

*Ym*,*n*1 þ

*Mm*,*n*2 *NmNn*2

� �

*Xm*,*n*2

*Mm*,*n*2 *NmNn*2

*Ym*,*n*2

, (10)

*:* (11)

X *p*

*k*¼0

*j*¼0

The right-hand sides of the linear Eq. (9) are easy to calculate, and the left-hand sides are directly measured experimentally using the Eq. (7) that gives the desired

As intensity of array of the beams depends only on the difference of phase between two pairs of modes, in computation of the initial phases we will assume that one of beams phase is given (let us say *β*0 ¼ 0). Now it is important to select the combination of moments of intensity *Jp*,*q* that so either the second or the third sum in Eq. (2) is considered in the calculation. It is natural to calculate not all differences of phase *βm*,*m*0, but only *β*0,*m*. However, we did not find such combinations *Jp*,*q* that could filter out only terms with *β*0,*m*. Generally speaking, to calculate the phase difference *βm*,*m*0 it is required to find the *M* equations for variables *Xm*,*m*0 ¼ *CmCm*0 cos *βm m*0 and *Ym*,*m*0 ¼ *CmCm*0 sin *βm*,*m*0 in Eq. (2) the second and third sum, number of which is equal to the number of 2-combinations of *N* elements. It turns out that the number of equations can be significantly decreased if using the moments *J*2ð Þ *p*þ1 ,1 for *Xm*,*m*0 variables and *J*1,2ð Þ *p*þ1 for *Ym:m*0 variables. The system of

*p j* � �

*m* ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *J*2,0 þ *J*0,2 p *:*

*J*2*j*,2ð Þ *p*�*j* , *p*≤ *N* � 2,

(9)

*m*¼0

*N* X�1 *m*¼0

values of the squared amplitudes *С*2

*J*2ð Þ *p*þ1 ,1 ¼ 1

22ð Þ *p*þ1

*J*1,2ð Þ *p*þ1 ¼ 1

0 *rm*þ*n*1,2þ1*G*2

**168**

22ð Þ *p*þ1

X *N*�1

*m*¼0,

X *N*�1

*m*¼0,

X *p*

*k*¼0

*(a) Model of experimental setup for real-time vortex measurement and value of the OAM, P—polarizer, FF light filter of space, SLM—space light modulator, L1, L2—spherical lenses with focal length fsh, BS—beam splitter, CL—cylindrical lens with a focal length* f*cyl, CCD1,2—CCD-color cameras; (b and c) the vortex spectrum of the beam array with N* ¼ *10: (b) С*2 mð Þ m , *(c) β*mð Þ m *the OAM: theory* ℓz ¼ *0:39, experiment* ℓ*z* ¼ *0:42; callout: the intensity distribution and corresponding diffraction grating.*

second CCD2 camera (Thorlabs DCC1645) then the second-order intensity moment *Jxy* was measured and the average OAM per photon was calculated in accordance with the technique described in detail in the paper [26].

Before proceeding to the current measurements, it was necessary to adjust the experimental setup. For this, a number of calibration measurements were carried out. The digitization of the intensity distribution *ℑ*ð Þ *x*, *y* required normalization of the transverse coordinates *x*, *y* in units of the Gaussian beam waist at the plane of the CCD camera. To measure the initial beam waist radius *w*0, the second-order intensity moments *J*12 were used in accordance with the beam quality measurement method [46]. In order to ensure calibration measurements, a special computer program was developed that assigned random values of the amplitudes and initial phases of the vortex modes in Eq. (1) due to the random-number generator. Then, an appropriate diffraction grating was formed on the liquid crystal element of the SLM modulator and the calibration beam was restored. A typical intensity distribution of the calibration beam is shown in the callout of **Figure 2(c)**. The corresponding spectra of the squared amplitudes j j *Сm* 2 and initial phases *βm* are presented in **Figure 2(b)** and **(c)**. Red and blue colors in **Figure 2(a)** and **(b)** present experimental and theoretical data, respectively. Note that the OAM measured in the second arm was ℓ*z* ¼ 4*:*2 and the first arm calculated according to Eq. (2) was ℓ*z* ¼ 4*:*5, while the theory gives ℓ*z* ¼ 4*:*9. The measurement error for *N* ¼ 10 beams did not exceed 3–4% for amplitudes and 5–6% for phases. The most interesting effect is observed in the region of OAM oscillations in **Figure 1(a)** with a characteristic perturbation of the holographic grating in **Figure 2(c)**. This effect is accompanied by a sharp restructuring of the vortex spectrum. For small perturbations of the singular beam, additional vortex modes arise only near the initial value of the topological charge *m* ¼ *M* (see [47]). However, in general case of the hologram perturbation, vortex modes with opposite signs of topological charges appear. Such a vortex spectrum reconstruction cannot be detected using the above approach. It is required to expand the measurement technique.

#### **3.2 Degenerate case and vortex avalanche**

We consider the case when the initial LG beam with a zero radial index *n* ¼ 0, *m* 6¼ 0 is subjected to local perturbation at the central region of the holographic grating. As a result of the perturbation, vortex modes with various types of topological charges appear in the beam. The method of the beam expansion in a series over the vortex modes depends on the type of the grating perturbation. We suppose (and then prove) that the chosen perturbation does not excite LG modes with radial indices *n* ¼ 0, but only modes with different topological charges *m* appear [39]. Then we write the complex amplitude of the perturbed beam in the form

$$\Psi(R,\rho,z) = \sum\_{m=-N}^{N} \frac{C\_m R^{|m|} e^{im\ \rho} G(R)}{M\_m \sigma^{|m|+1}(z)},\tag{12}$$

indices j j *m* � *n* ¼ 2 in Eq. (14) are imaginary, therefore the cross terms of the beam

and write a system of linear equations, the number of which is equal to the number of variables, as was done above for a non-degenerate vortex array. However, as shown by the assessed computations, as the indices *p* and *q* increase along with the 2*N*

of Eq. (14) cannot be closed. However, one can act otherwise. First, we make two

functions in the form of Hermite polynomials Eq. (14) as *Mp*,*0*ð Þ¼ *u Hp u=* ffiffi

It seems that for solving this problem it is reasonable to choose intensity moment

2 � � p . Next, one can use the orthogonality condition of the HG beams

*m*, the 2*N* � 6 cross terms *C*�*nС*�*m* and *C*�*nС*∓*m* are also added. The system

∞ð

ð ∞

cos *η*

2*n*þ1*=*2 *Nn*,*n*

*ξ*2*mL*2*m n*

*v* ffiffi 2 p

*ξ*2*nLn*

*ξ*2 2 � �*e* �*ξ*2

�*n* will be found, if we will make use the

*p*ð Þ *x* is Laguerre polynomials. Then,

*C*�*nC*∓*m*

*p*,*q*. Each equation contains 2*N*

� �*ℑ*ð Þ*η dη=J*00, (17)

� �*ℑ*ð Þ *x:y dx dy:*

*ξ*2 2 � �*e* �*ξ*2 4

4 , (16)

�*n* with the similar

*n* with only differ-

*=J*00, (18)

�∞

�∞

*N*

*n*¼0 *C*2 �*n*

ð Þ �1 *n*þ2*m*2*n*þ1*=*2 *Nm*,*n*

Thus, after the Fourier transform the distribution of intensity *ℑ*ð Þ *ξ*, *η* contains all the squared amplitudes, and, also, the cross terms with the difference j j *m* � j j *n* ¼ �4. The cross terms *ℑ*�*т*,*m* and *ℑm*,�*n* refer to vortices with positive and negative

ent coefficients. The expressions are received for the terms, that *ℑ*�*т*,�ð Þ *n*þ2*m* ð Þ*η* and

full basis of the Laguerre-Gaussian modes for expansion of the intensity moments:

ð Þ � *p*,2*q* ¼

X *N*

!

*m*, *n*¼1

The system of Eq. (18) enables us to find the values of the squared amplitudes

ð ∞

*ηL*2*q p*

ð Þ 2*m* þ *n* !ð Þ 1*=*2 *n n*!4�ð Þ 2*m*þ*n*þ1 *Nnm*

�

*m* and *N* � 3 cross terms *CmCn*. Increasing the value of the

*η*2 2

�∞

2 � � p and

(15)

intensity distribution with such indices difference disappear.

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

Fourier transforms of the intensity distribution *ℑ*ð Þ *x*, *y* in the form

� �*ℑ*ð Þ *x:y dx dy*, *ℑ*ð Þ¼ *η*

In particular, the terms *ℑ*�*m:*�*n* in the intensity distribution *ℑ* are

� �*ℑ*ð Þ *x:y dx dy* ¼ X

*M0*,*q*ð Þ¼ *v Hq v=* ffiffi

variables *C*2

*ℑ*ð Þ¼ *ξ*

*ℑ*�*m:*�*n*ð Þ¼ *ξ*

*J* ð Þ þ *p*,2*q* ¼

*J* ∓

**171**

*p*,*q* ¼ X *N*

*n*¼1 *C*2

terms for amplitudes *C*2

ð ∞

*ξL*2*q p*

�∞

∞ð

ð ∞

cos *ξ u*

ffiffi 2 p

cos *ξ u*

ffiffi 2 p

*C*�*nС*�ð Þ *n*þ2*m*

charges. Besides, the terms *ℑ*�*т*,�ð Þ *n*þ2*m* ð Þ*ξ* includes amplitudes *С*2

coefficients, and the member *ℑт*,ð Þ *n*þ2*m* ð Þ*ξ* includes amplitudes *С*2

� �*ℑ*ð Þ*ξ dξ=J*00, *J*

�∞

∞ð

∞ð

�∞

þ X *N*

> *m*, *n*¼0, j j *m*�*n* 6¼0, 2

*ℑт*,ð Þ *n*þ2*m* ð Þ*η* , if we swap the signs of the indices. The equations for squared amplitudes *С*2

> *ξ*2 2

where *J*00 is a total beam intensity and *Lq*

using [39], we find a system of linear equations

∓*n* þ ð Þ 2*q* þ *p* !ð Þ 1*=*2 *p p*!4�ð Þ 2*q*þ*p*þ1

via the measured values of the intensity moments *J*

�∞

�∞

where *Mm* is a normalized factor and *σ*ð Þ¼ *z* 1 � *i z=z*0.

It is important to notice the intensity moments *Jp:q* are degenerate with respect to the sign of the vortex topological charge *m* in the case of axial symmetry of the modes in Eq. (12). To remove degeneracy of the intensity moments, it is important to change symmetry of singular beams, but without destroying of structure beams so that the vortex modes with opposite topological charges in sign have distinct geometric contours. Such a requirement provides by the astigmatic transformation with a cylindrical lens of the Laguerre-Gauss beam, which was considered detail in [34] (the general theory of paraxial astigmatic transformations can be found in [41]). We suppose that the cylindrical lens with a focal length *f* is located at the plane *z* ¼ 0. A paraxial beam with complex amplitude (12) is projected at the lens input. The plane of the beam waist is matched with the plane of the lens. Besides, we restrict our study to zero initial phases of mode beams, i.e., amplitudes *Cm* are real values. Then, following [34], we write the complex amplitude of a combined beam at the wave diffraction zone as

$$\Psi\_{\rm CL}(x, y, z) = \sum\_{m = -N}^{N} A\_m H\_m(F) \exp\left(\Phi\right), \tag{13}$$

$$\begin{split} \text{where } A\_{m} &= -i\frac{x\_{0}}{x} \frac{1}{q q\_{0}} \left( -\frac{w\_{0}}{q\_{0}} \right)^{m} \left( \mathbf{1} - \frac{q\_{0}^{2}}{q^{2}} \right), \ \Phi = i\frac{k}{2x} \left( \mathbf{x}^{2} + \mathbf{y}^{2} \right) - \left( \frac{x\_{0}}{x w\_{0} q} \right)^{2} \mathbf{x}^{2} - \mathbf{1} \\ \left( \frac{x\_{0}}{x w\_{0} q\_{0}} \right)^{2} \mathbf{y}^{2}, \ F = \frac{x\_{0}}{w\_{0} x} \frac{\frac{q\_{0}}{q} \mathbf{x} - \frac{q}{q\_{0}} \mathbf{y}}{\sqrt{q^{2} - q\_{0}^{2}}}, \ q\_{0}^{2} = \mathbf{1} - i \frac{x\_{0}}{x}, q^{2} = \mathbf{1} + i \frac{x\_{0}}{x\_{1}}, z\_{1} = \frac{x f}{x - f} \dots \end{split}$$

Let us consider new normalized coordinates so that the beam axes are directed by an angle *π=*4 relative to the axes of the cylindrical lens: *u* ¼ ð Þ *x* þ *y =w*0, *v* ¼ ð Þ *x* � *y =w*0, we will consider field of paraxial beam in double focus plane *z* ¼ 2*f*, and also require that *z*0*=*2*f* ¼ 1. Then distribution of the beam intensity is written as

$$\mathcal{B} = |\Psi\_{CL}|^2 = \left\{ \sum\_{m,n=0}^{N} C\_{-n} C\_{-m} \frac{A\_{m,n}}{N\_{m,n}} H\_n \left( \frac{u}{\sqrt{2}} \right) H\_m \left( \frac{u}{\sqrt{2}} \right) + \sum\_{m,n=1}^{N} C\_n C\_m \frac{A\_{m,n}}{N\_{m,n}} H\_n \left( \frac{v}{\sqrt{2}} \right) H\_m \left( \frac{v}{\sqrt{2}} \right) \right. \tag{14}$$

$$+ \sum\_{n=0,\atop m=1}^{N} C\_{-n} C\_m \frac{A\_{m,n}}{N\_{m,n}} H\_n \left( \frac{u}{\sqrt{2}} \right) H\_m \left( \frac{v}{\sqrt{2}} \right) + \sum\_{n=1,\atop m=0}^{N} C\_n C\_{-m} \frac{A\_{m,n}}{N\_{m,n}} H\_n \left( \frac{v}{\sqrt{2}} \right) H\_m \left( \frac{u}{\sqrt{2}} \right) \right\} \exp\left( -\frac{u^2 + v^2}{2} \right), \tag{15}$$

where *Hn*ð Þ *x* is Hermite polynomials, *Am*,*n* ¼ �ð Þ1 j jþ*n* j j *m ei* j j� *m* j j *n* 4 *π*, *N*2 *n*,*m* ¼ *π* 4 *wo* 2 � �*n*þ*m*2*n*þ*m*þ1 *n*!*m*!. As can be seen from Eq. (14), the terms corresponding to optical vortices with positive and negative topological charges are only partially divided. The last two terms characterize the cross terms *C*�*nСm* and *CnС*�*m*. Besides, the first two terms also contain cross amplitudes *C*�*nС*�*m* and *CnСm*ð Þ j j *n* 6¼ j j *m* . On the other hand, our task is to measure only the squared amplitudes *C*2 *m*. It is important to note, that the factors *Am*,*n* ¼ �*An*,*m* ¼ *i* with the difference of the

### *Digital Sorting of Optical Vortices in Perturbed Singular Beams DOI: http://dx.doi.org/10.5772/intechopen.91419*

(and then prove) that the chosen perturbation does not excite LG modes with radial indices *n* ¼ 0, but only modes with different topological charges *m* appear [39]. Then we write the complex amplitude of the perturbed beam in the form

*CmR*j j *m ei m φG R*ð Þ

*Mmσ*j jþ *m* 1ð Þ*z* , (12)

*AmHm*ð Þ *F* exp ð Þ Φ , (13)

*z*�*f :*.

*Am*,*n Nm*,*n Hn v* ffiffi 2 p � � *Hm v* ffiffi 2 p

j j� *m* j j *n* 4 *π*, *N*2

9 >>>>=

>>>>;

*n*,*m* ¼

exp � *u*2 þ *v*2 2 � �

*m*. It is

(14)

*zw*0*q* � �2

*x*2 �

2*z x*2 þ *y*2 ð Þ� *z*0

*z*1 , *z*1 ¼ *z f*

*CnСm*

*N*

*m*¼�*N*

It is important to notice the intensity moments *Jp:q* are degenerate with respect to the sign of the vortex topological charge *m* in the case of axial symmetry of the modes in Eq. (12). To remove degeneracy of the intensity moments, it is important to change symmetry of singular beams, but without destroying of structure beams so that the vortex modes with opposite topological charges in sign have distinct geometric contours. Such a requirement provides by the astigmatic transformation with a cylindrical lens of the Laguerre-Gauss beam, which was considered detail in [34] (the general theory of paraxial astigmatic transformations can be found in [41]). We suppose that the cylindrical lens with a focal length *f* is located at the plane *z* ¼ 0. A paraxial beam with complex amplitude (12) is projected at the lens input. The plane of the beam waist is matched with the plane of the lens. Besides, we restrict our study to zero initial phases of mode beams, i.e., amplitudes *Cm* are real values. Then, following [34], we write the complex amplitude of a combined

Ψð Þ¼ *R*, *φ*, *z* X

where *Mm* is a normalized factor and *σ*ð Þ¼ *z* 1 � *i z=z*0.

*Nonlinear Optics - From Solitons to Similaritons*

Ψ*CL*ð Þ¼ *x*, *y*, *z* X

� *w*0 *q*0 � �*m*

> *Am*,*n Nm*,*n Hn u* ffiffi 2 p � � *Hm u* ffiffi 2 p � �

*N*

*m*¼�*N*

by an angle *π=*4 relative to the axes of the cylindrical lens: *u* ¼ ð Þ *x* þ *y =w*0, *v* ¼ ð Þ *x* � *y =w*0, we will consider field of paraxial beam in double focus plane *z* ¼ 2*f*, and also require that *z*0*=*2*f* ¼ 1. Then distribution of the beam intensity is written as

> þ X*N n*¼1, *m*¼0

where *Hn*ð Þ *x* is Hermite polynomials, *Am*,*n* ¼ �ð Þ1 j jþ*n* j j *m ei*

, Φ ¼ *i k*

*z* , *q*2 ¼ 1 þ *i z*0

þ X*N m*, *n*¼1

*Am*,*n Nm*,*n Hn v* ffiffi 2 p � � *Hm u* ffiffi 2 p � �

*n*!*m*!. As can be seen from Eq. (14), the terms corresponding to

( � �

*CnС*�*m*

optical vortices with positive and negative topological charges are only partially divided. The last two terms characterize the cross terms *C*�*nСm* and *CnС*�*m*. Besides, the first two terms also contain cross amplitudes *C*�*nС*�*m* and *CnСm*ð Þ j j *n* 6¼ j j *m* . On the other hand, our task is to measure only the squared amplitudes *C*2

important to note, that the factors *Am*,*n* ¼ �*An*,*m* ¼ *i* with the difference of the

Let us consider new normalized coordinates so that the beam axes are directed

1 � *q*2 0 *q*2 � �

0 ¼ 1 � *i z*0

beam at the wave diffraction zone as

*z* 1 *qq*0

where *Am* ¼ �*i z*0

*ℑ* ¼ j j Ψ*CL* 2 ¼ X*N*

*C*�*nСm*

� �*n*þ*m*2*n*þ*m*þ1

*y*2, *F* ¼ *z*0 *w*0*z i q*0 *q x*� *q q*0 *y* ffiffiffiffiffiffiffiffiffi *q*2�*q*2 0 p , *q*2

*m*, *n*¼0

*Am*,*n Nm*,*n Hn u* ffiffi 2 p � � *Hm v* ffiffi 2 p � �

*C*�*nС*�*m*

*z*0 *zw*0*q*0 � �2

þ X*N n*¼0, *m*¼1

*π* 4 *wo* 2

**170**

indices j j *m* � *n* ¼ 2 in Eq. (14) are imaginary, therefore the cross terms of the beam intensity distribution with such indices difference disappear.

It seems that for solving this problem it is reasonable to choose intensity moment functions in the form of Hermite polynomials Eq. (14) as *Mp*,*0*ð Þ¼ *u Hp u=* ffiffi 2 � � p and *M0*,*q*ð Þ¼ *v Hq v=* ffiffi 2 � � p . Next, one can use the orthogonality condition of the HG beams and write a system of linear equations, the number of which is equal to the number of variables, as was done above for a non-degenerate vortex array. However, as shown by the assessed computations, as the indices *p* and *q* increase along with the 2*N* variables *C*2 *m*, the 2*N* � 6 cross terms *C*�*nС*�*m* and *C*�*nС*∓*m* are also added. The system of Eq. (14) cannot be closed. However, one can act otherwise. First, we make two Fourier transforms of the intensity distribution *ℑ*ð Þ *x*, *y* in the form

$$\mathfrak{G}(\xi) = \int\_{-\infty}^{\infty} \int\_{-\infty}^{\infty} \cos \left( \xi \frac{u}{\sqrt{2}} \right) \mathfrak{G}(\varkappa; \mathfrak{y}) \, d\varkappa \, d\mathfrak{y}, \ \mathfrak{J}(\eta) = \int\_{-\infty}^{\infty} \int\_{-\infty}^{\infty} \cos \left( \eta \frac{v}{\sqrt{2}} \right) \mathfrak{G}(\varkappa; \mathfrak{y}) \, d\varkappa \, d\eta. \tag{15}$$

In particular, the terms *ℑ*�*m:*�*n* in the intensity distribution *ℑ* are

$$\begin{split} \mathfrak{I}\_{-m,-n}(\boldsymbol{\xi}) &= \int\_{-\infty}^{\infty} \int\_{-\infty}^{\infty} \cos\left(\boldsymbol{\xi} \, \frac{u}{\sqrt{2}}\right) \mathfrak{I}(\boldsymbol{x}, \boldsymbol{y}) \, d\boldsymbol{x} \, d\boldsymbol{y} = \sum\_{n=0}^{N} C\_{-n}^{2} \frac{2^{n+1/2}}{N\_{n,n}} \xi^{2n} L\_{n} \left(\frac{\xi^{2}}{2}\right) e^{-\frac{\xi^{2}}{4}} \\ &+ \sum\_{m,n=0,1}^{N} C\_{-n} C\_{-(n+2m)} \frac{(-1)^{n+2m} 2^{n+1/2}}{N\_{m,n}} \xi^{2m} L\_{n}^{2m} \left(\frac{\xi^{2}}{2}\right) e^{-\frac{\xi^{2}}{4}}, \end{split} \tag{16}$$

Thus, after the Fourier transform the distribution of intensity *ℑ*ð Þ *ξ*, *η* contains all the squared amplitudes, and, also, the cross terms with the difference j j *m* � j j *n* ¼ �4. The cross terms *ℑ*�*т*,*m* and *ℑm*,�*n* refer to vortices with positive and negative charges. Besides, the terms *ℑ*�*т*,�ð Þ *n*þ2*m* ð Þ*ξ* includes amplitudes *С*2 �*n* with the similar coefficients, and the member *ℑт*,ð Þ *n*þ2*m* ð Þ*ξ* includes amplitudes *С*2 *n* with only different coefficients. The expressions are received for the terms, that *ℑ*�*т*,�ð Þ *n*þ2*m* ð Þ*η* and *ℑт*,ð Þ *n*þ2*m* ð Þ*η* , if we swap the signs of the indices.

The equations for squared amplitudes *С*2 �*n* will be found, if we will make use the full basis of the Laguerre-Gaussian modes for expansion of the intensity moments:

$$J\_{p,2q}^{(+)} = \int\_{-\infty}^{\infty} \xi L\_p^{2q} \left(\frac{\xi^2}{2}\right) \Im(\xi) \, d\xi / J\_{00}, \ J\_{p,2q}^{(-)} = \int\_{-\infty}^{\infty} \eta L\_p^{2q} \left(\frac{\eta^2}{2}\right) \Im(\eta) \, d\eta / J\_{00},\tag{17}$$

where *J*00 is a total beam intensity and *Lq p*ð Þ *x* is Laguerre polynomials. Then, using [39], we find a system of linear equations

$$J\_{p,q}^{\mp} = \left(\sum\_{n=1}^{N} \mathcal{C}\_{\mp n}^{2} + \frac{(2q+p)!(1/2)^{p}}{p!4^{-(2q+p+1)}} \sum\_{m,n=1}^{N} \frac{(2m+n)!(1/2)^{n}}{n!4^{-(2m+n+1)}N\_{nm}} \mathcal{C}\_{\pm n} \mathcal{C}\_{\mp m}\right) / f\_{00}, \tag{18}$$

The system of Eq. (18) enables us to find the values of the squared amplitudes via the measured values of the intensity moments *J* � *p*,*q*. Each equation contains 2*N* terms for amplitudes *C*2 *m* and *N* � 3 cross terms *CmCn*. Increasing the value of the

indices *p* and *q* does not change the number of variables, so the system of Eq. (18) can be closed. Note that each of the Eq. (18) can be solved independently of each other and find the amplitudes *C*2 *m* and *C*2 �*m*. The total number of linear equations is relative to signs ð Þ � , so that the systems can be easy solved for finding the squared amplitudes in real time using modern computer software.

For measurements, the experimental setup shown in **Figure 2(a)** was used. However, in contrast to the non-degenerate case, the main experimental data were measured in the second arm, provided that the detection plane was located at the double focus region of the cylindrical lens, while the first arm was used for general adjustment [39].

We will implement the above technique for a new type of combined beams with conspicuous dips and bursts in the OAM spectrum similar to that shown in **Figure 1(a)** and **(b)**. Let us write the complex amplitude of the combined vortex beam in the form

$$\Psi(r,\rho,z) = \sin\left(p\pi\right)e^{i\,p\pi} \sum\_{m=-N}^{N} \binom{\cos\left(m\pi/2\right)}{\sin\left(m\pi/2\right)} \frac{\Omega^{|m|}}{|m|!(M-m+\delta p)} \frac{R^{|m|}e^{im\,\rho}}{\sigma(z)^{|m|+1}} \exp\left(-R^2\right),\tag{19}$$

where Ω is a scale parameter, *δp* stands for the perturbation of the holographic grating responsible for shaping the beam, *M* is a topological charge of the nonperturbed forked grating. Factors cosð Þ *mπ=*2 or sin ð Þ *mπ=*2 thin out the vortex spectrum remaining the modes with either even or odd *m*� indices. This makes the dips in the OAM spectrum more distinct. The complex amplitude (12) corresponds to the OAM in the form

$$\mathcal{E}\_x = \sum\_{m=-N}^{N} m \binom{\cos^2(m\pi/2)}{\sin^2(m\pi/2)} \frac{\mathfrak{L}^{2|m|}}{|m|^2 \left(M - m + \delta p\right)^2} / \sum\_{m=-N}^{N} \binom{\cos^2(m\pi/2)}{\sin^2(m\pi/2)} \frac{\mathfrak{L}^{2|m|}}{|m|^2 \left(M - m + \delta p\right)^2}. \tag{20}$$

**Figure 1(b)** and **(c)** displays the case of even vortices beams in the curves of ℓ*z*ð Þ *p* . The spectrum of OAM in **Figure 1(b)** characterizes the general view of resonant bursts and dips. The Ω parameter of scale admits displace along the axis of the resonant regions *p* ¼ *M* þ *δp*. **Figure 1(c)** (curves 1 and 2) describes their characteristic features. In the even vortex modes of integer topological charge are located resonant dips. Swaps of the cosine with the sine in Eq. (19) displaces along the axis *p* of the resonant regions by one. Increase of the topological charge *M* is accompanied by resonances compression, until their depth gets become approximately equal to *M* at Ω ¼ 18 and *p*< 6. For large parameters Ω > >20, the OAM is doubled ℓ*z* ¼ 2*M* in range of the initial topological charges *M*. As was shown in [47, 48], a perturbation *δp* causes a local distortion of the forked grating defect for the large values parameter *p* shown in **Figure 3(a)**.

charge from integer values *M* immediately causes a resonance dip in the OAM spectrum (**Figure 5(a)**). With increasing *M* the depth of the dip decreases while its

*18*, *M* ¼ *4; callouts: intensity distributions at the beam waist after spherical (sph) with focal length f* ¼ *0:1*m

*Forked defect perturbations of the holographic gratings with (a and b) M* ¼ *20, Ω* ¼ *2:1, (c and d) M* ¼ *4, Ω* ¼ *14: (а) δp* ¼ *10*�*3, (b) δp* ¼ *0:1, (c) δp* ¼ �*10*�*3 (d) δp* ¼ �*0:5, callouts: intensity distributions of the*

Absolutely other situation arises when the holographic grating with *M* >10 and a forked defect is subjected to a small perturbation. As shown in **Figure 3(c)**, even a weak perturbation *δp* ¼ �0*:*001 changes drastically the entire relief of the holographic grating although the intensity distribution. In the result, even small differences in the grating relief affects inevitably the shape of the vortex spectrum

*Dips and bursts in the OAM spectrum ℓz*ð Þ *p at various beam parameters Ω and M*, *p* ¼ *M* þ *δp: (a and b)*

*m*ð Þ *m for perturbed singular beams with (a) M* ¼ *20, Ω* ¼ *3:5 and (b) Ω* ¼

shape is distorted (see **Figure 5(b)**).

*and cylindrical (cyl) lens with f cyl* ¼ *0:5 m; red color—an initial mode.*

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

*M* ¼ *20, (c and d) M* ¼ *4; solid lines—theory, circles—experiment.*

**Figure 3.**

**Figure 4.**

**Figure 5.**

**173**

*Experimental vortex spectra C2*

*beams restored by the perturbed holograms.*

As the perturbation grows the area of the hologram defect increases (**Figure 3(b)**). Small variations in the holographic grating structure lead to a cardinal reconstruction of the vortex spectrum (**Figure 4(а)**). Weak perturbations *δp* � 10�2 of the grating cause an avalanche of optical vortices. The satellites appear in the OAM spectrum, the maxima of which has topological charges less than the initial *M*. The energy of the basic mode *M* flows partially into the spectral satellites, that are generated both in the region of positive *m* >0 and negative *m* >0 topological charges.

For large perturbations *δp* � 0*:*1, a vortex avalanche picks almost all the energy out of the basic mode. The vortex avalanche at small deviations *δp* of the topological

**Figure 3.**

indices *p* and *q* does not change the number of variables, so the system of Eq. (18) can be closed. Note that each of the Eq. (18) can be solved independently of each

relative to signs ð Þ � , so that the systems can be easy solved for finding the squared

For measurements, the experimental setup shown in **Figure 2(a)** was used. However, in contrast to the non-degenerate case, the main experimental data were measured in the second arm, provided that the detection plane was located at the double focus region of the cylindrical lens, while the first arm was used for general

We will implement the above technique for a new type of combined beams with conspicuous dips and bursts in the OAM spectrum similar to that shown in **Figure 1(a)** and **(b)**. Let us write the complex amplitude of the combined vortex

� � Ωj j *m*

where Ω is a scale parameter, *δp* stands for the perturbation of the holographic grating responsible for shaping the beam, *M* is a topological charge of the nonperturbed forked grating. Factors cosð Þ *mπ=*2 or sin ð Þ *mπ=*2 thin out the vortex spectrum remaining the modes with either even or odd *m*� indices. This makes the dips in the OAM spectrum more distinct. The complex amplitude (12) corresponds

2 ð Þ *M* � *m* þ *δp* 2 *=* X*N*

**Figure 1(b)** and **(c)** displays the case of even vortices beams in the curves of ℓ*z*ð Þ *p* . The spectrum of OAM in **Figure 1(b)** characterizes the general view of resonant bursts and dips. The Ω parameter of scale admits displace along the axis of the resonant regions *p* ¼ *M* þ *δp*. **Figure 1(c)** (curves 1 and 2) describes their characteristic features. In the even vortex modes of integer topological charge are located resonant dips. Swaps of the cosine with the sine in Eq. (19) displaces along the axis *p* of the resonant regions by one. Increase of the topological charge *M* is accompanied by resonances compression, until their depth gets become approximately equal to *M* at Ω ¼ 18 and *p*< 6. For large parameters Ω > >20, the OAM is doubled ℓ*z* ¼ 2*M* in range of the initial topological charges *M*. As was shown in [47, 48], a perturbation *δp* causes a local distortion of the forked grating defect for

As the perturbation grows the area of the hologram defect increases (**Figure 3(b)**). Small variations in the holographic grating structure lead to a cardinal reconstruction of the vortex spectrum (**Figure 4(а)**). Weak perturbations *δp* � 10�2 of the grating cause an avalanche of optical vortices. The satellites appear in the OAM spectrum, the maxima of which has topological charges less than the initial *M*. The energy of the basic mode *M* flows partially into the spectral satellites, that are generated both in the

For large perturbations *δp* � 0*:*1, a vortex avalanche picks almost all the energy out of the basic mode. The vortex avalanche at small deviations *δp* of the topological

*m*¼�*N*

j j *m* !ð Þ *M* � *m* þ *δp*

cos 2ð Þ *mπ=*2 sin 2

ð Þ *mπ=*2

� � Ω2j j *m*

j j *m* !

cosð Þ *mπ=*2 sin ð Þ *mπ=*2

�*m*. The total number of linear equations is

*R*j j *m ei m φ*

*σ*ð Þ*z* j jþ *m* 1 exp �*R*2 � �,

(19)

2 ð Þ *M* � *m* þ *δp* 2 *:*

(20)

*m* and *C*2

amplitudes in real time using modern computer software.

*i Pπ* X *N*

*m*¼�*N*

other and find the amplitudes *C*2

*Nonlinear Optics - From Solitons to Similaritons*

adjustment [39].

beam in the form

Ψð Þ¼ *r*, *φ*, *z* sin ð Þ *pπ e*

to the OAM in the form

*m* cos 2ð Þ *mπ=*2 sin 2

ð Þ *mπ=*2

� � Ω2j j *m*

the large values parameter *p* shown in **Figure 3(a)**.

region of positive *m* >0 and negative *m* >0 topological charges.

j j *m* !

ℓ*z* ¼ X*N m*¼�*N*

**172**

*Forked defect perturbations of the holographic gratings with (a and b) M* ¼ *20, Ω* ¼ *2:1, (c and d) M* ¼ *4, Ω* ¼ *14: (а) δp* ¼ *10*�*3, (b) δp* ¼ *0:1, (c) δp* ¼ �*10*�*3 (d) δp* ¼ �*0:5, callouts: intensity distributions of the beams restored by the perturbed holograms.*

#### **Figure 4.**

*Experimental vortex spectra C2 m*ð Þ *m for perturbed singular beams with (a) M* ¼ *20, Ω* ¼ *3:5 and (b) Ω* ¼ *18*, *M* ¼ *4; callouts: intensity distributions at the beam waist after spherical (sph) with focal length f* ¼ *0:1*m *and cylindrical (cyl) lens with f cyl* ¼ *0:5 m; red color—an initial mode.*

charge from integer values *M* immediately causes a resonance dip in the OAM spectrum (**Figure 5(a)**). With increasing *M* the depth of the dip decreases while its shape is distorted (see **Figure 5(b)**).

Absolutely other situation arises when the holographic grating with *M* >10 and a forked defect is subjected to a small perturbation. As shown in **Figure 3(c)**, even a weak perturbation *δp* ¼ �0*:*001 changes drastically the entire relief of the holographic grating although the intensity distribution. In the result, even small differences in the grating relief affects inevitably the shape of the vortex spectrum

#### **Figure 5.**

*Dips and bursts in the OAM spectrum ℓz*ð Þ *p at various beam parameters Ω and M*, *p* ¼ *M* þ *δp: (a and b) M* ¼ *20, (c and d) M* ¼ *4; solid lines—theory, circles—experiment.*

j j *Cm*ð Þ *m* 2 and Ω >10 induce the vortex avalanche (see **Figure 4(b)**). However, the avalanche type differs significantly from that in the previous case. Indeed, as well as in the above case, the perturbation provokes appearing resonance satellites with both positive and negative topological charges (TC). However, their values are always greater than those of the initial TC. In this case, energy is also pumped from the basic mode to the satellites but with higher topological charges under the condition of the resonance *m* ¼ *p*. In the OAM spectrum, it looks like the OAM burst at the resonance as shown in **Figure 5(c)** and **(d)** but the maxima ℓ*z* are at the region of fractional values *p*. As the scale parameter increases, the height of the flash decreases, and the resonance circuit broadens (**Figure 5(d)**).

**Figure 6(b)** and **(c)** shows a computer simulation (**Figure 6(b)**) and experiment (**Figure 6(c)**) of the intensity distribution *ℑ*ð Þ¼ *r*, *φ*, *z* ΨΨ∗ along the vortex beam length with *m* ¼ 15 modulated by the sector aperture with *α* ¼ *π=*4. It is important to note, that the cut out dark sector from the beam intensity distribution does not lead to its self-healing. When the beam propagates, the dark sector rotates synchronously. We showed, both theoretically and experimentally, that beam self-healing

We consider vortex beams spectra for topological charge that most clearly reflect properties of the sector perturbation. Computer simulation and experimental data of typical vortex spectra is shown in **Figure 7**. We revealed a clear maximum of

45°. The short part of the beam intensity redistributed symmetric at the neighboring vortices orders. But, if we increase the angle *α*, we see a violation of the symmetric distribution of intensity among the vortex modes. There formed a second maximum in the spectrum range of negative topological charges *n* <0. It is important to note, that the authors of Ref. [50], taking into consideration the optical uncertainty principle, also plotted the vortex spectrum for a topologically neutral beam *m* ¼ 0 and a beam with a small topological charge *m* ¼ 2 at the angle of aperture *α* ¼ 45°. The authors did not detect a maximum of spectra in the negative region of topological charges. As we presented above, the emergence of the second

*(a) Sketch of the hard-edged aperture (D) installed at the cylindrical lens plane (CL). Images LG and HG beams illustrate the astigmatic transformation of a single vortex beam. (b and c) Intensity distribution ℑ*ð Þ *r*, *φ*, *z of the vortex beam with m* ¼ *15 along length Z* ¼ *z=z0 perturbed by a hard-edged sector diaphragm with*

*m*,*n with charges m* ¼ *5 perturbed by the hard-edged aperture with the*

*m*,*n* for the initial topological charges *n* ¼ *m* ¼ 5 with *α* ¼

does not occur at any beam length *z* even at a very small sector angle *α*.

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

*4.1.1 The vortex spectrum*

vortex mode intensity *С*2

**Figure 6.**

**Figure 7.**

**175**

*α* ¼ *π=4: (b) theory, (c) experiment.*

*Distribution of vortex beams spectrum C*2

*sector adjacent angle β* ¼ *π* � *α, the solid line is spectrum outline.*

Further increasing the parameter *p* (Ω ¼ *const*) is accompanied by smoothing the curve ℓ*z*ð Þ *p* so that ℓ*z*jmax ! 2*p* except for the OAM dips at the integer values *m* ¼ *p*.

## **4. OAM, informational entropy and topological charge of truncated vortex beams**

A special case is represented by natural changes of the vortex spectrum due to external influences on the laser beam (diffraction by opaque obstacles, external interference in the beam, etc.). The simplest opaque obstacles are the sector, circular and annular apertures where we stop our attention using a simple technique of the vortex digital sorting.

#### **4.1 Sector aperture**

Let us consider propagation a scalar beam of Laguerre-Gaussian *LGm* 0 with a zero radial index *p* ¼ 0 and an azimuthal index (a topological charge) *m* through the rough regular sector aperture obstacle, that an angle shown in **Figure 1**. The edge of sector touches axis of the beam. The field of beam at the initial plane *z* ¼ 0 can be represented as [49].

$$\Psi\_m(r,\rho,a) = \left(\rho/w\right)^{|m|} e^{i m \cdot \rho} e^{-\rho^2/w^2} = r^{|m|} e^{i m \cdot \rho} e^{-r^2}, \quad a < \rho < 2\pi - a,\tag{21}$$

where *r* ¼ *ρ=w*. *w* is a beam waist radius at the plane *z* ¼ 0, *ρ* and *φ* represent polar coordinates. We rewrite the beam field (1) formed by the rigid-edges aperture with the half angle *α* as a sum of non-normalized Laguerre-Gauss beams *LGm p* in the form

$$\Psi\_m(r,\rho,a) = \sum\_{n=-\infty}^{\infty} \mathcal{C}\_{m,n}(a) L G\_0^n(r,\rho) = \sum\_{n=-\infty}^{\infty} \mathcal{C}\_{m,n}(a) r^{|n|} e^{in\cdot\rho} e^{-r^2},\tag{22}$$

where the beam amplitudes are

$$C\_{m,n}(a) = (-1)^{m-n} \Gamma\left(\frac{|m|+|n|}{2} + 1\right) \frac{\sin\left[ (m-n)(\pi-a) \right]}{m-n} / \left(\pi 2^{\frac{|m|-|n|}{2}} |n|!\right),\tag{23}$$

where Γð Þ *x* is a Gamma function. The terms in the series (23) with radial indices *p* 6¼ 0 disappear due to orthogonality of the LG modes.

The complex amplitude (22) can describe the spatial evolution of the perturbed beam if we replace *r* ! *r=σ* and multiply the sum by *σ* where *σ* ¼ 1 � *i z=z*0.

**Figure 6(b)** and **(c)** shows a computer simulation (**Figure 6(b)**) and experiment (**Figure 6(c)**) of the intensity distribution *ℑ*ð Þ¼ *r*, *φ*, *z* ΨΨ∗ along the vortex beam length with *m* ¼ 15 modulated by the sector aperture with *α* ¼ *π=*4. It is important to note, that the cut out dark sector from the beam intensity distribution does not lead to its self-healing. When the beam propagates, the dark sector rotates synchronously. We showed, both theoretically and experimentally, that beam self-healing does not occur at any beam length *z* even at a very small sector angle *α*.

## *4.1.1 The vortex spectrum*

j j *Cm*ð Þ *m*

*m* ¼ *p*.

**vortex beams**

the vortex digital sorting.

**4.1 Sector aperture**

represented as [49].

form

**174**

Ψ*m*ð Þ¼ *r*, *φ*, *α* ð Þ *ρ=w* j j *m e*

Ψ*m*ð Þ¼ *r*, *φ*, *α* X∞

*Сm*,*n*ð Þ¼ � *α* ð Þ1 *m*�*n*

where the beam amplitudes are

*n*¼�∞

Γ j j *m* þ j j *n*

*p* 6¼ 0 disappear due to orthogonality of the LG modes.

2 and Ω >10 induce the vortex avalanche (see **Figure 4(b)**). However, the avalanche type differs significantly from that in the previous case. Indeed, as well as in the above case, the perturbation provokes appearing resonance satellites with both positive and negative topological charges (TC). However, their values are always greater than those of the initial TC. In this case, energy is also pumped from the basic mode to the satellites but with higher topological charges under the condition of the resonance *m* ¼ *p*. In the OAM spectrum, it looks like the OAM burst at the resonance as shown in **Figure 5(c)** and **(d)** but the maxima ℓ*z* are at the region of fractional values *p*. As the scale parameter increases, the height of the flash

Further increasing the parameter *p* (Ω ¼ *const*) is accompanied by smoothing the curve ℓ*z*ð Þ *p* so that ℓ*z*jmax ! 2*p* except for the OAM dips at the integer values

**4. OAM, informational entropy and topological charge of truncated**

Let us consider propagation a scalar beam of Laguerre-Gaussian *LGm*

�*ρ*2*=w*2

*i m φe*

*Cm*,*n*ð Þ *α LGn*

2 þ 1

radial index *p* ¼ 0 and an azimuthal index (a topological charge) *m* through the rough regular sector aperture obstacle, that an angle shown in **Figure 1**. The edge of sector touches axis of the beam. The field of beam at the initial plane *z* ¼ 0 can be

> ¼ *r* j j *m e i m φe* �*r*2

where *r* ¼ *ρ=w*. *w* is a beam waist radius at the plane *z* ¼ 0, *ρ* and *φ* represent polar coordinates. We rewrite the beam field (1) formed by the rigid-edges aperture with the half angle *α* as a sum of non-normalized Laguerre-Gauss beams *LGm*

0ð Þ¼ *r*, *φ* X∞

� � sin ½ � ð Þ *m* � *n* ð Þ *π* � *α*

where Γð Þ *x* is a Gamma function. The terms in the series (23) with radial indices

The complex amplitude (22) can describe the spatial evolution of the perturbed

beam if we replace *r* ! *r=σ* and multiply the sum by *σ* where *σ* ¼ 1 � *i z=z*0.

*n*¼�∞

*Cm*,*n*ð Þ *α r*

*m* � *n = π* 2

j j *n e i n φe* �*r*2

A special case is represented by natural changes of the vortex spectrum due to external influences on the laser beam (diffraction by opaque obstacles, external interference in the beam, etc.). The simplest opaque obstacles are the sector, circular and annular apertures where we stop our attention using a simple technique of

0 with a zero

*p* in the

, (22)

, (23)

, *α*<*φ* <2*π* � *α*, (21)

j j� *m* j j *n* 2 j j *n* ! � �

decreases, and the resonance circuit broadens (**Figure 5(d)**).

*Nonlinear Optics - From Solitons to Similaritons*

We consider vortex beams spectra for topological charge that most clearly reflect properties of the sector perturbation. Computer simulation and experimental data of typical vortex spectra is shown in **Figure 7**. We revealed a clear maximum of vortex mode intensity *С*2 *m*,*n* for the initial topological charges *n* ¼ *m* ¼ 5 with *α* ¼ 45°. The short part of the beam intensity redistributed symmetric at the neighboring vortices orders. But, if we increase the angle *α*, we see a violation of the symmetric distribution of intensity among the vortex modes. There formed a second maximum in the spectrum range of negative topological charges *n* <0. It is important to note, that the authors of Ref. [50], taking into consideration the optical uncertainty principle, also plotted the vortex spectrum for a topologically neutral beam *m* ¼ 0 and a beam with a small topological charge *m* ¼ 2 at the angle of aperture *α* ¼ 45°. The authors did not detect a maximum of spectra in the negative region of topological charges. As we presented above, the emergence of the second

#### **Figure 6.**

*(a) Sketch of the hard-edged aperture (D) installed at the cylindrical lens plane (CL). Images LG and HG beams illustrate the astigmatic transformation of a single vortex beam. (b and c) Intensity distribution ℑ*ð Þ *r*, *φ*, *z of the vortex beam with m* ¼ *15 along length Z* ¼ *z=z0 perturbed by a hard-edged sector diaphragm with α* ¼ *π=4: (b) theory, (c) experiment.*

#### **Figure 7.**

*Distribution of vortex beams spectrum C*2 *m*,*n with charges m* ¼ *5 perturbed by the hard-edged aperture with the sector adjacent angle β* ¼ *π* � *α, the solid line is spectrum outline.*

intensity maximum is possible only at sufficiently large angles of overlap of beam *α* and also large topological charges *m* of the initial singular beam.

*4.1.4 The topological charge*

numbers *С*2

**Figure 10a**)

**177**

According to Berry [42], the topological charge of the vortex array is defined as the difference of the fluxes of vortex trajectories through the beam cross section

> 2 ð*π*

> > 0

Berry showed [42] that the vortex beam TC has always integer values when the spiral phase plate is perturbed equal to the integer value of the unperturbed plate step. We investigated the changes of the TC under the beam sector perturbations [52] and perturbation of the holographic grating [47]. Note that all further calculations are based on the requirement that the perturbation does not introduce changes in the mode phases. We performed a series of computer TC estimations of the perturbed beam (20) for various initial TC. The following restrictions were used. As the spectra of vortices in **Figure 7** show, the squared mode amplitudes quickly tend to zero as their TC increases. Therefore, we can restrict ourselves to a finite mode

*mn* � 10�3. The half-width beam radius can be estimated at the maxi-

replaced the limit *r* ! ∞ by *r* ! 3*rN*. The results of computer simulation are shown in **Figure 8(c)**. It can be seen from the figure that the total topological charge *M* at *α* ¼ 3*π:*4 remains equal to the initial TC of the unperturbed beam *M* ¼ *m*. Computer calculations for sector aperture angles *α* ¼ *π=*4, *π=*2, 30*π=*31 show similar results. Therefore, we assume that, in the general case of arbitrary angles, the initial TC will be preserved under sector perturbations. We also examined the constancy condition TC under the perturbation of the holographic grating considered above and made sure that this condition is strictly satisfied. Moreover, we traced the maxima displacement of the spectral satellites in the vortex avalanche in **Figure 4** [47] and convinced (see **Figure 8(d)**) that their TC cannot take fractional values.

The problem of the birth and annihilation of phase singularities has been considered as far back as at the beginning of the last century in connection with the peculiarities of light diffraction at the edges of the half-plane or lenses and microobjectives of telescopes and microscopes (see [53] and references therein). As a rule, the discussion came down to the technique of suppressing the corresponding aberrations. In this section, we focus on the digital vortex sorting after beam diffraction by the circular and annular apertures addressing the Shannon entropy

Note also that recently, special attention has been paid to studies on increasing

indices *n*, but a constant topological charge *m* [54]. Such modes are sorted using the holographic grating techniques and *a single phase screen* [55]. In this section, we will

*∂φ*Ψ*m=*Ψ*mdφ*

9 =

;*:* (25)

j j *m =*2 p . Therefore, in Eq. (25) we

*n* modes with various radial

0 beam at the *z* ¼ 0 plane with

2 � �, *R*<*r*<*R* þ *h:* (26)

taking into account vortex directions and "weights," and is calculated as

8 < :

2*π* Im lim*r*!∞

*M* ¼ 1

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

mum intensity of the initial vortex beam *rM* ¼ ffiffiffiffiffiffiffiffiffiffiffi

**4.2 Circular and annular apertures**

problem of the diffracted combined vortex beams.

the information capacity of optical channels due to *LGm*

We consider the perturbation of the vortex *LGm*

Ψ*m* ¼ *r*

j j *m e*

focus on the digital sorting of LG modes with various radial indices.

a ring aperture so that its complex amplitude is written in the form (see also

*i m φ* exp �*r*

#### *4.1.2 The orbital angular momentum*

The OAM per photon of a complex perturbed beam can be calculated in accordance with Eq. (23). The mode amplitudes are given by the normalized field Ψ*m* of the Laguerre-Gauss beams so that the squared amplitudes *С*2 *m*,*n* and *С*2 *m*,*n* are obeyed a simple relation *С*2 *m*,*n* ¼ 2�j j�*n* 2 j j *n* !*С*2 *m*,*n*. The changes of OAM ℓ*z*ð Þ *α*, *m* after increasing the aperture angle *α* (decreasing the adjacent angle *β* ¼ *π* � *α*) is presented by **Figure 7(a)** for topological charges *m* ¼ 5, *m* ¼ 10 and *m* ¼ 15 of the initial vortex beam.

The OAM is practically unchanged a wide range of angles and remains almost equal to the initial OAM, despite the rapid increase in the number of vortex states (see **Figure 8(a)**). After the second spectral maximum is formed in the negative region of topological charges **Figure 7**, there is a sharp decrease of the OAM. The OAM is equal to zero already at the angle *β* ≈2°.

#### *4.1.3 Informational entropy (Shannon's entropy)*

The normalized squared amplitude *С*2 *n* ∈ð Þ 0, 1 in the Eq. (23) can be treated as a conditional probability *P n*ð Þ *=m* of finding out a vortex beam in the state j i *n* among of 2*N* states (see, e.g., [51] and references therein). This approach to counting a number of vortex states can be used for the Shannon Entropy [51] and is written as

$$H\_{\rm I} = -\sum\_{n=0}^{N} P(n/m) \log\_2 P(n/m) = -\sum\_{n=0}^{N} \mathcal{C}\_n^2(a,m) \log\_2 \mathcal{C}\_n^2(a,m) > 0. \tag{24}$$

The Shannon entropy (24) characterizes the amount of uncertainty (randomness) that arises when a perturbation acts on a vortex beam. For example, in the case of the same amplitudes of the mode beams, the Shannon entropy of *HI* ¼ 2 means that to remove the uncertainty it is necessary to expend 1 bit of information.

In **Figure 8(b)** presents the dependence of the Shannon Entropy *HI* on the angle *α* for various of number topological charges *m*. In a broad range of aperture angles 0<*α*< 7*π=*8, the informational entropy *HI* increases equally for any topological charges *m* while the OAM ℓ*z* does not change (see **Figure 8(a)**) which corresponds to the same changes in the spectrum of vortex states in **Figure 7**.

#### **Figure 8.**

*Computer simulation of (a) the OAM* ℓzð Þ *α*, m *and (b) information entropy HI*ð Þ *α*, *m for the initial topological charges, m* ¼ *5*, *10*, *15; (c) transformation of the vortex TC Mof the perturbed beam via variation TC of the initial beam; (d) maximal TC mmax of the spectral vortex satellites via perturbation Ω of the perturbed holographic gratings; (solid lines); crosslets* ð Þ � *, circlets* ð Þ ⊙ *and squares (□) are correspondent experimental data.*

#### *4.1.4 The topological charge*

intensity maximum is possible only at sufficiently large angles of overlap of beam *α*

The OAM per photon of a complex perturbed beam can be calculated in accordance with Eq. (23). The mode amplitudes are given by the normalized field Ψ*m* of

ing the aperture angle *α* (decreasing the adjacent angle *β* ¼ *π* � *α*) is presented by **Figure 7(a)** for topological charges *m* ¼ 5, *m* ¼ 10 and *m* ¼ 15 of the initial vortex

The OAM is practically unchanged a wide range of angles and remains almost equal to the initial OAM, despite the rapid increase in the number of vortex states (see **Figure 8(a)**). After the second spectral maximum is formed in the negative region of topological charges **Figure 7**, there is a sharp decrease of the OAM. The

conditional probability *P n*ð Þ *=m* of finding out a vortex beam in the state j i *n* among of 2*N* states (see, e.g., [51] and references therein). This approach to counting a number of vortex states can be used for the Shannon Entropy [51] and is written as

*N*

*n*¼0 *C*2

The Shannon entropy (24) characterizes the amount of uncertainty (randomness) that arises when a perturbation acts on a vortex beam. For example, in the case of the same amplitudes of the mode beams, the Shannon entropy of *HI* ¼ 2 means that to remove the uncertainty it is necessary to expend 1 bit of information. In **Figure 8(b)** presents the dependence of the Shannon Entropy *HI* on the angle *α* for various of number topological charges *m*. In a broad range of aperture angles 0<*α*< 7*π=*8, the informational entropy *HI* increases equally for any topological charges *m* while the OAM ℓ*z* does not change (see **Figure 8(a)**) which corresponds

*m*,*n* and *С*2

*m*,*n*. The changes of OAM ℓ*z*ð Þ *α*, *m* after increas-

*n* ∈ð Þ 0, 1 in the Eq. (23) can be treated as a

*n*ð Þ *α*, *m* >0*:* (24)

*n*ð Þ *α*, *m* log 2*C*2

*m*,*n* are obeyed

and also large topological charges *m* of the initial singular beam.

the Laguerre-Gauss beams so that the squared amplitudes *С*2

*P n*ð Þ *=m* log 2*P n*ð Þ¼� *=m* X

to the same changes in the spectrum of vortex states in **Figure 7**.

*Computer simulation of (a) the OAM* ℓzð Þ *α*, m *and (b) information entropy HI*ð Þ *α*, *m for the initial topological charges, m* ¼ *5*, *10*, *15; (c) transformation of the vortex TC Mof the perturbed beam via variation TC of the initial beam; (d) maximal TC mmax of the spectral vortex satellites via perturbation Ω of the perturbed holographic gratings; (solid lines); crosslets* ð Þ � *, circlets* ð Þ ⊙ *and squares (□) are correspondent*

j j *n* !*С*2

*m*,*n* ¼ 2�j j�*n* 2

OAM is equal to zero already at the angle *β* ≈2°.

*4.1.3 Informational entropy (Shannon's entropy)*

The normalized squared amplitude *С*2

*4.1.2 The orbital angular momentum*

*Nonlinear Optics - From Solitons to Similaritons*

a simple relation *С*2

*HI* ¼ �X *N*

**Figure 8.**

**176**

*experimental data.*

*n*¼0

beam.

According to Berry [42], the topological charge of the vortex array is defined as the difference of the fluxes of vortex trajectories through the beam cross section taking into account vortex directions and "weights," and is calculated as

$$M = \frac{1}{2\pi} \text{Im} \left\{ \lim\_{r \to \infty} \left\{ \left. \partial\_{\rho} \Psi\_{m} / \Psi\_{m} d\rho \right. \right\} \right\}. \tag{25}$$

Berry showed [42] that the vortex beam TC has always integer values when the spiral phase plate is perturbed equal to the integer value of the unperturbed plate step. We investigated the changes of the TC under the beam sector perturbations [52] and perturbation of the holographic grating [47]. Note that all further calculations are based on the requirement that the perturbation does not introduce changes in the mode phases. We performed a series of computer TC estimations of the perturbed beam (20) for various initial TC. The following restrictions were used. As the spectra of vortices in **Figure 7** show, the squared mode amplitudes quickly tend to zero as their TC increases. Therefore, we can restrict ourselves to a finite mode numbers *С*2 *mn* � 10�3. The half-width beam radius can be estimated at the maximum intensity of the initial vortex beam *rM* ¼ ffiffiffiffiffiffiffiffiffiffiffi j j *m =*2 p . Therefore, in Eq. (25) we replaced the limit *r* ! ∞ by *r* ! 3*rN*. The results of computer simulation are shown in **Figure 8(c)**. It can be seen from the figure that the total topological charge *M* at *α* ¼ 3*π:*4 remains equal to the initial TC of the unperturbed beam *M* ¼ *m*. Computer calculations for sector aperture angles *α* ¼ *π=*4, *π=*2, 30*π=*31 show similar results. Therefore, we assume that, in the general case of arbitrary angles, the initial TC will be preserved under sector perturbations. We also examined the constancy condition TC under the perturbation of the holographic grating considered above and made sure that this condition is strictly satisfied. Moreover, we traced the maxima displacement of the spectral satellites in the vortex avalanche in **Figure 4** [47] and convinced (see **Figure 8(d)**) that their TC cannot take fractional values.

#### **4.2 Circular and annular apertures**

The problem of the birth and annihilation of phase singularities has been considered as far back as at the beginning of the last century in connection with the peculiarities of light diffraction at the edges of the half-plane or lenses and microobjectives of telescopes and microscopes (see [53] and references therein). As a rule, the discussion came down to the technique of suppressing the corresponding aberrations. In this section, we focus on the digital vortex sorting after beam diffraction by the circular and annular apertures addressing the Shannon entropy problem of the diffracted combined vortex beams.

Note also that recently, special attention has been paid to studies on increasing the information capacity of optical channels due to *LGm n* modes with various radial indices *n*, but a constant topological charge *m* [54]. Such modes are sorted using the holographic grating techniques and *a single phase screen* [55]. In this section, we will focus on the digital sorting of LG modes with various radial indices.

We consider the perturbation of the vortex *LGm* 0 beam at the *z* ¼ 0 plane with a ring aperture so that its complex amplitude is written in the form (see also **Figure 10a**)

$$\Psi\_m = r^{|m|} e^{i \, m \cdot p} \exp\left(-r^2\right), \ \ R < r < R + h. \tag{26}$$

If the axis of the annular aperture coincides with the beam axis then the perturbation excites only LG beams with the same TC (*m* ¼ *const*) but different radial indices *n* (this follows from the condition of orthogonally of LG beams). The perturbed beam field can be represented as an expansion over LG beams

$$\Psi\_m = \sum\_{n=0}^{\infty} \mathcal{C}\_{m,n} r^m L\_n^m \left( 2r^2 \right) e^{i \, m \varphi} e^{-r^2}. \tag{27}$$

by a annular aperture, which are the result of transformations in the vortex spectra in **Figure 9(a–c)**. The presented results show that any external interference in the beam immediately affects the uncertainty of the vortex beam state, the magnitude

Another interesting feature of the axial aperture action is manifested under the combined vortex beam perturbation consisting of two vortex beams with the same values but different TC signs (*m* and �*m*). If the amplitudes of the beams are the same, then the complex amplitude of such a combined beam is described by Eq. (28) when replacing the phase factor exp ð Þ! *imφ* 2 cos *mφ*. **Figure 10(b–d)** illustrates the intensity distribution of such perturbed beams at the measurement plane. Each of these beams receives the same perturbations, regardless of the sign of their TC. Even if the beam amplitudes are different, they receive the same amount of the vortex state uncertainty. However, this apparent indistinguishability of the modes can be easily detected experimentally due to opposite phase circulation of

*m*,*n and the Shannon entropy HI*ð Þ *h of the perturbed singular beams with different TC m:*

*(I) the circular aperture; (II) the annular aperture; (Id) the Shannon entropy HI via the aperture radius R0; (IIh) the Shannon entropy HI via the annular width, curves 1 and 3 correspond to the TC m* ¼ *0*, *m* ¼ *2*, *m* ¼ *4,*

*(a) Sketch of the annular aperture. (b–d) Intensity distributions of the singular beams with TC m* ¼ *5 and m* ¼ �*5 perturbed by the circular aperture with the radius (b) R0* ¼ *0:5, (c) R0* ¼ *0:1 and (d) the annular*

of which can be estimated by measuring informational entropy.

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

the fields, and the modes can be sorted out in different memory cells.

**Figure 9.**

**Figure 10.**

**179**

*aperture with R0* ¼ *0:5*, *h* ¼ *0:2.*

*The vortex spectra C2*

*respectively; callouts: corresponding intensity distributions.*

The LG vortex modes amplitudes restricted by the ring of thickness *h*, we find for the field difference the beams passing simultaneously through the circular aperture of radius *R* þ *h* and the opaque disk of radius *R*. As a result we obtain

$$\begin{split} C\_{m,n}(R) &= \int\_0^\infty \Psi\_m \left( L G\_n^m \right)^\* r dr \\ &= \left( 2(R+h)^2 \right)^{m+1} e^{-2(R+h)^2} \,\_1F\_1 \left( 1-n, m+2; -2(R+h)^2 \right) / (m+1)m! \\ &- \left( 2R^2 \right)^{m+1} e^{-2R^2} \,\_1F\_1 \left( 1-n, m+2; -2R^2 \right) / (m+1)m!, \end{split} \tag{28}$$

where 1*F*1 is a confluent hypergeometric function and we used the integrals from [56]. The expression (28) together with Eq. (27) allows covering three cases: (1) a circular aperture, *R* ¼ 0, *h* ¼ *R*0; (2) an opaque disk *R* ¼ *R*0, *h* ¼ 0; and (3) an annular aperture *R*, *h* 6¼ 0. An important property of mode amplitudes (28) is that the perturbation of a singular beam with a defined TC via these types of axial apertures does not excite vortex modes with other TC. This means that the OAM ℓ*z* does not change due to such a perturbation process. Does this mean that the perturbed vortex beam completely restores its initial properties during propagation, i.e., possesses the self-healing effect? We will peer into this process carefully.

Recall that under the action of axial perturbation, vortex modes with new topological charges do not appear in the perturbed beam. Therefore, the digital sorting of vortex modes in a perturbed beam can be carried out in accordance with Eq. (15) for the nondegenerate case (see Section 3.1). In this case, as a function of moments *Mp*,*q*ð Þ *r*, *φ* should choose the forms *Mp*,*q* ¼ sin *rLk p* 2*r*2 ð Þ or *Mp*,*q* ¼ cos *rLk p* 2*r*2 ð Þ while in the intensity distribution *ℑ*ð Þ *r*, *φ* use the complex amplitude (26). Variation of the indices *p* and *k* enables us to obtain a closed system of linear equations for the squared amplitudes *C*2 *m*,*n*.

The experimental results of measuring the vortex spectra are shown in **Figure 9**, where the average values of the squared amplitudes are plotted along the ordinate axis. A characteristic feature of the dependences *C*2 *m*,*n*ð Þ *n* is long spectral tails, which are omitted in the figures, but which make a significant contribution to the calculation of information entropy. Truncation of a topologically neutral beam (*m* ¼ 0) with a circular aperture in **Figure 9(Ia)** leads to overlapping many side rings, resulting in a wave-like form of the spectral tail. In the perturbed vortex beams shown in **Figure 9(Ibc)**, a broadening of the vortex spectrum and a decrease in the tail amplitudes are observed. Such characteristic features of the vortex spectra insert significant uncertainty into information carried by the vortex beam, which is represented as the dependence of information entropy *HI* on the aperture radius *R*0 in **Figure 9(Id)**. It is noteworthy that even small variations in the aperture radius lead to changes in the entropy *HI* that grow with increasing TC. Similar changes in the entropy *HI*ð Þ *h* shown in **Figure 9(IId)** occur when a vortex beam is perturbed

### *Digital Sorting of Optical Vortices in Perturbed Singular Beams DOI: http://dx.doi.org/10.5772/intechopen.91419*

by a annular aperture, which are the result of transformations in the vortex spectra in **Figure 9(a–c)**. The presented results show that any external interference in the beam immediately affects the uncertainty of the vortex beam state, the magnitude of which can be estimated by measuring informational entropy.

Another interesting feature of the axial aperture action is manifested under the combined vortex beam perturbation consisting of two vortex beams with the same values but different TC signs (*m* and �*m*). If the amplitudes of the beams are the same, then the complex amplitude of such a combined beam is described by Eq. (28) when replacing the phase factor exp ð Þ! *imφ* 2 cos *mφ*. **Figure 10(b–d)** illustrates the intensity distribution of such perturbed beams at the measurement plane. Each of these beams receives the same perturbations, regardless of the sign of their TC. Even if the beam amplitudes are different, they receive the same amount of the vortex state uncertainty. However, this apparent indistinguishability of the modes can be easily detected experimentally due to opposite phase circulation of the fields, and the modes can be sorted out in different memory cells.

#### **Figure 9.**

If the axis of the annular aperture coincides with the beam axis then the perturbation excites only LG beams with the same TC (*m* ¼ *const*) but different radial indices *n* (this follows from the condition of orthogonally of LG beams). The perturbed beam field can be represented as an expansion over LG beams

> *mLm n* 2*r* 2 � �*e*

The LG vortex modes amplitudes restricted by the ring of thickness *h*, we find for the field difference the beams passing simultaneously through the circular aperture of radius *R* þ *h* and the opaque disk of radius *R*. As a result we obtain

*i mφe* �*r*2

1*F*1 1 � *n*, *m* þ 2; �2ð Þ *R* þ *h* 2 � �

*p* 2*r*2 ð Þ or *Mp*,*q* ¼ cos *rLk*

*m*,*n*ð Þ *n* is long spectral tails, which

1*F*1 1 � *n*, *m* þ 2; �2*R*2 � �*=*ð Þ *m* þ 1 *m*!,

where 1*F*1 is a confluent hypergeometric function and we used the integrals from [56]. The expression (28) together with Eq. (27) allows covering three cases: (1) a circular aperture, *R* ¼ 0, *h* ¼ *R*0; (2) an opaque disk *R* ¼ *R*0, *h* ¼ 0; and (3) an annular aperture *R*, *h* 6¼ 0. An important property of mode amplitudes (28) is that the perturbation of a singular beam with a defined TC via these types of axial apertures does not excite vortex modes with other TC. This means that the OAM ℓ*z* does not change due to such a perturbation process. Does this mean that the perturbed vortex beam completely restores its initial properties during propagation, i.e., possesses the self-healing effect? We will peer into this process carefully.

Recall that under the action of axial perturbation, vortex modes with new topological charges do not appear in the perturbed beam. Therefore, the digital sorting of vortex modes in a perturbed beam can be carried out in accordance with Eq. (15) for the nondegenerate case (see Section 3.1). In this case, as a function of moments

while in the intensity distribution *ℑ*ð Þ *r*, *φ* use the complex amplitude (26). Variation of the indices *p* and *k* enables us to obtain a closed system of linear equations for the

The experimental results of measuring the vortex spectra are shown in **Figure 9**, where the average values of the squared amplitudes are plotted along the ordinate

are omitted in the figures, but which make a significant contribution to the calculation of information entropy. Truncation of a topologically neutral beam (*m* ¼ 0) with a circular aperture in **Figure 9(Ia)** leads to overlapping many side rings, resulting in a wave-like form of the spectral tail. In the perturbed vortex beams shown in **Figure 9(Ibc)**, a broadening of the vortex spectrum and a decrease in the tail amplitudes are observed. Such characteristic features of the vortex spectra insert significant uncertainty into information carried by the vortex beam, which is represented as the dependence of information entropy *HI* on the aperture radius *R*0 in **Figure 9(Id)**. It is noteworthy that even small variations in the aperture radius lead to changes in the entropy *HI* that grow with increasing TC. Similar changes in the entropy *HI*ð Þ *h* shown in **Figure 9(IId)** occur when a vortex beam is perturbed

*:* (27)

*=*ð Þ *m* þ 1 *m*!

*p* 2*r*2 ð Þ

(28)

Ψ*m* ¼ X∞

*r dr*

*e* �2ð Þ *R*þ*h* 2

*e* �2*R*2

*Mp*,*q*ð Þ *r*, *φ* should choose the forms *Mp*,*q* ¼ sin *rLk*

*m*,*n*.

axis. A characteristic feature of the dependences *C*2

*Сm*,*n*ð Þ¼ *R*

squared amplitudes *C*2

**178**

∞ð

Ψ*m LGm n* � � ∗

*Nonlinear Optics - From Solitons to Similaritons*

¼ 2ð Þ *R* þ *h* 2 � �*m*þ1

� 2*R*2 � �*m*þ1

0

*n*¼0

*Cm*,*nr*

*The vortex spectra C2 m*,*n and the Shannon entropy HI*ð Þ *h of the perturbed singular beams with different TC m: (I) the circular aperture; (II) the annular aperture; (Id) the Shannon entropy HI via the aperture radius R0; (IIh) the Shannon entropy HI via the annular width, curves 1 and 3 correspond to the TC m* ¼ *0*, *m* ¼ *2*, *m* ¼ *4, respectively; callouts: corresponding intensity distributions.*

#### **Figure 10.**

*(a) Sketch of the annular aperture. (b–d) Intensity distributions of the singular beams with TC m* ¼ *5 and m* ¼ �*5 perturbed by the circular aperture with the radius (b) R0* ¼ *0:5, (c) R0* ¼ *0:1 and (d) the annular aperture with R0* ¼ *0:5*, *h* ¼ *0:2.*

## **5. Conclusions**

We examined the technique of digital sorting of vortex modes that makes it possible to measure in real time the vortex spectrum (squared amplitudes and initial phases) including radial indices, OAM, and informational entropy of perturbed singular beams. The considered approach is based on the measurement of intensity moments of higher orders and a digital solution of a linear equations system that eliminates the cuts and gluing of the beam wavefront without losing information on the modes initial phases. Moreover, the digital vortex spectrum also enables us to restore the initial combined beam and, correcting parameters of the spectral modes, to improve its characteristics.

The digital approach has been tested on vortex beams free of wave defects perturbed both by local defects of holographic gratings responsible for the beam generation and by the sectorial, circular and annular aperture. We revealed that a local perturbation of the holographic grating near the central forked defect causes bursts and dips in the OAM spectrum. The depth and height of the spectral dips and bursts are controlled by the parameters of the holographic grating and can vary over a wide range. The perturbation inserted by the sector aperture is regulated by the sector angle. Over a wide range of sector angles, the beam OAM remains almost unchanged. However, when the sector angle is relatively large, so that most of the light flux is cut off by the aperture, the optical uncertainty principle begins to act, and the OAM sharply decreases to almost zero that is accompanied by a rapid growth of the Shannon entropy. At the same time, the beam topological charge remains unchanged for any sectorial perturbations. The axial perturbation via a circular and annular aperture does not change either the OAM or topological charge. However, a wide range of Laguerre-Gauss modes with the same topological charges but different radial indices leads to a rapid increase in information entropy as the pupil of the circular aperture or the ring thickness of the annular aperture decreases. This allows not only to estimate the noise level in the optical information transmission line, but also to record external interference in the information flow. We also note that the digital sorting of optical vortices in a perturbed light flux opens up broad prospects for its employment for medical express-diagnostics of skin diseases, since, for example, this allows us to detect slight changes in the vortex spectrum of a laser beam scattered by inflamed or dehydrated skin areas.

## **Acknowledgements**

The authors are grateful to E. Abramochkin (Samara Branch of the Lebedev Physical Institute, Russian Academy of Sciences, Samara, Russia) for a useful discussion of the mathematical approach. The reported study was funded by RFBR according to the research project № 19-29-01233.

**Author details**

**181**

Simferopol, Republic of Crimea, Russia

provided the original work is properly cited.

Alexander Volyar\*, Mikhail Bretsko, Yana Akimova and Yuriy Egorov

\*Address all correspondence to<EMAIL_ADDRESS>
*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

Physics and Technology Institute, V.I. Vernadsky Crimean Federal University,

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

*Digital Sorting of Optical Vortices in Perturbed Singular Beams DOI: http://dx.doi.org/10.5772/intechopen.91419*

**5. Conclusions**

to improve its characteristics.

*Nonlinear Optics - From Solitons to Similaritons*

**Acknowledgements**

**180**

according to the research project № 19-29-01233.

We examined the technique of digital sorting of vortex modes that makes it possible to measure in real time the vortex spectrum (squared amplitudes and initial phases) including radial indices, OAM, and informational entropy of perturbed singular beams. The considered approach is based on the measurement of intensity moments of higher orders and a digital solution of a linear equations system that eliminates the cuts and gluing of the beam wavefront without losing information on the modes initial phases. Moreover, the digital vortex spectrum also enables us to restore the initial combined beam and, correcting parameters of the spectral modes,

The digital approach has been tested on vortex beams free of wave defects perturbed both by local defects of holographic gratings responsible for the beam generation and by the sectorial, circular and annular aperture. We revealed that a local perturbation of the holographic grating near the central forked defect causes bursts and dips in the OAM spectrum. The depth and height of the spectral dips and bursts are controlled by the parameters of the holographic grating and can vary over a wide range. The perturbation inserted by the sector aperture is regulated by the sector angle. Over a wide range of sector angles, the beam OAM remains almost unchanged. However, when the sector angle is relatively large, so that most of the light flux is cut off by the aperture, the optical uncertainty principle begins to act, and the OAM sharply decreases to almost zero that is accompanied by a rapid growth of the Shannon entropy. At the same time, the beam topological charge remains unchanged for any sectorial perturbations. The axial perturbation via a circular and annular aperture does not change either the OAM or topological charge. However, a wide range of Laguerre-Gauss modes with the same topological charges but different radial indices leads to a rapid increase in information entropy as the pupil of the circular aperture or the ring thickness of the annular aperture decreases. This allows not only to estimate the noise level in the optical information transmission line, but also to record external interference in the information flow. We also note that the digital sorting of optical vortices in a perturbed light flux opens up broad prospects for its employment for medical express-diagnostics of skin diseases, since, for example, this allows us to detect slight changes in the vortex spectrum of a laser beam scattered by inflamed or dehydrated skin areas.

The authors are grateful to E. Abramochkin (Samara Branch of the Lebedev Physical Institute, Russian Academy of Sciences, Samara, Russia) for a useful discussion of the mathematical approach. The reported study was funded by RFBR

## **Author details**

Alexander Volyar\*, Mikhail Bretsko, Yana Akimova and Yuriy Egorov Physics and Technology Institute, V.I. Vernadsky Crimean Federal University, Simferopol, Republic of Crimea, Russia

\*Address all correspondence to<EMAIL_ADDRESS>
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Allen L, Beijersbergen MW, Spreeuw RJC, Woerdman JP. Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes. Physical Review A. 1992; **45**(11):8185-8189. DOI: 10.1103/ PhysRevA.45.8185

[2] Berry M. Nature's optics and our understanding of light. Contemporary Physics. 2015;**56**(1):2-16. DOI: 10.1080/ 00107514.2015.97162

[3] Mendenhall MJ, Nunez AS, Martin RK. Human skin detection in the visible and near infrared. Applied Optics. 2015;**51**(2):10559-10570. DOI: 10.1364/AO.54.010559

[4] Michel AP, Liakat S, Bors K, Gmachl CF. In vivo measurement of mid-infrared light scattering from human skin. Biomedical Optics Express. 2013;**4**(4):520-530. DOI: 10.1364/ BOE.4.000520

[5] Chorvat D, Chorvatova A. Multiwavelength fluorescence lifetime spectroscopy: A new approach to the study of endogenous fluorescence in living cells and tissues. Laser Physics Letters. 2009;**6**:175-193. DOI: 10.1002/ lapl.200810132

[6] Gbur GJ. Singular Optics. New York: CRC Press; 2017. DOI: 10.1201/<PHONE_NUMBER>260

[7] Lavery MPJ, Berkhout GCG, Courtial J, Padgett MJ. Measurement of the light orbital angular momentum spectrum using an optical geometric transformation. Journal of Optics. 2013; **13**:064006. DOI: 10.1088/2040-8978/ 13/6/064006

[8] Shen Y, Yang X, Qi R, Wan Z, Fu X, Gon M. Recent advances on tunable vortex beam devices for biomedical applications. BJSTR. 2018;**9**(3):2/5-4/5. DOI: 10.26717/BJSTR.2018.09.001801

[9] Otaka H, Shimakura H, Motoyoshi I. Perception of human skin conditions and image statistics. Journal of the Optical Society of America A. 2019; **36**(9):1609-1616. DOI: 10.1364/ JOSAA.36.001609

Perez-de-Tejada H, editor. Vortex Dynamics and Optical Vortices. Rijeka: IntechOpen; 2017. DOI: 10.5772/66343

*DOI: http://dx.doi.org/10.5772/intechopen.91419*

*Digital Sorting of Optical Vortices in Perturbed Singular Beams*

angular momentum states of light. Physical Review Letters. 2010;**105**: 153601. DOI: PhysRevLett.105.153601

[25] Volyar A, Bretsko M, Akimova Y, Egorov Y. Measurement of the vortex spectrum in a vortex-beam array without cuts and gluing of the

wavefront. Optics Letters. 2018;**43**(22): 5635-5638. DOI: 10.1364/OL.43.005635

Gopinath JT, Siements KE. Quantitative measurement of the orbital angular momentum of light with a single, stationary lens. Optics Letters. 2016;**41**: 5019-5022. DOI: 10.1364/OL.41.005019

[27] Alperin N, Siemens ME. Angular momentum of topologically structured darkness. Physical Review Letters. 2017; **119**:203902. DOI: PhysRevLett.119.

Porfirev AP. Methods for determining the orbital angular momentum of a laser beam. Computer Optics. 2019;**43**(1): 42-53. DOI: 10.18287/2412-6179-

[29] Soifer VA, Golub MA. Laser Beam Mode Selection by Computer-Generated Holograms. Boca Raton: CRC Press; 1994. p. 224. ISBN: 978-0-8493-2476-5

[30] Gibson G, Courtial J, Padgett MJ, Vasnetsov M, Pas'ko V, Barnett SM, et al. Free-space information transfer using light beams carrying orbital angular momentum. Optics Express.

[31] Berry MV. Paraxial beams of spinning light. In: Proceedings of the SPIE International Conference Singular Optics. Vol. 3487. 1998. pp. 6-11. DOI:

[32] D'errico A, D'amelio R, Piccirillo B, Cardano F, Marrucc L. Measuring the complex orbital angular momentum

[28] Kotlyar VV, Kovalev AA,

203902

2019-43-1-42-53

2004;**12**:5448-5456

10.1117/12.317704

spectrum and spatial mode

[26] Alperin SN, Niederiter RD,

[17] Segev M, Soljačić M, Dudley JM. Fractal optics and beyond. Nature Photonics. 2012;**6**(4):209-210. DOI:

[18] Sheppard CJR. Fractal model of light scattering in biological tissue and cells. Optics Letters. 2007;**32**(2):142-144. DOI: 10.1364/OL.32.000142

[20] Si K, Gong W, Sheppard CJR. Model for light scattering in biological tissue and cells based on random rough nonspherical particles. Applied Optics. 2009;**48**(6):1153-1157. DOI: 10.1364/

[21] Khonina SN, Kotlyar VV, Soifer VA, Paakkonen P, Turunen J. Measuring the light field orbital angular momentum using DOE. Optical Memory and Neural

[22] Khonina SN, Kotlyar VV, Soifer VA, Jefimovs K, Turunen J. Generation and selection of laser beams represented by a superposition of two angular harmonics. Journal of Modern Optics. 2004;**51**:

Networks. 2001;**10**(4):241-255

761-773. DOI: 10.1080/<PHONE_NUMBER>8235551

28067

**183**

[23] Khonina SN, Kazanskiy NL, Soifer VA. In: Yasin M, Harun SW, Arof H, editors. Optical Vortices in a Fiber: Mode Division Multiplexing and Multimode Self-Imaging, Recent Progress in Optical Fiber Research. Rijeka: IntechOpen; 2012. DOI: 10.5772/

[24] Berkhout GCG, Lavery MPJ, Courtial J, Beijersbergen MW,

Padgett MJ. Efficient sorting of orbital

10.1038/nphoton.2012.71

[19] Xu M, Alfano RR. Fractal mechanisms of light scattering in biological tissue and cells. Optics Letters. 2005;**30**(22):3051-3053. DOI:

10.1364/OL.30.003051

AO.48.001153

[10] Ney M, Abdulhalim I. Does human skin truly behave as an array of helical antennae in the millimeter and terahertz wave ranges? Optics Letters. 2010; **35**(19):3180-3182. DOI: 10.1364/ OL.35.003180

[11] Gbur G, Korotkova O. Angular spectrum representation for the propagation of arbitrary coherent and partially coherent beams through atmospheric turbulence. Journal of the Optical Society of America. A. 2007;**24**: 745-752. DOI: 10.1364/JOSAA.24.000745

[12] Soifer VA, Korotkova O, Khonina SN, Shchepakina EA. Vortex beams in turbulent media. Review Computer Optics. 2016;**40**(5):605-621. DOI: 10.18287/2412-6179-2016-40-5-605-624