Next Article in Journal
Universally Grasping Objects with Granular—Tendon Finger: Principle and Design
Next Article in Special Issue
Analytical Gas Sensing in the Terahertz Spectral Range
Previous Article in Journal
WPD-Enhanced Deep Graph Contrastive Learning Data Fusion for Fault Diagnosis of Rolling Bearing
Previous Article in Special Issue
Holographic Fabrication of 3D Moiré Photonic Crystals Using Circularly Polarized Laser Beams and a Spatial Light Modulator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities

by
Alexey A. Kovalev
1,2,
Victor V. Kotlyar
1,2 and
Anton G. Nalimov
1,2,*
1
Image Processing Systems Institute of the RAS—Branch of FSRC “Crystallography & Photonics” of the RAS, 151 Molodogvardeyskaya St., 443001 Samara, Russia
2
Technical Cybernetics Department, Samara National Research University, 34 Moskovskoe Shosse, 443086 Samara, Russia
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(7), 1470; https://doi.org/10.3390/mi14071470
Submission received: 20 June 2023 / Revised: 13 July 2023 / Accepted: 20 July 2023 / Published: 21 July 2023
(This article belongs to the Special Issue Photon-Driven Technologies: Exploring the Latest Advancements)

Abstract

:
Elements of micromachines can be driven by light, including structured light with phase and/or polarization singularities. We investigate here a paraxial vector Gaussian beam with an infinite number of polarization singularities residing evenly on a straight line. The intensity distribution is derived analytically and the polarization singularities are shown to exist only in the initial plane and in the far field. The azimuthal angle of the polarization singularities is shown to increase in the far field by π/2. We obtain the longitudinal component of the spin angular momentum (SAM) density and show that it is independent of the azimuthal angle of the polarization singularities. Upon propagation in free space, an infinite number of C-points is generated, where polarization is circular. We show that the SAM density distribution has a shape of four spots, two with left and two with right elliptic polarization. The distance to the transverse plane with the maximal SAM density decreases with decreasing distance between the polarization singularities in the initial plane. Generating such alternating areas with positive and negative SAM density, despite linear polarization in the initial plane, manifests the optical spin Hall effect. Application areas of the obtained results include designing micromachines with optically driven elements.

1. Introduction

Properties of light beams, and, in particular, optical vortices, can differ in near and far fields of diffraction. For instance, in [1], fractional-order optical vortices were studied in the near field. Such vortices contain chains of alternating ± first-order vortices, which disappear in the far field. As was demonstrated in [2], an optical vortex generated by a refractive spiral phase plate has an asymmetric shape in the Fresnel diffraction zone. Another work [3] investigated near-field diffraction of a Gaussian beam by fork gratings with different topological charges, and it was found that the generated optical vortices reside on spiral lines around the ± first diffraction orders. In far-field diffraction, transformation was found of the initially fractional topological charge [4,5]. In addition, in the far field (in the lens focus) of optical vortices, the spin Hall effect arises [6,7], i.e., alternating pairs of areas are generated in the focus with left and right circular polarization, despite linear polarization of the initial beam. This effect was discovered both for homogeneous linear polarization [8] and for inhomogeneous polarization (cylindrical) [9] of the initial field.
A natural generalization of a vortex light beam is a beam with several vortices. A seminal work with multiple vortices [10] investigated Gaussian beams with optical vortices located in the transverse plane arbitrarily. These beams are described by finite products with the number of multipliers equal to the number of vortices. Another work [11] describes propagation-invariant light fields with an arbitrary envelope analytical function in a closed form, whose zeros define positions of vortices in the beam. Based on [10], a light field can be constructed that has phase singularities residing on a circle [12]. In our paper [13], we investigated a similar field, but with polarization singularities on a circle. Recently we found that the spin Hall effect also arises in such fields [14], but in the Fresnel diffraction zone, rather than in the far field.
Besides the finite number of optical vortices, the approach from [11] allows for obtaining a light field with an infinite array of optical vortices, residing on a straight line [15]. Such fields have an infinite topological charge [16], can be generated by non-coaxial superposition of two tilted Gaussian beams [17], and are identified by density of the optical vortices, rather than by the topological charge, which can be measured interferometrically [15].
In this work, we study a vectorial Gaussian beam where, instead of phase singularities, an infinite number of polarization singularities reside on a straight line. The polarization singularities index (Poincaré–Hopf index) [18,19] of such a beam is shown to be also infinite. We found that the azimuthal angle of the polarization singularities [20] increases in the far field by π/2, i.e., initial radial polarization is converted to azimuthal and vice versa. It is demonstrated that when such a beam propagates in space, two pairs of areas are generated in the Fresnel zone with nonzero density of the longitudinal component of the spin angular momentum (SAM), despite linear polarization in the initial plane; i.e., the optical spin Hall effect arises.

