arxiv: v2 [physics.optics] 14 Feb 2018

Transcription

1 Goos-Hänchen and Imbert-Fedorov Shifts for Airy Beams Marco Ornigotti Institut für Physik, Universität Rostock, Albert-Einstein-Straße 23, Rostock, Germany (Dated: Compiled October 10, 2018) In this work, I present a full analytical theory for the Goos-Hänchen and Imbert-Fedorov shifts experienced by an Airy beam impinging on a dielectric surface. In particular, I will show how the decay parameter α associated with finite energy Airy beams is responsible for the occurence of giant angular shifts. A comparisons with the case of Gaussian beams is also discussed. arxiv: v2 [physics.optics] 14 Feb 2018 Accelerating beams, i.e., beams of light that propagate along curved trajectories in free space without any external potential present, attracted a lot of interest in the last years. The forefather of this new class of optical fields is undoubtledy the Airy beam, first introduced in quantum mechanics by Berry and Balazs in 1979 [1] and then brought to optics by by Siviloglou and co-workers in 2007 [2, 3]. This particular solution of the paraxial equation, in fact, propagates along a parabolic trajectory in free space, due to the fact that it experiences a nonzero tangential acceleration in its own comoving reference frame. Driven by the intriguing properties possessed by Airy beams, in the last decade many different classes and types of accelerating beams were proposed, each with unique features, such as the ability of propagating along curved [4, 5], spiraling [6] or arbitrary [7, 8] trajectories in space, route particles along non-conventional paths [9], or generate curved plasma channels [10], to name a few. Airy beams, in particular, were extensively studied [11], including a detailed analysis of their reflection from dielectric surfaces close to total internal reflection and Brewster incidence [12]. The occurrence of selective reflection was also reported [13]. Reflection of light upon an interface is a daily witnessed phenomenon, which can be easily understood and explained using Snell s law [14]. This law, however, holds exactly only when the light impinging upon a reflecting surface can be described within the geometrical optics approximation. When the wave properties of the impinging field cannot be neglected, like in the case of an optical beam for example, nonspecular reflection phenomena typically occur. Among these effects, the most famous are certaintly the so-called Goos-Hänchen [15 17] and Imbert-Fedorov [18 28] shifts, the former occurring in the plane of incidence and the latter in the plane perpendicular to it. These phenomena have been extensively studied in the past decades for a broad class of beam configurations [29 34] and interfaces [35 41]. Moroever, the analogy between beam shifts and the quantum mechanical weak measurements has also been proposed [42 46] and used as a way to get strongly amplified beam shifts [47, 48]. A comprehensive review on the beam shift phenomena can be found in Ref. [49]. For the case of Airy beams, intermittend Goos-Hänchen shift from nonlinear marco.ornigotti@uni-rostock.de FIG. 1. (Color online) Schematic representation of the geometry of the problem. The laboratory frame {ˆx, ŷ, ẑ} is attached to the reflecting surface (blue box), with the z-axis pointing downwards. The incident ({ˆx i, ŷ i, ẑ i}) and reflected ({ˆx r, ŷ r, ẑ r}) reference frames are instead shown by the green planes. θ is the angle of incidence. surfaces has been predicted and numerically verified in Ref. [50]. However, a comprehensive theory of Goos- Hänchen and Imbert-Fedorov shifts for Airy beams has not been developed yet. In this Letter, therefore, I will present a complete theory of Goos-Hänchen (GH) and Imbert-Fedorov (IF) shifts for Airy beams impinging upon a dielectric interface. In particular, I will show that for the case of finite energy Airy beams, giant angular GH and IF shifts occur, as a consequence of the finite width of the Airy beam spectrum. To begin with, let us consider a monochromatic, paraxial electric field impinging upon a dielectric surface, as depicted in Fig. 1. For later convenience, we can define three reference frames: the laboratory frame {ˆx, ŷ, ẑ} attached to the dielectric interface and oriented with the ẑ axis pointing towards the dielectric, and two auxiliary frames {ˆx k, ŷ k, ẑ k } attached to the incident (k = i) and the reflected (k = r) field, respectively. The impinging electric field can be then described, in the incident frame,

