## 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, diffraction-free propagation , self reconstruction , 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 . 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, 15]. The Airy beam is said to be free of diffractive spreading and freely accelerating in the absence of any external potential .

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 . 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,

$2i∂u(sx,sy,ξ)∂ξ+∇τ2u(sx,sy,ξ)=0$

where sx=x/x 0 and sy=y/y 0 are normalized transverse coordinates, ξ=z/k(x 0+y 0)2 is a normalized propagation distance, x 0 and y 0 are normalization constants, k=2πn0 and ∇2 τ is the second partial derivative in the transverse direction.

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

$u(sx,sy,ξ)=Ai(sx−(ξ2)2)Ai(sy−(ξ2)2)exp⁡[i(sxξ2+syξ2)−i(ξ36)],$

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 .

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, 6], 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(sx,sy,0)=Ai(sx)Ai(sy)exp[a(sx+sy)], where a<<1 and is a positive parameter that limits the infinite energy in the Airy tail, the electric field amplitude becomes

$u(sx,sy,ξ)=Ai(sx−(ξ2)2+iaξ)Ai(sy−(ξ2)2+iaξ)×$
$exp⁡[asx+asy−(aξ22)−i(ξ36)−i(a2ξ2)+iξ(sx+sy)2].$

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

$xm=λo2ξ216π2xo2andym=λo2ξ216π2yo2.$ Fig. 1. The intensity of the finite-energy Airy wave in the sx-sy 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, 11]. 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, 12, 13].

The Poynting vector is defined as :

$S→=(c4π)E→×B→,$

where c is the speed of light. Given a vector potential A⃗=$ε^$ u(sx,sy,ξ)exp[ik 2 ξ(x 0+y 0)], where $ε^$ is an arbitrary polarization and u(sx,sy,ξ) 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. , to calculate the time-averaged Poynting vector, <S⃗>. Assuming an -polarized field, <S⃗> becomes 

$c4π=c8π(E→×B*→+E*→×B→)$
$=c8π[iω(u∇τu*−u*∇τu)+2ωk∣u∣2ξ̂].$

The $ξ^$ term in the above equation is the energy flow in the $ξ^$-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 $ξ^$-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 sx-sy 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 sx=sy. This is in contrast to the locally varying direction of energy flow. Fig. 2. The numerically calculated Poynting vector in the sx-sy 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 sx or negative sy-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.  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 sx=sy everywhere. This leads to infinite energy propagation as discussed in section 1 above and in Ref. . 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 $ξ^$-axis is

$jξ(sx,sy,ξ)=r→×〈E→×B→〉ξ$
$=sx·Ssy+sy·Ssx.$

The numerically computed <S⃗> taken from Eq. (6) and shown in Fig. 2 is used to calculate the angular momentum in the $ξ^$-direction. Figure 3 (a-d) shows the $ξ^$-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 $ξ^$-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. This spiral imaging is similar to what Ref.  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)=12π∑mam(r,z)exp⁡(imϕ),$

where am(r,z)=1/(2π)1/22π 0 u(r,ϕ,z)exp(-imϕ) and the energy of each mode, m, is described by Cm=∫ 0|am(r,ϕ,z)|2 rdr. The power, or weight, of each angular momentum state for the arbitrary field u is given by 

$Pm=Cm∑−∞∞Cn.$

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.  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.

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]

### References

• View by:
• |
• |
• |

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]

#### 2007 (2)

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

#### 2000 (1)

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]

#### 1994 (1)

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

#### 1987 (2)

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

#### 1980 (1)

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

#### 1979 (1)

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

#### Allen, L.

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]

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

#### Babiker, M.

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

#### Balazs, N. L.

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

#### Berry, M. V.

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

#### Besieris, I. M.

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

#### Born, M.

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

#### Broky, J.

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

#### Christodoulides, D. N.

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

#### Dogariu, A.

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

#### Durnin, J.

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

#### Eberly, J. H.

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

#### Greenberger, D. M.

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

#### Jackson, J. D.

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

#### Miceli, J. J.

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

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]

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

#### Shaarawi, A. M.

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

#### Siviloglou, G. A.

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

#### Wolf, E.

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

#### Ziolkowski, R. W.

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

#### Am. J. Phys. (3)

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

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]

#### Opt. Commun. (1)

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]

#### Phys. Rev. Lett. (2)

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

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

#### Other (3)

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

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

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

### Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

### Figures (4)

Fig. 1. The intensity of the finite-energy Airy wave in the sx -sy 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).
Fig. 2. The numerically calculated Poynting vector in the sx -sy 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.
Fig. 3. The numerically calculated angular momentum in the $ξ ^$ -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.
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).

### Equations (12)

$2 i ∂ u ( s x , s y , ξ ) ∂ ξ + ∇ τ 2 u ( s x , s y , ξ ) = 0$
$u ( s x , s y , ξ ) = Ai ( s x − ( ξ 2 ) 2 ) Ai ( s y − ( ξ 2 ) 2 ) exp ⁡ [ i ( s x ξ 2 + s y ξ 2 ) − i ( ξ 3 6 ) ] ,$
$u ( s x , s y , ξ ) = Ai ( s x − ( ξ 2 ) 2 + i a ξ ) Ai ( s y − ( ξ 2 ) 2 + i a ξ ) ×$
$exp ⁡ [ a s x + a s y − ( a ξ 2 2 ) − i ( ξ 3 6 ) − i ( a 2 ξ 2 ) + i ξ ( s x + s y ) 2 ] .$
$x m = λ o 2 ξ 2 16 π 2 x o 2 and y m = λ o 2 ξ 2 16 π 2 y o 2 .$
$S → = ( c 4 π ) E → × B → ,$
$c 4 π < E → × B → > = c 8 π ( E → × B * → + E * → × B → )$
$= c 8 π [ i ω ( u ∇ τ u * − u * ∇ τ u ) + 2 ω k ∣ u ∣ 2 ξ ̂ ] .$
$j ξ ( s x , s y , ξ ) = r → × 〈 E → × B → 〉 ξ$
$= s x · S s y + s y · S s x .$
$u ( r , ϕ ; z ) = 1 2 π ∑ m a m ( r , z ) exp ⁡ ( im ϕ ) ,$
$P m = C m ∑ − ∞ ∞ C n .$