2. Paraxial Light Fields with an Infinite Number of Phase or Polarization Singularities

In [11], the following solution to the paraxial Helmholtz equation has been obtained (Equation (17) in [11]):
E r , φ , z = 1 q exp r 2 q w 0 2 f r e i φ q w 0 ,
where (r, φ, z) are the cylindrical coordinates, w0 is the waist radius of the Gaussian beam, q = 1 + iz/z0, and f(ξ) is an arbitrary entire analytical function. The field (1) does not change its intensity structure on propagation in space. It is only widened |q| = [1 + (z/z0)]1/2 times and rotated around the optical axis by an angle arg(q) = arctan(z/z0).
This general expression allows for obtaining a solution of the paraxial Helmholtz equation with infinite or a finite number of optical vortices. For instance, if f(ξ) = cos(w0ξ/α0), then the vortices reside evenly along a straight line [15]:
E r , φ , z = 1 q exp r 2 q w 0 2 cos r e i φ α 0 q .
This field is an example of light fields with an infinite topological charge [16]. In the initial plane of such a field, optical vortices reside in points with the Cartesian coordinates xp = α0(π/2 + πp), y = 0 with p being integer numbers.
It has long been known [21] that cylindrical polarization can be represented as a superposition of ± first-order optical vortices with opposite circular polarization. For the Jones vectors, such representation can be written as
cos φ + δ sin φ + δ = 1 2 exp i φ + i δ 1 i + 1 2 exp i φ i δ 1 i .
where δ is the azimuthal angle of cylindrical polarization (if δ = 0 or if δ = π/2, polarization is respectively radial or azimuthal) [20].
Then, if we use the same Jones vectors, but instead of the multipliers e and e we substitute the field (2) with the cosine argument re and re respectively, we construct a vector field with an infinite number of the polarization singularities:
E x , y , z = 1 2 q W 0 exp x 2 + y 2 q w 0 2 × exp i δ cos x + i y α 0 q 1 i + exp i δ cos x i y α 0 q 1 i ,
where (x, y) are the Cartesian coordinates in the transverse plane and W0 is a multiplier introduced for normalizing the beam energy (i.e., in order to make it equal to the unit). This multiplier can be obtained from an expression for the energy of the scalar field (2) [15]:
W 0 = π w 0 2 2 cosh w 0 2 2 α 0 2 .
Superposition (4) consists of two beams. Upon propagation in space, one of them is rotated clockwise and the other counterclockwise. Note that such a field (4) does not retain the intensity structure during propagation; that is, it is not structurally stable, because there are optical vortices with different signs in it.
For a compact description of such propagation, we introduce two rotated coordinate systems (Figure 1):
x ± = x cos ψ ± y sin ψ / α 0 q , y ± = y cos ψ x sin ψ / α 0 q ,
with ψ = arctan(z/z0) being the Gouy phase.
Then the complex amplitude (4) can be rewritten as
E x , y , z = 1 2 q W 0 exp x 2 + y 2 q w 0 2 × exp i δ cos x + + i y + 1 i + exp i δ cos x i y 1 i ,
or in a matrix form
E x , y , z = 1 2 q W 0 exp x 2 + y 2 q w 0 2 cos δ sin δ sin δ cos δ 1 1 i i cos x + + i y + cos x i y
where the matrix
S = 1 1 i i
converts phase singularities (optical vortices) into polarization singularities (radial polarization), while the matrix
R = cos δ sin δ sin δ cos δ
rotates the strength vectors by the azimuthal angle δ.
As seen from the matrix representation, the azimuthal angle δ of cylindrical polarization does not affect the intensity distribution in an arbitrary transverse plane.
In the initial plane, polarization is linear in each point. Adopting an approach from [22], we can derive the polarization singularities index (Poincaré–Hopf index) [18] of the field (4). It is equal to the topological charge of the following scalar complex field:
E c = E x + i E y = e i δ q W 0 exp x 2 + y 2 q w 0 2 cos x + i y α 0 q ,
In [15], the topological charge of such scalar fields was shown to be infinite, and therefore the Poincaré–Hopf index of the vector field (4) is also infinite.

3. Intensity Nulls of Light Fields with an Infinite Number of Polarization Singularities

Here we obtain the intensity nulls of the field (4). Since the determinants of both matrices are nonzero, for the zero intensity at some point, the following conditions should be fulfilled:
cos x + + i y + = 0 , cos x i y = 0 .
Both real and imaginary parts should be zero, thus we obtain
cos x + cosh y + = 0 , sin x + sinh y + = 0 , cos x cosh y = 0 , sin x sinh y = 0 .
The hyperbolic cosine cannot be zero. Thus, from the first and third equation in (13) we determine that cos x+ = cos x = 0. This means that sin y+ ≠ 0 and sin x ≠ 0 and, therefore,
cos x + = 0 , y + = 0 , cos x = 0 , y = 0 .
Since y+ = y = 0, we determine that y cos ψ = x sin ψ = 0.
In the initial plane, ψ = 0, thus the intensity nulls reside in points with the coordinates
x = α 0 π / 2 + π p , y = 0 .
In the far field, ψ → π/2, and the intensity nulls reside on the vertical axis in points with the coordinates
y = α 0 q π / 2 + π p , x = 0 .
For finite distances z, cos ψ ≠ 0 and sin ψ ≠ 0. Therefore, conditions (14) cannot be fulfilled and the field (4) does not have the intensity nulls.
Now we consider vicinities of the intensity nulls in the far field, i.e., points with the coordinates
x = ρ cos θ , y = α 0 q π / 2 + π p + ρ sin θ ,
with ρ << α0, w0. Since ψπ/2 in the far field, then in these points the rotated coordinates (6) read as
x ± = ± π / 2 + π p ± ρ sin θ / α 0 q , y ± = ρ cos θ / α 0 q .
Therefore, x ± ± i y ± ± π / 2 + π p i ρ e ± i θ / α 0 q , cos x ± ± i y ± ± 1 p i ρ e ± i θ / α 0 q , and the field amplitude is equal to:
E ρ cos θ , α 0 q π / 2 + π p + ρ sin θ , z > > z 0 1 p 1 2 q W 0 i ρ α 0 q exp α 0 2 q 2 q w 0 2 exp i δ + i θ 1 i exp i δ i θ 1 i = 1 p 1 q W 0 ρ α 0 q exp α 0 2 q 2 q w 0 2 cos θ + δ + π / 2 sin θ + δ + π / 2 .
This means that the azimuthal angle δ in the far field increases by π/2, i.e., radial polarization is converted to azimuthal and vice versa.

4. Intensity and Spin Angular Momentum Density Distribution of Light Fields with an Infinite Number of Polarization Singularities

The intensity distribution of the field (4) is given by (Appendix A):
I x , y , z = E x x , y , z 2 + E y x , y , z 2 = = 1 4 q 2 W 0 exp 2 x 2 + y 2 q 2 w 0 2 cos 2 x + + cos 2 x + cosh 2 y + + cosh 2 y .
In the same way, we can derive the distribution of the longitudinal component of the SAM density:
S z x , y , z = 1 4 q 2 W 0 exp 2 x 2 + y 2 w 0 2 q 2 cos 2 x + cosh 2 y cos 2 x + cosh 2 y + .
Hence, both the SAM density and intensity distributions are independent of the azimuthal angle δ of cylindrical polarization. Note that although the scalar field (2) has an infinite topological charge, and the vector field (4) has an infinite polarization singularity index, the density of the SAM (21) has a finite value at each point, since it has a Gaussian envelope (hyperbolic cosines in (21) depend linearly on the argument). The expressions (20) and (21) allow for obtaining the coordinates of C-points of the field (4). C-points are points where the tilt of the major axis of the polarization ellipse in the beam cross section is undefined [18]. That is, C-points are points with circular polarization. For instance, right circular polarization appears in points where Sz = I. In these points, cos 2 x + + cosh 2 y + = 0 and, therefore, cos 2x+ = –1 and cosh 2y+ = 1, i.e., y+ = 0 and x+ = (π/2)(2p + 1), with p being an integer number. Then, the coordinates of the C-points are
x RCP y RCP = π 2 α 0 q 1 + 2 p cos ψ sin ψ .
Similarly, left circular polarization appears in points with Sz = −I. Coordinates of these points are equal to
x LCP y LCP = π 2 α 0 q 1 + 2 p cos ψ sin ψ .
The dynamics of the C-points explain the destruction of the polarization singularities after the initial plane and their reconstruction in the far field (Figure 2). Due to the splitting of left and right circular polarization, the beam (4) acquires nonzero SAM density upon propagation, and appearing areas with alternating SAM density manifest about the spin Hall effect.
We failed to derive exact expressions for the points of maximal SAM density. However, expressions (20) and (21) are simplified when w0α0. Indeed, using identities for the sums and differences of trigonometric and hyperbolic functions, we obtain
I x , y , z = 1 2 q 2 W 0 exp 2 x 2 + y 2 q 2 w 0 2 × cos 2 x cos ψ α 0 q cos 2 y sin ψ α 0 q + cosh 2 x sin ψ α 0 q cosh 2 y cos ψ α 0 q ,
S z x , y , z = 1 2 q 2 W 0 exp 2 x 2 + y 2 w 0 2 q 2 × sin 2 x cos ψ α 0 q sin 2 y sin ψ α 0 q + sinh 2 x sin ψ α 0 q sinh 2 y cos ψ α 0 q .
Products of two trigonometric or hyperbolic functions can be represented as a sum of four exponents. Thus, Equations (20) and (21) contain eight exponential terms. The first four terms do not exceed the value (8|q|2W0)−1, whereas the other four terms describe off-axis Gaussian beams:
G B ± ± x , y , z = 1 8 q 2 W 0 exp 2 x 2 + y 2 w 0 2 q 2 exp ± 2 x sin ψ α 0 q ± 2 y cos ψ α 0 q = = 1 8 q 2 W 0 exp 2 w 0 2 q 2 x x c , ± 2 + y y c , ± 2 + w 0 2 2 α 0 2 ,
with
x c , ± = ± w 0 2 q 2 α 0 sin ψ , y c , ± = ± w 0 2 q 2 α 0 cos ψ .
These terms achieve values (8|q|2W0)−1exp[(w0/α0)2/2]. Thus, if w0 >> α0, the first four terms can be neglected. Then, the intensity and the SAM density (1) are equal to the intensity of the four Gaussian beams:
I x , y , z = G B + + + G B + G B + + G B + ,
S z x , y , z = G B + + + G B G B + G B + .
If these four beams are far enough from each other (i.e., |xc+xc−| >> w0|q|, |yc+yc−| >> w0|q|), then these terms almost do not affect each other and, obviously, the points with the maximal SAM density coincide with the points of maximal intensity. In this case, the dependence of the maximal SAM density on the propagation distance z is given by
max x , y   S z x , y , z 1 8 q 2 W 0 exp w 0 2 2 α 0 2 .
Hence, the SAM density decreases with the propagation distance from the initial plane in a similar law, as does the intensity in the Gaussian beam center. However, polarization is everywhere linear in the initial plane, i.e., the maximal SAM density is zero. This means that it at first increases near the initial plane and then, when the Gaussian beams are split, begins to decrease, i.e., the maximal SAM density is achieved in near-field diffraction.

5. Identification of Light Fields with an Infinite Number of Polarization Singularities

In our work [15], we studied an analogy between scalar fields with an infinite topological charge and conventional circularly symmetric optical vortices. Similarly, we can consider circular fields with cylindrical polarization and fields with an infinite number of polarization singularities. We suppose that a field with cylindrical polarization is composed of two circularly polarized single-ringed LG beams with opposite topological charges:
E r , φ , 0 = r w 0 m exp r 2 w 0 2 + i m φ 1 i + r w 0 m exp r 2 w 0 2 i m φ 1 i = 2 r w 0 m exp r 2 w 0 2 cos m φ sin m φ .
where (r, φ) are the polar coordinates.
Such a field can be easily identified using a polarizer. If it transmits only one polarization, then a multi-petal intensity distribution is obtained, which allows for determining the order of cylindrical polarization by counting the petals (Figure 3a,b).
Similarly, registering the intensity of only one transverse component of the field (4) allows for determining the density of polarization singularities. Indeed, in the initial plane, the field (4) can be written as
E x x , y , 0 = 1 2 W 0 exp w 0 2 4 α 0 2 × exp x 2 + y + w 0 2 / 2 α 0 2 w 0 2 + exp x 2 + y w 0 2 / 2 α 0 2 w 0 2 cos x α 0 .
Hence, the intensity distribution of the x-component has vertical zero-intensity lines, whose frequency allows for determining the density of polarization singularities (Figure 3c,d).

6. Numerical Simulation

Shown in Figure 4 are intensity distributions of two beams (4) (with a different distance between the polarization singularities) in the initial plane, in uniform and in logarithmic color scale, as well as polarization directions. The intensity distributions were computed as I(r, φ, 0) = |Ex(r, φ, 0)|2 + |Ey(r, φ, 0)|2 using Equation (4), whereas the logarithmic distribution was computed as ln(10−18 + I(r, φ, 0)/max I(r, φ, 0)), where the constant 10−18 was introduced for avoiding the logarithm of zero in points with zero intensity. Polarization directions were computed by the formula arg(Ex(r, φ, 0) + iEy(r, φ, 0)).
As seen in Figure 4, the intensity distribution has a shape of two light spots located symmetrically relative to the horizontal coordinate axis, and there are polarization singularities with radial polarization, residing periodically on this axis. Due to the low intensity, they are not visible, but they can be seen on the intensity distribution in logarithmic color scale.
Figure 5 illustrates the intensity and SAM density distributions of the vector beam from Figure 4c,d in several transverse planes.
According to Figure 5, both light spots split, each into two spots, one of which shifts to the left and the other shifts to the right. It is also seen that the maximal SAM density decreases upon propagation, which is consistent with Equation (30).
Figure 6 depicts the intensity and SAM density distributions of the vector beam from Figure 4a,b in several transverse planes. In contrast to Figure 4c,d, light spots in Figure 4a,b are closer to each other; thus, they do not split so fast upon propagation into spots with left and right circular polarization, compared to Figure 5. Therefore, maximal SAM density in Figure 5 decreases immediately from z = z0/4 till z = 5z0, whereas in Figure 6 it at first increases at distances up to z = z0/2 and then decreases. The dependence of the maximal SAM density on the propagation distance is illustrated in Figure 7.
Figure 7 reveals that the maximal SAM density is achieved closer to the initial plane, when the distance between the polarization singularities decreases. This effect has a physical explanation since, as seen from Equation (4), decreasing value α0 leads not only to moving the light spots away from each other, but also to increasing space frequency along the axis y. Therefore, with decreasing α0, each light spot splits into two spots with opposite circular polarization faster.

