Abstract
1 Introduction
The interplay between spin and orbital angular momenta of light beams results in complex polarization textures of light fields with optical properties important in imaging, metrology, and quantum technologies.1 For example, polarization variations appear in the structure of two-dimensional photonic spin-skyrmions at length-scales much smaller than the wavelength of light because, in contrast to the field and intensity variations, the polarization structure is not influenced by the diffraction of electromagnetic waves.2 Such polarization features often appear due to the spin–orbit interactions involving vector vortex beams and, in the case of evanescent fields, may be topologically protected by the optical spin-Hall effect.3 For three-dimensional (3D) free-space beams, such topological protection is not ensured and the polarization features may vary significantly upon beam propagation. In contrast, observations of robustness of polarization inhomogeneities in 3D structured light have been reported.4 Polarization singularities of optical fields and their relation to phase singularities5 has become an active field of research.6
Sign up for Advanced Photonics Nexus TOC. Get the latest issue of Advanced Photonics Nexus delivered right to you!Sign up now
In this paper, we show that the transverse size of certain polarization features of optical vortex beams is preserved independently of the diffraction of the beam. The effect is governed by the phase singularity in the cross section of the beam and arises due to the interplay of the longitudinal and transverse electromagnetic fields in the vector vortex. The important role of the longitudinal fields was previously emphasized in the context of absorption of twisted photons by atoms7,8 and optical vortex dichroism.9,10
While the studies of vortex beams are typically concerned with the properties of the light in the high-intensity ring, this work deals with the region within this ring around the center of the beam. The demonstrated effect can be attributed to the spin–orbit interaction and only appears in the case of antialigned spin and orbital angular momenta.
The paper is organized as follows. In Sec. 2, we derive analytic expressions for the longitudinal-to-transverse field ratio near the beam phase singularity and demonstrate its independence on the beam waist in a paraxial limit. In Sec. 3, we introduce a formalism for 3D optical field polarization and show, in an analytic model, that the transverse spatial profiles of the polarization features are independent of beam divergence due to diffraction and the beam focusing conditions but depend on the topological charge of the beam. Finally, in Sec. 4, we use full-wave numerical simulations beyond the paraxial limit to reveal the diffraction-independent polarization features and confirm the analytic results.
2 Optical Vortex Fields near the Phase Singularity
We initially consider a paraxial monochromatic Laguerre–Gauss beam with a topological charge and zero radial index, propagating in the -direction (see the Appendix C for a general case). The electric field components in the transverse () plane, with the position vector in cylindrical coordinates, are given as
For a conventional Gaussian beam [Eq. (1)] with , it follows that for the transverse field linearly polarized along the axis
For a vortex beam with , choosing appropriate expressions for in Eq. (1) for left-hand or right-hand circular (), linear, radial, or azimuthal polarization (see Appendix A) and keeping only lowest-power terms in —i.e., assuming —we obtain the following longitudinal-to-transverse field ratios:
3 Polarization of 3D Vortex Fields
An arbitrary complex 3D vector field can be expanded in terms of unit vectors in a Cartesian basis as . The same field can be represented in a helical basis as , where , , and . The polarization coherence matrix for electric optical fields, (an asterisk indicates a complex conjugation), is fully defined in terms of standard Stokes parameters (see Appendix B) only if the longitudinal component of the field is neglected. However, as shown in the previous section, the longitudinal field of optical vortices is not negligible and may even be dominant at certain regions across the wavefront. For this reason, a Stokes description becomes incomplete, and the formalism for field polarization has to include all three components of the field. The matrix is Hermitian by construction and can be fully defined by nine real parameters. Here, we will follow the convention previously adopted for description of polarization of spin-1 particles8,13
The corresponding vector and quadrupole polarization parameters, such as and , can be expressed in terms of the field amplitudes and as
In a limiting case of plane waves propagating in the direction, the polarization parameters defined above either become zero or reduce to Stokes parameters (Ref. 16 and Appendix B)
The ratio introduced in the previous section may be probed experimentally by measuring the polarization parameter . The independence of from the beam waist found in the previous section has immediate implications for which, as a result, maintains constant transverse spatial dimensions independently of beam divergence due to diffraction. Using Eq. (5) and the definitions from Eq. (8), -independent expressions can be obtained for for different beam polarizations
Figure 1.Polarization parameters for a focused Laguerre–Gaussian vortex beam with
The above results were obtained in a simplified analytical model for a paraxial optical vortex field. Next, we demonstrate the nondiffractive behavior of these polarization features using numerical simulations for a nonparaxial field.
4 Numerical Simulations and Discussion
We now outline a full-wave nonparaxial numerical approach. Any monochromatic electromagnetic field can be decomposed into a spectrum of plane-wave components with wave vectors lying on the -sphere of radius and hence , as follows:
For these exact field calculations, we employ the Laguerre–Gaussian vortex beam widely used in singular optics.11 The angular spectrum of a Laguerre–Gaussian vortex beam can be calculated by selecting the plane and performing a Fourier transform .20,21 The derivation of the plane wave polarization amplitudes from the paraxial Laguerre–Gaussian beam field is described in detail in Appendix C.
We can use Eq. (11) to calculate the full 3D electric field and plot required polarization parameters without any approximations. The first example is a linearly polarized nonparaxial vortex beam with and propagating along the positive axis (Fig. 1). The tight focusing creates a highly divergent beam. In this case, the , , and parameters are required to fully describe the polarization structure [please note that the only remaining nonzero polarization parameter is but its behavior in the plane is the same as that of in the plane; see Fig. 1(b)]. In contrast to the electric field, which diffracts naturally after a propagation distance of just a few wavelengths [Fig. 1(a)], all three polarization parameters in Figs. 1(b)–1(d) clearly show no divergence around the phase singularity that lies on the axis. The polarization structures remain invariant and extend far beyond the focal plane, in agreement with the analytical paraxial predictions for in Eq. (5), but numerically observed here beyond the paraxial approximation.
We now demonstrate the equivalent polarization properties of a circularly polarized vortex beam. Figure 2(a) shows the intensity of a circularly polarized vortex beam, similar to the previous case except with , so that the spin and orbital angular momenta are antialigned. The polarization parameters again reveal a nondivergent polarization structure around the beam axis [Figs. 2(b)–2(d)]. For a vortex beam with a beam waist of double the size and, therefore, weaker focusing, the nondiffractive polarization structure near the phase singularity remains unperturbed, with the (white) contour lying at in both Figs. 2(d) and 2(h), in good agreement with the analytical result from Eq. (10): . This structural invariance has been observed for all beam waists, independent of focusing.
Figure 2.Polarization parameters for a focused Laguerre–Gaussian vortex beam with
The main challenge in detecting this nondiffractive polarization property is the requirement to perform measurements in a region of space where the field intensity is weaker compared to its maximum. The ability to detect it is determined by the sensitivity of the detection apparatus, and can be mediated to some degree by the choice of wavelength, beam waist, how far the detection plane is from the focal plane, and what polarization structure is being investigated. Figure 3(a) shows the cross-sectional plots of the parameter for the linearly polarized vortex beam in Fig. 1 at different points along the axis. As before, we see nondiffractive behavior when near the beam centre, indicated by a horizontal dotted black line. Note how other features where at locations further away from the beam center are subjected to diffraction. The intensity of the normalized electric field in the same cross-sectional planes at the location of the nondiffracting polarization features, indicated by the vertical dotted black lines, is 37% of the peak intensity at the focal plane, and 2.8% at [Fig. 3(b)]. The polarization structure should, therefore, be easily detectable in the focal plane of the beam and measurable away from the focus. The intensity drop-off of a beam is dictated by the Rayleigh range, which is proportional to the square of the beam waist. However, increasing the beam waist reduces the intensity of the longitudinal field. The resulting optimization will depend on the measurement sensitivity and the desired application.
Figure 3.Cross sections of the linearly polarized vortex beams with
To experimentally verify the propagation-invariant polarization structures, a vortex beam will likely need to be focused using a lens with a defined numerical aperture (NA). The effect of a restricted NA was simulated by limiting the integration of the () plane in Eq. (11) from to . Figure 4(a) shows the nonzero component of the angular spectrum for a linearly polarized vortex beam with and , obtained using Eq. (16), which is equivalent to the back focal plane image of the beam. The phase singularity is clearly visible at . In Figs. 1–3, all the fields within the light line (indicated by a green line) are integrated using Eq. (11) to create the real-space field distribution. We now crop the angular spectrum down to a factor of with , as indicated by the red dotted line. Figure 4(b) shows the intensity profile of the beam after this NA restriction. When compared with the ideal beam in Fig. 1(b), the polarization parameter of the restricted beam in Fig. 4(c) reveals the same nondiffractive property near the phase singularity and only disturbances in the peripheral fields are observed. As the NA is reduced further (not shown), the beam waist widens but the behavior around the phase singularity is maintained.
Figure 4.(a) The angular spectrum of a linearly polarized (
One can proceed to add more imperfections or aberrations to the beam. A defect or a piece of dust on the focusing lens can perturb the beam. This can be approximated by deleting part of the angular spectrum. A lens can also introduce astigmatism to the beam. The effect can be roughly modeled by scaling and in the angular spectrum. Figure 4(d) shows the angular spectrum of the same vortex beam as Fig. 4(a) with a restricted NA of 0.5 but with these two additional perturbations applied. The angular spectrum is set to zero for and , and is transformed by , therefore, reciprocally stretching the beam in the direction. The nondiffractive nature of the polarization parameter near the phase singularity is preserved for such scattered focused beams with astigmatism [Fig. 4(f)]. We therefore conclude that the nondiffractive polarization structures within a vortex beam should be robust to a variety of experimental imperfections and experimentally observable in this respect.
When considering higher-order vortex beams, further complications can arise from beam imperfections, resulting in splitting the high-order vortex into multiple low-order vortices.22
5 Conclusions
We have studied the polarization of vector beams carrying optical angular momentum. We show the existence of polarization features within optical vortex beams that maintain constant transverse spatial dimensions independently of the beam divergence due to diffraction. The exact size of these vortex polarization structures is dictated by the presence of the longitudinal electric field in the beam, and such structures are expected for vortex beams of all topological charges. An analytical paraxial model predicts their presence in weakly focused beams and a numerical angular spectrum approach further extended this prediction to tightly focused beams, thereby proving applicability to all vortex beams. These polarization features are not affected by finite NAs and so should be experimentally measurable. It should be noted that the predicted nondiffractive polarization features have relatively small transverse dimensions of the order , centered on a low-intensity region of the optical vortex wavefront. Therefore, future measurements will require subwavelength resolution at low and increased sensitivity of the probe for larger values of .
The demonstrated effect allows one to pinpoint the position of a phase singularity with subwavelength accuracy independently of the size of a beam spot. This property may have useful applications in metrology, optical communications, optical networking, laser sensing, and radar operations.
6 Appendix A. Calculation of Polarization Parameters
In Table 1, we present analytic expressions for polarization parameters [Eq. (8)] calculated for different polarizations of optical vortex beams; is assumed. The transverse field is defined by Eq. (1), and the longitudinal field is obtained from combined with the paraxiality condition, Eq. (2), at radial positions near the beam axis (much smaller than the beam waist). As in the main text, we define the dimensionless radial parameter , which depends on the topological charge . The transverse polarization vector for different polarizations is given in terms of unit vectors in Cartesian or cylindrical coordinates.
Circular | Circular | Linear | Radial | Azimuthal | |
1 | 0 | 0 | 0 | ||
0 | 0 | ||||
0 | 0 | 0 | |||
0 | 0 | ||||
0 | 0 | 0 | |||
0 | 0 | 0 | |||
0 | 0 | 0 | 0 | ||
1 | 1 |
Table 1. Polarization parameters.
Let us demonstrate the derivation in the case of linearly polarized (along the axis) light, so that . The continuity equation , together with the paraxiality condition, gives us the longitudinal field from the transverse-field derivative: . After some algebra, the field derivative can be obtained as
The expressions for circular polarization () are simplified in cylindrical coordinates, for which radial components of the polarization vector and tensor are zero: , , , and .
7 Appendix B. Definition of Stokes Parameters
A polarization coherence matrix for two-dimensional electric fields is defined in terms of Stokes parameters as16
8 Appendix C. Decomposing the Angular Spectrum into a Polarization Basis
Here, we show how nonparaxial fields of a focused vortex beam are calculated using the angular spectrum approach. We start with the paraxial expression for a Laguerre–Gauss beam11
In the remainder of this section, we do not explicitly write the dependencies of the angular spectrum for ease of notation, but note that all the fields mentioned here are the spectra defined in the plane unless otherwise stated. All fields are assumed to be time-harmonic.
The angular spectrum of the total field can be represented in a Cartesian basis , which can be split into a transverse part and a longitudinal part . Similarly, in the polarization basis . Equating these, one can write the transverse part of the field as
9 Appendix D. Splitting of Higher-Order Vortices
The analytical model in Sec. 2 is concerned with ideal vortex beams of any order. However, when a vortex beam with is perturbed, the high-order vortex can split into multiple lower-order vortices.
The splitting of an vortex beam can be investigated with our nonparaxial angular spectrum method by employing the appropriate equations for and . The continuous integration of the angular spectrum is approximated by a summation of a finite number of plane waves and inevitably generates a numerically approximate beam, which tends toward the ideal case when the number of plane waves tends to infinity. Figures 5(a)–5(d) show the phase of the transverse electric field for a nonparaxial vortex beam linearly polarized along , constructed using the angular spectrum approach and integrating a finite number of plane waves. This integration introduces a small perturbation in the fields away from the ideal vortex beam and, therefore, promotes a splitting of the singularity into three distinct vortices. The separation distance among the three singularities is reduced by increasing the number of plane waves.
Figure 5.(a–d) The phase of
This split vortex beam can now be analyzed by simulating polarization parameters in Eq. (8). For an obtained, simulated with high , nondiffractive features in the polarization structure can be observed similar to what was previously seen for an vortex [cf. Figs. 5(e) and 1(c)]. When is high, the individual vortices are extremely close together, and a collective polarization structure (such as the white cylindrical contour) is present. When is reduced and the vortices separate slightly, the contour is warped [Fig. 5(f)]. With further reduction of , the singularities move far apart and do not exhibit a collective polarization structure; instead, three individual cylindrical contours are observed [Fig. 5(g)]. In all three cases, the polarization structures were found to be nondiffractive as they extend far beyond the divergent beam field intensity drop-off. In other words, an ideal vortex beam exhibits the polarization structure predicted in Sec. 2 with the analytical paraxial model, but this structure breaks down into three separate structures when the beam is imperfect, with each of the separate vortices carrying its own non-diffractive polarization structure.
Andrei Afanasev received his PhD in physics from Karazin Kharkiv National University in Ukraine. Currently, he is working as a Gus Weiss professor of physics and endowed chair in theoretical physics at the Department of Physics in George Washington University in Washington DC, United States. His research interests include quantum electrodynamics, nuclear physics, and quantum optics.
Jack J. Kingsley-Smith received his MSc and PhD degrees in physics from the Department of Physics at King’s College London. Currently, he is a postdoctoral research associate at the same institution and is working on optical forces and structured light.
Francisco J. Rodríguez-Fortuño received his BSc, MSc, and PhD degrees at Universitat Politecnica de Valencia, Spain, with research stays at University of Pennsylvania and King’s College London, where he obtained a permanent academic position in 2015, starting his research team. He was a PI of ERC Starting Grant PSINFONI and Co-I in EIC Pathfinder CHIRALFORCE. Currently, he is working as a reader at King’s College London Physics Department, researching on nanophotonics and new electromagnetic phenomena, with a focus on near-field effects.
Anatoly V. Zayats is chair in experimental physics and head of the photonics and nanotechnology at the Department of Physics, King’s College London. He is working as a codirector at London Centre for Nanotechnology and London Institute for Advanced Light Technologies. His research interests include nanophotonics and plasmonics, metamaterials and metasurfaces, electromagnetic field topology and optical spin-orbit coupling, and nonlinear and ultrafast optics.
References
[1] Y. Shen et al. Topological quasiparticles of light: optical skyrmions and beyond(2022).
[3] K. Y. Bliokh et al. Spin–orbit interactions of light. Nat. Photonics, 9, 796-808(2015).
[10] G. De Ninno et al. Photoelectric effect with a twist. Nat. Photonics, 14, 554-558(2020).
[11] B. E. A. Saleh, M. C. Teich. Fundamentals of Photonics(2007).
[14] D. A. Varshalovich, A. N. Moskalev, V. K. Khersonskii. Quantum Theory of Angular Momentum(1988).
[15] J. Eismann et al. Transverse spinning of unpolarized light. Nat. Photonics, 15, 156-161(2021).
[16] E. Collett. Field Guide to Polarization(2005).
Set citation alerts for the article
Please enter your email address