Negative force on free carriers in positive index nanoparticles

We theoretically demonstrate the reversal of optical forces on free charge carriers in positive refractive index nanostructures. Though optical momentum in positive refractive index materials is necessarily parallel to the local energy flow, reversal of optical momentum transfer can be accomplished by exploiting the geometry and size of subwavelength particles. Using the Mie scattering theory and separation of optical momentum transfers to the bound and free charges and currents, we have shown that metal nanoparticles can exhibit strong momentum transfer to free carriers opposite to the direction of incident electromagnetic waves. This can be explained for small particles in terms of a reversal of Poynting power inside the material resulting in a negative net force on free carriers in small particles. Two-dimensional simulations further illuminate this point by demonstrating the effect of incident wave polarization.

We theoretically demonstrate the reversal of optical forces on free charge carriers in positive refractive index nanostructures. Though optical momentum in positive refractive index materials is necessarily parallel to the local energy flow, reversal of optical momentum transfer can be accomplished by exploiting the geometry and size of subwavelength particles. Using the Mie scattering theory and separation of optical momentum transfers to the bound and free charges and currents, we have shown that metal nanoparticles can exhibit strong momentum transfer to free carriers opposite to the direction of incident electromagnetic waves. This can be explained for small particles in terms of a reversal of Poynting  The reversal of radiation pressure in negative index materials (NIMs) was first predicted by the seminal work of Veselago. 1 NIMs have received significant attention since their realization in 2001 2 due to the potential for interesting physics such as reversal of the Cerenkov effect, 3 perfect lensing, 4 and negative refraction. 5 Veselago's results predict in materials with simultaneously negative permeability and permittivity (µ < 0, < 0), light attraction instead of light pressure, and momentum transmitted into an NIM is antiparallel to the Poynting power. The radiation force on free carriers (F c ) in a material becomes reversed if the index of refraction is negative (n < 0) due to the reversal of wave momentum. 6 Maxwell first applied the electromagnetic wave theory to predict and calculate the radiation pressure of light. 7 Observation of light pressure on reflectors in vacuum validated this prediction. 8,9 After, detailed quantitative measurement of radiation pressure was presented by Poynting in 1905. 10 First observed in 1970, the photon drag effect results from the optical momentum transfer in semiconductors such as germanium or silicon. 11,12 The experiment reveals that when a photon of energy E is absorbed in a semiconductor of a refractive index n, the Minkowski momentum nE/c is transferred to free carriers, where c is the vacuum speed of light. From the momentum conservation, additional force on the bulk material must be directed toward the incident wave. Loudon et al. 13 provided a theoretical description of momentum transfer from electromagnetic wave to charge carriers and host semiconductor for limiting cases introducing the division of momentum transfers to a force on free carriers and a force on the host or bound carriers.
In 1954, Jones and Richards found that the radiation pressure exerted upon a submerged mirror is proportional to the index of reflection n of the surrounding medium by measuring the deflection of refractors in submerging fluids. 14 This experiment failed to determine whether the proportionality of the refraction index was the phase velocity refractive index n = c/v p = ck/ω or the group velocity refractive index n g = c/v g = c∂k/∂ω, where v p and v g are the phase and group velocities, respectively, ω is the frequency, and k is the wavenumber. Jones and Leslie repeated the original experiment with the use of a laser in 1978 and determined that the pressure upon a mirror submerged in a dielectric liquid is proportional to the phase velocity refractive index n. 15 In 2005, Campbell et al. presented a measurement of photon recoil momentum of single atoms in a dispersive medium consisting of a gas of ultra-cold atoms. 16 The experiment revealed the change in momentum of single atoms in the lowest momentum state caused by absorption of photons. It was observed using an atom interferometer that the recoil via optical momentum transferred to the atoms was directly proportional to the classical, macroscopic refractive index n of the gas. The experiment also shows that the presence of all surrounding atoms affect the single photon entering a dilute gas.
In 2015, Zhang et al. group demonstrated remarkable light induced propulsion and rotation of bulk graphene based material. 17 It is stated that the unique structure and properties of the graphene sponge make it not only absorb light but also can emit high velocity electrons which have net momentum towards the light source. This net electron momentum driven toward the laser source can propel the material along the propagation direction of the laser at a momentum larger than the incident photons. Although the incident optical momentum is positive in this case, the net momentum of ejected electrons towards the light source, which facilitates the propulsion, was experimentally shown to be negative. According to the previous discussion, the force on free carriers should be positive in a medium with positive real part of the refractive index n = n R + in I . In this experiment, the observed electron momentum is negative. However, the experiment has been presented with an energy analysis but description as to the mechanism for momentum transfer to the electrons yielding negative momentum. A recent thermal analysis questioned the physical mechanism responsible for the propulsion of the graphene sponge as radiometric forces may be responsible. 18 However, the debate is still ongoing as to the description of this phenomenon. 19 In this letter, we theoretically demonstrate that force on free carriers can be negative in a positive refractive index medium. It is shown that such negative momentum transfer to free carriers accompanies enhanced positive forces on the host material (i.e., bound carriers). We demonstrate this by using metal particles with (n R > 0). Our conclusions reveal that the negative momentum transfer to free carriers within positive refractive index media is highly dependent on the sub-wavelength geometry of the absorber. The justification for our application to metal particles is that it provides a proof-of-concept rooted in well-known material systems with analytical solutions. Whereas application to specific experimental conditions, as we point out, will require very detailed information about the sub-micron geometries involved. For example, we have calculated the optical forces on bound and free carriers for gold and silver nanospheres directly with complex permittivity (ω) = R + i I = | | exp(iφ ) and permeability µ(ω) = µ R + iµ I = | µ| exp(iφ µ ) in vacuum ( 0 , µ 0 ). The total time-averaged optical forceF =F b +F c due to e −iωt dependent harmonic excitation is decomposed into the force on bound carriersF b and force on free carriersF c based on the real and imaginary contributions to the complex permittivity and permeability. 13 The same approach has been previously been previously applied to model the momentum transfer to the isotropic left handed media 6 and absorbing Mie particles. 20 The momentum of light in media and the resulting forces is the subject of a longstanding debate, and attention must be given to the formulation of electrodynamics applied in the theoretical model. [21][22][23] We begin with the field-kinetic stress-energy-momentum tensor resulting from the Chu formulation, which is interpreted as providing the action of the field energy and momentum upon media. 24,25 The force density, 6,20 provides the field-kinetic momentum transfer to free carriers as the wave attenuates in the medium. By applying Poynting's theorem, the force density on free currents can be shown to be written as f c = −ẑ 1 2 Re n c∇ ·S , whereS =Ē ×H * is the complex Poynting power. 20 Thus, the direction of the force on free currents depends on the sign of the index of refraction n. Additionally, the force density on bound currents and charges, provides the remaining momentum transfer to the host material. The total force on a materialF and the total force on freeF c and boundF b carriers result from the integration of the force densities over the entire medium. The Chu force on a particle can be written as below, 6,20,26 where ρ e = −∇ ·P + ρ andJ e = ∂P ∂t +J. So for the non-magnetic particle, Chu force can be written as below,f From the definition (−∇ ·P)Ē + ∂P ∂t × µ 0H is the force on bound carriers and charges,f b . So the force on free carriers and charges becomesf Now from the Minkowski formula the force on a non-magnetic particle is 21,26 When we integrate outside the particle, we get total force but when we integrate just inside the homogeneous particle to find the force of free carriers and charges and in this case, ∇ becomes zero. So thef c = ρĒ +J × µ 0H for Minkowski formulation also gives the same formulation as Chu which is shown in Eq. (5) and this is consistent with the Lorentz force on free carriers. So, the calculated force densityf c is identical between the field-kinetic (Chu) and canonical (Minkowski) formulations, whereas the distribution off b differs. 21 However, the total values for bothF c andF b will be identical due to the necessity of resulting calculations forF being identical. 26 The connection of momentum transfer to bound and free currents is formalized here by applying the momentum conservation theorem via the Maxwell stress tensor. The time-average force on all currents and charges in a volume V enclosed by a surface of area A with area elementda can be determined by applying the divergence of the Maxwell (vacuum) stress tensor. For distinguishing the values ofF c andF b , we may apply the complex Minkowski stress tensor, to calculate the momentum transfer to free carriers enclosed by a surface of area A within a homogeneous medium,F Therefore, the force on free carriersF c is found by integrating the Minkowski stress tensor just inside the surface, and the total forceF is calculated by integrating the stress tensor outside the particle surface. 20 The force on the host (i.e., bound carriers) is determined, then, byF b =F −F c . The total fields presented herein are due to a plane waveĒ inc =xE 0 e ik 0 z incident from free space and are found using Mie theory. [26][27][28] Because the incident power and momentum areẑ directed, we consider the z components of the forces, which we denote simply as F, F b , and F c , where a negative force indicates the −ẑ direction. The dielectric constant of metal particles with radii less than 10 nm is size dependent because in this size range the radius of the particle is comparable with the conduction electron mean free path. 29 It is not our intent, or is this effect dependent on quantum size effects of the nanoparticles or nanorods. As such, we have removed mention of particles less than 10 nm from the calculation of optical force. Figure 1 shows the separation of total Lorentz force, F (black line) into F b (dashed line) and F c (double dashed line) for gold nanoparticles with size range, a = 10-30 nm. Gold nanoparticles have frequency dependent permittivity and which can be negative in some frequency range, but the index of refraction remains positive. Here gold has permittivity, = −5.6141 + 2.2556i at λ = 532 nm with a positive refractive index n = √ r µ r = 0.467 + 2.415i. 30,31 In this case, the F c gives negativity with a positive index of refraction for the gold nanoparticle. However, the F c is positive and F b is negative for absorbing particle with a similar size range of particle when the real part of permittivity is positive with a positive index of refraction (data not shown). So, enhancement of F b is possible for metal particles due to the negative force of free carrier and charges. For example, at a = 15 nm, F c = −1.6575 × 10 −19 nN and F = 9.1624 × 10 −19 nN, so F b = F − F c = 1.0820 × 10 −18 nN. We believe that even stronger F b is possible with a different parameter which can provide large propulsive force to the material. The negative F c depends on the particle's size and geometry. Double dashed line in the inset of Fig. 1 shows the zoomed out plot of F c for gold nanoparticles with a positive refractive index using the same parameters. The F c is minimum when the size of gold particle is 32 nm and becomes positive at 42 nm for gold. Similar results are also found for other materials (e.g., silver), although the results are omitted for brevity.
We further investigate the reason for the negativity of F c by plotting the real part of Poynting power in the x-z plane. According to Eq. (1), the F c depends on the real part of Poynting power. Figure 2 shows the direction of the real part of complex Poynting powerS for gold nanoparticle using the same parameters of Fig. 1 at a = 20 nm when the F c is negative. It is observed that the real part of S becomes negative at the inside of the material and gives negative F c . The derivation off c = −ẑ 1 2 Re n c∇ ·S from Ref. 20 assumes a plane wave solution as a uniform field amplitude. This condition is met in the Rayleigh regime as seen in Fig. 2, but not when the diameter is larger than λ/10. Furthermore, F c < 0 in the Rayleigh regime when the Poynting power in reversed. When the size parameter for the particle is larger than the Rayleigh regime, this effect is not observed. From Eq. (2), the F b depends on the imaginary part of Poynting power. It is also observed that the imaginary part ofS is strongly negative inside the material and gives strong F b when F c is negative. Negative F c is only found at the off-resonance condition. At the resonance condition of gold and silver, the force is stronger but no negative F c is observed (data not shown). We have also shown the frequency dependency of negative force on free carriers and charges, F c in Fig. 3. We have determined separation of total Lorentz F (solid line) into F b (dashed line) and F c (double dashed line) for a gold sphere of radii, a = 20 nm for the wavelength range, λ 0 = 200 − 700 nm. Wavelength dependent dielectric constant is used from tabular data. 30 A gold particle shows resonance at the scattering light, λ 0 = 481 nm when Re( ) = −2. 30,32 It is seen that the F c is negative for a wide frequency range due to the reversal of Poynting power,S inside the material when it crosses the resonance frequency. For the silver particle similar effects have also been found (data not shown). For completeness, that F c can be negative in metal nanostructures besides the spherical metal nanoparticles, we have also calculated the separation of total Lorentz F (solid line) into F b (dashed line) and F c (double dashed line) for the infinite gold cylinder using the same parameters of Fig. 1. We have used normalized force in the results of the infinite cylinder by using incident momentum, P 0 = S 0 /c = 1. Figures 4 and 5 show the incident wave polarization dependency of F c for both TE and TM modes. For TE mode of propagation, theĒ is polarized along the infinite axis of the cylinder, and there is no negative F c found as shown in Fig. 4. But for the case of TM mode of propagation, thē E is polarized along the radius of the cylinder, and the F c shows similar negativity as the spherical particle which is shown in Fig. 5. The size dependency as like spherical particle is also observed in an infinite cylinder.
In summary, the F c can be negative in positive index absorbing nanoparticles, which depends on the size of those nanoparticles. We also demonstrate the polarization effects with infinite cylinders, which indicate a geometric effect along the size effect for realizing negative optical momentum transfer to free carriers in positive refractive index media. If the free carriers come out from the particle, then the total force on the particle is just the F b , and when the F c is negative, by tuning  the parameters, polarization, and geometry of particles, useful strong F b can be achieved to gain propulsive force. 17 The reversal of F c is also can be an indirect explanation for the negative momentum of electron towards a light source and strong propulsive force to the host material. Whether the free carriers come out from the particle or remain in the particle depends on the momentum of incoming light which can be said as the photoelectric effect in metal. Negative F c results can be verified into a (near) vacuum to avoid the probable interference of air when the free carriers come out from the particle. 19 Although there is a lot of research on metal particle/positive index media, but to the best knowledge of authors, the negative force on free carriers is not measured for the positive index of media yet except for the graphene. Our results are theoretical, and we expect additional experimental results similar to the graphene those published by Zhang et al. in