OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 25 — Dec. 7, 2009
  • pp: 22432–22441
« Show journal navigation

Nonparaxial diffraction analysis of Airy and SAiry beams

Luis Carretero, Pablo Acebal, Salvador Blaya, Celia Garcia, Antonio Fimia, Roque Madrigal, and Angel Murciano  »View Author Affiliations


Optics Express, Vol. 17, Issue 25, pp. 22432-22441 (2009)
http://dx.doi.org/10.1364/OE.17.022432


View Full Text Article

Acrobat PDF (552 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We theoretically analyze Airy beams by solving the exact vectorial Helmholtz equation using boundary conditions at a diffraction aperture. As result, the diffracted beams are obtained in the whole space; thus, we demonstrate that the parabolic trajectories are larger than those previously reported, showing that the Airy beams start to form before the Fourier plane. We also demonstrate the possibility of using a new type of Airy beams (SAiry beams) with finite energy that can be generated at the focal plane of the lens due to diffraction by a circular aperture of a spherical wave modified by a cubic phase. The finite energy ensured by the principle of conservation of energy of a diffracted beam.

© 2009 Optical Society of America

1. Introduction

Since Airy optical beams were first observed [1

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

, 2

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

] intense research has been carried out on this finite-energy solution of the paraxial Helmholtz equation. For example, based on the characteristics of the intensity pattern of Airy light beams, optical micromanipulation has recently been demonstrated [3

3. J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using Airy wavepackets,” Nat. Photonics 2, 675–678 ( 2008). [CrossRef]

, 4

4. J. Baumgartl, G. Hannappel, D. J. Stevenson, D. Day, M. Gu, and K. Dholakia, “Optical redistribution of microparticles and cells between microcells,” Lab on a Chip 9, 1334–1336 ( 2009). [CrossRef] [PubMed]

]. So due to the possibility of removing microparticles and cells from a section of a sample chamber, in these papers new perspectives were opened in the field of redistribution of particles and cells between different media. Furthermore, using femtosecond Airy beams, curved plasma channels have been experimentally shown to be generated in air, and these are parabolic for the dominant intensity of the beam [5

5. P. Polynkin, M. Koleskik, J. V. Moloney, G. A. Siviloglou, and D. N. Christodoulides, “Curved plasma channels generation using ultraintense Airy beams,” Science 324, 229–232 ( 2009). [CrossRef] [PubMed]

]. Moreover, a new way of generating Airy beams using three-wave mixing processes in asymmetric nonlinear photonic crystals opens up new possibilities for all-optical switching and manipulation of Airy beams [6

6. T. Ellenbogen, N. Voloch-Bloch, A. Ganany-Padowicz, and A. Arie, “Nonlinear generation and manipulation of Airy beams,” Nat. Photonics 3, 395–398 ( 2009). [CrossRef]

]. In this sense, this methodology allows the creation of Airy beams at new wavelengths and high intensities, which is impossible with the conventional (linear) methods.

Other interesting properties of these beams have been described such as the ballistic dynamics of the beams akin to projectiles moving under the action of gravity [7

7. 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 self restoration properties (self-healing) after passing small obstacles or propagating in adverse environments [8

8. J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, “Self-healing properties of optical Airy beams,” Opt. Express 16, 12880–12891 ( 2008). [CrossRef] [PubMed]

], the transformation by paraxial optical systems (that can be described by four complex parameters) [9

9. M. A. Bandres and J. Gutierrez-Vega, “Airy-Gauss beams and their transformation by paraxial optical systems,” Opt. Express 15, 16719–16728 ( 2007). [CrossRef] [PubMed]

], or the evolution of the Poynting vector and angular momentum when the beam propagates through space [10

10. H. Sztul and R. Alfano, “The Poynting vector and angular momentum of Airy beams,” Opt. Express 16, 9411–9416 ( 2008). [CrossRef] [PubMed]

].

Fig. 1. Schematic diagram of the generation of Airy beams.

Thus, in this paper in contrast to other studies, we use the nonparaxial vectorial diffraction integrals, which are an exact solution of the Maxwell equations and the complete Helmholtz wave equation given by:

2Ex2+2Ey2+2Ez2+k2E=0
(1)

The exact vectorial Helmholtz equation is solved using the boundary conditions on a diffraction aperture situated at z=0 plane (see Fig. 1). Therefore, using this methodology several effects (that in the paraxial Helmholtz equation solution are not obtained) can be analyzed such as the propagation in the whole space. Moreover, in this work we propose the existence of a new family of Airy beams (SAiry beams) with finite energy that can be generated at the focal plane of the lens due to diffraction by a circular aperture of a spherical wave modified by a cubic phase. For the SAiry beams, the finite energy is ensured by principle of conservation energy of a diffracted beam.

2. Nonparaxial vectorial diffraction integrals

The solution to the Maxwell equations (and as consequence the Helmholtz equation [Eq. (1)]), valid through the half-space z>0 (see Fig. 1), based on the knowledge of boundary conditions Eo on the plane z=0:

EoxoyoO=Exxoyo0x̂+Eyxoyo0ŷ
(2)

is given by [11

11. R. K. Luneburg, Mathematical theory of optics (University California Press, Berkeley, California, 1964).

]:

Exx,yz=z2πExxoyo0ikR1R3exp(ikR)dxodyo
Eyxyz=z2πEyxoyo0ikR1R3exp(ikR)dxodyo
Ez(x,y,z)=xzExxyzyzEyxyz
12π(xoExxoyo0+yoEyxoyo0)=ikR1R3exp(ikR)dxodyo
(3)

where k is the wavenumber and R=[(x-xo)2+(y-yo)2+z 2]1/2. The propagation properties of vectorial beams can be studied using the Eqs. (2) and (3). Assuming that after the lens (for analysis of experimental set-up see for example references [1

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

, 3

3. J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using Airy wavepackets,” Nat. Photonics 2, 675–678 ( 2008). [CrossRef]

, 4

4. J. Baumgartl, G. Hannappel, D. J. Stevenson, D. Day, M. Gu, and K. Dholakia, “Optical redistribution of microparticles and cells between microcells,” Lab on a Chip 9, 1334–1336 ( 2009). [CrossRef] [PubMed]

]) we obtain an on-axis polarized, monochromatic spherical wave modified by an aperture function given by: t(xo,yo)={exp(iβ(xo3+yo3))exp(xo2+yo2wo2)if(xo,yo)D0else where a is the aperture radius and D={(xo,yo)\(x 2 o+y 2 o)≤ a 2}.

The field at the plane after the aperture may be represented as:

Eoxoyo0=t(xo,yo)exp(ik(xo2+yo2+zo2)12)(xo2+yo2+zo2)12[cos(α)x̂+sin(α)ŷ]
(4)

where zo=f2a2,, f is the focal length of the lens and α is the linearly polarized angle. In this work, we are going to analyze two types of diffraction patterns, the first one corresponds to those previously reported (see for example [1

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

]) in which β≠0 and wo=finite (Airy beam) and the second one which is characterized by β≠0 and wo→∞ (SAiry beam). In both cases we assumed diffraction by a finite aperture as can be seen in the transmittance function.

The intensity distribution at any point (x,y) of a plane z>0 is given by:

Ixyz=Exyz2=Exxyz2+Eyxyz2+Ezxyz2
(5)

The principle of conservation of energy in the solution to Maxwell’s equations ensures that total intensity must be equal and finite at all z-planes.

It is important to remark, that, from the analysis of Eq. (3), all the diffracted patterns on the half-space z>0 can be analyzed, including pre-focus (PRF) and post-focus (PF) regions (see Fig. 1). In this way, our method provides more complete information on the propagation characteristics of Airy beams which have not been previously analyzed due to the limitation of the paraxial Helmholtz equation (which can only be used at the PF region).

3. Numerical results and discussion

In order to analyze the propagation characteristics of the beams studied, we focus our attention on an electric field Eo(xo,yo,0) which is x̂ polarized, so α=0, and Eq. (3) can be reduced to:

Exxyz=z2πDEoxoyo0ikR1R3exp(ikR)dxodyo
Ez(x,y,z)=xzExxyz12π
DxoEoxoyo0ikR1R3exp(ikR)dxodyo
(6)

The parameters that we used in our numerical simulations are shown in Table 1 and were selected according to the experimental conditions mentioned in Ref. [1

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

].

Table 1. Values of parameters used in the numerical simulations

table-icon
View This Table

3.1. Airy beams

These numerical simulations of Airy beams can be compared with the experimental and theoretical results reported in [1

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

]. Thus, Figs. 2, 3 show the numerical simulations of diffraction patterns at different z planes in the PRF and PF zones for the Airy beam with parameters described in Table 1. The intensities have been normalized in such a way that the maximum value of intensity of the Airy beam at the focal plane is equal to 1.

As can be seen in Figs. 2(f)–2(h) an interesting result of the numerical simulations is that the Airy beam starts before the focal plane at the PRF region. This implies that the accelerating diffraction-free region that could be used in applications is 14 cm longer than predicted by the paraxial solution. These theoretical results justify the experimental measurements described in Ref. [5

5. P. Polynkin, M. Koleskik, J. V. Moloney, G. A. Siviloglou, and D. N. Christodoulides, “Curved plasma channels generation using ultraintense Airy beams,” Science 324, 229–232 ( 2009). [CrossRef] [PubMed]

], where plasma channels start before the Fourier plane. Furthermore, it is important to remark that the behavior of the Airy beam is not symmetrical to the focal plane, indeed at the PRF region the diffraction free distance of the Airy beam is nearly 60 % of the PF region. In order to compare the numerical results with those obtained analytically using the paraxial Helmholtz solutions, Fig. 4 shows the normalized intensities obtained at different Z=(z-f) planes of the PF region. As can be observed when Z<19 cm the intensity patterns are identical in both methods. However (since the condition of paraxial approximation 2Ez2<<k2E is not fulfilled), at Z>19 cm differences may be seen between the two method, such as the intensity distribution and the maximum intensity reached. However, these differences may greater when the a/f ratio tends to one, as was recently noted for spherical waves [1

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

].

Fig. 2. Normalized intensity distribution for Airy beams at different Z=(z-f) planes of PRF region.

In Fig. 5 the maximum values of intensity obtained at the PRF and PF regions are shown. As can be seen, the behavior of the maximum intensities is not symmetrical as regards the focal plane, and the slope of the growth zone in the PRF region is different to the decreasing slope in the PF region. Finally in Fig. 6 we show the parabolic trajectory of the Airy beam, where r=xmax2+ymax2 is represented as a function of the propagation distance Z, and (xmax,ymax) are the Cartesian coordinates of the maximum intensity at each plane z. As can be seen the parabolic trajectory is observed in the -0.2 m 0.2 m interval, and there is no acceleration outside this region. It is important to remark that outside the -0.15 m, 0.15 m z interval the maximum intensity is lower than 80% of the maximum obtained at the focal plane, so the effects of the peak intensity for several applications could may be negligible although the trajectory is still parabolic in the -0.2 m 0.2 m interval as mentioned above. For values of |Z|>0.2 m the parabolic trajectory disappears.

3.2. SAiry beams

In the second part of this study we present the result of the numerical simulations for the SAiry beams. The proposal of this new type of beams is based on a simple idea. Since the Gaussian amplitude modulation used is selected to ensure the finite energy at different propagation planes, this can also be fulfilled by using an aperture that limits the incident wave and solving Maxwell’s equations with these boundary conditions. Figures 7 and 8 show the numerical simulations of diffraction patterns at different Z planes in the PRF and PF zones for the SAiry beam parameters described in Table 1. Since we did not use a Gaussian beam the amplitude of transmittance is equal to 1. Moreover, the intensities were normalized using the same conditions as for Airy beam. As result of this normalization, the maximum intensities are higher than 1 at the focal plane. This is a result of energy conservation because the total intensity at Z=0 plane is higher if |t(xo,yo)|=1onD than if t(xo,yo)=exp(xo2+yo2wo2)onD..

Fig. 3. Normalized intensity distribution for Airy beams at different Z=(z-f) planes of PF region.

From these figures it may be deduced that accelerating beams can be obtained introducing only a cubic phase in a spherical wave. The differences in the diffraction patterns shown in Figs. 78 and 2–3 are exclusively due to the transmittance values (see Fig. 9) which modify the weight of the electric field Eo(xo,yo,o) in the integral 7. As can be seen the most significant difference between the transmittance that generates the Airy and that which generates the SAiry beams is that the high and low spatial frequencies have the same weight in SAiry beams, whereas the high spatial frequencies are nearly null for the Airy beams, and the contrast is lower for SAiry beams in the other of spatial frequencies.

As can be seen in Figs. 78 the region where the diffraction pattern of SAiry beams can be observed varies from -23 to 38 cm and in all cases the intensity of the main peak is higher than that the obtained in the same plane with an Airy beam as can be clearly observed in Fig. 5. Moreover it is important to point out that the behavior of the intensity is quite different between the two kinds of beams, for example, the maximum intensity is not obtained at the focal plane, but at a plane situated at -14 cm from the focus. Also the growth and decreasing slopes are very different in PRF and PF regions. Figure 6 shows the parabolic trajectory of the SAiry beam, which as can be seen is greater than that the observed for the Airy beam and approximately coincides in the -20, 20 cm Z region. Moreover, an observed drawback with SAiry beams is that the spread of the diffraction pattern is always wider than for Airy beams. In any case, from Figs. 5 and 6 it may be deduced that for some applications where the spread of the diffraction pattern is not important and the peak maximum and parabolic trajectory are the main parameters, SAiry beams are better than Airy beams.

Fig. 4. Normalized intensity distribution for Airy beams at different Z=(z-f) planes of the PF region numerically calculated (a),(c),(e),(g) and analytically obtained by means of the paraxial Helmholtz solutions [1],(b),(d),(f), (h).
Fig. 5. Normalized maximum intensities of Airy beams (blue) and SAiry Beams (red) at PRF and PF regions
Fig. 6. Transverse acceleration of Airy beams (blue) and SAiry Beams (red) at PRF and PF regions.
Fig. 7. Normalized intensity distribution for SAiry beams at different planes Z=(z-f) of PRF region.

4. Conclusions

We have analyzed Airy beams by solving the exact vectorial Helmholtz equation using the boundary conditions at a diffraction aperture and as a consequence, the diffracted beams were obtained in the whole space. The proposed method may be very useful if a/f parameters become close to unity or for the study of different types of linear polarized beams. We demonstrate theoretically that the parabolic trajectories are larger that those previously reported using the paraxial Helmholtz equation, showing that the Airy beams start to form before the Fourier plane. We also propose the possibility of using a new type of Airy beams (SAiry beams) with finite energy that can be generated at the focal plane of the lens due to diffraction by a circular aperture of a spherical wave modified by a cubic phase, and the finite energy is ensured by principle of conservation of energy of a diffracted beam. These SAiry beams may be a good choice for some applications where the spread of the diffraction pattern is not important but the peak maximum and the parabolic trajectory are the main parameters to be considered.

Fig. 8. Normalized intensity distribution for SAiry beams at different planes Z=(z-f) of PF region.
Fig. 9. Real (a)–(b) and Imaginary (c)–(d) part of t(xo,yo) for SAiry and Airy beams respectively

Acknowledgments

The authors acknowledge support from project FIS2009-11065 of Ministerio de Ciencia e Innovación of Spain.

References and links

1.

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

2.

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

3.

J. Baumgartl, M. Mazilu, and K. Dholakia, “Optically mediated particle clearing using Airy wavepackets,” Nat. Photonics 2, 675–678 ( 2008). [CrossRef]

4.

J. Baumgartl, G. Hannappel, D. J. Stevenson, D. Day, M. Gu, and K. Dholakia, “Optical redistribution of microparticles and cells between microcells,” Lab on a Chip 9, 1334–1336 ( 2009). [CrossRef] [PubMed]

5.

P. Polynkin, M. Koleskik, J. V. Moloney, G. A. Siviloglou, and D. N. Christodoulides, “Curved plasma channels generation using ultraintense Airy beams,” Science 324, 229–232 ( 2009). [CrossRef] [PubMed]

6.

T. Ellenbogen, N. Voloch-Bloch, A. Ganany-Padowicz, and A. Arie, “Nonlinear generation and manipulation of Airy beams,” Nat. Photonics 3, 395–398 ( 2009). [CrossRef]

7.

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

8.

J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, “Self-healing properties of optical Airy beams,” Opt. Express 16, 12880–12891 ( 2008). [CrossRef] [PubMed]

9.

M. A. Bandres and J. Gutierrez-Vega, “Airy-Gauss beams and their transformation by paraxial optical systems,” Opt. Express 15, 16719–16728 ( 2007). [CrossRef] [PubMed]

10.

H. Sztul and R. Alfano, “The Poynting vector and angular momentum of Airy beams,” Opt. Express 16, 9411–9416 ( 2008). [CrossRef] [PubMed]

11.

R. K. Luneburg, Mathematical theory of optics (University California Press, Berkeley, California, 1964).

12.

J. Guo, X. Zhao, and Y. Min, “The general integral expressions for on-axis nonparaxial vectorial spherical waves diffracted at a circular aperture,” Opt. Comm. 282, 1511–1515 ( 2009). [CrossRef]

OCIS Codes
(050.1940) Diffraction and gratings : Diffraction
(260.2110) Physical optics : Electromagnetic optics
(350.5500) Other areas of optics : Propagation

ToC Category:
Physical Optics

History
Original Manuscript: July 27, 2009
Manuscript Accepted: October 28, 2009
Published: November 23, 2009

Citation
Luis Carretero, Pablo Acebal, Salvador Blaya, Celia García, Antonio Fimia, Roque Madrigal, and Angel Murciano, "Nonparaxial diffraction analysis of Airy and SAiry beams," Opt. Express 17, 22432-22441 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-25-22432


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Observation of accelerating Airy Beams," Phys. Rev. Lett. 99, 213901 (2007). [CrossRef]
  2. G. A. Siviloglou and D. N. Christodoulides, "Accelerating finite energy Airy beams," Opt. Lett. 32, 979-981 (2007). [CrossRef] [PubMed]
  3. J. Baumgartl, M. Mazilu, and K. Dholakia, "Optically mediated particle clearing using Airy wavepackets," Nat. Photonics 2, 675-678 (2008). [CrossRef]
  4. J. Baumgartl, G. Hannappel, D. J. Stevenson, D. Day, M. Gu, and K. Dholakia, "Optical redistribution of microparticles and cells between microcells," Lab Chip 9, 1334-1336 (2009). [CrossRef] [PubMed]
  5. P. Polynkin, M. Koleskik, J. V. Moloney, G. A. Siviloglou, and D. N. Christodoulides, "Curved plasma channels generation using ultraintense Airy beams," Science 324, 229-232 (2009). [CrossRef] [PubMed]
  6. T. Ellenbogen, N. Voloch-Bloch, A. Ganany-Padowicz, and A. Arie, "Nonlinear generation and manipulation of Airy beams," Nat. Photonics 3, 395-398 (2009). [CrossRef]
  7. G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Ballistic dynamics of Airy beams," Opt. Lett. 33, 207-209 (2008). [CrossRef] [PubMed]
  8. J. Broky, G. A. Siviloglou, A. Dogariu, and D. N. Christodoulides, "Self-healing properties of optical Airy beams," Opt. Express 16, 12880-12891 (2008). [CrossRef] [PubMed]
  9. M. A. Bandres and J. Gutierrez-Vega, "Airy-Gauss beams and their transformation by paraxial optical systems," Opt. Express 15, 16719-16728 (2007). [CrossRef] [PubMed]
  10. H. Sztul and R. Alfano, "The Poynting vector and angular momentum of Airy beams," Opt. Express 16, 9411-9416 (2008). [CrossRef] [PubMed]
  11. R. K. Luneburg, Mathematical theory of optics (University California Press, Berkeley, California, 1964).
  12. J. Guo, X. Zhao, and Y. Min, "The general integral expressions for on-axis nonparaxial vectorial spherical waves diffracted at a circular aperture," Opt. Comm. 282, 1511-1515 (2009). [CrossRef]

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.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited