Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

Optik 114, No. 11 (2003) 504 Daomu Zhao504–508 et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems http...

155KB Sizes 0 Downloads 37 Views

Optik 114, No. 11 (2003) 504 Daomu Zhao504–508 et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems http://www.elsevier-deutschland.de/ijleo International Journal for Light and Electron Optics

Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems Daomu Zhao1, Haidan Mao1, Hongjie Liu1, 2, Feng Jing2, Qihua Zhu2, Xiaofeng Wei2 1 2

Department of Physics, Zhejiang University, Hangzhou 310027, China Research Center of Laser Fusion, CAEP, Mianyang 621900, China

Abstract: The apertured fractional Fourier transform (FRFT) system is introduced and applied to treat the propagation of Hermite-Gaussian beams. Based on the treatment that a rectangular function can be expanded into a approximate sum of complex Gaussian functions with finite numbers, the analytical expressions for the output field distribution of a Hermite-Gaussian beam through an apertured FRFT system are obtained and compared with those obtained from numerically integral calculation. The results show that our method can significantly improve the numerical calculation efficiency. Key words: Hermite-Gaussian beams – fractional Fourier transform – aperture – diffraction integral

1. Introduction The fractional Fourier transform (FRFT) is regarded as a generalization of the conventional Fourier transform. In the 1980s Namias [1], McBride and Kerr [2] took up the FRFT for the benefit of quantum mechanics, but did not gain much attention until its reintroduction to optics in 1993 by Mendlovic and Ozaktas [3, 4] and Lohmann [5]. Since then a lot of works have been done on its properties, optical implementations, and applications [6–20]. For instance, the FRFT has been used as a new method in signal processing, physical optics, laser optics, etc. From the optical point of view there exist two different approaches to obtain definitions of a FRFT. One definition is based on Hermite-Gaussian-modes and can be implemented optically by means of GRIN media [3, 4]. The other approach is based on rotating the Wigner distribution function of the input signal by a certain angle and can be implemented by means of bulk optical setups [5]. Both definitions are essentially equivalent [9]. The two kinds of optical setups for im-

plementing p order FRFT suggested by Lohmann are shown in fig. 1, respectively. Where fs is the standard focal length and f ¼ pp=2. Generally, the FRFT system without aperture had been extensively studied and used. The apertured FRFT system, to the best of our knowledge, has been paid little attention. Practically, aperture always exists in FRFT systems, for example, the finite size of lens. Therefore study on the behavior of light propagation through the apertured FRFT systems would be of practical interest. The paper is organized as follows. The theoretical analyses of a Hermite-Gaussian beam through the two types of apertured FRFT systems are given in section 2. Some numerical simulation comparisons by using the derived analytical formulae and the diffraction integral formulae are given in section 3. Finally, a simple conclusion is outlined in section 4.

2. Approximate analytical formulae for propagation of Hermite-Gaussian beams through apertured FRFT systems In this paper, we only consider the case of the finite size of lens in FRFT systems. For the sake of simplicity, we just consider the one-dimensional case, the extension to the two-dimensional one is straightforward. In case of the type I shown in fig. 1a, the whole system can be separated into two parts, the first is the free propagation from the input plane to z1 plane, and the second is the apertured propagation from z1 plane to the output plane because of the finite size of lens. For the first one we have  Eðx; z1 Þ ¼

Received 18 August 2003; accepted 5 November 2003. Correspondence to: D. Zhao Fax: ++86-571-88863887 E-mail: [email protected]

1  ð i 1=2 Eðx1 ; 0Þ ld 1

  ip 2 ðx1  2x1 x þ x2 Þ dx1  exp  ld

ð1Þ

0030-4026/03/114/11-504 $ 15.00/0

505

Daomu Zhao et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

f = f s / sin φ

d = f s tan

Fig. 1. The two optical setups for implementing FRFT. a) type I, b) type II.

0

φ

d = f s tan

2

f = f s / tan

φ

d = f s sin φ

Z1

Z

Z

0

i ld

1=2 ða

(b)

ing eqs. (1) and (5) into (4) and using the following integral formula

and for the second one we have 

Eðx; z1 Þ

1 ð

a



ip  exp  ld

   d 2 2 1 dx ; x  2xx2 þ x2 f ð2Þ

where a means the half-width of the lens aperture and the constant phase terms in eqs. (1) and (2) have been omitted which have no influence on the output intensity distribution. Introducing the aperture function  1; jxj  a Ap ðxÞ ¼ ð3Þ 0; jxj > a

"

# ðx  yÞ2 exp  Hm ðxÞ dx 2u

1

¼ ð2puÞ1=2 ð1  2uÞm=2 Hm ½yð1  2uÞ1=2  ;



i ld

1=2 1 ð 

Eðx2 ; zÞ ¼

 exp 

ip ld

 1

  d 2 dx : x  2xx2 þ x22 f ð4Þ

Generally, the aperture function can be expanded as the sum of complex Gaussian functions with finite numbers and is given by   N P Bn 2 An exp  2 x ; ð5Þ Ap ðxÞ ¼ a n¼1 where An and Bn are the expansion and Gaussian coefficients, respectively, which could be obtained by optimization-computation directly [21]. The Hermite-Gaussian beam is an important type of laser beams. Let us consider a Hermite-Gaussian beam pffiffiffi represented by exp ðx2 =W02 Þ Hm ð 2x=W0 Þ is incident on the input plane of the apertured FRFT system, where W0 is beam waist size of Gaussian part. Insert-

N P

 An

n¼1

1 Cn

1=2  1

B2 Cn

m=2

    ip p2 þ 2  exp  x2 ld l d2 Cn 2 "  1=2 # iBpx2 B2 1 ;  Hm ldCn Cn

Eðx; z1 Þ Ap ðxÞ

1

ð6Þ

after performing tedious integration, an analytical expression of the output field distribution in FRFT plane is obtained

then eq. (2) becomes Eðx2 ; zÞ ¼

2

2

(a)

Eðx2 ; zÞ ¼

φ

ð7Þ

where B¼

pffiffiffi  1=2 i 2p 2ld   1 ; 1 ip ld þ ipW02 þ ldW0 W02 ld

Cn ¼

Bn ip þ þ a2 ld

l2 d 2

p2 ip þ 1 ip ld þ W02 ld



 1

ð8Þ

 d : f (9)

In case of the type II shown in fig. 1b, there exist two apertures on the respective positions of the lenses, in fact the second aperture only truncates the output field distribution. According to Collins’ diffraction integral formula [22] and using the similar way as the type I, the analytical expression of the output field distribution of a Hermite-Gaussian beam in FRFT plane is ob-

506

Daomu Zhao et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

tained

Table 1. The coefficients An and Bn with N = 10.

Eðx2 ; zÞ ¼

N P



1=2 

m=2

1 2 1 2 Cn W0 Cn n¼1       ip f p2 1  exp  x22 þ 2 ld d l d2 Cn " pffiffiffi  1=2 # i 2 px2 2 1 2  Hm l dW0 Cn W0 Cn An

ð9Þ where Cn ¼

Bn 1 ip þ 2þ 2 a W0 ld

  d 1 : f

ð11Þ

3. Numerical simulation comparisons For comparing the results obtained by using the analytical expression with those by using the diffraction integral formula, some numerical simulations are done. In calculation, we choose the coefficients An and Bn of a rectangular function [21] as the table 1. Fig. 2 and fig. 3 show the relative intensity distributions in FRFT plane of a Hermite-Gaussian beam (m ¼ 1) passing through the two types of apertured FRFT systems, respectively, where we choose a=W0 ¼ 1. From the figures we find

N 1 2 3 4 5 6 7 8 9 10

An 11.428 þ 0.95175i 0.06002  0.08013i 4.2743  8.5562i 1.6576 þ 2.7015i 5.0418 þ 3.2488i 1.1227  0.68854i 1.0106  0.26955i 2.5974 þ 3.2202i 0.1484  0.31193i 0.2085  0.23851i

Bn 4.0697 þ 0.22726i 1.1531  20.933i 4.4608 þ 5.1268i 4.3521 þ 14.997i 4.5443 þ 10.003i 3.8478 þ 20.087i 2.528  10.31i 3.3197  4.8008i 1.9002  15.82i 2.634 þ 25.009i

that the simulation results obtained by using the analytical expression are in keeping with those by using the numerical integral calculation quite well except for some slight decinations outside the center part, and the two types of apertured FRFT systems are no longer equivalent. When p  0:5, the simulation results have enough high consistency, especially on the Fourier transform plane. It should be pointed out that the numerical calculation efficiency by using the analytical expression is over 100 times than that by using the numerical integral calculation. Meanwhile, it should also be pointed out that the calculation accuracy could be improved if the N is increased. Furthermore, after some numerical simulations we find that the diffraction effect of the lens aperture can be neglected when a=W0 > 4.

Fig. 2. Relative intensity distribution in the FRFT output plane of the type I with the incidence of a HermiteGaussian beam (m ¼ 1) by using the two methods. The solid lines mean the case by using the analytical formula, and the dotted lines denote the case by using the diffraction integral formula. a) p ¼ 1; b) p ¼ 0.5; c) p ¼ 0.2; d) p ¼ 0.1.

