Optical momentum and angular momentum in complex media: From the Abraham-Minkowski debate to unusual properties of surface plasmon-polaritons

We examine the momentum and angular-momentum (AM) properties of monochromatic optical fields in dispersive and inhomogeneous isotropic media, using the Abraham- and Minkowski-type approaches, as well as the kinetic (Poynting-like) and canonical (with separate spin and orbital degrees of freedom) pictures. While the kinetic Abraham-Poynting momentum describes the energy flux and the group velocity of the wave, the Minkowski-type quantities, with proper dispersion corrections, describe the actual momentum and angular momentum carried by the wave. The kinetic Minkowski-type momentum and AM densities agree with phenomenological results derived by Philbin. Using the canonical spin-orbital decomposition, previously used for free-space fields, we find the corresponding canonical momentum, spin and orbital AM of light in a dispersive inhomogeneous medium. These acquire a very natural form analogous to the Brillouin energy density and are valid for arbitrary structured fields. The general theory is applied to a non-trivial example of a surface plasmon-polariton (SPP) wave at a metal-vacuum interface. We show that the integral momentum of the SPP per particle corresponds to the SPP wave vector, and hence exceeds the momentum of a photon in the vacuum. We also provide the first accurate calculation of the transverse spin and orbital AM of the SPP. While the intrinsic orbital AM vanishes, the transverse spin can change its sign depending on the SPP frequency. Importantly, we present both macroscopic and microscopic calculations, thereby proving the validity of the general phenomenological results. The microscopic theory also predicts a transverse magnetization in the metal (i.e., a magnetic moment for the SPP) as well as the corresponding direct magnetization current, which explains the difference between the Abraham and Minkowski momenta.


Abraham and Minkowski momenta
The characterization of the momentum and angular momentum (AM) of light in continuous media is a long-standing problem, with the Abraham-Minkowski discussion in its center; see [1][2][3][4][5] for reviews. Although recently there were several works claiming the "resolution" of the Abraham-Minkowski controversy [5][6][7], debates on various aspects of optical momentum in media are still continuing. Naturally, different momenta can manifest in different types of problems or experiments. For example, one can investigate optical forces acting on the medium or on a small material probe in the medium [7][8][9][10]. On the other hand, one may wonder about the momentum carried by the wave per se and characterizing wave parameters such as the velocity of its propagation (phase or group) or wave vector [11][12][13][14][15][16][17][18][19]. In this work we mostly stick with the second approach. Throughout this paper we consider monochromatic waves with fixed frequency ω , lossless isotropic media described by permittivity and permeability ε and µ (which can depend on ω in the dispersive case), and cycle-averaged dynamical properties (energy, momentum, and AM) of waves.
To recall the basics of the problem, we start with the Poynting momentum density of monochromatic light in free space [20][21][22]: ) . (1.1) Hereafter, we use Gaussian units 1 with g = 8πω ( ) −1 , k 0 = ω / c , and all free-space quantities are marked by the subscript "0". In a non-dispersive medium, the Abraham and Minkowski momentum densities are given by [1][2][3][4][5] P A = P 0 , P M = ε µP 0 . (1.2) These two momenta are often interpreted as "kinetic" and "canonical" momenta of light, respectively [3][4][5]7,11,13,18]. In dispersive media, the Abraham momentum preserves its form, while the "canonical" momentum should be modified with dispersion-related terms [12][13][14]16,17]: (1.3) and the "naïve" Minkowski momentum (1.2) does not make physical sense. Hereafter, we mark by tilde all quantities modified by the presence of dispersion. For simplicity, we will refer to the momentum (1.3) as to the properly-modified Minkowski momentum in a dispersive medium. For the simplest optical fields -plane waves -in a transparent medium, the Abraham and Minkowski momenta "per photon" are reduced to the following simple form [3][4][5]7,[11][12][13][14][16][17][18]: where n p = εµ and n g = n p + ω dn p / dω are the phase and group refractive indices of the medium, respectively, v g = ∂ω / ∂k is the group velocity in the medium ( v g = c / n g ), and k is the wave vector in the medium, with magnitude k = n p k 0 . Equations (1.4) shed light on the physical meaning of the Abraham and Minkowski momenta, which can be associated with the group velocity and wave vector in the medium. This correspondence is very general [11,18]. For example, even for inhomogeneous waves in nontransparent inhomogeneous media, such as surface plasmon-polaritons (SPPs) at metal-vacuum interfaces [23], the group velocity is still determined by the integral value of the Poynting (=Abraham) momentum [24]. At the same time, the conservation of the wave momentum and momentum matching in various resonant problems involve Minkowski momentum, i.e., the wave vector in the medium. The most known example is Snell's law in the light refraction at planar interfaces [20,22]. Furthermore, more subtle spin and orbital Hall effects (transverse beam shifts) in the refraction at planar interfaces are intimately related to the conservation of the corresponding Minkowski AM based on the same wavevector k [25][26][27][28][29][30]. 1 For the conversion between Gaussian and SI units see Appendix 4 in [20]. This momentum describes only the orbital part of the AM, L 0 , while the spin part is provided by an independent intrinsic quantity S 0 [35,41,52,54,55,57]: (1.7) Recently, it was shown that the canonical quantities (1.6) and (1.7) are much more suitable for the description of the momentum and AM properties of free-space light than the kinetic Poynting characteristics (1.1) and (1.5). In particular, the optical force and torque on a small electricdipole particle or an atom are given by the electric parts of the canonical momentum (1.6) and spin AM (1.7), respectively [35,41,46,47,[58][59][60]. This makes canonical quantities directly measurable and immediately explaining numerous experiments involving spin/orbital AM [35,41,[61][62][63] and structured light fields [31][32][33][34][35]46,47,64]. Moreover, using canonical formalism enabled prediction and description of unusual phenomena, such as unusual transverse spin AM in evanescent and other structured fields [35,[41][42][43]46,[65][66][67][68] and super-momentum transfer higher than !k 0 per photon [31][32][33][34][35]. The canonical momentum density (1.6) can be written as a local expectation value of the quantum-mechanical momentum operator p = −i∇ , and hence can be associated with the local phase gradient or local wave vector k loc of the field [33,56]. This elucidates its canonical character, akin to Minkowski momentum (1.4). However, in contrast to Eq. (1.4), valid for a single plane wave, canonical momentum (1.6) describes the local phase gradient in an arbitrary structured field, which can consist of multiple plane waves propagating in different directions. In turn, the canonical spin AM density (1.7) describes the local ellipticity of the 3D polarization of an arbitrary structured field.
The relation between the kinetic (Poynting) momentum P 0 and canonical momentum P 0 in free space is given by the spin-orbital momentum decomposition [35,41,47,55,56]: Here the canonical momentum P 0 describes the orbital part (which determines the orbital AM L 0 ), while the spin momentum P 0 S is related to the spin AM S 0 via the nonlocal relation ∫ dV , valid for any localized fields vanishing at infinity. Importantly, the spin momentum vanishes for plane waves and does not contribute to the integral (expectation) value of the wave momentum for localized fields, so that the integral kinetic and canonical momenta coincide: Here, ... = ...
∫ dV , and hereafter it denotes suitable spatial integrals for localized fields.
In terms of relativistic field theory, the canonical momentum and spin densities (1.6) and (1.7) originate from the canonical energy-momentum and AM tensors, which are directly obtained from Noether's theorem applied to the electromagnetic field Lagrangian [36,44,45,54,55]. Two points should be emphasized here. First, the original form of these canonical tensors involves the gauge-dependent electromagnetic vector potential A . The standard procedure in this case is to consider only the "transverse" (i.e., divergence-free) gaugeinvariant part of this potential, A ⊥ , which for monochromatic fields is expressed via the wave electric field A ⊥ = −i c / ω ( )E [44,45,47,[48][49][50]52,54,55,60]. Second, the standard electromagnetic-field Lagrangian is not dual-symmetric with respect to the electric and magnetic fields. Due to this, it results in double electric-field parts of quantities (1.6) and (1.7), with no magnetic-field parts [36,44,47,60]. However, an alternative Lagrangian formalism, dualsymmetrized between electric and magnetic contributions [55,69,70], produces the symmetric quantities (1.6) and (1.7), more natural for free Maxwell fields [35,41,[54][55][56][57]59]. In this paper we employ the dual-symmetric formalism [55] and show that it is more consistent with the canonical optical properties in media than the dual-asymmetric (electric-biased) approach.
While the symmetrized (kinetic) energy-momentum tensor contains only the Poynting momentum density P 0 , the canonical energy-momentum tensor is non-symmetric, and contains both the Poynting vector P 0 , acting as the energy flux density, and the canonical momentum density P 0 . Remarkably, considering the energy-momentum tensor for electromagnetic waves in a medium, Dewar [11] and later Dodin and Fisch [18] found that the electromagnetic energymomentum tensor in a (non-dispersive) medium can be modified to the Minkowski form, where the Poynting energy flux and canonical momentum are substituted by the Abraham (Poynting) and Minkowski momenta, P A and P M , respectively. Schematically, these different forms of the energy-momentum tensors (EMTs) in free space and in media can be presented as follows:  Dewar (1977) [11] ! " ########## , (1.10) This provides a qualitative link between the Abraham-Minkowski and kinetic-canonical (in the relativistic field-theory sense) dilemmas.
Summarizing the above considerations, one should associate the Poynting-Abraham quantities with the energy flux and group velocity of the wave-packet propagation, while the canonical and Minkowski quantities are related to the momentum density carried by the wave and its wave-vector characteristics. At the same time, we emphasize that the kinetic-canonical dilemma between P 0 and P 0 in the sense of relativistic field theory originates from the separation of the spin and orbital degrees of freedom, while the Abraham-Minkowski dilemma between P A and P M is related to the separation of the medium and field contributions to the momentum. Therefore, one can consider the spin-orbital separation in both Abraham and Minkowski momenta in a medium, as well as Abraham and Minkowski forms of the kinetic and canonical (orbital) momenta in a medium, i.e., four types of momenta in the medium. In what follows, we use the "kinetic" and "canonical" characteristics in the field-theory sense of the spinorbital separation, also explicitly indicating the Abraham-and Minkowski-type quantities.

About this work
Here we aim to provide a complete Abraham-Minkowski and kinetic-canonical picture of optical momentum and AM in dispersive and inhomogeneous (but isotropic and lossless) media. For the reader's convenience, we summarize all the quantities under discussion in Table I, indicating their forms in free space, dispersion-free, and dispersive media. Our main emphasis in this work is on the Minkowski-type quantities, because these correspond to the actual wave momentum, spin, and AM in the medium, in contrast to the Abraham-type energy flux properties.
The paper is organized as follows. In this introductory Section 1 we provided a general overview of the problem. In Section 2, we discuss the general momentum and AM expressions listed in Table I and their properties. Then, in Section 3, we consider an explicit example of a SPP wave at a metal-vacuum interface.
This example (which to the best of our knowledge has never been considered in the Abraham-Minkowski context) provides a perfect test for the momentum and AM properties of structured optical fields in dispersive inhomogeneous media. The multiple advantages of the SPP system are as follows: (i) SPP waves are well studied and readily achievable experimentally; (ii) Even a single SPP wave is a structured field, for which the simplified plane-wave Eqs. (1.4) are not applicable; (iii) SPPs exist at interfaces, i.e., in essentially inhomogeneous media. (iv) SPPs exhibit nontrivial momentum and AM properties, including super-momentum [31][32][33][34][35] and transverse spin AM [35,[41][42][43]46,[65][66][67][68]; (v) Dispersion of the metal is crucial for the SPP properties.
Thus, SPPs provide both an accessible and highly nontrivial system to study optical momentum and AM.
In Section 3 we apply the general expressions from Table I to calculate the momentum and AM properties of SPPs. In particular, we present the first accurate calculations of the canonical and Minkowski momenta, as well as the transverse spin and orbital AM of a SPP. Notably, even the integral canonical or Minkowski momentum of SPP exceeds !k 0 per particle, and this offers the first example of the integral super-momentum (previously known only locally). Moreover, the intrinsic orbital AM of the SPP vanishes, whereas the integral transverse spin AM can change its sign depending on parameters.  [16,17], Table I. Four possible pictures of the optical momentum and angular momentum (AM) densities in free space, non-dispersive isotropic media, and dispersive (generally inhomogeneous) media. The Abraham-and Minkowski-type, kinetic and canonical (spin-orbital) quantities are shown. In all cases, the kinetic Abraham-Poynting momentum density P A = P 0 describes the energy flux and group velocity of the wave, whereas the canonical Minkowski-type momentum and spin densities ! P M and ! S M provide a clear and self-consistent picture of the momentum and angular momentum carried by the wave. In turn, the kinetic-Minkowski and canonical-Abraham quantities have less natural forms with cumbersome dispersive and gradient corrections (indicated as {dispers.} and {grad.} here).
Importantly, in Section 4, we provide microscopic calculations of the momentum and AM in SPPs. Taking into account both microscopic electromagnetic fields and the motion of free electrons in the metal, we obtain the Minkowski and canonical quantities previously introduced using macroscopic phenomenological considerations. This validates the use of these quantities for structured optical fields in dispersive and inhomogeneous media. Moreover, the microscopic theory predicts a transverse magnetization in the metal (i.e., a magnetic moment for the SPP) as well as the corresponding direct magnetization current, which corresponds to the difference between the Abraham and Minkowski-type momenta.
Finally, in Section 5, we briefly discuss issues related to the dual symmetry between electric-and magnetic-field contributions. We show that while integral electric and magnetic conributions to the momentum and spin are equal for localized fields in free space [52,55], this is not the case for localized fields in media. Most importantly, we find that the microscopic calculations of Section 4 are only compatible with the dual-symmetric (rather than electricbiased) forms of the canonical quantities.
We should also briefly mention preceding works, which considered some of the above aspects of SPPs. First, Nakamura [71] performed microscopic calculations of the transverse AM of SPPs. Although results of that work are erroneous in several aspects (calculation errors, mixing of the spin and orbital AM, etc.), its methodology inspired us to perform microscopic calculations presented in Section 4. Second, Kim and Wang aimed to calculated Abraham and "naïve" Minkowski (without dispersion terms) versions of the transverse spin AM in SPPs [72][73][74]. However, their results are also misleading. First, the definitions in [72][73][74] do not describe the Abraham-type spin, which was properly defined and calculated in [35,65], and which corresponds to an energy-flux property rather than the actual spin AM carried by the wave. Second, the "naïve" Minkowski expressions are not applicable to waves in dispersive media and lead to erroneous results. Thus, the first accurate calculation of the transverse spin and other "canonical" characteristics of SPPs are provided in the present work.

Momentum and angular momentum of light in dispersive inhomogeneous media
Throughout this paper we consider monochromatic electric and magnetic fields: (2.1) We now consider an isotropic lossless dispersive and inhomogeneous medium, which is characterized by the real frequency-dependent permittivity ε ω ,r ( ) and permeability µ ω ,r ( ) . In this case, the complex field amplitudes satisfy stationary Maxwell equations Note that these source-free equations are used in the decomposition of the Poynting momentum density into canonical and spin parts, Eq. (1.8).
The energy density of a monochromatic optical field in such a medium is described by the well-known Brillouin expression [20,22]: Describing the optical momentum density in a medium is a more sophisticated problem. On the one hand, the Abraham momentum P A preserves its Poynting-vector form in the medium, Eq. (1.2). By analogy with the canonical decomposition (1.8), one can decompose it into orbital and spin parts, This spin-orbital decomposition was introduced in [65] and was used in [35,46,80] because of its convenience in homogeneous media. However, in inhomogeneous media, the canonical momentum density (2.5) acquires cumbersome gradient terms. Moreover, the physical interpretation of the quantities (2.5) and (2.6) is not quite clear. Indeed, as we discussed above, the Abraham momentum P A should be associated with the energy flux density and group velocity, Eq. (1.10), rather than with the wave momentum density. Therefore, the Abraham-type quantities (2.5) and (2.6) correspond to the orbital and spin parts of the energy flux density, but cannot be regarded as canonical momentum and spin densities in the wave. In addition, at interfaces between media, the Abraham-type spin density (2.6) is discontinuous, and the corresponding gradient terms in Eq. (2.5) produce singular delta-function contributions to the canonical and spin momentum densities P A and P A S [65]. This makes the Abraham-type spinorbital decomposition imperfect. Note also that, similarly to the free-space Eqs. (1.8) and (1.9), the solenoidal spin part of the energy flux does not contribute to the plane-wave and integral characteristics: Therefore, in some plane-wave or integral calculations it could be more convenient to use Eq. (2.5) as the energy flux density.
To describe physically-meaningful momentum and AM densities in the optical field, one should use the Minkowski momentum. Its simple form (1.2) is not valid in the presence of dispersion, and several works discussed modifications of the Minkowski-type wave momentum in a dispersive medium [12][13][14]16,17]. The most general expression, suitable for structured wave fields was derived by Philbin [16,17] using the phenomenological Lagrangian formalism and Noether's theorem: Here, the first term is the "naïve" Minkowski momentum (1.2), while the second term describes the dispersion-related correction. For a plane wave in a transparent homogeneous medium, the Abraham and modified-Minkowski momenta P A and ! P M yield simplified Eqs. (1.4) [3][4][5]7,[11][12][13][14][16][17][18].
Since the Minkowski-type momentum (2.8) represents "canonical" wave-vector properties of the wave, it makes sense to find the spin-orbital decomposition, similar to Eq. (1.8), and introduce the corresponding canonical spin and orbital properties. In doing so, we apply the standard Poynting-vector decomposition (1.8) to the first (Minkowski) term in Eq. (2.8) and add the second dispersive term to the orbital part (because of its natural orbital form). This results in ! P M = ! P M + P M S : (2.10) Notably, the canonical (orbital) momentum density (2.9) has a nice form similar to the free-space momentum (1.6) with the ! ε and ! µ multipliers, exactly as in the Brillouin energy density (2.3).
Furthermore, the momentum (2.9) is free of cumbersome gradient terms, present in the canonical Abraham-type momentum (2.5). However, the quantity S M in Eq. (2.10) is the "naïve" Minkowski spin AM density, which lacks dispersive corrections. As we show below, this is not the canonical spin AM density of the wave. Notably, in the SPP example considered below, the quantity S M is continuous at the interface, and therefore the canonical and spin parts of the modified Minkowski momentum density, ! P M and P M S , are free of delta-function singularities.
This makes the Minkowski-type spin-orbital decomposition more appealing than the Abraham one. Akin to Eqs. (1.8), (1.9), and (2.7), the solenoidal spin momentum (2.10) vanishes for plane waves and does not contribute to the integral momentum values. Therefore, the integral values of the kinetic and canonical Minkowski-type momenta (2.8) and (2.9) coincide for localized fields: Thus, one can use either of these momenta in calculations of the integral or plane-wave properties.
To determine the canonical spin and orbital AM in a dispersive medium, we start with the kinetic (total) Minkowski-type AM. Again, the Minkowski-type analogue of Eq. (1.5) in a dispersive medium was found by Philbin and Allanson [17]: Importantly, the AM density (2.12) breaks the simple relation (1.5) between the kinetic momentum and AM densities, ! J M ≠ r × ! P M , and contains a dispersion-related correction of the spin-like local form. Substituting ! P M = ! P M + P M S with Eqs. (2.9) and (2.10) into Eq. (2.12) and using the nonlocal relation between the spin momentum and spin AM, we derive the canonical orbital and spin AM densities in a dispersive medium: Here the dispersion terms from Eq. (2.12) correct the "naïve" Minkowski spin density S M . Remarkably, Eqs. (2.13) have a very nice form, similar to the free-space Eqs. (1.7), but now with the same ! ε and ! µ multipliers as in the Brillouin energy density (2.3) and canonical momentum density (2.9). The integral values of the kinetic and canonical AM (2.12) and (2.13) coincide for localized fields: (2.14) Equations (2.3), (2.9), and (2.13) constitute a set of canonical characteristics of a monochromatic light in a dispersive and inhomogeneous medium. Importantly, these can be written as local expectation values of quantum-mechanical energy (ω ), momentum ( p = −i∇ ), and spin-1 ( Ŝ ) operators [35,[54][55][56]59]: (2.15) using the following wave-function: Exactly the same formalism for electromagnetic bi-linear forms (including Berry connection and other topological characteristics) in dispersive media was recently suggested in works by Silveirinha [81][82][83]. Thus, the above equations bring together approaches developed by (i) Philbin (kinetic Minkowski-type momentum and AM in dispersive media) [16,17] .15) suggests that the canonical form of the Minkowski-type momentum and AM densities is more suitable for the description of the optical momentum and AM than the previously-used kinetic Minkowski-type approach, Eqs. (2.8) and (2.12), and the Abraham-type quantities (2.5) and (2.6). For example, consider a polarized plane wave propagating in a homogeneous dispersive medium. All field components have the same phase factor exp i k ⋅ r ( ) , the electric and magnetic field amplitudes are related by E 2 / µ = H 2 / ε , whereas the ellipticity of the polarization can be characterized by the helicity σ , such that Im E * × E ( ) = σ E 2 k (and a similar equation for the magnetic field), where k = k / k characterizes the direction of the wave propagation. Using these simple properties, from Eqs. (2.3), (2.5), (2.6), (2.9) and (2.13), we readily obtain the ratios of the canonical Abraham-and Minkowski-type momentum and spin densities to the energy density of the wave: The first two of these equations exactly correspond to Eqs. (1.4), whereas the other two equations provide their counterparts for the spin AM (cf. [17,84]). In this manner, the Minkowski-type momentum ! P M and spin AM ! S M correspond to the values !k and !σ k per photon, as one would expect for photons, the Abraham-type momentum P A determines the group velocity (1.4), while the Abraham-type spin S A does not have a clear physical meaning.
It is important to note that it is the Minkowski-type wave momentum and AM that are conserved in media with the corresponding translational and rotational symmetries. First, this follows in the most general form from the results of [16,17], where the kinetic quantities (2.8) and (2.12) were derived from Noether's theorem. In view of Eqs. (2.11) and (2.14), this is also true for the canonical Minkowski-type quantities (2.9) and (2.13). Second, in the plane-wave Eqs. (2.17), the Minkowski-type momentum and spin exactly correspond to the tangentmomentum and normal-AM conservation laws in the wave refraction at an interface between two media [20,22,[25][26][27][28][29][30] (Snell's law and optical beam shifts).
In this Section we presented a macroscopic phenomenological introduction of these quantities. Below, considering SPPs at the vacuum-metal interface, we show that this macroscopic model is in exact agreement with microscopic calculations taking into account separate electron and field contributions. It should be also noted that in the absence of dispersion, ! ε = ε , ! µ = µ , and both kinetic and canonical characteristics discussed in this Section acquire simplified Minkowski forms, shown in Table I.

SPP fields and parameters.
We now consider an explicit example of a structured optical field in a dispersive and inhomogeneous medium: a surface plasmon-polariton (SPP) at the metal-vacuum interface [23]. The geometry of the problem is shown in Fig. 1(a), where the interface is the x = 0 plane with the vacuum in the x > 0 half-space (medium 1) and metal in the x < 0 half-space (medium 2), whereas the SPP wave propagates along the z -axis with the wavevector k p = k p z (hereafter, x , y , and z denote the unit vectors of the corresponding axes). The permittivity and permeability of the metal are given using the standard plasma model [23]: Here, ω p 2 = 4πn 0 e 2 / m is the plasma frequency, where n 0 is the volume density of free electrons in the metal, e < 0 is the electron charge, and m is the electron mass. Thus, the metal is a dispersive medium with ! ε = 1+ ω p 2 / ω 2 ≠ ε . SPPs are electromagnetic surface TM waves that exist at frequencies ω < ω p / 2 where ε < −1 [23]. The electric and magnetic fields of a single SPP wave can be written as [23,24,65]: where A is the field amplitude, whereas the wave number and spatial decay constants of the SPP field are The vacuum part of the SPP fields

Energy, energy flux, and momentum of SPPs.
Substituting Eqs. (3.2) and (3.3) into Eq. (2.3), we obtain the energy density distribution in the SPP field: Note that the distribution (3.4) is discontinuous at the interface x = 0 . Since SPP field is localized along the x coordinate, we can also calculate the integral "expectation value" of the SPP energy integrating ! W over x : In what follows, we express the integral momentum and AM characteristics of SPPs with respect to the energy (3.5), in order to highlight their values "per plasmon". From the energy-density distribution, we can also find the position of the center of energy along the x -axis: It follows from here that for ω < ω p /1.57 the energy centroid is located in the vacuum: x > 0 , while for ω p /1.57 < ω < ω p / 2 it moves into the metal: x < 0 . Dependences of the integral energy (3.5) and center-of-energy position (3.6) on the SPP frequency are shown in Fig. 2 Note that this flux is negative (backward) inside the metal, producing a vortex-like energy circulation in SPP wave packets [24,65,85]. Nonetheless, the integral energy flux P A is positive, and, in ratio to the energy (3.5), it determines the group velocity of SPPs [24,86]: which corresponds to the subluminal propagation of a SPP wave packet. This also means that the integral Abraham momentum is smaller than !k 0 "per plasmon". Although the plane-wave homogenous-medium equations (1.4) and (2.17) are not directly applicable to the structured SPP wave at the metal-vacuum interface, the group velocity (3.8) can be written in the form v g = c / n g eff , where we introduced an effective phase refractive index for SPPs, n p eff = k p / k 0 > 1, and the corresponding effective group index n g eff = n p eff + ω dn p eff / dω > 1 . This shows that the relations (1.4) between the Abraham energy flux, group velocity, and refractive indices are rather general and can be extended, using integral expectation values, to localized states in inhomogeneous media. The frequency dependence of the SPP group velocity (3.8) is depicted in Fig. 3(a).
The orbital and spin parts of the energy flux (3.7), P A and P A S = 1 2 ∇ × S A , Eqs. (2.5) and (2.6), have been analyzed for SPPs in [65], and we do not reproduce these here. We just recall that these parts have singular delta-function contributions at the interface x = 0 due to the gradient terms in Eq. (2.5) and discontinuity of the Abraham spin (2.6) S A (shown below). These singular contributions are crucial to satisfy Eq. (2.7). We now calculate the Minkowski-type momentum for the SPP fields. These have a more natural form in the canonical approach. Indeed, using Eq. (2.9), we readily obtain for the field (3.2) with the common exp ik p z ( ) phase factor: Note that the momentum density ! P M is positive both in the vacuum and in metal, in contrast to the energy flux (3.7). Moreover, it does not have delta-function singularities at the interface, in contrast to the canonical Abraham-type momentum P A [65]. Fig. 3. The group velocity (3.8) and integral Minkowski-type momentum (3.9) of a SPP versus frequency. While the group velocity is always sub-luminal, v g < c , the integral momentum of SPP corresponds to the super-momentum !k p > !k 0 "per plasmon".
Assuming that the energy is quantized as !ω "per plasmon", Eqs. (3.9) mean that the SPP carries super-momentum !k p > !k 0 per plasmon. Remarkably, so far the super-momentum was described only locally, in evanescent waves or near optical vortices [31][32][33][34][35]56,87,88]. At the same time, the integral momentum for any localized wave field in free space is always less than !k 0 per photon [33]. Equations (3.9) show an example of the integral super-momentum. Of course, one can say that a canonical momentum higher than !k 0 per photon appears in any medium with phase refractive index n p > 1 . Equations (3.9) can also be written in the form of Eqs. (1.4) and (2.17) with k = k p z and effective refractive index n p eff = k p / k 0 > 1. However, in the case of SPP waves, this effective refractive index originates from the inhomogeneous evanescent character of surface waves rather than from the high permittivity of the medium. Indeed, in the limit ω → ω p / 2 , we have n p eff → ∞ , while ε → −1. Note that although in some works the super-momentum was interpreted as a "superluminal group velocity" [33,88], our present analysis shows that it should rather be considered as a pure momentum property, while the group velocity is determined by the Poynting-Abraham energy flux and is always subluminal. The frequency dependence of the ratio of the SPP momentum and energy, Eq. (3.9), is shown in Fig. 3(b). According to Eq. (2.11), the kinetic form of the Minkowski-type momentum, Eq. (2.8), has the same integral value (3.9), but its local density does not exhibit the nice proportionality to the energy density as in Eq. (3.9): Moreover, this momentum density coincides with the Poynting-Abraham one in the vacuum. Therefore, for x > 0 , c ! P M / ! W = c P A / ! W = k 0 / k p < 1 , and the kinetic Minkowski-type momentum cannot explain the local super-momentum density in the vacuum evanescent field, c ! P M / ! W = k p / k 0 > 1, which is described by the canonical momentum (3.9) and is observed experimentally [31][32][33]. Comparing Eqs. (3.8)-(3.10), together with the singular character of the canonical Abraham-type momentum P A [65], confirms that the canonical picture is more natural for the description of the Minkowski-type wave momentum, while the kinetic approach is more suitable for the characterization of the Abraham-type energy fluxes, see Table I.

Spin and orbital AM of SPPs.
We are now in the position to determine the spin and orbital AM of SPPs. Akin to the momentum of SPPs, these should be described using the Minkowski-type canonical picture. However, for completeness and comparison with other approaches, we first calculate the Abraham-type spin density (2.6) and its integral value. With the SPP fields (3.2) and (3.3), we obtain This is the transverse helicity-independent spin, first described in [65] and now attracting considerable attention [35,[41][42][43]46,[66][67][68]. We wrote Eq. (3.11) using the κ 1,2 / k p factors to conform with the known results for the transverse spin in an evanescent wave in free space [35,[41][42][43]: (3.11) shows that the Abraham-type spin density S A is discontinuous at the interface x = 0 . As a result of this, the canonical and spin parts of the Abraham-Poynting energy flux have delta-function singularities [65], originating from the gradient terms in Eq. (2.5). Note also that the integral value S A y is always positive and is in agreement with calculations of Ref. [65] up to a factor of 2 missing there. [The missing factor of 2 in [65] originates from the improper application of the relation dx to the z -delocalized SPP wave, involving only the term ∂S Ay / ∂x under the integral. In fact, this relation is valid only for localized wave packets involving two terms: ∂S Ay / ∂x and ∂S Ay / ∂z .] At the same time, calculations of the integral Abraham spin "per particle" for surface Maxwell modes in [89] are not applicable in the SPP case because a dispersion-free model without proper Brillouin energy (2.3) was considered there. The frequency dependence of the Abraham-type spin (3.12) (in units of ! "per plasmon") is shown in Fig. 4(a). We now calculate the properly defined canonical Minkowski-type spin and orbital AM (2.13). Substituting Eqs. (3.2) and (3.3) into Eq. (2.13), we find The spin AM density (3.13) is directed oppositely in the vacuum and metal: S y < 0 for x < 0 .
This agrees with the opposite direction of rotation of the electric field E in the metal, Eq.
(3.14) Equation (3.14) provides the first accurate calculation of the transverse spin AM carried by a SPP. In contrast to the Abraham-like spin (3.12) considered before, its direction can vary depending on the SPP frequency. Namely, it is positive, ! S My > 0 , for ω < ω p / 3 and negative, ! S My < 0 , for ω p / 3 < ω < ω p / 2 , so that it vanishes at ω = ω p / 3 . The frequency dependence of the Minkowski-type spin (3.14) (in units of ! per plasmon) is shown in Fig. 4(a). Note that its absolute value never exceeds ! / 2 per plasmon; this is because of the pure-electric origin of the transverse spin, with no magnetic part. Interestingly, the critical zerospin value ω = ω p / 3 , or ε = −2 , corresponds to the elliptical x, z ( ) -polarizations of the electric field (3.2) with the axes ratios 2 and 1 / 2 (i.e., identical ellipticities but different orientations), and opposite directions of the rotation, in the vacuum and metal, respectively. The orbital AM density in the SPP field is determined by Eq. (2.13) and the canonical momentum density (3.9): ! L M = r × ! P M . This quantity is extrinsic, i.e., depends on the choice of the coordinate origin. However, the expectation value of the orbital AM can have both extrinsic and intrinsic contributions [41]. The intrinsic orbital AM is calculated with respect to the centroid of the energy density distribution. Using the x -shifted centroid (3.6) and z -directed momentum density (3.9), we calculate the y -directed intrinsic orbital AM of the SPP: Thus, the intrinsic orbital AM of SPP vanishes. This is because of the proportionality between the canonical momentum density (3.9) and energy density (3.4), which in turn determines the centroid (3.6). The vanishing intrinsic orbital AM (3.15) reflects the non-vortex character of the canonical momentum (3.9), in contrast to the circulating Abraham-type energy fluxes [24,65,85]. The extrinsic part of the orbital AM, calculated with respect to the interface x = 0 , can be written as

Microscopic fields and parameters of electron plasma.
The microscopic approach is based on the separation of the microscopic electromagnetic field ( E , H ) and charges/currents inside the medium. In our case, the metal can be described using the Bloch hydrodynamic model for electron plasma. In this model, the electron density is written as n r,t ( ) = n 0 + ! there is no effective medium, ε = µ = 1 ), but the Maxwell equations are modified by the presence of charges and currents [22]: These equations describe the influence of the electrons on the fields. The back action is described by the hydrodynamic equation for the electron gas [71]: This equation is a classical equation of motion of the electron in the electric field E (the Lorentz force from the magnetic field vanishes in the linear problem) with the additional quantum pressure term involving the coefficient where v F is the Fermi velocity of electrons. Our classical treatment of the SPP wave implies the limit β 2 → 0 . However, we cannot omit the last term in Eq. (4.2) from the beginning because it is crucial to satisfy the boundary conditions at the metal-vacuum discontinuity, x = 0 . Solving Eqs. (4.1) and (4.2) with standard boundary conditions at the metal-vacuum interface (continuity of E x , E z , H y , and vanishing of v x ) yields the microscopic electric and magnetic fields as well as electron plasma properties in the SPP wave. The magnetic field (both in the vacuum and in the metal) and electric field in the vacuum are still described by the macroscopic Eqs. (3.2) with parameters (3.1) and (3.3), while the electric field inside the metal becomes: Here γ 2 = k p 2 − ε ω 2 / β 2 , we still use ε as a parameter given by Eq. (3.1), and β 2 → 0 implies γ → ∞ . The electron density and velocity perturbations in SPPs are given by We can now consider the classical limit γ → ∞ . Since the γ -terms appear only in the metal half-space x < 0 , one can use the limiting transition γ exp γ x ( ) → δ x ( ), where the deltafunction describes surface effects at the interface. After doing so, the electric field (4.3) becomes equal to the macroscopic one, Eq. (3.2), while the electron density vanishes in the volume, ! n = 0 for x < 0 , and acquires a surface delta-function singularity: As we show below, this singularity is cancelled by another singularity in field gradients, and all dynamical properties of the SPP wave are determined by volume contributions in the metal and in the vacuum. Finally, the electron velocity (4.5) becomes proportional to the electric field (3.2) [which follows from Eq. (4.2) at β = 0 ]: Note that it is the vanishing of the electron density perturbation ! n in volume that makes the microscopic fields "transverse", i.e., divergence-free: ∇ ⋅ E = 0 . Because of this, we do not need to consider contributions of "longitudinal" (i.e., curl-free) fields to the energy, momentum, and AM [90].
In addition to the point-charge features of electrons, we will need their electric-dipole properties. Since velocity is a time derivative of the position of the electron, we can write the complex amplitude of the electron displacement as ! a = i ω ! v . From here, the complex amplitude of the density of the electron dipole moment is Substituting here Eq. (4.7), we can write the dipole-moment density (4.8) as Here α is the dipole polarizability of the electron gas, and the last equality shows that it is in perfect agreement with the macroscopic theory based on the permittivity ε [22]. Indeed, substituting the velocity (4.7) into Maxwell equations (4.1) for microscopic fields, we immediately obtain the source-free Maxwell equations (2.2) for macroscopic fields with permittivity ε .

Microscopic calculations of energy and momentum densities.
In the vacuum half-space x > 0 , the microscopic and macroscopic electromagnetic fields and their properties coincide, so we have to compare only the macroscopic and microscopic properties in the metal. Hereafter, we consider all quantities only in the x < 0 half-space.
The cycle-averaged energy density in the system of microscopic electromagnetic fields and electrons can be written as Here the first term is the microscopic-field energy [written as for free-space fields, Eq. (2.1)], and the second term is the kinetic energy of electrons. Note that the latter can also be presented as the energy of the dipole (4.8) and (4.9) in the electric field: Substituting fields (3.2) and velocity (4.7) into Eq. (4.10) results in the macroscopic energy density ! W , Eq. (3.4), at x < 0 . Thus, the microscopic and macroscopic calculations are in perfect agreement. The field and electron contributions in the metal are: The momentum density in the microscopic approach is also the sum of the field and electron contributions. The velocity of electrons in external fields is associated with the kinetic rather than canonical momentum of electrons, and, therefore, we consider the corresponding kinetic momentum density of the field. [It is worth noticing that the canonical electron momentum in the metal vanishes:  [16,17]: (4.12) The calculation of the electron contribution P mat requires a more sophisticated approach. Indeed, the simple expression n 0 m! v = i en 0 ω E (with oscillating, zero-average ! v ) for point electrons does not provide a meaningful result; instead of this, one has to consider an optical force acting on electric dipoles (4.8) and (4.9). We follow the formalism described in the review [4], Section 5.1 therein. Namely, we consider a long but finite wave packet instead of a monochromatic continuous wave. Afterwards, the length of the wave packet can be tend to infinity. Introducing slowlyvarying amplitudes E r ( ) → E r,t ( ) , H r ( ) → H r,t ( ) , etc., with the typical scale of the tvariations much larger than ω −1 , involves the corresponding narrow but finite frequency Fourier spectrum centered at ω . This produces the first-order Taylor-series correction to the relation (4.9) between the dipole moment and electric field: Next, we consider the cycle-averaged force density acting on the dipole moment (4.13) in an external electromagnetic field [4]: (4.14) Substituting here Eq. (4.13), after some transformations the force can be written as: where W mat = −α E 2 / 4 . The first term in Eq. (4.15) represents the gradient force, while the two terms subject to the time derivative should be associated with the momentum density carried by the electrons, i.e., P mat of Eq. (4.12). Expressing the dipole-moment density and polarizability via the electric field and ε , Eq. (4.9), we arrive at: (4.16) Here the second term is associated with the "Abraham force" [22]. Substituting the electron momentum density (4.16) into Eq. (4.12), results in the kinetic form of the Minkowski-type momentum density (2.8) (for µ = 1 ): Thus, using microscopic calculations for the SPP, we rigorously derived the kinetic momentum density (2.10) suggested previously from a phenomenological formalism [16,17]. We now trace the decomposition of the kinetic momentum (4.17) into the canonical (orbital) and spin parts (2.9) and (2.10). In principle, one substitutes the electron velocity ! v = ie / mω ( )E into the microscopic Maxwell equations (4.1), which results in the macroscopic Maxwell equations (2.2) with permittivity ε . Then, the decomposition becomes straightforward as described in Section 2. However, it is instructive to trace this decomposition at the microscopic level. For this purpose, we decompose the Poynting-like part of the kinetic momentum (4.17) into canonical and orbital parts using the microscopic Maxwell equations (4.1). The presence of sources, ! n and ! v , modifies this decomposition as compared to the freespace case (1.6)-(1.8): (4.18) Here we ascribed the two source-related terms to the canonical and spin parts such that the final result will coincide with the macroscopic equations. First, for the SPP fields (3.2) and velocity (4.7) the velocity-related term contains only the canonical-type contribution [because Im H * × H ( ) = 0 in SPPs]: Combining the canonical part of Eq. (4.18) with Eq. (4.19) and the second term in Eq. (4.17), we obtain the macroscopic canonical Minkowski-type momentum (2.9) (with µ = 1 ): Second, the curl of the free-space-like spin S 0 for the SPP fields (3.2) contains a deltafunction singularity at the metal-vacuum interface x = 0 . Remarkably, this singularity is exactly cancelled by the singularity (4.6) in the electron density distribution, so that the spin part in Eq. (4.18) becomes non-singular. This confirms the non-singular character of the spin-orbital decomposition in the Minkowski-type approach, in contrast to the Abraham one [65]. In the bulk, x < 0 , the spin part of Eq. (4.18) substituted in Eq. (4.17) immediately yields the corresponding part of the macroscopic Minkowski-type momentum involving the "naïve" Minkowski spin, Eq. (2.10): (4.21) Thus, we have obtained the macroscopic Minkowski-type momentum densities (2.8)-(2.10), both kinetic and canonical, using microscopic calculations in the metal with separated field and matter contributions.

Microscopic approach to the spin and orbital AM.
It might seem that the description of the AM quantities in dispersive media, given in Section 2, is somewhat "inconsistent" with the corresponding momentum quantities. Indeed, the Minkowski-type kinetic AM (2.12) is not simply determined by the corresponding kinetic momentum, r × ! P M , but contains additional dispersion-related terms. Furthermore, the spin momentum (2.10) is determined by the naïve Minkowski spin density S M , while the proper canonical spin AM ! S M , Eq. (2.13), differs from it in a dispersive medium. Importantly, the microscopic approach sheds light on these "inconsistencies".
Namely, the local motion of electrons provides an intrinsic contribution to the AM density [71], in fact, to its spin part. Using the electron displacement ! a , Eq. (4.8), and velocity ! v , Eq. (4.7), one can write this part of the AM density as: This term exactly describes the dispersion-related addition in the Minowski-type kinetic AM (2.12) for the SPP wave: Thus, akin to the momentum density (2.8) and (4.175), microscopic calculations justify the kinetic AM density (2.12), previously obtained by Philbin within a phenomenological approach [17]. Consider now the spin and orbital AM densities. The orbital AM density is straightforwardly determined by the canonical momentum (4.20): ! L M = r × ! P M , and we have already described its properties in Section 3. At the same time, the intrinsic electron contribution (4.22) elucidate the difference between the naïve Minkowski spin density (4.21) and canonical spin density (2.13): This justifies the use of the canonical spin AM in dispersive media and the transverse spin of a SPP calculated in Eqs. (3.13) and (3.14), Fig. 4(a). The dispersion-related contribution is absolutely crucial in the case of SPPs, because ε < 0 , ! ε = 2 − ε > 0 , and it changes both the magnitude and sign of the spin AM density in the metal.

Magnetization, magnetization current, and Abraham momentum.
The electron contribution (4.22) to the spin AM of a SPP corresponds to the microscopic circular motion of electrons in the SPP field. This microscopic orbital motion of the electrons produces multiple circulating currents, and, hence, the constant (non-oscillating) magnetization of the metal. Using the standard gyromagnetic ratio, we obtain the magnetization density in the metal: This equation exactly coincides with the results [91,92] obtained for the magnetization of plasma by electromagnetic radiation and the inverse Faraday effect [22]. For the SPP fields (3.2) and (3.3), we find Thus, the metal is magnetized along the positivey direction ( e < 0 ). The magnetization (4.26) means that the SPP, being a mixed photon-electron excitation, carries a non-zero magnetic moment. To characterize this magnetic moment "per plasmon", we calculate the integral magnetization (4.26): This corresponds to the magnetic moment µ = 2 −ε 1+ ε 2 µ B y per plasmon, where µ B = e ! / 2mc is the Bohr magneton. The absolute value of this magnetic moment grows from 0 to µ B as the SPP frequency ω changes from 0 to ω p / 2 . Moreover, the inhomogeneous magnetization (4.26) generates the corresponding magnetization electric current j magn = c∇ × M : This is a direct current which flows in the metal in the z direction, i.e., along the SPP propagation. It should be emphasized that the current (4.28) is obtained as a quadratic form of the SPP fields. Indeed, the linear-approximation current is determined by the electron velocity ! v and vanishes after cycle averaging.
We also note that the magnetization current is solenoidal (divergenceless). Therefore, it does not contribute to the charge transport and cannot be measured by an ammeter or voltmeter. Nonetheless, one can determine the electron velocity v magn = j magn / en 0 and momentum density P magn = mn 0 v magn = m / e ( ) j magn corresponding to the direct current (4.28). Using Eqs. (3.3), we write it as Thus, the magnetization momentum (4.29) is directed oppositely to the SPP propagation. Most significantly, it is exactly equal to the difference between the kinetic Abraham and Minkowskitype momentum densities in the metal, Eqs. (3.7) and (3.10): Equation (4.30) completes the microscopic picture and explains the origin of the difference between the Abraham and Minkowski momenta in the medium. Since this difference is produced by the direct magnetization current, it cannot be attributed to the wave (Minkowski-type) momentum but it does contribute to the energy flux (Abraham-Poynting momentum) and the group velocity of SPPs. Note that the Abraham momentum density can also be obtained as a pure microscopic-field momentum P A = P 0 , see Eq. (4.12) [19]. It follows from here that the total momentum of the metal vanishes in the problem under consideration: the electron contribution to the wave momentum is exactly cancelled by the magnetization-current momentum: P mat + P magn = 0 .
Interestingly, recent work [93] claimed that the magnetization of the medium should be twice as small as compared to Eq. (4.22) due to the contribution from a "drift current". If this were to be the case, then adding the corresponding magnetic current and momentum (4.29) to the Minkowski-type momentum would yield an average between the Minkowski and Abraham momenta in Eq. (4.30), similar to the suggestion by Mansuripur [6,94].

Helicity density and flux in a medium.
Here we briefly consider problems related to the dual symmetry between electric and magnetic fields in Maxwell equations [52,[54][55][56][57]69,70,75,76,79]. This symmetry is exact in freespace Maxwell equations, but the presence of electric charges and currents in matter breaks it. The dual symmetry corresponds to the conservation of electromagnetic helicity via Noether's theorem [55,57,69,70,75,76,[95][96][97][98]. Using the same phenomenological Lagrangian formalism as for the derivation of Minkowski-type momentum and AM densities (2.8) and (2.12), Philbin obtained the helicity density in a dispersive medium [77]: In a dispersion-free medium, ! ε = ε and ! µ = µ , and Eq. (5.1) yields the free-space result [55,57,59,69]. It should be emphasized that ! K is the helicity density, but not the chirality density (Lipkin's zilch); these quantities are simply proportional to each other only in free space [57,77,[99][100][101][102][103]. Similarly to Eqs. (2.17), for a plane wave in a homogeneous transparent medium, Eq. (5.1) yields This means that the helicity "per photon" in units of ! is limited in the medium by the −n p −1 ,n p −1 ( ) range. This is a strange result without clear physical meaning. Using the quantum wavefunction formalism (2.15) and (2.16) with the helicity operator [55,59] but only in the case of a non-dispersive medium. Therefore, the helicity density in a dispersive medium is still controversial and requires further investigation. In any case, for the SPP fields (3.2), the helicity vanishes [65] due to the orthogonality of the electric and magnetic fields. Note also that in free space the dual-symmetric spin AM density (1.7) determines the helicity flux density [55,57,69,70,100]. In a dispersive medium, calculations in Ref. [77] showed that the flux density of the helicity (5.1) is determined by the Abraham spin density S A , Eq. (2.6). In this case, the helicity flux density does not coincide with the proper Minkowskitype spin AM density (2.13) in a medium, and these are two different physical quantities. Even considering the Lipkin's chirality density, its flux becomes proportional to the "naïve" Minkowski spin density S M [77,102], which is different from the canonical spin density (2.13) in dispersive media.

Dual-symmetric and asymmetric quantities.
So far, we considered all definitions of the optical energy, momentum, and AM using forms symmetric with respect to the electric and magnetic fields (and, correspondingly, indices ε and µ ) [35,41,42,46,52,[54][55][56][57]59,65,69,70]. However, in standard electromagnetic field theory (or QED) the field Lagrangian is not dual-symmetric [36,44,55,69,70]. Due to this, the canonical momentum, spin, and orbital AM densities are often defined using dual-asymmetric field-theory expressions, which contain only the electric-field parts [36,44,[47][48][49][50]60,90]. In such "standard" formalism in free space, the energy density and kinetic Poynting momentum density are still given by the dual-symmetric expressions (1.1) and (2.1), whereas the canonical momentum, spin, and spin AM densities (1.6)-(1.8) become [54,55]: Here, we introduced the electric and magnetic parts of the momentum and spin densities (1.6) and (1.7): P 0 = P 0 e + P 0 m and S 0 = S 0 e + S 0 m . Adopting the definitions (5.4) and (5.5) would mean that the only the phase gradients of the electric (but not magnetic) field produce momentum and orbital AM, and only rotations of the electric (but not magnetic) field generate the spin AM. On the one hand, this is not satisfactory from a general physical perspective, where electromagnetic waves involve electric and magnetic fields on equal footing [52,55,57,56,[81][82][83]. On the other hand, the electric-biased quantities (5.4) and (5.5) can be useful in some practical problems considering interactions of fields with electric-dipole particles or atoms, which are not sensitive to magnetic fields [47][48][49]59,60]. However, for magnetic-dipole particles, the electric-biased quantities (5.4) and (5.5) do not make sense. Note that in free space, the dual-asymmetric densities (5.4) and (5.5) do not cause significant problems because: (i) the Poynting-vector decomposition still has the same form (1.8): P 0 = P 0 ′ + 1 2 ∇ × S 0 ′ , and (ii) the integral values of all these quantities for localized freespace fields coincide with the dual-symmetric definitions [52,55]: However, in a medium, which is dual-asymmetric ( ε ≠ µ ) in the generic case, the difference between the dual-symmetric and asymmetric definitions becomes crucial. We again consider the example of a SPP wave at a metal-vacuum interface. First, the electric and magnetic contributions to the energy density (2.3) and (3.4) for the macroscopic SPP fields (3.2) equal Obviously, these contributions are different and result in different integral electric and magnetic energy parts: Then, the electric and magnetic parts of the canonical Minkowski-type momentum (2.9) and (3.9) become proportional to the corresponding energy parts: From here and the difference in the integral electric and magnetic energies (5.8), it follows that using an electric-biased momentum density similar to Eq. (5.4), Thus, the dual-asymmetric definitions in a medium do not satisfy the convenient free-space relations (5.6). In the case of Eq. (5.10), this breaks the natural proportionality (3.9) and the momentum of the SPP becomes not equal to !k p per plasmon. Obviously, this is a physically unsatisfactory result. Therefore, we conclude that dual-symmetric definitions of momentum and AM densities are crucial for structured waves in inhomogeneous optical media. Moreover, the validity of the dual-symmetric (rather than electric-biased) formalism follows from the microscopic calculations of Section 4. Indeed, the dispersion-related material The transverse spin AM of SPPs is a quantity which is extremely sensitive to the duality. Indeed, for the metal-vacuum interface considered here, the transverse spin has a purely electric nature: only the electric field rotates in solutions (3.2), while the transverse magnetic field yields no contribution to the spin AM [65]. Using an electric-biased definition of the spin AM, similar to Eq. Overall, we conclude that only the dual-symmetric definitions of the momentum, spin, and orbital AM properties of optical fields in dispersive media provide physically consistent and satisfactory results.

Concluding remarks
We have examined momentum and angular-momentum (AM) properties of monochromatic optical fields in dispersive and inhomogeneous media. The two major problems that lie at the heart of this study are: (i) the Abraham-Minkowski dilemma and (ii) the canonical spin-orbital decomposition of optical momentum and AM. We have shown that, in principle, one can formulate four momentum-AM pictures using Abraham-type and Minkowski-type quantities in kinetic (i.e., Poynting-like, without spin-orbital separation) and canonical (i.e., with spinorbital separation) approaches. These four pictures are summarized in Table I. However, two of these sets of quantities are more physically meaningful.
First, the Abraham-Poynting kinetic momentum density (1.1) and (1.2), P A = P 0 , should be associated with the energy flux density rather than the momentum density. This quantity determines the energy transport and group velocity of a wave packet in the medium.
Second, the canonical Minkowski-type momentum density (2.9), ! P M , together with the corresponding spin and orbital AM densities (2.13), ! S M and ! L M = r × ! P M , provide a physicallymeaningful and self-consistent description of the momentum and AM of light in the medium. To the best of our knowledge, these quantities were derived for the first time in the present work, but these are consistent with several previously used approaches. On the one hand, the kinetic counterpart of this canonical picture, momentum (2.8), ! P M , and total AM (2.12), ! J M , exactly coincide with those obtained in the most general form by Philbin and Allanson [16,17]. On the other hand, our canonical characteristics (2.9) and (2.13) have a much more elegant form, exhibiting a pleasing similarity with the Brillouin energy density (2.3), ! W , in the medium. As a result, the energy, momentum, spin, and orbital AM densities in the medium can be written in a laconic unified form (2.15) using the corresponding quantum-mechanical operators and proper inner product modified by the ! ε and ! µ indices of the medium. This coincides with the general approach to electromagnetic bi-linear forms developed by Silveirinha [81][82][83]. We applied the above general theory to a surface plasmon-polariton (SPP) wave at a metalvacuum interface. This example provides a deep physical insight because it involves essentially inhomogeneous fields as well as an inhomogeneous, non-transparent, and dispersive medium. This is in sharp contrast to plane waves in homogeneous transparent media, which are considered in the majority of the Abraham-Minkowski studies. We have shown that in the nontrivial SPP field, the integral Abraham-Poynting momentum P A describes the group velocity v g = ∂ω / ∂k p < c of SPPs, Eq. (3.8), while the integral canonical momentum ! P M corresponds to the super-momentum !k p > !k 0 per plasmon, Eq. (3.9). This is the first example of a wave, which carries an integral momentum larger than that of a photon in vacuum, and this originates from the inhomogeneous-evanescent character of the surface wave rather than from the medium refractive index ( k p → ∞ at ε → −1 in the metal).
We have also provided the first accurate calculation of the transverse spin ! S M of a SPP, Eq. (3.14). The result differs considerably from previous calculations using the Abraham-type definition of the spin [65]. In particular, the integral transverse spin AM of a SPP can vanish (at ω = ω p / 3 ), change its sign and reach the value −! / 2 per plasmon, Fig. 4(a). In turn the intrinsic orbital AM (calculated with respect to the center of energy) of a SPP vanishes, Eq. (3.15). This agrees with the non-vortex character of the canonical momentum density (phase gradient) in the SPP field, and is in contrast to the Abraham-Poynting circulating energy flux [65]. Thus, the SPP example shows that the Abraham-type kinetic and Minkowski-type canonical properties provide intuitively clear and consistent description of complex optical fields in complex media. Importantly, we have also provided microscopic calculations of the momentum and AM densities for the SPP field, considering the microscopic electromagnetic fields and motion of free electrons in the metal. These calculations resulted precisely in the Minkowski-type momentum and AM densities, previously suggested from macroscopic phenomenological approaches [16,17]. This proves the validity of our approach and illuminates its physical origin. Importantly, the dispersion ε ω ( ) in the metal was absolutely crucial in the above calculations, affecting not only the magnitudes but also the signs of the dynamical quantities (because of ε < 0 and ! ε > 0 ). Using the microscopic theory, we have also predicted a transverse magnetization of the metal (the inverse Faraday effect) corresponding to the transverse spin of a SPP. This means that a SPP wave carries not only the spin AM but also the transverse magnetic moment, up to a Bohr magneton per plasmon. Furthermore, an inhomogeneous magnetization produces the direct magnetization current flowing along the metal surface in the SPP propagation direction. Remarkably, the momentum density corresponding to this current is exactly equal to the difference between the Abraham and Minkowski-type momenta in the metal, thereby providing one more "resolution" of this longstanding problem [1][2][3][4][5][6][7].
Finally, we briefly discussed the optical helicity and duality problems in dispersive media. The analysis of the SPP example leads us to conclude that the dual-symmetric description of the canonical momentum, spin and orbital AM, with symmetric electric-and magnetic-field contributions, is crucially important for the physically meaningful and consistent picture of these properties in dispersive media. In particular, we found that microscopic calculations, including dispersion-related material terms, are consistent with the dual-symmetric (rather than electricbiased) formalism.
Thus, the present study provides a complete analysis and description of the momentum and AM of light in dispersive and inhomogeneos (but isotropic and lossless) media. We have considered SPPs only as the simplest example of the application of our theory, where other approaches fail. Taking into account both material and structured-light properties is crucial in a variety of nanooptical and photonic systems, including photonic crystals, metamaterials, and optomechanical systems. Our theory provides an efficient toolbox for the description of dynamical properties of light in such systems. One of the main tasks for future studies is to extend this analysis to other classes of materials, including dissipation or gain and anisotropy. In particular, it is not clear if one can separate the spin and orbital degrees of freedom in anisotropic media. Close correspondence of our approach to some of the results of Refs. [18,[81][82][83]102] (dealing with quite general bi-anisotropic media), suggests that the analysis presented in this work can be extended to more complex cases.
Another important direction for future consideration is whether the fundamental wave characteristics introduced in this work can be observed in experiments. There were experiments measuring Minkowski-type momentum ( ∝ !k ) for plane waves in dispersive media [13,104,105]. In addition, the canonical momentum and spin AM densities in structured optical fields in free space are directly observable via the optical force and torque on dipole particles or atoms [31][32][33][34][35]41,46,47,[59][60][61][62][63]66]. Therefore, now it would be important to calculate and measure optical forces and torques on small particles in dispersive media (e.g., liquids or gases), and check if these are proportional to the canonical momentum and spin densities in dispersive media. Furthermore, it would be important to find an experimental setup where the integral super-momentum and transverse spin of SPPs, derived in this work, could be detected. So far, only local densities of these quantities were accessible via local free-space measurements using small particles or atoms [31][32][33][34][35]41,66,67,[106][107][108][109].
Note added. A short (letter) version of this work, highlighting its main results without derivations, has recently been published in [110]. After the submission of this work, we became aware of the recent preprint [111], where material corrections to the electromagnetic spin density, similar to those obtained in our work, were found in the form of a torque in the spin continuity equation.