7. Discussion

Usually, the spin Hall effect is considered in inhomogeneous media or in the presence of media interfaces, including metasurfaces. In most cases, the spin effect is associated with spin–orbit interaction. In this paper, we consider the propagation of a paraxial vortex laser vector beam in free space. The mechanism of appearance of the spin Hall effect in the considered beam is as follows. Beam (4) in the initial plane is a superposition of two vortex beams with left- and right-handed circular polarization, in which the centers of an infinite number of phase singularities are located at the same points on the horizontal axis, and have topological charges of different signs +1 and −1. The polarization of such a beam is linear in the initial plane, but as soon as the beam begins to propagate, optical vortices with left- and right-handed polarizations are separated and begin to rotate in different directions (vortices with a charge of −1 rotate clockwise, while those with a charge of +1 rotate counterclockwise). Therefore, regions with different spins are formed in the beam cross section, which leads to the Hall effect. The presence of an infinite number of screw dislocations in the initial beam is not fundamental; the number of vortices can be finite. Below, for comparison, we present an analysis of some works in which paraxial vortex beams were also considered, but with a finite number of polarization singularity points [14], with a fractional topological charge [23], and with vortex beams in a crystal [24].
In our recent work [14], we considered a paraxial laser beam with a finite number of polarization singularity points which are uniformly distributed over the beam cross section on a circle of some radius and in which the polarization is not defined. At all other points in the beam cross section in the initial plane, the polarization is linear. When such a beam propagates, each point of the polarization singularity in the beam cross section splits into two regions with elliptical polarization of a different sign, which move around the circle in different directions until they unite again at some distance and form points of the polarization singularity. The polarization becomes linear around these points. In this paper, we considered another paraxial beam, which is a superposition of optical vortices with left- and right-handed circular polarization and with an infinite number of phase singularities of different signs. In the initial plane, the centers of phase singularities for both beams coincide and are situated on the horizontal axis. In this case, the beam has a linear polarization at each point. When such a beam propagates in free space, then since the beams in superposition have topological charges of different signs, the centers of the phase singularities diverge and begin to rotate around the optical axis. In this case, the polarization in such parted centers became circular of different signs (C-points). In far-field diffraction, C-points of different signs coincide and form V-points, in which the polarization is undefined. Therefore, our previous work [14] is similar to this work in terms of the research method, but the considered beams are completely different.
In [23], a vector paraxial beam with an initial linear polarization and an initial fractional topological charge was also considered. The authors called this beam the erf-Gaussian beam. When it propagates in free space, a lot of optical vortices of different signs appear in the beam cross section, and the polarization near the intensity boundaries is elliptical of different signs. In fact, in [23], the authors showed that there is a spatial separation of the left- and right-handed circular polarizations in asymmetric paraxial beams with an initial linear polarization, which leads to the spin Hall effect.
In [24], a paraxial vortex laser beam propagating with a tilt to the axis of a birefringent crystal was considered. It was shown that two beams (ordinary and extraordinary) are formed in the crystal, each of which is a superposition of two optical vortices with topological charges +1 and −1 and with circular polarizations of different signs. It was also shown that in the cross section of such a beam there are points with circular polarization (C-points) of different signs. However, in these works [23,24] the SAM and OAM of the considered beams were not calculated.

8. Conclusions

