Introduction

Optical angular momentum is contributed to by physical orbital and spin parts1,2,3,4. The orbital component can be either extrinsic or intrinsic and is associated with the spatial structure of the light field. Orbital angular momentum (OAM) is most famously stirred into light by optical vortices, which have helical phase fronts twisted around a one-dimensional line singularity5,6,7,8. Spin angular momentum (SAM), meanwhile, develops as the electric and magnetic field vectors rotate during oscillation9, and is an intrinsic quantity, independent of co-ordinate origin. Both OAM and SAM arise from chiral structures in light, and as such, can access the chirality of matter, enhancing, pushing, twisting, and torquing with or without a parity bias10,11. Many commonly observed chiroptical interactions, such as circular dichroism, arise from the coupling of chiral matter to a photon’s spin state. Recent attention has also turned to OAM-dependent photon absorption and scattering, with a range of effects theorised12,13,14,15. Chiroptical effects involving OAM often hinge on the strong longitudinal field component introduced by beam focussing16,17. Unlike for photons with opposite spins, however, preferential absorption of photons with different OAM handedness does not occur in the dipole approximation for paraxial light18,19.

A well-known orbital and spin separation can be performed locally in the Poynting vector20 (kinetic momentum density, when divided by ({c}^{2})). The Poynting vector may be expressed as a sum of the orbital current, pointing in the direction of canonical momentum, and spin current, proportional to the virtual spin momentum. This is a physically meaningful decomposition, tying in to the famous Abraham-Minkowski dilemma21,22.

Surprisingly to most readers, a similar vector decomposition can be performed on the electromagnetic spin density ({bf{S}}=left(1/4omega right){text{Im}}{{epsilon }_{0}{{bf{E}}}^{{boldsymbol{* }}}times{bf{E}}+{mu }_{0}{{bf{H}}}^{{boldsymbol{* }}}times{bf{H}}}), which is split into, perhaps confusingly, orbital-like and spin-like contributions to the total electromagnetic spin. We name these two terms canonical ({{bf{s}}}_{{rm{c}}}) and Poynting spin ({{bf{s}}}_{{rm{p}}}), such that ({bf{S}}={{bf{s}}}_{{rm{c}}}+{{bf{s}}}_{{rm{p}}}). This little-known decomposition is the focus of our work—amidst discussion by Shi et al. who, without giving the complete decomposition explicitly, linked the two terms to longitudinal and transverse spin in some specific cases23,24,25, and Bekshaev, who discussed the validity of these links for a range of scenarios finding that they prove less tenuous with modifications26, no clear picture exists in the literature on the physical significance of the two terms in general, in the same way as for the decomposed Poynting vector. The decomposition as we present has been expressed previously (e.g., in ref.27, where the authors provided the decomposition in units of momentum in part to discuss chiral optical forces, but while neglecting recoil terms), but with almost no discussion on the function of terms as component parts of ({bf{S}}).

We therefore have two objectives: to clarify and understand the meaning of the orbit-like and spin-like components of ({bf{S}}), and to emphasise their value for the research community in interpreting electromagnetic SAM density. For example, by use of the Maxwell stress tensor, one understands that the canonical spin ({{bf{s}}}_{{rm{c}}}) is proportional to chiral pressure while Poynting spin manifests only in chiral recoil forces. Chiral pressure and therefore canonical spin is directly measurable, in a chiral analogy to the achiral and directly measurable radiation pressure developed by canonical momentum. While the two decomposed spin vectors are organised independently in general monochromatic light, a remarkable implication lies in optical vortices. We show that a linearly polarised vortex imparts a longitudinal component to canonical spin, resulting in longitudinal chiral pressure even in the absence of longitudinal total spin. This force can act even under the dipole approximation and should be present without strong focussing.

We generalise the spin decomposition to time-dependent fields and four-vector representation, unlocking a deeper understanding of the behaviour of the two terms. We also show that the ability to split spin into two terms appears to be a general feature of wave fields, similar to the decomposition of the Poynting vector, including those of the linearised theories of gravity28 and acoustics, which have distinctly different vector structures. Throughout, we use bold Latin letters ({bf{E}}) and ({bf{H}}) to represent time-harmonic complex electromagnetic field phasors, which are a function of position only [e.g., ({bf{E}}left({bf{r}}right)={{bf{E}}}_{0}exp left(i{bf{k}},{{cdot }},{bf{r}}right))], and scripted characters (boldsymbol{mathcal{E}}) and (boldsymbol{mathcal{H}}) denoting real, time-varying field vectors (e.g., ({boldsymbol{mathcal{E}}}left({bf{r}},tright)={rm{Re}}{{bf{E}}left({bf{r}}right)exp left(-iomega tright)})).

Result

Electromagnetic spin decomposition

We will first contextualise the decomposition of electromagnetic SAM density by recalling a well-known orbital and spin decomposition which exists in free-space and time-harmonic fields for the time-averaged Poynting vector,

$${bf{P}}=frac{1}{2}{text{Re}}left{{{bf{E}}}^{{boldsymbol{* }}}times{bf{H}}right}$$
(1)

which represents the flux of active power in the light field. The Poynting vector is separated into the orbital and spin currents20,22,

$$begin{array}{l}{bf{P}}={underbrace{displaystylefrac{c^2}{4omega}{rm{Im}}{{epsilon}_{0}{{bf{E}}}^{ast }cdot(nabla){bf{E}}+{mu}_{0}{{bf{H}}}^{ast }cdot(nabla){bf{H}}}}_{{mathbf{p}}_{mathrm{o}}}}+,{underbrace{displaystylefrac{c^2}{2}{nabla}timesoverbrace{frac{1}{4omega}{rm{Im}}{{epsilon}_{0}{{bf{E}}}^{ast }times{bf{E}}+{mu}_{0}{{bf{H}}}^{ast }times{bf{H}}}}^{bf{S}}}_{{mathbf{p}}_{mathrm{s}}}}end{array}$$
(2)

defined using the complex phasors ({bf{E}}) and ({bf{H}}). The inner product notation is ({bf{a}}cdot left(nabla right){bf{b}}={a}_{x}nabla {b}_{x}+{a}_{y}nabla {b}_{y}+{a}_{z}nabla {b}_{z}). When divided by ({c}^{2}), ({bf{P}}) has units of a momentum density and is termed the kinetic momentum density, while the orbital current ({{bf{p}}}_{{rm{o}}}) becomes the canonical momentum density, corresponding to the expectation value of linear momentum carried by photons at each point in space (unlike kinetic momentum, the canonical momentum is a directly measurable quantity). The net flow of the instantaneous Poynting vector into or out of a volume relates to the change in electromagnetic energy density over time, according to a continuity equation. An adjacent conserved quantity of light is its chirality, which has its own continuity equation and an associated chirality flux29,30,31. Chirality (chi) is one of an infinite hierarchy of conserved quantities in linear, non-dispersive media, to which helicity (h) also belongs9,31—helicity and chirality are often conflated because they (and their fluxes) are proportional by a factor of ({omega }^{2}) in monochromatic waves, though in general they are distinct. In monochromatic light, the flux of helicity corresponds to the time-averaged SAM density (often referred to simply as ‘spin’), given by,

$${bf{S}}=frac{1}{4omega }{text{Im}}left{{epsilon }_{0}{{bf{E}}}^{{boldsymbol{* }}}times{bf{E}}+{mu }_{0}{{bf{H}}}^{{boldsymbol{* }}}times{bf{H}}right}equiv {{bf{S}}}_{{rm{e}}}+{{bf{S}}}_{{rm{m}}}$$
(3)

Spin’s individual electric and magnetic contributions point in the normal direction to the electric and magnetic polarisation ellipses, drawn over time by the instantaneous vectors ({boldsymbol{mathcal{E}}}left({bf{r}},tright)) and ({boldsymbol{mathcal{H}}}left({bf{r}},tright)). Like the Poynting vector, and using the same procedure, spin ({bf{S}}) may be split into a sum of two vectors: the canonical spin ({{bf{s}}}_{{rm{c}}}) and the Poynting spin ({{bf{s}}}_{{rm{p}}}). This decomposition is our main focus and is, explicitly,

$$begin{array}{ll}{mathbf{S}}&=underbrace{displaystylefrac{1}{{4omega}^2}{rm{Re}left{{bf{E}}^{ast},{cdot},(nabla){bf{H}}-{bf{H}}^{ast}cdot(nabla){bf{E}}right}}}_{mathbf{s}_{rm{c}}}&+underbrace{displaystylefrac{1}{{2omega}^2}nablatimesoverbrace{frac{1}{2}{rm{Re}left{{bf{E}}^{ast}times{bf{H}}right}}}^{bf{P}}}_{mathbf{s}_{rm{p}}}end{array}$$
(4)

Any subsequent mention of ‘canonical spin’ in this work always refers to the first term, ({{bf{s}}}_{{rm{c}}}=(1/4{omega }^{2})text{Re}{{{bf{E}}}^{{boldsymbol{* }}}cdot left(nabla right){bf{H}}-{{bf{H}}}^{{boldsymbol{* }}}cdot left(nabla right){bf{E}}}), and ‘Poynting spin’ to the second, ({{bf{s}}}_{{rm{p}}}=(1/2{omega }^{2})nabla times {bf{P}}). Their sum ({bf{S}}) [Eq. (3)] we refer to as total spin. This naming scheme reflects that of the decomposed kinetic momentum, split into canonical momentum and spin momentum. The authors of ref. 23 proffer different expressions for the two decomposed terms which are not quite complete (we dissect their equations in the supplementary material).

Physical interpretation

Previous works23,24 have linked the canonical ({{bf{s}}}_{{rm{c}}}) and Poynting spin ({{bf{s}}}_{{rm{p}}}) vectors to the longitudinal and transverse spin of light, respectively. However, this is only true in a very limited number of cases such as linearly polarised evanescent waves. It cannot be a general result because well-defined longitudinal and transverse directions do not exist for many-plane wave or multiple-beam interference. In fact, it does not hold even when there is a well-defined longitudinal direction: a circularly polarised evanescent wave has both longitudinal and transverse spin components contained in the Poynting spin vector.

One way to gain an understanding of the physical significance of the canonical and Poynting spins is to study the interactions between light and chiral matter. Chiral light can exert preferential forces on enantiomers with opposite handedness, with direct proportionality to the light’s SAM. When light shines on a particle much smaller than the wavelength (Rayleigh regime) it gets polarised and acquires an electric ({{bf{p}}}_{{rm{d}}}) and a magnetic ({{bf{m}}}_{{rm{d}}}) dipole. In the linear regime, the dipole moments are proportional to the incident fields, ({{bf{p}}}_{{rm{d}}}={alpha }_{text{e}}varepsilon {bf{E}}+i{alpha }_{text{c}}{bf{H}}/c) and ({{bf{m}}}_{{rm{d}}}={alpha }_{text{m}}{bf{H}}-i{alpha }_{text{c}}{bf{E}}/eta), where ({alpha }_{text{e}}), ({alpha }_{text{m}}) and ({alpha }_{text{c}}) are the electric, magnetic, and chiral polarisabilities of the particle. The chiral polarisability ({alpha }_{text{c}}) is a pseudoscalar so it changes sign between the two enantiomers (mirror-reflected versions), and vanishes unless the matter is chiral. Under these assumptions, the illuminating light exerts a chiral optical force on the particle (which changes sign for the different enantiomers) given by32,33,34,35,36,37:

$$begin{array}{ll}{mathbf{F}}_{{rm{chiral}}}&=underbrace{omega nabla [{mathrm{Re}}({alpha }_{{rm{c}}})h]}_{rm{helicity}; rm{gradient}}+underbrace{2omega k[text{Im}({alpha }_{{rm{c}}}){{bf{s}}}_{{rm{c}}}]}_{rm{chiral}; rm{pressure}}&-underbrace{omega tfrac{{k}^{4}}{3pi }[{mathrm{Re}}({alpha }_{{rm{e}}}^{ast }{alpha }_{{rm{c}}}){{bf{s}}}_{{rm{e}}}+{mathrm{Re}}({alpha }_{{rm{m}}}^{ast }{alpha }_{{rm{c}}}){{bf{s}}}_{{rm{m}}}]}_{rm{spin}; rm{recoil}}end{array}$$
(5)

where (k=omega /c) is the wavenumber, (h=-{rm{Im}}left({{bf{E}}}^{* }cdot {bf{H}}right)/left(2omega cright)) is the cycle-averaged optical helicity density17,38, and ({{bf{S}}}_{{rm{e}}}) and ({{bf{S}}}_{{rm{m}}}) are the electric and magnetic parts of the total spin ({bf{S}}) from Eq. (3). One immediately sees that the canonical spin ({{bf{s}}}_{text{c}}), one of the two terms in the spin decomposition [Eq. (4)], appears in this force equation and is directly responsible for the chiral pressure (in this context, the product (k{{bf{s}}}_{text{c}}) is sometimes called the chiral momentum27), while the total spin appears in the spin recoil term. This chiral pressure is a chiral analogy to radiation pressure caused by canonical momentum and is transferred to matter through extinction. Poynting spin contributes (together with canonical spin) only to the higher-order recoil force caused by the unbalanced radiation pattern of an electric-magnetic dipole. This also implies that in an electromagnetic field whose SAM is pure Poynting spin (i.e., zero canonical spin), the particle will experience no chiral pressure, only helicity gradients and relatively weaker spin recoil terms of the force.

An alternative route to the physical meaning of the canonical and Poynting spins is found purely in the two terms’ mathematical expressions as shown in Eq. (4). The Poynting spin, being proportional to the curl of the Poynting vector, is easily interpreted as the vorticity in energy flow and in some sense couples the OAM of light to the total SAM density vector. The expression for the canonical spin is much harder to interpret initially but begins to unravel if we decompose a general electromagnetic field into circularly polarised plane waves. Any arbitrary electromagnetic field, no matter how complicated, can be expressed in momentum space by an angular spectrum (an infinite sum of plane waves of different wavevectors, weighted by an amplitude function). Using this property, we can further probe the canonical spin term in monochromatic light by separating electric and magnetic fields into two component fields of opposite helicities, indicated by the (+) and (-) subscripts,

$$begin{array}{ll}{mathbf{E}}left({mathbf{r}}right)&=iiint widetilde{{bf{E}}}left({bf{k}}right){{rm{e}}}^{i{bf{k}}cdot {bf{r}}}{{rm{d}}}^{3}k&=iiint left[{widetilde{E}}_{+}left({bf{k}}right){hat{{bf{e}}}}_{+}left({bf{k}}right)+{widetilde{E}}_{-}left({bf{k}}right){hat{{bf{e}}}}_{-}left({bf{k}}right)right]{{rm{e}}}^{i{bf{k}}cdot {bf{r}}}{{rm{d}}}^{3}k&={{bf{E}}}_{+}left({bf{r}}right)+{{bf{E}}}_{-}left({bf{r}}right)end{array}$$
(6)

