OSA's Digital Library

Optics Express

Optics Express

  • Editor: C. Martijn de Sterke
  • Vol. 17, Iss. 17 — Aug. 17, 2009
  • pp: 14825–14831
« Show journal navigation

Comparison of iterative angular spectrum and optimal rotation angle methods in designing beam-fanners

Seyyed H. Kazemi, Mir M. Mirsalehi, and Amir R. Attari  »View Author Affiliations


Optics Express, Vol. 17, Issue 17, pp. 14825-14831 (2009)
http://dx.doi.org/10.1364/OE.17.014825


View Full Text Article

Acrobat PDF (185 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We compare the iterative angular spectrum (IAS) and the optimal rotation angle (ORA) methods in designing two-dimensional finite aperture diffractive optical elements (FADOEs) used as beamfanners. The transfer functions of both methods are compared analytically in the spatial frequency domain. We have designed several structures of 1-to-4 and 1-to-6 beamfanners to investigate the differences in the performance of the beamfanners designed by ORA method for near field operation. Using the three-dimensional finite difference time-domain (3-D FDTD) method with perfect matched layer (PML) absorbing boundary condition (ABC), the diffraction efficiency is calculated for each designed FADOE and the corresponding values are compared.

© 2009 OSA

1. Introduction

In bidirectional methods, analysis and inverse-analysis are used iteratively to improve the DOE structure. Gerchberg–Saxton method [7

7. R. W. Gerchberg and W. O. Saxton, “Practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik (Stuttg.) 35, 237–250 (1972).

], iterative Fourier transform algorithm (IFTA) [8

8. F. Wyrowski and O. Bryngdahl, “Iterative Fourier-transform algorithm applied to computer holography,” J. Opt. Soc. Am. A 5(7), 1058–1065 (1988). [CrossRef]

], the method proposed by Di [9

9. F. Di, Y. Yingbai, J. Guofan, T. Qiaofeng, and H. Liu, “Rigorous electromagnetic design of finite-aperture diffractive optical elements by use of an iterative optimization algorithm,” J. Opt. Soc. Am. A 20(9), 1739–1746 (2003). [CrossRef]

], vectorial beam shaping by Jabbour [10

10. T. G. Jabbour and S. M. Kuebler, “Vectorial beam shaping,” Opt. Express 16(10), 7203–7213 (2008). [CrossRef] [PubMed]

] and iterative angular spectrum algorithm (IASA) [11

11. S. D. Mellin and G. P. Nordin, “Limits of scalar diffraction theory and an iterative angular spectrum algorithm for finite aperture diffractive optical element design,” Opt. Express 8, 706–722 (2001). [CrossRef]

] are various kinds of bidirectional methods. The Gerchberg–Saxton method and the IFTA method are used for far and near field cases. The method proposed by Di uses FDTD and thin element approximation (phase only approximation) as analysis and inverse-analysis methods, respectively and AS approach for wave propagation simulation. In the IASA, the optical field in the output plane is calculated using phase only approximation. This field is then propagated to the observation plane by AS approach.

In this paper, we compare the IAS and ORA methods in designing FADOEs as 2-D beamfanners. Since the original IAS method was presented for 1-D case, we have extended it for the 2-D case. The extended iterative angular spectrum (IAS) method is described in section 2 that is followed by a description of the ORA method. The transfer functions of these methods are compared in section 4. Both methods, as well as the FDTD method, are used to design and analyze 1-to-4 and 1-to-6 beamfanners in section 5 which is followed by conclusions.

2. Iterative angular spectrum method

The geometry of a 2-D finite aperture DOE is shown in Fig. 1
Fig. 1 Geometry of 2D finite aperture DOE and the two regions used in the IAS method.
. The input plane that includes the aperture is located in front of the DOE and the observation plane is located at z=z obs. It is assumed that a plane wave propagating along the z-axis interacts with the DOE at the input plane and it is expected to produce the desired pattern at the observation plane.

The space between the input and the observation plane is divided to two regions: DOE region and AS region. The DOE region extends from the input plane to the output plane (P0) and includes DOE, while the AS region extends from the output plane located at z=0 to the observation plane (P1). The thickness of DOE region is equal to the maximum etch depth (d max). The refraction index of the dielectric material used is n 2, while the surrounding medium is assumed to be air (n 1=1).

The transmission function of the DOE is
t(x,y)=τexp[jϕ(x,y)]Π(x/Lx,y/Ly),
(1)
where τ=2n 1/(n1+n 2) is Fresnel transmission coefficient, Π(x/Lx, y/Ly) is equal to 1 for any point inside a rectangular aperture of Lx×Ly, and equal to zero for other points. The phase function according to phase only approximation is φ(x,y)=k 0Δnd(x,y), where Δn is the difference between the two refractive indices, k 0=2π/λ, λ is the wavelength in free space and d(x,y) is the etch depth. We assume an incident plane wave, u inc, that is propagating normal to the aperture. The wave at the output plane of the first region is u 0=t(x,y)u inc. The angular spectrum of u 0, i.e., U 0, is defined as [4

4. J. W. Goodman, Introduction to Fourier Optics (Roberts & Company, Englewood, 2005).

]
U0(fx,fy)=P0u0(x,y)exp[j2π(fxx+fyy)]dxdy.
(2)
Using the fast Fourier transform (FFT) technique, U 0 can be computed. Assuming a square aperture, Lx=Ly=L, the distance between two adjacent points in spatial frequency domain is δf=δfx=δfy=1/L. To improve the resolution in the frequency domain, we can compute the FFT on an area larger than the DOE area. This is particularly useful if the desired intensity profile in the observation plane extends beyond the DOE area.

The field distribution in the observation plane, u 1, can be obtained by
u1(x,y)=U0exp[jk2zobs1(λ2fx)2(λ2fy)2]exp[j2π(fxx+fyy)]dfxdfy,
(3)
where k 2=n 2 k 0 and λ2=λ/n 2 are the wave-number and wavelength in the AS region, respectively [4

4. J. W. Goodman, Introduction to Fourier Optics (Roberts & Company, Englewood, 2005).

]. Using the inverse fast Fourier transform (IFFT), u 1 can be calculated as u 1=IFFT(HU 0), where
H(fx,fy)=exp[jk2zobs1(λ2fx)2(λ2fy)2]
(4)
is the angular spectrum transfer function. The maximum spatial frequency is f max=1/(2δ), where δ is the minimum feature size of the DOE. To consider all propagating waves, the condition f max>1/λ2, which is equivalent to δ <λ2/2, should be satisfied.

If the intensity at the observation plane satisfies the desired specification, the design process is terminated and the obtained DOE is used. Otherwise, the amplitude (not phase) of the field at the observation plane, u 1, is modified such that the difference between the obtained amplitude and the desired amplitude is reduced. The modified field distribution is then propagated backward to find the corresponding field distribution at the DOE output plane. Also, using the phase of the back going filed, the DOE structure is modified to satisfy the new field distribution. This modified DOE is quantized and used to repeat the above procedure.

3. Optimal rotation angle method

In this section, the ORA method is briefly described. For detailed information, the reader is referred to [5

5. J. Bengtsson, “Design of fan-out kinoforms in the entire scalar diffraction regime with an optimal-rotation-angle method,” Appl. Opt. 36(32), 8435–8444 (1997). [CrossRef]

] and [6

6. J. Stigwall and J. Bengtsson, “Design of array of diffractive optical elements with inter-element coherent fan-outs,” Opt. Express 12(23), 5675–5683 (2004). [CrossRef] [PubMed]

]. The geometry used in this method is shown in Fig. 2
Fig. 2 Geometry of DOE and the spots used in the ORA method.
. The output of the DOE is assumed to be on the left plane and it is divided to small cells. The output field at the center of the k-th cell is
u0,k=Akincexp[j(ϕkinc+ϕk)],
(5)
where Akincand ϕkincare the amplitude and phase of the incident field in the center of the k-th pixel, (xk , yk), respectively, and φk is the amount of phase modulation produced by the DOE. Considering the m-th spot in the observation plane, located at (xm , ym , z obs), the field at this spot can be represented as
u1,m=ku0,khk,m,
(6)
where hk,m is a complex value quantity that relates the field at the k-th cell of the DOE to the field at the m-th spot of the observation plane. Using the Helmholtz-Kirchhoff integral, it is obtained as [5

5. J. Bengtsson, “Design of fan-out kinoforms in the entire scalar diffraction regime with an optimal-rotation-angle method,” Appl. Opt. 36(32), 8435–8444 (1997). [CrossRef]

]
hk,m=14π[jk2zobs(jk21/rkmc)/rkmc]4exp(jk2rkmc)rkmcsin(k˜xa2)sin(k˜yb2)k˜xk˜y,
(7)
where rkmc=[(xk-xm)2+(yk-ym)2+z obs 2]½, k˜x=k2(xkxm)/rkmc, k˜y=k2(ykym)/rkmc, and a and b are the x and y dimensions of each DOE cell, respectively.

According to Bengtesson's ORA method, the intensities of the observation spots will be maximized if the phase of each DOE cell is changed by Δφk and this procedure is repeated until the phases of the cells are converged. For the formula of Δφk, the reader is reffered to [5

5. J. Bengtsson, “Design of fan-out kinoforms in the entire scalar diffraction regime with an optimal-rotation-angle method,” Appl. Opt. 36(32), 8435–8444 (1997). [CrossRef]

].

4. Transfer functions

The differences in the amplitude can be divided in two types. First, for zobs smaller than 10λ the frequencies near edge are amplified. Second, for z obs greater than 10λ, the spatial frequency bandwidth of the ORA transfer function becomes smaller (See Fig. 3
Fig. 3 The middle cross sections of the transfer functions for IAS (...) and ORA (___) methods at zobs=10λ, 30λ and 100λ
). The first difference produces high spatial frequency components in the DOE structures designed by ORA method, which can destructively affect the field distribution. Figure 4
Fig. 4 The middle cross sections of typical DOE profiles designed by (a) IAS and (b) ORA methods for zobs=10λ case.
shows the effect of these high spatial frequency components in two typical DOE profiles. The second difference eliminates high spatial frequencies that are important specially in cases that the spots are located at a large distance from each other compared with the DOE dimensions.

Therefore, the designs obtained by the IAS method for the near field cases are superior to the corresponding designs obtained by the ORA method. For a specific pattern, as the distance between the observation plane and the DOE is increased, the results obtained by the two methods become more similar.

5. Design and analysis

We investigate the effect of high spatial frequency components in the design of 1-to-4 and 1-to-6 beamfanners by IAS and ORA methods. The DOE is assumed to be a square of dimension D=10λ0 made of silicon with a refractive index of n 2=3.4 and λ0 is set to 5μm. The spatial quantization steps along the three directions (x, y and z), i.e., a, b, and δ, and the etch depth quantization step are chosen equal to λ0/25=.0.2μm. As mentioned in section 2, to improve the frequency resolution, we put additional zeroes around the above matrix to increase its dimension by a factor of 10. As a result, the spatial frequency resolution will be δf=1/(100λ0) and therefore, the transmission function is represented by a 2500×2500 complex-valued matrix.

