The Poynting vector and angular momentum of Airy beams

doi:10.1364/OE.16.009411

« Show journal navigation

The Poynting vector and angular momentum of Airy beams

H. I. Sztul and R. R. Alfano »View Author Affiliations

Optics Express, Vol. 16, Issue 13, pp. 9411-9416 (2008)
http://dx.doi.org/10.1364/OE.16.009411

View Full Text Article

Acrobat PDF (644 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Full Text
Article Info
References (15)
Cited By
Figures (4)
Metrics

Abstract

We analyze and describe the evolution of the Poynting vector and angular momentum of the Airy beam as it propagates through space. A numerical approach is used to show the Poynting vector follows the tangent line of the direction of propagation. A similar approach is used to show that while the total angular momentum of the Airy beam is zero, the angular momentum of the main intensity peak and the Airy “tail” are non-zero.

1. Introduction

Recently there has been a lot of attention paid to various interesting, non-trivial solutions to the paraxial wave equation. Some of these beams exhibit interesting features such as spin or orbital angular momentum[1

L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX, pp. 291–372 (1999). [CrossRef]

], diffraction-free propagation [2

J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]

], self reconstruction [3

J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651 (1987). [CrossRef]

], or acceleration. Depending on the coordinate system there are Hermite, Laguerre, Bessel, and Ince-Gaussian solutions to name just a few. Most recently there has been increasing attention paid to the Airy solution of the paraxial wave equation first described by Berry and Balazs in 1979 [4

M. V. Berry and N. L. Balazs, “Nonspreading wave packets,” Am. J. Phys. 47, 264 (1979). [CrossRef]

]. The demonstration of this solution was first given in the context of the Schrödinger equation describing a free particle. The 1-D and 2-D Airy solution was most recently shown experimentally by Siviloglou et. al. in optical beams of light [5

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]

, 15

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of Airy beams,” Opt. Lett. 33, 207–209 (2008). [CrossRef] [PubMed]

]. The Airy beam is said to be free of diffractive spreading and freely accelerating in the absence of any external potential [8

D. M. Greenberger, “Comment on ‘Non-Spreading Wave Packets’,” Am. J. Phys. 48, 256 (1980). [CrossRef]

In this paper the Poynting vector and angular momentum of the Airy beam as it propagates through space is investigated. First, classical electrodynamics will be used to numerically calculate the Poynting vector of the airy wave packet. Next, two methods are employed to analyze the angular momentum of the Airy beam. The first is calculating j⃗=⃗r×p⃗ while the second is calculating the angular momentum spectrum of this beam of light [9

L. Torner, J.P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13, 873 (2005). [CrossRef] [PubMed]

]. These methods provide complimentary information and insight as to how this beam travels through space.

2. Airy solution to the paraxial wave equation

The (2+1)D paraxial wave equation describes the wave propagation of the electric field ϕ as,

2 i \frac{\partial u (s_{x}, s_{y}, ξ)}{\partial ξ} + \nabla_{τ}^{2} u (s_{x}, s_{y}, ξ) = 0

(1)

where s_x =x/x ₀ and s_y =y/y ₀ are normalized transverse coordinates, ξ=z/k(x ₀+y ₀)² is a normalized propagation distance, x ₀ and y ₀ are normalization constants, k=2πn/λ₀ and ∇² _τ is the second partial derivative in the transverse direction.

The non-dispersive solution to this (2+1)D equation is [7

I. M. Besieris, A. M. Shaarawi, and R. W. Ziolkowski, “Nondispersive accelerating wave packets,” Am. J. Phys. 62, 519–521 (1994). [CrossRef]

u (s_{x}, s_{y}, ξ) = Ai (s_{x} - {(\frac{ξ}{2})}^{2}) Ai (s_{y} - {(\frac{ξ}{2})}^{2}) \exp [i (s_{x} \frac{ξ}{2} + s_{y} \frac{ξ}{2}) - i (\frac{ξ^{3}}{6})],

(2)

where Ai(x) is the Airy function and is shown in the inset at the center of Fig. 1. This type of beam can be produced experimentally with a specially designed diffraction grating [5

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]

In theory a non-dispersive beam of this sort would have infinite power. However in practice, this can not be the case because a beam can not propagate infinite power. As discussed in Refs. [5

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]

, 6

G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy Airy beams,” Opt. Lett. 32, 979 (2007). [CrossRef] [PubMed]

], to experimentally realize the Airy beam an initial condition must be employed to act as an exponential aperture function. Taking this initial condition as, u(s_x ,s_y ,0)=Ai(s_x )Ai(s_y )exp[a(s_x +s_y )], where a<<1 and is a positive parameter that limits the infinite energy in the Airy tail, the electric field amplitude becomes

u (s_{x}, s_{y}, ξ) = Ai (s_{x} - {(\frac{ξ}{2})}^{2} + i a ξ) Ai (s_{y} - {(\frac{ξ}{2})}^{2} + i a ξ) \times

\exp [a s_{x} + a s_{y} - (\frac{a ξ^{2}}{2}) - i (\frac{ξ^{3}}{6}) - i (\frac{a^{2} ξ}{2}) + i ξ \frac{(s_{x} + s_{y})}{2}] .

(3)

Figure 1(a-d) shows the intensity in the transverse plane at various ξ-positions when a=0.15 and x_o =y_o =0.1. All of the intensity peaks are changing position as ξ increases. The main intensity peak in the s_x -s_y plane of the Airy beam travels at 45°, along the line s_x =s_y , following the trajectories,

x_{m} = \frac{λ_{o}^{2} ξ^{2}}{16 π^{2} x_{o}^{2}} and y_{m} = \frac{λ_{o}^{2} ξ^{2}}{16 π^{2} y_{o}^{2}} .

(4)

Fig. 1. The intensity of the finite-energy Airy wave in the s_x -s_y plane, given by Eq. (3) with a=0.15 at ξ=0.025 (a), ξ=0.05 (b), ξ=0.075 (c), ξ=0.1 (d). The inset shows the behavior of Ai(x).

3. The Poynting vector of Airy beams

The rate of electromagnetic energy flow per unit area, or the Poynting vector, is a commonly known quantity in electrodynamics [10

J. D. Jackson, Classical Electrodynamics , (Wiley, New York, 1962).

, 11

M. Born and E. Wolf, Principles of Optics , 7th Edition, (Cambridge University Press, Cambridge, 1999).

]. This vector is routinely examined for plane waves but has received considerable attention in the literature with regard to Laguerre-Gaussian beams of light that have helical wavefronts [1

L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX, pp. 291–372 (1999). [CrossRef]

, 12

L. Allen and M. J. Padgett, “The Poynting vector in Laguerre-Gaussian beams and the interpretation of their angular momentum desnity,” Opt. Commun. 184, 67–71 (2000). [CrossRef]

, 13

H. I. Sztul and R. R. Alfano, “Double-slit interference with Laguerre-Gaussian beams,” Opt. Lett. 31, 999–1001 (2006). [CrossRef] [PubMed]

The Poynting vector is defined as [11

M. Born and E. Wolf, Principles of Optics , 7th Edition, (Cambridge University Press, Cambridge, 1999).

\vec{S} = (\frac{c}{4 π}) \vec{E} \times \vec{B},

(5)

where c is the speed of light. Given a vector potential A⃗= $\hat{ε}$ u(s_x ,s_y ,ξ)exp[ik ² ξ(x ₀+y ₀)], where $\hat{ε}$ is an arbitrary polarization and u(s_x ,s_y ,ξ) is the Airy field amplitude given by Eq. (3), we can use the E⃗ and B⃗-fields in the Lorenz gauge, as given by Ref. [1

L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX, pp. 291–372 (1999). [CrossRef]

], to calculate the time-averaged Poynting vector, <S⃗>. Assuming an x̂-polarized field, <S⃗> becomes [1

L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX, pp. 291–372 (1999). [CrossRef]

]

\frac{c}{4 π} < \vec{E} \times \vec{B} > = \frac{c}{8 π} (\vec{E} \times \vec{B^{*}} + \vec{E^{*}} \times \vec{B})

= \frac{c}{8 π} [i ω (u \nabla_{τ} u^{*} - u^{*} \nabla_{τ} u) + 2 ω k {∣ u ∣}^{2} \hat{ξ}] .

(6)

The $\hat{ξ}$ term in the above equation is the energy flow in the $\hat{ξ}$ -direction which is just proportional to the linear momentum density in that direction. This is typically the main contributing component of S⃗ in Gaussian optics. The first term is what we are really interested in here as it contributes a non-zero ŝ_x - and ŝ_y -component and an additional $\hat{ξ}$ -term to the Poynting vector.

Figure 2 (a-d) shows the numerically computed ŝ_x - and ŝ_y -components of the Poynting vector for a=0.15 at ξ=0.025, 0.050, 0.075, and 0.1, respectively. The direction and magnitude of the arrows (shown in red) correspond to the direction and magnitude of the energy flow in the transverse plane. The intensity of the Airy field is shown in the background of each frame to show the direction of the energy flow in relation to the peaks of the Airy beam. The flow of energy of the main peak at ξ>0 is consistently pointed at 45° relative to the s_x -s_y plane at all ξ-locations. In contrast, the direction of the energy flow for the Airy tails, or the peaks oriented along the horizontal or vertical axis approaches a direction perpendicular to that axis. The direction of the net energy flow is measured, however, to be constant and pointed in the direction that the main peak moves, i.e 45° or along the line s_x =s_y . This is in contrast to the locally varying direction of energy flow.

Fig. 2. The numerically calculated Poynting vector in the s_x -s_y plane, of the finite-energy Airy wave given by Eq. (6) with a=0.15 at ξ=0 (a), ξ=0.025 (b), ξ=0.050 (c), ξ=0.075 (d). The intensity of the Airy field is shown in the background of each frame.

It is interesting to note that in (a) of Fig. 2, the Poynting vector is initially pointing in the negative s_x or negative s_y -direction on each Airy tail; in (b) of the figure the direction starts to turn partially towards the direction of the energy flow of the main Airy peak (45°); and in (c) and (d) as the beam propagates further, the direction swings around even more towards 45°. The change in the Poynting vector can be used to explain why in Ref. [15

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of Airy beams,” Opt. Lett. 33, 207–209 (2008). [CrossRef] [PubMed]

] the optical Airy beam is said to be able to “ascend until it stalls due to downward acceleration” since S⃗∝P⃗ which points more along the tangential directions than the ξ-direction the further the beam propagates.

If the aperture function in Eq. (3) is ignored, i.e. setting a=0, then the Poynting vector would point along the line s_x =s_y everywhere. This leads to infinite energy propagation as discussed in section 1 above and in Ref. [5

G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]

]. The choice of aperture is partially responsible for the spatially varying Poynting vector yet does not seem to affect the calculations in the next section.

4. Angular momentum of the Airy beam

It is well known that as P⃗∝E⃗×B⃗ from which follows that angular momentum density about the $\hat{ξ}$ -axis is

j_{ξ} (s_{x}, s_{y}, ξ) = \vec{r} \times {〈 \vec{E} \times \vec{B} 〉}_{ξ}

= s_{x} \cdot S_{s_{y}} + s_{y} \cdot S_{s_{x}} .

(7)

The numerically computed <S⃗> taken from Eq. (6) and shown in Fig. 2 is used to calculate the angular momentum in the $\hat{ξ}$ -direction. Figure 3 (a-d) shows the $\hat{ξ}$ -component of the angular momentum density with a=0.15 at ξ=0.025, 0.050, 0.075, and 0.100, respectively. At ξ=0 the computed angular momentum is zero so it is not shown. In Fig. 3 reds are positive values (clockwise), blues are negative values (counter-clockwise), and green is zero. The non-discrete nature of these values will be discussed in the next section.

Fig. 3. The numerically calculated angular momentum in the $\hat{ξ}$ -direction of the finite-energy Airy wave given by Eq. (7) with a=0.15 at ξ=0.025 (a), ξ=0.050 (b), ξ=0.075 (c), ξ=0.100 (d). Reds are positive values, blues are negative values, and green is zero.

As the beam propagates, the net angular momentum about the ξ-axis is always zero. The spatial distribution of the angular momentum is changing however, and locally has non-zero values of angular momentum. Not only is the angular momentum changing in the Airy tails, but there are also changes to the angular momentum in the main Airy peak. This change of angular momentum is a torque that corresponds to the force present due to the changing linear momentum. Note that changing the axis about which the angular momentum is calculated would merely alter the angular momentum by a constant.

5. Angular momentum spectrum of Airy beams

Recently a new type of imaging was proposed that is based on the phase and spatial profile of the wavefront and is coined spiral imaging[9

L. Torner, J.P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13, 873 (2005). [CrossRef] [PubMed]

]. This spiral imaging is similar to what Ref. [14

M. V. Vasnetsov, J. P. Torres, D. V. Petrov, and L. Torner, “Observation of the orbital angular momentum spectrum of a light beam,” Opt. Lett. 28, 2285–2287 (2003). [CrossRef] [PubMed]

] refers to as the orbital angular momentum (OAM) spectrum of a beam of light that they show experimentally. One can use this technique to get a more complete picture of what the zero net angular momentum density with local non-discrete, non-zero values means.

Any optical beam can be decomposed into a superposition of angular harmonics in cylindrical coordinates written as

u (r, ϕ; z) = \frac{1}{\sqrt{2 π}} \sum_{m} a_{m} (r, z) \exp (im ϕ),

(8)

where a_m (r,z)=1/(2π)^1/2∫^2π ₀ u(r,ϕ,z)exp(-imϕ)dϕ and the energy of each mode, m, is described by C_m =∫^∞ ₀|a_m (r,ϕ,z)|² rdr. The power, or weight, of each angular momentum state for the arbitrary field u is given by [9

L. Torner, J.P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13, 873 (2005). [CrossRef] [PubMed]

]

P_{m} = \frac{C_{m}}{\sum_{- \infty}^{\infty} C_{n}} .

(9)

Figure 4(a-f) shows the angular momentum spectrum of the Airy beam where the field, u, is taken from Eq. (3), i.e. the weight of each spiral mode when the field is decomposed in these spiral harmonics, for ξ ranging from 0 to 0.125. Note that the sum of all of the modes at each ξ location is one and the net angular momentum is zero.

Fig. 4. The angular momentum spectrum of the finite-energy Airy wave given by eqs. (8, 9) with a=0.15 at ξ=0 (a), ξ=0.025 (b), ξ=0.05 (c), ξ=0.075 (d), ξ=0.100 (e), and ξ=0.125 (f).

The angular momentum, J_ξ , is shown to have non-discrete, non-integer values (positive and negative). Figure 4 shows that this Airy field, locally, has an integer sum of discrete values of orbital angular momentum while the total angular momentum is in fact zero.

6. Conclusions

In this paper we analyze the spatial evolution of the Airy solution to the paraxial wave equation and show that while momentum is changing, energy and momentum are conserved. Looking at the change of the magnitude and direction of the airy beam can give further insight into the dynamics of this class of solution to the paraxial wave equation. This technique can be further applied to Airy beams with different “launching angles” as described by Ref. [7

I. M. Besieris, A. M. Shaarawi, and R. W. Ziolkowski, “Nondispersive accelerating wave packets,” Am. J. Phys. 62, 519–521 (1994). [CrossRef]

] which discussed the ballistic dynamics of Airy beams. When the “launching angle” of the Airy beam is varied the quantity of interest would be the change in the direction of the energy flow.

We point out that the linear and angular momenta are changing as the Airy beam propagates in the ξ-direction which should have implications when analyzing the velocity of this field. Some form of the velocity, be it phase, energy, or signal velocity, should be changing as the beam propagates and should be investigated in future studies. These beams have promise for applications in optical trapping, imaging, and spectroscopy where a sample might interact with a changing momentum and spatially varying angular momentum.

Acknowledgments

This work was supported in part by organized research at CCNY, NASA URC - Center for Optical Sensing and Imaging (COSI) at CCNY (NASA Grant No.: NCC-1-03009), and DOD Center of Nanoscale Photonics at CCNY. HIS is grateful to Matthias Lenzner for helpful discussions.

References and links

1.	L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX, pp. 291–372 (1999). [CrossRef]
2.	J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58, 1499–1501 (1987). [CrossRef] [PubMed]
3.	J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4, 651 (1987). [CrossRef]
4.	M. V. Berry and N. L. Balazs, “Nonspreading wave packets,” Am. J. Phys. 47, 264 (1979). [CrossRef]
5.	G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Observation of Accelerating Airy Beams,” Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]
6.	G. A. Siviloglou and D. N. Christodoulides, “Accelerating finite energy Airy beams,” Opt. Lett. 32, 979 (2007). [CrossRef] [PubMed]
7.	I. M. Besieris, A. M. Shaarawi, and R. W. Ziolkowski, “Nondispersive accelerating wave packets,” Am. J. Phys. 62, 519–521 (1994). [CrossRef]
8.	D. M. Greenberger, “Comment on ‘Non-Spreading Wave Packets’,” Am. J. Phys. 48, 256 (1980). [CrossRef]
9.	L. Torner, J.P. Torres, and S. Carrasco, “Digital spiral imaging,” Opt. Express 13, 873 (2005). [CrossRef] [PubMed]
10.	J. D. Jackson, Classical Electrodynamics , (Wiley, New York, 1962).
11.	M. Born and E. Wolf, Principles of Optics , 7th Edition, (Cambridge University Press, Cambridge, 1999).
12.	L. Allen and M. J. Padgett, “The Poynting vector in Laguerre-Gaussian beams and the interpretation of their angular momentum desnity,” Opt. Commun. 184, 67–71 (2000). [CrossRef]
13.	H. I. Sztul and R. R. Alfano, “Double-slit interference with Laguerre-Gaussian beams,” Opt. Lett. 31, 999–1001 (2006). [CrossRef] [PubMed]
14.	M. V. Vasnetsov, J. P. Torres, D. V. Petrov, and L. Torner, “Observation of the orbital angular momentum spectrum of a light beam,” Opt. Lett. 28, 2285–2287 (2003). [CrossRef] [PubMed]
15.	G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, “Ballistic dynamics of Airy beams,” Opt. Lett. 33, 207–209 (2008). [CrossRef] [PubMed]

OCIS Codes
(260.2110) Physical optics : Electromagnetic optics
(350.5500) Other areas of optics : Propagation

ToC Category:
Physical Optics

History
Original Manuscript: April 17, 2008
Revised Manuscript: May 1, 2008
Manuscript Accepted: June 2, 2008
Published: June 11, 2008

Citation
H. I. Sztul and R. R. Alfano, "The Poynting vector and angular momentum of Airy beams," Opt. Express 16, 9411-9416 (2008)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-16-13-9411

Sort: Author | Year | Journal | Reset

References

L. Allen, M. J. Padgett, and M. Babiker, Progress in Optics XXXIX (Elsevier Science & Technology, 1999) pp. 291-372. [CrossRef]
J. Durnin, J. J. Miceli, and J. H. Eberly, "Diffraction-free beams," Phys. Rev. Lett. 58, 1499-1501 (1987). [CrossRef] [PubMed]
J. Durnin, "Exact solutions for nondiffracting beams. I. The scalar theory," J. Opt. Soc. Am. A 4, 651 (1987). [CrossRef]
M. V. Berry and N. L. Balazs, "Nonspreading wave packets," Am. J. Phys. 47, 264 (1979). [CrossRef]
G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Observation of Accelerating Airy Beams," Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]
G. A. Siviloglou and D. N. Christodoulides, "Accelerating finite energy Airy beams," Opt. Lett. 32, 979 (2007). [CrossRef] [PubMed]
I. M. Besieris, A. M. Shaarawi, and R. W. Ziolkowski, "Nondispersive accelerating wave packets," Am. J. Phys. 62, 519-521 (1994). [CrossRef]
D. M. Greenberger, "Comment on ???Non-Spreading Wave Packets???, " Am. J. Phys. 48, 256 (1980). [CrossRef]
L. Torner, J. P. Torres, and S. Carrasco, "Digital spiral imaging," Opt. Express 13, 873-881 (2005). [CrossRef] [PubMed]
J. D. Jackson, Classical Electrodynamics (Wiley, New York, 1962).
M. Born and E. Wolf, Principles of Optics, 7th Edition, (Cambridge University Press, Cambridge, 1999).
L. Allen and M. J. Padgett, "The Poynting vector in Laguerre-Gaussian beams and the interpretation of their angular momentum desnity," Opt. Commun. 184, 67-71 (2000). [CrossRef]
H. I. Sztul and R. R. Alfano, "Double-slit interference with Laguerre-Gaussian beams," Opt. Lett. 31, 999-1001 (2006). [CrossRef] [PubMed]
M. V. Vasnetsov, J. P. Torres, D. V. Petrov, and L. Torner, "Observation of the orbital angular momentum spectrum of a light beam," Opt. Lett. 28, 2285-2287 (2003). [CrossRef] [PubMed]
G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Ballistic dynamics of Airy beams," Opt. Lett. 33, 207-209 (2008). [CrossRef] [PubMed]

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.

Click here to see a list of articles that cite this paper