In this work, we have constructed a vector Gaussian beam with an infinite number of polarization singularities residing on a straight line. For such a beam, the intensity distribution was derived analytically, and it turned out that the polarization singularities appear only in the initial plane and in the far field. We found that the polarization singularities index (Poincaré–Hopf index) is infinite. After propagation from the initial plane to the far field, the azimuthal angle of polarization singularities increases by π/2, i.e., initial radial polarization is converted into azimuthal and vice versa. We obtained a distribution of the longitudinal component of the spin angular momentum density. Similar to the intensity distribution, it is independent of the azimuthal angle of polarization singularities. When such a vectorial field propagates in free space, an infinite number of C-points appears, where polarization is circular. The distance to the transverse plane with the maximal spin angular momentum density decreases with as the distance between the polarization singularities decreases in the initial plane. Generation of alternating areas with left and right circular polarization, despite linear polarization in the initial plane, manifests in the optical spin Hall effect. Application areas of the results obtained include designing micromachines for optically driving microscopic objects. The SAM causes particles to rotate around their centers of mass and engineering the SAM density distribution of the studied light field can allow simultaneous manipulation of an ensemble of four particles. Another application is optical information transmission, where the density of polarization singularities can be used for encoding the data.

Author Contributions

Conceptualization, A.A.K.; methodology, A.A.K.; software, A.A.K. and A.G.N.; validation, V.V.K.; formal analysis, V.V.K.; investigation, A.A.K. and V.V.K.; resources, V.V.K.; data curation, V.V.K.; writing—original draft preparation, A.A.K.; writing—review and editing, V.V.K.; visualization, A.A.K. and A.G.N.; supervision, A.A.K.; project administration, V.V.K.; funding acquisition, A.G.N. and V.V.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by the RUSSIAN SCIENCE FOUNDATION, grant number 23-12-00236 (theory), and the Ministry of Science and Higher Education of the Russian Federation within a government project of the FSRC “Crystallography and Photonics” RAS (simulation).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A. Derivation of the Intensity Distribution

Since the azimuthal angle δ of cylindrical polarization does not affect the intensity distribution, we suppose for simplicity that δ = 0. Then, the transverse field components are given by
E x x , y , z = 1 2 q W 0 exp x 2 + y 2 q w 0 2 cos x + + i y + + cos x i y , E y x , y , z = i 2 q W 0 exp x 2 + y 2 q w 0 2 cos x + + i y + cos x i y ,
and the intensity distribution is
I x , y , z = E x x , y , z 2 + E y x , y , z 2 = = 1 2 W 0 q 2 exp 2 x 2 + y 2 q 2 w 0 2 cos x + + i y + 2 + cos x i y 2 .
Using an identity cos x + i y = cos x cosh y i sin x sinh y , we obtain
I x , y , z = 1 2 W 0 q 2 exp 2 x 2 + y 2 q 2 w 0 2 × cos 2 x + cosh 2 y + + sin 2 x + sinh 2 y + + cos 2 x cosh 2 y + sin 2 x sinh 2 y .
Finally, using the formulae for trigonometric and hyperbolic functions with a double argument, we obtain an expression for the intensity distribution (20):
I x , y , z = 1 4 W 0 q 2 exp 2 x 2 + y 2 q 2 w 0 2 cos 2 x + + cos 2 x + cosh 2 y + + cosh 2 y .

References

  1. Matta, S.; Vayalamkuzhi, P.; Viswanathan, N.K. Study of fractional optical vortex beam in the near-field. Opt. Commun. 2020, 475, 126268. [Google Scholar] [CrossRef]
  2. Kotlyar, V.V.; Kovalev, A.A.; Stafeev, S.S.; Nalimov, A.G. An asymmetric optical vortex generated by a spiral refractive plate. J. Opt. 2013, 15, 025712. [Google Scholar] [CrossRef]
  3. Ikonnikov, D.A.; Myslivets, S.A.; Arkhipkin, V.G.; Vyunishev, A.M. Near-Field Evolution of Optical Vortices and Their Spatial Ordering behind a Fork-Shaped Grating. Photonics 2023, 10, 469. [Google Scholar] [CrossRef]
  4. Jesus-Silva, A.J.; Fonseca, E.J.S.; Hickmann, J.M. Study of the birth of a vortex at Fraunhofer zone. Opt. Lett. 2012, 37, 4552. [Google Scholar] [CrossRef]
  5. Wen, J.; Wang, L.; Yang, X.; Zhang, J.; Zhu, S. Vortex strength and beam propagation factor of fractional vortex beams. Opt. Express 2019, 27, 5893. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, S.Q.; Chen, S.Z.; Wen, S.C.; Luo, H.L. Photonic spin Hall effect: Fundamentals and emergent applications. Opto-Electron. Sci. 2022, 1, 220007. [Google Scholar] [CrossRef]
  7. Liu, S.; Qi, S.; Li, Y.; Wei, B.; Li, P.; Zhao, J. Controllable oscillated spin Hall effect of Bessel beam realized by liquid crystal Pancharatnam-Berry phase elements. Light Sci. Appl. 2022, 11, 219. [Google Scholar] [CrossRef]
  8. Kotlyar, V.V.; Kovalev, A.A.; Kozlova, E.S.; Telegin, A.M. Hall Effect at the Focus of an Optical Vortex with Linear Polarization. Micromachines 2023, 14, 788. [Google Scholar] [CrossRef]
  9. Kovalev, A.A.; Kotlyar, V.V. Spin Hall Effect of Double-Index Cylindrical Vector Beams in a Tight Focus. Micromachines 2023, 14, 494. [Google Scholar] [CrossRef]
  10. Indebetouw, G. Optical vortices and their propagation. J. Mod. Opt. 1993, 40, 73–87. [Google Scholar] [CrossRef]
  11. Abramochkin, E.; Volostnikov, V. Spiral-type beams. Opt. Commun. 1993, 102, 336–350. [Google Scholar] [CrossRef]
  12. Dennis, M.R. Rows of optical vortices from elliptically perturbing a high-order beam. Opt. Lett. 2006, 31, 1325–1327. [Google Scholar] [CrossRef] [Green Version]
  13. Kovalev, A.A.; Kotlyar, V.V. Tailoring polarization singularities in a Gaussian beam with locally linear polarization. Opt. Lett. 2018, 43, 3084–3087. [Google Scholar] [CrossRef]
  14. Kovalev, A.A.; Kotlyar, V.V.; Stafeev, S.S. Spin Hall Effect in the Paraxial Light Beams with Multiple Polarization Singularities. Micromachines 2023, 14, 777. [Google Scholar] [CrossRef]
  15. Kovalev, A.A.; Kotlyar, V.V. Propagation-invariant laser beams with an array of phase singularities. Phys. Rev. A 2021, 103, 063502. [Google Scholar] [CrossRef]
  16. Berry, M.V. Optical vortices evolving from helicoidal integer and fractional phase steps. J. Opt. A Pure Appl. Opt. 2004, 6, 259–268. [Google Scholar] [CrossRef]
  17. Vaity, P.; Aadhi, A.; Singh, R.P. Formation of optical vortices through superposition of two Gaussian beams. Appl. Opt. 2013, 52, 6652–6656. [Google Scholar] [CrossRef]
  18. Freund, I. Polarization singularity indices in Gaussian laser beams. Opt. Commun. 2002, 201, 251–270. [Google Scholar] [CrossRef]
  19. Wang, Q.; Tu, C.H.; Li, Y.N.; Wang, H.T. Polarization singularities: Progress, fundamental physics, and prospects. APL Photonics 2021, 6, 040901. [Google Scholar] [CrossRef]
  20. Zhan, Q. Cylindrical vector beams: From mathematical concepts to applications. Adv. Opt. Photon. 2009, 1, 1–57. [Google Scholar] [CrossRef]
  21. Tidwell, S.C.; Ford, D.H.; Kimura, W.D. Generating radially polarized beams interferometrically. Appl. Opt. 1990, 29, 2234–2239. [Google Scholar] [CrossRef] [PubMed]
  22. Kotlyar, V.V.; Kovalev, A.A.; Stafeev, S.S.; Nalimov, A.G.; Rasouli, S. Tightly focusing vector beams containing V-point polarization singularities. Opt. Laser Technol. 2022, 145, 107479. [Google Scholar] [CrossRef]
  23. Fadeyeva, T.; Alexeyev, C.; Rubass, A.; Volyar, A. Vector erf-Gaussian beams: Fractional optical vortices and asymmetric TE and TM modes. Opt. Lett. 2012, 37, 1397–1399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Fadeyeva, T.; Rubass, A.; Egorov, Y.; Volyar, A.; Swartzlander, G., Jr. Quadrefringence of optical vortices in a uniaxial crystal. J. Opt. Soc. Am. A 2008, 25, 1634. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Coordinate systems (x+, y+) and (x, y).