In designing 1-to-4 and 1-to-6 beamfanners by the ORA method, only 4 and 6 spots of equal intensities are defined at the observation plane, while in the IAS method, the desired intensity profile at the observation plane is used. We used a Gaussian shape as the desired intensity profile at the observation plane. The Gaussian functions were defined as exp{-[(x-xm)2+(y-ym)2]/(2w 0 2)]}, where (xm,ym) corresponds to the center of the m-th spot. In our simulations, we have used w 0= 0.14λ. This width is very small compared to the distance between the spot points (3λ to 8λ). Therefore, the two distribution profiles are approximately the same in numerical calculations.

To analyze the designed DOEs, we used a 3-D FDTD algorithm with unsplit step PML absorbing boundary condition. The Analyzed region is divided to the total field and the scattered filed regions, as shown in Fig. 5
Fig. 5 Cross section of the analyzed region.
. We have used the technique introduced by Sullivan for implementing the PML [12

12. D. M. Sullivan, “An unsplit step 3-D PML for use with the FDTD method,” IEEE Microwave and Guided Wave Letters 7(7), 184–186 (1997). [CrossRef]

, 13

13. D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, IEEE Press series on RF and microwave technology, 2000.

].

The mesh size for the FDTD simulation was chosen equal for all three directions, i.e., δxyz=δ. Also, δt=δ/(2c), where c is the velocity of light in free space, was chosen as the time step. The depths of PML and scattered field were set to 7δ. The incident field was assumed to be a TEM wave with a sharp Gaussian time varying amplitude, E xinc=exp{-(t-t d)2/(2w t 2)}, where t d and w t are the time delay and the spread parameter of the Gaussian pulse, respectively. Using the FDTD method, the distributions of the electric and magnetic fields were obtained for 200 time steps.