where ({hat{{bf{e}}}}_{pm }left({bf{k}}right)) represent the circularly polarised unit vectors for each plane wave with wavevector ({bf{k}}) (see, e.g., ref. 39). The magnetic field’s associated helicity components are obtained from Faraday’s law (nabla times {bf{E}}=iomega mu {bf{H}}) as,

$$begin{array}{ll}{mathbf{H}}left({bf{r}}right)&=frac{1}{eta }iiint frac{{bf{k}}}{k}times widetilde{{bf{E}}}left({bf{k}}right){{rm{e}}}^{i{bf{k}}cdot {bf{r}}}{{rm{d}}}^{3}k&=frac{1}{eta }iiint left[-i{widetilde{E}}_{+}{hat{{bf{e}}}}_{+}left({bf{k}}right)+i{widetilde{E}}_{-}{hat{{bf{e}}}}_{-}left({bf{k}}right)right]{{rm{e}}}^{i{bf{k}}cdot {bf{r}}}{{rm{d}}}^{3}k&={{bf{H}}}_{+}left({bf{r}}right)+{{bf{H}}}_{-}left({bf{r}}right)end{array}$$
(7)

where we used the property (left({bf{k}}/kright)times {hat{{bf{e}}}}_{pm }=mp i{hat{{bf{e}}}}_{pm }). Equations (6) and (7) show that the helicity-separated electric and magnetic fields are related by ({{bf{H}}}_{pm }left({bf{r}}right)=mp i{{bf{E}}}_{pm }left({bf{r}}right)/eta), characteristic of a circularly polarised plane wave—this is true not only in the spectral representation, but also in the spatial representation, for any arbitrary field (a property explored in depth in40). This allows us to substitute helicity-separated fields into many dual quantities, including canonical momentum as in41 and canonical spin, gaining further insight. Simply substituting ({bf{E}}={{bf{E}}}_{+}+{{bf{E}}}_{-}) and ({bf{H}}=-ileft({{bf{E}}}_{+}-{{bf{E}}}_{-}right)/eta) into the expression for orbital current [({{bf{p}}}_{{rm{o}}}) of Eq. (2)], we find, after some algebra,

$$begin{array}{ll}{{bf{p}}}_{{rm{o}}}&=frac{{c}^{2}}{2omega }{epsilon }_{0}text{Im}left{{{bf{E}}}_{+}^{* }cdot left(nabla right){{bf{E}}}_{+}+{{bf{E}}}_{-}^{* }cdot left(nabla right){{bf{E}}}_{-}right}&={c}^{2}left({{bf{p}}}_{+}+{{bf{p}}}_{-}right)end{array}$$
(8)

showing that the helicity segregation of the ({bf{E}}) and ({bf{H}}) fields translates to a separation of contributions to the field’s momentum ({bf{p}}={{bf{p}}}_{+}+{{bf{p}}}_{-}) by photons of positive and negative helicity. Note that orbital current is proportional to canonical momentum by ({{bf{p}}}_{{rm{o}}}={c}^{2}{bf{p}}); here, ({{bf{p}}}_{+}) and ({{bf{p}}}_{-}) are helicity-separated momentum densities. Making the same substitution in the expression for canonical spin ({{bf{s}}}_{{rm{c}}}) [first term of Eq. (4)] brings a clear physical interpretation,

$$begin{array}{ll}{{bf{s}}}_{{rm{c}}}&={displaystylefrac{1}{4{omega }^{2}}}{text{Re}}left{-displaystylefrac{i}{eta }left[{{bf{E}}}_{{boldsymbol{+}}}^{{boldsymbol{* }}}+{{bf{E}}}_{{boldsymbol{-}}}^{{boldsymbol{* }}}right]cdot left(nabla right)left[{{bf{E}}}_{{boldsymbol{+}}}-{{bf{E}}}_{{boldsymbol{-}}}right]+frac{i}{eta }left[{{bf{E}}}_{{boldsymbol{+}}}^{{boldsymbol{* }}}-{{bf{E}}}_{{boldsymbol{-}}}^{{boldsymbol{* }}}right]cdot left(nabla right)left[{{bf{E}}}_{{boldsymbol{+}}}+{{bf{E}}}_{{boldsymbol{-}}}right]right}&={displaystylefrac{1}{2k{omega }}}{{epsilon }}_{0}{rm{Im}}left{{{bf{E}}}_{+}^{* }cdot left(nabla right){{bf{E}}}_{+}-{{bf{E}}}_{-}^{* }cdot left(nabla right){{bf{E}}}_{-}right}&=displaystylefrac{1}{k}left({{bf{p}}}_{+}-{{bf{p}}}_{-}right)end{array}$$
(9)

which is that canonical spin is proportional to the difference in linear momentum densities carried by photons of oppositely signed helicity. Summarising, only with the presence of canonical spin can a direct chiral interaction occur between a field and a chiral particle whereas both canonical and Poynting spins cause chiral matter to experience recoil forces, in close analogy with the way that canonical momentum and spin momentum contribute to achiral radiation pressure and recoil force. This is further supported by Eq. (9) because photons of opposite helicity are absorbed or scattered in different amounts by chiral particles. While electric and magnetic total spins (({{bf{S}}}_{{rm{e}}}) and ({{bf{S}}}_{{rm{m}}})) are oriented normal to the ({bf{E}}) and ({bf{H}}) polarisation ellipse, there is no concrete geometric requirement for canonical and Poynting spins to point in any direction relative to polarisation or momentum. In a general 3D field, such as in Fig. 1, ({{bf{p}}}_{+}-{{bf{p}}}_{-},ne, {{bf{p}}}_{+}+{{bf{p}}}_{-}) and hence canonical spin is decoupled from any local longitudinal direction defined by canonical momentum.

Fig. 1
figure 1

Component parts of light’s total spin angular momentum density visualised at a single point (red circle) in a 3D electromagnetic field. The field is generated by interfering randomly polarised plane waves with random (mathbf{k}) vectors, bearing minimal spatial symmetry and allowing the spin decomposition to be depicted in its most generic form. An ({xtext{‘}ytext{‘}}) cut plane of the field is provided, on which two elements are plotted: the field’s non-uniform energy density, indicated by the colour, and streamlines of the Poynting vector’s ({xtext{‘}ytext{‘}}) projection. The co-ordinate axes ({xtext{‘}}), ({ytext{‘}}) and ({ztext{‘}}) are specially chosen so that ({ztext{‘}}) aligns with the Poynting spin vector, purely to allow the curling of the Poynting vector streamlines to be visualised on the ({xtext{‘}ytext{‘}}) plane. The blue and green ellipses represent the electric and magnetic polarisation ellipses respectively, and are oriented freely in 3D (partially submerged by the cut plane as a 3D visual aid). Total spin (red arrow) is equal to the sum of the electric ({mathbf{S}}_{text{e}}) and magnetic ({mathbf{S}}_{text{m}}) spins and the sum of canonical spin (yellow arrow) and Poynting spin (black arrow) which have no orientation restrictions and all point in different 3D directions. Only ({mathbf{S}}_{text{e}}) and ({mathbf{S}}_{text{m}}) hold a geometric significance in relation to polarisation by being orthogonal to the plane of the (mathbf{E}) and (mathbf{H}) ellipses (according to the right-hand rule) as illustrated separately below the main combined diagram. Dotted lines within the ellipses represent the ellipse major and minor axes (note that the ellipses are distorted as the 3D field is projected into a 2D image)

Full size image

Mathematical derivation

