Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Transverse spin and orbital angular momenta in rotating-wave structured light

Zhen-Lai Wang wangzhenlai@hbpu.edu.cn Center for Fundamental Physics and School of Mathematics and Physics, Hubei Polytechnic University, Huangshi 435003, China
(June 19, 2024)
Abstract

Within monochromatic optical fields, we demonstrate the rotating-wave structured light with wave vortex carrying an intrinsic transverse orbital angular momentum orthogonal to the propagation direction of light. Remarkably, we find that such a rotating-wave structured light reveals highly nontrivial features of transverse spin and orbital angular momentum densities. The normalized total angular momentum density is conserved universally, suggesting a mutual conversion of the intrinsic transverse spin and orbital angular momentum in free space. Despite such mutual conversion at local level, the integral intrinsic orbital angular momentum can be well defined with the topological charge of the vortex per photon. Moreover, the orientation of transverse spin density is governed by the direction of the Poynting momentum density, manifesting a spin-momentum-locking effect in free space.

I Introduction

Optical angular momentum manifests itself in two distinct forms [1]. One is spin angular momentum (SAM), which originates from the intrinsic polarization degrees of freedom [2]; the other is orbital angular momentum (OAM), which arises from the spatial phase degrees of freedom [3]. The former is purely intrinsic, or, more specifically, independent of the coordinate origin. The latter is allowed to make a clear distinction between the intrinsic and extrinsic ingredients [4]. SAM and OAM can be separately observable by the distinct effects upon test particles [5]. For the conventional light beams, both SAM and OAM are essentially longitudinal since they are aligned with the direction of propagation.

Nevertheless, there has been a growing interest in the transverse angular momentum of light [6]. Unlike the conventional and longitudinal angular momentum of light, the optical transverse angular momentum is perpendicular to the propagation direction of light. The transverse SAM (t-SAM) has been generated by evanescent waves [7, 8, 9], interference fields [10, 11] and focused beams [12, 13, 14]. The appearance of t-SAM requires the transverse inhomogeneity of the field intensity [6] and can be closely linked to the curl of the energy flow density [15, 16]. The t-SAM plays crucial roles in the intrinsic spin-momentum-locking effect [17, 18] and the intrinsic quantum spin Hall effect of light [19].

The extrinsic transverse OAM (t-OAM) is the most straightforward and trivial form of transverse angular momentum, which is analogous to the classical angular momentum [6]. Despite its trivial derivation, the extrinsic t-OAM can make remarkable physical effects including the spin-Hall and orbital-Hall effects of light [20, 21, 22, 23, 24, 25, 26, 27, 28, 29]. Recently, the intrinsic t-OAM is studied intensively in theory [30, 31, 32, 33] and in experiment [34, 35, 36, 37] with a new type of localized optical pulses characterized by spatiotemporal vortex structure. However, unlike the typical light beams with angular momentum, such pulses are fundamentally polychromatic as opposed to monochromatic.

Thus the intrinsic t-OAM stands for the last missing piece of transverse angular momentum within monochromatic optical fields. Recently, the intrinsic t-OAM of monochromatic optical fields is reported to be generated through the interference of two parallel and counter-propagating linearly-polarized focused beams, and to be observed directly by using an optomechanical sensor [38].

The motivation for this work is to find the intrinsic t-OAM in the monochromatic optical fields. In this work, we propose that the rotating-wave monochromatic structured light can carry the intrinsic t-OAM. In Sec. II, the rotating wave is introduced concisely, which is a monochromatic wave with wave vortex, propagating along the azimuthal direction and consisting of infinite plane waves. In Sec. III, we construct two modes of the rotating-wave structured light in the vector formulation of electromagnetic field. In Sec. IV, we show that the rotating-wave monochromatic structured light can possess an integral intrinsic t-OAM of well-defined value per photon and exhibits locally unusual t-SAM and t-OAM properties. Finally, the conclusion are presented.

II Rotating wave

Refer to caption
Refer to caption

(a)                               (b)

Figure 1: (a) The wave-vectors of the plane waves point radially outward from the origin and rotates on the kysubscript𝑘𝑦k_{y}italic_k start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT-axis to form a circle in the plane (kz,kx)subscript𝑘𝑧subscript𝑘𝑥(k_{z},k_{x})( italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) plane of 𝒌𝒌\bm{k}bold_italic_k space. (b) Surface plot of the rotating wave, calculated via the real part of the complex wavefunction in Eq. (1) with l=3𝑙3l=3italic_l = 3. The arrow denotes the wave’s direction of propagation.

Throughout this paper, let us restrict ourselves to the monochromatic wave field propagating in free space. The typical and simple wave field with OAM is Bessel beam. For simplicity, we can consider specifically a scalar Bessel beam propagating in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane. As illustrated in Fig. 1(a), such beam can be generated by superposing plane waves with the same frequency ω𝜔\omegaitalic_ω, wave vectors 𝒌𝒌\bm{k}bold_italic_k pointing radially outward from the origin and forming a circle in the (kz,kx)subscript𝑘𝑧subscript𝑘𝑥(k_{z},k_{x})( italic_k start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT , italic_k start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) plane, and with the corresponding azimuthal phase difference lϕ𝑙italic-ϕl\phiitalic_l italic_ϕ. Here l=0,±1,±2,𝑙0plus-or-minus1plus-or-minus2l=0,\pm 1,\pm 2,...italic_l = 0 , ± 1 , ± 2 , … is the topological charge of the vortex and ϕitalic-ϕ\phiitalic_ϕ is the azimuthal angle in 𝒌𝒌\bm{k}bold_italic_k space. Thus, the resulting scalar wave function can be written as the integral over ϕitalic-ϕ\phiitalic_ϕ from 00 to 2π2𝜋2\pi2 italic_π:

ψ(𝒓,t)𝜓𝒓𝑡\displaystyle\psi(\bm{r},t)italic_ψ ( bold_italic_r , italic_t ) =02πeikcosϕz+iksinϕx+ilϕiωtdϕabsentsuperscriptsubscript02𝜋superscript𝑒𝑖𝑘italic-ϕ𝑧𝑖𝑘italic-ϕ𝑥𝑖𝑙italic-ϕ𝑖𝜔𝑡differential-ditalic-ϕ\displaystyle=\int_{0}^{2\pi}e^{ik\cos\phi z+ik\sin\phi x+il\phi-i\omega t}{% \rm d}\phi= ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 italic_π end_POSTSUPERSCRIPT italic_e start_POSTSUPERSCRIPT italic_i italic_k roman_cos italic_ϕ italic_z + italic_i italic_k roman_sin italic_ϕ italic_x + italic_i italic_l italic_ϕ - italic_i italic_ω italic_t end_POSTSUPERSCRIPT roman_d italic_ϕ
=2πilJl(kρ)eilφiωt.absent2𝜋superscript𝑖𝑙subscript𝐽𝑙𝑘𝜌superscript𝑒𝑖𝑙𝜑𝑖𝜔𝑡\displaystyle=2\pi{i}^{l}J_{l}(k\rho)e^{il\varphi-i\omega t}.= 2 italic_π italic_i start_POSTSUPERSCRIPT italic_l end_POSTSUPERSCRIPT italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ( italic_k italic_ρ ) italic_e start_POSTSUPERSCRIPT italic_i italic_l italic_φ - italic_i italic_ω italic_t end_POSTSUPERSCRIPT . (1)

Here Jl(kρ)subscript𝐽𝑙𝑘𝜌J_{l}(k\rho)italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ( italic_k italic_ρ ) denotes the l𝑙litalic_lth-order Bessel function of the first kind, k=ω/c𝑘𝜔𝑐k=\omega/citalic_k = italic_ω / italic_c is the wave number and c𝑐citalic_c is the phase velocity of wave, and (ρ,φ)𝜌𝜑(\rho,\varphi)( italic_ρ , italic_φ ) indicate the polar coordinates in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane associated with z=ρcosφ𝑧𝜌𝜑z=\rho\cos\varphiitalic_z = italic_ρ roman_cos italic_φ and x=ρsinφ𝑥𝜌𝜑x=\rho\sin\varphiitalic_x = italic_ρ roman_sin italic_φ.

Note that Eq. (1) stands for a pure rotating wave [39]. The term eilφiωtsuperscript𝑒𝑖𝑙𝜑𝑖𝜔𝑡e^{il\varphi-i\omega t}italic_e start_POSTSUPERSCRIPT italic_i italic_l italic_φ - italic_i italic_ω italic_t end_POSTSUPERSCRIPT bears resemblance to the typical term eikxiωtsuperscript𝑒𝑖𝑘𝑥𝑖𝜔𝑡e^{ikx-i\omega t}italic_e start_POSTSUPERSCRIPT italic_i italic_k italic_x - italic_i italic_ω italic_t end_POSTSUPERSCRIPT in a traveling wave, enabling the wave to preserve its waveform while propagating along the azimuthal direction in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane. As an illustration, Fig. 1(b) shows the rotating wave for the mode l=3𝑙3l=3italic_l = 3. Unlike traveling waves, rotating waves are confined to specific modes due to their circular propagation. For mechanical waves, such as surface waves in water [39] and acoustic waves [40, 41, 42], rotating waves offer a clear understanding of angular momentum, a concept commonly encountered in quantum waves or electromagnetic waves [43, 44]. More interestingly, rotating waves can provide us with new insights into the intrinsic t-SAM and t-OAM of light, as we will discuss in the following sections.

III Rotating-wave structured light

We now consider the rotating-wave structured light. Light is a form of electromagnetic waves which manifest as vector fields, so there is a need for constructing vector beams of rotating wave. For simplicity and more concentration on the intrinsic t-OAM, we only consider the rotating-wave beam consisting entirely of linearly polarized plane waves. Due to the transversality of electromagnetic waves, the electric and magnetic fields 𝓔𝓔\bm{\mathcal{E}}bold_caligraphic_E, 𝓑𝓑\bm{\mathcal{B}}bold_caligraphic_B and the wave vector 𝒌𝒌\bm{k}bold_italic_k of each plane wave in the beam spectrum constitute a right-handed orthogonal triad. As a consequence, two basic modes of polarization should be considered in the problem. One mode is that 𝓑𝓑\bm{\mathcal{B}}bold_caligraphic_B is set along the positive y𝑦yitalic_y-direction and so 𝓔𝓔\bm{\mathcal{E}}bold_caligraphic_E is put in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane; see Fig. 2(a). Conversely, the other one is that 𝓔𝓔\bm{\mathcal{E}}bold_caligraphic_E is fixed along the negative y𝑦yitalic_y-direction and so 𝓑𝓑\bm{\mathcal{B}}bold_caligraphic_B is placed in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane; see Fig. 2(b).

Refer to caption
Refer to caption

(a)                                      (b)

Figure 2: Schematic diagrams of the two basic modes of polarization. (a) 𝓑𝓑\bm{\mathcal{B}}bold_caligraphic_B is polarized in the positive y𝑦yitalic_y-direction and then 𝓔𝓔\bm{\mathcal{E}}bold_caligraphic_E is located in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane. (b) 𝓔𝓔\bm{\mathcal{E}}bold_caligraphic_E points in the negative y𝑦yitalic_y direction and then 𝓑𝓑\bm{\mathcal{B}}bold_caligraphic_B lies in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane. Dots and crosses are used to represent a field coming out of or going into the plane of the paper, respectively.

With the help of Maxwell’s equations, the resulting electric and magnetic fields of the first mode can be given by

𝑩=ψ(𝒓,t)𝒆y,𝑬=ik×𝑩.formulae-sequence𝑩𝜓𝒓𝑡subscript𝒆𝑦𝑬𝑖𝑘bold-∇𝑩\bm{B}=\psi(\bm{r},t)\bm{e}_{y},~{}~{}~{}\bm{E}=\frac{i}{k}\bm{\nabla}\times% \bm{B}.bold_italic_B = italic_ψ ( bold_italic_r , italic_t ) bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , bold_italic_E = divide start_ARG italic_i end_ARG start_ARG italic_k end_ARG bold_∇ × bold_italic_B . (2)

Here 𝒆ysubscript𝒆𝑦\bm{e}_{y}bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT is the unit vector in the y𝑦yitalic_y-direction. For the convenience of future use and with a slight computation cost, the explicit components of the resulting electric field can be obtained

Exsubscript𝐸𝑥\displaystyle E_{x}italic_E start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT =il+1π[Jl1ei(l1)φJl+1ei(l+1)φ]eiωt,absentsuperscript𝑖𝑙1𝜋delimited-[]subscript𝐽𝑙1superscript𝑒𝑖𝑙1𝜑subscript𝐽𝑙1superscript𝑒𝑖𝑙1𝜑superscript𝑒𝑖𝜔𝑡\displaystyle=-i^{l+1}\pi\big{[}J_{l-1}e^{i(l-1)\varphi}-J_{l+1}e^{i(l+1)% \varphi}\big{]}e^{-i\omega t},= - italic_i start_POSTSUPERSCRIPT italic_l + 1 end_POSTSUPERSCRIPT italic_π [ italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_l - 1 ) italic_φ end_POSTSUPERSCRIPT - italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_l + 1 ) italic_φ end_POSTSUPERSCRIPT ] italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT , (3)
Ezsubscript𝐸𝑧\displaystyle E_{z}italic_E start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT =ilπ[Jl1ei(l1)φ+Jl+1ei(l+1)φ]eiωt.absentsuperscript𝑖𝑙𝜋delimited-[]subscript𝐽𝑙1superscript𝑒𝑖𝑙1𝜑subscript𝐽𝑙1superscript𝑒𝑖𝑙1𝜑superscript𝑒𝑖𝜔𝑡\displaystyle=-i^{l}\pi\big{[}J_{l-1}e^{i(l-1)\varphi}+J_{l+1}e^{i(l+1)\varphi% }\big{]}e^{-i\omega t}.= - italic_i start_POSTSUPERSCRIPT italic_l end_POSTSUPERSCRIPT italic_π [ italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_l - 1 ) italic_φ end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT italic_e start_POSTSUPERSCRIPT italic_i ( italic_l + 1 ) italic_φ end_POSTSUPERSCRIPT ] italic_e start_POSTSUPERSCRIPT - italic_i italic_ω italic_t end_POSTSUPERSCRIPT . (4)

For notational simplicity, from now on the argument kρ𝑘𝜌k\rhoitalic_k italic_ρ will be leave out from the Bessel functions.

In the similar way, the second mode can be represented by

𝑬=ψ(𝒓,t)𝒆y,𝑩=ik×𝑬.formulae-sequence𝑬𝜓𝒓𝑡subscript𝒆𝑦𝑩𝑖𝑘bold-∇𝑬\bm{E}=-\psi(\bm{r},t)\bm{e}_{y},~{}~{}~{}\bm{B}=-\frac{i}{k}\bm{\nabla}\times% \bm{E}.bold_italic_E = - italic_ψ ( bold_italic_r , italic_t ) bold_italic_e start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT , bold_italic_B = - divide start_ARG italic_i end_ARG start_ARG italic_k end_ARG bold_∇ × bold_italic_E . (5)

Clearly, the second mode is able to be fulfilled properly from the first one by the exchange 𝑬𝑩𝑬𝑩\bm{E}\to\bm{B}bold_italic_E → bold_italic_B and 𝑩𝑬𝑩𝑬\bm{B}\to-\bm{E}bold_italic_B → - bold_italic_E. In fact, such exchange is a typical transformation of electromagnetic duality in the free space [45, 46]. As a direct result of the transformation of electromagnetic duality, the case of the second mode can be derived from that of the first one. For this reason, we will concentrate on the first mode in the following section.

IV Spin and orbital angular momenta

In this core section, we turn now to analyze the t-SAM and t-OAM of the rotating-wave structured light. We start with the time-averaged densities of energy W𝑊Witalic_W, SAM 𝑺𝑺\bm{S}bold_italic_S and OAM 𝑳𝑳\bm{L}bold_italic_L in a monochromatic optical field with separate electric- and magnetic-field contributions

W𝑊\displaystyle Witalic_W =\displaystyle== We+Wm=gω2(|𝑬|2+|𝑩|2),superscript𝑊𝑒superscript𝑊𝑚𝑔𝜔2superscript𝑬2superscript𝑩2\displaystyle W^{e}+W^{m}=\frac{g\omega}{2}\big{(}|\bm{E}|^{2}+|\bm{B}|^{2}% \big{)},italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT + italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = divide start_ARG italic_g italic_ω end_ARG start_ARG 2 end_ARG ( | bold_italic_E | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + | bold_italic_B | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (6)
𝑺𝑺\displaystyle\bm{S}bold_italic_S =\displaystyle== 𝑺e+𝑺m=g2Im(𝑬×𝑬+𝑩×𝑩),superscript𝑺𝑒superscript𝑺𝑚𝑔2Imsuperscript𝑬𝑬superscript𝑩𝑩\displaystyle{\bm{S}}^{e}+{\bm{S}}^{m}=\frac{g}{2}{\rm Im\,}\big{(}\bm{E}^{*}% \times\bm{E}+\bm{B}^{*}\times\bm{B}\big{)},bold_italic_S start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT + bold_italic_S start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = divide start_ARG italic_g end_ARG start_ARG 2 end_ARG roman_Im ( bold_italic_E start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT × bold_italic_E + bold_italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT × bold_italic_B ) , (7)
𝑳𝑳\displaystyle\bm{L}bold_italic_L =\displaystyle== 𝑳e+𝑳m=g2Re(Ej𝓛^Ej+Bj𝓛^Bj).superscript𝑳𝑒superscript𝑳𝑚𝑔2Resubscriptsuperscript𝐸𝑗^𝓛subscript𝐸𝑗subscriptsuperscript𝐵𝑗^𝓛subscript𝐵𝑗\displaystyle{\bm{L}}^{e}+{\bm{L}}^{m}=\frac{g}{2}{\rm Re\,}\big{(}E^{*}_{j}% \hat{\bm{\mathcal{L}}}E_{j}+B^{*}_{j}\hat{\bm{\mathcal{L}}}B_{j}\big{)}.bold_italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT + bold_italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = divide start_ARG italic_g end_ARG start_ARG 2 end_ARG roman_Re ( italic_E start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT over^ start_ARG bold_caligraphic_L end_ARG italic_E start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_B start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT over^ start_ARG bold_caligraphic_L end_ARG italic_B start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) . (8)

Here g=(8πω)1𝑔superscript8𝜋𝜔1g=(8\pi\omega)^{-1}italic_g = ( 8 italic_π italic_ω ) start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT in Gaussian units, the OAM operator 𝓛^=𝒓×1i^𝓛𝒓1𝑖bold-∇\hat{\bm{\mathcal{L}}}=\bm{r}\times\frac{1}{i}\bm{\nabla}over^ start_ARG bold_caligraphic_L end_ARG = bold_italic_r × divide start_ARG 1 end_ARG start_ARG italic_i end_ARG bold_∇ and Einstein summation convention is used for the three-component vector.

Substituting Eqs. (2)-(4) into Eqs. (6)-(8), we can calculate the time-averaged densities of energy, SAM and OAM. The transverse angular momentum is the focus of attention in this study. Consequently, only the y𝑦yitalic_y-components of the SAM and OAM densities should be taken into consideration for the light propagating in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane. In fact, the SAM density is totally transverse because it is non-zero solely along the y𝑦yitalic_y-direction. Omitting the common factor gπ2𝑔superscript𝜋2g\pi^{2}italic_g italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and writing in explicit contributions from electric and magnetic fields, we obtain finally

Wesuperscript𝑊𝑒\displaystyle W^{e}italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT =ω(Jl12+Jl+12),Wm=2ωJl2,formulae-sequenceabsent𝜔superscriptsubscript𝐽𝑙12superscriptsubscript𝐽𝑙12superscript𝑊𝑚2𝜔superscriptsubscript𝐽𝑙2\displaystyle=\omega(J_{l-1}^{2}+J_{l+1}^{2}),~{}\qquad\qquad~{}W^{m}=2\omega J% _{l}^{2},= italic_ω ( italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = 2 italic_ω italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (9)
Syesubscriptsuperscript𝑆𝑒𝑦\displaystyle S^{e}_{y}italic_S start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT =Jl12Jl+12,Sym=0,formulae-sequenceabsentsuperscriptsubscript𝐽𝑙12superscriptsubscript𝐽𝑙12subscriptsuperscript𝑆𝑚𝑦0\displaystyle=J_{l-1}^{2}-J_{l+1}^{2},~{}~{}\qquad~{}~{}~{}~{}\qquad~{}S^{m}_{% y}=0,= italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_S start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 0 , (10)
Lyesubscriptsuperscript𝐿𝑒𝑦\displaystyle L^{e}_{y}italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT =(l1)Jl12+(l+1)Jl+12,Lym=2lJl2.formulae-sequenceabsent𝑙1superscriptsubscript𝐽𝑙12𝑙1superscriptsubscript𝐽𝑙12subscriptsuperscript𝐿𝑚𝑦2𝑙superscriptsubscript𝐽𝑙2\displaystyle=(l-1)J_{l-1}^{2}+(l+1)J_{l+1}^{2},~{}~{}L^{m}_{y}=2lJ_{l}^{2}.= ( italic_l - 1 ) italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_l + 1 ) italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = 2 italic_l italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (11)

From Eq. (10), we see that the t-SAM density is completely provided by the roatation of the electric field. The normalized electric SAM density ωSy/We𝜔subscript𝑆𝑦superscript𝑊𝑒\omega S_{y}/W^{e}italic_ω italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT can change sign from -- to +++ and its minimum and maximum values are 11-1- 1 and 1111, respectively; see Fig. 3(a). The rotation direction of the electric field then converts from the negative to the positive φ𝜑\varphiitalic_φ-direction relative to the y𝑦yitalic_y-axis. Namely, the electric field can shift from the left-handed to the right-handed circular polarization, and vice versa. For the magnetic sector, we can identify the relation Lym=lWm/ωsuperscriptsubscript𝐿𝑦𝑚𝑙superscript𝑊𝑚𝜔L_{y}^{m}=lW^{m}/\omegaitalic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT = italic_l italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT / italic_ω, and so we can confirm that the t-OAM coincides well with l𝑙Planck-constant-over-2-pil\hbaritalic_l roman_ℏ per photon as the quantization of photon’s energy ωPlanck-constant-over-2-pi𝜔\hbar\omegaroman_ℏ italic_ω. However, for the t-OAM density of the electric contribution, there does not exist a similarly direct correlation between Lyesubscriptsuperscript𝐿𝑒𝑦L^{e}_{y}italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT and Wesuperscript𝑊𝑒W^{e}italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT. Interestingly, from Eqs. (9)-(11) we can find the total transverse angular momentum density of the electric sector satisfies

𝒥ye=Lye+Sye=lωWe.subscriptsuperscript𝒥𝑒𝑦subscriptsuperscript𝐿𝑒𝑦subscriptsuperscript𝑆𝑒𝑦𝑙𝜔superscript𝑊𝑒\mathcal{J}^{e}_{y}=L^{e}_{y}+S^{e}_{y}=\frac{l}{\omega}W^{e}.caligraphic_J start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + italic_S start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = divide start_ARG italic_l end_ARG start_ARG italic_ω end_ARG italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT . (12)

Again, for the electric sector the total transverse angular momentum still corresponds to l𝑙Planck-constant-over-2-pil\hbaritalic_l roman_ℏ per photon. By Eq. (12), the normalized electric t-OAM density is then given by ωLy/We=lωSy/We𝜔subscript𝐿𝑦superscript𝑊𝑒𝑙𝜔subscript𝑆𝑦superscript𝑊𝑒\omega L_{y}/W^{e}=l-\omega S_{y}/W^{e}italic_ω italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT = italic_l - italic_ω italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT and so it can be regarded as varying with ωSy/We𝜔subscript𝑆𝑦superscript𝑊𝑒\omega S_{y}/W^{e}italic_ω italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT if l𝑙litalic_l is given, and vice versa. It achieves its minimum and maximum values of l1𝑙1l-1italic_l - 1 and l+1𝑙1l+1italic_l + 1; see Fig. 3(c).

In addition, the densities of energy, t-SAM and t-OAM display strong electric-magnetic asymmetry [6], especially the t-SAM density. Not coincidentally, the electric-magnetic asymmetry can be found in the second mode (5) as well. However, in such case the electric field will make no contribution to the t-SAM density and so only the magnetic field will generate the t-SAM density. In particular, if l=0𝑙0l=0italic_l = 0, the optical structure is so trivial that the vortex does not occur and both the densities of SAM and OAM vanish.

Refer to caption

(a)                             (b) Refer to caption (c)                             (d)

Figure 3: Spatial distributions of the normalized t-SAM and t-OAM. (a) ωSye/We𝜔superscriptsubscript𝑆𝑦𝑒superscript𝑊𝑒\omega S_{y}^{e}/W^{e}italic_ω italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT, calculated via Eqs. (9) and (10). (b) ωSy/W𝜔subscript𝑆𝑦𝑊\omega S_{y}/Witalic_ω italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W, Eqs. (13) and (14). (c) ωLye/We𝜔superscriptsubscript𝐿𝑦𝑒superscript𝑊𝑒\omega L_{y}^{e}/W^{e}italic_ω italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT / italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT, Eqs. (9) and (11). (d) ωLy/W𝜔subscript𝐿𝑦𝑊\omega L_{y}/Witalic_ω italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT / italic_W, Eqs. (13) and (15). The parameter is l=3𝑙3l=3italic_l = 3.

Evidently, the densities of energy, t-SAM and t-OAM contributed by the electric and magnetic fields together can be written as

W𝑊\displaystyle Witalic_W =ω(Jl12+2Jl2+Jl+12),absent𝜔superscriptsubscript𝐽𝑙122superscriptsubscript𝐽𝑙2superscriptsubscript𝐽𝑙12\displaystyle=\omega(J_{l-1}^{2}+2J_{l}^{2}+J_{l+1}^{2}),= italic_ω ( italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (13)
Sysubscript𝑆𝑦\displaystyle S_{y}italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT =Jl12Jl+12,absentsuperscriptsubscript𝐽𝑙12superscriptsubscript𝐽𝑙12\displaystyle=J_{l-1}^{2}-J_{l+1}^{2},= italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (14)
Lysubscript𝐿𝑦\displaystyle L_{y}italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT =(l1)Jl12+2lJl2+(l+1)Jl+12.absent𝑙1superscriptsubscript𝐽𝑙122𝑙superscriptsubscript𝐽𝑙2𝑙1superscriptsubscript𝐽𝑙12\displaystyle=(l-1)J_{l-1}^{2}+2lJ_{l}^{2}+(l+1)J_{l+1}^{2}.= ( italic_l - 1 ) italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + 2 italic_l italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + ( italic_l + 1 ) italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (15)

The normalized t-SAM and t-OAM density are depicted in Fig. 3(b) and (d). From Eqs. (13)-(15), we obtain easily the total transverse angular momentum density

𝒥y=Ly+Sy=lωW.subscript𝒥𝑦subscript𝐿𝑦subscript𝑆𝑦𝑙𝜔𝑊\mathcal{J}_{y}=L_{y}+S_{y}=\frac{l}{\omega}W.caligraphic_J start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT + italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT = divide start_ARG italic_l end_ARG start_ARG italic_ω end_ARG italic_W . (16)

Once again, the total transverse angular momentum corresponds to l𝑙Planck-constant-over-2-pil\hbaritalic_l roman_ℏ per photon. It is natural to interpret Eqs. (12) and  (16) as representing conservation of the transverse angular momentum. Such conservation manifests a new type of mutual conversion of the t-SAM and t-OAM without light-matter interaction [18]. Integrating Eq. (16) over the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane, we obtain the normalized integral values of total transverse angular momentum

𝒥yW=lω.delimited-⟨⟩subscript𝒥𝑦delimited-⟨⟩𝑊𝑙𝜔\frac{\langle\mathcal{J}_{y}\rangle}{\langle W\rangle}=\frac{l}{\omega}.divide start_ARG ⟨ caligraphic_J start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ end_ARG start_ARG ⟨ italic_W ⟩ end_ARG = divide start_ARG italic_l end_ARG start_ARG italic_ω end_ARG . (17)

Here delimited-⟨⟩\langle\cdots\rangle⟨ ⋯ ⟩ means the integration over the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane, =dzdx=002πρdρdφdelimited-⟨⟩subscriptsuperscriptsubscriptsuperscriptdifferential-d𝑧differential-d𝑥subscriptsuperscript0subscriptsuperscript2𝜋0𝜌differential-d𝜌differential-d𝜑\langle\cdots\rangle=\int^{\infty}_{-\infty}\int^{\infty}_{-\infty}\cdots{\rm d% }z{\rm d}x=\int^{\infty}_{0}\int^{2\pi}_{0}\cdots\rho{\rm d}\rho{\rm d}\varphi⟨ ⋯ ⟩ = ∫ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT ∫ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - ∞ end_POSTSUBSCRIPT ⋯ roman_d italic_z roman_d italic_x = ∫ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ∫ start_POSTSUPERSCRIPT 2 italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ⋯ italic_ρ roman_d italic_ρ roman_d italic_φ. As seen from Eq. (17), the integral total transverse angular momentum can be defined well with l𝑙Planck-constant-over-2-pil\hbaritalic_l roman_ℏ per photon. More specifically, such integral transverse angular momentum merely stands for the intrinsic t-OAM, which will be elucidated in the remainder of this section.

Refer to caption

(a)                              (b) Refer to caption (c)

Figure 4: Spin-momentum locking for the rotating-wave structured light. (a)-(b) The spatial distributions of the Poynting momentum density with l=±3𝑙plus-or-minus3l=\pm 3italic_l = ± 3, respectively, as given by Eq (20). Arrows are used to indicate the direction of the Poynting momentum density. (c) The spatial distributions of the t-SAM Sysubscript𝑆𝑦S_{y}italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT along the z𝑧zitalic_z-axis with l=±3𝑙plus-or-minus3l=\pm 3italic_l = ± 3, respectively, Eq (14).

It is worth further looking into the t-SAM. Following Refs. [15] and [16], the connection between the t-SAM density and the Poynting momentum density for an arbitrary surface evanescent wave can be established as follows:

𝑺=12k2×𝑷,𝑺12superscript𝑘2bold-∇𝑷\bm{S}=\frac{1}{2k^{2}}\bm{\nabla}\times\bm{P},bold_italic_S = divide start_ARG 1 end_ARG start_ARG 2 italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG bold_∇ × bold_italic_P , (18)

where the Poynting momentum density takes the form

𝑷=gkRe(𝑬×𝑩).𝑷𝑔𝑘Resuperscript𝑬𝑩\bm{P}=gk{\rm Re\,}\big{(}{\bm{E}}^{*}\times{\bm{B}}\big{)}.bold_italic_P = italic_g italic_k roman_Re ( bold_italic_E start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT × bold_italic_B ) . (19)

Plugging Eqs. (2)-(4) into Eq. (19), the Poynting momentum density is given by

𝑷=4lρJl2cosφ𝒆x4lρJl2sinφ𝒆z=4lρJl2𝒆φ.𝑷4𝑙𝜌superscriptsubscript𝐽𝑙2𝜑subscript𝒆𝑥4𝑙𝜌superscriptsubscript𝐽𝑙2𝜑subscript𝒆𝑧4𝑙𝜌superscriptsubscript𝐽𝑙2subscript𝒆𝜑\bm{P}=\frac{4l}{\rho}J_{l}^{2}\cos\varphi\bm{e}_{x}-\frac{4l}{\rho}J_{l}^{2}% \sin\varphi\bm{e}_{z}=\frac{4l}{\rho}J_{l}^{2}\bm{e}_{\varphi}.bold_italic_P = divide start_ARG 4 italic_l end_ARG start_ARG italic_ρ end_ARG italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_cos italic_φ bold_italic_e start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT - divide start_ARG 4 italic_l end_ARG start_ARG italic_ρ end_ARG italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT roman_sin italic_φ bold_italic_e start_POSTSUBSCRIPT italic_z end_POSTSUBSCRIPT = divide start_ARG 4 italic_l end_ARG start_ARG italic_ρ end_ARG italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT bold_italic_e start_POSTSUBSCRIPT italic_φ end_POSTSUBSCRIPT . (20)

Here the common factor gπ2𝑔superscript𝜋2g\pi^{2}italic_g italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT has been omitted. From Eq. (20), we find that the Poynting momentum density is directed along purely the azimuthal direction in the (z,x)𝑧𝑥(z,x)( italic_z , italic_x ) plane, as illustrated in Fig. 4(a)-(b). Therefore, the y𝑦yitalic_y-components of SAM and OAM are indeed transverse because of their orientations perpendicular to the Poynting momentum density as well as the propagation direction.

With elementary algebra, we can verify that the connection (18) still holds for the t-SAM density (14) and the Poynting momentum (19) despite non-evanescent wave. Such connection offers a guideline for understanding the t-SAM from the perspective of the transverse gradient of Poynting momentum. The t-SAM here is the interference result of plane waves with different wave-vector directions in the same plane, even taking place in the very simple interference systems of two [10] or three [11] plane waves. In our instance here, the t-SAM density reversals the local orientations universally once the Poynting momentum changes from the positive to the negative φ𝜑\varphiitalic_φ-direction, as illustrated in Fig. 4(c). Thus, the direction of the t-SAM density is determined by the direction of Poynting momentum density, which can be derived clearly from the connection (18). From the connection between the t-SAM density and the Poynting momentum, we show that the rotating-wave structured light reveals an intrinsic spin-momentum-locking effect [17, 18] without light-matter interaction.

Let us coming back to the integral transverse angular momentum. In general, the t-SAM of helicity independence exists locally, and yet its integral value vanishes [6]. The integral t-SAM in Eq. (14) can be verified to go to exactly zero:

Sydelimited-⟨⟩subscript𝑆𝑦\displaystyle\langle S_{y}\rangle⟨ italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ =02πdφ0(Jl12Jl+12)ρdρabsentsubscriptsuperscript2𝜋0differential-d𝜑subscriptsuperscript0superscriptsubscript𝐽𝑙12superscriptsubscript𝐽𝑙12𝜌differential-d𝜌\displaystyle=\int^{2\pi}_{0}{\rm d}\varphi\int^{\infty}_{0}(J_{l-1}^{2}-J_{l+% 1}^{2})\rho{\rm d}\rho= ∫ start_POSTSUPERSCRIPT 2 italic_π end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT roman_d italic_φ ∫ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( italic_J start_POSTSUBSCRIPT italic_l - 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - italic_J start_POSTSUBSCRIPT italic_l + 1 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) italic_ρ roman_d italic_ρ
=8πlk20Jl(kρ)dJl(kρ)d(kρ)d(kρ)=0.absent8𝜋𝑙superscript𝑘2subscriptsuperscript0subscript𝐽𝑙𝑘𝜌dsubscript𝐽𝑙𝑘𝜌d𝑘𝜌d𝑘𝜌0\displaystyle=\frac{8\pi l}{k^{2}}\int^{\infty}_{0}J_{l}(k\rho)\frac{{\rm d}J_% {l}(k\rho)}{{\rm d}(k\rho)}{\rm d}(k\rho)=0.= divide start_ARG 8 italic_π italic_l end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ( italic_k italic_ρ ) divide start_ARG roman_d italic_J start_POSTSUBSCRIPT italic_l end_POSTSUBSCRIPT ( italic_k italic_ρ ) end_ARG start_ARG roman_d ( italic_k italic_ρ ) end_ARG roman_d ( italic_k italic_ρ ) = 0 . (21)

Then the integral total transverse angular momentum can be written as

𝒥y=Ly+Sy=Ly=lωW.delimited-⟨⟩subscript𝒥𝑦delimited-⟨⟩subscript𝐿𝑦delimited-⟨⟩subscript𝑆𝑦delimited-⟨⟩subscript𝐿𝑦𝑙𝜔delimited-⟨⟩𝑊\langle\mathcal{J}_{y}\rangle=\langle L_{y}\rangle+\langle S_{y}\rangle=% \langle L_{y}\rangle=\frac{l}{\omega}\langle W\rangle.⟨ caligraphic_J start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ = ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ + ⟨ italic_S start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ = ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ = divide start_ARG italic_l end_ARG start_ARG italic_ω end_ARG ⟨ italic_W ⟩ . (22)

Hence the integral transverse angular momentum belongs entirely to the t-OAM.

An alternate approach to derive the above result is giving by performing directly the integrations of energy Wedelimited-⟨⟩superscript𝑊𝑒\langle W^{e}\rangle⟨ italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT ⟩ and Wmdelimited-⟨⟩superscript𝑊𝑚\langle W^{m}\rangle⟨ italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⟩ and the t-OAM Lyedelimited-⟨⟩subscriptsuperscript𝐿𝑒𝑦\langle L^{e}_{y}\rangle⟨ italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ and Lymdelimited-⟨⟩subscriptsuperscript𝐿𝑚𝑦\langle L^{m}_{y}\rangle⟨ italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩. Using the closure equation of Bessel Functions of the first kind [47], we have

We=Wm,Lye=Lym.formulae-sequencedelimited-⟨⟩superscript𝑊𝑒delimited-⟨⟩superscript𝑊𝑚delimited-⟨⟩subscriptsuperscript𝐿𝑒𝑦delimited-⟨⟩subscriptsuperscript𝐿𝑚𝑦\langle W^{e}\rangle=\langle W^{m}\rangle,~{}~{}\langle L^{e}_{y}\rangle=% \langle L^{m}_{y}\rangle.⟨ italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT ⟩ = ⟨ italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⟩ , ⟨ italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ = ⟨ italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ . (23)

Immediately, we obtain [48]

LyW=LyeWe=LymWm=lω.delimited-⟨⟩subscript𝐿𝑦delimited-⟨⟩𝑊delimited-⟨⟩subscriptsuperscript𝐿𝑒𝑦delimited-⟨⟩superscript𝑊𝑒delimited-⟨⟩subscriptsuperscript𝐿𝑚𝑦delimited-⟨⟩superscript𝑊𝑚𝑙𝜔\frac{\langle L_{y}\rangle}{\langle W\rangle}=\frac{\langle L^{e}_{y}\rangle}{% \langle W^{e}\rangle}=\frac{\langle L^{m}_{y}\rangle}{\langle W^{m}\rangle}=% \frac{l}{\omega}.divide start_ARG ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ end_ARG start_ARG ⟨ italic_W ⟩ end_ARG = divide start_ARG ⟨ italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ end_ARG start_ARG ⟨ italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT ⟩ end_ARG = divide start_ARG ⟨ italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ end_ARG start_ARG ⟨ italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⟩ end_ARG = divide start_ARG italic_l end_ARG start_ARG italic_ω end_ARG . (24)

Thus, the relation (22) between the integral t-OAM and energy is recovered, and the integral t-SAM vanishes.

We are now in the position to interpret why the integral t-OAM is intrinsic. For a translation of the coordinate origin 𝒓𝒓+𝒓0𝒓𝒓subscript𝒓0\bm{r}\to\bm{r}+\bm{r}_{0}bold_italic_r → bold_italic_r + bold_italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, the integral t-OAM is changed by LyLy+Lyext=Ly+(𝒓0×𝑷)ydelimited-⟨⟩subscript𝐿𝑦delimited-⟨⟩subscript𝐿𝑦delimited-⟨⟩superscriptsubscript𝐿𝑦extdelimited-⟨⟩subscript𝐿𝑦subscriptsubscript𝒓0delimited-⟨⟩𝑷𝑦\langle L_{y}\rangle\to\langle L_{y}\rangle+\langle L_{y}^{\rm ext}\rangle=% \langle L_{y}\rangle+\big{(}\bm{r}_{0}\times\langle\bm{P}\rangle\big{)}_{y}⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ → ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ + ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ext end_POSTSUPERSCRIPT ⟩ = ⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ + ( bold_italic_r start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT × ⟨ bold_italic_P ⟩ ) start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT. From Eq. (20) we find 𝑷=0delimited-⟨⟩𝑷0\langle\bm{P}\rangle=0⟨ bold_italic_P ⟩ = 0 and so the extrinsic t-OAM Lyextdelimited-⟨⟩superscriptsubscript𝐿𝑦ext\langle L_{y}^{\rm ext}\rangle⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ext end_POSTSUPERSCRIPT ⟩ is equal to zero. The integral t-OAM is unchanged under the translation of the coordinate origin and therefore it is intrinsic.

V Discussion and Summary

In summary, we have shown that the rotating-wave structured light carries purely integral intrinsic t-OAM within monochromatic optical fields. Furthermore, the rotating-wave structured light with wave vortex, consisting of linearly-polarized plane waves, exhibits locally unique t-SAM and t-OAM features. The sum of the locally-normalized t-SAM and t-OAM is conserved universally, which manifests a strong coupling of t-SAM and t-OAM. The t-SAM and t-OAM variations are heavily dependent on each other. Such strong coupling uncovers a new type of mutual conversion of the intrinsic t-SAM and t-OAM without light-matter interaction. In spite of such locally mutual conversion, the integral t-SAM and t-OAM are rather robust, and the integral intrinsic transverse angular momentum belong purely to the intrinsic t-OAM.

Impressively, the t-SAM density of the rotating-wave structured light also can be directly linked to the curl of the Poynting momentum density even though it is not associated with evanescent waves. As a consequence, the orientation of the local t-SAM is made directly subject to the local direction of the Poynting momentum density, which can be recognized as an intrinsic spin-momentum-locking effect of the rotating-wave structured light in free space.

We look forward to the future experimental verification of such nontrivial properties of the spin and orbital angular momentum in the rotating-wave structured light. The rotating-wave structured light might offer high potential for applications in the optical manipulation of matter.

Acknowledgments

The author is grateful to Z. Xiong and Y. Liu for fruitful discussions. This work was funded by the Scientific Research Project of Hubei Polytechnic University (No. 20xjz02R).

References

  • [1] S. Franke-Arnold, Nat. Rev. Phys. 4, 361 (2022).
  • [2] J.H. Poynting, Proc. R. Soc. Lond. A 82, 560 (1909).
  • [3] L. Allen, M.W. Beijersbergen, R.J.C. Spreeuw, and J.P. Woerdman, Phys. Rev. A 45, 8185 (1992).
  • [4] A.T. O’Neil, I. MacVicar, L. Allen, and M.J. Padgett, Phys. Rev. Lett. 88, 053601 (2002).
  • [5] V. Garcés-Chávez, D. McGloin, M.J. Padgett, W. Dultz, H. Schmitzer, and K. Dholakia, Phys. Rev. Lett. 91, 093602 (2003).
  • [6] K.Y. Bliokh and F. Nori, Phys. Rep. 592, 1 (2015).
  • [7] K.Y. Bliokh and F. Nori, Phys. Rev. A 85, 061801 (2012).
  • [8] K.-Y. Kim, I.-M. Lee, J. Kim, J. Jung, and B. Lee, Phys. Rev. A 86, 063805 (2012).
  • [9] K.Y. Bliokh, A.Y. Bekshaev, and F. Nori, Nat. Commun. 5, 3300 (2014).
  • [10] A.Y. Bekshaev, K.Y. Bliokh, and F. Nori, Phys. Rev. X 5, 011039 (2015).
  • [11] Z.-L. Wang, D.-D. Lian, and X.-S. Chen, Opt. Express 26, 33712 (2018).
  • [12] N. Yang and A.E. Cohen, J. Phys. Chem. B 115, 5304 (2011).
  • [13] R. Mathevet and G.L.J.A. Rikken, Opt. Mater. Express 4, 2574 (2014).
  • [14] M. Neugebauer, T. Bauer, A. Aiello, and P. Banzer, Phys. Rev. Lett. 114, 063901 (2015).
  • [15] P. Shi, L. Du, C. Li, A.V. Zayats, and X. Yuan, Proc. Natl. Acad. Sci. U.S.A. 118, e2018816118 (2021).
  • [16] P. Shi, L. Du, A. Yang, X. Yin, X. Lei, and X. Yuan, Commun. Phys. 6, 283 (2023).
  • [17] A. Aiello, P. Banzer, M. Neugebauer, and G. Leuchs, Nat. Photonics 9, 789 (2015).
  • [18] P. Shi, L. Du, and X. Yuan, Nanophotonics 10, 3927 (2021).
  • [19] K.Y. Bliokh, D. Smirnova, and F. Nori, Science 348, 1448 (2015).
  • [20] M. Onoda, S. Murakami, and N. Nagaosa, Phys. Rev. Lett. 93, 083901 (2004).
  • [21] K.Y. Bliokh and Y.P. Bliokh, Phys. Rev. Lett. 96, 073903 (2006).
  • [22] K.Y. Bliokh, Phys. Rev. Lett. 97, 043901 (2006).
  • [23] O. Hosten and P. Kwiat, Science 319, 787 (2008).
  • [24] K.Y. Bliokh and A. Aiello, J. Opt. 15, 014001 (2013)
  • [25] A. Aiello, N. Lindlein, C. Marquardt, and G. Leuchs, Phys. Rev. Lett. 103, 100401 (2009).
  • [26] L.-J. Kong, S.-X. Qian, Z.-C. Ren, X.-L. Wang, and H.-T. Wang, Phys. Rev. A 85, 035804 (2012).
  • [27] J. Korger, A. Aiello, V. Chille, P. Banzer, C. Wittmann, N. Lindlein, C. Marquardt, and G. Leuchs, Phys. Rev. Lett. 112, 113902 (2014).
  • [28] X. Ling, X. Zhou, K. Huang, Y. Liu, C.-W. Qiu, H. Luo, and S. Wen, Rep. Prog. Phys. 80, 066401 (2017).
  • [29] Z.-L. Wang and X.-S. Chen, Phys. Rev. A 99, 063832 (2019).
  • [30] K.Y. Bliokh and F. Nori, Phys. Rev. A 86, 033824 (2012).
  • [31] K.Y. Bliokh, Phys. Rev. Lett. 126, 243601 (2021).
  • [32] S.W. Hancock, S. Zahedpour, and H.M. Milchberg, Phys. Rev. Lett. 127, 193901 (2021).
  • [33] K.Y. Bliokh, Phys. Rev. A 107, L031501 (2023).
  • [34] S.W. Hancock, S. Zahedpour, A. Goffin, and H.M. Milchberg, Optica 6, 1547 (2019).
  • [35] A. Chong, C. Wan, J. Chen, and Q. Zhan, Nat. Photonics 14, 350 (2020).
  • [36] G. Gui, N.J. Brooks, H.C. Kapteyn, M.M. Murnane, and C.T. Liao, Nat. Photonics 15, 608 (2021).
  • [37] Y. Fang, S. Lu, and Y. Liu, Phys. Rev. Lett. 127, 273901 (2021).
  • [38] Y. Hu, J.J. Kingsley-Smith, M. Nikkhou, J.A. Sabin, F.J. Rodríguez-Fortuño, X. Xu, and J. Millen, Nat. Commun. 14, 2638 (2023).
  • [39] P.H. Ceperley, Rotating Waves, Am. J. Phys. 60, 938 (1992).
  • [40] P.H. Ceperley and A. Koren, J. Acoust. Soc. Am. 91, 2331 (1992).
  • [41] A.O. Santillán and K. Volke-Sepúlveda, Am. J. Phys. 77, 209 (2009).
  • [42] Z. Hong, J. Zhang, and B.W. Drinkwater, Phys. Rev. Lett. 114, 214301 (2015).
  • [43] J.E. Velazco and P.H. Ceperley, IEEE Trans. Microwave Theory Techn. 41, 330 (1993).
  • [44] S. Chávez-Cerda, G.S. McDonald, and G.H.C. New, Opt. Commun. 123, 225 (1996).
  • [45] M.G. Calkin, Am. J. Phys. 33, 958 (1965).
  • [46] J.D. Jackson, Classical electrodynamics, 3rd ed. (Wiley, Hoboken, NJ, 1999).
  • [47] G.B. Arfken and H.J. Weber, Mathematical methods for physicists, 6th ed. (Elsevier Academic Press, New York, 2005).
  • [48] Despite We=Wmdelimited-⟨⟩superscript𝑊𝑒delimited-⟨⟩superscript𝑊𝑚\langle W^{e}\rangle=\langle W^{m}\rangle⟨ italic_W start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT ⟩ = ⟨ italic_W start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ⟩ and Lye=Lymdelimited-⟨⟩subscriptsuperscript𝐿𝑒𝑦delimited-⟨⟩subscriptsuperscript𝐿𝑚𝑦\langle L^{e}_{y}\rangle=\langle L^{m}_{y}\rangle⟨ italic_L start_POSTSUPERSCRIPT italic_e end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ = ⟨ italic_L start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ are infinite, they indeed can be rigorously defined by the Dirac δ𝛿\deltaitalic_δ-function, and the normalized value of integral t-OAM Ly/Wdelimited-⟨⟩subscript𝐿𝑦delimited-⟨⟩𝑊\langle L_{y}\rangle/\langle W\rangle⟨ italic_L start_POSTSUBSCRIPT italic_y end_POSTSUBSCRIPT ⟩ / ⟨ italic_W ⟩ is exactly l/ω𝑙𝜔l/\omegaitalic_l / italic_ω.