After the fields were found at P 0 plane by the FDTD method, the angular spectrum technique was used to obtain the fields at the observation plane, i.e., P 1. The distribution of light intensity, I, at P 1 is then calculated. The diffraction efficiency defined as
DE=WI(x,y)dxdyP1I(x,y)dxdy
(10)
is then computed, where W is the detectors window which consists of four (six) circular areas of radius equal to λ0, centered at the four (six) spots positions.

The diffraction efficiency (DE) obtained by FDTD method are provided in Table 1

Table 1. The diffraction efficiency (DE) obtained by FDTD method for designed 1-to-4 and 1-to-6 beamfanners.

table-icon
View This Table
for 1-to-4 and 1-to-6 beamfanners, where S4 and S6 are separation distance of adjacent spots in the 1-to-4 and 1-to-6 beamfanners, respectively. It can be seen, that the IAS method provides better designs than the ORA.

6. Conclusions

We compared the transfer functions of the ORA and IAS methods in the near field region. This comparison was made in spatial frequency domain. We showed that, there are two weak points in the ORA transfer functions including 1) stronger high spatial frequency components and 2) smaller frequency bandwidths. By choosing the distance between the observation plane and the FADOE plane (z obs) as an important variable, we found, that if z obs is less than or equal to 10λ, the high spatial frequency components appear stronger in the ORA method. This results in increasing the number of discontinuities in the DOE phase profile. Vice versa, by increasing z obs, the ORA frequency bandwidth becomes smaller. This puts a limit in the separating distance between the spots.

To investigate these concepts, we designed several 1-to-4 and 1-to-6 beamfanners using both methods and compared their diffraction efficiencies. In all designs, the IAS method outperforms the ORA method in the near field region.

Acknowledgments

The authors are thankful to Iran Telecommunication Research Center (ITRC) for supporting this research.

References and links

1.

J. Jiang and G. P. Nordin, “A rigorous unidirectional method for designing finite aperture diffractive optical elements,” Opt. Express 7(6), 237–242 (2000). [CrossRef] [PubMed]

2.

D. W. Prather, M. S. Mirotznik, and J. N. Mait, “Boundary integral methods applied to the analysis of diffractive optical elements,” J. Opt. Soc. Am. A 14(1), 34–43 (1997). [CrossRef]

3.

D. W. Prather, J. N. Mait, M. S. Mirotznik, and J. P. Collins, “Vector-based synthesis of finite aperiodic subwavelength diffractive optical elements,” J. Opt. Soc. Am. A 15(6), 1599–1607 (1998). [CrossRef]

4.

J. W. Goodman, Introduction to Fourier Optics (Roberts & Company, Englewood, 2005).

5.

J. Bengtsson, “Design of fan-out kinoforms in the entire scalar diffraction regime with an optimal-rotation-angle method,” Appl. Opt. 36(32), 8435–8444 (1997). [CrossRef]

6.

J. Stigwall and J. Bengtsson, “Design of array of diffractive optical elements with inter-element coherent fan-outs,” Opt. Express 12(23), 5675–5683 (2004). [CrossRef] [PubMed]

7.

R. W. Gerchberg and W. O. Saxton, “Practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik (Stuttg.) 35, 237–250 (1972).

8.

F. Wyrowski and O. Bryngdahl, “Iterative Fourier-transform algorithm applied to computer holography,” J. Opt. Soc. Am. A 5(7), 1058–1065 (1988). [CrossRef]

9.

F. Di, Y. Yingbai, J. Guofan, T. Qiaofeng, and H. Liu, “Rigorous electromagnetic design of finite-aperture diffractive optical elements by use of an iterative optimization algorithm,” J. Opt. Soc. Am. A 20(9), 1739–1746 (2003). [CrossRef]