In this section, we briefly lay out the mathematical steps taken to arrive at Eq. (4), before developing the expression into a more fundamental 4-vector description. For each of the canonical and Poynting spin terms, the 4-vector decomposition incorporates a time component that characterises how the two terms transform differently between reference frames.

3-vector decomposition

Maxwell’s equations enable both decompositions of the Poynting vector Eq. (2) and spin Eq. (4) by bringing the relevant vector from its usual representation into a form that can be separated into two terms using a vector identity,

$${bf{a}}times left(nabla times {bf{b}}right)={bf{a}}cdot left(nabla right){bf{b}}-left({bf{a}}cdot nabla right){bf{b}}$$
(10)

where ({bf{a}}) and ({bf{b}}) are arbitrary vectors, and the previously unseen notation on the right-hand side means (left({bf{a}}cdot nabla right){bf{b}}={a}_{x}{partial }_{x}{bf{b}}+{a}_{y}{partial }_{y}{bf{b}}+{a}_{z}{partial }_{z}{bf{b}}). In time-harmonic fields, the phasors ({bf{E}}) and ({bf{H}}) can be substituted for curls of their counterpart via Faraday’s and Ampere’s laws. Substituting for the un-conjugated phasors in the spin vector definition Eq. (3) gives an expression in the form of Eq. (10),

$${bf{S}}=frac{1}{4omega }{text{Im}}left{{epsilon }_{0}{{bf{E}}}^{{boldsymbol{* }}}times left(frac{i}{omega {epsilon }_{0}}nabla times {bf{H}}right)-{mu }_{0}{{bf{H}}}^{{boldsymbol{* }}}times left(frac{i}{omega {mu }_{0}}nabla times {bf{E}}right)right}$$
(11)

Applying the vector identity Eq. (10), we have,

$$begin{array}{ll}{mathbf{S}}&=frac{1}{4{omega }^{2}}text{Re}left{{{bf{E}}}^{{boldsymbol{* }}}cdot left(nabla right){bf{H}}-{{bf{H}}}^{{boldsymbol{* }}}cdot left(nabla right){bf{E}}right}&+frac{1}{4{omega }^{2}}text{Re}left{left({{bf{H}}}^{{boldsymbol{* }}}cdot nabla right){bf{E}}-left({{bf{E}}}^{{boldsymbol{* }}}cdot nabla right){bf{H}}right}end{array}$$
(12)

Gauss’ law in free-space ((nabla cdot {bf{E}}=0) and (nabla cdot {bf{H}}=0)), combined with a second vector identity (nabla times left({bf{a}}times {bf{b}}right)={bf{a}}left(nabla cdot {bf{b}}right)-{bf{b}}left(nabla cdot {bf{a}}right)+left({bf{b}}cdot nabla right){bf{a}}-left({bf{a}}cdot nabla right){bf{b}}), extracts the curl of the Poynting vector from the second (text{Re}{}) term of Eq. (12),

$$begin{array}{ll}frac{1}{4{omega }^{2}}{text{Re}}left{left({{bf{H}}}^{{boldsymbol{* }}}cdot nabla right){bf{E}}-left({{bf{E}}}^{{boldsymbol{* }}}cdot nabla right){bf{H}}right}&=frac{1}{2{omega }^{2}}nabla times frac{1}{2}{text{Re}}left{{{bf{E}}}^{{boldsymbol{* }}}times{bf{H}}right}&=mathbf{s}_{rm p}end{array}$$
(13)

completing the final step in obtaining the electromagnetic spin decomposition as we present in Eq. (4).

In addition to the monochromatic case, canonical and Poynting spin analogies for the flow of chirality may also be defined in polychromatic fields, expressed using the full time-dependent fields ({boldsymbol{mathcal{E}}}left({bf{r}},tright)) and ({boldsymbol{mathcal{H}}}left({bf{r}},tright)) and stemming from the instantaneous flow of chirality ({boldsymbol{mathcal{F}}}left({bf{r}},tright))29,30,31,42.

$${boldsymbol{mathcal{F}}}=frac{1}{2}left({boldsymbol{mathcal{E}}}times left(nabla times {boldsymbol{mathcal{H}}}right)-boldsymbol{mathcal{H}}times left(nabla times {boldsymbol{mathcal{E}}}right)right)$$
(14)

Chiral flow (boldsymbol{mathcal{F}}), unlike the flux of helicity in time-dependent fields, is defined without the vector potentials (boldsymbol{mathcal{A}}) and (boldsymbol{mathcal{C}}), which avoids gauge complications and contains a direct relation to the instantaneous Poynting vector (boldsymbol{mathcal{P}}=boldsymbol{mathcal{E}}times boldsymbol{mathcal{H}}). For monochromatic light, (boldsymbol{mathcal{F}}) is not time-dependent and becomes proportional to the field’s SAM density ({bf{S}}). Equation (14) is immediately presented in a form that can be broken into two vectors using the two identities detailed in this section, giving the extension of the time-harmonic spin decomposition to time-dependent fields,

$${boldsymbol{mathcal{F}}}=frac{1}{2}left[left({boldsymbol{mathcal{E}}}cdot left(nabla right){boldsymbol{mathcal{H}}}-{boldsymbol{mathcal{H}}}cdot left(nabla right){boldsymbol{mathcal{E}}}right)+nabla times left({boldsymbol{mathcal{E}}}times {boldsymbol{mathcal{H}}}right)right]$$
(15)

The above is a more general decomposition of chirality flow than Eq. (4), and is valid for polychromatic light, at every instant in time. It is worth stressing again that the flow of chirality (boldsymbol{mathcal{F}}) is a different quantity to the flow of helicity (SAM density) for general time-dependent fields. To obtain the helicity flow equivalent of Eq. (15) in the Coulomb gauge, one simply substitutes (boldsymbol{mathcal{E}}to boldsymbol{mathcal{A}}) and (boldsymbol{mathcal{H}}to boldsymbol{mathcal{C}}). Equation (15) assumed a source-free medium, such that (nabla cdot {boldsymbol{mathcal{E}}}=rho /epsilon =0). Dropping this assumption, Eq. (15) can be further generalised by adding a third term, ((rho /2epsilon )boldsymbol{mathcal{H}}), to the decomposition.

4-vector representation

In an attempt to extract new physical insight, we reproduce the above spin decomposition Eq. (4) in 4-vector notation. This notation is particularly elegant and useful for verifying the Lorentz covariance of quantities and equations. We will begin by introducing some common definitions. Electrodynamics is intrinsically relativistic, and several related physical quantities can be expressed via well-known 4-vectors, such as the 4-potential (grouping the scalar and vector potentials) and the 4-current (grouping the charge and current densities). 4-vector notation is well-known to provide an efficient way to formulate Maxwell’s equations. Indeed, the time-dependent Maxwell’s equations (using scripted vectors) in vacuum are a combination of four equations, two scalar equations (nabla cdot {boldsymbol{mathcal{E}}}=0,nabla cdot {boldsymbol{mathcal{H}}}=0) and two vector equations (nablatimes{boldsymbol{mathcal{E}}}=-{mu}_{0}partial {boldsymbol{mathcal{H}}}/partial t), (nablatimes{boldsymbol{mathcal{H}}}={epsilon}_{0}partial {boldsymbol{mathcal{E}}}/partial t). In 4-vector notation, however, we can write Maxwell’s equations as:

$$begin{array}{cc}{partial }_{alpha }{{mathcal{F}}}^{alpha beta }=0 & {epsilon }^{alpha beta gamma lambda }{partial }_{alpha }{{mathcal{F}}}_{beta gamma }=0end{array}$$
(16)

Throughout, the Greek indices (mu ,nu ,alpha ,beta)… run from (0) to (3), where the ({0}^{{th}}) component labels the time direction. Repeated indices are summed over. Note that the distinction between subscript/superscript placement of indices is important, as ({A}_{mu }={eta }_{mu nu }{A}^{nu }) where the so-called Minkowski metric ({eta }_{mu nu }={rm{diag}}left(-1,1,1,1right)). We also use Roman indices (i=1,2,3) to denote the three-dimensional space. In Eq. (16), ({epsilon }^{alpha beta gamma lambda }) is the Levi-Civita symbol43, while ({{mathcal{F}}}_{alpha beta }) is the field strength tensor that conveniently packages the electric and magnetic fields:

$${{mathcal{F}}}_{mu nu }=left(begin{array}{cccc}0 & -{{mathcal{E}}}_{x}/c & -{{mathcal{E}}}_{y}/c & -{{mathcal{E}}}_{z}/c {{mathcal{E}}}_{x}/c & 0 & -{mu }_{0}{{mathcal{H}}}_{z} & {mu }_{0}{{mathcal{H}}}_{y} {{mathcal{E}}}_{y}/c & {mu }_{0}{{mathcal{H}}}_{z} & 0 & -{mu }_{0}{{mathcal{H}}}_{x} {{mathcal{E}}}_{z}/c & -{mu }_{0}{{mathcal{H}}}_{y} & {mu }_{0}{{mathcal{H}}}_{x} & 0end{array}right)$$
(17)

The field strength tensor can be expressed as ({{mathcal{F}}}_{alpha beta }={partial }_{alpha }{{mathcal{A}}}_{beta }-{partial }_{beta }{{mathcal{A}}}_{alpha }), where ({{mathcal{A}}}_{alpha }) is the 4-vector potential, that encases the scalar potential (varphi) and the vector potential (boldsymbol{mathcal{A}}) and ({partial }_{mu }) is a generalisation of gradient:

$${{mathcal{A}}}^{mu }=left(begin{array}{c}{{mathcal{A}}}^{0} {{mathcal{A}}}^{i}end{array}right)=left(begin{array}{c}varphi /c {boldsymbol{mathcal{A}}}end{array}right),{partial }_{mu }=left(begin{array}{c}frac{partial }{cpartial t} {partial }_{i}end{array}right)=left(begin{array}{c}frac{partial }{cpartial t} nabla end{array}right)$$
(18)

Equation (16) is a striking example of the neatness of the tensorial notation: the physical content is the same, but the tensorial notation is clearer and exhibits Lorentz covariance. For time-harmonic fields, these tensors will also have a phasor representation (e.g., ({mathcal{F}}_{mu nu }left({bf{r}},tright)={text{Re}}{F_{mu nu }left(bf{r}right)expleft(-iomega tright)})), including the time components of the gradient, which for phasors becomes ({partial }_{0}=-iomega /c).

We now have the tools to write, using phasors, the time-averaged spin decomposition Eq. (4) in tensorial notation. The first step is to equivalently express the spin, defined in Eq. (3), using spatial indices to replicate the vector operations,

$${S}^{i}=displaystylefrac{1}{4omega }{epsilon }^{ijk}{text{Im}}left{epsilon _{0}E_{j}^{* }E_{k}+{mu }_{0}H_{j}^{* }H_{k}right}$$
(19)

Then, the equivalent of the proposed decomposition in Eq. (4) is:

$$begin{array}{ll}S^i&=underbrace{frac{1}{4{omega }^{2}}{rm{Re}}{{E}_{j}^{ast }{partial }_{i}{H}^{j}-{H}_{j}^{ast }{partial }_{i}{E}^{j}}}_{s_{rm c}^i}&quad+underbrace{frac{1}{2omega^{2}}epsilon^{ijk}{partial }_{j}frac{1}{2}{mathrm{Re}}{epsilon_{klm}E^{*l}H^{m}}}_{{s}_{rm{p}}^i}end{array}$$
(20)

The next step is to find a 4-vector, for which the spatial part would reduce to Eq. (19). It turns out that this 4-vector quantity is exactly the helicity density and flux, a 4-current density associated with conserved helicity44

$${S}^{mu }=frac{1}{4}{mathrm{Re}}left{A_{nu}^{*}G^{nu mu }+C_{nu}^{*}{F}^{nu mu}right}$$
(21)

where ({C}^{mu }={left(psi ,-iomega {mu }_{0}^{-1}left({boldsymbol{nabla }}times {bf{A}}right)right)}^{intercal}) is the magnetic equivalent to the 4-potential ({A}_{mu }), and ({G}_{mu nu }={partial }_{nu }{C}_{mu }-{partial }_{mu }{C}_{nu }) is the corresponding field strength tensor. Note that Eq. (21) is time-averaged and the fields are in their phasor representation. We include the instantaneous variant in the supplementary information. Working in the Coulomb gauge, i.e., the scalar and vector potentials are chosen such that (nabla cdot {bf{A}}=phi =0) and (nabla cdot {bf{C}}=psi =0), and after some algebra (see supplementary Mathematica file45), we find that the spatial part of the 4-vector ({S}^{mu }) reproduces Eq. (19), whilst the temporal part is the cycle-averaged helicity density (h=-{text{Im}}{{{bf{E}}}^{{boldsymbol{* }}},{{cdot }},{bf{H}}}/left(2omega cright)),

$$S^mu=displaystylefrac{1}{4omega}{rm{Im}}left{left(begin{array}{c}-{frac{2}{c}}{mathbf{E}}^*cdot{mathbf{H}}epsilon_{0}{mathbf{E}}^*times{mathbf{E}}+mu_{0}{mathbf{H}}^*times{mathbf{H}}end{array}right)right}=left(begin{array}{c}h {mathbf{S}}end{array}right)$$
(22)

This is consistent with the fact that (h) and ({bf{S}}) are density and flux associated with integrated helicity, the conserved quantity associated with the dual symmetry9,31 (hence (nabla cdot {bf{S}}=0) for free-space monochromatic fields). This four-vector current ({S}^{mu }) has a decomposition equivalent to Eq. (4),

$${S}^{mu }={s}_{{rm{c}}}^{mu }+{s}_{{rm{p}}}^{mu }$$
(23)

where the individual components are given by,

$$begin{array}{ll} & {displaystyle{s}_{{rm{c}}}^{mu }=frac{1}{4}{mathrm{Re}}left{{A}_{nu }^{* }left({partial }^{mu }{C}^{nu }right)-{C}_{nu }^{* }left({partial }^{mu }{A}^{nu }right)right}=frac{1}{4{omega }^{2}}}left(begin{array}{l}-2left(frac{omega }{c}right){rm{Im}}left{{{bf{E}}}^{* }cdot {bf{H}}right} {mathrm{Re}}left{{{bf{E}}}^{* }cdot left({boldsymbol{nabla}}right){bf{H}}-{{bf{H}}}^{*}cdot left({boldsymbol{nabla}}right){bf{E}}right}end{array}right)=left(begin{array}{c}h {{bf{s}}}_{{rm{c}}}end{array}right) & {displaystyle{s}_{{rm{p}}}^{mu }=frac{1}{4}{mathrm{Re}}left{{C}_{nu}^{*}left({partial }^{nu}{A}^{mu }right)-{A}_{nu }^{* }left({partial }^{nu }{C}^{mu }right)right}=frac{1}{4{omega }^{2}}}left(begin{array}{c}0 mathrm{Re}left{left({{bf{H}}}^{* }cdot nabla right){bf{E}}-left({{bf{E}}}^{* }cdot {boldsymbol{nabla }}right){bf{H}}right}end{array}right)=left(begin{array}{c}0 {{bf{s}}}_{{rm{p}}}end{array}right)end{array}$$
(24)

reproducing the spin decomposition in 4-vector notation. Note that the temporal part of the total spin, the helicity density, is carried only by the canonical 4-spin term while the time component of the Poynting 4-spin is zero.

Let us briefly comment on the choice of the Coulomb gauge. As can be seen from the supplementary information, the decomposition of ({S}^{mu }) is independent of the choice of frame and gauge. As already discussed in previous work (e.g., refs. 9,44), while the local density ({S}^{mu }) is not gauge independent, upon integration one gets integrated helicity, which is a gauge-independent quantity that is conserved in any scenario where there is electromagnetic dual symmetry. Furthermore, the integral depends only on the transverse components of the potentials, which makes the Coulomb gauge the most convenient gauge to work with, and in this gauge, the helicity flux (spatial component of ({S}^{mu })) also coincides with the optical spin density. Note, however, that the gauge fixed version of this decomposition is not invariant under Lorentz boosts unless we simultaneously change our potentials to the Coulomb gauge associated with the new boosted frame.

Spin decomposition examples

It is instructive to apply the spin decomposition in Eq. (2) to some simple examples of light: a general evanescent wave and three commonly known focused beams, namely a Gaussian beam, a radial/azimuthal beam ((l=0)), and an (l=1) vortex beam (all linearly polarised in the transverse plane). What differentiates each of these fields are their energy density and phase structures, due to one-dimensionally monotonic, doughnut or Gaussian real-space amplitude profiles, and the presence of OAM, two characteristics which reorient the Poynting vector throughout space and generate different amounts of Poynting spin. Poynting spin can emerge even when the instantaneous vectors (boldsymbol{mathcal{E}}) and (boldsymbol{mathcal{H}}) do not actually rotate, because its counterpart (canonical spin) is able to counterbalance the total SAM density of the field as necessary. Interplay between spin’s canonical and Poynting components creates counter-intuitive effects, such as spin-free chiral interaction forces.

Evanescent wave

The electric and magnetic field phasors of an evanescent wave of arbitrary polarisation, propagating in the (z) direction (({k}_{z} ,>, k)) and decaying along the (x) axis with decay constant (gamma =sqrt{{k}_{z}^{2}-{k}^{2}}), are,

$$begin{array}{ll}{bf{E}}&=left(begin{array}{c}{A}_{p}dfrac{{k}_{z}}{k} {A}_{s} -i{A}_{p}dfrac{gamma }{k}end{array}right){{rm{e}}}^{i{k}_{z}z-gamma x} displaystyle{bf{H}}&=dfrac{1}{eta }left(begin{array}{c}-{A}_{s}dfrac{{k}_{z}}{k} {A}_{p} i{A}_{s}dfrac{gamma }{k}end{array}right){{rm{e}}}^{i{k}_{z}z-gamma x}end{array}$$
(25)

where (eta =sqrt{{mu }_{0}/{epsilon }_{0}}). Choosing complex values for ({A}_{s}) and ({A}_{p}), which are the amplitudes of the evanescent wave’s TE (transverse electric) and TM (transverse magnetic) modes, controls the wave’s polarisation—a circularly polarised wave, for instance, has ({A}_{p}=pm i{A}_{s}). The energy density of the wave decays in the (x) direction and is given by,

$$W={epsilon }_{0}frac{{k}_{z}^{2}}{{k}^{2}}{{rm{e}}}^{-2gamma x}frac{1}{2}left({left|{A}_{s}right|}^{2}+{left|{A}_{p}right|}^{2}right)$$
(26)

Using our formulae, we can calculate the energy-normalised total spin of the evanescent wave, as well as its canonical and Poynting spins,

$$frac{{bf{S}}}{W}=frac{1}{omega {k}_{z}}left[gamma hat{{bf{y}}}+ksigma hat{{bf{z}}}right]$$
(27)
$$frac{{{bf{s}}}_{{rm{c}}}}{W}=frac{1}{omega {k}_{z}}left[frac{{k}_{z}^{2}}{k}sigma hat{{bf{z}}}right]$$
(28)
$$frac{{{bf{s}}}_{{rm{p}}}}{W}=frac{1}{omega {k}_{z}}left[gamma hat{{bf{y}}}-frac{{gamma }^{2}}{k}sigma hat{{bf{z}}}right]$$
(29)

The parameter (sigma =2{text{Im}}{A_{s}A_{p}^{* }}/({|A_{s}|}^{2}+{|A_{p}|}^{2})) is the degree of circular polarisation in the sense of a plane wave ((sigma =pm 1) for circular polarisation, (sigma =0) for linear polarisation). It is well-known that evanescent waves carry transverse spin independently of polarisation27. This property is accounted for by the (hat{{bf{y}}}) component of ({bf{S}}), which is unaffected by the relationship between ({A}_{s}) and ({A}_{p}) and, interestingly, is a product solely of Poynting spin as was discovered in ref. 23. Meanwhile, the evanescent wave acquires a longitudinal spin component if (sigma ,ne, 0), a component which is contributed to by both decomposed spins, and in different amounts. The over-generous canonical spin develops a larger (hat{{bf{z}}}) component than is physical for the total spin of the wave (({{bf{s}}}_{{rm{c}}}cdot hat{{bf{z}}} ,>, {bf{S}}cdot hat{{bf{z}}})). Compensating, the Poynting spin’s (hat{{bf{z}}}) component points backwards to ensure (({{bf{s}}}_{{rm{c}}}+{{bf{s}}}_{{rm{p}}})cdot hat{{bf{z}}}={bf{S}}cdot hat{{bf{z}}}) (this becomes clear after substituting ({gamma }^{2}={k}_{z}^{2}-{k}^{2})).

Decomposing an evanescent wave’s spin reveals a physical distinction between its transverse and longitudinal spin components. Since Poynting spin is responsible for the transverse component of ({bf{S}}), transverse chiral forces felt by an enantiomer arise as it recoils from its own radiation, rather than from a direct field interaction. This is consistent with previous calculations33 and the fact that perpendicular to the wavevector, there is no phase advance to twist the rotating field vectors into helices—the evanescent wave only acquires helicity if it carries longitudinal spin ((h=Wsigma /omega)). Canonical spin, purely longitudinal on the other hand, couples directly to chiral matter to produce a significantly stronger (relatively) preferential force.

Beams

Let us first consider a tightly focussed, linearly polarised Gaussian beam, whose energy density and polarisation in the focal plane is plotted in the top row of Fig. 2a. In the transverse ({xy}) plane, the projection of the electric and magnetic field vectors are linear ((x) and (y) polarised respectively), though due to the tight focusing both ({bf{E}}) and ({bf{H}}) have significant longitudinal (z) components, such that they are elliptically polarised in 3D and contribute a circulation of transverse spin. The Poynting vector magnitude is not constant across the face of the beam and therefore has a non-zero curl in the transverse plane (except in the local maximum at the beam centre). This, the total spin and the decomposed canonical and Poynting spin of the Gaussian beam are plotted in Fig. 2a’s column (descending). Both ({bf{E}}) and ({bf{H}}) fields have (z) components and discretely symmetric polarisation ellipses in 3D (after a rotation of one field by (pi /2)), which appears to reduce the canonical spin of the beam to zero, ({{bf{s}}}_{{rm{c}}}={bf{0}}). All transverse spin developed by focussing of a linearly polarised Gaussian beam stems from Poynting spin, ({bf{S}}={{bf{s}}}_{{rm{p}}}). A similar conclusion could be made for an evanescent wave, making it tempting to argue for a general association between transverse spin and Poynting spin. This notion is not supported by the following example, however.

Fig. 2
figure 2

Spin decomposition of non-paraxial beams. The beams, each of waist (1.5lambda) and separated in columns, are (a) a linearly polarised Gaussian beam (({{mathbf{E}}}^{T}{||}hat{{mathbf{x}}})), (b) an azimuthally polarised (({mathbf{E}}{||}hat{{boldsymbol{phi}}})), (l=0) doughnut beam, and (c) a linearly polarised (({{mathbf{E}}}^{T}{||}hat{{mathbf{x}}})) vortex beam with topological charge (l=1). The top row of (isometric) plots across each sub figure column shows the beam energy density in colour, as well as the Poynting vector and electric (blue) and magnetic (green) polarisation ellipses, which are elliptical due to a significant (z) field component. Subsequent rows are vector plots of each beam’s Poynting vector (mathbf{P}), total spin (mathbf{S}) Eq. (3), canonical spin ({mathbf{s}}_{text{c}}) and Poynting spin ({mathbf{s}}_{text{p}}) from Eq. (4), respectively. White arrows are projections of the corresponding vector into the ({xy}) plane, while the red arrows are projections of the vector onto longitudinal ({yz}) and ({xz}) cut planes. Within each beam, arrows in the three spin decomposition plots are drawn to a consistent scale. Each non-paraxial beam is generated using an angular spectrum integration method46. Defining the 3D vector vortex beam (c) is a difficult problem and the method we used produces a small and physical longitudinal spin (third row of c), contributed by the magnetic field, which would not be present in a (non-physical) perfect paraxial beam

Full size image

An (l=0) azimuthally polarised (in the electric field) beam and its spin decomposition are plotted in Fig. 2b. Considering dual spin (contributions from both ({bf{E}}) and ({bf{H}})), the decomposition properties of this beam are identical for a radially polarised beam, which has an azimuthal magnetic field. Despite tight focussing, the azimuthal electric field vector does not carry a longitudinal component, as is well-known, and is completely linearly polarised in 3D. This means that the beam’s wholly transverse spin is supplied by the magnetic field alone. Combining to give the total spin of the beam, notably both canonical and Poynting spins are non-zero and transverse—this contrasts with the Gaussian beam, whose electric and magnetic fields are both elliptical in 3D. Neither beam considered so far has twisted wavefronts carrying OAM, though due to their azimuthal transverse spin components, both the Gaussian and radial beams possess a certain chiral OAM. A chiral particle could access the spatial structure of the beam’s spin field ({bf{S}}), feeling an azimuthal force (clockwise or anticlockwise depending on the enantiomer) which would cause the particle to orbit the beam centre with no need for helical phase fronts.

Finally, we treat the (x)-polarised (l=1) vortex beam [Fig. 2c], whose OAM leads to the most surprising spin decomposition features of the three. Numerically simulating such a non-paraxial, focussed beam is a challenging task because there is no unanimous 3D definition of a vector vortex beam (only scalar and paraxial beams are well-defined). The beams treated in this section are all electric linearly polarised and generated numerically using an angular spectrum integration technique46 which, by convention, forces the electric field to exactly match the paraxial description of the beam in the focal plane. All unenforced field components (i.e., all magnetic field components and the longitudinal electric component throughout space, and the transverse electric components outside of the focal plane) are subsequently calculated via Maxwell’s equations. While the beams are, therefore, physically valid, some magnetic-biased polarisation features develop in the (x)-polarised vortex beam considered here. In particular, although the electric field remains perfectly linear, the magnetic field alone gains a slight ellipticity in the transverse plane, contributing to a small longitudinal total spin outside of the beam centre as has been shown to exist in focussed linearly polarised beams16. Longitudinal spin remains zero in the singularity, however. Breaking continuous rotational symmetry via (x) polarisation, both the vortex beam’s electric and magnetic fields possess a (z) component from focussing, in contrast to the azimuthal beam for which the electric field is perfectly linearly polarised in 3D. The beam has helical wavefronts so that the average phase gradient (local wavevector) has an azimuthal component, also inherited by the Poynting vector (shown in the second row of Fig. 2c; there is a small transverse component circulating in the ({xy}) plane). The Poynting spin ({{bf{s}}}_{{rm{p}}}propto nabla times {bf{P}}), therefore, acquires a longitudinal component (even at the centre of the beam), in outright defiance of the fact that both transverse electric and magnetic fields are nearly linearly polarised, and zero in the centre, as is visible in the lowest plot of Fig. 2c. To neutralise the longitudinal Poynting spin and suppress the total longitudinal spin, the (z) component of canonical spin ({{bf{s}}}_{{rm{c}}}) must be non-zero too, satisfying ({s}_{{rm{c}}z}=-{s}_{{rm{p}}z}) for ({S}_{z}={s}_{{rm{c}}z}+{s}_{{rm{p}}z}=0) in the centre of the beam, and ({S}_{z}={s}_{{rm{c}}z}+{s}_{{rm{p}}z}approx 0) elsewhere (see the plot second from bottom). In the middle of the beam, this counteracting canonical spin points in the negative (z) direction as a consequence of the beam’s vortex right-handedness, and would switch sign in a left-handed vortex. Returning now to the chiral force equation Eq. (5), we can infer that the beam’s non-zero canonical spin, proportional to chiral momentum, produces a longitudinal chiral force, even though near to the vortex centre the ({bf{E}}) and ({bf{H}}) fields have negligible ellipticity in the transverse plane. Perhaps counter-intuitively, the longitudinal chiral force is strongest in the centre of the beam where the electromagnetic field is virtually zero. Its vortex structure causes the curl of the Poynting vector to be maximal in the central singularity and hence canonical spin, which is equal and opposite along (hat{{bf{z}}}) (ensuring that (({{bf{s}}}_{{rm{c}}}+{{bf{s}}}_{{rm{p}}})cdot hat{{bf{z}}}={bf{S}}cdot hat{{bf{z}}}=0)), must also be at its maximum. This is another example of spatially chiral light coupling to chiral matter which, even under the dipole approximation and in a linearly polarised field, can feel a discriminatory force by absorbing or scattering photons carrying OAM of a certain handedness. Only recently has OAM-dependent chirality, able to couple even to chiral dipoles, been demonstrated in focussed beams, resulting from strong longitudinal field components17,47. We point out, however, that the chiral pressure force which we theorise at the centre of a linearly polarised vortex appears only to require that the Poynting vector has an azimuthal component (i.e., no requirement for longitudinal field components), and should still be present in a paraxial beam.

It is important to distinguish this longitudinal chiral force, appearing in the absence of longitudinal spin, from what may appear at first glance to be similar effects, particularly in linearly polarised evanescent waves, where lateral chiral sorting forces are understood to exist. Although it may be excited via a linearly polarised source still a TE (TM) polarised evanescent wave is elliptically polarised in ({bf{H}}) (({bf{E}})) and carries transverse spin parallel to the direction of the chiral sorting force33. This is in stark contrast to the canonical spin force which we predict in a linearly polarised vortex, pointing in the longitudinal (hat{{bf{z}}}) direction despite there being zero longitudinal spins in the beam centre. Chiral forces present in a linearly polarised evanescent wave and a vortex are, meanwhile, distinct in their interaction with matter. The lateral sorting force calculated by the authors of ref.33 arises from spin recoil, the final term in Eq. (5) (in particular from the contribution of Poynting spin in accordance with23), whereas the canonical spin in Fig. 2c leads to a longitudinal chiral radiation pressure force in the linearly polarised vortex. These chiral force features are summarised in Table 1.

Table 1 Comparison of chiral forces in a linearly polarised evanescent wave and the centre of a linearly polarised vortex and their relation to the spin decomposition
Full size table

Analogies in other wave fields

Kinetic, canonical and spin momenta analogous to the terms in Eq. (2) can be identified more generally in other wave fields48,49,50. This, and the ability to split the SAM density of light into two distinct terms poses another curiosity: what does a spin decomposition of another wave field look like, particularly if, unlike the electromagnetic field, its quanta are not spin-1? Treated in this section are acoustic and gravitational waves, which both depart from light’s spin-1 structure. Despite their increased complexity, linearised theories exist for each of these fields under certain conditions. In a perfect fluid, acoustic waves are linear and oscillate longitudinally (spin-0), while gravitational waves are tensorial in nature (spin-2) and combine linearly when their amplitudes are sufficiently low.

In formulating equivalent expressions for total SAM, canonical and Poynting spins in linearised gravity, we use the Maxwellian representation of gravity in the weak field limit28, i.e., waves propagating over flat spacetime at a large distance from their source. Gravitational waves detected on earth51,52 arrive within this limit, and share some of light’s characteristics; both fields have two polarisation degrees of freedom (in gravity these are the ({h}_{+}) and ({h}_{times }) polarisations) and are massless. Some of light’s strangest behaviour, such as in evanescent fields, have also been predicted with additional properties in gravitational waves53. Note that we take the expressions for the helicity and SAM density that are dual-symmetric (i.e., ({bf{S}}={{bf{S}}}_{{rm{e}}}+{{bf{S}}}_{{rm{m}}}) rather than ({bf{S}}=2{{bf{S}}}_{{rm{e}}})), the same as in the case of electromagnetism, despite the fact that in the case of gravity we have no experimental evidence favouring its physical relevance over the asymmetric version54. However, if we take the point of view of ref.9 that only the integrated helicity is a physically meaningful quantity, the two will be equivalent.

Table 2 summarises the differences and similarities between linearised acoustics, electromagnetism, and linearised gravity with a focus on the spin-related quantities. An acoustic wave field can be described by a scalar pressure field (P) and a vector velocity field ({bf{v}}) which, in linearised acoustic theory, share a Maxwell-like relation49,50. The derivation of the decomposition for the acoustic field is given in the supplementary material. A gravitational wave can be described using a metric perturbation ({h}_{{ij}}), which can be thought of as components of a three-by-three symmetric matrix. If we consider only (i)-th row or column of this matrix as a vector potential, ({{bf{A}}}^{i}={h}^{{ij}}{hat{{bf{e}}}}_{j}), then the Maxwellian representation of gravity can be written in vector notation; the derivation of the spin decomposition then becomes simply Eq. (12). The acoustic helicity density is zero, which is a feature of spin-0 fields, but interestingly the canonical spin, which we showed in Eq. (24) to be related to helicity, is also zero. One can see that electromagnetic and gravitational waves are spin-1 and spin-2 respectively by taking a circularly polarised wave. The helicity density of a circularly polarised wave will be one (two) times the energy density and canonical spin will be one (two) times the canonical momentum for an electromagnetic (gravitational) wave.

Table 2 Comparison of time-averaged energy, momentum and spin densities between time-harmonic waves in theories of linearised acoustics, electromagnetism and linearised gravity
Full size table

Discussion

We have discussed a known but overlooked decomposition which exists for the spin angular momentum density of light, in a similar way to the Poynting vector, which can be split into orbital and spin currents. Spin is decomposed into two terms, which we call the canonical spin ({{bf{s}}}_{{rm{c}}}) and Poynting spin ({{bf{s}}}_{{rm{p}}}). We have further expressed the decomposition for time-varying fields, as well as in four-vector notation. The canonical spin is proportional to chiral linear momentum which can be transferred to chiral enantiomers in a positive or negative direction, depending on the enantiomer handedness. The resulting force is a chiral analogy to the (achiral) force due to radiation pressure—both are relatively strong forces and can act on matter in the dipole approximation.

The mechanism for preferential photon absorption or scattering by chiral matter is normally associated with the photon’s spin state (i.e., circular dichroism), and less so to the sign of its OAM (although these interactions have received recent attention). We emphasised, however, that light is capable of exerting a chiral force that does not directly depend on the total SAM of light, rather, its canonical spin, which can take a non-zero value even in linearly polarised fields. This is due to the second term in the spin decomposition, the Poynting spin, which by definition (being the curl of the Poynting vector) depends strongly on light’s OAM: optical vortices twisting energy flow around the beam axis naturally generates non-zero Poynting spin, regardless of polarisation. We showed that in a linearly polarised optical vortex, the longitudinal components of canonical spin and Poynting spin must be non-zero (due to OAM) and equal and opposite to ensure total longitudinal spin is zero, resulting in a chiral force along the beam axis in the absence of longitudinal SAM. We believe it is the intuition provided by Eq. (4) and the relation it shows between two components of total spin and chiral interaction and recoil forces that is one of the main values of the spin decomposition (which also simplifies expressions for optical forces, showing clear analogies between canonical spin and canonical momentum, Poynting spin and spin momentum37).

Interestingly, this OAM-dependent chiral force is strongest in the centre of the vortex, where the electromagnetic energy density is minimal, along with achiral forces such as radiation pressure and gradient force. Given that it is possible to trap atoms and molecules (using a blue-detuned wavelength with respect to a strong resonance) in dark spots55,56,57, it should be possible to directly measure a chiral force due to canonical spin by structuring a field in such a way as to isolate canonical spin in a certain direction. One plausible way to achieve this would be to use counter-propagating vortices of opposite topological charge, in whose centre chiral particles could be trapped and respond to the longitudinal component of canonical spin. Chiral forces due to Poynting spin could be readily measured from a linearly polarised evanescent wave using a modified nano-cantilever which, as pointed out by ref.33, has already been used to measure non-chiral lateral forces in evanescent waves58.

The ability to decompose spin appears to be a general property of wave fields, and as we demonstrated, is straightforward to perform in linearised acoustics and gravity, whose quanta are considered spin-0 and spin-2 respectively. This fact could be of significant interest to a broader community, beyond optics.

Although the vast majority of this work is cast in the context of monochromatic light propagating in free-space as is necessary to obtain Eq. (4), the ‘monochromatic’ and ‘free-space’ criteria are not required to maintain that SAM density is composed of two separate contributions. Time-varying light’s chiral flow can be split into canonical and Poynting-equivalent terms as we demonstrated in Eq. (15). It is also intriguing to ask, in view of the decades-long debate of which of the kinetic or canonical momenta is the correct measure of linear momentum in macroscopic media, whether one can find analogy to the Abraham-Minkowski dilemma in the relationship between total spin and canonical spin.

Materials and methods

The vast majority of our discussion has been theoretical. Numerical data in Fig. 2 for each of the non-paraxial Gaussian, azimuthal and vortex beams were generated using the software package46. This software simulates non-paraxial beams simply by superposing plane waves according to an analytic function for the beam’s angular spectrum. The angular spectrum function is calculated by the software package in such a way as to force the electric field vector to match the paraxial description of the beam in the focal plane (at (z=0)). The longitudinal electric field component throughout the simulation domain, the transverse electric field outside of the focal plane and the entire magnetic field are calculated from the subsequent plane wave superposition. This yields a completely physical beam, a full solution to Maxwell’s equations because of their linearity.