arxiv: v1 [math.ap] 17 Jan 2011

Transcription

1 Quantifying momenta through the Fourier transform B. M. Rodríguez-Lara, Centre for Quantum Technologies, National University of Singapore, 2 Science Drive 3, Singapore arxiv:0.32v [math.ap] 7 Jan 20 Abstract Integral transforms arising from the separable solutions to the Helmholtz differential equation are presented. Pairs of these integral transforms are related via Plancherel theorem and, ultimately, any of these integral transforms may be calculated using only Fourier transforms. This result is used to evaluate the mean value of momenta associated to the symmetries of the reduced wave equation. As an explicit example, the orbital angular momenta of plane and elliptic-cylindrical waves is presented. cqtbmrl@nus.edu.sg

2 I. INTRODUCTION The orbital angular momentum (OAM) of light has become a fundamental concept in optics [, 2]. The fact that light may carry orbital angular momentum which is not related to its polarization state has been shown both theoretically [3] and experimentally [4]. The optical tweezers toolbox now includes experiments involving the transfer of OAM from paraxial beams [5, 6] and vector fields [7] to dielectric microscopic particles; also, the feasibility of transferring OAM from laser beams into cold atoms [8, 9], and even store it in matter ensembles [0], has been discussed and shown. Furthermore, a couple of dynamical variables of light, involving compositions of linear and angular momenta, have been theoretically identified related to propagation invariant vector fields with transversal elliptic [] and parabolic [2] symmetries. In the paraxial regime, the elliptic [3] and parabolic [4] beams are well known; they have been generated by diverse means [5 7] and the orbital angular momentum per photon of the elliptic beams discussed [8]. Proposed applications for these light fields carrying compositions of linear and angular momenta include the optical manipulation of atoms, e.g. thermal atomic cloud splitters [9] and atomic billiards [20], as well as the optical manipulation of dielectric microscopic particles in elliptic, hyperbolic and parabolic trajectories in optical tweezers. The transition from mere optical manipulation towards optical control of matter, and further uses in other fields, e.g. in quantum information protocols [2, 22], implies measuring the linear and orbital angular momenta, or compositions of these, of arbitrary light fields. Here, a simple schema to calculate momenta for scalar wave fields through the Fourier transform, based on the Hilbert space model for separation of variables of Helmholtz equation [23, 24], is presented for the benefit of advanced undergraduate and graduate students interested in the OAM of light and related topics. In the following section, a brief overview of the Helmholtz equation, its symmetries and separable solutions, is given; the separable solutions for Cartesian, circular-cylindrical and elliptic-cylindrical coordinate systems are shown as an example. In the third section, wave decompositions and integral transforms related to the Helmholtz equation are treated; Plancherel theorem, the keystone of the proposed schema, is presented in this section. Fourier and Circular-cylindrical wave transforms of the separable solutions presented in section 2 are calculated for the sake of clarity. Next, the actual schema to calculate 2

3 momenta through Fourier transforms is discussed and the mean value of the orbital angular momentum for Plane and Elliptical waves is calculated. Finally, a summary and plausible applications are discussed. II. SYMMETRIES OF THE HEMLHOLTZ EQUATION The Helmholtz differential equation, ( 2 +k 2) Ψ = 0, () where the symbol 2 is the Laplacian operator and k is the magnitude of the wave vector k = k(sinϑcosϕ,sinϑsinϕ,cosϑ), and ϑ and ϕ are the Euler angles, is separable in eleven coordinate systems [25 27]. Each of these coordinate system defines a three element set of operators that span a Lie algebra [23, 24]. These operators fulfil the eigenvalue equation O (i) j Ψ = o (i) j Ψ, (2) with i =,...,, and j =,2,3. By defining a characteristic label set given by the three { } eigenvalues corresponding to a given coordinate system, Λ i =, it is possible to label the set of functions that solve the reduced wave equation in the i-th coordinate system as {Ψ Λi }. A common element to all characteristic label sets is the square linear momentum that fulfils p 2 Ψ Λi 2 Ψ Λi = k 2 Ψ Λi ; i.e., the first element of each and every characteristic label set Λ j is the magnitude of the wave vector, o (i) = k. The four coordinate systems with axial symmetry, Cartesian (i = ), circular- (i = 2), elliptic- (i = 3), parabolic- (i = 4) cylindrical, share the z-component of the linear momentum as the second dynamical variable, p z Ψ Λi = (kcosθ)ψ Λi ; i.e., for these four coordinate systems the second element of each of the characteristic label set is the Euler angle θ, o (i) 2 = θ with i =,...,4. There are five coordinate systems that share the orbital angular momentum, l z Ψ Λi = nψ Λi, with l = r p; for these, the second element of each of the characteristic label set is the discrete orbital angular momentum number n, also known as topological charge. Finally, there are another thirteen operators that can be easily found in the literature [23, 24]. In optics, there is a particular interest in propagation invariant scalar wave fields. There exist four fundamental separable families with this characteristic; for the sake of brevity, we will present just the cases for plane, circular- and elliptic-cylindrical waves. 3 o (i) j

