arxiv: v2 [physics.optics] 10 Oct 2018

Size: px

Start display at page:

Download "arxiv: v2 [physics.optics] 10 Oct 2018"

Coral Williams
5 years ago
Views:

1 Dissipation Effect on Optical Force and Torque near Interfaces Daigo Oue 1, 1 Division of Frontier Materials Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka, Japan arxiv: v2 [physics.optics] 10 Oct 2018 (Dated: October 12, 2018) Abstract Fresnel-Snell law, which is one of the fundamental laws in optics and gives insights on the behaviour of light on interfaces, is violated if there exists dissipation in transmitted media. In order to overcome this problem, we extend the angle of refraction from real number to complex number. In this research, we use this complex-angle approach to analyse the behaviour of light on interface between lossy media and lossless media. We reveal that dissipation makes the wavenumber of the light exceed the maximum in the case of lossless interfaces. This is surprising because, in general, dielectric loss only changes the intensity profiles of light, so this excess wavenumber cannot be produced in the bulk even if there exists dielectric loss. Also, anomalous circular polarisation emerges with dissipation. The direction of the anomalous circular polarisation is transverse, whereas without dissipation the direction of circular polarisation has to be longitudinal. We also discuss assist of optical force by the excess wavenumber and generation of optical transverse torque by the anomalous circular polarisation. This novel state of light produced by dissipation will pave the way for new-generation optical trapping and manipulation. daigo.oue@gmail.com; 1

2 I. INTRODUCTION Optical force and optical linear momentum was theoretically revealed by Maxwell [1] and observed by Nichols and Hull [2, 3]. Also, existence of optical torque and optical angular momentum was theoretically proposed by Poynting [4], and experimentally confirmed by Beth [5]. Although these interesting facts that light can push and rotate objects attracted many scientist, there are not so much progress in the study of optical force and torque for a while, because there are no available high-intensity light source for generating large force or torque at that time. Around 1950, laser was invented in Bell laboratories, and this makes large impact on the research of optical force and torque. Ashkin, who also belonged to Bell laboratories at that time, demonstrated manipulation of microparticles by utilizing laser beam [6, 7]. After Ashkin s experiments, many interesting setups for optical manipulation have been proposed. One of them is utilizing optical near field such as evanescent waves and surface plasmon polarion for optical manipulation [8 11]. Near field has peculiar properties such as steep intensity gradient, large wavevector compared to propagating field, and transverse circular polarisation. To clarify the effects of those characteristics of near field on optical force and torque, many theoretical and experimental studies have been conducted[12 19]. However, the dissipation, i.e. dielectric or magnetic loss, effects on the field, especially near interfaces, have never been discussed before. Electromagnetic fields also decay due to dissipation (without any total internal reflection), and in dissipative media they can have near-field-like characteristics such as intensity gradient. In this paper, we use complex-angle approach to analyse the behaviour of optical field near interfaces with dissipation, and we also analyse how the dissipation contributes to optical force and torque. II. ELECTROMAGNETIC FIELD NEAR INTERFACES WITH DISSIPATION We utilize complex-angle apploach [14] for calculation of electromagnetic field near dissipative interfaces. FIG. 1 shows the situation we consider in this paper. We consider interfaces between dissipative media and non-dissipative media. Field is incident on the interface from the lower side, and we take θ 1 for the angle of incidence. In dissipative media, the refractive index must be complex value ( ñ 2 R C ). ñ 2 = n 2 +iκ 2. (1) 2

$Let the angle of incidence θ 1 and the angle of refraction θ 2.$

3 FIG. 1. Setup we discuss in this paper. We consider interface between non-dissipative media (n 1 ) and dissipative media (ñ 2 = n 2 +iκ 2 ). Let the angle of incidence θ 1 and the angle of refraction θ 2. Considering Snell law, n 1 sinθ 1 = ñ 2 sin θ 2 (2) In this equation, the left hand side is real (n 1,sinθ 1 R), while ñ 2 has the imaginary part. Thus, sin θ ( 2 must be complex sin θ 2 R C ), and the angle of refraction must be a complex value, θ 2 = θ 2 +iψ 2 (ψ < 0), (3) which satisfies the following simultaneous equations (modified Snell law): { Re(ñ2 sin θ 2 ) = n 1 sinθ 1, (4a) Im(ñ 2 sin θ 2 ) = 0. Once we obtain the angle of refraction, we can calculate the explicit expressions of the transmitted field by rotating +z-propagating plane wave (5) towards the directions of refraction. (4b) E = E 0 exp(ik r) (5) Here, E p 0 E 0 = E s,k = 0. (6) 0 ñ 2 k 0 k 0 is the wavenumber in vacuum. For the rotation, we use cosθ 0 sinθ R(θ) = (7) sinθ 0 cosθ 3

4 By substituting the complex-angle of refraction θ 2 into the rotation matrix, we get the complex-angle rotation matrix which gives the wavevector of electromagnetic field in dissipative media: k k = R( θ 2 )k = k 0 ( k +i η), (8) where n 2 sinθ 2 coshψ 2 κ 2 cosθ 2 sinhψ 2 k = 0, (9) n 2 cosθ 2 coshψ 2 +κ 2 sinθ 2 sinhψ 2 0 η = 0. (10) n 2 sinθ 2 sinhψ 2 +κ 2 cosθ 2 coshψ 2 The imaginary part k 0 η represents the decay, and the real part k 0 k the propagation. Note that, from the second equation of the modified Snell law, n 2 cosθ 2 sinhψ 2 +κ 2 sinθ 2 coshψ 2 = 0, (11) so the x-component of the imaginary part of the wavevector η x vanishes. FIG. 2 shows the iso-frequency curves (the curve drawn by the real part of the wavevector for fixed k 0 ). The purple curve is the curve with the presence of dissipation (ñ 2 = i), and the dashed black curve is without dissipation (ñ 2 = i). It can be said that the curve with dissipation covers the non-dissipative curve, which means that higher (spatial) frequency field is produced by dissipation and that this could enlarge optical force on particle in the field. We can also observe that the larger angle of incidnece we have, the more excess wavenumber is produced. Likewise, we can calculate the field vector: T p E p E 0 Ẽ0 = R( θ 2 ) T s E s 0 T p E p (cosθ 2 coshψ 2 isinθ 2 sinhψ 2 ) = T s E s (12) T p E p (sinθ 2 coshψ 2 +icosθ 2 sinhψ 2 ) 4

5 FIG. 2. Excess wavevector is produced by dissipation. Dashed curve is the iso-frequency curve without the presence of dissipation (ñ 2 = i), while purple curve is the iso-frequency with the presence of dissipation (ñ 2 = i). Here, T p and T s are Fresnel coefficients. Finally, we get the explicit representation of the refracted field near interface with dissipation: ( ) Ẽ = Ẽ0exp i k r = Ẽ0exp ( ik 0 k r ) exp( k0 η r). (13) In this paper, we are interested in p-polarisation since it has non trivial longitudinal field vector, so below this we set E s = 0 and for simplicities set E p = 1 and n 1 = 1, but it is straightforward to discuss general case (E s,e p 0, and n 1 1) using our approach. We can see that there is phase difference between the transverse component (z) of field vector and the longitudinal component (x). That is because the conditions of continuity for the tangential component of field and for the normal component are different, and this induces phase difference between the two components, and causes rotation of the field. To describe degree of circular polarisation and direction of the field rotation, we can use a psuedo-vector s = g 4ω Im(E E), (14) where g = 1/4π is a Gaussian-unit factor. For our field, we have s = g 2ω sinhψ 2coshψ 2 exp( 2k 0 η r)e y. (15) 5

6 FIG. 3. Anomalous circular polarisation (transverse y-component of s) emerges with dissipation (ñ 2 = i). Here, e y is the unit vector in the direction of y. We can observe if we flip the sign of θ 2, then from (4b) the sign of ψ 2 is also flipped, and thus the direction of the psuedo-vector s is flipped. This is one kind of spin-momentum locking [20], in which the direction of circular polarisation and the direction of propagation is tied to each other. In FIG. 3, we plot the transverse (y) component of the psuedo-vector. This implies that the field is rotating in the transverse direction and that we can exert transverse optical torque on particle in the field. III. OPTICAL FORCE AND TORQUE NEAR DISSIPATIVE INTERFACES For the calculation of optical force and torque, particle which is small enough compared to the wavelength of the field can be regarded as a point dipole. Time-averaged optical force and torque exerted on a point dipole is given by[13, 17, 21] F = (d e )E +(d m )H+ 1 d c e B 1 d c m D, (16) T = d e E +d m H. (17) Here, E and H are real fields associated with complex fields, E = Re(Ee iωt ) and H = Re(He iωt ). Please note that we here do not consider higher order contribution, which could cause electric-magnetic dipolar interaction force [21] and torque [22, 23], but focus on the fundermental order by the dipole approximation. These complex fields satisfy monochro- 6

7 FIG. 4. Excess-wavenumber-assisted optical force. Upper figure is the plot of scattering forces excerted on gold nanosphere (100 nm in diameter) in the case without dissipation (ñ 2 = 1.3+0i). Purple is F scattering,x, and green is F scattering,z. Lower figure is the case with dissipation (ñ 2 = i). Likewise, purple is scattering force in x direction, and green is z direction. In both cases, the excitation wavelength is 650 nm. matic Maxwell equations: E = 0, (18) H = 0, (19) E = i ω µh, c (20) H = i ω εe. c (21) Complex flux densities, D and B, are characterized by permittivity ε and permiability µ: D = εe, B = µh, and these give the real fields: D = Re(De iωt ), B = Re(Be iωt ). We also use complex dipole moments d e, d m to give real dipole moments: d e (r,t) = Re(d e (r)e iωt ), d m (r,t) = Re(d m (r)e iωt ). Electric and magnetic dipoles are characterized by electric polarizability α e and by magnetic polarizability α m : d e = α e E, d m = α m H. InFIG.4, wecomparescatteringforce(radiationforce)on100nmgoldnanosphereby650 nm excitation with or without dissipation. Upper figure is the case without any dissipation (ñ 2 = 1.3+0i), and the lower is with dissipation (ñ 2 = i). We can naively say that dissipation assists optical force in both x and z direction. At large θ 2, the scattering force reach 0. This is because the larger angle of refraction is, the smaller Fresnel coefficient is (When the angle of incidence is π/2, the transmission coefficient is zero). 7

8 FIG. 5. Dissipation assist effect. Solid curves are the optical force on 100 nm gold nanosphere with dissipation (ñ 2 = i), and dashed curves are the case without dissipation (ñ 2 = 1.3+0i). Purple curves are the x-components of the force, while green curves are the z-components. It is clear that in both x and z direction, the optical radiation force is assisted. FIG. 6. Transverse torque excerted by anomalous circular polarisation on 100 nm gold nanosphere. The excitation wavelength is 650 nm. Solid curve is y-component of the optical torque T y with dissipation (ñ 2 = i), while the dashed curve is without any dissipation (ñ 2 = 1.3+0i). It is obvious that there are non-zero transverse torque with dissipation, whereas no torque without dissipation since without dissipation the field does not rotate in transverse direction. To calarify the excess wavenumber effect, we plot optical force on particle per unit power in FIG. 5. We can confirm in both x and z directions optical force per intensity is assisted by dissipation. FIG. 6 shows transverse optical torque on the gold nanosphere induced by anomalous circular polarisation, which cannot be generated without dissipation. As the anomalous 8

9 circularpolarisationvanishesatθ 2 0,thetransverseopticaltorquevanishes. Thisconfirms that the anomalous polarisation causes the transverse torque. IV. CONCLUSION To sum up, we utilize complex-angle approach for calculating electromagnetic field near interfaces with dissipation, and revealed dissipation-induced extraordinary behaviours of the field: production of excess wavevector and generation of anomalous transverse circular polarisation. We also studied what kind of effects these extraordinary behaviours cause on optical force and torque. Excess wavevector assisted optical force, and anomalous circular polarisation generated transverse optical torque. Since the effects disscussed in this paper can be caused only by adding some dissipation, it is easy to employ in experiments, and they will increase degrees of freedom of optical trapping and manipulation. [1] J. C. Maxwell, Philosophical transactions of the Royal Society of London 155, 459 (1865). [2] E. F. Nichols and G. F. Hull, in Proceedings of the American Academy of Arts and Sciences, Vol. 38 (JSTOR, 1903) pp [3] E. F. Nichols and G. F. Hull, Physical Review (Series I) 17, 26 (1903). [4] J. H. Poynting et al., Proc. R. Soc. Lond. A 82, 560 (1909). [5] R. A. Beth, Physical Review 50, 115 (1936). [6] A. Ashkin, Physical review letters 24, 156 (1970). [7] A. Ashkin, J. M. Dziedzic, J. Bjorkholm, and S. Chu, Optics letters 11, 288 (1986). [8] S. Kawata and T. Sugiura, Optics letters 17, 772 (1992). [9] S. Kawata and T. Tani, Optics letters 21, 1768 (1996). [10] K. Okamoto and S. Kawata, Physical review letters 83, 4534 (1999). [11] G. Volpe, R. Quidant, G. Badenes, and D. Petrov, Physical review letters 96, (2006). [12] P. Chaumet and M. Nieto-Vesperinas, Physical Review B 61, (2000). [13] A. Y. Bekshaev, Journal of Optics 15, (2013). [14] A. Y. Bekshaev, K. Y. Bliokh, and F. Nori, Optics express 21, 7082 (2013). 9

10 [15] A. Canaguier-Durand, A. Cuche, C. Genet, and T. W. Ebbesen, Physical Review A 88, (2013). [16] A. Canaguier-Durand and C. Genet, Physical Review A 89 (2014), /physreva [17] K. Y. Bliokh, A. Y. Bekshaev, and F. Nori, Nature communications 5, 3300 (2014). [18] M. Antognozzi, C. Bermingham, R. Harniman, S. Simpson, J. Senior, R. Hayward, H. Hoerber, M. Dennis, A. Bekshaev, K. Bliokh, et al., Nature Physics 12, 731 (2016). [19] A. Y. Bekshaev, Journal of Optics 20, (2018). [20] T. Van Mechelen and Z. Jacob, Optica 3, 118 (2016). [21] M. Nieto-Vesperinas, J. Sáenz, R. Gómez-Medina, and L. Chantada, Optics express 18, (2010). [22] M. Nieto-Vesperinas, Optics letters 40, 3021 (2015). [23] M. Nieto-Vesperinas, Physical Review A 92, (2015). 10

Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Electro Dynamic

Electrodynamics I Final Exam - Part A - Closed Book KSU 2005/12/12 Name Electro Dynamic Instructions: Use SI units. Short answers! No derivations here, just state your responses clearly. 1. (2) Write an