2 2 by the following Fourier representation: E inc (r) = 2 d 2 Kê λ (U, V ; θ)a λ (U, V )e i(uxi+v Yi+W Zi), λ=1 (1) where d 2 K = dudv, λ is a polarization index (λ = 1 corresponds to p-polarization and λ = 2 corresponds to s-polarization). In the expression above, U, V and W = 1 U 2 + V 2 1 (U 2 + V 2 )/2 are the dimensionless components of the k-vector in the incident frame, which are connected to the k-vector in the laboratory frame by the relation k = k 0 (U ˆx i + V ŷ i + W ẑ i ) = k xˆx + k y ŷ + k z ẑ [49]. Moreover, A λ (U, V ) = α λ (U, V ; θ)a(u, V ) is the vector spectral amplitude, with A(U, V ) being the angular spectrum of the incident field and α λ (U, V ) the polarisation functions given by [51] α λ (U, V ) = ˆf ê λ (k), (2) with ê λ (k) being the polarisation vectors attached to the single plane wave components of the impinging beam, whose definition is given as ê 1 (k) = (ê 2 k])/ ê 2 k and ê 2 (k) = (ẑ k)/ ẑ k. The unit vector ˆf = a pˆx i +a s e iη ŷ i (with a p, a s, η R) describes the polarization state of the beam in the incident frame. Upon reflection, each plane wave component of the Fourier spectrum in (1) gets reflected according to Fresnel s law, i.e., ê λ (k)e ik r r λ (k)ê λ ( k)e i k r, (3) where r λ (k) are the Fresnel reflection coefficients [14] and k = k 2ẑ (ẑ k) is determined by the law of specular reflection [52]. The expression for the electric field E refl (r) in the reflected frame can be the obtained from (1) upon substituting ê λ (k) ê λ ( k) and k i k r = U ˆx + V ŷ + W ẑ, thus obtaining E refl (r) = 2 d 2 K E(U, V ; θ)e i( UXr+V Yr+W Zr). λ=1 (4) where E λ (U, V ; θ) = ê λ ( U, V ; π θ)r λ (U, V )α λ (U, V ; θ)a(u, V ) contains informations about both the reflection coefficients and the angular spectrum of the impinging beam. According to the standard theory of beam shifts [33, 53], to evaluate the GH and IF shifts one first needs to calculate the centroid of the intensity distribution of the reflected field, i.e., X = X r ˆx r + Y r ŷ r, where X r = d 2 R X r I(R, z) d2 R I(R, z) being d 2 R = dx r dy r, I(R, z) = E refl (R, z) E refl (R, z) is the intensity distribution of the reflected field in the reflected frame.by substituting (4) in the above expression, after some straightforward calculations we obtain the following expresison for the centroid X r = 1 { N Im d 2 K E (U, V ; θ) E(U, V ; θ) U Z d 2 K U } W E(U, V ; θ) 2, (5) where N = d 2 K E(U, V ; θ) 2, and d 2 K = dudv. An analogue expression (with the substitution / U / V and U/W V/W ) is valid for Y r. Then, the spatial ( ) and angular (Θ) GH and IF shifts have the following form: GH = X r z=0, Θ GH = X r, z (6a) IF = Y r z=0, Θ IF = Y r z. (6b) In order to express the GH and IF shifts for Airy beams in a compact form, it is useful to first recall the expressions for the GH and IF shifts for a Gaussian beam impinging upon a dielectric surface. These can be obtained by setting A(U, V ) in (1) to be a Gaussian function and performing the integrations in (5) by expanding the integrand functions in the numerator and denominator up to the first order in U and V (by virtue of the paraxial approximation) for the spatial shifts and the numerator up to the second order for the angular shifts. After some manipulations, we get the following well known expressions [33]: (g) GH = 1 k 0 ( w p φ p Θ (g) GH = 1 k 0 z R (g) IF = 1 k 0 { cot θ + w s ( ln R p w p ) φ s, (7a) ), (7b) + w s ln R s [ wp a 2 s + w s a 2 p cos η a p a s + 2 ]} w p w s sin(η + φ s φ p ), (7c) [ ] Θ (g) IF = 1 wp a 2 s w s a 2 p cos η cot θ, (7d) k 0 z R a p a s where the superscript (g) stands for Gaussian beam, k 0 = 2π/λ and z R is the Rayleigh range of the beam. In the above expresions, R λ and φ λ are the modulus and phase of the reflection coefficients, i.e., r λ = R λ e iφ λ [51] and w λ = a 2 λ R2 λ /(a2 pr 2 p + a 2 sr 2 s) is the fractional energy contained in each polarisation. For the case of Airy beams, instead, the angular spectrum appearing in (1) has the following form [1]: A(U, V ) = e i 3 (U 3 +V 3), (8) which corresponds, in the plane z = 0, to the 2D pure (infinite energy) Airy beam E(x, y, z = 0) = Ai(x)Ai(y). To calculate the beam shifts for such field distribution, therefore, one should use (8) to calculate the integrals

3 3 appearing in (5). In doing so, however, one needs to deal with integrals of the type I(n, m) du dv U n V m A(U, V ) 2, (9) where n, m N. With the choice of angular spectrum as the one given by (8), it is not difficult to show that I(n, m) =. This result prevents one to calculate the integrals appearing in (5), and, moreover, their ratio. To overcome this problem, we can introduce a regularisation parameter in the integrals I(n, m), in such a way that they give finite results, thus allowing a correct evaluation of (5). Then, by taking the limit of the regularised result when the regularisation parameter goes to zero, we can obtain the correct result for pure Airy beams. Conceptually, this regularisation technique is the same one employed for the calculation of Goos-Hänchen and Imbert-Fedorov shifts for Bessel beams [34]. A convenient regularisation condition for Airy beams, is the one given in Ref. [2], which corresponds to the so-called finite energy Airy beams. This regularisation consists in adding a proper exponential factor to the expression of the Airy beam, namely E(x, y, z = 0) = Ai(x)Ai(y) exp [α(x + y)]. The regularisation parameter α > 0 ensures a containment of the tails of the Airy function, thus allowing such beams to carry finite energy and be then physically realisable [2]. The angular spectrum of such a regularised beam can be calculated using the integral representation of the Airy function [54] Ai(x) = 1 2π and the identity [54] This gives the following result: dz e i ( z 3 3 +xz ), (10) dx e xy Ai(x) = e y3 3. (11) A(U, V ) = e α2 (U 2 +V 2 )+ α3 3 e i 3[(U 3 +V 3 ) 3α 2 (U+V )]. (12) In the limit α = 0, the above angular spectrum reduces to (8). Notice, moreover, that A(U, V ) 2 = exp [ 2α 2 (U 2 + V 2 )] exp [2α 3 /3]. Therefore, the integrals in (5) become Gaussian integrals. This suggests that the GH and IF shifts for regularised Airy beam will have a form similar to the one for Gaussian beams. The beam shifts for an Airy beam can be then calculated using the regularised angular spectrum given by (12) (which, then, corresponds to finite energy Airy beams) and then take the limit of the results obtained in this way when α 0 to obtain the correspondent result for pure (infinite energy) Airy beams. If we now calculate the integrals in (5) for the case of regularised Airy beams, we obtain the following result: GH = (g) GH Γ, IF = (g) IF + Γ, Θ GH = 1 2α 2 Θ(g) GH, Θ IF = 1 2α 2 Θ(g) IF, (13a) (13b) where Γ = α 2 /k 0. This is the main result of this Letter. Before going any further, it is worth commenting on the results pre- sented above. To understand them, one needs first to recall, that when calculating the beam shifts, both the numerator and the denominator of (5) are Taylor expanded in U and V up to the first order. This means that only terms proportional toi(0, 0), I(1, 0), and I(0, 1) will contribute to define the shifts. For the case of the spectrum in (8), moreover, since A(U, V ) 2 is a Gaussian function of U and V, I(1, 0) = I(0, 1) = 0, and therefore only terms proportional to I(0, 0) will contribute in determining the spatial and angular shifts in (13). This justifies the fact that the shifts for regularised Airy beams looks very similar to the ones for Gaussian beams. For the angular shifts, instead, the presence of the factor U/W = U/ 1 U 2 V 2 in (5) requires the integrals to be expanded up to second order in U and V. Therefore, the result will be proportional to I(2, 0)/I(0, 0) 1/α 2. To understand the presence of the Γ-factor in the spatial shifts, instead, one should note that when calculating the derivative in (5), the term A(U, V )/ U (or its V -counterpart) gives a contribution of the form ( 2α 2 U + iu 2 iα 2 ). Since only O(U 2, V 2 ) terms must be taken into account, and Gaussian integrals with odd powers vanish, only the term proportional to α 2 will survive, thus giving the extra Γ term in Eqs. (13a). Notice, that this term comes from the extra linear phase exp [ iα(u + V )] present in the spectrum because of the regularisation. As α 0, therefore, this contribution vanishes, as the angular spectrum for a pure Airy beam does not contain such a linear phase term. For infinite energy Airy beams, therefore, the spatial Goos-Hänchen and Imbert-Fedorov shifts resemble the ones for a Gaussian beam. For finite energy Airy beams, the decay factor α affects both the spatial and angular shifts. For the spatial shifts, it accounts for an additive quantity ±Γ. For the angular shifts, on the other hand, the decay factor appears as a multiplicative constant in the angular shifts, thus resulting in an amplification or suppression of the shift (with respect to the case of a Gaussian beam), depending on whether α < 1 or α > 1, respectively. In practical cases, according to Refs. [2, 3], α is chosen to be a small number, of the order of This choice is motivated by the necessity of having a rapidly enough decaying function, which ensures the Airy beam to carry finite energy. In this case, therefore, Γ 1 and the spatial shifts are the ones of a Gaussian beam. This is in accordance with the numerical results obtained in in Ref. [13], where the predicted spatial GH shift is essentially the one of a Gaussian beam [see for example Fig. 5(c) of Ref. [13]]. For the angular shifts, instead, we have Θ GH,IF /Θ (g) GH,IF = 1/2α Finite energy Airy beams thus display giant angular GH and IF shifts. The amplitude of such angular shifts can be then controlled by suitably tuning the decay paramenter α. As a last remark, it is also interesting to explore the limiting case of (13) when α 0, corresponding to

4 s polarization where ε λ = (1/2)[(R p) 2 + R p R p R 2 s + R p R p cot θ (R 2 p R 2 s) csc 2 θ] is the second-order correction to the total energy carried by the beam [51], and the prime indicates derivation with respect to θ p polarization FIG. 2. Amplification factor χ (pure) for a pure Airy beam for an air-glass interface (n glass = 1.5) for p- (blue, lower, solid line) and s-polarisation (red, upper, solid line). As it can be seen, while χ (pure) p < 1 for θ < θ B, χ s (pure) > 1 for every incident angle. The black, dashed vertical line indicates the position of the Brewster angle, which in this case is θ B 56, 3. GH and IF shifts for nondiffracting, infinite energy Airy beams. For the case of spatial shifts, this limit is immediate, as Γ 0 when α 0, thus meaning that the spatial shifts for infinite energy Airy beams are the same as for a Gaussian beam. For the angular shifts, on the other hand, the presence of α at the denominator makes the limit of Eqs. (13b) diverge. This is due to the fact that, while the integral in the numerator of the angular part of the shifts is expanded up to second order in U and V, the denominator is only expanded up to first order in U and V. A similar problem, for example, occurs as well when dealing with bean shifts close to Brewster incidence [51]. To solve this problem, it is then sufficient to calculate the normalisation factor appearing in (5) up to second order precision in U and V. By doing (g) so, we obtain Θ GH = χ Θ GH, and Θ (g) IF = χ Θ IF, (where Θ (g) GH,IF are obtained from Eqs. (7) by substituting w λ with w λ = a 2 λ R2 λ /[a2 p(rp+ε 2 p )+a 2 s(rs 2 +ε s )]). The second order corrected amplification factor χ is given by χ = a 2 p(r 2 p + ε p ) + a 2 s(r 2 s + ε s ) a 2 p(2r 2 pα 2 + ε p ) + a 2 s(2r 2 sα 2 + ε s ), (14) The expression above for the amplification factor represents the correct form to use in the vicinity of Brewster incidence, and allows us to correctly take the limit α 0. For linear polarisation, in particular, we get χ (pure) λ = 1 + Rλ 2/ε λ, (λ = {p, s}). In particular, if we consider p-polarisation, at the Brewster angle we have R p = 0 and we therefore obtain χ p (pure) = 1. As can be seen from Fig. 2, moreover, the amplification factor is greater than one for every incident angle only for s-polarisation (red, upper, solid line), while for p- polarisation it is greater than one only for θ > θ B. In conclusion, I have calculated GH and IF shifts for both finite and infinite energy Airy beams. For the former, the decay factor α is responsible for the occurrence of a giant angular GH and IF shift, which, for realistic values of α are predicted to be approximately a thousand times bigger than the correspondent shifts for a Gaussian beam. Moreover, finite energy Airy beams possess nonzero GH and IF shifts even when no total internal reflection takes place. This is in contrast with usual results for Gaussian beams [39 41], where GH 0 only close to total internal reflection. This residual spatial shift, however, is significantly smaller than the wavelength, as Γ α 2 λ λ for realistic values of α. For the case of infinite energy carrying Airy beams, instead, the angular shifts are enhanced by a quantity χ, which is independent on the form of the beam and selectively amplifies s-polarisation at all incident angles. The spatial shifts GH,IF, on the other hand, are unaffected by the particular structure of Airy beams and correspond to the ones for Gaussian beams. The author would like to thank A. Szameit and L. Ornigotti for stimulating discussions and careful proofreading of the manuscript. [1] M. V. Berry and N. L. Balazs, Am, J. Phys. 47, 274 (1979). [2] G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, Phys. Rev. Lett. 99, (2007). [3] G. A. Siviloglou, and D. N. Christodoulides, Opt. Lett. 32, 979 (2007). [4] M. A. Bandres Opt. Lett., 34, 3791 (2009). [5] I. Kaminer, R. Bekenstein, J. Nemirovsky, and M. Segev, Phys. Rev. Lett. 108, (2012). [6] C. Vetter, T. Eichelkraut, M. Ornigotti, and A. Szameit, Phys. Rev. Lett. 113, (2014). [7] E. Greenfield, M. Segev, W. Walasik, and O. Raz, Phys.Rev. Lett., 106, (2011). [8] J. Zhao, I. D. Chremmos, D. Song, D. N. Christodoulides, N. K. Efremidis, and Z. Chen, Scientific Reports, 5, (2015). [9] J. Baumgartl, M. Mazilu, and K. Dholakia, Nat. Photonics 2, 675 (2008).

5 5 [10] P. Polynkin, M. Kolesik, J. V. Moloney, G. A. Siviloglou, D. N. Christodoulides, Science 10, 5924 (2009). [11] Y. Hu, G. A. Siviloglou, P. Zhang, N. K. Efremedis, D. N. Christodoulides, and Z. Chen, Self-accelerating Airy Beams: Generation, Control, and Applications, in Nonlinear Photonics and Novel Optical Phenomena (Springer, 2012). [12] Y. Li, L. Li, Y. X. Lu, X. Zhao, K. Xu, Y. Zhang, and Y. Zhang, Opt. Express 21, 8311 (2013). [13] I. D. Chremmos, and N. K. Efremidis, J. Opt. Soc. Am. A 29, 861 (2012). [14] M. Born and E. Wolf, Principles of Optics 7th edn (New York: Cambridge University Press) [15] F. Goos and H. Hänchen, Ann. Phys. 1, 333 (1947). [16] K. Artmann, Ann. Phys., 2, 87 (1948). [17] M. McGuirk and C. K. Carniglia, J. Opt. Soc. Am. 67, 103 (1977). [18] F. I. Fedorov, Dokl. Akad. Nauk SSSR, 105, 465 (1955). [19] C. Imbert, Phys. Rev. D 5, 787 (1972). [20] F. Pillon, H. Gilles and S. Girard, Appl. Opt. 43, 1863 (2004). [21] H. Schilling, Ann. Phys. (Berlin), 16, 122 (1965). [22] M. A. Player, J. Phys. A: Math. Gen. 20, 3667 (1987). [23] V. G. Fedoseyev, J. Phys. A: Math. Gen. 21, 2045 (1988). [24] V. S. Liberman, and B. S. Zeldovich, Phys. Rev. A 46, 5199 (1992). [25] K. Y. Bliokh, and Y. P. Bliokh, Phys. Rev. Lett. 96, (2006). [26] K. Y. Bliokh, and Y. P. Bliokh, Phys. Rev. E 75, (2007). [27] A. Aiello, and J. P. Woerdman, Opt. Lett. 33, 1437 (2008). [28] K. Y. Bliokh, I. V. Shadrivov, and Y. S. Kivshar, Opt. Lett. 34, 389 (2009). [29] P. T. Leung, C. W. Chen, and H. P. Chiang, Opt. Commun. 276, 206 (2007). [30] M. Merano, A. Aiello, G. W. t Hooft, M. P. van Exter, E. R. Elier, and J. P. Woerdman, Opt. Expr. 15, (2007). [31] T. Tamir, J. Opt. Soc. Am. A 3, 558 (1986). [32] G. D. Landry, and T. A. Maldonado, Appl. Opt. 35, 5870 (1996). [33] M. Ornigotti, and A. Aiello, J. Opt. 15, (2013). [34] A. Aiello, and J. P. Woerdman, Opt. Lett. 36, 543 (2011). [35] S. Kozaki, and H. Sakurai, J. Opt. Soc. Am. 68, 508 (1978). [36] D. Golla, and S. Dutta Gupta, arxiv: v1. [37] M. Merano, N. Hermosa, J. P. Woerdman, and A. Aiello, Phys. Rev. A 82, (2010). [38] S. Grosche, M. Ornigotti, and A. Szameit, Opt. Express 23, (2015). [39] M. Araújo, S. De Leo, and G. Maia, Phys. Rev. A 93, (2016). [40] O. Santana, S. Carvalho, S. De Leo, and L. de Araujo, Opt. Lett. 41, 3884 (2016). [41] M. Araújo, S. De Leo, and G. Maia, Phys. Rev. A 95, (2017). [42] M. R. Dennis, and J. B. Götte, New J. Phys. 14, (2012). [43] F. Töppel, M. Ornigotti, and A. Aiello, New J. Phys. 15, (2013). [44] O. Hosten and P. Kwiat, Science 319, 787 (2008). [45] Y. Qin, Y. Li, X. Feng, Y. F. Xiao, H. Yong, and Q. Gang, Opt. Express 19, 9636 (2011). [46] Y. Gorodetski, K. Y. Bliokh, B. Stein, C. Genet, N. Shitrit, V. Kleiner, E. Hasman, and T. W. Ebbsen, Phys. Rev. Lett. 109, (2012). [47] G. Jayaswal, G. Mistura, and M. Merano, Opt. Lett. 38, 1232 (2013). [48] G. Jayaswal, G. Mistura, and M. Merano, Opt. Lett. 39, 2266 (2014). [49] K. Y. Bliokh, and A. Aiello, J. Opt. 15, (2013). [50] P. Chamorro-Posada, J. Sanchez-Curto, A. B. Aceves, and G. S. McDonald, Opt. Lett (2014). [51] A. Aiello and J. P. Woerdman, arxiv: v2. [52] R. F. Gragg, Am. J. Phys. 56, 1092 (1988). [53] A. Aiello, New J. Phys. 14, (2012). [54] O. Vallé, and M. Soares, Airy Functions and Applications to Physics (World Scientific, 2004).