4 A. Plane waves In the Cartesian coordinate system, x = (x,y,z), the separable solution of the Helmholtz differential equation are plane waves, Ψ Λ = (sinϑ) /2 e ık(xsinϑcosϕ+ysinϑsinϕ+zcosϑ), π ϕ < π, 0 ϑ < π, (3) The conserved operators are the linear momenta p y and p z, p y Ψ Λ = (ksinϑsinϕ)ψ Λ, (4) p z Ψ Λ = (kcosϑ)ψ Λ, (5) yielding a characteristic label set given by the magnitude of the wave vector and the two Euler angles, Λ = {k,ϑ,ϕ}. B. Circular-cylindrical waves For circular-cylindrical coordinates, x = (rcosφ,rsinφ,z), the separable solution of the Helmholtz differential equation are circular-cylinder waves, Ψ Λ2 = ı n (2πsinϑ) /2 J n (ksinϑ r)e ı(nφ+kzcosϑ), 0 ϑ < π, n = 0,±,±2,... (6) where the function J n ( ) is the Bessel function of the first kind. These scalar wave fields, so called Bessel waves, conserve the linear momenta p z and the orbital angular momenta l z, p z Ψ Λ2 = (kcosϑ)ψ Λ2, (7) l z Ψ Λ2 = nψ Λ2, (8) thusthecharacteristic labelsetisgivenbythemagnitudeofthewave vector, oneoftheeuler angles and the orbital angular momentum number, or topological charge, Λ 2 = {k,ϑ,n}, in that order. C. Elliptic-cylindrical waves In elliptic-cylindrical coordinates (x = f coshξcosη,y = f sinhξsinη, z), the separable solution of the Helmholtz wave equation are the so called even and odd Mathieu waves, in 4

5 that order, Ψc Λ3 = (sinϑ) /2 c n Ce(ξ,q)ce(η,q)e ıkzcosϑ, 0 ϑ < π, n = 0,2,... (9) Ψs Λ3 = (sinϑ) /2 s n Se(ξ,q)se(η,q)e ıkzcosϑ, 0 ϑ < π, n =,2,... (0) where the normalization coefficients are given by the expressions c 2n = c 2n+ = s 2n+ = s 2n+2 = A (2n) 0 ce 2n (0;q)ce 2n (π/2;q), () (2n+) qa (2n+) qb qb (2n+2) 2 ce 2n+ (0;q) v ce 2n+ (v;q) v=π/2, (2) v se 2n+ (v;q) v=0 se 2n+ (π/2;q), (3) v se 2n+2 (v;q) v=0 v se 2n+2 (v;q) v=π/2, (4) where ce n (se n ) and Ce n (Se n ) are the ordinary and extraordinary n-th order Mathieu functions of even (odd) parity, sometimes named angular and radial Mathieu functions, in that order. The separation constant is given by q = (fksinϑ/2) 2 with 2f being the inter-foci distance that defines a particular elliptic cylindrical coordinate system. Mathieu waves conserve, apart from the z-component of the linear momenta, a linear superposition of orbital angular momenta l z and linear angular momenta p x (l 2 z +f 2 p 2 x)ψc Λ3 = a n Ψc (3) {k,γ,n}, (5) (l 2 z +f2 p 2 x )Ψs Λ 3 = b n Ψs (3) {k,γ,n}, (6) such that the characteristic label set is given by the magnitude of the wave vector, one Euler angles related to the z component of the wave vector and the discrete order number, Λ 3 = {k,ϑ,n} in that order. The even (odd) eigenvalue a n (b n ) is fixed for a given combination of n and q. III. WAVE DECOMPOSITIONS AND INTEGRAL TRANSFORMS Each one of the eleven sets of separable functions, {Ψ Λi }, that solve the reduced wave equation forms an orthonormal basis, dx Ψ Λi (x)ψ Λ i (x) Ψ Λi Ψ Λ i = δ(λ i Λ i), (7) 5

6 where the bracket,,, notation is used as a shorthand for integration over configuration space, the upper bar notation,, stands for complex conjugation andthe notation δ( ) stands for a combination of Dirac and Kronecker deltas particular to each coordinate system. A Ψ Λi -wave decomposition of an arbitrary well behaved function Φ(x) can be defined for each of the eleven orthonormal basis, Φ(x) = dλ i φ(λ i )Ψ Λi (x), (8) alongside an integral transform, φ(λ i ) = dx Ψ Λi (x)φ(x) Ψ Λi Φ. (9) It is possible to show that each and every one of these eleven wave decompositions can be related to each other through their amplitudes, by substituting the corresponding wave decomposition, Eq.(8), of Ψ Λi (x) and Φ(x) into the integral transform given by Eq.(9), φ(λ i ) = dλ j ψλi (Λ j )φ(λ j ), (20) In shorthand, Ψ Λi Φ = ψ Λi φ Λj, (2) where double angle-brackets, Λj, have been used to denote integration over the phase space defined by the characteristic label set Λ j. The above expression, Eq.(2), is the well known Plancherel theorem[28 30]; it is particularly useful when both the integral transforms related to a Ψ Λj -wave decomposition of Φ and Ψ Λi are known. Unitarity of the eleven Helmholtz-related integral transforms follows straightforward from Eq.(2) yielding Parceval theorem, Φ Φ = φ φ Λi. (22) A. Plane wave decomposition and Fourier transform Take as an example the plane wave decomposition Φ(x) = dλ φ(λ )(sinϑ) /2 e ık x, (23) 6

7 related to the Fourier transform φ(k) = dx (sinϑ) /2 e ık x Φ(x). (24) The Fourier transform of a plane wave is given by the expression ψ Λ (Λ ) = (sinγ) /2 δ(ϕ ϕ )δ(ϑ ϑ ), (25) where the function δ( ) is the Dirac delta function. A circular-cylindrical wave has a Fourier transform ψ Λ 2 (Λ ) = (2πsinϑ ) /2 e ınϕ δ(θ ϑ ). (26) For elliptic-cylindrical waves, Fourier transforms are ψc Λ 3 (Λ ) = (πsinϑ) /2 ce n (ϕ,q)δ(ϑ ϑ ), (27) ψs Λ 3 (Λ ) = (πsinϑ) /2 se n (ϕ,q)δ(ϑ ϑ ). (28) B. Circular-cylindrical wave decomposition and related integral transform Now, take the circular cylindrical wave decomposition, also known as Bessel wave decomposition, Φ(x) = n π 0 sinϑdϑ φ(λ 2 )ı n (2πsinϑ) /2 J n (ksinϑ ρ)e ı(nϕ+kzcosϑ). (29) This circular wave decomposition defines a Bessel transform related to the Helmholtz equation; let us call it a Helmholtz-Bessel transform (not to be confused with the Fourier-Bessel, so called Hankel, transform) for the sake of having a short name for it, φ(λ 2 ) = Ψ Λ2 Φ = dx Φ(x)( ı) n (2πsinϑ) /2 J n (ksinϑ ρ)e ı(nϕ+kzcosϑ). (30) Sometimes, it can be bothersome to calculate the value of this expression. If the Fourier transformφ(λ ) of the function Φ(x) isknown, it may be simpler to calculate the coefficients of the Bessel wave decomposition, by using Plancherel theorem, Eq.(2), as π φ(λ 2 ) = (2π) /2 (sinϑ) /2 dϕ φ(k,ϑ,ϕ)e ınϕ. (3) 7 π

8 For example, the Helmholtz-Bessel transform of a plane wave, described by the label set Λ, is given by the expression, ψ Λ (Λ 2 ) = (2πsinϑ) /2 e ınϕ δ(ϑ ϑ ), (32) while that of a Bessel wave is given by, ψ Λ 2 (Λ 2 ) = (sinϑ) /2 δ(ϑ ϑ )δ n,n, (33) where the symbol δ, stands for Kronecker delta. And that of a Mathieu wave is ψc Λ 3 (Λ 2 ) = 2 /2 A (n ) 2m+(n mod2) (q )δ(ϑ ϑ )δ n,2m+(n mod2), (34) ψs Λ 3 (Λ 2 ) = 2 /2 B (n ) 2m+(n mod2) (q )δ(ϑ ϑ )δ n,2m+(n mod2), (35) with m = 0,,2,... and the coefficients A (i) j ( B (i) j ) are the j-th coefficients of the cosine (sine) expansion of the ordinary Mathieu functions of i-th order. IV. QUANTIFYING MOMENTA THROUGH THE FOURIER TRANSFORM It is possible to calculate the mean value of one of the aforementioned sixteen conserved operators O (i) j, as well as linear superpositions of them, for any given wave Ψ using its Ψ Λi -wave decomposition, Φ O (i) j Φ = Λ i o (i) j Ψ Λ i Φ 2. (36) Via Plancherel theorem, Eq.(2), it is possible to give this expression in terms of Fourier transforms, which may be already available and provide an easier path to calculate the overlaps, Φ O (i) j Φ = Λ i o (i) j ψ Λ i φ Λ 2, (37) where the phase space defined by the label set Λ = {k,ϑ,ϕ} yields ψ φ Λ = dk π dϕ π 0 π 0 sinϑdϑ ψ Λi (Λ )φ(λ ). (38) 8

9 A. Orbital angular momentum of plane waves It is simple to calculate the well-known null mean value of the orbital angular momentum l z for a plane wave, Λ = {k,ϑ,ϕ}, by using Eq.(37) and Eq.(3), Ψ Λ l z Ψ Λ = ñ ψ ψ Λ2 Λ Λ 2, (39) = ñ= ñ= ñ (sinϑ) /2 (2πsin ϑ) /2 e ınϑ δ( ϑ ϑ) 2 (40) = 0. (4) B. Orbital angular momentum of elliptic-cylindrical waves Again, by use of Eq.(37) and Eq.(34), it is possible to calculate the mean value of the angular momenta for any given elliptic-cylindrical wave as Ψc Λ3 l z Ψc Λ3 = = = Ψs Λ3 l z Ψs Λ3 = = = ñ= ñ ψ Λ2 ψc Λ3 Λ 2, (42) [m+(nmod2)/2] m=0 m=0 ñ= ( sinϑ/sin ϑ ) /2A (n) 2m+(nmod2) δ( ϑ ϑ) 2, [m+(nmod2)/2] A (n) 2m+(nmod2) δ(0) 2, (43) ñ ψ Λ2 ψs Λ3 Λ 2, (44) [m+(nmod2)/2] m=0 m=0 ( sinϑ/sin ϑ ) /2B (n) 2m+(nmod2) δ( ϑ ϑ) 2. [m+(nmod2)/2] B (n) 2m+(nmod2) δ(0) 2. (45) Note that the divergences in all these expressions are related to the fact that ideal scalar wave solutions to the differential Helmholtz equation have a divergent normalization, Eq.(7). V. FINAL REMARKS A schema to calculate momenta related to the symmetries of the Helmholtz differential equation has been presented. The procedure relies on the orthonormal bases and the sixteen 9

10 conserved operators derived from the eleven coordinate systems where the reduced wave equation accepts a separable solution. It has been shown that all the integral transforms, or wave overlaps, yielding the coefficients for the wave decompositions related to the aforementioned orthonormal basis can be calculated through Fourier transforms via Plancherel theorem. These overlaps are particularly useful to calculate the mean values of the conserved operators, and linear superposition of them, constructed from the three linear and three angular momenta. As an example, the mean value of the orbital angular momentum of Elliptic-cylindrical waves, so called Mathieu waves, has been presented. With further efforts the schema can be extended from these ideal scalar wave fields to real world beams; by using the orthonormal bases given by Gauss, Bessel-, Mathieu- [3] and Weber-Gauss [32] beams it may be possible to implement computational procedures to calculate momenta of captured experimental images. An extension of these procedures to non-paraxial optical vector fields is plausible. Moreover, integral transforms and decompositions are of use for those in the quantum optics fields as they define decompositions between boson modes; e.g. expressing creation/annihilation operators of Bessel modes in terms of creation/annihilation operators of Plane wave modes [33]. Furthermore, to give just an example from outside the field of optics, in quantum mechanics, the treatment of molecular torsional levels for water and similar molecules is related to Ince polynomials that are a solution to the Helmholtz differential equation in paraboloidal coordinates [34]. In these cases, the presented schema can help in calculating suitable combinations of linear and orbital angular momenta of the different torsional eigenstates of the molecules. Acknowledgments The author is grateful to Prof. Dr. Eugenio Ley Koo who posed the interesting question on how to calculate momenta of arbitrary wave functions through wave decompositions and gave valuable comments and suggestions through the writing stage of this manuscript. 0

11 References [] L. Allen, S. M. Barnett, and M. J. Padgett, Optical Angular Momentum (CRC Press, 2003). [2] D. L. Andrews, Structured Light and Its Applications: An Introduction to Phase-Structured Beams and Nanoscale Optical Forces (Elsevier, 2008). [3] R. A. Beth, Phys. Rev. 50, 5 (936). [4] L. Allen, M. W. Beijersbergen, R. J. C. Spreeuw, and J. P. Woerdman, Phys. Rev. A 45, 885 (992). [5] H. He, M. E. J. Friese, N. R. Heckenberg, and H. Rubinsztein-Dunlop, Phys. Rev. Lett. 75, 826 (995). [6] G. Milne, K. Dholakia, D. McGloin, K. Volke-Sepulveda, and P. Zemánek, Opt. Express 5, 3972 (2007). [7] X.-L. Wang, J. Chen, Y. Li, J. Ding, C.-S. Guo, and H.-T. Wang, Phys. Rev. Lett. 05, (200). [8] J. W. R. Tabosa and D. V. Petrov, Phys. Rev. Lett. 83, 4967 (999). [9] A. Picón, A. Benseny, J. Mompart, J. R. V. de Aldana, L. Plaja, G. F. Calvo, and L. Roso, New J. Phys. 2, (200). [0] D. Moretti, D. Felinto, and J. W. R. Tabosa, Phys. Rev. A 79, (2009). [] B. M. Rodríguez-Lara and R. Jáuregui, Phys. Rev. A 78, (2008). [2] B. M. Rodríguez-Lara and R. Jáuregui, Phys. Rev. A 79, (2009). [3] J. C. Gutiérrez-Vega, M. D. Iturbe-Castillo, and S. Chávez-Cerda, Opt. Lett. 5, 493 (2000). [4] M. A. Bandres, J. C. Gutiérrez-Vega, and S. Chávez-Cerda, Opt. Lett. 29, 44 (2004). [5] S. Chávez-Cerda, M. J. Padgett, I. Allison, G. H. C. New, J. C. Gutiérrez-Vega, A. T. O Neil, I. MacVicar, and J. Courtial, J. Opt. B: Quantum Semiclass. Opt. 4, S52 (2002). [6] M. B. Alvarez-Elizondo, R. Rodríguez-Masegosa, and J. C. Gutiérrez-Vega, Opt. Express 6, 8770 (2008). [7] R. J. Hernández-Hernández, R. A. Terborg, I. Ricardez-Vargas, and K. Volke-Sepúlveda, Appl. Opt. 49, 6903 (200). [8] C. López-Mariscal, M. A. Bandres, J. C. Gutiérrez-Vega, and S. Chávez-Cerda, Opt. Express

12 3, 2364 (2005). [9] B. M. Rodríguez-Lara and R. Jáuregui, Phys. Rev. A 80, 083R (2009). [20] R. pérez-pascual, B. M. Rodríguez-Lara and R. Jáuregui, arxiv: v [quant-ph] (200). [2] J. Leach, B. Jack, J. Romero, A. K. Jha, A. M. Yao, S. Franke-Arnold, D. G. Ireland, R. W. Boyd, S. M. Barnett, and M. J. Padgett, Science 329, 662 (200). [22] A. K. Jha, J. Leach, B. Jack, S.Franke-Arnold, S.M. Barnett, R. W. Boyd, andm. J.Padgett, Phys. Rev. Lett. 04, 0050 (200). [23] C. P. Boyer, E. G. Kalnins, and W. Miller, Nagoya Math. J. 60, 35 (976). [24] W. Miller, Symmetry and Separation of Variables, Encyclopedia of Mathematics and Its Applications (Cambridge University Press, 984). [25] L. P. Eisenhart, Phys. Rev. 45, 427 (934). [26] P. M. Morse and H. Feshbach, Methods of theoretical physics, vol. of International series in pure and applied physics (McGraw-Hill, 953). [27] F. M. Arscott and A. Darai, IMA J. Appl. Math. 27, 33 (98). [28] M. Plancherel, Rend. Circ. Mat. Palermo 30, 289 (90). [29] E. C. Titchmarsh, Proc. London Math. Soc. 23, 279 (925). [30] G. N. Watson, Proc. London Math. Soc. 35, 56 (933). [3] J. C. Gutiérrez-Vega and M. A. Bandres, J. Opt. Soc. Am. A 24, 25 (2007). [32] B. M. Rodríguez-Lara, J. Opt. Soc. Am. A 27, 327 (200). [33] I. Bialynicki-Birula and Z. Bialynicka-Birula, Opt. Commun. 264, 342 (2006). [34] L. F. Roncaratti and V. Aquilanti, Int. J. Quantum Chem. 0, 76 (200). 2