Abstract
The local propagation and the energy flux in structured optical fields are often associated with the Poynting vector. However, the local phase gradient (i.e., local wavevector) in monochromatic fields in free space is described by another fundamental quantity: the canonical momentum density. Distributions of the Poynting and canonical momentum densities can differ significantly from each other in structured fields. We examine the role of these quantities in the propagation and diffraction of structured optical fields, exemplified by various circularly polarized vortex beams carrying orbital angular momentum. We describe the canonical and Poynting momentum distributions in such beams, experimentally measure the local transverse momentum density by a Shack-Hartmann wavefront sensor, and investigate fine features of the diffraction of various vortex beams on a knife-edge aperture. In all cases, the measured local momentum density and local beam evolution are consistent with the canonical momentum distribution rather than the Poynting vector. Furthermore, we introduce the local angular velocity in vortex beams and determine the universal integral $\pi$ angle of azimuthal rotation in an arbitrary (yet circularly symmetric) propagating and diffracting vortex beam. Finally, we discuss the “supermomentum” and “backflow” effects; both of these phenomena are examples of superoscillations and are related to the properties of the canonical momentum. Our results reveal the profound role of the canonical momentum in the evolution of light and demonstrate the importance of distinguishing between it and the Poynting vector in structured light.
© 2024 Optica Publishing Group
1. INTRODUCTION
Structured (i.e., inhomogeneous) wave fields are ubiquitous in nature and modern wave-based devices. Such waves are attracting ever growing attention in optics [1] and other fields of wave physics [2]. For a proper understanding of structured waves and their potential applications, it is important to characterize their local dynamical and propagational properties. In this context, the momentum density, or energy flux density, or wave current becomes the chief quantity of interest. Ideally, this quantity should (i) be well defined theoretically, (ii) be measurable in experiment, (iii) determine the wave-induced force (momentum transfer) in local wave-matter interactions, and (iv) determine the propagational evolution of the wave field. However, as it often happens, instead of a single well-defined quantity there are several competing alternatives with different properties and specific advantages in different types of problems.
There is a vast literature on various aspects of electromagnetic momentum density, including numerous controversies and paradoxes [3–8]. Here we consider only one class of problems, namely, propagation of structured monochromatic light in free space. Most textbooks on optics and electromagnetism introduce the Poynting vector as the main and only quantity characterizing the momentum or energy-flux density. Typical pictures of the structure of an inhomogeneous optical field (e.g., focused [9], diffracted [10], or scattered [11]) show the electric-field intensity (energy-density) distribution and streamlines of the Poynting vector as if these were the fluid-like currents determining the propagational evolution of light.
More recently, an alternative canonical (also called “orbital”) momentum density was introduced for monochromatic structured optical fields [6–8]. Compared to the Poynting vector this quantity has some advantages and disadvantages. The main features of the Poynting momentum density are as follows:
- • it is universally well defined for arbitrary (including non-monochromatic and static) electromagnetic fields in any media (anisotropic, inhomogeneous, etc.);
- • it is not directly observable in standard optical experiments;
In turn, the main features of the canonical momentum density are as follows:
Note that the last item is closely related to the fluid-like picture of propagation of light (or quantum-mechanical waves) within the Madelung-Bohm treatment of wave equations [20,21]. The streamlines of the canonical momentum current correspond to the experimentally measurable Bohmian trajectories of photons [19,22].
Importantly, the differences between the Poynting and canonical momenta are only local; their integral values for localized fields coincide. Yet, they coincide locally for linearly polarized (i.e., spinless) fields. Still, these local differences can be crucial and cause a number of controversies, such as the circular-polarization-dependent components of the Poynting vector orthogonal to the wavevectors in the problem [8,23–26].
The dilemma between the Poynting-like (sometimes called “kinetic”) and canonical momentum densities is not unique for electromagnetic waves; it equally appears in quantum-mechanical [27,28], acoustic [29,30], elastic [31–33], and water-wave [34] fields. (Incidentally, the acoustic analogue of the Poynting vector for waves in solids and fluids was derived by Umov [35] 10 years before the seminal Poynting work [36].) For acoustic and water waves the canonical momentum is also known as “pseudomomentum” [37–39], and it determines the Stokes drift of the medium particles [34,40,41]. Thus, in this case the canonical momentum is directly connected to the observable mechanical motion and momentum density of the medium. Therefore, thorough understanding of the Poynting versus canonical momentum properties of optical fields is also important for other fields of wave physics.
Despite the great progress in the description of the Poynting versus canonical momentum dilemma in the past decade, it still causes active debates in the structured-waves community. Furthermore, many scientists outside of this community are still unfamiliar with the problems and alternatives to the Poynting-vector description of structured light. Motivated by this, here we revisit the momentum density and propagational evolution of structured light using one of the most familiar examples: vortex beams [7,12,13,42]. It is now well established that direct optomechanical methods using small probe particles or atoms measure local optical forces proportional to the canonical momentum density (while the contribution proportional to the Poynting vector is much weaker) [8,25,26,43–45]. Moreover, “quantum weak measurements” of “Bohmian trajectories” of photons also correspond to the canonical momentum density [19,22]. In contrast to optomechanical and standard weak-measurement methods, here we analyze more traditional optical techniques based on (A) the Shack-Hartmann wavefront sensor and (B) propagational diffraction phenomena.
The earlier studies of vortex beams using the Shack-Hartmann wavefront sensor [46,47] and diffraction of truncated beams [48–51] dealt with linearly polarized paraxial beams and hence could not sense the difference between the Poynting and canonical momentum densities. That is why all these works mentioned the Poynting vector as the measured quantity. Here we revisit these approaches for the case of circularly polarized vortex beams and show that in fact both the Shack-Hartmann method and diffraction evolution reveal the canonical momentum density. Furthermore, we analyze the streamlines of the canonical momentum density (i.e., Bohmian trajectories) and their relation to the diffraction phenomena: despite the claim that “the rotation angle of the trajectory is unrelated to the Gouy phase” [52], we show that the mean rotation angle over all the trajectories exactly corresponds to the Gouy phase in the case of the Laguerre-Gaussian (LG) beams and always produces the same $\pi$ rotation in the $z \in (- \infty ,\infty)$ range for any vortex beams. Finally, we briefly consider the relation of the Poynting-canonical dilemma to “supermomentum” (i.e., anomalously high local wavevector) near the vortex core [6,53] and “backflow” (i.e., local retrograde propagation) [54,55]. These phenomena are currently attracting great attention [56–61] and represent examples of superoscillations [19,62]. We argue that these phenomena are determined by the canonical momentum, while the Poynting vector cannot produce supermomentum, although can exhibit “false backflow.”
2. CANONICAL AND POYNTING MOMENTUM DENSITIES IN VORTEX BEAMS
The time-averaged energy density, Poynting momentum density, and canonical momentum density for monochromatic paraxial optical fields in free space can be written as [6–8,13,63]
We now consider the well-known example of the LG vortex beams with circular polarization. Assuming propagation along the $z$ axis and the zero radial index, the transverse components of their electric field can be written as [7,12,42]
Substituting Eq. (5) into Eqs. (1) and (3), and neglecting quadratic contribution from the longitudinal field ${E_z}$, we find the canonical momentum density and energy density distributions:
Figure 1(a) shows the distributions of the transverse components of the canonical and Poynting momentum densities, as well as of the energy density, in the focal $z = 0$ plane of the LG beams with $\ell = 1$ and $\sigma = \pm 1$. The $\sigma$-dependent difference between the azimuthal momentum components, Eqs. (8) and (9), depends on the radius $r$ and generally is not small: $|{\Pi _\varphi} - {P_\varphi}| \sim |{P_\varphi}|$. Moreover, for $\ell \sigma \lt 0$ and radii $r \gt {w_0}\sqrt {|\ell |}$, the directions of the canonical and Poynting azimuthal momentum densities are opposite: ${\rm sgn}({P_\varphi}) = \def\LDeqbreak{}- {\rm sgn}({\Pi _\varphi})$. Thus, the choice of the canonical or Poynting momentum densities crucially determines the picture of the local propagational evolution of structured light.
The streamlines of the momentum density, Eq. (8) or Eq. (9), are the curves whose tangent vectors constitute the corresponding momentum density field. These curves can be regarded as “Bohmian trajectories” of the light propagation [19,22,52,67], by analogy with the streamlines of the velocity field in fluid mechanics. Figures 1(b) and 1(c) show such streamlines $(r(z),\varphi (z))$ in the LG beams with $\ell = 1$ and $\sigma = \pm 1$ (see Appendix A for their analytical calculations). One can see that the canonical and Poynting streamlines differ dramatically from each other, except for the straight-line trajectory at the radial intensity maximum: $\partial W/\partial r = 0$, ${r_0} = {w_0}\sqrt {|\ell |/2}$. The difference between the canonical and Poynting momentum densities vanishes there, because, according to Eq. (4), it requires gradients of the spin density (proportional to the intensity for paraxial uniformly polarized beams). The rectilinear propagation of the intensity-maximum points in diffracting LG beams was also noticed in other studies [68,69]. Note also that in the $\ell \sigma \gt 0$ case, the interplay between the azimuthal and radial components of the Poynting vector results in globally rectilinear streamlines; see the middle panels in Figs. 1(b) and 1(c). The role of the momentum streamlines in the propagation and diffraction of vortex beams will be discussed in Section 4.
3. MEASUREMENTS OF THE TRANSVERSE MOMENTUM DENSITY BY SHACK-HARTMANN SENSOR
We performed experimental measurements of the radial distributions of the azimuthal component of the transverse momentum density in circularly polarized vortex beams using a Shack-Hartmann wavefront sensor. This is a direct purely optical method to measure the transverse components of the momentum density, without involving optomechanical effects or quantum weak measurements. It was previously applied to linearly polarized vortex beams [46], where there is no difference between the canonical and Poynting momentum densities, so that the results were associated with the Poynting vector. Here we show that such a technique in fact measures the canonical momentum density. We used LG beams with $\ell = \pm 1$ and $\sigma = \pm 1$ in the focal plane $z = 0$ where the radial component of the momentum density vanishes.
The experimental setup is shown in Fig. 2. The LG beams, Eq. (5), with $\ell = 0, \pm 1$ were produced using phase masks on a phase-only spatial light modulator (Holoeye Pluto 2.0 SLM), as shown in the inset A of Fig. 2. A continuous wave laser with the wavelength $\lambda = 780\;{\rm nm} $ (Thorlabs CLD1015) was expanded and reflected off the SLM. In order to simultaneously modulate phase and amplitude using a phase-only SLM, we adopted the technique described in [61], such that the desired field was obtained after filtering the first diffraction order [70]. Polarizer P1 was used to set the linear polarization of the incident beam. Then, the light reflected from the SLM was circularly polarized using a quarter wave plate (QWP) oriented at 45° with respect to P1. Polarizer P2 was used to determine the orientations of the fast and slow axes of the QWP, and it was removed after that. Also, to determine the circular-polarization helicity $\sigma$, a $Q$-plate (QP) of the order $q = 1/2$ (Thorlabs-WPV10L-780) was employed [71] and removed afterwards. This process allowed us to create LG beams with desired circular polarization.
Next, the SLM was imaged using lenses L3 and L4 onto the microlens array (ThorLabs-MLA-150-5C-M; each lens has a pitch of 150 µm and a focal length of 5.6 mm) that focused the beam onto the CMOS camera (mvBlueFOX-200wG, pixel size 6 µm). The inset B in Fig. 2 shows the spotfield generated on the CMOS when the mask shown in the inset A was displayed on the SLM. (The spotfield is a 2D-grid-like arrangement of focal spots generated by a Shack-Hartmann wavefront sensor. This pattern is formed as each lens of the microlens array focuses the corresponding small area of the incoming wavefront onto the corresponding region on the CMOS sensor.) In accordance with the principle of the Shack-Hartmann sensor [72], a reference spotfield is generated by creating the hologram corresponding to a wide Gaussian beam on the SLM and subsequently imaging the light reflected from the SLM plane onto the microlens array. While generating the spotfield of the reference (a wide Gaussian beam), the polarizations of the light incident on and reflected from the SLM were once again set according to the aforementioned procedure.
The data analysis was performed as follows. The displacements of the centroids of each spot in the spotfield generated by the LG beam were measured with respect to the corresponding spots in the spotfield generated by the reference. The Cartesian coordinates of these displacements were used to obtain the azimuthal component of the displacement [61]. Then, the azimuthal displacement of the $i$th spot centroid, $\Delta {\varphi _i}$, is divided by the focal length ${f_m}$ of each microlens and multiplied by the corresponding intensity ${I_i}$ of the spot. The resulting quantity $P_{\varphi i}^{{\exp}} = {I_i}\Delta {\varphi _i}/{f_m}$ provides the experimentally measured azimuthal component of the momentum density, which is in agreement with the canonical momentum, Eq. (8), ${P_\varphi} \propto I\ell /kr$ (see Fig. 3). We collected several frames of the same spotfield image. For each such frame, the measured intensity and canonical momentum density are grouped into bins corresponding to segments of radii ${r_j}$. The average over all the frames of intensity ${\langle I\rangle _j}$ and canonical momentum density ${\langle P_\varphi ^{{\exp}}\rangle _j}$ are obtained for each ${r_j}$.
The results of the measurements are shown in Fig. 3. The grey symbols correspond to the measured intensity ${\langle I\rangle _j}$ as a function of the radius $r$, while the solid gray curves correspond to the theoretical dependences $|{\Psi _\ell}(r{)|^2}$, Eq. (5), optimally fitted to the experimental measurements. The red and blue curves correspond to the theoretical distributions of the azimuthal components of the canonical momentum density ${P_\varphi}(r)$ and Poynting momentum ${\Pi _\varphi}(r)$, Eqs. (8) and (9), respectively. Evidently, the measured azimuthal momentum densities (red symbols) are consistent with the canonical rather than Poynting momentum. While the azimuthal canonical momentum is independent of the polarization helicity $\sigma$ and does not change its sign across the beam radius, the azimuthal component of the Poynting vector depends on $\sigma$ and alters its sign with $r$. Our experimental results also confirm recent theoretical arguments [73] that the Shack-Hartmann wavefront sensor measures the weak value of the canonical momentum operator, i.e., the canonical momentum density.
We finally note that for the experimental data shown in Fig. 3, the statistical error in the measured azimuthal momentum density (red points) ranges from $5 \times {10^{- 3}}$ to $6 \times {10^{- 2}}$ across radius $r$. The maximum statistical error is typically observed near the vortex core as the intensities of the detected spots are low in these regions, resulting in imprecise detection of the corresponding spot centroids. It can also be observed that these data points for vortex beams ($\ell \ne 0$) exhibit some systematic errors around $r/{w_0} = 0.5$ and deviate from the theoretical red curves. These radii correspond to high intensity gradients (as seen from the gray curves), where systematic errors may appear due to cross-talks between neighboring microlenses (see supplementary section of [61] for more details).
4. PROPAGATIONAL DYNAMICS OF TRUNCATED VORTEX BEAMS
Another method to reveal the local propagational evolution of vortex beams is to consider diffraction of such beams after a suitable (e.g., knife-edge) aperture. This method has been applied in optics [48–51] and later for electron vortex beams [74,75]. However, the optical case was considered for linearly polarized beams, and the observed evolution was attributed to the Poynting vector. We argue that in fact this evolution is determined by the canonical momentum density.
We first recall the known facts. A vortex beam truncated by a knife edge or a similar aperture, usually placed in the focal plane, experiences an azimuthal rotation upon the $z$-propagation, as shown in Fig. 4. The direction of this rotation is determined by ${\rm sgn}(\ell)$ and its magnitude for LG beams is given by the Gouy phase ${\Phi _G}(z)$ [48–50,74]. The latter fact hints at the connection between the azimuthal rotation and diffraction, although similar rotations have also been observed for non-diffracting Bessel beams [48]. However, despite the general association with the Poynting vector [48,49], an accurate quantitative relation between this rotational evolution and the momentum density in vortex beams remains somewhat elusive. Even for linearly polarized beams, Berry and McDonald remarked that the rotation of the streamlines (trajectories) determined by the momentum density in LG beams “is unrelated to the Gouy phase” [52].
A. Laguerre-Gaussian Beams
To characterize the rotational evolution of truncated vortex beams within the optical-current paradigm, we employ the method used in [75] for electron vortex beams. Namely, assuming that the local velocity of the propagation in a structured optical field is given by ${\textbf v} = {c^2}{\textbf P}/W$ [19,76], we introduce the local angular velocity of the beam rotation with respect to the $z$-coordinate (rather than time):
Note that such radial distribution of the azimuthal velocity ${v_\varphi} \propto 1/r$ corresponds to “irrotational” hydrodynamical vortices, with ${({\boldsymbol \nabla} \times {\textbf v})_z} = 0$ and constant normalized angular-momentum density ${l_z} \propto r{v_\varphi}$ [77]. To characterize the global rotation of the beam, we calculate the mean value of this angular velocity akin to a quantum-mechanical expectation value:Equations (10)–(12) establish the direct relation between the azimuthal momentum density and the Gouy-phase rotation of truncated LG beams. It appears in the integral expectation values rather than in the local momentum density (i.e., individual Bohmian trajectories or weak values). Remarkably, if we use the Poynting momentum density ${\Pi _\varphi}$, Eq. (9), instead of the canonical momentum ${P_\varphi}$, Eq. (8),
Thus, the global rotation cannot discriminate between the Poynting and canonical momentum densities. To do this, one has to measure the local $r$-dependent rotations in vortex beams. In [48], following [67], this was done by involving the concept of the “radial maximum trajectory” $r = {r_{{\max}}}(z)$ based on the $z$-evolution of the radial maximum of the intensity distribution and by measuring rotations of the intensity rings with different ${r_{{\max}}}$ in beams with multiple radial maxima (e.g., Bessel beams). However, this concept of a single radial trajectory is elusive (the momentum density streamlines generate a continuum of trajectories [52], Fig. 1), and we would like to deal with the continuous momentum density field.
B. Bessel and Other Vortex Beams
Note that Expressions (8) and (10) for the $\ell$-dependent azimuthal component of the canonical momentum and the corresponding local angular velocity are universal for vortex beams with the azimuthal dependence $\propto \exp (i\ell \varphi)$. Consider, for example, paraxial non-diffracting Bessel beams [78]:
where ${J_\ell}$ is the Bessel function of the first kind, $\kappa \ll k$ is the radial wavevector component, and ${k_z} = \sqrt {{k^2} - {\kappa ^2}}$. The canonical momentum density in paraxial Bessel beams, Eq. (15), in the linear approximation in $\kappa$, equalsIn contrast, the azimuthal component of the Poynting vector differs considerably for different types of vortex beams. For paraxial Bessel beams, we have ${E_z} \propto - i\sigma (\kappa /\sqrt 2 k){J_{\ell + \sigma}}(\kappa r)\def\LDeqbreak{}{e^{i(\ell + \sigma)\varphi}}$ [66], which results in the Poynting momentum [79]:
Figure 5 shows the transverse distributions of the canonical and Poynting momentum and energy densities, as well as the streamlines of the canonical and Poynting currents for circularly polarized Bessel beams with $\ell = 1$ and $\sigma = \pm 1$. All the streamlines are spirals in this case because of the absence of the radial momentum components (i.e., the absence of diffraction); see Appendix A. Similar to the LG case, the canonical and Poynting streamlines coincide at the radial intensity maxima, i.e., for $d{J_\ell}(\kappa r)/dr = 0$. Furthermore, akin to the LG-beam case, the azimuthal component of the Poynting current has $r$-regions with positive and negative ${\Pi _\varphi}$. Thus, if the Poynting momentum density were to determine the $z$-evolution of truncated vortex beams, we could observe local negative rotations in certain $r$-ranges. Moreover, such rotations would strongly depend on the circular polarization $\sigma$. Unfortunately, the experiments [48–50] were performed for linear polarizations, but it is well known that the Kirchhoff-Fresnel integral determining diffraction of paraxial fields [81] does not depend on the polarization. Therefore, we can conclude that the local ($r$-dependent) diffraction evolution of truncated vortex beams cannot be described by the Poynting momentum density.
Although the local $r$-dependent angular velocity, Eq. (10), determined by the canonical momentum is independent of the form of the vortex beam ${\Psi _\ell}$, the resulting $z$-evolution can depend on it. Hence, we can write the local angle of the truncated vortex rotation as
where $F(z)$ is a function depending on the beam profile ${\Psi _\ell}$, which provides a $z$-scaling of the same universal $r$-dependent rotation. Figure 6(a) shows that the suitably scaled deformation of the knife-edge line described by Eq. (18) agrees very well with the experimentally observed [48] rotation of multiple truncated rings in a Bessel beam.The theoretical model of [48], based on the “radial maximum trajectory”, used Expression (18) with $F(z) = {w^2}(z){\Phi _G}(z)/2$ involving the width and Gouy phase for the LG beam, Eq. (5). Apparently, for Bessel or other beams, this function should have a different form. Since Bessel beams are non-normalizable (the integral $\int_0^\infty J_\ell ^2(\kappa r) r {\rm d}r$ diverges), let us now consider normalizable and diffracting Bessel-Gauss beams [82]:
To understand this result we calculate the integral (over all rings) rotation angle for the Bessel-Gauss beam, Eq. (19), similar to Eqs. (10)–(12):
We have found analytically that the total integral rotation for the whole range $z \in (- \infty ,\infty)$ equals
for any paraxial vortex beams of the form ${\Psi _\ell} = \psi (r,z){e^{i\ell \varphi}}$ (see Appendix B). This result is more fundamental than the Gouy-phase Eqs. (11) and (12) limited to the LG beams. The fundamental global rotation, Eq. (21), can be understood within the geometrical-optics ray picture. Indeed, the far-field $z = - \infty$ and $z = \infty$ regions can be connected by straight geometrical-optics rays [52], and projection of any straight ray onto the transverse $(x,y)$-plane covers the $\pi$ range in the azimuthal angle $\varphi$, while the direction of the azimuthal-angle evolutions is determined by ${\rm sgn}(\ell)$. Thus, the integral Relation (21) provides the connection between the local canonical-momentum current and the global geometrical-optics ray picture.5. ON THE “SUPERMOMENTUM” AND “BACKFLOW”
Here we briefly discuss the canonical and Poynting currents in the context of the “supermomentum” [6,53] and “backflow” [54,55]. Both phenomena are currently attracting significant attention [56–61] and represent examples of superoscillations [62]. Namely, if the wavefield consists of plane waves (Fourier spectrum) with some component of the wavevectors limited by the range ${k_i} \in (0,k)$, then the spatial regions where the local phase gradient (i.e., the local wavevector or normalized momentum density $\omega {P_i}/W$) is larger than $k$ or less than zero correspond to the supermomentum and backflow, respectively [19]. In other terms, these conditions correspond to the superluminal or negative local velocity ${v_i} = {c^2}{P_i}/W$: ${v_i} \gt c$ and ${v_i} \lt 0$. In the context of vortex beams, we consider these phenomena for the azimuthal momentum component, or angular momentum.
A. Supermomentum Near the Vortex Core
For the azimuthal component of the canonical momentum density, Eq. (8), the supermomentum occurs in the subwavelength vicinity of the vortex center:
However, the azimuthal component of the Poynting vector in Eqs. (9) and (17) for ${\rm sgn}(\sigma \ell) = - 1$ exhibits the same $1/r$ divergence (when divided by the intensity) as the canonical momentum and hence contradicts the general condition, Eq. (23). This means that our paraxial-approximation equations become inapplicable near the vortex center, and the second-order terms become important in this region. Indeed, near the vortex core (e.g., for $r \ll {w_0}$ in LG beams), $\partial {\Psi _\ell}/\partial r \simeq (|\ell |/r){\Psi _\ell}$, and the longitudinal field component is estimated as $|{E_z}| \sim (|\ell |/kr)|{\Psi _\ell}|$. Then, its contribution becomes of the order of the main field, $|{E_z}| \sim |{\Psi _\ell}|$, for $r \lesssim |\ell |{k^{- 1}}$, i.e., exactly in the supermomentum area, Eq. (22).
Importantly, this longitudinal-field contribution includes spin-orbit interaction effects [65,66], and all the beam properties become strongly $\sigma$-dependent in the supermomentum area. To take these effects into account, we calculate the azimuthal components of the canonical and Poynting momentum, and also the energy density, for circularly polarized Bessel beams without using paraxial approximation. In doing so, we use the exact Bessel-beam solutions described in [66]. These are helicity eigenmodes, i.e., represent a superposition of plane waves with different wave vectors and the same circular polarization (helicity) $\sigma = \pm 1$. Owing to this, the magnetic field exactly satisfies Eq. (7): ${\textbf H} = - i\sigma {\textbf E}$, whereas the Cartesian components of the electric field can be written as
Figure 7 shows the normalized values $c{P_\varphi}/W$ and $c{\Pi _\varphi}/W$ as functions of $r$, both exact and in the paraxial approximation, for Bessel beams with $\ell = 1,2,3$ and $\sigma = \pm 1$. One can see a drastic difference between the paraxial and exact curves for ${\rm sgn}(\sigma \ell) = - 1$. In agreement with the general Relation (23), the exact Poynting vector never corresponds to supermomentum. Furthermore, the exact canonical momentum does not produce supermomentum in the case of ${\rm sgn}(\sigma \ell) = - 1$ and $|\ell | \le 2$. Under these conditions $W(0) \ne 0$, and, hence, the azimuthal supermomentum can be achieved only in the case when the energy density (including all the field components) exactly vanishes in the vortex center. One can say that it is the vanishing energy density, rather than the anomalous momentum density, that provides for the supermomentum.
B. Azimuthal Backflow
Although, the backflow is usually considered for Cartesian components of the wave current, recently it was shown that it can be equally considered for the local azimuthal phase gradient [61]. In this case, one has to consider the expansion of the wavefield in the orbital-angular-momentum (vortex) eigenmodes instead of plane waves. In particular, if the field represents a superposition of the modes with $\ell \gt 0$, then the negative local azimuthal phase gradient corresponds to the azimuthal backflow.
In this context, we note that the azimuthal component of the Poynting vector in a circularly polarized vortex beam with $\ell \gt 0$ and $\sigma = - 1$, Eqs. (9) and (17), becomes negative in certain $r$-regions; see Figs. 1 and 5. However, one should not interpret it as an azimuthal backflow [60]. Such vortex beams are orbital-angular-momentum eigenmodes with fixed eigenvalues $\ell$ (in the paraxial regime one can consider the orbital and spin parts of the angular momentum independently [12,13,65]), while superoscillations are interference phenomena that appear only in superpositions of several momentum or angular-momentum eigenmodes [62]. Moreover, as we have shown, the Poynting momentum density does not correspond to the local phase gradient, or wavevector, or propagation direction of light, and the actual optical current corresponds to the canonical momentum density. Its azimuthal component is always positive in a vortex beam with fixed $\ell \gt 0$, and the interference of at least two vortex beams is needed to achieve the true azimuthal backflow [61].
6. CONCLUSIONS
In this work, we aimed to overview the roles of the canonical and Poynting momentum densities (currents) in the propagation and diffraction of structured monochromatic optical beams. In doing so, we considered vortex beams of different kinds. While the Poynting vector is determined by the well-known vector product of the electric and magnetic fields, the canonical momentum density is associated with the wave phase gradient (local wavevector) averaged over all the field components and multiplied by the wave intensity. In this manner, the canonical momentum provides a direct link to the quantum-mechanical de Broglie momentum and probability current. It is well established in other works that local optomechanical momentum effects (i.e., optical forces on small particles or momentum exchange with atoms) are determined by the canonical momentum density rather than the Poynting vector [8,13,14,17,18]. Furthermore, quantum weak measurements of the local momentum also result in the canonical momentum density [6,19,22]. In this work, we have considered the canonical-Poynting dilemma in the context of purely optical propagation and diffraction phenomena.
In the paraxial regime, assuming uniformly polarized light, the canonical current is independent of the polarization, while the Poynting vector depends on the degree of circular polarization, i.e., spin. Previous optical approaches to the optical currents in the propagation and diffraction of structured light [46,48,50,52,67] associated it with the Poynting vector. However, these works dealt with linearly polarized (zero-spin) beams where the canonical and Poynting currents coincide with each other. We have revisited these approaches using circularly polarized beams and found that it is the canonical momentum density rather than the Poynting vector that should be associated with the observable phenomena in the propagation of light.
First, we have performed experimental measurements of the transverse momentum density in vortex LG beams via a Shack-Hartmann wavefront sensor. This yielded polarization-independent results consistent with the canonical current rather than the Poynting vector. Second, we have analyzed diffraction of different kinds of vortex beams at a knife-edge aperture placed in the focal plane. The truncated beams experience a $z$-dependent rotation that depends on the sign of the vortex charge $\ell$ and can be associated with the Gouy phase. We have shown that this $z$-dependent rotation can be explained by the optical current properly integrated over the transverse cross-section of the beam. Notably, this integral current is the same in the canonical and Poynting versions; only the local $r$-dependent rotation can discriminate between these. However, since the paraxial diffraction is polarization-independent, we can conclude that it must be locally defined by the canonical rather than Poynting momentum density. Remarkably, we have derived a theoretical relation for the total integral rotation over the whole $z \in (- \infty ,\infty)$ range in a circularly symmetric vortex beam with an arbitrary $r$-profile. The angle of this rotation always equals $\pi {\rm sgn}(\ell)$, which provides a global relation between the optical currents, Gouy phase, and the geometrical-optics ray picture. While optical currents describe local propagation properties of a structured wavefield, geometrical optics rays connect far-field regions in the field evolution.
Finally, we have discussed the “supermomentum” [6,19,53,56,57] and “backflow” [19,54,55,58–61] effects in the canonical and Poynting currents in vortex beams. The former phenomenon means anomalously high momentum density normalized by the energy density, and it appears in a subwavelength area near the vortex core. We have shown that it is only possible for the canonical momentum density and strongly depends on the nonparaxial spin-dependent corrections. In this regime, the energy and canonical momentum densities become polarization-dependent. The latter phenomenon means anomalously negative local optical current, and it should be also associated with the canonical momentum density in interfering modes. While the Poynting vector can have negative azimuthal components in vortex beams with $\ell \gt 0$, it should not be associated with the backflow, because the actual propagation is determined by the canonical momentum density. Interference of at least two vortex modes with different $\ell$ is required for the proper canonical azimuthal backflow [61].
We hope that this work sheds light on important peculiarities of optical propagation and the relevant physical characteristics for their description. The canonical momentum density provides a useful theoretical tool for understanding the free-space evolution of monochromatic light fields. It would be important to further explore and extend this framework beyond its current limitations: e.g., for polychromatic optical fields and for structured light in anisotropic media.
APPENDIX A: STREAMLINES OF THE CANONICAL AND POYNTING CURRENTS IN LG AND BESSEL BEAMS
The streamlines $(r(z),\varphi (z))$ of a vector field ${\textbf V}({\textbf r})$ are solutions of differential equations
In turn, for the Poynting momentum density Eq. (9), the azimuthal angle equations are modified as
For Bessel beams, Eqs. (15)–(17), the radial evolution is trivial: $dr/dz = 0$, $r(z) = {r_0}$, while the azimuthal equations for the canonical and Poynting streamlines yield
APPENDIX B: INTEGRAL ROTATION ANGLE FOR ARBITRARY VORTEX BEAMS
Here we calculate the integral rotation angle over the whole range $z \in (- \infty ,\infty)$ for an arbitrary paraxial vortex beam ${\Psi _\ell} = \psi (r,z){e^{i\ell \varphi}}$:
To evaluate this expression, we represent the real-space wave function ${\Psi _\ell}$ by the plane-wave Fourier integral in cylindrical coordinates:
Thus, the denominator in Eq. (B.1) is $z$-independent, and the $z$-integration can be performed in the numerator. Substituting there the Fourier representation Eq. (B.3), akin to Eq. (B.4), and using the relation
Finally, substituting Eqs. (B.4) and (B.5) into Eq. (B.1), we arrive at Eq. (21):
Funding
Ministerstwo Edukacji i Nauki (contract MEiN/2023/DIR/3797); Horizon 2020 Framework Programme (the ENSEMBLE3 project MAB/2020/14); Narodowe Centrum Nauki (grant 2022/47/B/ST7/03465); European Regional Development Fund (the FIRST TEAM project POIR.04.04.00-00-3004/17-00, the TEAM project TEAM/2016-3/29); Fundacja na rzecz Nauki Polskiej (the FIRST TEAM project POIR.04.04.00-00-3004/17-00, the TEAM project TEAM/2016-3/29).
Acknowledgment
We are grateful to Miguel A. Alonso for helpful discussions and to Filippo Cardano for his support with the Q-plates for our experiments.
Disclosures
The authors declare no conflicts of interest.
Data availability
All data needed to evaluate the conclusions in this paper are present in the paper.
REFERENCES
1. H. Rubinsztein-Dunlop, A. Forbes, M. V. Berry, et al., “Roadmap on structured light,” J. Opt. 19, 013001 (2016). [CrossRef]
2. K. Y. Bliokh, E. Karimi, M. J. Padgett, et al., “Roadmap on structured waves,” J. Opt. 25, 103001 (2023). [CrossRef]
3. I. Brevik, “Experiments in phenomenological electrodynamics and the electromagnetic energy-momentum tensor,” Phys. Rep. 52, 133 (1979). [CrossRef]
4. R. N. C. Pfeifer, T. A. Nieminen, N. R. Heckenberg, et al., “Colloquium: momentum of an electromagnetic wave in dielectric media,” Rev. Mod. Phys. 79, 1197 (2007). [CrossRef]
5. S. M. Barnett and R. Loudon, “The enigma of optical momentum in a medium,” Philos. Trans. R. Soc. A 368, 927 (2010). [CrossRef]
6. M. V. Berry, “Optical currents,” J. Opt. A 11, 094001 (2009). [CrossRef]
7. A. Bekshaev, K. Y. Bliokh, and M. Soskin, “Internal flows and energy circulation in light beams,” J. Opt. 13, 053001 (2011). [CrossRef]
8. K. Y. Bliokh, A. Y. Bekshaev, and F. Nori, “Extraordinary momentum and spin in evanescent waves,” Nat. Commun. 5, 3300 (2014). [CrossRef]
9. A. Boivin, J. Dow, and E. Wolf, “Energy flow in the neighborhood of the focus of a coherent beam,” J. Opt. Soc. Am. 57, 1171–1175 (1967). [CrossRef]
10. W. Braunbek and G. Laukien, “Features of refraction by a semi-plane,” Optik 9, 174–179 (1952).
11. M. V. Berry, “Lateral and transverse shifts in reflected dipole radiation,” Proc. R. Soc. A 467, 2500–2519 (2011). [CrossRef]
12. L. Allen, M. J. Padgett, and M. Babiker, “The orbital angular momentum of light,” Prog. Opt. 39, 291–372 (1999). [CrossRef]
13. K. Y. Bliokh and F. Nori, “Transverse and longitudinal angular momenta of light,” Phys. Rep. 592, 1–38 (2015). [CrossRef]
14. E. Leader, “The photon angular momentum controversy: resolution of a conflict between laser optics and particle physics,” Phys. Lett. B 756, 303–308 (2016). [CrossRef]
15. K. Y. Bliokh, A. Y. Bekshaev, and F. Nori, “Optical momentum, spin, and angular momentum in dispersive media,” Phys. Rev. Lett. 119, 073901 (2017). [CrossRef]
16. K. Y. Bliokh, A. Y. Bekshaev, and F. Nori, “Dual electromagnetism: helicity, spin, momentum and angular momentum,” New J. Phys. 15, 033026 (2013). [CrossRef]
17. A. Ashkin and J. P. Gordon, “Stability of radiation-pressure particle traps: an optical earnshaw theorem,” Opt. Lett. 8, 511–513 (1983). [CrossRef]
18. M. I. Marqués, “Beam configuration proposal to verify that scattering forces come from the orbital part of the Poynting vector,” Opt. Lett. 39, 5122–5125 (2014). [CrossRef]
19. K. Y. Bliokh, A. Y. Bekshaev, A. G. Kofman, et al., “Photon trajectories, anomalous velocities and weak measurements: a classical interpretation,” New J. Phys. 15, 073022 (2013). [CrossRef]
20. E. Madelung, “Quantentheorie in hydrodynamische form,” Z. Phys. 40, 322–326 (1927). [CrossRef]
21. D. Bohm, B. J. Hiley, and P. N. Kaloyerou, “An ontological basis for the quantum theory,” Phys. Rep. 144, 321–375 (1987). [CrossRef]
22. S. Kocsis, B. Braverman, S. Ravets, et al., “Observing the average trajectories of single photons in a two-slit interferometer,” Science 332, 1170–1173 (2011). [CrossRef]
23. F. I. Fedorov, “To the theory of total reflection,” Dokl. Akad. Nauk. SSSR 105, 465–468 (1955).
24. C. Imbert, “Calculation and experimental proof of the transverse shift induced by total internal reflection of a circularly polarized light beam,” Phys. Rev. D 5, 787–796 (1972). [CrossRef]
25. A. Y. Bekshaev, K. Y. Bliokh, and F. Nori, “Transverse spin and momentum in two-wave interference,” Phys. Rev. X 5, 011039 (2015). [CrossRef]
26. M. Antognozzi, C. R. Bermingham, R. L. Harniman, et al., “Direct measurements of the extraordinary optical momentum and transverse spin-dependent force using a nano-cantilever,” Nat. Phys. 12, 731 (2016). [CrossRef]
27. W. Gordon, “Der strom der diracschen elektronentheorie,” Z. Phys. 50, 630–632 (1928). [CrossRef]
28. E. Leader and C. Lorcé, “The angular momentum controversy: what’s it all about and does it matter?” Phys. Rep. 541, 163–248 (2014). [CrossRef]
29. C. Shi, R. Zhao, Y. Long, et al., “Observation of acoustic spin,” Natl. Sci. Rev. 6, 707–712 (2019). [CrossRef]
30. L. Burns, K. Y. Bliokh, F. Nori, et al., “Acoustic versus electromagnetic field theory: scalar, vector, spinor representations and the emergence of acoustic spin,” New J. Phys. 22, 053050 (2020). [CrossRef]
31. Y. Long, J. Ren, and H. Chen, “Intrinsic spin of elastic waves,” Proc. Natl. Acad. Sci. USA 115, 9951–9955 (2018). [CrossRef]
32. G. J. Chaplain, J. M. D. Ponti, and R. V. Craster, “Elastic orbital angular momentum,” Phys. Rev. Lett. 128, 064301 (2022). [CrossRef]
33. K. Y. Bliokh, “Elastic spin and orbital angular momenta,” Phys. Rev. Lett. 129, 204303 (2022). [CrossRef]
34. K. Y. Bliokh, H. Punzmann, H. Xia, et al., “Field theory spin and momentum in water waves,” Sci. Adv. 8, eabm1295 (2022). [CrossRef]
35. N. Umow, “Ein theorem über die wechselwirkungen in endlichen entfernungen,” Z. Math. Phys. 19, 97–114 (1874).
36. J. H. Poynting, “On the transfer of energy in the electromagnetic field,” Philos. Trans. R. Soc. London 175, 343–361 (1884). [CrossRef]
37. M. E. McIntyre, “On the ‘wave momentum’ myth,” J. Fluid Mech. 106, 331–347 (1981). [CrossRef]
38. R. Peierls, Surprises in Theoretical Physics (Princeton University, 1979).
39. R. Peierls, More Surprises in Theoretical Physics (Princeton University, 1991).
40. T. S. van den Bremer and Ø. Breivik, “Stokes drift,” Philos. Trans. R. Soc. A 376, 20170104 (2017). [CrossRef]
41. K. Y. Bliokh, Y. P. Bliokh, and F. Nori, “Ponderomotive forces, Stokes drift, and momentum in acoustic and electromagnetic waves,” Phys. Rev. A 106, L021503 (2022). [CrossRef]
42. A. Bekshaev, M. Soskin, and M. Vasnetsov, Paraxial Light Beams with Angular Momentum (Nova Science, 2008).
43. O. V. Angelsky, A. Y. Bekshaev, P. P. Maksimyak, et al., “Orbital rotation without orbital angular momentum: mechanical action of the spin part of the internal energy flow in light beams,” Opt. Express 20, 3563–3571 (2012). [CrossRef]
44. A. Yevick, D. B. Ruffner, and D. G. Grier, “Tractor beams in the Rayleigh limit,” Phys. Rev. A 93, 043807 (2016). [CrossRef]
45. A. Hayat, J. P. B. Mueller, and F. Capasso, “Lateral chirality-sorting optical forces,” Proc. Natl. Acad. Sci. USA 112, 13190–13194 (2015). [CrossRef]
46. J. Leach, S. Keen, M. J. Padgett, et al., “Direct measurement of the skew angle of the Poynting vector in a helically phased beam,” Opt. Express 14, 11919 (2006). [CrossRef]
47. K. Murphy, D. Burke, N. Devaney, et al., “Experimental detection of optical vortices with a Shack-Hartmann wavefront sensor,” Opt. Express 18, 15448 (2010). [CrossRef]
48. J. Arlt, “Handedness and azimuthal energy flow of optical vortex beams,” J. Mod. Opt. 50, 1573 (2003). [CrossRef]
49. J. Hamazaki, Y. Mineta, K. Oka, et al., “Direct observation of Gouy phase shift in a propagating optical vortex,” Opt. Express 14, 8382 (2006). [CrossRef]
50. H. X. Cui, X. L. Wang, B. Gu, et al., “Angular diffraction of an optical vortex induced by the Gouy phase,” J. Opt. 14, 055707 (2012). [CrossRef]
51. A. Y. Bekshaev, K. A. Mohammed, and I. A. Kurka, “Transverse energy circulation and the edge diffraction of an optical vortex beam,” Appl. Opt. 53, B27–B37 (2014). [CrossRef]
52. M. V. Berry and K. T. McDonald, “Exact and geometrical optics energy trajectories in twisted beams,” J. Opt. A 10, 035005 (2008). [CrossRef]
53. S. M. Barnett and M. V. Berry, “Superweak momentum transfer near optical vortices,” J. Opt. 15, 125701 (2013). [CrossRef]
54. A. J. Bracken and G. F. Melloy, “Probability backflow and a new dimensionless quantum number,” J. Phys. A: Math. Gen. 27, 2197–2211 (1994). [CrossRef]
55. M. V. Berry, “Quantum backflow, negative kinetic energy, and optical retro-propagation,” J. Phys. A: Math. Theor. 43, 415302 (2010). [CrossRef]
56. A. Afanasev, C. E. Carlson, and A. Mukherjee, “Recoil momentum effects in quantum processes induced by twisted photons,” Phys. Rev. Res. 3, 023097 (2021). [CrossRef]
57. I. P. Ivanov, B. Liu, and P. Zhang, “Observability of the superkick effect within a quantum-field-theoretical approach,” Phys. Rev. A 105, 013522 (2022). [CrossRef]
58. Y. Eliezer, T. Zacharias, and A. Bahabad, “Observation of optical backflow,” Optica 7, 72–76 (2020). [CrossRef]
59. A. Daniel, B. Ghosh, B. Gorzkowski, et al., “Demonstrating backflow in classical two beams’ interference,” New J. Phys. 24, 123011 (2022). [CrossRef]
60. I. Bialynicki-Birula, Z. Bialynicka-Birula, and S. Augustynowicz, “Backflow in relativistic wave equations,” J. Phys. A: Math. Theor. 55, 255702 (2022). [CrossRef]
61. B. Ghosh, A. Daniel, B. Gorzkowski, et al., “Azimuthal backflow in light carrying orbital angular momentum,” Optica 10, 1217–1222 (2023). [CrossRef]
62. M. Berry, N. Zheludev, Y. Aharonov, et al., “Roadmap on superoscillations,” J. Opt. 21, 053002 (2019). [CrossRef]
63. M. Born and E. Wolf, Principles of Optics (Cambridge University, 1999).
64. M. V. Berry, “Five momenta,” Eur. J. Phys. 34, 1337–1348 (2013). [CrossRef]
65. K. Y. Bliokh, F. J. Rodríguez-Fortuño, F. Nori, et al., “Spin-orbit interactions of light,” Nat. Photonics 9, 796–808 (2015). [CrossRef]
66. K. Y. Bliokh, M. A. Alonso, E. A. Ostrovskaya, et al., “Angular momenta and spin-orbit interaction of nonparaxial light in free space,” Phys. Rev. A 82, 063825 (2010). [CrossRef]
67. M. Padgett and L. Allen, “The Poynting vector in Laguerre-Gaussian laser modes,” Opt. Commun. 121, 36 (1995). [CrossRef]
68. A. Y. Bekshaev, M. S. Soskin, and M. V. Vasnetsov, “Centrifugal transformation of the transverse structure of freely propagating paraxial light beams,” Opt. Lett. 31, 694–696 (2006). [CrossRef]
69. E. G. Abramochkin and V. G. Volostnikov, “Spiral light beams,”Phys.–Usp. 47, 1177–1203 (2004). [CrossRef]
70. E. Bolduc, N. Bent, E. Santamato, et al., “Exact solution to simultaneous intensity and phase encryption with a single phase-only hologram,” Opt. Lett. 38, 3546–3549 (2013). [CrossRef]
71. A. Rubano, F. Cardano, B. Piccirillo, et al., “Q-plate technology: a progress review,” J. Opt. Soc. Am. B 36, D70–D87 (2019). [CrossRef]
72. F. Kong, M. C. Polo, and A. Lambert, “Centroid estimation for a Shack–Hartmann wavefront sensor based on stream processing,” Appl. Opt. 56, 6466–6475 (2017). [CrossRef]
73. Y. Zheng, M. Yang, Z.-H. Liu, et al., “Detecting momentum weak value: Shack–Hartmann versus a weak measurement wavefront sensor,” Opt. Lett. 46, 5352–5355 (2021). [CrossRef]
74. G. Guzzinati, P. Schattschneider, K. Y. Bliokh, et al., “Observation of the Larmor and Gouy rotations with electron vortex beams,” Phys. Rev. Lett. 110, 093601 (2013). [CrossRef]
75. P. Schattschneider, T. Schachinger, M. Stöger-Pollach, et al., “Imaging the dynamics of free-electron Landau states,” Nat. Commun. 5, 4586 (2014). [CrossRef]
76. J. Lekner, “Phase and transport velocities in particle and electromagnetic beams,” J. Opt. A 4, 491 (2002). [CrossRef]
77. A. Y. Bekshaev, M. S. Soskin, and M. V. Vasnetsov, “An optical vortex as a rotating body: mechanical features of a singular light beam,” J. Opt. A 6, S170–S174 (2004). [CrossRef]
78. D. McGloin and K. Dholakia, “Bessel beams: diffraction in a new light,” Contemp. Phys. 46, 15–28 (2005). [CrossRef]
79. K. Volke-Sepulveda, V. Garcés-Chávez, S. Chávez-Cerda, et al., “Orbital angular momentum of a high-order Bessel light beam,” J. Opt. B: Quantum Semiclass. Opt. 4, S82–S89 (2002). [CrossRef]
80. K. Y. Bliokh, M. R. Dennis, and F. Nori, “Relativistic electron vortex beams: angular momentum and spin-orbit interaction,” Phys. Rev. Lett. 107, 174802 (2011). [CrossRef]
81. Y. A. Anan’ev, Laser Resonators and the Beam Divergence Problem (CRC Press, 1992).
82. R. Borghi, M. Santarsiero, and M. A. Porras, “Nonparaxial Bessel-Gauss beams,” J. Opt. Soc. Am. A 18, 1618–1626 (2001). [CrossRef]
83. S. Huard and C. Imbert, “Measurement of exchanged momentum during interaction between surface-wave and moving atom,” Opt. Commun. 24, 185–189 (1978). [CrossRef]
84. T. Matsudo, Y. Takahara, H. Hori, et al., “Pseudomomentum transfer from evanescent waves to atoms measured by saturated absorption spectroscopy,” Opt. Commun. 145, 64–68 (1998). [CrossRef]