Figure 1. Coordinate systems (x+, y+) and (x, y).
Micromachines 14 01470 g001
Figure 2. Mechanism of destruction of the polarization singularities after the initial plane and of their reconstruction in the far field. In the initial plane, points with left and right circular polarization coincide (a). Then, on propagation, they rotate around the optical axis in different directions (b), and in the far field they merge again (c).
Figure 2. Mechanism of destruction of the polarization singularities after the initial plane and of their reconstruction in the far field. In the initial plane, points with left and right circular polarization coincide (a). Then, on propagation, they rotate around the optical axis in different directions (b), and in the far field they merge again (c).
Micromachines 14 01470 g002
Figure 3. Intensity distribution of the Ex component of two beams with cylindrical polarization (31) in the initial plane at m = 3 (a) and at m = 7 (b), as well of two beams with an infinite number of polarization singularities (4) at a different singularities density: α0 = w0/5 (c) and α0 = w0/10 (d).
Figure 3. Intensity distribution of the Ex component of two beams with cylindrical polarization (31) in the initial plane at m = 3 (a) and at m = 7 (b), as well of two beams with an infinite number of polarization singularities (4) at a different singularities density: α0 = w0/5 (c) and α0 = w0/10 (d).
Micromachines 14 01470 g003
Figure 4. Intensity distributions of two beams (4) in the initial plane (a,c) and the logarithm of this distribution (b,d), as well as the polarization directions (green arrows) for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, and distance between the polarization singularities πα0 = πw0/2 ≈ 1.57 mm (a,b) and πα0 = πw0/5 ≈ 0.63 mm (c,d).
Figure 4. Intensity distributions of two beams (4) in the initial plane (a,c) and the logarithm of this distribution (b,d), as well as the polarization directions (green arrows) for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, and distance between the polarization singularities πα0 = πw0/2 ≈ 1.57 mm (a,b) and πα0 = πw0/5 ≈ 0.63 mm (c,d).
Micromachines 14 01470 g004
Figure 5. Intensity (ae) and SAM density (fj) distributions of the beam (4) from Figure 4c,d in several transverse planes for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, distance between the polarization singularities πα0 = πw0/5 ≈ 0.63 mm, and propagation distances from the initial plane z = z0/4 (a,f), z = z0/2 (b,g), z = z0 (c,h), z = 2z0 (d,i), and z = 5z0 (e,j). White dots on the SAM density distributions indicate the positions of maxima, obtained by Equation (27). The numbers near the color scales denote the minimal and maximal values.
Figure 5. Intensity (ae) and SAM density (fj) distributions of the beam (4) from Figure 4c,d in several transverse planes for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, distance between the polarization singularities πα0 = πw0/5 ≈ 0.63 mm, and propagation distances from the initial plane z = z0/4 (a,f), z = z0/2 (b,g), z = z0 (c,h), z = 2z0 (d,i), and z = 5z0 (e,j). White dots on the SAM density distributions indicate the positions of maxima, obtained by Equation (27). The numbers near the color scales denote the minimal and maximal values.
Micromachines 14 01470 g005
Figure 6. Intensity (ae) and SAM density (fj) distributions of the beam (4) from Figure 4a,b in several transverse planes for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, distance between the polarization singularities πα0 = πw0/2 ≈ 1.57 mm, and propagation distances from the initial plane z = z0/4 (a,f), z = z0/2 (b,g), z = z0 (c,h), z = 2z0 (d,i), and z = 5z0 (e,j). White dots in the SAM density distributions denote the positions of maxima computed by Equation (27). The numbers near the color scales denote the minimal and maximal values.
Figure 6. Intensity (ae) and SAM density (fj) distributions of the beam (4) from Figure 4a,b in several transverse planes for the following parameters: wavelength λ = 532 nm, Gaussian beam waist radius w0 = 1 mm, distance between the polarization singularities πα0 = πw0/2 ≈ 1.57 mm, and propagation distances from the initial plane z = z0/4 (a,f), z = z0/2 (b,g), z = z0 (c,h), z = 2z0 (d,i), and z = 5z0 (e,j). White dots in the SAM density distributions denote the positions of maxima computed by Equation (27). The numbers near the color scales denote the minimal and maximal values.
Micromachines 14 01470 g006
Figure 7. Dependence of the maximal intensity and of the maximal SAM density on the propagation distance for the distance between the polarization singularities equal to πα0 = πw0/2 ≈ 1.57 mm (a), πα0 = πw0/5 ≈ 0.63 mm (b), and πα0 = πw0/10 ≈ 0.31 mm (c).
Figure 7. Dependence of the maximal intensity and of the maximal SAM density on the propagation distance for the distance between the polarization singularities equal to πα0 = πw0/2 ≈ 1.57 mm (a), πα0 = πw0/5 ≈ 0.63 mm (b), and πα0 = πw0/10 ≈ 0.31 mm (c).
Micromachines 14 01470 g007
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kovalev, A.A.; Kotlyar, V.V.; Nalimov, A.G. Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities. Micromachines 2023, 14, 1470. https://doi.org/10.3390/mi14071470

AMA Style

Kovalev AA, Kotlyar VV, Nalimov AG. Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities. Micromachines. 2023; 14(7):1470. https://doi.org/10.3390/mi14071470

Chicago/Turabian Style

Kovalev, Alexey A., Victor V. Kotlyar, and Anton G. Nalimov. 2023. "Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities" Micromachines 14, no. 7: 1470. https://doi.org/10.3390/mi14071470

APA Style

Kovalev, A. A., Kotlyar, V. V., & Nalimov, A. G. (2023). Spin Hall Effect in Paraxial Vectorial Light Beams with an Infinite Number of Polarization Singularities. Micromachines, 14(7), 1470. https://doi.org/10.3390/mi14071470

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop