A proposal to measure Belinfante’s curl of the spin optical force based on the Kerker conditions

The linear momentum of electromagnetic radiation is proportional to the Poynting vector. However, when light interacts with a dipole, the radiation pressure generated is no longer proportional to the Poynting vector: the so-called Belinfante’s momentum or equivalently, the curl of the spin density of the light field, must be considered. In this paper, we propose an optical configuration, based on two evanescent counter-propagating waves, perpendicularly polarized, capable of detecting Belinfante’s mechanical action. The two beams interact with a high-refractive-index particle like silicon. The direction of the radiation pressure exerted on the particle, proportional only to the curl of the spin density, depends on the electric and magnetic response of the particle and changes sign at the radiation wavelengths corresponding with the Kerker conditions.


Introduction
Radiation pressure is a direct manifestation of the existence of light linear momentum which is transferred to objects inducing scattering forces. This linear momentum density is calculated dividing the electromagnetic power per unit area by the speed of light squared. By considering the energy conservation in Maxwell equations, the power per unit area is given by the Poynting vector S = E × B/μ 0 deriving in a linear momentum given by < P >=< S > /c 2 [1]. There are some debates related to the electromagnetic linear momentum like, for instance, the hidden momentum, the electromagnetic mass and the Abraham-Minkowski controversy [2] . In this paper, we are going to deal with the existence of the so-called virtual Belinfante's linear momentum.
The linear momentum of light may be obtained also using classical field theories by the principle of least action. However, following this path, a different result is obtained [3]. The momentum obtained, the so-called canonical momentum < P o >, is given by < P > plus a new additional solenoidal contribution that we may denote as < −P s >. In principle, since this new term is divergence free, it has no influence in Maxwell's conservation laws. However, the existence of this new term implies a non-symmetric stress energy tensor and then, in order to establish a symmetric field theory, < −P s > is usually discarded [4] and light's linear momentum is considered as < P >=< P o > − < −P s >. That is the reason a e-mail: manuel.marques@uam.es (corresponding author) why this solenoidal term < −P s >, known as Belinfante's momentum, is usually considered to be "virtual." However, this term has a fundamental importance when dealing with the radiation pressure induced by light on electric and magnetic dipoles. Light forces on dipoles are described as the sum of the gradient force and the radiation pressure or scattering force. The gradient force is proportional to the real part of the polarizability and the intensity gradient of the field, while the scattering force is proportional to the imaginary part of the polarizability and the phase gradient of the radiation [5]. Several theoretical studies have confirmed that the scattering force on a small particle is not given by the Poynting vector [6,7], there is an extra contribution coming from the curl of the spin density of the light field [8] implying that the momentum transferred is really the orbital part of the Poynting vector [9][10][11][12][13][14]. Therefore, we may say that the momentum transferred to the particle from electromagnetic radiation is not given by < P > but by < P o >=< P > + < −P s >.
This mechanism for momentum transfer cannot be experimentally confirmed using homogeneously distributed radiation, where < P s >= 0 and < P >=< P o >. However, if we consider an elliptically polarized, inhomogeneous field, then the curl of the spin density of the light field is different from zero [8] and the scattering force is no longer given by the full Poynting vector. By measuring the mechanical action of this kind of radiation, with a Belinfante's momentum different from zero, we should be able to experimentally detect how the radiation pressure on a small particle is no longer proportional to the Poynting vector.
Several experimental setups have been proposed in order to detect the special mechanical action on a dipole due to the existence of a curl of the spin density different from zero. The main problem is that usually, the Poynting vector and Belinfante's momentum are parallel and then, it is very difficult to separate both effects. We may say that there are sophisticated optical configurations, in which the action of the Poynting vector and the Belinfante's momentum are non-parallel even within the dipolar regime [15]. The other option is to look for optical configurations with a null value of the Poynting vector. In these cases, the scattering force on a dipolar particle is only due to Belinfante's momentum. These kinds of approaches are based on optical lattices [16,17] and evanescent waves [18].
Recently, a different type of forces coming from the electric and magnetic dipole interaction [19] have been studied. It is possible to find optical configurations where these forces are proportional to the transverse part of the curl of the spin density of the electromagnetic radiation [20][21][22][23][24][25][26][27][28][29]. However, these transverse forces, affecting Mie particles but not Rayleigh dipoles, are not the forces that we are treating on our paper. The forces we are about to analyze come fully and solely from the transfer of the Belinfante's momentum to a unique dipole (they show up also in configurations with a null value of the Poynting vector), are obtained from the interference between the incident and scattered fields and act separately on electric or magnetic dipoles [19]. In opposition, the transverse forces previously reported (which are zero if E × B * = 0) are obtained from the scattered fields, come from the interaction between the electric and magnetic dipoles and, only for particular cases, happen to be proportional to the transverse part of the Belinfante's force. In addition, we may say that, spin dependent forces are not univocally given by the curl of the spin; it is also possible to find systems presenting interactions that depend on the spin of the electromagnetic radiation, like the chiral systems [30][31][32] or proportional to the gradient of the spin, like the magneto-optical particles [33].
In this paper, we propose and analyze a setup based on evanescent waves and on the magnetic response of dielectric particles with a large value of the refractive index . On the one hand, transverse magnetic (TM) evanescent waves have a curl of the electric spin of the light field different from zero, affecting electric dipoles, while transverse electric (TE) evanescent waves have a curl of the magnetic spin of the light field different from zero, affecting magnetic dipoles. On the other hand, high-refractive-index scatterers like silicon micro-particles behave like a magnetic or electric dipole depending on wavelength [34]. Then, if a silicon particle is illuminated with a combination of evanescent waves, the mechanical response is going to depend on wavelength and on the polarization (TE or TM) of the evanescent wave. In this paper, we use these properties to propose an optical configuration where the characteristic wavelength dependence of the radiation pressure can be only attributed to the existence of a curl of the spin density of the light field different from zero or to the "virtual" mechanical action emanating from Belinfante's momentum.