10.

T. G. Jabbour and S. M. Kuebler, “Vectorial beam shaping,” Opt. Express 16(10), 7203–7213 (2008). [CrossRef] [PubMed]

11.

S. D. Mellin and G. P. Nordin, “Limits of scalar diffraction theory and an iterative angular spectrum algorithm for finite aperture diffractive optical element design,” Opt. Express 8, 706–722 (2001). [CrossRef]

12.

D. M. Sullivan, “An unsplit step 3-D PML for use with the FDTD method,” IEEE Microwave and Guided Wave Letters 7(7), 184–186 (1997). [CrossRef]

13.

D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, IEEE Press series on RF and microwave technology, 2000.

OCIS Codes
(000.4430) General : Numerical approximation and analysis
(050.1220) Diffraction and gratings : Apertures
(050.1970) Diffraction and gratings : Diffractive optics

ToC Category:
Diffraction and Gratings

History
Original Manuscript: May 4, 2009
Revised Manuscript: July 1, 2009
Manuscript Accepted: July 24, 2009
Published: August 6, 2009

Citation
Seyyed H. Kazemi, Mir M. Mirsalehi, and Amir R. Attari, "Comparison of iterative angular spectrum and optimal rotation angle methods in designing beam-fanners," Opt. Express 17, 14825-14831 (2009)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-17-17-14825


Sort:  Author  |  Year  |  Journal  |  Reset  

References

  1. J. Jiang and G. P. Nordin, “A rigorous unidirectional method for designing finite aperture diffractive optical elements,” Opt. Express 7(6), 237–242 (2000). [CrossRef] [PubMed]
  2. D. W. Prather, M. S. Mirotznik, and J. N. Mait, “Boundary integral methods applied to the analysis of diffractive optical elements,” J. Opt. Soc. Am. A 14(1), 34–43 (1997). [CrossRef]
  3. D. W. Prather, J. N. Mait, M. S. Mirotznik, and J. P. Collins, “Vector-based synthesis of finite aperiodic subwavelength diffractive optical elements,” J. Opt. Soc. Am. A 15(6), 1599–1607 (1998). [CrossRef]
  4. J. W. Goodman, Introduction to Fourier Optics (Roberts & Company, Englewood, 2005).
  5. J. Bengtsson, “Design of fan-out kinoforms in the entire scalar diffraction regime with an optimal-rotation-angle method,” Appl. Opt. 36(32), 8435–8444 (1997). [CrossRef]
  6. J. Stigwall and J. Bengtsson, “Design of array of diffractive optical elements with inter-element coherent fan-outs,” Opt. Express 12(23), 5675–5683 (2004). [CrossRef] [PubMed]
  7. R. W. Gerchberg and W. O. Saxton, “Practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik (Stuttg.) 35, 237–250 (1972).
  8. F. Wyrowski and O. Bryngdahl, “Iterative Fourier-transform algorithm applied to computer holography,” J. Opt. Soc. Am. A 5(7), 1058–1065 (1988). [CrossRef]
  9. F. Di, Y. Yingbai, J. Guofan, T. Qiaofeng, and H. Liu, “Rigorous electromagnetic design of finite-aperture diffractive optical elements by use of an iterative optimization algorithm,” J. Opt. Soc. Am. A 20(9), 1739–1746 (2003). [CrossRef]
  10. T. G. Jabbour and S. M. Kuebler, “Vectorial beam shaping,” Opt. Express 16(10), 7203–7213 (2008). [CrossRef] [PubMed]
  11. S. D. Mellin and G. P. Nordin, “Limits of scalar diffraction theory and an iterative angular spectrum algorithm for finite aperture diffractive optical element design,” Opt. Express 8, 706–722 (2001). [CrossRef]
  12. D. M. Sullivan, “An unsplit step 3-D PML for use with the FDTD method,” IEEE Microwave and Guided Wave Letters 7(7), 184–186 (1997). [CrossRef]
  13. D. M. Sullivan, Electromagnetic Simulation Using the FDTD Method, IEEE Press series on RF and microwave technology, 2000.

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.

Figures

Fig. 1 Fig. 2 Fig. 3
 
Fig. 4 Fig. 5
 

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited