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]):
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]:
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
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
eiφ and
e−iφ we substitute the field (2) with the cosine argument
reiφ and
re−iφ respectively, we construct a vector field with an infinite number of the polarization singularities:
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]:
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):
with
ψ = arctan(
z/
z0) being the Gouy phase.
Then the complex amplitude (4) can be rewritten as
or in a matrix form
where the matrix
converts phase singularities (optical vortices) into polarization singularities (radial polarization), while the matrix
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:
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:
Both real and imaginary parts should be zero, thus we obtain
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,
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
In the far field,
ψ → π/2, and the intensity nulls reside on the vertical axis in points with the coordinates
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
with
<<
α0,
w0. Since
ψ →
π/2 in the far field, then in these points the rotated coordinates (6) read as
Therefore,
,
, and the field amplitude is equal to:
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):
In the same way, we can derive the distribution of the longitudinal component of the SAM density:
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,
and, therefore, cos 2
x+ = –1 and cosh 2
y+ = 1, i.e.,
y+ = 0 and
x+ = (
π/2)(2
p + 1), with
p being an integer number. Then, the coordinates of the C-points are
Similarly, left circular polarization appears in points with
Sz = −
I. Coordinates of these points are equal to
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
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:
with
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:
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
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:
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
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 = 5
z0, 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
and the intensity distribution is
Using an identity
, we obtain
Finally, using the formulae for trigonometric and hyperbolic functions with a double argument, we obtain an expression for the intensity distribution (20):
References
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Indebetouw, G. Optical vortices and their propagation. J. Mod. Opt. 1993, 40, 73–87. [Google Scholar] [CrossRef]
- Abramochkin, E.; Volostnikov, V. Spiral-type beams. Opt. Commun. 1993, 102, 336–350. [Google Scholar] [CrossRef]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- Freund, I. Polarization singularity indices in Gaussian laser beams. Opt. Commun. 2002, 201, 251–270. [Google Scholar] [CrossRef]
- 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]
- Zhan, Q. Cylindrical vector beams: From mathematical concepts to applications. Adv. Opt. Photon. 2009, 1, 1–57. [Google Scholar] [CrossRef]
- Tidwell, S.C.; Ford, D.H.; Kimura, W.D. Generating radially polarized beams interferometrically. Appl. Opt. 1990, 29, 2234–2239. [Google Scholar] [CrossRef] [PubMed]
- 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]
- 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]
- 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]
| 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. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).