The electromagnetic fields
We are going to consider two electromagnetic fields consisting of two evanescent waves in z > 0 (TE (or s) and TM (or p) polarized) propagating in x and with the same amplitude. The TE field is given by: (2) and the TM field by: with being a the attenuation vector, δ the phase difference between both waves, k the wave vector, n 2 = √ the refractive index of the medium where the evanescent wave propagates (z > 0) and E 0 the amplitude of the y component of the electric field.
Note that 2.2 The scattering forces Let us consider a particle that behaves like an electric and a magnetic dipole depending on wavelength. The general expression for the time-averaged scattering force on a dipole with electric and magnetic response is given at first-order approach by [19] Eur. Phys. J. Plus (2021) 136:185 being < S >= c 2 0 (E × B * )/2 the time-averaged Poynting vector. < L e s >= 0 (E * × E)/i2ω and < L m s >= (B * × B)/i2ωμ 0 are the time-averaged electric and magnetic spin densities of the light fields and α ' e , α ' m are the imaginary parts of the electric and magnetic polarizabilities.
For the particular field considered, since electric field and complex conjugate of the magnetic field are parallel (see Eq. 7), < S >= 0. On the other hand, the electric and magnetic spin densities are simply related by < L m s >= − < L e s >: Thus, the force is given by which can be explicitly written as being i the unit vector in the x direction. Note that this force is just proportional to the time-averaged value of Belinfante's spin linear momentum density: Since B and E * are parallel, all forces coming from higher-order electric and magnetic dipoles interaction (α e α * m ) are equal to zero [19]. In the z direction, we have also forces coming from the gradient of the evanescent field intensity that may be compensated with the optical trap holding our particle over the surface.
To sum up, if the momentum transferred to the particle comes from < P o >=< P > − < P s > and not from < P >, with this configuration, we should be able to detect that, although < S >= 0, the force in the x direction is not zero and depends on Belinfante's momentum and on the value of (α ' m − α ' e ).

The particle
As a probe, we propose to use a particle with tuneable electric and magnetic dipolar response depending on wavelength. For example, we can use a silicon particle, with a diameter of half a micron and radiation wavelengths ranging from 1.4 to 2 microns. For plane wave incidence, the total extinction cross section calculated from Mie theory, the extinction coming from the electric and magnetic dipoles and the sum of the contributions coming from the electric and magnetic dipoles are plotted in Fig. 1. Note how the particle can be considered just like an electric and a magnetic dipole and how, for wavelengths smaller than roughly (1700nm), the electric response dominates while, for larger wavelengths, the magnetic response is the one prevailing . At 1700nm responses (electric and magnetic) are the same and the scattering force should be zero (see equation 11).
A close analysis to the values of the polarizability reveals that, close to 1700nm, α m = α e and also, α m = −α e . Then, the wavelength at which a Belinfante's force sign change is expected, corresponds with the generalized second Kerker condition [35]. In addition, the other wavelength (close to 2000nm) with zero value of the Belinfante's force in this

Description of the experiment proposed
We consider a semi-spherical glass prism (n 1 = 1.55). The flat surface is immersed in water (n 2 = 1.33) and the evanescent wave in water is created by the total internal reflection in the glass-water interface. We consider the emission from a laser with a 5 mm spot diameter that is split into two beams with the same intensity (100 mW) using a polarized beam splitter (see Fig. 2). Each beam is focused using a lens with 70 mm focal distance. Both beams enter perpendicularly on the glass surface and, due to the reflection on the glass, the power transmitted to the prism decreases in a 0.95 factor per beam.
Since the fields (TE and TM) are induced by the total internal reflection of two plane wave with the same intensity (s-polarized and p-polarized, respectively), propagating in a medium with n 1 > n 2 , the following condition on the transmission coefficients must be fulfilled |t H ⊥ |/|t ⊥ | = n 2 /n 1 . This condition will hold as long as the incident angle (θ ) for both waves takes the value In our case, this angle corresponds roughly to θ = 67 o being 59 o the critical angle needed to obtain total internal reflection.
The attenuation and k vectors are obtained from the incident angle and the wavelength in vacuum (λ 0 ): Eur. Phys. J. Plus (2021) 136:185 Fig. 2 Proposed experimental setup. The beam coming from an infrared non-polarized laser is split into two by a polarized beam splitter (dark blue). Both beams, s-polarized (crosses) and p-polarized (arrows) are reflected on a mirror (light blue) and focused on the same position of a glass-water interface. Two evanescent waves, TE and TM, with the same intensity and counter-propagating are generated on the water surface, inducing scattering forces on a silicon particle (yellow circle) proportional to Belinfante's spin momentum and depending on the value of the electric and magnetic dipoles induced on the particle In order to obtain this final incident angle (θ ), we use two tilted mirrors where the beams are reflected with the proper angle [37].

The radiation pressure
We have calculated, using equation 10, the force generated in the x direction on the half a micron diameter silicon particle as a function of wavelength. In order to ensure strong enough forces, we have considered that the spherical particle is trapped at a distance z = 100 nm from the glass surface. The expected measured force as a function of wavelength is shown in Fig. 3. For this calculation of the forces, we have considered that the power is the same for all wavelengths and we neglect possible corrections coming from the beam splitter performance. We consider the electric and magnetic polarizabilities inferred from Fig. 1. Corrections to the longitudinal scattering force felt by the electric dipole coming from the multiple scattering with the dielectric surface are only of the order of 3/1000 [38].
Note how the scattering force, coming solely from the curl of the spin density of the light field, i.e., from the spin or Belinfante's linear momentum, should go from positive to negative, being zero at the wavelengths at which the magnitudes of the electric and magnetic dipoles are equal. This is in clear disagreement with the response expected if the scattering Fig. 3 Force in the x direction versus wavelength (straight line for the analytic expression and circles for the results from full numerical simulations). Note how the force is negative when the electric dipole dominates and positive when the magnetic dipole prevails. The force is equal to zero close to 1700 nm where the contribution from both dipoles is equal force comes only from the Poynting vector, which should be negligible, remaining in any case with the same orientation for every wavelength.
In order to have an idea about the magnitude of the Belinfante's force, we have calculated, on our configuration, the ratio between the curl of the spin force (F B ) and the standard radiation pressure force proportional to the Poynting vector (F P ). For TE (TM) evanescent waves and for magnetic (electric) dipoles, the ratio between these two forces is given by which, in our case, corresponds to F B F P ∼ 0.3.

Full numerical simulation
In order to corroborate this analytic prediction, we have reproduced the experiment using full numerical simulations. By solving Maxwell equations in a mesh, we obtain the electromagnetic field of the two evanescent waves and the scattering of the particle. Then, forces on the probe are calculated by integrating, on a closed surface, the flux of the Maxwell stress energy tensor. The results are plotted in Fig. 3. Note how the outcome from the simulated experiment fits the behavior obtained from the analytic expression, demonstrating that scattering forces on the particle are due to the existence of a Belinfante's momentum, proportional to the curl of the spin density of the light field.

Conclusion
We have analyzed an optical configuration, based on the magnetic spin of TE evanescent waves and the electric spin of TM evanescent waves, capable of elucidating if the momentum transferred to a dipolar particle is the Poynting vector alone or the Poynting vector together with Belinfante's momentum. The measurement is based on the electric and magnetic response of silicon particles in the infrared region of the electromagnetic spectrum. If the momentum transferred to the dipolar particle is given only by the Poynting vector, the result of the proposed experiment should be a negligible force. However, if the canonical momentum is transferred, the result should be a force coming from the curl of the spin density of the light fields (Belinfante's momentum) depending on the difference between the magnetic and the electric response of the probe particle. In particular, for wavelengths between 1400nm and 1700nm the experiment should detect a negative displacement of the particle while, for wavelengths larger than 1700nm crossing the second Kerker condition, the particle should move in the positive direction. The full numerical simulations performed point toward the existence of a component of the radiation pressure coming from the "virtual" Belinfante's momentum.