Daomu Zhao et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

507

Fig. 3. Relative intensity distribution in the FRFT output plane of the type II with the incidence of a Hermite-Gaussian beam (m ¼ 1) by using the two methods. The others are the same as fig. 2.

4. Conclusions In conclusion, the apertured FRFT system is introduced and applied to treat the propagation of Hermite-Gaussian beams, and the analytical expressions for a Hermite-Gaussian beam through two types of apertured FRFT systems are obtained and compared with those by using the numerical integral calculation directly. It is shown that our method can significantly improve the numerical calculation efficiency. The results can be extended to the two-dimensional case, and the propagation of a Gaussian beam or a plane wave passing through the apertured FRFT systems can be regarded as the special cases. Acknowledgements. This work is supported by the National Natural Science Foundation of China (NSAF United Foundation), Grant No. 10276034.

References [1] Namias V: The fractional order Fourier transform and its applications to quantum mechanics. J. Inst. Math. Its Appl. 25 (1980) 241–265 [2] McBride AC, Kerr FH: On Namias fractional Fourier transform. IMA J. Appl. Math. 39 (1987) 159–175 [3] Mendlovic D, Ozaktas HM: Fractional Fourier transforms and their optical implementation: I. J. Opt. Soc. Am. A 10 (1993) 1875–1881

[4] Ozaktas HM, Mendlovic D: Fractional Fourier transforms and their optical implementation:II. J. Opt. Soc. Am. A 10 (1993) 2522–2531 [5] Lohmann AW: Image rotation, Wigner rotation, and the fractional Fourier transform. J. Opt. Soc. Am. A 10 (1993) 2181–2186 [6] Lohmann AW, Mendlovic D, Zalevsky Z: Fractional transformations in optics. In Wolf E (ed.): Progress in optics. Elsevier, Amsterdam, 1998 [7] Ozaktas HM, Zalevsky Z, Kutay MA: The Fractional Fourier Transform with applications in optics and signal processing. Wiley, New York 2000 [8] Torre A: The fractional Fourier transfomr and some of its applications to optics. In Wolf E (ed.): Progress in optics. Elsevier, Amsterdam 2002 [9] Mendlovic D, Ozaktas HM, Lohmann AW: Graded index fibers, Wigner-distribution functions, and the fractional Fourier transforms. Appl. Opt. 33 (1994) 6188–6193 [10] Pallat-Finet P: Fresnel diffraction and the fractionalorder Fourier transforms. Opt. Lett. 19 (1994) 1388– 1390 [11] Ozaktas HM, Mendlovic D: Fractional Fourier transforms as a tool for analyzing beam propagation and spherical mirror resonators. Opt. Lett. 19 (1994) 1678–1680 [12] Shih C: Optical interpretation of a complex-order Fourier transform. Opt. Lett. 20 (1995) 1178–1180 [13] Zalevsky Z, Mendlovic D, Dorsch RG: Gerchberg-Saxton algorithm applied in the fractional Fourier or the Fresnel domain. Opt. Lett. 21 (1996) 842–844 [14] Sahin A, Ozaktas HM, Mendlovic D: Optical implenations of two-dimensional fractional Fourier transform and linear canonical transform with arbitrary parameters. Appl. Opt. 37 (1998) 2130–2141 [15] Dragoman D, Dragoman M: Temporal implementation of Fourier-related transforms. Opt. Commun. 145 (1998) 33– 37

508

Daomu Zhao et al., Propagation of Hermite-Gaussian beams in apertured fractional Fourier transforming systems

[16] Zhao D: Multielement resonators and scaled fractional Fourier transforms. Opt. Commun. 168 (1999) 85–88 [17] Zhao D, Wang S: Effect of misalignment on optical fractional Fourier transforming systems. Opt. Commun. 198 (2001) 281–286 [18] Zhao D: A thick lens as a fractional Fourier transformer. Optik 111 (2000) 35–37 [19] Zhao D: Collins formula in frequency-domain described by fractional Fourier transform or fractional Hankel transform. Optik 111 (2000) 9–12

[20] Wang S, Zhao D: Matrix optics. Chapt. 2. CHEP-Springer, Beijing 2000 [21] Wen JJ, Breazeale MA: A diffraction beam field expressed as the superposition of Gaussian beams. J. Acoust. Soc. Am. 83 (1988) 1752–1756 [22] Collins SA: Lens-system diffraction integral written in terms of matrix optics. J. Opt. Soc. Am. 60 (1970) 1168– 1177