OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 14, Iss. 7 — Apr. 3, 2006
  • pp: 2888–2897
« Show journal navigation

Optical vortices with large orbital momentum: generation and interference

Anatoliy A. Savchenkov, Andrey B. Matsko, Ivan Grudinin, Ekaterina A. Savchenkova, Dmitry Strekalov, and Lute Maleki  »View Author Affiliations


Optics Express, Vol. 14, Issue 7, pp. 2888-2897 (2006)
http://dx.doi.org/10.1364/OE.14.002888


View Full Text Article

Acrobat PDF (405 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We demonstrate a method for generation of beams of light with large angular momenta. The method utilizes whispering gallery mode resonators that transform a plane electromagnetic wave into high order Bessel beams. Interference pattern among the beams as well as shadow pictures induced by the beams are observed and studied.

© 2006 Optical Society of America

1. Introduction

Optical vortices or light beams with nonzero angular momenta are interesting not only because of their underlying physics, but also because of their potential applications. Generation of such vortices and the study of their properties has recently attracted considerable attention. The beams have been realized with both passive and active schemes (see [1

1. L. Allen, M. J. Padgett, and M. Babiker, “The orbital angular momentum of light”, Progress in Optics 39, 291–372 (1999). [CrossRef]

, 2

2. L. Allen, “Introduction to the atoms and angular momentum of light special issue,” J. Opt. B 4, S1–S6 (2002). [CrossRef]

, 3

3. E. Santamato, “Photon orbital angular momentum: problems and perspectives,” Progress in Physics 52, 1141–1153 (2004).

, 4

4. D. McGloin and K. Dholakia, “Bessel beams: diffraction in a new light,” Contemp. Phys. 46, 15–28 (2005). [CrossRef]

] for review). Nonetheless, generation of beams with orbital angular momentum exceeding 104ħ per photon, though in principle possible, [5

5. J. Courtial, K. Dholakia, L. Allen, and M. J. Padgett, “Gaussian beams with very high orbital angular momentum,” Opt. Commun. 144, 210–213 (1997). [CrossRef]

], is not easily achievable experimentally. To our knowledge, the highest orbital angular momenta of photons observed thus far are approximately 25 ħ [5

5. J. Courtial, K. Dholakia, L. Allen, and M. J. Padgett, “Gaussian beams with very high orbital angular momentum,” Opt. Commun. 144, 210–213 (1997). [CrossRef]

] and 30ħ [6

6. J. E. Curtis, B. A. Koss, and D. G. Grier, “Dynamic holographic optical tweezers,” Opt. Commun. 207, 169–175 (2002). [CrossRef]

]. Beams with angular momenta up to 300ħ were recently demonstrated [7

7. S. Sundbeck, I. Gruzberg, and D. G. Grier, “Structure and scaling of helical modes of light,” Opt. Lett. 30, 477–479 (2005). [CrossRef] [PubMed]

]. In this paper we discuss another method for the generation of photons with large orbital momentum.

Fig. 1. Scheme of the experiment.

The feasibility of the applications of the optical waves possessing angular momenta directly depends on their propagation distance in the free space. A high order running Bessel wave that penetrates into a cylindrical waveguide from the resonator through evanescent field coupling cannot leave the waveguide. It runs until the end of the waveguide and reflects, because of the total internal reflection. It is possible to manipulate the wave inside the waveguide. However, the complete confinement of the waves in the waveguide reduces the spectrum of their applications. We used a tapered waveguide to release the Bessel beam into free space.

We have been able to generate Bessel beams with orders exceeding a thousand, and have studied interference of multiple beams with different angular momenta as well as pictures of shadow of an object placed in front of them. The shadow has properties that arise from the far field interference of the beams of different orders at different points at the exit of the waveguide. This interference experiment allowed us to estimate the propagation distance of the optical vortices. The distance drastically decreases with an increase of the value of the angular momentum. We explain the results theoretically and discuss the possibility of increasing the propagation distance. The experiments with shadows have shown that a shadow of a straight rod illuminated with a Bessel beam is not a straight line. This result is explained theoretically in the frame of geometrical, not wave, optics.

2. Experiment

We fabricated several low contrast fused silica WGM resonators being integral parts of tapered waveguides (Fig. 1). In one experiment, a segment of a fused silica multimode fiber rod is cleaved and the cleavage plane is polished. The diameter of the fiber is 1 mm. The other end of the rod is stretched into a conical shape in a hydrogen burner. This end is fused to form a WGM resonator having about 300 μm in diameter as an integral part of a cone. This cone starts at the WGM at diameter of about 250-280 μm and expands to 1 mm in diameter over 1 cm distance (Fig. 2). Another resonator we built has a 500 μm diameter, and the taper diameter changes from 450 μm to 3 mm over a 3 cm distance.

Fig. 2. A low contrast whispering gallery mode resonator for a Bessel beam generation. The resonator is created by the local increase of the waveguide radius. A tapered waveguide is used to generate a Bessel beam in the free space. Increasing the waveguide radius changes the wave vector of the light, guaranteeing escape of the light from the waveguide into free space. (A) A low contrast fused silica WGM resonator attached to the tapered waveguide. (B) Coupling light into the WGMs of the resonator using cleaved fiber. (C) Tapered fiber used to release generated Bessel beam into free space.
Fig. 3. (A) Far field interference pattern of multiple copropagating Bessel beams. (B) Far field shadow of a 250 μm-thick needle placed to the beam emerging the fiber taper.

The resonators were optically pumped with 650 nm light using an angle polished fiber coupler. We excited a family of WGMs that interacted with the modes of the tapered fiber. The fiber coupler is shown in (Fig. 2(B)). It is easy to see the glow on the resonator surface resulting from the scattering of the light in the WGMs, which abruptly disappears in the vicinity of the coupling point between the resonator and the tapered fiber. This glow shows the geometrical localization of light in the resonator.

As the WGMs of the resonator decay into the fiber, light propagates and exits the fiber into the free space as truncated Bessel beams. Because the ratio of the taper entrance over the exit diameters was small (rexit /rentrance = 6 or less) and the mode order was high (> 103), the propagation distance of Bessel beams in free space did not exceed ten millimeters. The beams spread out creating peculiar interference shapes in the far field region (Fig. 3(A)). This observation has a certain similarity with the interference pattern of zero and first order Laguerre-Gauss beams presented in [9

9. E. J. Galvez, P. R. Crawford, H. I. Sztul, M. J. Pysher, P. J. Haglin, and R. E. Williams, “Geometric phase associated with mode transformations of optical beams bearing orbital angular momentum,” Phys. Rev. Lett. 90, 203901 (2003). [CrossRef] [PubMed]

]. In what follows it is shown that in our case the interference occurs between multiple high-order truncated Bessel beams.

Aside from the interference we have observed peculiar shapes of the shadow of a straight, thin, opaque object parallel to the waveguide surface that crosses the symmetry axis of the beam possessing large angular momentum (Fig. 3(B)) (also see Fig. 9 for more details of the experiment). To generate the shadow we used a 250 μm needle as the object. The shadow is orthogonal to the object in the center of the interference pattern, and is parallel, but displaced, far from the center. The shape of the shadow depends on the distance of the object from the surface of the taper. The closer the wire is, the longer is the region of the orthogonal shadow. We have repeated the same experiment using 532 nm laser and a 25μm thick piece of a tungsten wire, and obtained similar result. The dynamics of the behavior of the shadow with the distance change is demonstrated in Fig. 4. In what follows we explain the result of the experiment theoretically.

Fig. 4. (0.7 Mb) Movie. A shadow of a straight wire illuminated with the light possessing large angular momentum. The distance between the wire and the waveguide surface changes from 12.5 mm to zero. Two mixed interference patterns come from two polarizations of the light.

3. Theory

3.1. Interference

The interference pattern (Fig. 3) could be explained if we consider the interference of multiple radially truncated high order Bessel waves with nearly the same diameter but different angular momenta. Light from an angle polished fiber coupler (Fig. 2(B)) excites many modes of the resonator with different quantum numbers and the same frequency. Those WGMs excite the Bessel waves in the taper. The waves have different propagation constants along the Z axis and, hence, different dispersion. The initially narrow angular distribution of light in the taper broadens up with the propagation distance.

Fig. 5. To the explanation of the influence of the dispersion of the Bessel waves on az-imuthal localization of light. Let us assume that the fiber coupler excites two Bessel waves shown with red and green curves in a cylindrical waveguide. The initially localized waves become azimuthally delocalized during the propagation. They cross the plane with coordinate Z at different points described by different azimuthal angles.

Let us consider a fiber coupler that excites Bessel waves with various angular moments in a cylindrical waveguide (Fig. 5). For the sake of simplicity we assume that the waveguide radius does not change and excitation of the modes is possible without a WGM resonator. The light leaving the fiber coupler is azimuthally localized. The localization disappears with distance because of the spatial dispersion of the excited modes. Results of an exact numerical simulation of the thought experiment (Fig. 5) are shown in (Fig. 6).

The longitudinal propagation constant kz for a Bessel wave can be found from equation k02 = kz2 + kl2 as a function of distance z and azimuthal index l:

kz(z)k02(lnr(z))2,
(1)

where k 0 = ω 0/c is the total wave vector,ω 0 is the frequency, n is the index of refraction of the material, kl is the transverse propagation constant, and r(z) is the radius.

The complex amplitude of a Bessel wave with azimuthal number l can be expressed as

f(r,z,ϕ,l)=A0lrexp[j(ϕ+ϕ0(l))l+j0zkz1z1ldz1]
(2)

at the distance z from the beginning of the fiber taper (we assume that the radius changes adiabatically so that beams with different azimuthal numbers do not interact), where ϕ 0 (l) and A 0(l,r) are the phase and amplitude of the wave at z = 0.

We assume that a group of Bessel beams with the azimuthal index distributed from l 0 to l 0 - Δl and the same frequency a>0 is created at the beginning of the fiber. The summation of all the waves inside the dispersion interval Δl gives:

F(r,z,ϕ,l0,Δl)=l=l0Δll0f(r,z,ϕ,l)=l0Δll=l0A0lrexp(j(ϕ+ϕ0(l))l+j0zk02(lnr(z1))2dz1).
(3)

F(r,z,ϕ, l 0l) is the complex field distribution in the taper.

Fig. 6. Azimuthal angle distribution of the real part of the optical field at different distances from the fiber coupler. The angular width of the distribution is increasing monotonically with the distance because of the different propagation constants of the Bessel waves supported by the waveguide and having different spacial distributions and the same frequency.

To obtain Fig. 6 we have considered a tapered optical fiber with refractive index n= 1.45, as well as initial and final radii r 0 =0.3 mm and rn = 1 mm, respectively. We assumed that all the waves initially have the same phase and amplitude (ϕ 0(l) ≡ π and A 0(l,r) ≡ 1). Fig. 6 picture shows angular dependence of [(F(r,z,ϕ(mod(2π)),l 0l) + c.c.)/2] for the given parameters and l = 3000, Δl = 1000.

The helix-shaped far field interference pattern, easily seen in Fig. 3(A), is created by the sum of all the waves at the exit of the taper. The length of the helical curve (the number of loops in the interference pattern) directly depends on the phase distribution shown in Fig. 6. For instance, if the taper is short enough, one will see the interference pattern as an unclosed helix. This is the case for a taper shorter than zn = 3 mm with the other parameters as used in the simulations resulting in Fig. 6. The pattern will have at least one loop for a longer taper (see, e.g., Fig. 7).

The distribution of the far field interference pattern is determined by the Fresnel-Kirchhoff diffraction integral

Bxy=CSFx'y'eik(xx'+yy')Hdx'dy',
(4)

where F(x,y) is given by Eq. (3) at the end of the taper crossection surface (here z, l and Δl are constants), C is a constant, S is the aperture area, and H is the distance between the screen and the aperture. A direct evaluation of the integral is complicated, but we can approximate this integral by assuming that amplitudes of the particular Bessel waves (A 0(l,r)) are nonzero only in a thin belt close to the circumference of the taper crossection. The radial thickness of the belt is approximately Δr = rn /l 2/3, as for high order WGMs. Splitting the circumference of the taper surface at small segments 1 ≫ Δϕλ/(2πΔr), we estimate the integral as

Bϕθ=C10ΔϕF(r=rnΔr2,z=znϕ+ϕ1,l0,Δl)ejk0rϕ1sinθk0rdϕ1,
(5)

where r = rn - Δr/2, C1 is a normalization parameter. To find the optimum step Δϕopt we first approximate the angle θ for the region of the localization of the interference pattern and take Δϕoptλ/[2π(sinθr]. For sinθ ≈ 0.22 and Δr ≈ 4.3 × 10-4 cm, which are the values for our taper, we get Δϕopt ≃ 0.11. The result of the optimum step calculation is shown in Fig. 7.

Fig. 7. Numerical simulation of the interference pattern obtained for the Bessel beams with l = 3000, Δl = 1000, wavelength λ = 650 nm, and zn = 20 mm (see text for the explanation).

3.2. Propagation

Let us roughly estimate the propagation distance LB for the generated Bessel beams. Because the beams can be considered as a superposition of plane waves propagating at the angle θ ≈ arcsin(kl (rn )/k 0) to the surface of the crossection of the taper, the propagation distance corresponds to LB = rn /arctg(θ), where kl (rn ) = kl (r 0)r 0/rn and at the very beginning of the taper kl (r 0) = k 0. Now the equations for the beam divergence θ and beam propagation length LB can be written as θ ≈ arcsin(r 0/rn ) and LBrn2/r 0, θ << 1. For parameters of our experiment r 0 = 0.5 mm, rn = 3 mm, and l = 7000 we get LB ≈ 18 mm. At this distance the Bessel beam decays significantly, and its phase structure is very disturbed, though its momentum is still preserved. Experimental data demonstrate a Bessel beam propagation length of the same order. In order to get a distance of a meter, the exit radius of the taper should be at least seven times larger. It is difficult to reproduce such a taper experimentally using our technique of fabrication.

Let us study the divergence of a single truncated Bessel beam numerically to demonstrate the nature of the beam decay. For the sake of simplicity of the computation we consider comparably low order Bessel beams here. We take a tapered fiber with ending radius rn = 10 μm as a source of light. According to Huygens the complex field amplitude can be described by the diffraction integral [11

11. L. Landau and E. Lifshitz, Classical Theory of Fields (Reed International Educational and Professional Publishing, Oxford, 1980).

]

B(x,y,z)=C2Sfejk0RRds,
(6)

where amplitude of optical field f = Jm (klr)e -jlϕ, Jm (klr) is the Bessel function of m-th order, C 2 is a scale parameter, R is the distance between point in the plane of end of the fiber and the point (x,y, z) where diffraction pattern is obtained.

In our case it is useful to employ cylindrical coordinates. The azimuth angle ϕ does not contain information on the divergence because our system is symmetrical. We can exclude it from the final consideration without any mathematical reduction. Now the diffraction integral looks like:

Brz=C20r002πJm(klr')Rejk0Rr''dr';R=(r'cosϕ'r)2+(r'sinϕ')2+z2
(7)

Please note that ϕ´ coordinate appeared under integral. After selection of certain angle where the optical field is to be calculated the function under the integral loose symmetry. This diffraction integral was computed numerically for m = 17 and l = 30 as well as for m = 10 and l = 15.

Fig. (8) represents results of simulations of the field intensity |B(r,z)|2 in the plane which is parallel to the fiber taper symmetry axis. The truncated Bessel beam does not change its shape but has radiative loss. This is seen as a straight Bessel beam and group of conically diverged beams in Fig. (8). Similar behavior was predicted for the zero-order Bessel-Gauss beam [10

10. F. Gori, G. Guattari, and C. Padovani, “Bessel-Gauss beams,” Opt. Commun. 64, 491–495 (1987). [CrossRef]

].

Fig. 8. Intensity pattern for the truncated Bessel beams leaving the fiber taper with 10 μm radius. The ordinate axis coincides with the symmetry axis of the taper. The pictures correspond to beams with numbers (A) m = 17 and l = 30, and (B) m = 10 and l = 15.

3.3. Shadow

We explain the shape of the shadow using principles of geometric optics. We assume that the Bessel beams are significantly truncated, such that ml - m and consider a one dimensional rod AA´ placed at the distance h from the fiber and H - h from the screen (Fig. 9). The surfaces of the fiber crossection and the screen are parallel to the rod. The point A of the rod is illuminated by the plane wave leaving point S on the fiber rim. Point B is the shadow of point A.

Fig. 9. Formation of the shadow (inset) of a straight rod AA´ illuminated with a Bessel beam.

Let us find angle γ (angle CFB), that determines the change of the position of the shadow (B) with respect the straight projection (C) of point A of the rod. To do this we need to know all the sides of triangle CFB.

Side CF is equal to the segment OD. The line segment SAB touches the surface of the fiber so angle OSD is equal to 90°. Segment AD is equal to h. Segment OS is equal to rn . Angle ADS is equal to 90°, and angle ASD is equal to α. Therefore

FC=rn2+(htgα)2.
(8)

Following the same logic we find

FB=rn2+(Htgα)2.
(9)

To find segment CB we note that angle ACB is equal to 90°, and angle CBA is equal to α, hence

CB=Hhtgα.
(10)

Using the theorem of cosines we get

cosγ=FC2+FB2CB22FCFBtgγ=rnHhtgαrn2+Hh(tgα)2.
(11)

Angle α can be directly found from the formula tgα = kz (rn )/kl (rn ). Finally, the distance from the shadow point (B) to the direct projection of the rod AA´ is equal to FBsiny.

Let us now assume that we generated a lot of Bessel beams with different kz (rn )/kl (rn ), as is done in the experiment. Each Bessel beam projects a point of the rod to different places of the screen. The complete shadow line is shown in the inset of Fig. (9). This result directly resembles our experimental observations.

4. Conclusion

Acknowledgments

The research described in this paper was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration and partial sponsorship from DARPA. A. Matsko appreciates useful discussions with V Ilchenko and I. Marinenko.

References and links

1.

L. Allen, M. J. Padgett, and M. Babiker, “The orbital angular momentum of light”, Progress in Optics 39, 291–372 (1999). [CrossRef]

2.

L. Allen, “Introduction to the atoms and angular momentum of light special issue,” J. Opt. B 4, S1–S6 (2002). [CrossRef]

3.

E. Santamato, “Photon orbital angular momentum: problems and perspectives,” Progress in Physics 52, 1141–1153 (2004).

4.

D. McGloin and K. Dholakia, “Bessel beams: diffraction in a new light,” Contemp. Phys. 46, 15–28 (2005). [CrossRef]

5.

J. Courtial, K. Dholakia, L. Allen, and M. J. Padgett, “Gaussian beams with very high orbital angular momentum,” Opt. Commun. 144, 210–213 (1997). [CrossRef]

6.

J. E. Curtis, B. A. Koss, and D. G. Grier, “Dynamic holographic optical tweezers,” Opt. Commun. 207, 169–175 (2002). [CrossRef]

7.

S. Sundbeck, I. Gruzberg, and D. G. Grier, “Structure and scaling of helical modes of light,” Opt. Lett. 30, 477–479 (2005). [CrossRef] [PubMed]

8.

A. B. Matsko, A. A. Savchenkov, D. Strekalov, and L. Maleki, “Whispering gallery resonators for studying orbital angular momentum of a photon,” Phys. Rev. Lett. 95, 143904 (2005). [CrossRef] [PubMed]

9.

E. J. Galvez, P. R. Crawford, H. I. Sztul, M. J. Pysher, P. J. Haglin, and R. E. Williams, “Geometric phase associated with mode transformations of optical beams bearing orbital angular momentum,” Phys. Rev. Lett. 90, 203901 (2003). [CrossRef] [PubMed]

10.

F. Gori, G. Guattari, and C. Padovani, “Bessel-Gauss beams,” Opt. Commun. 64, 491–495 (1987). [CrossRef]

11.

L. Landau and E. Lifshitz, Classical Theory of Fields (Reed International Educational and Professional Publishing, Oxford, 1980).

OCIS Codes
(140.3300) Lasers and laser optics : Laser beam shaping
(230.5750) Optical devices : Resonators
(260.1960) Physical optics : Diffraction theory

ToC Category:
Optical Devices

History
Original Manuscript: February 6, 2006
Revised Manuscript: March 28, 2006
Manuscript Accepted: March 29, 2006
Published: April 3, 2006

Citation
Anatoliy A. Savchenkov, Andrey B. Matsko, Ivan Grudinin, Ekaterina A. Savchenkova, Dmitry Strekalov, and Lute Maleki, "Optical vortices with large orbital momentum: generation and interference," Opt. Express 14, 2888-2897 (2006)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-14-7-2888


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. L. Allen, M. J. Padgett, and M. Babiker, "The orbital angular momentum of light," Progress in Optics 39, 291-372 (1999). [CrossRef]
  2. L. Allen, "Introduction to the atoms and angular momentum of light special issue," J. Opt. B 4, S1-S6 (2002). [CrossRef]
  3. E. Santamato, "Photon orbital angular momentum: problems and perspectives," Progress in Physics 52, 1141-1153 (2004).
  4. D. McGloin and K. Dholakia, "Bessel beams: diffraction in a new light," Contemp. Phys. 46, 15-28 (2005). [CrossRef]
  5. J. Courtial, K. Dholakia, L. Allen, and M. J. Padgett, "Gaussian beams with very high orbital angular momentum," Opt. Commun. 144, 210-213 (1997). [CrossRef]
  6. J. E. Curtis, B. A. Koss, and D. G. Grier, "Dynamic holographic optical tweezers," Opt. Commun. 207, 169-175 (2002). [CrossRef]
  7. S. Sundbeck, I. Gruzberg, and D. G. Grier, "Structure and scaling of helical modes of light," Opt. Lett. 30, 477-479 (2005). [CrossRef] [PubMed]
  8. A. B. Matsko, A. A. Savchenkov, D. Strekalov, and L. Maleki, "Whispering gallery resonators for studying orbital angular momentum of a photon," Phys. Rev. Lett. 95, 143904 (2005). [CrossRef] [PubMed]
  9. E. J. Galvez, P. R. Crawford, H. I. Sztul, M. J. Pysher, P. J. Haglin, and R. E. Williams, "Geometric phase associated with mode transformations of optical beams bearing orbital angular momentum," Phys. Rev. Lett. 90, 203901 (2003). [CrossRef] [PubMed]
  10. F. Gori, G. Guattari, and C. Padovani, "Bessel-Gauss beams," Opt. Commun. 64, 491-495 (1987). [CrossRef]
  11. L. Landau and E. Lifshitz, Classical Theory of Fields (Reed International Educational and Professional Publishing, Oxford, 1980).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Supplementary Material


» Media 1: MOV (708 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited